Гамильтоновы системы: различия между версиями
Timur23 (обсуждение | вклад) |
Timur23 (обсуждение | вклад) |
||
Строка 14: | Строка 14: | ||
\end{equation} | \end{equation} | ||
'''Определение 1.''' | '''Определение 1.''' | ||
− | Функция $$U(x)$$ называется '''потенциалом системы''', если $$ | + | Функция $$U(x)$$ называется '''потенциалом системы''', если $$\frac{dU(x)}{dx_i}=-f_i(x)$$. |
− | Далее предположим, что $$ | + | Далее предположим, что $$f_i(x)$$ — непрерывная функция, $$i = \overline{1, n}$$. Потенциал определяется с точностью до произвольной постоянной. |
'''Определение 2.''' | '''Определение 2.''' | ||
Строка 22: | Строка 22: | ||
\begin{equation} | \begin{equation} | ||
\label{eq2} | \label{eq2} | ||
− | H(x, p) = \frac{ | + | H(x, p) = \sum_{i=1}^{n} \frac{p_i^2}{2} + U(x) = const = C. |
\end{equation} | \end{equation} | ||
Действительно, непосредственно проверяется, что | Действительно, непосредственно проверяется, что | ||
\begin{gather*} | \begin{gather*} | ||
− | \frac{dH(x, p)}{dt} = \frac{\partial H}{\partial | + | \frac{dH(x, p)}{dt} = \sum_{i=1}^{n} \frac{\partial H}{\partial p_i}\dot{p}_i + \sum_{i=1}^{n} \frac{\partial H}{\partial x_i}\dot{x}_i = |
+ | \sum_{i=1}^{n} p_i f_i(x) + \sum_{i=1}^{n} -f_i(x) p_i \equiv 0 | ||
\end{gather*} | \end{gather*} | ||
Версия 23:33, 11 октября 2023
Определения
Рассмотрим многомерное уравнение Ньютона \begin{gather*} m \ddot{x} = F(x),\quad x(0) = x_0,\quad \dot{x}(0) = p_0,\quad x, p \in \mathbb{R}^n \end{gather*} Это уравнение эквивалентно системе \begin{equation} \label{eq1} \begin{cases} \dot{x}_i = p_i,\\ \dot{p}_i = f_i(x), \quad f_i(x) = \frac{1}{m}F_i(x). \end{cases}, \quad i = \overline{1, n} \end{equation} Определение 1. Функция $$U(x)$$ называется потенциалом системы, если $$\frac{dU(x)}{dx_i}=-f_i(x)$$.
Далее предположим, что $$f_i(x)$$ — непрерывная функция, $$i = \overline{1, n}$$. Потенциал определяется с точностью до произвольной постоянной.
Определение 2. Первый интеграл системы (\ref{eq1}) называется гамильтонианом и определяется равенством \begin{equation} \label{eq2} H(x, p) = \sum_{i=1}^{n} \frac{p_i^2}{2} + U(x) = const = C. \end{equation}
Действительно, непосредственно проверяется, что \begin{gather*} \frac{dH(x, p)}{dt} = \sum_{i=1}^{n} \frac{\partial H}{\partial p_i}\dot{p}_i + \sum_{i=1}^{n} \frac{\partial H}{\partial x_i}\dot{x}_i = \sum_{i=1}^{n} p_i f_i(x) + \sum_{i=1}^{n} -f_i(x) p_i \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}).
Определение 3. Пусть $$x, p \in \mathbb{R}^n$$. Гамильтоновой системой называется система вида \begin{equation} \label{eq3} \begin{cases} \frac{dx_i}{dt} = \frac{\partial H}{\partial p_i}\\ \frac{dp_i}{dt} = -\frac{\partial H}{\partial x_i} \end{cases}, \quad i = \overline{1, n} \end{equation}
Список литературы
- Абрамова В.В. "Лекции по динамическим системам и биоматематике", 2023.
- Братусь А.С., Новожилов А.С., Платонов А.П. "Динамические системы и модели биологии", 2011.