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

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

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

 

 

 

 

 

 

 

 

 

Окончание табл. 1

123452674869 3

 

 

 

 

 

 

2 296 3 6

 

 

 

 

 

4 96

 

 

8 86 39 96

 

 

 

 

 

 

 

 

 

 

 

 

1111

 

111

1

 

2 1

 

 

 

 

 

 

2 1

 

 

 

1232456 789 2328

 

11

1

1 1

 

 

 

2 1

1

 

67 53 232

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

11

1

2 1

1

12

1

1

111

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

11

1

 

1

 

 

11

1

 

1

 

 

 

 

 

2

1232456 789 2328

 

 

 

 

 

1

 

2

1

 

 

 

1

 

 

67 53 4 8

 

 

 

2 1

2

2 1

2

12

1

 

 

 

 

 

 

 

 

 

 

12

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

13

11

1

1

12

1

5

1

6

12

17

14

Рис. 1. Исходная структура САУ

Преобразованиевключаетнесколькоэтапов,показанныхнарис.2–5.

13

 

1

1

 

11

12

15

16

 

1

2

 

 

 

 

17

14 2316

Рис. 2. Перенос узла через сумматор

11

13

 

 

W5W6

 

 

 

W9 =

 

 

 

1 + W7W6W5

1

1

12

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W8 =

W4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W6

 

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3. Свертывание обратной связи и последовательного соединения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

111

 

 

112

 

 

 

 

 

 

 

 

 

1123

 

 

 

 

W10 =

 

W3W9

 

Wэкв

= W10W11

=

W1W3W5W6 + W2W3W5W6

 

W8W3W9

1 −W5W6W7 W3W4W5

1

 

 

 

 

 

 

 

 

 

W11 = W1 + W2

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4. Свертывание обрат7

 

Рис. 5. Свертывание последователь7

ной связи и параллельного

 

 

 

 

 

ного соединения

 

соединения

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

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

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

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

Контурыназываютсянекасающимися,еслионинеимеютобщихуз) лов.

Сигнальный граф однозначно соответствует структурной схеме.

12

Пусть X(s) и Y(s) – входная и выходная переменные системы. Тог) да для вычисления ПФ системы управления по ее графу можно вос) пользоваться формулой Мейсона:

 

 

N

 

 

X(s)

 

Pii

 

= W(s) =

i=1

 

,

 

 

Y(s)

 

где Pi i й путь от входа к выходу; N – количество путей; ∆– опреде) литель графа; ∆i – дополнительный множитель для пути.

Определитель графа получается по формуле:

K

M,Q

R,S,L

∆ =1− Lk +

LmLq

Lr LsLl +...,

k=1

m=1,q=1

r=1,s=1,l=1

K

где Lk – сумма коэффициентов передачи всех отдельных конту)

k=1

 

 

 

M,Q

 

ров;

LmLq – сумма произведений всех возможных комбинаций

m=1,q=1

 

 

 

R,S,L

из двух некасающихся контуров;

Lr LsLl – сумма произведе)

r =1,s=1,l=1

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

ров, касающихся этого пути.

Рассмотрим пример получения ПФ многоконтурной системы с ис) пользованием формулы Мейсона для структуры рис. 1, которой соот) ветствует граф, показанный на рис. 6.

От входа к выходу ведут два пути:

11

X(S)

12

13

14

15

2678

319 31

Рис. 6. Описание системы управления сигнальным графом

13

P1=W1W3W5W6.

P2= W2W3W5W6

В графе есть два контура:

L1=–W3W5W4

L2=–W5W6W7 .

Контур L1 касается контура L2, поэтому определитель графа вы) числяется по формуле:

∆ =1−(L1 + L2 ) .

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

 

 

1 = ∆2 =1.

 

 

 

 

 

Окончательно можно записать:

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Pii

 

W1W3W5W6 +W2W3W5W6

 

W(s) =

i=1

=

.

 

1−W W W +W W W

 

 

 

3

5

4

5

6

7

 

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

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

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

– series(w1,w2) – последовательное соединение динамических зве) ньев;

– parallel(w1,w2) – параллельное соединение динамических зве) ньев;

– feedback(w1,w2) – включение звена w2 в контур отрицательной обратной связи к w1;

– feedback(w1,w2) – включение звена w2 в контур отрицательной обратной связи звена w1;

– feedback(w1,w2,sign) – включение звена w2 в контур обратной связи звена w1 с указанием знака + или – (очевидно, feedback(w1,w2)= =feedback(w1,w2,)1));

Пример:

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

Transfer function: s + 2

14

——————

s^2 + 2 s + 2

>>w1=tf([1 2 3],[1 2 2])

Transfer function: s^2 + 2 s + 3

——————

s^2 + 2 s + 2

>>w2=series(w,w1)

Transfer function: s^3 + 4 s^2 + 7 s + 6

——————————————

s^4 + 4 s^3 + 8 s^2 + 8 s + 4

>>w3=parallel(w,w1)

Transfer function:

s^4 + 5 s^3 + 13 s^2 + 16 s + 10

————————————————

s^4 + 4 s^3 + 8 s^2 + 8 s + 4

>>w4=feedback(w,w1)

Transfer function:

s^3 + 4 s^2 + 6 s + 4

————————————————

s^4 + 5 s^3 + 12 s^2 + 15 s + 10

Дляпроверкиправильностипроведенныхпреобразованийнеобходимо собратьсхемуисходнойСАУисоответствующуюейэквивалентнуюсхемув MatLabSimulink.Задачасчитаетсярешенной,еслиприподаченавходобо) ихсхемодинаковыхтестовыхвоздействийнаблюдаютсяодинаковыевы) ходныесигналы.

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

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

вэквивалентнуюПФдвумяспособами:

непосредственноиспользуяправилатабл.1;

используяпредставлениеввидесигнальногографаиформулуМейсо)

на.

Вариантызаданийприведенывтабл.2.

ВкачествезвеньевW1–W3использоватьтиповыединамическиезвенья спараметрами,выбраннымивлабораторнойработе№1.

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

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

15

1

212

3

4

5

6

7

Таблица 2

Варианты структур САУ

3456789 45 9

22

 

 

1

 

21

1

23

 

 

1

2

 

2

3

 

 

 

 

 

 

 

 

22

 

 

 

 

1

 

21

1

23

 

1

1

 

 

2

 

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

1

 

21

1

23

 

1

1

2

 

2

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

1

 

21

1

23

 

1

1

2

 

2

1

3

 

 

 

 

 

 

 

 

22

 

 

 

1

 

 

1

23

 

1

 

 

 

21

 

 

1

2

 

2

1

3

 

 

 

 

 

 

 

сигнальный граф системы, описание путей, контуров и расчет ПФ по формуле Мейсона;

схемы экспериментов в Simulink MatLab и протокол команд MatLab.

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

16

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

ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ДИНАМИЧЕСКИХ ЗВЕНЬЕВ

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

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

Пусть входной сигнал системы имеет амплитуду a и частоту ω, т. е. описываетсяформулой

x = asin(ωt).

Выходной сигнал будет иметь амплитуду А1 и отличаться от вход) ного по фазе на величину ψ(фазовый сдвиг):

y = A1 sin(ωt +ψ).

Таким образом, можно рассчитать усиление по амплитуде

A = A1 . a

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

A(ω) – амплитудную частотную характеристику (АЧХ) и ψ(ω) фазо) вую частотную характеристику (ФЧХ).

Главное достоинство метода частотных характеристик заключа) ется в том, что АЧХ и ФЧХ объекта могут быть получены экспери) ментально. Для этого необходимо иметь генератор гармонических колебаний, который подключается к входу объекта, и измеритель) ную аппаратуру для измерения амплитуды и фазового сдвига колеба) ний на выходе объекта.

Частотные характеристики САУ могут быть получены по ее ПФ W(s). Для суждения о реакции звена на синусоидальный сигнал дос) таточно исследовать его реакцию на гармонический сигнал вида [1]

x(jω) = ejωt.

Тогда выходной сигнал

y(jω) = A(ω)ej(ωt(ω)),

17

и частотная ПФ

W(jω) = y(jω) = A(ω)ejψ(ω). x(jω)

Формально для получения частотной ПФ надо сделать в W(s) под) становку s = jω , и тогда полученная W(jω) является комплексным выражением, которое можно представить в виде:

W(jω) = a1(ω) + jb1(ω) . a2 (ω) + jb2(ω)

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

W(jω) =

a1(ω) + jb1(ω)

=

(a1(ω) + jb1(ω)) (a2 (ω) − jb2 (ω))

=

a (ω) + jb (ω)

 

 

 

 

 

 

 

 

 

 

 

 

 

(a (ω) + jb (ω)) (a (ω) − jb (ω))

 

 

 

 

2

 

2

 

1

 

 

2

 

 

 

2

 

 

2

 

 

 

=

a1(ω)a2 (ω) +b1(ω)b2 (ω)

+ j

a2 (ω)b1(ω) −a1(ω)b2 (ω)

=

 

 

 

 

 

 

 

 

 

 

 

 

a2

(ω) +b2

(ω)

 

 

 

 

 

 

a2

(ω) +b2

(ω)

 

2

2

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

=U(ω) + jV(ω) = A(ω) ejϕ(ω),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

+b2

 

 

A(ω) =

W(jω)

= U2 (ω) +V2

(ω) =

 

 

 

1

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

+b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V(ω)

 

b

 

 

 

 

b

 

ψ(ω) = arg(W(jω)) = arctg

 

 

= arctg

 

1

 

 

−arctg

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U(ω)

 

a1

 

 

a2

Графики функций U(ω) и V(ω) называют соответственно веще ственной и мнимой частотной характеристиками.

В практических расчетах удобно применять графики частотных характеристик, построенных в логарифмическом масштабе – лога рифмические частотные характеристики (ЛЧХ).

Логарифмическая амплитудная частотная характеристика

(ЛАЧХ) определяется следующим выражением:

L(ω) = 20 lgA(ω).

18

Логарифмической фазовой частотной характеристикой(ЛФЧХ) называется график зависимости ψ(ω), построенный в логарифмичес) ком масштабе частот.

Единицей L(ω) является децибел (дБ), а единицей логарифма часто) ты – декада. Декадой называют интервал частот, на котором частота изменяется в 10 раз. При изменении частоты в 10 раз говорят, что она измениласьнаоднудекаду.ОсьординатприпостроенииЛЧХпроводят через произвольную точку, а не через точкуω= 0. Частоте ω= 0 соответ) ствует бесконечно удаленная точка: lgω→–∞при ω→0.

Основное преимущество использования ЛЧХ заключается в том, что приближенные (асимптотические) ЛАЧХ типовых динамических звеньев изображаются отрезками прямых.

Пример. Построим ЛЧХ апериодического звена первого порядка. Передаточная функция звена

 

W (s)

=

 

k

 

.

 

 

 

Ts +1

 

 

 

Частотная передаточная функция

 

 

 

 

W(jω) =

k

 

 

=

k(1−Tjω)

,

 

 

 

 

 

 

 

Tjω+1

(Tω)2 +1

U =

k

 

,

V = −

 

kTω

.

(Tω)2

 

 

 

 

+1

 

 

 

 

(Tω)2 +1

Следовательно, АЧХ описывается формулой

A(ω) =

k

,

(Tω)2 +1

ФЧХ строится по формуле

ψ(ω) = −arctg(Tω).

ЛАЧХ апериодического звена 1)го порядка

L(ω) =20lgk −20lg (Tω)2 +1.

По этой формуле можно построить две асимптоты – прямые, к которым стремится ЛАЧХ при ω→ 0 и ω→ ∞. Так, при ω→ 0 второе слагаемое близко к нулю, и этот участок ЛАЧХ представляет собой горизонтальную прямую

L(ω) =20lgk.

19

При ω →∞получаем наклонную прямую:

L(ω) ≈ −20lg (Tω)2 +1 → −∞.

Для определения наклона этой прямой можно рассмотреть грани) цы декады:

ω= 1 и ω= 10.

T T

Изменение ЛАЧХ между этими точками:

 

10

2

 

 

1 2

 

L(ω) ≈ −20lg T

 

 

+1

+20lg

T

 

 

+1

≈ −20 (дБ/дек).

 

 

 

T

 

 

 

T

 

 

ЛЧХ часто называют диаграммами Боде.

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

Впакете MatLab ЛЧХ объекта, заданного с помощью ПФ, можно получить с командой bode.

 

 

 

43

 

 

53

8

 

 

 

 

 

9

3

8

7

 

6

 

5

 

4

 

153

3

 

2

 

1

143

123

 

3

149

3

9

8

91 3

2

1689

1673

 

 

6312

631

632

Рис. 1. Частотные характеристики динамических звеньев

20