Диссертация (785868), страница 28
Текст из файла (страница 28)
Рассмотрим поле электрических волн, примемH r 0 . Тогдасистема уравнений (А5) примет вид: E 0 sin E Er r sin E j r sin H r ErE r j rHr H j r sin Er sin H rH j rEr rH j rEr (А7)Из равенства нулю радиальной компоненты напряжённости магнитного поляследует соотношение: E sin E (А8)Используя приведенный в диссертационной работе алгоритм для решениясистемы уравнений (А7), удобно ввести вспомогательную функцию U. Всоотношении(А8)составляющиеEиEможновыразитьчерезвспомогательную функцию U:1 U r sin 1 U E r E (А9)259гдеUUrEE1 2Ur sin r1 2Ur r(А10)Подставляя (А10) в систему уравнений (А7), можно выразить компоненты H иH :j Ur j UH r sin H (А11)Подставив (А10) и (А11) в третье уравнение системы (А10), получаем выражениедля радиальной составляющей вектора напряжённости электрического поля: 3U Erk 2Ur 2 2U k 2U Err 2(А12)Вспомогательную функцию U можно получить, выразив компоненты поля черезU в четвёртом уравнении системы (А7): 2U U 1 2U r 2 sin k 2U sin 2 sin 2 rи придать вид: U 2U2U k 2U 0 sin 222 r sin 2 rr sin 1(А13)260Аналогично можно рассмотреть поле магнитных волн.
ПримемЕr 0 . Тогдасистема уравнений (2.2.5) примет вид: E j r sin H r sin E r sin E j r sin Hr rE jrHr H sin H H r r sin H j r sin E r Hr j rErH r (А14)Из равенства нулю радиальной компоненты напряжённости электрического поляследует соотношение:H sinH (А15)Используя алгоритм для решения системы уравнений (А14) удобно ввестивспомогательную функцию V.
В соотношении (А15) составляющие H и Hможно выразить через вспомогательную функцию V:1 V r sin 1 V H r Hгде V Vr(А16)261HH1 2Vr sin r1 2Vr r(А17)Подставляя (А17) в систему уравнений (А14), можно выразить компоненты E иE :j Vr j VE r sin E (А18 )Подставив выражения (А17) и (А18) в шестое уравнение системы (А14), получаемвыражение для радиальной составляющей вектора напряжённости магнитногополя:3VH rk 2Vr 2 2V k 2V H rr 2(А19)Вспомогательную функцию V можно получить, выразив компоненты поля черезV в первом уравнении системы (А14): 2V V 1 2V r 2 sin k 2V sin 2 sin 2 rи придать вид: V 2V2V k 2V 0 sin r sin 2 2 r 2r 2 sin 1(А20)262Функция U однозначно определяет электромагнитное поле электрического типа, афункция V-магнитного типа: 2V k 2V H r r 2 2U k 2U Er r 21 2UE r r21 U E r sin r Hr 0j U H r sin j U H r (А21а)21 V H r sin r Er 0j V E r sin j VE r (А21б)1 2VH r rПроведём интегрирование уравнения (А13) методом Фурье.
Интеграл уравненияпредставим в виде произведения трёх сомножителей:U Rr Подставив это выражение в уравнение (А7) и разделив на Rr , получаемуравнение в полных производных:1 d d 11 d 2 1 d 2Rr k 2 sin d r 2 sin 2 d 2Rr dr 2r 2 sin d (А22)1Так как в уравнении (А22) сумма трёх функций независимых переменных равнапостоянной величине, то его можно разбить на три уравнения аналогично.Умножаем уравнение (А22) на r 2 sin 2 :sin 1 d d 1 d 2 r 2 sin 2 d 2Rr 2 2 2 k r sin 0 sin d d d 2Rr dr 2и приравниваем полученное выражение постоянной величине m2:2631 d d r 2 sin 2 d 2 Rr 2 2 21 d 2 sin k r sin m2 sin d d Rr d 2dr 2(2.2.23)В результате получаются два уравнения:d 2 m 2 0d 2(А24)1 d d r 2 sin 2 d 2 Rr 2 2 2sin k r sin m 2 0 sin d d Rr dr 2(А25)Интеграл уравнения (А24) имеет вид: А sin m А cosm 12(А26)Для определения Rr и разделим уравнение (А25) на sin 2 :1 1 d d r 2 d 2 Rr 2 2m2k r 0 sin 2sin d d Rr dr 2sin (А27)Приравниваем уравнение (А27) к новой постоянной a2 и разбиваем его на двауравнения аналогично (А23):1 1 d d m 2r 2 d 2 Rr 2 2 k r а 2 sin 22sin d d sin Rr dr1 d d 2m 2 0 sin а sin d d sin 2 d 2 cos d 2m 2 а 022sinddsin d 2 Rr 2 а 2 k Rr 022drr(А28)(А29)264Уравнение (А29) приводится заменой переменной Rr kr F r к уравнениюцилиндрических функций (Бесселя):d 2 F r 1 dF r 2 4а 2 1 k F r 0r drdr 24r 2(А30)Введём новое обозначение для преобразования уравнения (А30):а 2 1(А31)Подставляя (А29) в (А30) получаем уравнение Бесселя для F(r):21 2 d 2 F r 1 dF r 2 k F r 0r drdr 2r2 (А32)Изменение поля по радиусу определяется цилиндрической функцией Бесселя,описывающейстоячуюволнуилицилиндрическойфункциейХанкеля,определяющей волну, расходящуюся от начала координат.krRr c kr J 112 krRr c kr H 21 1 2(А33)(А34)Решение уравнения (А32) удобно представить в виде сферических функцийБесселя или Ханкеля с помощью формул:j x x J2x 12265h2 x H22x 12x При большом значении r, сферические функции могут быть замененыасимптотическими выражениями :j x 1 1 cos x x2i 1 ix2h x exПодставляя асимптотические выражения в выражения (А33) и (А34), получаемрешение уравнения (А32) для больших значений r:2Rr c kr J kr сsin kr 11 kr 12 2Rr c kr H 21 kr 2kr e 1kr 2 1 2При большом значении аргумента, что соответствует условию данной задачи,асимптотические выражения функций Бесселя и Ханкеля с достаточнойточностью описывают эти функции.
Уравнение (А28) приводится к уравнениюЛежандра заменой переменной z cos и заменой постоянной разделения (А31):22 1 z 2 d z 2 z dz 1 mz 02 dz 2dz1 z Если m – целое число, то решение этого уравнения можно записать в видеприсоединённой функции Лежандра первого рода: c Pm cos 2(А35)266где Pm cos - присоединённая функция Лежандра первого рода, ν- степени и mпорядка, с2 - коэффициент.Интеграл волнового уравнения (А22) определяется произведением решений(А26), (А34) и (А35):2U c kr Hkrc Pm cos А sin m А cosm 1212 1 2(А36)Определение переменных (по другой терминологии постоянных) разделениядолжно соответствовать физическому процессу. При определении постояннойразделения по азимуту m поле должно изменяться с периодом 2.
Поэтому mпринимает только целочисленные значения m=0,1,2,3… Однако, индексы ν вобщем случае не будут целыми: они определяются из граничных условий:U 0иV0при Индексы для электрических волн определяются из уравнения:Pm cos 0Здесь через функцию Pm cos обозначена присоединённая функция Лежандра,принимающая при θ=0 конечные значения. Она выражается через функциюЛежандра P cos , удовлетворяющую условию:P 1 1267Поскольку ν не является целым числом, то функция P не будет полиномом,поэтому уравнению P m cos 0 удовлетворяет последовательность чисел m1, m2... mn , зависящих от угла γ.Первоначально определим общее решение, получающееся для волн типа Emn иHmn. Если возьмём число mn и соответствующую функцию U, то по формуле(А36) получим электрическую волну Emn. Если возьмём число mn исоответствующую функцию V, то получим магнитную волну Hmn.Компоненты поля типа Е определяются из уравнений (А21а).
Для записикомпонент векторов электрического и магнитного поля удобно принять обозначенияAmn c c А12 1Bc c А :mn 1 2 2Amn mcosm d Ekr H 2kr Pm cos ,mnr sin dr 1 mnmn 2Bmn msin m d kr Pm cos kr H 2r sin dr 1 mnmn 2 dА sin m d E mn kr H 2kr ,mn1rdr d mn 2 mcos P mn dB cosm d mn kr H 2kr1rdr d mn 2 mcos P mn268Amn mn mn 12Er ,mn sin m kr Hkr Pm cos 1 mnr2 mn 2Bmn mn mn 12mcosm kr Hkr Pcos 1 r2 mn mn2jkmAmnH ,mnr sin jkmBmnr sin 2sin m kr Hkr Pm cos 1 mn mn 2 jkAmnH ,mnr jkBmnr2cosm kr Hkr Pm cos 1 mn mn 2d2sin m kr Hkr1 d mn 2d2cosm kr Hkr1 d mn 2 mcos P mn mcos P mnН r, mn 0Обозначим индексы для магнитных волн μ, тогда из граничных условий следует:d mP cos 0dЗдесь присоединённая функция Лежандра также выражается через функциюЛежандра P cos , удовлетворяющую условию:P 1 1269Поскольку μ не является целым числом, то функция P не будет полиномом,поэтому уравнениюd mP cos 0 удовлетворяет последовательность чиселdm1, m2...mn , зависящих от угла γ.Следует также отметить, что индексы ν и μ могут быть и целочисленными внекоторых частных случаях.