%《數字信號處理教程——MATLAB釋義與實現》第七章例7.4.3程序hc743
% 用fir2函數生成限帶波形
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
N=input('N=');M=N-1;  % 給定序列長度
f=[0,0.45,0.5,1];   % 給定預期頻譜特性
m=[0,1,0,0];
h=fir2(M,f,m,boxcar(N));      % 生成脈沖響應序列
[H,w]=freqz(h,1);    % 計算脈沖序列的頻譜
subplot(3,1,1),plot(f,m)    % 畫出設定頻譜形狀
title('設定頻譜'),
xlabel(''),ylabel('m')
subplot(3,1,2),stem([0:N-1],h,'.')    % 畫出生成的脈沖序列形狀
title('用fir2函數生成的脈沖序列'),
xlabel(''),ylabel('h(n)')
subplot(3,1,3),plot(w/pi,abs(H))    % 畫出實際頻譜形狀
title('實際頻譜'),
xlabel('\omega/\pi'),ylabel('abs(H)')
set(gcf,'color','w')                % 置圖形背景色為白