24小时热门版块排行榜    

查看: 905  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cptll

银虫 (小有名气)

[求助] 偏微分方程组的求解 已有1人参与

在做毕业设计过程中遇到偏微分方程的求解问题,问问各位能否用Matlab解决。
方程组为:
(1)dN/dl=F*A*exp(-E/RT)*N
  (2)   dQ=cP*N*dT+h*dN;
   (3)dQ/dl=3.14a*D*(Tw-T);
其中Q、N、l、T是变量,其他已知,能否得到T随l的变化关系呢?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cptll

银虫 (小有名气)

引用回帖:
2楼: Originally posted by arising2010 at 2014-05-07 13:08:14
先算dQ/dl,先根据(2)可得:
dQ/dl= (cP*N*dT+h*dN)/dl=cp*N*dT/dl + h*dN/dl

  再将上式中把(1)的dN/dl 带入
dQ/dl=cp*N*dT/dl + h* F*A*exp(-E/RT)*N

所以  cp*N*dT/dl + h* F*A*exp(-E/RT)*N =3.14a ...

谢谢你~ 用你的方法计算了下,是没有解的,个人认为这是属于偏微分方程了,能否将偏微分方程和常微分方程联立求解呢?
3楼2014-05-07 19:24:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

arising2010

铜虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
cptll: 金币+5, ★★★★★最佳答案 2014-05-12 13:58:01
先算dQ/dl,先根据(2)可得:
dQ/dl= (cP*N*dT+h*dN)/dl=cp*N*dT/dl + h*dN/dl

  再将上式中把(1)的dN/dl 带入
dQ/dl=cp*N*dT/dl + h* F*A*exp(-E/RT)*N

所以  cp*N*dT/dl + h* F*A*exp(-E/RT)*N =3.14a*D*(Tw-T)

即   dT/dl =3.14a*D*(Tw-T)/(cp*N) -h* F*A*exp(-E/RT)/cp(如果N不为0的话,如果N可能为0,需要再讨论)

再联立 (1) : dN/dl=F*A*exp(-E/RT)*N

二元常微分方程组,可以直接利用matlab等软件解决
还需要初值:N和T在时刻为0时的初值
2楼2014-05-07 13:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arising2010

铜虫 (正式写手)

引用回帖:
3楼: Originally posted by cptll at 2014-05-07 19:24:16
谢谢你~ 用你的方法计算了下,是没有解的,个人认为这是属于偏微分方程了,能否将偏微分方程和常微分方程联立求解呢?...

没有解的原因很多:例如模型有问题的可能性很大。你模型的参数和初值设置不当,软件没有合理应用,都可能导致无解

你认为这是偏微分方程组?其实质就是个常微分方程组吧

[ 发自手机版 http://muchong.com/3g ]
4楼2014-05-07 22:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见