批处理最小二乘Matlab代码

%批处理最小二乘参数估计（LS）
clear all;

a=[1 -1.5 0.7]‘; b=[1 0.5]‘; d=3; %对象参数
na=length（a）-1; nb=length（b）-1; %na、nb为A、B阶次

L=500; %数据长度
uk=zeros（d+nb1）; %输入初值：uk（i）表示u（k-i）
yk=zeros（na1）; %输出初值
x1=1; x2=1; x3=1; x4=0; S=1; %移位寄存器初值、方波初值
xi=randn（L1）; %白噪声序列

theta=[a（2:na+1）;b]; %对象参数真值
for k=1:L
    phi（k:）=[-yk;uk（d:d+nb）]‘; %此处phi（k:）为行向量，便于组成phi矩阵
    y（k）=phi（k:）*theta+xi（k）;