poly_stable (MathScript RT Module Function)
Owning Class: linear systems
Requires: MathScript RT Module
Syntax
b = poly_stable(a)
Legacy Name: polystab
Description
Stabilizes a polynomial by ensuring that all roots of the polynomial a are located within the unit circle. If a root is located outside of the unit circle, poly_stable uses a circle-based geometric reflection to reflect the root inside of the unit circle.
Inputs
| Name | Description |
|---|---|
| a | Specifies a polynomial. |
Outputs
| Name | Description |
|---|---|
| b | Returns a stabilized polynomial. |
Examples
a = rand(1, 6);
b = poly_stable(a);
ra = roots(a);
rb = roots(b);
plot(real(ra), imag(ra), 'ro', real(rb), imag(rb), 'b*')