找回密码
立即注册
搜索
热搜: Java Python Linux Go
发回帖 发新帖

4738

积分

0

好友

616

主题
发表于 2 小时前 | 查看: 6| 回复: 0

L²空间中傅里叶级数的正交投影与阶梯函数示例幻灯片

这张图的核心思想是:

把一个复杂的周期函数,当成“无限维空间里的向量”,再用有限个正弦波和余弦波去做最接近它的近似。

这个“最接近的近似”,就是函数在三角函数空间上的正交投影;投影所得的系数,就是傅里叶系数

可以把它想成音乐分析:

  • 原函数 $f(x)$:一段复杂声音;
  • 三角函数基:不同频率的纯音;
  • 傅里叶系数:每种纯音需要放多大音量;
  • 傅里叶级数:把这些纯音重新叠加,恢复原来的声音。

一、图片最上方:三个函数集合的关系

图中研究的是以 $2\pi$ 为周期的函数。最上面给出的包含关系是:

$$L^2(\mathbb{R},\mathbb{C}) \subseteq L^1(\mathbb{R},\mathbb{C})$$

1. $L^1$:绝对值的面积有限

通俗地说:

在一个周期内,函数曲线与横轴之间的总面积不能无限大。

2. $L^2$:平方后的面积有限

物理上,$|f(x)|^2$ 常常代表能量密度,所以 $L^2$ 函数也常被理解为:

一个周期内总能量有限的函数。

在有限区间上,平方可积通常比绝对可积要求更强,因此有上面的包含关系。

3. 三角多项式:有限项的正弦、余弦之和

它指有限项的正弦、余弦之和,例如:

$$a_0 + \sum_{n=1}^{N} \big(a_n \cos(nx) + b_n \sin(nx)\big)$$

这种函数显然也属于 $L^2$,因此它是 $L^2$ 周期函数的一个子集。

这条包含关系表达的是:

三角多项式只是 $L^2$ 周期函数中的一小部分,但我们希望用这一小部分去逼近更复杂的函数。


二、函数也可以有“内积”

图中定义函数之间的内积为:

$$\langle f, g \rangle = \int_{-\infty}^{+\infty} \overline{f(x)}\, g(x)\, dx$$

如果 $f$$g$ 都是实函数,复共轭符号可以忽略,于是简化为:

$$\langle f, g \rangle = \int_{-\pi}^{\pi} f(x) g(x)\, dx$$

它与普通向量点积非常相似:

  • 普通向量的点积:衡量两个箭头方向有多相似;
  • 函数内积:衡量两条函数曲线整体上有多相似。

当内积为零时

就说 $f$$g$ 正交。这不是说两条曲线在图上垂直,而是说:

它们在整个区间上的正负贡献恰好相互抵消。

例如 $\sin x$$\cos x$ 在一个周期上的内积为零,所以它们正交。


三、傅里叶分析使用的正交归一系统

图中定义:

$$B_k = \left\{ x \mapsto \frac{1}{\sqrt{\pi}},\ x \mapsto \cos(kx),\ x \mapsto \sin(kx) \right\}$$

严格地说,每一项都是一个函数。这个系统一共有 $N = 2k + 1$ 个函数。

为什么常数函数写成 $\frac{1}{\sqrt{\pi}}$

因为这里希望每个基函数的“长度”都是 1。函数长度由范数

$$\|f\| = \sqrt{\langle f, f \rangle}$$

定义。对于常数函数 $c$,若要求 $\|c\| = 1$,则:

$$\int_{-\pi}^{\pi} c^2\, dx = 2\pi c^2 = 1 \Rightarrow c = \frac{1}{\sqrt{\pi}}$$

同样可以验证,$\cos(kx)$$\sin(kx)$ 的范数也都是 1,而不同频率之间的内积都是零。所以这些函数彼此正交,而且长度都是 1。这就是图中的 ONS,即标准正交系统


四、为什么“正交”以后特别容易求傅里叶系数?

把这组基函数写成:

$$B_k = \{ h_1, h_2, \dots, h_N \},\quad N = 2k+1$$

它们张成一个有限维空间 $U_N$$U_N$ 中的任意函数都可以写成:

$$\sum_{i=1}^{N} a_i h_i$$

例如 $N = 5$ 时,$U_N$ 中只允许使用:

  • 常数项;
  • 频率为 1、2 的余弦波;
  • 频率为 1、2 的正弦波。

五、图中的几何图:傅里叶部分和就是正交投影

图中绿色平面代表 $U_N$,蓝色箭头代表原函数 $f$。虽然真实的函数空间是无限维的,但我们仍然可以用普通三维几何来想象:

  • $f$ 是空间中的一个向量;
  • $U_N$ 是一个“平面”;
  • $f$ 垂直投影到 $U_N$
  • 投影结果记为 $\hat{x}_N(f)$

公式是:

$$\hat{x}_N(f) = \sum_{i=1}^{N} a_i h_i$$

其中 $a_i = \langle h_i, f \rangle$ 就是 $f$ 在第 $i$ 个基函数方向上的“投影长度”。它们正是傅里叶系数。

为什么说它是“最好的近似”?

在所有 $U_N$ 中的函数 $g$ 里,$\hat{x}_N(f)$$f$$L^2$ 距离最小:

$$\| f - \hat{x}_N(f) \| \le \| f - g \|, \quad \forall g \in U_N$$

这里的距离是:

$$\| f - g \| = \left( \int_{-\pi}^{\pi} |f(x) - g(x)|^2 dx \right)^{1/2}$$

所以“最好”指的不是每一个点都误差最小,而是:

在整个区间上,误差平方的平均值最小。

就像照片压缩时,不可能保证每个像素都完全相同,但可以让总体误差尽可能小。


六、傅里叶部分和的具体公式

图中定义:

$$F_k(f)(x) = \frac{\bar{a}_0}{2} + \sum_{n=1}^{\infty} \left[ a_n \cos(nx) + b_n \sin(nx) \right]$$

这叫做第 $k$ 阶傅里叶部分和。它只使用频率 $1, 2, \dots, k$ 的正弦和余弦,频率更高的波暂时不使用。

1. 常数系数

展开式中的常数项是:

$$\bar{a}_0 = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x)\, dx$$

所以它与通常写法中的 $\frac{a_0}{2}$ 对应。

2. 余弦系数

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

它表示:原函数中含有多少频率为 $n$ 的余弦波。

3. 正弦系数

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

它表示:原函数中含有多少频率为 $n$ 的正弦波。


七、图中的例子:半边为 1,半边为 0 的阶梯函数

定义:

$$f(x) = \begin{cases} 1, & x \in (0, \pi) \\ 0, & x \in [\pi, 2\pi] \end{cases}$$

然后把它按照 $2\pi$ 周期不断重复。它看起来像一个周期性的开关:

  • 左半个周期:开,数值为 1;
  • 右半个周期:关,数值为 0。

八、计算常数项

因为 $f(x)$ 只在 $(0, \pi)$ 上等于 1,所以:

$$\bar{a}_0 = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x)\, dx = \frac{1}{\pi} \int_{0}^{\pi} 1\, dx = \frac{1}{\pi} \cdot \pi = 1$$

因此展开式中的常数项是:

$$\frac{\bar{a}_0}{2} = \frac{1}{2}$$

这也很好理解:函数一半时间等于 1,一半时间等于 0,所以平均值是 $\frac{1}{2}$


九、为什么所有余弦系数都是零?

计算余弦系数:

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

积分得:

$$a_n = \frac{1}{\pi} \left[ \frac{\sin(nx)}{n} \right]_{0}^{\pi} = 0$$

因为 $\sin(n\pi) = 0$ 对所有整数 $n$ 都成立。因此这个函数的傅里叶级数里没有任何非恒定余弦项。

更深一层的直观解释

把平均值 $\frac{1}{2}$ 减掉,得到:

$$g(x) = f(x) - \frac{1}{2}$$

这个 $g(x)$ 是奇函数,而余弦是偶函数。奇函数乘偶函数仍然是奇函数,在对称区间上积分等于零,所以所有余弦系数消失。


十、计算正弦系数

因为 $f(x)$ 只在 $(0, \pi)$ 上等于 1,所以:

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

计算积分:

$$b_n = \frac{1}{\pi} \left[ -\frac{\cos(nx)}{n} \right]_{0}^{\pi} = \frac{1}{\pi} \left( -\frac{\cos(n\pi)}{n} + \frac{\cos(0)}{n} \right) = \frac{1}{\pi} \left( -\frac{(-1)^n}{n} + \frac{1}{n} \right)$$

现在分两种情况。

$n$ 为偶数

此时 $(-1)^n = 1$,所以:

$$b_n = \frac{1}{\pi} \left( -\frac{1}{n} + \frac{1}{n} \right) = 0$$

偶数频率全部消失。

$n$ 为奇数

此时 $(-1)^n = -1$,所以:

$$b_n = \frac{1}{\pi} \left( -\frac{-1}{n} + \frac{1}{n} \right) = \frac{1}{\pi} \cdot \frac{2}{n} = \frac{2}{\pi n}$$

因此:

$$b_n = \begin{cases} 0, & n \text{ 为偶数} \\ \dfrac{2}{\pi n}, & n \text{ 为奇数} \end{cases}$$

这说明:

偶数频率全部消失,只剩下第 1、3、5、7……个奇数频率。


十一、最终得到的傅里叶级数

把各项系数代回去,得到:

$$f(x) \sim \frac{1}{2} + \frac{2}{\pi} \sin x + \frac{2}{3\pi} \sin(3x) + \frac{2}{5\pi} \sin(5x) + \cdots$$

更紧凑地写成:

$$f(x) \sim \frac{1}{2} + \frac{2}{\pi} \sum_{m=0}^{\infty} \frac{\sin\big((2m+1)x\big)}{2m+1}$$

这里的“$\sim$”表示这是 $f$ 的傅里叶级数,而不是说有限项时已经完全相等。


十二、它是怎样逐步逼近阶梯函数的?

只取第一项:

$$f(x) \approx \frac{1}{2} + \frac{2}{\pi} \sin x$$

这只是一条比较圆滑的波,离方形阶梯还较远。

加入第 3 次谐波:

$$f(x) \approx \frac{1}{2} + \frac{2}{\pi} \sin x + \frac{2}{3\pi} \sin(3x)$$

顶部和底部开始变平,跳跃位置开始变陡。

再加入第 5 次谐波:

$$f(x) \approx \frac{1}{2} + \frac{2}{\pi} \sin x + \frac{2}{3\pi} \sin(3x) + \frac{2}{5\pi} \sin(5x)$$

继续加入第 7、9、11 次谐波,曲线会越来越像阶梯函数。

可以把它理解为:

低频波负责搭出大轮廓,高频波负责补充尖锐的边缘。


十三、最下面的图:为什么跳跃处会有波纹?

图中蓝色曲线在平坦区域已经接近 1 和 0,但在跳跃附近出现了明显的振荡和超调。这就是吉布斯现象

在跳跃点附近:

  • 曲线会超过上方的 1;
  • 也可能低于下方的 0;
  • 加入更多项后,振荡区域会越来越窄;
  • 但最大超调高度不会完全消失。

超调量大约是跳跃高度的 $9\%$。这里跳跃高度为 1,所以曲线可能达到约 $1.09$,或下降到约 $-0.09$


十四、在跳跃点,级数到底收敛到哪里?

$x = \pi$ 处,函数左极限为 $1$,右极限为 $0$。傅里叶级数在跳跃点收敛到左右极限的平均值:

$$\frac{1 + 0}{2} = \frac{1}{2}$$

所以即使原定义规定 $f(\pi) = 0$,傅里叶级数在 $x = \pi$ 仍然趋向 $\frac{1}{2}$

同样,由于函数是周期延拓的,在 $x = 0$ 也存在跳跃,级数同样收敛到 $\frac{1}{2}$


十五、图片中一个需要准确理解的地方

图片标题写的是 “Fourier Transform”,但这一页真正讲的是:

傅里叶级数,而不是通常意义上定义在整个实轴上的傅里叶变换。

两者的区别是:

  • 傅里叶级数:处理周期函数,频率是离散的 $n$
  • 傅里叶变换:通常处理非周期函数,频率是连续的。

此外,图中 $F_k(f)(x)$ 是第 $k$ 个傅里叶部分和。完整的傅里叶级数应理解为部分和序列 $F_k(f)(x)$ 或者对应的无限三角级数。


十六、整张图的逻辑主线

整张图可以浓缩成下面这条路线:

先选择正交归一的三角函数系统,再通过内积计算每个方向上的投影长度,然后把这些投影重新组合,得到由有限个正弦波和余弦波组成的最佳近似。

随着 $k$ 增大,部分和会逐渐逼近原函数。至少在这个分段常数例子中:

  • 连续点处趋向原函数;
  • 跳跃点处趋向左右极限的平均值;
  • 跳跃附近出现吉布斯振荡。

一句话总结:

傅里叶级数就是把函数投影到不同频率的正弦、余弦方向上,再把这些频率分量重新叠加起来。




上一篇:AI Native SDLC 实战指南:Anthropic 六阶段从 intent.md 到闭环
下一篇:GPUIX 发布:React 与 Vue 直连 GPUI,前端桌面新路线?
您需要登录后才可以回帖 登录 | 立即注册

手机版|小黑屋|网站地图|云栈社区 ( 苏ICP备2022046150号-2 )

GMT+8, 2026-8-28 03:51 , Processed in 0.797957 second(s), 40 queries , Gzip On.

Powered by Discuz! X3.5

© 2025-2026 云栈社区.

快速回复 返回顶部 返回列表