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

Lections_V

.pdf
Скачиваний:
16
Добавлен:
20.03.2015
Размер:
1 Мб
Скачать
☆

Уравнение 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{i(ωt - kr )}

 

 

 

 

 

ï E

 

 

 

 

 

í

r

r

 

 

 

rr

 

(8.1)

 

 

 

ï H

= H

0

exp{i(ωt - kr )}

 

 

 

 

 

î

 

 

 

 

 

 

 

 

 

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

εμ =

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

 

 

 

ìE

= E

 

 

 

 

 

 

 

 

 

í

1τ

 

 

2τ

 

 

 

(8.2)

 

 

 

îH1τ = H2τ

 

 

 

 

 

 

 

 

 

 

 

 

 

Е

 

 

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{i(ω it − k i r )}

, - падающая

 

 

 

E r

= E 0r

 

 

 

 

 

r

 

 

 

 

exp{ i(ω r 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{i(ωit − k ir )}+ Er

exp{i(ωrt − k rr )}= Et

exp{i(ωtt − k tr

 

0τ

 

 

 

 

 

 

0τ

 

0τ

 

 

 

 

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

 

ï

ω

 

 

= ω

 

= ω

 

ì

 

 

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 а), получим:

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