sode3.m

来自「非线性动力系统中分叉图用来定性的表示系统行为的演化」· M 代码 · 共 20 行

global A ;
for A = 2.082:-0.005:2.06;
[ T ,X] = ode45 ('ode1',[0,1500 ] , [ 0.05 ;0.05 ;0.05 ]) ;
data = X(:,1);
n = length (data) ;
N = round (n/2) ;
a = data (N - 2) ;
b = data (N - 1) ;
for i = N:n
if b>=a & b>=data(i) ;
plot (A ,b) ;hold on ;
axis on ;
axis ( [2.01 2.09 4 9 ]) ;
else
end
a = b ;
b = data (i) ;
end
end