![](/user_photo/_userpic.png)
книги из ГПНТБ / Востров М.В. Основы авиационной автоматики
.pdfВторая |
разность |
|
|
А2 |
= А хга+ 1 \ х п ~ ( X ll+ 2 -^/!+l) |
( 4 1 |
■*„) |
и окончательно |
х„. |
(3.110) |
|
|
Д2хл = А-л+2- 2 х л+1 + |
Аналогично можно определить и разность порядка г:
Т
Е ( - 1)*С*ж„+,_*. (3.111) fe-0
Коэффициенты С* являются биномиальными коэффи циентами.
х(Ьтл)
х'' |
|
v |
At=0 |
|
|
/ |
|
|
|
||
|
|
Т ^ -Г |
■птп |
||
|
|
2Тд |
|||
|
|
ЗТЛ |
|
|
|
Т п 1 |
|
\ |
At < 0 |
|
|
ф |
|
|
|
|
|
|
■ |
- J ________J— |
■— пТл |
||
__I___ |
|||||
Тп |
|
2ТП |
ЗТП |
|
|
f(^nM |
|
|
|
|
|
т"гI |
|
At >0 |
|
||
Y.1 -~*tt |
- |
I |
1 |
■пТп |
|
|
I |
: : . ^ |
|||
|
|
Ф и г. 3.55 |
|
|
Конечные разности являются до некоторой степени аналога ми дифференциалов или производных у непрерывных функций.
Аналогами интегралов непрерывных функций являются сум мы дискретных функций:
я—1
(3.112)
<=о
Легко понять, |
что первая |
разность суммы равна самой дис |
|
кретной функции |
|
п |
я —1 |
|
|
||
АКя = К я + 1 ~ > 'я |
= ^ л : | - |
(3.113; |
|
|
|
г=о |
«-о |
При исследовании импульсных систем приходится иметь де ло с уравнениями, связывающими дискретную функцию и ее разности. Такого рода уравнения называются уравнениями в конечных разностях или просто разностными уравнениями.
160
Существуют две формы записи разностных уравнений:
Первая форма:
a r А'г х п + a 'r _ i Д'’-1 х п + ' ' ' + a i ’ ^ х п + а о х п —
= Ьп' Д-Л+ Ь-т_хд»-'Л + ••.■+Ь{ Д/„ + V/» • (3.114)
Если к уравнению (3.114) применить соотношение (3.111), то .мы получили вторую форму разностного уравнения
Хи+г “Ь CLr —1 xn-f-r—i -i------ V a-о Хп ~ b mfn+m ~Ь ‘ ’ ‘ ^0 f п‘ (3.115)
Классический метод решения разностных уравнений во мно гом аналогичен классическому методу решения дифференциаль ных уравнений. 1
Если разностное уравнение — линейное, с постоянными коэф фициентами, то составляется характеристическое уравнение од нородного разностного уравнения, отыскиваются его корни и формируется решение. Далее отыскивается решение неоднород ного уравнения. Для определения неизвестных .коэффициентов учитываются начальные условия.
2, Д и ск р ет н о е п р е о б р а зо в а н и е Л а п л а с а (^ -п р е о б р а зо в а н и е )
Перейдем теперь к изложению сущности дискретного пре образования Лапласа, которое, с одной стороны, позволяет от носительно просто получить решение разностных уравнений, а с другой — является основным математическим аппаратом тео рии импульсных систем.
Дискретное преобразование Лапласа является функциональ ным. преобразованием и определяется соотношением:
о©
Х (р) = S х(п)е-~рпГ* . |
(3.116) |
п~ О |
|
Данное соотношение отличается от известного выражения для непрерывного преобразования Лапласа тем, что непрерыв ное время t заменено дискретным временем пТп и интеграл —
суммой.
Дискретное преобразование обычно применяется .в форме, предложенной профессором Ципкиным Я. 3.:
СО
X (q) = Yi |
*(л)е-*я, |
(3.117) |
7 1 = 0 |
|
|
(где q = рТ„ — комплексное |
число, называемое |
параметром |
преобразования) и в форме так называемого г-преобразования, широко используемого в иностранной литературе:
оо |
|
X (z) = S x { n ) z - \ |
(3.118) |
n=0 |
|
где 2 = ерГп; ДДг) называют при этом 2 -преобразованием дис кретной функции х(п).
11. Изд. № 3912 |
161 |
Дискретное преобразование обозначают иногда и так:
Ьж[х[п)} или 1Л[.«„], |
а также |
Lz [x{n)\ |
и Lz [xn\. |
|||
Соотношения, аналогичные |
(3.117) и (3.118), |
можно написать |
||||
и для смещенных дискретных функций: |
|
|
||||
x(q, |
е )= |
Yix{n, |
s)e |
qn; |
|
|
|
|
я=0 |
|
|
|
|
|
|
со |
|
|
|
|
X{ Z, |
е ) = |
L |
Х ( Я , |
e ) z ~ n . |
(3.119) |
и—О
Приведенные зависимости дают связь между дискретной функцией и изображением и потому называются прямым дис кретным преобразованием Лапласа.
Обратное преобразование представляется в таком виде
Х(п). |
1 |
X (z ) z n~l dz. |
(3.120) |
|
2*У
Здесь с — произвольный замкнутый контур, внутри которо
го лежат все полюса функции X(z).
В качестве примера определим 2 -преобразования некоторых
простейших функций:
1) пусть Х(п) — 1 (п):
* (* ) = E z-" = (l + 2 - 1 + 2 - 2 + ..•)
л=0
2)пусть х(пТп) = пТп
оо
X ( z ) = Тв V nz~n= Тп [ z ~ l + 2z~2-f - - - ) = = 7*п
л=*0
3) пусть ^(яГп)
X (г)= 2
—лпТ„
е _ я л Г п
(2 - 1)2
-аг„
л - 0
Таблицы 2 -преобразований наиболее часто встречающихся
функций обычно приводятся в литературе.
В данном пособии таблицы преобразований даны в прило жении 2.
162
Для дискретного преобразования Лапласа имеют место ана логичные основные правила и теоремы, что и для непрерывно
го преобразования Лапласа ДО]. |
теорема об изображении |
|
Большое значение для нас имеет |
||
смещенных функций (теорема сдвига) |
и теорема об изображе |
|
нии разностей. |
|
|
Пусть нам необходимо найти изображение смещенной функ |
||
ции x ( n - \ - k ) |
или x n+k. |
|
|
во |
|
По определению 1 д \x„+k] = X x n+k z~n. |
||
Сделаем |
П-—б |
|
подстановку п k — г, тогда |
|
L, [xn+k]= |
X |
|
г- |
|
|
|
|
|
г~т |
|
|||||
|
|
|
|
л=*0 |
|
|
|
|
k—1 |
r ~ k |
|
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
x r z~r — X х , z ~ r |
|
|
|
|
|||||||
Окончательно |
|
|
г=О |
|
|
|
|
Г—О |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
k - i |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ln\.xn-\-k\ |
— zk |
^■д |
|
|
X |
XrZ |
|
|
|
(3.121) |
|||||
|
|
|
|
|
|
|
|
|
|
r - 0 |
|
|
|
|
|
|
Если X q = X] = |
Л'2 |
= |
. . . = |
A'a _ ] |
= |
0, |
TO |
|
|
|
|
|
|
|||
|
|
|
|
■A Ая-Гй] ~ |
^ |
Ai Anl- |
|
|
|
|
(3.122) |
|||||
Аналогичное соотношение |
|
можно указать |
и |
|
для функции |
|||||||||||
Для к = I |
|
LA[xn_k\ = z ~ k LA\xn\. |
|
|
|
|
(3.123) |
|||||||||
Ад К +il == |
|
|
(хп) |
х0]. |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||
Полученные |
соотношения дают |
возможность |
написать изо |
|||||||||||||
бражение |
разностей. |
Так для |
кхп~= х п+1— хп |
|
|
|
||||||||||
А Л А * я ] = |
А А я А |
' я ) - |
* о |
] ~ |
А |
л |
( - к |
я ) = - |
( 2 - |
1 ) |
А |
( * |
я ) |
~ |
z x 0. (3.124) |
|
Для разности |
любого |
порядка: |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
т—1 |
|
|
|
|
|
|
|
|
Ад [Ая,^1 |
= (2 - |
1)я Ад(*я) - |
|
г"» |
X |
(z - |
l)" -1-- |
Д‘- |
(0). (3.125) |
|||||||
|
|
|
|
|
|
|
|
|
г-о |
|
|
Д°ап(0)= Д'хп(1)= |
||||
Здесь нужно считать Д°а (0) = а 0;при |
||||||||||||||||
= ... = 0, |
соотношение |
(3.125) |
упрощается |
|
|
|
|
|
||||||||
|
|
L J A mx n] = ( z - |
\)”Ч л [хя). |
|
|
|
(3.126) |
Зная основные теоремы дискретного преобразования Лапла са, мы можем разностные уравнения записать в форме г-пре- образования.
11* |
163 |
Применяя соотношение (3.126) к уравнению (3.114), получим первую форму разностного уравнения в 2-преобразованиях:
[ a / ( z — 1 ) 4 |
l)r_1 + |
----- b a / ( z — 1) + a0'\ |
X{z) = |
||||||
- |
[bm’( z - \ ) m + |
Ът_л (z - |
1 |
+ |
• • • + |
b ' (z - |
1) + b0] F (z). |
||
|
|
|
|
|
|
|
|
|
(3.127) |
|
Аналогично, мы можем соотношение (3.122) |
применить к |
|||||||
уравнению |
(3.115) |
и получить вторую форму разностного урав |
|||||||
нения в 2 -преобразованиях: |
|
|
|
|
|
|
|||
|
|
\arz r + ar_1z r~l -|------- b |
-f |
a0\X (z) = |
|
||||
|
= [bmzm |
------- \-bt z + b0]F(z). |
(3.128) |
||||||
из |
Легко видеть, что уравнение (3.128) |
может |
быть |
получено |
|||||
(3.127) |
непосредственно. |
|
|
|
|
|
|
||
|
3. |
Д и н а м и к а р а зо м к н у т ы х |
и м п ульсн ы х |
си стем |
|
||||
|
П е р е д а т о ч н ы е фу н к ц и и р а з о м к н у т ы х |
||||||||
|
|
|
и м п у л ь с н ы х с и с т е м |
|
|
||||
|
Динамика разомкнутых |
импульсных |
систем |
предполагает |
изучение поведения непрерывной части под воздействием посту пающих на ее вход импульсов.
Введем ряд определений и понятий, облегчающих исследова ние динамики как разомкнутых, так и замкнутых систем.
Так |
как |
разомкнутая |
система может быть |
представлена |
|||
в виде |
последовательного |
соединения |
8-импульсного |
элемента |
|||
и приведенной НЧ, то представляет |
интерес отыскания реакции |
||||||
непрерывной части на последовательность 8-импульсов. |
|||||||
Как известно, реакция |
НЧ на |
8-импульс есть весовая или |
|||||
импульсная |
переходная |
функция |
системы |
Hb{t), |
причем |
Hb(p )~ W { p ) .
Фиг. 3.56
Последовательность дискретных значений весовой функции Hb{t) в моменты времени пТп называется дискретной пере ходной функцией, или дискретной весовой функцией системы.
Дискретная весовая функция записывается так:
Нь{пТа) = д а п(л) = {«/„, w x,• • •, wn}. |
(3.129) |
Числа последовательности w0, wь..., «^называются коэффи циентами веса или весовыми коэффициентами (фиг. 3.56).
164
По аналогии с непрерывными системами
со
LAHb{nTn)]= W (z) = Е тпг~*. |
(3.130) |
п —0 |
|
z-преобразование дискретной весовой функции называется дис
кретной передаточной функцией |
разомкнутой |
системы и обо |
||||||||||
значается W(z). |
|
|
|
|
|
|
|
|
|
|
|
|
Определим передаточные функции для простейших элемен |
||||||||||||
тарных звеньев: |
|
|
|
|
|
|
, |
|
|
|
|
|
а) Интегрирующее звено |
W(p) — S - |
|
|
|
||||||||
Для |
интегрирующего звена |
|
|
Р |
k- |
1 ( / ) , |
следовательно, |
|||||
Н ь ( t ) = |
||||||||||||
Ht (nTa) = wn(n) - {k, k, •••}. |
|
|
|
|
|
|
|
|||||
Отсюда |
|
|
|
оо |
|
|
|
|
|
|
|
|
|
Г (г ) = |
k Y . z - n = |
— |
|
. |
|
(3.131) |
|||||
|
|
|
|
|
|
|
|
г — 1 |
|
|
|
|
б) |
Инерционное звено W (р) = |
|
у |
• |
|
|
||||||
Для |
этого звена |
|
|
_ |
Ь _ |
- |
~ |
и |
значит |
|
|
|
Я 8(^)> |
Т |
|
|
|
|
|||||||
|
|
|
|
^ |
с |
|
|
|
|
|
|
|
|
|
|
лГг |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
у |
( 2 т ) " |
(где |
г т = |
е |
|
||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
\^ (г ) |
|
■ У 7 2' Т - |
k |
|
Z |
|
(3.132) |
||||
|
|
Г |
|
Z r |
|
|||||||
|
|
|
|
Ц |
г |
) |
|
2 |
|
|
||
|
|
|
|
л=»0 |
|
|
|
|
|
|
|
|
Реакция НЧ на последовательность |
S-функций определяется |
|||||||||||
на основании соотношения, называемого сверткой |
двух дис |
|||||||||||
кретных функций ['1]. |
|
|
|
|
|
|
|
|
|
|
||
Сверткой двух дискретных функций называется следующая |
||||||||||||
зависимость: |
|
п |
|
|
|
|
|
|
|
|
||
|
х [пТа] = |
£ |
w [(/г - |
m) Тп] g [m Тп\, |
(3.133) |
|||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
при этом Щ[(п—tn) |
7'n]g'[m7'n] = |
0 для |
всех (п—т) |
< 0 . |
||||||||
Выражение для |
свертки двух |
дискретных функций можно |
||||||||||
записать и так: |
|
П |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
|
л:[лГп] = |
Е |
w [тТп\ g {(п — т) Тп\. |
(3.134) |
||||||||
|
|
|
т —0 |
|
|
|
|
|
|
|
|
|
Можно дать и другое определение дискретной передаточной |
||||||||||||
функции разомкнутой |
системы. |
|
|
|
|
|
|
|
165
Пусть |
на |
систему (фиг. 3.57) действует последовательность |
|||
g ['и?'.,]- |
Выходная величина системы х \п Т п\. |
|
|||
Тогда, |
следовательно, (3.133) можно записать |
|
|||
|
|
*о = |
®’о £ 0; |
|
|
|
|
*1 = |
w o g i + |
щ go; |
|
|
|
*2 = |
Щ ^2 + |
Wl g 1 + ^2 £<) ? |
|
|
|
*3 = W0 gi + |
£2 + w2 gi + W3 g0. |
|
|
Умножая каждую строку соответственно на 1, z !, z' |
2, и т. д. |
||||
и суммируя, |
после преобразования получим |
|
|||
Отсюда |
|
X(z) = W(z)G{z) . |
(3.135) |
||
|
|
|
|
||
|
|
|
W(z) = ^ - ^ - . |
(3.136) |
О (г)
Дискретная передаточная функция разомкнутой системы равна отношению z-преобразований выходной и входной ве личин.
Дискретная передаточная функция существенно зависит не только от структуры непрерывной части, но и от типа импульс ного элемента, точнее от фор мы импульса, формируемого
х |
им. |
В эт ом с л у ч а е |
W (г ) о п р е |
|
|
||
|
д е л я е т с я т а к ж е из со о т н о ш е н и я |
||
|
(3.130), т ол ь к о w n |
о п р е д е л я е т |
|
ся и г. 3.57 |
ся |
д л я п р и в е д е н н о й н е п р е р ы в |
|
|
ной |
части . |
|
Для примера определим W (г) интегрирующего и интерцион-
ного звеньев для случая, когда ИЭ формирует прямоугольные импульсы.
а) Интегрирующее звено W (р) = |
k |
|
|
|||||
— . |
|
|
||||||
Передаточная |
функция |
приведенной НЧ |
|
|
||||
W 0(P) = |
( 1 |
- е |
-рт, |
, 1 |
k |
_ |
|
Q~PTa |
п -------- |
Р2 |
|||||||
|
|
|
|
Р |
Р |
Р1 |
|
|
Находим Нг (t) — kt — k (t — Та) = kTn . |
|
|
||||||
Следовательно, |
wn = k T„ |
и значит |
|
|
||||
|
W(z) = kTn Y z - " |
= - ^ ^ |
. |
(3.137) |
||||
|
|
|
|
/1=1 |
2 — |
1 |
|
|
|
|
|
|
|
|
|
166
Суммирование здесь ведется от п ~ |
1, |
поскольку W0 —■0. |
|
||||||||||
б) Инерционное звено W (р) = |
|
|
k |
|
|
|
|
|
|
||||
Т р + 1 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
|||
Передаточная функция |
W 0[p) = (\ — е |
|
р п) |
Тр + |
■, или |
||||||||
|
|
|
|
kT |
|
|
|
|
|
Р |
\ |
||
W 0(p) = |
(1 — е~рТп) |
----- |
|
|
|
|
|
|
|
|
|
|
|
Тр + |
1 , |
|
|
|
|
|
|
|
|
||||
Как |
,Р |
|
|
|
|
Hb{t) |
и |
далее |
|||||
и в предыдущем случае, |
определяем |
||||||||||||
{пТп) = wn — k{\ — zT) z nT~l |
|
( я = |
1, 2, |
|
•••)• |
|
|
|
|
||||
Окончательно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W { z ) = |
k{l ~ |
Zt) . |
|
|
|
|
|
(3.138) |
||||
Можно для определения W (г) |
использовать |
уже |
готовые |
||||||||||
таблицы изображений |
(приложение 2). |
|
|
|
|
|
|
||||||
Так, например, для |
случая а) запишем W0(p) так: |
|
|
||||||||||
|
W 0(p) = z - |
1 |
1 |
|
|
|
|
|
|
||||
|
|
|
|
|
2 |
|
Р2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Далее по таблице находим--------: -* kTB |
2 |
|
|
|
|
||||||||
|
|
|
|
|
|||||||||
Окончательно |
|
|
Р2 |
|
|
( г - 1 ) 2 * |
|
|
|||||
|
|
|
|
|
|
|
кТп |
|
|
|
|||
|
W(z) = - — - |
k T n — -— |
|
|
|
|
|||||||
|
|
z - |
1 |
‘ |
|
|
|||||||
|
|
2 |
|
( 2 - |
1)2 |
|
|
|
|||||
Заметим одну важную особенность определения дискретных |
|||||||||||||
передаточных функций |
разомкнутой системы. Если |
|
|
||||||||||
|
W{p) = W 1(p)W2(p), |
|
|
|
|
|
|
||||||
то W(z) ф Wi (z) W2 (z) |
и |
т. |
д. |
|
|
|
|
|
|
|
|
Последнее равенство удовлетворяется только в том случае, если звенья разделены импульсными элементами, работающими синхронно (фиг. 3.58).
Фиг. 3.58
Если непрерывная часть импульсной системы не является элементарным звеном, а описывается передаточной функцией вида
w ( j = bmpm+ B m^ p m~4-------- |
h&ii + bp ^ |
anpn + an_iP*-1H------ |
{-axp + a0 |
167
то для нахождения дискретной передаточной функции необходи мо W (р) представить в виде суммы элементарных дробей
П
< 3 1 3 |
9 > |
<=1 |
|
Когда все p t действительны, W (р) можно представить в ви |
|
де параллельного соединения инерционных звеньев. При |
p t |
комплексных отдельные звенья будут колебательными. Дискрет ная передаточная функция каждого элементарного звена мо жет быть взята из таблицы (приложение 3).
Дискретная передаточная функция параллельного соедине ния звеньев определяется как сумма дискретных передаточных функций звеньев.
Ч а с т о т н ы е х а р а к т е р и с т и к и р а з о м к н у т ы х и м п у л ь е н ы х с и с т е м
Так же как для непрерывных систем, частотные характери стики импульсных систем характеризуют установившуюся реак цию на гармонический входной сигнал.
Гармонический входной сигнал будем далее брать в комп лексной форме
|
i l ^ n j = e ;’ra,!7% |
(3.140) |
или, обозначая |
u>T„ = v, |
|
|
^[п] = Д'1П. |
(3.141) |
Изображение |
входного гармонического |
сигнала |
|
LA g («)] = — Z— r - |
(3.142) |
|
г — e ' v |
|
Тогда на основании (3.135) изображение вынужденной со ставляющей выходной величины
X(z) = W ( z ) ---- z---- . |
(3.143) |
г — e-iv |
|
Откуда, сама вынужденная составляющая |
|
Щ п ) = W { d ' ) T A |
(3.144) |
где W — W (г) \z—txы — дискретная амплитудно-фазовая характеристика звена или системы.
Далее АФХ дискретной разомкнутой системы будем обозна чать просто W (уу).
168
АФХ мы можем представить в показательной и алгебраичеекой формах
(3.145)
W (yv) = и (v) + j V (v).
В отличие от непрерывных систем, дискретная АФХ является
периодической функцией частоты v. Период v= 2^^==
\ тп
Вследствие симметрии годографа относительно вещественной оси, его вычеркивают для частот v от 0 до тг.
Рассмотрим частотные характеристики простейших звеньев: а) Интегрирующее звено с 8 -импульсным элементом:
kz |
или |
|
|
kepTn |
|
W(z) = |
W (р) |
i |
|||
z - 1 |
|
|
|
ер7п _ |
|
.Отсюда, амплитудно-фазовая характеристика |
|
||||
W{jv) = |
|
|
ft(cosv+ /sln v) |
(3.146) |
|
|
|
|
|
||
е'1— 1 |
(cos v — 1) -j-j sin v |
||||
После преобразования, |
получаем |
|
|||
,| F . . . . |
k |
y |
. |
k , V |
(3.147) |
i r O v ) = Y - |
|
y ctg |
б) Интегрирующее звено с ИЭ, формирующим прямоуголь ные импульсы:
W (z) |
kl' |
|
W(p) |
|
kT„ |
|
|
||
|
1 |
|
|
|
|
||||
Z |
|
|
|
|
Q |
P l П |
|
|
|
Значит АФХ |
|
|
|
|
|
|
|
|
|
W(p) |
|
|
kTn |
|
|
|
|
||
|
(cosv — 1) + /sin |
|
|
||||||
|
|
|
|
||||||
Окончательно |
|
|
|
|
|
|
|
|
|
WU' ) = |
kTn |
. |
j |
kT„ |
, |
v |
|
(3.148) |
|
---- 5 |
|
2 |
—ctg — |
|
|||||
|
|
2 |
|
|
|
2 |
|
|
|
Дискретные АФХ для интегрирующего звена |
показаны на |
||||||||
фиг. (3.59,а и б). |
|
|
|
|
|
|
|
|
|
Следует отметить |
|
характерное |
положение |
этих |
двух |
||||
характеристик относительно |
АФХ |
|
непрерывного |
звена |
W ( n = * - j -
169