一种提高电力系统闭环辨识准确性的激励信号优化方法
【专利摘要】本发明公开了一种提高电力系统闭环辨识准确性的激励信号优化方法，属于电力系统辨识与控制【技术领域】。该方法采用伪随机码二进制信号或正弦信号作为初始激励信号，辨识得到初始辨识模型；基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型，根据第K-1步的辨识模型求出优化激励信号功率谱密度并转化为有限维FIR作为激励信号激励闭环电力系统，得出第K步的辨识模型及该模型相应的控制器Wk；通过第K步的辨识模型及该模型相应的控制器Wk判断该第K步的辨识模型是否满足闭环电力系统鲁棒稳定性，若是将得到的有限维FIR作为最终的优化激励信号。通过本发明方法可以得到较为准确的电力系统辨识模型，从而进一步地解决电力系统辨识模型的准确性问题。
【专利说明】一种提高电力系统闭环辨识准确性的激励信号优化方法
【技术领域】
[0001]本发明属于电力系统控制与辨识【技术领域】；特别涉及一种提高电力系统闭环辨识准确性的激励信号优化方法，该方法基于自适应离散卡尔曼滤波的优化方法得到电力系统激励信号进行优化后的优化电力系统模型，从而提高电力系统闭环辨识的准确性。
【背景技术】
[0002]系统辨识通常是指利用系统己知的和观测到的信息，构造该系统的数学模型的理论和方法，也就是按照一定的准则，在一组模型类中，选择一个与数据拟合得最好的模型。系统辨识理论和应用发展从20世纪60年代到90年代初逐渐走向成熟的关键时期。常用的系统辨识方法有最小二乘法，辅助变量，极大似然，随机逼近，预报误差等，除此之外，还有扩展的辨识方法和求解算法。
[0003]系统辨识的输入信号设计问题最早在1960年左右开始研究，在70年代该问题的研究逐渐深入，成为系统辨识领域研究的重要课题之一。常用的优化假定输入辨识信号是由白噪声通过有限维FIR滤波得到，将最优化问题转化为一个线性矩阵不等式问题，通过求解线性矩阵不等式得到有限维FIR滤波器的相关系数，再设计有限维FIR滤波器。闭环电力系统条件下的激励信号优化主要采用是迭代方法进行优化信号，考虑控制器和闭环电力系统鲁棒性稳定性，求取激励信号的优化能量谱，从而得到满足闭环系统性能的系统模型和相应的激励信号。
[0004]此外，卡尔曼滤波算法已广泛应用于信号优化和数据处理之中，它是从系统的测量参数中恢复一些动态参数的最主要方法之一，在使用卡尔曼滤波算法时会产生过程噪声协方差矩阵Q和测量噪声协方差矩阵R。自适应滤波基本思想就是在利用观测数据进行滤波的同时，不断地对未知的或不确知的系统模型参数和噪声统计特性和状态增益阵进行在线估计或修正，实现滤波器设计参数的在线改进，以缩小实际的滤波误差，提高滤波的精度，获得状态变量的最佳估计值。
【发明内容】
[0005]本发明的目的是针对现有技术的不足，提供一种提高电力系统闭环辨识准确性的激励信号优化方法，该方法基于自适应离散卡尔曼滤波的方法得到对电力系统激励信号进行优化后的电力系统优化模型，从而解决由于不完全的先验信息和不准确估计时间的因素给激励信号优化设计带来的难度，最终提高了电力系统闭环辨识的准确性。
[0006]本发明提出的提高电力系统闭环辨识准确性的激励信号的优化方法，其特征在于，该方法包括:选择一个接近真实电力系统的初始模型Gtl及其相应的初始控制器W。；
[0007]I)采用伪随机码二进制信号或正弦信号作为初始激励信号&辨识得到初始辨识
模型(4，成)；
[0008]2)基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型，然后根据第K-1步的辨识模型J求出优化激励信号功率谱密度Φ；1:0)，其中r为第K步的激励信号，K为1、2、".Κ-1、Κ、Κ+Ρ..;
[0009]3)将步骤2)中的功率谱密度转化为有限维FIR，以有限维FIR为激励信号激励闭环电力系统，得出第K步的辨识模型及该模型相应的控制器Wk;
[0010]4)通过第K步的辨识模型(it, ΑΛ.)及该模型相应的控制器Wk判断该第K步的辨
识模型^.，Αλ.:)是否满足闭环电力系统鲁棒稳定性，若是则转步骤5)，否则令Κ=Κ+1转步骤2)；
[0011]5)将步骤4)得到的有限维FIR作为最终的优化激励信号。
[0012]本发明的特点及有益效果:
[0013]本发明提出了一种提高电力系统闭环辨识准确性的激励信号优化设计方法，结合一种自适应离散卡尔曼滤波方法，得到电力系统激励信号优化设计后的电力系统优化模型，而且使过程噪声Q和测量噪声R阵在线适应当前的高动态、高机动以及环境的多变性，为解决电力系统辨识准确性问题提出了一种新的思路。本发明的理论和实施例都表明该方法优化了电力系统激励信号并得到了优化后的电力系统模型，并最终达到提高电力系统闭环辨识准确性的目的，使之更加适用于电力系统辨识与控制中。
【专利附图】
【附图说明】
[0014]图1是本发明的提 高电力系统闭环辨识准确性的激励信号优化方法流程图；
[0015]图2是本发明方法中的自适应离散卡尔曼滤波的方法具体框图。
【具体实施方式】
[0016]本发明提供了一种提高电力系统闭环辨识准确性的激励信号优化设计方法。下面是参照附图及实施例对本方法进行详细的说明。
[0017]本发明提出的一种自适应离散卡尔曼滤波的激励信号优化方法，如图1所示，该方法包括:选择一个接近真实电力系统的初始模型Gtl及其相应的初始控制器W。；
[0018]I)采用伪随机码二进制信号或正弦信号作为初始激励信号IV使用迭代辨识的方
法得到电力系统初始辨识模型((K,);可选择的常用的电力系统激励信号的类型有很多，例如:脉冲信号、阶跃信号、伪随机码二进制信号、正弦信号、余弦信号等。
[0019]2)基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型,Jk J ,
然后根据第K-1步的辨识模型求出优化激励信号功率谱密度,其中r为第K步的激励信号，K为1、2、".Κ-1、Κ、Κ+Ρ..;
[0020]步骤2)中基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型(Α ,太,)具体如图2所示，包括以下步骤:[0021]21)首先定义真实电力系统模型的最优参数GiQO和故(10的离散状态向量表达式(I):
[0022]χΜ) = [(,,α) HAk) L (Vvi⑷ Hs (/C)]'⑴
[0023]式(I) φ Gi(k) Hi(k)表示尚待估计的电力系统参数，以GiQO和_)建立电力系统转移方程如式(2)所示:
[0024]χ , (Ji + O - ' xJ (々.)+ (众)⑵
[0025]式(2)中，A是2(N-1)X2(N_1)的单位转移矩阵；A(幻是过程噪声(均值为零、协方差矩阵为Qj(k)的高斯白噪声序列)，k= 1，2，K
[0026]式(2)中，％(幻的数学期望为0，即:
[0027]Ι?\ω:{1?) \ = O
[0028]%(々)和A(k + τ)τ的数学期望为Qj (k) g δ ( τ )，如式(3)所示:
[0029]E^mj (k), ajj (I + r)r} = Qj (k).S(t)⑶
[0030]δ (τ)为单位脉冲函数，若电力系统的测量值用N-1维的向量7；#)表示，则有:
[0031]yj(k) = H{k)-x!{k) + vi(/<)⑷
[0032]式(4)中，H(k)是(N-1) X 2 (N-1)常值矩阵，且N_1维向量巧(幻是测量噪声(均值为零、协方差矩阵为Rj (k)的高斯白噪声序列)，且巧#)和巧认)互不相关，k= 1，2，K，式(4)中，巧W和+ 的数学期望`为Rj(k).δ (τ)，即
[0033]/:: {Vj V, (k + τ? } = Ri (k).δ(τ)(5)
[0034]5#)和％(& + ^.)的数学期望为0，即:
[0035]E j v.(k),ω.(k + τ)1} = O
_1 Tyj k O O OLO O "
O O I Tlj k 0L00
，
[0036]H(k)=_I^k表示第j个状态值；
00 0 L O O I Trik
_?--κ _
[0037]22)利用卡尔曼滤波方法中的测量更新的公式(6)- (10)求出第K_1步的辨识模型(4.'太的值:
【权利要求】
1.一种提高闭环电力系统辨识准确性的激励信号的优化方法，其特征在于，该方法包括:选择'~*个接近真实电力系统的初始t旲型Gci及其相应的初始控制器Wci ； .1)釆用伪随机码二进制信号或正弦信号作为初始激励信号r。，辨识得到初始辨识模型(?0,?0)； .2)基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型,?，然后根据第K-1步的辨识模型J求出优化激励信号功率谱密度Φ--?),其中r为第K步的激励信号，K为1、2、".Κ-1、Κ、Κ+Ρ..; .3)将步骤2)中的功率谱密度Φ?:(6>)转化为有限维FIR，以有限维FIR为激励信号激励闭环电力系统，得出第K步的辨识模型及该模型相应的控制器Wk; .4)通过第K步的辨识模型H)及该模型相应的控制器Wk判断该第K步的辨识模型/}v )是否满足闭环电力系统鲁棒稳定性，若是则转步骤5)，否则令K=K+1转步骤2); .5)将步骤4)得到的有限维FIR作为最终的优化激励信号。
2.如权利要求1所述的方法，所述步骤2)基于自适应离散卡尔曼滤波的优化方法求出第K-1步的辨识模型,),具体包括以下步骤: .21)首先定义真实电力系统模型的最优参数GiGO和故(10的离散状态向量表达式:
3.如权利要求1所述方法，其特征在于，所述步骤4)中通过电力系统辨识得到的模型(?Α_)及相应的控制器Wk判断电力系统模型(?Α_)是否满足电力系统闭环电力系统鲁棒稳定性具体为:若用控制器Wk控制真实系统(?,?)是稳定的，且满足I I Δ I I A 1，其中，A = T(Gu-0)0 '为加权相对模型误差，T = GoKl+Goff)-1为电力系统补灵敏度函数；根据鲁棒控制理论，将I I ΔI I I转化为Var(A) ( 1，Var(A)表示方差。
【文档编号】G05B13/00GK103488078SQ201310303686
【公开日】2014年1月1日 申请日期:2013年7月18日 优先权日:2013年7月18日
【发明者】陆超, 柳勇军, 于淼 申请人:清华大学, 南方电网科学研究院有限责任公司, 北京建筑大学