数字信号处理-第四次作业new.docVIP

下载本文档

3
0
约9.51千字
约 17页
2016-10-05 发布于江西
举报
版权申诉

数字信号处理-第四次作业new.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数字信号处理-第四次作业new.doc

数字信号处理作业四 2.1 x=[1,1,1,0,0,0,0]; xp=[x,x,x,x,x,x,x]; y=fft(x); Xp=[y,y,y,y,y,y,y]; subplot(3,1,1);stem([0:length(x)-1],x);ylabel(x(n));grid on; subplot(3,1,2);stem([0:length(xp)-1],xp);ylabel(xp(n));grid on; subplot(3,1,3);stem([0:length(Xp)-1],Xp);ylabel(Xp(k));grid 2.2 x=[1,1,1,1,0,0,0]; y=[0,0,0,0,1,1,1]; xp=[x,x,x,x,x,x,x]; yp=[y,y,y,y,y,y,y]; X=fft(x); Xp=[X,X,X,X,X,X,X]; Y=fft(y); Yp=[Y,Y,Y,Y,Y,Y,Y]; F=Xp.*Yp; f=[ifft(F),ifft(F),ifft(F),ifft(F),ifft(F),ifft(F),ifft(F)]; subplot(4,1,1);stem([0:length(xp)-1],xp);ylabel(xp(n));grid on; subplot(4,1,2);stem([0:length(yp)-1],yp);ylabel(yp(n));grid on; subplot(4,1,3);stem([0:length(f)-1],f);ylabel(xp(n)*yp(n));grid on; subplot(4,1,4);stem([0:length(F)-1],F);ylabel(Xp(k).*Yp(k));grid on; 2.3 n1=0;n2=5; n=n1:n2; nt=length(n); k0=0;k1=1;k2=2;k3=3; nk0=abs(k0-n1)+1;nk1=abs(k1-n1)+1;nk2=abs(k2-n1)+1;nk3=abs(k3-n1)+1; x=zeros(1,nt); x(nk0)=1;x(nk1)=2;x(nk2)=3;x(nk3)=4; xp=[x,x,x,x,x,x]; n=0:length(xp)-1; subplot(3,1,1);stem(n-2,xp);ylabel(x1(n));grid on; subplot(3,1,2);stem(2-n,xp);ylabel(x2(n));grid on; subplot(3,1,3);stem(-n,xp);ylabel(x3(n) );grid on; 2.4 取xn=[1,2,0,0] N=4 x=[1,2,0,0]; X1=[ifft(x),ifft(x),ifft(x),ifft(x)]; x2=[1,2,0,0,1,2,0,0]; X2=[ifft(x2),ifft(x2),ifft(x2),ifft(x2)]; subplot(2,1,1);stem(0:length(X1)-1,X1);ylabel(X1);grid on;subplot(2,1,2);stem(0:length(X2)-1,X2);ylabel(X2);grid on; 由图可以看出X2(k)=X1(2k) 2.5 (a): N = 32; n = 0:N-1; xn = [1,zeros(1,29)]; xk = fftshift(fft(xn)); subplot(3,1,1);stem([-pi+2*pi/length(xn)/2:2*pi/length(xn):pi-2*pi/length(xn)/2]/pi,abs(xn));xlabel(rad),ylabel(|xn|);grid on;title(N = 32); subplot(3,1,2);stem([-pi+2*pi/length(xk)/2:2*pi/length(xk):pi-2*pi/length(xk)/2]/pi,abs(xk));xlabel(rad),ylabel(|xk|);grid on; subplot(3,1,3);stem([-pi+2*pi/length(xk)/2:2*pi/length(xk):pi-2*pi/length(xk)/2]/pi,angle(xk));xlabel(rad),ylabel(degree);grid on; N = 31; n = 0:N-1; xn = [1,zeros(1,29)]; xk = fftshift(fft(xn)); figure; su