Способ наблюдения за поверхностью на базе многоканальной бортовой рлс

Способ наблюдения за поверхностью на базе многоканальной бортовой рлс

Реферат

Изобретение относится к радиолокации, а именно к радиолокационным системам наблюдения за поверхностью (и объектами на поверхности) на базе бортовой РЛС, работающей в режиме "реального луча" (РЛ) с многоканальной приемной системой, где многоканальность достигается или наличием большого числа пространственно разнесенных приемных элементов типа фазированной антенной решетки (ФАР) или за счет частотного (фазового) сканирования излучаемого сигнала [1].

При наблюдении бортовой моноимпульсной РЛС за наземными радиоконтрастными объектами в режиме РЛ осуществляется построчное сканирование лучом РЛС заданного участка поверхности путем последовательного смещения луча по азимуту на малую часть ширины диаграммы направленности антенны (ДН) на уровне 0,5 мощности с последующей алгоритмической обработкой принятых сигналов, прошедших амплитудное детектирование в элементах разрешения дальности с целью формирования двумерного радиоизображения поверхности в координатах дальность - азимут с повышенной разрешающей способностью по азимуту [2, 3].

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

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

Наиболее близким по технической сущности является способ синтезирования новых элементов разрешения по азимуту в режиме РЛ [3], который заключается в следующем. Повышение разрешающей способности с расширением зоны обзора РЛС по азимуту и формирование матрицы двумерного РИ поверхности в координатах дальность - азимут достигается за счет быстрого (электронного или механического) переключения (смещения) луча РЛС по азимуту (по j) на величину (2n+1)-й части ширины ДН размером в 2n+1 элементов дискретизации и обработки амплитуд отраженных сигналов РЛС, полученных при разных положениях луча на выходе суммарного и разностного каналов после амплитудного детектирования в элементах разрешения дальности, которая заключается в следующем.

1. Амплитуды y_q(i,j+j₁), , n₁≥n, полученные в суммарном (q=1) и разностном (q=2) каналах РЛС при 2n₁+1 j₁-x положениях луча (относительно центрального j-го направления по азимуту) в i-x элементах разрешения дальности суммируются с весами h_q(j₁), , которые вычисляются заранее по определенной методике. Результатом такой обработки являются оценки , , амплитуд сигнала отражения от поверхности в i-x элементах дальности, соответствующие центральному элементу дискретизации ДН при j-м положении луча:

где Q - число используемых каналов измерения (q=2); I - число элементов дальности в матрице изображения, соответствующих зоне обзора по дальности.

2. Для расширения зоны обзора по азимуту увеличивается число j₁-x сканирований луча по азимуту относительно j-го центрального направления: (N₁>n₁) и соответственно увеличивается число измерений: , . Это приводит к появлению 2l+1 оценок (l=N₁-n₁):

3. Совокупность оценок амплитуд x(i,j) сигналов, отраженных от соответствующих i,j-x элементов поверхности, представляет матрицу А двумерного амплитудного РИ поверхности в зоне обзора размером в I элементов разрешения по дальности и 2l+1 синтезированных элементов разрешения (дискретизации) по азимуту, размеры которых в 2n+1 раз меньше ширины ДН.

Однако такой способ обладает следующими недостатками.

1. Сканирование (смещение) луча в зоне обзора для движущейся (установленной на носителе) РЛС приводит к независимости и случайности фаз сигналов, отраженных от одних и тех же пространственных элементов дискретизации и распределенных по равномерному закону на [0,2π]. Устранение влияния случайности фазы с помощью амплитудного детектирования после прохождения отраженных сигналов тракта первичной обработки и фазового детектирования в квадратурных каналах приводит к существенному увеличению уровня помех и снижению отношения сигнал-шум. Следствием этого является невысокая разрешающая способность изображения в режиме РЛ.

2. Использование данных двух каналов с разными характеристиками ДН (суммарного и разностного) в моноимпульсных РЛС дает возможность одновременно принимать сигналы, отраженные от одних и тех же элементов дискретизации поверхности, и обрабатывать эти сигналы после прохождения квадратурных каналов без амплитудного детектирования, что существенно снижает уровень шумов. Однако число каналов в моноимпульсных РЛС значительно меньше числа оцениваемых параметров поля отражения, что не позволяет достичь необходимой точности оценивания.

3. Для получения начальной оценки амплитуды сигнала отражения в центральном элементе дискретизации требуется накопление 2n₁+1 измерений при 2n₁+1 положениях луча. При непрерывном обзоре поверхности последующие оценки находятся последовательно при каждом новом положении луча. Однако при разрывном обзоре (в разных угловых направлениях) требуется первоначальное накопление измерений, что заметно снижает быстродействие в случае использования механических антенных систем.

Технический результат направлен на повышение разрешающей способности по азимуту в элементах разрешения дальности с расширением зоны обзора РЛС по азимуту и увеличение точности и быстродействия оценивания амплитуд сигналов в синтезированных элементах разрешения азимута.

Технический результат предлагаемого технического решения достигается тем, что способ наблюдения за поверхностью на базе многоканальной бортовой РЛС в режиме РЛ с электронным сканированием заключается в формировании матрицы двумерного радиоизображения поверхности в координатах дальность - азимут, при этом за счет быстрого электронного переключения луча РЛС последовательно смещают луч по азимуту на величину (2n+1)-й части ширины ДН размером в 2n+1 элементов дискретизации на уровне 0,5 мощности и обрабатывают полученные при каждом j-м положении луча в i-x элементах разрешения дальности амплитуды отраженного сигнала, отличающийся тем, что при обработке измеряют амплитуды , отраженного сигнала в квадратурных каналах фазового детектирования (С - косинусном и S - синусном) одновременно в каждом q-м приемном канале антенной системы, состоящей из большого числа Q (Q≥2n+1) разнесенных по фазе приемных элементов, при этом измерения , , , суммируют с весами , , , найденными заранее, тем самым оценивают косинусную х^C(i,j) и синусную х^S(i,j) составляющие амплитуды x(i,j) отраженного сигнала, соответствующего центру j-го луча (центральному элементу дискретизации):

затем полученные оценки возводят в квадрат, суммируют и извлекают корень, тем самым вычисляют оценки амплитуд отраженного сигнала в i-x элементах дальности и j-м синтезированном элементе разрешения азимута:

указанные операции повторяют для всех j-x положений луча по азимуту в зоне обзора и получают матрицу А оценок амплитуд , представляющую двумерное радиоизображение поверхности с повышенным разрешением по азимуту.

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

1. Луч РЛС последовательно смещают по азимуту (по j) на величину (2n+1)-й части ширины ДН размером в 2n+1 элементов дискретизации на уровне 0,5 мощности. Антенная система состоит из большого числа Q (Q≥2n+1) измерительных каналов - приемных элементов, разнесенных по фазе принимаемого сигнала [1].

2. При каждом j-м положении луча в i-x элементах разрешения дальности измеряют амплитуды отраженного сигнала , в квадратурных каналах фазового детектирования (С - косинусном и S - синусном) одновременно в каждом q-м приемном канале .

3. Результаты измерений , , , в каждом i-м элементе дальности суммируют с весами , , , найденными заранее по определенной методике, тем самым оценивают косинусную x^C(i,j) и синусную х^S(i,j) составляющие амплитуды x(i,j) отраженного сигнала, соответствующего центру луча (ДН):

4. Вычисляют оценки амплитуд отраженного сигнала в j-м синтезированном элементе разрешения по формуле:

5. Указанные операции повторяют для всех j-x положений луча по азимуту в зоне обзора и тем самым получают матрицу А оценок амплитуд , , , представляющую амплитудное изображение поверхности в координатах дальность - азимут с повышенной точностью оценивания амплитуд в элементах дискретизации азимута и соответственно повышенным разрешением по азимуту.

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

где Q - число приемных каналов; - сигнал в квадратурных каналах фазового детектирования с измеряемой амплитудой s_q(t) и измеряемой фазой ψ_q(t); - нормированные комплексные коэффициенты ДН q-го канала, характеризующие интенсивность прихода сигналов от j-го углового направления относительно центрального направления; - полезная составляющая сигнала с амплитудой x_j(t), несущей информацию о поле отражения, и фазой ϕ_j(t); Δϕ_q(j) - известный фазовый сдвиг при приеме отраженного сигнала с j-го углового направления q-м приемным элементом; - комплексный гауссовский белый шум, действительная ξ_q(t) и мнимая η_q(t) составляющие которого распределены по нормальному закону с нулевым математическим ожиданием и дисперсией . Амплитуды x_j(t) и фазы ϕ_j(t) в общем случае случайны по j-м элементам дискретизации, а также на множестве положений антенны и их статистические характеристики определены. Случайность фазы обусловлена как движением носителя РЛС, так и тем, что длина волны излучения (например, 8 мм), отражающейся в данном i-м элементе разрешения дальности, меньше размера этого элемента (например, 1 м). Представим (1) в виде

где

Сигнал в (2) , где , , содержит действительную и мнимую составляющие:

где ,

, .

Выражение (3) представляет систему 2Q уравнений с 2N неизвестными и , N=2n+1, причем

После стробирования сигнала в i-x элементах разрешения дальности на промежутке [t_μ-1, t_μ] получается следующая общая модель измерения в q-м канале в i-м элементе дальности при j-м положении луча:

которая используется для восстановления искомого поля X={x(i,j)} на множестве интегральных (суммарных) измерений Y_q={y_q(i,j)}, , зашумленных помехами ξ_q и η_q. Так как корреляцией сигналов в соседних i-x стробах дальности можно пренебречь, то обработка измерений ведется независимо в i-x элементах разрешения дальности.

Для многоканальной антенной системы с Q излучающими и Q приемными элементами модель (1) принимает вид:

где - коэффициенты ДН, характеризующие интенсивность отраженного сигнала в q-м приемном элементе при q₁-м излучающем элементе; , - фазовый сдвиг сигнала от q₁-го излучателя, отраженного в j-м угловом направлении и принятого q₁-м приемным элементом.

После замены в (5) получается выражение (1).

Отношение сигнал-шум в модели (3) можно оценить, представив искомую амплитуду x_j(t) в виде суммы детерминированной (средней) составляющей и случайного отклонения Δx_j(t):

С учетом некоррелированности случайных составляющих и равномерности распределения ϕ_j(t) на [0,2π] отношение сигнал-шум по мощности для косинусной (или аналогично синусной) составляющей модели (4) найдется:

где M - символ математического ожидания; и - дисперсии случайных величин Δx_j(t) и ξ_q(t). Если принять приближенно , для 2n+1=5-7 (k=0,58-0,26), то отношение сигнал-шум по мощности составит , . С учетом осреднения на множестве L повторений измерений отношение сигнал-шум будет равно , что значительно больше, чем для модели амплитудного детектирования:

где отношение сигнал-шум при тех же условиях составляет примерно независимо от мощности полезного сигнала . Дальнейшее увеличение отношения сигнал-шум осуществляется в процессе алгоритмической обработки (3)-(4) за счет избыточного числа каналов измерения: Q>2n+1. Соответственно точность оценивания в многоканальной системе с моделью измерения (3), (4) при раздельной обработке составляющих отраженного сигнала в квадратурных каналах выше, чем в моноимпульсной с моделью (6).

Оптимальное оценивание сводится к следующему. Выражение (2) представляет систему Q линейных уравнений с помехами относительно N=2n+1 неизвестных или в матрично-комплексной форме:

где - Q-вектор комплексных измерений ; - Q×N-матрица комплексных коэффициентов ДН ; - N-вектор комплексных параметров поля отражения ; - Q-вектор комплексных гауссовских помех .

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

Практически удобно работать с действительными выражениями (3)-(4), которые также представляются в матричной форме:

где Y - 2Q-вектор действительных измерений и ; А - (2Q)×(2N)-матрица действительных коэффициентов ДН и ; Х - 2N-вектор действительных параметров поля отражения и , подлежащих оцениванию; Р - 2Q-вектор помех ξ_q и η_q.

Матричная запись (7) в случае некоррелированных помех Р позволяет находить стандартные МНК-оценки 2N-вектора X:

где - матрица весовых коэффициентов; δ - параметр регуляризации, необходимый для обращения плохо обусловленной матрицы А^T А, который с позиции статистической регуляризации [5] для некоррелированных полей имеет смысл отношения дисперсий:

Точность оценивания (8) характеризуется корреляционной матрицей К_ΔX ошибок оценивания . При этом наибольшая точность при малом числе каналов Q (Q≥N) достигается для тех составляющих вектора , которые соответствуют центру j-го луча (j₁=0). Эти составляющие и вычисляются по формулам:

где , - весовые коэффициенты центральной строки матрицы Н, соответствующие наименьшей дисперсии ошибки оценивания (в общем случае зависящие от j-го положения луча вследствие возможного изменения формы ДН при электронном сканировании), и используются для вычисления оценки амплитуды центрального элемента дискретизации азимута в каждом i-м элементе дальности:

В качестве примера при расчете матрицы К_ΔX примем:

,

где d - расстояние между приемными элементами; λ - длина волны; θ_q(j) - угол отклонения луча (от нормали к антенне); Δθ - ширина ДН в градусах (например: Δθ=1° при λ=0,008 м).

При значениях Δθ=1°, d=λ, 2n+1=3 и Q=3 после обращения А^TА с параметром δ=0,1 получается 6×6-матрица К_ΔX, диагональные элементы которой представляют дисперсии (D) оценок и , причем наименьшие их значения получаются при j₁=0:.

Для получения амплитудного изображения в зоне обзора (на множестве элементов дискретизации по углам) осуществляется сканирование луча (электронное или механическое) со смещением на один элемент дискретизации по азимуту и многократно повторяется оценивание (9), (10). При избыточном числе каналов (Q>>2n+1) точность оценивания увеличивается и отпадает необходимость поэлементного сканирования луча. В этом случае для расширения зоны обзора осуществляется сканирование со смещением по азимуту на ширину ДН, а в векторе оценок (8) используются все компоненты.

Предложенный способ позволяет без увеличения энергетических затрат станции повысить разрешающую способность многоканальной РЛС по азимуту в режиме РЛ за счет увеличения точности оценивания параметров поля отражения с расширением зоны обзора и сформировать на основе (9), (10) матрицу радиоизображения поверхности в виде совокупности оценок , , , амплитуд сигналов, отраженных от соответствующих i,j-x элементов дискретизации поверхности, которая позволяет наблюдать на экране индикатора объекты на поверхности в условиях отсутствия оптической видимости с более высоким разрешением по сравнению с известными способами обзора реальным лучом.

Источники информации

1. Воскресенский Д.И. Антенны с обработкой сигнала: Учеб. пособие для вузов. - М.: САЙНС-ПРЕСС, 2002. 80 с.

2. Пат. RU 2249832 С1. Способ наблюдения за поверхностью / В.К.Клочко, Г.Н.Колодько, В.И.Мойбенко, А.А.Ермаков. МПК: G 01 S 13/02, H 01 Q 21/00. Приоритет 02.09.2003. Опубл.: 10.04. 2005. Бюл. №10.

3. Пат. RU 2256193 С1. Способ наблюдения за поверхностью и воздушной обстановкой / В.К.Клочко, Г.Н.Колодько, В.И.Мойбенко, А.А.Ермаков. МПК: G 01 S 13/02. Приоритет 08.12.2003. Опубл.: 10.07. 2005. Бюл. №19.

4. Радиолокационные станции с цифровым синтезированием апертуры антенны / В.Н.Антипов, В.Т.Горяинов, А.Н.Кулин, Е.Ф.Толстов и др. Под ред. В.Т.Горяинова. М.: Радио и связь, 1988. С.13-14.

5. Василенко Г.И., Тараторин A.M. Восстановление изображений. М.: Радио и связь, 1986. С.76-82.

Способ наблюдения за поверхностью на базе многоканальной бортовой РЛС в режиме реального луча с электронным сканированием, заключающийся в формировании матрицы двумерного радиоизображения поверхности в координатах дальность - азимут, при этом за счет быстрого электронного переключения луча РЛС смещают луч по азимуту на величину (2n+1)-й части ширины диаграммы направленности антенны (ДН) размером в 2n+1 элементов дискретизации на уровне 0,5 мощности и обрабатывают полученные при каждом j-м положении луча в i-x элементах разрешения дальности амплитуды отраженного сигнала, отличающийся тем, что при обработке измеряют амплитуды , отраженного сигнала в квадратурных каналах фазового детектирования (С - косинусном и S - синусном) одновременно в каждом q-м приемном канале антенной системы, состоящей из большого числа Q (Q≥2n+1) разнесенных по фазе приемных элементов, при этом измерения , , суммируют с весами , , , найденными заранее, тем самым оценивают косинусную и синусную составляющие амплитуды отраженного сигнала, соответствующие центру j-го луча (центральному элементу дискретизации ДН)

затем полученные оценки возводят в квадрат, суммируют и извлекают корень, тем самым вычисляют оценки амплитуд отраженного сигнала в i-x элементах дальности и j-м синтезированном элементе разрешения азимута

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