| 查看: 3662 | 回复: 2 | ||
chnzhn木虫 (小有名气)
|
[求助]
ODE45能解非线性方程组吗？
|
|
今有一非线性微分方程组，欲求解，下列做法怎么求解不了，请教各位是为何？ function Ad_gene [t,y]=ode45(@hyperliu2,[0:0.01:10],[20;20;30;20;1;1]); figure; plot(t,y); function dy=hyperliu2(t,y) dy=zeros(6,1); dy(1)=y(5)*(y(2)-y(1))+y(2)*y(3); dy(2)=25*y(1)-y(1)*y(3)-y(2)-y(4); dy(3)=y(1)*y(2)-8/3*y(3); dy(4)=y(6)*y(1)+0.5*y(2)*y(3)+y(4); dy(5)=-y(1)+1; dy(6)=-2*y(4); 以上编程不知错在哪里，运行很久都不出结果，是否进入死循环呢？ |
回复此楼» 本帖已获得的红花（最新10朵）
chnzhn» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有26人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解非线性方程组
已经有16人回复
- 解线性代数方程组
已经有10人回复
- 【课件】非线性方程组课件（华东理工大学）
已经有26人回复
- 一个微分方程组，求解
已经有18人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
chnzhn: 金币+10, ★★★很有帮助, 你说的很好，与初始值、步长都有关，谢谢了 2012-04-06 09:35:26
xiegangmai: 金币+1, 谢谢应助！ 2012-04-06 10:13:44
感谢参与，应助指数 +1
chnzhn: 金币+10, ★★★很有帮助, 你说的很好，与初始值、步长都有关，谢谢了 2012-04-06 09:35:26
xiegangmai: 金币+1, 谢谢应助！ 2012-04-06 10:13:44
|
1. ode45可以解非线性方程组。 2. 计算时间很长有多长呢？ [t,y]=ode45(@hyperliu2,[0:0.01:10],[20;20;30;20;1;1]); 这句改成 [t,y]=ode45(@hyperliu2,[0,10],[20;20;30;20;1;1]); 试试。因为ode45默认自动步长的，会比你指定的步长更有效率。 3. 计算慢的原因可能是方程存在刚性，试试换成ode15s求解。 |
2楼2012-04-06 09:28:41
chnzhn
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3409.1
- 散金: 10
- 帖子: 163
- 在线: 108.9小时
- 虫号: 162699
- 注册: 2006-01-10
- 性别: GG
- 专业: 机械动力学
3楼2012-04-06 09:33:17
50