第5章 CHAPTER 5 取样积分与数字式平均 对于淹没在噪声中的正弦信号的幅度和相位,可以利用第4章中介绍的锁定放大器进行检测 。但是如果需要恢复淹没在噪声中的脉冲波形,则锁定放大器是无能为力的。脉冲波形的快速上升沿和快速下降沿包含丰富的高次谐波分量,锁定放大器输出级的低通滤波 器会滤除这些高频分量,导致脉冲波形的畸变。对于这类信号的测量,必须使用其他的有效 方法,取样积分与数字式平均就是这样的方法。 早在20世纪50年代,国外的科学家就提出了取样积分的概念和原理。1962年,加利福尼亚大 学劳伦茨实验室的Klein用电子技术实现了取样积分,并命名为BOXCAR积分器。为 了恢复淹没于噪声中的快速变化的微弱信号,必须把每个信号周期分成若干个时间间隔,间 隔的大小取决于恢复信号所要求的精度。然后对这些时间间隔的信号进行取样,并将各周期 中处于相同位置(对于信号周期起点具有相同的延时)的取样进行积分或平均。积分过程常 用模拟电路实现,称之为取样积分; 平均过程常通过计算机以数字处理的方式实现,称之为 数字式平均。 多年来,取样积分在物理、化学、生物医学、核磁共振等领域得到了广泛的应用,对于恢 复淹没在噪声中的周期或似周期脉冲波形卓有成效,例如,生物医学中的血流、脑电或心电 信号的波形测量,发光物质受激后所发出的荧光波形测量,核磁共振信号测量等,并研制出 多种测量仪器。对于非周期的慢变信号,常用调制或斩波的方式人为赋予其一定的周期性, 之后再进行取样积分或数字式平均处理。随着集成电路技术和微型计算机技术的发展,以微 型计算机为核心的数字式信号平均器的应用越来越广泛。 5.1取样积分的基本原理 取样积分包括取样和积分两个连续的过程,其基本原理示于图51。周期为T的被测信号s(t )叠加了干扰噪声n(t),可测信号x(t)=s(t)+n(t)经过放大输入到取样开关K。r(t)是与被测 信号同频的参考信号,也可以是被测信号本身。触发电路根据参考信号波形的情况(例如幅 度或上升速率)形成触发信号,触发信号经过延时后,生成一定宽度Tg的 取样脉冲,控制取样开关K的开闭,完成对输入信号x(t)的取样。 图51取样积分基本原理 取样积分的工作方式可分为单点式和多点式两大类。单点式取样在每个信号周期内只取样和 积分一次,而多点式取样在每个信号周期内对信号取样多次,并利用多个积分器对各点取样 分别进行积分。单点式电路相对简单一些,但是对被测信号的利用率低,需要经过很多信号 周期才能得到测量结果; 与此相反,多点式电路相对复杂一些,对被测信号的利用率高,经 过不太多的信号周期就可以得到测量结果。 单点式取样又可以分为定点式和扫描式两种工作方式。定点式工作方式是反复取样被测信号 波形上某个特定时刻点的幅度,例如被测波形的最大点或距离过零点某个固定延时点的幅度 。扫描工作方式虽然也是每个周期取样一次,但是取样点 沿着被测波形周期从前向后逐次移动,这可以用于恢复和记录被测信号的波形。 门积分器是取样积分器的核心,它的特性对于系统的整体特性具有决定性作用。门积分器 不同于一般的积分器,由于取样门的作用,在开关K的控制下,积分仅在取样时间内进行, 其余时间积分结果处于保持状态。根据实现电路的不同,图51中的积分器可以分为线性门 积分器和指数式门积分器。 5.1.1线性门积分 线性门积分器是由线性积分电路附加电子开关组合而成。为了理解线性门积分器的工作原理 ,首先分析图52所示普通线性积分电路的工作过程。 因为图52中左边放大器A1的负输入端为虚地,而且放大器输入端的输入阻抗可近似为无穷大,当输入电压 为ui(t)时,流过输入电阻R的电流i(t)等于流过电容C的电流,即 i(t)=ui(t)R=-Cdvo(t)dt (51) 由式(51)解得 vo(t)=-1RC∫t0ui(t′)dt′+vo0(52) 式中,vo0为A1输出vo(t)的初始电压,积分时间常数为Tc=RC。 当输入电压ui(t)是幅度为Vi的阶跃电压,而且初始电压vo0=0时, 由式(52)解得, 积分器阶跃响应输出为 uo(t)=-vo(t)=Vit/(RC)(53) 式(53)所表示的uo(t)与t之间的线性关系如图53中的点画线所示。 图52线性积分电路 图53线性门积分器的阶跃响应 线性门积分电路示于图54,图中,x(t)为被测信号,它包含有用信号s(t)和噪声n(t),s( t)是周期或似周期信号。r(t)是参考信号,由它触发取样脉冲产生电路,在被测信号周期中 的指定部位产生宽度为Tg的取样脉冲,在Tg期间使电子开关K闭合,以对 被测信号取样。 图54线性门积分器电路 设r(t)的周期为T,取样门闭合时间宽度为Tg,在取样门K的控制下,在r(t)的每个周 期内开关K只在Tg时段内闭合,这时输入电压x(t)经电阻R对C进行积分; 其余时 段开关断开,相当于输入电阻R=∞,电容C两端的电压保持不变。 这时的阶跃响应如图53中的折线所示,该折线可以用一条斜率取决于Tg/T的直线来近似。可以看出,由于取 样开关K的作用,积分的有效时间常数Te>Tc。 由于开关K的开闭作用,门积分器的等效积分电阻为 R(t)=R,门接通时∞,门断开时(54) 设开关闭合的占空因子为Δ=Tg/T,则平均积分电阻可以近似为R/Δ,幅度为Vi 的阶跃响应近似为 uo(t)≈VitΔ/(RC)=VitTg/(TRC)(55) 门积分的等效时间常数为 Te≈Tc/Δ=RC/Δ=RCT/Tg(56) 可见Tg越窄,等效时间常数Te越大。 由于线性门积分电路的输出幅度受到运算放大器线性工作范围的限制,所以比较适用于信号 幅度较小的场合。如果信号幅度较大,为数不多的若干次取样积分就有可能使运算放大器进 入非线性区,导致测量误差,在这种情况下只能使用指数式门积分器。 5.1.2指数式门积分 指数式门积分电路由普通的RC指数式积分器和采样电子开关K串联而成,如图55(a )所示。 图55(a)中的x(t)为被测信号,r(t)为参考输入。由r(t)触发取样脉冲产生电路,以 在被测信号周期的指定部位产生宽度为Tg的取样脉冲,控制电子开关闭合,对被测信 号进行取样。指数式门积分器电路的阶跃响应曲线示于图55(b),图中的点画线是开 关K始终闭合情况下的阶跃响应曲线; 当开关K以周期T、闭合时间宽度Tg周 期性地通断时,电路的阶跃响应曲线如图中的折线所示,这是一种台阶式的指数曲线,其平 均值用虚线示出。 图55指数式门积分器电路及其阶跃响应 (1) 在开关K始终闭合的情况下,图55(a)就是一个普通的RC积分电路,其输入输 出电压关系可以表示为 x(t)=RCduo(t)dt+uo(t)(57) 如果x(t)是幅度为Vi的阶跃电压,而且uo(t)的初始电压为0,则由式(5 7)可解得阶跃响应输出电压为 uo(t)=Vi(1-e-tRC)(58) uo(t)与t之间的关系曲线如图55(b)中的点画线所示,这是一条指数曲线。由 式(58)可以计算出,uo(t)由0上升到0.632Vi所需要的时间为Tc=R C,这就是RC积分电路的时间常数。 (2) 如果r(t)的周期为T,取样门闭合时间宽度为Tg,在取样门K的控制下,只在 Tg时段内输入电压x(t)经电阻R对C进行积分; 其余时段开关断开,输出电压uo (t)保持不变。当x(t)是一个幅度为Vi的阶跃电压时,如果初始条件为零,则阶跃响 应曲线变成图55(b)中的折线。所以,由于取样门K的开关作用,当开关闭合时积 分电阻为R,开关断开时,积分电阻为∞,设开关闭合的占空因子为Δ=Tg/T,则平均积 分电阻为R/Δ。阶跃响应可近似表示为 uo(t)=Vi(1-e-tΔRC)(59) 式(59)表示的uo(t)与t之间的关系如图55(b)中的虚线所示,该虚线可以看 成是图中折线的平均值。在这种情况下,令uo(t)由0上升到0.632Vi所需的等 效时间常数为Te,则 Te=Tc/Δ=TcT/Tg=RCT/Tg(510) 可见,取样门开关的作用使得积分的时间常数加长了很多。 与线性门积分相比,指数式门积分有利有弊。由图55(b)可以看出,随着取样次数的 增加,每个取样使积分输出上升的值逐渐减少。经过5倍的Te后接近稳定值,此后的 取样对积分输出影响很小,因此不会因为积分时间太长而过载。另一方面,当积分时间大于 2Te后,每次取样使得积分结果变化很小,而且会越来越小,即积分作用降低。这意 味着太长地增加测量时间是没有意义的,因为在2Te之后,继续采样积分对提高信噪 比作用不大。 相比之下,对于线性门积分,信噪比的改善会随着积分时间的增加而增加,它不受电路等效 时间常数的限制,只受电路工作线性范围的制约。所以在信号幅度较小的情况下,采用线性 门积分更为有利。而在信号幅度较大时,为了防止电路进入非线性区导致测量误差,必须采 用指数式门积分器。所以,在具体的门积分应用中,要根据实际检测情况和要求选择合适的 门积分方式。 5.2指数式门积分器分析 取样积分的关键部件是门积分器,取样积分抑制噪声的能力及其他一些重要的性能指标也主 要取决于门积分器的性能。在5.1节介绍的两种门积分器中,指数式门积分器不易因输入x( t)幅度过大而过载,抵御大幅度噪声的能力较好,所以得到了广泛的应用。 图55(a)所示的指数式门积分器可以分解为两级电路相串联而成: 前一级是取样门电 路,它以周期T、脉冲宽度Tg对输入信号进行取样; 后一级是RC积分电路。下面分别 分析它们的传输特性。 5.2.1取样过程频域分析 取样过程就是利用取样脉冲序列p(t)从被测连续信号x(t)中“抽取”一系列的离散样值,如图56所示,取样电路输出xs(t)可以看作是取样脉冲 图56取样门与取样脉冲序列 序列p(t)与连续信号x(t)的乘积,即 xs(t)=p(t)x(t)(511) 设取样脉冲序列p(t)的幅度为A,周期为T,脉冲宽度为Tg, 这样的周期性脉冲序 列可以展开为傅里叶 级数的形式 p(t)=∑∞n=-∞cnexp(jnωst)(512) 式中,ωs=2π/T为取样脉冲p(t)的基波角频率,复数集cn称为p(t)的频谱 cn=1T∫T/2-T/2p(t)exp(jnωst)d t=ATgTsin(nωsTg/2)nωs Tg/2(513) 令Δ=Tg/T,Δ为p(t)的占空系数,考虑到用p(t)控制取样开关的情况相当于p(t)的 幅度A=1,而且ωs=2π/T,则由式(513)可得 cn=Δsin(nπΔ)nπΔ(514) 将式(514)代入式(512)得 p(t)=∑∞n=-∞Δsin(nπΔ)nπΔ exp(jnωst)(515) 对式(515)进行傅里叶变换得p(t)的频谱为 P(ω)=∑∞n=-∞Δsin(nπΔ)nπΔ δ(ω-nωs)(516) 式中,ωs=2π/T为取样脉冲的角频率。可见,P(ω)的图形是包络线为取样 函数(sample function)sin(nπΔ)/(nπΔ)的一系列的冲激函数,如 图57(a)所示。取样函数在某些文献中又称之为sinc函数。 图57矩形脉冲取样信号的频谱变换 设输入信号x(t)的频谱为X(ω),如图57(b)所示。根据傅里叶变换的性质,式(51 1)所表示的p(t)和x(t)的相乘过程在频域表现为两者频谱的卷积,即 Xs(ω)=X(ω)*P(ω)(517) 式中,Xs(ω)表示取样信号xs(t)的频谱。将式(516)代入式(517),可 得 Xs(ω)=X(ω)*∑∞n=-∞Δsin(nπΔ)nπΔδ(ω-nωs) =∑∞n=-∞Δsin(nπΔ)nπΔX(ω- nωs)(518) 或者将式(515)代入式(511),可得 xs(t)=∑∞n=-∞x(t)Δsin(nπΔ)nπ Δexp(jnωst)(519) 取式(519)两侧的傅里叶变换,也可以得到式(518)的结果。 由式(518)可见,取样过程的作用是将输入信号x(t)的频谱X(ω)平移到nωs各点 (n为-∞~+∞的整数),再分别乘以相应的cn值。 图57说明了上述变换过程。图57(a)表示的是取样脉冲序列p(t)及其频谱P(ω),图 57(b)表示的是时域信号x(t)及其频谱X(ω),p(t)与x(t)相乘得到图57(c)中 的取样信号xs(t),X(ω)与P(ω)卷积的结果得到图57(c)中取样信号的频谱 Xs(ω),Xs(ω)是幅度按取样函数包络分布的离散频谱。 5.2.2指数式门积分器的输出特性 任何波形的周期信号都可以表示为三角函数的组合,考虑输入被测信号中频率为ω的单一频率正弦信号分量为 x(t)=xmcos[ω(t-τ)](520) 式中,xm是该频率分量的幅度。参考文献[12]中给出了下列几个有用的结论。 1. 当ω=ωs=2π/T时 这意味着x(t)的频率ω等于取样脉冲频率,这时指数式门积分器的稳态输出为 uo(t)=xmsin(πΔ)πΔcos(ωs τ)1-exp-ΔRCt(521) 式中,Δ=Tg/T为取样脉冲p(t)的占空系数。 由式(521)可以得出以下结果: (1) 式中的指数项说明,取样积分器的输出是沿着指数曲线逐渐积累的过程,时间常数为T e=RC/Δ,图55(b)也说明了这一点。取样脉冲宽度Tg越窄,占空系数Δ 越小,则所需积分时间越长。 (2) 当改变延迟时间τ时,输出按cos(ωsτ)的规律变化; 当τ从0逐渐变化到T 时,输出显示出一个完整周期的正弦波。在ω=ωs情况下,式(521)中的τ可以 看成是取样脉冲相对于被测信号x(t)起始点的延时。这说明,逐渐改变取样点相对于被测信 号x(t)起始点的时间,可以从噪声中恢复出被测信号的波形,这正是扫描式取样积分器的工 作原理。 (3) 取样积分器稳态输出时的衰减系数取决于sin(πΔ)/(πΔ),若要求衰减系 数小于3dB,则要求 sin(πΔ)πΔ<12 由上式解得πΔ<1.392,由Δ=Tg/T可得 Tg<0.4431T(522) 式(522)说明,为了使取样积分的稳态值衰减不多,取样脉冲宽度Tg要足够窄,如 果Tg>0.4431T,则衰减系数有可能大于3dB。 2. 当ω=nωs时 这意味着x(t)的频率ω等于取样脉冲p(t)的某次谐波频率,这时指数式门积分器的稳态输出为 uo(t)=xmsin(nπΔ)nπΔcos(nω sτ)1-exp-ΔRCt(523) 当n=1时,式(523)就退化为基波情况下的式(521)。由式(523)可得如下结论。 (1) 当改变延迟时间τ时,输出按cos(nωsτ)的规律变化,这说明通过改变延迟 时间τ可以恢复被测信号的任何高次谐波分量的波形。 (2) 式(523)右边中括弧中的项说明,输出信号中的高次谐波分量也是按指数规律逐渐积 累的过程,积分的时间常数也是Te=RC/Δ。 (3) 为了使恢复的被测信号的n次谐波分量衰减系数小于3dB,要求 sin(nπΔ)nπΔ<12 由上式解得nπΔ<1.392,由Δ=Tg/T可得 Tg<0.4431T/n(524) 式(524)说明,取样脉冲宽度越窄,输出信号的分辨率越高。要想使恢复信号的n次谐波分 量衰减系数小于3dB,取样脉冲的宽度必须小于n次谐波周期的0.4431倍。 3. 当ω=nωs+Δω时 这意味着x(t)=xmcos[ω(t-τ)]的频率ω偏离取样脉冲的n次谐波nωs一个Δω,这时指数式门积分器输出的稳态振幅为 uom=xm·sin(nπΔ)nπΔ·1 1+(ΔωRC/Δ)2(525) 式(525)除以x(t)的振幅xm,就可以得到在ω=nωs附近的幅度响应 |Hn(ω)|=sin(nπΔ)nπΔ·11+( ΔωRC/Δ)2(526) 式中,Δω=ω-nωs。 式(526)说明,经过取样积分,输入被测信号在各谐波处要经过一阶带通滤波,带宽取 决于等效时间常数Te=RC/Δ。而且,各次谐波处的幅度按取样函数sin(nπ Δ)/(nπΔ)分布。 当取样脉冲的占空比Δ<<1时,有ΔωRC/Δ>>1,nπΔ<<1,那么式(526)可以简化 为 |Hn(ω)|=ΔΔωRC(527) 将式(526)推广到各次谐波的情况,将各次谐波处的幅度响应求和就得到指数式取样 积分器总的幅度响应|H(ω)| |H(ω)|=∑∞n=-∞sin(nπΔ)nπΔ·11+[(ω-nωs)RC/Δ]2(528) 根据式(528),当Δ=Tg/T=0.2时,取样门电路的幅频响应|H(f)|示于图5 8。可以看出,|H(f)|是一个幅度服从取样函数规律的离散频域窗。 图58取样门电路的输出特性 由式(528)和图58可以得出以下几点结论。 (1) 在取样频率fs的各次谐波处的带宽随积分时间常数RC的增加而减少,随占空比Δ= Tg/T的减少而减少。也就是说,其带宽取决于Δ/(2πRC)。 (2) 在fs的各次谐波处的通带幅度服从取样函数sin(nπΔ)/(nπΔ)的 规律。 5.2.3指数式门积分的信噪改善比 1. 污染噪声为白噪声时的信噪改善比 由图55(b)的响应曲线可知,对于时间常数为Tc的指数式积分器,当有效积分 时间接近5Tc时,积分效果变得不明显。设积分次数为N,每次取样积分时间为T g,则有效积分时间为NTg。为了使积分作用有效,应该保证 NTg<5Tc(529) 在本章后面的数字累加部分将会看到,对于污染噪声是白噪声的情况,取样累积N次所 能实现的信噪改善比(有效值)为 SNIR=N 将式(529)代入上式得 SNIR<5Tc/Tg(530) 式(530)只是给出了信噪改善比的大致范围。用数学分析的方法确定准确的信噪改善 比SNIR比较繁琐,根据文献[15]中的分析,对于污染噪声是白噪声的情况,指数式门积分 可以达到的信噪改善比(有效值)为 SNIR=2Tc/Tg(531) 2. 污染噪声为有色噪声时的信噪改善比 对于污染噪声为有色噪声的情况,例如白噪声通过RC低通滤波器后的输出x(t)就是一种有色噪 声,根据式(165a),其自相关函数可以表示为 Rx(τ)=Pnexp(-|ατ|)(532) 式中,Pn=Rx(0)为噪声的功率,α=1/(RC)为有色噪声的相关函数指数因子,则指数式门积分可以 达到的信噪改善比SNIR(有效值)为 SNIR≈2Tc/Tg1+2e-αT (533) 式中,T为取样周期。对于白噪声的情况,相当于α=∞,式(533)可简化为式(531)。 5.3取样积分器的工作方式 取样积分器的工作方式可分为定点式和扫描式两种,一般将这两种工作方式组合在同一仪器 中,由用户选择使用哪种工作方式。定点工作方式用于检测信号波形上某一特定位置的幅度 ,而扫描工作方式用于恢复和记录被测信号的波形。 5.3.1定点工作方式 在定点工作方式中,参考触发信号与输入被测信号保持同步,经过延时后产生固定宽度为T g的门控信号,这样取样积分总是在被测信号周期的固定部位进行。定点工作方式 比较简单,适用于检测处理周期信号或似周期信号固定部位的幅度,例如接收斩波光的光电 倍增管(PMT)的输出电流,心电图一定部位(例如R波)的幅度等。 图59所示为定点取样积分电路原理方框图,它由信号通道、参考通道和门积分器组成。 信号通道中的前置放大器为宽带低噪声放大器,用于将叠加了噪声的微弱被测信号x(t)放大 到合适的幅度。参考通道由触发电路、延时电路和取样脉冲宽度形成电路组成。参考信号可 以是与被测信号相关的信号(例如交流电桥测量电路的交流电源),也可以是被测信号本身 。当参考信号的一定特征(例如幅度或变化率)达到一定数值时,则产生触发信号,触发信号 经过延时后触发门控电路,以形成宽度为Tg的取样脉冲,在被测信号周期中的固定部 位进行取样和积分。延时电路的延时量可调,以便调整取样的部位。对取样积分输出信号u o(t)可以进一步放大,以便于观测或记录。 图59定点取样积分原理示意图 定点工作方式中的各点波形示于图510。图中,参考触发信号经过一定时间的延迟Td 之后,电路产生宽度为Tg的取样脉冲,对被测信号的固定部位进行取样积分。可以 看出,在每个信号周期内,取样积分只进行了一次,在Tg期间取样并积分,而在其他 时间开关K断开,保持积分结果,输出信号呈现阶梯式累积的波形。经过多个周期的取 样积分,输出信号趋向于被测信号取样点处的平均值。 图510定点工作方式取样积分的各点波形 在定点工作方式中,因为取样点相对于信号起始时刻的延时是固定的,取样脉冲宽度Tg也保持不变,所以取样总是在被测信号距离原点为固定延时的某个小时段重复进行,积分 得到的结果是该时段的多次累加积分值。利用信号的确定性和噪声的随机性,重复取样积分 的结果将使信噪比得以改善。 图511所示是一种比较特殊的定点差值取样积分电路[12],图中的触发整形电路 和延时电路用于产生相对于参考信号原点固定延时Td的取样脉冲,其宽度Tg由 取样脉冲宽度控制电路设定。电阻R、电容C和A2组成积分器,放大器A1和A3组成差值积分电 路。被测信号经前置放大与上次取样积分结果的分压值相比较,在A1的输出端得到差值信 号,该差值信号被送到取样门进行定点取样,再经积分器积分得到输出信号。在电子开关 K接通期间,积分器对A1输出进行积分; 在电子开关K打开期间,由于运算放大器A 2的输入阻抗很高,积分器保持上次的积分结果。图中的R1、R2以及A3组成反馈支 路,A1将当前的前置放大输出与上次取样积分的输出进行比较,输出给取样门K的电 压为 u1(t)=A1A0x(t)-R2A3R1+R2·uo(t)(534) 所以这是一种差值取样积分,其工作原理类似于密勒积分器。 图511定点差值取样积分电路 5.3.2扫描工作方式 定点式取样积分器只能用于测量周期或似周期信号固定部位的电压,却不能用于恢复被测信 号的整个波形。在取样积分器的扫描工作方式中,取样点距离波形原点的延时量被逐渐延长 。随着一个个信号周期的到来,取样点沿着信号周期波形从前向后进行扫描,从而恢复被噪 声污染的波形。 扫描式取样积分器的结构方框图示于图512,图中的慢扫描电路用于产生覆盖很多个信号 周期的锯齿波,其宽度为Ts; 时基电路用于产生覆盖被测信号周期中需要测量部分的 锯齿波,其宽度为TB。比较器电路对两个锯齿波进行比较,从而产生逐渐增加的延时 ,这样就可以在被测信号的逐个周期中从前向后延时取样,以便实现对原信号的逐点恢复。 门控电路用于产生宽度为Tg的取样脉冲。 图512扫描式取样积分结构框图 图513示出扫描式取样积分器的各点波形,(a)为被测信号x(t)波形; (b)为 由参考信号r(t)产生的触发信号波形; (c)为由慢扫描电路产生的长周期Ts锯 齿波(虚线),以及覆盖被测信号需要测量部分的时基TB锯齿波(实线); 比较器根 据两个锯齿波的相交点产生延时脉冲,如图513(d)所示,其上升沿相对于触发脉冲 的延时逐渐增加; 由延时脉冲触发门控电路产生逐次后移的取样脉冲,如图513(e) 所示。可以看出,随着被测信号周期的逐个到来,取样点在信号周期中的位置从前向后逐次 移动,在经历了很多个信号周期后(为简单计,图513只画出5个周期),由取样值的包络 图513扫描式取样积分各点波形 线可以显现被测信号的波形,不过周期比原信号长了很多倍,如图513(f)所示。因 此,利用XY记录仪可以记录显示被测信号的波形。把记录仪的X输入端连接到慢扫描 锯齿波输出,Y输入端连接到取样积分输出信号端就能实现这种记录。 设取样脉冲相对于信号周期起始点的延迟量在每个信号周期中增加Δt,如果时基锯齿波宽度为T B,慢扫描锯齿波宽度为Ts,时基锯齿波与慢扫描锯齿波的幅度相同,则根据图5 13(c)中的几何关系,可得 Δt=TBTTs-TB(535) 式中,T为信号周期。考虑到Ts>>TB,式(535)可简化为 Δt≈TBTTs(536) 扫描式取样积分的工作过程是一种移动平均式的积分,如图514所示。图中,x(t)是 被测信号波形,取样积分器对Tg时段内(虚线框内)的x(t)进行积分, 图514扫描式移动取样积分示意图 得到一个u o(t)输出值,在信号的下一个周期虚线框向右移动一个小小的时段Δt,再次进行积分 。重复上述过程直到扫描完要测量的时段,就像积分框沿着信号周期向前移动一样,所以这 种积分方式又叫作运转积分(running integration),或称Boxcar积分。 任何信号都可以分解成不同频率分量的组合,设x(t)的频率为ω的分量为 xi(t)=Vmexp(jωt)(537) 式中,Vm为该频率分量的幅度,则在时刻t的积分平均结果为 uo(t)=1Tg∫tt-TgVmexp(jω t′)dt′ =Vmexp(jωt′)jωTgt t-Tg=Vmexp(jωt)jωTg[1-exp(- jωTg)] =xi(t)jωTg[1-cos(ωTg)+jsin(ω Tg)] =xi(t)exp(-jωTg/2)sin(ωTg/2)ωTg /2(538) 对比式(537)和式(538)可得,取样积分的频率响应为 |H(ω)|=|uo||xi|=sin(ωTg/2 )ωTg/2(539) 上述结论也可以由下面的数学推导得出。对于图514所示移动取样积分过程,时刻t积分器的输出为 uo(t)=1Tg∫tt-Tgx(u)du 式中的Tg为固定积分时间。积分器在时刻t-τ的输出为 uo(t-τ)=1Tg∫t-τt-τ-Tgx(v)dv 积分输出uo(t)的自相关函数为 Ruo(τ)=E[uo(t)uo(t-τ)]=E1Tg∫tt-Tgx(u)du·1Tg∫t-τt-τ-Tgx(v)dv =1T2g∫tt-Tgdu·∫t-τt-τ-TgdvE[x(u)x(v)] =1T2g∫tt-Tgdu·∫t-τt-τ-TgdvRx(u-v)(539a) 根据式(137),有 Rx(u-v)=12π∫∞-∞Sx(ω)ejω(u-v)dω 将上式代入式(539a),并改变积分顺序,得 Ruo(τ)=1T2g·12π∫∞-∞Sx(ω)dω∫tt-Tgejωudu·∫t-τt-τ-Tge-jωvdv(539b) 上式中,对u的积分为 ∫tt-Tgejωudu=ejωt(1-e-jωTg)jω 对v的积分为 ∫t-τt-τ-Tge-jωvdv=e-jω(t-τ)(1-ejωTg)-jω 它们的乘积为 ejωτsin2(ωTg/2)(ω/2)2 代入式(539b),得 Ruo(τ)=12π∫∞-∞Sx(ω)sin2(ωTg/2)(ωTg/2)2ejωτdω 根据式(137)所表示的自相关函数与功率谱密度的关系,由上式的被积函数可以看出,uo(t)的功率谱密度为 Suo(ω)=Sx(ω)sin2(ωTg/2)(ωTg/2)2 与式(160)对比,可知Boxcar积分的频率响应函数为 图515线性取样积分的幅频响应曲线 |H(ω)|=sin(ωTg/2)ωTg/2(539c) 上式与式(539)的结果相同。据此可画出取样积分器的幅频响应曲线,如图515所示。 由时基电路产生的TB锯齿波的起始点及斜率都可以根据需要进行调节,这有利于选择 恢复波形的位置和宽度。例如,图516(a)中的TB覆盖了被测波形的大部分,而图 516(b)中TB只覆盖了被测波形的中间部分。 图516TB变化与波形恢复的关系 在实际仪器中,一般都把图59电路和图512电路组合在一起,由使用者选择定点或扫描 工作方式,如图517所示。当开关K1打到上部时为单点工作方式; 当K1打到 下部时为扫描工作方式。 图517具有单点和扫描两种工作方式的取样积分器结构框图 5.4取样积分器的参数选择及应用 5.4.1取样积分器的参数选择 1. 取样脉冲宽度Tg的选择 取样脉冲宽度Tg不能选得太宽,否则会造成信号中高频分量的损失,使得恢复的信号 失真。下面以正弦信号为例,说明取样脉冲宽度Tg的选择原则。 对于图518所示的正弦波x(t)=Vmsinωt,以取样脉冲宽度Tg对x(t)取样 ,设取样脉冲的中心时刻为t0,则信号取样后的输出电压为 x(t)=Vmsinωt,t0-Tg2≤t≤t0+T g2 经积分器积分后的输出为 uo(t0)=∫t0+Tg2t0-T g2Vmsinωtdt=2VmωsinωTg2sin(ωt0)(540) 图518正弦波定点取样 式(540)是对任何ω都适合的结果。当频率ω很低时,ωTg→0,sin(ωT g/2)≈ωTg/2,式(540)可以近似为 uo(t0)≈2VmωωTg2sin (ωt0)(541) 对比式(540)与式(541)可知,当频率较高时,因为sin(ωTg/2)<ωTg/2,积分输出电压会下降,所以会引起信号中高频分量的损失,损失程度可表示 为 A=uo|ωuo|ω→0=sin(ωTg /2)ωTg/2(542) 式(542)与式(539c)的结果相同。式(542)说明,取样积分对被测信号高频分量的 衰减系数A与fTg相关, 图519取样积分衰减系数A与 fTg的关系 根据式(542)可画出A与fTg的关系,如图519所示 。若要求取样积分对被恢复信号的最高频率fc的衰减不大于3dB,即要求A≥1/2,则由式(542)得 fcTg≤0.42 或 Tg≤0.42/fc(543) 可见,希望恢复的信号频率越高,要求取样脉冲宽度Tg越窄。式(543)与式 (522)及式(524)对Tg的要求基本一致。 如果要恢复的信号波形包含很陡的上升沿或下降沿,则取样脉冲宽度Tg必须很窄。根 据脉冲技术和宽带放大器的分析,脉冲上升时间tr和频带宽度的关系为[12] fc≈0.35/tr(544) 由式(544)可见,信号波形的上升沿或下降沿越陡,则信号的最高频率fc越高,Tg应该越窄。但是Tg也不能选得太窄 ,从后面的分析可知,Tg越窄,测量时间就越长。式(543)和式(544)可以用 作取样脉冲宽度Tg的选择原则之一,实际应用中还要与测量时间综合权衡 。 2. 时基锯齿波宽度TB的选择 时基锯齿波的起始点及斜率都可以根据需要进行调节,TB的范围取决于被测信号周期 中需要恢复的区段长度,考虑到各种不确定因素,选择TB时,要在测量区段的两端都留 有一定的余地。 3. 积分器时间常数Tc=RC的选择 根据取样积分器的工作原理,对于指数式取样积分器,在每个取样脉冲作用期间,取样开关 闭合对积分电容充电,充电时间为Tg; 而在两次取样脉冲间隔期间,电容电压保持不 变,如图55(b)所示。对于指数式门积分器,当N次取样总的积分时间NT g接近5倍的积分器时间常数时,信号累积速度减慢,信噪比改善很少。在两倍的时间常数 内,信噪比改善得比较明显。根据式(531)可得 Tc=(SNIR)2Tg2(545) 而对于线性门积分器,信噪比的改善不受积分器时间常数的限制,仅受电路的动态范围限制 。取样积分次数N越大,信噪改善比SNIR越大,两者之间的关系为 SNIR=N(546) 这时就要根据取样脉冲宽度、测量范围、慢扫描测量时间的要求,综合考虑来确定积分时间 常数。 4. 慢扫描时间Ts的选择 设取样脉冲Tg相对于信号周期起始点的延迟量在每个信号周期中增加Δt,在许多个被测信 号周期内,Tg相对于被测波 图520取样脉冲相对于被测信号 周期起始点移动示意图 形移动的情况示于图520。被测信号的每一点只有在取 样脉冲宽度Tg内才能被取样,而Tg又以时间间隔Δt在跳跃移动。所以,对于 被测信号的任何一点,被取样次数为 N=TgΔt(547) 将式(536)代入式(547)得 N=TgTsTBT(548) 将式(548)代入式(546)得 SNIR=TgTsTBT(549) 对于线性门积分器,给定所要求的SNIR后,由式(549)可以估计出需要的测量时间Ts Ts≥(SNIR)2TBTTg(550) 在指数式门积分器中,为了使电容充电充分,须使总的积分时间 NTg 比积分器 的时间常数Tc大很多。当NTg=5Tc时,积分器充电值与稳定值之间的误 差为0.67%。若要求总的积分时间 NTg≥5Tc 则将式(548)代入上式得 T2gTsTBT≥5Tc 或 Ts≥5TBTcTT2g(551) 若按式(545)选择积分器时间常数Tc,则 Ts≥2.5TBT(SNIR)2Tg(552) 由式(551)和式(552)可见,波形恢复测量的总时间Ts正比于要恢复波形 宽度TB和信号周期T,这是比较容易理解的。要达到的信噪改善比SNIR对测量时间影 响很大,二者为平方关系。 取样脉冲宽度Tg的选择要考虑两方面的因素,即既要满足被测信号频率分辨率的要求 (见式(543)),又要考虑测量时间的问题。如果所选Tg太窄导致测量时间Ts 太长,那么电容的漏电和放大器的漂移也会引起测量误差,所以需要综合权衡两方面的因素 。 例51被测信号周期T=10ms,测量范围TB=2ms,要求信噪改善比SNIR=10,选TB=3ms以覆盖测量区,Tg=100μs,试求指数式取样积分器的参数Tc和Ts。 解: 根据式(545)及所要求的SNIR求出积分器时间常数Tc如下 Tc=(SNIR)2Tg2=5(ms) 再根据式(551)计算总的测量时间 Ts≥5TBTcTT2g=75(s) 综上所述,线性取样积分器和指数式取样积分器的参数选择过程分别如图521(a) 和图521(b)所示。 图521取样积分器参数选择过程框图 5.4.2基线取样与双通道取样积分器 1. 基线取样 利用取样积分的方法来改善被测信号的信噪比是以时间为代价的。因为在被测信号的每个周 期只取样一次,为了改善信噪比,必须对很多个信号周期进行取样积分,这就对系统的稳定 性提出了很高的要求。当被测信号周期较长或要求的SNIR较大时,必须进行长时间的测量, 才能达到要求。而在长时间的测量过程中,由于电容漏电、放大器的零点和增益变化、其他 元器件的温度漂移和时间漂移、激励源的起伏等因素,被测信号的零点基线已经发生了变化 ,称之为基线漂移。 图522基线漂移示意 图522示出没有漂移的信号波形和漂移后的信号波形。可以看出,消除基线漂移对于信号 的精确测量是十分重要的。采用低漂移器件,对器件进行老化筛选,改进电路设计等措施对 改善漂移都有一定的作用,但实际实现起来往往具有一定的困难,改善效果也很有限。而基 线取样方法对于改善取样积分测量的漂移是行之有效的。 基线取样的基本原理是: 在每个信号周期内先取样一次信号的有效成分,再取样一次信号的 基线,两者相减得到扣除漂移后信号成分的有效幅值。 具有基线取样补偿的取样积分器的结构框图示于图523,它包括两个相同的取样积分器, 由时基控制实现在不同时刻对信号和基线的交替选通,通过积分器1和积分器2分别对信号和 基线的取样值进行积分,两者相减就在输出中消除了基线漂移。 图523具有基线取样补偿的取样积分器结构框图 图524示出基线取样的工作波形。图524(a)所示为被测斩波信号波形,其低电压处 的起伏反映了背景噪声和基线漂移。在时基控制电路中,参考输入在每个信号周期产生一个 触发脉冲,如图524(b)所示; 由触发脉冲触发产生一个覆盖被测信号周期的锯齿波, 如图524(c)所示; 利用比较器将锯齿波与交替的固定电压V1和V2相比较,从而 在时刻t1和t2分别产生取样脉冲A和B,如图524(d)所示; 用取样脉冲 A和B分别对基线和有用信号进行取样,并分别积分。长时间的A点和B点交替 取样积分的结果为,积分器1输出为信号的累积平均值,积分器2输出为基线的累积平均值, 两者相减就实现了对基线漂移的补偿。 图524基线取样工作原理 如果V1为固定电压,V2用慢扫描锯齿波Ts来代替,则上述交替取样积分过程对于 基线为定点取样积分,而对于信号为扫描式取样积分,这样的取样积分可以用来恢复被测信 号所需部分的波形,又克服了基线漂移的不利影响。 2. 双通道取样积分器 双通道取样积分器系统由两路基线取样积分器组成,如图525(a)所示。它利用两路取 样积分器分别对被测信号和作为标准的另一路信号进行消除基线的取样积分,将被测量与标 准量对比,以消除由信号基线以及激励源起伏造成的误差,使测量更为准确和可靠。 图525双通道取样积分器 信号通道对被测信号及其基线进行取样积分,以消除信号基线漂移的影响,输出为A。标准通 道对用作激励源的标准信号及其基线进行取样积分,消除标准通道基线漂移的影响,输出为 B。之后进行A/B的运算,对测试结果进行归一化处理,从而克服激励源起伏引起的测量误差 。图525(b)示出对信号和标准分别取样的波形。为了实现各取样点的准确定时,确定 取样脉冲宽度等功能,系统还必须包括时基和延时控制部件,在图中未示出。 例如在激光测量中,当测量时间较长时,激光光源的强度可能发生变化,从而导致测量误差 。在这种情况下可以利用A通道作基本测量,利用B通道检测光源的强度,从而补偿光源强度 变化带来的影响。 除了常用的A/B运算外,功能插件还可以选择A输出、B输出、A-B、A×B、lg(A/B)等功能,给测试带来很大方便。 5.4.3多点取样积分器系统 上面介绍的定点式和扫描式取样积分器都属于单点式取样,即每个信号周期内只取样积分一 次,所以取样效率很低,需要经过很多信号周期才能得到测量结果。当被测信号重复频率较 低时,必然导致检测时间太长。尤其是当被测信号包含快速跳变沿时,为了恢复被测信号的高频 分量,根据式(543),取样脉冲宽度Tg必须很窄。这样一来,在扫描式工作方式中 ,为了恢复被测信号波形,由式(551)计算出的测量时间Ts会很长,电容漏电和电 路漂移都会造成测量结果的误差。 例如,对于10Hz的被测信号,若选Tg=0.1ms,要求SNIR=100,全周期恢 复波形则要求TB=T=0.1s,那么按式(552)计算出的测量时间Ts≈29 天,这当然是很不现实的。 多点式取样积分器是在每个被测信号周期内取样多点,其电路相对要复杂一些,但是对被测 信号的利用率高,经过不太多的信号周期就可以得到测量结果。图526示出多点式取样积 分器的电路结构,它包括多个电子开关和积分电容,相当于多个单点取样积分器组合在一起 。在取样时间控制逻辑的控制下,各个电子开关轮流闭合,在每个取样脉冲Tg期间只 有一个开关闭合,M个开关轮流闭合一次,扫描完一个被测信号周期中要测量的时段。当下 一个信号周期到来时,重复上述工作过程。这样就可以在信号周期的不同时段对不同的电容 充电,达到多点取样积分的目的。 图526多点取样积分器电路结构 多点取样积分器工作波形示于图527,图中,T为信号周期,Tg为取样脉冲宽度,每 个信号周期只画出了3个取样间隔和3个相应的积分器。注意,图526中轮流闭合的电子开关具有两个 作用: 一个作用是在信号周期的不同取样时段开关切换到不同的电容进行积分; 另一个作用 是使输出放大器的输入端也轮流切换到不同电容的积分电压上,也就是实现图527中轮流 输出选择的功能。因此,输出放大器A2的输出波形就反映了扫描式取样积分的结果。经过 多个信号周期的扫描式取样积分,输出波形的信噪比将得以改善,所实现的SNIR正比于信号 周期数i的平方根。 图527多点取样积分器工作波形示意图 图526电路正常工作的关键是,在信号周期中需要测量的时段内,取样时间控制逻辑电路 依次给各个电子开关输出取样脉冲,其余电路比较容易实现。这种多点取样积分器成本 低、取样效率高、工作速度快。但是,由于每个取样点都需要一套单独的电子开关和积分电容, 所以每个信号周期中的取样点数不可能太多,一般为50~100点。此外因为电容存在漏电问 题,所以信号电压保持时间有限,很难实现对低频信号的恢复。 5.4.4取样积分器应用实例 取样积分器已广 泛应用于物理、化学、生物学以及工业检测技术等许多领域,国内外已经研制开发出若干种 实用的取样积分检测仪器。这里只列举几种比较典型的应用实例。 1. 材料的光学特性检测 为了检测某种材料的光学特性,例如光学材料对某种特定波长的单色光的吸收特性,可以利 用马达带动等距开槽的斩光盘,使经过斩波的单色光透射过被测材料样品,如图528所示 。光电倍增管PMT将透射过被测材料的斩波光转换为电信号,其中包含杂散光、漏电流 、暗电流等各种噪声。利用定点式取样积分提高输出信号的信噪比,并利用基线取样补偿基 线的漂移,就可以根据输出信号判别被测材料的吸收特性。在斩波盘两侧装设小灯和光敏器 件,光敏器件输出的脉动信号可以用作取样积分的参考信号。 图528利用斩波光检测材料的光吸收特性 2. 霍尔效应测量 对于如图529所示的半导体薄片,若在它的两端通以控制电流I,并在薄片的垂直方向施加 磁感应强度为B的磁场,则在磁场和电流的垂直方向上将产生电动势VH, 图529霍尔效应原理图 这种现象称 为霍尔效应(Hall effect)。 霍尔效应的产生是由于半导体中的运动电荷受洛伦兹(Lorentz)力作用的结果。设在N型半 导体薄片的控制电流端通以电流I,则半导体中电子的运动方向与电流方向相反。如果在 垂直于半导体薄片的方向上施加磁场B,洛伦兹力的作用会使电子向一边偏转,见图529中的虚 线箭头方向。这会导致该边的电子堆积,而另一边则积累正电荷,于是产生电场。该电场阻 止运动电子的继续偏转,当电场作用力与洛伦兹力相等时,电子的积累就达到动态平衡。这 时,在薄片积累了电荷的两个端面之间建立的电场称为霍尔电场,相应的电位差称为霍尔电 势VH,当I与B垂直时,VH的大小为 VH=RHIB/d(V) 式中,RH为霍尔常数,取决于材料的物理性质,m3/C; I为控制电流,A; B为磁感应强度,T; d为霍尔元件的厚度,m。由上式可知,对于给定的霍尔元件,如果给定I和B中的任一项,就可以由VH测出另一 项,而且检测的特性曲线是线性的。在实际应用中常使用霍尔元件检测位移或磁场。 图530所示是利用霍尔元件检测磁场强度,由于干扰磁场的存在,霍尔元件的输出电势常 常夹杂着明显的噪声。利用取样积分器抑制噪声,可以使输出信噪比得以改善。霍尔元件输 出电势中还包含了不等位电势、寄生直流电势、 图530利用取样积分器检测霍尔电势 感应电势、温度误差等因素,这些因素都会 影响测量结果的精确度。如果霍尔元件的激励电流I是频率为f的方波,这相当于对磁场强度 进行斩波测量,那么前面介绍的基线取样方法可以用来补偿检测元件固有的偏差以及基线的 漂移。 3. 利用超声波检测材料特性 超声波是频率高于20kHz的机械振动波。当超声波入射到某种物体中时,部分超声波被 反射,部分超声波被吸收,部分超声波透射过该物体。当超声发射源与反射物体之间存在相 对运动时,接收到的超声频率会发生偏移,即出现多普勒效应。超声波的这些性质可以用来 检测多种物理量,例如距离、流速、料位等,还可以用来探测构件中的裂缝、气泡或其他缺 陷,称之为超声波探伤。此外,超声波在医学诊断领域也获得广泛应用。 超声波发送器和接收器一般都是由压电材料构成,在受到电脉冲激励时,超声波传感器会发送一 串短时段逐渐衰减的超声波,频率等于其共振频率。如果激励电压是频率等于其共振频率的 持续正弦波,则超声波传感器会发送持续的超声波; 如果超声波传感器接收到等于其共振频率的 超声波,它又可以把超声波转换为电信号。 利用超声波检测材料特性的原理框图示于图531。图中的脉冲发生器产生不断重复的短时 脉冲,该短时脉冲经驱动电路激励超声波发送探头发送断续的但波形重复的超声波。超声波接收 探头接收透射过被测材料的超声波,并将其转换为比较微弱的电信号,电信号经放大后再由 取样积分器提高信噪比。通过对被测材料声波速度特性和声波衰减特性的研究,可以揭示材 料的弹性系数和材料中压力、张力的关系。由脉冲发生器输出的重复脉冲是超声波的激励源 ,它可以用作取样积分的同步参考信号。 图531超声波检测材料特性 4. 荧光光谱测量 在荧光光谱测量中,通过测量荧光的衰减特性可以了解其发光机理。图532所示是利用双 通道取样积分器测量纳秒荧光光谱的结构框图。可调染料激光器光源经分光镜分成两路光束 A和B,光束A照射样品激发荧光,经单色仪用光电倍增管检测; 光束B用于对激光源强度进行 检测。因为测量时间长,光源强度会发生变化,利用双通道取样积分器作A/B运算,可以补 偿光源强度的变化,消除光源波动的影响。 图532双通道取样积分器测量荧光光谱 5.5数字式平均 随着集成电路技术和计算机技术的日益发展,数字处理方法得到越来越广泛的应用。与单点 取样积分器系统相比,多点数字式平均方法在每个信号周期内取样多次,信号利用率高,利 用数字式累加代替模拟电路积分,并利用数字式存储器存储处理结果,没有漏电和漂移问题 。这些特点使得数字式多点平均方法得到广泛应用,计算机的普及也为这种推广应用提供了条件 。 模拟式取样积分和数字式多点平均的特点对比示于表51。 表51取样积分和数字式多点平均方法的特点比较 方法SNIR取样效率频率保持时间取样脉冲宽度Tg 取样积分2Tc/TgTg/T×100%适于高频 差(电容)窄(分辨好) 数字式平均N(N为累加次数)100% 适于低频好(存储器)不太窄(分辨差) 从表51可见,取样积分和数字式平均各有特点,取样积分的取样门可以做得很窄,适用于 恢复高频信号; 数字式平均必须用A/D转换器将模拟信号转换为数字信号。因为模数转 换需要时间,取样脉冲宽度不容易做得很窄,所以更适用于低频和中频信号。但是,随 着微电子技术的发展,取样保持电路和A/D转换器芯片的工作速度越来越快,计算机的 运行速度也在不断提高。这样一来,数字式平均的取样脉冲宽度也可以做得很窄,也 能适应恢复高频信号的需要,所以数字式平均得到了越来越广泛的应用。 5.5.1数字式平均的原理及实现 在微弱信号检测领域中,常常遇到有限时段的信号重复出现许多次,相邻出现的信号之间的 时间间隔可能是固定值,也可能是变化量。这样的信号称之为重复信号。当各次重复信号之 间的时间间隔为固定值时,该信号就是周期信号。如果这样的信号被噪声污染,通过对多次 重复的信号进行数字式平均,就可以改善其信噪比。 数字式平均的工作过程为: 由采样保持器对被测信号进行取样,再由A/D转换器将被测 信号取样值变换为数字量,并将其存储在寄存器或存储器中; 累加平均的运算过程由微处理 单元(MPU)或数字信号处理器(DSP)完成,运算结果存储在寄存器或存储器中, 并可由D/A转换器输出相应的模拟量。 图533所示为周期信号的取样和数字式平均运算过程,被测信号的周期为T,在每个周期的 起始处触发取样过程,每个周期内均匀取样M次,取样时间间隔为Δt。对比图533与 图527可见,两者的工作过程非常类似,只不过实现方法不同。 图533数字式平均的取样和运算过程 对于第j道取样信号,数字式平均的运算过程可以表示为 A(j)=1N∑N-1i=0x(tj+iT),j=1,2,3,…,M(553) 式中,T为被测信号周期,tj是第一次对第j道信号取样的时刻。对于j=1,2,3,…,M,按 式(553)分别计算出各种j所对应的数字平均值A(j),并经过D/A转换器依次输出A(j ),就可得到平均后的被测信号波形。因为被测信号为确定性的信号,所以多次平均后仍然 为信号本身; 而干扰噪声为随机噪声,多次平均后其有效值会大为减少,从而提高信噪比。 对于周期不同的重复信号,运用数字式平均的关键是如何确定每段信号的起始点,换言 之,取样过程要与信号的出现保持同步。如果重复信号是由某个其他信号源激励产生的,即 使该激励信号是非周期或不规则的,那么也可以利用该激励信号作为每次重复开始取样的同 步信号。如果无法从其他信号源或先验知识确定每段信号的起始点,则只能由测量信号本身 来确定各次重复的起始点,例如利用被测信号的幅度或斜率来触发取样过程。对于噪声污染 严重 的信号,上述确定起始点的方法会遇到困难,必须采用比较复杂的检测和处理方法,这往往 需要一定的时间。为了在确定信号起点的过程中不丢失信号的前段数据,需要在信号通道中 设置延时环节或记录装置,这样才能在确定信号的起始点后(这需要一段时间的信号),再 从起始点处开始取样; 也可以在确定信号起点的过程中就对信号进行取样和存储,确定起始 点后,再确定这些取样数据的取舍及序号。 数字式平均器的功能结构框图示于图534。被测信号经预处理和取样保持(S/H),由 模数转换器(ADC)将其数字化后送入计算机进行累加,累加结果存入存储器。如果开 关K处于位置B,则无反馈,当下一次取样值被数字化后,与同一序号j的存储值相 累加,并将结果存入同一地址,这就完成了一步式(553)所表示的线性累加平均。如果 开关K处于位置A,上次的累加存储值还要经过DAC,变换成模拟量反馈到S/H ,计算机分别对信号值和反馈值进行取样,经运算后进行累加,这样可以实现下面将要介 绍的归一化平均或指数式平均。上述所有与时间相关的操作都由定时及控制电路来同步,以 保证各步操作协调一致。平均后的数字式信号波形可以经DAC作模拟量显示,也可以直 接进行数字量显示。 图534数字式平均器的功能结构框图 图534所表示的各种功能可以用普通计算机附加必要的接口插板来实现,也可以用单片机 附加相应的接口电路开发成专用的检测仪器,还可以开发成嵌入式微机系统与其他设备联合 使用。 5.5.2数字式平均的信噪改善比 设被测信号为 x(t)=s(t)+n(t)(554) 式中,s(t)为有用信号,n(t)为干扰噪声。设每个信号周期中的取样道数为j=1,2,…,M,取样间隔为Δt,重复次数为i=0,1,2,…,N-1(见图533)。第j道第i次的取样值 可以表示为 x(ti+jΔt)=s(ti+jΔt)+n(ti+jΔt)(555) 式中,ti是第i个取样周期中开始取样的时刻。因为s(t)是确定性信号,所以对于不同的 采样周期i,第j道的取样值基本相同,可以用sj来表示; 而噪声n(t)是随机的,其数值既 取决于i,又取决于j,所以式(555)可简记为 xij=sj+nij(556) 数字式累加平均的计算过程可以表示为 Aj=1N∑N-1i=0xij(557) 下面根据不同的情况,分析用式(557)所表示的数字式平均可以达到的信噪改善比SNIR。 1. 平均次数N→∞时的SNIR 当N→∞时,xij的平均结果为其数学期望值,由式(556)可得 E[xij]=E[sj]+E[nij] 对于确定性信号sj,其数学期望值E[sj]为信号本身sj; 对零均值噪声nij, 其数学期望值为零,即 E[nij]=0 所以 E[xij]=sj(558) 也就是说,当平均次数N→∞时,平均器输出信号的信噪比SNRo=sj/0=∞,达到 的信噪改善比SNIR也为无穷大。换言之,通过增加平均次数N可使SNIR达到任何希望的值。 2. 平均次数N有限时的SNIR 1) nij为高斯分布白噪声 设高斯分布零均值白噪声nij的有效值(均方根值)为σn,则对单次取样xij= sj+nij,其有用信号数值为sj,则平均处理之前的信噪比(有效值之比)为 SNRi=sj/σn(559) N次累加后的结果为 ∑N-1i=0xij=∑N-1i=0sj+∑N-1 i=0nij(560) 因为sj为确定性信号,N次累加后幅度会增加N倍。而噪声nij的幅度是随机的 ,累加的过程不会是简单的幅度相加,只能从其统计量的角度来考虑。取样累加后噪声的均 方值为 n2ij=E[n0j+n1j+…+n(N-1)j]2=E ∑N-1i=0n2ij+2E∑N-2i=0∑N-1m=i+1nijnmj(561) 式中右边的第1项表示噪声的各次取样值平方和的数学期望值,第2项表示噪声在不同时刻的 取样值两两相乘之和的数学期望值。 只要信号周期T足够大,则不同时刻的噪声取样值nij与nmj(i≠m) 互不相关, 其乘积的数学期望值为零,式(561)右边第2项为零,则 n2ij=E∑N-1i=0n2ij=Nσ2n(562) 累加后噪声的有效值为 σno=n2ij1/2=Nσ2n =Nσn(563) 累加后信号的电压值为 ∑N-1i=0sj=Nsj(564) 累加后输出信号的信噪比(有效值之比)为 SNRo=NsjNσn=Ns jσn(565) 由式(559)和式(565)可得 SNIR=SNRoSNRi=N(566) 式(566)说明,当污染噪声为白噪声时,N次不同时刻取样值的累积平均可以使信噪比 改善N倍,这就是常说的N法则。在取样积分的分析过程中, 已经多次使用过这一结论。 2) nij为高斯分布有色噪声 对于零均值有色噪声,式(561)中右边第2项所表示的噪声的不同时刻取样值两两相乘之 和的数学期望值不为零,设|i-m|=k,考虑到 E[nijn(i+k)j]=Rn(k)(567) 式中,Rn(k)表示噪声的自相关函数,Rn(0)=σ2n表示噪声的功率。取样累加后噪声的 均方值为 n2ij=E[n0j+n1j+…+n(N-1)j]2=E ∑N-1i=0n2ij+2E ∑N-2i=0∑N-1m=i+1nijnmj =NRn(0)+2∑N-1k=1(N-k)Rn(k)(568) 平均后噪声的均方值为 σ2no=n2ij/N2=[Rn(0)/N]1+2 N∑N-1k=1(N-k)ρn(k)(569) 式中 ρn(k)=Rn(k)/Rn(0)(570) 为噪声n(t)的归一化自相关函数。考虑到Rn(0)=σ2n,平均后输出信号的信噪比(有效值之比)为 SNRo=Sj/σno=sjσ2nN·1+2N∑N-1k=1(N-k)ρn(k)(571) 由式(559)和式(571)可得 SNIR=SNRoSNRi=N1+2N∑N-1k=1(N-k)ρn(k)(572) 当污染噪声n(t)为白噪声时,对于k≠0,其归一化自相关函数ρn(k)=0,式(572)就简化为式(566)。对比式(5 72)与式(566)可知,在污染噪声为有色噪声的情况下,数字式平均所达到的信噪改 善比要低于白噪声的情况。 5.5.3数字式平均的频域分析 1. 用傅里叶变换方法分析 被测信号x(t)的第i个取样x(t-iT)可以看成是x(t)与δ(t-iT)的卷积,即 x(t-iT)=x(t)*δ(t-iT)(573) 式中,T为取样周期。数字式平均的结果为 A(t)=1N∑N-1i=0x(t-iT)(574) 将式(573)代入式(574)得 A(t)=1N∑N-1i=0[x(t)*δ(t-iT)] =x(t)*1N∑N-1i=0δ(t-iT)(575) 如果将平均器看作一个线性系统,其输入信号是x(t),输出信号是平均结果A(t),系统 的冲激响应函数为h(t),频率响应函数为H(jω),如图535所示,则 A(t)=x(t)*h(t)(576) 图535平均器系统表示 对比式(575)与式(576)可知,平均器的冲激响应函数为 h(t)=1N∑N-1i=0δ(t-iT)(577) 对式(577)作傅里叶变换,得平均器的频率响应函数H(jω)为 H(jω)=∫∞-∞h(t)e-jωtdt =∫∞-∞1N∑N-1i=0δ(t-iT)e -jωtdt=1N∑N-1i=0e -jωT(578) 式(578)右边的级数是长度有限的几何级数,求和结果可以表示为 H(jω)=1N1-e-jNωT1-e-j ωT(579) 根据欧拉公式,得 H(jω)=1Nsin(NωT/2)sin(ωT/2)e-j(N-1)ωT/2(580) 相应的幅频响应特性为 |H(jω)|=1Nsin(NωT/2)sin(ω T/2)(581) 根据极限定理有 limx→0sinmxsinnx=mn 将上式关系应用于式(581)可知,当ωT=0或ωT为2π的整数倍时,|H(jω)| =1。 根据式(581),图536示出对于几种不同N的|H(jω)|图形。 由式(581)和图536所示图形可见,当平均次数N=1时,其传输特性为全通滤波器; 当N>1时,平均器 的 传输特性表现为梳状滤波器,每个齿是一个带通滤波器,各齿的中心频率位于fT为整数处, |H(jω)|的最大值为1。N值越大,各齿所表示的带通滤波器的带宽越窄。 图536平均器幅频响应 令| H(jω)|=1/2,可得平均器的-3dB带宽B为 B≈0.886NT(582) 可见,平均器的带宽反比于总测量时间NT。例如,若被测信号周期为T=1ms,平均次数 N=106,则总测量时间为NT=16.7min,带宽B=8.86×10-4Hz。所以, 时域平均也是频域的窄带化技术。 2. 用Z变换方法分析 用序列号i(i=0,1,2,…,N-1)代替式(574)中的t-iT,平均器计算过程为 A(t)=[x(i)+x(i-1)+x(i-2)+…+x(i-N+1)]/N(583) 利用Z变换方法,由式(583)可得平均器的离散传递函数为 H(z)=A(z)x(z)=1N(1+z-1+z-2+…+z -N+1)(584) 令z=ejωT,T为取样间隔,得平均器的稳态频率响应为 H(ejωT)=1N(1+e-jωT+e-j2ωT+…+e-j(N-1)ωT )(585) 对式(585)中的等比级数求和,可得 H(ejωT)=1N 1-e-jωNT1- e-jωT(586) 这与式(579)的结果是一致的。 5.5.4数字式平均算法 1. 线性累加平均算法 式(574)所表示的累加平均过程是一种线性累加平均过程,每个取样数据在累加中的权 重都一样。为简单计,用序号n来表示时刻t-iT,n=i=0,1,2,…,N-1,则式(574)可以 改写为 A(t)=1N∑N-1n=0x(n)(587) 这是一种批量算法,采集完N个数据后,再由计算机计算其平均值。这种算法的缺点是计算 量较大,需要做N次累加和一次除法才能得到一个平均结果,所以获得结果的频次较低。 根据前面的分析,对于噪声为白噪声的情况,平均过程能够实现的信噪改善比为 SNIR=N 2. 递推式平均算法 为了增加获得平均结果的频次,可以在每次取样数据到来时利用上次的平均结果做更新运 算 ,以获得新的平均结果。用A(N-1)表示时刻N-1的前N-1个数据的平均结果,A(N)表示时刻N 的平均结果,x(N)表示时刻N的取样值,式(587)可改写为 A(N)=1N∑Nn=1x(n)=N-1N 1N-1∑N-1n=1x(n)+1N x(N) =N-1NA(N-1)+x(N)N(588) 式(588)所表示的差分方程的运算过程示于图537。 图537递推式平均运算过程框图 利用这种递推式平均算法,当每个取样数据到来后,可以利用新数据对上次的平均结果进行 更新,这样相对于每个取样数据都会得到一个平均结果。随着一个个取样数据的到来,平 均 结果的信噪比越来越高,被测信号的波形逐渐清晰。针对几种不同的平均次数N,图538示 出这种逐渐清晰的被测波形。 由式(588)可得 A(N)=A(N-1)+x(N)-A(N-1)N(589) 可见,每次递推的过程都是对上次的运算结果附加一个修正量,修正量的大小取决于新的取 样数据与上次平均结果的差值以及平均次数N。随着时间的推移,平均次数N越来越大,式( 589)中右边的第2项所表示的修正量会越来越小,也就是新数据的作用越来越小。数字电路和 计 算机中的数据都有一定的字长和范围,当N大到一定程度后,该修正量会趋向于零,此后继 续取样和递推都不会对信噪比的改善起作用,平均结果稳定不变。如果被测信号波形发生 了变化,平均结果也不能跟踪这种变化,所以该算法不适于对时变信号进行处理。 图538递推式平均输出波形 3. 指数加权平均算法 1) 算法 如果把式(589)中的N用较大的固定正数α代替,可得 A(N)=α-1αA(N-1)+1αx(N)(590) 令β=(α-1)/α,若α>>1,则有0<β<1且接近于1,则式(590)可改写为 A(N)=βA(N-1)+(1-β)x(N)(591) 这就是指数加权数字式平均的算法,它也是在每次取样数据到来时,根据新数据对上次的平 均结果进行修正,得到本次的平均结果。参数β决定了递推更新过程中新数据和原平均结果 各起多大作用,算法的特性对β的依赖性很大。 相应的运算过程方框图示于图539。 2) 对取样数据的指数加权 将式(591)展开,得 A(N)=(1-β)∑Nn=1βN-nx(n)(592) 由式(592)可见,平均过程是把每个取样数据乘以一个指数函数,再进行累加,所以这是 一种指数加权平均,数据的序号n越大,权越重。因此,在平均结果中,新数据比老数据起 的作用要大,最新的数据(n=N)权重为1。图540示出n=N和n=N+1时的加权函数βN- n的图形。 图539指数加权数字式平均流程框图 图540加权函数βN-n的图形 3) 传递函数 对式(591)两边做Z变换,得 A(z)=βz-1A(z)+(1-β)x(z) 由此可得,指数加权数字式平均的离散传递函数为 H(z)=A(z)x(z)=1-β1-βz-1(593) 令式(593)中的z=ejωT,T为取样间隔,得指数加权平均器的稳态频率响应为 H(ejωT)=1-β1-βe-jωT(594) 其幅频响应为 |H(ejωT)|=1-β1+β2-2βcosωT (595) 根据式(595)画出几种不同β情况下的|H(ejωT)|图形,如图541所 示。可见,当β=0时,平均器为全通滤波器; 当0<β<1时,由x(k)到A(k)的传输过程为一 阶低通滤波器,其带宽取决于β,β越接近于1,带宽越窄。这种算法的特点是: 可以 用于对时变信号进行低通滤波处理,而且算法简单,易于实现。 图541指数加权平均器的幅频响应 4. 移动平均算法 数字式移动平均算法的工作原理类似于图514所示的扫描式移动取样积分,下面以常用的几种移动平均算法为例说明其原理。 1) 三点移动平均 三点移动平均算法为 A(n)=0.25[x(n-1)+2x(n)+x(n+1)](596) 对式(596)两边做Z变换,得 A(z)=0.25 x(z) [z-1+2+z](597) 由此可得,三点移动平均算法的传递函数为 H(z)=A(z)/ x(z)=0.25(z-1+2+z)(598) 令式(598)中的z=ejωT,得稳态频率响应为 H(ejωT)=0.25(2+e-jωT+ejωT)=0.25(2+2cosωT)=cos2(ωT/2)(599) 该函数在ωT=0~π区间(相应于f/fs=0~0.5,fs为采样频率)是单调减函数,因此具有低通特性,对采样数据具有平滑功能。式(596)所表示的算法又叫作海宁滤波器。 2) 五点移动平均 五点移动平均算法为 A(n)=[x(n-2)+x(n-1)+2x(n)+x(n+1)+x(n+2)]/6 其传递函数为 H(z)=(z-2+z-1+2+z+z2)/6 令z=ejωT,得稳态频率响应为 H(ejωT)=[cos(2ωT)+cos(ωT)+1]/3(5100) 3) 七点移动平均 七点移动平均算法为 A(n)=[ x(n-3)+x(n-2)+x(n-1)+2x(n)+ x(n+1)+x(n+2)+x(n+3)]/8 其传递函数为 H(z)=( z-3+z-2+z-1+2+z+z2+z3)/8 令z=ejωT,得稳态频率响应为 H(ejωT)=[cos(3ωT)+cos(2ωT)+cos(ωT)+1]/4(5101) 相应于式(599)、式(5100)和式(5101)的频率响应曲线示于图542。可见,移动平均的点数越多,通带的带宽越窄,对信号中高频分量的平滑作用越明显,但是计算工作量也越大。 图542三点、五点和七点移动平均的幅频响应 上述平均算法都具有低通特性。利用数字信号处理的不同算法,还可以实现高通、带通和带阻等滤波功能。 数字式运算的计算装置都有一定的字长,计算过程中数值的量化必然会有舍入误差,关于这种误差导致的算法频响偏差,以及由此产生的“死带”问题,读者可查阅《数字信号处理》的有关章节。 在实际应用中,数字式平均算法一般都是以微型计算机为中心,可以实现多种平均模式及 其他 数字信号处理功能; 而取样积分具有快速取样的优点,如果将其与微型计算机组合在一起, 就可以兼有两种方式的优点。在这种新型的时域平均系统中,一般都是把高速取样部分做成 几 种不同指标的组件,可以根据需要进行选择。例如,在某种已经产品化的以微型机为核心的 平均系统中,一种高速取样头的取样脉冲宽度为25ps~1ns,而另一种高速取样头的 取样脉冲宽度为2ns~2ms,使用者还可以对取样单元的参数进行设定。 利用时域平均方法提高信噪比,可以把有用信号从噪声中提取出来,取样积分还可以实现很 高 的分辨率,这种方法已经在很多领域中得到广泛的应用,例如用于荧光衰减测量、激光探测 实验 、激光诊断、深能级瞬变光谱、脉冲超声波检测裂缝或缺陷、脉冲核磁共振、噪声分析、拉曼 光谱测量、光纤特性测量、光吸收研究、电子自旋共振、脑电波测定、振动分析、电化学研 究、光电检测、地震信号与生物物理学信号采集与处理等。 习题 51指数式取样门积分器工作于定点模式,积分器时间常数TC=RC,取样门宽度为Tg。如果污染噪声为白噪声,要求功率信噪改善比SNIRP>100,试求Tg选择范围。 52污染信号s(t)的噪声n(t)为白噪声,经过带宽为B的滤波器后再累积平均。对于给定的 SNIR,须使平均重复次数≥N才能满足要求。试证明N∝B。 53指数加权平均算法如下: A(n)=0.96 A(n-1)+0.04x(n) 由于计算装置字长有限,运算结果都四舍五入为整数。设x(n)输入值都是100,A(n)初始值为(a)85,(b)115。找出这两种情况下A(n)输出的稳态值,以及这种输入数据条件下的“死带(dead band,数据进入此带不能再更新)”宽度。 54加权系数相同的12点移动平均算法为 A(n)=112[ x(n)+x(n-1)+…+x(n-11)] (1) 写出其递推公式。 (2) 画出其阶跃响应。 (3) 在0<f<1/(2T)(T为取样周期)频率范围内画出其幅频响应。 55数字滤波器传递函数为 H(z)=zz+0.9 图P56 在z平面画出其零点极点分布图,给出其输入和输出关系的递推公式,在-2π/T<ω<2π/T(T为取样周期)范围内画出其幅频响应曲线,说明滤波器功能。 56对于图P56所示变形的移动累加电路,y(n)=x(n)-2Ax(n-1)+x(n-2)。若A=1,试求其幅频响应函数|H(f)|。若A<1且可调,试求其幅频响应函数|H(f)|。讨论上述两种情况下电路实现的功能。