详细平稳解

2024-01-07 00:12

文章标签 详细平稳

本文主要是介绍详细平稳解，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

1.详细平衡

定义：一个在高斯白噪声激励下的动力学系统在状态空间中如果用如下运动方程描述：
$\frac{d}{dt}\mathbf{X}_{j}$ = $f_{j}(\mathbf{X})$ + $\sum_{l=1}^{m}g_{jl}(\mathbf{X})W_{l}(t)$ ,j=1,2,…n
它一般不属于平稳势类。为满足FPK方程，可将一阶导数矩分成如下两个部分：
$a_{j}(\mathbf{x})=a_{j}^{R}(\mathbf{x})+a_{j}^{I}(\mathbf{x})$
其中， $a_{j}^{R}(\mathbf{x})$ 为可逆分量， $a_{j}^{I}(\mathbf{x})$ 为不可逆分量。不可逆部分和阻尼力相应，可逆部分与惯性力和恢复力相应 (对于拉格朗日提法和哈密顿提法也是如此)

补充：维纳过程和高斯白噪声之间的关系

考虑方程 $\frac{dX(t)}{dt}=W(t),X(0)=0$ 其中， $W (t)$ 是谱密度为 $K$ 的高斯白噪声，即 $E[W(t)]=0，E[W(t)W(t+\tau)]=2\pi K\delta(\tau).$
按维纳过程的定义， $X (t)$ 是维纳过程，则 $\frac{dB(t)}{dt}=W(t).(*)$ 而维纳过程的强度 $\sigma^{2}$ 与高斯白噪声的谱密度之间的关系为 $\sigma^{2}=2\pi K.$
这里的(*）式只是一种形式上的关系，因为维纳过程 $B (t)$ 在 $L_{2}$ 意义上不可微。
作为最简单的马尔可夫扩散过程，维纳过程 $B (t)$ 可以通过随机微分方程用于构造其他马尔可夫过程，一个标量马尔可夫扩散过程可由下式产生： $dX(t)=m(X,t)dt+\sigma(X,t)dB(t)$ ，其中 $B (t)$ 是单位维纳过程，即
$E[B(t_{1})B(t_{2})]=min(t_{1},t_{2}),E[dB(t_{1})dB(t_{2})]=\left\{\begin{matrix} 0, &t_{1}\neq t_{2} \\ dt, &t_{1}=t_{2}=t \end{matrix}\right.$

2.以下举例说明：

2.1外激单自由度系统

考虑系统 $\mathbf{\ddot{X} }+h(\Lambda)\mathbf{\dot{X}}+u(\mathbf{X})=W(t)$ ,式中 $u(\mathbf{X})$ 是恢复力， $W (t)$ 是谱密度为 $K$ 的高斯白噪声， $\Lambda$ 是系统的总能量，也就是 $\Lambda=\frac{1}{2}\mathbf{\dot{X}}^{2}+ \int_{0}^{X}u(z)dz.$
自然地，我们可以知道阻尼力为 $\mathbf{\dot{X}}$ 对应部分，惯性力和恢复力为 $u(\mathbf{X})$ 。

首先，我们得到系统对应的伊藤方程(令 $\mathbf{X_{1}}=\mathbf{X}$ ， $\mathbf{X_{2}}=\mathbf{\dot{X}}$ ):
$d\mathbf{X_{1}}=\mathbf{X_{2}}dt\\ d\mathbf{X_{2}}=-[h(\Lambda)\mathbf{X_{2}}+u(\mathbf{X_{1}})]dt+\sqrt{2\pi K}dB(t).$
对于这里 $d B (t) = W (t) d t$ ，由于把普通随机微分方程转化成了伊藤随机微分方程，故按对于的规则，得到上式； $\sigma(X,t)dB(t)=\sqrt{2\pi K}dB(t)$ .

由此得到一、二阶导数矩：
$a_{1}=x_{2},a_{2}=-h(\lambda)x_{2}-u(x_{1}),b_{11}=b_{12}=b_{21}=0,b_{22}=2\pi K$
其中可逆部分：
$a_{1}^{R}=x_{2},a_{2}^{R}=-u(x_{1})$
不可逆部分：
$a_{1}^{I}=0,a_{2}^{I}=-h(\lambda)x_{2}$
将上式带入 $a_{j}^{I}=\frac{1}{2}\sum_{k=1}^{n}[\frac{\partial }{\partial x_{k}}b_{jk}(x)-b_{jk}(x)\frac{\partial \phi}{\partial x_{k}}],\\ \sum_{j=1}^{n}\frac{\partial }{\partial x_{j}}a_{j}^{R}(x)=\sum_{j=1}^{n}a_{j}^{R}(x)\frac{\partial \phi}{\partial x_{j}}.$
（若系统属于详细平衡类，则存在一个 $\phi(x)$ 满足所有这些方程。)
得到， $\pi K \frac{\partial\phi}{\partial x_{2}}=-h(\lambda)x_{2},(1)\\ x_{2}\frac{\partial\phi}{\partial x_{1}}=u(x_{1})\frac{\partial \phi}{\partial x_{2}}.(2)$
解(1)式：
请添加图片描述
带入(2)式，发现 $g(x_{1})$ 必为常数，又对于一个系统属于平稳势类则其平稳概率密度表示为： $p(x)=Cexp[-\phi(x)]$ 其中，C为归一化常数，于是
$p(x_{1},x_{2})=Cexp[-\frac{1}{\pi K}\int_{0}^{\lambda} h(z)dz],\lambda=\frac{1}{2}x_{2}^{2}+\int_{0}^{x_{1}}u(z)dz.$

故对于受外激作用的单自由度系统属于详细平衡类.
特别地，对于线性阻尼力情形: $h(\Lambda)\dot{X}=\alpha \dot{X}$ , $p(x_{1},x_{2})=Cexp\{-\frac{\alpha}{\pi K}[\int_{0}^{x_{1}}u(z)dz+\frac{1}{2}x_{2}^{2}]\}.$

2.2同受外激和参激的单自由度系统

补充：外激（施加在系统上的外力或外部扰动)；参激(系统内部的反馈信号或参数输入)。
$\ddot{X}+(\alpha+\beta X^{2})\dot{X}+\omega_{0}^{2}X=XW_{1}(t)+W_{2}(t)$
其中 $W_{1}(t)$ 与 $W_{2}(t)$ 是谱密度分别为 $K_{11}$ 与 $K_{22}$ 的独立高斯白噪声，以 $X_{1}$ 记 $X$ ， $X_{2}$ 记 $\dot{X}$ ，则对于的伊藤方程：
$dX_{1}=X_{2}dt\\ dX_{2}=-[(\alpha+\beta X_{1}^{2})X_{2}+\omega_{0}^{2}X_{1}]dt+\sqrt{2 \pi (K_{11}X_{1}^{2}+K_{22})}dB(t)$