Гамильтоновы системы: различия между версиями
Timur23 (обсуждение | вклад) |
Timur23 (обсуждение | вклад) |
||
Строка 4: | Строка 4: | ||
Функция $$V(x_1, x_2,\dots, x_n)$$ называется '''первым интегралом''', если $$V(x_1, x_2,\dots, x_n) = const$$ или если производная в силу системы функции $$V(x_1, x_2,\dots, x_n)$$ равна нулю. | Функция $$V(x_1, x_2,\dots, x_n)$$ называется '''первым интегралом''', если $$V(x_1, x_2,\dots, x_n) = const$$ или если производная в силу системы функции $$V(x_1, x_2,\dots, x_n)$$ равна нулю. | ||
− | Рассмотрим | + | Рассмотрим систему уравнений |
\begin{gather*} | \begin{gather*} | ||
− | m \ddot{x} = F(x),\quad x(0) = x_0,\quad \dot{x}(0) = p_0, \quad F(x) = \big(f_1(x), f_2(x), | + | \begin{cases} |
+ | \dot{x}_1 = f_1(x_1, x_2),\\ | ||
+ | \dot{x}_2 = f_2(x_1, x_2). | ||
+ | \end{cases} | ||
+ | \end{gather*} | ||
+ | Уравнение для первого интеграла этой системы будет иметь вид: | ||
+ | \begin{gather*} | ||
+ | \frac{dx_1}{dx_2} = \frac{f_1(x_1, x_2)}{f_2(x_1, x_2)} | ||
+ | \end{gather*} | ||
+ | или | ||
+ | \begin{gather*} | ||
+ | \frac{dx_1}{f_1(x_1, x_2)} = \frac{dx_2}{f_2(x_1, x_2)}. | ||
+ | \end{gather*} | ||
+ | |||
+ | Если система зависит от $$n$$ переменных, то уравнение для первого интеграла будет иметь вид: | ||
+ | \begin{gather*} | ||
+ | \frac{dx_1}{f_1(x_1, x_2)} = \frac{dx_2}{f_2(x_1, x_2)} = \dots = \frac{dx_n}{f_n(x_1, x_2)}. | ||
+ | \end{gather*} | ||
+ | В этом случае можно получить $$(n-1)$$ первый интеграл. | ||
+ | |||
+ | '''Пример''' | ||
+ | |||
+ | Рассмотрим систему | ||
+ | \begin{gather*} | ||
+ | \begin{cases} | ||
+ | \dot{x}_1 = - x_2,\\ | ||
+ | \dot{x}_2 = x_1. | ||
+ | \end{cases} | ||
+ | \end{gather*} | ||
+ | |||
+ | Найдем первый интеграл этой системы: | ||
+ | \begin{gather*} | ||
+ | - \frac{dx_1}{x_2} = \frac{dx_2}{x_1} | ||
+ | \end{gather*} | ||
+ | \begin{gather*} | ||
+ | - x_1 dx_1 = x_2 dx_2 | ||
+ | \end{gather*} | ||
+ | \begin{gather*} | ||
+ | - x_1^2 = x_2^2 = const | ||
+ | \end{gather*} | ||
+ | Откуда | ||
+ | \begin{gather*} | ||
+ | V(x_1, x_2) = x_1^2 + x_2^2 | ||
+ | \end{gather*} | ||
+ | Найдем производную в силу системы функции $$V(x_1, x_2)$$: | ||
+ | \begin{gather*} | ||
+ | \dot{V}(x_1, x_2) = 2 x_1 \dot{x}_1 + 2 x_2 \dot{x}_2 = - 2 x_1 x_2 + 2 x_2 x_1 \equiv 0. | ||
+ | \end{gather*} | ||
+ | Убедились в том, что найденная функция $$V(x_1, x_2)$$ является первым интегралом рассматриваемой системы. | ||
+ | |||
+ | Рассмотрим уравнение Ньютона (x $$\in \mathbb{R}^n$$): | ||
+ | \begin{gather*} | ||
+ | m \ddot{x} = F(x),\quad x(0) = x_0,\quad \dot{x}(0) = p_0, \quad F(x) = \big(f_1(x), f_2(x), \dots, f_n(x)\big). | ||
\end{gather*} | \end{gather*} | ||
Это уравнение эквивалентно системе | Это уравнение эквивалентно системе | ||
Строка 15: | Строка 67: | ||
\dot{p}_i = f_i(x), \quad f_i(x) = \frac{1}{m}F_i(x). | \dot{p}_i = f_i(x), \quad f_i(x) = \frac{1}{m}F_i(x). | ||
\end{cases}, \quad | \end{cases}, \quad | ||
− | i = \overline{1, | + | i = \overline{1, n} |
\end{equation} | \end{equation} | ||
'''Определение 2.''' | '''Определение 2.''' | ||
Система '''потенциальна''', если существует функция $$U(x)$$, называемая '''потенциалом системы''' такая, что | Система '''потенциальна''', если существует функция $$U(x)$$, называемая '''потенциалом системы''' такая, что | ||
\begin{gather*} | \begin{gather*} | ||
− | \frac{dU(x)}{dx_i}=-f_i(x), \quad i = \overline{1, | + | \frac{dU(x)}{dx_i}=-f_i(x), \quad i = \overline{1, n}. |
\end{gather*} | \end{gather*} | ||
− | Потенциал определяется с точностью до произвольной постоянной. Далее предположим, что $$f_i(x)$$ — непрерывная функция, $$i = \overline{1, | + | Потенциал определяется с точностью до произвольной постоянной. Далее предположим, что $$f_i(x)$$ — непрерывная функция, $$i = \overline{1, n}$$. |
'''Определение 3.''' | '''Определение 3.''' | ||
Строка 29: | Строка 81: | ||
\begin{equation} | \begin{equation} | ||
\label{eq2} | \label{eq2} | ||
− | H(x, p) = \sum_{i=1}^{ | + | 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} = \sum_{i=1}^{ | + | \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}^{ | + | \sum_{i=1}^{n} p_i f_i(x) + \sum_{i=1}^{n} -f_i(x) p_i \equiv 0 |
\end{gather*} | \end{gather*} | ||
Строка 41: | Строка 93: | ||
Для того, чтобы отыскать значение постоянной $$C$$, которая отвечает движению с начальными данными $$x(0) = x_0, p(0) = p_0$$ достаточно вычислить | Для того, чтобы отыскать значение постоянной $$C$$, которая отвечает движению с начальными данными $$x(0) = x_0, p(0) = p_0$$ достаточно вычислить | ||
\begin{gather*} | \begin{gather*} | ||
− | H(x_0, p_0) = \sum_{i=1}^{ | + | H(x_0, p_0) = \sum_{i=1}^{n} \frac{p_{0, i}^2}{2} + U(x_0). |
\end{gather*} | \end{gather*} | ||
Версия 00:27, 28 ноября 2023
Определения
Определение 1. Функция $$V(x_1, x_2,\dots, x_n)$$ называется первым интегралом, если $$V(x_1, x_2,\dots, x_n) = const$$ или если производная в силу системы функции $$V(x_1, x_2,\dots, x_n)$$ равна нулю.
Рассмотрим систему уравнений \begin{gather*} \begin{cases} \dot{x}_1 = f_1(x_1, x_2),\\ \dot{x}_2 = f_2(x_1, x_2). \end{cases} \end{gather*} Уравнение для первого интеграла этой системы будет иметь вид: \begin{gather*} \frac{dx_1}{dx_2} = \frac{f_1(x_1, x_2)}{f_2(x_1, x_2)} \end{gather*} или \begin{gather*} \frac{dx_1}{f_1(x_1, x_2)} = \frac{dx_2}{f_2(x_1, x_2)}. \end{gather*}
Если система зависит от $$n$$ переменных, то уравнение для первого интеграла будет иметь вид: \begin{gather*} \frac{dx_1}{f_1(x_1, x_2)} = \frac{dx_2}{f_2(x_1, x_2)} = \dots = \frac{dx_n}{f_n(x_1, x_2)}. \end{gather*} В этом случае можно получить $$(n-1)$$ первый интеграл.
Пример
Рассмотрим систему \begin{gather*} \begin{cases} \dot{x}_1 = - x_2,\\ \dot{x}_2 = x_1. \end{cases} \end{gather*}
Найдем первый интеграл этой системы: \begin{gather*} - \frac{dx_1}{x_2} = \frac{dx_2}{x_1} \end{gather*} \begin{gather*} - x_1 dx_1 = x_2 dx_2 \end{gather*} \begin{gather*} - x_1^2 = x_2^2 = const \end{gather*} Откуда \begin{gather*} V(x_1, x_2) = x_1^2 + x_2^2 \end{gather*} Найдем производную в силу системы функции $$V(x_1, x_2)$$: \begin{gather*} \dot{V}(x_1, x_2) = 2 x_1 \dot{x}_1 + 2 x_2 \dot{x}_2 = - 2 x_1 x_2 + 2 x_2 x_1 \equiv 0. \end{gather*} Убедились в том, что найденная функция $$V(x_1, x_2)$$ является первым интегралом рассматриваемой системы.
Рассмотрим уравнение Ньютона (x $$\in \mathbb{R}^n$$): \begin{gather*} m \ddot{x} = F(x),\quad x(0) = x_0,\quad \dot{x}(0) = p_0, \quad F(x) = \big(f_1(x), f_2(x), \dots, f_n(x)\big). \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} Определение 2. Система потенциальна, если существует функция $$U(x)$$, называемая потенциалом системы такая, что \begin{gather*} \frac{dU(x)}{dx_i}=-f_i(x), \quad i = \overline{1, n}. \end{gather*}
Потенциал определяется с точностью до произвольной постоянной. Далее предположим, что $$f_i(x)$$ — непрерывная функция, $$i = \overline{1, n}$$.
Определение 3. Первый интеграл системы (\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) = \sum_{i=1}^{n} \frac{p_{0, i}^2}{2} + U(x_0). \end{gather*}
Определение 4. Пусть $$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}
Примеры
Пример 1
Рассмотрим следующий пример гамильтоновой системы: \begin{gather*} \frac{dx}{dt} = p, \\ \frac{dp}{dt} = f(x) \end{gather*}
Список литературы
- Абрамова В.В. "Лекции по динамическим системам и биоматематике", 2023.
- Братусь А.С., Новожилов А.С., Платонов А.П. "Динамические системы и модели биологии", 2011.