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

vostrikov

.pdf
Скачиваний:
133
Добавлен:
23.01.2019
Размер:
12.99 Mб
Скачать

2.4. Импульсная переходная функция

31

2.4. ИМПУЛЬСНАЯ ПЕРЕХОДНАЯ ФУНКЦИЯ

Эта характеристика также используется для описания одноканальных объектов вида (2.5).

Импульсная переходная функция (характеристика) g(t) пред-

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

1)

(t

)

0,

t

,

2)

(t

)d

1.

(2.7)

,

t

;

 

 

 

 

 

 

 

 

u u

t

y y

g(tt ) )

t

Рис. 2.5. Пример импульсной переходной функции системы

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

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

t

y(t)

g(t )u( )d .

(2.8)

0

32 Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

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

 

t

 

g( )d .

(2.9)

g(t) h(t), h(t)

 

0

 

Уравнения (2.9) позволяют при одной известной характеристике определить вторую.

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

2.5. ПЕРЕХОДНАЯ МАТРИЦА

Данная динамическая характеристика применяется для описания многоканальных систем вида (2.1) и (2.2).

x

Ax,

x

Rn ,

n m .

(2.10)

 

 

 

Rm ,

y

Cx,

y

 

 

Переходная матрица представляет собой решение матричного дифференциального уравнения

A , dim (t) n n

(2.11)

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

 

1

0

0

(0) I , где I

0

1

0

 

 

.

 

0

0

1

Она обладает следующими свойствами:

1) det (t) 0

для любого t 0, ),

(2.12)

2) 1(t) (t).

2.5. Переходная матрица

33

Зная переходную матрицу, можно вычислить реакцию системы x Ax Bu

на произвольное входное воздействие при любых начальных условиях x(0) по выражению

t

x(t)

(t)x(0)

( )Bu( )d .

(2.13)

0

Здесь первое слагаемое описывает свободную составляющую движения, второе – вынужденную. Соотношение для выходных переменных следующее:

 

t

 

 

y(t) C (t)x(0)

C

( )Bu( )d .

(2.14)

 

0

 

 

Если система имеет нулевые начальные условия x(0)

0 , то выра-

жение (2.14) принимает вид

 

 

 

t

 

 

 

y(t) G(

)u(

)d .

(2.15)

0

 

 

 

Матрица G(t, ) называется матричной импульсной переходной функцией. Каждая ее компонента представляет собой импульсную переходную функцию gij ( ) , которая является реакцией i-го выхода

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

G( ) C ( )B.

(2.16)

Для многоканальных систем может быть определена также мат-

ричная переходная характеристика

t

H ( ) G( )d .

(2.17)

0

34 Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

Для линейных систем с постоянными параметрами переходная мат-

рица (t) представляет собой матричную экспоненту

 

(t) eAt I At

1

 

A2t2

,

(2.18)

2!

 

 

 

 

где dimeAt n n .

С учетом (2.18) выражения (2.13) и (2.14) можно записать следующим образом:

 

 

t

 

 

x(t)

eAt x(0)

eA(t

) Bu( )d ,

(2.19)

 

 

0

 

 

 

 

t

 

 

y(t)

CeAt x(0)

CeA(t

) Bu( )d .

(2.20)

 

 

0

 

 

В этом случае матричная импульсная переходная функция линейной системы с постоянными коэффициентами может быть найдена по соотношению

G(t, ) CeAt B.

(2.21)

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

2.6. ПЕРЕДАТОЧНАЯ ФУНКЦИЯ

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

p dtd ,

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

35

что позволяет записывать дифференциальные уравнения как алгебраические и вводить новую динамическую характеристику – передаточную функцию. Этот способ был предложен английским ученым Хевисайдом в 1895 г., позднее он был строго обоснован аппаратом интегральных преобразований Лапласа и Карсона [11] в предположении нулевых начальных условий.

Рассмотрим этот переход для многоканальных систем общего вида

x

Ax Bu,

x Rn , u Rm ,

y

Cx,

y R

m

, n m.

 

 

 

Запишем уравнение состояния в операторной форме px Ax Bu,

что позволяет определить вектор состояния

x

( pI

A) 1 Bu

(2.22)

и выходные переменные системы

 

 

y

C( pI

A) 1 Bu.

(2.23)

Матрица взаимосвязи между выходными переменными и управляющими воздействиями в выражении (2.23) называется матричной

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

 

 

 

W ( p) C( pI A) 1 B.

 

(2.24)

Она имеет размерность m m :

 

 

 

 

 

W11( p) W1m ( p)

 

 

W ( p)

 

 

 

,

(2.25)

 

Wm1( p) Wmm ( p)

 

 

где Wij ( p) yi u j скалярные

передаточные

функции,

которые

представляют собой отношение выходной величины к входной в символической форме при нулевых начальных условиях i 1, m, j 1, m .

36 Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

Собственными передаточными функциями i-го канала называ-

ются компоненты передаточной матрицы Wi i ( p) yi ui , которые на-

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

Как известно, обратная матрица ( pI

A) 1

 

может быть найдена по

выражению

 

 

 

 

 

( pI A)

1

( pI

A)

,

(2.26)

 

 

 

det( pI

A)

 

 

 

 

где ( pI A) – присоединенная матрица. Как следует из (2.26), все

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

A( p) det( pI A) pn a pn 1

a ,

n

1

который называется характеристическим полиномом и имеет n-й по-

рядок.

Если теперь характеристический полином приравнять нулю, то получим характеристическое уравнение системы

A( p) det( pI A) 0.

(2.27)

Уравнение (2.27) имеет n корней, которые называются полюсами системы p1,..., pn .

ПРИМЕР 2.6

Определить передаточную матрицу для объекта

 

 

 

x

Ax Bu,

x R2 , u R2 ,

 

 

 

y

Cx,

 

 

y

R

2

,

 

 

 

 

 

 

 

 

где A

0

1

, B

2

0

, C

1

1 .

 

 

1

2

 

0

1

 

0

2

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

37

Воспользуемся выражением для передаточной матрицы (2.24) и найдем

предварительно обратную матрицу (2.26). Здесь

 

 

 

 

 

 

 

 

p

 

 

1

 

 

 

 

 

 

 

 

pI

A

1

( p

2) .

 

 

 

 

 

 

Присоединенная матрица имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

( pI A)

 

p

 

 

1

 

 

,

 

 

 

 

1

 

p

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а характеристический полином –

 

 

 

 

 

 

 

 

 

 

 

 

 

 

det( pI

A)

 

p2

2 p

1.

 

 

 

 

 

Определим теперь обратную матрицу

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( pI A) 1

p2

2 p 1 p2

2 p 1 ,

 

 

 

 

1

 

 

 

 

 

p

2

 

 

 

 

 

p2

p2

1

 

p2

2 p

1

 

 

и передаточную матрицу объекта в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2( p

1)

 

 

 

 

p 1

 

 

 

 

 

 

 

 

 

 

W ( p) C( pI A) 1 B

 

p2

2 p 1 p2

2 p 1 .

 

 

 

 

4

 

 

 

 

 

2( p

2)

 

 

 

 

p2

2 p

1

 

p2

2 p

1

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

Чаще всего передаточные функции применяются для описания одноканальных систем вида (2.5)

y(n) a y(n 1)

a y

b u(m)

b u, n m .

(2.28)

n

1

m

0

 

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

38 Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

 

 

y

 

b

 

pm

b p

b

 

 

 

W ( p)

 

 

m

 

 

1

0

,

(2.29)

 

u

 

pn

an pn 1

a2 p a1

 

 

 

 

 

где A( p) pn

a pn 1

a p

 

a

характеристический полином.

 

n

 

2

 

1

 

 

 

 

 

Его корни,

p1, p2 , , pn , называются полюсами, а корни поли-

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

N

n1, n2 , , nm

, называют-

ся нулями системы.

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

 

 

 

 

 

d

m

pm

d p

1

 

 

 

 

 

 

W ( p) k

 

 

1

 

,

(2.30)

 

 

c pn

c p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

k b0 a1

– коэффициент

 

усиления; c j

 

a j 1 a1 ,

j 1, (n 1);

 

 

 

 

 

 

 

 

 

 

 

 

 

cn

1 a1 , di

bi b0 , i 1, m .

 

 

 

 

 

 

 

 

Передаточную матрицу (передаточную функцию) можно также определить с помощью изображений Лапласа или Карсона–Хевисайда. Если подвергнуть одному из этих преобразований обе части дифференциального уравнения и найти соотношения между входными и выходными величинами при нулевых начальных условиях, то получим ту же самую передаточную матрицу (2.24) или функцию (2.29).

Все динамические характеристики объекта взаимосвязаны: получив одну из них, можно определить все остальные. Мы рассмотрели переход от дифференциальных уравнений к передаточным функциям с помощью оператора дифференцирования p. Используя этот оператор, несложно сделать и обратный переход от передаточной функции к символической форме записи дифференциального уравнения, а затем к стандартному описанию объекта в форме (2.3) или (2.5).

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

 

t

y(t)

g(t )u( )d .

0

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

39

Подвергнем его преобразованию Лапласа [2, 9, 12]:

t

L[ y(t)] L g(t )u( )d

0

и получим соотношение y(s)

g(s)u(s) , из которого определим g(s) в

виде

 

 

 

 

g(s)

 

y(s)

W (s).

(2.31)

 

 

 

u(s)

 

 

 

 

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

ПРИМЕР 2.7

Определить передаточную функцию, нули и полюса для объекта, модель которого задана уравнением

y 6y 5y 2u 12u .

Запишем исходное уравнение объекта в операторной форме с помощью оператора дифференцирования p

 

( p2 6 p 5) y

(2 p

12)u .

 

 

Определим теперь передаточную функцию

 

 

 

 

 

W ( p)

y

2 p

12

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

u

p2

 

 

 

 

 

 

 

6 p

 

5

 

 

Характеристическое уравнение объекта имеет вид

 

 

 

A( p) p2

6 p

5

 

0.

 

 

Передаточная функция содержит два полюса ( p1

5 , p2

1 ) и один

нуль n1

6 .

 

 

 

 

 

 

 

 

 

 

ПРИМЕР 2.8

Определить передаточную функцию двигателя постоянного тока с независимым возбуждением (см. рис. 2.2).

40

Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ

Дифференциальное уравнение двигателя получено в примере 2.4 и

имеет вид

 

 

 

 

 

TяTм y

Tм y

y ku kм (Tя p 1)M .

Будем

полагать, что

возмущающее

воздействие отсутствует, т. е.

M 0 . Запишем это уравнение в символической форме с помощью опера-

тора дифференцирования p

 

 

 

T T p2 y T py

y ku

 

 

я м

м

 

или, рассматривая его как алгебраическое:

T T p2

T p 1 y ku .

я м

м

Определим теперь передаточную функцию двигателя постоянного тока с независимым возбуждением

W ( p)

 

k

 

.

 

 

 

T T p2

 

 

 

 

T p 1

 

я м

 

м

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

ными или комплексно-сопряженными.

2.7. МОДАЛЬНЫЕ ХАРАКТЕРИСТИКИ

Модальные характеристики соответствуют свободной составляющей движения системы (2.1) или, другими словами, отражают свойства автономной системы (2.10)

x Ax ,

x

Rn .

(2.32)

Будем искать ее решение в виде экспоненты

 

x(t)

e

t ,

(2.33)

где e t – скалярная экспонента,

x(0) – вектор начальных условий.

Подставляя решение (2.33) в исходное уравнение (2.32), после преобразований получим

[ I A]

0.

(2.34)

Соседние файлы в предмете Основы Теории Управления