Prony算法在电力系统负荷动态模型辨识中的应用研究.docVIP

电　力　系　统　自　动　化 1997年1月　　　　　　　　　　　　　　第21卷　第1期 AutomationofElectricPowerSystems 21 Prony算法在电力系统负荷动态模型辨识中 的应用研究 苟　北　陈　陈　提兆旭 (上海交通大学电力工程系　200240) 摘　要　提出了电力系统负荷动态模型辨识的Prony算法,给出了该算法所需要的Prony模型,并对传统Prony算法进行了发展有效性和准确性。 关键词　负荷动态模型辨识　Prony模型　Prony0　引言 意义[1]。精确计算[2]。统辨识过程。目前,所有的系统辨识方法可分为两大类:频域法和时域法。其中时域法有两种应用较为广泛的方法:①元件法;②现场测量法。元件法主要是计及系统中每个负荷元件的数学模型而获得系统总的负荷模型。在系统构造较复杂的情况下,该方法存在明显困难;而现场测量法则是采用从各个母线所监测到的实际数据来确定负荷模型的结构参数。该方法直接通过实际负荷行为来确定负荷模型。根据负荷模型的建立方法又可分为:机理模型[3]和输入输出模型。 本文所采用的Prony算法属于现场测量法范畴。它能通过给定的输入信号直接确定负荷模型的参数。Prony算法在电力系统中的应用研究目前尚局限于响应信号分析和控制设计,它在电力系统负荷动态模型辨识方面的应用研究还刚刚起步。其原因是应用Prony算法需要一个较复杂的指数形式的模型(Prony模型),而负荷模型转化为相应的指数形式存在一定困难。本文选择了一个被广泛应用的负荷模型。该模型是用电压和有功功率(或无功功率)的传递函数描述的,通过适当的变换,可将该传递函数从频域形式变换成时域形式,从而获得Prony算法所需要的指数形式。因为在最小方差意义下,Prony算法可以获得对负荷动态特性曲线的最佳拟合[4],所以本文采用该算法来辨识以传递函 1996-01-24收稿。 。文献[5]对应用Prony算法和最小二乘法辨识负荷动态模型的结果进行了对比分析,发现Prony算法不但比后者具有更好的辨识精确度,而且对于振荡的负荷动态特性,具有更强的捕捉能力。然而,传统的Prony算法只能应用于输入为冲击函数的情形。在文献[5]的基础上,本文发展了Prony算法,在输入为任意函数的情况下,得出了简单明了的计算方法,并应用两个算例分别对这两种输入情况进行了数值检验,结果表明Prony算法对于电力系统负荷动态模型辨识具有很高的精确度。 1　Prony算法简介 1795年Prony提出了Prony算法。随着当今计 算机技术的发展,该算法在电力系统中得到了广泛的重视。Prony算法可以通过给定输入信号下的响应直接估计系统的频率、振荡幅值和相对相位。正是由于这一特点,该算法在电力系统的响应信号分析中已获得初步的研究成果并显示出其良好的应用 6] 前景[3、。 n阶SISO动态系统 X(t)=AX(t)+Bu(t) (111)(112) y(t)=CX(t) 假定由于某种扰动使系统初始状态达到 X(t0)=X0 之后,系统输入u(t)保持为零,则方程(111)的解为: n 　　　X(t)= ∑R i=1 i X0exp(Κit)(2) 其中　Ri为n×n阶矩阵;Κ。i为矩阵A的特征值 将(2)式代入(112)式中: ? 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 22 n 电　力　系　统　自　动　化 y(t)= ∑ i=1 diexp(Κit)(3) 或Z D=Y(7) (3)式中的n可能与(2)式中的n有所不同,这主要 由矩阵C所决定。 在采样时间tk,(3)式可变换为: ^y(k)= n 辨识目标就是要估计所有的参数di和Κi,如果 已知zi,则根据式(5),Κi即可相应得到,而zi必定是某个n次方程的根: n-1n-2n (8)z-(a1z+a2z+…+anz0)=0形成向量A1: A1=-an,-an-1,…,-a 1,1,0,…,则 A1Y=A1ZD=0 ∑d i=1 i zi k (4)(5) zi=exp(Κi?t) (9)(10) 对于每一个采样时间tk,(4)式可改写成矩阵形式: z1z1 10 z2z2 1 …zn…zn …zn 01 d 1d2 ^y(1) = ^y(2) z1 N-1 z2 N-1 d(6) 在实际测量中,采样信号包含有噪声信号,因 此,(10)式应改写成: (11)1nen)。 (),可得到方程(12): ……… N- 1 ^y(N-1)^y(ny(n+1) ^y(n)^y(n ^y(n= ^yn^y(n-1)^y(n)