СПОСОБ ВЫДЕЛЕНИЯ СПЕКТРАЛЬНЫХ ОТСЧЕТОВ В МНОГОКАНАЛЬНОЙ ДОПЛЕРОВСКОЙ РЛС
METHOD OF DISTINGUISHING SPECTRAL SAMPLES IN MULTICHANNEL DOPPLER RADAR

Изобретение относится к радиолокации, а именно к бортовым импульсно-доплеровским радиолокационным станциям (РЛС), работающим в режиме узкополосной доплеровской фильтрации [1] и предназначенным для наблюдения за наземными или воздушными объектами. Традиционно в таких РЛС используется один пространственно-измерительный канал первичной обработки принимаемых антенной отраженных сигналов. На выходе тракта первичной обработки в каждом элементе разрешения дальности и доплеровской частоты получается комплексное измерение, амплитуда которого свидетельствует о наличии или отсутствии объекта отражения в пространственном объеме элемента разрешения.
Однако по двум координатам - дальность и доплеровская частоты невозможно определить пространственные координаты объекта. Существуют способы [2, 3] измерения пространственных координат элементов земной поверхности или объектов на поверхности, основанные на использовании нескольких пространственно-измерительных (в дальнейшем - измерительных) каналов первичной обработки по числу пространственно разнесенных приемных элементов антенной системы (АС). Запаздывание или опережение по фазе сигналов, принимаемых в приемных элементах АС, дает информацию о пространственном положении наземных или воздушных объектов.
Рассмотрим в качестве прототипа способ измерения координат элементов земной поверхности (или объектов на поверхности) в бортовой четырехканальной доплеровской РЛС [2]. Способ заключается в формировании на заданном промежутке времени синтезирования радиолокационного изображения участка земной поверхности в виде совокупности комплексных амплитуд (измерений) сигналов отражения в элементах разрешения дальности и доплеровской частоты одновременно в четырех измерительных каналах, определении тех частот, на которых амплитуда сигнала превышает порог обнаружения, и последующей обработке совокупности полученных измерений, приводящей к определению пространственных координат отражающих элементов земной поверхности.
К аналогичным действиям приводит и способ [3] формирования трехмерного изображения земной поверхности в бортовой четырехканальной доплеровской РЛС, который также можно рассматривать в качестве прототипа. Отличие [2] и [3] заключается лишь в разных правилах обработки четырех комплексных амплитуд. Определение доплеровских частот в [2] и [3], на которых амплитуда отраженного сигнала превышает порог обнаружения одновременно в четырех каналах, сопровождается выделением спектральных отсчетов (комплексных измерений), соответствующих этим частотам.
Рассмотренные способы-прототипы обладают следующим недостатком. При определении доплеровских частот, на которых амплитуда сигнала превышает порог обнаружения одновременно во всех измерительных каналах, предполагается, что флуктуация доплеровской частоты в каналах отсутствует. В действительности из-за колебаний доплеровской частоты при приеме отраженных сигналов в пространственно разнесенных элементах АС происходит «растекание» частот в соседних элементах разрешения частоты. Это приводит к появлению ложных спектральных отсчетов в доплеровских спектрах и изменению их амплитуды. Возникает проблема выделения в нескольких измерительных каналах спектральных отсчетов, соответствующих одним и тем же отражающим элементам - объектам, в условиях колебания частоты.
Технический результат направлен на устранение указанного недостатка.
Технический результат предлагаемого технического решения достигается применением способа выделения спектральных отсчетов в многоканальной доплеровской радиолокационной станции (РЛС), заключающегося в выделении в доплеровских спектрах измерительных каналов РЛС спектральных отсчетов (комплексных измерений), соответствующих разным объектам отражения, и обработке полученных измерений, приводящей к определению пространственных координат объектов, отличающийся тем, что при выделении спектральных отсчетов в доплеровских спектрах измерительных каналов РЛС, в условиях колебания доплеровской частоты, формируют непересекающиеся группы спектральных отсчетов, взятых по одному из каждого спектра, и выбирают по критерию наименьшей суммы квадратов отклонения параметров групп относительно их среднего значения наиболее правдоподобные группы, соответствующие разным объектам отражения, выделенные отсчеты используют при определении пространственных координат объектов. Алгоритмически способ заключается в следующем.
q
q
q
q
q
1
q
1. В каждом q-м измерительном канале (, Q - число каналов) амплитуды доплеровских спектров ,  (N - количество частот в спектре) сравниваются с порогом р обнаружения полезного сигнала. Если , то формируется k-й элемент  множества спектральных отсчетов в q-м канале Ω, , (m≤N, m - число таких отсчетов). В результате образуются множества Ω, . Параметры элементов множеств Ω запоминаются как координаты точек в 3-мерном пространстве: x(k,q)=j(k) - доплеровская частота,  - амплитуда спектрального отсчета,  - его фаза (аргумент комплексного измерения).
1
2. На первом шаге q=1 рассматривается множество Ω. Каждый k-й элемент  этого множества дает начало ρ-й группы, для нее устанавливаются начальные значения показателя правдоподобия и оценок координат центра ρ-й группы:
I(ρ,1)=0, , .
Номера k-x элементов запоминаются в массиве K(ρ, 1)=р, .
q
3. Для последующих множеств Ω (q=2, 3, …, Q) выполняются следующие операции.
q
1
q
3.1. Каждой ρ-й группе , полученной на предыдущем (q-1)-м шаге, ставятся в соответствие k-е элементы множества Ω. Выбору подлежат лишь те элементы, частота которых x (k,q)=j(k) удовлетворяет условию:
,
то есть попадает в доверительный интервал ρ-й группы, сформированной на предыдущем шаге. Интервал строится относительно экстраполированного значения частоты ρ-й группы на q-й шаг . Величина отклонения Δj выбирается с учетом возможных колебаний частоты (например: Δj=1 или 2).
L
L
Для таких k-x элементов вычисляется показатель правдоподобия, имеющий смысл суммы квадратов отклонений L-x координат точек относительно их среднего значения, взятых с определенными весами μ (μ>0), и вычисляемый по рекуррентной формуле
.
L
L
Веса μ (μ>0) выбираются эмпирически.
k
q
k
k
q
q
k
3.2. Показатель I сравнивается с порогом α, который выбирается с учетом статистических свойств величины I. Если I≤α, то k-й элемент множества Ω присоединяется к ρ-й группе и она получает свое продолжение с новым порядковым номером g на q-м шаге (g=1, 2, …). Для g-й группы запоминается показатель качества: I(g,q)=I, уточняются оценки средних значений L-x координат по рекуррентной формуле
, ,
и запоминается номер присоединенного элемента: K(g,q)=k.
q
3.3. Если к ρ-й группе, полученной на предыдущем (q-1)-м шаге, не присоединен ни один из элементов множества Ω, то эта группа считается неперспективной и не получает дальнейшего продолжения.
q
q
3.4. После выполнения операций п.п. 3.1-3.3 формируется n наиболее перспективных групп с оценками параметров , , , показателями правдоподобия I(g,q) и номерами K(g,q) присоединенных к g-м группам элементов множества Ω.
p
1
2
1
2
4. Для последнего номера q=Q после выполнения операций раздела 3 выделяются n наилучших непересекающихся ρ-x групп с наименьшими показателями I(g,Q), . Выделение групп производится последовательно в порядке увеличения показателей I(g,Q),I(g,Q), … где g - номер первой выделенной группы, g - номер второй и т.д. После выделения ρ-й группы (ρ=1, 2, …) из оставшихся групп исключаются те, номера элементов которых совпадают с номерами элементов выделенной группы. Для этого используется массив запомненных номеров K(g,q), , . Допускается минимально количество π одинаковых элементов групп (например, π=1).
Операции п. 4 заканчиваются после выделения очередной ρ-й группы по признаку отсутствия оставшихся групп, при этом оценкой n является последний номер ρ: n=ρ.
ρ
ρ
1
2
Q
5. Спектральные отсчеты , , , k=K(g,q), выделенных элементов множеств Ω, Ω, …, Ω передаются на алгоритмы оценивания координат.
Замечание. Возможно упрощение. Для этого вводятся ограничения:
q
k
а) к ρ-й группе, сформированной на (q-1)-м шаге, можно отнести только один k-й элемент множества Ω, доставляющий показателю I наименьшее значение;
б) выбранный k-й элемент не может быть привязан к другой g-й группе.
q
Это исключает возможность размножения вариантов, но нарушает оптимальность выбора на последнем Q-м шаге в смысле критерия правдоподобия. Также отпадает необходимость проверки пересечений выбранных групп, так как каждый элемент множества Ω в этом случае не может войти в состав более чем одной группы.
Выводы. Предложенный способ позволяет выделять спектральные отсчеты в доплеровских спектрах, полученных в нескольких измерительных каналах, в условиях колебания доплеровской частоты. Выделенные спектральные отсчеты, отнесенные к разным объектам отражения, подвергаются операциям определения пространственных координат объектов. Способ может быть использован в существующих доплеровских системах наблюдения за наземными и воздушными объектами.
Источники информации
1. Кондратенков Г.С., Фролов А.Ю. Радиовидение. Радиолокационные системы дистанционного зондирования Земли: учеб. пособие для вузов / под ред. Г.С. Кондратенкова. М.: Радиотехника, 2005, 368 с.
2. Патент RU 2534224 C1. Способ измерения координат элементов земной поверхности в бортовой четырехканальной доплеровской РЛС.
3. Патент RU 2572357. Способ формирования трехмерного изображения земной поверхности в бортовой четырехканальной доплеровской РЛС.