一本书上的分岔图程序

如枫

下面这个是一本书上的分岔图程序，大家看看这是用什么方法画分岔图的？怎么好像没有去掉瞬态解？


% Chapter 14 - Poincare Maps and Nonautonomous Systems in the Plane.
% Programs_14g - Bifurcation diagram for the Duffing equation.
% Copyright Birkhauser 2004. Stephen Lynch.
% Make sure Programs_14f is in your directory.
clear
global Gamma;
Max=120;step=0.001;interval=Max*step;a=1;b=0;
% Ramp the amplitude up.
for n=1:Max
    Gamma=step*n;
    [t,x]=ode45('Programs_14f',[0:(2*pi/1.25):(4*pi/1.25)],[a,b]);
    a=x(2,1);
    b=x(2,2);
    rup(n)=sqrt((x(2,1))^2+(x(2,2))^2);
end
% Ramp the amplitude down.
for n=1:Max
    Gamma=interval-step*n;
    [t,x]=ode45('Programs_14f',[0:(2*pi/1.25):(4*pi/1.25)],[a,b]);
    a=x(2,1);
    b=x(2,2);
    rdown(n)=sqrt((x(2,1))^2+(x(2,2))^2);
end
hold on
rr=step:step:interval;
plot(rr,rup)
plot(interval-rr,rdown)
hold off
fsize=15;
axis([0 .12 0 2])
xlabel('\Gamma','FontSize',fsize)
ylabel('r','FontSize',fsize)
title('Bifurcation Diagram of the Duffing System. Bistable Region.')
% End of Programs_14g.

hsfy919

第一个循环就是去掉瞬态，把第一个循环计算的结果付给第二个循环，作为初值。

lixu2007a

看这种程序，要对照公式看，

飞翼

hsfy919 发表于 2013-2-5 17:48

                               
登录/注册后可看大图

第一个循环就是去掉瞬态，把第一个循环计算的结果付给第二个循环，作为初值。

那第二个循环的瞬态就不用去了？偶尔看到不太明白

hsfy919

飞翼 发表于 2013-7-9 17:13

                               
登录/注册后可看大图

那第二个循环的瞬态就不用去了？偶尔看到不太明白

第二个循环就没有瞬态了

lihaitao123

hsfy919 发表于 2013-7-12 10:54

                               
登录/注册后可看大图

第二个循环就没有瞬态了

求解释一下，他是积分时间都是从[0:(2*pi/1.25):(4*pi/1.25)]如何去掉顺瞬态的？