本发明专利技术公开了一种卡尔曼滤波传感器信息融合的故障检测方法，包括以下几个步骤：步骤一：根据卡尔曼滤波理论，建立线性动态系统的状态方程和观测方程；步骤二：根据步骤一得出的观测方程，利用最小二乘方法获取状态估计和相应的均方误差阵、新息序列；步骤三：利用已知的新息序列，得到不同的渠道归一化的新息序列；并且组成m通道平行传感器的创新矩阵；步骤四：根据步骤三所得的创新矩阵，获取创新矩阵的谱范数和谱范数的均值；步骤五：对卡尔曼滤波传感器的信息融合进行故障检测；本发明专利技术采用数理统计和区间估计的方法，简化了复杂的计算，极大地提高了故障检测速度。
【技术实现步骤摘要】
本专利技术涉及一种基于标准化创新矩阵的卡尔曼滤波传感器信息融合故障检测方法，属于卡尔曼滤波器信息融合的故障检测
技术介绍
智能化、高精度、高可靠性是未来飞行器对导航系统的要求。信息融合技术可满足未来飞行器对导航系统的这些要求。在导航和控制领域中，不可避免要对卡尔曼滤波器进行有效的测试。基于此目的已经产生了一些算法并且也产生了相应测试卡尔曼滤波测试的算法技术，这样不仅能够确保故障定位和检测，还可以进行估计修正。目前，有很多卡尔曼滤波器测试的算法，应用哲学算法可以对卡尔曼滤波器不同特征标志进行故障检测。尽管卡尔曼滤波检测算法种类繁多，但是迄今为止，多通道卡尔曼滤波器修正的监测和诊断的问题还没有得到解决。在航空航天和海军导航系统中，单一的导航系统已经无法满足系统的要求，组合导航系统越来越受到高度的重视，GPS, DGPS, GL0NASS和INS系统都通过卡尔曼滤波集成不同的组合。在多传感器融合的复杂系统中，即同时从多个来源的动态系统状态向量来获取信息，联邦或平行卡尔曼滤波器用来整合不同的导航系统够能到得到和满意的效果。Kalman滤波器信息融合。系统参数和状态多通道估计算法已成熟并用于估计一个动态系统的一个数学模型，以及多个测量通道的量测。在卡尔曼滤波器的这些模型中，有效数据的联合处理，可以改善状态向量的估计精度和数据处理的可靠性。这些算法对检测Kalman滤波器信息融合的应用而言，每个估计信道都需要它自己的“故障检测器”，这样就需要很大的计算量。所以，考虑到多信道估计程序需要大量的计算，采用先前的技术实现测试Kalman滤波器信息融合就不是一个简单的问题。因此，有必要开发用于测试的多通道估计程序的简单算法，以执行实时测试滤波器没有故障情况下参量变化的先验信息。
技术实现思路
本专利技术的目的是为了简化在GPS/INS/DVL组合导航系统中卡尔曼滤波器的传感器信息融合故障检测方法，提出，本专利技术能够简单快速检测卡尔曼滤波器的传感器信息融合情况，在得到了标准化创新矩阵谱范数后，再对其进行数理统计和区间估计，从而对卡尔曼滤波器的传感器信息融合进行简单快速的故障检测。，包括以下几个步骤:步骤一:根据卡尔曼滤波理论，建立线性动态系统的状态方程和观测方程；步骤二:根据步骤一得出的观测方程，利用最小二乘方法获取状态估计和相应的均方误差阵、新息序列；步骤三:利用已知的新息序列，得到不同的渠道归一化的新息序列；并且组成m通道平行传感器的创新矩阵；步骤四:根据步骤三所得的创新矩阵，获取创新矩阵的谱范数和谱范数的均值；步骤五:对卡尔曼滤波传感器的信息融合进行故障检测；本专利技术的优点在于:(I)本专利技术采用平行算法同时联立所有测量通道组合为矩阵，这样不需要对每个通道分别进行检测，降低了系统的复杂程度；(2)本专利技术采用数理统计和区间估计的方法，简化了复杂的计算，极大地提高了故障检测速度。附图说明图1是GPS/INS/DVL组合导航系统的容错结构图；图2是卡尔曼 滤波器的传感器信息融合图；图3是本专利技术的方法流程图。具体实施例方式下面将结合附图和实施例对本专利技术作进一步的详细说明。如图1所示，在GPS/INS/DVL组合导航系统中(即卫星/惯导/多普勒计程仪组合导航系统)，认为惯导系统除了其固有的漂移之外不发生其它故障。当任意子系统发生故障时，对故障进行有效的检测和滤波处理。当故障不可消除时，将有故障传感器隔离，当故障消除后，自动进行系统恢复。如图2所示，在卡尔曼滤波器中有m个传感器通道，要同时对m个通道的传感器进行信息融合。本专利技术提出。基于引进的标准规范化创新矩阵的光谱范数数学期望的统计。这种方法允许数学期望和实时的创新序列方差的同时测试，并且不需要先验信息在故障时刻变化值的统计特征。使用这种方法，可对信息融合的卡尔曼滤波器进行故障检测。本专利技术是一种基于标准化创新矩阵的卡尔曼滤波传感器信息融合的故障检测方法，流程如图3所示，包括以下几个步骤:步骤一:根据卡尔曼滤波理论，建立线性动态系统的状态方程和观测方程；线性动态系统的状态方程为:χ(k+Ι) = Φ (k+1, k)x(k) +G(k+1, k) ω (k)(I)式中:x(k+l)是系统k+1时刻的η维状态向量，x(k)是系统η维状态向量，Φ (k+1, k)是系统nXn维转移矩阵；G(k+l，k)是nXn维系统噪声驱动阵，ω (k)是η维系统激励噪声序列。k表示卡尔曼滤波的更新时刻，状态x(k)是由m个测量通道所组成，所以x(k)的量测满足以下关系，第i个通道的观测方程:Zi(Ii) = Hi (k) χ (k) + V j (k)(2)式中:Zi(k)是第i个测量通道的η维测量向量^(k)是系统第i个通道的nXn维的测量矩阵；ViQO第i个通道的η维量测噪声序列。第i个通道的量测噪声v(k)的均值是零方差为Rii (k)，相关阵为/和,⑷<(./)] = RlMmD，各通道相互独立，V i O)是第 ο k^ji个通道j时刻的量测噪声序列，S(kj) = \:。步骤二:根据步骤一得出的观测方程，利用最小二乘方法获取状态估计和相应的均方误差阵、新息序列；系统的状态向量可以通过卡尔曼滤波传感器信息融合来估计，可以表示为:本文档来自技高网...
【技术保护点】
一种卡尔曼滤波传感器信息融合的故障检测方法，包括以下几个步骤：步骤一：建立线性动态系统的状态方程和观测方程；线性动态系统的状态方程为：x(k+1)＝Φ(k+1,k)x(k)+G(k+1,k)ω(k)???（1）式中：x(k+1)是系统k+1时刻的n维状态向量，x(k)是系统n维状态向量，Φ(k+1,k)是系统n×n维转移矩阵，G(k+1,k)是n×n维系统噪声驱动阵，ω(k)是n维系统激励噪声序列，k表示卡尔曼滤波的更新时刻，状态x(k)是由m个测量通道所组成，第i个通道的观测方程：zi(k)＝Hi(k)x(k)+νi(k)???（2）式中：zi(k)是第i个测量通道的n维测量向量，Hi(k)是系统第i个通道的n×n维的测量矩阵，νi(k)第i个通道的n维量测噪声序列，第i个通道的量测噪声ν(k)的均值是零方差为Rii(k)，相关阵为(k)δ(kj)，各通道相互独立，νi(j)是第i个通道j时刻的量测噪声序列，δ(kj)=0k≠j1k=j;步骤二：根据步骤一得出的观测方程，利用最小二乘方法获取状态估计和相应的均方误差阵、新息序列；系统的状态向量通过卡尔曼滤波传感器信息融合来估计，表示为：x^(k/k)=x^(k/k-1)+Σi=1mP(k/k)HiT(k)Rii-1(k)Δi(k)---(3)式中：是系统的状态估计，P(kk)表示估计均方误差的相关矩阵，Δi(k)表示i通道的新息序列，是状态一步预测：x^(k/k-1)=Φ(k,k-1)x^(k-1)---(4)Δi(k)是i通道的新息序列：Δi(k)=zi(k)-Hi(k)x^(k/k-1)---(5)估计均方误差的相关矩阵：P-1(k/k)=P-1(k/k-1)+Σi=1mHiT(k)R11-1(k)Hi(k)---(6)一步预测均方误差的相关矩阵：P(k/k?1)＝Φ(k,k?1)P(k?1/k?1)ΦT(k,k?1)+G(k,k?1)Q(k?1)GT(k,k?1)???（7）步骤三：利用已知的新息序列，得到不同的渠道归一化的新息序列，并且组成m通道平行传感器的创新矩阵；对于正常运作的卡尔曼滤波传感器信息融合，不同的通道归一化的新息序列为：Δ~i(k)=[Hi(k)P(k/k-1)HiT(k)+Rii(k)]-1/2Δi(k)---(8)其中，是第i个通道的归一化新息序列，服从N(0,1)分布；建立两种假设：γ0?Kalman滤波器工作在正常状态，γ1?估计系统发生故障；设定A是一个n×m矩阵且为m通道Kalman滤波器规范化创新矩阵，n≥2,m≥2，其中，各列为不同通道的规范化创新向量，创新矩阵为：A(k)=[Δ~1(k),Δ~2(k),...Δ~m(k)]---(9)其中，A(k)是k时刻的m通道平行传感器信息融合kalman滤波器的规范化创新矩阵，是各通道k时刻的归一化新息序列；对于假设检验γ0和γ1，矩阵AT(k)A(k)服从白噪声分布；步骤四：根据步骤三所得的创新矩阵，获取创新矩阵的谱范数和谱范数的均值；实数矩阵A(k)的谱范数||?||2为：||A(k)||2≡max{(λi[AT(k)A(k)])1/2}---(10)式中：λi[AT(k)A(k)]是AT(k)A(k)的特征值；A(k)的谱范数均值为：E{||A(k)||2}≈||A(k)||‾2=1kΣj=1k||A(j)||2---(11)假设ri和aj代表规范化创新矩阵A的行和列，最大行?列范数为：μ≡max[||ri||2,||aj||2]???（12）其中||ri||和||aj||是向量ri和aj的范数，并且有下列不等式成立：E{μ}≤E{||A||2}≤[max(n,m)]1/2E{μ}???（13）把E{μ}用它的下界来替换：所以公式（13）改为：式中f为未知函数，并且有在n＝m→∞时函数f渐近为2且f的区间是（1,2），所以定义f的估计值为2，针对以上所述，E{||A(k)||2}有以下定义：σmax(n,m)≤E{||A||2}≤2σmax(n,m)---(16)表达式（16）表现了随机矩阵A标准差和它的二阶范数之间的关系；通过卡尔曼滤波器的归一化新息矩阵的元素aij是否服从N(0,1)分布来确定，并要满足不等式（16）；不满足不等式（16）则说明元素ai...
【技术特征摘要】
1.一种卡尔曼滤波传感器信息融合的故障检测方法，包括以下几个步骤: 步...
【专利技术属性】
技术研发人员：沈锋，宋丽杰，张桂贤，陈潇，李平敏，刘海峰，李强，徐定杰，宋金阳，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。