Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
1787.pdf
Скачиваний:
42
Добавлен:
07.01.2021
Размер:
1.81 Mб
Скачать

Теперь, используя полученные коэффициенты, можно записать окончательное линейное дифференциальное уравнение по (2.6):

5 y 2 y 3 x 2 x 0.

(2.13)

2.4. Преобразования Лапласа

Исследование АСР существенно упрощается при использовании прикладных математических методов операционного исчисления. Например, функционирование некоторой системы описывается уравнением вида

a

d2y

a

dy

a y b

dx

b x,

(2.14)

 

 

 

0

1 dt

2

0

И

 

dt2

dt

1

 

где х и у – входная и выходная величины. Если в данное уравнение

 

 

Д

 

вместо x(t) и y(t) подставить функции X(p) и Y(p) комплексной пере-

менной величины p, такие, что:

 

 

 

 

А

 

 

X(p) x(t)e ptdt

и Y(p) y(t)e ptdt,

(2.15)

0

б

 

0

 

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

 

и dt

 

 

 

линейному алгебраическому уравнению

 

 

 

a0

p2 Y(p) + a1 p Y(p) + a2 Y(p) = b0 p X(p) + b1 X(p);

 

(2.16)

 

С

p

d

.

 

 

 

(2.17)

 

 

 

 

 

 

Переход от одной модели к другой достаточно прост и заключа-

 

 

 

 

 

dn

n

 

 

ется в замене знаков дифференциалов

 

на операторы p

,

знаков

dtn

 

 

 

 

 

 

 

 

интегралов ...dt на множители 1 , а самих x(t) и y(t) – изображения- p

ми X(p) и Y(p).

Такой переход от дифференциального к алгебраическому урав-

нению называется прямым преобразованием Лапласа, формулы

(2.15) – соответственно формулами прямого преобразования Лапласа,

аполученное уравнение – операторным уравнением или уравнением

воператорной форме.

19

Преобразование Лапласа обозначается как

 

L x(t) X p .

(2.18)

При этом x(t) называют оригиналом, а X(p) – изображением.

Для обратного перехода от операторного уравнения к функциям от времени используется метод обратного преобразования Лапласа.

Общая формула обратного преобразования Лапласа

x(t) L 1 X(p)

1

 

 

X(p)eptdt.

(2.19)

 

 

2

 

Эта формула достаточно сложна, поэтому были разработаны специальные таблицы (табл. 2.2), в которые сведены наиболее часто встречающиеся изображения функций X(p) и их оригиналы x(t). Данные формулы могут использоваться как в прямом, так и в обратном направлении. Они позволяют отказаться от использования формул

(2.15) и (2.19).

 

 

 

 

И

 

 

 

 

 

 

 

 

 

Таблица 2.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формулы преобразований Лапласа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

Функция

 

 

x(t) ↔ X(p)

идеальное запаздывание

А

(t ) e p

 

 

 

 

 

1(t)(t )

e p

 

единичная функция с запаздыван ем

 

 

 

 

 

 

 

 

 

 

p

 

б

 

 

 

 

 

 

 

 

 

 

 

 

дельта-функц я

 

 

 

1

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

степенная функция времени

 

 

tn

 

n!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pn 1

 

единичная функция

 

 

1(t)

1

 

 

 

 

 

 

 

p

С

 

 

 

 

 

 

 

 

«ступенька скорости»

 

 

1(t)t

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

p

 

квадрат времени

 

 

 

1(t) t2

2

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

p

 

куб времени

 

 

 

1(t) t3

6

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

p

 

экспоненциальное затухание

 

 

Me t

 

 

 

M

 

 

 

 

p

 

 

 

 

 

 

умножение на = const

 

 

x(t) ∙X(p)

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

Функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) ↔ X(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

сложение функций

 

 

 

 

 

 

 

 

xi(t)

Xi(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dnx(t)

 

 

 

n

 

 

 

 

 

 

 

дифференцирование функции

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

X(p)

 

 

 

 

 

 

 

 

dtn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

X(p)

 

 

 

интегрирование функции

 

 

 

 

 

 

 

 

 

x(t)dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

синус функции

 

 

 

 

 

 

 

 

 

 

sin( t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

косинус функции

 

 

 

 

 

 

 

 

 

cos( t)

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для примера рассмотрим решение дифференциального уравне-

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ния с использованием преобразований Лапласа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2 y

5

dy

6y 2

dx

12x.

 

 

 

 

 

 

 

 

(2.20)

 

 

dt2

 

 

А

dt

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Допустим, входной сигнал имеет форму единичного ступенча-

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

того воздействия, т. е. x(t) = 1(t). Тогда изображение входного сигнала

X(p) = 1/р.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Производим прео разование исходного дифференциального

уравнения по Лапласу подставляем в него X(p) и Y(p):

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.21)

 

p2Y(p) + 5∙pY(p) + 6∙Y(p) = 2∙pX(p) + 12∙X(p),

 

 

p2Y(pи) + 5∙pY(p) + 6∙Y(p) = 2∙p

1

 

+ 12∙

1

,

 

 

 

 

 

(2.22)

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

Y(p)∙(p3 + 5∙p2 + 6∙p) = 2∙p + 12.

 

 

 

 

 

 

 

 

(2.23)

Определяется выражение для Y:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y p

 

 

2p 12

.

 

 

 

 

 

 

 

 

 

 

 

 

(2.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p3 5p2

 

6p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оригинал полученной функции отсутствует в таблице оригиналов и изображений. Для решения задачи его поиска дробь разбивается на сумму простых дробей с учетом того, что знаменатель может быть представлен в виде p(p + 2)(p + 3):

Y p

2p 12

 

 

2p 12

 

2

 

4

 

2

. (2.25)

p3 5p2

 

p(p 2)(p 3)

 

 

 

 

6p

 

p

p 2

p 3

 

 

 

 

21

 

 

 

 

 

 

 

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

y(t) = 2 – 4∙e-2t + 2∙e-3t.

(2.26)

Таким образом, описание линейных непрерывных систем может быть представлено через преобразование Лапласа в операторной форме, что наиболее часто применяется в ТАУ.

2.5. Передаточная функция

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

Передаточной функцией называется отношение изображения выходного сигнала Y(p) к изображению входного воздействия X(p) при нулевых начальных условиях:

 

 

 

 

 

 

 

W(p)

 

Y(p)

.

И

 

 

 

 

(2.27)

 

 

 

 

 

 

 

 

X(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Передаточная функция являетсяДдробно-рациональной функци-

ей комплексной переменной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

m 1

 

 

 

 

 

m 2

 

 

 

 

 

 

 

 

 

 

 

B(p)

 

b0 p

 

А

b2 p

 

... bm

 

 

W(p)

 

 

 

b1 p

 

 

 

,

(2.28)

A(p)

 

 

n

 

 

n 1

 

 

 

 

 

 

 

 

a

p

 

a p

 

a

2

p

n 2

...

 

a

n

 

 

 

 

 

 

 

б0

1

 

 

 

 

 

 

 

 

 

 

 

 

где B(p) = b0 pm + b1 pm-1

+ b2 pm-2 + … + bm – полином числителя,

 

А(p) = a0 p

n

+ a1 p

n-1иn-2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a2 p

+ … + an – полином знаменателя.

 

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

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

порядком полинома знаменателя n. Из (2.27) следует, что изображение выходного сигнала можно найти как

Y(p) = W(p)∙X(p).

(2.29)

Для линейных систем с входным X и возмущающим воздействием F (рис. 2.1) можно применить принцип наложения (суперпозиции)

ивыделить следующие два случая:

сигнал F(p) = 0, тогда A(p) Y(p) B(p) X(p);

сигнал X(p) = 0, тогда A(p) Y(p) Q(p) F(p).

Тогда для такой АСР, имеющей входы по управлению и по возмущению, можно определить две передаточные функции

22

WX

(p)

 

Y(p)

 

B(p)

;

(2.30)

 

 

 

 

 

 

 

 

X(p) A(p)

 

WF

(p)

Y(p)

 

 

Q(p)

.

(2.31)

F(p)

 

 

 

 

 

 

A(p)

 

Уравнение (2.30) представляет передаточную функцию по управлению, а выражение (2.31) представляет передаточную функцию по возмущению. Общая передаточная функция такой системы является суммой (2.30) и (2.31).

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

Для линейных систем при нулевых начальных условиях возможен взаимный однозначный переход между операторной и частотной передаточными функциями в непрерывном времени. Для таких процессов преобразование Лапласа переходит в преобразование Фурье, т.е. производится замена оператора Лапласа p на частотный ком-

комплексную параметрическуюАфункциюД, модуль которой равен отношению амплитуды выходного сигнала к амплитуде входного сигнала А( ), а аргумент – сдвигу фаз между выходным и входным сигналом φ( ).Частотная передаточная функция W(jω), её модуль и аргумент φ( ) являются частотными характер стиками системы (см. п. 2.10).

плексный оператор jω.

И

 

Частотная передаточная функция W(jω) представляет собой

Частотную передаточнуюбфункцию, как и любое комплексное

число, можно представ ть в показательной или алгебраической форме

и

 

с вещественной U(ω) и мнимой V(ω) частями:

 

W( j ) A( )ej ( ) U( ) jV( );

(2.32)

С

 

 

 

 

 

A( )

 

W( j )

 

U2( ) V2( );

(2.33)

 

 

( ) argW( j ) arctg

V( )

.

(2.34)

 

 

 

 

 

 

U( )

 

Например, для апериодического звена (см. п. 2.11) в передаточную функцию вместо оператора p подставляется выражение jω. Затем, чтобы выделить вещественную U(ω) и мнимую V(ω) части, нужно освободиться от мнимости в знаменателе, умножив числитель и знаменатель функции на сопряженный комплекс:

23

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]