[发明专利]一种惯性器件线性系统误差系数确定方法

权利要求书

1.一种惯性器件线性系统误差系数确定方法，其特征在于，包括如下步骤：

S1、建立惯性器件线性系统方程；

S2、对S1所述惯性器件线性系统方程中的结构矩阵进行相关性检验，获得相关性检验后惯性器件线性系统方程；

S3、对S2所述相关性检验后惯性器件线性系统方程中的结构矩阵进行显著性检验，获得显著性检验后惯性器件线性系统方程；

S4、利用最小二乘法计算显著性检验后惯性器件线性系统方程中的未知参数向量估计值；所述未知参数向量估计值即为惯性器件线性系统的误差系数；

S1中所述的惯性器件线性系统方程为：

式中，为随机变量的观测向量；C为结构矩阵；为未知参数向量；ε为n维随机误差向量；

对S1所述惯性器件线性系统方程中的结构矩阵进行相关性检验包括如下步骤：

S21、设定相关系数临界值ρ_LJ、相关性简化次数p、相关性结构矩阵C_p、第一转换矩阵Q_A、第二转换矩阵Q_X；其中p＝0；C_p＝C；Q_A和Q_X均为单位矩阵；Q_A＝I_m；Q_X＝I_m；

S22、计算获得相关性结构矩阵C_p的列数为q；

S23、计算获得相关性结构矩阵C_p中最相关的两列C_ip与C_jp，同时计算获得C_ip与C_jp的相关系数为ρ_ij；其中ip和jp均为序号；

S24、比较相关系数临界值的绝对值|ρ_LJ|与S23中获得的C_ip与C_jp的相关系数的绝对值|ρ_ij|；如果|ρ_ij||ρ_LJ|，则将C_p作为相关性检验后惯性器件线性系统方程中的结构矩阵C′，转入步骤S28，否则，使相关性简化次数p的值增加1，转入S25；

S25、根据计算比例系数r_ji；

S26、根据S22中获得的列数为q，建立第一单位矩阵I_q；首先将第一单位矩阵I_q中的第jp列删除，获得第一降维矩阵Q_Ap；然后将第一单位矩阵I_q中第ip行加上第jp行的r_ji倍之和作为第一单位矩阵I_q的第ip行，然后将第一单位矩阵I_q的第jp行删除，获得第二降维矩阵Q_Xp；

S27、将C_p×Q_Ap作为相关性结构矩阵C_p、将Q_A×Q_Ap作为第一转换矩阵Q_A、将Q_Xp×Q_X作为第二转换矩阵Q_X；然后转入步骤S22；

S28、惯性器件线性系统方程中的结构矩阵进行相关性检验结束；

S23中C_ip与C_jp的相关系数ρ_ij的计算方法为：

式中，l为第一序数；

对S2所述相关性检验后惯性器件线性系统方程中的结构矩阵进行显著性检验包括如下步骤：

S31、设定显著性简化次数u、显著性结构矩阵C′_u、系数显著性临界值F_LJ；其中u＝0；C′_u的取值为相关性检验后惯性器件线性系统方程中的相关性结构矩阵C′；

S32、采用最小二乘法对进行求解，其中

为随机变量的观测向量；为未知参数组合向量，当u为0时，为为相关性检验后的未知参数向量；ε为n维随机误差向量；

S33、计算中的各分量的显著性水平，获得显著性水平最小的变量x_uj，以及x_uj对应的显著性水平F_uj；

S34、如果F_ujF_LJ，则将C′_u作为显著性检验后的结构矩阵C″，转入S35；否则将u的值增加1，x_uj＝0，将C′_u-1中的第j列删除后作为显著性结构矩阵C′_u；然后转入S32；

S35、惯性器件线性系统方程中的结构矩阵进行显著性检验结束。