Модификация метода Прони при приеме сигналов векторно-скалярной антенной

Размещено на http://www.allbest.ru/

Ростовский государственный строительный университет,

Научно-исследовательский институт физики Южного федерального университета,

Модификация метода Прони при приеме сигналов векторно-скалярной антенной

Т.Н. Ларина

Г.М. Глебова

Е.В. Винник

г. Ростов-на-Дону

Хорошо известный параметрический метод Прони является методом восстановления квазиполинома по конечному числу его значений на равномерной сетке временных или пространственных отсчетов [1-3]. Метод Прони использует представление наблюдаемого процесса в виде комплексного экспоненциального ряда. Метод позволяет по отсчетам сигнала найти параметры этих комплексных экспонент, что, в свою очередь, дает возможность записать выражение для спектральной или пространственной плотности исследуемого сигнала. Широкое применение метода Прони стало возможным только в последнее время, поскольку он существенно не линеен и требует больших вычислительных затрат. В прикладных задачах гидроакустики метод применяется как для оценки спектра принимаемых сигналов, так и для определения угловых координат сигналов от локальных источников. Модификация метода для определения параметров коррелированных сигналов или нормальных волн, образующих акустическое поле источника в волноводе предложена в работе [4].

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

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

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

,

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

,,

здесь - мощность сигнала на приемнике давления, - волновое число, - - координаты источника в полярной системе координат, - пеленг источника отсчитывается от оси Х. В работе [5] теоретически и экспериментально показано, что представление сигналов в виде (3) справедливо как при распространении сигналов в свободном пространстве, так и в волноводе.

Для гауссовых сигналов и шумов с нулевым математическим ожиданием статистика измерений полностью определяется матрицей ковариаций, которая рассчитывается для заданной модели сигналов и помех как K=U•U* , символ «*» означает эрмитово сопряжение. Матрица К размером имеет блочно-диагональный вид:

,

и представляет собой сумму сигнальной матрицы - и матрицы помех -

.

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

Метод Прони является «быстрым» методом решения системы уравнений следующего вида:

,

где -- неизвестные комплексные величины (2•L<M). Неизвестные находятся как корни полинома

,

коэффициенты которого удовлетворяют системе линейных уравнений

После нахождения величин ; значения их подставляют в (6) и решают полученную линейную систему уравнений относительно.

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

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

.

Введем обозначения:

и

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

.

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

.

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

Литература

1. Prony G.R.B. Essai experemenal et analytique: sur les lois de la dilitabilite de fluids elastques et sur celles de la force expanslve de la vapeure de l'eau et la vapeure de l'alkool, a differentes temperatures. //J. de L'Ecole Polytechnique. -1795. - T.1. -24-76.

2. Марпл С.П. Цифровой спектральный анализ и его приложения. -М. -Мир. -1990.

3. Backer H.P. Cjmparison of FFT and Prony algorithm for bearing estimation of narrow-band signals in a realistic ocean environment. // - JASA. - Mar. - 1977. - V.61. - P. 756-762.

4. Гительсон В.С., Глебова Г.М., Кузнецов Г.Н. Определение параметров коррелированных сигналов с использованием метода Прони. // Фкустический журнал. - 1988. - Т.XXXIV. - 1. - 170-172/

5. Гордиенко В.А. Векторно-фазовые методы в акустике. -М.: -ФИЗМАТЛИТ, -2007. - 480 с.

Размещено на Allbest.ru