傅里叶变换对评价正弦信号源波形失真的应用
摘要:介绍了一种直接利用经典傅里叶变换技术评价正弦信号源波形总失真度的 过程和方法。借助于周期的精确测量技术,实现基波参数的精确测量。用时域能量方式计算失真,其内涵全 面,包括了谐波、杂波、噪声等全 部影响,并对测量系统本身的影响进行了补偿同时,对不同条件下的失真度测量误差进行了分析, 用一组仿真数据的实验结果, 验证了正确性及可行性。所述方法可用于正弦信号波形失真度的精确 测量和计量校准,尤其适用于大失真度的测量。
关键词: 计量学; 失真度; 傅里叶变换;正弦波;评价
1.引言
对于正弦信号源产生的波形来说,以H1 表示其基波有效值幅度, H i 表示其含有的谐波有效值幅度( i = 2,..., M ) , { EMBED Equation.KSEE3 \* MERGEFORMAT |Hn表示其含有的噪声及非谐波
分量(不包括直流分量) 的有效值幅度。此时,该信号波形含有失真和噪声(其实噪声也是一种失真)。则该信号波形的总失真度 TD 定义为
在超低频、高频及射频范围内, FFT 法是较普遍采用的正弦波形失真度的数字化测量方法。通过使用 FFT 获得的基波和各次谐波分量, 代入式(1)中,获得正弦波形的总谐波失真度。该方法的优点是简单快捷,尤其适合大失真下的测量;弱点和局限性也非常明显:①要求测量序列的长度为2的整数次幂;②要求序列含有整数个信号周期, 即满足同步采样条件;③当同步采样条件不满足时, FFT法本身的”
栅栏效应”、“泄漏”等将给测量结果带来误差;尽管有准同步采样和补偿性方在使用,其高次谐波的测量准确度仍无法完全令人满意; ④使用 FFT 法只能对有限次谐波( 如19次以下,57次以下等) 进行计算,无法包括所有谐波分量和非谐波分量,也未包含噪声分量;⑤测量设备本身的量化误差对测量结果的影响未予以考虑。所有这些因素, 均限制了失真度的测量准确度。本文所述内容, 将主要讨论一种基于傅里叶变换的失真度测量方法, 其特点是可实现精确测量, 并能避免上述弱点。
2.周期信号的傅里叶分解
一个周期函数 y(t) = y(t+T) , 只要满足狄利克雷条件( 在一个周期中有有限个极值,并且或者处处连续, 或者有有限个第 1 类间断点) , 则可用傅里叶级数来表示该函数[1]。
其中,Δt为信号采样间隔; 为每个信号周期内含有的采样点数。
求取( m = 1, 2,...)并用它们来描述 y(t)的过程,就是周期信号y(t)的傅里叶分解过程,也是其谐波分析过程或频谱分析过程。
通过信号y(t)的离散采样值, ( i = 0,1,..., N - 1) , 准确获得参数 , 一直是人们所关心的问题,业已证明[2] ,只要满足对第 m 次谐波在一个信号周期内能采样两个以上信号点, 且N个采样点中恰好采集了整数个信号周期, 则:
这个过程及其条件和方法被称为同步采样。在多数情况下,同步采样也特指信号周期是采样间隔的整数倍的采样状态。实际上, N个采样点恰好含有整数个信号周期的情况是很难严格实现的,因而式
(9)~(11)式直接用于谐波测量分析会产生测量误差。由此引出了一系列以迭代收敛运算为特征的准同步采样方法[3~5]和已知信号周期为条件的补偿性方法[6,7]均被用于谐波分析。这些方法共同的特点是运算比较复杂。本文将试图从傅里叶分解的最基本关系出发,寻找出一种简单、精确的周期信号谐波测量分析方法, 并用于正弦波形的失真度测量。
3测量原理及方法
设正弦信号源产生的波形 函数是周期函数y(t)。由式(5)可见,当N个采样点中不是恰好含有整数个信号周期时,同步采样条件不被满足。式(10) 和式( 11) 实际上是在对周期为的信号进行谐波分析, 而不是对周期为T的信号进行谐波分析!以前者代替后者,将产生测量误差;另一方面,它也不满足式( 6)~(8)式的条件约定,直接使用式
(10)和式(11)仍将带来测量误差。鉴于上述原因,人们实际上可以在进行谐波分析之前, 首先对被分析信号 y(t)的周期T进行精确测量, 并获得精确的(通常不是整数) 。然后,通过判定过均值点截取首尾皆处于均值点附近、含有整数个信号周期的 N 个信号采集
点,(i=0,1...,N-1), 并进行谐波分析,获得信号 y(t)的各次谐波参数。其过程如下:
(1) 以采样间隔对周期信号 y(t)进行波形采样, 获得采样序
列 ,(i=0,1...,N-1);
(2) 在上述样本中, 截取含有整数个信号周期的N个样本(误差极限为) 点 ,(i=0,1...,N-1); 用式(9)计算波形均值:
;
同时获得零均值序列
,(i=0,1...,n-1);
(3) 使用周期精确测量方法[8]获得信号周期T,计算;
(4) 通过判定过零点,从序列中截取首尾值皆处于零附近, 并含有整数个信号周期的信号样本点,(i=0,1...,N-1), 以代替 N , 按式(12) ~ 式(13) 计算获得周期信号y(t)的谐波分量参数 ,(m=1,2...) ( 5) 进而通过式(4) 、式(5) 获得相应的幅度和相位 ;
( 6) 通过下式计算获得残差序列:
残差有效值为
波形数据的总失真度 TDZ 为
若信号峰峰值与波形测量系统量程之比为 G,测量系统有效位数为 BD , 则测量系统量化 误差带来的失真度 为:
4 误差分析
由于使用了周期精确测量技术, 由周期带来的谐波分析误差可
以忽略不计; 由式( 16) 的全微分, 得失真度 TD 的误差为
其中,有效位数BD作为修正因子,其误差约为0.1Bits,可参见有关文献从校准中获得[10]。
设信号 y(t) 在过均值点附近连续并导数存在且接近于常值, 约为或可以估算为
设,当不采取从均值点处开始截取整数个信号周期时, 执行式
(9) 进行均值 的估算, 给带来的最大误差约为
当采取从均值点 处开始截取整数个信号周期时, 执行式(9) 进行均值的估算, 给带来的最大误差约为
在上述方法中, 由于采样序列 N个点中不恰好含有整数个信号周期(周期长度为),在式(12)式(13)中给累加和所造成的最大误差约为,给和造成的最大误差约为
给造成的最大误差为
由于和具有正交性,和可以认为不相关,该因素给带来的最大误差约为
式(4)可以导出相应的相位误差:
由式(17)可得:
一般说来, 对于具体的周期信号y(t),是确定不变的, 从式(25) 可见, 增大序列样本数 N 可以降低失真度分析的误差, 降低采样间隔( 增高采样速率) 也可以降低失真度分析的误差, 采取从过均值点截取序列周期并经过剔除直流分量后进行分析, 同样将减小失真度分析误差。
当信号y(t)在过均值点附近不连续且出现巨大跳变时,如方波和锯齿波时,设跳变增量为(方波或锯齿波峰值),则上述方法中由于采样序列N个点中不恰好含有整个信号周期(长度为),在式(12)、式(13)中给累加和造成的最大误差约为,给和造成的最大误差约为
此时,从式(25)可见, 增大序列样本数 N 可以降低失真度分析误差, 降低采样间隔对降低失真度分析误差没有明显作用。
通常,每个信号周期至少有 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) 使用了信号周期的精确测量技术, 以代替 N 进行分解运算, 减小了谐波频率定值误差以及所带来的运算误差; (2) 在函数序列的均值点附近开始截取近整数个信号周期的数据 N ,执行谐波分 析, 减少了 非同步采样带来的运算误差和方法误差; (3) 从序列中减去均值分量后进行谐波参数运算, 减少了对各谐波参数带来的影响。采用这些措施后, 并未出现对测量过程的特别要求, 不要求同步采样状态, 实际上是处于非整周期采样状态, 使得上述方法具有较广泛的适应性和简便易行性, 因此, 可望在计量测试和精确测量中有良好的应用前景。
参考文献
[1]潘士先. 电路分析(下册) [M] . 北京: 国防工业出版社,1979, 54- 55.
[2]何天祥. 周期信号各参数的离散采样测量法[J] . 航空计测技术, 1997,17(3):3- 4.
[3]戴先中. 准同步采样及其在非正弦功率测量中的应用[J] . 仪器仪表学报, 1984, 5( 4) : 390- 396.
[4]戴先中, 唐统一. 周期信号谐波分析的一种新方法[J] .仪器仪表学报, 1989,10(3): 248- 255.
[5]戴先中. 进一步提高准同步采 样谐波分析 准确度的两种方法
[J] .仪器仪表学报,1992,13(4):350- 357.
[6]潘文. 准同步采样补偿方法及其误差估计[J] . 仪器仪表学报, 1990, 11(2): 192- 199.
[7]张建秋, 陶然, 沈毅, 邸丽敏. 非整周期采样应用于周期信号的谐波分析[J] . 仪器仪表学报, 1995,16( 1) : 1- 7.
[8]梁志国, 孙瞡宇. 信号周期的一种数字化测量方法.仪器仪表学报( 增刊),2003,24(3):195- 198.
[9]梁志国, 朱济杰, 孙瞡宇. 正 弦信号 总失 真度的 一种 精 确评价方法[J]. 计量学报, 2003, 24( 2) : 144- 148.
[10]梁志国. 动态有效位数的测量不确定度[J].工业计量, 2002, 12( 6) : 46- 49