第四章 相关检测.

1、微 弱 信 号 检 测第四章第四章 相关检测相关检测 相关检测技术是基于信号和噪声的统计特性进行检相关检测技术是基于信号和噪声的统计特性进行检测的，相关函数是两个时域信号相似性的一种度量。测的，相关函数是两个时域信号相似性的一种度量。4-1 4-1 相关检测概述相关检测概述一、相关检测技术与相敏检测器对比一、相关检测技术与相敏检测器对比x(t)y(t)LPF(a)相敏检测器结构图相敏检测器结构图x(t)y(t)(b)互相关检测器结构图互相关检测器结构图Rxy(t)微 弱 信 号 检 测二、相关检测技术应用二、相关检测技术应用从噪声中提取信号从噪声中提取信号确定信号的不同时刻取值一般都具较强的相

2、关性；而确定信号的不同时刻取值一般都具较强的相关性；而对于干扰噪声，因其随机性较强，不同时刻取值的相关对于干扰噪声，因其随机性较强，不同时刻取值的相关性差。性差。渡越时间检测渡越时间检测如两路随机信号具有延时关系，利用相关函数在该延如两路随机信号具有延时关系，利用相关函数在该延时值处取得最大值的特性，则可以由互相关函数的峰值时值处取得最大值的特性，则可以由互相关函数的峰值位置测量出该延时值的大小。位置测量出该延时值的大小。速度检测速度检测如两个测点的距离为确定值，检测出目标物通过这段如两个测点的距离为确定值，检测出目标物通过这段距离所需的时间，就测出了目标物的运动速度。距离所需的时间，就测出了

3、目标物的运动速度。微 弱 信 号 检 测距离检测距离检测如某种对象的运动速度已知，那么测出它在两点之间如某种对象的运动速度已知，那么测出它在两点之间的渡越时间，就可以计算出这两点之间的距离。的渡越时间，就可以计算出这两点之间的距离。系统动态特性识别系统动态特性识别又叫做系统辩识。又叫做系统辩识。其他应用其他应用如气体色谱分析、光子相关分析、火焰燃烧情况检测、如气体色谱分析、光子相关分析、火焰燃烧情况检测、天文想象观测和生物医学应用等。天文想象观测和生物医学应用等。三、相关检测发展三、相关检测发展19531953年，贝尔实验室，磁带记录仪技术来实现相关检测；年，贝尔实验室，磁带记录仪技术来实现相

4、关检测；19611961年，年，WeinrebWeinreb提出了利用自相关法从随机噪声中提取周提出了利用自相关法从随机噪声中提取周期信号的理论；期信号的理论；19661966年，年，Van FleckVan Fleck研究了用过零时刻相关法实现极性相关研究了用过零时刻相关法实现极性相关运算；运算；微 弱 信 号 检 测19691969年，英国年，英国BeckBeck教授确立了通过用相关法检测自然流动教授确立了通过用相关法检测自然流动噪声渡越时间来测定流速的基本理论；互相关流速仪发展；噪声渡越时间来测定流速的基本理论；互相关流速仪发展；19691969年，年，HPHP公司的公司的HP3721A

5、HP3721A相关仪问世，数字电路技术；相关仪问世，数字电路技术；19841984年，年，EgauEgau将极性相关应用于天文研究；将极性相关应用于天文研究；通用和专用相关仪的研发方面，通用和专用相关仪的研发方面，19721972年，用年，用PMOSPMOS技术实现技术实现溢出式极性峰点检测技术；此后专业仪表公司研制了多种通溢出式极性峰点检测技术；此后专业仪表公司研制了多种通用相关仪；用相关仪；19871987年，年，BeckBeck教授开发出实用的相关流速仪；教授开发出实用的相关流速仪；19841984年，年，VLSIVLSI相关仪问世；同年代英国的相关仪问世；同年代英国的KentKent公

6、司开发出公司开发出相关检漏仪；相关检漏仪；目前研究还在深入，同时扩展到光学信号等领域；目前研究还在深入，同时扩展到光学信号等领域；微 弱 信 号 检 测4-2 4-2 相关函数的实际运算及误差分析相关函数的实际运算及误差分析一、相关函数的实际运算一、相关函数的实际运算1.1.模拟积分方式模拟积分方式 对于平稳的随机信号对于平稳的随机信号x(tx(t) )和和y(ty(t) )，其自相关函数和互相，其自相关函数和互相关函数在实际积分运算时是在有限时间关函数在实际积分运算时是在有限时间T T内计算相关函数的估内计算相关函数的估计值，即：计值，即：dt)t (x) t (xT1)(RT0 xdt)t

7、 (x) t ( yT1)(RT0 xy积分时间越长，估计值越接近真实值。积分时间越长，估计值越接近真实值。微 弱 信 号 检 测2.2.数字累加方式数字累加方式 将被测信号将被测信号x(tx(t) )和和y(ty(t) )取样并模数转换，得离散的数字取样并模数转换，得离散的数字信号信号x(nx(n) )和和y(ny(n) )，则可用累加平均的方法实现积分运算：，则可用累加平均的方法实现积分运算：1N0nx)kn(x)n(xN1)k(R1N0nxy)kn(x)n( yN1)k(RN表示累加平均的次数，表示累加平均的次数，k为延时序号为延时序号3.3.实际相关器分类实际相关器分类 从构成原理和工

8、作方式从构成原理和工作方式模拟式相关器模拟式相关器数字式相关器数字式相关器-极性相关器是其特例极性相关器是其特例微 弱 信 号 检 测-继电器相关器是其特例继电器相关器是其特例混合式相关器混合式相关器修正的混合式相关器修正的混合式相关器二、运算误差分析二、运算误差分析1.1.估计值的方差估计值的方差 -以互相关函数为例进行分析以互相关函数为例进行分析对互相关函数估计值两边求数学期望得对互相关函数估计值两边求数学期望得: :dt)t (x) t ( yET1)(RET0 xy)(Rdt)(RT1xyT0 xy由式知，尽管由式知，尽管T有限，有限，Rxy() )是是R Rxyxy() )的无偏估计

9、。的无偏估计。微 弱 信 号 检 测估计值的均方误差为估计值的均方误差为: :2xyxyxy)(R)(R(E)(Rvar对于高斯分布零均值限带白噪声对于高斯分布零均值限带白噪声x(tx(t) )和和y(ty(t) )，若其带宽为，若其带宽为B B，则可以证明则可以证明: :2xyyxxy)(R) 0 (R) 0 (RBT21)(Rvar当当R Rxyxy()0)0 时，时， R Rxyxy() )估计值的归一化均方误差为估计值的归一化均方误差为: : )(11BT21)(R)(Rvar22xyxy2xy微 弱 信 号 检 测21) 0 (R) 0 (R)(R)(yyxxxyxyxyxy() )

10、为为x(tx(t) )和和y(t)y(t)的归一化相关函数的归一化相关函数: : 若若xyxy()=0.5)=0.5，B=100HZB=100HZ，要求，要求5%5%，则应使，则应使T10ST10S。 当信号带宽较窄时，需要较长的积分时间，这是相关当信号带宽较窄时，需要较长的积分时间，这是相关测量系统的主要缺点。测量系统的主要缺点。2.R2.Rxyxy() )估计值的归一化均方根误差估计值的归一化均方根误差)()(1BT21)(R)(Rvarxy2xyxyxy微 弱 信 号 检 测一般情况下一般情况下xyxy() )1/31/3，故，故 3.R3.Rxyxy() )估计值的信噪比估计值的信噪比

11、BT2)(1xy)(Rvar)(RESNRxyxy定义为)(R)(RExyxy有)(Rvar)(RSNRxyxy得微 弱 信 号 检 测)(1)(BT21SNR2xyxy一般情况下一般情况下xyxy() )1/31/3，故，故 BT2)(SNRxy4.4.数字相关量噪声导致的数字相关量噪声导致的SNRSNR退化退化退化系数定义为：退化系数定义为： SNRSNRD数字相关的模拟相关的微 弱 信 号 检 测4-34-3相关函数算法及实现相关函数算法及实现随着技术发展，当今的相关检测设备多采用数字式运算。随着技术发展，当今的相关检测设备多采用数字式运算。1N0nxy)kn(x)n( yN1)k(R，

12、k=0k=0，1 1，2 2，M-1M-1矩阵表示为：矩阵表示为：即：即：) 1N( y) 1 ( y) 0 ( y)MN( x)M2 ( x)M1 ( x) 2N( x) 0 ( x) 1( x) 1N( x) 1 ( x) 0 ( xN1) 1M(R) 1 (R) 0 (R) k (Rxyxyxyxy微 弱 信 号 检 测所有数据采集完毕后计算；所有数据采集完毕后计算；两种计算方法：两种计算方法：)M2 ( x) 0 ( x) 1 ( x) 1 ( yN1)M1 ( x) 1( x) 0 ( x) 0 ( yN1) 1M(R) 1 (R) 0 (R) k (Rxyxyxyxy边采集边计算；

13、边采集边计算；)MN( x) 2N( x) 1N( x) 1N( yN1微 弱 信 号 检 测一、递推算法一、递推算法根据上次相关函数的计算结果，当下一个取样数据到来时，根据上次相关函数的计算结果，当下一个取样数据到来时，对原有相关函数的计算结果进行更新，从而得到新的相关对原有相关函数的计算结果进行更新，从而得到新的相关憨数值。憨数值。)n( y )kn(x1N1)k(RN0nNxy)N( y)kN(x1N1)n( y )kn(x1N11N0n)N( y)kN(x1N1R1NN1Nxy微 弱 信 号 检 测随着取样数的增加，计算精度不断提高；随着取样数的增加，计算精度不断提高；递推算法特点：递

14、推算法特点：N N值越大，新数据作用越小。值越大，新数据作用越小。用固定值用固定值代替代替N/N+1N/N+1，得指数加权递推算法：，得指数加权递推算法：)N( y)kN(x)1 (R)k(R1NxyNxy可以跟踪时变的可以跟踪时变的R Rxyxy(k(k) )，越小，跟踪能力越强；越小，跟踪能力越强；10 , 算法具有一阶算法具有一阶LPFLPF特性，其带宽取决于特性，其带宽取决于，越接近越接近于于1 1，带宽越窄；，带宽越窄；指数加权算法特点：指数加权算法特点：算法简单，容易实现。算法简单，容易实现。N1nnNNxy)N( y)kN(x)1 ()k(R微 弱 信 号 检 测二、继电式相关算

15、法二、继电式相关算法在继电式相关算法中，一路输入信号为模拟量形式，而另在继电式相关算法中，一路输入信号为模拟量形式，而另一路输入信号被量化为一路输入信号被量化为1bit，即只取其正负符号。，即只取其正负符号。1. 1. 算法：算法：dt)t (xsgn) t ( yT1)(RT0 xy0 x, 10 x, 1xsgn其中：其中：在实际应用中，在零点在实际应用中，在零点附近设计一小的回差。附近设计一小的回差。 文献资料证明，继电式相关函数与原相关函数之间文献资料证明，继电式相关函数与原相关函数之间的关系为：的关系为：) 0 (R)(R2)(Rxxyxy微 弱 信 号 检 测. . 模拟积分继电式

16、相关的实现方法：模拟积分继电式相关的实现方法：单级继电器式相关检测运算电路单级继电器式相关检测运算电路x(t)y(t)-1Rxy()sgnx(t)K利用移位寄存器实现符号函数的延时利用移位寄存器实现符号函数的延时x(t)sgnx(t)移位寄存器移位寄存器f1mM微 弱 信 号 检 测多级继电器式相关检测运算电路多级继电器式相关检测运算电路y(t)-1Rxy()x(t)sgnx(t)移位寄存器移位寄存器f电子开关阵列电子开关阵列扫描多路开关扫描多路开关微 弱 信 号 检 测三、极性相关算法三、极性相关算法在极性相关算法中，两路输入信号均被量化为在极性相关算法中，两路输入信号均被量化为1bit，即

17、只，即只取其正负符号。取其正负符号。1. 1. 算法：算法：dt)t (xsgn)t ( ysgnT1)(RT0 xy 其中其中sgny(tsgny(t)和和sgnX(tsgnX(t-)分别表示分别表示y(ty(t) )和和x(tx(t-) )的符号函数。的符号函数。3. 3. 数字累加平均实现积分平均运算：数字累加平均实现积分平均运算： )kn(xsgn)n( yN1)k(R1N0nxy )kn(xsgn)n( ysgnN1)k(R1N0nxy数字累加平均算法：数字累加平均算法：微 弱 信 号 检 测. . 电路实现电路实现sgnysgny-1-1-1-1-1-1sgnxsgnx-1-1+1+1+1+1+1+1+1+1sgnxsgnx量化值量化值sgnysgny量化值量化值0 00 00 00