已知微分方程
0)1(
2
2
2
���� x
dt
dx
x
dt
xd
�
,讨论极限环.
在 matlab 环境下运行
clc
n=2000;
h=0.05;
x10=0.2;
x20=0.4;
x1=[];
x2=[];
%mu=0.05;
for mu=0.01:0.01:0.1
for i=1:n
x11=x10+x20*h;
x21=x20+(mu*(1-x10*x10)*x20-x10)*h;
x10=x11;
x20=x21;
x1=[x1 x11];
x2=[x2 x21];
end
end
plot(x1,x2)
hold on
m=1000;
x10=2;
x20=4;
x1=[];
x2=[];
%mu=0.05;
for mu=0.01:0.01:0.1
for i=1:m
x11=x10+x20*h;
x21=x20+(mu*(1-x10*x10)*x20-x10)*h;
x10=x11;
x20=x21;
x1=[x1 x11];
x2=[x2 x21];
end
end
plot(x1,x2,'r')
评论7