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);   % 理想脈沖響應