董世勇 发表于 2007-9-4 08:27

求助duffing方程的分岔图程序

求救!!duffing方程的分岔图程序???

[ 本帖最后由 xinyuxf 于 2007-9-4 10:58 编辑 ]

无水1324 发表于 2007-9-4 09:31

回复 #1 董世勇 的帖子

在论坛中搜索有关分岔的程序。将m程序修改为自己的方程就可以了

liliangbiao 发表于 2007-9-5 18:07

% 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);

moneylancy 发表于 2010-11-3 18:35

{:{28}:}厉害的，先拜读一下。

chnwangnannan 发表于 2011-5-5 18:42

好好学习一下，谢谢

lihaitao123 发表于 2011-5-6 16:43

回复 2 # 无水1324 的帖子

无水老师，怎么上传超过500k的图片

kangarooli 发表于 2011-5-6 16:57

回复 3 # liliangbiao 的帖子

我一直有个疑问，这个因为wd=1，所以无所谓，如果不是的话，step是不是应该为step=period/100啊！还请指教

zhangyb303 发表于 2012-3-4 10:45

xiexie,分享 不错

亦殊心 发表于 2012-3-7 15:31

感谢分享，可是为什么不出来图呢？
页: [1]
查看完整版本: 求助duffing方程的分岔图程序