传热学作业数值计算数值计算matlab程序内容:>> tw1=10; % 赋初值tw2=20;c=1.5;p2=20;p1=c*p2;L2=40;L1=c*L2;deltaX=L2/p2;a=p2+1;b=p1+1;ti=ones(a,b)*5;m1=ones(a,b);m1(a,2:b-1)=zeros(1,b-2);m1(2:a,1)=zeros(a-1,1);m1(2:a,b)=zeros(a-1,1);m1(1,:)=ones(1,b)*2;k=0;max1=1.0;tn=ti;while(max1>1e-6)max1=0;k=k+1;for i=1:1:afor j=1:1:bm=m1(i,j);n=ti(i,j);switch mcase 0tn(i,j)=tw1;case 1tn(i,j)=0.25*(tn(i,j+1)+tn(i,j-1)+tn(i+1,j)+tn(i-1,j));case 2tn(i,j)=tw1+tw2*sin(pi*(j-1)/(b-1));ender=abs(tn(i,j)-n);if er>max1max1=er;endendendti=tn;endktimax1t2=ones(a,b); %求解析温度场for i=a:-1:1for j=1:1:by=deltaX*(a-i);x=deltaX*(j-1);t2(i,j)=tw1+tw2*sin(pi*x/L1)*(sinh(pi*y/L1))/(sinh(pi*L2/L1));endendt2迭代次数k =706数值解温度场ti数值解每次迭代的最大误差max1 =9.8531e-07解析温度场t2取第11行的解析解和数值解的点曲线为第11行的解析解的直线,散点为其数值解的点第11行的误差=[数值解(11行) –解析解(11行)]/解析解数值温度场图像解析温度场图像数值解与解析解的误差数值计算matlab程序内容:>> tw1=10;tw2=20;c=1.5;p2=20;p1=c*p2;L2=20;deltaX=L2/p2;L1=c*L2;a=p2+1;b=p1+1;ti=ones(a,b)*5;m1=ones(a,b);m1(a,2:b-1)=zeros(1,b-2);m1(2:a,1)=zeros(a-1,1);m1(2:a,b)=zeros(a-1,1);m1(1,:)=ones(1,b)*2;k=0;max1=1.0;tn=ti;while(max1>1e-6)max1=0;k=k+1;for i=1:1:afor j=1:1:bm=m1(i,j);n=ti(i,j);switch mcase 0tn(i,j)=tw1;case 1tn(i,j)=0.25*(tn(i,j+1)+tn(i,j-1)+tn(i+1,j)+tn(i-1,j));case 2tn(i,j)=tw2;ender=abs(tn(i,j)-n);if er>max1max1=er;endendendti=tn;endktimax1tx=ones(a,b);for i=1:1:afor j=1:1:by=(a-i)*deltaX;x=(j-1)*deltaX;m=sym('m');g=(((-1)^(m+1)+1)/m)*sin(m*pi*x/L1)*sinh(m*pi*y/L1)/sinh(m*pi*L2/L1); h=symsum(g,m,1,100);tx(i,j)=2*h*(tw2-tw1)/pi+tw1;endendtx迭代次数k = 695数值解温度场ti数值解每次迭代的最大误差max1 =9.8243e-07解析温度场tx =取第11行的解析解和数值解的点曲线为第11行的解析解的直线,散点为其数值解的点第11行的误差=[数值解(11行) –解析解(11行)]/解析解数值温度场ti图像:解析温度场tx图像:数值解与解析解的误差数值计算matlab程序内容:>> t0=90;tf=10;L=10;c=0.25;p2=20;p1=p2/c;B=c*L;d=0.5*B;h=10;a=p2+1;b=p1+1;deltaX=B/p2;lambda=160;Bi=h*deltaX/lambda;ti=ones(a,b)*10;m1=ones(a,b)*3;m1(2:a-1,1)=zeros(a-2,1);m1(a,2:b-1)=ones(1,b-2);m1(1,2:b-1)=ones(1,b-2)*6;m1(2:a-1,b)=ones(a-2,1)*2;m1(1,b)=ones(1,1)*4;m1(a,b)=ones(1,1)*5;m1(1,1)=7;m1(a,1)=8;tn=ti;max1=1.0;k=0;while ( max1>1e-6)k=k+1;max1=0;for i=1:1:afor j=1:1:bm=m1(i,j);n=tn(i,j);switch mcase 0tn(i,j)=t0;case 1tn(i,j)=(2*tn(i-1,j)+tn(i,j-1)+tn(i,j+1)-4*tf)/(4+2*Bi)+tf;case 2tn(i,j)=(2*tn(i,j-1)+tn(i-1,j)+tn(i+1,j)-4*tf)/(4+2*Bi)+tf;case 3tn(i,j)=0.25*(tn(i,j-1)+tn(i,j+1)+tn(i-1,j)+tn(i+1,j));case 4tn(i,j)=(tn(i,j-1)+tn(i+1,j)-2*tf)/(2*Bi+2)+tf;case 5tn(i,j)=(tn(i,j-1)+tn(i-1,j)-2*tf)/(2*Bi+2)+tf;case 6tn(i,j)=(2*tn(i+1,j)+tn(i,j-1)+tn(i,j+1)-4*tf)/(4+2*Bi)+tf;case 7tn(i,j)=t0;case 8tn(i,j)=t0;ender=abs(tn(i,j)-n);if er>max1max1=er;endendendti=tn;endktita=ones(a,b);Bi1=h*d/lambda;sbi=sqrt(Bi1);for i=1:1:afor j=1:1:bif i>(a+1)/2y=-(i-(a+1)/2)*deltaX;else y=((a+1)/2-i)*deltaX;endx=deltaX*(j-1);ta(i,j)=(cosh(sbi*(L-x)/d)+sbi*sinh(sbi*(L-x)/d))*(t0-tf)/(cosh(sbi*L/d)+sbi*sinh(sbi*L/d))+tf;endendta迭代次数k =1461数值解温度场ti解析温度场ta取第11行的解析解和数值解的点曲线为第11行的解析解的直线,散点为其数值解的点第11行的误差=[数值解(11行) –解析解(11行)]/解析解图像如下数值温度场图像解析温度场图像数值解与解析解的误差数值计算matlab程序内容:>> tw=10;L2=15;c=0.75;L1=L2/c;p2=24 ;p1=p2/c;deltaX=2*L2/p2;a=p2+1;b=p1+1;lambda=16;qv0=24;ti=ones(a,b)*5;m1=ones(a,b);m1(1,:)=zeros(1,b);m1(2:a,b)=zeros(a-1,1);m1(2:a,1)=zeros(a-1,1);m1(a,2:b-1)=zeros(1,b-2);tn=ti;max1=1.0;k=0;while(max1>1e-6)max1=0;k=k+1;for i=1:1:afor j=1:1:bm=m1(i,j);n=tn(i,j);switch mcase 0tn(i,j)=tw;case 1tn(i,j)=0.25*(tn(i-1,j)+tn(i+1,j)+tn(i,j-1)+tn(i,j+1)+qv0*(deltaX^2)/lambda);ender=abs(tn(i,j)-n);if er>max1max1=er;endendendti=tn;endtik;tx=ones(a,b);for i=1:1:afor j=1:1:bif i>(a+1)/2y=-(i-(a+1)/2)*deltaX;elsey=((a+1)/2-i)*deltaX;endif j>(b+1)/2x=(j-(b+1)/2)*deltaX;elsex=-((b+1)/2-j)*deltaX;endm=sym('m');xi=(2*m-1)*pi/2;g=((-1)^m)/(xi^3)*(cosh(xi*y/L1)/cosh(xi*L2/L1))*cos(xi*x/L1); h=symsum(g,m,1,100);tx(i,j)=2*qv0*L1^2/lambda*h+qv0*(L1^2-x^2)/(2*lambda)+tw; endendtx数值温度场ti解析温度场tx取第13行的解析解和数值解的点曲线为第13行的解析解的直线,散点为其数值解的点第13行的误差=[数值解(13行) –解析解(13行)]/解析解数值温度场图像解析温度场图像数值解与解析解的误差。