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

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

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

ховкой снизу, — для измерения С2. В центре отмечен расположен­ ный между ними интервал интегрирования в устройствах разде­ ления каналов решающего устройства приемника. При таких ве­ личинах сдвигов в установившемся режиме последний интервал будет располагаться в центре посылки, границы которой отмече­ ны на рис. 5.2 моментами sT и (s+1) Т (здесь предполагается, что защитный шнтервал АГ0<Л 7’).

Рис. 5.2. Временная диаграмма ра-

Рис. 5.3. Входной преобразователь УС по

боты УС по минимуму переход-

квадрату группового сигнала:

НЫХ помех

/ — входной сигнал; 2 — преобразованный сигнал

Амплитуды Ci и С2 свободных колебаний КФ измеряются ам­ плитудными детекторами АД4 и АД2. Разность Сi—С2 использует­ ся для управления фс. Если эта разность отрицательна, то фс уве­ личивается (например, путем добавления импульса в УС с дис­ кретным управлением), если же разность положительна, то фс уменьшается.

Возможны различные модификации УС рис. 5.1. Например, ис­ пользуя не только знак, но и величину разности Сi— Съ можно построить многопозиционное замкнутое УС с дискретным управ­ лением. На базе ВП рассмотренного УС возможно также построе­ ние разомкнутого УС. Сигналом, подаваемым на вход ВИРУ, в этом случае может служить последовательность отсчетов величин Сь измеряемых с различными сдвигами друг относительно друга.

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

зации.

Предложенный в [75] способ синхронизации по квадрату груп­ пового сигнала основан на том, что результат интегрирования на интервале (Д, ti-{-T0) переменной составляющей квадрата группо­ вого сигнала (или другого преобразования сигнала нелинейным устройством с четной характеристикой) равен нулю, если этот ин­ тервал оказался внутри посылки, и тем больше отличается от ну­ ля, чем большая часть соседней посылки попала в интервал ин­ тегрирования [56].

Квадрат 'грунтового 'сигнала, имеющего 'вид (5Л) (с тем, одна­ ко, отличием, что информация передается по всем каналам), со­ держит постоянную составляющую, равную, как показано в [56],

130

L

собой

0,5 £ a2i, и переменную составляющую, представляющую

г=1

 

сумму слагаемых типа

 

хц (t) = 0,5 ata, {cos [(t — j) Qt + <p( (t) — фу (/)1

 

cos [2 ©j + (i -f- j) Йt + ф,• (t) + Ф/(01-

(5.6)

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

Как и в рассмотренном выше УС по минимуму переходных по­ мех, для построения разомкнутого УС можно воспользоваться по­ следовательностью отсчетов интеграла, вычисленных на сдвину­ тых друг относительно друга интервалах. Для построения зам­ кнутого УС достаточно двух отсчетов, вычисляемых на интерва­ лах, сдвинутых, например, как показано на рис. 5.2. Разность ме­ жду этими отсчетами используется для формирования команд о добавлении или вычитании импульсов в замкнутом УС с дискрет­ ным управлением (рис. 5.3). Вместо (квадраторов Кв, показанных на рис. 5.3, можно использовать устройство для вычисления аб­ солютной величины. Интеграторы И4 и И2 работают со сдвигами в соответствии с диаграммой рис. 5.2. Полоса режекции ФВЧ дол­ жна быть весьма узкой, но достаточной для подавления «постоян­ ной составляющей» группового сигнала, которая может медленно изменяться под влиянием, например, замираний сигнала.

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

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

5*

131

Отметим два возможных направления построения адаптивных УС. Первое направление ориентировано на борьбу с искажения­ ми группового сигнала в целом. Адаптация в этом случае может заключаться в прекращении подстройки фс на те интервалы вре­ мени, в течение которых сигнал настолько искажен, что результа­ ты измерений фс не могут быть надежными. На этих интервалах времени синхронизм обеспечивается за счет стабильности ЗГ пе­ редатчика и приемника (см. § 3.5). Сигнал о прекращении под­ стройки фс вырабатывается специальным устройством измерения степени искажения группового сигнала, которое может «обслужи­ вать» не только УС, но и другие узлы приемника (такой измери­ тель можно построить, например, на основе предложенного в [31, 403] для систем с ФРМ алгоритма измерения отклонений приня­ той разности фаз от ближайшего возможного варианта передан­ ной .разности ф.аз).

Ясно, что указанное направление построения адаптивных УС применимо во всех рассмотренных УС при наличии измерителя степени искажения сигнала.

Второе направление построения адаптивных УС предполагает борьбу с искажениями отдельных участков спектра сигнала. Реа­ лизация этого направления требует, во-первых, наличия в прием­ нике устройств измерения степени искажения спектральных ком­ понент сигнала. Во-вторых, необходимо, чтобы входной преобра­ зователь, формирующий команды о подстройке фс, «учитывал» различные спектральные компоненты сигнала с различными весо­ выми коэффициентами. Такой учет возможен, если ВП состоит из частных измерителей, каждый из которых обладает известной ча­ стотной избирательностью, и сумматора. Возможен, например, следующий алгоритм адаптивного УС.

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

Рассмотренное ниже УС по модулю вектора сигнала [31, 34, 35] представляет, в частности, интерес, благодаря возможности использования его в качестве адаптивного.

Синхронизация по модулю вектора сигнала. Алгоритм синхро­

низации аппаратуры «МС-5»

основан

на

использовании величин

С, определенных ф-лами (5.5)

и (5.3)

с той разницей, что на вход

корреляторов, вычисляющих проекции *)

X, Y (или на вход ком-

1) Интеграл от произведения двух функций можно .рассматривать как их ска­ лярное произведение, т. е. как проекцию одной функции на другую [79, 56].

132

мутируемого согласованного фильтра) подается весь групповой сигнал и поэтому величины проекций Л', Y (или модуля двумер­ ного вектора с такими проекциями) зависят не только от переход­ ной помехи, но и от ^-канального сигнала, выделяемого корреля­ торами. Поэтому по величинам проекций X, Y или по величине модуля С нельзя судить об интенсивности переходных помех. От вредного, в этом смысле, влияния k-ro канального сигнала можно избавиться вычитанием двух модулей, измеренных со сдвигом, равным длительности посылки. Действительно, если интервал ин­ тегрирования попал внутрь посылки, то выполняются условия ор­ тогональности и величины С„ и Сп-и измеренные на n-й и (п1)-й посылках, обусловлены только амплитудой k-ro каналь­ ного сигнала и равны друг другу. Разность С„—Сп-ь при этом равна 0 .

Таким образом, характеристикой переходной помехи может

служить величина

JC„— С„_,|.

(5.7)

v =

Для измерения величины

модуля С в устройстве

рис. 5.4, как

и в УС рис. 5.1, использованы КФ и АД. Величины С преобразу-

Рис. 5.4. Вариант схемы вычисления величины V

ются в V с помощью элемента памяти (ЭП), запоминающего на одну посылку постоянное напряжение С, вычитающего устройства

иустройства АВ для вычисления абсолютной величины. Измеритель V обладает определенной частотной избиратель­

ностью.

Как следует из

приложения 2,

«энергетический вклад»

/-го канального сигнала в величину V примерно пропорционален

sin2{nX(i—k)]/(i—k)2, где

)J 0 — длительность части соседней

по­

сылки,

попавшей внутрь интервала интегрирования («вклад» k-ro

канала

пропорционален л 2К2/2). Заметно

подавляются только

ка­

налы, для которых 2h(ik)>A, и

 

 

если необходимо, чтобы на величи­

 

 

ну V

существенно

влияли только

 

 

по 4— 6 соседних с

k-yi канала (по

 

 

2—3 с обеих сторон), то следует

 

 

параметры УС выбирать так, что­

 

 

бы

в

установившемся

режиме

 

 

Л>0 ,'2 .

построения

замкнутого УС

 

 

Для

 

 

необходимо, аналогично предыду­

 

 

щим

алгоритмам,

воспользоваться

 

 

парой

 

измерителен

V, работающих

 

 

со сдвигом, как показано на рис. 5.2.

Для борьбы с селективными зами­ Рис. 5.5. Входной преобразователь раниями следует взять несколько УС по модулю вектора сигнала

133

таких inaip. Схема входного |гпрео1бра'3ователя УС, содержащего т ша|р измерителей V, показана и,а рте. 5.5. Заметим, что пары из­ мерителей могут различаться не только номером канала, на ко­ торый «настроен» К.Ф, но и входными сигналами, которые могут подаваться с разных ветвей разнесения.

5.2. Расчет УС по минимуму переходных помех

Коэффициенты уравнения ФоккераПланка. В УС с двузначным управлением по минимуму переходных помех (рис. 5.1) решение о добавлении импульса при­ нимается на основании сравнения двух модулей — Ci и С2, — определяемых проекциями переходной помехи (5.5). Найдем совместное распределение этих

величин.

Так как проекции переходной помехи распределены по нормальному закону, имеют нулевые средние (см. приложение 2) и разноименные проекции взаимонезавюсимы, то совместное распределение четырех 'проекций Х\, Хг, Yь Уг можно представить в виде

 

 

1

 

 

1

*1 + 0?

И * ь *г. Ук Уг) =

(2 я Oi а2)г (1 — г2)

ехр

+

 

2 (1 - г*)1

 

+

*2 + У2

г *1*2

Ух Уа

 

 

-----:----- 2

а1 ст2

 

Вычисляя

интеграл

а,

 

 

 

 

 

 

 

 

2я 2я

 

 

 

 

 

J

J Pi р2 w (Pi cos Oj, pacos02, pxSinGj, p2 sin 0S) d 0! d 02,

о

о

 

 

 

 

 

определяющий двумерное распределение величин Ci, С2 [87], видим, что это распределение

wc (Pi.

Pi Рг

• ехр

 

1

РГ

Р г

X

Р г ) =

— г2

 

2 ( 1 — Г=) I 02 +

а2

 

а?о| ( 1

 

 

 

 

 

 

 

 

' P

i

Р г

 

 

(5.8)

 

 

X

/

ощ Oa (1

—г2)

 

 

где 1о(х) — функция Бесселя мнимого аргумента.

 

что С2> С |, т. е.

Вероятность

добавления

импульса равна

вероятности того,

 

 

 

<зо ио

 

 

 

 

 

 

<7i= J j* wC(Px, Рг) dptdpi-

О р,

Для вычисления этого интеграла следует, во-первых, перейти к полярной системе координат

Щ/ 1

 

°2V '

= 2

Sina,

 

----------------------------------P i

= г cos a,

-

Р

 

 

 

и, во-вторых, заменить функцию Бесселя ее интегральным представлением [41, 149]. Поменяв затем порядок интегрирования и вычислив с помощью таблиц [41] внутренний интеграл (по г), получим

<7i =

n-

I

 

sin2 a du d a

1

(1 — г sin 2 a cos и)г

 

 

.

о, о

 

 

arctg

 

134

Этот двойной интеграл после преобразований также последовательно сво­ дится к табличным [45] сначала по и, а затем — по а. В результате находим

„2

J2

СТ2

0[

4i =

(5.9)

V [ о? + а |) 2 — 4г2 о2 о!

Отсюда в соответствии с (3.45а) получаем выражение для коэффициента а0(х) уравнения относительно плотности вероятности фс

2

2

° 2

°1

а0 (<р) =

(5.10)

У ( о2 + o f)2 — 4 г2 а2 of

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

Вид функций 0i (ф), 02(<р), г(<р) зависит от величины сдвига АТ между ин­

тервалами интегрирования, на которых измеряются переходные помехи (рис. 5.2). Оптимальная величина АТ, обеспечивающая, например, наибольший наклон функ­ ции ао(ср) в окрестности ф » 0 , зависит, в свою очередь, от уровня помех в ка­ нале связи и при слабых помехах близка к величине защитного интервала АТ0. Ограничимся в дальнейшем этим случаем, т. е. примем АТ=АТо.

Из приложения 2 следует, что моменты распределения проекций определя­ ются величинами Я|, Яг, представляющими собой нормированные к Т0 длитель­ ности последующих посылок, попавших в интервал интегрирования. На основа­ нии (П2.10) и (П2.17) можно записать

 

 

QTp

 

 

1

(<= 1 .2)

 

 

 

2

 

 

2ft2

 

 

 

 

 

 

(5.11)

 

 

Q T о

 

 

 

1 — Л

 

г 01 о2

kmin (1

kmin)

 

 

2

2ft2

 

 

 

 

 

 

 

где kmax

и k m i n — соответственно

большая и меньшая из двух

величин — Xi

и Яг, Q =0,5a27'0 — энергия сигнала

на

интервале

интегрирования;

ft2 — отноше­

ние этой энергии к спектральной плотности помех; Л^АГо/TY

 

Для

того чтобы

представить

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

уравнения Фоккера—Планка

в функции от фс, необходимо установить зависимости ЯДф) и Яг(ф). Установим их только для 0 ^ —ф ^ л , так как можно показать, что

Оо(ф) = — М — Ф). Ь ( Ф ) = Ь ( — ф)-

(5 - 12)

Фазу синхросигнала при АТ=АТ0 можно отсчитывать по положению

грани­

цы посылок относительно конца интервала измерения величины С2. Из рис. 5.6

видно, что по характеру зависимости Я(ф)

можно

различать два интервала зна­

чений фс — 0s£ —ф ^2яЛ /(1+ Л )

и 2яЛ /(1+Л ) С -^ р ^ я , — которым на рисунке

соответствуют два возможных значения фс — фа и

фь. Нетрудно проверить, что

при —я ^ ф ^ О

 

 

 

0

 

 

2 я Л '

 

(0<- Ф < 1 + Л ,

Ях — kmin

 

 

 

(5.13)

 

 

I 2 я Л

 

 

 

< — ф < я

 

 

U + А

Я 2 — kmax

(1 Ч* Л)

(0 <

— ф < Я ).

 

 

 

135

Подставляя |(5ЛЗ) .в (5.11), а затем полученные соотношения в '(5.10), на­ ходим с учетом (3.456)') выражения для коэффициентов а0(ф) и А(<р) :

- Ф а 0 - Фа ) Фл (1 л ) + А2

- 0 , 5

 

 

 

Ф л (! - Ф л ) +

 

2 — Л - 0 , 5

0 < ф < 1

2 я Л '

 

 

+ : Л2

+ А )•

 

(5.14)

а0 (ср) =

 

 

 

 

 

_ / Л ( 1 + Л -

2 ф л ) ( 1 - Л

+ -^ -)

05 ^(ФЛ

Л )(2’ +

 

 

 

.

2 — Л 1 — 0 .5

/ 2 я Л

\

+ Л - Ф л ) + (

1 - Ф л )(2Ф л -

Л) +

- ^

- |

( Т ^ Х < Ф < Я)

 

Ъ(ф) = 1

а\ (ф),

 

 

(5.15)

где фл =ф(1 +Л )/2я — отклонение синхросигнала от идеального положения, из­ меренное в долях длительности интервала интегрирования. Вид функции ао(ф) при двух значениях Л представлен на рис. 5.7.

Рис. 5.6. К определению зависимости

Рис. 5.7. Функция ао(ф)

Л(ф)

Формулы (5.12), (5.14), (5.15) определяют коэффициенты уравнения относи­ тельно плотности вероятности фс и на основании результатов гл. 3 позволяют найти основные характеристики рассматриваемого УС.

Моменты распределения фс в установившемся режиме. Для нахождения математического ожидания и дисперсии фс восполь­ зуемся ф-лами (3.43) и (3.44). Так как взятая с отрицательным 31наком гтроиз'водна'я функции flo(<p) ири ср=0|равна

V (1 + Л)

(5.16)

— Оо(°) V X (2 — Л)

’*)

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

136

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

2 я V А (2 -

Л ) .

N

(5.17)

Фо =

Л2 (1 +Л)

 

Так, при Л2 = 10, Л =0,2

А7= 100,

6Ш=4 0~ 3 величина ф0составит

0,3, т. е. примерно 5% длительности посылки, или четверть защит­ ного интервала.

Величина Ь(0), как видно из (5.15), равна 1, поэтому из (3.44) находим

,

2 л2 / Л (2 — Л)

(5.18)

а ч> —

Wft2 (l + Л)

 

При указанных выше условиях величина среднеквадратичного отклонения фс равна 0,3.

Формулы (5.17) и (5.18), на первый взгляд, противоречат ин­ туитивному представлению о поведении моментов распределения фс при Л-+0, когда разница между величинами С4 и С2 уменьша­ ется и, следовательно, уменьшается информация о направлении подстройки фс. Необходимо, однако, иметь в виду, что приближен­ ные ф-лы (3.43) и (3.44) справедливы лишь при условии, что практически все значения фс принадлежат интервалу, где приме­ нима линейная аппроксимация функции Ою(ф). Этот интервал в данном случае принципиально не может превосходить интервала

| q>| <2лЛ /(1+ Л ),

вне которого вид функции а0(<р) отличается от

использованного

при получении (5.16). Таким образом, при Л-»-0,

с одной стороны,

растет крутизна Оо(ф) в окрестности точки ф= 0,

с другой стороны, уменьшается ширина этой окрестности. Поэтому

ф-лы

(5.17)

и (5.18) тем точнее, чем больше Л и Л. На практике

Л > 0 , 1

-т -0 ,5 ,

и эти формулы можно считать точными при h> 1 .

Время достижения синхронизма. Определим время достижения

синхронизма

Sm как наибольшее (при разных начальных значе­

ниях фс) значение математического ожидания времени первого достижения области синхронизма. Для нахождения Sm можно во­ спользоваться приближенной ф-лой (3.94). Для этого необходимо указать границы области синхронизма (—фь ф0 и узловые точки

..., X j . В качестве ф! естественно принять фазу, соответствую­ щую сдвигу на половину длительности защитного интервала. Тог­

да совпадающая с ф4 начальная узловая точка

 

*о = ф1 = яЛ/(1 + д ).

(5.19)

В этой точке функция а0(ф) принимает значение

 

Ое(*о) = — (2

(5 20)

В качестве единственной |П1ромеж'уточ1ной узловой

точки возь­

мем

 

Xj = 2 k A/(1 + Л),

(5.21)

137

в которой изменяется вид функциональной зависимости Оо(ф)- В этой точке

 

 

 

 

Л + ^ ) ( 1 - Л +

2 — Л \ 1-0.5

 

1(*i) = — (1 — Л)

-

: Л Л2 ■ Г

(5.22)

Последняя узловая точка хг совпадает с границей <р3 опреде­

ленной в § 3.5. Так

как

производная

функция А0(ф) при <р = я

Л ; ( п ) = - ^ а » =

 

+ Л)2[(

 

1

 

 

- 0 . 5

(1

■ А + —

1

— Л + : ft2 ■)

 

 

 

 

 

(5.23)

 

 

 

 

 

 

 

 

 

и, кроме того, В (л) =4n2N~2, то на основании

(3.92) и (5Л4)

 

 

 

2 я

 

 

 

 

2 —Л

Х 2 — Фз — я

1

 

 

■/(1-A+v)(1-A+

 

 

ЛГЛ(1 +

Л)2

1/ V"

 

 

 

'

ft2

 

 

 

 

 

 

 

 

 

(5.24)

 

 

 

а0 (х0) «

я/TV.

 

 

(5.25)

Подставляя ф-лы

(5.19) —(5.25)

в

(3.94),

 

найдем выражение

для численного расчета времени достижения синхронизма. Если

при

.подстановке

учесть,

что

—х2 » я —хи

 

а0(хг) а0(хi) «

« \ aQ(Xi) |,

то это выражение примет вид

 

 

 

 

 

 

 

Sm =

 

N

 

 

.ln M*i)_|+

 

 

 

 

 

+

Л)

<h(*i) — а„(деь)|

а0(лсо)

 

 

 

1 —Л

 

2(1

 

 

+

In

N

 

 

2 я Л

 

 

 

 

 

 

 

l“o (*i)l

я | а0 (^M

+

Л

 

 

 

 

 

(5.26)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При отсутствии помех (h

-оо)

время достижения

синхронизма

 

S*

 

N

 

Л +

 

N

 

 

1 — л

(5.27)

 

2(1

+

Л)

(1 — A) In-----1- 2 яЛ

1 +

Л

 

 

 

 

 

 

 

 

 

 

По формуле

(5.26)

при А=100, А =0,2

на

рис.

5.8 построена

зависимость S m(h). Как видим, при изменении h от

1 до оо вели­

чина времени достижения синхронизма уменьшается более чем

втрое. Интересно сопоставить 'влияние

отдельных слагаемых

в

квадратных скобках в >(5.26) на

величину Sm. Для этого

«а

рис. 5.8 приведены три кривые Si,

S&

S3, соответствующие трем

указанным слагаемым. Из сопоставления кривых можно заклю­ чить, что влияние слагаемого Si незначительно и, как показывают численные расчеты, при /V>50 и небольших защитных интервалах (Л<0,5) им можно пренебречь. Тогда выражение для наибольше­ го времени достижения синхронизма принимает вид

Sm =

N

1 — Л In -

 

N

+

2 я Л X

 

2 (1

+ Л)

L{ |До (*i)l

п К

(*i)|

 

1 + Л

 

х / ( ‘ -

л + £ ) ( 1 -

л + 2- ^ ) -

(5.28)

 

где величина a0(Xi) определена ф-лой (5.22).

138

Рис. 5.8. Зависимость времени достижения синхро­ низма от h

Вероятность срыва синхронизма. Вероят­ ность срыва синхронизма РС1 можно найти

по ф-ласм (3.100), (3.97),

которые для УС с

дискретным управлением можно предста­

вить в форме

 

 

 

С 1

b ( 0 ) ~ V \ Г (0)1 г (Л) exp

|Ц«)|

 

 

 

N

 

л

 

 

 

 

 

 

(5.29)

где

I (<р)

=

j a (z) b 1 (2) dz.

 

 

 

 

 

о

 

 

 

 

Таким образом, для нахождения Рсi не­

обходимо

вычислить интеграл /(<р) при

Ф = я и

вторые производные этого

интегра­

ла при ср = я и ф= 0.

Учитывая, что а0(я) =

= а0(0) = 0, &(0 ) = 6 (я) =

1 и что

 

<РЦФ)

 

d а0(Ф) =

ао(ф) ь (Ф) ~ ь' <Ф)«»(Я>)

dфа

dф Ь(Ф)

 

Ь2 (Ф)

 

нетрудно видеть, что

 

 

 

 

 

I" (0) = а’0(0),

Г (л) = а’0(л).

(5.30)

 

Для

вычисления

/(я)

можно

воспользоваться какими-нибудь

методами численного интегрирования. Достаточно простую оцен­

ку порядка величины Рсi

можно дать, если учесть, что для УС

с двузначным управлением

| сво(ф) | < 1 и

МФ1 =

а° (-Т>— ~ до(ф),

ь (ф)

1 а\ (ф)

Так как функция Яо(ф) близка к «треугольной» с вершиной в точке ф= 2 яА (1 + А )-1, то и подынтегральная функция а0(ф)Ь~Чф) близка к «треугольной». Поэтому величину интеграла можно оце­ пить площадью треугольника

Пя)

 

2яА

^_1 / 2яЛ \ _

_я_

1д0 (*i) |

1

+ А

\ 1 + я /

2

1 — Oq (X j )

 

Выражение для оценки порядка величины Pci принимает тог­ да вид

Pei

I «о(0) I ао (rt) ехР

N_

1 «о (<i) 1

(5.31)

2

1 — a§(*i)

 

 

 

 

Расчеты

показывают, что при N = 100,

Л =0,2, И= \

порядок

вероятности срыва составляет 10- ®, а при h = 5 он составляет К)-1', т. е. Ра практически равна нулю.

139

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