gmatrix.m

关于系统辨识的一些仿真程序

function HL=gmatrix(z,x1,x2)
b=length(z);
x3=x1+x2-1;
a=b-x1+1;
u=mseries(16);
HL=zeros(a,x3);
for j=1:(x1-1)
    for i=1:a
        HL(i,j)=-z(x1-j+i-1);
    end
end
for j=x1:x3
    for i=1:a
        HL(i,j)=u(x2+x1-j+i-1);
    end
end