经典的卡尔曼滤波算法.pdf

-.自适应卡尔曼滤波卡尔曼滤波发散的原因如果卡尔曼滤波是稳定的，随着滤波的推进，卡尔曼滤波估计的精度应该越来越高，滤波误差方差阵也应趋于稳定值或有界值。但在实际应用中，随着量测值数目的增加，由于估计误差的均值和估计误差协方差可能越来越大，使滤波逐渐失去准确估计的作用，这种现象称为卡尔曼滤波发散。引起滤波器发散的主要原因有两点：（1）描述系统动力学特性的数学模型和噪声估计模型不准确，不能直接真实地反映物理过程，使得模型与获得的量测值不匹配而导致滤波发散。这种由于模型建立过于粗糙或失真所引起的发散称为滤波发散。（2）由于卡尔曼滤波是递推过程，随着滤波步数的增加，舍入误差将逐渐积累。如果计算机字长不够长，这种积累误差很有可能使估计误差方差阵失去非负定性甚至失去对称性，使滤波增益矩阵逐渐失去合适的加权作用而导致发散。这种由于计算舍入误差所引起的发散称为计算发散。针对上述卡尔曼滤波发散的原因，目前已经出现了几种有效抑制滤波发散的方法，常用的有衰减记忆滤波、限定记忆滤波、扩充状态滤波、有限下界滤波、平方根滤波、和自适应滤波等。这些方法本质上都是以牺牲滤波器的最优性为代价来抑制滤波发散，也就是说，多数都是次优滤波方法。自适应滤波在很多实际系统中，系统过程噪声方差矩阵Q和量测误差方差阵R事先是不知道的，有时甚至连状态转移矩阵或量测矩阵H也不能确切建立。如果所.-可修遍--.建立的模型与实际模型不符可能回引起滤波发散。自适应滤波就是这样一种具有抑制滤波发散作用的滤波方法。在滤波过程中，自适应滤波一方面利用量测值修正预测值，同时也对未知的或不确切的系统模型参数和噪声统计参数进行估计修正。自适应滤波的方法很多，包括贝叶斯法、极大似然法、相关法与协方差匹配法，其中最基本也是最重要的是相关法，而相关法可分为输出相关法和新息相关法。在这里只讨论系统模型参数已知，而噪声统计参数Q和R未知情况下的自适应滤波。由于Q和R等参数最终是通过增益矩阵K影响滤波值的，因此进行自适应滤波时，也可以不去估计Q和R等参数而直接根据量测数据调整K就可以了。输出相关法自适应滤波的基本途径就是根据量测数据估计出输出函数序列{C}，再由{C}推算出最佳增益矩阵K，使得增益矩阵K不断地与实际量测数kk据{C}相适应。k.Sage-Husa自适应卡尔曼滤波是在利用量测数据进行递推滤波时，通过时变噪声估计估值器，实时估计和修正系统噪声和量测噪声的统计特性，从而达到降低系统模型误差、抑制滤波发散提高哦滤波精度的目的。xxwkk,k1k1kzHxvkkkk.-可修遍--.E(w)q,E(wwT)QkkkkkkjE(v)r,E(vvT)RkkkkkkjE(wvT)0kjSage-Husa自适应卡尔曼滤波算法可描述为xˆxˆK~zkk,k1kkxˆxˆqˆk,k1k,k1k1k1~zzHxˆrˆkkkk,k1kKPHT[HPHTRˆ]1kk,k1kkk,k1kkPPTQˆk,k1k,k1k1k,k1k1P(IKH)Pkkkk,k1.其中，rˆ、Rˆ、qˆ和Qˆ由以下时变噪声统计估值器获得：kKkkrˆ(1d)rˆd(zHxˆ)k1kkkk1k1k1,k~~Rˆ(1d)Rˆd(zzTHPHT)k1kkkk1k1k1k1,kk1qˆ(1d)qˆd(xxˆ)k1kkkk1k1,kk~~.Qˆ(1d)Qˆd(KzzTKTPPT)k1kkkk1k1k1k1k1k1,kkk1,k1b式中：d，0b1为遗忘因子。k1bk1如果系统状态变量的维数比较高，而Sage-Husa自适应滤波算法中又增加了对系统噪声统计特性的计算，计算量将大大增加，实时性也将难以得到保证。除此之外，对于阶次较高的系统，Sage-Husa自适应滤波算法中R和Q的在线kk估计有时会由于计算发散失去半正定性和正定性而出现滤波发散现象，此时Sage-Husa自适应滤波算法的稳定性和收敛性不能完全保证。.-可修遍--.基于极大似然准则的自适应卡尔曼滤波，通过系统状态方差阵和量测噪声方差阵实时估计系统噪声统计特性的变化，以保证滤波器更好地适应这种变化。极大似然估计从系统量测量出现概率最大的角度估计，其特点是不仅考虑新息的变化，而且考虑新息协方差矩阵C的变化。它的量测噪声协方差矩阵Rˆ和系统噪vk声协方差矩阵Qˆ为：RˆCˆHPHTkvkkk,k1k1kQˆxxTPTPkNiikk,k1k1k,k1ikN1xxˆxˆKvkkk,k1kk1kCvvTvkNiiikN1式中：vzzˆ，N为平滑窗口的宽度。