автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Разработка алгоритмов анализа систем синхронизации при воздействии гармонических и шумовых помех

кандидата технических наук
Рязанова, Мария Алексеевна
город
Москва
год
2011
специальность ВАК РФ
05.13.01
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Разработка алгоритмов анализа систем синхронизации при воздействии гармонических и шумовых помех»

Автореферат диссертации по теме "Разработка алгоритмов анализа систем синхронизации при воздействии гармонических и шумовых помех"

На правах рукописи УДК 621.396

Рязанова Мария Алексеевна

РАЗРАБОТКА АЛГОРИТМОВ АНАЛИЗА СИСТЕМ СИНХРОНИЗАЦИИ ПРИ ВОЗДЕЙСТВИИ ГАРМОНИЧЕСКИХ И ШУМОВЫХ ПОМЕХ

Специальность: 05.13.01 - Системный анализ, управление и обработка информации (в технических системах).

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата технических наук

МОСКВА. 2011

2 8 АПР 2011

4844598

Работа выполнена в Московском государственном техническом университете им. Н.Э. Баумана

Научный руководитель -

Научный консультант:

Шахтарин Борис Ильич,

доктор технических наук, профессор, заслуженный деятель науки и техники РФ, лауреат Государственной премии СССР к.т.н., доц. Чернышев Ю.Н.

Официальные оппоненты: д.т.н., проф.

Медведев Владимир Степанович к.ф.-м.н., Чернояров Олег Вячеславович

Ведущая организация: ФГУП «НПП «Дельта»

Защита состоится «17» мая 2011 г. в 14:30 часов на заседании диссертационного Совета Д 212.141.02 при Московском государственном техническом университете им. Н.Э. Баумана, по адресу: 107005, Москва, 2-ая Бауманская, д.5

С диссертацией можно ознакомиться в библиотеке МГТУ им. Н.Э. Баумана Автореферат разослан «11» апреля 2011 г.

Учёный секретарь диссертационного Совета

к.т.н., доц. Муратов И.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Современные системы и устройства управления, навигации и радиолокации, а также информационно-измерительные комплексы широко применяют системы синхронизации (СС). Эти системы выполняют многочисленные задачи; к числу наиболее распространенных из них относятся: синхронизация и демодуляция двоичных символов цифровой информации, измерение частоты и фазы сигналов, тактовая синхронизация, синтез сложных радиотехнических сигналов, синтез сетки высокостабильных частот, стабилизация частот генераторов различных диапазонов, когерентная демодуляция аналоговых и цифровых сигналов с частотной и фазовой модуляцией.

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

Особое значение при проектировании и создании СС имеют проблемы оценки устойчивости этих систем при воздействии помех различного рода. При этом статистические моменты фазовой и частотной ошибок слежения (среднее и среднеквадратическое значения) не дают полной информации о поведении СС. Поскольку СС - существенно нелинейная система, то в ряде случаев необходимо знание плотностей распределения вероятностей (IIPB) переменных состояния. Особенностью СС с рядом других систем (не фазовых) является существование множества устойчивых состояний равновесия, а в отдельных предельных случаях и устойчивых периодических движений 1-го и 2-го рода, что еще более усложняет картину при действии шумов.

Исследования СС существенно усложняются, если наряду с шумовой (флуктуационной) помехой на СС воздействует узкополосная негауссовская помеха. Здесь особенно важным является исследование СС при воздействии гармонических помех. К таким помехам относятся помехи соседних каналов в системах радиосвязи, помехи многолучевости, помехи космическим навигационным системам (GPS и ГЛОНАСС), организованные помехи радиолокационным станциям.

Основным методом исследования СС при наличии гауссовского шума на входе в настоящее время является метод марковских случайных процессов. Пионерами использования этих методов применительно к СС являются Р.Л.Стратонович и В.И.Тихонов. Существенный вклад в развитие теории синхронизации при наличии гауссовского шума на входе сделали Б.И. Шахтарин, В. Линдсей, А. Витерби, Дж. Холмс, В.Д. Шалфеев, H.H. Удалов, В.Н. Белых, В.Н. Кулешов, В.Д. Разевиг, В.В. Шахгильдян, А. Вайнберг и другие. В результате к настоящему времени теория непрерывных СС при наличии шумовых воздействий в основном завершена. В то же время теория дискретных систем синхронизации (ДСС) даже при шумовых воздействиях

разработана весьма незначительно. И если в непрерывных СС достаточно подробно методом марковских процессов исследованы системы 1, 2-го порядков, применительно к ДСС - это в основном касается систем только 1-го порядка. Такие исследования выполнены М.И. Жодзишским, В.Н. Кулешовым, В.В. Шахгильдяном, Б.И. Шахтариным, В.Н. Белых, В.П. Сизовым, Дж. Холмсом, X. Осборн, С. Гуптой.

Исследованиям дискретных СС в условиях даже простейших узкополосных помех посвящено ограниченное число работ. Одной из первых работ, где рассмотрены комбинированные воздействия на непрерывную СС (сигнал + шум + гармоническая помеха) является статья Журавлева А.Г.

Подобную ситуацию можно объяснить следующими причинами. Во-первых, представляет собой достаточно серьёзную проблему переход к марковским моделям, не существует общей методики перехода; ситуация значительно усложняется в условиях узкополосных воздействий. Во-вторых, необходимо обеспечить строгий переход от марковской модели к векторному уравнению Колмогорова - Чепмена, корректно подстроив условную плотность вероятности перехода; сложность вызвана периодическим характером нелинейности. В-третьих, усложняется задача о среднем времени до срыва слежения. Даже наличие простого сигнала без помехи, но с изменяющейся частотой, существенно усложняет решение задачи о срыве.

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

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

Задачи, решаемые в диссертации. Для достижения поставленной цели в диссертации решаются следующие основные задачи:

1. Разработка математических моделей ряда непрерывных и дискретных СС с учетом воздействия на них как шумовых, так и гармонических воздействий.

2. Разработка алгоритмов анализа статистических характеристик СС при комбинированном воздействии шумовых и гармонических помех.

3. Разработка методики решения интегрального уравнения Колмогорова-Чепмена приближёнными методами, описывающего состояние СС.

4. Исследование ДСС первого порядка приближённым методом Галёркина.

5. Анализ плотности распределения вероятности (ПРВ) фазовой ошибки в непрерывной системе при комбинированном воздействии помех приближённым методом.

6. Анализ двухконтурной СС (схемы Костаса) и получение ее статистических характеристик.

Положения, выносимые на защиту.

1. Марковские модели ряда непрерывных и дискретных СС при комбинированных воздействиях в виде смеси полезного колебания, детерминированной помехи и широкополосного гауссовского шума.

2. Методика приближенного исследования методом гармонического баланса и усреднения непрерывных СС при комбинированном воздействии.

3. Разработанные разностные схемы уравнений ФПК и уравнений Понтрягина, позволяющие численно находить ПРВ сигнала рассогласования в переходном режиме, а также среднее время и вероятность до срыва слежения.

4. Алгоритмы анализа статистических и динамических характеристик ДСС 1-го и ДСС 2-го порядков.

5. Сравнительные характеристики ДСС 1-го и ДСС 2-го порядков.

6. Результаты исследований статистических характеристик двухконтурной СС (схемы Костаса).

Научная новизна результатов.

1. Разработана разностная схема уравнения Фоккера-Планка-Колмогорова (ФПК) и уравнение Понтрягина, позволяющие получать ПРВ сигнала рассогласования в переходном режиме, а также численно находить среднее время и вероятность до срыва синхронизации.

2. Получены статистические характеристики фазового рассогласования ДСС 1 и 2-го порядка, дисперсии фазовой ошибки слежения при использовании метода Галёркина и численно-аналитических методов.

3. За счёт применения имитационной модели СС и разработанных методик, в отношении ряда систем получены уточняющие данные по сравнению с известными результатами.

4. Получены области применения приближённых методов при анализе дискретных СС 1-го порядка.

5. Приближенным методом гармонического баланса получены динамические характеристики непрерывных СС и на их основе и с использованием метода усреднения исследованы статистические характеристики этих систем.

6. Разработана модель двухконтурной схемы синхронизации в виде схемы Костаса при комбинированном воздействии и найдена дисперсия сигнала фазового рассогласования при двух моделях сигнальной формы.

Общая методика исследований. Разрабатываемые в диссертации методы анализа динамических и статистических характеристик непрерывных и дискретных СС базируются на общих положениях качественных методов теории систем автоматического управления и теории разностных уравнений, теории марковских процессов и цепей.

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

Разработанные методы и алгоритмы анализа статистических характеристик непрерывных и дискретных, в том числе и цифровых, СС ориентированы на использование персональных компьютеров.

Практическая ценность диссертации

1. В диссертации разработаны имитационная модель СС и методики исследования, позволяющие определить основные статистические характеристики различных непрерывных и дискретных СС. Разработаны алгоритмы для расчёта статистических характеристик; созданные автором программы апробированы в МГТУ им. Н.Э. Баумана, Институте криптографии, связи и информатики Академии ФСБ России, в С.-Петербургском Государственном Университете аэрокосмического приборостроения, а также в ОАО "Концерн "Созвездие".

2. На основе разработанных методик и алгоритмов автором создано программное обеспечение для анализа статистических характеристик различных непрерывных и дискретных систем. Разработанные программы позволяют оптимизировать параметры фильтра в цепи управления с целью обеспечения заданных статистических свойств непрерывных и дискретных СС при воздействии полезного сигнала и помехи.

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

Внедрение результатов работы осуществлено в:

1. НИР «Перспектива - СЧ» (ОАО «Концерн Созвездие»).

2. НИР [Синхронизация в радиосвязи и навигации, 2007 г.,

ГР 01200710182], на кафедре «Автономные информационные и управляющие системы» МГТУ им Н.Э. Баумана.

3. Учебные процессы.

3.1. Лабораторные работы на кафедре «Автономные информационные и управляющие системы» МГТУ им Н.Э. Баумана.

3.2. Учебное пособие [6].

3.3. Институт криптографии, связи и информатики (ИКСИ) Академии ФСБ.

3.4. С.-Петербургский Государственный Университет аэрокосмического приборостроения.

Апробация работы. Основные результаты работы докладывались на научно-технической конференции РНТОР и ЭС им. A.C. Попова, 2008, вып. №63, с.268-270., на семинарах кафедры СМ-5 МГТУ им. Н.Э.Баумана, а также на научных семинарах кафедры «Компьютерной математики и программирования» ГУАП (С.-Петербург, 2009 и 2010 г.).

Публикации. Основные результаты диссертации изложены в 9 работах, из них 1 отчет по НИР и 4 работы опубликовано в научных изданиях, входящих в Перечень ВАК.

Объём и структура диссертации.

Диссертация состоит из введения, пяти разделов, заключения, списка литературы (80 наименований), приложения и изложена на 163 листах машинописного текста, включая 48 листов иллюстраций.

Основное содержание работы.

Во введении дана общая характеристика диссертации: обоснована актуальность темы, сформулирована цель работы и поставлены задачи исследования, кратко изложено содержание работы, приведены основные положения, выносимые на защиту.

В первой главе выводятся основные соотношения для моделей непрерывной и дискретной СС в форме ДУ и разностных уравнений (РУ), учитывающих комбинированные воздействия в форме аддитивной смеси сигнала, шума и гармонических помех.

Структурная схема непрерывной системы синхронизации (ФАП) изображена на рис. 1.

Рис. 1

На вход ФАП поступает аддитивная смесь (рис. 1)

"«(')= "*(')+ "«(')+ «о(0. (!)

полезного сигнала (эталонного колебания) кэ(<), гармонической помехи и„(/) и полосового (белого) гауссовского шума ГБШ я0 (/).

Эталонное колебание имеет вид

м3(0 = у/2Авт(й)01 + + в0) где А - среднеквадратическое значение напряжения полезного сигнала, сод -частота сигнала, совпадающая с начальной частотой управляемого генератора (УГ), С10м в0 - начальные расстройки по частоте и фазе между эталонным колебанием и сигналом на выходе УГ.

Гармоническая помеха принимается в виде

К

"п (0 = ХУ 2 sin (йу + iV + 6»),

/=i

где N- число спектральных составляющих, попадающих в полосу пропускания линейного тракта, предшествующего ФАП; А, - среднеквадратическое значение напряжения /-ой спектральной составляющей; Q, и 6> начальные расстройки по частоте и фазе /'-ой составляющей относительно частоты и фазы на входе УГ. Полосовой гауссовский шум имеет вид

n0{t) = ^2nc{t) cos co0t + -j2ns{t)s\n(Dat, где nc (/) и ns (i) - независимые квадратурные гауссовские случайные процессы. Напряжение на выходе УГ имеет вид ujt) = yf2Aecos(a>0t + <p0(t)).

В результате преобразований получено стохастическое дифференциальное уравнение (ДУ) в символической форме

КА

рх = у- F(p)

N J

sin х + ^е, sin (х + ДС1, + Д 0,)+ — n„ (f)

•»I А

(2)

где у = П0/ КА; е. = А, / А; ДО, =П,- П0; Щ =в- в0; / = 1;Л'; Г(р) передаточная функция ФНЧ.

Низкочастотный шум пш (г) в пределах полосы пропускания ФАП является гауссовским шумом, который моделируется гауссовым белым шумом (ГБШ) с двусторонней спектральной плотностью Л'0/2 (Вт/Гц).

Помимо модели СС рис.1 в диссертации используется также двухконтурная СС (схема Костаса) (рис.2), на выход которой подается комбинированное воздействии в виде (1).

'10

£Г

л(')

-ч*>

ФНЧ

Лк, со»Ф(1)

ФНЧ

М<>

УГ KomvontA 4»№р, т

г, С)

Рис.2

Дискретная модель СС имеет вид дискретного аналога схемы Костаса.

Разностное уравнение РУ дискретной модели имеет вид

х(к + ¡) = ШыТ0 +х(к) + в/к + 1)-ес(к)-КцКф(:)АдАцА{*т(х(к)) +

+ ^ + Z 4 stoMk) + ß,kTo +e,(k)-ejk))}

A ¡A

В первой главе представлена разработанная автором имитационная БптшНпк - модель системы СС (рис.3), которая является в частности средством

визуализации процессов СС (см. рис.4). &

Ref 0se2 2*0*1 еб (rod/s)

iL

Band-Limited White Noisel

N

В

Reí Osc1

2*pMe6.2*p*1e3(r»i/s)

S

Ref Osc3 2*pi*1 e6 (rad/s)

Яе( 05с 2"р'1е6 (ГМЛ*)

IS

Power Spec) ral Oensrtyl

i î

IT

и

Power Spectral Density

1 Linear or analog Phase-Locked

^OjJ

Power Spectral Scope Voltage-Controlled Derv$i1y4 Oscillator

»"I—^ш]

"Чз

Power Spectral Density

delta_phase

To Workspace

ŒD

Power Spectral Density

Рис. 3

Управляющее напряжение УГ:

Рис. 4

Вторая глава посвящена исследованию непрерывной СС при воздействии единичной детерминированной гармонической помехи при отсутствии шума. Использован метод гармонического баланса для получения динамических характеристик СС. Найдены критические значения параметров колебания. В этом случае моделью СС служит ДУ при ¡=1; пш(1:)=0. 1 ¿¿с

-¡=— = у- /^[бш* + £5Ш(д: + АШ + Д6>)1. (4)

К-уБ (к

Предполагаемое решение ДУ (4) аппроксимируется двумя слагаемыми

•Х(/) = Л:0 + СОБ(ДП^ + М + ^) = Х0 +Х,СО5(0 4-у/), (5)

где в = в(1) = ДС2г + Ав,, а параметры колебания х0, .¡с,, у/ подлежат определению.

В предположении малой амплитуды когда выполняются условия

втх« $тлг0 + х1 со5х0соз(в + совдгя созл'0 - д:, Б1пд:0 со8(0 + у),

после подстановки (5) в ДУ (4) и последующей процедуры гармонического баланса находятся искомые параметры колебания

х,2=г2Д-'М2(ДП), (6)

где

Гг 12

АН

М (Шу^ + Ят (Р(АПу)

+ С052 (Р(Шуу\с1, х\й соэР(АП)

(7)

М( 0) 2 Л/(ДП) а также уравнение относительно фазы у/.

В соотношениях (6) и (7) М(ДО) = |/г(/й))|; />(ДС!) = а^(7й)); с1 = —.

КА

Зависимость амплитуды х1 первой гармоники предполагаемого решения (5) изображена на рис. 5 для случая системы первого порядка при у = 0.5,

е2 =0,2; 0,5; 1; 1,5 значения А = где Вь - шумовая полоса ФАП.

В1

Штриховыми линиями обозначены асимптотические значения амплитуды х,

_4£ Л '

Xi

Рис.5

Найдены критические значения отношения помеха/сигнал (ОПС) £ = £к и амплитуды хх = xxt, характеризующих срыв синхронизации.

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

В третьей главе рассматривается статистическая динамика непрерывной СС на основе метода гармонического баланса, разработанного в предыдущей главе, но с учетом всех трех составляющих комбинированного воздействия (сигнала, помехи и шума), так и с использованием других методов получения уравнения Фоккера-Планка-Колмогорова (ФПК).

Получены ПРВ сигнала фазового рассогласования в зависимости от отношения сигнал/шум (ОСШ) и помеха/сигнал (ОПС), а также от величины динамических параметров СС (частотных расстроек, параметров ФНЧ). Исследуются срывы слежения в непрерывной СС при комбинированных воздействиях.

Сначала рассматривается получение уравнения ФПК на основе аппроксимации решения стохастического ДУ (2) в форме (5), но в предположении медленного изменения первого слагаемого (5) л;0=х0(/). При этом показано, что стохастическое ДУ СС принимает вид 1 dx

——r- = y-F{p){\J0{x{) + eJx{x,) sin ]sin x0 + KA at

1 (8) +£J x(c,) COS V, COS x0 + — пш (t)}.

A

гдеУ0(л:,) и У,(x,) - функции Бесселя порядка 0 и 1 соответственно. В случае СС первого порядка найдено уравнение ФПК вида

dW 8 г -1 1 d2W

~дГ=(Х| )sin*~ 7+£j> (Xi 7&J-' (9)

где W=W(x, t) - искомая ПРВ; х = х0 (t) - марковский случайный процесс; < л\ г - ОСШ на выходе СС.

д1¥

Решение уравнения (9) в стационарном режиме (когда-=0) имеет вид

<Э/

1¥(х) = Ае"-рС(х) } е-^^и,

(10)

где А - постоянная интегрирования, V = уг\ у =у-eJ<l(xtly, р = ^й(х[);

С(и) = 1-созк; А = А(у,р') = С(у,р)1(\-е-2ю'у, С(у,р)=-

где (р) - модифицированная функция Бесселя чисто мнимого порядка, причем

, . ,2 ¿А/ту

\1,ЛР) I =

ЛУ

/02(р)+

(~1)Ч2(Р)

«2+у2

•гАя-у

XV

В результате ПРВ W(x) находится в виде

Щх) =

2л-Лт

(у сое пх - п вт их)/„ (р)

(П)

, (12)

где 1п (р) - модифицированная функция Бесселя п-го порядка.

Исследована сходимость знакопеременных рядов (11) и (12). Предложена разностная схема уравнений ФПК (9).

Результаты расчетов ПРВ \У(х) приведены на рис.6 при с!=2; е=0,8; р=1, 4, 16; у=0, а также на рис. 7 при р= 4; е =0,6; с1=-1,5 и различных значений у.

н'

1М>

О. 1

Рис.6

Рис. 7

Исследованы характеристики срыва слежения СС. Предложено использовать уравнение Понтрягина в форме

+ о,

(13)

сЬ^ ей

где 9|(х) = Л(х) + £5'т(дг + Д0); - нормированное среднее время

до срыва слежения. Получена разностная схемы уравнения Понтрягина (13).

Характеристики срыва слежения в форме зависимости нормированной величины ус среднего времени до срыва слежения от параметров СС изображены на рис. 8 при у=0,4; <1=1,5; рис. 9 у=0,4; <1=-1,5, когда

2лЧту Ус=—=— 7

Рис. 8

Рис.9

Делается вывод о различных свойствах СС в зависимости от знака частотного рассогласования при с!>0 и ё<0, а именно: при <1>0 увеличение ОПС £ приводит к росту (при любых ОСШ г) среднего времени до срыва слежения, при <1<0 имеется обратная тенденция, что подтверждает соответствующий вывод, сделанный в главе 1. Таким образом, наличие мощной помехи за пределами полосы синхронизации улучшает слежение, если сигнал и помеха имеют одинаковую по знаку расстройки от частоты собственных колебаний управляемого генератора и ухудшает слежение, если значения этих расстроек противоположные.

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

Для сопоставления с приведенным методом приводятся результаты вычисления ПРВ другими методами (помеха прицельная). Результаты расчета ПРВ решением уравнения ФПК получены на рис. 10 - 11 (рис. 10 при у=0; г=12 и на рис. 11 при ^=0,4; г=16.

-г?; # г/.' г

Рис.10 Рис.11

В четвертой главе рассматривается статистическая динамика двухконтурной схемы Костаса; приводится вывод стохастического ДУ - модели схемы Костаса. Получены соотношения для дисперсии фазового рассогласования при комбинированном воздействии.

Получены нормированные значения дисперсии сигнала ошибки при двух сигнальных функциях: а) двоичный сигнал без возврата к нулю (N117), и б) при манчестерском коде.

В пятой главе разрабатывается приближенный метод анализа статистических характеристик дискретной СС (ДСС) на основе марковской модели шума (уравнение Колмогорова-Чепмена) и приближенного метода Галёркина. Получены ПРВ фазового рассогласования как при отсутствии помехи (Приложение 1), так и при наличии комбинированного воздействия (аддитивная смесь сигнала, помехи и шума). В случае ДСС первого порядка стохастическое разностное уравнение (РУ) ДСС первого порядка является частным случаем РУ (3).

Показано, что переходная ПРВ имеет вид [5]

Ч(*\г) = -4—ехр{~[х- : + Г /2а2} л/ ¿по

где к1(:) = (1+А1)ят(2)-р, Кх - амплитуда гармонической помехи.

При использовании метода Галёркина в качестве системы {ц/„(х)} ортонормированных функций используется система тригонометрических функций [5,7]

\jijx )}= {/; ¡т( х);со.ч( х);$т( 2х ),соз( 2х )....}. После подстановки предполагаемого решения в уравнение Колмогорова-Чепмена получены в виде ряда ПРВ \?^(х) сигнала рассогласования

/7=0

Результаты представлены в виде графиков ПРВ \У(х) рис. 13 (при числе слагаемых N=30; А|=0,8; г=уаг; Т0=0,5; у= 0; <1=0).

-2 4^5 1=2 пБ

Рис. 12 Рис. 13

Также проводится анализ статистических характеристик ДСС второго порядка. На основе марковской модели шума получены ПРВ сигнала рассогласования и исследован процесс срыва синхронизации.

Стохастическое (РУ) в данном случае является частным случаем РУ (3). Показано, что система двух РУ, характеризующих ДСС 2-го порядка приводит к зависимости нового состояния систем не только от предыдущего состояния, но

и от более раннего, то есть исходная система ДУ не описывает марковскую последовательность.

Осуществлено преобразование координат, в результате которого получена система РУ, описывающая марковскую последовательность, двумерная ПРВ которой подчиняется векторному уравнению Колмогорова-Чепмена.

Численным методом найдено решение этого уравнения [8]. Результаты вычисления ПРВ фазовой ошибки представлены в виде графиков. Делается вывод о согласованности решений уравнения ФПК для непрерывной СС и для ДСС первого и второго порядков.

\У(х)

Рис. 14

В Приложении 1 [5] представлены результаты анализа ДСС 1-го порядка, полученные методом Галёркина при отсутствии гармонической помехи.

Основные результаты и выводы

1. Создана имитационная модель СС при комбинированном воздействии в виде аддитивной смеси сигнала, помехи и гауссовского шума.

2. Разработаны разностные схемы уравнений ФПК и Понтрягина, позволяющие вычислять ПРВ сигнала рассогласования в переходном режиме и находить характеристики срыва слежения при любой нелинейной характеристики фазового дискриминатора.

3. На основе приближенного метода Галёркина исследована ДСС как при отсутствии гармонической помехи, так и при наличии комбинированного воздействия на ДСС.

4. Разработаны алгоритмы и программы, позволяющие анализировать динамические и статистические характеристики непрерывных СС.

5. Показано, как в детерминированном, так и в статистическом случаях, что воздействие гармонической помехи по-разному влияет на динамику СС: подверженность СС помехе увеличивается, если частоты помехи и сигнала находятся по разные стороны от частоты покоя управляемого генератора.

6. Получены статистические характеристики (ПРВ сигнала рассогласования, среднее время до срыва слежения, среднее значение частотного рассогласования) систем синхронизации как непрерывных, так и дискретных. Рассмотрены две разновидности помех: прицельная и отстроенная по частоте.

7. Исследована двухконтурная СС (схема Костаса). Найдена дисперсия ошибки при двух формах сигнальной функции.

Основные результаты диссертации представлены в следующих работах:

1. Шахтарин Б.И., Иванов А.А, Рязанова М.А. Частотная и фазовая синхронизация с коррекцией импульсной характеристики канала передачи в OFDM - системе // Научный вестник МГТУ ГА. Радиофизика и радиотехника. -2008. -№126. - С.107-117.

2. Иванов A.A., Быков A.A., Рязанова М.А.. Статистический анализ цифровых систем синхронизации // Успехи современной радиотехники. - 2008. -№2.-С. 68-76.

3. Рязанова М.А., Митин C.B.. Схема Костаса в условиях наличия активных помех и аддитивного шума // Научная сессия, посвященная Дню радио. - 2008. — Вып. № 63. - С.457-459.

4. Рязанова М.А. Исследование схемы Костаса в условиях наличия непрерывного источника помех и аддитивного шума // Сборник научных работ лауреатов конкурса им. Е.Р. Дашковой. - 2008. - Вып. №2 . - С. 50-60.

5. Рязанова М.А. Сравнительный анализ приближенных методов для получения статистических характеристик дискретных систем синхронизации первого порядка // Сборник научных работ лауреатов конкурса им. Е.Р. Дашковой. -2008. - Вып. №2 . - С. 61-71.

6. Синхронизация в радиосвязи и радионавигации / Шахтарин Б.И., Иванов A.A., Рязанова М.А. [и др.] - М.:Гелиос АРВ, 2007 - 256 с.

7. Иванов A.A., Рязанова М.А., Кровяков И.И. Анализ бесфильтровой дискретной системы фазовой автоподстройки при наличии нормального белого шума // Научный вестник МГТУ ГА. Радиофизика и радиотехника. - 2007. -№117.-С. 137-148.

8. Рязанова М.А., Иванов A.A., Быков A.A.. Статистический анализ дискретной системы синхронизации 2-го порядка в условиях комбинированных воздействий // Научный вестник МГТУ ГА. Сер. Радиофизика и радиотехника. - 2007. - №117. -С. 160-168.

9. Синхронизация в радиосвязи и радионавигации: Отчет по НИР по теме СМ-2-285/МГТУ НИИ СМ; Рук. Темы Шахтарин Б.И. Исп.: Рязанова М.А. [и др.]: ГР№ 01.200710182,- М.:2007,- 285 с.

Подписано в печать: 09.04. И

Объем: 1,5 усл.п л. Тираж: 100 экч. Заказ № 785 Отпечатано в типографии «Реглет» 1 19526, г Москва, пр-т Вернадского,39 (495) 363-78-90; www.reglet.ru

Оглавление автор диссертации — кандидата технических наук Рязанова, Мария Алексеевна

Введение.

Глава 1. Анализ моделей систем синхронизации (СС).

1.1. Модель непрерывной системы синхронизации (СС).

1.2. Математические модели дискретных СС при наличии помех.

1.3. Выводы по первой главе.

Глава 2. Непрерывная система синхронизации при наличии гармонической помехи и отсутствии аддитивных шумов.

2.1. Введение.

2.2. Модель СС.

2.3. Обзор некоторых предыдущих работ.

2.4. Анализ случая 1: ^^ ^.

2.5. Выводы по второй главе.

Глава 3. Работа непрерывной СС при наличии помехи и аддитивного гауссовского белого шума (ГБШ).

3.1. Введение.

3.2. Анализ случая

3.3. Анализ характеристик СС - случай

3.4. Выводы по третьей главе.

Глава 4. Исследование схемы Костаса при наличии непрерывной помехи и аддитивного шума.

4.1. Введение.

4.2. Модель схемы Костаса при наличии помехи на входе.

4.3. Статистические характеристики возмущений.

4.4. Характеристики фазовой ошибки.

4.5. Частные случаи сигнальной формы.

4.6. Выводы по четвертой главе.

Глава 5. Статистические характеристики дискретной системы синхронизации 1-го порядка и 2-го порядка при наличии гармонической помехи и шума на входе.

5.1. Метод Галёркина.

5.2. Метод замены интегрального ядра.

5.3. Вид уравнения Колмогорова — Чепмена для случая отсутствия гармонической помехи.

5.4. Вид уравнения Колмогорова - Чепмена для случая наличия гармонической помехи.

5.5. Сравнение результатов, полученных различными методами.

5.6. Выводы по пятой главе.

Введение 2011 год, диссертация по информатике, вычислительной технике и управлению, Рязанова, Мария Алексеевна

Актуальность темы

Современные системы и устройства управления, навигации и радиолокации, а также информационно-измерительные комплексы широко применяют системы синхронизации (СС). Эти системы выполняют многочисленные задачи; к числу наиболее распространенных из них относятся: синхронизация и демодуляция двоичных символов цифровой информации, измерение частоты и фазы сигналов, тактовая синхронизация, синтез сложных радиотехнических сигналов, синтез сетки высокостабильных частот, стабилизация частот генераторов различных диапазонов, когерентная демодуляция аналоговых и цифровых сигналов с частотной и фазовой модуляцией.

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

Особое значение при проектировании и создании СС имеют проблемы оценки устойчивости этих систем при воздействии помех различного рода. При этом статистические моменты фазовой и частотной ошибок слежения (среднее и среднеквадратическое значения) не дают полной информации о поведении СС. Поскольку СС — существенно нелинейная система, то в ряде случаев необходимо знание плотностей распределения вероятностей (ПРВ) переменных состояния. Особенностью СС с рядом других систем (не фазовых) является существование множества устойчивых состояний равновесия, а в отдельных предельных случаях и устойчивых периодических движений 1-го и 2-го рода, что еще более усложняет картину при действии шумов.

Исследования СС существенно усложняются, если наряду с шумовой (флуктуационной) помехой на СС воздействует узкополосная негауссовская помеха. Здесь особенно важным является исследование СС при воздействии гармонических помех. К таким помехам относятся помехи соседних каналов в системах радиосвязи, помехи многолучевости, помехи космическим навигационным системам (GPS и ГЛОНАСС), организованные помехи радиолокационным станциям.

Основным методом исследования СС при наличии гауссовского шума на входе в настоящее время является метод марковских случайных процессов. Пионерами использования этих методов применительно к СС являются Р.Л.Стратонович и В.И.Тихонов. Существенный вклад в развитие теории синхронизации при наличии гауссовского шума на входе сделали Б.И. Шахтарин, В. Линдсей, А. Витерби, Дж. Холмс, В.Д. Шалфеев, H.H. Удалов, В.Н. Белых, В.Н. Кулешов, В.Д. Разевиг, В.В. Шахгильдян, А. Вайнберг и другие. В результате к настоящему времени теория непрерывных СС при наличии шумовых воздействий в основном завершена. В то же время теория дискретных систем синхронизации (ДСС) даже при шумовых воздействиях разработана весьма незначительно. И если в непрерывных СС достаточно подробно методом марковских процессов исследованы системы 1, 2-го порядков, применительно к ДСС - это в основном касается систем только 1 -го порядка. Такие исследования выполнены М.И. Жодзишским, В.Н. Кулешовым, В.В. Шахгильдяном, Б.И. Шахтариным, В.Н. Белых, В.П. Сизовым, Дж. Холмсом, X. Осборн, С. Гуптой.

Исследованиям дискретных СС в условиях даже простейших узкополосных помех посвящено ограниченное число работ. Одной из первых работ, где рассмотрены комбинированные воздействия на непрерывную СС (сигнал + шум + гармоническая помеха) является статья Журавлева А.Г.

Подобную ситуацию можно объяснить следующими причинами. Во-первых, представляет собой достаточно серьёзную проблему переход к марковским моделям, не существует общей методики перехода; ситуация значительно усложняется в условиях узкополосных воздействий. Во-вторых, необходимо обеспечить строгий переход от марковской модели к векторному уравнению Колмогорова - Чепмена, корректно подстроив условную плотность вероятности перехода; сложность вызвана периодическим характером нелинейности. В-третьих, усложняется задача о среднем времени до срыва слежения. Даже наличие простого сигнала без помехи, но с изменяющейся частотой, существенно усложняет решение задачи о срыве.

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

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

Задачи, решаемые в диссертации. Для достижения поставленной цели в диссертации решаются следующие основные задачи:

1. Разработка математических моделей ряда непрерывных и дискретных СС с учетом воздействия на них как шумовых, так и детерминированных гармонических воздействий.

2. Разработка алгоритмов анализа статистических характеристик СС при комбинированном воздействии шумовых и организованных помех.

3. Разработка методики решения интегрального уравнения Колмогорова-Чепмена приближёнными методами, описывающего состояние СС приближенными методами.

4. Исследование ДСС первого порядка приближённым методом Галёркина.

5. Анализ плотности распределения вероятности (ПРВ) фазовой ошибки в непрерывной системе при комбинированном воздействии помех приближённым методом.

6. Анализ двухконтурной СС (схемы Костаса) и получение ее статистических характеристик.

Общая методика исследований. Разрабатываемые в диссертации методы анализа динамических и статистических характеристик непрерывных и дискретных СС базируются на общих положениях качественных методов теории систем автоматического управления и теории разностных уравнений, теории марковских процессов и цепей.

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

Разработанные методы и алгоритмы анализа статистических характеристик непрерывных и дискретных, в том числе и цифровых, СС ориентированы на использование персональных компьютеров.

Научная новизна результатов.

1. Разработана разностная схема уравнения Фоккера-Планка-Колмогорова (ФПК) и уравнение Понтрягина, позволяющие получать ПРВ сигнала рассогласования в переходном режиме, а также численно находить среднее время и вероятность до срыва синхронизации.

2. Получены статистические характеристики фазового рассогласования ДСС 1 и 2-го порядка, дисперсии фазовой ошибки слежения при использовании метода Галёркина и численно-аналитических методов.

3. За счёт применения имитационной модели СС и разработанных методик, в отношении ряда систем получены уточняющие данные по сравнению с известными результатами.

4. Получены области применения приближённых методов при анализе дискретных СС 1-го порядка. Изучено влияние размерности уравнения на результаты вычислений.

5. Приближенным методом гармонического баланса получены динамические характеристики непрерывных СС и на его основе и с использованием метода усреднения исследованы статистические характеристики этих систем.

-86. Разработана модель схемы Костаса при комбинированном воздействии и найдена дисперсия сигнала фазового рассогласования при двух моделях сигнальной формы.

Практическая ценность диссертации

1. В диссертации разработаны имитационная модель СС и методики исследования, позволяющие определить основные статистические характеристики различных непрерывных и дискретных СС. Разработаны алгоритмы для расчёта статистических характеристик; созданные автором программы апробированы в МГТУ им. Н.Э. Баумана, Институте криптографии, связи и информатики Академии ФСБ России, в С.-Петербургском Государственном Университете космического приборостроения, а также в ОАО "Концерн "Созвездие".

2. На основе разработанных методик и алгоритмов автором создано программное обеспечение для анализа статистических характеристик различных непрерывных и дискретных систем. Разработанные программы позволяют оптимизировать параметры фильтра в цепи управления с целью обеспечения заданных статистических свойств непрерывных и дискретных СС при воздействии полезного сигнала и помехи.

3. Предложенные и развитые в диссертации методики и разрабатываемые на их основе алгоритмы и программы могут быть использованы в научно-исследовательских и опытно-конструкторских работах для анализа статистических свойств непрерывных дискретных СС и синтеза дискретных СС различного назначения.

Положения, выносимые на защиту.

1. Марковские модели ряда непрерывных и дискретных моделей СС при комбинированных воздействиях в виде смеси полезного колебания, детерминированной помехи и широкополосного гауссовского шума.

2. Методика приближенного исследования методом гармонического баланса и усреднения непрерывных СС при комбинированном воздействии.

3. Разработанные схемы уравнений ФПК и уравнений Понтрягина, позволяющие численно находить ПРВ сигнала рассогласования в переходном режиме, а также среднее время и вероятность до срыва слежения.

4. Алгоритмы анализа статистических и динамических характеристик ДСС 1 -го и ДСС 2-го порядков.

5. Сравнительные характеристики ДСС 1-го и ДСС 2-го порядков.

6. Результаты исследований статистических характеристик двухконтурной СС (схемы Костаса).

Объём и структура диссертации.

Диссертация состоит из введения, пяти разделов, заключения, списка

Заключение диссертация на тему "Разработка алгоритмов анализа систем синхронизации при воздействии гармонических и шумовых помех"

Выводы и заключение

1. Создана имитационная модель СС при комбинированном воздействии в виде аддитивной смеси сигнала, помехи и гауссовского шума.

2. Разработаны разностные схемы уравнений ФПК и Понтрягина, позволяющие вычислять ПРВ сигнала рассогласования в переходном режиме и находить характеристики срыва слежения при любой нелинейной характеристики фазового дискриминатора.

3. На основе приближенного метода Галёркина исследована ДСС как при отсутствии гармонической помехи, так и при наличии комбинированного воздействия на ДСС.

4. Разработаны алгоритмы и программы, позволяющие анализировать динамические и статистические характеристики непрерывных СС.

5. Показано, как в детерминированном, так и в статистическом случаях, что воздействие гармонической помехи по-разному влияет на динамику СС: подверженность СС помехе увеличивается, если частоты помехи и сигнала находятся по разные стороны от частоты покоя управляемого генератора.

6. Получены статистические характеристики (ПРВ сигнала рассогласования, среднее время до срыва слежения, среднее значение частотного рассогласования) систем синхронизации как непрерывных, так и дискретных. Рассмотрены две разновидности помех: прицельная и отстроенная по частоте.

7. Исследована двухконтурная СС (схема Костаса). Найдена дисперсия ошибки при двух формах сигнальной функции.

Библиография Рязанова, Мария Алексеевна, диссертация по теме Системный анализ, управление и обработка информации (по отраслям)

1. Шахтарин Б.И. Статистическая динамика систем синхронизации.- М.: Радио и связь, 1998. 488 с.

2. Шахтарин Б.И. Анализ систем синхронизации при наличии помех. -М.: ИПРЖР, 1996. 252 с.

3. Тузов Г.И. Выделение и обработка информации в доплеровских системах. — М.: Советское радио, 1967. — 256 с.

4. Рыжков А.В., Попов В.Н. Синтезаторы частот в технике радиосвязи.- М.: Радио и связь, 1991. 264 с.

5. Журавлев В.И. Поиск и синхронизация в широкополосных системах.- М.: Радио и связь, 1986. 240 с.

6. Линдсей В. Системы синхронизации в связи и управлении: Пер. с англ. / Под ред. Ю.Н. Бакаева и М.В. Капранова. М.: Сов. Радио, 1978. - 600 с.

7. Витерби Э.Д. Принципы когерентной связи. — М.: Советское радио. 1970.-350 с.

8. Best R. Е. Phase-locked loops: design, simulation, and applications. N.Y.r Fifth ed., 2003.-421 p.

9. Crawford J.A. Frequency Synthesizer Design Handbook. — Norwood: Artech. House. Inc., 1994.-435 p.

10. Crawford J.A. Advance phase-lock techniques L.: Artech House., 2008. -510pp.

11. Manassewitsch V. Frequency Synthesizers. Theory and Design. New York: Third Edition., 1987. - 611 p.

12. Bahai A.R.S., Saltzberg B.R., Ergen M. Multi carrier digital communications. Theory and applications of OFDM. -N.Y.: Springer, 2004. - 41 lp.

13. Meyer H. Moeneclaey M. Fettle digital communication receivers. Synchronization, channel estimation end signal processing. N.Y.: Willey, 1998.- 827 pp.

14. Беллами Д.К. Цифровая телефония, 3-е изд. М.:Экотрендз, 2004.- 639с.- 12415. Скляр Б. Цифровая связь. Теоретические основы и практическое применение. 2-е изд. М.: Вильяме, 2003, - 1104с.

15. Прокис Д. Цифровая связь. М.: Радио и связь, 2000. — 800с.

16. Феер К. Беспроводная цифровая связь: Пер. с англ. / Под ред. В.И. Журавлева. М.: Радио и связь. - 2000. - 520 с.

17. Г.И. Тузов и др. Помехозащищенность радиосистем со сложными сигналами, Под ред. Тузова Г.И. — М.: Связь, 1985. -279 с.

18. Тихонов В.И. Влияние шумов на работу схемы фазовой автоподстройки частоты // Автоматика и телемеханика. — 1959. — №9. — С. 11881196.

19. Тихонов В.И. Работа фазовой автоподстройки частоты при наличии шумов // Автоматика и телемеханика. 1960. - №3. - С.301-309.

20. Витерби А. Исследование динамики систем фазовой автоподстройки частоты в присутствии шумов с помощью уравнения Фоккера-Планка // ТИИЭР. 1963. - Т. 51, — №12. - С. 1704-1722.

21. Обрезков Г.В., Разевиг В.Г. Методы анализа срыва слежения. М: Советское радио, 1972. - 502 с.

22. Тихонов В.И., Харисов В.Н. Статистический анализ и синтез радиотехнических устройств и систем: Учебное пособие. — М.: Радио и связь. 1991.-608 с.

23. Bruno F. Tracking performance and loss of lock of carrier loop due to presence of a spoofed spread of spectrum signal // Proc., 1973, Symp. Spread Spectrum Communications. San Diego. (CA.), 1973. - v.l. - P. 71-75.

24. Blanchard A. Interference in phase locked loops // IEEE Trans. On aerospace and electronic systems. - 1974. - V. AE S-10, N 5, - P. 686-697.

25. Sue M.K. Effects of CW interference on the carrier tracking loop of the deep space network // IEEE Trans on electromagnetic compatibility. -1984. -V. EMC 26, N 3, -P. 120-128.

26. Endo T. Suzuki T. Influence of interfering second target on a certain amplitude comparison monopulse radar system // Proceedings of IEEE Global

27. Telecommunications conference Atlanta, GA. USA 26-29 November 1984. -P. 224229.

28. Kantak A. Davarian F. Performance of PLL in the presence of a CW interferer // Proceedings of IEEE Global Telecommunications conference Atlanta, GA. USA 26-29 November 1984. -P. 230-236.

29. Levitt B.K. Carrier tracking loop performance in the presence of a strong CW interference // IEEE Trans. On communications, 1981. V. COM -29, N6, - P. 911-916.

30. Olenberger C.F. Effects of automatic gain control on phase locked loop behavior in the presence of interference // IEEE Trans. On aerospace and electronic systems, 1976.-V. AES-12, N 11, - P. 803-805.

31. Kliger I.E. Olenberger C.F. Phase lock loop jump phenomenon in the presence of two signals // IEEE Trans. On aerospace and electronic systems, 1976. -V. AES-12, N 1,-P. 55-63.

32. Kliger I.E. Olenberger C.F. Multiple target effects on monopulse signal processing // IEEE Trans. On aerospace and electronic systems, 1975. -V. AES -11, N5,-P. 795-804.

33. Ziemer R.E. Perturbation analysis of the effects of CW interference in Costas loops // Rec. Nat. Telecommun. Conf. Houston. TX, Dec. 4-6, 1972. P. 20G-1-20G-6.

34. Britt C.L. Palmer D. Effects of CW interference on narrowband second -order phase locked loop // IEEE Trans. 1967, v. AES -3, N1, P. 123-135.

35. Yoon C.Y., Lindsey W.C. Phase locked loop performance in the presence of CW interference and additive noise // IEEE Trans. On communications 1982. -V. COM -30, N10. - P. 2305-2311.

36. Hasan P. Cycle-slip statistics of a first order phase locked loop in the presence of additive noise and co channel interference // European trans, on Telecommunications (ETT), 1994. -V.5, N 5, - P. 619-623.

37. Hasan P. Multiple co channel interference effects in a first order phase -locked loop // ETT, 1994. V. 5, N3, - P. 319-326.- 12638. Karsi M.F., Lindsey W.C. Effects of CW interference on carrier tracking // Proc. IEEE Conf. 1994, P. 301-305.

38. Sarcar B.C. Phase error dynamics of a first order phase-locked loop in the presence of co-channel tone interference and additive noise // IEEE Trans., 1900. — V. COM-38, № 7, P. 962-965.

39. Башмаков M.B. Статистические характеристики дискретных СФС в "условиях комбинированных воздействий. Кандидатская диссертация — М.: 2001,-216с.

40. Simon М.К., Lindsey W.C. Optimum performance of Costas receivers containing soft bandpass limiters. // IEEE Trans., 1977. V. COM-25, № 8, - P. 822831.

41. Lindsey W.C., Simon M.K. Optimum performance of suppressed carrier receivers with Costas loop tracking. // IEEE Trans., 1977, v. COM-25, № 2, P. 215227.

42. Шахтарин Б.И. Случайные процессы в радиотехнике, т.1. Линейные преобразования. 4-е изд. перераб. и доп. М.: Горячая линия - Телеком. 2010, 520 с.

43. Karsi M.F., Lindsey W.C. Effects of CW interference on phase-locked performance // IEEE Trans., 2000, v. COM-48, № 5, p. 886-896.

44. Прокис Дж. Цифровая связь. М.: Радио и связь. 2000, - 600 с.

45. Sarcar B.C. On the joint statistics of amplitude and phase of a signal with co-channel interference.// Proc. IEEE, 1988, v. 76, №3, P. 298-299.

46. Шахтарин Б.И., Иванов A.A., Рязанова M.A. Частотная и фазовая синхронизация с коррекцией импульсной характеристики канала передачи в OFDM системе // Научный вестник МГТУ ГА. Радиофизика и радиотехника. -2008. - №126. - С.107-117.

47. Иванов А.А., Быков А.А., Рязанова М.А. Статистический анализ цифровых систем синхронизации // Успехи современной радиотехники. 2008. — №2. С. 68-76.

48. Рязанова М.А., Митин С.В. Схема Костаса в условиях наличияактивных помех и аддитивного шума // Научная сессия, посвященная Дню радио. 2008. - Вып. № 63. - С.457-459.

49. Рязанова М.А. Исследование схемы Костаса в условиях наличия непрерывного источника помех и аддитивного шума // Сборник научных работ лауреатов конкурса им. Е.Р. Дашковой.- 2008.- Вып. №2 .- С. 50-60.

50. Mannella R. Computer experiment in nonlinear stochastic physics // Noise in nonlinear dynamical systems, 1983, v.3, ed. F. Moss, P.V.E. McClintock.

51. Mannella R., Palleschi V. Fast and precise algorithm for computer simulation of stochastic differential algorithms // Phisical Rewiev A., 1989, v.40, №6, pp. 2281-2286.

52. Сизых В.В. Статистическая динамика систем синхронизации. Докт. диссертация,-М.: 1998.

53. Башмаков М.В., Казаков JI.H., Шахтарин Б.И. Анализ срыва слежения дискретных систем фазовой синхронизации в условиях подвижных поглощающих границ // Радиотехника и электроника, 2002, №12, с. 1435-1449.

54. Башмаков М.В., Казаков JI.H., Шахтарин Б.И. Статистические характеристики дискретных СФС в условиях комбинированных воздействий. // Электротехника и электроника, 2002, №10, с. 1187-1197.

55. Шахтарин Б.И. Спектральный анализ фазовой автоподстройки частоты при наличии шума // Автоматика и телемеханика, 1997, №2, с. 118-129.

56. Журавлев А.Г. Работа системы фазовой автоподстройки частоты при гармонических помехах // Радиотехника, 1963, т. 18, №9, с. 38-46.

57. Ватсон Г.Н. Теория бесселевых функций. М.: Иностр. Лит. 1949.

58. Справочник по специальным функциям / Под ред. М. Абрамовича и Н. Стиган. М.: Наука, 1979, - 830с.

59. Статистические характеристики ДСС 1-го порядка при отсутствии помехина входе

60. Плотность распределения вероятностей (ПРВ) сигнала рассогласования ДСС 1-го порядка может быть получена методом Галёркина. В общем случае выражение для ПРВ будет иметь вид

61. WN(x)~ J qx(х|z)WN(z)dz = О, (П.5.2)тсгде х е (—я,я);дх(х \ z) = + 2яп | z) переходная ПРВ, приведённая кинтервалу (—я, я);д(х I z) = (1 / л/2^с7)ехр{-х- z + T0h(z)f / 2сг2. (П.5.3)

62. Для произвольной функции /(х) на этом основании можно записать соотношение

63. КхШАх) 2)1Г„Ю(Ь1фс = О (П.5.4)-л -л

64. По условию система = 0,1,2.) полная, поэтому справедливо00следующие разложение /(х) = ^1/тУ/т(х)-> подставив этот ряд в (П.5.4),т=0получим1. Е/,. = (П.5.5)т=О —я:1. Отсюда следуетл лт = 0,1. (П.5.6)-л

65. Выберем коэффициенты cn(N) так, что были равны нулю первые из интегралов (П.5.3). Такой подход соответствует нулевой проекции невязкилeN{x) = WN(x)-¡д{(х\ z)WN(z)dz-лна подпространство (N+1) функций у/т(х) . Запишем (П.5.6) в видел л л

66. J WN(z) J gi(x\z)y/m(x)dxdz = J WN(x)y/m(x)dx; m = 0,1. (П.5.7)л —л —л

67. С учётом ортогональности функций у/т(х) в правой части получаем1. Л со Л

68. J WN(х)¥т(x)dx = (AO jV„(х)у/т (x)dx = утст (АО, (П.5.8)-л -ллгде ym=\y/l(x)dx.л1. Обозначим интеграллчх(х\ = 40) = (q,vm\ (П.5.9)-лгде (q,yfm) скалярное произведение функций qx(х | z) и у/т(х). В результате получаем сумму1. П. 5.10)

69. Таким образом, коэффициенты cm(N)(m = l,N) должны определяться из решения системы линейных уравненийi^amnlym)cn{N) = cm(Ny, m = Qjj (П.5.11)п=О

70. Возьмём в качестве системы ортогональных функций {ут(х)} системутригонометрических функций {у/т(х)} = {l;sin(x);cos(x);sin(2x);cos(2x).} вэтом случаеcos(mx) при i чётномsin(mx) при i нечётном1. П.5.12)2 п при т = 071 при тф 0

71. Используя преобразования получаем следующие выражения для вычисления a{j, i,j номера столбцов и строк в матрице.

72. Введём следующие обозначения:ill при i чётном7 + 2)/2 при i нечётном1. Ут =т = <1. П.5.13)п = {j / 2 при j чётному + 2) / 2 при у нечётном

73. А = е-т2°2'2 cos (тТ0/3) В = e^*2'1 sin (тТф), (п.5.14)где ^ = (П.5.15)

74. Зп(х) функция Бесселя 1-го рода я-го порядка. Выражение (П.5.11) представим в матричном виде 1-А.См=В„1. П.5.16)1. П.5.17)- единичная матрица размером Л^х N, А матрица с элементами а у / л1. Л =1. Г« п а21 \ «0.л л л12 а22 аа1. Л Л л1 п Л п )1. П.5.18)

75. СМ,ВМ вектор столбец коэффициентов,

76. Сы =(С\(Ю с2{М) с„(Ю)Т (П.5.19)1. Вм =1. С а10 а2р ссмо ^ ■21. П.5.20)2 л1 2 л1 2л'

77. Для оценки быстроты сходимости ряда П.1 рассмотрим частный случай при Р = 0. Вычислим коэффициенты СЛ, при следующих параметрах Т0 =1, n = 15,р = уаг. В таблице 5.5 приведены вычисленные коэффициенты. Таблица 5.5р =7с1Ь р = 10 6Ь р=14 аь

78. С/ 0,296199894 0,302726245 0,307120213с2 0,238670644 0,26040611 0,275855716

79. С3 0,166540619 0,202608262 0,230660365с4 0,10064443 0,142585031 0,17954789с5 0,052681382 0,090763105 0,130109191с6 0,023887724 0,052260379 0,087771632с7 0,009384126 0,027219107 0,055122023

80. С8 0,003194198 0,012824003 0,032225998

81. С9 0,000942153 0,005465551 0,017540624

82. Сю 0,000240829 0,002107241 0,008887571

83. Си 5,33527Е-05 0,00073498 0,004192273

84. С12 1.02445Е-05 0,000231915 0,001841025с13 1,70503Е-06 6,62036Е-05 0,000752683

85. С14 2,45981Е-07 1,7098Е-05 0,000286486

86. Си 3,07618Е-08 3,99508Е-06 0,000101515

87. На рисунках 5.25 5.30 показаны графики ПРВ Щх) рассчитанные при различных параметрах

88. Рис. 5.25. ПРВ фазовой ошибки ДСС 1-го порядка при N=30, Тв=0.01 с, Р=0,p^r=var дБ.0.41. W(x)035 0.3 0.25 0.2 0.15 0.1 0.05 0

89. Рис. 5.26. ПРВ фазовой ошибки ДСС 1-го порядка при N=30, Т0=0.01 с,b=/3=var, г=1 дБ.1. W(x) °'5 0.450.40.350.30.250.20.150.10.050

90. Рис. 5.27. ПРВ фазовой ошибки ДСС 1-го порядка при N=30, T0=var с, р=0,г=1 дБ.08 W(x) 0.70.6