wangu 发表于 2012-7-16 16:40

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

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

wangu 发表于 2012-7-16 16:41

回复 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=;
f=0.28;
=ode45('duffing',,x0) ;
plot(x(7000:10:end,1),x(7000:10:end,2),'k.','markersize',2);
xlabel('x');ylabel('xdot');title('poincare');

wangu 发表于 2012-7-17 20:33

回复 1 # wangu 的帖子

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

gghhjj 发表于 2012-7-19 16:15

本帖最后由 gghhjj 于 2012-7-19 16:16 编辑

看不懂你的程序，我简单的写了一个
function df=dafen(t,x,force)
df=;


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

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

回复 4 # gghhjj 的帖子

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

gghhjj 发表于 2012-7-20 10:07

wangu 发表于 2012-7-19 19:30 static/image/common/back.gif
回复 4 # gghhjj 的帖子

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

是的，看你需要
页: [1]
查看完整版本: 请看看这个poincare图怎么回事？