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

UptoLike

ВУЗ: 

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

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

Рубрика: 

19
j=J;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = 1;
b(n) = Tenv;
j=1;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = -1;
a(n,n+I) = 1;
end
end
for s=2:S-1
for j=1:J
i=1;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = -1;
a(n,n+1) = 1;
i=I;
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) = -1;
a(n,n-1) = 1;
end
end
% Решение СЛАУ
TT=b/a';
T=zeros(I,J,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
% Вывод графиков по слоям по оси Z
for s=1:S
figure
mesh(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))))
max(max(max(T)))])
pause (0.1)
end
% Вывод 3D графика
figure
slice(Y, X, Z, T, [Y(fix(3*J/4))], [X(fix(3*I/4))],
[Z(fix(S/2))], length(X))
grid on
colorbar
xlabel(
'y, m')
ylabel('x, m')
zlabel('z, m')
title('T, K')

Страницы