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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Конспект лекций

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

Лекция №1. Модели объектов управления

Математические (аналитические) модели широко используются при изучении разнообразных процессов в науке и технике, экономике, в сфере социальных явлений. Известный шведский ученый Леннарт Льюнг пишет в своей книге «Идентификация систем»: «Формирование моделей на основе результатов наблюдений и исследование их свойств - вот, по существу, основное содержание науки. Модели («гипотезы», «законы природы», «парадигмы» и т.п.) могут быть более или менее формализованными, но все обладают той главной особенностью, что связывают наблюдения в некую общую картину». А вот высказывание известного специалиста в области автоматического управления Р. Калмана: «Теория управления не занимается исследованием реального мира, а лишь математическими моделями определенных аспектов реального мира».

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

Необходимо помнить, что объект моделирования, в соответствии с терминологией Канта, есть «вещь в себе» - он не доступен для непосредственного познания. Нам доступны лишь его отдельные фрагменты, полученные с помощью ощущений (измерений), которые можно сложить в единую картину только путем синтеза модели. Модель - это «вещь для нас», с ее помощью мы можем прогнозировать поведение объекта и затем сравнить результаты прогноза с реальной действительностью. Поскольку модель всегда проще объекта, прогноз и реальность никогда идеально не совпадают между собой. Обычно создание модели представляет собой итерационный процесс. Вначале создается простая, но неточная модель. Затем она совершенствуется.

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

Рис. 1.1. Система с одним входом и одним выходом

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

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

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

Любой системе присущи два основных атрибута: целостность и структурированность.

Целостность (единство) означает, что система отделена от внешней среды; среда может оказывать на нее действие (акцию) через входы и воспринимать отклик (реакцию) на эти действия через выходы.

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

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

Функционирование системы - это процесс, разворачивающийся во времени, т. е. множества возможных входов и выходов U, Y - это множества функций времени. Система называется функциональной (определенной), если каждой входной функции u(t) соответствует единственная выходная функция y(t). В противном случае система называется неопределенной. Неопределенность обычно возникает из-за неполноты информации о внешних условиях работы системы. Важным свойством, присушим реальным системам, является причинность. Она означает, что если входные функции и совпадают при , т.е. при , то соответствующие выходные функции удовлетворяют условию , т.е. «настоящее не зависит от будущего при заданном прошлом».

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

Аналитическими моделями различных объектов реального мира являются математические уравнения.

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

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

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

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

Пример 1.1. Построим аналитическую модель электронной схемы, изображенной на рис. 1.2.

Рис. 1.2. -цепь

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

,(1.1)

где - величина сопротивления резистора.

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

,(1.2)

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

Приравнивая правые части уравнений (1.1) и (1.2), получим

.(1.3)

Уравнение (1.3) равносильно линейному дифференциальному уравнению первого порядка

,(1.4)

где , а .

Уравнение (1.4) является аналитической моделью электрической схемы, изображенной на рис. 1.2. Если задать входной сигнал как известную функцию времени, то, решив уравнение (1.4) при заданных начальных условиях, получим выходной сигнал как функцию времени. Другими словами, уравнение (1.4) позволяет нам моделировать реакцию системы на различные входные сигналы. Поэтому оно и называется моделью. Произведение , имеющее размерность времени, является параметром данной системы.

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

,(1.5)

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

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

.(1.6)

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

.(1.7)

График функции, описываемой уравнением (1.7) для случая, когда , показан на рис. 1.3.

График показывает, что в ответ на ступенчатый входной сигнал система переходит в новое состояние не сразу, а спустя некоторое время. Приблизительно оно равно: (точнее, спустя промежуток времени выход системы составляет 63 % от величины выходного сигнала в установившемся режиме).

Рис. 1.3. Переходная функция апериодического звена

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

,(1.8)

где - модуль и - аргумент частотной функции.

Амплитудно-частотная (АЧХ) и фазово-частотная (ФЧХ) характеристики данной системы показаны на рис. 1.4.

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



Рис. 1.4. АЧХ и ФЧХ апериодического звена

Пример 1.2. Построим аналитическую модель электронной схемы, изображенной на рис. 1.5.

Рис. 1.5. -цепь

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

,(1.9)

, (1.10)

где - ток, протекающий через резистор,

- ток, протекающий через конденсатор. На основании (1.9) и (1.10) можем записать

.

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

.(1.11)

В общепринятых обозначениях

.(1.12)

Уравнение в изображениях: - позволит найти передаточную функцию

.(1.13)

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

,(1.14)

где - частота собственных колебаний системы,

- коэффициент демпфирования системы. Сравнивая формулы (1.13) и (1.14), получим:

, .

На рис. 1.6 показаны реакции линейной системы второго порядка на ступенчатый сигнал при различных коэффициентах демпфирования (; ; ). График получен в системе MATLAB.

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

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

Пример 1.3. Построим аналитическую модель механической системы, изображенной на рис. 1.7.



Рис. 1.7. Механическая система перемещения

Предполагаем, что движение возможно только в одном направлении - вдоль оси y(t). Никакое движение в поперечном направлении не допускается. На тело массы М действует три силы: внешняя сила, сила трения и сила упругости. Согласно принципу д'Аламбера, к ним нужно добавить силу инерции и результат приравнять нулю. В результате получим дифференциальное уравнение:

,(1.15)

где В - коэффициент трения,

К - коэффициент упругости.

Это дифференциальное уравнение является моделью механической системы, изображенной на рис. 1.7. Сравнив его с уравнением (1.12), можно заметить, что они отличаются только коэффициентами. Коэффициенту LC соответствует коэффициент М, коэффициенту - коэффициент В, единице соответствует коэффициент К.

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

.(1.16)

Частота собственных колебаний механической системы равна , а коэффициент демпфирования: .

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

Лекция №2. Классификация аналитических моделей

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

Математическая модель системы называется статической, если значение выхода y(t) зависит от значения входа u(t) только в тот же момент времени t. Символически это свойство записывается так:

.(2.1)

Уравнение (2.1) не дифференциальное, а алгебраическое. Мы можем рассматривать алгебраическое уравнение как частный случай дифференциального, в котором все производные неизвестной функции (выходного сигнала системы y(t)) равны нулю. В простейшем случае линейной системы:

, (2.2)

где - коэффициент усиления. Кроме явных функциональных зависимостей, статические модели могут задаваться неявно, в виде уравнения или системы: . Это уравнение должно быть однозначно разрешимо относительно y(t). Статическими моделями пользуются, когда в рамках поставленной задачи инерционностью и «памятью» реальной системы можно пренебречь. Это возможно при выполнении ряда условий:

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

2. входы меняются медленно, т.е. , где - время между изменениями входных воздействий;

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

В динамических, моделях значение y(t) может зависеть от всего прошлого (предыстории) входного процесса: .

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

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

Пример 2.1. Система, моделью которой служит линейное уравнение с постоянными коэффициентами, например следующее:

,(2.3)

- является стационарной. А модель следующего вида:

,(2.4)

- является нестационарной.

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

Лекция №3. Линейные непрерывные модели

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

,(3.1)

где - комплексный параметр преобразования Лапласа, - коэффициенты (вещественные или комплексные).

Иногда в системе имеет место задержка сигнала на величину (так называемое - транспортное запаздывание). Тогда в уравнение (3.1) нужно добавить сомножитель :

(3.2)

Конечное множество чисел: () называется полюсами, а множество () - нулями системы (3.1). Система является устойчивой, если действительная часть всех полюсов отрицательна. В противном случае - система неустойчивая. Комплексные нули и полюса образуют пары комплексно сопряженных чисел.

Пример 3.1. Если система имеет полюс: , то она также имеет и полюс: . Подобных пар может быть несколько. Если система имеет только комплексные полюса (или/и нули), то их количество четно.

Символ допускает различную трактовку. Либо это оператор дифференцирования: , тогда

, (3.3)

где - входной и - выходной сигналы системы.

Либо - это комплексная переменная (параметр преобразования Лапласа), тогда

,(3.4)

где и - преобразования Лапласа входного и выходного - сигналов соответственно.

Числитель и знаменатель передаточной функции (3.1) можно представить в виде полиномов степеней и :

, (3.5а)

, (3.5б)

где - коэффициенты, зависящие от полюсов,

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

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

,(3.6)

где - статический коэффициент усиления,

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

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

,(3.7)

где - статический коэффициент усиления,

- постоянная времени,

- частота собственных колебаний системы,

- коэффициент демпфирования.

В среде программирования «MATLAB» передаточные функции задаются с помощью специальной команды. С помощью «MATLAB» легко решаются многие задачи, связанные с моделированием систем.

Пример 3.4. Предположим, нам необходимо найти нули и полюса передаточной функции

.

В «MATLAB» эта задача решается с помощью следующей программы (тут же приводится и результат):

G=tf([4 7], [5 0 -2 1 3]);

zero(G)

ans =

-1.7500

pole(G)

ans =

0.6992 + 0.6004i

0.6992 - 0.6004i

-0.6992 + 0.4663i

-0.6992 - 0.4663i

Как можно видеть два полюса рассматриваемой системы имеют положительные действительные части. Значит, такая система неустойчива.

Уравнение (3.3) с учетом (3.1) можно записать следующим образом

.

Или по другому

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

(3.8)

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

(3.9)

Решение уравнений (3.8) и (3.9) имеет следующий вид

,

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

Пример 3.5. Дифференциальное уравнение линейной системы первого порядка имеет вид:

. (3.10)

Начальные условия предполагаются нулевыми: .

Пример 3.6. Дифференциальное уравнение линейной системы второго порядка имеет вид:

, (3.11)

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

Уравнение (3.6) можно переписать так

.(3.12)

Оригиналом уравнения (3.10) является свертка

,(3.13)

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

Пример 3.7. Импульсная весовая функция линейной системы первого порядка имеет вид:

.(3.14)

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

.(3.15)

Пример 3.8. Импульсная весовая функция линейной системы второго порядка имеет вид:

.(3.16)

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

(3.17)

В среде программирования «MATLAB» реакция системы с передаточной функцией на ступенчатый входной сигнал рассчитывается с помощью команды step(G).

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

.

В командном окне наберем программу:

G=tf([1], [1 0.1 1]);

step(G)

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

В результате получим график, показанный на рис. 3.1.

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

Поскольку в уравнениях (3.6) и (3.12) - понимается как комплексная переменная, то можем записать

,

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

,(3.18)

где , .

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

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

,(3.19а)

.(3.19б)

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

,(3.20а)

.(3.20б)

В «MATLAB» диаграмму Боде системы с передаточной функцией G(s) рассчитывают с помощью команды bode(G).

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

В результате получим график, показанный на рис. 3.2.

Рис. 3.2. Диаграмма Боде системы второго порядка

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

Заменив в формуле (2.12) на , получим следующее выражение

, (3.21)

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

Можно записать и так

,(3.22)

где , , , - амплитудные и фазовые спектры входного и выходного сигналов. Из выражения (3.22) следует

,(3.23а)

.(3.23б)

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

,.

С учетом (3.23а) можем записать

.(3.24)

Частотная характеристика является одной из возможных форм описания линейной системы (3.1). Она однозначно определяет все ее свойства.

Лекция №4. Линейные модели в пространстве состояний

В середине прошлого столетия в широкое употребление вошла еще одна форма линейной модели: модель в пространстве состояний (State Space model). Такая модель описывает связь между входными и выходными сигналами в виде системы дифференциальных уравнений первого порядка посредством введения вспомогательного вектора состояния .

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

Такую модель достаточно просто перевести в модель в пространстве состояний. Рассмотрим уравнение (3.8).

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

.(4.1)

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

(4.2а)

.(4.2б)

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

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

,(4.3а)

.(4.3б)

Здесь - квадратная матрица размером (); - матрица размером (), - размерность вектора входных сигналов; - матрица размером (), - размерность вектора выходных сигналов; - матрица размером ().

Выражения (4.3а) и (4.3б) являются математической моделью линейной системы в пространстве состояний. Уравнение (4.3а) называется уравнением состояния, а уравнение (4.3б) - уравнением наблюдения.

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

,,,.

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

,,,.

Корни уравнения

(4.4)

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

,(4.5)

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

Пример 4.3. Для системы первого порядка единственное собственное число равно (величина, обратная времени, имеет физический смысл частоты). Для системы второго порядка , где - частота собственных колебаний системы, - коэффициент демпфирования, собственные числа этой матрицы равны полюсам системы

,

.

Собственные векторы равны

,.

В «MATLAB» предусмотрена возможность определения передаточной функции по модели, заданной в пространстве состояний.

Пример 4.4. Задана модель в пространстве состояний, имеющая следующие матрицы:

,,, .

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

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

.(4.6)

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

A=[-3 -0.5; 4 0];

B=[0.5; 0];

C=[0 0.5];

h=ss(A, B, C, 0);

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



То же самое можно представить в матричной форме

,(4.7)

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

,(4.8)

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

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

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

.

Применив преобразование подобия (4.8), получим:

.

Решая это уравнение относительно , получим модель системы с вектором состояния .

.

Уравнение для выходной переменной

,

с учетом преобразования подобия (4.8) принимает вид

.

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

,(4.9а)

,(4.9б)

где , , , .

Пример 4.5. Задана система в пространстве состояний со следующими матрицами:

,,.

Необходимо определить матрицы , если матрица преобразования базиса равна:

.

Эту задачу можно решить в системе «MATLAB» следующим образом. Соответствующая программа и результат ее выполнения имеют вид:

Av=inv(P)*A*P

Bv=inv(P)*B

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

Лекция №5. Дискретные модели

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



Рис. 5.1. Цифровая система управления

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

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

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

. (5.1)

Порядок уравнения определяет «глубину памяти» системы.

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

Решить разностное уравнение означает найти последовательность . Такую последовательность называют решетчатой функцией. Существует три основных метода решения линейных разностных уравнений с постоянными коэффициентами. Первый (классический) метод состоит в нахождении общего и частного решений подобно тому, как это делается при классическом решении линейных дифференциальных уравнений. Этот метод мы рассматривать не будем. Второй метод является рекуррентным; он используется при решении разностных уравнений с помощью цифрового компьютера. Мы рассмотрим его на примере.

Пример 5.1. Получим решение следующего разностного уравнения

Причем, . Решения для можно получить, положив сначала в разностном уравнении , затем , затем и т.д. В результате получим

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

mkminus1=0; ekminus1=0; ek=1;

for k=0:5

mk=ek-ekminus1-mkminus1;

[k ek mk]

mkminus1=mk;

ekminus1=ek;

ek=1-ek;

end

В этой программе ekminus1 соответствует значению , ek - значению , mkminus1 - значению , а mk - значению .

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

.

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

.

Тогда дифференциальное уравнение приближенно примет вид:

.

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

.

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

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

.(5.2)

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

(5.3)

где - параметр -преобразования, - параметр преобразования Лапласа,

- шаг дискретизации,

- изображение входного сигнала,

- изображение выходного сигнала.

-преобразование основано на теореме операционного исчисления о запаздывании. Если , то .

Уравнение (5.3) можно переписать следующим образом

.(5.4)

Поскольку известно, то можно найти, применив обратное -преобразование к выражению (5.4).

Пример 5.2. Рассмотрим разностное уравнение из предыдущего примера

Найдем -преобразование этого уравнения

.

Отсюда следует

.

Изображение входного сигнала можно представить в виде

.

Решетчатая функция равна коэффициентам полученного ряда

.

В программе имитационного моделирования «Simulink», которая является частью языка технического программирования «MATLAB», модель дискретной системы задается в виде рациональной передаточной функции

,(5.4)

где - коэффициенты (вещественные или комплексные).

Конечное множество чисел: () называется полюсами, а множество () - нулями системы (5.4). Полюса (и нули) могут быть действительными, либо комплексными. В последнем случае они образуют пару комплексно-сопряженных чисел. Если система устойчивая, то модули всех ее полюсов меньше единицы. В противном случае - система неустойчивая.

Пример 5.3. Дискретная система первого порядка (инерционное звено) имеет передаточную функцию

,(5.5)

где и - коэффициенты ( - полюс системы).

Пример 5.4. Дискретная система второго порядка имеет передаточную функцию

,(5.6)

где и - полюса системы, , .

Пример 5.5. Построим в «MATLAB» модель дискретной системы второго порядка, показанной на рис. 5.2. На рис. 5.3 приведена реакция этой системы на ступенчатый входной сигнал.

Рис. 5.2. Устойчивая дискретная система второго порядка

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

Пример 5.6. Построим в «MATLAB» модель дискретной системы второго порядка, показанной на рис. 5.4. На рис. 5.5 приведена реакция этой системы на ступенчатый входной сигнал.



Рис. 5.4. Неустойчивая дискретная система второго порядка

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

Дискретная система, также как и непрерывная, может быть представлена в пространстве состояний:

- уравнение состояния;

- уравнение наблюдения, где

· - входной сигнал;

· - выходной сигнал;

· - вектор состояний;

· A, B, C, D - параметрические матрицы.

Пример 5.7. Система первого порядка может быть описана такими параметрами:

(5.4)

Пример 5.8. Система второго порядка может иметь следующие матрицы:

;;; (5.5)

Лекция №6. Нелинейные модели

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

Пример 6.1. Уравнение движения вертикально стартующей ракеты записывается в виде

,(6.1)

где - высота ракеты над точкой старта,

- мгновенная масса ракеты,

- ускорение свободного падения,

- положительные коэффициенты.

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

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

,(6.2а)

,(6.2б)

где - -мерный вектор состояния,

- -мерный вектор управления,

- -мерный вектор выходных переменных,

и - нелинейные вектор-функции. Зависимость этих функций от отражает действие возмущений. Как и в системе (4.2) уравнение (6.2а) называется уравнением состояния, а (6.2б) - уравнением наблюдения.

Пример 6.2. Преобразуем уравнение (6.1) в модель в пространстве состояний. Для этого введем замену переменных: , , , . В результате получим

(6.3а)

.(6.3б)

Здесь (6.3а) - векторное уравнение состояния, а (6.3б) - уравнение наблюдения. Как можно видеть, не все уравнения этой системы нелинейные.

Динамические нелинейные системы с гладкими нелинейностями описываются с помощью рядов Вольтерра (Вито Volterra (3.5.1860, Анкона, - 11.10.1940, Рим) итальянский математик, профессор университетов в Пизе, Турине, Риме, член национальной академии деи Линчеи в Риме). Построение рядов Вольтерра заключается в следующем. Представим нелинейную систему в виде двух последовательно соединенных блоков - линейного и нелинейного (рис. 6.1).

Рис. 6.1. Представление нелинейной системы в виде последовательного соединения линейного и нелинейного звеньев

Выходной сигнал линейного звена может быть представлен в виде интеграла свертки

.(6.4)

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

.(6.5)

Подставляя (6.4) в (6.5), приходим к выражению

,(6.6)

которое и определяет ряд Вольтерра. В уравнении (6.6) используются обозначения: ; и т.д.

Пример 6.3. Преобразуем уравнение (6.1) в ряд Вольтерра. Для простоты перейдем от неизвестной функции к функции , имеющей смысл скорости движения ракеты.

.(6.7)

Определив функцию , мы можем найти и функцию по формуле

.(6.8)

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

(6.9)

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

.(6.10)

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

.(6.11)

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

(6.12)

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

Пример 6.4. Преобразуем уравнение (6.7) в конечно-разностное уравнение. Для этого представим производную в виде разности:

,(6.13)

где - шаг дискретизации. Получим конечно-разностное уравнение

.(6.14)

Поскольку , можем записать:

.(6.15)

Используя начальное условие и задав закон изменения входного управляющего сигнала , получим решение конечно-разностного уравнения (6.14) путем организации итерационного процесса:

На рис. 6.2. показаны результаты численного решения уравнения (6.7) с использованием конечно-разностных уравнений (6.14) и (6.15). Как можно видеть, реакция системы на ступенчатое входное воздействие качественно подобна реакции хорошо демпфированной линейной системы второго порядка. Сразу после запуска ускорение ракеты растет. Однако при увеличении скорости возрастает также сила сопротивления воздуха, поэтому при постоянном управляющем воздействии ускорение достигает максимального значения и затем начинает уменьшаться. Как известно, для выхода на околоземную орбиту ракете необходимо достичь первой космической скорости. Оптимальная стратегия управления должна заключаться в том, чтобы достичь этой скорости без превышения определенного уровня ускорения, которое может губительно сказаться на здоровье космонавтов. Расход топлива при этом должен быть минимальным.

Рис. 6.2. Результаты численного решения нелинейного уравнения (6.7) при ступенчатом входном воздействии

Необходимо отметить, что модель (6.7) не учитывает падения давления воздуха с набором высоты. Поэтому она верно описывает движение ракеты только в приповерхностном слое атмосферы.

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

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

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

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

Рис. 6.3. Нелинейные звенья

На рис. 6.3 приведены примеры типовых статических нелинейных звеньев, используемые в программе «Simulink» языка технических вычислений «MatLab». Coulomb & Viscour Friction (Колебания Коломбо-Вискура) - расщепляет сигнал на ряд сигналов, колеблющихся в определенном диапазоне; Dead Zone (Мертвая зона) - заменяет значение входного сигнала, лежащее в заданном диапазоне, нулем; Hit Crossing (Обнаружено пересечение) - позволяет идентифицировать момент времени, когда входной сигнал «пересекает» некоторое значение: при появлении такой ситуации на выходе блока формируется единичный импульсный сигнал; Quantizer (Квантование сигнала) - преобразует непрерывный сигнал в дискретный; Rate Limiter (Ограничитель уровня) - имеет переходный процесс в виде прямой линии (если амплитуда синусоидального сигнала превысит определенное значение, сигнал на выходе приобретает треугольную форму); Relay (Реле) - имитирует работу реле с гистерезисом, если входной сигнал превысит некоторое пороговое значение, то на выходе блока формируется «разрешающий» сигнал; Saturation (Насыщение) - выходной сигнал не может превысить заданный уровень.

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

Этот метод применим к системам, которые содержат только одну нелинейность. Система такой конфигурации изображена на рис. 6.4.

Рис. 6.4. Вариант нелинейной системы

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

.(6.16)

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

.(6.17)

Коэффициенты ряда Фурье определяются выражениями

,(6.18)

,(6.19)

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

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

.(6.20)

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

.(6.21)

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

;.

Теперь можем записать

.(6.22)

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

.(6.23)

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

.(6.24)

Отношение изображений имитирует «передаточную функцию» нелинейного звена.

.(6.25)

Имитация «частотной функции» имеет следующий вид:

.(6.26)

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

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

Лекция №7. Цель управления. Методы синтеза регуляторов

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

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

Известный специалист в области теории автоматического управления Р. Калман дает следующее определение цели управления: «Цель управления состоит в том, чтобы изменить динамику поведения физической системы в соответствии с желаниями человека».

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

Регулятор проектируют под конкретный объект. Чтобы спроектировать регулятор, необходимо иметь модель объекта (например, в виде передаточной функции). Известны три способа синтеза регулятора: 1) синтез с помощью корневого годографа; 2) синтез в частотной области; 3) синтез с помощью современной теории управления (метод назначения полюсов).

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

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

Рис. 7.1. Замкнутая система управления

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

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

Для изучения метода корневого годографа будем рассматривать систему общего вида, изображенную на рис. 7.2, где - передаточная функция объекта управления, H(s) - передаточная функция обратной связи, - передаточная функция регулятора.

Рис. 7.2. Система для изучения метода корневого годографа

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

.(7.1)

Характеристическое уравнение соответственно равно

.(7.2)

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

В общем случае передаточная функция разомкнутой системы равна: . Корневой годограф можно построить с помощью программы MATLAB. Для этого нужно набрать команду: rlocus(F(s)), где («root locus» - по английски означает «корневой годограф»).

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

.

Решение. Запустив программу MATLAB, набираем следующие команды:

F=tf([0.2], [1 1 0.4 0.25]);

rlocus(F)

Результат показан на рис. 7.3.



Рис. 7.3. Корневой годограф для системы из примера

Из рисунка видно, что характеристическое уравнение замкнутой системы имеет три корня: один вещественный и два комплексно-сопряженных. Помеченные крестиками точки соответствую значению параметра K=0. Как можно видеть, при увеличении параметра комплексные полюса пересекают мнимую ось, так что действительные части их становятся положительными, а замкнутая система - неустойчивой. Таким образом, чтобы система оставалась устойчивой, значение K должно быть небольшим.

Заметим, что корневой годограф всегда симметричен относительно вещественной оси.

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

.(7.3)

Система (7.3) является системой второго порядка с частотой собственных колебаний и коэффициентом демпфирования, равным нулю. Ее полюса равны: . Поскольку , то это чисто мнимые полюса. Корневой годограф замкнутой системы показан на рис. 7.4.

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

Рис. 7.4. Корневой годограф системы с П-регулятором

Чтобы получить работоспособную систему, воспользуемся ПД-регулятором: . Обычно говорят, что в такой системе использована обратная связь по скорости, или по производной. Тогда разомкнутая система имеет вид