Гамильтоновы системы: различия между версиями

Материал из sawiki
Перейти к навигации Перейти к поиску
Строка 34: Строка 34:
 
H(x_0, p_0) = \frac{p_0^2}{2} + U(x_0).
 
H(x_0, p_0) = \frac{p_0^2}{2} + U(x_0).
 
\end{gather*}
 
\end{gather*}
 +
 +
Положения равновесия системы представляют пары чисел $$(x^∗, 0)$$, таких, что $$f(x^∗) = −U'(x^∗) = 0$$, т.е. $$x^∗$$ является критической точкой потенциала $$U(x)$$. Если задан потенциал системы, то фазовые траектории системы могут быть получены непосредственно из равенства (\ref{eq2}).

Версия 01:14, 10 октября 2023

Определения

Рассмотрим одномерное уравнение Ньютона \begin{gather*} m \ddot{x} = F(x),\quad x(0) = x_0,\quad \dot{x}(0) = p_0 \end{gather*} Это уравнение эквивалентно системе \begin{equation} \label{eq1} \begin{cases} \dot{x} = p,\\ \dot{p} = f(x), \quad f(x) = \frac{1}{m}F(x). \end{cases} \end{equation} Определение 1. Функция $$U(x)$$ называется потенциалом системы, если $$U'(x)=-f(x)$$.

Далее предположим, что $$f(x)$$ — непрерывная функция. Потенциал определяется с точностью до произвольной постоянной.

Определение 2. Первый интеграл системы (\ref{eq1}) называется гамильтонианом и определяется равенством \begin{equation} \label{eq2} H(x, p) = \frac{p^2}{2} + U(x) = const = C. \end{equation}

Действительно, непосредственно проверяется, что \begin{gather*} \frac{dH(x, p)}{dt} = \frac{\partial H}{\partial x}\dot{x} + \frac{\partial H}{\partial p}\dot{p} = -f(x)p + pf(x) \equiv 0 \end{gather*}

Гамильтониан задает полную энергию системы (\ref{eq1}), которая остается неизменной в течение всего времени эволюции системы. Для того, чтобы отыскать значение постоянной $$C$$, которая отвечает движению с начальными данными $$x(0) = x_0, p(0) = p_0$$ достаточно вычислить \begin{gather*} H(x_0, p_0) = \frac{p_0^2}{2} + U(x_0). \end{gather*}

Положения равновесия системы представляют пары чисел $$(x^∗, 0)$$, таких, что $$f(x^∗) = −U'(x^∗) = 0$$, т.е. $$x^∗$$ является критической точкой потенциала $$U(x)$$. Если задан потенциал системы, то фазовые траектории системы могут быть получены непосредственно из равенства (\ref{eq2}).