| 查看: 1651 | 回复: 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 |
回复此楼» 猜你喜欢
2026年国自然面上资助率
已经有24人回复
-
影响面上的因素
已经有11人回复
-
金属材料论文润色/翻译怎么收费?
已经有240人回复
-
filecode
已经有11人回复
-
给基金委的建议
已经有1人回复
-
听说今天filecode变了
已经有47人回复
-
filecode变化情况
已经有22人回复
-
filecode与中标关系的预测
已经有3人回复
-
关于filecode
已经有7人回复
-
好奇怪的filecode
已经有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
100