自相关系数

您所在的位置:网站首页 样本长度怎么求 自相关系数

自相关系数

2024-01-13 18:07| 来源: 网络整理| 查看: 265

1、介绍 相关函数是描述信号X(s),Y(t)(这两个信号可以是随机的,也可以是确定的)在任意两个不同时刻s、t的取值之间的相关程度。两个信号之间的相似性大小用相关系数来衡量。定义: 在这里插入图片描述

称为变量 X 和 Y 的相关系数。若相关系数 = 0,则称 X与Y 不相关。相关系数越大,相关性越大,但肯定小于或者等于1.。相关函数分为自相关和互相关。下面一一介绍 自相关函数是描述随机信号 x(t) 在任意不同时刻 t1,t2的取值之间的相关程度。定义式: 在这里插入图片描述

主要性质如下:

(1)自相关函数为偶函数,其图形对称于纵轴。 (2)当s=t 时,自相关函数具有最大值,且等于信号的均方值,即 (3)周期信号的自相关函数仍为同频率的周期信号。 自相关是互相关的一种特殊情况.。互相关函数是描述随机信号 x(t)、y(t) 在任意两个不同时刻s,t的取值之间的相关程度,其定义为: 在这里插入图片描述

对于连续函数,有定义: 在这里插入图片描述

对于离散的,有定义:

在这里插入图片描述

从定义式中可以看到,互相关函数和卷积运算类似,也是两个序列滑动相乘,但是区别在于:互相关的两个序列都不翻转,直接滑动相乘,求和;卷积的其中一个序列需要先翻转,然后滑动相乘,求和。所以,f(t)和g(t) 做相关等于 f*(-t) 与 g(t) 做卷积。 在图象处理中,自相关和互相关函数的定义如下:设原函数是f(t),则自相关函数定义为 R(u)=f(t)f(-t),其中表示卷积;设两个函数分别是f(t)和g(t),则互相关函数定义为R(u)=f(t)*g(-t),它反映的是两个函数在不同的相对位置上互相匹配的程度。

2、互相关公式 自相关和互相关表示两个的线性度。越接近正负1线性度越好。若接近0表示线性关系不好,但有可能是其他关系,如平方等。 相关的计算公式其实是为了让协方差有量纲: 在这里插入图片描述 也就是下面公式: 在这里插入图片描述 3、自相关 协方差是会受到单位的影响的,而相关系数就是消除了量纲的影响。 这里讲的自相关系数可以说是根据最原始的定义引伸出来的。 其实自相关系数可以这么理解:把一列数据按照滞后数拆成两列数据,在对这两列数据做类似相关系数的操作。 看一个例子:

在这里插入图片描述

这组数据是求滞后数为2的自相关系数,则变成求{x1,x2,…,x8}和{x3,x4,…,x10}两者的“相关系数”,相关系数打引号是因为这个相关系数的公式和以往的有点不一样。下面看一下公式的对比:

在这里插入图片描述 (注意:上面公式左边大的求和号只是对分子求和) 要注意的是在计算自相管系数的时候 是使用的总体的均值, 可以看到他们除了 取得不一样, 几乎就是一样的。 所以,我们可以这么理解自相关系数, 她就是用来表达一组数据前后数据 (自己和自己) 的相关性的。

4、matlab中计算自相关系数 比如矩阵A = [1 2 3 4] ;xcorr(A) = 4 11 20 30 20 11 4 自相关函数是信号间隔的函数,间隔有正负间隔,所以n个长度的信号,有2n-1个自相关函数值,分别描述的是不同信号间隔的相似程度,并且该2n-1个自相关函数值关于n对称。 比如,上面的矩阵,最后得到7个结果,其中第4个是自己和自己相乘,最后相加的结果,值最大11+22+33+44 = 30 。而第三个和第五个分别是间隔正负1的结果也就是12+23+34 = 20,21+32+43 = 20 。第二个和第六个分别是间隔正负2也就是13+24=11,31+42 = 11。第一个和第七个分别是间隔正负3也就是14 = 4 ,41=4 。 而matlab中autocorr和xcorr有什么不同呢?autocorr是对序列减去均值后做的自相关,最后又进行了归一化。而且由于自相关本身是偶函数,autocorr只是取了以中点n为起始的后面n个序列。

clear ; %使用autocorr函数 A = [1 2 3 4] ; n = length(A) ; [ACF,lags,bounds] = autocorr(A,n-1) ; subplot(2,1,1) ; plot(lags(1:end),ACF(1:end)) ; title(‘autocorr求自相关’) ; %使用xcorr函数 B = A - mean(A) ;%减掉均值 [c,lags] = xcorr(B) ; d = c ./ c(n) ;%归一化 subplot(2,1,2) ; plot(lags(n:end),d(n:end)) ;%取中点n为起始的后面n个序列 title(‘xcorr求自相关’) ;

得到的函数图像一样,如图所示: 在这里插入图片描述


【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3