pith. sign in

arxiv: 2604.13842 · v1 · submitted 2026-04-15 · 📡 eess.SY · cs.SY

Frequency Response of Nonlinear Systems: Notions, Analysis, and Graphical Representation

Pith reviewed 2026-05-10 12:45 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords nonlinear systemsfrequency responsephasor formdistortion functionBode diagramssteady-state analysiscontrol design
0
0 comments X

The pith

Nonlinear systems under periodic excitations have a frequency response defined as a phasor with gain, phase and distortion functions.

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

The paper shows that frequency response remains definable for nonlinear systems driven by nonlinear periodic inputs through a complex phasor representation of the steady-state output. Alongside gain and phase, a distortion function is introduced to quantify how the output deviates from a pure sinusoid. This matters because it creates a complete picture of nonlinear effects in the frequency domain. If the claim holds, engineers gain graphical tools to plot these quantities against both frequency and amplitude, supporting analysis and design of systems that must handle nonlinear behavior.

Core claim

For systems under nonlinear periodic excitations, the frequency response is defined as a complex-valued function in phasor form. This is completed by gain and phase functions plus a distortion function that quantifies the alteration introduced by the system in the steady-state output, allowing the full characterization to be shown in diagrams over input frequency and amplitude.

What carries the argument

the phasor representation of the steady-state output completed by a distortion function

If this is right

  • Gain, phase, and distortion can be plotted as functions of both input frequency and amplitude.
  • The resulting diagrams support performance analysis of nonlinear systems.
  • The loop-shaping design problem can be stated directly for nonlinear systems using the new representation.

Where Pith is reading between the lines

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

  • The distortion function may serve as a practical metric to decide when a nonlinear system is close enough to linear for standard design tools to apply.
  • This approach could extend frequency-domain methods to systems with harmonic-rich inputs without needing separate time-domain analysis.
  • Numerical computation of the three functions for specific nonlinear models would allow direct comparison of different system designs in the same format.

Load-bearing premise

That the steady-state output of a nonlinear system under periodic inputs settles into a pattern fully captured by one complex phasor plus a separate distortion term.

What would settle it

A concrete nonlinear system driven by a periodic input whose steady-state output waveform cannot be summarized by a single phasor value and one distortion metric would disprove the proposed characterization.

Figures

Figures reproduced from arXiv: 2604.13842 by Alessio Moreschini, Matteo Scandella.

Figure 5
Figure 5. Figure 5: Specifically, Figure 5a shows that the ω-gain is consistent with that of the linearized system until the effect of saturation becomes visible and dominates the linear effect. This is consistent with the shape of the ω-radius in Figure 5c, which proves that, at high amplitudes, the system bends the shape of the input, achieving high nonlinearity at high frequencies. Finally, Figure 5b confirms that the phas… view at source ↗
read the original abstract

The invariance principle, through which the steady-state behavior of nonlinear systems was introduced by Isidori and Byrnes, is leveraged in this article to bring forth a unifying characterization of the frequency response of nonlinear systems. We show that, for systems under nonlinear periodic excitations, the frequency response can still be defined as a complex-valued function in a phasor form. However, together with suitable notions of gain and phase functions, we show the existence of another function that completes the frequency response and allows quantifying the distortion introduced by the system in the steady-state output. This nonlinear characterization enabled the representation over input frequency and amplitude of the gain, phase, and distortion produced by the system, via a nonlinear enhancement of the Bode diagrams. This graphical representation of the frequency response is well-suited to performance analysis of a nonlinear system and, furthermore, allows for the formulation of the loop-shaping problem for nonlinear systems.

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 paper leverages the Isidori-Byrnes invariance principle to characterize the frequency response of nonlinear systems under nonlinear periodic excitations. It defines this response as a complex-valued phasor function, accompanied by gain and phase functions, and introduces a distortion function to quantify nonlinearity-induced effects in the steady-state output. This framework enables nonlinear Bode diagrams that represent gain, phase, and distortion over ranges of input frequency and amplitude, supporting performance analysis and the formulation of loop-shaping problems for nonlinear systems.

Significance. If the central derivations hold with the stated conditions, the work offers a principled extension of linear frequency-domain tools to nonlinear systems. The distortion function and resulting graphical representations provide concrete tools for analyzing steady-state behavior without full time-domain simulation, and the loop-shaping formulation could influence nonlinear controller design. The grounding in the established invariance principle is a clear strength, as is the explicit construction of amplitude- and frequency-dependent plots.

major comments (2)
  1. [Section 4] The main theorem (Section 4): the existence of the phasor representation and distortion function is asserted to follow directly from the invariance principle applied to the steady-state output, but the manuscript must state the precise hypotheses on the system class (e.g., relative degree, stability of the zero dynamics, and the form of the periodic excitation) that guarantee a unique attractive periodic response; without these, the claim that the phasor and distortion are well-defined for general nonlinear periodic inputs is not yet load-bearing.
  2. [Section 4] Definition of the distortion function (Eq. (12) or equivalent): the paper should demonstrate that this function is independent of the choice of coordinates on the invariant manifold and reduces to zero for linear systems; the current sketch leaves open whether the distortion term is uniquely determined or could be absorbed into a redefinition of the phasor.
minor comments (2)
  1. [Section 3] Notation: the symbols for the gain, phase, and distortion functions should be introduced with a single consistent table or list early in the paper to avoid later confusion when they appear in the nonlinear Bode plots.
  2. [Abstract] The abstract states the existence result but does not mention the required system hypotheses; adding one sentence on the class of systems would improve readability for readers outside the immediate subfield.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment, the recommendation for minor revision, and the constructive comments on Section 4. We address each major comment below and will incorporate the requested clarifications and explicit statements into the revised manuscript.

read point-by-point responses
  1. Referee: [Section 4] The main theorem (Section 4): the existence of the phasor representation and distortion function is asserted to follow directly from the invariance principle applied to the steady-state output, but the manuscript must state the precise hypotheses on the system class (e.g., relative degree, stability of the zero dynamics, and the form of the periodic excitation) that guarantee a unique attractive periodic response; without these, the claim that the phasor and distortion are well-defined for general nonlinear periodic inputs is not yet load-bearing.

    Authors: We appreciate the referee highlighting the need for explicit hypotheses. The derivations in Section 4 rely on the standard conditions from the Isidori-Byrnes invariance principle under which a unique attractive periodic response exists. In the revised manuscript we will state these hypotheses explicitly at the beginning of Section 4: the plant has well-defined relative degree, the zero dynamics are asymptotically stable, and the periodic excitation is generated by a neutrally stable exosystem. These conditions guarantee the existence of a unique invariant manifold and therefore well-defined phasor and distortion functions. We will also add a direct reference to the relevant theorem in Isidori and Byrnes. revision: yes

  2. Referee: [Section 4] Definition of the distortion function (Eq. (12) or equivalent): the paper should demonstrate that this function is independent of the choice of coordinates on the invariant manifold and reduces to zero for linear systems; the current sketch leaves open whether the distortion term is uniquely determined or could be absorbed into a redefinition of the phasor.

    Authors: We agree that the distortion function requires an explicit demonstration of coordinate independence and uniqueness. The distortion is defined intrinsically as the L2-norm of the residual between the steady-state output waveform and its fundamental harmonic component (extracted via the phasor). Because this definition operates directly on the scalar output signal, it is independent of any particular coordinate chart on the invariant manifold. For linear systems the steady-state output is purely sinusoidal at the driving frequency, so the residual vanishes and the distortion is identically zero. In the revision we will add a short remark (or appendix paragraph) showing that the phasor is the unique Fourier coefficient of the fundamental frequency on the manifold; the distortion is then the orthogonal complement in L2 and cannot be absorbed into a redefinition of the phasor. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's derivation begins from the external Isidori-Byrnes invariance principle (cited as prior literature, not self-authored) to establish existence of a unique attractive steady-state periodic output under nonlinear periodic forcing. From this foundation the authors introduce phasor-form frequency response together with explicit gain, phase, and distortion functions, then construct the associated nonlinear Bode representation. No equation or claim reduces by construction to a fitted parameter, self-definition, or load-bearing self-citation; the central steps remain independent once the cited invariance result is granted. The construction is therefore self-contained against an external benchmark.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The characterization depends on the invariance principle as background and on the assumption that periodic excitations admit a well-defined steady-state phasor representation; no free parameters or new invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption The invariance principle of Isidori and Byrnes applies to the steady-state behavior of the nonlinear system under periodic excitation.
    Invoked in the abstract as the foundation for defining the frequency response.

pith-pipeline@v0.9.0 · 5453 in / 1332 out tokens · 19094 ms · 2026-05-10T12:45:35.460338+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

68 extracted references · 68 canonical work pages

  1. [1]

    Bode-like diagrams illustrating the performance at vari- ous angular frequencyϖand amplitudea u are reported in Figure 5

    In this example, we consider the case withγ= 5,k on = 1, koff = 0.5,p tot = 20,¯u= 100.5,K= 0.05andK u = 0.1. Bode-like diagrams illustrating the performance at vari- ous angular frequencyϖand amplitudea u are reported in Figure 5. Specifically, Figure 5a shows that theω-gain is consistent with that of the linearized system until the effect of saturation ...

  2. [2]

    The development of frequency-response methods in automatic control [perspectives],

    A. MacFarlane, “The development of frequency-response methods in automatic control [perspectives],”IEEE Trans. Autom. Control, vol. 24, no. 2, pp. 250–265, 1979

  3. [3]

    Baillieul and T

    J. Baillieul and T. Samad, Eds.,Encyclopedia of Systems and Control. Springer, 2021

  4. [4]

    Gain and phase: Decentralized stability conditions for power electronics-dominated power systems,

    L. Huang, D. Wang, X. Wang, H. Xin, P. Ju, K. H. Johansson, and F. D¨orfler, “Gain and phase: Decentralized stability conditions for power electronics-dominated power systems,”IEEE Trans. Power Syst., vol. 39, no. 6, pp. 7240–7256, 2024

  5. [5]

    Interharmonic power–a new concept for power system oscillation source location,

    W. Xu, J. Yong, H. J. Marquez, and C. Li, “Interharmonic power–a new concept for power system oscillation source location,”IEEE Trans. Power Syst., 2025

  6. [6]

    M. J. Gibbard, P. Pourbeik, and D. J. V owles,Small-Signal Stability, Control and Dynamic Performance of Power Systems. Univ. Adelaide Press, 2015

  7. [7]

    Power system stability with a high penetration of inverter-based resources,

    Y . Gu and T. C. Green, “Power system stability with a high penetration of inverter-based resources,”Proc. IEEE, vol. 111, no. 7, pp. 832–853, 2022

  8. [8]

    Hard limits and performance tradeoffs in a class of antithetic integral feedback networks,

    N. Olsman, A.-A. Baetica, F. Xiao, Y . P. Leong, R. M. Murray, and J. C. Doyle, “Hard limits and performance tradeoffs in a class of antithetic integral feedback networks,”Cell Systems, vol. 9, no. 1, pp. 49–63, 2019

  9. [9]

    Cybergenetics: Theory and applications of genetic control systems,

    M. H. Khammash, “Cybergenetics: Theory and applications of genetic control systems,”Proc. IEEE, vol. 110, no. 5, pp. 631–658, 2022

  10. [10]

    Del Vecchio and R

    D. Del Vecchio and R. M. Murray,Biomolecular Feedback Systems. Princeton Univ. Press, 2015

  11. [11]

    Keener and J

    J. Keener and J. Sneyd,Mathematical Physiology. Springer, 2025

  12. [12]

    C. R. Paul,Analysis of Multiconductor Transmission Lines. John Wiley & Sons, 2007

  13. [13]

    Complex brain networks: graph theoretical analysis of structural and functional systems,

    E. Bullmore and O. Sporns, “Complex brain networks: graph theoretical analysis of structural and functional systems,”Nat. Rev. Neurosci., vol. 10, no. 3, pp. 186–198, 2009

  14. [14]

    Frequency-domain bounds for the multiconductor telegrapher’s equa- tion,

    D. Selvaratnam, A. Moreschini, A. Das, T. Parisini, and H. Sandberg, “Frequency-domain bounds for the multiconductor telegrapher’s equa- tion,”arXiv preprint arXiv:2504.01599, 2025

  15. [15]

    Memristors on ‘edge of chaos’,

    L. O. Chua, “Memristors on ‘edge of chaos’,”Nat. Rev. Electr. Eng., vol. 1, no. 9, pp. 614–627, 2024

  16. [16]

    Graph signal processing: Overview, challenges, and ap- plications,

    A. Ortega, P. Frossard, J. Kova ˇcevi´c, J. M. F. Moura, and P. Van- dergheynst, “Graph signal processing: Overview, challenges, and ap- plications,”Proc. IEEE, vol. 106, no. 5, pp. 808–828, 2018

  17. [17]

    Boashash,Time-Frequency Signal Analysis and Processing: A Com- prehensive Reference

    B. Boashash,Time-Frequency Signal Analysis and Processing: A Com- prehensive Reference. Acad. Press, 2015

  18. [18]

    Learning graphs from data: A signal representation perspective,

    X. Dong, D. Thanou, M. Rabbat, and P. Frossard, “Learning graphs from data: A signal representation perspective,”IEEE Signal Process. Mag., vol. 36, no. 3, pp. 44–63, 2019

  19. [19]

    Skogestad and I

    S. Skogestad and I. Postlethwaite,Multivariable Feedback Control: Analysis and Design. John Wiley & Sons, 2005

  20. [20]

    G. F. Franklin, J. D. Powell, and A. Emami-Naeini,Feedback Control of Dynamic Systems. Pearson, 2025

  21. [21]

    J. C. Doyle, B. A. Francis, and A. R. Tannenbaum,Feedback Control Theory. Courier Corp., 2013

  22. [22]

    H. W. Bode,Network Analysis and Feedback Amplifier Design. D. Van Nostrand Co., Inc., 1945

  23. [23]

    Regeneration theory,

    H. Nyquist, “Regeneration theory,”Bell Syst. Tech. J., vol. 11, no. 1, pp. 126–147, 1932

  24. [24]

    Feedback and optimal sensitivity: Model reference trans- formations, multiplicative seminorms, and approximate inverses,

    G. Zames, “Feedback and optimal sensitivity: Model reference trans- formations, multiplicative seminorms, and approximate inverses,”IEEE Trans. Autom. Control, vol. 26, no. 2, pp. 301–320, 2003

  25. [25]

    Feedback stabilization of linear dynamical plants with uncertainty in the gain factor,

    A. Tannenbaum, “Feedback stabilization of linear dynamical plants with uncertainty in the gain factor,”Int. J. Control, vol. 32, no. 1, pp. 1–16, 1980. 16

  26. [26]

    A loop-shaping design procedure using H∞ synthesis,

    D. McFarlane and K. Glover, “A loop-shaping design procedure using H∞ synthesis,”IEEE Trans. Autom. Control, vol. 37, no. 6, pp. 759– 769, 1992

  27. [27]

    Kramers–Kronig, Bode, and the meaning of zero,

    J. Bechhoefer, “Kramers–Kronig, Bode, and the meaning of zero,”Am. J. Phys., vol. 79, no. 10, pp. 1053–1059, 2011

  28. [28]

    Output regulation of nonlinear systems,

    A. Isidori and C. I. Byrnes, “Output regulation of nonlinear systems,” IEEE Trans. Autom. Control, vol. 35, no. 2, pp. 131–140, 1990

  29. [29]

    Isidori,Nonlinear Control Systems

    A. Isidori,Nonlinear Control Systems. Springer, 1995

  30. [30]

    Steady-state behaviors in nonlinear systems with an application to robust disturbance rejection,

    A. Isidori and C. I. Byrnes, “Steady-state behaviors in nonlinear systems with an application to robust disturbance rejection,”Annu. Rev. Control, vol. 32, no. 1, pp. 1–16, 2008

  31. [31]

    Isidori,Lectures in Feedback Design for Multivariable Systems

    A. Isidori,Lectures in Feedback Design for Multivariable Systems. Springer, 2017

  32. [32]

    Continuous nonlinear systems,

    D. A. George, “Continuous nonlinear systems,” Massachusetts Institute of Technology. Research Laboratory of Electronics, Tech. Rep., 1959

  33. [33]

    The analysis of nonlinear systems in the frequency domain using nonlinear output frequency response functions,

    R. S. Bayma, Y . Zhu, and Z.-Q. Lang, “The analysis of nonlinear systems in the frequency domain using nonlinear output frequency response functions,”Automatica, vol. 94, pp. 452–457, 2018

  34. [34]

    Fading memory and the problem of approximat- ing nonlinear operators with V olterra series,

    S. Boyd and L. Chua, “Fading memory and the problem of approximat- ing nonlinear operators with V olterra series,”IEEE Trans. Circuits Syst., vol. 32, no. 11, pp. 1150–1161, 2003

  35. [35]

    Disturbance attenuation andH ∞-control via measurement feedback in nonlinear systems,

    A. Isidori and A. Astolfi, “Disturbance attenuation andH ∞-control via measurement feedback in nonlinear systems,”IEEE Trans. Autom. Control, vol. 37, no. 9, pp. 1283–1293, 1992

  36. [36]

    Model reduction by moment matching for linear and nonlinear systems,

    A. Astolfi, “Model reduction by moment matching for linear and nonlinear systems,”IEEE Trans. Autom. Control, vol. 55, no. 10, pp. 2321–2336, 2010

  37. [37]

    Frequency response functions for nonlinear convergent systems,

    A. Pavlov, N. van de Wouw, and H. Nijmeijer, “Frequency response functions for nonlinear convergent systems,”IEEE Trans. Autom. Con- trol, vol. 52, no. 6, pp. 1159–1165, 2007

  38. [38]

    Frequency domain analysis of nonlinear systems: general theory,

    L. O. Chua and C.-Y . Ng, “Frequency domain analysis of nonlinear systems: general theory,”IEE J. Electron. Circuits Syst., vol. 3, pp. 165– 185, 1979

  39. [39]

    Analysing non-linear systems in the fre- quency domain–II. the phase response,

    S. Billings and H. Zhang, “Analysing non-linear systems in the fre- quency domain–II. the phase response,”Mech. Syst. Signal Process., vol. 8, no. 1, pp. 45–62, 1994

  40. [40]

    W. J. Rugh,Nonlinear System Theory. Johns Hopkins Univ. Press, 1981

  41. [41]

    On the input-output stability of time-varying nonlinear feedback systems–Part II: Conditions involving circles in the frequency plane and sector nonlinearities,

    G. Zames, “On the input-output stability of time-varying nonlinear feedback systems–Part II: Conditions involving circles in the frequency plane and sector nonlinearities,”IEEE Trans. Autom. Control, vol. 11, no. 3, pp. 465–476, 1966

  42. [42]

    Phase of nonlinear systems,

    C. Chen, D. Zhao, W. Chen, S. Z. Khong, and L. Qiu, “Phase of nonlinear systems,”arXiv preprint arXiv:2012.00692, 2021

  43. [43]

    The singular angle of nonlinear systems,

    C. Chen, D. Zhao, and S. Z. Khong, “The singular angle of nonlinear systems,”Automatica, vol. 181, p. 112515, 2025

  44. [44]

    F. W. King,Hilbert Transforms. Cambridge Univ. Press,, 2009

  45. [45]

    Graphical nonlinear system analysis,

    T. Chaffey, F. Forni, and R. Sepulchre, “Graphical nonlinear system analysis,”IEEE Trans. Autom. Control, vol. 68, no. 10, pp. 6067–6081, 2023

  46. [46]

    Scaled relative graphs: Nonexpan- sive operators via 2D Euclidean geometry,

    E. K. Ryu, R. Hannah, and W. Yin, “Scaled relative graphs: Nonexpan- sive operators via 2D Euclidean geometry,”Math. Program., vol. 194, no. 1, pp. 569–619, 2022

  47. [47]

    On the frequency response and loop shaping for nonlinear systems,

    A. Moreschini, M. Scandella, and A. Astolfi, “On the frequency response and loop shaping for nonlinear systems,” inIFAC WC 2026, (Accepted)

  48. [48]

    Some asymptotic methods for strongly nonlinear equations,

    J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” Int. J. Mod. Phys. B, vol. 20, no. 10, pp. 1141–1199, 2006

  49. [49]

    Cveticanin,Strong Nonlinear Oscillators: Analytical Solutions

    L. Cveticanin,Strong Nonlinear Oscillators: Analytical Solutions. Springer, 2018

  50. [50]

    Carr,Applications of Centre Manifold Theory

    J. Carr,Applications of Centre Manifold Theory. Springer, 2012

  51. [51]

    Hartman,Ordinary Differential Equations, 2nd ed

    P. Hartman,Ordinary Differential Equations, 2nd ed. Birkh ¨auser, 1982

  52. [52]

    Amplitude response and square wave describ- ing functions,

    T. Chaffey and F. Forni, “Amplitude response and square wave describ- ing functions,”Eur. J. Control., p. 101310, 2025

  53. [53]

    Oscillator with fraction order restoring force,

    L. Cveticanin, “Oscillator with fraction order restoring force,”J. Sound Vib., vol. 320, no. 4, pp. 1064–1077, 2009

  54. [54]

    F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark,NIST Handbook of Mathematical Functions. Cambridge Univ. Press, 2010

  55. [55]

    Dissipative dynamical systems part I: General theory,

    J. C. Willems, “Dissipative dynamical systems part I: General theory,” Arch. Ration. Mech. Anal., vol. 45, no. 5, pp. 321–351, 1972

  56. [56]

    Dissipative dynamical systems: Basic input- output and state properties,

    D. J. Hill and P. J. Moylan, “Dissipative dynamical systems: Basic input- output and state properties,”J. Franklin Inst., vol. 309, no. 5, pp. 327– 357, 1980

  57. [57]

    Brogliato, R

    B. Brogliato, R. Lozano, B. Maschke, and O. Egeland,Dissipative Systems Analysis and Control: Theory and Applications. Springer, 2020

  58. [58]

    H. K. Khalil,Nonlinear Control. Pearson Education, 2015

  59. [59]

    Moment matching by kernel-based learning,

    A. Moreschini, M. Scandella, A. Astolfi, and T. Parisini, “Moment matching by kernel-based learning,”IEEE Trans. Autom. Control, vol. 71, no. 4, pp. 2123–2138, 2026

  60. [60]

    Small phase theorem,

    W. Chen, D. Wang, and L. Qiu, “Small phase theorem,” inEncycl. Syst. Control, 2021, pp. 2082–2086

  61. [61]

    A generalized passivity theory over abstract time domains,

    A. Moreschini, M. Bin, A. Astolfi, and T. Parisini, “A generalized passivity theory over abstract time domains,”IEEE Trans. Autom. Control, vol. 70, no. 1, pp. 2–17, 2025

  62. [62]

    Eigenvalues and poles of nonlinear systems: A geometric approach,

    A. Padoan and A. Astolfi, “Eigenvalues and poles of nonlinear systems: A geometric approach,” inProc. 56th IEEE Conf. Decis. Control (CDC), 2017, pp. 2575–2580

  63. [63]

    Closed-loop interpolation by moment matching for linear and nonlinear systems,

    A. Moreschini and A. Astolfi, “Closed-loop interpolation by moment matching for linear and nonlinear systems,”IEEE Trans. Autom. Control, vol. 70, no. 5, pp. 2918–2933, 2025

  64. [64]

    Interconnection-based model order reduc- tion - a survey,

    G. Scarciotti and A. Astolfi, “Interconnection-based model order reduc- tion - a survey,”Eur. J. Control, vol. 75, p. 100929, 2024

  65. [65]

    The matrix equationAX+XB=C,

    V . Ku ˇcera, “The matrix equationAX+XB=C,”SIAM J. Appl. Math., vol. 26, no. 1, pp. 15–25, 1974

  66. [66]

    Controllability, observability and the solution ofAX−XB=C,

    E. de Souza and S. Bhattacharyya, “Controllability, observability and the solution ofAX−XB=C,”Linear Algebra Appl., vol. 39, pp. 167–188, 1981

  67. [67]

    The polynomial solution to the Sylvester matrix equation,

    Q. Hu and D. Cheng, “The polynomial solution to the Sylvester matrix equation,”Appl. Math. Lett., vol. 19, no. 9, pp. 859–864, 2006

  68. [68]

    I. M. Gel’fand and G. E. Shilov,Generalized Functions, Volume 1: Properties and Operations. AMS Chelsea Publishing, 1964. APPENDIX Uniqueness of the solution of the Sylvester equation (4) is ensured by the fact thatAandSdo not share eigenvalues [64]. Controllability of the pair(A, B)and observability of the pair (S, L)imply that the solution is full rank ...