北京理工大学信号与系统实验实验报告.doc

《北京理工大学信号与系统实验实验报告.doc》由会员分享，可在线阅读，更多相关《北京理工大学信号与系统实验实验报告.doc（55页珍藏版）》请在三一办公上搜索。

1、信号与系统实验报告 姓名：肖枫 学号：1120111431 班号：05611102 专业：信息对抗技术学院：信息与电子学院 实验1 信号的时域描述与运算一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MATLAB实现方法。 3. 利用MATLAB分析常用信号，加深对信号时域特性的理解。二、实验原理与方法 1. 连续时间信号的MATLAB表示连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。从严格意义上来说，MATLAB并不能处

2、理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： t=0:0.01:10; x=sin(t);利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： x=sin(t); ezplot(X

3、);利用ezplot(x)命令可以绘制上述信号的时域波形 常用的信号产生函数 函数名 功能 函数名 功能heaviside单位阶跃函数rectpuls门函数sin正弦函数tripuls三角脉冲函数cos余弦函数square周期方波sincsinc函数sawtooth周期锯齿波或三角波exp指数函数 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔

4、相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量和采样值向量表示的连续时间信号，其微分可以通过下式求得 其中表示采样间隔。MATLAB中用diff函数来计算差分。 连续时间信号的定积分可以由MATLAB的qud函数实现，调用格式为 quad (function_name,a,b)其中，function_name为被积函数名，a、b为积分区间。 对于符号对象表示的连续时间信号，MATLAB提供了diff函数和quad函数分别用于求微分和积分

5、。 3.离散时间信号的MATLAB表示离散时间信号仅在一些离散时刻有定义。在MATLAB中离散时间信号需要使用两个向量来表示，其中一个向量用于表示离散的时间点，另一个向量表示在这些时间点上的值。例如对于如下时间信号 采用MATLAB可以表示如下： n=-3:4; x=-3 2 -1 2 1 -1 2 3; stem(n,x,filled); xlabel(n); title(x(n);Stem函数用于绘制离散时间信号波形，为了与我们表示离散时间信号的习惯相同，在绘图时一般需要添加filled选项，以绘制实心的杆状图形。上述命令绘制的信号时域波形如图3所示。 4.离散时间信号的时域运算 离散时间

6、信号的相加相乘是将两个信号对应的时间点上的值相加或相乘，可以直接使用算术运算的运算符“+”和“*”来计算。离散时间信号的位移，则可看作是将表示时间的向量平移，而表示对应时间点上的值的向量不变。离散时间信号的反转，则可以看作是将表示时间的向量和表示对应时间点上的值的向量以零点为基准点，一纵轴为对称轴反折，向量的反折可以利用MATLAB的fliplr函数实现。三、实验内容（1） 利用MATLAB绘制下列连续时间信号波形。 MATLAB 程序如下：syms t;x=1-exp(-0.5*t)*heaviside(t);ezplot(x,0,3)波形图如下：MATLAB 程序如下：syms t;x=c

7、os(pi*t)*(heaviside(t)-heaviside(t-2);ezplot(x,0,3)波形图如下： MATLAB 程序如下：syms t;x=sqrt(t2)*(-0.5)*cos(pi*t)*(heaviside(t+2)-heaviside(t-2);ezplot(x,0,3)波形图如下：MATLAB 程序如下：syms t;x=-exp(-t)*sin(2*pi*t)*(heaviside(t)-heaviside(t-3);ezplot(x,0,3)波形图如下：（2） 利用MATLAB绘制下列离散时间信号波形 MATLAB 程序如下：n=-5:5;x=1 1 1 1 1

8、 1 1 1 1 1 1;n1=n+3;subplot(121);stem(n,x,filled);xlabel(n);title(x(n);subplot(122);stem(n1,x,filled);xlabel(n);title(x(n-3);波形图如下：MATLAB 程序如下：n=-5:5;x=power(-1/2,n).*heaviside(n);stem(n,x,filled);xlabel(n);title(x(n);波形图如下：MATLAB 程序如下：n=-3:7;x=n.*(heaviside(n)-heaviside(n-5);stem(n,x,filled);xlabel

9、(n);title(x(n);波形图如下：MATLAB 程序如下：n=-5:8;x=sin(n*pi/2).*heaviside(n);stem(n,x,filled);xlabel(n);title(x(n);波形图如下：（）利用MATLAB生成并绘制连续周期矩形波信号，要求周期为，峰值为，显示三个周期的波形。利用MATLAB 提供的square 函数来生成方波信号MATLAB 程序如下：t=0:0.001:6a=3;p=3;y=a*square(p*t);plot(t,y)axis(t(1)-1 t(end)+1 -(a+1) (a+1)波形图如下：（4） 已知信号，及信号，用MATLAB

10、绘出下列信号的波形。由教材上的波形可知为一个三角波的右半部分。MATLAB 程序如下：syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4);x2=sin(2*pi*t);ezplot(x1+x2,-2,5);波形图如下：MATLAB 程序如下：syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4);x2=sin(2*pi*t);ezplot(x1*x2,-2,5);波形图如下：MATLAB 程序如下：syms t;x1(t)=(-t+4)*(heaviside(t)-heaviside(t-4);x2(t)=sin(2*

11、pi*t);ezplot(x1(-t)+x1(t),-4,5);波形图如下：MATLAB 程序如下：syms t;x1(t)=(-t+4)*(heaviside(t)-heaviside(t-4);x2(t)=sin(2*pi*t);ezplot(x1(t-1)+x2(t-1)*x2(t),-2,5);波形图如下：(5) 已知离散时间信号，用MATLAB绘出的波形。 由教材上的波形可知MATLAB 程序如下：n=-3:3;x=0 1 2 3 3 3 3;stem(n,x,filled);xlabel(n);title(x(n);波形图如下：MATLAB 程序如下：n=-3:3;x=0 1 2

12、3 3 3 3;n1=-fliplr(n);x1=fliplr(x);%subplot(122)stem(n1,x1,filled);xlabel(n);title(x(-n);波形图如下：MATLAB 程序如下：x=0 1 2 3 3 3 3;n1=n-2;x1=fliplr(x);%subplot(122)stem(n1,x,filled);xlabel(n);title(x(n+2);波形图如下：MATLAB 程序如下：n=-3:3;x=0 1 2 3 3 3 3;n1=n+2;x1=fliplr(x);%subplot(122)stem(n1,x,filled);xlabel(n);t

13、itle(x(n-2);波形图如下：（6） 用MATLAB编程绘制下列信号的时域波形，观察信号是否为周期信号？若是周期信号，周期是多少？若不是周期信号，请说明原因。MATLAB 程序如下：t=-5:0.01:5; x=1+cos(pi/4*t-pi/3)+2*cos(pi/2*t-pi/4)+cos(2*pi*t); plot(t,x); xlabel(t); title(x(t);波形图如下：该信号是周期信号，周期TMATLAB 程序如下： t=-5:0.01:5; x=sin(t)+2*sin(pi*t); plot(t,x); xlabel(t); title(x(t);波形图如下：MA

14、TLAB 程序如下： n=-3:7;x=2+3.*sin(2.*n.*pi./3-pi./8);stem(n,x,filled);xlabel(n);title(x(n);波形图如下： 该信号是周期信号，周期T3MATLAB 程序如下：n=-3:7;x=cos(n.*pi./6)+sin(n.*pi./3)+cos(n.*pi./2);stem(n,x,filled);xlabel(n);title(x(n);波形图如下：该信号是周期信号，周期T=12四、 实验心得体会 通过此次实验我掌握了MATLAB的基本使用及其在信号分析上的应用，达到了实验目的，在实验的过程中我体会到了MATLAB功能的

15、强大，提高了我对信号与系统这门课程的兴趣。 实验3 信号频域分析一、 实验目的 1.深入理解信号频谱的概念，掌握信号的频域分析方法。 2.观察典型周期信号和非周期信号的频谱，掌握其频谱特性。二、 实验原理与方法 1.连续周期信号的频谱分析如果周期信号满足狄里赫利条件，就可以展开为傅里叶级数形式，即 式中，表示基波周期，为基波频率，表示任一个基波周期内的积分。式（1）和式（2）定义为周期信号复指数形式的傅里叶级数，系数称为的傅里叶系数。周期信号的傅里叶级数还可以由三角函数的线性组合来表示，即 式（3）中同频率的正弦项和余弦项可以合并，从而得到三角函数形式的傅里叶级数，即 可见，任何满足狄里赫利条

16、件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。一般来说周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中经常采用有限项级数来替代，所选项数越多就越逼近原信号。 2.连续非周期信号的频谱分析对于非周期连续时间信号，吸纳后的傅里叶变换和傅里叶逆变换定义为 式（7）和式（8）把信号的时域特性和频域特性联系起来，确立了非周期信号和频谱之间的关系。采用MATLAB可以方便地求取非周期连续时间信号的傅里叶变换，这里我们介绍常用的集中方法。1） 符号运算法MATLAB的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数，fourier函数和ifourier函数，基

17、本调用格式为 X=fourier(x) X=ifourier(X)默认的时域变量为t，频域变量为w。2） 数值积分法除了采用符号运算的方法外，我们还可以利用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析，quad函数是一个用来计算数值积分的函数。利用quad函数可以计算非周期连续时间信号的频谱。Quad函数的一般调用格式为：y=quad(fun,a,b)y=quad(fun,a,b,TOL,TRACE,p1,p2,)其中fun指定被积函数，可以采用inline命令来创建，也可以通过传递函数句柄的形式来指定，a、b表示定积分的下限和上限，TOL表示允许的相对或绝对积分误差

18、，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“p1,p2，”表示被积函数出时间t之外所需的其他额外输入参数。3）数值近似法我们还可以利用MATLAB的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换可以由式（9）近似计算 当为时限信号，且足够小，式（9）可以演变为 而式（10）中求和部分又可以表示成一个行向量和一个列向量的乘积式（11）可以很方便地利用MATLAB实现。3.离散周期时间信号的频域分析基波周期为N的周期序列可以用N个成谐波关系的复指数序列的加权和表示，即 这里k=表示求和仅需包括一个周期内的N项，周期序列在一个周期

19、内的求和与起点无关。将周期序列表示成式（12）的形式，成为离散傅里叶级数，而系数则称为离散傅里叶系数。离散傅里叶系数可以由式（13）确定。 傅里叶系数也称为的频谱系数，而且可以证明是以N为周期的离散频率序列。这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期N与傅里叶系数的乘积来表示周期离散时间信号的频谱，即 可以利用MATLAB提供的函数fft用来计算，调用格式为 该函数返回一个周期内的值，其中x表示一个周期内的样本值。4.离散非周期时间信号的频域分析非周期序列可以表示成一组复指数序列的连续和 其中 式（16）称为的离散时间傅里叶变换，式（15）和式（16）确立了非周期离

20、散时间信号及其离散时间傅里叶变换之间的关系。是连续频率的函数，称为频谱函数，且是周期的连续频率函数，其周期为。可见，非周期离散时间函数对应于频域中是一个连续的周期的频率函数。对于有限长的时间序列，式（16）可以表示为式（17）可以方便地利用MATLAB实现。三、实验内容（1）已知周期矩形脉冲信号。教材上的波形可知信号由式（3）和式（4）计算得 故的傅里叶级数为 利用MATLAB绘出由前N次谐波合成的信号波形，观察随着N的变化合成信号波形的变化规律；MATLAB程序如下： t=-3:0.01:3;N=input(N=);A=input(A=);T=input(T=);c=input(c=);x=

21、A*c/T;for n=1:Nx=x+2*A/(n*pi)*sin(n*pi*c/T)*cos(2*pi*n*t);endplot(t,x);xlabel(Time(sec)title(N= num2str(N)当N=10，A=1，T=1，c=0.5时，图形如下当N=20，A=1，T=1，c=0.5时，图形如下当N=40，A=1，T=1，c=0.5时，图形如下当N=80，A=1，T=1，c=0.5时，图形如下由述4个波形图可知，随着N的增大即选取的傅里叶级数的项数增加，合成波形越来越接近原来的矩形脉冲信号。利用MATLAB绘出周期矩形脉冲信号的频谱，观察参数T和变化时对频谱波形的影响。由式（1

22、）和式（2）计算得 MATLAB程序如下：N=input(N=);A=input(A=);T=input(T=);i=input(c=);n1=-N:-1;c1=A./n1./pi.*sin(n1.*pi.*i./T);c0=A.*i./T;n2=1:N;c2=A./n2./pi.*sin(n2.*pi.*i./T);cn=c1 c0 c2;n=-N:N;subplot(211);stem(n,abs(cn),filled);xlabel(omega/omega_0);title(Magnitude of ck);subplot(212);stem(n,angle(cn),filled);xl

23、abel(omega/omega_0);title(Phase of ck);当N=20，A=1，T=4，c=0.5时，图形如下当N=20，A=1，T=2，c=0.5时，图形如下当N=20，A=1，T=2，c=1时，图形如下由以上3个频谱波形图可知，频谱波形与的值有关，对于不同的T和，当比值相同时频谱波形图相同；当比值不同时，比值越小，频谱包络形状趋于收敛，过零点越少，谱线越密。思考题解答：1） 将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总

24、跳变值的9%。这种现象称为吉伯斯现象。产生吉伯斯现象的原因是当一个信号通过某一系统时，如果这个信号不是连续时间函数，则由于一般物理系统对信号高频分量都有衰减作用，从而产生吉伯斯现象。2） 由实验结果可知周期性矩形脉冲信号的频谱是离散的，谱线间隔为,即各次谐波仅存在于基频的整数倍上；而且，谱线的长度随着谐波次数的增高趋于收敛。至于频谱收敛规律，以及谐波含量则由信号波形决定。因此，离散型、谐波性和收敛性是周期信号的共同特点。3） 通常把包含主要谐波分量的这段频率范围称为矩形脉冲信号的有效频带宽度，简称为有效频宽。即 ，T为脉冲宽度。故有效频宽与信号的时域宽度成反比，实验结果与这相符。4） 频谱波形

25、与的值有关，对于不同的T和，当比值相同时频谱波形图相同；当比值不同时，比值越小，频谱包络形状趋于收敛，过零点越少，谱线越密。（2） 已知矩形脉冲信号。求该信号的傅里叶变换由教材上的波形可知用MATLAB函数库里的fourier函数来求的傅里叶变换MATLAB程序如下：syms t;A=input(A=);c=input(c=);x=A*heaviside(t+c/2)-heaviside(t-c/2); %x(t)X=fourier(x) collect(X) 当A=1；=1；ans = (2*sin(w/2)/w=Asinc（/2）所以的傅里叶变换为利用MATLAB绘出矩形脉冲信号的频谱，观

26、察矩形脉冲宽度变化时对频谱波形的影响。MATLAB程序如下：syms t;A=input(A=);c=input(c=);x=A*(heaviside(t+c/2)-heaviside(t-c/2);X=fourier(x);ezplot(abs(X),-8*pi,8*pi);grid on;xlabel(omega);ylabel(Magnitude);title(|X(omega)|);当A=1；=1；当A=1；=0.5；当A=2；=0.5；当A=1；=2；由以上4个频谱波形图可知，当矩形脉冲宽度增大时，信号占有频带减小，即信号的占有频带与脉冲宽度成反比。让矩形脉冲的面积始终等于1，改变矩

27、形脉冲宽度，观察矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。MATLAB程序如下：syms t;i=input(i=);A=1/i;x=A.*(heaviside(t+i/2)-heaviside(t-i/2);X=fourier(x);subplot(211);ezplot(x);xlabel(Time(seconds);title(x(t);subplot(212);ezplot(abs(X),-6*pi,6*pi);grid on;xlabel(omega); ylabel(Magnitude);波形图如下：随着脉冲宽度减小，时域波形的幅值增加，信号占有的频带宽度增加。思考题解答

28、：1） 相同点：它们的有效频带宽度都是与脉冲宽度成反比。 不同点：周期矩形脉冲信号的频谱是离散的，而矩形脉冲信号的频谱是连续的2） 矩形脉冲信号的有效频带宽度与时域宽度成反比，当，脉冲面积始终等于1时，其频谱会无限趋近于高度为1的一条直线。（3） 已知周期方波序列。由教材上的波形可知MATLAB程序如下：N1=input(N1=);N=input(N=);n=0:N-1;x=1 1 1 0 0 0 0 1 1;X=fft(x);subplot(211);stem(n,x,filled);xlabel(n);title(x(n);subplot(212);stem(n,X,filled);xla

29、bel(k);title(X(k) , N= num2str(N);波形图如下：思考题解答：） 相同点：周期序列与连续周期信号的频谱都是离散的，且都有收敛性和谐波性。不同点：连续周期信号在一个周期内要用无限多项级数来表示，而周期序列用有限项级数就可以表示。） 随着方波序列占空比的减小，其频谱的谱线越来越密集，谱线高度会增加。（）已知一矩形脉冲序列利用MATLAB绘制周期方波序列的频谱波形，改变矩形脉冲序列的宽度，观察频谱波形的变化趋势。MATLAB程序如下：w = -pi:0.01:pi;N = input(N=);n = -N:N;x = ones(size(n);X = x*exp(-j*

30、n*w);subplot(211);stem(n,x,filled);xlabel(n);title(x(n);subplot(212);plot(w/pi,X);xlabel(omega/pi);title(X(ejomega) , N= num2str(N);波形图如下：思考题解答：随着矩形脉冲序列宽度的增加，其频谱的有效频带宽度减小。其宽度与频谱的有效频带宽度成反比。四、 实验心得体会 通过此次实验我掌握了用MATLAB对连续时间信号和离散时间信号进行傅里叶变换的方法，信号的频域分析方法对于信号处理是十分有用而基本的，对我的知识掌握起到了很大作用。 实验5 连续时间系统的复频域分析一、

31、实验目的1. 掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB实现方法。2. 学习和掌握连续时间系统系统函数的定义及复频域分析方法。3. 掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、 实验原理与方法1. 拉普拉斯变换 连续时间信号的拉普拉斯变换定义为 拉普拉斯反变换定义为 在MATLAB中，可以采用符号数学工具箱的laplace函数和ilaplace函数进行拉氏变换和反拉氏变换。L=laplace(F)符号表达式F的拉氏变换，F中时间变量为t，返回变量为s的结果表达式。L=laplace(F,t)用t替换结果中的变量s。F=ilaplace(L)以s为变量的符号表达式

32、L的拉氏反变换，返回时间变量为t的结果表达式。F=ilaplace(L,x)用x替换结果中的变量t。除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下：当 X (s)为有理分式时，它可以表示为两个多项式之比： 式（3）可以用部分分式法展成一下形式 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。利用 MATLAB 的residue 函数可以将 X (s)展成式（1-2）所示的部分分式展开式，该函数的调用格式为：r,p,k = residue(b,a) 其中b、a 为分子和分母多项式系数向量，r、p、k 分别为上述展开式中的部分分式系数、极点和

33、直项多项式系数。2. 连续时间系统的系统函数 连续时间系统的系统函数是系统单位冲激响应的拉氏变换 此外，连续时间系统的系统函数还可以由系统输入和系统输出信号的拉氏变换之比得到 单位冲激响应反映了系统的固有性质，而从复频域反映了系统的固有性质。由式（6）描述的连续时间系统，其系统函数为s的有理函数 3.连续时间系统的零极点分析系统的零点是指式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统函数的值无穷大。通常将系统函数的零极点绘在s平面上，零点用表示，极点用表示，这样得到的图形称为零极点的分布图。由零极点的定义可知，零点和极点分别指式（7）的分子多项式和分母

34、多项式的根。利用MATLAB求多项式的根可以通过函数roots来实现，该函数的调用格式为：r=roots(c) c为多项式的系数向量，返回值r为多项式的根向量。分别对式（7）的分子多项式和分母多项式求根即可得到零极点。此外，在MATLAB中还提供了更简便的方法来求取零极点和绘制系统函数的零极点分布图，即利用pzmap函数，该函数的调用格式为：pzmap(sys)绘出由系统模型sys描述的系统的零极点分布图。p,z=pzmap(sys) 这种调用方法返回极点和零点，而不绘出零极点分布图。其中sys为系统传函模型，由t命令sys=tf(b,a)实现，b、a为传递函数的分子多项式和分母多项式的系数向

35、量。MATLAB还为用户提供了两个专用函数tf2zp和zp2tf来实现系统传递函数模型和零极点增益模型的转换，其调用格式为： z,p,k=tf2zp(b,a) b,a=zp2tf(z,p,k)其中b、a为传递函数的分子多项式和分母多项式的系数向量，返回值z为零点列向量，p为极点列向量，k为系统函数零极点形式的增益。三、实验内容（1）已知系统的冲激响应,输入信号，是采用复频域的方法求解系统的响应，编写MATLAB程序实现。使用卷积定理求解，先分别求的拉氏变换然后根据式（6）求出输出，最后对进行拉普拉斯反变换即可得到系统的响应。MATLAB程序如下：syms t;h=heaviside(t)-he

36、aviside(t-2);x=heaviside(t);H=laplace(h);X=laplace(x);F=X.*H;f=ilaplace(F);f=ilaplace(F)f = t - heaviside(t - 2)*(t - 2)故系统响应为（2） 已知因果连续时间系统的系统函数分别如下，试采用MATLAB画出其零极点分布图，求解系统的冲激响应并判断系统是否稳定。MATLAB程序及运行结果如下:b=1;a=1 2 2 1;sys=tf(b,a)p,z=pzmap(sys)sys = 1 - s3 + 2 s2 + 2 s + 1 Continuous-time transfer fu

37、nction.p = -1.0000 -0.5000 + 0.8660i -0.5000 - 0.8660iz = Empty matrix: 0-by-1系统没有零点，极点为p=-1，-0.50.866i；另：pzmap(sys)r,p,k=residue(b,a) 得：r = 1.0000 -0.5000 - 0.2887i -0.5000 + 0.2887ip = -1.0000 -0.5000 + 0.8660i -0.5000 - 0.8660ik = 故系统冲激响应为：h(t)=(exp(-1)+(-0.5-0.2887*i).*exp(-0.5+0.866*i)+(-0.5+0.

38、2887*i).*exp(-0.5-0.866*i).*heaviside(t)频率响应：b=1;a=1 2 2 1;H w=freqs(b,a);subplot(211);plot(w,abs(H);xlabel(omega(rad/s);ylabel(Magnitude);title(|X(jomega)|);grid on;subplot(212);plot(w,angle(H);xlabel(omega(rad/s);ylable(Phase);title(phi(omgea);由于该因果系统的所有极点都位于S 平面的左半平面，所以系统是稳定的。MATLAB程序及运行结果如下:b=1

39、0 1;a=1 2 -3 3 3 2;sys=tf(b,a);p,z=pzmap(sys)p = -3.1704 0.9669 + 0.9540i 0.9669 - 0.9540i -0.3817 + 0.4430i -0.3817 - 0.4430iz = 0 + 1.0000i 0 - 1.0000i另：pzmap(sys)则有系统的零点为z=i，极点为p=-3.1704，0.96690.9540i，-0.38170.4430i，r,p,k=residue(b,a)得：r = 0.0769 -0.0300 - 0.0881i -0.0300 + 0.0881i -0.0085 - 0.1436i -0.0085 + 0.1436ip = -3.1704 0.9669 + 0.9540i 0.9669 - 0.9540i -0.3817 + 0.4430i -0.3817 - 0.4430ik = 故系统冲激响应为：h(t)=(0.0769.*exp(-3.1704.*t)+(-0.03-0.0881.*i).*e