MATLAB课例子程序.docVIP

%% MATLAB Program 4-1 num=[1 3]; den=[1 5 20 16 0]; r=rlocus(num,den); plot(r,-); v=[-6 6 -6 6]; axis(v); axis(square); gtext(o) gtext(x) gtext(x) gtext(x) gtext(x) %% MATLAB Program 4-2 num=[1]; den=[1 1.1 10.3 5 0]; r=rlocus(num,den); plot(r,o) v=[-6 6 -6 6];axis(v) grid title(Root-Locus Plot of G(s)=K/[s(s+0.5)(s^2+0.6s+10)]) xlabel(Real Axis) ylabel(Imag Axis) %% MATLAB Program 4-3 num=[1]; den=[1 1.1 10.3 5 0]; K1=0:0.2:20; K2=20:0.1:30; K3=30:5:1000; K=[K1 K2 K3]; r=rlocus(num,den,K); plot(r,o) v=[-4 4 -4 4];axis(v) grid title(Root-Locus Plot of G(s)=K/[s(s+0.5)s^2+0.6s+10)]) xlabel(Real Axis) ylabel(Imag Axis) %% MATLAB Program 4-4 A=[0 1 0;0 0 1;-160 -56 -14]; B=[0;1;-14]; C=[1 0 0]; D=[0]; K=0:0.1:400; r=rlocus(A,B,C,D,K); plot(r,-); v=[-20 20 -20 20]; axis(v) grid gtext(o) gtext(x) gtext(x) gtext(x) %% MATLAB Program 4-5 num=[1]; den=[1 3 2 0]; numa=[1]; dena=[1 3 3 1]; K1=0:0.1:0.3; K2=0.3:0.005:0.5; K3=0.5:0.5:10; K4=10:5:100; K=[K1 K2 K3 K4]; r=rlocus(num,den,K); a=rlocus(numa,dena,K); y=[r a]; plot(y,-) v=[-4 4 -4 4];axis(v) grid title(Root-locus Plot of G(s)=K/[s(s+1)(s+2)]and Asymptotes) %% MATLAB Program 4-6 num=[1]; den=[1 3 2 0]; numa=[1]; dena=[1 3 3 1]; K1=0:0.1:0.3; K2=0.3:0.005:0.5; K3=0.5:0.5:10; K4=10:5:100; K=[K1 K2 K3 K4]; r=rlocus(num,den,K); a=rlocus(numa,dena,K); plot(r,o) a=rlocus(numa,dena,K); v=[-4 4 -4 4];axis(v) hold %current plot held plot(a,-) grid title(Root-Locus Plot of G(s)=K/[s(s+1)(s+2)] and Asympototes) xlabel(Real Axis) ylabel(Real Axis) %% MATLAB Program 4-7 num=[1 4 4]; den=[1 10 29 40 100]; r=rlocus(num,den); plot(r,o) v=[-8 4 -6 6];axis(v);axis(square) hold current plot held plot(r,-) grid title(Root-Locus Plot of G(s)=(s+2)^2/[(s^2+4)(s+5)^2]) xlabel(Real Axis) ylabel(Imag Axis) %% MATLAB Program 4-8 num=[1 0 25 0]; den=[1 0 404 0 1600]; K=0:0.4:1000; r=rlocus(num,den,K); plot(r,o); v=[-30 20 -25 25];axis(v);axis(square) grid title(Root-locus Plot of G(s)=K(s^2+25)s/(s^4+404s^2+1600)) xlabel(Real Axis