CC BY
107
УДК 537.876:551.51
К. В. Скобельцын, Г. М. Тептин, А. Р. Низамеев, Е. С. Нефедьев, И. Р. Низамеев
ИССЛЕДОВАНИЕ РАСПРОСТРАНЕНИЯ ДЕЦИМЕТРОВЫХ РАДИОВОЛН МЕТОДОМ ИМИТАЦИИ
Ключевые слова: имитационное моделирование, дециметровые радиоволны.
В представленной работе описана имитационная модель распространения радиоволн в неоднородной среде, которая основана на методе трассировки радиолучей с использованием полигональной модели представления среды распространения.
Keywords: imitation model, UHF radio waves.
The implementation model of radio waves propagation in inhomogeneous media which based on the method of tracing radio waves with used polygonal model is described.
Имитационная модель распространения дециметровых радиоволн
Введение. Одной из важнейших задач радиофизики является прогнозирование
распространения радиоволн в неоднородных средах. Вследствие невозможности решить эту задачу аналитически, было разработано множество приближенных эмпирических моделей, которые имели ограниченный диапазон применения [1,2].
Современные технологии позволяют производить численные эксперименты, которые дают возможность учитывать множество факторов при моделировании процесса.
В статье освещается численный метод моделирования распространения радиоволн под названием «трассировка лучей». Данный метод успешно применяется в оптике [3] и акустике [4,5]. Трассировка лучей позволяет учитывать неограниченное количество факторов, но предъявляет высокие требования к вычислительной мощности. Также трассировка лучей производится при помощи оптического приближения движения лучей, но данный недостаток можно компенсировать при помощи использования метода Монте-Карло [2,6]. Проведено множество исследований, которые доказывают работоспособность трассировки лучей для прогноза распространения радиоволн [7,8]. Описанные ниже методы представления среды и распространяющейся волны позволяют ассимилировать результаты других моделей или экспериментальных данных.
Основной сложностью реализации трассировки лучей оставались слабые вычислительные способности компьютеров, но современные технологии позволяют проводить расчеты трассировки большого количества лучей и хранения результатов.
Представление среды. При моделировании используется полигональное представление среды. В трехмерном моделировании полигональное представление является наиболее распространенным методом, так как позволяет моделировать формы предметов любой сложности [9]. Каждый трехмерный объект представляет сетку "склеенных" между собой элементарных фигур - полигонов. Размеры полигонов не ограничены, поэтому возможна экономия памяти при моделировании путем применения полигонов переменного размера. Также применение
полигональной модели позволяет производить расчеты различных масштабов. В представленной модели мы используем треугольники в качестве фигур для полигонов. На рис. 1 продемонстрированы примеры использования полигонов для представления трехмерных поверхностей.
Рис. 1 - Представление ландшафта при помощи набора треугольных полигонов
Выбор в пользу полигональной модели объясняется тем, что все больше исследований проводится для систем компьютерного зрения и распознавания образов, которые также применяют полигональную модель [10], что позволяет ассимилировать результаты систем анализа среды и получить более точные входные данные модели. Объекты моделируемой среды, составленные из полигонов, должны быть замкнутыми, так как внутри каждого объекта определены свои радиофизические параметры. В данной статье мы применяем результаты исследования индекса рефракции дециметровых радиоволн, проведенные в [11]. В качестве параметра, характеризующего вещество, мы использовали индекс рефракции.
Представление радиоволны. Фронт волны определяется набором точек, которые имеют начальные координаты, совпадающие с координатами источника сигнала. Траектория движения каждой точки подчиняется законам геометрической оптики и может
интерпретироваться как луч радиоволны. Каждый луч проходит свой путь, взаимодействуя с
преградами во время распространения. Совокупность точек с одинаковым временем генерации является волновым фронтом.
Движение точки волнового фронта дискретизировано, и определяется следующими параметрическими выражениями:
Хп = Хп-1 + ¡п--^1
Уп = Уп-1 + , гп = Ип-1 + кп-1Ь1
(1)
где Хп - координата X точки волнового фронта при итерации п, уп - координата У точки волнового фронта при итерации п, гп - координата 7 точки волнового фронта при итерации п, ¡п-1 - составляющая вектора следования вдоль оси абсцисс, VI - составляющая вектора следования вдоль оси ординат, кп-1 -составляющая вектора следования вдоль оси аппликат, ^ - это шаг времени модели.
При взаимодействии с препятствием луч может отразиться или отклониться от первоначального направления движения: изменятся параметры вектора следования луча ¡, ^ к. Критерии выбора происходящего явления и его реализация будут описаны ниже. Точность выходных данных модели зависит от количества моделируемых частиц. Количество частиц в моделируемой среде ограничено вычислительной мощностью и объемом хранилища данных.
Источник сигнала. В модели источник сигнала испускает лучи в разных направлениях в соответствии с заданной для него диаграммой направленности. Количество источников сигнала не ограничено. Количество лучей в каждом направлении пропорционально мощности сигнала источника в заданном направлении. На рис. 2 показана трассировка лучей точечного источника со сферической диаграммой направленности в однородном пространстве.
Рис. 2 - Излучение точечного источника со сферической диаграммой направленности в однородном пространстве
На рис. 2 можно видеть, как фронт волны равномерно распространяется по всем заданным направлениям.
Распространение лучей в пространстве. При генерации моделируемой среды в модели формируется набор элементов, состоящих из неограниченного
количества полигонов каждый. Каждый элемент проходит предобработку: для каждого полигона вычисляется нормаль и направляется от центра элемента в бесконечность. Нормаль вычисляется при помощи векторного произведения векторов, которые сформированы из двух сторон треугольника. Пусть вершины полигона описываются точками А(Хд,Уа,7д), Б(Хв,Ув,7в), (Хс,Ус,7с), тогда мы можем получить векторы:
а1 = (Хв - Ха, Уб - Уа, 7б - 7а) а2 = (Хс - Ха, Ус - Уа, 7с - 7а)
тогда вектор нормали будет равен:
п = ([а-| х а2]) / п
После вычисления нормали, формируется уравнение плоскости полигона вида:
Ах + Бу + с-г + й = 0, где й = - (аХА + ЬУа + о7а).
(2)
После вычисления уравнения плоскости для полигона, данные о полигоне отправляются в базу данных. При выполнении моделирования все элементы и их полигоны копируются из базы данных в оперативную память вычислительной станции.
Вывод: полигональная модель
представления среды позволяет формировать среду любой сложности для трассировки лучей. Данный подход дополняет способы трассировки радиолучей, опубликованных в источниках [1-3].
Взаимодействие луча со средой. При распространении луч пересекает множество элементов моделируемой среды. Каждый элемент состоит из полигонов, то взаимодействие происходит при пересечении луча с полигоном. Так как распространение луча описывается системой формул (1), то пересечение луча и полигона вычисляется следующим образом: выполняется поиск ближайшего полигона, который находится на пути следования волнового фронта. Поиск выполняется следующим образом: подставим значения координат из системы уравнений (1) в (2) и получим выражение, которое выглядит следующим образом:
А' (Хм-1+г' гстолкн )+ В' ^УN-1+ ^' г столкн )+ (3)
+ °-1+ к''столкн)+ С = 0
Решая уравнение (3) относительно 1столкн, получаем следующее выражение:
А' хм-1 + В' Ум-1 + С' ZN-1 + п
А' г + В' ] + С' к
(4)
Если 1с-
меньше нуля или его значение
равно бесконечности, то пересечения с плоскостью полигона нет. Если пересечение существует, то
г
столкн
выполняется проверка принадлежности точки полигону.
После проверки всех полигонов выбирается ближайший, а после этого сравнивается 1столкн с шагом ^ и если 1столкн < d1, то модель считает, что столкновение произошло. После этого происходит отражение или преломление.
Преломление лучей подчиняется закону Снеллиуса. Так как при решении задачи распространения луча мы используем векторную формулировку, то более удобным и эффективным является следующее выражение, которое представляет упрощенное выражение из [13,14]:
r 2:
n 2
rv
+1 -1
(л • n)->
где у2 - вектор преломленного луча, п - вектор
падающего луча, п - вектор нормали к поверхности, п - коэффициент преломления среды 1, п2 -коэффициент преломления среды 2.
Данная формулировка позволяет производить расчет преломления наиболее эффективным способом, так как не содержит даже тригонометрических функций. Как известно, при прохождении луча через границу раздела двух сред возможно как преломление, так и отражение. В данной модели учитывается явление полного внутреннего отражения, которое происходит, если выполняется следующее условие, которое выводится из закона Снелиуса:
тогда
П1 •sin 01 = П2 •sin 02'
sin 02 = П •sin 01 • П2
В случаях, когда п-|-8ш01 > п2, преломление луча невозможно, и происходит полное внутреннее отражение, реализуемое по принципу, представленному ниже. Формулировка отражения для радиолучей формулируется следующим образом:
Г 2 = Т\- 2(Л'п)п, где у2 - вектор отраженного луча, - вектор
падающего луча п - вектор нормали к плоскости полигона.
Пример использования
Для того, чтобы можно было наглядно показать работу трассировки, мы использовали модель атмосферы, которая представлена в виде элементов, имеющих формы параллелепипеда, находящихся друг над другом и имеющих значения коэффициента преломления, который бы повторял вертикальный профиль коэффициента преломления реальной атмосферы, представленный в работах [11,12]. Сформировав набор полигонов, мы провели трассировку лучей из точечного источника с
полусферической диаграммой направленности в указанной синтезированной среде, расположенного на высоте 20 м над идеально ровной поверхностью Земли. Вследствие неоднородности структуры атмосферы волновой фронт радиосигнала с поверхности Земли искажается. Трассировка лучей позволяет определить, как происходит это искажение. На рис. 3 мы можем наблюдать искажения фронта радиоволны при полусферической диаграмме направленности. Для наглядности рекомендуем рис. 3 сравнить с рис. 1.
Рис. 3 - Искажение волнового фронта радиоволны в неоднородной атмосфере
Данные рисунки несут лишь качественный характер, так как для получения количественных данных требуется больше информации о топографии, антропогенных факторах и др.
Заключение
Разработанная модель позволяет производить моделирование распространения радиоволн с заданной точностью. Применение полигональной модели представления среды позволяет проводить моделирование в среде любой геометрической формы. На представленном примере показана возможность ассимиляции данных из внешних источников.
Литература
3
Tetsuro, I. Ray-tracing Method for Estimating Radio Propagation Using Genetic Algorithm / I. Tetsuro // NTT DoCoMo Technical Journal. - 2007. - V.9. - No.3. - P.20-27.
Molisch, A. Wireless Communications / Wiley, 2011. -882p.
Chan, Y-C. A hybrid ray-tracing and radiosity method for calculating radiation transport and illuminance distribution in spaces with Venetian blinds / Y-C. Chan, A. Tzempelikos // Solar Energy. - 2012. - V.86.- P.3109-3124.
4. Smith, A. Numerical ray tracing in media involving continuous and discrete refractiVe boundaries / A. Smith, M. Goldberg, E.S.K. Liu // Ultrasonic Imaging. - 1980. - V.2. -P.291-301.
5. Schmitz, V. 3D ray tracing in austenite materials / V. Schmitz, F. Walte, S.V. Chakhlov // NDT & E International. - 1999. - V.32. - P.201-213.
6. Edward, R. Freniere Edge diffraction in Monte Carlo ray tracing / R. Edward, G. Groot Gregory, R. A. Hassler // Optical Design and Analysis Software - 1999. - V.3780. -P.151-157.
7. David, I.L., Indoor Radio Channel Propagation Modelling by Ray Tracing Techniques / Ph. D. Thesis, Diss., University of Edinburgh, 1994.
8. McKown, J. Ray tracing as a design tool for radio networks / J. W. McKown, R.L., Jr. Hamilton // Network, IEEE. - 1991. -V.5. - P.27-30.
9. Smith, C., On Vertex-Vertex Systems and Their Use in Geometric and Biological Modelling / Ph. D. dissertation, University of Calgary, 2006.
10. D.T. Lee, D. Computing the visibility polygon from an edge / D.T. Lee, A.K. Lin // Computer Vision, Graphics, and Image Processing. - 1986. -V.34. - No.34. - P.1-19.
11. Низамеев, А.Р. Сравнительный анализ индекса рефракции дециметровых радиоволн по радиоизмерениям ГЛОНАСС и с помощью радиозонда / А.Р. Низамеев, Г.М.Тептин, А.П. Шлычков // Научно-технические ведомости СПбГПУ. Физико-математические науки. -2012. - Т.158. - №4. - P.144-150.
12. Низамеев, А.Р. Трехмерная структура индекса рефракции радиоволн в тропосфере по измерениям сети приемных
станций GPS-ГЛОНАСС / А.Р. Низамеев, Е.С. Нефедьев, И.Р. Низамеев, Г.М. Тептин // Вестник Казанского технологического университета. - 2012. -№9. - С.27-31.
13. Schlick, C. An Inexpensive BRDF Model for Physically-Based Rendering / C. Schlick // Comput. Graph. Forum -1994. - P.233-246.
14. Низамеев, И.Р. Определение эффективной поверхности наноструктурированного платинового катализатора мембранно-электродных блоков топливного элемента при помощи потенциодинамического метода / И.Р. Низамеев, Е.С. Нефедьев, И.Э. Исмаев, М.К. Кадиров // Вестник Казанского технологического университета. - 2012. -Т.15 - №3. - С.126-130.
© К. В. Скобельцын - асп. каф. радиоастрономии КФУ, [email protected]; Г. М. Тептин - д-р физ.-мат. наук, проф. той же кафедры, [email protected]; А. Р. Низамеев - асп. той же кафедры, [email protected]; Е. С. Нефедьев - д-р хим. наук, зав. каф. физики КНИТУ, [email protected]; И. Р. Низамеев - канд. хим. наук, м.н.с. лаб. ЭХС ИОФХ им. А.Е. Арбузова КазНЦ РАН, доц. каф. физики КНИТУ, [email protected].
