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).
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)
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.
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.