
- •Глава 1. Оптимизация статических объектов 6
- •Глава 2. Основы линейного программирования 22
- •Глава 3. Способы описания динамических систем 34
- •Глава 4. Применение вариационных методов для поиска
- •Глава 1. Оптимизация статических объектов
- •1.1. Понятие статических и динамических объектов
- •1.2. Задача нелинейного программирования
- •Задачи на условный экстремум, неопределенные множители Лагранжа
- •1.4. Пример применения метода неопределенных множителей Лагранжа для поиска наибольших значений функций
- •1.5. Общая постановка задачи линейного оценивания параметров
- •Глава 2. Основы линейного программирования
- •2.1. Постановка задачи
- •2.2. Основные геометрические фигуры в линейном программировании
- •Экстремальные точки
- •2.4. Основные теоремы об экстремальных точках
- •2.5. Симплексный метод решения задач линейного программирования
- •2.6. Учет ограничений типа неравенств
- •2.7. Поиск начальной экстремальной точки
- •Глава 3. Способы описания динамических систем
- •3.1. Передаточные функции
- •3.2. Описание в форме Коши
- •3.3. Управляемость, наблюдаемость, стабилизируемость, обнаруживаемость
- •3.4. Понятие фильтра и общая задача регулирования
- •Глава 4. Применение вариационных методов для поиска оптимальных управлений
- •4.1. Понятие линейного пространства
- •4.2. Функционал и его вариация
- •4.3. Вычисление вариации функционала
- •4.4. Задача Эйлера
- •4.5. Применение уравнения Эйлера для поиска оптимального закона управления
- •4.6. Уравнение Эйлера – Пуассона и его применение
- •4.8. Неопределенные множители Лагранжа в вариационном исчислении
- •Глава 5. Вариационные задачи с подвижными границами
- •5.1.Основные виды задач с подвижными границами
- •5.2. Скольжение граничных точек по заданным траекториям
- •Глава 6. Принцип максимума
- •Постановка задачи поиска оптимального управления
- •Пояснения к получению принципа максимума
- •6.3. Динамическое программирование
- •Примеры применения динамического программирования
- •Глава 7. Аналитическое конструирование регуляторов (акор)
- •7.1. Постановка задачи
- •7.2. Решение задачи акор
- •7.3. Уравнение Риккати
- •7.4. Общие свойства решения уравнения Риккати
- •7.5. Способы решения уравнения Риккати
- •7.7. Метод диагонализации
- •Глава 8. Случайные процессы в системах управления
- •8.1. Описание случайных процессов
- •8.2. Стохастические дифференциальные уравнения
- •8.3. Прогнозирование (оценивание) значений случайных величин с использованием закона распределения
- •8.4. Линейное оценивание значений случайных величин
- •Глава 9. Фильтр калмана
- •9.1. Постановка задачи
- •9.2. Основные принципы получения формул дискретного фильтра Калмана
- •9.3. Получение формул фильтра Калмана
Глава 3. Способы описания динамических систем
3.1. Передаточные функции
Динамические системы – это такие системы, в которых имеется запаздывание передачи информации. Они обычно описываются с применением аппарата дифференциальных уравнений. Однако практика проектирования систем управления (например, выбор корректирующих устройств) показала эффективность применения других способов описания. Наиболее распространенными способа-ми описания динамических систем являются:
дифференциальные уравнения высоких порядков; частотные характеристики (амплитудные, фазовые, амплитудно-фазовые);
переходные характеристики (реакция на ступенчатое воз-действие);
импульсные переходные характеристики (реакция системы на δ – функцию);
передаточные функции;
структурные схемы линейных систем, включающие интег-рирующие, суммирующие, усилительные звенья;
описание систем в форме Коши.
34
Наиболее полными являются два последних способа описания. Поясним это на примере недостатков передаточной функции ли-нейной системы, когда появляется возможность сокращения оди-наковых множителей числителя и знаменателя.
При появлении в передаточной функции динамической системы одинаковых полиномов в числителе и знаменателе возникает во-прос о возможности сокращения таких полиномов. Простейшие примеры показывают, что при необоснованном сокращении поли-номов могут быть потеряны важные свойства рассматриваемой ди-намической системы. Пусть передаточная функция имеет вид
w( s) = ux ((ss )) = TsTs −−11 , где s – аргумент передаточной функции, x –
выходная величина, u – управляющее воздействие, T – постоян-ная, характеризующая динамические свойства объекта. Сокраще-ние числителя и знаменателя приведет к результату x( s ) = y ( s) , и,
казалось бы, можно сделать вывод, что выходной сигнал совпадает с управляющим воздействием. Однако, потеряно главное свойство, заключающееся в том, что объект является динамическим. В дан-ном примере, кроме того, объект является неустойчивым.
Известно, что передаточная функция wЗ ( s) |
замкнутой системы |
|
||||||
находится согласно правилу |
w |
(s) = |
|
|
w( s) |
, |
где w( s) – переда- |
|
|
|
|
|
|||||
|
З |
|
1 |
+ w( s) |
|
|
||
|
|
|
|
|
точная функция разомкнутой системы. Если w( s) является отноше-
нием двух полиномов w( s) = |
P(s) |
, то |
w (s) = |
|
P(s ) / Q (s) |
и воз- |
|
|
|
|
|||||
|
Q (s) |
З |
1+ P(s ) / Q (s) |
|
|||
|
|
|
никает вопрос о возможности сокращения числителя и знаменателя на Q ( s) .
Пусть объект описывается дифференциальным уравнением
a |
|
d n x |
+ |
... + a |
dx |
+ a x = b |
d mu |
+ ... + b u |
|
||
n dt n |
|
dtm |
|
||||||||
|
|
1 dt |
0 |
m |
0 |
|
|||||
функция имеет вид |
w( s) = |
|
a s n + ... + a |
|
|||||||
|
n |
0 |
|
||||||||
|
b s m + ... +b |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
m |
0 |
|
, и его передаточная
= QP(( ss)) . Введение единич-
ной обратной связи (замыкание объекта) означает, что справедливы
35
соотношения |
x( t ) = a |
d n x |
+ ... + a |
dx |
+ a x = b |
d mε |
+ ... + b ε, |
|
|
|
|
|
|
||||||
|
n dt n |
1 dt |
0 |
m dtm |
0 |
|
ε(t ) = u ( t ) − x (t ) . Исключая ε(t ) из последних соотношений, полу-
чим дифференциальное уравнение, описывающее замкнутый объ-ект:
|
a |
|
d n x |
|
+ ... + a |
|
d m +1 x |
|
+ ( a |
+ b ) |
d m x |
+ ... + ( a |
+ b )x = |
|
|||||||||||||||||||||||||||||||||||||||||
|
n dtn |
|
|
|
dtm +1 |
dtm |
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
m +1 |
m |
m |
|
|
|
0 |
|
|
0 |
|
|
||||||||||||||||||||||||||||||||||||||
|
= b |
d mu |
+ ... + b u, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
m dtm |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
которому |
|
|
|
|
соответствует |
|
|
передаточная |
функция |
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
b s m |
+ ... + b |
|
|
|
|
|
. |
Эта |
передаточная |
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
m |
|
|
|
0 |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
a s n |
+ ...a |
m +1 |
s m +1 + ( a |
m |
+ b ) s m + ... + ( a |
+b |
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
n |
|
|
|
|
|
|
m |
|
|
0 |
|
0) |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
функция также получается по формуле |
w (s) = |
|
P(s ) / Q (s) |
, где |
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
З |
1+ P(s ) / Q (s) |
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
необходимо произвести сокращение числителя и знаменателя на
Q ( s) . |
|
|
В |
частности, |
|
если |
w( s) = |
K |
|
, |
то |
|
|||||||
|
|
|
Ts −1 |
|
|||||||||||||||
|
|
|
K / (Ts −1) |
|
|
|
K |
|
|
|
|
|
|
|
|||||
w( s) = |
|
|
|
= |
|
; |
произведено сокращение |
на |
|
||||||||||
1 |
|
|
Ts + K −1 |
|
|||||||||||||||
|
+ K / (Ts − 1) |
|
|
|
|
|
|
|
|
(Ts −1) .
Все вопросы, связанные с возможностью сокращения нулей и полюсов в передаточных функциях следует решать с помощью анализа дифференциальных уравнений.
Другой пример. Пусть имеется инерционный объект, который описывается передаточной функцией w( s) = T0 sK+1 , где T0 –
большая постоянная времени (например, характеризует инерцион-ность электрической печи). Чтобы ускорить процесс разогрева, применим последовательное корректирующее звено с целью до-биться новой передаточной функции объекта передаточной функ-
цией w( s) = Tk sK+1 , где Tk значительно меньше T0 . Такое корректи-
36
рующее звено с передаточной функцией w( s) = (T0 s +1) можно реа-
Tk s +1
лизовать с помощью усилителя с большим коэффициентом усиле-ния, охваченного единичной отрицательной обратной связью с пе-
редаточной функцией w |
( s) = |
K y |
|
. После сокращения одинако- |
|
|
|
||||
y |
|
T0 s +1 |
|
|
|
|
|
|
|
вых множителей получим передаточную функцию системы в виде
w (s) = |
K |
|
, где T = |
T0 |
|
. Однако такая коррекция опасна, |
|
|
|
|
|
|
|||||
c |
Tk s +1 |
k |
K y +1 |
|
|
|||
|
|
|
|
поскольку теперь система имеет второй порядок, и на выходе кор-ректирующего устройства будет наблюдаться большой выброс, что останется незаметным при сокращении одинаковых множителей.