Критерий Рауса-Гурвица.
Пусть Р =- произвольная квадратичная матрица порядка n
Назовем главными минором детерминант (определитель) матрицы:
Пусть L(p) = а0∙р(n) + a1×р(n-1) + …+ an-1×р + an - произвольный многочлен с действительными (вещественными) коэффициентами
Q = - порядка n
Многочлен L(p) устойчив, когда все главные миноры положительные
k + j –отрицательные > n, то коэффициенты считаются = 0
§ Геометрический критерий устойчивости (Критерий Михайлова).
Допустим, есть однородное линейное уравнение:
ao·y(n) + a1·y(n-1) + … + an-1·y' + an·y = 0
введем характеристическое уравнение:
¦ (p) = ao·pn + a1·pn-1 + … + an-1·p + an = 0
вместо р введем р = i·w, i =
¦ (p) = ¦ (i·w) = u(w) + i·v(w), где
u(w) = an – an-2·w2 + an-n·wn – …
(*)
v(w) = an-1·w – an-3·w3 – an-5·w5 – …
¦ (i·w) при заданных значениях w и (*) можно изобразить на комплексной плоскости:
годограф Михайлова - ¥ < w < +¥
По виду годографа Михайлова можно судить об устойчивости системы.
Обозначим корни многочлена ¦ (p) через р1, р2,…, рn. Тогда можно записать:
|
|
¦ (p) = ao·(р – р1)·(р – р2)·…·(р – р2)
Тогда ¦ (p) = ¦ (i·w) = ao·(i·w – р1)·(i·w – р2)·…·(i·w – р2)
Комплексная плоскость корней:
S = a + i·w - вещественная часть
Подсчитаем угол j, на который повернется вектор ¦ (i·w) при изменении - ¥ < w < +¥.
При перемножении комплексных чисел их аргументы складываются, т.е. j = j1 + j2 +…+ jn
¦ (p) - многочлен с m корней с положительными вещественными частями
n – m корней слева от мнимой оси
j = (n – m)·p + m·(-p) = (n-2m)p
Замечание: u(w) – четная функция, т.е. можно изучать поведение корней только 0 ≤ w < +¥,
а - ¥ < w ≤ 0 – зеркальное отражение.
Тогда j = (n – 2·m)· при 0 ≤ w < +¥,
ao·y(n) + a1·y(n-1) + … + an-1·y' + an·y = 0
Чтобы система была устойчивой, необходимо и достаточно, чтобы все корни уравнения ¦ (p) = ¦ (i·w) = 0 имели отрицательные вещественные части. Т.е. для устойчивой системы m = 0.