Методика решения задач механики - Беляшкин А.Г.
Скачать (прямая ссылка):
уравнений возмущенного движения.
32
ЛИНЕЙНЫЕ ГАМИЛЬТОНОВЫ СИСТЕМЫ
[ГЛ. 2
§ 2. Нормальная форма автономной системы линейных
гамильтоновых уравнений в случае простых
чисто мнимых корней характеристического уравнении
Продолжим исследование линейной системы дифференциальных уравнений (1.1).
Предположим, что соответствующее ей характеристическое уравнение (1.3)
имеет только простые чисто мнимые корни. Обозначим их через = iak, Хп+к =
-iak (к = 1,2,. . . . . . ., п). Знаки вещественных величин ак пока не
фиксируем. Они будут определены ниже. Найдем вещественное линейное
каноническое преобразование (/ = 1, 2, . . ., 2п), приводя-
щее систему (1.1) к нормальной форме. Нормальной формой системы уравнений
(1.1) мы называем такую систему дифференциальных уравнений, которой
соответствует функция Гамильтона, равная алгебраической сумме
гамильтонианов п линейных, не связанных между собой осцилляторов:
П
Н= ?Ч(2/2* +2&*). (2.1)
к=1
Решение задачи о получении нормальной формы линейной системы
(1.1) необходимо при исследовании устойчивости нелинейных уравнений
возмущенного движения, при анализе нелинейных колебаний, при построении
приближенных решений нелинейных гамильтоновых систем, где в качестве
первого приближения берется обычно решение линейной задачи. Поэтому
крайне желательно выбирать такие координаты, в которых решение линейной
задачи записывалось бы наиболее просто. Простейшей вещественной формой
уравнений (1.1) и будет нормальная форма.
Общая , задача об алгебраических свойствах линейных систем гамильтоновых
дифференциальных уравнений исследована достаточно подробно [15, 28, 49,
90, 98, 109, 145,149, 154, 157, 179-183]. Для систем с постоянными
коэффициентами в работах [15, 28, 98] получены конструктивные методы
нормализации. Мы рассмотрим задачу получения нормальной формы иначе, чем
в упомянутых работах [15, 28, 98], и получим алгоритм нормализации,
который будет весьма простым, так как его применение сводится только к
нахождению собственных векторов матрицы IH.
Введем обозначение ут = (уг, уп+1, . . ., у2п). Тогда,
учитывая (2.1), получим, что нормальная форма линейной системы (1.1)
запишется в виде следующей гамильтоновой системы уравнений:
^=1Н*у, (2.2)
где Н* - вещественная диагональная матрица, диагональные
НОРМАЛЬНАЯ ФОРМА АВТОНОМНОЙ СИСТЕМЫ
33
элементы которой определены равенствами h*}с = &*+ s.n+k =
(к = 1, 2, . . п). Переход от переменных х к переменным у
зададим с помощью матрицы А в виде равенства
х = Ау. (2.3)
Из (1.1) и (2.2)-(2.3) получаем, что матрица А должна удовлетворять
следующему матричному уравнению:
AIH* = IHA. (2.4)
Кроме того, для каноничности преобразования (2.3) матрицу А ищем
симплектической [16], т. е. она должна удовлетворять еще одному
матричному уравнению
ATIA = I. (2.5)
Решение матричного уравнения (2.4) не единственно. Чтобы найти
нормализующее преобразование, надо из бесчисленного множества решений
матричного уравнения (2.4) выбрать хотя бы одно вещественное,
удовлетворяющее уравнению (2.5).
Решение А уравнения (2.4) будем искать в виде А = ВС, где матрица С
определена равенством
11 iEn Е
С==||-гЕп Е"
(2.6)
Подставив в уравнение (2.4) вместо А его выражение через В и С, получим
следующее уравнение для нахождения матрицы В:
BD = IHB, (2.7)
где D - диагональная форма матрицы IH. Для ее диагональных элементов
имеют место равенства dkli = - dn+kt n+& = iak (ft = 1,2,..., n). Таким
образом, теперь надо найти матрицу В, приводящую матрицу IH исходной
системы уравнений (1.1) к диагональной форме. Она строится следующим
образом [17]. Ее столбцами должны быть собственные векторы матрицы IH.
Именно, пусть т-й столбец матрицы В будет собственным вектором ет,
соответствующим собственному числу iam, а (п + т)-й столбец будет
вектором en+m, соответствующим собственному числу ^m+n ИТ jyi (jit == 1,
2, . . ., /l).
Собственные векторы определяются с точностью до постоянного множителя.
Примем этот множитель вещественным и одинаковым для векторов ет и еп+т.
Кроме того, соответствующие компоненты этих векторов выберем комплексно
сопряженными. Такой выбор собственных векторов обеспечивает
вещественность матрицы А. Произвольные множители собственных векторов
определяются из условия их нормировки, которое ниже будет получено из
условия (2.5) каноничности преобразования (2.3).
2 А. П. Маркеев
34
ЛИНЕЙНЫЕ ГАМИЛЬТОНОВЫ. СИСТЕМЫ
1ГЛ. 2
Подставив выражение А = ВС в уравнение (2.5), получим
CT'B'IBC = I. (2.8)
Рассмотрим подробно матрицу ВТ1В, которую для краткости обозначим через
F. Элемент fn этой матрицы, как нетрудно проверить, равняется скалярному
произведению векторов ек и 1ег:
/кг = (ек> 1ег)- (2-9)
Но так как для любых двух векторов а и Ь справедливо равенство
(а, 1Ь) = - (1а, Ь), (2.10)
то отсюда следует, что матрица F - кососимметрическая. Рас-
смотрим дальше структуру матрицы F. Докажем, что /&г = 0, если | I - к\ф
п.
Для доказательства рассмотрим очевидное равенство
