专利名称：基于svd数字滤波器的电力系统谐波分析方法
技术领域：
本发明涉及数字信号处理中的滤波器设计和电カ系统中的谐波分析领域，具体涉及一种基于SVD (singular value decomposition,奇异值分解)数字滤波器的电カ系统谐波分析技术。
背景技术：
随着电カ电子技术的快速发展，电カ系统的谐波污染日益严重，电能质量问题受到高度重视，因此有必要对电カ系统中的谐波分量进行提取以进行后续分析。一方面，发电厂和变电站的高压电能计量关系到发电、送电、用电三方的利益，这就对电能计量的准确度提出了更高的要求；另ー方面，对电网电能质量和使用情况的在线监测对于电网的安全稳 定运行和合理收费具有重要意义，可以根据监测结果进行相关分析，減少大量现场校验エ作。目前，对于含有高次谐波情况下的电信号的分析，主要有基于FIR(FiniteImpulse Response,有限冲击响应)数字滤波算法、傅立叶变换法、小波变换法、Prony算法等。基于FIR数字滤波器的电カ系统各次谐波分析及相应的电能參数计算是这些方法中原理简单且较为有效的ー种算法。然而该方法具有硬件执行复杂度(一般用算法硬件实现时所需乘法器个数来衡量)较高的缺点。基于FIR数字滤波器的电カ系统各次谐波分析算法通过设计若干符合指标要求的带通FIR滤波器滤除电カ系统信号中“不想要的信号”从而将需要分析的谐波信号提取出来以达到进行后续分析的目的。为了能从原信号中准确可靠的提取需要进行分析的谐波信号，所设计的带通FIR滤波器应满足窄通带、陡过渡带的幅度频率响应特性。因此，基于FIR数字滤波器的电カ系统各次谐波分析算法可归结为窄通帯、陡过渡带带通FIR滤波器的设计。采用传统的Parks-McClellan最优数字滤波器设计算法可以很方便的设计满足要求的窄通带、陡过渡带带通FIR滤波器。然而，该滤波器设计算法所设计的滤波器具有很高的阶数。较高的滤波器阶数意味着较高的硬件执行复杂度，而且硬件系统工作时的功耗将会较高。基于滤波器系数类周期性的低通FIR滤波器设计技术(外插法)，具有硬件执行复杂度较低的优点。然而外插法仅适用于低通FIR滤波器设计，且外插法中各參数的选取都是通过观察原型滤波器系数特点获得，这是因为低通FIR滤波器系数具有类周期性，而基于FIR滤波器的谐波分析算法所需要的窄通帯、陡过渡带带通FIR滤波器系数却不具有类周期性的特点。

发明内容
本发明所要解决的技术问题是为了減少基于FIR数字滤波器的谐波分析算法的硬件执行复杂度，提出一种基于SVD的窄通带、陡过渡带带通FIR滤波器设计算法。
本发明为解决上述技术问题采用以下技术方案一种基于SVD数字滤波器的电カ系统谐波分析方法，将奇异值分解应用到数字滤波器设计中，提取需要进行分析的电カ谐波分量；所述数字滤波器为窄通带、陡过渡带带通FIR滤波器，具体步骤如下步骤001，首先，设计得到窄通带、陡过渡带原型带通FIR滤波器系数h (η)，该原型带通FIR滤波器的长度为(2 · Nf+1)，其中I η I彡Nf, Nf为任意正整数； 步骤002，选取參数M、d和R，使得所述原型带通FIR滤波器系数h (η)组合成具有(R+1)个列向量的矩阵Η，η彡(Μ+1);其中，參数M为任意非负整数，參数d为任意正整数，參数 R= [(Nf-M)/d-1]
权利要求
1.一种基于SVD数字滤波器的电力系统谐波分析方法，其特征在于将奇异值分解应用到数字滤波器设计中，提取需要进行分析的电力谐波分量；所述数字滤波器为窄通带、陡过渡带带通FIR滤波器，具体步骤如下 步骤OOl，首先，设计得到窄通带、陡过渡带原型带通FIR滤波器系数h (η)，该原型带通FIR滤波器的长度为(2*NF+1)，其中|n|彡NF，Nf为任意正整数； 步骤002，选取参数M、d和R，使得所述原型带通FIR滤波器系数h(n)组合成具有(R+1)个列向量的矩阵H, η > (Μ+1);其中，参数M为任意非负整数,参数d为任意正整数,参数
2.根据权利要求I所述的一种基于SVD数字滤波器的电力系统谐波分析方法，其特征在于算法中的参数M、d和R按如下规则设置尽可能使得矩阵H的行向量的个数与列向量的个数接近。
3.根据权利要求I所述的一种基于SVD数字滤波器的电力系统谐波分析方法，其特征在于所述步骤OOl中设计得到窄通带、陡过渡带原型带通FIR滤波器系数h (η)的具体步骤如下 步骤Α、采集连续的电力系统信号，将其转换为数字信号X (η)，设信号采集时的采样频率为;然后根据需要确定从采集到的信号中提取谐波分量;步骤B，根据采样频率fsample和谐波分量も以及实际情况，确定窄通带、陡过渡带带通FIR滤波器的通带范围[fpl fp2]、阻带范围[O fsl]和[fs2 O. 5 · fsampJ，以及通带最大波纹值rp、阻带最大衰减值rs ； 步骤C，根据步骤B所述的采样频率fsample、带通FIR滤波器的通带范围和阻带范围，获得归ー 化带通FIR滤波器的通带范围为[2 · (fpl/fsample)2 · (fp2/fsample)],阻带范围为[O 2 · (fsl/fsample)]和[2 · (fs2/fsample) I]； 步骤D，根据步骤C归ー化后的带通FIR滤波器的通带范围、阻带范围以及通带最大波纹值rp、阻带最大衰减值rs,采用经典的Parks-McClellan最优化FIR滤波器设计方法设计得到窄通带、陡过渡带带通FIR滤波器系数h (η)。
全文摘要
本发明公开了一种基于SVD数字滤波器的电力系统谐波分析方法。本发明首次将奇异值分解应用到窄通带、陡过渡带带通FIR滤波器设计中，首先通过传统方法设计得到窄通带、陡过渡带原型带通FIR滤波器系数，然后选取合适的参数、和，使得原型带通FIR滤波器系数组合成具有个列向量的矩阵；之后对矩阵进行奇异值分解，获得原型带通FIR滤波器系数的近似表达式，根据这个近似表达式构造FIR滤波器，即可提取需要进行分析的谐波分量。试验表明本发明能够节省窄通带、陡过渡带带通FIR滤波器硬件执行所需的乘法器个数，从而使得基于FIR数字滤波器的电力系统谐波分析算法的硬件执行复杂度降低。
文档编号G01R23/167GK102721870SQ20121022001
公开日2012年10月10日 申请日期2012年6月29日 优先权日2012年6月29日
发明者何金陵, 崔恒志, 王浩, 王翀, 赵俊峰 申请人:江苏省电力公司信息通信分公司