Не могли бы вы рассказать мне быстрый и точный метод вычисления BesselK(mu,z), BesselI(mu,z), где mu — действительное число?
Алгоритм вычисления модифицированной функции Бесселя первого рода
Ответы (2)
Существует несколько способов представления функций Бесселя, и ни один из них не является простым математически или в закрытой форме. Тем не менее, эта статья, по-видимому, представляет собой недавний подход к их эффективной оценке:
http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.120.6055
Кажется, это лучший способ, чем метод старой школы 1950-х годов:
http://www.ams.org/journals/mcom/1959-13-066/S0025-5718-1959-0105794-5/S0025-5718-1959-0105794-5.pdf
У меня не было времени, чтобы всесторонне и всесторонне протестировать ее самостоятельно, кроме того, что я увидел, что она работает, но библиотека Colt для Java включает функцию Бесселя целочисленного порядка в пакете cern.jet.math.
https://dst.lbl.gov/ACSSoftware/colt/
Я также очень поддерживаю цифровую библиотеку математических функций: