急求：用matlab来画出Lorenz的关于参数r的分叉图程序

lirong75 · 发表于 2006-4-12 22:06

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

急求：用matlab来画出Lorenz的关于参数r的分叉图程序 老板要求做出来，很着急 谢谢各位大虾了

jumpwolf · 发表于 2006-4-13 12:12

clear all global r t0=[0 150];%积分时间 y0=[1,0,0];
%bifurcation for r=20:0.05:30 %r的变化精度 [t,y]=ode45('Lorenz',t0,y0); [Xmax]=getmax(y);
plot(r,Xmax,'b','markersize',1) hold on clear Xmax end xlabel('r') ylabel('Xmax') 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-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; 仅供参考！要想图漂亮一点，可以把积分时间和r的变化精度增大，但是计算时间会变长。

zhaolei · 发表于 2006-4-15 11:02

??? Undefined function or variable 'getmax'.怎么回出现这样的错误那

jumpwolf · 发表于 2006-4-15 12:30

getmax和 Lorenz是两个函数，要另外写两个M文件。

lirong75 · 发表于 2006-4-17 21:54

谢谢二楼的朋友

账号		自动登录	找回密码
密码			我要加入

[分形与混沌] 急求：用matlab来画出Lorenz的关于参数r的分叉图程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。