哈工程传热学数值计算大作业(2)

二维稳态导热数值计算,matlab

第一题：

a=linspace(0,0.6,121); t1=[60+20*sin(pi*a/0.6)]; t2=repmat(60,[80 121]); s=[t1;t2]; %构造矩阵

for k=1:10000000 %理论最大迭代次数，想多大就设置多大 S=s;

for j=2:120 for i=2:80

S(i,j)=0.25*(S(i-1,j)+S(i+1,j)+S(i,j-1)+S(i,j+1)); end end

if norm(S-s)<0.0001

break; %如果符合精度要求，提前结束迭代 else s=S; end end

S %输出数值解

数值解数据量太大，这里就不打印出来，只画出温度分布。

画出温度分布： figure(1)

xx=linspace(0,0.6,121); yy=linspace(0.4,0,81); [x,y]=meshgrid(xx,yy); surf(x,y,S)

axis([0 0.6 0 0.4 60 80]) grid on xlabel('L1') ylabel('L2')

zlabel('t（温度）')