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

vostrikov

.pdf
Скачиваний:
133
Добавлен:
23.01.2019
Размер:
12.99 Mб
Скачать

6.4. Частотный метод синтеза

 

 

 

 

 

 

 

 

191

 

 

 

 

 

L, дБ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( °)

 

 

 

 

 

 

 

 

 

 

 

 

L0

 

 

 

 

 

 

 

 

 

 

 

 

 

L0 (ω)

 

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

LK (

)

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2

L

 

 

+40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

lg

, дБ

 

 

 

 

 

 

 

 

 

 

 

L*(

)

 

 

 

 

Рис. 6.14. Логарифмические амплитудно-частотные характеристики

 

 

 

 

 

 

к примеру 6.7

 

 

 

 

 

σ, %

 

 

 

tn

 

LD,LдБ,

 

 

 

 

 

 

s,

%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

80

 

 

tn

 

 

 

25

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

%

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

n

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

n

10

 

 

 

 

 

 

20

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Pm ax

0

 

 

 

 

 

%

 

 

 

 

 

 

 

 

 

 

σ, , дБ

 

1,0

1,2

1,4

1,6

1,8

 

20

25

30

35

40

45

50

 

 

Рис. 6.15. Номограммы для определения параметров желаемой ЛАЧХ

192 Глава 6. СИНТЕЗ ЛИНЕЙНЫХ СИСТЕМ

Определим теперь

 

 

) и запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lк (

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22 p2

 

 

2 p 1

 

 

 

3 p 1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wк

( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 p

 

 

 

1 42 p2

 

 

 

4 p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где lg

1

1;

 

 

lg

1

 

0;

 

lg

1

 

 

 

0, 7;

lg

1

 

1,5 .

 

Отсюда

следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

1

10 с,

2

 

1 с,

3

 

0, 2 с,

4

 

0, 03 с . Окончательно запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p2

 

 

p

1)(0, 2 p 1)

 

 

,

 

 

 

Wк ( p) kкWк ( p) 50

(10 p 1)(0, 0009 p2

 

0, 03 p

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wк ( p) 5556

 

 

 

0, 2 p3

1, 2 p2

1, 2 p 1

.

 

 

 

 

 

 

 

 

p3

 

 

 

33, 4 p2

1114, 4 p

111,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

5556

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

p

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

33,4

1114,4

111,1

Рис. 6.16. Пример схемной реализации регулятора

Схемная реализация полученной передаточной функции корректирующего звена, соответствующая второму каноническому представлению (см. подразд. 3.6.2), показана на рис. 6.16.

6.5. Модальный метод синтеза

193

6.5. МОДАЛЬНЫЙ МЕТОД СИНТЕЗА

6.5.1. ОСНОВНЫЕ ПОНЯТИЯ

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

Рассмотрим основные соотношения метода для случая, когда математическая модель объекта управления представлена в переменных состояния

x Ax Bu M (t),

x Rn ,

(6.41)

 

u Rm , y Rm .

y Cx,

 

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

оценок переходных процессов (6.5): t

t * и

* , от которых мож-

п

п

 

но перейти к желаемому распределению корней на комплексной плоскости. На основе выбранных корней формируется желаемое (эталонное) характеристическое уравнение замкнутой системы

pn c pn 1

c p

c 0 .

(6.42)

n

2

1

 

Метод синтеза предполагает организацию «пропорционального»

закона управления

 

u Kx,

(6.43)

где K – матрица назначаемых коэффициентов.

После подстановки алгоритма управления (6.43) в уравнения объекта (6.41) получают уравнения замкнутой системы

x ( A BK ) x M (t),

(6.44)

y Cx

 

 

 

и записывают ее характеристическое уравнение

 

det[ pI (A BK)] pn a (K) pn 1

a (K) p

a (K) 0. (6.45)

n

2

1

194

Глава 6. СИНТЕЗ ЛИНЕЙНЫХ СИСТЕМ

Неизвестные коэффициенты матрицы регулятора необходимо определить таким образом, чтобы качество работы замкнутой системы соответствовало заданным оценкам. С этой целью приравнивают характеристическое уравнение замкнутой системы (6.45) желаемому (6.42) и получают соотношения для расчета элементов матрицы K в виде

 

 

 

 

 

ai (K ) ci ,

i 1, n .

(6.46)

Поскольку в общем случае зависимость ai (K) может быть нели-

нейной, найти коэффициенты матрицы K по выражению (6.46) не всегда удается даже для одноканального объекта, уравнения которого предварительно записывают в канонической форме.

Часто линейный одноканальный объект удобнее описывать с помощью передаточной функции, поэтому обсудим далее операторную методику модального метода синтеза, предложенную на кафедре автоматики НГТУ.

6.5.2. ПОСТАНОВКА ЗАДАЧИ СИНТЕЗА ДЛЯ ОДНОКАНАЛЬНОГО ОБЪЕКТА

Рассмотрим объект управления, поведение которого описывает передаточная функция

 

B( p)

 

b

pm

b pm 1

b

 

W

 

 

m 1

 

m

1

,

(6.47)

 

 

 

 

 

 

0

A( p)

 

pn

an pn 1 a1

 

 

 

 

где m n .

Модальный метод синтеза обеспечивает заданную реакцию системы на начальные условия, которая определяется корнями характеристического уравнения. Корни в свою очередь выбираются на основе требований, предъявляемых к динамике в виде условий (6.5). Кроме этого необходимо, чтобы в статике выполнялось условие (6.4), т. е.

lim y(t) v с точностью 0

*0 .

t

 

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

из двух составляющих: последовательного звена Ws ( p) на входе и звена с передаточной функцией Wd ( p) в цепи локальной обратной

6.5. Модальный метод синтеза

195

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

Звено прямого канала с передаточной функцией Ws ( p) будем называть корректором статики, а звено с передаточной функцией Wd ( p) –

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

 

 

 

 

 

 

 

 

M

v

 

u

 

 

 

y

 

 

Ws ( p)

 

W0

( p)

 

 

 

 

 

 

 

 

Wd ( p)

Рис. 6.17. Расчетная структурная схема для модального метода синтеза

Рассмотрим последовательно этапы модального метода синтеза системы регулирования.

6.5.3. ВЫБОР КОРРЕКТОРА СТАТИКИ

Для обеспечения условия статики (6.4) при произвольном возмуще-

нии M (t) , т. е. выполнения свойства

lim y(t)

v , предлагается в каче-

 

t

 

 

стве корректирующего звена Ws ( p) использовать интегратор

W ( p)

ks

,

(6.48)

 

s

p

 

 

 

 

 

где ks – коэффициент усиления регулятора; его численное значение

будет определено позже.

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

y

pM W0 ( p) ksv

 

p pW0 ( p) Wd ( p) W0 ( p) ks .

(6.49)

196

Глава 6. СИНТЕЗ ЛИНЕЙНЫХ СИСТЕМ

Поскольку в статике передаточные функции W0 ( p) и Wd ( p)

«вы-

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

v .

Таким образом, использование корректора статики Ws ( p) вида (6.48) делает систему астатической, и условие (6.4) можно обеспечить

сошибкой 0 0 .

6.5.4.РАСЧЕТ КОРРЕКТОРА ДИНАМИКИ

Вкачестве корректора динамики предлагается выбирать звено со следующей передаточной функцией:

 

D ( p)

 

d

n

pn 1 d

 

 

Wd ( p)

 

 

 

 

1

,

(6.50)

B ( p)

 

bm 1 pm

bm pm 1

b1

 

 

 

 

где B( p) – полином числителя передаточной функции объекта W0 ( p) , а D( p) – введенный расчетный полином с неизвестными коэффициен-

тами di , i 1, n .

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

Первоначально определим характеристическое уравнение системы, структурная схема которой приведена на рис. 6.17:

1 W0 ( p)Wd ( p) W0 ( p)Ws ( p) 0 .

(6.51)

С учетом (6.47), (6.48) и (6.50) уравнение (6.51) принимает вид

 

pA( p) pD( p) ks B( p) 0 ,

 

причем его порядок равен ( n 1 ).

Подставляя вместо A( p) , D( p) и B( p) их выражения, получим

действительное характеристическое уравнение замкнутой системы в следующей форме:

pn 1 (a

d

n

) pn

(a

d

k b ) p

k b

0 . (6.52)

n

 

 

1

1

s 2

s 1

 

Теперь на основе требований к качеству переходных

процессов

(заданного перерегулирования

и быстродействия tn* ) сформируем

6.5. Модальный метод синтеза

197

желаемое характеристическое уравнение того же порядка, что и (6.52). Для его конструирования используем корневые оценки переходных процессов, с помощью которых получим эталонное распределение корней на комплексной плоскости (см. подразд. 5.5.2).

Предварительно определим границу расположения заданных корней системы. Она зависит от заданного времени переходного процесса tn* и приближенно может быть найдена по соотношению

 

3

.

(6.53)

 

 

 

t n*

 

Заданное перерегулирование

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

ной плоскости, внутри которого должны располагаться заданные корни (рис. 6.18). С этой целью по соотношению

* 100exp( / * )

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

ным» размахом:

 

 

 

* .

Эталонные корни

*, ,

*

могут выби-

 

1

n 1

 

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

ходного процесса и больше потребуется ресурс управления объекта. Поэтому рекоменду-

ется выбирать корни i , i 1, n , достаточно

близкие друг к другу и правой границе области расположения корней, а затем сформировать желаемое уравнение следующим образом:

(6.54)

ImI

Re

Re

Риc. 6.18. К определению области расположения корней

C( p) ( p

*) ( p

*

) 0 .

(6.55)

 

1

n 1

 

 

198 Глава 6. СИНТЕЗ ЛИНЕЙНЫХ СИСТЕМ

Характеристическое уравнение (6.55) запишем в стандартном виде

C( p) pn 1 c

pn c 0 .

(6.56)

n 1

1

 

Приравнивая коэффициенты при соответствующих степенях оператора p желаемого (6.56) и действительного (6.52) характеристических уравнений системы, запишем выражения для определения неизвестных параметров регулятора:

c1

ksb1,

 

 

 

c2

a1

d1

b2ks ,

(6.57)

 

 

cn 1

an

 

dn .

 

Полученные из (6.57) расчетные соотношения имеют вид

 

 

c1

 

 

 

 

 

 

 

(6.58)

k

s

, d

i

c

a

k b , i 1, n .

 

 

b1

i 1

i

s i

 

 

 

 

 

 

 

 

 

 

Таким образом, мы определили параметры передаточных функций Ws ( p) и Wd ( p) регулятора, обеспечивающего в системе требуемые свойства в статике и динамике.

ПРИМЕР 6.8

Поведение одноканального объекта описывает передаточная функция

5

W0 ( p) p2 3 p 1 .

Требуется синтезировать систему, в которой качество процессов будет

отвечать следующим требованиям: tn 3 с;

0,

0

0.

 

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

В качестве корректора статики используем интегрирующее звено с передаточной функцией Ws ( p) k p , что гарантирует нулевую статическую

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

6.5. Модальный метод синтеза

 

 

199

W ( p)

d1 p

d0

.

 

 

d

5

 

 

 

 

 

 

Здесь k, d1, d0 – неизвестные коэффициенты регулятора, которые требу-

ется определить.

 

 

 

Используя структурные преобразования, запишем характеристическое

уравнение замкнутой системы (см. рис. 6.17)

 

A( p) p3 (3 d ) p2

(d

0

1) p 5k 0.

1

 

 

 

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

расстоянии не ближе

 

3 t n*

1 от мнимой оси. В результате выберем

следующие корни:

 

 

 

 

 

 

 

*

2,

*

2, 5,

*

3 .

 

1

2

3

 

 

 

 

В соответствии с (6.55) получим желаемое характеристическое уравнение

 

C( p)

p3 7, 5 p2

18, 5 p 15

0.

Запишем расчетные соотношения (6.57):

 

 

 

 

 

 

3

d1

7,5,

 

 

 

 

 

 

 

 

d0

1

18,5,

 

 

 

 

 

 

 

5k

15.

 

 

 

 

 

Отсюда находим параметры d1

4,5,

d0

19,5,

k

3 . Следовательно, пе-

редаточные функции регулятора имеют вид

 

 

W ( p)

3

,

W ( p)

 

4,5 p

19,5

 

0,9 p 3,9 .

 

 

 

 

 

 

s

 

p

d

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

6.5.5. РЕАЛИЗАЦИЯ РЕГУЛЯТОРА

Рассмотрим возможность реализации регулятора, рассчитанного модальным методом. Корректор статики с передаточной функцией Ws ( p) , представляющий собой обычный интегратор, не вызывает за-

труднений. Остановимся подробнее на реализации звена обратной связи с передаточной функцией Wd ( p) .

200

Глава 6. СИНТЕЗ ЛИНЕЙНЫХ СИСТЕМ

Поскольку для реальных объектов управления степень полинома числителя передаточной функции W0 ( p) обычно меньше степени полинома ее знаменателя ( m n ), корректор динамики

 

D ( p)

 

d

n

pn 1 d

 

Wd ( p)

 

 

 

 

1

,

B ( p)

 

bm 1 pm

bm pm 1

b1

 

 

 

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

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

состоит из модели Wm ( p) (с выходом yˆ ) и стабилизирующей добавки L( p) (рис. 6.19). Его называют фильтром Калмана–Бьюсси или параллельным фильтром.

u

 

 

 

yy

 

 

W00((pр)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L(рp)

В(р)

А(р) y

Рис. 6.19. Схема подключения фильтра

Здесь передаточная функция параллельной модели Wm ( p) W0 ( p) . Назначение стабилизирующей добавки L( p) – «сводить» к нулю разницу между выходом объекта у и выходом модели yˆ .

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

B( p)

u

B( p)

[u L( p) ] ,

 

 

A( p)

A( p)

 

которое после преобразований принимает вид

A( p) B( p)L( p) 0 .

Соседние файлы в предмете Основы Теории Управления