%《數字信號處理教程——MATLAB釋義與實現》第四章例5.5.1程序hc551
% 理想低通濾波器的脈沖響應及因果性
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
wc=1;n0=6;n=[-10:30]+1e-5;          % 給出偏移的時間序列
h=sin(n-n0)./(n-n0)/pi;        % 計算偏移6的理想低通濾波器的脈沖響應序列
subplot(2,1,1),stem(n,h)        % 畫出脈沖序列
title('理想因果低通濾波器的脈沖響應n0=6')
line([-10,35],[0,0])
line([0,0],[-0.1,0.4])
h1=sin(n-15)./(n-15)/pi;    % 計算偏移15的理想低通濾波器的脈沖響應序列
subplot(2,1,2),stem(n,h1)      % 畫出脈沖序列
title('理想因果低通濾波器的脈沖響應n0=15')
line([-10,35],[0,0])
line([0,0],[-0.1,0.4])
set(gcf,'color','w');                % 置圖形背景色為白