На правах рукописи
АНЮТИН АЛЕКСАНДР ПАВЛОВИЧ
РАЗВИТИЕ АСИМПТОТИЧЕСКИХ И ЧИСЛЕННЫХ МЕТОДОВ МОДЕЛИРОВАНИЯ ДИФРАКЦИОННЫХ ПОЛЕЙ СИГНАЛОВ В СРЕДАХ С ДИСПЕРСИЕЙ
Специальность: 01.04.03 - радиофизика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора физико-математических наук
Москва - 2008
Работа выполнена на кафедре Информационных технологий и естественно научных дисциплин Российского Нового Университета, г. Москва
Официальные оппоненты:
доктор физико-математических наук, профессор, заслуженный деятель науки и техники РФ, лауреат государственной премии РФ, зав. кафедрой Физфака МГУ Куницын Вячеслав Евгеньевич доктор физико-математических наук, профессор, заслуженный деятель науки и техники РФ зав. кафедрой МГТУГА Козлов Анатолий Иванович доктор физико-математических наук, профессор нач. отдела ИЗМИР РАН Черкашин Юрий Николаевич Ведущая организация :Ц Московский государственный институт радиотехники, электроники и автоматики (технический университет МИРЭА)
Защита диссертации состоится "___Ф__________2008 г. в 16 часов на заседании диссертационного совета Д 212.156.06 при Московском физико-техническом институте по адресу: 117 393, г. Москва, ул. Профсоюзная, д. 84/32, корп. В-2.
Отзывы направлять по адресу:141 700, г. Долгопрудный, Московской области, Институтский переулок, д. 9, МФТИ, учёному секретарю
С диссертацией можно ознакомиться в научно-технической библиотеке МФТИ.
Автореферат разослан "___"______ 2008 г.
Ученый секретарь диссертационного совета Д 212.156.к.т.н, доцент Н.П. Чубинский
Актуальность темы. Представляемая диссертация посвящена разработке асимптотических и численных методов описания искажений негармонических и гармонических полей в линейных средах в случаях, когда необходимо учитывать искажения частотного или пространственного (углового) спектра сигнала. Интерес к подобным задачам в последнее время связан с интенсивным развитием техники экспериментальных исследований в различных областях прикладной радиофизики, развитием средств вычислительной техники и математики, а также потребностями и внутренней логикой развития теоретических исследований, поскольку во многих приложениях радиофизики уже нельзя считать волновые процессы узкополосными, а свойства среды и устройств, осуществляющих прием (обработку) сигналов не влияющими на искажения частотного или пространственного (углового) спектра сигнала.
Так, например, анализ и интерпретация данных, получаемых современными цифровыми ЛЧМ ионозондами или цифровыми ионозондами типа "Сойка-6000", "Авгур"; разработка их дальнейших модификаций; расширение информативности действующих ионозондов и станций вертикального (ВЗ) и наклонного (НЗ) зондирования ионосферы Земли; увеличение потенциала станций загоризонтной радиолокации; ряд задач антенной техники и волоконной оптики, связанных с рассеянием, фокусировкой и канализацией энергии радиоволн в естественных средах и искусственных средах (метаматериалах); диагностика плазменных образований различной природы приводят к необходимости разработки теоретического аппарата, адекватного возросшим требованиям и запросам практики и теории, связанным с учетом нестационарности, способов излучения и обработки сигналов, частотной дисперсии и неоднородности среды, а так же преобразованием углового (пространственного) спектра сигналов вследствии сложной геометрии рассевающего тела (поверхности) или особых свойств искусственно созданной среды - метаматериала. Отметим, что активное освоение техникой СВЧ все более высокочастотного диапазона волн, создание новых СВЧ устройств, использование широкополосных импульсов и искусственных сред (метаматериалов) для создания новых приборов также связано с необходимостью учета дисперсии среды, широкополосности сигналов и дифракционных эффектов, обусловленных как особенностями геометрии рассеивающих тел (поверхностей), так и типом среды компактных рассеивателей.
Как известно, первые теоретические исследования в области распространения нестационарных волн в диспергирующей среде (плазме) восходят к Зоммерфельду А. и Бриллюэну (1914 г.). Эти и последующие работы, наиболее значительные из которых выполнены Гинзбургом В.Л., Бреховских Л.М., Вайнштейном Л.А., Жекулиным Л.А., Гершманом Б.Н., Денисовым Н.Г., Блиохом П.В., использовали представление решения в виде интеграла Фурье-Лапласа. Анализ такого интегрального представления проводился, как правило, приближенно в предположении об узкополосности сигнала и охватывал ограниченный круг сред и сигналов.
Начиная с шестидесятых годов 20 века стал развиваться более общий подход к подобным задачам, основанный на асимптотическом решении соответствующих волновых нестационарных уравнений. Такой подход получил в дальнейшем название метода пространственно-временной геометрической оптики (ПВГО или ВГО), и развивался в работах Люиса Р., Фелзена Л., Вакмана Д.Е., Кравцова Ю.А., Островского Л.А., Степанова И.О., Бабича В.М, Булдырева В.С., Молоткова И.А. Метод ПВГО использовал понятие пространственно временного луча и давал неравномерную одночленную асимптотику поля сигнала в диспергирующей среде. Ограничения метода ПВГО связаны с тем, что он не описывал явлений пространственно-временной дифракции, проявляющиеся в искажении первоначальной формы сигнала, и которые присущи широкополосным сигналам конечной длительности.
В 1965 г, вышла работа Люиса Р., в которой была предпринята попытка перенести идеи геометрической теории дифракции (ГТД) гармонических сигналов (общеизвестные работы Келлера Д,, Кинбера Б.Е и Боровикова В.А.) на одномерный случай нестационарного уравнения Кляйна-Гордона и, тем самым, была сделана первая попытка описать явления временной дифракции сигналов с использованием лучевых представлений. В дальнейшем такой подход развивался нами, Фелсеном Л, Орловым Ю.И., Лукиным Д.С. и Крюковским А.С. и получил название временной геометрической теории дифракции (ВГТД).
Следует отметить альтернативный подход к построению асимптотик полей сигналов в неоднородных средах, основанный на методе Маслова В.П. Использование на практике этого метода, как и представление решения в виде интеграла Фурье-Лапласа, приводил к сложной проблеме численного интегрирования быстроосциллирующих функций. Однако в настоящее время, благодаря работам по специальным функциям катастроф (СВК), выполненным в научной группе под руководством Лукина Д.С. и Крюковского А.С., эту проблему в значительной мере можно считать решенной.
Не менее важной и актуальной с практической и теоретической точек зрения является задача определения параметров неоднородности диспергирующей среды (плазмы) по наблюдаемым искажениям первоначальных временных, частотных и пространственных характеристик сигналов, расчет характеристик широкополосных сигналов в случае многолучевости и каустик лучей ВГО и ВГТД, а также создание эффективных строгих методов решения задач дифракции в квазиоптической области с учетом сложной геометрии и среды компактных рассеивателей, приводящей к взаимодействию различных типов лучей и каустик.
Разработка эффективных асимптотических и численных методов решения упомянутых выше прямых и обратных задач представляется актуальной не только с точки зрения радиофизики, но и с общенаучной точки зрения, поскольку, вследствие универсальности математического аппарата, разработанные в диссертации подходы применимы к волнам и средам различной физической природы.
Цель работы состоит:
- в создании и разработке метода пространственно-временной геометрической теории дифракции (ВГТД), позволяющей в наглядных физических образах интерпретировать, предсказывать и в областях своей применимости рассчитывать дисперсионные искажения широкополосных сигналов;
- в разработке равномерных модификаций асимптотической теории искажений сигналов конечной длительности в диспергирующих средах, позволяющих рассчитывать особенности пространственно временной фокусировки широкополосных сигналов в условиях многолучевости, каустик лучей ВГО и ВГТД;
- в создании и разработке асимптотического подхода к решению обратной задачи - определение параметров модели неоднородности плазмы по искажению первоначальной временной формы огибающей сигнала ВЗ (ВГТД интерферометрия) или его спектра:
- в создании и разработке методов расчета полей, рассеянных ограниченными идеально проводящими и диэлектрическими телами сложной формы или компактными рассевателями из метаматериалов в квазиоптической области частот в случаях образования каустик (фокусов) лучей ГО и ГТД различного типа.
Научная новизна работы заключается в том, что впервые:
- разработана пространственно-временная геометрическая теория дифракции (ВГТД);
- проведено детальное исследование пространственно-временных лучевых картин лучей ВГО и ВГТД, позволявшее наглядно объяснять и предсказать особенности временных искажений широкополосных сигналов конечной длительности связанные:
а/ с затягиванием и асимметрией переднего и заднего фронтов сигнала;
б/ с существованием и развитием осцилляций на фронтах сигнала и его вершине;
в/ "лепестковым" характером изменения амплитуды сигнала для случая, когда его несущая частота близка к критической частоте плазмы, а также уменьшение числа и амплитуды лепестков при обработке сигнала преемником с ограниченной частотной полосой пропускания (сначала эти эффекты были нами предсказаны, а затем экспериментально обнаружены нами совместно с Кольцовым В.В на цифровом ионозонде);
г/ с отличием в эффективной длительности обыкновенного и необыкновенного сигналов ВЗ при его отражении от окрестности максимума F-области ионосферы;
д / со степенным характером убывания поля сигнала за каустикой ВГО лучей;
е /с уменьшением уровня осцилляций огибающей сигнала при наличии амплитудной модуляции;
ж / с особенностью фокусировки широкополосных сигналов на каустиках лучей ВГО;
- указан способ построения специальных функций, описывающих дисперсионные искажения сигналов в областях неприменимости ВГТД;
- получены новые равномерные по положению точки наблюдения асимптотические формулы для поля сигнала в одномерном и многомерном случаях, позволяющие оценить границы применимости метода ВГТД и определить поле сигнала:
а / вблизи мирового конуса будущего (случай трехмерной плазмы);
б / вблизи фронта сигнала (случай трехмерной неоднородной диспергирующей среды);
в / вблизи "оборванных" сложных каустик лучей ВГО и ВГТД;
г / при обработке сигнала фильтром с конечной частотной полосой пропускания;
д / при идеальной фокусировке (когда все лучи ВГО сходятся в одну пространственно-временную точку) сигнала конечной длительности;
- разработаны эффективные алгоритмы расчета характеристик полей сигналов ВЗ и НЗ в пределах одного скачка. учитывающие неоднородность; потери анизотропию ионосферы; кривизну фронта волны и диаграмму направленности (ДН) излучающей антенны; конечную длительность, закон модуляции излученного сигнала и способ его обработки приемным устройством;
- предложен и разработан способ асимптотического решения обратной задача - определение параметров модели неоднородности плазмы по искажению первоначальной временной формы огибающей сигнала ВЗ (ВГТД интерферометрия), частотной (ЧХ) характеристике сигнала ВЗ, текущего частотного спектра сигнала НЗ или пространственной структуре поля сигнала НЗ в окрестности границы мертвой зоны;
- предложен метод продолженных граничных условий - универсальный метод решения задач рассеяния гармонических волн ограниченными телами (экранами) сложной геометрии, сводящий исходную задачу дифракции к решению интегральных уравнений (ИУ) Фредгольма первого рода с гладким ядром и на основе использования вейвлетных базисов разработан эффективный алгоритм решения задач рассеяния полей идеально проводящими и диэлектрическими телами (экранами) сложной геометрии и больших электрических размеров (совместно с Кюркчаном А.Г.) ;
- в строгой постановке получено решение и исследована задача пространственной фокусировки волн плоской линзой Веселаго конечных размеров с учетом потерь и значений показателя преломления метематериала и задача рассеяния гауссова пучка призмой из метаматериала. Установлено, что особенности лучевых картин и структура каустик для такой линзы определяются значениями модуля показателя преломления. Показана принципиальная ограниченность разрешающей способности таких линз.
Достоверность положений, результатов, выводов является следствием: использования моделей, учитывающих основные процессы, происходящие при излучении, приеме и распространении сигналов в диспергирующей среде, а так же рассеянием волн на ограниченных телах, использовании адекватного математического аппарата - уравнений Максвелла, согласия теоретических результатов с результатами измерений.
Научная и практическая значимость работы заключается в разработке асимптотических методов теории взаимодействия сигналов с неоднородной диспергирующей средой и развитии численных методов решения задач дифракции ограниченными телами в строгой постановке, которые позволили:
- уточнить ряд результатов традиционной теории гармонических полей и широкополосных сигналов в условиях образования сложных каустик, геометрии и среды компактного рассеивателя;
- обобщить результаты традиционной асимптотической теории гармонических полей на негармонические поля;
- поставить и решить ряд новых важных для практики прямых и обратных задач;
- создать и развить методы расчета полей, рассеянных идеально проводящими и диэлектрическими телами (экранами) сложной геометрии и больших электрических размеров в строгой постановке и контролируемой точностью.
К числу последних относятся: анализ и интерпретация временной структуры сигналов вертикального и наклонного зондирования ионосферы с учетом способа излучения и обработки в приемном устройстве; определение параметров модели неоднородности ионосферной плазмы по искажениям пространственной структуре КВ сигнала вблизи мертвой зоны;
ВГТД интерферометрия; расчет диаграмм рассеяния полей вогнуто-выпуклых идеально проводящих цилиндров и экранов в условиях образования волн шепчущей галереи, их фокусировки на каустиках и учетом взаимодействия с краями экрана; исследование структуры поля в областях фокусировки линзы Веселаго конечных размеров с учетом потерь и значений показателя преломления метаматериала линзы.
Полученные в работе результаты могут найти приложение в других разделах радиофизики, физики, а так же таких, как антенная техника, физика плазмы, акустика, физика космоса, теория упругости, гидромеханика.
Реализация результатов работы.
Полученные в диссертации результаты, алгоритмы и разработанные численные методы были использованы при проведении научно-исследовательских работ: "Скала", "Интерференция", выполненных по постановлению директивных органов, а так же в рамках грантов РФФИ : N 00-02-17639, N 02-02-06129, N 03-02-16336, N 06-02-16804, N 06-02-16483.
Результаты численного моделирования искажений сигналов были использованы в ИРЭ АН УССР, ИЗМИРАН СССР, НИИДАР, при исследованиях поля гармонического и негармонического сигнала наклонного зондирования вблизи границы мертвой зоны и вертикального зондирования ионосферы.
Результаты диссертации использовались в прочитанных автором курсах лекций:
"Асимптотические методы технической электродинамики", "Геометрическая оптика неоднородных сред и ее обобщенияФ. Получено положительное решение на заявку на изобретение:
"Способ определения параметров максимума электронной концентрации ионосферы".
Апробация -работы. Результаты диссертации докладывались на VП, VШ и IX Всесоюзных симпозиумах по дифракции и распространению волн (Ростов-на-Дону, 1977г.;
ьвов, 1981 г.; Телави, 1985 г. на ХII, ХIII и Х1V Всесоюзных конференциях по распространению волн (Горький, 1981 г.; Ленинград, 1984 г.; Алма-Ата, 1987 г.); на Ш Всесоюзной конференции по взаимодействию электромагнитных излучений с плазмой (Алма-Ата, 1982 г.);
Всесоюзной конференции "Машинное проектирование устройств и систем СВЧ (Тбилиси, 1979г.); Всесоюзной конференции НТОРЭС им. А.С Попова (Москва, 1976 г.); X и XII Всесоюзной школе-семинаре по дифракции и распространению волн (Москва-Волоколамск, 1984 г., Москва 2001 г.); Всесоюзных совещаниях по проблеме "Распространение декаметровых радиоволн в ионосфере" (Звенигород, 1978 г.; Троицк, 1988 г.); Всесоюзном семинаре "Исследование неоднородностей ионосферы радиофизическими методами (Душанбе, 19г.); XIV International Conference on Applied Electromagnetics and Communication, Croatia, Dubrovnik, 2001 y.; VI International Conference on Electromagnetic and Light Scattering by Nonspherical Particles, University of Florida, Gainesville, USA 2002 y.; Seville, Spain, 2005 y.; XX Всероссийской Конференции по Распространению Радиоволн, Н. Новгород, 2002 г. ; XXVII General Assembly of URSI, Maastricht, 2002 y.; WSEAS Int. Conf. on Microwaves, Antennas and Radar Systems, Spain, 2002 y.; WSEAS Int. Conference on Applied Mathematics and Computer Science (AMCOS 2002, AMCOS 2004), Rio de Janeiro, Brazil, 2002 y., 2004 y.; International Symposium on Antennas, France, Nice, 2002 y.; 2002 Asia-Pacific Microwave Conference (APMC 2002), Japan, Kyoto, 2002 y., 2003 y.; MATA 2003 Conference, Cancun, Mexico, 2003 y.; IV ISAAC Congress, Toronto, Canada, 2003; International Conference on Electromagnetics in Advanced Applications (ICEAA03), Torino, Italy, 2003 y.; XIV, XV and XVI Simposio Internacional de Matematicos Aplicados a las Ciencias, San Jose, Costa Rica, 2004 y., 2006 y., 2008 y.; International Symposium on Electromagnetic Theory, Pisa, Italy, 2004 (2004 URSI EMTS); International Symposium on Antennas and Propagation, Sendai, Japan, 2004 y. (2004 URSI ISAPТ04); 4th International Symposium on Information and Communication Technologies (ISICT05), Cape Town International Convention Center, Cape Town, South Africa, 2005 ; XXI Всероссийская научная конференция по Распространению радиоволн, Йошкар-Ола, 2005.
Полученные в диссертация результаты неоднократно докладывались и обсуждались на семинарах в РосНОУ, ИРЭ АН СССР, РА АН УССР, ИЗМИРАН СССР и Общероссийском семинаре УМатематическое моделирование волновых процессовФ.
В 1981, 1984 и 1985 гг. материалы диссертации включались в список важнейших достижений по секции 5а Научного Совета по комплексной проблеме "Распространение радиоволн" АН СССР.
Публикации. Основные научные результаты диссертации отражены в 24 статьях ( 3 - в Докладах Академии Наук, 6 - в Изв. ВУЗов, сер. Радиофизика, 10 - в Радиотехника и Электроника, 3 - в Электромагнитные волны и электронные системы, 2 - в Геомагнетизм и Аэрономия ), из них 15 работ в соавторстве. Кроме того, 7 статьей опубликованы в иностранных журналах : 4 - JQSRT, 2 - Revista de Matematica: Teoria y Aplicaciones, 1 - WSEAS Transactions on Systems. Остальные 35 работы опубликованы в тематических сборниках Академии Наук СССР, изданиях международных и отечественных симпозиумов и конференций.
В работах, выполненных в соавторстве, автору настоящей диссертации принадлежит равноценное участие в постановка задач и разработке методов их решения, получении и обсуждении результатов.
I. Основные положения диссертации.
1.1. Построена и развита асимптотическая теория распространения негармонических полей в плавнонеоднородной диспергирующей среде (пространственно-временная геометрическая теория дифракции (ВГТД)), позволяющая эффективно рассчитывать сложные явления пространственно-временной дифракции сигналов, а также их интерпретировать в наглядных физических образах. Исследованы особенности лучевых семейств и каустик, образованных лучами ВГО или ВГТД. Благодаря этому, теория дает возможность адекватно ставить и исследовать процессы распространения и приема сложных широкополосных сигналов в присутствии неоднородной диспергирующей среды (плазмы).
1.2. Предложен метод ВГТД-интерферометрии - способ определения параметров модели неоднородности плазменной среды по характерным дисперсионным искажениям первоначальной формы сигнала или его частотного спектра. Это позволяет решать задачи прогноза модели профиля неоднородности плазмы и ее потерь по наблюдаемым экспериментально характеристикам отраженных от плазмы сигналов.
1.3. Предложен метод продолженных граничных условий - универсальный метод решения задач рассеяния гармонических волн ограниченными телами (экранами) сложной геометрии, сводящий исходную задачу дифракции к решению интегральных уравнений (ИУ) Фредгольма первого рода с гладким ядром (совместно с Кюркчаном А.Г.). Впервые для решения таких ИУ были использованы вейвлетные базисы (Хаара и Батла-Лемарье) и на их основе разработан эффективный алгоритм решения таких ИУ в квазиоптической области частот, позволивший впервые в строгой постановке численно решить задачу рассеяния полей вогнуто-выпуклых идеально проводящим цилиндром и экраном в условиях образования волн шепчущей галереи, их фокусировкой и взаимодействием с краями экрана;
1.4. Разработан модифицированный метод дискретных источников (совместно с Кюркчаном А.Г.) и на его основе впервые в строгой постановке получено решение задачи дифракции цилиндрических волн плоской линзой Веселаго конечных размеров с учетом потерь и значений показателя преломления метематериала и исследованы особенности фокусировки таких линз. Установлено, что особенности лучевых картин в линзе Веселаго связаны с образованием каустики с одной точкой возврата и двумя точками прекращения, пространственное расположение которой зависит от значений модуля показателя преломления. Показана принципиальная ограниченность разрешающей способности таких линз.
Структура и объем диссертации. Диссертация состоит из введения. пяти разделов, заключения и списка литературы из 224 наименований. Объем работы - 369 страниц. Основная часть содержит 342 страниц машинописного текста, 155 рисунка и 2 таблиц.
СОДЕРЖАНИЕ РАБОТЫ
.
Во введении обосновывается необходимость развития асимптотической теории анализа и синтеза искажений сигналов в диспергирующих средах; сформулирована цель исследования; обоснована его актуальность, научная новизна, практическая ценность, реализация результатов и основные научные результаты; дан краткий обзор результатов исследований других авторов; приведены основные положения, выносимые на защиту.
Первый раздел диссертации посвящен обобщению и развития идей ГТД на нестационарные процессы в диспергирующей среде (плазме).
В п. 1.1 рассматривается асимптотика задачи Коши с осциллирующими начальными данными для нестационарного волнового уравнения:
p (r ) 1 2U U -- 2 U =, (1) c2 t2 c описывающего процесс распространения радиосигнала U (r,t)в плазме. Асимптотическое решение уравнения (1) при >>1 ищется в виде разложения:
-m U (r,t) = Umgo exp{igo}+ Umkr exp{ikr}]. (2) (i) [ m=Umgo Umkr Получены рекуррентные уравнения для амплитуд, и уравнения эйконала для go kr фазовых функций,. Для нулевого приближения ( m = 0) установлена связь начальных данных с начальными условиями для уравнений переноса и эйконала и построены их решения. Показано, что в главном ( m = 0) члене асимптотики (2) первое слагаемое соответствует приближению временной геометрической оптики (ВГО), а второе слагаемое может быть интерпретировано как вклад пространственно временных краевых лучей ВГТД. В заключение пункта сформулированы четыре универсальных положения асимптотической временной геометрической теории дифракции (ВГТД), основанной на двухчленном представлении главного члена асимптотики поля сигнала в произвольной диспергирующей среде.
В п. 1.2 путем асимптотического вычисления вкладов критических точек различного типа (стационарных точек, полюсов, точек ветвления, точек разрыва интегранды) однократного интеграла Фурье:
U (z,t) = S0(z,)T () exp[it - iz () / c]d, (3) 2 - описывающего распространение радиосигнала в диспергирующей среде, характеризуемой относительной действительной диэлектрической проницаемостью (), установлено, что вкладу в асимптотику поля сигнала каждого типа критической точки соответствует вклад соответствующего ВГО-луча, краевого ВГТД-луча или вершинного ВГТД-луча. В (3) S0() - первоначальный спектр радиосигнала, T () - частотный коэффициент передачи приемного устройства. Получены выражения для дифракционных коэффициентов краевых и вершинных лучей ВГТД. Показано, что при слиянии критических точек разного типа возникают интегральные представления специальных функций, пригодных для описания поля сигнала в областях неприменимости ВГТД.
В п. 1.3 получено приближение ВГТД для поля видеосигнала в однородной плазме без потерь. Введено понятие о "ближней", "промежуточной" и "дальней" зонах видеосигнала.
Получены выражения полей краевых ВГТД-лучей высших порядков и найдены выражения для дифракционных коэффициентов таких лучей.
Второй раздел диссертации посвящен построению равномерных асимптотик поля широкополосного сигнала в среде с дисперсией, справедливые в областях неприменимости одномерной ВГТД - в окрестностях оборванных каустик ВГО-лучей, каустик ВГТД-лучей и областей полутени. Равномерные асимптотики находятся путем асимптотического вычисления соответствующих интегральных представлений поля сигнала.
В п. 2.1 получена равномерная асимптотика поля широкополосного радиосигнала в произвольной диспергирующей среде (3), когда точка наблюдения оказывается расположенной вблизи границы свет-тень лучей ВГО или краевых ВГТД-лучей (фронтов сигнала).
Этот случай возникает тогда, когда в интегральном представлении поля сигнала точка стационарной фазы расположена вблизи полюса интегранды или в окрестности края области интегрирования. Показано, что в обоих случаях поле сигнала выражается через интеграл Френеля F(- ). Получены формулы, связывающие коэффициенты равномерной асимптотики с амплитудами и фазами лучей ВГО и ВГТД. Кроме того, показана регулярность найденных коэффициентов равномерной асимптотики.
В п. 2.2 построена равномерная асимптотика поля радиосигнала в диспергирующей среде, принимаемого фильтром с полосовой частотной характеристикой. При этом учитывается, что несущая частота сигнала может быть произвольно расположена относительно локальной частоты принимаемого сигнала и Участоты срезаФ (граничной частоты) частотной характеристики приемника. Показано, что равномерная асимптотика поля сигнала в этом случае выражается через обобщенный интеграл Френеля G( ,i ), интегральную 2 показательную функцию E1(i[ - ]) и интеграл Френеля F(- ). Установлена связь коэффициентов такой асимптотики с амплитудами и фазами лучей ВГО и ВГТД.
Показана регулярность коэффициентов найденной асимптотики.
В п. 2.3 получена равномерная асимптотика поля сигнала в одномерной диспергирующей среде для случая, имеет место простая временная фокусировка (временное сжатие) сигнала конечной длительности. В этом случае интегральное представление имеет две стационарные точки, произвольно (в том числе и близко) расположенные относительно друг друга и края области интегрирования. Физически, такая ситуация соответствует случаю расположения точки наблюдения вблизи границы свет-тень ВГО-лучей и простой каустики ВГО-лучей (каустика ВГО-лучей имеет точку прекращения). Показано, что равномерная асимптотика поля сигнала при этом выражается через неполную функцию Эйри:
v(,) = exp[i(u /3 - u)du и ее частные производные. Получены формулы, связывающие коэффициенты равномерной асимптотики с амплитудами и фазами двух ВГО-лучей и одного ВГТД-луча. Показана регулярность найденных коэффициентов в формулах равномерной асимптотики.
В п. 2.4 построена равномерная асимптотика для поля радиосигнала в одномерной диспергирующей среде, у которой частотная зависимость групповой скорости имеет изолированную точку экстремума. В этом случае поле сигнала выражается частотным интегралом Фурье, у которого две стационарные точки произвольно расположены относительно полюса интегранды. Физически, при этом точка наблюдения расположена в окрестности как границы свет-тень ВГО-лучей так и каустики краевых лучей ВГТД. Показано, что формула равномерной асимптотики в этом случае выражается через функцию Эйри ее производную и интеграл Эйри-Френеля:
exp{i( / 3-} V (,) = d - -.
Установлена связь коэффициентов такой асимптотики с амплитудами и фазами ВГО-луча и двух краевых ВГТД-лучей. Показана регулярность найденных коэффициентов в формулах равномерной асимптотики.
В п. 2.5 предложены асимптотические формулы равномерных интегральных представлений (ИП) для поля сигнала в одномерной диспергирующей среде. Такие интегральные представления физически представляют собой суперпозицию "временных" или "частотных" мод и получаются из строгого двукратного интегрального Фурье представления решения задачи путем асимптотического вычисления одного из интегралов либо :методом стационарной фазы, либо одной из модификаций метода стационарной фазы (например, см. п.2.1-2.4).
Удобство таких ИП связано с тем, что в ряде практически важных случаев при вычислении структуры сигналов не требуют предварительного анализа картин лучей ВГО и ВГТД. На основе аппроксимации амплитудной функции интегранды с помощью вэйвлетов Умексиканская шляпа или гауссоида Ф получены соотношения одномерного метода временных гауссовых пучков, удобного для численных расчетов сигналов с плавным законом изменения первоначальной амплитуды сигнала (или его частотного спектра).
В третьем разделе диссертации на основе метода эталонных функций получены равномерные асимптотики поля сигнала, распространяющегося в произвольной трехмерной плавнонеоднородной диспергирующей среде (скалярный случай).
В п. 3.1 построена равномерная асимптотика уравнения (I) для поля сигнала в неоднородной плазме для точек наблюдения расположенных вблизи мирового конуса будущего (предвестника сигнала). Асимптотика поля сигнала в этом случае выражается через функцию Бесселя Jm( ) и ее производную :
. (4) U(r,t) = 2ei{AJm()+ BJm()} Для коэффициентов и аргументов асимптотики (4) получена система дифференциальных уравнений в частных производных. Показано, что величины A, B,, в (4) выражается через амплитуды и фазы дух краевых лучей ВГТД, а индекс функций Бесселя "m" определяется характером изменения огибающей сигнала в момент его включения. Кроме того, для описания поля сигнала в средах со сложной дисперсионной характеристикой предложена новая специальная функция волновых катастроф - многомерная функция Эйри-Бесселя F(bn):
dy = 1 b-1 yn+1 F(bn) = exp -i + Pn-1( y) + 2 yn +1 - (5) -=-1[ f1( y) f2(y)] = [ f1( y)]-1[ f2( y)]d ;
где 1 1 b0 -1 f1(y) = exp [ ] -i - b1ydy = -i2 J0(2 | b0 | b1 )H (-b0) 2 yy - , dy = V (bq ) 1 yn+1 -1 f2(y) = exp [ ] -i -b y + Pn-2(y) + 2 n +1 - bn - коэффициенты полинома Pn-1(y); Ф-1[ f ] - обратное преобразование Фурье, а знак , как обычно, означает свертку от произведения функции Бесселя J0(2 | b0 | b1 )H (-b0) ( H (-b0) - функция единичного скачка) на многомерную функцию Эйри V (bq ).
В п. 3.2 получена равномерная асимптотика скалярного четырехмерного волнового уравнения для сигнала в произвольной диспергирующей среде, справедливая в случае образования сложной каустики ВГО лучей с одной точкой обрыва (возникновения). Решение имеет следующий вид:
M = U (r,t) = -1/ M +1{AI(l,) + i -S / M +1BS I(S,) S S =(6) -i-M / M +1C I(l,)}exp[i ] где I(l,) неполная многомерная функция Эйри:
l M - M +1 I(l,)= expi + l d.
M +1 l l= Показано, что коэффициенты А, A, BS,C и аргументы ,l, в асимптотическом выражении (5) выражается через амплитуды и фазы основных и одного краевого луча, проходящих в точку наблюдения. С использованием принципа асимптотической эквивалентности указана связь коэффициентов и аргументов равномерной асимптотики (6) с амплитудами и фазами лучей ВГО и краевого луча ВГТД.
В п. 3.3 построена равномерная асимптотика четырехмерного волнового уравнения для поля сигнала в произвольной среде с дисперсией, справедливая в случае, когда сложная каустика краевых лучей ВГТД располагается в окрестности границы свет-тень лучей ВГО.
Асимптотика поля сигнала имеет вид:
M = U(r,t) = -1/(M +1){AV (l,) +i -S /(M +1)BS V (S,) + S S=(7) +i-M /(M +1)C V (l,)}exp[i];
где многомерный интеграл Эйри-Френеля V (l,) :
M +1 M -1 l d V (l,) = expi + l M +1 l - l=- и ее частные производные. С помощью техники метода эталонных функций получены уравнения для неизвестных коэффициентов и аргументов асимптотики (7). На основе принципа асимптотической эквивалентности проиллюстрирована связь коэффициентов и аргументов асимптотики с амплитудами и фазами краевых лучей ВГТД и луча ВГО.
В п. 3.4 получена равномерная асимптотика скалярного четырехмерного волнового уравнения для поля сигнала в произвольной диспергирующей среде справедливая в случае, когда границы свет-тень лучей ВГО и ВГТД близки друг к другу. В качестве эталонных функций здесь используется обобщенный интеграл Френеля и интеграл Френеля. Решение имеет вид:
U(r,t)=Aexp{iВГО}{1-G(- , )}+ BC, (8) + exp{iкр}F( ) + exp{iв} G(- , ) A, B,C , F(- ) где, - неизвестные функции, - интеграл Френеля, - обобщенный интеграл Френеля. Неизвестные коэффициенты и аргументы асимптотики (8) определены с помощью метода асимптотической эквивалентности и выражены через амплитуды и фазы лучей ВГО и ВГТД (краевых и вершинных). Показано, что границы применимости ВГТД определяются двумя (а не одним) неравенствами, являющимися следствием перехода (8) в неравномерную асимптотику - ВГТД приближение.
В четвертом разделе приводятся и обсуждаются результаты численного и асимптотического исследования одномерных задач взаимодействия радиосигналов и видеосигналов со средой с плазменным законом дисперсии.
В п. 4.1 приводятся и анализируются результаты численных расчетов искажений радиосигналов с амплитудной и частотной модуляцией, распространяющихся в однородной плазме с потерями, полученные на основе численного интегрирования двукратного интеграла Фурье. Кроме того, показано, что неравномерная асимптотика двукратного интеграла Фурье для поля сигнала в такой среде совпадает с приближением ВГТД. Приведены результаты расчетов временной структуры огибающей и фазы АМ-ЧМ радиосигналов для различных значений пройденной дистанции z, времени нарастания фронтов сигнала TФ и его длительности T, скорости изменения девиации локальной частоты ЛЧМ сигнала в окрестности плазменной частоты среды и ширины частотной полосы приемника. Показано, что уменьшение длительности фронтов сигнала, как и общей длительности сигнала T, приводит к уменьшению значения параметра zгр, определяющей пространственное положение границы применимости метода ВГО.
В п. 4.2 обсуждаются полученные результаты расчетов искажений различных видеосигналов, распространяющихся в однородной плазме с потерями. Алгоритмы расчета поля видеосигнала основаны на использовании метода сеток, вычислении асимптотики соответствующих интегралов Фурье и представлении поля сигнала в виде ряда Неймана по функциям Бесселя. Показано, что характер изменения поля видеосигнала в момент его включения (выключения) определяет поведение (осцилляции) "хвоста" сигнала. Введено понятие ближней ( = Pz / c <1), промежуточной ( 1) и дальней ( >>1) зонах видеосигнала. Кроме того, обнаружено, что если точка наблюдения находится в ближней зоне видеосигнала, то потери среды уменьшают влияние плазмы на искажения первоначальной формы видеосигнала, а в дальней зоне они приводят к уменьшению, как числа осцилляции сигнала, так и их максимальных значений. Установлено, что при условии > 3 приближение ВГТД с точностью не хуже 5% описывает искажения видеосигнала с первоначальной формой в виде единичной ступеньки и отсутствии потерь (см. Рис.1). Пример расчета, иллюстрирующий влияние разрывов производной первоначальной формы видеосигнала (краевых ВГТД-лучей первого порядка) представлен на Рис.Рис.1 Искажения единичной ступеньки Рис.2 Искажения трапециидального ( = 0, = pz / c ); сигнала: ( = 0, = pz / c =1);
Строгое решение - непрерывная кривая (1); Приближение ВГО - пунктир (3);
Приближение ВГТД - точки (2).
Кроме того, показано, что учет внешнего постоянного магнитного поля приводит к развитию сложных разномасштабных осцилляции обыкновенной и необыкновенной компонент видеосигнала.
В п. 4.3 построено ВГТД приближение для поля плоского радиосигнала, отраженного от неоднородной плоскослоистой плазмы. Полученные выражения для амплитуд и фаз лучей ВГО и ВГТД конкретизированы на монотонных эталонных и немонотонных законах изменения квадрата плазменной частоты от одной пространственной координаты. Установлено существование временного аналога пространственного луча Педерсона.
В п. 4.4 приводятся результаты расчета и интерпретируются искажения плоских АМ-ЧМ радиосигналов, отраженных от плазменных слоев, квадрат плазменной частоты P(z) которых изменяется по линейному и экспоненциальному законам. Тенденции в изменении характера искажений огибающей фазы отраженного АМ-ЧМ сигнала, предсказанные на основе ВГТД анализа, проиллюстрированы примерами строгих численных расчетов. Алгоритм расчета структуры отраженных сигналов основан на использовании метода быстрого преобразования Фурье и созданных программ расчета функции Эйри и гамма функции произвольного комплексного аргумента. Показано, что одночленная локальная асимптотика поля сигнала может приводить к значительным погрешностям при расчете формы отраженного сигнала, а учет малых потерь среды вызывает дополнительное несимметричное уменьшение амплитуд лучей ВГО и ВГТД без нарушения их взаимных фазовых соотношений.
В п.4.5 анализируются результаты расчета и интерпретируются искажения плоских АМЧМ радиосигналов, отраженных от плазменного слоя Эпштейна и полу бесконечного параболического плазменного слоя. Строгие решения таких задач использованы для нахождения асимптотики поля отраженного сигнала и для строгих численных расчетов. На основе разработанных алгоритмов подробно исследованы искажения радиосигналов с первоначальной трапециидальной огибающей, ЛЧМ и AM радиосигналы, несущая частота 0 которых близка к критической частоте m плазменного слоя. Пример рассчитанной структуры огибающей ЛЧМ сигнала, излученного нитью магнитного тока и отраженного от параболического плазменного слоя при его фокусировке на простой временной каустике ВГО-лучей представлен на Рис. 3. Параметры слоя и сигнала характеризовались следующими значениями: высота слоя z0 = 200км, полутолщина слоя h =100км, критическая частота слоя fm = 7МГц, средняя несуm щая частота ЛЧМ сигнала f0 =10МГц, длительность ЛЧМ сигнала T = 100 мкс, девиация частоты ЛЧМ сигнала =16%, угол падения =58[град].
|U/Um| t[сек] Рис.3 Нормированная огибающая ЛЧМ (=-16%;=58[град]) Примеры расчета огибающей радиосигнала с первоначальной формой огибающей в виде трапеции длительности Т =100 мкс и Т = 50МKс, временем нарастания фронта сигнала Тф =10МKс, несущей частотой f0 = 7.03 МГц которая близка к критической частоте fm = 7.04 Мгц плазменного слоя и принятого приемником, частотна характеристика K( f ) которого изображена в верхней части рисунка, иллюстрирует рис.4. Из рисунка следует, что конечность частотной полосы пропускания = 30КГц приемника ограничивает число "лепестков" огибающей сигнала.
Рис.4 Огибающая радиосигнала на выходе приемника с неограниченной и ограниченной полосой пропускания 2 (T = 70Мкс., /0 =10-4, кривая 1 - 2 =, кривая 2 - 2= 30КГц ).
При малых длительностях TФ фронтов радиосигнала (T /TФ 0.2, где T - общая длительность сигнала) и несущей частоте близкой к критической частоте плазмы отраженный сигнал имеет "лепестковую" структуру. При этом, конечность частотной полосы пропускания приемника приводит к ограничению числа "лепестков" - эффекту временного окна. Сопоставление результатов проведенных теоретических расчетов с экспериментальными данными, полученными ВВ Кольцовым на цифровом ионозонде УСойка 6000Ф в эксперименте по вертикальному радиозондированию невозмущенной F-области ионосферы показывает их хорошее совпадение друг с другом. Кроме того, теоретически и экспериментально показано, что учет внешнего постоянного магнитного поля (в продольном приближении) приводит к уменьшению в 2,3 раза максимального расплывания необыкновенной компоненты сигнала по сравнению с обыкновенной компонентой. Исследована зависимость искажений сигналов с амплитудной модуляцией от типа неоднородности среды и параметров амплитудной модуляции.
В п.4.6 приводятся результаты исследования искажений широкополосных видеосигналов, отраженных от слоистой среды с потерями. Получены строгие аналитические формулы, позволяющие интерпретировать временную структуру отраженных сигналов. Разработан алгоритм расчета искажений временной структуры сигналов отраженных от трехслойной диспергирующей среды с потерями, моделирующей ситуацию зондирования поверхности Земли или моря, покрытых слоем снега или льда. Исследованы особенности искажений отраженных сигналов, в том числе связанные с наличием угла Брюстера. Так на Рис. 5а,5б представлены формы U (t) отраженных видеосигналов длительности 100 наносекунд, двух поляризаций (1- Е поляризованная плоская волна, 2- Н поляризованная плоская волна) при двух углах падения на трехслойную систему воздух - лед - сырая Земля.
U(t) U(t) t t а ) угол падения 30 [град] б) угол падения равен углу Брюстера Рис.5 Временная структура отраженного видеосигнала В п.4.7 предложен и реализован метод ВГТД интерферометрии, позволяющий по характерным дисперсионным искажениям огибающей радиосигналов определять параметры модели диспергирующей среды. Приводится и обсуждается примеры использования этого метода к задаче определения параметров параболической модели невозмущенной F области ионосферы для радиосигналов с частотой близкой к критической частоте такой области и задаче зондирования среды ЛЧМ сигналом при образовании каустики отраженных ВГОлучей.
В пятом разделе приводятся результаты расчетов гармонических и негармонических полей сигналов для точек наблюдения, расположенных вблизи каустик и фокуса ГО лучей (многомерные задачи). В основе разработанных алгоритмов лежит метод быстрого преобразования Фурье, созданная программа вычислений интеграла от быстро осциллирующей функции, и решение интегральных уравнений Фредгольма первого рода с гладким ядром с использованием вэйвлетных базисов.
В п. 5.1 обсуждаются результаты расчетов структуры поля КВ гармонического сигнала вблизи невозмущенной границы зоны молчания, полученные на основе строгих формул, приближения геометрической оптики и локальной одночленной каустической асимптотики.
В качестве модели F области ионосферы использовался полу бесконечный параболический плазменный слой. Сравнение результатов расчетов с данными эксперимента показало высокую точность такой аппроксимации, необходимость учета потерь среды и формы диаграммы направленности излучающей антенны. Пример сравнения результатов проведенных расчетов по полученным формулам (без учета потерь среды - верхняя кривая рисунка, с учетом потерь среды и данные эксперимента - нижняя кривая) с результатами экспериментального измерения поля вблизи невозмущенной границы зоны молчания (полученные научной группой Ю.М. Ямпольского) представлен на рис. 6.
Рис.6 Структура поля вблизи границы мертвой зоны:
кривая 1 - расчет при = 0, кривая 2 - расчет при /0 =10-4.
В п. 5.2 приводится численное и асимптотическое решение негармонической задачи п. 5.1. При этом полагается, что точка наблюдения располагается в окрестности границы зоны молчания, образованной на несущей частоте сигнала. В результате асимптотического анализа строгого решения получены условия, при которых исходная многомерная задача может быть сведена к решению более простых - двухмерных по пространственным координатам или по пространственной и временной координате задач. Показано, что структура текущего спектра сигнала носит осциллирующий характер, зависящий от длительности сигнала и близости точки наблюдения к каустике, а его временная форма и дисперсионное расплывание существенно зависит от положения точки наблюдения относительно границы мертвой зоны.
В п. 5.3 приводятся и обсуждаются результаты экспериментов и теоретического исследования поля KB гармонического сигнала и радиосигнала в условиях возмущения мертвой зоны. В качестве модели возмущенной F-области ионосферы использовался плазменный слой с двумя параболическими законами изменения относительной диэлектрической проницаемости (z) :
при z < z при (z, f ) = (z, f ) z0 < z zгр (z, f ) + 2(z, f ) при z > zгр где 2 2 2 fm1 fm1 z - zm1 fm 2 fm 2 z - zm 1(z, f ) = 1 - + ( )2 2 (z, f ) = + ( ),.
2 2 2 f f ym1 f f ym Показано, что небольшое (~5%) возмущение основного параболического слоя приводит к качественному изменению структуры мертвой зоны - появлению дополнительных двух точек выхода каустики на поверхность Земли. На основе сравнения полученных результатов расчетов и данных экспериментов показано, что такая аппроксимация возмущений области ионосферы позволяет интерпретировать наблюдаемые в эксперименте типы искажений полей KB гармонических сигналов вблизи зоны молчания. Учет потерь среды, характерных для F области ионосферы, не вызывает нарушения взаимных фазовых соотношений полей интерферирующих лучей, а сказывается лишь на их амплитудах. Установлено, что для точек наблюдения, расположенных вблизи зоны молчания, уменьшение длительности радиосигнала не приводит к увеличению разрешения "мод" (лучей) сигнала. Примеры расчета структуры гармонического поля и радиосигнала с гауссовой огибающей в окрестности возмущенной границы зоны молчания с учетом потерь среды (кривая 2 Рис. 7) и без учета потерь среды (кривая 1 Рис. 7) представлен на рис.7 и рис. 8 соответственно. При расчетах полагалось, что параметры слоя и сигнала были равны: f = 9.996 МГц, fm1 = 7.2 МГц, z0 =140 км, ym1 =1км, zгр =170 км, zm1 = zm2 = 240 км для рис.6; для рис. 7 - полуширина сигнала = 50 МКс, fm2 = 7.3 МГц, zгр =195 км, x = 638.5 км - непрерывная кривая, x = 636.6 км - пунктирная кривая и x = 634.25 км- штрихпунктирная кривая (остальные данные совпадают с данными рис.6) и на оси времени все три сигнала приведены к времени прихода главного максимума огибающей сигнала tm.
Рис.7 Структура поля вблизи возмущенной границы мертвой зоны ( fm2 = 7.68МГц ).
Кривая 1 - = 0, кривая 2 - / f0 =10-4.
Рис.8 Структура огибающей радиосигнала вблизи возмущенной границы мертвой зоны непрерывная кривая - x = 638.5 км, пунктирная кривая - x = 636.6 км, штрихпунктирная кривая - x = 634.25 км.
В п.5.4 приводятся и обсуждаются результаты численных расчетов диаграмм рассеяния (диаграмм направленности) для плоских и цилиндрических Е (или Н) поляризованных волн двухмерными вогнуто-выпуклыми идеально проводящими незамкнутыми одиночными экранами и системой, состоящей из двух экранов (Удвух зеркальнаяФ антенна). Расчет диаграмм рассеяния для таких случаев (диаграмм направленности) осуществлялся с помощью разработанного алгоритма, в основе которого лежат предложенный метод продвинутых граничных условий (в действительном и комплексном вариантах) и использование вэйвлетов Хаара и Батла-Лемарье в качестве системы базисных функций при численном решении интегральных уравнений (ИУ) Фредгольма первого рода с гладким ядром. Показано, что использование вэйвлетов Батла-Лемарье в качестве системы базисных функций приводит к снижению точности решения граничной задачи на один-два порядка по сравнению с использование вэйвлетов Хаара. Пример расчета нормированной диаграммы рассеяния для экрана, контур которой представлен на рис.9 а, представлен на рис. 9 б (параметры контура: kA = 60, kT = 40, kxk = 40, kxn =-45 ; координаты источника: kx0 = -40, ky0 = -59 ). В этом случае, как известно, каустика ГО-лучей имеет точку, уходящую на бесконечность.
y x Источник цилиндриче ской волны а) б) Рис. 9 Геометрия контура - а); структура нормированной диаграммы рассеяния - б);
Подробно исследовано изменение структуры диаграммы рассеяния в зависимости от положения источника, его поляризации и формы экрана в условиях, когда образуются волны шепчущей галереи. В случае системы, состоящей из двух параболических экранов (малого и большого) - двухзеркальной антенны, расчет диаграммы рассеяния показал существенное влияние малого зеркала на структуру диаграммы направленности в области ее первых боковых лепестков.
В п.5.5 приводятся результаты расчетов диаграмм рассеяния Е и Н поляризованных плоских и цилиндрических волн вогнуто-выпуклым замкнутым контуром больших электрических размеров (до kD 1500, где k - волновое число свободного пространства, D - максимальный размер рассеивающей области) в условиях образования волн шепчущей галереи. В основе полученных результатов лежит использование модифицированного метода дискретных источников. Исследовано влияние типа падающей волны, вида ее поляризации и места расположения источника на структуру рассеянного поля.
В п.5.6 приводятся результаты исследования диаграмм рассеяния и пространственной фокусировки поля цилиндрической волны линзой Веселаго конечных размеров при различных значениях модуля показателя преломления и потерь метаматериала линзы, а так же диаграммы рассеяния гауссова пучка призмой из метаматериала. Расчеты выполнены на основе использования модифицированного метода дискретных источников. Показано существование предельных значений для разрешающей способности такой линзы. Дана геометрооптическая интерпретация результатов расчетов поля в областях его фокусировки и отклонения максимума диаграммы рассеяния призмой. Пример расчета структуры поля в линзе и области фокусировки для случаев показателя преломления nr = -1.0001. nr =-0.5 и nr = -представлены соответственно на рис.10-12. Зависимость разрешающей способности от высоты линзы, представленной на рис. 13 показывает, что сверхразрешения у таких линз не наблюдается.
Рис.10 Структура ближнего поля идеальной Рис.11 Структура ближнего поля линзы линзы Веселаго ( nr =-1.0001- i0.001) Веселаго ( nr = -0.5 - i0.001) kD 1,2,2kb Рис.12 Структура ближнего поля линзы Рис.13 Разрешающая способность линзы Веселаго ( nr =-2.0 - i0.001) Веселаго (1- nr = -0.5; 2- nr -1;3- nr = -2 ) В заключение раздела приведены результаты строгих расчетов диаграммы рассеяния поля гауссова пучка призмой из метаматериала, подтверждающие прогноз В.Г. Веселаго об отклонении максимума рассеянного поля в левую сторону от нормали к поверхности скошенной части призмы.
Основные результаты диссертации.
Сформулированы четыре универсальных постулата пространственно временной геометрической теории дифракции (ВГТД) и получена асимптотика задачи Коши для трехмерного нестационарного волнового уравнения поля радиосигнала в плазме.
Исследованы особенности лучевых семейств, образованных лучами ВГО и ВГТД. Дана классификация критических точек интегральных представлений Фурье для поля сигнала в среде с дисперсией и установлена их связь с лучами ВГО и ВГТД, что позволило сформулировать путь построения эталонных функций, описывающих поведение поля сигнала вблизи сложных каустик лучей ВГО и ВГТД.
Получены равномерные асимптотические формулы для поля сигнала в областях неприменимости ВГТД - вблизи кустик лучей ВГО и ВГТД.
Предложен и разработан асимптотический подход к решению задачи определения параметров модели неоднородности плазмы и ее потерь по искажениям временной формы сигнала ВЗ, частотной характеристике сигнала ВЗ, пространственной структуре поля гармонического сигнала НЗ вблизи границы мертвой зоны.
Разработан комплекс алгоритмов расчета характеристик полей и сигналов, взаимодействующих с плазмой с учетом ее неоднородности, потерь, анизотропии, способа обработки сигнала в приемнике и ДН излучающей антенны.
Предложен метод продолженных граничных условий - универсальный метод решения задач рассеяния гармонических волн ограниченными телами (экранами) сложной геометрии, сводящий исходную задачу дифракции к решению интегральных уравнений (ИУ) Фредгольма первого рода с гладким ядром (совместно с Кюркчаном А.Г.) и на основе вейвлетных базисов разработан эффективный алгоритм решения таких ИУ в квазиоптической области частот. На основе такого алгоритма впервые решена задача рассеяния полей вогнутовыпуклых идеально проводящим цилиндром и экраном в условиях образования волн шепчущей галереи, их фокусировкой и взаимодействием с краями экрана;
Впервые в строгой постановке получено и решение задачи фокусировки волн плоской линзой Веселаго конечных размеров с учетом потерь и значений показателя преломления метематериала. Установлено, что особенности лучевых картин и структура каустики для такой линзы зависят от значений модуля показателя преломления. Исследована структура поля как в ближней зоне линзы - в окрестности точек фокусировки и каустик внутри и вне линзы, так и в ее дальней зоне. Показана принципиальная ограниченность разрешающей способности таких линз.
Таким образом, в работе разработана теория временной геометрической теории дифракции (ВГТД), позволяющая интерпретировать, предсказывать и рассчитывать дисперсионные искажения сигналов конечной длительности; предложен и разработан универсальный метод решения задач рассеяния волн на компактных телах и экранах сложной формы - метод продвинутых граничных условий, и исследованы новые явления и процессы, связанные с искажениями частотных и пространственных спектров широкополосных сигналов в средах с частотной и пространственной дисперсией и искусственных средах (метаматериале). В целом, совокупность предложенных и разработанных автором методов можно квалифицировать как новое крупное достижение в развитии важного научного направления: создание и развитие асимптотической теории, описывающей дисперсионные искажения сигналов конечной длительности в диспергирующих средах, а так же изучение распространения и дифракции сигналов в средах с частотной и пространственной дисперсией.
Основные публикации по теме диссертации:
1. Анютин А.П., Распространение плоского и цилиндрического радио импульса в плоско- слоистой плазме, / А.П. Анютин // Труды МЭИ, сер. Электроника и радиотехника, - 1974, - вып. 192, - С. 15-19.
2. Анютин А.П., Асимптотическая теория распространения модулированных радиосигналов в неоднородной диспергирующей среде, / А.П. Анютин, Ю.И. Орлов // Труды МЭИ, вып.
194, - 1974, - С. 99-102.
3. Анютин А.П., Равномерная асимптотика поля прямоугольных радиоимпульсов в магнии- тоактивной плазме, /А.П. Анютин, Ю.И. Орлов// Изв. ВУЗов, сер. Радиофизика, - 1976, - Т. 19, N 3, - C. 335-342.
4. Орлов Ю.И., О модификации лучевого метода в теории распространения радиосигналов в неоднородных диспергирующих средах, / Орлов Ю.И., Анютин А.П. // Изв. ВУЗов, сер.
Радиофизика, - 1974, - Т.17, N 9, - С.1369-135. Анютин А.П., Об особенностях искажении радиосигналов в неоднородном линейном слое плазмы, / А.П. Анютин, Ю.И. Орлов // Изв. ВУЗов, сер. Радиофизика, - 1976, - Т. 19, N 4, - С. 495-504.
6. Анютин А.П., Пространственно-временная геометрическая теория дифракции частотномодулированных радиосигналов в однородной диспергирующей среде, / А.П. Анютин, Ю.И. Орлов // Радиотехника и Электроника, - 1977, - Т. 22, N 10, - C. 2082-2090.
7. Анютин А.П., Асимптотическая теория распространения радиосигналов в неоднородной плазме, / А.П. Анютин // Распространение радиоволн в Ионосфере, - М.: ИЗМИРАН СССР, - 1978, - C. 19-27.
8. Анютин А.П., ВГТД видеосигналов в однородной среде, /А.П. Анютин// Тезисы докладов ХШ Всесоюзной конференции по распространению радиоволн, - М.: Наука, - 1981, - Т. 2, - C. 251-254.
9. Анютин А.П., Об одном подходе к решению задач нестационарной диагностики плазмы, /А.П. Анютин// Распространение декаметровых волн.- М.: ИЗМИРАН СССР, -1977, - C. 93-96.
10. Анютин А.П., Об отражении АМ-ЧМ сигналов от неоднородной плазмы, / А.П. Анютин // Распространеие радиоволн в ионоcфере, - М.: ИЗМИРАН СССР, - 1978, - С. 29-36.
11. Анютин А.П., Об особенностях искажений AM и ЧИ сигналов, отраженных от слоя Эпштейна, / А.П. Анютин // Изв. ВУЗов, сер. Радиофизика, - 1979, - Т. 22, N 6, - C. 703-710.
12. Анютин А.П., О расчете и интерпретации искажений сигналов, отраженных от неоднородной плазмы, В кн.: Вопросы распространения радиоволн в высоких и средних широтах, / А.П. Анютин // М.: ИЗМИРАН СССР, - 1979, - C. 153-163.
13. Анютин А.П., О влияние неоднородности плазмы на характер искажений АМ-ЧМ сигна- лов, / А.П. Анютин // Изв. ВУЗов, сер. Радиофизика,- 1980, - Т. 23, N 5, - C. 523-528.
14. Анютин А.П., О лучевом описании искажений радиосигналов в плоско слоистой плазме, / А.П. Анютин // Геомагнетизм и Аэрономия, - 1980, - Т. 20, N 3, - C. 555-557.
15. Анютин А.П., Некоторые обобщения и приложения BГTД, / А.П. Анютин, Б.Е. Кинбер // Волны и дифракция, - М.: ИРЭ АН СССР, - 1981, - Т. 3, - С. 276-279.
16. Анютин А.П., О расчете и интерпретации искажений сигналов в ионосферных волноводах, / А.П. Анютин, А.А. Апуневич // Тезисы докладов Ш Всесоюзной конференции по взаимодействию электромагнитных излучений с плазмой, - Алма-Ата, - 1982, - С. 155156.
17. Анютин А.П., Искажения сигналов в окрестности пространственной каустики, / А.П.
Анютин // В кн.: Тезисы докладов Х1У Всесоюзной конференции по распространению радиоволн, - М.: Наука, -1984, - Т. 2, - С. 247-248.
18. Анютин А.П., Искажения сигналов в окрестности каустики, / А.П.Анютин // Геомагнетизм и Аэрономия, - 1985, - Т. 25, N 6, - С. 935-940.
19. Анютин А.П., Дисперсионное расплывание сигналов с частотой близкой к критической частоте плазмы, / А.П. Анютин, В.В.Кольцов // Волны и дифракция,- М.: Наука,-1985, - Т. 2, - С. 478-481.
20. Анютин А.П., Численное и асимптотическое моделирование импульсных сигналов на выходе приемника при вертикальном радиозондировании ионосферы, / А.П. Анютин // Цифровые ионозонды и их применение, - М.: ИЗМИРАН СССР, - 1985, - С. 120-127.
21. Анютин А.П.,ВГТД интерферометрия, /А.П. Анютин И.О.Порохов// Волны и дифракция, - М.: АН СССР, - 1985, - Т. 1, - С.50-53.
22. Анютин А.П., Нестационарная диагностика параметров невозмущенной F области ионосферы, / А.П. Анютин, И.О.Порохов // Цифровые ионозонды и их применение, - М.: ИЗМИРАН СССР, - 1985,- С.128-132.
23. Анютин А.П., Равномерные модификации метода ВГТД в случае произвольной диспергирующей среды и каустик ВГО и ВГТД лучей, / А.П. Анютин // Дифракция и распространение волн, - М.: МФТИ, - 1985, - С. 32-36.
24. Анютин А.П., Об отражении сигналов от однородного полупространства магнитоактивной плазмы с потерями, / А.П. Анютин // Распространение декаметровых радиоволн, - М.:
ИЗМИР АН СССР, - 1979, - С. 97-99.
25. Анютин А.П., Апуневич А.А., О распространении сигналов в диспергирующей среде, / А.П. Анютин, А.А Апуневич // В кн.: Машинное проектирование устройств и систем СВЧ., - М.: МИРЭА, - 1982, - С. 205-213.
26. Анютин А.П., О расчете сигналов, взаимодействующих с однородным плазменным слоем, / А.П. Анютин // В кн.: Распространение декаметровых радиоволн, - М.: ИЗМИР АН СССР, - 1979, - С. 160-167.
27. Анютин А.П., О численном моделировании дисперсионных искажений сигналов в волноводах, / А.П. Анютин // В кн.: Краткие тезисы докладов Всесоюзной конференций "Машинное проектирование устройств и систем СВЧ, - Тбилиси, - 1979, - С. 264-265.
28. Анютин А.П., О влиянии потерь на распространение сигналов в плазме, / А.П. Анютин, А.А Апуневич // В кн.: Машинное проектирование устройств и систем СВЧ, - М.: МИРЭА, - 1979, - С. 228-234.
29. Анютин А.П., К вопросу о распространении сигналов в однородной диспергирующей среде, / А.П. Анютин, А.Н. Щербицкий // В кн.: Машинное проектирование устройств и систем СВЧ, - М.: МИРЭА, - 1979, - С. 225-227.
30. Анютин А.П., О возможности определения поглощения в отражающей области ионосферы по измерению поля вблизи зоны молчания, / А.П. Анютин, В.Г. Галушко, Ю.М. Ямпольский // Изв. ВУЗов, сер. Радиофизика, - 1985, - Т. 28, N 2, - С. 247-249.
31. Анютин А.П., ВГТД интерферометрия ионосферной плазмы, / А.П.Анютин, И.О. Порохов // Всесоюзный семинар НТО им. Попова "Исследование неоднородностей ионосферы радиофизическими методами", Душанбе, -1986, - С. 123-132. Анютин А.П., О расчете поля сигнала, отраженного возмущенной окрестностью максимума F области ионосферы, / А.П. Анютин // В кн.: Дифракция и распространение волн в неоднородных средах, - М.: МФТИ, - 1987, - С. 57-72.
33. Анютин А.П., Искажения амплитудно-модулированных радиоимпульсов вблизи критической частоты F области ионосферы, / А.П. Анютин, В.В. Кольцов // В кн.: Распространение радиоволн в ионосфере, - М.: ИЗМИРАН СССР, - 1986, - С. 69-73.
34. Анютин А.П., Решение задач дифракции методом вспомогательных токов с использованием вэйвлетных базисов, / А.П. Анютин, А.Г. Кюркчан // Радиотехника и Электроника, - 2002, - Т.47, N 6, - С. 688-693.
35. Кюркчан А.Г., Метод продолженных граничных условий и вейвлеты, / А.Г. Кюркчан, А.П. Анютин // Доклады Академии Наук, - 2002, - Т. 385, N 3, - С. 309-313.
36. Кюркчан А.Г., Вейвлеты и метод вспомогательных токов, / А.Г. Кюркчан, А.П. Анютин // Доклады Академии Наук, - 2002, - Т. 383, N 5, - С. 612-616.
37. Анютин А.П., О новой модификации метода дискретных источников, / А.П. Анютин, А.Г.
Кюркчан, С.А. Минаев // Радиотехника и Электроника, - 2002, - Т.47, N 8, - С. 955-960.
38. Anioutine A.P., Application wavelets technique for solving an integral equations of the method of auxiliary currents, / A.P. Anioutine, A.G. Kyurkchan // JQSRT, - 2003, - Vol. 79 - 80, - pp. 495-508.
39. Anioutine A.P., About universal modification to the method of discrete sources and its application, / A.P. Anioutine, A.G. Kyurkchan, S.A. Minaev // JQSRT, - 2003, - Vol. 79 - 80, - pp. 509-520.
40. Анютин А.П., ММДИ в задачах рассеяния цилиндрическими структурами больших электрических размеров, / А.П. Анютин, С.А. Минаев // Радиотехника и электроника, - 2003, - Т. 48, N 11, - С. 1343-1348.
41. Анютин А.П., Об отражении широкополосных сигналов плоскослоистой диспергирующей средой с потерями, / А.П. Анютин // Электромагнитные волны и электронные системы, - 2004, - Т.9, N 1, - C. 27-32.
42. Анютин А.П., Решение задач теории дифракции и антенн с использованием метода продолженных граничных условий и техники вейвлетов, / А.П. Анютин, А.Г. Кюркчан // Радиотехника и Электроника, - 2004, - Т. 49, N 1, - C. 15-23.
43. Анютин А.П., О расчете рассеяния волн шепчущей галереи, / А.П. Анютин, В.И. Стасевич // Доклады Академии Наук, - 2004, - Т.399, N 4, - С. 1-4.
44. Анютин А.П., О рассеянии Е и Н поляризованных волн вогнуто-выпуклыми экранами, / А.П. Анютин, В.И. Стасевич // Радиотехника и Электроника, - 2004, - Т. 49, N 11, - С. 1338-1343.
45. Anyutin A.P., Scattering of E polarized whispering-gallery mode from concave boundary, / A.P. Anyutin and V.I. Stasevich // Programa y Resumenes, XIV Simposio Internacional de Matematicos Aplicados a las Ciencias, San Jose, (17-20 Febrero), Costa Rica, - 2004, - pp. 18-19.
46. Anyutin A.P., MMDS and 2D scattering problem by cylindrical structure with piece-wise boundary, / A.P. Anyutin and V.I. Stasevich // Revista de Matematica: Teoria y Aplicaciones, - 2006, - Vol. 13, N1, - pp.37-47.
47. Anyutin A.P., Scattering of E polarized whispering-gallery mode from concave boundary, / A.P. Anyutin and V.I. Stasevich // Revista de Matematica: Teoria y Aplicaciones, - 2005, - Vol. 12, N1-2, - pp. 121-128.
48. Anyutin A.P., Scattering of whispering-gallery mode from concave-convex boundary, / A.P.
Anyutin and V.I. Stasevich // Proceedings of International Symposium on Electromagnetic Theory, - 2004, - Vol. 1, - pp. 543-549. Anyutin A.P., Scattering of E and H polarized waves by covered cylindrical structures, / A.P.
Anyutin and V.I. Stasevich // Proceedings of ELS-NSP 8th International Conference on Electromagnetic and Light Scattering by Non-spherical Particles, - 2005, - pp.5-8.
50. Anioutine A.P., About 3D solution of diffraction problems by MMDS, / A.P. Anioutine, A.G.
Kyurkchan, S.A. Minaev // Proceedings of ELS-NSP 8th International Conference on Electromagnetic and Light Scattering by Non-spherical Particles, - 2005, - pp.9-12.
51. Anyutin A.P., About scattering problems for 2D and 3D bodies with wedges in kD>>1 region, / A.P. Anyutin // Programa y Resumenes, XV Simposio Internacional de Matematicos Aplicados a las Ciencias, - 2006, - pp.28-29.
52. Anyutin A.P., On a new>
53. Anyutin A.P. About focusing problem for Veselago lens with finite size, /A.P. Anyutin // XV1 Simposio Internacional de Metodos Matematicos Aplicados a las Ciencias, San Jose, 19-22 Febrero, Costa Rica, - 2008, - pp. 31-32.
54. Anioutine A.P., Wavelets in 2D and 3D scattering and propagation problems, / A.P. Anioutine, A.G. Kyurkchan // Proceedings of XXVII-th General Assembly of URSI, Maastricht, - 2002.
55. Anioutine A.P., Application of the wavelets in 2D scattering problems, / A.P. Anyutin and V.I.
Stasevich // WSEAS Transactions on Systems, - 2003, - Vol. 2, Issue 1, - pp. 74-76.
56. Анютин А.П., О рассеянии Е и Н поляризованных волн вогнуто-выпуклым цилиндром больших электрических размеров, / А.П. Анютин, В.И. Стасевич // Радиотехника и Электроника, - 2004, - Т. 49, N 12, - С. 1421-1426.
57. Анютин А.П., Рассеяние цилиндрических Е поляризованных волн многоугольными цилиндрическими структурами, / А.П. Анютин, В.И. Стасевич // Радиотехника и Электроника, - 2006, - Т.50, N 12, - С. 1334-1341.
58. Anyutin A.P., Scattering of E and H polarized waves by covered cylindrical structures, / A.P.
Anyutin and V.I. Stasevich // JQSRT, - 2006, - Vol. 100, N1-3, - pp.16-25.
59. Анютин А.П., Пространственно-временная геометрическая теория дифракции сигналов в диспергирующей среде и ее некоторые обобщения., /А.П. Анютин // Электромагнитные волны и электронные системы, - 2007, - Т.12, N 9, - С. 26-38.
60. Anioutine A.P., About 3D solution of diffraction problems by MMDS, / A.P. Anioutine, A.G.
Kyurkchan, S.A. Manenkov, S.A. Minaev // JQSRT, - 2006,- Vol.100, N 1-3, - pp.26-36.
61. Анютин А.П., Моделирование характеристик зеркальных антенн на основе метода продолженных граничных условий, /А.П.Анютин, А.Г.Кюркчан, Н.И.Смирнова, // Электромагнитные волны и электронные системы, - 2007, - Т. 12, N 8, - С. 63-70.
62. Анютин А.П., О фокусировке и рассеянии поля цилиндрической волны линзой Веселаго конечных размеров, /А.П. Анютин // Радиотехника и Электроника, - 2008, - Т. 53, N 4, - С.413-422.
63. Анютин А.П., О рассеянии волн призмой из метаматериала с потерями, /А.П. Анютин // Радиотехника и Электроника, - 2008, - Т. 53, N 5, - С 570-574.