СПОСОБ ФОРМИРОВАНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ЗЕМНОЙ ПОВЕРХНОСТИ И ВОЗДУШНОЙ ОБСТАНОВКИ С ПОМОЩЬЮ АНТЕННОЙ РЕШЕТКИ Российский патент 2015 года по МПК G01S13/89 

Описание патента на изобретение RU2539558C1

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

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

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

Наиболее близким по технической сущности является способ повышения разрешающей способности фазированной антенной решетки (ФАР) бортовой станции [4], заключающийся в формировании ДН N-канальной ФАР с использованием рассчитанных комплексных весовых коэффициентов, отличающийся тем, что при приеме отраженного от земной поверхности зондирующего сигнала с направления θ, φ угла места θ и азимута φ комплексные сигналы y ˙ k , принятые в k-x каналах ФАР ( k = 1 , N ¯ , N - число каналов), суммируют с весовыми коэффициентами рассчитанными методом восстановления n угловых составляющих поля отражения от земной поверхности в пределах ширины ДН ФАР (n<N), и осуществляют оценивание j-х комплексных амплитуд поля отражения по формуле

затем берут модули оценок и получают амплитудное изображение элементов поверхности, наблюдаемых в пределах ширины ДН с повышенным в n раз угловым разрешением.

Суммирование в формуле (1) ведется по всем j-м элементам линейной или плоской ФАР. Для плоской ФАР, ДН которой обладает свойством разделения по угловым переменным, формула (1), так же как в [3], раскрывается в виде экономичной двухэтапной процедуры (восстановление матрицы изображения вначале по i-м строкам и затем по j-м столбцам):

где m и n - число элементов дискретизации по угловым координатам θ и φ (m<M, n<N); y ˙ q , k - комплексные измерения в q,k-x каналах ФАР и - элементы матриц весовых коэффициентов WA и WB.

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

1. Для рассмотренной в [4] системы координат ДН ФАР не разделяется по угловым переменным (азимуту и углу места). Следовательно, заявленная двухэтапная процедура в такой системе координат не применима.

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

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

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

Технический результат предлагаемого технического решения достигается тем, что способ формирования трехмерного изображения земной поверхности и воздушной обстановки с помощью антенной решетки в виде M×N-матрицы приемных элементов заключается в том, что при наблюдении за земной поверхностью или воздушной обстановкой с помощью сканирующего по углам луча антенны в зоне обзора измеряют при каждом положении луча в q,k-x измерительных каналах значения комплексных амплитуд принимаемых сигналов y ˙ q , k , формируют из этих амплитуд M×N-матрицу измерений Y с элементами и далее обрабатывают,

отличающийся тем, что при сканировании последовательно смещают луч на ширину диаграммы направленности антенны (на уровне 0,5 мощности) построчно в зоне обзора, комплексные амплитуды измеряют в каждом l-м элементе разрешения дальности ( l = 1 , L ¯ , L - число таких элементов), запоминают их в M×N-матрице измерений Y и обрабатывают в антенной системе координат, при этом умножают матрицу Y слева на рассчитанную заранее m×M-матрицу весовых коэффициентов WA и получают m×N-матрицу S восстановленного по строкам изображения: S=WAY, далее матрицу S умножают справа на рассчитанную заранее N×n-матрицу весовых коэффициентов WB и получают m×n-матрицу U восстановленного по строкам и столбцам изображения с повышенным в mn раз (m<N, n<N) угловым разрешением: U=SWB, затем вычисляют модули комплексных элементов матрицы повторяют операции для всех l-х элементов дальности и запоминают координаты тех i,j,l-x пространственных элементов, амплитуды которых превышают порог γ0 обнаружения полезного сигнала: u(i,j,l)>γ0, в результате получают трехмерное изображение объектов отражения в виде совокупности их прямоугольных координат xj,yj,zl в зоне видимости антенны, затем повторяют операции для всех положений луча и получают трехмерное изображение в полной зоне обзора.

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

1. Луч антенны с помощью электронного управления сканирует зону обзора, смещаясь по азимуту и углу места построчно на ширину ДН (на уровне 0,5 мощности).

2. После посылки зондирующих импульсов отраженные сигналы принимаются q,k-ми элементами АР (одиночными или объединенными в группы), которые расположены в плоскости АР в матричном порядке по q-м строкам и k-м столбцам (, М и N - число элементов по строке и столбцу).

3. Принятые сигналы поступают в q, k-e приемные каналы , проходят в них тракт первичной обработки (стробирование по дальности, фазовое детектирование, фильтрацию на низкой частоте [5]), в результате которой на видеочастоте выделяются комплексные амплитуды принимаемых сигналов в каждом l-м элементе разрешения по дальности ( l = 1 , N ¯ , L - число таких элементов).

4. Комплексные амплитуды помещаются в М×N-матрицу измерений Y в каждом l-м элементе дальности. Дальнейшая обработка осуществляется в l-х элементах дальности независимо и параллельно в антенной системе угловых координат φx, θy.

5. Матрица Y умножается слева на рассчитанную заранее m×М-матрицу весовых коэффициентов WA. Получается m×N-матрица S восстановленного по строкам изображения: S=WAY.

6. Матрица S умножается справа на рассчитанную заранее N×n-матрицу весовых коэффициентов WB. Получается m×n-матрица U восстановленного по строкам и столбцам изображения с повышенным в mn раз (m<N, n<N) угловым разрешением: U=SWB.

7. Вычисляются модули комплексных элементов матрицы и запоминаются номера i,j амплитуд u(i,j,l), которые превышают порог обнаружения полезного сигнала γ0: u(i,j,l)>γ0 в l-м элементе дальности.

8. На множестве l-х элементов дальности для запомненных номеров i,j,l, имеющих угловые координаты θy(i), φx(j) и дальность r(l), вычисляются прямоугольные координаты x=φx(j)r(l), y=θy(i)r(l), z=r(l), совокупность которых представляет трехмерное изображение объектов отражения в зоне видимости АР.

9. Для всех положений сканирующего луча антенны получается трехмерное изображение в виде совокупности точек с координатами x, y, z в полной зоне обзора.

10. Полученные координаты пересчитываются в глобальную систему координат с целью пространственно-временной обработки изображений.

Расчетная часть

Наблюдение земной поверхности или воздушной обстановки ведется в антенной прямоугольной системе координат o,x,y,z, ось oz совпадает с центральным направлением излучения и приема отраженных сигналов. В антенной системе также используется система координат r,φxy, где r - наклонная дальность; φx - угол между осью oz и геометрической проекцией луча РЛС на горизонтальную плоскость o,x,z; θy - угол между осью oz и геометрической проекцией луча на вертикальную плоскость o,y,z.

На фигуре 1 показано следующее: заштрихованная часть - плоскость антенны; o,x,y,z - прямоугольные антенные координаты; φ, θ - сферические угловые координаты (азимут и угол места); φx=φ - азимут и θy - угол места в проекции на координатную плоскость o,y,z; - вектор луча отражения от точки M; α - угол между вектором и вектором скорости движения носителя РЛС.

Антенна выполнена в форме плоской М×N-решетки с q,k-ми элементами, расположенными в матричном порядке работающими на излучение и прием отраженных сигналов. Центры приемных элементов антенны расположены в точках М(xk,yq) с координатами xk, yq.

Комплексная модель сигналов представлена двойной суммой:

где

где - комплексная амплитуда принимаемого сигнала в l-м элементе разрешения дальности с углового направления θy(i), φx(j); множитель i под знаком экспоненты - мнимая единица; ξl - изменение фазы сигнала при приеме в l-м элементе дальности; u(i,j,l) - амплитуда сигнала отражения в i,j-м угловом направлении и l-м элементе дальности, характеризующая отражающую способность i,j,l-го пространственного элемента разрешения (ненулевая для элемента отражения); ε - неизвестное изменение фазы в l-м элементе дальности, которое в l-м элементе можно считать неизменным, а на множестве значений l - случайным, равномерно распределенным на [0,2π]; λ - длина волны; - комплексная нормированная ДН q,k-го приемного элемента АР, амплитудная часть которой g(i,j) - известная функция θy(i) и φx(j) (µ - известный коэффициент); Δ - ширина круговой ДН антенны; δq,kq,k(i,j) - известное опережение или запаздывание по фазе сигнала, принимаемого с i,j-го углового направления в q,k-м приемном элементе по сравнению с центром антенны; - комплексный гауссовский белый шум с нулевым математическим ожиданием и дисперсией действительной и мнимой части.

Суммирование в (3) ведется по i,j-м элементам дискретизации угломерной области {(θ,φ)} в пределах ДНА.

Расчет δq,kq,k(i,j) в угловых антенных координатах θy(i), φx(j) сводится к следующему. Для плоского фронта волны вектор нормали или

По формуле расстояния от точки M(xq,yk) центра q,k-го элемента антенны до плоскости волны, проходящей через центр антенны, получаем

Для ширины ДН Δ порядка 10-20° (в радианах это малая величина по сравнению с единицей), формулу (4) можно упростить:

Тогда комплексная ДН каждого q,k-го элемента АР разделяется по угловым переменным φx и θy следующим образом:

где

С учетом (6) модель (3) принимает вид повторной суммы:

или в матричной форме:

где Y - М×N-матрица измерений

А - М×m-матрица коэффициентов

U - m×n-матрица искомых амплитуд;

В - m×N-матрица коэффициентов записанная в транспонированном виде по отношению к матрице A; P - М×N-матрица помех

Оптимальную оценку матрицы U находим по критерию минимума следа матрицы (Y-AUB)H(Y-AUB) по формуле [3, 4]:

где H - символ операции транспонирования и сопряжения; I - единичная матрица; δ - параметр регуляризации.

Операции (9) удобно выполнять в виде двухэтапной процедуры:

где - матрицы комплексных весовых коэффициентов. Для круговой ДН (m=n): В=AH.

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

Модули элементов комплексной матрицы оценок , превышающие порог обнаружения полезного сигнала, по совокупности значений i, j дают амплитудное изображение объектов отражения в l-м элементе дальности в зоне видимости АР. На множестве значений l получаем трехмерное изображение поверхности в виде совокупности точек с координатами x=φx(j)r(l), y=θy(i)r(l), z=r(l) в антенной прямоугольной системе.

Точность оценок амплитуд при заданном числе m, n элементов дискретизации угломерной области {(θ,φ)} определяется числом и взаимным расположением элементов антенной решетки и соответственно структурой матриц A и B. Ковариационную матрицу Coν ошибок оценивания на каждом этапе оценивания (10) приближенно можно определить по формуле . Наоборот, при заданном предельном значении СКО ошибки оценивания находятся числа m и n, определяющие угловую разрешающую способность трехмерного изображения.

Литература

1. Монзинго Р.А., Миллер Т.У. Адаптивные антенные решетки: Введение в теорию / пер. с англ. - М.: Радио и связь, 1986. 448 с.

2. Активные фазированные антенные решетки / Под ред. Д.И. Воскресенского и А.И. Канащенкова.- М.: Радиотехника. 2004. 487 с.

3. Патент RU 2379705 C2.

4. Патент RU 2464680 C1.

5. Радиолокационные станции с цифровым синтезированием апертуры антенны / А.Н. Антипов, В.Т. Горяинов, А.Н. Кулин и др.; под ред. В.Т. Горяинова.- М.: Радио и связь, 1988. 304 с.

Похожие патенты RU2539558C1

название год авторы номер документа
СПОСОБ ИЗМЕРЕНИЯ КООРДИНАТ ЭЛЕМЕНТОВ ЗЕМНОЙ ПОВЕРХНОСТИ В БОРТОВОЙ ЧЕТЫРЕХКАНАЛЬНОЙ ДОПЛЕРОВСКОЙ РЛС 2013
  • Клочко Владимир Константинович
  • Нгуен Чунг Тхык
RU2534224C1
СПОСОБ ФОРМИРОВАНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ПОВЕРХНОСТИ НА БАЗЕ БОРТОВОГО РАДИОТЕПЛОЛОКАТОРА 2008
RU2368918C1
СПОСОБ НАБЛЮДЕНИЯ ЗА ПОВЕРХНОСТЬЮ И ВОЗДУШНОЙ ОБСТАНОВКОЙ НА БАЗЕ МНОГОКАНАЛЬНОЙ БОРТОВОЙ РЛС 2006
  • Клочко Владимир Константинович
RU2316787C1
СПОСОБ ФОРМИРОВАНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ЗЕМНОЙ ПОВЕРХНОСТИ В БОРТОВОЙ ДОПЛЕРОВСКОЙ РЛС С ЛИНЕЙНОЙ АНТЕННОЙ РЕШЕТКОЙ 2014
  • Клочко Владимир Константинович
RU2569843C1
СПОСОБ НАБЛЮДЕНИЯ ЗА ПОВЕРХНОСТЬЮ НА БАЗЕ МНОГОКАНАЛЬНОЙ БОРТОВОЙ РЛС 2006
  • Клочко Владимир Константинович
RU2316786C1
СПОСОБ ФОРМИРОВАНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ЗЕМНОЙ ПОВЕРХНОСТИ В БОРТОВОЙ ЧЕТЫРЕХКАНАЛЬНОЙ ДОПЛЕРОВСКОЙ РЛС 2014
  • Клочко Владимир Константинович
RU2572357C1
СПОСОБ ПОВЫШЕНИЯ РАЗРЕШАЮЩЕЙ СПОСОБНОСТИ ФАЗИРОВАННОЙ АНТЕННОЙ РЕШЕТКИ БОРТОВОЙ СТАНЦИИ 2011
  • Клочко Владимир Константинович
RU2464680C1
СПОСОБ НАБЛЮДЕНИЯ ЗА ПОВЕРХНОСТЬЮ И ВОЗДУШНОЙ ОБСТАНОВКОЙ НА БАЗЕ БОРТОВОЙ РЛС 2005
  • Клочко Владимир Константинович
RU2284548C1
СПОСОБ ИЗМЕРЕНИЯ ВЫСОТЫ ОБЪЕКТОВ НА ПОВЕРХНОСТИ ПРИ ПОЛУЧЕНИИ ТРЕХМЕРНОГО РАДИОЛОКАЦИОННОГО ИЗОБРАЖЕНИЯ ПОВЕРХНОСТИ С ОБЪЕКТАМИ НА БАЗЕ БОРТОВОЙ РЛС МАЛОВЫСОТНОГО ПОЛЕТА 2005
  • Клочко Владимир Константинович
RU2300780C1
СПОСОБ ПОЛУЧЕНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ПОВЕРХНОСТИ ПО ДАННЫМ БОРТОВОЙ ИМПУЛЬСНО-ДОПЛЕРОВСКОЙ РЛС МАЛОВЫСОТНОГО ПОЛЕТА 2005
  • Клочко Владимир Константинович
RU2299448C2

Реферат патента 2015 года СПОСОБ ФОРМИРОВАНИЯ ТРЕХМЕРНОГО ИЗОБРАЖЕНИЯ ЗЕМНОЙ ПОВЕРХНОСТИ И ВОЗДУШНОЙ ОБСТАНОВКИ С ПОМОЩЬЮ АНТЕННОЙ РЕШЕТКИ

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

Формула изобретения RU 2 539 558 C1

Способ формирования трехмерного изображения земной поверхности и воздушной обстановки с помощью антенной решетки в виде M×N-матрицы приемных элементов, заключающийся в том, что при наблюдении за земной поверхностью или воздушной обстановкой с помощью сканирующего по углам луча антенны в зоне обзора измеряют при каждом положении луча в q, k-x измерительных каналах ( q = 1 , M ¯ , k = 1 , N ¯ ) значения комплексных амплитуд принимаемых сигналов y ˙ q , k , формируют из этих амплитуд M×N-матрицу измерений Y с элементами y ˙ q , k , q = 1 , M ¯ , k = 1 , N ¯ и далее обрабатывают,
отличающийся тем, что при сканировании последовательно смещают луч на ширину диаграммы направленности антенны (на уровне 0,5 мощности) построчно в зоне обзора, комплексные амплитуды y ˙ q , k ( l ) измеряют в каждом l-м элементе разрешения дальности ( l = 1 , L ¯ , L - число таких элементов), запоминают их в M×N-матрице измерений Y и обрабатывают в антенной системе координат, при этом умножают матрицу Y слева на рассчитанную заранее m×M-матрицу весовых коэффициентов WA и получают m×N-матрицу S восстановленного по строкам изображения: S=WAY, далее матрицу S умножают справа на рассчитанную заранее N×n-матрицу весовых коэффициентов WB и получают m×n-матрицу U восстановленного по строкам и столбцам изображения с повышенным в mn раз (m<N, n<N) угловым разрешением: U=SWB, затем вычисляют модули u ( i , j , l ) = | u ˙ ( i , j , l ) | i, j-х комплексных элементов матрицы U ( i = 1 , m ¯ , j = 1 , n ¯ ) , повторяют операции для всех l-х элементов дальности ( l = 1 , L ¯ ) и запоминают координаты тех i, j, l-x пространственных элементов, амплитуды которых превышают порог γ0 обнаружения полезного сигнала: u(i,j,l)>γ0, в результате получают трехмерное изображение объектов отражения в виде совокупности их прямоугольных координат xj, yj, zl в зоне видимости антенны, затем повторяют операции для всех положений луча и получают трехмерное изображение в полной зоне обзора.

Документы, цитированные в отчете о поиске Патент 2015 года RU2539558C1

СПОСОБ ПОВЫШЕНИЯ РАЗРЕШАЮЩЕЙ СПОСОБНОСТИ ФАЗИРОВАННОЙ АНТЕННОЙ РЕШЕТКИ БОРТОВОЙ СТАНЦИИ 2011
  • Клочко Владимир Константинович
RU2464680C1
СПОСОБ ПОЛУЧЕНИЯ ТРЕХМЕРНОГО РАДИОЛОКАЦИОННОГО ИЗОБРАЖЕНИЯ ВРАЩАЮЩЕГОСЯ ПО КУРСУ, ТАНГАЖУ И КРЕНУ ОБЪЕКТА ПРИ МНОГОЧАСТОТНОМ ИМПУЛЬСНОМ ЗОНДИРОВАНИИ 2006
  • Блиновский Александр Михайлович
  • Мезенцев Павел Викторович
  • Половинкин Леонид Петрович
  • Половинкин Алексей Леонидович
RU2327190C1
Передняя вилка мотоцикла 1959
  • Игнатов Г.В.
SU124820A1
US 6563451 B1, 13.05.2003
EP 1895472 A2, 05.03.2008
Способ получения этилового эфира 9-фенилгидразоно-6-метил-4-оксо-6,7,8,9-тетрагидро-4н-пиридо/1,2-а/-пиримидин-3-карбоновой кислоты 1980
  • Иштван Хермец
  • Золтан Месарош
  • Тибор Брайнинг
  • Шандор Вираг
  • Лелле Вашвари
  • Агнеш Хорват
  • Габор Надь
  • Аттила Манди
  • Тамаш Сюч
  • Иштван Биттер
  • Дьюла Шебештьен
SU978730A3

RU 2 539 558 C1

Авторы

Клочко Владимир Константинович

Усачев Алексей Николаевич

Даты

2015-01-20Публикация

2013-11-05Подача