Модель тракта прослушивания гидроакустических сигналов

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

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

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

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

Процесс преобразования акустической энергии в электрическую (и наоборот) выполняют подводные электроакустические приёмники и излучатели, входящие в состав антенны, и называемые гидроакустическими преобразователями (ГАП).

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

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

С начала развития подводного флота тракт прослушивания являлся неотъемлемым и единственным средством обнаружения и классификации сигналов. Первоначально этот тракт состоял из одиночных приемников, которые размещались с левого и правого борта, либо где-то в носу судна. Это были одиночные приемники и все сигналы, которые поступали на вход этих приемников, оператор имел возможность прослушать. И по разности прихода сигнала с разных гидрофонов оператор определял, с какого борта идет цель. Тракт прослушивания существует столько же, сколько существует и подводное плавание, так как есть необходимость прослушивать сигналы, обнаруживать, откуда они пришли, и определять, что это за цель, противник либо безопасный объект. Несмотря на наличие в современных ГАС режима объективной классификации, значение тракта прослушивания так же велико. Гидроакустики во время поиска объекта хотят иметь возможность самостоятельно прослушивать сигнал и чаще всего именно они и оценивают параметры цели.

Таким образом, во всех ТЗ на ГАС и ГАК обязательно присутствует требование о наличии тракта прослушивания сигналов и помех с любого направления в заданном секторе обзора ГАС.

К числу основных задач, решаемых гидроакустическими средствами подводных лодок (ПЛ) при освещении окружающей обстановки, относится освещение подводной, надводной и воздушной обстановки в интересах самообороны ПЛ. Для наблюдения за морскими целями ПЛ должна обладать мощным гидроакустическим комплексом, работающим преимущественно в пассивных режимах (шумопеленгование, обнаружение гидроакустических сигналов - ОГС).

При наблюдении за ПЛ ПЛО возникает дуэльная ситуация “ПЛ против ПЛ”, выигрыш в которой определяется соотношением комплекса параметров противоборствующих ПЛ и их гидроакустических комплексов (ГАК): акустическая шумность ПЛ, уровень корабельных акустических помех работе собственного ГАК, энергетический потенциал ГАК, совершенство цифрового вычислительного комплекса ГАК, в том числе алгоритмов и программного обеспечения.

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

Целью дипломного проектирования является разработка тракта прослушивания для ГАС обнаружения гидроакустических сигналов по их шумоизлучению в звуковом диапазоне частот. Тракт обнаружения проектируется для работы в трех частотных диапазонах, рассчитанных под обнаружение целей на различных дальностях, при этом необходимо обеспечить наилучшие условия для прослушивания сигнала цели для каждого из трех частотных диапазонов. Поскольку, как известно, человеческое ухо наилучшим образом воспринимает частоты в полосе от ~300 Гц до 3-4 кГц, возникает задача понижать частоту прослушивания, то есть нужно выделить (вырезать) требуемую частотную полосу и перенести ее (гетеродинировать) в область частот, наиболее комфортную для оператора.

В данной работе требуется:

Разработать структуру тракта прослушивания гидроакустических сигналов на выходе сформированного пространственного канала (канала наблюдения) в тракте шумопеленгования с использованием многоэлементной антенной решетки;

разработать программный макет тракта прослушивания;

установить взаимосвязь основных параметров тракта прослушивания с базовыми параметрами тракта ШП;

с использованием программного макета выбрать параметры тракта прослушивания применительно к заданным в ТЗ на проект условиям.

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

Сектор обзора - ±45° относительно нормали к антенной решетке.

Частота дискретизации входных выборок fd=24000 Гц.

Частотные диапазоны прослушивания: I – (1-2.5) кГц;

II – (2-5) кГц;

III – (4-8) кГц.

Полоса пропускания усилителя тракта прослушивания (0.3-4.5) кГц.

1 Место тракта прослушивания в структуре режима ШП типовой ГАС

Обобщенная структурная схема ГАС шумопеленгования представлена на рисунке 1.

Рис.1 Структурная схема ГАС шумопеленгования аналогового типа

антенное устройство, содержащее несколько гидрофонов.

предварительный широкополосный усилитель.

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

основные усилители и частотные фильтры

детекторы с интеграторами

пороговые устройства

устройство измерения угловых координат

устройство вторичной обработки информации

Непрерывные электрические сигналы от элементов АР поступают на вход аппаратуры предварительной обработки (АПО), основными функциональными узлами которой являются предварительные усилители (ПУ) и диапазонные фильтры (ДФ). С точки зрения системного построения ГАС, АПО разбивается на ряд отдельных трактов, при этом каждой антенной решетке или группе однородных антенных решеток соответствует свой тракт. В отдельных случаях в состав АПО включают устройства формирования характеристик направленности на базе аналоговых или дискретно-аналоговых линий задержки, то есть элементы, относящиеся к первичной обработке.

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

В этом случае дискретные по пространству (принятые дискретной АР) процессы с выхода АПО подвергаются дискретизации с помощью аналого-цифровых преобразователей (АЦП) и поступают на вход ЦВК, где реализуются алгоритмы первичной и вторичной цифровой обработки сигналов.

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

- пространственную фильтрацию сигналов (формирование ХН антенной решетки);

- временную обработку сигналов (спектральный анализ и оценивание).

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

Сопровождение обнаруженных целей по углу с точным измерением текущего пеленга и ВИП в интересах целеуказания оружию осуществляется системой АСЦ.

Система вторичной обработки оперирует уже не с сигналами, принятыми антенной ГАС, а с формализованной информацией о фактах обнаружения сигналов, превысивших порог, а также о параметрах сигналов. К числу основных задач, решаемых аппаратурой вторичной обработки информации, относятся: классификация целей, определение их координат и параметров движения, формирование трасс движения целей и их анализ. В системах вторичной обработки накапливается информация от всех подсистем ГАК, в связи с чем возникает задача комплексной обработки данных. Аппаратура вторичной обработки должна подготовить данные для индикации их на дисплеях системы отображения, регистрации, документирования и управления (СОРДиУ).

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

При классификации в режиме ШП процедура начинается при малых отношениях сигнал/помеха и продолжается в течение всех последующих этапов наблюдения за обстановкой. По мере сближения с целью и увеличения ОСП появляется возможность отнести обнаруженную цель к одному из следующих классов: торговое судно, военный корабль, ПЛ, торпеда, сформированная средствами ГПД ложная цель; шумы портов, гаваней, плавучих платформ; морские животные.

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

-большой корабль или малый;

-есть кавитация или нет;

-шум дизеля или иного механизма;

-число оборотов винтов и число лопастей винта;

-имеет ли место перекладка рулей;

-наблюдается ли пуск оружия;

-не обусловлен ли наблюдаемый сигнал использованием средств ГПД;

-не имеет ли обнаруженный сигнал биологическую или сейсмическую природу;

-как ведет себя цель;

-насколько велик уровень сигнала и в скольких ХН он наблюдается.

Кроме того, на основании наблюдаемого значения ВИП (величина изменения пеленга) оператор может определить дополнительно:

-близко либо далеко находится цель;

-закономерно ли себя ведет пеленг на цель;

может ли цель идти так быстро или так медленно.

Структурная схема режима ШП типовой ГАС цифрового типа представлена на рисунке 2.

Рис 2. Структура режима ШП типовой ГАС цифрового типа

Тракт АСЦ (тракт автоматического сопровождения цели) содержит:

- блок формирования двух диаграмм направленности по половинам апертурного окна;

- Дискриминатор угла, вырабатывающий оценку рассогласования между углом наведения тракта АСЦ и углом прихода сигнала;

- Сглаживающий фильтр с поправками на крен и дифферент;

- блок формирования суммарного канала в частотной области.

Информация с выхода тракта АСЦ поступает в тракт прослушивания и в систему классификации.

В тракте одновременного секторного обзора для обнаружения широкополосных сигналов (ОСО – ШП) выполняется:

- формирование веера каналов наблюдения (пространственных каналов – ПК);

- квадратичное детектирование;

- формирование частотных диапазонов;

- накопления по времени;

- предындикаторная обработка.

Информация с одного из выходов системы формирования веера каналов наблюдения тракта ОСО-ШП также может поступать в тракт прослушивания.

Таким образом, оператор может выбирать, какой сигнал слушать: тракта АСЦ или тракта ОСО-ШП.

Общая структура обработки сигнала в тракте прослушивания представлена на рисунке 3

Рис. 3 Общая структура обработки сигнала в тракте прослушивания

-Аналого-цифровой преобразователь

-Быстрое преобразование Фурье

-Формирование Характеристик направленности

-Сдвиг полосы частот

-Обратное быстрое преобразование Фурье

-Стыковка реализаций

-Цифро-аналоговое преобразование

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

2 Формирование канала наблюдения в частотной области

Задачей формирования канала наблюдения является сбор энергии по большой апертуре, с компенсацией задержки, возникшей на элементах АР при прохождении сигнала в среде. Это формирование можно делать и во временной области, но для получения хороших характеристик такого формирования частота дискретизации входных выборок должна быть очень большой (fd>>fв). Поэтому чаще процедуру формирования пространственного канала выполняют в частотной области, где частота дискретизации fd может быть лишь в 2.5-3 раза выше верхней частоты полосы обработки.

Основана эта процедура на следующем свойстве преобразования Фурье.

Пусть сигнал имеет непрерывный спектр , тогда задержанный сигнал на время t сигнал S1(t)=имеет спектр , который на каждой частоте f отличается от спектра исходного сигнала лишь фазой, которая пропорциональна задержке t и частоте f:

.

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

,

получив М-кратную копию спектра исходного сигнала.

Расчет задержек прихода фронта волны сигнала, приходящего с направления, определяемого углами (j,q) в сферической системе координат, на элементы АР возникающих проводим по формуле :

Рис. 4 Координаты элементов АР

Где – направляющая косинуса нормали к фронту волны.

Зависит от плоского фронта сигнала и координат приемных элементов АР.

По исходным данным мы строим плотный веер ХН для сигнала приходящего с известного нам направления.

Реализация алгоритма обнаружения широкополосного сигнала предполагает фазирование АР в широкой полосе частот. На практике ширина полосы обработки иногда составляет 2-3 октавы и более. Направление фазирования осей ХН принадлежат секторы обзора (α–r, αr) , в котором формируется 2R+1 пространственных канала с шагом по углу

Рис.5 Веер ХН на линейной АР

Алгоритм фазирования АР на направлении ,, и алгоритм формирования ХН имеет вид :

В пространственно - частотной области:

где выходные величины вычисляются на всех частотах из полосы обработки .

Далее мы будем работать с одним направлением в секторе обзора, с которого у оператора может возникнуть желание прослушать сигнал.

3 Факторы влияющие на восстановление сигнала

3.1 Перекрытие входных выборок в тракте прослушивания

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

Таким образом, при реализации ФХН в частотной области происходит циклический сдвиг реализаций на отдельных приемных элементах друг относительно друга. В результате после сдвига реализации на время τm и перевода сигнала во временную область в начале задержанного фрагмента появятся отсчеты из его конца или, наоборот, отсчеты из начала реализации окажутся в конце ее (в зависимости от знака вносимой задержки). Таким образом, как внутри задержанной реализации, так и на стыках соседних по времени фрагментов возникает разрыв фазы (восстановленного) сигнала. Поэтому в суммарном сигнале (на выходе процедуры ФХН) часть отчетов на краях реализации оказывается некорректной. Число этих отчетов напрямую зависит от максимальной временной задержки, вносимой при формировании ПК. Иллюстрация этого эффекта приведена на рис.6, где изображены реализации тонального сигнала (красная кривая), его задержанная во временной области на величину τ копия (черная линия) и она же, после внесения задержки τ в частотной области и перехода к временному представлению (зеленый пунктир).

Рис. 6 Иллюстрация эффекта, возникающего при восстановлении сигнала после ФХН в частотной области

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

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

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

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

Рис. 7а Исходный сигнал во временной области

Рис. 7б Спектр исходного сигнала

Рис. 8а Стыковка восстановленных во временную области фрагментов сигнала на выходе сформированного канала (без отбрасывания некорректных отсчетов)

Рис. 8б Спектр восстановленного сигнала (стыковка без исключения некорректных отсчетов)

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

Для линейной АР величина перекрытия может быть вычислена по формуле:

N = tмах*fd, где

N – Величина перекрытия

tмах=(М-1)*d*sin(j)/c – максимальная задержка,

fd – частота дискретизации,

М – число приемных элементов линейной АР,

d – Межэлементное расстояние,

j - граница сектора обзора (от нормали к антенной решетке),

с – скорость звука в воде.

Для заданных в ТЗ на дипломную работу параметрах эта величина составляет: ~ 25%

3.2 Сужение полосы частот восстановленного сигнала

Приведенные в предыдущем разделе рисунки соответствует восстановлению сигнала во всей полосе (|ΔF|=fd/2). Реальные сигналы на входе тракта ШП после усиления и фильтрации имеют ограниченную полосу, тракт прослушивания (по заданию) предназначен для прослушивания сигналов в трех частотных диапазонах (шириной 1.5, 3 и 4 кГц), что также влияет на качество восстановленного сигнала. На рисунке 9 изображены входной сигнал (черная пунктирная линия) и восстановленный после вырезания требуемой полосы прямоугольным окном в частотной области (синяя кривая).

Рис. 9 Иллюстрация эффекта, возникающего при вырезании полосы фильтром с прямоугольной ЧХ

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

На рисунке 10 приведена импульсная переходная характеристика (ИПХ) фильтра с прямоугольной частотной характеристикой.

Рис. 10 ИПХ фильтра с прямоугольной частотной характеристикой

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

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

На следующем рисунке 11 приведены ИПХ некоторых фильтров (частотных окон): красная линия – окно Хана, черная – окно Хэмминга, зеленая – Кайзера.

На рисунке 12 представлены соответствующие амплитудно-частотные характеристики этих окон: красная линия – окно Хана, черная – окно Хэмминга, зеленая – Кайзера.

Рис. 11 ИПХ фильтров Ханна, Хэмминга, Кайзера

Рис. 12 Вид частотных окон Ханна, Хэмминга и Кайзера

Все представленные окна имеют общий параметр – размер окна L, который равен 27 частотным отсчетам, что составляет примерно 5 % от полосы принимаемого сигнала.

На рисунке 13 приведен результат восстановления сигнала при применении частотного окна Хэмминга:

Чем больше n, тем большая доля энергии, сосредоточенная в главном лепестке спектра (и тем шире этот главный лепесток), и тем меньше уровень боковых лепестков. Черная пунктирная линия – исходный сигнал, синяя – восстановленный сигнал при применении частотного окна Хэмминга.

Рис. 13 Иллюстрация эффекта, возникающего при вырезании полосы фильтром Хэмминга

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

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

3.3 Частотный сдвиг при понижении полосы частот прослушивания

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

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

Рис.14а. Спектр восстановленного сигнала после сдвига его полосы в частотной области на 2Δf. Перекрытие входных выборок 50%

Рис. 14б. Спектр восстановленного сигнала после сдвига его полосы в частотной области на 3Δf. Перекрытие входных выборок 50%

Рис 15а. Состыкованная реализация восстановленного сигнала после сдвига его полосы в частотной области на 4Δf. Перекрытие входных выборок 50%

Рис 15б. Состыкованная реализация восстановленного сигнала после сдвига его полосы в частотной области на 3Δf. Перекрытие входных выборок 50%

Рис. 16а Спектр восстановленного сигнала после сдвига его полосы в частотной области на 4Δf. Перекрытие входных выборок 25%

Рис. 16б Спектр восстановленного сигнала после сдвига его полосы в частотной области на 4Δf. Перекрытие входных выборок 17%

Рис 17а. Состыкованная реализация восстановленного сигнала после сдвига его полосы в частотной области на 4Δf. Перекрытие входных выборок 25%

Рис 17б. Состыкованная реализация восстановленного сигнала после сдвига его полосы в частотной области на 5Δf. Перекрытие входных выборок 25%

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

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

Таблица 1 . Взаимосвязь минимальной полосы частотного сдвига от величины перекрытия входных выборок

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

- обязательное исключение части отчетов в восстановленной реализации сигнала

из дальнейшей обработки;

- применение частотного фильтра с ЧХ, отличной от прямоугольной;

- расширение полосы обработки.

При этом величина перекрытия входных выборок должна быть больше суммы максимальной задержки сигнала на элементах АР и длительности ИПХ частотного фильтра. Других ограничений на величину перекрытия не накладывается.

Традиционно информация, используемая в тракте шумопеленгования, берется с перекрытием входных выборок на величину η. Это реализуется в интересах следующих задач и трактов гидроакустического комплекса:

тракты автоматического сопровождения цели (АСЦ) и выработки классификационной информации по ним;

тракт обнаружения дискретных составляющих сигнала во всем секторе обзора;

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

Однако для тракта прослушивания сигналов и помех, в котором для удобства оператора реализована (в цифровой области) возможность понизить частоту прослушивания, оказалось, что величина перекрытия входных выборок должна быть делителем числа n - количество точек БПФ.

Поскольку для процедуры БПФ обычно выбирают число n, являющееся степенью 2: n=2x ,

где x – целое положительное число, то и величина перекрытия должна быть тоже степенью 2, так как других простых делителей у n нет.

4 Программный макет тракта прослушивания

4.1 Структурная схема алгоритма обработки в тракте прослушивания

Структурная схема алгоритма обработки в тракте прослушивания цифровой ГАС с формированием каналов наблюдения в частотной области приведена на рисунке 18.

Рис. 18 Структурная схема алгоритма обработки в тракте прослушивания

1) АР – антенная решетка, M – количество ее приемных элементов;

2) АПО – фильтрация сигнала, предварительное усиление;

3) АЦП – аналого-цифровое преобразование с частотой дискретизации 24000Гц, процессов, принятых на М приемных элементах антенны; получаем дискретизированные по пространству и по времени выборки входного поля;

4) Набор реализаций с перекрытием 25 % точек;

5) БПФ – n-точечное быстрое преобразование Фурье по всем М каналам. Частотное разрешение системы df=fd/n;

6) ФПК – формирование пространственного канала: осуществляется в частотной области умножением на фазирующие коэффициенты и суммированием выходов задержанных реализаций. Получаем один выход в полосе от 1 до 8 кГц, границам полосы частот соответствуют номера частотных отсчетов Кн=(1000/df)=21 и Кв=(8000/df+0.5)=171.

7) Вырезание полосы частот, соответствующей заданному оператором номеру частотного диапазона Nd (от 1 до 3):

I чд - от Кн= 21 до Кв=53,

II чд - от Кн= 42 до Кв=106,

III чд - от Кн= 85 до Кв=171;

8) Сдвиг полосы в область от 0.3 кГц (для удобства оператора); сдвиг осуществляется на величину:

4Δf = 187 Гц - в I чд,

32Δf = 1497 Гц – во II чд,

76Δf = 3556 Гц – в III чд;

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

K=53 – в I чд,

K=81 – во II чд,

K=99 – в III чд ;

10) Восстановление сигнала во временную область с использованием процедуры ОБПФ на те же самые n точек;

11) Отбрасывание некорректных отсчетов – по n/8 точек в начале и в конце реализации;

12) Стыковка реализаций;

13) Цифро-аналоговое преобразование (ЦАП) с частотой дискретизации 24000кГц;

14) Вывод результатов на динамик или