Xlabel('t axis'); % название осей
ylabel('x axis');
zlabel('{\theta}(t,x)');
View(-53,23); % угол разворота графика
colorbar; % цветной индикатор значений
==============================================================
script
% параметры
N=1000; % количество цленов суммы ряда
a=0.318; % коэффициент теплопроводности [см^2/c]
alfa=0.0162; %
D=3.27 ; % тепловая активность [(Дж/см^2 К)^2/с]
x=0 % [см]
x =
0
% Характерные формы импульсов для пластин разной толщины
l=0.0010; % толщина пластины, [см]
% первая характерная форма
tp=1e+3; % период колебаний, [c]
ti=tp/10; % длительность импульса, [c]
t=0:ti/100:2*tp; % задаём временной интервал от 0 до tp в долях ti, [c]
st1 = num2str(ti); st2 = num2str(tp); st3 = num2str(l);
st = strcat('ti =', st1, ' s, tp =', st2, ' s, l =', st3, ' sm');
Hsur=plot(t,normalize_v(teta12(N,x,l,a,alfa,D,t,ti,tp))); title(st);
Hsur=plot(t,normalize_v(teta12(N,x,l,a,alfa,D,t,ti,tp))); title(st);
set(Hsur, 'LineWidth', 2);
==============================================================
N=1000;
a=0.0162;
alfa=0.318;
D=3.2;
tp=1e-5;
ti=tp/2;
l=0.01;
x=0;
t=ti;
y=0:l/103:l;
x=y.*y;
teta=0;
for k=1:N
omega=2*pi/tp;
gamma=ti/tp;
ksik=pi*k*gamma;
ck1=(sin(ksik)/(pi*k))*exp(-i*ksik);
h=sqrt(omega*D)/alfa;
hk1 = exp(-x*sqrt(i*k*omega/a))/(1+h*sqrt(i*k));
teta = teta + real((ck1)*(hk1).*exp(i*omega*k*t));
end;
teta=(2/alfa)*teta;
plot(teta);
xlabel('x^2 axis');
ylabel('teta axis');