《Matlab自相关函数和互相关函数的计算和作图》由会员分享,可在线阅读,更多相关《Matlab自相关函数和互相关函数的计算和作图(4页珍藏版)》请在金锄头文库上搜索。
1、自相关函数(Autocorrelation function,缩写 ACF)是信号处理、时间序列分析中常用的数学工具,反映了同一序列在不同时刻的取值之间的相关程度。自相关函数在不同的领域,定义不完全等效。在某些领域,自相关函数等同于自协方差(autocovariance)。信号处理在信息分析中,通常将自相关函数称之为自协方差方程。 用来描述信息在不同时间 的,信息函数值的相关性。,其中“*”是卷积算符, 为取共轭自相关函数的性质以下以一维自相关函数为例说明其性质,多维的情况可方便地从一维情况推广得到。 对称性:从定义显然可以看出 R(i) = R(i)。连续型自相关函数为偶函数当 f 为实函数
2、时,有: 当 f 是复函数时,该自相关函数是厄米函数,满足: 其中星号表示共轭。 连续型实自相关函数的峰值在原点取得,即对于任何延时 ,均有 。该结论可直接有 柯西-施瓦茨不等式得到。离散型自相关函数亦有此结论。 周期函数的自相关函数是具有与原函数相同周期的函数。 两个相互无关的函数(即对于所有 ,两函数的互相关均为 0)之和的自相关函数等于各自自相关函数之和。 由于自相关函数是一种特殊的互相关函数,所以它具有后者的所有性质。 连续时间白噪声信号的自相关函数是一个 函数,在除 = 0 之外的所有点均为 0。 维纳- 辛钦定理 (Wiener Khinchin theorem)表明,自相关函数和
3、 功率谱密度函数是一对傅里叶变换对: 实值、对称的自相关函数具有实对称的变换函数,因此此时 维纳-辛钦定理中的复指数项可以写成如下的余弦形式:白噪声的自相关函数为 函数:自相关函数和偏相关函数的问题在时间序列分析的研究中,首先是判别时间序列的稳定性,如果时间序列是平稳的就可以计算这些数据的自相关函数和偏相关函数。如果自相关函数是拖尾的,偏相关函数是截尾的,那麽数据符合 AR(P)模型。如果自相关函数是截尾的,偏相关函数是拖尾的,那麽数据复合 MA( Q )模型如果自相关函数和偏相关函数都是拖尾的,那麽数据复合 ARMA( P,Q )模型。自相关函数和互相关函数的 matlab 计算和作图 1.
4、 首先说说自相关和互相关的概念。这个是信号分析里的概念,他们分别表示的是两个时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度,即互相关函数是描述随机信号 x(t),y(t)在任意两个不同时刻 t1,t2 的取值之间的相关程度,自相关函数是描述随机信号 x(t)在任意两个不同时刻 t1,t2 的取值之间的相关程度。互相关函数给出了在频域内两个信号是否相关的一个判断指标,把两测点之间信号的互谱与各自的自谱联系了起来。它能用来确定输出信号有多大程度来自输入信号,对修正测量中接入噪声源而产生的误差非常有效.事实上,在图象处理中,自相关和互相关函数的定义如下:设原函数是f(t),则自相
5、关函数定义为 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 中,
6、求解 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);p
7、lot(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 高度相关