Статистические характеристики дискретных систем фазовой синхронизации Текст научной статьи по специальности «Математика»

научное издание мгту им. н. э. баумана

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

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

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

Статистические характеристики дискретных систем

фазовой синхронизации

# 08, август 2012

Б01: 10.7463/0812.0475421

Ковальчук А. А.

УДК: 621.396.662

Россия, МГТУ им. Н.Э. Баумана nastia kov-k@rambler.ru

Введение

Теория стохастических систем фазовой автоподстройки (ФАП) начала формироваться работами Р.Л. Стратоновича [1] и В.И. Тихонова. Одним из основных методов исследования является метод марковских и полумарковских случайных процессов и цепей, разработанные А.Н. Колмогоровым и Л.С. Понтрягиным [2]; применяются интегральное уравнение Колмогорова-Чепмена и уравнение среднего времени до срыва синхронизации (работы Ч. Цзы [3], А. Вайнберг и В. Ли [4], В.И. Витюцкого и В.Н. Сердюкова [5]). В указанных работах найдены плотности распределения вероятности (ПРВ) координат, в меньшей степени исследованы характеристики срыва слежения, причем рассмотрены детерминированные системы, и не проводится анализ стохастической системы. Так же среди систем ФАП особое внимание уделяется непрерывным системам [6], а статистические характеристики дискретных ФАП до сих пор исследованы в меньшей степени, отчасти в силу многообразия дискретных систем.

В данной статье приводится анализ стохастической дискретной ФАП при синусоидальной характеристики детектора. Приводятся зависимости стационарных ПРВ при нулевой частотной расстройке, рассчитанные методом Галеркина. Приведен приближенный метод вычисления ПРВ на основе метода Галёркина. Сопоставлены результаты точных и приближенных методов. Делается вывод о том, что система алгебраических уравнений, полученная методом расщепления ядра интегрального уравнения, имеет тот же вид, что и система, полученная методом Галёркина. Приводятся графики нестационарной ПРВ при равномерном начальном распределении. Исследованы моменты времени до срыва слежения в дискретной ФАП. Приводятся зависимости среднего значения времени до срыва слежения

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

1. Анализ стохастической дискретной ФАП при синусоидальной характеристике детектора

Рассмотрим статистическую динамику дискретной ФАП при синусоидальной характеристике детектора g (x) = sin x. Разностное уравнение ФАП в этом случае принимает вид

Хь

= Хк _ T0 (sin Хк _Р) + ñk, (1)

Ч+1

причем среднее значение равно

Е (Хк+1| Хк ) = Е (Хк+1| *) = * - то (sin *-Р),

где 7д - период, в - частотная растройка, Пк - гауссова последовательность. Уравнение Колмогорова-Чепмена запишем в виде [7]

•п

)-п

где

Wm+1 =\Я_яЧ1 ( Х^ )Wm ( Z ) dz, (2)

(x | z) = ^ q (x + 2nn | z); q (x | z) = [ 1/\¡2na2 ] exp{_[x _ z + T) (sin z _в)]2/2<г2}.

W(x) = Пq(x | z)W(z)dz.

qi\

n=-<ю

В стационарном режиме уравнение Колмогорова-Чепмена переходит в интегральное уравнение Фредгольма второго рода

•п

>-п

Предположим, что, как и в случае непрерывной системы при в = 0, сигнал рассогласования в стационарном режиме распределен по закону Тихонова [8] при r = р, где р = r0 (1- T0 / 2). Тогда

1 1 ю

W (x) =-— ерcosx = — + У an cos nx, (3)

V > 2nIо (v) 2п П=1 n

где an = nI0 (р)In (р); In (р) - функция Бесселя.

Рассмотрим влияние интервала дискретизации сигнала на плотность распределения вероятности (ПРВ) W ( x ). Вычисление интеграла в правой части уравнения Колмогорова-

Чепмена проведем при T0 = 0,5 и T0 = 1. Для сравнения с [7] составим перечень значений 10.7463/0812.0475421 264

отношения сигнал/шум г = = / А2, где <7^ = 72 / 7. Пусть Т = 0,5; тогда значениям аг = 3; 1; 0,6; 0,2 соответствуют р = 1; 3; 5; 15 и г = 4/3; 4; 20/3; 20. При 7 = 1 для тех же значений а1 имеем р = 1/3 (1); 1 (2); 5/3 (3); 5 (4) и г = 2/3; 2; 10/3; 10 (рис. 1).

Щх) Щх)

0.3

4

/ 3

V 2

У 1

1 Т

/ ] V

/ /А

1 з

1 у2

1Ж / 1

■ и / V /1 \

\ /4

1 ^ \

<1/2

/Л\____ 4 и/1

у

х/тг

1

о 0.5 'х/л 0 °-5 1 х/л

а) 7=0,1 б) 7=0.5

Щх)

Рисунок 1 - Стационарная ПРВ Ж (х) при в = 0 и различных значениях 70

Результаты сравнения видны на рис.1 а-в, где приведены зависимости стационарных ПРВ Ж от х при в = 0 и значениях 7, равных 0,1; 0,5; 1. Сплошными линиями обозначены

результаты расчета методом Галеркина, штриховыми линиями - ПРВ W ( x ) (3) непрерывной системы.

2. Вычисление статистических характеристик ФАП приближенным методом

Рассмотрим приближенный метод вычисления ПРВ Ж (х) на основе метода Галёркина [8]. Пусть

N

ж(х)« wN (х), wN (х) = x сй (n)¥п (х),

п =0

где { у/п(х) } (п=0,1...) - полная система ортогональных на интервале (—п,п) функций. Точное значение ПРВ находим в форме предельного соотношения

W (х )= lim WN (х).

N ^да

Тогда

да 1 да

W (х) = Z CnVn (х) = 2- + Z (An C0s пх + Bn sin nX) (4)

п=0 2П п=0

где cn = lim cn (N), An, Bn - расчетные коэффициенты.

N ^да

Приведенные соотношения для коэффициентов ряда (4) могут быть получены методом замены ядра интегрального уравнения Колмогорова-Чепмена на вырожденное

разложение ядра qi (х | z) в двойной ряд Фурье.

Рассмотрим анализ нестационарного режима методом Галёркина. Плотность распределения вероятности

N

Wk ( х) = Wk ( х 1 х0 )= lim Z dn (k )Vn (X )

N ^да n=0

где х0 - начальное значение фазового рассогласования.

Согласно уравнения Колмогорова-Чепмена и условию W (х) = &(х — х0 ) имеем

да

W1 (х | х0 ) = q(х | z — х0 ). Отсюда W1 (х) = Z dn (l)¥n (х).

n=0

Подставляя ряд для функции qi (х | z) в интегральное уравнение (2), получаем выражение

Wm +1 (x )=X Ak (m+1)cos kx+Bk (m+1)sin kx,

k=0

где

ж

Ак (m + 1) = X ank £Wm ( z) c0s nzdz + Cnk ¡п Wm ( z) sin nzdz, n=0

ж

Bk ( m + 1)= X bnk ¡U_nWm ( z ) c0s nzdz + dnk ¡U_nWm ( z ) sin nzdz. (5)

n=0 П _

В частном случае в = 0 выражение для стационарной ПРВ сигнала рассогласования в дискретной системе принимает вид

1 ж

W(x) = — + X An cosnx. (6)

2 n =1

Ограничимся первым слагаемым ряда (6). Тогда

А = - e~°2/2 JT0

1 _ ь 2 ' п 1-(2/То)

Результаты расчета ПРВ ж(х) приведены сплошными линиями на рис. 1 а-в. На рис. 2 и рис. 4 изображены графики при в Ф 0 . На рис. 2 изображены графики ПРВ при в = 0,4, р = 4; Т0 = 0(1); Т0 = 0,5 (2); Т = 0,1 (3); штриховой линией отмечена ПРВ для непрерывной системы. На рис. 3 изображены графики нестационарной ПРВ при равномерном начальном распределении Ж^)(Х) = 1/2 п, когда в = 0, Р = 4, Т = 1; /=1

(1), /=2 (2), /=4 (3), /=10 (4); штриховой линией отмечена ПРВ для непрерывной системы. На рис. 4 а,б принято, что р = 0,3 (1); р = 1 (2); р = 3 (3); р = 5 (4). На рис. 4 а: в = 0,5; Т0 = 0,1; на рис. 4 б: в = 0,5; Т0 = 0,5.

\У(х)

0.5

-0.5

0.3

Рисунок 2 - Стационарная ПРВ Щх) при /5 = 0,4; р= 4

¡¥(х)

3\ V* А\ /1

\ 4 1 1 / 1 1

г 1 1

/ ? // л 1 \

-,/4

ь* Ж3

. / 2 М

щ \\ 1

/ ! § ,/'// ■ / л \\ ^

- - - / у ' V чЛО^:- - -

! х/тг

0.3

1 х/я

Рисунок 3 - Нестационарная ПРВ при равномерном начальном распределении

\¥(х)

0Ух)

0.8

0.6

0.4

0.2

/V-

п

а I ! 2 1

^Л \\ \

0.5

0.6

6.4

1

Л/ 3

/ \

/ 1 (у 2 1

-05

0.5

1 х/л

-0.5

0.5

а) 7,\ = 0,1

б)Го = 0,5

Рисунок 4 - ПРВ Ж(х) при р Ф 0

1 х/я

3. Моменты времени до срыва слежения в дискретной ФАП

Для вероятности срыва слежения на г-м шаге можно записать [9] в виде

Сг (х0 ) = | q (х | х0) di-1 (х) dx; dl (х0 ) = 1 q (х | х0) сХ. (7)

Определим к-й начальный момент числа шагов до срыва синхронизации (Б = 2п)

т,

ж ж

(хо) = Ё^ (хо) =1 +Я(* I хо )\ Ё[гкй -1 (х)

/ =1 I/ = 2

-1

йх.

Получаем рекуррентные соотношения для начальных моментов числа шагов до срыва синхронизации

ть

ж

(хо) =1 + ¡-5ткЯ(х 1 хо)йх + ЁС'ктк-/ (хо )

(8)

/ = 2

п-1

где Ск - биномиальный коэффициент, тп (*о ) = тп (*о ) — 1 — Ё С'^к — (Хо ).

/=1

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

На рис. 5 приведены зависимости среднего значения времени ус = т^о до срыва

слежения при Хо = о, в = о и различных значениях р и То: (1) - То =о; (2) - То =о,1; (3) -

То =о,5; (4) - То =1 (штриховой линией отмечены данные из [1о]).

Рисунок 5 - Среднее значение времени до срыва слежения ус (р) при *о = о; в = о

Заключение

Проведен анализ стохастической дискретной ФАП при синусоидальной характеристики детектора. Приведены зависимости стационарных ПРВ Ж (х) при в = о , рассчитанные методом Галеркина. Рассмотрен приближенный метод вычисления ПРВ Ж(х)

на основе метода Галёркина. Сделан вывод о том, что система алгебраических уравнений, полученная методом расщепления ядра интегрального уравнения, имеет тот же вид, что и система, полученная методом Галёркина. Приведены графики нестационарной ПРВ при

равномерном начальном распределении Wq (x). Исследованы моменты времени до срыва

слежения в дискретной ФАП. Приведены зависимости среднего значения времени до срыва

слежения yc от р.

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

1. Стратонович Р.Л. Избранные вопросы теории флуктуации в радиотехнике. М.: Советское радио, 1961. 210 с.

2. Андронов В.А., Витт А.А., Понтрягин Л.С. О статистическом рассмотрении динамических систем // Журнал Экспериментальной и Теоретической Физики (ЖЭТФ). 1933. Т. 3, № 3. С. 165-180.

3. Chie C.M. Mathematical analogies between first-order digital and analog phase-locked loops // IEEE Trans. on Communications. 1978. Vol. COM-26, № 6. P. 860-865. DOI: 10.1109/TCOM.1978.1094148

4. Weinberg A., Liu B. Discrete time analysis of nonuniform sampling first- and second-order phase-locked loops // IEEE Trans. on Communications ,1974. Vol. COM-22, № 2. P. 123-137. DOI: 10.1109/TCOM.1974.1092168

5. Битюцкий В.И., Сердюков П.Н. Оценка времени до срыва синхронизма в импульсной системе ФАПЧ // Радиотехника. 1973. № 8. С. 95-97.

6. Ковальчук А. А., Сидоркина Ю. А., Рязанова М. А. Воздействие на систему синхронизации гармонических помех и шума // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2011. № 3. Режим доступа: http://technomag.edu.ru/doc/169529.html (дата обращения 19.10.2012).

7. Канторович Л.В., Крылов В.И. Приближенные методы высшего анализа. М.-Л.: Гостехиздат, 1949. 695 с.

8. Шахтарин Б.И. Анализ кусочно-линейных систем с фазовым регулированием. М.: Машиностроение, 1990. 192 с.

9. Шахтарин Б.И., Сизых В.В., Курочка Б.Я. Исследование статистических характеристик дискретных ФАС первого порядка // Вестник МГТУ. Сер. Приборостроение. 1992. №3. С. 89-110.

10. Первозванский А.А. Случайные процессы в нелинейных автоматических системах. М.: Наука, 1962. 352 с.

scientific periodical of the raijman ms tu

SCIENCE and EDUCATION

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

electronic scientific and technical journal

Statistical characteristics of discrete systems of phase locking

# 08, August 2012

DOI: 10.7463/0812.0475421

Kovalchuk A.A.

Russia, Bauman Moscow State Technical University

nastia kov-k@rambler.ru

This paper provides statistical characteristics of discrete systems of first-order phase locking in the absence of noise, in contrast to most other studies which consider deterministic systems. One of the main research methods is the method of Markov and semi-Markov random processes and circuits developed by A.N. Kolmogorov and L.S. Pontryagin. The theory of stochastic systems of phase locking first emerged in works by R.L. Stratonovich and V.I. Tikhonov. Among the phase locking systems the main focus is on continuous systems whereas statistical characteristics of discrete phase locking have been studied to a lesser extent, partly because of diversity of discrete systems. The author gives an analysis of stochastic discrete phase locking with sinusoidal characteristics of the detector. The results of exact and approximate methods were compared for the first time. As a result, the main statistical characteristics were obtained: probability density of phase mismatch in the stationary and nonstationary modes, time left till tracking loss in pulse phase systems.

Publications with keywords:tracking failure, discrete systems of phase synchronization, a sinusoidal response detector, the probability density of the phase mismatch, the Galerkin method, the failure of the tracking method of Markov random processes

Publications with words:tracking failure, discrete systems of phase synchronization, a sinusoidal response detector, the probability density of the phase mismatch, the Galerkin method, the failure of the tracking method of Markov random processes

References

1. Stratonovich R.L. Izbrannye voprosy teorii fluktuatsii v radiotekhnike [Selected questions of the theory of fluctuations in radio engineering]. Moscow, Sovetskoe radio, 1961. 210 p.

2. Andronov V.A., Vitt A.A., Pontriagin L.S. O statisticheskom rassmotrenii dinamicheskikh system [The statistical analysis of dynamic systems]. Zhurnal Eksperimental'noi i Teoreticheskoi Fiziki (ZhETF) [Journal of Experimental and Theoretical Physics), 1933, vol. 3, no. 3, pp. 165-180.

3. Chie C.M. Mathematical analogies between first-order digital and analog phase-locked loops. IEEE Trans. on Communications, 1978, vol. COM-26, no. 6, pp. 860-865. DOI: 10.1109/TCOM.1978.1094148

4. Weinberg A., Liu B. Discrete time analysis of nonuniform sampling first- and second- order phase-locked loops. IEEE Trans. on Communications, 1974, vol. COM-22, no. 2, pp. 123-137. . DOI: 10.1109/TCQM.1974.1092168

5. Bitiutskii V.I., Serdiukov P.N. Otsenka vremeni do sryva sinkhronizma v impul'snoi sisteme FAPCh [Estimate of the time to stalling of synchronism in pulsed PLL]. Radiotekhnika, 1973, no. 8, pp. 95-97.

6. Koval'chuk A. A., Sidorkina Iu. A., Riazanova M. A. Vozdeistvie na sistemu sinkhronizatsii garmonicheskikh pomekh i shuma [Affecting system of synchronization of harmonic hindrances and noise]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2011, no. 3. Available at: http://technomag.edu.ru/doc/169529.html , accessed 19.10.2012.

7. Kantorovich L.V., Krylov V.I. Priblizhennye metody vysshego analiza [Approximate methods of higher analysis]. Moscow-Leningrad, Gostekhizdat, 1949. 695 p.

8. Shakhtarin B.I. Analiz kusochno-lineinykh sistem s fazovym regulirovaniem [Analysis of piecewise linear systems with phase regulation]. Moscow, Mashinostroenie, 1990. 192 p.

9. Shakhtarin B.I., Sizykh V.V., Kurochka B.Ia. Issledovanie statisticheskikh kharakteristik diskretnykh FAS pervogo poriadka [Study of statistical characteristics of the discrete phase automatic systems (FAS) of the first order]. Vestnik MGTU im. N.E. Baumana. Ser. Priborostroenie [Herald of the Bauman MSTU. Ser. Instrument Engineering], 1992, no. 3, pp. 89-110.

10. Pervozvanskii A. A. Sluchainye protsessy v nelineinykh avtomaticheskikh sistemakh [The random processes in nonlinear automatic systems]. Moscow, Nauka, 1962. 352 p.