ch1_1离散信号与系统-北京交通大学陈后金教授信号处理教案解析.pptVIP

例: 已知x[k]={2, 1, -2, 1;}，y[k]={-1, 2, 1, -1; k=0,1,2,3}， 试计算互相关函数rxy[n] 和ryx[n]，以及自相关函数rx[n]。 解：根据序列的相关运算定义可得 相关在信号检测中的应用（检测信号） N=500; n=0:N-1; ns = randn(1,N); k=0:N-1; x0=k.*(0.95.^k); x1=k*0; x2=k*0; x=k*0; for k=200:N x1(k)=(k-200).*(0.9.^(k-200)); end for k=300:N x2(k)=(k-300).*(0.85.^(k-300)); end x=x0+x1+x2; y=x+ns; subplot(321);stem(n,x(n+1),.);title(Signal x);grid on subplot(323);stem(n,ns(n+1),.);title(Noise ns);grid on subplot(325);stem(n,y(n+1),.);title(y=x+ns);grid on [Ryy,n] = xcorr(y,y); subplot(322);stem(n,Ryy/N,r.);title(Ryy);grid on [Rxx,n] = xcorr(x,x); subplot(324);stem(n,Rxx/N,r.);title(Rxx);grid on [Rnn,n] = xcorr(ns,ns); subplot(326);stem(n,Rnn/N,r.);title(Rnn);grid on 相关在信号检测中的应用 相关在信号检测中的应用 由y的波形很难判断信号y 中是否存在周期性信号x。 由相关可判断y 中是否存在周期信号x及其周期。 N=500; ns = randn(1,N); k=0:N-1; n=0:N-1; x= 2*sin(0.01*pi*k); y=x+ns; subplot(321);stem(n,x(n+1),.);title(Signal x); grid on subplot(323);stem(n,ns(n+1),.);title(Noise ns); grid on subplot(325);stem(n,y(n+1),.);title(y=x+ns); grid on [Ryy,n] = xcorr(y,y); subplot(322);stem(n,Ryy/N,.);title(Ryy); grid on [Rxx,n] = xcorr(x,x); subplot(324);stem(n,Rxx/N,.);title(Rxx); grid on [Rnn,n] = xcorr(ns,ns); subplot(326);stem(n,Rnn/N,.);title(Rnn); grid on 相关在信号检测中的应用 相关在信号检测中的应用（回声检测） 相关在信号检测中的应用 load mtlb; x=mtlb; M=length(x); a=0.9; N=3000; x0=[x,zeros(1,N)]; xN=[zeros(1,N),x]; y=x0+a.*xN; yy=xcorr(y,y); subplot(2,2,1); stem(x0,.); title(Signal x[k]));grid on subplot(2,2,3); stem(y,.); title(y[k]=x[k] + 0.9*x[k-N]); grid on [Rxx,n] = xcorr(x0,x0); subplot(2, 2, 2); stem(n,Rxx,r.); title(Rxx[n]); grid on [Ryy,n] = xcorr(y,y); subplot(2, 2, 4); stem(n,Ryy,r.); title(Ryy[n]); grid on 相关在信号检测中的应用 load mtlb; x=mtlb; M=length(x); a=0.9; N=3000; x0=[x,zeros(1,N)]; xN=[zeros(1,N),x]; y=x0+a.*xN; yy=xcorr(y,y); subplot(2,2,1); stem(x0,.); title(Signal x[k]));grid on subplot(2,2,3); stem(y,.); title(y[k]=x[k] + 0.9*x[k-N]); grid on [Rxx,n] = xcorr(x0,x0); subplot(2, 2, 2); stem(n,Rxx,r.);