examp53.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 17 行

clc,echo on 
%EXAMPLE 53 
n=-20:20;                       % Create a DT index 
Fc=0.2; h=2*Fc*sinc(2*Fc*n);    % Create impulse response 
a=[zeros(1,20) 1 zeros(1,20)];  % Denominator array of DTFT 
f=(-200:199)/400;w=2*pi*f;      % Frequency array 
HI=(abs(f)<=0.2);               % Ideal Filter Spectrum,  $F_{C}=0.2$ 
wb=1-abs(n)/20;                 % Bartlett window 
wv=0.5+0.5*cos(pi*n/20);        % von Hann window 
wh=0.54+0.46*cos(pi*n/20);      % Hamming window 
HB=abs(freqz(h.*wb,a,w));       % Spectrum of Bartlett windowed h[n] 
HV=abs(freqz(h.*wv,a,w));       % Spectrum of von Hann windowed h[n] 
HH=abs(freqz(h.*wh,a,w));       % Spectrum of Hamming windowed h[n] 
plot(f,HI,f,HB,f,HV,f,HH),grid  % Overplot
echo off %end of example