hd=ideal_lp(wc,N);

function hd=ideal_lp(wc,N);
alpha=(N-1)/2;		%Ӧλ
n=[0:1:(N-1)];
m=n-alpha+eps;		%epsȷĸǷ
hd=sin(wc*m)./(pi*m);