| 查看: 5033 | 回复: 1 | |||
zhaohuiyuan金虫 (著名写手)
|
[求助]
一维热传导偏微分方程的数值解法求助
|
|
各位虫友们，你们好！由于我对编程软件运用的不是很熟悉，希望得到各位的帮助，我编出来的程序怎么都运行不出来 不知道是为什么，希望各位高手能指点一下，我要求解的方程是（见图片），如果哪位能帮上我，让我的程序能运行出结果来的话，我将用所有的金币来酬谢，若能给出有用意见的也行，我将酌情给金币，谢谢！ 程序如下 function heat_conduction() %一维齐次热传导方程 options={'空间半径ro','空间点数N' ,'时间点数M','密度p','比热c','导热系数s',‘'对流传热系数h',}; topic='seting'; lines=1; def={'0.032','100','1000','1714.8','1.48','0.535','61.27'}; b=inputdlg(options,topic,lines,def);%创建输入对话 ro=eval(b{1});%eval执行字符串 N=eval(b{2}); M=eval(b{3}); p=eval(b{4}); c=eval(b{5}); s=eval(b{6}); h=eval(b{7}); %*************************************************** Fo=s/(p*c); or=ro/N;%空间步长 r=0  r:ro;r=r'; ot=Fo*or^2/a;%时间步长 tm=M*ot;%热传导的总时间tm t=0 t:tm;t=t'; %计算初值和边值 T=zeros(N+1,M+1); Ti=init_fun(r);%fun为目标函数的表达式字符串或MATLAB自定义函数的函数柄 To=border_funo(t); Te=border_fune(t); T(:,1)=Ti;%t=0时的初始温度 T(1,  =To;%r=0时的温度值T(N+1, =Te;%r=r0时的温度%用差分法求出温度T与半径r、时间t的关系 for k=1:M m=2; while m<=N&&Fo<=(ro/(2*r+or)) T(m,k+1)=(Fo+(or*Fo)/ro)*(T(m+1,k)+(1-(or*Fo)/r+2*Fo)*T(m,k))+Fo*T(m-1,k); m=m+1; end; end; %设置立体网格 for i=1:M+1 X(:,i)=r; end; for j=1:N+1 Y(j, =t;end mesh(X,Y,T); view([1 -1 1]); xlabel('r'); ylabel('t'); zlabel('T'); function y=init_fun(r)%初值条件t=0 y=14+t.*0 return function y=border_funo(t)%r=0的边界条件 for k=1:M m=1; while Fo<=0.5 T(m+1,k+1)=(1-2Fo)*T(m+1,k)+2Fo*T(m,k); end; end; return function y=border_fune(t)%r=ro的边界条件 Tf=40+t.*0; A=(h*ot)/(p*c*or); for k=1:M m=N; while (1-A-Fo)>=0 T(m+1,k+1)=A*Tf+(1-A-Fo)*T(m+1,k))+Fo*T(m,k); end; end; return  需求解的方程 |
回复此楼» 收录本帖的淘帖专辑推荐
 好帖子 |
» 本帖已获得的红花（最新10朵）
zhaohuiyuan» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 跟大家请教一个一维热传导方程建立的问题（高分求助）
已经有3人回复
- matlab 求解三元方程组
已经有9人回复
- 求助，计算方程最大值
已经有5人回复
- 中心差分求解二维热传导的matlab程序
已经有12人回复
- erfc的求解问题。谢谢。
已经有5人回复
- 求助-帮我用mathematics计算下下面的方程组
已经有13人回复
- 问个热传导问题的计算
已经有15人回复
- 求助matlab编程：一维热传导问题
已经有5人回复
- 数值求解N-S方程一般用什么编程啊？？？
已经有9人回复
- 求解，热传递方程
已经有6人回复
- 利用有限元方法求解Maxwell方程组、流体力学方程以及热传导方程的软件
已经有10人回复
- 哪位高手能用matlab编程求解一维瞬态热传导方程，先谢谢啦
已经有7人回复
- 【求助】紧束缚模型与薛定谔方程的区别与联系
已经有13人回复
- 【资源】数值解 光纤中 非线性薛定谔方程 软件
已经有27人回复
- 【求助】一道博士生面试题：薛定谔方程和热传导方程
已经有10人回复
- 【求助】解含时薛定谔方程的问题
已经有5人回复
- 【求助】O3的薛定谔方程表达式，及每一项的物理意义
已经有4人回复
zhaohuiyuan
金虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 1368.9
- 散金: 3625
- 红花: 21
- 帖子: 1135
- 在线: 412.2小时
- 虫号: 2005943
- 注册: 2012-09-16
- 性别: MM
- 专业: 工程热物理相关交叉领域
2楼2014-03-06 18:06:57
