title: Equazione di Binet
text: # Equazione di Binet
  
  L'equazione di Binet, dovuta a Jacques Philippe Marie Binet, fornisce la forma di una forza centrale data la traiettoria orbitale in coordinate polari. L'equazione può anche essere usata per ricavare la forma di un'orbita data una certa forza, ma ciò solitamente comporta la risoluzione di un'equazione differenziale ordinaria non-lineare del secondo ordine. In caso di moto circolare intorno al centro della forza una soluzione unica è impossibile.
  
  ## Equazione
  
  La forma di un'orbita è spesso facilmente descritta in termini della distanza relativa $r$ come funzione dell'angolo $\theta$. Per l'equazione di Binet, la forma orbitale è invece descritta dal suo reciproco $u=1/r$ in funzione di $\theta$. Si definisca il momento angolare specifico $h=L/m$ dove $L$ è il momento angolare e $m$ la massa. L'equazione di Binet, ricavata nella prossima sezione, dà la forza in termini della funzione $u(\theta )$:
  
      $F({u}^{-1})=-mh^{2}u^{2}\left({\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u\right).$
  
  ## Dimostrazione
  
  La seconda legge di Newton per una forza puramente centrale è
  
      $F(r)=m({\ddot {r}}-r{\dot {\theta }}^{2}).$
  
  Per la conservazione del momento angolare segue che
  
      $r^{2}{\dot {\theta }}=h={\text{costante}}.$
  
  Le derivate di $r$ rispetto al tempo possono essere riscritte come derivate di $u=1/r$ rispetto all'angolo:
  
      ${\begin{aligned}&{\frac {\mathrm {d} u}{\mathrm {d} \theta }}={\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {1}{r}}\right){\frac {\mathrm {d} t}{\mathrm {d} \theta }}=-{\frac {\dot {r}}{r^{2}{\dot {\theta }}}}=-{\frac {\dot {r}}{h}}\\&{\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}=-{\frac {1}{h}}{\frac {\mathrm {d} {\dot {r}}}{\mathrm {d} t}}{\frac {\mathrm {d} t}{\mathrm {d} \theta }}=-{\frac {\ddot {r}}{h{\dot {\theta }}}}=-{\frac {\ddot {r}}{h^{2}u^{2}}}\\\end{aligned}}$
  
  Mettendo insieme i risultati precedenti, si arriva a
  
      $F=m({\ddot {r}}-r{\dot {\theta }}^{2})=-m\left(h^{2}u^{2}{\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+h^{2}u^{3}\right)=-mh^{2}u^{2}\left({\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u\right)$
  
  ## Esempi
  
  ### Problema di Keplero
  
  Per calcolare l'orbita data da una forza che soddisfa la legge dell'inverso del quadrato si risolve l'equazione di Binet ottenendo:
  
      $-ku^{2}=-mh^{2}u^{2}\left({\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u\right)$
      ${\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u={\frac {k}{mh^{2}}}\equiv {\text{costante}}>0.$
  
  Se l'angolo $\theta$ è misurato dal periapside, allora la soluzione generale per l'orbita espressa in coordinate polari (reciproche) è
  
      $lu=1+\varepsilon \cos \theta .$
  
  L'equazione di cui sopra descrive sezioni coniche, con $l$ il semilato retto (pari a $h^{2}/\mu =h^{2}m/k$) e $\varepsilon$ l'eccentricità orbitale.
  
  L'equazione relativistica ricavata per le coordinate di Schwarzschild è
  
      ${\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u={\frac {r_{s}c^{2}}{2h^{2}}}+{\frac {3r_{s}}{2}}u^{2}$
  
  dove $c$ è la velocità della luce e $r_{s}$ è il raggio di Schwarzschild. Per la metrica di Reissner-Nordström si otterrà
  
      ${\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u={\frac {r_{s}c^{2}}{2h^{2}}}+{\frac {3r_{s}}{2}}u^{2}-{\frac {GQ^{2}}{4\pi \varepsilon _{0}c^{4}}}\left({\frac {c^{2}}{h^{2}}}u+2u^{3}\right)$
  
  dove $Q$ è la carica elettrica e $\varepsilon _{0}$ è la costante dielettrica del vuoto.
  
  ### Problema di Keplero inverso
  
  Che tipo di forza produce un'orbita ellittica (o più in generale una sezione conica) intorno a un fuoco dell'ellisse?
  
  Facendo due volte la derivata dell'equazione di Binet per un'ellisse si ha
  
      $l\,{\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}=-\varepsilon \cos \theta .$
  
  La legge della forza è perciò
  
      $F=-mh^{2}u^{2}\left({\frac {-\varepsilon \cos \theta }{l}}+{\frac {1+\varepsilon \cos \theta }{l}}\right)=-{\frac {mh^{2}u^{2}}{l}}=-{\frac {mh^{2}}{lr^{2}}},$
  
  che è la legge dell'inverso del quadrato come previsto. Facendo combaciare i valori orbitali $h^{2}/l$ a valori fisici come $GM$ o $k_{e}q_{1}q_{2}/m$ si riproduce rispettivamente la legge di gravitazione universale o la legge di Coulomb.
  
  