zerophase - Zero-phase response of digital filter

Syntax

Description

[Hr,w] = zerophase(b,a) returns the zero-phase response Hr, and the frequency vector w (in radians/sample) at which Hr is computed, given a filter defined by numerator b and denominator a. For FIR filters where a=1, you can omit the value a from the command. The zero-phase response is evaluated at 512 equally spaced points on the upper half of the unit circle.

The zero-phase response, Hr(w), is related to the frequency response, H(w) by

where H(e^jω) is the frequency response, Hr(ω) is the zero-phase response and [[PHI1]](w) is the continuous phase.

Note   The zero-phase response is always real, but it is not the equivalent of the magnitude response. The former can be negative while the latter cannot be negative.

[Hr,w] = zerophase(b,a,nfft) returns the zero-phase response Hr and frequency vector w (radians/sample), using nfft frequency points on the upper half of the unit circle.

[Hr,w] = zerophase(b,a,nfft,'whole') returns the zero-phase response Hr and frequency vector w (radians/sample), using nfft frequency points around the whole unit circle.

[Hr,w] = zerophase(b,a,w) returns the zero-phase response Hr and frequency vector w (radians/sample) at frequencies in vector w.

[Hr,f] = zerophase(...,fs) returns the zero-phase response Hr and frequency vector f (Hz), using the sampling frequency fs (in Hz), to determine the frequency vector f (in Hz) at which Hr is computed.

[Hr,w,phi] = zerophase(...) returns the zero-phase response Hr, frequency vector w (rad/sample), and the continuous phase component, phi. (Note that this quantity is not equivalent to the phase response of the filter when the zero-phase response is negative.)

zerophase(...) with no output arguments, plots the zero-phase response versus frequency.

Examples

Example 1

Plot the zero-phase response of a constrained least squares FIR filter:

b=fircls1(54,.3,.02,.008);
zerophase(b);

Example 2

Plot the zero-phase response of an elliptic filter:

[b,a]=ellip(10,.5,20,.4);
zerophase(b,a,512,'whole');

See Also

freqs | freqz | fvtool | grpdelay | invfreqz | phasedelay | phasez

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS