А, чтобы другим читателям было понятно о чем у нас идет речь, я просто приведу данные по нескольким вычислительным экспериментам из этой статьи, где явно видно, что действие на прямом пути не является минимальным, т.е. ПНД не соблюдается. Пусть
Ну что ж, давайте и в самом деле проведём парочку численных экспериментов.
Я так и не смог вычленить из вашего поста полное описание данных для расчётов, поэтому мои данные могут не в точности соответствовать вашим расчётам, но если вы скажете, что на что поменять - я поменяю.
Итак, рассматривается движение пробного заряда в центральном кулоновском поле неподвижного заряда в начале координат. За единицу заряда принимаем величину неподвижного заряда (
); за единицу массы - массу пробного заряда (
); за единицу длины - величину проекции радиус-вектора начального положения заряда на ось
(
); единица времени подбирается так, чтобы коэффициент в законе Кулона был
, она получается равной
.
Исходные данные. Величина пробного заряда
. Начальное положение
, начальная скорость
. Время в пути
.
Вначале считаем истинную траекторию.
(Код Mathematica)
Код:
q = -1
T = 7.32
s = NDSolve[{
x''[t] == q*x[t]/(x[t]^2 + y[t]^2)^(3/2),
y''[t] == q*y[t]/(x[t]^2 + y[t]^2)^(3/2),
x[0] == -5,
y[0] == 1,
x'[0] == 0.6667,
y'[0] == -0.4
},
{x, y},
{t, 0, T}
]
xt[t_] = First[x[t] /. s]
yt[t_] = First[y[t] /. s]
ParametricPlot[{xt[t], yt[t]}, {t, 0, T}, Frame -> True, ImageSize -> Medium]
xend = xt[T]
yend = yt[T]
Вот что получается:
То есть довольно похоже на траекторию на вашем скриншоте. Конечная точка траектории имеет координаты
.
Теперь давайте проверим, действительно ли этой траектории соответсвует минимум действия. Для этого подсчитаем действие для этой траектории и для нескольких других, начинающихся и кончающихся в тех же точках, но отличающихся от истинной. Я взял несколько разных синусоидальных отклонений с амплитудой
(Код Mathematica)
Код:
d1[t_] := 0.1 * Sin[t * Pi/ T]
d2[t_] := 0.1 * Sin[-t * Pi/ T]
d3[t_] := 0.1 * Sin[2*t * Pi/ T]
d4[t_] := 0.1 * Sin[3*t * Pi/ T]
d5[t_] := 0.1 * Sin[-3*t * Pi/ T]
и получил
(Код Mathematica)
Код:
S[x_, y_] := NIntegrate[(x'[t]^2 + y'[t]^2) / 2 - q/Sqrt[x[t]^2 + y[t]^2], {t, 0, T}]
Sdelta[x_, y_, xdelta_, ydelta_] := (
xsumm[t_] = x[t] + xdelta[t];
ysumm[t_] = y[t] + ydelta[t];
S[xsumm, ysumm]
)
St = S[xt, yt]
St1 = Sdelta[xt, yt, d1, d1]
St2 = Sdelta[xt, yt, d2, d2]
St3 = Sdelta[xt, yt, d3, d4]
St4 = Sdelta[xt, yt, d5, d3]
следующие значения для действия:
для истинной трактории
, для близких к истинной -
,
,
и
. Мы видим, что
для истинной траектории действие минимально.
Теперь проделаем то же самое не для истинной траектории, а для равномерного движения по прямой из начальной точки в конечную.
(Код Mathematica)
Код:
xn[t_] = xt[0] + ((xt[T] - xt[0]) / T) * t
yn[t_] = yt[0] + ((yt[T] - yt[0]) / T) * t
Sn = S[xn, yn]
Sn1 = Sdelta[xn, yn, d1, d1]
Sn2 = Sdelta[xn, yn, d2, d2]
Sn3 = Sdelta[xn, yn, d3, d4]
Sn4 = Sdelta[xn, yn, d5, d3]
Мы получим следующие значения действия:
для движения по прямой
, а для вариаций около неё
,
,
и
.
Видим, что
для прямой траектории действие не минимально (
), как и должно быть в соответствии с принципом наименьшего действия, ведь движение по прямой не является истинной траекторией в данной задаче.
Таким образом, численный эсперимент подтверждает принцип наименьшего действия, а вам,
ser, надо искать у себя ошибку.