m0nkey

简单的一维随机游走问题的解析解及其在电路问题中的应用

其实应该更准确的称为“First Exit Time”问题：我们考虑到一个粒子，其初始位置在X=0处。每个时刻其向左或向右走\Delta_X距离，其中\Delta_X服从均值为\mu，标准差为\sigma的正态分布，同时在坐标轴上，在-\frac{l}{2}和\frac{l}{2}处各有两个吸收壁，当粒子碰到吸收壁时，即消逝。求随着时间增加，粒子仍然未被吸收的概率密度分布。

这类问题在电路中也可能会产生应用：考虑一个伪差分的压控振荡器（VCO），其相位差会随着时间的增加而增加，同样也会受到噪声影响，当相位差超过一定范围，就可以触发一个鉴相器输出结果；或者考虑一个积分器对电压或电流采样，每次累加的电压服从一个正态分布，此时若对溢出的时间进行统计，则可以获得采样的电压的信息。

First Exit Time（FET）问题通常可以使用Fokker-Planck Equation方程来描述，这里以一个VCO比较器为例。

考虑一个伪差分VCO比较器，后面接一个鉴相器。设l表示鉴相器的死区大小，正或负的输入电压V_\mathrm{in}使累积的相位差分别向+l/2或-l/2漂移（drift）。
设\phi(t)表示\mathrm{X}和\mathrm{Y}之间的瞬态相位差，随机变量\Phi_t的概率密度p(\phi, t)的Fokker-Planck方程为

(1)   \begin{equation*} \frac {\partial} {\partial t }p(\phi, t) = -\mu\frac {\partial} {\partial \phi }p(\phi, t)+\frac{\sigma}{2} \frac {\partial^2} {\partial \phi^2 }p(\phi, t), \end{equation*}

为了简化问题的分析，我们假设鉴相器是连续时间工作的，即一旦\phi(t)超过[-l/2,+l/2]，振荡会立即停止，由此可以得到方程的边界条件为

(2)   \begin{equation*} p(\phi, 0) = \delta(\phi), \end{equation*}

(3)   \begin{equation*} p\left(-\frac{l}{2}, t\right) = p\left(\frac{l}{2}, t\right) = 0. \end{equation*}

其中 \delta(\phi) 是迪拉克Delta函数。下面给出方程(1)在边界条件(2)-(3)的解法。

一维定边界噪声时间无关的Fokker-Planck的解析解法

为了方便表示，我们将方程(1)和边界条件(2)-(3)重写为以下形式

(4)   \begin{equation*} \begin{dcases} f_t = -\mu f_x + \frac{\sigma}{2}f_{xx} (t > 0, -\frac{l}{2} < x < \frac{l}{2})\\ f|_{x=-\frac{l}{2}}=f|_{x=\frac{l}{2}} = 0\\ f|_{t=0}=\delta (x) \end{dcases} \end{equation*}

设u(x, t) = f(x-\frac{l}{2}, t)，有

(5)   \begin{equation*} \begin{dcases} u_t = -\mu u_x + \frac{\sigma}{2}u_{xx} (t > 0, 0 < x < l)\\ u|_{x=0}=u|_{x=l} = 0\\ u|_{t=0}=\delta (x-\frac{l}{2}) \end{dcases} \end{equation*}

分离变量

(6)   \begin{equation*} \frac{\frac{\sigma}{2}X^{''}-\mu X^{'}}{X}=\frac{Y^{'}}{Y}=-\lambda \end{equation*}

(7)   \begin{equation*} \begin{dcases} \frac{\sigma}{2}X^{''}-\mu X^{'}+\lambda X = 0\\ X(0)=X(l)=0 \end{dcases} \end{equation*}

可得

(8)   \begin{equation*} \begin{dcases} \frac{\sigma}{2}r^2-\mu r+\lambda = 0\\ r_{1,2}=\frac{u \pm \sqrt{\Delta}}{\sigma}, \Delta = \mu ^2-2\sigma \lambda \end{dcases} \end{equation*}

若\Delta > 0, \lambda < \frac{\mu ^2}{2\sigma}, r_1 = \frac{\mu+\sqrt{\Delta}}{\sigma}, r_2 = \frac{\mu-\sqrt{\Delta}}{\sigma}

(9)   \begin{equation*} X = C_1 e^{r_1 x}+C_2 e^{r_2 x} \end{equation*}

(10)   \begin{equation*} \begin{dcases} C_1+C_2=0\\ C_1e^{wr_1l}+C_2e^{er_2l}=0 \end{dcases} \end{equation*}

解得C_1=C_2=0，平凡解。
若\Delta = 0, \lambda = \frac{\mu ^2}{2\sigma}, r_1 = r_2 = \frac{\mu}{\sigma}

(11)   \begin{equation*} X = C_1 e^{r x}+C_2xe^{rx} \end{equation*}

解得C_1=C_2=0，平凡解。
若\Delta < 0, \lambda > \frac{\mu ^2}{2\sigma}, r_{1,2}=\frac{\mu}{\sigma}\pm i\frac{\sqrt{-\Delta}}{\sigma}, \Delta = \mu^2-2\sigma \lambda

(12)   \begin{equation*} \begin{dcases} X = e^{\frac{\mu}{\sigma}x}(C_1\cos \frac{-\Delta}{\sigma}x +C_2\sin \frac{\Delta}{\sigma}x)\\ X(0)=C_1=0\\ X(l)=C_2e^{\frac{\mu}{\sigma}l}\sin \frac{\sqrt{\Delta}}{\sigma}l \end{dcases} \end{equation*}

(13)   \begin{equation*} \begin{dcases} \frac{\sqrt{-\Delta}}{\sigma}=\frac{n\pi}{l}\\ 2\sigma\lambda = \mu^2+\frac{n^2\pi^2\sigma^2}{l^2} \end{dcases} \end{equation*}

解得

(14)   \begin{equation*} \lambda_n=\frac{\mu^2l^2+n^2\pi^2\sigma^2}{2\sigma l^2} (n=1, 2, ...) \end{equation*}

(15)   \begin{equation*} X_n(x)=C^{'}_{n}e^{\frac{\mu}{\sigma}x}\sin \frac{n\pi}{l}x \end{equation*}

(16)   \begin{equation*} \begin{dcases} Y^{'}_{n}+\lambda_n Y_n=0\\ Y_n=A_ne^{-\lambda_n t} \end{dcases} \end{equation*}

(17)   \begin{equation*} u(x, t)=\sum_{n=1}^{\infty} C_ne^{-\lambda_n t}e^{\frac{\mu}{\sigma}x}\sin \frac{n\pi}{l}x \end{equation*}

由边界条件

(18)   \begin{equation*} u(x, 0)=\sum_{n=1}^{\infty} C_ne^{\frac{\mu}{\sigma}x}\sin \frac{n\pi}{l}x=\delta(x-\frac{l}{2}) \end{equation*}

(19)   \begin{equation*} \sum_{n=1}^{\infty} C_n\sin \frac{n\pi}{l}x=e^{-\frac{\mu}{\sigma}x}\delta(x-\frac{l}{2}) \end{equation*}

利用固有函数系的正交性可得

(20)   \begin{equation*} \begin{aligned} C_n&=\frac{2}{l}\int\limits_{0}^{l}e^{-\frac{\mu}{\sigma}x}\delta(x-\frac{l}{2})\sin\frac{n\pi}{x} dx\\ &=\frac{2}{l}e^{-\frac{\mu l}{2\sigma}}\sin\frac{n\pi}{2} \end{aligned} \end{equation*}

因此，函数的解为

(21)   \begin{equation*} u(x, t)=\sum_{n=1}^{\infty}\sin\frac{n\pi}{2}\frac{2}{l}e^\frac{-\mu l}{2\sigma}e^{-\lambda_n t}e^{\frac{\mu}{\sigma}x}\sin \frac{n\pi}{l}x \end{equation*}

即

(22)   \begin{equation*} \begin{aligned} f(x,t) & =u(x+\frac{l}{2}, t)\\ &=\sum_{n=1}^{\infty}\sin\frac{n\pi}{2}\frac{2}{l}e^\frac{-\mu l}{2\sigma}e^{-\lambda_n t}e^{\frac{\mu}{\sigma}x+\frac{l}{2}}\sin \frac{n\pi}{l}(x+\frac{l}{2})\\ &=\frac{1}{l}e^{\frac{\mu}{\sigma}x}\sum_{n=1}^{\infty}e^{-\lambda(n)t}(1-(-1)^n)\cos(\frac{n\pi}{l}x) \end{aligned} \end{equation*}

其中\lambda_n=\frac{\mu^2l^2+n^2\pi^2\sigma^2}{2\sigma l^2} (n=1, 2, ...)

一个利用概率密度函数的例子：简单积分、亚稳态深度

由上一小节，我们可以得到方程(1)的解为

(23)   \begin{equation*} p(\phi,t)=\frac{2}{l}e^{\frac{\mu}{\sigma}\phi}\sum_{n=1}^{\infty}e^{-\lambda(n)t}\sin^2\left(\frac{n\pi}{2}\right)\cos\left(\frac{n\pi}{l}\phi\right), \end{equation*}

其中

(24)   \begin{equation*} \lambda(n)=\frac{\mu^2l^2+n^2\pi^2\sigma^2}{2\sigma l^2}~(n=1, 2, ...). \end{equation*}

公式(23)即描述了在存在噪声的情况下基于闭环VCO的比较器的瞬态行为。

我们可以定义当比较器在一定时间后仍然未能做出决定为一个亚稳态（Metastability）事件，即此时压控振荡器仍在振荡。
设P(k)为k个周期后发生亚稳态的概率，则即对上面的概率密度函数进行积分，有

(25)   \begin{equation*} \begin{split} P(k)=&\int_{-l/2}^{l/2}p\left(\phi, \frac{k}{f_0}\right) d\phi\\ =&4\pi\sigma^2\cosh\left(\frac{\mu l}{2\sigma}\right)\sum_{n=1}^{\infty}\left[\frac{ne^{-\lambda\left(n\right)\frac{k}{f_0}}\sin\left(\frac{n\pi}{2}\right)}{\mu^2 l^2+n^2\pi^2\sigma^2}\right]. \end{split} \end{equation*}

Fokker-Planck方程更直觉的理解

可以看出，上面的解析解还是比较复杂的，为了方便直觉理解，我们重新在没有边界条件(3)的情况下解方程(1)，那么得到的结果就非常直接明了了：

(26)   \begin{equation*} p_{\mathrm{s}}(\phi,t)=\frac{1}{\sqrt{2\pi t\sigma}}\exp\left(-\frac{\left(\phi-t\mu\right)^2}{2t\sigma}\right), \end{equation*}

这个结果，即高斯分布的概率密度函数，均值t\mu和方差t\sigma随着时间线性增加。这其实就是一个一维维纳过程（Wiener Process）。

延伸思考

  • 在添加了吸收壁的情况下，函数的概率密度函数就从维纳过程变成了复杂的无穷级数。在工程应用中，一定条件下可以用维纳过程来减少复杂度，帮助理解。
  • 如果吸收壁是可移动的，系统的模型会更复杂，但是会不会有更好的效果：更低的亚稳态概率，或者更低的噪声。

相关阅读

发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注