examp67.m

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

clc,echo on 
%EXAMPLE 67 
Fc=0.2; n=-7:7; lp1=2*Fc*sinc(2*n*Fc);          % LPF cutoff=Fc 
hw=lp1 .* window('blac', 15);                   % Windowed response
tfplot('z',hw,1,[0 0.5],0,1);pause(2)
hhp=udelta(n)-hw;                               % LP2HP transformation
tfplot('z',hhp,1,[0 0.5],0,1);pause(2)
S=100; fp=20;fs=30;                             % BP specifications
F0=0.5*(fp+fs)/S; Fc2=0.5*(fs-fp)/S;            % Find F0 and 
n=-9.5:9.5; lp2=2*Fc*sinc(2*n*Fc2);             % LPF cutoff=Fc2 
hw2=lp2 .* window('vonh', 20);hw2=hw2/sum(hw2); % Windowed normalized response 
hbp2=2*cos(2*pi*F0*n) .* hw2;                   % LP2BP transformation 
tfplot('z',hbp2,1,[0 0.5],0,1);
echo off %end of example