山东大学第三四次数字信号处理实验报告.docxVIP

数字信号处理实验报告实验三、DFT和FFT1、实验目的：（1）掌握DFT/FFT及其性质（2）掌握采用DFT/FFT做信号频谱分析的方法（3）掌握利用DFT/FFT做序列的圆周卷积和线性卷积的方法2、实验内容（1）用Matlab编程实现pp167 习题3.6。在同一幅图上表示你的结果。实验程序：clc;clear all;close all;b1=[2,1,4,2,3];[H1,w1]=freqz(b1,1,whole); %用freqz函数求x（n）的DTFT变换subplot(2,1,1);plot(w1,abs(H1));hold on;y1=fft(b1,5);%用fft函数求x（n）的DFT变换n1=0:4;k1=2*pi*n1/5;stem(k1,abs(y1));title(对x(n)进行DTFT和DFT变换)b2=[2,1,4,2,3 0 0 0];[H2,w2]=freqz(b2,1,whole);subplot(2,1,2);plot(w2,abs(H2));hold on;y2=fft(b2,8);n2=0:7;k2=2*pi*n2/8;stem(k2,abs(y2))title(对补零后的x(n)进行DTFT和DFT变换)实验结果：X(k)等于X(e^(jw))中w=2*π*k/5，并且对x(n)补零后的DFT抽样点比之前更多。（2）用Matlab编程实现pp168 习题3.11。画图表示你的结果。实验程序：clc;clear all;close all;n=0:71;xn=cos(pi*n/6)+5*cos(pi*n/3)+4*sin(pi*n/7);y=fft(xn,72);%对x(n)进行72点DFTstem(n,abs(y));title(对x(n)做72点DFT);实验结果：x(n)的周期为84，对序列进行72点截断不能得到周期序列，进行频谱分析时，其频谱的周期延拓不是周期序列，所以会产生频谱泄露。（3）用Matlab编程实现pp170习题3.27。画图表示你的结果。实验程序：clc;clear all;close all;n=0:10;xn=cos(0.48*pi*n)+cos(0.52*pi*n);y=fft(xn,11); %对x(n)进行72点DFTsubplot(311)stem(n,abs(y));title(11点DFT);n1=0:10;n2=0:100;xn1=cos(0.48*pi*n1)+cos(0.52*pi*n1);xn2=[xn1 zeros(1,90)];%在x(n)后补上90个零y1=fft(xn2,101);%对补零后的x2(n)进行101点DFTsubplot(312)stem(n2,abs(y1));title(对x(n)补零后的101点DFT);n3=0:100;xn3=cos(0.48*pi*n3)+cos(0.52*pi*n3);y3=fft(xn3,101);%对101个数据的x3(n)进行101点DFT subplot(313)stem(n3,abs(y3));title(101点DFT);实验结果：时域抽样点数不变的情况下，在有效数据后增加零值点可以减小栅栏效应，提高DFT的计算分辨率，但不能提高频率分辨率。只有增加数据有效长度才能提高频率分辨率。（4）用Matlab编程实现pp208 习题4.15。画图表示你的结果。实验程序：clc;clear all;close all;t=0:0.0001:0.02;xt=cos(2*pi*300*t)+cos(2*pi*450*t)+sin(2*pi*1200*t)+sin(2*pi*2500*t);subplot(211)plot(t,xt) %画出连续时间信号x（t）title(连续时间信号x（t）)fs=2000;Ts=1/fs;t1=0:Ts:1;xt1=cos(2*pi*300*t1)+cos(2*pi*450*t1)+sin(2*pi*1200*t1)+sin(2*pi*2500*t1);load coef415; %调用低通滤波器[b,a]=sos2tf(SOS,G);% SOS和G是以直接II型结构给出的系数，需要通过调用函数 [B, A]=sos2tf(SOS,G)，转化为多项式形式的系数yn=filter(b,a,xt1); %用filter函数实现滤波[Yn,w]=freqz(yn,1,512,whole,fs);%求频率响应subplot(212)stem(w,abs(Yn),.)title(连续时间信号x（t）经过抽样和低通滤波后包含的频率分量)理想低通滤波器设计：实验结果：由图可知连续时间信号x（t）经过抽样和截止频率为800Hz的理想低筒低通滤波器后包含的频率分量由300