Question

随机模拟方法

Answer 1

2.2.1 随机模拟的直接算法

对随机问题（带有场域随机性或时域随机性）分析的直接模拟方法的基本步骤是：

（1）建立描述和刻画系统行为功能的确定性分析模型，并确定其求解方法。

（2）分析和确认基本的随机变量（随机场）及其分布函数。采用蒙特卡罗方法产生随机数（随机样本）。

（3）根据所产生的随机样本，按确定性分析方法求解所模拟问题的输出量（系统反应）。

（4）计算分析系统反应量的样本反应估计值，如样本反应均值，样本反应方差及样本反应谱密度估计等。

2.2.2 Neumann展开算法

对许多工程问题的分析计算最终都归结为计算求解下列形式的线性方程组：

地下水系统随机模拟与管理

式中：K——受随机变量影响的系统结构整体刚度矩阵；

F——由边界和系统外部条件确定的列向量；

H——系统状态反应列向量。

由于矩阵K一般具有对称正定性，所以可用Cholesky分解法求解上述方程，即取：

地下水系统随机模拟与管理

而H可由下列两步求解过程给出：

地下水系统随机模拟与管理

在随机模拟的直接算法中，每产生一次随机样本结构，就要进行一次（2.14）式的分解运算，而为了获得更加精确的随机样本的统计量，必然会有大量的随机样本产生，所以其计算工作量非常之大，而利用Neumann展开思想，可以在全部模拟计算过程中只进行Cholwsky分解，从而使计算工作量大大降低。为此设：

地下水系统随机模拟与管理

式中：K₀——均值参数结构所对应的系统整体结构矩阵；

ΔK——样本结构关于均值参数结构的偏差部分。

显然：

地下水系统随机模拟与管理

由式（2.13）及式（2.17）有：

令：

地下水系统随机模拟与管理

则由Neuman展开公式有：

地下水系统随机模拟与管理

将式（2.21）代入式（2.19）有：

地下水系统随机模拟与管理

显然可得下列递推公式：

地下水系统随机模拟与管理

将式（2.20）代入式（2.23）并稍作变换即有：

地下水系统随机模拟与管理

于是，一旦由K₀ H₀=F（2.25）

求出H₀，即可用递推公式依此求得H₁，H₂，…，H_n。代入（2.22）便得样本反应量H。

原则上，Neumann展开算法只是随机模拟算法实施过程中为节省运算工作量而采取的一项技术措施，而对于随机模拟的思想则未做任何改进。

2.2.3 摄动算法

对于确定性物理问题的控制方程可以表示为带有小参数的方程。

地下水系统随机模拟与管理

式中：L——一般的线性算子；

H——所研究物理问题的解，H一般可表示为H=H（X，ε）；

X——控制自变量；

ε——小参数。

一般地说，方程（2.26）所描述的问题往往不能精确地解出，但根据H为X和ε的函数且ε是小参数的特点，可以用ε的一个渐近展开式来表示H，即：

地下水系统随机模拟与管理

式中H_i（x）与ε无关。

将式（2.27）代入式（2.26），并将ε的同次幂系数合并起来可得：

地下水系统随机模拟与管理

式中：L₀，L₁，L₂，…——空间H中的线性算子；

h——关于x的实函数，可根据具体情况给出它们的形式。

由于方程（2.28）对所有的ε都必须成立，又因为ε的序列是线性无关的，故ε各次幂前面的系数项必须自动为零，即有：

地下水系统随机模拟与管理

式（2.29）构成了关于H_i（x）的递推方程组。

根据边初值条件可依次求得序列H₀，H₁，H₂，…从而代入式（2.27）可得H（x，ε）的近似解。

将上述关于确定性物理问题的摄动求解思想推广到带有随机参数的问题中来，就构成了随机参数摄动问题，为此，设所考虑问题的随机微分算子方程为：

地下水系统随机模拟与管理

式中：L——算子符号；

x——自变量；

ξ——某一给定分布的随机变量；

Y——一随机函数，可表示为Y=Y（x，ξ）。

随机变量 可转化为用标准随机变量表示的形式：

地下水系统随机模拟与管理

式中：ξ₀——随机变量ξ的均值；

ξ_r——随机变量ξ的均方差；

b——均值为零，方差为1的标准化随机变量。

将式（2.31）代入式（2.30）有：

地下水系统随机模拟与管理

利用随机函数的幂级数展开式可将解Y（x，ξ）展开为关于随机变量b的级数：

地下水系统随机模拟与管理

由于-b=0，故上式可表示为：

地下水系统随机模拟与管理

由于 y 未知，所以系数式等也是未知的，但可将上式等价地写为：

地下水系统随机模拟与管理

显然，上式中的系数U_i（x）与b无关，为一确定性函数。

将式（2.34）代入式（2.32）并经适当的运算后将b的不同次幂系数项合并起来，可得：

地下水系统随机模拟与管理

式中：L₀，L₁，L₂，…，L_n——确定性算子；

h——关于x的实函数。

由于随机变量b具有任意性，因此，式（2.35）成立的充分条件是各系数项皆为零，由此可得：

地下水系统随机模拟与管理

上述方程组为一组确定性算子方程，当给定边界条件和初值条件以后，便可依次求出解U₀，U₁，U₂，…回代方程（2.13）后即可得到y（x，ξ）的形式解答，而解答的均值与方差分别为：

地下水系统随机模拟与管理