前回から引き続き $K$ が正方形の場合を考えていきましょう。今回の関数のクラスは，自乗可積分、絶対可積分、微分連続を満足する場合です。フーリエ変換が使えるわけですね。

\begin{align} \mathscr{F}[k\ast f]=\mathscr{F}k\mathscr{F}f \end{align}

$\mathscr{F}k$ で割ってから逆フーリエ変換を実行すれば

\begin{align} \mathscr{F}^{-1}\bigg[\dfrac{\mathscr{F}[k\ast f]}{\mathscr{F}k}\bigg]=f \end{align}

となりそうですね。逆フーリエ変換の中身の分母が零となってしまう場合には困ってしまいますが、分母 $\mathscr{F}k$ の零点の測度が零であれば、左辺はルベーグ積分として計算できます。

今、$k$ は正方形の定義関数なので、そのフーリエ像は以前計算した通り $\mathrm{sinc}$ 関数の積になります。零点の測度が零であることが確認できますので、上の計算式は正しく計算できます。これで反転公式が導出できましたね。

ただし、$\mathscr{F}k$ の零点に近い所では非常に小さい値で $\mathscr{F}[k\ast f]$ を割らなければなりません。理論上は $\mathscr{F}[k\ast f]$ の値も同様に小さいのですが、実用上ではあらゆる理由により $\mathscr{F}[k\ast f]$ は小さくなってくれないかもしれません。結果として $\mathscr{F}k$ の零点近傍では理論からの誤差が大きく増幅されてしまうリスクがあります。

ポンペイユ変換は畳込みと違い回転運動も許可されています。$R$ を回転行列として $k(Rx)$ のフーリエ像を計算してみましょう。

\begin{align} \int_{\mathbb{R}^2}k(Rx)e^{-2\pi i\langle x,y\rangle}\,dS=\int_{\mathbb{R}^2}k(Rx)e^{-2\pi i\langle Rx,Ry\rangle}\,dS=\mathscr{F}k(Ry) \end{align}

フーリエ像も同様にして回転させればよいことが分かります。零点だった周波数も、回転させれば零点から外すことができます。零点の近傍では、このように回転させて零点を回避してから計算することで、誤差増幅のリスク回避が期待できます。

今回はここまで。次回は $f$ が微分連続な場合を考えます。(written by J.Nishigata)

研究4- フーリエ変換を用いたポンペイユ変換の反転公式