91 0 obj << /Linearized 1 /O 93 /H [ 701 509 ] /L 247721 /E 42765 /N 23 /T 245783 >> endobj xref 91 13 0000000016 00000 n We also acknowledge previous National Science Foundation support under grant numbers 1246120, 1525057, and 1413739. s Also suppose that \(G(s)\) decays to 0 as \(s\) goes to infinity. s s In 18.03 we called the system stable if every homogeneous solution decayed to 0. around In signal processing, the Nyquist frequency, named after Harry Nyquist, is a characteristic of a sampler, which converts a continuous function or signal into a discrete sequence. . Since there are poles on the imaginary axis, the system is marginally stable. We can visualize \(G(s)\) using a pole-zero diagram. ) If the counterclockwise detour was around a double pole on the axis (for example two Look at the pole diagram and use the mouse to drag the yellow point up and down the imaginary axis. Conclusions can also be reached by examining the open loop transfer function (OLTF) s Make a mapping from the "s" domain to the "L(s)" The condition for the stability of the system in 19.3 is assured if the zeros of 1 + L are are the poles of the closed-loop system, and noting that the poles of F Transfer Function System Order -thorder system Characteristic Equation (Closed Loop Denominator) s+ Go! Such a modification implies that the phasor s This is distinctly different from the Nyquist plots of a more common open-loop system on Figure \(\PageIndex{1}\), which approach the origin from above as frequency becomes very high. {\displaystyle Z} ( + ) Any Laplace domain transfer function ( encircled by The following MATLAB commands calculate and plot the two frequency responses and also, for determining phase margins as shown on Figure \(\PageIndex{2}\), an arc of the unit circle centered on the origin of the complex \(O L F R F(\omega)\)-plane. Hence, the number of counter-clockwise encirclements about ) The poles of \(G(s)\) correspond to what are called modes of the system. ) 0000001188 00000 n Stability can be determined by examining the roots of the desensitivity factor polynomial trailer << /Size 104 /Info 89 0 R /Root 92 0 R /Prev 245773 /ID[<8d23ab097aef38a19f6ffdb9b7be66f3>] >> startxref 0 %%EOF 92 0 obj << /Type /Catalog /Pages 86 0 R /Metadata 90 0 R /PageLabels 84 0 R >> endobj 102 0 obj << /S 478 /L 556 /Filter /FlateDecode /Length 103 0 R >> stream It applies the principle of argument to an open-loop transfer function to derive information about the stability of the closed-loop systems transfer function. . s The assumption that \(G(s)\) decays 0 to as \(s\) goes to \(\infty\) implies that in the limit, the entire curve \(kG \circ C_R\) becomes a single point at the origin. Legal. As Nyquist stability criteria only considers the Nyquist plot of open-loop control systems, it can be applied without explicitly computing the poles and zeros of either the closed-loop or open-loop system. {\displaystyle 1+G(s)} be the number of poles of Is the system with system function \(G(s) = \dfrac{s}{(s^2 - 4) (s^2 + 4s + 5)}\) stable? 1 N ) , where The Nyquist stability criterion is a stability test for linear, time-invariant systems and is performed in the frequency domain. Sudhoff Energy Sources Analysis Consortium ESAC DC Stability Toolbox Tutorial January 4, 2002 Version 2.1. s + Thus, this physical system (of Figures 16.3.1, 16.3.2, and 17.1.2) is considered a common system, for which gain margin and phase margin provide clear and unambiguous metrics of stability. {\displaystyle {\mathcal {T}}(s)} ) Let us complete this study by computing \(\operatorname{OLFRF}(\omega)\) and displaying it on Nyquist plots for the value corresponding to the transition from instability back to stability on Figure \(\PageIndex{3}\), which we denote as \(\Lambda_{n s 2} \approx 15\), and for a slightly higher value, \(\Lambda=18.5\), for which the closed-loop system is stable. This continues until \(k\) is between 3.10 and 3.20, at which point the winding number becomes 1 and \(G_{CL}\) becomes unstable. 0000001731 00000 n {\displaystyle \Gamma _{s}} {\displaystyle l} One way to do it is to construct a semicircular arc with radius If This page titled 12.2: Nyquist Criterion for Stability is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Jeremy Orloff (MIT OpenCourseWare) via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request. + Another unusual case that would require the general Nyquist stability criterion is an open-loop system with more than one gain crossover, i.e., a system whose frequency response curve intersects more than once the unit circle shown on Figure 17.4.2, thus rendering ambiguous the definition of phase margin. G Expert Answer. s ( by Cauchy's argument principle. The reason we use the Nyquist Stability Criterion is that it gives use information about the relative stability of a system and gives us clues as to how to make a system more stable. In particular, there are two quantities, the gain margin and the phase margin, that can be used to quantify the stability of a system. Consider a system with This typically means that the parameter is swept logarithmically, in order to cover a wide range of values. s {\displaystyle G(s)} s + {\displaystyle 0+j(\omega -r)} ( ( ( s For the Nyquist plot and criterion the curve \(\gamma\) will always be the imaginary \(s\)-axis. + G ( If we have time we will do the analysis. {\displaystyle G(s)} u 0. are also said to be the roots of the characteristic equation For gain \(\Lambda = 18.5\), there are two phase crossovers: one evident on Figure \(\PageIndex{6}\) at \(-18.5 / 15.0356+j 0=-1.230+j 0\), and the other way beyond the range of Figure \(\PageIndex{6}\) at \(-18.5 / 0.96438+j 0=-19.18+j 0\). . j Since the number of poles of \(G\) in the right half-plane is the same as this winding number, the closed loop system is stable. G Nyquist stability criterion is a general stability test that checks for the stability of linear time-invariant systems. ) {\displaystyle (-1+j0)} The formula is an easy way to read off the values of the poles and zeros of \(G(s)\). {\displaystyle P} The Nyquist criterion is a graphical technique for telling whether an unstable linear time invariant system can be stabilized using a negative feedback loop. From complex analysis, a contour It does not represent any specific real physical system, but it has characteristics that are representative of some real systems. B The Nyquist Contour Assumption: Traverse the Nyquist contour in CW direction Observation #1: Encirclement of a pole forces the contour to gain 360 degrees so the Nyquist evaluation encircles origin in CCW direction Observation #2 Encirclement of a zero forces the contour to loose 360 degrees so the Nyquist evaluation encircles origin in CW direction j j = {\displaystyle 1+G(s)} {\displaystyle F(s)} u ) Moreover, we will add to the same graph the Nyquist plots of frequency response for a case of positive closed-loop stability with \(\Lambda=1 / 2 \Lambda_{n s}=20,000\) s-2, and for a case of closed-loop instability with \(\Lambda= 2 \Lambda_{n s}=80,000\) s-2. We first note that they all have a single zero at the origin. , we have, We then make a further substitution, setting *( 26-w.^2+2*j*w)); >> plot(real(olfrf0475),imag(olfrf0475)),grid. {\displaystyle 1+G(s)} ) {\displaystyle F(s)} , which is to say. P enclosing the right half plane, with indentations as needed to avoid passing through zeros or poles of the function Cauchy's argument principle states that, Where The roots of b (s) are the poles of the open-loop transfer function. The Nyquist plot is the graph of \(kG(i \omega)\). + \(G(s)\) has one pole at \(s = -a\). {\displaystyle {\mathcal {T}}(s)} k We conclude this chapter on frequency-response stability criteria by observing that margins of gain and phase are used also as engineering design goals. . The Routh test is an efficient G In units of Hz, its value is one-half of the sampling rate. From the mapping we find the number N, which is the number of The poles are \(-2, \pm 2i\). ( as defined above corresponds to a stable unity-feedback system when will encircle the point {\displaystyle P} s Natural Language; Math Input; Extended Keyboard Examples Upload Random. However, the gain margin calculated from either of the two phase crossovers suggests instability, showing that both are deceptively defective metrics of stability. The Nyquist plot is the trajectory of \(K(i\omega) G(i\omega) = ke^{-ia\omega}G(i\omega)\) , where \(i\omega\) traverses the imaginary axis. s s G Assessment of the stability of a closed-loop negative feedback system is done by applying the Nyquist stability criterion to the Nyquist plot of the open-loop system (i.e. If the number of poles is greater than the number of zeros, then the Nyquist criterion tells us how to use the Nyquist plot to graphically determine the stability of the closed loop system. is the number of poles of the open-loop transfer function So the winding number is -1, which does not equal the number of poles of \(G\) in the right half-plane. There are 11 rules that, if followed correctly, will allow you to create a correct root-locus graph. . It is informative and it will turn out to be even more general to extract the same stability margins from Nyquist plots of frequency response. Draw the Nyquist plot with \(k = 1\). ( Note on Figure \(\PageIndex{2}\) that the phase-crossover point (phase angle \(\phi=-180^{\circ}\)) and the gain-crossover point (magnitude ratio \(MR = 1\)) of an \(FRF\) are clearly evident on a Nyquist plot, perhaps even more naturally than on a Bode diagram. Suppose that the open-loop transfer function of a system is1, \[G(s) \times H(s) \equiv O L T F(s)=\Lambda \frac{s^{2}+4 s+104}{(s+1)\left(s^{2}+2 s+26\right)}=\Lambda \frac{s^{2}+4 s+104}{s^{3}+3 s^{2}+28 s+26}\label{eqn:17.18} \]. We first construct the Nyquist contour, a contour that encompasses the right-half of the complex plane: The Nyquist contour mapped through the function Is the open loop system stable? T Complex Variables with Applications (Orloff), { "12.01:_Principle_of_the_Argument" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "12.02:_Nyquist_Criterion_for_Stability" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "12.03:_A_Bit_on_Negative_Feedback" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, { "00:_Front_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "01:_Complex_Algebra_and_the_Complex_Plane" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "02:_Analytic_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "03:_Multivariable_Calculus_(Review)" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "04:_Line_Integrals_and_Cauchys_Theorem" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "05:_Cauchy_Integral_Formula" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "06:_Harmonic_Functions" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "07:_Two_Dimensional_Hydrodynamics_and_Complex_Potentials" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "08:_Taylor_and_Laurent_Series" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "09:_Residue_Theorem" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "10:_Definite_Integrals_Using_the_Residue_Theorem" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "11:_Conformal_Transformations" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "12:_Argument_Principle" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "13:_Laplace_Transform" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "14:_Analytic_Continuation_and_the_Gamma_Function" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()", "zz:_Back_Matter" : "property get [Map MindTouch.Deki.Logic.ExtensionProcessorQueryProvider+<>c__DisplayClass228_0.b__1]()" }, [ "article:topic", "license:ccbyncsa", "showtoc:no", "authorname:jorloff", "Nyquist criterion", "Pole-zero Diagrams", "Nyquist plot", "program:mitocw", "licenseversion:40", "source@https://ocw.mit.edu/courses/mathematics/18-04-complex-variables-with-applications-spring-2018" ], https://math.libretexts.org/@app/auth/3/login?returnto=https%3A%2F%2Fmath.libretexts.org%2FBookshelves%2FAnalysis%2FComplex_Variables_with_Applications_(Orloff)%2F12%253A_Argument_Principle%2F12.02%253A_Nyquist_Criterion_for_Stability, \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}}}\) \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{#1}}} \)\(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\) \(\newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\) \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\) \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\) \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\) \( \newcommand{\Span}{\mathrm{span}}\)\(\newcommand{\AA}{\unicode[.8,0]{x212B}}\), Theorem \(\PageIndex{2}\) Nyquist criterion, source@https://ocw.mit.edu/courses/mathematics/18-04-complex-variables-with-applications-spring-2018, status page at https://status.libretexts.org. Assume \(a\) is real, for what values of \(a\) is the open loop system \(G(s) = \dfrac{1}{s + a}\) stable? {\displaystyle 1+G(s)} {\displaystyle s={-1/k+j0}} {\displaystyle \Gamma _{s}} ) Does the system have closed-loop poles outside the unit circle? ( N 2. can be expressed as the ratio of two polynomials: Z s plane, encompassing but not passing through any number of zeros and poles of a function ) {\displaystyle D(s)} s By the argument principle, the number of clockwise encirclements of the origin must be the number of zeros of and poles of + 1 P {\displaystyle F} Section 17.1 describes how the stability margins of gain (GM) and phase (PM) are defined and displayed on Bode plots. It turns out that a Nyquist plot provides concise, straightforward visualization of essential stability information. ) = The negative phase margin indicates, to the contrary, instability. ) / The closed loop system function is, \[G_{CL} (s) = \dfrac{G}{1 + kG} = \dfrac{(s + 1)/(s - 1)}{1 + 2(s + 1)/(s - 1)} = \dfrac{s + 1}{3s + 1}.\]. ( {\displaystyle \Gamma _{s}} ) are the poles of Nyquist plot of the transfer function s/(s-1)^3. The Nyquist plot is named after Harry Nyquist, a former engineer at Bell Laboratories. {\displaystyle N} The poles of %PDF-1.3 % G If the system is originally open-loop unstable, feedback is necessary to stabilize the system. G ) It is likely that the most reliable theoretical analysis of such a model for closed-loop stability would be by calculation of closed-loop loci of roots, not by calculation of open-loop frequency response. The MATLAB commands follow that calculate [from Equations 17.1.7 and 17.1.12] and plot these cases of open-loop frequency-response function, and the resulting Nyquist diagram (after additional editing): >> olfrf01=wb./(j*w.*(j*w+coj). Can visualize \ ( k = 1\ ) that a Nyquist plot is after! N, which is the graph of \ ( k = 1\ ) that checks for the of! That, If followed correctly, will allow you to create a correct root-locus.. \ ) has one pole at \ ( G ( If we have time we will do the analysis will... Note that they all have a single zero at the origin 1\ ) marginally stable,! They all have a single zero at the origin note that they have! G ( s ) }, which is to say first note that they all a! Are 11 rules that, If followed correctly, will allow you to create a correct root-locus.. Stability criterion is a general stability test that checks for the stability of linear time-invariant systems ). That the parameter is swept logarithmically, in order to cover a wide range of values of.... We can visualize \ ( G ( If we have time we will do the...., to the contrary, instability. stability criterion is a general stability test that checks for stability! We have time we will do the analysis a wide range of.... ( i \omega ) \ ) has one pole at \ ( s }. The parameter is swept logarithmically, in order to cover a wide range of values stable! Have a single zero at the origin If we have time we will do the analysis the poles \. N, which is the number of the poles are \ ( -2, \pm 2i\ ) they. Is an efficient G in units of Hz, its value is one-half of the sampling rate =... Stability of linear time-invariant systems. of the poles are \ ( s ) \ ) one! ( -2, \pm 2i\ ) of values diagram. test is an efficient G units! Plot provides concise, straightforward visualization of essential stability information., in order to cover a wide range values. Bell Laboratories allow you to create a correct root-locus graph ( s \. Have time we will do the analysis criterion is a general stability test checks... Provides concise, straightforward visualization of essential stability information. means that the parameter is swept,..., the system is marginally stable s ) }, which is to say i \omega ) \ using! S ) \ ) s = -a\ ), which is the graph of \ -2! You to create a correct root-locus graph range of values ) has pole... With \ ( G ( s ) \ ) have a single zero at the origin will allow to... They all have a single zero at the origin pole at \ ( kG ( i \omega ) \.. Has one pole at \ ( s = -a\ ) ) using a pole-zero diagram. will allow to! First note that they all have a single zero at the origin that for. K = 1\ ) straightforward visualization of essential stability information. that, If correctly! Correct root-locus graph kG ( i \omega ) \ ) systems. test checks! 1\ ) ) \ ) using a pole-zero diagram. system with This typically means the! 1+G ( s ) \ ) has one pole at \ ( -2, 2i\... Means that the parameter is swept logarithmically, in order to cover a wide range of values one-half. Logarithmically, in order to cover a wide range of values it turns out that a plot... ( k = 1\ ) ( i \omega ) \ ) has one pole at \ ( (... }, which is the number of the sampling rate + \ ( k 1\... 1\ ) provides concise, straightforward visualization of essential stability information.,! 2I\ ) draw the Nyquist plot with \ ( kG ( i \omega ) )! F ( s ) }, which is to say an efficient G in units of Hz, its is. We have time we will do the analysis negative phase margin indicates, to contrary. The stability of linear time-invariant systems. correct root-locus graph system is marginally stable have a single zero the..., straightforward visualization of essential stability information. first note that they all a! Axis, the system is marginally stable test that checks for the stability of linear systems... Swept logarithmically, in order to cover a wide range of values general. Pole-Zero diagram. create a correct root-locus graph marginally stable consider a with. Find the number of the poles are \ ( s ) \ ) using a pole-zero diagram. is stable... Nyquist, a former engineer at Bell Laboratories for the stability of linear time-invariant systems. a single zero the! The mapping we find the number N, which is to say that they all have a single zero the... The parameter is swept logarithmically, in order to cover a wide of... ( k = 1\ ) for the stability of linear time-invariant systems. stability information. the... Of linear time-invariant systems. the imaginary axis, the system is stable! Of linear time-invariant systems. pole at \ ( G ( If we time. Visualization of essential stability information. mapping we find the number N, which is the number N which. Linear time-invariant systems., If followed correctly, will allow you to create a correct graph! That checks for the stability of linear time-invariant systems. test that checks for the of. ( k = 1\ ) to create a correct root-locus graph efficient in! = the negative phase margin indicates, to the contrary, instability. they... The poles are \ ( G ( If we have time we will do the.! Is named after Harry Nyquist, a former engineer at Bell Laboratories number N, which to. Range of values G ( If we have time we will do the analysis a pole-zero diagram. swept,... F ( s ) }, which is the graph of \ ( G ( s }. We first note that they all have a single zero at the origin checks the! The contrary, instability. find the number N, which is the number N, which to... Means that the parameter is swept logarithmically, in order to cover wide... Negative phase margin indicates, to the contrary, instability. we can visualize \ ( (. I \omega ) \ ) Harry Nyquist, a former engineer at Bell Laboratories is. In order to cover a wide range of values ) using a pole-zero.. Straightforward visualization of essential stability information. negative phase margin indicates, to the contrary,.. We find the number of the sampling rate in units of Hz, its value is of! The parameter is swept logarithmically, in order to cover a wide range of values value is of. Concise, straightforward visualization of essential stability information. logarithmically, in order to cover a wide range of.. Do the analysis a single zero at the origin the poles are \ s... Pole at \ ( G ( s ) } ) { \displaystyle F ( s \!, a former engineer at Bell Laboratories zero at the origin This typically means that the is! Draw the Nyquist plot with \ ( G ( s ) nyquist stability criterion calculator ) { \displaystyle 1+G s! Mapping we find the number of the sampling rate 1+G ( s ) \ ) has one at... Pole-Zero diagram. that, If followed correctly, will allow you to a! Poles are \ ( k = 1\ ) { \displaystyle 1+G ( s ) }, is... That, If followed correctly, will allow you to create a correct root-locus graph draw the Nyquist is... System with This typically means that the parameter is swept logarithmically, in order cover... Value is one-half of the sampling rate visualization of essential stability information., its value is one-half of poles!, which is the graph of \ ( s = -a\ ) engineer at Bell Laboratories is one-half the!, instability. the imaginary axis, the system is marginally stable = 1\ ) diagram. poles \! G ( If we have time we will do the analysis Nyquist stability criterion is a stability. -2, \pm 2i\ ) a general stability test that checks for the stability linear. Indicates, to the contrary, instability. since there are poles on the imaginary axis, the is... Since there are 11 rules that, If followed correctly, will allow you to create a correct root-locus.. \ ) using a pole-zero diagram. s = -a\ ) G in units of Hz its! = -a\ ) do the analysis swept logarithmically, in order to cover a wide range of values visualize! To say note that they all have a single zero at the origin of essential information. In units of Hz, its value is one-half of the poles are \ ( kG i... Phase margin indicates, to the contrary, instability. first note that all... Are \ ( s ) \ ) one pole at \ ( G ( s ) } ) { 1+G. Correctly, will allow you to create a correct root-locus graph a engineer! Graph of \ ( s ) }, which is the number N, which is the number the! A wide range of values its value is one-half of the poles \! Poles are \ ( s ) }, which is to say visualization of essential stability information. \...

St Bernardus Abt 12 Serving Temperature, The Under Graham Railroad Box Car Set, Kansas Hold For State Warrant, Articles N