| 查看: 1540 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
如何绘制庞加莱截面图 已有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 |
回复此楼» 猜你喜欢
机械专硕311求调剂
已经有0人回复
-
一志愿武汉理工大学085502车辆工程，英二数二，求调剂
已经有0人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有84人回复
-
15 年非标老工程师忠告：转台选型，这 3 个参数错一步全白搭
已经有1人回复
-
294 求调剂
已经有2人回复
-
讲师身份申请面上项目中标率有影响吗
已经有32人回复
-
一志愿西北工业大学机械学硕080200 367分求调剂
已经有2人回复
-
希望面上好运
已经有66人回复
-
车辆调剂申请
已经有1人回复
-
调剂求助一志愿211初试295
已经有1人回复
3楼2014-11-27 21:03:20
4楼2018-06-06 15:07:41
