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

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

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

188

 

Глава 5

 

 

где функция yl.xif)

определена в виде

 

 

 

ч

 

 

 

 

I

 

(!)

 

Ц .х ( П =

1-\

■= 1

(5.100)

 

Syy (/)

 

Syy(f)

 

Величина y\.x(f) называется функцией множественной когерент­ ности выходного процесса у{і) и всех входов xt(t), £ = 1, 2 , q. Эта функция удовлетворяет неравенству

0 < УІ.х (/) < 1 при всех /.

(5.101)

Неравенство (5.101) может быть получено непосредственно из формулы (5.100) в предположении, что Szz(f) > 0 и SZZ(J) ^ Syy(f), что вполне справедливо для рассматриваемой модели. Для дока­ зательства неравенства SZZ(J) ^ Syy(f) достаточно показать, что процессы z(t) и Хі(1) некоррелированы при всех £. В этом случае из формулы (5.95) следует, что функция Syy(f) представляет со­ бой сумму функции SZZ(J) с некоторой (всегда неотрицательной) суммой спектральных плотностей, а все функции взаимной спек­ тральной плотности равны нулю. По определению взаимная кор­ реляционная функция^процессов z{t) и xt(t) есть

/?гі (т) = М [г (ф :г(£ + т)],

где z(t) определяется уравнением (5.95):

ч' °°

(0 = у (/) — Г hj (ß)Xj (t— ß)dß.

/-i.oj;

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

R2i (т)=М [у(()х, (t + т)] _ 2 J hj (ß)M [xt (t + x)xj (t— ß)]dß =

/ “ 1 0

= [Rui W - 2 ] hj (ß)/?„ ( ^ T - ß ) r f ß . ' l

/-»o

Согласно уравнению (5.97),

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

189

Следовательно,

(5.102)

£*»(*)=о,

что и требовалось доказать.

Функция множественной когерентности, определенная фор­ мулой (5.100), включает в себя как частные случаи функции ІЬбычной и частной когерентности, определенные ранее. Так, функцию частной когерентности можно рассматривать как функ­ цию множественной когерентности для остаточных (условных)

случайных

величин. Рассмотрим формулу (5.100) для

случая

4 = 1 :

 

 

%

? ;.,(/)= " '( ' ’У 1 .

(б.іоз)

Преобразование Фурье уравнения (5.97) в данном случае имеет вид

5іЛ/) = ^ (/)5 и (Л -

(5.104)

Исключая из равенства (5.104) функцию Яг(/), получаем

(5.105)

а это выражение совпадает с выражением для функции обычной когерентности. Таким образом, функция множественной когерент­ ности содержит функцию обычной когерентности как частный случай.

^ В общем случае частотные характеристики Яг(/), входящие в фі^шулу (5.100), могут быть определены из решения системы линейных уравнений

я

 

Siy (/) = 2 Ні

(/>»

і = 1 > 2, • • • • ? .

(5.106)

получаемых при помощи преобразования Фурье

обеих частей

равенства (5.97).

 

 

 

Рассмотрим теперь идеальный случай, когда остаточный шум

z(t) равен

нулю. Согласно формуле (5.42), функция Suy(f), со­

держащаяся

в уравнении

(5.100),

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

виде

яя

Syy ( / ) = 2

2 н ‘

Wsu (/).

(5.107)

Подстановка в формулу

(5.100) выражения, комплексно сопря­

же ого выражению (5.106),

приводит к результату

 

Q Я

 

 

2

Ц я і(/)я 7 № <•(/)

 

УІ.АП

 

Syy (f)

1.

(5.108)

 

 

 

 

190

Глава 5

поскольку Syy(f) = SyU(/), а знаменатель дроби в (5.108) ес£$ именно Sty if). Таким образом, при отсутствии помех и строго!! линейности системы со многими входными процессами и одним выходом функция множественной когерентности равна единице, Очевидно, что этот результат можно получить и более непосредст венным путем, просто полагая в формуле (5.100) Szz(f) = 0.

Более детальный анализ функций множественной когерент ности, а также функций частной когерентности и других типов этой функции содержится в работах [20] и [26]. Вопросы практи ческой интерпретации функций множественной когерентности рас­ смотрены в гл. 6 этой книги.

&

Р и с . 5.12. Функция множественной когерентности между тремя вход­ ными процессами и выходным процессом для системы, изображенной на рис. 5.10.

Пример 5.8. Измерение функции множественной когерент­ ности. Рассмотрим систему с тремя входными процессами и од­ ним выходным процессом, обсуждавшуюся ранее в примере 5.7 и показанную на рис. 5.10. Функция множественной когерент­ ности выходного процесса y[t) и всех трех входов хг(/), і = 1, 2 , 3, изображена на рис. 5.12. Как видно из рисунка, числовые значе­ ния функции множественной когерентности весьма велики в<| всем диапазоне частот (свыше 0,95), но с ростом частоты проис­ ходит некоторое их уменьшение. Единственной причиной того, что функция множественной когерентности не равна единице во всем диапазоне частот, является присутствие инструментальнаго шума в процессе на выходе системы. Поскольку спектральная плотность инструментальных помех на всех частотах постоянна^ а процессы хг( 0 подвергаются низкочастотной фильтрации, относи* тельный вклад инструментального шума в суммарный сигнал на'

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

191

Я^іходе системы возрастает с частотой. Именно на это обстоятель­ ство и указывает уменьшение числовых значений функции мно­ жественной когерентности с увеличением частоты.

5.4.4.Матричная форма результатов

}Приведенные выше формулы могут быть выражены более сжа­ то в матричной форме, что позволяет, кроме того, более просто получить некоторые новые результаты. Определим вначале <7-мер-

ный вектор входных процессов

^ X (0 = [*1 (fl, Х2 (fl...........Х„(і)]. (5.109)

Определим также g-мерный вектор частотных

характеристик

Н (/) = [Я Х(/), Я а ( / ) , . . . , Я ,(/)].

(5.110)

Далее определим g-мерный вектор взаимных спектральных плот­ ностей выходного процесса у(!) и входов лД0 :

(/) - [Sly (/), S 2, ( / ) , . . . , s w (/)!,

(5.111)

где

 

SiB(f)=Sxiy(f), i = 1, 2Д . . , q.

(5.112)

Определим, наконец, матрицу размерности q X q взаимных спек­ тральных плотностей всех входных процессов xt(t):

'Snif)

S12(f). .

. S lq([f

 

S« (/) = Sn (/)

S„(f). .

. S2, (/)

(5.113)

Л і ( /)

4 (/) ■ ■

 

 

где

 

 

 

Su ( n = S xlxj(f),

i , j = \ , 2 , . . . , q .

(5.114)

Соотношение (5.42) может быть записано в матричной форме:

Syy(f) = W{f)Sxx( f W ’ {f),

(5.115)

^де Н*'(/) — комплексно-сопряженный транспонированный век­ тор. В развернутом виде равенство (5.115) записывается так:

М /) = [Яі(/), я , (fl...........Hq(f)\

X

 

 

'Sn(f)

Sxa(fl. .

. S lq(f)~

■яг (/)■

 

S2i (/)

5ла (/) . .

. S2q(/)

ЯИ/)

(5.116)

ß qi (/)

s q2(f). .

. s qq([)_

_ял/)_

 

192 Глава 5

Вектор-столбец в правой части есть вектор, комплексно-сопря женный с транспонированным вектором Н (/):

Щ(П ~

 

Н*'(/)=

н т

(5.117)

 

 

_ Щ(П

 

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

равенстве (5.115) функция Syy(f) — скалярная

характеристика,

а другие

функции — векторы.

Равенство

(5.115) служит также правильным представлением равенству (5.43) для некоррелированных входов; в этом случае матрицй- SXX(J) диагональна. Следует, однако, подчеркнуть, что этот ре­ зультат справедлив только для идеального случая отсутствия постороннего шума.

Система (5.44) может быть записана в матричной форме для

і= 1» 2 , ..., q:

s;,(/)= s« (/)H '(/).

или в развернутом виде

S l y ( f )

^11 (/)

^12 (/) • • • Ѵ /Г

~ H i

( f r

$ 2 У (/)

 

= ^21 (/)

S 2o (/) . '

■ S 2 q ( f )

н

т

 

 

 

 

 

 

 

S gy ( f )

_

Л і(/)

$,»(/).

 

н

т

-

(5.118)

(5.1ід^

где векторы-столбцы в правой части суть транспонированные век­ торы-строки Sxy(f) и H(f). Уравнение в матричной форме (5.118) может быть решено относительно транспонированного векторастроки Н'(/) при измеренных или известных Sxy(f) и Sxx(f). Это выражение представляет собой, очевидно, систему q совместных линейных уравнений, решение которой имеет вид

H'(/)=s;i-(/)SV(/),

(5.120)

где SХЩ~ ) —матрица, обратимая матрице Sxx(f) . Уравнение

(5.120)

определяет каждую функцию #Д /) как функцию взаимных спект­ ров Siy(f) входных и выходного процессов и взаимных спектррв Sjj(f) входных процессов. Уравнение справедливо независимо от того, коррелированы или некоррелированы входные процессы. і

Для получения других соотношений, касающихся функций? множественной когерентности, определим пополненную спект-'

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

193

фальную матрицу размерности (<7 + 1) X {q + 1) для входных 'процессов xt(t) и выходного процесса y(t) в виде

'Syy(f)

Syl(f). .

. S yq(f)~

 

Suj (/)

s u (/). .

,S lq(f)

(5.121)

SУ х х (/)= s iy(f)

S21(/). .

. 52<7 (/)

Л у(/)

(/) ■ •

(/) _

 

Покажем теперь, что в идеальном случае отсутствия шума де­ терминант этой матрицы равен нулю для всех /.

^ Согласно равенству (5.44), элементы Siy(f) первого столбца матрицы Syxx, начиная со второй строки и ниже, могут быть пред­ ставлены в виде

s iy(f)= S Я, (/)S„(/),

1. 2 , . . . , q.

(5 .122)

/=1

 

 

Таким образом, каждый член S ^ ) есть линейная комбинация из элементов данной строки. Согласно формуле (5.42), которая спра­ ведлива только для идеального случая отсутствия шума, остав­ шийся член Syy(f), принадлежащий первой строке и первому столб­ цу, можно записать в виде

^ (/) = ! ]

2

я ; (f)Hj (n s u (/) = 2 я ; (/)stlf(/) =

(/)Syi (/).

/ = і

/ = 1

і = і

/ = і

 

 

 

(5.123)

Последние два равенства получаются в результате подстановки выражения (5.44) в (5.42) с учетом равенства Syiß ) = S*yy{f). Таким образом, функция Syy(J) есть линейная комбинация, состав­ ленная из элементов Syi(f) первой строки матрицы' Syxx(f). Итак, в случае отсутствия шума матрица Syxx(f) обладает при всех / следующим свойством: элементы ее первого столбца представляют собой линейные комбинации из соответствующих элементов дру­ ги х столбцов. Согласно хорошо известной теореме (см., напри­ мер, [32]), детерминант этой матрицы равен нулю.

Рассмотрим теперь более общий случай присутствия остаточ­ ного шума, обсуждавшийся ранее в подразд. 5.4.3. Как следует ^формулы (5.99), элемент Syy(f) выражается в виде

* Syy(n = i}'ßi(f)Syi(n + S22(f), (5.124)

13*—2244

194 Глава 5

что не совпадает с равенством (5.123). Однако элементы Siy(fy£ определяемые формулой (5.106), не меняются, поскольку условие (5.106) должно оставаться справедливым независимо от наличи или отсутствия шума. Следовательно, в этом общем случае, с учетом

полученного выше результата при отсутствии шума,

детерминант

I Syxx(f) I

матрицы Syxx(J) определяется как

 

 

 

\Syxx(f)\=Szz(f)\Sxx(t)\,

(5.125)

где I Sxx(f) I — детерминант матрицы Sxx(f). Очевидно, что

при

Szz(/)i=

0 детерминант матрицы Syxx(f) обращается

в нуль,

как

это и должно быть при отсутствии шума.

 

 

Соотношение (5.125) приводит к очень простому выражении^ для функции множественной когерентности, определенной ране*?' формулой (5.100), через детерминанты матриц Syxx(f) и Sxx(f)

t h r r -

Для проверки этой формулы заметим, что при q = 1

' l S y y ( f ) Syx (/)'

А * ( / ) .

Sу х X А Л / )

Здесь

\Syxx\= S xx(f)Syy( f t - \ S xy(f)\* и \Sxs(f).\=Sxx(/).

<5Л26>

(5.127)

Произведя подстановку в соотношение (5.126) и сокращая общие

члены, получим выражение

\ s x u U ) \ 2

і

у і . А П =

(5.123W

$хх (f)Syy U)

 

'Ф\

совпадающее с функцией обычной когерентности.

 

Другое эквивалентное представление функции множественной

когерентности процесса y(t) и всех процессов хг( 0

дается форму­

лой

 

(5.129)

?А(/) = 1-ІАЛ /АУ(/)Г\

где Syy(f) — первый элемент главной диагонали матрицы S^.(/), обратной матрице Syxx(f). Этот результат соответствует формуле (5.126), поскольку обратная матрица S получается делением матрицы, транспонированной и присоединенной матрице Syxx(f), на ее детерминант | Syxx(J) |. Матрица, присоединенная матрице Syxxif), находится путем подстановки вместо элементов S^if) их/ алгебраических дополнений Ai}(f):

 

Аyy{f) Аіу(П А2,(/).

qy{f)

 

Ауг(П А ц (/)

А21(/).

.

ql(f)

ss .(/)

 

 

 

 

V 1

J y x x I Aj/2 (/)

А12 (/)

А22 (/).

 

АЛ/)

 

 

 

 

 

4

 

_ А Л /)

А ? (/)

А'lqЛ1/).

.

. А Л /)_

 

 

 

 

 

W'

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

195

Следовательно,

 

 

Суу /С) __

^UV (/)______I Яг* (/) I

(5.130)

ч>

|SyXJC(f)| - | S ^ ( / ) | •

 

причем, согласно равенству (5.121), Ayy(f) = |S ^(/)|. Таким об­ разом, равенства (5.126) и (5.129) дают два различных, но эквива­ лентных способа вычисления определяемой в общем случае форму­ лой (5.100) функции множественной когерентности.

Упражнения

1- Покажите, каким образом можно вычислить для изображен­ ной на рис. 5.1 модели функции Gx(f) и H(f) при известных спект­ ральных плотностях Gy(J) и GXy{f) и с учетом ограничений, налагае­ мых на соответствующие формулы.

2 . Убедитесь в справедливости результатов, приведенных в примерах 5.1 и 5.2.

3. Убедитесь в справедливости результатов, приведенных в примерах 5.3 и 5.4.

4. Рассмотрите модель, изображенную на рис. 5.3, где

G„ ( / ) = 5, Gm( / ) = 0 .

^Найдите следующие величины:

а) G0(f) и Gm(f);

б)

y'Uif) и y l y i f ) .

5.

В случае системы со многими входами, когда входные про

цессы имеют отличные от нуля средние значения и не коррелиро­ ванъ! между собой, формула (5.38) принимает вид

при і = /,

при і Ф /.

Покажите, что в этом случае формулы (5.43) и (5.48) запишутся --в виде

Sy (/) = £ I Н і (/) 12S,- (/) + 23 2 Hi № (0)№-s (/).

1 /=1

(/ +/)

Siy(f)=Hl (f)Si (f)+ J ] я у(0)р^6(/).

/-1

cf*n

13»

196

 

 

Глава 5

6 .

Рассмотрите модель с двумя входными процессами (рис. 5.7Jf

где

 

 

 

 

 

" iW -F T T T ’

H^ = i h r -

Пусть

входные процессы удовлетворяют условиям

 

 

JR1 (x) = 36(x),

Ga (/)= 12, 012(/) = 8 .

Найдите

следующие величины:

а) Yi2</);

 

б) Я„(х) и G,(/);

3

в)

 

иGiö(/);

г)

7. Используя формулу (5.84), вычислите функцию частной когерентности для системы, определенной в упражнении 6 .

8 . Выведите формулы (5.86) и (5.88).

9.Используя формулу (5.100), вычислите функцию множест­ венной когерентности для системы, определенной в упражнении 6 ,

впредположении, что спектральная плотность постороннего шума z{t) имеет вид Sz(f) — 19/(25 + /2)-

10.Убедитесь в справедливости соотношения (5.125) для слу­

чая, рассмотренного в упражнении 9.

Я

I

СТАТИСТИЧЕСКИЕ ОШИБКИ ПРИ АНАЛИЗЕ СЛУЧАЙНЫХ ПРОЦЕССОВ

^ Как отмечено в гл. 4, точно определить характеристики слу­ чайны х величин по выборочным данным невозможно. По выборке конечной длины можно найти лишь оценки интересующих иссле­ дователя параметров. В гл. 4 приведены сведения о точности определения ряда основных параметров случайного процесса по выборке объема N независимых измеренных значений этого про­ цесса. В настоящей главе рассматривается вопрос о точности оценок параметров для непрерывной реализации процесса дли­ ной Т. Предполагается, что анализируемые реализации принад­ лежат стационарному (эргодическому) случайному процессу с не­ прерывным временем. Рассмотрены оценки как основных пара­ метров процесса, так и частотной характеристики и функций коге­ рентности, полученных на основе оценок спектральной плотности для систем с одним и многими процессами на входе. Здесь анали­ зируются только те ошибки, которые обусловлены статистиче­

ск о й изменчивостью исследуемых процессов. Сведения о других ошибках, связанных с получением исходных данных и их обра­ боткой, приведены в гл. 7.

6.1. Понятие о статистических ошибках

Точность оценки некоторого параметра случайного процесса, полученной на основании выборки, характеризуется, как указано в разд. 4.1, средним значением квадрата ошибки

(с. з. к. о.):

С. з. к.

о.=М [(Ф —Ф)2],

(6.1)

где CD— оценка параметра

Ф.

 

Равенство (6.1) можно представить в виде М[(Ф—Ф)2] = М[(Ф—М [Ф ]+М [Ф ]—Ф)2] =

= М [(Ф—М [Ф])2] + 2М [(Ф—М [Ф])(М [Ф]—Ф)1 + М [(М [Ф]—Ф)*Ь

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