| 查看: 1514 | 回复: 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 |
回复此楼» 猜你喜欢
英国巴斯大学 University of Bath 博士后及博士招聘
已经有5人回复
-
求调剂收留
已经有1人回复
-
机械工程论文润色/翻译怎么收费?
已经有243人回复
-
085500机械专硕 331分求调剂
已经有1人回复
-
380，求调剂
已经有1人回复
-
机械专硕311求调剂
已经有0人回复
-
学硕080200 315分 求调剂 2026年安徽省优秀毕业生
已经有0人回复
-
机械专硕求调剂
已经有0人回复
 |
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
