DFT及其性质的MATLAB实现

1、实验报告实验题目： DFT及其性质的MATLAB实现 姓名：崔金浩 院别：电气工程学院 专业：电子信息班级： 2 0 1 4级 三 班 指导教师： 许 春 和 时间： 2016 年 11 月 07 日电气工程学院实验六 DFT及其性质的MATLAB实现一、实验目的：（1）掌握DFT的定义和给定任意序列的DFT的MATLAB实现方法。（2）掌握DFT的主要性质和MATLAB环境下对性质的仿真程序。（3）更好的理解DFT及其性质在实际工程中的应用。二、实验仪器：PC机一台 MATLAB软件三、实验用MATLAB函数介绍： 可能涉及的函数abs( ); %求函数模值 angle ( ); %求函数相

2、角real ( ); %求函数的实部imag ( ); %求函数的虚部max( ); %求求函数的最大值length( ); %确定序列的长度四、实验内容：1.已知序列为，求此序列的DFT及IDFT并画出原信号及其DFT的X(k) 和IDFTX(k)的图像。xn=1,2,4,5,6,7,8,10;N=length(xn); n=0:N-1;k=0:N-1; subplot(3,1,1); stem(xn,.);%绘出xn xlabel(n); ylabel(|x(n)|);title(xn); Xk=xn*(exp(-1i*2*pi/N).(n*k);subplot(3,1,2); stem(

3、k,abs(Xk),.);%绘出xn的DFT xlabel(k);ylabel(|X(k)|); title(xn的DFT); subplot(3,1,3); stem(xn,.);%绘出xn xlabel(n);ylabel(x(n); title(xn的IDFT);2.已知，。（1）求再由的N点DFT获得。（2）由求，再求.用图形分别表示以上结果，将两种方法得到的进行比较。xn1=0,1,3,4,5,0; xn2=1,2,1,2,1,5; yn=2*xn1+3*xn2; N=length(yn); n=0:N-1;k=0:N-1; Xk=yn*(exp(-1i*2*pi/N).(n*k);

4、figure(1);stem(k,abs(Xk),.);xlabel(k);ylabel(|X(k)|);grid on;title(Y(k)的DFT); Xk1=xn1*(exp(-1i*2*pi/N).(n*k);figure(2);stem(k,abs(Xk1),.); xlabel(k);ylabel(|X1(k)|);grid on;title(X1(k)的DFT); Xk2=xn2*(exp(-1i*2*pi/N).(n*k); figure(3);stem(k,abs(Xk2),.);xlabel(k);ylabel(|X2(k)|);grid on;title(X2(k)的DF

5、T); Yk=2*Xk1+3*Xk2; figure(4);stem(k,abs(Yk),.); xlabel(k); ylabel(|Y(k)|); grid on; title(Y(k)的DFT);3. 已知序列为，设为圆周移位右移2位得到的新序列，（1）求出，再由的N点的DFT获得。 （2）先求的N点DFTX(k)，再由时域移位定理获得的N点的DFT.用图形分别表示以上结果，将、进行比较。xn=2,3,4,5,1,7,9;yn=(circshift(xn,2);figure(1);stem(k,abs(yn),.)xlabel(n);ylabel(y);grid on;title(yn);N=length(yn);n=0:N-1;k=0:N-1;Y1k=yn*(exp(-1i*2*pi/N).(n*k);figure(2);stem(k,abs(Y1k),.);%y1n的DFTxlabel(k);ylabel(|Y(k)|);grid on;title(Y1(k)的DFT);Xk=fft(xn);figure(3);stem(k,abs(Xk),.);%xn的DFT xlabel(k