求助duffing方程的分岔图程序
求救!!duffing方程的分岔图程序???[ 本帖最后由 xinyuxf 于 2007-9-4 10:58 编辑 ]
回复 #1 董世勇 的帖子
在论坛中搜索有关分岔的程序。将m程序修改为自己的方程就可以了 % Author: Thomas Lee% E-mail: lixf1979@126.com % Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
function dx=duffing(t,X)
global F wd;
r=0.168;
x=X(1);
y=X(2);
psi=X(3);
dx=zeros(3,1);
dx(1)=y;
dx(2)=-r*y+1/2*x*(1-x^2)+F*sin(psi);
dx(3)=wd;
clear;
global F wd;
wd=1.0;
range=;
period=2*pi/wd; %
k=0;
YY2=[];
step=2*pi/100;%步长。
for F=range
y0=;
F
k=k+1;
% discard the first 60 periodic data;
%除去前面60个周期的数据,并将最后的结果作为下一次积分的初值
tspan=;
=ode45(@duffing,tspan,y0);
y0=Y(end,:);
j=1;
for i=60:200
tspan=;
=ode45(@duffing,tspan,y0);
YY1(k,j)=Y(end,1); % get the omega data from every period end
j=j+1; %取出每一个周期内的第一个解的最后一个值。
y0=Y(end,:);
end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k.','markersize',1); {:{28}:}厉害的,先拜读一下。 好好学习一下,谢谢 回复 2 # 无水1324 的帖子
无水老师,怎么上传超过500k的图片 回复 3 # liliangbiao 的帖子
我一直有个疑问,这个因为wd=1,所以无所谓,如果不是的话,step是不是应该为step=period/100啊!还请指教 xiexie,分享 不错 感谢分享,可是为什么不出来图呢?
页:
[1]