
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf168 |
Глава 5 |
ранее. Как следует из формул (3.83), функция когерентности пр4 всех / удовлетворяет условию
О < ?!„(/) < 1 - |
(5.24) |
Для того чтобы исключить влияние дельта-функций в начале координат, следует, прежде чем пользоваться двумя последними соотношениями, центрировать подвергаемые анализу данные. Заметим, что функция когерентности аналогична квадрату нор мированной взаимной корреляционной функции р^„(т), определен ной формулами (3.52) или (3.54).
Для линейной системы с постоянными параметрами справед ливы равенства (5.8) и (5.11), подстановка которых в_(5.23) дае?
Ѵ5у(/) |
\Н (О ГС?(Л |
(5.25) |
с , U ) \ h ( f ) \ 2Gx w — u |
Следовательно, в идеальном случае линейной системы с постоян ными параметрами и одним входом при полном отсутствии помех на входе и выходе функция когерентности будет равна единице.
|
|
|
u(t) Система |
№ |
y(t) |
|
|
т - |
► 0 ------ ►*© |
|
|
Mt) |
|
Р и с . |
б.ЗЛЛинейная |
система с одним |
входом |
при наличии инструмен- |
||
|
|
|
тального |
шума. |
|
|
Если |
процессы |
x(t) |
и 2/(0 полностью |
независимы, то функция |
когерентности будет равна нулю. Если функция когерентности оказывается больше нуля и меньше единицы, то возможны три
случая: |
|
содержат посторонний |
шум; |
|
а) |
результаты измерения |
|||
б) |
система, |
связывающая |
процессы x(t) и у(і), |
нелинейна; |
в) |
выход 2/(0 |
определяется |
не только входом х(0 , |
но и други |
ми входами.
При анализе линейных систем функция когерентности yly(f) позволяет выделить ту часть среднего значения квадрата выход ного процесса 2/(0 , которая на частоте/определяется входным про цессом x{t). И обратно, величина [1 — y%(f)l служит мерой той части среднего значения квадрата процесса 2/(0 , которая на ча
стоте/ не зависит от процесса x(t). с- |
А. |
Пример 5.5. Влияние постороннего шума. Предположим, |
что |
измеряемые процессы на входе х{1) и выходе 2/(0 включают в себя истинные сигналы u(t) и v(t) и некоррелированные шумовые со ставляющие п(і) и m(t) соответственно, как показано на рис. 5 .3 .
Процессы на входе и выходе физических систем |
169 |
Измеряемые на входе и выходе процессы определяются выраже ниями
x(t)=u(t) + n (О, y(t)=v (t) + m(t).
Измеренные значения спектральной плотности при идеальных ^условиях связаны соотношениями
Gx (n = Gu(n + Gn(f),
Gu(f) = G0(f) + Gm(n,
GxU {f)=GU0(/).
'JCftf
Р и с . 5.4. Самолет, летящий в турбулентном воздушном потоке.
В х о д н о й п р о ц е с с дг(0 — с к о р о с т ь в е р т и к а л ь н ы х п о р ы в о в в е т р а ; в ы х о д н о й п р о ц е с с
y(J) |
— в е р т и к а л ь н о е у с к о р е н и е с а м о л е т а . |
|
В данной задаче искомая функция когерентности
|
/£\_ |
I Gw (Л 12 |
|
(5.26) |
|
УиЛ[) |
ои (f)Gv (I) |
|
|
|
|
|
||
Измеренное же значение функции когерентности составит |
|
|||
У%(}) = |
\Gxy (f) I; |
\GUV{D\' |
|
|
G x U ) G y Щ ~ [üu (/) + Gn (/)]IGC(/) + Gm (/)1 |
|
|||
|
Y?,o ff) |
< Ѵ Ш . |
(5.27) |
|
1 + (NJGJ + (N2/G2) + (W1/G1)(/V2/G2) |
||||
где |
Ni = Gn(f)> |
G±= Gu(/), |
|
|
|
|
|
||
|
N ^ G m{f), |
Gt ~G0{f). |
|
|
Зависимость (5.27) позволяет оценить вероятные значения функ ции когерентности, когда в процессах на входе и выходе присутст вует некоррелированный посторонний шум. Заметим, что, по скольку величины Ni и іѴ2 положительны, измеренное значение функции когерентности будет теоретически меньше искомого. Отметим также', что вклад других некоррелированных входных дЛоцессов, помимо х({), проявится в выходном процессе как влия ние некоррелированного постороннего шума.
Пример 5.6. Измерение функции когерентности. Рассмотрим самолет, летящий в турбулентном воздушном потоке (рис. 5.4). Пусть входным процессом x(t) будет скорость вертикальной со
170 |
Глава 5 |
ставляющей порывов ветра, измеряемая датчиком, установлен ным в носовой части самолета; процесс на выходе у(1) есть верти кальное ускорение самолета в единицах g, измеряемое акселеро метром, который находится в центре тяжести самолета. Пример спектральной плотности и функции когерентности для реальных данных такого типа представлен на рис. 5.5. Функции спектраль ной плотности вычислялись для диапазона частот от 0,1 до 4,0 Гі^ при разрешающей способности 0,05 Гц и длине реализации
10 мин.
Р и с . 5.5а. Функция когерентности между скоростью порывов ветра результирующим ускорением самолета.
Как видно из рис. 5.5а, входной процесс (скорость вертикаль ных порывов ветра) и выходной процесс (вертикальное ускорение самолета) обнаруживают довольно значительную когерентность, достигающую 0,8—0,9 в диапазоне частот примерно от 0,3 до 2,0 Гц. На более высоких и более низких частотах функция коге рентности близка к нулю. На низких частотах вертикальное у корение самолета определяется в основном работой пилотирующего устройства, а не влиянием атмосферной турбулентности. П оэт^у уменьшение функции когерентности на этих частотах отражает влияние на выходной процесс других входных процессов, помимо х((). Как показано на рис. 5.56, на более высоких частотах спект'' ральная плотность процесса на выходе y{t) резко уменьшается в
Процессы на входе и выходе физических систем |
171 |
^зультате того, что частотная характеристика самолета близка
кхарактеристике низкочастотного фильтра, а спектральная
плотность процесса на входе уменьшается с ростом частоты.
Р и с . |
5.5б. Спектральные плотности скорости^порывов ветра и резуль- |
- |
тирующего ускорения самолета. |
tдругой стороны, уровень шума в датчиках и регистрирующих
. тройствах обычно не убывает с ростом частоты. Поэтому убы вание функции когерентности на сравнительно высоких частотах объясняется, по-видимому, влиянием постороннего аппаратур ного шума. На этом пример 5.6 заканчивается.
172 |
Глава 5 |
При использовании функции когерентности в задаче оценив вания частотной характеристики линейной системы можно рас сматривать эту функцию как отношение квадратов двух различ ных мер амплитудной частотной характеристики системы. Со гласно равенству (5.8), одна из этих мер есть
|Я ( / ) |? = 0 л 7 г |
(5,28) |
Другая мера, согласно'формуле (5.11)," есть
\Н (I)\a2 = |
(5.29) |
Отношение этих функций определяет функцию когерентности
I Схи {{) 12 |
Уху ( / ) • |
(5.30) |
|
\Н U)\?~GxU)GyU) |
|||
|
|
На практике измеряемые значения функции (5.30) находятся между нулем и единицей. Оценки амплитудной частотной ха рактеристики (5.28), получаемые по спектральным плотностям входного и выходного процессов, будут смещены всегда, кроме случая yly(j) — 1. Однако оценка амплитудной частотной характе ристики (5.29), получаемая по спектру входа и взаимному спект ру входа и выхода, будет смещена только при наличии посторон него шума на входе; при наличии постороннего шума на выходе оценка смещена не будет (при очень длинных реализациях). В частности, формула (5.29) позволяет получить несмещенную оценку амплитудных частотных характеристик для системы с(£ многими некоррелированными входами. Дальнейшее обсуждение этих вопросов проводится в гл. 6 , где показано, что точность несмещенной оценки частотной характеристики возрастает по мере приближения функции когерентности к единице.
5.3. Линейные системы со многими входами
Рассмотрим теперь линейную систему с постоянными^параметрами, на вход которой поступает несколько стационарных случайных процессов1). Пусть на вход системы поступает q про цессов и наблюдается один выходной процесс. По соображениям удобства здесь будут использованы двусторонние спектральные плотности S(f), а не односторонние измеряемые функции G(f). Все полученные в этом разделе результаты будут справедливы также и при замене 5 на соответствующие G. д )*
*) Процессы должны быть не только стационарными, но и стационарно связанными, т. е. их смешанные моменты также должны удовлетворять усло виям стационарности.— Прим, перев.
Процессы на входе и выходе физических систем |
173 |
£.5.3.1. Автокорреляционные и спектральные соотношения
Рассмотрим q линейных систем с постоянными параметрами, на вход которых поступают q входных процессов с гладкими свойствами і = 1 , 2 , 3, q, а на выходе измеряется один
I
Р и с . 5.6. Линейная система со многими входами.
процесс |
y(t) (рис. 5.6). Этот выход |
можно |
рассматривать как |
||
^ сумму q |
выходных процессов |
г/ДО. |
£ = 1, 2 , 3, ..., |
q: |
|
|
= |
|
|
(5.31) |
|
|
|
t~i |
|
|
|
где yt(t) |
определяется как та |
часть |
выхода, |
которая |
создается |
£-м входом при равенстве нулю других входов.
Пусть ht(х) — весовая функция системы, на вход которой по ступает процесс Хі(і). Тогда, как следует из формулы (5.1),
(5.32)
,и процесс на выходе системы есть
hi (т)хг (t—т)с£т. |
(5.33) |
174 |
Глава 5 |
Предположим далее, что x^t) — выборочные функции раз-і личных стационарных случайных процессов со средними зна чениями рг. Здесь независимо от t
|і£= М [*,(/)!, |
(5.34) |
где каждая функция xt(t) определена на соответствующем выбо рочном пространстве. Имея в виду линейность оператора мате матического ожидания, найдем математическое ожидание функ ции y{t):
Ру=М[у(()] = М |
ч |
г |
2 |
(t)jc, (/ — x)rfx |
.‘•=і о
(5.35)
Равенство (5.35) справедливо при любых рг. Если каждое р, = О, то Pj, = 0. При практических расчетах следует прежде всего цен трировать анализируемые реализации. Это значительно упро щает окончательные формулы и их интерпретацию.
Можно также |
найти автокорреляционную функцию |
Ry(т). |
|||
В предположениистационарности процесса получим |
|
||||
Ry(т) = М \у (t)y (t + т)] = М [ 2 Уі (0 S у( (t + т)] |
= |
|
|||
|
|
VU |
|
|
|
= М £ |
£ |
— ЪШ ^ |
+ |
fi)d7) |
|
1 |
І |
о |
|
|
|
2 SJ |
hM)mxi{t - s w +t- -n)№ ч = |
|
|||
|
i |
|
|
|
|
|
|
hi (l)h j (Tj)R j j ( l — q + |
x)dU -q. |
(5.36) |
|
В соотношении (5.36) функция Ru(т) определена как |
|
||||
|
Я ^ ) = Ѵ / ( Т)==М*ДС*(* М < + 'С))- |
(5-37> |
Соотношение (5.36) выполняется и при коррелированных входах. Если предположить, что все входы взаимно некоррелированно и средние значения равны нулю, то
(т) при i= j,
(5.38)
при І Ф /.
Процессы на входе и выходе физических систем |
175 |
этом случае формула (5.36) сводится к виду
<7 |
00 |
|
Ry(*)= è |
f J a*(щ т лг—п+ѵтч |
(5.39) |
|
1"n |
|
Найдем теперь соотношения для спектральных плотностей. Двусторонняя спектральная плотность S^/) стационарного слу чайного процесса определяется как преобразование Фурье авто корреляционной функции Ry(т). Поэтому преобразование Фурье обеих частей равенства (5.36) позволяет получить спектральную плотность выходного процесса у ( і ) в виде
|
|
|
00 |
|
|
|
S ,( /) = |
§Ry(r)<r-WdT= |
|
|
|||
|
|
со |
|
|
|
|
SS I I h ( Щ (rijR ij ( l — 7J+x)dU x\)dx . |
(5.40) |
|||||
/ |
/ |
о |
|
|
|
|
Для упрощения формулы |
(5.40) введем в нее |
множитель |
||||
е- 2л/7 (£—Т)) е2л/7 (ё—Т]) = 1 . |
Тогда |
|
|
|
||
|
|
ІЯ h t (%)e2n‘^ h j (~q)е~2я^ |
x |
|
||
i |
J |
—со |
0 |
|
|
|
X R tj (5— 3j + |
г)e~ :ni! (S - л + VdM-qdx. |
(5.41) |
||||
Произведем теперь замену |
переменных t — t, — .7 + x |
и d t = |
||||
— d x . Выделяя отдельныемножители, |
получим следующий общий |
|||||
результат для коррелированных входов: |
|
|
||||
|
|
СО |
СО |
|
|
|
Sy ( / ) = 2 S fК { 1 ) < м т I h j W e - ^ d - q X |
|
|||||
1 |
і |
о |
о |
|
|
|
Rij(i)e-2n'!tdt = 2 2 |
mWHjinSijif). |
(5.42) |
||||
|
|
|
1=1 / = 1 |
|
|
|
В формуле (5.42) функции S;j(/) представляют собой взаимные спектральные плотности входов х,(0 и Xj(t), т. е. Su(f) = Sx[x.(f).
Вместо Sy(f) и Si}(f) можно ввести односторонние спектральные плотности Gy(f) и Guff), поскольку общий множитель, равный 2,
176 |
Глава 5 |
сокращается. В случае |
взаимно некоррелированных входов/ |
при которых автокорреляционная функция определяется соот ношением (5.38), формула (5.42) принимает вид
(5.43)
г=і
5.3.2.Взаимные корреляционные и взаимные спектральные соотношения
Соответствующие рис. 5.6 соотношения для взаимных спект
ров получаются |
путем вычисления взаимной спектральной плот |
ности Sxiy(f) = |
Siy(f) выхода y(t) с одним из входов, скажем |
с Xi(t). Общий |
результат для коррелированных входов имеет |
вид |
|
|
(5-44) |
|
/=і |
Вывод этого соотношения аналогичен выводу предыдущего соот
ношения |
для Sy(f). Для каждого |
і = 1, 2, ..., q |
|
СО |
|
siy(f)= |
§ R iy(r)erWhch = |
|
|
e~mihj{l) |
у |
|
е-2я/Цх-& Я0.(т — |
- ^ № = S # ;.(/)S;,.(/). |
(5.45) |
/=і |
|
Здесь использовано выражение для взаимной корреляционной функции Riy(т), имеющее вид
Riy(f)=M [xt (І)у (t +■T)]=M X;i ( 0 2 |
f |
+ |
i |
0 |
|
00
= 2 f hi ® ^ ®xj V +
Процессы на входе а выходе физических систем |
177 |
f При взаимно некоррелированных входах и нулевых средних
[соотношение (5.38)] формула (5.46) сводится к виду
СО |
|
а д т) = | М Е ) Я , ( т - Е № |
(5-47) |
о |
|
Следовательно, соотношение для взаимного спектра |
входного |
и выходного процессов есть |
|
S iy([) = Hi (f)Si (f). . |
(5.48) |
Формула (5.48) представляет собой полезный и интересный результат, поскольку из него следует, что частотьая характери стика тракта системы, соответствующего входу xt(t), может быть измерена при помощи взаимных спектров независимо от влияния других входов при условии, что входы взаимно некоррелированы. Этот результат решает проблему измерения отдельных частотных характеристик при некоррелированных входах. Статистическая точность измерения определяется функцией когерентности. Этот вопрос рассматривается в гл.6 .
5.3.3.Система с двумя входами
,Рассмотрим линейную систему с двумя входами (рис. 5.7)
ихг(1), вообще говоря коррелированными. Пусть по данным
оспектральных плотностях требуется определить частотные ха рактеристики системы. Формула (5.44) в этом случае имеет вид
S1y(f) = H1(f)S11(f) + H2([)S12(f),
Sa,(/)= tfi(/)S al(/) + tfa(/)SM(/).
Решения для # г(/) и # 2(/) |
при у1 2 (f) Ф 1 |
имеют вид |
|||
Sly(f) |
s^tnSzyffv |
||||
•S22 (f)Sly ( / ) J |
|||||
|
|
||||
H l ^ = |
Su \f) [ 1 - ѵ Ы / |
Л |
|
||
Szy (f) |
S 21 (f)Siy |
(f ) |
|||
1 |
. J |
( f ) |
|||
НЛП-- |
|
( f ) $ 2 |
|||
S22 (/Я1 — V12 (/)] |
|
||||
* |
|
||||
|
|
|
|
||
где |
|
|
|
|
|
??«(/) |
I Sia (fl 12 |
|
|
||
■Su {f)S22 (f)‘ |
|
(5.49)
(5.50)
(5.51)
12—2244