examp49.m

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

clc,echo on 
%EXAMPLE 49 
n=0:8; x=[3 3 3 2 2 2 1 1 1]; 
F=(0:199)/200; W=2*pi*F;
X=freqz(x,[zeros(1,4) 1 zeros(1,4)],W); 
Y=freqz(x,[1 zeros(1,8)],W);
plot(F,abs(X));grid,pause(2)
plot(F,abs(Y));grid,pause(2)
plot(F,180*angle(X)/pi);grid,pause(2)
plot(F,180*angle(Y)/pi);grid,pause(2)
phi=unwrap(angle(Y)-angle(X));grid,plot(F,phi)
echo off %end of example