信号与系统课程设计报告

课程设计任务书课程名称信号与系统课程设计题目连续LTI系统的分析专业通信工程班级通信1102学号姓名起止日期2013.06.26设计地点计通楼304指导教师职称分数辽宁石油化工大学计算机与通信工程学院一、课程设计的内容和要求用MATLAB软件设计连续LTI系统的分析。参考题目如下：1、已知信号f1(t)=u(t)-u(t-2)，f2(t)=cos(2πt)用MATLAB绘制f1(t)+f2(t)和f1(t)×f2(t)的波形。2、已知某连续系统的微分方程2y(t)y(t)8y(t)f(t)f(t)etu(t)试用MATLAB绘制出该系统的冲激响应和阶跃响应的波形。3、编程求解正弦函数f(t)sin(t)、取样函数f(t)=Sa(t)、门函数)1()1()(2tututg、双边指数函数te2信号的傅立叶变换式。4、已知某连续系统的的系统函数为H(s)(s23s2)/(8s42s33s25)，试用MATLAB求出该系统的零极点，画出零极点分布图5、设有两个稳定的LTI系统，可分别由下列微分来描述:'()3()3()3''()4'()()''()5()ytytftytytytftft试分别画出它们的系统频率响应的幅值和相位特性曲线。6、信号fl(t)和f2(t)如图所示。tt12111f2(t)f1(t)取t=0：0.005:2.5，计算信号f(t)=f1(t)+f2(t)cos(50t)的值并画出波形；而一可实现的实际系统的H(w)为H(w)=443223410()26.131()3.4142*10()2.6131*10()10jwjwjwjw用freqs函数画出H(w)的幅度和相位曲线。用lsim函数求出信号f1(t)和、f2(t)通过系统H(w)的零状态响应y1(t)和y2(t),并根据理论知识解释所得的结果。二、对课程设计成果的要求（包括课程设计说明书（论文）、图纸、实物样品等）1.课程设计要求符合《课程设计说明书规范》2.图表格式要规范。3上机运行结果，应满足课程设计任务书的基本要求，并有明确的结论。三、课程设计进度计划1）知识回顾和课程设计安排，集中进行，半天2）查阅文献，分散进行，1天3)数字信号处理设计并实施，两个半天4）整理文档，编写课程设计说明书，两天系（教研室）主任审查意见系主任（签字）：年月日摘要信号与系统课程设计是在学习完信号与系统以及复变函数与积分变换等课程的基础上设立的一门实践课。通过本次设计让学生学会使用MATLAB软件，同时学会建立m文件，并且在m文件下建立MATLAB函数和语句以及对其进行调试、运行、函数图像处理保存。进一步加深对信号的各种运算的理解，提高学生的综合分析问题和解决问题的能力，为后续课程和相关专业研究奠定基础AbstractSignalandsystemcourseisdesignedafterstudyingsignalsandsystems,andcomplexfunctionandintegraltransformcourseestablishedonthebasisofapracticallesson.ThroughthedesigntoletthestudentslearnhowtousetheMATLABsoftware,atthesametimelearnhowtobuildmfile,andtheMATLABfunctionandstatementinthe.mfileandcarriesonthedebuggingandrunning,thefunctionofimageprocessing.Furtherdeepentheunderstandingofthevariousoperationsofthesignal,improvingthestudents'comprehensiveproblemanalysisandproblemsolvingskills,laythefoundationforfollow-upcoursesandrelatedresearch目录一、中英文摘要二、正文2.1设计题目2.2试验环境2.3实验目的2.4课程设计原理及说明2.5课程设计的过程及调试2.6设计的收获及结论2.7主要参考文献一、设计题目1、已知信号f1(t)=u(t)-u(t-2)，f2(t)=cos(2πt)用MATLAB绘制f1(t)+f2(t)和f1(t)×f2(t)的波形。2、已知某连续系统的微分方程2y(t)y(t)8y(t)f(t)f(t)etu(t)试用MATLAB绘制出该系统的冲激响应和阶跃响应的波形。3、编程求解正弦函数f(t)sin(t)、取样函数f(t)=Sa(t)、门函数)1()1()(2tututg、双边指数函数te2信号的傅立叶变换式。4、已知某连续系统的的系统函数为H(s)(s23s2)/(8s42s33s25)，试用MATLAB求出该系统的零极点，画出零极点分布图5、设有两个稳定的LTI系统，可分别由下列微分来描述:'()3()3()3''()4'()()''()5()ytytftytytytftft试分别画出它们的系统频率响应的幅值和相位特性曲线。6、信号fl(t)和f2(t)如图所示。tt12111f2(t)f1(t)取t=0：0.005:2.5，计算信号f(t)=f1(t)+f2(t)cos(50t)的值并画出波形；而一可实现的实际系统的H(w)为H(w)=443223410()26.131()3.4142*10()2.6131*10()10jwjwjwjw用freqs函数画出H(w)的幅度和相位曲线。用lsim函数求出信号f1(t)和、f2(t)通过系统H(w)的零状态响应y1(t)和y2(t),并根据理论知识解释所得的结果。二、试验环境计算机、MATLAB软件三、实验目的信号与系统课程设计是在学习完基础上设立的一门实践课。通过本次设计让学生学会使用MATLAB软件，同时学会建立m文件，并且在m文件下建立MATLAB函数和语句以及对其进行调试、运行、函数图像处理保存。更深层次的学习信号与系统以及复变函数与积分变换等课程，掌握信号的的各种处理及计算，同时提高学生的综合分析问题和解决问题的能力，为后续课程和相关专业研究奠定基础四、课程设计原理及说明1、信号的MATLAB表示（1）用MATLAB实现函数f(t)=Sa(t)，并绘制f(t)波形。抽样函数Sa(t)在MATLAB中用sinc函数表示，定义为sinc(t)=sin(pt)/pt其调用形式为y=sinc(t)clearall;t=-5:0.1:5;f=sinc(t);%产生抽样函数plot(t,f);grid;xlabel('t');ylabel('Sa(t)')（2）正弦函数f(t)=Ksin(wt+a)正弦信号和余弦信号在MATLAB中分别用sin和cos表示，其调用形式为clearall;t=-8:0.01:8;k=2;w=1;a=pi/6;f=k*sin(w*t+a);plot(t,f);grid;xlabel('t');ylabel('f(t)')（3）单边指数函数f(t)=-atKe指数信号Keat在MATLAB中可用exp表示，其调用形式为t=0:0.001:10;k=1;a=2;f=k*exp(-a*t);plot(t,f);grid;xlabel('t');ylabel('f(t)')（4）单位冲激信号clearall;t0=0;tf=5;dt=0.01;t1=1;t=[t0:dt:tf];st=length(t);n1=floor((t1-t0)/dt);x1=zeros(1,st);x1(n1)=1/dt;stairs(t,x1),gridon;axis([0,5,0,22])（5）单位阶跃信号利用MATLAB的可视化功能，在work目录下创建函数Heaviside()，可方便地绘出单位阶跃信号波形。需要创建的Heaviside()函数代码如下：functionf=Heaviside(t)f=(t0);调用该函数，并运用plot命令。%利用plot命令绘制单位阶跃信号波形clearall;t=-1:0.01:3;f=Heaviside(t);plot(t,f);axis([-1,3,-0.2,1.2]);set(gcf,'color','w');title('单位阶跃信号')（6）矩形脉冲信号矩形脉冲信号在MATLAB中用rectpuls()函数表示，其调用形式为用以产生一个幅度为1，宽度为width以t为对称的矩形波。clearall;t=0:0.001:4;T=1;ft=rectpuls(t-2*T,2*T);%产生一个以t=t-2T为中心，宽度为2T的矩形脉冲plot(t,ft);gridon;axis([04-0.51.5])（7）三角波脉冲信号三角波脉冲信号在MATLAB中用tripuls()函数表示，其调用形式为用以产生一个最大幅度为1，宽度为width的三角波。函数值的非零范围为（-width/2，width/2）。t=-3:0.001:3;%t决定函数的横坐标范围ft=tripuls(t,4,0.5);%产生一个最大幅度为1，宽度为4，斜度为0.5的三角波plot(t,ft);gridon;axis([-33-0.51.5])（7）用MATLAB实现函数，并绘制x(t)波形。clearall;t=0:0.1:30;x=exp(-0.1*t).*sin(2/3*t);axis([030-11]);plot(t,x);grid;xlabel('time(second)');ylabel('x(t)')2、连续时间系统零状态响应的求解MATLAB工具箱提供了一个用于求解零初始条件下微分方程数值解的函数lsim()，其调用形式为该调用格式对向量t定义的时间范围内，绘制LTI系统的时域波形，同时绘制出系统的激励信号对应的时域波形。其中，t表示计算系统响应的抽样点向量，f是系统输入信号向量，sys是LTI系统模型。在求解微分方程时，LTI的模型sys要借助MATLAB中的tf()函数来获得，其调用形式为其中，b和a分别为微分方程右端和左端各项系数的向量。例如，对于下列微分方程：解：（1）ts=0;te=5;dt=0.01;num=[1];den=[12100];sys=tf(num,den);%调用LTI系统模型的函数t=ts:dt:te;f=10*sin(2*pi*t);y=lsim(sys,f,t);%求零初始条件微分方程数值解plot(t,y);xlabel('t(sec)');ylabel('y(t)');grid（2）clearall;ts=0;te=5;dt=0.01;num=[10];den=[12100];sys=tf(num,den);t=ts:dt:te;y=impulse(sys,t);%求解冲激响应的函数plot(t,y);grid;xlabel('time(sec)');ylabel('h(t)')3、求解冲激响应可用MATLAB工具箱提供的impulse()函数，求解阶跃响应可用step()函数。其调用形式分别为y=step(sys,t)y=impulse(sys,t)其中，t表示计算系统响应的抽样点向量，sys是LTI系统模型。4、连续时间系统的零、极点分析与MATLAB实现利用MATLAB绘制连续系统零、极点分布图如果连续系统的系统函数已知，则可利用MATLAB提供的函数tf、pole、zero、pzmap可以方便地求出系统函数的零、极点，并绘出其零、极点分布图。首先，根据系统函数分子和分母多项式的系数，调用tf函数生成系统函数对象（tf对象）。tf调用格式为其中输入参量num为系统函数分子多项式系数构成的行向量，den为系统函数分母多项式系