example1.m

是关于系统辨识的一篇论文和本人对文中所提论文方法的Matlab程序验证

% 算例1
clf;
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);
th=arx(z,[2 2 1]);
%
imp=[1;zeros(19,1)];
irth=idsim(imp,th);
%
subplot(211)
plot(ir,'--');
hold on
plot(irth,'r')
legend('相关分析','ARX 建模')
title('impulse responses')
%
subplot(212)
plot(cumsum(ir),'--');
hold on
plot(cumsum(irth),'r')
legend('相关分析','ARX 建模')
title('step responses')