ideal_lp.m

新版Matlab 7辅助信号处理技术与应用

function hd=ideal_lp(Wc,N)
% compute the ideal lowpass 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=sin(Wc*m)./(pi*m);