Патент на изобретение №2372627

(54) СПОСОБ ПОЛУЧЕНИЯ ДВУМЕРНОГО РАДИОЛОКАЦИОННОГО ИЗОБРАЖЕНИЯ ОБЪЕКТА В БОЛЬШОМ ДИАПАЗОНЕ ИЗМЕНЕНИЯ ВЕЛИЧИН ЭФФЕКТИВНЫХ ПЛОЩАДЕЙ РАССЕИВАНИЯ ЛОКАЛЬНЫХ ЦЕНТРОВ ПРИ МНОГОЧАСТОТНОМ ИМПУЛЬСНОМ ЗОНДИРОВАНИИ

(57) Реферат:

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

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

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

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

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

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

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

_отр – скачок фазы, возникающий при отражении от объекта.

Модель исследуемого объекта – совокупность локальных рассеивающих центров (РЦ), у которых эффективная площадь рассеивания (ЭПР) при величине , изменяющейся в относительно небольшой полосе ||0,15 _ср, где _ср – средняя в полосе перестройки длина волны, а также координаты в связанной с объектом системе отсчета являются постоянными.

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

r, r_d, полученным для каждого i-го РЦ, строят РЛИ цели.

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

координата r_d при малых секторах (₂–₁) может быть определена только с большой погрешностью; при узкополосном зондировании разрешение РЦ является плохим; имеют место значительные трудности при идентификации образов РЦ, полученных в портретах при различных углах визирования цели. Что же касается важного для практики параметра – величины ЭПР отдельных РЦ, то в рассмотренном способе оценка его при малых секторах углов (₂–₁) является грубой.

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

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

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

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

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

Расширение этого сектора теоретически должно приводить к повышению разрешающей способности и точности оценок [4]. Однако в рассматриваемом способе этого не происходит из-за того, что координаты r, r_d изменяются при вращении объекта, а процесс синтезирования дает их усредненную оценку. Изображение РЦ расфокусируется, падает разрешающая способность – вследствие этого способ не может обеспечить высокой точности определения координат и ЭПР РЦ, особенно в большом диапазоне изменения ЭПР РЦ.

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

Способ решает задачу получения двумерного РЛИ объекта, обеспечивающего повышение разрешающей способности и точности определения координат и ЭПР РЦ в широком диапазоне изменения их ЭПР.

Требования к диапазону разрешаемых ЭПР РЦ в настоящее время составляет 20 дБ и более.

Технический результат, состоящий в повышении разрешающей способности в широком диапазоне изменения ЭПР РЦ и повышении точности оценки координат и ЭПР РЦ, достигается благодаря измерению с повышенной точностью частоты f сигнала и угла курса объекта, являющихся параметрами процесса синтезирования РЛИ, в значительной мере определяющими качество РЛИ; реализации процесса синтезирования в виде свертки функции комплексных огибающих отраженных сигналов со специальной опорной функцией для повышения фокусировки РЛИ; использовании криволинейных координат f_x, f_z для устранения нелинейного влияния угла на процесс синтезирования; оценивании в процессе синтезирования координат РЦ в связной системе отсчета (где они постоянны) для повышения точности оценивания; уточнении полученных из РЛИ оценок координат РЦ за счет использования информации об энергии комплексных огибающих; использовании уточненных координат РЦ и полученных из РЛИ оценок ЭПР РЦ для коррекции матрицы комплексных огибающих с целью изъятия из нее информации о группе РЦ с большими ЭПР, что позволяет перенести диапазон разрешаемых при синтезировании РЦ в сторону РЦ с малыми ЭПР; повторной реализации процесса синтезирования со скорректированной матрицей комплексных огибающих для получения РЛИ с малыми значениями ЭПР РЦ.

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

измеряют частоту f(t_nm) зондирующих импульсов, момент времени t_nm запоминают, где n – номер шага перестройки частоты, m – номер цикла перестройки (номер повторной перестройки), измеряют в земной системе отсчета координаты фазового центра антенны РЛС x₀, у₀, z₀ и координаты x_м(t_nm), y_м(t_nm), z_м(t_nm) точки М на объекте, выбранной в качестве центра синтезирования, измеряют относительно земной системы отсчета угол курса (t_nm) связанной с объектом системы отсчета с началом в точке М, вычисляют пространственные частоты

f_x(t_nm)=-2f(t_nm)/c·sin((t_nm)-₀),

где с – скорость света,

₀ – угол между линией, соединяющей фазовый центр антенны РЛС с центром синтезирования, и осью z земной системы отсчета,

фазы измеренных значений комплексных огибающих Ф(t_nm) отраженных сигналов корректируют на величину 4f(t_nm)/c·R₀(t_nm),

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

полученные после коррекции значения комплексных огибающих Ф^*(t_nm) размещают в элементы матрицы T(k,l), двумерные номера (k,l), которых определяют по формулам

где int [] – функция определения целой части числа,

d – шаг элементов матрицы Т по координатам f_x, f_z,

матрицу T(k,l) подвергают двумерному дискретному преобразованию Фурье

и по полученной двумерной матрице синтезированных откликов g(x,z) определяют двумерное радиолокационное изображение объекта, задают величину порога разрешения по величине ЭПР рассеивающих центров и для элементов матрицы g(x,z), у которых квадрат модуля |g(x,z)|² превышает порог , определяют оценки координат x_i^*, z_i^*, и оценки ЭПР _i^* выделенных РЦ, координаты x_i^*, z_i^* уточняют путем поиска в окрестности каждого выделенного РЦ координат x_i^**, z_i^**, дающих минимум по переменным х, z функции энергии

где

R(x,z) – расстояние от фазового центра антенны РЛС до точки объекта с координатами x, z из окрестности выделенного i-го РЦ,

и по полученным оценкам ЭПР i^* и уточненным значениям координат x_i^**, z_i^** выделенных РЦ определяют корректирующие поправки для комплексных огибающих

скорректированную матрицу комплексных огибающих со значениями ее элементов

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

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

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

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

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

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

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

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

Оператор синтезирования РЛИ, являющийся инструментом оценки неизвестных параметров (координат и ЭПР РЦ), при использовании метода максимального правдоподобия должен представлять собой интеграл свертки комплексных огибающих Ф(q), зарегистрированных в координатах вектора q, изменяющегося в интервале q, с опорной функцией F(q, x):

где вектор x является совокупностью неизвестных координат и ЭПР РЦ.

Опорная функция в соответствии с [5] является комплексной огибающей отраженного сигнала фиктивного РЦ, координаты x, y, z которого в связной системе отсчета произвольны:

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

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

где Х_з, Z_з и Х_св, Z_св – оси координат земной и связанной систем отсчета.

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

где

где ₀ – угол ориентации линии визирования (линии, соединяющей фазовый центр антенны РЛС с центром синтезирования на объекте) относительно оси z земной системы отсчета.

При R₀/Rl соотношение (10) преобразуется к виду:

где f_x=-2/sin(–₀), f_z=2/cos(–₀) – пространственные частоты.

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

выражается в виде:

где Ф^*(f_x, f_z) – значения Ф(t_nm), размещенные в области переменных f_x, f_z.

В выражении (16) третий сомножитель в процессе синтезирования (при изменении , ) изменяется мало.

Действительно, при R₀=900 м, =3 см в диапазоне изменения , равном ±0,15, для РЦ, удаленных от центра синтезирования на 1 м (т.е. при (х²+у²+z²)^0,5=1 м), величина

₃=4/·(х²+у²+z²)/(2R₀) изменяется не более чем на ±2° от ее среднего значения (при =0).

А для РЦ с удалением от центра синтезирования на 3 м изменение ₃ составляет ±18°.

Поскольку ошибки измерения фаз комплексных огибающих Ф^* на практике укладываются в интервале ±30°, то влиянием на процесс синтезирования третьего сомножителя в (16) можно пренебречь при указанных выше значениях дальности РЛС – объект R₀ и размера объекта (х²+у²+z²)^0,5.

И тогда оператор синтезирования РЛИ (16) становится двумерным преобразованием Фурье

Измерения Ф^* производят при дискретных значениях f, , вследствие этого преобразование (17) должно быть представлено в дискретной форме. Для этого переменные интегрирования f_x, f_z заменяют номерами элементов двумерной матрицы Т

где int [] – функция определения целой части числа,

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

Поскольку диапазон изменения координат f_x, f_z одинаков и равен ±2f_max/c, где f_max – максимальное значение частоты в полосе F ее перестройки, то величина d может быть выбрана одинаковой для f_x и f_z.

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

Опыт моделирования показывает, что размерность К матрицы Т в области переменных f_x, f_z является в 2÷4 раза большим, чем в области переменных f, .

При дискретизации переменных f_x, f_z преобразование (17) в дискретной форме выражается в виде суммы

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

Если размерность матрицы Т удовлетворяет условию К=2^P, где p – целое число, то реализация (19) выполняется с помощью быстрого двумерного дискретного преобразования Фурье.

После реализации алгоритма синтезирования (19) в полученном РЛИ выделенные по превышению порога элементы матрицы синтезированного отклика содержат информацию о координатах x, z и об ЭПР рассеивающих центров .

Координаты выделенных РЦ определяются по двумерному номеру отклика в матрице синтезированных откликов, а величины ЭПР РЦ – по квадрату модуля отклика, который применительно к способу-прототипу пропорционален величине , полосе перестройки зондирующих сигналов F и интервалу времени синтезирования [6].

Моделирование показывает, что при достижимом в настоящее время уровне значений ошибок измерения амплитуды (/A0,5÷1 дБ) и фазы (30°) комплексных огибающих отраженных от объекта сигналов в диапазоне перестройки частоты зондирующих импульсов от 8500 до 11500 МГц, ошибок измерения координат центра синтезирования _х, _z2 мм и угла курса <3′ синтезирование в секторе =16° обеспечивает получение оценок координат x, z и ЭПР РЦ с ошибками

_х, _z5 см, A/А=1 дБ.

Однако это выполняется только для тех РЦ, чьи величины ЭПР находятся в диапазоне 0÷-12 дБ от уровня максимального значения ЭПР этих РЦ.

Вероятность правильного выделения (идентификации) РЦ с малыми значениями ЭПР (-12÷-20 дБ) падает от 1 при =-12 дБ до 0 при =-20 дБ.

Для расширения диапазона ЭПР правильно идентифицируемых РЦ и нахождения оценок их координат и ЭПР в предлагаемом способе информация об идентифицированных в процессе синтезирования РЦ с большими ЭПР (в диапазоне 0÷-12 дБ) изымается из исходной матрицы Т комплексных огибающих и производится повторное синтезирование с использованием скорректированной матрицы комплексных огибающих.

Однако эффективность такой коррекции существенно зависит от точности определения оценок координат РЦ.

Моделирование показывает, что используемое здесь преобразование Фурье комплексных огибающих Ф^* не реализует потенциальную точность, определяемую информацией матрицы Т. По крайней мере, точность оценок координат x, z идентифицированных (выделенных) РЦ может быть увеличена более, чем на порядок.

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

Задается порог разрешения по ЭПР для РЦ с большими значениями ЭПР. Выделяется РЦ, чьи ЭПР превысили порог, определяются координаты выделенных РЦ x_i, z_i и их ЭПР _i. Пусть количество выделенных РЦ составило I, для каждого выделенного i-го РЦ () производится уточнение координат. В окрестности точки x_i, z_i создается область G_i(x, z) с центром в точке x_i, z_i и радиусом =5÷10 см. В области G_i задается множество точек x, z с удалением друг от друга на расстоянии, меньшем ·10^-2.

Каждая из точек x, z области G проверяется на возможность нахождения в ней выделенного РЦ.

Для этого определяются комплексные огибающие Ф_i(t_nm), создаваемые фиктивным РЦ с координатами x, z и ЭПР _i для имеющегося множества измеренных значений f(t_nm), (t_nm):

где R(x, z) – дальность от фазового центра антенны РЛС до точки (x, z) области G.

Комплексные огибающие отраженных сигналов от объекта, в котором удален РЦ с параметрами x, z, _i, находятся как

Энергия процесса, представленного всей совокупностью комплексных огибающих Ф_i(t_nm) при n= и m=, есть:

Величина Э_i(, x, z) находится для всех точек x, z области G_i(x, z). Точка x, z, в которой достигается min Э_i(, x, z) по переменным x, z, берется в качестве уточненной оценки местоположения выделенного i-го РЦ:

Реализация описанной процедуры позволяет уменьшить ошибки определения координат выделенных РЦ до 1-2 мм.

Установленные координаты i-ых выделенных РЦ xi, z_i и оценки их ЭПР _i используется для удаления информации о выделенных РЦ из матрицы комплексных огибающих.

Действительно, вклад i-го РЦ с координатами x_i, z_i и ЭПР _i в комплексную огибающую отраженного сигнала на частоте зондирования f_n в момент времени t_nm есть:

где R(x_i, z_i) – расстояние от фазового центра антенны РЛС до точки объекта с координатами x_i, z_i в момент времени t_nm.

Это выражение в соответствии с (14) преобразуется к виду

Если известны оценки x_i, z_i, _i для , то при изъятии из матрицы Т информации об i-х РЦ элементы матрицы Т корректируется на величину

Скорректированные значения комплексных огибающих есть:

Выражения (26), (27) в области переменных k, 1 запишутся как

Реализация процесса синтезирования РЛИ по данным скорректированных значений комплексных огибающих (29) производится с помощью преобразования (19).

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

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

Задавалась модель объекта в виде совокупности неподвижных относительно связанной системы отсчета РЦ.

Зондирующие сигналы РЛС – импульсы с периодом повторения 20 мкс. Несущая частота сигнала изменялась от импульса к импульсу с шагом f=3000/511 МГц в полосе частот от 8500 до 11500 МГц.

Объект вращался по курсу со скоростью 12°/с.

Сектор синтезирования =16°.

РЦ в количестве 30 располагались тремя изолированными группами с удалением соседних РЦ на 40 см.

Величины ЭПР РЦ (или амплитуды отраженных от РЦ сигналов, равные ) задавались равными 0,1 м²; 0,01 м (-10 дБ относительно 0,1 м²) и 0,001 м ² (-20 дБ).

В каждой группе имелись РЦ с большими, средними и малыми значениями ЭПР.

На фиг.2 приведено двумерное РЛИ, полученное при использовании исходной (некомпенсированной) матрицы комплексных огибающих отраженных сигналов, а на фиг.3 – двумерное РЛИ, полученное после компенсации информации о РЦ с ЭПР, имеющими 0 и – 10 дБ. Максимальные значения ошибок оценок координат и ЭПР компенсируемых РЦ составляли х, z2 мм, A/А=1 дБ.

В РЛИ на фиг.2 РЦ с малыми значениями ЭПР (=-20 дБ) теряются на фоне изображения.

В РЛИ на фиг.3 указанные РЦ становятся доступными для идентификации.

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

Использованные источники информации

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

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

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

4. Кондратенков Г.С., Фролов А.Ю. Радиовидение. М.: Радиотехника, 2005.

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

6. Митрофанов Д.Г. Формирование двумерного радиолокационного изображения цели с траекторными нестабильностями полета. Радиотехника и электроника. Том 47, 7, 2002.

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

Способ получения двумерного радиолокационного изображения объекта в большом диапазоне изменения величин эффективных площадей рассеивания локальных рассеивающих центров при многочастотном импульсном зондировании, включающий излучение импульсов с изменением несущей частоты от импульса к импульсу с шагом f полосе частот F, прием отраженных сигналов, измерение комплексных огибающих Ф отраженных сигналов и запоминание их в течение времени синтезирования в угловом секторе , образование матрицы комплексных огибающих и преобразование ее с помощью быстрого двухмерного преобразования Фурье в матрицу синтезированных откликов, определение величины порога по уровню первых боковых лепестков наиболее интенсивного отклика, сравнение величин откликов с порогом для выделения превышающих порог элементов матрицы, совокупность которых принимают за радиолокационное изображение объекта, в котором по положению выделенных откликов в матрице синтезированных откликов определяют значения координат, а по величинам квадратов модулей откликов определяют значения эффективных площадей рассеивания рассеивающих центров, отличающийся тем, что измеряют частоту f(t_nm) зондирующих импульсов, момент времени измерения t_nm запоминают, где n – номер шага перестройки частоты в полосе частот F, m – номер повторного цикла перестройки, измеряют в земной системе отсчета координаты фазового центра антенны радиолокационной станции (РЛС) x₀, y₀, z₀ и координаты x_м(t_nm), y_м(t_nm), z_м(t_nm) точки М на объекте, выбранной в качестве центра синтезирования, измеряют относительно земной системы отсчета угол курса
(t_nm) связанной с объектом системы отсчета с началом в точке М, вычисляют пространственные частоты
f_x(t_nm)=-2f(t_nm)/c·sin((t_nm)-₀),
f_z(t_nm)=2f(t_nm)/c·cos((t_nm)-₀),
где с – скорость света,
₀ – угол между линией, соединяющей фазовый центр антенны РЛС с центром синтезирования, и осью z земной системы отсчета, фазы измеренных значений комплексных огибающих Ф(t_nm) отраженных сигналов корректируют на величину (4f(t_nm))/c·R₀(t_nm), где R₀(t_nm) – расстояние от фазового центра антенны РЛС до точки синтезирования М, полученные после коррекции значения комплексных огибающих Ф*(t_nm) размещают в элементы матрицы T(k, l), двумерные номера которых (k, l) определяют по формулам
k=int[f_x(t_nm)/d],
l=int[f_z(t_nm)/d],
где int [] – функция определения целой части числа,
d – шаг элементов матрицы Т по координатам f_x, f_z, матрицу T(k, l) подвергают двумерному дискретному преобразованию Фурье
g(x,z)=Ф*(k, l)exp[j2(kx+lz)]
и по полученной двумерной матрице синтезированных откликов g(x, z) определяют двумерное радиолокационное изображение объекта, задают величину порога разрешения по величине эффективной площади рассеяния рассеивающих центров и для элементов матрицы g(x, z), у которых квадрат модуля |g(x,z)|² превышает порог , определяют оценки координат x_i^*, z_i^*, и оценки эффективных площадей рассеивания _i* выделенных рассеивающих центров, координаты x_i*, z_i* уточняют путем поиска в окрестности каждого выделенного рассеивающего центра координат x_i**, z_i**, дающих минимум по переменным x, z функции энергии

где
R(x, z) – расстояние от фазового центра антенны РЛС до точки объекта с координатами x, z из окрестности выделенного i-го рассеивающего центра, и по полученным оценкам эффективных площадей рассеивания _i* и уточненным значениям координат х_i^**, z_i^** выделенных рассеивающих центров определяют корректирующие поправки для комплексных огибающих

скорректированную матрицу комплексных огибающих со значениями ее элементов
Ф^**(k,l)=Ф^*(k,l)-Ф_кор(k,l)
подвергают двумерному дискретному преобразованию Фурье
g^*(x,z)=Ф^**(k,l)ехр[j2(kx+lz)]
и по полученной двумерной матрице g^*(x,z) синтезированных откликов определяют двумерное радиолокационное изображение объекта для рассеивающих центров с величинами эффективных площадей рассеивания, меньшими, чем заданный порог .

РИСУНКИ