getrz.m

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

function rz=getrz(c,x1,x2,u)
  rz=zeros(1,16);
  for k=x1:16
    for i=2:x1
      rz(k)=rz(k)-c(i-1)*rz(k-i+1);
    end
    for j=1:x2
      rz(k)=rz(k)+c(x1+j-1)*u(k-j+1);
    end
  end