63931 (Анализ систем автоматического управления)

Министерство образования Республики Беларусь

Белорусский Государственный Университет Информатики и Радиоэлектроники

Факультет Информационных Технологий и Управления

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

Курсовая работа

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

Тема: Анализ систем автоматического управления

Выполнил:Проверил:

Студент гр. 722401Кузнецов А. П.

Царик А.С.

Минск 2010

Содержание

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

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

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

Литература

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

Задание:

1. Найти передаточную функцию разомкнутой системы W(s) и передаточную функцию замкнутой системы Ф(s), ;

2. Построить область устойчивости системы в плоскости общего коэффициента передачи К = К₁К₂К₃ и постоянной времени Т₂ при заданных значения Т₁ и Т₃. Найти граничное значение при заданном значении Т₂, при котором система выходит на границу устойчивости.

3. Построить графики логарифмических амплитудной и фазовой частотных характеристик L(w) и φ(w) при значении коэффициента передачи K=0,7K’.

4. Оценить запасы устойчивости по модулю ∆L и фазе ∆ φ, величину ошибки по скорости е_ск при v(t) = v₁t и f= 0, время переходного процесса t_p и перерегулирование σ в исходной системе при K=0,7K’.

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

6. Вычислить в скорректированной системе переходный процесс на выходе y(t) при подаче на вход единичной ступенчатой функции v(t)=1(t)( f= 0). Найти t_p, σ по переходному процессу и сравнить их с требуемым по заданию.

Исходные данные:

Структура исследуемой замкнутой линейной непрерывной САУ представлена на рис. 1.1, где v(t)- управляющее воздействие, (f)- возмущающее воздействие, е(t)- сигнал ошибки, y(t)- выходной сигнал. Значения параметров Т₁ Т₂, Т₃ заданы в табл. 1. Размерность Т₁ Т₂, Т₃ в секундах, общий коэффициент передачи К = К₁К₂К₃ имеет размерность 1/с, в табл. 1 заданы также желаемые показатели качества системы: максимальная ошибка по скорости е_ск при скачке по скорости v(t) = v₁t и f= 0, время переходного процесса t_п.п в секундах, и перерегулирование у в процентах.

Таблица 1. Структура исследуемой замкнутой линейной непрерывной САУ

Рисунок 1.1

Выполнение:

1. Требуемые передаточные функции находят с использованием правил структурных преобразований. Коротко сформулируем основные правила:

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

• Передаточные функции параллельно соединенных звеньев складываются.

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

(по условию)

Передаточная функция разомкнутой системы W(s) = Y(s)/U(s) при f= 0, e = u (т.е. разомкнута главная обратная связь) определится выражением:

где обозначим К = К₁К₂К₃,

0,03135

1,12127

5,223

Главная передаточная функция или передаточная функция замкнутой системы при f = 0:

Передаточная функция по ошибке при f= 0, которая позволяет выразить ошибку e(t) в системе при известном входном воздействии:

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

2. . Передаточная функция разомкнутой исходной системы имеет вид W(s) = K/sL(s), где L(s) = (T₁s+1)(T₂s +1)(T₃s+1). Характеристическое уравнение замкнутой системы будет D(s) = K+L(s)s = b₀s⁴ +b_}s³ +b₂s² +b₃s + b₄ =0, где при заданных из таблицы исходных данных числовых значениях Т₁ и Т₃ коэффициенты b_j будут зависеть от параметров К и Т₂. Применение критерия Гурвица к характеристическому уравнению четвертого порядка дает следующие условия устойчивости:

b₃(b₁b₂-b₀b₃)-b₄b₁² > 0, b, > 0, i = 0,...,4.

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

К = К₁К₂К₃

b₀= =0,165 =с₀

5,033 с₀

b₃=1 b₄=K

Выразим К через параметр Т₂:

Зависимость К(Т₂) приведена на рис. 1.2

Рис.1.2

K_гр=K_T2=0.19=4,633

3. Полагая К = 0.7К_ГР, записываем аналитическое выражение для φ(w)= argW(jw), L(w) = 20lg|W(jw)| из W(s) при s = jw.

К=0.7К_гр= 3,243

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

где

тогда:

где

Строим графики логарифмических характеристик разомкнутой системы, с помощью MATLAB (оператор bode или margin) Рис. 1.3 а.

Рис. 1.3а

Строим график АФЧХ с помощью MATLAB (оператор nyquist) рис. 1.3 б для разомкнутой системы.

Рис 1.3 б

Запасы устойчивости по модулю и фазе определяются по логарифмическим характеристикам (см. рис. 1.3 а): на частоте среза w_с определяется запас по фазе —∆φ, а запас по амплитуде ∆L - на частоте при которой φ(w) = -180. Таким образом, ∆L≈0. 1дБ, ∆φ≈ 0°, что является недостаточным.

4. Величина ошибки по скорости определяется как e_ск=V₁/K. Для ориентировочной оценки t_пп и σ следует построить переходной процесс h(t) (оператор step в MATLAB) при v(t) = 1[t] и по нему определить t_пп и σ.

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

D(s)y(t)=Kv(t). Если D(s)=b₀s⁴+b₁s³+b₂s²+b₃s+b₄=0, ,то уравнение состояния имеет вид

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

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

где: х - вектор состояния; v- вектор входа; у - вектор выхода.

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

sys=ss(A,B,C,D)

В результате под именем sys получаем ss-объект с числовыми характеристиками в виде четверки матриц {А, В, С, D}, которые должны иметь согласованные размеры. Матрицу D в данном случае полагаем равной 0.

Для построения переходного процесса h(t) воспользуемся оператором step в MATLAB.

Реализация функций имеет вид:

sys=ss([0 1 0 0;0 0 1 0;0 0 0 1;-b4/b0 -b3/b0 -b2/b0 -b1/b0],[0 0 0 K/b0]', eye(4), zeros(4,1))

a =

x1 x2 x3 x4

x1 0 1 0 0

x2 0 0 1 0

x3 0 0 0 1

x4 -104.6 -32.26 -168.5 -36.16

b =

u1

x1 0

x2 0

x3 0

x4 104.6

c =

x1 x2 x3 x4

y1 1 0 0 0

y2 0 1 0 0

y3 0 0 1 0

y4 0 0 0 1

d =

u1

y1 0

y2 0

y3 0

y4 0

Continuous-time model.

>> step(sys)

В результате получим графики представленные на рис. 1.4. Нас будетинтересовать Out(l). Величина ошибки по скорости определяется как:

е_ск=V₁/K = 1,4/3,243 = 0,432>е_скзад = 0,04.

Для ориентировочной оценки t_nn и о следует построить переходной процесс h{t) (оператор step в MATLAB) при v(t)=1(t) и по нему определить t_пп и σ. Эти величины из графика Out(l) определяются следующим образом:

Время переходного процесса определяется с учетом следующих соотношений: ε_уст=v(t)/(l+K), где v(t)=l[t], а К=3,243 - общий коэффициент передачи разомкнутой системы. Тогда е_уст= 1/(1+3,243)=0,236 и следовательно t_пп из графика Out(l) t_пп ≈50с > t_ппзад = 2.5с.

Рис 1.4

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

5. Если исходная система не удовлетворяет заданным показателям качества, ее следует скорректировать. В случае применения частотных методов синтеза коррекции строится желаемая ЛАЧХ L_ж(w). В низкочастотной части желаемой ЛАЧХ при сохранении порядка астатизма (наличие интегратора 1/s в системе) требуемый коэффициент усиления выбирается из соотношения Kz=v₁/e_ск=1,4 / 0.04 = 35. На частоте среза желательно иметь наклон ЛАЧХ -20 дБ/дек с протяженностью этого участка не менее одной декады. Далее среднечастотная часть ЛАЧХ сопрягается с низкочастотной отрезком прямой с наклоном -40(если необходимо -60) дБ/дек, а высокочастотная часть желаемой и исходной ЛАЧХ по возможности должны совпадать.

Учет требований качества переходного процесса: t_пп и σ, запасов устойчивости учитываются при формировании среднечастотной области L_ж(w). Здесь можно воспользоваться графиком (рис. 1.5).

Рис 1.5

По графику рис. 1.5 для заданных значений у и t_nn находим w_п, и затем из соотношения w_c = (0.6 - 0.9) w_п, частоту среза w_c.