p10_9.m

来自「数字信号处理实验指导书MATLAB版程序电子版」· M 代码 · 共 20 行

% Program P10_9
% Design of L-th Band FIR Filter Using the 
% Windowed Fourier Series Approach
%
clf;
K = 11;
n = -K:K;
% Generate the truncated impulse response of  
% the ideal lowpass filter
b = sinc(n/2)/2;
% Generate the window sequence
win = hamming(23);
% Generate the coefficients of the windowed filter
fil = b.*win';
c = fil/sum(fil);
% Plot the gain response of the windowed filter
[h,w] = freqz(c,1,256);
g = 20*log10(abs(h));
plot(w/pi,g);axis([0 1 -90 10]); grid
xlabel('\omega/