clear
clf??
%定義預測長度?和控制長度
N2=13;????????????????????%?預測長度??（分別取5，6，7）
Nu=5;????????????????????%控制長度??（分別取1，23）
gama=0.7;??????????????????%制加權系數
lam=0.998;???????????????%遺忘因子
P=（1e+6）*eye（N2+1）;??????%P的初始值
%參數初始值
uuu=0;yyy=0;
uu=zeros（N21）;u=zeros（Nu1）;
yy=zeros（N21）;y1=zeros（N21）;
theta=zeros（N2+11）;theta（11）=1;theta（N2+11）=1;

%產生周期為100，時間為T，幅值為1的方波信號的給定值
?T=600;[yr0t]=gensig（‘square‘100T1）;
?nm=length（t）;??????????%確定循環次數
?for?ij=2:nm
?????yr=yr0（ij）+1;??????%產生周期為100，時間為T，幅值為在1和2之間的方波信號
%根據系統模型，計算k時刻的輸出值y（k）
?????y=1.5*yy（N21）-0.5*yy（N2-11）+0.5*uu（N2-41）;
%保存k時刻以及以前的N1個輸出值y（k）y（k-1）...y（k-N2）以供模型運算
?????for?i=1:N2-1
?????????yy（i1）=yy（i+11）;
?????end
?????yy（N21）=y;
?????yyy=[yyy;y];??????%保存各k時刻的nm個輸出量以便繪圖
%根據最小二乘公式，由y（k）計算F的各元素的值f0