1.3 为什么需要噪声:逃离伪极小值#

玻尔兹曼分布提供了对平衡态的概率描述,但它本身并没有解释为什么必须在神经动力学中引入随机性。要理解噪声的关键作用,我们需要先考察没有噪声时的情形:确定性的 Hopfield 网络,以及它众所周知的陷入 伪极小值 的倾向。

Hopfield 网络:确定性下降#

Hopfield 网络是一类具有对称权重 wij=wjiw_{ij} = w_{ji} 且无自连接 wii=0w_{ii} = 0 的循环神经网络。其状态通过异步更新演化:每个时间步选择一个单元 ii,并由其总输入的符号决定新状态:

xisign(jiwijxj+bi)x_i \leftarrow \text{sign}\left( \sum_{j \neq i} w_{ij} x_j + b_i \right)

该更新规则具有一个关键性质:它永远不会增加网络能量。为了看到这一点,考虑单元 iixix_i 翻转到 xix_i' 时的能量变化:

ΔE=E(x)E(x)=(xixi)(jiwijxj+bi)\Delta E = E(\mathbf{x}') - E(\mathbf{x}) = -(x_i' - x_i)\left( \sum_{j \neq i} w_{ij} x_j + b_i \right)

如果输入符号与新状态一致,能量会下降;如果与旧状态一致,能量保持不变。因此,该动力学在能量景观上实现了严格的 梯度下降:网络状态不可逆地向最近的局部极小值滑落,并会无限期停留在那里。

作为能量极小值的记忆#

Hopfield 网络最初被提出作为 联想记忆 模型。通过 Hebbian 学习规则存储一组模式 {ξ1,ξ2,,ξP}\{\boldsymbol{\xi}^1, \boldsymbol{\xi}^2, \ldots, \boldsymbol{\xi}^P\}

wij=1Nμ=1Pξiμξjμw_{ij} = \frac{1}{N} \sum_{\mu=1}^P \xi_i^\mu \xi_j^\mu

每个存储模式都会成为能量函数的一个局部极小值。当以某个记忆的损坏版本或部分版本作为初始状态时,确定性动力学会自然向下流动并恢复完整模式。这就是 内容寻址记忆:检索由与存储模式的相似性驱动,而不是由显式地址驱动。

伪极小值的困扰#

然而,Hebbian 学习规则并不完美。权重矩阵 W\mathbf{W} 不仅包含来自存储模式的预期贡献,也包含不同记忆之间重叠导致的非预期 串扰 项。这些串扰项会在能量景观中创造额外局部极小值:这些状态在确定性动力学下是稳定的,却并不对应任何存储模式。这样的状态称为 伪极小值(或伪记忆)。

伪极小值有几种常见形式:

  • 混合状态:奇数个存储模式的线性组合(例如 sign(ξ1+ξ2+ξ3) \text{sign}(\xi^1 + \xi^2 + \xi^3))常常会形成稳定吸引子。

  • 自旋反转状态:对于每个存储模式,它的全局取反(所有自旋翻转)也可能成为能量极小值。

  • 随机谷地:即使是随机权重矩阵,也会仅因淬火无序而产生大量局部极小值。

随着存储模式数 PP 相对于神经元数 NN 增加,伪极小值密度会灾难性增长。当超过约 0.14N0.14 N 的临界 存储容量 后,网络会经历相变:目标记忆不再稳定,能量景观由伪吸引子主导,检索完全失败。

引入噪声:随机动力学与热涨落#

玻尔兹曼机通过引入 热噪声 解决这一困境。它不再把神经元确定性地设为输入符号,而是按概率采样状态:

P(xi=1)=σ(jwijxj+biT)P(x_i = 1) = \sigma\left( \frac{\sum_j w_{ij} x_j + b_i}{T} \right)

其中 σ(z)=1/(1+exp(z))\sigma(z) = 1/(1+\exp(-z)) 是 sigmoid 函数,T>0T > 0温度 参数。

在有限温度下,系统转移到更高能量状态的概率不为零。具体来说,仅在单元 ii 上不同的两个状态,其概率比满足:

P(xi=1)P(xi=0)=exp(ΔEiT)\frac{P(x_i = 1)}{P(x_i = 0)} = \exp\left( \frac{\Delta E_i}{T} \right)

其中 ΔEi=Exi=0Exi=1\Delta E_i = E_{x_i=0} - E_{x_i=1} 是单元 ii 打开时的能量下降量。如果 ΔEi\Delta E_i 为正(打开会降低能量),概率比偏向打开状态。关键是,即使 ΔEi\Delta E_i 为负(打开会增加能量),打开的概率仍然大于零。热涨落会偶尔把系统推向 上坡 方向。

这种攀越能垒的能力,正是逃离伪极小值的关键。困在浅层局部极小值中的网络,经过一连串不利的热涨落后,可能越过周围能垒并下降到更深、更有利的吸引域。噪声提供了确定性动力学所缺少的 探索 机制。

统计物理基础:作为热平衡的噪声#

引入噪声并不只是算法技巧;它直接来自 1.2 玻尔兹曼分布与平衡 中推导的 正则系综。当系统与温度为 TT 的热源发生热接触时,占据任意微观状态的概率由玻尔兹曼分布给出。单个神经元的随机更新规则,正是从玻尔兹曼分布推导出的条件概率 P(xi=1xi)P(x_i = 1 \mid \mathbf{x}_{-i})

P(xi=1xi)=exp(E(xi=1,xi)/T)exp(E(xi=1,xi)/T)+exp(E(xi=0,xi)/T)=σ(jwijxj+biT)P(x_i = 1 \mid \mathbf{x}_{-i}) = \frac{\exp(-E(x_i=1, \mathbf{x}_{-i})/T)}{\exp(-E(x_i=1, \mathbf{x}_{-i})/T) + \exp(-E(x_i=0, \mathbf{x}_{-i})/T)} = \sigma\left( \frac{\sum_j w_{ij} x_j + b_i}{T} \right)

因此,噪声 不是 启发式补丁,而是系统与环境处于热平衡时的物理表现。向上爬升的概率恰好由玻尔兹曼因子比值决定,从而保证系统满足细致平衡并收敛到正确的平衡分布。

模拟退火:从探索到利用#

不过,仅有噪声并不是万能药。如果温度始终很高,网络永远不会稳定下来,而是在能量景观中无休止游走,近似均匀地采样各种构型。为了得到有意义的低能状态,必须随时间逐渐降低温度,这一过程称为 模拟退火

退火调度通常从高温 ThighT_{\text{high}} 开始,此时系统容易逃离局部极小值并探索能量景观的全局结构。随后按照冷却计划缓慢降低温度,常见的是几何降温:

Tk+1=αTk,0<α<1T_{k+1} = \alpha T_k, \quad 0 < \alpha < 1

T0T \to 0 时,随机动力学接近确定性 Hopfield 极限,网络冻结到最近的局部极小值。如果退火足够缓慢,系统将以高概率收敛到 全局最小值,这一点由马尔可夫链蒙特卡洛(MCMC)理论形式化保证。

实践中,由于低温下混合时间会指数级变慢(临界慢化),大型网络几乎无法实现真正的全局优化。即便如此,模拟退火仍然是寻找深层低能状态的强大启发式方法,这些状态对应模型学习分布下的高概率构型。

噪声在学习中的双重作用#

在玻尔兹曼机学习语境中,噪声还有第二个同样重要的作用。学习算法依赖于比较网络运行的两个阶段:

  • 正相(固定):可见单元固定为某个训练样本,隐藏单元允许涨落。噪声保证隐藏单元探索与可见数据一致的所有构型,从而为依赖数据的统计量提供无偏估计。

  • 负相(自由运行):整个网络在没有外部输入的情况下自由运行。噪声驱动系统走向自身的平衡分布,生成反映模型当前“信念”的样本。

如果没有随机性,网络在两个阶段都会简单坍缩到最近的确定性吸引子,既无法捕捉隐藏原因的完整分布,也无法生成多样化样本。噪声不是需要消除的缺陷,而是使玻尔兹曼机能够表示和学习复杂概率分布的必要特性。

小结:从确定性陷阱到随机自由#

确定性 Hopfield 网络表明,对称循环网络可以把记忆存储为能量极小值,但其贪婪下降动力学使它脆弱且容易陷入伪状态。玻尔兹曼机继承了基于能量的架构,但用 受温度控制的随机采样 取代了确定性更新。

这一单一修改带来两种根本能力:

  1. 推断时逃离局部极小值,使网络能够为模糊输入找到全局一致的解释。

  2. 学习时探索完整状态空间,在配分函数不可处理的情况下提供梯度估计所需的统计样本。

后续章节将形式化玻尔兹曼机使用的能量函数,并介绍使随机采样在计算上可行的具体架构,尤其是受限玻尔兹曼机(RBM)。