- •СОДЕРЖАНИЕ
- •ВВЕДЕНИЕ
- •1. ОБЩИЕ СВЕДЕНИЯ
- •1.1. Основные понятия
- •1.2. Классификация импульсных систем
- •1.3. Примеры дискретных систем
- •2. МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ДИСКРЕТНЫХ СИСТЕМ
- •2.1. Понятие систем с дискретным временем
- •2.2. Решетчатые функции и разностные уравнения
- •2.3. Дискретное преобразование Лапласа и его свойства
- •2.4. Решение разностных уравнений
- •2.6. Представление данных в импульсной системе
- •2.7. Частотные свойства дискретного преобразования
- •2.8. Восстановление данных
- •2.9. Импульсная передаточная функция разомкнутой системы
- •2.10. Импульсная передаточная функция замкнутой системы
- •2.11. Процессы в импульсных системах
- •3. АНАЛИЗ УСТОЙЧИВОСТИ ДИСКРЕТНЫХ СИСТЕМ
- •3.1. Условия устойчивости
- •3.2. Билинейное преобразование
- •3.3. Критерий Рауса-Гурвица
- •3.4. Критерий Найквиста
- •3.5. Логарифмический критерий Найквиста
- •3.6 Критерий Михайлова
- •4. ТОЧНОСТЬ ДИСКРЕТНЫХ СИСТЕМ
- •5.СИНТЕЗ ДИСКРЕТНЫХ СИСТЕМ
- •5.3. Синтез регулятора с отставанием по фазе
- •5.4 Синтез регулятора с опережением по фазе
- •5.5. Цифровые ПИД-регуляторы
- •5.6. Особенности реализации дискретной коррекции
- •6. МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ЛИНЕЙНЫХ ДИСКРЕТНЫХ СИСТЕМ В ПРОСТРАНСТВЕ СОСТОЯНИЯ
- •6.1. Уравнения состояния дискретных систем и схемы моделирования
- •6.2. Решение уравнений состояния
- •6.3. Основные формы уравнений состояния импульсных систем
- •6.4. Преобразование уравнений состояния
- •6.5 Управляемость и наблюдаемость дискретных систем
- •НЕЛИНЕЙНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ
- •1. ОБЩИЕ СВЕДЕНИЯ О НЕЛИНЕЙНЫХ СИСТЕМАХ
- •1.1. Основные понятия и особенности нелинейных систем
- •1.2. Характеристики типовых нелинейностей и их соединений
- •1.3. Классификация и примеры нелинейных систем
- •1.4. Методы исследования нелинейных систем
- •2. МЕТОД ФАЗОВОГО ПРОСТРАНСТВА
- •2.1. Фазовое пространство
- •2.2. Методы построения фазовых портретов
- •2.3. Анализ нелинейной системы с насыщением
- •3. АНАЛИЗ НЕЛИНЕЙНЫХ СИСТЕМ НА ОСНОВЕ МЕТОДОВ А.М. ЛЯПУНОВА И В.М. ПОПОВА
- •3.3. Анализ устойчивости нелинейных систем методом В. М. Попова
- •4. МЕТОД ГАРМОНИЧЕСКОЙ ЛИНЕАРИЗАЦИИ
- •4.1. Сущность метода гармонической линеаризации
- •4.2. Определение амплитуды и частоты автоколебаний
- •4.2.1. Аналитический метод
- •4.2.2. Графоаналитический метод
- •5. ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ В НЕЛИНЕЙНЫХ СИСТЕМАХ
- •5.1. Определение параметров и устойчивости вынужденных колебаний (задача Дуффинга)
- •5.2. Метод эллипса
- •5.3. Метод Гольдфарба
- •6. ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ НЕЛИНЕЙНЫХ СИСТЕМ
- •6.1. Построение переходного процесса методом фазовой плоскости
- •6.3. Компьютерное моделирование переходного процесса
- •7. КОРРЕКЦИЯ НЕЛИНЕЙНЫХ СИСТЕМ
- •7.1. Коррекция нелинейной системы с помощью обратной связи
- •7.2. Синтез компенсационных моделей
- •7.3. Метод вибрационной линеаризации
- •ОПТИМАЛЬНЫЕ И АДАПТИВНЫЕ СИСТЕМЫ
- •1. ОБЩИЕ СВЕДЕНИЯ
- •2. ОПТИМИЗАЦИЯ ПАРАМЕТРОВ
- •3. ОЦЕНКА ТОЧНОСТИ НАСТРОЙКИ ОПТИМАЛЬНЫХ ПАРАМЕТРОВ
- •4. ВЛИЯНИЕ РАЗБРОСА ПАРАМЕТРОВ НА ДИНАМИКУ СИСТЕМЫ
- •5. УЧЕТ ФИЗИЧЕСКИХ ОГРАНИЧЕНИЙ
- •6. СИНТЕЗ ОПТИМАЛЬНЫХ РЕГУЛЯТОРОВ ПО КВАДРАТИЧНОМУ КРИТЕРИЮ КАЧЕСТВА
- •7. МЕТОДИКА СИНТЕЗА ОПТИМАЛЬНОЙ ПЕРЕДАТОЧНОЙ ФУНКЦИИ РЕГУЛЯТОРА
- •8. ВЫБОР ВЕСОВЫХ КОЭФФИЦИЕНТОВ
- •10. АДАПТИВНЫЕ СИСТЕМЫ. ПРИНЦИПЫ ПОСТРОЕНИЯ
- •СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
2.3. Анализ нелинейной системы с насыщением
Пусть структурная схема НСАУ имеет вид, приведенный на рис. 2.17.
v = 0 |
e |
U (e) |
y |
K1 |
|||
+ |
å |
|
|
|
- |
|
|
Рис. 2.17. Структурная схема НСАУ с насыщением
В режиме свободного движения v(t) º 0 , и, следовательно, e(t ) = -y (t ).
Передаточная функция линейной части системы |
Y (s ) |
|
= |
K1 |
. Диф- |
|||||||||||||||||||||
U (s ) |
s (T s +1) |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
||
ференциальное нелинейное уравнение системы T p2 y (t ) + py (t ) = K U (e). |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|||
Заменяем |
|
p на |
d |
|
и y (t ) на - е (t), получим T |
d 2e(t ) |
+ |
de(t ) |
= -K U (e ), |
|||||||||||||||||
|
dt |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
dt2 |
|
|
dt |
1 |
|||||
ìUm , |
|
e ³ a, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ï |
Um |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где U (e )= íï |
е, |
-a < e < a , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
ï а |
|
, |
e £ -a. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ï-U |
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Нелинейная система в данном случае будет описываться системой урав- |
||||||||||||||||||||||||||
нений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ìT |
d 2e |
+ |
de |
= -K U |
m |
, |
|
е ³а, |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
ï |
1 |
dt |
2 |
|
dt |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
ï |
|
|
|
|
U |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
ï |
|
d 2e |
|
de |
|
m |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
T |
|
|
+ |
|
= -K |
|
|
e, |
- а < е < a , |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
í |
1 dt 2 |
|
dt |
1 |
a |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
ïT |
d 2e |
+ |
de |
= K U |
m |
, |
|
|
|
е £ -а. |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
ï |
1 |
dt |
2 |
|
dt |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Разделив правую и левую части уравнений наa, |
и обозначив K1 |
U m |
= K , |
||||||||||||||||||
|
|||||||||||||||||||||
запишем систему уравнений в виде: |
|
|
|
|
|
|
|
|
|
a |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
ìT |
d 2 |
æ e ö |
+ |
d æ e ö |
= -K , |
|
|
е |
³1, |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
2 |
ç ÷ |
ç ÷ |
|
|
|
|
|
|||||||||||||
ï 1 |
dt |
|
|
|
|
|
а |
|
|
|
|
||||||||||
ï |
|
è a ø |
|
dt è a ø |
|
|
|
|
|
|
|
|
|||||||||
ï |
d 2 |
æ e ö |
|
d æ e ö |
æ e ö |
|
|
е |
|
|
|
||||||||||
íT1 |
|
|
ç |
|
÷ |
+ |
|
ç |
|
÷ |
= -K ç |
|
÷ |
, -1 < |
|
<1, |
|
|
|||
dt |
2 |
|
|
|
|
а |
|
|
|||||||||||||
ï |
|
è a ø |
|
dt è a ø |
è a ø |
|
|
|
|
|
|||||||||||
ï |
d 2 |
æ e ö |
|
d æ e ö |
|
|
е |
|
|
|
|
|
|
||||||||
ïT1 |
|
|
ç |
|
÷ |
+ |
|
ç |
|
÷ |
= K , |
|
|
|
£ -1. |
|
|
|
|
||
dt |
2 |
|
|
|
а |
|
|
|
|
||||||||||||
î |
|
è a ø |
|
dt è a ø |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
89 |
Разделив правую и левую части на K и обозначив e = x , получим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
ì |
|
|
d |
2 |
x |
|
|
|
dx |
|
|
|
|
||
ïT K |
|
|
|
+ |
= -1, |
х ³1, |
|
||||||||
d (Kt 2) |
|
|
|||||||||||||
ï |
1 |
|
|
d (Kt ) |
|
|
|
||||||||
ï |
|
|
d 2 x |
|
|
|
|
dx |
|
|
|
|
|||
ï |
|
|
|
|
|
|
|
|
|
|
|||||
íT K |
|
|
|
|
|
|
+ |
|
|
= -x, |
-1 < х <1, |
|
|||
d (Kt 2) |
|
|
|
||||||||||||
ï |
1 |
|
|
d (Kt ) |
|
|
|
||||||||
ï |
|
|
d |
2 |
x |
|
|
|
dx |
|
|
|
|
||
ïT K |
|
|
|
+ |
=1, |
х £ -1. |
|
||||||||
|
|
|
|
|
|
|
|
||||||||
ï |
1 |
|
|
|
|
|
2 |
|
|
d (Kt ) |
|
|
|
||
î |
|
d (Kt ) |
|
|
|
|
|
|
|
|
|||||
Введем новую переменную t = Kt , тогда |
|
|
|||||||||||||
|
|
|
d |
2 x |
|
|
dx |
ì-1, |
х ³1, |
|
|||||
|
|
|
|
|
ï |
|
|
|
(2.11) |
||||||
|
T K |
|
|
|
|
+ |
|
|
|
= í-х, |
-1 |
< х < +1, |
|||
|
|
|
|
|
|
|
|
||||||||
|
1 |
|
dt2 |
|
|
dt |
ï |
|
х £ -1. |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
î1, |
|
|
Таким образом, математическая модель нелинейной системы с насыщением описывается системой дифференциальных уравнений в безразмерных величинах.
Исследуем динамику системы (2.11) методом изоклин. Заменим |
d |
= s , |
|||||||||||
d t |
|||||||||||||
тогда структурная схема принимает вид, изображенный на рис. 2.11. |
|
||||||||||||
|
|
||||||||||||
v = 0 |
|
e |
|
|
|
|
|
|
|
|
|
|
|
å |
|
|
|
|
U (e) |
1 |
|
|
|
||||
|
|
|
|
|
|
s (KT1s +1) |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.18. Структурная схема НСАУ, описанная системой (2.11)
Составим уравнения изоклин для двух случаев.
1. Для линейной области ( -1 < x < +1) уравнение (2.11) приводится к виду
|
|
T K |
d 2 x |
+ |
|
dx |
|
+ x = 0 , |
-1 < x <1. |
||||||
|
|
d t2 |
|
dt |
|||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
||||
Обозначаем |
dx |
= x& , тогда |
d 2 x |
= x& |
dx& |
, |
T Kx& |
dx& |
+ x& |
+ x = 0 . Для построения |
|||||
dt |
dt2 |
dx |
dx |
||||||||||||
|
|
|
|
|
|
1 |
|
|
изоклин обозначим |
dx& |
= Ni и подставляем в предыдущее уравнение. Оконча- |
|
dx |
|||
|
|
||
тельно получим T1KNi x& + x& + x = 0 , откуда |
x& = - |
|
x |
|
. |
(2.12) |
1 |
+ KT N |
|
|||
|
|
1 |
i |
|
90
2. Для нелинейных областей (вывод |
уравнений аналогичен |
предыдуще- |
|||||
му) |
|
|
|
|
|
|
|
x& = - |
1 |
|
|
, |
x ³1; |
|
|
|
|
|
|
||||
|
|
1 + KT1 Ni |
, |
(2.13) |
|||
|
|
1 |
|
|
|
||
x& = |
|
|
, x £ -1. |
|
|||
1 + KT N |
|
|
|||||
|
|
1 |
i |
|
|
Для упрощения расчетов принимаем KT1 =1. Результат построения фазового портрета системы представлен на рис. 2.19.
N = 0 |
x& |
N = -2 |
|
N |
= -1 |
||
|
|||
N =1 |
|
|
|
|
|
N = -¥ |
|
|
|
g |
|
|
|
N =1 |
|
N = -2 |
|
N = 0 |
|
|
|
Рис. 2.19. Фазовый портрет НСАУ с насыщением
Анализ фазовых траекторий показывает, что насыщение элементов системы приводит к увеличению длительности переходного процесса, так как скорость движения на различных участках траектории ограничена, с другой стороны в системе уменьшается перерегулирование.
2.4.Анализ нелинейной системы с насыщением
иместной ОС по скорости
Пусть структурная схема системы имеет вид, приведенный на рис. 2.20.
v = 0 |
e |
e1 |
Um |
U (e1 ) |
|
å |
å |
|
-Um |
|
|
|
|
|
|
|
|
|
KОС |
K1 |
ω |
1 |
y |
T1 p +1 |
|
p |
|
Рис. 2.20. Структурная схема НСАУ с насыщением и местной обратной связью по скорости
91
|
Согласно схеме |
E1 ( p) = E ( p) - KОС pY ( p), где E ( p) = V ( p) - Y ( p) . |
||||||
|
В |
режиме |
свободного |
движенияV ( p) º 0 , E ( p) = -Y ( p) , и |
||||
E |
( p) = 1 + K |
ОС |
p |
] |
E ( р) и исходную систему можно представить в соответствии |
|||
1 |
[ |
|
|
|
|
|
с рис. 2.21.
v |
e |
e1 |
Um |
|
1 |
+ KОС p |
-Um |
|
|
|
U (e1) |
K1 |
|
y |
|
p(T1 p +1) |
|
|
|
|||
|
|
|
|
Рис. 2.21. Преобразованная структурная схема НСАУ
Перейдем по рассмотренному в пункте2.3 алгоритму к структурной схеме, описываемой безразмерными переменными (рис. 2.22).
v = 0 |
e |
e1 |
1 |
|
|
1+ KОСKs |
-1 |
|
|
|
U (e1 ) |
1 |
|
|
x |
|
s(KT1s +1) |
|
|
|
|
|
|
|
|
Рис. 2.22. Преобразованная структурная схема НСАУ
Запишем дифференциальные уравнения НСАУ с обратной связью по пе-
редаточной функции |
|
|
X (s) |
= |
|
|
|
1 |
|
|
|
|
. Учитывая, что X (s) = -E (s) и |
|||||||||||||
|
|
|
|
(s )ù |
|
s (KT s +1) |
||||||||||||||||||||
|
|
|
U éE |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
d |
|
|
|
|
ë 1 |
û |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
s = |
, получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
d t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
ì1, |
|
|
|
|
|
|
|
|
|
|
|
|
e1 (t) ³1, |
|
||
|
|
|
d 2e(t ) |
|
de(t ) |
ï |
|
|
|
|
|
|
|
de(t ), |
|
|
||||||||||
|
|
-KT |
- |
= íïe(t )+ KK |
oc |
|
-1 < e (t )<1, |
(2.14) |
||||||||||||||||||
|
|
|
|
d t |
|
|
||||||||||||||||||||
|
|
1 dt2 |
|
ï |
|
|
|
|
|
|
dt |
|
|
1 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
ï-1, |
|
|
|
|
|
|
|
|
|
|
|
|
e1 (t )£1. |
|
||
|
|
|
|
|
|
|
|
|
|
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Уравнения для построения фазовых траекторий методом изоклин в -ли |
|||||||||||||||||||||||||
нейной области: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
x& = - |
|
|
|
e(t) |
|
|
|
|
, |
|
|
(2.15) |
||||||
|
|
|
|
|
|
|
|
1 + KT N |
i |
+ KK |
oc |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
в областях насыщения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
x& = - |
|
e(t) |
|
|
. |
|
|
|
|
(2.16) |
||||||
|
|
|
|
|
|
|
|
|
|
|
+ KT N |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
i |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||
92 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из полученных уравнений видно, что стыковку изоклин необходимо осуществлять на линиях переключения, описываемых уравнениями:
|
|
|
|
|
e(t) + KKOC x& =1, |
|
(2.17) |
|||
|
|
|
|
|
e(t )+ KKOC x& = -1. |
|||||
|
|
|
|
|
|
|||||
Пусть |
KKOC =1, KT1 =1. |
Некоторые |
результаты |
расчетов: N =1, |
||||||
x& = - |
e(t) |
; |
N = -1, x& = -e(t); |
N = -2,5 , |
x& = 2e(t); |
N = -3 , x& = e(t); |
||||
|
||||||||||
3 |
|
|
e(t) |
|
|
e(t) |
|
|
|
|
N = -0,5 , |
x& = - |
; N = 0 , x& = - |
. |
|
|
|||||
|
|
|
|
|||||||
|
|
|
1,5 |
|
2 |
|
|
|
Фазовый портрет приведен на рис. 2.23.
Сопоставляя результаты исследования динамики системы с ОС по скорости с результатами исследования без ОС можно сделать вывод: коэффициент усиления в цепи ОС по скорости KOC можно подобрать так, что в переходном процессе будет отсутствовать перерегулирование, а время переходного процесса будет минимальным.
Подобные результаты были положены в основу построения позиционных систем управления, которые широко распространены в современной технике.
x& = de dt
N =1
0
e
N =1
N = -1
Рис. 2.23. Фазовый портрет НСАУ с насыщением и местной обратной связью по скорости
|
|
|
1 |
|
|
|
|
1 |
|
|
||||||
1 |
|
|
|
KT s2 |
|
1 |
|
|
KT |
|
|
|||||
Учитывая, что |
|
= |
1 |
|
|
= |
|
1 |
|
, построим структурную |
||||||
s (KT1s +1) |
|
1 + |
1 |
|
|
s2 |
|
1 + |
1 |
|||||||
|
|
KT1s |
KT1s |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
схему позиционной системы (рис. 2.24):
93
U Д = K Д jД |
|
|
У |
U |
* |
x& |
|
У2 |
|
|
|
xв |
|
|
|
1 |
|
|
|
|
|||
|
|
|
1 |
|
|
å |
|
1 |
|
å |
|
K Д |
+ |
å |
- 1 |
1 |
+ |
- 1 |
|
1 |
|||
|
|
- 1 |
|
- 1 |
+ |
- |
|
- |
- |
|
UC = KC jC |
U x& |
Ускорение |
||
|
|
|
|
|
|
|
Cкорость |
KOC |
|
|
K у |
|
|
|
|
Положение |
|
|
|
|
|
|
|
|
|
|
|
|
|
У3 |
x = &x& |
|
x2 = x& |
|
1 |
3 |
|
1 |
x |
|
1 |
|||
KT1 |
|
s |
s |
|
1 |
|
|
|
|
s |
|
|
|
|
Рис. 2.24. Структурная схема позиционной системы
Пусть в некоторый момент времени на вход подано воздействиеxв достаточно большой величины. Учитывая, что система находилась в состоянии покоя в начальный момент времени, сигналы в цепях ОС по положению, по скорости и ускорению равны нулю. В этом случае усилители У1, У2 оказываются в насыщении, У3 выполняет роль регулятора ускорения. По мере движения системы сигнал ОС по скорости U x& возрастает и когда достигает определенного уровня, У2 выходит из насыщения. Затем на очень коротком отрезке времени сигнал на выходе У2 оказывается равным нулю, следовательно ускорение равно нулю, а скорость x2 достигла своего максимального значения. При дальнейшем движении возрастает сигнал U Д и достигает значения, при кото-
ром У1 выходит из насыщения. В этот момент U* x& <U x& и усилитель У2 снова оказывается в насыщении, но с обратным знаком. Система переходит в режим торможения. Скорость движения снижается и когда ошибка рассогласования на входе У1 станет равной нулю, то и скорость равна нулю, система пришла в точку позиционирования (рис. 2.25).
x
0 |
|
x& |
t |
0 |
|
&x& |
t |
0 |
|
|
t |
Рис. 2.25. Графики изменения во времени выходной координаты x , её скорости x& и ускорения &&x
94
2.5.Исследование динамики релейной системы
сидеальным двигателем
Для улучшения динамических свойств систем в контур управления часто включают элементы, обладающие релейными характеристиками. Рассмотрим структуру простейшей системы (рис. 2.26).
v(t ) e(t ) x&2 =U x&1 = x2 x = x1
Рис. 2.26. Структура релейной системы
На схеме x = |
y |
|
, s =1× p, t = t , |
U = sign e - sign x1 , |
v º 0 . Запишем |
||||||||
Umax K |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
уравнения согласно структурной схеме: |
|
|
|
|
|||||||||
|
|
|
|
|
dx1 |
= x2 , |
(2.18) |
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
dt |
|
|
|
|
||||
|
|
|
|
dx2 |
=U , |
(2.19) |
|||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
dt |
|
|
|
|
||||
|
|
|
U = -sign x1 . |
(2.20) |
|||||||||
Интегрируя (2.18) и (2.19), находим: |
|
|
|
|
|||||||||
|
|
|
x2 (t ) = x2 (0) +Ut , |
(2.21) |
|||||||||
x1 |
(t )= x1 |
(0 )+ x2 (0 )t + |
U t 2 |
|
x1 (0 )= x2 (0 )= 0. |
(2.22) |
|||||||
|
|||||||||||||
|
|||||||||||||
|
|
|
2 |
|
|
|
Полагая x1 (0)= x2 (0)= 0 , построим в соответствии с (2.20 – 2.22) графики
U (t ), x1 (t ) и x2 (t ) (рис. 2.27). |
|
|
|
|
|
|
|
|
|
|
|
|
Отметим, что площадь, ограниченная кривой x2 = |
dx |
и осью абсцисс, ве- |
||||||||||
dt |
||||||||||||
|
tпп dx |
|
|
|
|
|
|
|
||||
|
|
(t |
|
|
)=1. |
|
|
|||||
|
ò dt |
dt = x |
|
|
|
|
||||||
личина постоянная и равна единице, так как |
|
1 |
п.п. |
|
|
|||||||
|
|
|
1 |
|
|
|
|
|
|
|||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
Иная кривая, не может проходить выше x2 |
в силу условия |
|
x2 |
|
=U £1, при |
|||||||
|
|
|||||||||||
этом x2max – максимальная допустимая скорость. |
|
|
|
|
|
|
|
|
|
|
95
x 1 (t )
0 x2 (t )
x2 max
x2ср
0
U (t )
0
·
tп.п. = Tmin
tп
tп
tп
Рис. 2.27
Поясним физический смысл рис. 2.27. В начальный момент времени при U = 1 к двигателю прикладывается напряжение, при котором он развивает полный момент и начинает разгоняться. В момент tп прикладывается обратное напряжение и двигатель начинает тормозить. Если правильно рассчитать tп, то процесс будет апериодическим, без ошибок, с минимальными затратами по времени. Приведенные характеристики являются идеализированными, так как трудно анализировать момент окончания переходного процесса и обеспечить равенство U = 0 и x2 = 0.
|
Докажем, что tп.п. |
является минимальным. Действительно, из уравнения |
||||||||||||||||||||
x |
= x t |
, где x |
= 1, |
x |
|
= |
1 |
|
x |
|
, находим t |
= |
|
2 |
. |
|
|
|||||
|
|
|
x2 max |
|
||||||||||||||||||
1k |
2cp п.п. |
1k |
|
2cp |
|
2 2max |
|
|
п.п. |
|
|
|
|
|
||||||||
|
Допустим теперь, что x2(tп) < x2max. Тогда t |
|
= |
|
2 |
|
> |
|
2 |
. |
|
|||||||||||
|
п.п. |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x2 (tп |
) |
|
x2 max |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Следовательно, t |
|
= |
2 |
|
|
|
действительно |
|
|
является |
минимальным |
||||||||||
|
п.п. |
x2 max |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( tп.п. = Tmin ), а такая система – оптимальной по быстродействию.
Дальнейшее исследование динамики релейной системы продолжим при помощи метода изоклин на фазовой плоскости. Разделим (2.18) на (2.19):
96
|
|
dx 1 |
|
|
x 2 (t ) |
|
U = -sign x1 , следовательно dx1 |
|
1 |
|
(t )dx2 |
|
|||||
|
|
dt |
|
|
, где |
|
и после инте- |
||||||||||
|
|
|
= |
= U (t )x2 |
|||||||||||||
|
dx 2 |
|
U (t ) |
|
|||||||||||||
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
грирования |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
x |
(t )= |
1 |
x2 |
(t )+ C , |
|
|
|
|
(2.23) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
1 |
|
2U 2 |
|
|
|
|
|
|
где С – постоянная интегрирования.
По этому уравнению построим семейство парабол для разных значений
Ci и U(t) = ±1 (рис. 2.28).
При v(t) = 1(t), С = 1. Через точку x1k = C = 1 проводим параболу при U = –1. Из рис. 2.28 видно, что после перемещения системы из начального состояния x1(0) = x2(0) = 0 в конечное x1k = 1, x2k = 0, движение не заканчивается, так как в системе возникают автоколебания.
U = -1 |
x2 max |
x2 = x&1 |
|
|
|
C = 0 C > 0 |
|||
C = 0 C > 0 |
C < 0 |
|||
|
||||
C < 0 |
|
A |
|
1
x2k |
x1k |
v(t ) = é1(t )ù x1 |
|
|
|
ë |
û |
U =1
Рис. 2.28. Фазовый портрет релейной системы
Из теории линейных систем известно, что динамические свойства можно улучшить, если в сигнал управления ввести дополнительный сигнал, пропорциональный первой производной (т.е. ввести обратную связь по скорости). Используем этот метод для того, чтобы исключить автоколебания. Структурная схема системы примет вид (рис. 2.29).
v = 0 |
|
|
e |
|
|
e1 |
|
|
|
|
|
|
U = x&2 |
1 |
|
x2 = x&1 |
1 |
|
x = x1 |
||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
KОС
Рис. 2.29. Структура релейной системы с обратной связью по скорости
97
Этой структуре соответствуют следующие соотношения:
e1 = e - KOC x2 = -(x1 + KOC x2 ), |
(2.24) |
|
U = sign (e1 )= -sign (x1 + KOC x2 ). |
||
|
||
а уравнение линии переключения имеет вид: x1 + KOC x2 = 0 , откуда |
|
|
x1 = -KOC x2 . |
(2.25) |
Фазовый портрет системы представлен на рис. 2.30.
Изменяя значение KOC , можно подобрать его таким, чтобы парабола при U = –1 проходила через начало координат (точка А на линии переключении при начальном условии (x0 , 0) ), тогда переходный процесс в системе будет монотонный, апериодический, а время переходного процесса– минимальным, дру-
гими |
словами |
режим |
работы |
системы |
будет |
оптимальным |
|
быстродействию (рис. 2.30). |
|
|
|
|
|||
|
|
|
|
x2 |
|
|
|
|
|
U = +1 |
|
U = -1 |
|
|
|
|
|
|
|
А |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
x |
нач |
x0 |
x |
|
|
|
|
|
|
1 |
|
|
x1 = -KOC x2
Рис. 2.30. Фазовый портрет линейной системы с обратной связью
2 |
|
|
|
|
|
|
|
|
x2 |
1 |
||
C |
|
|
|
|
|
|
|
|
|
|||
|
g |
|
|
|
|
|
|
|
|
|
|
|
|
A g |
|
|
|
|
|
|
|
|
|
|
|
|
|
g |
|
D |
x3 |
|||||||
|
E |
|
g |
|
0 |
|
|
|||||
|
|
|
||||||||||
|
|
|||||||||||
|
|
|||||||||||
x1 |
x2 |
|
|
|
|
|
|
|
|
|
x1 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
g |
|
|||||
|
|
|
|
|
|
|
|
|
|
В |
||
|
|
|
|
|
|
|
|
H |
g |
G
Рис. 2.31. Скользящий режим в релейной системе с обратной связью
98
В том случае, когда в точке переключения реле угол наклона линии переключения становится равным наклону или меньше угла наклона касательной к фазовой траектории, по которой движется изображающая точка после переключения реле, возникают условия существования скользящего режима
(рис. 2.31).
Скользящим режимом называется режим работы релейной системы, характеризующийся колебательным движением изображающей точки вдоль линии переключения с высокой частотой и малой амплитудой.
Пусть начальное состояние системы задано точкой (x1, 0) (рис. 2. 31), от
которой изображающая точка перемещается по фазовой траектории типа1 до встречи с линией переключения в точкеС. В этой точке происходит переключение реле и изображающая точка перемещается по фазовой траектории типа2 до точки D. В точке D происходит переключение реле в другую сторону, после чего изображающая точка будет перемещаться по фазовой траектории типа1. Как только снова увеличится результирующий сигнал обратной связи, произойдет переключение реле и изображающая точка будет перемещаться по фазовой траектории типа 2 и т. д. Таким образом, изображающая точка, достигнув точки D, непрерывно перемещается к началу координат, как бы «скользя» вдоль линии переключения.
Как видно из рис. 2.31, скользящий режим возникает на участке АВ. При начальном положении изображающей точки(x2, 0) после ее прихода по траектории типа 1 в точку на линии переключенияЕ сразу начинается скользящий режим. При начальном положении изображающей точки(x3, 0) скользящий режим имеет место после переключения реле в точкеН, когда изображающая точка скользит по линии переключения в четвертом квадранте.
Быстрое движение по линии скольжения обуславливает следующие особенности, характерные для нелинейных систем: возможность получения с использованием релейного регулятора гораздо меньшего времени переходного процесса, чем, например, при использовании стандартного ПИД–регулятора; возможность получения практически конечного времени переходного процесса (времени достижения заданного состояния).
Из уравнения (2.18) с учетом (2.25) следует, что |
dx1 |
= - |
1 |
x |
. Это значит, |
|
|
||||
|
dt |
KОС |
1 |
|
|
|
|
|
что движение по линии переключения соответствует системе первого порядка, процессы не зависят от параметров прямой цепи, а определяются только величиной KОС .
Отметим, что на практике всегда реализуется режим, близкий к скользящему, отличающийся от истинно скользящего конечной частотой переключения. В самом деле, релейный элемент не может переключаться с бесконечной частотой вне зависимости от его реализации: аппаратной (реле), электронной (электронная ключевая схема) или программной.
99