dir2fs.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 20 行

function [C,B,A]=dir2fs(h)
M=length(h);
H=fft(h,M);
magH=abs(H);phaH=angle(H)';
if(M==2*floor(M/2))
  L=M/2-1;
  A1=[1,-1,0;1,1,0];
  C1=[real(H(1)),real(H(L+2))];
else
  L=(M-1)/2;
  A1=[1,-1,0];
  C1=[real(H(1))];
end
K=[1:L]';
B=zeros(L,2);A=ones(L,3);
A(1:L,2)=-2*cos(2*pi*K/M);A=[A;A1];
B(1:L,1)=cos(phaH(2:L+1));
B(1:L,2)=-cos(phaH(2:L+1));
C=[2*magH(2:L+1),C1]';