ОГЛАВЛЕНИЕ
Введение.........................................................................3
Глава 1. Влияние поляризации на отражение волн..................................12
1.1. Варианты задания ограниченного слоя, при которых уравнение для падающей плоской волны ТМ или ТЕ поляризации сводится к 11 У.............................12
1.2. Исследование решения задачи о падении плоской 7М-волны на переходный
слой............................................................................22
1.3. Точное аналитическое решение задачи распространения волны горизонтальной поляризации через переходный слой с переменным масштабом........................28
1.4. Выводы.....................................................................63
Глава 2. Отражение плоских волн горизонтальной поляризации от плазменных слоёв
с переменным масштабом изменения плазменной частоты.............................64
2.1. Строгое решение задачи о поле волны горизонтальной поляризации
для плазменного слоя с переменным масштабом.....................................64
2.2. Исследование частотной зависимости модуля коэффициента отражения горизонтально поляризованной плоской волны от симметричного плазменного
слоя............................................................................76
2.3. Вид плазменного переходного слоя, для которого решение уравнения для волны горизонтальной поляризации выражается через вырожденные ГГФ.....................83
2.4. Выводы.....................................................................86
Глава 3. Распространение 7Л/-волны в плазменном слое с максимумом электронной концентрации при малых потерях .................................................87
3.1. Задача о падении ТМ-волны на симметричный плазменный слой...................90
3.2. Обсуждение полученных результатов..........................................99
Заключение.....................................................................101
Библиографический список использованной литературы.............................102
Приложение 1...................................................................108
Приложение 2. Трансформация электромагнитной волны в нестационарной плазме.. 109
3
ВВЕДЕНИЕ
Задачи распространения волн в неоднородных средах являются предметом исследования разных разделов физики (механики, квантовой механики, радиофизики, огггики). В большинстве случаев математические модели изучаемых явлений сводятся к линейным обыкновенным дифференциальным уравнениям (ЛОДУ) второго порядка, которые для многих задач теории колебаний оказываются схожими и можно говорить об общих способах решения этих уравнений. Изучению этих уравнений посвящено огромное число работ. По теории распространения электромагнитных волн в детерминированных неоднородных средах имеется целый ряд монографий [1-13] и журнальных статей [14-31], а также диссертаций [32-34], в которых рассмотрены точные и приближенные решения многих задач распространения волн. Интерес к этим задачам возник в начале XX века, после того как эксперименты Кеннели в Америке и Хевисайда в Европе доказали возможность отражения волн от ионосферы [35].
Распространение радиоволн в изотропной среде зависит от дисперсионных свойств и степени неоднородности среды, а также от вида поляризации волны. При рассмотрении слоисто-неоднородных сред (плоскослоистых, цилиндрически-слоистых, сферически-слоистых) уравнения Максвелла распадаются на две независимые пары, которые соответствуют двум видам поляризации (вертикальной и горизонтальной). Как известно [7], простейшим типом электромагнитной волны является плоская волна. Сферическая или цилиндрическая волна на большом расстоянии от источника может в силу малого искривления участков фронта рассматриваться в ограниченной области пространства как плоская. Это позволяет в ряде случаев при изучении отражения волн от слоёв пользоваться выражениями для коэффициентов отражения плоских волн. Линейно поляризованная волна (распространяющаяся в плоскослоистой среде се(г)), у которой вектор электрического поля лежит в плоскости распространения, проходящей через прямую (г) и направление распространения падающей волны, называется вертикально поляризованной (или ТМ), а с вектором, перпендикулярным плоскости, - горизонтально поляризованной (или 7Е).
Изучение отражения волн от неоднородностей среды является одной из широкого круга задач, касающихся теории распространения волн и энергетического расчета радиолиний. В случае плоскослоистой среды се(г) и гармонической зависимости поля от времени в результате Фурье-преобразования поля по двум другим пространственным переменным из уравнений Максвелла (с учётом линейности уравнений), являющихся (в дифференциальной форме) уравнениями в частных производных, для комплексных амплитуд полей получаются ЛОДУ [4]. Рассмотрение задач, для которых известны точные анали-
4
тические решения (ТАР) (т е. решения, выраженные через специальные функции математической физики) этих уравнений, имеет весьма большое значение. Дело в том, что вытекающие из точных решений выражения для коэффициентов отражения и прохождения волн раскрывают важные закономерности в зависимостях эпгих коэффициентов от частоты волны, угла падения волны на слой, параметров слоя, а также от поляризации волны. Радиофизика, как наука, имеет дело с определенными математическими моделями. Наибольший интерес представляет задание зависимости е от расстояния г посредством аналитических функций, для которых известны ТАР уравнений Максвелла. Такие зависимости е(г) при рассмотрении других качественно сходных моделей среды можно рассматривать в качестве эталонных [20,25].
В задачах радиофизики, описываемых дифференциальными уравнениями, параметры задачи, входящие в коэффициенты этих уравнений, не могут быть измерены точно, и они могут изменяться под влиянием различных возмущающих факторов. В связи с этим большое значение имеет теорема [36] об условиях (Липшица), при которых решение ЛОДУ непрерывно зависит от параметров, входящих в коэффициенты уравнения. Условия этой теоремы в практических задачах обычно выполняются, за исключением случаев, связанных с некоторыми предельными переходами по параметру либо параметрам [37]. Так, например, в случае, когда дифференциальное уравнение имеет малый параметр при старшей производной при частоте волны о —► оо, модуль коэффициента отражения в пределе может терпеть скачкообразное изменение. Возможна также ситуация, когда решение, являющееся непрерывным по каждому из параметров, не является непрерывным по их совокупности. В этом случае значение коэффициентов отражения и прохождения волны можег зависеть от порядка предельного перехода по параметрам. Учёт одного из них, как правило, меняет структуру особых точек дифференциального уравнения. Такая ситуация имеет место в задаче об экранировании ТМ - волны в симметричном плазменном слое с малыми потерями и на частоте волны, близкой к максимальной плазменной частоте слоя [25, 33]. При этом решение зависит от порядка стремления параметров задачи -эффективной частоты столкновений и угла падения волны - к нулю. Реально эти предельные значения никогда не достигаются. Значения этих параметров могут находиться в окрестности их предельных значений, и представляет тггерес исследование коэффициентов отражения и прохождения в более широкой области изменения параметров.
Хорошо известно, что впервые законы отражения и преломления волн в электромагнитной теории света были угаданы Френелем в 1823 году. Затем этими вопросами занимался Рэлей, Жамен, Коши. Ряд вопросов о согласовании теории и эксперимента, а также
5
обсуждение возможных причин, ответственных за некоторое их несогласование полно освещён в лекциях [38] в разделе, посвящённом некоторым вопросам теории колебаний.
Поскольку механика, а вместе с ней и теория дифференциального и интегрального исчисления появилась значительно раньше радиофизики, то нет ничего удивительного в том, что многие модели неоднородных сред в радиофизике были заимствованы из механики. ЛОДУ второго порядка описывают в радиофизике, например, распространение электромагнитных полей, а в механике - колебания с сосредоточенными или распределенными параметрами и одной степенью свободы. Разумеется, аналогия между механикой и радиофизикой прослеживается и при рассмотрении задач, математические модели которых даются в виде систем ЛОДУ.
Теория линейных дифференциальных уравнений второго порядка с аналитическими коэффициентами, которая берёт начало в XIX веке, привела к получению всех важнейших специальных функций [39-42]. Важную роль среди них занимают уравнения, все особые точки которых регулярны [43]. Гипергеометрическое уравнение (ГГУ) является одним из таких наиболее изученных уравнений с тремя регулярными особыми точками. Коэффициенты в функциональных соотношениях (связывающих линейно-независимые решения ГГУ в различных интервалах изменения независимой переменной) выражаются через гамма-функцию [44], теория которой в основном была завершена к концу XIX - к началу XX века. По теории ЛОДУ существует обширная литература [45-60]. Такой большой интерес к ним связан с их важностью в прикладной математике и отсутствием метода нахождения общего решения произвольного ЛОДУ второго порядка. Вместе с тем в течение последних тридцати лет теория дифференциальных уравнений (прежде всего по таким вопросам, как устойчивость решений, применение методов функционального анализа, разработка новых методов приближенного решения уравнений) развивалась настолько интенсивно, что претендующее на полноту сопоставление с литературой заняло бы более сотни страниц списка литературы (заметим, что в книге [61] библиография занимает 140 страниц, а с момента её выхода появился целый ряд монографий и статей)
Одной из самых первых (обзорных) работ об интегрировании ЛОДУ второго порядка с переменными коэффициентами, которое встречается в теории колебаний, является работа [57]. В ней рассмотрены функции /(х), для которых дифференциальное уравнение
Ув+Ро/(х)У = 0 > (1)
(где Р0 - параметр), посредством замен зависимой и независимой переменных сводится а) к уравнению с постоянными коэффициентами, б) к уравнению Бесселя и в) к ГГУ.
атити тш .
/(*> = ~~77~. :—гг- • <2>
6
I Д + огех р(|--—-----5-)
В случае а) /(*)-__. вслучаеб) т = ——-±Щ±Ш-,га,
к,т,п,д,а,/?-произвольные величины. В 157] утверждалось (без доказательства), что в случае в)
а + Ьх + сх2 + дЬ:3 (Л + тх + пх2)3
Однако, как показано в [47], уравнение Римана, которое дробно-линейным преобразованием независимой переменной и линейной заменой зависимой переменной сводится к
___ _ „ ... . а + Ьх + сх2 ,
ГГУ, имеет канонический вид (1), где /(х) =-------—, так что эта функция совпа-
(к +тх + пх )
дает с (2) только при с/ =0.
Решения уравнения (1) можно применить в задаче о нормальном падении ТЕ волны на слой, функция диэлектрической проницаемости е которого зависит от безразмерного расстояния х по закону /(х). Некоторые из этих зависимостей е (х)= /(х) одними из первых были рассмотрены в радиофизической литературе в задачах о распространении плоских электромагнитных волн [23].
Наиболее полный обзор функций е(г), для которых решения уравнений Максвелла для волны ТМ поляризации имеют ТАР в терминах ГТФ и вырожденных ГГФ, сделан в работе [24]. В интервале г е (-со,-ко) почти все (кроме двух вариантов) рассмотренные в ней зависимости являются либо неограниченными, либо одно из их предельных значений при г —> -оо (или г -> +оо ) равно нулю. На этих двух вариантах е(г) следует подробнее остановиться. Одна из этих функций задаётся выражением
е(2) = Ш2(аг + Ь), (3)
где К, а, Ь- постоянные. Для этой модели среды в [15, 22] рассмотрены ТАР уравнения и для ТЕ волны. Другая зависимость £:(г), для которой было найдено ТАР для 7М поля, задавалась в неявном виде е = е(С00). однако, об ограниченности этой функции £г(^(г)) в работе [24] не упоминалось. В ней не рассматривался и вопрос о возможном виде (или видах) слоя при той или иной взаимно однозначной связи между переменными £ и г . Оказывается (это рассмотрено в диссертации), что при определенной такой связи эта функция будет описывать переходный слой.
Что касается систематического анализа функций, для которых известны ТАР уравнения для ТЕ волны, то этот вопрос был рассмотрен в [15,22]. В работах [14, 15, 22-24] приведены зависимости е(г), для которых ЛОДУ для полей горизонтальной и верти-
кальной поляризации сводятся к изученным специальным функциям. Эти работы имеют радиофизическую специфику, в том числе связанную и с различными вариантами задания среды распространения волны (в какой-либо выбранной системе координат). В общем же математическом плане книги [49, 56] охватывают более широкие виды линейных дифференциальных уравнений второю порядка с переменными коэффициентами, которые заменами переменных сводятся к уравнениям Бесселя, Лежандра, Уиттекера, Матье, Ламе и гипергеометрическому. В них проведено более полное по сравнению с [22-24] исследование тех случаев, когда ЛОДУ второго порядка приводятся к изученным уравнениям. Однако не все результаты исследований, проведённых в этих работах, можно приложить к радиофизическим задачам. Это связано с тем, что в ЛОДУ для полей параметры задачи (частота волны, угол падения волны на слой, либо производная Шварца (при рассмотрении ТМ- волны)) входят в уравнение определенным образом, что заметно сужает применение многих уравнений, рассмотренных в [49], к теории распространения электромагнитных волн. Кроме того, в [49] проанализированы только линейные преобразования зависимой переменной. Заметим, что некоторые результаты этой работы были получены в [56] методом факторизации, обобщенным на случай переменных коэффициентов.
Обсудим коротко приближенные методы решения задач о падении плоской волны какой-либо поляризации на слой, необходимость привлечения которых связана с отсутствием известных ТАР уравнений для полей. Все эти приближенные методы [60-72] связаны с наличием малого параметра задачи. Это методы: частичных отражений (частный случай метода последовательных приближений) [1], ВКБ [62, 64], метод геометрической оптики [67], метод фазовых интегралов [9,62]. Границы применимости метода геометрической оптики были расширены [73]. Здесь следует упомянуть предложенный в 1952 году метод осреднения функциональных поправок [74], применяемый для приближенного решения интегро-дифференциальных уравнений. Ссыпок на его применение в задачах радиофизики автор не нашёл, но этот метод, возможно, будет применён и при решении задач о распространении электромагнитных волн в неоднородных средах. Большая библиография в [74] указывает на его различные (математические) видоизменения и обобщения. В случае среды, характеризующейся бесконечно дифференцируемой функцией диэлектрической проницаемости, а также при отсутствии точек поворота на вещественной оси модуль коэффициента отражения при стремлении частоты падающей волны со —> <х> экспоненциально мал [9]. Таким образом, он не может быть получен в рамках лучевого метода, где для него получается степенное убывание с ростом частоты. Приближенные методы решения задачи распространения волн тесно связаны с математической теорией устойчивости [75] и с теорией асимптотического (в редких случаях сходящегося) разложе-
8
ния (в монографии [64] решения уравнений с малым параметром при старшей производной представлены в виде сходящихся рядов). Укажем также на метод эталонного уравнения (МЭУ) (который в отечественной физико-математической литературе связывают с работами А. А. Дородницына, а в западной математической литературе (вероятно учитывая историю возникновения, а не саму идею метода, безусловно принадлежащую А. А. Дородницыну, более известный как преобразование Лангера) [60,70]. Этот метод приближенного решения ЛОДУ при наличии большого параметра задачи нашёл широкое применение в задачах распространения электромагнитных волн в неоднородных (изотропных и анизотропных) средах (см. работу [20], а также большое число работ Г. И. Макарова, на которые в ней приведены ссылки). Этот метод является одним из эффективных способов нахождения асимптотического поведения специальных функций математической физики.
С математической точки зрения возможна ситуация, когда слой настолько медленно выходит (при z —> ±оо) на вакуум (или другую однородную среду), что полученное решение при z -> ±00 не будет иметь вид плоских волн. Математическая формулировка условий, при которых это имеет место, дана в работе [61]. Два варианта такого слоя приведены в диссертации (один вариант - в первой главе, а другой - во второй). Указанные выше асимптотические методы часто приводят к решению в виде расходящегося ряда.
В радиофизической литературе интерес к точно решаемым задачам о распространении волны в среде никогда не пропадал. Трудность этих задач в том, что основной метод интегрирования дифференциальных уравнений - это введение удобных замен зависимой и независимой переменных, преобразующих уравнение к простейшему виду, но для нахождения этих замен нет общего правила. Заранее не ясно, к какому уравнению целесообразнее свести исходное дифференциальное уравнение. При этом количество и характер его особых точек не всегда может навести на успешное преобразование. В математической литературе линейно-дифференциальные и другие более сложные (нелинейные) преобразования зависимой переменной (которые могут изменить количество особых точек уравнения) рассмотрены не так подробно, как линейные преобразования. При рассмотрении нелинейных преобразований зависимой переменной тем более нет общего правила.
Из обзора вышеперечисленных работ можно сделать вывод, что в классе аналитических функций e(z), вещественных и непрерывных при действительных значениях г, а также принимающих произвольные конечные значения при z -* ±оо, существует только одна зависимостьc(z), допускающая ТАР уравнения для 77С-волны - это слой Эпштейна [1, 14]. Однако даже для простейшего случая переходного слоя Эпштейна решение уравнения для ТМ- волны не получено. С помощью замен независимой и зависимой перемен-
9
ных его можно свести к уравнению с четырьмя особыми точками [58, 59], теория которого к настоящему времени еще не завершена. Функциональные соотношения, связывающие асимгпотики линейно-независимых решений произвольного уравнения с четырьмя правильными особыми точками на сегодняшний день не известны. Полученные в [59] решения такого уравнения (в случае вещественного коэффициента уравнения) в виде очень громоздких рядов, сходимость которых имеет место только в отдельных областях, не дают возможности получить простые (т.е. в конечном виде) выражения для коэффициентов отражения и прохождения волны даже в случае вещественного слоя.
В литературе исследовано все ещё недостаточное количество точно решаемых моделей, чтобы провести какую либо классификацию закономерностей коэффициента отражения, поэтому представляет интерес дальнейшее изучение закономерностей отражения волн при рассмотрении некоторых новых моделей среды. С целью их изучения предприняты исследования, результаты которых изложены ниже.
В диссертационной работе представлены некоторые новые ТАР задач о распространении плоских гармонических волн в безграничных изотропных плоскослоистых средах.
В математическом отношении эти задачи сводятся к решению дифференциальных уравнений с соответствующими условиями на бесконечности [7, 17].
В диссертации рассмотрены такие зависимости диэлектрической проницаемости е от расстояния г, для которых ТАР уравнений Максвелла для гармонических полей горизонтально и (или) вертикально поляризованных волн сводятся к гипергеометрическим и вырожденным гипергеометрическим функциям. В диссертации впервые, исходя из ТАР уравнений Максвелла для ТМ и ТЕ- волны, рассматривается задача о влиянии поляризации на особенности отражения волн для положительной и ограниченной функции е(г) специального вида, моделирующего переходный слой. Для него автором найдены коэффициенты отражения R и прохождения Т плоских гармонических волн и исследованы зависимости коэффициентов отражения от частоты волны и угла падения волны на слой.
В первой главе диссертации для этого слоя (который в общем случае задается неявно) найдены выражения для коэффициентов отражения и прохождения ТМ- волны. Затем (для этого же слоя) автором получено ТАР уравнения и для ТЕ- волны, что позволяет исследовать влияние поляризации на особенности отражения волн от такого слоя при произвольных значениях частоты волны и угла падения волны на слой. При получении ТАР уравнения для ТЕ- волны применяется линейно-дифференциальное преобразование зависимой переменной уравнения с целью понижения количества его особых точек с четырёх до трёх регулярных особых точек и сведения уравнения (для новой функции) к ГГУ.
- Київ+380960830922