matlab作业题

num=[-1 -2];
den=conv([1 3],[1 2 2]);
r=rlocus(num,den);
plot(r)
v=[-6 6 -2 2];axis(v)
grid
title('Root-Locus Plot')
xlabel('Real Axis')
ylabel('Imag Axis')

num=[5];
den=[1 16 56 36];
a=2:10;
r=rlocus(num,den,a);
plot(r,'o')
%v=[-12 12 -12 12];axis(v)
grid
title('Root-Locus Plot')
xlabel('Real Axis')
ylabel('Imag Axis')

num=[1 1];
den=conv([1 -1 0],[1 4]);
rlocus(num,den);
r=rlocus(num,den);
[K,r]=rlocfind(num,den)
sysg=tf(K*num,den);
sys=feedback(sysg,[1]);
step(sys)

w=logspace(0,1);wn=0.7;zeta=[0.1 1.0 2.0];
num=[wn.^2];
for k=1:length(zeta)
    den=[1 2*zeta(k)*wn wn.^2];
    bode(num,den,w);
    hold on
end
grid;
title('Bode Plot')
hold off
    

num=[3.6];
den=[1 3 5];
sys=tf(num,den);
w=logspace(-1,2,100);
bode(sys,w);
[Gm,pm,wcp,wcg]=margin(sys);
GmdB=20*log10(Gm);
[Gm,pm,wcp,wcg]
[mag,phase,w]=bode(sys,w);
[Mp,k]=max(mag);
resonant_peak=20*log10(Mp);
resonant_frequency=w(k);

close all;clc;clear all;
num=[0.0167 1];
den=conv([0.03 1 0],conv([0.0025 1],[0.001 1]));
sys1=tf(num,den);
K=[5 500 800 3000];
for n=1:4
    sys=sys1*K(n);
    w=logspace(-1,2,100);
    figure(1);hold on;
    bode(sys,w)
    [Gm,pm,wcp,wcg]=margin(sys)
    GmdB=20*log10(Gm);
    [Gm,pm,wcp,wcg]
end