做分岔图时,程序中这个语句是什么含义(含程序)
本帖最后由 蓝云天 于 2011-9-8 21:31 编辑各位大哥大姐,做分岔图的时候,ode45后面括号里面有一项的含义搞不清楚,已经用红色标出来了,哪位指教一下里面都是什么含义?r应该是参数,其他的数值是什么意思呢?谢谢了?
clear all
t0=;%积分时间
%bifurcation
for r=linspace(1,500,1000);%r的变化精度
=ode45('Lorenz',t0,); =getmax(y(:,1));
plot(r,Xmax,'b','markersize',1)
hold on
clear Xmax
end
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-2,1)-b^2/(4*c);
j=j+1;
end
end 三维的,初值,其它是参数赋值。对应着lorenz系统应该能看出来!! 自己晕了,确实是,呵呵,刚才看了一下lorenz函数把系数也当状态向量处理了,我说怎么那么多初值。好久不用都忘了怎么回事了?谢谢lihaitao123 热心帮忙哈
页:
[1]