CC BY
24
1К1 Инженерный вестник Дона. №4 (2015) Н| ivdon.ru/ru/magazine/arcliive/n4p2y2015/3490
Оптимизация метода измерения угла фазового сдвига между двумя квазигармоническими сигналами
В.К. Игнатьев, Д.А. Станкевич Волгоградский государственный университет
Аннотация: В статье предложена модификация ранее описанного метода, позволяющая сократить количество вычислительных операций и повысить точность оценивания угла фазового сдвига, которая приближается к пределу, установленному неравенством Рао-Крамера. Также в статье приведены результаты численного моделирования систематических и случайных погрешностей параметрического оценивания угла фазового сдвига предложенным методом.
Ключевые слова: угол фазового сдвига, квазигармонический сигнал, неравенство Рао-Крамера, параметрический метод.
Введение
Измерение угла фазового сдвига необходимо в задачах радиолокации и навигации, радиоастрономии, радиофизике и в других прикладных областях. Точность измерения угла фазового сдвига существенно зависит от наличия аддитивных и мультипликативных помех, скорости и диапазона изменений огибающей и частоты сигналов. Задача измерения угла фазового сдвига давно исследуется, и разработано множество способов ее решения [1 - 3]. В работах [4 - 7] описаны методы оценивания угла фазового сдвига между двумя квазигармоническими сигналами. Особенностями методов является их устойчивость к частотной и амплитудной модуляции сигналов.
Описание метода
В методе из работы [7] используются комбинации перекрестных произведений отсчетов двух квазигармонических сигналов х^) = +
Фо) + П1(0 и Х2(0 = й^п(0(О) + П2(0:
В1 (() = х1 (г - 4Д)х2 (г) - х1 (г)х2 (г - 4Д), В2 (() = х1 (( - 3Д)х2 (( - Д) - х1 (( - Д)х2 ((- 3Д), В3 (() = х1 (( - 3Д)х2 (( - Д) - х1 (( - 2Д)х2 (( - 2Д), В4 (() = х1 (( - Д)х2 (( - 3Д) - х1 (( - 2Д)х2 (( - 2Д),
где обозначено Пь П2 - шумовые составляющие первого и второго сигнала, соответственно, А - некоторый временной интервал, такой, что << 1,
= Показано, что из этих комбинаций можно построить оценку угла фазового сдвига между двумя сигналами ф0:
1в(ф0 )«± Вз - В
2В2 - В, Пч
2ВТГВ1+х0, (1)
Вз + В4^
где %0 - флуктуационная составляющая, которая определяется шумовыми составляющими ^(0 и ^2(0 сигналов, а также скоростью изменения амплитуды и мгновенных значений частот этих сигналов. В выражении (1) необходимо использовать знак «+», если сумма в числителе отрицательна и «-» - если положительна.
Оценка (1) угла сдвига фазы квазигармонических сигналов получается по пяти отсчетам, что при наличии шума приведет к значительному смещению и большой неопределенности оценки. Аналитическое исследование и численное моделирование показывают, что х0 осциллирует с частотой ю, а среднее значение на интервале Т >> 2л/ю близко к нулю. В работе [7] для уменьшения влияния шумов предлагается использовать метод наименьших квадратов совместно с цифровым полосовым фильтром, однако при малом отношении сигнал/шум применение фильтра приводит к значительному смещению оценки при приемлемой ее неопределенности.
Детальный анализ показал, что каждая из комбинаций В1-В4 содержит, кроме полезной постоянной составляющей, которая используется для формирования оценки, аддитивные шумовые компоненты в виде перекрестных произведений отсчетов ц^) и ^2(0 и гармонических функций со случайной амплитудой. Отсюда следует, что для увеличения точности метода необходимо использовать средние значения величин В1-В4:
Щ - В} (2)
2Вг + В/ К)
Iljfl Инженерный вестник Дона. №4 (2015) Н| ivdon.ru/ru/magazine/arcliive/n4p2y2015/3490
В этом случае, если отсчеты шумовой составляющей сигнала дельта-коррелированы, аддитивные шумовые компоненты будут в среднем равны нулю.
Результаты
Наиболее простым критерием качества оценок параметров сигналов при априорной неопределенности этих параметров являются неравенства Рао-Крамера. Сравним погрешности описанного метода и методов [4], [7] с пределом Рао-Крамера, для этого рассмотрим задачу определения параметров двух синхронно дискретизированных сигналов (t = nAt, n = 0,..., N - 1) с одинаковыми частотами f и фазовым сдвигом ф0:
x1 [n] = a1 sinilnfn + ф0) +1 [n] x2 [n] = a2 sin(2nfn) + r 2 [n]
Будем также предполагать, что амплитуды сигналов можно считать постоянными на интервале наблюдения, тогда функция распределения для пары последовательностей из N не коррелирующих между собой отсчетов сигналов имеет вид:
P 1
P=т2-) exP
N-1 N-1
X (x1 [n] - a1 sin(2^fn + ф0 ))2 X (x2 [n] - a2 sin(2^/n))2
2a12 2a 22
Здесь o1, o2 - дисперсии нормальных шумовых последовательностей n1[n] и n2[n] соответственно. Коэффициенты информационной матрицы Фишера [8], рассчитанные по этой функции распределения отсчетов, имеют вид:
/а2 \ i N-1 / -2 \ i N-1
B^ =Ч-Гт ln P) = — X sin2 (2nf + Ф0), B^ =---7 ln P) = — X sin2 (nfn),
\da1 / — n=0 \-a2 / a2 „=0
/-2 , n\ 4a12n2 N-1 22( ) 4a22n2 N-1 2 2( , )
B/f =4-^7lnP = Xn2cos2(2П/„ + Ф0) + —Xn2cos2(2n>),
\ — / — n=1 - 2 n=1
/— 2 \ a n N-1
Bf, = - f- ln P = -П X n2 sin ((nfn + 2ф0),
,-foa1 / a
2
1 n=1
n=0
n=0
:
I?
В/аг \я/Яа.
а п ы~х Я2 \ а N-1
1п Р > = а2ПIП2 з1п/ Вфоа1 = - -Я— 1п Р = I *ш(4я> + 2ф0 ),
/_Я1
В/ф0 \Я/Эф0
1пр) = ^Iпсоз2(2п/ + ф0), Вф0ф0 = -/ 1пР ) = %Iсоз2(2п/п + ф0),
/ О, п=1 \Яф0 / а, п=0
2
а2 п=1
2 N-1
-ф0-а1
-2
2а1 п=0
2 N-1
2
1 п=1
-ф0
2
1 п=0
Неравенство Рао-Крамера [8, 9] для /-го (аь а2, /, ф0 - 1...4) оцениваемого параметра имеет вид:
а,2 >
В аа 0 Ва М В ф0а1
0 В агаг В/аг 0
В /а\ В/аг Ваа В ф0
В ф0а 0 В /ф0 В ф0ф0
Численное моделирование для оценки точности метода (2) было проведено на дискретизированных сигналах, вычисления производились с использованием арифметики с плавающей запятой двойной точности. Оценки угла фазового сдвига строились по N = 100001 отсчетам при различном отношении сигнал/шум, / = 0,1, ф = 0,15. Точность, а именно среднеквадратичное отклонение оценки от среднего значения о и модуль среднего смещения от заданного значения Ь, оптимизированного метода сравнивалась с точностью методов [4] и [7].
При этом перед вычислениями оценок по методам [4] и [7] сигналы обрабатывались цифровыми полосовыми фильтрами с центральной частой /, построенными по методу частотной выборки [10] с полосой 210-/0. Для минимизации влияния переходного процесса, вызванного применением фильтра, из начала выборки профильтрованного сигнала отбрасывалось число отсчетов, равное длительности импульсной характеристики фильтра.
Численное моделирование методов проводилось по 50 реализациям сигналов, результаты приведены в таблице 1. По результатам численного моделирования можно сделать вывод, что предложенный метод не уступает в точности описанным ранее методам при большом отношении сигнал шум и превосходит их - при малом.
1К1 Инженерный вестник Дона. №4 (2015) Н| ivdon.ru/ru/magazine/arcliive/n4p2y2015/3490
Стоит отметить, что методу требуется значительно меньше вычислительных операций, поскольку вычисление квадратного корня и арктангенса требуется производить один раз на выбранном интервале наблюдения, и методу не требуется предварительная фильтрация сигналов.
Таблица № 1
Граница Рао-Крамера, СКО оценки о и модуль среднего смещения от заданного значения Ь, оптимизированного метода (2) и методов [4], [6]
Сигнал / шум, дБ Метод [4] и ПФ Метод [7] и ПФ Метод (2) Граница Рао -Крамера
b о b о b о
80 10-7 8 10-7 10-7 8 10-7 10-7 610-7 5,6 10-7
60 410-7 710-6 5 10-7 710-6 3 10-7 610-6 5,610-6
40 5 10-6 710-5 10-6 710-5 5 10-7 610-5 5,610-5
20 810-4 710-4 410-4 710-4 410-5 610-4 5,610-4
0 610-2 210-2 410-2 10-2 610-4 10-2 5,610-3
Заключение
В статье предложена оптимизация ранее описанного метода. Модификация метода позволяет сократить количество вычислительных операций и повысить точность оценивания угла фазового сдвига, которая приближается к пределу, установленному неравенством Рао-Крамера.
Теоретическая часть исследования выполнена при поддержке гранта РФФИ № 15-47-0229715. Численное моделирование исследования выполнено за счет гранта Российского научного фонда № 15-19-00028.
Литература
1. Massimo, C. 6/12-channel synchronous digital phasemeter for ultrastable signal characterization and use // Frequency Control Symposium & the European Frequency and Time Forum (FCS), 2015 Joint Conference of the IEEE International.
2. Kawagoe, J., Kawasaki, T. A New Precision Digital Phase Meter and Its Simple Calibration Method // IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, VOL. 59, NO. 2, FEBRUARY 2010.
3. Дюкин A^., Медведев С.Ю., Мишагин K.r. Перспективный цифровой частотный компаратор/анализатор фазовых шумов // Вестник метролога. 2011. № 2. С. 18 - 20.
4. Игнатьев В.К, Никитин AB., Бернардо-Сапрыкин В.Х., Орлов A.A. Измерение разности фаз квазигармонических сигналов в реальном времени // Наука и образование. 2013. № 7. URL: http://dx.doi.org/10.7463/0713.0588392 (дата обращения: 21.11.2015).
5. Игнатьев В.К, Никитин AB., дочкин A^. Параметрический метод измерения разности фаз квазигармонических сигналов // Инженерный вестник Дона. 2013. № 3. URL: http://www.ivdon.ru/uploads/article/pdf/IVD_6_Kvochkin.pdf_1749.pdf (дата обращения 21.11.2015).
6. Игнатьев В.К, Станкевич ДА. Aппаратно-программный комплекс для параметрического анализа сигналов в задачах технической диагностики // Инженерный вестник Дона. 2013. № 3. URL: http://ivdon.ru/uploads/article/pdf/IVD_106_ignatjev.pdf_1843.pdf (дата обращения 21.11.2015).
7. Игнатьев В.К, Никитин A3., Юшанов С.В. Измерение фазового сдвига квазигармонических сигналов // Вычислительные методы и программирование. 2013. Т.14. С. 424 - 431.
8. Левин Б.Р. Теоретические основы статистической радиотехники. М.: Радио и связь. 1989. 656 с.
9. Kay S.M. Fundamentals of Statistical Signal Processing: Estimation Theory. NJ: Prentice Hall, Upper Sadle River, NJ, 1998. 595 p.
10. Сергиенко А.Б. Цифровая обработка сигналов. СПб.: Питер. 2002. 608 с.
References
1. Massimo, C. 6/12-channel synchronous digital phasemeter for ultrastable signal characterization and use. Frequency Control Symposium & the European Frequency and Time Forum (FCS), 2015 Joint Conference of the IEEE International.
2. Kawagoe, J., Kawasaki, T. A New Precision Digital Phase Meter and Its Simple Calibration Method. IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, VOL. 59, NO. 2, FEBRUARY 2010.
3. Dyukin A.B., Medvedev S.Yu., Mishagin K.G. Vestnik metrologa. 2011. № 2. pp. 18 - 20.
4. Ignat'ev V.K., Nikitin A.V., Bernardo-Saprykin V.Kh., Orlov A.A. Nauka i obrazovanie. 2013. № 7. URL: http://dx.doi.org/10.7463/0713.0588392 (data obrashcheniya: 21.11.2015).
5. Ignat'ev V.K., Nikitin A.V., Kvochkin A.I. Inzenernyj vestnik Dona (Rus). 2013. № 3. URL: ivdon.ru/uploads/article/pdf/IVD_6_Kvochkin.pdf_1749.pdf (data obrashcheniya 21.11.2015).
6. Ignat'ev V.K., Stankevich D.A. Inzenernyj vestnik Dona (Rus). 2013. № 3. URL:ivdon.ru/uploads/article/pdf/IVD_106_ignatjev.pdf_1843.pdf (data obrashcheniya 21.11.2015).
7. Ignat'ev V.K., Nikitin A.V., Yushanov S.V. Vychislitel'nye metody i programmirovanie. 2013. T.14. pp. 424 - 431.
8. Levin B.R. Teoreticheskie osnovy statisticheskoy radiotekhniki [Theoretical Foundations of Statistical Radio Engineering]. M.: Radio i svyaz'. 1989. 656 p.
9. Kay S.M. Fundamentals of Statistical Signal Processing: Estimation Theory. NJ: Prentice Hall, Upper Sadle River, NJ, 1998. 595 p.
J
10. Sergienko A.B. Tsifrovaya obrabotka signalov [Digital signal processing]. SPb: Piter. 2002. 608 p.
