
vostrikov
.pdf
11.3. Метод локализации |
381 |
Отметим, что все нелинейности и действующие на объект внешние возмущения, отраженные в функциях f (t, y) и b(t, y) , подавляются
большим коэффициентом k , численные значения которого рекомендуется выбирать по соотношению
bmin |
k (20 100) . |
(11.24) |
В этом случае точность обеспечения требуемых свойств (11.23) будет составлять
(0,05 0,01) f (t, y) ,
что соответствует обычной в практике 5%-ной зоне ошибок.
Полученный эффект объясняется структурной локализацией возмущений, что иллюстрирует схема системы, показанная на рис. 11.3.
|
|
|
|
|
|
t |
|||
|
|
|
|
|
|
f |
|||
|
|
t |
y |
|
|
y(0) |
|||
v |
|
u |
|
|
|
|
|
|
y |
F |
|
k |
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 11.3. Структурная иллюстрация метода
Как видим, система двухконтурная: внешний контур образован обычной обратной связью по выходной переменной, а внутренний – обратной связью по производной. Именно в нем локализовано влияние функций f (t, y) и b(t, y) , которое подавляется коэффициентом уси-
ления k .
Отметим, что внутренний контур является безынерционным, так как не содержит инерционных звеньев.
Для реализации закона (11.21) необходимо убедиться в том, что численные значения управления в любой момент времени не будут превышать имеющегося ресурса объекта. С этой целью исследуем управляющее воздействие в замкнутой системе, подставляя в (11.21) вместо y правую часть уравнения (11.19):
u k F ( ) f ( ) b( )u .

382 |
|
Глава 11. СИНТЕЗ НЕЛИНЕЙНЫХ СИСТЕМ |
||
После несложных преобразований получим |
|
|||
u |
k |
F ( ) |
f ( ) . |
|
|
|
|||
1 b( )k |
|
|||
В асимптотике при k |
управляющее воздействие в замкнутой |
|||
системе принимает вид |
|
|
|
|
u |
b 1( ) |
F ( ) |
f ( ) . |
(11.25) |
Отсюда следуют несколько важных выводов.
1.Асимптотическое управление (11.25) совпадает с управляющим воздействием (11.16), соответствующим точному решению поставленной задачи синтеза. Следовательно, основанный на методе локализации закон (11.21) представляет собой неявную реализацию «точного» управления (11.16).
2.Управляющее воздействие в замкнутой системе остается конечным даже при бесконечном коэффициенте усиления.
3.Выражение (11.25) позволяет определить максимальное значение управления в наихудшей ситуации, когда все функции достигают своих предельных значений
u |
max |
b 1 |
|
F |
|
f |
max |
|
. |
|
|
|
|||||||
|
min |
|
max |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
Больше этого значения управление в замкнутой системе быть не может. Если выполняется условие
u |
|
b 1 |
|
F |
|
f |
|
|
|
|
|
|
|
max |
|
|
max |
|
U |
|
, |
(11.26) |
|||||
|
min |
|
max |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то в ней можно обеспечить желаемые процессы (11.20).
11.3.2. ОЦЕНКА ВЛИЯНИЯ ПОМЕХ ИЗМЕРЕНИЯ
Выходная переменная физической системы измеряется датчиком, который всегда имеет некоторую помеху измерения. Как известно, дифференцирование усиливает влияние высокочастотных помех, поэтому в системах с законом управления (11.21) важно оценить их влияние.

11.3. Метод локализации |
|
|
383 |
Полагаем, что выходная переменная датчика имеет вид |
|
||
y |
y |
h(t) , |
(11.27) |
где h(t) – помеха измерения. |
|
|
|
Закон управления (11.21) в этом случае принимает форму |
|
||
u k |
F( y |
|
(11.28) |
,v) y . |
С учетом (11.27) и (11.19) запишем управление в замкнутой системе
u k F ( ) f ( ) b( )u h
или после преобразований
uk F( ) f ( ) h . (11.29)
1b( )k
Впределе при k асимптотическое управление будет следующим:
u b 1( ) |
F( ) |
|
f ( ) h , |
(11.30) |
где слагаемое |
|
|
|
|
u |
b |
1 |
|
(11.31) |
|
( )h |
|||
h |
|
|
|
|
определяет «вклад» помехи в управление. Сравнение (11.25) и (11.30) показывает, что при наличии помехи измерения может требоваться больший ресурс управления для объекта, чем без нее.
Выражение (11.30) следует иметь в виду при выборе датчика.
11.3.3. ДИФФЕРЕНЦИРУЮЩИЙ ФИЛЬТР
При практической реализации закона управления (11.21) необходимо оценивать производную y . С этой целью в систему следует доба-
вить специальное устройство, реализованное на интеграторах, которое будем называть дифференцирующим фильтром [8]. Он представляет собой динамическую систему первого порядка
|
yˆ y |
(11.32) |
yˆ |

384 |
Глава 11. СИНТЕЗ НЕЛИНЕЙНЫХ СИСТЕМ |
||
или второго |
|
|
|
2 |
|
yˆ y , |
(11.33) |
yˆ |
2d yˆ |
выбираемую в зависимости от частотного состава и уровня помех измерения. Здесь yˆ – измеренное значение выходной переменной; –
малый параметр, отражающий инерционность фильтра.
Структурная схема дифференцирующего фильтра второго порядка показана на рис. 11.4.
yˆ x
y |
|
|
|
|
yˆ |
x |
|
1 |
|
1 |
x |
|
|
a
2d2
a3
Рис. 11.4. Структурная интерпретация фильтра второго порядка
Поскольку дифференцирующий фильтр является линейным звеном, можно записать его передаточную функцию
Wф ( p) |
yˆ |
1 |
, |
(11.34) |
|
|
|
|
|||
y |
|
D( p) |
|||
|
|
|
|
где D( p) – характеристический полином фильтра, его будем назы-
вать также «фильтрующим» полиномом.
Покажем, что это устройство действительно позволяет получить оценку производной, выражение для которой с учетом (11.34) представим в виде
|
p |
|
yˆ |
|
y . |
D( p) |
Так как p есть оператор дифференцирования, то py y и
|
1 |
y . |
yˆ |
|
|
D( p) |

11.3. Метод локализации |
|
385 |
|
При |
0 получим |
|
y , т. е. оценка производной совпадает с ее |
yˆ |
точным значением. Таким образом, для реализации закона управления (11.21) необходимо использовать дифференцирующий фильтр с малой инерционностью. На практике достаточно, чтобы инерционность фильтра была на порядок меньше, чем основные процессы в системе.
Заметим, что наличие таких звеньев приводит к возникновению в системе разнотемповых процессов, причем для ее работоспособности быстрые движения должны быть устойчивыми.
11.3.4.АНАЛИЗ ВЛИЯНИЯ МАЛЫХ ИНЕРЦИОННОСТЕЙ
Спомощью метода разделения движений выделим медленные и быстрые процессы в системе с дифференцирующим фильтром (рис. 11.5).
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
t |
|
xy |
|
|
|
|
|
|
|
|||
v |
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
x |
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
||||
F |
|
|
|
|
k |
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
W |
(p) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
ф( p) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ф |
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 11.5. Структурная схема системы с фильтром
При использовании фильтра первого порядка уравнения замкнутой системы имеют вид
y f ( ) b( )k F ( ) yˆ ,
(11.35)
yˆ y yˆ.
Поскольку первое уравнение содержит производную в правой части, необходимо сначала представить описание в стандартной форме. С этой целью введем новую переменную
z 1( y yˆ)

386 |
Глава 11. СИНТЕЗ НЕЛИНЕЙНЫХ СИСТЕМ |
и получим следующую математическую модель системы:
y f ( ) b( )k F ( ) z ,
(11.36)
z f ( ) b( )k F ( ) z z.
Выделим подсистему быстрых движений
y const, |
z f ( ) b( )k F ( ) z z. |
Запишем характеристическое уравнение подсистемы быстрых движе-
ний (рис. 11.6)
|
|
|
|
|
|
|
|
p |
bk |
1 0 . |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
|||||
|
F |
|
|
|
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
|
|
y |
|
|
F |
|
|
|
u |
|
|
|
|
|
|
|
|
|
|
x& |
||||||
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yˆ
€
x& WWфф((pp))
Рис. 11.6. Структурная схема контура быстрых движений
В общем случае с учетом D( p) оно принимает вид
D( p) bk 0 . |
(11.37) |
Как видим, контур быстрых движений линейный, поэтому для анализа устойчивости можно применять известные критерии устойчивости (см. главу 4). Если используется фильтр первого или второго порядка, то этот контур будет устойчив при любых положительных значениях b( ) .
Выделим теперь медленные движения, полагая в (11.36) |
0 , |
|
y |
f ( ) b( )k F ( ) z , f ( ) b( )k F ( ) z z. |
|
Так как f ( ) |
b( )k F ( ) z y , то уравнения подсистемы медленных |
движений можно записать в форме


388 |
|
|
|
|
|
|
|
Глава 11. СИНТЕЗ НЕЛИНЕЙНЫХ СИСТЕМ |
|||||||||||||||
|
|
|
|
|
|
|
|
|
f |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
F |
|
|
|
k |
u |
b |
|
|
|
|
|
|
|
|
|
|
y(yn) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H(n) |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1
D( p)
Рис. 11.8. Контур быстрых движений с учетом помехи измерения
Так как этот контур линейный и для него справедлив принцип суперпозиции, то определим реакцию u на воздействие h(t) :
u |
Wф ( p)k |
h(t) , |
|
|
||
|
|
|
|
|
||
1 Wф |
( p)kb |
|
|
|
||
которую можно представить в виде |
|
|
|
|
|
|
u Wк ( p)h |
|
|
kp |
h , |
(11.39) |
|
|
|
|
||||
|
D( |
p) kb |
||||
|
|
|
|
|
где Wк ( p) – передаточная функция контура.
|
Рассмотрим теперь влияние полинома |
||||||
L, дБ |
D( p) , при этом выделим два варианта. |
||||||
|
1. В |
случае |
дифференцирующего |
||||
+ 20, дБ/дек. |
фильтра первого порядка D( p) |
p |
1 . |
||||
При этом асимптотическая логарифми- |
|||||||
|
|||||||
lgω |
ческая амплитудно-частотная характери- |
||||||
стика, соответствующая Wк ( p) , имеет |
|||||||
|
|||||||
|
вид, показанный на рис. 11.9. |
|
|
||||
Рис. 11.9. ЛАЧХ контура |
Как видим, высокочастотная помеха |
||||||
будет проходить с постоянным коэффи- |
|||||||
с фильтром первого порядка |
|||||||
циентом усиления. |
|
|
|
|
|||
|
|
|
|
|
|||
При использовании фильтра второго порядка D( p) |
2 p2 |
2d p |
1 |
соответствующая асимптотическая логарифмическая амплитудночастотная характеристика представлена на рис. 11.10. Очевидно, что в этом случае высокочастотная помеха будет отфильтровываться контуром.

11.3. Метод локализации |
|
|
389 |
|||
2. Таким образом, порядок |
|
L, дБ |
|
|||
D(μ p) |
полинома |
должен быть |
|
|
|
|
выше порядка требуемого диф- |
|
|
– 20 дБ/дек. |
|||
ференцирования (хотя бы на |
+ 20 дБ/дек. |
|||||
|
||||||
единицу). Только при этом ус- |
|
|
lg ω |
|||
ловии в контуре быстрых дви- |
|
|
||||
|
|
|
||||
жений будет подавляться высо- |
|
|
|
|||
кочастотная помеха измерения. |
Рис. 11.10. ЛАЧХ контура с фильтром |
|||||
Выбирать конкретные чис- |
||||||
|
второго порядка |
|||||
ленные |
значения |
параметров |
|
|||
|
|
|
||||
фильтрующего полинома необходимо с учетом условий разделимости |
||||||
движений в системе. На практике можно пользоваться следующими |
||||||
соотношениями: |
|
|
|
|
||
|
|
0,1T |
* , |
|
(11.40) |
|
|
|
d (0,5 0,7), |
|
|||
|
|
|
|
|||
где T * – желаемая постоянная времени, tn* |
3T * . |
|
||||
|
11.3.6. СИСТЕМЫ ПРОИЗВОЛЬНОГО ПОРЯДКА |
Обсудим задачу синтеза системы управления на основе метода локализации для объекта, математическая модель которого имеет вид уравнения произвольного порядка
|
|
y(n) |
|
f (t, y, , y(n 1) ) |
b(t, y, , y(n 1) )u , |
(11.41) |
|||||
|
fmax , |
|
b( ) |
|
bmax |
|
|
|
|
|
|
где |
f ( ) |
|
|
и b(t, y) |
0 |
y |
y , t [0; |
) . |
|||
Желаемая динамика задается эталонным уравнением n-го порядка |
|||||||||||
|
|
|
|
|
|
y(n) |
F( y, y,..., y(n 1) ,v) . |
|
(11.42) |
||
Закон управления следует формировать в виде |
|
|
|||||||||
|
|
|
|
|
|
u |
k F ( ) |
y(n) . |
|
|
(11.43) |
