⑥p1=1+4j,p2=1-4j
z = [] p = [1+4j,1-4j] k = [1]
[b,a] = zp2tf(z,p,k) sys = tf(b,a) pzmap(sys) impulse(sys)
答:由程序执行结果可以看出,在无零点的情况下: 当极点唯一且在原点时,h(t)为常数;
当极点唯一且是负实数时,h(t)为递减的指数函数; 当极点唯一且是正实数时,h(t)为递增的指数函数; 当H(s)有两个互为共轭的极点时,h(t)有sint因子;
当H(s)有两个互为共轭的极点且他们位于右半平面时,h(t)还有?????因子; 当H(s)有两个互为共轭的极点且他们位于左半平面时,h(t)还有?????因子。
(4)已知连续时间系统的系统函数分别如下: ①?? ?? =②?? ?? =③?? ?? =
????+???????????+????+????????+????+????????+????+????
上述三个系统具有相同的极点,只是零点不同,试用MATLAB分别绘制系统的零极点分布图及相应冲激响应的时域波形,观察并分析系统函数零点位置对冲激响应时域特性的影响。 ①?? ?? =
??????+????+????
MATLAB程序如下: a = [1 2 17] b = [1] sys = tf(b,a) subplot(211) pzmap(sys) subplot(212) impulse(b,a) 程序执行结果如下:
②?? ?? =
??+??
????+????+????
MATLAB程序如下: a = [1 2 17] b = [1 8] sys = tf(b,a) subplot(211) pzmap(sys) subplot(212) impulse(b,a) 程序执行结果如下:
③?? ?? =
?????????+????+????
MATLAB程序如下: a = [1 2 17] b = [1 -8] sys = tf(b,a) subplot(211) pzmap(sys) subplot(212) impulse(b,a) 程序执行结果如下:
由程序执行结果看出,当极点不变时,零点分布只影响系统时域响应的幅度和相位,对时域响应模式没有影响。
不会改变是衰减振荡还是增长振荡。
四,心得体会
MATLAB在拉普拉斯变换处又一次化繁为简,简化了繁杂的计算,奖结果直观的呈现在了我的眼前。