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

Лабораторные / UP_Nelineynye_sistemy_upravlenia

.pdf
Скачиваний:
20
Добавлен:
16.05.2021
Размер:
2.6 Mб
Скачать

s

f(s)

M

N

L

0

sц

s0

s

 

 

 

Рис. 2.12. Графическое изображение функции последования

Возьмем теперь исходную точку s0 левее sц и тем же способом проследим ход точечного преобразования, как показано стрелками на рис. 2.12. Очевидно, что этот процесс тоже сходится к тому же предельному циклу sц . Следовательно, мы имеем устойчивый предельный цикл (автоколебания).

Отсюда условие устойчивости предельного цикла имеет вид:

ds

1.

 

 

 

 

ds s s

 

 

 

ц

 

В противном случае, изображенном на рис. 2.13, а, получается неустойчивый предельный цикл. На других графиках рис. 2.13 показаны: б – случай двух предельных циклов, где один – устойчивый, а второй – неустойчивый; в – случай расходящихся колебаний; г – случай затухающих колебаний.

Такого типа графики называются диаграммами точечного преобразова-

ния. Изображение хода точечного преобразования на такой диаграмме эквивалентно сопряжению начальных и конечных условий соседних участков в методе припасовывания. Но производится это специальным и довольно простым геометрическим построением.

31

s

s

0

s

ц

s

0

s

ц1

s

ц2

s

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

б

 

 

 

s

 

 

 

s

 

 

 

 

0

s

0

s

 

 

 

в

 

г

Рис. 2.13. Нарушение условия устойчивости предельного цикла:

а неустойчивый предельный цикл; б два предельных цикла, из которых один – устойчивый, а второй – неустойчивый; в расходящиеся колебания;

г затухающие колебания

Основным в методе является нахождение функции последования s f s на основе решения уравнений динамики системы (2.36). Найти эту функцию в явной форме не всегда легко. Поэтому иногда используют представление функции последования в параметрической форме.

Рассмотрим данный метод на примере генератора незатухающих прямоугольных сигналов, схема которого показана на рис. 2.14.

Данный генератор относится к классу систем с положительной обратной связью, в которых при определенных условиях могут возникать незатухающие колебания.

32

Uвых

L

iL

C

+

R

iC

Eп

 

 

 

iк

VT

Рис. 2.14. Генератор незатухающих сигналов

Учитывая связь тока коллектора транзистора VT с токами, протекающими через вторичную обмотку трансформатора и конденсатор, которые в свою очередь образуют колебательный RLC-контур, составим систему уравнений этого контура:

i

 

i

L

i ,

 

 

к

 

 

 

C

 

u

 

i

 

 

R L diL .

(2.39)

 

 

C

 

 

L

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

С другой стороны, можно записать

uC

1

iC dt .

(2.40)

C

 

 

 

Тогда с учетом (2.39) и (2.40) уравнение колебательного контура можно переписать так:

iL R L

diL

 

1

iC dt .

(2.41)

 

 

 

dt

C

 

 

Дифференцирование обеих частей (2.41) и дальнейшая подстановка первого уравнения (2.39) даст:

LC

d 2iL

RC

diL

i

L

i .

(2.42)

 

 

 

dt 2

 

dt

к

 

 

 

 

 

 

 

 

33

 

 

 

 

При этом ток коллектора в данной схеме определяется падением напряжения на переходе база-эмиттер:

 

 

 

 

 

 

 

 

 

 

u

БЭ

M

diL

.

 

 

 

 

 

(2.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем полагать, что ток коллектора изменяется по идеальной характери-

стике, определяемой (2.43):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iк

I , uБЭ 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.44)

 

 

 

 

 

 

 

 

 

 

 

0, uБЭ 0,

 

 

 

 

 

 

где I Eп

R .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда уравнение колебательного контура (2.42) с учетом (2.43), (2.44) и

замены iL x можно переписать в рассмотренном ранее виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

I , x 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2hx 02 x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, x 0,

 

 

 

 

 

 

где 2h R L , 2

LC 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

Фазовые траектории

в верхней

и

нижней

 

 

 

 

 

 

 

 

полуплоскостях представляют собой спирали,

 

 

 

 

 

 

x

 

скручивающиеся к точке

x I в верхней полу-

 

b

 

 

 

 

 

плоскости и к началу координат в нижней. Рас-

 

 

 

 

 

x0 = a

 

 

 

 

 

 

 

смотрим

 

движение,

начинающееся

в

точке

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0 a на оси абсцисс, как показано на рис. 2.15.

Рис. 2.15. Фазовые траектории

 

 

Решение уравнения для фазовой траекто-

 

колебательного контура

 

 

 

 

 

рии, выходящей из этой точки и проходящей в

 

 

 

 

 

 

 

 

нижней полуплоскости, будет иметь вид:

 

 

 

 

 

 

 

 

 

x e ht c

cos t c

2

sin t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

cos t

 

 

 

sin t e ht c sin t c

 

cos t .

 

y he ht

c

2

2

 

 

 

1

0

 

 

 

0

 

 

 

 

 

 

1

0

0

0

 

0

 

Подставляя в (2.45) начальные условия x0 x 0 a , y0

 

 

 

(2.45)

 

y 0 0 , по-

лучим значение постоянных интегрирования c1 и c2 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

 

 

 

 

 

 

 

 

(2.46)

 

 

 

 

 

 

 

 

 

 

 

c

 

hx0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

Учитывая (2.46), уравнения движения из (2.45) окончательно могут быть переписаны как:

 

 

 

 

 

 

 

 

 

h

 

 

 

 

ht

 

 

 

 

 

 

x x0e

 

 

cos 0t

 

 

sin

0t ,

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

(2.47)

y x

 

 

h

 

 

 

 

e

ht

 

2

sin t .

 

 

 

 

0

 

 

2

0

 

 

0

 

 

 

 

 

 

0

 

 

 

 

Для продолжения траектории в верхней полуплоскости необходимо знать значение абсциссы изображающей точки при пересечении с осью абсцисс. Из второго уравнения (2.47) определим момент времени, в который происходит пересечение с осью абсцисс из условия y 0 sin 0t . Таким образом 0t или t 0 . Подставляя получившееся значение времени t 0 в первое уравнение (2.47), получим

 

 

 

 

 

h

 

 

 

 

x0e

h 0

 

 

 

 

x

0

 

cos

0

sin

 

 

 

 

 

или

 

 

 

 

 

e h 0

 

 

x

 

 

x

0

b.

(2.48)

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Можно показать, что движение в верхней полуплоскости будет проис-

ходить в течение времени 0

t 2 0 .

При этом начальные условия

x0 x 0 b , y0 y 0 0 . Тогда можно показать, что уравнение движения имеет вид

h t h

x I e 0

0

По истечении времени цисс справа в точке

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b I sin 0 t

 

 

b I cos 0 t

0

.

 

0

 

 

t 2 0 траектория снова пересечет ось абс-

 

2

 

 

h 0

b I I c .

 

 

 

e

(2.49)

x

0

 

 

 

 

 

 

 

 

Введя дополнительное обозначение e h 0 , из (2.48) и (2.49) можно получить систему уравнений, связывающую точки пересечения фазовыми траекториями оси абсцисс:

35

c b 1 I ,

(2.50)

 

b a.

 

Система (2.50) определяет процесс преобразования точки а на положительной полуоси абсцисс в точку с на этой же полуоси. Если с < a, то траектории будут скручиваться, и генератор будет совершать затухающие колебания. Если с > a, то траектории будут раскручиваться, а колебания – нарастать. Если с = a, то колебания становятся незатухающими и траектория превращается в замкнутый цикл. Уравнения (2.50), называемые уравнениями точечного преобразования, позволяют найти параметры этого цикла. Так, положив с = a, решим систему (2.50) относительно a и b:

 

I

 

 

a

 

 

,

 

1

 

 

 

(2.51)

 

I

 

 

b

.

 

 

 

 

1

 

 

 

 

 

Если 1, то решение (2.51) единственное. Таким образом, найденный замкнутый цикл является изолированным. Изолированную замкнутую траекторию называют предельным циклом. Предельный цикл окружен навивающимися на него или скручивающимися с него траекториями. Если в результате малого отклонения от цикла в любом направлении изображающая точка попадает на траекторию, неограниченно приближающуюся к циклу, то цикл устойчив. Устойчивый предельный цикл соответствует устойчивым колебаниям, называемым автоколебаниями.

Пусть связь между точками с и a определяется из (2.50) как c T a . Для рассматриваемого случая

c T a 2a 1 I .

(2.52)

Если в результате возмущения точка сместится с цикла, то по истечении

периода получим приращения a и c , связанные зависимостью

 

 

c

dT

 

 

lim

 

 

.

 

 

 

 

a

da

 

 

Соответственно, если dTda 1, то c a , т. е. начальное отклонение уменьшается и предельный цикл устойчивый, в противном случае – неустойчивый.

Для (2.52) T a 2 e2h 0 1, т. е. цикл устойчивый.

36

2.3.4. Метод припасовывания (сшивания) граничных условий

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

Полученные линейные уравнения реша-

 

 

 

 

R( p)

 

x

 

 

 

 

W

 

ются соответственно по участкам статической

 

 

 

 

Q( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

нелинейности обычным или

операционным

y

 

 

 

 

 

 

 

 

 

 

методами. Полученные решения содержат по-

 

 

M

 

 

 

 

 

-M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

стоянные интегрирования,

зависящие от

Рис. 2.16. Структурная схема

начальных значений движения. Для нахожде-

 

 

нелинейной системы

 

 

 

 

 

 

 

 

 

 

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

ния последующего участка, т. е. участки решения припасовываются друг к

другу.

y

 

 

 

 

 

Рассмотрим использование метода припа-

M

 

 

 

совывания на примере схемы, приведенной на

 

 

T/2

T

t

рис. 2.16 и состоящей из линейной части и

 

 

 

 

 

 

 

 

 

двухпозиционного реле в обратной связи.

 

M

 

 

 

 

 

 

x

 

 

 

Допустим, что в приведенной системе

 

 

 

 

 

 

 

 

существуют симметричные автоколебания, вид

 

 

 

 

 

которых показан на рис. 2.17. Требуется дока-

 

 

 

 

t

 

 

 

 

зать их существование и определить период

 

 

 

 

 

автоколебаний.

 

 

 

 

 

Рассмотрим движение системы на началь-

 

 

 

 

 

 

 

Рис. 2.17. Переменные

 

ном интервале времени при 0 t T 2 . С уче-

 

 

 

 

 

нелинейной системы

 

 

том отрицательной обратной связи решение системы x t

запишется как ре-

акция на воздействие y M :

 

n

 

xI t C0 Cie pit ,

(2.53)

i 1

 

где Ci – неизвестные постоянные интегрирования.

 

37

На втором интервале времени T2 t T решение xII t можно найти как алгебраическую сумму двух решений:

продолжения решения xI t , считая, что воздействие y M на входе не снимается;

решения, полученного при подаче на вход линейной части ступенчатого воздействия y 2M , компенсирующего первое решение.

В результате получим решение

n

 

 

xII t C0 Cie pi t T 2 x2

t ,

(2.54)

i 1

где третье слагаемое x2 t определяется согласно формуле Хевисайда

 

t 2M

R 0

n

R pi

p t

 

x2

 

 

2M

 

 

e i .

(2.55)

Q 0

 

 

 

 

i 1 piQ pi

 

 

Суммарный процесс на втором (2.55) можно записать так:

x t C n C e pi t T 2

II 0 i i 1

интервале времени с учетом (2.54) и

 

R 0

n

R pi

p t

2M

 

 

2M

 

 

e i . (2.56)

Q 0

 

 

 

i 1 piQ pi

 

В силу симметрии характеристики двухпозиционного реле, если автоколебания существуют, то решение на втором интервале должно быть равно решению на первом, взятому с обратным знаком, т. е.

 

 

 

 

 

 

 

 

xII t xI t .

 

 

 

 

 

 

 

 

 

(2.57)

Из (2.53), (2.56) и (2.57) следует:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

p

t T 2

 

 

 

 

R 0

 

 

n

 

 

R pi

p t

 

 

n

p t

 

 

 

 

 

 

 

 

 

 

 

 

 

C0 Cie i

 

 

2M

 

 

2M

 

 

 

 

 

e i

 

C0 Cie i

 

 

Q 0

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

i 1 piQ pi

 

 

 

 

 

i 1

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R 0

 

n

 

 

piT 2

 

1

2MR pi

 

 

 

pit

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2C0 2M

 

 

 

 

Ci e

 

 

 

 

 

 

e

 

0 . (2.58)

Q 0

 

piQ pi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку равенство (2.58) должно быть справедливым для любых значений t, то необходимо порознь приравнять нулю свободный член и множи-

тели при epit. Тогда получим систему значений постоянных интегрирования

38

(x)
Рис. 2.18. Структурная схема нелинейной системы
W(p)
x

 

 

M

R 0

 

 

 

 

 

C0

 

 

,

 

 

 

 

Q 0

 

 

 

 

 

 

 

 

 

 

 

(2.59)

 

 

 

 

2MR pi

 

C

 

 

 

, i 1 n.

 

 

 

 

 

i

 

piQ pi e

piT

2

1

 

 

 

 

 

Таким образом, все постоянные интегрирования выражены через известные коэффициенты и полюсы передаточной функции и через неизвестный период колебаний T. Для нахождения периода автоколебаний составим уравнение периодов, для чего используем условия переключения: на границах интервалов x проходит через нуль. С учетом этого, а также выражения (2.53) для конца первого интервала в момент времени t T2 справедливо равенство

 

 

 

 

n

 

 

 

 

 

 

C0 Cie piT

2 0 .

 

(2.60)

 

 

 

 

i 1

 

 

 

Подставив C0 и Ci из (2.59)

в (2.60),

получим уравнение

периодов

 

R 0

n

2MR p e piT 2

 

 

M

 

 

 

i

 

0 .

(2.61)

 

 

 

 

 

Q 0

i 1 piQ pi e piT 2 1

 

 

Уравнение (2.61) содержит только одно неизвестное T. После нахождения периода автоколебаний необходимо сделать проверку выполнения условий переключения. Если окажется, что период автоколебаний является вещественным числом, то автоколебания существует, если же комплексным, то вероятно, что колебания искомой формы отсутствуют.

2.3.5. Метод гармонической линеаризации

Метод гармонической линеаризации относится к приближенным аналитическим методам исследования нелинейных систем. Метод гармонической линеаризации предназначен прежде всего для исследования периодических (автоколебательных) процессов в нелинейных системах, однако возможно применение метода и для исследования колебательных переходных процессов, а также для решения широкого круга других задач.

Для рассмотрения метода гармониче- y ской линеаризации на рис. 2.18 приведем структурную схему нелинейной системы с одним нелинейным блоком в обратной связи.

39

H W j
y cx f x .

В основе метода гармонической линеаризации лежат две гипотезы:

гипотеза о порождающей системе;

гипотеза фильтра.

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

(2.62)

При обращении малого параметра в нуль уравнение (2.62) вырождается в порождающее линейное.

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

должно быть выполнено неравенство H H k , где k 2, 3, 4 ..., т. е. коэффициент передачи для основной частоты значительно превосходит коэффициент передачи для высших частот.

Уравнения, описывающие приведенную систему (рис. 2.18), имеют следующий вид:

x yW p ,

(2.63)

y x ,

(2.64)

где W p – передаточная функция линейной части системы от входа y к выходу x ; x R – выход линейной части системы; y R – поступающий на вход линейной части выход нелинейной части со статической характеристикой x .

Пусть в системе, задаваемой уравнениями (2.63) и (2.64), имеет место периодический процесс с некоторой частотой и периодом T 2 . Прежде всего интересно определение характеристик этого процесса (ампли-

40

Соседние файлы в папке Лабораторные