完整word版MATLAB编辑一维热传导方程的模拟程序.doc

(完整word版)MATLAB编辑一维热传导方程的模拟程序 (完整word版)MATLAB编辑一维热传导方程的模拟程序 PAGE / NUMPAGES (完整word版)MATLAB编辑一维热传导方程的模拟程序 求解以下热传导问题： 2T 1 T 0 z L z2 0 t T z,0 1 z2 T 0, t 1, T L,t 0 L 1, 1 程序： function heat_conduction() % 一维齐次热传导方程 options={ 空间杆长 L, 空间点数 N , 时间点数 M,  扩散系数 alfa,  稳固条 件的值 lambda( 取值一定小于 0.5),}; topic=seting; lines=1; def={1,100,1000,1,0.5}; h=inputdlg(options,topic,lines,def); L=eval(h{1}); N=eval(h{2}); M=eval(h{3}); alfa=eval(h{4}); lambda=eval(h{5});%lambda 的值一定小于 %*************************************************** h=L/N;%空间步长 z=0:h:L; z=z; tao=lambda*h^2/alfa;% 时间步长 tm=M*tao;%热传导的总时间 tm t=0:tao:tm; t=t; %计算初值和边值 T=zeros(N+1,M+1); Ti=init_fun(z); To=border_funo(t); Te=border_fune(t); T(:,1)=Ti; T(1,:)=To; T(N+1,:)=Te; %用差分法求出温度 T与杆长 L、时间 t 的关系 for k=1:M m=2; while m=N T(m,k+1)=lambda*(T(m+1,k)+T(m-1,k))+(-2*lambda+1)*T(m,k); m=m+1; end; end; %设置立体网格 for i=1:M+1 X(:,i)=z; end; for j=1:N+1 Y(j,:)=t; end mesh(X,Y,T); view([1 -1 1]); xlabel(Z); ylabel(t); zlabel(T); function y=init_fun(z)% y=1-z.^2; return  初值条件 function y=border_funo(t)%z=0 y=1+t.*0; return  的界限条件 function y=border_fune(t)%z=L y=t*.0; return  的界限条件 运转状况： 按“ run”运转时，弹出窗口 将图中有关数据改正为： 点击图框中的“ OK”，在“ command window”中输出结果为：