李文俊;肖辉;江维;曾林俊
【摘 要】非同步采样是造成频谱泄露和栅栏效应的根本原因,而频谱泄露和栅栏效应是影响DFT变换谐波测量精度的主要原因.因此采样方法的准确性至关重要.针对传统采样法上存在的误差问题.为此,提出了一种改进后的的双速率同步采样法,该方法在使用传统的采样间隔补偿方法抑制周期信号产生的周期误差后,对补偿后残留下来的误差进行分析,比较残留误差与定时器分辨率Td的大小关系,从而做进一步的补偿,再用均根方值算法进行处理,这样就更有效的减小了误差.最后,仿真结果验证了文中所提方法的有效性. 【期刊名称】《电测与仪表》 【年(卷),期】2018(055)022 【总页数】5页(P96-99,109)
【关键词】非同步采样;频谱泄露;栅栏效应;双速率同步采样法;残留误差 【作 者】李文俊;肖辉;江维;曾林俊
【作者单位】长沙理工大学电气与信息工程学院,长沙410114;长沙理工大学电气与信息工程学院,长沙410114;长沙理工大学电气与信息工程学院,长沙410114;长沙理工大学电气与信息工程学院,长沙410114 【正文语种】中 文 【中图分类】TM93
0 引 言
随着经济的发展和科学的进步，越来越多的电力电子设备、电弧放电性负载、气体放电灯以及其他非线性用电负载的投入使用，致使大量的谐波和间谐波电流涌入到各级电网，引起电网的电压正弦波形严重畸变，造成发电设备损坏，严重的危及了电力系统的安全稳定运行。因此，谐波和间谐波的检测有着重要的意义。 在信号的频谱分析中，快速傅里叶变换[1-3](FFT)是应用中最为广泛的一种方法。但它要求对信号进行整周期的采样，才能够得到准确的频谱，否则就会出现频谱泄露和栅栏现象[4]。加窗插值[5-7]的算法虽然可以很好的解决这个问题，但是窗函数的选择上较为复杂。为了更方便的减小频谱泄露和栅栏现象的影响，则更需要在信号采样的方法上进行改进。
通常信号的获取采用的是同步采样法，所谓同步采样法[8]指的是：假设交流电信号的周期为T，采样点数为N，采样间隔为t，严格的满足于条件T=N×t。但是由于间谐波频率的未知性与随机性，在实际中，是很难满足上述公式的条件。只能够近似的满足。同步采样按同步方法来分类，分别是：硬件同步法和软件同步法[9]。硬件同步法电路复杂，成本较高。所以我们常常选用软件同步法 。软件同步法是当前数字信号处理中常用的一种采样法，经过多年来的探索和发展，新的采样方法也层出不穷。但都得首先测出信号的周期T，在根据采样点数为N的条件下求得采样间隔t，t=T/N，但是这种测量方法存在着技术上的困难，主要是受定时器分辨率Td的影响，由于定时器的分辨率的有限性，定时器给出的采样间隔与理论上的采样间隔之间存在着截断误差[10-11]。一个采样点上的误差，积累N个点，即积累一个周期，就会产生所谓的周期误差。为了较好的减少周期误差对测量带来的影响，双速率同步采样法应运而生，并且在某些领域上得到了重用[12-14]。文献[15]率先提出了双速率同步采样法。基于这种方法，文章在此之上进行了改进，
使得周期误差[16]变得更小，测量结果更为精确。并通过Matlab仿真分析，验证了文章所提检测方法的准确性。 1 双速率同步采样法原理及其DFT模型
假设有一模拟信号x(t)，其周期为T。采样点数为N，则理想的采样间隔ts=T/N，这就实现了同步采样。可实际采样中却受到定时器的分辨率Td的影响，所以实际的采样间隔应该是一个关于Td的公式，即Ts=Td×INT[T/NTd], INT为取整符号。此时产生的截断误差显然则一个周期上产生的周期误差应为： (1)
令M=INT[e1/Td],M 但是采用的双速率同步采样法前后运用的是不 同的采样间隔，所以其离散傅里叶变换应该为: (4) 若令S=Ts/Td,可将式(4)简化为： (5) 这种采样方法由于先后用了2种采样速率，所以称为双速率同步采样法。利用这种采样方法所积累的周期误差在一个定时器分辨率Td以内，因此，它极大的抑制了周期误差。 2 改进后的双速率同步采样法 上面讨论的周期误差为e1，M=INT[e1/Td]，此时设一个参数H， H=MOD{e1/Td}，MOD为取余符号，H (2)当2×H≤Td时，则代表残留误差H已经不大于0.5Td，此时可以不进行补偿。 对残留误差进行分析，可以将原本的周期误差在一个Td的范围以内缩小到0.5Td的范围内。不需要用到昂贵的硬件材料和复杂的软件程序的情况下，在某种程度上就使得周期误差减少了一半。均方根值也称为有效值，它的计算方法是先平方，再平均，然后开方。通过改进后的双速率同步采样后，已经求得较为精确的谐波间谐波频谱，再利用均方根值算法，进行处理，使得测量结果更为精确。 为证明上述结论的有效性，通过实验得到了三个仿真图，分别是采用了一般的FFT算法、传统的双速率同步采样法和改进后的双速率同步采样法对幅值为10 V，频率为50 Hz的基波进行了检测。 如上述的频谱图所示，图1和图2都存在一定程度上的频谱泄露，而图3采用的 是改进后的双速率同步采样法，频谱泄露问题得到了很大程度上的解决，测量的精度明显的得到了提高，仿真图进一步的证实了补偿的作用。 图1 FFT算法频谱泄漏示意图Fig.1 Schematic diagram of general FFT algorithm spectrum leakage 图2 双速率同步采样法的频谱泄漏示意图Fig.2 Schematic diagram of spectral leakage of the two-speed synchronous sampling method 图3 改进后的双速率同步采样法的频谱泄漏示意图Fig.3 The improved dual-rate synchronous sampling method for the spectral leakage diagram 3 实验数据及分析 对于工频信号，设定谐波、间谐波的参数如表1所示，用MATLAB仿真得到的数据来验证文中所提方法的有效性。 表1 谐波、间谐波的参数设置Tab.1 Parameters setting of harmonics and interharmonics波型频率/Hz幅值/V间谐波252.28基波50380间谐波752.28谐波15019间谐波1751.9谐波25015.2 分别采用三种方法对谐波、间谐波进行MATLAB仿真计算。 (1)方法1：一般的FFT算法； (2)方法2：双速率的同步采样法； (3)方法3：改进后的双速率同步采样法。 表2、表3分别是三种采样方法在频率和幅值上的检测结果。对比分析可知，一般的FFT算法无法检测靠近基波的间谐波，基波和谐波的频点会将间谐波的频点淹没，而没有办法检测到间谐波的估计值。显然，一般的FFT算法对于含有间谐波的检测精度值太低。而双速率同步采样法在谐波、间谐波的测量精度上相对于的FFT算法有了较大的提高，因为双速率采样法相对于单速率的采样法能够在一定的程度上减少周期误差。而文中采用的改进后的双速率同步采样法在原有的双速率同 步采样法的补偿后又进行了一定程度上的补偿，仿真结果表明，该方法具有比方法1和方法2更高的检测精度。 表2 三种方法的频率估计结果Tab.2 Frequency estimation of the three methods波型设定值频率/Hz一般的FFT算法估计值/Hz偏差/%双速率同步采样法估计值/Hz偏差/%改进后的双速率采样法估计值/Hz偏差/%间谐波 25010024.52-1.9225.030.12基波5051.954.90649.28-1.4450.080.16间谐波75010075.230.30975.110.14谐波150152.201.432151.751.167150.020.013间谐波175176.300.725175,860.591175.140.08谐波250252.300.935251.630.652250.020.008 表3 三种方法的幅值估计结果Tab.3 Amplitude estimation result of the three methods波型设定值频率/Hz一般的FFT算法估计值/Hz偏差/%双速率同步采样法估计值/Hz偏差/%改进后的双速率采样法估计值/Hz偏差/%间谐波2.2801002.2910.4562.2820.087基波380377,800-0.588381.4500.381380.250.065间谐波2.2801002.2850.2582.2810.043谐波1918.320-3.56618.842-0.30518.997-0.015间谐波1.91.515-20.261.883-0.8941.8980.105谐波15.213.490-11.20015.8355.49315.2120.078 4 结束语 文章从对谐波、间谐波的采样方法这个角度上出发，为了提高采样精度，在传统的采样方法上提出了改进后的的双速率同步采样法，这种方法能够有效的抑制采样周期误差并且将误差控制在0.5Td以内，使得采样数据更精确，经仿真测试是有效的。它与普通的的单速率采样方法相比，很好的解决了频谱泄露带来的影响，比起传统的双速率同步采样法，精确度上又有了一定程度的提高，却只仅仅增加了不太复杂的软件处理程序，因而更显其实用性。由于只需要采样一个周期，更突出了该方法的实时性。故文中方法不失为一种很好的软件同步采样法。 参 考 文 献 【相关文献】 [1] 薛蕙, 基于FFT的高精度谐波检测算法[J]. 中国电机工程学报, 2002, 42(12): 106-110. Xue Hui. High precision harmonic detection algorithm based on FFT[J]. Journal of the Chinese electrical engineering project, 2002, 42(12): 106-110. [2] 邱宽民, 赵胜凯. DFT与FFT在实际应用时的性能比较[J]. 北方交通大学学报, 2000, 24(5): 60-62. Qiu Kuanmin, Zhao Shengkai. Performance com-parison between DFT and FFT in actual application[J]. Journal of Northern Jiaotong University, 2000, 24(5): 60-62. [3] 胡广书. 数字信号处理-理论、算法与实现[M]. 北京: 清华大学出版社(第二版), 2003. Hu Guangshu. Digital signal processing-theory, algorithm and implementation[M]. Beijing: tsinghua university press (2nd edition), 2003. [4] Hui Xue. A novel algorithm for harmonic measurement in Power system. 2002. Proceed-ings.Power on 2002[J]. International Conference on Power system.Technology. Oct. 2002: 438-442. [5] 祁才君, 王小海, 基于插值FFT 算法的间谐波参数估计[J]. 电工技术学报, 2003, 18(1): 92-95. Qi Caijun, Wang Xiaohai. Inter harmonics estimation based on interpolation FFT algorithm[J]. Transactions of china Electrotechnical Society, 2003, 18(1): 92-95. [6] 卿柏元, 滕召胜, 高云鹏, 等. 基于 Nuttall 窗双谱线插值FFT的电力谐波分析方法[J]. 中国电机工程学报, 2008, 28(25): 153-158. Qing Baiyuan, Teng Zhaosheng, Gao Yunpeng, et al. An approach for electrical harmonic analysis based on Nuttall window double-spectrum-line interpolation FFT[J]. Proceedings of the CSEE, 2008, 28(25): 153-158. [7] 翟晓军, 周波. 一种改进的插值FFT谐波分析算法[J]. 中国电机工程学报, 2016, 36(11): 2952-2958. Zhai Xiaojun, Zhou Bo. An Improved Interpolated FFT Algorithm for Harmonic Analysis[J]. Proceedings of the CSEE, 2016, 36(11): 2952-2958. [8] 戴先中. 准同步采样及其在非正弦功率测量中的应用[J]. 仪器仪表学报, 1984, (11): 390-396. Dai Xianzhong. Quasi-synchronous sampling and its application in non-sinusoidal power measurement[J]. Instrument journal, 1984, (11): 390-396. [9] 马宏忠. 软件实现同步采样的误差分析[J]. 电工技术学报, 1996, (11): 43-47. Ma Hongqi. The software realizes the error analysis of synchronous sampling[J]. 1996, (11): 43-47. [10]周军, 李孝文, 盛艳. 双速率同步采样法在电力系统谐波测量中的应用[J]. 计量学报, 1999, 20(2): 151 155. Zhon Jun, Li Xiaowen, Sheng Yan. Application of double-speed synchronous sampling in the power system harmonic measurement[J]. Acta Metrologica Sinica, 1999, 20(2): 151 155. [11]徐震. 谐波分析仪中采样截断误差的分析及其抑制[J]. 电测与仪表, 1990, 27(2): 3-5. Xu Zhen. Analysis and suppression of sampling truncation error in harmonic analyzer[J]. Electric Measurement & Instrument, 1990, 27(2): 3-5. [12]何丕雁. 一种改进的正弦拟合时基失真估计算法[J]. 系统工程与电子技术, 2003, 25(3): 359-361. He Peiyan. An improved algorithm for the distortion estimation of a sinusoidal composite[J]. Systems engineering and electronic technology, 2003, 25(3): 359-361. [13]董伯藩, 蔡菲娜, 任国海, 等. 非正弦波电参数综合测量仪的研制[J]. 浙江大学学报: 工学版, 1998, 32(3): 279-285. Dong Bofan, Cai Feina, Ren Guohai, et al. The precise instrument of electrical parameters of non-sinu-soidal wave-forms[J]. Journal of Zhejiang University: Engineering Science, 1998, 32(3): 279-285. [14]王跃科. 现代动态测试技术[M]. 北京: 国防工业出版社, 2003. [15]方伟林, 王立功. 双速率同步采样法在交流测量中的应用[J]. 电测与仪表, 1997, 34(4): 21-23. Fang Weilin, Wang Ligong. Application of double-speed synchronous sampling in alternating current meas-urement[J]. Electrical Measurement & Instrumentation, 1997, 34(4): 21-23. [16]Xi Jiangtao, Chicheiro J F. A new algorithm for improving the accuracy of periodic signal analysis[J]. IEEE Trans,1996, IM-45(4): 827-830. 因篇幅问题不能全部显示，请点此查看更多更全内容