互相关函数是信号分析里的概念,表示的是两个
时间序列之间的相关程度,即描述信号x(t),y(t)在任意两个不同时刻t1,t2的取值之间的相关程度。描述两个不同的信号之间的
相关性时,这两个信号可以是
随机信号,也可以是
确知信号。
对于
离散信号公式表示为R(n)=(1/N)∑[x(m)y(m+n)]其中m从0到N-1变化。
特殊地,若
离散信号为二进制信号,互相关函数应表示为R(n)=(A-D)/(A+D)其中A、D分别为x序列与循环移n位后的y序列之间相同的
码元数和不同的码元数。
由自相关的性质可知,函数的峰值在原点处,因此的峰值出现在点处,也就是函数和的
偏移量。
这个是信号分析里的概念,他们分别表示的是两个
时间序列之间和同一个时间序列在任意两个不同时刻的取值之间的相关程度。
自相关函数是描述
随机信号X(t)在任意两个不同时刻t1,t2的取值之间的相关程度;互相关函数给出了在
频域内两个信号是否相关的一个判断指标,把两
测点之间信号的互谱与各自的自谱联系了起来。它能用来确定
输出信号有多大程度来自
输入信号,对修正测量中接入
噪声源而产生的误差非常有效。
事实上,在
图象处理中,自相关和互相关函数的定义如下:设
原函数是f(t),则自相关函数定义为R(u)=f(t)*f(-t),其中*表示
卷积;设两个函数分别是f(t)和g(t),则互相关函数定义为R(u)=f(t)*g(-t),它反映的是两个函数在不同的
相对位置上互相匹配的程度。
上面代码是求自相关函数并作图,对于互相关函数,稍微修改一下就可以了,即把[a,b]=xcorr(x,'unbiased');改为[a,b]=xcorr(x,y,'unbiased');便可。
2. 实现过程:
在Matlab中,求解xcorr的过程事实上是利用Fourier变换中的
卷积定理进行的,即R(u)=ifft(fft(f)×fft(g)),其中×表示乘法(注:此公式仅表示形式计算,并非实际计算所用的公式)。当然也可以直接采用
卷积进行计算,但是结果会与xcorr的不同。事实上,两者既然有定理保证,那么结果一定是相同的,只是没有用对公式而已。下面是检验两者结果相同的代码: