function [Hr,w,c,L]=hr_type3(h);
%線性相位濾波器類型3，計算濾波器振幅響應
%Hr振幅響應；w在[0,2pi]之間計算Hr的500個頻率點；
%c為3型低通濾波器系數；L為Hr的階次；h為3型低通濾波器的脈沖響應。
M=length(h);L=(M-1)/2;
c=[ 2*h(L+1:-1:1)];   %1乘(L+1)行向量
n=[0:1:L];                %(L+1)乘1列向量
w=[0:1:500]'*2*pi/500;
Hr=sin(w*n)*c';