DFT及其性质的MATLAB实现.docVIP

实验报告 实验题目： 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 ( ); %求函数相角 real ( ); %求函数的实部 imag ( ); %求函数的虚部 max( ); %求求函数的最大值 length( ); %确定序列的长度 四、实验内容： 1.已知序列为，求此序列的DFT及IDFT并画出原信号及其DFT的X(k) 和IDFT[X(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(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); 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)的DFT); 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的DFT xlabel(k); ylabel(|Y(k)|); grid on; title(Y1(k)的DFT); Xk=fft(xn); figure(3); stem(k,abs(Xk),.);%xn的DFT xlabel(k); ylabel(|X(k)|); grid on; title(Y2(k)的DFT); 五、思考题： 以上实验内容若采用快速傅里叶变换能否实现，应使用的函数是什么。 批语： 时间：