(完整word版)信号与系统matlab实验习题4连续系统的零极点分布与频响特性的关系.doc

连续系统零极点分布与频响特性的关系 （答案在下方） 连续系统零极点图的画法 pzplot()函数可用来绘制连续系统的零极点图，具体用法如下： pzplot(SYS) 计算 LTI 模型 SYS 的零点和极点并绘制零、极点图。其中 SYS 的 产生可以采用传递函数法： SYS=tf(b,a)，b 和 a 分别为系统函数的分子多项式和分母多项式系数矩阵。 例 1：系统函数为 H s 2s2 5s 12 ，完成一幅适当标注的零极点图。 4s2 14s 20 s3 解： MATLAB 代码如下： b=[2 5 12]; a=([1 4 14 20]); SYS=tf(b,a); pzplot(SYS); 系统的零极点图如图 1 所示。 Pole-Zero Map 3 2 s 1 x A y r 0 a n g a m I -1 -2 -3 -2 -1.5 -1 -0.5 0 Real Axis 图 1 例 1 系统的零极点图 由系统的零极点分布决定系统的频率响应特性 稳定的连续时间系统的频率响应特性可以由系统函数得到 m m s z j j ω zj H j ω H s s jω K jn 1 s j ω K jn 1 (1) s pi j ω pi i 1 i 1 令分子中每一项 j ω zj N j ejψj , 分母中每一项 j ω pi M i ej θi ，则 jψ jψ j ψ H j ω K N1 e 1 N2 e 2 L Nm e m (2) jθ jθ j θ M 1 e 1 M 2 e 2 L M n e n H j ω K N1 N 2 L N m (3) M 1M 2 L M n ω ψ1 ψ2 L ψm θ1 θ2 L θn (4) 分析频率响应特性的方法： 1. H j ω H s s jω ，带入数值，得到 H j ω ～ω的分布； 根据零极点图中零极点的分布，用几何的方法定性判断系统的频率响应特性； 对模拟系统， MATLAB 信号处理工具箱提供了 freqs()函数是用来求取模拟滤波器的频率响应。具体函数的用法是： H=freqs(B, A, W) 返回模拟滤波器的频率响应 H j ω ，即复数频率响应矩阵 H，其中的 B、A 是已知系统的传递函数模型， W 为已选定的频率点。如果直 接使用 freqs(B, A, W) ，则直接绘制幅度响应和相位响应曲线，不返回任何值。 如果缺省频率范围 W，函数将自动选用一组 200 个频率点范围。 freqs(B,A,N) 选择 N 个频率点。 B,A 矩阵的写法如下 b 1 snb 1 b 2 snb 2 H s a 2 sna a 1 sna 1 2  b nb L a na 例如： (a) s 0.5 B [1,0.5]; A [1,3,2]; H s 3s 2 s2 (b) H s 1 B [1]; A [1,3,2]; 3s s2 2 (c) H s s B [1 0]; A [1,3,2]; 3s s2 2 注意： (c)的 B 矩阵由于缺少常数项，故补了一个 0。 例 2：已知 C 0，利用 MATLAB 的 freqs 函数分析如下系统的频率响应特性。 (1) H 1 s C ; (2) H 2 s s ; s C s C 解： MATLAB 代码如下： close all w_C=100; B1=[w_C];A1=[1 w_C];SYS1=tf(B1,A1); B2=[1 0];A2=[1 w_C];SYS2=tf(B2,A2); pzplot(SYS1);axis([-2.1*w_C 2.1*w_C -w_C w_C]); figure;pzplot(SYS2);axis([-2.1*w_C 2.1*w_C -w_C w_C]); figure;freqs(B1,A1); figure;freqs(B2,A2); MATLAB 绘制的零、极点图和频率响应特性曲线如图 2 所示。 e 100 d u n g Pole-Zero Map a M 10-1 100 1) 1 2 3 s 10 10 10 d 50 Frequency (rad/s) n ) o c s e 0 s e s 0 r i g x e A d -50 y ( r e s i -50 g a a -100 m h I P 101 102 103 -100 -100 0 100 200 -200 Real Axis (seconds -1 )  Frequency (rad/s) (a) (b) e 100 d u -2 n 10 g Pole-Zero Map a 100 M 10-4 0 2 -