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

Формула изобретения

Способ получения трехмерного радиолокационного изображения вращающегося по курсу, тангажу и крену объекта при многочастотном импульсном зондировании, заключающийся в излучении импульсных сигналов с несущей частотой, изменяющейся от импульса к импульсу с шагом f в полосе F, приема отраженных сигналов, измерении комплексных огибающих отраженных сигналов Ф, запоминании их в течении времени синтезирования в угловом секторе , поворота объекта по курсу, образовании матрицы комплексных огибающих и преобразовании ее с помощью быстрого двумерного преобразования Фурье в матрицу синтезированного отклика, определении величины порога по уровню первых боковых лепестков наиболее интенсивного откликом, сравнении значений элементов матрицы с порогом для выделения превышающих порог элементов матрицы, совокупность которых принимают за радиолокационное изображение объекта, отличающийся тем, что при больших секторах угла синтезирования по курсу при наличии вращения объекта по тангажу и крену, измеряют углы курса (t), тангажа (t) и крена (t) связанной системы отсчета с началом в точке М на оси вращения объекта относительно земной системы отсчета, измеряют координаты точки М и точки 0 фазового центра антенны радиолокационной станции (РЛС) в земной системе отсчета, координаты вектора преобразуют в координаты связанной системы отсчета

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

R₀ - расстояние от точки 0 фазового центра антенны РЛС до точки М, являющейся центром синтезирования, производят коррекцию фазы в измеренных значениях комплексных огибающих Ф отраженного сигнала на величину где , с - скорость света, f - частота зондирующего сигнала, определяют пространственные частоты

,

,

при значениях , , , соответствующих моменту времени измерения Ф, скорректированные по фазе значения комплексных огибающих отраженного сигнала размещают в элементы матрицы Т, трехмерные номера которых (l, m, n) определяют по формулам

где d - шаг матрицы Т по координатам f _x, f_y, f_z,

матрицу Т подвергают трехмерному дискретному преобразованию Фурье

и по полученной трехмерной матрице синтезированного отклика определяют трехмерное радиолокационное изображение объекта.

Описание изобретения к патенту

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

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

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

При зондирующем сигнале

где:S - зондирующий сигнал;

А - амплитуда сигнала;

f - несущая частота сигнала;

₀ - начальная фаза сигнала.

Комплексная огибающая отраженного сигнала есть

где:а - амплитуда отраженного сигнала,

R - расстояние от РЛС до объекта,

- длина волны зондирующего сигнала,

_отр - фаза, возникающая при отражении от объекта.

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

Известен триангуляционный способ построения двумерного РЛИ цели в РЛС сопровождения с инверсным синтезированием апертуры [1], в котором при узкополосном импульсном зондировании цели измеряют и запоминают комплексные огибающие отраженных сигналов, их подвергают одномерному преобразованию Фурье и формируют одномерный доплеровский портрет цели путем измерения доплеровской частоты i-го Рц и пересчета ее в линейное удаление Рц от линии визирования цели r в направлении линии, перпендикулярной плоскости, проходящей через линию визирования и ось относительного вращения цели при ее сопровождении, величины r определяют для двух ракурсов сопровождения цели, разделенных по углу визирования интервалом =0,5°, по данным r ₁, r ₂, зафиксированным при углах ₁, ₂, определяют координату r _d в направлении линии визирования, по координатам r , r_d, полученным для каждого i-го Рц, строят РЛИ цели.

Указанному способу присущи существенные недостатки:

координата r_d при малых углах может быть определена только с большой погрешностью,

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

имеют место трудности при идентификации образов различных Рц, полученных в портретах при разных углах визирования цели.

Известен также способ получения трехмерного изображения объекта [2], в котором измеряют в пределах полусферы объекта комплексный коэффициент рассеяния объекта являющийся отношением комплексных амплитуд падающей и отраженной волн - волновой вектор отраженной волны, - вектор нормали к поверхности объекта), определяют объединенный коэффициент отражения

где ^* - число, комплексно сопряженное к ,

и с помощью трехмерного преобразования Фурье определяют характеристическую функцию объекта

равную единице внутри поверхности объекта и нулю - вне ее, функцию подвергают преобразованию Радона для определения трехмерного изображения объекта.

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

Известен [3] способ построения двумерного радиолокационного изображения прямолинейно летящей цели при многочастотном узкополосном зондировании, основанный на излучении импульсных сигналов, приеме отраженных сигналов и накоплении их в течение рассчитываемого интервала синтезирования Т_с, отличающийся тем, что несущую частоту зондирующих импульсов изменяют от импульса к импульсу по линейному закону в диапазоне частот от f _o до f_o+ F с шагом F/2^N, частоту повторения импульсов Т_И выбирают таковой, что величина Т _и·2^N на порядок меньше времени корреляции траекторных нестабильностей полета целей (составляющих 25-100 мс), накопление отраженных импульсов в объеме 2 ^N·2^N производят в 2 ^N этапов с интервалом между этапами Т_с /2^N, составляют матрицу с 2 ^N строками и 2^N столбцами, элементами которой являются амплитуды и фазы отраженных сигналов, данные матрицы подвергают двумерному быстрому преобразованию Фурье, полученную двумерно-спектральную матрицу синтезированного отклика преобразуют в графическое матричное изображение цели, для чего определяют уровень первых боковых лепестков отклика наиболее интенсивного рассеивателя цели, принимают данный уровень за пороговую величину, сравнивают с ней значения двумерно-спектральной матрицы и в случае превышения порога выделяют соответствующий элемент в матричном поле 2^N·2 ^N, а совокупность всех выделенных элементов принимают за радиолокационное изображение цели.

Описанный способ взят в качестве прототипа.

Указанный способ обеспечивает получение двумерного портрета распределения Рц, характеризуемого интенсивностью отражения Рц и его координатами в радиальном r _d и азимутальном r направлениях. При этом координаты r _d являются результатом преобразования во временную область спектра отраженного сигнала в полосе частот F, размещенного в столбце (столбцах) матрицы. Отсчет величины r_d производится от центра синтезирования по линии визирования.

Величина r_d по существу является проекцией на линию визирования расстояния Рц до центра синтезирования, выбираемого на оси вращения объекта.

Другая координата r , получающаяся в результате преобразования Фурье строки (строк) матрицы, образованной после преобразования Фурье столбцов, является расстоянием Рц до плоскости, проходящей через линию визирования и ось вращения [4].

Недостатком способа является то, что при синтезировании РЛИ расширение сектора углов поворота объекта сверх =10-15° (с целью повышения точности спектрального разрешения Рц) невозможно из-за того, что координаты r _d, r , являются переменными величинами, изменяющимися при вращении объекта. Кроме того, если к процессу формирования РЛИ подходить с позиции метода максимального правдоподобия, используемого для решения задач обнаружения или измерения параметров сигналов, то оператор синтезирования должен являться интегралом свертки комплексной огибающей отраженного сигнала (q), зарегистрированной в координатах вектора q, изменяющегося в интервале q, и опорной (или фокусирующей) функции F(q, x),

где вектор х является совокупностью параметров, характеризующих интенсивность рассеяния и местоположения Рц.

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

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

Однако углы , являются полезными при оценке третьей координаты Рц - высоты.

Наличие движений объекта по тангажу и крену, а также учет того факта, что линия визирования РЛС - объект в реальных условиях на полигоне располагается под углом 2-3° к плоскости, перпендикулярной оси вращения объекта, создают предпосылки для получения в процессе синтезирования некоторого разрешения по высоте Рц, т.е. появляется возможность получения трехмерного РЛИ объекта.

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

Способ решает задачу построения РЛИ, в котором имеет место разрешение высоты Рц (при наличии движения объекта по тангажу и крену), т.е. фактически решается задача построение трехмерного РЛИ.

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

Кроме того, благодаря введению специальной опорной функции, осуществляющей фокусировку координат Рц, и использованию криволинейных координат f_x, f_y, f _z, реализуется процесс синтезирования при любом (0-2 ) секторе углов поворота объекта по курсу при наличии движений объекта по тангажу и крену, вследствие чего обеспечивается повышение точности оценки и разрешения координат Рц.

В предлагаемом способе оценка координат Рц производится в связанной системе отсчета, где эти координаты постоянны, и это создает дополнительные предпосылки для повышения точности оценки.

Для достижения указанного технического результата в способе получения трехмерного РЛИ, включающем излучение зондирующих импульсов с изменением несущей частоты от импульса к импульсу с шагом f в полосе частот F, прием отраженных сигналов, измерение комплексных огибающих отраженных сигналов (t) и запоминание их в течение времени синтезирования в угловом секторе , образование матрицы комплексных огибающих и преобразование ее с помощью быстрого двумерного преобразования Фурье в матрицу синтезированного отклика, определение величины порога по уровню первых боковых лепестков наиболее интенсивного отклика, сравнение значений элементов матрицы с порогом для выделения превышающих порог элементов матрицы, совокупность которых принимают за РЛИ объекта, измеряют углы курса (t), тангажа (t) и крена (t) связанной с объектом системы отсчета с началом в точке М на оси вращения объекта относительно земной системы отсчета, измеряют в земной системе отсчета координаты вектора соединяющего точку М с фазовым центром антенны РЛС (точка О), преобразуют их в координаты связанной системы отсчета

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

производят регистрацию комплексных огибающих отраженных сигналов (t) в координатах пространственных частот

при значениях f, , , соответствующих моментам измерения , корректируют фазу измеренных значений комплексных огибающих на величину получая величины и их размещают в элементы матрицы Т(l, m, n), трехмерные номера которых (l, m, n) определяют по формулам

где int() - функция взятия целой части числа,

d - шаг элементов матрицы Т по координатам f _x, f_y, f_z.

Матрицу Т подвергают трехмерному дискретному преобразованию Фурье

по полученной трехмерной матрице синтезированного отклика определяют трехмерное радиолокационное изображение объекта.

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

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

В предлагаемом способе значения размещаются в матрице в варианте двумерного РЛИ по правилу При этом значения , полученные при f=const и различных углах , будут размещены в элементы матрицы, расположенные по окружности, т.к. Кроме того, в предлагаемом способе при измерениях не требуется постоянный шаг изменения (или времени).

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

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

Процесс синтезирования РЛИ с позиции метода максимального правдоподобия выражается интегралом свертки (1), где фокусирующая функция в соответствии с [5] есть

R₀ - расстояние от фазового центра РЛС (точка О) до центра синтезирования (точка М),

R - расстояние от фазового центра РЛС до Рц, координаты которого неизвестны.

Взаимное положение точек О, М, Рц показано на фигуре 1.

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

Вектор в земной системе отсчета есть

где: x_m, y _m, z_m - координаты точки М;

x₀, y₀, z ₀ - координаты точки О.

В связанной системе отсчета вектор с координатами x_в,y _в,z_в выразится как

где A( , , ) является известной матрицей, определяющей переход из земной в связанную систему отсчета [6];

A ^T( , , ) - транспонированная матрица.

При использовании для образования углов , , последовательности правых поворотов:

- поворот вокруг оси у_земн,

- вокруг оси z₁, получившейся после поворота на угол ,

- вокруг оси x₂ (после , ).

Матрица A( , , ) выражается в виде

A^T( , , ) - транспонированная матрица.

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

Откуда:

Вектор

В связанной системе отсчета

где (x, y, z) - координаты Рц в связанной системе отсчета.

И следовательно, фокусирующая функция (13), определяемая с точностью до постоянной фазы, может быть представлена в виде:

где

f_x, f_y , f_z - пространственные частоты.

Фокусирующая функция (23) выражена через пространственные частоты f_x, f_y, f _z, являющиеся функциями и углов , , отсчитываемых в двух взаимно перпендикулярных направлениях связанной системы отсчета и представляющих собой сложные функции углов ориентации объекта , , и угла , определяющего наклон линии визирования РЛС по отношению к плоскости, перпендикулярной оси вращения объекта по углу :

Для практического использования оператора синтезирования (1) комплексные огибающие отраженного сигнала должны быть зафиксированы в координатах f_x, f_y, f _z.

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

И тогда оператор синтезирования (1) выразится как

Измерения производят при дискретных значениях f, , , . Вследствие этого преобразование (5) должно быть представлено в дискретной форме. Для этого переменные интегрирования f _x, f_y, f_z заменяют номерами элементов трехмерной матрицы Т

где d - шаг построения элементов матрицы.

Поскольку диапазон изменения координат f_x , f_y, f_z одинаков и равен величина d должна быть одинаковой для f _x, f_y, f_z.

Величина d выбирается из условия обеспечения размещения величин в различные элементы матрицы при различных параметрах f, , , . При этом величина d будет зависеть от того, насколько изменяются f, , , от измерения к измерению.

На практике выбор величины d целесообразно проводить методом проб: задают размер матрицы N³, размещают измерения в элементы матрицы и фиксируют события повторного размещения в элементы матрицы. Если повторы имели место, то размер матрицы увеличивают до (2N)³ и т.д.

Преобразование (5) в дискретной форме выражается суммой

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

При N=2^k реализация (6) выполняется с помощью быстрого дискретного преобразования Фурье.

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

Действительно, при , 0 матрица A( , , ) приобретает вид

И тогда:

Обозначив получим

Угол ₀ определяет начальную ориентацию линии визирования объекта относительно земной системы отчета.

При малых значениях угла : sin 0, cos 1

Оператор синтезирования приобретает вид:

Приведение этого оператора к дискретному двумерному преобразованию Фурье аналогично описанному для случая трехмерного варианта.

Построение РЛИ, основанное на использовании предлагаемого способа, проведено методом моделирования.

Задавалась модель объекта в виде совокупности неподвижных относительно связанной системы отсчета Р_ц. Зондирующие сигналы РЛС - импульсы с периодом повторения 20 мкс. Несущая частота сигнала изменялась от импульса к импульсу с шагом 30 мГц в полосе от 2500 до 3500 МГц. Объект вращался по курсу со скоростью 12°/с. Линия визирования (РЛС - объект)составляла =3° и 20° по отношению к плоскости, перпендикулярной оси вращения объекта.

Сектор синтезирования =90°. Движение объекта по тангажу задавалось в виде гармонического колебания с амплитудой 5° и количеством периодов колебаний 2. Количество заданных Рц - 7. Амплитуды отраженных сигналов полагались одинаковыми для каждого Рц, равными 1. На фигурах 2 и 3 приводится трехмерное изображение объекта с проекциями синтезированных откликов Рц на координатные плоскости при углах визирования =3° и =20°. При этом сферами показаны положения рассеивающих центров с истинными координатами - в центре сферы. Точками в районе сферы отмечены области, где амплитуда сигнала в синтезированном отклике превышает половину максимального значения.

Значения амплитуд отклика на трехмерном изображении не показаны.

Полученные в откликах значения амплитуд были близкими к единице.

Параметры, при которых получено изображение:

сектор угла курса : 90°,

амплитуда тангажа: А =5°,

количество периодов колебаний тангажа: 2,

вертикальный угол визирования объекта =3°.

На фигурах 4 и 5 приводится двумерное РЛИ объекта при секторе угла : 10° и двумерное РЛИ объекта при секторе угла : 90°.

Здесь на плоскости (x, z) показаны координаты синтезированных откликов Рц, а в направлении, перпендикулярном плоскости (x, z), отложены амплитуды откликов.

Литература

1. Триангуляционный способ построения двумерного РЛИ цели в РЛС сопровождения с инверсным синтезированием апертуры. Патент RU 2099742 С1, кл. G 01 S 13/89.

2. Астанин Л.Ю., Костылев А.А. Основы сверхширокополосных радиолокационных измерений. М.: Радио и связь. 1989.

3. Способ построения двумерного РЛИ прямолинейно летящей цели при многочастотном узкополосном зондировании. Патент RU 2099743 С1, кл. G 01 S 13/89.

4. Зиновьев Ю.С., Пасмуров А.Я. Методы обращенного синтезирования апертуры в радиолокации с помощью узкополосных сигналов. Зарубежная радиоэлектроника. №3, 1985.

5. I.Fortuny. An Efficient 3-D Near-Field ISAR Algorithm. IEEE Transactions on Aerospace and Electronic Systems. Vol.34, №4, October 1998.

6. Г.Корн, Т.Корн. Справочник по математике. М.: Наука, 1973.