[論文レビュー] The asymptotic and numerical inversion of the Marcum $Q-$function

Amparo Gil, Javier Segura+1 more|arXiv (Cornell University)|Apr 1, 2014

Statistical Distribution Estimation and Applications参考文献 10被引用数 7

ひとこと要約

本稿では、大規模な $\mu$ に対して一般化された Marcum $Q$-関数およびその補完関数 $P$-関数を逆算するための、新しい漸近的かつ数値的手法を提示する。漸近的展開と単調性・凸性の性質を活用することで、ニュートン法やセカント法の収束を信頼できるものにしている。主な貢献は、数値計算により検証された、効率的かつ高精度な逆算アルゴリズムであり、特にレーダーおよび統計的応用において顕著に効果的である。

ABSTRACT

The generalized Marcum functions appear in problems of technical and scientific areas such as, for example, radar detection and communications. In mathematical statistics and probability theory these functions are called the noncentral gamma or the noncentral chi-squared cumulative distribution functions. In this paper we describe a new asymptotic method for inverting the generalized Marcum $Q-$function and for the complementary Marcum $P-$function. Also, we show how monotonicity and convexity properties of these functions can be used to find initial values for reliable Newton or secant methods to invert the function. We present details of numerical computations that show the reliability of the asymptotic approximations.

研究の動機と目的

• 大規模な $\mu$ に対して一般化された Marcum $Q$-関数およびその補完関数 $P$-関数を効率的かつ信頼性高く逆算するための手法を開発すること。
• 単調性および凸性の性質と漸近的展開を活用し、反復的逆算手法のための堅牢な初期推定値を生成すること。
• パrameter $\mu$, $x$, $y$ の広い範囲にわたって、Marcum 関数の逆数を数値的に安定して計算するアルゴリズムを提供すること。
• Marcum 関数と補完誤差関数との間の関係を確立し、漸近的近似を改善すること。

提案手法

• 鞍点法を用いた経路積分表現を用いて、$Q_\mu(\mu x, \mu y)$ の漸近的展開を導出し、補完誤差関数を含む表現に帰着させる。
• 勾配降下路を虚軸に写像するために変換 $s^2 = 2(\phi(t) - \phi(t_0))$ を用い、$s$ に関する級数展開を可能にする。
• 関数 $g(s) = (f(s) - f(\zeta))/(s - \zeta)$ のテイラー展開から得られる係数 $d_n(\zeta)$ を用いて、剰余項の級数展開を構築する。
• 補完誤差関数 $\mathrm{erfc}(\zeta\sqrt{\mu/2})$ を $Q_\mu(\mu x, \mu y)$ の一次近似として用いる。
• $Q_\mu(x,y)$ の $y$ に関する単調性および凸性に基づいた初期値を用いて、ニュートン法またはセカント法を適用し、高速かつ信頼性の高い収束を保証する。
• $\zeta$ が小さい場合には級数展開による数値的チェックを実施し、$\zeta$ が大きい場合には方程式 (7.12) からの導関数情報を利用したニュートン反復に切り替える。

実験結果

リサーチクエスチョン

• RQ1 大規模な $\mu$ に対して、漸近的近似を用いて一般化された Marcum $Q$-関数の逆数をどのように効率的に計算できるか。
• RQ2 $Q\_\mu(x,y)$ の単調性および凸性の性質が、反復的逆算手法の初期推定値を効果的に生成する上で果たす役割は何か。
• RQ3 経路積分表現から出発し、$Q_\mu(\mu x, \mu y)$ の漸近的展開をどのように導出し、補完誤差関数を含む計算可能な級数に変換できるか。
• RQ4 漸近的展開における剰余項の構造は何か。その係数 $d_n(\zeta)$ はどのように計算されるか。
• RQ5 提案手法は、Marcum 関数の既存の数値的逆算技術と比較して、精度および効率の面でどのように優れているか。

主な発見

• 関数 $Q_\mu(\mu x, \mu y)$ の漸近的近似は、鞍点解析から得られる $\zeta$ を用いた $\mathrm{erfc}(\zeta\sqrt{\mu/2})$ に基づき、強力な一次近似項を提供する。
• 漸近的級数における係数 $d_0(\zeta)$ は明示的に計算され、$d_0(\zeta) = \frac{1 - f(0)}{\zeta}$ と表され、$f(0)$ は $s=0$ における $dt/ds$ の導関数から得られ、閉形式の表現が得られる。
• 小さい $\zeta$ に対しては、級数展開 (7.5) がニュートン反復の代替手段として信頼性があり、漸近的領域における安定性を保証する。
• 方程式 (7.1) の関数 $g(y)$ は凸であり、$y = x+1$ で値が 0 になるため、収束が保証されるニュートン法の適用が可能である。
• 数値計算により、漸近的近似の信頼性が確認され、特に本手法が最も効果を発揮する大規模 $\mu$ の領域で顕著である。
• 本手法は漸近的解析と反復的改良を組み合わせており、$\mu \geq 1$, $x > 0$, $y > 0$ のパrameter 空間全域で高い精度と頑健性を達成している。