Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
![]()
|
(7.32) на примере преобразования Лапласа. Перейдем в левой и правой частях (7.32) к изображениям Лапласа. При этом оператор дифференцирования в полиномах D (р) и N (р) заменяется на комплексную величину , а вместо оригиналов хи) появляются их изображения X (р) и F (р. В результате получаем где D 0 (р) обозначает сумму всех членов, содержащих начальные условия. Отсюда находится изображение искомой величины (7.33) Последнее выражение требует некоторых пояснений в связи с различными возможными трактовками понятия начальных условий. Интегральное преобразование Лапласа (7.17) следует, вообще говоря, записать в более строгом виде (при заменена р (7.34) Это дает возможность введения двух несколько отличающихся понятий преобразования Лапласа (и соответственно преобразования Карсона — Хевисайда). 1. Преобразование Лапласа по начальным условиям справа. Если в выражении (7.34) нижний предел интегрирования стремится к нулю, оставаясь положительным (а >0), тов изображении производной (7.26) следует брать начальные условия прите. для момента времени, который будет сразу после приложения к системе внешних воздействий. В этом случае Для использования последней формулы необходимо знание начальных условий справа, что оказывается не всегда удобными требует расчета по формулам § 7.3. Заметим, что даже в тех случаях, когда до приложения воздействия система находилась в покое, начальные условия справа могут быть ненулевыми и полином D 0 (р, как правило, отличен от нуля. Кроме того, если рассматриваемая функция времени f(t) имеет при t = 0 особенности типа -функции, то это обстоятельство не будет учтено в найденном изображении. Так, например, изображение самой -функции и ее производных оказывается при этом равным нулю 2. Преобразование Лапласа по начальным условиям слева. Если в формуле (7.34) нижний предел интегрирования стремится к нулю, оставаясь отрицательным (а < 0), тов выражении для изображения производной (7.26) следует брать начальные условия прите. для момента времени, который будет непосредственно предшествовать моменту приложения воздействия. Такие начальные условия называются также предначалъными. В этом случае Расчет получается более простым, так как предначальные условия должны быть известны всегда и никаких дополнительных операций здесь не требуется. В частном случае, когда до приложения воздействия система находилась в покое, предначальные условия нулевые и выражение (7.33) приобретает вид (7.35) Только это выражение и позволяет строго сформулировать понятие передаточной функции ТУ(р) как отношение изображений входной и выходной величин при нулевых предначальных условиях. Кроме того, преобразование Лапласа в случае, когда нижний предел интегрирования стремится к нулю, оставаясь отрицательным (а < 0), позволяет учитывать наличие в рассматриваемой функции при 2 = 0 особенностей типа -функции. Так, например, изображение единичной -функции оказывается равным единице а изображение ее производной го порядка Влияние особенностей f (t) и ее первых m производных, где m — порядок полинома р) , на изображение р) в этом случае и проявляется в виде автоматического учета начальных условий, которые будут иметь место справа (при t = + 0) в самом изображении N (р) Р (р) без введения дополнительного члена D 0 (р) при нулевых предначальных условиях или без его изменения при ненулевых предначальных условиях. В связи с этим функция иногда называется также функцией начальных условий. В дальнейшем изложении под преобразованием Лапласа будет пониматься именно этот случай (а < 0). Зная изображение искомой величины X (р) в виде (7.33) или (7.35), можно найти оригинал х (t). Это и будет решением исходного дифференциального уравнения (7.32). Для отыскания оригиналах) по его изображению X (р) можно пользоваться таблицами изображений и существующими теоремами, в частности теоремой разложения, которая устанавливает следующее. Если изображение Лапласа имеет вид отношения двух многочленов (7.36) то при отсутствии нулевых корней знаменателя (7.37) где р k — некратные корни знаменателя (7.36). Если знаменатель изображения Лапласа (7.36) имеет нулевой корень (p 0 = 0), то изображение надо представить в виде (7.38) Тогда оригинал может быть найден по формуле (7.39) Аналогичным образом теорема разложения может быть записана и для преобразований Карсона — Хевисайда. Так, например, если изображение искомой величины может быть представлено в виде отношения двух полиномов (7.40) то при отсутствии нулевых и кратных корней знаменателя оригинал будет определяться выражением (7.41) Это выражение полностью совпадает с формулой (7.39), так как изображения Лапласа (7.38) и Карсона — Хевисайда (7.40) отличаются на множитель р. Использование изображений часто называют также операторным методом, хотя в действительности операторному методу, разработанному Хеви-сайдом, оказывается полностью аналогичным использование только преобразований Карсона — Хевисайда (7.21) и (7.22). Метод использования изображений обладает тем преимуществом, что в нем полностью сохраняется лишь одна операция — вычисление корней характеристического уравнения (знаменателя изображения. Что касается определения произвольных постоянных интегрирования, то эта операция отпадает, потому что начальные условия автоматически учитываются в процессе решения с самого начала (при нахождении изображения искомой величины. Поэтому этот метод оказывается удобными его часто применяют в задачах теории регулирования. Практически важной для отыскания оригинала решения является еще теорема свертывания. Она гласит следующее. Если изображение представляет собой произведение (7.42) то оригинал выражается формулой (7.43) где т представляет собой вспомогательное время интегрирования. В частности, пусть для некоторой системы с передаточной функцией р) известна реакция на единичную импульсную функцию , представляющую собой функцию веса и связанную р) преобразованием -Лапласа Если на вход этой системы поступает некоторая функция времени f(t), изображение которой р, то изображение выходной величины будет Тогда функция времени на выходе может быть найдена по интегралу свертывания (7.43), который совпадает с интегралом Дюамеля (4.9): (7.44) л ^ о о Если входная функция определена только для положительного времени прикладывается на вход в момент времени t = 0), то функция отлична от нуля только при . В этом случае верхний предел интеграла в формуле (7.44) может быть заменен на бесконечность иона приобретает вид (7.44') § 7.5. Использование вещественных частотных характеристик Опишем метод приближенного построения кривой переходного процесса в автоматической системе (при воздействиях в виде скачка и импульса) по заданной вещественной частотной характеристике замкнутой системы, разработанный В. В. Солодовниковым в 1948 году [121]. Этот способ полезен тогда, когда расчет системы ведется с самого начала, частотными методами. Он совершенно необходим, если известны уравнения не всех звеньев системы, а часть из них задается экспериментально снятыми частотными характеристиками. На основании интеграла Фурье (7.16) оригинал искомой величины может быть представлен в виде (7.45) где изображение Фурье искомой функции времени ха частотное изображение искомой величины, полученное из изображения Карсона — Хевисайда (р) подстановкой р = jω. Однако использовать интегральную зависимость (7.45) можно только в том случае, когда все полюсы функции X (jω) лежат в левой полуплоскости. Тогда интегрирование может вестись по мнимой оси. Это значит, что для преобразования Лапласа (7.18) абсцисса абсолютной сходимости сир. В действительности изображение Фурье X (jω) даже для устойчивой системы, когда все полюсы передаточной функции системы лежат в левой полуплоскости, может иметь полюсы на мнимой оси за счет входного воздействия. Так, например, пусть передаточная функция системы имеет вид причем аи Полюсы этой передаточной функции лежат в левой полуплоскости. Если на вход поступает сигнал типа единичной ступенчатой функции x 1 (t) = 1(t), изображение которого по Лапласу равно , то изображение выходной величины будет Это изображение имеет однократный полюс вначале координат (р = 0). Если на вход системы поступает сигнал типа , изображение которого , то изображение выходной величины будет иметь вначале координат двукратный полюс р = р = 0). В связи с этим для использования интегральной зависимости (7.45) необходимо отделить от изображения Фурье искомой функции времени члены, содержащие полюсы на мнимой оси. Рассмотрим частный случай, когда изображение Карсона — Хевисайда не имеет полюсов на мнимой оси. К этому случаю сводится например, задача нахождения переходной функции в устойчивой системе, если даны ее передаточная функция Фр, не имеющая полюсов на мнимой оси, и входное воздействие типа единичной ступенчатой функции . Тогда изображение по Лапласу выходной величины будет и соответственно Тогда оказывается, что частотное изображение совпадает с частотной передаточной функцией замкнутой системы Фа вещественная и мнимая части в формуле (7.46) совпадают с вещественной Р (ω) и мнимой S(ω) частотными характеристиками замкнутой системы. К аналогичному результату можно прийти, если рассматривать реакцию системы на скачок внешнего возмущения. Тогда вещественная и мнимая части в формуле (7.46) будут совпадать с вещественной и мнимой частями частотной передаточной функции по возмущению К тому же частному случаю могут сводиться и другие задачи исследования переходных процессов в системах регулирования, например нахождение' ошибки системы при приложении скачкообразного внешнего возмущения, нахождение функции веса системы и др. В этом случае существует ограниченное установившееся значение искомой функции времени х (t), которое можно получить, подставляя в значение р = 0. Учитывая, что , получаем . Тогда подынтегральная функция (7.45) может иметь однократный полюс вначале координат. Его можно устранить, рассматривая не саму величину ха разность , которой соответствует разность изображений . В результате приходим к следующей интегральной зависимости (7.47) Используем формулу Эйлера Подставляя последнее выражение в (7.47), используя формулу (7.46) и отбрасывая мнимую часть, которая должна быть равной нулю, так как функциях) является, конечно, вещественной, получаем (7.48) Подынтегральное выражение представляет -собой четную функцию частоты. Поэтому интегрирование по всем частотам можно заменить интегрированием только по положительным частотам, а затем удвоить результат. Так как тов результате имеем (7.49) Если принять нулевые начальные условия, то до приложения внешнего воздействия при t<0) . Заменив в (7.49) время t на —t, получим (7.50) Совместное решение (7.49) и (7.50) дает два выражения для нахождения искомой функции времени (7.51) (7.52) причем Таким образом, можно отыскать оригинал х (t) по известной вещественной или известной мнимой частям частотного изображения . Обычно для этих целей используется вещественная часть изображения Если входное воздействие представляет собой единичный скачок, то, как указывалось выше, частотное изображение совпадает с частотной передаточной функцией замкнутой системы . Тогда в формулы (7.51) и (7.52) будут входить вещественная и мнимая части частотной передаточной функции замкнутой системы . Следовательно, в этом случае для построения переходного процесса, который будет представлять собой переходную функцию системы h(t), необходимо в формуле (7.52) положить — вещественная характеристика системы. В результате получим (7.53) Аналогичным образом, при нахождении реакции системы на единичный скачок возмущающего воздействия необходимо использовать вещественную часть частотной передаточной функции по возмущению. В дальнейшем изложении будем иметь ввиду случай, определяемый формулой (7.53), хотя методика построения переходного процесса остается единой и для общего случая (7.52). Интегрирование выражения (7.53) представляет большие трудности. Поэтому обычно используется приближенное решение задачи. Для этой цели вводится понятие типовой единичной трапецеидальной вещественной характеристики, (рис. 7.3). Единичная трапеция имеет высоту, равную единице и частоту среза , также равную единице, точнее, 1 сек -1 Единичная трапеция характеризуется частотой излома, которая может быть задана в виде коэффициента наклона трапеции Для единичных трапеций с различными коэффициентами наклона по выражению (7.53) может быть вычислен оригинал, те. функция времени. Эта функция получила название функции. В настоящее время составлены подробные таблицы функции для различных коэффициентов наклона, лежащих в пределах 0 < χ < 1 (см. приложение 1). По такой таблице для каждого коэффициента наклона единичной трапеции может быть построена функция времени — безразмерное время, соответствующее единичной трапецеидальной характеристике. Метод построения кривой переходного процесса заключается в том, что построенную вещественную характеристику исследуемой системы (рис. 7.4) разбивают наряд трапеций, заменяя приближенно кривые линии прямолинейными отрезками так, чтобы при сложении ординат всех трапеций получилась исходная характеристика. Затем для каждой трапеции определяется коэффициент наклона. При известном коэффициенте наклона по таблицам могут быть построены функции для каждой трапеции. Кривая переходного процесса может быть получена суммированием построенных h- функций с учетом правил масштабов. Правила масштабов заключаются в следующем. 1. Перед сложением ординаты каждой функции необходимо умножить на высоту соответствующей трапеции (см. рис. 7.4), так как функция построена для трапеции, имеющей единичную высоту. При этом необходимо учитывать знак высоты, считая высоту положительной для трапеций, расположенных выше абсцисс. 3. Перед сложением необходимо изменить масштаб времени каждой функции, так как функции построены для единичной трапеции, имеющей частоту среза = 1 сек -1 Изменение масштаба времени делается в соответствии с теоремой подобия (табл. 7.2). Действительное время равно времени приведенному в таблице функций, деленному на частоту среза соответствующей трапецеидальной характеристики При нахождении реакции системы на единичную импульсную входную функцию, те. функции веса , можно пользоваться общей формулой (7.52). При этом должно быть вещественной частью частотного изображения искомой функции , где 1 представляет собой изображение единичной импульсной функции . Однако можно преобразовать формулу (7.53) так, что и при нахождении функции веса можно будет исходить из вещественной частотной характеристики замкнутой системы Р (ω). Для этой цели продифференцируем выражение (7.53) повремени) Если разбить исходную вещественную характеристику на трапецеидальные характеристики (рис. 7.4), то аналогично построению переходной функции выражение (7.54) можно представить в виде где n — число трапеций, на которые разбита вещественная характеристика Р. Можно показать, что это выражение приводится к виду (7.55) где введены обозначения Следовательно, в данном случае искомая функция времени приближенно определяется простым подсчетом ее ординат по формуле (7.55) для разных t и последующим построением по точкам. Для облегчения подсчетов можно воспользоваться готовыми таблицами значений зш а/а, которые имеются в справочниках. В заключение заметим, что при построении кривой переходного процесса по трапецеидальным частотным характеристикам наибольшие ошибки получаются в начальной части кривой, так как отбрасываемый хвост вещественной частотной характеристики замкнутой системы влияет главным образом именно на начальную часть кривой переходного процесса. Кроме изложенного здесь частотного метода В. В. Солодовникова существует еще предложенный А. А. Вороновым [28] аналогичный способ построения кривых переходного процесса по треугольным частотным характеристикам. § 7.6. Использование вычислительных машин За последнее время для исследования систем автоматического регулирования и управления ив частности, для построения переходных процессов стали широко применяться вычислительные машины непрерывного и дискретного действия. Наибольшее применение находят вычислительные машины непрерывного действия, относящиеся к классу моделирующих установок электронного или электромеханического типа. Удобство моделирующих вычислительных машин заключается в том, что физическому процессу, протекающему в исследуемой системе регулирования, соответствует протекание в вычислительной машине (модели) некоторого другого аналогового процесса, описываемого теми же дифференциальными уравнениями, что и исходный процесс. Это позволяет изучать процессы в системах регулирования наиболее наглядно, так как каждой обобщенной координате в исследуемой системе соответствует некоторая переменная в вычислительной машине, например электрическое напряжение, ток (в электронной модели) или угол поворота (в электромеханической модели. Моделирующие вычислительные машины позволяют моделировать как всю систему в целом, таки отдельные ее части. Так, например, часто вычислительная машина используется для моделирования объекта регулирования, например самолета, корабля, паровой турбины, двигателя внутреннего сгорания и т. пасам регулятор может быть реальным. При сопряжении реального регулятора с объектом, в качестве которого выступает модель, получается замкнутая система регулирования, которая может быть исследована еще до того, как будет построен сам объект. Вычислительные машины целесообразно использовать для исследования обыкновенных линейных систем в тех случаях, когда последние описываются дифференциальными уравнениями сравнительно высокого порядка и их аналитическое исследование становится малоэффективным. Однако наибольшее значение имеют вычислительные машины при исследовании линейных систем с переменными параметрами и нелинейных систем, поскольку для этих случаев пока еще мало разработано приемлемых для практики методов, а иногда аналитические методы вообще отсутствуют. Точность моделирующих вычислительных машин обычно не превосходит нескольких процентов. В большинстве случаев этого оказывается достаточно для целей практики. Получение точности в десятые доли процента и выше связано со значительным |