[论文解读] Nonparametric estimation of mark's distribution of an exponential Shot-noise process

Paul Ilhe, Éric Moulines+2 more|arXiv (Cornell University)|Jun 26, 2015

Stochastic processes and statistical mechanics参考文献 28被引用 10

一句话总结

本文提出了一种非参数估计器，用于核谱学中伽马光子能量分布的估计，利用间接的散粒噪声观测的特征函数。通过建立观测信号特征函数与标记分布特征函数之间的新型反演公式，该方法即使在重叠光子信号使单个能量难以分辨的堆积条件下，也能实现以多项式速率的均匀收敛。

ABSTRACT

In this paper, we consider a nonlinear inverse problem occurring in nuclear science. Gamma rays randomly hit a semiconductor detector which produces an impulse response of electric current. Because the sampling period of the measured current is larger than the mean inter arrival time of photons, the impulse responses associated to different gamma rays can overlap: this phenomenon is known as pileup. In this work, it is assumed that the impulse response is an exponentially decaying function. We propose a novel method to infer the distribution of gamma photon energies from the indirect measurements obtained from the detector. This technique is based on a formula linking the characteristic function of the photon density to a function involving the characteristic function and its derivative of the observations. We establish that our estimator converges to the mark density in uniform norm at a logarithmic rate. A limited Monte-Carlo experiment is provided to support our findings.

研究动机与目标

• 解决从核谱学中间接、噪声大且重叠的测量中估计伽马光子能量概率密度函数（PDF）的逆问题。
• 克服堆积问题带来的挑战，即多个光子到达导致半导体探测器中脉冲响应重叠，使标准的一对一能量估计不可行。
• 开发一种非参数方法，不假设标记分布的参数形式，从而在真实世界核数据中实现灵活推断。
• 在现实采样和信号条件下，建立所提估计器的理论收敛速率。

提出的方法

• 该方法利用观测到的散粒噪声过程的特征函数，该函数源自 Lévy-Khintchine 表示，将观测数据与潜在的标记分布联系起来。
• 推导出一个关键的反演公式，将标记密度的特征函数表示为观测信号特征函数及其导数的函数。
• 通过用独立同分布观测的理论特征函数替换其经验对应物，构建估计器。
• 理论分析依赖于观测过程特征函数衰减的界，这些界在标记密度的矩和支撑条件之下推导得出。
• 利用函数不等式和有界密度卷积的下界，在一致范数下建立收敛速率。
• 该方法避免了对无穷 Lévy 测度的直接估计，而是通过特征函数反演专注于标记密度。

实验结果

研究问题

• RQ1 当探测器信号是指数衰减脉冲的叠加时，能否为伽马光子能量分布构建一个非参数估计器？
• RQ2 在堆积条件下，如何从观测到的散粒噪声过程的特征函数中恢复标记密度的特征函数？
• RQ3 所提非参数估计器在一致范数下的收敛速率是多少？
• RQ4 在标记密度具有有界支撑和矩条件等现实假设下，该方法表现如何？

主要发现

• 所提非参数估计器在一致范数下以多项式速率收敛至真实的标记密度，具体上界由样本大小和底层密度的光滑性决定。
• 在标记密度在紧区间上远离零的假设下，结合额外的矩和正则性条件，建立了收敛速率。
• 观测信号的特征函数以 (1 + |u|)^{-λ/α} 的速率衰减，这为收敛证明提供了关键的下界。
• 该方法在无需估计无穷 Lévy 测度的情况下实现了统一一致性，从而避开了经典 Lévy 过程推断中的主要障碍。
• 蒙特卡洛实验支持了理论发现，展示了该估计器在有限样本设置下的可行性和稳定性。
• 理论结果基于指数脉冲响应的假设，这在实际的 HPGe 探测器中很常见，增强了该方法的实际相关性。