Uma massa m>0 em uma mola com constante elástica k>0 satisfaz a equação diferencial m u^{\prime \prime}+k u=0, onde u(t) é o deslocamento da massa no instante t a partitr de sua posição de equilíbrio.
(a) Sejam x=u e y=u^{\prime}. Mostre que o sistema resultante é:
(b) Encontre os autovalores/autovetores da matriz de coeficientes e escreva a solução geral do sistema.
(c) Esboce as trajetórias do sistema usando o plano de fase e seu campo de inclinações.
(d) Esboce as soluções x, y em função de t num único gráfico e interprete fisicamente este resultado.
Passo 1
O enunciado dá a equação de um oscilador massa-mola
m\,u'' + k\,u = 0
com m>0 e k>0, onde u(t) é o deslocamento em relação ao equilíbrio. Vamos definir as variáveis de estado
x = u,\quad y = u',
de modo a transformar essa EDO de ordem 2 num sistema linear de ordem 1. Em seguida usaremos teoria de autovalores/autovetores para resolver, depois analisaremos o plano de fase e, por fim, veremos o comportamento temporal de x(t) e y(t).
Para \lambda = i\omega um autovetor é \displaystyle v_1 = \begin{pmatrix}1\i\omega\end{pmatrix},
e para \lambda = -i\omega temos \displaystyle v_2 = \begin{pmatrix}1\\\ -\,i\omega\end{pmatrix}.
Assim a solução geral complexa é
E = \frac12\,m\,y^2 + \frac12\,k\,x^2 = \text{constante}
define elipses (ou círculos, se escalonarmos y por \sqrt{m}). O campo de inclinações mostra vetores tangentes que percorrem essas órbitas no sentido anti-horário, caracterizando oscilações periódicas sem amortecimento.
Passo 5
Traçando x(t) e y(t) num mesmo gráfico, vemos que são funções sinuoidais defasadas de \pi/2: quando x(t) está no máximo (\cos máximo), y(t) passa por zero, e vice-versa. Fisicamente, x(t) é o deslocamento e y(t) a velocidade da massa; a energia alterna entre potencial (x máximo, y=0) e cinética (y máximo, x=0), com frequência \omega=\sqrt{k/m}.
Resposta
(a) \displaystyle\binom{x'}{y'}=\begin{pmatrix}0 & 1\\[3pt]-\tfrac{k}{m} & 0\end{pmatrix}\binom{x}{y}.
(b) Os autovalores são \lambda=\pm i\sqrt{k/m}; autovetores \begin{pmatrix}1\\pm i\sqrt{k/m}\end{pmatrix}. A solução geral real é
(c) As trajetórias no plano de fase são elipses (ou círculos após escala) centradas na origem, correspondendo a órbitas periódicas de energia constante, com campo de inclinações dy/dx=-(k/m)(x/y).
(d) Num único gráfico, x(t) e y(t) são senos/cossenos defasados \pi/2, mostrando que a velocidade lidera o deslocamento em \pi/2. Fisicamente, é a oscilação harmônica simples onde energia potencial e cinética se alternam sem perda.]