试验原理有限长序列xn的离散傅里叶变换DFT.PPT

一、实验目的 1、掌握离散傅立叶变换（DFT）的计算方法及意义； 2、掌握实数序列的DFT系数的对称特点； 3、学习利用DFT计算程序计算IDFT的方法； 4、利用FFT对信号进行频谱分析。 5、分析DFT计算模拟信号中出现的问题 二、实验仪器 微型计算机 三、实验原理 有限长序列x(n)的离散傅里叶变换（DFT）： 三、实验原理 例1：求序列 的12点离散傅里叶变换（DFT） 三、实验原理 求上述频域序列X(K)的12点离散傅里叶反变换（IDFT） 三、实验原理 计算离散傅立叶变换矩阵(DFT) w=dftmtx(n) 其中n为采样点，返回W阵 DFT正变换： DFT反变换： 三、实验原理 例2 若x(n)=8(0.4)n是一个N=20的有限长序列，利用MATLAB计算它的DFT，并画出图形。 三、实验原理 有限长序列x(n)的快速傅里叶变换（FFT）： 三、实验原理 三、实验原理 例3 如果 是一个16点的有限序列，用MATLAB求其DFT 的结果，并画出其结果图。 三、实验原理 N=16; n=0:1:N-1; xn=sin(n*pi/8)+sin(n*pi/4); k=0:1:N-1; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.^nk; Xk=xn*WNnk; subplot(3,1,1); stem(n,xn); % x(n) subplot(3,1,2); stem(k,abs(Xk)); % X(k) WNnk1=WN.^(-nk); xn1=(Xk*WNnk1)./N; subplot(3,1,3); % X(k)的反变换 stem(n,real(xn1)); 三、实验原理 N=16; n=0:1:N-1; xn=sin(n*pi/8)+sin(n*pi/4); Xk=fft(xn,N); subplot(3,1,1); stem(n,xn); subplot(3,1,2); k=n; stem(k,abs(Xk)); n1=ifft(Xk, N); subplot(3,1,3); stem(n,real(xn1)); 例4 对连续的单一频率周期信号 按采样频率 采样，截取长度N分别选N=20和N=16，观察其DFT结果的幅度谱。 k=8; n1=[0:1:19]; xa1=sin(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) ; title( 20点序列); xlabel(t/T);ylabel(x(n)); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) ; title( 20点DFT); xlabel(k);ylabel(X(k)); n2=[0:1:15]; xa2=sin(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) ; title(16点序列); xlabel(t/T);ylabel(x(n)); xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4) stem(n2,xk2);title(16点DFT); xlabel(k);ylabel(X(k)); 三、实验原理 例5、利用FFT从受噪声污染的信号x(t)中鉴别出有用的信号。 （1）首先产生以下信号： x=sin(2*pi*100*t)+sin(2*pi*200*t) ; （2）给信号叠加随机信号（利用rand函数）; （3）利用FFT提取信号. 设采样频率1000Hz。采样周期0.001s。 t=0:0.001:1; %采样周期为0.001s,即采样频率为1000Hz; x=sin(2*pi*100*t)+sin(2*pi*200*t)+rand(size(t)); %产生受噪声污染的正弦波信号； subplot(2,1,1); plot(x(1:50)); %画出时域内的信号； Y=fft(x,512); %对X进行512点的傅立叶变换； f=1000*(0:256)/512; %设置频率轴（横轴）坐标，1000为采样频率； subplot(2,1,2)，plot(f,Y(1:257)); %画出频域内的信号 运行结果： 四、实验内容 1. 将例5中的最后的语句改为f=1000*(0:512)/512; subplot