基于MATLAB的序列DFT实现

基于MATLAB的序列DFT实现

ID:47357880

大小:180.50 KB

页数:15页

时间:2019-08-05

基于MATLAB的序列DFT实现_第1页
基于MATLAB的序列DFT实现_第2页
基于MATLAB的序列DFT实现_第3页
基于MATLAB的序列DFT实现_第4页
基于MATLAB的序列DFT实现_第5页
资源描述:

《基于MATLAB的序列DFT实现》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、课程设计课程设计名称:数字信号处理课程设计专业班级:电子信息工程0604学生姓名:学号:20064300406指导教师:课程设计时间:2009年6月14日数字信号处理专业课程设计任务书学生姓名专业班级电信0604学号20064300406题目基于MATLAB的序列DFT实现课题性质其他课题来源自拟课题指导教师同组姓名主要内容根据已学的知识并结合MATLAB对下面三个序列计算在N=16和N=8时序列的N点DFT,并画出幅频特性曲线,最后用DFT理论解释为何两种N值下的DFT结果差别如此之大?复正弦序列余弦序列正弦序列任务要求1、写出设计原理和设计思路,画出程序流程图2、用MATLAB

2、画出序列的N点DFT及其幅频特性曲线3、用MATLAB画出N点DFT的幅频特性曲线参考文献1.程佩青著,《数字信号处理教程》,清华大学出版社,20012.SanjitK.Mitra著,孙洪,余翔宇译,《数字信号处理实验指导书(MATLAB版)》,电子工业出版社,2005年1月3.郭仕剑等,《MATLAB7.x数字信号处理》,人民邮电出版社,2006年审查意见指导教师签字:教研室主任签字:年月日说明:本表由指导教师填写,由教研室主任审核后下达给选题学生,装订在设计(论文)首页一设计内容与技术要求1.设计内容:设计题目:基于MATLAT的序列DFT实现内容要求:根据已学的知识并结合MA

3、TLAB对下面三个序列计算在N=16和N=8时序列的N点DFT,并画出幅频特性曲线,最后用DFT理论解释为何两种N值下的DFT结果差别如此之大?复正弦序列余弦序列正弦序列2技术要求:1.写出设计原理和设计思路,画出程序流程图2.用MATLAB画出序列的N点DFT及其幅频特性曲线3.用MATLAB画出N点DFT的幅频特性曲线二 设计原理及设计思路离散傅里叶变换(DFT),是连续傅里叶变换在时域和频域上都离散的形式,将时域信号的采样变换为在离散时间傅里叶变换(DTFT)频域的采样。在形式上,变换两端(时域和频域上)的序列是有限长的,而实际上这两组序列都应当被认为是离散周期信号的主值序列

4、。即使对有限长的离散信号作DFT,也应当将其看作经过周期延拓成为周期信号再作变换。时域中N点序列x[n]的DFT定义为X[k]=∑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)

5、;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(

6、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.

7、对于复正弦序列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('

8、X1(k)

9、')title('16点的DFT[x1(n)]')Subplot(2,2,3);stem(k,abs(X2k),'.');axis([0,20,0,20]);ylabe

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。