Toward Robust Stability Margin Evaluation in Continuous and Discrete Systems

Research article
DOI:
https://doi.org/10.23670/IRJ.2023.129.22
Issue: № 3 (129), 2023
Suggested:
08.01.2023
Accepted:
22.02.2023
Published:
17.03.2023
811
1
XML
PDF

Abstract

The work presents the problem of determining the robust stability margin, i.e. the stability margin in conditions of uncertainty in setting all or part of parameters of continuous and discrete dynamic systems. A procedure for determining the maximum stability margin, based on the introduction of the parameter determining the value of the stability margin into the characteristic equation in explicit form, is suggested. When determining the robust stability of reduced polynomials, the Routh–Hurwitz theorem for continuous dynamical systems and the Korsakov criterion for discrete dynamical systems, as well as the basic statements of the theory of robust stability, which allows to take into account the initial parametric uncertainty, are used. Examples of finding the maximum robust stability margin in concrete continuous dynamical systems with the second and third order characteristic polynomials are presented.

1. Введение

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

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

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

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

Актуальность указанного подхода связана прежде всего с задачами, где речь идет не только о качественных исследованиях, но и получении некоторых количественных характеристик интервально неопределенных как непрерывных, так и дискретных систем. В частности с задачей получения оценок длительности переходных процессов в интервально неопределенных системах.

2. Запас устойчивости непрерывных систем и метод Рауса-Гурвица

Будем характеризовать запас устойчивости величиной img, для которой корни семейства полиномов

img
(1)

удовлетворяют условию img. Для того чтобы свести задачу определения img к проблеме Гурвица перейдем на плоскости корней к новому значению img, соответствующему сдвигу мнимой оси на величину img влево. В этом случае полином (1) переходит в полином

img
(2)

коэффициенты которого определяются из тождества

img
(3)

и равны

img
(4)

Известно

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

img
(5)

при условии, что при img. Следовательно, для определения максимального значения параметра img необходимо анализировать выражение img при увеличении img, начиная со значения img. Представим при img выражение для img в виде разности

img
(6)

где img – члены правой части (5), содержащие img в четных степенях, включая ноль, а img – члены, содержащие img в нечетных степенях. Тогда очевидно, что

img
(6)

где

img
(7)
img
(8)

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

Представим img в виде разности положительных членов

img
(9)

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

img
(10)

Здесь черта у img и img переносится однотипно на все входящие в них img. Тогда выполнения неравенства img достаточно для утверждения, что img и, следовательно, при img комплексно-сопряженные корни семейства полиномов (1) не могут находиться на прямой img.

Для оценки величины img на основе оценки знакоположительности img можно, скажем, считать, что img. Тогда условие знакопостоянства img можно установить на основе анализа величины img, где

img
(11)

а черта у img и img переносится однотипно на все входящие в них img и параметр img.

3. Алгоритм оценивания максимального значения робастного запаса устойчивости и его применение к непрерывным системам с характеристическими полиномами второй и третьей степени

Как следует из приведенных выше рассуждений, при img достаточными условиями робастной устойчивости полинома (2) являются

img
(12)

Предположим, что при img для полинома img в плоскости каких-либо двух выбранных параметров img при интервальности других построена область робастной устойчивости

. Для оптимизации качества переходного процесса в системе необходимо определить такое img, чтобы выполнялись соотношения (12). Так как img и img являются полиномами относительно img, процесс определения img можно провести следующим образом:

1) Находим аналитически или с помощью известных алгоритмических процедур корни уравнений img и img. Обозначим их через img и img соответственно.

2) Образуем последовательность действительных положительных чисел img из этих корней.

3) Определяем img.

4) Строим область робастной устойчивости

в плоскости параметров img при интервальности всех остальных параметров.

1. Пусть характеристическое уравнение img имеет вид

img
(13)

После введения величины img будем иметь

img
(14)

Для определения img имеем два уравнения img

img
(15)

Из (15) определяем

img
(16)

Предположим, что img. Тогда достаточная граница робастной устойчивости, оптимизирующая качество переходного процесса в системе, будет определяться из уравнения

img
(17)
в котором выполняется достаточное условие робастной устойчивости полинома img, в качестве параметров для выделения области робастной устойчивости выбраны img, а параметр img остался интервальным. Тогда область робастного запаса устойчивости, оптимизирующего качество переходного процесса будет иметь вид, представленный на рисунке 1. Следует отметить, что если при img областью робастной устойчивости был весь квадрант img, то в данном случае ей является только часть этого квадранта.
Область робастного запаса устойчивости полинома второй степени

Рисунок 1 - Область робастного запаса устойчивости полинома второй степени

Аналогично можно выделить область робастного запаса устойчивости в случае, когда img и img.

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

.

2. Предположим теперь, что характеристический полином img имеет третью степень, т.е. характеристическое уравнение будет иметь вид

img
(18)

Такая ситуация имеет место в математической модели автоколебаний в основной стойке шасси самолета в случае бесконечной жесткости на кручение без учета момента инерции ориентирующейся части стойки относительно оси стойки

. Полином img при этом будет иметь вид

img
(19)

где

img
(20)

Уравнения для определения img имеют вид

img
(21)
img
(22)

В безразмерных переменных img для системы из

могут изменяться в пределах img. Тогда, если предположить, что минимальный из img определяется корнем первого уравнения, то, согласно (20), (21), img. А граница, определяющая область робастного запаса устойчивости, задается уравнением

img
(23)

где в качестве параметров для выделения области робастной устойчивости выбраны img, а остальные параметры остались интервальными. Уравнение границы (23) можно записать в виде

img
(24)

а если учесть ограничения на параметры, получим

img
(25)

В плоскости выбранных параметров img график функции (25) будет иметь вид, приведенный на рисунке 2. Из (25) и рисунка 2 видно, что область робастного запаса устойчивости будет располагаться внутри приведенной кривой.

Область робастного запаса устойчивости полинома третьей степени

Рисунок 2 - Область робастного запаса устойчивости полинома третьей степени

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

Запас устойчивости дискретных динамических систем может быть охарактеризован параметром img таким, что img, где img – корни характеристического полинома (1). Поскольку устойчивость (1) имеет место при img, задача обеспечения заданного запаса устойчивости (или определение минимально возможного img) является более общей задачей по отношению к определению условий устойчивости характеристического полинома. Параметр img можно ввести формально в (1) за счет замены img. В этом случае при любом img корни полинома должны лежать в единичном круге, т.е. так, как это имеет место при исследовании устойчивости. При таком подходе

img
(26)

где

img
(27)

При использовании критерия Корсакова как для проверки устойчивости img относительно единичного круга при img, так и для нахождения возможного значения img такого, что при img img устойчив, можно ограничиться анализом условий невыхода в пространстве img за границы области устойчивости. В этом случае расположение в сторону области устойчивости от границ img соответствует выполнению неравенств img и img, где img – матрицы Корсакова полинома img.

Проверка условия img на семействе очевидна, т.к для этого необходимо и достаточно, чтобы выполнялось неравенство img, где

img
(28)

Для анализа условия img достаточно, чтобы имело место img, где

img
(29)

Процедура определения img при попадании на границы img сводится к определению максимального корня img, при котором обращается в ноль правая часть (28) или правая часть (29).

При рассмотрении выполнимости условия img с квадратными матрицами img-го порядка вида

img
(30)

после замены img основная сложность состоит в нахождении выражения для img. Так в случае img

img
(31)

а значит img, где

img
(32)

и условие img является достаточным для выполнения img. А в случае img

img
(33)

т.е. представляет собой уже полином от img достаточно высокой степени, и получение оценки img уже является затруднительным, хотя и возможным

img
(34)

И получение достаточных оценок тем сложнее, чем больше img.

5. Заключение

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

Article metrics

Views:811
Downloads:1
Views
Total:
Views:811