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

книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов

.pdf
Скачиваний:
71
Добавлен:
25.10.2023
Размер:
16.48 Mб
Скачать

178

Глава 5

В частном случае

равенства у?2(/) нулю функции S i2(f) и S2l(fft

также равны нулю, и формулы (5.50) и (5.51) обращаются в обыч­ ные соотношения

Я і(/)Н Ы т Ь

(5-52)

я 2( / ) Ң Ц ^ - .

(5-53)

Р и с . 5.7. Линейная система с двумя входами.

Случай уЬ(Л = 1 следует рассмотреть особо. Равенство еди­ нице функции когерентности между процессами лу(/) и х2(() озна­ чает полную линейную зависимость. Следовательно, нужно рас­ сматривать уже несколько иную линейную систему, как это по­ казано на рис. 5.8. В этом случае выходной процесс у(() образует­ ся только одним входным процессом xx(t), но двумя различными

Р и с . 5.8. Пример полностью когерентных входных процессов.

способами. Поэтому процессы хх{1) и у(і) связаны только одной частотной характеристикой

H(f) = Hx{f)+H2(f)H3{f).

В общем случае, когда у\2if) ф 0 или 1, знаменатель в формуле (5.50) представляет собой спектр остаточного процесса (Дл:^/)), равного разности процесса jjc^/)} и линейного м. и. к.-прогнала процесса |хі(0} по (х2(0}- Аналогично числитель в формуле (5.50) есть взаимный спектр остаточных процессов (At/(/)] и. (ДХ](0), где (Ау(/)) — остаточный процесс, равный разности процесса \у(і)} и линейного м. н. к.-прогноза процесса {«/(Of"

Процессы на входе и выходе физических систем

179

ПРО (х2(0 ). Эти утверждения базируются на аналогии с основными "понятиями регрессионного анализа, рассмотренными ранее в гл. 4, а также на соображениях, изложенных в следующем разделе.

В общем случае системы с двумя коррелированными входами функции Sbß ) и S2y(f) определяются формулами (5.49), а функ­ ция S uu(f) — формулой (5.42). В частности,

SyV( /) = I Нг(/) I *su (/)+m (f)H2(/)su (/)+

+ m

(f ) H ! ( f ) S n (/) +

I # 2 (/) IZS 22 (/).

(5.54)

Функции обычной

когерентности

выхода с каждым

из входов

^шеют вид

 

 

 

 

..а / п

I H i ( f ) S l l ( f ) + H 2 ( f ) S l t ( f ) \ *

 

У'«и>—

SuUiSyyU)

 

Ѵ2 /ГЧ_ ■ | g i

( » S a (f) +

( f ) S W (Л [ 8

(5.55)

 

 

 

 

S 2 2 ( f ) S y y ( f )

5.4. Функции частной и множественной когерентности

Определения остаточных случайных величин и функций част­ ной когерентности, базирующиеся на критерии среднего квадрата ошибки, приведены в подразд. 5.4.1 и 5.4.2 для частного случая линейной системы с двумя входными процессами. Обобщения этих понятий основываются на том же критерии. Общие определе­ ния функций множественной когерентности даны в^подразд. 5.4.3

и5.4.4.

5.4.1.Остаточные случайные величины

Рассмотрим два действительных стационарных случайных про­ цесса (х(/)) и [y(t)). При фиксированном t величины x(t) и y(t) являются случайными в рассматриваемом выборочном простран­ стве. Для того чтобы избавиться от необходимости центрировать эти величины, положим, что средние значения равны нулю. Ли­

нейный прогноз y(t) величины y(t) по x{t) может

быть запи­

сан в виде

 

у ( 0 = I' h0 {х)х (t— x)dx.

(5.56)

У

 

При необходимости положить равным —оо. весовая функция. При

нижний предел интегрирования можно Здесь h0(x) есть подлежащая определению прогнозировании по методу наименьших

12*

180 Глава 5

квадратов, описанному в подразд. 4.8.2, функция h0(х) выби^

рается таким образом, чтобы средний квадрат ошибки

"

е 2= М [ {

* / ( / ) - £

( / ) П = М

(x)x(t— т)гіт)

 

 

со

со

 

 

=Яуу (0) - 2

J h0 m

xv(x)dx + j j

h0( T )h0 (v)Rxx ( T - v ) d v i T

(5.57)

принимал минимальное значение по всей совокупности возмож­ ных функций /г0(т). Это условие удовлетворяется, если функция h0(x) обладает свойством

де2

0 .

(5.58)

dh0

(т)

 

 

Отсюда следует свертка

 

 

 

СО

 

 

 

W = j К

т х х ( T - v ) d v .

(5.59)

о

 

 

 

Выполнив преобразование Фурье обеих частей этого равенства, получим эквивалентное соотношение

S xy(j) = H 0(f)Sx x (f).

(5.60)

Следует отметить сходство соотношений (5.60) и (5.48).

.-if

Остаточная (условная) случайная величина Ду(І), которую получают в результате вычитания из величины у(() м. н. к.-прог­ ноза величины y{t) по x{t), определяется в виде

СО

 

дг/(0 = «/(0—y(t)=y(t) — j h0(x)x(t~x)dx.

(5.61)

о

 

Ее автокорреляционная функция

 

Я&уАу (т) = М[Д г/(/)Дг/(/ + т)].

(5.62)

Отсюда после вычитания и упрощений получаем

 

СО

 

Яауьу (т) = Ryy (г)—j h0 (v)Rux(т—v)dv.

(5.63)

о

jjè*

Преобразование Фурье последнего уравнения дает равенство

Säyby ( f ) = S y y ( f ) - H 0(f)Syx(f).

(5.64)

Процессы на входе и выходе физических систем

181

Используя формулу (5.60), перепишем этб равенство в виде

S w y ( f ) = S y y ( n - Sxyl^x y{i f — = Syy(f) П ІЛ /)]. (5-ß5)

где

I Sxy ( / ) I

2

 

(5.66)

Уху ( / )

(/)

S X X { l ) S y y

 

есть функция обычной когерентности,

 

связывающая процессы

х(0 и у{1). Величина 5Д4,4г/ (/) в формуле (5.65) называется остаточ­ ной (условной) спектральной плотностью процесса Ay(t)1). Для функции SAyby(f) принято обозначение

Sbyby(n=Sy„.x (f),

(5.67)

указывающее на ее связь с м. н. к.-прогнозом величины у(() по x(t).

Рассмотрим снова случай линейной системы с двумя входами (рис. 5.7), причем входы хх(/) и x2(t) могут быть коррелированьи

Пусть у{1) —- линейный м. н.

к.-прогноз процесса у(і)

по x2(t)

оо

 

 

 

у (t)— j* с (х)х2

(t%)dx.

(5.68)

о

 

 

 

Аналогично x^t) — линейный

м.

н. к. -прогноз процесса Хі(і)

по х2(і)

 

 

 

ОО

 

 

 

хх (t) ^ Ь (х)х2

(t —т)dx.

(5.69)

о

 

 

 

Весовые функции с(т) и Ь{т) определяются подобно тому, как это сделано в соотношениях (5.56) — (5.60). Таким образом, они должны удовлетворять равенствам

S*„(/)=C(/)Si2(/).

(5.70)

Sn(f)=B(f)S22(f).

(5.71)

Определим теперь остаточные случайные величины

 

Ax1(t) = x1(t)— x1(t),

(5.72)

Ay (t) = y(t)—y((),

(5.73)

11 Правильнее называть эту величину остаточной (условной) спектраль­ ной плотностью процесса у(і) либо спектральной плотностью остаточного

»роцесса Ду{і).— Прим, перев.

182

Глава 5

вычитая из хг(/) и y(t) их м. и. к.-прогнозы по х2((). Соответствую­ щие этим величинам остаточные спектральные плотности обозна­ чаются в дальнейшем как

S ЬхіЬхі (/) = S n .2 (/)1

S&yby (/) = Syy.2 (/)>

SAxiby (/) = Sly-2 (/).

Из формулы (5.65) следует, что

S u .2 (/) = S U (/) [1 vf2 (/)],

S y y . 2 ( f ) = Syy( f ) [ l - y yAf)h3

где

712 W

\Si,(f) Г

SuMSnUy

r 2

1S y2 (f ) 12

Ъ А П

S y y ( f ) S 2 2 ( f ) -

1(5.74)

(5.75)

(5.76)

(5.77)

(5.78)

Вывод формулы для функции Slj/-2 (/) несколько сложнее. Согласно уравнениям (5.72) и (5.73), остаточная (условная)

взаимная корреляционная функция1) величин Ахх(1) и Дг/( 0 есть

Riy.t (т) = м lA*i (t)by (t + Т)І-

(5.79)

После подстановки и упрощений получаем

 

00

 

R i y . 2 W = R i y W - f C(ѵ)Яг2 (x—v)dv.

(5.80>7

b

 

Преобразование Фурье обеих частей последнего уравнения дает равенство

Sly.2(f)=Sly( n - C ( n S 12(n,

(5.81)

которое с учетом равенства (5.70) определяет остаточную (услов­ ную) взаимную спектральную плотность1) процессов ДXi(t) и

Ду(0

$1у.2

(/) — Siy (/)

1

__Sl2 ( f ) S 2y ( f )

(5.82)

 

■S22 U )$ ly (f )

 

 

 

 

Сопоставляя уравнения (5.75) и (5.82) с определением частотной характеристики Ях(/) [формула (5.50)1, получаем интересный результат

Яі(/) =

Siy.2 (f)

( 5 . 8 3 )

■Sil.2 (/)

 

 

J) См. примечание на стр. 181.— Прим, перев.

 

Процессы на входе и выходе физических систем

183

5.4.2.

Функции частной когерентности

 

 

Функция

частной когерентности,

связывающая

процессы

x:(t)

и у{1)

при условии, что вход х2(/)

исключен из хг{{) и у(і)

при

помощи линейного м. н. к.-прогноза, определяется как функ­

ция

обычной когерентности между процессами Длу^) и Дy(t):

 

УІул(П= ■ S ii - 2 ( / )Syy.i (!)

(5.84)

Функция частной когерентности удовлетворяет обычному нера­ венству

о < Т?у.2 (/) < 1.

(5.85)

Формулы (5.84) и (5.85) представляют собой прямые обобще­ ния формул (5.23) и (5.24). На основе приведенных ранее соотно­ шений можно показать, что в идеальном случае в предположении, что HX(J) Ф 0 и y\zif) Ф 1, справедливы равенства

S ly.2 (f)=H1 (f)Sn (/) [1 - у 1 2 (/)],

(5.86)

S„.2(/)=Sn(/)[l-yM/)l.

(5-87)

Syy* (/) = I (/) 12SU(/) [ 1- v b (/)].

(5.88)

При этом функция частной когерентности, определяемая форму­ лой (5.84), принимает вид

Yii/.2 (/) = 1 при всех /.

(5.89)

Из сказанного должно быть ясно, что функция частной когерент­ ности позволяет обнаружить линейною зависимость между про­ цессами Л*,(/) и Л//(/) даже в том случае, если функция обычной когерентности между процессами x,if) и у(і) такой зависимости не показывает. Более того, этот результат остается справедливым независимо от наличия или отсутствия корреляции между двумя входными процессами.

Функции обычной когерентности, связывающие процессы хх (t) и y(t) и процессы х2(і) и у{і), имеют вид

•Ьц \f)Syy (!)

(5.90)

 

\ s 2y ( f )\2

(5.91)

■$22 (i)Sy (!)

 

^Эти две функции когерентности мсгут быть равны единице или

быть меньше единицы в зависимости от рассматриваемой ситуа­ ции. Функции обычной когерентности могут оказаться как оши­ бочно высокими, так и ошибочно низкими. Рассмотрим конкрет­ ные примеры.

184

Глава 5

В частном случае, когда входные процессы некоррелированы, функции 5 12(/) = у12(/) = 0. При этом функции обычной когерент­ ности принимают вид

*.2

m ___________ \ H l

( F) \

zS n ( f ) _______

У і у

Ч ) -

I Ні

I *Sll

іп +

I н2(Л 12S22 (/)

..а

/ « _____________ 1Я2 (Л I »sw tf)_________

V2J, (.л—

! Ні (п I *Su (/) +

I я 2 сл 12522 (Л

Таким образом, в случае некореллированных входов

У%(П + У2у(/) = !•

(5.92)

(5.93)

(5.94)

а*(і>

Р и с . 5.9. Пример ошибочно высокой когерентности.

Это равенство показывает, что если одна из функций обычной ко­ герентности не равна нулю, то другая обязательно будет меньше единицы. Заметим, что при вычислении функции обычной коге­ рентности у2y(f) между процессами хД/) и у(І) влияние процесса'/ х2Д) проявляется в создании шума в выходном процессе (Д/).4- Этим объясняется ошибочно низкое значение функции обычной когерентности, маскирующее существование строго линейной за­ висимости между процессами хД^) и y(t). Однако, как следует из равенства (5.89), при вычислении функции частной когерент­ ности, связывающей процессы xx(t) и у((), влияние процесса х2Д) исключается как из хД/), так и из y(t), что позволяет получить истинное значение когерентности, равное единице.

Пример ошибочно высокой когерентности показан на рис. 5.9. Пусть значение функции когерентности между двумя процессами xi(t) и у{() оказалось близко к единице. Казалось бы, это дает основание считать, что эти функции можно рассматривать как вход и выход некоторой линейной системы. Но предположим, что существует и третий процесс х2(/), когерентный с хД/) и также дающий вклад в у(і) после прохождения через линейную систему-Д

В этом случае высокая когерентность

между хДО и y(t) может

всего лишь отражать

тот факт, что

когерентность между хДЛ

и х2(Л также высока,

а процесс хг{і) связан через некоторую ли­

нейную систему с y{t). В действительности же процессы Хі(Л

Процессы на входе и выходе физических шстем

185

ки y(t) могут быть не связаны никакой физической системой. Если в этом случае построить функцию частной когерентности между хг(і) и 2/(0 , то она, вероятно, окажется очень малой, близкой к нулю. Однако функции обычной и частной когерентности между

х2( 0 и у(і) будут близки к единице.

Пример 5.7. Ошибочно высокая когерентность. Рассмотрим линейную систему с тремя входными процессами при наличии инструментального шума на выходе; тракты системы представляют собой низкочастотные ÄC-фильтры (рис. 5.10). Заметим, что

Р и с . 5.10.

Линейная система с тремя

входным!

процессами щи і

 

личин инструментального

шума на

выходе.

И м о . z,(2), z,(0

н 2 ,(0 — четыре реализации независимых гауссовских процессов с иуде-

 

выми средними значениями и единичными дисперсиями.

входной процесс х2( 0 коррелирован с двумя другими входами Xi(t) и x3(t). Какова когерентность процесса y(t) с процессом х2(і), прошедшим через соответствующий тракт системы?

Для ответа на этот вопрос были вычислены функции обычной и частной когерентности между входным процессом x2(t) и вы­ ходным процессом y(t) для диапазона частот от 0 Гц примерно до 1000 Гц с разрешающей способностью 32 Гц и при длине реали­ зации, равной 1 с. Результаты вычислений представлены на рис. 5.11. Как видно из рис. 5.11, а, функция обычной когерент­ ности процессов 2/(0 и хг( 0 весьма велика во всем диапазоне ча­ стот и принимает значения почти от единицы на низких частотах до 0,97 на более высоких частотах. Однако это в значительной мере объясняется тем, что процессы хг(і) и х;(і) когерентны с про­ цессом х2( 0 и также дают непосредственный вклад в выходной процесс 2/(0 после преобразования их в соответствующих трак­

186

Глава 5

тах системы. Если вычислить функцию частной когерентности процессов хг(/) и у{1), то, как видно из рис. 5.11, б, ее значения оказываются несколько меньшими (до 0,75).

Р и с . 5.11. Функции когерентности между входным процессом x2{f) и выходным процессом y(t) для системы, изображенной на рис. 5.10.

u — функция обычной когерентности; б — функция частной когерентности.

5.4.3. Функции множественной когерентности

Процесс у({) на выходе линейной системы со многими входами и при наличии аддитивного шума (рис. 5.6) описывается равенст­ вом

ч

Z (0 —у (0 — hi (x)xi —т)бт,

(5.95^

' - ‘ о

Здесь г{()"— остаточный шум, определяющий отклонение модели от линейной и неизвестные внешние помехи (если таковые име-

Процессы на входе и выходе физических систем

187

^ются). • Величина z{t) сходна с остаточной случайной

величиной

Дy(t) в формуле (5.61) в том смысле, что она указывает, насколь­ ко точно может быть предсказан выходной процесс y(t) только по линейным преобразованиям входов х;(/), і — 1 , 2 , ..., q, в при­ сутствии шума. Заметим, что процессы y(t) в соотношениях (5.95) и (5.33) не идентичны.

Автокорреляционная функция остаточного стационарного слу­

чайного процесса

(z(/)|

 

определяется в виде

 

 

R zz (Т)=М [Z (t)z (/ + т)]= Л>9у

Г h t ( ß ) R ffl ( T - ß ) d ß -

 

 

 

q ■° °

i-1 £

 

 

 

 

 

 

 

 

 

—2

 

f hi(a)Riy(%+a)da +

 

 

 

 

i- l o

 

 

 

 

+ 2q

q

f

00fh

( a ) h j f ë ) R , j ( x + а —

ß ) d ß d a .

(5.96)

 

*-i / - 1V

 

 

 

Но,

как показано в формуле (5.46), весовые функции fy(ß)

долж­

ны

удовлетворять

условию

 

 

 

 

Я « Л * )= І]Т м Р )Я н (‘'-Р Д О -

'

(5-97)

 

 

 

 

/- 1

о

 

 

і^Это условие вытекает, как и при выводе уравнения (5.59), из требования такого выбора функций A; (ß), которые минимизи­ руют средний квадрат отклонения реального выходного процесса у{1) от его линейной аппроксимации только процессами xt(t).

Оно непосредственно следует также и из уравнения (5.95) в пред- | положении, что остаточный шум z{() некоррелирован с процессами хг(/), а средние квадраты ошибок не учитываются. И наоборот,

из формулы (5.97) следует, что процессы z(t) и xt(i) некоррелированы. Таким образом, подстановка выражения (5.97) в равен­ ство (5.96) дает

(5.98)

Выполняя преобразование Фурье обеих частей последнего ра­ венства, получаем выражение для спектральной плотности оста- ■Я5чных случайных величин

5й (Л =5 га( / ) - 2 Hi (f)Syi ( f ) = s y y (f) U - y U m ,

(5.99)

t=l

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