资源描述:
《实验五 用DFT(FFT)对时域离散信号进行频谱分》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、07电信35号林烟泉(1)、分别以变换区间N=8,16,32,对x1(n)=R4(n)进行DFT(FFT),画出相应的幅频特性曲线;%DFT的MATLB计算xn=[1111];%输入时域序列向量xn=R4(n)Xk8=fft(xn,8);%计算xn的8点DFTXk16=fft(xn,16);%计算xn的16点DFTXk32=fft(xn,32);%计算xn的32点DFTk=0:7;wk=2*k/8;%产生8点DFT对应的采样点频率(关于π归一化值)subplot(3,2,1);stem(wk,abs(Xk8),'linewidth',2
2、);%绘制8点DFT的幅频特性图title('(a)8点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度');gridsubplot(3,2,2);stem(wk,angle(Xk8),'linewidth',2);%绘制8点DFT的相频特性图title('(b)8点DFT的相频特性图');gridon;xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])k=0:15;wk=2*k/16;%产生16点DFT对应的采样点频率(关于π归一化值)subplot(3,2,3);s
3、tem(wk,abs(Xk16),'linewidth',2);%绘制16点DFT的幅频特性图title('(c)16点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度');gridsubplot(3,2,4);stem(wk,angle(Xk16),'linewidth',2);%绘制16点DFT的相频特性图title('(d)16点DFT的相频特性图');gridon;xlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])k=0:31;wk=2*k/32;%产生32点D
4、FT对应的采样点频率(关于π归一化值)subplot(3,2,5);stem(wk,abs(Xk32),'linewidth',2);%绘制32点DFT的幅频特性图title('(e)32点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度');gridsubplot(3,2,6);stem(wk,angle(Xk32),'linewidth',2);%绘制32点DFT的相频特性图title('(f)32点DFT的相频特性图');gridon;xlabel('ω/π');ylabel('相位');axis([0,2
5、,-3.5,3.5]);图形输出如下所示:(2)、分别以变换区间N=8,16,32,对x2(n),x3(n)分别进行DFT(FFT),画出相应的幅频特性曲线;<1>、xn=[12344321];%输入时域序列向量x2Xk8=fft(xn,8);%计算xn的8点DFTXk16=fft(xn,16);%计算xn的16点DFTXk32=fft(xn,32);%计算xn的32点DFT%以下为绘图部分k=0:7;wk=2*k/8;%产生8点DFT对应的采样点频率(关于π归一化值)subplot(3,2,1);stem(wk,abs(Xk8),'l
6、inewidth',2);%绘制8点DFT的幅频特性图title('(a)8点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度');gridsubplot(3,2,2);stem(wk,angle(Xk8),'linewidth',2);%绘制8点DFT的相频特性图title('(b)8点DFT的相频特性图');gridxlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5])k=0:15;wk=2*k/16;%产生16点DFT对应的采样点频率(关于π归一化值)subplot(
7、3,2,3);stem(wk,abs(Xk16),'linewidth',2);%绘制16点DFT的幅频特性图title('(c)16点DFT的幅频特性图');xlabel('ω/π');ylabel('幅度');gridsubplot(3,2,4);stem(wk,angle(Xk16),'linewidth',2);%绘制16点DFT的相频特性图407电信35号林烟泉07电信35号林烟泉title('(d)16点DFT的相频特性图');gridxlabel('ω/π');ylabel('相位');axis([0,2,-3.5,3.5
8、])k=0:31;wk=2*k/32;%产生32点DFT对应的采样点频率(关于π归一化值)subplot(3,2,5);stem(wk,abs(Xk32),'linewidth',2);%绘制32点DFT的