Воздействие гармонической помехи на фазовую автоподстройку второго порядка Текст научной статьи по специальности «Физика»

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. II. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408

электронный научно-техническим журнал

Воздействие гармонической помехи на фазовую

автоподстройку второго порядка

# 05, май 2013

Б01: 10.7463/0513.0551477

Асланов Т. Г.

УДК 621.396

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

Все эти системы работают в условиях воздействия помех [1-3].

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

Помехоустойчивости систем синхронизации посвящен ряд работ [4, 5, 7 и др.].

В данной статье путем анализа фазовых портретов ФАП получены уравнения захвата за сигнал и за гармоническую помеху.

1. Влияние гармонической помехи на систему ФАП второго порядка

Функциональная схема ФАП второго порядка приведена на рис. 1,а, где ив - смесь сигнала и гармонической помехи, ид - напряжение на выходе фазового детектора (ФД), ир - напряжение на выходе фильтра, иг - сигнал управляемого генератора (УГ). Данная схема отличается от схемы первого порядка наличием пропорционально-интегрирующего фильтра (ПИФ) на выходе фазового детектора приведенного на рис. 1,б.

На вход системы воздействует аддитивная смесь сигнала и гармонической помехи. Дифференциальное уравнение ФАП в символической форме имеет вид [6]

Россия, МГТУ им. Н.Э. Баумана tabasik@gmail.com

Введение

(1)

где в = Лп/Ас - отношение помеха/сигнал; Ас, Ап - соответственно амплитуды колебаний сигнала и помехи; - время; О - полоса синхронизации ФАП; Од = ®с — ®0 -расстройка частоты сигнала ас и частоты УГ а0; АО = (Оп — 0)с - разность частот помехи и сигнала а = - оператор дифференцирования, К (а) - передаточная

функция ПИФ.

Дифференциальное уравнение (1) может быть записано в форме

ах = Од — ОК (а){вт х + вБт [ х — (Од — Оп ) ^ ]},

(2)

где шп = ёфп/дХу; фп - фаза колебания помехи детектируемая ФД в момент времени Введем обозначения

у = х-(Од — Оп)Iу; Фу (х,у) = Бтх + вБту.

(3)

а) б)

Рис. 1 Функциональная схема ФАП второго порядка

Преобразуем передаточную функцию ПИФ к виду

К (а) = (1 + атфа)/(с0 + аТфа) = а + ^Д сд + аТфа),

(4)

—2

где ^ = 1 + атфа; ад =Отф, Тф - постоянная времени ФНЧ; сд = 1V 0, сд = 0

соответствует вырожденному ПИФ.

Тогда второе слагаемое в (1) можно записать в виде

К (ст)Ф1 ( X, у ) = аф>1 ( х, у)-г; г = -

£

С0 + тф°

ф1 ( х, у ).

(5)

При условии (3) уравнению (2) можно сопоставить систему из трех дифференциальных уравнений первого порядка

= Г (w),

(6)

где w = (х, у, г )Т; Г - некоторая вектор-функция переменных х, у, г. Третье уравнение системы сразу получаем из (5)

шг=-ао [сог+^ф1 (х, у)],

где а0 = у Отф.

По уравнению (3) находим

!=I-(О-Оп).

(7)

где Р = Оо/О; вп = Од/О.

С учетом (3) по уравнению (2) получаем

— = О0-О

(

а +--

V 1 + тф°у

ф1 ( ^ у).

Отсюда

Шх = в-аФ1 ( х, у )--^-Ф1 ( х, у ) = в-аФ1 ( х, у )-г.

ш Со + ТфО

Обозначим [7]

в—* = А (г),

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

ёх = А (г)—аФ1 (х у)

и с учетом (7) второе дифференциальное уравнение

ёУ = Вд (г) —аф (х, у); Вд (г) = рп — г.

Итак, система дифференциальных уравнений (6) окончательно принимает вид [7]

ё: = А (*)—аф1 (х у);

ёу = В0 ( г ) — аф1 ( х, у );

йг

(8)

йг

—ао [с0 г + ^ф1 (х, у )],

Рассмотрим частный случай ФАП первого порядка. Тогда при ад ^ ^ по третьему уравнению системы получим

г = (1 — а )ФХ ( х, у ).

В результате по системе из трех дифференциальных уравнений (8) находим систему из двух дифференциальных уравнений

I=в—ф- (х, у);

йу = Рп — аф1 (х, у).

/0 0 Введем параметр Х = а/ ао, и обозначим а^ао. Тогда по (8) получим систему

дифференциальных уравнений

— = Р + г -аХ( sin х + в sin у); Шг

— = в + Лр + г -аХ ( sin х + в sin у); Шг

Шг

Шг

-а0 [сог + (1 - асо )(sin х + в sin у)],

(9)

где Лр = рП-р.

Процедуру усреднения [8] можно применять в случае, если а «1, тогда а « 1, а X > 1 и сохраняет постоянное значение. В этом случае переменные х и у изменяются быстро, а переменная г - медленно.

Будем искать решение системы дифференциальных уравнений (9) в виде

0

х = хо +ах^ +а х0 +...;

< у = уо +ау1 + а0у0 +...; 0

г = го + аг^ +а г0 +...,

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

г зависят от времени явно: хо = хо (г), уо = уо (г), го = го (г) и неявно: х^ = х^ (хо,уо, го), у = Л (^ Уо, го ), г1 = г1 Уо, го ) и тд.

Пусть также выполняются равенства

— = Хо + аХ +а0 Х0 +...; Шг о 1 0

Шу = Уо +аУл +а2У0 +...; Шг о 1 0

— = а2у +а020 +..., Шг

т.е. правые части дифференциальных уравнений, описывающих ФАП, так же могут быть представлены в виде функциональных рядов с коэффициентами, которые есть не что иное, как степени малого параметра а. Будем учитывать слагаемые функциональных рядов с нулевого по второе включительно и пренебрегать остальными (нулевое слагаемое функционального ряда для йг/йг, очевидно, равно нулю). Заметим также, что по предположению функции X),Х^,Xо,Уд,У\... зависят от времени только неявно: х0 = Хо (гд), Уд = Уд (гд), 2 = 2Х (гд) ,...

Получим соотношения, связывающие друг с другом переменные х0, у0, г0, х1, у1, г1, х2, у2, г2 и Х0, У0, Х1, У1, 21, Х2, У2, Для этого продифференцируем выражения для х, у, г по времени г:

— = (Х0 + аХх + а2 X2 + ...) + а

дх-

дх,

1 (Х0 +аХ1 + ...) +

(У) +ау + ...) + д1 (а2: + ...)

дуо дг0

+ а

дх2 дх0

дх2

(X) (Уо +...)

дуо

+ ...;

— = (уо + аУх + а2У2 +...) + а

ду1 дх0

(Х0 + аХх + ...) +

д1 (Уо +аУ: + № +...) дуд дго

+ а

£ (Х0 +->+£ (У0 +■•■>

+ ...;

— = (а2т + а2 22 + ...) + а

йг \ 1 2 /

дг дх

1 (Х0 +аХ1 + ...) +

■д1 (У) +аУ + (а2х + ...)

дуо дг0

+ а

дг2 дхо

дг2

(Х0 + -)+дТ (У0 +■) дуо

+ ....

С другой стороны, справедливы дифференциальные уравнения, задающие модель ФАП. Подставим в их правые части соотношений для х, у, г и разложим функции б1пх и вшу в ряды Тейлора в окрестности точек х0 и у0:

— = (р + го +аг! +а0 г0 + ...)-аХ^т хо +8sin уо )-Шг ^ '

-а0Х( cos хо х! + в cos уо у1) + •••; = (в + Лр + го +аг1 +а0 г0 + ...)-аХ^т хо +вsin уо )--а0Х( cos хо х1 +вcos уо у^) + •••;

Шу Шг

Шг Шг

а[со (го +аг1 +...) + (1 -асо)(sinхо +вsinуо)

+ (1 -асо )а( cos хо х1 +вcos уо у1)] + •••.

Приравнивая в правых частях выражений для Шх/Шг, Шу/Шг и Шг/Шг слагаемые, содержащие в качестве сомножителей одинаковые степени параметра а, получаем

Хо = го +р;

Х1 + Хо + Уо = г1 -Х(sinхо +вмПуо);

сх,

о

дуо

Х0 +|х1 Х1 +|х1 У +дх121 +дх0 Хо +дх0 Уо = г0-Х( cos хо х1 уо у1);

ЙКо дуо дго дхо дуо

Уо = го + Р + ЛР;

У 1 Хо + ду1 Уо = г1 -Х^т хо +вsin уо );

сх.

о

дуо

У0 + !1 Х1 +|у1 У + ^ 2Х +ду° Хо +ду° Уо = г0-Х( га хо хх +вcos уо у ); дхо суо дго дхо дуо

21 Хо + ^ Уо =-[сого +(1 -асо)(sinхо +вsinуо)]

<

2 0 Х1 + У1 21 + Х о Уо = -[со го +(1 - асо )(cos хо х1 +в sin уо ух

дхо дуо & о дхо дуо

Теперь необходимо учесть, что если применяется процедура усреднения, то функции х1, у1, г1, хг, уг, г0 периодические (с периодом 2п) по переменным хо и уо, либо зависят от одной из этих переменных. Следовательно,

1()2П 1о ^ х1 ( х0, уо , г0 ) йх0йу0 = 0

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

Проинтегрируем полученные выше уравнения связи функций х0, у0, г0, х1, у1, г1, х2, у2, г2 и Х0, У0, 21, Х1, У1, 21, Х2, У2, 22 в пределах [0;2п] по х0, у0, разделив полученные

равенства на 4п2 . Тогда ясно, что

X о = гд + Р;

XI = 0; X 2 = —XI;

Уо = г о +Р + Ар;

У1 = 0; У2 = —Х1;

[21 = сог о;

I2 2 = —(1 — ас0 )1,

где I = |02п |02п ( СОБ хд х: + в СОБ уо у1 ) йхд йуд.

4п

С учетом полученных равенств уравнения связи упрощаются:

^ (го + Р) + ^"1 (го +Р + Ар) = гх — Х(бш хд + вмп уд ); ОХ) ОУо

^ (го +Р) + |У1 (го +Р + Ар) = г1 — Х(бш хо + вмп Уо ); дхо °Уо

дг дг

^ (го +Р) + тг1 (го +Р + АР) = —(1 — ас) )( мп хо + вбш уд ); дх0 ду0

дх дх дх

—2 (г0 + Р) ++ —2 (г0 +Р + АР) = с0 г0 + г 2 — ^(сОБ х0 х1 + в СОв Уо У1 — 1); дхо °Уо дг0

(го + Р)+(го + Р + АР) = со г0 + г2 — х(соб хо х1 + в СОБ УоУ1 — I);

дх0 дУо дг0

дг дг дг

—2 (г 0 + Р) + —2 (г0 +Р + АР) = с0 г 0 V1 — [с0 г1 +(1 — ас0 ХСО8 хо х1 + в СОВ Уо У1 —1)]

дх0 дУо дг 0

Исключая из трех первых уравнений функцию г1, убеждаемся, что

^ (го +Р)2 (го + Р + АР)2 + 2 (го +Р)( го + Р + АР) =

дх)2 ду2 дходУо

= — X СОБ хд (г) + Р) —Хв СОБ Уо (г) +Р +АР) —(1 — ас) )(зШ хд + В б1П Уо ).

Аналогичное уравнение справедливо для функции у1.

Чтобы I ^ о, необходимо, чтобы функция х1=х1(хо, уо, го) не зависела от уо или по крайней мере содержала слагаемое, не зависящее от уо: х^ = х^ (хо, г о ) +... Тогда умножая последнее дифференциальное уравнение в частных производных на cosх0, интегрируя в пределах [0;0п] по хо и уо, и, наконец, разделив результат на 4п0, получаем

^о * ^ хох1 (хо, уо , го )Шхошуо = |о0П х1 (хо, го)cos хоШхо = 0(0.

Аналогично, чтобы I ^ о, необходимо, чтобы функция у1=у1(хо, уо, го) не зависела от хо или по крайней мере содержала слагаемое, не зависящее от хо: у^ = у^ (хо, го) +... Тогда умножая последнее дифференциальное уравнение в частных производных на cosy0, интегрируя в пределах [0;0 п] по х0 и у0, и, наконец, разделив результат на 4п0, получаем

4п

^ |00п |00п cos уо у (хо, уо, го )Шхо Шуо = 0П |00п Л (уо, го)cos уошуо =-—

2 ( го +Р + ЛР)'

Отсюда имеем

1 = Х

0

г

го +Р го +р + Лр

Таким образом, найдены явные выражения через го для функций Хо, Уо, 21, доказано, что Х1=0, У1=0. Теперь с учетом явного выражения интеграла I через г0 получаем соотношения для функций Х0, У0, 20:

0 (

Х 0 =-Х1 = -—

1 в

- +

0 ^

го +Р го +р + Лр

0 (

У0 =-Х1 = -—

1 в

- +

0 ^

го +Р го +р + Лр

0

0

X

=-(! - ас0 )1 = -(1 - ас0 ) 2

1

z0 +р z0 +р + Ар

Следовательно, можно записать систему дифференциальных уравнений для функций х0, у0, 20, правые части которой согласно сделанным предположениям зависят только от функций Х0, Уо, XI, У1, Z1 и Х2, У2, Z2:

йх0 X

—- = В + ¿0--а

йг 0 2

(

+в ^0 +Р + АР

^=р+ар+¿02 а

йг 0 2

1

¿0 +в ¿0 +Р + АР

йг° = -ас0 ¿0 -(1 - ас0 )—

йг

(

¿0 +в ¿0 +Р + АР

(10)

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

При с0=1 (невырожденный ПИФ) нельзя провести аналитическое интегрирование последнего уравнения усредненной системы (10). Однако при больших значениях постоянной времени ПИФ можно пренебречь величиной с0 и положить с0=0, что соответствует вырожденному ПИФ. Рассмотрим последнее дифференциальное уравнение системы (10) при этом условии:

йг

0=-—а2

(

\

¿0 +в ¿0 +р + Ар

Разделяя переменные и производя интегрирование, находим

2 ( 2 АО ^/-А

г0

+ г0

Р +

^АР

V 1+

(АР)2

1п

1 + в2

АР

го + Р + —^ 1 + в2

+ С = —(1 + в2 )Ха2 г/ 2,

где С - произвольная постоянная. Если положить

V = — (г0 +Р), АР v 0 '

то получим

V2 , ,, ^ Л ?л3 ха2 1

2 а в2 (АР)2

2в

- + ^ — 1п|1 + V| ) + С = —(1 + в2 )3

где С - новая произвольная постоянная, связанная с С соотношением

_ С С = С — 1п

в2

АР

1 + в2

Р(1 + в2

АР

1 + -

Р(1 + в2 2АРв2

Введем также новое время

т = (1 + в 2 Г — а 2—т^—- г.

2 в 2 (АР)2

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

VI

2в:

- + (V — 1п|1 + V)) + С

= —т;

V = в2 (1 + V)

й т

V

V + (1 + в2

Зависимости V (т) могут быть заданы только неявно. Функция —т^) имеет два минимума в точках V = 0 и V = —(1 + в2), а также вертикальную асимптоту в точке

2

V = -1. На рис. 2 приведены графики функций -т(¥) при С = 0 и различных в: 8=0.6 (рис. 2, а), 8=1 (рис. 2, б) и 8=1.4 (рис. 2, в).

Результаты представлены на рис. 3, а-в. Для различных значений С (4; 0; -4) заметим также, что из выражения для производной dV|dт следует, что функция V (т)

возрастает при -го< V <(1+ 82 ), -1 < V < 0 и убывает при (1+ 82 )<V <-1, 0 <V <го.

Направление изменения функции V (т) с ростом т показано на рис. 3 стрелками.

В результате можно сделать вывод, что, задав начальный момент времени Т0 и начальное значение переменной V (т), равное V), можно однозначно определить решение дифференциального уравнения для переменной V (т), соответствующее данным начальным условиям. Процесс поиска такого решения показан на рис. 3 штрихпунктирными линиями. Если V <-1, то процесс изменения переменной V (т)

заканчивается в точке V = -(1+ 82 ) или же ¿0 = -0-Ар, при этом, как следует из второго

дифференциального уравнения усредненной системы, установившееся значение переменной ¿0 компенсирует первоначальную частотную расстройку между помехой и свободными колебаниями генератора, входящего в состав ФАП, т.е. происходит захват за помеху. Если ^ >-1, то процесс изменения переменной V (т) заканчивается в точке

V = 0 (или ¿0 = -Р). При этом, как следует из первого дифференциального уравнения усредненной системы, установившееся значение переменной ¿0 компенсирует первоначальную частотную расстройку между сигналом и свободными колебаниями генератора, т.е. происходит захват за сигнал. Задав начальное значение переменной ¿0 (г),

(0)

равное ¿0 , из неравенства

V» = ^ (40) + э)< 1

-2-'-'-'-'-1-'-'-'-'-

-5 -4 -3 -2 -1 0 1 2 3 4 V

в)

Рис. 2. Графики зависимости функции —т от V при различных значениях 8

можно определить границы области захвата за помеху на плоскости параметров (Р, АР) :

Р<— г00) — -^ при АР> 0; 1 + в 2

Р > — г)0)--при АР > 0.

1 + в 2

Аналогично определяются границы области захвата за сигнал:

Р > —г)0) —АР2 при АР>0; 1 + в2

Р<—г)0)—АР2 при АР <0. 1 + в2

/ / / / Xза сигнал XXXX / 5

ХХХХХХХХХХХХхХХХ/,

С:>С;>Сз Захват за помеху

-20

-15

/// Захват X /за сигнал /

////////////^ УХ^к/

С1>С2>С3 У су Захват уХ за помеху ^ С-^Х

-10

-20

-15

-10

а)

V

3 2 1

б)

жщ Х/^Зах-ват X/ у/ /.за сигна1 /-, ///%/

ххххххххх^ххх/^ у/22//,

С:>С;>С3 Ч Захват _ т за помеху /

-14 -12 -10

4

в)

Рис. 3 Графики зависимости функции -т от V при различных значениях 8 и С

На рис. 4 приведена бифуркационная диаграмма ФАП второго порядка при наличии гармонической помехи на плоскости параметров (в, Ав), если начальное

значение переменной ¿0 (г) равно ¿00).

Из полученной диаграммы следует, что режим работы системы определяется не только частотными расстройками сигнала и помехи (параметры в и ДР), но и начальным значением переменной ¿0. Это говорит о том, что на режим захвата оказывает влияние предыстория работы системы.

Рис. 4. Бифуркационная диаграмма ФАП второго порядка

Из диаграммы также следует, если система находиться в режиме захвата за сигнал г)0) =Р^, то перейти в режим захвата за помеху она не может ни при каких 8 и Др.

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

Проведенное методом Рунге-Кутты четвертого порядка численное решение системы (9) показало, что при определенной интенсивности помехи происходит ее захват. Поиск критических значений Др проводился при условии, что система осуществляет

слежение за сигналом ^ г)0) = Р^ и при отсутствии сигнальной частотной расстройки

(Р=0). На рис. 5 квадратными метками показаны критические значения Др, при которых происходит смена режима захвата.

При малых значениях Др (~0,1) и близких амплитудах сигнала и помехи (8~1) в системе возникают асинхронные движения (заштрихованная область на рис. 5). ФАП не может осуществить захвата ни помехи, ни сигнала, а осуществляет синхронизацию на некоторой промежуточной частоте.

С учетом полученных результатов бифуркационная диаграмма при 8=2 принимает вид рис. 6. Здесь в отличие от рис. 4 имеется участок, на котором происходит захват

помехи (при 8>1) при любых значениях Р и г.

(0) 0 .

Е/|\

я

о 0;2 0;4 0=6 0:8 др

Рис. 5 График зависимости 8 от в

Рис. 6 Бифуркационная диаграмма ФАП

Характер движений, возникающих в системе при различной интенсивности помехи, иллюстрируют рис. 7-8, где изображены проекции фазовой траектории на плоскости (х,у), (х,г), (у,г). Рис. 7 соответствует режиму захвата за сигнал (Дв=0,1; 8=0,6; а=0.8; а=0.16), рис. 8 - режиму захвата за помеху (Дв=0,1; 8=1,4 а=0.8; а=0.16).

а) б) в)

Рис. 7. Проекция фазовой траектории на плоскости в режиме захвата сигнала

а) б) в)

Рис. 8. Проекция фазовой траектории на плоскости в режиме захвата помехи

Используемый численный метод решения системы дифференциальных уравнений позволяет получить зависимости, аналогичные рис. 5, в случае невырожденного ПИФ в кольце ФАП (с0=1). Параметры системы аналогичны случаю вырожденного ПИФ (а=0,05; а=0,01). Сравнивая рис. 9 с рис. 5, можно заметить, что характер зависимости практически не изменился, отличие состоит лишь в размерах области неустойчивой синхронизации при малых Дв=0,1 и 8~1. Это подтверждает, что при больших Тф можно без больших потерь рассматривать случай вырожденного ПИФ в кольце ФАП.

%

В 6 4 2-

О

Рис. 9. График зависимости 8 от в

Рассмотрим влияние полосы пропускания ПИФ на процессы захвата сигнала или помехи. На рис. 9 треугольными метками показаны результаты численных расчетов при а=0,01; а=0,005, что соответствует увеличению постоянной времени ПИФ в пять раз по сравнению с предыдущем случаем. Из полученных зависимостей следует, что увеличение постоянной времени фильтра (уменьшение полосы пропускания) повышает помехоустойчивость системы. Однако, при малых Дв достаточно незначительного

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

В заключение рассмотрим определение времени переходного процесса, т.е. промежутка времени, за которой ФАП компенсирует первоначальную частотную расстройку и произойдет захват за сигнал либо за помеху.

Пусть V) > 1 или, что то же самое, (1 + 82 )(г)0) + р)/Лр > -1. Тогда конец процесса регулирования частоты - это достижение точки V=0 и захвата за сигнал. Зависимость V (т), соответствующая такому процессу, соединяет точки (то,^) и (тк,^к). При этом тк -то =Дт - длительность переходного процесса; Vk = 0, тд, V) и т^, Vk - соответственно начальные и конечные значения переменных т и V. Так как обе точки лежат на одной и той же кривой, то

V?

2в:

VI

2в:

- + - 1п 1 + V)|) + С = -то;

- + - 1п 1 + Vk\) + С = -тк.

Вычитая из первого равенства второе, имеем длительность переходного процесса при захвате за сигнал

Лт = + (Vo- 1п|1 + V)|); V) >-1.

2в

Аналогично из равенств

У2

2в2

2в2

■ + (Vo - 1п 1 + V)|) + С = -то;

■ + (Vk - 1п 1 + Vk|) + С

(1+82 )2

+ С -21п8-11 + 82 ) = -

28 2

"тк,

вытекающих из того факта, что при V) <-1 и захвате за помеху зависимость V(т), соответствующая процессу регулирования частоты, соединяет точки (то,Р0) и (тк,Рк),

причем Vk =-(1+ 82), может быть получена длительность переходного процесса при

захвате за помеху:

Лт = 1Р + 8) + 1п

1 + V0

; V) <-1.

Приложение

Найдем критические значения первой гармоники х1=х1к и отношения помеха/сигнал е=£к, которые определяют условия срыва синхронизации:

Бт X) =±1; соб X) = 0.

(П1)

При учете второго условия по уравнению

х2 = 8 2М 2 32

Х1 = ^ , ^2

8 2 Т 2 8 3 0

+ Р2

й

СОБ X) БШ Р

2 31 +

^ йХу

М

СОБ Р

М

находим х2 = М 21й2 . Отсюда следует

82 = & VМ 2.

При учете первого равенства в (П.1 ), находим

Х 2 = Х1к "

+1

\м о ;

2М

й СОБ Р

Подставляя это значение х^ в (П.1) окончательно получим

2

8

8 к =

2й

у А

М соб Р

■ +1

Vм 0 ;

2й

М соб Р

7

V М 0

+ й

Зависимость 8^ = / (й) при а=0.8 и а = 6.25 изображена на рис. П.1, где прямая

1 получена при ФАП 1-го порядка и в=0; прямая 2 при невырожденном и вырожденном фильтре и в=0; 3 - ФАП 1-го порядка в=0.8; 4 - невырожденный и вырожденный фильтр и в=0.8;

Рисунок П.1

По рис. П. 1 замечаем, что при Р>0 и ё>0 помехоустойчивость ФАП выше, чем Р>0,

а ё<0.

Можно показать, что условиями захвата за сигнал при действии гармонической помехи является [5] условие

или условие

\р + 8\< 1; 8-в< 1,

-(1 + е)<Р< 1 -8; 8-1 <в< 1 + 8.

С другой стороны, по [9] находим

8 < 1 - в; 8> 1+ Д

или

s-1 <в< 1 -s; 1 -S<P<S-1.

Причем во втором случае границы захвата по сигналу совпадают с [4] несмотря на то, что получены [4,5] и [9] на основе разных критериев.

Достаточным условием захвата за помеху является условие [5]

вж_ + 1 < 1; 1

s s s s

< 1,

или

-(1 + s) < вп < s — 1; 1 — s < вn <s +1.

При этом левая граница в данном случае совпадает с [4].

Заключение

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

Список литературы

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

2. Meyr H., Ascheid G. Synchronization in digital communications frequency - locked loops, and amplitude control. N.Y.: Wiley, 1990. 510 p.

3. Stephens D.R. Phase - locked loops -for Wireless communications. Digital, analog and implementations. 2nd ed. N.Y.: Kluwer, 2002. 421 p.

4. Шахтарин Б.И. Анализ систем синхронизации методом усреднения. М.: Радио и связь, 1999. 495 с.

5. Nakagawa M. Effects of interfering signals in phase-locked loops // Frequentz. 1978. Vol. 32, no. 5. P. 146-153.

6. Шахтарин Б.И. Квазигармонический метод и его применение к анализу нелинейных фазовых систем. М.: Энергоатомиздат, 1987. 192 с.

7. Быховский М.А. Влияние помехи на процессы захвата в системе фазовой автоподстройки частоты // Радиотехника и электроника. 1987. № 10. С. 2131-2141.

8. Richman D. Color-carrier reference phase synchronization accuracy in NTSC color television // Proc. IRE. 1954. Vol. 42, no. 1. P. 106-133.

9. Karsi M.F., Lindsey W.C. Effects of CW interference on phase-locked performance // IEEE Trans. 2000. Vol. COM-48, no. 5. P. 886-896.

SCIENTIFIC PERIODICAL OF THE BAUMAN MS TU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Influence of harmonic noise on the second order phase-lock loop

# 05, May 2013

DOI: 10.7463/0513.0551477

Aslanov T.G.

Bauman Moscow State Technical University, 105005, Moscow, Russian Federation

tabasik@gmail.com

Influence of harmonic noise on phase-lock loop (PLL) leads to distortion of received information and a loss of PLL stability. This article analyzes effects of harmonic interference on synchronization processes in the second order PLL system. Analysis of signal grabbing and interference capturing cases was carried out in this work; projections of a phase path onto planes, which illustrate the specified modes, were also obtained. By analyzing phase planes, equations of capturing signal and interference were obtained. Dependence diagrams of the signal-to-noise ratio and frequency error are presented along with a bifurcation diagram. Application compares critical values of the PLL parameters obtained by different methods.

Publications with keywords: error performance, harmonic interference, the relation a noise/signal

Publications with words: error performance, harmonic interference, the relation a noise/signal

References

1. Shakhtarin B.I. Analiz sistem sinkhronizatsiipri nalichiipomekh [Analysis of the systems of synchronization at presence of noise]. Moscow, IPRZhR, 1996. 252 p.

2. Meyr H., Ascheid G. Synchronization in digital communications frequency - locked loops, and amplitude control. N.Y., Wiley, 1990. 510 p.

3. Stephens D.R. Phase - locked loops -for Wireless communications. Digital, analog and implementations. 2nd ed. N.Y., Kluwer, 2002. 421 p.

4. Shakhtarin B.I. Analiz sistem sinkhronizatsii metodom usredneniia [Analysis of the systems of synchronization by method of averaging]. Moscow, Radio i sviaz', 1999. 495 p.

5. Nakagawa M. Effects of interfering signals in phase-locked loops. Frequentz, 1978, vol. 32, no. 5, pp. 146-153.

6. Shakhtarin B.I. Kvazigarmonicheskii metod i egoprimenenie к analizu nelineinykh fazovykh system [Quasiharmonic method and its application to the analysis of non-linear phase systems]. Moscow, Energoatomizdat, 1987. 192 p.

7. Bykhovskii M.A. Vliianie pomekhi na protsessy zakhvata v sisteme fazovoi avtopodstroiki chastoty [The influence of noise on the capture processes in the system of phase-locked loop]. Radiotekhnika i elektronika, 1987, no. 10, pp. 2131-2141.

8. Richman D. Color-carrier reference phase synchronization accuracy in NTSC color television. Proc. IRE, 1954, vol. 42, no. 1, pp. 106-133.

9. Karsi M.F., Lindsey W.C. Effects of CW interference on phase-locked performance. IEEE Trans.Commun., 2000, vol. COM-48, no. 5, pp. 886-896.