clc;
clear;
close all;
warning off;
n=10000;
a=3:0.0001:4;
%这里取了0.001，比原题的更精细
len=length(a);
a=reshape(a,len,1);
sum=zeros(len,1);
unit=ones(len,1);
x=unit*0.1;
for i=1:n
y=a.*(unit-2*x);
sum=sum+log(abs(y));
x=a.*x.*(unit-x);
end
lamuda=sum/5000;
%根据评论修改过
plot(a,lamuda)
axis([3.6 4 -5 1])
grid on xlabel('u')
ylabel('Lyapunov指数\lambda')
title('指数Lyapunov 指数的关系曲线')
【MATLAB|【指数Lyapunov】Lyapunov指数的MATLAB仿真】
文章图片
推荐阅读
- MATLAB|模版匹配定位跟踪
- MATLAB|【智能驾驶】基于合成雷达的传感器融合的智能驾驶MATLAB仿真
- 其他|如何将本地文件上传到Gitee
- 嵌入式|2021-07-08 使用RT-Thread Studio程序调试过程中的要点记录（芯片使用STM32F407ZE）
- Matlab|Matlab实现小世界网络生成及其分析
- 学习笔记|C++ STL概述
- 图像算法|非局部均值滤波算法（NL-means）
- dfs|dfs bfs连通区域算法 matlab,【算法】图论(一) —— 基本图算法(BFS/DFS/强连通分量)...
- 算法优化|回归预测 | MATLAB实现PSO-RBF多输入单输出