ymsncx 发表于 2008-6-2 15:09

为什么我的分岔图出不来

function dx=yang(t,x)
global A ;
dx=zeros(3,1);
dx(1)=6.3*x(2)-1.8*x(1)-A*(abs(x(1)-1)-abs(x(1)+1));
dx(2)=0.7*x(1)-0.7*x(2)+x(3);
dx(3)=-7.0*x(2);

function chua
clear;
global A
range=;
period=2*pi;
k=0;
YY1=[];
Step=2*pi/100;
for A=range
x0=;
A
k=k+1;   %出去前面60个周期的数据,并将最后的结果作为下一次积分的初值。
tspan=;
=ode45('yang',tspan ,x0);
y0=Y(end,:);
j=1;
for i=60:200
    tspan=;
    ode45('yang',tspan ,x0);
    YY1(k, j)=Y(end,1);
    j=j+1;    %取出每一个周期内的第一个解的最后一个值。
    y0=Y(end,:);
   end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k','markersize',1);
各位请都给点意见吧

咕噜噜 发表于 2008-6-2 19:07

ode45('yang',tspan ,x0);
这里有问题
不过你的程序有问题吧

function chua
clear;
global A
range=;
period=2*pi;
k=0;
YY1=[];
Step=2*pi/100;
for A=range
x0=;
A
k=k+1;   %出去前面60个周期的数据,并将最后的结果作为下一次积分的初值。
tspan=;
=ode45('yang',tspan ,x0);
y0=Y(end,:);
j=1;
for i=60:200
    tspan=;
    ode45('yang',tspan ,x0);
    YY1(k, j)=Y(end,1);
    j=j+1;    %取出每一个周期内的第一个解的最后一个值。
    y0=Y(end,:);
   end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k','markersize',1);
为什么出现两次 =ode45('yang',tspan ,x0);

ymsncx 发表于 2008-6-2 19:17

我是按照liliangbiao提供的Duffing系统的分岔图的Matlab实现中的程序编的,实现不了
页: [1]
查看完整版本: 为什么我的分岔图出不来