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

Vedyakov_articles / Операторный метод анализа и синтеза ЛСУ

.pdf
Скачиваний:
27
Добавлен:
21.03.2016
Размер:
1.12 Mб
Скачать

Y (s) = L[ y(t)]

к изображению входной переменной

U (s) = L[u(t)]

при нулевых начальных условиях

y(0) = py(0) =K = pn1 y(0) = u(0) = pu(0) =K= pm1 y(0) .

Напомним, что преобразование Лапласа

F(s) = L[ f (t)] = f (t)est dt

0

ставит в соответствие каждой функции f (t) (оригиналу), для которой несобствен-

ный интеграл сходится, единственную функцию F(s) (изображение) комплексной переменной s . Использование преобразования Лапласа для изучения дифференциальных уравнений основывается на утверждении:

L[ pk f (t)] = sk F(s) ,

если равны нулю значения f (t) и ее производных вплоть до ( k 1 )–ой при t = 0 .

Применяя преобразование Лапласа к правой и левой частям уравнения (2.1) при

нулевых начальных условиях, получим

 

a(s)Y (s) = b(s)U (s) ,

 

и, следовательно, передаточная функция -

 

W (s) =

b(s)

.

(2.3)

 

 

a(s)

 

Передаточная функция у которой полиномы a(s) и b(s) не являются взаимно простыми называется вырожденной передаточной функцией. При выводе передаточной функции W (s) системы необходимо следить за тем, чтобы не произошло со-

кращения ее числителя и знаменателя. При выполнении последнего условия, знаменатель a(s) передаточной функции называется характеристическим полиномом

системы. Устойчивость объекта управления определяется корнями характеристического уравнения a(s) = 0 . Для асимптотической устойчивости необходимо и доста-

точно отрицательность вещественных частей всех корней. Любой полином, все кор-

11

ни которого имеют строго отрицательную вещественную часть, называется устой-

чивым полиномом.

Корни знаменателя передаточной функции называются полюсами, а корни числителя - нулями. Передаточная функция, у которой вещественные части всех нулей отрицательны называется минимально-фазовой. Для существования решения задачи управления достаточно, чтобы передаточная функция была минимально-фазовой. Если W (s) не является минимально-фазовой передаточной функцией, то решение

задачи управления может существовать только при условии, что наибольшее общее кратное полиномов a(s) и b(s) - устойчивый полином.

Если W (s) является минимально-фазовой передаточной функцией, то решение

задачи будем искать в два этапа. Для этого представим сигнал управления в виде

 

 

 

 

u =W0 ( p)(u1 + u2 ) ,

(2.4)

где

W0

( p) =

b0 ( p)

- дробно – рациональная функция. На первом этапе определим

a0 ( p)

 

 

 

 

 

алгоритм формирования переменной u1 для того, чтобы решить задачу стабилиза-

ции, т.е. обеспечить выполнение предельного соотношения lim y(t) = 0 при любых

t→∞

начальных условиях. На втором этапе, выбором сигнала управления обеспечим ре-

шение задачи слежения, т.е. выполнение неравенства | y(t) y (t) |≤ γ, где y(t) -

функция времени такая, что

lim[ y(t) y(t)] = 0 ,

t→∞

y (t) - заранее неизвестная переменная величина (задающее воздействие), а γ > 0 -

некоторое заданное постоянное число. Отметим, что определенная подобным обра-

зом функция y(t) называется установившейся реакцией системы.

В заключении первого этапа курсовой работы строятся ЛАХ и ЛФХ передаточной функции W (s) объекта управления. Напомним, что логарифмической амплитудно-

частотной характеристикой, соответствующей передаточной функции W (s) , назы-

вается график функции

L(ω) = 20lg | W ( jω) |

12

от логарифма lgω, где ω - действительная переменная, которая называется часто-

той. Функция W ( jω) , которую получают из передаточной функции W (s) при под-

становке в нее s = jω, называется частотной передаточной функцией. Принято измерять значение ЛАХ в децибелах (дБ), а значение наклона ЛАХ – в децибелах на декаду. Декадой называют интервал на котором частота ω меняется в 10 раз. Логарифмической фазочастотной характеристикой называется график функции

ϕ(ω) = argW ( jω)

от логарифма lgω.

2.2. Решение задачи стабилизации

Для решения задачи стабилизации выберем алгоритм формирования сигнала u1 в

виде

 

 

 

 

 

 

 

u1 = −W1( p) y ,

 

 

(2.5)

где W

( p) =

b1( p)

. Подставляя (2.4), (2.5) в уравнение (2.2) и разрешая его относи-

 

 

1

 

a1( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тельно выходной переменной найдем уравнение замкнутой системы

 

 

 

 

 

 

 

 

 

y =Wy ( p)u2 ,

 

 

(2.6)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wy ( p) =

 

 

 

W ( p)W0 ( p)

 

=

 

b( p)b0 ( p)a1( p)

 

.

(2.7)

 

1

+W ( p)W0 ( p)W1

( p)

a( p)a0

( p)a1( p) +b( p)b0

( p)b1( p)

 

 

 

 

 

 

Для объектов управления с минимально-фазовыми передаточными функциями, по-

линомы a0 ( p), a1( p), b0 ( p), b1( p)

выбирают из условия

 

 

 

 

 

 

b( p)b0 ( p)a1( p)

 

=

 

 

1

,

(2.8)

 

a( p)a ( p)a ( p) +b( p)b ( p)b ( p)

a

y

( p)

0

1

0

1

 

 

 

 

 

где ay ( p) - произвольный устойчивый полином степени n m . Например, полином ay ( p) можно выбрать следующим образом

ay ( p) = ( p +1)nm ,

(2.9)

13

а многочлены a0 ( p), a1( p), b0 ( p), b1( p) -

a0 ( p) = a0,n1 pn1 + a0,n2 pn2 +K+ a0,0 ,

a1( p) = a1,n1 pn1 + a1,n2 pn2 +K+ a1,0 ,

b0 ( p) = b0,n1 pn1 +b0,n2 pn2 +K+b0,0 ,

b1( p) = b1,n1 pn1 +b1,n2 pn2 +K+b1,0 ,

где 4(n 1) неизвестных коэффициентов этих полиномов находятся из тождеств

 

 

 

 

a ( p) = b ( p) = ( p +1)n1 ,

 

 

 

 

 

 

(2.10)

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a( p)a ( p) +b( p)b ( p) = b( p)( p +1)2nm1 .

 

 

 

 

(2.11)

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, принимая во внимание

u2

 

 

 

 

 

 

 

 

 

 

 

y

(2.7), (2.9)-(2.11), передаточная функ-

 

 

 

 

 

W0 ( p)

 

 

W ( p)

 

 

 

 

 

 

 

 

 

 

 

ция объекта

управления

охваченного

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

отрицательной обратной связью по вы-

 

 

 

 

 

 

 

 

 

 

W1 ( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ходу будет иметь вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

1

 

 

 

 

 

 

Wy (s) =

 

 

, (2.12)

 

 

 

 

 

 

системы управления

 

 

(s

+1)nm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Структурная схема системы управления представлена на рис.2.2.

В заключении данного этапа строятся ЛАХ и ЛФХ передаточной функции Wy (s) .

2.3. Синтез следящей системы управления

С учетом того, что измерения проводятся с помехами, управляющее воздействие объектом (2.6) строится в форме:

 

 

 

 

 

 

u2

~

(2.13)

 

 

 

 

 

 

=W2 ( p)ε,

 

W2

( p) =

b2

( p)

 

~

 

 

где

 

 

и

ε = y * (t)

y f (t) = ε − f (t) . Структурная

схема системы

a2

( p)

 

 

 

 

 

 

 

(2.2), (2.4), (2.5), (2.13) приведена на рис.2.3.

14

Задача этого этапа со-

 

u2

 

W0 ( p)

u

W ( p)

y

y *

стоит в выборе такой переда-

 

 

 

точной

 

функции

регулятора

 

 

 

u1

 

 

 

 

ε

 

 

 

 

 

 

 

 

 

 

W2 (s) , которая обеспечит бли-

 

 

 

 

 

 

W1 ( p)

 

 

зость

ошибки

управления

 

 

 

 

 

 

 

~

f

ε(t) = y

 

(t) y(t)

к

нулю

и

 

 

 

 

 

 

W2 ( p)

ε

 

 

 

 

 

 

 

 

 

 

 

требуемые

характеристики

 

 

 

Рис. 2.3. Структурная схема сле-

 

замкнутой

системы

управле-

 

 

 

дящей системы управления

 

 

 

 

 

 

 

 

 

 

ния. Для решения этой задачи воспользуемся методом динамической компенсации.

Используя уравнения (2.2), (2.4), (2.5) и (2.13) найдем выражение для ошибки

управления

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε(t) = [1 W

( p)]y (t) +W

( p) f (t) ,

 

 

 

(2.14)

 

 

 

 

 

 

 

 

з

 

 

 

з

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wз( p) =

 

Wp ( p)

 

,

Wp ( p) =Wy ( p)W2 ( p) .

 

(2.15)

 

 

 

1

+Wp

( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция Wз(s) называется передаточной функцией замкнутой системы, а функция

Wp (s) - передаточной функцией разомкнутой системы.

 

 

 

Выбор передаточной функции регулятора определяет вид передаточной функции

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

том, чтобы предъявить требования к самой Wз(s) , к примеру, потребовать выполне-

ние условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (s)

 

 

 

 

 

 

 

 

 

W (s) =W (s) =

p

,

(2.16)

 

 

 

 

 

 

 

1 +Wp (s)

 

 

 

 

 

 

 

з

 

з

 

 

 

 

 

b

(s)

 

 

 

bp (s)

 

 

 

 

где W

з

(s) =

з

 

и W

p

(s) =

 

- соответственно, желаемые передаточные функ-

aз (s)

ap (s)

 

 

 

 

 

 

 

 

ции замкнутой и разомкнутой системы, которые формируются на основе требований

ксистеме управления. Преобразуем требования, предъявленные в курсовой работе,

кзамкнутой системе в условия на передаточную функцию замкнутой системы.

15

Требование к устойчивости замкнутой системы (a) означает, что все корни харак-

теристического полинома

 

(s) =[a(s)a0 (s)a1 (s) + b(s)b0 (s)b1 (s)]a2 (s) + a1 (s)b(s)b0 (s)b2 (s)

(2.17)

замкнутой системы (2.2), (2.4), (2.5), (2.13) должны иметь строго отрицательную вещественную часть.

Установившаяся ошибка, вызванная воздействием (1.1), равна

 

 

 

 

 

 

d

 

 

 

ε

(t) = [1 W

(0)](a

+ a t) +

[1 W (s)]

a .

(2.18)

 

 

з

0

1

ds

з

1

 

 

 

 

 

 

 

s=0

 

 

 

 

 

 

 

 

Ошибка может быть ограниченной при a1

0 , только если выполнено условие ас-

татизма

 

 

 

 

 

 

 

 

 

 

 

Wз (0) =1

 

 

(2.19)

или

 

 

 

 

 

 

 

 

 

 

 

Wp1(0) = 0 ,

 

 

(2.20)

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

Wp (s) =

kp

~

~

 

 

 

Wp (s) , Wp (0)

=1,

(2.21)

s

 

 

 

 

 

где kp - коэффициент усиления разомкнутого контура. Тогда используя выражения

(2.15) и (2.21), найдем

 

 

 

 

 

 

d

 

 

= k1

 

 

[1 W (s)]

 

(2.22)

 

 

 

 

ds

 

з

s=0

p

 

 

 

 

 

 

и, следовательно,

 

 

 

 

 

 

 

ε

(t) = k1a .

(2.23)

 

 

 

 

p

1

 

Таким образом, для выполнения требования (1.2) необходимо, чтобы коэффициент удовлетворял условию

k

p

a1

= δ1.

(2.24)

ε

 

 

1

 

 

 

1

 

 

Установившаяся ошибка, вызванная воздействием (1.3), является гармонической функцией с амплитудой

16

εy =

 

1 Wз ( jω)

 

a y ,

(2.25)

 

 

где ω- частота воздействия. Таким образом, для удовлетворения требования (1.4) необходимо, чтобы

 

1 W

( jω)

 

εy

= δ

y

, ω ≤ ω .

(2.26)

 

 

 

 

 

 

з

 

 

 

a y

 

y

 

 

 

 

 

 

 

 

 

 

Установившаяся ошибка вызванная воздействием (1.5), является гармонической функцией с амплитудой

 

ε f =

 

Wз ( jω)

 

a f .

(2.27)

 

 

Таким образом, для удовлетворения требования (1.7) необходимо, чтобы

 

 

Wз ( jω)

 

εf

= δ f , ω ≥ ωf .

(2.28)

 

 

 

 

 

 

 

 

 

 

 

a f

 

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

Явные зависимости L(ω) , ϕ(ω) достаточно сложны. Поэтому часто ограничива-

ются построением асимптотических логарифмических частотных характеристик. Для примера рассмотрим произвольную минимально-фазовую передаточную функцию с вещественными нулями и полюсами

 

 

 

 

 

−ν

(Tis +1)

 

 

W (s) = ks

i

,

(2.29)

 

′′

 

 

(Tis +1)

 

 

i

где k - положительный коэффициент усиления, Ti, Ti′′ - положительные постоян-

ные времени,. ν - натуральное число (порядок астатизма). Тогда соответствующие ЛАХ и ЛФХ даются формулами

2

′′

2

+1]

,

(2.30)

L(ω) = 20lg k 20νlgω+ 20lg[(Tiω)

 

+1] 20lg[(Ti ω)

 

i

 

i

 

 

 

 

17

ϕ(ω) = ν

π

′′

(2.31)

2

+ arctg(Tiω) arctg(Ti ω) .

 

i

i

 

Асимптотическая ЛАХ есть кусочно-линейная функция, получаемая заменой в

(29) членов 20 lg[(Tω) 2 +1] на

0, ω ≤ T1 ,

20lgTω, ω ≥ 1 ,

T

где частота ω = T1 называется сопрягающей.

Асимптотическая ЛФХ есть кусочно-постоянная функция, получаемая заменой в

(30) членов arctg(Tω) на

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

0,

ω ≤

 

 

 

,

 

T

 

 

 

 

 

 

 

π

 

 

 

1

 

 

 

, ω ≥

 

.

 

 

T

2

 

 

 

Нетрудно заметить, что для систем с передаточной функцией вида (2.29) по асимптотической ЛАХ можно восстановить асимптотическую ЛФХ и саму передаточную функцию. Данное свойство характерно для любой передаточной функции, не имеющей нулей и полюсов в правой полуплоскости.

Преобразуем требования к замкнутой системе в ограничения на свойства логарифмической амплитудной характеристики разомкнутой системы. Условие (2.21) на отработку линейно растущего задающего воздействия сводятся к ограничению на

поведение Wp ( jω) при низких частотах, близких ω = 0 :

 

Wp ( jω)

kp

 

(2.32)

jω

 

 

 

 

 

 

и, следовательно,

 

Lp (ω) = 20lg

 

Wp ( jω)

 

= 20lg kp 20lg ω.

(2.33)

 

 

Иначе говоря, низкочастотная асимптота Lp (ω) должна иметь наклон -20 дБ/дек,

причем в силу (2.24) ее уровень определяется условием

18

 

20 lg k p ≥ −20 lg δ1 .

(2.34)

Условия (2.26), (2.28) при достаточно малых δy , δ f можно заменить на

 

 

W

p

( jω)

 

 

≥ δ1 , ω ≤ ω ,

(2.35)

 

 

 

 

 

 

 

 

 

 

 

y

y

 

 

Wp ( jω)

 

≤ δf

, ω ≥ ωf

(2.36)

 

 

и, следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

Lp (ω) = 20lg

 

Wp ( jω)

 

≥ −20lg δy , ω ≤ ωy ,

(2.37)

 

 

Lp (ω) = 20lg

 

Wp ( jω)

 

20lg δf , ω ≥ ωf .

(2.38)

 

 

Условие (2.37) задает ограничение на поведение логарифмической амплитудной частотной характеристики в области низких частот, а условие (2.38) - в высокочастотной области.

Поведение ЛАХ в области средних частот определяет запасы устойчивости по фазе и амплитуде и в значительной мере качество системы в переходном режиме (время переходного процесса и перерегулирование). Для обеспечения приемлемых запасов устойчивости наклон ЛАХ на частоте ωc такой, что Lp (ωc ) = 0 обычно выбира-

ется равным -20 дБ/дек, причем длительность этого участка должна быть не менее декады, что соответствует изменению частоты в 10 раз. Частота ωc называется час-

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

Итак, все требования, которые были сформулированы в задании курсовой роботы, сведены к ограничениям на допустимое поведение ЛАХ разомкнутой системы (см. рис. 2.4, где заштрихованы границы зон, в которые не может заходить ЛАХ). Оста-

ется подобрать передаточную функцию разомкнутого контура Wp (s) , для которой эти ограничения выполнены. Если Wp (s) минимально-фазовая передаточная функ-

ция, то вначале строят асимптотическую ЛАХ, удовлетворяющую всем ограничени-

ем, а затем по ней находят саму Wp (s) .

Приведем в готовом виде сводку допустимых передаточных функций разомкнутой системы:

19

Lp (ω)

 

 

 

20 lg k p

-20 дБ/дек

 

 

 

 

 

20lg δy

 

 

 

 

-20 дБ/дек

 

lg δ1

20 lg δf

lg ωy lg ωc

lg ωf

lg ω

 

 

 

Рис. 2.4. Ограничения на допустимое поведение ЛАХ

1) Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

≤ δ1

δy

.

(2.39)

 

δ

 

ω

 

ω

 

 

f

f

 

y

 

 

 

 

 

 

 

 

 

Тогда передаточная функция

 

 

 

 

 

 

 

 

 

 

 

Wp (s) =

kp

, kp = δf ωf

(2.40)

 

 

 

 

 

 

 

s

 

 

 

 

 

удовлетворяет всем ограничениям. При тех же условиях допустима, но требует меньшего усиления на частотах ω ≥ ωy и обеспечивает лучшее подавление помех передаточная функция вида

 

 

 

Wp (s) =

 

 

 

 

kp (T2s +1)

 

 

 

,

 

(2.41)

 

 

 

 

s(T1s +1)(T3s +1)

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δy

 

 

 

 

 

 

 

 

 

 

 

 

 

k

p

= δ1

, T =

 

 

 

 

 

 

, T =

 

 

 

µδy

 

 

, T 0.1T ,

(2.42)

 

1

1

 

δ1 ωy2

2

 

 

ωy

 

 

 

 

3

2

 

а параметр µ, характеризующий расположение частоты среза ωc

= µ / T2 , может вы-

бираться в пределах 2÷4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ1

 

δy

 

, δ1 <

 

 

 

 

1

 

 

 

.

 

 

(2.43)

 

 

 

 

ω

y

δ

f

ω

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20