多维傅里叶变换性质与计算
题目
问题7.
(a) 证明多维傅里叶变换具有与问题1中相同的性质(参见小节5.2.5)。
(b) 证明如果多维函数 f f f 具有旋转对称性(即对所有正交变换 Q Q Q,有 f(Qx)=f(x) f(Qx) = f(x) f(Qx)=f(x)),那么其傅里叶变换 f^ \hat{f} f^ 也具有旋转对称性(反之亦然)。
注:等价地,f f f 具有旋转对称性当且仅当 f(x) f(x) f(x) 仅依赖于 ∣x∣ |x| ∣x∣。
问题8. 求以下多维函数的傅里叶变换:
(a) f(x)={1∣x∣≤a,0∣x∣>a; f(x) = \begin{cases} 1 & |x| \leq a, \\ 0 & |x| > a; \end{cases} f(x)={10∣x∣≤a,∣x∣>a;
(b) f(x)={a−∣x∣∣x∣≤a,0∣x∣>a; f(x) = \begin{cases} a - |x| & |x| \leq a, \\ 0 & |x| > a; \end{cases} f(x)={a−∣x∣0∣x∣≤a,∣x∣>a;
© f(x)={(a−∣x∣)2∣x∣≤a,0∣x∣>a; f(x) = \begin{cases} (a - |x|)^2 & |x| \leq a, \\ 0 & |x| > a; \end{cases} f(x)={(a−∣x∣)20∣x∣≤a,∣x∣>a;
(d) f(x)={a2−∣x∣2∣x∣≤a,0∣x∣>a; f(x) = \begin{cases} a^2 - |x|^2 & |x| \leq a, \\ 0 & |x| > a; \end{cases} f(x)={a2−∣x∣20∣x∣≤a,∣x∣>a;
(e) f(x)=e−α∣x∣; f(x) = e^{-\alpha |x|}; f(x)=e−α∣x∣;
(f) f(x)=∣x∣e−α∣x∣; f(x) = |x| e^{-\alpha |x|}; f(x)=∣x∣e−α∣x∣;
(g) f(x)=∣x∣2e−α∣x∣. f(x) = |x|^2 e^{-\alpha |x|}. f(x)=∣x∣2e−α∣x∣.
提示:使用问题7(b),观察到我们只需要计算 f^(0,…,0,k) \hat{f}(0, \ldots, 0, k) f^(0,…,0,k)(即频率空间中沿最后一个坐标轴的值),并使用适当的坐标系(如 n=2 n = 2 n=2 时用极坐标,n=3 n = 3 n=3 时用球坐标等)。注意:此问题可针对 n=2 n = 2 n=2、n=3 n = 3 n=3 或任意 n≥2 n \geq 2 n≥2 求解。
解决问题
问题7(a):证明多维傅里叶变换的性质
多维傅里叶变换在 Rn \mathbb{R}^n Rn 上定义为:
f^(ξ)=∫Rnf(x)e−2πix⋅ξdx,其中x⋅ξ=∑j=1nxjξj. \hat{f}(\xi) = \int_{\mathbb{R}^n} f(x) e^{-2\pi i x \cdot \xi} dx, \quad \text{其中} \quad x \cdot \xi = \sum_{j=1}^n x_j \xi_j. f^(ξ)=∫Rnf(x)e−2πix⋅ξdx,其中x⋅ξ=j=1∑nxjξj.
问题1(未提供)假设涉及一维傅里叶变换的标准性质,如线性、位移、调制、缩放、共轭、卷积等。这些性质在多维情况下类似成立,证明方法类似一维,通过积分变换直接验证。以下是关键性质及简要证明:
-
线性:F{af+bg}=aF{f}+bF{g} \mathcal{F}\{a f + b g\} = a \mathcal{F}\{f\} + b \mathcal{F}\{g\} F{af+bg}=aF{f}+bF{g}
证明:
F{af+bg}(ξ)=∫Rn[af(x)+bg(x)]e−2πix⋅ξdx=a∫Rnf(x)e−2πix⋅ξdx+b∫Rng(x)e−2πix⋅ξdx=af^(ξ)+bg^(ξ). \mathcal{F}\{a f + b g\}(\xi) = \int_{\mathbb{R}^n} [a f(x) + b g(x)] e^{-2\pi i x \cdot \xi} dx = a \int_{\mathbb{R}^n} f(x) e^{-2\pi i x \cdot \xi} dx + b \int_{\mathbb{R}^n} g(x) e^{-2\pi i x \cdot \xi} dx = a \hat{f}(\xi) + b \hat{g}(\xi). F{af+bg}(ξ)=∫Rn[af(x)+bg(x)]e−2πix⋅ξdx=a∫Rnf(x)e−2πix⋅ξdx+b∫Rng(x)e−2πix⋅ξdx=af^(ξ)+bg^(ξ). -
位移:若 g(x)=f(x−x0) g(x) = f(x - x_0) g(x)=f(x−x0),则 g^(ξ)=e−2πix0⋅ξf^(ξ) \hat{g}(\xi) = e^{-2\pi i x_0 \cdot \xi} \hat{f}(\xi) g^(ξ)=e−2πix0⋅ξf^(ξ)
证明:
令 y=x−x0 y = x - x_0 y=x−x0,则 dx=dy dx = dy dx=dy,
g^(ξ)=∫Rnf(x−x0)e−2πix⋅ξdx=∫Rnf(y)e−2πi(y+x0)⋅ξdy=e−2πix0⋅ξ∫Rnf(y)e−2πiy⋅ξdy=e−2πix0⋅ξf^(ξ). \hat{g}(\xi) = \int_{\mathbb{R}^n} f(x - x_0) e^{-2\pi i x \cdot \xi} dx = \int_{\mathbb{R}^n} f(y) e^{-2\pi i (y + x_0) \cdot \xi} dy = e^{-2\pi i x_0 \cdot \xi} \int_{\mathbb{R}^n} f(y) e^{-2\pi i y \cdot \xi} dy = e^{-2\pi i x_0 \cdot \xi} \hat{f}(\xi). g^(ξ)=∫Rnf(x−x0)e−2πix⋅ξdx=∫Rnf(y)e−2πi(y+x0)⋅ξdy=e−2πix0⋅ξ∫Rnf(y)e−2πiy⋅ξdy=e−2πix0⋅ξf^(ξ). -
调制:若 g(x)=e2πix⋅ξ0f(x) g(x) = e^{2\pi i x \cdot \xi_0} f(x) g(x)=e2πix⋅ξ0f(x),则 g^(ξ)=f^(ξ−ξ0) \hat{g}(\xi) = \hat{f}(\xi - \xi_0) g^(ξ)=f^(ξ−ξ0)
证明:
g^(ξ)=∫Rne2πix⋅ξ0f(x)e−2πix⋅ξdx=∫Rnf(x)e−2πix⋅(ξ−ξ0)dx=f^(ξ−ξ0). \hat{g}(\xi) = \int_{\mathbb{R}^n} e^{2\pi i x \cdot \xi_0} f(x) e^{-2\pi i x \cdot \xi} dx = \int_{\mathbb{R}^n} f(x) e^{-2\pi i x \cdot (\xi - \xi_0)} dx = \hat{f}(\xi - \xi_0). g^(ξ)=∫Rne2πix⋅ξ0f(x)e−2πix⋅ξdx=∫Rnf(x)e−2πix⋅(ξ−ξ0)dx=f^(ξ−ξ0). -
缩放:若 g(x)=f(ax) g(x) = f(a x) g(x)=f(ax)(a≠0 a \neq 0 a=0 实数),则 g^(ξ)=∣a∣−nf^(a−1ξ) \hat{g}(\xi) = |a|^{-n} \hat{f}(a^{-1} \xi) g^(ξ)=∣a∣−nf^(a−1ξ)
证明:
令 y=ax y = a x y=ax,则 dx=∣a∣−ndy dx = |a|^{-n} dy dx=∣a∣−ndy,
g^(ξ)=∫Rnf(ax)e−2πix⋅ξdx=∫Rnf(y)e−2πi(a−1y)⋅ξ∣a∣−ndy=∣a∣−n∫Rnf(y)e−2πiy⋅(a−1ξ)dy=∣a∣−nf^(a−1ξ). \hat{g}(\xi) = \int_{\mathbb{R}^n} f(a x) e^{-2\pi i x \cdot \xi} dx = \int_{\mathbb{R}^n} f(y) e^{-2\pi i (a^{-1} y) \cdot \xi} |a|^{-n} dy = |a|^{-n} \int_{\mathbb{R}^n} f(y) e^{-2\pi i y \cdot (a^{-1} \xi)} dy = |a|^{-n} \hat{f}(a^{-1} \xi). g^(ξ)=∫Rnf(ax)e−2πix⋅ξdx=∫Rnf(y)e−2πi(a−1y)⋅ξ∣a∣−ndy=∣a∣−n∫Rnf(y)e−2πiy⋅(a−1ξ)dy=∣a∣−nf^(a−1ξ). -
卷积:F{f∗g}=F{f}F{g} \mathcal{F}\{f * g\} = \mathcal{F}\{f\} \mathcal{F}\{g\} F{f∗g}=F{f}F{g},其中 (f∗g)(x)=∫Rnf(y)g(x−y)dy (f * g)(x) = \int_{\mathbb{R}^n} f(y) g(x - y) dy (f∗g)(x)=∫Rnf(y)g(x−y)dy
证明:
F{f∗g}(ξ)=∫Rn(∫Rnf(y)g(x−y)dy)e−2πix⋅ξdx=∫Rnf(y)e−2πiy⋅ξ(∫Rng(z)e−2πiz⋅ξdz)dy=f^(ξ)g^(ξ), \mathcal{F}\{f * g\}(\xi) = \int_{\mathbb{R}^n} \left( \int_{\mathbb{R}^n} f(y) g(x - y) dy \right) e^{-2\pi i x \cdot \xi} dx = \int_{\mathbb{R}^n} f(y) e^{-2\pi i y \cdot \xi} \left( \int_{\mathbb{R}^n} g(z) e^{-2\pi i z \cdot \xi} dz \right) dy = \hat{f}(\xi) \hat{g}(\xi), F{f∗g}(ξ)=∫Rn(∫Rnf(y)g(x−y)dy)e−2πix⋅ξdx=∫Rnf(y)e−2πiy⋅ξ(∫Rng(z)e−2πiz⋅ξdz)dy=f^(ξ)g^(ξ),
其中令 z=x−y z = x - y z=x−y.
这些性质与一维情况类似,证明依赖于积分的基本操作(如变量替换、Fubini 定理)。小节 5.2.5 应提供更多细节,但此处已覆盖主要内容。
问题7(b):证明旋转对称性的等价性
旋转对称性指 f(Qx)=f(x) f(Qx) = f(x) f(Qx)=f(x) 对所有正交变换 Q Q Q(即 QTQ=I Q^T Q = I QTQ=I,detQ=±1\det Q = \pm 1detQ=±1) 成立。等价地,f(x) f(x) f(x) 仅依赖于 ∣x∣ |x| ∣x∣(即径向函数)。需证:f f f 旋转对称当且仅当 f^ \hat{f} f^ 旋转对称。
证明:
设 g(x)=f(Qx) g(x) = f(Qx) g(x)=f(Qx)。其傅里叶变换为:
g^(ξ)=∫Rnf(Qx)e−2πix⋅ξdx. \hat{g}(\xi) = \int_{\mathbb{R}^n} f(Qx) e^{-2\pi i x \cdot \xi} dx. g^(ξ)=∫Rnf(Qx)e−2πix⋅ξdx.
令 y=Qx y = Qx y=Qx,则 x=QTy x = Q^T y x=QTy(因为 Q Q Q 正交,Q−1=QT Q^{-1} = Q^T Q−1=QT),且 dx=dy dx = dy dx=dy(雅可比行列式为 ∣detQ∣=1 |\det Q| = 1 ∣detQ∣=1)。代入得:
g^(ξ)=∫Rnf(y)e−2πi(QTy)⋅ξdy=∫Rnf(y)e−2πiy⋅(Qξ)dy=f^(Qξ), \hat{g}(\xi) = \int_{\mathbb{R}^n} f(y) e^{-2\pi i (Q^T y) \cdot \xi} dy = \int_{\mathbb{R}^n} f(y) e^{-2\pi i y \cdot (Q \xi)} dy = \hat{f}(Q \xi), g^(ξ)=∫Rnf(y)e−2πi(QTy