dir2fs.m

新版Matlab 7辅助信号处理技术与应用

function [C,B,A]=dir2fs(h);
%变h(n)值形式为频率取样形式
%[C,B,A]=dir2fs(h)
%C=并联部分的增益列向量
%B=按列排列的分子系数
%A=按列排列的分母系数
%h=FIR滤波器的冲激响应

%计算FIR滤波器冲激响应的频率采样
M=length(h);
H=fft(h,M);
magH=abs(H);
phaH=angle(H)';
%判断M的奇偶性
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)-(2*pi*k/M));
C=[2*magH(2:L+1),C1]';