这样的分岔和poincare能说明什么?

yh0247

这是一个分岔图对应的poincare映射


这四幅是另一组分岔图
谢谢

[ 本帖最后由 mjhzhjg 于 2007-6-15 13:53 编辑 ]

咕噜噜

图呢？用附件传上来

[ 本帖最后由 mjhzhjg 于 2007-6-15 13:54 编辑 ]

yh0247

[ 本帖最后由 mjhzhjg 于 2007-6-15 13:52 编辑 ]

yh0247

[ 本帖最后由 mjhzhjg 于 2007-6-15 13:52 编辑 ]

无水1324

应该是一个稳定的周期解，你要去掉不稳定的阶段，画图

[ 本帖最后由 无水1324 于 2007-6-15 10:54 编辑 ]

咕噜噜

应该是有阻尼的拟周期解，这是从图上看出来的
不过你最好还是能把你的系统微分方程贴上来如果方便的话

yh0247

无水1324

注意参数的选取，关于这种方程论坛上、书上都是以此为例的，你可以看看！程序也有很多

yh0247

是这样的，只不过我选取的参数都是固定的了，所以才会出来这样的图，但我不是搞非线性的，所以对这些周期，拟周期的搞不懂。不管怎么截取，图形基本都是这样的，总有一些离散点什么的

[ 本帖最后由 无水1324 于 2007-6-15 12:40 编辑 ]

无水1324

你把程序放上来，我帮你修改一下！

yh0247

clear all

n=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=[x0_1;x0_2];
  
    % options=odeset('RelTol',1e-3,'AbsTol',1e-4);
    [t,x]=ode45('poin',[10:T/10:t_final],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');

yh0247

如果我画关于激励P的分岔图，是不是把w（n）换成P（n）就行了，同时T=2*pi/omega （固有频率）

无水1324

对的

无水1324

你的每周期取点太少了至少也要取到200，你才取10

yh0247

其实我有些是取的很多的，出来的图形形式是一样的