数字信号处理-课堂作业-频响特性及fft.doc

下载文档 降价啦

22
0
约 8页
2017-09-23 发布于河南
举报
版权申诉
保障服务

数字信号处理-课堂作业-频响特性及fft.doc

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

1. 程序如下： function[Hk]=cy(Hk,n) T=input(T=); Tp=input(Tp=); N=Tp/T; t=-4:T:3.75; t=t+eps; hn=sin(pi*t)./t/pi; hn1=[hn zeros(1,N-length(hn))]; k=[0:1:N-1]; n=[0:1:N-1]; Wn=exp(-j*2*pi/N); kn=k*n; Wn_kn=Wn.^kn; Hk1=hn1*Wn_kn; Hk=T*Hk1; f=(0:N/2-1); Hk2=abs(Hk(1:N/2)); plot(f,Hk2); 键入T=0.25 Tp=8 得键入T=0.25 Tp=16 得键入T=0.25 Tp=32 得 2 G=input(G=); n=1:G; xn=sin(0.48*pi*n)+sin(0.52*pi*n); subplot(3,1,1); stem(n,xn); N=length(xn); T=1/32; k=[1:1:N]; m=[1:1:N]; Wn=exp(-j*2*pi/N); kn=k*m; Wn_kn=Wn.^kn; Xk=T*xn*Wn_kn; f=(0:N/2-1); subplot(3,1,2); stem(f,abs(Xk(1:N/2))); 键入G=10 得键入G=100 得 n=1:10; xn=sin(0.48*pi*n)+sin(0.52*pi*n); xn1=[xn zeros(1,90)]; N=length(xn)+90; T=1/32; k=[1:1:N]; m=[1:1:N]; Wn=exp(-j*2*pi/N); kn=k*m; Wn_kn=Wn.^kn; Xk=T*xn1*Wn_kn; f=(0:N/2-1); G=1:100; subplot(3,1,1); stem(G,xn1); subplot(3,1,2); stem(f,abs(Xk(1:N/2))); 3 G=input(G=); n=1:G; xn=0.15*sin(0.125*pi*n)+sin(0.125*pi*2*n)-0.1*sin(0.125*pi*3*n); subplot(3,1,1); stem(n,xn); N=length(xn); T=1/32; k=[1:1:N]; m=[1:1:N]; Wn=exp(-j*2*pi/N); kn=k*m; Wn_kn=Wn.^kn; Xk=T*xn*Wn_kn; f=(0:N/2-1); subplot(3,1,2); stem(f,abs(Xk(1:N/2))); 键入G=16 键入G=24 键入G=32 补16个0 n=1:16; xn=0.15*sin(0.125*pi*n)+sin(0.125*pi*2*n)-0.1*sin(0.125*pi*3*n); xn1=[xn zeros(1,16)]; N=length(xn)+16; T=1/32; k=[1:1:N]; m=[1:1:N]; Wn=exp(-j*2*pi/N); kn=k*m; Wn_kn=Wn.^kn; Xk=T*xn1*Wn_kn; f=(0:N/2-1); G=1:32; subplot(3,1,1); stem(G,xn1); subplot(3,1,2); stem(f,abs(Xk(1:N/2)));