%批處理最小二乘參數估計（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）;