马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
n=10; % the number of the mesh
h=10; % the heat transfer coefficient
k=1;
dx=1/n;
dy=1/n;
T_old=zeros(n+1,n+1)+500;
T_new=T_old;
eps=1;
iter=0;
while (eps>0.001)
iter=iter+1;
for j=2:n
i=1;
T_new(i,j)=1/(2*(h*dx/k+2))*(2*T_old(i+1,j)+T_new(i,j-1)+T_old(i,j+1)+2*h*dx/k*300);
end
for i=2:n
for j=2:n
T_new(i,j)=1/4*(T_old(i+1,j)+T_old(i,j+1)+T_new(i-1,j)+T_new(i,j-1));
end
end
ErrorTOT=T_old-T_new;
eps=max(max(ErrorTOT));
T_old=T_new;
end
[XI YI]=meshgrid(0:dx:1,0:dy:1);
figure(1)
contourf(XI,YI,T_new,30);
caxis([300 500]);
colorbar('vert')
title('Temp Distribution at S-S')
% 不行用while,怎么替换谢谢 |