Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975. Теория систем автоматического управления. В.А. Бесекерский, Е.П.. 3 Примеры непрерывных автоматических систем
 Скачать 25.93 Mb.
|
|
k M j ) — значение частотной передаточной функции, полученное из Ф (z, ε ) подстановкой Аналогичным образом по передаточной функции Ф (z, ε ) может быть получена для установившегося режима выходная величина у [n, ε ]. Более простой метод заключается в следующем. Рассмотрим, например, задающее воздействие g[n], представляющее собой периодическую последовательность, изображение которой (15.113) где G 0 (z) — изображение g[n] на интервале 0 - М. Пусть рассматриваемая последовательность действует на входе системы с передаточной функцией Ф. Тогда изображение выходной величины можно представить в виде суммы изображений переходной составляющей У (z) и установившегося периодического режима У (z). Первая составляющая определяется полюсами функции Фи стечением времени затухает, так как система предполагается устойчивой. Периодическая составляющая на выходе может быть представлена в виде где У (z) — изображение у [n] на интервале 0 - М в установившемся режиме, которое и является искомой величиной, коэффициенты а, . . ., а должны быть определены при разложении на сумму дробно-рациональных функций Y 0 (z) и У (z). Для этой цели могут использоваться известные методы, например теорема разложения. Так, если степень У (z) равна степени У (z) и У (z) = У (z), то (15.178) где ) ,..., 1 ( M q z q = — корни уравнениям. Однако пользоваться этой формулой при М > 2 практически неудобно. При М >1 удобнее найти переходную составляющую У (z), а затем периодическую У (z) = У (z) — У (z). Далее можно найти (15.179) Так как полюсы Ф (z) известны, то отыскание переходной составляющей не представляет труда. Так, например, если степень числителя У (z) меньше степени знаменателя и полюсы Ф (z) не кратные, то где z i (i =1,2, ...,l) — полюсы Ф. Если степень числителя У (z) одинакова со степенью знаменателя, но числитель имеет общий множитель г, то можно записать Тогда Другие возможные случаи — см. § 15.2, п. 13. Если входное воздействие представляет собой симметричную периодическую последовательность с полупериодом N, изображение которой дается формулой (15.114), то аналогичная зависимость для изображения периодической последовательности выходной величины на интервале 0 - N будет (15.180) где У (z) — переходная составляющая, определяемая полюсами Ф (z). Пример Рассмотрим входную последовательность в виде прямоугольной волны (рис. 15.10, внос полупериодом N = 3 и систему с передаточной функцией Изображение периодической последовательности на входе (15.114) Найдем периодический режим на выходе. В соответствии с (15.180), учитывая, что Ф (z) имеет единственный полюс z 1 = 0, имеем Отсюда следует, что в установившемся периодическом режиме на выходе, если совместить начало положительного полупериода с началом отсчета, будет у [0] = 1 + ау у [2] = 1 - а. В следующем полупериоде будет у [3] = -у [0] = - (1 + ау у [5] = у [1] = -(1 - аи т. д. § 15.5. Случайные процессы в импульсных системах Введем понятие случайной решетчатой функции f[n], которую можно образовать из непрерывной случайной функции f(z) ее дискретизацией. В этом случае она будет определена в дискретные моменты времени t= Т. Будем рассматривать стационарные процессы, когда вероятностные характеристики не зависят от времени. Среднее значение решетчатого случайного стационарного процесса (15.181) или на основании эргодического свойства (15.182) где w(f[n]) — одномерная плотность вероятности. Для центрированных процессов среднее значение равно нулю. Введем понятие корреляционной функции (15.183) Аналогично главе 11 можно сформулировать основные свойства корреляционной функции. 1. Для случая m = 0 (15.184) 2. При m = 0 корреляционная функция достигает наибольшего значения (15.185) 3. Корреляционная функция является четной (15.186) При наличии двух случайных процессов f 1 [n] и f 2 [n] можно ввести понятие взаимной корреляционной функции (10.187) Свойства ее схожи со свойствами взаимной корреляционной функции для непрерывных процессов. Введем понятие спектральной плотности случайного стационарного решетчатого процесса как двустороннего преобразования корреляционной функции (15.188) где Т — нормирующий множитель, равный периоду дискретности, а F(z) представляет собой преобразование корреляционной функции R[m]. Нормирующий множитель Т введен в (15.188) для того, чтобы сделать физическую размерность спектральной плотности дискретного случайного процесса равной размерности спектральной плотности непрерывного процесса и сохранить ее физический смысл. Аналогично непрерывному случаю можно ввести понятие спектральной плотности как функции круговой частоты (15.189) или при учете четности (15.190) Наконец, можно определить спектральную плотность как функцию абсолютной псевдочастоты. Для этого в формуле (15.188) необходимо перейти к преобразованию, используя подстановку (15.163), а затем перейти к псевдочастоте посредством подстановки k T j w 2 = . В результате получим (15.191) Аналогичным образом может быть определена взаимная спектральная плотность двух процессов. Заметим, что все приведенные формулы могут быть записаны и для случая 0 ≠ ε , тогда рассматривается случайная решетчатая функция f[n, ε ], корреляционная функция R[m, ε ], спектральные плотности S (z, ε ), S (w, ε ) и S* ( λ Я, ε ). Основное свойство спектральной плотности, как ив непрерывном случае, заключается в том, что интеграл от нее по всем частотам дает средний квадрат случайной величины. Можно показать [136], что в дискретном случае соответствующая формула имеет вид (15.192) Так как имеют место равенства то формула (15.192) может быть записана в виде (15.193) Выражение (15.193) обычно является более удобным для расчетов по сравнению с (15.192), так как позволяет использовать таблицы интегралов-(см. приложение 2). Типовые случайные стационарные процессы. Если для функции f(z), представляющей собой центрированную помеху, эффективное время корреляции (15.194) меньше периода дискретности, T < ∆ τ , то такой процесс может быть представлен как дискретный белый шум с корреляционной функцией (15.195) где R[0] = D — дисперсия, а δ 0 [m] — единичная решетчатая импульсная функция (15.32), равная единице при m = 0 и равная нулю при 0 ≠ m . Этому белому шуму соответствует спектральная плотность (15.196) Если эффективное время корреляции T > ∆ τ , то корреляционная функция R [m] может быть получена из соответствующей корреляционной функции непрерывного процесса R( τ ) заменой mT = τ . Спектральная плотность может быть получена использованием формул (15.188) — (15491). В табл. 15.2 приведены некоторые типовые дискретные стационарные случайные процессы. Прохождение сигнала через линейную систему. Пусть на входе линейного звена с известной дискретной передаточной функцией W(z) действует случайная функциях, для которой известны корреляционная функция R 1 [m] и спектральная плотность 1 S (w) или S*( λ ). Тогда для выходной, величины х, аналогично непрерывному случаю, можно найти спектральную плотность умножением спектральной плотности входного сигнала на квадрат модуля частотной передаточной функции (15.197) Интегрирование спектральной плотности по всем частотам в соответствии си) позволяет найти средний квадрат выходной величины ] [ 2 2 n x . Это позволяет для замкнутой импульсной системы производить расчеты, аналогичные изложенным в § 11.8. Так, например, пусть в схеме, изображенной на рис. 15.11, на входе действуют полезный сигнал g(t) и помеха n(t), не коррелированные между собой. Обозначим их спектральные плотности ) ( * λ g S и ) ( * λ n S . Тогда спектральная плотность ошибки (15.198) где Фи Ф (j λ ) — частотные передаточные функции замкнутой системы и замкнутой системы по ошибке. Интегрирование (15.198) по всем частотам в соответствии сдает средний квадрат ошибки (15.199) Подобным же образом могут быть найдены расчетные формулы и для других возможных случаев (см. § 11.8). РАЗДЕЛ IV НЕЛИНЕЙНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ ГЛАВА 16 СОСТАВЛЕНИЕ УРАВНЕНИЙ НЕЛИНЕЙНЫХ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ § 16.1. Общие понятия Нелинейной системой автоматического регулирования называется такая система, которая содержит хотя бы одно звено, описываемое нелинейным уравнением. Перечислим виды нелинейных звеньев 1) звено релейного типа (рис. 1.12); 2) звено с кусочно-линейной характеристикой (рис. 1.10, д и др 3) звено с криволинейной характеристикой любого очертания 4) звено, уравнение которого содержит произведение переменных или их производных и другие их комбинации 5) нелинейное звено с запаздыванием, причем запаздывание понимается в смысле § 14,1, а нелинейность может иметь любой вид 6) нелинейное импульсное звено 7) логическое звено 8) звенья, описываемые кусочно-линейными дифференциальными уравнениями, в том числе переменная структура. Различают статические и динамические нелинейности. Первые представляются в виде нелинейных статических характеристика вторые — в виде нелинейных дифференциальных уравнений. Общий метод составления уравнений для нелинейных систем состоит в следующем. Сначала по правилам § как делалось в главе 5, производится линеаризация уравнений всех звеньев системы, для которых это допустимо, кроме существенно нелинейных звеньев (чаще всего одного-двух). Затем составляются уравнения этих последних звеньев со всеми допустимыми упрощениями их характеристик. В результате получается система обыкновенных линейных уравнений, к которым добавляется одно-два (иногда более) нелинейных. В соответствии с этим обобщенную структурную схему любой нелинейной системы автоматического регулирования в случае одного нелинейного звена можно представить в виде риса, где линейная часть может иметь структуру любой сложности (с обратными связями и т. п, как, например, на рис. 16,1, били в. В случае двух нелинейных звеньев могут быть разные комбинации, в зависимости оттого, в какие цепи системы они входят (см, например, рис. 16.2). Часто при исследовании нелинейных систем автоматического регулирования удается выделить нелинейность так, чтобы она описывалась непосредственно зависимостью между выходной и входной величинами (16.1) которая может иметь любую форму (релейного типа, кусочно-линейного или криволинейного. Но иногда, как будет показано в следующих параграфах, не удается этого сделать и приходится исследовать нелинейные дифференциальные зависимости вида (16.3) Встречаются и более сложные случаи, когда обе величины (входная и выходная) оказываются под знаком нелинейной функции раздельно (16.4) или же вместе (16.5) Разделим все нелинейные системы регулирования на два больших класса. 1. К первому классу нелинейных систем отнесем такие, в которых уравнение нелинейного звена приводится к любому из видов (16.1) — (16.3), те. когда под знаком нелинейной функции стоит только входная величина (и ее производные) либо только выходная величина (и ее производные).При этом имеется ввиду, что схема системы в целом может быть приведена к виду рис. 16.1 с одним нелинейным звеном. К этому классу сводится, например, также случай с двумя нелинейными звеньями, указанный на рис. 16.2, в, так как там они могут быть объединены водно нелинейное звено. Сюда же относится и случай, показанный на рис. 16.2, г, где имеются два нелинейных звена (если их уравнения содержат под знаком нелинейности только входную величину х, например, вида (16.1) или (16.2)). 2. Второй класс нелинейных систем включает системы с любым числом нелинейных звеньев, когда под знаки нелинейных функций входят различные переменные, связанные между собой линейной передаточной функцией. Так будет в случае системы с одним нелинейным звеном вида (16.4) или (16.5), а также в системе с двумя нелинейными звеньями (риса или г, если в первом из них под знак нелинейности входит входная величина, а во втором — выходная. Система же рис. 16.2, б относится ко второму классу, если под знаки нелинейностей входят в обоих звеньях либо только входные, либо только выходные величины нелинейных звеньев. Ко второму классу нелинейных систем относятся также системы с двумя и более нелинейностями, в уравнениях которых под знаки нелинейных функций входят разные переменные, связанные между собой нелинейными дифференциальными уравнениями (те. связанные через линейные части и нелинейные звенья. К таким системам относятся, например, система на риса, если в ее уравнениях под знаками нелинейных функций находятся входные (или выходные) величины обоих нелинейных звеньев, и многие другие системы. Системы с логическими устройствами относятся обычно к нелинейным системам второго класса. Заметим, что во всех случаях, когда под знак нелинейной функции входит какая-либо линейная комбинация разных переменных, их следует обозначать одной буквой, а данную линейную комбинацию учесть при составлении общего уравнения линейной части системы. Это бывает, например, в тех случаях, когда на вход нелинейного звена подаются производные или включается обратная связь. Так, если для рис. 16.1, б то, обозначая можно привести уравнение нелинейного звена к виду (16.1). Из всех уравнений линейных звеньев, а также добавочных линейных выражений типа (16.6), получаемых при выделении нелинейности, составляется общее уравнение линейной части системы (16.7) где Q (р) и р) — операторные многочлены) или передаточная функция линейной части системы (16.8) Составление уравнений будет проиллюстрировано ниже на примерах. Процессы в нелинейных системах автоматического регулирования имеют целый ряд весьма существенных особенностей, которые не встречаются в линейных системах. Благодаря этим существенным особенностям даже вопрос об устойчивости системы становится здесь более сложным. Кроме структуры системы и значений ее параметров для устойчивости того или иного установившегося процесса имеют значение здесь, в отличие от линейных систем, также и начальные условия. Возможен новый вид установившегося процесса _ автоколебания, те. устойчивые собственные колебания с постоянной амплитудой при отсутствии внешних колебательных воздействий. Когда в системе возникают автоколебания, то установившееся состояние, соответствующее постоянному значению регулируемой величины, часто становится невозможным. Следовательно, в общем случае на плоскости параметров системы могут быть не два вида областей (устойчивости и неустойчивости, как в линейных системах, а больше 1) область устойчивости равновесного состояния с постоянным значением регулируемой величины 2) область устойчивых автоколебаний 3) область неустойчивости системы 4) области, соответствующие другим, более сложным случаям. Если процессы в системе имеют вид, указанный на риса, то равновесное состояние х = 0) неустойчиво. В том случае, когда оба указанных на риса колебания в переходных процессах стремятся к одной и той же амплитуде и к одной и той же частоте, система будет обладать устойчивыми автоколебаниями с амплитудой а. На рис. 16.3, бив показаны случаи, когда равновесное состояние (х = 0) системы устойчиво в малом, те. при начальных условиях, не выводящих отклонения в переходном процессе за определенную величину аи неустойчиво в большом, те. при начальных условиях, выводящих отклонение в переходном процессе за пределы величины а. Здесь граничным процессом является неустойчивый периодический процесс собственного движения системы с амплитудой а (переходные процессы расходятся от него в обе стороны. На рис. 16.3, г показан случай трех возможных установившихся состояний 1) равновесное состояние (х = 0), 2) колебания с постоянной амплитудой а 3) колебания с постоянной амплитудой а. При этом колебания с амплитудой а неустойчивы. В результате система будет устойчива в малом по отношению к равновесному состоянию ха в большом система будет обладать устойчивыми автоколебаниями с амплитудой а 2 Пример. Для иллюстрации особенностей нелинейной системы исследуем переходный процесс и автоколебания в релейной системе автоматического регулирования температуры, изображенной на рис. 1.35. Для этого составим сначала уравнение регулируемого объекта и регулятора. Пусть регулируемый объект представляет собой некоторую камеру. Учитывая инерционность процесса нагрева и охлаждения, запишем уравнение регулируемого объекта в виде где θ — отклонение температуры, ф — отклонение регулирующего органа, f(t) — внешние возмущения. При отклонении температуры θ появляется ток в диагонали моста того или иного направления (рис. 1.35) и замыкается тот или иной контакт реле 3, включающего постоянное напряжение в ту или иную обмотку возбуждения 4 электродвигателя 5. Приняв во внимание некоторое отставание в этом процессе включения, получим релейную характеристику вида г) рис. 1.36. Далее, считая, что ток I пропорционален отклонению температуры объекта θ , а скорость dt d ϕ отклонения регулирующего органа 6 пропорциональна напряжению на обмотках возбуждения электродвигателя, можно в данном случав выходной величиной для указанной релейной характеристики считать прямо dt d ϕ , а входной — θ (риса. Следовательно, уравнение регулятора запишется здесь следующим образом (16.10) (16.11) Рассмотрим два произвольных участка переходного процесса (при f(t) = 0) в данной системе (участки АВ и В на рис. 16.4, б. На участке АВ уравнение регулятора согласно рис. 16.4, в будет с dt d + = ϕ Дифференцируя (16.9) пои подставляя туда +с, получаем при f(t) =0 следующее уравнение системы регулирования на участке АВ: (16.12) а на участке В (16.13) Решение уравнения (16.12) будет (16.14) откуда получаем (16.15) Условимся для простоты отсчитывать время t от начала участка АВ (риса. Тогда начальные условия будут где А пока неизвестно. Используя начальные условия, находим произвольные постоянные для уравнения (16.15): (16.16) Аналогично для участка В согласно (16.13), отсчитывая время I тоже от начала этого участка (рис. 16.5, б, получим решение |