如果你学过量子力学,一定记得这条公式:把波函数从位置表象换到动量表象,要做一次傅里叶变换。教科书上写得很轻松,一道积分就过去了。但很多人第一次算完会愣住:位置和动量,一个描述“在哪里”,一个描述“有多快”,凭什么它们之间的关系恰好是傅里叶变换?换一种积分变换不行吗?
这个问题我当年也纠结了很久。后来把正则对易关系、辛几何、群表示论这几块拼图凑齐,才真正看明白:位置与动量是傅里叶变换对,不是数学上的巧合,而是量子力学结构性决定的必然结果。这次我准备把这个链条从头到尾拉一遍,内容包括:
- 位置表象里的动量算符为什么是微分算子;
- 正则对易关系如何“逼”出傅里叶变换;
- 相空间辛结构下的几何图像,以及最近通信领域很热的逆辛有限傅里叶变换(ISFFT);
- 弯曲时空中这个关系会发生什么变化;
- 最后聊一个意外的主场:深度学习里的“视觉傅里叶变换”。
这篇文章适合学过一点量子力学、但一直没想透这个问题的朋友,也适合想把信号处理和量子力学两个视角打通的人。我尽量把每步推导背后的物理动机说清楚,不会只甩公式。
1. 从一道常微分方程说起
1.1 动量算符为什么是微分算子
在量子力学里,一个粒子在某个时刻的状态用波函数 ψ(x) 描述,|ψ(x)|² 代表在 x 处找到粒子的概率密度。这是位置表象,x 是自变量。
那动量怎么表示?一个朴素的想法是:粒子的动量决定了它运动的方向和快慢,而一个“匀速运动”的自由粒子,波函数必然带有相位因子 e^{i(kx - ωt)}。波长越长,动量越小;波长越短,动量越大。德布罗意关系把这一步定量化:
p = ℏk
这里 ℏ 是约化普朗克常数。问题是,如何在波函数上“提取”这个动量信息?
关键洞察来自平移对称性。假设把整个坐标轴向右平移 a,波函数变成 ψ(x - a)。如果波函数是相位因子 e^{ikx},平移后变成 e^{ik(x-a)} = e^{-ika} e^{ikx},多了一个整体相位。反过来看,这个相位因子随 x 变化的快慢,正是由 k 决定的。
如果把平移无穷小化,令 a = ε,泰勒展开:
ψ(x - ε) ≈ ψ(x) - ε ψ'(x) = (1 - ε d/dx) ψ(x)
平移的生成元就是 d/dx。既然动量是平移的“源头”,那么动量算符正比于这个生成元。加上量纲匹配后:
p̂ = -iℏ d/dx
这就是位置表象里动量算符为什么是微分算子的原因。负号不是随手加的,它保证了 e^{ikx} 是 p̂ 的本征态:
p̂ e^{ikx} = -iℏ d/dx e^{ikx} = ℏk e^{ikx} = p e^
所以,动量算符在位置表象里写成微分算子,本质上是“平移对称性”和“德布罗意关系”共同作用的结果。
1.2 平面波是位置与动量之间的“翻译官”
现在要做的是:给定一个任意波函数 ψ(x),如何读出它里面各种动量成分的分布?
因为平面波 e^{ipx/ℏ} 是 p̂ 的本征函数,一个自然策略是把 ψ(x) 按平面波展开:
ψ(x) = ∫ ψ̃(p) e^{ipx/ℏ} dp
系数 ψ̃(p) 就代表动量本征值 p 的成分有多大。这个展开式本身就是一个积分变换——把函数 ψ(x) 变成系数函数 ψ̃(p)。
求出这个系数的过程恰好的傅里叶变换。利用平面波的正交归一性:
∫ e^{ipx/ℏ} e^{-ip'x/ℏ} dx = 2πℏ δ(p - p')
两边乘上 e^{-ip'x/ℏ} 再积分,就得到:
ψ̃(p) = (1/√(2πℏ)) ∫ ψ(x) e^{-ipx/ℏ} dx
这就是教科书里的傅里叶变换公式。也就是说,从位置表象换到动量表象,做傅里叶变换这件事,本质上是“把波函数在平面波基底上展开”。
这一步让我印象很深:位置和动量之间的关系,不是人为规定的“配对”,而是由动量算符作为平移生成元这一身份自然推导出来的。动量本征态必须是指数函数 e^{ipx/ℏ},而把任意函数分解到指数函数上,就是傅里叶分析的本职工作。
需要模型API调用? 免费领10W Token,多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 正则对易关系:傅里叶变换的真正源头
2.1 对易关系包含了量子力学的全部“新意”
上面推导看起来非常顺畅,但有一个问题:为什么动量算符偏偏是 -iℏ d/dx,而不是别的算子?
更深一层的原因藏在正则对易关系里:
[x̂, p̂] = x̂p̂ - p̂x̂ = iℏ
这个关系是量子力学的核心公设之一。它到底在说什么?
如果把 x̂ 看成“测量位置”的操作,p̂ 看成“测量动量”的操作,那么对易关系告诉我们:这两个操作不能同时对角化。换句话说,不存在一个量子态,让位置和动量都有确定值。这正是不确定性原理的根源。
但它的意义远不止“测量限制”。对易关系 [x̂, p̂] = iℏ 还决定了一个群结构——海森堡群。位置平移和动量平移通过这个关系耦合在一起,形成了量子力学独特的代数结构。
可以这样理解:在经典力学里,x 和 p 是相互独立的坐标,它们可以同时取确定值。到了量子层面,x̂ 和 p̂ 不对易,意味着这两个自由度“纠缠”在一起,无法分开处理。
2.2 海森堡群告诉你为什么非得是傅里叶变换
这部分我会讲得稍微抽象一点,但这是理解“为什么必须是傅里叶变换”的关键。
考虑两个操作:
- 位置平移:e^{-ia p̂/ℏ},把波函数沿 x 方向平移 a;
- 动量平移:e^{ib x̂/ℏ},给波函数乘上一个相位 e^{ibx/ℏ},等价于在动量空间平移 b。
这两个操作不对易。根据 Baker-Campbell-Hausdorff 公式和 [x̂, p̂] = iℏ,可以得到:
e^{-ia p̂/ℏ} e^{ib x̂/ℏ} = e^{iab/ℏ} e^{ib x̂/ℏ} e^
多出来的相位因子 e^{iab/ℏ} 不是无关紧要的东西,它定义了一个群的中心扩张。这个群叫海森堡群,量子力学的整个表象理论都建立在它的表示论之上。
这里有一个深刻的定理:Stone-von Neumann 定理。它说,满足正则对易关系的不可约酉表示,在酉等价意义下是唯一的。也就是说,“位置表示”和“动量表示”看起来不同,但本质上是同一个表示的不同表象,它们之间必然存在一个酉变换把它们联系起来。
这个酉变换是什么?正是傅里叶变换。
换句话说,只要接受了正则对易关系 [x̂, p̂] = iℏ,位置表象和动量表象之间的变换形式就被唯一确定了。你不能换一个小波变换、拉普拉斯变换,因为它们改变了对易关系的表示形式。傅里叶变换不仅是选择之一,而是唯一合法选择。
这一节内容是我觉得最有说服力的论证。很多教材把傅里叶变换当成纯粹的数学工具摆在量子力学边上,没有点明它就是量子力学代数结构的化身。
2.3 不确定性原理与信号处理的时宽-带宽积是同一件事
如果你学过信号处理,肯定知道一个经典结论:信号的时域宽度和频域宽度不能同时任意小。高斯脉冲在时域越短,频谱就越宽。这个“时宽-带宽积”不等式为:
Δt · Δf ≥ 1/(4π)
这个公式和量子力学里的不确定性原理:
Δx · Δp ≥ ℏ/2
是同一个数学事实,只是换了符号。把 x → t,p → ℏk → ℏ ω/v,就完全对应上了。
我在学这两门课的时候,最震撼的时刻就是发现它们用的是同一套数学。傅里叶变换的尺度性质决定了:一个函数和它的傅里叶变换不能同时紧致。位置分布越窄,动量分布必然越宽。
高斯波包是这个不等式取等号的情形。位置波函数是高斯函数时,傅里叶变换仍然是高斯函数,而且两者的标准差乘积恰好取到最小值。这也是为什么相干态、最小不确定态在量子光学里那么重要——它们是“最接近经典”的量子态。
所以要理解不确定性原理,与其背诵海森堡显微镜的实验思想实验,不如直接看傅里叶变换的数学性质:不确定性是波的本质属性,不是你测量手法不够好。
3. 相空间的视角:辛结构与逆辛有限傅里叶变换
3.1 相空间旋转:把傅里叶变换看成九十度转动
前面讲的是代数层面的必然性,这一节换个几何视角。
在经典力学里,系统的状态可以用相空间 (x, p) 中的一个点来描述。相空间不是普通空间,它上面有一个自然的几何结构叫辛结构,用微分形式表示出来就是 dp ∧ dx。辛结构的特殊之处在于:所有保持辛结构的坐标变换,都对应着一个合法的经典正则变换。
傅里叶变换在这个框架里恰好扮演一个特殊角色:它是相空间里的一个 π/2 旋转。经过傅里叶变换,坐标 (x, p) 变成 (-p, x),位置和动量互换位置(带一个负号)。这个变换保持 dp ∧ dx 不变,所以它是一个辛变换。
把这个观点量子化,就得到相空间量子力学的语言。一个量子态不再只由 ψ(x) 或 ψ̃(p) 单独刻画,而是由整个相空间上的 Wigner 函数 W(x, p) 来描述。Wigner 函数不是普通的概率分布,因为它在非经典区域可以取负值,但它的边缘分布是有意义的:
∫ W(x, p) dp = |ψ(x)|²
∫ W(x, p) dx = |ψ̃(p)|²
注意到这两条边缘分布分别由波函数和它的傅里叶变换给出,所以 Wigner 函数本身也携带了位置和动量之间的“傅里叶对偶”信息。
用这个视角看,位置和动量互相对偶,不是一种奇怪的特殊关系,而是辛几何的自然结论:它们本来就是相空间里一对共轭的规范化坐标。
3.2 有限辛空间里的新热点:ISFFT
前面对论的都是连续相空间。但在实际工程和有限系统中,位置和动量常常只取有限个离散值。比如,通信系统里的时延(τ)和多普勒频移(ν)就是一组共轭变量,它们定义在一个 N × M 的离散网格上。
这时候,连续傅里叶变换要换成有限离散版本。但“有限”带来一个新问题:普通离散傅里叶变换(DFT)是一位时间序列的变换,而时延-多普勒域是二维的,且 τ 和 ν 之间存在类似 x 和 p 的辛耦合。
处理这个问题的工具叫辛有限傅里叶变换(Symplectic Finite Fourier Transform,SFFT),它的逆变换叫逆辛有限傅里叶变换(ISFFT)。公式长这样:
X_tf[n,m] = (1/√(NM)) Σ_{n'=0}^{N-1} Σ_{m'=0}^{M-1} X_dd[n',m'] e^
对应的逆变换:
X_dd[n',m'] = (1/√(NM)) Σ_{n=0}^{N-1} Σ_{m=0}^{M-1} X_tf[n,m] e^
注意看指数里的符号:一个维度是正号,另一个维度是负号。这个符号差异正是辛结构的体现——它保证了时延和多普勒这一对变量在变换时保持共轭关系,而不是像普通二维傅里叶那样两个维度完全对称。
这个 ISFFT 是 OTFS(正交时频空间)调制技术的基础。OTFS 是 6G 通信里的一个热门候选方案,它的核心思想就是:先把数据符号放到时延-多普勒域,用 ISFFT 映射到时频域,再通过传统 OFDM 发送。多普勒频移在时延-多普勒域里是静态的,信道估计和均衡都会好做很多。
我第一次看到 OTFS 的调制框图时,瞬间就明白为什么它要叫“逆辛”傅里叶变换了——这本质上就是在离散的相空间里做坐标变换,和量子力学里位置-动量变换是同一个数学结构的工程应用。
4. 到了弯曲时空,傅里叶变换还成立吗
4.1 弯曲时空没有全局平面波
前面讨论的都是平直时空里的量子力学。但标题里提到了时空弯曲,这部分来回答一个更有挑战的问题:在弯曲时空中,位置和动量还是傅里叶变换对吗?
先说结论:不再是,至少不是全局意义上的傅里叶变换对。
原因在于,平面波 e^{ipx/ℏ} 的定义高度依赖平移对称性。e^{ipx/ℏ} 被定义为动量算符的本征函数,而动量算符是空间平移的生成元。在平直时空中,空间各处平移对称性完全一致,所以平面波可以覆盖整个空间。
但弯曲时空没有这个性质。引力会带来潮汐效应,不同位置的“方向”和“距离”定义会变,空间不再有全局的平移对称性。你无法在整个弯曲时空上用单一的一簇平面波展开波函数,因为“动量守恒”不再成立。
更技术化的说法是:弯曲时空里没有全局的 Killing 矢量场,所以不存在全局意义下的能量-动量守恒,也就不存在全局的傅里叶共轭变量对。
4.2 局域惯性系里的急救方案
没有全局傅里叶变换,物理学家怎么办?
标准做法是“化整为零”。在弯曲时空的每一个点附近,可以建立一个局域惯性系(黎曼法坐标),在这个小范围内时空近似平直,傅里叶变换照常使用。然后在不同点之间,通过联络(平行移动)把局域动量连接起来。
这就像在地球表面上处理问题:整个球面不是平直的,但在你脚下的一小块区域,可以当成平面来处理。你可以在这一小块上做傅里叶分析,但想用同一个平面坐标系描述整个球面,就会出错。
量子场论在弯曲时空中的很多计算,比如霍金辐射的推导,就是用这种局域模式展开的办法。在黑洞视界附近建立一组局域平面波,经过引力传播到无穷远后,这组模式会被拉伸和变形,最后在无穷远处看起来像是带有热谱的粒子分布。
这个过程中,局域的傅里叶变换依然成立,但“全局的”位置-动量对偶被时空曲率破坏了。
4.3 广义不确定原理:最小长度修正
再往前走一步,量子引力研究提出了更激进的修改。如果时空在普朗克尺度下有颗粒性或最小长度,那么位置的本征态不可能无限精确地定域,位置和动量的对易关系也会被修正。
这种修正常常写成广义不确定原理(GUP)的形式:
Δx Δp ≥ (ℏ/2) [1 + β (Δp)²]
等价地,有效对易关系变成:
[x̂, p̂] = iℏ (1 + β p̂²)
注意,这里的 β 通常正比于 (ℓ_P/ℏ)²,ℓ_P 是普朗克长度。当 Δp 远小于普朗克尺度关联的动量时,修正项可以忽略,回到标准傅里叶变换;但当动量接近普朗克尺度时,修正变得重要。
这个修正带来的后果是:动量本征态对应的波函数不再是指数平面波,而是一个被修正了的“变形平面波”。位置和动量之间的变换不再是我前面写的标准傅里叶积分,而是一种变形傅里叶变换。
这是不是意味着“位置与动量是傅里叶变换对”这句话在极端条件下会被修改?是的。这正是很多量子引力现象学工作里讨论的内容:给出一个 GUP 模型,算出对应的变形傅里叶变换核,然后预言高能实验里可能观测到的偏差。
不过要泼一盆冷水:到目前还没有任何实验证据支持具体的 GUP 模型。但这确实是理论推演的合理方向,也说明“位置-动量互为傅里叶变换对”这个命题,在更深层的物理里是有继续深挖空间的。
5. 傅里叶变换的现代应用视角:视觉傅里叶变换
5.1 卷积神经网络的频域本质
聊完物理,再说一个最近几年很热的词:视觉傅里叶变换。
你可能在 AI 相关文章里见过这个词,它指的不是某种新算法,而是对卷积神经网络(CNN)的另一层理解:卷积操作本质上就是傅里叶域里的滤波。
卷积定理告诉我们:两个函数卷积的傅里叶变换,等于各自傅里叶变换的逐点乘积。换言之:
卷积在空间域 = 相乘在频率域
CNN 的卷积层在图像上滑动一个小卷积核,本质上就是在对图像做局部的频率滤波。卷积核的权重不同,滤波的效果就不同——有些提取高频边缘,有些提取低频平滑信息。
更直观的证据是,把训练好的 CNN 第一层卷积核可视化,你会看到很多类似 Gabor 滤波器的图案:正弦条纹乘以高斯窗。Gabor 滤波器在信号处理里就是“带定位功能的傅里叶基”,同时提取空间位置和频率信息,和量子力学里的相干态有异曲同工之妙。
所以现在的视觉研究里,经常把 CNN 前几层理解为“局部傅里叶分析器”,把后面几层理解为“从局部频率信息中组合出语义结构”的过程。
5.2 频率偏置与注意力机制
除了 CNN,傅里叶视角还能解释深度学习中一些有意思的现象。
一个是“谱偏置”(spectral bias)。很多研究观察到,神经网络训练时先学到的往往是函数的低频分量,高频细节后学。这种偏置在傅里叶域里看非常清晰:网络的输出误差在训练初期主要集中在高频段。低频先收敛,高频后收敛,这和傅里叶基函数的“光滑度排序”完全一致。
另一个是注意力机制(Attention)的频域解释。Transformer 里的自注意力层会对特征图做全局加权,这个加权在频域里等于对频谱做一次非均匀调制。有人分析过,多头注意力的每个头可以看成在不同频率波段上做筛选,一个头关注全局低频结构,另一个头关注局部高频细节。从“视觉傅里叶变换”的角度看,Transformer 就是在用一种更灵活的方式做频率分解。
有意思的是,近年来也有模型直接“把参数放在频率域”。比如 Fourier Neural Operator(FNO),它直接在频域里做神经网络层的变换,用来求解偏微分方程。这说明傅里叶变换作为底层工具,正在从信号处理、量子力学一路渗透到最前沿的机器学习架构里。
6. 几个容易被绕进去的坑
6.1 误区速查表
| 误区 | 正解 |
|---|---|
| 位置和动量互为傅里叶对是量子力学特有的 | 信号处理里的时频对偶早就存在,量子力学只是给它加上了 ℏ 量纲 |
| 不确定性原理是因为测量扰动了粒子 | 不定性根源是波函数本身的傅里叶宽度,与测量手段无关 |
| 傅里叶变换是可以随意选择的其他积分变换之一 | 正则对易关系固定后,Stone-von Neumann 定理决定了它在表示论上的唯一性 |
| 弯曲时空里没有傅里叶变换 | 严格说是没有全局傅里叶变换,局域惯性系内仍然成立 |
| “视觉傅里叶变换”是一种新算法 | 它是对 CNN/Transformer 频域行为的一种解释性视角 |
这个表格里的误区,我自己几乎全部踩过一遍。特别是第一条:我最初在信号处理课上学过傅里叶变换,后来学量子力学时完全没有意识到这就是同一个数学对象。直到把两个学科串联起来,才知道知识是可以这样迁移的。
6.2 给学习者的实操建议
如果你想亲手验证“位置和动量是傅里叶变换对”,我强烈建议做一个简单的数值实验:造一个高斯波包,把它傅里叶变换到动量空间,看看两边宽度乘积。
用 Python 的 numpy 就能搞定,核心思路就这几行:
python复制import numpy as np
sigma = 1.0
x = np.linspace(-30, 30, 60000)
dx = x[1] - x[0]
# 高斯波包:位置空间的概率幅
psi = np.exp(-x**2 / (2 * sigma**2))
# 傅里叶变换到动量空间
psi_k = np.fft.fftshift(np.fft.fft(psi))
# 动量坐标(忽略常数因子)
k = np.fft.fftshift(np.fft.fftfreq(len(x), d=dx))
跑完之后,分别算位置分布的宽度和动量分布的宽度,你会发现两者的乘积约等于 1/2(在合适的量纲约定下)。这个实验比我当年背十遍不确定性公式都管用。
如果你想深入理论层面,建议按这个顺序读:
- 先看 Griffiths《量子力学概论》里“动量空间波函数”那一节,理解基本积分;
- 再看 Sakurai《现代量子力学》关于海森堡群和 Stone-von Neumann 定理的讨论,理解唯一性;
- 有兴趣的话,翻一翻 OTFS 相关的综述文章,看 ISFFT 在工程里怎么用;
- 最后读一些弯曲时空量子场论的入门讲义,理解为什么全局傅里叶变换会失效。
这个路径不会太陡,每一级需要的数学工具都是上一级的自然延伸。
最后再分享一个个人体会
我真正把位置和动量理解为傅里叶变换对,不是在看量子力学教材时,而是在同时学过信号处理和量子场论之后。两边用的是完全一样的数学,讲的是不同物理,但一旦打通,很多公式就不再需要死记:你只需要知道对易关系、知道平面波是平移本征态、知道傅里叶变换是把函数拆成平面波叠加,剩下的自然就会推导出来。
如果你看完这篇文章,下次再看到傅里叶变换,能下意识想到它背后的对易关系和辛结构,那这篇东西就没白写。位置与动量之间这道“看似普通的积分”,其实藏着量子力学最深的秘密:一个体系里所有关于不确定性的故事,都写在那个小小的 [x̂, p̂] = iℏ 里。
