Laplacian 算子的 Fourier 变换

在一些优化模型中含有 偏导项与数据项，对其求最小化。一种方法是利用 Eular Lagrange 变换，求解出关于未知变量的解。 还有一种方法是利用 快速Fourier 变换来求解。
如L0 滤波中，对子问题 $S$ 的解，公式 (7)

最小化目标函数： $\sum_{p}\{(S_p-I_p)^2+\beta((\partial_xS_p-h_p)^2+(\partial_yS_p-v_p)^2) \}$ （1）

上式的解用FFT 表示为：
$S=\mathcal{F}^{-1}\{\frac{\mathcal{F(I)}+\beta(\mathcal{F(\partial_x)^*\mathcal{F(h)+\mathcal{F}(\partial_x)^*\mathcal{F}(v))}})}{\mathcal{F(1)+\beta(\mathcal{F}(\partial_x)^*\mathcal{F}(\partial_x))+\mathcal{F}(\partial_y)^*\mathcal{F}(\partial_y)}} \}$

这个Fourier 等式是如何推导出来的 ？

推导过程中主要涉及两个主要内容
（1） Eular lagrange 变换。
（2）Laplacian 算子的 FFT 表达。

首先以 1 维函数为例推导 Laplacian 算子的 FFT 表达：
记函数$f(x)$的 FFT 为$\mathcal{F}(u)$, 有

$f(x)\Leftrightarrow \mathcal{F}(u)$,

平移函数的FFT,

$f(x-n_0)\Leftrightarrow exp(-j2\pi un_0)\mathcal{F}(u)$

则前向差分算子 $\nabla_xf=f(x+1)-f(x)$ 的FFT为：

$f(x+1)-f(x)\Leftrightarrow(exp(j2\pi u)-1)\mathcal{F}(u)$

后向差分算子的FFT为

$f(x)-f(x-1)\Leftrightarrow(1-exp(-j2\pi u))\mathcal{F}(u)$

则 Laplacian 算子 $\triangle^2f(x)=f(x+1)-2f(x)+f(x-1)$ 的FFT 变换为

$\triangle^2f(x)\Leftrightarrow (exp(j2\pi u)-2+exp(-j2\pi u))\mathcal{F}(u)$
$=-((exp(j2\pi u)-1)\dot(exp(-j2\pi u)-1)\mathcal{F}(u)$
$=-((exp(j2\pi u)-1)\dot(exp(j2\pi u)-1)^*\mathcal{F}(u)$

这时得到了Laplacian 算子与梯度算子的 FFT 关系
$\mathcal{F}(\triangle^2f(x))\Leftrightarrow-\mathcal{F}(\nabla_xf)*\mathcal{F}(\nabla_xf)^*\mathcal{F}(u)$
其中$^*$表示复共轭。

下面看Eular_Lagrange 等式
（1）的Eular _Lagrange 等式表达为

$2(S_p-I_p)-2\beta((\partial_x(\partial_xS_p-h_p)+\partial_y(\partial_yS_p-v_p))=0$

即：
$S-\beta(\partial_{xx}S+\partial_{yy}S)=I-\beta(\partial_xh+\partial_yv)$

对等上式两边进行FFT有

$\mathcal{F}(S)+\beta（\mathcal{F}(\partial_x) \mathcal{F}(\partial_x)^*+\mathcal{F}(\partial_y)\mathcal{F}(\partial_y)^* )\mathcal{F}(S)=\mathcal{F}(I)-\beta(\mathcal{F}(\partial_x)\mathcal{F}(h)+\mathcal{F}(\partial_y)\mathcal{F}(v))$

整理得
$\mathcal{F}(S)=\frac{\mathcal{F}(I)-\beta(\mathcal{F}(\partial_x)\mathcal{F}(h)+\mathcal{F}(\partial_x)\mathcal{F}(v))}{1+\beta（\mathcal{F}(\partial_x) \mathcal{F}(\partial_x)^*+\mathcal{F}(\partial_y)\mathcal{F}(\partial_y)^* )}$