function [Hr,w,b,L] = Hr_Type2(h);
% 計算第二種低通濾波器設計的振幅響應Hr(w) 
% -----------------------------------------------------------
%  [Hr,w,b,L] = Hr_Type2(h)
%  Hr = 振幅響應
%  w = 在[0 pi] 區間內計算Hr 的頻率點
%  b = 第二種低通濾波器的系數
%  L = Hr的階次
%  h = 第二種低通濾波器的頻率響應
% 
 M = length(h);
 L = M/2;
 b = 2*[h(L:-1:1)];
 n = [1:1:L]; n = n-0.5;
 w = [0:1:500]'*pi/500;
Hr = cos(w*n)*b';