数字信号处理课程设计典型序列频谱分析.doc

各专业全套优秀毕业设计图纸 数字信号处理课程设计报告 题 目： 典型序列频谱分析 学 院： 班 级： 学 号： 姓 名： 联系方式： 指导教师： 报告成绩： 2014年 12月24日引 言1课程设计目的1课程设计要求2设计内容3一、三种序列的时域波形图31.1单位采样序列31.2实指数序列31.3矩形序列4二、 三种典型序列的傅里叶变换52.1傅里叶变换原理5三、 三种典型序列的幅度谱与相位谱53.1单位采样序列的幅度谱及相位谱53.2实指数序列的幅度谱及相位谱63.3矩形序列的幅度谱及相位谱7四、 三种典型序列的时移84.1时移原理84.2单位采样信号的时移及时移后的频谱图94.3实指数序列的时移及时移后的频谱图104.4矩形序列的时移及时移后的频谱图11五、 三种典型序列的频移125.1频移的原理125.2单位采样信号的频移及频移后的特性图125.3实指数序列的频移及频移后的特性图135.4矩形序列的频移及频移后的特性图14六、 结果分析156.1单位采样序列分析156.2实指数序列分析166.3矩形序列分析16七、 周期序列的设计167.1周期序列时域波形167.2周期序列的ft177.3周期序列的dfs187.4周期序列dfs与ft比较分析20心得体会21参考文献21引 言数字信号处理的主要研究对象是数字信号，且是采用运算的方法达到处理的目的的，因此，其实现方法，基本上分成两种实现方法，即软件和硬件实现方法。软件实现方法指的是按照原理和算法，自己编写程序或者采用现成的程序在通用计算机上实现，硬件实现指的是按照具体的要求和算法，设计硬件结构图，用乘法器加法器延时器、控制器、存储器以及输入输出接口部件实现的一种方法。显然前者灵活，只要改变程序中的有关参数，但是运算速度慢，一般达不到实时处理，因此，这种方法适合于科研和教学。后者运算速度快，可以达到实时处理要求，但是不灵活。目前dsp芯片已进入市场，且正在高速发展，速度高，体积小，性能优良，价格也在不断下降。可以说，用dsp芯片实现数字信号处理，正在变成工程技术领域的主要方法。用合适的dsp芯片，配有合适的芯片语言及任务要求的软件，来实现信号处理功能无疑是一种最佳的数字信号处理系统。课程设计目的1. 熟悉离散信号和系统的时域特性。2. 掌握数字信号处理的基本概念，基本理论和基本方法。3. 掌握序列快速傅里叶变换方法。4. 学会matlab的使用，掌握matlab的程序设计方法。课程设计要求1.1对于三种典型序列-单位采样序列、实指数序列、矩形序列要求： （1）画出以上序列的时域波形图； （2）求出以上序列的傅里叶变换； （3）画出以上序列的幅度谱及相位谱，并对相关结果予以理论分析； （4）对以上序列分别进行时移，画出时移后序列的频谱图，验证傅里叶变换的时移性质； （5）对以上序列的频谱分别进行频移，求出频移后频谱所对应的序列，并画出序列的时域波形图，验证傅里叶变换的频移性质。1.2自行设计一个周期序列要求： （1）画出周期序列的时域波形图； （2）求周期序列的dfs，并画出幅度特性曲线； （3）求周期序列的ft，并画出幅频特性曲线； （4）比较dfs和ft的结果，从中可以得出什么结论。 设计内容一、三种序列的时域波形图1.1单位采样序列公式程序：n=-25:25;y=zeros(1,25),1,zeros(1,25);stem(n,y)axis(-25,25,0,2);title(单位采样序列时域波形)特点： 单位采样序列也称为单位脉冲序列，仅在n=0时，数值才为1，其它时候取值全是0.它类似于模拟信号和系统中的单位冲激函数，但是不同的是在t=0时，取值无穷大，时取值为零，对时间t的积分为1。1.2实指数序列公式 a为实数特点当0a1时，该函数是单调递增函数，称为发散序列。程序n=0:40;a=1.2;x=power(a,n);stem(x,fill);xlabel(时间(n);ylabel(幅度x(n);title(实指数序列时域波形);1.3矩形序列公式式子中n为矩形序列的长度程序n0=0;n1=-10;n2=10;n3=50; n=n0:n3; x=(n=n0)&(n=n0)&(n0时每一个样本向右移动，称为x(n)的延时序列；当k=n0+n4)&(n=n0+n4)&(nlength(x) n=0:n-1; x=x zeros(1,n-length(x);endk=0:n-1;wn=exp(-j*2*pi/n);nk=n*k;wnnk=wn.nk;xk=x*wnnk;end利用dfs计算实现代码如下：fs=25;%设定采样频率n=256;n=0:n-1;t=n/fs;f0=10;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);figure(1);subplot(231);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(正弦信号y=2*pi*10t时域波形);x=fft(t,1024);xk=dfs(n,t,8);figure(2);plot(-1023:2048)/2048*16,abs(x) abs(x) abs(x),-);hold on;stem(-8:15,abs(xk2) abs(xk2) abs(xk2),linewidth,2);grid;set(gcf,color,w);图137.4周期序列dfs与ft比较分析 ft是傅立叶变换，它主要用于分析连续非周期信号，由于信号是非周期的，它必包含了各种频率的信号，所以具有时域连续非周期对应频域连续非周期的特点。离散的信号为周期序列时，严格的讲，傅立叶变换是不存在的，因为它不满足信号序列绝对级数和收敛（绝对可和）这一傅立叶变换的充要条件，但是采用dfs（离散傅立叶级数）这一分析工具仍然可以对其进行傅立叶分析。采用类似dfs的分析方法可以把有限长非周期序列假设为一无限长周期序列的一个主直周期，即对有限长非周期序列进行周期延拓，延拓后的序列完全可以采用dfs进行处理，即采用复指数基频序列和此有限长时间序列取相关，得出每个主值在各频率上的频谱分量以表示出这个“主值周期”的频谱信息。心得体会本次课程设计是关于三种典型序列的谱分析及特性，要求学生掌握典型序列的特性，能熟练使用matlab软件，学会运用matlab编写一些简单的程序，用图形表示出来，能把数字信号处理的知识学以至用，运用理论对典型序列的频谱加以分析，达到理论知识实践化。matlab主要用于数值运算， 可利用为数众多的附加工具箱，它也适合不同领域的应用，例如控制系统设计与分析、图像处理、信号处理与通讯、金融建模和分析等。它还是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。在通常情况下，可以用它来代替底层编程语言。在计算要求相同的情况下，使用matlab的编程工作量会大大减少。matlab的这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。通过这次的课程设计，使我对数字信号中的三种典型序列的谱分析及相关特性更加理解，从实验结果的分析中，对信号时移与频移的特性有了进一步的认识；另外，在这次课程设计中还加深了我对matlab软件