Тогда функция вида (31) допускает такую интерпретацию. Пусть фиксированному начальному моменту времени t = t0 отвечает неко- торая поверхность уровня
ϕ0(х, y, z) = ωt0.
С течением времени эта поверхность как бы перемещается в про- странстве, занимая последовательно положения с различными, уве- личивающимися значениями ϕi (x, y, z) = ω ti. Такую бегущую по- верхность уровня принято называть волновым фронтом. Если в на- чальный момент времени t = t0 на поверхности волнового фронта за- фиксировать некоторую точку с координатами x0, y0, z0, то при дви- жении фронта волны эта точка опишет некоторую кривую, называе- мую лучом. Радиус-вектор
r(t) = xi + yj + zk
текущей точки на луче является функцией времени, и нашей бли- жайшей целью будет составить дифференциальные уравнения для определения этой функции. Очевидно, что скорость точки на луче
v(t) = dr = dx i + dy j + dz k dt dt dt dt
должна во все время движения оставаться коллинеарной вектору градиента функции ϕ (х, y, z), вычисленного в данной точке. Это по- ложение можно воспринимать и как определение луча.
Положим g = grad ϕ.
Тогда
dr = λg, dt
или в координатной форме
|
dx |
= λ |
∂ϕ , |
||
|
|
|
|
||
|
dt |
∂x |
|||
dy |
= λ |
∂ϕ , |
|||
|
|
||||
|
dt |
∂y |
|||
|
|
|
|
|
(34) |
dz |
= λ |
∂ϕ . |
|||
|
|||||
dt |
∂z |
||||
Коэффициент пропорциональности λ (х, y, z) находится следую- щим образом. Из данных выше определений волнового фронта и лу- ча вытекает тождество
21
ϕ ( x(t), y(t), z(t)) = ωt,
дифференцируя обе части которого по времени, получим
∂ϕ dx + ∂ϕ dy + ∂ϕ dz = ω
.
∂x dt ∂y dt ∂z dt
Вносим сюда вместо компонент скорости их выражения из урав- нений луча (34). И, с учетом уравнения эйканала (32), найдем
λ = c2 .
ω
Теперь нетрудно вычислить квадрат элемента длины луча ds2 = dr2 = dx2 + dy2 + dz2 .
Принимая во внимание (32), (34) и (35), имеем ds2 = c2 dt2 .
Соответственно, получим величину лучевой скорости
v = c = ds . dt
(35)
(36)
Вводя в систему уравнения луча (34) в качестве переменной дли- ну дуги s, можно представить ее и в такой форме:
dx |
= |
|
|
c |
|
∂ϕ , |
|
|||
ds |
|
|
|
|
||||||
|
|
|
|
ω ∂x |
|
|||||
dy |
|
= |
c |
∂ϕ , |
(37) |
|||||
|
|
|
|
|||||||
ds |
|
|
|
|
ω ∂y |
|
||||
dz |
= |
c |
∂ϕ . |
|
||||||
|
|
|
||||||||
ds |
|
|
|
|
ω ∂z |
|
||||
Системы уравнений луча (34) или (37) содержат в себе, вообще говоря, некую неопределенность в виде неизвестных функций
∂ϕ , |
∂ϕ , |
∂ϕ , так как относительно функции ϕ (x, y, z) известно лишь |
∂x |
∂y |
∂z |
то, что она должна быть одним из решений уравнения эйканала (32). Но по значению только модуля градиента скалярное поле однозначно не восстанавливается. Однако систему уравнений луча возможно преобразовать и во вполне определенную систему, включающую только заданные функции.
22
Запишем следующие тождества:
|
∂g2 |
|
= 2 |
∂ϕ ∂2ϕ + 2 ∂ϕ |
|
|
∂2ϕ |
+ 2 |
∂ϕ |
∂2ϕ |
, |
|||||||||||||||||||||||||||||||||||||
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
∂x ∂x2 |
|
|
∂y ∂x∂y |
|
|
∂z ∂x∂z |
|
|
||||||||||||||||||||||||||||||||||||||
|
∂g2 |
|
= 2 |
∂ϕ |
|
|
∂2ϕ |
+ 2 |
∂ϕ ∂2ϕ + 2 |
∂ϕ |
∂2ϕ |
, |
||||||||||||||||||||||||||||||||||||
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
∂x ∂x∂y |
|
|
|
∂y ∂y2 |
|
|
∂z ∂y∂z |
|
|
|||||||||||||||||||||||||||||||||||||
|
∂g2 |
|
= 2 |
∂ϕ |
|
|
∂2ϕ |
+ 2 |
∂ϕ |
|
|
∂2ϕ |
|
+ 2 ∂ϕ ∂2ϕ |
, |
|||||||||||||||||||||||||||||||||
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
∂x ∂x∂z |
|
|
|
∂y ∂y∂z |
|
|
∂z ∂z2 |
|
|
|||||||||||||||||||||||||||||||||||||
а также |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
d |
∂ϕ = |
∂2ϕ |
dx |
+ |
|
∂2ϕ |
|
dy |
+ |
|
∂2ϕ |
|
|
|
|
dz |
, |
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
dt ∂x |
∂x2 dt |
∂x∂y dt |
|
∂x∂z dt |
|
|
|||||||||||||||||||||||||||||||||||||||||
|
d |
∂ϕ = |
∂2ϕ |
|
|
dx |
+ |
|
∂2ϕ |
dy |
+ |
|
∂2ϕ |
|
|
dz |
, |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
dt ∂y |
∂x∂y dt |
|
|
∂y2 dt |
|
∂y∂z dt |
|
|
|||||||||||||||||||||||||||||||||||||||
|
d |
∂ϕ = |
∂2ϕ |
|
dx |
+ |
|
∂2ϕ |
|
dy |
+ ∂2ϕ |
dz |
. |
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
dt ∂z |
∂x∂z dt |
|
|
∂y∂z dt |
|
∂z2 dt |
|
|
|||||||||||||||||||||||||||||||||||||||
Отсюда, с учетом уравнений луча (34), |
||||||||||||||||||||||||||||||||||||||||||||||||
|
d |
∂ϕ = |
c2 |
|
∂ |
|
g2 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
dt ∂x |
2ω ∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
d |
∂ϕ = |
c2 |
|
∂ |
g2 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
dt ∂y |
2ω ∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
d |
∂ϕ = |
c2 |
|
∂ |
g2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
dt ∂z |
2ω ∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
или, заменяя здесь квадрат градиента его значением из уравнения эйканала (32),
d |
∂ϕ = −ω |
∂ |
|
ln c, |
|
|
|
|
|
|
|||
dt ∂x |
∂x |
|
||||
d |
∂ϕ = −ω |
∂ |
ln c, |
(38) |
||
|
|
|
||||
dt ∂y |
∂y |
|
||||
d |
∂ϕ = −ω |
∂ |
ln c. |
|
||
|
|
|
||||
dt ∂z |
∂z |
|
||||
Таким образом, вводя в рассмотрение три новые неизвестные функции времени
23
l = |
∂ϕ , m = |
∂ϕ , n = |
∂ϕ |
|
∂x |
∂y |
∂z |
и прибавляя к системе из трех уравнений (34) еще три равенства (38), получим для определения хода луча систему из шести уравнений первого порядка, у которых правые части составлены из заданных функций. А именно
|
|
|
dx |
|
= |
|
c2 |
|
l, |
|
||||||||
|
|
|
dt |
ω |
|
|
|
|
|
|
||||||||
|
|
|
dy |
|
= |
c2 |
|
m, |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
dt |
ω |
|
|
|
|
|
|||||||||
|
|
|
dz |
= |
c2 |
|
n, |
|
||||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
dt |
ω |
|
|
|
|
|
|||||||||
dl |
|
= −ω |
|
∂ |
ln c, |
|
||||||||||||
|
|
|
|
|||||||||||||||
dt |
∂x |
|
|
|
|
|
||||||||||||
dm |
= −ω |
∂ |
ln c, |
(39) |
||||||||||||||
|
|
|
|
|||||||||||||||
dt |
|
|
∂y |
|
||||||||||||||
dn |
= −ω |
|
∂ |
ln c. |
|
|||||||||||||
|
|
|
|
|||||||||||||||
dt |
|
∂z |
|
|||||||||||||||
В качестве текущей лучевой координаты может быть выбран лю- бой удобный параметр. Система (39) приобретет еще более компакт- ную форму, если вместо времени t ввести лучевой параметр σ, опре- делив его дифференциал так:
dσ = c2 dt.
ω
При этом переменная σ имеет размерность площади ([σ] = м2). Теперь система (39) принимает вид
dx = l, dσ
dy = m, dσ
dz = n, dσ
24
dl |
|
= ω2 |
|
∂ |
|
|
|
1 |
|
, |
|||
dσ |
|
|
|
|
|
|
|||||||
|
2 ∂x c2 |
|
|
||||||||||
dm |
= ω2 |
∂ |
|
|
1 |
|
|
, |
|||||
|
|
|
|
|
|
||||||||
dσ |
|
2 ∂y c2 |
|
|
|||||||||
dn |
|
= ω2 |
∂ |
|
1 |
. |
|||||||
dσ |
|
|
|
||||||||||
|
2 ∂z c2 |
|
|
||||||||||
Из этой системы, конечно, можно исключить три неизвестные функции l, m, n и получить систему из трех уравнений второго по- рядка:
d 2 x = ω2 dσ2 2
d 2 y = ω2 dσ2 2
d 2 z = ω2 dσ2 2
∂1
∂x c2 ,
∂1
∂y c2 ,
∂1
∂z c2 .
Если возвратиться к первоначальному аргументу t, означающему время, то получится такая система трех уравнений второго порядка:
d |
2 |
x |
|
= −c |
∂c + |
2 |
|
|
|
dx |
|
|
∂c |
dx |
+ |
∂c |
dy |
+ |
∂c dz |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
, |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
dt2 |
∂x c dt |
∂x dt |
∂y dt |
∂z dt |
|||||||||||||||||||||||||||||||||
d |
2 |
y |
= −c |
∂c + |
2 |
|
|
|
dy |
|
|
∂c |
dx |
+ |
∂c |
dy |
+ |
∂c dz |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
, |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
dt2 |
∂y c dt |
∂x dt |
∂y dt |
∂z dt |
|||||||||||||||||||||||||||||||||
d |
2 |
z |
= −c |
∂c + |
2 |
|
|
|
dz |
|
|
∂c |
dx |
+ |
∂c |
dy |
+ |
∂c dz |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
, |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
dt2 |
∂z c dt |
∂x dt |
∂y dt |
∂z dt |
|||||||||||||||||||||||||||||||||
которую, очень удобно представить, в векторной форме:
d 2 r |
= −c gradc + |
2 dr |
dr |
|
|||||
|
|
|
|
gradc |
|
|
. |
(40) |
|
dt2 |
c |
dt |
dt |
||||||
25