基于奇异值分解的Lyapunov指数计算（2）
Chen_euler.m
% Chen图形(欧拉方法)
%
% Author：yujunjie
% Author's email: ustb03-07@yahoo.com.cn
%
clear
h=0.004;a=35;b=3;c=28;
x=10;y=10;z=20;
Y=[];
for i=1:5000
x1=x+h*a*(y-x);
y1=y+h*((c-a)*x-x*z+c*y);
z1=z+h*(x*y-b*z);
x=x1;y=y1;z=z1;
Y(i,:)=[x y z];
end
plot3(Y(:,1),Y(:,2),Y(:,3));
chen_le_eu.m
% 奇异值分解求Lyapunov法
% 微分chen系统
%
% Author：yujunjie
% Author's email: ustb03-07@yahoo.com.cn
%
h=0.002;a=35;b=3;c=28;
x=10;y=10;z=20;
V=eye(3);
S=V;b1=0;
k=6000;
for i=1:k
x1=x+h*a*(y-x);
y1=y+h*((c-a)*x-x*z+c*y);
z1=z+h*(x*y-b*z);
x=x1;y=y1;z=z1;
J=[-a a 0
c-a-z c -x
y x -b];
J=eye(3)+h*J;
B=J*V*S;
[V,S,U]=svd(B);
a_max=max(diag(S));
S=(1/a_max)*S;
b1=b1+log(a_max);
end
Lyapunov=(log(diag(S))+b1)/(k*h)
>> chen_le_eu
Lyapunov =
2.1531
0.0051
-11.8983
6 comments:
你好啊！看了你的很多文章，看你学的这么好，好羡慕你啊！我也在学混沌，现在在做硕士论文，是基于混沌的图像加密，我对这个不是很熟，不知道可不可以和你交流一下，我的QQ是215753800
你好，我看到你的QQ留言了，谢谢你跟我联系，不知道你的邮箱是？？
ustb03-07@yahoo.com.cn
老师,我看过你的Lorenz的matlab制作,并且演示了.很好!我能请教你关于一些matlab的制作!
0.05x-0.01y+x*y-0.003294>0,
0.000072012-0.01y* x+0.002118y-0.0054174x>0,
以上是一个不等式方程组,如何用matlab程序做出一个二维网格线性图!本来对这个就很陌生,查资料尝试过一些方法,但不行!
我的e-mail:gago999@126.com
QQ:30809190
谢谢!
我不是老师，另外没查到什么是二维网格线性图。所以抱歉帮不了忙。
Post a Comment