ex2_7a.m

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

% ex2_7
% 用P.R.B.S.伪随机双位信号作为控制输入信号，
% 对上述过程进行系统辨识，获得了输出y，并存放在bjdata文件中。
% 根据输入输出数据建立BJ模型
figure('pos',[10,10,300,500],'color','w');
load bjdata; % 输入输出数据
axes('pos',[.1,.7,.8,.2]);
t=length(y);
[x,xx]=stairs(u);
plot(x,2*xx,1:t,y,'k');
grid;
title('原始数据');
legend('输入','输出');
m=pem([y u],[2,2,2,2,2,1])
e=0.02*randn(t,1);
y1=idsim([u e],m);
axes('pos',[.1,.4,.8,.2]);
plot(1:t,y,'b--',1:100,y1,'k');
grid;
title('实际和估计输出比较');
legend('实际','估计');
axes('pos',[.1,.1,.8,.2]);
err=max(abs(y1-y))
plot(1:t,y1-y,'k');
grid;
title('误差');