laogao 发表于 2006-10-25 18:57

求助:分叉程序问题

clear all
global r
t0=;%积分时间
y0=;

%bifurcation
for r=20:0.05:30 %r的变化精度
=ode45('Lorenz',t0,y0);
=getmax(y);

plot(r,Xmax,'b','markersize',1)
hold on
clear Xmax
end
xlabel('r')
ylabel('Xmax')


function = getmax(y)
a=length(y);
j=1;
for i=(a-1)/2:a

b=(y(i,1)-y(i-2,1))/2;
c=(y(i,1)+y(i-2,1))/2-y(i-1,1);

if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
Xmax(j)=y(i-1,1);
j=j+1;
elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
Xmax(j)=y(i-1,1)-b^2/(4*c);
j=j+1;
end
end

function dy = Lorenz(t,y)
global r
dy=zeros(3,1);
dy(1)=-10*(y(1)-y(2));
dy(2)=-y(1)*y(3)+r*y(1)-y(2);
dy(3)=y(1)*y(2)-8*y(3)/3;

这个程序中取点的语句看不懂
b=(y(i,1)-y(i-2,1))/2;
c=(y(i,1)+y(i-2,1))/2-y(i-1,1);

if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
Xmax(j)=y(i-1,1);
j=j+1;
elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
Xmax(j)=y(i-1,1)-b^2/(4*c);
就是这几个语句.请大家指点.谢谢

liujie_hch 发表于 2006-11-19 17:06

请多情清秋解释一下,好象是他发的.

呵呵: 请多情清秋解释一下,好象是他发的.
页: [1]
查看完整版本: 求助:分叉程序问题