Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Электроника и схемотехника, учебное пособие, Щ.А.С., О.А.Г

..pdf
Скачиваний:
18
Добавлен:
13.10.2024
Размер:
17.78 Mб
Скачать

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

 

 

 

где Т – период автоколебаний РГН.

392

Задачу оптимального быстродействия сформулируем следующим образом. В фа-

зовом пространстве U системы (8.127) даны две точки u0 и u1 . Среди всех допустимых

управлений g(t) , переводящих фазовую точку из положения u0 в положение u1 , найти такое, для которого функционал (8.129) принимает наименьшее возможное значение. Эта задача с закрепленными концами обращается в задачу с подвижными концами, если требуется оптимальным образом перевести фазовую точку из заранее незаданного положения

u0 U 0 в некоторое положение u1 U1 .

Для решения задачи синтеза вырожденной автоколебательной системы необходимо оптимальное управление определить в функции координат фазовой точки g(t) f (u) .

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

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

ния (8.127) и на существовании в фазовом пространстве области их притяжения Z. При этом вопросы устойчивости найденных решений - оптимальных предельных циклов и оп-

ределения границы области притяжения Z остались открытыми/ В то же время эти вопросы являются важными, так как определяют возможности физической реализации оптимальных управлений. Устойчивость периодических решений в смысле А.М.Ляпунова [17] определяется поведением фазовых траекторий вблизи предельного цикла [21].

Результат решения системы уравнений (8.127) и (8.128) можно представить в такой форме

u(t)

 

k (t)

u(t)

 

 

1

r

,

(8.130)

 

 

 

 

 

 

 

где скаляр k (t) - есть решение уравнения Риккати [17]:

k(t) 2 1 k(t)

k 2 (t)

q

 

r

(8.131)

 

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

393

с граничным условием

k(t1 )

f

.

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрирование выражения (8.131) даѐт

 

 

 

 

 

(

 

 

)

 

( f / r

 

1

)

exp[2

(t t1 )]

1

1

 

 

 

 

f / r

 

 

 

 

k(t)

 

 

 

 

 

 

 

 

 

1

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f / r

1

 

 

 

 

 

 

 

1

 

 

 

exp[2

(t

t1 )]

 

 

 

 

 

 

 

 

 

 

f / r

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где параметр

 

определяется выражением

 

 

 

 

 

 

q

2 .

 

 

r

1

 

 

 

(8.132)

(8.133)

Задавая соотношение между неопределѐнными коэффициентами в функционале

(8.128) в виде

 

f / r

1

1

,

 

 

 

 

 

 

 

f / r

1

 

 

 

 

 

 

(8.134)

 

 

 

 

формулу (8.130) можно переписать

 

 

 

 

u(t)

th[ (t t1 )] u(t) .

(8.135)

 

 

 

 

Принимая начальное значение напряжения на конденсаторе С (рис. 8.23) равным u0 , решение дифференциального уравнения (8.135) можно записать как

u(t) u

 

ch(

t1 )

 

 

 

0 ch[

(t t

1

)] .

(8.136)

 

 

 

 

 

 

 

Решая уравнение (8.136), а затем подставляя результат в формулу (8.135), получим

u(t)

1

u 2

ch

2 ( t )

 

 

 

 

 

 

 

 

u02

 

1

.

(8.137)

Последнее выражение зависят от параметра

 

. Этот параметр определяет динамиче-

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 .

395

Учитывая, что s может принимать дискретные значения (0 или 1), а значение 1 обыч-

но мало, то для моментов времени t 0 и t t1 при изменении s значение p1 остаѐтся практически постоянным.

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

 

 

 

 

 

H

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.142)

 

 

 

 

 

s

 

 

 

Момент смены состояний РГН должен происходить при

 

 

 

 

t

 

 

1

Arth

 

1

.

 

 

 

(8.143)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этот момент напряжение на конденсаторе С (рис. 8.23) достигает значения

 

u1

 

u(t1 )

 

 

 

u0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2 .

(8.144)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

Учитывая, что в режиме диссипации разряд конденсатора С (рис. 8.23) (рис. 3.1) будет происходить через резистор R , то нетрудно получить выражение, связывающее период автоколебаний с определенными ранее параметрами автоколебания РГН

u1

u0

exp[ 1

(T t1 )].

 

 

 

(8.145)

Совместное решение последних уравнений (8.144) и (8.145) позволяет определить период автоколебаний релаксационного генератора

T

1

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

2

 

 

 

 

 

1

 

 

 

 

 

 

 

1 Arth 1 .

(8.146)

Форма импульсов на конденсаторе С (рис. 8.23) , состоящая из двух экспонент (заряд и разряд конденсатора С), на экспериментальном РГН практически полностью совпадает с вычисленными кривыми.

8.7. Определение функции цепи по заданному модулю

Рассмотрим обратную операцию: построение рациональной по одной еѐ вещественной или мнимой части или по еѐ модулю или углу. Из предыдущего материала известно, что вещественная и мнимая части аналитических функций связаны между собой извест-

принадлежат функции F(s), а полюсы правой по-

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). Распределение

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, б, в показаны два способа приближенной реализации этой идеальной функции. Первый способ называется максимально гладкой аппроксимаци-

ей или аппроксимацией Баттерворта, а второй; способ - равноволновой аппроксимацией

или аппроксимацией. Чебышева. Максимально гладкая аппроксимация монотонна в полосе пропускания и в полосе задерживания, причем максимальная ошибка получается вблизи границы полосы пропускания. Аппроксимация Чебышѐва является осциллирующей в полосе пропускания, причѐм пики волн одинаковы. Благодаря этому ошибка распределена в полосе пропускания более равномерно.

398

Рис. 8.24

Аналитические выражения этих функций с точностью до постоянного множителя имеют вид:

F ( j )

F ( j )

2

 

1

 

 

 

 

 

 

 

 

- аппроксимация Баттерворта;

(8.149)

1

2n

2

 

1

 

 

 

 

 

 

 

 

 

 

1

2T 2

( ) - аппроксимация Чебышѐва;

(8.150)

 

 

n

 

 

 

где - малое число, определяющее амплитуду волн, а ω = 1 соответствует границе поло-

сы пропускания. Функция

Tn (

)

представляет

полином Чебышева.-,

 

 

 

 

 

 

 

 

 

 

 

 

Tn (s / j)

ch[n

Arcch(s / j)]

,

 

 

(8.151)

 

 

 

 

 

 

 

 

 

 

который приводится после подстановки s

j

к виду

 

 

 

 

 

 

Tn ( ) cos(n arccos )

при

 

 

1

;

(8.152)

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдѐм передаточную функцию F(s), если известен квадрат ее модуля на оси j. Рассмотрим сначала реакцию фильтра Баттерворта. Сначала в (8.149) осуществим следую-

2

s 2

щую замену:

и получим:

G(s) F (s)F ( s)

1

.

.

 

1 ( 1)n s 2n

(8.153)

Эта функция не имеет конечных нулей, поэтому требуется разложить лишь знамена-

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)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последние выражения можно объединить, произведя обычные операции

400

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

k

 

 

 

 

 

 

k

 

 

1.

 

sh

2 1

Arsh

1

 

 

ch

2 1

Arsh

1

 

(8.161)

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это уравнение эллипса в плоскости s. Горизонтальная ось эллипса совпадает с осью

j

и будет короче вертикальной (ζ), так как гиперболический косинус вещественного

 

переменного всегда больше гиперболического синуса.

Рассмотрим реакцию фильтра Чебышѐва, когда заданы параметры: 0,1; n=4. Передаточная функция F(s) в этом случае принимает следующий вид:

F (s)

1

1

 

 

 

 

 

(s 2 0,641s 1,534 )(s 2 1,519 s 0,823)

 

s 4 2,16s3 3,31s 2 2,86s 1,26 . (8.162)

 

 

Найденная функция F(s) имеет вид, представленный на рис. 8.24, в.

9. ЦИФРОВЫЕ ФИЛЬТРЫ

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

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

Первый расчѐтный алгоритм основывается на использовании свѐртки двух функций, одна из которых является передаточной функцией, а вторая – входным сигналом. Рассмотрим сигнал х(пТ). Для простоты рассмотрения ограничим х(пТ) условием, что он равен нулю при n 0, так что входной сигнал состоит из последовательности чисел х(0), х(Т), х(2Т) и т. д. Постоянная Т является интервалом дискретизации, поэтому при таком обозначении неявно подразумевается, что х(пТ) может быть получено путем взятия отсчетов у непрерывного сигнала x(t) через равные промежутки времени Т, т. е. в моменты времени 0, Т, 2Т и т. д.

Линейный цифровой фильтр можно определить с помощью принципа суперпозиции следующим образом:

 

n

n

y(nT )

h(mT ) x(nT mT )

h(nT mT ) x(mT ) ,

m 0

m 0