精通信号处理设计小Tips(8):检测淹没在噪声中的信号

精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第1张

  本文作者maxfiner,毕业于西安电子科技大学,拥有信号与信息处理专业硕士学位。maxfiner曾供职于华为通信技术公司无线通信部门,拥有多年的工程项目研发经验,同时兼备算法理论研究,仿真验证,以及对应的硬件设计实现能力;具备通信物理层开发设计各个方面的实战经验...

  相关函数的应用很广,比如噪声中信号的检测,信号中隐含周期性的检测,信号时延长度的测量等等。这一节专门讨论利用自相关函数检测淹没在噪声下的周期性信号。

  一个信号,由于噪声的影响,从波形上看,已经面目全非,杂乱不堪。如下图所示:

  

精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第2张

  基于信号的自相关累积可以检测是否有周期信号的存在,这是基于几个前提:

  第一,通常把噪声都看作加性的白噪声对待,因为这种噪声从理论上便于分析和处理。而白噪声的自相关函数有着非常独特的个性,就是自相关函数精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第3张只在精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第4张时有一定的数值,在精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第5张则为零值。这是白噪声本身的完全随机性决定的。也正是基于这个性质,我们才有可能从强噪声背景中提取出我们感兴趣的信号。

  第二,周期信号的自相关函数,仍是周期函数,并且周期和原函数周期相同。这从自相关函数的定义可以很容易推导出来。

  第三,白噪声和周期信号完全不相关,白噪声的完全随机性导致它和任何函数都不相关恋。它们之间的互相关函数可认为是零。这是很多理论推导的一个前提。

  基于自相关函数的定义,我们可以从公式上做一下简单的推导,看看添加了噪声后的信号的自相关函数到底是什么样子。假定我们的信号为精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第6张,白噪声信号为精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第7张,则添加了白噪声的信号为精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第8张

  那么,基于上面的三个前提,我们可以把精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第9张的自相关函数精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第10张表示成如下的形式。公式推导看上去很麻烦,其实一步一步都有严格的规律和理由,认真观察和分析,可以增进对自相关函数和互相关函数的理解,也可以加强公式推导的基础能力。

精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第11张

也就是说,精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第12张

  理论公式非常严密但不直观。还是找个实际例子验证下吧。构造一个正弦波周期信号和白噪声信号,看看它们叠加后的信号的自相关函数到底是个什么样子。

精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第13张精通信号处理设计小Tips(8):检测淹没在噪声中的信号,第14张

  自相关函数在0位置是一个相对较大的数值,这是白噪声自身相关累积导致的。在非零位置,则完全由周期信号来决定了。

  相应的matlab仿真验证代码如下,

  fs = 5e3;

  n = 0:1/fs:1;

  len = length(n);

  freq = 100;

  s = sin(2*pi*freq*n);

  n = randn(1,len);

  x = s + n;

  rr = xcorr(x,‘unbiased’);

  figure,subplot(121),plot(x);TItle(‘ÐźÅʱÓòµÄ²¨ÐÎ’);

  subplot(122),plot(rr);TItle(‘ÐźŵÄ×ÔÏà¹Øº¯ÊýµÄ²¨ÐÎ’);

  值得说明的是,相关函数采用matlab工具箱提供的函数xcorr,它有两种计算方式,一种叫“biased”,即有偏的。一种叫“unbiased”,即无偏的。这是由于实际计算相关函数时,实际数据的长度总是有限的,那么随着相关函数中m的增大,相关累积求和的样点数是逐渐减小的,当采用上面提到的自相关函数计算公式时,统一都除以N,导致相关函数随着m的增大会线性减,这就是有偏的计算方式,意思是计算值的统计上的均值和实际值之间并不一致,存在一定的偏差。这种计算方式导致当m的绝对值接近实际采集信号长度N时,误差会变得越来越大。若是把m的增大导致累积样点减少这个因素考虑进去,不再统一除以N,而是除以N-m,则相关函数不再会产生线性减小的问题,这就对应代码中的无偏的计算方式。

 

下期开讲——精通信号处理设计小TIps(9):估算信号在模拟通道的延时,敬请关注!

声明:电子发烧友网版权所有,谢绝转载!
 

往期回顾

精通信号处理设计小TIps(1):信号和信息
精通信号处理设计小Tips(2):数学的作用
精通信号处理设计小Tips(3):必须掌握的三大基石
精通信号处理设计小Tips(4):最频繁使用的几个信号
精通信号处理设计小Tips(5):三个应用广泛的数学概念
精通信号处理设计小Tips(6):卷积是怎么得到的?
精通信号处理设计小Tips(7):应用极其广泛的相关

 

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/dianzi/2437451.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-08-03
下一篇 2022-08-03

发表评论

登录后才能评论

评论列表(0条)

保存