For AI agents: the complete documentation index is available at https://shijin.space/llms.txt, the full documentation bundle is available at https://shijin.space/llms-full.txt, and this page is available as Markdown at https://shijin.space/knowledge/fourier-transform-and-poisson-summation/index.md.

10 傅里叶变换与 Poisson 求和公式

整理日期:2026-10-16。
前置笔记:02 基本区域、格行列式与正交化、07 对偶格与对偶基、09 高斯函数与离散高斯分布。

设 Λ=L(B)⊆Rn\Lambda=\mathcal L(B)\subseteq\mathbb R^n 为满秩格,Λ∗=B−TZn\Lambda^*=B^{-\mathsf T}\mathbb Z^n。

以下可取 ff 为 Schwartz 函数,即函数及其各阶导数都比任意负幂更快地衰减。高斯函数满足这一条件,后面的求和、积分和傅里叶展开均可使用。

1. 傅里叶变换

采用如下约定:

f^(y)=∫Rnf(x)e−2πi⟨x,y⟩ dx,y∈Rn.\boxed{ \widehat f(y) = \int_{\mathbb R^n} f(x)e^{-2\pi i\langle x,y\rangle}\,dx, \qquad y\in\mathbb R^n. }

其中 i2=−1i^2=-1,yy 表示频率。特别地,

f^(0)=∫Rnf(x) dx,\widehat f(0)=\int_{\mathbb R^n}f(x)\,dx,

即零频值等于函数的总积分。

后面用到两个性质:

g(x)=f(x−c)⟹g^(y)=e−2πi⟨c,y⟩f^(y);g(x)=f(x-c) \quad\Longrightarrow\quad \widehat g(y) = e^{-2\pi i\langle c,y\rangle}\widehat f(y); h(x)=f(x/s),s>0⟹h^(y)=snf^(sy).h(x)=f(x/s),\quad s>0 \quad\Longrightarrow\quad \widehat h(y)=s^n\widehat f(sy).

平移引入相位因子;空间中的伸缩对应频率尺度的反向变化。上述约定与性质见 Regev《Lecture 9: Fourier Transform》,Fall 2004,PDF 第 3 页,Definition 3 与 Theorem 3。

2. 高斯函数的傅里叶变换

沿用第九篇的记号:

ρ(x)=e−π∥x∥2,ρs(x)=e−π∥x∥2/s2.\rho(x)=e^{-\pi\|x\|^2}, \qquad \rho_s(x)=e^{-\pi\|x\|^2/s^2}.

高斯函数的一个关键性质是

ρ^=ρ.\widehat\rho=\rho.

一维情形可以用简短推导说明:对 ρ^(y)\widehat\rho(y) 求导,再分部积分,得到

ddyρ^(y)=−2πy ρ^(y),ρ^(0)=1.\frac{d}{dy}\widehat\rho(y) = -2\pi y\,\widehat\rho(y), \qquad \widehat\rho(0)=1.

因此 ρ^(y)=e−πy2\widehat\rho(y)=e^{-\pi y^2}。利用各坐标的乘积分解,可以推广到 nn 维。

再应用伸缩性质,得到

ρs^(y)=sne−πs2∥y∥2=snρ1/s(y).\boxed{ \widehat{\rho_s}(y) = s^n e^{-\pi s^2\|y\|^2} = s^n\rho_{1/s}(y). }

所以,高斯函数经过傅里叶变换后仍是高斯函数,尺度由 ss 变为 1/s1/s,并带有系数 sns^n。该结论见 Regev,PDF 第 3 页,Example 3。

对于第九篇的归一化密度 fs,0=ρs/snf_{s,0}=\rho_s/s^n,相应地有 f^s,0=ρ1/s\widehat f_{s,0}=\rho_{1/s}。

3. 为什么频率落在对偶格上

考虑一个基本振荡函数

χy(t)=e2πi⟨t,y⟩.\chi_y(t)=e^{2\pi i\langle t,y\rangle}.

要求它沿任意格向量 v∈Λv\in\Lambda 平移后保持不变,即

χy(t+v)=χy(t),\chi_y(t+v)=\chi_y(t),

等价于

e2πi⟨v,y⟩=1⟺⟨v,y⟩∈Z.e^{2\pi i\langle v,y\rangle}=1 \quad\Longleftrightarrow\quad \langle v,y\rangle\in\mathbb Z.

这个条件对所有 v∈Λv\in\Lambda 成立,恰好意味着 y∈Λ∗y\in\Lambda^*。

因此,与格周期相容的频率,正好组成对偶格。

对于一般函数 ff,其傅里叶变换仍定义在整个 Rn\mathbb R^n 上;这里的离散频率结构来自函数的格周期性。

4. Poisson 求和公式

对满秩格 Λ\Lambda,Poisson 求和公式为

∑x∈Λf(x)=1det⁡(Λ)∑y∈Λ∗f^(y).\boxed{ \sum_{x\in\Lambda}f(x) = \frac{1}{\det(\Lambda)} \sum_{y\in\Lambda^*}\widehat f(y). }

它把原格上的函数值求和,转化为对偶格上的傅里叶变换值求和。系数也可写为 det⁡(Λ∗)=1/det⁡(Λ)\det(\Lambda^*)=1/\det(\Lambda)。见 Micciancio《The Gaussians Distribution》,Fall 2019,§4,PDF 第 7–8 页,Theorem 5 与 Corollary 6。

证明的关键是将 ff 周期化:

F(t)=∑x∈Λf(t+x).F(t)=\sum_{x\in\Lambda}f(t+x).

此时 F(t+v)=F(t)F(t+v)=F(t),所以它的傅里叶级数以 Λ∗\Lambda^* 为频率集合。对 y∈Λ∗y\in\Lambda^*,记对应系数为 aya_y,则

ay=1det⁡(Λ)∫P(B)F(t)e−2πi⟨t,y⟩ dt=1det⁡(Λ)∫Rnf(u)e−2πi⟨u,y⟩ du=f^(y)det⁡(Λ).\begin{aligned} a_y &= \frac{1}{\det(\Lambda)} \int_{\mathcal P(B)} F(t)e^{-2\pi i\langle t,y\rangle}\,dt\\ &= \frac{1}{\det(\Lambda)} \int_{\mathbb R^n} f(u)e^{-2\pi i\langle u,y\rangle}\,du\\ &= \frac{\widehat f(y)}{\det(\Lambda)}. \end{aligned}

第二步使用了两个事实:基本区域的格平移覆盖整个空间;对偶条件保证平移格向量时相位不变。

于是

F(t)=1det⁡(Λ)∑y∈Λ∗f^(y)e2πi⟨t,y⟩.F(t) = \frac{1}{\det(\Lambda)} \sum_{y\in\Lambda^*} \widehat f(y)e^{2\pi i\langle t,y\rangle}.

令 t=0t=0,便得到 Poisson 求和公式。

5. 应用于格上的高斯总质量

将 f=ρsf=\rho_s 代入,得到

ρs(Λ)=sndet⁡(Λ)ρ1/s(Λ∗).\boxed{ \rho_s(\Lambda) = \frac{s^n}{\det(\Lambda)} \rho_{1/s}(\Lambda^*). }

这里 ρs(A)=∑x∈Aρs(x)\rho_s(A)=\sum_{x\in A}\rho_s(x)。将对偶格中的零向量单独取出,可写成

ρs(Λ)=sndet⁡(Λ)(1+∑y∈Λ∗y≠0e−πs2∥y∥2).\rho_s(\Lambda) = \frac{s^n}{\det(\Lambda)} \left( 1+ \sum_{\substack{y\in\Lambda^*\\y\ne0}} e^{-\pi s^2\|y\|^2} \right).

这是一个精确等式。零频项贡献 sn/det⁡(Λ)s^n/\det(\Lambda),其余贡献由非零对偶向量控制。

当非零项的总和足够小时,就有

ρs(Λ)≈sndet⁡(Λ).\rho_s(\Lambda)\approx\frac{s^n}{\det(\Lambda)}.

这个主项也符合体积直觉:高斯函数的总积分为 sns^n,格点的平均密度为 1/det⁡(Λ)1/\det(\Lambda)。

6. 平移中心与下一篇的联系

第九篇定义了 ρs,c(x)=ρs(x−c)\rho_{s,c}(x)=\rho_s(x-c)。利用平移性质,

ρs,c(Λ)=sndet⁡(Λ)∑y∈Λ∗ρ1/s(y)e−2πi⟨c,y⟩.\rho_{s,c}(\Lambda) = \frac{s^n}{\det(\Lambda)} \sum_{y\in\Lambda^*} \rho_{1/s}(y)e^{-2\pi i\langle c,y\rangle}.

去掉零频项,再利用相位因子的模为 11,得到

∣det⁡(Λ)snρs,c(Λ)−1∣≤ρ1/s ⁣(Λ∗∖{0}).\boxed{ \left| \frac{\det(\Lambda)}{s^n}\rho_{s,c}(\Lambda)-1 \right| \le \rho_{1/s}\!\left(\Lambda^*\setminus\{0\}\right). }

右侧与中心 cc 无关。因此,只要它不超过 ε\varepsilon,所有中心对应的高斯总质量相对于 sn/det⁡(Λ)s^n/\det(\Lambda) 的误差都不超过 ε\varepsilon。

这种展开与估计见 Regev 的最坏情况到平均情况归约讲义,§2,PDF 第 3–4 页,Lemma 5 的证明。

下一篇将通过平滑参数描述:为了让非零对偶向量的高斯总质量足够小,尺度 ss 需要达到多大。

引用此笔记