Оптимальные процессы во времени приведены на рис. 8.25,6 и в. Графики •/(т) характеризуют закон изменения вращающего момента (или напряжения) двигателя.
Из рассмотренного следует, что для получения минималь ного времени отработки начального отклонения в системе не обходимо сначала осуществить разгон х с максимально допу стимым положительным ускорением (максимальным положи тельным моментом), а затем при максимальном отрицательном по знаку ускорении (максимальном отрицательном моменте) произвести торможение этой величины.
Указанный способ обеспечения минимального времени регу лирования в системе является следствием более общего прин ципа максимума Понтрягина Л. С., из которого вытекает [1], что для получения максимального быстродействия при любом порядке линейного дифференциального уравнения объекта все гда необходимо использовать релейное управление.
Если порядок уравнения объекта равен двум, то минималь ное время регулирования получается при предельном отклоне нии регулирующего органа в одну сторону и переключении его затем в противоположном направлении также на максимально возможную величину (одно переключение реле). В более об щем случае, если порядок уравнения объекта равен п, то для получения максимального быстродействия (при отрицательных вещественных корнях характеристического уравнения) за вре мя переходного процесса регулирующий орган должен п—1 раз
|
|
|
|
|
|
|
|
менять |
знак |
отклонения. |
Последнее |
положение |
доказано |
Фельдбаумом А. А. [5]. |
что |
оптимальный процесс при С |
Из рис. 8.25,п |
видно, |
- const |
может |
быть только |
при какой-либо одной |
группе на |
чальных |
условий |
(например, |
когда все |
начальные |
положения |
изображающей точки лежат па отрезке траектории хоцщ). При начальном положении на д'озиіз и на х^пц процессы уже не будут оптимальными. В первом случае процессы будут с пере регулированием, во втором — будет скользящий режим.
Для того чтобы при любых начальных отклонениях переход ный процесс был оптимальным, линией переключения должна быть сама фазовая траектория отп\р (и соответственно оп'т\'р' — в противоположном квадранте).
Непрямая линия переключения означает нелинейный харак тер воздействия корректирующего устройства.
Уравнение системы стабилизации углового положения кос мического летательного аппарата с нелинейной корректирую щей цепью будет иметь вид:
Зс + F[x + /(х )] = 0. |
(8.35) |
Функция f нелинейного корректирующего сигнала зависит от скорости X и свойств реле. Следовательно, для получения опти
мального процесса релейная система должна снабжаться спе циальным вычислительным устройством, образующим функ цию f. Для уравнения (8.35) линией переключения ріщпо будет парабола
X + — з/2 = 0.
2
Следовательно, искомая функция |
|
f[x) = |
(jr)2signA\ |
(3.36^ |
Множитель sign.v указывает, что квадратичная функция про изводной должна быть нечетной функцией (для образования линии переключения р'т\'п'о).
Для релейной следящей системы нелинейная функция име ет еще более сложный вид. Однако фазовая траектория, на ко торой должно происходить переключение, не сильно отличается от прямой. Мирясь с некоторым отступлением процессов от оп
тимальных, линию переключения можно |
взять в виде прямой, |
и оптимальное значение порции сигнала |
от тахогенератора бу |
дет при С= 0,3 [1]. |
|
§ 8.4. МЕТОД ГАРМОНИЧЕСКОЙ ЛИНЕАРИЗАЦИИ
Метод фазовой плоскости дает достаточно точную картину движения нелинейной системы, однако он применим лишь для систем второго порядка. Для исследования автоколебаний в системах, линейная часть которых описывается уравнениями более высоких порядков, удобно использовать метод гармониче
ской линеаризации, разработанный |
академиками Крыловым |
Н. М. и Боголюбовым Н. Н. |
Л. С., распространившим |
Метод развит далее Гольдфарбом |
его на системы автоматического регулирования для исследо вания устойчивости и автоколебаний нелинейных систем, и По повым Е. П., осуществившим его дальнейшее обобщение и раз витие (в частности для исследования качества регулирования нелинейных систем). Метод гармонической линеаризации по зволяет исследовать устойчивость систем и исследовать уста новившиеся процессы, т. е. определить амплитуду и частоту автоколебаний и получить зависимость основных параметров автоколебаний от параметров системы. Метод является прибли женным, но довольно широко применяется, так как по своему математическому аппарату является весьма простым и, кроме того, обеспечивает достаточно высокую для практики точность. Чтобы пояснить сущность метода гармонической линеаризации, представим нелинейную систему в виде двух частей: нелиней ной и линейной (рис. 8.26,а).
Применяя метод гармонической линеаризации, мы полага ем, что автоколебания в системе существуют в синусоидальной форме
X = А sin t,
где А и и) — искомые амплитуда и частота автоколебаний. Эти автоколебания, действуя на нелинейную часть, создадут на вы ходе периодический несинусоидальный сигнал, действующий на
Нелинейная У, |
Линейная X |
часть |
часть |
о) |
|
б)
Рис . 8.26. К пояснению метода гармонической линеаризации:
а — представление нелинейной системы в виде нелинейной и линейной частей; б —
структурная схема нелинейной системы с од ним нелинейным элементом
входе линейной части системы. Однако линейная часть систе мы, как правило, представляет собой фильтр низких частот, хорошо пропускающий низкие частоты. Вследствие этого перио дический сигнал у на выходе линейной части системы превра щается в сигнал X, близкий к синусоидальному.
Подавление линейной частью высокочастотных гармоник по зволяет при анализе автоколебаний в нелинейной системе огра ничиться в первом приближении учетом лишь первой гармони ки периодической несинусоидальной кривой у.
Этот метод позволяет получить тем более точные результа ты, чем ближе линейная часть по своим «подавляющим свой ствам» к фильтру низких частот.
Структурную схему нелинейной системы с одним нелиней ным элементом представим в виде рис. 8.26,6.
Запишем уравнение линейной части системы |
|
X — W (D) у = Q ( £ > ) |
ba + b^DA-------b bmDm |
У, (8.37) |
|
|
P {D) У = |
c0 + ^ ö - j ----- cnDn |
где D |
d |
символ дифференцирования. |
|
|
|
dt
Нелинейная часть системы может быть представлена нели нейной функцией
|
|
|
y = F [ - x ) |
= |
- F{x). |
|
(8.38) |
Знак |
минус |
обусловлен |
учетом |
уравнения |
рассогласования |
е = л;вх — X |
= — X , |
так как |
входной сигнал |
|
полагаем рав |
ным нулю (хвх — 0). |
|
сигнал |
|
л;= уі sin оі ^ |
с |
неизвестными |
Пусть гармонический |
ш |
нам амплитудой А п частотой |
автоколебаний |
действует |
на |
вход |
нелинейного |
элемента. В |
этом случае его |
выходная |
ве |
личина будет периодической функцией времени |
|
|
|
|
|
у = — F(x) = — F(i4 sin <о£). |
|
(8.39) |
Полагая, что F (х) — однозначная нечетная функция х, разло жим ее в ряд Фурье и ограничимся лишь первым членом ряда
|
у = |
— F [А sin u>f) == |
— i]) (А) sin id*;, |
(8.40) |
где |
4>(Л) = — f |
F (А sin ш t) s i n |
d[y>t) |
— коэффициент ря- |
|
J . |
|
|
|
|
да |
Фурье или амплитуда первой |
гармоники выходной |
величи |
ны нелинейного элемента. ф(/1) зависит от вида характеристи ки нелинейного элемента и амплитуды сигнала А, действующе го на его входе.
Рассмотрим некоторые типы нелинейных характеристик и со ответствующие нм зависимости ф(4).
Если на вход идеального нелинейного элемента (рис. 8.27,а) действует гармонический сигнал а' = і4 зіпшЛ то на его выходе
^ y ip(f) ^ |
F(Asincöt) |
I |
\ |
I |
L |
|
и |
ъ |
|
І |
” |
S) |
|
ФЩ |
|
|
|
|
|
■4В. |
О) |
0 |
|
1К |
|
А |
|
|
Ь) |
Рис . 8.27. Идеальный релейный элемент: |
а — характеристика и |
зависимость |
х = А sin coif; |
б — периодический несинусоидальный сигнал ре
лейного элемента |
F{A sin о> t) |
и закон |
изменения |
амплитуды |
первой |
гармоники |
ф(.Л) во |
времени; |
в |
— зависимость ф(4) |
|
получим 'Последовательность |
знакопеременных |
импульсов |
(рис. 8.27,6). Амплитуда первой гармоники такой |
периодиче |
ской функции постоянна и не |
зависит от А (рис. |
8.27,в). Она |
равна: |
|
|
|
45 |
АI ) |
ФИ) = |
(8 |
ТС |
|
Для нелинейного элемента с зоной нечувствительности (рис. 8.28,а) выходной сигнал представляет собой серию прямоуголь ных импульсов, следующих с паузой друг за другом (рис. 8.28,6). Амплитуда первой гармоники в этом случае зависит от амплитуды входного сигнала:
45 |
(8.42) |
ф(А) = — V А 2- я2- |
itА |
|
Характеристика ф(А) приведена на рис. 8.28,в.
/Ф(А |
|
/ F(AsLTi(Ot) |
1 Л |
/ |
/ ' |
' |
" Т |
|
|
|
|
Р и с. 8.28. Релейный элемент с зоной нечувствитель ности:
а — характеристика н зависимость х = А sin ш і\ б — периодический несииусоидальный сигнал релейного элемента F (А sin ші) и закон изменения амплиту ды первой гармоники ф (Л> во времени; в — зави
симость ф
Для получения основных зависимостей, с помощью которых могут быть определены А и ш, вводится понятие о передаточ ной функции нелинейного элемента.
Если входной сигнал нелинейного элемента запишем в комплексной форме
x = A e Jmt, |
(8.43) |
26* |
403 |
а амплитуду первой гармоники нелинейного элемента предста вим как
3> = ф ( Л ) е Уо,<, |
( 8 . 4 4 'і |
то отношение указанных комплексных величин и будет назы ваться передаточной функцией нелинейного элемента
|
V |
Ф(Л) |
|
|
|
|
= |
. |
( 8 . 4 5 ) |
или |
X |
А |
|
|
|
|
|
|
|
У = |
Л |
|
( 8 . 4 6 ) |
|
|
|
|
Понятие |
передаточной функции |
^ „(4 ) |
позволяет |
заменить |
данное нелинейное звено для |
периодического режима и по |
стоянной |
амплитуды А некото*рым линейным эвеном |
с коэффи |
циентом усиления W H(A), так как нелинейное уравнение (8.33) заменяется линейным (8.46). Такой процесс сведения нелиней ной зависимости к линейной носит название гармонической ли неаризации нелинейных зависимостей. Это название вытекает из разложения нелинейных колебаний на гармонические состав ляющие (гармоники) и пренебрежения высшими гармониками на выходе нелинейного элемента и линейной части системы. В результате гармонической линеаризации мы получаем своеоб разное линейное звено, коэффициент усиления которого ІУДА) зависит от амплитуды входного сигнала А.
Передаточные функции нелинейных звеньев с однозначными характеристиками (см. рис. 8.2,а, б, в) представляют собой дей
ствительные числа, так как они не |
создают фазовых сдвигов |
при прохождении гармонических сигналов и ф(А) |
для них чис |
ло вещественное. |
|
характеристиками |
(типа |
|
Для звеньев с неоднозначными |
рис. 8.3,а, б, в) гармонический выходной сигнал |
нелинейного |
элемента (первая гармоника) имеет |
сдвиг |
фазы по отношению |
к сигналу на входе |
|
|
|
|
|
|
у=Ъ(А) еЛш,- 9,Д)1 = |
ф(А)еМ/|) |
= ф (/А) e/W , |
(8.47) |
где |
ф(/А) — комплексная |
амплитуда первой гармоники; |
|
ср (А) — сдвиг фаз между выходным и входным |
сигна |
|
лами. |
|
|
|
|
|
та |
Следовательно, передаточная функция нелинейного элемен |
с неоднозначной характеристикой |
есть |
число |
кташлекеное |
|
|
Фи А) |
|
|
(8.48) |
|
W H( j A ) = Z - |
А |
|
|
|
|
X |
|
|
|
Как и для нелинейных элементов с однозначными |
характери |
стиками, выражение (8.48) |
можно |
рассматривать |
как |
комп |
лексный коэффициент усиления некоторой линейной системы и, следовательно, в данном случае имеет место гармоническая линеаризация нелинейной зависимости.
Гармоническая линеаризация нелинейных элементов дает возможность исследовать устойчивость нелинейных систем, вы явить автоколебательные режимы, оценить их устойчивость.
Понятие о передаточной функции нелинейного элемента по зволяет представить структурную схему в виде последователь ного соединения W H(jA) и W (jus) (рис. 8.29).
Р и с. 8.29. К пояснению метода гармо нической. линеаризации. Представление нелинейной системы в виде последова тельного соединения Wn{jA) и W (/шj
По правилу встречно-параллельного соединения звеньев пе редаточная функция для замкнутой системы есть выражение
Ф ( і А М) |
W » U A ) W U * ) |
W ( j A , * ) |
(8.49) |
|
\ + W B ( j A ) W ( M |
. 1 + IF ( / Л , со) ’ |
а характеристическое уравнение замкнутой системы имеет вид:
|
I + WUA, Ш) = 0, |
|
где |
W {j А, со) = W H(JA) W{Ja>) — передаточная |
функция ра |
|
зомкнутого контура нелиней |
В |
ной системы. |
амплитуды |
общем случае W (j A , со) является функцией |
входных колебаний и частоты. При постоянном значении А и изменении со от 0 до оо вектор W (/А , ш) вычертит на комплекс ной плоскости годограф, аналогичный амплитудно-фазовой ха рактеристике разомкнутой системы W (/со).
Семейство годографов W{jA, со) |
некоторой системы для раз |
личных значений А |
приведено на |
рис. 8.30. |
Если |
годограф |
W{jA,<ü) при .всех А |
не охватывает |
точку — |
1,/0, |
то система |
устойчива, и автоколебаний в ней не возникает. В этом случае
при любых А модуль вектора |
| W{jA, < o j | < 1 . |
Здесь |
— ча |
стота, |
при |
которой со (со, Л) = |
— іг. Физически |
это |
будет |
озна |
чать, |
что |
если в результате действия каких-либо |
возмущений |
в системе возникнут колебания, то после снятия этих возму
щений колебания будут затухать. |
| W [JA, <о_)| > 1, |
Если при некоторых значениях А на частоте |
то система неустойчива, колебания в ней |
будут нарастать. |
II, наконец, если j W (jA, <uj | = I, система находится па гра нице устойчивости, в системе будут автоколебания. Это проис ходит при определенных значениях А и ш, которые и являют ся искомыми амплитудой и частотой автоколебаний Аа и ша.
Следует заметить, что при однозначной характеристике не линейного элемента равенство t f ( с о ) г = — л может быть исполь зовано для определения возможной частоты автоколебаний, а
|
|
|
|
|
|
|
|
|
|
|
|
второе |
соотношение! W ( jA , ш_) | = 1 |
позволяет определить ам |
|
плитуду. Вследствие того, что вектор |
W (JA, со) |
имеет |
зависи |
|
|
|
мость |
не только |
от частоты, |
но |
и от |
|
|
|
амплитуды, |
может иметься |
несколько |
|
|
|
значений А и со, |
при |
которых |
годо |
|
|
|
граф |
W JA , <о) |
проходит |
через |
точ |
|
|
|
ку — 1,/0. |
|
|
|
|
|
|
|
|
|
Это значит, что в системе возмож |
|
|
|
но несколько периодических |
режимов |
|
|
|
или автоколебаний с различными ча |
|
|
|
стотами и амплитудами. Часть из этих |
|
|
|
■режимов действительно может наблю |
|
|
|
даться |
в системе, а |
другая |
часть не |
|
|
|
может. Наблюдаться могут только ус |
|
|
|
тойчивые периодические режимы |
или |
|
|
|
автоколебания. Другая часть периоди |
|
|
|
ческих |
режимов |
является |
неустойчи |
|
|
|
вой и наблюдаться в реальной систе |
|
Рис . |
8.30. Годографы |
ме не |
может. Возникает таким обра |
|
зом задача |
выделить из всех периоди |
|
W [jA, ш)для трех значений |
|
|
|
ческих |
решений |
уравнения |
движения |
системы устойчивые решения, соответствующие устойчивым пе
риодическим движениям |
или автоколебаниям |
системы. |
При |
ближенное рассмотрение этого вопроса будет |
показано |
ниже, |
на примерах |
|
|
|
|
|
|
|
сто |
Для определения амплитуды и частоты автоколебаний вме |
построения |
годографа |
W (jA, |
ш) |
обычно |
используют |
дру |
гой |
графический |
прием, |
связанный |
с |
решением уравнения |
|
|
1 + |
W„ (JA) W (у'ш) |
= 0. |
|
(8.5Ü) |
Это уравнение можно представить в виде |
|
|
|
|
— Ен JA) = W (jui), |
|
(8.51) |
где |
E„(jA)==---- ------ |
—обратная |
передаточная функция не- |
|
WJJA) |
|
линейного элемента. |
|
|
Если на комплексной плоскости построить амплитудно-фа зовую характеристику линейной части системы W (у'ш) и об ратную отрицательную передаточную функцию нелинейного элемента — Ен(/А), то пересечение характеристик будет ука
зывать на наличие в системе автоколебаний. Пересечение ха рактеристик соответствует нахождению системы на границе ус тойчивости или прохождению годографа W ( j A , < s > ) через точку
- і,уо. |
пересечения характеристик дадут искомые значения |
Точки |
ша и Л„. |
Частота автоколебаний при этом отсчитывается по |
W (/со ), |
а амплитуда — по характеристике — Е н (/ А). |
Следует заметить, что если в системе имеются интегрирую щие звенья, то для получения области, охватываемой годогра фом W (/со), его приходится дополнять дугой бесконечно боль шого радиуса. В соответствии с определением передаточной функции нелинейного элемента его обратная передаточная функ- *
|
ция есть отношение |
А |
|
|
1 |
(8.52) |
|
W j j A ) |
У (JA ) |
|
|
Поскольку W„(/'A) и En(jA) являются функциями комплекс ного аргумента, их можно представить в алгебраической и по казательной формах:
W„ (JA) = U„ (А) + jV„ (А) = |
W„ (А) ем*); |
Ен(./ А ) == S„{A) + |
jQ H(А) = |
Ен (А) |
. |
Рассмотрим передаточные |
функции |
некоторых |
нелинейных |
звеньев.
З в е н о с и д е а л ь н о й р е л е й н о й х а р а к т е р и с т и
кой (см. рис. 8.2,а). |
|
|
|
|
|
|
ф ( 4 ) = — ; |
W„VA)=W„{A) = - ^ ; |
? И ) = |
0; |
ТС |
|
|
|
|
7ГА |
|
|
|
|
|
4В |
|
|
|
(8.53) |
|
|
|
|
|
|
|
Амплитудные зависимости |
W H(Л) |
и |
Ен (А) приведены на |
рис. 8.31. График |
Ен(Л) |
построен |
на |
комплексной плоскости. |
|
|
|
|
|
|
/®н |
|
|
|
|
|
~Й |
ш |
й=0 |
|
|
|
|
|
|
flz |
/?7 |
|
|
|
|
|
|
|
V |
|
Рис. 8.31. |
Идеальный релейный элемент: |
и |
а — графики амплитудной характеристики WH(A) |
обратной |
амплитудной |
характеристики |
Ен (А)\ б |
— |
годограф обратной отрицательной передаточной функ ции нелинейного элемента £„ (]А)
Однако он построен со знаком минус ( — EH(JA)), что’бы удоб нее было его использовать при расчетах, связанных с опреде
лением амплитуды |
и частоты |
автоколебаний. |
Как |
видно из |
рис. 8.31,6, |
годограф |
— £ и (/А) |
расположен |
вдоль |
отрица |
тельной вещественной оси, так |
как — En {jA) |
в данном случае |
не содержит мнимой части. |
|
|
|
|
З в е н о |
с |
з о н о й |
н е ч у в с т в и т е л ь н о с т и |
(см. рис. |
8.2,6). |
|
|
|
|
|
|
|
|
Ф(А) = ~ |
V A 2 - |
а2; W H(JA)= W H{A) = |
У Ж ^ Г ; |
|
■кА |
|
|
|
|
кА-1 |
|
Ф(А) = 0; |
|
|
|
- л 2 |
I |
(8.54) |
ЕкUА) — £■„ (А) = -—------ = = = = .. |
п ' |
|
н |
|
4 |
’ 4В У А2 — а2 |
|
На рис. 8.32 приведены амплитудные характеристики и годо граф — Е„ (JA ). Условно две сливающиеся ветви показаны раздельно. В действительности они обе располагаются на отри цательной вещественной оси.
а |
|
Рис. 8.32. Релейный элемент |
с зоной |
нечувствительности: |
— амплитудная |
характеристика |
WH(Л); |
б — обратная амплитуд |
ная |
характеристика £„ (Л); |
в — годограф |
обратной |
отрицательной |
|
|
|
передаточной функции нелинейного элемента ЕИ(/А) |
|
Н е л и н е й н ы й э л е м е н т с н е о д н о з н а ч н о й х а р а к |
т е р и с т и к о й |
(см. рис. 8.3,6). |
|
|
|
|
|
Если на вход такого элемента подать синусоидальный сиг |
нал |
(рис. |
8.33), то |
выходная величина |
будет серией |
знакопе |
ременных |
прямоугольных |
импульсов (рис. |
8.33,а). Изменение |
знака |
выходной величины |
будет |
опаздывать |
на |
угол |
<р(А) = |
|
|
|
а |
|
|
|
|
|
|
входной вели |
= arc sin — по отношению к изменению знака |
чины. |
|
А |
|
|
|
|
|
|
|
|
С учетом этого сдвига |
|
|
|
|
|
|
|
|
|
|
|
Ф(У'Л)= — |
|
|
|
|
(8.55) |
|
|
|
|
|
|
TZ |
|
|
|
|
|
где |
|
<р (Al = |
arcsin —— |
фазовый сдвиг. |
|
|
|
|
|
|
|
|
А |
|
|
|
|
|
|