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

книги из ГПНТБ / Востров М.В. Основы авиационной автоматики

.pdf
Скачиваний:
13
Добавлен:
29.10.2023
Размер:
7.14 Mб
Скачать

Вторая

разность

 

 

А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

Соседние файлы в папке книги из ГПНТБ