北京交通大学(数字信号处理研究生课程)ch离散信号与系统概要.ppt

下载文档 降价啦

6
0
约8.98千字
约 56页
2017-05-14 发布于湖北
举报
版权申诉
保障服务

北京交通大学(数字信号处理研究生课程)ch离散信号与系统概要.ppt

1、本文档共56页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

北京交通大学(数字信号处理研究生课程)ch离散信号与系统概要

例: 已知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.);