请看看这个poincare图怎么回事？

wangu

我刚接触非线性不久，想把刘延柱的《非线性振动》中的p255页对应的duffing方程的分岔图，相图和pioncare图画出来。当激振力f为0.28时应为P-4运动，但我画出来的确是这样的

wangu

回复 1 # wangu 的帖子

对应的程序
function dx=duffing(t,X)
global w f;
c=0.3;
x=X(1);
y=X(2);
psi=X(3);
dx=zeros(3,1);
dx(1)=y;
dx(2)=-c*y+x-x^3+f*cos(psi);
dx(3)=w;


主程序
function duffing_poincare_f
clear;
global f w;
w=1.2;
T=2*pi/w;
x0=[0,0.001,0];
f=0.28;
[t,x]=ode45('duffing',[0:T/10:10000*T],x0) ;
plot(x(7000:10:end,1),x(7000:10:end,2),'k.','markersize',2);
xlabel('x');ylabel('xdot');title('poincare');

wangu

回复 1 # wangu 的帖子

先问一下，我的程序对不？

gghhjj

看不懂你的程序，我简单的写了一个
function df=dafen(t,x,force)
df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];


clc;
clear;
force=0.28;
options=odeset('RelTol',1e-7);
[t,x]=ode45(@dafen,[0:0.01:500],[0,0],options,force);
plot(x(10000:end,1),x(10000:end,2))

wangu

回复 4 # gghhjj 的帖子

谢谢，大概明白点了，跟频闪法的分岔图程序类似，但是有个问题，这个时间的积分区间是自己随便选的吗？[0:0.01:500]

gghhjj

wangu 发表于 2012-7-19 19:30

                               
登录/注册后可看大图

回复 4 # gghhjj 的帖子

谢谢，大概明白点了，跟频闪法的分岔图程序类似，但是有个问题，这个时间的积分区 ...

是的，看你需要