close

Вход

Забыли?

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

?

Численный метод определения неоднородной комплексной диэлектрической проницаемости плоской поверхности объектов по поляризационной структуре поля отраженной электромагнитной волны.

код для вставкиСкачать
НАУЧНЫЙ ВЕСТНИК МГТУ ГА
2012
№ 179
УДК 621.396
ЧИСЛЕННЫЙ МЕТОД ОПРЕДЕЛЕНИЯ НЕОДНОРОДНОЙ КОМПЛЕКСНОЙ
ДИЭЛЕКТРИЧЕСКОЙ ПРОНИЦАЕМОСТИ ПЛОСКОЙ ПОВЕРХНОСТИ
ОБЪЕКТОВ ПО ПОЛЯРИЗАЦИОННОЙ СТРУКТУРЕ ПОЛЯ
ОТРАЖЕННОЙ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ
А.И. КОЗЛОВ, В.Ю. МАСЛОВ
Предлагается алгоритм численного метода определения электрических параметров неоднородных объектов на
основе информации об изменении поляризационной структуры отраженной волны. Метод использует лучевые
представления полей рассеяния, базирующиеся на принципе Гюйгенса-Френеля.
Ключевые слова: определение комплексной диэлектрической проницаемости, поляризация электромагнитных волн.
Существенным практическим приложением поляризационных свойств радиоволн является
определение электрических параметров (диэлектрической проницаемости ε0 и проводимости σ0)
материальных сред на основе информации об изменении поляризационной структуры рассеянной волны с эллиптической или круговой поляризацией [1].
Разработка приближенных численных методов радиополяризационного анализа рассеянных электромагнитных волн на поверхностях, проводимость которых конечна, с целью определения электрических параметров ε0 и σ0 представляет существенный практический интерес.
Для объекта, размеры и главные радиусы кривизны поверхности которого значительно
больше длины падающей волны λ, характер и степень деполяризации зависят от электрических
параметров объекта ε0 и σ0, угла ее облучения. Так как поляризация падающей волны, ее длина
λ и угол облучения являются величинами известными, параметры ε0 и σ0 можно определить,
зная поляризацию отраженной волны.
Вопрос о возможности определения электрических параметров ε0 и σ0 объекта поляризационным методом впервые был рассмотрен в работе [2]. Там же получено выражение для определения комплексной диэлектрической проницаемости ε0 однородного объекта, имеющего плоскую поверхность.
В работе [3] предложен базирующийся на принципе Гюйгенса-Френеля алгоритм метода
численного решения обратной задачи рассеяния на объекте, который позволяет найти комплексную функцию рассеяния объекта. Используя этот метод, можно определить неоднородные
значения электрических параметров ε0 и σ0 плоской поверхности исследуемого объекта.
Пусть на плоскую поверхность (рис. 1) с конечной проводимостью под некоторым углом θ
r
(исключая значения угла θ = 0 и θ = π 2 ) падает монохроматическая волна E п , например, с
круговой поляризацией. Электрические свойства поверхности характеризуются комплексной
диэлектрической проницаемостью
ε 0 = ε 0 + j 60λσ 0 ,
(1)
где ε 0 - диэлектрическая проницаемость; σ 0 - проводимость.
Полагая, что длина волны λ значительно меньше линейных размеров поверхности, и пренебрегая дифракцией на границах поверхности (краевым эффектом), можно записать выражеr
ние для волны E о , отраженной в зеркальном направлении, в виде суммы двух ортогональных
r
r
линейно поляризованных компонент E1о и E 2о
r
r
r
r
r
r
E o = E1o + E 2o = i cos θ + k sin θ F 1 (ε 0 ,θ ) + ij F 2 (ε 0 ,θ ) eikρ ,
(2)
[(
)
]
141
Численный метод определения неоднородной комплексной …
где F 1 (ε 0 ,θ ) и F 2 (ε 0 ,θ ) - коэффициенты Френеля, которые равны:
F 1 (ε 0 , θ ) =
ε 0 cos θ − ε 0 − sin 2 θ
- для случая параллельного поля;
ε 0 cos θ + ε 0 − sin 2 θ
F 2 (ε 0 , θ ) =
cos θ − ε 0 − sin 2 θ
cos θ + ε 0 − sin 2 θ
- для случая поперечного поля.
(3)
(4)
Разрешая соотношения (3) и (4) относительно ε 0 и произведя необходимые преобразования,
получается выражение для комплексной диэлектрической проницаемости в виде [2]

2

 1 − ip 
 tg 2θ  sin 2 θ ,

 1 + ip 
ε 0 = 1 + 

(5)
где p = E 1о = − i F 1 , p = p .
E 2о
(6)
F2
z
e1mn
e2mn
Е1о
Е1п
Е2о
Е2п
ρ
Плоскость
наблюдения
θ
θ
k
i Поверхность объекта
x
j
y
Рис. 1. Отражение поляризованной волны
от плоской поверхности с конечной проводимостью
Если на поверхность объекта (рис. 1) падает локально плоская электромагнитная волна с
круговой поляризацией, то в соответствии с принципом Гюйгенса - Френеля в дискретных
точках наблюдения (m,n), комплексная амплитуда напряжённости электрического поля отраr
r
r
женной волны e m, n = e 1m, n + e 2 m, n может быть вычислена по формуле
e v m, n =
∫
D( x, y )
F v x, y
eikρ
ρ
dx dy ,
(7)
где v = 1, 2 ; F v x , y – коэффициенты Френеля; ρ – расстояние между рассеивающей точкой на
поверхности объекта и точкой, где определяется величина отраженного поля; k = 2π λ – волновое число; λ – длина волны; D(x,y) – контур объекта. В каждой точке (m,n) плоскости, расположенной перпендикулярно направлению распространения отраженной волны (рис. 1), измеряются ортогональные составляющие отраженного поля e1 m, n и e 2 m, n . Расстояние между со-
142
А.И. Козлов, В.Ю. Маслов
седними точками равно d. Массивы измеренных составляющих отраженного поля e1 m , n и e 2 m, n
имеют размер N × N .
Нахождение распределения коэффициентов Френеля F 1 x , y и F 2 x , y на поверхности объекта
сводится к решению двух независимых систем линейных алгебраических уравнений вида [3]
ΦP1 = G1 и ΦP2 = G 2 .
(8)
Матрицы G v (v = 1, 2) размером N × 1 образуются из массивов значений комплексных со2
ставляющих отраженного поля e1 m, n , e 2 m, n и имеют вид
(
Gv = ev1,1 ev1, 2 L ev1, N
Матрица Φ размером N 2 × N 2
 φ 1,1
φ 11,,12
 1,1
 φ 1,1
φ 11,, 22
 1, 2
M
 M
1, 2
 φ 1,1
φ 1, N
Φ =  11,,1N
φ
φ 12,,21
2
,
1
 1,1
 φ 2, 2 φ 12,,22

M
 M
1,1
1, 2
φ
φ N ,N
 N ,N
ikρmx ,,yn
имеет элементы φ m,n = e
x, y
(9)
N ,N
L φ 1,1 
N ,N
L φ 1, 2 

M 
N ,N
L φ 1, N 

N ,N
L φ 2,1 

N ,N
L φ 2, 2 

M 
L φ NN ,, NN 

(10)
L
φ 11,,1N
φ 12,,11
φ 12,,12
L
φ 11,, 2N
φ 12,,21
φ 12,,22
M
M
M
L
φ 11,, NN
φ 12,,N1
φ 12,,N2
φ 12,,N1
φ 22,,11
φ 22,,12
φ 12,,N2
φ 22,,12
φ 2N, 2, N
L
L
L
M
M
M
φ 1N, N, N
φ 2N,1, N
φ 2N, 2, N
, где ρ mx,, yn - расстояние между точками (x,y) и (m,n).
Искомые матрицы коэффициентов Френеля имеют вид
(
Pv = p v1,1
)T
ev 2,1 ev 2,2 L ev N , N .
p v1, 2 L p v 1, N
p v 2,1
)
p v 2, 2 L p v N , N T .
(11)
Эти матрицы затем следует преобразовать в двумерные матрицы размером N × N
 p

 v 1,1 L p v 1, N 
(12)
P0 v =  M
M .


 p v N ,1 L p v N , N 


Используя формулы (6) и (5) для известного угла θ из массивов коэффициентов Френеля
(12) рассчитываются массивы комплексных диэлектрических проницаемостей ε 0v m, n .
Рис. 2. Массив распределения модуля комплексной амплитуды
напряжённости электрического поля отраженной волны e1 m, n (N = 40)
Численный метод определения неоднородной комплексной …
143
Рис. 3. Массив распределения модуля комплексной амплитуды
напряжённости электрического поля отраженной волны e 2 m, n (N = 40)
На рис. 2 и рис. 3 изображены массивы значений модуля комплексной амплитуды напряжённости ортогональных составляющих отраженного поля e1 m, n и e 2 m, n .
Рис. 4. Массив распределения диэлектрической проницаемости ε 0
по поверхности объекта (N = 40)
Рис. 5. Массив распределения диэлектрической проводимости σ 0
по поверхности объекта (N = 40)
144
А.И. Козлов, В.Ю. Маслов
Результаты решения систем уравнений (8) и последующего расчета по формулам (12), (6) и (5)
позволяют определить значения диэлектрической проницаемости ε 0 и проводимости σ 0 в различных точках поверхности исследуемого объекта.
Массивы распределения диэлектрической проницаемости ε 0 и проводимости σ 0 , характеризующие электрические свойства поверхности исследуемого объекта, изображены на рис. 4 и
рис. 5.
ЛИТЕРАТУРА
1. Козлов А.И., Логвин А.И., Сарычев В.А. Поляризация радиоволн. Поляризационная структура радиолокационных сигналов. - М.: Радиотехника, 2005.
2. Поздняк С.И., Мелитицкий В.А. Введение в статистическую теорию поляризации. - М.: Сов. радио, 1974.
3. Козлов А.И., Маслов В.Ю. Численный метод решения трехмерной обратной задачи рассеяния электромагнитных волн на препятствии // Научный Вестник МГТУ ГА, серия Радиофизика и радиотехника. - 2012. - № 180.
NUMERICAL METHOD OF THE DETERMINATION COMPLEX PERMEABILITY TO FLAT SURFACE
OBJECT ON POLARIZATION STRUCTURE OF THE FIELD REFLECTED ELECTROMAGNETIC WAVE
Kozlov A.I., Maslov V.Ju.
The algorithm of the numerical method of the determination dielectric parameter object is offered on base of information on change polarization structures reflected waves. The method uses the beam presentations by flap of the dissipation, basing on Gygens-Frenel principle.
Key words: the determination of complex permeability, the polarization of the electromagnetic waves.
Сведения об авторах
Козлов Анатолий Иванович, 1939 г.р., окончил МФТИ (1962), заслуженный деятель науки и техники
РФ, академик Российской академии транспорта и Международной академии информатизации, доктор физико-математических наук, профессор, автор более 300 научных работ, область научных интересов – радиофизика, радиолокация, радиополяриметрия, дистанционное зондирование окружающей среды.
Маслов Виктор Юрьевич, 1945 г.р., окончил МГУ им М.В. Ломоносова (1968), доктор технических
наук, профессор МГТУРЭА, автор более 70 научных работ, область научных интересов - радиофизика,
радиолокация, радиополяриметрия.
1/--страниц
Пожаловаться на содержимое документа