MATLAB编辑一维热传导方程的模拟程序.docxVIP

求解下列热传导问题: 2T 1 T 门 c , 2 0 0 z L z t T乙0 1 z2 T 0,t 1, T L,t 0 L 1, 1 程序： fun ctio n heat_c on ductio n() % 一维齐次热传导方程 稳定条options={空间杆长L,空间点数N,时间点数M,扩散系数alfa, 件的值lambda(取值必须小于0.5),}; 稳定条 topic=set in g; lin es=1; def={1,100,1000,1,0.5}; h=in putdlg(opti on s,topic, lin es,def); L=eval(h{1}); N=eval(h{2}); M=eval(h{3}); alfa=eval(h{4}); lambda=eval(h{5});%lambda 的值必须小于 0.5 o%*************************************************** h=L/N;%空间步长 z=0:h:L; z=z; tao=lambda*hA2/alfa;% 时间步长 tm=M*tao;%热传导的总时间tm t=0:tao:tm; t=t; %计算初值和边值 T=zeros(N+1,M+1); Ti=i nit_fu n(z); To=border_fu no (t); Te=border_fu ne