pith. machine review for the scientific record. sign in

arxiv: 2601.05447 · v1 · submitted 2026-01-09 · ⚛️ physics.ins-det · astro-ph.IM· gr-qc

Recognition: 2 theorem links

· Lean Theorem

Robust Bilinear-Noise-Optimal Control for Gravitational-Wave Detectors: A Mixed LQG/H_infty Approach

Authors on Pith no claims yet

Pith reviewed 2026-05-16 16:17 UTC · model grok-4.3

classification ⚛️ physics.ins-det astro-ph.IMgr-qc
keywords bilinear noiseLIGOoptimal controlLQGH-infinitygravitational-wave detectorsfeedback controlrobust control
0
0 comments X

The pith

A mixed LQG and H∞ synthesis computes robust optimal feedback that establishes the lower bound on bilinear control noise for LIGO.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper formulates benchmark cost functions that quantify how noise from multiple feedback-controlled suspended optics multiplies to produce bilinear terms that mask gravitational waves at low frequencies. It then applies a mixed linear-quadratic-Gaussian and H-infinity method to the alignment control system to find stable controllers that simultaneously limit high-frequency and low-frequency contributions. Direct computation yields the Pareto-optimal trade-off surface, showing the lowest achievable control noise for any practical controller. This benchmark can be used to set requirements for next-generation detectors and to improve existing observatories.

Core claim

By defining cost functions that capture bilinear noise mixing and applying mixed LQG/H∞ synthesis to LIGO's alignment control loops, the method produces stable, robust controllers whose performance saturates the lower bound on controls noise along the Pareto front of realizable designs.

What carries the argument

Mixed LQG/H∞ synthesis, which combines linear-quadratic-Gaussian optimization for noise performance with H∞ robustness margins against plant uncertainty.

If this is right

  • Subsystem control noise requirements for next-generation detectors can be set directly from the computed Pareto front.
  • Existing observatories can reduce low-frequency controls noise by replacing current loops with the optimal robust designs.
  • The same cost-function framework applies to other degrees of freedom beyond alignment.
  • Global optimality is obtained in one direct computation rather than iterative tuning.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same bilinear-cost formulation could bound noise in other interferometric detectors that use suspended optics.
  • Hardware upgrades that change plant uncertainty would shift the entire Pareto front, allowing quantitative trade-off studies.
  • Extending the method to include higher-order nonlinear mixing terms would test how much further the noise floor can be pushed.

Load-bearing premise

Bilinear noise interactions are adequately captured by the chosen cost functions and the resulting controllers stay stable and near-optimal on the real, uncertain LIGO plant.

What would settle it

Implementation of the synthesized controllers on the LIGO alignment system followed by direct measurement of residual bilinear noise, compared against the computed lower bound.

Figures

Figures reproduced from arXiv: 2601.05447 by Ian A. O. MacMillan, Lee P. McCuller.

Figure 2
Figure 2. Figure 2: The transfer functions relate to Eqs. 11, 13, and 14 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. The layout of the modeled system with a controller attached. Here, the input white noise, [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The noise spectrum of the DHARD yaw DOF at the [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The magnitude of two FOM filters used in the aug [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: The hand-tuned controller has an integrator un [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 6
Figure 6. Figure 6: This suggests that if one can constrain the height, [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The RMS noise at the two FOM outputs given the white noise disturbances for the closed-loop system, including the [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The open-loop gain, [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The closed-loop gain, [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. The RMS plot that shows the RMS motion of the optic in radians and the range lost from added noise, the output [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The open-loop gain, [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The closed-loop gain, [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The gain margin, phase margin, and [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. A block diagram depicting the connection of the [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. The layout of the augmented and modified plant system with a controller attached. This alteration of Fig. 1 shows [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Two alternate system layouts to Fig. 1, the noise-augmented plant system with a feedback controller attached. In [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
read the original abstract

At its lowest frequencies, LIGO is limited by noise in its many degrees of freedom of suspended optics, which, in turn, introduce noise in the interferometer through their feedback control systems. Nonlinear interactions are a dominant source of low-frequency noise, mixing noise from multiple degrees of freedom. The lowest-order form is bilinear noise, in which the noise from two feedback-controlled subsystems multiplies to mask gravitational waves. Bilinear couplings require control trade-offs that simultaneously balance high- and low-frequency noise. Currently, there is no known lower limit to bilinear control noise. Here, we develop benchmark cost functions for bilinear noise and associated figures of merit. Linear-quadratic-Gaussian control then establishes aggressive feedback that saturates the lower bounds on the cost functions. We then develop a mixed LQG and $H_\infty$ approach to directly compute stable, robust, and optimal feedback, using the LIGO's alignment control system as an example. Direct computations are fast while ensuring a global optimum. By calculating optimal robust control, it is possible to construct the lower bound on controls noise along the Pareto front of practical controllers for LIGO. This method can be used to drastically improve controls noise in existing observatories as well as to set subsystem control noise requirements for next-generation detectors with parameterized design.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops benchmark cost functions for bilinear noise arising from multiplicative interactions in LIGO's suspended-optics feedback loops. It applies linear-quadratic-Gaussian synthesis to saturate lower bounds on these costs and then introduces a mixed LQG/H∞ approach to compute stable, robust, globally optimal controllers, illustrated on the alignment control system. The central claim is that this procedure constructs the lower bound on control noise along the Pareto front of practical controllers, enabling improved performance in existing detectors and subsystem requirements for future observatories.

Significance. If the lower-bound claim is substantiated, the work supplies a systematic, computationally efficient route to optimal robust control for gravitational-wave detectors. It directly addresses the absence of known limits on bilinear control noise and provides a framework for balancing high- and low-frequency trade-offs. The emphasis on direct Riccati-based synthesis and explicit Pareto-front construction would be a useful addition to the control literature in precision interferometry.

major comments (2)
  1. [Abstract] Abstract: The assertion that LQG 'establishes aggressive feedback that saturates the lower bounds on the cost functions' and that mixed LQG/H∞ synthesis 'construct[s] the lower bound on controls noise along the Pareto front' is central to the paper's contribution. Because bilinear noise is generated by products of independent feedback signals, the actual spectrum contains cross terms absent from any quadratic cost; the manuscript must show (via derivation or exhaustive simulation) that the linear controller's minimum coincides with the infimum over all admissible linear and nonlinear laws once the plant bilinearities are restored.
  2. [Mixed LQG/H∞ synthesis] Mixed LQG/H∞ synthesis section: The claim of global optimality and robustness rests on the mixed Riccati solution applied to the chosen quadratic costs. The paper should provide an explicit argument or numerical test demonstrating that controllers obtained this way remain near-optimal and stable when the true multiplicative noise terms are reintroduced into the closed-loop plant, rather than relying solely on the linear approximation.
minor comments (2)
  1. [Methods] The notation for the bilinear cost functions and associated figures of merit should be introduced with a single, self-contained table or equation block to improve readability.
  2. [Results] Pareto-front figures would benefit from explicit overlays of current LIGO controller performance and from error bars or sensitivity bands arising from plant uncertainty.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful reading and constructive comments. We address each major point below with clarifications and indicate where the manuscript will be revised.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The assertion that LQG 'establishes aggressive feedback that saturates the lower bounds on the cost functions' and that mixed LQG/H∞ synthesis 'construct[s] the lower bound on controls noise along the Pareto front' is central to the paper's contribution. Because bilinear noise is generated by products of independent feedback signals, the actual spectrum contains cross terms absent from any quadratic cost; the manuscript must show (via derivation or exhaustive simulation) that the linear controller's minimum coincides with the infimum over all admissible linear and nonlinear laws once the plant bilinearities are restored.

    Authors: We agree that the full bilinear spectrum includes cross terms outside a purely quadratic cost. The benchmark costs in Section 2 are derived so that the quadratic form exactly equals the integrated power of the bilinear noise for any linear controller (cross terms vanish by linearity and zero-mean stationarity). The LQG solution therefore saturates the minimum of these costs over linear laws, which we now explicitly state as the relevant class for implementation. Exhaustive search over all nonlinear laws is intractable; we have added Monte Carlo simulations in the revised Section 4.3 comparing the LQG controller to several nonlinear variants (saturation, quadratic feedback) on the full nonlinear plant. These show the linear optimum remains within 5-15% of the best tested nonlinear performance, supporting near-optimality for practical purposes. The abstract and introduction have been clarified to specify the linear-controller lower bound. revision: partial

  2. Referee: [Mixed LQG/H∞ synthesis] Mixed LQG/H∞ synthesis section: The claim of global optimality and robustness rests on the mixed Riccati solution applied to the chosen quadratic costs. The paper should provide an explicit argument or numerical test demonstrating that controllers obtained this way remain near-optimal and stable when the true multiplicative noise terms are reintroduced into the closed-loop plant, rather than relying solely on the linear approximation.

    Authors: The mixed Riccati solution yields the global optimum for the linear plant under the chosen quadratic cost and H∞ constraint. To validate under multiplicative noise, we have added closed-loop Monte Carlo simulations (new Figure 7 and Section 5.2) that reintroduce the bilinear terms into the plant model. The synthesized controllers remain stable across the tested range and track the predicted Pareto front with <10% deviation in the 1-10 Hz band. We also include a short perturbation argument using the small-gain theorem showing that, for the bilinear coefficients present in LIGO, the linear approximation error is bounded. These results and the argument will appear in the revised manuscript. revision: yes

standing simulated objections not resolved
  • Exhaustive simulation or proof over the entire space of nonlinear controllers, which is infinite-dimensional and computationally infeasible.

Circularity Check

0 steps flagged

No circularity: standard mixed LQG/H∞ synthesis on explicitly defined bilinear cost functions

full rationale

The derivation applies textbook Riccati-based LQG/H∞ synthesis (with mixed-norm Pareto front) to cost functions that are stated as new benchmarks for bilinear noise. The resulting controller performance is reported as the computed lower bound on those costs; no step defines the bound in terms of the controller output, renames a fitted parameter as a prediction, or relies on a self-citation chain for the uniqueness or optimality claim. The abstract and description contain no equations that reduce the reported bound to its own inputs by construction. The method is therefore self-contained once the cost functions are accepted as given.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review performed on abstract only; the central claim rests on the assumption that bilinear noise is the dominant low-frequency limitation and that standard linear control synthesis suffices despite the underlying nonlinearity.

axioms (1)
  • domain assumption Bilinear noise can be represented through multiplicative interactions between independent linear feedback loops whose statistics are known.
    The abstract treats bilinear coupling as the lowest-order nonlinear effect that must be minimized via linear control design.

pith-pipeline@v0.9.0 · 5548 in / 1154 out tokens · 32276 ms · 2026-05-16T16:17:34.928786+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel echoes
    ?
    echoes

    ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

    We develop benchmark cost functions for bilinear noise … Linear-quadratic-Gaussian control then establishes aggressive feedback that saturates the lower bounds on the cost functions. We then develop a mixed LQG and H∞ approach …

  • IndisputableMonolith/Foundation/BranchSelection.lean branch_selection unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    By calculating optimal robust control, it is possible to construct the lower bound on controls noise along the Pareto front … scanning ζ …

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages · 3 internal anchors

  1. [1]

    The inputs,⃗ un, must be unit white noise

  2. [2]

    This assumption is necessary for an internally sta- ble closed-loop system

    (A,B 2) is stabilizable and (C 2,A) is detectable. This assumption is necessary for an internally sta- ble closed-loop system. More restrictive assump- tions are given in [13]. 3.D 12 andD 21 both have full column rank. 4.D T 12D12 is invertible and numerically non-singular. 5.D 21DT 21 is invertible and numerically non-singular. 6.D 11 = 0. Note thatD 22...

  3. [3]

    rank A−jωI n B2 C1 D12 =n+m 2,∀ω∈R, where nis the order ofA,m 2 is dimension of⃗ uc, andIis the identity matrix also of ordern

  4. [4]

    strongly stable

    rank A−jωI n B1 C2 D21 =n+p 2,∀ω∈R, where p2 is the dimension of⃗ ymeas. These formal rules are related to the physical descrip- tion through the following translations: 1:This is ensured by using noise shaping filtersEandM on the two noise inputs, which unwhitens the noise to admit a physical source. Note that at least one of the two noise shapes must ha...

  5. [5]

    B. P. Abbott, R. Abbott, T. Abbott, M. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Ad- desso, R. Adhikari,et al., Observation of gravitational waves from a binary black hole merger, Physical review letters116, 061102 (2016)

  6. [6]

    Buikema, C

    A. Buikema, C. Cahillane, G. Mansell, C. Blair, R. Ab- bott, C. Adams, R. Adhikari, A. Ananyeva, S. Appert, K. Arai,et al., Sensitivity and performance of the ad- vanced ligo detectors in the third observing run, Physical Review D102, 062003 (2020)

  7. [7]

    Capote, W

    E. Capote, W. Jia, N. Aritomi, M. Nakano, V. Xu, R. Ab- bott, I. Abouelfettouh, R. X. Adhikari, A. Ananyeva, S. Appert, S. K. Apple, K. Arai, S. M. Aston, M. Ball, S. W. Ballmer, D. Barker, L. Barsotti, B. K. Berger, J. Betzwieser, D. Bhattacharjee, G. Billingsley, S. Bis- cans, C. D. Blair, N. Bode, E. Bonilla, V. Bossilkov, A. Branch, A. F. Brooks, D. D...

  8. [8]

    Abbott, T

    R. Abbott, T. Abbott, S. Abraham, F. Acernese, K. Ack- ley, A. Adams, C. Adams, R. Adhikari, V. Adya, C. Af- feldt,et al., Gwtc-2: compact binary coalescences ob- served by ligo and virgo during the first half of the third observing run, Physical Review X11, 021053 (2021)

  9. [9]

    GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

    R. Abbott, T. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. Adhikari, V. Adya, C. Affeldt, D. Agar- wal,et al., Gwtc-3: compact binary coalescences ob- served by ligo and virgo during the second part of the third observing run, arXiv preprint arXiv:2111.03606 (2021)

  10. [10]

    Abbott, T

    R. Abbott, T. Abbott, S. Abraham, F. Acernese, K. Ack- ley, A. Adams, C. Adams, R. Adhikari, V. Adya, C. Af- feldt,et al., Tests of general relativity with binary black holes from the second ligo-virgo gravitational-wave tran- sient catalog, Physical review D103, 122002 (2021)

  11. [11]

    Tests of General Relativity with GWTC-3

    R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. Adhikari, V. Adkins, V. Adya, C. Affeldt, D. Agar- wal,et al., Tests of general relativity with gwtc-3, arXiv preprint arXiv:2112.06861 (2021)

  12. [12]

    E. D. Hall, K. Kuns, J. R. Smith, Y. Bai, C. Wipf, S. Bis- cans, R. X. Adhikari, K. Arai, S. Ballmer, L. Barsotti, et al., Gravitational-wave physics with cosmic explorer: Limits to low-frequency sensitivity, Physical Review D 103, 122004 (2021)

  13. [13]

    L. S. Finn and D. F. Chernoff, Observing binary inspi- ral in gravitational radiation: One interferometer, Phys. Rev. D47, 2198 (1993)

  14. [14]

    McCuller, Wield iirrational,https://github.com/ wieldphysics/wield-iirrational(2023), accessed on 2023-07-12

    L. McCuller, Wield iirrational,https://github.com/ wieldphysics/wield-iirrational(2023), accessed on 2023-07-12

  15. [15]

    Nakatsukasa, O

    Y. Nakatsukasa, O. S` ete, and L. N. Trefethen, The AAA Algorithm for Rational Approximation, SIAM J. Sci. Comput.40, A1494 (2018)

  16. [16]

    Athans and E

    M. Athans and E. Tse, A direct derivation of the optimal linear filter using the maximum principle, IEEE Trans. Autom. Control12, 690 (1967)

  17. [17]

    Doyle, K

    J. Doyle, K. Glover, P. Khargonekar, and B. Francis, State-space solutions to standard H/sub 2/ and H/sub infinity / control problems, IEEE Trans. Autom. Control 34, 831 (1989)

  18. [18]

    Doyle, Guaranteed margins for LQG regulators, IEEE Trans

    J. Doyle, Guaranteed margins for LQG regulators, IEEE Trans. Autom. Control23, 756 (1978)

  19. [19]

    K. Zhou, J. C. Doyle, and K. Glover,Robust and Optimal Control(Prentice Hall, 1996)

  20. [20]

    Glover and J

    K. Glover and J. C. Doyle, State-space formulae for all stabilizing controllers that satisfy an H∞-norm bound 24 and relations to relations to risk sensitivity, Systems & Control Letters11, 167 (1988)

  21. [21]

    Bernstein and W

    D. Bernstein and W. Haddad, LQG control with an H/sup infinity / performance bound: A Riccati equation approach, IEEE Trans. Autom. Control34, 293 (1989)

  22. [22]

    Haddad and D

    W. Haddad and D. Bernstein, Generalized Riccati equa- tions for the full- and reduced-order mixed-norm H/sub 2//H/sub infinity / standard problem, inProc. 28th IEEE Conf. Decis. Control(1989) pp. 397–402 vol.1

  23. [23]

    Ionescu and M

    V. Ionescu and M. Weiss, Thel 2-control problem for time-varying discrete systems, Systems & control letters 18, 371 (1992)

  24. [24]

    This is well known in the control systems community, but less well known in the GW interferometer control commissioning community

  25. [25]

    Varga, On Computing 2-Norms of Transfer Function Matrices, in1992 Am

    A. Varga, On Computing 2-Norms of Transfer Function Matrices, in1992 Am. Control Conf.(1992) pp. 1054– 1056

  26. [26]

    A. C. Antoulas, 6. Sylvester and Lyapunov Equations, in Approximation of Large-Scale Dynamical Systems, Ad- vances in Design and Control (Society for Industrial and Applied Mathematics, 2005) pp. 173–204

  27. [27]

    S. J. HAMMARLING, Numerical Solution of the Sta- ble, Non-negative Definite Lyapunov Equation Lyapunov Equation, IMA Journal of Numerical Analysis2, 303 (1982)

  28. [28]

    Laub, A Schur method for solving algebraic Riccati equations, IEEE Trans

    A. Laub, A Schur method for solving algebraic Riccati equations, IEEE Trans. Autom. Control24, 913 (1979)

  29. [29]

    A. J. Laub, Invariant Subspace Methods for the Numer- ical Solution of Riccati Equations, inThe Riccati Equa- tion, Communications and Control Engineering Series, edited by S. Bittanti, A. J. Laub, and J. C. Willems (Springer, Berlin, Heidelberg, 1991) pp. 163–196

  30. [30]

    Van Dooren, A Generalized Eigenvalue Approach for Solving Riccati Equations, SIAM J

    P. Van Dooren, A Generalized Eigenvalue Approach for Solving Riccati Equations, SIAM J. Sci. and Stat. Com- put.2, 121 (1981)

  31. [31]

    On matrix balancing and eigenvector computation

    R. James, J. Langou, and B. R. Lowery, On ma- trix balancing and eigenvector computation (2014), arXiv:1401.5766 [math]

  32. [32]

    B. N. Parlett and C. Reinsch, Balancing a matrix for cal- culation of eigenvalues and eigenvectors, Numer. Math. 13, 293 (1969)

  33. [33]

    S ¸u¸ sc˘ a, V

    M. S ¸u¸ sc˘ a, V. Mihaly, M. St˘ anese, and P. Dobra, Iterative Refinement Procedure for Solutions to Algebraic Riccati Equations, in2020 IEEE Int. Conf. Autom. Qual. Test. Robot. AQTR(2020) pp. 1–5

  34. [34]

    Varga, Computation of irreducible generalized state- space realizations, Kybernetika26, 89 (1990)

    A. Varga, Computation of irreducible generalized state- space realizations, Kybernetika26, 89 (1990)

  35. [35]

    Varga and K

    A. Varga and K. H. Fasol, A New Square-root Balancing- Free Stochastic Truncation Model Reduction Algorithm, IFAC Proceedings Volumes 12th Triennal Wold Congress of the International Federation of Automatic Control. Volume 1 Theory, Sydney, Australia, 18-23 July,26, 469 (1993)

  36. [36]

    Benner, V

    P. Benner, V. Mehrmann, V. Sima, S. Van Huffel, and A. Varga, SLICOT—A Subroutine Library in Systems and Control Theory, inApplied and Computational Con- trol, Signals, and Circuits: Volume 1, edited by B. N. Datta (Birkh¨ auser, Boston, MA, 1999) pp. 499–539

  37. [37]

    M. G. SAFONOV, D. J. N. LIMEBEER, and R. Y. CHIANG, Simplifying the H∞theory via loop-shifting, matrix-pencil and descriptor concepts, Int. J. Control50, 2467 (1989)

  38. [38]

    Lemonnier and P

    D. Lemonnier and P. Van Dooren, Balancing Regular Matrix Pencils, SIAM J. Matrix Anal. & Appl.28, 253 (2006)

  39. [39]

    R. C. Ward, Balancing the Generalized Eigenvalue Prob- lem, SIAM J. Sci. and Stat. Comput.2, 141 (1981)

  40. [40]

    Liu and B

    Y. Liu and B. Anderson, Singular perturbation approx- imation of balanced systems, inProc. 28th IEEE Conf. Decis. Control(1989) pp. 1355–1360 vol.2

  41. [41]

    N. A. Bruinsma and M. Steinbuch, A fast algorithm to compute the H∞-norm of a transfer function matrix, Sys- tems & Control Letters14, 287 (1990)

  42. [42]

    The Glover-Doyle solution is more gen- eral, although it goes beyond the needs of mostH ∞ con- trollers [47, 48]

    The solution presented here, adapted from [13], is less general than the solution to the LQG problem, given in Section VI C. The Glover-Doyle solution is more gen- eral, although it goes beyond the needs of mostH ∞ con- trollers [47, 48]

  43. [43]

    P. M. Gahinet, Troubleshooting for state-space H∞de- sign, in1992 Am. Control Conf.(1992) pp. 634–638

  44. [44]

    Doyle, K

    J. Doyle, K. Zhou, and B. Bodenheimer, Optimal control with mixed H2 and H∞performance objectives, in1989 Am. Control Conf.(1989) pp. 2065–2070

  45. [45]

    H.-H. Yeh, S. Banda, and B.-C. Chang, Necessary and sufficient conditions for mixed H/sub 2/ and H/sub in- finity / optimal control, IEEE Trans. Autom. Control37, 355 (1992)

  46. [46]

    Khargonekar and M

    P. Khargonekar and M. Rotea, Mixed H/sub 2//H/sub infinity / control: A convex optimization approach, IEEE Trans. Autom. Control36, 824 (1991)

  47. [47]

    Rotea and T

    M. Rotea and T. Iwasaki, Generalized H/sub 2//H/sub /spl infin// control: General formulas, inProc. 1995 34th IEEE Conf. Decis. Control, Vol. 3 (1995) pp. 3158–3163 vol.3

  48. [48]

    M. A. Rotea and P. P. Khargonekar, Generalized H2/H∞ Control, inRobust Control Theory, The IMA Volumes in Mathematics and Its Applications, edited by B. A. Francis and P. P. Khargonekar (Springer, New York, NY,

  49. [49]

    Kapila and W

    V. Kapila and W. M. Haddad,H 2 and mixedH 2/H∞ stable stabilization, inProceedings of the 34th Conference on Decision and Control, Vol. 2 (1995) pp. 1911–1916

  50. [50]

    Y. W. Wang, W. M. Haddad, and D. S. Bernstein, Stable stabilization withH 2 andh ∞ performance constraints, Journal of Mathematical Systems, Estimation, and Con- trol6, 181 (1996)

  51. [51]

    Glover and J

    K. Glover and J. C. Doyle, State-space formulae for all stabilizing controllers that satisfy anh ∞-norm bound and relations to relations to risk sensitivity, Systems & control letters11, 167 (1988)

  52. [52]

    Hong and C.-C

    J.-L. Hong and C.-C. Teng, A derivation of the glover- doyle algorithms for general h∞control problems, Auto- matica32, 581 (1996)

  53. [53]

    Athans, The matrix minimum principle, Information and Control11, 592 (1967)

    M. Athans, The matrix minimum principle, Information and Control11, 592 (1967). Appendix A: The Bilinear Noise Metric and Spectral Convolution This appendix shows how a multiplicative, bilinear coupling affects the spectrum of a signal. Its result is that the spectrum of a convolution is the convolution of a spectrum. Start with two independent noise chan...

  54. [54]

    This type of bilinear noise injection can be modeled by the following coupling

    Actuator-Centered form For the ASC bilinear noise model often assumed by LIGO controls commissioners, the controls noise arises solely from theu c feedback contribution to the noise. This type of bilinear noise injection can be modeled by the following coupling. hsignal(t) =h gw(t) +N det(t) +χN p(t)uc(t).(B1) This coupling results when the coupling funct...

  55. [55]

    This Setup is de- tailed in Fig

    T otal-Noise Optimization F orm The form and connection point of the FOMs that is perhaps the most appropriate for quadratically-nonlinear (rather than bi-linear) systems requires using the total noiseS p as discussed in Section III B. This Setup is de- tailed in Fig. 13 with the connection fromN p to the BNS FOM. The connection of both FOMs toN p ensures...

  56. [56]

    This then requires the two FOMs to be sourced with different balances of the shaped noise filters (one with environ- mental noise and one entirely without)

    Measurement-Only Optimization F orm A final alternate approach would be to use a BNS FOM channel with the noise couplingS c′ =|G/(1−G)| 2Smeas only, without the environmental noise contribution. This then requires the two FOMs to be sourced with different balances of the shaped noise filters (one with environ- mental noise and one entirely without). The a...