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

权利要求书

1.一种多维多变量非高斯空间随机场模拟方法，其特征在于，包括如下步骤：

(A)通过参数的相关距离，建立参数之间的目标多维非高斯功率谱矩阵，给定一个目标多维非高斯功率谱矩阵与参数的分布特性；

(B)假定一个初始的多维高斯功率谱矩阵为目标多维非高斯功率谱矩阵；

(C)由上述多维高斯功率谱矩阵求解多维非高斯功率谱矩阵；

(D)比较多维非高斯功率谱矩阵与目标多维非高斯功率谱矩阵之间的误差；

(E)对目标多维非高斯功率谱矩阵、多维高斯功率谱矩阵以及多维非高斯功率谱矩阵进行分解；

(F)对分解后矩阵的元素采用指数型迭代公式进行迭代；

(G)通过对迭代后得到的元素求解新的多维高斯功率谱矩阵，并对矩阵元素进行标准化；

(H)对各个变量进行随机排序；

重复步骤(C)至步骤(H)，直到误差矩阵的平均误差大于前一次迭代的平均误差，得到潜在多维高斯功率谱矩阵；

(I)使用迭代得到的潜在多维高斯功率谱矩阵模拟多维多变量高斯空间随机场，再通过非线性变换，得到所要模拟的多维多变量非高斯空间随机场。

2.根据权利要求1所述的多维多变量非高斯空间随机场模拟方法，其特征在于：步骤(C)中，采用多维傅里叶变化公式，由所述多维高斯功率谱矩阵求解得到高斯互相关函数再根据下述公式求解得到非高斯互相关函数：

RNjk(ξ)=∫-∞+∞∫-∞+∞Fyj-1[Φ(gj1)]Fyk-1[Φ(gk2)]φ[gj1,gk2;ρjk(ξ)]dgj1dgk2]]>

其中，φ[g_j1,g_k2；ρ_jk(ξ)]是第j个元素与第k个元素的在相关系数ρ_jk(ξ)下的联合正态分布函数，表示第j个函数概率分布函数的逆函数；

由所述非高斯互相关函数求解得到多维非高斯功率谱矩阵。