| 查看: 3613 | 回复: 17 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
高手进来用matlab求解常微分方程，十分感谢！
|
||
|
本人数学建模型后，得到一常微分方程，相对较复杂，手算不出来，听说matlab很强大，据说能解常微分方程，但本人是一只小菜鸟啊，小菜鸟，所以把本人的常微分方程写出来让大家看看，高手们帮帮忙。成功求解后还有重赏。。。提前感谢下 dy/dx=0.002*A*y1/y3+0.002*B*y*y2/y3; (A,B是两个常数) y1=4957.64*ln[1.04*(2.529-y*y3)/(0.5*y-y*y3); y2=89.24*(2.749-0.473*y-0.145*y*y3)/(0.5-y3); y3=0.0000006*x^2-0.0007803*x; |
回复此楼» 本帖已获得的红花（最新10朵）
科大贝贝» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
17楼2013-05-09 09:03:37
2楼2012-04-26 21:09:06
科大贝贝
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3868.8
- 散金: 5
- 红花: 1
- 帖子: 302
- 在线: 152.6小时
- 虫号: 825964
- 注册: 2009-08-11
- 性别: GG
- 专业: 环境化工
3楼2012-04-26 21:10:32
zhaohu_727
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7116.9
- 红花: 2
- 帖子: 776
- 在线: 495.7小时
- 虫号: 809332
- 注册: 2009-07-15
- 专业: 资源化工
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
soliton923: 金币+1, 谢谢参与讨论~~~ 2012-04-26 22:35:03
傻妞丫头: 金币+5, ★有帮助, 输入程序里面好像运行不出来，有错误，而且y1中的对数是log以e为底的，我再看看吧，先给你5个金币啊。如果运行出来，肯定全给金币给你的 2012-04-27 09:02:49
傻妞丫头: 金币+10, ★★★很有帮助, 高手，不知道咋回事，这个程序，我这边运行不了，你能不能再帮忙看看？添加初始条件y(0)=0。拜托拜托了。最好能将程序写好了给我。谢谢 2012-04-27 14:36:19
感谢参与，应助指数 +1
soliton923: 金币+1, 谢谢参与讨论~~~ 2012-04-26 22:35:03
傻妞丫头: 金币+5, ★有帮助, 输入程序里面好像运行不出来，有错误，而且y1中的对数是log以e为底的，我再看看吧，先给你5个金币啊。如果运行出来，肯定全给金币给你的 2012-04-27 09:02:49
傻妞丫头: 金币+10, ★★★很有帮助, 高手，不知道咋回事，这个程序，我这边运行不了，你能不能再帮忙看看？添加初始条件y(0)=0。拜托拜托了。最好能将程序写好了给我。谢谢 2012-04-27 14:36:19
|
syms y1 y2 y3 y x A B; y3=0.0000006*x^2-0.0007803*x; y1=4957.64*log10(1.04*(2.529-y*y3)/(0.5*y-y*y3)); y2=89.24*(2.749-0.473*y-0.145*y*y3)/(0.5-y3); >> y=dsolve(['Dx-(0.002*A*y1/y3+0.002*B*y*y2/y3)=0']) y = C2 + (t*(A*y1 + B*y*y2))/(500*y3) |
4楼2012-04-26 21:59:40
