matlab 二阶系统 伯德图 阶跃响应

    目的是为了分析不同阻尼比对阶跃响应的影响，以及观察二阶系统的频率响应。

   二阶系统传递函数如下：

                              

    系统框图：

                         

  （1）  系统阶跃响应，不同阻尼比对应不同的输出。具体参考自动控制原理。当a=1，无超调的到达，可以作为滤波器，也是自抗扰中的跟踪微分器原理。

             

    （2）伯德图，是一个低通滤波器。

             

附实现代码：

clc;
clear;

%% G(s) = wn^2 / (s^2 + 2*wn*a + wn^2)

wn = 100; %fre
a = 1;  %damping ratio
num = [0 0 wn^2];
den = [1 2*wn*a, wn^2];
sys1 = tf(num,den);

a = 0.5;  %damping ratio
num = [0 0 wn^2];
den = [1 2*wn*a, wn^2];
sys2 = tf(num,den);

%% plot
bode(sys1)
grid 
xlabel('频率f');title('   ');

figure
step(sys1,sys2,1)
grid on;
xlabel('时间t');title('   ');