Pêndulos Estocásticos

Grupo : Gustavo H. Guesser, Joshua L. Kipper, Marcos Pasa.

Pêndulo Simples

Equação de movimento

Primeiramente vamos inserir ruído em um pêndulo simples, que é constituído de uma barra de comprimento $l$, sem massa e rígida que contém uma massa $m$ pontual em sua ponta, conforme ilustrado na figura a seguir.

Esquema de um pêndulo simples em um campo gravitacional constante.

Considerando que o pêndulo está sob o efeito da gravidade e se encontra submerso em um fluido viscoso (como o ar), tal que a força de resistência que atua na massa é $-2b{\vec {v}}$, a equação de movimento é dada por:

${\ddot {\theta }}(t)=-{\frac {2b}{m}}{\dot {\theta }}-{\frac {g}{l}}sen(\theta )$

Vamos supor que existe uma força ruidosa atuando em Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m} , tal que sua componente tangencial (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle F_{r}(t)} ) pode ser modelada por um ruído branco gaussiano Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi(t)} da seguinte forma

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle F_r(t) = m \alpha \xi(t) }

em que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} é a intensidade do ruído. Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi(t)} é caracterizado pelas seguintes propriedades:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle \xi(t) \rangle = 0,~~ \langle \xi(t_2)\xi(t_1) \rangle = \delta(t_2 -t_1) }

Adicionando essa nova força nas equações de movimento, ficamos com

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta(t) = -\frac{2b}{m}\dot \theta - \frac{g}{l}sen(\theta) + \frac{\alpha}{l}\xi(t) }

A partir de agora, por questão de simplicidade, vamos supor que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle g = l = 1 } , então

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta(t) = -2b\dot \theta - sen(\theta) + \alpha\xi(t) }

Método de integração

Vamos montar um método para integrar o sistema no tempo. Primeiramente vamos dividir a equação em duas equações diferencias de primeira ordem, introduzindo a variável Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega := \dot \theta } , então ficamos com o seguinte sistema

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} \dot \theta &= \omega \\ \dot \omega &= -2b \omega - sen(\theta) + \alpha\xi(t) \end{aligned} }

que pode ser escrito na forma diferencial

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} d\theta &= \omega dt \\ d\omega &= (-2b \omega - sen(\theta))dt + \alpha\xi(t)dt \end{aligned} }

mas Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi(t)dt} é o incremento do processo de Wiener (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle W(t) = \int_0^t \xi(t')dt' } ), então

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} d\theta &= \omega dt \\ d\omega &= (-2b\omega - sen(\theta))dt + \alpha dW(t) \end{aligned} }

Discretizando o tempo e lembrando que a densidade de probabilidade de transição de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle W(t)} para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle W(t+\Delta t)} tem desvio padrão igual a Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sqrt{\Delta t}}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} \theta_{j+1} &= \theta_{j} + \omega_j \Delta t \\ \omega_{j+1} &= \omega_{j} + (-2b \omega_j - sen(\theta_j))\Delta t + \alpha {R_G}_j \sqrt{\Delta t} \end{aligned} }

em que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle R_G} é uma amostra de uma distribuição gaussiana com média 0 e variância 1, e o método de Euler foi utilizado para a parte determinística da equação.

Nas próximas seções será analisado a energia do sistema, e como o método de Euler não é muito bom para preservar a energia de sistemas conservativos, será utilizado o método preditor corretor (com adição do método de Heun para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta } ) para a parte determinística da equação, que consiste nos seguintes passos:

• Calcular um theta intermediário:

   Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle  \theta^{(2)}_{j+1} = \theta_j + \omega_j \Delta t }

• Com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta^{(2)}_{j+1} } calcular um theta médio e utilizá-lo para obter um omega intermediário:

   Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle  \begin{aligned} \bar \theta_j &= \frac{\theta^{(2)}_{j+1} + \theta_j}{2} \\ \omega^{(2)}_{j+1} &= \omega_j + f(\bar \theta_j, \omega_j, {R_G}_j) \end{aligned} }

Em que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle f} é a expressão do método de Euler visto logo acima.

• Recalcular theta utilizando um omega intermediário

   Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle  \begin{aligned} \bar \omega_j &= \frac{\omega^{(2)}_{j+1} + \omega_j}{2} \\ \theta_{j+1} &= \theta_j + \bar \omega_j \Delta t  \end{aligned} }

• Recalcular omega com um theta intermediário atualizado

   Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle  \begin{aligned} \bar \theta^{(2)}_j &= \frac{\theta_{j+1} + \theta_j}{2} \\ \omega_{j+1} &= \omega_j + f(\bar \theta^{(2)}_j, \bar \omega_j, {R_G}_j) \end{aligned} }

OBS: No cálculo de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega^{(2)}_{j+1}} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_{j+1}} foi utilizado o mesmo Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle {R_G}_j} .

Energia (Sem amortecimento)

Logo após terminar a implementação do método numérico, rapidamente notamos que a adição do ruído gera um aumento na energia mecânica do pêndulo (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E} ), vamos explorar esse fenômeno, sempre utilizando Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta = 0} . Para ilustrar esse efeito, segue uma animação do pêndulo partindo do repouso na configuração de equilíbrio estável (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta(0) = 0,~ \omega(0)=0} ) com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha = 0.1}

Pêndulo partindo do repouso com ruído.

Para realizar uma exploração quantitativa, o seguinte procedimento foi feito para vários valores de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha } :

• Utilizando Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta t = 0.01} , integrar o sistema até Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle t_f = 100} , calculando e armazenando a energia em cada passo temporal. Repetir essa integração 700 vezes para fazer médias temporais. Como exemplo, segue os dados obtidos da energia em função do tempo (com a média temporal feita) para um determinado Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha } utilizado

Energia média em função tempo.

O gráfico nos indica que o ruído gera uma potência média sobre o pêndulo.

• Realizar um ajuste linear nos dados Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle E \rangle_t \times t} para obter o coeficiente angular, que corresponde a potência média gerada pelo ruído (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar P} ).

Com as simulações executadas, foi realizado o gráfico Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar P \times \alpha} . Notamos que os dados se alinham em linha reta com os eixos em escala logarítmica, ou seja, os mesmos seguem uma lei de potência Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar P = ae^{b \alpha}} , então foi realizado outro ajuste linear para encontra o expoente Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle b} , a figura a seguir ilustra os dados e os resultados do ajuste:

Potência em função do ruído (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} ). O painel da esquerda possui eixos em escala linear e o da direita em escala logarítmica.

Portanto, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar P} aumenta, aproximadamente, de forma quadrática com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} . Por fim, gostaríamos de mostrar que mesmo para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} muito pequeno, ainda existe energia sendo injetada no sistema com taxa constante, e isso não é um artefato da simulação. Para tal, foram realizadas 700 simulações com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha=5\cdot10^{-4}} e os resultados foram comparados com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha = 0}

Energia média para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} muito pequeno comparado com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} nulo.

Energia (Com amortecimento)

Até o momento, o amortecimento foi negligenciado. Vamos, então, introduzi-lo e rodar várias simulações (700, neste caso) e observar como a média temporal da energia evolui ao longo do tempo, assim como foi feita na seção anterior. Utilizando Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta =0.1} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha = 0.485} foi obtido o seguinte resultado

Energia mecânica média em função do tempo com amortecimento.

claramente o comportamento neste caso é diferente do observado sem amortecimento, agora a anergia aumenta até um certo valor e permanece nele. Para explorar este novo fenômeno, os seguintes passos foram feitos para cada valor de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta = 0.01, 0.1, 0.2 } :

• Para diversos valores de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} , executar 700 simulações até a energia estabilizar, salvando a média da energia entre as simulações.
• Para cada conjunto de dados gerados por um determinado Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} , selecionar um intervalo de tempo onde a energia está estabilizada e calcular a sua média (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle E \rangle} ).

Produzindo o gráfico de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle E \rangle \times \alpha} obtemos

Energia estabilizada média em função de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} para diferentes valores de amortecimento. Os painéis da esquerda possui eixos em escala linear e os da direita em escala logarítmica.

as linhas vermelhas são os melhores ajustas de leis de potência na forma Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle E \rangle = a \alpha^b } . Para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta=0.1} os dados utilizados no ajuste foram apenas até Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha = 0.6} (indicado pela reta preta vertical no gráfico), pois após esse limite, a lei de potência deixa de ser um ótimo ajuste. É chamativo o fato de todos os coeficientes, independente de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta} , serem aproximadamente 2.

Pêndulo invertido

O próximo pêndulo a ser considerado é um pêndulo invertido, que possui um potencial harmônico em seu ângulo, ou seja, um dos termos de sua energia potential é Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{k}{2}\theta^2 } , sendo que agora Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta} é zero quando a haste está apontando para cima, conforme ilustrado na figura a seguir

Esquema do pêndulo invertido com movimento vertical livre na base.

Ainda, a base do pêndulo é livre para movimentar-se na direção vertical, é justamente nesse local onde será adicionado uma força ruidosa. A equação de movimento neste caso é

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta = \frac{-2b}{m}\dot \theta + \frac{1}{l}(g + \ddot u)sen(\theta) - \frac{k}{m}\theta }

O primeiro termo vem da resistência do ar, o segundo se origina da gravidade e do deslocamento de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle u} e o último provém da "mola" em Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta } . Vamos supor que existe um vínculo mantendo a base fixa, mas o local onde o pêndulo é fixado pode se movimentar de forma aleatória na direção vertical, supondo que tal movimento pode ser modelado por ruído branco gaussiano, segue que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot u(t) = \alpha \xi(t)} . Introduzindo a variável Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega := \dot \theta} , ficamos com os seguintes sistemas de equações na forma diferencial

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} d\theta &= \omega dt \\ d\omega &= (\frac{-2b}{m}\omega + \frac{g}{l}sen(\theta) - \frac{k}{m}\theta )dt + sen(\theta) \frac{\alpha}{l}\overbrace{dW(t)}^{\xi(t)dt} \end{aligned} }

Note que agora o ruído é multiplicativo, em contraste com o ruído aditivo dos pêndulos anteriores, para lidar com esta complicação, no momento da integração vamos utilizar um Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta} médio no argumento do seno que multiplica Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle dW}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{aligned} \theta_{j+1} &= \theta_j + \omega_j \Delta t \\ \bar \theta_j &= (\theta_{j+1} + \theta_{j} ) / 2 \\ \omega_{j+1} &= \omega_j + (\frac{-2b}{m}\omega_j + \frac{g}{l}sen(\bar \theta_j) - \frac{k}{m} \bar \theta_j ) \Delta t + sen(\bar \theta_j) \frac{\alpha}{l} {R_G}_j \sqrt{\Delta t} \end{aligned} }

Esse tipo de pêndulo é de grande interesse em algumas áreas, como na engenharia estrutural, pois uma coluna comprimida pode ser modelada como um pêndulo invertido com uma mola na base, o ruído da base pode representar um terremoto.

Uma questão interessante neste modelo é sobre a estabilidade da configuração Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta_0 = \omega_0 = 0 } quando é adicionado ruído e amortecimento. Se Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle k} é grande o suficiente, sem ruído, o equilíbrio é estável, mas é possível sair dessa configuração adicionando ruído, cujo valor limite vai depender do amortecimento. Podemos demostrar isso integrando o sistema, partindo da condição inicial Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta_0 = 10^{-4}, \omega_0 = 0 } , com valores de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha} muito próximos:

os seguintes valores foram utilzados

• g = l = 1
• k = 1.1

No gráfico da esquerda, o ângulo oscilou um pouco e permaneceu em 0, já no outro, eventualmente, o ângulo explodiu. Abaixo segue uma animação dessa situação, mas com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha=1} para o pêndulo sair do equilíbrio mais rápido

Pêndulo invertido partindo do equilíbrio com ruído.

Pêndulo Duplo Estocástico

O pêndulo duplo estocástico é um sistema dinâmico que combina a complexidade intínseca do pêndulo duplo com a introdução de elementos de aleatoriedade ou incerteza, tornando o comportamento do sistema ainda mais imprevisível e caótico. O pêndulo duplo em si é um exemplo clássico de um sistema caótico, onde pequenas variações nas condições iniciais podem resultar em trajetórias drasticamente diferentes. Quando um termo estocástico é adicionado, por exemplo, na forma de uma força externa aleatória ou de flutuações nos parâmetros do sistema, a análise e a previsão do movimento se tornam desafiadoras.

Pêndulo Duplo não estocástico.

Equação de movimento

O pêndulo duplo consiste em dois pêndulos acoplados, onde o segundo pêndulo está suspenso na extremidade do primeiro. As equações de movimento para o pêndulo duplo sem termos estocásticos podem ser derivadas utilizando as equações de Lagrange, considerando as coordenadas angulares Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta_1 } e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta_2 } como as variáveis generalizadas e são dadas por:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_1 = \frac{-g(2m_1 + m_2)\sin{\theta_1} - m_2 g(\theta_1 - 2 \theta_2) - 2 \sin{(\theta_1 - \theta_2)}m_2(\dot {\theta_2}^2 l_2 + \dot {\theta_1}^2 l_1 \cos{(\theta_1 - \theta_2)})}{l_1 (2 m_1 + m_2 - m_2 \cos{(2 \theta_1 - 2 \theta_2)})} }

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_2 = \frac{2 \sin{(\theta_1 - \theta_2)} (\dot {\theta_1}^2 l_1 (m_1 + m_2) + g (m_1 + m_2) \cos{\theta_1} + \dot {\theta_2}^2 l_2 m_2 \cos{(\theta_1 - \theta_2)})}{l_2 (2 m_1 + m_2 - m_2 \cos{(2 \theta_1 - 2 \theta_2)})} }

mantendo Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_{1} } , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m_{2} } , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle l } , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle g } , fixos, defina :

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle f_{1}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) := \frac{-g(2m_1 + m_2)\sin{\theta_1} - m_2 g(\theta_1 - 2 \theta_2) - 2 \sin{(\theta_1 - \theta_2)}m_2(\dot {\theta_2}^2 l_2 + \dot {\theta_1}^2 l_1 \cos{(\theta_1 - \theta_2)})}{l_1 (2 m_1 + m_2 - m_2 \cos{(2 \theta_1 - 2 \theta_2)})} }

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle f_{2}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) := \frac{2 \sin{(\theta_1 - \theta_2)} (\dot {\theta_1}^2 l_1 (m_1 + m_2) + g (m_1 + m_2) \cos{\theta_1} + \dot {\theta_2}^2 l_2 m_2 \cos{(\theta_1 - \theta_2)})}{l_2 (2 m_1 + m_2 - m_2 \cos{(2 \theta_1 - 2 \theta_2)})} }

assim as equações dinâmicas ficam escritas de maneira mais compacta:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_1 = f_{1}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) }

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_2 = f_{2}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) }

suponha que haja um força externa ruidosa, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle F_r = \sigma \xi (t) } , onde Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma } além de dar a intensidade do ruído branco, terá consigo todas as constantes agrupadas. Assim as equações de movimento se tornam:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_1 = f_{1}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) + \sigma_{1} \xi_{1} (t) }

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ddot \theta_2 = f_{2}(\theta_{1}, \dot {\theta_{1}}, \theta_{2}, \dot {\theta_{1}}) + \sigma_{2} \xi_{2} (t) }

Retrato de fase

Simulação de um pêndulo duplo estocástico com ruído branco em Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta_1 } .

Retrato de fase.