DFT性质及DFT应用的研究

1、DFT性质及DFT应用的研究易长城生物医学工程专业2014级本科指导老师：王燕 李霞摘要：为了掌握DFT正变换及反变换之间的计算关系；观察和分析 DFS 应用中存在的混叠、泄漏、栅栏效应等问题和改善途径、方法；了解 DFT 实部、虚部之对应关系；验证 DFT 的若干性质；学习 Matlab 中子程序的编写和调用；运用matlab编程，研究DFT的性质和以及可能在应用出现的问题；通过上述方法，实现了DFT函数的编程和调用，详细探究了DFT的在时域和K域的能量守恒关系、复数序列的圆周相关定理，观察了DFT在采样频率不够时的混叠现象、没加窗之前的泄露现象以及栅栏效应。关键词： DFT matlab

2、时域 K域 Abstract：In order to understand the relationship between DFT positive transform and inverse transformation, we can observe and analyze the problems such as aliasing, leakage and fence effect in DFS application, and improve the way and method. We understand the correspondence between DFT real p

3、art and imaginary part. The use of matlab programming, research the nature of DFT and may be in the application of the problem; through the above method to achieve the DFT function programming and call, a detailed study of the DFT in the Time domain and K domain, the algebraic phenomenon of DFT in t

4、he case of insufficient sampling frequency, the leakage phenomenon and the fence effect before windowing are observed.Keywords：DFT matlab time domain K domain1.前言：了掌握DFT正变换及反变换之间的计算关系；观察和分析 DFS 应用中存在的混叠、泄漏、栅栏效应等问题和改善途径、方法；了解 DFT 实部、虚部之对应关系；验证 DFT 的若干性质；学习 Matlab 中子程序的编写和调用。应用前景及研究现状：MATLAB 的应用范围非常广，

5、包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用MATLAB函数集）扩展了MATLAB 环境，以解决这些应用领域内特定类型的问题。2.理论和方法2.1理论DFT的定义离散傅里叶变换（Discrete Fourier Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式，将信号的时域采样变换为其DTFT的频域采样。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作其周期延拓的变换。公式如下Xk=

6、DFTxn=n=0N-1x(n)WNnk,0kN-1xn=IDFTXk=1Nk=0N-1X(k)WN-nk, 0nN-1在用 DFT 计算连续信号时可能出现，频谱混叠；频谱泄露；栅栏效应等问题。MATLABMATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境

7、中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。2.2方法编写一个对xn=sin(n)n(-1)n,n=0,1,2,.,31进行DFT正变换和逆变换的程序；序列 x(n)的N点DTFT的物理意义是对X(ejw)的在0，2上进行 N 点等间隔采样。xn=R4(n) Xejw=DTFTxn=1-e-j4w1-e-jw,绘出相频曲线和幅频曲线，令w=2k/N，分别计算并图示x(n)的8点、16点DFT；验证能量守恒关系S1=n=0N-1x(n)

8、y*(n) 时域能量；S2=1Nk=0N-1XkY*k K域能量；其中x(n)=sin(n/8)+i*cos(2n) n=0,1,31 ym=DFTx(n) m=0,1,31混叠现象衰减正弦信号 x(t)=1+sin(7f0t)cos(2f0t), f0=120Hz, fs=200Hz,采样点数 N=64，绘出相频曲线和幅频曲线，观察混叠现象；泄露现象x(t)= cos(2f 0 t),f 0 =60Hz,fs=200Hz,采样点数 N=64；哈明窗函数：w(n)=0.54-0.46cos(2n/N),n=0,1,2,N-1 分别图示x(n)的频谱，窗函数频谱以及加窗后的频谱(x(t)*w(n

9、)即为加窗后函数)栅栏效应内容1的x(n)取07记为x(n)，并补零至L点，L分别为16,32,64；计算并显示x(n)的频谱X(k),与x(n)的频谱X(k)比较，观察补零的效果。验证圆周相关定理若Rxx=XkY*k;则有rxy=n=0N-1xny*n-mNRN(m),其中x(n),y(n)均为复数序列。3结果n=0:1:31;k=0:1:31;N=32;xn=sin(n).*n.*(-1).n;subplot(3,2,1);plot(xn);title('ԭͼ');subplot(3,2,2);stem(xn);title(

10、9;ԭͼ');WN=exp(-j*2*pi/N);nk1=n'*k;WNnk=WN.nk1;Xk=xn*WNnk;subplot(3,2,3);plot(abs(Xk);title('DFT');subplot(3,2,4);stem(abs(Xk);title('DFT');nk2=n'*k;WNnk=WN.(-nk2);xn1=(Xk*WNnk)/N;subplot(3,2,5);plot(real(xn1);title('IDFT')subplot(3,2,6);stem(

11、real(xn1);title('IDFT')图 1N1=1024;x1=(0:2*pi/1024:2*pi);X1=(1-exp(-1i*4*x1)./(1-exp(-1i*x1);N2=8;for i=1:N2; X2(i)=X1(i*N1/N2); x2(i)=2*pi*i/N2;endN3=16;for i=1:N3; X3(i)=X1(i*N1/N3); x3(i)=2*pi*i/N3;endsubplot(3,2,1),plot(x1,abs(X1),title('X1µÄ·ùƵÇ

12、50;Ïߣ¬N=1024');subplot(3,2,2),plot(x1,angle(X1),title('X1µÄÏàƵÇúÏߣ¬N=1024');subplot(3,2,3),stem(x2,abs(X2),title('X2µÄ·ùƵÇúÏߣ¬N=8');sub

13、plot(3,2,4),stem(x2,angle(X2),title('X2µÄÏàƵÇúÏߣ¬N=8');subplot(3,2,5),stem(x3,abs(X3),title('X3µÄ·ùƵÇúÏߣ¬N=16');subplot(3,2,6),stem(x3,angle(X3),title('X3&#

14、181;ÄÏàƵÇúÏߣ¬N=16'); 图 2clc;clear all;n=0:1:31;xn=sin(n/8)+j*cos(2*n);yn=dft(xn,32);yn1=conj(yn);for n=0:1:31 s1=sum(xn.*yn1);endXk=dft(xn,32);yn=Xk;图 3Yk=dft(yn,32);Yk1=conj(Yk);for k=0:1:31; s2=sum(Xk.*Yk1)/32;endclc,clear all;fs=200;f0=

15、120;N=64;Ts=1/fs;n=1:N;xt=(1+sin(7*pi*f0*n*Ts).*cos(2*pi*f0*n*Ts);Xk=dft(xt,64);subplot(2,2,1);plot(abs(Xk);title('·ùƵ')subplot(2,2,2);stem(abs(Xk);title('·ùƵ')subplot(2,2,3);plot(angle(Xk);title('ÏàƵ')subplot(2,2

16、,4);stem(angle(Xk);title('ÏàƵ')图 4clc,clear all;N=64;fs=200;f0=60;Ts=1/fs;n=1:N;x=cos(2*pi*f0*n*Ts); Xk=dft(x,64);subplot(3,2,1)plot(abs(Xk);title('Ô­º¯ÊýµÄƵÆ×') subplot(3,2,2);stem(abs(Xk);title('&

17、#212;­º¯ÊýµÄƵÆ×') n=0:(N-1);wn=0.54-0.46*cos(2*pi*n/N); Wn=dft(wn,64);subplot(3,2,3);plot(abs(Wn);title('´°º¯ÊýµÄƵÆ×')subplot(3,2,4);stem(abs(Wn);title('´°&#

18、186;¯ÊýµÄƵÆ×')x2=wn.*x;X2=dft(x2,64);subplot(3,2,5);plot(abs(X2);title('¼Ó´°ºóµÄƵÆ×')subplot(3,2,6);stem(abs(X2);title('¼Ó´°ºóµÄÆ

19、81;Æ×')图 5clc,clear all;n=0:7;N=8;xn1=sin(n).*n.*(-1).n;Xn1=dft(xn1,8);subplot(4,1,1);stem(abs(Xn1);title('N=8');n=0:15;N=16;xn2=sin(n).*n.*(-1).n;Xn2=dft(xn2,16);subplot(4,1,2);stem(abs(Xn2);title('N=16');n=0:31;N=32;xn3=sin(n).*n.*(-1).n;Xn3=dft(xn3,32);subplot(4,1,3);stem(abs(Xn3);title('N=32');n=0:63;N=64;xn4=sin(n).*n.*(-1).n;Xn4=dft(xn4,64);subplot(4,1,4);stem(abs(Xn4);title('N=64');图 6clc,clear all;x=1+1i,1-1i,2+1i,2-1i,3+1i,3-1i,4+1i,4-1i;y=1-1i