exam5_8.m

来自「TI的DSP C55X的应用程序」· M 代码 · 共 15 行

%
% exam5_8.m - M-file for Example 5-8
%
omegac = 0.4;	 % cutoff frequency
L = 61;		 % filter order L = 61
M = (L-1)/2;     % M
l = 0:2*M;       % coefficient index l
h = omegac/pi*sinc(omegac*(l-M)/pi);  % compute coefficients
wn = 0.54 -0.46*cos(2*pi*l/(L-1));    % hamming window
hwn = h.*wn;     % windowing
omega = -pi:2*pi/200:pi;   % frequency range
Hd = freqz(hwn,1,omega);   % frequency response  
plot(omega,abs(Hd)), title('Magnitude Response'),...
    xlabel('Frequency'), ylabel('Magnitude'), grid;
axis([-3.14159 3.14159 0 1.2]);