matlab程序 梯度校正参数辨识程序

%FLch4GAeg1
clear
u=[-1-1-1-1111-111-1-11-11-1-1-1-11];
y=[0-2-6-7-7-3573-153-5-31-11-5-7-7];
%画出u和y图形
figure（1） subplot（211） stem（u） subplot（212） stem（y） hold on
k=1:20
plot（ky）
%给出初始值
h1=[-100]‘;h2=[-1-10]‘;g=[000]‘;I=[100;01/20;001/4];
h=[h1h2zeros（316）];
%计算输入样本数据h（k）
for k=3:18
    h（:k）=[u（k）u（k-1）u（k-2）]‘;
end
%计算出权矩阵R（k）和估计值g
for k=1:18
a=h（1k）^2+（h（2k）^2）/2+（h（3k）^2）/4;%按照式（4.45）开始计算权矩阵
a1=1/a;
R=a1*I;%按照式（4.45）计算出权矩阵
g（:k