Математические основы теории автоматического управления

5. - время достижения первого максимума.

6. - время нарастания переходного процесса, время от начала переходного процесса до момента первого пересечения графиком линии установившегося значения.

7. - декремент затухания, равный отношению модулей двух смежных перегулирований -

.

Перечисленные выше показатели могут быть дополнены и другими, если этого требуют специфические технические задания на разработку или исследование систем управления.

Устойчивые переходные процессы, возникающие при ступенчатом воздействии, принято делить на три группы:

1. Монотонные процессы. Такие процессы, где первая производная выходной величины по времени не меняет знак.

2. Апериодические процессы. Здесь производная меняет знак не более одного раза.

3. Колебательные процессы. Производная меняет свой знак периодически

На рис. 6 показан примерный вид колебательного, апериодического и монотонного процессов.

Рис. 6

Требования к качеству переходного процесса могут быть представлены графически. Они сводятся к требованию, чтобы отклонение регулируемой величины при ступенчатом воздействии не выходило за границы некоторой области, изображенной на диаграмме качества процесса. Простейшая диаграмма показана на рис. 7.

Рис. 7

Уточненная диаграмма, в которой накладываются ограничения на скорость нарастания сигнала, показана на рис. 8. На уточненной диаграмме обозначает время запаздывания, равное отрезку времени, заключенному между моментами приложения скачкообразного сигнала () и моментом времени, при котором выходная величина достигает половины установившегося значения.

Рис. 8

3. Контрольные вопросы и задачи

1. В каких системах автоматического управления важен характер и параметры протекания переходных процессов?

2. Перечислите основные области применения критериев качества.

3. Какие показатели качества относят к прямым показателям?

4. Какие показатели (критерии) качества относят к косвенным критериям?

5. Как определить время регулирования по графику переходного процесса?

6. Как определить перегулирование по графику переходного процесса?

7. Как определить время нарастания по графику переходного процесса?

8. Какие переходные процессы относят к апериодическим процессам?

9. Какие переходные процессы относят к монотонным процессам?

Лекция 17

1. Оценка качества при гармонических воздействиях

Качество переходных процессов в системах и объектах при гармонических воздействиях оценивают по частотным характеристикам, снятым экспериментально или рассчитанным по параметрам математической модели. В этом случае получают частотные критерии качества переходных процессов, которые относят к косвенным критериям.

Частотные критерии определяют по одной из частотных характеристик: амплитудно-фазовой, амплитудной, вещественной или логарифмической амплитудной. Рассмотрим в качестве основы определение косвенных показателей качества по амплитудной частотной характеристике (АЧХ) систем.

1. - показатель колебательности, определяемый как отношение максимального значения АЧХ к ее значению при

Рис. 1

Показатель колебательности характеризует склонность систем или объектов к колебательности. Чем выше показатель колебательности, тем более колебательна система, то есть менее качественна.

При переходная характеристика системы неколебательная, а АЧХ системы имеет примерный вид, показанный на рис. 2.

Рис. 2

Если , то говорят о незатухающих колебаниях переходной характеристики, а АЧХ системы имеет вид, показанный на рис. 3.

Рис. 3

Считается допустимым, если .

2. - резонансная частота системы или объекта, при которой АЧХ имеет максимум.

При гармонические сигналы проходят через систему с наибольшим усилением.

3. Полоса пропускания системы управления - это интервал частот от до частоты , в котором выполняется условие

.

Рис. 4

Полоса пропускания систем не должна быть очень широкой, чтобы не проходили высокочастотные помехи. С другой стороны, чем выше частота сигналов, которые пропускает система, тем выше ее быстродействие. В электроприводах полосу пропускания задают в герцах, для систем управления скоростью вращения двигателя полоса пропускания, в зависимости от типа двигателя и преобразователя энергии, находится в диапазоне Гц,

,

при этом накладывается дополнительное ограничение о том, что фазочастотная характеристика в полосе пропускания не должна опускаться ниже .

4. - частота среза, при которой АЧХ системы принимает значение, равное единице

,

как это показано на рис. 5.

Рис. 5

косвенно характеризует длительность переходного процесса, так же как и . Время регулирования системы обратно пропорционально частоте среза -

,

если переходный процесс имеет колебания, то можно установить связь между частотой среза т временем достижения первого максимума -

.

Полоса среза, как характеристика быстродействия систем получила распространение, так как она легко определяется на логарифмических частотных характеристиках (см. рис. 6).

Рис. 6

2. Интегральные оценки качества

Интегральными оценками качества переходного процесса систем управления называют интегралы по времени от некоторых функций переходного процесса изменения ошибки регулирования.

Рассмотрим скалярную линейную систему, показанную на рис. 7. На вход управления системы поступает ступенчатый сигнал с амплитудой , система предназначена для стабилизации заданного значения регулируемой переменной , на систему так же может действовать скалярное ступенчатое возмущение стремящееся снизить значение регулируемой переменной.

Рис. 7

Примерный вид графиков переходных процессов регулируемой величины и ошибки регулирования -

,

показаны соответственно на рис. 8 и 9.

Рис. 8

Рис. 9

Основные области применения интегральных оценок в теории автоматического управления:

1. Общая оценка быстроты затухания и величины отклонения регулируемой величины в совокупности, без определения того и другого в отдельности.

2. Выбор при синтезе параметров систем, обеспечивающих оптимальность переходного процесса с точки зрения достижения минимума интегральных оценок.

Простейшей интегральной оценкой может служить линейная интегральная оценка следующего вида -

.

Геометрическая интерпретация этого интеграла представляет собой площадь под кривой , как это показано на рис. 10 для переходных процессов изменения ошибки при управлении и возмущении.

Рис. 10

Если система управления устойчива и обладает свойством астатизма, тогда

,

а интеграл стремится к конечному значению, равному площади под кривой . Параметры системы управления стремятся выбирать таким образом, чтобы добиться минимума , при этом идеальный переходный процесс будет стремиться к идеальной ступенчатой форма.

3. Контрольные вопросы и задачи

1. Как определить показатель колебательности по АЧХ системы?

2. Как определить частоту среза и полосу пропускания по АЧХ?

3. Как соотносятся частота среза системы и ее быстродействие?

4. Дайте определение линейной интегральной оценке, укажаите ее достоинства и недостатки.

5. Почему для интегральных оценок переходного процесса используют график изменения ошибки регулирования?

6. По графику АЧХ системы

определить показатель колебательности и частоту среза системы. Ответ: Показатель колебательности , частота среза .

7. По графикам АЧХ двух систем: САУ1 - , САУ2 - ,



определить систему управления, переходные процессы которой имеют большую колебательность. Ответ: Большую колебательность переходных процессов имеет система САУ2.

8. По графикам АЧХ двух систем: САУ1 - , САУ2 - ,

определить систему управления, которая имеет большее быстродействие. Ответ: Большую колебательность переходных процессов имеет система САУ1.

Лекция 18

1. Вычисление линейных интегральных оценок

Рассмотрим проблему вычисления интеграла линейной интегральной оценки. Можно сначала решить аналитически дифференциальные уравнения, описывающие систему, долее определить ошибку регулирования, затем подставить выражение для ошибки в интеграл линейной оценки и, взяв его, получить выражение для .

Но можно поступить и иначе.

Пусть свободное движение ошибки регулирования системы описывается уравнением

(1)

Проинтегрируем это уравнение -

После интегрирования получаем -

(2)

Подстановки верхнего предела дают члены следующего вида -

 (3)

так как все производные ошибки в установившемся режиме обращаются в ноль.

Подстановки нижнего предела дают члены вида -

(4)

которые являются начальными условиями уравнения (1).

Подставив (3) и (4) в (2), получим

(5)

А так как

,

окончательно получаем

(6)

Решая (6) относительно , получим выражение для вычисления линейной интегральной ошибки -

 (7)

Теперь мы может определить по коэффициентам характеристического уравнения системы и начальным условиям переходного процесса ошибки.

Для синтеза систем, определения параметров минимизирующих , следует воспользоваться обычными методами исследования функций на экстремум. Следовательно, если мы хотим определить параметр системы, на пример, параметр , обеспечивающий , необходимо решить относительно параметра следующее уравнение -

.

Рассмотрим несколько примеров использования линейной интегральной оценки.

Пример

Система имеет характеристическое уравнение

(8)

Определим выражение для , если начальные условия имеют вид -

.

Определим значение параметра , при котором интегральная оценка имеет минимум.

Решение

Обозначим -

.

Используем для нахождения выражение (7) -

(9)

Из рассмотрения (9) получаем, что в этом случае не имеет экстремума, а меньшее значение интегральной ошибки мы будем получать при меньшем значении . Действительно, ведь уравнение (8) является характеристическим уравнением апериодического звена, параметр - это постоянная времени. Переходный процесс для двух разных постоянных времени будет иметь вид, показанный на рис. 1.

Рис. 1

Пример

Система имеет характеристическое уравнение

.

Определим выражение для , если начальные условия имеют вид -

.

Определим значение параметра , при котором интегральная оценка имеет минимум.

Решение

Обозначим -

.

Используем для нахождения выражение (7) -

.

Если , то процессы монотонные, обеспечивается при наименьших и . Если , то уменьшение коэффициента затухания уменьшает линейную интегральную оценку, но это приводит к ухудшению переходного процесса, повышению его колебательности.

При колебательных процессах в системах линейная интегральная оценка дает значительную погрешность. При этом минимум оценки может соответствовать процессу с большим числом колебаний со значительной амплитудой, малым быстродействием, так как, по сути, в оценке происходит сложение положительных и отрицательных областей площади под интегральной кривой. Это иллюстрируют рис. 2 и 3, показывая два процесса, которые могут иметь одно и то же значение линейной интегральной оценки.



Рис. 2

Рис. 3

И так как форма переходного процесса при анализе системы автоматического управления часто заранее неизвестна, то применять линейные интегральные оценки на практике нецелесообразно.

Можно попытаться использовать интеграл от модуля ошибки следующего вида -

(10)

На рис. 4 показан примерный вид кривых изменения ошибки и ее модуля. Но аналитическое вычисление интеграла от модуля ошибки по математической модели системы оказалось весьма громоздким, поэтому эта оценка широкого распространения не получила.

Рис. 4

2. Квадратичная интегральная оценка

В большинстве случаев, при возможности возникновения в системе колебательного переходного процесса, используют квадратичную интегральную оценку, которая имеет следующий вид -

(11)

Оценка не зависит от знака отклонений ошибки, а значит и от формы переходного процесса, монотонный, апериодический или колебательный характер он будет иметь. На рис. 5 и 6 показан примерный вид кривых изменения ошибки и квадрата ошибки.

Рис. 5

Рис. 6

Рассмотрим процедуру вычисления квадратичной оценки по математической модели системы. Система управления представляется в виде, показанном на рис. 7.



Рис. 7

Изображение по Лапласу сигнала на выходе системы имеет вид -

(12)

где - изображение по Лапласу единичной ступенчатой функции - входного сигнала системы.

Для системы автоматического управления, математическая модель которой приведена к виду (12), интегральная квадратичная ошибка определяется по следующему выражению -

(13)

где (14)

в все элементы с индексами меньше 0 и больше заменяются 0.

Определители в (13), где , получаются заменой в определителе (14) ()-го столбца столбцом следующего вида -

.

Коэффициенты в выражении (13) определяются следующим образом -

(15)

при определении коэффициенты, индексы которых меньше 0 и больше , заменяются 0.

3. Контрольные вопросы и задачи

1. Какие параметры математической модели объекта требуются для вычисления линейной интегральной оценки?

2. Почему нельзя использовать линейную интегральную оценку в случае колебательного характера переходных процессов?

3. Какие интегральные оценки целесообразно использовать в том случае если в системе возможно наличие колебательных переходных процессов?

4. Дайте определение квадратичной интегральной оценке переходного процесса.

5. При минимизации квадратичной оценки, к какому виду стремится переходный процесс?

6. Какие параметры математической модели объекта требуются для вычисления квадратичной интегральной оценки?

7. Объект управления описывается передаточной функцией - . Вычислите линейную интегральную оценку переходного процесса при начальном значении ошибки . Ответ: Линейная интегральная оценка .

8. Объект управления описывается передаточной функцией - . Вычислите линейную интегральную оценку переходного процесса при начальном значении ошибки . Ответ: Линейная интегральная оценка .

Лекция 19

1. Квадратичная интегральная оценка с учетом производной

Недостатком квадратичной интегральной оценки , как и предыдущих оценок, является то, что при минимизации оценки не накладываются ограничения на форму переходного процесса. На пример, показанные на рис. 1 графики - (а, б, в) могут иметь одинаковые значения существенно при этом отличаясь по форме переходного процесса.

Рис. 1

Кроме того, часто оказывается, что выбранные по параметры системы приводят к существенно колебательному процессу, большим производным из-за стремления приблизить процесс к идеальному скачку.

Поэтому используют еще один вид интегрально квадратичной оценки, в которой ограничение накладывается не только на величину отклонения , но и на скорость его изменения . Эта оценка имеет следующий вид -

 (1)

где - некоторая постоянная времени.

Разницу между оценками и можно представить графически, как это показано на рис. 2.

Рис. 2

То есть оптимизированный по переходный процесс стремиться к идеальному скачку, а оптимизированный по - к кривой экспоненциального вида, которая описывается следующим выражением -

.

Докажем последнее утверждение. Для этого проанализируем выражение (1).

,

с учетом того, что

,

Получаем

(2)

С учетом того, что последнее слагаемое в (2) является величиной постоянной -

,

квадратичная оценка будет иметь минимум при

(3)

Решение дифференциального уравнения (3) имеет вид -

,

а если перейти от ошибок к выходным переменным, то получим -

,

что и требовалось доказать.

Следовательно, выбирая параметры системы по , можно приблизить переходный процесс к экспоненте с заданной постоянной времени , тем самым вводится ограничение на скорость нарастания выходной величины .

Методика определения может быть аналогичной методике определения , рассмотренной выше, если представить квадратичную оценку с учетом производной в следующем виде -

,

где определяется по формулам для , но с учетом того, что порядок числителя - увеличивается на 1.

В теории автоматического управления используют квадратичные оценки с производными более высокого порядка (до ) для более точного задания желаемой формы переходного процесса, естественно, что при этом усложняется и процесс вычисления оценок.

2. Вычисление квадратичных интегральных оценок

Рассмотрим вычисление и использование квадратичных ошибок на примере.

Пример

В системе управления с передаточной функцией -

,

зададим :

· из условия ,

· из условия ,

и сравним переходные процессы для двух этих случаев.

Решение

Получим выражение для . Для этого преобразуем передаточную функцию системы к заданному виду

,

тогда получим

(4)

Выражение для принимает вид -

(5)

Определим компоненты (5) по параметра передаточной функции системы (4).

(6)

Для нахождения определим (), при ,

,

Заменим в выражении (6) для первый столбец столбцом вида

.

Тогда получаем

.

Определим -

.

После подстановки полученных компонент в (5) получаем выражение для квадратичной интегральной оценки.

(5)

Найдем выражение для частной производной по от выражения (5)

,

приравнивая полученное выражение к нулю получаем уравнение для нахождения оптимального значения .

.

В результате получаем оптимизированное по квадратичной оценке значение -

(6)

Передаточная функция системы при примет вид -

.

На рис. 3 покажем вид переходного процесса системы при единичном ступенчатом воздействии и оптимизированным по параметром.

Рис. 3

Таким образом, имеем следующие показатели качества переходного процесса,

(7)

Определим по отработанной выше методике для -

,

выражение для берем из предыдущего случая -

.

Определим теперь . Передаточная функция системы для этого случая имеет вид -

,

тогда получим

(8)

Выражение для принимает вид -

(9)

Определим компоненты (9) по параметра передаточной функции системы (8).

(10)

Определим коэффициенты -

.

не определяем, так как . Для нахождения определим (), при ,

,

Заменим в выражении (10) для второй столбец столбцом вида

.

Тогда получаем

.

После подстановки полученных компонент в (9) получаем выражение для квадратичной интегральной оценки.

(11)

Окончательно получаем

(12)

Найдем выражение для частной производной по от выражения (12)

,

приравнивая полученное выражение к нулю получаем уравнение для нахождения оптимального значения .

.

В результате получаем оптимизированное по квадратичной оценке с учетом производной значение -

(13)

Полагаем для определенности , тогда

.

Передаточная функция системы при примет вид -

.

На рис. 3 покажем вид переходного процесса системы при единичном ступенчатом воздействии и оптимизированным по параметром.

Рис. 4

Таким образом, имеем следующие показатели качества переходного процесса,

(14)

Сравнивая переходные процессы, видим, что при оптимизации по квадратичной оценке с учетом производной () получили существенно меньшие значения перерегулирования и быстродействия, при более плавном нарастании переменной.

3. Контрольные вопросы и задачи

1. Дайте определение квадратичной интегральной оценки с учетом производной, поясните ее компоненты.

2. К какому виду стремиться переходный процесс при минимизации интегральной квадратичной оценки с учетом производной?

3. Как вычисляют квадратичную интегральную оценку с учетом производной?

4. Вычислите интегральную квадратичную оценку переходного процесса в системе с передаточной функцией - , если на вход системе подается единичная ступенчатая функция. Ответ: Интегральная квадратичная оценка .

5. Вычислите интегральную квадратичную с учетом производной оценку переходного процесса в системе с передаточной функцией - , если на вход системе подается единичная ступенчатая функция, а постоянная времени оценки . Ответ: Интегральная квадратичная оценка .

6. Определите параметр регулятора системы управления, обеспечивающий минимум квадратичной оценки



Ответ: Параметр пропорционально-интегрального регулятора .

Лекция 20

1. Корневые критерии качества переходных процессов

Эта группа критериев основана на оценке качества переходных процессов по значениям полюсов и нулей передаточной функции системы между интересующими нас входами и выходами системы.

Как известна, переходная характеристика системы может быть определена следующим образом -

(1)

где - корни характеристического уравнения системы

.

Очевидно, что на характер переходного процесса оказывает влияние и числитель и знаменатель передаточной функции. Но, в большинстве случаев, при анализе систем по реакции на управляющее воздействие, не имеет корней, то есть передаточная функция не имеет нулей. Тогда характер переходного процесса можно оценить только по полюсам передаточной функции, подвергая тем самым анализу корни характеристического уравнения системы -

(2)

В случае приближенной оценки качества по корням характеристического уравнения на комплексной плоскости выделяют область расположения корней, границы которой задаются по требованиям к качеству процессов, как это показано на рис. 1.

Рис. 1

Границы области, показанной на рис. 1, задаются следующими параметрами:

· - критерий длительности переходного процесса,

· - колебательность переходного процесса, определяется по ,

· - максимальное удаление корня от мнимой оси.

Рассмотрим эти параметры.

Критерий длительности определяется как расстояние от мнимой оси до ближайшего действительного корня или ближайшей пары комплексно сопряженных корней.

Выясним, действительно ли этот параметр характеризует длительность переходного процесса? Возможны два случая расположения корней на границе области.

1. Пусть ближайшим к мнимой оси, то есть лежащий на границе области, будет действительный корень -

,

тогда соответствующая ему компонента переходного процесса, в соответствии с (1) будет иметь вид -

(3)

где - коэффициент разложения (1).

2. Если ближайшей к мнимой оси будет комплексно-сопряженная пара корней -

,

тогда соответствующая им компонента переходного процесса, в соответствии с (1) будет иметь вид -

(4)

где - частота колебаний.

Из (3) и (4) мы видим, что время затухание компоненты определяет сомножитель -

,

где - величина минимального действительного корня или минимальной действительной части корней, - соответствующая , наибольшая постоянная времени. Таким образом, можно считать, что переходный процесс системы завершится не раньше, чем затухнет компонента . Следовательно, определяет длительность переходного процесса, будучи величиной, обратно пропорциональной времени регулирования. Зная , мы можем оценить время регулирования или переходного процесса по следующему соотношению -

,

где - половина ширины области, при попадании в которую переходной процесс считается завершенным. Если , а крайний корень действительный, то имеем -

.

Критерий колебательности определяется по углу следующим образом -

.

где - соответственно действительная и мнимая части комплексно сопряженной пары корней расположенных на границе области (см. рис. 1). При увеличении возрастает колебательность системы.

Дальнюю от мнимой оси границу области , определяют корни, оказывающие предельно малое влияние на переходный процесс.

При прочих равных условиях от системы требую увеличения и снижения .

В качестве примера влияния расположения корней на характер переходных процессов покажем графики, представленные на рис. 2 и 3.

Рис. 2

Рис. 3

Если передаточная функция системы имеет нули, то оценка качества системы только по полюсам может дать существенную погрешность.

Чтобы пояснить характер влияния нулей на качество переходных процессов, представим систему следующим образом, как это показано на рис. 4.



Рис. 4

Конкретизируем задачу, пусть

,

а имеет вид, показанный на рис. 5. При этом рассмотрим два варианта графика:

· ,

· .

Из рассмотрения рис. 5 можно сделать вывод, что члены с положительными коэффициентами приводят к повышению колебательности и быстродействия, а отрицательные коэффициенты затягивают переходный процесс.

Рис. 5

В тех случаях, когда требуется получить желаемый вид переходного процесса, используют методы, основанные на связи коэффициентов характеристического уравнения системы или его корней с видом переходного процесса, с заданными динамическими показателями.

Рассмотрим характеристическое уравнение вида -

(5)

Преобразуем (5)

(6)

По формулам Виета определяется как сумма всех корней уравнения, - сумма произведений всех пар корней, - сумма произведений всех троек корней и т. д., а определяется как произведение всех корней уравнения -

.

Теперь, если мы сможем задать расположение корней на комплексной плоскости, исходя из требований качества динамики, то по формулам Виета можно найти значения коэффициентов характеристического уравнения, которые связаны с параметрами системы.

Обратим особое внимание на коэффициент , чем больше , то, при прочих равных условиях, больше действительные части корней, следовательно, быстрее переходный процесс. Если корни действительные и кратные, тогда -

.

Обозначим

(7)

где носит название среднегеометрического корня характеристического уравнения.

Тогда уравнение (6) с учетом (7) имеет вид -

На комплексной площади расположения корней характеристического уравнения определяет точку на действительной оси - геометрический центр всех корней системы, а коэффициенты определяют взаимное расположение корней. При этом легко показать, что определяют кривую переходного процесса в относительном времени , а величина определяет масштаб времени для этого процесса.

На практике рассмотренный выше подход используют следующим образом:

1. Для конкретной системы определяют требуемый вид переходного процесса.

2. Для обеспечения заданных требований выбирают из имеющихся в справочной литературе предварительно рассчитанные значения коэффициентов характеристического уравнения, тем самым выбирается "желаемое" характеристическое уравнение -

 (8)

3. Определяют характеристическое уравнение по структуре и параметрам системы -

(9)

4. где - коэффициенты, функционально связанные с параметрами системы.

5. Получают систему алгебраических уравнений, приравняв коэффициенты уравнений (8) и (9) при одинаковых степенях оператора Лапласа -

(10)

6. Решают систему (10) относительно изменяемых параметров системы (параметров регуляторов), что позволяет определить параметры, обеспечивающие заданный вид и качество переходного процесса.

Описанный выше алгоритм часто называют методом стандартных коэффициентов или стандартного расположения корней характеристического уравнения системы управления. Рассмотрим в качестве иллюстрации два стандартных расположения корней, которые наиболее распространенны в системах управления электромеханическими приводами различных установок.

2. Биномиальное распределение корней

Биномиальное распределение корней используют для обеспечения заданного быстродействия при монотонности переходных процессов. Стандартное биномиальное характеристическое уравнение имеет вид -

В этом случае имеем кратных действительных корней с отрицательной действительной частью, равной . Вид переходных процессов для от 1 до 4 показан на рис. 6. Характеристические уравнения для этих случаев имеют вид -

Распределение Баттерворта

Корректным является сопоставление системы автоматического управления и идеальным фильтром низкой частоты (ФНЧ), когда для полосы пропускания системы (НЧ) требуют максимальной горизонтальности ЛАЧХ, что обеспечивает пропускание без искажений сигналов управления. Для диапазона высоких частот (ВЧ) требуют максимального подавления сигнала, так как это диапазон сигналов помех. Рис. 7 иллюстрирует приближение желаемой характеристики системы к характеристике "идеального" фильтра низкой частоты.

Распределение корней по Баттерворту обеспечивает компромисс между этими требованиями, достигая высокой равномерности в полосе пропускания НЧ при приемлемой крутизне характеристики в полосе подавления ВЧ.

Рис. 6

Рис. 7

При этом корни характеристического уравнения располагаются на комплексной плоскости, на окружности с радиусом и угловым расстоянием между корнями - , симметрично относительно действительной оси, как это показано на рис. 8.

Рис. 8

Вид переходных процессов для от 1 до 4 показан на рис. 9.

Рис. 9

Характеристические уравнения и параметры переходного процесса для этих случаев имеют вид -



Сравнение переходных характеристик показывает, что распределение Баттерворта обеспечивает более высокое, чем биномиальное распределение, быстродействие с малым перерегулированием и колебательностью.

3. Контрольные вопросы и задачи

1. Как объяснить влияние на переходные процессы корней характеристического уравнения?

2. Какую компоненту переходного процесса дает отрицательный действительный корень характеристического уравнения?

3. Какие компоненты переходного процесса дают комплексно сопряженные корни характеристического уравнения?

4. Что определяют корни характеристического уравнения ближе всего расположенные к мнимой оси комплексной плоскости?

5. Как связана с быстродействием системы величина среднегеометрического корня характеристического уравнения?

6. Какое влияние оказывает на переходный процесс нули передаточной функции?

7. В каких случаях следует использовать на настройки системы биномиальное распределение корней характеристического уравнения?

8. В каких случаях следует использовать на настройки системы распределение корней характеристического уравнения Баттерворта?

9. Определите коэффициенты характеристического уравнения с биномиальным распределением корней для системы управления третьего порядка, если требуемое время регулирования . Ответ: Желаемое характеристическое уравнение имеет вид - .

10. Определите коэффициенты характеристического уравнения с распределением корней по Баттерворту для системы управления четвертого порядка, если требуемое время регулирования . Ответ: Желаемое характеристическое уравнение имеет вид - .

Лекция 21

1. Определения и задачи идентификации математических моделей

Идентификация динамических объектов в общем случае состоит в определении их структуры и параметров по наблюдаемым данным - входному воздействию и выходным величинам.

В этом случае объект (элемент системы, объект управления, элемент технологического процесса и т. п.) представляет собой "черный ящик". Исследователю необходимо, подвергая объект внешним воздействиям и анализируя его реакции, получить математическую модель (описание его структуры и параметров), то есть превратить "черный ящик" в "белый ящик", добиться его "информационной прозрачности". Графически процесс идентификации иллюстрирует рис. 1.

Рис. 1

Важным моментом этого процесса является выбор точек приложения внешних воздействий и сбор информации о реакциях объекта, то есть размещение управляющих устройств и датчиковых систем.

Решается при идентификации объектов и более простая (относительно простая) задача, это задача идентификации параметров, когда заранее известна структура математической модели объекта, но не известны ее параметры. В этом случае говорят о переходе от "серого ящика" к "белому ящику". Графически процесс идентификации параметров иллюстрирует рис. 2.

Рис. 2

Задача идентификации параметров может либо входить компонентом в общую задачу идентификации объекта, либо решаться самостоятельно.

Рассмотрим на обобщенной структуре и процедуре процесса идентификации основы подхода к решению задач идентификации. Обобщенная структура процесса идентификации показана на рис. 3.

Обобщенная процедура идентификации

1. Классификация объекта.

2. Выбор для определенного класса объекта настраиваемую модель, то есть модель, структуру и параметры которой можно менять в процессе идентификации.

3. Выбрать критерий (оценку) качества идентификации, характеризующий в виде функционала доступных для наблюдения переменных отличие модели и объекта.

4. Выбрать алгоритм идентификации (механизм настройки модели), обеспечивающий сходимость процесса идентификации, минимум критерия качества идентификации.



Рис. 3

Методы идентификации принято разделять на две группы:

· активная идентификация - идентификация вне контура управления,

· пассивная идентификация - идентификация в контуре управления.

Активная идентификация

В этом случае объект исследования выводится из условий нормальной окружающей среды (нормальный режим эксплуатации, номинальные параметры рабочего режима и т. п.). Исследования проводятся в специализированных лабораторных условиях, как это показано на рис. 4. На входы объекта (рабочие и дополнительные) подаются тестовые сигналы специального вида. Это могут быть:

· ступенчатые и импульсные временные сигналы,

· гармонические сигналы,

· случайные воздействия с заданными параметрами.

Активную идентификацию используют при разработке новых технологий применительно к действующим промышленным объектам, в изучении новых явлений, в первоначальной разработке математической модели.

Пассивная идентификация

При пассивной идентификации объект функционирует в контуре управления, находится в процессе нормальной эксплуатации. На его входы поступают только естественные сигналы управления.

Пассивную идентификацию используют для уточнения математической модели, для слежения за изменениями в объекте. Информация оперативно используется в системе управления объектом, процесс такой идентификации иллюстрируется рис. 5.

Рис. 4

Рис. 5

Кроме перечисленных групп методов реализуются и системы идентификации смешанного типа, когда объект не выводится из нормального режима эксплуатации, но к управляющим сигналам добавляются тестовые воздействия, позволяющие идентифицировать объект, не ухудшая качества основного процесса управления.

Более подробно рассмотрим активную идентификацию.

Активная идентификация объектов управления может производиться как во временной области, так и в частотной области. При этом В каждой области используют собственные алгоритмы и методы идентификации.

При активной идентификации в большинстве случаев используют полученные в результате экспериментов характеристики:

· частотные характеристики (АФЧХ, ЛАЧХ, ЛФЧХ и др.),

· временные характеристики (ступенчатое изменение задания, "узкий" импульс задания и др.).

Рассмотрим в качестве примера один из подходов решения задачи идентификации структуры и параметров объекта в частотной области. Ограничим рассмотрение объектом с одним входом и одним выходом.

Мы знаем, что если имеется математическая модель такого объекта в виде передаточной функции -

(1)

то это соответствует наличию полной информации о структуре и параметрах объекта, всех его характеристиках.

Преобразуем передаточную функцию (1) к полюсно-нулевому представлению, форме Боде -

(2)

где - коэффициент усиления (), - соответственно полюсы и нули передаточной функции.

Если среди корней () встречаются комплексно сопряженные пары корней, то разложение (2) необходимо дополнить сомножителями следующего типа -

.

Предполагая для простоты изложения отсутствие комплексно сопряженных корней, можно преобразовать (2) к следующему виду -

(3)

где .

По выражению для передаточной функции в виде (3) получим частотную характеристику объекта -

,

ЛАЧХ -

(4)

ЛФЧХ -

 (5)

С другой стороны, нам известно, что ЛАЧХ и ЛФЧХ динамических звеньев с передаточными функциями -

имеет вид, показанный на рис. 6, так как звенья являются соответственно форсирующим и апериодическим динамическими звеньями первого порядка.

Рис. 6

Исходя из изложенного материала, можно предложить следующую процедуру активной идентификации структуры и параметров линейной системы с одним входом и одним выходом:

1. В процессе эксперимента с объекта снимается частотная характеристика в виде ЛАЧХ и ЛФЧХ.

2. Полученная экспериментально ЛАЧХ аппроксимируется кусочно-линейной криво1 - набором отрезков (асимптот) с целочисленным наклоном кратным .

3. По наклону асимптот и частотам сопряжения асимптот определяется передаточная функция объекта в виде произведения передаточных функций соответствующих асимптотам элементарных динамических звеньев (апериодических и форсирующих).

При наличии в полученной ЛАЧХ и ЛФЧХ признаков звеньев второго порядка, то есть асимптот с наклоном кратным , необходимо ввести такие звенья в модель.

Колебательное звено с передаточной функцией -

,

имеет ЛАЧХ и ЛФЧХ, показанные на рис. 7.

Рис. 7

Форсирующее звено второго порядка с передаточной функцией

,

имеет ЛАЧХ (ЛФЧХ) симметричные показанным на рис. 7 характеристикам колебательного звена относительно оси частот.

Рассмотрим пример идентификации по рассмотренной процедуре.

Пример

По экспериментально полученной ЛАЧХ объекта определить передаточную функцию.

Решение

Аппроксимируем экспериментальную ЛАЧХ набором асимптот, как это показано на рис. 8.

Рассмотрим теперь участки аппроксимированной ЛАЧХ, на которых наклон не изменяется:

1. .

На этом интервале виду асимптоты соответствует ЛАЧХ интегрирующего звена, его передаточная функция -

.

Этому звену соответствует следующее выражение ЛАЧХ -

.

Используем последнее выражение для определение , подставив значение характеристики при частоте

.



Рис. 8

2. .

На этом интервале наклон асимптоты возрос на 20 дБ/дек, что соответствует добавлению апериодического звена первого порядка с передаточной функцией -

,

где постоянная времени определяется по точке сопряжения асимптот -

.

3. .

На этом интервале наклон асимптоты уменьшился на 20 дБ/дек, что соответствует добавлению форсирующего звена первого порядка с передаточной функцией -

,

где постоянная времени определяется по точке сопряжения асимптот -

.

4. .

На этом интервале наклон асимптоты возрос на 20 дБ/дек, что соответствует добавлению апериодического звена первого порядка с передаточной функцией -

,

где постоянная времени определяется по точке сопряжения асимптот -

.

Перемножая полученные передаточные функции, получим передаточную функцию объекта -

.

2. Контрольные вопросы и задачи

1. Дайте определение процессу идентификации объекта, идентификации параметров.

2. Какие методы относят к методам пассивной идентификации?

3. Какие методы относят к методам активной идентификации?

4. Поясните процедуру идентификации параметров объекта управления по его экспериментальной ЛАЧХ.

5. Асимптотическая ЛАЧХ объекта имеет вид -

6.

Определите передаточную функцию объекта. Ответ: .

7. Асимптотическая ЛАЧХ объекта имеет вид -



Определите передаточную функцию объекта. Ответ: .

8. Асимптотическая и точная ЛАЧХ динамического звена имеет вид -

Определите передаточную функцию звена. Ответ: .

Размещено на Allbest.ru