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

权利要求书

1.一种电力系统谐波阻抗计算方法，其特征包括下列步骤：

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

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

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

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

3.根据权利要求2所述的一种电力系统谐波阻抗计算方法，其特征在于，所述求解系统谐波阻抗回归系数，采用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(|f_i|)； med(|f_i|)是|f_i|的中位数；由式，计算基于Ramsay影响函数的赋权函数，其中a为调整Ramsay影响函数形状的调整因子；

d)取赋权函数的值ω_i (i=1,2,3…m)，记权重对角阵为W=diag(ω_i),(i=1,2,…,m,m=0.95n)；

e)由加权复最小二乘计算公式：；计算出赋权的回归参数，即λ^ω=ZshU·sh]]>；

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

g)计算Δλ₁和Δλ₂模值|Δλ₁|和|Δλ₂|，并判断max{|Δλ₁|,|Δλ₂|}＜ε₀是否成立，若成立，则执行步骤h；否则，令=，返回步骤b；

h)由即可取出谐波阻抗Z_sh。