实验八由系统函数零极点分布决定时域频域特性

1、实验八由系统函数零、极点分布决定 时域、频域特性实验目的用MATLAB模拟系统函数零、极点分布决定时域、频域特性的过程实验原理计算零、极点可以用 roots函数，若参数为传统函数 H(s)的分子多项式系数 b,则得到 零点；若为分母多项式系数 a,则得到极点。 MATLAB还提供了 zero(sys)和pole(sys)函 数直接计算零、极点，其中 sys表示系统传递函数。另外，【p,z】 =pzmap(sys)函数也具有 计算极点p和零点z的功能。不带返回值的 pzmap(sys)则绘制出系统的零、极点分布图。零极点和传递函数的多项式系数一样，可以作为LTI系统的描述方法。MATLAB提供了

2、( b,a) =zp2tf(z,p,k)和【z,p,k =tf2zp(b,a)两个函数用于在上述两种描述方法之间进 行转换，其中k为用零、极点表示传函数时的系统增益。实验过程1.当F(s)极点(一阶)落于图所示s屏幕图中的各方框所处位置时，画出对应的f 波形填入方框中。解：clear all , close all , clc;t = 0:.1:40'figure, id = 1;for omega = .5:-.25:0for sigma = -.06:.03:.06p = sigma + j*omega;if omega = 0p = p;p' endb a = zp2tf

3、(,p,1);subplot(3,5,id);impulse(b,a,t);set(gca, 'YLim' ,-20,20);endend输出结果如下图所示，可见随着极点从虚轴左侧移动到右侧，其冲激响应有衰减变为 发散；随着极点由实轴向上、下两侧移动，冲激响应由单调变化转为振荡变化，且振荡 周期逐渐减小。Figure 1File Edit View Insert Tools Desktop .Vmdow Help。仃“3氐叫口国国用 Fnzld 与 五专£lpln*d 与impulEe Response ImpuLse Respofise Impulse Respon

4、se Impukse Response Impulse ResponseO 2040 O 2040 D 20 4C 020400Tine (sec) Time (&ec) Time sec) Time (sec) Time (sec)Impjlse Response Impukse Response Impulse Response Impulse Response Impulse ResponseTime (sec) Time (s&c) Tkne sec)Time eecjTime (sec)impulse Response impulse Respcnae impulse

5、 Response Response Impulse Response02040020 4Q 0204002040020 40Time (sec) Time (sec) Time tsec) Time Uec) line (sec)MATLA4口识点（9）-虚数单位i,j以及常量恢复和数学公式中经常以i、j作为下标一样，包才Fortran和C语言在内的打部分机器语言教科书都习惯的把这两个字符作为循环变量使用，但在MATLAB i和j被默认为定义为表示虚数单位的常量。虽然用户还可以将其自定义成任意类型变量或者数据结构，但这样的很容易带来混乱。请看下面的命令和输出：Clear alla=1+2*i

6、a=1.0000+2.0000ii=1;b=1+2*ib=3Sqrt(-1)ans=0+1.0000iC=1+2iC=1.0000+2.0000i可见MATLA除了把i定义为表示虚数单位的数值，在解释表达式时还将其理解为表示虚数单位的符号，所以把i和j被用户定义使用并不影响复数的表达。但为了避免混淆，还是建议调整自己的编程习惯，不再使用i和j作为循环变量。MATLAB提供了 freqs(b,a)函数用来绘制系统的频率响应，包括幅度响应和相位响应， 其中b和a分别对应传递函数的分子和分母多项式系数。如果将调用方式改为 H=freqs(b,a ,omega),则不绘图输出，而是计算抽样点omega

7、处得频响并传递到 H中。2.若H(s)的零极点分布如图，试讨论它们分别是哪种滤波网络(低通、高通、带通、 带阻)。解：clear all, close all, clc;data = struct('title','(a)','(b)','(c)','(d)','(e)','(f)','(g)','(h)','zeros',口,0,0;0,卜0.5,0,1.2j;-1.2j,0;0,1.2j;-1.2j,'poles',

8、-2;-1,-2;-1,-2;-1,-2;-1,卜1+j;-1-j,-1+j;-1-j,-1+j;-1-j,j;-j);omega = 0.01:0.01:6'figure;for id = 1:8b,a = zp2tf(data(id).zeros,data(id).poles,1);H = freqs(b,a,omega);subplot(4,2,id);plot(omega,abs(H);set(gca,'YScale','log','FontSize',16);title(data(id).title);xlabel('omega');ylabel('H(omega');end输出结果如下所示，可见八个系统分别是：低通、带通、高通、带通、