ideal_lp.m

来自「基于matlab的数字信号处理经典实验代码及所有完整实验报告」· M 代码 · 共 12 行

function hd=ideal_lp(wc,M);
% Ideal Lowpass filter computation 
% ---------------------
% [hd]=ideal_lp(wc,M)
% hd=ideal impulse response between 0 to M-1
% wc=cutoff frequency in radians
% M=length of the ideal filter
%
alpha=(M-1)/2;
n=[0:(M-1)];
m=n-alpha+eps;
hd=sin(wc*m)./(pi*m);