examp48.m

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

clc,echo on 
%EXAMPLE 48 
F=(0:399)/800; W=2*pi*F;                       % Frequency array 
h1=freqz([1 0],[1 -0.9],W);                    % Evaluate H(F) for filter 1 
h2=freqz([1 0],[1 0.9],W);                     % For filter 2
h3=freqz([0 -0.9],[1 -0.9],W);                 % For filter 3 
h4=freqz([1 0 0],[1 0.4 0.8],W);               % For filter 4 
h5=freqz([1 0.1 1],[1 0.1 0.8],W);             % For filter 5 
h6=freqz([1 zeros(1,7) -1],[1 zeros(1,8)],W);  % For filter 6 
plot(F,abs(h1));grid,pause(2)                  % Plot magnitude for filter 1 
plot(F,angle(h1));grid,pause(2)                % And phase (radians) 
plot(F,abs(h2));grid,pause(2)                  % Plot magnitude for filter 2 
plot(F,angle(h2));grid,pause(2)                % And phase (radians) 
plot(F,abs(h3));grid,pause(2)                  % Plot magnitude for filter 3 
plot(F,angle(h3));grid,pause(2)                % And phase (radians) 
plot(F,abs(h4));grid,pause(2)                  % Plot magnitude for filter 4 
plot(F,angle(h4));grid,pause(2)                % And phase (radians) 
plot(F,abs(h5));grid,pause(2)                  % Plot magnitude for filter 5 
plot(F,angle(h5));grid,pause(2)                % And phase (radians) 
plot(F,abs(h6));grid,pause(2)                  % Plot magnitude for filter 6 
plot(F,angle(h6));grid                         % And phase (radians) 
echo off %end of example