非线性电路-Lorenz方程的Matlab求解

1、精选优质文档-倾情为你奉上题目：Lorenz方程的Matlab求解姓名：Webster-jie学号：xxxx班级：硕4019日期：2014年12月Lorenz方程的Matlab求解硕4019班 xxx xxxx 1963年美国麻省理工学院的气象学家E. Lorenz通过对对流实验的研究，得到了第一个表现奇怪吸引子的连续动力系统，该系统描述了从水桶底部加热时，桶内液体的运动情况，加热时，底部的液体越来越热，并开始逐渐上升，产生对流，当提供足够的热量并保持不变时，对流便会以不规则的和湍流的方式运动。通过对该动力学模型进行数值计算发现了一个由非线性微分方程组描述的著名的Lorenz方程，这

2、就是混沌现象的第一个奇怪吸引子Lorenz 吸引子。由于在天气、对流、斜波等现象及水轮机、发电机、激光机等真实物理系统中发现，Lorenz方程可以作为许多现实混沌运动的精确模型，因此对Lorenz方程的特性的研究受到许多学者的关注。一、Lorenz方程二、源程序clearh=0.005;%欧拉法，步长取0.005a=10;b=8/3;u=100;x=20;y=20;z=50;%起始点选为（20,20,50）Y=;for i=1:8000 x1=x+h*a*(y-x); y1=y+h*(u*x-x*z-y); z1=z+h*(x*y-b*z); x=x1; y=y1; z=z1; Y(

3、i,:)=x y z;endplot3(Y(:,1),Y(:,2),Y(:,3);figure(2)%x与t（i=8000）的关系plot(Y(:,1)figure(3)%y与t的关系plot(Y(:,2)figure(4)%z与t的关系plot(Y(:,3)三、程序运行结果及分析：图1（1）当0<u<1时，平衡点是稳定的，例如u=0.8时，奇点O（0,0,0）处三个特征值为（-10.8151，-0.1849，-2.6667），三个特征值均为负值，因此，全部轨道在时趋于零。如图2所示：图2（2）当时，例如当u=1.5时，系统有三个平衡点，在奇点O（0,0,0）处三个特征值分别为（-

4、11.4372,0.4372，-2.6667），一个特征值是正值，另外两个是负值，故平衡点O（0,0,0）是鞍点，因此平衡点O（0,0,0）不稳定；另外两个平衡点为P+，P-。图3分析图3，随着i从时（即随着时间变化），轨道从初始点（20,20,50）趋于平衡点零，之后离开平衡点。当时，平衡点O（0,0,0）是不稳定的，奇点O（0,0,0）处三个特征值有一个为正，两个为负，即为（-21.86,10.86，-2.6667），平衡点O（0,0,0）是鞍点，从而奇点O的稳定流形是一曲面，不稳定流形为一曲线。当时，不稳定流形最终螺旋地趋于与之同一侧的平衡点P+或P-，由图可以看出，曲线最终趋于与初始点

5、（20,20,50）同一侧的P+点。当u从1逐渐增大时，这种螺旋圈逐步增大。当时，不稳定流形刚好无限趋于原点O，即出现同宿轨。如图4所示：图4当时，不稳定流形将绕到另一侧，最终趋于与之异侧的平衡点P+或P-，由图5可看出，不稳定流形绕到了初始点（20,20,50）不同侧的P-，如下图所示。图5可见，是一个同宿分岔点。当时，奇点P+和P-变成不稳定的，也就是说u在24.74之后系统进入“混沌区”，对混沌区内的参量u值有时定常态是混沌解，即吸引子是奇怪吸引的，有时也可能出现稳定的周期解，如图6，这时轨线P+转几圈后被甩到P-附近转几圈，又回到P+附近，而且刚好头尾接上，成为闭曲线。图6当u很大时，取时，轨线绕P+、P-各半圈，形成一个空间闭轨道，在x-z或y-z平面上投影类似横“8”字形。如图7所示：图7四、结论 Lorenz方程是非线性微分方程，没有解析解，只有数值解。通过对本课程的学习，了解了非线性系统Lorenz方程在u取不同值情况下呈现出来的貌似无序但又遵循一定规律的复杂动力学行为。随着人类对混沌理论的不断探索和认识，在