Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
Скачать 25.93 Mb.
|
W p W p W p W p ± = (5.59) или ) ( ) ( ) ( ) ( ) ( ) ( ) ( 2 1 2 1 2 1 p R p R p Q p Q p Q p R p W p ± = (5.60) Здесь знак минус относится к положительной, а знак плюс — к отрицательной обратной связи. Обратные связи будут рассмотрены подробно в главе, посвященной методам улучшения динамических свойств системы автоматического регулирования. При использовании динамических звеньев обычно наиболее просто находится передаточная функция разомкнутой системы (рис. 5.1). Затем по формулам, приведенным в § 5.2, легко находятся все уравнения системы автоматического регулирования. При анализе системы автоматического регулирования необходимо составить ее так называемую структурную схему, представляющую собой совокупность динамических звеньев со связями между звеньями. Такая структурная схема часто является весьма простой и ее составление не представляет особого труда. Однако в некоторых случаях составление структурной схемы сопряжено с большими трудностями и может быть сделано только на основании детального анализа исходных дифференциальных уравнений системы регулирования. В этом случае структурная схема не облегчает нахождения основных уравнений системы однако ив этом случае она остается весьма ценной, так как на ней в наглядной форме представлены все узлы исследуемой системы и все существующие между ними связи. Это может оказаться полезным во всех дальнейших исследованиях. На рис. 5.6 в качестве примера приведена структурная схема разомкнутой системы регулирования в том случае, когда цепь регулирования представляет собой простую цепь последовательно включенных звеньев. В этом случае передаточная функция разомкнутой системы W(p) = W 0 (р (р (р (р) (5.61) Здесь W 0 (р, W 1 (р, W 2 (р) и W 3 (р) представляют собой заданные передаточные функции объекта регулирования и отдельных звеньев, входящих в систему регулирования. Нетрудно видеть, что для нахождения передаточной функции разомкнутой системы можно разомкнуть систему необязательно так, как это показано на риса в произвольном месте. На рис. 5.7 изображен более сложный пример системы автоматического регулирования. Передаточная функция разомкнутой системы в этом случае (5.62) Ив этом случае для нахождения передаточной функции разомкнутой системы можно разомкнуть систему в другом месте, например в точках ас или А. Для рассмотренных на рис. 5.6 и 5.7 систем, зная передаточную функцию разомкнутой системы W(p), легко найти по формулами) дифференциальные уравнения для регулируемой величины и ошибки, записанные в символической форме где g(t) задающее воздействие. На рис. 5.8 изображена структурная схема системы стабилизации. В этом случае задающее воздействие g(t) = const представляет собой настройку регулятора. Определив передаточную функцию разомкнутой системы (5.63) можно по формулами) получить символические записи дифференциальных уравнений для регулируемой величины и ошибки где f(t) — возмущение, действующее на объекта (р) — передаточная функция регулируемого объекта по возмущению. В тех случаях, когда структурная схема оказывается сложной и содержит многоразличных перекрестных связей, можно попытаться ее упростить и свести к простейшему виду, например к изображенной на рис. 5.6. Преобразование структурных схем линейных систем делается на основе некоторых правил, которые даны в табл. 5.2. На рис. 5.9 изображены этапы упрощения сложной структурной схемы на основе приведенных выше правил. При упрощении введены дополнительные передаточные функции, определяемые выражениями Полученная в результате преобразования схема (рис. 5.9, в) уже относится к простейшим. Использование графов. Подобно структурным схемам графы прохождения сигналов используются для наглядного изображения математических зависимостей в системах регулирования. Графом (рис. 5.10, б) называется множество вершин и ребер. Каждому ребру соответствуют две вершины — начало и конец ребра. Вершине и ребру могут быть сопоставлены или некоторые величины, или операторы, например передаточные функции. Основные свойства графов прохождения сигналов следующие. 1. Каждая вершина, отмеченная на графе кружком или точкой, соответствует некоторой переменной (координате) рассматриваемой системы. 2. Каждое ребро графа, изображаемое в виде линии со стрелкой, указывающей направление прохождения сигнала, имеет вершину-начало (входную величину) и вершину-конец (выходную величину. Если из вершины выходит несколько ребер, то все они имеют одинаковую входную величину. 3. Выходная величина ребра получается как результат преобразования, осуществляемого соответствующим ребру оператором, входной величины ребра. 4. Если к одной вершине подходит несколько ребер, то величина, соответствующая этой вершине, получается алгебраическим суммированием выходных величин этих ребер. Между структурной схемой и графом прохождения сигналов имеется прямое соответствие прямоугольник структурной схемы соответствует ребру, а линия передачи сигнала — вершине графа. На рис. 5.10 для сравнения изображены одновременно структурная схема (аи граф прохождения сигналов (б) одной и той же системы. Правила преобразования графов подобны правилам преобразования структурных схем линейных систем. Эти правила изображены на рис. 5.11 в виде исходных (первый столбец) и эквивалентных (второй столбец) схем. В дальнейшем изложении будут использоваться более удобные структурные схемы. § 5.5. Многомерные системы регулирования К многомерным относятся системы управления и регулирования, имеющие несколько регулируемых величину. Это имеет место во многих современных сложных системах. К ним относятся, например, системы регулирования напряжения и частоты синхронных генераторов, системы управления подвижных объектов, многие системы регулирования технологических процессов и др. Многомерная система предполагает наличие многомерного объекта управления рис. 5.12), который характеризуется существованием нескольких видов (точек приложения управляющих и возмущающих воздействий) и нескольких выходов, определяемых регулируемыми величинами. Многомерный объект описывается системой уравнений, которую удобно представлять в матричной форме. Введем одностолбцовую мерную матрицу регулируемых величин (5.64) одностолбцовую мерную матрицу управляющих величин (5.65) и одностолбцовую мерную матрицу возмущающих воздействий (5.66) Здесь штрихом обозначена операция транспонирования матрицы. Если регулируемые величины имеют одинаковую физическую размерность и могут трактоваться как проекции некоторого вектора на оси координат, матрица-столбец может отождествляться с этим вектором. Тогда можно говорить о векторе регулируемых величин. Если регулируемые величины имеют разную физическую размерность, то переход от матрицы-столбца к вектору в принципе может быть сделан ив этом случае, если ввести в матрицу-столбец весовые коэффициенты, уравнивающие размерности отдельных составляющих. Однако такой переход не является единственным, а имеет бесчисленное количество вариантов. Аналогичным образом при равенстве физических размерностей отдельных составляющих матриц-столбцов управляющих величин и возмущений может быть введен вектор управления и вектор возмущения. При разных физических размерностях отдельных составляющих матриц-столбцов переход к вектору возможен, ноне будет единственным. Линеаризованные уравнения движения многомерного объекта могут быть записаны в матричном виде (5.67) Здесь введена квадратная матрица операторных коэффициентов (5.68) и прямоугольные матрицы операторных коэффициентов (5.69) (5.70) Если в выражениях (5.64) — (5.70) перейти к изображениям Лапласа при нулевых начальных условиях, то матричное уравнение (5.67) может быть записано для изображений в следующем виде (5.71) Здесь р, р) и р) — матрицы-столбцы изображений регулируемых величин, управляющих величин и возмущений. В уравнение (5.71) входят также квадратная матрица n m ij p q p Q × = ) ( ) ( и прямоугольные матрицы k m ij p r p R × = ) ( ) ( и Если матрица р) неособая, те. определитель 0 ) ( ≠ p Q , то, умножив левую и правую части (5.71) слева на обратную матрицу Q -1 (р, получим (5.72) Здесь введены матрицы передаточных функций объекта для управляющих величин (5.73) и для возмущении (5.74) В (5,74) символом m m ij p Q p Q × = ) ( ) ( обозначена матрица, присоединенная для матрицы р, а Q ij - алгебраическое дополнение определителя ) ( Формулы (5.72) — (5.74) позволяют получить связь между регулируемыми величинами и управляющими и возмущающими воздействиями. Так, например, если m= 3, k = 2 и l = 0, то из (5.72) и (5.73) можно получить для изображений (5.75) Если в матрице передаточных функций объекта (5.73) или (5.74) для каждого элемента матрицы (частной передаточной функции) найти обратное преобразование Лапласа (оригинал, то будет получена так называемая матрица Коши (матрица весовых функций. Запишем ее, например, для управляющих воздействий (5.76) Если в момент времени t= 0 на все входы поступают управляющие воздействия u i (t), где i= 1, 2, . . ., k, то изменение и регулируемой величины может быть записано посредством интеграла Дюамеля — Карсона (4.9) на основании принципа суперпозиции На рис. 5.13 изображена условная структурная схема замкнутой многомерной системы регулирования. На схеме все указанные символы соответствуют матрицам g(t) — задающих воздействий, у) — регулируемых величин, х (t) — ошибок для каждой регулируемой величины, u(t) — управляющих воздействий, f(t) — возмущений, о (р) — передаточных функций для управлений, W f (р) — передаточных функций для возмущений. Кроме того, введена прямоугольная матрица передаточных функций регулирующего устройства m k ij рег p k p W × = ) ( ) ( которая определяет используемые законы регулирования. Она дает связь между изображениями управляющих величин и ошибок (5.77) Уравнения многомерной системы (рис. 5.13) могут быть получены действиями, аналогичными одномерному случаю (§ 5.2). Матрица передаточных функций разомкнутой по всем каналам системы (5.78) Характеристическая матрица системы представляет собой квадратную матрицу размером т х т (5.79) Здесь I — единичная матрица размером тахт, те. квадратная матрица, у которой все элементы главной диагонали равны единице, а остальные — нулю. Характеристическое уравнение системы получается приравниванием нулю определителя характеристической матрицы (5.80) Заметим, что в случае, когда многомерная система представляет совокупность т независимых одномерных систем, характеристическая матрица будет диагональной и определитель системы тогда равен произведению частных определителей каждой из систем, те, В этом случае общее характеристическое уравнение распадается на т независимых характеристических уравнений 0 ) ( = p D i , i=1,2,... m. Матрицы передаточных функций замкнутой системы, замкнутой системы по ошибке и замкнутой системы по возмущениям при условии, что матрица р) неособая, что означает независимость исходных дифференциальных уравнений, могут быть определены из выражений (5.81) (5.82) (5.83) Здесь m m ij p D p D × = ) ( ) ( — матрица, присоединенная для матрицы р, ар алгебраическое дополнение определителя Полученные выражения для матриц передаточных функций замкнутой системы позволяют использовать формулы, аналогичные формулам § 5.2, но записанные уже для матриц-столбцов ошибок и регулируемых величин. Так, например, для матрицы изображений ошибок имеем (5.84) На рис. 5.14 изображены для иллюстрации некоторые структурные схемы двумерных систем регулирования. Схема на риса соответствует так называемому сепаратному регулированию объекта с двумя входами и двумя выходами. Матрица передаточных функций регулирующего устройства в этом случае получается диагональной. Матрицу изображений управляющих величин для этого случая можно представить в виде (5.85) Схемы на рис. 5.14, б ив соответствуют комбинированному регулированию. В этом случае (5.86) Исходные дифференциальные уравнения многомерной системы регулирования могут быть также представлены в форме Коши в матричной записи (5.87) В этих выражениях — матрица-столбец фазовых координат системы, n— порядок дифференциального уравнения, —матрица-столбец регулируемых величин, — матрица-столбец управляющих величин, — матрица-столбец возмущающих и задающих воздействий, — квадратная матрица коэффициентов, и — прямоугольные матрицы коэффициентов. Величины x i (i = 1, 2, . . ., n) представляют собой некоторые абстрактные величины, задание которых полностью определяет текущее состояние системы- Эти величины называются фазовыми координатами системы. Состояние системы может быть также отождествлено с положением изображающей точки в мерном пространстве, которое носит название пространства состояния. При переходе к изображениями совместном решении система уравнений может быть приведена, например, к виду (5.84). Характеристическое уравнение, соответствующее системе (5.87), имеет вид (5.88) где I — единичная матрица n x n. Векторная запись исходных уравнений. Введем в рассмотрение мерное векторное пространство состояния, которое определяется базисом (набором векторов) так, что с матрицей-столбцом фазовых координат может быть отождествлен вектор состояния (5.89) Длины векторов базиса играют роль весовых коэффициентов в переходе от матрицы-столбца фазовых координат к вектору состояния. Заметим, что в общем случае, когда рассматриваются абсолютные, а не относительные значения фазовых координат, их физические размерности не совпадают и длины векторов базиса не могут считаться единичными. Аналогичным образом может быть введено векторное пространство управления, возмущения и выходных величин. При введении векторов исходные уравнения системы могут быть записаны в векторной форме (5.90) Фазовые координаты, а также составляющие управления, возмущения и выходных координат могут быть получены как проекции соответствующих векторов на оси, определяющие выбранные векторные пространства. При использовании относительных (безразмерных) величин в качестве базиса может приниматься совокупность ортогональных векторов единичной длины, те. обычное n- мерное евклидово пространство. § 5.6. Управляемость и наблюдаемость Рассмотрим мерное пространство состояния X, в котором каждому состоянию системы соответствует некоторое положение изображающей точки, определяемое значениями фазовых координат x i (i = 1, . . ., n). Пусть в пространстве состояния X заданы два множества и Рассматриваемая система будет управляемой, если существует такое управление , определенное наконечном интервале времени , которое переводит изображающую точку в пространстве X из подобласти Г в подобласть Г 2 Можно сузить определение управляемости и понимать под ней возможность перевода изображающей точки из любой области пространства состояния X в начало координат. Система будет полностью управляемой, если каждое состояние управляемо в этом смысле. От пространства состояния X перейдем к другому пространству посредством неособого преобразования , причем , где R — матрица коэффициентов n x n. Тогда вместо (5.87) будем иметь (5.91) Здесь использованы преобразованные матрицы коэффициентов , , Введение новых фазовых координат посредством неособого преобразования приводит к эквивалентным системам различной структуры. При некотором преобразовании может оказаться, что часть управляющих величин не входит в некоторые дифференциальные уравнения (5.91) или часть фазовых координат не участвует в формировании выхода у. В первом случае система будет не полностью управляемой, а во втором — не полностью наблюдаемой. В случае не полностью управляемой системы ее исходные уравнения (5.87) могут быть представлены в виде (5.92) Ох. Это иллюстрирует рис. 5.15. Набор фазовых координат х соответствует управляемой части фазовых координата набор ж — неуправляемой части. Р. Калманом [50] был доказан критерий управляемости, который гласит, что размерность V управляемой части системы, то есть порядок первой группы уравнений (5.91), совпадает с рангом матрицы (5.93) При v = n система полностью управляема, при 0 < v < n не полностью управляема и при v = 0 полностью неуправляема. На риса изображен простейший пример. Если рассматривать выходную величину у (t) при ненулевых начальных условиях, то можно записать (5.94) где определяются начальными условиями до приложения входного сигнала u 1 (t), ау в (t) — вынужденная составляющая. Система устойчива при аи с. Если начальные условия до приложения u 1 (t) были нулевыми, то поведение системы может быть рассчитано по передаточной функции (5.95) В этом случае по интегралу Дюамеля — Карсона (5.96) Как следует из выражений (5.95) и (5.96), система во втором случае-описывается дифференциальным уравнением не третьего, а второго порядка. Система будет устойчивой даже при а. < 0. Рассмотренная система будет не полностью управляемой. В ней оказывается, что n = 3, а v = 2. При введении второй составляющей управления u 2 (t) система оказывается полностью управляемой, и ей будет соответствовать матрица-строка передаточных функций по управлению В случае не полностью наблюдаемой системы ее уравнения могут быть, представлены в виде (5.97), Эти уравнения отличаются от (5.87) тем, что фазовые координаты группы x 2 не входят нив выражения для у и u, нив первое уравнение, куда входят только фазовые координаты группы x 1 . Группа фазовых координат x 2 относится к ненаблюдаемым. Это иллюстрирует рис. 5.17. Р. Калманом [50] показано, что порядок первой группы уравнений v совпадает с рангом матрицы (5.98) При v = n система полностью наблюдаема, при 0 < v , управляемую и наблюдаемую часть x 2 , неуправляемую и ненаблюдаемую часть x 3 и неуправляемую, но наблюдаемую часть Исходные уравнения системы (5.87) в этом случае можно для самого, общего случая записать следующим образом (5.99) Левая часть характеристического уравнения системы в этом случае содержит четыре сомножителя Управляемость и наблюдаемость системы в изложенном смысле не всегда совпадает с практическими представлениями. Даже если какая-либо фазовая координата и может быть вычислена по доступным для измерения выходным величинам, обработка измеренных величин может быть, во-первых, сложной и, во-вторых, она может быть затруднена наличием помех. Поэтому практически наблюдаемыми координатами обычно считаются те из них, которые могут быть непосредственно измерены датчиками различных типов. |