4-11.m

Matlab数字信号处理与应用源代码

%例程4－11   Ⅰ、Ⅱ型理想低通滤波器的单位冲激响应 计算

function hd=ideal_lp1(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);