
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf188 |
|
Глава 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 + М [(М [Ф]—Ф)*Ь