CC BY
35
УДК 621.396.96
ФОРМИРОВАНИЕ ПАРАМЕТРОВ ПОЛЯРИЗАЦИОННОГО СПЕКЛА ПРИ КОГЕРЕНТНОМ РАССЕЯНИИ КАК ИНТЕРФЕРЕНЦИОННЫЙ ПРОЦЕСС
А.И. КОЗЛОВ, В.Н. ТАТАРИНОВ, С.В. ТАТАРИНОВ, Н.Н. КРИВИН
По заказу редакционной коллегии
В настоящей работе продемонстрировано, что формирование поляризационно-энергетического спекла при рассеянии волн многоточечными сложными радарными объектами (МСРО) представляет собой интерференционный процесс. Такое представление позволяет определить как автокорреляционные функции параметров Стокса рассеянного поля, так и пространственный спектр мощности МСРО.
Ключевые слова: поляризация радиоволн, спекл, рассеяние, интерференция.
Введение
Поляризационно-энергетические параметры электромагнитного поля при когерентном рассеянии волн МСРО (спекл) тесно связаны с процессом формирования когерентного изображения этого объекта. Настоящая работа посвящена анализу углового распределения поляризационно-энергетического спекла как интерференционного процесса с теоретическим и экспериментальным исследованием автокорреляционных функций и пространственных энергетических спектров многоточечных случайных радиолокационных объектов.
1. Формирование углового распределения поляризационно-энергетических параметров как интерференционный процесс
Используя модель процесса рассеяния электромагнитных волн сложным объектом в виде интеграла Стрэттона-Чу [1], мы можем представить поле, рассеянное таким объектом, как сумму волн, рассеянных элементарными рассеивателями («блестящими» точками), составляющими сложный объект. Соответствие этого подхода физике рассеяния подтверждают рис. 1а, б. Так, на рис. 1 а изображен военный корабль, а на рис. 1 б представлен результат моделирования процесса рассеяния этим объектом с использованием метода масштабного моделирования [2]. На данном рисунке четко видны блестящие точки, обусловленные как элементами искусственного объекта (корабль), так и подстилающей поверхностью.
а б
Рис. 1:
а - фотография военного корабля; б - результат масштабного моделирования
Рассмотрение результатов моделирования позволяет построить достаточно простую, но репрезентативную модель когерентного рассеяния волн сложным искусственным объектом. Как известно [1-5], при рассеянии точечным объектом, обладающим матрицей рассеяния 5^ (і, к = 1,2), поле
в дальней зоне определяется как Е$(0) = exp(-у^кЛ^) $I Е0 exp(]2кХ70)/л/4яЛ0, где Х1 - координаты элементарных отражателей, а 0 - позиционный угол объекта (рис. 2).
fxv] L X. X,
ч Ч \ \ V ¡¡К- \ ^ б\ * с \ 'V\V \ ‘А X ц \\ 1 /- / ^^ —м
Рис. 2. Г еометрия рассеяния многоточечным объектом
Рассмотрим теперь процесс формирования поляризационно-энергетических параметров в дальней зоне при рассеянии многоточечным МСРО как интерференционный процесс, полагая без уменьшения общности, что элементы объекта расположены на одной прямой (рис. 2). В качестве примера найдем поле, рассеянное 4-точечным объектом, в дальней зоне (точка Q )
для случая совпадающих линейных поляризаций передачи и приема
4
ES (0) = - exp( •j2kRo)Eо exp(- j2kXi0). Определим
V4pR0 i=1
мгновенное распределение мощности
ю
рассеянного поля в пространстве как функцию угла 0
2
P(0) = \Es (0)|2 = A
4
exp( - j2kXi0)
i=1
4
4 4
cos(2kdmn0), (n > m). (1)
i=1
m=1n=2
Таким образом, мгновенное распределение мощности рассеянного поля в зависимости от угла 0 формируется в виде суммы ЭПР элементарных рассеивателей плюс 6 пространственных гармонических колебаний. Нетрудно видеть, что каждое из гармонических колебаний обусловлено эффектом интерференции между волнами, рассеянными одной из пар элементарных рассеивателей объекта. Число этих пар определяется биномиальным коэффициентом М!
, где М - общее число рассеивателей, а N = 2 есть число рассеивателей в ин-
rN _ wy =
'M N!(M - N)!
терференционной паре. Для случая M = 4 имеем C4 = 6 . Угловая функция отклика рассматриваемого объекта включает 6 пространственных гармоник, как это следует из выражения (1), определяемых величинами dmn = Xm - Xn, m = 1,3; n = 2,4; n > m, которые представляют собой пространственный разнос элементов интерференционных пар. Пространственные гармоники y]smsn cos(2kdmn0) характеризуются полной фазой y(0) = 2kd0 = 2p2d0/1, производная которой —dy = 2d = есть пространственная частота, имеющая размерность Rad-1. Период
2p d0 1
гармоники TSp = 1/fSp =1 /2d имеет размерность Rad.
Таким образом, полная мощность, рассеянная МСРО, есть сумма интерференционных картин, формируемых набором элементарных двухточечных интерферометров. Тогда выражение (1) можно переписать в виде
M C
p(0) = X + 2XVs sk cos(2kdik 0),
m =1
где C = CM - биномиальный коэффициент.
Рассмотрим теперь вопрос о формировании параметров Стокса поля, рассеянного МСРО, как интерференционного процесса. В работе [3] было показано, что параметры Стокса поля, рассеянного МСРО, имеют вид
So(0) = SS + So + 2yfss\[Щcos(2k/0+0,5hi), S3(0) = S3 + S 3 + 2т[#!тЩ,{цЬ sin(2k/0 - 0,5^2).
Отсюда следует, что пространственные гармоники cos(2k/0 ± h) имеют амплитуды
■yJsO\jNab и tJs3-\[s3■sjDab . Здесь величины Nab, Dab представляют собой близость
(удаленность) состояний поляризации элементов одной из интерференционных пар, составляющих МСРО.
В этом случае угловое распределение параметров Стокса МСРО примет вид суммы обобщенных законов интерференции, формируемых элементарными двухточечными интерферометрами
So(0) = X So" + 2^ SoiSokJNk cos(2kdik 0 + h*),
m=1 1
M C
Ss(0) = X S3‘ + 2XJSA^ cos(2kd,k 0+hfi),
т=1 1
где С = СМ - биномиальный коэффициент.
Поскольку как амплитуды, так и фазы пространственных гармоник случайны, дальнейший анализ должен быть статистическим.
2. Автокорреляционные функции параметров Стокса. Пространственный спектр
Определим теоретический вид АКФ (автокорреляционной функции) параметра Стокса 5з(0) углового распределения рассеянного поля. Прежде всего, мы удалим нестационарную
м
составляющую ^ 53т , используя метод усреднения скользящим окном. После устранения нет = 1
стационарности и нормирования запишем случайную функцию 5з (0) в виде
вз(0)=£л/в;ксо8(2кё1к0+л1к) . (2)
1
Полагая случайную функцию (2) стационарной, найдем её автокорреляционную функцию как В8(Д0) = М{5(01)• 5(02)}=М\ ^ со8[2М^0 + ^]со8[2Ы£(0 + А0) + 'п]1 (3)
[^=1 Ь=1 [
Здесь амплитуды л/Ь и фазы ^ есть случайные величины, характеризуемые двумерной
плотностью распределения Щ2(л/Ь, ^) ; М - оператор усреднения; Д0=01 -02 . Полагая, что случайные амплитуды и фазы независимы, запишем двумерную плотность вероятностей произведением одномерных функций Щ2 (4Ь, ^) = Щ (4Ь)Щ (^) .
Условие ортогональности J cos(WSjp 0) cos(W^P 0)d0 :
Перепишем соотношение (з) в виде
В5 (Д0) = М|^=(л/Д^) соъ\1Ык0 + ^]соб[2Ы^ (0 + Д0) + л]|. (4)
Учитывая, что реализация 5з(0) есть функция случайных аргументов л/Ь и 0, найдем АКФ этой реализации как среднее значение функции двух случайных величин
¥ ¥
M{у(ХЪ x2)} = J JУ(ХЪ x2)W2(ХЪ x2)dx1dx2.
— ¥ —¥
Используя данное выражение, найдем АКФ (3) как
C ¥ ¥
BS ( D0) = X J J (iDN f cos[2kdN 0 + h] cos [2kdN (0 +D0) + h]W2(VD, h)d (4D )dh. (5)
N = 1 — ¥ — ¥
Для вычисления двойного интеграла в выражении (5) воспользуемся условием, полагая также, что случайная фаза распределена равномерно на интервале (—p;p), т.е. W(h) = 1/2p.
Распределение случайной амплитуды VD во всех случаях будет односторонним. Представим
C
(5) в виде BS (D0) = X (Л + / 2 + /3).
N=1
Вычисление первого интеграла дает , ¥ p
Ix = 0,5 J J (DN) cos(2kdN A0)W (VD )d(VD )dh = 0,5(DN) cos(2kdN Л0),
о -p
2p
где (Pn) есть среднее значение поляризационной дистанции, которое представляет собой среднее по статистическому ансамблю случайных значений (Pn) для всех пространственных гармоник fSp = 2dN /1.
Второй и третий интегралы равны нулю.
Таким образом, теоретическая форма АКФ углового распределения третьего параметра Стокса имеет вид
BS (Л0) = £ cos(2kdNЛ0). (6)
n=1 2
Каждое из слагаемых суммы (6) есть АКФ отдельной пространственной гармоники Sn (0), имеющей случайную амплитуду ypN и случайную фазу hN
(Pn)
BSN (Л0) = cos(2kdN Л0). (7)
Нетрудно видеть, что АКФ случайной реализации параметра Стокса есть сумма индивидуальных АКФ всех пространственных гармоник, которые формируют случайную реализацию
C
Bs (Л0) = X Bsn (Л0).
N=1
Пространственный спектр мощности может быть найден путем преобразования Фурье над АКФ (5). Спектр мощности изолированной пространственной гармоники может быть получен преобразованием Фурье соотношения (7)
P(Wsp ) = J Bsn (Л0) exp(- jWspЛ00d(Л0) = 0,5(Pn - WNp ) + d(WSp + WNp )],
где О5р = 2п/$р = 2%2ё /1 есть пространственная частота. Спектральные линии расположены на
расстояниях ± О^р от начала координат. Эта пространственная частота отвечает пространственному разносу между интерференционной парой рассеивателей. Интенсивность спектральной линии определяется поляризационной дистанцией между состояниями поляризации элементов пары. Полный пространственный спектр случайного углового отклика имеет вид
Р(0 5р) = 0,5 £ Ь )[б(-П N) +6(0 &)].
N=1
Такая форма обусловлена дискретной структурой МСРО.
3. Результаты экспериментов
Результаты экспериментального исследования АКФ и пространственного спектра поляризационно-угловой функции отклика вращающегося МСРО (гусеничный вездеход) приведены на
рис. 3, 4. На рис. 3 изображены АКФ в угловом интервале ±200 к борту объекта (пунктир) и к корме объекта (сплошная линия). Измерения для этих направлений позволяют определить различия как в АКФ, так и в пространственных спектрах, обусловленные различной протяженностью объекта по длине и ширине. Кроме того, из рис. 4 видно, что пространственный спектр имеет двухмодовую форму. Это подтверждает ранее сделанное [4] теоретическое предположение о принципе эквивалентности «в среднем».
Рис. 3. АКФ поляризационного отклика МСРО
Рис. 4. Средние спектры мощности МСРО
Заключение
Результаты исследования подтверждают, что формирование поляризационно-энергетических параметров поля при рассеянии МСРО представляет собой интерференционный процесс.
ЛИТЕРАТУРА
1. Stratton J.A., Chu L.J. Diffraction theory of electromagnetic waves. Physical Review, v. 56, pp. 308-316.
2. Штагер Е.А. Методы расчёта радиолокационных характеристик объектов, находящихся вблизи неровной земной или морской поверхности // Зарубежная радиоэлектроника. - 1994. - № 4, 5. - С. 22-40.
3. Козлов А.И., Татаринов В.Н., Татаринов С.В., Кривин Н.Н. Рассеяние электромагнитных волн двухточечными радиолокационными объектами // Статья в настоящем Вестнике.
4. Tatarinov V.N., Tatarinov S.V., L.P. Ligthart The use of an equivalence principle «on the average» for a statement of definition of a random complex radar object model. Proc. of the MIKON’2000, Wrozlaw, Poland. Vol. 2, pp. 12-17.
5. Козлов А.И., Сарычев В.А., Логвин А.И. Поляризация радиоволн. - М.: Радиотехника. - 2006. - Т. 1.
POLARIZATION SPECKLE PARAMETERS FORMATION AS THE INTERFERENCE PROCESS AT THE COHERENT SCATTERING
Kozlov A.I., Tatarinov V.N., Tatarinov S.V., Krivin N.N.
In this paper is demonstrated that polarization-energetical speckle formation at the complex multipoint radar objects scattering is an interference process. Such point of view allows to determine Stocks parameter autocorrelation function of scattered waves and space power complex multipoint radar objects spectrum.
Key words: polarization, scattering, speckle, radio wave, interferention.
Сведения об авторах
Козлов Анатолий Иванович, 1939 г.р., окончил МФТИ (1962), заслуженный деятель науки и техники РФ, академик Академии транспорта РФ и Международной академии информатизации, профессор, доктор физико-математических наук, Соросовский профессор, заведующий кафедрой технической эксплуатации радиоэлектронных систем воздушного транспорта МГТУ ГА, автор более 300 научных работ, область научных интересов - радиофизика, радиополяриметрия, радиолокация.
Татаринов Виктор Николаевич, 1941 г.р., окончил ТУСУР (1964), доктор технических наук, профессор, действительный член Академии электромагнетизма (Массачусетс, США), заведующий кафедрой конструирования и производства радиоаппаратуры ТУСУР, автор более 200 научных работ, область научных интересов - теория когерентности и поляризации электромагнитного поля, статистическая радиофизика, рассеяние волн сложными объектами, поляризационная радиолокация.
Татаринов Сергей Викторович, 1969 г.р., окончил ТУСУР (1994), кандидат технических наук, доцент кафедры конструирования и производства РЭА ТУСУР, автор более 70 научных работ, область научных интересов - статистическая теория поляризации при рассеянии волн сложными объектами.
Кривин Николай Николаевич, 1985 г.р., окончил ТУСУР (2007), аспирант кафедры конструирования и производства РЭА ТУСУР, автор 7 научных работ, область научных интересов - теория поляризационного контраста малоразмерных объектов на подстилающей поверхности.
