% 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('誤差');