%例程4－10    Ⅳ型濾波器的幅頻函數

function [Hr,w,d,L]=hr_type4(h);         
%計算所設計的4型濾波器的振幅響應
%Hr=振幅響應
%b=4型濾波器的系數
%L=Hr的階次
%h=4型濾波器的單位沖擊響應

M=length(h);
L=M/2;
d= 2*[h(L:-1:1)];
n=[1:1:L];
n=n-0.5;
w=[0:1:500]'*2*pi/500;
Hr=sin(w*n)*d';