1、数字信号课设 电话拨号燕山大学课 程 设 计 说 明 书题目： 移位与卷积运算 学院（系）： 电气工程学院 年级专业： 学 号： 学生姓名： 指导教师： 教师职称： 燕山大学课程设计（论文）任务书院（系）：电气工程学院 基层教学单位：仪器科学与工程系 学 号学生姓名专业（班级）设计题目移位与卷积运算设计技术参数其循环移位序列为。M取2和-2两种情况，N分别取10,20,30。，。设计要求（1） 绘出x(n)及其移位与循环移位序列波形（2） 绘出x1、x2波形，并完成有限长序列x1和x2的卷积与循环卷积并显示结果，并验证线性卷积与循环卷积的关系（circonv等）工作量数字信号处理方面资料MAT
2、LAB方面资料工作计划收集消化资料、学习MATLAB软件，进行相关参数计算，编写仿真程序、调试参考资料数字信号处理基础及matlab实现，北京希望电子出版社，周辉 董正宏编数字信号处理的matlab实现，科学出版社，万永革编；matlab7.x数字信号处理，人民邮电出版社，郭仕剑 王宝顺 贺志国编；matlab程序设计及在信号处理中的应用，西安交大出版社，聂祥飞，王海宝，谭泽富编；matlab5,x应用与技巧，科学出版社，蒙以正，柳成茂编；信号处理原理及应用，机械工业出版社，谢平，王娜，林洪彬编。指导教师签字基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。年
3、月 日 目 录摘 要 第3页第一章 概述第4页第二章 离散时间信号卷积的实现第5页第三章 仿真中用到的命令介绍第10页第四章 仿真程序与仿真图第11页第五章 线性卷积与循环卷积的验证与比较第15页课程设计总结 第17页参考文献 第18页评审意见表 第19页 摘 要数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。因此在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数转换器实现。而数字信号处理的输出经常也要变换到模拟域，这是通过数模转换器实现的。数字信号处理的算法需要利用计算机或专用处理设备如数字信号处理器（DSP）和专用集成电路（ASIC）等。数字信号处理技术及设备
4、具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等突出优点，这些都是模拟信号处理技术与设备所无法比拟的。 数字信号处理的核心算法是离散傅里叶变换(DFT)，是DFT使信号在数字域和频域都实现了离散化，从而可以用通用计算机处理离散信号。而使数字信号处理从理论走向实用的是快速傅里叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。广义来说，数字信号处理是研究用数字方法对信号进行分析、变换、滤波、检测、调制、解调以及快速算法的一门技术学科。但很多人认为：数字信号处理主要是研究有关数字滤波技术、离散变换快速算法和谱分析方法。随着数字电路与系
5、统技术以及计算机技术的发展，数字信号处理技术也相应地得到发展，其应用领域十分广泛。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在
6、数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其它编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 第一章 概 述数字信号处理就是对时间变量均为离散的信号进行变换、加工和处理的技术。数字信号处理主要研究的是数字系统（离散系统），核心算法是离散傅里叶变换(DFT)，是DFT使信号在数字域和频域都实现了离散化，从而可以用通用计算机处理离散信号。而使数字信号处理从理论走向实用的是快速傅里叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的
7、发展。由于客观世界中大量存在的是模拟信号，因此，在工程应用中常常是用“数字系统”处理模拟信号。数字信号处理技术是一门理论性、技术性都很强的科学，涉及知识面非常广，涉及信号分析、离散系统分析、综合实现技术等。数字信号处理的过程包含变换、运算和识别；对象既有确定信号由于随机信号；处理方法涉及时域、频域和复频域；实现方法不同于模拟处理系统，可以分为软件实现和硬件实现。软件实现是指按照原理和算法，在通用计算机上编程实现对数字信号的处理；硬件实现指的是按照具体要求，单独设计硬件结构。 本次课设主要研究线性卷积与循环卷积的计算过程，并设计程序用matlab实现卷积计算的过程。第二章 离散时间信号卷积的实现
8、在科学技术飞速发展的今天，现代通信技术正在蓬勃发展，如何快速准确、方便、易行地进行远距离信息传递，这时我们人类追求的目标。要实现这一目标必须进行离散信号卷积的计算。一、 线性卷积1、 定义及计算步骤设系统输入x（n）=（n），系统输出y（n）的初始状态为零，定义这种条件下系统输出为系统的单位取样相应，用h（n）表示，换句话说，单位取样相应即是系统对于（n）的零状态相应。用公式表示为h（n）=T（n） 设系统的输入用x（n）表示，按照式（3.2.3）表示成单位采样序列移位加权和为x（n）= 那么系统输出为X(n)=T根据线性系统的叠加性质y（n）=又根据时不变性质y（n）=式中的符号“*”代表卷
9、积运算，由式可见，这类卷积运算中有翻转、位移、相乘和相加，因此称这类卷积为线性卷积，简称卷积。式表示线性是不变系统的输出等与输入序列和该系统的单位取样响应的卷积。只要知道系统的单位取样响应，由公式可知，对于任意输入x（n）都可以求出系统的输出。计算离散信号卷积时，最常见的方法就是图解法，通过对信号h（n）的翻转、平移、乘积、取和等步骤完成计算。下面介绍图解法计算线性卷积的具体步骤：首先应该将x（n）和h(n)中的变量n用变量m替换，得到x（m）和h（m）。依据离散信号卷积运算过程（包括翻转，平移、乘积及取和），然后根据信号h（m）的图形作其翻转信号h（-m）的图形和平移信号h（n-m）的图形，
10、当n0时，序列右移；当n=maxN1,N2，则一般称此式所表示的运算为x1（n）和x2（n）的循环卷积，还可记为线性卷积是求离散系统响应的主要方法之一，许多重要应用都建立在这一理论基础上，如卷积滤波等。循环卷积，又称圆周卷积，它的计算速度要远远快于线性卷积，如果选择圆周卷积的长度L=N1+N2-1，则可以用圆周卷积取代线性卷积。循环卷积的计算过程为：先将周期化，形成，再反转形成，去主值序列移位n，得到，通常称之为的循环反转。将的循环反转序列移位n，得到。当n0时，序列右移；n=N+M-1时，以L为周期进行周期延拓才无混叠现象。此时取主值序列显然满足。由此证明循环卷积等于线性卷积的条件是L=N+
11、M-1。如果L=N+M-1，则可用DFT（FFT）计算线性卷积，计算框图如下：第三章 仿真中用到的命令介绍第一，conv(x1,x2)是用来实现卷积运算的，对x1序列和x2序列进行卷积，其中x1序列长为N，x2序列长为M,输出的结果的长度为L=N+M-1。第二，figure(n)，为当前工作空间的第n个图像。第三，subplot(a,b,c)函数，对图像空间进行子图划分，a表示行数，b表示列数，c表示子图序号。第四，stem(n,m)函数，画离散点的波形。第五，title(x1序列)，是规定子图名称，在图正中显示序列名称。 xlabel(n)和ylabel(幅度)，规定图的横纵坐标。第六，le
12、ngth(x)，取序列x的长度，x1=x1 zeros(1,N-length(x1)，若x1序列不够N则在x1后补零。 第四章 仿真程序与仿真图1、circonvt函数的源程序（均以 .m文件保存）function y= circonvt( x1,x2,N )%计算x1与x2的循环卷积，长均为N，即y(n)=sum(x1(m)*x2(n-m)modN)% 检查x1的长度 Detailed explanation goes hereif length(x1)N error(N必须=x2的长度)end%检查x2的长度if length(x2)N error(N必须=x2的长度)endx1=x1 z
13、eros(1,N-length(x1);%若长度不够补零x2=x2 zeros(1,N-length(x2);%若长度不够补零m=0:1:N-1;x2=x2(mod(-m,N)+1);H=zeros(N,N);for n=1:1:N H(n,:)=cirshfit(x2,n-1,N);endy=x1*H;end2、cirshfit函数的源程序function y= cirshfit( x,m,N )%将长度为N的序列x循环移m位,即y(n)=x(n-m)mod N)if length(x)N error(N必须=x的长度)endx=x zeros(1,N-length(x);n=0:1:N-1
14、;n=mod(n-m,N);y=x(n+1);end3、线性卷积程序n=0:9;x1=(0.8).n;x2=(0.6).n; figure(1),subplot(3,1,1);stem(x1);title(x1序列) xlabel(n);ylabel(幅度); figure(1),subplot(3,1,2);stem(x2);title(x2序列) xlabel(n);ylabel(幅度); y=conv(x1,x2)y = 1.0000 1.4000 1.4800 1.4000 1.2496 1.0774 0.9086 0.7549 0.6207 0.5066 0.2979 0.1739
15、0.1005 0.0572 0.0318 0.0171 0.0087 0.0039 0.0014 figure(1),subplot(3,1,3);stem(y);title(x1,x2线性卷积) xlabel(n);ylabel(幅度);线性卷积的仿真波形图：4、循环移位程序n=0:9;x1=(0.8).n;x2=(0.6).n; %生成x1,x2序列 figure(1),subplot(3,2,1);stem(n,x1);title(x1序列) xlabel(n);ylabel(幅度); %生成标题为x1序列的子图 figure(1),subplot(3,2,2);stem(n,x2);t
16、itle(x2始序列) xlabel(n);ylabel(幅度);N=10;x1=x1 zeros(1,N-length(x1); x2=x2 zeros(1,N-length(x2); xn1=cirshfit(x1,2,10)xn1 =0.1678 0.1342 1.0000 0.8000 0.6400 0.5120 0.4096 0.3277 0.2621 0.2097figure(1),subplot(3,2,3);stem(n,xn1);title(x1循环移位序列,m=2)xlabel(n);ylabel(幅度); xn2=cirshfit(x2,2,10)xn2 = 0.0168
17、 0.0101 1.0000 0.6000 0.3600 0.2160 0.1296 0.0778 0.0467 0.0280figure(1),subplot(3,2,4);stem(n,xn2);title(x2循环移位序列,m=2)xlabel(n);ylabel(幅度);xn1=cirshfit(x1,-2,10)xn1 = 0.6400 0.5120 0.4096 0.3277 0.2621 0.2097 0.1678 0.1342 1.0000 0.8000figure(1),subplot(3,2,5);stem(n,xn1);title(x1循环移位序列,m=-2) xlabe
18、l(n);ylabel(幅度);xn2=cirshfit(x2,-2,10)xn2 = 0.3600 0.2160 0.1296 0.0778 0.0467 0.0280 0.0168 0.0101 1.0000 0.6000figure(1),subplot(3,2,6);stem(n,xn2);title(x2循环移位序列,m=-2) xlabel(n);ylabel(幅度);循环移位的仿真波形图：5、循环卷积的程序n=0:9;x1=(0.8).n;x2=(0.6).n; %生成x1,x2序列N=10;x1=x1 zeros(1,N-length(x1); %把x1补为N-1的有限长序列x
19、2=x2 zeros(1,N-length(x2); %把x1补为N-1的有限长序列y=circonvt(x1,x2,N) %对x1和x2求循环卷积y = 1.2979 1.5739 1.5805 1.4572 1.2814 1.0946 0.9173 0.7588 0.6221 0.5066 %循环卷积序列figure(1),subplot(3,1,1);stem(y);title(N=10时，x1与x2的循环卷积序列)xlabel(n);ylabel(幅度);循环卷积的仿真波形图第五章 验证与比较线性卷积的计算结果为：y = 1.0000 1.4000 1.4800 1.4000 1.24
20、96 1.0774 0.9086 0.7549 0.6207 0.5066 0.2979 0.1739 0.1005 0.0572 0.0318 0.0171 0.0087 0.0039 0.0014循环卷积的结果为y = 1.0000 1.4000 1.4800 1.4000 1.2496 1.0774 0.9086 0.7549 0.6207 0.5066 0.2979 0.1739 0.1005 0.0572 0.0318 0.0171 0.0087 0.0039 0.0014 0 0 0 0 0 0 0 0 0 0 0通过对比线性卷积与循环卷积的计算结果，可以看出对于循环卷积当N取值不
21、同时结果不同，所以只有当循环卷积长度L=N+M-1时，循环卷积的计算结果才与线性卷积结果相同，循环卷积的周期延拓无混叠线性。因此线性卷积与循环卷积的关系得以验证。课程设计心得通过一周的课设，对数字信号处理和matlab有了全新的认识。通过课程设计我学习了很多，也收获了很多。这次我第一次接触matlab软件，从大一到大三一直都比较困惑，不知道自己所学的专业有什么意义，经过这次课设之后，我知道我们可从事的工作非常广泛，也产生了更加深刻的认识，对此我也产生了浓厚的兴趣。本次设计，首先针对题目进行分析，将所涉及的波形，频谱及相关函数做了研究，大体上能够把握了设计的流程以及思路。再通过查阅相关资料，能对
22、相关的知识做正确的记录，以便随时查看。在问题的分析阶段中，就原始信号的频率和载波信号的频率做了比较，确定了具体的方案后，在针对MATLAB中的有关画图处理函数进行学习和分析，这样就提高了学习的针对性，同时节约了设计的时间。在设计过程中，也遇到了许多的困难。由于一开始对知识的掌握程度不够，所以课程设计曾一度停滞不前，但在后面的学习及讨论过程中，最终解决了问题。此外在分析所设计的图中，根据相关的通信原理知识可以对结果作出判断，这样就提高了自己的相关知识，同时加深了对MATLAB的运用。由于本次设计运用了不同的知识，这样我就可以更好的将不同科目的知识进行联系学习，对牢靠的学习有着巨大的支持。这一个星
23、期的数字信号处理课程设计，通过老师们的讲解和自己的动手不但使我熟悉了数字信号课程设计的基本思想和基础知识，初步掌握了其应用软件matlab的使用，而且更为深入的体会了数字电路在现代高科技信息产业领域中的重要地位，当遇到了自己无法解决的困难与问题的时候，还请教老师给予指导和帮助。本着认真的态度，我回到宿舍后又对原理图进行了修改，并改变了相应的模块图，波形图，工作量非常大，但我认为这是自己必须要做的。这次实训给我最深的印象就是扩大自己的知识面，了解更多与本专业有关的科技信息，与时代共同进步，才能在将来成为有用的科技人才。通过简单的电路设计，提高了我得独立思考能力，并延伸了我在课堂上学到的知识，作为
24、一名工科学生，这样的数字信号处理课程设计是非常有意义的。在课程设计中，通过思考、发问，自己解惑并动手设计的过程，将课本上的知识学以致用，不仅强化了自己原有的知识体系，并扩展了自己的思维。经过这次课程设计，我相信在以后的课程设计中我们会吸取经验教训，做出更好的设计来。参考文献1、数字信号处理基础及matlab实现，北京希望电子出版社，周辉 董正宏编2、数字信号处理的matlab实现，科学出版社，万永革编；3、matlab7.x数字信号处理，人民邮电出版社，郭仕剑 王宝顺 贺志国编；4、matlab程序设计及在信号处理中的应用，西安交大出版社，聂祥飞，王海宝，谭泽富编；5、matlab5,x应用与技巧，科学出版社，蒙以正，柳成茂编；6、信号处理原理及应用，机械工业出版社，谢平，王娜，林洪彬编。燕山大学课程设计评审意见表指导教师评语：成绩： 指导教师： 年 月 日答辩小组评语：成绩： 评阅人： 年 月 日课程设计总成绩：答辩小组成员签字：年 月 日
