View in NBViewer Open in binder Download notebook View source


7.5. Hamiltonianos

Obtendo o Hamiltoniano no caso de um grau de liberdade

\[ \begin{cases} \displaystyle \dot q = \frac{\partial H}{\partial p}(q,p), \\ \displaystyle \dot p = - \frac{\partial H}{\partial q}(q,p), \end{cases} \] \[ \frac{\mathrm{d}}{\mathrm{d} t} H(q(t),p(t)) = \frac{\partial H}{\partial q} \dot q + \frac{\partial H}{\partial p}\dot p = \frac{\partial H}{\partial q} \left( \frac{\partial H}{\partial p}\right) + \frac{\partial H}{\partial p}\left(- \frac{\partial H}{\partial q}\right) = 0. \]

Forma compacta de se escrever o sistema Hamiltoniano

\[ \nabla H(q,p) = \left( \frac{\partial H}{\partial q}(q,p), \frac{\partial H}{\partial p}(q,p) \right), \qquad J = \left[ \begin{matrix} 0 & 1 \\ -1 & 0 \end{matrix} \right]. \] \[ (\dot q,\dot p) = J \nabla H(q,p). \]

Hamiltoniano a partir da integral de Jacobi

\[ h(q,\dot q) = \dot q \frac{\partial L}{\partial \dot q}(q,\dot q)- L(q,\dot q). \] \[ p = \frac{\partial L}{\partial \dot q}(q,\dot q) \] \[ \dot q = \dot q(q,p). \] \[ H(q,p) = p\dot q(q,p) - L(q,\dot q(q,p)). \] \[ \begin{align*} \frac{\partial H}{\partial p}(q,p) & = \dot q(q,p) + p\frac{\partial \dot q}{\partial p}(q,p) - \frac{\partial L}{\partial \dot q}(q,\dot q(q,p))\frac{\partial \dot q}{\partial p}(q,p) \\ & = \dot q + p\frac{\partial \dot q}{\partial p}(q,p) - p\frac{\partial \dot q}{\partial p}(q,p) \\ & = \dot q. \end{align*} \] \[ \begin{align*} \frac{\partial H}{\partial q}(q,p) & = p\frac{\partial \dot q}{\partial q}(q,p) - \frac{\partial L}{\partial q}(q,\dot q(q,p)) - \frac{\partial L}{\partial \dot q}(q,\dot q(q,p))\frac{\partial \dot q}{\partial q}(q,p) \\ & = p\frac{\partial \dot q}{\partial q}(q,p) - \frac{\partial L}{\partial q}(q,\dot q(q,p)) - p\frac{\partial \dot q}{\partial q}(q,p) \\ & = - \frac{\partial L}{\partial q}(q,\dot q(q,p)) \\ & = - \frac{\mathrm{d}}{\mathrm{d} t} \frac{\partial L}{\partial \dot q}(q,\dot q(q,p)) \\ & = - \frac{\mathrm{d}}{\mathrm{d} t} p \\ & = - \dot p. \end{align*} \] \[ \begin{cases} \displaystyle \dot q = \frac{\partial H}{\partial p}(q,p), \\ \displaystyle \dot p = - \frac{\partial H}{\partial q}(q,p). \end{cases} \]

Exemplo do pêndulo planar

\[ L(\theta,\dot\theta) = \frac{1}{2}m\ell^2\dot\theta^2 + mg\ell\cos\theta. \]

caso é o momento angular do sistema, e é dado por

\[ \psi = \frac{\partial L}{\partial \dot\theta}(\theta,\dot\theta) = m\ell^2\dot\theta \] \[ \dot\theta = \dot\theta(\psi) = \frac{1}{m\ell^2}\psi. \] \[ h(\theta,\dot\theta) = \dot\theta\psi - L(\theta,\dot\theta). \] \[ H(\theta,\psi) = \frac{1}{2m\ell^2}\psi^2 -mg\ell\cos\theta. \] \[ \begin{cases} \displaystyle \dot \theta = \frac{1}{m\ell^2}\psi, \\ \displaystyle \dot \psi = -mg\ell\sin\theta. \end{cases} \]

Obtendo o Hamiltoniano no caso de vários graus de liberdade

\[ L(\mathbf{q},\mathbf{\dot q}). \] \[ h(\mathbf{q},\mathbf{\dot q}) = \mathbf{\dot q}\cdot\mathbf{p} - L(\mathbf{q},\mathbf{\dot q}), \qquad \mathbf{p} = \nabla_{\mathbf{\dot q}} L(\mathbf{q},\mathbf{\dot q}) \] \[ \mathbf{\dot q} = \mathbf{\dot q}(\mathbf{q},\mathbf{p}). \] \[ H(\mathbf{q},\mathbf{p}) =\mathbf{p}\cdot\mathbf{\dot q}(\mathbf{q},\mathbf{p}) - L(\mathbf{q},V(\mathbf{q},\mathbf{p}). \] \[ \begin{cases} \displaystyle\mathbf{\dot q} = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p}), \\ \displaystyle \mathbf{\dot p} = - \frac{\partial H}{\partial \mathbf{q}}(\mathbf{q},\mathbf{p}), \end{cases} \] \[ \begin{align*} \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p}) & = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p}) = \left(\frac{\partial H}{\partial p_1}(\mathbf{q},\mathbf{p}), \ldots, \frac{\partial H}{\partial p_d}(\mathbf{q},\mathbf{p})\right), \\ \frac{\partial H}{\partial \mathbf{q}}(\mathbf{q},\mathbf{p}) & = \frac{\partial H}{\partial\mathbf{q}}(\mathbf{q},\mathbf{p}) = \left(\frac{\partial H}{\partial q_1}(\mathbf{q},\mathbf{p}), \ldots, \frac{\partial H}{\partial q_d}(\mathbf{q},\mathbf{p})\right). \end{align*} \]

Obtendo o Lagrangiano a partir do Hamiltoniano

\[ H(\mathbf{q},\mathbf{p}) = \mathbf{\dot q}\cdot\mathbf{p} - L(\mathbf{q},\mathbf{\dot q}), \qquad \mathbf{\dot q} = \mathbf{\dot q}(\mathbf{q},\mathbf{p}), \] \[\mathbf{p} = \nabla_{\mathbf{\dot q}} L(\mathbf{q},\mathbf{\dot q}). \] \[ L(\mathbf{q},\mathbf{\dot q}) = \mathbf{\dot q}\cdot\mathbf{p} - H(\mathbf{q},\mathbf{p}), \qquad\mathbf{p}=\mathbf{P}(\mathbf{q},\mathbf{\dot q}) \] \[ \mathbf{\dot p} = - \frac{\partial H}{\partial\mathbf{q}}(\mathbf{q},\mathbf{p}). \]

Obtendo o Lagrangiano do Hamiltoniano no caso do pêndulo planar

\[ H(\theta,\psi) = \frac{1}{2m\ell^2}\psi^2 - mg\ell\cos\theta. \] \[ \dot\theta = \frac{\partial H}{\partial \psi}(\theta,\psi) = \frac{1}{m\ell^2}\psi, \] \[ \psi = m\ell^2\dot\theta. \] \[ \begin{align*} L(\theta,\dot\theta) & = \dot\theta\psi - H(\theta,\psi) \\ & = m\ell^2 \dot\theta^2 - \left(\frac{1}{2m\ell^2}(m\ell^2\dot\theta)^2 - mg\ell\cos\theta\right) \\ & = \frac{1}{2}m\ell^2\dot\theta^2 + mg\ell\cos\theta, \end{align*} \]

Sistemas não-autônomos

\[ L(\mathbf{q},\mathbf{\dot q},t). \] \[ \begin{cases} \displaystyle\mathbf{\dot q} = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p},t), \\ \displaystyle \mathbf{\dot p} = - \frac{\partial H}{\partial \mathbf{q}}(\mathbf{q},\mathbf{p},t), \end{cases} \] \[ H(\mathbf{q},\mathbf{p},t). \] \[ \frac{\mathrm{d} }{\mathrm{d} t} H(\mathbf{q}(t),\mathbf{p}(t),t) = \frac{\partial H}{\partial t}(\mathbf{q}(t),\mathbf{p}(t),t). \]

autônomo.

Transformada de Legendre

\[ \mathbf{q}=q\in \mathbb{R}, \qquad \mathbf{\dot q} = \dot q, \qquad \mathbf{p} = p = \frac{\partial L}{\partial \dot q}(q,\dot q). \]

A transformada de Legendre relacionando Lagrangiano e Hamiltoniano

\[ H(p) = pq - L(\dot q), \qquad \text{ onde } \dot q = \dot q(p) \text{ é a inversa de } p=\frac{\partial L(\dot q)}{\partial \dot q}. \] \[ g^*(s) = sr - g(r), \qquad \text{ onde } r=r(s) \text{ é a inversa de } s=g'(r). \] \[ g^*(s) = \max_{r\in \mathbb{R}} \left\{ sr - g(r) \right\}. \]

Reescrevendo a transformada de Legendre

\[ g^*(s) = \max_{r\in \mathbb{R}} \left\{ sr - g(r) \right\}. \] \[ \frac{\partial}{\partial r} \left\{ sr - g(r) \right\} = s - g'(r) = 0, \] \[ \frac{g(r)}{|r|} \rightarrow \infty, \qquad \text{quando } |r|\rightarrow \infty. \]

Exercícios

  1. Ache o Hamiltoniano e as equações de Hamilton nos seguintes casos:

    1. Corpo em queda livre, com o seguinte Lagrangiano, onde \(m,g>0\):

\[ L(h,\dot h) = \frac{1}{2}m\dot h^2 + mgh. \]

1. Corpo deslizando em um plano inclinado, com o seguinte Lagrangiano, onde \(m,g>0\), \(\alpha\in \mathbb{R}\):

\[ L(x,\dot x) = \frac{1}{2}m\dot x^2\sec^2\alpha + mg x\tan\alpha, \]

1. Sistema massa-mola, com o seguinte Lagrangiano, onde \(m,k>0\):

\[ L(x,\dot x) = \frac{1}{2} m \dot x^2 + \frac{1}{2} k x^2, \]
  1. Ache o Hamiltoniano no caso de um pêndulo girante, cujo Lagrangiano é

\[ L(\varphi,\dot\varphi) = \frac{1}{2}m\ell^2(\dot\varphi^2 + \omega\sin^2\varphi) + mg\ell\cos\varphi, \]

onde \(m,g,\ell>0\), \(\omega\in \mathbb{R}\), e verifique que as equações de Hamilton desse sistema tomam a forma

\[ \begin{cases} \displaystyle \dot \varphi = \frac{1}{m\ell^2}\psi, \\ \displaystyle \dot \psi = 2\omega^2\sin\varphi\cos\varphi --mg\ell\sin\varphi. \end{cases} \]

Exercícios - parte 1

  1. Analogamente ao caso bidimensional, ache uma matriz antisimétrica \(\mathbf{J}\in \mathbb{R}^{2d \times 2d}\) tal que um sistema Hamiltoniano em \(\mathbb{R}^{2d}\) pode ser escrito na forma

\[ (\mathbf{\dot q}, \mathbf{\dot p}) = \mathbf{J} \boldsymbol{\nabla} H(\mathbf{q},\mathbf{p}), \]

onde \(\boldsymbol{\nabla} H(\mathbf{q},\mathbf{p}) = (\boldsymbol{\nabla}_\mathbf{q} H(\mathbf{q},\mathbf{p}),\boldsymbol{\nabla}_\mathbf{p} H(\mathbf{q},\mathbf{p}))\) é o gradiente de \(H(\mathbf{q},\mathbf{p})\).

  1. Verifique, nesse caso de vários graus de liberdade, que as soluções da equações de Euler-Lagrange de um Lagrangiano \(L(\mathbf{q},\mathbf{\dot q})\) satisfazem as equações de Hamilton

\[ \begin{cases} \displaystyle\mathbf{\dot q} = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p}), \\ \displaystyle \mathbf{\dot p} = - \frac{\partial H}{\partial \mathbf{q}}(\mathbf{q},\mathbf{p}), \end{cases} \]

para o Hamiltoniano \(H(\mathbf{q},\mathbf{p})\) correspondente.

  1. Considere o problema de três corpos planar restrito circular, em que um corpo celeste \(L\) com massa \(M_L>0\) gira em movimento circular com freqüência constante \(\omega>0\) em torno de um outro corpo celeste \(T\) de massa \(M_T>0\), e um satélite de massa \(m\ll M_L,M_T\) muito menor que a dos dois corpos se movimenta sob a atração gravitacional desses corpos. Considerando um referencial girante centrado em \(T\) e acompanhando a rotação de \(L\) em torno de \(T\), o Lagrangiano para o movimento do satélite pode ser escrito na forma polar

\[ L(r,\theta,\dot r,\dot \theta) = \frac{1}{2} m \dot r^2 + \frac{1}{2} m r^2 (\omega+\dot\theta) + G \frac{mM_T}{r} + G\frac{mM_L}{\sqrt{r^2+R^2 -2rR}}. \]

Obtenha o Hamiltoniano desse sistema e verifique que o sistema hamiltoniano associado pode ser escrito na forma

\[ \begin{cases} \displaystyle \dot r = \frac{p_r}{m}, \\ \displaystyle \dot \theta = \frac{p_\theta}{mr^2} - \omega, \\ \displaystyle \dot p_r = \frac{p_\theta}{r} - G\frac{mM_T}{r^2} - 2G(r-R)\frac{mM_L}{\sqrt{r^2+R^2 -2rR}}, \\ \displaystyle \dot p_\theta = 0, \end{cases} \]

onde \(p_r\) e \(p_\theta\) são os momentos conjugados às coordenadas generalizadas \(r\) e \(\theta\), respectivamente.

  1. Ache o Hamiltoniano e as equações de Hamilton de um problema

de dois pêndulos co-planares com a mesma massa \(m>0\) e hastes de mesmo comprimento \(\ell>0\).

Exercícios - parte 2

  1. Nos problemas das seções anteriores, faça o caminho inverso, obtendo o Lagrangiano a partir do Hamiltoniano, conforme descrito acima.

Exercícios - parte 3

  1. Mostre, nesse caso não-autônomo, que o Hamiltoniano obtido da relação

\[ H(\mathbf{q},\mathbf{p},t) = \mathbf{\dot q}\cdot\mathbf{p} - L(\mathbf{q},\mathbf{\dot q},t), \qquad \mathbf{\dot q} =\mathbf{\dot q}(\mathbf{q},\mathbf{p},t), \]

onde \(\mathbf{\dot q} = \mathbf{\dot q}(\mathbf{q},\mathbf{p},t)\) é solução de

\[\mathbf{p} = \nabla_{\mathbf{\dot q}} L(\mathbf{q},\mathbf{\dot q},t), \]

leva as equações de Euler-Lagrange no sistema hamiltoniano

\[ \begin{cases} \displaystyle\mathbf{\dot q} = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p},t), \\ \displaystyle \mathbf{\dot p} = - \frac{\partial H}{\partial \mathbf{q}}(\mathbf{q},\mathbf{p},t), \end{cases} \]
  1. Mostre que nesse caso não-autônomo o Hamiltoniano satisfaz

\[ \frac{\mathrm{d} }{\mathrm{d} t} H(\mathbf{q}(t),\mathbf{p}(t),t) = \frac{\partial H}{\partial t}(\mathbf{q}(t),\mathbf{p}(t),t). \]
  1. Mostre o Lagrangiano \(L(\mathbf{q},\mathbf{\dot q},t)\) associado a um Hamiltoniano

\(H(\mathbf{q},\mathbf{p},t)\) pode ser recuperado através das relações

\[ L(\mathbf{q},\mathbf{\dot q}) = \mathbf{\dot q}\cdot\mathbf{p} - H(\mathbf{q},\mathbf{p}), \qquad\mathbf{p}=\mathbf{P}(\mathbf{q},\mathbf{\dot q}) \]

onde \(\mathbf{p} = \mathbf{P}(\mathbf{q},\mathbf{\dot q})\) é solução da primeira equação do sistema Hamiltoniano,

\[ \mathbf{\dot q} = \frac{\partial H}{\partial\mathbf{p}}(\mathbf{q},\mathbf{p}). \]
  1. Considere o Lagrangiano

\[ L(\theta,\dot\theta,t) = m\ell\omega\dot\theta\cos\omega t\cos\theta + \frac{1}{2}m\ell^2\dot\theta^2 + mg\ell\cos\theta \]

associado a um pêndulo planar cuja haste fixa está sujeita a um movimento periódico com freqüência \(\omega\in \mathbb{R}\). Obtenha o Hamiltoniano correspondente.

Exercícios - parte 4

  1. Seja \(g:\mathbb{R}\rightarrow \mathbb{R}\) contínua com crescimento superlinear e considere a sua transformada de Legendre

\[ g^*(s)= \max_{r\in \mathbb{R}}\{sr - g(r)\}. \]

a) Mostre que \(g^*:\mathbb{R}\rightarrow \mathbb{R}\) está bem definida. b) Mostre que \(g^*\) é contínua em \(\mathbb{R}\). c) Mostre que \(g^*\) é coerciva. d) Mostre que se \(g\) é convexa, então \(g^*\) é convexa. e) Mostre que a segunda transformada \(g^{**}\) de Legendre de \(g\) está bem definida, onde

\[ g^{**}(r) = \max_{s\in \mathbb{R}} \{rs-g^*(s)\}, \]

f) Mostre que se \(g\) é estritamente convexa, então a segunda transformada de Legendre de \(g\) coincide com a função original, ou seja, \(g^{**} =g\). g) Suponha que \(g\) seja continuamente diferenciável e com \(g': \mathbb{R}\rightarrow \mathbb{R}\) invertível. Mostre que \(g^*\) pode ser escrito na forma

\[ g^*(s) = sr(s) - g(r(s)), \]

para todo \(s\in \mathbb{R}\), onde \(r(s)\) é o único número real \(r\) tal que \(s=g'(r)\). De outra forma,

\[ g^*(s) = s(g')^{-1}(s) - g((g')^{-1}(s)). \]

h) Seja \(g:\mathbb{R}\rightarrow \mathbb{R}\) duas vezes continuamente diferenciável e tal que a sua derivada \(g':\mathbb{R}\rightarrow \mathbb{R}\) é bijetiva. Para \(s\in \mathbb{R}\), considere \(h:\mathbb{R}\rightarrow\mathbb{R}\) dada por

\[ h(s) = -g(0) + \int_{g'(0)}^s (g')^{-1}(\xi)\;\mathrm{d}\xi. \]

Usando a substituição de variáveis \(\xi=g'(\eta)\) e integração por partes mostre que \(h(s)=rs-g(r)\) onde \(r\) é tal que \(s=g'(r)\), i.e. \(h\) é a transformada de Legendre de \(g\).