Main Content

Polynomial Stability Test

Use Schur-Cohn algorithm to determine whether all roots of input polynomial are inside unit circle

Library

Math Functions / Polynomial Functions

dsppolyfun

  • Polynomial Stability Test block

Description

The Polynomial Stability Test block uses the Schur-Cohn algorithm to determine whether all roots of a polynomial are within the unit circle.

y = all(abs(roots(u)) < 1)								% Equivalent MATLAB code

Each column of the M-by-N input matrix u contains M coefficients from a distinct polynomial,

f(x)=u1xM1+u2xM2++uM

arranged in order of descending exponents, u1, u2, ..., uM. The polynomial has order M-1 and positive integer exponents.

Inputs to the block represent the polynomial coefficients as shown in the previous equation. The block always treats length-M unoriented vector input as an M-by-1 matrix.

The output is a 1-by-N matrix with each column containing the value 1 or 0. The value 1 indicates that the polynomial in the corresponding column of the input is stable; that is, the magnitudes of all solutions to f(x) = 0 are less than 1. The value 0 indicates that the polynomial in the corresponding column of the input might be unstable; that is, the magnitude of at least one solution to f(x) = 0 is greater than or equal to 1.

Applications

This block is most commonly used to check the pole locations of the denominator polynomial, A(z), of a transfer function, H(z).

H(z)=B(z)A(z)=b1+b2z1++bmz(m1)a1+a2z1++anz(n1)

The poles are the n-1 roots of the denominator polynomial, A(z). When any poles are located outside the unit circle, the transfer function H(z) is unstable. As is typical in DSP applications, the transfer function above is specified in descending powers of z-1 rather than z.

Supported Data Types

  • Double-precision floating point

  • Single-precision floating point

  • Boolean — Block outputs are always Boolean.

See Also

Least Squares Polynomial FitDSP System Toolbox
Polynomial EvaluationDSP System Toolbox
polyfitMATLAB

Extended Capabilities

Version History

Introduced before R2006a