3. Если уравнение (8.2) не имеет корней с положительной вещественной частью, но имеется часть корней с нулевой вещественной частью, то может иметь место как устойчивость (не асимптотическая), так и неустойчивость.

Таким образом, вопрос об устойчивости нулевого решения системы (8.1) сводится к исследованию характера корней уравнения (8.2). Раскрывая определитель (8.2), получим уравнение

/(Х) = Х« + а,Х«-+...+ а„ = 0. (8.3)

Из коэффициентов полинома f(k) составим матрицу 0 0 0 0. 10 0. 1 .

Дв ui ui ai

В записи этой матрицы полагаем а„ = 0, если тп. Рассмотрим определители

а ац

Для того чтобы все корни уравнения (8.3) имели отрицательные вещественные части, необходимо и достаточно выполнение неравенств

\>0 при k==\, 2, ..., я.

(8.4)

Это утверждение носит название теоремы Гурвица; доказательство теоремы Гурвица имеется, например, в монографии [26]. Условия (8.4) часто называют условиями Рауза- Гурвица.

Для уравнения второй степени

X«-f-aiX + a2 = 0

условия Рауза-Гурвица имеют вид ajO, или, что то же самое, а, 0, 0.

а, О

линейные системы

Для уравнения третьей степени условия Рауза-Гурвица запишутся в виде

а,>0,

>0, а,

а, 1 аз а.з

очевидно, эти условия эквивалентны условиям

а, > О, а > О, aia.2 > aj.

В дальнейшем нам понадобятся некоторые простейшие сведения из линейной алгебры.

Рассмотрим квадратичную форму t)== biXiXji, где

i, ft = 1

й;;, = bi. Обозначая через В матрицу, составленную из коэффициентов bi, и через х вектор с проекциями х,, ... , х„, можем записать

v = {x, Вх),

где {х, Вх) означает скалярное произведение векторов л: и Вх.

Если вектор х является функцией времени t, то и форма v также будет функцией t и имеет место следующее правило дифференцирования:

(8.5)

Если матрица В не симметрична, то для векторов хну имеет место соотношение

{X, Ву) = {В*х, у).

(8.6)

где В* - транспонированная матрица.

Заметим еще следующий хорошо известный факт: если матрица Т-неособая матрица, т. е. матрица, определитель которой отличен от нуля, то матрицы В и ТВТ (подобные матрицы) имеют одинаковые собственные числа.

Запишем систему (8.1) в матричной форме

dt -

(8.7)

2 Е. А. Барбашин

и найдем производную квадратичной формы v = (x, Вх) в силу этой системы. Имеем согласно (8.5)

=(Ах, Вх)-\- (х, В Ах).

Используя соотношение (8.6), получим

=(х, [А*В + ВА]х). (8.8)

Потребуем теперь, чтобы форма v удовлетворяла уравнению

S = (8.9)

гае w = {x, Сх) - заданная квадратичная форма. Сравнивая (8.8) и (8.9), получим матричное уравнение для определения матрицы В

А*В-{-ВА = С. (8.10)

Уравнение (8.10) позволяет найти матрицу формы v по заданной матрице формы w. Исследование уравнения (8.10) представляет большой интерес, так как это уравнение позволяет найти функцию Ляпунова в виде квадратичной формы по заданной производной. Уравнение (8.10) ставит в соответствие всякой матрице В матрицу С, причем это соответствие линейно. Таким образом, в пространстве квадратичных матриц л-го порядка может быть определен линейный оператор F{B) = A*B-\-BA. Задача о разрешимости уравнения (8.10) сводится к задаче определения оператора, обратного оператору F, так как B - F~C. Так как оператор F действует в конечномерном пространстве (я измерений), то для существования обратного оператора необходимо и достаточно, чтобы среди собственных чисел оператора F не было нулевого.

Теорема 8.1. Если корни характеристического уравнения системы (8.7) таковы, что сумма i-j-i, не обращается в нуль ни при каких I, k, то, какова бы ни была наперед заданная квадратичная форма w, существует единственная квадратичная форма v, удовлетворяющая уравнению (8.9).

Докажем теорему. По определению, собственным числом оператора F называется число [х, такое, что уравнение F(B) =