论坛 › 数学与非线性科学 › 请高手看看这分岔程序

earth 发表于 2007-10-30 12:49

请高手看看这分岔程序

clear all
global c
zhang=[];
M=;
counter=1;
for counter=1:length(M)
    c=M(counter);
    i=2;
xmax=0;
xmaxold=0;
frmdata=[];
error=0;
tspan=;
var=1;
y0=;
=ode45(@chaoLu2,tspan,y0);
while i < size(x,1)
       if x(i-1,var) < x(i,var) & x(i+1,var) <= x(i,var)
            xmax=x(i,var);
            if xmaxold ~= 0
                frmdata=;
            end
                xmaxold=xmax;
       end
      i=i+1;
end
r= length(frmdata)-20:length(frmdata);
fradata1(1,r)=frmdata(r,1);
zhang=;%这个地方老是显示错误，这是错误信息{All rows in the bracketed expression must have the same
number of columns.}%
end

plot(M,zhang,'k.','markersize',1);
xlabel(sprintf('c'));
ylabel(sprintf('x'))

function dx=chaoLu2(t,X)

x=X(1);

y=X(2);

z=X(3);
w=X(4);

global c;
dx=zeros(4,1);

dx(1)=-35*(x-y-yz);

dx(2)=-x*z+w+14*y;

dx(3)=-3*z+y*x;
dx(4)=-c*x
请高手能给我指导下，谢谢！

咕噜噜 发表于 2007-10-30 13:45

搜索论坛，用最大值法求分叉图
你这个就是最大值法

无水1324 发表于 2007-10-31 11:51

回复 #1 earth 的帖子

是程序出了问题还是什么？

查看完整版本: 请高手看看这分岔程序