[论文解读] The generalized Marcum $Q-$function: an orthogonal polynomial approach

Szilárd András, Árpád Baricz+1 more|arXiv (Cornell University)|Oct 16, 2010

Advanced Wireless Communication Techniques参考文献 53被引用 85

一句话总结

本文提出了一种针对任意正实数阶 ν 的广义马尔库姆 Q-函数的新型幂级数表示，采用广义拉盖尔多项式。该方法确保了绝对收敛性，并提供了收敛速度的截断误差分析，为现有计算方法提供了一种数值稳定且高效的替代方案，尤其适用于标准实现无法处理的非整数 ν 值的情况。

ABSTRACT

A novel power series representation of the generalized Marcum $Q-$function of positive order involving generalized Laguerre polynomials is presented. The absolute convergence of the proposed power series expansion is showed, together with a convergence speed analysis by means of truncation error. A brief review of related studies and some numerical results are also provided.

研究动机与目标

• 开发一种针对任意正实数阶 ν 的广义马尔库姆 Q-函数的数值稳定且高效的级数表示。
• 克服现有算法中 ν 仅限于正整数的局限，特别是在非中心卡方分布和衰落信道等应用中。
• 提供收敛性分析及截断误差界，以支持实际实现。
• 通过正交多项式展开，将 Pent 的一阶结果推广至一般实数阶马尔库姆 Q-函数。
• 为标准软件（如 MATLAB 的 marcumq）在非整数 ν 情况下失效时，提供一种可行的马尔库姆 Q-函数计算替代方案。

提出的方法

• 利用广义拉盖尔多项式 L_n^(ν-1)(a) 作为基函数，推导广义马尔库姆 Q-函数的新级数展开。
• 通过拉盖尔多项式生成函数恒等式，对修正贝塞尔函数 I_{ν-1}(at) 进行变换，建立该表示。
• 应用马尔库姆 Q-函数的积分表示，并在绝对收敛条件下交换求和与积分顺序。
• 使用生成函数恒等式：∑_{n≥0} [L_n^(α)(x)/L_n^(α)(0)] * (-1)^n z^n / n! = Γ(α+1) e^{-z} (xz)^{-α/2} I_α(2√(xz))，其中 α = ν-1。
• 推导出关键级数形式：Q_ν(a,b) = e^{-a} ∑_{n≥0} (-1)^n L_n^(ν-1)(a) / Γ(ν+n) * ∫_b^∞ z^{n+ν-1} e^{-z} dz，得到表达式 (2.7)。
• 通过与 MATLAB 的 marcumq 函数比较结果，对整数 ν 进行数值验证，并提供 marcumq 不可用的 ν = 7.7 时的数值结果。

实验结果

研究问题

• RQ1 能否利用广义拉盖尔多项式，为任意实数阶 ν > 0 的广义马尔库姆 Q-函数构造幂级数表示？
• RQ2 该级数的收敛行为如何依赖于参数 a、b 和 ν？截断误差行为如何？
• RQ3 与现有方法相比，该级数在数值稳定性与效率方面是否更优，尤其是在非整数 ν 的情况下？
• RQ4 与标准级数形式 (2.4) 相比，该表示在结构与计算实用性方面有何差异？
• RQ5 当标准软件（如 ν = 7.7 时）失效时，该方法能否可靠计算马尔库姆 Q-函数？

主要发现

• 所提出的级数表示 (2.7) 对所有 a, ν > 0 和 b ≥ 0 均绝对收敛。
• 该方法可准确计算非整数 ν（如 ν = 7.7）下的广义马尔库姆 Q-函数，此时 MATLAB 的 marcumq 函数不可用。
• 数值结果表明，对于 ν = 1, 3, 5，该方法与 MATLAB 的 marcumq 在多种 a 和 b 取值下完全一致。
• 截断误差分析提供了收敛速度的定量度量，支持带有误差界的实际实现。
• 该方法允许对固定 a 值预先计算拉盖尔多项式项，从而在多个 b 值上实现快速评估，显著提升计算效率。
• 由于在不同变量上展开，该级数表示与标准形式 (2.4) 显著不同，为分析与计算提供了全新的路径。