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

книги из ГПНТБ / Гинзбург, В. В. Теория синхронизации демодуляторов

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

ч

тожиые» ее изменения с падки зрения фс, так как изменения фс в - течение посылки весьма незначительны и одну посылку можно считать очень 'маленьким интервалом времени. Для нахождения ; вероятностных характеристик фс заметим, что вероятностные свя- ' зи между значениями %(t) (или величинами kt) распространяют­ ся на сравнительно небольшое число посылок (значительно мень­ шее N в УС с дискретным управлением). Так, если последователь­ но передаваемые информационные символы взаимно независимы, искажения сигнала в канале связи обусловлены сравнительно ши­ рокополосной помехой, а инерционность ВП не превосходит дли­ тельности посылки, то величины £(7) и £(7+Т) (или £г- и &i+i) ; практически независимы. Таким образом, интервал корреляции приращений фс значительно меньше постоянной времени УС и мо­ жно считать [27, 88, 119, 125, 133, 147], что плотность 'вероятности фс, подчиняющейся ур-иию (3.2), удовлетворяет уравнению Фок- кера-Планка (называемому также диффузионным уравнением и уравнением А. Н. Колмогорова)

£ - - £ [ * w -.1 + t £ [ * w . , ] .

(3-3)

где Wy =w v (х, s\y) — условная плотность 1ве(роятности значений

фс в точке

в момент времени s при условии, что <р (0) = г/;

А(х) и В(х)

— коэффициенты, представляющие собой математиче­

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

Решение ур-ния (3.3), являющееся плотностью вероятности, должно, очевидно, быть неотрицательным и удовлетворять условию нормировки, т. е.

П

(3.4)

%(х, S I у) > О, j шф(х, s\y)d x= 1,

—Я

 

и, кроме того, удовлетворять граничному и начальному условиям;

В>„(— л. s | у) = Шф(л, s\y)\ wv (x, 0\у) = б(х— у). (3.5)

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

В установившемся режиме плотность вероятности фс не зависит от начальных условий и неизменна во времени, т. е. юф(х, s|*/) =

==до ф ( Х ) , а ее производная по времени равна нулю. При этом из

уравнения в частных производных

(3.3) получаем

обыкновенное

дифференциальное уравнение второго порядка

 

-j- [А (х) % (*)]=

^ юфW1,

(3-6)

где доф (х) — не зависящая от начального значения фс плотность вероятности фс.

60

Коэффициенты уравнения для плотности вероятности фс. Как отмечалось выше, фс на нескольких посылках остается практиче­ ски неизменной в том смысле, что изменения статистических харак­ теристик приращения фс незначительны. Поэтому при изучении статистических характеристик приращений можно считать цепь об­ ратной 'связи в УС разомкнутой, а величину фс — фиксированной. При фиксированной фс процесс \(t) является периодически ста­ ционарным, причем его статистические характеристики зависят от значения ср. Для того чтобы подчеркнуть это обстоятельство бу­ дем вместо \(t) писать | ф(7). Таким образом, без учета расстрой­

ки тактовых частот приращение фс

H-AsT

Дф = J 1Ф(0 ^ .

а математическое ожидание приращения

< Л ф > = 5 < ц , (*)><#. i

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

<

Лф> = A s j <

(0 > dt.

 

 

 

0

 

 

Так как коэффициент Л(ф)

ур-ния (3.3) по определению равен

 

A (m) =

Пт < Д < Р >

 

 

Т

д з - о Д s

 

то с учетом расстройки тактовых частот

 

 

Л(ф) =

2л6ш+

ДДф),

(3.7)

т

— значение коэффициента А (ф)

при ра-

где Л0(ф) = J < 5Ф(0 ^ >

венстве тактовых частот, совпадающее с точностью до множителя 1/Г с постоянной составляющей математического ожидания

< £ ф(*)>■

Для устройства синхронизации с дискретным управлением, где ^■(i) представляется в виде последовательности 6-функций,

Л(ф) = ^ - а ( ф ) = y [б ш^ + а0 (ф)] ,

(3.7а)

гДе ао(ф) = <Л,->, т. е. равно математическому ожиданию числа Добавленных на посылке импульсов при условии, что фс равна ф.

61

Найдем теперь коэффициент 5(ф ), который по определению ра­

вен

В (ср) = Игл < Аф2>

As-*0 Д s

О

где Дф=Аф—<Д ф > — флуктуирующая часть приращения фс.. Из (3.1) имеем

t+AsT t+AsT

q

0

< Дф2 > — j

j

< 1 Ф(z)

(zi)> dztdz,

t

i

 

 

откуда после замены переменной z ченных пределов интегрирования во < т < tz-\-AsT на бесконечные1) — ка интегрирования находим

на x = z iz, изменения полу­ внутреннем интеграле t—z < о о < т < о о и изменения поряд­

оо t-{-As

<Дф2> = J j K\(z, x)dzdx,

--00 t

причем во внутреннем интеграле теперь интегрирование выполня­ ется по 2 .

Подынтегральное выражение в этом интеграле представляет со­ бой автокорреляционную функцию периодически стационарного процесса (t) и является периодическим по z. Поэтому, восполь­

зовавшись рассуждениями, аналогичными привлеченным при вы­ воде (3.7), можно записать

О

w

j

<Дф2>== Д s

J

| К(. (t, т) did т,

 

-оОО

откуда, принимая во внимание, что внутренний интеграл, совпадаю­ щий с точностью до коэффициента 1 с постоянной составляющей разложения автокорреляционной функции в ряд Фурье, является четкой функцией т, находим выражение для второго коэффициента ур-ния (3.3)

В(Ф)=

2 J j* /С£ (t, т) dtdx = j j Ki(t, т) dxdt,

(3.8)

 

о о

 

где Ki (t, т) = <

| ф (0 1ф(t + т)> .

 

*) Такое изменение не приводит к существенной погрешности, поскольку зна­ чения т, при которых автокорреляционная функция процесса §ф (t) существенно

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

62

Для УС с дискретным управлением из (3.8) можно получить,

что

В (ф) = (4n2//V2) Ь(ф),

(3.8а)

О

где Ъ(ф) = D (k) + 2 £ Kk{s),

D(k) — дисперсия числа добавленных в течение посылки импуль­ сов; Kk(s) — автокорреляционная функция числа импульсов, до­ бавленных на двух сдвинутых на s посылках (при фс, равной <р).

Если величины k на соседних посылках независимы, то

(3.86)

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

Пусть ф = 0 — некоторое наилучшее значение фс, при котором, например, обеспечивается наименьшая вероятность ошибки. Есте­ ственно потребовать, чтобы в установившемся режиме наиболее ве­ роятные значения фс группировались около нулевого значения, т. е. чтобы точка ф = 0 была точкой «притяжения», точкой устойчивого (в среднем) равновесия. Целесообразно также, чтобы в переход­ ном режиме изменение фс в направлении этой точки происходили по возможности быстрее.

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

в точке ф = 0 и иметь в ней корень, т. е.

 

 

' =

0

при ф =

О,

 

А (ф) > 0

при ф <

О,

(3.9)

. <

0

при ф>

0.

 

Для того чтобы точка ф = 0

была единственной точкой устойчи­

вого равновесия, функция А (ф)

не должна иметь других корней в

точках, где функция 'убывает.

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

(3.10)

т. е. А (ф) должна в точке ф = я иметь второй корень и возрастать. Из соображений симметрии естественно также потребовать, что­

бы А (ф) была нечетной функцией относительно точки ф = 0.

63

Функция В (<р), представляющая собой дисперсию приращений фс, очевидно, неотрицательна. Флуктуации этих приращений дол­

жны быть примерно одинаковы при положительных и отрицатель­

ных значениях фс. Поэтому функция В (ф)

близка к четной.

Вреальных условиях из-за искажений сигнала в канале связи,

атакже под влиянием аппаратурных ошибок функция А (<р) может несколько отклоняться от нечетной, а В (<р) — от четной. Однако эти отклонения должны быть сравнительно небольшими.

Для приближенного решения ур-ния (3.3) необходимо выбрать удобное приближенное представление функций .4 (<р) и В (<р). По­ скольку эти функции заданы на интервале (—л, я), естественно попытаться представить их отрезком ряда Фурье. При таком пред­ ставлении коэффициенты при нечетных членах разложения функ­ ции А (ф) «а основании изложенного выше должны существенно превосходить по абсолютным значениям коэффициенты при четных

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

В разложении функции В(ф) наибольшим из коэффициентов по абсолютной величине является постоянная составляющая. Четные

члены разложения должны превосходить

нечетные.

Представление функций Л(ф)

и В(ф)

в виде отрезков рядов

Фурье «физически» представляется наиболее оправданным. Вместе с тем, такое представление неудобно для нахождения плотности ве­

роятности фс в установившемся режиме. Действительно, представ­

ление А (ф)

и В (ф) с помощью ряда Фурье предполагает, что плот­

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

Обычно область наиболее вероятных значений фс сравнительно узка и ширина ее имеет порядок коэффициента передачи по петле обратной связи, например порядок 2n/N для УС с дискретным уп­ равлением. Поэтому требуемое для точного представления число гармоник в разложении плотности вероятности должно быть по­ рядка N, т. е. достигающим в отдельных УС нескольких сотен или даже тысяч. Поэтому представление плотности вероятности фс в виде ряда Фурье может оказаться громоздким и неудобным.

Так как в установившемся режиме большие значения плотности вероятности сконцентрированы в узкой окрестности точки ф = 0, то

вид функции

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

определяется поведением

функций А (ф)

и В (ф)

только в этой окрестности.

Функции А (ф)

и В(ф) в

окрестности точки ф=0

обычно хорошо

аппроксимиру­

ются степенными полиномами-

 

 

 

 

 

Л(ф)

= А0 — Лхф +

Аг ф2 А,ф*

(3.11)

 

 

В (ф ) = В0+ Bt ф +

В 2 фа.

(3.12)

64

Члены нечетного порядка выражения (3.11) взяты с отрица­ тельным знаком для того, чтобы коэффициент Л4 был положи­ тельным. Это придаст несколько более удобный вид полученным

вследующем параграфе формулам.

Всилу перечисленных условий, которым должны удовлетво­

рять функции А (ф)

и В(ф),

коэффициенты полиномов

(3.11),

(3.12) подчиняются соотношениям:

 

д > о ,

в0> 0; 14,/А |, IА / А I, I Bj/B0 \<zi.

(злз)

Если функции А

(ф)

и В (ф)

непрерывны вместе со своими про­

изводными соответствующих порядков, то коэффициенты Л„ В< можно определять с помощью разложения в ряд Тейлора. Иногда их удается определить с помощью разложения в ряд по ортого­ нальным полиномам или другим путем.

Для последующего важно отметить еще соотношение между

порядками малости коэффициентов Л(ф)

и В(ф).

Приращения фс за одну посылку, как указывалось выше, име­ ют порядок коэффициента передачи по петле обратной связи, что для УС с дискретным управлением составляет 2n/N. Коэффициент Л(ф), т. е. математическое ожидание приращений имеет, очевид­ но, тот же порядок (см. (3.7а)]. Коэффициент В(ф), представляю­ щий собой дисперсию приращений, по порядку величины совпа­ дает с квадратом коэффициента передачи [см. (3.8а)]. Таким об­ разом, коэффициент В ( ф ) — величина меньшего порядка по срав­ нению с Л (ф).

3.3. Распределение фазы синхросигнала в установившемся режиме

Решить ур-ние (3.3) в сколько-нибудь общем виде не удается, однако с его помощью можно найти все наиболее важные харак­ теристики УС. В данном параграфе на основе (3.3) определена плотность вероятности фс в установившемся режиме (при s-*~oo), которая для многих технических задач с исчерпывающей полно­ той характеризует УС.

Установившемуся режиму соответствует, как отмечалось выше, частный случай ур-ния (3.3) в виде линейного дифференциаль­ ного уравнения с переменными коэффициентами (3.6). Извест­ ное 'решение этого уравнения [88, 119, 125] в 'виде общего инте­ грала приведено несколько ниже. Такое решение, однако, не всег­ да удобно для инженерных расчетов. Найдем поэтому сначала

приближенное решение.

Будем искать

при­

Приближенное определение моментов фс.

ближенное

решение в виде отрезка ряда

Грамма-Шарлье

[88,

125]

 

 

 

 

 

 

(ЗЛ4)

3 -6 5

65

 

 

где Hk(x) — полином Эрмита 6-го порядка относительно веса ехр (—х2/2), а величины ф0, сгф, уз и у4 — соответственно мате­

матическое ожидание, дисперсия и коэффициенты асимметрии и эксцесса распределения фс.

Как известно, параметры %>, аф , уз и у4 выражаются через кумулянты и,- распределения, представляющие собой коэффици­ енты разложения в ряд Тейлора кумулянтной функции

Ч'ф0 и) = In 0ф (i со) » Ху i (0 +

-i- x,(i<o)*+

к3(i to)3 +

x4 (i (o)4,

 

 

 

(3.15)

где

 

 

 

Jl

exp^©*)©^*)^*:

(3.16)

0^(1©) = |

—Я

— характеристическая функция (преобразование Фурье плот­ ности вероятности).

В разложениях (3.14) и (3.15) слагаемые пятого и более высо­ ких порядков приняты равными нулю.

Связь между коэффициентами разложений (3.14) и (3.15) да­ ется соотношениями:

= Фо- *2 = стф- *з = *4 = • (3-17)

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

представление функций А

(ф)

и В (<р) ib виде

(3.11) и (3.12). Тог­

да, подставляя в ур-ние

(3.6)

полиномы (3.U)

и (3.12), умножая

его части на exp (icojc) и интегрируя их на интервале (—я, я), по­ лучим, учитывая правила дифференцирования оригинала и изо­ бражения преобразования Фурье,

А00ф (i со) — А 9; (i со) + Aj 0' (i со) А30ф" (i со) +

 

+ 0,51 © [В00Ф(1 ©) + в J 0; (1 ©) +

в20; (1 ©)] =

о.

Так как

0ф (ico) = exipTr(p(i©), то

это

уравнение

приводится к

уравнению относительно кумулянтной функции

 

 

 

А - a y ; о о) +

а { К

(i to)]2+

a to)} -

-

А {[y ; (i to)]3 +

ЗЧГф (i

со)

(i ©) +

4 7 (!©)} +

+ If-{ B Q+ B, ^ ;(i ©) + A pp;(i to)]2+ вг y ;( i ©)} = o.

Подставляя сюда выражение (3.15) для

i(ico) и приравни­

вая нулю суммы коэффициентов при одинаковых

степенях

(L©), получим следующую систему нелинейных

алгебраических

уравнений относительно фо, <тф , хз и Xk.

 

 

А — А ф0 + Агф2—АдфЗ + о ;(А — ЗА Фо) — А*з =

(3.18)

66

—стф(Л—2Л2(Ро+ ЗЛ3ф2) + х*(Л2—ЗЛ3ф0)—

— Л (ЗОф +

щ) + у

(Во +

Вгф0 +

В2 Фо +

В2СТф) =

0; (3.19)

- у (А- 24фо+ ЗА,ф§) + а;(Л2-

ЗЛзФо)+ у (Л- зл,ф0)-

—У А °1 *3 + у

[Оф(Si+ 2В2ф0)+ В2х3] = 0;

 

(3.20)

4(gА

2Л2ф0-JЗЛ3фд) -jо*х3(Л2

2Л3ф0)

 

 

- А, (2» .^ + о; +

Y »!) + Y [ т - <в*+ 2В* > +

в» ( « i

= 0.

 

 

 

 

 

 

(3-21)

В каждом из этих уравнений

в соответствии с (3.13)

основную

роль играют слагаемые с коэффициентами Л4 и В0. Поэтому при приближенном решении целесообразно определять математичес­ кое ожидание фо из ур-ния (3.18), дисперсию о£ из ур-ния (3.19),

а кумулянты третьего и четвертого порядков — соответственно из ур-ний (3.20) и (3.21).

Будем искать решение системы уравнений, используя метод последовательных приближений Ньютона, в соответствии с кото­ рым i-e приближение а» 'корня а функции f(x) можно получить из (i—1)-то приближения по формуле

а<= - [/ ( «f-i)]/ [/' ( <Vi)] ' (3-22)

В нулевом приближении все искомые величины положим рав­

ными

нулю. Тогда

из ур-ний (3.18)

и (3.19)

находим величины

фо и

в первом приближении

 

 

 

 

 

 

 

 

 

 

 

 

 

Фо =

Л0/Лх-,

 

 

 

 

 

 

(3.23)

 

 

 

 

о2 =

B J 2 A , .

 

 

 

 

 

 

(3.24)

Первое приближение для кумулянта третьего порядка х» будем искать с уче­

том первых приближений (3.23) и (3.24) для математического ожидания

и дис­

персии. Использование (3.22) дает

 

 

 

 

 

 

 

 

 

 

 

 

А*

о АоАз .

 

.

g

 

d®_

 

 

 

 

Bl

Al

А>

Be

 

 

Bp

At

 

 

(3.25)

 

2

4

?

 

А \А4 з

4

2

 

1

ВрА/

д

,

 

 

I

 

 

 

j

^

+ О

-I

I

О

4

" Г о

\ 4

?

 

 

 

 

4

?

А]

 

 

 

а \

2

 

В

силу (3.13),

а также учитывая,

что

 

малы

по

сравнению

с 4i,

знаме­

натель в последнем сомножителе близок к единице и коэффициент асимметрии

Vs

Bi

ВгА0

\

— + 2-----

(3.26)

 

Bp

B p A i

I

При 4 д = 4 г = Bi =0

коэффициент асимметрии,

как « следовало ожидать.

Равен нулю.

 

 

 

3*

67

 

 

Из ур-ния

(3.21) с такими же допущениями находим

 

 

 

И. = 0,76 ( BgM?) М з /Л т -В ,/А ,) .

 

(3.27)

Величина

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

эксцесса, определяемая с помощью

(3.17),

равна

 

yi =

Z(Aa/A l ) ( A 3/A 1- B 2/B0).

 

(3.28)

При Лз = Вг=0 коэффициент эксцесса равен нулю. Итак, если

в достаточно

большой окрестности точки

х = 0, точнее — в окрестности корня функции

А(х),

функция А(х)

линейна, а В(х) постоянна, то искомая плотность вероятности нор­

мальна [61, 119, 133].

Поправки, даваемые вторым приближением, невелики при выполнении усло­ вий (3.13). Так, выражения для математического ожидания и дисперсии во вто­

ром приближении после некоторых упрощений принимают вид:

 

 

Ао

л2

( .

АрА3

 

 

Вр

А%

АрА3

•Ро =

. _

I

л

+

 

л

+ .9

 

2Аг

*

 

(3.29)

 

^1

 

 

 

At

 

 

 

Л?

 

 

Во +

 

Л„

Ва

Л]

 

 

~ А3

Bi

 

 

 

Вг —— +

Л

 

+

 

 

 

 

 

 

 

Л1

 

 

 

 

 

(3.30)

 

 

 

 

 

 

А2

 

 

 

 

 

 

 

 

 

^0

л

,

 

ч в _

6в

Аз

 

 

Ai — 2Лj

л о

 

 

 

 

о

1

 

2 Г

2

6В °

Аг

 

 

 

 

 

 

А \

h

 

 

 

 

 

 

 

 

 

 

 

 

Формулы (3.23),

(3.24),

так же как и (3.25)— (3.30),

можно уточнить, заме­

нив (3.11), (3.12)

разложениями по степеням (ф—ф<и) в окрестности точки Ф01,

в качестве которой используется первое или второе приближение величины мате­ матического ожидания (3.23) или (3.29).

При очень малых значениях отношения В(ф)/Л(ф), например при очень больших N в УС с дискретным управлением, математическое ожидание фс мож­

но определять из условия

 

Л (фо) = 0, Л' (ф0) < 0.

(3.31)

Дисперсия при этом асимптотически равна

 

° Ф ------- В (ф0) /Л ' (фо).

(3.32)

Из ф-лы (3.31) следует, что необходимым условием существования уста­ новившегося (стационарного) распределения фс является существование корня функции Л(ф). Если же функция Л(ф) не меняет знака на интервале (—я,л), что может иметь место, например, при большой расстройке тактовых частот передатчика и приемника и сильных помехах или при слишком большой инер­ ционности УС, то установившегося распределения фс не существует.

Формулы (3.26) и (3.28) позволяют оценить близость распределения фс к нормальному. Такую оценку можно получить, например, найдя поправку к зна­

чению функции плотности вероятности

при х = 3 а ^ (т.

е. на

границе

области,

которой принадлежат примерно 99,7%

значений

фс),

обусловленную

отличием

от нуля коэффициента эксцесса уч-

 

как видно из

(3.14),

 

Относительная величина поправки составляет,

 

~Я4(3 )= -|у - 30= 74.1,25.

Для того чтобы относительная поправка не превышала 10%, необходимо, чтобы у‘<0,1- Допустимую величину коэффициента эксцесса можно выразить через дисперсию фс и оценить таким образом допустимую дисперсию, при кото­ рой распределение фс незначительно отличается от нормального.

На основании (3.28) и (3.24)

у 4 « б3/ОАхф — ВВ0)2. /

68

Сомножитель в скобках на практике меньше единицы. Поэтому можно счи-

О

тать, что -у4< 6сгф и для выполнения неравенства -у*< 0,1, достаточно, чтобы вы­

полнялось неравенство 6оф < 0,1, т. е.

 

аф < 0 ,1 ,

(3.33)

что составляет примерно 2% длительности посылки.

котором

Соотношение (3.33) с некоторым запасом отражает условие, при

распределение фс близко к нормальному. Поэтому можно считать, что фс рас­ пределена по нормальному закону, если среднеквадратичное отклонение фс не превосходит 0,1-ь0,2 рад, т. е. 2-М% длительности посылки. При этом по правилу «трех сигма», почти все значения фс будут находиться в интервале, составляющем 10<-20% длительности посылки. Для многих практических задач такую точность можно считать достаточной. Уточнение закона распределения фс получено ниже на основе общего интеграла ур-иия (3.6) [88, 119, 125].

Плотность вероятности фазы синхросигнала. Интегрируя левую и правую

части (3.6), находим

 

 

 

 

 

А (х) шф W + Сх = y

[В (х) шф (х)],

 

(3.34)

где Ci — произвольная постоянная.

 

 

 

 

Решение ур-ния (3.34) имеет вид [18]

 

X

 

 

 

exp (—

 

 

 

%,(*)

1+

2

]

(3.35)

boo"

jexp [A (z) d

 

 

*

о

 

 

 

 

 

 

где А (х) = 2 I dz; Cj — произвольная постоянная,

о B(z)

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

Для нахождения Ci заметим, что правая часть в (3.34) не содержит постоян­ ной составляющей. Поэтому постоянная Ci представляет собой среднее значение произведения — Л(х)шф (х), т. е.

я

Ci = — ^ j А (х) шф (дг) dx.

—я

Для получения приближенного значения Ci заменим в подынтегральном вы­ ражении А(х) полиномом (3.11), а плотность вероятности шф(х) — нормальной

плотностью вероятности. Заменим, кроме того, пределы интегрирования на бес­ конечные. Тогда

Ci =

^ [— А0 + Аг Фо — Аг ( о* + ф§) +

Л,ф0 ( Ф^ + За*)] .

Поставив

сюда приближенные выражения для

математического ожидания

и дисперсии из (С.23), (3.24), находим

 

Заметим, что если функция Л(ф) — нечетная, а В(ф) — четная, то

 

Сх = 0.

(3.36а)

69

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