to можем считать, что на выходе • нелинейного элемента (рис. 13.12) будет только первая гармоника из всего спектра гармоник
ряда Фурье, т. е. |
получим |
|
х =*■•А |
sin w / -4- Ь(Л) cos ш / 1, |
(13.16) |
где коэффициенты для -первой гармоники ряда, Фурье равны (для упрощения индексы «1» у g(A) и Ь{А) опущены):
2-к
g(A) — z^ | F(x1) sin idtdu>t\
О
(13.17)
2 ~
Ь(А) = ~д\^ F{x i) cos со tdo) t ,
a F (xx) — нелинейная функция, соответствующая релейной ха рактеристике нелинейного элемента.
Запишем изображения по Лапласу для входной и выходной величин нелинейного элемента:
х л р У - а р*-f- е>-
Х { р ) - л S W ^ s + b i A ) ^ ^
р* +
Следовательно, передаточная функция нелинейного элемента рав на '
|
т = |
Х(р) |
g(A) + b(A) £ |
(13.18) |
|
х л р ) |
|
|
|
Из выражения (13Л8) видно, что передаточная функция нелиней ного элемента в общем виде зависит от амплитуды входного сиг нала А и частоты со. Наиболее распространены случаи, когда пе редаточная функция нелинейного элемента зависит только от амп литуды- Л входного сигнала и не зависит от частоты. Заменив п выражении (13.18) р на /ш и .произведя необходимые сокращения, получим выражение, которое не зависит от частоты w . Значит, комплексная передаточная функция данного нелинейного элемента зависит только от амплитуды А входного сигнала. Комплексная
передаточная функция нелинейного элемента, которую теперь обозначаем через J(j А), равна
J<JA) = g { A ) + j b { A \ |
(13.19) |
Отсюда .можем записать уравнение для нелинейного |
элемента в |
следующем виде: |
|
x — J( jA)xv |
(13.20) |
При неизменных значениях амплитуды и частоты входного сигна
ла ./( / Л) = |
const и связь (13,20) оказывается линейной. |
Замена |
фактической нелинейной связи между входом и выхо |
дом нелинейного элемента линейной связью (13.20) носит назва ние гармонической линеаризации. Гармонической она называется потому, что выходная величина нелинейного элемента заменяется первой гармонической составляющей (первой гармоникой).
Определим передаточную функцию |
нелинейного элемента, |
имеющего |
однозначную |
релейную характеристику (без |
гистере |
зисной петли), |
при подаче на вход элемента входного |
сигнала |
(13.14). |
На |
выходе |
получим сигнал |
прямоугольной |
формы |
F{X\), |
из которого выделяем первую гармонику. Чтобы по |
лучить аналитическую зависимость первой гармоники, нужно вы числить коэффициенты ,?(Л) и й(Л), для чего необходимо знать пределы интегрирования. С этой целью на рис. 13.13 приве
ден порядок получения выходного сигнала |
х = |
Р(хг). Из рисун |
ка видно, что |
|
|
|
|
|
|
0 |
при |
j3> |
w > 0 |
|
|
|
+ с |
при |
( " — |
? ) > " > t > э |
|
|
0 |
при |
( т : + ^ > Ю* > 0 :- Э ) |
(13.21) |
— с |
при |
(2- - | 1 ) > ю * > ( « |
+ ?) |
|
0 |
при |
2 - > w г* > |
(2 т; |
р). |
|
Подставляя значения функции Р(х^) |
в формулу |
(13.17), можно |
вычислить коэффициенты передаточной функции нелинейного эле мента.
Так как функция /*’(^i) нечетная, то Ь{А) — 0 и передаточ
ная функция |
нелинейного элемента становится действительным |
числом, |
- . |
Коэффициент g)A) определяется следующим образом:
|
2~ |
|
3 |
|
- - Г 1 я-Ср 2 - - Э |
|
2 - |
^ ( Л ) = |
—^-г I |
F ( х , ) si n w /'rfw ^ = |
. |
|
+ I + 4~| | + |
|
|
я А . ] |
|
я Л |
|
|
|
|
|
О |
|
|
р |
я — р я р ф |
|
2 я — 'р |
|
|
|
|
|
2 |
я - р |
|
|
|
|
я Л |
О + с j sin w tci a>i 4 - 0 - |
с |
|
J |
sin w t dxо t -f |
0 |
|
g(A)= |
c o |
s p = ’ я Л \ 1 ~ |
5in2 p ■ |
|
|
|
|
|
|
|
|
|
Рис. |
13.13. |
Из рис. |
13.13 видно, |
что |
а := Л sin К поэтому sin(3 = -^-. Отсюда |
|
e W |
^ |
V ' A |
( 1 3 . 2 2 ) |
|
* |
Несколько более сложными получаются.коэффициенты £(Л) |
и b (Л) |
при входном сигнале х шфй. |
Аналогичным образом |
вычисляются |
коэффициенты ^(<4) и |
Ь(А) для других нелинейных элементов1. |
нелинейной сис |
Перейдем к определению автоколебаний в |
теме, изображенной на рис. |
13.12. Будем |
считать, |
что линейная |
часть системы в разомкнутом состоянии является устойчивой. За
менив р на |
усе |
в передаточной функции линейной части W(p) = |
Х.,(Р) |
> |
, |
= — |
получим выражение для амплитудно-фазовои ха |
рактеристики разомкнутой системы W(ju>), иногда называемой комплексной передаточной функцией. Поэтому выходная величина
линейной части х2 будет |
равна |
|
|
|
Xs = W{J<*)x. |
|
(13.23) |
Подставляя |
в (13.23) значение х из (13.20), |
получим |
|
х., — W (у со) J (у Л) х г. |
' |
(13.24) |
Значит, |
W ( j w ) J ( j A ) |
является комплексной |
передаточной |
функцией разомкнутой линеаризованной нелинейной системы. Как известно, для определения устойчивости системы, в том
числе и определения незатухающих периодических колебаний, не обходимо составить характеристическое уравнение замкнутой сис темы. На этом основании для нахождения автоколебаний в нели нейной системе составим уравнение собственных колебаний лине аризованной нелинейной системы, по своему виду напоминающее характеристическое уравнение замкнутой системы. Уравнение
собственных колебаний системы |
будет |
|
W ( у ю) У ( у Д ) + |
1 = 0 . |
( 1 3 . 2 5 ) |
Аналитическое решение уравнения (13.25) позволяет определить частоту о и амплитуду А автоколебаний. Для этого необходимо из комплексного числа (левой части уравнения 13.25) выделить вещественную и мнимую части, приравнять их раздельно нулю и
решить |
полученные |
уравнения. |
|
|
|
|
|
|
Однако значительно проще уравнение (13.25) решается гра |
фически. |
|
Для этого |
перепишем |
его так: |
|
|
|
|
|
|
Щум) = |
|
1 __ |
|
(13.26) |
|
|
|
|
J ( j A ) ‘ |
|
|
|
|
|
|
|
|
|
|
|
1 Э ти |
коэф ф ициенты |
вычислены и |
для |
ка ж д о го нел инейного |
элемента |
им е |
ю тся готовы е |
вы раж ени я |
передаточной |
ф ункц ии |
J (j А) , которы е |
м о ж н о |
найти |
в п р и л о ж е н и и |
3 кн и ги |
А . |
А . |
К р а с о в с ко го и |
Г. С. |
П оспелова [4 ] |
. П о это м у при |
и н ж ен ерны х |
расчетах |
м о ж н о |
пользоваться |
готовы м и передаточны м и ф ункц иям и |
нелинейны х |
элементов, т. е, |
нет необходим ости |
в ы ч и с л я т ь . интегралы (13.17). |
|
где |
|
1 |
часто |
обозначают через |
г U Л) |
и называют |
|
J{jA) |
|
|
амплитудной |
характеристикой |
нелинейного |
|
|
|
|
|
|
|
|
элемента. |
|
|
(13.25) по- |
|
|
Некоторые авторы преобразовывают уравнение |
|
другому |
и получают |
уравнение |
|
|
|
|
|
|
|
- |
W U “) = J U A ) - |
■ |
° 3'27) |
|
|
Для |
графического решения |
уравнения |
(13.26) |
необходимо |
построить на одной комплексной плоскости амплитудно-фазовую
характеристику W (у ы), |
изменяя частоту |
о ) от 0 |
до о © , |
и отрица |
тельную амплитудную характеристику |
z(jA), |
изменяя |
амплиту |
ду А от 0 до о о . Тогда, |
если эти характеристики |
пересекутся, то |
это значит, что в данной системе могут существовать автоколеба ния (система будет находиться на границе устойчивости по отно шению к первой гармонике автоколебаний). Если эти обе кривые не пересекутся, то в системе невозможны автоколебания.
На рис. 13.14 построены кривые W(J ю) и — z(y А) для слу чая, когда они пересекаются в двух точках: 7 и 2. Значит, в систе ме могут возникать автоколебания с частотой Wj и амплитудой
.4], соответствующие точке пересечения 7, причем частота ац бе рется из кривой IF(/'со), а амплитуда А\ — из кривой г (у Л). Могут также возникать автоколебания с частотой о>2 и амплиту дой А2, соответствующие точке пересечения 2.
Теперь необходимо выяснить, являются ли эти автоколебания устойчивыми? Для этого воспользуемся следующим критерием устойчивости, базирующимся на критерии устойчивости Найк виста:
Если амплитудная характеристика нелинейного элемен
|
|
|
|
|
|
|
|
та — z (j А), |
пересекая АФХ |
W (j со) линейной части |
систе |
мы, выходит из области, в которой характеристика |
— z ( j А) |
охватывается |
характеристикой |
W (y'w), то |
точка их |
пере |
сечения |
соответствует устойчивым автоколебаниям |
(точка 2 |
на рис. 13.14); если характеристика |
—z { j А) |
входит |
в об |
ласть, где характеристика — г (у А) |
охватывается |
характе |
ристикой |
W( jd)), то точка |
их пересечения |
соответствует |
неустойчивым |
автоколебаниям |
(точка 1 на рис. 13.14). |
|
Поясним |
несколько этот критерий |
устойчивости. |
Допустим, |
что в системе возникли автоколебания с амплитудой А\ и частотой Wj. Пусть по каким-либо причинам амплитуда А\ уменьшилась. Тогда автоколебания срываются (они по критерию в точке 1 неус тойчивы) и их амплитуда и частота уменьшаются до нуля. Это значит, что замкнутая система находится в устойчивом состоянии (в соответствии с критерием устойчивости Найквиста) и в ней отсутствуют автоколебания. Теперь допустим, что амплитуда А\ увеличилась, а так как автоколебания в точке 1 неустойчивы, то амплитуда и частота автоколебаний увеличиваются до Л2 и со,, т. е. устанавливаются устойчивые автоколебания, соответствующие точке 2.
Пусть в системе установились автоколебания с амплитудой
.42 и частотой со2. По какой-либо причине амплитуда А2 уменьши лась. Тогда амплитуда и частота снова увеличиваются до значе ния А2 и со2 , т. е. снова устанавливаются устойчивые автоколеба ния. То же самое происходит и при увеличении амплитуды А2, т. е. снова в системе устанавливаются устойчивые автоколебания, со ответствующие точке 2.
Если амплитуда устойчивых автоколебаний оказалась боль шой по величине, тогда нелинейная система считается практиче ски неустойчивой. Для того чтобы система стала практически.ус тойчивой, надо в систему ввести корректирующее устройство, ко торое уменьшит амплитуду автоколебаний или совсем ее сведет к нулю, т. .е. приведет к ликвидации автоколебаний в системе. Кор ректирующие устройства, вводимые в нелинейные системы, позво ляют улучшить качество переходного процесса в системе.
Метод гармонического баланса, предложенный Е. П. Поло вы.ч, дает возможность исследовать как устойчивость системы и автоколебаний, так п качество процессов регулирования нелиней ной системы.
П Р И Л О Ж Е Н И Е 2
ПРЕОБРАЗОВАНИЕ ЛАПЛАСА И ЕГО ПРИЛОЖЕНИЕ К РЕШЕНИЮ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИИ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ
§1 . Общие сведения о преобразовании Лапласа
Вданном курсе основное внимание уделяется линейным систе мам, которые описываются обыкновенными линейными дифферен циальными уравнениями с постоянными коэффициентами. /Для ре шения таких дифференциальных уравнений в теории автоматиче ского регулирования широко применяется операционный . метод на базе одностороннего преобразования Лапласа из операционного исчисления.
Преобразование |
Лапласа позволяет |
перейти от |
дифферен |
циальных |
уравнений |
относительно интересующей нас |
переменной |
(функции |
времени t) |
к алгебраическим |
уравнениям относительно |
некоторой, функции комплексной переменной р. Применительно к системам автоматического регулирования преобразование Лапла са позволяет перейти от дифференциальной связи между входной и выходной величинами к алгебраической связи между их изобра жениями и учесть начальные условия системы. После того как ал гебраическое уравнение будет решено и будет получено выраже ние для изображения выходной величины через известную вход ную величину, совершается обратный переход от изображений к реальным функциям времени.
Функция / ( 0 называется оригиналом, а функция |
/. ]/(/)] — |
преобразованием Лапласа или изображением функции |
/(£). За |
пишем это так: |
|
£ ! / ( / ) ! - > / (О или L[f{t)\ = № , |
|
где знак /->- или = не определяет равенства изображения и ори гинала, а означает их соответствие, определяемое выражением
О С
прямого преобразования Лапласа L \f (t)\ - j / (/) е" pt d t .
О
В некоторых случаях прямое преобразование Лапласа запи |
сывают |
L [/’(01 |
= Л ( р ) а обратное преобразование Лапласа— |
f(t) = L-i\F(p)\. |
|
по оригиналу составляет прямую |
Определение изображения |
задачу |
теории преобразования |
Лапласа. Определение оригинала |
по изображению составляет обратную задачу теории преобразова
ния |
Лапласа. |
|
|
|
|
н правила опера |
|
Кратко сформулируем основные теоремы |
ционного исчисления на базе преобразования |
Лапласа, которые |
используются в теории автоматического регулирования. |
|
1. Определение изображения суммы нескольких функций по |
|
Если Z L1 |
(О!Н-/;(0; |
М.г(01“*/а(0;-.; |
Г/»(01-М»(О, |
изображениям |
этих |
функций. |
|
|
то |
i | / ](0 + / a ( 0 + |
. - . + / B(0 |
l = |
^ f/i( 0 1 + ^ l / 2 (0 1 - f ..'-+ ^ [/» (0 1 - |
|
2. Определение |
изображения |
производной |
от заданной функ |
Тогда |
L |
dfjt) 1 |
р М / ( 0 |
1 - |
/ ( 0 |
) . |
|
|
|
|
|
at |
|
|
|
|
|
|
Если |
/ (0 ) = 0 , |
то |
d f { t ) |
|
/’/-[/(0 1 ; |
|
|
d t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Д Д О |
P-L [/(01 |
Pf( 0)4- |
dm |
|
|
|
|
tW |
|
dt |
Если |
/ ( 0 ) |
d m . |
0 , |
t o |
L |
|
Pn [ f ( t ) \ . |
d t |
|
d t 2 |
|
|
|
|
|
|
|
|
|
Аналогично н для производных более высоких порядков. |
то |
Если |
/.[/(01> /(0. |
|
|
от определенного интеграла. |
|
3. |
Определение |
изображения |
|
|
|
|
|
|
' О |
|
|
|
|
|
|
|
|
|
L |
\ f { t ) d t |
1 / ( 0 1 • |
4. Теорема запаздывания. Теорема указывает, каким образом изменяется изображение при смещении оригинала вдоль оси t на постоянную величину т0.
Если
М Я 0 1 Д / ( 0 >
то
е_ХР/. 1/ ( 0 1 - > /(^ -- 0 .
5.Теорема умножения изображений, или теорема свертки двух функций,