Изучая метрические геодезические Шварцшильда, можно легко прийти к следующему дифференциальному уравнению
Я попытался вычислить первый интеграл с помощью Mathematica, и результат, который я получил, абсурден, что дает сложное собственное время. Вот код, который я использовал:
Integrate[1/Sqrt[C^2 - (1 - 2*G*M/r)], {r, 2*G*M, 2*G*M + h}]
И вывод:
ConditionalExpression[(G M (2 Sqrt[C^2 - C^4] - I Log[2] - I Log[((-I + 2 I C^2 + 2 Sqrt[C^2 - C^4]) G M)/Sqrt[1 - C^2]]))/(1 - C^2)^(3/2) + (-Sqrt[1 - C^2] (h + 2 G M) Sqrt[C^2 - h/(h + 2 G M)] + I G M Log[2] + I G M Log[(-I h + I C^2 h - I G M + 2 I C^2 G M + Sqrt[1 - C^2] h Sqrt[C^2 - h/(h + 2 G M)] + 2 Sqrt[1 - C^2] G M Sqrt[C^2 - h/(h + 2 G M)])/Sqrt[1 - C^2]])/(1 - C^2)^(3/2), ((G M)/h != 0 && Re[(GM)/h] >= 0) ||Re[(G M)/h] < -(1/2) || (G M)/h \[NotElement] Reals]
Что мне здесь не хватает? Спасибо
Я тоже не совсем понимаю, как это "легко увидеть" из того.
Однако вы можете решить интеграл для собственного времени, сначала заметив, что для
у нас (с
)
что следует из включения
в
и рассматривая падающую частицу, которая покоится в
, значение
.
Это вы можете использовать, чтобы написать свое первое уравнение как
Отсюда интеграл по собственному времени будет равен
которые вы можете решить с помощью Mathematica или также вводя параметризацию
с .
результат, который я получаю, абсурден, давая комплексу правильное время.
В вашем результате есть I
, что я предполагаю, это обозначение Mathematica для
, но это не означает, что результат сложный. Для некоторых методов интеграции характерно получение результатов, которые выглядят сложными, но на самом деле оказываются реальными, когда вы их оцениваете.
Вы излишне усложняете свои результаты, принимая и как параметры. При выполнении такого рода интеграла нужно преобразовать как можно больше переменных в безразмерные величины, что обычно приводит к устранению ненужных констант, подобных этим. Здесь вы хотите перейти к переменной . Не вводите его в систему компьютерной алгебры (CAS), не выполнив эту подготовку. CAS глупы и производят сложный вывод, поэтому дайте ему шанс побороться, выполнив ту же настройку, которую вы обычно делаете, работая без CAS.
Я использую CAS maxima с открытым исходным кодом, а не mathematica. У него нет проблем с получением результата без явного в нем.
$ maxima -q --batch-string="assume(h>0 and (c^2-1)<0 and c>0); integrate(-(c^2-(1-1/x))^(-1/2),x,1+h,1);"
(%i1) assume(h > 0 and c^2-1 < 0 and c > 0)
2
(%o1) [h > 0, c < 1, c > 0]
(%i2) integrate(-(c^2-(1-1/r))^((-1)/2),r,1+h,1)
2 2 2
2 sqrt(1 - c ) sqrt(h + 1) sqrt((c - 1) h + c )
(%o2) (sqrt(1 - c ) atan(----------------------------------------------)
2 2
(c - 1) h + c - 1
2 2 2 4 2
+ (c - 1) sqrt(h + 1) sqrt((c - 1) h + c ))/(c - 2 c + 1)
2
2 c sqrt(1 - c ) 3
sqrt(1 - c ) atan(--------------) + c - c
2
c - 1
- ------------------------------------------
4 2
c - 2 c + 1
Возможно, вы захотите попробовать настроить его таким же образом в mathematica. Безусловно, после замены переменных интерпретировать будет проще и легче.
На самом деле нет никакой гарантии, что этот результат окажется реальным только потому, что нет в нем. В нем действительно есть один квадратный корень, который выглядит так, как будто он может быть мнимым. То, как мы знаем, что оно должно быть реальным, состоит в том, что это определенный интеграл от действительного подынтегрального выражения. Если вы хотите проверить, действительно ли это правда, попробуйте просто использовать случайные числа для и . Если результат окажется реальным, то это не совпадение. Если получается сложно, то может иметь место одно из следующего: (1) вы закодировали его неправильно, (2) есть ошибка в CAS или (3) есть проблема с ответвлениями для арктангенса.
Триаттикус
безопасная сфера
Аксионоподобные частицы
Аксионоподобные частицы
ПрофРоб