iir_dirtime.m

来自「the file contain many matlab signal sour」· M 代码 · 共 18 行

function [a,b]=IIR_dirTime(hd,M,N)
hd=hd(1:M+N+1));
h=zeros(M+N+1,N+1);
for n1=1:M+N+1
    for n2=1:N+1
        if n1>n2
         h(n1,n2)=hd(n1-n2+1);
     end
 end
end
h1=h(1:M+1,:);
h2=h((M+2,M+N+1,:);
h11=h2(:,1);
h3=h2(:,2:N+1);
a=-h3\h11;
a=[1,a'];
b=(h1*a')';