You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I'm not aware of any plans along these lines, but it would be great if you wanted to submit some code implementing this functionality. Information on how to contribute is here.
We do compute what this paper calls the S-based disk margin in stability_margins; it's called "Stability margin" in our docs.
Slycot already has a function to compute the transfer function matrix L-infinity norm (which is the H-infinity norm for stable systems), so the hard part to implement disk-margins is done.
The SISO disk margins are straightforward, and I've got prototype code to find the disk margin, and the associated destabilising complex perturbation and destabilising LTI system.
The MIMO version of disk margins needs structured singular values, or at least I'm pretty sure it does. SLICOT has AB13MD, which computes an upper bound on structure singular values, and if we want MIMO disk margins, I guess we'd start there. For now, one can get conservative upper bounds on (two of) these using linfnorm on the sensitivity and complementary sensitivity transfer functions.
Activity
murrayrm commentedon Apr 21, 2022
I'm not aware of any plans along these lines, but it would be great if you wanted to submit some code implementing this functionality. Information on how to contribute is here.
sylvain-bougnoux commentedon Apr 22, 2022
I wish I could do it, currently I am not familiar enough with the algorithm. FTTB, one can find an implementation here (in julia)
https://github.com/JuliaControl/RobustAndOptimalControl.jl/blob/master/src/diskmargin.jl
roryyorke commentedon Apr 24, 2022
We do compute what this paper calls the S-based disk margin in
stability_margins; it's called "Stability margin" in our docs.Slycot already has a function to compute the transfer function matrix L-infinity norm (which is the H-infinity norm for stable systems), so the hard part to implement disk-margins is done.
roryyorke commentedon Apr 27, 2022
See also the published form of the disk margin paper, and Matlab's diskmargin command. One of the paper's authors, the late Andrew Packard, was an author of Matlab's Robust Control Toolbox.
The SISO disk margins are straightforward, and I've got prototype code to find the disk margin, and the associated destabilising complex perturbation and destabilising LTI system.
The MIMO version of disk margins needs structured singular values, or at least I'm pretty sure it does. SLICOT has AB13MD, which computes an upper bound on structure singular values, and if we want MIMO disk margins, I guess we'd start there. For now, one can get conservative upper bounds on (two of) these using
linfnormon the sensitivity and complementary sensitivity transfer functions.