信号与系统课程设计(MATLAB)

1、信号与系统课程设计开设方案信号与系统课程设计是一门实践环节必修课。其目的是培养学生正确的设计思想，理论联系实际的科学态度，严肃认真、实事求是的科学态度和勇于探索的创新精神。培养学生综合运用所学信号与系统及数字信号处理的知识、分析和解决工程技术问题的能力，为毕业设计打下基础。一、课程基本要求课程设计对象：电子信息工程专业本科班课程设计学时：24二、时间安排开设时间任务目标第8周星期日指导老师给学生布置课程设计的任务及要求。学生完成组队，选题工作。第9周星期日学生根据设计任务及要求，查阅资料，熟悉相关的理论知识。第10周星期日学生设计并确定方案第11周星期日学生编写并调试程序第12周星期日学生撰写

2、课程设计报告第13周星期日课程设计答辩，评分三、设计题目、任务及要求 课程设计可参考以下6个题目，学生可在此范围内选择。设计题目一、连续信号卷积积分的MATLAB实现设计题目二、连续信号频域分析的MATLAB实现设计题目三、基于MATLAB的系统零极点及频率响应分析设计题目四、LTI系统零状态响应的MATLAB仿真设计题目五、在MATLAB环境下实现连续信号复频域分析设计题目六、基于MATLAB的离散信号与系统的Z域分析设计题目七信号的产生与时域运算设计题目八 信号频谱分析设计题目九 线性连续时间系统的分析设计题目十 离散时间信号与系统（详见：课程设计指导书）四、设计报告及书写内容要求 课程设

3、计任务完成后，每位同学必须独立书写一份课程设计报告，注意：不得抄袭他人的报告（或给他人抄袭），一旦发现，成绩为零分。课程设计报告的内容应包括以下四个部分：（1）设计题目，要求。（2）软件设计调试过程及仿真结果。（3）总结：包括课程设计过程中的学习体会与收获。五、考核方式指导老师负责验收设计结果，并结合学生的工作态度、实际动手能力、创新精神和设计报告等进行综合考评，并按优秀、良好、中等、及格和不及格五个等级给出每位同学的课程设计成绩。具体考核标准包含以下几个部分：（1）平时出勤（占10%）。（2）程序设计调试（占30%）。（3）设计报告（占30%）。（4）答辩效果（占30%）。 课程设计指导书设

4、计题目一  连续信号卷积积分的MATLAB实现一、 设计目的1  熟悉卷积积分的定义和性质2  了解卷积积分在系统分析中的应用3  熟悉用MATLAB实现卷积的方法二、 基本原理信号的卷积是数学上的一种积分运算，两个信号的卷积定义为：信号的卷积运算在系统分析中主要用于求解系统的零状态响应。一般情况，卷积积分的运算比较困难，但在MATLAB中则变得十分简单，MATLAB中是利用conv函数来实现卷积的。conv功能：实现二个函数 和 的卷积。格式： 说明：      表示二个函数，    

5、;   表示两个函数的卷积结果。例题：已知两信号                                    求卷积           &

6、#160;    MATLAB程序如下：t1=1:0.01:2;f1=ones(size(t1);               %高度为一的门函数，时间从t=1到t=2t2=2:0.01:3;f2=ones(size(t2);               %高度为一的门函数，时间从t=2到

7、t=3g=conv(f1,f2);                  %对f1和f2进行卷积t3=3:0.01:5;subplot(3,1,1);plot(t1,f1);grid; %画f1的波形subplot(3,1,2),plot(t2,f2);grid  %画f2的波形subplot(3,1,3),plot(t3,g);grid   %画g的波形三、设计内容1  已知两信号

8、0;                        求卷积     比较此题与例题。2 已知两信号                      

9、;             求卷积           四、 设计报告要求1  简述设计目的和原理。2  理论上计算信号的卷积积分。3  记录程序运行结果。4  收获与建议。设计题目二  连续信号频域分析的MATLAB实现一、设计目的1  熟悉傅里叶变换的性质2  熟悉常见信号的傅里叶变换3  了解傅里叶变换

10、的MATLAB实现方法二、 基本原理傅里叶变换是信号分析的最重要的内容之一。从已知信号f (t)求出相应的频谱函数F(j)的数学表示为：f (t)的傅里叶变换存在的充分条件是f (t)在无限区间内绝对可积，即f (t)满足下式：    但上式并非傅里叶变换存在的必要条件。在引入广义函数概念之后，使一些不满足绝对可积条件的函数也能进行傅里叶变换。    傅里叶反变换的定义为：       在这一部分的学习中，大家都体会到了这种数学运算的麻烦。在MATLAB语言中有专门对信号进行

11、正反傅里叶变换的语句，使得傅里叶变换很容易在MATLAB中实现。在MATLAB中实现傅里叶变换的方法有两种，一种是利用MATLAB中的Symbolic Math Toolbox提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换，另一种是傅里叶变换的数值计算实现法。下面分别介绍这两种实现方法的原理。1、直接调用专用函数法   在MATLAB中实现傅里叶变换的函数为：l         F=fourier( f )     对f(t)进行傅里叶变换，其结果为

12、F(w)l         F=fourier(f,v)   对f(t)进行傅里叶变换，其结果为F(v)l         F=fourier( f,u,v )  对f(u)进行傅里叶变换，其结果为F(v)傅里叶反变换l         f=ifourier( F )     对F(w)进

13、行傅里叶反变换，其结果为f(x)l         f=ifourier(F,U)   对F(w)进行傅里叶反变换，其结果为f(u)l         f=ifourier( F,v,u )  对F(v)进行傅里叶反变换，其结果为f(u)    由于MATLAB中函数类型非常丰富，要想了解函数的意义和用法，可以用mhelp命令。如在命令窗口键入：mhelp fourier回车，

14、则会得到fourier的意义和用法。 注意： 在调用函数fourier( )及ifourier( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。对fourier( )中的f及ifourier( )中的F也要用符号定义符sym将其说明为符号表达式。具体方法参见第一部分第四章第三节采用fourier( )及fourier( )得到的返回函数，仍然为符号表达式。在对其作图时要用ezplot( )函数，而不能用plot()函数。fourier( )及fourier( )函数的应用有很多局限性，如果在返回函数中含有()等函数

15、，则ezplot( )函数也无法作出图来。另外，在用fourier( )函数对某些信号进行变换时，其返回函数如果包含一些不能直接表达的式子，则此时当然也就无法作图了。这是fourier( )函数的一个局限。另一个局限是在很多场合，尽管原时间信号f(t)是连续的，但却不能表示成符号表达式，此时只能应用下面介绍的数值计算法来进行傅氏变换了，当然，大多数情况下，用数值计算法所求的频谱函数只是一种近似值。例：求门函数f(t)=(t+1)-(t-1)的傅里叶变换，并画出幅度频谱图MATLAB程序如下：syms t w       &#

16、160;                       %定义两个符号变量t,wGt=sym('Heaviside(t+1)-Heaviside(t-1)');     %产生门宽为2的门函数Fw=fourier(Gt,t,w);        

17、              %对门函数作傅氏变换求FFw=maple('convert',Fw,'piecewise');        %数据类型转换FFP=abs(FFw);                

18、;         %求振幅频谱ezplot(FFP,-10*pi 10*pi);grid;            %绘制函数图形，并加网格axis(-10*pi 10*pi 0 2.2)                 %限定坐标轴范围运行结果：例

19、：求函数的傅里叶反变换f(t)    MATLAB程序如下：syms t w                             %定义两个符号变量t,wFw=sym('1/(1+w2)');      

20、60;          %定义频谱函数ft=ifourier(Fw,w,t);                   %对频谱函数进行傅氏反变换运行结果：ft = 1/2*exp(-t)*Heaviside(t)+1/2*exp(t)*Heaviside(-t)2、傅里叶变换的数值计算实现法严格说来，如果不使用symbolic工

21、具箱，是不能分析连续时间信号的。采用数值计算方法实现连续时间信号的傅里叶变换，实质上只是借助于MATLAB的强大数值计算功能，特别是其强大的矩阵运算能力而进行的一种近似计算。傅里叶变换的数值计算实现法的原理如下：对于连续时间信号f(t)，其傅里叶变换为：其中为取样间隔，如果f(t)是时限信号，或者当大于某个给定值时，f(t)的值已经衰减得很厉害，可以近似地看成是时限信号，则上式中的n取值就是有限的，假定为N，有： 若对频率变量进行取样，得：   通常取：，其中是要取的频率范围，或信号的频带宽度。采用MATLAB实现上式时，其要点是要生成f(t)的N个样本值f（n）的向量，以

22、及向量，两向量的内积（即两矩阵的乘积），结果即完成上式的傅里叶变换的数值计算。注意：时间取样间隔的确定，其依据是必须小于奈奎斯特（Nyquist）取样间隔。如果f(t)不是严格的带限信号，则可以根据实际计算的精度要求来确定一个适当的频率为信号的带宽。       例：用数值计算法实现上面门函数f(t)=(t+1)-(t-1)的傅里叶变换，并画出幅度频谱图    分析： 该信号的频谱为，其第一个过零点频率为，一般将此频率认为是信号的带宽。但考虑到的形状（为抽样函数），假如将精度提高到该值的50倍，即取，则据

23、此确定的Nyquist取样间隔为：MATLAB程序如下：R=0.02;                                %取样间隔=0.02t=-2:R:2;          

24、;                     % t为从-2到2，间隔为0.02的行向量%有201个样本点ft=zeros(1,50),ones(1,101),zeros(1,50);      % 产生f(t)的样值矩阵（即f(t)的样本值%组成的行向量）W1=10*pi;      &#

25、160;                       %取要计算的频率范围M=500; k=0:M; w=k*W1/M;           %频域采样数为M, w为频率正半轴的采样点Fw=ft*exp(-j*t'*w)*R;    

26、;            %求傅氏变换             FRw=abs(Fw);                       %取振幅W=-f

27、liplr(w),w(2:501) ;              %形成负半轴和正半轴的2M+1个频率点WFW=fliplr(FRw),FRw(2:501);         %形成对应于2M+1个频率点的值Subplot(2,1,1) ; plot(t,ft) ;grid;          %画

28、出原时间函数f(t)的波形，并加网格xlabel('t') ; ylabel('f(t)');                %坐标轴标注title('f(t)=u(t+1)-u(t-1)');               %文本标注subplot(2,1,2)

29、 ; plot(W,FW) ;grid;        %画出振幅频谱的波形,并加网格xlabel ('W') ; ylabel ('F(W)');            %坐标轴标注title('f(t)的振幅频谱图');            &#

30、160;  %文本标注运行结果如下：三、 设计内容1、 编程实现下列信号的幅度频谱    a.  求出f1(t)=(2t+1)-(2t-1)的频谱函数F1(j)，请将它与上面门宽为2的门函数f(t)=(t+1)-(t-1)的频谱进行比较，观察两者的特点，说明两者的关系。b.三角脉冲      f2(t)=         c.单边指数信号f3(t)=d.高斯信号    

31、0; f4(t)=         2、 利用ifourier( ) 函数求下列频谱函数的傅氏反变换a.  b.  四、 基本要求1、熟悉常见信号的频谱2、了解MATLAB语言中的函数：fourier , ifourier , abs ,syms ,ezplot。五、设计报告要求1、  傅里叶变换的原理及其数值计算实现法的理论依据2、  记录设计的波形图，写出程序清单3、  收获与建议设计题目三、基于MATLAB的系统零极点及频率响应分析一、  设计目的1 

32、; 掌握系统函数零极点的定义2  零极点与频率响应的关系3  极点与系统稳定性的关系4  状态方程与系统函数的关系5  在MATLAB中实现系统函数与状态方程间的转换二、  基本原理1  原理描述连续系统的系统函数H(s)的一般表示形式为：其对应的零极点形式的系统函数为： 共有n个极点：p1，p2，pn和m个零点：z1，z2，zm。把零极点画在S平面中得到的图称为零极点图，人们可以通过零极点分布判断系统的特性。当系统的极点处在S的左半平面时系统稳定；处在虚轴上的单阶极点系统稳定；处在S的右半平面的极点及处在虚轴上的高阶极点，系统是不稳

33、定的。描述系统除了可以用系统函数和零极图以外，还可以用状态方程。对应上述用系统函数H (s) 描述的系统，其状态方程可用相变量状态方程和对角线变量状态方程描述，形式分别为相变量状态方程：输入方程为 对角线变量方程 输出方程 矩阵中的p为系统函数的极点，k为部分分式展开中的系数，即          上述状态方程和输出方程均可表示为               

34、;             A.B.C.D分别表示对应的矩阵，上述二种表示中D=0。系统在频域中的特性可以用频域中的系统函数表示              H(j)是复函数，可表示为称为幅频特性， 称为相频特性。MATLAB语言提供了系统函数，零极点和状态方程之间的相互转换语句，也提供了得到系统频率特性的语句：tf2zp：从系统函数的一般形式求出其

35、零点和极点。zp2tf：从零极点求出系统函数的一般式。ss2zp：从状态方程式求系统的零极点。zp2ss：从零极点求系统的状态方程。freqs：由H(s)的一般形式求其幅频特性和相频特性。2例题     已知系统函数 ，求其零极点图。MATLAB程序如下：num = 1  0.5  2；           %分子系数，按降幂顺序排列。den = 1  0.4  1；    &

36、#160;        %分母系数，按降幂顺序排列。z，p = tf2zp（num，den）；     %求零点z和极点pzplane (z,p)                    %作出零极点图运行结果如下：     已知系统和状态方程和输出方程

37、60; 求其系统的零极点。MATLAB程序如下：         A = 1，0；1，-3；         B = 1，0；         C = - ，1；         D = 0       

38、60; z，p =  ss2zp (A,B,C,D)       %求出零极点         zplane (z，p)                  %画出零极点图     已知系统的传递函数为 ，求其频率特性。MATLAB程序如下：

39、num = 0.2  0.3  1；           den = 1  0.4  1； w =logspace (-1，1)；   %频率范围freqs(num，den，w)    %画出频率响应曲线运行结果如下：三、  设计内容1  已知下列系统函数H (s)或状态方程，求其零极点，并画出零极点图。     状态方程：输出方程：y = 4 

40、5  1X2  已知下列系统函数H (s)，求其频率特性。 3.       已知系统函数H (s)，求其频率特性和零极图。四、   设计报告要求1  简述设计目的及基本原理2  计算系统的零极点并与设计结果进行比较3  记录频率特性曲线4  收获与建议设计题目四   LTI系统零状态响应的MATLAB仿真一、 设计目的1、  熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2、  熟悉离散时间系统的单位冲激

41、响应、阶跃响应的意义及求解方法3、  熟悉连续（离散）时间系统在任意信号激励下响应的求解方法4、  熟悉应用MATLAB实现求解系统响应的方法二、  基本原理1、连续时间系统对于连续的LTI系统，当系统输入为e(t)，输出为r(t)，则输入与输出之间满足如下的线性常系数微分方程：当系统输入为单位冲激信号(t)时产生的零状态响应称为系统的单位冲激响应，用h(t)表示。若输入为单位阶跃信号(t)时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。系统的单位冲激响应h(t)包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关

42、。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应h(t)对我们进行连续系统的分析具有非常重要的意义。在MATLAB中有专门用于求解连续系统冲激响应和阶跃响应， 并绘制其时域波形的函数impulse( ) 和step( )。如果系统输入为e(t)，冲激响应为h(t)，系统的零状态响应为r(t)，则有：若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。在MATLAB中，应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信

43、号作用下的响应。lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。以上各函数的调用格式如下： impulse( ) 函数    函数impulse( )将绘制出由向量a和b所表示的连续系统在指定时间范围内的单位冲激响应h(t)的时域波形图，并能求出指定时间范围内冲激响应的数值解。l         impulse(b,a)     以默认方式绘出由向量a和b所定义的连续系统的冲激

44、响应的时域波形。l         impulse(b,a ,t0)        绘出由向量a和b所定义的连续系统在0 t0时间范围内冲激响应的时域波形。l         impulse(b,a,t1:p:t2)      绘出由向量a和b所定义的连续系统在t1 t2时间范围内，并且以时间间隔p均匀取样的冲激响应的

45、时域波形。l         y=impulse(b,a,t1:p:t2)     只求出由向量a和b所定义的连续系统在t1 t2时间范围内，并且以时间间隔p均匀取样的冲激响应的数值解，但不绘出其相应波形。    step( ) 函数函数step( )将绘制出由向量a和b所表示的连续系统的阶跃响应，在指定的时间范围内的波形图，并且求出数值解。和impulse( )函数一样，step( )也有如下四种调用格式：step( b,a)step(b,a,

46、t0)step(b,a,t1:p:t2)y=step(b,a,t1:p:t2) 上述调用格式的功能和impulse( )函数完全相同，所不同只是所绘制（求解）的是系统的阶跃响应g(t)，而不是冲激响应h(t)。     lsim( )函数根据系统有无初始状态，lsim( )函数有如下两种调用格式：系统无初态时，调用lsim( )函数可求出系统的零状态响应，其格式如下：l         lsim(b,a,x,t)    绘出由向量a和b所定义的连续系统在输

47、入为x和t所定义的信号时，系统零状态响应的时域仿真波形，且时间范围与输入信号相同。其中x和t是表示输入信号的行向量，t为表示输入信号时间范围的向量，x则是输入信号对应于向量t所定义的时间点上的取样值。l         y=lsim(b,a,x,t)   与前面的impulse 和step函数类似，该调用格式并不绘制出系统的零状态响应曲线，而只是求出与向量t定义的时间范围相一致的系统零状态响应的数值解。系统有初始状态时，调用lsim( )函数可求出系统的全响应，格式如下：l  

48、;       lsim(A,B,C,D,e,t,X0)     绘出由系数矩阵A,B,C,D所定义的连续时间系统在输入为e和t所定义的信号时，系统输出函数的全响应的时域仿真波形。t为表示输入信号时间范围的向量，e则是输入信号e(t)对应于向量t所定义的时间点上的取样值，X0表示系统状态变量X=x1,x2,.xn'在t=0时刻的初值。l         Y,X= lsim(A,B,C,D,e,t,X0)

49、60; 不绘出全响应波形，而只是求出与向量t定义的时间范围相一致的系统输出向量Y的全响应以及状态变量X的数值解。显然，函数lsim( )对系统响应进行仿真的效果取决于向量t的时间间隔的密集程度，t的取样时间间隔越小则响应曲线越光滑，仿真效果也越好。说明：（1）当系统有初始状态时，若使用lsim( )函数求系统的全响应，就要使用系统的状态空间描述法，即首先要根据系统给定的方式，写出描述系统的状态方程和输出方程。假如系统原来给定的是微分方程或系统函数，则可用相变量法或对角线变量等方法写出系统的状态方程和输出方程。其转换原理如前面设计四所述。（2）显然利用lsim( )函数不仅可以分析单输入单输出系

50、统，还可以分析复杂的多输入多输出系统。例题1： 若某连续系统的输入为e(t)，输出为r(t)，系统的微分方程为：求该系统的单位冲激响应h(t)及其单位阶跃响应g(t)。若 求出系统的零状态响应r(t)    求冲激响应及阶跃响应的MATLAB程序：a=1  5  6;b=3  2;subplot(2,1,1), impulse(b,a,4)subplot(2,1,2), step(b,a,4)运行结果如下：    求零状态响应的MATLAB程序：a=1  5  6;b=3 

51、2;p1=0.01;            %定义取样时间间隔为0.01t1=0:p1:5;           %定义时间范围x1=exp(-2*t1);       %定义输入信号lsim(b,a,x1,t1),       %对取样间隔为0.01时系统响应进行仿真

52、hold on;             %保持图形窗口以便能在同一窗口中绘制多条曲线p2=0.5;             %定义取样间隔为0.5t2=0:p2:5;           %定义时间范围x2=exp(-2*t2);  &#

53、160;    %定义输入信号lsim(b,a,x2,t2), hold off   %对取样间隔为0.5时系统响应进行仿真并解除保持运行结果如下：  例题2  已知一个过阻尼二阶系统的状态方程和输出方程分别为：                r(t)=0  1X(t)  若系统初始状态为X(0)=4  -5T   

54、求系统在作用下的全响应。求全响应程序如下：A=0  1 ; -2  -3 ;B=0  2'C=0  1;D=0;X0=4  -5'                       %定义系统初始状态t=0: 0.01:10;        

55、;               E=3*exp(-4*t).*ones(size(t)'        %定义系统激励信号r , x=lsim(A,B,C,D,E,t,X0);          %求出系统全响应的数值解plot(t,r)    

56、60;                      %绘制系统全响应波形运行结果如下:2、离散时间系统LTI离散系统中，其输入和输出的关系由差分方程描述：         （前向差分方程）        （后向差分方程）当系统的输入为单位序列

57、(k)时产生的零状态响应称为系统的单位函数响应，用h(k)表示。当输入为(k)时产生的零状态响应称为系统的单位阶跃应，记为：g(k)，如下图所示。  如果系统输入为e(k)，冲激响应为h(k)，系统的零状态响应为y(k)，则有：与连续系统的单位冲激响应h(t)相类似，离散系统的单位函数响应h(k)也包含了系统的固有特性，与输入序列无关。我们只要知道了系统的单位函数响应，即可求得系统在不同激励信号作用下产生的响应。因此，求解系统的单位函数响应h(k)对我们进行离散系统的分析也同样具有非常重要的意义。MATLAB中为用户提供了专门用于求解离散系统单位函数响应， 并绘制其时域波形的函数im

58、pz( )。同样也提供了求离散系统响应的专用函数filter( )，该函数能求出由差分方程所描述的离散系统在指定时间范围内的输入序列作用时，产生的响应序列的数值解。当系统初值不为零时，可以使用dlsim( )函数求出离散系统的全响应，其调用方法与前面连续系统的lsim( )函数相似。另外，求解离散系统阶跃响应可以通过如下两种方法实现：一种是直接调用专用函数dstep( ),其调用方法与求解连续系统阶跃响应的专用函数step( )的调用方法相似；另一种方法是利用求解离散系统零状态响应的专用函数filter( )，只要将其中的激励信号看成是单位阶跃信号(k)即可。函数的调用格式分别如下：.impz

59、( )函数l  impz(b,a)     以默认方式绘出由向量a和b所定义的离散系统单位函数响应的时域波形。l   impz(b,a,n)     绘出由向量a和b所定义的离散系统在0 n （n必须为整数）的离散时间范围内单位函数响应的时域波形。l impz(b,a,n1:n2)      绘出由向量a和b所定义的离散系统在n1 n2 （n1、n2必须为整数）的离散时间范围内单位函数响应的时域波形。l 

60、y=impz(b,a,n1:n2)     求出由向量a和b所定义的离散系统在n1 n2 （n1、n2必须为整数）的离散时间范围内单位函数响应的数值解，但不绘出波形。 filter( ) 函数l filter(b,a,x)     其中a和b与前面相同，x是包含输入序列非零样值点的的行向量。此命令将求出系统在与x的取样时间点相同的输出序列样值。例题：已知描述离散系统的差分方程为：          且已知系统输入序列

61、为 求出系统的单位函数响应h(k)在-3 10离散时间范围内响应波形。  求出系统零状态响应在0 15区间上的样值； 并画出输入序列的时域波形以及系统零状态响应的波形   求系统的单位函数响应的MATLAB程序：      a=1,-0.25,0.5;        b=1,1,0;                im

62、pz(b,a,-3:10)           %绘出单位函数响应在-3 10区间上的波形      运行结果如下：求零状态响应的MATLAB程序：a=1,-0.25,0.5;b=1,1,0k=0:15;              %定义输入序列取值范围x=(1/2).k;    

63、;       %定义输入序列表达式y=filter(b,a,x)         %求解零状态响应样值subplot(2,1,1)stem(k,x)             %绘制输入序列的波形    subplot(2,1,2)stem(k,y)     

64、        %绘制零状态响应的波形运行结果如下：y =  Columns 1 through 5     1.0000    1.7500    0.6875   -0.3281   -0.2383  Columns 6 through 10     0.1982    0.2156   -0.

65、0218   -0.1015   -0.0086  Columns 11 through 15     0.0515    0.0187   -0.0204   -0.0141    0.0069  Column 16     0.0088三、 设计内容1、已知描述系统的微分方程和激励信号e(t) 分别如下，试用解析方法求系统的单位冲激响应h(t)和零状态响应r(t)，并用MA

66、TLAB绘出系统单位冲激响应和系统零状态响应的波形，验证结果是否相同。                                               

67、60;如下图所示的电路中，已知（），（H），且两电感上初始电流分别为（A）， (A)，如果以电阻上电压作为系统输出，请求出系统在激励  (V)作用下的全响应。2、请用MATLAB分别求出下列差分方程所描述的离散系统，在020时间范围内的单位函数响应、阶跃响应和系统零状态响应的数值解，并绘出其波形。另外，请将理论值与MATLAB仿真结果在对应点上的值作一比较，并说出两者的区别和产生误差的原因。 y(k)+2y(k-1)+ y(k-2)=e(k)            

68、      y(k+2)-0.7y(+1)+0.1 y(k)=7e(k+2)-2e(k+1)              一个带通滤波器可由下列差份方程描述：y(k)+0.81y(k-2)=e(k)- e(k-2)  其中 e(k)为系统输入， y(k)为系统输出。请求出当激励为e(k)=10+10cos(k/2)+10cos(k)(k)   时滤波器的稳态输出。四、 基本要求 

69、   1、熟悉系统响应的求解方法   2、了解MATLAB语言中关于系统分析的各个函数如：impulse、step、lsim、impz、filter等函数的调用方法：五、 设计报告要求1、理论上计算出系统的单位冲激响应/单位函数响应、阶跃响应、零状态响应、全响应的表达式，并写出解题过程。2、记录仿真结果（包括数据和波形）。      3、写出程序清单。     4、收获与建议设计题目五  在MATLAB环境下实现连续信号复频域分析一、设计目的1、

70、熟悉拉普拉斯变换的原理及性质2、熟悉常见信号的拉氏变换3、了解正/反拉氏变换的MATLAB实现方法4、了解利用MATLAB绘制三维曲面图的方法5、了解信号的零极点分布对信号拉氏变换曲面图的影响6、了解连续信号的拉氏变换与傅氏变换的关系二、 基本原理拉普拉斯变换是分析连续时间信号的重要手段。对于当t 时信号的幅值不衰减的时间信号，即在f(t)不满足绝对可积的条件时，其傅里叶变换可能不存在，但此时可以用拉氏变换法来分析它们。连续时间信号f(t)的单边拉普拉斯变换F(s)的定义为：拉氏反变换的定义为：         显然，上式中

71、F(s)是复变量s的复变函数，为了便于理解和分析F(s)随s的变化规律，我们将F(s)写成模及相位的形式：其中， 为复信号F(s)的模，而 为F(s)的相位。由于复变量s=+j，如果以为横坐标（实轴），j为纵坐标（虚轴），这样，复变量s就成为一个复平面，我们称之为s平面。从三维几何空间的角度来看， 和 分别对应着复平面上的两个曲面，如果绘出它们的三维曲面图，则就可以直观地分析连续信号的拉氏变换F(s)随复变量s的变化情况，在MATLAB语言中有专门对信号进行正反拉氏变换的函数，并且利用 MATLAB的三维绘图功能很容易画出漂亮的三维曲面图。   在MATLAB中实现拉氏变换

72、的函数为：l         F=laplace( f )       对f(t)进行拉氏变换，其结果为F(s)l         F=laplace (f,v)   对f(t)进行拉氏变换，其结果为F(v)l         F=laplace ( f,u,v)  &

73、#160; 对f(u)进行拉氏变换，其结果为F(v)拉氏反变换l         f=ilaplace ( F )      对F(s)进行拉氏反变换，其结果为f(t)l         f=ilaplace(F,u)   对F(w)进行拉氏反变换，其结果为f(u)l       

74、0; f=ilaplace(F,v,u )    对F(v)进行拉氏反变换，其结果为f(u) 注意： 在调用函数laplace( )及ilaplace( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。对laplace( )中的f及ilaplace( )中的F也要用符号定义符sym将其说明为符号表达式。具体方法参见第一部分第四章第三节。例：求出连续时间信号 的拉氏变换式，并画出图形求函数拉氏变换程序如下：syms t s      

75、60;              %定义符号变量ft=sym('sin(t)*Heaviside(t)');    %定义时间函数f(t)的表达式Fs=laplace(ft)                %求f(t)的拉氏变换式F(s)运行结果：Fs = 1/(s2+1)

76、绘制拉氏变换三维曲面图的方法有2种：方法一syms x y ss=x+i*y;                    %产生复变量sFFs=1/(s2+1);               %将F(s)表示成复变函数形式FFss=abs(FFs);  

77、             %求出F(s)的模ezmesh(FFss);                %画出拉氏变换的网格曲面图ezsurf(FFss);              

78、0;  %画出带阴影效果的三维曲面图colormap(hsv);               %设置图形中多条曲线的颜色顺序 方法二figure(2)                    %打开另一个图形窗口x1=-5: 0.1:5;  

79、               %设置s平面的横坐标范围y1=-5: 0.1: 5;                %设置s平面的纵坐标范围x,y=meshgrid(x1,y1);        %产生矩阵s=x+i*y; 

80、0;                  %产生矩阵s来表示所绘制曲面图的复平面区域，%其中矩阵s包含了复平面-6<<6，-6<j<6范围内%以间隔0.01取样的所有样点fs=1./(s.*s+1);               %计算拉氏变换在复平面上的样点值ffs

81、=abs(fs);                 %求幅值mesh(x,y,ffs);               %绘制拉氏变换的三维网格曲面图surf(x,y,ffs);          &#

82、160;     %绘制带阴影效果的三维曲面图axis(-5,5,-5,5,0,8);          %设置坐标显示范围colormap(hsv);              %设置图形中多条曲线的颜色顺序说明：从拉普拉斯变换的三维曲面图中可以看出，曲面图上有象山峰一样突出的尖峰，这些峰值点在s平面的对应点就是信号拉氏变换的极点位置。

83、而曲面图上的谷点则对应着拉氏变换的零点位置。因此，信号拉氏变换的零极点位置决定了其曲面图上峰点和谷点位置。 例：求出函数 的拉氏反变换式    MATLAB程序如下： syms t s                     %定义符号变量Fs =sym('1/(1+s2)');     &#

84、160;    %定义F(s)的表达式ft=ilaplace(Fs)               %求F(s)的拉氏反变换式f(t)运行结果：       ft= sin(t)注意: 在MATLAB中，求拉氏反变换的函数ilaplace()，在默认情况下是指拉氏右变换，其运行结果是单边函数。如例中的运行结果为ft= sin(t)，实际上是指ft= sin(t)e(t)

85、三、 设计内容1、  求出下列函数的拉氏变换式，并用MATLAB绘制拉氏变换在s平面的三维曲面图 2、已知信号的拉氏变换如下，请用MATLAB画出其三维曲面图，观察其图形特点，说出函数零极点位置与其对应曲面图的关系，并且求出它们所对应的原时间函数f (t)，       3、已知连续时间信号 ，请分别求出该信号的拉氏变换 及其傅里叶变换 ，并用MATLAB绘出 的曲面图及振幅频谱 的波形，观察 的曲面图在虚轴上的刨面图，并将它与信号的振幅频谱曲线进行比较，分析两者的对应关系。四、  基本要求1、  熟悉信号的拉氏变换2、&#

86、160; 了解MATLAB中的有关函数的调用方法五、  设计报告要求1、  写出拉氏正反变换的原理2、  理论上计算出信号的拉氏正/反变换表达式，并写出解题过程3、  写出相关的程序清单4、  记录设计波形5、 获与建议设计题目六  基于MATLAB的离散信号与系统的Z域分析一、设计目的1、熟悉离散信号Z变换的原理及性质2、熟悉常见信号的Z变换3、了解正/反Z变换的MATLAB实现方法4、了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5、了解利用MATLAB实现离散系统的频率特性分析的方法二、基本原理1、正/

87、反Z变换Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号为：理想抽样信号的双边拉普拉斯变换Fd (s)为：若令 ， ，  那么的双边拉普拉斯变换Fd (s)为：则离散信号f(k)的Z变换定义为：  从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号的拉氏变换Fd (s)之间存在以下关系：               

88、 同理，可以推出离散信号f(k)的Z变换F(z)和它对应的理想抽样信号的傅里叶变换之间的关系为：            如果已知信号的Z变换F(z)，要求出所对应的原离散序列f(k)，就需要进行反Z变换，反Z变换的定义为：   其中，C为包围的所有极点的闭合积分路线。在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans( ) 和itrans( )。其调用格式分别如下：l       &

89、#160; F=ztrans( f )      对f(n)进行Z变换，其结果为F(z)l         F=ztrans(f,v)   对f(n)进行Z变换，其结果为F(v)l         F=ztrans(f,u,v)    对f(u)进行Z变换，其结果为F(v)l         f=itrans ( F )      对F(z)进行Z反变换，其结果为f(n