[发明专利]一种多维多变量非高斯空间随机场模拟方法

摘要

本发明公开了一种多维多变量非高斯空间随机场模拟方法，包括如下步骤首先给定一个目标多维非高斯功率谱矩阵与参数的分布特性，并假定一个初始的多维高斯功率谱矩阵，求解多维非高斯功率谱矩阵，比较其与目标功率谱矩阵的误差；其次对目标多维功率谱矩阵、多维高斯功率谱矩阵以及多维非高斯功率谱矩阵进行分解，并对分解后矩阵的元素进行迭代；求解新的多维高斯功率谱矩阵；对各变量进行随机排序，重复迭代，使精度满足要求；最后由迭代得到的高斯功率谱矩阵模拟高斯空间随机场，再进一步得到多维多变量非高斯空间随机场。本发明模拟的多维多变量非高斯空间随机场精度高、迭代收敛性强，且可以与FFT联合使用，从而保证了模拟效率，适于推广应用。

主权项

一种多维多变量非高斯空间随机场模拟方法，其特征在于，包括如下步骤：(A)通过参数的相关距离，建立参数之间的目标多维非高斯功率谱矩阵，给定一个目标多维非高斯功率谱矩阵与参数的分布特性；(B)假定一个初始的多维高斯功率谱矩阵为目标多维非高斯功率谱矩阵；(C)由上述多维高斯功率谱矩阵求解多维非高斯功率谱矩阵；(D)比较多维非高斯功率谱矩阵与目标多维非高斯功率谱矩阵之间的误差；(E)对目标多维非高斯功率谱矩阵、多维高斯功率谱矩阵以及多维非高斯功率谱矩阵进行分解；(F)对分解后矩阵的元素采用指数型迭代公式进行迭代；(G)通过对迭代后得到的元素求解新的多维高斯功率谱矩阵，并对矩阵元素进行标准化；(H)对各个变量进行随机排序；重复步骤(C)至步骤(H)，直到误差矩阵的平均误差大于前一次迭代的平均误差，得到潜在多维高斯功率谱矩阵；(I)使用迭代得到的潜在多维高斯功率谱矩阵模拟多维多变量高斯空间随机场，再通过非线性变换，得到所要模拟的多维多变量非高斯空间随机场。