close

Вход

Забыли?

вход по аккаунту

?

Адаптивная оптимизация характеристик активных поляризационных локаторов.

код для вставкиСкачать
3
РАДИОЭЛЕКТРОНИКА
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
Ошлаков В. Г., 2Илюшин Я. А.
1
Институт оптики атмосферы им. В.Е. Зуева, Сибирское отделение Российской академии наук, http://www.iao.ru
634021 Томск, Российская Федерация
1
Московский государственный университет им. М.В. Ломоносова, http://www.phys.msu.ru
119991 Москва, Российская Федерация
Поступила в редакцию 05.11.2012, после доработки 28.05.2013
2
Рассмотрены задачи адаптивной оптимизации характеристик активного локатора на основе
известных данных по поляризационным свойствам рассеяния цели и подстилающей поверхности.
Дано единое описание всех элементов схемы локатора в терминах матриц Мюллера и параметров
Стокса оптического излучения. Представлена структурная схема адаптивного однопозиционного
поляризационного локатора, а также схема поляриметра. Рассмотрена задача оптимальной селекции
целей с известными поляризационными характеристиками рассеяния на фоне помех, а также задача
идентификации одиночных и групповых целей на фоне подстилающей поверхности с известными
характеристиками рассеяния. Проведены расчеты матриц рассеяния пресной воды и гидрозоля на
длине волны 440 нм и график контрастности гидрозоля на фоне рэлеевского рассеяния пресной воды.
Проведены расчеты матриц рассеяния атмосферного аэрозоля и контраста на фоне рэлеевского
рассеяния воздуха в приземном слое атмосферы на длинах волн 300 и 550 нм. Показана возможность
повышения контраста цели с использованием информации о поляризационных свойствах рассеяния
цели и окружающего фона.
Ключевые слова: контраст цели, подстилающая поверхность, зондирующий сигнал, поляризационные
характеристики, вектор Стокса, матрицы Мюллера, гидрозоль, аэрозоль
УДК 535.36+530.1
содержание
1. введение (3)
2. адаптивный поляризационный локатор (4)
3. численный эксперимент (11)
3.1. зондирование гидрозоля (12)
3.2. зондирование аэрозолей (12)
4. заключение (13)
литература (13)
1. ВВЕДЕНИЕ
Актуальная задача повышения эффективности
активных локаторов решается, как правило, на
путях уменьшения вероятности ложной тревоги
и вероятности пропуска цели, что в значительной
мере определяется отношением сигнал/шум
на входе решающего устройства и выбором
пространства признаков цели. Так в работе [1]
показано, что пространственная, временная,
частотная, фазовая и поляризационная селекции
полезного сигнала на фоне помех повышают
отношение сигнал/шум. В работе [2] для защиты
приема от сильных сигналов, отраженных от
гидрометеоров, рассмотрена поляризационная
селекция. В настоящей работе предлагается
использовать для этих целей поляризационные
характеристики всех элементов схемы активного
локатора.
Цели могут располагаться в объеме и
на поверхности. Назовем все, что создает
пассивные помехи при локации, подстилающей
поверхностью, а активные и пассивные помехи
фоном. Диаграмма направленности (ДН)
антенны определяет телесный угол, в пределах
которого
распространяется
зондирующий
импульс, характеризуемый несущей частотой f0
и шириной спектра Δf. Если Δf/f << 1, то такой
сигнал является квазимонохроматическим.
В дальнейшем будем рассматривать локацию
квазимонохроматическими
импульсами.
При
локации целей, находящихся в объеме, разрешающий
объем с одной стороны ограничен шириной
ДН антенны, а с другой стороны по дальности
величиной cτp/2, где c – скорость распространения
зондирующего импульса, а τp – его длительность.
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
4 ОШЛАКОВ В.Г., ИЛЮШИН Я.А.
При локации целей, находящихся на
поверхности,
ограничение
разрешаемой
площадки поверхности по дальности в отличие
от зондирования целей в объеме, становится
равным примерно cτp/2cosγ, где γ – угол между
осью ДН антенны и зондируемой поверхностью.
В каждый данный момент времени сигнал,
действующий на вход приемника, является
результатом сложения сигналов, отраженных
от элементарных отражателей и от цели,
расположенных случайным образом в пределах
разрешаемого объема или разрешаемой
площадки. Рассматриваемая нами локация
занимает релеевскую, квазиоптическую и
резонансную области, а габаритные размеры
цели ограничиваются условием: при дискретном
сканировании ДН антенны цель должна
вписываться в размеры разрешающего объема
или разрешаемой площади.
Если определять электромагнитное поле
вектором напряженности электрического поля,

то электромагнитное поле E (t ) в любой точке
приемной антенны, в общем случае, является
векторной суммой полей, обусловленных
рассеянием зондирующего импульса целью


Et (t ) , подстилающей поверхностью Es (t ) и

поля внешних источников Eext (t ) . Эти поля
могут быть описаны в двухкомпонентных
обозначениях [3], в частности, в декартовой
системе координат XYZ в виде



E (t ) = Ex(t) i + Ey(t) j = (Etx(t) + Esx(t) +


Eextx(t)) i + (Ety(t) + Esy(t) + Eexty(t)) j ,
(1)
где Ex(t), Etx(t), Esx(t), Eextx(t) и Ey(t), Ety(t), Esy(t),



Eexty(t) – составляющие полей E (t ) , Et (t ) , Es (t )

 
и Eext (t ) по осям Х и Y соответственно; i , j –
орты осей соответственно Х и Y.
Представление (1) позволяет описывать

поляризацию электромагнитного поля E (t ) .
Как известно [4], в оптимальном приемнике
детерминированного сигнала в шумах все
компоненты
его
спектра
суммируются
арифметически в некоторый момент времени
и в результате в этот момент достигается
максимальное превышение сигнала над шумом.
Аналогично приемник поляризованного

электромагнитного поля E (t ) можно снабдить
поляризационным
устройством,
которое
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
РАДИОЭЛЕКТРОНИКА
осуществляет
алгебраическое
сложение
компонент Ex(t) и Ey(t), умноженных на
коэффициенты ax и ay соответственно, которые
определяют поляризационные характеристики
приемника электромагнитного поля.
Если E(t) – комплексный сигнал, то сигнал
на выходе поляризационного устройства E∑(t)
определяется в виде
E∑(t) = Re{axEx(t)} + Re{ayEy(t)}
(2)
где ax, ay – комплексные числа; Re{...} –
действительная часть {...}.
Всегда возможны такие ax и ay, что в (2)
компоненты Re{axEsx(t)}, Re{ayEsy(t)}, Re{axEextx(t)},
Re{ayEexty(t)} в результате арифметического
вычитания полностью компенсируют друг
друга, в то время как компоненты Re{axEtx(t)} и
Re{ayEty(t)} дадут величину отличную от нуля.
Назовем адаптивным поляризационным локатором
(АПЛ) активный локатор, который позволяет
ослабить фоновый сигнал по отношению к сигналу
цели за счет создания оптимальных поляризационных
характеристик зондирующего сигнала и приемника
[5, 6, 7] при известных матрицах рассеяния цели и
подстилающей поверхности цели.
Матрица
рассеяния
подстилающей
поверхности
цели
и
оптимальные
поляризационные характеристики приемника,
а также контраст цели образуют признаковое
пространство цели.
Далее рассмотрим один из вариантов АПЛ [8].
2. АДАПТИВНЫЙ ПОЛЯРИЗАЦИОННЫЙ
ЛОКАТОР
Каждая
радиолокационная
цель
может
быть представлена некоторым
оператором

(оператором рассеяния) A , отображающим
пространство параметров падающей на цель
волны в пространство параметров отраженной
волны. Представление цели оператором
рассеяния математически корректно ставит
задачу оптимизации локатора [3].
Зондирующий сигнал и фон характеризуются

своими векторами Стокса, соответственно S0
= (I0Q0U0V0)T и Sb = (IbQbUbVb)T, где T – знак
транспонирования, а цель и подстилающая
поверхность характеризуются своими матрицами
рассеяния, соответственно [9, 10]
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
РАДИОЭЛЕКТРОНИКА
Sc
t/
S0
Sф
4
τ, α
θ, δ
1
β, γ
5
A
2
7
S
8
θ, δ
τ, α
M
Y
Z
Х
0
6
β0, γ0
Z
∆
t′′
3
Рис. 1. Структурная схема АПЛ: 1 – источник
излучения, 2 – измеритель контраста, 3 – компьютер, 4 –
приемопередающая антенна, 5 – поляриметр, 6 – основание,
7 – индикатор контраста, 8 – поляризационный фильтр.
A = (aij), i,j = 1, 4 ,
M = (mij), i,j = 1, 4 .
Структурная схема одного из вариантов
АПЛ (моностатическая схема зондирования,
вектор Стокса зондирующего сигнала и вектор
Стокса излучения, падающего на приемник,
определяются в системе XYZ) показана на рис. 1.
Источник излучения 1 содержит: генератор
линейно-поляризованного
излучения,
вращатель плоскости поляризации и линейное
фазосдвигающее устройство, и генерирует
электромагнитную волну, вектор напряженности

электрического поля E которой представим
ортогональными компонентами, направленными


вдоль осей X и Y, соответственно Ex и E y в
виде


Ex = Excosωt i ,


E y = Eycos(ωt + δ) j ,
где Ex и Ey – амплитуды компонент соответственно


Ex и E y ; ω – круговая частота генерируемой
электромагнитной волны; t – время; δ – фазовый


сдвиг компонент Ex и E y .
Излучение
источника
полностью
поляризовано и его поляризация описывается
параметрами:
Если излучение нормировано, то I0 = 1.
Зондирующий импульс излучается генератором
1 при поступлении импульса t'.
Сначала примем, что антенна 4 и
пространство, в котором распространяются
излучаемая и отраженная от цели и
подстилающей поверхности электромагнитные
волны, не изменяют их поляризацию. Ниже
будет рассмотрен случай, когда это положение
не выполняется.
Отраженный от цели, а также фоновый
сигналы, которые описываются векторами

Стокса St = (ItQtUtVt)T и Sb = (IbQbUbVb)T, через
антенну 4 поступают на поляриметр 5 (рис. 2),
который
измеряет параметры вектора Стокса

T
=
(IQUV)
поступившего излучения. При
S
наличии одновременно излучения, отраженного

от цели и излучения фона, вектор Стокса S
определяется
в виде



T
S = (IQUV) = ( St + Sb ) =
= (It + Ib, Qt + Qb, Ut + Ub, Vt + Vb)T.
Поляриметр 5 содержит: вращающийся
фазовый элемент (оптимальный фазовый
сдвиг, который должен создавать вращающийся
фазовый
элемент,
составляет
131,795°),
поляризатор,
измеритель
интенсивности,
процессор и контроллер [11, 12].
Антенна 4 закреплена на основании 6,
которое может вращаться вокруг двух взаимноортогональных осей с помощью управляемого
I0 = Ex2 + E y2 ; Q0 = I0cos2θ;
U0 = I0sin2θcosδ; V0 = I0sin2θsinδ, (3) Рис. 2. Структурная схема поляриметра: 1 –
где θ – угол поворота вращателя плоскости вращающийся фазовый элемент, 2 – поляризатор,
3 – измеритель интенсивности, 4 – процессор, 5 –
поляризации.
электрический привод, 6 – контроллер.
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
5
6
ОШЛАКОВ В.Г., ИЛЮШИН Я.А.
электропривода. Кроме того, при некотором
положении основания 6 компьютер 3, изменяя
положение ДН антенны 4, сканирует ею в некотором телесном угле Ω с наперед заданным
шагом изменения углов β и γ, соответственно
по азимуту и углу места, в диапазоне Δβ
и Δγ (рис. 3). Шаг изменения углов β и γ
выбирается меньше ширины ДН антенны 4
по уровню 0,7. Сканируемый телесный угол
Ω разбивается на Z точек. Будем считать, что
перемещение ДН производится построчно
слева направо и сверху вниз. Если положение
ДН соответствует точке с номером Zi, азимут
и угол места которой соответствуют βi и γi, то
будем обозначать это положение в виде ДН
{βi, γi}.
 Измерение параметров вектора Стокса
S сигнала, поступившего с антенны 4,
производится для каждой точки Zj телесного
угла Ω в поляриметре 5 следующим образом:
при поступлении нового числа Zj контроллер с
помощью электрического привода устанавливает
быструю ось вращающегося фазового элемента
относительно оси X под следующими углами:
ξ1 = 38,54°; ξ2 = 75,14°; ξ3 = 105,38°; ξ4 = 141,857°
[11, 12].
Процессор 4 запоминает данные измерителя
интенсивности 3 в момент поступления импульса
t″. При каждом установленном значении ξ
достаточно поступления одного импульса t″.
Так как зондирующий импульс излучается
генератором 1 при поступлении импульса t', то
временной интервал между импульсами t' и t″
определяет расстояние АПЛ до зондируемого
пространства. Далее процессор по этим четырем
РАДИОЭЛЕКТРОНИКА
измеренным величинам интенсивности, путем
решения системы четырех линейных
уравнений,

определяет параметры Стокса S = (IQUV)T для
точки пространства, соответствующей временному
интервалу между импульсами t' и t″ [11, 12].
Изменением временного интервала между
импульсами t' и t″ определяют значения
векторов Стокса сигнала, отраженного от точек
пространства по всей дальности.
Поляризационный фильтр 8 может быть
реализован на базе микропроцессора и
преобразует
параметры I, Q , U, V вектора Стокса

,
поступающего
с поляриметра 5, так как это
S
выполняют последовательно-фазовый элемент
(с параметрами: α – угол между направлением
его быстрой оси и осью X; τ – сдвиг фаз
ортогональных компонент) и поляризатор,
плоскость пропускания которого совпадает с
осью X.
Электрические
сигналы
на
выходе
поляризационного фильтра 8 пропорциональны
параметрам вектора Стокса:


S s = Pr Pr′ S = (IsQsYsVs)T, (4)
где Pr и Pr′ – матрицы Мюллера поляризатора и
фазового элемента соответственно.
При облучении цели
зондирующим сигналом

с вектором Стокса S0 = (I0Q0U0V0)T на приемник
попадает
отраженный сигнал с вектором Стокса

St , который можно записать в виде


St = (ItQtUtVt)T = CA S0 (5)
где C – коэффициент пропорциональности,
определяемый расстоянием до зондируемого
объекта и параметрами приемопередающей
системы локатора.
На выходе поляризационного фильтра ему

соответствует вектор Стокса Sts = ( I ts QtsU tsVt s )T .
Параметр I ts можно записать
в виде скалярного


произведения векторов St и P в виде:
 
(6)
I ts = ( St , P)
где
T
  1 + cosτ 1 − cosτ
1 − cosτ

+
P=
cos 2α ,
sin 4α , − sin τ sin 2α  .
1,
2
2
2



Вектор Стокса Sb излучения, отраженного
Рис. 3. Следы ДН антенны на сферической поверхности подстилающей поверхностью, определяется в
с центром в точке расположения АПЛ.
виде
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
РАДИОЭЛЕКТРОНИКА


Sb = (IbQbUbVb)T = CM S0 .
На выходе поляризационного фильтра 8 ему

соответствует вектор Стокса Sb = ( I bs QbsU bsVbs )T .
Параметр I bs можно записать
в виде скалярного


произведения векторов Sb и P в виде:
 
(7)
I bs = ( Sb , P).
В общем случае электрический сигнал
на выходе поляризационного фильтра 8
s
пропорционален I=
I ts + I bs .
s
It
Отношение s = N характеризует контраст


Ib
цели на фоне помех.
Оптимальные S0 и P
обеспечивают максимум N.
Параметры Q0, U0, V0, обеспечивающие
максимум
N при условии полной поляризации

S0 , определяются параметрами приемника τ и α
фазового элемента поляризационного фильтра
и матрицей рассеяния цели A и не зависят от

Sb [5, 6]. При любых τ, α и A вектор Стокса S0
при условии полной поляризации и нормировки
имеет вид:



A3
A2
A4
S0 = 1,
,
,
2
2
2
2
2
2
2

A +A +A
A +A +A
A + A2 + A2
2
3
4
2
3
4
2
3
4

где
A=
2 a12 + a22
+ a32
1 + cos τ
1 – cos τ
cos 4α +
+ a22
2
2
1 + cos τ
1 – cos τ
cos 4α +
+ a23
2
2
1 – cos τ
sin 4α – a43 sin τ sin 2α,
2
A=
4 a14 + a24
+ a34
(8)
1 – cos τ
sin 4α – a42 sin τ sin 2α,
2
A=
3 a13 + a23
+ a33
T


 ,


1 + cos τ
1 – cos τ
cos 4α +
+ a24
2
2
1 – cos τ
sin 4α – a44 sin τ sin 2α.
2
Таким образом, если известна только матрица
рассеяния цели A, то максимум максиморум N
можно получить практически путем перебора
с наперед заданным шагом параметров τ и
α в диапазоне 0 ≤ τ ≤ 2π рад, 0 ≤ α ≤ π рад с
одновременным
изменением
поляризации
источника излучения 1 согласно (8), которое
удовлетворяется при
1
θ A = arccos
2
−1
±1
1
, δ A =arcsin
.
2
A32 + A42
A32
1+
1+
A42
A22
(9)
Если
фон
создается
отражением
зондирующего сигнала от подстилающей
поверхности цели с известной
матрицей

рассеяния M, то вектор Стокса S0 , при котором
N достигает максимума (при условии полной
поляризации и нормировки зондирующего
импульса), имеет вид



M3
M2
M4
1, −
S0 =
,−
,−

2
2
2
2
2
2
2
M +M +M
M +M +M
M + M2 + M2

2
3
4
2
3
4
2
3
4

где
M=
2 m12 + m22
+ m32
1 + cos τ
1 – cos τ
cos 4α +
+ m22
2
2
1 – cos τ
sin 4α – m42 sin τ sin 2α,
2
M=
3 m13 + m23
+ m33
T

 ,(10)



1 + cos τ
1 – cos τ
cos 4α +
+ m23
2
2
1 – cos τ
sin 4α – m43 sin τ sin 2α,
2
M=
4 m14 + m24
1 + cos τ
1 – cos τ
cos 4α +
+ m24
2
2
1 – cos τ
sin 4α – m44 sin τ sin 2α.
2
Если
известна
матрица
рассеяния
подстилающей поверхности цели M, но
неизвестна матрица рассеяния цели A, то
операция нахождения максимума максиморума
величины контраста N любой цели может
быть выполнена практически путем перебора
с наперед заданным шагом параметров τ и
α в диапазоне 0 ≤ τ ≤ 2π рад, 0 ≤ α ≤ π рад с
одновременным изменением вектора Стокса S0
источника излучения 1 согласно (10), которое
удовлетворяется при
+ m34
1
θ M =arccos
2
±1
1
−1
, δ M =arcsin
.
2
2
2
M 32
M3 + M 4
1+
1+
M 42
M 22
(11)
Параметры τ и α, при которых достигается
максимальная
величина
контраста
цели
Nmax, определяются матрицей рассеяния
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
7
8
РАДИОЭЛЕКТРОНИКА
ОШЛАКОВ В.Г., ИЛЮШИН Я.А.
подстилающей поверхности цели M и матрицей
рассеяния цели A, которая является признаковой
характеристикой цели. Поэтому признаковое
пространство цели образуют параметры τ̂ и
α̂ , при которых достигается максимальное
значение контраста цели N̂ max , матрица
рассеяния подстилающей поверхности цели
M и максимальное значение контраста цели
N̂ max . Матрица A определяется путем решения
обратной задачи, используя τ̂ , α̂ и N̂ max .
Рассмотрим работу измерителя контраста
2. Информационное поле каждой точки i
сканируемого пространства Ω содержит данные
о параметрах этой точки: i – номер точки
(соответствует угловому положению ДН), время
импульса t″, азимут этой точки βi и ее угол места
γi, которые поступают с компьютера 3, а также
значение IS, поступившее с поляризационного
фильтра 8.
Алгоритм работы измерителя контраста 2
показан на рис. 4.
На 1-, 2-, 3-, 4-м шагах формируется массив
данных о Z точках сканируемого пространства
Ω, который содержит Z информационных
полей. После того, как сформирован этот
массив данных, выполняется 6-й шаг, на котором
вычисляется контраст N каждой точки, начиная
с 1-й. Например, для j-й точки с координатами
{βj,γj}, удаленной от АПЛ на расстоянии,
соответствующем t″, вычисляется значение
{
}
s
параметра I t′′ β j , γ j , соответствующее фону
в этой точке, по формуле:
χ
{
I ts′′ β j , γ j
I ts′′
}
∑ Its′′ {βmj , γ mj }
m =1
=
{
χ
{
I ts′′
}
{
{
}
{
}
}
,
(13)
}
β j , γ j – значение параметра I вектора
где

Стокса S s при положении ДН {βj,γj} антенны,
соответствующее времени t″.
7-й шаг. Если при некоторых значениях
параметров βj = βˆ j , γj = γˆ j и t″ = tˆ′′
удовлетворяется неравенство
Nˆ tˆ′′ βˆ j , γˆ j ≥ N th ,
(14)
где Nth – пороговое отношение сигнал/фон,
которое определяется вероятностью ложной
тревоги и вероятностью пропуска цели, то
осуществляется 8-й шаг – на индикатор контраста
7 и компьютер 3 выдается электрический сигнал,
пропорциональный Nˆ tˆ′′ βˆ j , γˆ j . На экране
индикатора контраста 7 в точке, соответствующей
положению ДН βˆ j , γˆ j антенны, высвечивается
{
S
}
{
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
(12)
βmj , γ mj – значение параметра IS при
где
m-м соседнем положении ДН {βmj,γmj} антенны,
относительно положения ДН {βj,γj} антенны,
соответствующее времени t″. Положение
ДН {βmj,γmj} является соседним относительно
её положения {βj,γj}, если оно находится в
некоторой его окрестности; χ – число соседних
положений ДН антенны (для хорошего
усреднения необходимо брать χ >> 1).
Далее вычисляется значение контраста
N{βj,γj} точки с координатами {βj,γj} по формуле:
I ts′′ β j , γ j
N t′′ β j , γ j =
I ts′′ β j , γ j
Рис. 4. Структурная схема программы измерения
контраста точек, равноудаленных от АПЛ.
,
{
}
}
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
РАДИОЭЛЕКТРОНИКА
точка
с
яркостью,
пропорциональной
ˆ
ˆ
N tˆ′′ β j , γˆ j . Чем выше яркость точки на
экране, тем больше вероятность наличия
цели в направлении ДН βˆ j , γˆ j антенны на
расстоянии, определяемом временным интервалом
t″ = tˆ′′ . Значение контраста вычисляется для
каждой точке сканируемого пространства.
Если исследуемый объект (цель) дисперсно
распределен в среде, которая является
подстилающей поверхностью с известной
матрицей рассеяния M, то Nt″{βj, γj} определяется
по формуле:
{
}
{
{
N t′′ β j , γ j
}
{
}
}
I ts′′ β j , γ j
,
= s
I bt′′
(15)
где I bs – определяется по формуле (7).
Формулу (15) можно использовать при
зондировании одиночных целей, расположенных
на известной подстилающей поверхности.
Рассмотрим работу АПЛ в двух случаях:
1. Известна матрица рассеяния цели A и
неизвестна матрица рассеяния подстилающей
поверхности цели M.
2. Известна матрица рассеяния M подстилающей
поверхности цели, но неизвестна матрица
рассеяния цели A.
Параметры τ и α в поляризационном фильтре
8 принимают дискретные значения в диапазоне
0 ≤ τ ≤ 2π рад, 0 ≤ α ≤ π рад и устанавливаются с
компьютера 3.
Рассмотрим работу АПЛ в первом случае
(решается задача оптимальной селекции на фоне
помех целей, обладающих матрицей рассеяния
A).
Сначала устанавливаются минимальные
значения параметров τ и α из указанного выше
диапазона. В источнике излучения 1 вращатель
плоскости
поляризации
устанавливает
плоскость поляризации генератора линейно
поляризованного излучения на угол θA
1
θ A = arccos
2
−1
1+
A32 + A42
A22
,
а линейное фазосдвигающее устройство
осуществляет
сдвиг
фаз
ортогональных
компонент на угол δA,
1
δ M =arcsin
2
±1
M2
1+ 3
M 42
.
соответствующие этим τ и α согласно (9).
Поляризация
импульсов
источника
излучения 1 канала
 описывается согласно (3)
вектором Стокса S0

S0 = (I0, I0cos2θA, I0sin2θAcosδA, I0sin2θAsinδA)T.
Компьютер 3, изменяя положение ДН
антенны 4, дискретно сканирует пространство
в некотором телесном угле Ω при неизменных
параметрах τ, α, θA, δA. Излучение, отраженное
от цели и фоновое или только фоновое,
принимается в каждом j-м положении ДН
{βj,γj} антенной 4 и направляется вполяриметр
5, где измеряется вектор Стокса S = (IQUV)T
этого излучения, соответствующее моменту
поступления импульса t″. Электрические
сигналы, пропорциональные параметрам I ts ,
подаются на компьютер 3 и поляризационный
фильтр 8, где преобразуются виртуальным
фазовым элементом и поляризатором согласно

T
(4) в вектор Стокса Sts′′ = I ts′′ , Qts′′ , 0, 0
при
этих, установленных в них, параметрах τ и α.
Электрический сигнал, пропорциональный
параметру I ts , поступает в измеритель
контраста 2. После того, как ДН антенны
4 в процессе сканирования пройдет все Z
точек сканируемого угла Ω, в измерителе
контраста 2 будет определен контраст Nt″{βj,γj}
каждой этой точки по формуле (13). Далее
устанавливаются в поляризационном фильтре
8 большие значения τ и α, а в источнике
излучения соответствующие им значения
параметров θA и δA, и описанный выше
процесс повторяется. Процесс изменения τ и α
продолжается до достижения ими значений 2π
рад и π рад, соответственно. Если в некоторой
точке βˆ , γˆ на расстоянии, соответствующем
(
)
( )
временному интервалу tˆ′′ – t ′ , находится цель
с матрицей рассеяния A, то найдутся τ̂ и α̂
такие, при которых контраст цели Nˆ tˆ′′ βˆ j , γˆ j
{
}
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
9
10
РАДИОЭЛЕКТРОНИКА
ОШЛАКОВ В.Г., ИЛЮШИН Я.А.
будет на данном фоне максимальный и больше
Nth. В результате на индикаторе контраста 7
в точке, соответствующей положению ДН
βˆ j , γˆ j , будет яркая точка. Информация о
яркой точке ДН βˆ j , γˆ j антенны и величине
контраста Nˆ tˆ′′ βˆ j , γˆ j , времени tˆ′′ подается в
компьютер 3.
Рассмотрим работу АПЛ во втором случае.
В этом случае АПЛ решает задачу
идентификации объекта в точке, отличающейся
от других точек сканируемого пространства
характеристикой отражать электромагнитные
волны.
Рассмотрим работу АПЛ в случае, когда целью
является группа дисперсных рассеивателей или
одиночная цель.
Сначала устанавливаются минимальные
значения параметров τ и α из указанного
выше диапазона их изменений. В источнике
излучения 1 вращатель плоскости поляризации
устанавливает плоскость поляризации генератора
линейно-поляризованного излучения на угол θM
{
}
{
1
θ M =arccos
2
{
}
}
−1
,
M 32 + M 42
1+
M 22
а линейное фазосдвигающее устройство
осуществляет
сдвиг
фаз
ортогональных
компонент на угол δM, согласно (11)
1
δ M =arcsin
2
±1
M2
1+ 3
M 42
.
Поляризация
импульсов
источника
излучения 1 описывается согласно (3) вектором
Стокса

S0 = (I0, I0cos2θM, I0sin2θMcosδM, I0sin2θMsinδM)T.
Компьютер 3 изменяет положение ДН
антенны 4 и дискретно сканирует пространство
в телесном угле Ω при неизменных τ, α, θM, δM.
В каждом j-м положении ДН {βj,γj} антенны 4
излучение, отраженное от цели и фоновое или
только фоновое, принимается ей, направляется

в поляриметр 5 и измеряется вектор Стокса S
= (IQUV)T этого излучения, соответствующее
моменту поступления импульса t″.
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
Электрические сигналы, пропорциональные
параметрам IQUV, подаются на компьютер 3 и
поляризационный фильтр 8, где преобразуются
виртуальным
фазовым
элементом
и
поляризатором согласно (4) в вектор Стокса

T
Sts′′ = I ts′′ , Qts′′ , 0, 0
при этих, установленных τ
и α.
Электрический сигнал, пропорциональный
параметру I ts , поступает в измеритель контраста
2. После того, как ДН антенны 4 в процессе
сканирования пройдет все Z точек сканируемого
угла Ω, в измерителе контраста 2 будет
определен контраст Nt″{βj,γj} каждой j-й точки
положения ДН {βj,γj} антенны по формуле (13)
в случае одиночной цели и по формуле (15), в
случае, когда цель является группой дисперсных
рассеивателей.
Далее устанавливаются в поляризационном
фильтре 8 большие значения τ и α, а в
источнике излучения соответствующие им
значения параметров θM и δM, и описанный выше
процесс повторяется. Процесс изменения τ и α
продолжается до достижения ими значений 2π
рад и π рад соответственно. Если в некоторой
точке βˆ , γˆ на расстоянии, соответствующем
временному интервалу tˆ′′ – t ′ , находится цель
с характеристикой отражения, отличающейся
от характеристики отражения подстилающей
поверхности, то найдутся τ̂ и α̂ согласно (10),
при которых контраст цели Nˆ tˆ′′ βˆ j , γˆ j будет
на данном фоне максимальный и больше Nth. В
результате на индикаторе контраста 7 в точке,
соответствующей положению ДН βˆ , γˆ , будет
яркая точка. Информация о яркой точке (азимут
ˆ
β̂ , угол места γˆ и величина контраста Nˆ tˆ′′ β , γˆ
) подается в компьютер 3. Параметры τ̂ и α̂ , при
которых достигается максимальная величина
контраста цели N, определяется матрицей
рассеяния подстилающей поверхности М и
матрицей рассеяния цели А, которая является
признаковой характеристикой цели. Поэтому
признаковое пространство цели образуют
параметры τ̂ и α̂ , при которых достигается
максимальное значение контраста цели N̂ ,
матрица рассеяния подстилающей поверхности
М и значение контраста N̂ max .
(
)
( )
{
}
( )
{ }
РАДИОЭЛЕКТРОНИКА
Если антенна изменяет поляризацию
излучения источника излучения 1, а также
поляризацию
принимаемого
излучения,
то такие изменения можно представить
соответствующими
матрицами
Мюллера
антенны
[13] и записать


S0ex = T1 S0 ,


где S0 и S0ex – вектора Стокса излучения
источника излучения 1 и выходного излучения
антенны 4 соответственно; Т1– матрица Мюллера


Sex = T2 S , где S и
антенны
4
на
передачу,
и

Sex – вектора Стокса излучения, принимаемого
антенной 4, из пространства и излучения на ее
выходе соответственно; Т2 – матрица Мюллера
антенны 4 на прием.
Аналогично, если пространство изменяет
поляризацию распространяющегося в нем
излучения,

 то можно записать
S2 =T S1 , 
где S1 и S2 – вектора Стокса излучения в
пространстве в начальной и конечной точках
распространения соответственно; Т – матрица
Мюллера пространства.
Отнесем поляризационные изменения
антенны 4 и пространства к матрице
рассеяния цели А и матрице рассеяния
подстилающей поверхности цели М путем
введения приведенной матрицы рассеяния
цели Ar и приведенной матрицы рассеяния
подстилающей поверхности цели Mr, которые
определим в виде
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
дисперсных рассеивателей с матрицей
рассеяния М, которая принимается за
подстилающую поверхность, создающую
фоновый отраженный сигнал.
Обозначим вектор Стокса
излучения

зондирующего импульса S0 = (I0Q0U0V0)T,
параметры I0, Q0, U0, V0 определяются согласно
(10). Вектор Стокса S = (IQUV)T излучения,
отраженного от зондируемого объема, равен


S = (IQUV)T = C(A + M) S0 = (It + Ib, Qt +
Qb, Ut + Ub, Vt + Vb)T,
где параметры It, Qt, Ut, Vt обусловлены группой
зондируемых дисперсных рассеивателей, а
параметры Ib, Qb, Ub, Vb обусловлены группой
дисперсных рассеивателей, создающих
фон.

Параметры вектора Стокса Sb = (IbQbUbVb)T на
выходе поляризационного фильтра 8 определяются


по формуле S s = Pr Pr′ S . Параметр IS можно
записать в виде:
IS = I ts + I bs ,
где I ts – интенсивность полезного сигнала;
I bs – интенсивность фонового
сигнала.

Если вектор Стокса S0 определяется по
формуле (10), то [3, 4]:
I bs
= M1 − M 22 + M 32 + M 42 .
C
Рассмотрим отношение
Is
I ts
=+
1 N,
(16)
I bs
I bs
где N – контраст зондируемой группы дисперсных
Ar = T1TATT1
рассеивателей.
и
Если изменять параметры поляризационного
Mr = T1TMTT1.
Теперь считая, что антенна 4 и пространство фильтра τ и α в диапазоне 0 ≤ τ ≤ 2π рад, 0
не изменяют поляризацию излучения, можно ≤ α ≤ π рад с одновременным изменением
импульса
использовать формулы (8)–(11), подставив в них вектора Стокса S0 зондирующего
S
и нормировать сигнал I с поляризационного
A вместо А и M вместо М.
r
=+
1
r
фильтра на величину M1 − M 22 + M 32 + M 42 , то
3. ЧИСЛЕННЫЙ ЭКСПЕРИМЕНТ
 Is 
Описание рассеяния электромагнитных волн
C 1 + ts 
она, согласно (16), будет равна  I  и, при
b 

с помощью матриц рассеяния экспериментально
отсутствии зондируемой группы дисперсных
подтверждено в оптическом диапазоне, поэтому
рассеивателей не будет изменяться.
покажем эффективность работы АПЛ в
Приводимые ниже расчеты на основе
оптическом диапазоне.
экспериментальных данных показывают, что при
Рассмотрим
зондирование
группы
изменении τ и α в указанных выше диапазонах с
дисперсных рассеивателей с матрицей
одновременным изменением вектора Стокса S0
рассеяния А, смешанной с другой группой
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
11
12
РАДИОЭЛЕКТРОНИКА
ОШЛАКОВ В.Г., ИЛЮШИН Я.А.
зондирующего импульса, найдутся такие τ и α,
при которых контраст N зондируемой группы
дисперсных рассеивателей будет максимальным
при ее присутствии в зондируемом объеме.
Контраст N зондируемой группы дисперсных
рассеивателей на выходе поляризационного
фильтра рассчитывался по формуле
I ts
1
N
= =
×
I bs M1 − M 22 + M 32 + M 42

 1 + cos τ 1 − cos τ
 
cos 4α  + 
+
 I t + Qt 
2
2

 ,
×
1 − cos τ


sin 4α − Vt sin τ sin 2α 
 +U t
2

где

 It   a11 a12 a13 a14  
  
 − M2
 Qt  =  a21 a22 a23 a24  
 U t   a31 a32 a33 a34   − M
2
  
 

V
a
a
a
a
 t   41 42 43 44 
 − M2

1


M 22 + M 32 + M 42 

,
2
2
2
M 2 + M3 + M 4 

M 22 + M 32 + M 42 

Рис. 5. Графики контрастности: а – контраст
гидрозоля, λ = 440 нм; б – контраст плотной дымки, λ
= 550 нм; в – контраст плотной дымки, λ = 300 нм;
г – контраст дымки, λ = 550 нм.
что в четырех точках контраст гидрозоля имеет
здесь a11 – коэффициент направленного максимум с N = 11.8. По значениям τ̂ и α̂ , при
светорассеяния, который отвечает условию которых N максимально, можно решать обратную
задачу – подбор элементов матрицы рассеяния
σ = ∫ a11 ( I, I 0 ) dω ′, ( I, I ) – угол между
гидрозоля A, т.е. определение типа гидрозоля и
0
4π
направлением облучения I 0 и направлением его массовой концентрации.
рассеяния I , dω′ – элемент телесного угла.
3.2. ЗОНДИРОВАНИЕ АЭРОЗОЛЕЙ
Результаты численного эксперимента по
3.1. ЗОНДИРОВАНИЕ ГИДРОЗОЛЯ
Наличие солей в воде приводит к таким зондированию атмосферы для двух вариантов
малым изменениям показателя преломления, оптической погоды (дымка и плотная дымка)
что молекулярное рассеяние на флуктуациях приведены в таблице 2. Значения элементов
плотности для пресной и соленой воды матриц релеевского рассеяния воздуха M у
поверхности Земли и аэрозоля А получены на
практически остается постоянным [14].
Полученные в численном расчете значения длинах волн λ = 300 и λ = 550 нм. Размерность
–1
–1
элементов матриц релеевского рассеяния воды элементов (км ср ) [16].
Таблица 2
M и гидрозоля A (размерность м–1ср–1) на длине
Матрицы рассеяния воздуха M и аэрозоля A
волны λ = 440 нм представлены в таблице 1.
на длинах волн λ = 550 и λ = 300 нм
Таблица 1
Матрицы рассеяния воды M и гидрозоля A
на длине волны λ = 440 нм
M
m∙10
-3
A
a∙10-3
m11
m22
m33
m44
0.387
0.352
-0.352
-0.352
a11
a22
a33
a44
1.213
1.024
-1.024
-0.8
На рис. 5а приведен расчетный график
контрастности гидрозоля на фоне релеевского
рассеяния пресной воды, на котором видно,
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
M
m22
m33
λ = 550 0.137
0.133
-0.133 -0.133 m∙10-2
λ = 300 0.17
0.165
-0.165 -0.165 m∙10-1
a22
a33
a44
д ы м - λ = 550 0.44
ка
λ = 300 0.8
0.44
-0.3
-0.3
a∙10-2
0.56
-0.56
-0.56
a∙10-2
плотная
λ = 550 0.14
0.14
-0.066 -0.066 -
дымка
λ = 300 0.61
0.61
-0.29
воздух
A
m11
a11
m44
-0.29
a∙10-1
РАДИОЭЛЕКТРОНИКА
На рис. 5б-г приведены разрезы по плоскости
симметрии расчетных графиков контрастности
аэрозолей на фоне релеевского рассеяния чистой
атмосферы. Вершины графиков имеют форму
двух овалов, значение контрастности во всех
точках овалов одинаковое. При значениях τ̂ и
α̂ , соответствующих вершины графиков (а их
большое количество) достигается максимальный
контраст аэрозоля на фоне релеевского рассеяния
воздуха.
Анализ
графиков
показывает,
что
наибольший контраст можно получить на длине
волны λ = 550 нм, поэтому данная длина волны
является предпочтительной при зондировании в
данном случае.
При зондировании АПЛ позволяет получать
графики контрастности и, используя их, решать
обратную задачу: подбор элементов матрицы
рассеяния аэрозоля А, которая дает аналогичный
график контрастности, а также оценивать
массовую концентрацию аэрозоля.
4. ЗАКЛЮЧЕНИЕ
В АПЛ используются поляризационные
свойства цели и ее подстилающей поверхности
для повышения контрастности цели. Повышение
контрастности цели улучшает точностные
характеристики локатора и расширяет его
возможности: позволяет решать обратную
задачу – определять матрицу рассеяния цели и
производить ее идентификацию.
БЛАГОДАРНОСТЬ
Авторы благодарят Д.Д. Малых за составление
программы при проведении численного
эксперимента и выполнение расчетов.
ЛИТЕРАТУРА
1. Свистов ВИ. Радиолокационные сигналы и их
обработка. М., Сов.Радио, 1977, с. 447.
2. Давыдов ПС, Сосновский АА, Хаймович ИА.
Авиационная радиолокация. Справочник. Под.
ред. Давыдова ПС. М., Транспорт, 1984, 116 с.
3. Варганов МЕ, Зиновьев ЮС, Астанин ЛЮ.
Под ред. Тучкова ЛТ. Радиолокационные
характеристики летательных аппаратов. М.,
Радио и связь, 1985, 236 с.
4. Вакман ДЕ. Сложные сигналы и принцип
неопределенности в радиолокации. М., Советское
радио, 1965, 304 с.
АДАПТИВНАЯ ОПТИМИЗАЦИЯ ХАРАКТЕРИСТИК
АКТИВНЫХ ПОЛЯРИЗАЦИОННЫХ ЛОКАТОРОВ
5. Oshlakov VG, Eremina T. The increase of the
efficiency of research tools of active radio and
optical sensing. Proc. SPIE Europe Remote Sens.,
6744B-90, Florence, Italy, 2007.
6. Ошлаков ВГ, Ошлаков ВК, Еремина
ТA.
Оптимальные
поляризационные
характеристики приборов, регистрирующих
рассеянное излучение. Оптика атмосферы и
океана, 2007, 20(7):635–640.
7. Oshlakov VG, Andreev MI, Malykh DD. An
adaptive lidar. Proc. SPIE, Remote Sens., 74730N,
Berlin, Germany, 2009.
8. Ошлаков ВГ. Способ комплексной локации
цели. Патент РФ 2416108С1: G01S13/95 от
25.08.2009, опубл. 10.04.2011. Бюл. № 10, 2011.
9. Bohren CF, Huffman DR. Absorption and Scattering
of Light by Small Particles. New York, Willey, 1983.
10.Van de Hulst HC. Light scattering by small particles.
New York, Willey, 1957.
11.Ошлаков ВГ. Оптимальный измеритель
матрицы рассеяния. Оптика атмосферы и
океана, 1992, 5(11):1190–1193.
12.Ошлаков ВГ, Барков ЮГ. Численный анализ
аппаратной матрицы поляризационного
измерителя. Оптика атмосферы и океана, 2002,
15(7):635–640.
13.Есепкина НА, Корольков ДВ, Парийский
ЮН. Радиотелескопы и радиометры. М., Наука,
1973, 415 с.
14.Иванов АП. Физические основы гидрооптики.
Минск, Наука и техника, 1978, 504 с.
15.Kenneth JV, Edward SF. Measurement of the
Muller Matrix for Ocean Water. Appl. Opt., 1984,
23(23):4427-4439.
16.Горчаков ГИ. Матрица рассеяния света
приземным воздухом. Изв. АН СССР. Физика
атмосферы и океана, 1966, 2(6):595–605.
Ошлаков Виктор Григорьевич
к.ф.-м.н., с.н.с.
Институт оптики атмосферы им. В.Е. Зуева СО РАН
1, ул. акад. Зуева, г. Томск 634021, Россия
oshlakov@iao.ru
Илюшин Ярослав Александрович
к.ф.-м.н., доцент, чл.-корр. РАЕН
Московский государственный университет им. М.В.
Ломоносова
1/2, Ленинские горы, Москва 119991, Россия
+7 495 939 3252, ilyushin@physics.msu.ru.
РЭНСИТ | 2013 | ТОМ 5 | НОМЕР 2
13
14
RADIOELECTRONCS
ADAPTIVE OPTIMIZATION OF THE ACTIVE POLARIZATION
RADARS' CHARACTERISTICS
Oshlakov V. G.
Zuev Institute of Atmosheric Optics, Siberian Branch of Russian Academy of Sciences, http://www.iao.ru
634021 Tomsk, Russian Federation
oshlakov@iao.ru
Ilyushin Ya. A.
Lomonosov Moscow State University, Faculty of Physics, http://www.phys.msu.ru,
119991 Moscow, Russian Federation,
ilyushin@physics.msu.ru.
In the presented paper, the problem of maximization of the radar target contrast due to distinctions in the
scattering properties of the target and background is considered. Unified description of all the constructive elements of the locator in terms of Mueller matrices and Stokes vectors is formulated. The structural
diagram of the monostatic adaptive polarization locator is presented, as well as the polarimeter principal
scheme. The problem of optimal selection of the target with know polarimetric properties from the noise
is considered, as well as the problem of identification of single and group targets on the background with
known scattering properties. Scattering matrices of fresh wate and hydrosol at 440 nm wavelength are evaluated, and relative radar contrast of the hydrosol in the fresh water is estimated. The same problem is posed
and solved for aerosol in surrounding air at two wavelength 300 and 550 nm. The possibility of enhacement
of radar target contrast using available information about polarization scattering properties of the target
and background is shown.
Keywords: radar target contrast, underlying surface, probe signal, polarization characteristics, Stokes vector,
Mueller matrices, hydrosol, aerosol.
UDC 535.36+530.1
Bibliography - 16 references
RENSIT, 2013, 5(2):3-14
REFERENCES
1. Svistov VI. Radiolokatsionnye signaly i ikh obrabotka [Radar
signals and their processing]. Moscow, Sov.Radio Publ.,
1977, 447 p.
2. Davydov PS, Sosnovsky AA, Khaymovich IA.
Aviatsionnaya radiolokatsiya. Spravochnik. [Aviation radar.
Handbook.]. Davydov PS (ed.). Moscow, Transport
Publ., 1984, 116 p.
3. Varganov ME, Zinoviev YuS, Astanin LYu. Ed. Tuchkov
LT. Radar characteristics of aircraft. Moscow, Radio i svyaz'
Publ., 1985, 236 p.
4. Vakman DE. Complex signals and the uncertainty principle in
radar. Moscow, Sov. radio Publ., 1965, 304 p.
5. Oshlakov VG, Eremina T. The increase of the efficiency
of research tools of active radio and optical sensing.
Proc. SPIE Europe Remote Sens., 6744B-90, Florence,
Italy, 2007.
6. Oshlakov VG, Oshlakov VK, Eremina TA. Optimal
polarization characteristics of the devices, recording
the scattered radiation. Optika atmosfery i okeana, 2007,
20(7):635-640 (in Russ.).
7. Oshlakov VG, Andreev MI, Malykh DD. An adaptive
lidar. Proc. SPIE, Remote Sens., 74730N, Berlin, Germany,
2009.
2 НОМЕР | ТОМ 5 | 2013 | РЭНСИТ
Received 05.11.2012, revised 28.05.2013
8. Oshlakov VG. Sposob kompleksnoy lokatsii tseli [Method for
Complex Terget Location]. Patent RF no. 2416108C1:
G01S13/95. Date of Fil. 25.08.2009, Bull.11, 2011.
9. Bohren CF, Huffman DR. Absorption and Scattering of
Light by Small Particles. New York, Willey, 1983.
10. Van de Hulst HC. Light scattering by small particles. New
York, Willey, 1957.
11. Oshlakov VG. Optimum measuring the scattering
matrix. Optika atmosfery i okeana, 1992, 5(11):1190-1193
(in Russ.).
12. Oshlakov VG, Barkov YuG. Numerical analysis of the
instrumental matrix of polarization meter. Optika atmosfery
i okeana, 2002, 15(7):635-640 (in Russ.).
13. Esepkina NA, Korol'kov DV, Pariyskiy YuN.
Radioteleskopy i radiometry [Radio telescopes and
radiometers]. Moscow, Nauka Publ, 1973, 415 p.
14. Ivanov AP. Fizicheskie osnovy gidrooptiki [Physical Basis
of Hydrooptics]. Minsk, Nauka i tekhnika Publ., 1978,
504 p.
15. Kenneth JV, Edward SF. Measurement of the
Muller Matrix for Ocean Water. Appl. Opt., 1984,
23(23):4427-4439.
16. Gorchakov GI. Matrix light scattering by surface air. Izv.
AN SSSR. Fizika atmosfery i okeana, 1966, 2(6):595-605
(in Russ.).
Документ
Категория
Без категории
Просмотров
3
Размер файла
2 520 Кб
Теги
оптимизация, адаптивных, характеристика, поляризационных, локатором, активный
1/--страниц
Пожаловаться на содержимое документа