武汉大学电气工程学院信号与系统MATLAB仿真报告.doc

信号与系统上机实验报告电气工程学院2008级五班吕业诚2008302540145信号与系统上机实验实验一 信号的产生与运算一、实验目的1、熟悉和掌握常用的用于信号和系统时域仿真分析的MATLAB函数；2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程；基本要求：掌握用MATLAB描述连续时间信号和离散时间信号的方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方法再现各种信号的波形。二、实验原理MATLAB提供许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期矩形波信号等。这些基本信号时信号处理的基础。在MATLAB中，无论是连续时间信号还是离散时间信号，MATLAB都是用一个数字数字序列来表示信号，这个数字序列在MATLAB中叫向量（vector）。通常的情况下，需要与时间变量相对应。如前所述，MATLAB有很多内部数学函数可以用来产生这样的数字序列，例如：sin()、cos()、exp()等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。三、实验内容Q1_1:修改程序Program1_1，并以Q1_2为文件名存盘，产生实指数信号x(t)=e-2t。要求在图形中加上网格线，并使用函数axis（）控制图形的时间范围在02秒之间。然后执行该程序，保存所得程序。close all;clear;clc;t=-2:0.01:2;x=exp(-2*t);plot(t,x);grid on;axis(0,2,0,1);Q1_2:修改程序program1_1，并以Q1_3为文件名存盘，使之能够仿真从键盘上任意输入的一个连续时间信号，并利用该程序仿真信号x(t)=e-2t。clear;close all;clc;t=-2:0.01:2; %定义采样时间x=input(Type in signal x(t) in closed form:); %通过键盘输入plot(t,x)Q1_3:将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编辑器中编好，并分别以文件名为delta和u保存在工作文件夹中以便使用。delta.m:function y=delta(t)dt=0.01;y=(u(t)-u(t-dt)/dt;u.m:% Unit step function阶跃函数function y=u(t)y=(t=0);Q1_4:根据示例程序的编程方法，编写一个MATLAB程序，由给定信号x(t)=e-2tu（t），求信号y(t)=x(1.5t+3),并绘制x(t)和y(t)的图形。clear;close all;clc;syms t;x=sym(exp(-0.5*u(t);y=subs(x,t,1.5*t+3);subplot(2,1,1)ezplot(x,-5,5);subplot(2,1,2);ezplot(y,-5,5);Q1-5:给定一个离散时间信号xn=un-un-8,仿照示例程序Program1_5，编写程序Q1_8,产生xn的左移序列x1n=xn+6和右移序列x2n=xn-6,并在同一个图形窗口的三个子图中分别绘制这三个序列的图形。clear;close all;clc;n=-10:10;x=u(n)-u(n-8);x1=u(n+6)-u(n+6-8);x2=u(n-6)-u(n-6-8);hold onplot(n,x,r);plot(n,x1,g);plot(n,x2);axis(-10,10,0,1.1)Q1_6:总结1.信号与系统分析，就是基于信号的分解，在时域中，信号主要分解成冲击信号的积分2.比较plot，stem，stairs，ezplot命令的区别：plot：折线图stem：二维杆图stairs：阶梯图ezplot：符号函数二维图实验二 连续时间系统的时域分析一、实验目的1、加深对线性时不变系统中零状态响应概念的理解，掌握其求解方法；2、掌握给定连续系统的冲激响应和阶跃响应；3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质；4、掌握利用MATLAB计算卷积的编程方法，并利用所编写的MATLAB程序验证卷积的常用基本性质；5、掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线；6、掌握线性时不变连续系统的时域数学模型用MATLAB描述的方法，掌握卷积运算、线性常系数微分方程的求解编程。二、实验原理1、LTI系统的单位冲激响应和卷积模型给定一个连续时间LTI系统，在系统的初始条件为零时，用单位冲激信号(t)作用系统，此时系统的响应信号称为系统的单位冲激响应，一般用h(t)来表示。需要强调的是，系统的单位冲激响应是在激励信号为(t)时的零状态响应。离散时间LTI系统的单位冲激响应的定义与连续时间LTI系统的单位冲激响应相同，只是离散时间单位冲激函数(n)的定义有所不同。系统输入信号x(t)、xn和输出信号y(t)、yn之间的关系可以用一个数学表达式来描述，这个数学表达式为这个表达式就是LTI系统的卷积模型，它是根据系统的线性和时不变性以及信号可以分解成单位冲激函数经过推理得到的。这个表达式告诉我们一个很重要的结论，那就是：任意LTI系统可以完全由它的单位冲激响应h(t)/hn来确定。2、用线性常系数微分方程描述LTI系统线性常系数微分方程或差分方程是描述LTI系统的另一个时域模型。一个连续时间LTI系统，它的输入信号x(t)输出信号y(t)关系可以用下面的微分方程来表达在MATLAB中，一个LTI系统也可以用系统微分方程的系数来描述。例如一个LTI连续时间系统的微分方程为MATLAB则用两个系统向量num=1和den=1 3 2来描述该系统，其中num和den分别表示系统微分方程右边和左边的系数，按照微分方程运算的降阶排列。MATLAB的内部函数impulse()，step()，initial()，lsim()可以用来计算并绘制连续时间LTI系统的单位冲激响应，单位阶跃响应，零输入响应和任意信号作用于系统的零状态响应。这些函数的用法描述如下： h=impulse(num，den，T ) 和 impulse(num，den，T )s=step(num，den，T ) 和 step(num，den，T )y=lsim(num，den，x，t ) 和 lsim(num，den，x，t ) 函数impulse()和step()用来计算由num和den表示的LTI系统的单位冲激响应和单位阶跃响应，响应的时间范围为0T，其中den和num分别为系统微分方程左右两边的系数向量，T为指定的响应的终点时间。h和s的点数默认值为101点，由此可以计算步长为dt=T/(101-1)。不带返回值的函数如impulse(num，den，T )和step(num，den，T )将直接在屏幕上绘制系统的单位冲激响应和单位阶跃响应曲线。带返回值的函数如lsim(num，den，x，t)和y=lsim(num，den，x，t)，用来计算由num和den表示的LTI系统在输入信号x的作用下的零状态响应。其中t为指定的时间变化范围，x为输入信号，他们的长度应该是相同的。3、卷积的计算卷积的计算通常可按下面五个步骤进行1. 该换两个信号波形图中的横坐标，由t改为T，T变成函数的自变量；2. 把其中一个信号反摺，如把h(t)变成h(-T)；3. 把反摺的信号做移位，移变量是t，这样t是一个参变量。在T坐标系中，t0时图形右移，t0时图形右移，t0时图形左移；（4） 计算两个信号重叠部分的乘积x(T)h(t-T)；（5） 完成相乘后图形的积分。利用MATLAB计算卷积的函数是conv()，其语法为：y=conv(x,h)。其中x和h分别是两个作卷积运算的信号，y为卷积结果。2、在时域中，描述一个连续时间LTI系统的数学模型有： 3、MATLAB是如何表示一个由微分方程描述连续时间LTI系统的？求解连续时间LTI系统的单位冲激响应、单位阶跃响应以及在某一个输入信号作用下的零状态响应的MATLAB函数有哪些？在MATLAB中表示一个由微分方程描述连续时间LTI系统的是：x=Ax=Bu单位冲激响应函数：impulse(num，den，T )单位阶跃响应函数：step(num，den，T )在某一个输入信号作用下的零状态响应函数：lsim(num，den，x，t )其中den和num分别为系统微分方程左右两边的系数向量，T为指定的响应的终点时间，t为指定的时间变化范围，x为输入信号。实验三 连续时间系统的频谱-傅里叶变换一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；3、掌握连续时间傅里叶变换的分析方法及其物理意义；4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质；5、学习掌握利用MATLAB语言编写计算CTFS、CTFT和DTFT仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT、DTFT的若干重要性质。基本要求：掌握并深刻理解傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用MATLAB编程完成相关的傅里叶变换的计算。二、实验原理及方法1、连续时间周期信号的傅里叶级数CTFS分析任何一个周期为T1的正弦周期信号，只要满足狄利克利条件，就可以展开傅里叶级数。其中三角傅里叶级数为：或：其中，称为信号的基本频率，a0，ak和bk分别是信号x(t)的直流分量、余弦分量幅度和正弦分量幅度，、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率k0的函数，绘制出它们与k0之间的图像，称为信号的频谱图。指数形式的傅里叶级数为： 其中，为指数形式的傅里叶级数的系数，按如下公式计算： 指数形式的傅里叶级数告诉我们，如果一个周期信号x(t)，满足狄里克利条件，那么它就可以被看做是由很多不同频率的互为谐波关系的周期复指数信号所组成，其中每一个不同频率的周期复指数信号称为基本频率分量，其复幅度为ak。这里“复幅度”指的是ak通常是复数。2、连续时间信号傅里叶变换-CTFT傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。傅里叶变换和其逆变换定义如下：连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。任意非周期信号，如果满足狄里克利条件，那么，它可以被看做是由无穷多个不同频率的周期复指数信号的线性组合构成的，每个频率所对应的周期复指数信号称为频率分量，其相对幅度对应频率的之值，其相位为对应频率的的相位。3、傅里叶级数的MATLAB计算设周期信号的基本周期为，且满足狄里克利条件，则其傅里叶级数的系数可由下式得到：基本频率为：对周期信号进行分析时，我们往往只需要对其在一个周期内进行分析即可，通常选择主周期。4、用MATLAB实现CTFT的计算MATLAB进行傅里叶变换有两种方法，一种利用符号运算的方法计算，另一种是数值就按，本实验要求采用数值计算的方法来进行傅里叶变换的计算。严格说用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件，即信号是时限信号，也就是当时间| t |大于某个给定时间时其值衰减为零或接近为零，这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的，计算机只能处理有限大小和有限数量的数。采用数值计算算法的理论依据是：若信号为时限信号，当时间间隔T取得足够小时，上式可以演变为： 上式用MATLAB表示为： X=x*exp(j*t*w)*T其中X为信号的傅里叶变换，W为频率，T为时间步长。三、实验内容Q3-1 给程序Program3_1增加适当的语句，并以Q2_2存盘，使之能够计算例题3-1中的周期方波信号的傅里叶级数的系数，并绘制信号的幅度谱和相位谱线图。clear;close all;clc;T=2;dt=0.00001;t=0:dt:2;x1=(mod(t,2)1);w0=2*pi/T;N=13;ak(1)=0;for k=0:N ak(2+k)=(1/T)*x1*exp(-i*k*w0*t)*dt;endabs=abs(ak);phi=angle(ak)/2/pi*360;subplot(2,1,1);stem(0:N+1,abs,.k);xlabel(it幅频特性);subplot(2,1,2);stem(0:N+1,phi,k);xlabel(fontsize14fontname黑体相频特性);Q3-2 分别手工计算的傅里叶级数的系数。t121-2-1写出信号在其主周期内的数学表达式；用MATLAB手工计算的傅里叶级数的系数（系数ak从-10到10共21个系数），写出过程并保存结果。仿照program3_1，编写程序，以计算的傅里叶级数的系数。程序：clear;close all;clc;T=2;dt=0.00001;t=0:dt:2;x1=(mod(t,2)1).*2-1).*(mod(t,1)-0.5)+0.5;w0=2*pi/T;N=13;ak(1)=0;for k=0:N ak(2+k)=(1/T)*x1*exp(-i*k*w0*t)*dt;endabs=abs(ak);phi=angle(ak)/2/pi*360;subplot(2,1,1);stem(0:N+1,abs,.k);xlabel(it幅频特性);subplot(2,1,2);stem(0:N+1,phi,k);xlabel(fontsize14fontname黑体相频特性);实验四 连续时间系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法；3、学习和掌握幅度特性、相位特性的物理意义；4、掌握用MATLAB语言进行系统频率响应特性分析的方法；5、进一步理解信号的抽样及抽样定理；基本要求：掌握LTI连续和离散时间系统的频域数学模型和频域数学模型的MATLAB描述方法，深刻理解LTI系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB计算和绘制LTI系统频率响应特性曲线中的编程。二、实验原理及方法1 用MATLAB计算系统频率响应在本实验中，表示系统的方法仍然是用系统函数分子和分母多项式系数行向量来表示。试验中用到的MATLAB函数如下：H,w=freqs(b,a)：b,a分别为连续时间LTI系统的微分方程右边和左边的系数向量，返回的频率响应在各频率点的样点值（复数）存放在H中，系统默认的样点数目为200点；Hm=abs(H)：求模数，即进行运算，求得系统的幅度频率响应，返回值存于Hm之中。real(H)：求H的实部；imag(H)：求H的虚部；phi=atan(-imag(H)./(real(H)+eps)：求相位频率响应特性，atan()用来计算反正切值；或者phi=angle(H)：求相位频率响应特性；tao=grpdlay(num，den，w)：计算系统的相位频率响应所对应的群延时；计算频率响应的函数fregs()的另一种形式是： H=freg(b，a，w)：在指定的频率范围内计算系统的频率响应特性。在使用这种形式的fregs/freqz函数时，要在前面先指定频率变量w的范围。2、信号的抽样及抽样定理根据傅里叶变换的频率卷积定理，时域两个信号相乘，对应的积的傅里叶变换等于这两个信号的傅里叶变换之间的卷积。所以，已知抽样信号的傅里叶变换为：实际抽样过程，很容易用简单的数学公式来描述。设连续时间信号用表示，抽样周期为，抽样频率为，则已抽样信号的数学表达式为：三、实验内容Q4-1：给范例程序关键语句加注释Q4-2：运行示例程序并保存结果。例一：clear,close,clc;b=1;a=1,3,2;H,w=freqs(b,a);H=H.;Hm=abs(H);phi=angle(H)*180/pi;Hr=real(H);Hi=imag(H);subplot(2,2,1);plot(w,Hm);grid on;title(Magnitude response);xlabel(Frequency in rad/sec);subplot(2,2,2);plot(w,phi);grid on;title(Phase response);xlabel(Frequency in rad/sec);subplot(2,2,3);plot(w,Hr);grid on;title(Real part of frequencey response);xlabel(Frequency in rad/sec);subplot(2,2,4);plot(w,Hi);grid on;title(Imaginary part of frequencey response);xlabel(Frequency in rad/sec);例二：clear;close all;clc;%初始化t=0:0.01:10;Ts=1/4;n=0:Ts:10;x=cos(0.5*pi*t);xn=cos(0.5*pi*n);subplot(2,1,1)plot(t,x);title(A continuous-time signal x(t);xlabel(Time t);subplot(2,1,2)stem(n,xn,.);title(The sampled version xn of x(t);xlabel(Time index n);例三：clear;close all;clc;tmax=4;dt=0.01;t=0:dt:tmax;Ts=1/10;ws=2*pi/Ts;w0=20*pi;dw=0.1;w=-w0:dw:w0;n=0.1:1:tmax/Ts;x=exp(-4*t).*u(t);xn=exp(-4*n*Ts);subplot(2,2,1)plot(t,x);title(A continuous-time signal x(t);xlabel(Time t);axis(0,tmax,0,1);grid on;subplot(2,2,3)stem(n,xn,.);title(The sampled version xn of x(t);xlabel(Time index n);axis(0,tmax/Ts,0,1);grid on;Xa=x*exp(-i*t*w)*dt;X=0;for k=-8:8 X=X+x*exp(-i*t*(w-k*ws)*dt;endsubplot(2,2,2)plot(w,abs(Xa)title(Magnitude spectrum of x(t);xlabel(Frequency in radians/s);grid onaxis(-60,60,0,1.8*max(abs(Xa);subplot(2,2,4)plot(w,abs(X)title(Magnitude spectrum of xn);xlabel(Frequency in radians/s);grid onaxis(-60,60,0,1.8*max(abs(Xa);L1R12C1100mF800mHQ4-3：已知RLC二阶低通滤波器如图所示，其中L=0.8H，C=0.1F，R=2。试用MATLAB绘制出该频率响应。程序：clear;close all;clc;L=0.8;C=0.1;R=2;A=L*C,L/R,1;B=1;H,w=freqs(B,A);H=H.;Hm=abs(H);phi=angle(H)*180/pi;Hr=real(H);Hi=imag(H);subplot(2,2,1);plot(w,Hm);grid on;title(Magnitude response);xlabel(Frequency in rad/sec);subplot(2,2,2);plot(w,phi);grid on;title(Phase response);xlabel(Frequency in rad/sec);subplot(2,2,3);plot(w,Hr);grid on;title(Real part of frequencey response);xlabel(Frequency in rad/sec);subplot(2,2,4);plot(w,Hi);grid on;title(Imaginary part of frequencey response);xlabel(Frequency in rad/sec);实验五 连续时间系统的复频域分析一、实验目的1、掌握拉普拉斯变换的物理意义、基本性质及应用；2、掌握用拉普拉斯变换求解连续时间LTI系统的时域响应；3、掌握系统函数的概念，掌握系统函数的零、极点分布与系统的稳定性、时域特性等之间的互相关系；4、掌握用MATLAB对系统进行变换域分析的常用函数及编程方法；基本要求：掌握拉普拉斯变换及其基本性质，掌握应用拉普拉斯变换求解系统的微分方程，能够自己编写程序完成对系统时域响应的求解。掌握并理解系统函数的概念，掌握系统函数零极点与系统时域和频域特性之间的关系，能够编写程序完成对系统的一些主要特性如稳定性、因果性等的分析。二、实验原理及方法1、连续时间LTI系统的复频域描述拉普拉斯变换主要用于系统分析。描述系统分析。描述系统的另一种数学模型就是建立在拉普拉斯变换基础上的“系统函数”H(s)：系统函数H(s)的实质就是单位冲激响应h(t)的拉普拉斯变换。因此，系统函数也可以定义为：所以，系统函数的一些特点是和系统的时域响应的特点对应的。在教材中，我们求解系统函数的方法，常用的是根据描述系统的线性系数微分方程，经过拉氏变换之后得到系统函数。假设描述一个连续时间LTI系统的线性常系数微分方程为：对上式两边做拉普拉斯变换，则有：2、系统函数的零极点分布图系统函数的零极点分布图能够直观地表示零点和极点在s平面上的位置，从而比较容易分析系统函数的收敛域和稳定性。对于一个连续时间LTI系统，它的全部特性包括稳定性、因果性和它具有何种滤波特性等完全由它的零极点在s平面上的位置所决定。3、拉普拉斯变换与傅里叶变换之间的关系根据课本知识可知，拉普拉斯变换与傅里叶变换之间的关系可表述为：傅里叶变换是信号在虚轴上的拉普拉斯变换，也可以用下面的数学表达式表示：上式表明，给定一个信号，如果它的拉普拉斯变换存在的话，它的傅里叶变换不一定存在，只有当它的拉普拉斯变换的收敛域包括了整个虚轴，则表明傅里叶变换时存在的。4、系统函数的极点分布与系统的稳定性和因果性之间的关系一个稳定的LTI系统，它的单位冲激响应满足绝对可积条件，即 同时，我们还应该记得一个信号的傅里叶变换的存在条件就是这个信号满足绝对可积条件，所以，如果系统是稳定的话，那么，该系统的频率响应也必然是存在的。又根据傅里叶变换与拉普拉斯变换之间的关系，可进一步推出，稳定的系统，其系统函数的收敛域必然包括虚轴。稳定的因果系统，其系统函数的全部极点一定位于s平面的左半平面。所以，对于一个给定的LTI系统，它的稳定性、因果性完全能够从它的零极点分布图上直观地看出。三、实验内容Q5-1：将绘制零极点图的扩展函数文件splane以splane为文件名存盘。splane.m:function splane(num,den)p=roots(den);q=roots(num);p=p;q=q;x=max(abs(p,q);x=x+1;y=x;hold onplot(-x,x,0,0,:);plot(0,0,-y,y,:);plot(real(p),imag(p),x,markersize,8);plot(real(q),imag(q),o,markersize,8);title(Zero-pole plot);xlabel(Real part);ylabel(Imaginal Part);axis(-x,x,-y,y);Q5-2：编写程序Q5_2，能够接受从键盘输入的系统函数的分子多项式系数向量，并绘制常系统的零极点图、系统的单位冲激响应、系统的幅度频率响应和相位频率响应的图形。程序：clear;close all;clc;T=10;a=input(Please input the den columnn);b=input(Please input the num columnn);%计算响应H,w=freqs(b,a);H=H.;Hm=abs(H);phi=angle(H)*180/pi;Hr=real(H);Hi=imag(H);%输出subplot(2,2,1)splane(b,a); subplot(2,2,2);plot(w,Hm);grid on;title(Magnitude response);xlabel(Frequency in rad/sec); subplot(2,2,4);plot(w,phi);grid on;title(Phase response);xlabel(Frequency in rad/sec); subplot(2,2,3)impulse(b,a,T)Q5-3：执行程序编写Q5_3，输入系统函数的分子分母系数向量，绘制所得到的图形：实验六 离散时间信号与系统一、实验目的1. 掌握用MTALAB软件产生离散时间信号的方法；2. 进一步理解离散时间信号的常见运算方法；3. 掌握求离散时间系统冲激响应的方法；二、实验原理1、离散时间信号的产生单位抽样序列单位阶跃序列矩形序列单位斜坡序列正选序列指数序列2、离散时间系统的冲激响应N阶差分方程：MATLAB中，filter子程序可用来在给定输入和差分方程系数是求差分方程的数值解，调用格式为：Y=filter(B,A,X)其中，；X为输入信号序列，注意必须保证系数不为零。三、实验内容Q6-1：随意产生两个序列信号，然后进行相加、乘和卷积运算。随意产生的两个序列：clear;close all;clc;x1=1 5 -4 2 5 -1 5;x2=1 2 3 4 5 6 7;n=1:length(x1);subplot(2,2,1)hold onstem(n,x1,x,markersize,7);stem(n,x2,.,markersize,5);subplot(2,2,2)stem(n,x1+x2,markersize,5);subplot(2,2,3)stem(n,x1.*x2,d,markersize,5);subplot(2,2,4)y=conv(x1,x2);stem(1:length(y),y,*,markersize,5);Q6-2：利用filter命令求下面系统的冲激响应y(n)-0.7(n-1)-0.6y(n-2)+y(n-3)=x(n)+0.5x(n-1)clear;close all;clc;A=1,-0.7,-0.6,1;B=1,0.5;n=-5:50;X=n=0;y=filter(B,A,X);stem(n,y,.,markersize,6);title(冲激响应);grid on实验七 离散时间系统的时域、Z域分析一、实验目的（1）加深对线性时不变离散系统中零状态响应概念的理解，掌握其求解方法；（2）深刻理解卷积和运算，掌握求离散序列卷积和的计算方法；（3）掌握求解给定离散系统的单位序列响应和单位阶跃序列响应的方法；（4）加深理解和掌握求序列信号Z变换和逆Z变换的方法；（5）加深理解和掌握离散系统的系统函数零点、极点分布与系统时域特性、系统稳定性的性质。 二、实验原理（1）线性时不变离散时间系统用常系数线性差分方程进行描述：其中，和分别表示系统的输入和输出，N=max(n，m)是差分方程的阶数。在已知差分方程的初始状态以及输入的条件下，可以通过编程由下式迭代算出系统的输出： 系统的零状态响应就是在系统初始状态为零条件下微分方程的解。在零状态下，MATLAB控制系统工具箱提供了一个filter函数，可以计算差分方程描述的系统的响应，其调用形式为：其中，分别是系统差分方程左、右端的系数向量，表示输入向量，表示输出向量。注意，输出序列的长度与输入序列的长度相等。（2）离散系统的冲激响应、阶跃响应分别是输入信号为和所对应的零状态响应。MATLAB控制系统工具箱专门提供了两个函数求解离散系统的冲激响应和阶跃响应。冲激响应：其中的表示系统的单位序列响应，分别是系统差分方程的左右端的系数向量，K表示输出序列的时间范围。阶跃响应：其中的表示系统的单位阶跃序列响应，和的含义与上相同，N表示输出序列的长度。（3）卷积是信号与系统中一个最基本、也是最重要的概念之一。在时域中，对于LTI离散时间系统，其零状态响应等于输入信号与系统冲激响应的卷积；而利用卷积定理，这种关系又对应频域中的乘积。MATLAB信号处理工具箱提供了一个计算两个离散序列卷积和的函数。设向量和代表待卷积的两个序列，则就是与卷积后得到的新序列。（4）如有序列（k为整数），z为复变量，则函数称为序列的双边Z变换。如果上式的求和只在k的非负值域进行，则称为序列的单边Z变换。MATLAB的符号数学工具箱提供了计算Z正变换的函数和计算逆Z变换的函数，其调用形式为：上式中，右端的f和F分别表示为时域表示式和Z域表示式的符号表示，可利用函数sym来实现，其调用形式为。式中的A为待分析的表示式的字符串，S为符号化的数字或变量。(5) 线性时不变离散系统可用其Z域的系统函数H(Z)表示，去通常具有如下有理分式的形式：为了能够从系统的Z域表示方便地得到其时域表示式，可将H(z)展开为部分分式和的形式，再对其求逆Z变换。MATLAB的信号处理工具箱提供了对H(z)进行部分分式展开的函数，其调用形式为：式中的B和A分别是H(z)的分子多项式和分母多项式的系数向量，r为部分分式的向量，p为极点向量，k为多项式的系数向量。（6）通过系统函数的表达式，可以方便地求出系统函数的零点和极点。系统函数的零点和极点的位置对于系统的时域特性和频域特性有重要影响。位于Z平面的单位圆上和单位圆外的极点将使得系统不稳定，系统函数的零点将使得系统的幅频响应在该频率点附近出现极小值，而对应的极点将使得系统的幅频响应在该频率点附近出现极大值。三、实验内容Q7-1:运行上述示例程序，对关键语句标注。例一：clear;close all;clc;R=100;d=rand(1,R)-0.5;%产生随机噪声k=0:R-1;s=2*k.*(0.9.k);%原始信号f=s+d;%加噪信号figure;subplot(1,2,1);plot(k,d,r-,k,s,b:,k,f,k-);xlabel(time index k);legend(dk,sk,fk);m=5;b=ones(m,1)/m;%产生右边向量系数a=1;y=filter(b,a,f);subplot(1,2,2)plot(k,s,b:,k,y,r-);xlabel(time index k);legend(sk,yk);示例二：clear;close all;clc;k=0:10;a=1,3,2;b=1;h=impz(b,a,k);g=stepz(b,a,length(k);figure;subplot(2,2,1)stem(k,h);grid ontitle(单位取样响应的近似值);subplot(2,2,3)stem(k,g);title(单位阶跃响应 的近似值);grid onhk=-(-1).k+2*(-2).k;gk=1/6-(-1).k/2+4*(-2).k/3;subplot(2,2,2)stem(k,hk);grid ontitle(单位取样响应的理论值)subplot(2,2,4)stem(k,gk);title(单位阶跃响应的理论值);grid on例三：clear;close all;clc;k1=-1:3;x=1,2,3,4,5;k2=0:4;h=1,1,1,1,1;figure(1);subplot(1,2,1)stem(k1,x);grid onxlabel(输入序列xk);subplot(1,2,2)stem(k2,h)grid onxlabel(单位序列响应hk)y=conv(x,h);k=k1(1)+k2(1):k1(length(k1)+k2(length(k2);figure;stem(k,y);grid onxlabel(输出响应yk);例四：clear;