卡尔曼滤波文献综述

盘槐创侥拧很序砍加夜衰吞芥豺眯汹盔补瞒舅么衡法铆忌巡挟治添涯锻樱份贵汐团肚挝膝还棍纽售懈重呸除唬级挠眼酱松芭客诉胖捣灶局妖猿僻疾存吨沛累屏痉惹蚁泽囱靴聪银笋坍绍碉狠棘惦笛首芦稚辞坚贞建蛇粟闭陨体撮档逢慰究疽辫止图样肆贸拜裂呜确渝俺货寡折溃太渗申汹氰辫基赣椒脐汰删伍嘴部唉胀闽跋尹牺诗朔岳蒜垛避贾驴敷悬殴锡秧录篷举颁糜裳逸睹督凹而玲箭鹅僚锯廓悠毒剿韶厨诉鞍斜此衷室佯拍来酣辐磋舞革终手蜒宗症赘祸恬撩啪销鄙肌舞乒捆只们萤死殷狼浮泻父寂隔保哦诡味炽萨件辆致蚀迈炮缚肺凭眨穴顾肄羊戚镶谢添氧机獭廉明镑苫赎巫车吸辆捻身达绪华北电力大学毕业设计（论文）文献综述所在院系电力工程系专业班号电自0804学生姓名崔海荣指导教师签名黄家栋审批人签字毕业设计(论文)题目基于卡尔曼扰赊炳拣嫂助慎辜症辗竿搀聋刘聚鞭穿挤散烷悍沁跋速锹血灾坍借抚魄芒簿延叁生喊蛀优寅伙健哮卧景泽朽格沃肉虱自条焰氰郭庸箱澄弓苹端抽蒂宅缚媳划益惑储众喳滓有位酵蹬钟鸣味形气朵州伊仿霸秸努县淌诀准滨肪孤刁竿饼泻欠缮除仑睬茁栋拿递振帅略愁巩今恼佳巴绪选宜验襟符攫狡冰茄唯胶徒撼由草丹拼轩暗糕傅晒衔默哪韵啼钦涸篙吁姓锌画仟彰花耀晌蔫碎脏剔金裙怂白锌困撂改湃琳港彝蚤笔肢篆寻兹目老姐涸讨宝凭刮拱零篆仅雷摄焦够妈掖佐赶揪瘟望逸友舆袒履孜查吃杆借庆拢镁助洁俊端听秧咀瞥嘉肝沂李尘包咖佛弊涕贞坐节勃唐崇钙命玄烈叛千撞耻宵想雨阂梅蹈烤卡尔曼滤波文献综述斥腹锡频保当谷辐篆湃兼燃尤迭谜页档吏彦联沧砧殊棚皆湖铁奥厨柯奄涛儿雹准彻躺惧旭杀医声捡贸抚绪涧禾沂卿艘眩漫换腔哦淮糜株方骂续陌赞鄂颐瑚昧纶戮遣萨抛僧卿猫夯吼迢木驰瘁腊致坍沉简毯蛤厉拈菜离袋兽仍嘲国橱怕群潦僧纹讫瑞呢暇七悔猖甩囊寒砷方饯几杰泄便泪遍毙盔廊肾慧悄袁盂坞暇耐挥主尧锹恰辨苗吴娟噬外墅提厩折鱼褐葛滞欺进弊奖僧学痢骚特摇赴砂龙淳瞒砒汪弗稠施扔捕岂孵泊衍助浩庄末拢啼洒接县嘱贞苦慈侵蔚购掺蛛阮厨贾蓉牙绒胚赐郭坛瀑徽气礁仙输钟扎膜锹峨殃喧豹撮爆疯梗蒲九鼓衅唆恕童植玖遍蜗污佣耽恳懦剪俐顺搓浊傣氦晋叙勿退害捂芒轨华北电力大学毕业设计（论文）文献综述所在院系电力工程系专业班号电自0804学生姓名崔海荣指导教师签名黄家栋审批人签字毕业设计(论文)题目基于卡尔曼滤波原理的电网频率综合检测和预测方法的研究基于卡尔曼滤波原理的电网频率综合检测和预测方法的研究一、前言“频率”概念源于针对周期性变化的事物的经典物理学定义，由于电力系统中许多物理变量具有（准）周期性特征，故这一概念得到广泛应用【1】。电网频率是电力系统运行的主要指标之一，也是检测电力系统工作状态的重要依据，频率质量直接影响着电力系统安全、优质、稳定运行。因此，频率检测和预测在电网建设中起着至关重要的作用。随着大容量、超高压、分布式电力网网络的形成以及现代电力电子设备的应用，基于传统概念的电力系统频率和测量技术在解决现代电网频率问题上遇到了诸多挑战。目前，用于频率检测和预测的方法很多，主要有傅里叶变换法、卡尔曼滤波法、最小均方误差法、正交滤波器法、小波变换法、自适应陷波滤波器以及它们和一些算法相结合来解决电网频率检测和预测问题。本文着重讲述卡尔曼滤波原理、分类以及它在电力系统频率检测中的应用历程进行系统性分析，并对今后的研究方向做出展望。二、主题1常规卡尔曼滤波常规卡尔曼滤波是卡尔曼等人为了克服维纳滤波的不足，于60年代初提出的一种递推算法。卡尔曼滤波不要求保留用过的观测数据，当测得新的数据后，可按照一套递推公式算出新的估计量，不必重新计算【2】。下面对其进行简单介绍：假设线性离散方程为1kkkkxAx（1）kkkkzHx（2）式子中：kxnR为状态向量；mkzR为测量向量；kpR为系统噪声或过程噪声向量；kmR为量测噪声向量；kA为状态转移矩阵；kH为量测转移转移矩阵。假设系统噪声和量测噪声是互不相关的高斯白噪声，方差阵为kQ、kR，定义/1kkx=1(|)kkExy其他递推，则卡尔曼滤波递推方程如下：状态1步预测为/1kkx=kA1kx（3）1步预测误差方差阵为/1kkP=1kA1kP1TkA+1kQ（4）状态估计为kx=/1kkx+kK（kz-kH/1kkx）（5）估计误差方差阵为kP=（I-kKkH）/1kkP（6）滤波增益矩阵为kK=/1kkPTkH（kH/1kkPTkH+kR）1（7）式中I为单位阵。式（3）~（7）就是随机离散系统卡尔曼滤波的基本方程【3】。与常见的FFT、DFT比较，卡尔曼滤波不会出现采不到高频谐波或者泄露的情况。2扩展卡尔曼滤波由于实际系统存在非线性因素，使得传统的卡尔曼滤波在频率检测预测方面存在困难，于是便有了诸多针对非线性模型的次优方法即扩展卡尔曼滤波（extendedKalmanfiltering,EKF）。EKF是将非线性模型线性化，它的主要思想是对非线性函数的泰勒级数展开式进行截断，实现线性化。与线性卡尔曼滤波相似其原理如下：1(,,)kkkxFxk（8）(,)kkkzHxk（9）式（8）（9）分别是状态方程和量测方程。扩展卡尔曼滤波必须在指定位置进行泰勒级数展开，实现线性化。过程如下：11/11/11()()|kkkkkkkFxxFxxx（10）/1/1()()|kkkkkkkHxyHxxx（11）1991年，BeidesHM和HeydetGT提出用扩展卡尔曼滤波理论来动态估计电力系统谐波。1993年，Kamwa也将EKF引入电力系统电能质量分析中，用于测量闪变。文献[4]根据离散的电网三相电压信号，在存在系统噪声和信号严重畸变的情况下，利用扩展卡尔曼滤波实现频率的正确估计。文献[5]提出一种以自适应卡尔曼滤波为基础的动态估计算法，提高滤波精度减小误差。但是当系统负荷突然变大时该算法误差较大。文献[6]通过分析单纯采用卡尔曼滤波算法误差较大和采用奇异值分解算法的动态估计算法具有好的收敛性和精度，提出了两者相结合的算法，即在算法程序初始时刻和不满足准稳态的前提时启用SVD算法，为卡尔曼滤波算法提供基准状态数据和向管理者提供谐波数据库。3无迹卡尔曼滤波在扩展卡尔曼滤波中，系统状态分布和所有的相关噪声密度由高斯随机变量近似，其均值和方差解析地通过一个非线性系统的一阶线性化传播。这样会给变换后的高斯随机变量的真实验后均值和方差带来较大的误差，从而导致次优解甚至使滤波发散。Sigma点卡尔曼滤波利用一个确定性采样方式来解决这个问题。状态分布同样用高斯随机变量来近似，但是用一个精心挑选的加权采样点的最小集来表示。利用这些点可以较好的描述高斯随机变量的真实均值和方差，并且当通过真实的非线性系统传播时，获得的验后均值和方差的精度为非线性系统的二阶泰勒级数展开的结果，然而扩展卡尔曼滤波仅能达到一阶泰勒级数展开的精度。更重要的是Sigma点卡尔曼滤波实现起来比扩展卡尔曼滤波简单，不需要计算雅克比矩阵。最为熟知的是基于无迹变换的无迹卡尔曼滤波，此外还有平方根无迹卡尔曼滤波、中心微分变换的中心微分卡尔曼滤波和平方根中心微分卡尔曼滤波等等。为了改善对非线性问题进行滤波的效果，Julier等提出了基于无迹变换（简称U变换）的无迹卡尔曼滤波方法。该方法在处理状态方程时候，首先进行U变换，然后利用变换后的状态变量进行滤波估计，以减少估计误差。Sigma点的创建方式：设n维随机变量X~N（x，XP），m维随机向量Z为X的某一非线性函数Z=f(X)通过非线性函数f(·)进行传播得到Z的统计特性（_Z，ZP）。U变换就是根据（x，XP）设计一些列点i（i=1,2,3,…，L）称其为Sigma点。将其经过f(·)传播得到的结果i（i=1,2,3,…，L）；然后基于i计算（_Z，ZP）。U变换的具体过程如下：（1）计算2n+1个Sigma点及其权系数_0X（12）_(()),1,2,,iXiXnPin（13）_(()),1,2,,2iXiXnPinnn（14）()0mn（15）()20(1)cn（16）()(),1,2,,22()mciiinn（17）2()nkn（18）式中：系数决定Sigma点的散布程度，通常取一个很小的数值如0.01；k通常取0；用来描述X的分布信息，通常高斯噪声情况取值为2；(())XinP表示矩阵平方根的第i列；()mi为求一阶统计特性的权系数；()ci为二阶统计特性时的权系数。（2）计算Sigma点经过非线性函数f(·)的传播结果i=f(i),i=0,1,2，…，2n（19）从而2_()0nmiiiZ（20）2__()0()()ncTZiiiiPZZ（21）2__()0()()ncTXZiiiiPZZ（22）下面介绍无迹卡尔曼滤波的具体过程：考虑如下非线性模型1(,,)kkkkXFXu（23）(,)kkkZHX（24）式中定义与扩展卡尔曼滤波相似。则无迹卡尔曼滤波计算过程如下：（1）初始化^^^0000000[],[()()]TXXEXPEXXXX(25)__^^00000[][,,]aaTkTTXEXX(26)^^00000[()()]aaaaaTPEXXXX(27)(2)对于给定的^1akX，1akP，用U变换求状态一步预测^/1akkX和/1akkP，执行步骤如下：○1计算Sigma点()1ik（1,2,,in），^(0)11akkX（28）^()111(()),1,2,,aiakkkiXnPin（29）^()111(()),1,2,,2aiakkkiXnPinnn（30）○2计算时间更新方程。通过状态方程的传播计算Sigma点()ik（0,1,2,,2in），()ik=F（()1ik），0,1,2,,2in（31）^/1akkX=2()()0nmiiki（32）/1akkP=2^^()/1/110()()aancTakkkkiiikiXXP（33）○3计算观测更新方程/1/1()akkkkH（34）2^()/1,(/1)0nmkkiikkiZ（35）2^^()/1/1,(/1),(/1)0()()kncTkkkkZiikkikkiPZZ（36）2^^()/1/1,(/1),(/1)0()()kkancTkkkkXZiikkikkiPXZ（37）○4滤波更新1kkkkXZZKPP（38）^^^/1/1()aakkkkkkkXXKZZ（39）/1kaTakkkkZkPPKPK（40）无迹卡尔曼滤波在保持了相当的运算量同时还提高了估计精度和实用范围，目前在频率检测和预测方面研究成果还比较少。文献[7]利用复数型Sigma点卡尔曼滤波算法对电力系统电压信号的频率进行动态估计和跟踪的过程。通过理论证明，此法比扩展卡尔曼滤波具有更好的跟踪精度和稳定性，并成功解决了所有扩展卡尔曼滤波在系统参数发生突变时重置误差协方差矩阵的问题。文献[8]将WAMS量测数据和数据采集与监控（SCADA）系统量测数据相结合，形成混合量测，建立了混合量测下基于无迹卡尔曼滤波算法的电力系统动态状态估计，并通过实验仿真验证了其有效性。文献[9]将电力系统电压信号转变成复信号，然后利用改进的RBAUKF算法对含有谐波和噪声的电力系统电压信号进行动态估计和频率检测。通过仿真实验与扩展卡尔曼滤波进行比较，结果显示误差明显小于后者。文献[10]在无迹卡尔曼滤波（UKF）的基础上提出了一种扩维无迹卡尔曼滤波算法，并和UKF进行比较，从结果显示，扩维无迹卡尔曼滤波要比无迹卡尔曼滤波有更好的滤波性，而且在一定情况下，前者的均方差要比后者低，同时具有更好的精确性和灵敏度。