1、信号检测与估计理论信号检测与估计理论章信号波形的检测章信号波形的检测引言引言l研究目的:研究目的:l根据性能指标要求,设计与环境相匹配的接收根据性能指标要求,设计与环境相匹配的接收机(检测系统),机(检测系统),l从噪声污染的接收信号中提取有用的信号;从噪声污染的接收信号中提取有用的信号;l在噪声干扰背景中区别不同特性、不同参量的在噪声干扰背景中区别不同特性、不同参量的信号。信号。l理论基础:理论基础:l假设检验和似然比检验假设检验和似然比检验l最佳检测的判决方式(判决表示式)最佳检测的判决方式(判决表示式)l检测系统的结构检测系统的结构l检测性能分析检测性能分析l最佳波形设计最佳波形设计引言
2、引言l简化的信号流程模型简化的信号流程模型二元数字通信系统波形检测模型二元数字通信系统波形检测模型把随机过程把随机过程 用正交级数表示,进行统计描述,再应用信号的用正交级数表示,进行统计描述,再应用信号的统计检测理论来处理信号波形的检测问题。统计检测理论来处理信号波形的检测问题。From Steven M. Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入二元信号的模型:二元信号的模型:01: 0,1,.,1: 0,1,.,1Hx nw nnNHx ns nw nnN10;( );p x HL xp x H若似然比超过门限,即则检测器判断为H1成立。From Steven M.
3、 Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入检测信号似然函数:01.1xxx Nx121202211;exp 22NNnp x Hx ns n120202211;exp 22NNnp x Hx n 11222001exp 2NNnnLx ns nx nxFrom Steven M. Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入两边取对数,并转换可以根据不同的准则,确定门限值。 11222001ln ln2NNnnlLx ns nx n xx112220011 ln2NNnnx n s ns n 1122001 ln 2NNnnTx n s ns nx
4、10 NnTx n s nxFrom Steven M. Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入上式表明:检验统计量根据信号的值对数据样本进行加权。上式表明:检验统计量根据信号的值对数据样本进行加权。大的信号样本采用大较大的加权。大的信号样本采用大较大的加权。把接收到的数据和信号的仿形品进行相关运算。把接收到的数据和信号的仿形品进行相关运算。检测器称为相关器或仿形检测器称为相关器或仿形-相关器。相关器。 10 NnTx n s nxFrom Steven M. Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入用FIR(Finite Impulse Re
5、sponse)滤波器模型:0 nky nh nk x k =10,1.,-1h ns Nnn=N 0 x nN -1假定信号只在,上是非零的。令冲击响应为信号的镜像, 0,1.,h nn=N -1冲击响应当是非零的。0 1 nky ns Nnkx k 1011 NknNy Ns k x k时刻,From Steven M. Kayl检测器的建立检测器的建立匹配滤波器的引入匹配滤波器的引入用FIR滤波器模型来表达:将信号相对将信号相对n=0反转,再向右移反转,再向右移N-1个样本得到匹配滤波器的冲击响应。个样本得到匹配滤波器的冲击响应。From Steven M. Kayl检测器的建立检测器的建
6、立匹配滤波器的引入匹配滤波器的引入WGN中已知的确定性信号的检测问题,利用中已知的确定性信号的检测问题,利用NP准则与最大准则与最大SNR准则都准则都可以导出匹配滤波器。可以导出匹配滤波器。对于非高斯噪声,匹配滤波器输出的信噪比最大。对于非高斯噪声,匹配滤波器输出的信噪比最大。匹配滤波器的概念匹配滤波器的概念l通信、雷达等电子信息系统的接收机模型通信、雷达等电子信息系统的接收机模型对信号进行加工、处理,使信噪比最大检测判断(例如与门限值进行比较判断)若线性时不变滤波器输入的信号是确知信号,噪声是加性平稳噪声,若线性时不变滤波器输入的信号是确知信号,噪声是加性平稳噪声,则在输入功率信噪比一定的条
7、件下,使输出功率信噪比最大的滤波器,则在输入功率信噪比一定的条件下,使输出功率信噪比最大的滤波器,是一个与输入信号相匹配的最佳滤波器,即匹配滤波器(是一个与输入信号相匹配的最佳滤波器,即匹配滤波器(MF)。)。匹配滤波器的设计匹配滤波器的设计图4.3 线性滤波器匹配滤波器的设计匹配滤波器的设计l输出噪声的功率谱密度输出噪声的功率谱密度设滤波器输出信号设滤波器输出信号 在在 时刻出现峰值,有时刻出现峰值,有匹配滤波器的设计匹配滤波器的设计l输出功率信噪比输出功率信噪比利用利用Schwarz不等式,满足式(不等式,满足式(4.2.12) , 等号成立。等号成立。匹配滤波器的设计匹配滤波器的设计令令
8、由由有有当当式(式(4.2.16)中的等号成立。)中的等号成立。匹配滤波器的设计匹配滤波器的设计噪声为有色噪声时,广义滤波器:噪声为有色噪声时,广义滤波器:当滤波器输入为白噪声时,当滤波器输入为白噪声时, ,有有匹配滤波器的主要特点匹配滤波器的主要特点1. 匹配滤波器的脉冲响应与匹配滤波器的脉冲响应与 时刻的选择时刻的选择图4.4 匹配滤波器的脉冲响应特性在在 时刻,输入信号时刻,输入信号 s(t) 已全部送入滤波器,已全部送入滤波器,因此,因此, 至少要选择在输入信号至少要选择在输入信号 s(t) 的末尾。的末尾。匹配滤波器的主要特点匹配滤波器的主要特点2. 匹配滤波器的输出功率信噪比匹配滤
9、波器的输出功率信噪比若输入信号若输入信号 s(t) 的能量为的能量为 ,白噪声白噪声 n(t) 的功率谱密度为的功率谱密度为 ,匹配滤波器的输出功率信噪比为匹配滤波器的输出功率信噪比为与输入信号与输入信号 s(t) 的能量的能量 有关,与波形无关。有关,与波形无关。 22102sEst dtSdR匹配滤波器的主要特点匹配滤波器的主要特点3. 匹配滤波器的适应性匹配滤波器的适应性对振幅和时延参量具有适应性;对振幅和时延参量具有适应性;对频移信号不具有适应性。对频移信号不具有适应性。匹配滤波器的主要特点匹配滤波器的主要特点4. 匹配滤波器与相关器的关系匹配滤波器与相关器的关系对平稳输入信号对平稳输
10、入信号 ,自相关器的输出为,自相关器的输出为图4.5 自相关器匹配滤波器的主要特点匹配滤波器的主要特点图4.6 互相关器匹配滤波器的主要特点匹配滤波器的主要特点图4.7 输入为正弦信号时,相关器和匹配滤波器的输出波形在 时刻,匹配滤波器的输出与相关器的输出信号相等。5. 匹配滤波器输出的频谱函数与输入信号频谱函数的关系匹配滤波器输出的频谱函数与输入信号频谱函数的关系输入信号s(t)的频谱函数 模的平方,称为 s(t) 的能量频谱。 S随机过程的正交级数展开随机过程的正交级数展开l完备的正交函数集及确知信号完备的正交函数集及确知信号 s(t) 的正交级数展开的正交级数展开在(0,T)时间内满足式
11、(4.3.1),则函数集 构成相互正交的函数集。若不存在另一个函数 g(t),使 ,则正交函数集 是完备的正交函数集。展开系数随机过程的正交级数展开随机过程的正交级数展开接收信号接收信号 用正交级数展开用正交级数展开随机过程 x(t) 完全由其展开系数确定。 随机过程的卡亨南洛维展开:随机过程的卡亨南洛维展开:根据噪声干扰的特点,正确选择正交函数集根据噪声干扰的特点,正确选择正交函数集 ,以使展开系数,以使展开系数 之间是互不相关的随机变量。之间是互不相关的随机变量。随机过程的卡亨南随机过程的卡亨南-洛维展开洛维展开展开系数的均值:展开系数的均值:若要求展开系数若要求展开系数 互不相关,有:互
12、不相关,有:kx随机过程的卡亨南随机过程的卡亨南-洛维展开洛维展开正交函数集每个函数正交函数集每个函数 需满足积分方程:需满足积分方程:根据平稳噪声 n(t) 的自相关函数 ,求解上式的积分方程,得到特征函数 ,作为正交函数集 的坐标,对平稳随机过程 进行展开,展开系数 之间是互不相关的。核函数核函数特征函数特征函数特征值特征值白噪声情况下正交函数集的任意性白噪声情况下正交函数集的任意性任意取正交函数集任意取正交函数集 ,x(t) 的展开系数协方差的展开系数协方差当当 时,协方差时,协方差 。n(t) 是白噪声的条件下,取任意正交函数集是白噪声的条件下,取任意正交函数集 对平稳随机过程对平稳随
13、机过程 x(t)进行展开,展开系数进行展开,展开系数 之间都是互不相关的。之间都是互不相关的。参量信号随机过程的正交级数展开参量信号随机过程的正交级数展开把参量信号把参量信号 看作以看作以 为条件的信号,有为条件的信号,有其中,有展开系数互不相关 ,应满足高斯白噪声中确知信号波形的检测高斯白噪声中确知信号波形的检测主要内容:主要内容:简单二元信号波形的检测简单二元信号波形的检测一般二元信号波形的检测一般二元信号波形的检测M元信号波形的检测元信号波形的检测高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测技术路线技术路线分析信号模型分析信号模型推导信号状态的判决表示式推导信号状态的
14、判决表示式设计检测系统设计检测系统分析检测性能分析检测性能研究最佳波形设计研究最佳波形设计简单二元信号波形的检测简单二元信号波形的检测1. 信号模型信号模型接收信号中的信号分量 是能量为的确知信号,n(t) 是均值为零,功率谱密度为 的高斯白噪声。2. 判决表示式判决表示式用正交级数展开式 表示接收信号简单二元信号波形的检测简单二元信号波形的检测2. 判决表示式判决表示式第一步,用正交级数展开式 表示接收信号假设 下,有第二步,简单二元信号波形的检测简单二元信号波形的检测2. 判决表示式判决表示式假设 下,有假设 下和假设 下,展开系数的概率密度函数为:第二步,取前N项,构成似然比检验。简单二
15、元信号的波形检测简单二元信号的波形检测2. 判决表示式判决表示式简单二元信号的波形检测简单二元信号的波形检测2. 判决表示式判决表示式第二步,第三步,取 的极限,将离散判决表示式变成连续形式的判决表示式。简单二元信号的波形检测简单二元信号的波形检测3. 检测系统的结构检测系统的结构图4.8相关检测系统结构(相关接收机)图4.9匹配滤波器检测系统结构简单二元信号的波形检测简单二元信号的波形检测4. 检测性能分析检测性能分析检验统计量检验统计量 在假设在假设H0或假设或假设H1下,都是高斯随机变量。下,都是高斯随机变量。通过分析两种假设下的均值和方差,计算判决概率,通过分析两种假设下的均值和方差,
16、计算判决概率,并据此分析检测性能。并据此分析检测性能。 可以得到, , ,简单二元信号的波形检测简单二元信号的波形检测偏移系数:偏移系数:简单二元信号的波形检测简单二元信号的波形检测5. 最佳信号波形设计最佳信号波形设计在高斯白噪声条件下,简单二元确知信号波形的检测性能在高斯白噪声条件下,简单二元确知信号波形的检测性能由偏移系数由偏移系数d2决定,决定,d2取决于信号的能量取决于信号的能量Es,与信号波形无关。,与信号波形无关。图图4.10 接收机工作特性接收机工作特性 图图4.11检测概率检测概率PD与参数与参数d的关系的关系 简单二元信号的波形检测简单二元信号的波形检测6. 充分统计量的分
17、析方法充分统计量的分析方法第一个坐标函数选择为确知信号的归一化函数第一个展开系数两种假设下其余展开系数在两种假设下因此,展开系数因此,展开系数 x1是充分统计量。是充分统计量。参考线性代数的几何意义page133/174简单二元信号的波形检测简单二元信号的波形检测6. 充分统计量的分析方法充分统计量的分析方法利用充分统计量利用充分统计量 x1构造似然比检验构造似然比检验x1 是高斯随机变量,有是高斯随机变量,有返回返回一般二元信号波形的检测一般二元信号波形的检测1. 信号模型信号模型2. 判决表示式判决表示式用正交级数展开系数表示接收信号:一般二元信号波形的检测一般二元信号波形的检测2. 判决
18、表示式判决表示式取展开系数的前N项一般二元信号波形的检测一般二元信号波形的检测2. 判决表示式判决表示式一般二元信号波形的检测一般二元信号波形的检测3. 检测系统的结构检测系统的结构图4.12双路相关检测系统结构 图4.13双路匹配滤波器检测系统结构 一般二元信号波形的检测一般二元信号波形的检测4. 检测性能分析检测性能分析定义两个信号波形的相关系数为 一般二元信号波形的检测一般二元信号波形的检测4. 检测性能分析检测性能分析一般二元信号波形的检测一般二元信号波形的检测5. 最佳信号波形设计最佳信号波形设计信号的检测性能随偏移系数信号的检测性能随偏移系数d2的增加而增大。的增加而增大。因此,信
19、号的检测性能与以下两个因素有关。因此,信号的检测性能与以下两个因素有关。(1)信号相对于噪声的能量大小;)信号相对于噪声的能量大小;(2)信号之间的波形差异。)信号之间的波形差异。两个信号反相:两个信号反相: , ,两个信号正交:两个信号正交: ,两个信号满足:两个信号满足: , , 。一般二元信号波形的检测一般二元信号波形的检测6. 充分统计量的分析方法充分统计量的分析方法(1) 选择正交函数集中的第一个坐标函数为:(2) 根据Gram-Schmidt正交化方法,构造第二个坐标函数为:一般二元信号波形的检测一般二元信号波形的检测6. 充分统计量的分析方法充分统计量的分析方法由x1与x2构成的
20、二维矢量是充分统计量。x1和x2都是高斯随机变量,且相互统计独立。(4.4.53)(4.4.54)(4.4.55)一般二元信号波形的检测一般二元信号波形的检测6. 充分统计量的分析方法充分统计量的分析方法一般二元信号波形的检测一般二元信号波形的检测6. 充分统计量的分析方法充分统计量的分析方法一般二元信号波形的检测一般二元信号波形的检测6. 充分统计量的分析方法充分统计量的分析方法一般二元信号波形的检测一般二元信号波形的检测例题例题4.4.1图4.14 超越方程 的解 TTddtan一般二元信号波形的检测一般二元信号波形的检测7. 二元信号波形检测归纳二元信号波形检测归纳(1)对于高斯白噪声背
21、景的接收信号,进行正交展开的函数集 可以任意选择,展开系数xk是相互统计独立的高斯随机变量。 采用格拉姆-施密特正交化方法生成的正交函数集,可以获得有限维 且与假设Hj有关的充分统计量。(2)检验统计量 是高斯分布,因此判决概率P(HiIHj)完全由偏移 系数d2决定,即有效功率信噪比决定。 对于简单二元信号波形检测, 对于一般二元信号波形检测, 因此, 的信号互补关系是最佳的波形。 )(txlsENd02201012202ssssEEEENd)()(10tsts一般二元信号波形的检测一般二元信号波形的检测7. 二元信号波形检测归纳二元信号波形检测归纳(3)采用充分统计量分析方法的判决表示式:
22、图4.15 判决域划分示意图 一般二元信号波形的检测一般二元信号波形的检测7. 二元信号波形检测归纳二元信号波形检测归纳(3)分界线:直线的斜率:原信号差矢量的斜率:有:判决域分界线是垂直于信号间连线的直线!一般二元信号波形的检测一般二元信号波形的检测7. 二元信号波形检测归纳二元信号波形检测归纳(4)若二元信号假设的先验概率相等,采用最小平均错误概率准则, 则判决域分界线满足:分界线是两信号连线的垂直平分线。若进一步,两信号的能量相等,有分界线是两信号连线的垂直平分线,并通过判决域的原点。一般二元信号波形的检测一般二元信号波形的检测7. 二元信号波形检测归纳二元信号波形检测归纳 1f t 2
23、ft x t返回返回M元信号波形的检测元信号波形的检测系统每次发送系统每次发送M个可能信号的一个,个可能信号的一个,接收信号接收信号x(t)后,判断后,判断M个可能信号的哪一个。个可能信号的哪一个。 M元信号的检测元信号的检测 模式识别中的分类或辨识模式识别中的分类或辨识一般采用最小平均错误概率准则或一般贝叶斯准则。一般采用最小平均错误概率准则或一般贝叶斯准则。(4. 4. 73)M元信号波形的检测元信号波形的检测(1)根据格拉姆施密特正交化方法:正交函数集 中, 的坐标函数不必设计。考虑到M个信号中,可能有N个是线性不相关的,只需构造N个正交函数集。M元信号波形的检测元信号波形的检测(2)转
24、换为N维随机矢量问题M元信号波形的检测元信号波形的检测(3)平均错误概率最小的准则等价为最大后验概率准则由于 是N维联合高斯概率密度函数,有判决假设Hi成立。两边取对数;注意:C为对角阵M元信号波形的检测元信号波形的检测(3)平均错误概率最小的准则等价为最大后验概率准则若各假设的先验概率相等,有判决假设Hi成立。即M元信号波形的检测元信号波形的检测例 4.4.3图4.16四元信号检测判决域划分 图4.17四元信号检测系统结构 返回返回高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测加性高斯有色噪声背景中二元确知信号波形检测加性高斯有色噪声背景中二元确知信号波形检测1 卡亨南-洛
25、维展开法 根据噪声的自相关函数 选择合适的正交函数集 。2 将接收信号 先通过白化滤波器,其输出端噪声变为白噪声。nr tu 1,2,.kftk x t有色噪声:噪声过程有色噪声:噪声过程 的功率谱密度在频域上的分布不均匀。的功率谱密度在频域上的分布不均匀。一般采用高斯功率谱密度的模型。一般采用高斯功率谱密度的模型。 n t白色包含了所有的颜色,因此白噪声的特点就是包含各种噪声。白噪声定义为在无限频率范围内功率谱密度为常数的信号,这就意味着还存在其它“颜色”的噪声,即,其功率谱密度函数不平坦。大多数的音频噪声,如移动汽车的噪声,电钻的噪声,周围人们走路的噪声等等,其频谱主要都是非白色低频段频谱
26、。而且,通过信道的白噪声受信道频率的影响而变为有色的。高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测1 信号模型及其统计特性信号模型及其统计特性 均值为零,自相关函数为 的高斯有色噪声。 n tnr tu2 信号检测的判决表示式信号检测的判决表示式展开系数:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式展开系数之间的协方差函数:目的要得到正交函数集 的坐标函数 ,当 时,协方差函数 ,展开式中各展开系数 之间互不相关。 1,2,.kftk kftjk,0jkCov xxkx高斯有色噪声中确知信号波形的检测高斯有色噪
27、声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式取前 N 个展开系数构造 N 维矢量的似然比检验:再取 的极限,求信号波形下的判决表示式。N 高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式N维随机矢量 分别求解在两个假设下的概率密度函数:Nx高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式由前N项构成的似然比函数为:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式由前N项构成的似然比函数为:N 高
28、斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的判决表示式信号检测的判决表示式高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测2 信号检测的等效判决表示式信号检测的等效判决表示式若 ,即高斯白噪声环境02nNr tutu高斯白噪声环境下的结果是高斯有色噪声结果的特例。高斯白噪声环境下的结果是高斯有色噪声结果的特例。高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测3 检测系统的结构检测系统的结构双路相关器检测系统结构 高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测4 检测性能分析检测性能分析检验统计量:均值与方差:高斯
29、有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测4 检测性能分析检测性能分析确定函数 与 :0( )g t1( )g t带入上页均值函数,有:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测4 检测性能分析检测性能分析统计量的方差:统计量的方差:统计量的均值:统计量的均值:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测4 检测性能分析检测性能分析偏移系数偏移系数 :检测性能随着检测性能随着 的增大而增大。的增大而增大。 2d高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测4 检测性能分析检测性能分析高斯白噪声背景下:高斯白噪声背景下
30、: 1102gts tN 0002gtstN变为:变为:对比:对比:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测5 最佳信号波形设计最佳信号波形设计在一定条件下,寻找使在一定条件下,寻找使 达到最大的信号波形。达到最大的信号波形。2l构造辅助函数:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测5 最佳信号波形设计最佳信号波形设计若 与 分别代表 与 的最佳波形:求偏微分,并令偏微分的值为0:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测5 最佳信号波形设计最佳信号波形设计上式说明需满足条件:信号间的波形相关系数 。1 接下来确定信号的函数
31、形式:高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测5 最佳信号波形设计最佳信号波形设计等号左边的积分:有:信号 是上述齐次积分方程的特征函数,响应的特征值为 。 1y t高斯白噪声情况下, ,只要满足 。高斯有色噪声中确知信号波形的检测高斯有色噪声中确知信号波形的检测5 最佳信号波形设计最佳信号波形设计当约束条件为:信号能量之和为常数,积分方程中最小的特征值 对应的特征函数为最佳信号 ,同时 。 1y tmin 01yty t 01yty t 返回返回高斯白噪声中随机参量信号波形的检测高斯白噪声中随机参量信号波形的检测一般二元信号波形检测(带随机参量)一般二元信号波形检测(带随机参量)1 利用参量的先验概率密度函数,通过统计平均方法转换为 无条件概率密度函数。高斯白噪声中随机参量信号波形的检测高斯白噪声中随机参量信号波形的检测一般二元信号波形检测(带随机参量)一般二元信号波形检测(带随机参量)1 利用参量的先验概率密度函数,通过统计平均方法转换为 无条件概率密度函数。2 求得参量的最大似然估计量,替代信号中的随机参量或未知参量。