一维稳态导热matlab程序求解

《一维稳态导热matlab程序求解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用文档文案大全实用文档文案大全实用文档文案大全实用文档文案大全实用文档文案大全实用文档设定N=100，a=3,b=40,T0=8,s=20x-4TMatlab源程序代码：length=2;N=100;iterN=0;eps=1.0e-4;deltQ=2;T2=(1:N);T2(1)=8;T2(N)=1;dx=length/(N-2);while(deltQ>=eps&&iterN<50)T=T2;fori=1:Nk(i)=3*T(i)^3;sc(i)=20*(i-1)/N;sp(i)=4;endk(N/2+1:N)=40*T(N/

2、2+1:N).^3;fori=2:N-1ifi==2a(i-1)=-2*k(i)/dx;c(i-1)=-2*k(i)*k(i+1)/dx/(k(i)+k(i+1));b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;d(i-1)=sc(i)*dx-a(i-1)*T(i-1);elseifi==N-1;a(i-1)=-2*k(i)*k(i-1)/dx/(k(i)+k(i-1));文案大全实用文档c(i-1)=-2*k(i+1)/dx;b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;d(i-1)=sc(i)*d

3、x-c(i-1);elsea(i-1)=-2*k(i)*k(i-1)/dx/(k(i)+k(i-1));c(i-1)=-2*k(i)*k(i+1)/dx/(k(i)+k(i+1));b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;d(i-1)=sc(i-1)*dx;endendendifabs(b(1))<=abs(c(1))

4、

5、abs(b(N-2)<=abs(c(N-2)))warndla('²»Âú×ã×·¸Ï·¨Ìõ¼þ','¾¯¸æ£¡');endfori=1:N-2ifabs(b(i))

6、bs(a(i)))warndla('²»Âú×ã×·¸Ï·¨Ìõ¼þ','¾¯¸æ£¡');endendfori=1:N-2ifi==1u(1)=b(1);elsel(i)=a(i)/u(i-1);u(i)=b(i)-l(i)*c(i-1);endendfori=1:N-2文案大全实用文档ifi==1y(1)=d(1);elsey(i)=d(i)-l(i)*y(i-1);endendforj=N-2:-1:1ifj==N-2x(N-2)=y(N-2)/u(N-2);elsex(j)=(y(j)-c(j)*x(j+1))/u(j);e

7、ndenditerN=iterN+1;T2(2:N-1)=x(1:N-2);deltQ=max(abs(T2-T));endfori=1:10Q(i,1:10)=T2((1+(i-1)*10):i*10);enddisp(Q);上面程序有乱码，附截图一份：文案大全实用文档文案大全实用文档文案大全