Programación de noviembre
Como el método Nicolis, el método Benettin, el método Wolf, el método Jacobia, etc.
Te dejo mi programa anterior para calcular el índice de Lie del mapeo de Enon para tu referencia:
- Calcular el índice de Laypunov máximo del atractor de Enon.
¿Está claro? all;clc
a = 0.9:0.001:1.4;k=length(a);
b = 0.3p = 600
¿Para qué? n=1：k
? ¿para qué? m=2: p
x(1,n)= 0.4; y(1,n)= 0.6;
x(m,n)=1 b*y(m -1, n)-a(n)*x(m-1, n)^2;
y(m, n)=x(m-1, n);
? Fin
Fin
¿Para qué? R=1: k calcula la matriz jacobiana
? ¿para qué? h=2: p
A{1, r} =;
A{h, r}=*A{h-1, r}; Note la multiplicación de matrices de celdas; orden.
? Fin
Fin
¿Para qué? T=1: k calcula el índice de mentira máximo.
vv(：,t)=eig(A{p,t}); v = max(ABS(vv));
le 1 = 1/p * log( v);
Fin
plot(a, LE1, ' k '); Abierto;
plot(a, 0, 'k:');
Eje([a(1),a(k),-1?1]);
p>
xlabel(' a '); ylabel(' le 1 '); título ("Índice de vida máxima");