РАДИОТЕХНИКА И ЭЛЕКТРОНИКА, 2015, том 60, № 7, с. 731-741
ТЕОРИЯ РАДИОТЕХНИЧЕСКИХ ЦЕПЕЙ
УДК 621.373
АЛГЕБРАИЧЕСКИЕ КРИТЕРИИ УСТОЙЧИВОСТИ КОМПЛЕКСНЫХ ПОЛИНОМОВ И ИХ ПРИМЕНЕНИЕ В РАДИОЭЛЕКТРОНИКЕ © 2015 г. В. М. Богачев
Национальный исследовательский университет "Московский энергетический институт" Российская Федерация, 111250, Москва, ул. Красноказарменная, 14 E-mail: BogachevVM@mpei.ru Поступила в редакцию 18.07.2013 г.
Предложен новый подход к обоснованию ряда алгебраических критериев устойчивости и теорем о локализации корней комплексного полинома на комплексной плоскости. Решение основано на теоремах о связи комплексных полиномов с методами синтеза реактивных и квазиреактивных двухполюсников. Рассмотрены аналоги табличного критерия Рауса и детерминантный критерий Эрми-та—Гурвица, включая критические случаи. Показано, что непревзойденным по простоте и точности вычислений является алгоритм Рауса. Применение критериев проиллюстрировано на моделях реальных динамических систем: двух- и трехконтурных генераторов с близкими собственными частотами.
DOI: 10.7868/S0033849415060029
ВВЕДЕНИЕ
Алгебраические критерии устойчивости полиномов с комплексными коэффициентами разработаны Эрмитом (1856 г.) и адаптированы для инженерных приложений Гурвицем (1895 г.). Наиболее известны детерминантный критерий Эрмита— Гурвица [1—3] и его иннорное представление, предложенное Э. Джури [4]. Однако для решения технических задач упомянутые критерии, как правило, не применяют, поскольку полученные строгими методами характеристические уравнения линеаризованных динамических систем имеют вещественные коэффициенты. В то же время комплексные характеристические уравнения получаются при анализе методом укороченных операторных уравнений [5—8] устойчивости положений равновесия и периодических режимов узкополосных динамических систем [8—15], при исследовании стационарных режимов автоматических систем [16] и в ряде других случаев, связанных со специальным преобразованием исходных уравнений системы или изменением области локализации корней [2—4]. Область приложения рассматриваемых методов существенно расширилась в связи с открытием нового направления радиоэлектроники — проектирования "комплексных" цифровых и активных ^С-фильтров, для которых точные характеристические уравнения являются комплексными [17—19]. Таким образом, перед критериями устойчивости и локализации корней комплексных полиномов открыта широкая область технических приложений.
Цель работы — привести новое, базирующееся на теории цепей, обоснование группы алгебраических критериев устойчивости для комплексных полиномов, проиллюстрировать их применение на моделях реальных радиоэлектронных устройств, а также изучение эффективности рассматриваемых критериев, включая критические случаи.
1. КРИТЕРИЙ УСТОЙЧИВОСТИ ТИПА РАУСА-КАУЭРА
Критерий основан на разложении ассоциированной с комплексным характеристическим полиномом H(p) обобщенной реактивной функции в цепную дробь по методу типа Кауэра. Аналогичный критерий, известный в теории цепей для вещественных полиномов, условимся называть критерием Рауса-Кауэра, а его обобщение — критерием типа Рауса—Кауэра.
По определению цепь с входным иммитансом W(p) является реактивной, если W(p) — рациональная неприводимая функция с комплексными (в общем случае) коэффициентами и выполняются два условия: цепь устойчива при подключении единичного резистора и Re W (j ю) = 0 при всех ю. Цепь, удовлетворяющую второму, но не удовлетворяющую первому условию, назовем квазиреактивной.
Для доказательства критерия используем теорему о связи комплексных полиномов Гурвица с обобщенными реактивными функциями [11].
Xi/j X/J Xn/j
—ПИ— -пп— — ГШ—
C1 = = Ci = Cn
Форма 1А
(а)
L: Ln
Форма 1В
Форма 1С
(б)
Ci C:
HI"" Т-II'"
JBift jBi 0 jBn
Cn
jX1
•-ГГП-
Li
C2 JX3
HI—tm-i
JB2
Форма 2А Форма 2В Форма 2С
Рис. 1. Канонические формы типа Кауэра: аналоги 1-й и 2-й форм Кауэра.
Ln
Теорема 1. Комплексный полином Н(р) тогда и только тогда асимптотически устойчив (т.е. является полиномом Гурвица), когда ассоциированный с H(jю) = Нв(ю) + у'Нм(ю) иммитанс W(jю) = = Нв(ю)/у'Нм(ю) является обобщенной реактивной функцией, причем вещественная и мнимая части полинома Н(jю) не имеют общего множителя, т.е. дробь Нв/ jHM не может быть сокращена.
Относительно теоремы сделаем следующие замечания.
Если вещественная и мнимая части полинома Н (j ю) имеют общий множитель, то он может быть только вещественной функцией частоты: F (ю). В этом случае Н (j ю) можно представить в виде произведения двух полиномов: Н(jю) = Н(jra)F(ra), где Н(jю) — гурвицева часть полинома порядка m< n, F(ю) — вещественный множитель порядка v = n — m.
Вещественный полином, в частности F (ю), может иметь только вещественные и комплексно -сопряженные корни. Вещественный по переменной ю корень располагается на мнимой оси p-плоскости, а каждой комплексно-сопряженной паре = а ± ур соответствует пара корней pk i = ±р - ja, один из которых расположен в левой, а другой — в правой части плоскости p. Таким образом, для решения задачи о типе корней полинома F(p) достаточно определить число вещественных корней Ув (с учетом кратности) полинома F (ю), остальные распределятся поровну слева и справа от мнимой оси: L = R = (v - VJ/2. Полином Н(p) устойчив, но не асимптотически, если Н(р) —
полином Гурвица и все корни Р (ю) вещественные простые, иначе полином Н(р) неустойчив.
Рассмотрим сначала комплексный полином Н(р), не имеющий общего вещественного множителя. Проведем замену переменной р = у ю и, разделив Н (у ю) на вещественную и мнимую части, получим импеданс
Z(jЮ) = - aom" + a1®" +••• + an
Н м j(bo®" + ¿1®" + • • ■ + bn)
(1)
Для проверки реактивности импеданса 2 (у ю) можно воспользоваться одним из методов синтеза обобщенных реактивных цепей 1-й или 2-й формы типа Кауэра, каждая имеет три разновидности: А, В и С [11, 15]. Соответствующие цепи показаны на рис. 1. Они содержат "обычные" реактивные элементы Ь, С(- и "псевдореактансы" уХк, ]В1 — частотно независимые реактивные элементы.
Используем форму 1А и проследим процедуру синтеза шаг за шагом. Математически процедура соответствует записи 2 (у ю) в виде цепной дроби.
На первом шаге, выполнив в соответствии с алгоритмом Евклида деление числителя на знаменатель: Нв/ уНм, выделим псевдореактанс Х1 = а0/Ь0 и запишем результат в виде суммы 2 (у ю) = Х1/у + 21(у ю). Здесь 21(ую) — "остаток" от деления Н„/]Нм :
п-1 , п-2 ,
7(. ч_ао® + а1ю +... + ап-1 21(у ®) — ■
j(b0&" + ¿1ю
n—1
• + bn)
где ak = ak+1 — X1bk+1 — коэффициенты числителя (k = 0,...,n - 1).
На втором шаге обращаем "остаток" и вновь выполняем деление. В результате выделяем емкостный элемент С1 = а0/Ь0, тогда можем записать
2 М = ^ + >
у уш С! + 1/ 22
где
Z2(j Ю) =
n-1 , n-2 ,
а 0ю + а^ + ... + а n-i yXp0®n-1 + pi®"-2 +^ + P„-i)'
(2)
Рк = рк+1 — С$к+1 — коэффициенты, определяемые по схеме, аналогичной ак.
Выражение для Z2 имеет тот же вид, что и (1), но его порядок на единицу ниже. В регулярном случае, когда все коэффициенты с нулевыми индексами не равны нулю, разложение можно довести до конца и представить 2 (у ю) в форме цепной дроби:
2 (ую) = -у'Х 1 + 1
+
j ю Ci + -
1
(3)
-jX 2 +■
1
j ЮС2 + ... + ll(-jXn + 1/j юС„)
Выражению (3) соответствует лестничная цепь (рис. 1а), содержащая псевдореактансы Xjj и емкостные элементы Ck. Для пассивной (реактивной) цепи знаки X могут быть любыми, но все величины Ck, к = 1,2,..., п, должны быть положительные. Согласно доказанной в [15] теореме разложение (3) справедливо и для квазиреактивных цепей с тем лишь отличием, что теперь часть элементов Ck становится отрицательной и их число определяет число правых корней характеристического полинома. Сформулируем полученный результат как теорему типа Рауса—Кауэра.
Теорема 2. Комплексный полином H(p) порядка п асимптотически устойчив тогда и только тогда, когда все коэффициенты Ck (к = 1,2,..., n) в разложении ассоциированного с Н(ую) иммитанса W(ую) = = Нв(ю)/jHH(®) в цепную дробь вида (3) положительные. Число правых корней равно числу отрицательных коэффициентов Ck. Знаки псевдореактансов X на устойчивость цепи не влияют.
Теорема верна для регулярного случая, когда все элементы с нулевым индексом (a0,b0, а0,..) отличны от нуля и число коэффициентов Ck равно порядку полинома. Нарушение этих условий принято классифицировать как критические или особые случаи 1-го и 2-го типов [1—4].
В особом случае типа 1, как и в схеме Рауса для вещественного полинома, нулевой коэффициент (например, у0) можно заменить малой величиной s и продолжить разложение. Однако такой подход может не дать правильного решения, если оба усло-
вия регулярности нарушаются одновременно [1]. Мы предлагаем новый более радикальный метод, применимый как для вещественных, так и комплексных полиномов. Его идея заключается в смене формы разложения импеданса 2 (у ю) в цепную дробь, начиная со строки, предшествующей появлению нулевого коэффициента, и возможно на более ранней стадии. Для вещественного полинома это означает переход от 1-й формы Кауэ-ра ко 2-й, а для комплексного полинома — от формы 1А к одной из оставшихся (подходящих) форм типа Кауэра (см. рис. 1). Эта вариативность метода Кауэра становится очевидной, если учесть, что "остаточный" иммитанс на каждом шаге реализации является реактивной или квазиреактивной цепью, и для ее реализации, кроме критического, остается несколько (возможно, один) реализуемых вариантов. Обсудим этот прием в Приложении; заметим, что он применим и для табличной формулировки критерия.
2. ТАБЛИЧНЫЙ КРИТЕРИЙ ТИПА РАУСА
Оформим вычисление коэффициентов X,, Ск при разложении импеданса 2 (у ю) в цепную дробь (3) в виде табл. 1. Назовем ее таблицей типа Рауса. В первые две строки таблицы входят коэффициенты вещественной и мнимой частей полинома Н (у ю), остальные заполняются по схеме Рауса для вещественных полиномов — по правилу перекрестного перемножения. В дополнительном левом столбце рассчитываются элементы X, Ск цепи, синтезированной по форме 1А. Признаком физической реализуемости цепи (т.е. у
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.