Frequency Response of Nonlinear Systems: Notions, Analysis, and Graphical Representation
Pith reviewed 2026-05-10 12:45 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
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
axioms (1)
- domain assumption The invariance principle of Isidori and Byrnes applies to the steady-state behavior of the nonlinear system under periodic excitation.
Reference graph
Works this paper leans on
-
[1]
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]
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
work page 1979
-
[3]
J. Baillieul and T. Samad, Eds.,Encyclopedia of Systems and Control. Springer, 2021
work page 2021
-
[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
work page 2024
-
[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
work page 2025
-
[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
work page 2015
-
[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
work page 2022
-
[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
work page 2019
-
[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
work page 2022
-
[10]
D. Del Vecchio and R. M. Murray,Biomolecular Feedback Systems. Princeton Univ. Press, 2015
work page 2015
- [11]
-
[12]
C. R. Paul,Analysis of Multiconductor Transmission Lines. John Wiley & Sons, 2007
work page 2007
-
[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
work page 2009
-
[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]
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
work page 2024
-
[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
work page 2018
-
[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
work page 2015
-
[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
work page 2019
-
[19]
S. Skogestad and I. Postlethwaite,Multivariable Feedback Control: Analysis and Design. John Wiley & Sons, 2005
work page 2005
-
[20]
G. F. Franklin, J. D. Powell, and A. Emami-Naeini,Feedback Control of Dynamic Systems. Pearson, 2025
work page 2025
-
[21]
J. C. Doyle, B. A. Francis, and A. R. Tannenbaum,Feedback Control Theory. Courier Corp., 2013
work page 2013
-
[22]
H. W. Bode,Network Analysis and Feedback Amplifier Design. D. Van Nostrand Co., Inc., 1945
work page 1945
-
[23]
H. Nyquist, “Regeneration theory,”Bell Syst. Tech. J., vol. 11, no. 1, pp. 126–147, 1932
work page 1932
-
[24]
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
work page 2003
-
[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
work page 1980
-
[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
work page 1992
-
[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
work page 2011
-
[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
work page 1990
-
[29]
Isidori,Nonlinear Control Systems
A. Isidori,Nonlinear Control Systems. Springer, 1995
work page 1995
-
[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
work page 2008
-
[31]
Isidori,Lectures in Feedback Design for Multivariable Systems
A. Isidori,Lectures in Feedback Design for Multivariable Systems. Springer, 2017
work page 2017
-
[32]
D. A. George, “Continuous nonlinear systems,” Massachusetts Institute of Technology. Research Laboratory of Electronics, Tech. Rep., 1959
work page 1959
-
[33]
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
work page 2018
-
[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
work page 2003
-
[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
work page 1992
-
[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
work page 2010
-
[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
work page 2007
-
[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
work page 1979
-
[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
work page 1994
-
[40]
W. J. Rugh,Nonlinear System Theory. Johns Hopkins Univ. Press, 1981
work page 1981
-
[41]
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
work page 1966
-
[42]
C. Chen, D. Zhao, W. Chen, S. Z. Khong, and L. Qiu, “Phase of nonlinear systems,”arXiv preprint arXiv:2012.00692, 2021
-
[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
work page 2025
-
[44]
F. W. King,Hilbert Transforms. Cambridge Univ. Press,, 2009
work page 2009
-
[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
work page 2023
-
[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
work page 2022
-
[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)
work page 2026
-
[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
work page 2006
-
[49]
Cveticanin,Strong Nonlinear Oscillators: Analytical Solutions
L. Cveticanin,Strong Nonlinear Oscillators: Analytical Solutions. Springer, 2018
work page 2018
-
[50]
Carr,Applications of Centre Manifold Theory
J. Carr,Applications of Centre Manifold Theory. Springer, 2012
work page 2012
-
[51]
Hartman,Ordinary Differential Equations, 2nd ed
P. Hartman,Ordinary Differential Equations, 2nd ed. Birkh ¨auser, 1982
work page 1982
-
[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
work page 2025
-
[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
work page 2009
-
[54]
F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark,NIST Handbook of Mathematical Functions. Cambridge Univ. Press, 2010
work page 2010
-
[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
work page 1972
-
[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
work page 1980
-
[57]
B. Brogliato, R. Lozano, B. Maschke, and O. Egeland,Dissipative Systems Analysis and Control: Theory and Applications. Springer, 2020
work page 2020
-
[58]
H. K. Khalil,Nonlinear Control. Pearson Education, 2015
work page 2015
-
[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
work page 2026
-
[60]
W. Chen, D. Wang, and L. Qiu, “Small phase theorem,” inEncycl. Syst. Control, 2021, pp. 2082–2086
work page 2021
-
[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
work page 2025
-
[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
work page 2017
-
[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
work page 2025
-
[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
work page 2024
-
[65]
V . Ku ˇcera, “The matrix equationAX+XB=C,”SIAM J. Appl. Math., vol. 26, no. 1, pp. 15–25, 1974
work page 1974
-
[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
work page 1981
-
[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
work page 2006
-
[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 ...
work page 1964
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.