dir2fs.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 23 行

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]'