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

下载本文档

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

文档简介

1、MATLAB编辑一维热传 导方程的模拟程序 -CAL-FENGHAI-(2020YEAR-YICAI)INGBIAN 求解下列热传导问题: -1=0(。5厶) a dt r(z,o)=i-z2 r(07)=i, t(lj)=o L = l, a = 1 程序： function heat_conduction() %一维齐次热传导方程 options彳空间杆长L：空间点数NJ时间点数扩散系数alfa稳定条件的值 lambdaf取值必须小于 topic=,seting,; lin es=l; def=l,；100,z,1000,z,l,/; h=in putdlg(opti on s,topic

2、,li nesdef); L=eval(hl); N=eval(h2); M=eval(h3); alfa=eval(h4); lambda=eval(h5);%lambda 的值必须小于 h=L/N;%空间步长 z=O:h:L; z=z； tao=lambda*hA2/alfa;% 时间步长 tm=M*tao;%热传导的总时间tm t=0:tao:tm; t=f; %计算初值和边值 T=zeros(N+l,M+l); Ti=init_fu n(z); To=border_funo(t); Te=border_fu ne(t); T(:,l)=Ti; T(l,：)=To; T(N+l,:)=T

3、e; %用差分法求出温度T与杆长L、时间t的关系 for k=l:M m=2; while m=N T(m,k+l)=lambda* (T(m+l/k)+T(m-l/k)+(-2*lambda+l)*T(m,k); m=m+l; end; end; %设置立体网格 for i=l:M+l X(：,i)=z； end; for j=l:N+l Y(j,：)=t； end mesh(XYT); view(l -1 1); xlabel(Z); ylabelCf); zlabel(T); function y=init_fun(z)%初值条件 y=l-z.A2; return function y=border_funo(t)%z=0的边界条件 y=l+t.*O; return function y=border_fune(t)%z=L的边界条件 y=t*.O; return 运行情况: 按run运行时，弹出

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论