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

книги из ГПНТБ / Марченко Б.Г. Метод стохастических интегральных представлений и его приложения в радиотехнике

.pdf
Скачиваний:
14
Добавлен:
24.10.2023
Размер:
9 Mб
Скачать

Выражение (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 (2+ 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

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