ideal_bs.m

matlab exprience,M文件,实验代码

function hd=ideal_bs(Wcl,Wch,N)
%compute the ideal bandstop filter unit pulse respondense hd(n)
%wcl:low cutoff frequency
%wch:high cutoff frequency
%N:window length
%hd:unit pulse respondense

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