基于matlab的lorenz系统仿真研究精品.pdfVIP

基于 Matlab的 Lorenz 系统仿真研究 摘要：本文利用matlab这一数学工具对Lorenz系统进行了研究。首先使用matlab 分析求解Lorenz方程，利用matlab 的绘图功能，直观地观察了Lorenz 混沌吸引子的三维图形，并简单观察了Lorenz混沌系统对初值的敏感性； 然后对Lorenz系统进行仿真，比较分析在不同参数下的Lorenz系统仿真 结果；最后验证了通过添加反馈控制的方式，可以使Lorenz方程不稳定的 平衡点成为稳定的平衡点。 关键词：Lorenz系统；matlab；混沌系统 1.引言 Lorenz方程是由美国著名的气象学家Lorenz在1963年为研究气候变化，通过 对对流实验的研究，建立的三个确定性一阶非线性微分方程。这三个方程是混沌 领域的经典方程，Lorenz系统也是第一个表现奇怪吸引子的连续动力系统，具有 着举足轻重的作用。Lorenz方程的表达式如下： 其中，σ、μ、b为正实常数。 本文利用matlab这一数学工具，对Lorenz系统进行了研究，得到了仿真结果， 加深了对Lorenz系统的认识。 2.matlab求解Lorenz方程并绘图 首先建立m文件“Lorenz.m”来定义Lorenz方程，固定σ 10，μ30，b 8/3， 程序如下所示： function dx Lorenz(t,x) dx [-10*(x (1)-x (2));30*x (1)-x (2)-x (1)*x (3);x (1)*x (2)-2.6667*x (3)]; end 然后利用ode45命令来求解Lorenz方程并绘制图形，初值取x y z 0.1。程 序如下所示： clf x0 [0.1,0.1,0.1]; [t,x] ode45(Lorenz,[0,100],x0); subplot(2,2,1) plot(x(:,1),x(:,3)) title((a)) subplot(2,2,2) plot(x(:,2),x(:,3)) title((b)) subplot(2,2,3) plot(x(:,1),x(:,2)) title((c)) subplot(2,2,4) plot3(x(:,1),x(:,2),x(:,3)) title((d)) 运行上述程序，可得到如下波形： (a) (b) 60 60 40 40 20 20 0 0 -20 -10 0 10 20 30 -40 -20 0 20 40 (c) (d) 40 50 20 40 30 0 20