
Teorija_Polja
.pdfÒåìà 6
Постоянное магнитное поле
Будем считать, что заряды совершают финитное движение в ограниченной области пространства; будем также считать, что мы интересуемся только усредненными по времени характеристиками системы. В рамках такой модели можно считать, что постоянные токи порождают не зависящее от времени магнитное поле . Опуская значки усреднения, запишем соответствующие уравнения (с. 74):
div H = 0, rot H = |
4π |
j |
|
|
|
|
|||
|
|
|
|
|
|||||
|
|
|
c |
|
|
|
|
||
Вводя векторный потенциал A: H = rot A, получим |
|
|
|
||||||
rot [rotA] ≡ grad divA − A = |
4π |
j |
|
|
|
|
4π |
j. (6.1) |
|
c |
|
A = |
|||||||
|
|
||||||||
divA = 0 (кулоновская калибровка) |
|
− c |
|||||||
|
|
|
|
|
|
|
|
|
Для каждой компоненты вектора A получается уравнение Пуассона ви-
да (5.1), решение которого уже было найдено (с. 86). По аналогии с (5.3) мы сразу можем записать:
|
1 |
Z |
j(r ) |
|
|||
A(r) = |
|
|
|
0 |
dV 0. |
(6.2) |
|
|
c |
|r − r0| |
В случае токов, порожденных движением дискретных зарядов, имеем:
j(r) = |
|
eαvα δ(r − rα) = A(r) = |
1 |
|
|
|
eαvα |
|
(6.3) |
|||||
α |
|
|
|
α |
|
|
|
|
|
|
||||
|
c |
|
| |
r |
− |
rα |
| |
|||||||
|
X |
|
|
|
|
|
X |
|
|
|
||||
|
|
~ |
~ |
|
|
|
|
|
|
|
~ |
|
|
|
С помощью тождества: rot(a · b) ≡ a · rot(b) + [grad a × b], ãäå a ñêà-
лярная, ~
b векторная функция, находим выражение для напряженности
91

Тема 6. Постоянное магнитное поле
магнитного поля, порожденного токами:
H(r) = rotr |
c |
Z |
| |
rj(r0r0 |
| |
dV 0 |
= c |
Z |
[j(r0)r× (rr0− r0)] dV 0 |
|||
|
1 |
|
|
|
) |
|
|
|
1 |
|
|
|
закон Био-Савара.
92
(6.4)
Также, как и в случае электростатики, можно строить мультипольное разложение, полагая, что точка наблюдения удалена на расстояние, гораздо большее размеров системы. При
этом следует считать, что все макроскопические величины являются результатом усреднения по времени, гораздо большему, чем характерный период финитного движения зарядов (если движение строго периодич- но, можно усреднять по периоду).
Запишем мультипольное разложение для векторного потенциала A, ограничиваясь первым и вторым (дипольным) членом:
1 |
|
|
|
|
|
|
e v |
|
(1) |
1 |
|
|
e v |
α |
|
1 |
|
|
|
|
|
rα |
1 |
|
|
||||||||||||||||||
A(r) = |
|
|
|
|
|
|
|
|
|
|
α α |
|
' |
|
|
|
|
|
α |
|
− |
|
α |
eαvα |
· r |
|
|
= |
|||||||||||||||
c |
|
|
α |
| |
r |
|
rα |
|
c |
|
α |
|
r |
| |
|
c |
r |
|
|||||||||||||||||||||||||
|
|
|
|
X |
|
|
− |
|
| |
|
|
|
X |
| |
|
|
|
|
|
X |
|
|
|
|
|
|
|
| |
| |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
" |
|
|
|
|
|
|
|
|
|
# |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 1 |
|
|
|
eαrα |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
(2) |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
eαvα (rα · r) |
= |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
dt |
|
c |
|
|
|
|
r |
|
c |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
α |
|
|
r 3 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
|
|
|
| |
| |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X | |
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
(2) |
|
d |
|
1 |
|
|
|
|
e r |
|
1 |
|
|
|
e v (r |
|
r) |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
= |
|
|
" |
|
|
|
|
α |
α |
+ |
|
|
|
|
|
α α |
|
α · |
|
|
# + |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
dt |
c |
|
α |
r |
| |
2c |
|
α |
|
r |
|
3 |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
| |
|
|
|
|
|
X |
|
|
| | |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
1 |
|
|
|
|
|
eαvα {vα (rα · r) − rα (vα · r)} |
|
|
(3) |
[m × r] |
|
|
|
|
|
||||||||||||||||||||||||||
+ |
|
|
|
|
|
|
|
= |
|
. |
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
2c |
|
|
|
|
|
|
|
|
|
|
|
|
|
r 3 |
|
|
|
|
|
|
|
|
|
|
|
|
r 3 |
|
|
|
||||||||||||
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
| | |
|
|
|
|
|
|
|
|
|
|
|
|
|
| | |
|
|
|
|
|
|
||||
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь сделаны следующие преобразования:
1.выполнено мультипольное разложение с сохранением дипольных членов;
2.в этом равенстве сделано преобразование вида uv˙ = (uv˙ ) − uv˙ ;
3.здесь были отброшены полные производные по времени, поскольку при финитном движении полные производные обращаются в нуль при усреднении по времени.

Тема 6. Постоянное магнитное поле |
93 |
Наконец, окончательный результат выражен через магнитный момент
m:
1 |
X |
eα [rα × vα]. |
|
|||
m = |
|
(6.5) |
||||
2 c |
|
α |
||||
Итак, в дипольном приближении |
|
|
|
|
||
|
A = |
[m × r] |
. |
(6.6) |
||
|
|
|
|
|r|3 |
|
Если у всех частиц отношение заряда к массе одинаково, то, полагая
в нерелятивистском пределе (v/c 1) p = mv, можно получить простую формулу, связывающую магнитный момент m и механический момент M (момент импульса):
1 |
X |
eα [rα × vα] = |
e |
|
X |
m [rα × vα] = |
|||||
m = |
|
|
|
||||||||
2 c |
α |
2 mc |
|
α |
|||||||
|
|
= |
|
e |
X |
[rα × pα] = |
e |
|
M. |
||
|
|
|
|
|
|
||||||
|
|
2 mc |
α |
2 mc |
Как и в случае электростатического поля, можно рассмотреть силовое взаимодействие системы зарядов с магнитным полем и выполнить мультипольное разложение.
Существенным фактором, определяющим движение, является |
ìî- |
||||||||||||||
ìåíò ñèë K, действующий на магнитный диполь: |
|
|
|||||||||||||
|
|
1 |
X |
|
|
|
|
1 |
X |
|
|
||||
FЛоренц = |
|
|
eα [vα × H] ≡ |
|
|
fα |
|
||||||||
c |
α |
c |
α |
|
|||||||||||
|
X |
|
|
|
|
1 |
X |
|
|
|
|
|
|
||
K = |
[rα × fα] = |
|
|
eα [rα × [vα × H]] . |
(6.7) |
||||||||||
α |
c |
|
α |
Можно произвести преобразования, при которых из полученного для K выражения будут выделены члены, представляющие собой полные
производные по времени; при усреднении по времени они обращаются в нуль, в результате чего остается простая формула:
K = [m × H]. |
(6.8) |
Запишем уравнение движения магнитного момента, помещенного в магнитное поле. Для этого следует выразить магнитный момент через

Тема 6. Постоянное магнитное поле |
94 |
механический и применить уравнения механики для механического момента M:
˙ |
˙ |
|
M = K m˙ = γ M = γ K = γ [m × H] |
|
|
m˙ |
= γ [m × H], |
(6.9) |
e
2 mc гиромагнитное отношение.
Решая получившееся уравнение (6.9), можно показать, что оно описывает ларморовскую прецессию вектора m движение, при кото-
ром конец вектора m описывает окружность, а сам вектор заметает поверхность конуса, ось которого направлена вдоль магнитного поля H.
Òåìà 7
Волновые решения. Поле движущихся зарядов
7.1Электромагнитные волны в свободном пространстве
Важный тип решений уравнений Максвелла ненулевые решения, в отсутствие зарядов и токов, которые описывают волны электромагнитного поля. Можно считать, что волновые решения соответствуют самоподдерживающемуся процессу изменение со временем поля одного типа служит причиной появления другого поля и наоборот.
Запишем уравнения Максвелла в отсутствие зарядов и токов:
rot
rot
|
|
|
1 |
∂H |
|
|
||||
E = − |
|
|
|
|
|
, |
div H = 0, |
|||
c |
|
∂t |
||||||||
H = |
1 |
∂E |
, |
|
div E = 0. |
|||||
c |
|
|
|
|||||||
|
|
∂t |
|
|
Применяя формулу преобразования двойного векторного произведения, имеем:
hi
rot rot A ≡ r × [r × A] |
= r (r · A) − (r · r)A = grad (div A) − A, |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 ∂ |
|
|
1 ∂2 |
|
||||||||||||
rot rot E = grad ( |
div E) − |
E = − |
|
|
|
|
|
|
rot H = − |
|
|
|
|
|
E, |
|||||||||||||
c |
∂t |
c2 |
∂t2 |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|( |
|
{z } |
− H = c ∂t |
E = −c2 ∂t2 H |
|
|||||||||||||||||||||||
|
|
|
=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 ∂2 |
|
||||||||
rot rot H = grad div H) |
|
|
1 ∂ |
|
, |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
rot |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
| |
|
|
{z |
|
} |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
95
Тема 7. Волновые решения. Поле движущихся зарядов |
96 |
|||||||||
откуда следуют волновые уравнения: |
|
|||||||||
E − |
1 ∂2 |
H − |
1 ∂2 |
|
||||||
|
|
|
E = 0, |
|
|
|
H = 0. |
(7.1) |
||
c2 |
∂t2 |
c2 |
∂t2 |
Воспользуемся известным математическим результатом: с помощью метода Д'Аламбера можно представить решение однородного волнового уравнения в виде суперпозиции встречных бегущих волн:
∂2f |
= c |
2 ∂2f |
= f(x, t) = Φ(x ± ct). |
(7.2) |
|
∂t2 |
|
∂x2 |
Роль функции Φ может играть произвольная дважды дифференцируе-
мая функция. Если подставить Φ в виде разложения в ряд (интеграл)
Фурье, можно получить решение волнового уравнения в виде суперпозиции гармонических (синусоидальных) бегущих волн. В физике важны именно такие решения, т.к. часто синусоидально-осциллирующий ток порождает именно синусоидальную волну. Кроме того, изучив свойства одиночной синусоидальной волны, мы можем потом собрать из таких волн произвольное решение. Синусоидальную волну можно записать как
f(x, t) = f0 sin(ωt ± kx + ϕ0) ω = k · c.
Можно ввести обобщение на трехмерный случай:
Re |
|
i (ωt k r) |
|
|
|
|
2 |
|
n i (ωt k r)o |
||
|
|
|
E = −k2 · E n Re ei (ωt−k·r) , |
||||||||
E = E n cos(ωt − k · r) = |
|
||||||||||
= E n e |
− · |
o |
E¨ = |
− |
ω |
|
· |
E n Re e |
− · . |
||
|
n |
|
|
|
|
|
n |
o |
¨
Подставляя E è E в волновое уравнение, мы убеждаемся, что равен-
ство имеет место при
c2k2 = ω2.
Мы видим, что если выбрать волновой вектор k произвольным об-
разом и положить ω = ±|k| c, то записанное выражение для E всегда
будет решением волнового уравнения.
Множество пространственных точек, которым отвечает постоянное значение фазы волны, называется волновым фронтом. В данном приме-
ре уравнение волнового фронта имеет вид: k·r = const. Нетрудно видеть,

Тема 7. Волновые решения. Поле движущихся зарядов |
97 |
что здесь записано уравнение плоскости, перпендикулярной волновому вектору k. Поэтому такое решение называют плоской волной.
Аналогичным образом можно записать решение волнового уравнения для магнитного поля: H = H q Re ei (ωt−k·r) .
Согласуем параметры полей:
rot E = |
1 |
∂H |
|
|
rot E ≡ [r × E] = −i [k × n] E Re nei (ωt k·r)o , |
||||||
|
−c |
∂t |
|
1 ∂H = |
1 i ω H q Re ei (ωt k·r) . |
||||||
|
|
|
|
|
|
|
|
|
|
n |
o |
|
|
|
|
|
|
|
|
−c ∂t |
−c |
Приравнивая эти выражения, имеем: −[k ×n] E = −1c ω H q = |k|H q,
èëè: |
k × n E = ±q H. |
|
||||
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
Учитывая, что |n| = |q| = 1, видим, что |
|
|||||
q = |
n × k |
, n = ± q × k |
, E = H. |
|||
|
|
k |
k |
|
Таким образом, в свободной электромагнитной волне:
1) направления E è H взаимно-перпендикулярны, а также перпенди-
кулярны волновому вектору k, указывающему направление распростра-
нения волны (поперечная волна); 2) амплитуды электрической и магнитной составляющей равны:
|E| = |H|.
7.2Поле движущегося заряда
Комбинируя уравнения Максвелла, мы получили выше волновые уравнения непосредственно для напряженностей полей E и H. Однако
уравнения такого же типа можно записать и для потенциалов, выбрав лоренцеву калибровку (3.18). Имеем
|
|
|
∂Ak |
∂Ai |
|
∂F ik |
|
|
|
4π |
|
|
∂Ak |
|
|
|
|||||||||||
F ik ≡ |
|
|
|
|
− |
|
|
|
, |
|
|
|
= − |
|
|
ji, |
|
|
|
= 0 = |
|
||||||
|
∂xi |
|
∂xk |
∂xk |
c |
∂xk |
|
||||||||||||||||||||
∂F ik |
∂ |
|
∂Ak |
|
∂Ai |
|
|
|
∂ |
|
|
∂Ak |
|
|
|
∂2Ai |
|
π |
|
||||||||
|
= |
|
|
|
|
|
− |
|
= |
|
|
|
|
|
− |
|
= − |
4 |
ji. |
||||||||
∂xk |
∂xk |
|
∂xi |
|
∂xk |
∂xi |
∂xk |
∂xk∂xk |
c |

Тема 7. Волновые решения. Поле движущихся зарядов |
|
|
|
|
98 |
|||||||||||||
Таким образом, |
|
Ai = 4π ji |
|
|
A − c12 ∂t2 |
= − c |
j, |
|||||||||||
|
∂2Ai |
|
|
|||||||||||||||
|
∂xk∂xk |
≡ |
|
c |
|
|
|
|
|
|
∂2A |
|
|
4π |
||||
|
|
|
|
1 ∂2ϕ |
|
|
|
|
(7.3) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ |
− c2 ∂t2 |
= |
− |
4π ρ. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если поле порождается точечными зарядами (как частный случай одним зарядом), то ρ(r, t) = e(t) δ(r), т.е. распределение зарядовой плот-
ности сингулярно. В этом случае решения волнового уравнения будут иметь вид решения для пустого пространства, дополненного так, чтобы учесть наличие зарядов.
При наличии единственного заряда в точке r = 0 имеется сферическая
симметрия. Симметричное решение удобно искать, вводя сферические координаты:
y = r sin θ sin ϕ |
|
|
f = 1 ∂ |
r2 ∂f + |
1 |
|
∂ |
sin θ∂f |
+ |
1 |
|
|
∂2f |
|||||||||||||
x = r sin θ cos ϕ |
|
r2 ∂r |
|
∂r |
r2 sin θ ∂θ |
|
∂θ |
|
r2 sin2 |
θ ∂ϕ2 |
||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
z = r cos θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=0, |
(не зависит от углов) |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
| |
|
|
|
|
|
{z |
|
|
} |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
если функция f сферически симметрична
В итоге волновое уравнение для сферически-симметричного потенциала ϕ примет вид:
|
1 |
|
∂2ϕ |
|
1 |
|
∂ |
|
|
∂ϕ |
|
|
|
|
1 |
|
∂2ϕ |
|
||||||
ϕ − |
|
|
|
= |
|
|
|
r2 |
|
|
− |
|
|
|
|
= 0 |
||||||||
c2 |
∂t2 |
r2 |
∂r |
∂r |
c2 |
∂t2 |
||||||||||||||||||
|
|
|
|
χ(r, t) |
|
|
|
|
∂2χ |
|
1 ∂2χ |
|
|
|
|
|||||||||
ϕ = |
|
|
|
|
|
|
|
|
− |
|
|
|
|
= 0 |
|
|||||||||
|
r |
|
∂r2 |
c2 |
∂t2 |
|
Простая замена переменных позволила нам свести уравнение к однородному волновому. Однако, поскольку роль пространственной переменной
теперь играет радиальная координата r, уравнение описывает сфери-
ческие волны; максимумы и минимумы таких волн лежат на концентрических сферах. Формальное решение Д'Аламбера имеет вид:
χ(r, t) = χ1(t − r/c) + χ2(t + r/c)
и описывает суперпозицию двух волн: волны, убегающей на бесконеч- ность (запаздывающее решение) и волны, приходящей из бесконечности

Тема 7. Волновые решения. Поле движущихся зарядов |
99 |
и сходящейся к началу координат (опережающее решение). Поскольку принципу причинности удовлетворяет только запаздывающее решение, имеем:
χ = χ(t − r/c).
Заметим, что при r → 0 поведение членов, входящих в волновое уравнение для ϕ, таково: ϕ r−3 ∂2ϕ/∂t2 r−1. Смысл этого
неравенства прост: в малой окрестности точки r = 0, где находится то-
чечный заряд, эффектами запаздывания, связанными с конечностью скорости света, можно пренебречь . Формально это сводится к отбра- сыванию члена ∂2ϕ/∂t2; волновое уравнение в этом случае превраща-
ется в уравнение Пуассона (5.1), решение которого нам известно это кулоновский потенциал: ϕ 1/r. Таким образом, χ(t) = e(t), è
ϕ(r, t) = e(t − r/c). r
Можно сказать, что найден обобщенный закон Кулона, учитывающий
конечность скорости распространения электромагнитного взаимодействия.
Далее, исходя из выражений (4.6) для плотности системы зарядов, а также для потенциала системы зарядов (5.3, 6.2), нетрудно угадать
интегральное соотношение, связывающее ρ è ϕ. Сразу можно записать и
аналогичное соотношение для j è A (в нем нужно добавить множитель
1/c):
ϕ(r, t) = |
Z |
ρ(r0, t − |r − r0|/c) |
dV 0 |
+ ϕ |
, |
||||||
|
Z |
|
|r − r0| |
|
0 |
|
|||||
A(r, t) = |
|
1 |
j(r0 |
, t − |r − r0|/c) |
dV 0 |
+ A . |
|||||
|
c |
|
|
|
|||||||
|
|
|
|r − r0| |
|
|
0 |
(7.4)
(7.5)

Тема 7. Волновые решения. Поле движущихся зарядов |
100 |
7.3Поле точечного заряда, движущегося по заданной траектории
Пусть заряд движется по траектории, которая задана в параметри- ческой форме уравнением:
r0(t) = i · x0(t) + j · y0(t) + k · z0(t)
(параметром служит время t). В этом случае, если величина заряда фик-
сирована, плотность распределения заряда может быть задана формулой:
ρ(r, t) = e δ(r − r0(t)).
Подставим это в выражение (7.4), предварительно введя интегрирование по времени и δ-функцию:
|
ρ(r0 |
, t − 1c |r − r0|) |
|
ρ(r0, τ) |
1 |
|
|||||||||||||||
ϕ(r, t) = Z |
|
|
|
|
|
|
dV 0 = ZZ |
|
|
|
|
|
|
δ τ − t + |
|
|
|r − r0| dV 0dτ |
||||
|
|
|
|r − r0| |
|
|r − r0| |
c |
|||||||||||||||
|
|
|
|
e δ |
(r0 |
|
r (τ)) |
|
1 |
|
|
|
|
|
|||||||
|
= |
|
|
− 0 |
δ τ − t + |
|
|
|r − r0| dV 0dτ = |
|||||||||||||
|
|
|
|
r |
r0 |
c |
|||||||||||||||
|
|
ZZ |
|
| |
− |
| |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
||||
|
|
= e Z |
|
|
|
δ τ − t + |
|
|r − r0(τ)| dτ. |
|
||||||||||||
|
|
|
|r − r0(τ)| |
c |
|
Далее учтем, что
Z |
f(τ) δ g(τ) dτ = |
f(τ ) |
|
, ãäå g(τ0) = 0. |
|||||
|g0(τ )| |
|||||||||
В нашем случае: |
g0(τ) = ∂τ τ − t + c|r − r0(τ)| |
= |
|||||||
|
|||||||||
|
|
∂ |
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
1 ∂ q
=1 + c ∂τ (x − x0(τ))2 + (y − y0(τ))2 + (z − z0(τ))2 = 1 −2(x − x0) x˙ 0 − 2(y − y0) y˙0 − 2(z − z0) z˙0
=1 + c 2q(x − x0(τ))2 + (y − y0(τ))2 + (z − z0(τ))2 =
= 1 − |
1 |
|
r − r0(τ) · v0(τ) |
. |
|
|
|||
|
c |
|r − r0(τ)| |