DOI: 10.12731/2227-930X-2016-2-79-87
ПРИНЯТИЕ РЕШЕНИЯ О ФОРМЕ ДИФРАКЦИОННЫХ СТРУКТУР В ЗАДАЧЕ РАССЕЯНИЯ РАДИОВОЛН
Преображенский А.П., Чопоров О.Н.
В работе рассматривается задача рассеяния электромагнитных волн на различных дифракционных структурах. Решение задачи рассеяния осуществляется на основе метода интегральных уравнений. По диаграммам обратного рассеяния при разных частотах падающей волны принимается решение о форме исследуемого объекта.
Ключевые слова: дифракция; интегральное уравнение; рассеяние радиоволн; радиосвязь; дифракционная структура; метод моментов; радиоизлучение.
THE DECISION OF FORM FOR DIFFRACTIVE STRUCTURES IN THE PROBLEM OF SCATTERING OF RADIO WAVES
Preobrazhensky A.P., Choporov O.N.
This paper considers the problem of scattering of electromagnetic waves in different diffraction structures. The solution of the scattering problem is based on the method of integral equations. On diagrams of backscattering at various frequencies of the incident wave, the decision about the form of the object is carried out.
Keywords: diffraction; integral equation; scattering of radio waves; telecommunication; diffraction structure; method of moments; radioirradiation.
Для того, чтобы проектировать и оценивать эффективность современных радиолокационных систем требуется знать априорные характеристики рассеяния исследуемых объектов [1-4]. Вследствие развития техники и способов обработки радиолокационных сигналов возникают возможности для увеличения числа и точности измеряемых параметров поля, которое рассеивается от объектов при рассмотрении процессов радиолокации. Следует отметить, что актуальной является задача, связанная с разработкой новых подходов и алгоритмов, направленных на обнаружение и селекцию радиолокационных сигналов [5-8].
Информацию, связанную с состояниемпринимаемых электромагнитных волн можно использовать для того, чтобы обнаруживать и проводить селекцию радиосигналов по соответствующим признакам [7].
Идентификация объектов в радиолокационной сфере может осуществляться по спектрам отраженных сигналов и значениям эффективной поверхности рассеяния (ЭПР). По удаленным (точечным) рассеивателям ЭПР рассматривается как мера их видимости для радаров. Так как точных решений задач определения ЭПР для произвольных тел в общем случае не найдено, были разработаны несколько способов ([9-11]) для того, чтобы ее приближенным образом решать: аналитические, численные и комбинированные.
Когда размеры тел являются сопоставимыми с длиной падающей плоской монохромной волны 1 - то мы говорим о зоне резонансного рассеивания здесь эффективными являются метод конечных элементов (МКЭ) в частотной области и метод моментов (МоМ).
В данной работе мы рассматриваем задачу об определении формы дифракционной структуры на основе информации о угловых зависимостях ЭПР. Примеры рассматриваемых структур приведены на рис. 1.
Рассмотрим решение задачи в рамках метода интегральных уравнений. Запишем интегральное уравнение Фредгольма второго рода для электрического тока [12] при учете граничных условий для поверхностей идеально проводящих объектов (рис. 1)
1/г; = 2п X НгУг; +—п X Д/г; X роТвсЬ!, (1)
2п й
где О = ехр(-]кг)/т - является трехмерной функцией Грина, относящейся к свободному пространству, она будет являться решением уравнения Гельмгольца когда 5-образный источник; ^ - обозначает поверхность объектов; п - показывает внешнюю нормаль к поверхности объекта для точки наблюдения;
в) г)
Рис. 1
= [п х Н - является поверхностной плотностью эквивалентного электрического тока; Н\г) = хЯ^ + уЩ + 2Н[ - является вектором первичной плоской электромагнитной волны, которая вертикально поляризована.
Для того, чтобы определить неизвестные поверхностные электрические токи на поверхности анализируемых структур проводится решение интегрального уравнения (1) на основе метода моментов [12] таким образом. Определяем систему базисных функций, относительно которых происходит разложение электрических токов на поверхности структур. Для нашего случая в качестве базисных функций мы выбираем кусочно-постоянные функции. Потом требуется ввести систему пробных функций. Пробными функциями в нашем случае были взяты 5-функции Дирака. Применение такой системы с точки зрения физики показывает, что выполнение граничных условий по магнитному полю происходит не для всех поверхностей объектов S, а только для дискретного числа точек. Анализ показывает [12], что когда выбирается достаточное количество точек, по которым происходит удовлетворение граничных условий, мы можем получить высокую точность решений, когда осуществляется моделирование рассеяния электромагнитных волн на различных объектах.
Как результат применения метода моментов интегральное уравнение (1) превращается в систему уравнений:
(2)
здесь - являются компонентами плотности поверхност-
ных электрических токов.
В указанной выше системе для вычисления блоков в матрице и используются следующие выражения:
1
и* Л X"
и,. и„ и, =
и* 3 _ г _ к
= ТГ" \((Пу)т(8Га(1у )тп + (Пг)т(ё^ЛшЖ - 5„
2 ' '
1
1
(и^)^ =-— КПу)т(gradx)nmdsn; (С/^ = -— \(n2)m(gradx)lmdsn;
(Uyx)„m = -^-\(nx)m(grad' ^ds'^ (U ^ =~ ¡(njjgrad' 2ns 2ns
(uyy)mn = — \((nx)Jgradx)mn + (nz)m(gradz)mn)dsn -5
mn
271 S
(uzx)mn =--!:r\(nx)m(grad'z)mndsn, (Uzy)tm = -^-\(ny)m(grad'z)nmdsn,
Zn S ZJli
(V*)wrn ^-!:-\((nx)m(grad'x)mn + (Пу )rn(Srad'y )тп№п ~ ^mn
2 ns у у
где m, n = 1, ..., N, где N - количество точек разбиения поверхности объектов.
Здесь dmn - является символом Кронекера,
grad Gmn = -гт exp(-jkrmn ) =
^ (4)
У grad, )„ + Кgrady)mn + k(grad'z ^
где rmn = - рассчитывается в виде единичного вектора, направил |
ленного из точки источника в точку наблюдения [12].
Для вектора-столбца свободных членов записываются такие выражения:
(RJm =2((ny)JH[)m-(nJJH^U, (Ry)m =-2 ((nJJHU-fnJJH'JJ; (5)
(RJm =2((пх)т(Н'у)т ~(пу)т(Щ)т). После того, как решена приведенная система уравнений (2) происходит вычисление рассеянного электромагнитного поля для дальней зоны, которое связано с определенным электрическим током JS(r) таким образом[12]:
Н- (г) = J^exp(-jkr) х г, (6)
4 tv „
где k - является волновым числом;
=4-71-10"'Гн/м, е0 =8,85-10"12Ф/м - представляют собой магнитную и электрическую постоянную;
г - является радиус-вектором точки источника;
г - является радиус-вектором точки наблюдения по дальней зоне.
Проводится расчет ЭПР объекта на основе такой формулы [13]
Ит 2|нмГг/ о = 4кг-Ь (7)
|н'Уг;|
Определение формы объекта может быть полезным в различных автоматизированных системах [14, 15]. Для того, чтобы принять решение о форме объекта, мы будем рассматривать случаи, когда сначала падает электромагнитная волна с частотой £[, а затем с частотой причем их значение отличается на порядок.
Характеристика рассеяния объекта, приведенного на рис 1а, в обоих случаях будет иметь характерный максимум при нулевом угле наблюдения.
Объекты, приведенные на рис 1б и 1г, по характеристикам рассеяния отличаются тем, что во втором случае их значения отличаются на большие величины, что позволяет принять решение, что объектом является полая прямоугольная структура. При относительно небольших изменениях (не более 5-6 дБ) для сектора углов 0°-20° принимается решение, что рассматривается структура, приведенная на рис. 1б.
Для того, чтобы принять решение о том, что мы рассматриваем объект, приведенный на рис. 1в, необходимо провести анализ изменений характеристик рассеяния в секторе углов наблюдения 0°-30°. Если значения изменяются не более, чем на 2-3 дБ, то мы принимаем решение.
Вывод
В работе рассмотрены вопросы принятия решения о форме объектов на основе анализа характеристик рассеяния при различных частотах падающей волны.
Список литературы
1. Головинов С О., Хромых А.А. Проблемы управления системами мобильной связи // Вестник Воронежского института высоких технологий. 2012. № 9. С. 13-14.
2. Казаков Е.Н. Разработка и программная реализации алгоритма оценки уровня сигнала в сети wi-fi // Моделирование, оптимизация и информационные технологии. 2016. № 1. С. 13.
3. Пекшев Г.А., Скляр А.Г. Некоторые свойства лучевых методов, используемых для анализа распространения электромагнитных волн // В мире научных открытий. 2015. № 12. С. 17-23.
4. Пекшев Г.А., Скляр А.Г. Проблемы управления рассеянными электромагнитными полями // В мире научных открытий. 2015. № 12. С. 24-30.
5. Максимова А.А., Юрочкин А.Г. Методы исследования характеристик рассеяния электромагнитных волн объектами // Вестник Воронежского института высоких технологий. 2016. № 16. С. 53-56.
6. Канавин С.В., Лукьянов А.С. Перспективы применения систем мобильного широкополосного доступа в сетях подвижной радиосвязи на основе стандартов mobile WIMAX и LTE // Вестник Воронежского института высоких технологий. 2016. № 16. С. 79-82.
7. Кульнева Е.Ю., Гащенко И.А. О характеристиках, влияющих на моделирование радиотехнических устройств // Современные наукоемкие технологии. 2014. № 5-2. С. 50.
8. Глотова Т.В. Применение гибридного метода для расчета характеристик рассеяния объектов над шероховатой поверхностью // Моделирование, оптимизация и информационные технологии. 2016. № 1. С. 11.
9. Баранов А.В. Проблемы функционирования mesh-сетей // Вестник Воронежского института высоких технологий. 2012. № 9. С.49-50.
10. Милошенко О.В. Методы оценки характеристик распространения радиоволн в системах подвижной радиосвязи // Вестник Воронежского института высоких технологий. 2012. № 9. С. 60-62.
11. Ерасов С.В. Оптимизационные процессы в электродинамических задачах // Вестник Воронежского института высоких технологий. 2013. № 10. С. 20-26.
12. Щербатых С.С. Метод интегральных уравнений как основной способ анализа в САПР антенн // Моделирование, оптимизация и информационные технологии. 2016. № 1. С. 10.
13. Справочник по радиолокации / Под ред. Сколника М. М.: Сов. радио, 1976. Т.1. 455 с.
14. Мишин Я.А. О системах автоматизированного проектирования в беспроводных сетях // Вестник Воронежского института высоких технологий. 2013. № 10. С. 153-156.
15. Болучевская О.А., Горбенко О.Н. Свойства методов оценки характеристик рассеяния электромагнитных волн // Моделирование, оптимизация и информационные технологии. 2013. № 3. С. 4.
References
1. Golovinov S.O., Khromykh A.A. Vestnik Voronezhskogo instituta vysokikh tekhnologiy. 2012. № 9, рр. 13-14.
2. Kazakov E.N. Modelirovanie, optimizatsiya i informatsionnye tekhnologii. 2016. № 1. P. 13.
3. Pekshev G.A., Sklyar A.G. Vmire nauchnykh otkrytiy. 2015. № 12, рр.17-23.
4. Pekshev G.A., Sklyar A.G. Vmire nauchnykh otkrytiy. 2015. № 12, рр. 24-30.
5. Maksimova A.A., Yurochkin A.G. Vestnik Voronezhskogo instituta vysokikh tekhnologiy. 2016. № 16, рр. 53-56.
6. Kanavin S.V., Luk'yanov A.S. Vestnik Voronezhskogo instituta vysokikh tekhnologiy. 2016. № 16, рр. 79-82.
7. Kul'neva E.Yu., Gashchenko I.A. Sovremennye naukoemkie tekh-nologii. 2014. № 5-2. P. 50.
8. Glotova T.V. Modelirovanie, optimizatsiya i informatsionnye tekh-nologii. 2016. № 1. P. 11.
9. Baranov A.V. Vestnik Voronezhskogo instituta vysokikh tekh-nologiy. 2012. № 9, рр. 49-50.
10.Miloshenko O.V. Vestnik Voronezhskogo instituta vysokikh tekh-nologiy. 2012. № 9, рр. 60-62.
11.Erasov S.V. Vestnik Voronezhskogo instituta vysokikh tekhnologiy. 2013. № 10, рр. 20-26.
12.Shcherbatykh S.S. Modelirovanie, optimizatsiya i informatsionnye
tekhnologii. 2016. № 1. P. 10. 13.Spravochnikpo radiolokatsii [Reference radar] / Skolnik M. (ed.). M.: Sov. ra-dio, 1976. t.1. 455 p.
14.Mishin Ya.A. Vestnik Voronezhskogo instituta vysokikh tekhnologiy. 2013. № 10, рр. 153-156.
15.Boluchevskaya O.A., Gorbenko O.N. Modelirovanie, optimizatsiya i informatsionnye tekhnologii. 2013. № 3. P. 4.
данные об авторах
Преображенский А.П., чопоров О.Н.
Воронежский институт высоких технологий
ул. Ленина, 73А, Воронеж, Воронежская обл., 394043,
Российская Федерация
komkovvivt@ya.ru
DATA ABOUT THE AUTHORS Preobrazhensky A.P., Choporov O.N.
Voronezh Institute of High Technologies 73A, Lenin Str., Voronezh region, 394043, Russian Federation
komkovvivt@ya.ru