
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf178 |
Глава 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