ch4_1.m

来自「清华大学《matlab 控制系统应用与实例》第一部分M文件的源码」· M 代码 · 共 16 行

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')