长沙理工大学数字信号处理实验三为下届学妹做贡献.doc

《长沙理工大学数字信号处理实验三为下届学妹做贡献.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、3-1t=0:0.001:1;x=sin(2*pi*100*t)+sin(2*pi*200*t)+1.5*rand(1,length(t));subplot(2,1,1);plot(x(1:50));y=fft(x,512);f=1000*(0:256)/512;subplot(2,1,2);plot(f,y(1:257));3-2N=input('输入序列的长度=');M=input('输入离散傅立叶变换长度=');f1=20;f2=20.5;f3=40;fs=100;t=0:1:N-1;x=sin(2*t*pi*f1/fs)+sin(2

2、*t*pi*f2/fs)+sin(2*t*pi*f3/fs);X=fft(x,M);stem(t,x);title('原始时域序列');xlabel('时间序号n');ylabel('振幅');pause;subplot(2,1,1);k=0:1:M-1;stem(k,abs(X));title('DFT抽样点的幅度');xlabel('频率序号k');ylabel('幅度');subplot(2,1,2);stem(k,angle(X));title('DFT抽样点的相位');xlabel('频率序号k');ylabel('相位');(1

3、)(2)(3)3-3（1）M=input('请输入M=');n1=0:M/2;n2=M/2:M;x1=n1;x2=M-n2;n=[n1,n2];x=[x1,x2];stem(n,x);（2）M=40;F=64;n1=0:M/2;n2=M/2:M;x1=n1;x2=M-n2;n=[n1,n2];x=[x1,x2];X=fft(x,F);stem(n,x);title('原始时域序列');xlabel('时间序号n');ylabel('振幅');pause;subplot(2,1,1);k=0:1:F-1;stem(k,abs(X));titl

4、e('DFT抽样点的幅度');xlabel('频率序号k');ylabel('幅度');subplot(2,1,2);stem(k,angle(X));title('DFT抽样点的相位');xlabel('频率序号k');ylabel('相位');(3)M=input('输入M=');F=input('输入离散傅立叶变换长度=');n1=0:M/2;n2=M/2:M;x1=n1;x2=M-n2;n=[n1,n2];x=[x1,x2];X=fft(x,F);X1=zeros(1,32);fork=1:32X1(k)=X(2*k);end;su

5、bplot(2,1,1);k=0:31;stem(k,abs(X1));title('DFT抽样点的幅度');xlabel('频率序号k');ylabel('幅度');subplot(2,1,2);stem(k,angle(X1));title('DFT抽样点的相位');xlabel('频率序号k');ylabel('相位');(4)M=40;F=64;n1=0:M/2;n2=M/2:M;x1=nx2=M-n2;n=[n1,n2];x=[x1,x2];X=fft(x,F);X1=zeros(1,32);fork=1:32X1(k)=X(2*

6、k);end;t=0:31;s1=(X1*exp(j*2*pi/32).^(t'*t))/32;stem(t,s1);sn=ifft(X1,32);stem(0:31,sn);(5)n=0:32*31;sn=s1(mod(n,32)+1);stem(n,sn);3-4h=fir1(30,125/500,boxcar(31));fs=1000;f1=100;f2=200;m=1024;y=fft(h,1024);fy=fs*(0:1023)/1024;subplot(2,1,1);plot(fy,abs(y));axis([0,500,0,1

7、.5]);w=exp(-j*2*pi*(f2-f1)/(m*fs));a=exp(j*2*pi*f1/fs);z=czt(h,m,w,a);fz=(f2-f1)*(0:1023)/1024+f1;subplot(2,1,2);plot(fz,abs(z));