时域信号的MATLAB仿真系统

1、xxxxxxxxX号分析与处理课程设计时域信号的MATLA昉真系统设计学生姓名x学号x所在学院通信工程学院专业名称通信工程班级2012级广播电视班指导教师x成绩四川师范大学成都学院二。一四年十二月四川师范大学成都学院课程设计报告课程设计任务书学生姓名x学生学号x学生专业通信工程广播电视方向学生班级广播电视班指导教师x职称讲师发题日期2014年10月30日完成日期2014年12月1日设计题目时域信号的MATLAB仿真系统设计设计目的：1,熟悉MATLA颤件平台2.掌握MATLA褊程方法、常用语句和可视化绘图技术3,熟悉和掌握常用的用于信号和系统时域仿真分析的MATLA函数具体任务及要求：对MAT

2、LA进行研究，了解MATLA歆件在信号处理中常用的几种功能。对于连续系统，了解信号描述及运算，连续时域系统分析和连续时域系统的仿真。对于离散系统，了解线性离散时间信号的基本运算，线性离散时域系统的响应，单位取样响应以及线性离散时间信号的卷积和运算。课程设计进度安排：1丁P内容安排时间1收集资料，确定课题的实施方案2014.11.1之前2软件设计2014.11.1-2014.12.13整理资料，编写程序报告2014.12,1-2014.12.204完成定稿，打印报告2014.12.15之前课程设文献：1管致中、夏恭恪、孟桥：信号与线性系统，高等教育出版社，2004.8,P63-P812郑君里、应

3、启衍、杨为理：信号与系统，高等教育出版社，2000.4,P93-P993上美、图四全：数字信号处理，四安电子科技大学出版社，2001.8,P45-P664程佩青：数字信号处理教程，清华大学出版社，2007.6,P92-P1135王沐然：MATLA的科学计算，电子工业出版社，2003.7,P281-P316指导教师签字院长审核签字I四川师范大学成都学院课程设计报告时域信号的MATLA昉真系统设计内容摘要：信号与系统是通信和电子信息类专业的核心基础课，其中的概念和分析方法广泛应用于通信、自动控制、信号与信息处理、电路域系统等领域。MATLABMatrixLABoratoy的缩写，它是一种以矩阵运算

4、为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。MATLA是美国MathWorks公司开发的大型工程计算软件，其顶尖的数值计算功能、强大的图像可视化功能、简洁易学的编程语言、可交互式的集成环境，深受工程技术人员欢迎。由于其内部包含规模极大、覆盖面极广的工具箱，因此应用范围十分广泛。利用其强大的信号处理工具箱可以生成信号，而且还可计算系统的响应，并可完成对连续时间系统的时域、频域及复频域的分析和仿真。MATLABT应用于连续系统中零输入状态响应的数值求解和符号求解，冲击响应和阶跃响应，以及卷积积分的求解和仿真。可运用于离散系统的时域分析和响应、单位取样响应、卷积运算、信号抽样的运算和

5、仿真。关键词：时域信号MATLAB仿真II四川师范大学成都学院课程设计报告ThetimedomainsignalofMATLABsimulationsystemdesignAbstract:Signalandsystemisthecoreofthebasiccourseofcommunicationandelectronicinformationspecialty,domainconceptandtheanalysismethodwidelyusedTongxin,automaticcontrol,signalandinformationprocessing,thecircuitdomains

6、ystemetc.MATLABistheabbreviationofMatrixLABoratoy,whichisamatrixbasedinteractiveprogramminglanguage,specificallytotheneedsforscientificandengineeringcalculationanddrawing.MATLABisalargeprojectdevelopmentAmericancalculationsoftwareMathofWorkscompany,theleadingfunctionofnumericalcalculation,powerfulim

7、agevisualizationfunction,simpleandeasytolearnprogramminglanguage,interactiveintegrationenvironment,bytheengineeringandtechnicalpersonnelwelcome.Duetoitsinternalcontainsscalegreat,widecoverageofthetoolbox,thereforetheapplicationrangeisveryextensive.Cangenerateasignalusingitspowerfulsignalprocessingto

8、olbox,butalsocanresponsecalculationsystem,andcancompletetheanalysisandSimulationoftimedomain,frequencydomainandcomplexfrequencydomainforcontinuoustimesystems.NumericalandsymbolicsolutionofMATLABcanbeappliedtocontinuoussystemzeroinputresponse,impulseresponseandstepresponse,andtheconvolutionintegralso

9、lutionandsimulation.Canbeusedinthetimedomainanalysisofdiscretesystemsandresponse,unitsamplingresponse,convolutioncomputationandSimulationofsignalsampling.Keywords：ThetimedomainsignalMATLABSimulationiii四川师范大学成都学院课程设计报告前言11 MATLAB简介21.1 MATLAB的特点和优势41.2 MATLAB在信号处理上的应用42信号分析51.1 信号的描述及运算51.1.1 信号的时移、反

10、褶和尺度变换51.1.2 积分运算71.2 信号时域分析81.2.1 冲激响应和阶跃响应81.2.2 零状态响应、零输入响应和完全响应93 时域信号分析及MATLA防真113.1 连续线性系统的时域分析及仿真113.2 离散线性系统的时域分析及仿真153.2.1 线性离散时间信号的基本运算153.2.2 零状态响应153.2.3 单位取样响应173.2.4 信号抽样183.2.5 离散系统仿真204 结束语22参考文献24IV四川师范大学成都学院课程设计报告时域信号的MATLA昉真系统设计前言信号与系统课程是一门实用性较强、涉及面较广的专业基础课，是通信、电子信息类专业本科生的必修课程。也是许

11、多通信、电子信息类专业硕士研究生入学必考课程，其中的概念和分析方法广泛应用于通信、自动控制、信号与信息处理、电路域系统等领域。该课程是将学生从电路分析的知识领域引入信号处理与传输领域的关键性课程。对后续专业课起着承上启下的作用。在信号与系统中，时域分析法是系统三大分析法之一，在系统分析中具有重要意义。这不仅是因为时域分析简单直观，概念容易理解，而且，还有一些非线性系统在限定范围与指定条件下，遵从线性特性的规律，另一方面，线性系统时域分析方法已经形成了完整的严密的体系，日趋完善和成熟。时域分析法直接分析时间变量的函数，研究系统的时间响应特性。这种方法的优点是概念清楚。对于输入、输出描述的数学模型

12、，可以利用经典法解常系数线性方程或差分方程；对于状态变量描述的数学模型，则需解矩形方程，在线性系统时域分析方法中，卷积方法最受重视，它的优点表现在许多方面。课题主要应用MATLAB勺数值运算和图形功能对线性时域系统进行分析，并利用其仿真功能对线性时域系统进行分析。MATLA前身为美国NewMexico大学计算机主任CleveMoler博士利用业余时间为方便同学使用EISPAC陆口LINPACKS程序编写得FORTRA犊口程序，并且使用MATrix和LABoratory两个单词的前三个字母组成的名字MATLAB初期的MATLAB业版本由SteveBangert主持开发编译解释程序，SteveKl

13、eiman完成图形功能的设计，JohnLittle和CleveMoler主持开发了各类数学分析的子模块；随后经过几年的快速发展，MATLA的mathworks公司在1993年突出了MATLAB的4.0版本；1995年推出了MATLAB4.2c(forWin3.x)版本；1995年推出了MATLA5.0版本；2000年10月推出了MATLA6.0版本；随后推出了MATLAB6.1、MATLAB6.5MATLAB7.0matlab7.1和MATLABR2006aK本。其功能越来越强大，应用领域也越来越广，成为当今应用比较广泛的工具软件之一。MATLAB自产生之日起就具有方便的数据可视化功能以将向量

14、和矩阵用图形表现出来并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视四川师范大学成都学院课程设计报告MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求应用广泛的模块集合工具箱。MATLABf巴工具箱延伸到了科学研究和工程应用的诸多领域诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统开发、定点仿真、DS叫通讯、电力系统仿真等都在工具箱Toolbox家族中有了自己的一席之

15、地。MATLAB1一个高级的矩阵阵列语言。它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序M文件后再一起运行语法特征。其与C+胡言极为相似而且更加简单更加符合科技人员对数学表达式的书写格式。使之更利于非计算机专业的科技人员使用。而且这种语言可移植性好、可拓展性极强这也是MATLABt够深入到科学研究及工程计算各个领域的重要原因。当前，科技飞速发展，新技术不断涌现，例如超宽带通信、人工神经网络信号处理、盲信号处理以及量子计算机与量子通信等，都可能成为今后实用的技术，相关的知识将成为“信号与系统”课程

16、中新的知识点，并可能对整个课程体系产生重大的影响。在今后必须重视跟踪新技术发展趋势。1 MATLAB简介MATLA%美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB口Simulink两大部分。MATLAB!矩阵实验室MatrixLaboratory的简称，和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAEM以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图

17、像处理、信号检测、金融建模设计与分析等领域。MATLAB勺基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLA来解算问题要比用C、FORTRAN语言完相同的事情简四川师范大学成都学院课程设计报告捷得多。并且mathwork也吸收了像Maple等软件的优点，使MATLAEK为一个强大的数学软件。在新的版本中也加入了对C、FORTRANC+、JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAEH数库中方便自己以后调用。此外许多的MATLA要好者都编写了一些经典的程序，用户可以直接进行下载就可以用。20世纪70年代，美国新墨西哥大学计算机科学系

18、主任CleveMoler为了减轻学生编程的负担，用FORTRAN写了最早的MATLAB1984年由Little、Moler、SteveBangert合作成立了的MathWorks公司正式把MATLABt向市场。至U20世纪90年代，MATLABB成为国际控制界的标准计算软件。历经升级，到2001年已经有了6.0版。现在MATLAB.1、6.5、7.0、7.X都已相继面世。早期的MATLAB在DO斜境下运行，1990年推出了Windows版本。1993年Mathworks公司又推出了MATLAB勺微机版本，充分支持在MicrosoftWindows界面下的编程。随着版本的升级，内容不断扩充，它的

19、功能越来越强大。特别是在系统仿真和实时运行等方面有很多新进展，更扩大了它的应用前景，是各种科学计算软件中使用频率最高的软件。MATLAEfe括拥有数百个内部函数的主包和三十几种工具包。工具包又可以分为功能性工具包和学科工具包。功能工具包用来扩充MATLAB勺符号计算、可视化建模仿真、文字处理及实时控制等功能。学科工具包是专业性比较强的工具包，控制工具包、信号处理工具包、通信工具包等都属于此类。1993年出现SIMULINK这是基于框图的仿真平台。SIMULINKS接在MATLAB环境上，以MATLAB勺强大计算功能为基础，以直观的模块框图进行仿真和计算。SIMULINK提供了各种仿真工具，尤其

20、是它不断扩展的、内容丰富的模块库，为系统的仿真提供了极大便利。在SIMULINK平台上，拖拉和连接典型模块就可以绘制仿真对象的模型框图，并对模型进行仿真。在SIMULINK台上，仿真模型的可读性很强，这就避免了在MATLA弱口使用MATLA峰令和函数仿真时需要熟悉记忆大量M函数的麻烦。对广大工程技术人员来说，这无疑是最好的福音。现在的MATLAB同时捆绑了SIMULINK的版本也在不断升级。从1993年的MATLAB4.0ZSIMULINK1.Ofi至U2001年的MATLAB6.1/SIMULINK4.1®,2002年即推出了MATLAB6.5/SIMULINK5.0版。MATLA

21、史经不再是单纯的“矩阵实验室”了，它已经成为一个高级计算和仿真平台。四川师范大学成都学院课程设计报告1.1 MATLAB的特点和优势特点：此高级语言可用于技术计算此开发环境可对代码、文件和数据进行管理交互式工具可以按迭代的方式探查、设计及求解问题数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数值积分二维和三维图形函数可用于可视化数据各种工具可用于构建自定义的图形用户界面各种函数可将基于MATLAB勺算法与外部应用程序和语言（如C、C+、Fortran、Java、COM以及MicrosoftExcel）集成优势： 友好的工作平台和编程环境 简单易用的程序语言 强大的科学计算机数据处理能

22、力 出色的图形处理功能 应用广泛的模块集合工具箱 实用的程序接口和发布平台 应用软件开发（包括用户界面）1.2MATLAB在信号处理上的应用MATLAB1一个包含大量计算算法的集合。拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C+,在计算要求相同的情况下，使用MATLAB勺编程工作量会大大减少。MATLAB勺这些函数集包括从最简单最基本的函数到诸如矩、特征向量、快速傅里叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性

23、方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅里叶变换和数据的统计分析、工程中的优化问题、复数的各种运算、三角函数和其他初等数学运算、多维数组操作及其建模动态仿真。正是因四川师范大学成都学院课程设计报告为MATLAIW有如此强大的科学计算机数据处理的功能，其在信号处理上就更能体现优势。2信号分析2.1 信号的描述及运算2.1.1 信号的时移、反褶和尺度变换信号的时移、反褶和尺度变换是针对自变量而言的，其数学表达式与波形变化之间存在一定的变化规律。信号f(t)的时移就是将信号数学表达式中的自变量t用t士to替换，其中to为正实数。因此，波形的时移变换是将原来的f(t)波形在时间轴上向

24、左或向右移动。已知信号f(t)的波形如下图所示，试用MATLA瑜令画出f(t-2)、f(3t)、f(-1)、f(-3t-2)的波形图。2-1f(t)波形图其实现的MATLAB(代码如下：clearall;t=2:0.01:4;ftl=functl(t-2);subplot(2,2,1);piot(t,ftl);xlabel('(a)f(t-2)');gridon;axis(-24-0.52);四川师范大学成都学院课程设计报告ft2=functl(3*t);subplot(2,2,2);plot(t,ft2);xlabel('(b)f(3t)');gridon;a

25、xis(-24-0.52);ft3=functl(-t);subplot(2,2,3);plot(t,ft3);xlabel('(c)f(-t)');gridon;axis(-24-0.52);ft4=functl(-3*(t+2/3);subplot(2,2,4);plot(t,ft4);xlabel('(d)f(-3t-2)');gridon;axis(-24-0.52);、f(3t)、f(-t)、f(-3t-2)的波形图图2-2f(t-2)四川师范大学成都学院课程设计报告2.1.2 积分运算用MATLA脸令求函数：y1=sin(x),y2=sinxlnx关

26、于x的一阶导数。其实现的MATLA觊代码如下：clearall;symsaxy1y2y1=sin(a*xA2);y2=x*sin*log(x);dy1=2*a*x*cos(a*xA2)dy2=sin(x)+log(x)*sin(x)+x*cos(x)log(x)_2.1.3 信号的相加与相乘运算信号的相加与相乘是指在同一时刻信号取值的相加与相乘。因此，MATLAB对于时间信号的相加与相乘都是是基于向量的点运算。已知f1(t)=sinCt,f2(t)=sin8Qt,试用MATLA琬令绘出f1(t)+f2(t)和f1(t)f2(t)的波形图，其中f=£=1Hz。2二其实现的MAYLAB?

27、序代码如下：clearall;f=1;t=0:0.01:3/f;f1=sin(2*pi*f*t);f2=sin(2*pi*8*f*t);subplot(2,1,1);plot(t,f1+1,':',t,f1-1,':',t,f1+f2);gridon;xlabel('(a)f1(t)+f2(t)');subplot(2,1,2);plot(t,f1,':',t,-f1,':',t,f1.*f2);gridon;xlabel('(b)f1(t)+f2(t)');四川师范大学成都学院课程设计报告H曲)图

28、2-3.(t)+f2(t)和.(t)f2(t)的波形图2.2信号时域分析2.2.1 冲激响应和阶跃响应在连续线性时域系统中，冲激响应和阶跃响应是系统特性的描述，对它们的分析是线性系统中极为重要的问题。输入为单位冲激响应函数6(t)所引起的零状态响应称为单位冲激响应，简称为冲激响应，用h(t)表示；输入为单位阶跃函数u(t)所引起的零状态响应称为单位阶跃响应，简称为阶跃响应，用g(t)表示。在MATLAB,对于连续时域系统的冲激响应和阶跃响应的数值解，可分别利用控制系统工具提供的函数impulse和step来求解。其调用格式如下：y=impulse(sys,t)y=step(sys,t)其中，t

29、表示计算系统响应的时间抽样点向量，sys表示线性时域系统模型。已知一连续时域系统的微分方程如下：y"(t)2y,(t)32y(t)=f(t)16f(t)(2-1)四川师范大学成都学院课程设计报告下面用MATLA瑜令绘出0MtM4范围内系统的冲激响应h(t)和阶跃响应g(t)。其实现的MATLA默码如下：clearall;t=0:0.001:4;sys=tf(1,16,1,2,32);h=impulse(sys,t);g=step(sys,t);subplot(2,1,1);plot(t,h);gridon;xlabel('时间/s');ylabel('h(t)

30、');title('冲激响应);subplot(2,1,2);plot(t,g);gridon;xlabel('时间/s');ylabel('h(t)');title('阶跃响应);其仿真波形如图：图2-4冲激响应h(t)和阶跃响应g(t)2.2.2 零状态响应、零输入响应和完全响应前面已述及线性连续系统可用常系数微分方程来描述。该系统的完全响应由零输入响应和零状态响应两部分组成MATLAB?号工具多!提供了dsolve函数,四川师范大学成都学院课程设计报告可实现常系数微分方程的符号求解，其调用格式如下：Dsolve('eq1&#

31、39;,'eq2','condl','cond2',v')其中，参数eq1、eq2、，表示各微分方程，它与MATLA符号表达式的输入基本相同，微分或导数的输入是用D1y、D2y、D3y、，表示y的一阶导数y,、二阶导数y.、三阶导数y,”、，；参数cond1、cond2、，表示各初始条件或起始条件；参数v表示自变量，默认为变量to可利用dsole函数来求解系统微分方程的零输入响应和零状态响应，进而求解完全响应。下面用MATLA瑜令求解彳分方程y”(t)+3y,(t)+2y(t)=x,(t)+3x(t),当输入x(t)=e'tu(t

32、),起始条件为y(0_)=1、y,(0_)=2时，系统的零输入响应，零状态响应及完全响应。求得零输入响应和零状态响应后，完全响应则为两者之和。其实现的MATLA毓序代码如下：clearall;eq=，D2y+3*Dy+2*y=0，;cond='y(0)=1,Dy(0)=2'yzi=dsolve(eq,cond);yzi=simplify(yzi)yzi=-3*exp(-2*t)+4*exp(-t)eq1=Dy+3*Dy+2*y=Dx+3*x;eq2=x=exp(-3*t)*Heaviside(t)'cond=y(-0.001)=0,Dy(-0.001)=0'yz

33、s=dsolve(eq1,eq2,cond);yzs=simplify(yzs.y)yzs=1/4*exp(-1/2*t)*Heaviside(t)yt=simplify(yzi+yzs)yt=-3*exp(-2*t)+4*exp(-t)+1/4*exp(-1/2*t)*Heaviside(t)subplot(3,1,1);ezplot(yzi,0,8);gridon;title('零输入响应);10四川师范大学成都学院课程设计报告subplot(3,1,2);ezplot(yzs,0,8);gridon;title('零状态响应);subplot(3,1,3);ezplot(

34、yt,0,8);gridon;title('完全响应);其波形如图所示：图2-5零输入响应，零状态响应及完全响应3时域信号分析及MATLAB!3.1连续线性系统的时域分析及仿真S2(1)已知连续系统的系统函数H(S)=2s，绘制系统函数的零极点S22s5分布图和计算系统函数的零极点位置和冲激响应时域波形。MATLA取现程序如下：%xh06碰察和分析系统零极点对系统时域特性的影响%绘制系统函数的零极点分布图和计算系统函数的零极点位置11四川师范大学成都学院课程设计报告b=12;a=125;H=tf(b,a);p,z=pzmap(H)subplot(2111);pzmap(H);subpl

35、ot(212);impulse(b,a);%xh07法制系统(a)的幅频响应曲线z=0;%p=-50-100;f1=0;%f2=80;k=0.1;%subplot(221)pz_plxy(f1,f2,f3,k,p,z);%法制系统(b)的幅频响应曲线Z=i*150-i*150;%P=-100+i*100-100-i*100;f1=0;%f2=80;k=0.1;%subplot(222)pz_plxy(f1,f2,k,p,z);%法制系统(c)的幅频响应曲线z=20+i*4020-i*40;%p=-20+i*40-20-i*40;f1=0;%定义系统（a）零点行向量定义幅频响应曲线频率范围定义频

36、率采样间隔绘制系统幅频响应曲线定义系统（b）零点行向量定义幅频响应曲线频率范围定义频率采样间隔绘制系统幅频响应曲线定义系统（c）零点行向量定义幅频响应曲线频率范围12定义频率采样间隔绘制系统幅频响应曲线定义系统（d）零点行向量定义频率采样间隔绘制系统幅频响应曲线四川师范大学成都学院课程设计报告f2=80;k=0.1%subplot(223)pz_plxy(f1,f2,k,p,z);%绘制系统(d)的幅频响应曲线z=;%p=-100；f1=0;f2=80;k=0.1;%subplot(224)pz_plxy(f1,f2,k,p,z);%仿真结果：p=-1.0000+2.0000i-1.0000-

37、2.0000iz=-2Tiima(aec).i(b和*mAXIS图3-1零极点图13四川师范大学成都学院课程设计报告颊军单位H2)小挽幅颊响在朗城图3-2四系统幅频特性曲线2连续时间LTI系统，其微分万程"D+1电g+25y(t)=%),若dt2出出f(t)=5sin3t),绘制系统的幅度响应特性、相位响应特性曲线图，利用系统仿真绘制稳定响应波形图。%应用freqs函数求频率响应并作图b=1,0；a=1,1,25;w=0:0.1:2*pi;H=freqs(b,a,w);HM=abs(H)HP=angle(H)subplot(211);plot(w/pi,HW);xlabel('

38、;omega单位：pi');title('幅度特性');axis(02-0.11.1);subplot(212);plot(w/pi,HP*180/pi);title('相位特性);14四川师范大学成都学院课程设计报告axis(02-100200);%End最后的仿真框图为：图3-3系统仿真框图图3-4系统输入输出波形3.2离散线性系统的时域分析及仿真3.2.1 线性离散时间信号的基本运算对离散时间序列实行基本运算可得到新的序列，这些基本运算主要包括加、减、乘、除、移位和反折等。两个序列的加减乘除是对应离散样点值的加减乘除,因此，可通过MATLAB（乘和点除、序

39、列移位和反折来实现，与连续时间信号处理方法基本一样。3.2.2 零状态响应离散线性时域系统可用下式表示，即15四川师范大学成都学院课程设计报告nm£%y(k-j)=£bmJf(ki),a0=1(2-2)j=0i=0MATLA中的函数filter可对上式差分方程在指定时间范围内的输入序列所产生的响应进行求解。函数filter的调用格式如下：y=filter(b,a,x)(2-3)其中x为输入的离散序列；y为输出的离散序列；y的长度和x的长度一样；b和a分别为差分方程右端与左端的系数向量。一离散线性时域系统的差分方程如下：3y(n)-4y(n-1)2y(n-2)=x(n)2x(

40、n-1)(2-4)下面用MATLA脸令绘出当?殷励信号为x(n)=-!u(n)时，该系统的零状态响应。2其实现的MATLA毓序代码如下：clearall;a=3-42;b=12;n=0:30;x=(1/2).An;y=filter(b,a,x);stem(n,y,'f川');gridon;xlabel('n');title('系统零状态响应y(n)');其仿真波形如下图：16四川师范大学成都学院课程设计报告图3-5零状态响应3.2.3 单位取样响应系统的单位取样响应定义为系统在仃(t)激励下系统的零状态响应，用h(n)表示。MATLABt解单位取

41、样响应可利用函数filter,另一种求单位取样响应的方法是利用控制系统工具提供的函数impz来实现。有一离散线性系统的差分方程如下：3y(n)-4y(n-1)2y(n-2)=x(n)2x(n-1)(2-5)下面利用MATLAB勺impz函数绘出该系统的单位取样响应。其实现的程序代码如下：clearall;a=3-42;b=12;n=0:30;impz(b,a,30);gridon;title('系统单位取样响应h(n)');其仿真波形如下图：17四川师范大学成都学院课程设计报告图3-6单位取样响应3.2.4 信号抽样信号抽样是连续时间信号分析向离散时间信号分析、连续信号处理向数

42、字信号处理的第一步广泛应用于实际的各类系统中。信号抽样也称为取样或采样，就是利用抽样脉冲序列p(t)从连续信号f(t)中抽取一系列的离散样值，通过抽样过程得到的离散样值信号称为抽样信号，用fs(t)表示。从数学上讲，抽样过程就是抽样脉冲p(t)和原连续信号f(t)相乘的过程，即：fs(t)=f(t)p(t)(2-6)因此，可以用傅里叶变换的频域卷积性质来求抽样信号fs(t)的频谱。信号在时域被抽样后，它的频谱是原连续信号的频谱以抽样和角频率为间隔周期的延拓,即信号在时域抽样或离散化，相当于频域周期化。升余弦脉冲信号为f(t)尊cos-(2-7)卜面用MATLA琳程实现该信号经冲激脉冲抽样后得到

43、的抽样信号fs(t)及其频1谱。参数，E=1,£=n则£代)=(1+cost)。当米样间隔Ts=1时，其实现的2MATLA毓序代码如下:18四川师范大学成都学院课程设计报告clearall;ts=1;dt=0.1;t1=-4:dt:4;ft=(1+cos(t1)/2).*(uCT(t1+pi)-uCT(t1-pi);subplot(2,2,1);plot(t1,ft);gridon;axis(-44-0.11.1);xlabel('Time(sec)');ylabel('f(t)');title('升余弦脉冲信号);N=500;k=-

44、N:1:N;W=pi*k/(N*dt);Fw=dt*ft*exp(-j*t1'*W);subplot(2,2,2);plot(W,abs(Fw);gridon;axis(-1010-0.21.1*pi);title('升余弦脉冲信号的频谱);t2=-4:Ts:4;fst=(1+cos(t2)/2).*(uCT(t2+pi)-uCT(t2-pi);subplot(2,2,3);plot(t1,ft,':');holdon;stem(t2,fst);gridon;axis(-44-0.11.1);xlabel('Time(sec)');ylabel(

45、'fs(t)');title('抽样后的信号)；19holdoff;Fsw=Ts*fst*exp(-j*t2'*W);subplot(2,2,4);plot(W,abs(Fsw);gridon;axis(-1010-0.21.1*pi);title('抽样信号的频谱);其波形如下图：升余弦脉冲僖号-F-一一，'T'1-十e:/*:娶0;J%nJQ1V-I-11-1"1*J-2024Time(se)施样后的信号：JC：用?1*b4,1i卬dA1*下05.f|=；升余弦脉冲信号的SS语.;一-10-50610抽样偌号的频谱，齐k，八

46、12飞,-t”-F'四川师范大学成都学院课程设计报告I10-6Q5103图3-7信号抽样波形图3.2.5离散系统仿真(1)已知某滤波器的差分方程:y(n)=x(n)+1.4x(n-1)+0.3x(n2)0.3y(n1)+0.2y(n-2)，用直接n型结构实现该滤波器，求冲激响应。20四川师范大学成都学院课程设计报告图3-8离散系统的仿真图MATLA毓序如下:n=0:20;subplot(212);stem(n,yout2,.');gridon;%直接口型冲激响应1614i11*114*111I1114j111i1I1iji.Iii卜A1idi*IL*r/施*1|111Iii卜t1*4iiiftI&111i1iV11»l11ii4»|I*i1i«-*k4J'-a4-bi«««=«4-L-V«->aJii*|*11iifII1111>l41*<11lt1*1*&l