信号与系统实验指导书(matlab软件仿真)

1、信号与系统实验指导书（MATLAB仿真） # 目录TOC o 1-5 h z HYPERLINK l bookmark4实验一MATLAB基本应用2 HYPERLINK l bookmark22实验二信号的时域表示7 HYPERLINK l bookmark28实验三连续信号卷积11 HYPERLINK l bookmark52实验四典型周期信号的频谱表示18 HYPERLINK l bookmark78实验五傅立叶变换性质研究23 HYPERLINK l bookmark82实验六离散信号分析26 HYPERLINK l bookmark88实验七离散系统的Z域分析29 HYPERLINK

2、l bookmark198Matlab相关符号及函数说明37 实验一MATLAB基本应用头验目的：学习mATLAB的基本用法，了解MATLAB的目录结构和基本功能以及MATLAB在信号与系统中的应用。二头验内容：例一已知x的取值范围，画出y=sin(x)的图型。参考程序：x=0:0.05:4*pi;y=sin(x);plot(y)例二计算y二sin(冗/5)+4cos(冗/4)例三已知z取值范围，x=sin(z);y=cos(z);画三维图形。z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)xlabel(x)ylabel(y)zlabel(z)例四已

3、知x的取值范围，用subplot函数绘图。参考程序：x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(2,2,1),plot(x,y1),title(sin(x)subplot(2,2,2),plot(x,y2),title(1.5*cos(x)subplot(2,2,3),plot(x,y3),title(sin(2*x)subplot(2,2,4),plot(x,y4),title(5*cos(2*x)连续信号的MATLAB表示1、指数信号：指数信号Aeat在MATLAB中可用exp函数表示，其调用形式

4、为:y=A*exp(a*t)(例取A=1,a=-0.4)参考程序：A=l;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);gridon;注：gridon是一个函数，表示在画图的时候添加网格线。2、正弦信号：正弦信号Acos(wOt+)和Asin(wOt+)分别由函数cos和sin表示，其调用形式为：A*cos(wOt+phi);A*sin(wOt+phi)(例取A=1，w0=2冗，=冗/6)参考程序：A=1;w0=2*pi;phi=pi/6;t=0:0.001:8;ft=A*sin(w0*t+phi);plot(t,ft);gridon;3、抽样函数：抽

5、样函数Sa(t)在MATLAB中用sinc函数表示，其定义为：sinc(t)=sin(冗t)/(冗t)其调用形式为：y=sinc(t)参考程序：t=-3*pi:pi/100:3*pi;ft=sinc(t/pi);plot(t,ft);gridon;4、矩形脉冲信号:在MATLAB中用rectpuls函数来表示,其调用形式为:y=rectpuls(t,width),用以产生一个幅值为1,宽度为width，相对于t=0点左右对称的矩形波信号，该函数的横坐标范围由向量t决定，是以t=0为中心向左右各展开width/2的范围，width的默认值为1。例：以t=2T(即t-2,T=0)为对称中心的矩形脉

6、冲信号的MATLAB源程序如下：(取T=1)t=0:0.001:4;T=1;ft=rectpuls(t-2*T，3*T);plot(t，ft);gridon;axis(-0.54-0.11.5);5、周期性矩形波(方波)信号在MATLAB中用square函数来表示，其调用形式为：y=square(t,DUTY)用以产生一个周期为2冗、幅值为1的周期性方波信号，其中的DUTY参数表示占空比，即在信号的一个周期中正值所占的百分比。例如频率为30Hz的周期性方波信号的MATLAB参考程序如下：t=-0.0625:0.0001:0.0625;y=square(2*pi*30*t,75);plot(t,

7、y);axis(-0.06250.0625-1.51.5);gridon;三、实验练习：1、编写代码实现如下图矩形脉冲信号（利用rectpuls函数）t=-2:0.001:2;T=1;ft=rectpuls(t,2*T);plot(t,ft);gridon;axis(-22-0.51.5);2、编写代码实现一个频率50Hz,占空比为50%的周期矩形波（方波）信号（利用square函数）。t=-0.065:0.0001:0.065;y=square(2*pi*50*t,50);plot(t,y);axis(-0.06250.0625-1.51.5);gridon 实验二信号的时域表示实验目的：利

8、用MATLAB实现信号的时域表示以及图形表示题目一：连续信号的MATLAB描述设计要求：列出单位冲激函数、单位阶跃函数、复指数函数的MATLAB表达式。一、建模：1、单位冲激函数(t)无法直接用MATLAB描述，可以把它看作是宽度为(程序中用dt表示)，幅度为1/的矩形脉冲，即x1(t)二A(t-t1)=1/A11tt1+A0其余表示在t=t1处的冲激。2、单位阶跃函数：在t=t1处跃升的阶跃可写为u(t-t1)。定义为x2(t)=u(t-t1)=1t1lhnu=0;nh=lu-lh;elseiflulhnh=0;nu=lh-lu;elsenu=0;lh=0;enddt=input(输入时间间

9、隔dt=(例如0.5)lt=lmax;u=zeros(1,lt),uls,zeros(1,nu),zeros(1,lt);t1=(-lt+1:2*lt)*dt;h=zeros(1,2*lt),hls,zeros(1,nh);hf二fliplr(h);%反褶y=zeros(1,3*lt);fork=0:2*ltp=zeros(1,k),hf(1:end-k);%h(k-t)y1=u.*p*dt;yk=sum(y1);y(k+lt+1)=yk;subplot(4,1,1);stairs(t1,u)axis(-lt*dt,2*lt*dt,min(u),max(u),holdonylabel(u(t)

10、subplot(4,1,2);stairs(t1,p)axis(-lt*dt,2*lt*dt,min(p),max(p)ylabel(h(k-t)subplot(4,1,3);stairs(t1,y1)axis(-lt*dt,2*lt*dt,min(y1),max(y1)+eps)ylabel(s=u*h(k-t)subplot(4,1,4);stem(k*dt,yk)axis(-lt*dt,2*lt*dt,floor(min(y)+eps),ceil(max(y+eps)holdon,ylabel(y(k)=sum(s)*dt)ifk=round(0.8*lt)disp(暂停，按任意键继续)

11、，pauseelsepause(1),end运行结果如图3-15图3-1程序运行结果程序2：u=input(输入u数组u=);h=input(输入h数组h=);dt=input(输入时间间隔dt=);y=conv(u,h);%卷积和plot(dt*(l:length(y)-l),y),gridon;连续信号的卷积卷积积分运算实际上可用信号的分段求和来实现。用MATLAB实现连续信号f1(t)与f1(t)卷积的过程如下：1.将连续信号f1(t)与f1(t)以时间间隔A进行取样，得到离散序列f1(kA)与f1(kA)； 构造与f1(kA)和f1(kA)相对应的时间向量k1和k2(注意，此时时间序号

12、向量k1和k2的元素不再是整数，而是取样时间间隔A的整数倍的时间间隔点)；调用conv()函数计算卷积积分f(t)的近似向量f(nA)；构造f(n)对应的时间向量k。根据以上步骤我们自己可以比较容易写出求连续信号卷积的函数程序sconv()，该函数实现的功能为：在计算连续信号的卷积的同时绘出信号及卷积结果的波形，程序如下：functionf,k=sconv(f1,f2,k1,k2,p)%计算连续信号卷积积分f(t)=f1(t)*f2(t)%f:卷积积分f(t)对应的非零样值向量%k:f(t)的对应时间向量%f1:f1(t)的非零样值向量%f2:f2(t)的非零样值向量%k1:f1(t)的对应时

13、间向量%k2:f2(t)的对应时间向量%p:取样时间间隔f=conv(f1,f2);%计算需要序列f1与f2的卷积和ff=f*p;%将f进行采样k0=k1(1)+k2(1);%计算序列f非零样值的起点位置k3=length(f1)+length(f2)-2;%计算卷积和f的非零样值的宽度k=k0:p:k3*p;%确定卷积和f非零样值的时间向量subplot(2,2,1)plot(k1,f1)%在子图1绘f1(t)时域波形图title(f1(t)xlabel(t)ylabel(f1(t)subplot(2,2,2)plot(k2,f2)%在子图2绘f2(t)时域波形图title(f2(t)xla

14、bel(t)ylabel(f2(t)subplot(2,2,3)plot(k,f)%画卷积f(t)的时域波形title(f(t)=f1(t)*f2(t)xlabel(t)ylabel(f(t)给定下面两个连续信号，我们调用sconv()函数来求两连续信号的卷积。f1(t)=2U(t)-U(t-4)f2(t)=U(t)-U(t-2)MATLAB命令如下：p=0.01;k1=0:p:4;f1=2*ones(1,length(k1)k2=0:p:2;f2=ones(1,length(k2)%f1,f2两个信号的向量表示f,k=sconv(f1,f2,k1,k2,p)%子函数调用运行程序，绘制卷积结果

15、波形如图3.7所示。练习题：调用sconv()函数实现两连续信号的卷积。f1(t)=U(t)-U(t-2)f2(t)=U(t)-U(t-2)I*nt)f纱P123图3.7连续时间信号的卷积题目二：LTI系统的零状态响应设计要求：设二阶连续系统，其特性可用常微分方程表示：d2y/dt2+2（dy/dt）+8y=u求其冲激响应。若输入为u=3t+cos（0.11）,求其零状态响应。参考程序：clf,cleara=input（多项式分母系数向量a=）;b=input（多项式分子系数向量b=）;t=input（输入时间序列t=0:dt:tf）;u=input（输入序列u=）;tf=t(end);dt=

16、tf/(length(t)-1);r,p,k二residue(b,a);%完成有理多项式的部分分式展开h=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);subplot(2,1,1),plot(t,h);grid;y=conv(u,h)*dt;subplot(2,1,2);plot(t,y(1:length(t);grid运行该程序，取a=l,2,8,b=l,t=0:0.1:5及u=3*t+cos(0.1*t),所得结果如图3-2。 # # # #图3-2冲激响应和卷积法求输出 式4.3式4.1实验四典型周期信号的频谱表示实验目的：用MATLAB分析周期矩形脉冲、三角波脉冲

17、的频谱题目一：周期信号的频谱设计要求：周期电流、电压(统称其为信号)f(t)可展开为直流与各次谐波之和，即,AVf(t)，o+Acos(kt+Q)2kmkk,1式中=2k/T是基波角频率，T为周期。A,a2+b2kmkka=JT/2f(t)cos(kt)dt,k=0,1,2,.kT_t/2b=Jt/2f(t)sin(kt)dt,k=0,1,2,.kT_T/2周期信号的有效值定义为A，1JTf(t)2dt式4.2T0若用各谐波有效值1A:则表示为V2km丿A，A2+1A20V2km丿k,1全波整流电压Us(t)的波形如图4-1所示，用傅立叶级数可求得4U14Um,Am兀kmk21兀k2,4,6,

18、.k1,3,5,.A0kmU(t)二ms,1丄124n21n1cos2n0t式4.4可写出其展开式为（它只含直流和偶次谐波，令k=2n）若Um=100V,频率f=50Hz,（相应的T=0.02s,l=100,rad/s）,分别用式（4.2）和式（4.3）计算其有效值Usl和Us2（取至十次谐波），并求Us2的误差。参考程序：clear,formatcompact%数据显示的一种样式;空格比较少，数字排列比较紧凑Um=100;T=0.02;w=2*pi*5N=input(取的谐波次数N=);t=linspace(-T/2,T/2);dt=T/99;%linspace产生-T/2,T/2间100个

19、行矢量u=Um*abs(sin(w*t);fork=0:Na(k+1)=trapz(u.*cos(k*w*t)*dt/T*2;%trapz求数值积分b(k+1)=trapz(u.*sin(k*w*t)*dt/T*2;A(k+l)=sqrt(a(k+l)2+b(k+l)2);end0:N,A(1)/2,A(2:end)stem(0:N,a(1)/2,A(2:end)%绘制离散函数图形 #cleartf=10; %方法一：按傅立叶分析定义计算Usll=sqrt(trapz(u.2)*dt/T)%式4.2Usl2二sqrt(A(l)2/4+sum(A(2:end).2/2)%式4.3%方法二：按推导

20、出的全波傅立叶分量公式计算Us21=Um*sqrt(trapz(sin(w*t).2)*dt/T)%式4.2Us22=4*Um/pi*sqrt(0.52+0.5*sum(l./(4*l:3.2-l).2)%式4.3e=(Us21-Us22)/Us21运行程序，按提示输入。取得谐波次数N=10全波信号的波形图如图4-1所示，全波信号的各谐波分量如图4-2所示图4-1全波信号的波形图图4-2全波信号的各谐波分量题目二：非周期信号（方波）的频谱分析计要求：如图4-3a的矩形脉冲信号，求其在=-40rad/s40rad/s区间的频谱。参考程序：N=input(N=);%时间分隔的点数dt=10/N;t

21、=1:N*dt;f=ones(1,N/2),zeros(1,N/2);wf=input(wf=);%频谱宽度Nf=input(Nf=);%频谱点数w1=linspace(0,wf,Nf);dw=wf/(Nf-1);F1=f*exp(-j*t*w1)*dt;w=-fliplr(w1),w1(2:Nf);F=fliplr(F1),F1(2:Nf);subplot(1,2,1),plot(t,f,linewidth,1.5),gridsubplot(1,2,2),plot(w,abs(F),linewidth,1.5),grid程序运行结果：取时间分隔的点数N=256,需求的频谱宽度wf=40,需求

22、的频谱点数Nf=64,得出图4-3b10.90.80.70.60.50.40.30.20.10图4-3a矩形脉冲信号图4-3b程序运行结果 # #cleartf=10; # # #cleartf=10; #若取时间分隔的点数N=64，需求的频谱宽度wf=40,需求的频谱点数Nf=256,得出图4-4。 # cleartf=10; #图4-4程序运行结果 cleartf=10; #实验五傅立叶变换性质研究实验目的：学习用MATLAB的Fourier变换函数，验证Fourier变换的一些性质题目一：方波分解为多次正弦波之和设计要求：如图5-1所示的周期性方波，其傅立叶级数为：f(t)=4/冗sin

23、t+l/3(sin3t)+(l/(2k-l)sin(2k-l)t+k=1,2,用MATLAB演示谐波合成情况。参考程序：t=0:0.01:2*pi;y=sin(t);plot(t,y),figure(gcf),pausey=sin(t)+sin(3*t)/3;plot(t,y),pausey=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9;plot(t,y)y=zeros(10,max(size(t);x=zeros(size(t);fork=1:2:19x=x+sin(k*t)/k;y(k+1)/2,:)=x;endpause,figur

24、e(1),plot(t,y(1:9,:),gridline(0,pi+0.5,pi/4,pi/4)text(pi+0.5,pi/4,pi/4)halft二ceil(length(t)/2);pause,figure(2),mesh(t(l:halft),l:10,y(:,l:halft)程序运行结果如图5-2-pi2pi3pi-10.50图5-1周期性方波图5-2程序运行结果题目二：周期信号的滤波设计要求：如图5-3滤波电路，已知L=400mH,C=1OF,R1=2OO0.如激励电压us(t)为全波整流信号，Um=100V,1=100冗rad/s,求负载R两端的直流和各次谐波(它只含偶次谐波)

25、分量。Lr图5-3滤波电路参考程序：clear,formatcompactL=0.4;C=10e-6;R=200;Um=100;wl=100*pi;N=input(需分析的谐波次数2N=(键入偶数)；n=l:N/2;w=eps,2*n*wl;Us=4*Um/pi*0.5,T./(4*n.2-l);z1=j*w*L;z2=l./(j*w*C);z3=R;z23二z2.*z3./(z2+z3)UR=Us.*z23./(zl+z23)disp(谐波次数谐波幅度谐趣波相移(度)disp(2*0,n,abs(UR),angle(UR)*180/pi)程序运行结果：根据程序提示：需分析的谐波次数2N=（键

26、入偶数），如键入10后，得出结果。实验六离散信号分析实验目的：学习用MATLAB对离散信号分析的方法和编程方法。题目一：离散信号的MATLAB表述设计要求：编写MATLAB程序来产生下列基本脉冲序列：单位脉冲序列，起点nO，终点nf,在ns处有一单位脉冲(nOnsnf)。单位阶跃序列，起点n0,终点nf,在ns前为0,在ns后为1(n0nsnf)。复指数序列。参考程序：clear,n0=0;nf=10;ns=3;n1=n0:nf;x1=zeros(1,ns-n0),1,zeros(1,nf-ns);n2=n0:nf;x2=zeros(1,ns-n0),ones(1,nf-ns+1);n3=n0

27、:nf;x3=exp(-0.2+0.5j)*n3);subplot(2,2,1),stem(n1,x1);title(单位脉冲序列)subplot(2,2,3),stem(n2,x2);title(单位阶跃序列)subplot(2,2,2),stem(n3,real(x3);line(0,10,0,0)title(复指数序列)，ylabel(实部)subplot(2,2,4),stem(n3,imag(x3);line(0,10,0,0),ylabel(虚部)程序运行结果如图7-1 # #cleartf=10; # # #cleartf=10; #单悝胁沖序列复指敢序列单位阶氐庠列 # #cl

28、eartf=10; # #cleartf=10; #图7-1程序运行结果题目二：差分方程的通用递推程序设计要求：描述线性时不变离散系统的差分方程为:ay(n)+ay(nl)+ay(nn+l)=bu(n)+bu(nl)+bu(nn+l)12naa12nbb编写解上述方程的通用程序。参考程序：a=input(差分方程左端的系数向量a=a(1),.a(na)=);b=input(差分方程右端的系数向量b=b(1),.b(na)=);u=input(输入信号序列u=);nu=length(u);y=filter(b,a,u);stem(y),gridon执行此程序，输入a=1,0.1,0.15,-0.

29、225;b=3,7,1u=exp(0.1*1:20)程序运行结果如图7-2图7-2程序运行结果 # cleartf=10; #实验七离散系统的Z域分析 # #cleartf=10; # cleartf=10; #、目的(1)掌握利用MATLAB绘制系统零极点图的方法(2)掌握离散时间系统的零极点分析方法(3)掌握逆Z变换概念及MATLAB实现方法二、离散系统零极点线性时不变离散系统可用线性常系数差分方程描述，即,ay(n-i)=ii0,bx(n-j)jj07-1)其中y(k)为系统的输出序列，x(k)为输入序列。将式(7-1)两边进行Z变换的H(z)=Y(Z),bz-j=j=B(2X(z厂，a

30、z-丽ii07-2)将式(7-2)因式分解后有：叫z-q)H(z)=Cj1-叫z-p)i7-3)i1其中C为常数，q(j=1,2，,M)为H的M个零点，p(i1,2,N)为H的NJi个极点。系统函数H(z)的零极点分布完全诀定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性：系统单位样值响应h(n)的时域特性;离散系统的稳定性;1+ 1+ #cleartf=10; #离散系统的频率特性；三、离散系统零极点图及零极点分析1零极点图的绘制设离散系统的系统函数

31、为H(z)则系统的零极点可用MATLAB的多项式求根函数roots()来实现，调用格式为：p=roots(A)其中A为待求多项式的系数构成的行矩阵，返回向量P则是包含多项式所有根的列向量。31如多项式为B(Z)z2+Z+，则求该多项式根的MATLAB命令为：48A=13/41/8;P=roots(A)运行结果为：P=-0.5000-0.2500需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分母多项式均按z的降幕次序排列；另一种是分子、分母多项式均按Z，1的升幕次序排列。这两种方式在构造多项式系数向量时稍有不同。(1)H(z)按z的降幕次序排列：系数向量一定要由多项式最高

32、次幕开始，一直到常数项,缺项要用0补齐；如z4+3z3+2z2+2z+1其分子、分母多项式系数向量分别为A=1020、B=13221。(2)H(z)按z，i的升幕次序排列：分子和分母多项式系数向量的维数一定要相同，不足的要用0补齐，否则z0的零点或极点就可能被漏掉。如1+2z-1112z4z_21+ #1+ #cleartf=10; #1+ #1+ #cleartf=10; #其分子、分母多项式系数向量分别为A=120、B=11/21/4。 # cleartf=10; # cleartf=10; #用roots()求得H(z)的零极点后，就可以用plot()函数绘制出系统的零极点图。下面是求系

33、统零极点，并绘制其零极点图的MATLAB实用函数ljdt(),同时还绘制了单位圆。functionljdt(A,B)%Thefunctiontodrawthepole-zerodiagramfordiscretesystemp=roots(A);%求系统极点q=roots(B);%求系统零点p=p;%将极点列向量转置为行向量q=q;%将零点列向量转置为行向量x=max(abs(pq1);%确定纵坐标范围x=x+0.1;y=x;%确定横坐标范围clfholdonaxis(-xx-yy)w=0:pi/300:2*pi;t=exp(i*w);plot(t)axis(square)plot(-xx,0

34、0)plot(00,-yy)text(0.1,x,jImz)text(y,1/10,Rez)plot(real(p),imag(p),x)plot(real(q),imag(q),o)%用来清除图形的命令%确定坐标轴显示范围%复数定义exp(i*x)=cos(x)+i*sin(x)%画单位园%将当前坐标系图形设置为方形。横轴及纵轴比例是1：1%画横坐标轴%画纵坐标轴%画极点%画零点title(pole-zerodiagramfordiscretesystem)%标注标题holdoff例1：绘制如下系统函数的零极点(1)H(z)，3z35z2+10z10.5z1z3一3z2+7z一5H(z)，1

35、+3z一1+8z一2解：MATLAB命令如下A=1-37-5;B=3-5100;ljdt(A,B)绘制的零极点图如图7-1(a)所示。A=13/41/8;B=1-0.50;ljdt(A,B)绘制的零极点图如图7-1(b)所示。(a)(b)图7-1离散系统的零极点图2离散系统零极点分析离散系统零极点分布与系统稳定性信号与系统课程已讲到离散系统稳定的条件为：时域条件：离散系统稳定的充要条件为h(n)，即系统单位样值响应绝n=，对可和；Z域条件：离散系统稳定的充要条件为系统函数H(z)的所有极点均位于Z平面的单位圆内。对于三阶以下的低阶系统，可以利用求根公式求出系统函数的极点，从而判断系统的稳定性，

36、但对于高阶系统，手工求解则显得十分困难，这时可以利用MATLAB来实现。实现方法是调用前述的函数ljdt()绘出系统的零极点图，然后根据极点的位置判断系统的稳定性。例2：系统函数如例1所示，判断两个系统的稳定性。解：由例1绘出的零极点图可以看出两个系统的稳定性分别为：第(1)个系统不稳定；第(2)个系统稳定。零极点分布与系统单位样值时域特性的关系从信号与系统课程中已经得知，离散系统的系统函数H(z)与单位样值响应h(n)是一对Z变换对；因而，H(z)必然包含了h(n)的固有特性。离散系统的系统函数可以写成,(z-q)H(z)Cji(7-4),(z-p)ii1若系统的N个极点均为单极点，可将H(

37、z)进行部分分式展开为：H(z)=迓-(7-5)z-P由Z逆变换得：h(n)=迓k(p)nu(n)(7-6)iii1从式(7-5)和(7-6)可以看出离散系统单位样值响应h(n)的时域特性完全由系统函数H(z)的极点位置决定。从信号与系统的学习中已经得出如下规律:H(z)位于Z平面单位圆内的极点决定了h(n)随时间衰减的信号分量;H(z)位于Z平面单位圆上的一阶极点决定了h(n)的稳定信号分量;H(z)位于Z平面单位圆外的极点或单位圆上高于一阶的极点决定了h(n)的随时间增长的信号分量；下面以例子证明上述规律的正确性：例3：已知如下系统的系统函数H(z)，试用MATLAB分析系统单位样值响应h

38、(n)的时域特性。H(z)=丄，单位圆上的一阶实极点；z-1H(z)3-，单位圆上的一阶共轭极点；z2-2zcos()+1 #cleartf=10; #H(z)k，单位圆上的二阶实极点;4)单位圆内的一阶实极点；5)H(z)1(z-0.5)2单位圆内的二阶实极点；6)单位圆外的一阶实极点； # #cleartf=10; # # cleartf=10; #解：利用MATLAB提供的函数impz()绘制离散系统单位样值响应波形，impz()基本调用方式为(其他方式，请读者参看MATLAB帮助)：impz(b,a,N),其中，b为系统函数分子多项式的系数向量，a为系统函数分母多项式的系数向量，N为产

39、生序列的长度；需要注意的是，b和a的维数应相同，不足用0补齐，例如H(z)1(Z-1)2的b=001，a=1-21。下面是求解个系统单位样值响应的MATLAB命令:a=1-1;b=01;impz(b,a,10)运行结果如图7-2(a)所示。a=1-2*cos(pi/8)1;b=001;impz(b,a,50)运行结果如图7-2(b)所示。a=1-21;b=010;impz(b,a,10)运行结果如图7-2(c)所示。(4)a=1-0.8;b=01;impz(b,a,10)运行结果如图7-2(d)所示。(5)a=1-10.25;b=001;impz(b,a,10)运行结果如图7-2(e)所示。(6)a=1-1.2;b=01;impz(b,a,10)运行结果如图7-2(f)所示。_y6d._Jnn-堂naI音UL02d6Bn焙可niR忖功(b)(c)(d)