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

u_lectures

.pdf
Скачиваний:
52
Добавлен:
22.03.2016
Размер:
2.16 Mб
Скачать

 

 

 

 

 

 

223

 

 

 

 

 

W

(a) =

X 2

 

 

(27.15)

 

 

Н

 

 

 

 

 

X 1

 

 

 

 

 

 

 

 

 

W

( jω) = −

 

X 1

,

(27.16)

 

Л

 

 

 

 

 

 

 

 

 

X 2

 

 

где X 2 и X 1 - изображения Фурье для синусоидальных сигналов. Перемножим (27.15) и (27.16) и получим

WH (a)WЛ ( jω) =WР(a, jω) = −1,

(27.17)

при ω = ωП , a = aП .

Выражение (27.17) представляет необходимое условие для существования автоколебаний в системе и называется условием гармонического баланса.

Уравнение (27.17) согласуется с условием нахождения линейной системы на границе устойчивости. По критерию Найквиста в линейной системе будут автоколебания, если годограф WР(jω) пройдет через точку (1, j0),то есть при некоторой частоте ω = ωп, WР(jωп) = 1.

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

27.5.Графоаналитический метод определения параметров автоколебательного режима и его устойчивости

Практически пользоваться выражением (27.17) неудобно, так как требуется подобрать два параметра, поэтому задачу решают графически.

Для этого запишем (27.17) следующим образом

WЛ ( jω) = −

1

= −M Н (a),

(27.18)

WН (a)

 

 

 

где: –MН(a) – обратный комплексный коэффициент передачи НЗ, который равен

M Н (a) = −

1

 

.

(27.19)

q(a) +

 

jq (a)

 

 

Уравнение (27.18) можно решить графически. Для этого нужно построить на комплексной плоскости годографы линейной части системы

224

WЛ ( jω) =U (ω) + jV (ω)

иобратного комплексного коэффициента передачи НЗ

M H (a) =U (a) + jV (a) .

Вточке пересечения годографов (рис. 27.13)

Im

Re

a

aп

ω

 

 

ωп

–Mн(a) Wл(jω)

Рис. 27.13

выполняется условие (27.18) и могут существовать автоколебания.

Параметры их определяются следующим образом: значение частоты ωп определяют из годографа линейной части WЛ(jω), а значение амплитуды aп определяется из годографа

MН(a).

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

Анализ устойчивости системы и устойчивости предельного цикла.

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

Aл(ω) e j[ϕл (ω)+ϕн (a)] = −1.

Aм(а)

Здесь Ам(a) = |Mн(a)| – модуль передачи нелинейного элемента,

ϕH (a) = argWH (a) .

Уравнение (27.20) можно представить в виде двух уравнений

Aл(ω) = 1 , Aм(а)

φл(ω) + φн(а) = – (2m+1) π ; где m = 0, ±1, ±2,…

(27.20)

(27.21)

Допустим, получили следующее расположение годографов (рис. 27.14).

D0 Im

D0

Im

 

Re

 

 

Re

D1

1`→ ← 1``

D1

Mн(a) Wл(jω)

1

2

Wл(jω)

 

–Mн(a) a

 

 

Рис. 27.14

 

Рис.

27.15

 

 

225

Пусть весь годограф нелинейного звена лежит в области D0, то есть

условие (27.20) не выполняется, так как Aл(ω) <1, значит |Wp(a, jω)| < 1 и в

Aм(а)

соответствии с критерием Найквиста система устойчива и в ней будут затухающие колебания.

Аналогично можно показать, что если весь годограф –Mн(а) лежит в области D1, то система будет неустойчива и в ней будут существовать расходящиеся колебания.

Если годографы пересекаются (рис. 27.15) , то в системе возможны автоколебания в точках 1 и 2. Определим, где автоколебания устойчивы?

В точке 1– амплитуда – ап; в точке 1`– амплитуда – п + ∆a); в точке 1``– амплитуда – ( ап – ∆a).

Пусть существуют автоколебания в точке 1 с параметрами ап и ωп. Допустим, произошло увеличение амплитуды ап + ∆a (точка 1`). Рабочая точка находится в области D0, то есть система устойчива и амплитуда колебаний уменьшится до ап.

Пусть амплитуда колебаний уменьшается, то есть ап – ∆a (точка 1``). Рабочая точка находится в области D1, то есть система неустойчива и амплитуда колебания возрастет до ап.

Следовательно, точка 1 соответствует устойчивому предельному циклу.

Рассмотрим точку 2.

Пусть а понижается, то есть ап – ∆a. В этом случае система устойчива и колебания затухающие. Следовательно, точка 2 соответствует неустойчивому предельному циклу.

На основе проведенного анализа сформулируем правило определения устойчивости автоколебаний: автоколебания устойчивы, если годограф

WЛ(jω) не охватывает точку на годографе

–MH(a),

в которой a = aп + a,

a > 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Метод гармонического баланса.

k=10c-1

xвх=0

x1

 

 

 

 

 

 

x2

 

 

 

X

 

 

 

 

 

 

 

 

 

b=0.5,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К/Р

 

 

вых

 

 

 

 

 

 

 

 

 

 

 

c=10.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 27.16

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

Комплексный коэффициент передачи нелинейного звена будет равен

Wн(а)=q(a)+jq`(a),

где

 

 

 

 

 

 

 

 

226

q(a) =

4C

1

b2

,

a b,

πa

a2

 

 

 

 

 

 

q`(a) = −

4Сb

,

 

 

a b.

 

 

πa2

 

 

 

 

 

Тогда обратный комплексный коэффициент передачи примет вид

Mн

(a) = −

 

 

1

 

= −

 

 

 

1

 

 

= −

πa

1

b2

+ j

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (a)

 

 

 

 

 

2

 

a

 

 

 

 

4C

 

 

b

b

4C

 

 

a

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

Построим его

 

 

 

 

πa

( 1a2

j a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) a = b = 0.5

πa

j b

= − j 3.14 0.5

 

 

 

 

 

 

 

Mн(a) = −

= − j0.04

 

 

 

 

 

2) a>b

 

4C

 

b

 

 

 

 

4C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

н(a) = −

πa

 

1

b2

j

πa b

 

 

 

 

 

4C

 

a

2

4C a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соnst

Годограф линейной части строим по выражению

Wл ( jω) =

k

= − j

k

= − j

10

 

jω

 

ω

 

ω

 

ω = 100

 

Wл(jω) = -j0.1

ω = 250

 

Wл(jω) = -j0.04

ω = 500

 

Wл(jω) = -j0.02

ω = 1000

 

Wл(jω) = -j0.01

Годографы

построены на рис. 27.17.

 

Im

Re

 

ω=1000

-0.01

 

ω=500

 

WЛ(jω)

 

 

а

ω=250

-0.03

-0.04

 

 

 

Н(а)

a=b=0.5

 

Рис. 27.17

Автоколебания возможны и имеют параметры:

a = b = 0.5; ω =250 рад/с .

28. МЕТОД ЛЯПУНОВА

28.1. Понятие о знакоопределенных, знакопостоянных и знакопеременных функциях

 

 

227

 

 

Пусть

имеется

функция

нескольких

переменных

V =V (x1, x2 ,...., xn ) .

 

 

 

При n = 2

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

при n = 3 – в

трехмерном

пространстве.

В каждой точке

пространства

функция V имеет определенное значение.

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

x1 = x2 =…= xn = 0,

V (0,0,....,0) = 0.

Кроме того, V должна быть непрерывна в некоторой области вокруг начала координат.

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

Для n = 2 такая функция может иметь вид

V = x12 + x22 .

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

Пример для n = 3

V = (x1 + x2 )2 + cx32 .

Обращается в нуль на прямой x2 = -x1 и x3 = 0.

Функция V называется знакопеременной, если в некоторой области вокруг начала координат она может иметь любой знак

Для n = 2 V = x1 + x2 .

28.2. Функция Ляпунова и её производная по времени

Нелинейная система может быть описана системой нелинейных уравнений:

dxdt1 = F1(x1, x2 ,.., xn ) dxdt2 = F2 (x1, x2 ,.., xn )

K

dxdtn = Fn (x1, x2 ,.., xn )

где x1 , x2 ,..., xn - переменные состояния системы.

F1 = F2 =... = Fn = 0 при

x1 = x2 =... = xn = 0

Любая функция

V =V (x1, x2 ,...., xn )

(28.1)

(28.2)

(28.3)

228

называется функцией Ляпунова, если в качестве переменных

x1 , x2 ,..., xn взяты переменные системы (28.1).

 

 

Производная по времени функции Ляпунова имеет вид

 

 

dV

=

V

dx1

+

V

dx2

+... +

 

V

dxn

 

(28.4)

 

dt

x

 

 

 

x

 

 

 

 

 

dt

 

 

 

x

2

 

dt

 

 

 

n

dt

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

dxi

 

 

 

Подставим в (28.4) выражения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dV

=

V

F +

V

 

F +... +

V

 

F

 

(28.5)

 

dt

x

 

 

 

 

 

 

 

1

 

 

x

2

 

 

2

 

x

n

 

n

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

=W (x , x

2

,..., x

n

)

 

 

 

 

 

 

 

 

 

 

(28.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

 

 

является

функцией отклонений.

Причем,

 

 

 

 

 

 

 

 

 

 

 

 

dt

W = 0 при x1 = x2 =... = xn = 0.

 

 

согласно

 

свойства

 

(28.2),

Поэтому

к этой

функции можно применить понятие знакоопределенной,

знакопостоянной и

знакопеременной функции.

 

 

 

 

 

 

 

 

 

 

 

Теорема Ляпунова. Если при

 

 

заданных в форме

(28.1) уравнениях

нелинейной системы n-го порядка можно подобрать такую знакоопределенную функцию V =V (x1, x2 ,...., xn ) , чтобы её производная по

времени W (x1, x2 ,..., xn )

 

 

 

 

также

была

знакоопределенной

или

знакопостоянной, но имела знак, противоположный знаку функции V, то

данная система устойчива.

 

 

 

Пример. Пусть задана нелинейная САУ, которая описывается системой

уравнений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx1

 

 

= F (x

, x

2

, x

3

)

 

 

 

 

 

 

 

 

 

dt

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

= F

(x

 

, x

2

 

, x

3

)

 

 

(28.7)

 

 

 

 

 

 

 

 

dt

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx3

= F (x , x

2

, x

3

)

 

 

 

 

 

 

 

 

 

dt

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выберем Функцию Ляпунова в виде:

 

 

 

V = а2 x2

+b2 x2

+ c2 x2

 

 

(28.8)

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

V > 0 при любых xi,

а,b,c - произвольные вещественные числа.

Будем придавать этой функции постоянные значения 0,С12,…., причем 0<С12… Тогда получим

 

 

 

 

 

229

V = а2 x2

+b2 x2

+ c2 x2

= 0,

 

 

1

2

3

 

 

 

V = а2 x2

+b2 x2

+ c2 x2

= C ,

(28.9)

1

2

3

1

 

 

V = а2 x2

+b2 x2

+ c2 x2

= C

.

 

1

2

3

2

 

 

Х3

MС2 С1

Х2

Х1x1

Рис. 28.1

 

Возьмем производную от V по времени

 

 

 

 

 

 

dV = 2a2 x F (x , x , x ) + 2b2 x F (x , x , x ) + 2c2 x F (x , x , x )

(28.10)

dt

1

1

 

1

2

3

2

2

1

2

3

3

3

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если полученная функция W (x1, x2 ,..., xn ) является знакоопределенной

отрицательной, то есть

 

 

 

 

 

 

 

 

 

 

 

 

W =

dV

 

< 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то V убывает во всех точках пространства, кроме начала координат,

 

V

 

 

 

 

 

 

 

 

Таким образом, при любых начальных

 

 

 

 

 

 

 

условиях изображающая точка M будет двигаться в

 

 

 

 

 

 

 

 

сторону уменьшения функции V, то есть будет

 

 

 

 

 

 

 

 

пересекать эллипсоиды и стремиться к началу

 

 

 

 

 

 

 

 

координат. А это означает, что все отклонения

 

 

 

 

 

 

 

t

системы x1 ,x2 ,x3 уменьшаются, то есть

система

 

Рис. 28.2

 

 

устойчива.

 

 

 

 

 

 

 

 

 

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

Однако, несмотря на это, метод применяется не только для исследования нелинейных систем, но и при анализе и синтезе адаптивных законов управления. Адаптивные системы управления относятся к классу нелинейных систем.

28.3. Формулировка теоремы Ляпунова в векторно-матричной форме

230

В общем случае, если используется математическое описание САУ в форме уравнений состояния, то функция Ляпунова часто имеет вид квадратичной формы, которая включает в себя члены второй степени по xi и xj

n

n

 

V (x1, x2 ,....xn ) = ∑∑pij xi xj ,

(28.11)

i=1

j=1

 

где: x - вектор состояния, pi,j коэффициенты.

В векторно-матричной форме

V (x , x

,..., x

n

) = xT Px

(28.12)

1 2

 

 

 

где xT – вектор-строка, полученная транспонированием вектора x;

P- квадратная определенно положительная симметрическая матрица типа n×n, элементы которой постоянны, причем P = PT .

Критерий Сильвестра.

Квадратичная форма (28.12) положительно определенная, если все “n” определителей |P1|,|P2|,…,|Pn| - положительны.

Пример :

V (x , x

 

) = [x x

 

] p11

p12

x1

 

= [x x

] p11x1 + p12 x2

 

= p x2

+ p x x

 

1

2

1

2

p

p

22

x

2

 

1

2 p x

+ p x

2

 

11 1

12 1

2

 

 

 

 

21

 

 

 

 

21 1

22

 

 

 

 

+ p12 x1x2 + p22 x22 = p11x12 + p22 x12 + 2 p12 x1x2

так как р12 = р21.

Определение устойчивости линейной системы.

Дифференцирование функции Ляпунова в виде (28.12).

Пусть имеем линейную автономную систему, которая описывается уравнением состояния

x = Ax

 

 

 

 

 

 

 

 

 

 

 

(28.13)

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда сопряженная система будет иметь вид

 

 

 

&T

= x

T

A

T

 

 

 

 

 

 

 

 

(28.14)

x

 

 

 

 

 

 

 

 

 

 

Пусть V( x )=xTPx

- функция Ляпунова.

 

 

Тогда найдем производную функции Ляпунова по времени

 

&

 

 

T

Px + x

T

Px = x

T T

T

PAx = x

T

T

(28.15)

V (x) = x

 

 

A Px + x

 

 

(A P + PA)x

 

 

&

 

 

 

 

 

&

 

 

 

 

 

 

 

Введем обозначение AT P +PA= - Q.

 

 

( 28.16)

Тогда если матрица Q положительно определена, то система (28.13) асимптотически устойчива, так как

231

W =V&(x) = −xT Qx < 0

(28.17)

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

29.УСЛОВИЯ АБСОЛЮТНОЙ УСТОЙЧИВОСТИ

НЕЛИНЕЙНЫХ СИСТЕМ

29.1. Понятие абсолютной устойчивости

Пусть имеем нелинейную систему с одной нелинейной функцией

x2=F(x1)

(29.1)

Пусть на характеристику НЗ наложены ограничения

 

0

F(x1 )

< kн

(29.2)

 

 

x1

 

232

Условие (29.2) графически изображено на рис. 29.1 и показывает, что нелинейная функция должна находиться внутри заштрихованного сектора. Примером такой функции могут служить характеристики АЦП и ЦАП.

x2

F(x1)

Условие

абсолютной

 

устойчивости гласит, что при любой

 

 

 

аrctg kн

нелинейности,

удовлетворяющей

 

 

ограничению (29.2), нелинейная

 

x1

система будет устойчива.

Рис. 29.1

29.2. Частотный критерий абсолютной устойчивости В. М. Попова

Пусть нелинейная система содержит однозначную нелинейность и в ней выделены линейная часть и нелинейное звено (рис. 29.2).

x1

 

x2

F(x1)

 

 

 

 

 

Линейная часть

Рис. 29.2

НЗ описывается уравнением (29.1) и удовлетворяет ограничениям

(29.2).

Пусть линейная часть имеет передаточную функцию (29.3) и

описывается уравнением (29.4):

 

W ( p) =

X1

( p)

= −

R( p)

 

(29.3)

X2

( p)

Q( p)

 

 

 

Q( p)X1 ( p) = −R( p)X2 ( p),

(29.4)

где Q(p) и R(p) – полиномы степени n и m

 

Q(p)=a0pn + a1pn-1 +…+ an

(29.5)

R(p)=b0pm + b1pm-1 +…+ bm

 

На линейную часть накладываются следующие ограничения: Характеристическое уравнение линейной части Q(p) = 0 имеет все

простые корни с отрицательной вещественной частью и не более двух

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]