hr_type.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function [Hr,w,f,L]=hr_type(h);
%线性相位滤波器，计算滤波器振幅响应
%Hr振幅响应；w在[0,pi]之间计算Hr的500个频率点；
%a为低通滤波器系数；L为Hr的阶次；
%h为低通滤波器的脉冲响应，N为脉冲响应长度。
N=length(h);
if h(1)==h(N)
    M=rem(N,2);
    if  M==1
        [Hr,w,a,L]=hr_type1(h);
        f=a;
    else
        [Hr,w,b,L]=hr_type2(h);
        f=b;
    end
else
    M=rem(N,2);
    if  M==1
        [Hr,w,c,L]=hr_type3(h);
        f=c;
    else
        [Hr,w,d,L]=hr_type4(h);
        f=d;
    end 
end