Skip to content

定理Theorem

Whittaker–Shannon采样定理

Whittaker-Shannon sampling theorem · Nyquist-Shannon sampling theorem · 带限采样定理 · Shannon sampling theorem

从有限频带的正交展开恢复整实线函数,核清采样步长、连续代表、能量等式和一致误差,并区分平方可和噪声与逐样本有界噪声。

一串格点值通常不能确定格点之间的曲线。采样定理增加了一份实质输入:函数没有超出指定区间的频率。这个输入限制整条函数,而不是说手头的离散频谱图看起来没有高频。

形式陈述 ​

先为点采样选定代表 ​

沿Fourier变换的规范,

f^(ξ)=∫Rf(t)e−2πitξdt.

固定B>0,记

PWB={f∈L2(R):f^(ξ)=0 对几乎处处 |ξ|>B}.

这是本页的带限函数类。L2变换最初是等价类,但频带有限,Cauchy–Schwarz给 ‖f^‖1≤2B‖f^‖2。因此选定连续代表

(1)f(t)=∫−BBf^(ξ)e2πitξdξ.

本页所有样本都取这个代表的值。不能另改一个格点的值,再声称L2频谱没变所以采样定理仍适用。

定义 sincx=sin⁡(πx)/(πx),在x=0时取1。若步长满足

(2)0<h≤12B,

则每个 f∈PWB 都有

(3)f(t)=∑k∈Zf(kh)sinc(th−k).

对称部分和在L2及整个实线的一致范数中收敛。每个固定t的级数还绝对收敛。能量与尾误差满足

(4)‖f‖22=h∑k∈Z|f(kh)|2,(5)‖f−SNf‖22=h∑|k|>N|f(kh)|2,‖f−SNf‖∞≤(∑|k|>N|f(kh)|2)1/2.

这里 SNf 为式(3)中|k|≤N的有限和。h的因子只出现在L2等式,不能把两个误差对象混写。

唯一性依赖全部样本和带宽 ​

在式(2)的带限类内,两份函数若所有f(kh)相同,就处处相同。这个结论使用双边无限样本;有限样本需要式(5)所需的尾信息或其他先验约束。定理没有从有限输入免费恢复无限精度。

直觉

图中黑点是右侧两函数共有的样本。这里只画有限区间,右图函数在全部格点不可区分的证明见本页式(9)及终点;不能从有限图形单独推断带宽。

在频域把样本认成系数 ​

令 Ih=[−1/(2h),1/(2h)],将 F=f^ 在本来没有频率的部分补零。I_h的长度为1/h,函数

ek(ξ)=he−2πikhξ,k∈Z,

由Fourier级数的完备指数基组成L2(I_h)的正交规范基。本站内积对第一变量共轭线性,因此

⟨ek,F⟩=h∫IhF(ξ)e2πikhξdξ=hf(kh).

频域展开是

(6)F(ξ)=h∑k∈Zf(kh)e−2πikhξ于 L2(Ih).

这里并没有预先假定时间域的采样级数成立;是把已存在的频谱在有限区间展开。

先反演有限和,再取极限 ​

每个有限项的逆变换为

(7)h∫Ihe2πi(t−kh)ξdξ=sinc(th−k).

所以频域正交投影的逆变换正是S_Nf。Plancherel酉性把L2收敛和正交尾能量原样搬回实线,得到式(4)和式(5)第一式。这里始终先积分有限和,不把任意L2级数当作绝对可积二重和。

对频域尾项 RN=F−h∑|k|≤Nf(kh)e−2πikhξ,其支撑仍在I_h。于是

|∫IhRN(ξ)e2πitξdξ|≤h−1/2‖RN‖2=(∑|k|>N|f(kh)|2)1/2.

右侧不含t且趋零,给全实线一致收敛。逐点采样之所以在这里能由L2信息控制,关键正是尾项始终留在一个共同的有限频带。

绝对收敛不靠样本绝对可和 ​

对任意实x,在L2([−1/2,1/2])中将单位范数函数 e−2πixξ 展开到整数指数基。Parseval给

(8)∑k∈Z|sinc(x−k)|2=1.

与式(4)合用Cauchy–Schwarz,

∑k|f(kh)sinc(t/h−k)|≤(∑k|f(kh)|2)1/2.

即使样本不属于ℓ1,仍能得到这里的绝对收敛。不能反过来据此声称单独的 ∑k|sinc(x−k)| 也收敛。

例子与边界

一份有无限尾的简单信号 ​

取f(t)=sinc(t)。它是矩形频谱 1[−1/2,1/2] 的逆变换,属于PW_{1/2}且范数平方为1;它不属于L1。h=1时样本仅f(0)=1,其他整数点为零,式(3)精确退化成一个sinc项。

这个例子解释为何本页不应强加f∈L1,也不应把有限个非零样本理解成函数具有有限时间支撑。sinc在样本间仍有无限远的振荡尾。

低于采样速率要求时,真的会丢掉函数 ​

若h>1/(2B),可取 0<ε≤2B−1/h,定义

(9)g(t)=sin⁡(πt/h)sinc(εt).

sinc(εt)属于L2,乘有界正弦后仍属于L2。它的矩形频谱被正弦拆成中心在 ±1/(2h) 的两份平移,每份半宽ε/2,故仍在[−B,B]中。g并非零函数,却有g(kh)=0对全部整数k成立。因此零信号与g完全不可区分,任何只读这些样本的恢复算法都不可能对二者同时正确。

在临界h=1/(2B)处,裸正弦sin(πt/h)也有零样本,但它不平方可积,不能用来反驳本页定理。频带端点上的单个点在L2频谱中为零测集;点质量频谱属于另一种函数或分布合同。

只给有限样本会遗漏什么 ​

在临界带宽1/(2h)下,函数sinc(t/h−(N+1))在全部|k|≤N的采样点取零,却在t=(N+1)h处为1。因此这些有限零样本不能认证全函数为零。

如果已知未读样本的平方和上界,式(5)才给有限证书;否则单纯增加一个“已采N点”的标签并未补足缺失信息。

推论与应用

平方可和噪声具有精确尺度 ​

若样本扰动序列 (ak)∈ℓ2,令

e(t)=∑kaksinc(t/h−k).

同一正交展开给

(10)‖e‖2=h‖(ak)‖ℓ2,‖e‖∞≤‖(ak)‖ℓ2.

这里重建误差属于PW_{1/(2h)};若原信号严格过采样,任意扰动未必还具有原来更窄的带宽。能量稳定性没有悄悄替噪声加上那个约束。

每点误差小,不等于整个噪声序列能量小 ​

固定η>0和观察位置t/h=1/2。对|k|≤N令

ak=ηsgnsinc(1/2−k),

其余项为零。每个样本误差至多η,但这组噪声在观察点的输出为

(11)e(h/2)=ηπ∑k=−NN1|k−1/2|.

调和和随N无界增长,所以不存在对所有N都有效的常数C,使未加窗sinc重建总满足 ‖e‖∞≤Cη。这与式(10)不矛盾:这些噪声序列的ℓ2范数也随样本数增长。

严格余带允许改用平方衰减的重建核,在保持准确重建的同时给出有界样本噪声及有限观察窗的显式预算。更密采样本身不改变式(11)所用核的调和尾,必须把新的核和误差合同一起写出。

与有限周期FFT的边界 ​

谱混叠在有限周期网格上把整数模态按网格数折回;本页是在非周期整实线上,以连续频带和双边无限格点建立对应。FFT可以参与实际近似,但有限窗口、周期延拓及核截断都还需要各自的误差分析。

参考资料
关系图谱11 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系