matlab 绘 洛伦兹系统 3D相图
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;
endplot3(x(:,1),x(:,2),x(:,3));
grid on
龙格库塔方法
其中
% ============================================================
% -- 洛伦兹方程的四阶龙格库塔函数
% ============================================================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
洛伦兹方程
%%******************************************
% -- 洛伦兹方程
%%******************************************A = 10;
B = 28;
C = 8/3;function g=f1(x,y,z)
g = A*(y - x);
endfunction g=f2(x,y,z)
g = B*x - y - x*z;
endfunction g=f3(x,y,z)
g = x*y - C*z;
end
matlab 绘 洛伦兹系统 3D相图相关推荐
- 使用matlab 仿真洛伦兹方程
参考 matlab 绘 洛伦兹系统 3D相图_颹蕭蕭的博客-CSDN博客_matlab三维相图 代码如下: % 洛伦兹系统仿真 clear all; % 洛伦兹系统 % % ************* ...
- matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程
1. 洛伦兹方程求解 本文说明用Matlab工具箱求解洛伦兹方程的过程,并给出吸引子的三维动态图象.洛伦兹方程如下: (1)这是一个自洽的方程组,求解过程如下: (1) 建立自定义函数 functio ...
- matlab中洛伦兹拟合,基于MATLAB洛伦兹线型非线性拟合算法实现
[1] Yao Hua. Research on Remote Sensing of Methane Based Tunable Diode Laser Absorption Spectroscopy ...
- R语言求解 常微分方程 绘制洛伦兹系统
library(deSolve) # ODE 常微分方程 library(scatterplot3d) # 三维曲线图library(shape) library(ReacTran) # PDE 偏微 ...
- 用matlab画多普勒加宽线性函数,洛伦兹线性函数
对于 CEST MRI 图像中的任一个像素,利 用伪佛克脱线型(PVP)代替洛伦兹作为拟合函数, 使拟合方法适合在更大的饱和功率和组织中磁化 转移含量更大的情况下拟合,...... 3.3谱线加宽和线 ...
- c语言模拟洛伦兹吸引子,洛伦兹吸引子相图的简易实现.pdf
洛伦兹吸引子相图的简易实现.pdf 洛伦兹吸引子相图的简易实现 晏志武 ( 临沧市临翔区第一中学 云南 临沧 ) ( 收稿日期 ) 摘要 用平面动画软件 绘制出高质量的洛伦兹吸引子相图 关键词 洛伦兹 ...
- 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 数据处理(二)非线性拟合——洛伦兹拟合(Lorentz fit)
文章目录 一.拟合示例 二.单峰洛伦兹 2.1 洛伦兹线型函数表达式与物理含义 2.2 lsqcurvefit非线性拟合 2.3 代码实现 三.双峰洛伦兹 3.1 洛伦兹线型函数表达式与物理含义 3. ...
- matlab做基尼曲线,计算基尼系数和matplotlib绘制洛伦兹曲线
基尼系数和洛伦兹曲线,在表示数据的不平均方面特别是财富的不平均上被广泛应用.但是目前在python里面并没有找到很好的可以直接绘制洛伦兹曲线的函数,由于目前项目用到,也就在实际应用中使用到,就把如何使 ...
最新文章
- 使用Hyperledger Ursa简化区块链安全性
- 如何布局文章标题才更吸引搜索引擎注意?
- 操作系统-进程同步与互斥易混点
- 汇编语言:利用栈的特性对内存数据进行逆置
- HTML+CSS制作Windows启动加载动画
- python初学者怎么入门-初学者如何学习Python?掌握这17个实用小技巧快速入门!...
- webstorm设置字体没有反应
- 月薪30k的PHP架构师的成长路线图1.0!
- CTF 每日一题 Day39 达芬奇密码
- iOS各种证书之间详解关于Certificate、Provisioning Profile、App ID的介绍及其之间的关系
- java.lang.IllegalStateException: Only fullscreen opaque activities can request orientation
- C++学习8——菜鸟教程自学
- 对于按字寻址和按字节寻址的理解
- ruby rails + grape + sidekiq 项目实践
- Html help安装步骤,Doxgen+Graphiz+htmlhelp配置
- 多因素身份认证 (MFA) 插件:手机验证码认证因素配置流程
- 程序员欣宸的文章分类汇总
- 用计算机里可以加50度电,一台电脑开一天用多少度电?台式机节电节能建议和措施...
- multisim 11
- 大学期间如何自学编程?