1.一种适用于背景谐波电压变化的谐波责任计算方法，其步骤为：A、数据采集

采样得到时刻n的母线即公共连接点的谐波电压v(n)及馈线k上的谐波电流Ik(n)，将n＝1～N时刻的谐波电压v(n)组成母线谐波电压向量V，V＝[v(1),v(2),…v(n)…Tv(N)]；同时，将馈线k上n＝1～N时刻的谐波电流ik(n)组成馈线k的谐波电流向量Ik，TIk＝[ik(1),ik(2),…ik(n)…ik(N)]，再将所有馈线上的谐波电流向量Ik组成馈线谐波电流矩阵I，I＝[I1,I2,…Ik…,IK]，k为馈线的序号，K为馈线的条数，T为矩阵的转置；

B、数据预处理及谐波责任计算

B1、背景谐波电压向量的计算：根据母线谐波电压向量V、馈线谐波电流矩阵I，利用主导波动量法求得系统谐波阻抗Zs，再求得背景谐波电压向量V0，V0＝[v0(1),v0(2),…Tv0(n)…v0(N)]＝V-ZsIpcc，其中Ipcc为公共连接点电流向量，等于各条馈线电流向量Ik之和，B2、采用均值漂移算法对背景谐波电压向量进行分段处理：B21、随机选取背景谐波电压向量V0中的一个元素作为中心v′ 0；

B22、选出满足 的所有元素，计算这些元素的均值 其中h为分段参数，其取值范围为5％～15％；

B23、如 则将这些元素的均值 作为新的中心，返回步骤B22；如则将这些元素按时刻n从小到大排列，组成背景谐波电压向量段V0,j，V0,j＝[v0,j1,v0,j2...v0,jl...v0,jL]，j为背景向量段的序号，j＝1～J，jl为背景谐波电压向量段V0,j元素的序号，jL为第j个向量段的元素个数；返回步骤B21，直至所有数据全部提取；

B24、将第j段背景谐波电压向量段V0,j中的背景谐波电压v0,l所对应的母线谐波电压v(n)，组成第j段谐波电压向量段V(j)；将第j段背景谐波电压向量段V0,j中的背景谐波电压v0,l所对应的各条馈线的谐波电流ik(n)，组成各条馈线的谐波电流向量段Ik(j)，进而组成谐波电流子矩阵I(j)＝[I1(j),I2(j),…Ik(j)…IK(j)]；

B3、谐波责任计算：根据提取的谐波电压向量段V(j)，谐波电流子矩阵I(j)，利用偏最小二乘法得到馈线k的第j段谐波责任uk(j)；加权求和得到1～N时刻的测试时间段内馈线k的谐波责任uk， 其中ω(j)为第j段的权重，。

2.如权利1所述的一种适用于背景谐波电压变化的谐波责任计算方法，其特征在于：所述的步骤B1中计算系统谐波阻抗的主导波动量方法的具体操作是：B11、计算公共连接点时刻n的谐波电流波动量Δipcc(n)，Δipcc(n)＝ipcc(n)-ipcc(n-1) (1)再计算出1-N时刻的谐波电流波动量的均值 和方差

选出满足下式的谐波电流波动量Δipcc，组成计算谐波阻抗的谐波电流波动量向量ΔIpcc，其中α为奈尔系数，取值为1～1.5；

B12、将式(1)～(4)中的Δipcc(n)替换为谐波电压波动量Δv(n)，ipcc(n)替换为谐波电压v(n)，即可得到组成计算谐波阻抗的谐波电压波动量向量ΔV；

B13、计算得出系统谐波阻抗Zs，