|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
一个例子,代码如下:——用的是gghhjj提出的方法
%%%方程定义
function df=dafen(t,x,flag,force)
df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];
%%%Poincare截面绘制
clear
ff=0.222;
options=odeset('RelTol',1e-7);
tt=3*pi/2;
[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff);
plot(x(100:end,1),x(100:end,2),'.')
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')
在绘图的时候需要取点,我做了三组,一个是从100点开始取、从500点开始取,还有从100点开始取,请问各位,那个截面是最合理的,选点有没有什么要求注意的啊??
请大家一起讨论帮忙!
还有一个问题——用这样的方法作出来的截面图感觉就是将相图进行点离散,那向很多书上做出的那些很漂亮的Poincare截面是如何作出来的呢?还请大家指教啊! |
-
这个是从第100点开始取点
-
这个是从第500点开始取点
-
这个是从第1000点开始取点
|