RU2492498C2 - Способ определения направления на гидроакустический маяк-ответчик по горизонтальному и вертикальному углу - Google Patents

Abstract

Изобретение относится к области гидроакустики. Сущность: способ включает излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F_s, вычисляют взаимную корреляционную функцию сигналов x₁ и x₂ и для каждой пары преобразователей вычисляют разность моментов времени прихода плоского фронта волны к преобразователям. Затем для каждой плоскости вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z' соответствующей плоскости. Технический результат: уменьшение вычислительных операций для определения направления на гидроакустический маяк-ответсик, уменьшение количества преобразователей в антенных решетках при применении в качестве ответного сигнала маяка-ответчика сложного навигационного сигнала с фазовой манипуляцией или с линейной частотной модуляцией. 5 ил.

Description

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

Известен способ обнаружения шумящих в море объектов в фиксированном частотном диапазоне, при котором принимают шумовые сигналы в горизонтальной и вертикальной плоскостях многоэлементной антенной решеткой гидролокатора и осуществляют первичную обработку, для чего преобразуют в цифровую форму напряжения шумовых сигналов антенной решетки, выполняют преобразование Фурье отсчетов напряжений шумовых сигналов антенной решетки, вычисляют для каждого из полученных частотных отсчетов амплитудные и фазовые коэффициенты синфазного сложения напряжений сигналов антенной решетки, суммируют выходные напряжения сигналов антенной решетки с постоянными весами, равными произведению амплитудных и фазовых коэффициентов, чем образуют пространственные каналы наблюдения в горизонтальной и вертикальной плоскостях, квадрируют и осуществляют вторичную обработку на каждом цикле обзора, для чего суммируют по всем частотным отсчетам выходные напряжения образованных пространственных каналов в фиксированном частотном диапазоне, усредняют по времени, центрируют и нормируют шумовые сигналы к помехе, осуществляют наблюдение на каждом цикле обзора полученных отметок принятых шумовых сигналов и принимают решение об обнаружении путем сравнения с пороговым значением отношения сигнал-помеха, отличающийся тем, что на каждом цикле обзора до квадрирования для каждого частотного отсчета формируют адаптивные пространственные каналы наблюдения, каждый из которых образован, по крайней мере, тремя смежными пространственными каналами в горизонтальной или в вертикальной плоскости, для чего формируют взаимные спектры мощности между шумовыми сигналами пространственных каналов, участвующих в формировании адаптивного пространственного канала наблюдения, накапливают взаимные спектры мощности шумовых сигналов для заданного времени накопления, составляют матрицу накопленных взаимных спектров мощности шумовых сигналов и осуществляют ортогональное преобразование матрицы, вычисляют векторы фазовых коэффициентов синфазного сложения сигналов пространственных каналов, участвующих в формировании адаптивного пространственного канала наблюдения, вычисляют вектор выходных напряжений пространственных каналов, участвующих в формировании адаптивного пространственного канала путем решения векторно-матричного уравнения для ортогонально преобразованной матрицы накопленных взаимных спектров мощности и вектора фазовых коэффициентов синфазного сложения сигналов, квадрирование производят, вычисляя отклик полученного адаптивного пространственного канала наблюдения, равного обратной величине суммы квадратов элементов вектора выходных напряжений, а вторичную обработку осуществляют для выходных напряжений адаптивных пространственных каналов наблюдения (см. описание изобретения к патенту РФ №2339050, МПК G01S 3/80, G01S 15/04, публикация 20.11.2008).

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

Известен способ адаптивной пространственной фильтрации сигналов, заключающийся в приеме сигналов в N пространственно разнесенных точках, где N≥2, расчете весовых коэффициентов, взвешенном суммировании принятых сигналов, в котором генерируют первый и второй опорные сигналы с различными частотами, разделяют каждый из них на N равноамплитудных составляющих, сдвигают по фазе каждую составляющую на предварительно заданную величину фазы, после чего i-й принятый сигнал, где i=1, 2, … N, и i-e составляющие первого и второго опорных сигналов суммируют, причем i-й весовой коэффициент рассчитывают по формуле

$$W_i(t+\Delta t)=W_i(t)-\mu [V_{yi}-Y(t)Z^{*}{}_i(t)]$$

,

где W_i(t) - текущее значение i-го весового коэффициента;

W_yi - значение i-го управляющего сигнала;

Δt - интервал временной дискретизации;

Y(t) - взвешенно суммированный сигнал;

Z_i(t) - суммарный сигнал от сложения i-го принятого сигнала и i-x составляющих первого и второго опорных сигналов;

* - обозначение операции комплексного сопряжения;

µ - коэффициент усиления (передачи), а частоты первого ω₁ и второго ω₂ опорных сигналов выбирают из условия ω₁<(ω₀-Δω), ω₂>(ω₀+Δω), где ω₀ и Δω соответственно несущая частота и ширина спектра принимаемого полезного сигнала (см. описание изобретения к патенту РФ №2141706, МПК H01Q 21/00, публикация 20.11.1999).

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

Известен способ формирования диаграмм направленности N-элементной линейной эквидистантной антенной решетки, заключающийся в том, что сигналы с элементов антенной решетки усиливают, ограничивают в полосе частот и осуществляют преобразование Фурье каждого сигнала, осуществляют инвертирование преобразованного сигнала первого элемента и его попарное суммирование с преобразованными сигналами последующих элементов, для каждого из полученных (N-1) сигналов измеряют спектральную плотность мощности, умножают их на вещественные весовые коэффициенты A_n, которые определяют из соотношения

A_n=sin[nπ(N-1)/N]/n,

где n=1, (N-1),

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

$$F(\omega )={\displaystyle \sum _{n=1}^{N=i}{A}_n\left[1-\cos \left(n\omega \frac{d\sin \alpha }{c}\right)\right]}$$

,

где d - расстояние между соседними элементами;

c - скорость распространения волновых процессов в среде;

α - угол между направлением главного лепестка диаграммы направленности и нормалью к плоскости апертуры (см. описание изобретения к авторскому свидетельству СССР №1327026, МПК G01S 3/80, H01Q 3/26, публикация 30.07.1987).

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

Известен способ обзора пространства гидроакустической системой, включающий прием акустического сигнала в водной среде с помощью неподвижной гидроакустической антенны, преобразование акустического сигнала в электрические колебания, усиление и частотную фильтрацию этих колебаний, формирование N направлений пространственной фильтрации в пространстве обзора, сдвинутых относительно друг друга на угол $$\alpha =\frac{{\beta }^{\circ }}{N}$$

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

$$\frac{G_1(f)}{G_2(f)}=\frac{R^2(f,{\theta }_x)}{R^2(f,{\theta }_x+d)}$$

,

где G₁(f), G₂(f) энергетические спектры, соответствующие соседним направлениям пространственной фильтрации;

R(Θ_x) любая из характеристик направленности статического веера, причем Θ_x определяют для каждой из частот спектрального анализа в рабочем диапазоне,

после чего усредняют по ансамблю полученные данные, при этом после сопоставления энергетических спектров G₁(f) и G₂(f) подавляют сигналы в направлениях с меньшими значениями энергетических спектров на всех частотах (см. описание изобретения к патенту РФ №2047278, МПК G01S 7/52, H04R 1/44, публикация 27.10.1995).

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

Известен способ подледного приема сигналов спутниковых навигационных систем при нахождении подводного объекта на горизонте плавания с использованием гидроакустического канала передачи навигационной информации, включающий ввод в лед антенны приемника спутниковой навигационной системы в месте приледнения данной антенны к нижней кромке льда, прием этой антенной сигналов космических аппаратов, измерение приемником параметров этих сигналов и по полученным данным определение вычислительным путем искомых поправок к геодезическим счислимым координатам и курсу, вырабатываемых бортовой инерциальной навигационной системой подводного объекта, осуществление путем учета данных поправок коррекции бортовой инерциальной навигационной системы, в котором доставляют от горизонта плавания подводного объекта до места приледнения к нижней кромки льда антенну с приемником сигналов спутниковой навигационной системы и с механизмом внедрения ее в лед, двухканальный маяк-ответчик, который имеет разные коды запроса и разную частоту излучаемого гидроакустического сигнала, измеритель скорости распространения звука в воде и гидростатический датчик, внедряют антенну с приемником сигналов спутниковой навигационной системы в лед в месте ее приледнения к нижней кромке льда с помощью механизма внедрения до момента выхода данной антенны изо льда в воздушную среду, а на подводном объекте используют двухканальную гидроакустическую систему с двумя гидроакустическими приемоизлучающими антеннами, одну приемоизлучающую гидроакустическую антенну устанавливают в заданном месте в корме, а другую гидроакустическую приемоизлучающую антенну устанавливают в заданном месте в носу подводного объекта, при этом две данные антенны располагают вдоль направления его диаметральной плоскости или под известным углом к ней, определяют геодезические координаты маяка-ответчика приемником сигналов спутниковой навигационной системы, в период подъема до нижней кромки льда антенны с приемником сигналов спутниковой навигационной системы и их погружения до подводного объекта измеряют гидростатическое давление воды гидростатическим датчиком, а также скорость звука в воде измерителем скорости звука, а на подводном объекте синхронно излучают в направлении маяка-ответчика сигналы «запроса» одной и второй гидроакустическими приемоизлучающими антеннами и принимают переизлученные сигналы маяка-ответчика данными антеннами, по измеренным параметрам которых определяют вычислительным путем наклонные расстояния от подводного объекта до маяка-ответчика и пеленги на маяк-ответчик из двух заданных мест установки данных двух гидроакустических приемоизлучающих антенн двухканальной навигационной гидроакустической системы (см. описание изобретения к патенту РФ №2398316, МПК H01Q 1/04, публикация 27.08.2009).

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

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

Сущность заявляемого изобретения заключается в следующем.

Способ определения направления на гидроакустический маяк-ответчик по горизонтальному и вертикальному углу, включающий излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F_s, вычисляют взаимную корреляционную функцию сигналов x₁ и x₂, принятых преобразователями каждой пары по формуле:

$$r_{x_1{x}_2}\left(k+k_{max}\right)=\frac{1}{N-k_{max}}{\displaystyle \sum _{i=k_{max}+1}^{N-\left|k\right|}{x}_1\left(i\right)x_2\left(i-k\right)}$$

где k - аргумент взаимной корреляционной функции соответствующий сдвигу в отсчетах двух сигналов x₁ и x₂, изменяющийся в диапазоне - k_max≤k≤k_max;

$$k_{\max }=\frac{d}{c}\cdot F_s$$

- количество отсчетов максимально возможной с задержки сигнала на преобразователях одной пары;

N - количество отсчетов анализируемого фрагмента сигнала;

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

$$\tau =\left(k_{\max }+1-\arg \max \left(r_{x_1{x}_2}\right)\right)/Fs$$

,

где $$\arg \max \left(r_{x_1{x}_2}\right)$$

- позиция максимума взаимной корреляционной функции в отсчетах,

затем для каждой плоскости по найденным относительным задержкам для двух пар преобразователей и известным координатам акустических центров преобразователей вычисляют горизонтальный A′ и вертикальный E′ углы направления на источник сигнала в системе координат O′X′Y′Z′ соответствующей плоскости:

$$A^{\text{'}}=\pi +arctg\left(\frac{X_1{H}_2-X_2{H}_1}{Y_2{H}_1}\right)$$

,

$$E_1^{\text{'}}=\arccos \left(\frac{\sqrt{{\left(X_1{H}_2-X_2{H}_1\right)}^2+Y_2^2{H}_1^2}}{x_1{y}_2}\right)$$

,

$$E_2^{\text{'}}=2\pi -\arccos \left(\frac{\sqrt{{\left(X_1{H}_2-X_2{H}_1\right)}^2+Y_2^2{H}_1^2}}{X_1{Y}_2}\right)$$

,

где X₁=x₁-x₃ - разность абсцисс акустических центров преобразователей 1 и 3, образующих плоскость;

X₂=x₂-x₄ - разность абсцисс акустических центров преобразователей 2 и 4;

H₁=h₁-h₃, где h₁ и h₃ - расстояние, проходимое плоским фронтом волны от преобразователей 1 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;

H₂=h₂-h₃, где h₂ и h₃ - расстояния, проходимые плоским фронтом волны от преобразователей 2 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;

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

$${\nu }_{m1}^{\text{'}}=-\left[\begin{array}{c}\cos \left(E_{m1}\right)\cos \left(A\right)\\ \cos \left(E_{m1}\right)\sin \left(A\right)\\ \sin \left(E_{m1}\right)\end{array}\right]$$

и $${\nu }_{m2}^{\text{'}}=-\left[\begin{array}{c}\cos \left(E_{m2}\right)\cos \left(A\right)\\ \cos \left(E_{m2}\right)\sin \left(A\right)\\ \sin \left(E_{m2}\right)\end{array}\right]$$

,

где $${\nu }_{m1}^{\text{'}}$$

и $${\nu }_{m2}^{\text{'}}$$

- единичные векторы для плоскости m, направлений на маяк-ответчик (A, E₁) и (A, E₂) в системах координат плоскостей преобразователей $$O^{\text{'}}{X}^{\text{'}}{Y}^{\text{'}}{Z}_m^{\text{'}}$$

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

$$\nu \left(A,E\right)={\left[\begin{array}{c}{\cos }_{O^{\text{'}}{X}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Y}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Z}^{\text{'}}}{}^T\end{array}\right]}^T\cdot \nu \text{'}\left(A,E\right)$$

где cos_O′X′ - направляющий косинус оси O'X' в системе координат OXYZ;

cos_O′Y′ - направляющий косинус оси O'Y' в системе координат OXYZ;

cos_O′Z′ - направляющий косинус оси O'Z' в системе координат OXYZ;

после вычисления векторов $${\nu }_{m1}^{\text{'}}$$

и $${\nu }_{m2}^{\text{'}}$$

для всех m плоскостей находят истинный вектор направления $$\overline{\nu }$$

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

$$A=\arctan \left(\frac{y_{\overline{\nu }}}{x_{\overline{\nu }}}\right)$$

$$E=\arcsin \left(z_{\overline{\nu }}\right)$$

,

где $$x_{\overline{\nu }}$$

, $$y_{\overline{\nu }}$$

, $$z_{\overline{\nu }}$$

- прямоугольные координаты усредненного вектора направления на источник сигнала;

это позволяет определить значения горизонтального A и вертикального E углов, задающих направление на гидроакустический маяк-ответчик.

Сущность заявляемого изобретения поясняется чертежами: где на фиг.1 показана антенная решетка со своей системой координат OXYZ и фронт волны акустического сигнала маяка-ответчика:

S₁÷S₇, S_i - акустические центры преобразователей; α₀ - плоскость фронта волны принимаемого антенной сигнала в момент времени t=0, α_i - в момент времени t=t_i, когда фронт достигает преобразователя S_i; ν - нормаль к плоскости фронта волны, r_i - радиус-вектор преобразователя S_i, h_i - расстояние между фронтами волны в момент времени t=t_i и t=0; A и E - горизонтальный и вертикальный углы направления на источник сигнала - гидроакустический маяк-ответчик;

на фиг.2 показано разбиение антенной решетки на плоскости:

S₁÷S₇ - акустические центры преобразователей; ν - нормаль к плоскости фронта волны - вектор, задающий направление прихода сигнала; A и E - горизонтальный и вертикальный углы направления на источник сигнала - гидроакустический маяк-ответчик; OXYZ - система координат антенной решетки, O'X'Y'Z' - система координат, связанная с плоскостью преобразователей S₂S₃S₇.;

на фиг.3 показана плоскость со связанной системой координат O'X'Y'Z', образованная тремя преобразователями S₁, S₂, S₃, и ее расположение относительно фронта волны в момент времени t=0; t₁ и t₂ - моменты времени, когда фронт волны достигает преобразователи S₁ и S₂ соответственно, h₁ и h₂ - расстояния, измеренные относительно нулевого положения, проходимые фронтом волны при достижении преобразователей S₁ и S₂ соответственно;

на фиг.4 показана плоскость со связанной системой координат O'X'Y'Z', образованная двумя парами преобразователей S₁, S₃ и S₂, S₄, и расположение ее относительно фронта волны в момент времени t=0; t₁, t₂, t₃, t₄ - моменты времени, когда фронт волны достигают преобразователи S₁, S₂, S₃ и S₄ соответственно; h₁, h₂, h₃, h₄ - расстояния, измеренные относительно нулевого положения, проходимые фронтом волны при достижении преобразователей S₁, S₂, S₃ и S₄ соответственно;

на фиг.5 показаны координаты $$(x_1^{\text{'}},{\text{ y}}_{\text{1}}^{\text{'}})$$

, $$(x_2^{\text{'}},{\text{ y}}_{\text{2}}^{\text{'}})$$

, $$(x_3^{\text{'}},{\text{ y}}_{\text{3}}^{\text{'}})$$

, $$(x_4^{\text{'}},{\text{ y}}_{\text{4}}^{\text{'}})$$

преобразователей S₁, S₂, S₃ и S₄ в системе координат O′X′Y′, связанной с плоскостью преобразователей.

Заявляемый способ осуществляют следующим образом.

С пространственной антенной решеткой связывают трехмерную прямоугольную систему координат OXYZ, ось OZ которой направлена вниз для удобства, так как вертикальный угол под водой отсчитывают от плоскости OXY вниз. Радиус-вектор произвольного преобразователя S_i в системе координат OXYZ обозначим как r_i=[x_i y_i z_i]^T. Расположение преобразователей и системы координат показано на фиг.1.

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

Направление на источник сигнала задают горизонтальным и вертикальным углами - вектором [A, E] в сферической системе координат с началом координат в точке O. Горизонтальный угол A - угол между осью OX и проекцией вектора направления на плоскость OXY, отсчитывают от оси OX по часовой стрелке, если смотреть на плоскость OXY сверху. Вертикальный угол E - угол между вектором направления на источник и его проекцией на плоскость OXY и отсчитывают от плоскости OXY вниз (см. фиг.1).

Плоскости фронта волны α₀ и α_i, схематически показаны на фиг.1 эллипсами. Направление на источник сигнала есть нормаль к плоскости фронта волны - вектор ν на фиг.1. Благодаря пространственной удаленности преобразователей, фронт волны сигнала от источника излучения достигнет преобразователей в разные моменты времени. Для расчетов необходимы только относительные времена прихода сигнала к преобразователям, поэтому для удобства расчетов за t=0 принимают момент времени, когда плоский фронт волны достигнет точки O - начала системы координат антенной решетки. Фронт волны достигает i-й преобразователь в момент времени t=t_i, причем, если преобразователь находится ближе к источнику чем т.O, то знак у t_i будет отрицательный, а если дальше - положительный.

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

С каждой плоскостью связывают свою прямоугольную систему координат (на фиг.2 показана система координат O'X'Y'Z', связанная с плоскостью преобразователей S₂S₃S₇) и вычисляют горизонтальный и вертикальный углы направления на источник сигнала в системе координат соответствующей плоскости преобразователей. Для этого преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары. В случае трех преобразователей один из преобразователей будет общим для обоих пар. Для каждой пары преобразователей вычисляют разность моментов времени прихода фронта волны к преобразователям - относительная задержка. Для этого сигнал оцифровывают с частотой дискретизации F_s, вычисляют взаимную корреляционную функцию сигналов x₁ и x₂, принятых преобразователями каждой пары по формуле:

$$r_{x_1{x}_2}\left(k+k_{max}\right)=\frac{1}{N-k_{max}}{\displaystyle \sum _{i=k_{max}+1}^{N-\left|k\right|}{x}_1\left(i\right)x_2\left(i-k\right)}$$

где k - аргумент взаимной корреляционной функции соответствующий сдвигу в отсчетах двух сигналов x₁ и x₂, изменяющийся в диапазоне - k_max≤k≤k_max;

$$k_{\max }=\frac{d}{c}\cdot F_s$$

- количество отсчетов максимально возможной с задержки сигнала на преобразователях одной пары;

N - количество отсчетов анализируемого фрагмента сигнала;

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

$$\tau =\left(k_{\max }+1-\arg \max \left(r_{x_1{x}_2}\right)\right)/Fs$$

,

где $$\arg \max \left(r_{x_1{x}_2}\right)$$

- позиция максимума взаимной корреляционной функции в отсчетах,

Для каждой плоскости по найденным относительным задержкам для двух пар преобразователей и известным координатам акустических центров преобразователей вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z', связанной с данной плоскостью преобразователей.

Для этого зададим плоскость фронта волны ее нормалью ν' (см. фиг.1). В сферической системе координат вектор ν определяют тройкой чисел (|ν|, A', E'), где |ν|, - длина вектора. С учетом того, что вектор ν′ направлен от плоскости α_i (см. фиг.1) к точке O' (в соответствии с направлением прихода сигнала), его прямоугольные координаты определяют следующим образом:

$$x=-\left|\nu \text{'}\right|\cos (E)\cos (A)$$

,

$$y=-\left|\nu \text{'}\right|\cos (E)\sin (A)$$

,

$$z=-\left|\nu \text{'}\right|\sin (E)$$

.

Что можно записать в векторной форме:

$$\nu \text{'}(A\text{'},E\text{'})=-\left|\nu \text{'}\right|\left[\begin{array}{c}\cos \left(E\text{'}\right)\cos \left(A\text{'}\right)\\ \cos \left(E\text{'}\right)\sin \left(A\text{'}\right)\\ \sin \left(E\text{'}\right)\end{array}\right]$$

Так как нас интересует лишь направление вектора, то примем |ν|=1.

Тогда

$$\nu \text{'}(A\text{'},E\text{'})=-\left[\begin{array}{c}\cos \left(E\text{'}\right)\cos \left(A\text{'}\right)\\ \cos \left(E\text{'}\right)\sin \left(A\text{'}\right)\\ \sin \left(E\text{'}\right)\end{array}\right]$$

.

Рассмотрим фронт волны в двух положениях: в момент времени t=0, когда он находится в начале координат, и в момент времени t=t_i, когда он проходит преобразователь S_i (см. фиг.1). Обозначим расстояние между фронтами в эти моменты времени как h_i, которое определяют как скалярное произведение нормали к плоскости ν' и радиус-вектора r_i преобразователя S_i.

h_i=ν'(A', E')·r_i=-(x_icos(E')cos(A')+y_icos(E')sin(A')+z_isin(E'))

Сначала рассмотрим случай, когда две базы образованы тремя преобразователями S₁, S₂, S₃, т.е. один преобразователь общий, пусть это будет S₃. Свяжем с этими преобразователями систему координат O'X'Y'Z' таким образом, чтобы точка O' совпадала с S₃, ось О'Х' совпадала с вектором S₃S₁, а база S₃S₂ произвольно лежала в плоскости O'X'Y'Z'. Схему такого расположения преобразователей см. на фиг.3.

При анализе сигналов, приходящих на преобразователи возможно измерение только относительных задержек прихода сигнала на соответствующих преобразователях, в данном случае это S₃S₁ и S₃S₂. Зная скорость распространения волн, можно вычислить путь, который проходит фронт волны между преобразователями, соответственно Δh₁₃ и Δh₂₃.

Составим систему уравнений для трех преобразователей:

x₁cos(E')cos(A')+y₁cos(E')sin(A')+z₁sin(E')=-h₁

x₂cos(E')cos(A')+y₂cos(E')sin(A')+z₂sin(E')=-h₂

x₃cos(E')cos(A')+y₃cos(E')sin(A')+z₃sin(E')=-h₃

Δh₁₃=h₁-h₃

Δh₂₃=h₂-h₃

Так как за точку, в которой t=0, принято начало координат, в котором находится преобразователь S₃, то h₃=0, x₃=y₃=z₃=0. Координата по оси O'Z' у всех преобразователей равна 0, у S₁ также равна нулю координата по оси O'Y'.

Учитывая все вышеизложенные условия, систему упрощают до следующего вида:

x₁cos(E')cos(A')=-Δh₁₃

x₂cos(E')cos(A')+y₂cos(E')sin(A')=-Δh₂₃

Для простоты обозначим Δh₃₁=H₁, Δh₃₂=H₂ и окончательно получим:

x₁cos(E')cos(A')=-H₁

x₂cos(E')cos(A')+y₂cos(E')sin(A')=-H₂

Рассмотрим теперь случай, когда плоскость образована пересечением прямых, соответствующих двум базам преобразователей S₄S₂ и S₃S₁ (см. фиг.4). Свяжем с преобразователями систему координат O'X'Y'Z' таким образом, чтобы точка O' лежала на пересечении прямых S₄S₂ и S₃S₁, ось О'Х' совпадала с вектором S₃S₁, а база S₄S₂ произвольно лежала в плоскости O'X'Y'Z' (см. фиг.4). По относительным задержкам сигналов можно найти Δh₁₃ и Δh₂₄.

Составим систему уравнений для данного случая:

x₁cos(E')cos(A')+y₁cos(E')sin(A')+z₁sin(E')=-h₁

x₂cos(E')cos(A')+y₂cos(E')sin(A')+z₂sin(E')=-h₂

x₃cos(E')cos(A')+y₃cos(E')sin(A')+z₃sin(E')=-h₃

x₄cos(E')cos(A')+y₄cos(E')sin(A')+z₄sin(E')=-h₄

Δh₁₃=h₁-h₃

Δh₂₄=h₂-h₄

Учитывая, что координата по оси O'Z' у всех преобразователей равна 0, у S₁ и S₃ также равны нулю координаты по оси O'Y', упрощаем систему:

x₁cos(E')cos(A')=-h₁

x₂cos(E')cos(A')+y₂cos(E')sin(A')=-h₂

x₃cos(E′)cos(A′)=-h₃

x₄cos(E′)cos(A′)+y₄cos(E′)sin(A′)=-h₄

Δh₁₃=h₁-h₃

Δh₂₄=h₂-h₄

Вычитая третье уравнение из первого, а четвертое из второго, и подставляя в них Δh₁₃, Δh₂₄, получаем:

(x₁-x₃)cos(E′)cos(A′)=-Δh₁₃

(x₂-x₄)cos(E′)cos(A′)+(y₂-y₄)cos(E′)sin(A′)=-Δh₂₄

Произведем замену: x₁-x₃=X₁, x₂-x₄=X₂, Δh₃₁=H₁, Δh₃₂=H₂, в результате которой получим окончательно:

x₁cos(E′)cos(A′)=-H₁

x₂cos(E′)cos(A′)+Y₂cos(E′)sin(A′)=-H₂

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

Полученная система тригонометрических уравнений дает два решения относительно A' и E':

$$\{\begin{array}{l}A\text{'}=\pi +arctg(\frac{x_1{H}_2-x_2{H}_1}{y_2{H}_1})\\ E_1^{\text{'}}=arccos(\frac{\sqrt{{\left(x_1{H}_2-x_2{H}_1\right)}^2+y_2^2{H}_1^2}}{x_1{y}_2})\end{array}$$

$$\{\begin{array}{l}A\text{'}=\pi +arctg(\frac{x_1{H}_2-x_2{H}_1}{y_2{H}_1})\\ E_2^{\text{'}}=2\pi -arccos(\frac{\sqrt{{\left(x_1{H}_2-x_2{H}_1\right)}^2+y_2^2{H}_1^2}}{x_1{y}_2})\end{array}$$

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

Для того чтобы по направлению, полученному в системе координат плоскости преобразователей, определить вектор направления в системе координат антенны, необходимо выполнить поворот вектора направления в соответствии с направляющими косинусами осей системы координат O'X'Y'Z' в системе OXYZ. Вектор направления на источник получим по горизонтальному A' и вертикальному E' углам по формуле:

$$\nu \text{'}\left(A\text{'},E\text{'}\right)=-\left[\begin{array}{c}\cos \left(E\text{'}\right)\cos \left(A\text{'}\right)\\ \cos \left(E\text{'}\right)\sin \left(A\text{'}\right)\\ \sin \left(E\text{'}\right)\end{array}\right]$$

Направляющие косинусы системы координат O'X'Y'Z' получим, исходя из вышеперечисленных условий привязки этой системы координат к двум парам преобразователей. Рассмотрим теперь сразу случай четырех преобразователей в одной плоскости, так как случай трех преобразователей аналогичен. Система координат и преобразователи изображены на фиг.5. Для начала нам необходимо определить точку пересечения прямых, образованных базами преобразователей. Это можно сделать, решив систему уравнений для прямых S₄S₂ и S₃S₁ (см. фиг.5):

$$\frac{x-x_1}{x_3-x_1}=\frac{y-y_1}{y_3-y_1}=\frac{z-z_1}{z_3-z_1}$$

$$\frac{x-x_2}{x_4-x_2}=\frac{y-y_2}{y_4-y_2}=\frac{z-z_2}{z_4-z_2}$$

Таким образом, получаем координаты (x_o, y_o, z_o) точки O' - начало координат для координатной системы O'X'Y'Z'.

Далее по двум точкам из множества S₄S₂S₃S₁ и точке O' нужно определить направляющие косинусы осей координат системы координат O'X'Y'Z'. Мы условились, что одна из прямых, соединяющих преобразователи, будет осью O'X'. Пусть ось O'X' проходит через точку S₁, как и в предыдущих случаях.

Определим длины сторон треугольника O'S₁S₂:

$$\left|O\text{'}{S}_1\right|=\sqrt{{\left(x_1-x_0\right)}^2+{\left(y_1-y_0\right)}^2+{\left(z_1-z_0\right)}^2}$$

$$\left|O\text{'}{S}_2\right|=\sqrt{{\left(x_2-x_0\right)}^2+{\left(y_2-y_0\right)}^2+{\left(z_2-z_0\right)}^2}$$

$$\left|S_1{S}_2\right|=\sqrt{{\left(x_2-x_1\right)}^2+{\left(y_2-y_1\right)}^2+{\left(z_2-z_1\right)}^2}$$

Вектор направляющих косинусов оси О′X′ определяют просто:

$${\cos }_{O^{\text{'}}{X}^{\text{'}}}=\left[\begin{array}{c}\left(x_1-x_0\right)/\left|O\text{'}{S}_1\right|\\ \left(y_1-y_0\right)/\left|O\text{'}{S}_1\right|\\ \left(z_1-z_0\right)/\left|O\text{'}{S}_1\right|\end{array}\right]=\left(S_1-O\text{'}\right)/\left|O\text{'}{S}_1\right|$$

Определим направляющие косинусы оси O'Y':

$${\cos }_{O^{\text{'}}{Y}^{\text{'}}}=\left[\begin{array}{c}\left(x_R-x_0\right)/\left|O\text{'}R\right|\\ \left(y_R-y_0\right)/\left|O\text{'}R\right|\\ \left(z_R-z_0\right)/\left|O\text{'}R\right|\end{array}\right]=\left(R-O\text{'}\right)/\left|O\text{'}R\right|$$

Для этого необходимо определить координаты точки R в системе координат OXYZ. Точка R является проекцией S₂ на ось O'Y'.

Прямая RS₂ параллельна оси O'X', а значит, их направляющие косинусы совпадают.

Вектор направляющих косинусов прямой RS₂ вычисляют как:

$${\cos }_{R{S}_2}=\left[\begin{array}{c}\left(x_2-x_R\right)/\left|R{S}_2\right|\\ \left(y_2-y_R\right)/\left|R{S}_2\right|\\ \left(z_2-z_R\right)/\left|R{S}_2\right|\end{array}\right]=\left(S_2-R\right)/\left|R{S}_2\right|$$

Поскольку прямая ОХ параллельна RS₂, то cos_OX=cos_RS2, соответственно:

(S₂-R)/|RS₂|=cos_O′X′.

Откуда получаем координатный вектор точки R:

R=S₂-|RS₂|=cos_O′X′.

Отрезок |O'R| определяем из прямоугольного треугольника O′PS₂:

$$\left|O\text{'}R\right|=\sqrt{{\left|O\text{'}{S}_2\right|}^2-{\left|O\text{'}P\right|}^2}$$

.

Отрезок |O'P| можно определить следующим образом:

|O'P|=|O′S₂|cos(α).

По теореме косинусов:

$$\cos \left(\alpha \right)\frac{{\left|O\text{'}{S}_1\right|}^2+{\left|O\text{'}{S}_2\right|}^2-{\left|S_1{S}_2\right|}^2}{2\left|O\text{'}{S}_1\right|\left|O\text{'}{S}_2\right|}$$

Подставляя данное выражение в предыдущее, получаем:

$$\left|O\text{'}P\right|=\left|O\text{'}{S}_2\right|\frac{{\left|O\text{'}{S}_1\right|}^2+{\left|O\text{'}{S}_2\right|}^2-{\left|S_1{S}_2\right|}^2}{2\left|O\text{'}{S}_1\right|\left|O\text{'}{S}_2\right|}$$

Упрощая, получаем:

$$\left|O\text{'}P\right|=\frac{{\left|O\text{'}{S}_1\right|}^2+{\left|O\text{'}{S}_2\right|}^2-{\left|S_1{S}_2\right|}^2}{2\left|O\text{'}{S}_1\right|}$$

Подставляя полученное выражение для |O'P| в выражение для вычисления значения |O′R|, окончательно получаем длину отрезка O'R:

$$\left|O\text{'}R\right|=\sqrt{{\left|O\text{'}{S}_2\right|}^2-\frac{{\left({\left|O\text{'}{S}_1\right|}^2+{\left|O\text{'}{S}_2\right|}^2-{\left|S_1{S}_2\right|}^2\right)}^2}{4{\left|O\text{'}{S}_1\right|}^2}}$$

Итак, зная значение координат точки R и длину отрезка O'R, получим направляющие косинусы оси O'Y'.

Направляющие косинусы оси O'Z' определяют, исходя из того, что эта ось перпендикулярна к плоскости O'X'Y', а значит, ее направляющие косинусы равны направляющим косинусам нормали к этой плоскости.

Зная уравнение плоскости:

A·x+B·y+C·z+D=0,

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

$${\cos }_{O^{\text{'}}{Z}^{\text{'}}}=\left[\begin{array}{c}\frac{A}{A^2+B^2+C^2}\\ \frac{B}{A^2+B^2+C^2}\\ \frac{C}{A^2+B^2+C^2}\end{array}\right]$$

В нашем случае, проще задать уравнение плоскости по 3-м точкам:

O'(x₀, y₀, z₀), S₁(x₁, y₁, z₁), S₂(x₂, y₂, z₂), тогда A, B, C задают следующими определителями:

$$A=\left|\begin{array}{ccc}{y}_0& z_0& 1\\ y_1& z_1& 1\\ y_2& z_2& 1\end{array}\right|$$

$$B=\left|\begin{array}{ccc}{z}_0& x_0& 1\\ z_1& x_1& 1\\ z_2& x_2& 1\end{array}\right|$$

$$C=\left|\begin{array}{ccc}{x}_0& y_0& 1\\ x_1& y_1& 1\\ x_2& y_2& 1\end{array}\right|$$

Зная направляющие косинусы всех осей системы координат O'X'Y'Z', можно определить матрицу вращения T для перехода из координатной системы O'X'Y'Z' в OXYZ.

$$T=\left[\begin{array}{c}{\cos }_{O^{\text{'}}{X}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Y}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Z}^{\text{'}}}{}^T\end{array}\right]$$

Теперь, используя данную матрицу вращения, мы можем перевести вектор направления на источник сигнала ν' из координатной системы O'X'Y'Z' в OXYZ следующим образом (в матричном виде):

ν(A, E)=T^T·ν(A, E)'.

После вычисления векторов $${\nu }_{m1}^{\text{'}}$$

и $${\nu }_{m2}^{\text{'}}$$

для всех m плоскостей антенной решетки находят истинный вектор направления $$\overline{\nu }$$

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

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

$$A=\arctan \left(\frac{y_{\overline{\nu }}}{x_{\overline{\nu }}}\right)$$

,

$$E=\arcsin \left(z_{\overline{\nu }}\right)$$

,

где $$x_{\overline{\nu }}$$

, $$y_{\overline{\nu }}$$

, $$z_{\overline{\nu }}$$

- прямоугольные координаты усредненного вектора направления на источник сигнала,

в результате чего получают значения горизонтального A и вертикального E углов, задающих направление на гидроакустический маяк-ответчик.

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

Claims (1)

1. Способ определения направления на гидроакустический маяк-ответчик по горизонтальному и вертикальному углу, включающий излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, отличающийся тем, что пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F_s, вычисляют взаимную корреляционную функцию сигналов x₁ и x₂, принятых преобразователями каждой пары по формуле:
   $$r_{x_1{x}_2}\left(k+k_{max}\right)=\frac{1}{N-k_{max}}{\displaystyle \sum _{i=k_{max}+1}^{N-\left|k\right|}{x}_1\left(i\right)x_2\left(i-k\right)},$$

   

   
   где k - аргумент взаимной корреляционной функции, соответствующий сдвигу в отсчетах двух сигналов x₁ и x₂, изменяющийся в диапазоне - k_max≤k≤k_max;
   $$k_{\max }=\frac{d}{c}\cdot F_s$$

   

   - количество отсчетов максимально возможной задержки с сигнала на преобразователях одной пары;
   N - количество отсчетов анализируемого фрагмента сигнала;
   и для каждой пары преобразователей вычисляют разность моментов времени прихода плоского фронта волны к преобразователям:
   $$\tau =\left(k_{\max }+1-\arg \max \left(r_{x_1{x}_2}\right)\right)/Fs$$

   

   ,
   где $$\arg \max \left(r_{x_1{x}_2}\right)$$

   

   - позиция максимума взаимной корреляционной функции в отсчетах,
   затем для каждой плоскости по найденным относительным задержкам для двух пар преобразователей и известным координатам акустических центров преобразователей вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z' соответствующей плоскости:
   $$A^{\text{'}}=\pi +arctg\left(\frac{X_1{H}_2-X_2{H}_1}{Y_2{H}_1}\right),$$

   

   
   $$E_1^{\text{'}}=\arccos \left(\frac{\sqrt{{\left(X_1{H}_2-X_2{H}_1\right)}^2+Y_2^2{H}_1^2}}{x_1{y}_2}\right),$$

   

   
   $$E_2^{\text{'}}=2\pi -\arccos \left(\frac{\sqrt{{\left(X_1{H}_2-X_2{H}_1\right)}^2+Y_2^2{H}_1^2}}{X_1{Y}_2}\right),$$

   

   
   где X₁=x₁-x₃ - разность абсцисс акустических центров преобразователей 1 и 3, образующих плоскость;
   Х₂=x₂-x₄ - разность абсцисс акустических центров преобразователей 2 и 4;
   H₁=h₁-h₃, где h₁ и h₃ - расстояние, проходимое плоским фронтом волны от преобразователей 1 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;
   H₂=h₂-h₃, где h₂ и h₃ - расстояния, проходимые плоским фронтом волны от преобразователей 2 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;
   затем для каждого направления, задаваемого значениями горизонтального и вертикального углов, вычисляют единичный вектор в прямоугольной системе координат, связанной с каждой из плоскостей:
   $${\nu }_{m1}^{\text{'}}=-\left[\begin{array}{c}\cos \left(E_{m1}\right)\cos \left(A\right)\\ \cos \left(E_{m1}\right)\sin \left(A\right)\\ \sin \left(E_{m1}\right)\end{array}\right]$$

   

   и $${\nu }_{m2}^{\text{'}}=-\left[\begin{array}{c}\cos \left(E_{m2}\right)\cos \left(A\right)\\ \cos \left(E_{m2}\right)\sin \left(A\right)\\ \sin \left(E_{m2}\right)\end{array}\right],$$

   

   
   где $${\nu }_{m1}^{\text{'}}$$

   

   и $${\nu }_{m2}^{\text{'}}$$

   

   - единичные векторы для плоскости m, направлений на маяк-ответчик (A, E₁) и (A, E₂) в системах координат плоскостей преобразователей $$O^{\text{'}}{X}^{\text{'}}{Y}^{\text{'}}{Z}_m^{\text{'}}$$

   

   ,
   и выполняют перевод прямоугольных координат найденных векторов в систему координат, связанную с антенной решеткой, по формуле:
   $$\nu \left(A,E\right)={\left[\begin{array}{c}{\cos }_{O^{\text{'}}{X}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Y}^{\text{'}}}{}^T\\ {\cos }_{O^{\text{'}}{Z}^{\text{'}}}{}^T\end{array}\right]}^T\cdot \nu \text{'}\left(A,E\right),$$

   

   
   где cos_O'X' - направляющий косинус оси O'X' в системе координат OXYZ;
   cos_O'Y' - направляющий косинус оси O'Y' в системе координат OXYZ;
   cos_O'Z' - направляющий косинус оси O'Z' в системе координат OXYZ;
   после вычисления векторов $${\nu }_{m1}^{\text{'}}$$

   

   и $${\nu }_{m2}^{\text{'}}$$

   

   для всех m плоскостей находят истинный вектор направления $$\overline{\nu }$$

   

   на источник сигнала путем усреднения тех векторов, расстояния между которыми, вычисленные по Евклидовой норме, минимальны, затем координаты найденного таким образом вектора переводят в сферическую систему координат, связанную с антенной решеткой, с помощью формул:
   $$A=\arctan \left(\frac{y_{\overline{\nu }}}{x_{\overline{\nu }}}\right),$$

   

   
   $$E=\arcsin \left(z_{\overline{\nu }}\right)$$

   

   ;
   где $$x_{\overline{\nu }}$$

   

   , $$y_{\overline{\nu }}$$

   

   , $$z_{\overline{\nu }}$$

   

   - прямоугольные координаты усредненного вектора направления на источник сигнала;
   в результате чего получают значения горизонтального A и вертикального E углов, задающих направление на гидроакустический маяк-ответчик.

Images