自动控制原理 matlab仿真

1. 闭环反馈

clear;clc;
num1=6;den1=[1 2 3 1];
num2=1;den2=[1 2];
gs=tf(num1,den1);hs=tf(num2,den2);
gshs=series(gs,hs);
figure(1)
nyquist(gshs);
figure(2)
g=feedback(gs,hs);
step(g)

2 阶阻尼系统

clear;clc;
wn=10;
for zeta=[0.1 0.5 1 2]
    num=wn^2;
    den=[1 2*zeta*wn wn^2];
    figure(1)
    bode(num,den)
    hold on
end
zeta=0.707;

for wn=[1 5 10 20]
    num=wn^2;
    den=[1 2*zeta*wn wn^2];
    figure(2)
    bode(num,den)
    hold on
end

3. 典型2阶系统

g(s)=(wn^2)/(s^2+2*wn*s+wn^2)

clear;clc;
wn=10;
for zeta=[0.1 0.5 1 2]
    num=wn^2;
    den=[1 2*zeta*wn wn^2];