| 查看: 871 | 回复: 3 | |||
[交流]
【求助】急问，各位师哥师姐帮小弟检查一下matlab分岔程序吧！！小弟万分感谢！！！ 已有2人参与
|
|
我想用Matlab软件得到下面微分方程组的分岔图，哪位高手能告诉我分岔程序如何写才能画出分岔图，谢谢！ 我的公式：（公式粘不上，在附件中也有） function dx=fun5(t,x,a) dx(1)=0.000088+((0.051253+1.74*x(1)+2981*x(1)^2)/(1+157*x(1)+8.8*x(1)^2+12.67*x(3)+0.05*x(3)^2+48.5*x(1)*x(3)))-0.198*x(1); dx(2)=((0.014+259*x(1)+172*x(1)^2)/(1+157*x(1)+49*x(1)^2+25.133867*x(3)+0.0134*x(3)^2+1.17*x(1)*x(3)))-0.075*x(2); dx(3)=0.324+((264*x(1)+763*x(1)^2+1.89*x(2)+358*x(2)^2)/(1+x(1)+0.055*x(1)^2+0.224*x(2)+0.3*x(2)^2+8.4*x(1)*x(2)))+0.133*x(4)+0.043*x(5)-1.538*x(3); dx(4)=a*x(3)-0.1332*x(4); dx(5)=0.0913*x(3)-0.04365*x(5); dx=[dx(1);dx(2);dx(3);dx(4);dx(5)];即如附件所示，其中a 的变化范围 0.99~1.077 时要画出x(4)和a 之间的分岔图，初始值 x0=[0.01;0.05;12.5;72;19]。 以下我的程序： a=0.99:0.001:1.075; hold on for j=1:length(a) [t,x]=ode45('fun5',[0,8000],[0.01;0.05;12.5;72;19],[],a(j)); plot(a(j),x(:,4)) hold on end 可是运行不出来，小弟是新手，不知道错了哪里，希望各位高手帮我一下，小弟真的是急死了！！！！  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大家有没有那种idea理论推导挺合理，但matlab仿真仿不出来的情况啊？
已经有6人回复
- 希望明年公派德国，现在研二，求各位师兄师姐指导！
已经有13人回复
- 求助格氏反应，小弟感谢！！！
已经有9人回复
- 德国公派携签，资金证明急问
已经有10人回复
- 小弟初学生物信息学，完全摸不着头脑，求启蒙
已经有7人回复
- 关于XPS的几个问题，小弟先再次谢过了
已经有8人回复
- 求助申请文部省奖学金的师兄师姐们帮小弟解解惑呀
已经有11人回复
- 关于XPS的几个问题，求指导！（小弟金币实在寒颤，真心希望高手解答）
已经有12人回复
- 小弟新手，谁帮我分析一下这个交流阻抗图啊
已经有8人回复
- 中科院的师哥师姐、、、帮忙给讲解一下呗
已经有7人回复
- 小弟想请教一个Origin作图问题
已经有16人回复
- 哪位师哥师姐能告诉小弟NH3BH3（硼烷铵）在哪买啊
已经有4人回复
- 【求助】各位师哥师姐帮小弟检查一下复杂参数回归的程序吧！！小弟万分感谢！！！
已经有3人回复
xunyuqi
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2209.9
- 散金: 68
- 红花: 3
- 帖子: 304
- 在线: 242.8小时
- 虫号: 487564
- 注册: 2007-12-28
- 性别: GG
- 专业: 管理复杂性研究
我也同求此问题
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
据说有奇点的话需要忽略，否则matlab程序运行不过去，我也正愁如何画呢，如果解决了也请告知我，多谢 ，另将我的问题贴上http://blog.sina.com.cn/s/blog_67454ac70100t3kh.html |
2楼2011-04-17 00:32:58
lzwhw2000
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1143.1
- 散金: 100
- 帖子: 21
- 在线: 32.3小时
- 虫号: 2008082
- 注册: 2012-09-17
- 性别: GG
- 专业: 常微分方程与动力系统
3楼2013-02-02 23:07:26
lzwhw2000
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1143.1
- 散金: 100
- 帖子: 21
- 在线: 32.3小时
- 虫号: 2008082
- 注册: 2012-09-17
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2013-02-02 23:09:31
19