基于单片机的谐波分析精度讨论及定点实现.docVIP

基于单片机的谐波分析精度讨论及定点实现 　　摘要：谐波分析是电能质量测试中的重要环节，准确的测试为谐波的补偿提供了重要的依据。该文着重研究了不同步采样下谐波分析的理论，给出了C语言的仿真，并比较了不同窗函数的特点。对FFT算法做了改进，将浮点的运算改为定点的运算，使其更能符合项目的特点。 　　关键词：谐波；窗函数；A/D采样；FFT；定点 　　中图分类号：TP391文献标识码：A文章编号：1009-3044(2008)34-1782-03 　　Discussion about the Accuracy of Harmonic Analysis and Sentinel Achieve 　　ZHANG Peng, HE Qiang, LI Fei-fei 　　(School of Electrical an Information Engineering, Sichuan University, Chengdu 610065, China) 　　Abstract: Harmonic analysis is power quality test of an important part of accurate test for harmonic compensation provided an important basis. This article focuses on the non-synchronous sampling Harmonic Analysis under the theory, given the C language simulation, and compare the different characteristics of the window function. FFT algorithm to do the improvements, the floating-point operations will be designated the operator, to make it more consistent with the characteristics of the project. 　　Key words: harmonic; window function; sampling；FFT; sentinel 　　 　　1 引言 　　谐波是由于电网中非线性的负载的增多而引起的，随着电网中谐波地增加已严重影响了电能的质量，因此谐波分析成为电能质量测试的重要环节。 　　现有的测量方法大都采用的FFT算法，生产出了便携式的谐波分析仪。本文基于凌阳75F2413A型单片机，对谐波分析中的A/D采样精度问题及对于FFT算法的定点实现问题进行了分析。 　　2 谐波分析理论 　　2.1 频谱泄露与栅栏效应 　　假设输入的信号为单一频率的正弦信号，即 　　 ■ (1) 　　对其进行等间隔的采样，采样的频率为fs，则采样后的序列为 　　 ■ (2) 　　令■，则 ■，其频谱为■(3) 　　在频域为一冲击串。对输入序列进行加窗，即在时域相乘，那么在频域就是卷积，任何函数与冲击函数的卷积等效为对原函数进行移位，所以对矩形窗来说，其频域函数为 　　 ■(4) 　　其零点在■ ，■。 　　所以对单频率的正弦函数加矩形窗就是对函数(4)进行移位。对于DFT而言，由于其周期为2π，所以我们只讨论在[0，2π]的范围。由于在θ0的范围内也会有一个函数(4)，所以会对正半周的采样产生误差，这就是所谓的频谱泄露。当存在其他次数的谐波时，各谐波之间相互干扰，以至于幅度较高的谐波淹没幅度较低的谐波，造成分析上的误差。从上面的分析可以知道，函数(4)在■ 处的值为零，所以如果我们采样的频率是信号频率的整数倍的话，刚好在采样的位置函数(4)的值为0，所以存在频谱泄露不一定会产生频谱分析误差。但在实际的时候我们的采样频率很难做到是周期的整数倍。 　　对输入序列做DFT变换则相当与对傅立叶变换X(ejω)在[0，2π]上的N点等间隔采样，而这N点的采样值不一定就刚好在函数(4)的峰值，也就是谐波频率的位置，就相当于透过栅栏观赏风景，看不到全貌，这就是所谓的栅栏效应。 　　2.2 窗函数 　　对输入的信号直接采样，相当与对信号加矩形窗，而矩形窗的频谱泄露较大，会有很大的旁瓣干扰，所以要选用性能更好的窗函数。比较常用的窗有汉宁(Hanning)窗，哈明(Hamming)窗，布莱克曼(Blackman)窗，布莱克曼-哈里斯(Blackman-Harris)窗，几种窗的基本参数[1]如表1。 　　用Matlab的Wintool工具进行分析，得到下面图1的频谱图形。 　　从表格1和图形1我们