| 查看: 1655 | 回复: 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 |
回复此楼» 猜你喜欢
山东交通学院航运学院硕士研究生还有大量名额，欢迎工科学院报考
已经有0人回复
-
本科机器人工程，总分296，共产党员，希望调剂B区大学，考研科目085500
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有198人回复
-
一志愿重庆大学机械学硕080200 313分，求二轮调剂
已经有7人回复
-
沈阳高校招收机械、控制类等相关专业研究生（调剂）
已经有0人回复
-
陶瓷销盘实验销材料求助
已经有2人回复
-
C/C-SiC复合材料材料性能问题
已经有5人回复
-
【26年入学】红外半导体/电子信息/嵌入式视觉/诚招博士生 光电教育部重点实验室
已经有0人回复
-
【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院
已经有9人回复
4楼2018-06-06 15:07:41
luck0110
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 745.3
- 帖子: 28
- 在线: 23.7小时
- 虫号: 1992625
- 注册: 2012-09-11
- 性别: GG
- 专业: 机械测试理论与技术
3楼2014-11-27 21:03:20
20