Как я могу реализовать сферическую функцию ханкеля первого рода с помощью scipy/numpy или sympy? [закрытый]
Я знал, что в scipy нет встроенного sph_hankel1, тогда я хочу знать, как это правильно реализовать?
Дополнительно: просто покажите мне одну правильную реализацию sph_hankel1 либо с использованием Scipy, либо Sympy.
2 ответов:
Что касается реализации его в SymPy, здесь есть руководство по реализации специальных функций . Мы хотели бы вытянуть запросы на хорошо известные специальные функции.
Для численных процедур они реализованы вmpmath . Это выглядит как он просто использует определение напрямую.
Хотя было бы неплохо, если бы существовала существующая процедура для вычисления сферических функций Ханкеля (как и для обычных функций Ханкеля), они являются просто (сложной) линейной комбинацией сферических функций Бесселя первого и второго рода, поэтому их можно легко вычислить из существующих процедур. Поскольку функции Ханкеля сложны и зависят от того, как вы их применяете, может оказаться полезным переписать ваше выражение в терминах функций Бесселя из первый и второй вид, то есть вполне реальные величины, особенно если ваш конечный результат реален.
Comments