CC BY
115
УДК 396.96.6
ОПРЕДЕЛЕНИЕ ПОЛЯ И ХАРАКТЕРИСТИК РАССЕЯНИЯ РАДИОЛОКАЦИОННЫХ ОБЪЕКТОВ В СВЕРХШИРОКОПОЛОСНОМ И МОНОХРОМАТИЧЕСКОМ СЛУЧАЯХ НА ОСНОВЕ ПЛАНАРНЫХ
ИЗМЕРЕНИЙ
В.С. КАЛАШНИКОВ, А.В. КИРПАНЕВ, А.И. КОЗЛОВ, Е.Е. НЕЧАЕВ,
С.В. ВАСЕНДИН, М.Ю. ПОНОМАРЕВ
По заказу главной редакции
В статье приведены методы определения поля и характеристик рассеяния на основе планарных измерений суммарного поля, обусловленного облучающим зондом и рассеиваемым объектом. С помощью измерения на двух плоских поверхностях определяется амплитудный спектр рассеянного поля. Определяются компоненты обобщенной матрицы рассеяния для случая, когда облучающее поле обладает линейной поляризацией.
Ключевые слова: электромагнитное поле, вторичное поле, эффективная площадь рассеяния, матрица рассеяния, планарные измерения, передающий зонд, приемный зонд, поляризация.
Эффективная площадь рассеяния (ЭПР) является, как известно, одной из основных радиолокационных характеристик, от которой зависит дальность обнаружения целей, эффективность их сопровождения и наведения на цель. Еще более информативной характеристикой является обобщенная матрица рассеяния, характеризующая также явление деполяризации поля [1]. В общем случае для определения ЭПР исследуемого объекта (ИО) необходимо знать рассеянное им электромагнитное поле (ЭМП). Падающее на ИО ЭМП, как правило, считается известным. Обычно ИО по отношению к передатчику и приемнику находится в дальней зоне, т.е. падающее ЭМП представляет собой плоскую волну определенной поляризации. Для реальных ИО в СВЧ диапазоне дальняя зона может располагаться на 50, 100 и более метров, что делает затруднительным организацию измерений в закрытых помещениях.
В закрытых экранированных помещениях (компактных полигонах) наибольшее распространение получили методы определения характеристик рассеяния с помощью коллиматорных измерений [2]. Применение коллиматорного метода дает возможность непосредственно определять ЭПР ИО по результатам измерений, а также исследовать характер отраженных от объектов сигналов во временной области. Однако коллиматор позволяет определять только моноста-тические характеристики рассеяния ИО, а в некоторых случаях, для решения задач распознавания целей, необходимо знать бистатические характеристики рассеяния ИО. Кроме того, в большинстве случаев, коллиматор представляет собой зеркальную антенну с вынесенным облучателем, что приводит к принципиальной невозможности создания идеального плоского фронта в зоне размещения ИО. Диаграмма направленности сверхширокополосного облучателя на некоторых частотных областях может существенно отличаться от “столообразной”, что также негативно влияет на равномерность амплитудного распределения в рабочей области коллиматора.
Полные характеристики рассеяния ИО могут быть определены на основе амплифазометри-ческих измерений в ближней зоне. Изложение общих принципов определения характеристик рассеяния по измерениям в ближней зоне можно найти в [3; 4].
Рассмотрим случай, когда ИО размещается на двухкоординатном поворотном столе, облучается неподвижной антенной, а результирующее поле измеряется на двух плоских поверхностях в области между облучающей антенной и ИО. Первичное поле, падающее на ИО, обладает линейной поляризацией, причем плоскость поляризации перпендикулярна вертикальной стойке сканера, перемещающей приемный зонд. Конструкция стойки и крепления зонда должны обес-
печить минимальное затенение первичного поля. Предполагается, что в качестве зонда используется антенна, близкая по своим свойствам к электрическому диполю. В обобщенном виде подобная схема измерения приведена на рис. 1, где обозначены: ИО - исследуемый объект; 1 -поворотный стол; 2 - радиопоглощающий материал; ИЗ - излучающий зонд; ПЗ - приемный зонд; Г - генератор; ПУ - приемное устройство. Поворачивая ИО на двухкоординатном поворотном столе и для каждого положения ИО проводя серию планарных измерений, можно приближенно определить независимые компоненты матрицы рассеяния.
Рис. 1. Обобщенная схема измерения
Направление падения фронта поля на ИО можно определить вектором kx = kx ex + ky ey,
где k = sin J0cos jo , kyo = sin J0 sin ф0. Применяя принцип суперпозиции однородных плоских волн [5], рассеянное ИО электрическое поле в сверхширокополосной интерпретации можно определить следующим образом
і і
Е*(г,X) = к21 11В(к', о)е-]к^ешхёк'ёк^ю , (1)
о -1 -1
где г = хех + уеу + 2вг — радиус-вектор в системе координат 0ХУ2, связанной с ИО; X - время; к = 0 ; о - круговая частота; с - скорость света; к' = к'хех + к'уеу; к'х и к'у — направляющие косинусы однородных плоских волн рассеянного поля (нормированные к к пространственные частоты); В(кт', о) — векторный амплитудный спектр однородных плоских волн рассеянного поля. Так как в свободной среде йг\Е = 0, то справедливо равенство
кХ (к>) + к'уВу (к» + кВ (ко) = 0. (2)
Из (2) следует, что только две из трех компонент векторного амплитудного спектра однородных плоских волн полностью определяют рассеянное ИО поле. Учитывая последнее и используя принцип стационарной фазы [5], независимые составляющие спектра можно определить следующими равенствами
В,(к',ш) - і—к, к2е(Ях(к,к, , ш)Ах (£ , ш)+Яу(к',к. , ш)А„(£, ш)),
Г
Ву (к X, ш)-і 2р к,о к2 е’ (яу (к, к0, ш) А, (к 0, ш)+яу (к, к 0, ш) Ау (к 0, ш)),
(3)
где г — расстояние между ИО и ИЗ (рис. 1); Ах (кх , о) и Ау (кх , о) — независимые компоненты векторного амплитудного спектра однородных плоских волн падающего на ИО поля; 8^ (к', кх , о), ^у (к', кХо , о), 8-%. (к', кх , о), £у (к', кХо , о) — искомые компоненты матрицы рассеяния.
Для определения составляющих матрицы рассеяния необходимо выполнить планарные измерения составляющих напряженности электрического поля (рис. 1) для двух взаимноортогональных положений ИЗ “(1)” и “(2)” при заданном к' , соответствующим углам поворота $о и фо .
Рассмотрим один из методов, позволяющих учесть влияние первичного поля при проведении измерений с помощью зонда, имеющего широкую ДН. Предположим, что зонд по своим направленным свойствам близок к электрическому диполю. В этом случае необходимо использовать методику, позволяющую по планарным измерениям суммарного поля, выделить вторичное поле (рассеянного ИО) на фоне первичного (поля ИЗ) [6]. Применение этой методики с учетом свойств линейной поляризации ИЗ позволяет получить для элементов матрицы рассеяния
8/ (/ = х, у; ] = х, у ), имеющей в общем случае вид
{V, оуЛ Я х Я х
Ях ЯУ
VЯу Яу1
следующие значения
Я, (к, кх 0, ш) = П(к 0, ш) Л’(кX, к о, ш)/ А_х,|(кX 0, ш), я Xу (к, к 0, ш) = П(к0, ш)./™ (к, к 0, «>)/ А,1 (к0, ш), яу (к,к 0, ш) = П(к о, ш) 42> (к,к 0, ш)/ А<2'' (к 0, ш) Я у (к,к 0, ш) = й(к 0, ш) / <2) (к,к 0, ш)/ А(2) (к, , ш),
(4)
ге
ікг
/Х‘>(кX,к., ш) - (е
- (е - і'кк,,01
/«(к, к,, ш, ,02 ) - еі™ (к, к,, ш, ,01 »А',
/ учк, к 0, ш) - (е -/1'к'г01 /'‘»(к, к0, ш, ,02) - е-т,‘і2 іу"»(к, к0, ш, ,и»/4', /х2|(к,к, ш) - (е-/‘к,'г0‘ іх2)(к, к„, ш, ,02) - е-/“>02 іх2,(к, к„, ш, ,м»/4'
/ <2)(к X, к X 0, ш) - (е - і“,,01 і у-’(к X, kх 0, ш, ,02) - е - ’“',02 і !г’(к X, kх 0, ш, ,ш))/ 4'.
- („-’кк'г01 Т(2)
-’кк,,02 т(2)і
(5)
г
/«(к, к
,ш, ,0
^ 2 хтах утах к к/ к/
г г г ^2,(1)(у, у, ,01,(02), ¿, кх())еІkkxX+іккууе“іш *ёхёуЖ, (6а)
г
хтіп утіп
Т_
к к 1 2 хшах утах к ~
I<y)(kТ,~0,—^01,(02)) = — І І I Еу’(1)(хy,^01,(02),кх0УккхХ+ккуУе“•/Ш 1ахауж, (бб)
ТхУ
__лшт /шт
2
Т_
к к 1 2 хтах ушах к к/
/Х2)(~ , кт 0’ — г01,(02)) = 2Р І І І ЕХР)( х’ У’ ^01,(02) , *■, ~т0 )еМхХ+іккуУе“ •/Ш ^ЫуЖ’ (бв)
_Т хшіп Ушіп 2 Т_
к к 1 2 хшахушах
IУ2)(кТ ,кт0’ —, ^01,(02)) = — I I I Е
к к Л Шал »' шал к ~ ~
(кТ,кт0,—,г01,(02)) = — І І І Еу,(2)(x,У,г01,(02),t,кт0Уккхх+ккуУе~•/Ю 1йхйуй1. (бг)
Т
хшіп У шіп
д' = j Юшах 5іп(—шах~(z01 _г02)), ~ = л/1 _~х2 _~У2 , Юшах - наибольшая частота рассматри-рс с ''
ваемого спектра, Т - период обработки, который, как правило, равен длительности “окна” ПЗ.
У У
В подынтегральных выражениях системы (б) составляющие Еху (х, у, z01) и Еху (х, у, z02) измеряются соответственно на плоскостях 701 и 702 (рис. 1). При этом z02 _ z04 = 1шіп /4 , 1 шіп -длина волны, соответствующая —шах. Таким образом, для каждого положения ИО, определяемого $0 и Ф0 (кТ0 ), необходимо выполнять 4 серии измерений (по 2 на плоскостях 701 и 702).
Угловой сектор а, в котором при заданных $0 и Ф0 определяются компоненты матрицы рассеяния, определяется равенством
Ь _ в
а = аг^-----------------------------------------------, (7)
2 zo
где Ь _ размер стороны плоскости измерения; В _ наибольший линейный размер ИО, перпендикулярный направлению прихода волнового фронта поля ИЗ; zo _ расстояние от ПЗ до ИО.
Таким образом, предлагаемая методика позволяет определить как моно-, так и бистатиче-ские характеристики рассеяния. Вместе с тем, так как необходимо выполнение условия
2 В 2
г <-----, то применение метода возможно только для малогабаритных ИО, наибольший ли-
1 шіп
нейный размер В которых не превышает 10Ішіп. Для крупногабаритных ИО (В >10Ішіп) необходимо обеспечить облучение ИО плоским фронтом, создаваемым, например, зеркальным коллиматором. В этом случае, проводя планарные измерения в области пространства между коллиматором и ИО и используя изложенную выше методику, можно также определять как моно-, так и бистатические характеристики рассеяния. Однако в этом случае вертикальное перемещение зонда сканером должно быть реализовано таким образом, чтобы удавалось избежать существенных искажений поля коллиматора.
Рассмотрим далее монохроматический случай, когда используются два зонда, обладающих возможностью независимого перемещения по плоской поверхности в области ближней зоны ИО (рис. 2). Обобщенная схема такого сканера показана на рис. 2. Математические основы такого способа определения характеристик рассеяния изложены в [3]. Ниже приведены соотношения для случая планарного зондирования и с учетом влияния на приемный зонд поля передающего зонда.
Примем следующие условия: направленные свойства приемного зонда не учитываются, т.к. предполагается, что в качестве зонда используется электрический диполь; зонд, работающий на
2
прием, измеряет результирующее поле, обусловленное рассеянным полем ИО и падающим на ИО полем передающего зонда; направленные свойства передающего зонда считаются известными. С передающим зондом связана правая система координат. На рис. 2 обозначено: (х0у0) -координаты приемного зонда в системе координат 0ХУ2, связанной с ИО; (х1,у1) - координаты передающего зонда в системе 0ХУ2. В системе координат, связанной с ИО, напряженность электрического поля, измеряемого зондом, определяется равенством
И
у
1 1
( 1 Іт ► і і ¡7і
1 х ] 1 1 1 і і І Г
N о _ ы" о
х
ИЗ
(хь
Ґ1 п П
! (х0,
где
Рис. 2. Обобщенная схема двухзондового сканера
Е = Е
1 1
+ Е Р = к 2 ІІ Б* (к, х1, У,)е _И‘хх^кку
х_ ДкУу_
$к х$к у +
_1 _1 11
+ к 2 ІІ Л(~Т, х1, У1)е _ккх
x_jkk'yУ_jkkZzг!к> к
йк'хйк'у.
_1 _1
Бх(1)(~Т,7) = (е~М^01 іхі)(к5т, 7, z02) _ е~]кк'^02іх2)(~т,7,z01))/А',
БУ(1) (~Т, 7) = (е --^о11 У1 (кт, 7, z 02) _ е _кк^021 У2) (к, 7, z 01))/ А',
(1)
V
х шах У шах
(кТ,Z01(02)) = І І ^х1^7ъ ^^У^0 + '/ккуУ0dxodyo.
хшіп ушіп
к х шах У шах
1 У1)(кт, rl, z01(02)) = І І R, z01(02))е
хшіп Ушіп хшах У шах
jkkХxo + іккУУ0
dxodyo.
І?)(К, 7Ъ z01(02)) = І І Ех2)( 'Іь 7Ь z01(02))ejkkxХ0 +МуУ0 ^0^0.
хшіп Ушіп хшах ушах
1 Ї\К, 7Ь ^1(02)) = І І Е(у)(7Ь z01(02))ejkkxХ0 + МуУ0 ^0^0.
хшіп Ушіп
К= Кех + к;еу, г0 = х0ех + У0еу, 71 = х1ех + У1еу.
(8)
(9)
(10а)
(10б)
(10в)
(10г)
у
Амплитудный спектр однородных плоских волн рассеянного ИО поля зависит от искомой матрицы рассеяния и амплитудного спектра падающего на ИО поля и определяется для положений “(1)”-“(2)” ИЗ следующими системами
11
Б
(1)
(~т,Г) = к2 ІІ(^(~т,~т)Ах1}(~т,71) + ^(~т,~т)¿ун-т,Rl)dk'xdk'y,
_1_1
11
бх2) (~т, А) = к2 ІІ (¿х (~т, ~т) А х2)(~т, Г) + ¿у (~т, ~т) Ау2 (к т, А^к^ку,
_1 _1
(11)
в';>(кт, а) = к2 ІІ (¿у (кт, кт) Ат(кт, А)+¿у (к'т, кт) А»^, А^к,.,
_1 _1
Б^2)(кГ, А) = к2 І І(^У (кт, кт) Ах2)(кт, 71) + ¿УУ -Д) А^2)(кт, АМ~Д,.
(12)
Подставляя в (11)-(12) А(кт,Г) = А (к т)е ]ккхХ ]ккуУ1+}ккгг\, Где А(~) - амплитудный спектр однородных плоских волн ИЗ и применяя двумерное преобразование Фурье, получим следующие системы
¿х (кт, кт) лх1}(кт) + ^ (кт, кт) АУ\кт) = ^\К, кт),
¿х (~т, ~т) ах2)(~т)+¿х (~т, ~т) аУ2)(~т) = 42)(кт, ~т),
(13)
(кт, кт) лххкт)+¿у (кт, кт) лу\кт) = ^УЧК, кт), ¿У (~т, ~т) лх2)(кт) + (~т, кт) л;у2)(кт) = 3 (у)(~т, кт),
(14)
где
^х1 (кт, -т ) =
_ ]ккгг1 х шах У шах к ~ ~
е ГГ г>(1)(к/ А )еІ хх1 + }]ку у1
т' /1/2 2 4к р
І І в™]', ає
йх1ф 1.
хшіп шіп
(15а)
е_]кк2І1 хшах ушах
3х2)(*т',-т) = ‘—ГТ І 1 вх2)(]т,А)еДк'х‘*МуУ1 0x^1-
4к р2 : :
хшіп У шіп
(15б)
к к _ Ікк2г1 х шах У шах к
3У1)(к,кт) = —^ І І БУ\к',А)еккхХ1 + іккуУ1
4к 2 р2
0x1^1.
хшіп Ушіп
(15в)
_ Ікк71л хшах ушах к ~ ~
зУ2|(*т,-т)=І ш в<2|(кт,а)‘ 1+ уУ1Огі^і.
4к 2 р 2
х шіп у шіп
(15 г)
1 1
1_1
Из систем (13)-(14) определяются независимые компоненты матрицы рассеяния
s; (k'„ k,) = (^(k’J,) - jf(k;, k,) Ay«(kz))/д (16а)
s; (~t, ~)=( jf(~;, ~,) Af(~,) - jf(~;, ~,) a<«(!,))/д (16б)
SX (~t,~) = (Jy;Ц(~'Д) Af(~,) - j;2,(~;, ~) A“)(~))/д (16в)
sy (~;, ) = ( j® (~;, к,) a™ (~,) - j; 1 (~;, ) а;\ї , »/ д, (16г)
где А = A® (~) AУ2 (~) - 41} (~t) аХ2) (~ ).
Если ИЗ обладает линейной поляризацией, то компоненты матрицы рассеяния определяются следующим образом
SX (~, ~) = У" (~, ~) / А<'> (~,). (17а)
Sy (~, ~) = 42)(~, ~)/А<'>(~). (176)
SX (~, ~) = ./“(I,', ~)/А<'>(~). (17е)
Sy (~, ~t) = ~t)/А<2)(~,). (17г)
Таким образом, для определения независимых компонент матрицы рассеяния при использовании в качестве ПЗ антенны, близкой по направленным свойствам к электрическому диполю, необходимо для каждого положения ИЗ (xiyi) при определенной ориентации его плоскости поляризации (1 или 2) выполнять 2 серии планарных измерений на плоскостях, расстояние между которыми не должно быть кратным 1/2. Если предположить, что количество положений ПЗ и ИЗ в плоскости одинаково и равно N, то полное количество измерений, необходимое для
определения характеристик рассеяния, составит 4N , что осложняет практическую реализацию метода [1]. Последнее побуждает в качестве инструмента реализации изложенной идеи использовать плоскую решетку или независимо перемещаемые по горизонтали линейные вертикальные решетки [7]. Очевидным преимуществом применения плоской решетки является быстрота получения первичной информации о рассеянном поле, а к недостаткам можно отнести следующие: сложность обеспечения сверхширокополосных измерений и идентичности приемопередающих каналов; взаимное влияние между облучателями; высокая стоимость. Приведенные в работе методы измерения и обработки позволяют избежать этих недостатков ценой продолжительных по времени измерений.
ЛИТЕРАТУРА
1. Козлов А.И., Логвин А.И., Сарычев В.А. Поляризационная структура радиолокационных сигналов / Поляризация радиоволн. - М.: Радиотехника, 2005. - Т. 1.
2. Балабуха Н.П., Зубов А.С., Солосин В.С. Компактный полигон для измерения рассеивающих свойств объектов и параметров антенн // Антенны, 2008. - Вып.6 (133). - С. 59-80.
3. Вилькоцкий М.А., Гринчук А.П., Громыко А.В., Михасенко А.Г. Разработка голографического метода измерения характеристик рассеяния объекта // Тезисы докл. V Всесоюз. конф. “Метрологическое обеспечение антенных измерений”. - Ереван, 1990. - С. 277-280.
4. Barry J.C., Charles E. Ryan Jr. Near-field scattering measurements for determiniig complex target RCS // IEEE Trans. Ant. and Prop. 1989. vol. AP-37. N5. P. 576-585.
5. Kerns D.M. Plane-Wave Scattering-Matrix Theory of Antennas and Antenna-Antenna Interactions. Washington: U.S. Government Printing Office, NBS Monograph 162, 1981.
6. Кирпанев А.В. Определение электромагнитных полей излучающих систем с учетом влияния отраженных полей и направленных свойств зонда при планарных измерениях // Вопросы радиоэлектроники, серия РЛТ. - 1995. - Вып. 1. - С. 72-76.
7. Кирпанев А.В. Определение характеристик рассеяния объектов по измерениям в ближней зоне // Физика волновых процессов и радиотехнические системы. - Т. 7. - № 3. - 2004. - С.45-48.
DETERMINATION OF THE FIELD AND CHARACTERISTICS OF THE DISSIPATION RADAR OBJECT IN SUPERBROADBAND AND MONOCHROMATIC EVENTS ON BASE OF THE PLANAR
MEASUREMENTS
Kalashnikov V.S., Kirpanjev A.V., Kozlov A.N., Nechaev E.E., Vasendin S.V., Ponomarjev M.U.
Methods of the determination of the field and characteristics of the dissipation are brought in article. The base of the methods form the planar measurements of the total field. This field presents itself amount of the field, created by radiator and floor of the dissipation. The Measurements were conducted on two flat surfaces. Amplitude spectrum of the absent-minded field was defined. The elements of the generalised matrix of the scattering are defined for linear polarization.
Key words: electromagnetic field, secondary electromagnetic field, secondary field, efficient dissipation surface, matrix of the dissipation, planar measurements, transmitting probe, receiving probe, polarization.
Сведения об авторах
Калашников Вадим Серафимович, 1940 г.р., окончил ЛИАП (1963), доктор технических наук, профессор, начальник отдела ОАО “ВНИИРА”, автор более 60 научных работ, область научных интересов - антенны и устройства СВЧ, исследование антенных обтекателей и укрытий, антенные измерения.
Кирпанев Алексей Владимирович, 1962 г., окончил ЛИАП (1986), доктор технических наук, старший научный сотрудник, начальник сектора ОАО “ВНИИРА”, автор 77 научных работ, область научных интересов - электромагнитная совместимость, антенные измерения, исследование радиолокационных характеристик рассеяния.
Козлов Анатолий Иванович, 1939 г.р., окончил МФТИ (1962), заслуженный деятель науки и техники РФ, академик Академии транспорта РФ и Международной академии информатизации, профессор, доктор физико-математических наук, Соросовский профессор, заведующий кафедрой технической эксплуатации радиоэлектронных систем воздушного транспорта МГТУ Г А, автор более 300 научных работ, область научных интересов - радиофизика, радиополяриметрия, радиолокация.
Нечаев Евгений Евгеньевич, 1952 г., окончил МАИ (1974), профессор, доктор технических наук, проректор по научной работе, заведующий кафедрой УВД МГТУ ГА, автор свыше 100 научных работ, область научных интересов - радионавигация, антенная техника, управление воздушным движением.
Васендин Сергей Васильевич, 1986 г., окончил ЛЭТИ (2010), аспирант ОАО “ВНИИРА”, автор 2 научных работ, область научных интересов - обработка и кодирование сигналов, системный анализ, исследование радиолокационных характеристик рассеяния.
Пономарев Максим Юрьевич, 1980 г., окончил СПбГУ (2005), инженер 2 категории ОАО “ВНИИРА”, автор 5 научных работ, область научных интересов - обработка сигналов, электромагнитная совместимость, антенные измерения, исследование характеристик рассеяния.
