CC BY
85
УДК 621.396.98
Радиоэлектроника и системы связи
ТРАССИРОВКА ПРЕЛОМЛЕННЫХ ЛУЧЕЙ Д.В. Асотов, Б.В. Матвеев, В.Б. Авдеев, В.А. Середа
В статье рассматриваются вопросы трассировки лучей при расчёте электромагнитного поля в лучевом приближении. Решаемая задача - трассировка преломлённых лучей. Предложен алгоритм трассировки двухмерной задачи, сводящий решение к решению уравнений четвёртой степени. Разобраны методы Кардано и Декарта-Эйлера для решения уравнений 3й и 4й степени. Получены аналитические формулы для расчёта. Приведён способ сведения трехмерной задачи к двумерной
Ключевые слова: трассировка лучей, геометрическая оптика, преломление, передача
Введение
Решения задач распространения радиоволн актуальны как в инженерной практике, так и в научной деятельности. Существующие программные продукты, в частности Wireless Insite [1] и FEKO [2], используют лучевые методы. Они основаны на приближениях геометрической оптике (ГО) и геометрической теории дифракции. Постулаты ГО обобщены в [3, с. 11-14, 4]. Основное преимущество ГО перед прочими численными методами заключается в уменьшении необходимых вычислений, а следовательно снижению затрат времени на моделирование.
В общем случае моделирование методом ГО можно разбить на два условно независимых этапа: трассировка лучей и расчёт параметров поля (напряжённости, мощности и т. д.) основываясь на известных траекториях. В специальной литературе основной акцент сделан на определение параметров поля, а вопрос трассировки лучей практически не рассматривается. Однако для моделирования поля, особенно для автоматизированного моделирования с использованием ГО [1-2], трассировка лучей необходима. Поэтому данная теория также получила развитие.
Одна из актуальных задач в данной области -трассировка преломлённых лучей. В известных системах автоматизированного моделирования задача трассировки преломленных лучей решается асимптотически. Так в программе Wireless Insite преломление лучей в тонких слоях (в частности в стенах зданий) не учитывается. Траекторией проходящего через стену луча - прямая линия. Необходимые поправки в расчёты вводятся заменой коэффициента преломления. В программе FEKO закон преломления электромагнитных волн, в частности преломления лучей в слоистом диэлектрике упрощается [5]. Аналитическое решение данной задачи приведено в работе [7]. В данной статье будет рассмотрено аналитическое решение, альтернативное ранее опубликованному, но
целью будет получение алгоритма наиболее подходящего для автоматизированного вычисления.
Алгоритм трассировки
Задача трассировки лучей на границе раздела 2х сред формулируется следующим образом: определить координаты точки преломления при известных электрических и магнитных параметрах сред, координатах передающей и приёмной антенны относительно границы раздела сред. Для описания границ будем использовать плоские поверхности. В геометрической двухмерной интерпретации задача сводится к следующему: определитьХт (илиХк) при известных Х, Нт, Нк и коэффициенте преломления п12 (см. рис. а).
sin 0 n2
■=— = n1 n
sin0
sin 0 sin 0J
% = xtJxï+h
XD
x^x2 + H2
(1)
(2)
Приравнивая правые части (1) и (2), избавляясь от радикалов, приводя к общему знаменателю и производя замену Хк = X - Хт получаем уравнение 4 степени:
mX4 - 2XmXT + (H + X2m)x2T -- 2 XH2 n22 Xt + H X n22 = 0
(3)
где т = п^2 -1 , Н = Н2Тп122 - Н\ .Кроме того,
большой интерес представляет задача трассировки лучей при прохождении через тонкий диэлектрический слой (рис. 1 б), в частности стены и другие строительные конструкции, тогда уравнение (2) преобразуется в
Асотов Дмитрий Валериевич - ВГТУ, аспирант, e-mail: [email protected], тел. (473) 243-76-65 Матвеев Борис Васильевич - ВГТУ, канд. техн. наук, доцент, тел. 8-960-138-45-61
Авдеев Владимир Борисович - ВГТУ, д-р техн. наук, профессор, тел. 8-903-65-45-520
Середа Владимир Александрович - ВГТУ, студент, e-mail: [email protected], тел. 8-920-444-64-15
sin 0 sin 0,
X^Xl + H = XrtJXS + H
xj X2 + H2
xs^x1 + HR
(4)
в результате преобразований системы уравнений (1) и (4) и осуществления замены Х8 = X - Хт - Хя и Хя = Хт• Нк/Нт получаем уравнение
R
шк2Х\ - 2ХшкХ'т + (Н1 + X2ш)х2т -
(5)
где к = 1 +—-
Н'
Н =(Нт + нК )2 - н
б)
Преломление луча а) на границе раздела, б) 2х сред в тонком слое
Уравнения (3) и (5) отличаются лишь коэффициентами при переменной Хт и имеют аналитические решения: либо решение Феррари, либо решение Декарта-Эйлера [6, с. 48-49].В общем виде приведённое уравнение 4й степени записывается как
У4 + аУ3 + ЬУ2 + еУ + а = 0.
(6)
Тогда решение уравнения получается из решения вспомогательного уравнения
У3 + а'У 2 + Ь'У + е' = 0,
<7 <7 <7
коэффициенты которого определяются по формулам
а =■
3а2 - Ь 16 2
За4 а 2Ь + ае ^ = 256 Ив Т ;
, ,,, 1 ( а3 аЬ , е = а Ь — I----+ е | .
8 I 8 2
(8)
(9)
(10)
Для получения решения уравнения (6) достаточно найти один корень уравнения (7). Наиболее простой способ - использовать формулу Кардано [5, с. 47-48]:
р = Ь'-
9 = е+-
а
2а'3 а' Ь'
27 3
У =3
9 К
+3
- — + , 1-Р1 +¡—I +
а
7
(11) (12)
(13)
Р I + ¡9
Используя данное решение, корни уравнения (6) определяются как
У + а' У = ± -±у
2
У
_<7_
2
а' _ а - 4аЬ + 8е а
(14)
16Л/27У9
+ а
4
Из четырёх полученных корней необходимо выбрать один корень, удовлетворяющий физическому смыслу задачи, т.е. лежащий в пределах от 0 до X.
Вспомогательные алгоритмы
Для сравнения уравнений (3) и (4) с уравнением четвёртой степени в общем виде (6) приведём его коэффициенты (см. таблицу).
Чтобы свести произвольную трехмерную задачу к изложенной выше двухмерной необходимо вычислить величины Нт, НК, X (Н5 для второго случая) на основе координат точек Тх (хТх, уТх, гТх) и Кх (хКх, уКх, 2р.), и коэффициентов уравнения (уравнений) плоскости (плоскостей). Уравнение плоскости удобно представить в виде
ах ■ х + а ■ у + аг ■ г + а = 0.
(15)
где {ах, ау, а2} = а - единичные нормальный
вектор к рассматриваемой плоскости, а - расстояние от плоскости до начала системы координат (с положительным или отрицательным знаком в зависимо-
НТ ХП2кХт + н2 х 2и2 = 0
2
3
2
2
сти от расположения плоскости). Координаты точек пересечения плоскостей и перпендикуляров к ним, опущенных из заданной точки
x' > f a a a
x У z
У = a У -a x 0
v az 0 -a
Y1
(
- d a ■ x - a ■ y
y x s
\
(16)
С помощью данной формулы можно определить координаты точек Тх и Кх' (проекций точек Тх и Кх на плоскость отражения). Расстояния между точками Тх и Кх'вычисляется по формуле
X = TX Rx' =
= V(XTX' - XRx' )2 + (Утх' - УшУ + (ZTx' - ZRf ) 2
(17)
Коэффициенты приведённых уравнений в задачах трассировки преломлённых лучей
Коэффициент Задача
а Б
a -2 X -2 X / h
b (H + X2 m)/m (H1 + X2 m)/mh2
c - 2НгТ Xnlz / m - НгТ Xn22 / mh
d H2 X X2 / m H T X 2n IJ mh2
TxTx
Та же формула позволяет вычислить HT = и HR =
RxRX
Для решения второй задачи
(см. рис. б) также необходимо вычислить Н8 = - а2\, а точку Кх и расстояние НК определить для второй
плоскости с нормальным вектором {а2х, а2 , а2г} = а2
и расстоянием а2 до начала системы координат.
Заключение
Полученное решение позволяет вычислить координаты точек преломления, следовательно, решать задачу трассировки преломлённых лучей. Его пре-
имущество перед другими известными решениями отсутствие векторных уравнений, которые тяжело поддаются автоматизации.
Данный алгоритм, может быть использован совместно с аналогичными алгоритмами [8-9] для трассировки лучей в системах автоматизированного моделирования процессов распространения радиоволн.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (Грант РФФИ №13-08-97538-р-центр-а.)
Литература
1. Официальный сайт компании Remcom. {URL: www. remcom. com}.
2. Официальный сайт программы FEKO. {URL: www.feko.info}.
3. Боровиков В.А., Кинбер Б.Е. Геометрическая теория дифракции. - М.: Связь, 1978 - 248 с.
4. C. Balanis. Advanced Engineering Electromagnetics. New York: Wiley, 1989.
5. Методы геометрической оптики. Учебные материалы компании «Родник» {URL: http://www.rodnik.ru}.
6. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров под ред. И.Г. Араманович, Л.Я. Цлаф. М.: «Наука», 1970 - 720 с.
7. Панычев А.И. Дубинская И.В. Синтез лучевой траектории проникновения сигналов WLAN в смежные помещения // Известия ЮФУ. Технические науки. - 2013. - № 5 (142). - С. 116-122.
8. Панычев А.И. Алгоритм трехмерной трассировки радиоволн локальной беспроводной сети // Известия ЮФУ. Технические науки. - 2012. - № 11 (136). С. 31-41.
9. Авдеев В.Б. Моделирование процессов распространения радиоволн в условиях городской застройки [Текст] / В.Б. Авдеев, Д.В. Асотов, В.А. Сладких, Б.В. Матвеев // Вестник Воронежского государственного технического университета - 2013. -Т.9. - № 2. - С.4-7.
Воронежский государственный технический университет
REFRACTION RAYTRASING ALGORITHM
x ; i az ■x - ax ■ Z ,
D.V. Asotov, B.V. Matveev, V.B. Avdeev, V.A. Sereda
In article ray tracing questions are considered in case of calculation of an electromagnetic field in ray approximation. The solved task is trace of refracted beams. The algorithm of trace of the two-dimensional task, bringing together the decision to the solution of biquadratic equations is offered. Methods Cardano and Descartes-Euler for the solution of the equations 3h and 4h levels are sorted. Analytical formulas for calculation are received. The method of convergence of the three-dimensional task is brought to the sorted two-dimensional task
Key words :ray tracing, geometry optic, refraction
