P=
-1.5000+1.3229i -1.5000-1.3229i
需要注意的是，系数向量A的元素一定要由多项式的最高幂次开始直到常数项，缺项要补零。例如：
A(S)?S6?3S4?2S2?S?4
应表示为： A=[1 0 3 0 2 1 -4]
用roots（）函数求得系统函数H（S）的零极点后，就可以用plot命令在复平面上绘制系统函数的零极点图，方法是在零点位置标以符号“O”，而在极点位置标以“X”。
三．实验内容
1、通过u（t）-u（t-2）的拉氏变换曲面图，观察频域和复频域的关系。
S2?42、通过F(S)?3,画出零极点图。
S?5S2?10S?50
四．实验设备
PC机一台能全速运行MATLAB
五．实验程序
1、u（t）-u（t-2）的拉氏变换曲面图 a=0:0.1:5; b=-20:0.13:20; [a b]=meshgrid(a,b); c=a+i*b;
c=(1-exp(-2*c))./c; c=abs(c);
31
mesh(a,b,c);
F(S)?S22、?4S3?5S2?10S?50曲面图 a=-6:0.1:6; b=-5:0.1:5; [a b]=meshgrid(a,b); s=a+i*b; d=s.^2-4;
e=s.^3-5*s.^2+10*s-50; c=d./e; c=abs(c); mesh(a,b,c);
S23、F(S)??4S3?5S2?10S?50零极点图
a=[1 -5 10 -50]; b=[1 0 -4]; p=roots(a); z=roots(b); p=p'; z=z';
x=max(abs([p,z])); x=x+0.1; y=x; hold on axis([-x,x,-y,y]); plot([-x,x],[0,0]); plot([0,0],[-y,y]); plot(real(p),imag(p),'x');
32
plot(real(z),imag(z),'o'); grid off
33
34