四参数正弦波曲线拟合的快速算法.pdf

四参数正弦波曲线拟合的快速算法梁志国,张大治,孙王景宇,李新良(中国一航北京长城计量测试技术研究所计量与校准国防科技重点实验室,北京1 0 0 0 9 5)摘要:介绍了一种基于信号频率估计方法,利用三参数正弦波拟合算法实现的四参数最小二乘正弦波拟合算法,特点是:1)绝对收敛;2)不需要对参数进行预估计;3)算法简洁;4)不需要迭代过程;5)运算速度高。仿真实验验证了方法的有效性、正确性和局限性。关键词:正弦波;曲线拟合;收敛性;评价中图分类号:O2 4 1.5文献标识码:A文章编号:1 0 0 2-6 0 6 1(2 0 0 6)0 1-0 0 0 4-0 4AF a s t A r i t h me t i cMe t h o do f F o u r-p a r a me t e rS i n eWa v eC u r v e-f i tL I ANGZ h i-g u o,Z HANGD a-z h i,S UN J i n g-y u,L I X i n-l i a n g(Na t i o n a l L a b o r a t o r yo f Me t r o l o g ya n dC a l i b r a t i o no f C h a n g c h e n gI n s t i t u t eo f Me t r o l o g y&Me a s u r e me n t,B e i j i n g1 0 0 0 9 5,C h i n a)A b s t r a c t:I nt h i sp a p e r,af o u r-p a r a me t e rl e a s ts q u a r es i n ew a v ec u r v e-f i tme t h o db a s e do nt h ef r e q u e n c ye s t i ma t i o nme t h o da n dr e a l i z e db yu s i n gt h et h r e e-p a r a me t e rs i n ew a v ec u r v e-f i tme t h o di si n t r o d u c e d.T h e f e a t u r e s o f i t a r e:1)a b s o l u t e c o n v e r g e n c e;2)w i t h o u t t h e p r e-e s t i ma t i o no f p a r a me t e r s;3)s i mp l e a n de a s yt ob e r e a l i z e d;4)n oi t e r a t i o n;5)h i g hs p e e dc a l c u l a t i o nc o u r s e.T h e s i mu l a t i o nr e s u l t sh a v ep r o v e dt h ec o r r e c t n e s s,r i g h t n e s sa n dt h ed e f e c t o f t h i sme t h o d.K e yw o r d s:s i n u s o i d a l w a v e;c u r v e-f i t;c o n v e r g e n c e;e v a l u a t i o n收稿日期:2 0 0 5-0 7-1 8作者简介:梁志国(1 9 6 2-),男,研究员,硕士,研究方向为数字化仪器与系统校准评价、模式识别。0概述正弦波曲线的四个独立参数分别为:幅度、频率、相位和直流分量。由正弦波形的等间距采样序列获得其四参数拟合正弦曲线函数,是一种基本信号处理方法,在许多场合下获得了应用,如评价数据采集系统的有效位数、采集速率、交流增益、通道间延迟、触发特性等 1,在调制信号的数字化解调和失真度测量中,也有应用 2-5。针对四参数正弦波曲线拟合算法,已有众多学者作过大量研究工作 6-1 6。对于这类算法,人们比较关注其以下几方面的特性 6:相对精度;绝对精度;效率;收敛性;运行时间;残差形式;鲁棒性。通过上述研究,人们发现,已知频率,关于幅度、相位和直流分量的三参数正弦波曲线拟合是一种闭合的线性过程,绝对收敛。而四参数正弦波曲线拟合则不然,尚无确切的数学公式可直接计算获得拟合参数,多数已知的方法都属非线性迭代拟合过程,若拟合初始值距离目标值太远,则很容易导致迭代过程发散或收敛到局部最优点而不是总体最优点上,致使拟合结果错误。另外,多数四参数正弦波拟合算法的共同缺点是拟合过程需要大量运算时间,从而影响了算法的效率和实时性应用。因此,提高运算效率和缩短运行时间,也一直是四参数正弦波曲线拟合的目标之一。一种显而易见的方法是使用组合方法与技术来达到目的,即不是同时使用四参数拟合方法获得最终结果,而是通过其它方法首先对信号的频率进行估计,在此基础上,再使用已知信号频率的三参数正弦拟合算法进行拟合,最终获得正弦波形的四个参数 1 2,1 6。当正弦信号采样序列是等间隔采样序列时,该拟合过程将不需要叠代运算,从而具有速度快、过程简洁的特点,但精度通常较四参数直接拟合差。本文下面的过程具有这样的特点。首先,将详细介绍一种三参数正弦波曲线拟合算法;然后介绍一种正弦频率估计方法;最后,介绍由它们的组合获得的四参数正弦波曲线拟合算法。1三参数正弦波曲线拟合法 1 5 设理想正弦信号为42 0 0 6年第2 6卷第1期理论与实践y(t)=A0c o s(2 f t)+B0s i n(2 f t)+D0=C0c o s(2 f t+0)+D0(1)数据记录序列为时刻t1,t2,tn的采集样本y1,y2,yn,采集速率v已知,采样间隔为 t,ti=i t=i/v,(i=1,n),数字角频率=2 f/v,则式(1)可表示成下列离散形式y(i)=A0c o s(i)+B0(s i n(i)+D0=C0c o s(i+0)+D0(2)三参数正弦波曲线拟合过程,即为输入信号的数字角频率已知,选取或寻找A,B,D,使下式所述残差平方和最小:=ni=1 yi-Ac o s(i)-B s i n(i)-D2(3)则,参数A,B,D即为A0,B0,D0的最小二乘拟合值。为寻找出A,B,D,构造矩阵G=c o s(H1)s i n(H1)1c o s(H2)s i n(H2)1IJKLc o s(Hn)s i n(Hn)1MIJ=y1y2yKLnNO0=IJKLABD则式(3)可用矩阵表示为=()=(M-G O0)P(M-G O0)当式(3)最小时可得O0的最小二乘解OQ0:OQ0=(GPG)-1(GPM)拟合函数为yQ(i)=Ac o s(i)+B s i n(i)+D其幅度和相位表达形式yQ(i)=C c o s(i+)+D其中C=A2+BR2=S T c U S n(-BA)N AV0S T c U S n(-BAWXY)+N AZ0拟合残差为i=yi-Ac o s(i)-B s i n(i)-D拟合残差有效值为 =RnN其中=ni=12i=ni=1(yi-yQ(i)2由于这是一种闭合算法,因而收敛是肯定的。将上述三参数正弦波曲线拟合方法与信号频率估计方法相结合,可获得绝对收敛的四参数正弦波曲线拟合方法。四参数正弦波曲线拟合的快速收敛算法 _正弦信号序列频率的估计 _ 关于正弦波序列频率的估计方法,有很多文献做了讨论 1 2,1 a3 1,根据精度b速度b实时性和信号质量的不同,可以采取不同的方法。当式(1)所述正弦信号被执行均匀采样后获得的等间隔采样序列为y1,y2,yn,其表述的函数关系如式(2)所示。令c(i)=C0c o s(i+0),i=1,2,n(d)若测量点yi的误差为ei,f=2 c o s ,则有yi=y(i)+ei=c(i)+D+ei(g)c(i)+c(i-2)=(2 c o s )h c(i-1)=fh c(i-1)(i)y(i)-D+y(i-2)-D=f(y(i-1)-D)()yi-D-ei+yi-2-D-ei-2=f(yi-1-D-ei-1)(j)yi+1-D-ei+1+yi-1-D-ei-1=f(yi-D-ei)(k)令li=yi-yi-1(1 0)mi=ei-ei-1(1 1)则由式(j)和(k)得li+1+li-1-f li=mi+1+mi-1-f mi(1 2)由于可以认为ei为随机误差,故式(1 2)右侧也可以认为是随机误差。通过选取f使得n=n-1i=3(li+1+li-1-f li)2=oi n(1 3)令p np f=0,则得f=lnln-1+l3l2+2n-1i=dlili-1n-1i=3l2i(1 d)=S T c c o s(f/2)(1 g)即为数字角频率的最小二乘估计值。四参数正弦波曲线拟合的快速算法将第2 1节和第1节的运算过程结合起来,即构成了组合形式的四参数正弦波最小二乘拟合算法。从上述过程可见,两种算法均只涉及使用了加b减b乘b除的四则运算,属于不需要经过叠代的闭合算法,因而没有收敛性问题,即过程总是收敛的。q仿真实验验证对于四参数正弦波曲线拟合算法的较全面评价需HgH计 测 技 术理论与实践要使用蒙特卡罗仿真方法,或使用文献 3 2 所述方法首先获得其指标以最终确定优劣。限于篇幅,这里只选择在一组特定条件下的不同算法的结果比较来间接考察其相对优劣。设 定,仿 真 信 号 标 称 幅 度 为4V、标 称 频 率 为6 2 5 4 3 2 1Hz、标称相位0r a d、标称直流分量0V。仿真测量系统量程-5 5V、采样速率4G S a/s。当仿真A/D位数变化时,其用本文上述方法获得的四参数拟合结果如表1所示。表2为相同条件下使用文献 4 所述的频率搜索法获得的拟合结果,表3为相同条件下使用Ne wt o n迭代法获得的拟合结果 6。表1本文所述方法四参数正弦波拟合结果(A/D)/b i t s幅度/V频率/Hz相位/r a d直流分量/V有效位数/b i t s信噪比/d B噪信比1 60.4 4 5 1 0 7 6 36 5 5 5 2 9 4-0.3 9 36.6 9 1 0-30.0 3 9 82.0 48.9 21 73.4 4 8 1 4 6 0 26 3 3 2 7 2 3-0.9 1 42.8 6 1 0-31.0 37.9 90.5 8 01 83.9 5 5 0 9 7 1 76 2 7 4 8 0 6-0.2 3 22.5 0 1 0-52.8 81 9.1 30.1 4 01 93.9 9 6 3 5 7 7 76 2 5 9 1 7 3-0.0 4 7 4-5.1 1 1 0-54.9 53 1.5 90.0 3 3 12 03.9 9 9 6 5 8 2 76 2 5 5 2 2 0-0.0 0 0 7 8 7-1.4 3 1 0-57.3 84 6.2 20.0 0 6 1 32 13.9 9 9 9 2 6 5 76 2 5 4 5 5 30.0 0 7 0 9-4.6 3 1 0-69.3 45 8.0 10.0 0 1 5 82 23.9 9 9 9 8 7 0 16 2 5 4 3 6 50.0 0 9 3 1-1.3 9 1 0-61 1.7 67 2.5 70.0 0 0 2 9 52 34.0 0 0 0 0 4 5 66 2 5 4 3 0 60.0 1 0 0-1.9 1 1 0-71 3.2 58 1.5 50.0 0 0 1 0 52 43.9 9 9 9 9 1 7 66 2 5 4 3 4 90.0 0 9 5 0-4.8 9 1 0-71 2.4 17 6.4 90.0 0 0 1 8 8表2频率搜索迭代方法四参数正弦波拟合结果(A/D)/b i t s幅度/V频率/Hz相位/r a d直流分量/V有效位数/b i t s信噪比/d B噪信比33.8 1 5 3 2 2 6 36 2 5 4 3 8 29.1 7 1 0-3-0.5 0 03.3 42 1.9 00.1 0 643.9 3 4 5 6 1 3 86 2 5 4 3 4 99.5 2 1 0-3-0.2 5 04.3 02 7.7 15.2 5 1 0-253.9 7 6 8 3 4 2 16 2 5 4 3 3 09.7 1 1 0-3-0.1 2 55.2 93 3.6 72.6 2 1 0-263.9 9 1 7 9 2 1 06 2 5 4 3 2 59.7 6 1 0-3-6.2 5 1 0-26.2 93 9.6 91.3 0 1 0-273.9 9 7 0 9 2 6 06 2 5 4 3 2 19.8 1 1 0-3-3.1 3 1 0-27.3 04 5.7 46.4 9 1 0-383.9 9 8 9 6 2 0 66 2 5 4 3 2 39.8 1 1 0-3-1.5 7 1 0-28.3 05 1.7 93.2 3 1 0-393.9 9 9 6 2 3 8 76 2 5 4 3 2 29.8 1 1 0-3-7.8 3 1 0-39.3 15 7.8 31.6 1 1 0-31 03.9 9 9 8 7 0 4 96 2 5 4 3 2 29.8 1 1 0-3-3.9 1 1 0-31 0.3 16 3.8 68.0 4 1 0-41 13.9 9 9 9 5 4 8 46 2 5 4 3 2 19.8 2 1 0-3-1.9 5 1 0-31 1.3 16 9.8 94.0 2 1 0-41 23.9 9 9 9 8 7 3 16 2 5 4 3 2 19.8 3 1 0-3-9.7 6 1 0-41 2.3 27 5.9 42.0 0 1 0-41 33.9 9 9 9 9 7 1 46 2 5 4 3 2 19.8 3 1 0-3-4.9 0 1 0-41 3.3 18 1.9 31.0 0 1 0-41 44.0 0 0 0 0 1 9 86 2 5 4 3 1 59.9 0 1 0-3-2.4 5 1 0-42 6.5 11 6 1.3 71.0 7 1 0-81 54.0 0 0 0 0 0 3 86 2 5 4 3 2 19.8 2 1 0-3-1.2 2 1 0-41 5.8 59 7.2 41.7 2 1 0-51 64.0 0 0 0 0 0 2 06 2 5 4 3 2 19.8 2 1 0-3-6.0 9 1 0-51 6.2 29 9.4 31.3 4 1 0-51 73.9 9 9 9 9 9 9 76 2 5 4 3 2 19.8 2 1 0-3-3.0 5 1 0-51 5.8 09 6.9 01.7 9 1 0-51 84.0 0 0 0 0 0 0 66 2 5 4 3 2 19.8 2 1 0-3-1.5 2 1 0-51 5.5 29 5.2 22.1 7 1 0-51 94.0 0 0 0 0 0 0 16 2 5 4 3 2 19.8 2 1 0-3-7.5 2 1 0-61 6.5 21 0 1.2 51.0 9 1 0-52 04.0 0 0 0 0 0 0 26 2 5 4 3 2 19.8 2 1 0-3-3.7 3 1 0-61 7.8 01 0 8.9 34.4 9 1 0-62 13.9 9 9 9 9 9 4 46 2 5 4 3 2 39.8 0 1 0-3-1.8 3 1 0-62 6.0 41 5 8.5 51.4 8 1 0-82 23.9 9 9 9 9 9 1 66 2 5 4 3 2 49.7 9 1 0-3-8.9 3 1 0-72 6.9 71 6 4.1 97.7 4 1 0-92 33.9 9 9 9 9 9 4 16 2 5 4 3 2 39.8 0 1 0-3-4.0 7 1 0-72 5.8 31 5 7.2 91.7 1 1 0-82 43.9 9 9 9 9 9 4 26 2 5 4 3 2 39.8 0 1 0-3-1.7 3 1 0-72 6.5 51 6 1.6 41.0 4 1 0-862 0 0 6年第2 6卷第1期理论与实践表3 Ne wt o n迭代法四参数正弦波拟合结果(A/D)/b i t s幅度/V频率/Hz相位/r a d直流分量/V有效位数/b i t s信噪比/d B噪信比33.8 1 3 5 0 9 5 6 76 2 5 4 3 8 4-2.9 0 1 0-4-0.5 0 03.3 42 1.9 00.1 0 643.9 3 4 2 1 0 3 4 56 2 5 4 3 4 9-3.0 3 1 0-4-0.2 5 04.3 02 7.7 15.2 5 1 0-253.9 7 6 7 3 5 0 5 56 2 5 4 3 3 0-1.1 1 1 0-4-0.1 2 55.2 93 3.6 72.6 2 1 0-263.9 9 1 7 7 8 9 0 36 2 5 4 3 2 5-6.9 6 1 0-5-6.2 5 1 0-26.2 93 9.6 90.0 1 3 073.9 9 7 0 8 8 0 7 86 2 5 4 3 2 1-9.4 3 1 0-6-3.1 3 1 0-27.3 04 5.7 46.4 9 1 0-383.9 9 8 9 6 2 0 7 56 2 5 4 3 2 3-1.9 2 1 0-5-1.5 7 1 0-28.3 05 1.7 93.2 3 1 0-393.9 9 9 6 2 3 4 9 46 2 5 4 3 2 2-1.4 4 1 0-5-7.8 3 1 0-39.3 15 7.8 31.6 1 1 0-31 03.9 9 9 8 6 9 6 4 56 2 5 4 3 2 2-1.6 1 1 0-5-3.9 1 1 0-31 0.3 16 3.8 78.0 3 1 0-41 13.9 9 9 9 5 3 9 8 46 2 5 4 3 2 1-2.9 5 1 0-6-1.9 5 1 0-31 1.3 16 9.9 14.0 1 1 0-41 23.9 9 9 9 8 6 4 3 96 2 5 4 3 2 15.8 1 1 0-6-9.7 6 1 0-41 2.3 27 5.9 52.0 0 1 0-41 33.9 9 9 9 9 6 2 7 56 2 5 4 3 2 12.4 4 1 0-6-4.9 0 1 0-41 3.3 28 1.9 71.0 0 1 0-41 43.9 9 9 9 9 9 3 0 66 2 5 4 3 2 14.7 1 1 0-7-2.4 5 1 0-41 4.3 28 7.9 95.0 0 1 0-51 54.0 0 0 0 0 0 3 6 86 2 5 4 3 2 1-1.5 1 1 0-7-1.2 2 1 0-41 5.3 29 4.0 32.4 9 1 0-51 64.0 0 0 0 0 0 1 9 66 2 5 4 3 2 11.3 5 1 0-7-6.0 9 1 0-51 6.3 21 0 0.0 71.2 4 1 0-51 73.9 9 9 9 9 9 9 7 16 2 5 4 3 2 12.6 5 1 0-8-3.0 5 1 0-51 7.3 31 0 6.1 06.2 1 1 0-61 84.0 0 0 0 0 0 0 5 26 2 5 4 3 2 1-1.7 5 1 0-8-1.5 2 1 0-51 8.3 21 1 2.0 83.1 2 1 0-61 94.0 0 0 0 0 0 0 0 76 2 5 4 3 2 12.8 3 1 0-8-7.5 2 1 0-61 9.3 11 1 8.0 61.5 7 1 0-62 04.0 0 0 0 0 0 0 2 16 2 5 4 3 2 14.3 9 1 0-9-3.7 3 1 0-62 0.3 01 2 4.0 47.8 8 1 0-72 14.0 0 0 0 0 0 0 1 56 2 5 4 3 2 18.9 3 1 0-1 0-1.8 1 1 0-62 1.2 91 2 9.9 43.9 9 1 0-72 23.9 9 9 9 9 9 9 9 66 2 5 4 3 2 1-3.1 4 1 0-9-8.5 8 1 0-72 2.1 51 3 5.1 62.1 9 1 0-72 34.0 0 0 0 0 0 0 0 26 2 5 4 3 2 1-1.7 6 1 0-9-3.8 3 1 0-72 2.8 11 3 9.1 11.3 9 1 0-72 43.9 9 9 9 9 9 9 9 46 2 5 4 3 2 1-1.4 2 1 0-1 0-1.4 9 1 0-72 3.1 51 4 1.1 61.1 0 1 0-7从上述仿真过程及结果可见,本文所述方法的运算速度最快,也最为便捷,没有收敛性问题,但对数据要求较高,且拟合精度最差,当A/D位数低到1 6位时,拟合结果已经非常差了,原因在于本文方法的频率估计精度受数据质量影响非常严重;频率搜索方法的运算 速度最慢,但收敛性较好,精度中等;而Ne wt o n法的精度最高,速度中等,但有时拟合结果不收敛。4结束语综上所述可见,用本文所述方法进行四参数正弦波拟合,具有如下特点:1)不存在过程的收敛问题,因而总是可以获得一个估计结果;2)不需要对参数进行预估计,算法简洁;3)不需要迭代过程,所有结果都可以一次运算获得,因而具有较高的运算速度;4)量化噪声和随机误差对本文所述频率估计算法影响非常大,因此其最适合于测量数据比较精确的序列模型估计;5)估计过程对序列所含信号周期数没有要求,故可以依据精确的局部周期数据对参数模型进行估计,此为其重要特点;6)如果对测量序列进行一下预处理,使其近似含有整数个信号周期,将可以进一步降低估计结果的不确定度。7)如果选取其它效果更加优良的信号频率估计方法,将可以获得更加有价值的正弦波拟合算法。参 考 文 献 1 梁志国,孙王景宇.正弦波模型化测量方法及应用 J .航空计测技术,2 0 0 1,2 1(6):3-7.2 梁志国,朱济杰.调相信号的一种数字化精确解调方法 J .计量学报,2 0 0 3,2 4(4):3 2 5-3 2 9.3 梁志国,孙王景宇.调幅信号的数字化解调 J .测控技术,2 0 0 4,2 3(1 0):1 7-2 0.4 梁志国,孙王景宇.调频信号的数字化解调 J .测试技术学报,2 0 0 5,1 9(2):1 9 0-1 9 4.5 梁志国,朱济杰,孙王景宇.正弦信号总失真度的一种精确评价方法 J .计量学报,2 0 0 3,2 4(2):1 4 4-1 4 8.6 J o h n Ku f f e l,T e r e n c e R Mc C o mb,R y s z a r d Ma l e ws k i.C o mp a r a t i v eE v a l u a t i o no f C o mp u t e rMe t h o d sf o rC a l c u l a t i n gt h eB e s t-F i t S i n u s o i dt ot h eD i g i t a l R e c o r do f aHi g h-P u r i t yS i n eWa v e J .I E E ET r a n s a c t i o n so nI n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 8 7,3 6(3):4 1 8-4 2 2.7 T e r e n c eR Mc C o mb,J o h nKu f f e la n dB r y a n tCL eR o u x,AC o mp a r a t i v e E v a l u a t i o n o fS o me P r a c t i c a lAl g o r i t h ms Us e d i n t h eE f f e c t i v eB i t sT e s to fWa v e f o r m R e c o r d e r s J .I E E ET r a n s a c t i o n so nI n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 8 9,3 8(1):3 7-4 2.8 J e n gYC.Hi g hp r e c i s i o ns i n u s o i d a l f r e q u e n c ye s t i ma t o rb a s e do n we i g h t e d l e a s t s q u a r e me t h o d J.I E E E T r a n s a c t i o n s o nI n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 8 7,3 6(2):1 2 4-1 2 7.(下转第4 7页)7计 测 技 术理论与实践小了定焦带来的测量不确定度。另外采用双刀口扫描可进行器件调平判断和环境振动、噪声的判断,以便及时发现问题,设法消除其对测量不确定度的影响。参 考 文 献 1 麦伟麟.光学传递函数及数理基础 M.北京:国防工业出版社,1 9 7 9.2 王庆有.C C D应用技术 M.天津:天津大学出版社,1 9 9 1.3 Mi k eMa r c h y wk a,D e n n i sG S o c k e r.Mo d u l a t i o n t r a n s f e rf u n c t i o nme a s u r e me n tt e c h n i q u ef o rs ma l l-p i x e ld e t e c t o r s J .Op t i c a lE n g i n e e r i n g,1 9 9 3,3 2(2):9 5-1 0 0.4 Mi l e s N We r n i c k,G Mi c h a e lMo r r i s.E f f e c to fs p a t i a lc o h e r e n c eo nk n i f e-e d g eme a s u r e me n t so fd e t e c t o rmo d u l a t i o nt r a n s f e rf u n c t i o n J .Op l i c a l E n g i n e e r i n g,1 9 9 1,3 0(9):1 3 9 4-1 3 9 8.5 G B/T 4 3 1 5.2-1 9 9 8,光学传递函数测量导则 S .6 J o h n C F e l t z,Mo h a mma d A Ka r i m.Mo d u l a t i o n t r a n s f e rf u n c t i o no fc h a r g e-c o u p l e dd e v i c e s J .Ap p l i c a t i o nOp t i c a l,1 9 9 0,2 7 (2):7 1 7-7 2 2.(上接第7页)9 _i h-C h y u n J e n q,h i l i p B C r o s a y.S i n e wa v e a r a me t e rE s t i ma t i o nAl g o r i t h m wi t hAp p l i c a t i o nt oWa v e f o r m D i g i t i z e rE f f e c t i v eB i t s Me a s u r e me n t J .b E E E T r a n s a c t i o n so n b n s t r u me n t a t i o n a n dMe a s u r e me n t,1 9 8 8,3 7(4):5 2 9-5 3 2.1 0 F e l i c e C e n n a mo,a s q u a l e D a p o n t e,Ma r i o S a v a s t a n o.D y n a mi c T e s t i n g a n d D i a g n o s t i c s o fS i g n a lAn a l y z e r s J .b E E ET r a n s a c t i o n so nb n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 9 2,4 1(6):8 4 0-8 4 4.1 1 G i a q u n i n t oN,T r o t t a A.F a s t a n da c c u r a t e AD Ct e s t i n gv i a a ne n h a n c e d s i n e wa v e f i t t i n g a l g o r i t h m J .b E E E t r a n s a c t i o n s o nb n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 9 7,4 6(4):1 0 2 0-1 0 2 5.1 2 c h a n gJd,c e i n mi n,f e i a o,S J i n we i.S i n ewa v ef i ta l g o r i t h m a a s e do nt o t a l l e a s t-s q u a r e sme t h o dwi t ha p p l i c a t i o nt oAD Ce f f e c t i v ea i t sme a s u r e me n t J .b E E E t r a n s a c t i o n so nb n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 9 7,4 6(4):1 0 2 6-1 0 3 0.1 3 e t e rfa n d e l.r o p e r t i e s o ft h e b E E E-S T D-1 0 5 7 F o u r-a r a me t e r S i n e Wa v e F i t Al g o r i t h m J .b E E E T r a n s a c t i o n s o nb n s t r u me n t a t i o na n dMe a s u r e me n t,2 0 0 0,4 9(6):1 1 8 9-1 1 9 3.1 4 田社平,王坚,颜德田,等.基于遗传算法的正弦波信号参数提取方法 J .计量技术,2 0 0 5,(5):3-5.1 5 b E E ES t d1 0 5 7-1 9 9 4,b E E Es t a n d a r df o rd i g i t i z i n gwa v e f o r mr e c o r d e r s S .1 6 B e r t o c c oM,Na r d u z z iC.S i n e-f i tv e r s u sd i s c r e t eF o u r i e rt r a n s f o r m-a a s e da l g o r i t h m i nS Ng t e s t i n go fwa v e f o r m d i g i t i z e r s J .b E E Et r a n s a c t i o n so nb n s t r u me n t a t i o na n dMe a s u r e me n t,1 9 9 7,4 6(2):4 4 5-4 4 8.1 7 h a d k eA G,T h o r pJS,Ad a mi a kM.A n e w me a s u r e me n tt e c h n i q u ef o rt r a c k i n gv o l t a g ep h a s o r s,l o c a l s y s t e m f r e q u e n c ya n dr a t eo f c h a n g e o f f r e q u e n c y J .b E E Et r a n s a c t i o n s o n o we r Ap p a r a t.S y s t.,1 9 8 3,1 0 2:1 0 2 5-1 0 3 8.1 8 G i r g i sA A,fwa n gT h D.Op t i ma lE s t i ma t i o no fv o l t a g ep h a s o r sa n d f r e q u e n c y d e v i a t i o n u s i n g l i n e a ra n d n o n l i n e a rKa l ma nf i l t e r i n g J .b E E Et r a n s a c t i o n so n o we rAp p a r a t.S y s t.,1 9 8 4,1 0 3:2 9 4 3-2 9 4 9.1 9 S a c h d e vM S,G i r a yM M.A l e a s ts q u a r et e c h n i q u ef o rd e t e r mi n i n gp o we rs y s t e m f r e q u e n c y J .b E E E t r a n s a c t i o n so n o we rAp p a r a t.S y s t.,1 9 8 5,1 0 4:4 3 7-4 4 3.2 0 S a c h d e v MS,G i r a y MM.Of f-n o mi n a l f r e q u e n c yme a s u r e me n t i ne l e c t r i cp o we r s y s t e ms J .b E E Et r a n s a c t i o n so n o we rD e l i v e r y.1 9 8 9,4:1 5 7 3-1 5 7 8.2 1 G i r g i sA A,e t e r s o nW h.Ad a p t i v ee s t i ma t i o no fp o we rf r e q u e n c yd e v i a t i o na n di t sr a t eo fc h a n g ef o rc a l c u l a t i n gs u d d e np o we rs y s t e m o v e r l o a d s J .b E E E t r a n s a c t i o n so n o we rD e l i v e r y,1 9 9 0,5:5 8 5-5 9 4.2 2 T e r z i i ajj,D i r i cM B,Ko v a c e v i cBD.jo l t a g ep h a s o ra n dl o c a ls y s t e m f r e q u e n c ye s t i ma t i o nu s i n gNe wt o n-t y p ea l g o r i t h ms J .b E E Et r a n s a c t i o n so n o we r D e l i v e r y,1 9 9 4,9:1 3 6