
- •СОДЕРЖАНИЕ
- •ВВЕДЕНИЕ
- •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. АДАПТИВНЫЕ СИСТЕМЫ. ПРИНЦИПЫ ПОСТРОЕНИЯ
- •СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

6. ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ НЕЛИНЕЙНЫХ СИСТЕМ
Рассмотренные до сих пор переходные процессы в нелинейных системах часто носили качественный характер, поскольку время в явном виде не фигурировало. Для количественной оценки качества переходного процесса необходимо построить зависимость x = f (t) .
6.1. Построение переходного процесса методом фазовой плоскости
Пусть фазовая траектория системы имеет вид, приведенный на рис. 6.1, а. Возьмем начальную точку M 0 (x0 , x&0 ) . Разобьем ось x на участки (жела-
тельно равные) x1 - x0 , x2 - x1, x3 - x2 и т.д. Если интервалы xi - xi-1 малы, то можно считать, что на интервале скорость изменения координаты неизменна и определяется средней величиной
x& |
= |
x&i + x&i-1 |
= |
Dxi |
= |
xi |
- xi-1 |
. |
(6.1) |
|||||
|
|
|
|
|
|
|
|
|||||||
ср.i |
|
|
2 |
|
|
Dti |
Dti |
|
||||||
|
|
|
|
|
|
|||||||||
Отсюда |
|
|
|
Dxi |
|
|
|
xi - xi-1 |
|
|
|
|||
|
Dt = |
|
= 2 |
, |
(6.2) |
|||||||||
|
|
|
||||||||||||
|
|
i |
|
x&ср.i |
|
|
x&i + x&i -1 |
|
||||||
|
|
|
|
|
|
|
||||||||
|
|
|
ti +1 = ti + Dti. |
|
|
|
(6.3) |
Естественно, отрезки Dti в общем случае будут различными по величине. Точность построения зависит от величины интервалов разбиения оси на участки. Кривая переходного процесса x(t) в данном случае представлена на
рис. 6.1, б.
x&
x&8
x6 |
x7 |
x3 x2 x1 x0 |
x |
|
x5 |
x4 |
|
|
x8 |
|
x&0 |
M 0 |
x&5 |
|
x&1 x&4 |
|
x&3 |
|
x&2 |
|
x
x0 g
x1 g
x2 g
x3 g
x4 |
0 |
|
|
g |
|
|
|
|
x8g |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
x5 |
|
1 |
2 |
3 |
4 |
5 |
6 |
7 t |
|||||||
|
g |
gx |
|||||||||||||
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
x6 g
Рис. 6.1. Построение переходного процесса по фазовой траектории: а – фазовая траектория, б – кривая x = f (t)
117

6.2. Построение переходного процесса методом гармонической линеаризации
Рассмотрим колебательное звено второго порядка. В линейном звене свободные движения (при v(t ) = 0 ) соответствуют характеристическому уравне-
нию T 2l2 + 2zT l +1 = 0 , корни которого l = g ± jw . Эти движения определя- |
|||||||||
ются выражением |
|
|
|
|
1,2 |
|
|
||
|
|
|
|
|
|
|
|||
|
|
|
|
|
x = A egt sin (wt + j), |
|
(6.4) |
||
|
|
|
0 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
1 - x2 |
|
|
|
||
причем g = |
|
= const, w = |
|
|
|
= const , так как |
показатель |
колебательности |
|
T |
|
|
|||||||
x = const . |
|
|
|
T |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
В нелинейном звене g = g(t ) и w = w(t ), поэтому |
|
||||||||
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
òg(t )dt |
æ t |
ö |
|
|
|
|
|
x = A0e0 |
sin ç òw(t)dt + j÷. |
(6.5) |
||||
|
|
|
|
|
|
|
è 0 |
ø |
|
В связи с этим коэффициенты гармонической линеаризации в общем случае будут функциями A, w, g. Выражение (4.15) с учетом этого примет вид:
ìR e [D (jw, A, g, П )] = 0, |
(6.6) |
í |
|
îIm [D (jw, A, g, П )] = 0. |
|
По этим выражениям строят диаграммы качества нелинейных колебаний. Это семейства линий g = const и w = const на плоскости (A, П), где П – параметр, например, коэффициент усиления системы.
Для этого построения из первого уравнения (6.6) выделяют w = f1(A, g, П) и подставляют во второе, получают А = f2 (П, g) . Придавая g постоянные значения, строят семейство, показанное на рис. 6.2, а. Затем, подставляя А = f2 (П, g) в w = f1 (A, g, П) , исключают g и получают А = f3 (П, w) - рис. 6.2, б.
Отметим, что для линейных САУ линии на рис. 6.2 вертикальны.
Из рис. 6.2, а следует, что при g = 0 можно найти амплитуду автоколебаний (для точки C – амплитуда AC ). Слева от заштрихованной линии– затухающие колебания, справа – расходящиеся.
Если зафиксировать параметр П, то переходному процессу соответствует движение изображающей точки по вертикали.
Если параметр П соответствует точке L (затухающие колебания), то переходный процесс будет от начальной точки M0 с амплитудой A0 до A ® 0.
Если же П соответствует точке E, то характер процесса зависит от выбора начальных условий. Если это точка M02, то колебания расходятся до автоколе-
118

баний. Если это точка M03, то колебания затухают до автоколебаний. В обоих |
||||||||||||||||||
случаях амплитуда автоколебаний будет A4. |
|
|
|
|
|
|
||||||||||||
|
|
При этом частота изменяется в соответствии с графикомA = f3 (П,w) |
по |
|||||||||||||||
вертикали. На рис. 6.3, б линии w = 0 соответствуют монотонные переходные |
||||||||||||||||||
процессы. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
a |
|
A |
|
g < 0 g = 0 |
б |
A |
w = 0 |
|
w0 |
|
|
||||
|
|
|
|
|
A |
M |
g > 0 |
|
|
|
|
|||||||
|
|
|
|
|
01 g0 M 03 |
|
|
|
|
w |
|
|
||||||
|
|
|
|
A1 |
0 |
|
г |
|
|
|
|
A |
|
|
|
1w2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
w3 |
|
||
|
|
|
|
|
|
|
|
|
|
C |
A2 |
|
|
|
|
|
||
|
|
|
AC A |
g |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
2 |
|
|
|
|
|
A3 |
|
|
|
|
|
|
|
|
|
|
|
A3 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
g |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
A |
g4 |
|
|
|
|
|
A4 |
|
|
|
w4 |
|
||
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
L |
|
E |
M |
02 |
П |
|
П = const |
П |
|
|||
|
|
|
|
|
П = const |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 6.2 |
|
|
|
|
|
|
|
|
На основании этих рассуждений можно приближенно оценить переход- |
||||||||||||||||
ной процесс, определив при П = const огибающую и частоту колебаний. |
|
|||||||||||||||||
|
|
Из |
графика |
на |
рис. 6.2, а определяют T (A ) = - |
1 |
, |
а из рис.6.2, б – |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
i |
i |
gi |
(Ai ) |
|
|
|
|
|
|
2p . |
|
|
|
|
|
|
|
|
|
|
|
||||
t |
(A ) = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
i |
i |
gi (Ai |
) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
Огибающую можно построить, используя свойство экспоненты, заклю- |
||||||||||||||||
чающееся |
в |
том, что |
|
длина |
подкасательной |
равна |
постоянной |
времени |
||||||||||
(рис. 6.3). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A2 |
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
142443 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
T0 (A0 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 42 43 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
T1 (A1 ) |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 6.3 |
|
|
|
|
|
|
T (A ) = - 1 |
; T (A ) = - 1 |
; T (A ) = - 1 ; и т. д. |
|
|
|
|
|
|
||||||||||
0 |
0 |
|
g0 |
1 |
1 |
|
g1 |
|
2 |
2 |
g2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
119 |

|
|
|
Аналогично оценивается и период колебаний: |
t |
(A ) = |
2p |
; |
t (A ) = |
2p |
; |
|||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
0 |
0 |
w0 |
|
1 |
1 |
w1 |
|
|
|
|
|
2p |
|
|
|
|
|
|
|
||
t |
|
(A ) = |
и т.д. |
|
|
|
|
|
|
|
|
||
2 |
|
|
|
|
|
|
|
|
|
||||
|
2 |
|
w2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Окончательно зависимость A = f (t) имеет вид рис. 6.4.
A
A @ A4 , w @ w3
0
2p/w0 |
4p/w1 |
6p/w2 |
t |
8p/w3 10p/w4 12p/w5 |
|||
T0 (A0 ) |
T1(A1) |
T2 (A2 ) |
T3 (A3 ) |
Рис. 6.4
6.3. Компьютерное моделирование переходного процесса
При расчетах переходных процессов рассмотренными методами возникают трудности, связанные с громоздкими построениями и невысокой точностью. Кроме того, наличие в системе нескольких нелинейностей усложняет решение задачи. Существенные ограничения накладывает порядок системы и необходимость вариации величины и характера внешних воздействий. От этих недостатков в значительной степени свободен метод компьютерного моделирования процессов.
Рассмотрим систему, структура которой приведена в примере5.1 на рис. 5.3, а ( M C учитывать не будем). Выполним исследование системы при подаче на вход гармонического воздействия. Будем считать, что усилитель K1 имеет характеристику типа «насыщение».
Схема моделирования системы в SIMULINK изображена на рис. 6.5. Нелинейный элемент представлен звеном, имеющим коэффициент
усиления в линейной зоне, равный единице. При моделировании исследовано влияние величины "a" нелинейного элемента на характер ошибки при гармоническом воздействии и переходной процесс системы при ступенчатом воздействии.
При подаче на вход системы сигнала v(t ) = 0,18sin 0, 25t осциллограммы выходного сигнала изображены на рис. 6.6.
120

|
1 |
|
|
|
v |
|
|
|
НЭ |
+ |
x |
|
a |
|
|
|
|||
+ |
|
-a |
||
|
|
|
a |
|
|
|
K1 |
= 23000 |
-a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
.052s+1 |
|
1 |
|
.0047 |
|
||||||
|
|
|
s |
|
|
|
.12s+1 |
|
|
|||
|
.003s+1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 6.5. Схема моделирования нелинейной системы в SIMULINK
Рис. 6.6, а соответствует значению "a" =10 , при котором НЭ не входит в
насыщение, при этом ошибка x |
= 4 ×10-4 . Рис. 6.6, б соответствует значению |
|
max |
|
|
"a" = 8 , при котором начинает сказываться |
насыщение. Выходной сигнал от- |
|
личен от синусоиды, ошибка x |
= 7,5 ×10-3 ; рис. 6.6, в соответствует "a" = 6 – |
|
max |
|
|
НЭ входит в глубокое насыщение. |
|
|
а |
б |
в |
Рис. 6.6
При подаче на вход системы сигнала1(t) при значении "a" = 25000 , нелинейный элемент не входит в насыщение, переходной процесс и ошибка соответствуют линейной системе рис. 6.7, а. При a =10000 , нелинейный элемент входит в насыщение; переходной процесс затягивается рис. 6.7, б. При дальнейшем уменьшении "a" НЭ входит в глубокое насыщение, система становится «вялой», t p возрастает, перерегулирование уменьшается рис. 6.7, в.
а |
б |
в |
Рис. 6.7
121