matlab 计算雅可比矩阵的lorenz李雅普图

1、内容简介

略
629-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

clear
clc
close all
x0_span = -5:0.01:5;
Z = [];
for i = 1:length(x0_span)ily(i,:) = ly_spectra(x0_span(i));
endfigure
plot(ly(:,[2 3 4 5]),'DisplayName','ly')
figure
plot(sum(ly,2))

function ly = ly_spectra(x1)
% clear
% clc
T = 13.5;
dt = 0.01;
a = 3.39;
b = 1.95;
c = 9.36;
r = 0.42;
g = 1.42;
A = 0.0676;
B = 0.3682;
N=T/dt;
x2=0;x3=0;x4=0; x5=0;
[t,x] = ode45('func',[0:1:T],[x1;x2;x3;x4;x5]);
n=length(t);
x1=x(n,1); x2=x(n,2); x3=x(n,3); x4=x(n,4); x5=x(n,5);
[t,x] = ode45('func',[0:dt:T],[x1;x2;x3;x4;x5]);JN = eye(5);
w = eye(5);
J = eye(5);
for k=1:Nx1 = x(k,1);x2 = x(k,2);x3 = x(k,3);x4 = x(k,4);x5 = x(k,5);W1 = A+B*abs(x4);W2 = A+B*abs(x5);J = (eye(5)+[-c*r-c/W2 c c 0 0; -a a*g a*g 0 0; -b b*g b*g-b*W1 0 0;0 0 1 0 0;1/W2 0 0 0 0]*dt);[~,w] = qr(J*w);w(:,1) = w(:,1)/norm(w(:,1));w(:,2) = w(:,2)/norm(w(:,2));w(:,3) = w(:,3)/norm(w(:,3));w(:,4) = w(:,4)/norm(w(:,4));w(:,5) = w(:,5)/norm(w(:,5));
end

4、参考论文

略

matlab 计算雅可比矩阵的lorenz李雅普图相关推荐

用matlab计算雅可比矩阵
如何用matlab轻松算出雅可比矩阵? 举例: 已知函数:f1=-x1+x2 f2=x1-x2-x2^3 求: 具体步骤: 1.matlab里定义2个变量 2.输入f的两个表达式 3.直接调用jaco ...
matlab画复变函数,科学网—复数复变函数的Matlab计算与绘图 - 周铁戈的博文
复数复变函数的Matlab计算与绘图周铁戈复数的表示存在两种表示方法,一种是代数式,一种是指数式,在Matlab中的方式如下: >> z=1+2i #代数式,1 ...
达摩院李雅亮：大规模预训练模型的压缩和蒸馏
作者 | 李雅亮博士阿里巴巴来源 | DataFunTalk 导读:本次分享的主题为大规模预训练模型的压缩和蒸馏,主要是从自动机器学习的角度,介绍大规模预训练模型的压缩和蒸馏.将介绍阿里巴巴达摩院 ...
MATLAB计算气象水文要素年内分配指数
MATLAB计算气象水文要素年内分配指数降水/径流等水文要素年内分配指标 1 不均匀系数nuniformity coefficient (Cn) 1.1 原理 1.2 MATLAB代码 1.3 案例 ...
matlab 计算变异系数,[转载]用Matlab编的计算结构可靠指标的改进一次二阶矩法程序（验算点法）...
题目:编制改进一次二阶矩法计算可靠指标的程序,并给出算例,要求提供源程序,算法语言不限. 选取的算例为:z=g(x,y)=x*y-1140,其中x,y服从正态分布,μx=38,Vx=0.1, μy=3 ...
如何用Matlab建立信用风险模型,基于Matlab计算的KMV模型在商业银行信用风险管理中的实践应用.pdf...
科技论坛 2014.6 基于Matlab计算的KMV模型在商业银行信用风险管理中的应用李园 (天津大学管理与经济学部,天津,300072) 摘要:社会经济的不断发展,金融行业也处于不断发展之中, ...
【数字信号处理】相关函数应用 ( 使用 matlab 计算相关函数 )
文章目录一.相关函数应用场景 1.生成高斯白噪声 2.信噪比 SNR 3.根据信噪比 SNR 求信号幅度 4.产生单载波信号及最终信号 5.求自相关函数及功率 6.matlab 完整代码一.相关函 ...
Matlab计算基站覆盖面积示例
站间距用D表示:小区半径为 R: 单扇区覆盖面积一般作为一个正六边形:D和R示意如下: 计算基站覆盖面积,一般算三个正六边形:下图填充部分: 公式如下: D一般算作R的1.5倍: 用matlab计算 ...
P2 Matlab计算基础-《Matlab/Simulink与控制系统仿真》程序指令总结
上一篇回到目录下一篇 <Matlab/Simulink与控制系统仿真>程序指令总结 Matlab_Simulink_BookExample 2. Matlab 计算基础表2.1 Ma ...

matlab 计算雅可比矩阵的lorenz李雅普图

matlab 计算雅可比矩阵的lorenz李雅普图相关推荐

最新文章

热门文章