SciPy

scipy.signal.unique_roots

scipy.signal.unique_roots(p, tol=0.001, rtype='min')[source]

Determine unique roots and their multiplicities from a list of roots.

Parameters:

p : array_like

The list of roots.

tol : float, optional

The tolerance for two roots to be considered equal. Default is 1e-3.

rtype : {‘max’, ‘min, ‘avg’}, optional

How to determine the returned root if multiple roots are within tol of each other.

  • ‘max’: pick the maximum of those roots.
  • ‘min’: pick the minimum of those roots.
  • ‘avg’: take the average of those roots.
Returns:

pout : ndarray

The list of unique roots, sorted from low to high.

mult : ndarray

The multiplicity of each root.

Notes

This utility function is not specific to roots but can be used for any sequence of values for which uniqueness and multiplicity has to be determined. For a more general routine, see numpy.unique.

Examples

>>> from scipy import signal
>>> vals = [0, 1.3, 1.31, 2.8, 1.25, 2.2, 10.3]
>>> uniq, mult = signal.unique_roots(vals, tol=2e-2, rtype='avg')

Check which roots have multiplicity larger than 1:

>>> uniq[mult > 1]
array([ 1.305])

Previous topic

scipy.signal.remez

Next topic

scipy.signal.residue