
- •Глава 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.2. Описание в форме Коши
Описание системы в форме Коши учитывает полную структуру динамической системы, позволяет вести анализ всех меняющихся в ней сигналов и имеет вид
dxdt1 = f1 ( x1 , x2 ,..., xn , u),
.......................................,
.......................................,
dxdtn = fn ( x1 , x2 ,..., xn , u).
Величины x1 , x 2 ,..., xn называются переменными состояния.
Для линейных систем форма Коши задается с помощью матриц
A(t ), B (t) :
x = A(t )x + B (t )u ,
где x – вектор с составляющими x1 , x 2 ,..., xn ; x – символ производ-
ной.
Если система является стационарной, то матрицы A(t ), B (t) яв-ляются постоянными (не зависят от времени).
37
Поскольку с помощью передаточной функции могут быть най-дены частотные и переходные характеристики, то последние спо-собы описания также являются неполными.
3.3. Управляемость, наблюдаемость, стабилизируемость, обнаруживаемость
Рассмотренные ранее примеры заставляют искать особенности движения систем даже при полном их описании [6,7,8,9]. Рассмот-рим линейную динамическую систему второго порядка, которой соответствуют уравнения
dxdt1 = x1 + x2 +u, dxdt2 + x2 = 0.
Ясно, что переменная x2 движется независимо от управляющего
воздействия u . Это пример неуправляемой системы, в которой не удается с помощью управляющего воздействия заставить систему перейти в наперед заданное состояние.
Изучим общий характер движения неуправляемых систем. Для этого используем теорему Гамильтона – Кели: «Всякая квадратная матрица удовлетворяет своему характеристическому уравнению».
Поясним теорему на примере. Пусть имеется матрица |
A = |
1;2 |
|
|
|
. |
|
||
|
|
3;4 |
|
|
Найдем ее характеристическое уравнение | A − λE |= 0 |
(где |
E |
– |
|
единичная матрица), которое в данном примере является квадрат-
ным: λ2 − 5 λ − 2 = 0 . Матричное выражение A 2 − 5 A − 2E можно представлять как матрицу, состоящую из нулевых элементов. Дей-
ствительно, |
7;10 |
|
, поэтому |
A 2 − 5 A − 2 E = 0M , где 0M – |
|
A2 = |
|
|
|||
|
15;22 |
|
|
|
матричный нуль.
38
Пусть стационарная динамическая система находится в состоя-
нии покоя, т.е. x = 0 . Тогда x(t ) = ∫t |
e A ( t −τ ) Bu (τ)dτ . |
|
|
|
|
||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
Раскладывая |
в |
|
ряд |
|
|
экспоненту, |
|
получим |
|
||||
x(t ) = Bα |
0 |
+ ABα |
+ A2Bα |
2 |
+ ... + AnBα |
n |
+ R , где |
R – |
остаток, за- |
|
|||
|
1 |
|
|
|
|
|
|
|
|
|
|||
висящий от старших степеней матрицы A , αk = ∫t |
(t −τ)k |
u (τ)dτ . |
|
||||||||||
|
|
||||||||||||
|
|
|
|
|
|
|
|
0 |
|
k! |
|
|
|
Согласно теореме Гамильтона – Кели степень |
An |
является ли- |
|
||||||||||
нейной комбинацией степеней Ai , |
где i ≤ n −1 . |
Поэтому степени |
|
||||||||||
Aj , где j |
больше n , также являются линейными комбинациями |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
n−1 |
|
||
степеней |
|
Ai . В итоге справедливо равенство x(t ) = ∑ci (t )BAi для |
|
i=0
любого момента времени t . Таким образом, в любой момент вре-мени вектор x(t) является линейной комбинацией столбцов матри-
цы Y = (B ; BA; ...; BAn−1 ). Размерность пространства, в котором на-
ходится вектор x , совпадает с числом линейно независимых столбцов матрицы Y (ее называют матрицей управляемости), т.е. с ее рангом. Если матрица Y является вырожденной, то размерность рассматриваемого пространства меньше – n и вектор состояния не может находиться в любой наперед заданной точке n – мерного пространства. Если неуправляемая часть движения (оно является свободным, связанным с начальными условиями) стремится к ну-лю, то система называется стабилизируемой.
Аналогичные представления рассматриваются и при наблюде-нии за переменными состояния. В приведенном примере с электри-ческой печкой не наблюдалась переменная состояния на выходе корректирующего звена. Система называется наблюдаемой, если при наблюдаемых (измеряемых) сигналах и известном управляю-щем воздействии можно, не решая дифференциальных уравнений, восстановить в каждый момент времени значения всех переменных состояния.
39
Рассмотрим пример. Пусть система описывается дифференци-альными уравнениями:
x1 = a11 x1 + a12 x2 +b1u, x2 = a 21 x1 + a 22 x2 +b2 u,
и измеряется величина y = c1 x1 + c2 x2 . Попробуем, не решая диффе-
ренциальных уравнений, восстановить значения переменных со-стояния. Для этого продифференцируем наблюдаемый сигнал и получим систему
y = c1 x1 + c2 x2 , y = c1 x1 + c2 x2 ,
которую с учетом дифференциальных уравнений можно записать в виде
c1 x1 + c2 x2 = y,
( c1 a11 + c2 a21 ) x1 + ( c1 a12 + c2 a22 ) = y − ( c1b1 + c2 b2 ) u.
Решить последнюю систему относительно x1 |
и x2 |
можно, если |
|
|||
ее определитель отличен от нуля. Матрица N этой системы может |
|
|||||
быть записана следующим образом: |
|
T |
|
, где |
cT – строка |
|
N = c |
T |
|
|
|||
|
|
|
|
|
|
|
|
c |
|
A |
|
|
|
( c1 , c2 ) . Поэтому, если матрица N вырождена, то восстановить x1
x2 невозможно.
общем случае, когда линейная динамическая система имеет порядок n , восстановление всех переменных состояния возможно,
|
|
|
T |
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
если матрица |
c T A |
|
является невырожденной. Матрицу |
N |
|
|||
N = |
... |
|
|
|
||||
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
||
|
|
|
T |
A |
n−1 |
|
|
|
|
c |
|
|
|
|
|
называют матрицей наблюдаемости, которая используется при из-мерении одного сигнала y . Сигнал y вовсе не обязательно совпа-
дает с регулируемым выходным сигналом. Для обеспечения на-блюдаемости системы следует правильно выбирать измеряемые
40
величины. Если ненаблюдаемая часть движения системы стремится к нулю, то систему называют обнаруживаемой.
Упражнение. Показать, что две различные системы, описывае-
мые уравнениями x = |
A1 x + B1u и |
x = A2 x + B2 u , где |
A1 |
0;1 |
|
, |
|
|||||||||
= |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
−1; |
−2 |
|
|
B1 |
0 |
|
, |
A2 |
|
−2;1 |
и |
B2 |
1 |
, |
в первой системе наблюдаемый |
|
||||
= |
|
= |
|
= |
|
|||||||||||
|
1 |
|
|
|
|
−1;0 |
|
|
1 |
|
|
|
|
|
|
|
выходной сигнал y1 = x1 + x2 , а во второй – y 2 = x1 , имеют одинако-
вые передаточные функции, но одна из них неуправляемая, другая – ненаблюдаемая.
Построить структурные схемы систем и убедиться в их разли-чии.