firminphase (Digital Filter Design Toolkit, MathScript Function)

Owning Class: singlerate

Syntax

b2 = firminphase(b1)

Description

Calculates the minimum phase spectral factor of a linear phase, finite impulse response (FIR) filter. The resulting spectral factor is also an FIR filter whose zeroes correspond to the zeroes of the original linear phase FIR filter inside or on the unit circle. If a zero is on the unit circle, the zero must be an even-multiplicity zero. In other words, the zero must occur an even number of times. The magnitude response of the spectral factor is the square root of that of the original FIR filter.

Examples

Inputs

Name	Description
b1	Specifies the coefficients of a linear phase filter with a nonnegative zero-phase response. b1 is a real vector.

Outputs

Name	Description
b2	Returns the minimum phase FIR spectral factor of b1. b2 is a real vector.

Examples

bmin = fircband(12, [0, 0.4, 0.5, 1], [1, 1, 0, 0], [1, 1], 'ww', 'minphase');bmax = fliplr(bmin);
b1 = conv(bmin, bmax);
b2 = firminphase(b1);
figure;
zplane(b1, 1);
figure;
zplane(b2, 1);