| 查看: 609 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】matlab常微分方程组求解
|
|||
|
用以下程序进行计算，数据结果却是一正一负交替， 但是方程预期的解是不可能存在负值的，思考过可能存在的问题，初值，方程组都没问题，请问matlab自身的解法有可能造成这种结果吗？如果有，可以告诉我怎么解决吗? -------------------------------------------------------------------------------------- 部分程序： Ki=108.85;Kp=108.85;Ktm=0.012;Ktal=0.12;Kth=4.8;Kd=0.001; M=0.21；H=0.002 y0=zero(1,2*n+2) y0(1,1)=0.0000045 y0(1,2)=0.00045 [t,y]=ode45(@mwd,tspan,y0) function dydt=mwd(t,y) global Ki Kp Ktm Ktal Kth Kd M H n ju=0; for j=3: (n+2) ju=ju+y(j); end mwd(1)=(Kth*H^0.5+Ktal*y(2)+Ktm*M+Kd)*ju-Ki*y(1)*M-Kd*y(1); mwd(2)=-Ktal*y(2)*ju; mwd(3)=Ki*y(1)*M-Kp*y(3)*M-(Kth*H^0.5+Ktal*y(2)+Ktm*M+Kd)*y(3); for k=4: (n+2) mwd(k)=Kp*y(k-1)*M-Kp*y(k)*M-(Kth*H^0.5+Ktal*y(2)+Ktm*M+Kd)*y(k); end for k=3: (n+2) mwd(k+n)=(Kth*H^0.5+Ktal*y(2)+Ktm*M+Kd)*y(k); end dydt=mwd( : ); [ Last edited by gg0520 on 2009-9-18 at 15:37 ] |
回复此楼» 猜你喜欢
实验室接单子
已经有6人回复
-
假如你的研究生提出不合理要求
已经有11人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
backstroke9680
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 578.4
- 散金: 1800
- 帖子: 374
- 在线: 100.7小时
- 虫号: 720194
- 注册: 2009-03-11
- 性别: GG
- 专业: 无机非金属类光电信息与功
3楼2009-09-28 19:34:57
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
把问题说清楚一些！ 2楼2009-09-27 22:13:18
4楼2009-09-28 21:21:01
sui2066
木虫 (职业作家)
- 应助: 5 (幼儿园)
- 金币: 6188.4
- 散金: 8
- 红花: 3
- 帖子: 3600
- 在线: 809小时
- 虫号: 382693
- 注册: 2007-05-25
- 性别: GG
- 专业: 聚合物共混与复合材料
5楼2009-09-28 21:22:14