| 查看: 593 | 回复: 0 | ||
apoorone铁虫 (小有名气)
|
[求助]
求助matlab
|
|
https://wenku.baidu.com/view/95c28019ff00bed5b9f31d18.html 这是我下载的一维热传导方程，我想将这个方程改成分段函数，即当0<=l<=0.5时，alfa=1，当0.5<=L<=1,alfa=2,想请问怎么修改 我是菜鸟，刚开始学，还是不知怎么改，希望大神帮忙 这是，他的程序 有限差分法求出温度T与杆长L、时间t的关系 求解下列热传导问题： 程序： 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的值必须小于0.5 %*************************************************** 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 [ Last edited by apoorone on 2017-5-9 at 14:21 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有31人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有72人回复
-
FileCode能看出啥？
已经有29人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
