clc;
clear;
load?data.txt;
u=data（:1）;%input
y=data（:2）;%output
n=2;
N=..;??%N為數據長度－n
r=2;

for?i=1:N????????%用最小二乘法初步估計a1a2b1b2
????Y（i1）=y（i+n）;
????Si（i1）=y（i+1）;
????Si（i2）=y（i）;
????Si（i3）=u（i+1）;
????Si（i4）=u（i）;
end
sita?=inv（Si‘*Si）*Si‘*Y;
sitaN=sita;

%用廣義最小二乘法求sita
deta?=?1.0;
num=0;
while（deta>0.001）
????num=num+1;