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

心剑

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

无水1324

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

心剑

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

octopussheng

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

心剑

LET的源码在哪里有啊？

octopussheng

用google搜索一下吧！

心剑

多谢谢你了！

syjgf619

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

czy

回复 5 # 心剑 的帖子

找到Let 源码了吗？我有

gghhjj

LET工具包

gghhjj

基于RHR算法的求Lyapunov程序

1. Lya1=[];Lya2=[];Lya3=[];
   
2. V=eye(3);
   
3. S=V;b1=0;
   
4. a=0.4;c=0.2;gama=3.5;
   
5. b=4.0;
   
6. h=0.01;
   
7. x(1)=0.1;y(1)=0;z(1)=0;n=0;
   
8. while z<=200
   
9. n=n+1;
   
10. k1=h*y(n);
    
11. m1=h*(-sin(x(n))-a*y(n)+b*cos(gama*z(n)).*sin(x(n))+c);
    
12. k2=h*(y(n)+m1/2);
    
13. 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);
    
14. k3=h*(y(n)+m2/2);
    
15. 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);
    
16. k4=h*(y(n)+m3);
    
17. m4=h*(-sin(x(n)+k3)-a*(y(n)+m3)+b*cos(gama*(z(n)+h)).*sin(x(n)+k3)+c);
    
18. x(n+1)=x(n)+(k1+2*k2+2*k3+k4)/6;
    
19. y(n+1)=y(n)+(m1+2*m2+2*m3+m4)/6;
    
20. z(n+1)=n*h;
    
21. J = [0 1 0;
    
22. 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));
    
23. 0 0 0];
    
24. J=eye(3)+h*J;
    
25. B=J*V*S;
    
26. [V,S,U]=svd(B);
    
27. a_max=max(diag(S));
    
28. S=(1/a_max)*S;
    
29. b1=b1+log(a_max);
    
30. Lyapunov1=(log(diag(S))+b1)/(n*h);
    
31. Lya1=[Lya1,Lyapunov1(1,:)];
    
32. Lya2=[Lya2,Lyapunov1(2,:)];
    
33. Lya3=[Lya3,Lyapunov1(3,:)];
    
34. 
    
35. end
    
36. 
    
37. Lyapunov1
    
38. n=1:20001;
    
39. plot(n,Lya1,'k',n,Lya2,'k',n,Lya3,'k')
    
40. %grid on
    
41. axis([0,30001,-0.8,0.5])
    
42. title('Lyapunov exponents of Warship')
    
43. xlabel('n'),ylabel('LCE')

复制代码

zhao522

回复 10 # gghhjj 的帖子

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

ranling

回复 12 # zhao522 的帖子

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

shaoshuyi

混沌都不知道怎么搞

shaoshuyi

duffing系统lyapunov程序咋编？