【摘要】 本发明公开了一种赛德尔式递推贝叶斯方法及其在状态估计中的应用，本发明是通过提取不同模式的特征量，实现不同模式后验概率的递推估计计算，在递推计算中，赛德尔式递推贝叶斯方法利用在本次递推中刚刚得到的后验概率来代替先验概率以实现后验概率的计算。本方法可以快速、精确地进行模式识别，对噪声具有很强的鲁棒性。基于赛德尔式递推贝叶斯方法的电网拓扑错误辨识方法或变压器分接头位置估计方法可以快速、精确地获得正确的电网运行模式或正确的变压器分接头位置，从而可以大幅度提高状态估计的合格率，进一步促进整个能量管理系统的实用化。 【专利类型】发明授权 【申请人】海南电网公司; 清华大学 【申请人类型】企业,学校 【申请人地址】570203 海南省海口市海府路34号 【申请人地区】中国 【申请人城市】海口市 【申请人区县】美兰区 【申请号】CN201210253525.7 【申请日】2012-07-23 【申请年份】2012 【公开公告号】CN102738794B 【公开公告日】2014-12-10 【公开公告年份】2014 【授权公告号】CN102738794B 【授权公告日】2014-12-10 【授权公告年份】2014.0 【IPC分类号】H02J3/00; G06K9/62 【发明人】黄良毅; 刘锋; 陈艳波; 何光宇; 梅生伟; 付艳兰 【主权项内容】1.基于赛德尔式递推贝叶斯估计的电网拓扑错误辨识方法，其具体的计算步骤如下： 步骤(1)初始化 形成所有可能的电网拓扑模式；令所有待辨识的电网拓扑模式的后验概率都相等，即令 其中“…”代表 i的取值范围为1≤i≤N， 为第i个电网拓扑模式的后验概率，N为电网拓扑模式的总数目；设置递推计数器k＝1； 步骤(2)递推计算以获得正确的电网拓扑模式 步骤(2.1)判断所有电网拓扑模式所对应的后验概率的最大值是否小于门槛值，即判断 是否成立，是则转步骤(2.2)；否则转步骤(3)；threshold为门槛值； 步骤(2.2)输入真实量测、虚拟量测以及伪量测，这些量测量形成量测方程为 z＝h(x)+τ 其中，z∈Rm为量测矢量，包括节点注入有功量测、节点注入无功量测、支路有功量测、支路无功量测以及节点电压量测；这些量测量一般取真实量测量，若真实量测量缺失的话，则用伪量测来代替，伪量测采用在上一次状态估计的结果叠加10％正态分布的随机数模式；同时在量测矢量中应加入0注入虚拟量测；x∈Rn为包括节点电压幅值和相角在内的状态变量；h(x)为状态矢量到量测矢量的非线性映射；τ～N(0, R)为量测误差； 为量测误差方差矩阵， 根据以上量测方程，通过以下修正方程进行迭代，直到状态估计收敛，可获得状态矢量估计值 xk+1＝xk+G(xk)HT(xk)R-1(z-h(xk)) 其中， 为雅可比矩阵，G(xk)＝(HT(xk)R-1H(xk))-1； 根据下式，分别独立计算所有电网拓扑模式下的残差矢量 以及误差方差矩阵Cf, i 根据下式，分别独立计算所有电网拓扑模式下的误差方差矩阵Cf, i Cf, i中的 从[HT R-1 H]-1对角元素的对应项取得，Cf, i中的 各项可从矩阵H[HT R-1 H]-1HT对角元素的对应项取得； 当进行电网拓扑错误辨识时，以上所进行的状态估计是全网状态估计； 步骤(2.3)利用以下公式计算所有电网拓扑结构所对应的初始后验概率； 其中， 为第k次递推时模式i所对应的残差矢量，其中i的取值范围为1≤i≤N； 为第k次递推时模式i的初始后验概率， 为第k-1次递推后模式i的后验概率； 为第k次递推中模式i对应残差的条件概率，即似然函数数值； 中求和号的意思是j的取值范围为1≤j≤N且j

未经允许不得转载：http://www.zhongzhencnc.com/1791779726.html