微弱信号检测第四章 相关检测 NEW

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
0
1 T
T
0 R xy ()dt R xy ()
由式知,尽管T有限,Rxy(τ)是Rxy(τ)的无偏估计。
微 弱信号检测
估计值的均方误差为:
varR~
xy
()
E(R~
xy
()
R
xy
())
2
对于高斯分布零均值限带白噪声x(t)和y(t),若其带宽为B,
则可以证明:
varR~ xy()
1 2BT
1 2
若ρxy(τ)=0.5,B=100HZ,要求ε<5%,则应使T>10S。
当信号带宽较窄时,需要较长的积分时间,这是相关 测量系统的主要缺点。
2.Rxy(τ)估计值的归一化均方根误差
varR~ xy () 1
1 2xy ()
R xy ()
2BT xy ()
微 弱信号检测
1
一般情况下ρxy(τ)<1/3,故
~ Rxy (k)
~
R xy (1)
R~ xy (M 1)
1 N
x(0)
x(1)
x(1 M)
x(1) x(0)
x(2 M)
x(N 1) y(0)
x(N
2)
y(1)
x(N M)y(N 1)
微 弱信号检测
两种计算方法:①所有数据采集完毕后计算;
②边采集边计算;
~
R
xy
⑦通用和专用相关仪的研发方面,1972年,用PMOS技术实 现溢出式极性峰点检测技术;此后专业仪表公司研制了多种 通用相关仪; ⑧1987年,Beck教授开发出实用的相关流速仪;
⑨1984年,VLSI相关仪问世;同年代英国的Kent公司开发 出相关检漏仪;
⑩目前研究还在深入,同时扩展到光学信号等领域;
T0
积分时间越长,估计值越接近真实值。
微 弱信号检测
2.数字累加方式
将被测信号x(t)和y(t)取样并模数转换,得离散的数字信 号x(n)和y(n),则可用累加平均的方法实现积分运算:
~
Rx (k)
1
N1
x(n)x(n k)
N n0
~
R xy (k)
1
N1
y(n)x(n k)
N n0
N表示累加平均的次数,k为延时序号
微 弱信号检测
④距离检测 ◆如某种对象的运动速度已知,那么测出它在两点之间的 渡越时间,就可以计算出这两点之间的距离。
⑤系统动态特性识别 ◆又叫做系统辩识。
⑥其他应用 ◆如气体色谱分析、光子相关分析、火焰燃烧情况检测、 天文想象观测和生物医学应用等。
三、相关检测发展
①1953年,贝尔实验室,磁带记录仪技术来实现相关检测;
4.数字相关量噪声导致的SNR退化
退化系数定义为:D
模拟相关的 数字相关的
SNR SNR
微 弱信号检测
§4-3相关函数算法及实现
随着技术发展,当今的相关检测设备多采用数字式运算。
~
即:R xy (k)
1
N1
y(n)x(n k)
,k=0,1,2,…,M-
N n0
1
矩阵表示为:
~
R xy (0)
文献资料证明,继电式相关函数与原相关函数之间的
关系为:
~
R xy ()
2 Rxy ()
Rx (0)
微 弱信号检测
2. 模拟积分继电式相关的实现方法:
y(t)
K
∫ -1
~ Rxy(τ)
x(t)
sgn[x(t)] τ
单级继电器式相关检测运算电路
f
移位寄存器 x(t)
sgn[x(t)]
1 …m … M 利用移位寄存器实现符号函数的延时
~
R xy (0)
(k)
~
R xy (1)
R~ xy (M 1)
1 N
x(0)
y(0)
x(1)
x(1 M)
1 N
x(1)
y(1)
x(0)
x(2 M)
x(N 1)
1
y(N 1)
x(N 2)
N
x(N M)
微 弱信号检测
一、递推算法
根据上次相关函数的计算结果,当下一个取样数据到来时, 对原有相关函数的计算结果进行更新,从而得到新的相关 憨数值。
若两路离散输入信号x(n)和y(n)它们的离散傅里叶变换 式分别为X(m)和Y(m),即:
N 1
X (m) x(n) exp( j2nm / N )
n1
N 1
Y (m) y(n) exp( j2nm / N )
则:
n1
~
DF [ R xy (k )]
DF [
1
N 1
y(n)x(n
k)]
④极性相关只能用于对时延和速度等和幅度无关的测量
微 弱信号检测
四、其他相关算法
1. 修正的极性相关算法:
通过叠加符合一定条件的伪随机信号,可以消除极性相关 的非性偏差。
x(t)
n1(t) y(t)
Ρ’xy(τ) 极性相关
n2(t) 修正的极性相关器原理
微 弱信号检测
若x(t)和y(t)为有界的随 机函数,n1(t)和n2(t)互相 独立、均匀分布且和分别对 x(t)和y(t)独立。满足 |x(t)|≤max|n1(t)|=A和 |y(t)|≤max|n2(t)|=A的条 件下,文献资料证明,得到 的修正极性相关函数:
1.
算法:
~
R xy ()
1
T
sgn[y(t)]sgn[x(t )]dt
T0
其中sgn[y(t)]和sgn[X(t-τ)]分别表示y(t)和x(t-τ)的符
号函数。
~
数字累加平均算法:R xy (k)
1
N1
sgn[y(n)]sgn[x(n k)]
N n0
微 弱信号检测
2. 电路实现 两符号函数乘积结果
②1961年,Weinreb提出了利用自相关法从随机噪声中提取 周期信号的理论; ③1966年,Van Fleck研究了用过零时刻相关法实现极性相 关运算;
微 弱信号检测
④1969年,HP公司的HP3721A相关仪问世,数字电路技术; ⑤1969年,英国Beck教授确立了通过用相关法检测自然流动 噪声渡越时间来测定流速的基本理论;互相关流速仪发展; ⑥1984年,Egau将极性相关应用于天文研究;
微 弱信号检测
f
移位寄存器
x(t)
sgn[x(t)]

y(t) 电子开关阵列
-1
∫∫


扫描多路开关 ~ Rxy(τ)
多级继电器式相关检测运算电路
微 弱信号检测
3. 数字累加平均实现积分平均运算:
~
R xy (k)
1
N1
y(n)sgn[x(n k)]
N n0
三、极性相关算法
在极性相关算法中,两路输入信号均被量化为1bit,即只 取其正负符号。
文献资料证明,极性相关函数与原相关函数之间的关
系为:
~
R xy
()
2
arc sin
R xy () Rx (0)R y (0)
2
arcsin[xy ()]
极性相关的特点:
①极性相关函数R’’xy(τ)是有偏估计,取值【-1,+1】;
② R’’xy(τ)中只有输入信号的符号信息,没有幅度信息;
③极性相关函数与归一化相关函数呈单调的反正弦关系;
X~(m)Y (m)
N n0
N
微 弱信号检测
从而:
~
Rxy(n)]
DF 1[X~(m)Y (m)]
1
N1 X~(m)Y (m) exp(
j2nm / N)
N n1
x(t) S/H
ADC
FFT 求共轭
相乘
Rxy(τ) FFT
y(t) S/H ADC FFT
基于FFT的相关运算过程
微 弱信号检测
微 弱信号检测
原理:对相关函数进行微分,能获得延时跟踪环的调整信号。 其微分值有正有负,但在相关函数的峰点处,它总为零, 两侧的符号相反。
Rxy(τ)
问题:计算工作量更大。
τ dRxy(τ)/dτ
解决方法:计理论分析证明,先对一路输入信号进行微分, 再将其与另一路信号进行相关处理,得到的就是相关函数的 微分。
①从噪声中提取信号 ◆确定信号的不同时刻取值一般都具较强的相关性;而对 于干扰噪声,因其随机性较强,不同时刻取值的相关性 差。
②渡越时间检测
◆如两路随机信号具有延时关系,利用相关函数在该延时 值处取得最大值的特性,则可以由互相关函数的峰值位 置测量出该延时值的大小。
③速度检测 ◆如两个测点的距离为确定值,检测出目标物通过这段距 离所需的时间,就测出了目标物的运动速度。
微 弱信号检测
基于这种原理的两种相关函数峰点跟踪的实现方案框图:
x(t) 延时线
y(t) d/dt


(a)
VCO
f=K/τ
x(t)
x(t-τ) 延时线 -

VCO
f=K/τ
y(t) d/dt
(b) 相关函数峰点跟踪的两种实现方案
R x (0)R y (0) R xy ())2
当Rxy(τ)≠0 时, Rxy(τ)估计值的归一化均方误差为:
2
varR~ xy ()
R
2 xy
()
1 2BT
1
1
2
xy
()
微 弱信号检测
ρxy(τ)为x(t)和y(t)的归一化相关函数:
xy()
R xy ()
R xx
(0)R yy(0)
微 弱信号检测
§4-2 相关函数的实际运算及误差分析
一、相关函数的实际运算
1.模拟积分方式
对于平稳的随机信号x(t)和y(t),其自相关函数和互相关
函数在实际积分运算时是在有限时间T内计算相关函数的估计
值,即:
~
1T
Rx () x(t)x(t )dt T0
~
Rxy ()
1
T
y(t)x(t )dt
x(t)
n1(t) y(t)
Ρ’xy(τ) 极性相关
n2(t) 修正的极性相关器原理
~
Rxy( )
1 A2
Rxy( )
Rx (0)Ry (0)
1 A2
xy( )
①极性相关函数与归一化相关函数为线性关系;
②人为引入随机噪声,为求精度,需更长的积分时间。
微 弱信号检测
2. 基于快速FFT的相关算法:
FFT是时域信号和频域信号相互转换的工具,是计算相关 函数的一种有效方法。
~N
~ N1
R xy (k) R xy (1 )x(N k)y(N) ,0 1
~N
N
R xy (k) (1 ) Nnx(N k)y(N)
指数加权算法特点: n 1
①算法具有一阶LPF特性,其带宽取决于β,β越接近于1, 带宽越窄;
②可以跟踪时变的Rxy(k),β越小,跟踪能力越强;
§4-4相关函数峰点跟踪
在相关检测的许多应用中,不要求相关函数的数值大 小,只求与其峰点位置的延时值,如相关法测速、超声测 距、雷达测距和泄漏点定位等。
在相关函数的计算中,为消除噪声的不利影响,需较 长的积分时间,故相关检测仪器的响应速度往往较慢。为 解决这问题常利用相关函数峰点跟踪系统。
相关函数峰点跟踪系统----一种闭环跟踪系统,不是通 过反复计算所有延时范围内的相关函数来得到其峰点所在 位置的延时,而是随着峰点位置的变化自动调整延时,。
xy () 2BT
3.Rxy(τ)估计值的信噪比
定义为 SNR
ER~ xy ()
varR~ xy ()
有 ER~ xy () R xy ()
得 SNR R xy () varR~ xy ()
微 弱信号检测
SNR 1 2BTxy () 1 2xy ()
一般情况下ρxy(τ)<1/3,故 SNR xy () 2BT
同或门真值表
sgn[y]
sgn[x] -1
+
sgn[x] sgn[y]量化值
量化值 0
1
-1
+ -11
0
பைடு நூலகம்
1
0
+
-11 +
1
0
1
1
x(t)
延 时= τ
1
~
x(t)
∫ Rxy(τ)

N个计数脉冲
时= τ
加减计 +/-数器
y(t) (a)模拟积分式
y(t) 清零
(b)数字累加式
微 弱信号检测
3. 估计值的偏差
③算法简单,容易实现。
微 弱信号检测
二、继电式相关算法
在继电式相关算法中,一路输入信号为模拟量形式,而另 一路输入信号被量化为1bit,即只取其正负符号。
1. 算法:
~
R xy ()
1
T
y(t)sgn[x(t )]dt
T0
其中:
sgn[x]
1, 1,
x0 x0
在实际应用中,在零点 附近设计一小的回差。
~ N
R xy (k)
1
N
x(n k)y(n)
N 1 n0
1
N1
x(n k)y(n)
1
x(N k)y(N)
N 1 n0
N 1
N
~ N1
R xy
1
x(N k)y(N)
N1 N1
微 弱信号检测
递推算法特点: ①随着取样数的增加,计算精度不断提高;
②N值越大,新数据作用越小。
用固定值β代替N/N+1,得指数加权递推算法:
3.实际相关器分类 —从构成原理和工作方式
①模拟式相关器
②数字式相关器 --极性相关器是其特例
微 弱信号检测
③混合式相关器 --继电器相关器是其特例
④修正的混合式相关器
二、运算误差分析
1.估计值的方差 --以互相关函数为例进行分析
对互相关函数估计值两边求数学期望得:
ER~
xy
()
1 T
T Ey(t)x(t )dt
微 弱信号检测
第四章 相关检测
§4-1 相关检测概述
相关检测技术是基于信号和噪声的统计特性进行检 测的,相关函数是两个时域信号相似性的一种度量。
一、相关检测技术与相敏检测器对比
x(t)
x(t)
LPF
Rxy(t) ∫
y(t) (a)相敏检测器结构图
y(t) τ
(b)互相关检测器结构图
微 弱信号检测
二、相关检测技术应用
相关文档
最新文档