[发明专利]一种牵引供电系统-区域电网并联谐波谐振频率识别方法

权利要求书

1.一种牵引供电系统-区域电网并联谐波谐振频率识别方法，其特征在于，包括，

步骤1：选择有且只有一辆机车运行的时间段，获取有效采样率大于r*H*50Hz的牵引变压器高压侧区域电网任意一母线的A相、B相和C相中任意一相的电压波形K，以及机车网侧高压电路的电压波形P；其中，H为要分析的最高谐波次数，r为最高谐波次数的倍数；

步骤2：对电压波形K和电压波形P分别进行快速傅里叶变换，得到：K的基波电压和谐波电压的方均根值时间序列，P的基波电压和谐波电压的方均根值时间序列；进而得到：K的各次谐波电压含有率时间序列P的各次谐波电压含有率时间序列其中，h为谐波次数，h＝2,...,H；

步骤3：对任一h次谐波，采用窗宽为L、滑动步长为L/2的滑动时间窗，将和分别分割为N个子时间段，获取对应的N个谐波电压含有率段；使用二元线性回归方程拟合得到谐波电压含有率段中K与P的h次谐波电压含有率之间的关系Y_i＝λ_iX_i+b_i+ε_i；其中，i为谐波电压含有率段的序号，i＝1,…,N，Y_i为K的谐波，X_i为P的谐波，λ_i为回归直线的斜率即h次谐波从P到K的谐波放大倍数，b_i为截距，ε_i为回归残差；

步骤4：对谐波电压含有率段i中K与P的h次谐波电压含有率进行线性关系显著性检验得到线性度F_i：

其中，是P的第i个谐波电压含有率段中第j个值，Y_i^j是K的第i个谐波电压含有率段中第j个值；

步骤5：选择随着F_i的增加具有一致聚集效应的λ_i，将这些λ_i的平均值作为h次谐波从P到K的谐波放大倍数；如放大倍数大于1，则判定该牵引供电系统存在h次的并联谐波谐振频率。

2.如权利要求1所述的一种牵引供电系统-区域电网并联谐波谐振频率识别方法，其特征在于，步骤5中，选择随着F_i的增加具有一致聚集效应的λ_i，其方法为：

5.1将N个F_i从大到小排序；令m为需选择的λ_i的数量与N的百分比；

5.2计算m％*N，向下取整得到M；

5.3计算前M个F_i的方差Var1与N个F_i的方差Var2的比值即

5.4若小于阈值Δ，则判定所选择的λ_i随着F_i的增加具有一致聚集效应；否则，令m＝m-1返回步骤5.2重新计算，直到满足小于阈值Δ。

3.如权利要求2所述的一种牵引供电系统-区域电网并联谐波谐振频率识别方法，其特征在于，所述5.1中，还包括，令阈值Δ＝m％。