《DFT及其性质的MATLAB实现.docx》由会员分享，可在线阅读，更多相关《DFT及其性质的MATLAB实现.docx（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、实验报告+实验题目：DFT及其性质的MATLAB实现 姓名：崔金浩院别:电气工程学院专业:电子信息班级：2014级三班指导教师:许春和时间： 2016 年 U 月 07 日电气工程学院实验六DFT及其性质的MATLAB实现一、实验目的：(1)掌握DFT的定义和给定任意序列的DFT的MATLAB实现方法。(2)掌握DFT的主要性质和MATLAB环境下对性质的仿真程序。(3)更好的理解DFT及其性质在实际工程中的应用。二、实验仪器：PC机一台 MATLAB软件三、实验用MATLAB函数介绍:可能涉及的函数abs(); %求函数模值angle (); %求函数相角real (); %求函数的实部im
2、ag();%求函数的虚部max(); %求求函数的最大值length(); %确定序列的长度四、实验内容：1.序列为15) = 1,2,4,5,6,7,8,10,求此序列的DFT及IDFT并画出原信号x()及其DFT的X(k)和 IDFTX(k)的图像。xn=1,2,4,5,6,7,8/0;xn1011V4(N=length(xn);g 5-In=0:N-l;k=0:N-l;X, n LtILttsubplot(3,l,l);U12345678stem(xn,绘出 xnnxn 的 DFTxlabel(*n);50r:ylabel(|x(n)r);title(xn);Xk=xn*(exp(-1
3、i*2*pi/N).A(n*k);V01234567subplot(3,l,2);K xn的 IDFTstem(k,abs(Xk),绘出 xn 的 DFT10 111xlabel(k);S 5-tXylabcl(,|X(k)rj);titlefxn 的 DFT);U,i:12345678nsubplot。, 1,3);stem(xn,Y);%绘出 xnxlabel(nr);ylabel(x(n),);titleCxn 的 IDFT);2 .石() =0,1,3,4,5, x2 (n) = 1,2,1,2,1,5。(1)求)，() =2xI () + 3x2(n)再由 y()的 N 点 DFT
4、获得 Y(k)。(2)由芭(江 / 求 X，Xk),再求 丫(幻= 2X(K) + 3X2(K).用图形分别表示以上结果，将两种方法得到的y(Q进行比拟。xn 1=0,1,3,4,5,0;xn2=l,2,l,2,l,5;yn=2*xnl+3*xn2;N=length(yn);n=0:N-1;k=0:N-l;Xk=yn*(exp(-1 i*2*pi/N).A(n*k);figure(l);stem(k,abs(Xk),7);xlabel(k);ylabel(,|X(k)r);grid on;title(，Y(k)的 DFT);Xk l=xn 1 *(cxp(-l i*2*pi/N).A(n*k)
5、;figure(2);stem(k,abs(Xkl xlabel(k);ylabel(|Xl(k)r);grid on;title(Xl(k)的 DFT);Xk2=xn2*(exp(-1 i*2*pi/N).A(nf*k);figure(3);stem(k,abs(Xk2),.);xlabel(k);ylabel(lX2(k)n);grid on;title(X2(k)的 DFT);Yk=2*Xkl+3*Xk2;figure(4);stcm(k,abs(Yk),.);xlabel(k);ylabel(lY(k)n);grid on;X1(k)的 DFTtitle(Y(k)的 DFT);3 .序
6、列为x() = 2,3,451,7,9,设y()为x()圆周移位右移2位得到的新序列,(1)求出y(n)，再由)，()的N点的DFT获得K(幻。(2)先求x()的 N 点 DFTX(k),再由时域移位定理获得y()的N点的DFTr2a).用图形分别表示以上结果，将工伙)、xn=2,3,4,5,l,7,9;yn=(circshift(xn,2);figure(l);stem(k,abs(yn),.)xlabelCn);ylabel(y);grid on;title(yn);N=length(yn);n=0:N-l;k=0:N-l;Y1 k=yn*(exp(-1 i*2*pi/N).A(n*k); figure ；stem(k,abs(Y1 k),.);%y 1 n 的 DFT xlabel(k);ylabel(|Y(k)r);grid on;titleCYl(k)的 DFT);Xk=fft(xn);figure(3);stem(k,abs(Xk),7);%xn 的 DFTxlabel(k);ylabel(lX(k)r);五、思考题:以上实验内容假设采用快速傅里叶变换能否实现,应使用的函数是什么。批语:时间: