[发明专利]一种部分模糊度解算方法

摘要

本发明公开了一种部分模糊度解算方法，包含以下步骤：A、给定模糊度浮点解和协方差阵，对协方差阵进行下三角Cholesky分解，获得单位下三角矩阵和对角矩阵，本发明采用对角矩阵D顺序作为模糊度的精度顺序；采用Bootstrap成功率作为第一重约束，保证模糊度固定成功率足够高；采用FFRT检验作为第二重约束，保证模糊度固定解的可靠性；采用基线精度增益作为第三重约束，保证模糊度子集能够获得较高的基线解算精度；采用模糊度双频一致性检验作为最后一重约束，进一步保证固定解的正确性。相比于现有的子集选取方法，本发明可以更好地保证固定解的可靠性，同时能够获得足够高的基线解算精度，具有较好的实用价值。

主权项

1.一种部分模糊度解算方法，其特征在于，包含以下步骤：A、给定模糊度浮点解和协方差阵对进行下三角Cholesky分解，获得单位下三角矩阵L和对角矩阵D；B、采用LAMBDA算法进行整数变换，获得整数变换后的D，并对D阵对角元素作为模糊度剔除的对象；c、求取当前模糊度的Bootstrap成功率P_{s_ib}，若P_{s_ib}大于预先设置的阈值P₀，则进入下一步；否则，按照D阵顺序依次剔除模糊度分量直至满足要求，如果剔除后模糊度子集小于3则直接退出循环返回模糊度浮点解；其中，P_{s_ib}计算公式为：式中，表示连乘运算符；n为模糊度维数；Φ(·)为标准正态分布的累积函数运算符，d_i为对角矩阵D的对角线元素；D、采用基于固定失败率的Ratio检验方法求取检验阈值c，并采用SEVB搜索算法枚举出当前子集的候选模糊度，计算当前模糊度子集的Ratio值，若Ratio值大于阈值c，则进入步骤E；否则，按照对角矩阵D的顺序依次剔除模糊度分量直至满足要求，如果剔除后模糊度子集小于3则直接退出循环返回模糊度浮点解；其中，Ratio值的计算公式为：式中，||·||表示范数运算符；和分别为整数变换后模糊度子集的浮点解和协方差阵；和分别为模糊度子集二次型的最优解和次优解；E、计算当前模糊度子集的精度增益函数Gain_p，若Gain_p大于预先设置的阈值g₀，则进入步骤F；否则，根据Gain_p≥g₀反向求取子集z_Gain，返回步骤D；F、对当前模糊度子集进行双频一致性检验，若通过该检验则得到模糊度子集的固定解；否则直接采用模糊度全集的浮点解。