使用matlab 仿真洛伦兹方程
参考
matlab 绘 洛伦兹系统 3D相图_颹蕭蕭的博客-CSDN博客_matlab三维相图
代码如下:
% 洛伦兹系统仿真
clear all;
% 洛伦兹系统
%
% ****************************************************
% * 参数初始化
% ****************************************************
%
x0 = 0.1;
y0 = 0.1;
z0 = 0.1;
n = 1e4;
x = zeros(n, 3);
x(1,:) = [x0, y0, z0];
% ****************************************************
% * 调用龙格库塔迭代公式进行迭代
% ****************************************************
%
for k=2:n
[dx,dy,dz] = dxdt_Lorenz(x(k-1,1),x(k-1,2),x(k-1,3));
x(k,1) = x(k-1,1) + dx;
x(k,2) = x(k-1,2) + dy;
x(k,3) = x(k-1,3) + dz;
end
plot3(x(:,1),x(:,2),x(:,3));
grid on
%%******************************************
% -- 洛伦兹方程
%%******************************************
function g=f1(x,y,z)
A = 10;
g = A*(y - x);
end
function g=f2(x,y,z)
B = 28;
g = B*x - y - x*z;
end
function g=f3(x,y,z)
C = 8/3;
g = x*y - C*z;
end
% ============================================================
% -- 洛伦兹方程的四阶龙格库塔函数
% ============================================================
function [dx,dy,dz] = dxdt_Lorenz(x,y,z)
h = 1e-2; %步长
K1 = f1(x,y,z);
K2 = f1(x + h*K1/2,y + h*K1/2,z + h*K1/2);
K3 = f1(x + h*K2/2,y + h*K2/2,z + h*K2/2);
K4 = f1(x + h*K3,y + h*K3, z + h*K3);
L1 = f2(x,y,z);
L2 = f2(x + h*L1/2,y + h*L1/2,z + h*L1/2);
L3 = f2(x + h*L2/2,y + h*L2/2,z + h*L2/2);
L4 = f2(x + h*L3,y + h*L3, z + h*L3);
M1 = f3(x,y,z);
M2 = f3(x + h*M1/2,y + h*M1/2,z + h*M1/2);
M3 = f3(x + h*M2/2,y + h*M2/2,z + h*M2/2);
M4 = f3(x + h*M3,y + h*M3, z + h*M3);
dx = (K1 + 2*K2 + 2*K3 + K4)*h/6;
dy = (L1 + 2*L2 + 2*L3 + L4)*h/6;
dz = (M1 + 2*M2 + 2*M3 + M4)*h/6;
end
我的理解:
-----
----
结果:x y z在空间上的轨迹。每个点代表一个时刻。
使用matlab 仿真洛伦兹方程相关推荐
- matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程
1. 洛伦兹方程求解 本文说明用Matlab工具箱求解洛伦兹方程的过程,并给出吸引子的三维动态图象.洛伦兹方程如下: (1)这是一个自洽的方程组,求解过程如下: (1) 建立自定义函数 functio ...
- matlab 洛伦兹方程,求解洛伦兹方程的MATLAB程序
求解洛伦兹方程的MATLAB程序 clear;a=0;b=100;h1=0.0001;h2=0.1;ya=5;2;10;sigma=10;gamma=28;rho=8/3;z=f2(sigma,gam ...
- matlab 绘 洛伦兹系统 3D相图
matlab 代码 clear all;% % **************************************************** % * 参数初始化 % *********** ...
- 洛伦兹方程c语言,龙格库塔在C为洛伦兹方程
我想计算洛伦茨系统使用龙格库塔方法,但我找不到我的代码有错误的地方.当我运行它时,系统进入静态点,我应该获得一个蝴蝶(洛伦兹吸引子).我认为这是Runge Kutta方法部分的'for'循环中的一些东 ...
- C语言 最小二乘 向量旋转 欧拉方法求洛伦兹方程
最小二乘 This question will develop a set of functions to least square fit the linear model
- matlab中洛伦兹拟合,基于MATLAB洛伦兹线型非线性拟合算法实现
[1] Yao Hua. Research on Remote Sensing of Methane Based Tunable Diode Laser Absorption Spectroscopy ...
- 用matlab画多普勒加宽线性函数,洛伦兹线性函数
对于 CEST MRI 图像中的任一个像素,利 用伪佛克脱线型(PVP)代替洛伦兹作为拟合函数, 使拟合方法适合在更大的饱和功率和组织中磁化 转移含量更大的情况下拟合,...... 3.3谱线加宽和线 ...
- 【Scipy高级计算】(2) 常微分方程、洛伦兹吸引子,附python完整代码
大家好,在上一篇博文中,我介绍了如何使用Scipy库计算定积分和二重积分,感兴趣的可以看一下:https://blog.csdn.net/dgvv4/article/details/124226759 ...
- MATLAB 数据处理(二)非线性拟合——洛伦兹拟合(Lorentz fit)
文章目录 一.拟合示例 二.单峰洛伦兹 2.1 洛伦兹线型函数表达式与物理含义 2.2 lsqcurvefit非线性拟合 2.3 代码实现 三.双峰洛伦兹 3.1 洛伦兹线型函数表达式与物理含义 3. ...
最新文章
- 构造函数不能为虚/重载函数总结
- 科学计算机怎么算坐标,用科学计算器fx-5800算坐标怎么按
- React state和props使用场景
- 1000道Python题库系列分享19(81道填空判断题)
- Python的web相关及Django简介
- 关于windows2008重新启动需要重新设置分辨率的问题
- 飞翔的小鸟(FlyBird)游戏C语言编程(含撞柱子)
- 巨佬Jake Wharton谈Android对Java 8的支持
- 内存取证-volatility工具的使用 (史上更全教程,更全命令)
- 使用redis存储全球IP库
- mysql 数据转移历史表_mysql 历史数据表迁移方案
- 用STC12C5A60S2的P4接口点亮LED
- jfinal获取mysql_jfinal 调用mysql结果集-问答-阿里云开发者社区-阿里云
- ADAMS-Simulink联合仿真-零基础(一)
- Caused by: java.lang.IllegalArgumentException的解决方法
- 2021 Python入门资料汇总
- 来,看一个真实的用户分析案例!
- CSS———移动端案例分享(美图公司首页)
- 风口猪炒股指标操作复盘 2021-01-20
- 百度广告产品系统级测试技术演进