Модели динамики в пространстве состояний

Для описания динамики стационарных систем в классической теории управления используются дифференциальные уравнения вида:

(1.1)

где константы;

y, U- выходные и входные переменные соответственно.

Другой формой описания системы управления является использование передаточных функций :

(1.2)

где Y(s),U(s) - изображения по Лапласу выходной и входной координаты соответственно.

Однако поведение системы во времени характеризуется не только выходными величинами, но и выходными величинами всех составляющих систем элементов. Комплекс таких величин, полностью характеризующих состояние, порождает пространство состояний системы, т.е. характеризует поведение во времени всех динамических элементов, составляющих систему.

Пространством состояний называется метрическое пространство, каждый элемент которого полностью определяет состояние рассматриваемой системы. Состояние понимается в данный момент времени, т.е. как мгновенное состояние. В большинстве случаев достаточно использовать евклидово пространство состояний с евклидовой метрикой.

Метрическим пространством называется множество X, в котором задано расстояние между каждыми двумя элементами xÎX и yÎY в виде действительной функции r(x,y), удовлетворяющей трем аксиомам:

1. r(x,y)=0 тогда и только тогда, когда x=y;

2. r(x,y)=r(y,x) (аксиома симметрии);

3. r(x,z)£r(x,y)+r(y,z) (аксиома треугольника).

Пусть n - порядок системы, равный числу динамических элементов, тогда движение системы характеризуется траекторией движения n – мерного вектора в пространстве . Таким образом, за координаты вектора состояния принимаются выходные сигналы динамических объектов, обладающих свойством накапливать кинетическую или потенциальную энергию, вещество, положение и т.д. В качестве примера таких координат могут служить угол поворота вала электродвигателя, напряжение на конденсаторе и др. Многие элементы не являются идеальными интеграторами и должны рассматриваться с учетом их диссипативных свойств.

Примерами таких элементов могут служить электрическая цепь с индуктивностью и омическим сопротивлением, механическая система с массой и сопротивлением в виде вязких сил трения и др.

При этом дифференциальные уравнения системы представляются в нормальной форме Коши. Если возмущающие воздействия являются детерминированными, то систему можно представить в виде известных функций времени и уравнение указанной системы записать в виде

(1.3)

где пространство управлений.

В более сложном и общем случае рассматриваемая область пространства состояний разделена на подобласти, для каждой из которых задана своя непрерывная функция, удовлетворяющая указанному условию в области определения. Поведение системы на границах подобластей доопределяется, т.е. определяется дополнительными условиями, например, условиями непрерывности.

Если система описывается уравнением (1.3), то процесс в этой системе описывается решением уравнения (1.3) при определенном начальном условии и определенном управлении .

Таким образом, системой в нормальной форме Коши называют систему дифференциальных уравнений первого порядка, разрешенных относительно производных. В общем случае стационарную систему, описываемую линейными дифференциальными уравнениями, можно представить в виде

(1.4)

где постоянные коэффициенты;

возмущающие воздействия.

В матричной форме указанную систему записывают в следующем виде

(1.5)

Матрицы-столбцы называются векторами. Вектор X называют фазовым вектором или вектором состояния, а его координаты x₁,..., x_n - фазовыми координатами. Вектор U называют вектором управления или просто управлением, а его координаты u₁ ,..., u_m - параметрами управления.

Вектор F называют вектором возмущения или просто возмущением, а его j-ая координата - j-м возмущающим воздействием.

Таким образом, движение системы, описываемой уравнением n–го порядка, характеризуется траекторией движения n - мерного вектора в пространстве .

Для нестационарных систем уравнение состояния содержит матрицы как функции времени

(1.6)

При наличии аддитивного белого шума уравнение состояния имеет вид

(1.7)

где V₀(t) - случайный процесс в виде белого шума:

(1.8)

Q₀(t) - матричная функция времени, называемая матрицей интенсивностей.

В случае Q₀=const она называется матрицей спектральных плотностей белого шума.

1.1.2.Примеры вывода уравнений состояния

Уравнение состояния электродвигателя манипулятора

Рассмотрим электродвигатель постоянного тока, используемый в качестве исполнительного элемента в манипуляторах.

Уравнение электрической цепи:

(1.9)

где U - напряжение сети; e_g - ЭДС электродвигателя, i _я - ток якоря, L_я - индуктивность якоря.

Уравнение механической цепи:

(1.10)

(1.11)

где - момент на валу,

M_c - момент сопротивления,

J - момент инерции,

W - угловая скорость вращения электродвигателя,

b - коэффициент.

Момент электродвигателя:

(1.12)

где K_m - коэффициент.

ЭДС электродвигателя:

(1.13)

где K_e - коэффициент.

После подстановки (1.11-1.13) в (1.9) и (1.10) получим:

, (1.14)

. (1.15)

Из (1.14) и (1.15) следует:

(1.16)

где

Уравнение состояния электродвигателя:

. (1.17)

Если обозначить

получим

(1.18)

где .

Структурная схема электродвигателя приведена на рис. 1.1.

Уравнение движения манипулятора

Пусть задан манипулятор с цилиндрической рабочей зоной (рис.1.2).

Для вывода уравнения воспользуемся уравнениями Лагранжа второго рода:

(1.19)

где_i - обобщенные координаты,

Q_i - обобщенные силы,

E_k, E_n - кинетическая и потенциальная энергии системы.

Рассмотрим три основные степени подвижности по координатам r,j,z.

Размеры ориентирующего узла и объекта в схвате считаем существенно малыми по сравнению с остальными линейными размерами манипулятора. Тогда ориентирующий узел можно представить как твердое тело, что существенно упрощает задачу. Считаем все элементы абсолютно твердыми телами. Рука считается однородным стержнем длиной l и массой m_1. Схват с грузом представляет собой точечную массу m. Масса частей, участвующих только в поступательном движении по координатам z и во вращательном движении по координатам j, равна m₂. Соответствующий момент инерции, приведенный к вертикальной оси, равен J.

Приведенный момент инерции частей, участвующих во вращательном движении, равен J₁. Вращающиеся части приводов по координатам r и z учитываются дополнительными приведенными массами и .

Кинетическая энергия манипулятора представляет сумму кинетических энергий отдельных частей.

(1.20)

Потенциальная энергия

(1.21)

где - ускорение свободного падения.

Выбирая обобщенные координаты g₁=r, g₂=j, g­₃=z и производя необходимые преобразования в соответствии с (1.19),получим уравнение движения в матричной форме:

(1.22)

где

вектор ускорений,

вектор, обусловленный силами инерции от взаимного влияния движений по координатам r и j,

- компоненты вектора, обусловленные соответственно центробежной и кориолисовой силами инерции,

вектор сил тяжести,

вектор обобщенных сил,

F_r, F_z – силы, развиваемые приводами по координатам r и z,

M - момент, развиваемый приводом по координате j.

Уравнение состояния интегро-диффенцирующего звена

Принципиальная схема приведена на рис 1.3.а

 

В качестве переменных состояния выбираем напряжение на конденсаторах, характеризующее накопленное количество электричества.

Получим

.

Так как

то после подстановки и соответствующих преобразований получим такое уравнение состояния:

. (1.23)

Выходная величина y=U₂=U₁-U_C1.

Структурная схема приведена на рис 1.3б.

Уравнение состояния смесительного бака

На рис. 1.4. приведена схема смесительного бака, в котором необходимо поддерживать заданный уровень и концентрацию раствора.

Уравнение баланса масс для бака имеет вид:

Мгновенный расход выходного потока F(t) зависит от высоты жидкости h:

где k - константа.

При постоянной площади сечения бака S

Тогда уравнения баланса примут вид:

, (1.24)

(1.25)

Установившееся состояние при F₁₀=const, F₂₀=const, F₀=const, V₀=const, C₀=const:

При небольших отклонениях от установившегося состояния:

где m₁, m₂ - рассматриваются как входные переменные, а x₁, x₂ - переменные состояния.

Линеаризация (1.24) и (1.25) приводит к уравнениям

(1.26)

(1.27)

Введем параметр - время заполнения бака:

. (1.28)

Исключив x(t) из (1.27), получим линеаризованное дифференциальное уравнение:

(1.29)

где

Если определить выходные переменные в виде:

то уравнение (1.29) можно дополнить линеаризованным уравнением выходной переменной

где

1.1.3.Общее решение дифференциальных уравнений

Решение дифференциальных уравнений в матричной форме выполняется с использованием фундаментальной матрицы.

Уравнение (1.5) называется неоднородным.

Однородное уравнение имеет вид

(1.30)

Пусть образует n линейно - независимых решений этого уравнения.

Любую такую систему называют фундаментальной системой решений уравнения (1.30). Матрица, составленная из этой системы решений, называется фундаментальной матрицей уравнений (1.6) и (1.30):

(1.31)

Если при t=t₀ фундаментальная матрица обращается в единичную, то она называется нормированной. Используя произвольную фундаментальную матрицу Ф(t),нормированную матрицу X(t,t₀) можно представить в виде:

где Е - единичная матрица.

Нормированная фундаментальная матрица обладает следующими свойствами:

.

Рассмотрим методику решения дифференциальных уравнений в матричной форме. Пусть динамика системы описывается уравнением в скалярной форме:

(1.32)

где x - неизвестный скаляр, а - константа.

Решение этого уравнения известно:

. (1.33)

Разложение (1.33) в ряд Тейлора приводит к следующему выражению

. (1.34)

Для определения характеристического уравнения системы (1.32) используются преобразование Лапласа уравнения (1.32):

(1.35)

где p - оператор Лапласа.

По аналогии с рассмотренным случаем находят характеристическое уравнение системы при векторных переменных. Допустим, в однородном матричном уравнении начальные условия имеют вид X(0)=E.

Разложим X(t) в ряд Тейлора в окрестности X(0):

(1.36)

Очевидно

С учетом (1.15) получим

(1.37)

Следовательно, решение однородного уравнения (1.6) имеет вид:

(1.38)

Характеристическое уравнение получим путем преобразования по Лапласу (1.38).

(1.39)

откуда характеристическое уравнение:

(1.40)

С помощью нормированной фундаментальной матрицы решение неоднородного уравнения (1.5) при всех t и X(t₀)=X_­0 можно представить в виде соотношения

(1.41)

которое называется формулой Коши.

Если матрица А постоянная, то матрица X(t,t₀) зависит только от разности (t-t₀) и имеет вид

. (1.42)

Матричная функция называется экспоненциальной матрицей или матричным экспоненциалом и определяется суммой ряда:

(1.43)

Уравнение, сопряженное с (1.30), имеет вид:

(1.44)

Если Z(t,t₀) – нормированная фундаментальная матрица этого уравнения, т.е. то формулу Коши представляют в виде

(1.45)

При решении задач оптимального управления удобно пользоваться также операторным методом решения систем дифференциальных уравнений.

Этот метод состоит в том, что посредством интегрального преобразования Лапласа от дифференциального уравнения переходят к вспомогательному алгебраическому уравнению. Затем находят решение преобразованного уравнения, после чего при помощи обратного преобразования получают решения дифференциального уравнения:

.

Пусть требуется найти решение линейных дифференциальных уравнений с постоянными коэффициентами.

(1.46)

где .

, (1.47)

где постоянные.

Преобразование Лапласа приводит к уравнению

(1.48)

где , или сокращено

. (1.49)

Решением полученного уравнения имеет вид

(1.50)

Обратное преобразование только для второго слагаемого дает решение дифференциального уравнения с нулевыми начальными условиями.

Пример 1.1.

Тогда

Решение вспомогательного уравнения имеет вид

Для вычисления обратного преобразования раскладываем правую часть на простейшие дроби:

При помощи таблицы преобразования Лапласа по этому изображению находим оригинал:

.

Операторный метод удобно использовать также для решения линейных систем с постоянными коэффициентами. Пусть дана задача Коши:

,

.

при условии det(a_ik)¹0. Преобразование Лапласа приводит к системе

.

Из этой вспомогательной системы находим Y_k (p), производим обратное преобразование и получим решение.

Пример 1.2. Найти общее решение системы:

.

Пусть y₁ (0)=с, y₂(0)=с₂.

После преобразования получим

откуда

Выполнив обратное преобразование, получим общее решение

1.1.4.Преобразование дифференциальных уравнений

Преобразование дифференциальных уравнений к нормальной форме

Дифференциальное уравнение, разрешимое относительно старшей производной, всегда можно привести к нормальной системе. Рассмотрим линейное уравнение n - го порядка с постоянными коэффициентами:

, (1.51)

или в символической форме

(1.52)

Такая система эквивалентна линейной системе n уравнений первого порядка, собственные значения которых равны корням полинома:

. (1.53)

Пусть собственные значения системы - различные действительные числа.

Предположим y(t) - решение уравнения (1.51) и определим n - мерную функцию x(t), примем

(1.54)

x_i(t) - компоненты вектора x(t).

Отсюда имеем:

(1.55)

В векторной форме эта система имеет вид:

. (1.56)

или сокращенно:

(1.57)

С другой стороны, если x(t) – решение уравнения (1.57), то функция y=x₁(t) – решение уравнения (1.51), т.е. полученная система n уравнений первого порядка эквивалентна одному уравнению n - го порядка.

Чтобы убедиться в том, что собственные значения матрицы А системы (1.56) являются корнями полинома (1.53), найдем собственные значения матрицы А:

. (1.58)

Детерминант матрицы (1.58) после раскрытия по минорам

(1.59)

или

(1.60)

где является собственным значением матрицы А в том и только в том случае, когда l - корень полинома (1.53), т.е.

Следовательно, уравнения (1.51) и (1.56) эквивалентны.

Пример 1.3. Рассмотрим систему:

Полагая, получим систему:

Матрица .

Собственные значения матрицы - корни полинома:

Допустим собственные значения матрицы А - действительные и различные числа l₁,l₂,....,l_n и матрица собственных значений:

. (1.61)

Существует невырожденная матрица Р, такая, что:

(1.62)

где Р - матрица Вандермонда

. (1.63)

Детерминант матрицы Р равен

т.е. Р - невырожденная матрица, если все l_i различны.

Представление в пространстве состояний при содержании только полюсов.

В этом случае дифференциальное уравнение системы в символической форме имеет вид

(1.64)

где - скалярный выход, - скалярный вход.

Передаточная функция

(1.65)

В общем случае корни характеристического уравнения комплексные будут полюсами , поэтому

Начальные условия системы (всего ) - . По аналогии с (1.54) определяем переменных

(1.66)

Из (1.66) следует, что

(1.67)

Из уравнения (1.64) следует

(1.68)

Следовательно

(1.69)

откуда вытекает, что удовлетворяет системе дифференциальных уравнений, которые в матричной форме имеют вид

, (1.70)

или

(1.71)

Вектор представляет собой переменную состояния системы (1.64), поскольку знание и полностью определяют и выход . Если и известны, то единственное решение уравнения (1.71) с в качестве исходной точки и управляющей функцией выражается соотношением:

(1.72)

где - фундаментальная матрица уравнения (1.66).

Следовательно, полностью определяется и . Уравнениями состояния системы являются:

(1.73)

Для выяснения физического смысла понятия состояния системы рассмотрим модель решения уравнения (1.70) (рис. 1.5).

Как видно по рис.1.5, переменные состояние X₁(t),X₂(t),…,X_n(t) соответствуют выходным сигналам интеграторов. Постоянные коэффициенты а₀, а₁,…, а_n являются коэффициентами усиления в цепях каналов обратной связи.

Так как собственные значения матрицы А являются корнями полинома

. (1.74)

То они совпадают с полосами передаточной функции W(p) (1.65).

Собственные значения являются инвариантным свойством линейного преобразования и матрица, связанная с линейным преобразованием, зависит от базисов векторных пространств, поэтому полюса передаточной функции (собственные значения матрицы) являются инвариантными величинами. При этом переменные состояния и матрицы могут быть различными.

Случай, когда передаточная функция системы имеет полюса и нули

Пусть объект с одним входом и одним выходом описывается линейным дифференциальным уравнением с постоянными коэффициентами

(1.75)

Передаточная функция

(1.76)

Здесь принято, что число нулей равно числу полюсов. В противном случае можно выбирать тогда будем иметь m нулей и n полюсов.

Обозначим через p₁, p₂, …, p_n комплексные корни полинома знаменателя и s₁, s₂, …,s_n - комплексные корни полинома числителя, т.е. p₁, p₂, …, p_n – полюса, а s₁, s₂, …,s_n – нули W(p). Тогда можно написать

(1.77)

Для того, чтобы решить дифференциальное уравнение (1.75), необходимо знать 2n начальных условий: y(o), y(o),…,y^(n-1)(0) и U(0), U(0),…, U^(n-1)(0). Найдем вектор состояния. Пусть X(t) – n – мерный вектор с компонентами x₁(t), x₂(t),…, x_n(t), заданными соотношениями

(1.78)

где h₀, h₁,…, h_n-1 - постоянные, которые необходимо найти.

Соотношения (1.78) можно записать в виде

(1.79)

Найдем дифференциальное уравнение, которому удовлетворяет x_i(t). Из (1.78) заключаем, что

(1.80)

Аналогично, можно получить

(1.81)

так как

(1.82)

Дифференцирование выражения (1.79) для i=n дает:

(1.83)

Дифференциальное уравнение (1.79) дает соотношение

(1.84)

Из выражения (1.82) следует, что

(1.85)

Следовательно

(1.86)

Путем непосредственных вычислений можно получить

(1.87)

Подставив это выражение в уравнение (1.86) и сгруппировав члены, получим

(1.88)

Выберем k таким образом, чтобы не зависела от всех производных U(t). Это можно сделать, положив n₀=b₀,

для (1.89)

Подставив выражение (1.89) в уравнение (1.88), получим дифференциальное уравнение

(1.90)

где

(1.91)

Полученные дифференциальные уравнения относительно x_i(t) можно записать в векторной форме

, (1.92)

или сокращенно

(1.93)

В развернутой форме вычисление h_i, i=1,n имеет вид

h₀=b_0,

h₁=b_n-1-h₀a_n-1,

h₂=b_n-2-h₀a_n-2-h₁a_n-1 и т.д. (1.94)

Вектор X(t) представляет вектор состояния системы (1.75), так как знание X(t₀) и U(t₀,t) полностью определяет решение дифференциального уравнения (1.93) на интервале [t₀,t]. Поскольку выход y(t) равен

Y(t)=x₁(t)+h₀U(t), (1.95)

то y(t₀,t) полностью определяется величинами X(t₀) и U(t₀,t).

Случай, когда n>m

Предположим, что в уравнении (1.75)

b_n=b_n-1=…=b_m+1=0; b_m¹0, (1.96)

(1.97)

Выполним преобразования

(1.98)

где

, (1.99)

можно рассматривать как вектор состояния системы с передаточной функцией (1.97). Дифференциальное уравнение имеет вид

. (1.100)

Пример 1.4.

Дана система

для которой n=3, a₀=0, a₁=2, a₂=3, b₃=0, b₂=1, b₁=7, b₀=12.

Из уравнения (1.91) находим h₀=0, h₁=1, h₂=4, h₃=-2.

Фазовые координаты определяются соотношением (1.90):

Дифференциальное уравнение

.

Преобразование уравнений состояния

При анализе и синтезе систем управления возникает необходимость преобразования линейных дифференциальных уравнений с постоянными коэффициентами.

Пусть уравнение системы

(1.101)

Определим новую переменную состояния

, (1.102)

где T-невырожденная nxn матрица. Подставляя (1.102) в (1.101), получим

откуда

(1.103)

Тогда

(1.104)

где

Матрицы A и Q подобны, следовательно, имеют одинаковые собственные значения, поэтому вектор Z(t) можно рассматривать как вектор состояния системы.

. (1.105)

Таким образом, по известному вектору Z(t) можно найти вектор X(t).

Заменой переменных можно преобразовать модель системы управления. Это возможно, если собственными значениями матрицы системы (1.70) являются различные действительные числа. Как было показано выше (см 1.62), в этом случае можно отыскать невырожденную матрицу Р, которая является матрицей Вандермонда, такую, что

. (1.106)

Если определить вектор состояния, положив

.

То Z(t) будет удовлетворять дифференциальному уравнению

(1.107)

Обозначим элементы P^-1 через r_ij, i,j=1,2,…,n. Так как вектор B задан в виде B=[0 0 … 0 b₀]^T, то уравнение (1.107) можно записать:

(1.108)

где -собственные значения матрицы системы, r_in-элементы последнего (n-го) столбца матрицы P^-1.

На рис. 1.6 приведена модель системы (1.107). Компоненты вектора состояния Z_i(t) представляют собой выходные сигналы n интеграторов. Выход y(t)=x₁(t) равен сумме Z_i(t), т.е.

(1.109)

так как первая строка матрицы Вандермонда состоит целиком из единиц.

Числа в (1.108) можно вычислить непосредственно, проинвертировав матрицу P. Более рациональным способом можно получить эти числа следующим способом. Так как полюса W(p) расположены при p=l_i, можно написать

(1.110)

где представляет собой вычет W(p) в i-м полюсе l_i. Так как =l_i различны, имеем

.

Предположим, что в начале система находится в покое (начальные условия нулевые) и что приложен вход U(t), тогда

(1.111)

При этих же условиях из выражения (1.107) находим

(1.112)

Так как , то

(1.113)

Сравнивая (1.111) с (1.113), получим

(1.114)

т.е. элементы последнего столбца матрицы равны вычетам W(p).

Преобразование подобия позволило заменить модель системы (рис.1.5) моделью (рис.1.6), коэффициенты усиления которой равны вычетам передаточной функции W(p) в полюсах соответственно.

Устойчивость систем управления в пространстве состояний

В общем случае система устойчива, если евклидова норма ||x(t)|| остается ограниченной при для любого решения x(t) системы. Систему называют абсолютно устойчивой в том и только в том случае, если

, (1.115)

для любого решения системы x(t). Если система неустойчива, то существует такое решение , при котором

(1.116)

Пусть собственные значения матрицы А будут

где - действительные числа, .

Тогда критерий устойчивости гласит: система устойчива, если

- для всех i в случае, когда корни – действительные числа;

- для всех k в случае кратных корней вида .

Система абсолютно устойчива, если для любого i.

Линеаризация уравнений систем

В большинстве случаев системы описываются нелинейными дифференциальными уравнениями. В общем случае уравнение состояния имеет вид

. (1.117)

Нелинейную систему можно линеаризировать. При этом принимают следующие допущения:

- отклонения переменных достаточно ма