matlab 已知离散时间lti系统的单位脉冲响应 求h(z),北京交通大学(数字信号处理研究生课程)ch1_1离散信号与系统分析.ppt...
北京交通大学(数字信号处理研究生课程)ch1_1离散信号与系统分析
例: 已知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[