离散LSI系统分析

1、 信号与系统 分析实验报告 实验项目名称：离散线性时不变系统分析； 连续时间系统分析所属课程名称：信号与系统实验教程实 验 类 型 ：验证型指 导 教 师 ：实 验 日 期 ：2013.06.04班 级 ：学 号 ：姓 名 ：离散线性时不变系统分析一、实验目的1. 掌握离散线性时不变系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB求解方法。2.掌握离散线性时不变系统的频域分析方法；3.掌握离散线性时不变系统的复频域分析方法；4. 掌握离散线性时不变系统的零极点分布与系统特性的关系。二、实验原理及方法1. 离散线性时不变系统的时域分析描述一个N阶线性时不变离散时间系统的数学模型是线

2、性常系统差分方程，N阶线性时不变离散系统的差分方程一般形式为 （2.1）也可用系统函数来表示 （2.2）系统函数反映了系统响应和激励间的关系。一旦上式中，的数据确定了，系统的性质也就确定了。特别注意必须进行归一化处理，即。对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位序列或单位阶跃序列的线性叠加，把这些单元激励信号分别加于系统求其响应，然后把这些响应叠加，即可得到复杂信号作用于系统的零状态响应。因此，求解系统的单位序列响应和单位阶跃响应尤为重要。由图2-1可以看出一个离散LSI系统响应与激励的关系。图2-1 离散LSI系统响应与激励的关系(1) 单位序列响应（单位响应）单位响应

3、是指离散线性时不变系统在单位序列激励下的零状态响应，因此满足线性常系数差分方程(2.1)及零初始状态，即 ， (2.3)按照定义，它也可表示为 (2.4)对于离散线性时不变系统，若其输入信号为，单位响应为，则其零状态响应为 (2.5)可见，能够刻画和表征系统的固有特性，与何种激励无关。一旦知道了系统的单位响应，就可求得系统对任何输入信号所产生的零状态响应。MATLAB提供了专门用于求离散系统冲激响应的函数impz()，其调用格式有h,n=impz(b,a)求解离散系统的单位响应，其中，；h,n=impz(b,a,N) 求解离散系统的单位响应，采样点数由N确定，；impz(b,a) ：在当前窗口

4、，用stem(n,h)绘出图形。（2）单位阶跃响应 单位阶跃响应是指离散离散时不变系统在单位阶跃序列激励下的零状态响应，它可以表示为 (2.6)上式表明，离散线性时不变系统的单位阶跃响应是单位响应的累加和，系统的单位阶跃响应和系统的单位响应之间有着确定的关系，因此，单位阶跃响应也能完全刻画和表征一个线性时不变系统。MATLAB提供了专门用于求离散系统单位阶跃响应的函数stepz( )，其调用格式有s,n=stepz(b,a) ：求解离散系统的单位阶跃响应，其中，；s,n=stepz(b,a,N) ：求解离散系统的单位阶跃响应，采样点数由N确定，；stepz(b,a) ：在当前窗口，用stem(

5、n,s)绘出图形。（3）任意激励下的零状态响应已经知道，离散线性时不变系统可用常系数线性差分方程(2.1)式来描述，Matlab提供的函数dlsim( )能对上述差分方程描述的离散线性时不变系统的响应进行仿真，该函数不仅能绘制指定时间范围内的系统响应波形图，而且还能求出系统响应的数值解。其调用格式有dlsim(b,a, x) ：求解输入序列为x的零状态响应需要特别强调的是，Matlab总是把由分子和分母多项式表示的任何系统都当作是因果系统。所以，利用impz (b,a)，stepz(b,a)，dlsim(b,a,x)函数求得的响应总是因果信号。同时，卷积和也是线性时不变系统求解零状态响应的重要

6、工具之一。假设系统的输入信号为，单位响应为，则系统的零状态响应可由（2.5）式求解。Matlab提供了专门用于求离散系统卷积和的函数conv( )，其调用格式有y=conv(x,h) ：求解序列x，h的卷积和，若序列x的长度为n1，序列h的长度为n2，卷积和y的长度为n1+n2-1。这一点需要特别注意，否则，作图时容易造成横纵坐标长度不匹配。（4）带初始状态的任意激励下的全响应任意激励下的离散线性时不变系统的全响应为零输入响应和零状态响应之和，表示为 （2.7）Matlab提供了用于求离散系统全响应的函数filter( )，其调用格式有y=filter( b,a,x) ：求解零状态响应；y=f

7、ilter( b,a,x,zi) ：求解初始条件为zi的系统的全响应，zi向量的长度为max(length(a),length(b)-1，返回值为系统的全响应。z = filtic(b,a,y,x)：将初始状态转换为初始条件，其中，；z = filtic(b,a,):将初始状态转换为初始条件,2 离散线性时不变系统的复频域（Z域）分析（1）利用Z变换解差分方程在前面图2-1中表示了离散系统的响应与激励的关系，由图可知，系统的响应既可以用时域的方法求解，也可以用Z域的方法求解。当已知系统输入序列的Z变换，系统函数时，系统响应序列的Z变换可由求出。Matlab提供了用于求序列Z变换和Z反变换的函数

8、，其调用格式有X=ztrans(x)：求无限长序列x的Z变换，返回Z变换的表达式，注意这里x，X都是符号表达式；x=iztrans(X)：求X（z）的Z反变换x(n)，返回Z反变换的表达式，注意这里x，X都是符号表达式；r,p,c=residuez(b,a)：把b(z)/a(z)展开成部分分式；b,a=residuez(r,p,c):根据部分分式的r、p、c数组，返回有理多项式。（2）系统的零极点分布与系统因果性和稳定性的关系因果系统的单位响应一定满足当时，=0，那么其系统函数的收敛域一定包含点，即点不是极点，极点分布在某个圆的圆内，收敛域是圆外区域。系统稳定要求 ，对照z变换定义，系统稳定要

9、求收敛域包含单位圆。如果系统因果且稳定，收敛域包含点和单位圆，那么收敛域可表示为： (2.8)MATLAB提供了用于求系统零极点的函数，其调用格式有roots()：利用多项式求根函数来确定系统函数的零极点位置；roots(a)：求极点位置，a为系统函数分母多项式所构成的系数向量；roots(b)：求零点位置，b为系统函数分子多项式所构成的系数向量；zplane(b,a)：绘制由行向量b和a构成的系统函数的零极点分布图；zplane(z,p)：绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。（3）系统的零极点分布与系统频率响应的关系将式(2.2)因式分解，得到 (2.9)式中，是

10、的零点，是其极点。参数影响频率响应的幅度大小，影响系统特性的是零点和极点的分布。下面采用几何方法研究系统零极点分布对系统频率特性的影响。将式(2.9)的分子、分母同乘以，得到： (2.10)假设系统稳定，将代入上式，得到频率响应 (2.11)设，由式(2.11)得到 (2.12)在平面上，用一根由零点指向单位圆（）上任一点B的向量表示，同样用一根由极点指向单位圆（）上任一点B的向量来表示，和分别称为零点矢量和极点矢量，用极坐标表示为： ，。将和表示式代入式(2.12)，得到 (2.13) (2.14)系统或者信号的频率特性由式(2.13)和式(2.14)确定。按照式(2.13)，知道零极点的分

11、布后，可以很容易地确定零极点位置对系统特性的影响。当点转到极点附近时，极点矢量长度最短，因而幅度特性可能出现峰值，且极点愈靠近单位圆，极点矢量长度愈短，峰值愈高愈尖锐。如果极点在单位圆上，则幅度特性为，系统不稳定。对于零点，情况相反，当点转到零点附近时，零点矢量长度变短，幅度特性将出现谷值，且零点愈靠近单位圆，谷值愈接近零。当零点在单位圆上时，谷值为零。综上所述，极点位置主要影响频响的峰值位置及尖锐程度，零点位置主要影响频响的谷值位置及形状。Matlab提供了专门用于求离散系统频率响应的函数freqz(),其调用格式如下：H,w = freqz(b,a,n)：返回数字系统的n点频率值（复数），

12、这n个点均匀地分布在0，上，系统默认的采样点数目为512点；H,f = freqz(b,a,n,Fs)：用于对在0，Fs/2上等间隔采样n点，采样点频率及相应的频响值分别存放在f和H中。H = freqz(b,a,w)：用于对在0，2上进行采样，采样频率点由w指定。H = freqz(b,a,f,Fs)：用于对在0，Fs上进行采样，采样频率点由f指定。freqz(b,a,n)：用于在当前图形窗口中绘制幅频和相频特性曲线。下面介绍几个实用的函数：mag=abs(H):求解系统的绝对幅频响应；db=20*log10(mag+eps)/max(mag)：求解系统的相对幅频响应；pha=angle(H

13、): 求解系统的相频响应；grd=grpdelay(b,a,w): 求解系统的群延迟。三、实验内容1.设某LTI系统的，输入，求系统输出:主程序如下：clear;clf;clc;n=-5:20; %设定一个n的观察范围h=delta(n-5);x=stepseq(0,-5,20)-stepseq(5,-5,20);y,ny=conv_m(x,n,h,n)subplot(2,2,1);stem(n,x);title('xn');xlabel('n');axis(-5,20,0,1.2);grid on;subplot(2,2,2);stem(n,h);title(

14、'hn');xlabel('n');axis(-5,20,0,1.2);grid on;subplot(2,2,3);stem(ny,y);title('yn');xlabel('n');axis(-5,20,0,1.2);grid on;delta.mfunction y=delta(x)y=(x=0);stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2)if nargin=3%nargin(number of arguments input)；“ =”表示不相等 disp('Us

15、age: Y=stepseq(n0,n1,n2)'); return;else if(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=n2')%n0<n1或者n0>n2或者n1>n2endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程序如下：functiony,ny=conv_m(x,nx,h,nh)%Modified convolution routine for signal processing%y,ny=conv_m(x

16、,nx,h,nh)%y=convolution result%ny=support of y%x=first signal on support nx%nx=support of x%h=second signal on support nh%nh=support of hif nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1); %nys beginingnye=nx(length(x)+nh(length(h); %nys endny=nyb:nye; %ny仅仅为了计算一下结果y对应

17、的横坐标范围y=conv(x,h); %Convolution and polynomial multiplication;conv为MATLAB的固有函数%conv_m函数比conv函数多用了用于表示横坐标范围的ny、nx、nh,因为这里假定二维坐标范围nx和nh可能从负数开始，所以要重新计算y的横坐标范围ny；如果nx和nh均为0:N，则可以直接得出ny为0:2N(正如conv函数中那样）；在conv_m和conv函数中，x和h的横坐标范围都可以实不相等的思考题：若,试写出与的关系，并对MATLAB的仿真结果。解：因为函数是延时器，所以yn=xn* =x-。2.设某线性时不变系统的hn=

18、输入信号为xn=求输出：a ： b：（a）主程序如下：n=-5:20;u1=stepseq(0,-5,20);u2=stepseq(6,-5,20);%u1=un；u2=un-6%input xnx=u1-u2;%impulse response hnh=n.*x;subplot(3,1,1);stem(n,x);axis(-5 20 0 2);title('Input Sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-5 20 0 6);title('Inpulse Response'

19、);ylabel('hn'); %output responsey1,ny=conv_m(x,n,h,n);%conv_m为自定义求卷积的函数subplot(3,1,3);stem(ny,y1);title('Output Sequence');xlabel('n');ylabel('y_1n');（b）主程序如下：n=-10:20;u1=stepseq(0,-10,20);u2=stepseq(6,-10,20);%u1=un；u2=un-6%input xnx=u1-u2;u3=stepseq(-5,-10,20);u4=st

20、epseq(1,-10,20);%u3=un+5；u4=un-1x1=u3-u4;%impulse rsponse hn+5h=(n+5).*x1;subplot(3,1,1);stem(n,x);axis(-10 20 0 2);title('Input sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-10 20 0 6);title('Inpulse Response');ylabel('hn+5');%output responsey2,ny=conv_m(x,n

21、,h,n);subplot(3,1,3);stem(ny,y2);title('Output Sequence');xlabel('n');ylabel('y_2n');stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2&

22、#39;)endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程序如下：functiony,cy=conv_m(x,nx,h,nh)if nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1);nye=nx(length(x)+nh(length(x);ny=nyb:nye;y=conv(x,h);3.设,输入，求系统输出:n=-5:50;u1=stepseq(0,-5,50);u2=stepseq(10,-5,50);%u1=un;u2=un-10%input xnx=u

23、1-u2;%impluse response hnh=(0.9).n).*u1;subplot(3,1,1);stem(n,x);axis(-5 20 0 2);title('Input Sequence');ylabel('xn');subplot(3,1,2);stem(n,h);axis(-5 20 0 6);title('Inpulse Response');ylabel('hn');%output responsey,ny=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y);title(&

24、#39;Output Sequence');xlabel('n');ylabel('yn');stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2')endn=n1:n2;x=(n-n0)>=0;conv_m.m的源程

25、序如下：functiony,cy=conv_m(x,nx,h,nh)if nargin=4 disp('Usage:Y=conv_m(x,nx,h,nh)'); return;end;nyb=nx(1)+nh(1);nye=nx(length(x)+nh(length(x);ny=nyb:nye;y=conv(x,h);思考题：在此问题中为无限长序列，在计算卷积和时应如何处理？解：取尽可能多的包含能量分布集中的具有代表性的点4某离散时间线性时不变系统的差分方程为：a.计算时的系统冲激响应 b.计算时的系统阶跃响应a=1,-0.8;b=1;n=-10:50;x=impseq(0,

26、-10,50);h=filter(b,a,x);%MATLAB固有函数filter表示系统方程ay=bx；x是冲击信号，则h就是冲击响应subplot(2,1,1);stem(n,h);title('Inpulse Response');xlabel('n');ylabel('hn');x=stepseq(0,-10,50);s=filter(b,a,x); %x是阶跃信号，则s就是阶跃响应subplot(2,1,2);stem(n,s);title('Step Response');xlabel('n');yla

27、bel('sn');impseq.m的源程序如下：functionx,n=impseq(n0,n1,n2)%Generates xn=delta(n-n0)；n1<=n<=n2,n1<=n0<=n2%x,n=impseq(n0,n1,n2)；if nargin=3 disp('Usage:Y=impseq(n0,n1,n2)');return;elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=n2')endn

28、=n1:n2;x=(n-n0)=0;stepseq.m的源程序如下：function x,n=stepseq(n0,n1,n2);if nargin =3 disp('Usage:Y=stepseq(n0,n1,n2)');elseif(n0<n1)|(n0>n2)|(n1>n2) error('arguments must satisfy n1<=n0<=2')endn=n1:n2;x=(n-n0)>=0; 连续时间系统分析1.设有两个稳定的LTI系统，分别可由下列微分方程来描述：；。请分别画出它们的系统频率响应的幅值和相位特性曲线。lab41a.m的源程序如下：a=1 3; % y's coeff. in the system equationb=3; % x's coeff.in the system equationfreqs(b,a); % freqs为MATLAB中求连续时间系统拉普拉斯变换域（s域）的频率响应的固有函数Lab41b.m的源程序如下：a=3 4 1; % 系数矩阵b=1 0 5; % 系数矩阵freqs(b,a); % freqs为MATLAB中求连续时间系统拉普拉斯变换域（s域）的频率响应的固有函数思考题：求出系统a的频率响应表达式；根