Lyapunov计算程序.docVIP

Lyapunov计算程序 2008-05-22 22:48 % 奇异值分解求Lyapunov法 % 微分lorenz系统 % dx/dt = SIGMA*(y - x) % dy/dt = BETA*x - y -x*z % dz/dt= x*y - GAMA*z % In this demo, SIGMA = 16,BETA = 45.92, GAMA = 4 % 理论值: % LE1 = 1.5, LE2 = 0.00, LE3 = -22.5 sigma=16;beta=45.92;gama=4; x(1)=20;y(1)=10;z(1)=50;V=eye(3);S=V;I=V;b=0; h=0.005;k=6000; for i=1:k k11=sigma*(-x(i)+y(i)); k21=beta*x(i)-x(i)*z(i)-y(i); k31=x(i)*y(i)-gama*z(i); k12=sigma*(-(x(i)+0.5*h*k11)+(y(i)+0.5*h*k21)); k22=beta*(x(i)+0.5*h*k11)-(x(i)+0.5*h*k11)... *(z(