function [Hr,w,d,L]=hr_type4(h);
%線性相位濾波器類型4，計(jì)算濾波器振幅響應(yīng)
%Hr振幅響應(yīng)；w在[0,2pi]之間計(jì)算Hr的500個(gè)頻率點(diǎn)；
%d為4型低通濾波器系數(shù)；L為Hr的階次；h為4型低通濾波器的脈沖響應(yīng)。
M=length(h);L=M/2;
d=2*[h(L:-1:1)];            %1乘(L+1)行向量
n=[1:1:L];n=n-0.5;          %(L+1)乘1列向量
w=[0:1:500]'*2*pi/500;
Hr=sin(w*n)*d';