Методы математической физики. Куликова И.В - 24 стр.

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

Составители: 

Рубрика: 

24
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
i=I;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
end
for i=2:I-1
j=1;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
j=J;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
end
end
% Решение СЛАУ
TT=b/a';
T=zeros(I,J,S);
Ttime=zeros(1,S);
for s=1:S
for i=1:I
for j=1:J
T(i,j,s)=TT(i+(j-1)*I+(s-1)*I*J);
end
end
end
for s=1:S
figure
meshc(Y,X,T(:,:,s))
grid on
xlabel('y, m')
ylabel('x, m')
zlabel('T, K')
axis([min(Y) max(Y) min(X) max(X) (min(min(min(T)))-10)
max(max(max(T)))])
pause (0.1)
M(s)=getframe;
end
figure
movie(M,5,3)
for s=1:S
Ttime(s)=T(ceil(I/2),ceil(J/2),s);
end
figure
plot(time,Ttime)
grid on
xlabel('t, c')
ylabel('T, K')

Страницы