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

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

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

путем фильтрации из белого шума. Этот случай соответствует случаю сигнала с постоянной амплитудой частоты wH, промодулированного

по фазе моногармоническим сигналом

частоты ш0 и амплитуды U,

а также стационарным гауссовым шумом.

 

 

 

 

 

 

 

 

Выражение для смещенной корреляционной функции согласно

(475) в этом случае принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

АІг\

—/і(0) I M s )

 

2U sin

)nS

.

/ ,

,

0)fjS

\

 

 

 

33 (s, t) =

о

 

 

 

+

 

 

е"~' cos

 

sin (co0^ 4----1—J — co„s

+I еh( s)

cos

2U cos CO0S

COS

-j

 

 

 

 

ti>Ht

+

2<p

 

.

(499)

 

 

 

 

 

I—j -f- coHs -f~ 2

 

 

 

 

 

 

 

 

 

 

 

 

Корреляционная функ­

 

 

 

 

 

 

 

 

 

ция выражения (465) соглас­

 

 

 

 

 

 

 

 

 

но (477)

получается

в виде

 

 

 

 

 

 

 

 

 

B(s, t) =

 

О

л(0))і ~h(s)

 

 

 

 

 

 

 

 

 

- f - e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1] cos

2(7 sin

(B„S

X

 

 

 

 

 

 

 

 

 

 

Xsin |^co0( 4-

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

■A(s)

-1] X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X cos

2U cos -M - X

 

 

 

 

 

 

 

 

 

 

X cos

 

+

M„s

+

 

 

 

 

 

 

 

 

 

-j"- coHs -j- 2сон/ -|—2ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(500)

 

 

 

 

 

 

 

 

 

Величина

 

смещения

 

 

 

 

 

 

 

 

 

корреляционной

функции

 

 

 

 

 

 

 

 

 

определяется согласно(478)

 

 

 

 

 

 

 

 

 

следующим образом:

 

 

АВ (s, () = Лс<з-Л<0> cos (сoHt +

U cos со0( +

cp) X

 

 

 

 

 

 

 

X

cos [сон (t -f- s) -(- U cos COQ(/ -f- s) - Ь

Ф I

 

 

 

(001)

Воспользовавшись выражением (499) с учетом

 

(292),

получим

усредненную

 

по

времени смещенную

корреляционную

функцию

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

__0_ g-A(0)+A(s) ^ 2 и s i n

 

j c o s Ö)hS

 

 

 

 

п р и

сон ф 0 ,

3 3 0 ( s ) =

АA l

 

е _- h т

|eA(s,/o ^

sin

 

_

e - H

s ) J o

( 2 U cos

 

 

 

(502)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

o)„ = 0,

150

где J0 (X) — функция Бесселя действительного аргумента нулевого порядка (бесселева функция первого рода). Усредненная по времени корреляционная функция согласно (500) и (293) получается в виде

е~щ0) [еЛ(!!>— 1] cos(oH J 0 | 2t/sin— -j

при

сон Ф 0,

В0 (s) = ( _ |L

е- т j[e*(», __ !j j o (2U sin i l l ) _

 

(503)

—' ie~h(s)— 1] J0(2U cos - ^ - j |

при

coH= 0.

Следует отметить, что обычно при модуляции на практике выполняется условие ©н to„, поэтому в (502) и (503) выражения,

приведенные в предположении, что сон Ф 0, представляются ин­ тересными как с теоретической, так и с практической точки зрения, а выражения, приведенные в предположении, что сон = 0, инте­ ресны с чисто теоретической точки зрения.

Ча рис. 20, а и б приводятся семейства усредненных по времени смещенных корреляционных функций (}30 (s), значения которых вы­

числялись по формуле (502) при А0

1,

0 =

1, <ÜH = 0,

ш0 = 0,

± 10, ±30, ±50, s £

[0; 0,3] и а =

1

(а), а

— 2 (б).

функций

На рис. 20, в и г

приводятся

аналогичные

семейства

$ 0(s), полученные путем расчета по формуле (502)

приЛ0=1, 0 = 0,5,

(он = 0, ш0 = 0,±10, ±30, ±50, s £ 10; 0,31 и а

= 1 (б), а = 2 (г).

151

Из этих графиков и выражения (502) видно, что при сон = 0 значения ^30 (s) удовлетворяют условию

А е -* 1'’■[еА(!% т - e- A(s,l< ?J0(S)< 4 - е~т ^

- e-hi%(2U)].

 

(504)

При этом график верхней границы (верхняя «огибающая» се­

мейства % (s) при о)н = 0) в (504) на всех рисунках

совпадает

с

(s), построенной для случая со0 = 0, т. е. при Ü)0 =

0 справа

в

(504) получено строгое равенство (график нижней границы — ниж­ няя «огибающая» — в (504) на всех рис. 20 выполнен пунктирной линией).

Из рисунков видно, что с уменьшением U, или с увеличением а, сужаются границы значений % (s), определенные соотношением (504), при условии неизменности остальных параметров, входящих в (504). Отсюда следует практический вывод: расширение частотной полосы шума в определенных пределах можно компенсировать со­ ответствующим увеличением индекса фазовой модуляции Абсцис­

сы точек пересечения

кривых семейства корреляционных функ­

ций (ß0 (s) при изменении а и U остаются неизменными

Заметим, что на рис.

20 приводятся значения

(s) только при

s >- 0, а при s <; 0 график получается симметричным отображением относительно оси ординат. В точке s = 0 функция ))3(| (s) недифферен­ цируема, т. е. соответствующий случайный процесс также не диф­ ференцируем.

152

Если ввести обозначения

 

 

x = sfB,

6/ = J a - ,

af = ~ - ,

 

Г

 

Iн

Ан

(457)

где J ---- 2^-----частота, измеренная в герцах, то из выражения

получаем

 

 

 

 

 

а выражение (502) при сон

0 может быть представлено в виде

ФоМ

=

Л2

 

 

(505)

- у - ехр [/„/І/ (х) — fHhf (0)] / 0[2Usin (лхб/)] cos 2ях.

На рис.

21 приводится семейство кривых для функции ^30 (s /н),

определенной

согласно (505),

при А0 = 1,

U = 1, /н = 100,

6/ =

= 0, ±0,05,

±0,1,

±0,15, sfH£ [0, 3] и а

= 2. Пунктиром приво­

дится семейство огибающих для соответствующих кривых.

 

Введем обозначение

 

 

 

 

А2

 

 

(506)

Фо(*) = —

exp lfHhf ( X ) fjif (0)] J0[2U sin (лхб/)].

Из выражения (505) видно, что

 

 

 

 

I Фо (*) I < Фо ДО-

 

 

Таким образом, функция

(х) является «огибающей» для ф,, (х)

График этой функции уже частично приводился на рис. 21 при по­

строении семейства кривых 3?0 (s/H).

 

 

 

для функции

На рис. 22 и 23 приводятся семейства кривых

Фо (S/J> определенной согласно (506),

при

Л0 =

1,

ос = 2, Щ= 0,

п.0,05,

±0,10,

±0,15, sfu £ [0, 12] и U =

1, /н =

100

(рис. 22, а),

U = 1,

/н = 1000

(рис. 22, б), и = 0,5,

/н = 100

(рис.

23, a), U =

= 0,5, /н = 1000 (рис. 23, б).

Сопоставление последних четырех семейств значений огибаю­ щих (sfH) при сі)н Ф 0 показывает, что с ростом несущей частоты скорость убывания огибающих уменьшается, а с убыванием U их «размах» становится меньше. Полученные выше корреляционные

153

функции и соответствующие им графики имеют непосредственное отношение к практике, так как характеризуют поведение смещенных корреляционных функций реальных технических систем, исполь­ зующих фазовую модуляцию в условиях гауссовой помехи, и могут быть применены при вычислении различных спектральных (на­ пример, выражения (324), (325) и (326)) и энергетических характе­ ристик такого сигнала.

Аналогичные графики можно построить и для корреляционной функции В0 (s), определенной согласно (503). По характеру они по-

 

 

добны

приведенным выше

 

 

для Д0 (s), только более

 

 

быстро

убывают с

ростом

 

 

s. Их отличия легко прос­

 

 

ледить

из

сопоставления

а

 

самих исходных выражений

 

(502) и (503).

 

 

 

 

 

 

 

 

 

 

На

этом

рассмотрение

 

 

примеров анализа сигналов

 

 

с фазовой модуляцией

за­

°'30 1 2 3 4 5 6

7 8 9 10 11 sfH

канчивается,

хотя

без

осо­

бых

затруднений,

исполь­

Рис.

23.

зуя

метод

интегральных

представлений, их рассмот-

 

 

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

Подобным образом можно получить выражение для корреляцион­ ной функции отклика вида (465) для случая частотной модуляции при наличии шума, если считать, что в (465) функции N (t) и v (t) определяются согласно (467) и (468) соответственно. При этом выражения (475), (477) и (478) внешне остаются неизменными, только в них N (t) и h (t, t + s) определяются с учетом (467) и (468). Однако окончательные выражения, рассмотренные в п. 1—5 настоящего параграфа, в этом случае будут несколько отличны­ ми. Для сравнения приведем выражение (s), аналогичное (502), но полученное для случая частотной модуляции:

■л,(0)+мцу

о

2

CO S (öHS

 

 

 

 

 

где

оо t

при сон Ф 0,

при сон = 0, (507)

hi (s) = И

ф (*)dx 1

Ф (у) dydt’

(508)

Оо

о

 

 

154

Ф (() — ядро линейного гауссового процесса, являющегося ад­ дитивной помехой при частотной модуляции, удовлетворяющее

условию J ф (x)dx £ Lz (0, оо).

о

В заключение этого параграфа отметим, что если ф — случай­ ная величина, имеющая равномерное распределение на интервале [—я, я] и независимая от v(f), то вычисление корреляционной функ­ ции от процесса (465) сводится к вычислению математического ожи­ дания от процесса вида (465). Этот случай вычисления корреляцион­ ной функции с применением методов, отличных от изложенных здесь, рассмотрен в [73, стр. 2111.

Отметим также, что в случае нестационарных процессов доволь­ но основательный их анализ в рамках энергетической теории мож­ но производить, пользуясь только корреляционными функциями, не переходя в частотную область и не применяя какого-либо из многочисленных определений «мгновенного спектра». Это наглядно видно из приведенного выше.

§ 5. Случай смешанной модуляции (амплитудная и угловая)

Этот параграф посвящен анализу в рамках корреляционной теории отклика, полученного на выходе устройства, осуществля­ ющего одновременно модуляцию по амплитуде и по фазе (АФМ) или по амплитуде и частоте (АЧМ) моногармонического несущего сигнала с помощью полезных неслучайных сигналов и гауссовых линейных, в общем случае нестационарных, процессов. При этом рассматривается случай зависимых гауссовых компонент, а это зна­ чит, что полоса пропускания фильтра в канале модуляции по ам­ плитуде и полоса пропускания фильтра в канале модуляции по фазе или частоте пересекаются (имеют общий участок).

Так как случай АЧМ с теоретической точки зрения не вносит ничего нового по сравнению с АФМ, то подробно будет анализи­ роваться только АФМ, а по поводу АЧМ ограничимся некоторыми замечаниями сравнительного характера.

Процесс на выходе устройства, осуществляющего одновремен­ но амплитудную и фазовую модуляции, может быть записан в виде

У (0 = [Л (t) + Ѵу(01cos [М (t) + ц2(t) + ф],

(509)

где А0 (t) — некоторая неслучайная действительная неотрицатель­ ная функция; N (() — некоторая ограниченная на любом конечном интервале неслучайная функция, которая для случая модуляции по фазе определяется согласно (466), а для соответствующего случая с частотной модуляцией — согласно (467); ѵг (/) и ѵ2 ( t) — ком­ поненты гауссового случайного вектора Ѵ2 (/), определенные со­ гласно (73), с матрицей вторых моментов, определенной выражением (364), при условии, что Му, (t)vz (t) Ф 0, так как в противном

155

случае задача может быть сведена к рассмотренной в предыдущем параграфе; ср — постоянная начальная фаза.

Остановимся сначала на определении математического ожи­ дания процесса (509). Пользуясь результатами § 6 гл. II, со­ вершаем переход от компонент вектора Ѵ2 (/) к компонентам век­

тора Ѵ2 (t) согласно соотношению (366), которое в рассматриваемом случае принимает вид

ѵЛі) = ап (() цг(/),

(510)

Щ(0 = «21 (0 «1 (0 + «22 (0 Ѵ2 (t)

при условии, что элементы матрицы А (t) = ||a,7(/)|f определяются согласно (369). С учетом (510) выражение (509) принимает вид

у (0 = [А0 (t) + аи (0 ѵг(0] cos [N (t) + a2l (t) vx(/) + a.,2(0 v2 (t) + ф].

(511)

Процесс (511) можно рассматривать как отклик некоторой нелинейной системы, имеющей два входа и характеризующейся изменяющейся во времени динамической характеристикой вида

Gt (xv х2) = [А0 (t) + ап (0 Ху\ cos [N (t) -f a21(t) хг + a22(t) xt + <p]

(512)

на воздействия

Хг (t) = yj (t) и x2(t) = v2 (t).

(513)

Так как процесс (511) удовлетворяет условию (352), то (512) можно представить рядом Фурье — Эрмита, коэффициенты которо­ го, согласно (356), получаются в виде

 

cmim2 (t) =

exp {— ~Y [«22 (t) + Ö2i (Ol) X

X

Л0(t) a%\ (t) cos

N (t) + ф + {tn1-f- m2) -~- + « n( 0 X

X

1(/) — a™iI+1 (0J sin ^ ( 0 + Ф + (Щ + т2) ~ ]. (514)

Из последнего выражения согласно (247) с учетом (376) и (379) получаем

Му(і) = с 00(0 =

= expf— -^-h22{t, t)\{A0{t) cos [N{t) + y\ — h21(t, t)sm [N (t) -f cp]}.

(515)

Откуда, согласно (363),

AB(s, t) = exp {---- L [fc22(t, t) + h22(t + s, t + s)]} 'A0(t)cos[N(t) +

156

+ ф] — h21 0 sin [N (t) + ф]} {А0 + s) cos [N (t +

s) -f- ф] —

— /і2і (t + s, t -j- s) sin \N (t -f- s) -f- ф]}.

(516)

Из выражений (515) и (516) при Л0 (і) == А0 и при условии неза­ висимости компонент, т. е. при h21 (t, t) = 0, получаем (470) и (478) соответственно.

Остановимся теперь на получении выражения для смещенной корреляционной функции. Пользуясь результатами § 6 гл II, совершаем переход от компонент вектора Ѵ4 (t, s), определенного

согласно (372), к ортонормированным компонентам вектора Ѵ4(t, s) согласно соотношения (384), которое в рассматриваемом случае

принимает вид (для краткости записи вместо Vj (t, s) везде пишем ѵ/)

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

ѵіФ

 

(t,

 

= агѵ13

 

 

 

 

1

 

 

 

 

_i_

 

 

(t +

s) = h

2 (t, s) a21 (t,

s) vx+

h22

(t, s) v2 =

+ b2v2, (517)

 

_L

 

 

 

 

 

 

 

 

j_

v2 (t) =

h *

(t,

s) ß31(t,

s) Oj + h21

(t, s) a32 (t,

s) v2+

(i, s) v3=

 

 

 

 

 

= h vi +

/ 2

^ 2 +

ftP»,

 

 

 

 

 

1

 

 

 

 

 

1

 

 

t»2(/ +

s) =

ft *

(/,

s) a41(f,

s) u4 +

7t22 (f,

s) a42 (/,

s) n2 +

1

 

 

 

 

 

1

 

 

 

 

 

+ Й33

(f,

s) a43 (f,

s) ü3 +

h 2

(t, s)vi =;d1ü1 + d2v2-Jr d3v3 + divi,

где элементы матрицы А (t, s) и М (^, s) определяются согласно (385) и (386) соответственно, а сокращенные обозначения а/, bj, fj, dj понятны из (517).

Для процесса (509) с учетом (517) выписываем выражение (388)

в виде

 

 

 

 

У (t) y(t + s) = [А0(0 + а л і И 0(/ +

s) +

M i +

Ь2оя] cos [N (f) + Ф +

+ fivi + h u 2 + /зуз1 cos [N (t + s) +

ф +

dlv1+

d2v2+ d3v3+

d4t>J.

 

 

 

 

(518)

Введем обозначения

 

 

 

 

N t = N(t) + N(t + s), NT = N(f) — N(t + s).

(519)

Процесс (518) можно рассматривать как отклик некоторой не­ линейной системы с динамической характеристикой

Gt,s ( д ? х3>х3, Х4) —•

= [Л0 (t) + ßjXj [А0 (t + s) 4- b1xl + b2x2] {cos [УѴ+ + 2cp +

157

 

+

(/i +

di) xi +

(/2 + d2) x2 +(dg + / 3 ) x3+ d4x4]-f-

 

 

+

cos [/V/

 

 

+

(fi d4) xl -f- (/2 d2) x2+ (/3

d3) x3d4x4]}(520)

на независимые нормированные воздействия

 

 

 

 

 

 

 

 

 

х4 (t) =

vx (t, s), x2 (i) = v2 (t, s),

 

 

 

 

 

 

 

 

x3(0 =

va (t, s),

xt (t) = v4 (t, s).

 

 

 

 

Смещенную корреляционную функцию процесса (509) будем

искать, исходя

из выражения

(389). Для

этого

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

приведем соотношение [7 ]

 

 

 

 

 

 

 

сю

 

 

 

 

 

 

 

 

 

 

 

 

j cos (ах +

ф) ехр ^---- j dx — К 2 я exp ^-------- ^ -|со вф .

(521)

 

—сю

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируя последнее соотношение по ф, а затем по а , по­

лучаем

 

 

 

 

 

 

 

 

 

 

 

 

I sin (ах +

ф) ехр

 

 

dx — ] / 2я ехр

 

 

(522)

 

 

 

 

 

 

 

 

 

 

 

 

 

сю

 

 

 

 

 

 

 

 

 

 

 

 

j X cos ( а х +

ф) ехр (---- dx =

У 2 л а ехр ^----------sin ф.

 

 

—сю

 

 

 

 

 

 

 

 

 

 

 

Дифференцируя второе соотношение из (522) сначала по ф, а

затем по а,

получаем

 

 

 

 

 

 

 

 

сю

 

 

 

 

 

 

 

 

 

 

 

 

j X sin (ах +

ф) ехр (---- —х

\ dx — V 2 n a ехр (---- cos ф,

 

 

сю

 

 

 

 

 

 

 

 

'

 

 

(523)

 

 

 

 

 

 

 

 

 

 

 

 

j

*2 COS ( a x +

ф) exp (-----dx =

Ѵ 2 л ( 1 a 2) exp (---------- cos ф.

Согласно (389) с учетом (356) имеем

 

 

 

 

 

 

 

 

 

 

53(8,0 =

 

 

 

 

— (2я) I

1

1

j*

(*i> х2>Хд, х4) dx4dx2dxgdx4 =

(524)

где

Gt,s (•) определяется согласно (520),

 

 

 

 

 

со

оо

оо

сю

 

 

 

 

 

 

 

h = §

]"

I

J

Мо (0 A if +

s) + [0іА>

"Ь s) +

b4A0(01 xx +

 

—0000—0000

 

 

 

 

 

 

 

 

“f"

"t“

^ 2 ^ 0

(^) -^2 “ h

 

C O S / 1 “f~

2<p -{ -

(/^ - f -

 

“ {*

+

(/2 + d2) x2+

(/3 + d3) x3+

d4x4]exp (---- Щ xfj dx4dx2dx^x4,

158

а выражение для определения / 2 получается из выражения для

/ 4

путем формальной замены всех

dj на —djy / = 1, 4 ,

и

N t +

2ср

на NT-

 

(523),

получаем

 

 

 

Воспользовавшись (521), (522) и

 

 

 

Л = (2л)2exp I---- [(/х + dx)2+

(/2+

4 )2+

(/з + 4 )2 4

dll} X

 

X {А0 (t) Ай (t 4 s) cos (Nt 4

2ф) — (fx 4

dj) [a4A0 (t 4

s) 4

 

4 M o (0] sin (Nt 4 2<p) + axbx [1 — (/x + dt)2]cos (Nt 4

2<p) —

 

b2 (f2+ d2) A0 (t) sin (N t +

2<p) — axb2 (fx + dy) (f2+

d2) x

 

X cos(A^ +2ф)}.

 

 

(525)

С учетом (387) имеем (для краткости аргументы опускаем)

 

(fi ± di)2+ (/2 ± d2)24 (/з ± da)24

= hn ati 4 h22a t 4 h334

 

+

Ьищі +

h22a42 + h33ah 4 hu ± 2 (/iu a31a41 4

^22a32a42 "Ь ^33а4з) —

 

 

 

= h22(t, t) 4

h22(t +

s, 1 4

s) ±

 

2/t22 (i, 1 4

s),

(526)

где знак «4 » относится к выражению для Іх,

а знак

«—»—к выра­

жению для

/2.

 

 

 

 

 

 

 

 

 

 

 

 

С учетом (385) и (386) получаем

 

 

 

 

 

 

 

 

 

(fi ± dj) 1Л„ (t 4

s) 4 M o (t)] 4 b2(f2±

d2) A0(i) —

 

 

 

 

1

 

 

1

 

 

 

1

 

 

 

 

 

 

=

 

( a 31 ±

a 4 l ) [^ 1 1 A 0 (^ +

S ) +

h - t

a 2 1 ^ 0 ( 0 1 +

 

 

J_

 

_L

 

 

 

 

 

 

 

 

 

 

4

^22

(^22a32 ±

^22 ß42) ^O (0 =

1^12 (0 0 І

^12 (0 ^ +

S)1 A0(t 4

S) +

 

 

 

4

[h21(0 ^ +

s) dz h12

(t 4

 

s. 1

4

s)] A0 (t).

(527)

Аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aА

(1 — (/i ± ^i)2] — М 2(/j 4

^i) (h + d2) =

 

 

 

 

 

 

 

 

 

 

'

1

_L

 

 

 

 

=

^11^21 fl

 

(П31 і

п41)21

/і (|

h22

(hw П31 zh

 

 

 

 

 

1

1

 

I

 

 

 

 

 

 

 

 

 

 

±

au ) (/i22 a32 ±

^22

a4a) =

Ьц (t, t 4

s) zb

 

 

zb \hX2(0 t) -F hl2 (t,

t 4

s)] [h21 (t 4

s, 1 4

 

s) 4

/г24 (t, t 4 s)].

(528)

Если теперь в выражениях (526), (527) и (528) взять верхние зна­ ки и подставить эти выражения в правую часть (525), то получим выражение для Іъ а если проделать то же, взяв нижние знаки, то получим выражение для / 2. Подставляя полученные таким образом выражения для Іхи / 2 в (524), получаем окончательный вид искомого выражения для смещенной корреляционной функции

(s, f) “ —- exp I---- 2~ [h22 (t, f) 4 b22 (t 4 s,14 s)]j X

159

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