samp9_3.m

数字信号处理的Matlab实现的随书光盘

%Samp9_3
clf;N=1000; Fs=500; %数据长度和采样频率
n=0:N-1;t=n/Fs;      %时间序列
Lag=200;        %最大延迟单位
x=sin(2*pi*10*t);    %周期信号x
y=0.5*sin(2*pi*10*t+90*pi/180);   %与x有90o相移的信号y
[c,lags]=xcorr(x,y,Lag,'unbiased');   %求无偏互相关
subplot(2,1,1);plot(t,x,'r');   %绘制x信号
hold on;plot(t,y,':');     %在同一幅图中绘y信号
legend('x信号', 'y信号')
xlabel('时间/s');ylabel('x(t) y(t)');
title('原始信号');grid on;
hold off
subplot(2,1,2),plot(lags/Fs,c,'r');   %绘制x,y的互相关
xlabel('时间/s');ylabel('Rxy(t)');
title(' 信号x和y的相关');grid on