Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

лабораторная работа / metodicheskie_ukazaniya_po_vypolneniyu_laboratornyh_rabot_po

.pdf
Скачиваний:
34
Добавлен:
22.02.2014
Размер:
663.11 Кб
Скачать

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ

Государственное образовательное учреждение высшего профессионального образования

САНКТ)ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ АЭРОКОСМИЧЕСКОГО ПРИБОРОСТРОЕНИЯ

ТЕОРИЯ

АВТОМАТИЧЕСКОГО

УПРАВЛЕНИЯ

Методические указания

квыполнению лабораторных работ

1–9

Санкт)Петербург

2006

Составители: М. В. Бураков, Т. Г. Полякова, А. В. Подзорова Рецензент С. А. Гусев

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

Предназначены для студентов специальности 220201 «Управле) ние и информатика в технических системах», изучающих дисцип) лину «Теория автоматического управления».

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

Редактор А. М. Смирнова Верстальщик C. Б. Мацапура

Подписано к печати 11.10.05. Формат 60×84 1/16. Бумага офсетная. Печать офсетная. Уч. )изд. л. 3,2. Тираж 100 экз. Заказ №

Редакционно)издательский отдел Отдел электронных публикаций и библиографии

Отдел оперативной полиграфии ГУАП

190000, Санкт)Петербург, ул. Б. Морская, 67

© ГОУ ВПО «СПбГУАП», 2006

2

ПРЕДИСЛОВИЕ

Настоящие указания составляют учебно)методическую базу для выполнения лабораторных работ по курсу «Теория автоматического управления» (ТАУ) студентами направления «Автоматизация и уп) равление».

Цель работы: развить и закрепить у студентов навыки практичес) кого анализа и проектирования систем управления.

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

В процессе выполнения лабораторных работ студенты должны исследовать такие вопросы, как:

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

исследование устойчивости замкнутых систем;

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

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

Выполнение лабораторных работ предполагает использование популярного мощного пакета моделирования MatLab с расширения) ми Control System Toolbox и Simulink. Каждая работа содержит опи) сание команд MatLab, которые могут быть использованы при выпол) нении заданий.

При подготовке и выполнении лабораторных работ рекомендует) ся использовать классические учебные пособия по теории автомати) ческого управления [1, 2], а также [3–6].

Для углубленного знакомства с возможностями системы MatLab могут быть рекомендованы учебные пособия [7–9].

3

Лабораторная работа №1

ДИНАМИЧЕСКИЕ ЗВЕНЬЯ И ИХ ХАРАКТЕРИСТИКИ ВО ВРЕМЕННОЙ ОБЛАСТИ

1. Методические указания

Рассмотрим систему автоматического управления (САУ), описы) ваемую линейным дифференциальным уравнением вида:

a

dny(t)

+ a

dn−1y(t)

+ ...+ a

dy(t)

+ a y(t) =

 

 

dtn

dtn−1

 

 

 

n

 

n−1

 

1 dt

 

0

 

= bm

 

dmu(t)

+ bm−1

dm−1u(t)

+ ...b1

du(t)

+ b0u(t),

(1)

 

dtm

dtm−1

dt

 

где u(t) – входной процесс; y(t) – выходной процесс; ai, bj, – постоян) ные коэффициенты; n, m (n >= m) – постоянные числа.

Если ввести обозначение p для оператора дифференцирования p = d , то можно записать (1) в операторной форме:

dt

 

 

 

 

 

 

(a pn + a

pn−1 + ...+ a p + a )y(t) =

 

n

n−1

1

0

 

= (bm pm + bm−1 pm−1 + ...+ b p + b )u(t),

(2)

 

 

1

0

откуда получается:

 

 

 

 

 

 

 

y(t)

=

B( p)

= W( p),

 

 

 

u(t)

 

 

 

 

 

A( p)

 

 

где A(p) и B(p) – полиномы из формулы (2).

Выражение (2) по виду совпадает с определением передаточной функции (ПФ) как отношения преобразования по Лапласу выход) ной переменной к преобразованию по Лапласу входной переменной при нулевых начальных условиях:

y(s)

=

B(s)

= W(s) ,

(3)

u(s)

 

 

A(s)

 

где s – комплексная переменная.

Комплексные числа, являющиеся корнями многочлена B(s), на) зываются нулями передаточной функции, а корни многочлена A(s) –

полюсами.

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

4

Таблица

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

¹

 

Название звена

 

 

 

 

 

ПФ звена

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Интегрирующее

 

 

 

 

 

W(s) =

K

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

Дифференцирующее

 

 

 

 

 

W(s) = Ks

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

Усилительное

 

 

 

 

 

W(s) = K

 

(безынерционное)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

Апериодическое 1)го

 

 

 

 

 

W(s) =

 

 

 

K

 

порядка (инерционное)

 

 

 

 

 

Ts

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Апериодическое 2)го

 

 

 

W(s) =

 

K

 

 

 

 

 

 

 

; T1 2T2

5

 

порядка (все корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2s2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вещественные)

 

 

 

 

 

+ T s + 1

 

 

 

 

 

2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

Kолебательное*

 

 

 

W(s) =

 

 

K

 

 

 

 

 

 

 

; T1 < 2T2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2s2 + T s + 1

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

Kонсервативное

 

 

 

 

W(s) =

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ts2

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрирующее

 

 

 

 

 

 

 

 

 

 

 

 

 

K

8

 

с запаздыванием

 

 

 

W(s) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(реальное

 

 

 

 

s(Ts+1)

 

 

 

 

 

 

 

 

 

 

 

 

интегрирующее)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцирующее

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

с запаздыванием

 

 

 

 

 

W(s) =

Ks

 

 

 

 

(реальное

 

 

 

 

 

 

 

 

 

 

 

Ts + 1

 

 

дифференцирующее)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

Форсирующее

 

 

 

W(s) = K(Ts+1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

Изодромное

 

 

 

W(s) =

K(Ts +1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* Часто и с п ользуется описание

ко лебательного звена в виде:

 

 

 

 

 

 

 

 

 

 

 

W(s) =

 

K

 

;

T =T ,

 

ξ =

 

T1

 

 

 

 

 

 

 

 

 

2T2 .

 

 

 

T2s2

+2ξTs +1

2

 

 

 

 

 

 

 

 

5

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

Единичный скачок 1(t) определяется условиями:

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

1 при t >0.

Реакция САУ на единичный скачок называется переходной функ цией системы и обозначается h(t). При неединичном ступенчатом воз) действии g(t)=N1(t), где N = const, в соответствии с принципом супер) позиции выходная реакция системы будет

y(t) = Nh(t) .

Импульсная функция δ(t) определяется условиями:

∞ при t =0,

δ(t) = ≠0 при t 0.

Очевидно:

δ(t) =1′(t).

Реакция САУ на импульсную функцию называется импульсной переходной функцией системы (функцией веса) и обозначается w(t).

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

t

h(t) = w(τ) dτ.

0

2.Использование пакета MatLab

Впакете MatLab имеется два основных варианта для исследова) ния передаточных функций и моделирования САУ:

– использование команд пакета расширения Control System Toolbox;

– использование пакета Simulink.

Control System Toolbox [8, 9] предназначен для работы с LTI)моде) лями (Linear Time Invariant Models – линейные модели с постоянны) ми параметрами) систем управления.

Команда, создающая LTI)систему c одним входом и одним выхо) дом в виде передаточной функции, имеет следующий синтаксис:

TF([bm, …, b1, b0], [an, …, a1, a0]),

6

где bm, …, b1, b0 и an, …, a1, a0 значения коэффициентов полиномов В и A в (3).

Например, если требуется описать ПФ вида

W =

 

s + 1

 

 

2s2

+ 8s + 5

 

и узнать значения ее нулей и полюсов, то нужно ввести в окне команд MatLab следующие команды:

>>w=tf([1 1],[2 8 5])

>>zero(w)

>>pole(w)

Исследовать реакцию LTI)модели на типовые входные воздействия можно с помощью команд

>>step(w)

>>impulse(w)

Можно получить на одном графике реакцию сразу нескольких ди) намических звеньев, если использовать команды вида:

>>step(w,w1,w2)

>>impulse(w, w1,w2)

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

>>step(w, w1, w2, t),

где t – время моделирования в секундах.

На рис. 1 показан пример моделирования динамики колебатель) ного звена при различных параметрах:

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

379

 

 

 

 

 

 

 

 

 

 

 

378

 

 

 

 

 

 

 

 

 

 

9

376

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

6

374

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

179

 

 

 

 

 

 

 

 

 

 

 

178

 

 

 

 

 

 

 

 

 

 

 

176

 

 

 

 

 

 

 

 

 

 

 

174

 

 

 

 

 

 

 

 

 

 

 

1 1

2

31

32

41

42

51

52

61

62

21

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1. Исследование реакции колебательного звена

7

>>w=tf([1],[2 0.3 1]);

>>w1=tf([1],[2 0.5 1]);

>>w2=tf([1],[2 0.1 1]);

>>step(w,w1,w2,50).

ВSimulink MatLab ПФ можно описать с помощью блока Transfer fcn в разделе библиотеки Continuous. Для подачи типовых воздей) ствий надо использовать блок Step из раздела Sources. Импульсную переходную характеристику звена можно получить, подавая на вход

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

3.Задание на лабораторную работу

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

Отчет по лабораторной работе должен содержать:

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

– экспериментально полученные характеристики при вариации параметров каждого звена;

– выводы, обобщающие проделанные эксперименты по каждому звену.

8

Лабораторная работа №2

ИССЛЕДОВАНИЕ ПРЕОБРАЗОВАНИЙ СТРУКТУРНЫХ СХЕМ

1. Методические указания

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

Структурной схемой называется схема САУ, изображенная в виде соединения ПФ составляющих ее звеньев.

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

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

Структурная схема может быть составлена по уравнению системы

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

Различные способы преобразования структурных схем облегчают определение ПФ сложных САУ и дают возможность привести много) контурную систему к эквивалентной ей одноконтурной схеме.

Преобразование структурной схемы должно осуществляться на основании правил. Правила преобразования структурных схем мож) но найти в справочной литературе [1, 2], основные из них приведены

втабл. 1.

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

9

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1

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

 

 

 

 

 

123452674869 3

 

 

 

 

2 296 3 6

 

 

 

 

 

 

 

 

 

 

4 96

 

 

 

8 86 39 96

 

 

 

 

 

 

 

 

 

 

 

 

 

12345627893

11

12

 

 

 

11

 

1

 

 

1

 

 

 

 

 

3 2753 8

 

 

 

 

 

2

 

3

 

 

 

 

2

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

3 98389

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

141112611

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

111

 

 

 

 

 

1

 

 

1

1

 

 

 

 

12345627893

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

747 3 8

21

 

12

 

 

 

31

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 98389

 

111

 

 

 

141151256511

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

1

 

1

1

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

12345627893

 

211

 

 

 

 

3

 

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4758 2 9

 

 

 

 

 

 

 

 

 

 

W1

 

 

 

 

 

 

1

2

 

 

 

 

W =

 

 

 

 

 

 

 

 

 

 

 

± W1W2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

21

11

 

 

3 1

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

3438 7 343

2

1 1

 

3

1

 

 

 

 

 

 

111

71

 

 

238 2 343

 

71

 

 

 

 

 

= 1

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

W

 

 

 

 

 

 

 

 

11

 

 

 

 

 

1

 

 

11

31

 

 

 

3438 7 343

21

 

 

 

31

 

21

 

 

 

 

238 87 7

 

 

 

 

 

 

 

 

11

3

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

1

 

 

 

211

 

11

 

 

3

 

1

3438 75 47

 

 

 

 

 

 

 

 

 

 

 

 

343 238 2 343

 

 

211

 

1

31

 

22 1

 

 

2111

 

 

 

 

22 1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2111

 

 

31

 

3438 75 47

 

1

 

 

 

 

 

 

 

 

 

 

 

2

1

1

 

1

 

3

1

 

 

 

 

 

 

 

 

 

343 238 87 7

 

 

 

 

 

 

11

 

 

 

 

1

211

 

 

 

 

 

 

 

 

22 1

 

 

 

22

1

 

 

 

W1 =

1

 

 

 

 

 

 

 

 

 

 

 

 

W

 

 

 

 

 

 

3438 4

1

13 1

 

 

 

 

 

 

12

 

 

 

 

 

 

1

2 9 343 238

 

 

2

 

 

112

 

 

1

 

 

 

 

 

21

1

 

 

11

111

 

 

3

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

21

11

 

 

 

 

 

3

 

 

10