Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Turn off warn_infinite for phase_crossover_frequencies() calculation #1106

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 2 commits into from
Jan 30, 2025
Merged

Turn off warn_infinite for phase_crossover_frequencies() calculation #1106

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
6fcf741
turn off warn_infinite for phase_crossover_frequencies() calculation
murrayrm Jan 29, 2025
d9e44d2
address @slivingston review comments
murrayrm Jan 30, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 9 additions & 4 deletions control/margins.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -81,11 +81,16 @@ def _poly_iw_sqr(pol_iw):

def _poly_iw_real_crossing(num_iw, den_iw, epsw):
# Return w where imag(H(iw)) == 0

# Compute the imaginary part of H = (num.r + j num.i)/(den.r + j den.i)
test_w = np.polysub(np.polymul(num_iw.imag, den_iw.real),
np.polymul(num_iw.real, den_iw.imag))

# Find the real-valued w > 0 where imag(H(iw)) = 0
w = np.roots(test_w)
w = np.real(w[np.isreal(w)])
w = w[w >= epsw]

return w


Expand Down Expand Up @@ -471,7 +476,7 @@ def phase_crossover_frequencies(sys):
omega : ndarray
1d array of (non-negative) frequencies where Nyquist plot
intersects the real axis
gain : ndarray
gains : ndarray
1d array of corresponding gains

Examples
Expand All @@ -493,13 +498,13 @@ def phase_crossover_frequencies(sys):
omega = _poly_iw_real_crossing(num_iw, den_iw, 0.)

# using real() to avoid rounding errors and results like 1+0j
gain = np.real(evalfr(sys, 1J * omega))
gains = np.real(sys(omega * 1j, warn_infinite=False))
else:
zargs = _poly_z_invz(sys)
z, omega = _poly_z_real_crossing(*zargs, epsw=0.)
gain = np.real(evalfr(sys, z))
gains = np.real(sys(z, warn_infinite=False))

return omega, gain
return omega, gains


def margin(*args):
Expand Down
13 changes: 6 additions & 7 deletions control/tests/margin_test.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,12 +12,9 @@
from numpy import inf, nan
from numpy.testing import assert_allclose

from control.frdata import FrequencyResponseData
from control.margins import (margin, phase_crossover_frequencies,
stability_margins)
from control.statesp import StateSpace
from control.xferfcn import TransferFunction
from control.exception import ControlMIMONotImplemented
from control import ControlMIMONotImplemented, FrequencyResponseData, \
StateSpace, TransferFunction, margin, phase_crossover_frequencies, \
stability_margins

s = TransferFunction.s

Expand Down Expand Up @@ -111,15 +108,17 @@ def test_margin_3input(tsys):
out = margin((mag, phase*180/np.pi, omega_))
assert_allclose(out, np.array(refout)[[0, 1, 3, 4]], atol=1.5e-3)


@pytest.mark.parametrize(
'tfargs, omega_ref, gain_ref',
[(([1], [1, 2, 3, 4]), [1.7325, 0.], [-0.5, 0.25]),
(([1], [1, 1]), [0.], [1.]),
(([2], [1, 3, 3, 1]), [1.732, 0.], [-0.25, 2.]),
((np.array([3, 11, 3]) * 1e-4, [1., -2.7145, 2.4562, -0.7408], .1),
[1.6235, 0.], [-0.28598, 1.88889]),
(([200.0], [1.0, 21.0, 20.0, 0.0]),
[4.47213595, 0], [-0.47619048, inf]),
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Warnings do not cause the test to fail, so this case does not provide a regression test. To fix this, add the decorator @pytest.mark.filterwarnings("error") as done in another test in margin_test.py.

])
@pytest.mark.filterwarnings("error")
def test_phase_crossover_frequencies(tfargs, omega_ref, gain_ref):
"""Test phase_crossover_frequencies() function"""
sys = TransferFunction(*tfargs)
Expand Down
Loading