getval.m

来自「关于系统辨识的一些仿真程序」· M 代码 · 共 19 行

function [z,A,B]=getval(num,den)
  z=zeros(1,16);
  C=tf(num,den);
  D=c2d(C,0.01);
  [B,A]=tfdata(D,'v');
  u=mseries(16);
  for k=length(A):16
    for i=2:length(A)
      z(k)=z(k)-A(i)*z(k-i+1);
    end
    for j=1:length(B)
      z(k)=z(k)+B(j)*u(k-j+1);
    end
  end