|

楼主 |
发表于 2011-2-25 19:49
|
显示全部楼层
回复 2 # hsfy919 的帖子
主任,你帮我看看我的庞加莱程序有没有问题?
clear;
options=odeset('RelTol',1e-6,'AbsTol',1e-6);
f=2.9;
c=160
tt=1/f;
y0=[0,0,0,0,0,0,0,0];
[t,y]=ode45('biyelunwen21',[0:tt/100:1500*tt],y0,options,c) ;
plot(y(140000:100:end,1),y(140000:100:end,2),'r.')
hold on
plot(y(140000:1:end,1),y(140000:1:end,2))
激励频率是2.9,这样的方法做庞加莱图合适不?主要是这一块,[0:tt/100:1500*tt]
我一个周期取100个点,间隔100个点取点画庞加莱
我现在担心我的程序是不是错了,是不是对我这个系统部合适。图出的这种结果有点矛盾啊 |
|