tshpro.m

来自「两个关于互相关函数的matlan程序。第一个求互相关」· M 代码 · 共 18 行

%两个sinc信号有0.2s的时移，用互相关函数计算时移的大小

clf
N=1000;n=0:N-1;Fs=500;t=n/Fs;
Lag=200;   %最大延迟单位数
x1=90*sinc(pi*(n-0.1*Fs));  %第一个原始信号延迟0.1s
y1=50*sinc(pi*(n-0.3*Fs));  %延迟0.3s
[c,lags]=xcorr(x1,y1,Lag,'unbiased');    %计算两个函数的互相关
subplot(2,1,1),plot(t,x1,'r');        %第一个信号
hold on;plot(t,y1,'b:');        %在同一幅图中绘制第二个信号
legend('信号 x','信号 y');     %绘制图例
xlabel('时间/s');ylabel('x(t) y(t)')
title('信号x和y');hold off
subplot(2,1,2),plot(lags/Fs,c,'r'); %绘制互相关信号
xlabel('时间/s');ylabel('Rxy(t)');
title('信号x和y的相关');