-300-16
-14-12-10-8Real Axis
-6-4-20
程序（3）画零极点图 num=[50];
den=conv([1,0],conv([1,10],[3,1])); sys=tf(num,den) sys1=feedback(sys,1) pzmap(sys1) gird 图（3）
Pole-Zero Map
1.5
1
0.5
Imaginary Axis
-0.5
-1
-1.5-12
-10-8-6Real Axis
-4-20
程序（4）
计算相角裕量和幅值裕量 num=[50];
den=conv([1,0],conv([1,10],[3,1])); sys=tf(num,den) sys1=feedback(sys,1)
[gm,pm,wcg,wcp]=margin(sys) 运行结果 Transfer function: 50 --------------------- 3 s^3 + 31 s^2 + 10 s
Transfer function: 50
-------------------------- 3 s^3 + 31 s^2 + 10 s + 50 gm =
2.0667 pm =
7.5615 wcg =
1.8257 wcp =
1.2645 程序（5） 绘制阶跃响应曲线 num=[50];
den=conv([1,0],conv([1,10],[3,1])); sys=tf(num,den) sys1=feedback(sys,1) step(sys1)
图（5）
Step Response
21.8
1.61.41.2
Amplitude
10.80.60.40.20
0102030
Time (sec)
40506070
【7.4】 程序如下： num=[300];
den=conv([1,0,0],conv([0.2,1],[0.02,1])); sys=tf(num,den) margin(sys) grid
波特图如下：
Bode Diagram
Gm = Inf , Pm = -78 deg (at 11 rad/sec)
10050
Magnitude (dB)Phase (deg)
0-50-100-150-180-225
-270-315-360
10
-1
10
10
1
10
2
10
3
Frequency (rad/sec)
【9.3】程序： A=[-2 2 -1;0 -2 0;1 -4 0]; B=[0;0;1]; C=[1,0,0]; D=0; M=ctrb(A,B) m=rank(M) if m==3;
disp('系统可控') else
disp('系统不可控') end
N=obsv(A,C) n=rank(N) if n==3;
disp('系统可观')
else
disp('系统不可观') end
sys=ss(A,B,C,D)
[num,den]=ss2tf(A,B,C,D) sys1=tf(num,den) [z,p,k]=ss2zp(A,B,C,D) 运行结果： M =
0 -1 2 0 0 0 1 0 -1 m = 2
系统不可控 N =
1 0 0 -2 2 -1 3 -4 2 n =