%《數字信號處理教程——MATLAB釋義與實現》第七章演示程序fgp731
% 理想脈沖序列與矩形窗的時域相乘和頻域卷積
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
clear
wc=1;N=15;alpha=(N-1)/2;                        % 給出理想頻率特性參數
ns=-19;nf=29;n=[ns:nf]+1e-10;                   % 給出序列長度區間
hd=sin(n-alpha)./(n-alpha)/pi;                  % 理想脈沖序列
subplot(3,2,1),stem(n,hd,'.')                   % 畫脈沖序列圖
win=[zeros(1,0-ns),ones(1,N),zeros(1,nf-N+1)];  % 給出矩形窗函數序列
subplot(3,2,3),stem(n,win,'.')                  % 畫矩形窗圖        
w=[-1:0.01:1]*pi+1e-10;                         % 給出頻率自變量向量
h=seqmult(hd,n,win,n);                   % 理想脈沖序列與矩形窗相乘，得到實際的脈沖響應
Win=sin(0.5*N*w)./(0.5*N*w);                    % 矩形窗函數的頻譜
N1=round(1/0.01/pi);
Hd=[zeros(1,100-N1),ones(1,2*N1+1),zeros(1,100-N1)];    % 理想濾波器頻譜
H=conv(Hd,Win);                          % 理想濾波器頻譜與矩形窗頻譜卷積，得到實際頻譜
subplot(3,2,5),stem(n,h,'.')                    % 繪實際的脈沖響應圖
subplot(3,2,2),plot([-pi,-1,-1,1,1,pi],[0,0,1,1,0,0]) % 繪理想濾波器頻譜圖
subplot(3,2,4),plot(w,Win);                     % 繪矩形窗函數的頻譜圖
subplot(3,2,6),plot(w,H(floor(0.5*length(w):1.5*length(w)-1)));  % 繪實際濾波器頻譜圖
set(gcf,'color','w')                        % 置圖形背景色為白