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

无水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=432;
hold on
for n=1:length(w)
    T=2*pi/w(n);
    t_final=400*T;
    x0=[x0_1;x0_2];

    [t,x]=ode45('poin',[1:T/200:t_final],x0,[],w(n));
    %plot(w(n),x(1000:1:end,1),'k.','markersize',5),hold on;
end
grid
plot(x(1000:1:end,1),x(1000:1:end,2),'k.','markersize',5),hold on;

xlabel('位移 Displacement a/m');
ylabel('速度 Speed V/m/s');

%*****************************************
%poin 函数
%**************
function xdot=poin(t,x,flag,w)
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;
xdot=[x(2);
    -mu*x(2)-beta*(x(1))^3-omega^2*x(1)-alpha*(x(1))^2+F*cos(w*t)];

无水1324

改为：plot(w(n),x(1000:200:end,1),'k.','markersize',5),hold on

200个点是一个周期，一个周期取一个点，所以是x(1000:200:end,1),

无水1324

w=1:432;% w要变化，当画分叉图的时候
hold on
for n=1:length(w)
    T=2*pi/w(n);
    t_final=400*T;
    x0=[x0_1;x0_2];
    [t,x]=ode45('poin',[1:T/200:t_final],x0,[],w(n));
    plot(w(n),x(70000:200:end,1),'k.','markersize',5),hold on;
end

无水1324

仔细看一下别人程序是怎么用的

yh0247

:loveliness: 以前自己的疑惑也解开了，谢谢无水：）真是专业阿！！

无水1324

解决了，以后多回来看看，多帮助其他人，自己也可以多学习点东西！

yh0247

恩，我正在重画，等画好了再请教。
我经常来，只是大部分问题我确实不知道
所以都是看大家的回答，一知半解
才导致一直画不出图来

sssssxxxxx921

你怎么取了４００个周期啊，你取周期时试算了吗　　还是有经验啊

无水1324

这个可以试算一下，会减少不少时间，但是你试算的点最好分散点

gghhjj

原帖由 sssssxxxxx921 于 2007-6-22 17:37 发表

                               
登录/注册后可看大图

你怎么取了４００个周期啊，你取周期时试算了吗　　还是有经验啊

这个显然是试算的，简单算几个频率下的就可以了

无水1324

综合考虑各试算点，取稳定时间比较长的哪个