基于三次样条函数的加blackman-harris窗插值fft算法 blackman-harris window interpolated fft algorithm based on cubic spline function.pdfVIP

第29卷第2期 电力自动化设备 V01．29No．2 ElectricPowerAutomation Feh．2009 2009年2月 Equipment o 基于三次样条函数的 许珉，刘凌波 (郑州大学电气工程学院，河南郑州450002) 式和复振幅的插值修正函数过于复杂。提出利用三次样条函数逼近其频率修正系数的7次多项式和复振幅 的修正函数，采用三次样条插值函数的有效形式计算频率修正系数和复振幅的修正系数，将插值FFT算法的 频率修正系数曲线分为lo段，给定n个等间距插值点，并将复振幅修正系数曲线以频率修正系数间隔O．1 分为10段，给定11个等间距插值点，分别构造出频率修正系数和复振幅修正系数的快速计算公式。公式简 单，计算量小，程序实现方便，实时性好，并且在分段处连续，分段处的计算值为精确值。仿真结果表明，该算 法计算所得幅值误差小于0．Ol％，频率误差小于0．003Hz，相位误差小于l％。 关键词：三次样条；BIackman—harris窗；插值FFT算法；频谱泄漏；谐波；电力系统 930．1 中图分类号：TM 文献标识码：A 不能实现真正意义上的同步采样。用插值FFT算 0引言 法[¨]和其他改进算法睁8]是一种有效方法。 用快速傅里叶变换(FFT)进行电力系统谐波计 加窗插值FFT算法睁14|，是一种异步采样方 算时，必须将周期电压和电流信号截断，截取一个或 法，它是以固定不变的采样频率对信号进行采样，利 多个整周期，这必然产生截断效应，即频谱泄漏和谱 用窗函数的频谱特性修正泄漏的频谱来得到信号的 实际频谱值，可以减小异步采样引起的误差。但是， 间干扰[1吨]。由于电网频率不断波动和FFT存在栅 栏效应，当采样频率不等于电网频率的整数倍时，截 截断信号时使用的窗函数不同，其频谱泄漏和谱间 取的信号长度不是整周期的，计算得到的是泄漏的 干扰对计算结果的影响也不同。频谱泄漏的主瓣可 频谱值，使得测量的谐波幅值、相位和频率有误差， 能淹没主谱附近的谐波分量，影响是短范围的rlo]。 尤其是相位误差很大，无法满足准确的谐波测量要 旁瓣的影响是长范围的n引，它相邻的基波和高次谐 波分量影响较大。选择频谱主瓣窄的窗函数可减小 求。采用硬件锁相环电路和测频调整定时器时间常 短范围的影响，选择频谱旁瓣小的窗函数，可减小频 数实现同步采样是常用方法，但由于硬件锁相环电 谱旁瓣的影响。但是，减小了频谱旁瓣的窗函数， 路内有低通滤波环节，且对基波实现同步采样而对 其频谱主瓣较宽，实际应用中常靠增加测量时间来 间谐波不一定能实现同步采样，因此对于频率快速 减小泄漏频谱主瓣的影响，同时也可以减小旁瓣的 变化的信号和需要测量信号中的分数次谐波时，并 影响，但这又影响了测量的实时性。常见的加窗插 收稿日期：2007—08—21；修回日期：2008一06一03值FFr算法有加矩形窗插值FFT、加汉宁窗插值 Currentcontrol ofmatrixconverter