傅里叶级数

三角函数系

正交函数系

  • 正交向量

    $$ \mathbf{a} \cdot \mathbf{b} = \sum_{i=1}^{n} a_i b_i = 0 $$

我们知道对于一个向量空间,其中任意向量可由一组 线性无关的向量线性表出(极大无关组,向量空间的一组基),对于一组基${\mathbf{v}_1, \mathbf{v}_2, \dots, \mathbf{v}_n}$ ,任一向量$\mathbf{x}$,有

$$ \mathbf{x} = \sum_{i=1}^{n} k_i \mathbf{v_i} $$

正交基就是要求一组基中的向量,两两正交,这样可以更好的求出系数,$ \mathbf{x} \cdot \mathbf{v_m} = \sum_{i=1}^{n} k_i \mathbf{v_i} \cdot \mathbf{v_m} = k_m \mathbf{v_m} \mathbf{v_m} $

  • 正交函数

深入理解正交函数 - 知乎

我们可以把函数看成一个无穷维向量

$$ \begin{equation*} f(x) = \lim_{\Delta x \to 0} \begin{pmatrix} f(a), & f(a+\Delta x), & f(a+2\Delta x), & \ldots, & f(a+(M-1)\Delta x) \end{pmatrix} \qquad x \in [a,b] \end{equation*} $$

无穷维空间(如函数空间)中,我们可以类似地把一个函数视为在某种基底上的线性组合。假设我们考虑一个函数 ( f ),例如定义在某个区间上的实值连续函数。我们可以将函数表示为某种基(例如正交基)上的线性组合:

$$ \mathbf{f(x)} = \sum_{i=1}^\infty c_i \mathbf{v(x)} $$

那么类似$\mathbf{a} \cdot \mathbf{b} = \sum_{i=1}^{n} a_i b_i = 0$两个函数正交:

$$ \langle f, g \rangle = \int_{a}^{b} f(x)g(x)dx =0 $$$$ \langle f, v_m \rangle = \int_{a}^{b} f(x)v_m(x)dx = \int_{a}^{b} \sum_{i=1}^\infty c_i v_i(x) v_m(x)dx = \int_{a}^{b} c_m v_m(x) v_m(x)dx = c_m \langle v_m, v_m \rangle $$

三角函数系

对于函数集合

$$ C = \left\{ 1, \sin x, \cos x, \sin 2x, \cos 2x, \ldots, \sin nx, \cos nx \mid x \in (-\pi, \pi) \right\} $$

称之为三角函数系

  • 正交性

    对于$\sin nx, \cos nx$,有

$$ \int_{-\pi}^{\pi} \sin nx \sin mx dx = 0, \qquad \int_{-\pi}^{\pi} \cos nx \, dx = 0 $$

由$\sin(a + b) = \sin a \cos b + \cos a \sin b$ $和\cos(a + b) = \cos a \cos b - \sin a \sin b$

$$ \int_{-\pi}^{\pi} \sin nx \sin mx dx = \int_{-\pi}^{\pi} \frac{1}{2}[\cos(m-n)x - \cos(m+n)x] dx = 0 $$$$ \int_{-\pi}^{\pi} \cos nx \cos mx dx = \int_{-\pi}^{\pi} \frac{1}{2}[\cos(m-n)x + \cos(m+n)x] dx = 0 $$$$ \int_{-\pi}^{\pi} \sin nx \cos mx dx = \int_{-\pi}^{\pi} \frac{1}{2}[\sin(n+m)x - \sin(n-m)x] dx = 0 $$
  • 特别地

    $$ \int_{-\pi}^{\pi} \sin nx \sin nx dx = \int_{-\pi}^{\pi} \frac{1}{2}[\cos(n-n)x - \cos(n+n)x] dx = \int_{-\pi}^{\pi} \frac{1}{2} (1 - \cos 2nx) dx = \pi $$$$ \int_{-\pi}^{\pi} \cos nx \cos nx dx = \int_{-\pi}^{\pi} \frac{1}{2}[\cos(m-n)x + \cos(n+n)x] dx = \int_{-\pi}^{\pi} \frac{1}{2} (1 + \cos 2nx) dx = \pi $$

傅里叶级数

周期为$2\pi$的傅里叶级数

对于一个以 $2\pi$ 为周期的函数 $f(x) = f(x + 2\pi)$,可以将其展开为傅里叶级数:

$$ f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} (a_n \cos nx + b_n \sin nx) $$

根据三角函数的正交性:

对于 $1 \cdot f(x)$,仅 $\frac{a_0}{2}$ 项的积分结果不为 0:

$$ \int_{-\pi}^{\pi} f(x) \, dx = \frac{1}{2} a_0 \int_{-\pi}^{\pi} dx = \pi a_0 $$

对于 $\cos nx \cdot f(x)$,仅 $a_n \cos nx$ 项的积分结果不为 0:

$$ \int_{-\pi}^{\pi} \cos nx \cdot f(x) \, dx = a_n \int_{-\pi}^{\pi} \cos^2 nx \, dx = \pi a_n $$

对于 $\sin nx \cdot f(x)$,仅 $b_n \sin nx$ 项的积分结果不为 0:

$$ \int_{-\pi}^{\pi} \sin nx \cdot f(x) \, dx = b_n \int_{-\pi}^{\pi} \sin^2 nx \, dx = \pi b_n $$

因此傅里叶级数中的系数满足:

$$ a_0 = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x) \, dx $$$$ a_n = \frac{1}{\pi} \int_{-\pi}^{\pi} \cos nx \cdot f(x) \, dx $$$$ b_n = \frac{1}{\pi} \int_{-\pi}^{\pi} \sin nx \cdot f(x) \, dx $$

周期为 $2L$ 的傅里叶级数

对于一个以 $2L$ 为周期的函数 $f(t) = f(t + 2L)$,可以使用变量代换 $t = \frac{L}{\pi}x$ 将其转变为一个以 $2\pi$ 为周期的函数:

$$ g(x) = f(t) = f( \frac{L}{\pi}x ) \quad \Rightarrow \quad g(x + 2\pi) = f\left(\frac{L}{\pi}(x + 2\pi)\right) $$

因此将 $x = \frac{\pi}{L}t $ 代入 $g(x)$ 的傅里叶级数中,即可得到$f(t)$ 的傅里叶级数。

定义 $\omega_0 = \frac{\pi}{L}$,三角函数换元后有:

$$ \sin nx = \sin n\omega_0 x, \qquad \cos nx = \cos n\omega_0 x $$

积分换元后有:

$$ \int_{-\pi}^{\pi} g(x) \, dx = \frac{\pi}{L} \int_{-L}^{L} f(t) \, dt $$

因此对于周期为 $2L$ 的傅里叶级数有:

$$ \begin{aligned} f(t) &= \frac{a_0}{2} + \sum_{n=1}^{\infty} \left[ a_n \cos(n\omega_0 t) + b_n \sin(n\omega_0 t) \right] \\[6pt] a_0 &= \frac{1}{L} \int_{-L}^{L} f(t) \, dt \\[6pt] a_n &= \frac{1}{L} \int_{-L}^{L} \cos(n\omega_0 t) \cdot f(t) \, dt \\[6pt] b_n &= \frac{1}{L} \int_{-L}^{L} \sin(n\omega_0 t) \cdot f(t) \, dt \end{aligned} $$

周期信号的傅里叶级数(工程常用形式)

在工程中,信号为从 0 开始、以 T 为周期的函数,因此定义变换(其中 $\omega_0$ 为基频):

$$ T = 2L, \qquad \omega_0 = \frac{\pi}{L} = \frac{2\pi}{T} $$

易得对于周期为 T 的函数:

$$ \int_{-L}^{L} f(t) \, dt = \int_{0}^{2L} f(t) \, dt = \int_{0}^{T} f(t) \, dt $$

由此可以得到工程中常用的傅里叶级数形式:

$$ \begin{aligned} f(t) &= \frac{a_0}{2} + \sum_{n=1}^{\infty} \left[ a_n \cos(n\omega_0 t) + b_n \sin(n\omega_0 t) \right] \\[6pt] a_0 &= \frac{2}{T} \int_{0}^{T} f(t) \, dt \\[6pt] a_n &= \frac{2}{T} \int_{0}^{T} \cos(n\omega_0 t) \cdot f(t) \, dt \\[6pt] b_n &= \frac{2}{T} \int_{0}^{T} \sin(n\omega_0 t) \cdot f(t) \, dt \end{aligned} $$

傅里叶级数的复数形式

傅里叶级数的复数形式

根据欧拉公式有:

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

据此以及三角函数的奇偶性,可将三角函数表示为复数形式:

$$ \sin\theta = \frac{1}{2i}(e^{i\theta} - e^{-i\theta}) = -\frac{i}{2}(e^{i\theta} - e^{-i\theta}) $$$$ \cos\theta = \frac{1}{2}(e^{i\theta} + e^{-i\theta}) $$

复数形式的展开将复数形式的三角函数代入傅里叶级数:

$$ f(t) = \frac{a_0}{2} + \sum_{n=1}^{\infty} \left[ a_n \cos(n\omega_0 t) + b_n \sin(n\omega_0 t) \right] $$$$ = \frac{a_0}{2} + \sum_{n=1}^{\infty} \left[ \frac{a_n - i b_n}{2} e^{i n \omega_0 t} + \frac{a_n + i b_n}{2} e^{-i n \omega_0 t} \right] $$$$ = \frac{a_0}{2} + \sum_{n=1}^{\infty} \frac{a_n - i b_n}{2} e^{i n \omega_0 t} + \sum_{n=-\infty}^{-1} \frac{a_n + i b_n}{2} e^{i n \omega_0 t} $$$$ = \sum_{n=-\infty}^{\infty} C_n e^{i n \omega_0 t} $$

现推导系数 $C_n$: 当 $ n = 0 $ 时:

$$ C_0 = \frac{a_0}{2} = \frac{1}{T} \int_0^T f(t) \, dt $$

当 $ n > 0 $ 时:

$$ \begin{aligned} C_n &= \frac{a_n - i b_n}{2} \\ &= \frac{1}{T} \int_0^T f(t) \cdot \left[ \cos(n\omega_0 t) - i \sin(n\omega_0 t) \right] dt \\ &= \frac{1}{T} \int_0^T f(t) \cdot e^{-i n \omega_0 t} \, dt \end{aligned} $$

当 $ n < 0 $ 时:

$$ \begin{aligned} C_n &= \frac{a_{-n} + i b_{-n}}{2} \\ &= \frac{1}{T} \int_0^T f(t) \cdot \left[ \cos(-n\omega_0 t) + i \sin(-n\omega_0 t) \right] dt \\ &= \frac{1}{T} \int_0^T f(t) \cdot e^{-i n \omega_0 t} \, dt \end{aligned} $$

故:

$$ C_n = \frac{1}{T} \int_0^T f(t) \cdot e^{-i n \omega_0 t} \, dt \qquad n \in Z $$
  • 周期函数傅里叶级数的特点

根据傅里叶级数的复数表达:

$$ f(t) = \sum_{n=-\infty}^{\infty} C_n e^{i n \omega_0 t} $$

在复平面上,傅里叶级数的每一项中:

  1. $ e^{i n \omega_0 t} $体现了一个以 $ n \omega_0 $为角速度旋转的向量
  2. $ |C_n| $ 确定了这个向量的幅值
  3. $ \angle C_n $ 确定了这个向量在 $ t = 0 $ 的相位角

对于复数 Cn=a+ib(其中 a 是实部,b 是虚部):

  • ∣Cn∣:模长(幅值),$∣C_n∣=\sqrt {a^2+b^2} $

  • ∠Cn:辐角(相位角),$∠C_n=arg(C_n)=arctan(\frac{b}{a})$ (需根据象限调整)

因此周期函数的傅里叶本质上是将周期函数分解为频率为基频的整数倍 nω0​ 的无数个周期性旋转 即周期函数的频域以 Δω=ω0 为间隔

$$ f(t+T)=f(t) \to e^{iw(t+T)}=e^{iwt} \to e^{iwT}=1 \to wT=2n\pi \to w=n \cdot \frac{2 \pi}{T} \\ w=n w_0 \to \Delta ω=ω0 $$

非周期函数的傅里叶变换

对于非周期函数 $f(t)$ 其周期 $T\to\infty$, 基频 $\omega_0 = \frac{2 \pi}{T}\to 0$ 因此函数的频率分布之间的间隔 $\omega_0=\Delta\omega\to 0$, 此时函数的频域连续分布根据以下变换

$$ n \omega_0 \to \omega \qquad T= \frac{2 \pi}{\omega _0} ;\; \frac{1}{T}= \frac{\Delta \omega}{2 \pi} \to \frac{\omega}{2 \pi} $$

带入 $f(t)$ 的傅里叶级数

$$ \begin{aligned} f(t) &= \sum_{n=-\infty}^{\infty} C_n e^{in\omega_0 t} \\ &= \frac{1}{T} \sum_{n=-\infty}^{\infty} \left( \int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) e^{-in\omega_0 t} dt \right) e^{in\omega_0 t} \\ &= \frac{d\omega}{2\pi} \int_{-\infty}^{\infty} \left( \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt \right) e^{i\omega t} \\ &= \frac{1}{2\pi} \int_{-\infty}^{\infty} \left( \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt \right) e^{i\omega t} d\omega \end{aligned} $$

此时$ C_n $ 可以被表示成一个连续的函数$ F(\omega) $,因此有

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

其中

$ F(\omega) = \mathcal{F}[f(t)] $ 为傅里叶正变换

$ f(t) = \mathcal{F}^{-1}[F(\omega)] $ 为傅里叶逆变换

根据推导可得, 傅里叶变换得到的 $F(\omega)$相当于$ C_n=F(\omega)d \omega$, 因此 $F(\omega)$ 并不代表信号中

率为 $\omega$的成分, 而是在频率范围$ \omega \sim \omega + d\omega $内的密度

因此 F(ω) 体现的是信号在 ω 的幅值密度而非幅值 (与傅里叶级数做区分)


傅里叶变换 | Tonyddg’s Noteverse

傅里叶变换和卷积

卷积

$$ f_1(x) * f_2(x) = \int_{-\infty}^{\infty} f_1(\xi) f_2(x - \xi) \, d\xi $$

卷积的本质及物理意义(全面理解卷积)

设两时域信号 $ f(t), g(t) $,对于卷积有:

$$ f(t) * g(t) = \int_{-\infty}^{\infty} f(\tau) g(t - \tau) \, d\tau $$

那么其傅氏变换为

$$ \begin{aligned} \mathcal{F}[f(t) * g(t)] &= \int_{-\infty}^{\infty} \left[ \int_{-\infty}^{\infty} f(\tau) g(t - \tau) \, d\tau \right] e^{-j\omega t} \, dt \\ &= \int_{-\infty}^{\infty} f(\tau) \left[ \int_{-\infty}^{\infty} g(t - \tau) e^{-j\omega t} \, dt \right] d\tau \\ \text{交换次序} \\ &= \int_{-\infty}^{\infty} f(\tau) e^{-j\omega \tau} \, d\tau \left[ \int_{-\infty}^{\infty} g(t - \tau) e^{-j\omega (t - \tau)} \, d(t - \tau) \right] \\ &= \mathcal{F}[f(\tau)] \, \mathcal{F}[g(t - \tau)] \\ &= F_1(\omega) F_2(\omega) \end{aligned} $$
Licensed under CC BY-NC-SA 4.0
使用 Hugo 构建
主题 StackJimmy 设计