exa100800_xcov.m

来自「功率谱估计的一些matlab例程」· M 代码 · 共 29 行

%------------------------------------------------------------------------
% exa100800_xcov.m,  to test 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);
 pn=var(randn(1,1000));
 ps=a1^2;
 SNR(1:N)=10*log(ps/pn);
 subplot(3,1,1);
 plot(x(1:N/2));grid on;
 subplot(3,1,2);
 plot(y(1:N/2));grid on;
 
 % 求这两个数据向量的协方差函数；
 %covxy=xcov(x,y,Mlag,'biased');
 %subplot(3,1,3);
 %plot((-Mlag:1:Mlag),covxy); 
 subplot(3,1,3);
 plot(SNR(1:N),x(1:N/2));