fb.m

设计一个原型滤波器

%-----Filter Bank Design for Subband Adaptive Microphine Arrays----%
%      h=fb(M1,D1,L,wp)                               %
%      h为原型分析滤波器,                 %
%      M1为子带分解个数,D1为抽取个数,L为分析综合滤波器的长度              %
%      wp带通滤波器的截止频率
function h=fb(M1,D1,L,wp)
M=M1;%子带个数
D=D1;%过采样率
Lh=L;%分析滤波器长度(阶数)
Th=(Lh)/2;%分析滤波器组的群延时
wp=wp;%带通滤波器的分离频率点,这里我们可以选多个进行试验
%%%%%%%%%%%%%%%%%%以下是分析理想原型滤波器的设计%%%%%%%%%%%%%%%%%%%%%%
for i=1:Lh
    for p=1:Lh
%-------------------------以下是C的求法
        k=-Lh/D:Lh/D;
        n=find(p-i==k*D);%只有一个成立
        if n%为真
           fai=D-1;
        else
           fai=-1;
        end
        
        if i==p
           C(i,p)=fai/D;
        else
           C(i,p)=fai*sin((p-i)*pi/D)/((p-i)*pi);
        end
%-------------------------以下是A的求法 
        if p==i
            A(i,p)=1;
        else
            A(i,p)=sin(wp*(p-i))/(wp*(p-i));
        end    
    end    
%-------------------------以下是b的求法
    if Th==i
        b(i,1)=1;
    else
        b(i,1)=sin(wp*(Th-i))/(wp*(Th-i));%96*1的列向量
    end
end
%%-----------------------以下是求理想分析原型滤波器
h=inv(A+C)*b;%96*1的列向量