ex2_5.m

华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 1. 安装MATLAB软件。安装时

% ex2_5
% 过程输入是伪随机双位信号，建立过程的ARX模型
figure('pos',[10,10,300,500],'color','w');
load arxdata; % 输入输出数据文件名为arxdata
axes('pos',[.1,.7,.78,.2]);
t=length(y);
[x,xx]=stairs(u);
plot(x,2*xx,1:t,y,'k');
grid;
title('原始数据');
legend('输入','输出');
m=arx([y u],[3,3,1]);
[A,B]=polydata(m)
m1=idpoly(A(1:3),B(1:3));
e=0.02*randn(t,1);
y1=idsim([u e],m1);
axes('pos',[.1,.4,.78,.2]);
plot(1:t,y,'b--',1:t,y1,'k');
grid;
title('实际和估计输出比较');
legend('实际','估计');
axes('pos',[.1,.1,.78,.2]);
plot(1:t,y1-y);
grid;
title('误差');
;