std::cyl_bessel_i, std::cyl_bessel_if, std::cyl_bessel_il

来自cppreference.com
 
 
 
 
在标头 <cmath> 定义
(1)
float       cyl_bessel_i ( float nu, float x );

double      cyl_bessel_i ( double nu, double x );

long double cyl_bessel_i ( long double nu, long double x );
(C++17 起)
(C++23 前)
/* 浮点类型 */ cyl_bessel_i( /* 浮点类型 */ nu, /* 浮点类型 */ x );
(C++23 起)
float       cyl_bessel_if( float nu, float x );
(2) (C++17 起)
long double cyl_bessel_il( long double nu, long double x );
(3) (C++17 起)
在标头 <cmath> 定义
template< class Arithmetic1, class Arithmetic2 >
/* 公共浮点类型 */ cyl_bessel_i( Arithmetic1 nu, Arithmetic2 x );
(A) (C++17 起)
1-3) 计算 nux常规修正贝塞尔函数标准库提供所有以无 cv 限定的浮点类型作为参数 nux 的类型的 std::cyl_bessel_i 重载。 (C++23 起)
A) 为算术类型的所有其他组合提供额外重载。

参数

nu - 函数的阶数
x - 函数的参数

返回值

如果没有发生错误,那么返回 nux 的常规修正柱贝塞尔函数,即 I
nu
(x) = Σ
k=0
(x/2)nu+2k
k!Γ(nu+k+1)
(对于 x≥0)。

错误处理

可能报告 math_errhandling 中指定的错误。

  • 如果参数是 NaN,那么返回 NaN 且不报告定义域错误。
  • 如果 nu>=128,那么行为由实现定义。

注解

不支持 C++17,但支持 ISO 29124:2010 的实现在定义了 __STDCPP_MATH_SPEC_FUNCS__ 为至少 201003L 的值,且用户在包含任何标准库头文件前定义了 __STDCPP_WANT_MATH_SPEC_FUNCS__ 时也会提供此函数。

不支持 ISO 29124:2010 但支持 TR 19768:2007 (TR1) 的实现也会在标头 tr1/cmath 及命名空间 std::tr1 中提供此函数。

此函数的一种实现参考 boost.math

额外重载不需要以 (A) 的形式提供。它们只需要能够对它们的第一个实参 num1 和第二个实参 num2 满足以下要求:

  • 如果 num1num2 具有 long double 类型,那么 std::cyl_bessel_i(num1, num2)std::cyl_bessel_i(static_cast<long double>(num1),
                      static_cast<long double>(num2))
    的效果相同。
  • 否则,如果 num1 和/或 num2 具有 double 或整数类型,那么 std::cyl_bessel_i(num1, num2)std::cyl_bessel_i(static_cast<double>(num1),
                      static_cast<double>(num2))
    的效果相同。
  • 否则,如果 num1num2 具有 float 类型,那么 std::cyl_bessel_i(num1, num2)std::cyl_bessel_i(static_cast<float>(num1),
                      static_cast<float>(num2))
    的效果相同。
(C++23 前)

如果 num1num2 具有算术类型,那么 std::cyl_bessel_i(num1, num2)std::cyl_bessel_i(static_cast</* 公共浮点类型 */>(num1),
                  static_cast</* 公共浮点类型 */>(num2))
的效果相同,其中 /* 公共浮点类型 */num1num2 的类型中浮点转换等级浮点转换子等级最高的浮点类型,整数类型的实参被视为具有与 double 相等的浮点转换等级。

如果不存在等级和子等级最高的浮点类型,那么在重载决议时不会从提供的重载中产生可用的候选。

(C++23 起)

示例

#include <cmath>
#include <iostream>
 
int main()
{
    // 对 nu == 0 点检查
    double x = 1.2345;
    std::cout << "I_0(" << x << ") = " << std::cyl_bessel_i(0, x) << '\n';
 
    // 对 I_0 的级数展开 
    double fct = 1;
    double sum = 0;
    for (int k = 0; k < 5; fct *= ++k)
    {
        sum += std::pow(x / 2, 2 * k) / std::pow(fct, 2);
        std::cout << "sum = " << sum << '\n';
    }        
}

输出:

I_0(1.2345) = 1.41886
sum = 1
sum = 1.381
sum = 1.41729
sum = 1.41882
sum = 1.41886

外部链接

Weisstein, Eric W. “第一类修正贝塞尔函数”来自 MathWorld--A Wolfram Web Resource。

参阅

(第一类)圆柱贝塞尔函数
(函数)