Enunciado

x ' = - y - z ,       y ' = x + a y ,       z ' = b + z x - c

em que a, b e c são parâmetros positivos, é conhecido como o sistema de Rössler. 20 É um sistema relativamente simples que consiste em duas equações lineares e uma terceira equação com uma única não linearidade quadrática. Nos problemas de 11 a 15, pedimos que você faça algumas investigações numéricas sobre esse sistema com o objetivo de explorar sua propriedade de duplicação de período. Para simplificar, faça a   =   0,25 , b   =   0,5 e deixe c   >   0 permanecer arbitrário.

11. (a) Mostre que não existem pontos críticos quando c   < 0.5 , existe um único ponto crítico quando c   = 0.5 e existem dois pontos críticos quando c   > 0.5 .

(b) Encontre o(s) ponto(s) crítico(s) e determine os autovalores correspondentes da matriz jacobiana quando c   = 0.5 e quando c   =   1 .

(c) Como você acha que será o comportamento das trajetórias do sistema para c   =   1 ? Desenhe a trajetória que começa na origem. Ela se comporta da maneira que você esperava?

(d) Escolha um ou dois outros pontos iniciais e desenhe as trajetórias correspondentes. Esses gráficos estão de acordo com suas expectativas?

Passo 1

L e t r a   A :

A primeira coisa que devemos aqui é substituir as derivadas por 0 , já que estamos buscando os valores dos pontos críticos:

x ' = - y - z ,       y ' = x + 0,25 y ,       z ' = 0,5 + z x - c     →

0 = - y - z ,       0 = x + 0,25 y ,       0 = 0,5 + z x - c     →

Da segunda equação...

x = - 0,25 y  

E da primeira equação, tiramos que...

z = - y

Jogando esses dois últimos resultados na última equação...

0 = 0,5 + z x - c   →

0 = 0,5 - y - 0,25 y - c     →

Poxa! Uma equação do segundo grau:

0,25 y 2 + c y + 0,5 = 0

Sabemos que ela só tem então solução real se Δ for maior ou igual a 0 :

Δ = c 2 - 4 ⋅ 0,25 ⋅ 0,5

Portanto...

c 2 - 4 ⋅ 0,25 ⋅ 0,5 = 0   →

c 2 = 0,5   →

c = 0,5

Então se c < 0,5 , implica no Δ < 0 e, portanto não teremos ponto crítico nenhum. De fato com c = 0,5 teremos um único valor para y e, portanto, para x e z também, logo um único ponto crítico nesse caso.

Agora se c > 0,5 , então teremos dois valores possíveis para y , cada um gera um par de x e z , ou seja, teremos então dois pontos críticos!

Passo 2

L e t r a   B :

PRIMEIRO PARA c = 0,5 :

0,25 y 2 + c y + 0,5 = 0

Sabemos que esse caso é o que nos gera um delta igual a 0 :

Δ = c 2 - 4 ⋅ 0,25 ⋅ 0,5 = 0

Logo, o único valor real de y será...

y = - c 2 ⋅ 0,25 = - 0,5 0,5

Então do passo anterior, podemos tirar as informações de x e z ...

x = - 0,25 y           →               x = 0,5 2

z = - y       →               z = 0,5 0,5

O ponto crítico no caso será...

0,5 2 ,   - 0,5 0,5 ,   0,5 0,5

E quanto a matriz Jacobiana? De acordo com a nossa teoria podemos tirar essa informação da matriz:

J = F x F y F z G x G y G z H x H y H z

Onde...

x ' = - y - z       →           F = - y - z

y ' = x + 0,25 y       →               G = x + 0,25 y

z ' = b + z x - c           →             H = 0,5 + z x - 0,5

Calculando essas derivadas parciais então...

J = 0 - 1 - 1 1 0,25 0 z 0 x - 0,5   →

Logo, o valor do nosso Jacobiano será no nosso único ponto crítico no caso será...

J 0,5 2 ,   - 0,5 0,5 ,   0,5 0,5 = 0 - 1 - 1 1 0,25 0 0,5 0,5 0 0,5 2 - 0,5     →

Simplificando um pouquinho...

J 0,5 2 ,   - 0,5 0,5 ,   0,5 0,5 = 0 - 1 - 1 1 0,25 0 0,5 0,5 0 - 1 2 0,5    

Para extrair as informações de autovalores dessa matriz, podemos usar o auxilio computacional, já que o problema nos permite tal acesso . Aqui, vou utilizar o WOLFRAMALPHA, afinal ele é uma plataforma muito simples e eficaz. Basta você escrever a nossa matriz lá na barrinha de fórmula do Site e dar um enter, ganhando:

λ = 0         e           λ = - 0,0517767 ± 1.52419 i

Sacou? Então temos um autovalor nulo e outros dois imaginários!

Passo 3

PRIMEIRO PARA c = 1 :

Refazendo as continhas agora para c = 1 ...

0,25 y 2 + c y + 0,5 = 0   →

0,25 y 2 + y + 0,5 = 0   →

Δ = 1 - 4 ⋅ 0,25 ⋅ 0,5 = 0,5

y = - 1 ± 0,5 2 ⋅ 0,25           →               y 1 = - 1 + 0,5 0,5 y 2 = - 1 - 0,5 0,5

Dessa vez obtivemos dois valores para y , então é de se esperar dois valores diferentes para x e z ...

x = - 0,25 y           →                 x 1 = - - 1 + 0,5 2 x 2 = - - 1 - 0,5 2

z = - y       →               z 1 = - - 1 + 0,5 0,5 z 2 = - - 1 - 0,5 0,5

Os dois pontos críticos nesse caso então serão...

- - 1 + 0,5 2 ,   - 1 + 0,5 0,5 ,   - - 1 + 0,5 0,5             e                 - - 1 - 0,5 2 ,   - 1 - 0,5 0,5 ,   - - 1 - 0,5 0,5

O jacobiano desse primeiro ponto será então...

J = 0 - 1 - 1 1 0,25 0 z 0 x - 0,5   →

J - - 1 + 0,5 2 ,   - 1 + 0,5 0,5 ,   - - 1 + 0,5 0,5 = 0 - 1 - 1 1 0,25 0 - - 1 + 0,5 0,5 0 - - 1 + 0,5 2 - 0,5     →

Jogando no Wolfram para obter seus autovalores ganharemos...

λ = - 0.287838         e           λ = - 0,0 114112 ± 1.19955 i

MAS CALMA, meu Jovem! Temos um segundo ponto crítico e que certamente precisa ser avaliado também!

J - - 1 - 0,5 2 ,   - 1 - 0,5 0,5 ,   - - 1 - 0,5 0,5 = 0 - 1 - 1 1 0,25 0 - - 1 - 0,5 0,5 0 - - 1 - 0,5 2 - 0,5     →

Calculando os seus autovalores por sua vez obteremos...

λ = 0.226637         e           λ = 0.0849048 ± 2.09884   i

Passo 4

L e t r a   C :

Bem, como no caso de c = 1 , obtivemos dois pontos críticos diferentes, onde para um deles...

- - 1 + 0,5 2 ,   - 1 + 0,5 0,5 ,   - - 1 + 0,5 0,5

Tivemos autovalores imaginários com a parte real negativa e a outra parte real positiva, é de se esperar que pontos próximos deste ponto referenciado aí em cima sejam assintoticamente estáveis, enquanto para regiões mais próximas do outro ponto seja instável.

No caso da origem, por ela estar mais próximo de:

- - 1 + 0,5 2 ,   - 1 + 0,5 0,5 ,   - - 1 + 0,5 0,5

É esperado que ela entre em estabilidade assintótica!

De fato, usando o mesmo código que utilizamos no exercício 6, mas adaptando para este exercício dessa forma aqui:

Onde a função de definição da EDO foi feita aqui separadamente:

Bem, mas plotando a solução para o ponto inicial de 0 ,   0 ,   0 ...

Viu só! De fato a nossa observação está certa! A solução tende a estabilidade ;)

Passo 5

L e t r a   D :

Para um ponto inicial 1 ,   1 ,   1 ...

Temos a convergência.

Para um ponto - 5 ,   - 5 ,   - 5 ...

Temos uma divergência, bem como previmos!

Resposta

Letra A: Demonstração.

Letra B:

P a r a   c = 0.5       →         0,5 2 ,   - 0,5 0,5 ,   0,5 0,5

P a r a   c = 1       →

- - 1 + 0,5 2 ,   - 1 + 0,5 0,5 ,   - - 1 + 0,5 0,5               e                   - - 1 - 0,5 2 ,   - 1 - 0,5 0,5 ,   - - 1 - 0,5 0,5

Letra C:

Letra D: Ponto 1 ,   1 ,   1 :

Ponto - 5 ,   - 5 ,   - 5 :

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

1.(a) Mostre que os autovalores do sistema linear (8), válidos perto da origem, são dados pela Eq. (10).(b) Determine os autovetores associados.(c) Determine os autovalores e autovetores do sistema (8
7. Para r=28 , faça as projeções nos planos xy e xz, respectivamente, da trajetória que começa no ponto ( 5 ,5 ,5 ) . Os gráficos são iguais aos das Figuras 9.8.5 e 9.8.6?
2. (a) Mostre que a aproximação linear válida perto do ponto crítico P 2 é dada pela Eq. (11).(b) Mostre que os autovalores do sistema (11) satisfazem a Eq. (12).(c) Para r=28 , resolva a Eq. (12) e d
Em cada um dos problemas de 6 a 10, faça as investigações pedidas sobre as equações de Lorenz.6. Para r=28 , faça o gráfico de x em função de t para os casos ilustrados nas Figuras 9.8.2 e 9.8.3. Seus
4. Use a função de Liapunov V x ,y ,z = x 2 + σ y 2 + σ z 2 para mostrar que a origem é um ponto crítico global assintoticamente estável para as equações de Lorenz (1) se r1 .