ex2_6.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 24 行

% ex2_6
% 过程输入是随机双位信号，建立过程的ARMAX模型
figure('pos',[10,10,300,500],'color','w');
load armaxdata; % 输入输出数据
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=armax([y u],[2,2,2,1])
e=0.02*randn(t,1);
y1=idsim([u e],m);
axes('pos',[.1,.4,.78,.2]);
plot(1:t,y,'b--',1:100,y1,'k');
grid;
title('实际和估计输出比较');
legend('实际','估计');
axes('pos',[.1,.1,.78,.2]);
plot(1:t,y1-y);
grid;
title('误差');