книги из ГПНТБ / Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи
.pdfмного меньше B\ Полагая |
в уравнении |
(4.66) /г(іі)=0 |
и учитывая (4.68), находим |
|
|
а\Ч—"о— |
о„ — 2а |
(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 |
|||||
г. |
4Д 0 |
р. |
4Д |
|
|
|
|
где 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