最小二乘法辨识一阶系统-lsq.m

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;