Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
![]()
|
t ). В связи с этим функцию веса можно изобразить в виде некоторой поверхности (рис. 13.2). Эта поверхность переходит в плоскость t0 ϑ при t< ϑ . Границе перехода поверхности в плоскость соответствует биссектриса t = ϑ . Это обстоятельство объясняется тем, что в реальных системах реакция не может появиться ранее приложения на входе системы импульса. Поэтому при t= ϑ функция веса должна быть тождественно равна нулю. Сечение поверхности весовой функции вертикальной плоскостью, параллельной оси t риса, дает весовую функцию для фиксированного момента приложения единичного импульса на входе системы ( ϑ = const). Эта функция называется нормальной весовой функцией системы с переменными параметрами (13.4) Она является параметрической функцией, так как в нее входит фиксированный параметр ϑ = const;. Нормальная весовая функция может быть сделана зависящей от аргумента ϑ τ − = подстановкой. В результате получаем функцию (13.5) Сечение поверхности весовой функции вертикальной плоскостью, параллельной оси и, дает кривую, образованную ординатами семейства нормальных весовых функций для фиксированного значения времени t=const (рис. 13.2, б. Эта кривая может быть получена путем обработки семейства нормальных весовых функций, построенных для различных моментов приложения единичного входного импульса (рис. 13.3). Получающуюся зависимость будем называть сопряженной функцией веса ) , ( ϑ ϑ ω − t , t=const. (13.6) Она также является параметрической функцией, так как содержит параметр t=const. Сопряженная функция веса является функцией смещения , но может быть представлена также как функция аргумента ϑ θ − = рис. 13.2, б, называемого реверс-смещением, поскольку θ отсчитывается от точки ϑ =t в сторону, противоположную смещению ϑ . Это осуществляется подстановкой в сопряженную весовую функцию значения θ ϑ − = t при t=const. В результате получаем (13.7) Проиллюстрируем все сказанное примером. Пусть функция веса системы с переменными параметрами имеет вид Зафиксировав смещение и положив, например, ϑ =ϑ 0 = const, получаем нормальную функцию веса или в другом виде, при переходе к аргументу ϑ τ − = t : Зафиксировав текущее время и положив, например, t=t 0 =const, получаем сопряженную функцию веса Перейдя к реверс-смещению ϑ θ − = t , имеем Заметим, что в системах с постоянными параметрами весовая функция является функцией только времени ϑ θ − = t и не зависит от момента приложения ϑ входного импульса. Рельеф функции веса (рис. 13.2) в этом случае получается цилиндрическим, а оба рассмотренных выше сечения (риса и б) совпадают по форме и отличаются только знаками аргументов. При переходе к реверс-смещению получаем полное совпадение двух функций веса — нормальной и сопряженной Пусть на систему (13.1) с функцией веса действует входной сигнал f(t). Элементарная реакция на выходе системы в произвольный момент времени t ϑ ≥ будет (13.8) Полный сигнал на выходе линейной, системы определяется как суперпозиция элементарных реакций интегрированием (13.8) в пределах от 0 до t: (13.9) Так как при >t функция веса равна нулю, то выражение (13.9) можно также записать в виде (13.10) Из двух последних выражений видно, что в интегральном уравнении связи между входной и выходной величинами используется сопряженная функция веса (13.6), те. разрез рельефа функции веса (рис. 13.2, б) вдоль аргумента ϑ . Если использовать реверс-смещение ϑ θ − = t , то интегральная связь (13.9) может быть представлена в виде интеграла свертки (13.11) Как уже отмечалось, в случае постоянства параметров системы функция веса зависит только от времени (t- ϑ В этом случае формула (13.11) переходит в интеграл свертки (7.44) § 13.2. Нахождение функции веса и построение переходных процессов Функция веса системы с переменными параметрами является исчерпывающей характеристикой этой системы, и нахождение ее важно последующим соображениям. Функция веса характеризует протекание временных процессов в системе регулирования, и по ее виду можно судить о качестве регулирования, аналогично тому, как это делалось для систем с постоянными параметрами (§ 8.4). По имеющейся функции веса можно определить время протекания переходного процесса, как характеристику быстродействия, и склонность системы к колебаниям. Кроме того, по имеющейся функции веса можно строить процесс на выходе системы регулирования при заданных входных воздействиях, не производя при этом каждый раз полного решения исходного уравнения (13.1). В соответствии с формулами (13.9) и (13.11) для этой цели необходимо иметь сопряженные функции веса. Ввиду сложности проблемы существующие методы позволяют пока решать задачу нахождения функции веса в численном виде. Только для систем регулирования, описываемых дифференциальными уравнениями первого и иногда второго порядков, удается решать задачу в общем виде. Поэтому в некоторых случаях приходится сложную систему с переменными параметрами приближенно сводить к более простой системе, движение которой описывается уравнением не выше второго порядка. Следует заметить, что большинство систем регулирования с переменными параметрами относится к так называемым квазистационарным системам, или системам, параметры которых меняются сравнительно медленно. В подобных системах коэффициенты дифференциального уравнения (13.1) мало меняются в течение времени переходного процесса, определяемого временем затухания нормальной функции веса. Дифференциальное уравнение первого порядка. В некоторых случаях для оценки вида переходных процессов системы с переменными параметрами ее уравнение приближенно можно свести к дифференциальному уравнению первого порядка (13.12) Это уравнение имеет аналитическое решение (13.13) где а С — постоянная интегрирования. Пусть, например, имеется уравнение (13.14) Определим для него семейство переходных характеристик ) , ( ) , ( υ υ υ t h t h = − . Для единичной ступенчатой функции при 0 ≠ υ уравнение можно записать в следующем виде Приведем его к виду (13.12): Далее получаем На основании формулы (13.13) получаем При нулевых начальных условиях (для υ = t ) должно быть 0 ) , 0 ( = υ h . Отсюда определяется постоянная интегрирования Окончательно получаем Дифференцируя последнее выражение по υ , можно получить функцию веса или вином виде Для дифференциального уравнения (13.12) можно сразу найти функцию веса из общего решения (13.13), если положить в (13.12) входной сигнал равным единичному смещенному импульсу ) ( ) ( υ δ − = t t Q . Проделав необходимые выкладки, получаем (13.15) где Распространим этот результат на более общий случай записи дифференциального уравнения в виде (13.16) Приведем его к виду (13.12): Положив ) ( ) ( υ δ − = t t f , получаем для функции веса решение в виде (13.18) где Рассмотрим снова в качестве примера уравнение (13.14). Приведем его к виду (13.17): Обратившись к формуле (13.18), находим и функцию веса что совпадает с полученным ранее выражением. Дифференциальное уравнение второго порядка. Рассмотрим случай, когда дифференциальное уравнение (13.1) сводится к уравнению второго порядка (13.19) При помощи подстановки (13.20) это уравнение приводится к виду (13.21) Здесь введено обозначение (13.22) При действии единичного импульса ( ) υ δ − = t t f ) ( для уравнения (13.21) получится решение ) , ( υ υ − = t z u , которое связано с весовой функцией ) , ( υ υ ω − = t z исходного уравнения (13.19) на основании формулы (13.20) соотношением Если же положить ( ) υ δ − = t t f ) ( , то для уравнения (13.21) будет получена весовая функция ( ) υ δ − = t t r ) ( , которая на основании (13.9) связана с решением ( ) υ δ − = t t z зависимостью Эта зависимость на основании свойства дельта-функции может быть представлена в виде (13.24) В результате из (13.23) и (13.24) получаем (13.25) Таким образом, для отыскания функции веса ) , ( υ υ ω − t необходимо предварительно решить уравнение (13.21), которое приобретает вид (13.26) с нулевыми начальными условиями 0 ) ( = t u и 0 ) ( = t u при υ = t . Полученную при решении весовую функцию ( ) υ υ , ) ( − = t r t u необходимо затем подставить в (13.25) и найти Решение уравнения (13.26) может быть произведено при помощи использования функций Бесселя [118]. Для этого функция F(t) должна быть аппроксимирована отрезками прямых линий, уравнение которых сводится к виду t b a i i + . Однако это решение является сравнительно сложным. Ограничимся рассмотрением так называемого аппроксимирующего решения, которое может применяться, если функция F(t) мало изменяется относительно своего среднего большого значения F cp (t) (рис. 13.4). Это решение называется аппроксимацией Бриллуина — Вентцеля — Крамера [118]. Рассмотрим однородное дифференциальное уравнение (13.27) Предположим теперь, что для некоторого однородного дифференциального уравнения второго порядка получено частное решение (13.28) где (13.29) Найдем дифференциальное уравнение, которому удовлетворяет решение (13.28). Продифференцировав его дважды и исключив промежуточные переменные, получаем (13.30) Сравнивая (13.30) и (13.27), видим, что выражение (13.28) будет частным решением уравнения (13.27), если выполняется тождество (13.31) Решение уравнения (13.31) и отыскание функции N(t) является сложной задачей вследствие наличия нелинейностей в (13.31). Однако может быть найдено приближенное решение (13.31) в виде ряда, если удовлетворяются неравенства Тогда решение (13.21) можно представить в виде (13.32) Подставляя этот ряд в (13.31), получаем формулы для определения членов ряда (13.33) Часто можно ограничиться только первым членом ряда (13.22), что будет справедливым, если функция F(t) изменяется медленно, оставаясь в среднем большой (риc. 13.4). Тогда (13.34) При выполнении условия F(t)> 0 в качестве второго частного решения можно взять комплексно-сопряженную величину (13.29) (13.35) Тогда можно показать, что решение уравнения (13.26) будет (13.36) или, после подстановки (13.28) и (13.35), (13.37) В предельном случае постоянства параметров const t F = Ω = 2 ) ( . Тогда t t S Ω = ) ( и υ υ Ω = ) ( S . В результате из формулы (13.37) можно получить функцию веса консервативного звена Для исходного дифференциального уравнения (13.19) на основании (13.25) и (13.37) получаем искомую функцию веса (13.38) Критерием медленности изменения функции F(t) и, следовательно, применимости полученного выражения может служить неравенство (13.39) которое получается из (13.31) и (13.34). Метод последовательных приближений. Рассмотрим уравнение (13.1): Ограничиваясь случаем квазистационарных систем и полагая, что коэффициенты меняются медленно, найдем функцию веса для этого уравнения. Переменные коэффициенты в левой части исходного уравнения представим в виде суммы постоянной и изменяющейся частей (13.40) где ) ( 0 υ i i a a = — переменный коэффициент, зафиксированный для момента- приложения входной величины Тогда исходное дифференциальное уравнение (13.1) можно представить в виде (13.41) где (13.42). (13.43), Поскольку мы предположили, что коэффициенты ) (t a i меняются медленно, то функция у мала по сравнению с левой частью (13.41). Эту функцию можно рассматривать как возмущение, и тогда к уравнению (13.41) можно применить метод последовательных приближений. В уравнении (13.41) можно перейти к изображениям по Лапласу. Тогда получим (13.44) Здесь введено обозначение (13-45) Решение уравнения (13.41) или (13.44) можно записать в виде ряда (13.46) Для получения первого приближениях зафиксируем переменные коэффициенты ) ( ) ( υ i i a t a = . Тогда первое приближение может быть найдено как решение дифференциального уравнения (13.47) Решение этого уравнения можно получить, используя обычные методы см. главу 7), в том числе путем нахождения оригинала, соответствующего изображению (13.44) при У(р)= 0: (13.48) Для получения второго приближения в правую часть (13.41) или (13.44) подставляется первое приближение х = ха в левую часть -х = х + х. Тогда получается уравнение с фиксированными коэффициентами для определения поправки (13.49) Это уравнение также может быть решено с использованием преобразования Лапласа посредством нахождения оригинала изображения где У 1 (р) — изображение у при подстановке в формулу (13.43) х = х. Повторяя этот процесс многократно, можно найти рекуррентное соотношение для определения го члена ряда (13.46): (13.50) Ряд (13.46) сходится тем быстрее, чем медленнее изменяются коэффициенты Рассмотренный метод может использоваться как для нахождения функции веса и переходной функции, таки для построения переходного процесса при любом известном воздействии f(t). Численно-графический метод. Численно-графический метод ДА. Башкирова [98] разработан также применительно к системам с переменными во времени параметрами, причем можно вводить любое переменное возмущающее или задающее воздействие и произвольные начальные условия. Неоднородные уравнения первого порядка с переменными коэффициентами. Пусть требуется построить решение уравнения с начальным условием х = х 0 при t = 0. Разделив его на уравнение к виду (13.51) где Уравнение (13.51) можно решать графически, если считать Т постоянными равным ) 2 ( t t T ∆ + внутри каждого интервала времени ) , ( t t t ∆ + , но различным для разных интервалов. Формула для решения в этом случае будет а процесс построения сводится к следующему. Наносим заданные кривые f(t) и Т (рис. 13.5). Из точки Е кривой f(t), взятой в середине первого интервала t ∆ , откладываем по горизонтали отрезок ) 2 ( t T EM ∆ = , величина которого берется равной ординате точки Н заданной кривой Т, те. тоже в середине первого интервала t ∆ . Полученная точка М соединяется прямой линией с заданной начальной точкой процесса А. В результате получается новая точка В искомой кривой х (t). Затем аналогично берется ордината точки I, откладывается в виде отрезка Ни проводится прямая В, дающая новую точку С решениях, и т. д. Неоднородные уравнения второго порядка с переменными коэффициентами. Требуется построить решение уравнения которое можно записать также в виде (13.52) где с начальными условиями х = х, 0 x x = при t = 0. Если в правой части <13.52) имеется операторное выражение, то предварительно производим вычисление правой части и сводим ее к F(t). Если обозначить х = Т (х, то уравнение (13.52) разобьется на два (13.53) где (13.54) а начальные условия будут Формулы для решения уравнений (13.53) согласно [64, 78] будут (13.55) где причем во второй из формул (13.55) значения x ∆ берутся со сдвигом на t ∆ вправо по сравнению с Отсюда вытекает следующее построение. Наносим заданные кривые T 1 (t) и T2(t), а также кривую Т, ординаты которой определяются по второй из формул (13.54). Они показаны на графике (риса. На другом графике наносим заданное f(t) (рис. 13.6, 6). На основании заданных начальных условий (см. выше) наносим на последнем графике точки хи в середине первого интервала (как в § 7.6) еще точку Ас ординатой (7.76), те. Из точки Ев середине первого интервала t ∆ на кривой f(t) откладываем вниз отрезок вниз, когда он положителен, и вверх, когда он отрицателен. При этом величина берется как ордината уже имеющейся точки А, величина ) 2 ( 3 t T ∆ берется из графика Та величинах из заданных начальных условий. Из полученной точки Е откладываем горизонтальный отрезок размер которого берется из графика T 1 (t). Точку М соединяем прямой линией сточкой, что дает новую точку Н кривой х при t = t ∆ . Из точки Н откладываем вниз отрезок равный ординате точки А. Из точки Н проводим горизонтальный отрезок размер которого берется из графика Т. Точку К соединяем сточкой, что дает новую точку В искомой кривой х (t) в середине второго интервала t ∆ . Опишем еще второй шаг интегрирования. Из точки F 1 кривой F(t) в середине второго интервала t ∆ откладывается вниз отрезок где х в — ордината точки B, полученной выше α tg — тангенс угла наклона прямой КА, проведенной ранее (он дает требуемое значение х. Откладываем отрезок и проводим прямую N 1 H 1 , получая при этом новую точку I 1 кривой х. Из точки I 1 откладываем вниз а затем вправо после чего проводим прямую В. Это дает новую точку С искомой кривой хит. д. Все описанные построения можно заменить числовыми расчетами. |