说明:该分析的基础是《Fundamentals of Computerized Tomography Image Reconstruction from Projections 2ndEd》书中(8.10)式及其上下文。
函数$\phi(u)$的希尔伯特变换$[\mathcal{H}\phi](v)$为
\begin{equation}[\mathcal{H}\phi](v)=-\frac{1}{\pi} \bigg(\int\limits_{-\infty}^{\infty}\frac{\phi(u)}{v-u}du\bigg)\end{equation}
其实就是$\phi$与$\rho$的卷积,即$\mathcal{H}\phi=\phi*\rho$,其中
\begin{equation}\rho(u)=-\Big(\frac{1}{\pi u}\Big)\end{equation}
称为希尔伯特卷积核函数。通过数值计算求(1)式的值较为困难,一种解决方法就是正则化卷积核函数,即用一族函数$\{\rho_A|A>0\}$来近似$\rho$,其中$A$是正整数,该正则化函数的目的是满足:$\phi$在任意合理的实数点$v$上都有
\begin{equation}\lim\limits_{A\to\infty}[\phi*\rho_A](v)=[\mathcal{H}\phi](v)\end{equation}
这就是希尔伯特变换的卷积近似表达式。现在的问题是如何确定$\rho_A$,由于很难找到可以直接逼近$\rho$的近似函数,所以从其傅里叶变换入手,寻找可行的近似表达。根据
\begin{equation}\mathcal{F}\big(\frac{1}{u}\big)(U)=-i\pi \mathrm{sgn}(U)\end{equation}