Инвариантность систем

1. Инвариантность систем Инвариант – отображение φ рассматриваемой совокупности М математических объектов, снабженной

фиксированным отношением
эквивалентности ρ, в другую совокупность N математических
объектов, постоянное на классах эквивалентности М по ρ.
Концепция инвариант является одной из важнейших в математике,
поскольку изучение инварианта непосредственно связано с задачами
классификации объектов того или иного типа.
Пример: ранг – инвариант квадратичной формы.
Проблема инвариантности.
Это проблема определения таких структур и параметров систем
управления, при которых влияние некоторых произвольно меняющихся
внешних воздействий и собственных параметров системы на
динамические характеристики процессов управления могут быть
частично или полностью компенсированы.

2. Более простая постановка – требуется сделать по возможности независимой ту или иную переменную (обобщенную координату) от

одного или нескольких внешних воздействий.
В 1938 г. – идея инвариантности была высказана Т.В.Щипановым, а
достаточные и необходимые условия сформулированы Н.Н.Лузиным.
Рассмотрим линейную стационарную систему с тремя
степенями свободы, состоящую из
- объекта регулирования с регулируемой координатой х1(t);
- регулятора с двумя обобщенными координатами х2(t) и х3(t).
а11x1(t)+ а12x2(t)+ а13x3(t)=f1(t) – возмущенное воздействие на объект
а21x1(t)+ а22x2(t)+ а23x3(t)=g(t) – управляющее воздействие
а31x1(t)+ а32x2(t)+ а33x3(t)=f3(t) – возмущение на регулятор
aij=mijp2+lijp+rij, где ; i,j=1, 2, 3
р = d/dt
Допустим, что функция удовлетворяет требованиям
оригинала, и перейдем от дифференциальных уравнений к уравнениям
алгебраическим с помощью преобразования Лапласа.

3. Cтруктура системы

-a32/a33
-a31/a33
f3
1/a33
-a13/a11
-a21/a22
x3
q
x2
-a23/a22
1/a22
f1
a13
1
a12
x1 f1
x2 ( ) x3 ( )
a11
a11
a11
a23
1
a21
x2 g
x1 (
) x3 (
)
a22
a22
a22
a31
a32
1
x1 f 3
x1 ( ) x2 (
)
a33
a33
a33
-a12/a11
1/a11
x1

4. Упорядоченность системы заключается в том, что порядковый номер уравнения соответствует номеру обобщенной координаты, для

которой это уравнение
составлено.
Поэтому элементы главной диагонали операторной матрицы (р)
представляют собой собственные (характеристические) операторы
каждой из обобщенных координат схемы.
Остальные операторы отражают воздействие одних обобщенных
координат на другие.
Найдем условие, при котором регулируемая координата x1(t) не будет
зависеть от внешнего воздействия f1(t).
От оригиналов перейдем к лапласовым изображениям переменных,
x1 (0), x 1 (0), x2 (0)...
учитывая начальные условия в виде
а11(S)x1(S)+ а12(S)x2(S)+ а13(S)x3(S)=F1(S)+r1(S)
а21(S)x1(S)+ а22(S)x2(S)+ а23(S)x3(S)=G(S)+r2(S)
а31(S)x1(S)+ а32(S)x2(S)+ а33(S)x3(S)=F3(S)+r3(S)
где r1,r2,r3 (S) определяются начальными условиями
и согласно исходному уравнению имеет вид
r1 (S ) m1[ x1 (0)s x 1 (0)] l11x1 (0) m2[ x2 (0)s x 2 (0)] ...

5. Решаем систему относительно х1(S)

a11 a12
F1 r1
a13
a21 a22 a23 x1 ( S ) G r2
a13
a22 a23
F3 r3 a23 a33
a31 a23 a33
Раскрывая
a12
определитель
третьего
порядка , запишем
a22 a23
a12 a13
a12 a13
a21
a31
a11
x1
a32 a33
a22 a23
a32 a33
a22 a23
a12 a13
a12 a13
( F1 r1 )
(G r2 )
( F r3 )
a32 a33
a32 a33
a22 a23
11
21
31

6. Условие инвариантности, при котором x1(t) не будет зависеть от возмущения f1(t) (полученное Щипачевым)

11 ( S )
a22
a23
a32
a33
0
Итак
(a21 21 a31 31 ) x1 (G r2 ) 21 ( F3 r3 ) 31
Перейдя к Лапласову
преобразованию для x1 ( S )
при выполнении условий инвариант. по f1 (t )
(G ( S ) r2 ( S )) 21 ( S ) ( F3 ( S ) r3 ( S )) 31 ( S )
х1 ( S )
a21 ( S ) 21 ( S ) a31 ( S ) 31 ( S )
видим, что оказывается исключается только компонента ,
вызываемая действием f1 (t )

7. Теперь положим, что g(t) = 0, f3(t) = 0 и все начальные условия нулевые. Рассмотрим только влияние f1(t) . Очевидно, что при

удовлетворении для
координаты x1(t) условий абсолютной инвариантности все члены вида
ai1(s)x1(s) = 0 при всех i ≠ 1, а это эквивалентно размыканию системы на
выходе элемента, поведение которого характеризуется координатой x1(t).
а11(S)x1(S)+ а12(S)x2(S)+ а13(S)x3(S)=F1(S)
0
+ а22(S)x2(S)+ а23(S)x3(S)= 0
0
+ а32(S)x2(S)+ а33(S)x3(S)= 0
Из уравнений следует, что
p (S )
11 ( S )
x1 ( S )
F1 ( S )
p( S )
- характеристический определитель разомкнутой системы по
координате xi(t)
p(S ) a11(S ) 11(S )
, т.е. при соблюдении условий абсолют. инвар.
( Δ11(S) = 0) матрица системы разомкнутой по x1(t) оказывается особой и
следовательно при f1(t) поведение x1(t) будет описываться особым решением
1
x1 ( S )
F1 ( S ) 0
a11 ( S )

8. Но тогда оказывается, что x1(S) зависит от F(S), т.е. реализация абс. инвар. оказалась невозможной. Это объясняется тем, что в

системе
имеется лишь один канал распространения сигнала от точки приложения
воздействия f1(t) к точке измерения x1(t) .
По аналогии рассмотрим реализацию абсолютной инвариант.
координаты x1 от f3(t).
31 (S )
a12
a13
a22
a23
а11(S)x1(S)+ а12(S)x2(S)+ а13(S)x3(S)=0
0 + а22(S)x2(S)+ а23(S)x3(S)= 0
0 + а32(S)x2(S)+ а33(S)x3(S)= F3(S)
x1 ( S )
31 ( S )
F3 ( S )
p( S )
Но при Δ31(S) = 0 характерист. определитель разомкнутой
системы Δp(S) не обращается в нуль и система уравнений не
становится особой.

9. Это происходит (как следует из схемы) благодаря тому, что между f3 и х1 есть два канала, и тогда возможна компенсация. Это

условие известно как принцип
двухканальности (Петров).
Это необходимое условие физической реализуемости абсолютной
инвариантности.
Критерии реализуемости условий абсолютной инвариантности:
Необходимое условие реализуемости инвариантности переменной хi(t) по
отношению к внешнему воздействию fj(t) является тождественное совпадение
множества решений исходной системы и системы разомкнутой на выходе
элемента хi(t) в точке измерения при выполнении условия инвариантности и
равенстве всех остальных воздействий и равенстве начальных условий.
Для достижения абс. инвар. некоторой координаты xi(t) относительно fi(t)
необходимо и достаточно, чтобы передаточная функция между точкой
приложения внешней силы и точкой измерения была тождественно равна 0,
когда все остальные воздействия отсутствуют и нулевые нач. условия.
Для получения абс. инвар. как в разомкнутом, так и в замкнутом состоянии
необходимо удовлетворить требование Δij(S) = 0 , а также потребовать, чтобы
Δ(S) ≠ 0 для разомкнутой и замкнутой системы.
Абс. инвар. не может быть реализована в том случае, если матрица системы
ДУ в разомкнутом состоянии становится особой после учета условий абс.
инвариантности.

10. Устойчивость движения Под устойчивостью функционирования сложной системы понимают способность системы сохранять требуемые

свойства в условиях
действий возмущений.
Рассматривая нелинейные системы вводят понятие устойчивости
«в малом», «в большом», «в целом».
Система устойчива «в малом», если констатируют лишь факт
наличия области устойчивости, но не определяют каким-либо образом
ее границы.
Систему называют устойчивой «в большом», когда определены границы
области устойчивости, т.е. определены границы области начальных
отклонений, при которых система возвращается в исходное
состояние.
Когда система возвращается в исходное состояние при любых
начальных условиях, систему называют устойчивой «в целом».
Устойчивость «в целом» для определенного класса нелинейностей
называют «абсолютной устойчивостью».

11. Постановка задачи. y1…yn - вещественные переменные, характеризующие состояние системы. Y1…Yn – известные функции,

Постановка задачи.
dy1
Y1 ( y1 , y2 ,..., yn , t )
dt
......................................
dyn
Yn ( y1 , y2 ,..., yn , t )
dt
(1)
y y (0)
y1…yn - вещественные переменные, характеризующие состояние
системы.
Y1…Yn – известные функции, удовлетворяющие условию существования
и единственности решения.

12. Невозмущенное движение – некоторое вполне определенное движение системы, подлежащее исследованию на устойчивость. Возмущению

подвергаются только начальные условия.
Невозмущенному движению системы отвечает определенное
частное решение ДУ (1)
y1 f1 (t ) ... y n f n (t )
(2)
получили такое реш ение.
При t t 0 начальноеусловие
y1 f1 ( t 0 ) ... y n f n ( t 0 )
(3)
Дадим начальным значениям некоторое приращение ε
y1 f1 (t 0 ) 1
.........................
y n f n (t 0 ) n
(4)
Движение системы, отвечающие измененным начальным условиям
(4), есть возмущенное движение, а ε1… ε n - возмущения.

13. Обозначим yj – возмущенное движение; fj - невозмущенное движение; xj - отклонение или вариация xj=yj(t) – fj(t) (j =1…n) (5)

Обозначим
yj – возмущенное движение;
fj - невозмущенное движение;
xj - отклонение или вариация xj=yj(t) – fj(t) (j =1…n)
Если х1=0,…,хn=0
то возмущенное движение совпадает с невозмущенным
х1
0
х2
(5)
(6) ,
Геометрически можно интерпретировать так:
совокупность отклонений в n–мерном
пространстве переменных x1…xn определяет
М(х1, х2, х3) точку М (изображающая точка).
В возмущенном движении при изменении в-н
x1…xn , М будет описывать некоторую траекторию.
х3
Невозмущенному движению xj = 0 отвечает
неподвижная точка – начало координат.
n
2
Мера отклонения: х х ... х x j (7)
2
1
2
2
2
n
j 1
При t = t0 xj=x0j=εj (j=1..n), т.е. начальные значения отклонений xoj
представляют возмущения системы
(8)

14. Определение устойчивости движения по Ляпунову. Если по любому положительному числу ε, как бы оно не было мало, можно найти

такое положительное число δ, что при всяких возмущениях x0j, удовлетворяющих
условию х02 j
(9) и при любых t 0 будет
x
определяться неравенство
(10),
то невозмущенное движение устойчиво, в противном случае нет.
2
j
Практически устойчивость данного невозмущенного движения означает,
что при достаточно малых начальных возмущениях, возмущенное движение
будет сколь угодно мало отличаться от невозмущенного.
Если же невозмущенное движение неустойчиво, то возмущенное
движение будет отходить от него, как бы малы ни были начальные возмущения.
Если невозмущенное движение устойчиво и при этом любое возмущенное
движение при достаточно малых начальных возмущениях стремится к
невозмущенному движению, т.е lim x 2j (t ) 0
(11),
t
то невозмущенное движение называется асимптотически устойчивым .

15. Рассмотрим сферу Выберем радиус √ε произвольно малым. Если движение устойчиво, то для этой сферы должна найтись другая сфера ,

Рассмотрим сферу
М0
●0
δ
М
ε
x
2
j
Выберем радиус √ε произвольно малым.
Если движение устойчиво, то для этой сферы
должна найтись другая сфера
x 2j ,
обладающая следующим свойством.
Изображающая точка М, начав свое движение из любого положения М0, лежащего
внутри или на поверхности сферы δ, при своем дальнейшем движении остается
всегда внутри сферы ε, никогда не достигая ее поверхности.
Если же невозмущенное движение неустойчиво, то хотя бы одна траектория
изображающей точки М с течением времени пересечет сферу ε изнутри наружу при
сколь угодно близком положении точки М0 к началу координат.
Геометрически это означает, что при асимптотической устойчивости
изображающая точка должна неограниченно стремится к началу координат,
не выходя из сферы ε.
В тех случаях, когда асимптотическая устойчивость имеет место при любых
возмущениях (не обязательно малых), невозмущенное движение называется
асимптотически устойчивым в целом.
Иногда устойчивость имеет место не при любых возмущениях, а при возмущениях,
подчиненных некоторым условиям. Такая устойчивость называется условной.

16. Особенности определения устойчивости по Ляпунову. 1. Возмущения накладываются только на начальные условия, что физически

говорит о том, что возмущенное движение происходит при тех же
источниках энергии, что и невозмущенное.
2. Устойчивость рассматривается на бесконечно большом интервале времени.
3. Возмущения предполагаются малыми.
Тем не менее, методы развитые Ляпуновым лежат в основе
исследования других видов устойчивости движения.
Существуют два метода Ляпунова:
1. Оценка устойчивости по приближенному решению – основан на
линеаризации.
2. Прямой метод Ляпунова – осуществляется через функцию Ляпунова.

17. Составим уравнения возмущенного движения yj(t) = fj(t) + xj(t) Подставим в уравнение (1) где - совокупность членов, зависящих

Составим уравнения возмущенного движения
yj(t) = fj(t) + xj(t)
Подставим в уравнение (1)
df j
dx j
Y j ( f1 x1 , f 2 x 2 ,..., f n x n , t )
dt
dt
Разложим правые части этих уравнений в ряды Тейлора по степеням
х j (предполагается, что фукнции ра злагаемые в ряды удовлетворяют
соответствующим требованиям)
Y j
Y j ( f1 , f 2 ,..., f n , t )
dt
dt
x1
df j
dx j
Y j
x1
0
x 2
Y j
x 2 ...
0
x n
x n X *j
0
где
X -*j совокупность членов, зависящих от отклонений xi в степени
выше первой. Учтем, что в невозмущенном движении функции fj(t) должны
удовлетворять уравнению (1), т.е.
df j
dt
Y j ( f1 , f 2 ,..., f n , t )

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

Тогда диф. уравнения возмущенного движения
Y j
xi X *j
dt
i 1 xi 0
dx j
n
j 1...n ( 12 )
a ji
коэф. a jk (
Y j
xk
)x 0 ( 13 )
в общем случае являются функциями времени, в частности могут быть
постоянными.
Если в уравнениях (12) отбросить члены X *j , то полученные при этом
уравнения называются уравнениями первого приближения.
dx j
dt
a j1 x1 a j 2 x2 ... a jn xn (14)
Ур-ия первого приближения во многих случаях дают верный ответ
на вопрос об устойчивости движения, но иногда заключение, которое
можно получить из этих приближенных уравнений ничего общего не
имеет с решением исходных уравнений

19. Пример. Пусть уравнения возмущенного движения имеют вид Умножим первое уравнение на х1, второе на х2 и сложим почленно оба

Пример.
Пусть уравнения возмущенного движения имеют вид
dx1
x2 x1 x12 x22
dt
dx2
x1 x2 x12 x22
dt
Умножим первое уравнение на х1, второе на х2 и сложим почленно оба ур-ия:
dx1
dx2
x1
x2
( x12 x22 )3 / 2 или
dt
dt
1 d 2
( x1 x22 ) ( x12 x22 ) 3 / 2
2 dt
Положим x12+x22=r2 , где r – расстояние от начала координат до изображающей
точки. После перехода к новой переменной r имеем

20. При α>0 r неограниченно возрастает при t→t0+1/αr0 [1-αr0(t-t0)=0 αr0t=1+αr0t0 t=t0+1/αr0 , т.е. при t=t0+1/αr0 знаменатель

1 d 2r
dr
r 3 или
r 2
2 dt
dt
r
t
r
t
dr
1
Интегрируем 2 dt r
t
r0
t0
r
r0
t0
1 1
(t t 0 ) r0 r rr0 (t t0 )
r r0
r0
r
1 r0 (t t0 )
При α>0 r неограниченно возрастает при t→t0+1/αr0
[1-αr0(t-t0)=0 αr0t=1+αr0t0 t=t0+1/αr0 , т.е. при t=t0+1/αr0 знаменатель
равен 0 и r-неограниченно возрастает], т.е. движение неустойчиво.
Это пример того, что одного предельного условия
асимптотической устойчивости недостаточно,
необходимо, чтобы движение было устойчивым.
lim r 2 0
t

21. При α<0 r монотонно убывает, стремясь к нулю при t→∞, т.е. движение асимптотически устойчиво. Рассмотрим теперь уравнения

При α<0 r монотонно убывает, стремясь к нулю при t→∞, т.е. движение
асимптотически устойчиво.
Рассмотрим теперь уравнения первого приближения, которые
получаются отбрасыванием членов порядка выше первого
dx1
x2
dt
dx2
x1
dt
для этих ур ий имеем
dr
0
dt
или , интегрируя, r r0
Это решение показывает, что изображающая точка М,
отвечающая уравнениям первого приближения, движется по окружности,
радиус которой равен начальному отклонению т.М от начала координат.
Т.о. из уравнений первого приближения следует устойчивость
невозмущенного движения при всех α (х1=х2=0),
но этот вывод не совпадает с результатами анализа исходных
уравнений.

22. Прямой метод исследования устойчивости Так называют второй метод Ляпунова, который позволяет судить об устойчивости

непосредственно по уравнениям возмущенного движения, не
прибегая к их интегрированию.
Уравнения возмущенного движения записываются в нормальной
форме Коши в отклонениях от невозмущенного движения
dxi/dt = Xi(x1,x2…xn) i=1…n.
Невозмущенному движению при этом соответствует тривиальное
решение, т.е. x1=x2=..xn=0.
Если Xi(x1…xn) представляют собой функции фазовых координат
(отклонений), непрерывные в некоторой области Rn, содержащей начало
координат и имеющей частные производные по всем аргументам, то для
анализа устойчивости невозмущенного движения могут быть использованы
специальные функции фазовых координат , называемые функциями
Ляпунова.
Прямой метод опирается на известную теорему Лагранжа, согласно
которой равновесие устойчиво, если в положении равновесия потенциальная
энергия системы минимальна.
Идея метода состоит в том, чтобы подобрать такую функцию фазовых
координат, которая бы в некотором смысле была аналогична потенциальной
энергии системы в состоянии покоя.

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

области Rn Σxj2≤ μ (15) ( μ- постоянное положит. число), удовлетворяющие
условию V(x1…xn)=0 при x1=x2=…xn=0 (16) и имеющие производные по всем
аргументам.
Цель состоит в том, чтобы, предполагая невозмущенное движение
устойчивым, попытаться подобрать такую функцию фазовых координат,
которая при любом движении системы уменьшалась, т.е dV/dt<0.
Если в окрестности начала координат функция V кроме нуля может
принимать значения только одного знака, то она называется
знакопостоянной (положительной или отрицательной)
Если знакопостоянная ф-я обращается в нуль только в том
случае, когда все x1…xn равны нулю, то ф-я V называется
знакоопределенной (определенно-положительная или определенноотрицательная)
Знакоопределенная функция имеет при x1=…xn=0 экстремум (min
для опред.-положит. функции и mах для опред.-отриц).
Знакопостоянная функция в начале координат экстремума не имеет, т.к. в
окрестности начала координат есть другие точки, в котoрых V=0.

24. Пусть V=V(x) непрерывна вместе с производными первого порядка: кроме того предположим, что V(x) знакоопределенная. Тогда при

x1=…xn=0 она будет иметь изолированный.
экстремум и все
V
(
) 0 0 j 1....n (17)
x j
Разложим V в ряд Маклорена по степеням x1…xn
V
1 n n 2V
V V (0) (
) 0 x j (
) xk x j ...
2 k 1 j 1 xk x j
j 1 x j
n
Учитывая (16) и (17), получим
1 n n 2V
V (
) xk x j ... (19)
2 k 1 j 1 x j xk
Обозначим
2V
rkj (
) 0 и rkj r jk или
x j xk
1 n n
V rkj xk x j ...
2 k 1 j 1
(20)
(18)

25. Т.о. разложение знакоопределенной функции V в ряд по степеням x1…xn не содержит членов первой степени, т.е. остается

квадратичная
форма:
1 n n
rkj xk x j
(21)
2 k 1 j 1
Пусть квадратичная форма принимает положительные значения и в нуль
обращается только при x1=…xn=0 .
Тогда вне зависимости от членов высшего порядка при достаточно
малых по модулю xj функция V будет принимать тоже положительные
значения и в нуль обращается только при x1=…xn=0 .
Если квадратичная форма (21) определенно – положительна, то и
функция будет определенно положительной.
Рассмотрим матрицу коэффициентов
квадратичной формы (21) и составим из нее
r11 r12 ... r1n
n главных диагональных миноров.
r
r21
r22
... r2 n
...
...
...
rn1
rn 2
... rnn
...
( 22)

26. В линейной алгебре доказывается следующий критерий Сильвестра: Для того, чтобы квадратичная форма была

определенноположительной, необходимо и достаточно, чтобы все главные
диагональные миноры матрицы ее коэффициентов были
положительны, т.е. Δ1>0, Δ2>0… Δn>0
1 c11 , 2
c11
c12
c21 c22
c11 ... c1n
... n ...
...
(23)
cn1 ... cnn
Для того, чтобы квадратичная форма была определенноотрицательной, необходимо и достаточно, чтобы имели место
неравенства Δ1<0, Δ2>0, Δ3<0 , т.е. oпределители должны
последовательно чередовать знак, причем знак Δ1=c11 д.б.
отрицательным.
Критерий Сильвестра для квадратичной части функции V является
достаточным (но не необходимым) условием определенной
положительности самой функции V.

27. Пример. Знакоопределенная функция Знакопостоянная функция V V x2 x2 x1 x1

Пример.
V ( x) 2 x12 4 x24
V ( x) x12 2 x1 x2 x22 ( x1 x2 ) 2
Знакоопределенная функция
Знакопостоянная функция
V
V
x2
x1
x2
x1
V
2V
1
x j
V ( x ) V0
j 1 x
x x
2
j
1
k
1
j
j k
0
n
n
0
1 n n
V
( x) r jk x j x k
2 j 1 k 1
квадратичная
форм а
n
2
V
x j x k ...
r jk
x j x k

28. Может оказаться, что разложение знакоопределенной функции V в ряд по степеням x1…xn начинается не с членов второго, а с членов

более высокого порядка. В этом
случае общих приемов исследования функции на знакоопределенность нет.
Некоторые св-ва функции V
Если функция V знакоопределенная, то поверхность V(x1…xn) - замкнута. Выберем на
поверхности V(x)=c произвольную точку М и вычислим в этой точке вектор gradV
e1…en - орты осей x1…xn.
gradV
V
V
V
e1
e2 ...
en
x1
x2
xn
(24)
Известно, что вектор gradV напрaвлен по нормали к поверхности V=c в точке М в
сторону возрастания функции V.
Из этого следует, что вектор gradV направлен во внешнюю часть поверхности V=c,
если функция V определена положительно и внутрь V=c, если V определена
отрицательно.
gradV
gradV
M
M
V=c>0
V=c<0

29. Одновременно с функцией V будем рассматривать ее полную производную по времени t, взятую в предположении, что переменные xj

удовлетворяют
диф.ур. возмущенного движения dxj/dt=xj (j=1…n)
(25)
Имеем
dV
V
V
V
V
x 1
x 2 ...
x n
dt
x1
x2
xn
( 26)
или с учетом (25)
V
V
V
V
X1
X 2 ...
Xn
x1
x2
xn
( 27)
Напомним, что величины Xj равны проекциям скорости U
изображающей точки М, а производные
V
x j
– проекциям gradV.
Поэтому правая часть равенства (27) равна скалярному
произведению векторов U и gradV , т.е.
V U gradV

30. Теоремы Ляпунова об устойчивости 1. Если при заданных уравнениях возмущенного движения системы можно найти такую

знакоопределенную функцию V, полная производная
которой является функцией знакопостоянной противоположного знака
по отношению к V или тождественно равна нулю, то невозмущенное
движение устойчиво.
2. Если при заданных уравнениях возмущенного движения системы
можно найти такую знакоопределенную функцию V, полная производная
которой по времени в силу этих уравнений является функцией
знакоопределенной противоположного знака по отношению к V, то
невозмущенное движение устойчиво асимптотически.
Рассмотренные теоремы Ляпунова дают достаточные условия
устойчивости, т.е. невыполнение этих условий не означает, что
невозмущенное движение неустойчиво. Полнота определения области
устойчивости в пр-ве фазовых координат зависит от выбора конкретной
функции.
Полученные условия устойчивости могут не охватывать всей
области устойчивости системы по параметрам.

31. В нелинейных системах, в отличие от линейных, возможны случаи, когда невозмущенное движение устойчиво при «малых» отклонениях

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

32. Пример Рассмотрим функцию V=1+sin2x1-cos(x1-x2) Разложим эту функцию по степеням x1 и x2 sin2x1=x12+…

cos(x1-x2)=1-1/2(x1-x2)2+…
V=1+x12-1+1/2(x1-x2)2+…= 1/2(3x12-2x1x2+x22)+…
Составим матрицу коэффициентов квадратичной части функции V (…
обозначены члены выше второй степени). По главной диагонали стоят
коэф. при квадратах переменных, элементы с12, с21 равны половине коэф. при
3 1
произв. х1х2)
1
1
Вычислим главный диагональный
минор
3 1
1 3, 2
1
1
2
т.к. все Δj>0 , то неравенство Сильвестра выполнено и рассматриваемая
функция V будет определенно-положительной.

33. Пример Вычислим полную производную Пусть функция соответствует уравнениям возмущенного движения т.е. производная V

Пример
Вычислим полную производную
V (3 x1 x2 ) x 1 ( x1 x2 ) x 2
V
x1
V
x2
Пусть функция соответствует уравнениям возмущенного движения
1
x1 x22
2
1
x2 3x2 x1 x2 x12 x2 x1 x22
2
1
2
a V (3 x1 2 x1 x2 x22 )
2
V 3 x14 2 x12 x2 2 x22
x 1 x2 x1 x2 x13
1 3 0, 2
3
1
1
2
5 0
т.е. производная V определенно-отриц. и движение асимптотически
устойчиво