1.首先说说自相关和互相关的概念。

    这个是信号分析里的概念,他们分别表示的是两个时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度,即互相关函数是描述随机信号x(t),y(t)在任意两个不同时刻t1,t2的取值之间的相关程度,自相关函数是描述随机信号x(t)在任意两个不同时刻t1,t2的取值之间的相关程度。

    自相关函数是描述随机信号X(t)在任意两个不同时刻t1,t2的取值之间的相关程度;互相关函数给出了在频域内两个信号是否相关的一个

判断指标,把两测点之间信号的互谱与各自的自谱联系了起来。它能用来确定输出信号有多大程度来自输入信号,对修正测量中接入噪声源而产生

的误差非常有效.

    事实上,在图象处理中,自相关和互相关函数的定义如下:设原函数是f(t),则自相关函数定义为R(u)=f(t)*f(-t),其中*表示卷积;设两个函数分别是f(t)和g(t),则互相关函数定义为R(u)=f(t)*g(-t),它反映的是两个函数在不同的相对位置上互相匹配的程度。

   通信系统中,自相关性决定多径干扰,互相关性决定了多址干扰。


那么,如何在matlab中实现这两个相关并用图像显示出来呢?

dt=.1;
t=[0:dt:100];
x=cos(t);
[a,b]=xcorr(x,'unbiased');
plot(b*dt,a)
上面代码是求自相关函数并作图,对于互相关函数,稍微修改一下就可以了,即把[a,b]=xcorr(x,'unbiased');改为[a,b]=xcorr(x,y,'unbiased');便可。
2. 实现过程:
     在Matalb中,求解xcorr的过程事实上是利用Fourier变换中的卷积定理进行的,即R(u)=ifft(fft(f)×fft(g)),其中×表示乘法,注:此公式仅表示形式计算,并非实际计算所用的公式。当然也可以直接采用卷积进行计算,但是结果会与xcorr的不同。事实上,两者既然有定理保证,那么结果一定是相同的,只是没有用对公式而已。下面是检验两者结果相同的代码:
dt=.1;
t=[0:dt:100];
x=3*sin(t);
y=cos(3*t);
subplot(3,1,1);
plot(t,x);
subplot(3,1,2);
plot(t,y);
[a,b]=xcorr(x,y);
subplot(3,1,3);
plot(b*dt,a);
yy=cos(3*fliplr(t)); % or use: yy=fliplr(y);
z=conv(x,yy);
pause;
subplot(3,1,3);
plot(b*dt,z,'r');
即在xcorr中不使用scaling。

3. 其他相关问题:
(1)相关程度与相关函数的取值有什么联系?

    相关系数只是一个比率,不是等单位量度,无什么单位名称,也不是相关的百分数,一般取小数点后两位来表示。相关系数的正负号只表

示相关的方向,绝对值表示相关的程度。因为不是等单位的度量,因而不能说相关系数0.7是0.35两倍,只能说相关系数为0.7的二列变量相关程度

比相关系数为0.35的二列变量相关程度更为密切和更高。也不能说相关系数从0.70到0.80与相关系数从0.30到0.40增加的程度一样大。
对于相关系数的大小所表示的意义目前在统计学界尚不一致,但通常按下是这样认为的:
相关系数     相关程度
0.00-±0.30   微相关
±0.30-±0.50  实相关
±0.50-±0.80  显著相关
±0.80-±1.00  高度相关

(2)matlab计算自相关函数autocorr和xcorr有什么不一样的?

    分别用这两个函数对同一个序列计算,为什么结果不太一样?因为xcorr是没有将均值减掉做的相关,autocorr则是减掉了均值的。而且,用离散信号做自相关时,信号截取长度(采样点N)不一样,自相关函数就不一样。
(3)xcorr是计算互相关函数,带有一个option的参数:
a=xcorr(x,y,'option')
option=baised时,是计算互相关函数的有偏估计;
option=unbaised时,是计算互相关函数的无偏估计;
option=coeff时,是计算归一化的互相关函数,即为互相关系数,在-1至1之间;
option=none,是缺省的情况。
所以想要计算互相关系数,可用'coeff'参数。
*************************************************************************
用这个xcorr函数作离散互相关运算时要注意,当x,y是不等长向量时,短的向量会自动填0与长的对齐,运算结果是行向量还是列向量就与x一样。
互相关运算计算的是x,y两组随机数据的相关程度,使用参数coeff时,结果就是互相关系数,在-1至1之间,否则结果不一定在这范围,有可能很大也有可能很小,这视乎x,y数据的大小,所以一般要计算两组数据的相关程度,一般选择coeff参数,对结果进行归一化。
所谓归一化简单理解就是将数据系列缩放到-1到1范围,正式的就是一种简化计算的方式,即将有量纲的表达式,经过变换,化为无量纲的表达式,成为纯量。变换式为X=(X实测--Xmin)/(Xmax-Xmin)。
一般来说选择归一化进行互相关运算后,得到结果绝对值越大,两组数据相关程度就越高。

数字信号处理中的自相关和互相关计算和物理意义(一)相关推荐

  1. 数字信号处理中的自相关和互相关计算和物理意义(二)

    在信号处理中,经常要研究两个信号的相似性,或者一个信号经过一段时间延迟后自身的相似性,以便实现信号检测.识别与提取等. 可用于研究信号相似性的方法称为相关,该方法的核心概念是相关函数和互相关函数. 1 ...

  2. 用MATLAB绘制国债NSS模型,Matlab在数字信号处理中的运用.ppt

    <Matlab在数字信号处理中的运用.ppt>由会员分享,可在线阅读,更多相关<Matlab在数字信号处理中的运用.ppt(68页珍藏版)>请在装配图网上搜索. 1.第七讲 M ...

  3. 数字信号处理中的 SIMD

    系列文章目录 数字信号处理中的 SIMD Neon intrinsics 简明教程 文章目录 系列文章目录 0. 1. 前言 2. SIMD 是什么 3. SIMD 伪代码示例 4. SIMD 是如何 ...

  4. 数字信号处理中卷积的计算

    数字信号处理的一条原则呢就是把信号分解成一个一个的脉冲信号,输入到系统之后得到输出响应,再把这些输出响应做一个线性的叠加就可以得到真是的响应了.这一点是非常重要的,不管是卷积还是傅里叶变换,本质就是这 ...

  5. 【数字信号处理】相关函数应用 ( 使用 matlab 计算相关函数 )

    文章目录 一.相关函数应用场景 1.生成高斯白噪声 2.信噪比 SNR 3.根据信噪比 SNR 求信号幅度 4.产生单载波信号及最终信号 5.求自相关函数及功率 6.matlab 完整代码 一.相关函 ...

  6. 序列的自相关和互相关计算

    -- Ref [1] [2] [3] ------------------------------------------------------------------------ 1.自相关和互相 ...

  7. 归一化数字角频率_数字信号处理中的各种频率

    在学习数字信号处理时,很多种频率很容易搞混淆,有模拟/数字/频率/角频率等等,也不是特别清楚不同频率之间的关系,希望这篇文件可以为各种频率来个了结. 4种频率及其数量关系 实际物理频率表示物理信号的真 ...

  8. c语言复数序列求自相关,序列的自相关和互相关计算

    -- Ref [1] [2] [3] ------------------------------------------------------------------------ 1.自相关和互相 ...

  9. 数字信号处理中的归一化频率

    4种频率及其数量关系 实际物理频率表示AD采集物理信号的频率,fs为采样频率,由奈奎斯特采样定理可以知道,fs必须≥信号最高频率的2倍才不会发生信号混叠,因此fs能采样到的信号最高频率为fs/2. 角 ...

最新文章

  1. QTableView中点击单元格弹出QComboBox
  2. 都在建议你不要直接使用 @Async 注解,为什么?
  3. uva111346Probability
  4. StringBuffer和StringBuilder使用方法比較
  5. ORA-12514: TNS:监听程序当前无法识别连接描述符中请(转)
  6. python 小海龟鼠标画图_Python小海龟画图
  7. 网络安全与渗透测试工具集合
  8. HTML DOM教程 37-HTML DOM Meta 对象
  9. 开源web应用防火墙 - Naxsi
  10. flutter实战1:完成一个有侧边栏的主界面
  11. zk和redis分布式锁比较
  12. 理解WebSocket心跳及重连机制以及加强版的reconnectingwebsocket.js实现
  13. 广州爱立信java笔试题_爱立信笔试经历
  14. ug草图精度怎么设置_UG怎么设置草图尺寸? UG草图尺寸的编辑方法
  15. 如何避开微信小程序的审核机制(实测有效)
  16. 深度linux双系统切换快捷键,苹果 Mac双系统如何切换?用Option键切换双系统的步骤分享...
  17. python远程主机强迫关闭了_[转] python 远程主机强迫关闭了一个现有的连接 socket 超时设置 errno 10054...
  18. 数学建模系列-模糊综合评价FCE
  19. 智能穿戴新方向华米科技血压手表相当于A类血压计申请标准
  20. C# windows服务部署到服务器上指令

热门文章

  1. 一款N-沟道耗尽型JFET晶体管 MPF102
  2. 利用pyBoard的实验来讨论部分MicroPython特性
  3. 航天智慧物流创意组-技术培训二期
  4. 2020春季学期作业提交统计处理
  5. 高频信号对LM386直流偏置的影响
  6. linux脚本中sed -i,Linux 中常用的sed命令
  7. java base64 压缩_在线等(急)Base64(带压缩)出问题
  8. 计算机桌面文件夹删除如何找回,电脑删除文件如何恢复 误操作的一剂后悔药...
  9. oracle的结构是什么,Oracle的软件结构是什么呢?
  10. 服务器论坛有哪些_SEO工作中,经常使用的无效外链有哪些?