基于MATLAB序列DFT实现.docVIP

课 程 设 计 课程设计名称：数字信号处理课程设计 专 业 班 级 ： 电子信息工程0604 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 课程设计时间： 2009年6月14日 数字信号处理 专业课程设计任务书 学生姓名 专业班级 学号 题 目 基于MATLAB的序列DFT实现 课题性质 其他 课题来源 自拟课题 指导教师 同组姓名 主要内容 根据已学的知识并结合MATLAB对下面三个序列计算在N＝16和N＝8时序列的N点DFT，并画出幅频特性曲线，最后用DFT理论解释为何两种N值下的DFT结果差别如此之大？ 复正弦序列 余弦序列 正弦序列 任务要求 写出设计原理和设计思路，画出程序流程图 用MATLABMATLAB画出N点DFT的幅频特性曲线 参考文献 1．程佩青著，《数字信号处理教程》，清华大学出版社，2001 2．Sanjit K. Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书（MATLAB版）》，电子工业出版社，2005年1月 3．郭仕剑等，《MATLAB 7.x数字信号处理》 设计题目：基于MATLAT的序列DFT实现 内容要求：根据已学的知识并结合MATLAB对下面三个序列计算在N＝16和N＝8时序列的N点DFT，并画出幅频特性曲线，最后用DFT理论解释为何两种N值下的DFT结果差别如此之大？ 复正弦序列 余弦序列 正弦序列 2技术要求： 1. 写出设计原理和设计思路，画出程序流程图 2. 用MATLABMATLAB画出N点DFT的幅频特性曲线 二　设计原理及设计思路 离散傅里叶变换（），是连续傅里叶变换在时域和频域上都离散的形式，将时域信号的采样变换为在离散时间傅里叶变换（DTFT）频域的采样。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作经过周期延拓成为周期信号再作变换。 ∑x[n]exp[-j*2pi*kn/N] 三　程序流程图 输入相对应的参数 N=16；N=8 分别调用序列复正弦序列，正弦序列，余弦序列 通过matlab仿真得出结果 根据幅频特性曲线对结果进行分析 四　程序源代码 设计过程： 1.对于余弦序列 k=16; n1=[0:1:15]; xa1=cos(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) xlabel(t/T);ylabel(x(n)); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) xlabel(k);ylabel(X(k)); n2=[0:1:7]; xa2=cos(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) xlabel(t/T);ylabel(x(n)); xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4) stem(n2,xk2) xlabel(k);ylabel(X(k)); 所得原序列和幅频特性曲线： 2.对于正弦序列 k=16; n1=[0:1:15]; xa1=sin(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) xlabel(t/T);ylabel(x(n)); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) xlabel(k);ylabel(X(k)); n2=[0:1:7]; xa2=sin(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) xlabel(t/T);ylabel(x(n)); xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4) stem(n2,xk2) xlabel(k);ylabel(X(k)); 所得原序列和幅频特性曲线： 3.对于复正弦序列 N=16;N1=8; n=0:N-1;k=0:N1-1; x1n=exp(j*pi*n/8); X1k=fft(x1n,N); X2k=fft(x1n,N1); Subplot(2,2,1);stem(n,abs(X1k),.);axis([0,20,0,20]);ylabel(|X1(k)|) title(16点的DFT[x1(n)]) Subplot(2,2,3);s