invztrans.m

来自「wav声音文件的通过多脉冲激励编码器」· M 代码 · 共 15 行

function hn=invztrans(lina,p)
syms z;
syms s;
syms n;
for i=1:p
    if i==1
        s=lina(i)*z^(-i);
    else 
        s=s+lina(i)*z^(-i);
    end
end
HZ=1/(1-s);

hn=iztrans(HZ,z,n);