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

Lections_V

.pdf
Скачиваний:
11
Добавлен:
20.03.2015
Размер:
1.21 Mб
Скачать

Уравнение 2t3 + t(1− 2y) − x = 0 есть уравнение нормали к параболе y = x2 − b в точке (t, t2). Если заменить b и а на y и x, то перейдем к уравнению (6.20).

Мы приходим к тривиальному выводу, что центр тяжести (а, b) в положении равновесия должен располагаться непосредственно над точкой касания. Пересечем по- верхность (6.20) рис.6.10 плоскостью b = b0. Получим кубическую кривую:

2t3 t(1− 2b0 ) − a = 0 .

Если рассматривать а как функцию t, то эта функция имеет максимум и минимум точку перегиба при b0 = ½. То есть поверхность имеет складку.

Уравнение (6.20) может иметь 1,2 и 3 решения. Число решений зависит от координат (а, b). Каждое решение отвечает положению равновесия качалки. Наша поверхность

 

 

 

 

 

Z

X

 

 

 

1

2

 

Y

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

5

 

 

 

 

 

Z `

 

 

1

2

3

4

5

 

 

a

t

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

Р и с . 6 . 10 . С б о р к а У и тн и в п р о с тр а н с тв е и в п р о е к ц и и н а п ло с ко с ть .

катастроф делит пространство (t, a, b) на две области: в области выше поверхности, мы имеем U(t) > 0, в другой U(t) < 0. Минимум потенциальной энергии достигается при тех значениях t, где функция U(t) меняет знак с минуса на плюс.

Когда U(t) = 0 имеем одно решение это минимум; когда решений три одно из них максимум, а два других минимума. Два решения соответствуют: в одной из точек по- верхности (6.20) имеет вертикальную касательную плоскость. Это решение ни макси- мум ни минимум, т.е. V`(t) не меняет знака, а второе решение минимум. Внутренняя часть складки дает только максимум. Найдем условие, при котором U(t) = 0 имеет кратный корень. Это условие состоит в том, что U(t) =U`(t) = 0, т.е.:

2t3 + t(1− 2b) + a = 0

6t2 + (1− 2b) = 0

(6.22)

 

Исключим t из системы двух уравнений:

 

27a2 = 2(2b −1)3

(6.23)

Отметим, что значения а = 0, b = ½ играют особую роль. Они дают трехкратный ко- рень. В этом случае U(t) = t3 и t = 0 – корень кратности 3. Кривая (6.23) – полукубиче- ская парабола на (a, b) – плоскости, т.е. в пространстве управления, а множество точек

этой кривой называется бифуркационным множеством функции V. Эта кривая отделя- ет точки, дающие одно решение уравнения U(t) = 0 от точек, дающих три решения. Если смотреть на поверхность (6.20) сверху (проекция поверхности на плоскость аb), то видимый контур, вдоль которого перегибается поверхность даст прямую, аналогич-

ную (2.23) рис.6.10.

Посмотрим на это с другой точки зрения.

Мы уже говорили, что (7.22 а) есть уравнение нормали к параболе b = a2 в точке (t, t2). Когда t меняется, а а и b – фиксированы, получаем семейство нормалей. Огибающая этого семейства эволюта, находится из уравнения (6.22). Т.е. уравнение эволюты есть (6.23).Для заданного положения центра тяжести (а, b) мы имеем уравнение (6.23), определяющее число положений устойчивого равновесия. Будем изменять координаты центра масс, и будем интересоваться не динамикой качалки, а состоянием устойчивого равновесия рис.6.11. В точке 1 – только одно положение устойчивого равновесия. В точке 2 и 3 имеется два положения равновесия. Пока положение равновесия медленно меняет положение параболы. Однако в точке 4 уже не возможно непрерывное измене- ние положения параболы. Малый сдвиг точки вызывает сильное изменение положений равновесия. Происходит катастрофа.

1

2

3

4

5

Р и с .6 .1 1 .

Вопрос №61

Геометрическая каустика точечного источника, помещенного в прозрачную среду. Понятие структурной устойчивости каустик.

Как уже было сказано, тонкая линза собирает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

все лучи в одной точке фокусе. Однако по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ложение этой точки неустойчиво относи-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тельно слабых смещений пучка лучей отно-

Рис. 6.12 . Изм енение аперту ры пу чка.

сительно оси. То есть поверхность лучей

 

 

 

 

 

 

 

 

фокальная коническая поверхность

 

 

 

 

 

 

 

 

 

неустойчива. Какая поверхность яв-

 

 

 

 

Z

 

 

 

 

ляется устойчивой? Эту задачу решал

 

 

 

 

 

 

 

P

М. Берри. Им была рассмотрена ус-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тойчивость относительно расширения

n 0 = 1

 

 

 

 

 

r

 

 

 

 

 

апертуры пучка (рис. 6.12).

 

 

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим некоторый светящийся

n

 

 

 

 

 

 

X

 

 

 

 

i

 

 

 

объект (рис. 6.13), расположенный

 

D

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

под водой. Пусть это будет точечный

 

 

 

 

S

 

 

 

 

источник света S, помещённый на

 

 

 

 

 

 

 

 

глубине D ниже горизонтальной по-

Рис . 6. 13 . Наб лю дени е све тя щ е го ся

верхности Z = 0, имеющий коэффи-

 

 

 

о б ъекта в ср е де с п о каза телем

циент преломления n. Рассмотрим

 

 

 

пр ело м лени я n .

двумерную задачу. В этом случае источник S порождает одно параметрическое семей-

ство лучей. Удобно в качестве параметра выбрать угол i. Наблюдатель находящийся в

точке Р увидит, что светящийся объект располагается на части виртуальной каустики,

которая является касательной к лучу, исходящему из точки Р. Очевидно, для того что-

бы найти каустику, сначала надо найти уравнение луча. Его легко определить из прин-

ципа Ферма.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, траектория, проходящая через точку Р и параметризуемая углом i есть оптиче-

ский путь Ц(i, x, z) : Φ(i, x, z) = nD seci + (z2 + (x - Dtg i)2 )1/ 2

(6.16)

уравнение Φ (i , x , z )

 

 

= const – уравнение луча. По условию огибающей:

 

∂Φ = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.17)

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

Φ(i, x, z)

 

 

 

D2

 

 

ì

 

 

 

 

 

(x - D tg i) 2

 

ü

(6.18)

 

=

 

i

ínsin i -

(z

2

 

1/ 2

ý = 0

 

 

i

 

 

 

 

 

cos

 

î

 

 

 

+ (x - Dtg i)

)

þ

 

Продифференцируем (6.18):

2Φ

 

 

+

¶ ¶Φ

dx = 0

,

(6.19)

i

2 di

x i

 

 

2Φ

 

 

 

 

 

Φ

 

 

 

 

 

 

 

 

 

 

 

так как:

=

 

0,

 

= 0

 

 

 

 

 

 

 

 

(6.20)

i

2

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует:

 

ncos3 i + D cos r

(1− n2 sin2 i) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

Определим координаты x и y из уравнений (6.18) и (6.20):

 

 

ì

= D(n

2

-1)tg

3

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

= -

D(1- n

2

sin

2

i)

3 / 2

 

 

 

 

 

 

 

 

(6.21)

 

ïz

 

 

 

 

 

 

 

 

 

 

 

 

 

î

 

 

 

 

 

 

ncos3 i

 

 

 

 

 

 

 

 

 

 

 

 

График уравнения (6.21) в двумерном случае (рис. 6.14) представляет собой casp -

(клюв).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

3

n

æ

+

D

ö

ü3 / 2

D

 

 

 

 

 

 

 

 

 

 

x = ±í

2

ç z

 

n

 

÷

ý

 

 

n2 -1

 

 

 

 

 

 

 

 

 

î

 

 

è

 

 

 

ø

þ

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эволюта это линия, на которой рас-

 

 

 

 

O

r

положен источник S, если наблюдать

 

 

 

 

 

X

за ним под разными углами r = (0;900).

 

 

-D/n

Малые шатания источника S не изме-

 

 

 

 

няет форму каустики. Поэтому в дву-

 

 

 

 

i

 

мерном случае эта каустика устойчива.

 

 

 

-D

 

Если n → 1, то все лучи пересекутся в

 

 

 

 

 

 

 

 

 

одной точке (в источнике).

 

 

 

 

 

 

 

 

 

S

 

 

Рассмотрим теперь трёхмерный

Рис. 6.14. Положение изображения источника S

случай. На первый взгляд кажется, что

в 3 мерном пространстве мы будем

 

 

 

 

 

накаустической поверхности.

иметь в конце каустики клюв, полу-

 

 

 

 

 

 

 

Рис. 6.15. Две ветви каустики – casp и фокальная линия.

ченный вращением клюва вокруг оси z. Посмотрим, так ли это? Забегая вперёд, скажу, что это предположение ошибочное!

Пусть теперь уравнение луча, приходящего в точку Р зависит от третьей коорди- наты y, и пусть направление распространения луча определяется ещё и азимутальным углом β в дополнении к полярному углу i.

Тогда вместо уравнения (6.16) получаем:

Φ(i, β, x, y, z) = nD seci +{z2 + (x Dtg i cos β )2 + (y D tg isin β )2}1/ 2

(6.23)

Очевидно, что теперь условия для огибающей поверхности будет дополнено ещё

одной производной:

∂Φ

= 0 .

(6.24)

 

 

∂β

 

В удобной форме из (6.23) находим:

 

ìx = (Dtg i + z tg r) cos β

(6.25)

í

 

 

îy = (Dtg i + z tg r)sin β

 

Можно доказать, что в общем случае условие для огибающей поверхности есть требо-

 

 

 

 

 

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

вание равенства нулю якобиана:

i

∂ β

 

= 0

(6.26)

x

y

 

 

 

 

 

 

 

 

 

 

i

∂ β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

x2 + y2

 

 

(D tg i + z tg r) = 0 .

 

 

 

(6.27)

 

sin i

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

Второй множитель в уравнении (6.27) задает уравнение поверхности вращения клюва,

а первый множитель x2 + y2 = 0 при х = у

= 0 и z < 0 – линию.

То есть каустика получает вторую ветвь, содержащую фокальную линию, прости-

рающуюся с вершины клюва к поверхности воды рис. 6.15. На этой линии не определе- но множество пересечений лучей. Напри- мер, внутри клюва в двухмерном случае пе- ресекалось три луча. Вне клюва два. Эта неопределённость связана с неопределён-

ностью азимутальной координаты β на этой линии. Поэтому на линии пересекается не- определённо много лучей с любой величиной β. Отсюда точка фокуса в трёхмерном случае оказывается более неустойчивой, чем в 2 двухмерном случае.

Вопрос № 23

Отражение и преломление света на границе двух диэлектрическихсред. Граничные условия.

Любое дифференциальное уравнение в част-

ных производных имеет бесконечное множество возможных решений. Чтобы из всех решений вы- брать единственное требуется задачу дополнить на- чальными и граничными условиями.

Рассмотрим формулировку граничных усло- вий для оптических задач. Найдем связь между на-

r r r r

пряжённостями E, H и индукциями D, H электри-

ческого и магнитного полей при прохождении света через границу раздела двух диэлектриков.

Определим граничные условия для нормаль- ных компонент поля.

Пусть имеем поверхность раздела (рис. 7.1)

r

двух сред, где n12 единичный вектор к поверх-

ности границы раздела. Заменим поверхность рез- кого раздела тонким переходным слоем, внутри ко- торого ε и μ быстро, но непрерывно изменяется.

Пусть толщина слоя δh. Внутри слоя построим не- большую цилиндрическую поверхность. Основани- ем цилиндра являются площадки дA1 и дA2 , парал- лельные поверхности. Запишем:

 

δA2

 

 

 

 

 

 

 

n2

 

 

 

 

 

 

 

 

n12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δh

 

 

 

 

 

 

 

 

 

 

 

 

ε2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

δA1

 

 

n1

 

 

 

 

 

ε1 1

 

 

 

 

Рис. 7.1.Квыводу граничныхусловий для нормальныхкомпонентВи D.

 

divB = 0 .

 

r

(7.1)

 

 

Поскольку во всём цилиндре вектор индукции B и его производные непрерывны, применим теорему Га-

усса к интегралу по объёму:

 

r

 

 

 

 

r

= 0 Þ

ò div B dV

= 0

 

 

òdivBdV

= ò B * n dS

(7.2)

 

V

 

 

V

S

 

 

Поскольку площадки δА1 и δА2 малы, будем считать, что индукция B на них принимает значения B1 и

B2 . Уравнение (7.2) примет вид:

 

 

 

 

 

 

r

r

 

 

 

 

 

B1n1δA1 + B2n2δA2 + вклад боковых стенок = 0

(7.3)

 

Устремим к нулю высоту цилиндра дh → 0 . Тогда:

 

 

 

 

r

r

 

 

 

 

 

(B1n1 + B2n2 A = 0 ,

 

 

 

(7.4)

 

r

 

 

 

 

r

для второй среды, то

если n1 внешний вектор нормали к площадке δА1

для первой среды, а n2

r

r

r

 

 

 

 

 

n1

= −n2

= n12 - нормали направлены в разные стороны.

 

 

 

 

 

r

= 0

 

 

 

Откуда:

 

n12 (B2 − B1 )

(7.5)

 

 

Запись (7.5) означает, что на границе раздела нормальная компонента вектора индукции магнитного поля

непрерывна. Подобным образом можно поступить при нахождении соотношения для D1

и D2 . В сис-

r

 

теме СГС уравнения запишутся: òdivDdV = ò D *ndS = 4π ò ρdV

(7.6)

V

При слиянии площадок δА1 и δА2 полный заряд остаётся

 

 

 

 

 

Q2

t

 

n12

 

 

 

 

 

 

 

P2

 

 

 

 

 

 

 

 

 

 

 

 

конечным а, следовательно, объёмная плотность стре-

ε2, μ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мится к бесконечности. Поэтому вместо объёмной плот-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q1

 

b

)

 

 

 

 

 

 

ности ρ следует ввести поверхностную плотность с из

ε1, μ1

P1

 

 

 

соотношения:

lim ò

)

 

 

с dV = ò с dA ,

Рис. 7.2.Квыводу граничныхусловийдля

дh → 0

(7.7)

 

тангенциальныхкомпонентЕиН.

таким образом для малой площадки δА имеем:

)

 

r

r

 

D1n1δA1 + D2 n2δA2 + вклад боковых

стенок = 4πρδА

(7.8)

Тогда при δh → 0 следует:

 

 

r

)

 

 

n12 (D2 D1) = 4πρ . (7.9)

 

 

)

При наличии на поверхности раздела поверхностного заряда с нормальная компонента вектора электри-

)

ческой индукции испытывает скачёк равный 4рс .

Найдём граничные условия для тангенциальных компонент.

)r

Заменим поверхность раздела переходным слоем прямоугольной площади, где b – единичный вектор к плоскости прямоугольника P1Q1P2Q2. Рассмотрим одно из уравнений Максвелла:

r

= − 1

B

 

 

 

 

 

 

 

 

 

rot E

 

 

 

 

 

(7.10)

 

 

 

c

t

 

 

 

 

 

 

 

 

 

Возьмём от (7.10) поверхностный интеграл:

 

 

 

 

 

 

)

 

 

 

)

 

 

 

 

 

 

r r

 

1

 

r r

 

 

 

 

 

ò rot E b dS

= −

ò

B b dS

 

 

 

 

(7.11)

c

t

 

 

 

 

S

 

 

S

 

 

 

 

 

 

Применим к (7.11) теорему Стокса:

 

 

 

 

 

)

 

 

)

 

 

 

r

r

 

 

 

 

 

r r

 

 

 

1

 

r r

 

ò rot

E b

dS

=

ò

E

d l

ò

&

 

 

c

B b dS ,

(7.12)

S

 

 

 

l

 

 

 

 

 

 

где точка над буквой означает дифференцирование по времени. Пусть P1Q1 = dS1, а P2Q2 = dS2 очень малы. Тогда на этих длинах E1 и E2 не меняются. Пусть ∂h стремится к нулю, тогда S → 0, и:

r )

r )

 

1

r r)

r

r

+ вклад от стенок = -

&

E1 t1

+ E2 t2

с

Bb S = 0

или: (E1 E2 )tr = 0 ,

 

 

(7.13)

тангенциальная компонента напряжённости электрического поля непрерывна на границе раздела диэлек- триков.

Аналогично доказывается, что:

r

H2 ) = 0 .

 

n(H1

(7.14)

Нормальная компонента магнитного поля при отсутствии токов на границе раздела непрерывна.

 

 

 

 

 

 

 

 

 

Вопрос № 24

 

 

 

 

 

 

Формулы Френеля.Законы отражения и преломления.

 

Пусть плоская электромагнитная волна падает на плоскую границу раздела двух диэлектриков:

 

ì

r

r

 

 

 

rr

 

 

 

 

 

 

= E0 exp{it - kr )}

 

 

 

 

 

ï E

 

 

 

 

 

í

r

r

 

 

 

rr

 

(8.1)

 

 

 

ï H

= H

0

exp{it - kr )}

 

 

 

 

 

î

 

 

 

 

 

 

 

 

 

Предположим, что μ = 1 и n =

εμ =

ε . Согласно граничным условиям:

 

 

 

ìE

= E

 

 

 

 

 

 

 

 

 

í

 

 

 

 

 

(8.2)

 

 

 

îH= H

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

 

Hi, Ei

Hr, Er

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

ϕ′

 

 

 

 

 

 

 

 

 

 

 

n 1

 

X

 

 

ε1

 

 

 

 

 

 

X

n 2

Ht, Et

 

 

 

 

 

 

 

 

 

 

 

 

 

ε2

 

 

 

 

Y

 

 

 

 

ψ

r

 

 

 

 

 

 

Z

 

 

kr

k x

 

 

 

 

 

 

 

 

 

Z

 

 

а

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

k z

 

 

Рис. 8.1. Прелом ление и отражение п лоской волны : (а) - в пространстве;

 

 

 

(б) – п лоскость па дения.

 

 

 

Ei

= E0i

 

 

 

 

 

r

 

 

 

 

exp{iit k i r )}

, - падающая

 

 

 

E r

= E 0r

 

 

 

 

 

r

 

 

 

 

exp{ ir t k r r )} , - отражённая

(8.3)

 

 

E t

= E 0t exp{

i t t

r

 

 

 

 

k t r )} , - преломлённая.

 

 

 

В соответствии с граничными условиями (8.2) для тангенциальной компоненты:

 

 

Ei

 

 

 

 

 

 

r

 

r

r

}(8.4)

 

exp{iit k ir )}+ Er

exp{irt k rr )}= Et

exp{itt k tr

 

 

 

 

 

 

 

 

 

 

 

 

В линейных оптических процессах на границе раздела всегда выполняются условия:

 

ï

ω

 

 

= ω

 

= ω

 

ì

 

 

i

 

r

 

t

í k xi

= k xr = k xt

 

ï

 

i

 

r

 

t

(8.5)

î k y

= k y

 

= k y

 

Пусть волновой вектор падающей волны лежит в плоскости см. рис. 8.1, тогда kiy = 0 . Это

условие выполняется, если пренебречь потоком через боковые поверхности (оптический эффект Маг- нуса), тогда:

kyr = kyt = 0 .

(8.6)

Постулируется, что луч падающий, отражённый и преломлённый лежат в одной плоскости с нормалью в плоскости падения (рис. 8.1, б). Это условие выполняется только приблизительно на пло- ских границах раздела двух сред и вообще не выполняется на кривых границах!

Разложим волновой вектор на составляющие:

i

= k

i

sinφ,

 

r

= k

r

 

 

 

t

= k

t

sinψ .

 

kx

 

kx

 

sinφ ,

kx

 

(8.7)

 

 

 

 

r

 

ω , откуда волновые векторы падающей, отраженной и преломленной

Но из (2.8) следует, что | k |=

волн:

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ki = ω

, kr =

ω

, kt

=

ω

,

 

 

 

 

(8.8)

 

v

 

 

v

 

 

 

 

v

2

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

где v1 и v2 скорости волны в средах с абсолютным показателем преломления n1 и n2. Тогда из (8.5) и (8.7) следует:

sinφ = sinφ′ = sinψ

 

v1

 

v1

 

 

.

(8.9)

 

 

v2

 

Выражение (8.9) характеризует два закона Снеллиуса: закон отражения (угол падения равен углу от- ражения)

ϕ = ϕ′, (8.10)

и закон преломления:

sinφ = v1

=

ε2

= n2

 

 

 

 

 

 

 

(8.11)

sinψ v2

 

ε1

 

 

 

n1

Соотношения между амплитудами падающей, отражённой и прошедшей волн. Формулы Френеля.

Отобразим компоненты электрического вектора для падающей, отраженной и преломленной волн на рис. 8.2. Разложим вектор напряжённости поля на две компоненты: параллельную и перпендикуляр-

ную составляющие: E i = E pi + ESi , E r = E pr + ESr , E t = E pt

+ ESt

,

(8.12)

 

 

Тогда компоненты поля из (8.12) можно представить в следующем виде:

 

 

Exi = Eip cosφ, Eiy = ESi , Eiz = −Epi sinφ ,

(8.13 а)

 

 

 

 

 

 

 

 

k zr

r

 

E iz

E ip

 

 

 

 

r

E zr

k r

 

 

 

 

 

 

E p

 

r

 

ϕ E i

 

 

 

E r

 

 

k x

 

 

 

 

 

E Sr = E ry

 

E Si = E iy

x

i

 

 

x

 

 

 

 

 

k x

 

ϕ

 

 

 

 

 

 

 

r

ϕ

 

 

 

 

 

 

i

 

 

 

X

 

 

 

 

 

 

 

 

 

 

k i

 

 

 

 

 

 

 

k z

 

E zt

 

 

 

 

 

 

 

 

ψ

E pt

 

 

 

 

 

 

E St

= E ty

E xt

 

 

 

 

 

 

 

 

 

 

k xt

 

 

 

 

 

Z

 

 

 

r

 

 

 

 

 

 

 

t

k t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k z

 

 

 

Рис. 8.2. Распо ло ж ение ко м по нент эле ктр ическо го век то р а Е в

 

падаю щ ей, о тр аженно й и пр ело м ленно й во лна х.

 

Exr

= −Epr cosφ, Ery = ESr , Ezr

= −Epr

sinφ ,

 

(8.13 б)

Ext

= Etp cosψ , Ety

= ESt , Ezt

= −Ept

sinψ .

 

(8.13 в)

Согласно граничным условиям:

 

 

 

 

 

 

 

Eix

+ Exr = Ext

, Eiy

+ Ery = Eiy ,

 

 

 

(8.14 а)

Hix

+ Hrx = H xt , Hiy

+ Hry

= Hty .

 

 

 

(8.14 б)

В соответствии с (8.13) и (8.14 а) получим систему уравнений:

(Eip Epr )cosφ = Etp cosψ ,

(8.15)

ESi + ESr = ESt

(8.16)

Согласно уравнению: μ H = ε ( k × E) с учетом того, что для прозрачной среды μ = 1, представим вектор магнитного поля через компоненты электрического поля:

 

 

 

r

r

r

 

 

 

 

 

 

 

r

 

 

i

j

k

 

 

 

 

kx

ky

kz

 

 

 

 

 

 

H = ε

,

(8.17)

 

 

 

Ez

Ey

Ez

 

 

 

 

 

 

 

 

 

 

 

 

 

r

: k ix

= k i sin

φ ;

k iy = 0;

k iz

= k i cos φ ,

k i

)

 

 

 

 

 

 

 

krr

: k rx

= k r sin

φ ;

k ry

= 0 ;

k zr

= − k r cos φ ,

где )

 

 

 

 

 

 

(8.18)

krt

: k xt

= k t sin ψ ;

k ty

= 0;

k zt

= k t cos ψ .

Для падающей волны вектор магнитного поля можно представить в виде:

 

 

 

 

r

r

r

 

 

 

 

r

 

 

 

i

j

k

= ε i

 

sin φ

 

cos φ

H i

 

0

 

 

 

 

E xi

E yi

E zi

Разложим вектор Нi по компонентам поля:

Hxi = −ε1Eyi cosφ, Hiy = ε1 (Exi cosφ − Ezi sinφ), Hiz = ε1Eiy sinφ .(8.19 а)

Соответственно, для отраженной волны вектор магнитного поля:

 

 

 

r

r

r

 

 

 

r

 

 

i

j

k

 

 

sinφ

 

cosφ

Hr = ε1

 

0

 

 

 

Exr

Eyr

Ezr

 

 

 

 

 

 

Разложим по компонентам магнитного поля:

Hxr =ε1Eyr cosφ, Hry =ε1(Exr cosφ −Ezr sinφ), Hrz =ε1Eyr sinφ. (8.19 б)

И для преломленной волны:

 

 

 

 

 

r

r

 

r

 

 

 

 

 

 

 

 

 

 

r

 

 

i

j

 

k

 

 

 

 

Ht = ε2

 

sinψ

0

cosψ

 

 

 

 

 

 

 

 

 

Ext

Eyt

 

Ezt

 

 

 

 

Представим по компонентам поля:

 

 

 

 

 

 

Hxt =−

 

Eyt cosψ, Hyt

=

 

(Ext cosψ −Ezt sinψ), Hzt =

 

Eyr sinψ . (8.19в)

ε2

ε2

ε2

Подставим компоненты поля (8.19) в уравнения (8.14 б) для компонент x поля, получим:

ε1 Eyi cosφ + ε1 Eyr cosφ = −ε2 Eyt cosψ .

Исходя из (8.13), приведем вышестоящее выражение к виду:

ε1

(ESi ESr ) cosφ =

ε2

ESt cosψ

(8.20 а)

Аналогично поступим для y компонент поля:

ε1 (Exi cosφ − Ezi sinφ) + ε1 (Exr cosφ − Ezr sinφ) = ε2 (Ext cosψ − Ezt sinψ ) .

Преобразуем это выражение, используя (8.13):

ε1 (Eip cos2 φ + Eip sin2 φ) + ε1 (Epr cos2 φ + Epr sin2 φ) = ε2 (Etp cos2ψ + Etp sin2ψ )

Упростим вышестоящее выражение:

ε1

(Epi + Epr ) =

ε2

Ept .

(8.20 б)

Подставим значение ESt в уравнение (8.20 а), получим:

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]