ideallp.m

来自「这是数字信号处理——matlab释义与实现一书的源代码」· M 代码 · 共 16 行

function hd = ideallp(wc,N);
%《数字信号处理教程——MATLAB释义与实现》第七章例7.3.2程序hc732
% 理想低通滤波器的脉冲响应子程序
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% hd = ideallp(wc,N)
% ----------------------------------------------
% hd = 点0 到 N-1之间的理想脉冲响应
% wc = 截止频率(弧度) 
% N = 理想滤波器的长度
%
tao = (N-1)/2;              % 理想脉冲响应的对称中心位置
n = [0: (N-1)];             % 设定脉冲响应长度
m = n - tao + eps;	        % 加一个小数以避免零作除数
hd = sin(wc*m) ./ (pi*m);   % 理想脉冲响应