MATLAB编辑一维热传导方程的模拟程序_第1页
MATLAB编辑一维热传导方程的模拟程序_第2页
MATLAB编辑一维热传导方程的模拟程序_第3页
MATLAB编辑一维热传导方程的模拟程序_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、求解下列热传导问题:Ft1s!:-；：t=02*T(z,0 )=1 _zT(0,t )=1, T(L,t )=0L =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',

2、'1','0.5'h=in putdlg(opti on s,topic, lin es,def);L=eval(h1);N=eval(h2);M=eval(h3);alfa=eval(h4);lambda=eval(h5);%lambda 的值必须小于 0.5o%*h=L/N;%空间步长z=0:h:L;z=z'tao=lambda*hA2/alfa;% 时间步长tm=M*tao;%热传导的总时间tmt=0:tao:tm;t=t'%计算初值和边值T=zeros(N+1,M+1);Ti=i nit_fu n(z);To=border_fu no (

3、t);Te=border_fu ne(t);T(:,1)=Ti;T(1,:)=To;T(N+1,:)=Te;%用差分法求出温度T与杆长L、时间t的关系for k=1:Mm=2;while m<=NT(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+1X(:,i)=z;end;for j=1:N+1Y(j,:)=t;endmesh(X,Y,T);view(1 -1 1);xlabel('Z');ylabel('t');zlabel('T');fun ctio n y=i nit_fun(z) %初值条件y=1-z.A2;returnfun ctio n y=border_fu no(t)%z=0的边界条件y=1+t.*0;returnfun cti on y=border_f un e(t)%z=L的边界条件y=t*.0;return运行情况:按“run”运行时，弹出窗口seting穏定条件的值滋砒加（取值必须小于05OK C

人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论