如枫 发表于 2013-2-1 20:38

一本书上的分岔图程序

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


% 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;
    =ode45('Programs_14f',,);
    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;
    =ode45('Programs_14f',,);
    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()
xlabel('\Gamma','FontSize',fsize)
ylabel('r','FontSize',fsize)
title('Bifurcation Diagram of the Duffing System. Bistable Region.')
% End of Programs_14g.

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

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

lixu2007a 发表于 2013-2-12 23:42

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

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

hsfy919 发表于 2013-2-5 17:48 static/image/common/back.gif
第一个循环就是去掉瞬态，把第一个循环计算的结果付给第二个循环，作为初值。

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

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

飞翼 发表于 2013-7-9 17:13 static/image/common/back.gif
那第二个循环的瞬态就不用去了？偶尔看到不太明白

第二个循环就没有瞬态了

lihaitao123 发表于 2013-7-12 16:45

hsfy919 发表于 2013-7-12 10:54 static/image/common/back.gif
第二个循环就没有瞬态了

求解释一下，他是积分时间都是从如何去掉顺瞬态的？
页: [1]
查看完整版本: 一本书上的分岔图程序