基于matlab的Lorenz系统模拟实验仿真.docVIP

基于matlab的Lorenz系统模拟实验仿真 2010年第17期SCIENCEamp;TECHNOLOGYINFORMATIONO本刊重稿.科技信息 基于matlab的Lorenz系统模拟实验仿真 赖宏慧陈澜祯 (赣南医学院江西赣州341000) 【摘要】非线性系统的研究难度一直较大,借助数学工具madab进行模拟实验是目前研究的趋势.本文选取Lorenz系统为实验模型,探 讨了采用rnatlab对Lorenz系统实验的方法.快速求解能够获得实验结果以及模拟实验仿真图,最后使用mat.1ab的动态建模环境simulik模拟 仿真出的系统相图与前面实验比较.实验证明了模拟实验的可行性和通用性. 【关键词】madab;混沌系统;模拟仿真;lorenz系统 非线性是～切动力学复杂性之源.自然界和现实生活中所有系统 都是非线性的.混沌是非线性动力学系统所特有的一种运动形式,是 自然界及人类社会中的一种普遍现象,它是在一个确定性系统中出现 的一种貌似不规则的,内在的随机性运动.其动力学模型多采用微分 方程.微分方程求解难度较大,随着计算机技术的迅速发展,普遍借助 计算机数值方法予以解决.Matlab是一套数值分析软件,从上个世纪 8O年代发展至今,已经成为国际控制界公认的标准计算软件,用于算 法开发,数据可视化,数据分析以及数值计算的高级技术计算语言和 交互式环境.在数值计算方面独占鳌头_l-31. 1963年,美国麻省理工大学的气象学家E.N.Lorenz在研究大气运 动时提出了着名的Lorenz混沌模型.开辟了混沌发展的新纪元. Lorenz系统无论从数学还是物理的角度来说都值得详细地研究.本文 选取Lorenz系统为实验模型.利用maflab对Lorenz方程计算分析.快 速求解得到实验结果以及模拟实验仿真图.实验证明Matlab进行 Lorenz方程的模拟实验快速有效,具有一定的通用性. 1matlab分析求解Lorenz方程 1963年,气象学家Lorenz发现了第一个结构简单的三维自治混 沌系统一renz系统,其动力学方程为: f(I)=a(I)一x(1)) {(【):cx(I)一x(【)z(1)一y(1)(1) 【i(t)=x(1)y(t)一bz(t) 当系统参数a=10,b=8/3,c=28时,Lorenz系统处于混沌状态. 接下来描述采用matlab进行求解的详细步骤: 1.1建立自定义函数用edit命令或file菜单新建一个自定义函数 名为lorenz_f的m文件.代码为 functiondx=lorenz_ f(t,x) a=1O;b=8/3;c=28; dx(1)-a(1)一x(2)); d,【(2)=cx(1)-x(2)-x(1)x(3); dx(3)=x(1)x(2)-bx(3); d)【=【d]((1);dx(2);dx(3)】; 】.2采用4—5阶Runge—Kutta方程,利用ode45命令求解.同上再新 建一个lorenz.m文件,输入程序代码为: clear%清除内存变量和函数 set(gcLColor,【111]); set(gef,qnvertHardcopy,,0ff3; tspan=]O:O.01:10o]; x0=[0.2;3;1.31; [t,x]=ode45(lorenz_f,tspan,xo); subplot(221) plot(x(:,1),x(:,2));%Y平面上的相图 xlabel(~itx);ylabel(kity; subplot(222) plot(x(:,1),x(:,3;%xg平面上的相图 xlabel(kitx);ylabel(~itz subplot(223) plot(x(:,2),x(:,3));%y-z平面上的相图 xlabel(~ty);ylabel(~itz1; subplot(224) plot(t,x(:,1),r,x(:,2),g,t,x(:,3),b);%红绿 18 %蓝分别对应x,y.z三个分量的时间序列图 1.3在MATLAB窗口中运行程序Lorenz.m.可以得到Lorenz混沌系 统在x—Y,x—z,y-z平面上的相图,以及系统三个变量x—y-z的时间序 列图 1.4动态显示吸引子的绘制过程,新建m文件lorenz1.m,输入代码为 tspan=[0:0.01:100]; xO=[O.2;3;1.31; [t,x]=ode45(loren』tspan,xo); axis@,50,-20,20,一30,30]) view(I-37.5,30】); holdon comet3(x(:,3),x(:,1),x(:,2));%显示吸引子绘 %制过程 zlabel(kitz);xlabel(adtx);ylabel(%ty1; 可以得到lor