c++ - c++ 中 "Matlab besselk(x,y,1)"的确切等价物是什么？

Question

我试过 boost::math::cyl_bessel_k(x,y) * exp(y)。在大多数情况下，这等于 Matlab 的缩放 besselk(x,y,1)。但在某些情况下(例如，x=1、y=2000)当 besselk(x,y)=0 和 boost::math::cyl_bessel_k(x,y)=0，Matlab 的缩放版本 besselk(x,y,1) 给我的不同值在 10^-3 左右变化。但是 boost::math::cyl_bessel_k(x,y) * exp(y) 返回 -nan。我想为 Matlab 的 besselk(x,y,1) 找到一个等效语句。我该如何处理？

Answer 1

我在 Boost 中没有看到任何满足您需要的东西(尽管您可以使用较低级别的函数自己实现它)。正如您所发现的，缩放贝塞尔函数并不是简单地通过乘以 exp(z) 来计算的。 GSL似乎已包含此功能，例如 gsl_sf_bessel_Knu_scaled .对于“完全等效”，您可以查看 the paper和 code by Amos ，例如 CBESK . Matlab 和 Octave 似乎都使用了这个实现。但是，代码是用 Fortran 语言编写的，因此您需要对其进行翻译或对其进行包装(this project 似乎已经这样做了，所以它可能会有用 - 还有其他的)。

您还可以使用 Matlab 的 Coder 和 codegen也输出一些东西。