×
10.09.2013
216.012.68e3

СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ

Вид РИД

Изобретение

Юридическая информация Свернуть Развернуть
Краткое описание РИД Свернуть Развернуть
Аннотация: Изобретение относится к области гидроакустики. Сущность: способ включает излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F, вычисляют взаимную корреляционную функцию сигналов x и x и для каждой пары преобразователей вычисляют разность моментов времени прихода плоского фронта волны к преобразователям. Затем для каждой плоскости вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z' соответствующей плоскости. Технический результат: уменьшение вычислительных операций для определения направления на гидроакустический маяк-ответсик, уменьшение количества преобразователей в антенных решетках при применении в качестве ответного сигнала маяка-ответчика сложного навигационного сигнала с фазовой манипуляцией или с линейной частотной модуляцией. 5 ил.
Основные результаты: Способ определения направления на гидроакустический маяк-ответчик по горизонтальному и вертикальному углу, включающий излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, отличающийся тем, что пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F, вычисляют взаимную корреляционную функцию сигналов x и x, принятых преобразователями каждой пары по формуле: где k - аргумент взаимной корреляционной функции, соответствующий сдвигу в отсчетах двух сигналов x и x, изменяющийся в диапазоне - k≤k≤k; - количество отсчетов максимально возможной задержки с сигнала на преобразователях одной пары;N - количество отсчетов анализируемого фрагмента сигнала;и для каждой пары преобразователей вычисляют разность моментов времени прихода плоского фронта волны к преобразователям: ,где - позиция максимума взаимной корреляционной функции в отсчетах,затем для каждой плоскости по найденным относительным задержкам для двух пар преобразователей и известным координатам акустических центров преобразователей вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z' соответствующей плоскости: где X=x-x - разность абсцисс акустических центров преобразователей 1 и 3, образующих плоскость;Х=x-x - разность абсцисс акустических центров преобразователей 2 и 4;H=h-h, где h и h - расстояние, проходимое плоским фронтом волны от преобразователей 1 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;H=h-h, где h и h - расстояния, проходимые плоским фронтом волны от преобразователей 2 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;затем для каждого направления, задаваемого значениями горизонтального и вертикального углов, вычисляют единичный вектор в прямоугольной системе координат, связанной с каждой из плоскостей: и где и - единичные векторы для плоскости m, направлений на маяк-ответчик (A, E) и (A, E) в системах координат плоскостей преобразователей ,и выполняют перевод прямоугольных координат найденных векторов в систему координат, связанную с антенной решеткой, по формуле: где cos - направляющий косинус оси O'X' в системе координат OXYZ;cos - направляющий косинус оси O'Y' в системе координат OXYZ;cos - направляющий косинус оси O'Z' в системе координат OXYZ;после вычисления векторов и для всех m плоскостей находят истинный вектор направления на источник сигнала путем усреднения тех векторов, расстояния между которыми, вычисленные по Евклидовой норме, минимальны, затем координаты найденного таким образом вектора переводят в сферическую систему координат, связанную с антенной решеткой, с помощью формул: ;где , , - прямоугольные координаты усредненного вектора направления на источник сигнала;в результате чего получают значения горизонтального A и вертикального E углов, задающих направление на гидроакустический маяк-ответчик.
Реферат Свернуть Развернуть

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

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

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

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

,

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

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

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

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

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

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

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

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

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

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

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

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

,

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

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

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

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

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

,

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

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

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

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

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

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

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

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

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

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

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

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

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

,

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

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

,

,

,

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

X2=x2-x4 - разность абсцисс акустических центров преобразователей 2 и 4;

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

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

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

и ,

где и - единичные векторы для плоскости m, направлений на маяк-ответчик (A, E1) и (A, E2) в системах координат плоскостей преобразователей

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

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

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

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

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

,

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

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

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

S1÷S7, Si - акустические центры преобразователей; α0 - плоскость фронта волны принимаемого антенной сигнала в момент времени t=0, αi - в момент времени t=ti, когда фронт достигает преобразователя Si; ν - нормаль к плоскости фронта волны, ri - радиус-вектор преобразователя Si, hi - расстояние между фронтами волны в момент времени t=ti и t=0; A и E - горизонтальный и вертикальный углы направления на источник сигнала - гидроакустический маяк-ответчик;

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

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

на фиг.3 показана плоскость со связанной системой координат O'X'Y'Z', образованная тремя преобразователями S1, S2, S3, и ее расположение относительно фронта волны в момент времени t=0; t1 и t2 - моменты времени, когда фронт волны достигает преобразователи S1 и S2 соответственно, h1 и h2 - расстояния, измеренные относительно нулевого положения, проходимые фронтом волны при достижении преобразователей S1 и S2 соответственно;

на фиг.4 показана плоскость со связанной системой координат O'X'Y'Z', образованная двумя парами преобразователей S1, S3 и S2, S4, и расположение ее относительно фронта волны в момент времени t=0; t1, t2, t3, t4 - моменты времени, когда фронт волны достигают преобразователи S1, S2, S3 и S4 соответственно; h1, h2, h3, h4 - расстояния, измеренные относительно нулевого положения, проходимые фронтом волны при достижении преобразователей S1, S2, S3 и S4 соответственно;

на фиг.5 показаны координаты , , , преобразователей S1, S2, S3 и S4 в системе координат O′X′Y′, связанной с плоскостью преобразователей.

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

С пространственной антенной решеткой связывают трехмерную прямоугольную систему координат OXYZ, ось OZ которой направлена вниз для удобства, так как вертикальный угол под водой отсчитывают от плоскости OXY вниз. Радиус-вектор произвольного преобразователя Si в системе координат OXYZ обозначим как ri=[xi yi zi]T. Расположение преобразователей и системы координат показано на фиг.1.

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

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

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

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

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

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

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

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

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

,

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

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

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

,

,

.

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

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

Тогда

.

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

hi=ν'(A', E')·ri=-(xicos(E')cos(A')+yicos(E')sin(A')+zisin(E'))

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

При анализе сигналов, приходящих на преобразователи возможно измерение только относительных задержек прихода сигнала на соответствующих преобразователях, в данном случае это S3S1 и S3S2. Зная скорость распространения волн, можно вычислить путь, который проходит фронт волны между преобразователями, соответственно Δh13 и Δh23.

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

x1cos(E')cos(A')+y1cos(E')sin(A')+z1sin(E')=-h1

x2cos(E')cos(A')+y2cos(E')sin(A')+z2sin(E')=-h2

x3cos(E')cos(A')+y3cos(E')sin(A')+z3sin(E')=-h3

Δh13=h1-h3

Δh23=h2-h3

Так как за точку, в которой t=0, принято начало координат, в котором находится преобразователь S3, то h3=0, x3=y3=z3=0. Координата по оси O'Z' у всех преобразователей равна 0, у S1 также равна нулю координата по оси O'Y'.

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

x1cos(E')cos(A')=-Δh13

x2cos(E')cos(A')+y2cos(E')sin(A')=-Δh23

Для простоты обозначим Δh31=H1, Δh32=H2 и окончательно получим:

x1cos(E')cos(A')=-H1

x2cos(E')cos(A')+y2cos(E')sin(A')=-H2

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

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

x1cos(E')cos(A')+y1cos(E')sin(A')+z1sin(E')=-h1

x2cos(E')cos(A')+y2cos(E')sin(A')+z2sin(E')=-h2

x3cos(E')cos(A')+y3cos(E')sin(A')+z3sin(E')=-h3

x4cos(E')cos(A')+y4cos(E')sin(A')+z4sin(E')=-h4

Δh13=h1-h3

Δh24=h2-h4

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

x1cos(E')cos(A')=-h1

x2cos(E')cos(A')+y2cos(E')sin(A')=-h2

x3cos(E′)cos(A′)=-h3

x4cos(E′)cos(A′)+y4cos(E′)sin(A′)=-h4

Δh13=h1-h3

Δh24=h2-h4

Вычитая третье уравнение из первого, а четвертое из второго, и подставляя в них Δh13, Δh24, получаем:

(x1-x3)cos(E′)cos(A′)=-Δh13

(x2-x4)cos(E′)cos(A′)+(y2-y4)cos(E′)sin(A′)=-Δh24

Произведем замену: x1-x3=X1, x2-x4=X2, Δh31=H1, Δh32=H2, в результате которой получим окончательно:

x1cos(E′)cos(A′)=-H1

x2cos(E′)cos(A′)+Y2cos(E′)sin(A′)=-H2

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

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

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

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

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

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

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

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

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

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

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

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

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

Поскольку прямая ОХ параллельна RS2, то cosOX=cosRS2, соответственно:

(S2-R)/|RS2|=cosO′X′.

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

R=S2-|RS2|=cosO′X′.

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

.

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

|O'P|=|O′S2|cos(α).

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

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

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

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

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

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

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

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

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

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

O'(x0, y0, z0), S1(x1, y1, z1), S2(x2, y2, z2), тогда A, B, C задают следующими определителями:

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

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

ν(A, E)=TT·ν(A, E)'.

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

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

,

,

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

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

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

Способ определения направления на гидроакустический маяк-ответчик по горизонтальному и вертикальному углу, включающий излучение судном сигнала запроса и прием пространственной антенной решеткой сигнала ответа от маяка-ответчика, отличающийся тем, что пространственную антенную решетку разбивают на набор плоскостей таким образом, чтобы каждая плоскость образовывалась акустическими центрами трех или четырех преобразователей, лежащих на двух взаимно пересекающихся прямых, с каждой плоскостью связывают свою прямоугольную систему координат и вычисляют горизонтальный и вертикальный углы направления на источник сигнала, для чего преобразователи, акустические центры которых образуют данную плоскость, разбивают на две пары, сигнал оцифровывают с частотой дискретизации F, вычисляют взаимную корреляционную функцию сигналов x и x, принятых преобразователями каждой пары по формуле: где k - аргумент взаимной корреляционной функции, соответствующий сдвигу в отсчетах двух сигналов x и x, изменяющийся в диапазоне - k≤k≤k; - количество отсчетов максимально возможной задержки с сигнала на преобразователях одной пары;N - количество отсчетов анализируемого фрагмента сигнала;и для каждой пары преобразователей вычисляют разность моментов времени прихода плоского фронта волны к преобразователям: ,где - позиция максимума взаимной корреляционной функции в отсчетах,затем для каждой плоскости по найденным относительным задержкам для двух пар преобразователей и известным координатам акустических центров преобразователей вычисляют горизонтальный A' и вертикальный E' углы направления на источник сигнала в системе координат O'X'Y'Z' соответствующей плоскости: где X=x-x - разность абсцисс акустических центров преобразователей 1 и 3, образующих плоскость;Х=x-x - разность абсцисс акустических центров преобразователей 2 и 4;H=h-h, где h и h - расстояние, проходимое плоским фронтом волны от преобразователей 1 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;H=h-h, где h и h - расстояния, проходимые плоским фронтом волны от преобразователей 2 и 3 до начала системы координат, связанной с рассматриваемой плоскостью;затем для каждого направления, задаваемого значениями горизонтального и вертикального углов, вычисляют единичный вектор в прямоугольной системе координат, связанной с каждой из плоскостей: и где и - единичные векторы для плоскости m, направлений на маяк-ответчик (A, E) и (A, E) в системах координат плоскостей преобразователей ,и выполняют перевод прямоугольных координат найденных векторов в систему координат, связанную с антенной решеткой, по формуле: где cos - направляющий косинус оси O'X' в системе координат OXYZ;cos - направляющий косинус оси O'Y' в системе координат OXYZ;cos - направляющий косинус оси O'Z' в системе координат OXYZ;после вычисления векторов и для всех m плоскостей находят истинный вектор направления на источник сигнала путем усреднения тех векторов, расстояния между которыми, вычисленные по Евклидовой норме, минимальны, затем координаты найденного таким образом вектора переводят в сферическую систему координат, связанную с антенной решеткой, с помощью формул: ;где , , - прямоугольные координаты усредненного вектора направления на источник сигнала;в результате чего получают значения горизонтального A и вертикального E углов, задающих направление на гидроакустический маяк-ответчик.
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
СПОСОБ ОПРЕДЕЛЕНИЯ НАПРАВЛЕНИЯ НА ГИДРОАКУСТИЧЕСКИЙ МАЯК-ОТВЕТЧИК ПО ГОРИЗОНТАЛЬНОМУ И ВЕРТИКАЛЬНОМУ УГЛУ
Источник поступления информации: Роспатент

Показаны записи 1-4 из 4.
10.05.2014
№216.012.c0e8

Сферическая гидроакустическая антенна

Использование: изобретение относится к гидроакустической технике. Сущность: антенна содержит тонкостенную полую сферическую оболочку, пьезоэлектрические преобразователи, опору для крепления антенны к носителю. Сферическая оболочка выполнена акустически прозрачной из пластика с равномерно...
Тип: Изобретение
Номер охранного документа: 0002515133
Дата охранного документа: 10.05.2014
10.05.2014
№216.012.c116

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

Изобретение относится к области гидроакустики. Сущность: в способе определения направления на гидроакустический маяк-ответчик в условиях многолучевого распространения навигационного сигнала определяют направление одновременно в горизонтальной и вертикальной плоскостях на гидроакустический...
Тип: Изобретение
Номер охранного документа: 0002515179
Дата охранного документа: 10.05.2014
27.11.2015
№216.013.942f

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

Изобретение относится к области радиотехники. Способ борьбы с гармонической помехой при автокорреляционном методе приема информации с использованием шумоподобных сигналов включает вычисление комплексных огибающих первого и второго периодов принимаемого сигнала, вычисление с помощью дискретного...
Тип: Изобретение
Номер охранного документа: 0002569554
Дата охранного документа: 27.11.2015
20.12.2015
№216.013.9b55

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

Изобретение относится к области гидроакустической связи. Технический результат - повышение помехоустойчивости передачи информации. На передающей стороне осуществляют формирование синхросигнала и исходного информационного сигнала, имеющих малый коэффициент взаимной корреляции, на приемной...
Тип: Изобретение
Номер охранного документа: 0002571390
Дата охранного документа: 20.12.2015
Показаны записи 1-10 из 11.
10.01.2013
№216.012.17e9

Продольно-фрезерный деревообрабатывающий станок

Изобретение относится к деревообрабатывающему оборудованию. Станок содержит станину, загрузочный и подающий механизмы, базовые направляющие, установленные в направлении обрабатывающей зоны станка. Обрабатывающая зона оснащена последовательно расположенными в направлении подачи заготовки...
Тип: Изобретение
Номер охранного документа: 0002471613
Дата охранного документа: 10.01.2013
27.08.2013
№216.012.644c

Способ отверждения карбинольного оптического клея "бальзамин-н"

Изобретение относится к отверждению карбинольного оптического клея, для использования в оптической, опто-электронной и электронной промышленности, и предназначено для склеивания оптических деталей из стекла и кристаллов со светоделительными покрытиями и без, светофильтров и поляроидов с...
Тип: Изобретение
Номер охранного документа: 0002491313
Дата охранного документа: 27.08.2013
10.08.2014
№216.012.e871

Трансформатор с трехфазной и круговой обмотками

Изобретение относится к электротехнике и может быть использовано в многофазных полупроводниковых преобразователях, в частности в выпрямителях, инверторах и преобразователях частоты. Технический результат состоит в повышении качества питающего и выходного напряжения управляемого преобразователя,...
Тип: Изобретение
Номер охранного документа: 0002525298
Дата охранного документа: 10.08.2014
10.07.2015
№216.013.60aa

Способ вибрационных испытаний радиоэлектронной аппаратуры

Изобретение относится к испытательной технике. Способ реализуют следующим образом. На испытуемое изделие воздействует гармоническая вибрация, воспроизводимая на вибростенде и являющаяся эквивалентной ударным воздействиям, возникающим при транспортировании изделия. Предварительно расчетом...
Тип: Изобретение
Номер охранного документа: 0002556287
Дата охранного документа: 10.07.2015
20.07.2015
№216.013.6483

Пьезокерамический материал

Изобретение относится к области производства пьезокерамических материалов, предназначенных для изготовления излучателей и приемников ультразвука, электромеханических преобразователей. Технический результат изобретения заключается в повышении температурной стабильности пьезокерамики....
Тип: Изобретение
Номер охранного документа: 0002557278
Дата охранного документа: 20.07.2015
20.11.2015
№216.013.8fc4

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

Изобретение относится к оптическим монтажным кабелям и способам их производства. Согласно способу оптическое волокно подают в зону технологической обработки, где на него наносят буферный слой. После прохождения обработанным изделием системы роликов на него наносят буферный слой и повив...
Тип: Изобретение
Номер охранного документа: 0002568420
Дата охранного документа: 20.11.2015
27.12.2015
№216.013.9e07

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

Изобретение относится к области определения местоположения подвижных подводных объектов технической природы и может быть использовано при поиске и обнаружении подводных аппаратов и платформ. Достигаемый технический результат - увеличение дальности, угла обзора, а также повышение скрытности...
Тип: Изобретение
Номер охранного документа: 0002572085
Дата охранного документа: 27.12.2015
27.08.2016
№216.015.4f3c

Антенный пост в виде башенно-мачтовой конструкции с пониженной радиозаметностью

Изобретение относится к радиотехнике. Особенностью заявленного антенного поста является то, что металлические валы через редукторы и электромагнитную муфту сцепления соединены с возвратными электродвигателями, крепящимися к нижним бимсам, радиопрозрачные тяги, обеспечивающие продвижение...
Тип: Изобретение
Номер охранного документа: 0002595325
Дата охранного документа: 27.08.2016
27.08.2016
№216.015.51ab

Устройство защиты от импульсных помех

Изобретение относится к радарным системам с защитой от активных импульсных непреднамеренных радиопомех (НРП) радиоэлектронных средств (РЭС), расположенных на одном объекте. Достигаемый технический результат - уменьшение потерь полезной информации в защищаемых от импульсных НРП приемных...
Тип: Изобретение
Номер охранного документа: 0002596148
Дата охранного документа: 27.08.2016
13.01.2017
№217.015.6b85

Астронавигационная система

Изобретение относится к высокоточным астроинерциальным навигационным системам для применения в составе пилотируемых и беспилотных летательных аппаратов. Астронавигационная система, установленная на летательном аппарате, содержит бесплатформенную инерциальную навигационную систему, включающую...
Тип: Изобретение
Номер охранного документа: 0002592715
Дата охранного документа: 27.07.2016
+ добавить свой РИД