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

ЦСУ_Егоров

.pdf
Скачиваний:
16
Добавлен:
12.02.2016
Размер:
2.4 Mб
Скачать

6 ОБРАБОТКА СИГНАЛОВ, ПРЕДСТАВЛЕННЫХ ЧИСЛОВЫМИ ПОСЛЕДОВАТЕЛЬНОСТЯМИ

6.1 Типовые воздействия

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

Рассмотрим математическое описание цифровых систем управле- ния применяемое для расчетов регуляторов и их моделирования.

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

Цифровой единичный импульс (рис. 6.1 а).

1 [n] = 1,

при n = 0

(6.1)

 

0

0, при n ¹ 0

 

Цифровой единичный импульс задержанный на

мер, n0 = 2 (рис.6.1, б).

 

 

 

 

 

1, при n = n

o

1

 

[n n

 

 

 

 

] =

 

 

0

 

0

0, при n ¹ n

 

 

 

 

 

o

n0 отсчетов, напри-

(6.2)

102

а)

 

 

10[n]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-3

 

-2

 

-1

 

0

 

1

 

2

 

3

 

 

n

б)

10

[n n0

]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-3

-2

-1

0

1

2

3

4

n

Рисунок 6.1 Цифровой единичный импульс

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

Единичный скачок 1[n] (рис. 6.2) задаётся следующим образом

1[n]= 1, при n ³ 0

 

0, при n < 0

 

Единичный скачок связан с единичным импульсом соотношением

1[n] =

 

10[n−m].

(6.6)

 

m=o

 

1[n]

n

Рисунок 6.2. Цифровой единичный скачок. Следующей важной последовательностью является синусоидальная

(рис. 6.3).

x[n] = А× sin(ωnT + ϕ ) n = 0,1... ,

(6.7)

где A- амплитуда, ϕ - начальная фаза.

 

103

x[n]

nT с

Рисунок 6.3. Дискретная синусоида x[nT ] = Asin(ωnT + ϕ )

Синусоидальная последовательность может быть заменена симво- лической записью в комплексном виде

x[n] = e j(ωnT ).

(6.8)

Символичность записи заключается в том, что на самом деле дис- кретная синусоида определяется мнимой частью приведенной комплекс- ной функции.

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

Рассмотрим числовую последовательность a[0], a[1], a[2]...a[n],

где а[ n] - значение n-го элемента. Такая последовательность опи- сывается равенством:

a[m]10

(6.9)

a[n] =

[n m].

m = −∞

 

В равенстве (6.9) используется задержка единичного импульса на m отсчётов и масштабирование его умножением на множитель a[ m] .

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

Эффективный математический метод описания дискретных функ- ций основывается на замене последовательности импульсов х[nT] на ре- шетчатую функцию, состоящую из идеальных импульсов. Такая решетча- тая функция физически не реализуема, она служит лишь удобной формой представления реальных импульсов.

Будем считать, что преобразование непрерывного сигнала в реаль- ный дискретный можно реализовать с помощью схемы (рис. 6.4).

104

Если период замыкания и размыкания ключа Т значительно боль- ший чем продолжительность времени положения ключа в замкнутом со- стоянии tз, то дискретная последовательность х[nT], имеющая площадь x[nT]tз можно приближенно заменить идеальными импульсами δ (t) той

же площади.

 

 

Т

 

 

 

x (t)

 

 

 

x [nT]

 

x (t)

 

a)

 

 

 

0.3

 

 

 

 

 

 

 

 

 

 

0.23

 

 

 

 

Yt()

0.15

 

 

 

 

 

0.075

 

 

 

 

 

0

0.5

1

1.5

2

 

0

 

 

 

t

 

 

 

 

 

б)

 

 

г)

Рисунок 6.4. Преобразование непрерывного сигнала в реальный дискретный а) n; б) непрерывный сигнал; г) реальный сигнал.

Последовательность идеальных импульсов описывается выражением:

 

x[δ (t)] = x[nT ]δ (t nT )

(6.10)

n=0

x [nT] и

В дальнейшем будем считать, что последовательности

x [δ (t)]примерно равны

 

x [nT] x[δ (t nT )] .

(6.11)

Преобразование Лапласа заключается в том, что функция веще- ственной переменной x (t) преобразуется в функцию F(p) комплексной пе- ременной p = a + jω . Исходная функция называется оригиналом, а пре-

образованная изображением.

105

F ( p) = L{f (t)}= f (t)ept dt

 

(6.12)

0

 

 

 

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

L{δ (t)}= δ (t)ept dt = 1 ,

 

(6.13)

0

 

 

 

а применение этого преобразования к импульсу, сдвинутому на nT

тактов дает

 

 

 

L{δ (t nT )}= δ (t nT )ept dt .

 

(6.14)

0

 

 

 

Введем новую переменную λ = t nT , тогда t = λ + nT . С уче-

том принятых обозначений выражение (6.14) примет вид

 

 

 

L{δ (t nT}= δ (λ)epλ enT dλ = epnT δ (λ)epλ dλ = epnT

(6.15)

0

 

0

 

С учетом допущения соотношений преобразование Лапласа реаль-

ной дискретной последовательности приобретает вид

 

 

 

 

L{x[nT ]}= x[nT ]epnT .

 

(6.16)

n=0

 

 

 

Такое преобразование называется дискретным преобразованием

Лапласа.

 

 

 

Обозначив Z(p) = e pt

получим Z-преобразование

 

 

 

 

Z{X [nT ]}= x[nT ]z

n .

 

(6.17)

n=0

6.3Математические модели линейных цифровых систем

Цифровая система (звено), с точки зрения её математического опи- сания, по существу является преобразованием входной последовательно- сти в выходную в соответствии с заданным алгоритмом. Обозначим вход- ную последовательность через x[n], а выходную через y[n]. Если на вход

системы подать некоторые входные последовательности x1[n], x 2 [n]...x n [n], то на выходе будут некоторые выходные последова-

тельности соответственно y1[n], y 2 [n]...y m [n], которые называются откли-

ками системы.

Линейной называется такая система, у которой при подаче на вход суммы двух и более последовательностей x1[n]+ x2[n], на выходе образу-

ется последовательность y1[n]+ y2[n], .

106

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

тельность y[n], то выходной последовательности x n - n 0 при любых n0 ,

соответствует последовательность y[n − n0 ].

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

В разностных уравнениях аналогом первой производной является либо первая прямая разность:

Dy[n] = y[n + 1] - y[n],

(6.18)

либо первая обратная разность

 

Ñy[n] = y[n] - y[n - 1].

(6.19)

Прямая разность определяется в момент времени t = nT по буду- щему значению решетчатой функции при t = (n + 1)T . Это можно сделать,

если будущее значение известно (числовая последовательность выборок сигнала на каком-то интервале времени записана в память машины), либо если это будущее значение можно вычислить по аналитическим зависимо- стям.

Обратная разность определяется

для момента времени

tn = nT по

прошлому значению решетчатой

функции в момент

времени

tn−1 = (n −1)T .

 

 

Аналогом второй производной непрерывной функции служат вто- рые разности прямая

D2 y[n] = Dy[n + 1] - Dy[n] = y[n + 2] - 2y[n + 1] + y[n],

(6.20)

и обратная

 

Ñ2 y[n] = Ñy[n] - Ñy[n - 1] = y[n] - 2y[n - 1] + y[n - 2].

(6.21)

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

Аналогом интеграла непрерывной функции в пределах выборки от 0 до n для решетчатой функции является неполная сумма

107

n−1

(6.22)

y[n] = y[m].

m=0

иполная сумма

n

(6.23)

y[n] = y[m]

m=0

Отличие выражений (6.22) и (6.23) состоит в том, что значение y[ n] при определении полной суммы также участвует в формировании ре-

зультата.

Динамические непрерывные объекты описываются дифференци- альными уравнениями, дискретные - разностными. Для однозначного ре- шения дифференциального уравнения необходимо знать начальные или конечные условия. Если система описывается уравнениями n-го порядка, необходимо знать n-1 начальных или конечных условий. В большинстве задач анализа систем управления начальные условия задаются в момент времени t=0.

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

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

amÑm y[n]+ am−1Ñm−1 y[n]+ ... + a0 y[n] = x[n],

(6.24)

где х[n] - заданная, а y[n] - искомая решётчатая функции. При x[n] = 0 уравнение становится однородным.

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

Так разностные уравнения первого порядка самого общего вида

y[n] = a1 y[n −1]+ b0 x[n]+ b1x[n −1],

(6.25)

можно реализовать с помощью схемы, изображенной на рис 6.5

108

Блок задержки, на этой схеме, реализованный программным путём, осуществляет задержку на один отсчёт. Аналогично может быт решено разностное уравнение второго порядка (рис 6.8).

y [n ] = − a1 y [n − 1]a 2 y [n − 2 ]+ b0 x [n ]+ b1 x [n − 1]+

+ b 2 x [n − 2 ].

(6.26)

 

 

 

 

 

 

 

 

 

 

b1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задержка

 

 

 

+

 

 

 

 

 

 

 

 

y[n]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x[n]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b0

y[n-1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-a1

 

 

 

 

Задержка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 6.5 Схема решения на цифровой машине разностного

 

 

 

 

 

 

 

 

 

уравнения первого порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y[n-2]

 

 

 

 

 

 

 

Задержка

 

 

 

 

Задержка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b0

 

 

 

 

 

 

 

 

 

-a2

 

 

y[n-1]

 

 

-a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х[n]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y[n]

 

 

 

b1

 

 

 

Х[n-1]

 

b2

 

 

 

 

 

 

Х[n-2]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задержка

 

 

 

Задержка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 6.6 Cхема решения на цифровой машине разностного уравнения второго порядка

Для примера рассмотрим решение разностного уравнения первого порядка.

a1Ñy[n]+ a0 y[n] = b x[n].

(6.27)

 

0

 

В записанном разностном уравнении перейдём к дискретам

 

a1(y[n]− y[n − 1]) + ao y[n] = box[n],

(6.28)

 

y[n]×(a1 + ao ) = -a1y[n -1]+ box[n] и окончательно

 

109

y[n] = −ry[n−1] + qx[n],

 

 

(6.29)

где

r =

 

a1

, q =

 

bo

.

 

 

 

 

+ ao

a1

 

 

a1

 

+ ao

В качестве входной последовательности x[n]возьмём единичный скачок и примем нулевые начальные условия: y[n − 1] = 0 . На (рис. 6.7)

приведена огибающая решетчатой функций полученной в результате ре- шения разностного уравнения.

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

Пусть дано разностное уравнение

aoy[n]+ aoy[n −1]+ ... + ay[n − m] = box[n]+ b1x[n −1]+ ... + bnx[n − 2] (6.30)

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

 

F (z){f

[n kn]}= Z knF (z), в результате получим

 

 

ao + a z

−1 + ... + a mz

−m Y(z) =

bo + b z−1

+ ... + bn z−n

X(z) .

(6.31)

 

1

 

 

 

1

 

 

 

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

Передаточная функция будет иметь вид:

 

b

o

+ b z−1 + ... + b

r

z−n

(6.32)

W(Z) =

 

 

1

 

 

 

 

.

 

a

o

+ a z

−1

+ ... + a

m

z−m

 

 

 

 

1

 

 

 

 

 

Умножив числитель и знаменатель выражения (6.36) на Z m получим выражение передаточной функции с положительными степенями числите- ля и знаменателя:

 

b

o

zm + b zm−1

+ ... + b

zm−n

 

W(Z) =

 

 

 

1

r

 

 

.

(6.33)

 

a

 

zm + a zm−1 + ... + a

 

 

 

o

m

 

 

 

 

 

1

 

 

 

Этой передаточной функции соответствует структурная схема при- веденная на (рис. 6.7).

110

b0

b1

b2

 

bn

 

z−1

 

 

 

z−1

 

 

z−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y[n]

X[n]

 

-a2

-a1

-am-1

-am

−1

z−1

z−1

z

Рисунок 6.7 Схема решения разностного уравнения

6.4 Импульсные и частотные характеристики цифровой системы

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

Линейная система с постоянными параметрами, импульсная харак- теристика которой является последовательностью конечной длины, назы- вается системой с конечной импульсной характеристикой. Пример им- пульсной характеристики { K[n] }, при n=1...64 конечной длины изображена на рис. 6.8.

K[n]

0.8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0.2

1

4

7

10

13

16

19

22

25

28

31

34

37

40

43

46

49

52

55

58

61

-0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.6.8 Импульсная характеристика конечной длины КИХ.

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

В отличие от определения, данного для КИХ систем, системой с бесконечной импульсной характеристикой (БИХ) называют системы, дли- ны импульсной характеристики которой не ограничиваются слева или справа.

111