Всё для Учёбы — студенческий файлообменник
1 монета
pdf

Студенческий документ № 000013 из ДГТУ (бывш. РИСХМ)

Кудряшев Сергей Борисович

Канд. техн. наук доцент каф. "Автоматизация

производственных процессов"

Введение

• Теория управления - наука, изучающая законы управления, принципы управления, принципы построения систем управления, а также методы анализа и синтеза систем управления.

• Автоматика - это отрасль науки и техники, охватывающая теорию и практику автоматического управления.

• Система - это целенаправленное множество взаимосвязанных элементов любой природы.

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

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

Необходимыми условиями управления являются наличие объекта управления и цели управления. Задача управления заключается в формировании требуемого для достижения цели закона управления.

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

• Объект управления - это часть окружающего нас мира, на которую можно воздействовать с определенной целью.

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

В зависимости от участия человека в процессе управления последнее подразделяется на ручное, автоматическое и автоматизированное.

Водяные часы Ктесибиоса

Центробежный регулятор Уатта (1769 г.)

Поплавковый регулятор уровня воды

Перспективы развития роботов и систем управления

Принципы управления

• Принцип разомкнутого управления;

• Принцип замкнутого цикла

(обратной связи);

• Комбинированный;

• Принцип адаптации.

Принцип разомкнутого управления

ОУ - объект управления; ЗУ - задающее устройство; R - регулятор.

Координаты (переменные) системы:

g(t) - задающее воздействие; y(t) - управляемая (регулируемая) величина; f(t) - возмущающее воздействие; u(t) - управляющее воздействие.

Принцип замкнутого цикла (обратной связи)

ОУ - объект управления; ЗУ - задающее устройство; ИПУ - измерительнопреобразовательное устройство;

СУ - сравнивающее устройство; R - регулятор.

Координаты (переменные) системы: g(t) - задающее воздействие; y(t) - управляемая (регулируемая) величина;

f(t) - возмущающее воздействие; x(t) - рассогласование (ошибка); u(t) - управляющее воздействие.

Комбинированный принцип

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

• КЦВ - компенсирующая цепь по возмущающему воздействию, позволяет скомпенсировать негативное влияние возмущающего воздействия на работу системы.

Принцип адаптации

Адаптивные системы разделяются на:

• экстремальные, • самонастраивающиеся;

• самоорганизующиеся.

Функциональная схема экстремальной системы

УАПЭ - устройство автоматического поиска экстремума, которое анализирует параметр объекта управления ?, определяющий его режим работы, и воздействует на задающее устройство с целью изменения задающего воздействия g(t) для обеспечения наивыгоднейшего режима работы объекта управления.

Классификация систем управления

1. По принципу действия:

а) разомкнутые системы;

б) замкнутые системы;

в) комбинированные системы;

г) адаптивные системы

2. По виду задающего воздействия g(t):

а) системы стабилизации, если g(t)=const;

б) системы программного управления, если g(t) - предварительно заданная функция времени;

в) следящие системы, если g(t) - случайная величина.

3. По математическому описанию:

а) линейные системы;

б) нелинейные системы.

4. По характеру передачи сигналов:

а) непрерывные системы, такие, у которых все координаты или переменные являются непрерывными функциями времени;

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

5. По реакции системы на входное воздействие:

а) детерминированные системы - это системы, отвечающие на один и тот же входной сигнал всегда одним и тем же вполне определенным выходным сигналом;

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

в) стационарные системы - это системы, реакция которых не зависит от момента времени подачи входного воздействия;

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

6. По виду используемой энергии:

а) электрические системы, обладают удобством и легкостью обработки и передачи информации;

б)пневматическиесистемы,используютэнергиюсжатогогазаи обеспечивают высокое быстродействие;

в)гидравлическиесистемы,используютэнергиюжидкостии обеспечивают высокую мощность;

г) электропневматические системы;

д) электрогидравлические системы.

7. По числу управляемых величин.

а) одномерные системы, имеют одну управляемую величину;

б) многомерные или многосвязные системы - это системы, имеющие много входов и выходов.

8. По точности в статическом режиме.

Статически и астатические

Функциональная схема типовой автоматической системы

МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ

ЭЛЕМЕНТОВ И СИСТЕМ УПРАВЛЕНИЯ

Поведение координат системы во времени называется динамикой системы.

Математическое описание элемента устанавливает связь во времени между его текущими значениями выходных y(t) и входных x(t) величин.

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

Привремени,стремящемсякбесконечностиt??, текущиекоординатыy(t)иx(t)принимаютпостоянные установившиеся значения и наступает статика элемента, которая описывается алгебраическими уравнениями.

Встатикевходныеивыходныевеличиныэлемента постоянные:x(?)=x0=const;y(?)=y0=const.Этипостоянные величиныназываютсяустановившимися.Апроцесс, соответствующийстатике,называетсяустановившимся процессом.

Линеаризация

Большинство реальных элементов имеют нелинейные характеристики и описываются нелинейными ДУ. Многие нелинейные элементы можно линеаризовать, то есть заменить нелинейные уравнения элемента приближенными линейными.

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

Текущее значение координат y и x

y(t) = y0 + ?y(t);

x(t) = x0 + ?x(t); (1)

где:y0, x0 - установившиеся значение,

?y, ?x - отклонения от установившихся значений. В рабочей точке ( x0, y0),заменим участок кривой касательной

y = yн + kx ,(2)

где: yн - постоянная величина, k - коэффициент, определяемый наклоном касательной к кривой в рабочей точке (x0, y0).

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

?y(t) = k ?x(t)

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

Формы записи линеаризованных уравнений

an y?n? ?an?1y?n?1? ???a1y?1? ?a0y?

?bmu?m? ?b?m?1?u?m?1? ???b1u?1? ?b0u?(3)

?crv?r? ?c?r?1?v?r?1? ???c1v?1? ?c0v

Это вход-выходное дифференциальное уравнение, в котором u,v является входами, а y(t) - выходом, a, b, c - коэффициенты, могут быть нулевыми.

При исследовании линейных систем широко применяются так называемые операторные методы. Из операторных методов наиболее известны:

1) Преобразование Лапласа:

? X?p?? L?x?t????x?t?e?p?tdt

0 2) Преобразование Фурье:

? X? j??? ? ? x?t?e? j??tdt

?? Применяем преобразование Лапласа с учетом нулевых начальных условий.

(an pn ? an?1pn?1 ??? a1p ? a0)?Y(p) ?

? (bm pm ?b?m?1? pm?1 ???b1p ?b0)?U(p)?(4.1)

?(cr pr ?c?r?1? pr?1 ???c1p ?c0)?V(p)

nmr (?ai pi )?Y(p) ? (?bi pi )?U(p) ? (?ci pi )?V(p) (4.2) i?0 ? i?0 ? i?0 ?

A(p)B(p)C(p)

A(p)Y(p)?B(p)U(p)?C(p)V(p)

при нулевом возмущающем воздействии:

A(p)Y(p)?B(p)U(p)Y(p) ? B(p) ?W(p)-функцияпередаточная

U(p) A(p)

Передаточной функцией объекта или системы называется отношение изображения по Лапласу выходной переменной к входной при нулевых начальных условиях.

Свойства передаточной функции динамического звена.

1) Передаточная функция линейного звена является дробно-рациональной, т.е. отношением полиномов к комплексной переменной Р (или S).

2) Для физической реализуемости (реальные объекты и системы) порядок полинома в числителе не должен превышать порядок полинома в знаменателе.

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

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

A(p) ?anpn ?an?1pn?1 ?.....?a1p?a0 ?an(p?p1)(p?p2).....(p?pn)

p1 ? 0p2 ??bp3,4 ???? j?

A(p) ? a4p(p ?b)(p2 ?2?p ??2 ??2)

vks A(p) ?an p??(p?pi )?(p2 ?2?i p??j2 ??j2)

i?1j?1

? - количество нулевых корней в полиноме? - количество вещественных корней

ks - количество пар комплексно - сопряженных корней

Произвольная форма записи

bmpm ?b(m?1) pm?1 ?.....?b1p?b0

W(p) ? a pn ?a(n?1) pn?1 ?.....?a1p?a0

n Каноническая форма записи

W (p) b0 p???bm*pm ?b(m?1)*pm?1 ?.....?b1*p?1

K?a0an*pn ?a(n?1)*pn?1 ?.....?a1*p?1

ai* ?aa0ibi* ?bb0iab00 ?k -(коэффициент передачи усиления при к>0)

? - количество нулевых корней в числителе? - количество нулевых корней в знаменателе

Последовательно-структурированная форма записи

v1ks1

?(Ti p?1)?(Tj 2p2 ?2?Tj p?1)

WК ПС (p) ?kp??? iv?21ksj?21

?(Ti p?1)?(Tj 2p2 ?2?Tj p?1)

i?1j?1

Ti, j- постоянные времени элементарных звеньев 1

Ti ? (?pi ) Tj ?

?? 1 22 ?j ??j

??j

22 ?j ??j ? - Коэффициент демпфирования

??v - характеризует порядок астатизма сложного динамического звена; если 0, то звено является статическим, если больше 0, то звено обладает астатизмом.

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

Приведенная форма записи:

bm p??? ppmn ??ab(m?1)**ppnm??11??..........??ab11**pp??ab00**

WПР(p) ? an(n?1)

bamn ?kПРbi* ?bbmiai* ?aani

Такая форма записи используется для перехода в матрично-векторную форму.

Условие реализуемости передаточной функции: m?n

Типовые входные воздействия

Характеристика звена - это его реакция на определенное входное воздействие.

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

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

1) дельта-функция ?(t);

2) единичная ступенчатая функция 1(t);3) гармоническая функция X0sin(?t). ?-функцией называется функция, равная нулю всюду, кроме начала координат, но притом так, что интеграл от нее по любому интервалу, содержащему нуль, равен единице, т.е.

?(t)= ??0 при t ? 0,Кроме того, ? при любом ?>0. ?? при t = 0. ???? (t)dt =1

Весовая функция w(t) представляет собой реакцию на выходе звена при подаче на его вход ?-функции.

Единичная ступенчатая функция удовлетворяет следующему условию

?0 при t ? 0, 1(t)= ?

?1 при t ? 0.

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

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

Если на вход линейного звена подать гармоническое воздействие x(t)=X0sin(?t),

гдеX0 - амплитуда,

? - угловая частота, имеющая размерность [рад/с] или [c-1 ], то, на выходе звена в установившемся режиме будет также гармоническая функция той же частоты, но, в общем случае, другой амплитуды Y0 и сдвинутая по фазе относительно входной величины на угол ?

y(t)=Y0sin(?t+?).

Связь между выходной гармоникой и входной устанавливается с помощью частотной передаточной функции звена W(j?).

Частотная передаточная функция W(j?), как видно, представляет собой комплексное число, которое можно записать как в полярной, так и декартовой системах координат:

W(j?) = A(?) ej?(?) = U(?) + jV(?), где А(?) - модуль или амплитуда частотной передаточной функции,

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

???? - аргумент или фаза частотной передаточной функции, показывает фазовый сдвиг выходной гармоники по отношению к входной на частоте ?

???? = arg W(j?);

U(?) - вещественная составляющая частотной передаточной функции U(?) = Re W(j?);

V(?) - мнимая составляющая частотной передаточной функции V(?) = Im W(j?).

22(?) = arctg V(?) A(?) =U (?)? V (?)?

U(?)

Частотные характеристики

Амплитудно-фазовая частотная характеристика (АФЧХ). Строится на комплексной плоскости и представляет собой геометрическое место конца вектора, соответствующего частотной передаточной функции W(j?) при изменении ? от 0 до ?.

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

где ?р - резонансная частота, т.е. частота, на которой амплитудная частотная

характеристика достигает максимума,

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

?п - частота пропускания, частота, на которой амплитудная частотная характеристика, уменьшаясь, принимает значение, равное 0,707, и при дальнейшем повышении частоты не увеличивается;

??п=2?п - полоса пропускания, диапазон частот гармонических колебаний, пропускаемых звеном без заметного ослабления.

Фазовая частотная характеристика (ФЧХ). Показывает фазовые сдвиги, вносимые звеном на различных частотах.

Вещественная частотная характеристика (ВЧХ). Мнимая частотная характеристика (МЧХ)

Логарифмические частотные характеристики (ЛЧХ).

L(?) = 20 lg A(?) = 20 lg?W(j?)????? = argW(j?)

Типовые динамические звенья

Пропорциональное (безинерционное) звено

Описывается алгебраическим уравнением y(t) = k·x(t).

Передаточная функция:

W(p) = k.

Амплитудно-фазовая частотная характеристика:

W(j?) = k, A(?) = k, ?(?) = 0.

Интегрирующее звено

kkjk W(p) ?? W( j?) ??

T ? pTj? ??T kk?

U(?) ? 0; V(?) ? ?;A(?) ?; ?(?) ? arctg(??) ? ?

?T?T2

k L(?) ? 20 lg ?T

При увеличении частоты на декаду оцениваем наклон ЛАЧХ: k1kk L(10??) ? 20 lg? 20 lg?20 lg? ?20?20 lg

10?T10?T?T Наклон ЛАЧХ: ?20дБ дк

ЛАЧХЛФЧХПереходной процесс

Дифференцирующее звено

W(p) ? kTp? W( j?) ? jkT?

U(?) ? 0; V(?) ? kT?;A(?) ? kT?; ?(?) ? arctg(?) ??

2 L(?) ? 20 lg(kTp)

При увеличении частоты на декаду оцениваем наклон ЛАЧХ:

L(10??) ? 20 lg10kT?? 20?20 lgkT?

Наклон ЛАЧХ: ?20дБ

дк ЛАЧХЛФЧХПереходной процесс

Апериодическое звено

W(p) ? Tpk?1? W( j?) ? Tj?k?1 ? (1? kj?(1T?)(j1??Tj)?T) ? k1???jk2?TT2

U(?) ? 1??k2T2 ; V(?) ? ?1?k??T2T2 ; ?(?) ? arctg(??T)

при? ??(?) ? ?45

С10; ?С ?сопрягающая частота

при при ????C

????C W( j?) ? k

k W( j?) ? j?T ЛАЧХ ЛФЧХ Переходной процесс T

L, дБ С Пропорционально-дифференцирующее звено

W(p) ? k(Tp?1) ? W( j?) ? kTj??k

U(?) ? k; V(?) ? kT?; ?(?) ? arctg(?T)

ЛАЧХЛФЧХПереходной процесс

Колебательное звено k

W(p)?T 2 p2 ?2?Tp?1 ? два комплексно сопряженных корня

ЛАЧХЛФЧХПереходной процесс

Пропорционально-дифференцирующее 2-го порядка

W(p)?k(T2p2 ?2?Tp?1)

ЛАЧХЛФЧХПереходной процесс

Правила преобразования структурных схема

Виды соединений звеньев

Всистемахуправлениявыделяюттривида соединений звеньев:

- последовательное

- параллельное

- с обратной связью (положительной или отрицательной).

Последовательное

Обозначим промежуточные выходы

W1(p) ? у1(р) , W2(p) ? у2(р) , W3(p) ? у(р) , х(р)у1(р)у2(р)

WЭ (p) ? у(р) ?W1(p)?W2(p)?W3(p) ? ?n Wi (p) х(р)i?1

где n - число последовательно включенных звеньев

Параллельное

ху у(р)?у1(р)?у2(р)?у3(р)

W1(p) ? у1(р) , W2(p) ? у2(р) , W3(p) ? у3(р) y(р)у(р)у(р)

WЭ (p) ? у(р) ? у1(р)?у2(р)?у3(р) ? х(р)х(р)

n ?W1(p)?W2(p)?W3(p) ??Wi (p)

i?1 где n - число параллельно включенных звеньев

Обратная связь

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

связи.

W1(p) ? у(р) , W2(p) ? у1(р) , е(p) ? х(р)? у(р) е(р)у(р)

WЭ(p) ? у(р) ? W1(р)?e(p)

х(р) e(р)? y1(p)

WЭ(p) ? 1?W1у(1p()р) ? 1?W1W(p1()p?W) 2(p) e(р)

где знак "+" соответствует ООС, знак "?" - ПОС.

W234 ?W2(W3 ?W4)

W2345 ?W2(W3 ?W4)

1?W2 ?W5 ?(W3 ?W4)

WY ?W1? W2(W3 ?W4)

1?W2 ?W5 ?(W3 ?W4)

Основные правила эквивалентных преобразований структурных схем

- перенос сумматора через звено - перенос точки съёма через звено

- перенос сумматора через сумматор

- перенос точки съёма через точку съёма

Сумматор через звено

у 1.А yA ? (x1 ? x2)?W ? x1?W ? x2 ?W

yA ?yB при

2. yB ? x1?W ? x2 ? A

A?W ху 2 yA ? x1?W ? x2yA ? (x1 ? x2 ? A)?W ? x1?W ? x2 ? A?W

1 yA ? yBi?eA? W

Точка съема через звено

уху 1. 11 y ? x?Wx ? y??A ?

WW у 2.

1 y ? x?Wx ? y??A ?W

W Сумматор через сумматор (точка съема через точку съема)

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

? ?

Аналогично для точек съема стоящих рядом

? ? WWЭ ? W1 ?(W2 ?W3)

Э Построение логарифмических частотных характеристик

(ЛЧХ) сложных динамических звеньев

Построение логарифмических частотных характеристик (ЛЧХ) сложных динамических звеньев.

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

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

W ПС (p) kp????W (p)

K?i i?1 Построение ЛАЧХ сложного динамического звена.

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

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

Построение ЛФЧХ сложного динамического звена.

Строится ЛФЧХ элементарных динамических звеньев, входящих в сложное по

отдельности. После этого они суммируются и формируют ЛФЧХ сложного звена.

При построении ЛАЧХ важно:

1) Правильно строить низкочастотную асимптоту;

2) При построении двигаться строго слева направо.

Правило построения низкочастотной асимптоты:

20*lg(kp???) ? 20*lg(k)?(???)20*lg( p)

Таким образом на оси частот откладывается точка соответствующая нулю, на оси

ординат откладывается 20*lg(k) и через эту точку проводится линия с наклоном 20(???).

Правило построения среднечастотной асимптоты:

Определяется сопрягающие частоты всех элементарных звеньев входящих в сложное;

они располагаются на оси частот в порядке возрастания, через эти частоты проводятся ординаты пунктиром. Построенная низкочастотная асимптота проводится до точки пересечения с первой ординатой - это будет первая точка поворота ЛАЧХ. Далее определяется звено, которое соответствует этой частоте и в соответствии с его характеристиками изменяется наклон ЛАЧХ сложного звена; она продолжается до точки пересечения со следующей ординатой.

Правило построения высокочастотной асимптоты:

Последняя построенная асимптота продолжается до бесконечности.

Пример.б)

W ( p) ?; 1) k=40; ?? 0; ?? 0

20lg(40) ? 32; 20(???) ? 0дб / дк

в) 2)а) W1(p) ?5p?1;

11 ?C1 ? ? ? 0.2

T 5 lg(0.2) ??0.7

W2(p) ?(КЗ) T? 4?2

1 ?C2 ? ? 0.5 2 lg(0.5) ??0.3

2?T ? 3; 4?? 3 ??? 0.75

1 W3(p) ?(АЗ)

0.2p?1 1 ?C3 ?? 5

0.2 lg(5) ? 0.7

Устойчивость систем. Критерии устойчивости.

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

WЗСZY ? WPWOYWЗСZE ?1H(p) ?1?WPWOY

1?WPWOY1?WPWOY

- характеристическая функция.

L(p) WPC(p)?WPWOY ?

S(p) где S(p) - характеристический полином

разомкнутой системы. S(p)

H(p) ??L(p)

S(p) WЗСZY ? L(p)/S(p) ?L(p)

L(p)S(p) ?L(p) 1? N(p) S(p)

В характеристической функции в числителе стоит характеристический полином замкнутой системы, а в знаменателе характеристический полином разомкнутой системы.

Понятие устойчивости.

Методы оценки устойчивости объектов и систем.

Корневой критерий (Ляпунова)

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

1 - устойчивые корни

2 - нулевой корень

3 - нейтральные корни

4 - неустойчивые корни

Качественная оценка характеристического полинома.

Необходимое условие устойчивости всех корней характеристического полинома:

1) Все корни характеристического полинома должны быть одного знака; 2) Характеристический полином должен быть полным, т.е. содержать все нулевые коэффициенты, т.е. для каждого аi ? 0 i ? 0,n

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

Все методы оценки устойчивости делятся на 2 основных класса: а) Алгебраические (с помощью алгебраических вычислений)

б) Частотные (используя частотные характеристики)

Алгебраические критерии:Раусса

Гурвица КритерийГурвица

В 1895 году новый критерий изобрел немецкий математик Гурвиц.

N(p) ?a0pn ?a1pn?1 ?.....?an?1p?an

Для оценки устойчивости также используется характеристический полином. На его основе составлен главный определитель Гурвица a1 a300Для проверки правильности

составления определителя,

?n ?a0 a2......анализируется главная диагональ, где должны стоять все коэффициенты

0 a1...0характеристического

полинома от а1 до аn

0 a0... an Критерий устойчивости Михайлова:

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

Пример: устойчивые

Пример: неустойчивые

Пример: на границе устойчивости

Модифицированный критерий Михайлова.

Обычно применяется при построении и анализе годографа Михайлова в математических пакетов. N( j?)

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

Критерий Найквиста: Для оценки устойчивости используют АФЧХ разомкнутой системы.

L(p)L(p)?S(p)

WPC(p) ? S(p) H(p)?1?WPC(p)?S(p)?0

S(p) - характеристический полином разомкнутой системы

L(p)?S(p)- характеристический полином замкнутой системы

nPC ?nЗС

Графическая трактовка критерия Найквиста.

1) Разомкнутая система является устойчивой

Для устойчивости замкнутой системы необходимо и остаточно, чтобы годограф АФЧХ не охватывал точку с

координатами (-1, j0). Для этого применяется эмпирическое правило правой штриховки: на годограф АФЧХ наносится штриховка справа по мере увеличения частоты, если точка (-1 , j0) попадает в область штриховки, значит годограф охватывает ее.

2)Разомкнутаясистема обладаетастатизмом

первого порядка 3) Разомкнутая система обладает астатизмом произвольного порядка (допустим третьего) q?1 1 - устойчива q?3

2 - неустойчива

Понятие запасов устойчивости замкнутой системы.

Критерий Найквиста имеет некоторое преимущество над другими, анализируя АФЧХ разомкнутой системы можно оценить количественно близость устойчивости системы к границе неустойчивости и наоборот. Она определяется с помощью специальных запасов устойчивости.

?m - запас по амплитуде или по модулю

?? - запас по фазе

Трактовка критерия Найквиста с использованием ЛЧХ.

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

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

?L -запаспоамплитуде(по

модулю) (в дб.)

?F - запас по фазе, в град. (рад.)

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

Связь между качеством системы и значением запасов - не жесткая, но величина запасов зависит от порядка разомкнутой системы и степени ее астатизма. Чем больше порядок астатизма и размерность разомкнутой системы, тем выше запас устойчивости.

Статические системы управления(?? 0)

Степень ХП ?F (град) ?L (дб) n?3 35-45 7-10 3 ?n? 6 45-60 10-15 n ? 6 >60 >15 Астатические первого порядка (? ?1)

Степень ХП ?F (град) ?L (дб) n? 3 45-60 10-15 3?n?6 60-90 15-20 n ? 6 >90 >20 Астатические второго порядка(?? 2)

Степень ХП ?F (град) ?L (дб) n? 3 60-90 15-20 3?n?6 90-120 20-25 n? 6 >120 >25 Оценка качества управления в линейных САУ.

Основные понятия и определения.

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

К качественным показателям относят такие оценки как:

- быстродействие;

- колебательность;

- астатизм ...

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

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

В теории управления применяют два основных подхода для оценки качества:

- Прямая оценка качества (осуществляется по переходному процессу);

-Косвенная оценка качества (осуществляется по другим характеристикам, например, частотным или корневым).

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

Вид переходного процесса зависит не только от динамических и

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

1) Переходной процесс с

перерегулированием, где ?

- величина перерегулирования.

tH - время нарастания переходного процесса; tMAX - время достижения первого максимума;

?? hMAX ? hУСТ ?100% hУСТ

В реальных системах ?? 35%

2) Апериодический переходной процесс

?доп - допустимая трубка регулирования (по сути - допустимое отклонение переходного процесса от hУСТ

?доп = 1% - 10%

Самые распространенные - 3% и 15% tp1,2- время регулирования переходного процесса

A1 Косвенный метод оценки качества

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

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

2 ?sin ? t

y(t)= h(t)= ? P(?) d?

? 0? где P(?) - вещественная частотная характеристика замкнутой системы;

Корневые оценки качества

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

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

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

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

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

Такая степень устойчивости называется колебательной.

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

Колебательность

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

Именно эта величина характеризует быстроту затухания колебаний за каждый период T=2?/?. Чем выше колебательность, тем слабее затухание колебаний в переходном процессе.

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

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

где ?з и ?з ? заданные значения степени устойчивости и

колебательности

Оценка точности работы систем.

Коэффициенты ошибок.

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

E(p)=Wg?(p)?G(p) =1+W(1 p)G(p)

Где Wg? - передаточная функция замкнутой системы по каналу ошибки, G(p) - изображение задающего воздействия.

Разложив выражение в ряд по возрастающей степени комплексной величины p

E(p) ? ???c0 ?c1p ? c22! p2 ?...????G(p)

C0, C1 - коэффициенты ошибок

? i 1 ?d Wg?(p)??Коэффициенты c0, c1 и c2 называются соответственно

ci ? i!??dpi??коэффициентами ошибки по положению, скорости и

??p = 0ускорения.

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

2 dg(t)d g(t)

E(t)= c0g(t)+c1?c22 ?...

dtdt Статическая система:с0?0, с1?0, с2?0.

Астатическая 1 порядка: с0=0, с1?0, с2?0.

Астатическая 2 порядка: с0=0, с1=0, с2?0. Астатическая 3 порядка: с0=0, с1=0, с2=0.

Построение областей устойчивости САУ в плоскости их параметров.

Для построения областей устойчивости применяется метод предложенный Неймарком и названным - метод D - разбиения.

Исходными данными для построения областей - это характеристический полином замкнутой системы.

Re[N(?,{?})]? Im[N(?,{?})] ? 0

?Re[N(?,{?})] ? 0

??{?} ?Im[N(?,{?})] ? 0

Однопараметрическое D-разбиение.

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

N(p,?)?S(p)??Q(p)?0

S(j?)??Q(j?)?0 S(j?) Re?(?)? jIm?(?)

???? Для определения области устойчивости САУ напараметрическуюкривуюнеобходимо нанести штриховку слева по мере увеличения частоты (по аналогии с корневой плоскостью).

Примечание: для построения параметрического годографа необязательно изменять частоту от ??до ??. Достаточно рассмотреть ?от 0 до ??,

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

Двухпараметрическое D-разбиение

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

N(p,?,?)??K(p)??L(p)?M(p)?0 ??K1(?)??L1(?)??M1(?)

p? j???K2(?)??L2(?)??M2(?)

? ????KL((jj??))??LK1(1?(?)?)?jLjK2(?2(?) ) ??KK12((??)) LL12((??))?1???MM12((??)) LL12((??));

?M(j?)?M1(?)? jM2(?)

?K(j?)??L(j?)?M(j?)?0 ?2?K1(?) ?M1(?)???1; ???2

K2(?) ?M2(?)?? Так как M2 L2 и K2 являются нечетными функциями, то и ?,?1 и ?2 также нечетные функции, а функции ? и ? будут четными. При построении параметрических кривых необходимо соблюдать следующее: по оси абсцисс располагается параметр ?, а по оси ординат - ?.

Рассмотрим три возможных случая решения систем уравнений.

Уравнения совместны и Уравнения несовместны и не имеют общее решениеимеют общее решение

?K(?k)??L(?k)?M(?k)?0-

- уравнение особой прямой

В большинстве практических случаев особые

прямые получаются при ??0 и ???.

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

1. (Для параметрической кривой) При возрастании ? от ?? до ?? граница

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

2. (Для особых прямых) Штриховка особых прямых, как правило, одинарная и и производится так, чтобы в точке сопряжения особой прямой и параметрической

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

После нанесения штриховки определяется область, претендующая на область устойчивост

Области, претендующие на

устойчивость - I и II. Область II после проверки будет действительно являться областью устойчивости.

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

Регуляторы. Законы регулирования.

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

Регуляторы подразделяются

1) по виду регулируемой величины: давление, температуры, переходного уровня и т.д.

2) по виду источника:

- прямого действия - не требуется дополнительной энергии (сливной бочок) - непрямого действия - требует дополнительной энергии (электрические, гидравлические, пневматические др.)

Типовые законы регулирования

Существует 5 типовых регуляторов. Среди этих регуляторов есть простые и сложные.

1) П - пропорциональный регулятор, W(p)=k (является простым, т.к. один параметр настройки - k). Является статическим регулятором.

2) И - интегрирующий регуляторW(p) ? 1

TИp

Так же является простым, но в отличие от первого он является астатическим.

1 3) ПИ - пропорционально-интегрирующий регуляторW(p) ?k?

TИp Является сложным и астатическим.

4) ПД - пропорционально-дифференцирующий регулятор WИ(p) ? k ?TД p

W P(p) ? kTp?T?Д1pT??TД- является сложным и статическим.

T ?0 5) ПИД - пропорционально-интегрально-дифференцирующий регулятор

1 W(p) ?k?TД p?- является сложным и астатическим

TИ p 1) П - пропорциональный регуляторW(p) ? k

Недостатки: Если объект статический, то П - регулятор не повысит порядок астатизма.

Является достаточно эффективным для астатических объектов.

2) И - интегрирующий регуляторW(p) ? 1

TИ p (применяется для статических систем)

Достоинства: Позволяет повысить порядок астатизма и прост в настройке;

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

3) ПИ - пропорционально-интегрирующий регулятор1kT p?1

Достоинства: Позволяет повысить порядок астатизма, при одновременном уменьшении

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

Применяется для статических систем и является самым распространенным в инженерной практике.

4) ПД - пропорционально-дифференцирующий регуляторW И (p) ?k?TД p

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

Недостатки: - Не повышает порядок

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

5) ПИД - пропорционально-интегрально-дифференцирующий регулятор

1 TИTД p2 ?kTИ p?1 W(p) ?k?TД p? ?

TИ pTИ p

Корректирующие устройства

Корректирующие устройства

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

Различают три вида основных корректирующих устройств.

1.Последовательные корректирующие устройства.

2.Параллельные корректирующие устройства 3.Обратные связи

1. Последовательные КУ

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

2. Параллельные КУ

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

3. Обратные связи

Вводятся в цепь регулятора и охватывают какие-либо его звенья. Могут быть положительными (ПОС) и отрицательными (ООС), кроме того - жесткими и гибкими.

Основнымивидамикорректирующихобратныхсвязей

являются:k

а) жесткая обратная связь WОС(p) = kОС;ос

б) инерционная жесткая обратная связь WОС(p) =T s+1 ос

в) гибкая обратная связь WОС(p) = kОС p;kосs

г) инерционная гибкая обратная связь WОС(p) = Ts+1 ос

Синтез последовательного КУ

Алгоритм синтеза

• Строится ЛАЧХ исходной системы;

• Строится желаемая ЛАЧХ;

• Получим ЛАЧХ КУ вычитая их желаемой характеристики исходную.

При построении желаемой ЛАЧХ выделяют три частотные области:

1. низких частот (НЧ); 2. средних частот (СЧ); 3. высоких частот (ВЧ)

Низкочастотная часть желаемой ЛАЧХ формируется из условия обеспечения требуемой точности работы системы1 управления в установившемся режиме, то есть из условия?? того,чтоустановившаясяошибкасистемынедолжна1?К превышать заданное значение.

Среднечастотныйучасток желаемой ЛАХ строится из условия обеспеченияосновныхпоказателей качествапереходногопроцессаперерегулированияивремени регулирования. Это достигается тем, что он имеет наклон ?20 дб/дек и пересекаетосьчастотначастоте среза ?с, которая определяется по номограммеВ.В.Солодовникова,

исходяиззаданныхзначений???

величиныперерегулирования??и?срж ?

времени регулирования tр.t

p Протяженность среднечастотного участка определяется номограммой, устанавливающей связь между показателями качества и запасами устойчивости. Так, например, для обеспечения ?=25% требуется запас устойчивости по модулю Lh=20 дб и по фазе ?=580. Среднечастотный участок проводится с наклоном ?20 дб/дек влево и вправо от частоты среза ?с до достижения модулей, равных Lh и ?Lh. После этого участки средних и низких частот сопрягаются прямой с наклоном ?40 или ?60 дб/дек, как удобнее.

Анализ систем в пространстве состояний

Современная теория управления основана на описании систем в

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

Понятие состояния является определяющим в современной теории

управления.

Под состоянием системы понимается минимально-необходимый набор

переменных величин системы x1,x2,...,xn, способных однозначно и единственным образом определить положение системы в любой момент времени t.

Совокупность переменных величин x1,x2,...,xn образует n-мерное

пространство состояний Rn. Вектор с компонентами x1,x2,...,xn называется вектором состояния.

Рассмотрим систему с m входами (u1,u2,...,um), r выходами (y1,y2,..., yr) и n переменными координатами (x1,x2,...,xn).

Линейную систему управления можно описать с помощью системы дифференциальных уравнений вида

??x?1 ? a11?x1 ? a12?x2?...? a1n ?xn ?b11?u1 ?b12 ?u2...?b1m ?um

?? ?x2 ? a21?x1 ? a22?x2?...? a2n ?xn ?b21?u1 ?b22?u2...?b2m ?um ?

?......... ?? ??xn ? an1?x1 ? an2 ?x2?...? ann ?xn ?bn1?u1 ?bn2 ?u2...?bnm ?um

Кроме того, для описания выходной переменной используем уравнение

y1 ? c11?x1 ?c12?x2 ?...?c1n ?xn ? d11?u1 ? d12?u2 ?...? d1m ?um y2 ? c21?x1 ?c22?x2 ?...?c2n ?xn ? d21?u1 ? d22?u2 ?...? d2m ?um

......

yr ? cr1?x1 ?cr2 ?x2 ?...?crn ?xn ? dr1?u1 ? dr2 ?u2 ?...?drm ?um

Рассматриваемаясистемаможетбытьопределенаследующейвекторноматричной формой

?X? =AX+BU ? ?Y=CX+DU

где X - вектор состояния системы, Y - вектор выходных управляемых величин, U вектор внешних воздействий (задающих и возмущающих),A,B,C,D - матрицы, а именно:

?x ??y ?? u ??aa...a?

?x1 ? ?y1 ? ? u1 ? ?a11 a12 ... a1n ? X = ? 2 ? Y = ? 2? U = ? 2 ? A = ? 21 22 2n?

? ... ?? ... ?? ... ?? ............ ?

?b 11 ?b B = ? 21

? ...

?b ?? n1 b 12 b 22

... b n2 ... ...

... ... b? 1m ?

b 2m?

... ? b? nm??n x m ?d

11 ?d D = ? 21

? ... ?d ?? r1 d

12 d 22

... d r2 ... ...

... ... d? 1m ?

d 2m? ... ? d?

rm ??r x m ? n1n2

?cc , 1112

?cc C = ? 2122 ? ......

?cc , ?? r1r2 ...

... ...

... nn?n x n c ?

1n ? c 2n? ... ? c ?

rn ??r x n ???xn ??????yr ??????um ?????aa...a??

A - системная матрица,

В - матрица управления,

С - матрица выхода по состоянию

(наблюдения),

D - матрица выхода по управлению (связи).

Структурная схема систему управления в пространстве состояний

Построение МВ ММ по структурной схеме.

Алгоритм построения:

Исходная математическая модель, представленная в графической

1) форме преобразуется к структуре, содержащей только элементарные звенья нулевого или первого порядков (ПЗ или ИЗ).

На выходах динамических звеньев первого порядка назначаются

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

Формируются сигналы на входах динамических звеньев, согласно

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

Составляется система операторных алгебраических уравнений, после

4) применения обратного преобразования Лапласа она преобразуется в систему дифференциальных уравнений в форме Коши, по которой строится МВ ММ.

ПримерДана структурная схема объекта управления

Необходимо построить МВ ММ

?px1 ?10(U?15x2)

??px2 ??0,1x2 ? 2x1

? ?px3 ? 2x2

??y? 5x1?x3 ? 0 -150 0?

A?? 2 -0,1 0?

?? ?? 0 2 0??

?x? ? Ax?Bu

? ?y?Cx ?10? B?? 0 ?C??5 0 1?

? ? ?? 0 ?? Построение вход-выходных моделей по

матрично-векторным

?x? ? Ax?Bu?px(p)? Ax(p)?Bu(p) ??y?Cx?Du??Y(p)?Cx(p)?Du(p)

(pE?A)x(p) ?Bu(p)

????? ? характеристическая матрица

det(pE?A) ?N(p) ? характеристический полином x(p)?(pE?A)?1Bu(p)

???????

WПМux ?(p)?n?m

Передаточная матрица (ПМ) по состоянию

Y(p) ?Cx(p)?Du(p)

Y(p)?C[WПМux (p)u(p)]?Du(p)

Y(p)??[CW??ПМux?(p?)??D?]u(p)

? WПМuy (p)?p?m

Передаточная матрица по выходу

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

Пример

?x1 ???1p p10?1?? ??U1 ??WU1X1(p)? 1pWU1X2 (p)? 5p

?? ? ??? ???x2??????5p p?1 2??? ???U2???WU2X1(p)? p10?1 WU2X2 (p)? p1?2

Свойства системы в пространстве состояний

К фундаментальным свойствам системы в ПС можно отнести: управляемость и наблюдаемость.

Управляемость и наблюдаемость, так же как и устойчивость, относятся к числу важнейших характеристик динамических систем.

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

Система называется управляемой, если для любого начального состояния X(0)?Rn существует управление U(t), переводящее ее за конечное время T в нулевое состояние X(T)=0.

Для оценки управляемости использую матрицу управляемости

Qу ?B AB A2B ... An-1B

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

Система управляема, если ранг матрицы управляемости равен размерности системы n. Если ранг меньше n система частично управляема, равен 0 - неуправляема полностью.

Система называется наблюдаемой, если по наблюдениям за выходным сигналом Y(t) в течение конечного времени T можно определить ее начальное состояние X(0).

Для оценки управляемости использую матрицу наблюдаемости

C CA

Qн ?CA2 ...... Необходимым и достаточным условием наблюдаемости системы является невырожденность матрицы наблюдаемости.

Система наблюдаема, если ранг матрицы наблюдаемости равен размерности системы n. Если ранг матрицы меньше n система частично наблюдаема, равен 0 - ненаблюдаема полностью. CAn?1

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

A= ??10 10 11?? B= ??10??C = ?001?

?1 1 0??????0???

?? Решение. Характеристический определитель матрицы A

p ?1 ?1 pE?A = ?1 p ?1 ? p(p2 ?1)?p?1?1?p? p3?3p?2

?1 ?1 p

p3?3p?2?0 ?1 Qу=?0

? ??0 0 1 1 2? 1?

? 1?? ?0 Qн=?0

? ??1 1

1 0 1? 1? ? 2??

Система неустойчива.

Матрица управляемости

, det Qу=1?1=0, следовательно, система неуправляема.

Матрица наблюдаемости

, det Qн=1?1=0, следовательно, система ненаблюдаема.

Основные формы представления МВ ММ и способы перехода к ним

Для исследования динамических систем в пространстве состояний в теории управления используется семь основных форм: 1) Основная нормальная форма (ОН)

? 0 1 0 ... 0

? 0 0 1 ... 0

? AOH ? ?.........................................

? ? 0 0 0 ... 1

? a a a ... a

??0? ??0? ?? ?

? BOH ? ?...? COH ? ?1 0 ... 0?

????...?? ??1 ?

?? 0 ? 1 ? 2? n?1?? ?

N(p) ? pn ?an?1pn?1 ?an?2 pn?2 ?...?a1p?a0

В ОН форме управляемой матрица В имеет стандартную форму, а С может быть произвольной.

В ОН форме наблюдаемой матрица С должна быть стандартной, а В может быть произвольной.

Для перехода к ОН форме используют преобразования на основе матриц управляемости или наблюдаемости.

2) Транспонированная ОН (ТОН)

?0 0 0 ... ?a0 ??1 ?

??1 0 0 ... ?a1 ????0??

AТОН ? ?.................................... ? BTOH ? ?...? CTOH ??0 0 ... 1?

??? ?

?0 0 0 ... ?an?2??...? ??0 0 0 ... ?an?1 ????0??

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

?CAn?1 ?an?1CAn?2 ?...?a1C?

?? ?............................................?

TOH? CA2 ?an?1CA?an?2C ?

а) Р1? ?? ? CA?an-1C?

??? C???

Исходный базис должен быть в ОН - форме.

б) Р2ТОН ?Qу?1 , где Qу - матрица управляемости

3) Дополнительная нормальная форма

??an?1 1 0 ... 0?

???an?2 0 1 ... 0??

AДОП ? ?............................?

?? ВДОП ?ВОН СДОП ?СОН

? ?a1 0 0 ... 1 ? ?? ?a0 0 0 ... 0??

4) Обратная нормальная форма

??an?1 ?an?2 ... ?a0?

? 1 0 ... 0 ? ??

AОБР ??................................... ? ? 0 0 . .. 0 ?

? 0 0 . .. 0 ??? ВОБР ? ВТОН СОБР ?СТОН

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

5) Каноническая диагональная форма (КД)

??1 0 ... 0 ??1 ??b1 ?

?? 0 ?2 ... 0 ????1 ????b2 ??

AКД ? ?.................... ?BКД ? ?...? ???....?? CКД ??b1 b2 ... bn?

???...??....? ?1 1 ... 1 ?

??? 00 00 ...... ?0n??????bn ??

??1 ???

где ?i - собственные числа матрицы А

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

? 1 1 1 ... 1 ? ??? ?2 ?3 ... ?n ?? 1

W ? ???12 ?22 ?32 ... ?2n ?? PКД ?W?1

?.......................................?

??n?1 ?n2?1 ?3n?1 ... ?nn?1??

? 1 Для перехода к диагональной форме записи необходимо сначала преобразовать к ОН форме. 6) Жарданово - каноническая форма

?A1 0 ... 0 ?

?? 0 A2 ... 0 ??

AЖК ? ??...................... ?? ВЖК ?ВКД СЖК ?СКД ,где Ai - клетки Жордана

? 0 0 ... 0 ? ?? 0 0 ... An??

7) Треугольная формаИз данных 7 форм, наибольшее

??1 a12 ... a1,n ?распространение получили три:

?? ? 0 ?2 ... a 2,n ?а) ОН

AТР ? ??.................... ...........? ВТР ?ВКД СТР ?СКДб) ТОН

? ? 0 0 ... a n-1,n-1 ?в) КД

?? 0 0 ... ?n ??

Понятие базиса в пространстве состояний

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

Форма записи матрично-векторных моделей не является единственной для каждых конкретных A,B,C, конкретные A,B,C соответствуют некоторому конкретному базису в пространстве состояний. Любое n-мерное линейное пространство имеет множество различных базисов.

?x? ? Ax?Bu

? ?y?Cx?Du Для перехода к новому базису существует некоторое преобразование, которое обозначается: .x* ? Px

Это преобразование должно удовлетворять следующим условиям: - p -

квадратная матрица и невырожденная [ ].rang(p) ? 0

? x?Р?1x*????Рy??1x?C*Р??1AxР* ??1xDu* ? Bu?????xy?*??CРAРР?1x?*1x?*Du?BРР

A* ?РAР?1B* ?РBC* ?Р?1CD* ?D

Общие закономерности преобразования МВ ММ к произвольно выбранному базису

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

1) Преобразование на основе управляемости. Необходимо, чтобы объект был полностью управляемым.

Qу ? [B | AB | A2B |...| AN?1B]

Qу* ? [B* | A*B* | (A*)2B* |...| (A*)N?1B*]

Qу*? [РB | РAР?1РB |...| (РAР?1)N?1РB]

Qу* ?PQу P?Qу*Qу?1

После определения преобразования обязательно необходимо сделать его проверку. Если преобразование из исходного базиса в новый не дает положительного результата, то преобразовать объект исходного базиса в заданный невозможно.

2) Преобразование на основе свойств наблюдаемости

Для его построения необходимо, чтобы объект был полностью наблюдаемым.

QH ? [CT | ATCT |...| (AT )n?1CT ]

Q*H ? [(C*)T | (A*)T (C*)T |...| ((A*)T )n?1(C*)T ]

P?(Q*HT)?1QHT

Наблюдатели (идентификаторы) состояния линейных динамических объектов и систем

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

x?н ? Aнxн ?BнU?FнY

Схема наблюдателя:

В теории управления применяется 2 типа наблюдателей состояний:

1) Наблюдатель состояния Калмана (так называемый полный наблюдатель).

2) Наблюдатель состояния Луенбергера (так называемый редуцированный наблюдатель).

Наиболее распространенным в инженерной практике является ДНС Калмана.

Условия асимптотической оценки наблюдателя состояния: ?1) BH ? B

? ?2) AH ? a ? FHC

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

Модель наблюдателя (Калмана):

x?н ? Axн ? BU ? FH (Y ?Yн)

Для обеспечения асимптотического затухания ошибки наблюдения необходимо, чтобы имела устойчивый спектр корней и чтобы характер затухания удовлетворял условиям асимптотической ошибки необходимо выполнение условия:min Re( pi H )? maxRe( pi )

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

? 0 1 0 ? 1) AOH ? ?? 0 0 1 ??

?? - a0 ?a1 ?a2 ??

? 0 1 0 ? AOHH ? ?? 0 0 1 ??

? - a0H ?a1H ?a2H ??

? ? f1 ?

FH ? ?? f2?? C ?[1 0 0]

?? f3?? ? f1 0 0?

FHC?A?AH COHн ?[1 0 0] Вывод: Основная нормальная форма записи является непригодной для формирования вектора обратных связей ДНС, т.к. невозможно определить все его ( FН ) коэффициенты. COH ? [1 0 0]

FHC ? ?? f2 0 0?? ? AOH ? A0HH

?? f3 0 0??

? 0 0 0?

? ??(aH ? a0 ) (a1H ? a1) (a2H ? a2 )??

0 2) ?0 0 ? a0?

AТОН ? ?1 0 ? a1 ??

? ??0 1 ? a2 ?? CTOH ??0 0 1? ?0 0 ?a0Н ?

?? ATOHH ? ?1 0 ?a1Н ?

?Н ? ??0 1 ?a2 ?? CTOHH ? ?0 0 1? ?0 0 f1 ?

F C ?0 0 f ?A ?0 0 a0H ? a0?

H?0 0 aHa ? A H ? ?2? ? ТOH ? TOH ? ?1 ? 1 ?

??0 0 f3????0 0 a2H ? a2?? fi ? aiH?1 ? ai?1 (i ?1,n)

aiH - коэффициент характеристического полинома наблюдателя.

Преимуществом использования ДНС при построении САУ по сравнению с классическими вход - выходными методами является то, что появляется возможность для формирования закона управления не только непосредственно изменяемых Y и входных воздействий U, но и всех наблюдаемых переменных состояния, оценивающих реальные переменные с точностью определенной выбранными параметрами асимптотической оценки.

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

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

Пример.

? 0 1 0 ? ?0? A ? ?? 0 0 1 ?? В???0??

?? -0,03 ?0,43 ?1,4 ????1??

p1 ??0.1 p2 ??0.3 p3 ??1

?0 0 -0.03?

ATOH ? P1TOH A(P1TOH )?1 ? ??1 0 -0.43??pTOH1

??0 1 -1.4 ??

?1? BTOH ? ??0??CTOH

??0??

С??1 0 0? ?0.43 1.4 1 ?

? ? 1.4 1 0 ?

?? ?? 1 0 0 ??

??0 0 1?

Выберем корни наблюдателя:p1,2,3 ??2

NH(p) ? (p?2)3 ? p3 ?6p2 ?12p?8

?0 0 ?8?a0 ? 0.03 a0H ? 8

ATOHH ? ??1 0 ?12??a1 ? 0.43 a1H ?12

??0 1 ? 6??a2 ?1.4 aH2 ? 6

f1 ? 8? 0.03 ? 7.97 ?7.97 ? f2 ?12 ? 0.43 ?11.57 FH ? ??11.57?? f3 ? 6 ?1.4 ? 4.6 ?? 4.6 ??

x?H ? AxH ?BU ?FH(Y

?1?

?YH)BTOH ? ??0??

??0?? Нелинейные системы

Нелинейной системой называется такая система, в состав которой входит хотя бы одно звено, описываемое нелинейным уравнением. Такое звено называется нелинейным звеном или нелинейным элементом.

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

(НЭ) и линейной части (ЛЧ), как показано на рисунке

Нелинейные звенья классифицируются по различным признакам.

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

Гладкие криволинейные характеристики: а - гистерезисная; б, в - усилительные

Кусочно-линейные характеристики: а - с насыщением;

б - с зоной нечувствительности; в - с насыщением и зоной нечувствительности; г - люфт

Релейные характеристики: а - идеальная; б - с зоной нечувствительности; в - гистерезисная

Особенности нелинейных систем.

Поведение нелинейных систем, при наличии существенных нелинейностей, значительно отличается от поведения их линейных моделей.

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

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

Система устойчива "в малом", если она устойчива при малых

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

3. Для нелинейных систем характерен режим незатухающих периодических колебаний с постоянной амплитудой и частотой

(автоколебаний), возникающий в системах при отсутствии периодических внешних воздействий.

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

Прямой метод Ляпунова

При использовании прямого метода Ляпунова (второй методой Ляпунова), исследуемая система описывается

дифференциальными уравнениями первого порядка, полагая, что они записаны для переходного процесса в отклонениях всех переменных xi (i = 1, 2, ... , n) от их значений в установившемся процессе при новых постоянных значениях возмущающего f = f0 и задающего g = g0 воздействий.

dxdti ? Fi(x1,x2,...,xn) (1)

где Fi - нелинейные функции произвольного вида, удовлетворяющие условию

F1 = F2 = ... = Fn = 0 при x1 = x2 = ... = xn = 0

ТеоремаЛяпуноваобустойчивостинелинейныхсистем

Если при заданных в форме (1) уравнениях системы n-го порядка можно подобрать такую знакоопределенную функцию Ляпунова V(x1, x2, ..., xn ), чтобы ее производная по времени W(x1, x2, ..., xn ), тоже была знакоопределенной (или знакопостоянной), но имела знак, противоположный знаку V, то данная система устойчива; при знакоопределенной функции W будет иметь место асимптотическая устойчивость.

ТеоремаЛяпуноваонеустойчивостинелинейныхсистем Если при заданных в форме (1) уравнениях системы n-го порядка производная по времени W(x1, x2, ..., xn) от какойнибудь функции Ляпунова V(x1, x2, ..., xn) окажется знакоопределенной, причем сама функция V в какой-нибудь области, примыкающей к началу координат, будет иметь знак, одинаковый со знаком производной W, то данная система неустойчива.

Замечания к теореме Ляпунова об устойчивости.

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

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

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

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

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

Пример. Исследовать устойчивость системы, заданной уравне

dx122dx222

??(x1?? x2)(1? ax1? bx2)??(x2?? x1)(1? ax1? bx2) dtdt где ?? ?? а, b - положительные постоянные числа.

Р е ш е н и е. Выбираем положительно-определенную функцию

; Ляпунова

V = ?x12 + ?x22 .

Находим производную от функции Ляпунова по времени

W ? dV ?dx1 2?x dx2 2?x? dt1 dt2 dt

W=?2?x1 (x1 ? ?x2 )(1 ? ax12 ? bx22) ?2?x2 (x2 + ?x1 )(1 ? ax12 ? bx22) = = ?2(1 ? ax12 ? bx22)( ?x12 + ?x22) .

Тогда W 0 или ax12 + bx22 0,

где k - предельное значение коэффициента передачи нелинейного элемента;

WЛЧ (j?) - амплитудно-фазовая частотная характеристика линейной части системы.

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

При наличии одного нулевого полюса требуется еще, чтобы

Im WЛЧ (j?) ? ?? при ? ? 0,

а при двух нулевых полюсах

Re WЛЧ(j?) ? ?? при ? ? 0, а Im WЛЧ(j?) 0.

Метод гармонической линеаризации

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

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

Коэффициенты гармонической линеаризации

В основе метода гармонической линеаризации нелинейностей (рисунок), описываемых уравнением yн = F(x), лежит предположение, что на вход нелинейного элемента подается гармоническое воздействие с частотой ? и амплитудой a, а из всего спектра выходного сигнала выделяется только первая гармоника yн1 = aн1 sin(?t + ?н1) при этом высшие гармоники отбрасываются и устанавливается связь между первой гармоникой выходного сигнала и входным гармоническим воздействием

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

Для нелинейных элементов с характеристикой yн = F(x) в результате разложения периодической функции F(x) в ряд Фурье при синусоидальных колебаниях на входе x = a sin?t получим выражение для первой гармоники сигнала на выходе yн1 = b1F sin ?t + a1F cos ?t

где b1F, a1F - коэффициенты разложения в ряд Фурье, которые определяются по формулам:

1 2?1 2?

a1F ? ?F(asin ?t )cos?t d?t. b1F ?? ?0 F(asin ?t )sin ?t d?t,

? 0 Так как x = a sin?t, то px = a? cos ?t (где p = d/dt).

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

p y ? (q ? q? )x н1?

гдеq = b1F/a, q? = a1F/a.

Уравнениеназываетсяуравнениемгармонической линеаризации,акоэффициентыqиq?-коэффициентами гармонической линеаризации.

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

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

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

WЭ(j?, a) = q + jq? = AЭ(?, a) e j?э(?, a)

где модуль и аргумент эквивалентного комплексного

коэффициента передачи связаны с коэффициентами гармонической линеаризации выражениями

AЭ(?, a) = [q(a,?)]2 ? [q?(a,?)]2 ;

?Э(?, a) = arctg[q?(a, ?)/q(a, ?)].

Исследование симметричных периодических режимов в нелинейных системах.

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

Частотный метод. В соответствии с критерием устойчивости Найквиста незатухающие колебания в линейной системе возникают в том случае, когда амплитудно-фазовая характеристика разомкнутой системы проходит через точку с координатами [?1, j0]. Данное условие является также условием существования автоколебаний в гармонически линеаризованный нелинейной системе, т.е.

Wн(j?, a) = ?1

Wн(j?, a) = Wлч(j?) WЭ(j?, a)=-1 Wлч(j?) =-1/ WЭ(j?, a)

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

Wлч(j?) и годографа обратной характеристики нелинейной части 1/WЭ(j?). Если эти годографы не пересекаются, то режим автоколебаний в исследуемой системе не существует.

Для устойчивости автоколебательного режима с частотой ?0 и амплитудой a0 требуется, чтобы точка на годографе нелинейной части, соответствующая увеличенной амплитуде a0+?a по сравнению со значением в точке пересечения годографов, не охватывалась годографом частотной характеристики линейной части системы и охватывалась точка, соответствующая уменьшенной амплитуде a0??a.

На рисунке дан пример расположения годографов для случая, когда в нелинейной системе существуют устойчивые автоколебания, так как a3 0

центр

2. a12>4a2, a1>0, a2>0

устойчивый узел

x 3. a20 седло

5. a120, a2>0x

особая точка -

устойчивый фокус

6. a12>4a2, a10 неустойчивый узел

Дискретные системы управления

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

К дискретным системам относятся импульсные, релейные и цифровые.

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

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

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

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

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

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

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

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

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

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

Математическое исследование дискретных систем существенно упрощается в том случае, когда все величины рассматриваются в дискретные равноотстоящие моменты времени.

Решетчатая функция времени x[nT], или в сокращенной записи x[n] - это математическая функция, значения которой определены в дискретные равноотстоящие друг от друга моменты времени t = nT, где n - целое положительное число 0, 1, 2 ..., а Т период дискретности.

Z-пребразованием решетчатой функции x[nT] называется функция комплексного аргумента z, определяемая выражением

? -n

X(z) = Z? x[nT] ?? ?x[nT]? z

n = 0 Преобразование, в котором z = epT, называется - "дискретное преобразование Лапласа".

Устойчивость дискретных систем

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

Линейная дискретная система называется устойчивой, если переходная составляющая процесса регулирования yп[n,?] затухает с течением времени.

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

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

Передаточные функции разомкнутых импульсных систем

Разомкнутая линейная импульсная система может быть схематически представлена в виде последовательного соединения импульсного элемента (ИЭ) и непрерывной части (НЧ) ).

Импульсный элемент преобразует задающее воздействие g(t) в последовательность импульсов x*, амплитуда которых пропорциональна входному непрерывному сигналу. Импульсная последовательность после прохождения через непрерывную часть вследствие сглаживающих свойств последней превращается в непрерывную величину на выходе y(t).

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

прямоугольным

Tимп =?Т, W(s) = L w

Частотные характеристики импульсных систем

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

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

получаются из их передаточных функций путем замены оператора z на ej?T. Так как частота ? входит в показатель степени числа e, то частотные характеристики являются периодическими функциями частоты, период изменения которых равен ???T. Следовательно, нельзя различить составляющие, частоты которых кратны частоте квантования импульсного элемента ?0 = 2?/Т.

На практике применяются частотные характеристики относительно абсолютной псевдочастоты ?. Переход к псевдочастоте основан на переходе от z-преобразования к w-преобразованию с помощью подстановки

1+w c последующей заменой комплексной переменной w на

z=абсолютную псевдочастоту w= j?T/2.

1?w При этом реальная частота ? и псевдочастота ? связаны соотношением

2 ?T ?? tg. T2

Описание дискретных систем в пространстве состояний

Современная теория дискретных систем, так же как и

непрерывных, базируется на описании их в пространстве состояний

?X(n+1) = ФX(n)+ Г U(n);

? ? Y(n)= CX(n)+ DU(n),

Ф ? T Ap Bdp eAT, Г ? ?e ?

0 Матрицы Ф и Г определяется разложением матричной экспоненты в ряд ?i

?АТ? Ф ? ?i?0i!, при i ?1 Ф ? E ? AT

Г ? BT

Точность импульсных систем

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

В замкнутой импульсной системе ошибка x, задающее воздействие g и возмущающее воздействие f связаны следующей зависимостью относительно z-изображений

X(z,?) = Xg(z,?) + Xf(z,?)

Выражение содержит z-изображения двух составляющих ошибки:

Xg(z,?) - от задающего и Xf(z,?) - от возмущающего воздействий.

Коэффициенты ошибок. Если задающее воздействие g(t) имеет произвольный вид, предельное значение ошибки вычисляется по

формулеcc

x[n]= c g[n]? c g?[n]? 2 g??[n]?...? k g(k)[n],

012!k! Коэффициентыошибокнаходятподискретнойпередаточной функции замкнутой импульсной системы по ошибке di

c ? i! Ф(z), iixgz?1

dz

Коррекция импульсных систем

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

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

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

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

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

k u[n]? ? w [n?i]x[i],

k i ? 0 где wk[n] - импульсная функция непрерывной части импульсного фильтра.

Наиболее просто импульсные корректирующие устройства реализуются с помощью импульсных RC-цепей.

РазличаюттриструктурыимпульсныхRC-цепей:

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

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

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

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

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

Передаточная функция исходного объекта имеет вид: k

W(p) ? pv(T2? p?1)(T3? p?1)(T45? p?1)

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

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

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

Задачаоценкиработысистемыприизменениипараметров составляющих ее элементов решается путем количественной оценки чувствительностисистемы.Дляэтоготребуетсяописатьсистему управления уравнениями в нормальной форме, т.е. dx i ? ai1x1+ ai2x2+ ... + ainxn + fi(1)

dt где n - порядок системы;

xi - координаты состояния системы;

fi - внешние воздействия, прикладываемое к системе;

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

Изменяющиеся со временем параметры элементов системы в процессе

эксплуатации и от разброса при изготовлении обозначим через ?j (j=1, 2, ... , m). Тогда уравнение системы (1) можно записать в виде

dx i ?? (x ,x , ... ,x ,?,? , ... ,? ,f ) (2) dt i 1 2 n 1 2 m i

Решение уравнений (2) определяет координаты системы: x1(t), x2(t), ... ,

xn(t), образующие исходное движение системы.

Пусть параметры ?j изменяются на малые величины ??j , тогда имеем

?~ ?? ????~ ?? ???

111mmm Рассматривая малые изменения параметров ?j (j=1, 2, ... , m), получим новые уравнения ; d~x

i ?? (~x ,~x , ... ,~x ,? + ?? ,? + ?? , ... ,? + ?? ,f )(3)

dti 1 2n 11 22mmi

Процесс в той же системе, но с измененными параметрами, определяемый решением уравнений (3), т.е.

~x (t),~x (t), ... ,~x(t)называется варьированным движением.

12n Возникшее различие в протекании процессов в системе за счет изменения параметров

?x (t)= x (t)-~x (t) называется дополнительным движением. i i i

При малых отклонениях ??j эта разность может быть определена следующим образом:

? x? x? x

?x (t) =i ???i ???...?i ??

i ? ? 1 ? ? 2 ? ? m

12m ? x Обозначимuij(t)= ? ?i(4) j

Тогда дополнительное движение будет

?xi(t)= ui1 ??1 ? ui 2 ??2 ? ...? uim ??mпри i=1, 2, ... , n.

Величины uij(t) , представляют собой функции чувствительности i-ой координаты системы по j-ому параметру.

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

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

Если показатели качества выражаются не функциями, а числами, то uij называются коэффициентами чувствительности.

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

1 ?W2 ?W4 ?W1

2 . Оценить устойчивость исходной замкнутой системы.

3 . Определить перерегулирование и время регулирования.

4 . Построить ЛАЧХ и ЛФЧХ исходной разомкнутой системы.

5 . Увеличить требования к качеству и реализовать их с помощью последовательной коррекции.

---------------

------------------------------------------------------------

---------------

------------------------------------------------------------

0?, дк0t, с

0?, дк0t, с

0?, дк0t, с

Показать полностью…
Рекомендуемые документы в приложении