傅里叶变换对评价正弦信号源波形失真的应用

傅里叶变换对评价正弦信号源波形失真的应用摘要：介绍了一种直接利用经典傅里叶变换技术评价正弦信号源波形总失真度的过程和方法。借助于周期的精确测量技术,实现基波参数的精确测量。用时域能量方式计算失真,其内涵全面,包括了谐波、杂波、噪声等全部影响,并对测量系统本身的影响进行了补偿同时,对不同条件下的失真度测量误差进行了分析,用一组仿真数据的实验结果,验证了正确性及可行性。所述方法可用于正弦信号波形失真度的精确测量和计量校准,尤其适用于大失真度的测量。关键词:计量学;失真度;傅里叶变换;正弦波;评价1.引言对于正弦信号源产生的波形来说,以H1表示其基波有效值幅度,Hi表示其含有的谐波有效值幅度(i=2,...,M),nH表示其含有的噪声及非谐波分量(不包括直流分量)的有效值幅度。此时,该信号波形含有失真和噪声(其实噪声也是一种失真)。则该信号波形的总失真度TD定义为）（112i2i2nHHHTDM在超低频、高频及射频范围内,FFT法是较普遍采用的正弦波形失真度的数字化测量方法。通过使用FFT获得的基波和各次谐波分量,代入式(1)中,获得正弦波形的总谐波失真度。该方法的优点是简单快捷,尤其适合大失真下的测量;弱点和局限性也非常明显:要求测量序列的长度为2的整数次幂;要求序列含有整数个信号周期,即满足同步采样条件;当同步采样条件不满足时,FFT法本身的”栅栏效应”、“泄漏”等将给测量结果带来误差;尽管有准同步采样和补偿性方在使用,其高次谐波的测量准确度仍无法完全令人满意;④使用FFT法只能对有限次谐波(如19次以下,57次以下等)进行计算,无法包括所有谐波分量和非谐波分量,也未包含噪声分量;⑤测量设备本身的量化误差对测量结果的影响未予以考虑。所有这些因素,均限制了失真度的测量准确度。本文所述内容,将主要讨论一种基于傅里叶变换的失真度测量方法,其特点是可实现精确测量,并能避免上述弱点。2.周期信号的傅里叶分解一个周期函数y(t)=y(t+T),只要满足狄利克雷条件(在一个周期中有有限个极值,并且或者处处连续,或者有有限个第1类间断点),则可用傅里叶级数来表示该函数[1]。)5(22)4(arg)3()2)(cos(sincos)(y02210110tNTABtgBACmCAmBmAAtmmmmmmmmtmmtmmtm其中,Δt为信号采样间隔;0N为每个信号周期内含有的采样点数。dtsinmty27dtcosmty26dtty1t22t22220TTmTTmTTTBTATA）（）（）（）（）（求取mm0BAA、、(m=1,2,...)并用它们来描述y(t)的过程,就是周期信号y(t)的傅里叶分解过程,也是其谐波分析过程或频谱分析过程。通过信号y(t)的离散采样值,(i=0,1,...,N-1),准确获得参数mm0BAA、、,一直是人们所关心的问题，业已证明[2],只要满足对第m次谐波在一个信号周期内能采样两个以上信号点,且N个采样点中恰好采集了整数个信号周期,则:)11(2sin)(2)10(2cos)(2)9()(110i10i10i0NmiityNBNmiityNAityNANmNmN这个过程及其条件和方法被称为同步采样。在多数情况下,同步采样也特指信号周期是采样间隔的整数倍的采样状态。实际上,N个采样点恰好含有整数个信号周期的情况是很难严格实现的,因而式(9)~(11)式直接用于谐波测量分析会产生测量误差。由此引出了一系列以迭代收敛运算为特征的准同步采样方法[3~5]和已知信号周期为条件的补偿性方法[6,7]均被用于谐波分析。这些方法共同的特点是运算比较复杂。本文将试图从傅里叶分解的最基本关系出发,寻找出一种简单、精确的周期信号谐波测量分析方法,并用于正弦波形的失真度测量。3测量原理及方法设正弦信号源产生的波形函数是周期函数y(t)。由式(5)可见,当N个采样点中不是恰好含有整数个信号周期时,同步采样条件不被满足。式(10)和式(11)实际上是在对周期为NTt的信号进行谐波分析,而不是对周期为T的信号进行谐波分析!以前者代替后者,将产生测量误差;另一方面,它也不满足式(6)~(8)式的条件约定,直接使用式(10)和式(11)仍将带来测量误差。鉴于上述原因,人们实际上可以在进行谐波分析之前,首先对被分析信号y(t)的周期T进行精确测量,并获得精确的0N(通常不是整数)。然后,通过判定过均值点截取首尾皆处于均值点附近、含有整数个信号周期的N个信号采集点）（ity,(i=0,1...,N-1),并进行谐波分析,获得信号y(t)的各次谐波参数。其过程如下:(1)以采样间隔t对周期信号y(t)进行波形采样,获得采样序列）（ity,(i=0,1...,N-1);(2)在上述样本中,截取含有整数个信号周期的N个样本(误差极限为1)点）（ity,(i=0,1...,N-1);用式(9)计算波形均值:10i0)(1NityNA;同时获得零均值序列00)(y)(yAitit,(i=0,1...,n-1);(3)使用周期精确测量方法[8]获得信号周期T,计算t0TN;(4)通过判定过零点,从序列中截取首尾值皆处于零附近,并含有整数个信号周期的信号样本点）（ity,(i=0,1...,N-1),以0N代替N,按式(12)~式(13)计算获得周期信号y(t)的谐波分量参数mmBA、,(m=1,2...)(5)进而通过式(4)、式(5)获得相应的幅度mC和相位m;(6)通过下式计算获得残差序列:)1...,1,0(;tmax2sin2cos)t(tmax01010NiiNiBNiAiyi）（）（残差有效值为)14()(11i2NitN波形数据的总失真度TDZ为)15(21CTDZ若信号峰峰值与波形测量系统量程之比为G,测量系统有效位数为BD,则测量系统量化误差带来的失真度为:)16(43222212BDCTD4误差分析由于使用了周期精确测量技术,由周期带来的谐波分析误差可以忽略不计;由式(16)的全微分,得失真度TD的误差为)17)(432l(22131221BDnCCCTDTDBD其中，有效位数BD作为修正因子，其误差BD约为0.1Bits，可参见有关文献从校准中获得[10]。设信号y(t)在过均值点附近连续并导数)(yt存在且接近于常值,约为)0(y或))1((ytN可以估算为)]1()0([5.02)1(ˆ)1(ˆ2)1()1())1((y2)1(ˆ)1(ˆ2)1()1()0(yNyyytNyNytNyNytNtyytyy设0)(maxyAtyp,当不采取从均值点处开始截取整数个信号周期时,执行式(9)进行均值0A的估算,给0A带来的最大误差约为)18(y00NAAp当采取从均值点处开始截取整数个信号周期时,执行式(9)进行均值0A的估算,给0A带来的最大误差约为)19(y000NANtAA在上述方法中,由于采样序列N个点中不恰好含有整数个信号周期(周期长度为0N),在式(12)式(13)中给累加和所造成的最大误差约为ty,给mA和mB造成的最大误差约为)21(22sin2)20(22cos20m0mtyNNmtyNBtyNNmtyNA给造成的最大误差为)22(22max2max2maxNN由于mA和mB具有正交性，mA和mB可以认为不相关，该因素给mC带来的最大误差mC约为)29(222mmtyNBACm式（4）可以导出相应的相位误差：)24)(cossin(1)coscossin(12mmmmmmmmmmmmBACBAA由式（17）可得：)25)(432lny22(22312212maxBDBDNCtNCTDTD一般说来,对于具体的周期信号y(t）,y是确定不变的,从式(25)可见,增大序列样本数N可以降低失真度分析的误差TD,降低采样间隔t(增高采样速率)也可以降低失真度分析的误差,采取从过均值点截取序列周期并经过剔除直流分量后进行分析,同样将减小失真度分析误差。当信号y(t）在过均值点附近不连续且出现巨大跳变时，如方波和锯齿波时，设跳变增量为py（方波或锯齿波峰值），则上述方法中由于采样序列N个点中不恰好含有整个信号周期（长度为0N），在式（12）、式（13）中给累加和造成的最大误差约为py，给mA和mB造成的最大误差约为)28(2)27(22sin2)26(22cos20m0mpmppppyNCyNNmyNByNNmyNA此时,从式(25)可见,增大序列样本数N可以降低失真度分析误差,降低采样间隔t对降低失真度分析误差没有明显作用。通常,每个信号周期至少有2个“过均值点”,多数时候这两个点处的导数有明显差异,由式(18)~式(25)可见,选取导数存在且导数绝对值最小的“过均值点”进行样本截取,然后执行谐波分析,将有利于获得良好的失真度测量结果。5.仿真实验验证用上述公式和结论在计算机上进行仿真,获得结果如表1所示。通常,未考虑对量化误差等带来的影响进行修正时,均以测量数据失真度TDZ代替输入信号失真度TD,实际上的输入信号失真度应该是TD。在已知测量系统动态有效位数BD的情况下,可对测量数据失真度TDZ进行修正,获得输入信号失真度的较准确评价值TD。表1验证了这一结论。从表1可见,APD位数对于失真度的评价具有较大影响,当其过低时,即使使用修正结果,也无法完全消除其影响,这表明实际应用中,在测量较小失真度和进行较精确的失真度测量时,需要选择动态有效位数足够大的测量系统。6结束语综上所述可见,本文所述失真度测量方法具有如下特点:(1)直接从经典傅里叶变换定义出发执行运算,不使用FFT方法,因而不要求测量序列长度为2的整数次幂;(2)不要求测量序列恰好含有整数个信号周期;即不要求满足同步采样条件。(3)只对基波分量参数进行分解变换,不涉及高次谐波,有利于提高测量准确度;(4)直接在时域使用信号与基波的残差来计算失真,属于软件陷波方法,其内涵包括了所有谐波、杂波及噪声失真,不存在方法原理上的缺陷;(5)通过已知的测量系统动态有效位数,可以对测量仪器本底噪声带来的影响进行修正,从而进一步提高测量准确度;(6)测量误差因素可以明确估计与控制。目前其它FFT分析方法尚未提供如此具体的结论。从上述仿真实验验证可见,本文所述方法用于正弦波形的失真度测量,具有足够高的测量准确度,尤其对于大失真度情况。这主要是由于上述方法采取了如下关键性措施减小谐波分析误差:(1)使用了信号周期的精确测量技术,以0N代替N进行分解运算,减小了谐波频率定值误差以及所带来的运算误差;(2)在函数序列的均值点附近开始截取近整数个信号周期的数据N,执行谐波分析,减少了非同步采样带来的运算误差和方法误差;(3)从序列中减去均值分量0A后进行谐波参数运算,减少了0A对各谐波参数带来的影响。采用这些措施后,并未出现对测量过程的特别要求,不要求同步采样状态,实际上是处于非整周期采样状态,使得上述方法具有较广泛的适应性和简便易行性,因此,可望在计量测试和精确测量中有良好的应用前景。参考文献[1]潘士先.电路分析(下册)[M].北京:国防工业出版社,1979,54-55.[2]何天祥.周期信号各参数的离散采样测量法[J].航空计测技术,1997,17(3):3-4.[3]戴先中.准同步采样及其在非正弦功率测量中的应用[J].仪器仪表学报,1984,5(4):390-396.[4