Nichols plot
The Nichols plot, also known as the Nichols chart, is a graphical tool in control systems engineering for analyzing the frequency response of open-loop transfer functions in feedback systems. It displays the magnitude of the transfer function G(j\omega) in decibels (dB) on the vertical axis against the phase angle \angle G(j\omega) in degrees on the horizontal axis, with frequency \omega serving as a parameter traced along the curve.[1] This representation reverses the coordinates of a standard Nyquist plot and uses a logarithmic scale for magnitude, enabling the overlay of contours for constant closed-loop magnitude (M-circles) and phase (N-circles) to assess stability and performance directly from open-loop data.[1][2] Named after American control engineer Nathaniel B. Nichols (1914–1997), the plot originated as a "decibel-phase-angle diagram" during World War II efforts to design reliable servomechanisms for radar and fire-control systems.[3][4] It was formally introduced in the 1947 book Theory of Servomechanisms by Hubert M. James, Nathaniel B. Nichols, and Ralph S. Phillips, part of the MIT Radiation Laboratory Series, where it extended earlier frequency-response methods by Hendrik Bode to facilitate practical stability analysis and parameter tuning.[3] Nichols, a key contributor to early control theory through his work at MIT and later at Raytheon, co-developed the technique to address challenges in optimizing servo performance under sinusoidal inputs, ensuring no encirclement of the critical point (-1, 0) in the complex plane for stability.[3][2] In practice, the Nichols plot is constructed by transferring points from a Bode plot—lifting magnitude and phase values at discrete frequencies onto the gain-phase plane—and is particularly effective for single-sheeted versions spanning phases from 0° to -360°.[1] It excels in determining gain and phase margins graphically: stability requires the plot to avoid the (0 dB, -180°) point, with margins read as horizontal and vertical distances from curve intersections near this critical location.[4][2] Compared to Nyquist diagrams, which use polar coordinates for encirclement counts, or Bode plots, which separate magnitude and phase, the Nichols plot combines their strengths by providing explicit frequency-domain insights into closed-loop behavior, such as resonance peaks and bandwidth, without complex-plane mapping.[1][2] This makes it valuable for controller design, including gain adjustments and applications in quantitative feedback theory (QFT), though modern software like MATLAB has reduced its manual use while preserving its interpretive power.[4][2]Definition and Background
Definition
The Nichols plot is a graphical representation in control engineering of the frequency response for the open-loop transfer function G(j\omega) of a linear time-invariant (LTI) system. It serves as a variant of the polar plot by positioning the magnitude $20 \log_{10} |G(j\omega)| in decibels (dB) along the vertical axis and the phase angle \angle G(j\omega) in degrees along the horizontal axis, with the latter typically ranging from 0° to -360° to capture phase lag.[1][5] This plotting method consolidates magnitude and phase data into one chart, allowing for efficient evaluation of stability margins and overall system performance in feedback control designs.[5] The logarithmic scaling of the magnitude axis accommodates broad dynamic ranges common in control systems, while the phase axis highlights shifts that indicate potential instability.[1] Originating in control theory, the Nichols plot originated as a tool for analyzing LTI systems and remains valuable for visualizing closed-loop responses through overlaid contours of constant magnitude and phase.[5][1]Historical Development
The Nichols plot was developed in the 1940s by Nathaniel B. Nichols, an American engineer working at the Massachusetts Institute of Technology (MIT), as part of efforts during World War II to analyze servo-mechanisms for military applications such as radar and fire control systems.[6] Nichols' work emerged from the MIT Radiation Laboratory, where interdisciplinary teams advanced feedback control theory to meet wartime demands for precise automatic control.[7] The plot's initial formal presentation appeared in 1947 within Chapter 4 of the influential book Theory of Servomechanisms, co-authored by Hubert M. James, Nathaniel B. Nichols, and Ralph S. Phillips as Volume 25 of the MIT Radiation Laboratory Series.[6] This publication integrated the Nichols plot into early feedback control research, building on foundational concepts like Hendrik Wade Bode's 1940s explorations of gain-phase relationships and Harry Nyquist's 1932 stability criterion, which highlighted limitations in polar and logarithmic representations that the new chart addressed.[8] In the 1950s and 1960s, the Nichols plot evolved alongside advances in control systems design, with the advent of digital computation facilitating easier generation and analysis of frequency responses, reducing reliance on manual graphical methods.[9] Its adoption grew in industries like aerospace for flight control systems and process control for chemical engineering, where it proved valuable for assessing stability margins in complex loops.[7] A key milestone came in the late 20th century with its integration into software tools such as MATLAB and Simulink, enabling automated plotting and interactive design within the Control System Toolbox.[10]Construction and Mathematical Foundation
Construction Process
To construct a Nichols plot, begin by obtaining the open-loop transfer function of the system, denoted as G(s).[11] The first step involves deriving the frequency response by substituting s = j\omega into G(s), resulting in G(j\omega), where \omega is the angular frequency in radians per second.[1] For a suitable range of frequencies, such as from 0.1 to 10 rad/s, compute the magnitude |G(j\omega)| and express it in decibels as $20 \log_{10} |G(j\omega)|.[11] Simultaneously, calculate the phase angle \angle G(j\omega) in degrees for the same frequency values.[11] Plot these computed values as points, with the phase angle on the horizontal x-axis (typically decreasing from 0° to -360°) and the magnitude in decibels on the vertical y-axis (often spanning -60 dB to +40 dB).[1] To form the complete locus curve, evaluate and connect points across multiple frequencies, commonly using logarithmic spacing to capture the system's behavior across decades efficiently.[1] In practice, software tools facilitate this process; for instance, MATLAB'snichols function automates the computation and plotting of the Nichols response from a transfer function model.[12]