[Interpretação e significado do espectro de Lyapunov:] Considere um sistema dinâmico tridimensional não linear genérico descrito por
\dot{\vec{x}}=\vec{F}(\vec{x}),
com \vec{x}=(x, y, z). O espectro de Lyapunov obtido numericamente para uma trajetória de longa duração é:
\lambda_{1}=0,125 ; \quad \lambda_{2}=0,000 ; \quad \lambda_{3}=-14,60
A) Explique o significado físico de cada expoente e o que esse conjunto de valores indica sobre o comportamento da trajetória no espaço de fase. B) Determine qualitativamente o tipo de atrator associado a esse espectro (ponto fixo, ciclo limite, toro ou atrator caótico). C) Calcule a dimensão de Kaplan-Yorke (dimensão de Lyapunov) do atrator e interprete o resultado. D) Descreva o que aconteceria com o espectro se o sistema fosse estabilizado por controle linear que eliminasse o comportamento caótico. [\mathbf{2 , 5}]

Passo 1

Enunciado e dados fornecidos (o que o problema pede e quais informações temos):
  • Temos um sistema autônomo tridimensional com estado \vec{x}=(x,y,z) e equações \dot{\vec{x}}=\vec{F}(\vec{x}).
  • O espectro de Lyapunov numérico para uma trajetória de longa duração é: \lambda_{1}=0{,}125;\qquad \lambda_{2}=0{,}000;\qquad \lambda_{3}=-14{,}60.
  • Perguntas que o enunciado faz:
    A) Explicar o significado físico de cada expoente e o que esse conjunto indica sobre o comportamento da trajetória no espaço de fase.
    B) Determinar qualitativamente o tipo de atrator associado a esse espectro (ponto fixo, ciclo limite, toro ou atrator caótico).
    C) Calcular a dimensão de Kaplan–Yorke (dimensão de Lyapunov) do atrator e interpretar o resultado.
    D) Descrever o que aconteceria com o espectro se o sistema fosse estabilizado por controle linear que eliminasse o comportamento caótico.

Passo 2

O que é um expoente de Lyapunov (aula rápida e intuitiva)
  • Um expoente de Lyapunov \lambda mede a taxa exponencial média de divergência (ou convergência) entre duas trajetórias infinitesimalmente próximas no espaço de fase. Mais precisamente, se uma pequena perturbação inicial \delta \vec{x}(0) evolui como \delta \vec{x}(t)\approx \delta \vec{x}(0)\,e^{\lambda t} (num sentido médio), então \lambda é a taxa no expoente.
  • Sinais e interpretações rápidas:
    • \lambda>0: separação exponencial -> sensibilidade a condições iniciais (componente instável).
    • \lambda=0: nenhum crescimento nem decaimento exponencial médio -> direção neutra (por exemplo, deslocamento temporal em sistemas autônomos).
    • \lambda<0: contração exponencial -> direção estável.
  • Num sistema de dimensão n obtemos n expoentes ordenados \lambda_1\ge\lambda_2\ge\cdots\ge\lambda_n. Eles descrevem como volumes infinitesimais evoluem: a soma \sum_i \lambda_i é a taxa média de variação logarítmica do volume de um pequeno paralelepípedo no espaço de fase.

Passo 3

Por que temos um expoente zero em sistemas autônomos (importante para interpretar \lambda_2=0{,}000)
  • Em sistemas autônomos (sem dependência explícita do tempo), se deslocarmos uma solução um pequeno tempo \delta t ao longo da mesma trajetória obtemos outra solução. Esse deslocamento temporal gera uma perturbação que não cresce nem decai exponencialmente — é uma direção neutra associada ao fluxo; portanto, costuma existir sempre pelo menos um expoente Lyapunov igual a zero.
  • Conclusão prática: o valor \lambda_{2}=0{,}000 provavelmente corresponde exatamente a essa direção tangente ao fluxo (é comum encontras um zero num sistema contínuo).

Passo 4

Interpretação individual dos três expoentes (ensinando como ler os números)
  • \lambda_{1}=0{,}125>0:
    • Significado: existe uma direção no espaço de fase onde perturbações crescem, em média, como e^{0{,}125 t}. Isso traduz sensibilidade a condições iniciais nessa direção — característica associada a comportamento caótico/estranho.
    • Intuição para prova: um expoente positivo é a assinatura matemática de divergência exponencial local entre trajetórias próximas; é o que dá imprevisibilidade a médio/longos prazos.
  • \lambda_{2}=0{,}000:
    • Significado: direção neutra. Em um sistema autônomo isso normalmente é a direção tangente à órbita (deslocamento no tempo), ou seja, mudar a fase temporal não altera a forma da solução.
    • Intuição: não confunda “zero numérico” com insignificância — aqui é esperado e tem interpretação clara.
  • \lambda_{3}=-14{,}60<0:
    • Significado: existe uma direção muito fortemente contraída, com taxas médias de decaimento muito rápidas (ordem e^{-14{,}60 t}). Isso “aperta” as trajetórias na direção transversal correspondente, comprimindo volumes.
    • Intuição: um expoente negativo grande em magnitude indica que o atrator é fino numa direção — forte contração transversal.

Passo 5

O que diz o conjunto completo \{\lambda_1,\lambda_2,\lambda_3\} sobre o comportamento global (volume, dissipatividade, caos)
  • Soma dos expoentes: \sum_{i=1}^3 \lambda_i = 0{,}125 + 0{,}000 - 14{,}60 = -14{,}475.
    • Interpretação: a soma é negativa, logo o sistema é dissipativo: volumes no espaço de fase encolhem, em média, com taxa exponencial \exp(-14{,}475\,t) (muito rápido).
  • Presença de expoente positivo:
    • Como \lambda_1>0, existe sensibilidade a condições iniciais -> há possibilidade de comportamento caótico (se houver um conjunto compacto invariant e recorrente associado). Em sistemas dissipativos com pelo menos um expoente positivo costuma-se obter um atrator estranho (caótico) de dimensão fracionária.
  • O zero:
    • Garante que a dinâmica não é simplesmente uma contração para um ponto em tempo estacionário (se fosse ponto fixo, todos os expoentes seriam negativos). O zero é compatível tanto com órbitas periódicas (um zero) quanto com fluxo caótico contínuo (também um zero por ser autônomo).
  • Conclusão qualitativa (prévia, para usar depois na escolha do tipo de atrator): combinação de um expoente positivo, um zero e um grande negativo é típica de um atrator caótico (estranho) em sistemas 3D dissipativos.

Passo 6

Decisão do tipo de atrator (razão prática entre as opções pedidas)
  • Vamos comparar com características dos candidatos:
    • Ponto fixo estável: todos os \lambda_i<0. Aqui temos \lambda_1>0, então não é ponto fixo estável.
    • Ciclo limite (órbita periódica): típico em sistemas contínuos tem um expoente zero (direção de fase) e os demais negativos (um negativo no caso 3D). Para um ciclo limite não haveria \lambda>0. Aqui temos \lambda_1>0, logo não é um ciclo limite estável.
    • Toro (toroidal/quasi-periódicidade de dimensão 2): um toro de duas frequências em um fluxo contínuo tende a apresentar dois expoentes nulos (um por fluxo, outro por fase invariante) e os restantes negativos; portanto esperaríamos algo como \lambda_1=0,\;\lambda_2=0,\;\lambda_3<0. Não é o caso.
    • Atrator caótico (estranho): caracteriza-se por pelo menos um \lambda>0, geralmente um zero (fluxo), e soma negativa se é dissipativo. Isto coincide com os nossos valores.
  • Portanto (justificando didaticamente): o conjunto dado é compatível com um atrator caótico (estranho) dissipativo em 3D.

Passo 7

Cálculo da dimensão de Kaplan–Yorke (dimensão de Lyapunov) — fórmula, significado e cálculo passo a passo
  • Fórmula de Kaplan–Yorke (explicação):
    • Ordene os expoentes decrescentemente \lambda_1\ge\lambda_2\ge\cdots.
    • Seja j o maior inteiro tal que \sum_{i=1}^{j}\lambda_i\ge 0.
    • A dimensão de Kaplan–Yorke é D_{KY} = j + \frac{\sum_{i=1}^{j}\lambda_i}{\big|\lambda_{j+1}\big|}.
    • Intuição: j conta o número inteiro de direções cuja soma de expansões ainda é não negativa; a fração acrescenta a porção “fracionária” necessária na direção seguinte para anular o crescimento — resultando numa dimensão fractal entre j e j+1.
  • Aplicando aos nossos expoentes:
    • Ordenados já estão: \lambda_1=0{,}125, \lambda_2=0{,}000, \lambda_3=-14{,}60.
    • Soma até j=1: \lambda_1=0{,}125>0 então j\ge 1.
    • Soma até j=2: \lambda_1+\lambda_2=0{,}125+0{,}000=0{,}125\ge 0 então j\ge 2.
    • Soma até j=3: \lambda_1+\lambda_2+\lambda_3=-14{,}475<0 então j=2.
    • Agora aplicar a fórmula com j=2: D_{KY} = 2 + \frac{\lambda_1+\lambda_2}{\big|\lambda_3\big|} = 2 + \frac{0{,}125}{14{,}60}.
    • Calcule a fração: \frac{0{,}125}{14{,}60}\approx 0{,}0085616.
    • Assim D_{KY}\approx 2{,}00856\quad(\text{aprox. }2{,}009).
  • Interpretação do valor:
    • A dimensão D_{KY} está apenas um pouco acima de 2: isso significa que o conjunto atraente ocupa essencialmente uma superfície (dimensão 2) no espaço de fase, mas com uma estrutura fractal muito fina que lhe dá uma dimensão fractal não inteira — ou seja, um “folheado” fractal de dimensão ligeiramente maior que 2.
    • Em termos físicos: existe expansão em uma direção (produzindo complexidade) e forte compressão em outra, resultando num conjunto que é estendido em duas direções mas comprimido na terceira, com pequenas folgas que geram fractalidade.

Passo 8

Efeito de um controle linear que elimina o comportamento caótico — como o espectro mudaria (intuição e relação com autovalores do linearizado)
  • O que significa “estabilizar por controle linear que elimina o comportamento caótico”:
    • Um controle linear aplicado tipicamente altera o campo \vec{F}(\vec{x}) para que o ponto de operação desejado (ponto fixo ou órbita periódica) se torne assintoticamente estável.
    • Localmente, os expoentes de Lyapunov são muito relacionados às partes reais dos autovalores da linearização (Jacobiano) no regime de operação.
  • Possíveis destinos e seus espectros qualitativos:
    • Se o controle leva o sistema a um ponto fixo estável: todos os expoentes tornam-se negativos. Assim \lambda_1<0,\;\lambda_2<0,\;\lambda_3<0. Não há expoente zero (a menos que haja uma degenerescência não típica); o sistema deixa de ser sensível às condições iniciais.
    • Se o controle estabiliza uma órbita periódica (ciclo limite) estável: o espectro terá um expoente zero (direção de fase) e os demais negativos; ou seja, \lambda_1=0,\;\lambda_2<0,\;\lambda_3<0 (ou ordenados adequadamente). Não haverá expoente positivo.
    • Em termos gerais: a principal mudança é que o expoente positivo \lambda_1 deixará de ser positivo — será empurrado para valores não positivos pelo controle. O controle linear age como “deslocar para a esquerda” as partes reais dos autovalores, e assim transforma \lambda_1 em \le 0.
  • Consequências sobre soma e dissipatividade:
    • A soma provavelmente continuará negativa (ou ficará ainda mais negativa se a contração for reforçada). A principal mudança qualitativa é a supressão da expansão exponencial que dava origem ao caos.
  • Observação prática para prova: ao projetar um controlador linear estável, você espera que os expoentes de Lyapunov do sistema controlado sejam todos negativos (se quiser um ponto fixo estável) ou um zero e restantes negativos (se quiser preservar uma órbita periódica estável).

Resposta

A) \lambda_{1}=0{,}125 indica uma direção instável com crescimento exponencial médio (sensibilidade a condições iniciais); \lambda_{2}=0{,}000 é a direção neutra (em sistemas autônomos corresponde à direção tangente ao fluxo, deslocamento temporal); \lambda_{3}=-14{,}60 indica forte contração numa direção transversal. O conjunto mostra expansão em uma direção, uma direção neutra e forte contração na terceira, e soma negativa dos expoentes ( -14{,}475), ou seja, sistema dissipativo com sensibilidade típica de comportamento caótico.
B) Atrator caótico (atrator estranho dissipativo) — pois existe \lambda>0 e a soma é negativa; não é ponto fixo, nem ciclo limite estável, nem toro.
C) Dimensão de Kaplan–Yorke: \displaystyle D_{KY}=2+\frac{0{,}125+0{,}000}{14{,}60}\approx 2{,}00856. Interpretação: dimensão fractal ligeiramente acima de 2, isto é, o atrator é essencialmente uma superfície comprimida com estrutura fractal fina.
D) Com controle linear que elimina o caos, o expoente positivo deixaria de existir: para um ponto fixo estabilizado todos os \lambda_i seriam negativos; para um ciclo limite estabilizado haveria um zero e os demais negativos. Em qualquer caso o \lambda_{1}>0 seria levado a \le 0, suprimindo a sensibilidade exponencial.

Exercícios relacionados