| 查看: 1650 | 回复: 3 | ||
luck0110金虫 (初入文坛)
|
[求助]
如何绘制庞加莱截面图 已有1人参与
|
|
照着别人的思路，绘制了达芬方程的仿真图，但是感觉不对； 程序如下： 首先建立了存放函数的M函数； function dx=mm_4(t,x) dx=[x(2);cos(t)-x(1)^3+x(1)-4*x(2)]; 然后又编写了画制庞加莱截面图的M函数。 function f tt=2*pi/1.2; [t,x]=ode45(@mm_4,[0:tt/100:80*tt],[1,1]); figure subplot(2,1,1) plot(t,x(:,1)); title('位移曲线'); xlabel('t');ylabel('x'); subplot(2,2,3) plot(x(2000:end,1),x(2000:end,2),'-') xlabel('x');ylabel('dx/dt'); title('相图,r=4,w=1,f=1'); subplot(2,2,4) i=2000:100:3000 ; plot(x(i,1),x(i,2),'*') title('庞加莱截面'); 通过修改r的值，得到不同的图形； 通过图片可以看出，r=4，相图为一个封闭图形，庞加莱界面应该是一个点才对；不知道错在哪个地方？ 拜托高手帮忙了~谢谢。  1.JPG  2.JPG |
回复此楼» 猜你喜欢
测面比电阻，600-850℃
已经有0人回复
-
北京航空航天大学航空发动机研究院2026年暑期学校通知
已经有2人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有108人回复
-
Progress in Materials Science 期刊投稿求助
已经有2人回复
-
【征稿通知】2027 IPPTC 国际石油石化技术会议征稿启动
已经有0人回复
-
含铝高熵合金 真空感应熔炼
已经有0人回复
-
材料博士申请
已经有21人回复
 |
2楼2014-11-27 21:02:53
luck0110
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 745.3
- 帖子: 28
- 在线: 23.7小时
- 虫号: 1992625
- 注册: 2012-09-11
- 性别: GG
- 专业: 机械测试理论与技术
3楼2014-11-27 21:03:20
4楼2018-06-06 15:07:41
100