К.Ю. Поляков - Теория автоматического управления для чайников (1054016), страница 7
Текст из файла (страница 7)
Это значит, что все его переменные состояния равнынулю и внутренняя энергия также нулевая.Если начальные условия ненулевые, то для построения сигнала выхода при любом входенужно использовать дифференциальные уравнения объекта или модель в пространстве состояний. Это значит, что переходная характеристика дает меньше информации, чем исходные уравнения.Пусть модель объекта задана дифференциальным уравнением первого порядка:dy (t )T+ y (t ) = k ⋅ x(t ) ,(16)dtгде k – безразмерный коэффициент, а T – некоторая постоянная, которая имеет размерностьвремени (измеряется в секундах). Найдем переходную характеристику этого звена.
Решая уравнение (16) при x(t ) = 1 ( t > 0 ), получаем⎛ t⎞y (t ) = k + C1 ⋅ exp⎜ − ⎟ ,⎝ T⎠где постоянная C1 должна определяться из начальных условий. Поскольку нас интересует переходная характеристика, начальные условия считаем нулевыми, то есть y (0) = 0 , что даетC1 = −k и поэтому⎡⎛ t ⎞⎤h(t ) = y (t ) = k ⎢1 − exp⎜ − ⎟⎥ .(17)⎝ T ⎠⎦⎣На рисунке показаны переходные характеристики (17) при различных значениях параметра T,который называется постоянной времени звена:yT = 0,5 ckT =1c0123456tВидно, что при увеличении T выход y медленнее достигает установившегося значения, равного k , то есть постоянная времени характеризует инерционность звена (16).
Чем больше постоянная времени, чем медленнее реагирует объект на управление и тем больше усилий нужно длятого, чтобы перевести его в новое состояние.Заметим, что ступенчатый сигнал легко получить на практике, поэтому переходную характеристику можно снять экспериментально.23© К.Ю. Поляков, 20083.4.
Импульсная характеристика (весовая функция)В качестве тестового сигнала можно, в принципе, использовать любой сигнал. Например,можно изучать реакцию системы на прямоугольный импульс. Вопрос в том, чтобы определитьнекоторый стандартный вид этого импульса. На рисунках а)-в) показаны три импульса, имеющих одинаковые площади. Для простоты будем считать, что эта площадь равна единице.а)б)в)г)1δ (t )t00ttt00Что будет, если мы будем уменьшать ширину импульса, сохраняя его площадь? Очевидно, что высота импульса будет расти и в пределе (когда ширина стремится к нулю) станет бесконечной. Таким образом, мы получили еще один классический тестовый сигнал – единичныйимпульс или дельта-функцию Дирака δ (t ) .
Это идеальный (невозможный в реальной жизни)сигнал, который равен нулю во всех точках, кроме t = 0 , где он уходит к бесконечность, причем его площадь (интеграл по всей оси времени) равен единице:∞⎧ ∞, t = 0δ (t ) = ⎨,∫ δ (t ) dt = 1.⎩ 0, t ≠ 0−∞Поскольку бесконечный импульс невозможно нарисовать, на графике он изображается стрелкой, высота которой равна единице (см. рисунок г).Иногда определяют дельта-функцию как производную от единичного ступенчатого сигнала 1(t ) . Действительно, эта производная равна нулю при всех значениях t , кроме нуля, где онаобращается в бесконечность.Реакция системы на единичный импульс (дельта-функцию) называется импульсной характеристикой и обозначается w(t):w(t)δ(t)δ(t)w(t)U00ttИмпульсная характеристика, так же, как и переходная характеристика, определяется при нулевых начальных условиях, то есть, объект должен находиться в состоянии покоя.Рассматривая дельта-функцию как предельный случай прямоугольного сигнала единичной площади, можно найти связь между переходной функцией и импульсной характеристикой.Пусть ширина прямоугольного импульса равна ε , а высота – 1 / ε .
Такой импульс можнопредставить в виде разности двух ступенчатых сигналов1x(t ) = [1(t ) − 1(t − ε )] ,εгде 1(t − ε ) – это единичный ступенчатый сигнал, который приходит в момент t = ε , то есть,смещен по времени на ε (см. рисунок далее).24© К.Ю. Поляков, 2008x(t )11ε1(t − ε )1(t )1εt0 εt0Так как для линейных систем справедлив принцип суперпозиции, сигнал на выходе будет равенразности реакций системы на входы 1(t ) и 1(t − ε ) , умноженной на коэффициент 1 / ε .
Учитывая, что реакция на сигнал 1(t ) – это переходная функция h(t ) , получаем1y (t ) = [h(t ) − h(t − ε )] .εПереходя к пределу при ε → 0 , наодим, что импульсная характеристикаh(t ) − h(t − ε ) dh(t )w(t ) = lim=,ε →0εdtкак оказывается, равна производной от переходной функции. Наоборот, переходная функция –это интеграл от импульсной характеристики на интервале от 0 до t:th(t ) = ∫ w(τ ) dτ .0Дифференцируя переходную характеристику (17) звена первого порядка, получаем соответствующую импульсную характеристику:d⎛ ⎡⎛ t⎞⎛ t ⎞⎤ ⎞ kw(t ) = ⎜⎜ k ⎢1 − exp⎜ − ⎟⎥ ⎟⎟ = exp⎜ − ⎟ .dt ⎝ ⎣⎝ T⎠⎝ T ⎠⎦ ⎠ TДругое название импульсной характеристики – весовая функция.
Это название связано стем, что для произвольного входного сигнала x(t ) выход системы y (t ) при нулевых начальныхусловиях вычисляется как интегралy (t ) =t∞−∞0∫ x(τ ) w(t − τ ) dτ = ∫ x(t − τ ) w(τ ) dτ .Здесь функция w(t ) как бы «взвешивает» входной сигнал x(t ) в подынтегральном выражении.Заметим, что импульсная характеристика дает неполную информацию об объекте, поскольку неучитывает ненулевые начальные условия.В отличие от ступенчатого сигнала, мгновенный импульс бесконечной величины невозможно получить на реальном устройстве, поэтому снять импульсную характеристику системы,строго говоря, экспериментально не удается.3.5. Передаточная функцияВы уже знаете, выходной сигнал системы можно представить как результат действия некоторого оператора на ее вход.
Для линейных моделей такой оператор можно записать следующим образом.Пусть модель объекта задана линейным дифференциальным уравнением второго порядка,связывающим вход x(t ) и выход y (t ) :d 2 y (t )dy (t )dx(t )b2+ b1+ b0 y (t ) = a1+ a0 x(t )2dtdtdtгде ai (i = 0,1) и bi (i = 0,1,2) – постоянные.(18)25© К.Ю. Поляков, 2008Введем оператор дифференцирования p =d, который действует на сигнал x(t ) по праdtdx(t ). Обратите внимание, что запись p x(t ) обозначает не умножение оператораdtp на x(t ) , а действие этого оператора, то есть дифференцирование x(t ) .Теперь запишем производные сигналов x(t ) и y (t ) по времени в операторной формевилу p x(t ) =dx(t )dy (t )d 2 y (t )= py (t ), &y&(t ) == p 2 y (t ), x& (t ) == px(t ) .2dtdtdtПодставляя эти выражения в (18), получим(19)b2 p 2 y (t ) + b1 py (t ) + b0 y (t ) = a1 px(t ) + a0 x(t ) .Можно формально вынести за скобки y (t ) в левой части равенства (19) и x(t ) в правой части:(b2 p 2 + b1 p + b0 ) y (t ) = (a1 p + a0 ) x(t ) .(20)y& (t ) =Левая часть (20) означает, что оператор b2 p 2 + b1 p + b0 действует на сигнал y (t ) , а в правой части оператор a1 p + a0 действует на сигнал x(t ) .
«Разделив» (условно, конечно) обе части (20) наоператор b2 p 2 + b1 p + b0 , связь выхода и входа можно записать в видеa1 p + a0y (t ) =x(t ) = W ( p) x(t ) ,(21)b2 p 2 + b1 p + b0где запись W ( p ) x(t ) означает не умножение, а действие сложного оператораa1 p + a0W ( p) =.(22)b2 p 2 + b1 p + b0на сигнал x(t ) . Иначе говоря, формула y (t ) = W ( p) x(t ) – это не что иное, как символическаязапись уравнения (18), которую удобно использовать.Функция W ( p) называется передаточной функцией объекта, который описываетсяуравнением (18). Она полностью описывает связи между выходом и входом объекта при нулевых начальных условиях, но не учитывает его внутреннее устройство.Часто передаточной функцией называют функцию W (λ ) , которая получается из (22) в результате замены оператора p на некоторую независимую переменную λ .
Эта фукнция представляет собой отношение двух полиномов (многочленов) от λ .Передаточная функция W (λ ) называется правильной, если степень ее числителя небольше, чем степень знаменателя; строго правильной, если степень числителя меньше степенизнаменателя; неправильной, если степень числителя больше, чем степень знаменателя.
Напри1λмер, функция– строго правильная и одновременно правильная;– правильная, но неλ +1λ +1λ2 + λ + 1– неправильстрого правильная (иногда такие функции называют биправильными), аλ +1ная.Нулями передаточной функции называются корни ее числителя, а полюсами – корниλ −1знаменателя. Например, функция W (λ ) = 2имеет нуль в точке λ = 1 и два полюса вλ + 3λ + 2точках λ = −1 и λ = −2 .3.6. Преобразование Лапласа3.6.1.
Что такое преобразование Лапласа?Одна из первых задач, которые были поставлены в теории управления – вычисление выхода системы при известном входе. Мы видели, что для ее решения нужно решать дифференци26© К.Ю. Поляков, 2008альные уравнения. Чтобы упростить процедуру, математики придумали преобразование, которое позволило заменить решение дифференциальных уравнений алгебраическими вычислениями, то есть, операциями с полиномами (многочленами) и рациональными функциями.Для функции f (t ) вводится преобразование Лапласа, которое обозначается как L{ f (t )} :∞F ( s ) = L{ f (t )} = ∫ f (t ) e −st dt .(23)0Функция F ( s) называется изображением для функции f (t ) (оригинала).