数字信号处理上机实验dft综述.docVIP

本 科 实 验 报 告 实验名称： 数字信号处理上机实验 作业1：用DFT分析周期序列的频谱 任务：设周期序列截取 N 点长得到 （1）N=10，做10点DFT，得到 X1(k); （2）N=10，做100点补零DFT，得到 X2(k); （3）N=100，做100点DFT，得到 X3(k)。 要求：针对以上三种情况，分别输出|X1(k)|、|X2(k)|、|X3(k)|的图形，并进行比较、分析和讨论。 程序： clear all; n=0:1000; xn=cos(pi*0.48*n)+cos(pi*0.52*n); Xk1=fft(xn(1:10),10); X1=abs(Xk1); subplot(3,1,1); stem(X1,.); xlabel(k); ylabel(|X1(k)|); title(N=10,10点DFT); Xk2=fft(xn(1:10),100); X2=abs(Xk2); subplot(3,1,2); stem(X2,.); xlabel(k); ylabel(|X2(k)|); title(N=10,100点补零DFT); Xk3=fft(xn(1:100),100); X3=abs(Xk3); subplot(3,1,3); stem(X3,.); xlabel(k); ylabel(|X3(k)|); title(N=100,100点DFT); 运行结果： 分析： 从幅度谱中我们可以明显看出，X1(k)的相邻谱线间隔大，栅栏效应明显，频率分辨率低。 X2（k）是在采样数据10点后面补零再作DFT，降低了X1(k)中存在的栅栏效应，相邻谱线间隔变小，但并没有提高频率分辨率，没有改变频谱函数的主瓣宽度，只是将X1(k)的已存在的频谱函数进行内插，显示出更多细节。 X3(k)与X1(k)相比，信号采样数据增加，幅度谱栅栏效应减小，频率分辨率提高。X2（k）与X3(k)对照，X3(k)相当于时域上乘以一个较长的矩形窗，幅度谱的主瓣与X2(k)相比减小?(主瓣宽度为)，频谱泄露明显减轻，频率分辨率提高,但计算复杂度增加。 作业2：用DFT对模拟信号进行谱分析 任务：设 用FFT分析其频谱结构，选择不同的截取长度Tp，观察截断效应，试用加窗的方法降低旁瓣。 程序： fs=400; n=0:1000; xn=0.3*cos(200*pi*n/fs)+0.5*sin(100*pi*n/fs)+cos(50*pi*n/fs); Vk1=fft(xn(1:16),2048); Vk1=abs(Vk1)./max(abs(Vk1)); subplot(3,1,1); plot((0:length(Vk1)-1)*fs/length(Vk1),abs(Vk1)); xlabel(HZ); ylabel(|V1(k)|); title(截取时间长度Tp=0.04s，N=16,2048点DFT); Vk2=fft(xn(1:64),2048); Vk2=abs(Vk2)./max(abs(Vk2)); subplot(3,1,2); plot((0:length(Vk2)-1)*fs/length(Vk2),abs(Vk2)); xlabel(HZ); ylabel(|V2(k)|); title(截取时间长度Tp=4*0.04s，N=64,2048点DFT); Vk3=fft(xn(1:256),2048); Vk3=abs(Vk3)./max(abs(Vk3)); subplot(3,1,3); plot((0:length(Vk3)-1)*fs/length(Vk3),abs(Vk3)); xlabel(HZ); ylabel(|V3(k)|); title(截取时间长度Tp=16*0.04s，N=256,2048点DFT); 运行结果： 程序： fs=400; n=1:1000; xn=0.3*cos(200*pi*n/fs)+0.5*sin(100*pi*n/fs)+cos(50*pi*n/fs); Vk1=fft(xn(1:16).*hamming(16),2048); subplot(3,1,1); plot((0:length(Vk1)-1)*fs/length(Vk1),abs(Vk1)); xlabel(HZ); ylabel(|V1(k)|); title(N=16,2048点DFT); Vk2=fft(xn(1:64).*hamming(64),2048); subplot(3,1,2); plot((0:length(Vk2)-1)*fs/length(Vk2),abs(Vk2)); xlabel(HZ);