matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程

1. 洛伦兹方程求解

本文说明用Matlab工具箱求解洛伦兹方程的过程，并给出吸引子的三维动态图象.洛伦兹方程如下：

(1)这是一个自洽的方程组，求解过程如下：

(1) 建立自定义函数

function

dy=Lorenz(t,y) %

y(1)=x,y(2)=y,y(3)=z

dy=zeros(3,1);

dy(1)=10*(-y(1)+y(2));

dy(2)=28*y(1)-y(2)-y(1)*y(3);

dy(3)=y(1)*y(2)-8*y(3)/3;

(2)用ode45命令求解

[t,y]=ode45(@Lorenz,[0,30],[12,2,9]);

subplot(221);

plot(t,y(:,1));

subplot(222);

plot(t,y(:,2))

subplot(223);

plot(t,y(:,3))

subplot(224);

plot3(y(:,1),y(:,2),y(:,3))

view([20 42]);

(3)求解结果

(4)动态显示吸引子的绘制过程

[t,y]=ode45(@Lorenz,[0 30],[12 2 9]);

clf;

axis([-20

20 -25 25 10 50]);

view([20

42]);

hold

on;

comet3(y(:,1),y(:,2),y(:,3));%显示吸引子的绘制过程

(5)生成动画

[t,]=ode45(@Lorenz,[0

30],[12 2 9]);

m=moviein(100);

axis([-20 20 -25 25 -10

50]);

shading flat;

h=plot3(y(:,1),y(:,2),y(:,3));

for j=1:100

rotate(h,[0 0 1],1.8); %沿Z轴旋转

axis([-20 20 -25 25 -10 50]);

shading flat;

m(:,j)=getframe;

[X,map]=getframe;

imwrite(X,map,['Lz'int2str(j)'.bmp'],'bmp');%写入bmp文件

end

movie(m,5) %循环播放5次

(6)验证蝴蝶效应

clf;

hold on;

[t,u]=ode45('Lorenz',[0 6],[12 2 9]);

plot(t,u(:,3),'color','r');

[t,v]=ode45('Lorenz',[0 6],[12 2.01 9]);

plot(t,v(:,3),'color','b');

[t,w]=ode45('Lorenz',[0 6],[12 1.99 9]);

plot(t,w(:,3),'color','k');

matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程相关推荐

matlab求二元函数极值算法_最优化计算与matlab实现（3）——进退法
参考资料 <精通MATLAB最优化计算(第二版)> 数值实现 Matlab 2019a 目录石中居士:最优化计算与Matlab实现--目录zhuanlan.zhihu.com 进退法 ...
亮度均匀性 matlab,求：亮度保持的夜景图像直方图均衡算法 matlab程序
该楼层疑似违规已被系统折叠隐藏此楼查看此楼求:亮度保持的夜景图像直方图均衡算法 matlab程序我是大四学生,最近在做一个论文,头疼死了,不知道这个论文的程序怎么写.这里是matlab论坛我想 ...
【Matlab 六自由度机器人】运动学正解（附MATLAB机器人正解完整代码）
[Matlab 六自由度机器人]求运动学正解往期回顾前言正文一.运动学正解 1. 齐次变换矩阵 2. 总变换二.代码实现 1. 定义各连杆参数 2. 齐次变换矩阵及总变换 3. 代码运行结果 ...
matlab求单自由度振动方程为,单自由度振动方程与Matlab/Simulink求解
1.问题引用1:质量-弹簧-阻尼系统引用2:模型推导 2.运动方程 Step1: 将微分方程最高阶变量移到等式左边式1 Step2: 为每一阶微分式选择状态变量,最高阶除外 2.1 ' '' . ...
matlab求曲线每点的斜率_如何用matlab求出图中各条直线的斜率
展开全部 1.瞬时斜率62616964757a686964616fe4b893e5b19e31333431373236:求导. 方法:d(k)=(u(k)-u(k-1))/T(k);d(0)=0; d ...
matlab求连个书叫大叔,科学计算与MATLAB语言尔雅答案
[单选题]避雷器与被保护设备(). A. 串联于电源侧 B. 并联于电源侧 C. 串联或并联于电源侧 D. 之间无连接 [单选题]排气式避雷器是()的另一种叫法. A. 保护间隙 B. 管型避雷器 C ...
【MATLAB】求二元一次方程所有可能的解 + 画图
MATLAB求二元一次方程所有可能的解 + 画图画图输出结果相关知识点示例函数:x.y均为正整数,x∈[0, 10], y∈[0, 100], z = 3 * y - x, 求z = 0时,x ...
matlab 定义一个有自变量的方程_常微分方程：（第四章）高阶微分方程
参考<常微分方程>第三版(王高雄) 常微分方程王高雄第四章高阶微分方程_哔哩哔哩 (゜-゜)つロ干杯~-bilibiliwww.bilibili.com 对于高阶微分方程,线性部分 ...
matlab求解不等式方程组解集,matlab求不等式方程组的可行解
求高手用matlab解一个三元的方程组我算的你这个解都是无穷大.具体过程如下:第一步:建立M文件函数myf.mfunctionf=myf(x)f(1)=600+x(1)*cos(10)+x(2)*s ...

matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程

matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程相关推荐

最新文章

热门文章