阶跃响应面积法辨识传递函数&matlab编程

 

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

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

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

大多数自平衡对象的y(t)可近似为  a_n*dⁿy(t)/dt+a_n-1 *d^n-1y(t)/dt+······+a₁*dy(t)/dt+y(t)=u(t)

现以三阶系统为例

将y(t)移至右边，在[0,t]上积分    a₃*d²y(t)/dt+a₂*dy(t)/dt+a₁y(t)=∫^t₀[1-y(t)]dt

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

所以当t趋于无穷时,a₁=∫^t₀[1-y(t)]dt;

同理可求a_n

matlab实现：

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