Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 409

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

292 ГЛ . 7, М ЕТОДЫ И ССЛЕДОВАН ИЯ ТОЧНОСТИ Л И Н Е Й Н Ы Х СИСТЕМ

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

УкИ = У (tk, со) = ф (гео) eiahT*= ¥ (еіаТп) е™кТ* (7.5.3)

(& = 0, + 1, ± 2 , ...) .

Подставляя это выражение в формулы (7.2.39) и (7.2.40) и прини­ мая во внимание, что в данном случае

?t= oo, a*= 0, Я2 = -1^п ,

G() = s£(со),

получим следующие формулы для

корреляционной функции

и дисперсии выходной переменной дискретной линейной системы:

2 я /Т п

ку (0і — ti) = ку {{h1)ТП)= J

^(со )|Ф (ш )ре{“ <л- г>Гп ^ ,

(7.5.4)

о

 

 

2я/Тп

 

 

Du = j

(со) I Ф (гео)

(7.5.5)

о

 

 

Интегрирование в формулах (7.5.4) и (7.5.5) производится вдоль отрезка мнимой оси на плоскости комплексной переменной

s от 0 до точки 2яі/Тп. При переходе к переменной z = е5Тп этот отрезок перейдет в единичную окружность С4 плоскости z. Поэтому

при замене переменных z = ег<оТп формулы (7.5.4) и (7.5.5) при­ мут вид

кѵ(mTn) = ± J а* (z) I ¥ (z) I2 zm-i dz,

(7.5.6)

Ci

 

D y = \ o i ( z ) \ ^ ( z ) \ ^ ,

(7.5.7)

Ci

 

где

а интегрирование производится вдоль единичной окружности. Если перейти в формулах (7.5.4)—(7.5.7) к переменной К =

= ѴІІ = (z — 1)li (z + 1) = (егшТп — i)/i (ешТп -f 1), то, учиты­ вая, что единичной окружности плоскости комплексной перемен­ ной z соответствует вся мнимая ось плоскости переменной ѵ,


 

§ 7.5. О П Р Е Д Е Л Е Н И Е УСТАНОВИВШ ЕЙСЯ ДИ СПЕРСИИ

293

получим

оо

 

 

 

 

 

 

 

 

 

 

 

к у ( т Т а ) =

J '?*W|Q(a)|2(^±g)mT^

=

 

 

—оо

оо

 

тп

 

 

 

 

 

 

 

 

 

= J ** ( V

I Q№ ) I22

(ü)2ft

, (7.5.9)

 

 

— оо

оо

f t = 0

 

 

 

 

 

 

 

 

 

 

D„= {

М ^ І Щ ^ І 2-

^ ,

(7.5.10)

где

 

 

_0°

 

 

 

 

 

 

 

 

<7-5Л1)

Формулы (7.5.9) и (7.5.10) в случае рациональных функций sx (Я) и Q (V) удобны тем, что интегралы в них могут быть вычислены по формуле (7.4.6).

Формулы (7.5.4)—(7.5.7), (7.5.9) и (7.5.10), как и формулы предыдущего параграфа, применимы только к устойчивым дискрет­ ным стационарным линейным системам для моментов времени, достаточно удаленных от момента начала работы системы.

П р и м е р 7.5.1. Найти установившуюся дисперсию выходной пере­ менной стационарной дискретной линейной системы, состоящей из импульс­ ного устройства, вырабатывающего прямоугольные импульсы с периодом повторения Тп, и апериодического звена, если входное возмущение пред­ ставляет собой стационарную случайную функцию, интервал корреляции которой меньше периода повторения импульсов Ти. В этом случае значения входного возмущения, действующие на систему, не коррелированъ! (т. е. вход­ ное случайное возмущение представляет собой импульсный белый шум с периодом повторения Тп). В примере 5.4.1 была найдена передаточная функция рассматриваемой системы:

У (z)= fczt Zl

1 , Zl= e- V r , ѵ = | а -

(7.5.12)

Спектральная плотность стационарной последовательности некоррелирован­ ных импульсов 4 (со) постоянна и равна TnD/2n, где D — дисперсия каждого из этих импульсов. Следовательно, как показывает формула (7.5.8), функция

а* (z) в данном случае постоянпа и равна D/2n. Подставляя это выражение спектральной плотности о% (z) и выражение (7.5.12) передаточной функции

системы

 

в формулу (7.5.7),

находим

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

системы:

 

 

 

_ D k 4 \

у

2 Г ___ dz__

 

 

 

 

 

(7.5.13)

 

 

 

 

 

і

( 1

} J

 

 

 

 

 

 

 

 

Формула

(7.5.5)

в этом

случае дает

Сі

 

 

 

 

 

 

 

 

 

2 П /Т и

 

 

 

 

_ ТиРкЧ\

 

у

{

da>

 

 

 

Dy =

(zrT-l)2

ia > T „

 

 

 

 

 

 

 

-*і I2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dkzz\

dcp

(7.5.14)

 

 

 

 

 

 

 

 

 

 

 

 

2n (zfV- l ) 2 j

1 —2zj cos ф+ 2і

 

 

 

 

 

 

 

Ü


294 гл. 7. М ЕТО Д Ы И ССЛЕДО ВА Н И Я ТО ЧН О С ТИ Л И Н Е Й Н Ы Х СИСТЕМ

К этому же виду приводится формула (7.5.10) заменой переменных z =

=

е

п =

егф. Выполняя интегрирование в

формуле

(7.5.14), полудим

 

 

 

Dy = Dk4\ {Z\1—Zj1)2

 

 

(7.5.15)

Если

воспользоваться формулой (7.5.10), то,

принимая

во внимание, что

в

данном

случае

 

 

 

 

 

Q(v) = v

1 -Г

 

~ .

D

(7.5.16)

 

( r ^ ) = fcz‘ (zr v- 1) iz : *i+ (l + *i)»*

S*()

я

 

 

цолучим

 

 

 

 

 

 

 

 

 

dk

 

 

(7.5.17)

 

 

 

I т т з zi + (l +

zi) ik

I2

 

 

 

 

Применив для вычисления интеграла формулу (7.4.6), снова получим форму­ лу (7.5.15).

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

Из теории случайных функций известно, что для нахождения интегрального канонического представления случайной функции X (t) в данном интервале (/0? h) изменения аргумента t достаточно найти три функции х (t, к), а (t, к), G (к), удовлетворяющие усло­

виям (см. [54], § 6.5 или [53], §

67)

 

«1

 

k')dt = ö(k k’),

 

j x(t,

k)a(t,

(7.6.1)

*0

 

 

 

J x(t,

k)a{t',

k)dk = 8 (t t'),

(7.6.2)

Ki

ti

 

 

 

 

 

г(і,к) = Щ£}<\

Kx (t, t') а (t', к) dt'.

(7.6.3)

 

to

 

 

Если такие функции найдены, то случайная функция

 

 

ц

 

 

V(k)= j

a (t, к) Х° (t) dt

(7.6.4)

 

<0

 

 

будет представлять собой белый шум, интенсивностью которого является функция G (к), а случайная функция X (t) выразится через белый шум V (к) интегральным каноническим представле­ нием

X ( t ) ^ m x ( t) + ^ V ( k )x (t, k)dk (t0< W i ) .

(7.6.5)

ii

 


§ 7.6. П РИ М ЕН Е Н И Е ТЕО РИ И Л И Н Е Й Н Ы Х СИСТЕМ

295

Корреляционная функция белого шума V (А,) определяется фор­ мулой

К в (Я, X') = G (Я) б (Я — Я').

(7.6.6)

В общей теории интегральных канонических представлений параметр Я может быть произвольной величиной и в соответствии с этим V (X) является белым шумом не в физическом, а в обобщен­ ном, математическом смысле. Однако в приложениях часто оказы­ вается возможным рассматривать случайную функцию времени как результат прохождения некоторого физического белого шума, представляющего собой случайную функцию времени, через некоторую линейную систему. В этом случае параметр X пред­ ставляет собой время, а функцию х (t, X) можно рассматривать как весовую функцию системы, формирующей данную случайную функцию времени из белого шума. Функция а (t, Я) в этом случае может рассматриваться как весовая функция обратной системы, преобразующей данную случайную функцию в белый шум. Эти соображения наводят на мысль попытаться применить для нахож­ дения интегральных канонических представлений случайных функций теорию линейных систем. Как будет видно из дальней­ шего, это во многих случаях приводит к положительным резуль­ татам.

В § 4.2 было показано, что весовые функции двух взаимно обратных физически возможных линейных систем в любом интер­ вале (г0> h) связаны соотношениями (4.2.10) и (4.2.11), которые можно переписать в других обозначениях в виде

н

J

н>(*, Я)ыг(Я', *).Л = 6(Я —Я')

(t0<X, X'< h),

(7.6.7)

to

 

 

 

ti

w (t, X) w~ (X, t') dX= 8 (t 1')

(t0*Ct, t' C ^).

 

j

(7.6.8)

to

Здесь мы произвольно расширили пределы интегрирования, имея в виду, что вследствие условия физической возможности подынте­ гральная функция в первой формуле отлична от нуля только при Я < Я ' в интервале (Я, Я'), а подынтегральная функция во второй формуле отлична от нуля только при t > t' в интервале (t', t). Соотношения (7.6.7) и (7.6.8) по существу не отличаются от (7.6.1) и (7.6.2). Поэтому, если известны весовые функции w и w~ двух взаимно обратных линейных систем, то, принимая во внимание, что эти весовые функции действительны, можно удовлетворить условиям (7.6.1) и (7.6.2), приняв

ж (t, X) = w (t, Я), a (t, Я) = W- (Я, t).

(7.6.9)

Таким образом, весовые функции двух взаимно обратных линей­ ных систем всегда удовлетворяют двум условиям (7.6.1) и (7.6.2)


296 Г Л . 7. М ЕТО Д Ы И С С Л Е Д О В А Н И Я ТО ЧН О СТИ Л И Н Е Й Н Ы Х СИСТЕМ

в любом интервале (t0, ti), если у одной из них первый аргумент рассматривать как независимую переменную, а второй — как параметр, а у другой наоборот. Следовательно, для нахождения интегрального канонического представления случайной функции времени X (^’достаточно среди всех известных пар весовых функ­ ций взаимно обратных физически возможных линейных систем найти такую пару, для которой удовлетворяется и условие (7.6.3) в некотором интервале (t0, £t). Критерием выполнения усло­ вия (7.6.3), как это видно из формул (7.6.3) и (7.6.9), может слу­ жить независимость от времени t выражения

w (t, X) I К х

w

dt>‘

to

 

 

Если это выражение при некоторых значениях tQи ^ не зависит от времени t, то, полагая при этих t0 и t\

ti

 

Kx (t,t')w-(X,t')dtr,

(7.6.10)

’ to

мы определим функцию G (Я), удовлетворяющую совместно с функ­

циями (7.6.9) условию (7.6.3).

функция

На основании формул

(7.6.4) и (7.6.9) случайная

 

«1

 

V (Я) =

j иг (Я, t) Х° (0 dt

(7.6.11)

 

*0

 

будет при этом белым шумом параметра Я, интенсивность которого равна G (Я). Случайная функция X (t) выразится через этот белый шум интегральным каноническим представлением

ti

 

X (t) = тх (t) + £ w (t, Я) V (Я) d%

(7.6.12)

to

Корреляционная функция случайной функции X (t) выразится соответствующим интегральным каноническим представлением *)

ti

 

 

Kx(t, t') — j G (Я) w (t, Я) w (t’, Я) d%

(t0*Ct,

(7.6.13)

io

*) Для вывода формулы (7.6.13) достаточно воспользоваться форму­ лой (7.2.6), расширив в ней пределы интегрирования и заменив верхние пределы обоих интегралов величиной U, предполагая, что t, t' < Ц. Тогда, принимая во внимание (7.6.6), получим

11 ti

Их(6 С)—j £ w(t, Я) w{t', X')G(X)6(X — X')dXdX'.

to to

Выполняя здесь интегрирование по X', мы и получим формулу (7.6.13).