马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 蓝云天 于 2011-9-8 21:31 编辑
各位大哥大姐,做分岔图的时候,ode45后面括号里面有一项的含义搞不清楚,已经用红色标出来了,哪位指教一下里面都是什么含义?r应该是参数,其他的数值是什么意思呢?谢谢了?
clear all
t0=[0 100];%积分时间
%bifurcation
for r=linspace(1,500,1000); %r的变化精度
[t,y]=ode45('Lorenz',t0,[1;1;1;16;r;4]); [Xmax]=getmax(y(:,1));
plot(r,Xmax,'b','markersize',1)
hold on
clear Xmax
end
function [Xmax] = 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 |