## factorize

Purpose

Produce a spectral factorization of a scalar polynomial.

Synopsis

h=factorize(g)

Description

This is a standard system theory function providing spectral factorization of scalar polynomials without zeros on the imaginary axis:

$$g(\omega^{2})=h(j\omega)\divideontimes h(-j\omega),$$

where h is a Hurwitz polynomial and h(0) > 0. The function is useful when handling complex frequency scaling factors obtained via curve ﬁtting, such as in iqc delay.

Example

>> factorize([1 0 4])