![](/user_photo/_userpic.png)
Электроника и схемотехника, учебное пособие, Щ.А.С., О.А.Г
..pdf![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs391x1.jpg)
391
чтобы выбранный функционал принимал минимальное значение. Здесь - множество состояний предельного цикла периодического решения уравнения (8.127).
Поставленную задачу оптимального управления для вырожденных автоколебательных систем можно будет считать решенной, если будут определены условия существования, период и число интервалов оптимального предельного цикла, определена его область притяжения, а также дано решение задачи синтеза, т. е. определено управление в функции
положения фазовой точки g(t) f (u) .
Если оптимальный предельный цикл существует в вырожденной автоколебательной системе, то он имеет два интервала знакопостоянства управления независимо от структуры системы.. Для выделения этих двух интервалов введем дискретную перемев-
шую s , принимающую значение «1» во время подпитки РГН от источника питания, и значение «0» -в.режиме диссипации электрической энергии.
Выбор оптимального закона управления во время подпитки РГН от источника питания сделать нетрудно. Любая система стремится к минимуму потребления мощности на управление и к минимальной разности между заданным режимом и действительным значением фазовой координаты. Учитывая, что наша автоколебательная система вырожденная, запишем искомый функционал в скалярной форме
|
|
1 |
t1 |
|
|
J |
|
f u 2 (t ) [q u 2 (t) r g 2 (t)]dt |
|
|
|
|
|
, |
(8.128) |
||
|
2 |
||||
|
1 |
1 |
|||
|
|
|
|
||
|
|
|
0 |
|
|
где f , q и r - в общем случае некоторые функции времени, но для упрощения будем их считать постоянными.
Как уже указывалось ранее, работа коммутирующего устройства, управляющего процессом заряда конденсатора, начинается в момент времени, равный нулю, и заканчивается
в момент времени t 1 . Этим обусловлен выбор пределов интегрирования в формуле
(8.128).
Найдя решение уравнения (8.128), определим оптимальную форму переходного процесса во время заряда конденсатора. Определение основных параметров предельного цикла (частоты, амплитуды, формы сигнала и т.д.) требует нахождения минимума функционала вида
|
T |
|
|
J 2 |
dt T |
, |
(8.129) |
|
0 |
||
|
|
|
где Т – период автоколебаний РГН.
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs392x1.jpg)
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs393x1.jpg)
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs394x1.jpg)
394
Определим параметры оптимального предельного цикла [21]. Для чего воспользуемся методом припасовывания (точечных отображений) L3IJ. Определим переходные процессы заряда-разряда конденсатора С (рис. 8.23), если рассматривать РГН как устройство, оптимальное по быстродействию, установившимся режимом которого будет предельный цикл вырожденной автоколебательной системы [21].
Для нахождения параметров оптимального по быстродействие переходного процесса воспользуемся принципом максимума Л.С.Понтрягина [18]. Используя дискретную переменную s, можно записать следующую систему уравнений:
du |
s |
u 1 |
u 2 |
ch 2 ( |
t1 ),0 t t1 , |
|
u 2 |
||||||
dt |
|
|
|
(8.138) |
||
|
|
0 |
|
|||
|
|
|
|
|
||
|
|
(1 s) ( |
1 u),t1 |
t T , |
где s – дискретная переменная.
Найдѐм минимум функционала (8.129) для вырожденной автоколебательной системы, описываемой, описываемой системой уравнений (8.138). Решением поставленной задачи будет закон изменения управляющего воздействия s во времени, который бы удовлетворяло ограничениям, накладываемым на параметры предельного цикла, и обеспечивал пе-
ревод системы (8.138) в область за минимальное время.
Если оптимальное управление s существует, то этому управлению должен соответствовать вектор сопряжѐнных переменных. Гамильтониан для рассматриваемой системы (8.138) представим в виде
H |
1 p1 |
u(t) |
, |
(8.139) |
|
|
|
где p1 - сопряжѐнная переменная, которую можно определить, решив дифференциальное уравнение
p1 |
H |
|
|
u . |
(8.140) |
||
|
После соответствующих подстановок и интегрирования уравнения (8.140), получим
|
|
|
|
ch2 [ (t |
t )] |
sh[ |
(t t )] |
|
|
||
p1 |
p10 |
exp[ 1 |
(1 s) t] |
|
1 |
|
1 |
|
|
|
|
ch2 ( |
t1 ) |
sh( |
t1 ) |
, |
(8.141) |
||||||
|
|
|
|
где p10 const .
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs395x1.jpg)
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs396x1.jpg)
396
ными интегральными формулами. Далее будем рассматривать функцию модуля, причѐм будем пользоваться квадратом модуля, а не самим модулем, так как это проще. Очевидно, если функция модуля дана, еѐ можно возвести в квадрат. Мы полагаем, что квадрат модуля задан в виде рациональной функции от s, а, не в виде кривой, для которой нужно найти приближенное выражение.
Допустим теперь, что дана рациональная функция от ω, скажем, G( j ) , удовлетворяющая необходимым условиям, при которых она является квадратом модуля функции цепи. Эти условия таковы: функция должна быть чѐтной функцией от ω, а степень числителя еѐ не выше степени знаменателя больше, чем на два, так как передаточная функция не может иметь больше одного простого полюса в бесконечности. Кроме того, если
G( j ) имеет конечные полюсы на оси j, они должны быть двукратными, так как полюсы
функции цепи на оси j |
должны быть простыми. Данную функцию можно представить |
|
||||||||
как |
|
|
|
|
|
|
|
|
|
|
|
G( j |
) |
|
F( j ) |
|
2 |
F( j ) F( |
j ) . |
(8.147) |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||
Необходимо определить функцию F ( j |
) , если известна передаточная функция G( j |
) |
||||||||
По существу можно рассматривать G как функцию комплексного переменного s, |
а G( j |
) |
||||||||
- как ее значение на оси |
j |
. Таким образом, можно написать соотношение |
|
|
||||||
|
|
|
|
G ( s ) = F ( s ) F ( - s ) , |
(8.148) |
|
||||
|
|
|
|
|
|
|
|
|
||
которое справедливо для всех значений s, а не только при s j . Однако G( j ) |
не равна |
|||||||||
квадрату модуля F(s), за исключением случая, когда s |
j . |
|
|
Будем разлагать числитель и знаменатель функции G ( s ) на множители. Некоторые из них нужно включить в F(s), а другие в F(-s). Сначала найдем, какое будет расположе-
ние нулей и полюсов функции G(s). Напомним, что G ( s ) есть отношение, двух четных многочленов.
Теперь возникает вопрос, как выбрать полюсы и нули функции F(s) из полюсов и ну-
лей функции G(s). Для полюсов ответ простой. Мы знаем, что F(s) должна быть регулярна в правой полуплоскости. Поскольку полюсы функции F(-s) равны полюсам функции F(s) с противоположным знаком, то функция F(-s) регулярна в левой полуплоскости. Следовательно, распределение полюсов функции G(s ) может быть однозначным : п о -
люсы G(s ) в левой полуплоскости
луплоскости принадлежат F(-s). Полюсы на оси j распределены равномерно. Они явля-
ются двукратными в функции G ( s ) и простыми в функциях F(s), и F(-s). Распределение
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs397x1.jpg)
397
нулей приведено на рис. 8.13. Такая симметрия распределения нулей называется квадратной симметрией. Говорят, что полюсы и нули передаточной функции G(s ) или квадра-
та модуля функции распределены симметрично по квадратам.
Что касается выбора нулей, то он не столь прост. Если F(s) не является входной функцией, то условие, что она не должна иметь нулей в правой полуплоскости, не обязательно. Поэтому не обязательно приписывать функции F(s) все нули функции
G ( s ) , расположенные в левой полуплоскости. Но если задано, что F(s) должна быть функцией минимальной фазы, то все нули функции G(s), расположенные в левой полу-
плоскости, следует приписать F(s). Если это условие не задано, то функции F(s) можно приписать некоторые нули правой полуплоскости. Однако надлежит помнить, что пара комплексных нулей должна быть объединена, так как F(s) должна быть действительной функцией от s. Итак, F(s) однозначно определяется функцией G(s) только в том случае, если F(s) должна быть функцией минимальной фазы.
Наконец, нули на оси j следует разделить поровну между F(s) и F(-s), так же как полюсы.
Разберем теперь для пояснения этого метода несколько примеров; имеющих практическое значение. Для обычных электрических фильтров требуются передаточные функции, имеющие совершенно постоянные модули на оси j в данном интервале частот, который называется полосой пропускания, и в точности равные нулю на остальной оси j , в так называемой полосе задерживания. Модуль рациональной функции на оси j не может иметь такую идеальную форму. Однако можно найти передаточные функции, модули которых на оси j так или иначе приближаются, к желаемой функции.
' Рассмотрим передаточную функцию идеального фильтра нижних частот, изображенную на рис. 8.24, а. На рис. 8.24, б, в показаны два способа приближенной реализации этой идеальной функции. Первый способ называется максимально гладкой аппроксимаци-
ей или аппроксимацией Баттерворта, а второй; способ - равноволновой аппроксимацией
или аппроксимацией. Чебышева. Максимально гладкая аппроксимация монотонна в полосе пропускания и в полосе задерживания, причем максимальная ошибка получается вблизи границы полосы пропускания. Аппроксимация Чебышѐва является осциллирующей в полосе пропускания, причѐм пики волн одинаковы. Благодаря этому ошибка распределена в полосе пропускания более равномерно.
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs398x1.jpg)
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs399x1.jpg)
399
тель. В данном случае это выполняется сравнительно просто. Нули знаменателей находятся из выражения
|
s 2n exp[( 2k |
1 |
n) |
] |
|
1 |
, |
|
|
|
|
|
(8.154) |
|
|
|
|
|
|
|
|
|
|
|
|
||
где знак минус соответствует чѐтному n. |
Извлекая корень степени 2n |
из обеих частей |
|||||||||||
уравнения (8.154), найдѐм полюсы функции G(s): |
|
|
|
|
|
|
|
|
|
|
|||
sk |
exp[ j(2k 1 |
n) |
/(2 |
n)] |
|
|
|
|
|
|
|
|
|
; |
|
k 1,2n |
. |
|
|
(8.155) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||
Передаточная функция имеет 2n полюсов и модуль каждого из них равен единице. |
|||||||||||||
Полюсы расположены равномерно по окружности единичного радиуса. При n |
|
4 получа- |
|||||||||||
ем расположение корней: |
|
|
|
|
|
|
|
|
|
|
|
|
|
s1 exp[ j(5 / 8)]; s2 |
exp[ j(7 / 8)]; s3 |
exp[ j(9 |
|
/ 8)]; s4 exp[ j(11 |
/ 8)] |
. |
(8.156) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
Передаточная функция F(s) в этом случае принимает следующий вид:
F (s) |
1 |
1 |
|
||
|
|
|
|
||
(s s1 )(s s2 )(s s3 )(s s4 ) |
1 2,613 s 3,414 s 2 2,613 s3 s 4 . (8.157) |
||||
|
Найденная функция F(s) имеет вид, представленный на рис. 8.24, б.
Рассмотрим реакцию, вызванную аппроксимацией Чебышѐва (8.150). Чтобы найти полюсов, заменим j на s и приравняем знаменатель нулю. Получим
Tn (s / j) |
ch[nArch(s / j)] |
j / |
. |
|
(8.158) |
|
|
|
|
|
|
||
Используем преобразование w |
Arch(s / j) |
x |
jy . Тогда получим |
|
||
s jch(w) |
jch(w) |
jch(x |
jy) , |
|
|
(8.159, а) |
Tn (s / j) ch(n w) ch[n (x jy )] |
j / |
. |
(8.159, б) |
|||
|
|
|
|
|
Если в последнем уравнении раскрыть ch(n·w) и приравнять вещественные и мнимые части этого выражения нулю, то получатся значения х и у, удовлетворяющие этому урав-
нению. Если их обозначить как sk |
k |
j |
k , то в результате получим |
|
|
|
||||||||||||
|
sh |
1 |
Arsh |
1 |
sin |
|
(2k 1) |
|
|
ch |
1 |
Arsh |
1 |
cos |
|
(2k |
1) |
|
k |
n |
|
|
2 n |
, |
k |
n |
|
|
2 n |
. (8.160) |
|||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Последние выражения можно объединить, произведя обычные операции
![](/html/93263/27/html_TzgGJ_mAH0.txFo/htmlconvd-JfEvCs400x1.jpg)