[发明专利]一种电力系统谐波阻抗的测量计算方法

权利要求书

1.一种电力系统谐波阻抗的测量计算方法，其特征在于，包括下列步骤:

步骤1：采集公共连接点的母线电压瞬时值和用户注入系统的电流瞬时值，通过傅里叶变换得到谐波电压和谐波电流相量数据序列，并建立谐波电压和谐波电流相量的复回归方程；

步骤2：对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，排除异常点对回归系数的影响；

步骤3：在初步筛选基础上，求解系统谐波阻抗回归系数；求出最优的赋权回归系数，并得到谐波阻抗的值。

2.根据权利要求1所述的电力系统谐波阻抗计算方法，其特征在于，对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，所述初步筛选采用基于相对复残差排序的方法，定义复残差列向量，e=[e₁,e₂,e₃…e_n]＇，其中，指没有经过数据预处理过程回归参数的估计值，Y为谐波电压相量数据序列组成的列向量，X为谐波电流相量数据序列组成的增广矩阵，为X的共轭转置矩阵；定义相对复残差，（i=1,2，…，n），e_j为第i个复残差值，Y_j为第i个谐波电压相量，n为列向量的个数；将e_ri由小到大排序，得到相对复残差序列{e_ri}，删除处于尾部5%相对复参差e_ri所对应的谐波电压和谐波电流数据，并将余下谐波电压和谐波电流相量记为Y_95%、X_95%。

3.根据权利要求1所述的电力系统谐波阻抗的测量计算方法，其特征在于，所述求解系统谐波阻抗回归系数，采用Ramsay影响函数赋权的反复加权复最小二乘的稳健回归方法；具体步骤为：

a）选择剩余的谐波电压X_95%和谐波电流Y_95%，根据式，计算回归参数；为X_95%的共轭转置矩阵；

b）由式，计算复残差f=[f₁,f₂,f₃,…,f_m]＇,即f={f_j}，（i=1,2,3…，m，m=0.95n）；

c）由式计算Ramsay影响函数自变量，其中，=1.5med（）；med（）是的中位数；由式，计算基于Ramsay影响函数的赋权函数，其中a为调整Ramsay影响函数形状的调整因子；

d）取赋权函数的值=（i=1,2,3…，m），记权重对角阵为W=diag（），（i=1,2,3…，m，m=0.95）；

e）由加权复最小二乘计算公式：；计算出赋权的回归参数，即；

f）计算，即计算相邻两次回归参数的差值；其中j代表回归计算循环次数，其初始值为0，当j=0时，即是步骤a求出的回归参数，而是步骤e求出的回归参数；和都是复数；

g）计算和模值和，并判断max{，}<是否成立，若成立，则执行步骤h；否则，令=，返回步骤b；

h）由即可求出谐波阻抗。