前言

    在信号处理和计算机视觉领域,我们谈论傅里叶变换时,往往第一反应是看它的幅度谱(Magnitude Spectrum)——它告诉我们信号里有哪些频率成分,能量有多大。然而,真正决定信号“骨架”和图像“几何结构”的,往往是那个经常被当成配角的相位谱(Phase Spectrum)。我们将从最直观的三角函数平移,一路过渡到复指数表示,彻底理清相位的物理本质。

直观视角:从三角函数的“左加右减”说起

    回到我们最熟悉的标准正弦波公式:

                                                     f(t)=Asin(\omega t + \phi)

    这里的 A 为该频率下的振幅(或强度),这里的 \phi 被称为初相(Initial Phase)。单纯看 \phi 这个角度可能有些抽象,为了看清它的物理本质,我们把角频率 \omega 提取出来,改写成如下形式: 

                                                    f(t)=Asin(\omega(t + \frac{\phi}{\omega}))

    根据函数平移的铁律(左加右减),这个波形相对于标准的 $\sin(\omega t)$,在时间轴上向左平移了 $\Delta t = \frac{\phi}{\omega}$ 个单位。所谓的相位,其实就是描述了正弦波左右平移情况。很显然,对于一个确定周期的正弦波,其最大可分辨平移距离为该正弦波的一个周期,如平移1.5个周期与平移0.5个周期看起来是一致的。

    如果信号的周期 $T = 1$,那么角频率 $\omega = \frac{2\pi}{T} = 2\pi$。带入上式:

                                                    f(t)=Asin(2\pi(t + \frac{\phi}{2\pi}))

    此时,相位 $\phi$ 与波形平移的周期数存在着极其纯粹的线性对应关系:

                                                    平移的周期数 = \frac{\phi}{2\pi}

  •    $\phi = 0$:未平移。
  •    $\phi = \frac{\pi}{2}$:向左平移了 \frac{1}{4} 个周期,正弦波变成了余弦波。
  •    $\phi = \pi$:向左平移了 $\frac{1}{2}$ 个周期,波峰和波谷刚好对调。
  •    $\phi = 2\pi$:向左平移了 $1$ 个完整周期,由于周期信号的重复性,波形与原波形完全重合。

    在实数三角函数的世界里,相位本质上就是描述了正弦信号在轴上相对平移了多少个周期。它不决定波的大小,只决定波在哪里发生。因此,在公式 f(t)=Asin(\omega t + \phi) 中,我们紧跟 \phi 这个参数,关注该参数是如何进一步被封装的,从而实现对相位的深度理解。

 数学飞跃:走向复指数表示

    虽然实数形式的正余弦组合很直观,但在处理复杂信号或多维信号时,三角函数的和差化度公式会让人抓狂。于是,数学家请出了欧拉公式:

                                                                 $e^{j\theta} = \cos\theta + j\sin\theta$

   在复指数形式的傅里叶级数中,一个周期为 1 的信号可以优雅地写为:

                                                                 $f(t) = \sum_{n=-\infty}^{\infty} c_n e^{j 2\pi n t}$

    这里的傅里叶系数 $c_n$ 绝大多数情况下是一个复数。我们可以把它写成极坐标形式:

                                                                $c_n = |c_n|e^{j\phi_n}$

    在三角函数中的相位信息 \phi 被封装到了复数 c_{n} 中,通过如下变换,你会发现:

                        $c_n e^{j 2\pi n t} = |c_n|e^{j\phi_n} \cdot e^{j 2\pi n t} = |c_n|e^{j(2\pi n t + \phi_n)}$

  • $|c_n|$(复数的模):就是这个频率分量的振幅

  • $\phi_n$(复数的辐角):就是这个频率分量的初相

复指数的精妙之处在于:它把原本空间域中复杂的“左右平移”转换成了复平面上向量的“角度旋转”。空间域的位移定理在复频域变成了极为简单的乘法:

                                                           $f(t - t_0) \longleftrightarrow c_n e^{-j 2\pi n t_0}$

原本在时域或空间域需要大费周章的平移操作,在频域只需要给每个频点的复数向量乘以一个旋转因子即可。

以上的讲解直接给出了傅里叶级数的复指数形式,这样看起来有些突兀,下面通过一些推导将傅里叶级数与复指数关联起来,并重点关注 \phi 是如何被封装起来的。

我们有如下实数域的波动方程:

                                              f(t)=Asin(\omega t + \phi)

根据欧拉公式(Euler's Formula),复指数与三角函数的关系为:

                                             $e^{jx} = \cos x + j \sin x$

                                             $e^{-jx} = \cos x - j \sin x$

两式相减,可以消去余弦项 $\cos x$,从而得到正弦项的复指数表示(这也叫欧拉公式的正弦变形)

                                             $\sin x = \frac{e^{jx} - e^{-jx}}{2j}$

我们将标准正弦波中的整体相位 $x = \omega t + \phi$ 代入上式形式中:

                                              $f(t) = A \cdot \left[ \frac{e^{j(\omega t + \phi)} - e^{-j(\omega t + \phi)}}{2j} \right]$

展开括号,将系数 $A$ 和分母 $2j$ 分配进去:

                                              $f(t) = \frac{A}{2j} e^{j(\omega t + \phi)} - \frac{A}{2j} e^{-j(\omega t + \phi)}$

利用指数的运算法则 $e^{a+b} = e^a \cdot e^b$,我们将随时间变化的项($\omega t$)与不随时间变化的初相项($\phi$)剥离开来:

                                               $f(t) = \left( \frac{A}{2j} e^{j\phi} \right) e^{j\omega t} - \left( \frac{A}{2j} e^{-j\phi} \right) e^{-j\omega t}$

为了让形式看起来更像标准的傅里叶级数项(即各项之间用“+”号连接),我们微调一下第二项的符号:

                                              $f(t) = \left( \frac{A}{2j} e^{j\phi} \right) e^{j\omega t} + \left( -\frac{A}{2j} e^{-j\phi} \right) e^{-j\omega t}$

此时,复指数形式的雏形已经出现了!它可以写成:

                                              $f(t) = c_1 e^{j\omega t} + c_{-1} e^{-j\omega t}$

目前的系数中还挂着分母上的虚数单位 $j$,这不够直观。在复数中,我们知道 $\frac{1}{j} = -j$,且根据极坐标表示:

  • $-j = e^{-j\frac{\pi}{2}}$ (在复平面上指向正下方 90度)
  • $j = e^{j\frac{\pi}{2}}$ (在复平面上指向正上方 90度)

 对于正频率项的系数 $c_1$

                 $c_1 = \frac{A}{2j} e^{j\phi} = -j \frac{A}{2} e^{j\phi} = \left( \frac{A}{2} e^{-j\frac{\pi}{2}} \right) e^{j\phi} = \frac{A}{2} e^{j\left(\phi - \frac{\pi}{2}\right)}$

 对于负频率项的系数 $c_{-1}$

                 $c_{-1} = -\frac{A}{2j} e^{-j\phi} = j \frac{A}{2} e^{-j\phi} = \left( \frac{A}{2} e^{j\frac{\pi}{2}} \right) e^{-j\phi} = \frac{A}{2} e^{-j\left(\phi - \frac{\pi}{2}\right)}$

 将清洗后的系数带回,我们得到了完美的复指数表达

                 $f(t) = \frac{A}{2} e^{j\left(\phi - \frac{\pi}{2}\right)} e^{j\omega t} + \frac{A}{2} e^{-j\left(\phi - \frac{\pi}{2}\right)} e^{-j\omega t}$   

同时我们关注到正弦波动方程中的初始相位 \phi 被封装到了e^{j(\phi - \frac{\pi}{2})} 和 e^{-j(\phi - \frac{\pi}{2})}两项中,这两项构成了我们之前讨论傅里叶级数中的 c_{n} 项。

完成上述推导后,可给出以下醍醐灌顶的物理直觉解释

  • 为什么振幅变成了 $\frac{A}{2}$?从最终公式可见,原本空间域里高为 $A$ 的波,在频域被拆成了正频率 $\omega$ 和负频率 $-\omega$ 两部分,每个分支的模长都变成了 $\frac{A}{2}$。这在物理上意味着:能量没有变,只是被平均分到了正、负两个旋转方向相反的复空间中。

  • 为什么相位平白无故减去了 $\frac{\pi}{2}$?因为复指数波 $e^{j\omega t}$ 的老家(实部)默认是余弦波$\cos$。我们一开始给出的信号是正弦波($\sin$),根据三角变换,$\sin(x) = \cos(x - \frac{\pi}{2})$。这个 $-\frac{\pi}{2}$ 极其精准地表达了正弦波天然滞后余弦波 $90^\circ$ 的物理事实。

  • 怎么理解“负频率” $-\omega t$?在实数世界里,根本没有“每秒转成负圈”的说法。但通过上面的推导,你会发现 $c_{-1}$ 恰好是 $c_1$ 的共轭复数( $c_{-1} = c_1^*$ )。因此,在处理实数信号的傅里叶变换时,我们总是得到一对幅度相等但旋转角度相反的系数,这里的旋转角度就是我们一直关注的相位。

概念的延续:非周期信号与连续傅里叶变换(FT)

有人会问:傅里叶级数探讨的是周期信号,那对于非周期信号,傅里叶变换(FT)中的相位概念还一致吗?

完全一致。 傅里叶变换可以被视为周期趋于无穷大($T \to \infty$)时的极限情况。此时,离散的频谱变成了连续谱:

                                               $F(s) = |F(s)|e^{j\theta(s)}$

这里的连续相位函数 $\theta(s)$ 依然遵循完全相同的物理定律:

  • 线性相位与时延:如果整个非周期信号在时间轴上平移了 $t_0$,其频域响应将整体乘以 $e^{-j2\pi s t_0},这意味着相位谱满足 $\theta(s) = -2\pi s t_0$。相位随频率 $s$ 呈严格的线性正比变化。
  • 色散现象(非线性相位):如果 $\theta(s)$ 随频率的变化不是线性的,说明不同频率的分量在空间/时间上传播的速度不同(延迟不同),这会导致信号在时域上发生脉冲展宽或畸变。

结语

    无论是周期信号的级数展开,还是非周期信号的连续变换,相位都是波形合成的“粘合剂”。振幅决定了素材有多丰富,而相位决定了这些素材在什么空间位置‘对齐集合’。 破坏了相位,整个信号的几何结构(如图像的边缘)将荡然无存。

Logo

openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构

更多推荐