ch4_1.m

这是书籍:Matlab控制系统与应用的源码,

A=[1 -0.5 0.7];B=[0 1 0.5];
th0=poly2th(A,B);
u=idinput(300, 'rbs');
y=idsim([u,randn(300,1)],th0);
z=[y u];
ir = cra(z);
m = arx(z,[2 2 1]);
imp = [1;zeros(20,1)];
irth = sim(m,imp);
subplot(211)
plot([ir irth])
title('impulse responses')
subplot(212)
plot([cumsum(ir),cumsum(irth)])
title('step responses')