这样的分岔和poincare能说明什么?
这是一个分岔图对应的poincare映射
这四幅是另一组分岔图
谢谢
[ 本帖最后由 mjhzhjg 于 2007-6-15 13:53 编辑 ] 图呢？用附件传上来
[ 本帖最后由 mjhzhjg 于 2007-6-15 13:54 编辑 ]
不好意思，不会贴图
[ 本帖最后由 mjhzhjg 于 2007-6-15 13:52 编辑 ]
第二组图
[ 本帖最后由 mjhzhjg 于 2007-6-15 13:52 编辑 ] 应该是一个稳定的周期解，你要去掉不稳定的阶段，画图
[ 本帖最后由 无水1324 于 2007-6-15 10:54 编辑 ] 应该是有阻尼的拟周期解，这是从图上看出来的
不过你最好还是能把你的系统微分方程贴上来如果方便的话
方程
注意参数的选取，关于这种方程论坛上、书上都是以此为例的，你可以看看！程序也有很多回复 #8 无水1324 的帖子
是这样的，只不过我选取的参数都是固定的了，所以才会出来这样的图，但我不是搞非线性的，所以对这些周期，拟周期的搞不懂。不管怎么截取，图形基本都是这样的，总有一些离散点什么的[ 本帖最后由 无水1324 于 2007-6-15 12:40 编辑 ]
回复 #9 yh0247 的帖子
你把程序放上来，我帮你修改一下！ clear alln=0.0005;
ro=2600;
E0=2.3*10^5;
E1=5.3*10^4;
E2=-1.6*10^4;
A=20*10^(-6);
pe=0.005;
l=0.1;
P=0.5;
T=10;
mu=n/ro/A;
alpha=(1-pe*T)*E1*pi^2/ro/l^3/3;
beta=3*(1-pe*T)*E2*pi^4/ro/l^4/64;
omega=sqrt((1-pe*T)*E0*pi^2/ro/l^2/4);
F=2*P/ro/A;
%w1=1/3*omega-0.2;
x0_1=0.00001;x0_2=0.00001;
w=0.3:0.005:2;
hold on
for n=1:length(w)
T=2*pi/w(n);
t_final=500*T;
x0=;
% options=odeset('RelTol',1e-3,'AbsTol',1e-4);
=ode45('poin',,x0,[],w(n));
plot(w(n),x(10:10:end,1),'k.','markersize',5),hold on;
end
grid
%xlabel('轴向力 Excitation P/N');%,'fontsize',14);
ylabel('位移 Displacement a/m');
%ylabel('速度 Speed V/m/s');
xlabel('频率 Frequency \Omega/rad/s');
%ylabel(' 位移 Displacement a/m');