| 查看: 3460 | 回复: 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朵）
科大贝贝» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有101人回复
-
时间戳又变了8-15
已经有26人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 线性常微分方程求解
已经有5人回复
- matlab解微分方程
已经有10人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】如何用matlab求解这个微分方程？【已解决】
已经有3人回复
zhaohu_727
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7116.9
- 红花: 2
- 帖子: 776
- 在线: 495.7小时
- 虫号: 809332
- 注册: 2009-07-15
- 专业: 资源化工
9楼2012-04-27 10:49:26
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
