心剑 发表于 2007-12-18 16:54

求duffing方程lyapunov的QR分解算法或RHR算法的程序。

求各位大哥大姐：
最近算duffing方程的lyapunov指数用定义法求解不够理想，有谁有QR分解，或RHR算法的程序。
咱先谢谢了。

无水1324 发表于 2007-12-19 10:18

回复 #1 心剑 的帖子

QR分解方法就是矩阵的一种计算方法吗？

心剑 发表于 2007-12-19 14:57

回复 #2 无水1324 的帖子

我看了李月写的一本书上说的求解lyapunov指数的方法，其中对标准QR分解算法，RHR算法以及RHR改进算法，看到得到数据与图还不错，但我是初学者看了一段时间还不知算法的具体意思。

octopussheng 发表于 2007-12-19 19:25

QR分解法，其实那个LET用的就是这种方法，需要看的话就看LET的源码吧！

心剑 发表于 2007-12-21 10:01

回复 #4 octopussheng 的帖子

LET的源码在哪里有啊？

octopussheng 发表于 2007-12-22 10:36

用google搜索一下吧！

心剑 发表于 2007-12-24 09:38

回复 #6 octopussheng 的帖子

多谢谢你了！

syjgf619 发表于 2010-4-10 13:30

急急急

求duffinf方程lyapunov的QR分解算法或RHR算法及RHR改进的程序

czy 发表于 2010-11-23 19:11

回复 5 # 心剑 的帖子

找到Let 源码了吗？我有

gghhjj 发表于 2010-11-26 07:07

LET工具包

gghhjj 发表于 2010-11-26 07:15

基于RHR算法的求Lyapunov程序

Lya1=[];Lya2=[];Lya3=[];
V=eye(3);
S=V;b1=0;
a=0.4;c=0.2;gama=3.5;
b=4.0;
h=0.01;
x(1)=0.1;y(1)=0;z(1)=0;n=0;
while z<=200
n=n+1;
k1=h*y(n);
m1=h*(-sin(x(n))-a*y(n)+b*cos(gama*z(n)).*sin(x(n))+c);
k2=h*(y(n)+m1/2);
m2=h*(-sin(x(n)+k1/2)-a*(y(n)+m1/2)+b*cos(gama*(z(n)+h/2)).*sin(x(n)+k1/2)+c);
k3=h*(y(n)+m2/2);
m3=h*(-sin(x(n)+k2/2)-a*(y(n)+m2/2)+b*cos(gama*(z(n)+h/2)).*sin(x(n)+k2/2)+c);
k4=h*(y(n)+m3);
m4=h*(-sin(x(n)+k3)-a*(y(n)+m3)+b*cos(gama*(z(n)+h)).*sin(x(n)+k3)+c);
x(n+1)=x(n)+(k1+2*k2+2*k3+k4)/6;
y(n+1)=y(n)+(m1+2*m2+2*m3+m4)/6;
z(n+1)=n*h;
J = [0 1 0;
b*cos(gama*z(n+1))*cos(x(n+1))-cos(x(n+1)) -a -b*gama*sin(gama*z(n+1))*sin(x(n+1));
0 0 0];
J=eye(3)+h*J;
B=J*V*S;
=svd(B);
a_max=max(diag(S));
S=(1/a_max)*S;
b1=b1+log(a_max);
Lyapunov1=(log(diag(S))+b1)/(n*h);
Lya1=;
Lya2=;
Lya3=;

end

Lyapunov1
n=1:20001;
plot(n,Lya1,'k',n,Lya2,'k',n,Lya3,'k')
%grid on
axis()
title('Lyapunov exponents of Warship')
xlabel('n'),ylabel('LCE')

zhao522 发表于 2011-3-18 19:16

回复 10 # gghhjj 的帖子

谢谢了 ，最近在学习Duffing混沌系统状态的判别，谢谢。

ranling 发表于 2011-3-21 09:28

回复 12 # zhao522 的帖子

可以联系我 qq 357809978 咱俩的方向差不多

shaoshuyi 发表于 2011-12-2 19:15

混沌都不知道怎么搞

shaoshuyi 发表于 2012-3-4 10:21

duffing系统lyapunov程序咋编？
页: [1]
查看完整版本: 求duffing方程lyapunov的QR分解算法或RHR算法的程序。