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

Конспект лекц ТАУ _1ч _11 л

.pdf
Скачиваний:
60
Добавлен:
12.03.2015
Размер:
1.78 Mб
Скачать

41

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

Uу( p) Lя p Rя Iя( p) cе ( p) ,

Jp ( p) cмIя( p) Mc ( p) .

Исключая промежуточную переменную Iя( p) , получим выражение для ( p):

 

 

( p)

 

cм

 

 

 

 

U

у( p)

 

Lя p Rя

Mc ( p) .

(20)

 

 

 

JL p2 JR p c c

 

 

JL p2

JR

p c c

 

 

 

 

 

я

я

 

м е

 

 

 

 

я

 

 

я

 

м е

 

 

 

Обычно индуктивность Lя

малая величина и ею можно пренебречь. Тогда вы-

ражение

 

(20) с

учетом

принятых

обозначений k

дв

1/ c ,

k f

R

/(c c ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е

дв

я

м е

T

Jk f

 

можно записать в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дв

дв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p) W

 

( p)U

у

( p) W f ( p)M

c

( p) ,

 

 

(21)

 

 

 

 

 

 

дв

 

 

 

 

 

 

дв

 

 

 

 

 

 

 

где передаточные функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

дв

 

 

, W f

 

k f

 

 

 

 

 

 

 

 

 

 

W ( p)

 

 

 

 

( p)

 

дв

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дв

 

Tдв p 1

дв

 

Tдв p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

теля и знаменателя не нулевой степени имеют значение единица. Удобство та-

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

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

Tдв kдвuу kдвf Mc ,

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

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

эффициентами усиления.

42

Если в числителе и знаменателе передаточной функции W ( p) в изобра-

жениях Лапласа имеются одинаковые множители или, иначе говоря, одинако-

вые нули и полюса, то при определении оригинала y(t) по изображению (19)

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

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

Таким образом, передаточная функция W (s) в символьном виде и в изо-

бражениях Лапласа W ( p) совпадают с точностью обозначения s p , если не проводятся сокращения в передаточной функции W ( p) .

Вопросы для самопроверки

1. Чем отличаются изображения операций дифференцирования и интегрирова-

ния оригинала?

2.Какой оригинал соответствует произведению изображений?

3.В чем особенность начального значения оригинала, определяемого по из-

вестному его изображению?

4. В каких случаях по известному изображению можно найти конечное значе-

ние оригинала?

5. При каких условиях решения оригиналов являются асимптотически устойчи-

выми?

6. При каких начальных условиях определяется передаточная функция в преоб-

разованиях Лапласа?

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

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

43

ЛЕКЦИЯ 4

Структура движения системы. Типовые входные сигналы. Временные характери-

стики. Частотные характеристики. Свойства АФЧХ. Логарифмические частотные ха-

рактеристики и их свойства.

1. Структура движения системы.

В задачу анализа динамики САУ входит исследование движения системы,

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

деляется формулой (3.17). Полагая для простоты

f (t) 0, получим

 

 

 

 

b p b

 

 

 

 

y( 0) p y( 0) a y( 0) bu( 0)

 

 

 

Y( p)

1

0

 

U( p)

 

1

1

 

.

(1)

 

 

 

p2 a p a

 

 

p2 a p a

 

 

 

 

 

1

0

 

 

 

 

1

 

0

 

 

 

Выражение (1) представим в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y( p) Yсв( p) Yвын( p),

 

 

 

 

 

(2)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Yсв( p)

y(0) p y(0) a1y(0)

 

 

 

(3)

 

 

 

 

p2 a p a

0

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

соответствует свободному движению yсв(t) , вызванному ненулевыми началь-

ными условиями выходной координаты y(t) ,

 

b p b

 

 

bu(0)

 

 

 

Y ( p)

1

0

 

U( p)

1

 

.

(4)

 

 

 

 

 

вын

p2 a p a

p2 a p a

 

 

 

1

0

 

1

0

 

 

соответствует вынужденному движению увын(t), вызванному входным воздей-

ствием u(t) с учетом его начальных условий.

Если изображение входного сигнала U ( p) m( p)/ d( p), то выражение (4)

с помощью формулы разложения (2.21) или (2.24) можно представить в виде

Y ( p)

m1( p)

 

 

m2( p)

 

b1u(0)

 

.

(5)

p2 a p a

 

 

 

вын

 

d( p)

p2 a p a

 

 

1

0

 

 

 

1

0

 

 

44

Выражение (5) можно записать следующим образом

Yвын( p) Yсоб( p) Yуст( p) ,

(6)

 

 

где

Yсоб( p)

m1( p)

 

 

b1u(0)

 

(7)

p2 a p a

p2 a p a

1

0

1

0

 

соответствует собственному движению yсоб(t) , вызванному входным воздейст-

вием, но независящего от его вида,

Y ( p)

m2 ( p)

(8)

уст d( p)

соответствует установившемуся движению yуст(t), зависящего только от вида входного воздействия (полюсов уравнения d( p) 0 ).

Таким образом, структура полного движения выхода y(t) ФЭ с учетом

обратного преобразования Лапласа можно представить следующим образом

 

 

 

y(t) yсв(t) yвын(t) yсв(t) yсоб(t) yуст(t) yпер(t) yуст(t) ,

 

(9)

 

 

 

 

 

 

 

 

 

где yпер(t) yсв(t) yсоб(t) –

переходное

движение

зависит от полюсов

pi ,

i

 

(n 2) полинома оператора входа

 

 

 

 

1,n

 

 

 

 

 

 

 

d( p) p

2 a p a ( p p )( p p ) .

 

 

 

 

 

1

0

1

2

 

 

 

 

Представление (9) справедливо для решения дифференциального уравне-

ния произвольного порядка.

 

 

 

 

 

 

 

 

Согласно свойству 2.1

решения

yсв(t) , yсоб(t)

затухают при Re pi

0 ,

i 1,n, поэтому спустя некоторое время, называемое временем установления переходного процесса ty (или временем регулирования tр ), будет выполняться условие y(t) yуст(t) при t ty для произвольного входа u(t) .

Если имеются полюса pi , у которых Re pi 0 , то согласно свойству 2.2 с

течением времени решения yсв(t) , yсоб(t) расходятся и определить yуст(t) по измерениям выхода y(t) для произвольного входа u(t) не удается. Однако для

45

известного входа u(t) в этом случае можно подобрать начальные условия y(0), y(0), u(0), при которых в числителе изображения (1) появятся нули, равные не-

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

(9) будут отсутствовать указанные полюса и с течением времени y(t) yуст(t)

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

го сигнала, а также параметров системы приводят к проявлению свойства неус-

тойчивости системы. Отсюда следует важный вывод:

Для определения свойств асимптотической устойчивости движения сис-

темы с помощью эксперимента необходимо учитывать произвольные на-

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

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

2. Типовые входные сигналы.

Одним из способов анализа движения системы на практике является по-

строение ее реакции выхода на типовые входные воздействия u(t) , при нуле-

вых начальных условиях (значениях входа, выхода и их производных). В каче-

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

Различают следующие типовые сигналы.

1) Единичный скачок 1(t) , который с физической точки зрения можно

рассматривать, как предел lim (t, ) 1(t) функции (t, ), указанной на рис. 1.

0

Единичный скачок 1(t) используется для имитации скачкообразного из-

менения команды или возмущения.

46 2) Импульсная функция (t) , являющаяся пределом прямоугольного им-

пульса (t, ) d (t, )/ dt с единичной площадью (рис. 2)

при 0, которая

определяется по формулам

 

 

 

 

 

 

 

 

 

 

 

при t 0,

 

(t)dt 1.

 

(10)

 

(t)

 

 

 

0 при t 0,

 

 

 

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

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

Установим связь единичного скачка 1(t) и импульсной функции (t) . Для

этого найдем производную от функции 1(t) :

 

 

 

 

 

d

0

при

t 0;

 

 

 

 

 

 

 

 

 

 

 

1(t)

при

t 0.

.

 

 

 

 

 

 

dt

lim (t, )

 

 

 

 

 

 

0

 

 

 

Поскольку интеграл

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

1(t)dt d1(t) 1( ) 1( ) 1 0 1,

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

то справедливо выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

1(t) (t) .

 

 

(11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью выражений (10) найдем изображение Лапласа для функции

(t) :

 

 

 

L (t)

(t)e ptdt (0)e p0dt

(t)dt 1.

0

0

 

Данный результат также можно получить с помощью выражения (11):

d

 

 

1

 

L{ (t)} L

 

1(t)

p

 

1( 0) 1 0 1.

 

p

dt

 

 

 

Функцию (t) удобно использовать для задания начальных условий сис-

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

47

ставления системы:

y ay , y(0) y0 ;

y ay y0 (t) , y(0) 0 ,

получим одно и тоже изображение Лапласа

Y ( p) y0 . p a

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

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

мощью (t) функций через входные воздействия.

3) Гармонический сигнал, определяемый по формуле u(t) um cos t или

u(t) um sin t , используется для имитации периодических воздействий, а так-

же может выступать в качестве эквивалентного сигнала,

если для сигнала u(t)

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

и ускоре-

ние umax .

В этом случае эквивалентный сигнал

можно

записать в

виде

u

э

(t) u

max

sin

t , где

 

э

u

max

/u

max

. Учитывая,

что

u

(t) u

 

2 sin t ,

 

 

э

 

 

 

 

 

 

 

э

 

max э

э

должно выполняться условие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

э

 

 

 

.

 

 

 

 

 

(12)

 

 

 

 

 

 

 

 

 

 

umax /umax

 

 

 

 

 

 

 

Для определения реакции системы y(t) при нулевых начальных условиях

на типовое входное воздействие u(t)

можно воспользоваться выражением

 

 

 

 

 

 

 

 

 

 

Y ( p) W ( p)U ( p),

 

 

 

 

 

(13)

где передаточная функция W ( p) системы в общем случае имеет вид (11).

 

 

 

Тогда при известном изображении U ( p) реакцию

y(t) можно найти с

помощью обратного преобразования Лапласа

 

 

 

 

 

 

 

 

 

 

 

 

y(t) L 1 Y ( p) L 1 W ( p)U ( p) .

 

 

 

 

(14)

Для облегчения определения оригинала y(t) изображение Y ( p) с помо-

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

вующие оригиналы и затем складываются.

48

3. Временные характеристики.

Реакция выхода y(t) при нулевых начальных условиях, построенная для входного сигнала 1(t) или (t) с помощью выражения (2.14), называется вре-

менной функцией, а ее график – временной характеристикой.

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

рим их свойства.

3.1. Переходная характеристика

Определение 1. Реакция выхода системы y(t) на единичный скачок 1(t)

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

рая обозначается символом h(t) .

Учитывая, что для единичного скачка U ( p) 1/ p , переходная функция

согласно (14) определяется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(t) L 1 H ( p)

L 1

W ( p)

.

 

(15)

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По известному

 

 

 

 

 

 

конечное значение ори-

изображению H ( p)

можно найти

гинала h( ) lim h(t), если он существует, по формуле (3.10):

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h( ) lim pH ( p).

 

 

 

 

 

(16)

 

 

 

 

 

p 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.2. Весовая характеристика

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение 2. Реакция выхода системы y(t)

на импульсную функцию

 

(t) при нулевых начальных условиях называется весовой функцией, кото-

 

 

рая обозначается символом w(t) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

Учитывая, что для u(t) (t)

изображение U ( p) 1, весовая функция со

гласно (14) определяется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(t) L 1 W ( p) .

 

 

 

 

 

 

(17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

49

Установим связь переходной и весовой функций. Согласно (15) имеем

H ( p) W ( p) . p

Данное выражение перепишем в виде

pH ( p) W ( p) .

Тогда согласно свойству изображения для производной получим

 

 

 

 

 

 

 

 

 

d

h(t) w(t) .

 

(18)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

Уравнение (18) лежит в основе экспериментального метода определения весовой характеристики по измеренной переходной характеристике h(t) . Затем

по известной весовой характеристике

w(t) можно найти реакцию систему на

произвольное входное воздействие u(t) с помощью теоремы о свертке (3.6):

 

t

t

 

y(t) L 1 W ( p)U ( p) w(t )u( )d w( )u(t )d .

(19)

0

0

 

Вычисление интеграла (19) можно проводить методами численного ин-

тегрирования.

По экспериментально найденной весовой характеристике можно постро-

ить приближенную передаточную функцию [14].

Пример 1. Рассмотрим систему с передаточной функцией

W ( p) k ,

T p 1

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

Из формулы (1) получим

 

1

W ( p)

 

 

1

k

 

 

1

1

 

T

 

 

h(t) L

L

kL

 

k

1(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

p(T p 1)

 

 

p

 

T p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C помощью формулы (3) найдем

w(t) L 1 W ( p)

L 1

k

 

 

 

1

t

 

 

 

ke T .

 

 

 

 

 

 

 

 

 

 

 

 

 

T p

1

 

 

 

1 t

e T .

50

Решение аналогичной задачи для передаточных функций W ( p) общего вида можно проводить с помощью вычислительного пакета MATLAB [13]. При этом для построения характеристики h(t) используется команда step, а для построения весовой характеристики команда impulse. С помощью фрагмента программы для принятых значений параметров k 10, T 0,5 :

k=10;T=0.5;W=tf(k,[T 1]);step(W);hold on; impulse(W);grid

на рис. 3 приведены графики указанных характеристик.

Рис. 3

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

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

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

a0 y(n) (t) a1y(n 1) (t) an y(t)

(20)

b0u(m) (t) b1u(m 1) (t) bmu(t),

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

жениях Лапласа

Y ( p) W ( p)U ( p),

(21)

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