СПОСОБ ПОДАВЛЕНИЯ ПОМЕХ

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

Известен способ подавления помех ("Применение цифровой обработки сигналов" под ред Э.Опенгейма, М, "Мир", 1980), в котором из принятых радиолокатором сигналов в N периодах повторения когерентной пачки формируют отсчеты комплексной огибающей, отсчеты в каждом элементе дальности обрабатывают в режекторном фильтре с весовыми коэффициентами W=(w₁, w₂, …, w_N)^T, заранее выбранными таким образом, чтобы обеспечить минимальное значение мощности помех на выходе режекторного фильтра. В данном способе обеспечивается подавление некоторых видов пассивных помех, например, помех от неподвижных МП. В этом случае используют биномиальные весовые коэффициенты, находящиеся из выражения .

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

Наиболее близким техническим решением является способ подавления помех (Д.И.Попов. Синтез цифровых адаптивных режекторных фильтров. Радиотехника, 1981 г., т.36, №10), в котором из принятых радиолокатором сигналов в N периодах повторения когерентной пачки формируют отсчеты комплексной огибающей, при обработке отсчетов в k-м элементе дальности оценивают матрицу ковариации помех по отсчетам в N_w элементах дальности, окружающих обрабатываемый k-й элемент, определяют весовые коэффициенты режекторного фильтра W=(w₁, w₂, …, w_N)^T из уравнения

где I - единичная матрица размером N×N, λ_min - минимальное собственное значение матрицы размером N×N, обрабатывают отсчеты k-го элемента дальности в данном режекторном фильтре. В указанном способе обеспечивается подавление разного вида помех, в том числе смеси пассивных помех от неподвижных МП и движущихся ПП При этом достигается минимальное значение мощности смеси помех на выходе режекторного фильтра.

Основной недостаток этого способа заключается в сложности его технической реализации, связанной с тем, что из-за изменения во времени матрицы ковариации смеси помех необходимо в каждом элементе дальности получать новое решение матричного уравнения (1). Это требует производить порядка ~N³ операций комплексного умножения за время, равное периоду дискретизации комплексной огибающей τ_d, что для типовых значений N>3 и τ_d<1 мкс вызывает значительные технические трудности.

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

Технический результат (решаемая задача) в предлагаемом способе достигается тем, что в известном способе подавления помех, представляющих собой смесь помех от неподвижных местных предметов (МП) и движущихся пассивных помех (ПП), в котором из принятых радиолокатором сигналов и помех в N периодах повторения когерентной пачки, где N≥3, формируют отсчеты комплексной огибающей x_i,k, где i - номер периода повторения когерентной пачки, 1≤i≤N, k - номер элемента дальности, формируют окно из N_w элементов дальности, окружающих k-й элемент, в котором оценивают корреляционные параметры помех, включая коэффициенты корреляции помех между соседними периодами повторения когерентной пачки, оценки корреляционных параметров помех используют для фильтрации сигналов на фоне смеси помех от МП и ПП, согласно изобретению отсчеты комплексной огибающей x_i,k подают на N каналов обработки, причем в j-м канале, где j=1, …, N, при j≥2 их вначале подают на режекторный фильтр подавления помех от неподвижных МП (j-1)-го порядка с биномиальными весовыми коэффициентами , где m=1, …j, затем из отсчетов z_i,k(j) на выходе указанного фильтра формируют упомянутое окно из N_w элементов дальности, окружающих k-й элемент, в котором оценивают упомянутый коэффициент корреляции помех между соседними j-м и j+1-м периодами повторения когерентной пачки где - соответственно оценки модуля и фазы коэффициента корреляции помех, отсчеты z_i,k (j) поворачивают на фазу в ρ-м периоде повторения когерентной пачки, где ρ=j, …, N, затем при j≤N-1 их подают на режекторный фильтр подавления помех от движущихся ПП с действительными коэффициентами w_2-n, где n=1, …, N-j+1, которые находят из предварительно рассчитанных для подавления ПП с нулевой фазой коэффициента корреляции, гауссовой формой спектра и разных значений модуля коэффициента корреляции ρ₀, заданных с точностью Δρ₀ в диапазоне от 0 до 1, так, чтобы расчетное значение ρ₀ отличалось от полученной оценки нe более чем на величину точности Δρ₀, на выходах режекторных фильтров подавления помех от движущихся ПП в N каналах обработки измеряют мощность помех в окне из N_p элементов дальности, окружающих k-й элемент, где N_p выбирается исходя из допустимого снижения подавления помех, по минимальной оценке мощности помех выбирают канал обработки, выполняющий упомянутую функцию фильтрации сигналов на фоне смеси помех от МП и ПП.

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

- отсчеты комплексной огибающей x_1,k подают на N каналов обработки, причем в j-м канале, где j=1,…,N, при j≥2 их вначале подают на режекторный фильтр подавления помех от неподвижных МП (j-1)-го порядка с биномиальными весовыми коэффициентами , где m=1,…j,

- из отсчетов z_1.k (j) на выходе режекторного фильтра подавления помех от неподвижных МП формируют упомянутое окно из N_w элементов дальности, окружающих k-й элемент, в котором оценивают упомянутый коэффициент корреляции помех между соседними j-м и j+1-м периодами повторения когерентной пачки , где соответственно оценки модуля и фазы коэффициента корреляции помех;

- отсчеты z_i,k(j) поворачивают на фазу в ρ-м периоде повторения когерентной пачки, где p=j, …,N, затем при j≤N-l их подают на режекторный фильтр подавления помех от движущихся ПП с действительными коэффициентами w_2n, где n=1, …, N - j+1,

- коэффициенты w_2n находят из предварительно рассчитанных для подавления ПП с нулевой фазой коэффициента корреляции, гауссовой формой спектра и разных значений модуля коэффициента корреляции ρ₀, заданных с точностью Δ_ρ0 в диапазоне от 0 до 1, так, чтобы расчетное значение ρ₀ отличалось от полученной оценки не более чем на величину точности Δ_ρ0;

- на выходах режекторных фильтров подавления помех от движущихся ПП в N каналах обработки измеряют мощность помех в окне из N_p элементов дальности, окружающих k-ый элемент, где N_p выбирается исходя из допустимого снижения подавления помех;

- по минимальной оценке мощности помех выбирают канал обработки, выполняющий упомянутую функцию фильтрации сигналов на фоне смеси помех от МП и ПП.

Применение всех новых признаков совместно с признаками прототипа позволит упростить техническую реализацию способа. В прототипе необходимо при обработке сигналов в каждом элементе дальности оценивать всю матрицу ковариации помех , состоящую в общем случае из различных элементов. Затем требуется найти решение уравнения (1), для чего необходимо использовать значительные вычислительные ресурсы В предлагаемом способе подавления помех обработка производится в N канальной системе, причем в j-м канале отсчеты комплексной огибающей вначале подают на режекторный фильтр (j-1)-го порядка с биномиальными коэффициентами, в котором производится подавление помех от неподвижных МП. Подавление движущейся ПП в этом режекторном фильтре не обеспечивается, для этого выполняют следующее. На выходе упомянутого режекторного фильтра подавления помех от неподвижных МП в окне дальности, окружающем k-й элемент дальности, измеряют коэффициент корреляции помех , далее отсчеты z_i,k поворачивают таким образом, чтобы компенсировать фазу движущейся ПП затем подают их на режекторный фильтр (N-j)-го порядка, весовые коэффициенты которого предварительно рассчитывают для подавления ПП с нулевой фазой коэффициента корреляции, гауссовой формой спектра и разных значений модуля коэффициента корреляции ρ₀. Из рассчитанных значений весовых коэффициентов выбирают такие, чтобы расчетное значение ρ₀ отличалось от полученной оценки не более чем на величину точности Δρ₀. В последнем упомянутом режекторном фильтре обеспечивается подавление движущейся ПП.

В зависимости от соотношения мощностей МП и движущейся ПП наибольшее подавление смеси помех обеспечивается в одном из каналов обработки, отличающихся порядком режекторного фильтра подавления МП - 0, …, N-1 и порядком режекторного фильтра подавления движущейся ПП - соответственно, N-1, …, 0. (Режекторный фильтр 0-го порядка представляет собой звено с единичным коэффициентом передачи на всех частотах). Выбор наилучшего канала обработки производится путем сравнения оценок мощности помех в окне дальности на их выходах. Для реализации предложенного способа требуются существенно меньшие вычислительные ресурсы.

Изобретение поясняется чертежами:

- Фиг.1 - Структурная схема устройства, реализующего предлагаемый способ подавления помех.

- Фиг.2 - Структурная схема устройства, реализующего адаптивный режекторный фильтр на Фиг.1.

- Фиг.3 - Структурная схема устройства, реализующего измеритель мощности на Фиг.1.

- Фиг.4 - Структурная схема устройства, реализующего блок определения канала с минимальным значением на Фиг.1.

- Фиг.5 - Зависимости коэффициента улучшения отношения сигнал/помеха от частоты Доплера движущейся ПП (прототип).

- Фиг.6 - Зависимости коэффициента улучшения отношения сигнал/помеха от частоты Доплера движущейся ПП (изобретение).

Техническая реализация предлагаемого способа подавления помех возможна на основе устройства, показанного на Фиг 1. Устройство содержит формирователь отсчетов комплексной огибающей 1, вход которого является входом устройства, (N-1) режекторных фильтров подавления МП 2, (N-1) адаптивных режекторных фильтров подавления ПП 3, N измерителей мощности 4, блок определения канала с минимальным значением 5, мультиплексор 6, выход которого является выходом устройства.

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

Техническая реализация адаптивного режекторного фильтра подавления ПП 3 на Фиг.1 возможна в виде устройства, представленного на Фиг.2. Устройство содержит первый умножитель 7, блок задержки на период 9, блок оценки коэффициента корреляции помех 10, входы которых соединены между собой и являются входом устройства, режекторный фильтр подавления ПП 8, выход которого является выходом устройства, регистр 11, второй умножитель 7, блок памяти весовых коэффициентов режекторного фильтра 12.

Техническая реализация измерителя мощности 4 на Фиг.1 возможна в виде устройства, представленного на Фиг.3. Устройство содержит вычислитель квадрата модуля 13, вход которого является входом устройства, (N_p+1)-отводный регистр 14, сумматор 15, делитель 16, выход которого является выходом устройства.

Техническая реализация блока определения канала с минимальным значением 5 на Фиг.1 возможна в виде устройства, представленного на Фиг.4. Устройство содержит схему выбора минимального значения 17, N входов которой являются входами устройства, N схем сравнения 18 и шифратор 19, выход которого является выходом устройства.

Устройство на Фиг.1 работает следующим образом. На формирователь отсчетов комплексной огибающей 1 поступают сигналы x(t) с выхода приемного устройства РЛС. С выхода формирователя отсчетов комплексной огибающей 1 отсчеты комплексной огибающей x_i,k поступают на N каналов обработки, причем j-й канал (j=1, …, N} состоит из последовательно соединенных режекторного фильтра подавления МП 2 (j-1)-го порядка, адаптивного режекторного фильтра подавления ПП 3 (N-j)-го порядка и измерителя мощности 4. Причем в первом канале обработки режекторный фильтр подавления МП 2 отсутствует и отсчеты комплексной огибающей с выхода формирователя отсчетов комплексной огибающей 1 подаются прямо на адаптивный режекторный фильтр подавления ПП 3, а в N-м канале обработки отсутствует адаптивный режекторный фильтр подавления ПП 3 и отсчеты z_i,k с выхода режекторного фильтра подавления МП 2 подаются прямо на измеритель мощности 4. Оценки мощности помех , полученные с помощью измерителей 4 на выходах адаптивных режекторных фильтров подавления ПП 3, подаются на блок определения канала с минимальным значением 5, с выхода которого снимается и подается на управляющий вход мультиплексора 6 номер канала обработки I_m, на выходе которого мощность помех имеет минимальное значение На N информационных входов мультиплексора 6 подаются сигналы u_k,j с выходов адаптивных режекторных фильтров подавления ПП 3 соответствующих каналов обработки. С выхода мультиплексора 6 снимается смесь сигналов и подавленных помех, полученная на выходе адаптивного режекторного фильтра подавления ПП 3 канала обработки j_m, в котором обеспечивается наибольшее подавление помех.

Формирователь отсчетов комплексной огибающей 1 представляет собой устройство, известное, например, из (Я.Д.Ширман, В.Н.Манжос. Теория и техника обработки радиолокационной информации на фоне помех, М., Радио и связь, 1981, 416 с.).

Режекторные фильтры подавления МП 2 представляют собой j-отводные трансверсальные фильтры, известные, например, из (Л.Рабинер, Б.Гоулд. Теория и применение цифровой обработки сигналов, М., "Мир", 1978). Весовые коэффициенты режекторных фильтров подавления МП 2 определяются из выражения , где m=1, …, j.

Устройство на Фиг.2 работает следующим образом. В (j+1)-м периоде повторения когерентной пачки отсчеты z_j+1,k(j) с выхода режекторного фильтра подавления МП 2 j-го канала обработки поступают на первый вход блока оценки коэффициента корреляции помех 10, вход блока задержки на период 9 и первый вход первого умножителя 7. Задержанные на период повторения когерентной пачки T отсчеты z_j,k(j) с выхода блока задержки на период 9 поступают на второй вход блока оценки коэффициента корреляции помех 10, в котором в окне из N_w элементов дальности, окружающих k-й элемент, оценивают коэффициент корреляции помех между j-м и j+1-м периодами повторения когерентной пачки . Оценка фазового вектора помех с первого выхода блока оценки коэффициента корреляции помех 10 подается на первый вход второго умножителя 7, на второй вход которого с выхода регистра 11 поступает хранящийся в нем вектор поворота отсчетов по фазе ρ-ом периоде повторения когерентной пачки , где p≥j. С выхода второго умножителя 7 новое значение вектора поворота отсчетов по фазе записывается в регистр 11 для использования в следующем (ρ+1)-м периоде повторения когерентной пачки. Вектор поворота отсчетов по фазе в ρ-м периоде повторения когерентной пачки с выхода регистра 11 подается на второй вход первого умножителя 7, с выхода которого отсчеты с нулевой после поворота доплеровской фазой ПП подаются на информационный вход режекторного фильтра подавления ПП 8. Оценка модуля коэффициента корреляции помех co второго выхода блока оценки коэффициента корреляции помех 10 подается на вход блока памяти весовых коэффициентов режекторного фильтра 12, с выхода которого весовые коэффициенты режекторного фильтра w_2,1, w_2,2, …, w_2,N-j+1, соответствующие значению , подаются на управляющий вход режекторного фильтра подавления ПП 8. В режекторном фильтре подавления ПП 8 происходит оптимальное подавление движущейся ПП с нулевой доплеровской фазой, гауссовым спектром и модулем коэффициента корреляции Смесь сигналов и подавленных помех с выхода режекторного фильтра подавления ПП 8 поступает на выход данного устройства.

В блоке оценки коэффициента корреляции помех 10, известного, например, из (Д.И.Попов. Синтез цифровых адаптивных режекторных фильтров. Радиотехника, 1981 г., т.36, №10), определяется оценка коэффициента корреляции между отсчетами комплексной огибающей j-го и (j+1)-го периодов пачки повторения когерентной пачки по следующей формуле для j-го канала обработки:

Режекторный фильтр подавления ПП 8 аналогичен режекторному фильтру подавления МП 2 с постоянными весовыми коэффициентами, реализация которого известна, за исключением наличия управляющего входа, по которому передаются значения используемых в данный момент времени весовых коэффициентов трансверсального фильтра w_2,1, w_2,2, …, w_2,N-j+1.

Устройство на Фиг.3 работает следующим образом. На вычислитель квадрата модуля 13 поступают сигналы u_k,j c выхода адаптивного режекгорного фильтра подавления ПП 3. Далее с выхода вычислителя квадрата модуля 13 значения |u_k,j|² записываются в (N_p+1)-отводный регистр 14, с отводов которого сигналы подаются на сумматор 15 и далее на первый вход делителя 16, на второй вход которого подается коэффициент N_p. Величина среднего значения квадратов модулей сигналов с выхода делителя 16 подается на выход данного устройства.

Устройство на Фиг.4 работает следующим образом. Оценки мощности помех с выхода измерителей мощности 4 поступают на 1,…,N входы схемы выбора минимального значения 17 и первые входы N схем сравнения 18. Минимальная оценка мощности помех с выхода схемы выбора минимального значения 17 подается на вторые входы N схем сравнения 18. На выходах схем сравнения 18 единичный уровень появляется при равенстве входных сигналов, то есть при , в противном случае формируется нулевой уровень. С выходов N схем сравнения 18 логические сигналы подаются на 1, 2, …, N входы шифратора 19, на выходе которого формируется число, равное номеру входа с единичным логическим уровнем, в данном случае это номер I_m.

Расчеты характеристик предлагаемого способа подавления помех проводились методом статистического моделирования в среде программы MatLab v.7. На Фиг.5 для прототипа показаны зависимости коэффициента улучшения отношения сигнал/помеха К_u от значения центральной частоты спектра движущейся ПП ƒ₀, умноженного на период повторения когерентной пачки Т (ƒ₀×T), для разных значений отношения мощностей МП и движущейся ПП G_mp. Число импульсов когерентной пачки N=4, ширина спектра движущейся ПП Δƒ=50 Гц. Аналогичные зависимости для изобретения показаны на Фиг.6. Значения коэффициента улучшения отношения сигнал/помеха К_u изобретения, характеризующего степень подавления помех, близки к значениям прототипа.

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