Xlabel('t');
figure;
======================================================================
% синусоидальный нагрев
tp=1e-6; % период колебаний, [c]
t=0:ti/100:12*tp; % задаём временной интервал от 0 до tp в долях ti, [c]
N=100; % количество членов суммы ряда
Hsur=plot(t,normalize_v(sin(t/tp)),'r',t,normalize_v(sin(t/tp+1/tp)),'b');
set(Hsur, 'LineWidth', 3);
ylabel('{\theta}(t,x), q(t)');
Xlabel('t');
% периодический импульсный нагрев
tp=1e-6; % период колебаний, [c]
ti=tp/2; % длительность импульса, [c]
t=0:ti/100:2*tp; % задаём временной интервал от 0 до tp в долях ti, [c]
N=1000; % количество членов суммы ряда
Hsur=plot(t,normalize_v(real(q1(N,t,ti,tp))),'r',t,normalize_v(teta12(N,x,l,a,alfa,D,t,ti,tp)),'b');
set(Hsur, 'LineWidth', 3);
ylabel('{\theta}(t,x), q(t)');
Xlabel('t');
figure;
script % демонстрационный файл функции teta12(t,x)
% =========================================================================
% Первая характерная форма
%
N=100; % количество членов суммы ряда
tp=1e-6; % период колебаний, [c]
ti=tp/2; % длительность импульса, [c]
l=0.001; % толщина пластины, [см]
a=0.0162; % коэффициент температуропроводности, [см^2/c]
alfa=0.318; % линеаризованный коэффициент теплопотерь, [Вт/см^2]
D=3.2; % тепловая активность [(Вт/см^2 K^2)*c]
t=0:ti/50:2*tp; % задаём временной интервал от 0 до 2tp в долях ti, [c]
% Строим поверхность teta=teta(t,x)
% для этого создаём матрицу поверхности - A
% далее заполняем её в цикле соответствующими значениями t и х
% первая характерная форма
tp=1e-6; % период колебаний, [c]
ti=tp/2; % длительность импульса, [c]
t=0:ti/100:2*tp; % задаём временной интервал от 0 до tp в долях ti, [c]
A = [];
i=1;
for x=0:l/20:l;
c = teta12(N,x,l,a,alfa,D,t,ti,tp);
A (i, :) = c; i=i+1;
end;
% задаём параметры и рисуем поверхность
Hsur=mesh(A); % рисуем поверхность
colormap('copper'); % цвет поверхности
set(Hsur, 'FaceColor', 'interp'); % способ окраски
set(Hsur, 'LineStyle', 'none');
set(Hsur, 'FaceAlpha', [0.8]); % полупрозрачность
title('periodical impulse | thin film | tp=1e-6 s, l=10 mkm'); % название графика
set(gca, 'Projection', 'perspective'); % способ проецирования
set(gca, 'Color', [0 1 1]); % цвет фона
Xlabel('t axis'); % название осей
ylabel('x axis');
zlabel('{\theta}(t,x)');
View(-53,23); % угол разворота графика
colorbar; % цветной индикатор значений
script % демонстрационный файл функции teta12(t,x)
% =========================================================================
% Вторая характерная форма
%
N=100; % количество членов суммы ряда
l=0.001; % толщина пластины, [см]
a=0.0162; % коэффициент температуропроводности, [см^2/c]
alfa=0.318; % линеаризованный коэффициент теплопотерь, [Вт/см^2]
D=3.2; % тепловая активность [(Вт/см^2 K^2)*c]
% вторая характерная форма
tp=1e-3; % период колебаний, [c]
ti=tp/2; % длительность импульса, [c]
t=0:ti/100:2*tp; % задаём временной интервал от 0 до tp в долях ti, [c]
% Строим поверхность teta=teta(t,x)
% для этого создаём матрицу поверхности - A
% далее заполняем её в цикле соответствующими значениями t и х
A = [];
i=1;
for x=0:l/20:l;
c = teta12(N,x,l,a,alfa,D,t,ti,tp);
A (i, :) = c; i=i+1;
end;
% задаём параметры и рисуем поверхность
Hsur=mesh(A); % рисуем поверхность
colormap('copper'); % цвет поверхности
set(Hsur, 'FaceColor', 'interp'); % способ окраски
set(Hsur, 'LineStyle', 'none');
set(Hsur, 'FaceAlpha', [0.8]); % полупрозрачность
title('periodical impulse | thin film | tp=1e-3 s, ti=tp/2 s, l=10 mkm'); % название графика
set(gca, 'Projection', 'perspective'); % способ проецирования
set(gca, 'Color', [0 1 1]); % цвет фона