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

книги / Теория автоматического управления. Дискретные системы

.pdf
Скачиваний:
3
Добавлен:
19.11.2023
Размер:
72.74 Mб
Скачать

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

цифровой системой.

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

Рис. 3.4. Структурная схема одноконтурной цифровой системы

Упрощенная структура разомкнутой цифровой системы будет иметь вид, приведенный на рис. 3.5. Цифровое управляющее устройство (ЦУУ), реализующее операции с цифровыми кодами, описывается передаточной функцией D(z). Линеаризованный коэффициент передачи ЦУУ совместно с АЦП и ЦАП определяется как

kЦ= kЭВМ kА kЦ = kЭВМ δЦ / δA.

Рис. 3.5. Структурная схема разомкнутой цифровой системы

Операторное уравнение разомкнутой системы в соответствии с рис. 3.5 будет иметь вид

Y (p)=WO (p)WЭ (p)kЦ X * (p)=WO (p) 1еpTp kЦ X * (p), (3.11)

применяя к которому Z-преобразование Лапласа, получим

41

 

 

 

Tp

 

 

 

 

 

 

Y (z)=Z

1e

 

W

(p)

 

k

 

X (z)=

 

 

 

 

Ц

 

p

O

 

 

 

 

 

 

 

 

 

 

(3.12)

 

 

 

 

 

 

 

 

 

=Z 1eTp WO (p) kЦkАD(z)E(z).

p

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

Wраз(z) = kЦ kАD(z)WОWЭ(z).

(3.13)

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

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

 

m

 

W (p)=

k(1+τj p)

 

j=1

.

n

 

p(1+Ti p)

 

 

i=1

 

Разлагая исходную передаточную функцию на простые дроби:

 

k

n

ci

 

 

W (p)=

+k

 

,

p

1+T

p

 

 

i=1

i

 

 

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

 

kTz

n

c z

 

T

 

T

W (z)=

 

 

+k

i

, di =e

i .

z 1

z di

 

i=1

 

 

42

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

Wз(z)=Wраз(z)/(1+ Wраз(z)).

(3.14)

Наряду с непрерывными сигналами задания и обратной связи (см. рис. 3.4) широко используется базовая структурная схема цифровых систем управления, приведенная на рис. 3.6, описываемая системой уравнений:

E* ( p) =V * ( p) kAY* ( p);

 

 

 

 

 

 

 

 

 

Y ( p) =k

Ц

W W * ( p) X * ( p);

(3.15)

 

 

 

 

Э

O

 

 

*

 

 

*

 

*

 

 

( p) =E

 

 

X

 

 

( p) WM ( p).

 

Рис. 3.6. Базовая структурная схема цифровой САУ

В уравнениях (3.15)–(3.20) передаточная функция ЦУУ обозначена WМ* (р).

Из системы уравнений получим изображения выходного сигнала Y*(p) и рассогласования E*(p) относительно изображения задания V*(p):

Y* ( p) =k

Ц

W W *

( p) W*

( p)[V * ( p)k Y* ( p)];

 

 

Э O

 

 

 

М

 

 

A

 

E* ( p) =V * ( p)k

A

k

Ц

W W* ( p) W*

( p) E* ( p).

(3.16)

 

 

 

 

Э

O

М

 

 

 

 

 

 

 

 

 

 

 

 

43

Из (3.16) находим дискретные передаточные функции системы в формах D- и Z-преобразования Лапласа по отношению к рассогласованию (ошибке):

W *

(p)=

E* (p)

 

=

 

 

1

 

 

 

 

 

 

;

 

 

 

1+k k W W *

(p)W

 

 

(p)

 

зе

 

V * (p)

 

*

 

 

 

 

 

 

 

 

 

 

 

A Ц

Э O

 

M

 

 

 

 

W

(z)=

 

 

 

 

1

 

 

,

 

 

 

 

1+k

k W W (z)W

(z)

 

 

 

 

зе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

Ц Э

O

M

 

 

 

 

 

и выходной переменной:

 

 

 

 

 

 

 

 

 

 

 

 

 

W

(p)=

Y * (p)

=

 

kA kЦWЭWO* (p)WM* (p)

;

 

 

1+k k W W

(p)W

 

 

(z)

зv

 

 

V * (p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A Ц

Э O

 

M

 

 

 

 

 

(z)=

k

k W W

(z)W

 

(z)

 

 

 

 

 

 

W

 

 

A

Ц Э O

M

 

 

.

 

 

 

 

1+k

k W W (z)W

(z)

 

 

 

 

 

зv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

Ц Э

O

M

 

 

 

 

 

(3.17)

(3.18)

Определим дискретную передаточную функцию системы по каналу возмущения, воздействующего на объект. Пусть на объект, состоящий из двух последовательно соединенных непрерывных частей, действует возмущение f(t), показанное на структурной схеме (рис. 3.7). Тогда ошибку, вызванную возмущением, можно определить по формуле

E( p) =−F( p)k W W ( p)W*

( p)k

Y* ( p).

(3.19)

Ц Э 1

М

A

 

 

Рис. 3.7. Структура системы с разделением объекта

44

Отсюда дискретное изображение выходной переменной

Y* ( p) =−FW* ( p)k k W WW* ( p)W*

( p)Y* ( p),

 

2

A

Ц Э 1 2

М

 

 

 

и в форме Z-преобразования

 

 

 

 

 

Y(z) =−FW (z)k k W WW (z)W (z)Y(z),

 

2

A

Ц Э 1 2

M

 

 

 

Y (z)=

 

FW (z)

 

 

.

(3.20)

 

2

 

 

1+k k W W W (z)W (z)

 

 

 

 

A Ц

Э 1 2

M

 

 

 

Параметры входного воздействия входят в (3.20) неявно, поэтому

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

но в подразд. 3.1), выход W1(p) должен быть функцией времени во все моменты времени, а не только вмоменты квантования.

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

АЦП. На рис. 3.7 это воздействие показано как f1(t)=W2(p)f(t) вместо действующего f(t).

Большинство реальных микропроцессорных систем управления являются аналого-цифровыми, поскольку приложенные к ним воздействия и ошибка управления соответствуют непрерывным функциям времени вследствие применения непрерывных фильтров

взадающих устройствах и датчиках выходного сигнала. Базовая структурная схема таких систем приведена на рис. 3.8.

Рис. 3.8. Структурная схема аналого-цифровой системы

45

Данная система, содержащая непрерывное динамическое звено с передаточной функцией WД(p) на входе, в разомкнутом состоянии соответствует схеме на рис. 3.3 и не имеет передаточной функции. Для определения передаточной функции такой системы структурную схему преобразуют к виду, приведенному на рис. 3.9.

Рис. 3.9. Преобразованная структурная схема аналого-цифровой системы

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

 

 

k

D(z)Z W

(p)W

(p)

 

 

W

(z)=W (z)

 

Ц

{ Э

НЧ

}

 

. (3.21)

1+k D(z)Z W (p)W (p) W (z)

p

Д

 

 

 

{ Э

 

}

 

 

 

 

Ц

 

НЧ

Д

 

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

e(t)=eV(t)+eν(t)+eνA(t)+eνЦ(t),

(3.22)

используемой в статистических методах исследований систем.

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

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

Сравнение структурных схем импульсной (см. рис. 3.2) и цифровой (см. рис. 3.5) систем управления в разомкнутом состоянии показывает разницу дискретных передаточных функций, состоя-

46

щую в том, что в уравнение (3.23) импульсной системы для цифровой системы добавлена передаточная функция микроЭВМ D(z), которая не вносит изменений в порядок построения частотных характеристик. Кроме того, в уравнение непрерывной части введено суммарное запаздывание, присущее импульсным элементам второго рода относительно элементов первого рода цифровых систем, время для обработки информации в ЦВМ, наличие звеньев чистого запаздывания в объекте и т.п.:

 

 

 

W

 

(z)

=(1z

1 )Z W

 

(p) ;

 

 

 

 

(3.23)

 

 

 

 

p

 

 

 

 

 

{

ПНЧ

 

}

 

 

 

 

 

 

(

 

)

 

δЦ

(

 

1 )

 

 

WO (p)

 

 

pτ

(

 

)

 

 

Wp

z

=

 

 

Z

 

 

e

D

z

.

(3.24)

 

 

 

 

1z

 

p

 

 

 

 

 

 

 

δA

 

 

 

 

 

 

 

 

 

 

 

 

Это позволяет проводить исследование дискретных систем управления по единым методам анализа и синтеза.

3.3. Передаточная функция цифрового управляющего устройства

Передаточная функция ЦУУ в соответствии с рис. 3.4 определяется как:

D(z)=

X (z)

=

b +b

z +...+b zs

 

 

 

s

s1

0

.

(3.25)

E(z)

a

+a

z +...+a zk

 

 

 

 

 

 

 

k

k1

0

 

 

Для физически реализуемых законов управления степень числителя W(z) должна быть не выше степени знаменателя, иначе для вычисления x[k] требуется знать будущие значения последовательности e[k]. Например, пусть x[k]=b1e[k+1]+b0e[k]. Такой закон управления нереализуем на практике, поскольку значение e[k+1] неизвестно в момент времени t0=kT.

Разделив числитель и знаменатель выражения (3.23) на zk, получим для предельного случая s=k:

47

D(z)=

X (z)

=

b +b z1

+...+b zk

 

 

0

1

k

.

E(z)

a

+a z1

+...+a zk

 

 

 

 

 

 

0

1

k

 

Представив уравнение (3.26) в виде

(a0+ a1z–1+…+ akzk)X(z)=(b0+ b1z–1+…+ bkz–k)E(z)

(3.26)

(3.27)

и перейдя к оригиналам слагаемых с учетом теоремы смещения на целое число тактов, получим разностное уравнение, связывающее решетчатые последовательности x[n] и e[n]:

a0 x[n]+a1x[n –1]+…+ak x[n k]=

(3.28)

=b0e[n] +b1e[n –1]+…+bk e[n k].

Из уравнения (3.26) получим алгоритм работы ЦВМ при a0=1:

x[n]= b0e[n]+ b1e[n –1]+…+ bk e[n k]

(3.29)

{a1x[n –1]+…+ ak x[n k]}.

В соответствии с (3.29) текущее значение выходного сигнала x[n] определяется с использованием текущего значения входного сигнала e[n] и предыдущих значений выходного и входного сигналов. Таким образом, ЦВМ соответствует некоторому фильтру с передаточной функцией D(z), свойства которого определяются видом передаточной функции. Передаточную функцию фильтра находят при синтезе цифровой системы с учетом требований по точности и динамическим свойствам.

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

Передаточная функция имеет вид П-регулятора

DП(p) = kП.

(3.30)

И-регулятор. Численное интегрирование выходной функции u(t) при нулевых начальных условиях по методу Эйлера:

48

du

u

=

u[n]u[n 1]

 

=kИe[n].

(3.31)

dt

t

T

 

 

 

 

Из уравнения (3.31) получим зависимость изображений выходной U(z) и входной E(z) переменных в форме Z-преобразования Лапласа: U(z) – z–1U(z) = kИTE(z), и передаточную функцию вида

D

(z)=

U (z)

=

k

И

Tz

.

(3.32)

E(z)

 

 

 

И

 

z 1

 

 

 

 

 

 

 

D-регулятор. При численном дифференцировании функции входного сигнала e(t) при нулевых начальных условиях методом

простой разности выходная функция будет при u(t)=kД dedt :

u[n]=k

e =k

 

e[n]e[n 1]

,

(3.33)

 

 

 

T

 

Д t

Д

 

 

 

откуда передаточная функция в форме Z-преобразования Лапласа:

D (z)=

kД

 

z 1

.

 

(3.34)

 

 

 

 

 

 

Д

T

 

z

 

 

 

 

 

 

 

 

 

49

4. ЧАСТОТНЫЕ СВОЙСТВА ДИСКРЕТНЫХ СИСТЕМ УПРАВЛЕНИЯ

4.1. Частотные характеристики идеального импульсного элемента

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

На входе элемента – непрерывный сигнал x(t), на выходе – дискретная последовательность, описываемая импульсной функцией x*(t), равной сумме модулированных и смещенных единичных импульсов.

Импульсная функция, представляющая собой периодическую функцию с интервалом T на всей оси времени ∞< t <∞, может быть представлена рядом Фурье в комплексной форме:

 

δ(t nT )= Cr e jrωИt ,

(4.1)

n=−∞

r=−∞

 

где коэффициенты разложения Cr =1/T для любого r, поэтому импульсная функция с учетом (2.4) может быть записана в виде

 

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

(4.2)

n=0

n=−∞

 

поскольку x(t) = 0 при t < 0 и δ(t nT) = 0 при t nT. После подстановки (4.1) получим

 

1

 

x* (t)= x(t)

e jrωиt .

(4.3)

 

 

T r=−∞

 

50