exa6_4.m

利用相关函数求信号功率谱,信号的自相关函数及不同信号的互相关函数

% exa060104_xcorr_xcov.m, for example 6.1.4  
% to test xcorr.m and xcov.m ;

clear;
N=256;
f=.1;
a1=5;
a2=3;
Mlag=N/4;

% 产生两个正弦加白噪声的数据；
x=a1*sin(2*pi*f*(0:N-1))+2*randn(1,N);
y=a2*sin(2*pi*f*(0:N-1))+randn(1,N);
subplot(2,2,1);
plot(x(1:N/2));
grid on;
title('原始信号1')
subplot(2,2,2);
plot(y(1:N/2));
grid on;
title('原始信号2')

% 求这两个数据向量的互相关函数；
corrxy=xcorr(x,y,Mlag,'unbiased');
subplot(2,2,3);
plot((-Mlag:1:Mlag),corrxy);
grid on;
title('互相关函数') 

% 求这两个数据向量的互协方差函数；
covxy=xcov(x,y,Mlag,'biased');
subplot(2,2,4);
plot((-Mlag:1:Mlag),covxy); 
title('互协方差函数')