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

книги из ГПНТБ / Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи

.pdf
Скачиваний:
21
Добавлен:
23.10.2023
Размер:
7.76 Mб
Скачать

много меньше B\ Полагая

в уравнении

(4.66) /г(іі)=0

и учитывая (4.68), находим

 

 

а—"о—

о„ —

(4.71)

Чтобы очертить границы применимости подобной замены, воспользуемся известным в статистической ра­ диотехнике результатом (см., например, [2]) : при воз­ действии белого шума на параллельный колебательный контур дисперсия шумового тока <^ в стационарном ре­ жиме определяется формулой

<=%

Из (4.71), (4.72) получаем условие малости шума n{t):

Подставляя (4.71) в (4.70) и имея в виду (4.63), окончательно получим

Х = К— '*sinx —Y,(0.

 

(4.73)

где

 

 

 

 

 

Ï , (<)=;? - *(* ) cos («* +

?„)

+

 

 

0

 

 

 

(4.74)

+ ^ ^ ( О С 0 5 К +

Ф0 + Ф(0].

 

Здесь уз — случайный

процесс

со

свойствами белого

шума

 

 

 

 

 

< Тз (t) > = 0. <

T, (t) Тз (' +

т) >

== % S (*). (4.75)

Соотношения (4.75) записаны на основе сопоставления процессов ya{t) и y%(t) (4.22) и учета статистических характеристик последнего

(4.32), (4.33).

2. Синхронизация генерато­ ра при наличии помех [6, 4].

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

Рис 4.4. обстоятельство позволяет пре-

140

небречь флюктуациями

синхронизирующего колебания

и считать его абсолютно

стабильным.

cos ®at

Пусть синхронизирующее

напряжение и с = Л с

подается в индуктивную

цепь

контура (рис. 4.4).

Необ­

ходимые обозначения физических величин указаны на схеме.

Для напряжения на контуре вместо (4.42)

имеем

или

LT\ +

Rt\ — "с = 4" j " (/ — -ц) dt

 

 

 

 

 

 

 

 

 

Если повторить

выкладки

 

(4.44) — (4.46)

'и учесть

обозначения

(4.47), то можно

получить

 

-

I

2

• / і

 

 

1

-2

 

•Чі + ^

Ч

^ ^ о і Л

1

-

- щ

\

 

 

+ ^ Ю - - ^ 8

І

П Ш

^

( 4 - 7 7 )

От соотношения (4.77) перейдем к уравнениям для амплитуды и фазы. Однако в отличие от случая § 4.3 этот переход производится относительно частоты син­ хронизации кос, а не через резонансную частоту контура

ЩІ(І)=А{І)СО$Ф(І),

Ф(*)=юс*+іфо+ф(*)-

Как и ранее, имеем

 

^ ( 0 =

- ° ^ (0 sin<5(f).

Поступая далее известным образом (§ 4.2, § 4.3), можно получить интересующее нас уравнение для флюктуирую­ щей фазы

? =

Д -

sin 9 +

НО cos

+ 90 + <Р (01.

(4- 78)

где

Д =

ш0 — ш с ;

Е=-Ц-Ас.

 

 

Исследования [6, 4] показывают, что синхронный ре­

жим работы генератора

возможен лишь при расстройках

 

 

д < - ! г = л -

 

С4-79)

14І

где Ac параметр, называемый шириной полосы син­ хронизации. С учетом (4.79) уравнение (4.78) перепи­ шем в виде

ç W - A s i t K p +

YiC).

(4.80)

где у* (Опроцесс со свойствами,

близкими к

белому

шуму.

 

 

3. Фазовая автоподстройка частоты под действием

шумов

[7]. Принципы действия

системы' фазовой

авто­

подстройки частоты (ФАП) поясним с помощью

рис. 4.5.

Гармонические

колебания

«о (О

основного генератора

 

f }

 

ОГ и uc(t)

синхронизи­

or

U ФД

Uytt}

руемого

 

генератора

СГ

подаются

на фазовый де­

 

Uottl

 

тектор

(ФД), на

выходе

 

Ucltl

которого

возникает

на­

 

 

 

пряжение

Uy(t),

завися­

 

cr

 

щее от разности фаз меж­

 

 

 

ду

uc(t)

 

и u0(t).

Напря­

 

 

 

жение Uy(i) поступает на

 

Рис. 4.5.

фильтр

 

нижних

частот

 

 

 

(ФНЧ),

выходное

напря­

жение с которого u(t) воздействуют на реактивную лампу (РЛ), которая изменяет частоту синхронизируемого ге­ нератора, приводя ее к совпадению с частотой основного генератора. Для практики особенно интересен случай, когда сигнал u0(t) поступает на фазовый детектор вместе с шумом 1(0- Это обычно имеет место при пространст­ венном разнесении основного генератора и остальных элементов ФАП.

Рассмотрим упрощенный

вариант

системы ФАП,

в которой фазовый детектор

заменен

перемножающим

устройством, а фильтр нижних частот обладает идеаль­ ной характеристикой (коэффициент передачи равен еди­ нице в области нижних частот и нулю при высоких частотах). При этом будем учитывать флюктуации ча­ стот основного и синхронизируемого генераторов. Пред­ положим также, что шум £ (0 можно представить в ви­ де узкополосного квазигармонического колебания со слу­ чайной амплитудой Е(і) и фазой Ѳ(0

•1(0 = £ ( O c o s 0 ( O = £ ( 0 c o s M + 6(0],

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

Іі^ (т) = a'R (т) COS аус.

142

Пренебрегая амплитудными флюктуациями, колебания "оСО и и с ( 0 можно записать в виде

« o ( 0 = ^ i c o s © i ( 0 = ^ i c o s [ c ù 0 f + 8i(0];

(4.81)

/ic(O=^2Sin02 (O = i 4 2 s i n H + e2 (0].

(4.82)

Здесь At, Аг — постоянные амплитуды, шо и <в— средние частоты, Ѳі(0 и 02(0—случайные фазы колебаний.

Запишем сумму напряжений основного генератора и - шума таким образом:

СО + 1 =А І cos Ѳ І (/) + Ес (0cos Ѳ,С0 —

 

— £ a ( O s i n 0 i ( O ,

 

где £ c ( 0 = £(*)cos(0—Ѳі), Es(,t) =E(4)sm{&—0i).

Мож­

но показать [8], что синусная És(t) и косинусная

Ec(t)

составляющие огибающей E(t) имеют нулевое матема­ тическое ожидание и распределены по нормальному

закону, а их функции

автокорреляции ^авт(т) и взаим­

ной корреляции kB3(x)

определяются соотношениями:

kam(x)=<Ea(t)Es(t+x)>

 

=

= <Ec(t)Ec{t+x)>

=

ozR{x);

kB3(x)

= <Es(t)Ec(t+x)>

=

=

<Ee(t)Ee(i+x)>

 

= 0.

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

щ, (t) = р ["о ( 0 + * (t)] ис (t) = 4

« а + £ =) t s i n X

X (Ѳ2 - Ѳ,) + sin (Ѳ, + Ѳ2)] - Es [cos (Ѳ2 - Ѳ.) -

- c o s f ö + e,)]}.

Фильтр низких частот пропустит лишь колебания разностной частоты, так что на вход реактивной лампы поступает напряжение

и (*) = - L цЛ, [(Лj + Е0) sin ? - Es cos.?],

(4.83)

<р = e2 —0i.

(4.84)

143

Пусть реактивная лампа управляет изменением частоты синхронизируемого генератора по линейному закону

 

Cù='CûcO—su(t),

 

где CÙCO средняя

частота

синхронизируемого

генерато­

ра при u(t)=0; s — крутизна линейного участка

модуля­

ционной характеристики реактивной лампы.

 

Из соотношений

(4.81),

(4.82), (4.84) получим

 

т = ѳ ,

-

ѳ , = к „ - «•>„) -

su (О-)- (Ѳ2

- Ѳ,),

 

или

 

 

 

 

 

 

 

 

 

 

 

<р =

Д0 --Д sin<p — -^-à(Ecsm

 

f —

 

 

Escos<p)-j-(4.85)

где

 

 

 

 

 

 

 

 

 

 

 

д о — шсо — <Ѵд

~2

psA1A2,ty

=

i)3 — 6,, Д0 — начальная

расстройка генераторов, Д—полоса

удержания.

 

Для того чтобы к уравнению

(4.85)

можно было

при­

менить

аппарат

марковских

процессов,

необходимо

удовлетворить условию [7]

 

 

 

 

 

 

 

 

H-s/U (Л, + £с )

^

I

Кф

 

к

*

 

и X . в р е м я

 

 

Л-

 

 

t(t) и

где

корреляции

для

процессов

<j> (t) соответственно. Неравенство (4.86),

характеризующее

флюктуации

£('0

и

ty(t)

как

быстрые,

справедливо

в ряде практических случаев [7]. С

учетом (4.86)

урав­

нение (4.85) можно записать в виде

 

 

 

 

 

 

 

? = ; д о - : Д 8 І п < р

+ Т6 (0.

 

(4.87)

где уб(0 —процесс, который может быть заменен экви­ валентным белым шумом.

Сопоставляя результаты анализа трех приведенных задач [см. уравнения (4.73), (4.80), (4.87)], приходим к выводу о том, что поведение фазы колебательной си­ стемы во всех случаях подчиняется одному и тому же стохастическому уравнению. Обращаясь для определен­ ности к уравнению (4.87) и полагая

<Y 5 Y6 (t + 1 ) > = 4 ^ 8 W = H T S M'

"14*4

запишем коэффициешы сноса и диффузии

 

 

 

Я, (?) = Д0

A sin ср, К2 (? ) = А .

(4.88)

Уравнение

Фоккера—I Іланка—Колмогорова для стацио­

нарного состояния (w(<?) = 0)

согласно

(4.88) имеет вид

 

д

rdw (у)

~(D0-Dsm<?)w(?)

=

0,

(4.89)

 

а

г с Of

г.

0

р.

 

 

 

 

где D. =

-j^-\ D

 

 

 

 

 

Вследствие

периодичности

коэффициентов

сноса

по <р

плотность

w(<p) будет также

периодична, т.

е. функция

w(y+2nk)

будет удовлетворять уравнению

(4.89)

при

любом k. При отсутствии граничных условий нельзя указать такого значения k, для которого в стационарном режиме поток вероятности был бы равен нулю. Поэтому '

стационарное решение уравнения (4.89) следует

искать

из условия G— const, а не G = 0, как это бывает

в боль­

шинстве случаев {см. выражение (3.129)].

 

 

Общее решение (4.89) записывается в

виде

 

ч>

 

 

w (<р) = С, ехр (£0<f -f- D cos <р) j exp (—

D cos f') df,

c,

 

(4.9Q)

где Ci и C2 постоянные интегрирования,

 

определяемые

из двух условий:

 

 

1) периодичности — w (<р + 2я&) = ai(q>) ;

2) нормировки на каждом периоде — j * tw(<p)d<p=l.

о

Вычисления [6,4] приводят к следующему результату:

 

 

9+2*

o»(<f>) =

-ç-exp(D0 -|-Z)cos<p)

J exp (—£>0?'—Z) cos ?')

 

 

 

(4.91)

где С =

tâ<T*Do

I / . О о (D)|2 ;

/ Ш о (D) —функция Бесселя мни­

мого индекса и мнимого аргумента.

При нулевой расстройке (До—0) плотность вероятно­

сти для

разности

фаз определяется выражением

10—186

145

и имеет симметричную формулу (рис. 4.6). Если шумы велики (Do мало),' то функция до(<р) близка к равно­ мерной.

Рис. 4.6.

Различные аспекты функционирования системы ФАП под воздействием шумов с помощью теории марковских процессов анализируются, например, в работах В. И. Ти­ хонова (7, 9—12], Э. Витерби [13].

4.5.Воздействие шума на детектор

сэкспоненциальной характеристикой

Разберем задачу [4] по определению плотности рас­ пределения напряжения на выходе экспоненциального детектора, состоящего из диода Д и параллельной

и

Mr J с " : Ü

Рис. 4.7.

цепочки RC (рис. 4.7). Характеристика нелинейного эле­ мента i=if\(u) выбрана экспоненциальной для того, что-

146

бы обойти трудности, связанные с использованием -не­ линейных зависимостей, не имеющих непрерывных про­ изводных [14—17].

Составим дифференциальное уравнение для выход­ ного напряжения г)(/). Имея в виду обозначения, ука­

занные на рис. 4.7, запишем

 

 

 

 

 

 

 

 

 

і = і1-\-і1,

 

У —

 

 

^^dt^LR,

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ ( 6 - , ) = СІ + - | - .

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

^ + Ж ^ І Г ^ - 7

! ) -

 

 

( 4 - 9 3 )

 

 

 

 

 

 

 

 

 

 

 

 

f

(Ç — 7j) =

ô e 0

( E _ , l )

;

a, b =

const.

(4.94)

Удобно ввести новую

переменную

 

 

 

 

 

 

 

 

 

 

х =

еа , ) ,

 

 

 

 

(4.95)

для которой уравнение (4.93)

примет вид

 

 

 

 

 

RCx +

л In X =

aWîe**.

 

(4.96)

Введем

обозначения:

 

 

 

 

 

 

 

 

 

 

 

 

a=\IRC,

 

k=>aabR,

 

 

(4.97)

 

 

x(f) = e o t - m ,

m = < e a

* >

*

(4.98)

Когда

\(t)—нормальный

процесс

с

нулевым

средним

значением,

величину

m легко

определить,

основываясь

на выражении для характеристической

функции

 

 

 

 

< е " ' * > = е

 

2

,

 

 

(4.99)

где и — вспомогательный

параметр,

а2 — дисперсия ве­

личины g.

 

 

 

 

 

 

 

 

 

 

 

Полагая

в (4.99)

и = —га,

 

 

 

 

 

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/п = < е ^ > =

е 2

.

 

 

 

С учетом (4.97),

(4.98)

уравнение

(4.96)

преобразуем

к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

— axlnx-\-mk-\-kx,(t).-

 

 

(4.100)

10*

147

При малом времени

корреляции

т к случайное воздейст­

 

вие K(t)

можно

представить как белый шум с

нулевым

 

средним и функцией

корреляции

 

 

 

 

 

 

< х ( 0 н ( * + т ) > = - ^ 8 ( - « ) ,

 

 

где

спектральная

плотность Na% определяется, как обычно,

 

по

соотношению

(3.61).

 

 

 

 

 

 

Стохастическому

уравнению

(4.100)

соответствует

 

уравнение Фоккера — Планка — Колмогорова

 

 

— 5 Г = w \ { m k - a x \ n x ) w { x ) } - \ — Ä

^ L -

(4.101)

 

Отыскание нестационарного решения (4.101) затрудни­

 

тельно. Что касается одномерного стационарного рас­

 

пределения, то оно легко находится по формуле

(3.132).

 

(jc) = C e x p | - - A - [ axa(lnx-±)-2mkx]}>

 

(4.102)

 

где С — нормировочная

постоянная.

 

 

 

 

От распределения

(4.102) необходимо

перейти к плот­

 

ности распределения для выходного напряжения ц, ко­

 

торое связано с x безынерционным нелинейным

преобра­

 

зованием (4.95). Используя известную формулу

 

 

окончательно получаем

 

 

 

 

 

W fa) =

Ca ехр J ~ -щ-

[ае 2 0 1 1 ~~т)-

2 m k &

a n ]

} •

4.6.Анализ работы накопителя

Внастоящем и следующем параграфах рассматри­

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

Разберем в этой связи [18] работу простейшего #С-накопителя, на выходе которого имеется пороговое устройство (рис. 4.8). Предполагается, что с момента ^о=0 на 'вход накопителя подается детерминированный сигнал s(rf) и белый шум n(t). На выходе фильтра на­ пряжение складывается из детерминированной b(t) и случайной x (t) составляющих и сравнивается в порого-

148

вом устройстве с напряжением Е. В момент равенства величин b(t) +x(t) —Е пороговое устройство срабаты­ вает. Задача состоит в определении вероятности сраба­ тывания порогового устройства за время t.

-CZ>

 

 

R

 

Пороговое

s(t)+n(t)

ßltl+xftj

истоойст-

 

во

 

 

Рис. 4.8.

Поставленная задача эквивалентна задаче о вычис­ лении вероятности появления выброса случайного про­ цесса x(t) над переменным порогом Un(t)=E—b(t). При отсутствии сигнала порог становится постоянным: uu(t)=E. Согласно (3.157), (3.149) вероятность Р(х0, t) превышения за время t случайным напряжением x(t) уровня ua(t) находится по формуле

 

 

 

 

 

 

«„О

 

 

P(x0,t)

= l -Q(x0,t)

=

\ -

j

q{x0,x,t)dx.

 

(4.103)

 

 

 

 

 

 

—oo

 

 

При

отсутствии поглощающей

границы

в точке х =

= un(t)

плотность

распределения

вероятностей

выходно­

го процесса

x(t)

характеризуется

нормированной плот­

ностью

вероятности

w(xQ,

х,

і),

которая

подчиняется

уравнению

Фоккера — Планка — Колмогорова

(3.55).

Если же существует поглощающая граница, то процесс

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

плотностью

q[xo, х, t),

удовлетворяющей тому же уравнению

(3.55),

 

dq

д

[X q}-

aW0

Ô*q

(4.104)

 

дх

 

 

 

 

 

но с граничными

условиями

типа (3.148):

 

 

д{хо,

ип,

-0=0.

 

(4.105)

Для общности будем полагать, что в начальный мо­ мент to на емкости - С существует начальное напряже­ ние хо, так что

q(x0)

X, U)

=Û(X—XQ).

Для решения сформулированной задачи воспользуем­ ся известным в математической физике методом- [19],

149

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