欢迎来到天天文库
浏览记录
ID：41152351
大小：183.16 KB
页数：5页
时间：2019-08-17
《一维热传导方程分离变量法与有限差分法Matlab解法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、模拟与仿真根据课上所学知识，我们有如下方程：uau−2=0,0<0txxuux=0,xxl=0,t>0=0=u=ϕ(),x0<2、^2/4*t));end;surf(x,t,u);xlabel('x'),ylabel('t'),zlabel('T');title('分离变量法（无穷）');disp(u);得到如图所示的热传导方程：有限差分法u=zeros(20,100);%t=1x=pi20行100列横坐标为x纵坐标为ts=(1/100)/(pi/20)^2;fprintf('稳定性系数S为:');disp(s);fori=1:20u(i,1)=i/20*pi;;end;forj=1:100u(1,j)=0;endforj=1:99fori=2:19u(i,j+1)=s*u(i+3、1,j)+(1-2*s)*u(i,j)+s*u(i-1,j);endendforj=1:100u(20,j)=u(19,j);end;disp(u);[x,t]=meshgrid(1:100,1:20);surf(x,t,u);xlabel('t'),ylabel('x'),zlabel('T');title('有限差分法解');我们得到如图所示的热传导方程：结论：比较可得由以上两种方法作出的三维图形基本相同，符合热传导的热量分布随时间和空间的变化规律第四题完成
2、^2/4*t));end;surf(x,t,u);xlabel('x'),ylabel('t'),zlabel('T');title('分离变量法（无穷）');disp(u);得到如图所示的热传导方程：有限差分法u=zeros(20,100);%t=1x=pi20行100列横坐标为x纵坐标为ts=(1/100)/(pi/20)^2;fprintf('稳定性系数S为:');disp(s);fori=1:20u(i,1)=i/20*pi;;end;forj=1:100u(1,j)=0;endforj=1:99fori=2:19u(i,j+1)=s*u(i+
3、1,j)+(1-2*s)*u(i,j)+s*u(i-1,j);endendforj=1:100u(20,j)=u(19,j);end;disp(u);[x,t]=meshgrid(1:100,1:20);surf(x,t,u);xlabel('t'),ylabel('x'),zlabel('T');title('有限差分法解');我们得到如图所示的热传导方程：结论：比较可得由以上两种方法作出的三维图形基本相同，符合热传导的热量分布随时间和空间的变化规律第四题完成
此文档下载收益归作者所有
