Bode plot
In electrical engineering and control theory, a Bode plot /ˈboʊdi/ is a graph of the frequency response of a system. It is usually a combination of a Bode magnitude plot, expressing the magnitude (usually in decibels) of the frequency response, and a Bode phase plot, expressing the phase shift.
As originally conceived by Hendrik Wade Bode in the 1930s, the plot is an asymptotic approximation of the frequency response, using straight line segments.^{[1]}
Overview
Among his several important contributions to circuit theory and control theory, engineer Hendrik Wade Bode, while working at Bell Labs in the 1930s, devised a simple but accurate method for graphing gain and phaseshift plots. These bear his name, Bode gain plot and Bode phase plot. "Bode" is often pronounced /ˈboʊdi/ BOHdee, although the Dutch pronunciation is Boduh. (Dutch: [ˈboːdə]).^{[2]}^{[3]}
Bode was faced with the problem of designing stable amplifiers with feedback for use in telephone networks. He developed the graphical design technique of the Bode plots to show the gain margin and phase margin required to maintain stability under variations in circuit characteristics caused during manufacture or during operation.^{[4]} The principles developed were applied to design problems of servomechanisms and other feedback control systems. The Bode plot is an example of analysis in the frequency domain.
Definition
The Bode plot for a linear, timeinvariant system with transfer function [math]\displaystyle{ H(s) }[/math] ([math]\displaystyle{ s }[/math] being the complex frequency in the Laplace domain) consists of a magnitude plot and a phase plot.
The Bode magnitude plot is the graph of the function [math]\displaystyle{ H(s=j \omega) }[/math] of frequency [math]\displaystyle{ \omega }[/math] (with [math]\displaystyle{ j }[/math] being the imaginary unit). The [math]\displaystyle{ \omega }[/math]axis of the magnitude plot is logarithmic and the magnitude is given in decibels, i.e., a value for the magnitude [math]\displaystyle{ H }[/math] is plotted on the axis at [math]\displaystyle{ 20 \log_{10} H }[/math].
The Bode phase plot is the graph of the phase, commonly expressed in degrees, of the transfer function [math]\displaystyle{ \arg \left( H(s =j \omega) \right) }[/math] as a function of [math]\displaystyle{ \omega }[/math]. The phase is plotted on the same logarithmic [math]\displaystyle{ \omega }[/math]axis as the magnitude plot, but the value for the phase is plotted on a linear vertical axis.
Frequency response
This section illustrates that a Bode plot is a visualization of the frequency response of a system.
Consider a linear, timeinvariant system with transfer function [math]\displaystyle{ H(s) }[/math]. Assume that the system is subject to a sinusoidal input with frequency [math]\displaystyle{ \omega }[/math],
 [math]\displaystyle{ u(t) = \sin (\omega t), }[/math]
that is applied persistently, i.e. from a time [math]\displaystyle{ \infty }[/math] to a time [math]\displaystyle{ t }[/math]. The response will be of the form
 [math]\displaystyle{ y(t) = y_0 \sin (\omega t + \varphi), }[/math]
i.e., also a sinusoidal signal with amplitude [math]\displaystyle{ y_0 }[/math] shifted by a phase [math]\displaystyle{ \varphi }[/math] with respect to the input.
It can be shown^{[5]} that the magnitude of the response is

[math]\displaystyle{ y_0 = H(\mathrm{j} \omega) }[/math]
(
)
and that the phase shift is

[math]\displaystyle{ \varphi = \arg H(\mathrm{j} \omega). }[/math]
(
)
In summary, subjected to an input with frequency [math]\displaystyle{ \omega }[/math], the system responds at the same frequency with an output that is amplified by a factor [math]\displaystyle{ H(\mathrm{j} \omega) }[/math] and phaseshifted by [math]\displaystyle{ \arg H(\mathrm{j} \omega) }[/math]. These quantities, thus, characterize the frequency response and are shown in the Bode plot.
Rules for handmade Bode plot
For many practical problems, the detailed Bode plots can be approximated with straightline segments that are asymptotes of the precise response. The effect of each of the terms of a multiple element transfer function can be approximated by a set of straight lines on a Bode plot. This allows a graphical solution of the overall frequency response function. Before widespread availability of digital computers, graphical methods were extensively used to reduce the need for tedious calculation; a graphical solution could be used to identify feasible ranges of parameters for a new design.
The premise of a Bode plot is that one can consider the log of a function in the form
 [math]\displaystyle{ f(x) = A \prod (x  c_n)^{a_n} }[/math]
as a sum of the logs of its zeros and poles:
 [math]\displaystyle{ \log(f(x)) = \log(A) + \sum a_n \log(x  c_n). }[/math]
This idea is used explicitly in the method for drawing phase diagrams. The method for drawing amplitude plots implicitly uses this idea, but since the log of the amplitude of each pole or zero always starts at zero and only has one asymptote change (the straight lines), the method can be simplified.
Straightline amplitude plot
Amplitude decibels is usually done using [math]\displaystyle{ \text{dB} = 20 \log_{10}(X) }[/math] to define decibels. Given a transfer function in the form
 [math]\displaystyle{ H(s) = A \prod \frac{(s  x_n)^{a_n}}{(s  y_n)^{b_n}}, }[/math]
where [math]\displaystyle{ x_n }[/math] and [math]\displaystyle{ y_n }[/math] are constants, [math]\displaystyle{ s = \mathrm{j}\omega }[/math], [math]\displaystyle{ a_n, b_n \gt 0 }[/math], and [math]\displaystyle{ H }[/math] is the transfer function:
 At every value of s where [math]\displaystyle{ \omega = x_n }[/math] (a zero), increase the slope of the line by [math]\displaystyle{ 20 a_n\ \text{dB} }[/math] per decade.
 At every value of s where [math]\displaystyle{ \omega = y_n }[/math] (a pole), decrease the slope of the line by [math]\displaystyle{ 20 b_n\ \text{dB} }[/math] per decade.
 The initial value of the graph depends on the boundaries. The initial point is found by putting the initial angular frequency [math]\displaystyle{ \omega }[/math] into the function and finding [math]\displaystyle{ H(\mathrm{j}\omega) }[/math].
 The initial slope of the function at the initial value depends on the number and order of zeros and poles that are at values below the initial value, and is found using the first two rules.
To handle irreducible 2ndorder polynomials, [math]\displaystyle{ ax^2 + bx + c }[/math] can, in many cases, be approximated as [math]\displaystyle{ (\sqrt{a}x + \sqrt{c})^2 }[/math].
Note that zeros and poles happen when [math]\displaystyle{ \omega }[/math] is equal to a certain [math]\displaystyle{ x_n }[/math] or [math]\displaystyle{ y_n }[/math]. This is because the function in question is the magnitude of [math]\displaystyle{ H(\mathrm{j}\omega) }[/math], and since it is a complex function, [math]\displaystyle{ H(\mathrm{j}\omega) = \sqrt{H \cdot H^*} }[/math]. Thus at any place where there is a zero or pole involving the term [math]\displaystyle{ (s + x_n) }[/math], the magnitude of that term is [math]\displaystyle{ \sqrt{(x_n + \mathrm{j}\omega)(x_n  \mathrm{j}\omega)} = \sqrt{x_n^2 + \omega^2} }[/math].
Corrected amplitude plot
To correct a straightline amplitude plot:
 At every zero, put a point [math]\displaystyle{ 3 a_n\ \text{dB} }[/math] above the line.
 At every pole, put a point [math]\displaystyle{ 3 b_n\ \text{dB} }[/math] below the line.
 Draw a smooth curve through those points using the straight lines as asymptotes (lines which the curve approaches).
Note that this correction method does not incorporate how to handle complex values of [math]\displaystyle{ x_n }[/math] or [math]\displaystyle{ y_n }[/math]. In the case of an irreducible polynomial, the best way to correct the plot is to actually calculate the magnitude of the transfer function at the pole or zero corresponding to the irreducible polynomial, and put that dot over or under the line at that pole or zero.
Straightline phase plot
Given a transfer function in the same form as above,
 [math]\displaystyle{ H(s) = A \prod \frac{(s  x_n)^{a_n}}{(s  y_n)^{b_n}}, }[/math]
the idea is to draw separate plots for each pole and zero, then add them up. The actual phase curve is given by
 [math]\displaystyle{ \varphi(s) = \arctan \frac{\operatorname{Im}[H(s)]}{\operatorname{Re}[H(s)]}. }[/math]
To draw the phase plot, for each pole and zero:
 If [math]\displaystyle{ A }[/math] is positive, start line (with zero slope) at 0°.
 If [math]\displaystyle{ A }[/math] is negative, start line (with zero slope) at −180°.
 If the sum of the number of unstable zeros and poles is odd, add 180° to that basis.
 At every [math]\displaystyle{ \omega = x_n }[/math] (for stable zeros [math]\displaystyle{ \operatorname{Re}(z) \lt 0 }[/math]), increase the slope by [math]\displaystyle{ 45 a_n }[/math] degrees per decade, beginning one decade before [math]\displaystyle{ \omega = x_n }[/math] (e.g., [math]\displaystyle{ x_n/10 }[/math]).
 At every [math]\displaystyle{ \omega = y_n }[/math] (for stable poles [math]\displaystyle{ \operatorname{Re}(p) \lt 0 }[/math]), decrease the slope by [math]\displaystyle{ 45 b_n }[/math] degrees per decade, beginning one decade before [math]\displaystyle{ \omega = y_n }[/math] (e.g., [math]\displaystyle{ y_n/10 }[/math]).
 "Unstable" (right halfplane) poles and zeros ([math]\displaystyle{ \operatorname{Re}(s) \gt 0 }[/math]) have opposite behavior.
 Flatten the slope again when the phase has changed by [math]\displaystyle{ 90 a_n }[/math] degrees (for a zero) or [math]\displaystyle{ 90 b_n }[/math] degrees (for a pole).
 After plotting one line for each pole or zero, add the lines together to obtain the final phase plot; that is, the final phase plot is the superposition of each earlier phase plot.
Example
To create a straightline plot for a firstorder (onepole) lowpass filter, one considers the normalized form of the transfer function in terms of the angular frequency:
 [math]\displaystyle{ H_{\text{lp}}(\mathrm{j} \omega) = \frac{1}{1 + \mathrm{j} \frac{\omega}{\omega_\text{c}}}. }[/math]
The Bode plot is shown in Figure 1(b) above, and construction of the straightline approximation is discussed next.
Magnitude plot
The magnitude (in decibels) of the transfer function above (normalized and converted to angularfrequency form), given by the decibel gain expression [math]\displaystyle{ A_\text{vdB} }[/math]:
 [math]\displaystyle{ \begin{align} A_\text{vdB} &= 20 \logH_{\text{lp}}(\mathrm{j}\omega) \\ &= 20 \log \frac{1}{\left 1 + \mathrm{j} \frac{\omega}{\omega_\text{c}} \right} \\ &= 20 \log \left 1 + \mathrm{j} \frac{\omega}{\omega_\text{c}} \right \\ &= 10 \log \left( 1 + \frac{\omega^2}{\omega_\text{c}^2} \right). \end{align} }[/math]
Then plotted versus input frequency [math]\displaystyle{ \omega }[/math] on a logarithmic scale, can be approximated by two lines, forming the asymptotic (approximate) magnitude Bode plot of the transfer function:
 The first line for angular frequencies below [math]\displaystyle{ \omega_\text{c} }[/math] is a horizontal line at 0 dB, since at low frequencies the [math]\displaystyle{ \omega/\omega_\text{c} }[/math] term is small and can be neglected, making the decibel gain equation above equal to zero.
 The second line for angular frequencies above [math]\displaystyle{ \omega_\text{c} }[/math] is a line with a slope of −20 dB per decade, since at high frequencies the [math]\displaystyle{ \omega/\omega_\text{c} }[/math] term dominates, and the decibel gain expression above simplifies to [math]\displaystyle{ 20 \log(\omega/\omega_\text{c}) }[/math], which is a straight line with a slope of −20 dB per decade.
These two lines meet at the corner frequency [math]\displaystyle{ \omega_\text{c} }[/math]. From the plot, it can be seen that for frequencies well below the corner frequency, the circuit has an attenuation of 0 dB, corresponding to a unity passband gain, i.e. the amplitude of the filter output equals the amplitude of the input. Frequencies above the corner frequency are attenuated – the higher the frequency, the higher the attenuation.
Phase plot
The phase Bode plot is obtained by plotting the phase angle of the transfer function given by
 [math]\displaystyle{ \arg H_{\text{lp}}(\mathrm{j} \omega) = \tan^{1}\frac{\omega}{\omega_\text{c}} }[/math]
versus [math]\displaystyle{ \omega }[/math], where [math]\displaystyle{ \omega }[/math] and [math]\displaystyle{ \omega_\text{c} }[/math] are the input and cutoff angular frequencies respectively. For input frequencies much lower than corner, the ratio [math]\displaystyle{ \omega/\omega_\text{c} }[/math] is small, and therefore the phase angle is close to zero. As the ratio increases, the absolute value of the phase increases and becomes −45° when [math]\displaystyle{ \omega = \omega_\text{c} }[/math]. As the ratio increases for input frequencies much greater than the corner frequency, the phase angle asymptotically approaches −90°. The frequency scale for the phase plot is logarithmic.
Normalized plot
The horizontal frequency axis, in both the magnitude and phase plots, can be replaced by the normalized (nondimensional) frequency ratio [math]\displaystyle{ \omega/\omega_\text{c} }[/math]. In such a case the plot is said to be normalized, and units of the frequencies are no longer used, since all input frequencies are now expressed as multiples of the cutoff frequency [math]\displaystyle{ \omega_\text{c} }[/math].
An example with zero and pole
Figures 25 further illustrate construction of Bode plots. This example with both a pole and a zero shows how to use superposition. To begin, the components are presented separately.
Figure 2 shows the Bode magnitude plot for a zero and a lowpass pole, and compares the two with the Bode straight line plots. The straightline plots are horizontal up to the pole (zero) location and then drop (rise) at 20 dB/decade. The second Figure 3 does the same for the phase. The phase plots are horizontal up to a frequency factor of ten below the pole (zero) location and then drop (rise) at 45°/decade until the frequency is ten times higher than the pole (zero) location. The plots then are again horizontal at higher frequencies at a final, total phase change of 90°.
Figure 4 and Figure 5 show how superposition (simple addition) of a pole and zero plot is done. The Bode straight line plots again are compared with the exact plots. The zero has been moved to higher frequency than the pole to make a more interesting example. Notice in Figure 4 that the 20 dB/decade drop of the pole is arrested by the 20 dB/decade rise of the zero resulting in a horizontal magnitude plot for frequencies above the zero location. Notice in Figure 5 in the phase plot that the straightline approximation is pretty approximate in the region where both pole and zero affect the phase. Notice also in Figure 5 that the range of frequencies where the phase changes in the straight line plot is limited to frequencies a factor of ten above and below the pole (zero) location. Where the phase of the pole and the zero both are present, the straightline phase plot is horizontal because the 45°/decade drop of the pole is arrested by the overlapping 45°/decade rise of the zero in the limited range of frequencies where both are active contributors to the phase.
Gain margin and phase margin
Bode plots are used to assess the stability of negativefeedback amplifiers by finding the gain and phase margins of an amplifier. The notion of gain and phase margin is based upon the gain expression for a negative feedback amplifier given by
 [math]\displaystyle{ A_\text{FB} = \frac{A_\text{OL}}{1 + \beta A_\text{OL}}, }[/math]
where A_{FB} is the gain of the amplifier with feedback (the closedloop gain), β is the feedback factor, and A_{OL} is the gain without feedback (the openloop gain). The gain A_{OL} is a complex function of frequency, with both magnitude and phase.^{[note 1]} Examination of this relation shows the possibility of infinite gain (interpreted as instability) if the product βA_{OL} = −1 (that is, the magnitude of βA_{OL} is unity and its phase is −180°, the socalled Barkhausen stability criterion). Bode plots are used to determine just how close an amplifier comes to satisfying this condition.
Key to this determination are two frequencies. The first, labeled here as f_{180}, is the frequency where the openloop gain flips sign. The second, labeled here f_{0 dB}, is the frequency where the magnitude of the product βA_{OL} = 1 = 0 dB. That is, frequency f_{180} is determined by the condition
 [math]\displaystyle{ \beta A_\text{OL}(f_{180}) = \beta A_\text{OL}(f_{180}) = \beta A_\text{OL}_{180}, }[/math]
where vertical bars denote the magnitude of a complex number, and frequency f_{0 dB} is determined by the condition
 [math]\displaystyle{ \beta A_\text{OL}(f_\text{0 dB}) = 1. }[/math]
One measure of proximity to instability is the gain margin. The Bode phase plot locates the frequency where the phase of βA_{OL} reaches −180°, denoted here as frequency f_{180}. Using this frequency, the Bode magnitude plot finds the magnitude of βA_{OL}. If βA_{OL}_{180} ≥ 1, the amplifier is unstable, as mentioned. If βA_{OL}_{180} < 1, instability does not occur, and the separation in dB of the magnitude of βA_{OL}_{180} from βA_{OL} = 1 is called the gain margin. Because a magnitude of 1 is 0 dB, the gain margin is simply one of the equivalent forms: [math]\displaystyle{ 20 \log_{10} \beta A_\text{OL}_{180} = 20 \log_{10} A_\text{OL}  20 \log_{10} \beta^{1} }[/math].
Another equivalent measure of proximity to instability is the phase margin. The Bode magnitude plot locates the frequency where the magnitude of βA_{OL} reaches unity, denoted here as frequency f_{0 dB}. Using this frequency, the Bode phase plot finds the phase of βA_{OL}. If the phase of βA_{OL}(f_{0 dB}) > −180°, the instability condition cannot be met at any frequency (because its magnitude is going to be < 1 when f = f_{180}), and the distance of the phase at f_{0 dB} in degrees above −180° is called the phase margin.
If a simple yes or no on the stability issue is all that is needed, the amplifier is stable if f_{0 dB} < f_{180}. This criterion is sufficient to predict stability only for amplifiers satisfying some restrictions on their pole and zero positions (minimum phase systems). Although these restrictions usually are met, if they are not, then another method must be used, such as the Nyquist plot.^{[6]}^{[7]} Optimal gain and phase margins may be computed using Nevanlinna–Pick interpolation theory.^{[8]}
Examples using Bode plots
Figures 6 and 7 illustrate the gain behavior and terminology. For a threepole amplifier, Figure 6 compares the Bode plot for the gain without feedback (the openloop gain) A_{OL} with the gain with feedback A_{FB} (the closedloop gain). See negative feedback amplifier for more detail.
In this example, A_{OL} = 100 dB at low frequencies, and 1 / β = 58 dB. At low frequencies, A_{FB} ≈ 58 dB as well.
Because the openloop gain A_{OL} is plotted and not the product β A_{OL}, the condition A_{OL} = 1 / β decides f_{0 dB}. The feedback gain at low frequencies and for large A_{OL} is A_{FB} ≈ 1 / β (look at the formula for the feedback gain at the beginning of this section for the case of large gain A_{OL}), so an equivalent way to find f_{0 dB} is to look where the feedback gain intersects the openloop gain. (Frequency f_{0 dB} is needed later to find the phase margin.)
Near this crossover of the two gains at f_{0 dB}, the Barkhausen criteria are almost satisfied in this example, and the feedback amplifier exhibits a massive peak in gain (it would be infinity if β A_{OL} = −1). Beyond the unity gain frequency f_{0 dB}, the openloop gain is sufficiently small that A_{FB} ≈ A_{OL} (examine the formula at the beginning of this section for the case of small A_{OL}).
Figure 7 shows the corresponding phase comparison: the phase of the feedback amplifier is nearly zero out to the frequency f_{180} where the openloop gain has a phase of −180°. In this vicinity, the phase of the feedback amplifier plunges abruptly downward to become almost the same as the phase of the openloop amplifier. (Recall, A_{FB} ≈ A_{OL} for small A_{OL}.)
Comparing the labeled points in Figure 6 and Figure 7, it is seen that the unity gain frequency f_{0 dB} and the phaseflip frequency f_{180} are very nearly equal in this amplifier, f_{180} ≈ f_{0 dB} ≈ 3.332 kHz, which means the gain margin and phase margin are nearly zero. The amplifier is borderline stable.
Figures 8 and 9 illustrate the gain margin and phase margin for a different amount of feedback β. The feedback factor is chosen smaller than in Figure 6 or 7, moving the condition  β A_{OL}  = 1 to lower frequency. In this example, 1 / β = 77 dB, and at low frequencies A_{FB} ≈ 77 dB as well.
Figure 8 shows the gain plot. From Figure 8, the intersection of 1 / β and A_{OL} occurs at f_{0 dB} = 1 kHz. Notice that the peak in the gain A_{FB} near f_{0 dB} is almost gone.^{[note 2]}^{[9]}
Figure 9 is the phase plot. Using the value of f_{0 dB} = 1 kHz found above from the magnitude plot of Figure 8, the openloop phase at f_{0 dB} is −135°, which is a phase margin of 45° above −180°.
Using Figure 9, for a phase of −180° the value of f_{180} = 3.332 kHz (the same result as found earlier, of course^{[note 3]}). The openloop gain from Figure 8 at f_{180} is 58 dB, and 1 / β = 77 dB, so the gain margin is 19 dB.
Stability is not the sole criterion for amplifier response, and in many applications a more stringent demand than stability is good step response. As a rule of thumb, good step response requires a phase margin of at least 45°, and often a margin of over 70° is advocated, particularly where component variation due to manufacturing tolerances is an issue.^{[9]} See also the discussion of phase margin in the step response article.
Bode plotter
The Bode plotter is an electronic instrument resembling an oscilloscope, which produces a Bode diagram, or a graph, of a circuit's voltage gain or phase shift plotted against frequency in a feedback control system or a filter. An example of this is shown in Figure 10. It is extremely useful for analyzing and testing filters and the stability of feedback control systems, through the measurement of corner (cutoff) frequencies and gain and phase margins.
This is identical to the function performed by a vector network analyzer, but the network analyzer is typically used at much higher frequencies.
For education and research purposes, plotting Bode diagrams for given transfer functions facilitates better understanding and getting faster results (see external links).
Related plots
Two related plots that display the same data in different coordinate systems are the Nyquist plot and the Nichols plot. These are parametric plots, with frequency as the input and magnitude and phase of the frequency response as the output. The Nyquist plot displays these in polar coordinates, with magnitude mapping to radius and phase to argument (angle). The Nichols plot displays these in rectangular coordinates, on the log scale.
Figure 11: A Nyquist plot.
Figure 12: A Nichols plot of the same response from Figure 11.
See also
 Analog signal processing
 Phase margin
 Bode's sensitivity integral
 Bode's magnitude (gain)–phase relation
 Electrochemical impedance spectroscopy
Notes
 ↑ Ordinarily, as frequency increases, the magnitude of the gain drops, and the phase becomes more negative, although these are only trends and may be reversed in particular frequency ranges. Unusual gain behavior can render the concepts of gain and phase margin inapplicable. Then other methods such as the Nyquist plot have to be used to assess stability.
 ↑ The critical amount of feedback where the peak in the gain just disappears altogether is the maximally flat or Butterworth design.
 ↑ The frequency where the openloop gain flips sign f_{180} does not change with a change in feedback factor; it is a property of the openloop gain. The value of the gain at f_{180} also does not change with a change in β. Therefore, we could use the previous values from Figures 6 and 7. However, for clarity the procedure is described using only Figures 8 and 9.
References
 ↑ R. K. Rao Yarlagadda (2010). Analog and Digital Signals and Systems. Springer Science & Business Media. p. 243. ISBN 9781441900340. https://archive.org/details/analogdigitalsig00yarl_849.
 ↑ Van Valkenburg, M. E. University of Illinois at UrbanaChampaign, "In memoriam: Hendrik W. Bode (19051982)", IEEE Transactions on Automatic Control, Vol. AC29, No 3., March 1984, pp. 193–194. Quote: "Something should be said about his name. To his colleagues at Bell Laboratories and the generations of engineers that have followed, the pronunciation is bohdee. The Bode family preferred that the original Dutch be used as bohdah."
 ↑ "Vertaling van postbode, NL>EN". mijnwoordenboek.nl. http://www.mijnwoordenboek.nl/vertalen.php?s1=&s2=NL+%3E+EN&s3=NL+%3E+EN&woord=postbode.
 ↑ David A. Mindell Between Human and Machine: Feedback, Control, and Computing Before Cybernetics JHU Press, 2004, ISBN:0801880572, pp. 127–131.
 ↑ Skogestad, Sigurd; Postlewaite, Ian (2005). Multivariable Feedback Control. Chichester, West Sussex, England: John Wiley & Sons, Ltd.. ISBN 047001167X.
 ↑ Thomas H. Lee (2004). The design of CMOS radiofrequency integrated circuits (2nd ed.). Cambridge UK: Cambridge University Press. pp. 451–453. ISBN 0521835399. http://worldcat.org/isbn/0521835399.
 ↑ William S. Levine (1996). The control handbook: the electrical engineering handbook series (2nd ed.). Boca Raton FL: CRC Press/IEEE Press. p. 163. ISBN 0849385709. https://books.google.com/books?id=2WQP5JGaJOgC&q=stability+%22minimum+phase%22&pg=RA1PA163.
 ↑ Invariance and Systems Theory: Algebraic and Geometric Aspects. New York, NY: SpringerVerlag. February 1981. ISBN 9783540105657.
 ↑ ^{9.0} ^{9.1} Willy M C Sansen (2006). Analog design essentials. Dordrecht, The Netherlands: Springer. pp. 157–163. ISBN 0387257462. http://worldcat.org/isbn/0387257462.
External links
 How to draw piecewise asymptotic Bode plots
 Gnuplot code for generating Bode plot: DINA4 printing template (pdf)
Original source: https://en.wikipedia.org/wiki/Bode plot.
Read more 