ideal_bs.m

FIR滤波器设计

function hd=ideal_bs(wc1,wc2,N)

%compute the ideal bandstop fiter unit pulse respondence hd(n)

%wcl: low cutoff frequency

%wch: high cutoff frequency

%N: window length

%hd: unit pulse respondence

alpha= (N-1)/2;

n=0:1:N-1;

m=n-alpha+eps;

hd=[sin(wc1*m)+sin(pi*m)-sin(wc2*m)]./(pi*m);