Второй метод Ляпунова.

Одним из наиболее эффективных способов исследования устойчивости является второй метод Ляпунова (часто его называют также прямым методом Ляпунова), который не требует нахождения решения дифференциального уравнения.

Метод основан на идее замены анализа решений нелинейных уравнений произвольного порядка на оценку свойств этих решений с помощью скалярного дифференциального неравенства. При этом теряется информация о виде решений, но приобретается простота анализа устойчивости, поскольку исследуется изменение «расстояния» в пространстве состояний от текущей точки системы до начала координат. В качестве такой оценки расстояния можно использовать скалярную функцию переменных состояния, которую обозначим через V(x).

Основная философия прямого метода Ляпунова является математическим расширением фундаментального физического наблюдения: если полная энергия механической (электрической) системы непрерывно рассеивается (убывает), то с течением времени система независимо от того является ли она линейной или нелинейной должна стремиться к состоянию равновесия. Отсюда мы можем судить об устойчивости системы путем исследования изменения единственной скалярной функции.

Пример. Рассмотрим нелинейную автономную управляемую систему масса-пружина-демфер, динамика которой описывается уравнением

, (*)

где - нелинейное рассеивание или демпфирование,

- нелинейная сила пружины.

Полная механическая энергия = кинетическая энергия + потенциальная энергия,

. (**)

Сопоставляя определения устойчивости и механической энергии, мы можем увидеть определенную связь между механической энергией и понятиями, описанными ранее:

· Нулевая энергия соответствует состоянию равновесия ( )

· Асимптотическая устойчивость влечет за собой сходимость механической энергии к нулю

· Неустойчивость связана с ростом механической энергии.

Эти связи говорят о том, значение скалярной величины, механической энергии непосредственно отражает длину вектора состояния, и кроме того, что свойство устойчивости системы можно охарактеризовать с помощью анализа изменения механической энергии.

Скорость изменения энергии в процессе движения системы, фактически движения изображающей точки, можно определить путем дифференцирования

первого равенства в выражении (**) и используя (*),

. (***)

Выражение (***) влечет за собой, что энергия системы, начиная с определенного исходного значения, рассевается до тех пор, пока масса не успокоится, т.е. .

Прямой метод Ляпунова основан на обобщении понятий, связанных с вышерассмотренной системой масса-пружина-демпфер, на более сложные системы.

Понятно, что в случае устойчивого состояния равновесия фазовые траектории системы

, , , (13)

с течением времени «стягиваются» к началу координат (рис. 8.2,а), при этом уменьшается расстояние от текущей точки до точки равновесия.

Для неустойчивого состояния равновесия фазовые траектории системы (13) расходятся (рис. 2,б), а расстояние увеличивается с течением времени.

Таким образом, суть второго метода Ляпунова сводится к оценке изменения некоторой функции переменных состояния системы V(x) , которую называют функцией Ляпунова, вдоль траекторий движения. Обсудим формальные свойства функций V(x), используемые в дальнейшем.

Рис. 2

Рассмотрим функции V(x) определенные и непрерывные в некоторой области пространства состояний D, содержащей начало координат и обладающие в этой области непрерывными частными производными по переменным х.

Функция V(x) называется положительно определенной в области D, если она положительна для любых значений переменных из этой области и обращается в нуль только в начале координат, т. е. выполняются свойства

(14)

Дадим геометрическую интерпретацию положительно определенной функции. Рассмотрим положительно определенную функцию V(x) двух переменных и . В трехмерном пространстве V(x) обычно является поверхностью подобной перевернутой чашке, как показано на рис. ниже. Нижняя точка чашки расположена в начале координат.

Двухразмерная геометрическая интерпретация может быть дана следующим образом. На плоскости с координатами и линии уровня обычно представляют овальные кривые, охватывающие начало координат, причем каждому овалу соответствует постоянное положительное значение . Эти овалы, часто называемые линиями уровня или линиями постоянного значения можно трактовать как проекции сечений чашки горизонтальными плоскостями на плоскость , как показано на рис. ниже.

Функция переменных состояния называется отрицательно определенной в области D, если она отрицательна для любых значений переменных из этой области и обращается в нуль только в начале координат.

Полной производной функции Ляпунова в силу уравнения системы (13) (производной функции Ляпунова вдоль траекторий движения системы) называется функция V(x), которая определяется следующим образом:

. (15)

Здесь - вектор-строка частных производных.

Полная производная функции Ляпунова представляет собой скалярное произведение вектор - строки и вектор - столбца f(x) и может быть представлена в развернутой форме

.

Следует отметить, что полная производная функции Ляпунова , так же как и сама функция Ляпунова V(x), тождественно обращается в нуль в начале координат, поскольку f(0) = 0.