Гамильтонова оптика

редактировать

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

Содержание
  • 1 Принцип Гамильтона
  • 2 Лагранжева оптика
    • 2.1 Принцип Ферма
    • 2.2 Эйлер- Уравнения Лагранжа
    • 2.3 Оптический момент
    • 2.4 Уравнения Гамильтона
  • 3 Применения
    • 3.1 Преломление и отражение
    • 3.2 Лучи и волновые фронты
    • 3.3 Фазовое пространство
    • 3.4 Сохранение внешнего пространства
    • 3.5 Оптика для формирования изображений и без формирования изображений
  • 4 Обобщения
    • 4.1 Общая параметризация лучей
    • 4.2 Обобщенные координаты
  • 5 См. Также
  • 6 Ссылки
Принцип Гамильтона

В физике, принцип Гамильтона гласит, что эволюция системы (q 1 (σ), ⋯, q N (σ)) {\ displaystyle \ left (q_ {1} \ left (\ sigma \ right), \ cdots, q_ {N} \ left (\ sigma \ right) \ right)}\ left (q_ {1} \ left (\ sigma \ right), \ cdots, q_ {N} \ left (\ sigma \ right) \ right) описывается N {\ displaystyle N}N обобщенными координатами между двумя заданными состояниями при двух заданных параметрах σ A и σ B - это стационарная точка (точка, где вариация равна нулю), action функционал или

δ S знак равно δ ∫ σ A σ BL (q 1, ⋯, q N, q ˙ 1, ⋯, q ˙ N, σ) d σ знак равно 0 {\ Displaystyle \ delta S = \ delta \ int _ {\ sigma _ {A}} ^ {\ sigma _ {B}} L \ left (q_ {1}, \ cdots, q_ {N}, {\ dot {q}} _ {1}, \ cdots, {\ dot { q}} _ {N}, \ sigma \ right) \, d \ sigma = 0}\ delta S = \ delta \ int _ {{\ sigma _ { {A}}}} ^ {{\ sigma _ {{B}}}} L \ left (q_ {1}, \ cdots, q_ {N}, {\ dot {q} } _ {1}, \ cdots, {\ dot {q}} _ {N}, \ sigma \ right) \, d \ sigma = 0

где q ˙ k = dqk / d σ {\ displaystyle {\ dot {q}} _ {k} = dq_ {k} / d \ sigma}{\ dot {q}} _ {k} = dq_ {k} / d \ sigma . Условие δ S = 0 {\ displaystyle \ delta S = 0}{\ displaystyle \ delta S = 0} действительно тогда и только тогда, когда выполняются уравнения Эйлера-Лагранжа

∂ L ∂ qk - dd σ ∂ L ∂ q ˙ К знак равно 0 {\ Displaystyle {\ frac {\ partial L} {\ partial q_ {k}}} - {\ frac {d} {d \ sigma}} {\ frac {\ partial L} {\ partial {\ точка {q}} _ {k}}} = 0}{\ frac {\ partial L} {\ partial q_ {k }}} - {\ frac {d} {d \ sigma}} {\ frac {\ partial L} {\ partial {\ dot q} _ {k}}} = 0

с k = 1, ⋯, N {\ displaystyle k = 1, \ cdots, N}k = 1, \ cdots, N .

Импульс определяется как

pk = ∂ L ∂ q ˙ k {\ displaystyle p_ {k} = {\ frac {\ partial L} {\ partial {\ dot {q}} _ {k}}}}p_ {k} = {\ frac {\ partial L} {\ partial {\ dot q} _ {k}}}

и метод Эйлера-Лагранжа уравнения затем можно переписать как

p ˙ k = ∂ L ∂ qk {\ displaystyle {\ dot {p}} _ {k} = {\ frac {\ partial L} {\ partial q_ {k}}}}{\ dot p} _ { k} = {\ frac {\ partial L} {\ partial q_ {k}}}

где p ˙ k = dpk / d σ {\ displaystyle {\ dot {p}} _ {k} = dp_ {k} / d \ sigma}{\ dot {p}} _ {k} = dp_ {k} / d \ sigma .

Другой подход к решению этой проблемы состоит при определении гамильтониана (взяв преобразование Лежандра из лагранжиана ) как

H = ∑ kq ˙ kpk - L {\ displaystyle H = \ sum _ {k} {{ \ dot {q}} _ {k}} p_ {k} -L}H = \ sum _ {k} {{\ dot q} _ {k}} p_ {k} -L

, для которого новый набор дифференциальных уравнений может можно получить, посмотрев, как полный дифференциал лагранжиана зависит от параметра σ, позиции qi {\ displaystyle q_ {i} \,}q_i \, и их производные q ˙ i {\ displaystyle {\ dot {q}} _ {i}}\ dot q_i относительно σ. Этот вывод такой же, как в гамильтоновой механике, только с заменой времени t на общий параметр σ. Эти дифференциальные уравнения представляют собой уравнения Гамильтона

∂ H ∂ q k = - p ˙ k, ∂ H ∂ p k = q ˙ k, ∂ H ∂ σ = - ∂ L ∂ σ. {\ displaystyle {\ frac {\ partial H} {\ partial q_ {k}}} = - {\ dot {p}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial p_ {k}}} = {\ dot {q}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial \ sigma}} = - {\ partial L \ over \ partial \ sigma} \,.}{\ frac {\ partial H} {\ partial q_ {k}}} = - {\ dot {p}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial p_ {k}}} = {\ dot {q}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial \ sigma}} = - {\ partial L \ over \ partial \ sigma} \,.

с k = 1, ⋯, N {\ displaystyle k = 1, \ cdots, N}k = 1, \ cdots, N . Уравнения Гамильтона - это дифференциальные уравнения первого порядка, а уравнения Эйлера-Лагранжа - второго порядка.

Лагранжева оптика

Общие результаты, представленные выше для принципа Гамильтона, могут быть применены к оптике. В 3D евклидовом пространстве обобщенные координаты теперь являются координатами евклидова пространства.

Принцип Ферма

Принцип Ферма гласит, что оптическая длина пути, пройденного при свете между двумя фиксированными точками, A и B, является неподвижной точкой. Это может быть максимум, минимум, константа или точка перегиба . В общем, когда свет распространяется, он движется в среде с переменным показателем преломления, который представляет собой скалярное поле положения в пространстве, то есть n = n (x 1, x 2, x 3) {\ displaystyle n = n \ left (x_ {1}, x_ {2}, x_ {3} \ right)}{\ displaystyle n = n \ left (x_ {1}, x_ {2}, x_ {3} \ right)} в 3D евклидовом пространстве. Предполагая, что теперь свет движется вдоль оси x 3, путь светового луча может быть параметризован как s = (x 1 (x 3), x 2 (x 3), x 3) {\ Displaystyle s = \ left (x_ {1} \ left (x_ {3} \ right), x_ {2} \ left (x_ {3} \ right), x_ {3} \ right)}{\ displaystyle s = \ left (x_ {1} \ left (x_ {3} \ right), x_ {2} \ left (x_ {3} \ right), x_ {3} \ right)} начиная с точки A = (x 1 (x 3 A), x 2 (x 3 A), x 3 A) {\ displaystyle \ mathbf {A} = \ left (x_ {1} \ left (x_ {3A} \ right), x_ {2} \ left (x_ {3A} \ right), x_ {3A} \ right)}{\ displaystyle \ mathbf {A} = \ left (x_ {1} \ left (x_ {3A} \ right), x_ {2} \ left (x_ {3A} \ right), x_ {3A} \ right)} и заканчивается в точке B = (x 1 (Икс 3 В), Икс 2 (Икс 3 В), Икс 3 В) {\ Displaystyle \ mathbf {B} = \ left (x_ {1} \ left (x_ {3B} \ right), x_ {2} \ left (x_ {3B} \ right), x_ {3B} \ right)}{\ displaystyle \ mathbf {B} = \ left (x_ {1} \ left (x_ {3B} \ right), x_ {2} \ left ( x_ {3B} \ right), x_ {3B} \ right)} . В этом случае, по сравнению с принципом Гамильтона выше, координаты x 1 {\ displaystyle x_ {1}}x_ {1} и x 2 {\ displaystyle x_ {2} }x_{2}взять на себя роль обобщенных координат qk {\ displaystyle q_ {k}}q_ {k} , а x 3 {\ displaystyle x_ {3}}x_ {3} играет роль параметра σ {\ displaystyle \ sigma}\ sigma , то есть параметра σ = x 3 и N = 2.

В контексте вариационного исчисления это можно записать как

δ S = δ ∫ AB nds = δ ∫ x 3 A x 3 B ndsdx 3 dx 3 = δ ∫ Икс 3 A Икс 3 BL (Икс 1, Икс 2, Икс ˙ 1, Икс ˙ 2, Икс 3) dx 3 = 0 {\ displaystyle \ delta S = \ delta \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} n \, ds = \ delta \ int _ {x_ {3A}} ^ {x_ {3B}} n {\ frac {ds} {dx_ {3}}} \, dx_ {3} = \ дельта \ int _ {x_ {3A}} ^ {x_ {3B}} L \ left (x_ {1}, x_ {2}, {\ dot {x}} _ {1}, {\ dot {x}} _ {2}, x_ {3} \ right) \, dx_ {3} = 0}\ delta S = \ delta \ int _ {{{\ mathbf {A}}}} ^ {{{\ mathbf {B}}}} n \, ds = \ delta \ int _ {{x _ {{3A}}}} ^ {{x _ {{3B}}}} n {\ frac {ds} {dx_ {3}}} \, dx_ {3} = \ delta \ int _ {{x _ {{3A}}}} ^ {{x _ {{3B}}}} L \ left (x_ {1}, x_ {2}, {\ dot {x}} _ {1}, {\ точка {x}} _ {2}, x_ {3} \ right) \, dx_ {3} = 0

где ds - бесконечно малое смещение вдоль луча, задаваемое формулой ds = dx 1 2 + dx 2 2 + dx 3 2 {\ displaystyle ds = {\ sqrt {dx_ {1} ^ {2} + dx_ {2} ^ {2} + dx_ {3} ^ {2}}}}ds = {\ sqrt {dx_ {1} ^ {2} + dx_ {2} ^ {2} + dx_ { 3} ^ {2}}} и

L знак равно ndsdx 3 знак равно N (Икс 1, Икс 2, Икс 3) 1 + Икс ˙ 1 2 + Икс ˙ 2 2 {\ Displaystyle L = n {\ frac {ds} {dx_ {3}}} = n \ left ( x_ {1}, x_ {2}, x_ {3} \ right) {\ sqrt {1 + {\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2}}}}L = n {\ frac {ds} {dx_ {3}}} = n \ left (x_ {1}, x_ {2}, x_ {3} \ right) {\ sqrt {1 + {\ dot {x}} _ {1} ^ {2} + { \ dot {x}} _ {2} ^ {2}}}

- оптический лагранжиан, а x ˙ k = dxk / dx 3 {\ displaystyle {\ dot {x}} _ {k} = dx_ {k} / dx_ {3} }\ dot {x} _k = dx_k / dx_3 .

длина оптического пути (OPL) определяется как

S = ∫ AB nd s = ∫ ABL dx 3 {\ Displaystyle S = \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} n \, ds = \ int _ {\ mathbf {A}} ^ {\ mathbf {B }} L \, dx_ {3}}S = \ int _ {{{\ mathbf {A}}}} ^ {{{\ mathbf {B}}}} n \, ds = \ int _ {{{\ mathbf {A}}}} ^ {{{{\ mathbf {B}}}} L \, dx_ {3}

где n - локальный показатель преломления как функция положения на пути между точками A и B.

Уравнения Эйлера-Лагранжа

Общие результаты, представленные выше для принципа Гамильтона, могут быть применены к оптике с использованием лагранжиана, определенного в принципе Ферма. Уравнения Эйлера-Лагранжа с параметром σ = x 3 и N = 2, примененные к принципу Ферма, приводят к

∂ L ∂ xk - ddx 3 ∂ L ∂ x ˙ k = 0 {\ displaystyle {\ frac {\ partial L} {\ partial x_ {k}}} - {\ frac {d} {dx_ {3}}} {\ frac {\ partial L} {\ partial {\ dot {x}} _ {k }}} = 0}{\ frac {\ partial L} {\ partial x_ {k}}} - {\ frac {d} {dx_ {3}}} { \ frac {\ partial L} {\ partial {\ dot x} _ {k}}} = 0

с k = 1,2 и где L - оптический лагранжиан, а x ˙ k = dxk / dx 3 {\ displaystyle {\ dot {x}} _ {k} = dx_ {k} / dx_ {3}}\ dot {x} _k = dx_k / dx_3 .

Оптический импульс

Оптический импульс определяется как

pk = ∂ L ∂ x ˙ k {\ displaystyle p_ {k} = {\ frac {\ частичное L} {\ partial {\ dot {x}} _ {k}}}}p_ {k} = {\ frac {\ partial L} {\ partial {\ dot x } _ {k}}}

и из определения оптического лагранжиана L = n 1 + x ˙ 1 2 + x ˙ 2 2 {\ стиль отображения L = п {\ sqrt {1 + {\ точка {х}} _ {1} ^ {2} + {\ точка {x}} _ {2} ^ {2}}}}L = n {\ sqrt {1 + {\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2}}} это выражение можно переписать как

pk = nx ˙ kx ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 = ndxkdx 1 2 + dx 2 2 + dx 3 2 = ndxkds {\ displaystyle p_ {k} = n {\ frac {{\ dot {x}} _ {k}} {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2 } + {\ dot {x}} _ {3} ^ {2}}}} = n {\ frac {dx_ {k}} {\ sqrt {dx_ {1) } ^ {2} + dx_ {2} ^ {2} + dx_ {3} ^ {2}}}} = n {\ frac {dx_ {k}} {ds}}}p_ {k} = n {\ frac {{\ dot {x}}} _ {k}} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ { 3} ^ {2}}}}} = n {\ frac {dx_ {k}} {{\ sqrt {dx_ {1} ^ {2} + dx_ {2} ^ {2} + dx_ {3} ^ { 2}}}}} = n {\ fr ac {dx_ {k}} {ds}}
Оптический импульс

или в векторной форме

p = ndsds = (p 1, p 2, p 3) {\ displaystyle \ mathbf {p} = n {\ frac {\ mathbf {ds}} {ds}} = \ left (p_ {1}, p_ {2}, p_ {3} \ right)}{\ mathbf {p}} = n {\ frac { {\ mathbf {ds}}} {ds}} = \ left (p_ {1}, p_ {2}, p_ {3} \ right)
= (n cos ⁡ α 1, n cos ⁡ α 2, n cos ⁡ α 3) = ne ^ {\ displaystyle = \ left ( n \ cos \ alpha _ {1}, n \ cos \ alpha _ {2}, n \ cos \ alpha _ {3} \ right) = n \ mathbf {\ hat {e}}}= \ left (n \ cos \ alpha _ {1}, n \ cos \ alpha _ {2}, n \ cos \ alpha _ {3} \ right) = n {\ mathbf {{\ hat {e}}}}

где e ^ {\ displaystyle \ mathbf {\ hat {e}}}\ mathbf {\ hat {e}} - это единичный вектор и углы α 1, α 2 и α 3 - углы p относительно оси x 1, x 2 и x 3 соответственно, как показано на рисунке «оптический момент». Следовательно, оптический импульс - это вектор norm

‖ p ‖ = p 1 2 + p 2 2 + p 3 2 = n {\ displaystyle \ | \ mathbf {p} \ | = {\ sqrt { p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2}}} = n}\ | {\ mathbf {p}} \ | = {\ sqrt {p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2}}} = n

, где n - показатель преломления, при котором p рассчитано. Вектор p указывает направление распространения света. Если свет распространяется в оптике с градиентным показателем преломления , путь светового луча изогнут, и вектор p касается светового луча.

Выражение для длины оптического пути также можно записать как функцию оптического импульса. Учитывая, что x ˙ 3 = dx 3 / dx 3 = 1 {\ displaystyle {\ dot {x}} _ {3} = dx_ {3} / dx_ {3} = 1}{\ dot { x}} _ {3} = dx_ {3} / dx_ {3} = 1 выражение для оптического лагранжиана можно переписать как

L = nx ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 = x ˙ 1 nx ˙ 1 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 + x ˙ 2 nx ˙ 2 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 + nx ˙ 3 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 {\ displaystyle L = n {\ sqrt { {\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}} = {\ dot {x}} _ {1} {\ frac {n {\ dot {x}} _ {1}} {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ точка {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}} + {\ dot {x}} _ {2} {\ frac {n { \ dot {x}} _ {2}} {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}} + {\ frac {n {\ dot {x}} _ {3}} {\ sqrt {{\ dot {x}} _ {1} ^ { 2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}}L = n {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2 } ^ {2} + {\ dot {x}} _ {3} ^ {2}}} = {\ dot {x}} _ {1} {\ frac {n {\ dot {x}} _ {1 }} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}} + {\ dot {x}} _ {2} {\ frac {n {\ dot {x}} _ {2}} {{\ sqrt {{\ dot {x}}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}} + {\ frac {n { \ dot {x}} _ {3}} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ точка {x}} _ {3} ^ {2}}}}}
= x ˙ 1 p 1 + x ˙ 2 п 2 + Икс ˙ 3 п 3 знак равно Икс ˙ 1 п 1 + Икс ˙ 2 п 2 + п 3 {\ displaystyle = {\ dot {x}} _ {1} p_ {1} + {\ dot {x} } _ {2} p_ {2} + {\ dot {x}} _ {3} p_ {3} = {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + p_ {3}}= {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + {\ точка {x}} _ {3} p_ {3} = {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + p_ {3}

и выражение для длины оптического пути:

S = ∫ L dx 3 = ∫ p ⋅ ds {\ displaystyle S = \ int L \, dx_ {3} = \ int \ mathbf {p} \ cdot d \ mathbf {s}}S = \ int L \, dx_ {3} = \ int {\ mathbf {p}} \ cdot d {\ mathbf {s}}

уравнения Гамильтона

Аналогично что происходит в гамильтоновой механике, а также в оптике гамильтониан определяется выражением, приведенным выше для N = 2, соответствующих функциям x 1 (x 3) {\ displaystyle x_ {1} \ left (x_ {3} \ right)}x_ {1 } \ left (x_ {3} \ right) и x 2 (x 3) {\ displaystyle x_ {2} \ left (x_ {3} \ right)}x_ {2} \ left (x_ {3} \ right) подлежит определению

H = x ˙ 1 p 1 + x ˙ 2 p 2 - L {\ displaystyle H = {\ dot {x}} _ {1} p_ {1} + {\ dot { x}} _ {2} p_ {2} -L}H = {\ dot {x}} _ {1} p_ {1} + {\ dot {x }} _ {2} p_ {2} -L

Сравнение этого выражения с L = x ˙ 1 p 1 + x ˙ 2 p 2 + p 3 {\ displaystyle L = {\ dot {x }} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + p_ {3}}L = {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + p_ {3 } для лагранжевых результатов в

H = - p 3 = - n 2 - p 1 2 - p 2 2 {\ displaystyle H = -p_ {3} = - {\ sqrt {n ^ {2} -p_ {1} ^ {2} -p_ {2} ^ { 2}}}}H = -p_ {3} = - {\ sqrt {n ^ {2} -p_ {1} ^ {2} -p_ {2} ^ {2}}}

И соответствующие уравнения Гамильтона с параметром σ = x 3 и k = 1,2, примененные к оптике:

∂ H ∂ xk = - p ˙ k, ∂ H ∂ пк знак равно Икс ˙ К {\ Displaystyle {\ fra c {\ partial H} {\ partial x_ {k}}} = - {\ dot {p}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial p_ {k}}} = {\ dot {x}} _ {k}}{\ frac {\ partial H} {\ partial x_ {k}}} = - {\ dot {p}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial p_ {k}}} = {\ dot {x}} _ {k}

с x ˙ k = dxk / dx 3 {\ displaystyle {\ dot {x}} _ {k} = dx_ {k} / dx_ { 3}}\ dot {x} _k = dx_k / dx_3 и p ˙ k = dpk / dx 3 {\ displaystyle {\ dot {p}} _ {k} = dp_ {k} / dx_ {3}}{\ dot { p}} _ {k} = dp_ {k} / dx_ {3} .

Приложения

Предполагается, что свет движется вдоль оси x 3, в принципе Гамильтона выше, координаты x 1 {\ displaystyle x_ {1}}x_ {1} и x 2 {\ displaystyle x_ {2}}x_{2}берут на себя роль обобщенных координат qk {\ displaystyle q_ {k}}q_ {k} а x 3 {\ displaystyle x_ {3}}x_ {3} играет роль параметра σ {\ displaystyle \ sigma}\ sigma , то есть параметра σ = x 3 и N = 2.

Преломление и отражение

Если плоскость x 1x2разделяет две среды с показателем преломления n A ниже и n B над ней, коэффициент преломления индекс задается функцией шага

n (x 3) = {n A if x 3 < 0 n B if x 3>0 {\ displaystyle n (x_ {3}) = {\ begin {cases} n_ {A} {\ t_dv {if}} x_ {3} <0\\n_{B}{\t_dv{if }}x_{3}>0 \\\ end {cases}}}n(x_{3})={\begin{cases}n_{A}{\t_dv{if }}x_{3}<0\\n_{B}{\t_dv{if }}x_{3}>0 \\\ end {cases}}

и из уравнений Гамильтона

∂ H ∂ xk = ∂ ∂ xkn (x 3) 2 - p 1 2 - p 2 2 = 0 {\ displaystyle {\ frac {\ partial H} {\ partial x_ {k}}} = - {\ frac {\ partial} {\ partial x_ {k}}} {\ sqrt {n (x_ {3}) ^ {2} -p_ {1} ^ {2} -p_ {2} ^ {2}}} = 0}{\ frac {\ partial H} {\ partial x_ {k}}} = - {\ frac {\ partial} {\ partial x_ {k}}} {\ sqrt {n (x_ {3}) ^ {2} -p_ {1} ^ {2} -p_ {2} ^ {2}}} = 0

и, следовательно, п ˙ К = 0 {\ displaystyle {\ dot {p}} _ {k} = 0}{\ dot {p}} _ {k} = 0 или pk = Constant {\ displaystyle p_ {k} = {\ text {Constant) }}}{\ displaystyle p_ {k} = {\ text {Constant}}} для k = 1,2.

Входящий световой луч имеет импульс pAдо рефракции (ниже плоскости x 1x2) и импульс pBпосле рефракции ция (над плоскостью x 1x2). Луч света составляет угол θ A с осью x 3 (нормаль к преломляющей поверхности) до преломления и угол θ B с осью x 3 после рефракции. Поскольку компоненты импульса p 1 и p 2 постоянны, только p 3 изменяется с p 3A на p 3B.

Преломление

Рисунок «преломление» показывает геометрию этого преломления, из которой d = ‖ p A ‖ sin ⁡ θ A = ‖ p B ‖ sin ⁡ θ B {\ displaystyle d = \ | \ mathbf {p} _ {A} \ | \ sin \ theta _ {A} = \ | \ mathbf {p} _ {B} \ | \ sin \ theta _ {B}}d = \ | {\ mathbf {p}} _ {A} \ | \ sin \ theta _ {A} = \ | {\ mathbf {p}} _ {B} \ | \ sin \ theta _ {B} . Поскольку ‖ p A ‖ = n A {\ displaystyle \ | \ mathbf {p} _ {A} \ | = n_ {A}}\ | {\ mathbf {p}} _ {A} \ | = n_ {A} и ‖ p B ‖ = n B {\ displaystyle \ | \ mathbf {p} _ {B} \ | = n_ {B}}\ | {\ mathbf {p}} _ {B } \ | = n_ {B} , это последнее выражение можно записать как

n A sin ⁡ θ A = n B sin ⁡ θ B {\ displaystyle n_ {A} \ sin \ theta _ {A} = n_ {B} \ sin \ theta _ {B}}{\ displaystyle n_ {A} \ sin \ theta _ {A} = n_ {B} \ sin \ theta _ {B}}

который является законом Снеллиуса преломления.

На рисунке «преломление» нормаль к преломляющей поверхности указывает в направлении оси x 3, а также вектора v = p A - p B {\ displaystyle \ mathbf {v } = \ mathbf {p} _ {A} - \ mathbf {p} _ {B}}{\ mathbf {v}} = {\ mathbf {p}} _ {A} - {\ mathbf {p}} _ {B} . Единичная нормаль n = v / ‖ v ‖ {\ displaystyle \ mathbf {n} = \ mathbf {v} / \ | \ mathbf {v} \ |}{\ mathbf {n}} = {\ mathbf {v}} / \ | {\ mathbf {v}} \ | к преломляющей поверхности может тогда быть полученным из импульсов входящих и выходящих лучей следующим образом:

n = p A - p B ‖ p A - p B ‖ = n A i - n B r ‖ n A i - n B r ‖ {\ displaystyle \ mathbf {n} = {\ frac {\ mathbf {p} _ {A} - \ mathbf {p} _ {B}} {\ | \ mathbf {p} _ {A} - \ mathbf {p} _ {B } \ |}} = {\ frac {n_ {A} \ mathbf {i} -n_ {B} \ mathbf {r}} {\ | n_ {A} \ mathbf {i} -n_ {B} \ mathbf { r} \ |}}}{\ mathbf {n}} = {\ frac {{\ mathbf {p}} _ {A} - {\ mathbf {p}} _ {B}} {\ | {\ mathbf {p}} _ {A} - {\ mathbf {p} } _ {B} \ |}} = {\ frac {n_ {A} {\ mathbf {i}} - n_ {B} {\ mathbf {r}}} {\ | n_ {A} {\ mathbf {i}} - n_ {B} {\ mathbf {r} } \ |}}

где i и r - единичные векторы в направлениях падающих и преломленных лучей. Кроме того, исходящий луч (в направлении p B {\ displaystyle \ mathbf {p} _ {B}}{\ mathbf { p}} _ {B} ) содержится в плоскости, определяемой входящим лучом (в направлении p A {\ displaystyle \ mathbf {p} _ {A}}{\ mathbf {p}} _ {A} ) и нормальный n {\ displaystyle \ mathbf {n}}\ mathbf {n} на поверхность.

Аналогичный аргумент может использоваться для отражения при выводе закона зеркального отражения, только теперь с n A=nB, что приводит к θ A=θB. Кроме того, если i и r являются единичными векторами в направлениях падающего и преломленного луча соответственно, соответствующая нормаль к поверхности задается тем же выражением, что и для преломления, только с п A=nB

N = я - р ‖ я - р ‖ {\ displaystyle \ mathbf {n} = {\ frac {\ mathbf {i} - \ mathbf {r}} {\ | \ mathbf {i} - \ mathbf {r} \ |}}}{\ mathbf {n}} = {\ frac {{\ mathbf {i}} - {\ mathbf {r}}} {\ | {\ mathbf {i}} - {\ mathbf {r}} \ |}}

В векторной форме, если i - единичный вектор, указывающий в направлении падающего луча, а n - единичная нормаль к поверхности, направление r преломленного луча определяется выражением:

r = n A n B i + (- (i ⋅ n) n A n B + Δ) n {\ displaystyle \ mathbf { r} = {\ frac {n_ {A}} {n_ {B}}} \ mathbf {i} + \ left (- \ left (\ mathbf {i} \ cdot \ mathbf {n} \ right) {\ frac {n_ {A}} {n_ {B}}} + {\ sqrt {\ Delta}} \ right) \ mathbf {n}}{\ mathbf {r}} = {\ frac {n_ {A}} {n_ {B}}} {\ mathbf {i}} + \ left (- \ left ({\ mathbf {i}} \ cdot { \ mathbf {n}} \ right) {\ frac {n_ {A}} {n_ {B}}} + {\ sqrt {\ Delta}} \ right) {\ mathbf {n}}

с

Δ = 1 - (n A n B) 2 ( 1 - (я ⋅ n) 2) {\ displaystyle \ Delta = 1- \ left ({\ frac {n_ {A}} {n_ {B}}} \ right) ^ {2} \ left (1- \ left (\ mathbf {i} \ cdot \ mathbf {n} \ right) ^ {2} \ right)}\ Delta = 1- \ left ({\ frac {n_ {A}} {n_ {B}) }} \ right) ^ {2} \ left (1- \ left ({\ mathbf {i}} \ cdot {\ mathbf {n}} \ right) ^ {2} \ right)

Если i·n<0 then -nследует использовать в вычислении ионы. Когда Δ < 0 {\displaystyle \Delta <0}\ Delta <0, свет испытывает полное внутреннее отражение, и выражение для отраженного луча является выражением отражения:

r = i - 2 (i ⋅ n) n {\ displaystyle \ mathbf {r} = \ mathbf {i} -2 \ left (\ mathbf {i} \ cdot \ mathbf {n} \ right) \ mathbf {n}}{\ displaystyle \ mathbf {r} = \ mathbf {i} -2 \ left (\ mathbf {i} \ cdot \ mathbf {n} \ right) \ mathbf {n}}

Лучи и волновые фронты

Из определения оптического пути длина S = ∫ L dx 3 {\ displaystyle S = \ int L \, dx_ {3}}{\ displaystyle S = \ int L \, dx_ {3}}

∂ S ∂ xk = ∫ ∂ L ∂ xkdx 3 = ∫ dpkdx 3 dx 3 = pk {\ displaystyle {\ frac {\ partial S} {\ partial x_ {k}}} = \ int {\ frac {\ partial L} {\ partial x_ {k}}} \, dx_ {3} = \ int {\ frac { dp_ {k}} {dx_ {3}}} \, dx_ {3} = p_ {k}}{ \ frac {\ partial S} {\ partial x_ {k}}} = \ int {\ frac {\ partial L} {\ partial x_ {k}}} \, dx_ {3} = \ int {\ frac {dp_ {k}} {dx_ {3}}} \, dx_ {3} = p_ {k}
Лучи и волновые фронты

с k = 1,2, где уравнения Эйлера-Лагранжа ∂ L / ∂ xk = dpk / dx 3 {\ displaystyle \ partial L / \ partial x_ {k} = dp_ {k} / dx_ {3}}\ partial L / \ partial x_ {k} = dp_ { k} / dx_ {3} с k = 1,2 были использованы. Кроме того, из последнего из уравнений Гамильтона ∂ H / ∂ x 3 = - ∂ L / ∂ x 3 {\ displaystyle \ partial H / \ partial x_ {3} = - \ partial L / \ partial x_ {3}}\ partial H / \ partial x_ {3} = - \ partial L / \ partial x_ {3} и из H = - p 3 {\ displaystyle H = -p_ {3}}{\ displaystyle H = -p_ {3}} выше

∂ S ∂ x 3 = ∫ ∂ L ∂ Икс 3 dx 3 знак равно ∫ dp 3 dx 3 dx 3 = p 3 {\ displaystyle {\ frac {\ partial S} {\ partial x_ {3}}} = \ int {\ frac {\ partial L} {\ partial x_ {3}}} \, dx_ {3} = \ int {\ frac {dp_ {3}} {dx_ {3}}} \, dx_ {3} = p_ {3}}{\ frac {\ partial S} {\ partial x_ {3}}} = \ int {\ frac {\ partial L} {\ partial x_ {3}}} \, dx_ {3} = \ int {\ frac {dp_ {3}} {dx_ {3}}} \, dx_ {3} = p_ {3}

объединяя уравнения для компоненты импульса p приводят к

p = ∇ S {\ displaystyle \ mathbf {p} = \ nabla S}{\ mathbf {p}} = \ nabla S

Поскольку p является вектором, касательным к световые лучи, поверхности S = ​​Constant должны быть перпендикулярны этим световым лучам. Эти поверхности называются волновыми фронтами. Рисунок «Лучи и волновые фронты» иллюстрирует эту взаимосвязь. Также показан оптический импульс p, касательный к световому лучу и перпендикулярный волновому фронту.

Векторное поле p = ∇ S {\ displaystyle \ mathbf {p} = \ nabla S}{\ mathbf {p}} = \ nabla S - это консервативное векторное поле. Теорема о градиенте затем может быть применена к длине оптического пути (как указано выше), в результате чего

S = ∫ AB p ⋅ ds = ∫ AB ∇ S ⋅ ds = S (В) - S (A) {\ Displaystyle S = \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} \ mathbf {p} \ cdot d \ mathbf {s} = \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} \ nabla S \ cdot d \ mathbf {s} = S (\ mathbf {B}) -S (\ mathbf {A})}S = \ int _ {{{\ mathbf {A}}}} ^ {{{\ mathbf {B} }}} {\ mathbf {p}} \ cdot d {\ mathbf {s}} = \ int _ {{{\ mathbf {A}}}} ^ {{{{\ mathbf {B}}}} \ nabla S \ cdot d {\ mathbf {s}} = S ({\ mathbf {B}}) - S ({\ mathbf {A}})

и оптический путь длина S, вычисленная вдоль кривой C между точками A и B, является функцией только ее конечных точек A и B, а не форма кривой между ними. В частности, если кривая замкнута, она начинается и заканчивается в одной и той же точке, или A=B, так что

S = ∮ ∇ S ⋅ ds = 0 {\ displaystyle S = \ oint \ nabla S \ cdot d \ mathbf {s} = 0}S = \ oint \ nabla S \ cdot d {\ mathbf {s}} = 0

Этот результат может быть применен к замкнутому пути ABCDA, как на рисунке «длина оптического пути»

S = ∫ AB p ⋅ ds + ∫ BC p ⋅ ds + ∫ CD п ⋅ DS + ∫ DA п ⋅ ds = 0 {\ displaystyle S = \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} \ mathbf {p} \ cdot d \ mathbf {s} + \ int _ {\ mathbf {B}} ^ {\ mathbf {C}} \ mathbf {p} \ cdot d \ mathbf {s} + \ int _ {\ mathbf {C}} ^ {\ mathbf {D} } \ mathbf {p} \ cdot d \ mathbf {s} + \ int _ {\ mathbf {D}} ^ {\ mathbf {A}} \ mathbf {p} \ cdot d \ mathbf {s} = 0}S = \ int _ {{{\ mathbf {A}}}} ^ {{{{\ mathbf {B}}}}} {\ mathbf {p}} \ cdot d { \ mathbf {s}} + \ int _ {{{\ mathbf {B}}}} ^ {{{\ mathbf {C}}}} {\ mathbf {p}} \ cdot d {\ mathbf {s}} + \ int _ {{{\ mathbf {C}}}} ^ {{{\ mathbf {D}}}} {\ mathbf {p}} \ cdot d {\ mathbf {s}} + \ int _ {{ {\ mathbf {D}}}} ^ {{{\ mathbf {A}}}} {\ mathbf {p}} \ cdot d {\ mathbf {s}} = 0
Длина оптического пути

для сегмента кривой AB оптический момент p перпендикулярен смещению d s вдоль кривой AB или p ⋅ ds = 0 {\ displaystyle \ mathbf {p} \ cdot d \ mathbf {s} = 0}{\ mathbf {p }} \ cdot d {\ mathbf {s}} = 0 . То же верно и для сегмента CD . Для сегмента BC оптический момент p имеет то же направление, что и смещение d s и p ⋅ ds = nds {\ displaystyle \ mathbf {p} \ cdot d \ mathbf {s} = nds}{\ mathbf {p}} \ cdot d {\ mathbf {s}} = nds . Для сегмента DA оптический момент p имеет направление, противоположное смещению d s и p ⋅ ds = - nds {\ displaystyle \ mathbf {p } \ cdot d \ mathbf {s} = -n \, ds}{\ mathbf {p}} \ cdot d {\ mathbf {s}} = - п \, дс . Однако при инвертировании направления интегрирования, так что интеграл берется от A до D, d s меняет направление и p ⋅ ds = nds { \ Displaystyle \ mathbf {p} \ cdot d \ mathbf {s} = n \, ds}{\ mathbf {p}} \ cdot d {\ mathbf {s}} = n \, ds . Исходя из этих соображений

∫ BC nds = ∫ AD nds {\ displaystyle \ int _ {\ mathbf {B}} ^ {\ mathbf {C}} n \, ds = \ int _ {\ mathbf {A}} ^ {\ mathbf {D}} n \, ds}\ int _ {{{\ mathbf {B}}}} ^ {{{\ mathbf {C}}}} n \, ds = \ int _ {{{\ mathbf {A}} }} ^ {{{\ mathbf {D}}}} n \, ds

или

SBC = SAD {\ displaystyle S _ {\ mathbf {BC}} = S _ {\ mathbf {AD}}}S _ {{\ mathbf {BC}}} = S _ {{\ mathbf {AD}}}

и оптический путь длина S BCмежду точками B и C вдоль соединяющего их луча такая же, как длина оптического пути S ADмежду точками A и D вдоль луча, соединяющего их. Длина оптического пути постоянна между фронтами волн.

Фазовое пространство

На рисунке «2D-фазовое пространство» вверху показаны некоторые световые лучи в двухмерном пространстве. Здесь x 2 = 0 и p 2 = 0, поэтому свет распространяется по плоскости x 1x3в направлениях увеличения значений x 3. В этом случае p 1 2 + p 3 2 = n 2 {\ displaystyle p_ {1} ^ {2} + p_ {3} ^ {2} = n ^ {2}}p_ {1} ^ {2 } + p_ {3} ^ {2} = n ^ {2} и направление светового луча полностью определяется компонентой импульса p 1 p = (p 1, p 3) {\ displaystyle \ mathbf {p} = (p_ {1}, p_ {3})}{\ mathbf {p}} = (p_ {1}, p_ {3}) , поскольку p 2 = 0. Если задано p 1, p 3 может быть вычислено (с учетом значения показателя преломления n), и поэтому p 1 достаточно для определения направления луч света. Показатель преломления среды, в которой проходит луч, определяется как ‖ p ‖ = n {\ displaystyle \ | \ mathbf {p} \ | = n}\ | { \ mathbf {p}} \ | = n .

2D-фазовое пространство

Например, луч r C пересекает ось x 1 в координате x B с оптическим импульсом pC, вершина которого находится на окружности радиуса n с центром в позиции x В. Координата x B и горизонтальная координата p 1C количества движения pCполностью определяют луч r C, когда он пересекает ось x 1. Затем этот луч может быть определен точкой r C = (x B,p1C) в пространстве x 1p1, как показано внизу рисунка. Пространство x 1p1называется фазовым пространством, и разные световые лучи могут быть представлены разными точками в этом пространстве.

Таким образом, луч r D, показанный вверху, представлен точкой r D в фазовом пространстве внизу. Все лучи, пересекающие ось x 1 с координатой x B между лучами r C и r D, представлены вертикальной линией, соединяющей точки r C и r D в фазовом пространстве. Соответственно, все лучи, пересекающие ось x 1 с координатой x A, находящиеся между лучами r A и r B, представлены вертикальной линией. точки соединения r A и r B в фазовом пространстве. В общем, все лучи, пересекающие ось x 1 между x L и x R, представлены объемом R в фазовом пространстве. Лучи на границе ∂R объема R называются краевыми лучами. Например, в позиции x A оси x 1 лучи r A и r B являются краевыми лучами, поскольку все остальные лучи находятся между этими двумя. (Луч, параллельный x1, не будет между двумя лучами, так как импульс не находится между двумя лучами)

.

В трехмерной геометрии оптический импульс задается как p = (p 1, p 2, p 3) {\ displaystyle \ mathbf {p} = (p_ {1}, p_ {2}, p_ {3})}{\ mathbf {p}} = (p_ {1}, p_ {2}, p_ {3}) с p 1 2 + p 2 2 + p 3 2 = n 2 {\ displaystyle p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2} = n ^ {2}}p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2} = n ^ {2} . Если p 1 и p 2 заданы, p 3 может быть вычислено (с учетом значения показателя преломления n) и, следовательно, p 1 и p 2 достаточно для определения направления светового луча. Затем луч, движущийся вдоль оси x 3, определяется точкой (x 1,x2) в плоскости x 1x2и направлением (p 1,p2). Затем он может быть определен точкой в ​​четырехмерном фазовом пространстве x1x2p1p2.

Сохранение внешнего вида

На рисунке «изменение объема» показан объем V, ограниченный областью A. Со временем, если граница A движется, объем V может меняться. В частности, бесконечно малая область dA с направленной наружу нормалью n движется со скоростью v.

Изменение объема

Это приводит к изменению объема d V = d A (v ⋅ n) dt {\ displaystyle dV = dA (\ mathbf {v} \ cdot \ mathbf {n}) dt}dV = dA ({\ mathbf {v}} \ cdot {\ mathbf {n}}) dt . Используя теорему Гаусса, изменение во времени полного объема V объема, движущегося в пространстве, равно

d V dt = ∫ A v ⋅ nd A = ∫ V ∇ ⋅ vd V {\ displaystyle { \ frac {dV} {dt}} = \ int _ {A} \ mathbf {v} \ cdot \ mathbf {n} \, dA = \ int _ {V} \ nabla \ cdot \ mathbf {v} \, dV }{\ frac {dV} {dt}} = \ int _ {A} {\ mathbf {v}} \ cdot {\ mathbf {n}} \, dA = \ int _ {V} \ nabla \ cdot {\ mathbf {v}} \, dV

Самый правый член - это объемный интеграл по объему V, а средний член - поверхностный интеграл по границе A объема V. Кроме того, v - это скорость, с которой движутся точки в V.

В оптике координата x 3 {\ displaystyle x_ {3}}x_ {3} играет роль времени. В фазовом пространстве луч света идентифицируется точкой (x 1, x 2, p 1, p 2) {\ displaystyle (x_ {1}, x_ {2}, p_ {1}, p_ {2}))}{\ displaystyle (x_ {1}, x_ {2}, p_ {1}, p_ {2})} который движется со "скоростью " v = (x ˙ 1, x ˙ 2, p ˙ 1, p ˙ 2) {\ displaystyle \ mathbf {v } = ({\ dot {x}} _ {1}, {\ dot {x}} _ {2}, {\ dot {p}} _ {1}, {\ dot {p}} _ {2})}{\ displaystyle \ mathbf {v} = ({\ dot {x}} _ {1 }, {\ dot {x}} _ {2}, {\ dot {p}} _ {1}, {\ dot {p}} _ {2})} , где точка представляет собой производную относительно x 3 {\ displaystyle x_ {3}}x_ {3} . Набор световых лучей, распространяющихся по dx 1 {\ displaystyle dx_ {1}}{\ displaystyle dx_ {1}} в координатах x 1 {\ displaystyle x_ {1}}x_ {1} , dx 2 {\ displaystyle dx_ {2}}{\ displaystyle dx_ {2}} в координате x 2 {\ displaystyle x_ {2}}x_{2}, dp 1 {\ displaystyle dp_ {1}}{\ displaystyle dp_ {1}} в координате p 1 {\ displaystyle p_ {1}}p_ {1} и dp 2 {\ displaystyle dp_ {2}}{\ displaystyle dp_ {2}} в координатах p 2 {\ displaystyle p_ {2}}p_ {2} занимает объем d V = dx 1 dx 2 dp 1 dp 2 {\ displaystyle dV = dx_ {1} dx_ {2} dp_ {1} dp_ {2}}{\ displaystyle dV = dx_ {1} dx_ {2} dp_ {1} dp_ {2}} в фазовом пространстве. В общем, большой набор лучей занимает большой объем V {\ displaystyle V}V в фазовом пространстве, к которому теорема Гаусса может быть применена

d V dx 3 = ∫ V ∇ ⋅ vd V {\ displaystyle {\ frac {dV} {dx_ {3}}} = \ int _ {V} \ nabla \ cdot \ mathbf {v} \, dV}{\ frac {dV} {dx_ {3}}} = \ int _ {V} \ nabla \ cdot {\ mathbf {v}} \, dV

и с использованием Уравнения Гамильтона

∇ ⋅ v = ∂ x ˙ 1 ∂ x 1 + ∂ x ˙ 2 ∂ x 2 + ∂ p ˙ 1 ∂ p 1 + ∂ p ˙ 2 ∂ p 2 = ∂ ∂ x 1 ∂ H ∂ p 1 + ∂ ∂ x 2 ∂ H ∂ p 2 - ∂ ∂ p 1 ∂ H ∂ x 1 - ∂ ∂ p 2 ∂ H ∂ x 2 = 0 {\ displaystyle \ nabla \ cdot \ mathbf {v} = {\ frac { \ partial {\ dot {x}} _ {1}} {\ partial x_ {1}}} + {\ frac {\ partial {\ dot {x}} _ {2}} {\ partial x_ {2}} } + {\ frac {\ partial {\ dot {p}} _ {1}} {\ partial p_ {1}}} + {\ frac {\ partial {\ dot {p}} _ {2}} {\ частичное p_ {2}}} = {\ frac {\ partial} {\ partial x_ {1}}} {\ frac {\ partial H} {\ partial p_ {1}}} + {\ frac {\ partial} { \ partial x_ {2}}} {\ frac {\ partial H} {\ partial p_ {2}}} - {\ frac {\ partial} {\ partial p_ {1}}} {\ frac {\ partial H} {\ partial x_ {1}}} - {\ frac {\ partial} {\ partial p_ {2}}} {\ frac {\ partial H} {\ partial x_ {2}}} = 0}\ nabla \ cdot {\ mathbf {v}} = {\ frac {\ partial { \ dot {x}} _ {1}} {\ partial x_ {1}}} + {\ frac {\ partial {\ dot {x}} _ {2}} {\ partial x_ {2}}} + { \ frac {\ partial {\ dot {p}} _ {1}} {\ partial p_ {1}}} + {\ frac {\ partial {\ dot {p}} _ {2}} {\ partial p_ { 2}}} = {\ frac {\ partial} {\ partial x_ {1}}} {\ frac {\ partial H} {\ partial p_ {1}}} + {\ frac {\ partial} {\ partial x_ { 2}}} {\ frac {\ partial H} {\ partial p_ {2}}} - {\ frac {\ partial} {\ partial p_ {1}}} {\ frac {\ partial H} {\ partial x_ {1}}} - {\ frac {\ partial} {\ partial p_ {2}}} {\ frac {\ partial H} {\ partial x_ {2}}} = 0

или d V / dx 3 = 0 {\ displaystyle dV / dx_ {3} = 0}{\ displaystyle dV / dx_ {3} = 0} и d V = dx 1 dx 2 dp 1 dp 2 = константа {\ displaystyle dV = dx_ {1} dx_ {2} dp_ {1} dp_ {2} = {\ text {Constant}}{\ displaystyle dV = dx_ {1} dx_ {2} dp_ {1} dp_ {2} = {\ text {Константа}}} , что означает, что объем фазового пространства сохраняется при движении света по оптической системе.

Объем, занимаемый набором лучей в фазовом пространстве, называется etendue, который сохраняется по мере прохождения световых лучей в оптической системе в направлении x 3. Это соответствует теореме Лиувилля, которая также применима к гамильтоновой механике.

Однако смысл теоремы Лиувилля в механике весьма отличается от теоремы о сохранении эталонности. Теорема Лиувилля является статистической по своей природе и относится к эволюции во времени ансамбля механических систем с идентичными свойствами, но с разными начальными условиями. Каждая система представлена ​​одной точкой в ​​фазовом пространстве, и теорема утверждает, что средняя плотность точек в фазовом пространстве постоянна во времени. Примером могут служить молекулы идеального классического газа, находящиеся в равновесии в контейнере. Каждая точка в фазовом пространстве, которая в этом примере имеет 2N измерений, где N - количество молекул, представляет собой одну из ансамбля идентичных контейнеров, ансамбля, достаточно большого, чтобы позволить получить статистическое среднее значение плотности репрезентативных точек. Теорема Лиувилля гласит, что если все контейнеры остаются в равновесии, средняя плотность точек остается постоянной.

Оптика с визуализацией и без визуализации

На рисунке «сохранение бесконечности» слева изображена двумерная оптическая система. в которой x 2 = 0 и p 2 = 0, поэтому свет распространяется по плоскости x 1x3в направлениях увеличения значений x 3.

Сохранение внешнего вида

Световые лучи, пересекающие входную апертуру оптики в точке x 1=xI, содержатся между краевыми лучами r A и r B, представленными вертикальной линия между точками r A и r B на фазовом пространстве входной апертуры (правый нижний угол рисунка). Все лучи, пересекающие входную апертуру, представлены в фазовом пространстве областью R I.

Кроме того, световые лучи, пересекающие выходную апертуру оптики в точке x 1=xO, содержатся между краевыми лучами r A и r B представлен вертикальной линией между точками r A и r B в фазовом пространстве выходной апертуры (правый, верхний угол рисунка). Все лучи, пересекающие выходную апертуру, представлены в фазовом пространстве областью R O.

Сохранение длины волны в оптической системе означает, что объем (или область в этом двумерном случае) в фазовом пространстве, занимаемый областью R I на входной апертуре должен быть таким же, как объем в фазовом пространстве, занимаемый R O на выходной апертуре.

В оптике формирования изображений все световые лучи, пересекающие входную апертуру в точке x 1=xI, перенаправляются ею к выходной апертуре в точке x 1=xO, где x I = mx O. Это обеспечивает формирование изображения входа на выходе с увеличением m. В фазовом пространстве это означает, что вертикальные линии в фазовом пространстве на входе преобразуются в вертикальные линии на выходе. Это было бы в случае вертикальной линии r ArBв R I, преобразованной в вертикальную линию r ArBв R O.

В оптике без визуализации, цель не состоит в том, чтобы сформировать изображение, а просто переносить весь свет из входной апертуры в выходную апертуру. Это достигается путем преобразования краевых лучей ∂R I из R I в краевые лучи ∂R O из R O. Это известно как принцип краевого луча.

Обобщения

Выше предполагалось, что свет движется вдоль оси x 3, в принципе Гамильтона выше., координаты x 1 {\ displaystyle x_ {1}}x_ {1} и x 2 {\ displaystyle x_ {2}}x_{2}берут на себя роль обобщенных координат qk {\displaystyle q_{k}}q_ {k} while x 3 {\displaystyle x_{3}}x_ {3} takes the role of parameter σ {\displaystyle \sigma }\ sigma , that is, parameter σ =x3and N=2. However, different parametrizations of the light rays are possible, as well as the use of generalized coordinates.

General ray parametrization

A more general situation can be considered in which the path of a light ray is parametrized as s = ( x 1 ( σ), x 2 ( σ), x 3 ( σ)) {\displaystyle s=\left(x_{1}\left(\sigma \right),x_{2}\left(\sigma \right),x_{3}\left(\sigma \right)\right)}{\ displaystyle s = \ left (x_ {1 } \ left (\ sigma \ right), x_ {2} \ left (\ sig ma \ right), x_ {3} \ left (\ sigma \ right) \ right)} in which σ is a general parameter. In this case, when compared to Hamilton's principleabove, coordinates x 1 {\displaystyle x_{1}}x_ {1} , x 2 {\displaystyle x_{2}}x_{2}and x 3 {\displaystyle x_{3}}x_ {3} take the role of the generalized coordinates q k {\displaystyle q_{k}}q_ {k} with N=3. Applying Hamilton's principleto optics in this case leads to

δ S = δ ∫ A B n d s = δ ∫ σ A σ B n d s d σ d σ {\displaystyle \delta S=\delta \int _{\mathbf {A} }^{\mathbf {B} }n\,ds=\delta \int _{\sigma _{A}}^{\sigma _{B}}n{\frac {ds}{d\sigma }}\,d\sigma }{\ displaystyle \ delta S = \ delta \ int _ {\ mathbf {A}} ^ {\ mathbf {B}} n \, ds = \ delta \ int _ {\ sigma _ {A}} ^ {\ sigma _ {B}} n {\ frac {ds} {d \ sigma}} \, d \ sigma}
= δ ∫ σ A σ B L ( x 1, x 2, x 3, x ˙ 1, x ˙ 2, x ˙ 3, σ) d σ = 0 {\displaystyle =\delta \int _{\sigma _{A}}^{\sigma _{B}}L\left(x_{1},x_{2},x_{3},{\dot {x}}_{1},{\dot {x}}_{2},{\dot {x}}_{3},\sigma \right)\,d\sigma =0}= \ delta \ int _ {{\ sigma _ {A}} } ^ {{\ sigma _ {B}}} L \ left (x_ {1}, x_ {2}, x_ {3}, {\ dot {x}} _ {1}, {\ dot {x}} _ {2}, {\ dot {x}} _ {3}, \ sigma \ right) \, d \ sigma = 0

where now L = n d s / d σ {\displaystyle L=nds/d\sigma }{\ displaystyle L = nds / d \ sigma} and x ˙ k = d x k / d σ {\displaystyle {\dot {x}}_{k}=dx_{k}/d\sigma }{\ dot {x} } _ {k} = dx_ {k} / d \ sigma and for which the Euler-Lagrange equations applied to this form of Fermat's principle result in

∂ L ∂ x k − d d σ ∂ L ∂ x ˙ k = 0 {\displaystyle {\frac {\partial L}{\partial x_{k}}}-{\frac {d}{d\sigma }}{\frac {\partial L}{\partial {\dot {x}}_{k}}}=0}{\ frac {\ partial L} {\ partial x_ {k} }} - {\ frac {d} {d \ sigma}} {\ frac {\ partial L} {\ partial {\ dot x} _ {k}}} = 0

with k=1,2,3 and where L is the optical Lagrangian. Also in this case the optical momentum is defined as

p k = ∂ L ∂ x ˙ k {\displaystyle p_{k}={\frac {\partial L}{\partial {\dot {x}}_{k}}}}p_ {k} = {\ frac {\ partial L} {\ partial {\ dot x } _ {k}}}

and the Hamiltonian P is defined by the expression given abovefor N=3 corresponding to functions x 1 ( σ) {\displaystyle x_{1}\left(\sigma \right)}x_ {1} \ left (\ sigma \ right) , x 2 ( σ) {\displaystyle x_{2}\left(\sigma \right)}x_ {2} \ left (\ sigma \ right) and x 3 ( σ) {\displaystyle x_{3}\left(\sigma \right)}x_ {3} \ left (\ sigma \ right) to be determined

P = x ˙ 1 p 1 + x ˙ 2 p 2 + x ˙ 3 p 3 − L {\displaystyle P={\dot {x}}_{1}p_{1}+{\dot {x}}_{2}p_{2}+{\dot {x}}_{3}p_{3}-L}P = {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + {\ dot {x}} _ {3} p_ {3} -L

And the corresponding Hamilton's equations with k=1,2,3 applied optics are

∂ H ∂ x k = − p ˙ k, ∂ H ∂ p k = x ˙ k {\displaystyle {\frac {\partial H}{\partial x_{k}}}=-{\dot {p}}_{k}\,,\quad {\frac {\partial H}{\partial p_{k}}}={\dot {x}}_{k}}{\ frac {\ partial H} {\ partial x_ {k}}} = - {\ dot {p}} _ {k} \,, \ quad {\ frac {\ partial H} {\ partial p_ {k}}} = {\ dot {x}} _ {k}

with x ˙ k = d x k / d σ {\displaystyle {\dot {x}}_{k}=dx_{k}/d\sigma }{\ dot {x} } _ {k} = dx_ {k} / d \ sigma and p ˙ k = d p k / d σ {\displaystyle {\dot {p}}_{k}= dp_{k}/d\sigma }{\ dot {p}} _ {k} = dp_ {k} / d \ sigma .

The optical Lagrangian is given by

L = n d s d σ = n ( x 1, x 2, x 3) x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 = L ( x 1, x 2, x 3, x ˙ 1, x ˙ 2, x ˙ 3) {\displaystyle L=n{\frac {ds}{d\sigma }}=n\left(x_{1},x_{2},x_{3}\right){\sqrt {{\dot {x}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {x}}_{3}^{2}}}=L\left(x_{1},x_{2},x_{3},{\dot {x}}_{1},{\dot {x}}_{2},{\dot {x}}_{3}\right)}L = n {\ frac {ds} {d \ sigma}} = n \ left (x_ {1}, x_ {2}, x_ {3} \ right) {\ sqrt {{\ точка {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}} = L \ слева (x_ {1}, x_ {2}, x_ {3}, {\ dot {x}} _ {1}, {\ dot {x}} _ {2}, {\ dot {x}} _ { 3} \ right)

and does not explicitly depend on parameter σ. For that reason not all solutions of the Euler-Lagrange equations will be possible light rays, since their derivation assumed an explicit dependence of L on σ which does not happen in optics.

The optical momentum components can be obtained from

p k = n x ˙ k x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 = n d x k d x 1 2 + d x 2 2 + d x 3 2 = n d x k d s {\displaystyle p_{k}=n{\frac {{\dot {x}}_{k}}{\sqrt {{\dot {x}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {x}}_{3}^{2}}}}=n{\frac {dx_{k}}{\sqrt {dx_{1}^{2}+dx_{2}^{2}+dx_{3}^{2}}}}=n{\frac {dx_{k}}{ds}}}p_ {k} = n {\ frac {{\ dot {x}}} _ {k}} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ { 3} ^ {2}}}}} = n {\ frac {dx_ {k}} {{\ sqrt {dx_ {1} ^ {2} + dx_ {2} ^ {2} + dx_ {3} ^ { 2}}}}} = n {\ fr ac {dx_ {k}} {ds}}

where x ˙ k = d x k / d σ {\displaystyle {\dot {x}}_{k}=dx_{k}/d\sigma }{\ dot {x} } _ {k} = dx_ {k} / d \ sigma . The expression for the Lagrangian can be rewritten as

L = n x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 = x ˙ 1 n x ˙ 1 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 + x ˙ 2 n x ˙ 2 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 + x ˙ 3 n x ˙ 3 x ˙ 1 2 + x ˙ 2 2 + x ˙ 3 2 {\displaystyle L=n{\sqrt {{\dot {x}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {x}}_{3}^{2}}}={\dot {x}}_{1}{\frac {n{\dot {x}}_{1}}{\sqrt {{\dot {x}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {x}}_{3}^{2}}}}+{\dot {x}}_{2}{\frac {n{\dot {x}}_{2}}{\sqrt {{\dot {x}}_{1}^{2}+{\dot {x}}_{2}^{2}+{\dot {x}}_{3}^{2}}}}+{\dot {x}}_{3}{\frac {n {\ dot {x}} _ {3}} {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}}L = n {\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2 }}} = {\ dot {x}} _ {1} {\ frac {n {\ dot {x}} _ {1}} {{\ sqrt {{\ dot {x}} _ {1} ^ { 2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}} + {\ dot {x}} _ {2} {\ frac {n {\ dot {x}} _ {2}} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}} + {\ dot {x}} _ {3} {\ frac {n {\ dot {x}} _ {3 }} {{\ sqrt {{\ dot {x}} _ {1} ^ {2} + {\ dot {x}} _ {2} ^ {2} + {\ dot {x}} _ {3} ^ {2}}}}}
= x ˙ 1 p 1 + x ˙ 2 p 2 + x ˙ 3 p 3 {\ displaystyle = {\ dot {x} } _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + {\ dot {x}} _ {3} p_ {3}}= {\ dot {x}} _ {1} p_ {1} + {\ dot {x}} _ {2} p_ {2} + {\ dot {x}} _ {3} p_ {3}

Сравнение этого выражения для L с этим для гамильтониана P можно сделать вывод, что

P = 0 {\ displaystyle P = 0}P = 0

Из выражений для компонентов pk {\ displaystyle p_ {k}}p_ {k} результатов оптического импульса

p 1 2 + p 2 2 + p 3 2 - n 2 (x 1, x 2, x 3) = 0 {\ displaystyle p_ {1} ^ {2} + p_ {2 } ^ {2} + p_ {3} ^ {2} -n ^ {2} \ left (x_ {1}, x_ {2}, x_ {3} \ right) = 0}p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2} -n ^ {2} \ left (x_ {1}, x_ {2}, x_ {3} \ right) = 0

Оптический гамильтониан выбрано как

P = p 1 2 + p 2 2 + p 3 2 - n 2 (x 1, x 2, x 3) = 0 {\ displaystyle P = p_ {1} ^ {2} + p_ {2 } ^ {2} + p_ {3} ^ {2} -n ^ {2} \ left (x_ {1}, x_ {2}, x_ {3} \ right) = 0}P = p_ {1} ^ {2} + p_ {2} ^ {2} + p_ {3} ^ {2} -n ^ {2} \ left (x_ {1}, x_ {2}, x_ {3} \ right) = 0

хотя другие варианты могут быть сделано. Уравнения Гамильтона с k = 1,2,3, определенные выше, вместе с P = 0 {\ displaystyle P = 0}P = 0 определяют возможные световые лучи.

Обобщенные координаты

Как и в гамильтоновой механике, также можно записать уравнения гамильтоновой оптики в терминах обобщенных координат (Q 1 (σ), Q 2 (σ), Q 3 (σ)) {\ Displaystyle \ влево (q_ {1} \ влево (\ сигма \ вправо), q_ {2} \ влево (\ сигма \ вправо), q_ {3} \ left (\ sigma \ right) \ right)}{\ displaystyle \ left (q_ {1} \ left (\ sigma \ right), q_ {2} \ left (\ sigma \ right), q_ {3} \ left (\ sigma \ right) \ right)} , обобщенные импульсы (u 1 (σ), u 2 (σ), u 3 (σ)) {\ displaystyle \ left (u_ {1} \ left (\ sigma \ right), u_ {2} \ left (\ sigma \ right), u_ {3} \ left (\ sigma \ right) \ right)}{\ displaystyle \ left (u_ {1} \ left (\ sigma \ right), u_ {2} \ left (\ sigma \ right), u_ {3} \ left (\ sigma \ right) \ right)} и гамильтониан P как

dq 1 d σ = ∂ P ∂ u 1 du 1 d σ = - ∂ P ∂ q 1 {\ displaystyle {\ frac {dq_ {1}} {d \ sigma}} = { \ frac {\ partial P} {\ partial u_ {1}}} \ quad \ quad {\ frac {du_ {1}} {d \ sigma}} = - {\ frac {\ partial P} {\ partial q_ { 1}}}}{\ frac {dq_ {1}} {d \ sigma}} = {\ frac {\ partial P} {\ partial u_ {1}}} \ quad \ quad {\ frac {du_ {1}} {d \ sigma }} = - {\ frac {\ partial P} {\ partial q_ {1}}}
dq 2 d σ = ∂ P ∂ u 2 du 2 d σ = - ∂ P ∂ q 2 {\ displaystyle {\ frac {dq_ {2}} {d \ sigma}} = {\ frac {\ partial P} {\ partial u_ {2}}} \ quad \ quad {\ frac {du_ {2}} {d \ sigma}} = - {\ frac {\ partial P} {\ partial q_ {2 }}}}{\ frac {dq_ {2 }} {d \ sigma}} = {\ frac {\ частичный P} {\ partial u_ {2}}} \ quad \ quad {\ frac {du_ {2}} {d \ sigma}} = - {\ frac {\ partial P} {\ partial q_ {2}}}
dq 3 d σ = ∂ P ∂ u 3 du 3 d σ знак равно - ∂ п ∂ q 3 {\ displaystyle {\ frac {dq_ {3}} {d \ sigma}} = {\ frac {\ partial P} {\ partial u_ {3}}} \ quad \ quad {\ frac {du_ {3}} {d \ sigma}} = - {\ frac {\ partial P} {\ partial q_ {3}}}}{\ frac {dq_ {3}} {d \ sigma}} = {\ frac {\ partial P} {\ partial u_ {3}}} \ quad \ quad {\ frac {d u_ {3}} {d \ sigma}} = - {\ frac {\ partial P} {\ partial q_ {3}}}
P = p ⋅ p - n 2 = 0 {\ displaystyle P = \ mathbf {p} \ cdot \ mathbf {p} -n ^ {2} = 0}P = {\ mathbf {p}} \ cdot {\ mathbf {p}} - n ^ {2 } = 0

, где оптический импульс определяется как

p = u 1 ∇ q 1 + u 2 ∇ q 2 + U 3 ∇ Q 3 {\ Displaystyle \ mathbf {p} = u_ {1} \ nabla q_ {1} + u_ {2} \ nabla q_ {2} + u_ {3} \ nabla q_ {3}}{\ mathbf {p}} = u_ {1} \ nabla q_ {1} + u_ {2} \ nabla q_ {2 } + u_ {3} \ nabla q_ {3}
= u 1 ‖ ∇ q 1 ‖ ∇ q 1 ‖ ∇ q 1 ‖ + u 2 ‖ ∇ q 2 ‖ ∇ q 2 ‖ ∇ q 2 ‖ + u 3 ‖ ∇ q 3 ‖ ∇ q 3 ‖ ∇ q 3 ‖ { \ Displaystyle = u_ {1} \ | \ nabla q_ {1} \ | {\ frac {\ nabla q_ {1}} {\ | \ nabla q_ {1} \ |}} + u_ {2} \ | \ nabla q_ {2} \ | {\ frac {\ nabla q_ {2}} {\ | \ nabla q_ {2} \ |}} + u_ {3} \ | \ nabla q_ {3} \ | {\ frac {\ набла q_ {3}} {\ | \ набла q_ {3} \ |}}}= u_ {1} \ | \ nabla q_ {1} \ | {\ frac {\ nabla q_ {1}} {\ | \ nabla q_ {1 } \ |}} + u_ {2} \ | \ nabla q_ {2} \ | {\ frac {\ nabla q_ {2}} {\ | \ nabla q_ {2} \ |}} + u_ {3} \ | \ nabla q_ {3} \ | {\ frac {\ nabla q_ {3}} {\ | \ nabla q_ {3} \ |}}
= u 1 a 1 e ^ 1 + u 2 a 2 e ^ 2 + u 3 a 3 e ^ 3 {\ displaystyle = u_ {1} a_ {1} \ mathbf {\ hat {e}} _ {1} + u_ {2} a_ {2} \ mathbf {\ hat {e}} _ {2} + u_ {3} a_ {3} \ mathbf {\ hat {e}} _ {3}}= u_ {1} a_ {1} {\ mathbf {{\ hat {e}}}} _ {1} + u_ {2} a_ {2} {\ mathbf {{\ hat {e}}}} _ {2} + u_ {3} a_ {3} {\ mathbf {{\ hat {e}}} } _ {3}

и e ^ 1 {\ displaystyle \ mathbf {\ hat {e}} _ {1}}{\ mathbf {{\ hat {e}}}} _ {1} , e ^ 2 {\ Displaystyle \ mathbf {\ hat {e}} _ {2}}{\ mathbf {{ \ Hat {e}}}} _ {2} и e ^ 3 {\ displaystyle \ mathbf {\ hat {e}} _ {3}}{\ mathbf {{\ hat {e}}}} _ {3} равны единичные векторы. Частный случай получается, когда эти векторы образуют ортонормированный базис , то есть все они перпендикулярны друг другу. В этом случае ukak / n {\ displaystyle u_ {k} a_ {k} / n}{\ displaystyle u_ {k} a_ {k} / n} - это косинус угла оптического импульса p {\ displaystyle \ mathbf {p }}\ mathbf {p} преобразовывает в единичный вектор e ^ k {\ displaystyle \ mathbf {\ hat {e}} _ {k}}{\ mathbf {{\ hat {e}}}} _ {k} .

См. Также
Ссылки
Последняя правка сделана 2021-05-22 11:54:26
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте