Matlab信号处理：自相关和互相关分析
Matlab信号处理：自相关和互相关分析
在信号处理中，自相关和互相关是相关分析非常重要的概念，它们能分析一个信号或两个信号在时间维度的相似性，在振动测试分析、雷达测距和声发射探伤得到了广泛的应用。自相关分析的研究对象为一个信号，互相关分析的研究对象为两个信号。
本文以互相关分析为例，流程如下图所示，构建了一个余弦信号，并人为设置滞后延迟20个采样点，构建了另一个余弦信号，通过计算了两个余弦信号的互相关函数，分析了两个信号的相似性，找到了两个信号之间的时间延迟，在工程实际中常用于进行故障源定位。此外，对比matlab的互相关函数xcorr，采用手动计算的方法，计算互相关函数并验证本文手工计算方法的准确性。
1.具体案例
互相关分析是分析信号间时间前后相似的重要手段，它在故障源定位中得到了广泛的应用。本文分析的仿真信号如下图所示：
其中，y1为原始的余弦信号，y2为延迟了20个采样点（准确来说是采样间隔）的余弦信号，采样频率为100，信号长度为100，即采集时间为1s，获得时域波形如下图所示：
从上图中，能明显发现y2相对于y1延迟了差不多0.2s，一个采样间隔是1/100=0.01s，20个就是0.20s。采用matlab的互相关计算函数xcorr函数，获得互相关函数如下图所示：
上图为相关函数的茎叶图，使用stem函数。从图中能发现，存在正相关和负相关的最大点，表明了不同延迟或超前时两个信号的相关性最强。下图为手动计算互相关函数和matlab的xcorr函数的结果对比图：
从上图中能定性地发现，手动计算结果基本和matlab的结果保持一致。同时，定量计算的平均偏差（平均偏差为两数差的绝对值的平均数）结果如下，二者的偏差为0。
紧接互相关函数结果进行分析，如下图所示：
以最大正相关的点为例，发现互相关函数最大的峰值点（-20，20.1535），20表示使两个信号之间的相关程度最大的时间间隔是20，负号表示超前还是滞后，使用matlab的xcorr函数，代码具体为[c,lags] = xcorr(y1,y2)，c为互相关函数的y轴，lags为互相关函数的x轴。因此可以断定，-20的意思为y2比y1滞后20个采样间隔。
如果是正号的意思就是超前，通过分析手动计算互相关函数过程也能发现，负号的意思为滞后。代码后续增加了一个if判断函数，来自动分析互相关函数结果，结果如上上图所示。综上分析，互相关函数分析结果与信号人为设置结果相一致，验证了本案例的正确性。
2.具体代码
主函数main.m，用于互相关函数分析：
%% 信号的自相关和互相关分析
clc
clear all
close all
%%
fs=100; %采样频率
L=100; %信号长度
t=(0:L-1)/fs; %时间序列
y1=0.7*cos(2*pi*2*t-pi/4); %信号y1
y2=circshift(y1,20);
y2(1:20)=0;
%画图
figure;
subplot(211);plot(t,y1,'b-');xlabel('时间(s)');ylabel('幅值');
legend('y1');
subplot(212);plot(t,y2,'b-');xlabel('时间(s)');ylabel('幅值');hold on;
plot(t(21),y2(21),'Color','r','Marker','pentagram');
legend('y2','y1的起始点')
%计算互相关
[c,lags] = xcorr(y1,y2);
figure;
stem(lags,c)
%手动计算互相关
y11=zeros(3*length(y1)-2,1);
y22=zeros(3*length(y2)-2,1);
y11(length(y1):2*length(y1)-1)=y1;
y22(1:length(y2))=y2;
for i=1:2*length(y1)-1
c1(i)=sum(y11.*y22);
y22=circshift(y22,1);
end
%matlab函数和手动计算函数之间对比分析
figure;
plot(lags,c,'b-*');hold on;
plot(lags,c1,'o');
legend('matlab的xcorr函数','手动计算')
fprintf('matlab的xcorr函数获得相关函数与手动计算之间的差为%f \n',sum(abs(c1-c))/length(c1));
%分析结果
[c_max,index]=max(c);
%相关函数
figure;plot(lags,c,'b-');hold on;
plot(lags(index),c_max,'Color','r','Marker','pentagram','MarkerSize',10);
legend('相关函数','最大正相关的延迟位置')
%信号的对应关系
figure;
subplot(211);plot(t,y1,'b-');xlabel('时间(s)');ylabel('幅值');
legend('y1');
subplot(212);plot(t,y2,'b-');xlabel('时间(s)');ylabel('幅值');hold on;
plot(t(abs(lags(index))+1),y2(abs(lags(index))+1),'Color','r','Marker','pentagram');
legend('y2','y1的起始点')
%输出分析结果
if lags(index)<0
fprintf('y2比y1滞后%d 个采样点 \n',abs(lags(index)));
else
fprintf('y2比y1超前%d 个采样点 \n',abs(lags(index)));
end
xcorr是matlab的自相关和互相关的函数，函数只有一个输入，请大家参考官方使用介绍。互相关分析代码为：[c,lags] = xcorr(y1,y2)，c为互相关函数的y轴，lags为互相关函数的x轴，c的结果中，如果最大峰值的横坐标(即lags(max(c)))为-20，则表示y2比y1滞后了20个采样间隔（y1比y2超前了20个采样间隔），如果为20，则表示y2比y1超前了20个采样间隔（y1比y2滞后了20个采样间隔）。
上述案例是采用互相关函数，分析两个信号间的相关程度。通过设置第二个信号相对于第一个信号的延迟，利用互相关函数，确定潜在的延迟数，对比设置参考，验证互相关函数分析信号间相关程度的可行性，其原理常用于故障源定位。
此外，为了更好地理解互相关函数的计算原理，手动计算复现了互相关函数，对比matlab的自带函数，验证了手动计算的可行性。