Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Маркузе, Д. Оптические волноводы

.pdf
Скачиваний:
45
Добавлен:
22.10.2023
Размер:
17.73 Mб
Скачать

450

Глава I)

искажением стенки, покидает волновод и исчезает в про­ странстве. Ниже в этом разделе будет рассчитана структу­ ра поля излучения в дальней зоне. При потере малой части мощности падающая мода практически ие изменяется и ее амплитуда уменьшается незначительно. Преобразованная в излучение мощность больше не направляется волново­ дом, а излучается в пространство. Это не имеет места для случая связи падающей моды с направляемой, поскольку переданная направляемой моде мощность снова взаимодей­ ствует с мощностью падающей моды. Однако в случае излучения такое взаимодействие очень слабое. Этот результат позволяет расширить диапазон применения теории возмущений. Используя формулу (9.3.19) неодно­ кратно к уменьшенному значению мощности падающей моды, получим на конце волновода длиной D мощность

p = p 0e - 2 a D '

(9.3.22)

Коэффициент затухания а определяется из формулы (9.3.19). Выражение (9.3.22) справедливо для малых значений 2аХ, где X — длина волны, и точность его увеличивается с ро­ стом угла, под которым излучение покидает волновод.

Из анализа структуры поля в дальней зоне получаем угол излучения [100]. Для коэффициента разложения в предположении (9.3.3) из формул (9.2.38), (9.3.5) п (9.3.0)

находим следующее выражение:

9е+)(Р. L ) — — («1 —

п \ ) X

• т, L

 

арк% cos я0Ф) cos аФ)еШЬ/ 2

X

81,1 Г т

'[/"яр (р 0Ф)

(р2 cos2 orf0 + а2si n2 ad0)

, (9.3.23)

 

 

где Г определяется формулой (9.3.16). Составляющая электрического поля излучения получается из (8.5.16) в виде

со

 

£ „ = j ?<+’ (Р, Ь )£ г/(р, z)dp.

(9.3.24)

о

 

Выражения (9.3.23) и (9.3.24) описывают поле излучения плоского волновода с синусоидальным изменением его толщины в точке пространства с координатами х, z. Излу-

Лерегцлярп ые Ни.ментр теснив вол поводы

451

чеппе вызывается нерегулярным участком длиной L пло­ ского волновода. Нас интересует только поле в дальней зоне, поэтому можно считать, что

L < z.

(9.3.25)

Поле моды излучения Е у (р) вне слоя волновода опреде­ ляется формулами (8.4.4), (8.4.9), (8.4.10) и (8.4.18).

Объединение этих выражений приводит к следующему результату:

р V 2сор0Р ^cos ad0 cos р ( | х | —dQ)

Я„(Р) =

— ^-sinadosnip d- T — <4)1 )J

V

.-iflz. (9.3.2G)

 

яр (р2 cos2 od0-j- a2 sin2 ad0)

Окончательно выражение для составляющей электриче­ ского поля излучения получается из формул (9.3.23), (9.3.24) и (9.3.2(1) в виде

гг_ оА-д У 2 ш р 0Р ('1Г — « !) cosy-orfo v,

 

 

-

-=

х

 

 

,rL'2p2 cos adо |^соз ad0 cos p (| x |d0)—^ sin ad0sin p (| x | d0) j

 

X

J0

 

 

P (p2 cos2 arf0-|-o2sin2 ad0)

X

 

 

 

 

sin Г A

 

 

 

(9.3.27)

 

 

X — гг—

dp.

 

 

 

 

Г

 

 

 

 

 

Интеграл в этом выражении достаточно сложен и не при­ ходится рассчитывать на получение общего решения. Одна­ ко можно получить его приближенное значение методом стационарной фазы [9] и определить, таким образом, поле в дальней зоне. Этот полезный приближенный метод вычис­ ления интегралов, имеющих быстро осциллирующие функ­ ции в подынтегральных выражениях, был рассмотрен в разд. 2.3. Предполагается, что х и z становятся больши­ ми по сравнению с L. Это заставляет функции, содержа­ щие х и z в аргументах, осциллировать намного быстрее функций с L и dQв аргументах. Главный вклад в интеграл

20*

452

Глава i)

имеем в области, где аргументы быстро осциллирующих функций принимают стационарные значения. Таким обра­ зом, нужно проанализировать выражения вида

 

Р* +

РZ

= / (р).

(9.3.28)

Используя

формулу (8.4.6),

получим

 

 

df

 

р

(9.3.29)

 

- + = £ ----п г

 

dp

 

р

 

и

ri2/

 

 

 

 

 

пЩ ,

(9.3.30)

 

dpi —

 

р“з Z-

Требование

стационарности

фазы dj/dp =

0 приводит

к условию

 

 

 

 

 

J js - = 4 = tg a ..

(9.3.31)

 

Ps

*

 

 

Разлагая / (р) в ряд Тейлора в окрестности стационарной точки, получаем

СО

j e -i(px+PO ф ~ g-*(Pe*+ P sz) х

О

оо

X j e i(nnho/2Ps)-(p-os)2 ^ (р — р 3) = "

— СО

= (1 + 0 J +

e- i(V+M). (9.3.32)

П 2^о

^

Для оценки интеграла была использована формула (2.3.23). В (9.3.27) необходимо было вычислить лишь интеграл (9.3.32). Разложим sin и cos по экспоненциальным функ­ циям. В разложении необходимо рассмотреть только ту его часть, которая содержит функции вида (9.3.28). Экспоненциальные функции с аргументами (pz — рх) дают незначительные вклады при х я z больше нуля. Множи­ тели при экспоненциальных функциях можно рас­ сматривать как постоянные в точке р = ps и их можно вынести из-под интеграла. Окоичательное интегрирование

Нерегулярные диэлектрические

волноводы

453

выполнено в (9.3.32). Таким

 

образом,

имеем

 

(1 г) k0Ps У Ps У 2шрqP {ill— nl) X

 

 

X cos Kpdocos as d0etPsCioe,rsL/2]

 

^

 

[ 23 У n ]/~ ^ Po^o ~r

j X

 

 

 

X (pscos asd0-\-iossin asdQ)j

 

 

 

■ г

L

-inzkoix sin a s+zcos a s)

 

sm I s

-s-

(9.3.33)

X ----- -----------------------

 

 

-y =

------------- ,

1 s

 

 

у

z

 

 

где

n2/c0since*

 

 

(9.3.34)

ps=

 

 

и

n2/cocosa*,

 

 

(9.3.35)

P* =

 

 

а Г* получается из (9.3.16) заменой p на (3*. Выражение для составляющей поля в дальней зопе указывает па не­ сколько интересных особенностей. Видно, что поле рас­ пространяется примерно так же, как плоская волна под углом as, под которым поле приходит в точку х, z из источника в области нерегулярного участка х). Для поля в дальней зоне источник ограниченной протяженно­ сти подобен точечному источнику, расположенному в на­ чале координат. Последний множитель выражения (9.3.33) представляет собой цилиндрическую волну. Это становится более очевидным, если ввести в соответствии с (9.3.31)

и

a:=

?-sinas

(9.3.36)

z=

rcosa*.

(9.3.37)

 

Тогда получается соотношение

 

г

-—inzhotx sin Gts+z cos a s)

„____ „-inzhor

/ P

s ----------zrp---------- = V n 2k0e— — . (9.3.38)

 

Vz

 

V r

Появление здесь вместо сферической волны цилиндриче­ ской связано с двумерной геометрией волновода. Наиболее

*) Из этого примера видно, что суперпозпцпя стоячих воли излучения в интегральном разложении дает распространяющуюся волну.

454 Глава 9

важным множителем в формуле (9.3.33) является

• г

L

si» Г8

-5-

G=

(9.3.39)

*S

Известно, что on является приближением дельта-функцпи для возрастающих значений L. Чтобы величина G/л была хорошей аппроксимацией дельта-функции, необходимо, чтобы L было много больше длины волны К. Это условие не противоречит условию (9.3.25), при котором было полу­ чено поле в дальней зоне. Множитель (9.3.39) определяет диаграмму направленности излучения участка плоского вол­ новода с синусоидальным возмущением стенки. Главный ле­

песток атой

диаграммы становится уже с увеличением L.

Направление

главного

лепестка в пространство задается

формулой

Г, =

0 —(р0_ р ’,) = 0.

(9.3.40)

Поэтому основное излучение происходит только в направ­ лении, определяемом из формулы (9.3.35):

соза, = -^ = ^ -.

(9.3.41)

"2*0

 

Излучение направлено перпендикулярно поверхности вол­ новода при 0 = ро п по касательной к его поверхности при 0 = Ро — »2^о- Другими словами, излучение перпендику­ лярно поверхности при ps = Он приближается к касатель­ ному с увеличением |3S. При ps ~ п2к0 излучение главного лепестка параллельно поверхности волновода. Большие значения постоянных распространения мод излучения не­ возможны. При ps > п2к0 получаем направляемые моды. Данное рассмотрение достаточно ясно показывает, каким образом поле излучения становится все более направлен­ ным вперед и как его направление достигает своего пре­ дела — направления, параллельного оси волновода, кото­ рое сливается с направляемой модой. При рассмотрении направления излучения лепестка диаграммы направлен­ ности и изменения этого направления предполагалось, что можно непрерывно перестраивать механическую частоту 0. Точность приближения в формуле (9.3.22) возрастает при увеличении угла a s, определенного в (9.3.41).

Читатель может удивиться, что обсуждению плоского волновода с синусоидальными искажениями стенки отве-

Нерегулярные диэлектрические нолиоводы.

455

депо так много места. К тому же совершенно невероятно, чтобы такой случай мог иметь место на практике, если не позаботиться специально о создании волновода с подобной геометрией. Но все недоразумения исчезают, если функцию, описывающую любое произвольное искажение стенки, раз­ ложить в ряд или интеграл Фурье по синусоидальным функ­ циям. Мы видели, что данное синусоидальное искажение стенки связывает падающую направляемую моду только с одной направляемой модой или только с очень узкой областью непрерывного спектра мод излучения. Таким образом, наше рассмотрение применимо к каждой фурьекомпоиенте разложения произвольной функции искаже­ ния стенки. При определении потерь на излучение для произвольной деформации стенки введем фурье-компоненту F2dQ вместо амплитуды а2 и будем интегрировать (9.3.19)

по интервалу (9.3.12).

При этом выражения

(9.2.38)

и (9.2.39) подставляются

в формулу (9.3.14). При

Г =

0

из (9.3.16) получим йр =

dQ. Величина | F |2

[ II

|2)

называется спектром мощности функции искажения степкн. Этот термин заимствован из анализа Фурье электри­ ческих сигналов. Квадрат абсолютного значения коэффи­ циента разложения Фурье изменяющегося во времени сигнала пропорционален плотности мощности при данной

частоте.

Ограничиваясь

для

простоты искажением стенки

только с одной стороны волновода, h (z ) = 0 ,

из формул

(9.2.38),

(9.2.39) и (9.3.14),

получим

 

 

 

нгЬо

 

 

 

 

1

АР = j

7 (p )|F (P o -P )|2Ldp,

(9.3.42)

 

L

Р

 

 

 

где

 

—П2&0

 

 

 

pA-;j cos2 x0d0

 

 

 

 

 

 

 

cos2 ad0

 

 

(9.3.43)

p2 cos2 ad0+

a2 sin2 ad0

 

 

 

 

 

Результаты, полученные при рассмотрении влияния строго синусоидального искажения стенки, помогают интуитивно предсказывать влияние более сложного искажения стенки,

456 Глава 9

которое всегда можно представить в виде интеграла Фурье. Если спектр функции искажения стенки известен, можно вычислить потери на излучение спомощью формулы (9.3.42). Обычно спектр мощности функции искажения стенки не из­ вестен, поэтохму необходимо изучить некоторые модели, чтобы попять и оценить ожидаемые потери. Вместо изу­ чения моделей с различными «спектрами мощности» мы в следующем разделе воспользуехмся статистическим под­ ходом.

С помощью теории возмущений было показано, что потери на излучение, испытываемые направляемой модой, можно вычислить даже для достаточно длинных участков волноводов. Теперь покажем, что аналогичный подход можно использовать для изучения взаимодействия между направляемыми модами. Нельзя, конечно, непосредственно распространить формулу для потерь (9.3.10) на (9.3.22). Потери по формуле (9.3.10) не зависят линейно от L, поэто­ му пз нее нельзя определить потери на единицу длины. Физические соображения, использованные для получения соотношения (9.3.22) пз (9.3.19), неприменимы в дан­ ном случае к направляемым модам, поскольку мощность, переходящая от падающей моды к паразитной моде, остает­ ся в волноводе и продолжает распространяться, взаимо­ действуя с падающей модой. Существуют определенные закономерности связи направляемых мод, вызванной стро­ го синусоидальным искажением стенки, которые можно использовать для получения решения, справедливого за пределами применимости теории возмущений. Выше мы видели [см. формулу (9.3.7)], что при синусоидальном иска­ жении стенки эффективно связаны только две направляе­ мые моды. Поэтому пренебрежем влиянием всех других направляемых мод п поля излучения и будем считать, что взаимодействуют только две направляемые моды, связан­ ные синусоидальным искажением стенки. Обозначая ампли­ тудные коэффициенты этих двух мод через с0и щ, из урав­ нения (9.2.6) получим

д2ер

2 ifo - ? g - = - c 0F00- c iF10

(9.3.44)

9z2

и

(9,3.45)

Нерегулярные диэлектрические волноводы

457

Эта система уравнений применима только в особом случае синусоидально изменяющейся функции искажения стенки. Здесь диагональные члены F00 и Fn изменяются синусои­ дально в зависимости от z. Они представляют собой сину­ соидально изменяющиеся поправки к постоянной распро­ странения для регулярного волновода, стенка которого имеет среднее положение. Так как рассматривается незна­ чительное отклонение стенки, то этими диагональными членами можно пренебречь. Предположим, что обе моды распространяются в положительном иаправленип осп z. В отсутствие связи коэффициенты с0и Cj являются постоян­ ными. Для слабой связи, которая рассматривается в дан­ ном случае, с0 и щ изменяются на длине волны очень незначительно. Поэтому можно пренебречь вторыми про­

изводными коэффициентов

с0 н с, по сравнению с много

большими произведениями

первых

производных на 2 р.

Пренебрегая

указанными

малыми членами,

получим

следующую

упрощенную

систему

уравнений

[сравните

с уравнениями связанных волн (2.7.15) и (2.7.16)]:

 

дс0

= —XCj

 

(9.3.46)

 

дг

 

 

 

 

 

 

 

(9.3.47)

Коэффициенты связи в (9.3.46) и (9.3.47) можно определить из формулы (9.2.7). Однако проще воспользоваться фор­ мулами (9.2.34), (9.3.8) и (9.3.9) при L = z. Учитывая, что коэффициент с0 равен единице в формуле (9.2.34), путем дифференцирования по z п сравнения с (9.3.47), получим

a k f i COS K Q d g COS X jd 0

. (9.3.48)

 

M o + ^r) ( M o + -?r

Члены типа cos xd0 можно исключить с помощью формулы (8.6.16). Коэффициент связи является вещественной вели­ чиной. Коэффициент связи для соотношения (9.3.46) полу­ чается заменой с0 на cLв (9.2.34). Как видно из формулы (9.3.48), коэффициент х инвариантен относительно пере­ становки индексов 0 и 1. Изменение в знаке в (9.3.46) можно получить, рассматривая коэффициенты Фурье

458

Глава 9

в формулах (9.2.35) и (9.3.8). При выводе равенства (9.3.8) мы пренебрегли членом вида

sm х

L

(9.3.49)

где х = 0 + (Ро — Pi), поскольку х конечная величина. Однако после перестановки индексов 0 н 1 этот отброшен­ ный член становится существенным, а членом, содержа­ щим х = 0 — (р0 — p j, можно пренебречь. Два остав­ ленных члена имеют противоположные знаки, поскольку они получены из разложения функции sin 0z на экспонен­ циальные компоненты. Данное рассмотрение объясняет различие знаков в формулах (9.3.46) н (9.3.47).

Исключая Ci из формул (9.3.46) и (9.3.47), получим

- ^ + х 2с0= 0

(9.3.50)

и аналогично

 

- g f + x 2Ci= 0.

(9.3.51)

Так как предполагалось, что с0 — амплитудный коэффи­ циент падающей моды, причем с0 = 1 при z — 0 , то имеем решение

с0 =

cos xz

(9.3.52)

и

sin xz.

(9.3.53)

С] =

Для малых значений xz из формулы (9.3.53) при z — L получаем

^ . = | Cl|* = x*La,

(9.3.54)

что хорошо согласуется с (9.3.10).

Решения (9.3.52) и (9.3.53) справедливы за пределами применимости теории возмущений. Они хорошо описывают взаимодействие мощностей двух направляемых мод, кото­ рые связаны посредством синусоидального искажения стенки. Теперь стало ясно, почему приближенные отно­ сительные потери мощности (9.3.10) пропорциональны L2, а ие L.

Две направляемые моды могут полностью обмениваться своей мощностью. Из формул (9.3.52) и (9.3.53) видно, чтр

Нерегулярные диэлектрические волноводы

459

мода с индексом 1 растет за счет моды с индексом 0. Поело прохождения расстояния

D =

(9.3.55)

падающая мода полностью передает свою мощность моде

с индексом 1 .

Теория взаимодействия двух мод, связанных посред­ ством синусоидального искажения стенки, не может быть распространена на случай произвольной функции иска­ жения стенки. Хотя только одна фурье-компонеита функ­ ции искажения стенки ответственна за связь двух мод, данное рассмотрение оказывается несправедливым, если синусоидальная функция (9.3.1) является одной из многих фурье-компонент полной функции искажения стенки. Существуют по крайней мере две причины, в силу которых простая теория связи неприменима в этом более сложном случае. Прежде всего отметим, что гораздо большее число мод связано с падающей направляемой модой. Моды теря­ ют мощность и на излучение за счет фурье-компонент, соответствующих непрерывному спектру, и из-за связи с другими направляемыми модами. Таким образом, в общем случае нельзя ограничиваться рассмотрением связи только двух направляемых мод, а необходимо принимать во вни­ мание и другие конкурирующие процессы. Однако, если даже потери па излучение были бы малы и кроме падаю­ щей существовала только одна направляемая мода, мы не смогли бы использовать простую теорию связи мод, поскольку направление передачи мощности зависит крити­ чески от фазового соотношения между двумя связанными модами и фазы синусоидальной функции изменения стенки. Фурье-компонеита является функцией z, так как она полу­ чается интегрированием в точке z, в которой рассматрива­ ются обе моды. Не только амплитуда а фурье-компоненты изменяет свое значение с увеличением z, но, что более важ­ но, и ее фаза не имеет постоянного значения, необходимого для непрерывного обмена энергией. Для случайной функ­ ции / (z) фаза фурье-компонеиты изменяется случайно как функция z. Это заставляет направление передачи мощ­ ности между двумя модами меняться случайным образом то р одну, то в другую сторону, так что полная передача

Соседние файлы в папке книги из ГПНТБ