傅里叶变换

傅立叶变换，表示能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。在不同的研究领域，傅立叶变换具有多种不同的变体形式，如连续傅立叶变换和离散傅立叶变换。本文梳理基本知识。概述Fourier transform或Transformée de Fourier有多个中文译名，常见的有“傅里叶变换”、“付立叶变换”、“傅立叶转换”、“傅氏转换”、“傅氏变换”、等等。傅

苍蓝儿

1438人浏览 · 2021-06-05 09:33:51

苍蓝儿 · 2021-06-05 09:33:51 发布

傅立叶变换，表示能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。在不同的研究领域，傅立叶变换具有多种不同的变体形式，如连续傅立叶变换和离散傅立叶变换。本文梳理基本知识。

概述

Fourier transform或Transformée de Fourier有多个中文译名，常见的有“傅里叶变换”、“付立叶变换”、“傅立叶转换”、“傅氏转换”、“傅氏变换”、等等。

傅立叶变换是一种分析信号的方法，它可分析信号的成分，也可用这些成分合成信号。所谓信号，从狭义上说可以认为是自然界中作为信息载体的各类波，一般来说简谐震动产生的正弦波是最常见的研究对象。许多波形可作为信号的成分，比如正弦波、方波、锯齿波等，傅立叶变换用正弦波作为信号的成分。
对于自然界存在的所有波，我们可以利用所谓的傅立叶级数展开法将它们分解为有限或无限个不同频率不同振幅的正弦、余弦波的集合

定义

连续傅里叶变换

$f (t)$ 是 $t$ 的周期函数，如果 $t$ 满足狄里赫莱条件：在一个以 $2 T$ 为周期内 $f (X)$ 连续或只有有限个第一类间断点，附 $f (x)$ 单调或可划分成有限个单调区间, 则 $\mathrm{x} )$ 以 $2 T$ 为周期的傅里叶级数收玫, 和函数 $(\mathrm{x})$ 也是以 $2 T$ 为周期的周期函数，且在这些间断点上，函数是有限值；在一个周期内具有有限个极值点；绝对可积。
则有傅里叶变换成立：

$F(\omega)=\mathcal{F}[f(t)]=\int_{-\infty}^{\infty} f(t) e^{-i w t} d t$

其中 $F (ω)$ 叫做 $f (t)$ 的象函数， $f (t)$ 叫做 $F(ω)F(\omega)$ 的象原函数。 $F(ω)\mathrm{F}(\omega)$ 是 $f(t)f(\mathrm{t})$ 的象。 ${f}({t})$ 是 $F(ω)\mathrm{F}(\omega)$ 原象。

傅立叶逆变换

$f(t)=\mathcal{F}^{-1}[F(\omega)]=\frac{1}{2 \pi} \int_{-\infty}^{\infty} F(\omega) e^{i w t} d \omega$

离散傅里叶变换

设 $x (n)$ 是一个长度为 $M$ 的有限长序列，则定义 $x (n)$ 的 $N$ 点离散傅里叶变换为

$X(k)=\operatorname{DFT}[x(n)]=\sum_{n=0}^{N-1} x(n) W_{N}^{k n} \quad k=0,1, \cdots, N-1$

$X (k)$ 的离散傅里叶逆变换(Inverse Discrete Fourier Transform, IDFT)为

$x(n)=\operatorname{IDFT}[X(k)]=\frac{1}{N} \sum_{k=0}^{N-1} X(k) W_{N}^{-k n} \quad n=0,1, \cdots, N-1$

式中, $WN=e−j2πN,NW_{N}=\mathrm{e}^{-j\frac{2\pi}{N}}, N$ 称为 $DFT\mathrm{DFT}$ 变换区间长度, $\geqslant M)$ 通常称上述两式为离散傅里叶变换对。

为了叙述简洁，常常用 ${DFT[x(n)]} _ { N}$ 和 $IDFT[X(k)]_{N}$ 分别表示 $N$ 点离散傅里叶变换和 $N$ 点离散傅里叶逆变换。

实际应用

给出一幅图像，我们求出图像中圆形的周期和相位

将图像极坐标变换

叠加减去均值得到时域信号：

离散傅里叶变换，计算模长

其中能量最大的就是信号的周期 12，与实际相符
计算频率为12的相位，得到 -10.279644688897708，与实际相符

参考代码

import matplotlib.pyplot as plt
import cv2
import numpy as np
from scipy.fftpack import fft


if __name__ == '__main__':

    image_path = 'test.png'
    image = cv2.imread(image_path)

    center = (np.array(image.shape[:2]) / 2).astype('int')[::-1].tolist()
    assert center[0] == center[1]
    radius = center[0]

    size_x = 200
    size_y = 200

    output_img = cv2.warpPolar(image, (size_x, size_y), tuple(center), radius, 12)

    output_img = output_img[:, :, 0]

    signal = np.sum(output_img, axis=1)

    signal = signal - signal.mean()

    scipy_fft_y = fft(signal)

    scipy_fft_y = scipy_fft_y[:len(scipy_fft_y)//2]

    # energy_fft
    abs_y = np.abs(scipy_fft_y)

    # phase_fft
    angle_y = np.angle(scipy_fft_y)

    cycle_num = np.argmax(abs_y)
    # cycle angle
    cycle_angle = 360 / cycle_num
    # phase
    phase = cycle_angle / 2 / np.pi * angle_y[cycle_num]

    print(f"周期 : {cycle_num}")
    print(f"角度 : {cycle_angle}")
    print(f"相位 : {phase}")

    plt.plot(abs_y)
    plt.show()

    pass

参考资料

https://baike.baidu.com/item/%E5%82%85%E9%87%8C%E5%8F%B6%E5%8F%98%E6%8D%A2/7119029?fr=aladdin
https://www.cnblogs.com/amanlikethis/p/3505248.html
https://www.zhihu.com/question/21040374
https://baike.baidu.com/item/%E6%97%B6%E5%9F%9F%E9%A2%91%E5%9F%9F/9399325?fr=aladdin