lorenzp.m

来自「这是一个混沌理论中重要的Lorenz映射的程序。通过曲线很容易看到奇妙的混沌现象」· M 代码 · 共 15 行

function dy=Lorenzp(t,y); 
% 定义Lorenz方程 
% dx/dt=-b*x+z*y 
% dy/dt=-σ(y-z) 
% dz/dt=r*y-z-x*y 
% 取σ=10,b=8/3;r=28; 
% Example: 
%    [t,y]=ode45('Lorenzp',[0,80],[0;0;1e-10]); 
%    u=y(:,1);v=y(:,2);w=y(:,3);plot3(u,v,w); 
% 
global xigma b r 
xigma=10;b=8/3;r=28; 
dy=[-b*y(1)+y(2)*y(3);-xigma*(y(2)-y(3));... 
       r*y(2)-y(3)-y(1)*y(2)];