книги из ГПНТБ / Марченко Б.Г. Метод стохастических интегральных представлений и его приложения в радиотехнике
.pdfВыражение (478) в приложениях имеет самостоятельный фи зический смысл, так как в рассматриваемом случае оно определяет мощность «постоянной» составляющей при условии, что процесс V (t) —стационарный
Сравнивая два приведенных выше способа, которыми было получено выражение (475), можно заметить следующее. Первый путь, использующий разложение (332), является универсальным, но не наиболее простым. Это особенно видно в тех случаях, когда ряды допускают замену эквивалентными выражениями, содержащими конечное число элементарных или специальных функций. Следова
тельно, в тех |
случаях, когда искомая корреляционная функция |
выражается в |
конечном виде через элементарные или специаль |
ные функции, |
удобнее второй путь анализа, использующий переход |
к новым компонентам случайного вектора. Однако в тех случаях, когда такие представления получить не удается, следует пользо ваться разложениями в ряд (314) или (332) и применять для постро ения графиков и таблиц корреляционных функций ЭВМ.
Заметим, что можно применить и третий путь получения кор
реляционной |
функции, |
воспользовавшись результатами § 6 гл. II |
|||
в одномерном |
случае, |
т. е. совершить переход |
к независимым |
||
компонентам |
при вычислении корреляционной |
функции. Этот |
|||
путь анализа |
в настоящем |
параграфе |
не рассматривается, так как |
||
о нем подробно будет идти |
речь в § 5 |
настоящей главы, из резуль |
|||
татов которого, как частный случай, можно будет получить некото рые формулы этого параграфа.
При постоянной амплитуде выражения (475) и (477) являются самыми общими для случая угловой модуляции в присутствии га уссовой помехи, так как в них учитываются стационарный и неста ционарный случаи гауссовых помех и при соответствующем опреде лении ядра ф (т, f) в V (і) и функции N (t), согласно (466) и (467), они оказываются пригодными как для случая фазовой, так и для случая частотной модуляции. Для более глубокого и наглядного раскрытия свойств этих функций рассмотрим ряд частных задач анализа сигнала вида (465), приводящих к корреляционным функ циям (475) и (477) и имеющих самостоятельное прикладное и теоре тическое значения.
1. Пусть в (465) функция N (t) == 0, а v (t) — линейный ста ционарный гауссов процесс, описанный в настоящей главе или по лученный путем фильтрации из белого шума, как это упоминалось в § 3 настоящей главы.
Тогда выражение |
для смещенной |
корреляционной |
функции |
(475) упрощается следующим образом: |
|
|
|
$ (s) = |
- у - е~т [eh(s) + |
е~Ш) cos 2<pJ, |
(479) |
а корреляционная функция процесса (465), согласно (477), в этом
140
случае определяется выражением
|
|
Я(5) = |
|
-^ -< Г Л(0> M s ) |
1 4 - [е |
h(s) — 1] cos 2ф) = |
|
|||||
= Лае‘ |
sh |
h(s) |
cos ф sh |
h |
(s) |
+ sin. „ ф ch, I —Y ~ |
(480) |
|||||
|
2 |
|
-А(0) |
|
|
2 |
|
2 |
h (s) |
|
||
|
|
|
|
|
|
|
|
|
|
|||
Величина смещения кор |
|
|
|
|
|
|||||||
реляционной |
функции |
в |
|
|
|
|
|
|||||
рассматриваемом |
|
случае |
|
|
|
|
|
|||||
равна |
|
|
|
|
|
|
|
|
|
|
|
|
АВ (s) = |
(s) — В (s) |
- |
|
|
|
|
|
|||||
= 4oe_MO)cos2? . (481)
Теперь предположим, что процесс V (t) получен в ре зультате фильтрации бело го гауссового шума при по мощи /?С-фильтра с посто янными параметрами, т. е. выражен формулой (456) с корреляционной функцией, определенной согласно (457)
сучетом, что к2 Iw (1)] =
=1. Рассмотрим для этого случая графики корреля
ционных функций, |
постро |
|
|
енные согласно полученным |
|
||
выше выражениям. |
|
|
|
На рис. 10, а приво |
|
||
дится |
семейство |
кривых |
|
B(s)/B (0), построенных с |
|
||
учетом (480) для случая #С-фильтра с корреляционной |
функцией, |
||
определенной согласно (457) при А0 = 1; ф = 0; а |
= 0,5; 1; 2; |
||
4; 7 |
и s £ [0, 2]. В этом случае (480) принимает вид |
|
|
|
|
а |
|
|
B(s) = |
2е~ 2 sh2(-^<T “|S|) , s g(—oo, oo). |
(482) |
Так как функция (482) симметрична относительно оси ординат, то на графиках рис. 10 приводятся значения В (s)/B (0) только для s >> 0. Функция (482) в нуле не дифференцируема, причем недиффе ренцируемость преобразованного случайного процесса (465) в рассматриваемом случае вызвана недифференцируемостью исход ного процесса ѵ (t).
Из семейства графиков рис. 10, а видно, что с увеличением а кор реляционная связь отдельных значений процесса (465), находящихся на фиксированном расстоянии s, монотонно ослабевает. Для рассмат риваемого случая дисперсия, согласно (482), определяется формулой
В (0) = 2е 2 sh2'(*)■ <х> 0.
141
Нетрудно видеть, что
= Г ~ sh (-J-) (ch -j- — sh-J-) > 0 при а > О,
следовательно, с изменением а от 0 до оо дисперсия монотонно воз растает от 0 до оо.
При ф = |
0 и А (s), определеннойасогласно (461), |
смещение кор |
||||||||
реляционной функции Aß (s) = |
е 2 , т. е. отлично |
от |
нуля. Это |
|||||||
|
|
объясняется тем, |
что дина |
|||||||
|
|
мическая характеристика в |
||||||||
|
|
этом |
случае |
является |
чет |
|||||
|
|
ной |
функцией. При |
ф |
== |
|||||
|
|
|
л |
смещение |
АВ (s) |
=з |
||||
|
|
= |
Т |
|
||||||
|
|
0, |
|
т. |
е. |
‘Р (s) = В |
(s). |
|||
|
|
Вэтом случае динамическая |
||||||||
|
|
характеристика |
является |
|||||||
|
|
нечетной функцией. Таким |
||||||||
|
|
образом, |
начальная |
фаза |
||||||
|
|
Ф |
влияет |
|
на |
значения |
||||
|
|
корреляционной |
функции |
|||||||
|
|
В (s), смещенной |
корреля |
|||||||
ционной функции ^3 (s) и на значения Aß (s). |
|
|
В (s)/ß (0) |
при |
||||||
На рис. |
10, б приводится |
семейство |
кривых |
|||||||
Ф = -у. В этом случае А0, s и а принимают такие же значения, как
и на рис. 10, а. |
(0), В (0)/$ (0) и Aß (0)/$ (0) |
|
На рис. |
11 приводятся графики |
|
как функций начальной фазы ф для |
двух случаев: а — 0,5 (сплош |
|
ные линии) |
и а = 1,0 (пунктирные линии), построенные согласно |
|
(479), (480) |
и (481) при Л0 = 1 и h (s), определенной согласно (457). |
|
Анализ этих графиков позволяет сделать вывод, что при изменении
пл
Фот 0 до — вследствие появления асимметрии динамической харак
теристики (471) мощность постоянной составляющей падает, а переменной — возрастает, т. е. происходит перераспределение мощ ностей.
2. Пусть в (465) функция N (t) = 0, а ѵ (і) — линейный неста ционарный гауссов процесс, полученный в результате прохождения белого шума через фильтр, параметры которого меняются со време нем, или нестационарный гауссов шум, описанный в § 3'настоящей главы.
Выражение (475) для смещенной корреляционной функции в этом случае принимает вид
А2 |
h(t>t)-j- h(t s,t s) |
|
Ло |
j { ^ ^ > + e -W+*>cos2q>}, |
|
V (S, 0 = -if-exp |
2 |
|
|
(483) |
|
|
|
142
а выражение для корреляционной функции (465) согласно (477):
B(s, = |
exp (-----L[h(t, f) + h(t + s,t + s)]J X |
|
X {ew 't+s) — 1 + [érft(M+s) — 1J cos 2ф}. |
(484) |
|
Величина смещения корреляционной функции в рассматрива емом случае определяется так:
АВ (s, f) = Ао ехр {— [h(t, t) + h(t -f s, t + s)]J cos2Ф. (485)
Следовательно, в рассматриваемом нестационарном случае ве личина смещения уже зависит от t и s, причем, как и в стационарном
случае, при фиксированных t я sc увеличением ф от 0 до ~ она убы
вает до нуля, оставаясь положительной.
Предположим, что ѵ (t) — линейный нестационарный процесс, полученный в результате фильтрации белого гауссового шума при помощи /?С-фильтра, у которого сопротивление меняется по линейному закону, т. е. процесс ѵ (t) представим в виде (460). Им пульсная переходная функция такого фильтра определяется со гласно (458), а ее вид при а — 1 приведен на рис. 8 . Корреляцион ная функция этого фильтра определяется согласно (461) с учетом, что х2 \w (1)] — 1. В частности, при ос = 1 и а = 2 она принимает вид (462) и (463), который использовался при построении приведен ных ниже графиков.
На рис. 12, а сплошными линиями приводится семейство кри вых В (s, /), построенных согласно (484) с учетом (462) при А(І = 1,
Ф =» 0, |
а |
= |
1, s £ [—2, 31 и |
t = 0,5; |
1; 2; |
3, т. е. для выражения |
|||
В (s, t) |
= |
2АІ |
ехр |
(t + s) 1(0 + П (t + s) |
sh2 |
KQiP + s) |
|
(486) |
|
2t (t + |
s) |
max (t, t + s) |
CO S2 ф . |
||||||
Из графиков видно, что с увеличением ^ дисперсия В (0, t) процес са (465) в рассматриваемом случае падает, а корреляционная связь между отдельными значениями процесса меняется более плавно по s. В (s, t) является недифференцируемой в нуле ввиду недифферен цируемости h (t, t + s), а следовательно, и преобразованный про цесс (465) также является недифференцируемым (по вероятности).
В момент і — 1 величина R (t) С рассматриваемого фильтра сов падает с параметром ^С-фильтрас импульсной переходной функцией (454), который рассматривался в предыдущем параграфе. Сравни
вая кривую В (s, 1) рис. 12, а с кривой а = 1, приведенной на
рис. 10, а, видим, что первая из них на рассматриваемом участке
[0, 1/2] изменяется немного медленнее, |
чем вторая. Зато |
дисперсия |
|
в первом случае (В (0, 1) = 0,199) в два с |
лишним раза |
выше, чем |
|
во втором (В (0) = 0,077), и только |
к |
моменту і = 2 дисперсии |
|
143
приблизительно выравниваются, а к моменту t = 3 уже
5 (0 , 3 ) < ß ( 0 ) .
Физически это можно объяснить так. Линейная система обладает памятью и средняя мощность на выходе системы с импульсной пере ходной функцией (458) к некоторому моменту времени / <С 2 за счет запаса энергии, полученного в предыдущие моменты времени, вы ше (так как полоса прозрачности фильтра нижних частоте течением
времени сужается непрерывно от оо до 0), чем средняя мощность на выходе такого же фильтра с импульсной переходной функцией (454) и неизменными во времени параметрами в установившемся режиме. Эти мощности выравниваются только к моменту / s 2 , в то время как равенство потребляемых мощностей наблюдалось в момент / = 1.
При отрицательных s функция В (s, /) равна 0 для | s| > |
/ и воз |
||
растает с ростом s от нуля при js| |
= |
/ до максимального значения |
|
при s = 0. Пунктиром на рис. |
12, |
а приводится такое |
же се |
мейство кривых, только для ф =
На рис. 12, б приводится семейство кривых В (s, /), аналогичных рис. 12, а, но для случая а = 2. При сравнении графиков рис. 12, а и 12, б видим, что во время перехода а от 1 до 2 при неизменности прочих параметров площадь под кривой В (s, /) увеличивается. Увеличивается также дисперсия, «интервал» корреляции, а макси мум смещается в сторону отрицательных значений.
144
Физически это можно объяснить следующим образом. Энерге тические компоненты процесса, претерпевающие затухание, кор релируются с аналогичными и более мощными, пришедшими в си стему ранее.
На рис. 12, в и 12, г приводятся семейства функции В (s, t),
определенной, согласно (484), при ф = j и а = 1 (в), а = 2 (г),
в предположении, что остальные параметры, входящие в эту фор мулу, такие же, как и для семейств, приведенных на рис. 12, а и 12, б. Нетрудно убедиться из сопоставления соответствующих
графиков рис. 12, что с увеличением ср от 0 до ТС соответствующие
значения В (s, t) возрастают в пределах, ограниченных соответству ющей сплошной кривой и пунктирной линией, приведенной на этих
графиках. |
13 приводится семейство для функции В (s, t), опре |
||
На рис. |
|||
деленной согласно (484), при ср = |
т. е. для случая, когда ^3 (s, t)- |
||
— В (s, t) |
и |
AB (s, t) = 0, и при |
t = 0,5; 1; 2; 3; 5; 10; 20 для |
случая а = |
1 |
(рис. 13, а) и случая а — 2 (рис. 13, б). Из сравнения |
|
семейства этих кривых замечаем, что с увеличением t корреляцион ная связь между отдельными значениями процессов, взятыми на фиксированном временном расстоянии, быстро возрастает, так что «излом» при s = 0 становится практически незаметен. Этот рисунок еще раз подтверждает выводы, сделанные при анализе рис. 12.
Ю 3 — 6 6 2 |
145 |
|
|
На рис. 14 представлены семейства кривых для функции В (s, t) |
||||||||
(484) |
с h(t, |
t 4- s) (461) как |
функции t, причем |
Ад = |
1, t £ [О, 3], |
||||
s = |
0; |
0,1; |
0,5; |
1; 5 для а = |
1 |
(а) и для а — 2 |
(б). |
Сплошными |
|
линиями даны |
кривые для ф = 0, а пунктирными — для ф = |
JT |
|||||||
Из |
этих рисунков видно, что с |
ростом £ дисперсия, определяемая |
|||||||
как |
В (0, f), убывает при а = 1 |
быстрее, чем при а = |
2, что |
уже |
|||||
подтверждалось при анализе предыдущих рисунков. |
|
|
|||||||
|
На рис. |
15 представлено семейство кривых для функции В (s, t) |
|||||||
(483) с h (t, t + |
s) (461) в зависимости от изменения s; при этом А0 = |
||||||||
Рис. 15.
— 1, s £ I—3, 31, t — 0,5; 1; 2; 3, для а — 1, (а) и для а = 2 (б). Сплошными линиями указаны кривые для ф = 0, а пунктирными —
для ф = Так как при ф = АВ (s, t) = 0, т. е. В (s, t) = 5) (s, (),
то пунктирные линии этих рисунков совпадают с пунктирными ли ниями рис. 12.
Так как нелинейные преобразователи с косинусоидальными динамическими характеристиками на практике встречаются до вольно редко, а сигналы, соответствующие рассмотренным выше
случаям, непосредственно в радиотехнике |
не применяются, |
то |
|
эти случаи представляют |
интерес больше с теоретической точки |
||
зрения — они раскрывают |
физическую суть |
рассматриваемого |
не |
линейного преобразования и описывают некоторые наблюдаемые при этом энергетические трансформации.
3.Пусть в (465) функция N (t) определена согласно (466) при
U — 0, т. е. N (t) = о)н^, а v (t) — линейный стационарный гаус сов процесс, как и в п. 1. Этот случай соответствует анализу сигнала с постоянной амплитудой частоты сон, промодулированного по фазе гауссовым шумом.
Выражение для смещенной корреляционной функции в этом случае согласно (475) имеет вид
— |
А2 |
coHs + 2ф)]. |
(487) |
h(0) [eh{s) cos coHs -J- é~h(s) cos (2мн/ |
U 6
Корреляционная функция выражения (465) согласно (477) получается в виде
B(s, t) = A e- ft<0> |
_ 1] coso)H + \e~m |
- 1] |
X |
|
X cos {2<dJ + (oHs + 2cp)}. |
|
|
(488) |
|
Величина смещения в рассматриваемом случае равна |
|
|||
ДВ (s, і) — Alé~h(0) cos (сон/ + ф) cos (сон/ + |
соHs + |
ф)- |
(489) |
|
Хотя в этом примере процесс v (t) стационарный, однако нали чие функции N (і) делает смещение зависимым как от s, так и от L
Усредненная по времени |
||||
корреляционная функция в |
||||
этом случае согласно |
(292) |
|||
и (293) записывается в виде |
||||
I— |
|
$о(*) = |
|
|
р—ft(0)+A(s) co s (О S |
||||
2 |
е |
|
Ф 0, |
(490) |
при сон |
|
|||
а д |
= |
- і г * - Л(0Ѵ |
(5>- |
|
— 1] coscoHsnpHcoHTfcO. (491) |
||||
Если предположить, что процесс v (t) определен согласно (456) |
||||
и имеет корреляционную функцию (457) при щ Iw (1)1 = 1, то (488) примет вид
^2 |
_ J L 2L е—exist |
|
В (s, f) = |
е 2 {[е 2 |
— l]coscoH + |
_ — е—exist
+ [е 2 — 1] cos (2(öHt + a)Hs + 2ф)}, (492)
а смещенная корреляционная функция в этом случае определится выражением
|
д2 |
_ _ 2 _ |
е—*lsl |
cos coHs + |
|
|
|
|
|
53 (s, t) = ~Y ~e 2 [e2 |
|
|
|
|
|||
|
_ JL e—“lsl |
cos (2CÜb£ + CÜ„S + |
2ф)]. |
|
(493) |
|||
|
+ e 2 |
|
||||||
На рис. |
16 приводится семейство кривых для функции В (s, |
f) |
||||||
как функции s, определенной согласно |
(492), построенных |
при |
||||||
А0 — 1, ф = |
0, (он — 10, s£ |
[—1,5, |
1,5], |
t — 0,5; 1; |
1,5 и а |
= |
1. |
|
С ростом I s| эта функция приближается к нулю. По |
аргументу |
t |
||||||
она является периодической. |
|
|
|
|
|
|
t) |
|
На рис. |
17 приводится семейство кривых для функции 53 (s, |
|||||||
как функции s, определенной согласно (493) |
при Л0 = 1, ф |
= |
0,, |
|||||
Ю* |
|
|
|
|
|
|
14Т |
|
©в = 10; |
sgt—1,5, 1,51, t = 0,5; 1; 1,5 и о = |
1. С ростом |s | эта |
функция приближается к периодической, определенной согласно |
||
(489), с А (0) = -у. По аргументу t она является периодической при |
||
всех фиксированных значениях s. |
|
|
4. |
Пусть в (465) функция N (t) определена согласно (466) при |
|
U = 0, |
а V (t) — линейный нестационарный |
гауссов процесс, |
определенный, как в п. 2 .
Тогда выражение для смещенной корреляционной функции
согласно (475) принимает вид |
|
|
Ф (s, t) = ~Y ~exp {— ~ [А {t, f) + h(t + s,t + |
s)]| [ew + s) cos coHs + |
|
+ e~h(t,t+s) cos (2<aHt -f- coHs + |
2cp)]. |
(494) |
Корреляционная функция процесса |
(465) |
в этом случае имеет |
||
вид |
|
|
|
|
в (S, t) = |
exp I-----L[h (7, і) + А {t + |
s, t + |
s)]| {[ew -t+s) — 1] ^ |
|
m cos соHs + [e~h(t’t+s) — 1] cos (2coH/ -f- coHs + 2cp)}. |
(495) |
|||
Величина смещения при этом определяется формулой |
|
|||
AB (s, f) = |
Лоexp I---- i- [A (t, t) + A (t + |
s, t + |
s)]J cos (coH* + |
qp) SJ |
|
igScos (coH^ -J- соHs -f- qp). |
|
(496) |
|
Заметим, что выражения, входящие в (484) и (495), а также аналогичные им, полученные для стационарного случая, входящие в (480) и (488) при t >■ 0 и s £ (—оо, оо) равномерно по s, удовлетво ряют условиям -elft««<+s)V— 1 > 0 и g - lft(M+s)l — 1 < 0. Поэтому
148
с учетом того, что модуль суммы не превышает суммы модулей, из (494) и (495) получаем
^ ( s , /)I == -4оexp I---- [h(t, t) + h(t + s, f + s)]}ch/і(*, t + |
s) (497) |
|||
и |
|
|
|
|
IB(s, 0 | < 4 |
exp {---- [h(t, t) + |
h(t + s,t + s)]}sh\h{t,t + s)\ |
||
при / > |
0 H S ^ |
(— oo, oo). |
|
(498) |
|
|
|||
Из |
последнего неравенства, |
в частности, следует, |
что если |
|
h ( t , t |
+ s) |
0, то IВ (s, t) J 0. |
|
|
Предположим теперь, что процесс ѵ (t) получен в результате фильтрации белого гауссового шума 7?С-фильтром с перемен ными параметрами и определен согласно (460) с корреляционной функцией (461). В этом случае (s, t) и В (s, f) имеют сравнительно сложный вид. На рис. 18, а приводятся сечения функции В (s, t),
построенные согласно (495) для фиксированных значений |
t = |
0,5; |
|||
1; 1,5 при А0 = 1 , (і)н = |
10, |
s£ [—1,5, 1,5] и при |
а — 1, |
ср |
= 0. |
Это же семейство, но при а |
= |
2, приводится на рис. |
18, б. |
Сравни |
|
вая эти графики с графиками, приведенными на рис. 16, замечаем
следующее. На рис. 18 при отрицательных s, если |s| |
■< г1,функция |
В Cs, t) — 0, причем она быстро увеличивает размах |
колебаний по |
мере приближения s к нулю, а затем этот размах убывает медленно. На рис. 16 размах убывает симметрично с ростом | s |. На положитель ной полуоси s для стационарного случая убывание размаха более
быстрое, |
чем в нестационарном случае. |
построенные |
||
|
На рис. 19, а приводятся сечения функции В (s, f), |
|||
согласно |
выражению (495) для фиксированных значений s — 0; |
|||
|
2л; 20я, соответствующих «гребням» рис. 18 при Л0 = 1, сон = |
|||
= |
10; / С [0, 3] и при а = 1, ср = О.Это же семейство, но при а = |
2, |
||
приводится на рис. 39, б. Пунктиром приводятся «огибающие», |
||||
которые являются верхней и нижней границами значений В (s, |
t) |
|||
при соответствующих s, независимо от частоты сон. |
согласно |
(466), |
||
а |
5. |
Пусть в (465) функция N (t) определяется |
||
V (() — линейный стационарный гауссов процесс, |
полученный |
|||
149
