Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике Рассмотрим систему $m$ дифференциальных уравнений первого порядка, которую запишем в векторной форме где под $x$ и $f(x)$ подразумеваются векторы-столбцы с элементами $x_{k}$ и $f_{k}(x)(k=1, \ldots, m)$. Пусть $x=x^{(0)}$ есть равновесное решение, и предположим, что $f_{k}$ являются не зависящими от $t$ регулярными функциями величин $x_{1}, \ldots, x_{m}$ в некоторой окрестности $x^{(0)}$. Без ограничения общности можно принять, что $x^{(0)}=0$, так что ряд Тейлора для функций $f_{k}(x)$ по переменным $x_{l}$ вблизи $x=0$ запишется в виде В векторных обозначениях это равенство можно записать короче: Сначала в правой части пренебрежем членами высших порядков и вместо системы уравнений (1) будем решать линейную систему Попытаемся использовать решения системы (2) для решения предыдущей нелинейной системы (1). Это не всегда удается, однако мы докажем, что если (1) есть система Гамильтона, то из периодического решения системы (2), вообще говоря, можно получить также периодическое решение системы (1). Для решения (2) используем известную алгебраическую теорему о нормальной форме квадратной матрицы. Собственными значениями $\mathfrak{A}$ будут $m$ корней $\lambda_{1}, \ldots, \lambda_{m}$ характеристического уравнения $m$-й степени причем $\mathfrak{E}$ есть единичная матрица порядка $m$. Допустим далее, что эти $m$ корней различны. Тогда по упомянутой теореме существует такая обратимая комплексная матрица $\mathfrak{C}$, что матрица будет диагональной с диагональными элементами $\lambda_{1}, \ldots, \lambda_{m}$. Рассмотрим также, в какой степени матрица $\mathfrak{C}$ определяется матрицей $\mathfrak{A}$. Если обозначить через $c^{(1)}, \ldots, c^{(m)}$ один из столбцов матрицы $\mathfrak{C}$, то матричное уравнение $\mathfrak{A} \mathfrak{C}=\mathfrak{C} \mathfrak{L}$ равносильно $m$ векторным уравнениям откуда получим Так как все характеристические корни $\lambda_{k}$ различны, то матрица имеет ранг $m-1$, и потому $c^{(k)}$ определяются из уравнений (5) с точностью до постоянного множителя $p_{k}$. С другой стороны, если положить с произвольными $p_{k} Если матрица $\mathfrak{A}$ действительна, то характеристическое уравнение (3) имеет действительные коэффициенты, и поэтому для каждого корня $\lambda_{k}$ найдется комплексно сопряженная величина $\bar{\lambda}_{k}$, которая также является корнем $\lambda_{l}$ уравнения (3). При этом каждому $k=1, \ldots, m$ будет соответствовать точно одно $l=l_{k}$ из ряда $1, \ldots, m$, так что $\bar{\lambda}_{k}=$ $=\lambda_{l}$; в частности, $k=l$ для действительного $\lambda_{k}$. Так как все $\lambda_{k}$ попарно различны, то соответствие между $k$ и $l_{k}$ взаимно-однозначно. Теперь в силу уравнений (5) также таким образом, со скалярными $\rho_{k}$. Так как имеем также $c^{(k)}=\overline{c^{(l)}} \rho_{l}$, и вследствие $|\mathfrak{C}| Так как вместо $c^{(k)}$ можно написать $c^{(k)} \rho_{l}^{1 / 2}$, то $\rho_{k}=1$, значит, $c^{(l)}=\overline{c^{(k)}}$. Это же можно увидеть непосредственно из уравнений (5). Но для дальнейшего применения к системам Гамильтона удобнее не пользоваться нормированием $\rho_{k}=1$. систему (2) можно преобразовать в систему Если $y_{1}, \ldots, y_{m}$ – элементы столбца $y$, то полное решение (8) запишется в виде с $m$ постоянными интегрирования $\alpha_{1}, \ldots, \alpha_{m}$. В частности, $y_{k}$ будет периодической функцией действительного переменного $t$, если $\lambda_{k}$ чисто мнимое. Для исследования условий вещественности используем равенство (6). Чтобы было действительным, должно быть также откуда в силу $|\mathfrak{C}| и произведем замену переменных аналитическую в окрестности $y_{1}-0, \ldots, y_{u}-0$ и переводящую $y-0$ в $x=0$. Пусть соответствующее разложение Тейлора будет и предположим, что $|\mathfrak{B}| Это доказывает, что собственные значения $\lambda_{1}, \ldots, \lambda_{m}$ остаются инвариантными при аналитических преобразованиях дифференциальных уравнений (1). Остановимся теперь на дифференциальных уравнениях (1), имеющих форму Гамильтона Чтобы эту систему также записать в векторной форме, обозначим через $w$ столбец из $2 n$ элементов $u_{1}, \ldots, u_{n}, v_{1}, \ldots, v_{n}$ и через $H_{w}-$ столбец соответствующих производных от $H$. Если обозначить опять через $\mathfrak{E}$ единичную матрицу порядка $n$ и положить то систему уравнений (11) можно записать в сокращенном виде Пусть функция Гамильтона $H=H(w)$ будет регулярной в окрестности $w=0$. Так как постоянный член ряда Тейлора для функции $H$ при $w=0$ для дифференциальных уравнений (11) не имеет значения, то его можно положить равным нулю. Если $w=0$ соответствует состоянию равновесия системы (11), то должны быть равными нулю и первые производные от $H$ при $w=0$. Таким образом, степенной ряд для $H$ начинается с членов второго порядка, и потому можно положить где через $\mathfrak{S}$ обозначена симметричная матрица порядка $2 n$ и через $w^{\prime}$ транспозиция соответствующей строки $w$. Тогда и уравнение (12) переходит в Поэтому нужно исследовать характеристический многочлен $p(\lambda)=$ $=|\lambda \mathfrak{E}-\mathfrak{J} \mathfrak{S}|$. Теперь $\mathfrak{J}^{\prime}=\mathfrak{J}^{-1}=-\mathfrak{J},|\mathfrak{J}|=1, \mathfrak{S}^{\prime}=\mathfrak{S}$, значит, таким образом, $p(\lambda)$ является четной функцией. Если $\lambda$ есть корень многочлена $p(\lambda)$, то корнем будет и $-\lambda$, и притом оба корня будут иметь одинаковую кратность. Если нуль также корень, то его кратность будет четной. Допустим опять, что все собственные значения являются простыми и, следовательно, не равны нулю; при соответствующем расположении их можно обозначить через $\lambda_{k}, \lambda_{k+n}=-\lambda_{k}(k=1, \ldots, n)$. Положим $\mathfrak{L}_{0}=\left[\lambda_{1}, \ldots, \lambda_{n}\right]$ и найдем тогда такую обратимую комплексную матрицу $\mathfrak{C}$ порядка $2 n$, что будет нормальной формой для $\mathfrak{J} \mathfrak{S}$. Переходом к транспонированной матрице получаем С другой стороны, матрица будет симметричной, так что Из уравнений (15) и (16) следует Если положить, кроме того, $\mathfrak{B}=\left(\mathfrak{J}^{-1} \mathfrak{C}^{\prime} \mathfrak{J}\right)^{-1}$, то $|\mathfrak{B}| Но из нашего прежнего определения $\mathfrak{C}$ следует $\mathfrak{C}=\mathfrak{B P}$ с обратимой диагональной матрицей $\mathfrak{P}$, которую мы представим в виде с двумя диагональными матрицами $\mathfrak{P}_{1}, \mathfrak{P}_{2}$ порядка $n$. Следовательно, Кроме того, матрица $\mathfrak{C}^{\prime} \mathfrak{J C}$ альтернированная, так как альтернированной будет $\mathfrak{J}$;ледовательно, $\mathfrak{P}_{1}=\mathfrak{P}_{2}$. Полагая тогда получим $\mathfrak{Q}^{\prime} \mathfrak{J Q}=\mathfrak{J P}$. Если обозначить, наконец, $\mathfrak{C Q}^{-1}$ опять через $\mathfrak{C}$, то уравнение (14) останется справедливым, и, кроме того, теперь следовательно, матрица $\mathfrak{C}$ является симплектической. Отсюда далее следует Будем понимать под $z$ столбец с $2 n$ элементами $x_{1}, \ldots, x_{n}$, $y_{1}, \ldots, y_{n}$ и рассмотрим линейную подстановку Так как функциональная матрица $\mathfrak{C}$ симплектическая, то подстановка (19) в соответствии с уравнением $(2 ; 20)$ будет каноническим преобразованием. Оно переводит систему Гамильтона (12) в где в силу уравнений (13) и (18) поэтому квадратичные члены $H$ будут иметь нормальную форму. на которую она справа может быть умножена. Соответствующим выбором матрицы $\mathfrak{R}$ можно упростить вид условий, при которых $z$ действительно. Прежде всего пусть собственное значение $\lambda_{k}(k \leqslant n)$ не является чисто мнимым, и $\lambda_{l}=\bar{\lambda}_{k}\left(l=l_{k}\right.$ ). Если при этом $l>n$, то $\lambda_{k} Если $\lambda_{k}$ не чисто мнимое, то уравнение (6) дает формулу но если $\lambda_{k}$ чисто мнимое, то $l_{k}=k+n$, и Если теперь опять заменить произведение матриц $\mathfrak{C}$ матрицей $\mathfrak{C}$, то первоначальные столбцы $c^{(k)}, c^{(k+n)}$ умножатся на скалярные множители $r_{k}, r_{k}^{-1}$. Если собственное значение $\lambda_{k}$ не чисто мнимое, то тогда в соответствии с уравнением (6) $\rho_{k}$ умножаются на множитель $r_{l} \bar{r}_{k}^{-1}$. В случае когда $\lambda_{k}$ не будет действительным и $k<l_{k}$, выбирая $r_{k}=\bar{\rho}_{k}$, $r_{l}=1$, получим $\rho_{k}=1$; тогда в силу уравнений $(7),(22)$ имеем также $\rho_{k+n}=1, \rho_{l}=1, \rho_{l+n}=1$. Если $\lambda_{k}$ действительно, значит, $k=l_{k}$, и, согласно уравнению (7), модуль числа $\rho_{k}$ равен 1 ; заменяя тогда $r_{k}$ на $\rho_{k}^{-1 / 2}$, получим сразу искомый результат. Наконец, если $\lambda_{k}$ чисто мнимое, то $\rho_{k}$ имеет положительный множитель $\left(r_{k} \bar{r}_{k}\right)^{-1}$, в то время как в соответствии с уравнениями (7), (23) из уравнения $\bar{\rho}_{k}=-\rho_{k}$ следует, что $\rho_{k}$ будет чисто мнимым. Поэтому в этом случае можно добиться выполнения равенства $\rho_{k}=\rho_{k+n}= \pm i$. Если переставить $\lambda_{k}$ и $\lambda_{k+n}$, то можно положить $c^{(k)}=c^{(k+n)}$ и $c^{(k+n)}=-c^{(k)}$, причем $\mathfrak{C}$ остается симплектической. Так как тогда $\rho_{k}$ заменяется на $-\rho_{k+n}$, то можно, следовательно, сделать $\rho_{k}=-i$. При этом $\rho_{k}=-i$ будет множителем для всех чисто мнимых собственных значений $\lambda_{k}$, в противном случае $\rho_{k}=1$. Условие вещественности $\mathfrak{C} z=\overline{\mathfrak{C}} \bar{z}$ совпадает с $\bar{z}_{k}=\rho_{k} z_{l_{k}}$, следовательно, $x_{l}=\bar{x}_{k}, y_{l}=\bar{y}_{k}$ для $\lambda_{l}=\bar{\lambda}_{k} Если рассматривать в правой части уравнения (20) только линейный член, то получится линейная система $\dot{z}=\mathfrak{L} z$, или, точнее, с общим решением где $\xi_{k}, \eta_{k}$ – постоянные интегрирования. Если одно из собственных значений $\lambda_{k}$ чисто мнимое, например $\lambda_{1}$, то решение есть периодическое решение линейной системы. При этом, чтобы $w$ было действительным, нужно взять $\eta_{1}=i \bar{\xi}_{1}$. Цель ближайших двух параграфов состоит в том, чтобы, исходя из этого решения, получить периодическое решение общей системы нелинейных уравнений Гамильтона (11). Что же касается не гамильтоновых дифференциальных уравнений (1), то может случиться, что они, кроме постоянных решений, не имеют никаких других периодических решений, в то время как соответствующая линейная система (2) может иметь периодическое решение, не являющееся константой. Это можно показать на примере системы с двумя неизвестными функциями $x, y$, где $g$ есть заданное натуральное число. Соответствующая линейная система $\dot{x}=-y, \dot{y}=x$ имеет только периодические решения, а именно $x=\alpha \cos t+\beta \sin t, y=\alpha \sin t-$ $-\beta \cos t$ с постоянными $\alpha, \beta$. Это будут концентрические окружности в плоскости $(x, y)$, которые пробегаются с постоянной угловой скоростью, равной единице, и центр этих окружностей находится в начале координат. С другой стороны, если записать данную систему (24) в полярных координатах с помощью подстановки $x=r \cos \varphi, y=r \sin \varphi$, то Для решений, отличных от тривиального решения $r=0$, получим с двумя постоянными интегрирования $a, b$. Эти решения суть спирали, следовательно они заведомо не будут периодическими. Данный пример, в частности, показывает, что свойство системы $\dot{x}=-y, \dot{y}=x$ иметь периодические решения исчезнет, если добавить в правые части члены более высоких степеней относительно $x$ и $y$; например, в нашем случае можно выбрать параметр $g=1,2, \ldots$ Однако заметим, что данная система не является канонической. Теперь нам нужно исследовать, как можно получить периодические решения системы Гамильтона, если их имеет соответствующая линейная система и если, кроме того, выполняется некоторое условие.
|
1 |
Оглавление
|