4-15.m

1.本书包含的实例所使用的源程序清单

%例程4－15   Ⅲ、Ⅳ型理想高通滤波器的单位冲激响应 计算

function hd=ideal_hp2(Wc,N)
%compute the ideal highpass fiter unit pulse respondence hd(n)
%wc: cutoff frequency
%N: window length
%hd: unit pulse respondence

alpha=(N-1)/2;
n=0:1:N-1;
m=n-alpha+eps;
hd=[cos(pi*m)-cos(Wc*m)]./(pi*m);