利用阶跃响应曲线对系统进行辨识，有切线法、两点法和面积法。

面积法的原理是是通过求取微分方程的系数来辨识传递函数：

首先求得放大系数，得到无因次阶跃响应y(t)

大多数自平衡对象的y(t)可近似为  an*dny(t)/dt+an-1 *dn-1y(t)/dt+······+a1*dy(t)/dt+y(t)=u(t)

现以三阶系统为例

将y(t)移至右边，在[0,t]上积分    a3*d2y(t)/dt+a2*dy(t)/dt+a1y(t)=∫t0[1-y(t)]dt

由初始条件可知，当t趋于无穷时，d2y(t)/dt=dy(t)/dt=0，y(t)=1

所以当t趋于无穷时,a1=∫t0[1-y(t)]dt;

同理可求an

matlab实现：

num=[1];den=[960 368 36 1];disp('原系统：')Gs=tf(num,den)     %原系统的传递函数；h=yout(:,1);t0=0.05;t=0:t0:300;