微弱信号检测——取样积分器

1、第四局部：取样积分器背景：实际信号：不都是正弦或方波,通常为一复杂的宽带函数,表现为淹没在噪声中的微弱周期短脉冲.例如：脑电、心电信号；核磁共振信号.实际要求：恢复淹没在噪声中的微弱周期短脉冲信号波形,并显示记录.锁定放大器无能为力.功能实现：用取样积分器,对信号屡次取样并加以平均,抑制混于信号中的噪声,得到完整的波形恢复.模拟取样积分器T微计算机组成的数字式多点信号平均器.第一节取样积分器的根本原理、定点式取样积分器：信号通道,参考通道,门积分器积分器类型：1、指数积分器：上图定点式取样积分器使用的RC积分器为指数积分器.微分方程：Uj(t)=RC叫(t)dt+u0(t)；微分方程解：阶跃信

2、号Eu作用下,Uo6=E(1-e/RC),Tc=RC为时间常数.作用：对信号进行积分和平均；阶跃信号作用下,当t=3RC时,到达输入的95%;当t=5RC时,到达输入的99%o2、线性积分器：定点式取样积分器还可使用线性积分器.微分方程：-c业式dtUi(t)RnCUo(t)1tt微分万程解：阶跃信号Eu作用下(电谷初始电压为0),u0(t)=-Ui(t)dt=-ERC0RCTc=RC为时间常数.3、积分器选择：线性门积分器的输出受到运算放大器线性工作范围的限制,较适用于信号幅度较小的场合.如果信号幅度较大,为数不多的假设干次取样积分就有可能使运算放大器进入非线性区,导致测量误差.在这种情况下

3、只能使用指数门积分器.定点式取样积分器工作过程间歇式充电,一个周期内只在Tg内进行,因此一次无法充至该点的信号值.r(t)为参考信号,与被测信号s(t)同频,周期为T;Tg称为门宽,周期也为T0每隔周期T取样一次,积分器(RC)对取样信号进行积分和平均.缺点：延时to需手动调解,因此要恢复全部波形需逐步改变取样脉冲的延时量,麻烦而且很难得到一个完整的波形.s(t)+r(t),二、扫描式取样积分器(Boxcar积分器)1、扫描式取样积分器组成2、扫描式取样积分器工作过程f为非常缓慢的逐渐延时的脉冲,对波形取样,可理解为使被测信号上的每一点依次取出上百次以至上万次样值,等效为定点取样积分,提升输出

4、信噪比.输出波形与被测脉冲波形完全一致,观测噪声得到抑制,但波形时间轴变化大大放慢,总测量时间达几分钟.缓慢变化的信号只能用记录仪来记录三、取样积分器抑制噪声原理n1n4n11.1一1.u0=x(t0kT)=一qs(t0kT)n(t0kT)nkzenk=snk=s1nJ一对白噪严,由于不同时刻噪严值不相关：一£n(t0+kT)电0,故nkzo1nlU0二一"s(tokT)=s(to)nk=0取样积分分为两个过程:1、周期信号的取样：满足取样定理(在这里要注意对取样定理的理解)2、取样信号的积累平均：n次取样平均后,信号与噪声值分别为:1nl1n41Us=_gs(t0+kT)

5、=s(t°);Un=_jzn2(t0+kT)=-n1=-=nk-0n.k-0nnn输入信噪比：SNR=s城；输出信噪比：SNRo=；信噪改善比:SNIR=剑风=4二n二n/、nSNR结论：增加积累次数,可更充分抑制噪声,提升信噪比,检测更微弱信号；积累次数增加,测量时间nT也同时增加；时间换取精度.第二节取样积分器的频响特性一、频响特性：对输入信号的一次取样,相当于：y(t)-x(t)-(t)=-j-x(t).(t-t)dtJ_oO对信号的n次取样累加,相当于：一y(t)=x、(t)x(t-T)、.(t)x(t-(n-1JT)、(t)二i-x(t)<(t-t)：(t-t-T)-

6、,It-t-(n-1)Tbdt(由于：x(t_T)、.(tt)dt=i-x(t)；(t-t-T)dt)_jQOJ_jqQ根据冲激响应的定义：y(t)=x(t).h(t)=：x(t)h(tt)df得取样积分器的冲激响应为：h(t)=、：(t)、：(t-T)、.It-(n-1)T1,幅频响应特性为: 取样积分器相当是一个梳状滤波器； 被测周期信号s(t)的谐波分量nfi可以无失真地通过； 观察噪声具有很宽的功率谱,除了nfi附近的功率谱分量可通过外,绝大多数噪声谱分量均被梳状滤波器抑制掉； n越大,梳状滤波器越理想,信噪比提升就越大.、梳状滤波器带宽1、梳状滤波器每个梳齿的带宽-3dB带宽：sin

7、(n2T)sin(T)1=、,可求出if=0.886/nT2例：设信号周期T=10ms,即频率f=100Hz,n=106,nT=104秒,相当Af=8.66父10Hz结论：取样积分器同样具有与锁定放大器类似的窄带滤波功能,从而可以充分排除噪声来提取信号.2、等效噪声带宽1由等效噪声带宽定义：吊=H0广H2codf,由于H侬是以取样脉冲基涉及各次谐波频率为中央的梳状滤波器,每周期内的滤波特性是相同的,因此先求01/T或0fi内的f:11T2fn12"2(.)df=n1H2(0)0二Tn2二,sinnx、211()dx=0sinxnT设信号的最高次谐波频率为fm,fm内包含的谐波数为fm

8、/f1,那么0fm的频带内的等效输出噪声带宽为：?为=上f1n3、信噪改善比结论：随着累加次数SNIR=Vnfnn的增加,等效噪声带宽减小,同时信噪改善比以新方式增加.第三节取样积分器的参数选择方法一、累加次数n：根据信噪改善比确定SNIR=M例：要求SNIR=100,那么n=104.对于极微弱信号检测时,n值要求更大.、取样脉冲宽度TgytAto_TLTg0.42*A对正弦信号取样,取样后的输出电压为：ut=UmSincot,to-Tg/2<t<to+Tg/2经RC积分后输出为：Uo=1°t2Umsincotdt=2Umsinsincetoto4.2当频率很低时,8Tg

9、To,贝JUo之巴口,叫sinto2TT当频率较图时,sin:二一,故输出电压下降,从而引起信号中高频分量的损失,损22失程度可用比值A表示:A_UoLsin®Tg/2"Uoo一工/2设信号的高端截止频率为fc,那么对该频率分量sin2二fJg/21/曰卬>,得至1:2二fcTg/22fcTg<o.42,取样门宽Tg满足TgM*2恢复总时间为:故取样门宽脉冲TgM呸gfc结论：Tg要根据被恢复信号的最高频率fc来确定；fc越高,那么Tg越小.取样门脉冲宽度Tg不能选择过宽,否那么会造成信号高频成分损失平滑作用,导致恢复的信号失真.三、积分器时间常数Tcn次积累时

10、,总的充电时间为nTg0由于是间歇式充电,电容上电压为阶梯状上升,而不是指数式上升.根据RC充电的指数形式,其有效的充电时间为5RC到达稳定值,故定点积分器积分器最多可充电次数为：nTg=5Tc=n=5Tc/Tg. TTg已确定后,信噪改善比为：SNIR3=向/Tg 已要求SNIR情况下,积分器的时间常数可按下式确定：Tc=nTg/5g 对某一时刻信号值的恢复时间为：Ts=nT=5TJ/Tg注意：对于指数门积分器,电容充电具有指数特性,充电接近饱和时即总积分时间接近5倍积分器时间常数,信号积累速度减慢,信噪比改良较少.针对白噪声情况,在2倍积分时问常数内,信噪比就改善的比拟明显,此时取Tc=n

11、Tg/2=SNIR2Tg/2四、波形恢复总时间Tss慢扫描时间,测量区域Tb内取样点数为Tb/Tg；单点恢复时间为Ts,波形五、取样积分器参数选择过程计算门脉冲宽度Tg_0.42/fc估计信号的最高频率分量fc第四节基线取样与双通道取样积分器一、基线取样基线漂移输出信号波形1、背景：取样积分器对信号按一周期移动取样,并屡次平均得到信号的恢复.因此一个信号需要长时间记录才能恢复,这就对系统的稳定性提出了苛刻的要求,否那么会引起基线漂移.基线不漂移输出信号波形 必须解决基线漂移问题 解决方法一：采用低漂移器件、老化挑选,改善系统漂移特性 解决方法二：采用基线取样方法改善漂移2、原理：一周期内先取样

12、一次信号的有效成份,再取样一次信号的基线,两者相减,即得到信号成份扣除漂移的有效幅值.基于假定一周期内基线漂移不变I参考信号取样脉冲AtB二、双通道取样积分器双通道取样积分器系统由两路基线取样积分器组成,利用两路取样积分器分别对被测信号和作为标准的另一路信号进行消除基线的取样积分,然后将被测量与标准量比照,以消除由信号基线以及鼓励源起伏造成的误差,使测量更为准确和可靠.例如：在激光测量中,当测量时间较长时,激光光源的强度可能发生变化,从而导致测量误差,在这种情况下,可利用A通道作根本测量,利用B通道检测光源的强度,A通道和B通道的输出作除法运算,从而补偿光源强度变化带来的影响.在运算功能上,除

13、了常用的A/B运算外,还可以采用A-B、A-B、lgA/B等运算功能.如果在B通道中附加延时环节,还可实现相关运算功能.三、快速信号取样对快速信号取样积分,关键要求小到几十ps的时间取样常用门积分器不易实现,必须将取样与保持分开实现累积所需要 取样电容C的数值不大,否那么会影响开门及关门时间 取样期间Ci上建立的电压用来继续对C2的积分 通过反应构成闭环,使系统具有工作范围扩大、线性度好、漂移小及稳定度高的特点第五节多点信号平均器 取样积分器优点：极高的分辨力和对快速脉冲信号的处理水平 取样积分器缺点：1、对低重复频率信号的恢复需要极长的时间.例SNIR=100累加1002次,10Hz波形,T

14、b=T,门宽Tg=100%,恢复时间Tss=11.6天；2、积分电容由于漏电而使保持时间有限；3、因零漂而使恢复的信号失真一、多点信号平均器：实时取样系统,在信号的一个周期内要取样多点,并一一对应地存贮在相应的存贮器内并进行累加和平均.等效于大量单点取样积分器在不同延时的情况下并联使用. 顺序闭合开关,对波形上各点顺序取样,在电容上进行累积; 累积完成后,再依顺序输出电容上的电压； 测量时间减少为原来的1/N,N为每一轮回取样点数； 实际上不应用电容组.平均次数1三、数字多点信号平均器：计算机技术,存贮器为数据存贮器1、原理128204832768工作方式：线性累加平均,归一化平均,指数平均开

15、关拨向1:无反应,完成线性累加平均开关拨向2:有反应,上一次存储值反应到S/H,分别对信号值和反应值取样,实现归一化平均或指数式平均2、用取样积分器提取信号举例X-Y记录仪yx3、用数字示波器完成多点信号平均器功能第六节噪声信号数据处理方法背景： 实验数据包含噪声,测量结果偏离真值,表现在曲线上为高频波动及毛刺； 经数据处理,使结果尽量趋近于被测物理量的期望值； 为获得其它参数,也可利用数据处理.一、递推式平均算法：增加获得平均结果的频次,每次取样数据到来时,利用上一次平均结果作更新运算,获得新的平均结果.A(N-1):时刻N-1前的N-1个数据平均结果1NN111N-11、A(N)=Xx(n

16、)-,一Zx(n)+x(N)-A(N-1)+x(N)(差分万程)NnJNN-1nJNNN随着新数据的逐渐到来,平均结果的信噪比越来越高,被测信号波形逐渐清楚 递推：A(N)=A(N-1).x(N)-A(N1)N 每次递推过程是对上一次的运算结果附加一个修正量 不适合对时变信号进行处理,时变信号可采用指数加权平均算法.、曲线的移动平滑(平均光滑)方法1、过程：平均T移动一点T平均T整个曲线数学表达式:y(n)=12M1MZk-JMx(n-k)2、一种移动平滑滤波器举例1二点平均：y(n)=-x(n-1)+2x(n)+x(n+1)41两边取Z变换：y(z)=z,x(z)+2x(z)+zx(z)41

17、传递函数：H(z)z2z411频响特性：H(e)=-e2e'(cos1)421五点平均：y(n)=-x(n-2)+x(n1)+2x(n)+x(n+1)+x(n+2)6j.1频响特性：H(ej)=-(cos2cos1)1七点干均：y二产逸+近0俨皿1)“2)"1频响特性：H(e)(cos3cos2cos1)43、频率响应特性曲线H(ej)4、总结：相当于数字低通滤波器,平均点数增多,频带越窄,对高频分量的平滑越明显三、多项式最小二乘平滑通过实验数据点,需要画出一条最正确拟合曲线,常用曲线拟合准那么为多项式最小二乘曲线存在波动和高频噪声“毛刺,波动误差都来自纵坐标,横坐标误差可忽

18、略不计1、原理五点数据移动多项式进行最小二乘平滑过程示意：首先最左端选五点测量数据,组成数据组,对此数据组中央点进行拟合；然后自左至右,在左边去掉一个测量数据,同时在右边增加一个测量数据,组成新数据组,再进行数据组中央点拟合；重复下去.每次拟合得到中心点的平滑数据,准确度较高.举例说明：常用二次和三次多项式对实验数据点进行最小二乘平滑,以二次多项式、5点数据平滑为例说明拟合过程测量数据值为x(n),n=k-2,k-1,k,k+1,k+2,x(k)对应x(n)中央点,即待拟合数据点拟合计算值为y(i)=a20+a2ii+a22i2,i=2,1,0,+1,+2,i=0对应y(i)中央点2k22k22拟合方均误差；(a20,a21,a20)-ly-x(n)1'(a20a21ia22i2)-x(n)i-2n=k-2i-2n-k-2-取偏导并令其为0,即-=0,-=0,-=0,止匕时君(a20,a21,a20)最小,实现最小二202122乘平滑.22kd25a20+0+10a22=£Zx(n)i=2一/2M20+10a2i+0=£Zix(n)i=2n*/2k也,一一一一一_.2.10a20+0+34a22=££ix(n)-n*N解方程得：a20=l-3x(k-2)12x(k-1)17x(k)12x(k1)-3x(k2)1351