Matlab实验报告书.doc

广东海洋大学学生实验报告书（学生用表） 实验名称 离散傅里叶变换及其快速算法 课程名称 数字信号处理 课程号 学院(系) 信息学院 专业 电子信息工程 班级 电子1122 学生姓名 熊祥伟 学号 201211211228 实验地点 钟海楼 实验日期 2015.5.11 11.5.1 实验目的 （1）理解DFT算法，并能用MATLAB实现DFT。 （2）加深对FFT的理解，体会DFT和FFT之间的关系。 熟悉应用FFT实现两个序列的线性卷积的方法。 11.5.2 实验原理 N点序列x（n）的DFT和IDFT分别定义为 若将DFT变换的定义写成矩阵形式，则得到X=A*x，其中DFT变换矩阵A为 可以用函数U=fft（u，N）和uifft （U,N）计算N点序列的DFT正反变换。 成绩 指导教师 日期 注:请用A4纸书写，不够另附纸。 第　　　页，共　　　页 广东海洋大学学生实验报告书（学生用表） 实验名称 离散傅里叶变换及其快速算法 课程名称 数字信号处理 课程号 学院(系) 信息学院 专业 电子信息工程 班级 电子1122 学生姓名 熊祥伟 学号 201211211228 实验地点 钟海楼 实验日期 2015.5.11 11.5.2 实验内容 1.离散傅里叶变换 （1）用MATLAB求N=16的有限序列的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(2,1,1) stem(n,xn); subplot(2,1,2) stem(k,abs(Xk)); 成绩 指导教师 日期 注:请用A4纸书写，不够另附纸。 第　　　页，共　　　页 广东海洋大学学生实验报告书（学生用表） 实验名称 离散傅里叶变换及其快速算法 课程名称 数字信号处理 课程号 学院(系) 信息学院 专业 电子信息工程 班级 电子1122 学生姓名 熊祥伟 学号 201211211228 实验地点 钟海楼 实验日期 2015.5.11 （2）矩阵序列（n）=R5（n），求N分别取8,32时的DFT,最后会出结果图形。 function[Xk]=dft(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.^nk; Xk=xn*WNnk; 调用上面函数解题 N=8; x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); magX=abs(X); phaX=angle(X)*180/pi; k=(0:length(magX)-1); subplot(2,2,1); stem(n,x); ylabel(x(n)); subplot(2,2,2); stem(k,magX); axis([0,10,0,5]); ylabel(|X(k)|); N=32; x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); magX=abs(X); phaX=angle(X)*180/pi; k=(0:length(magX)-1)*N/length(magX); subplot(2,2,3); stem(n,x); ylabel(x(n)); subplot(2,2,4); stem(k,magX); axis([0,32,0,5]); ylabel(|X(k)|); 广东海洋大学学生实验报告书（学生用表） 实验名称 离散傅里叶变换及其快速算法 课程名称 数字信号处理 课程号 学院(系) 信息学院 专业 电子信息工程 班级 电子1122 学生姓名 熊祥伟 学号 201211211228 实验地点 钟海楼 实验日期 2015.5.11 快速傅里叶变换 已知一个8点的时域非周期离散阶跃信号，其起点n1=0，终点在呢=7，在n0=4前为0，n0以后为1，用N=32点进行FFT变换，做其时域信号图及其信号频谱图。 n1=0;n0=4;n2=7;N=32; n=n1:n2; w=[(n-n0)=0]; subplot(2,1,1); stem(n,w); i=0:N-1; y=fft(w,N); aw=abs(y); subplot