Step response

The step response of a system in a given initial state consists of the time evolution of its outputs when its control inputs are Heaviside step functions. In electronic engineering and control theory, step response is the time behaviour of the outputs of a general system when its inputs change from zero to one in a very short time. The concept can be extended to the abstract mathematical notion of a dynamical system using an evolution parameter.

From a practical standpoint, knowing how the system responds to a sudden input is important because large and possibly fast deviations from the long term steady state may have extreme effects on the component itself and on other portions of the overall system dependent on this component. In addition, the overall system cannot act until the component's output settles down to some vicinity of its final state, delaying the overall system response. Formally, knowing the step response of a dynamical system gives information on the stability of such a system, and on its ability to reach one stationary state when starting from another.

Time domain "versus" frequency domain

Depending on the application, instead of frequency response, system performance may be specified in terms of parameters describing time-dependence of response. The step response can be described by the following quantities related to its time behavior,

*overshoot
*rise time
*settling time
*ringing

In the case of linear dynamic systems, much can be inferred about the system from these characteristics. Below the step response of a simple two-pole amplifier is presented, and some of these terms are illustrated.

tep response of feedback amplifiers

This section describes the step response of a simple negative feedback amplifier shown in Figure 1. The feedback amplifier consists of a main open-loop amplifier of gain "A"_OL and a feedback loop governed by a feedback factor β. This feedback amplifier is analyzed to determine how its step response depends upon the time constants governing the response of the main amplifier, and upon the amount of feedback used.

Analysis

A negative feedback amplifier has gain given by (see negative feedback amplifier):

:

where "A"_OL = open-loop gain, "A"_FB = closed-loop gain (the gain with negative feedback present) and β = feedback factor. The step response of such an amplifier is easily handled in the case that the open-loop gain has two poles (two time constants, τ₁, τ₂), that is, the open-loop gain is given by:

:$A\_\{OL\}\; =\; frac\; \{A\_0\}\; \{(1+j\; omega\; au\_1)\; (1\; +\; j\; omega\; au\_2)\}\; ,$

with zero-frequency gain "A"₀ and angular frequency ω = 2π"f", which leads to the closed-loop gain:

:   •  

The time dependence of the amplifier is easy to discover by switching variables to "s" = "j"ω, whereupon the gain becomes:

:$A\_\{FB\}\; =\; frac\; \{A\_0\}\; \{\; au\_1\; au\_2\; \}$   •  

The poles of this expression (that is, the zeros of the denominator) occur at:

:$2s\; =\; -\; left(\; frac\; \{1\}\; \{\; au\_1\}\; +\; frac\; \{1\}\; \{\; au\_2\}\; ight)$ ::::

which shows for large enough values of β"A"₀ the square root becomes the square root of a negative number, that is the square root becomes imaginary, and the pole positions are complex conjugate numbers, either "s"₊ or "s"₋; see Figure 2:

:$s\_\{pm\}\; =\; -\; ho\; pm\; j\; mu\; ,$

with

::$ho\; =\; frac\; \{1\}\{2\}\; left(\; frac\; \{1\}\; \{\; au\_1\}\; +\; frac\; \{1\}\; \{\; au\_2\}\; ight\; )\; ,$

and::Using polar coordinates with the magnitude of the radius to the roots given by |"s"| (Figure 2):

:$|\; s\; |\; =\; |s\_\{\; pm\; \}\; |\; =\; sqrt\{\; ho^2\; +mu^2\}\; ,$

and the angular coordinate φ is given by:

: $mathrm\; \{cos\}\; phi\; =\; frac\; \{\; ho\}\; \{\; |\; s\; |\; \}$   $mathrm\; \{sin\}\; phi\; =\; frac\; \{\; mu\}\; \{\; |\; s\; |\; \}\; .$Tables of Laplace transforms show that the time response of such a system is composed of combinations of the two functions:

::$e^\{-\; ho\; t\}\; mathrm\; \{sin\}\; (\; mu\; t)$ $quad$ and $quad$ $e^\{-\; ho\; t\}\; mathrm\; \{cos\}\; (\; mu\; t)\; ,$

which is to say, the solutions are damped oscillations in time. In particular, the unit step response of the system is: cite book
author=Benjamin C Kuo & Golnaraghi F
title=Automatic control systems
year= 2003
pages=p. 253
publisher=Wiley
edition=Eighth Edition
location=New York
isbn=0-471-13476-7
url=http://worldcat.org/isbn/0-471-13476-7] :$S(t)\; =\; 1\; -\; e^\{-\; ho\; t\}\; frac\; \{\; mathrm\; \{sin\}\; left(\; mu\; t\; +\; phi\; ight)\}\{\; mathrm\; \{sin\}(\; phi\; )\}\; .$

Notice that the damping of the response is set by ρ, that is, by the time constants of the open-loop amplifier. In contrast, the frequency of oscillation is set by μ, that is, by the feedback parameter through β"A"₀. Because ρ is a sum of reciprocals of time constants, it is interesting to notice that ρ is dominated by the "shorter" of the two.

Results

Figure 3 shows the time response to a unit step input for three values of the parameter μ. It can be seen that the frequency of oscillation increases with μ, but the oscillations are contained between the two asymptotes set by the exponentials [ 1 - exp (−ρt) ] and [ 1 + exp (−ρt) ] . These asymptotes are determined by ρ and therefore by the time constants of the open-loop amplifier, independent of feedback.

The phenomena of oscillation about final value is called ringing. The overshoot is the maximum swing above final value, and clearly increases with μ. Likewise, the undershoot is the minimum swing below final value, again increasing with μ. The settling time is the time for departures from final value to sink below some specified level, say 10% of final value.

The dependence of settling time upon μ is not obvious, and the approximation of a two-pole system probably is not accurate enough to make any real-world conclusions about feedback dependence of settling time. However, the asymptotes [ 1 - exp (−ρt) ] and [ 1 + exp (−ρt) ] clearly impact settling time, and they are controlled by the time constants of the open-loop amplifier, particularly the shorter of the two time constants. That suggests that a specification on settling time must be met by appropriate design of the open-loop amplifier.

The two major conclusions from this analysis are:
#Feedback controls the amplitude of oscillation about final value for a given open-loop amplifier and given values of open-loop time constants, τ₁ and τ₂.
#The open-loop amplifier decides settling time. It sets the time scale of Figure 3, and the faster the open-loop amplifier, the faster this time scale.

As an aside, it may be noted that real-world departures from this linear two-pole model occur due to two major complications: first, real amplifiers have more than two poles, as well as zeros; and second, real amplifiers are nonlinear, so their step response changes with signal amplitude.

Control of overshoot

How overshoot may be controlled by appropriate parameter choices is discussed next.

Using the equations above, the amount of overshoot can be found by differentiating the step response and finding its maximum value. The result for maximum step response "S"_max is: cite book
author=Benjamin C Kuo & Golnaraghi F
title=p. 259
isbn=0-471-13476-7
url=http://worldcat.org/isbn/0-471-13476-7]

:$S\_\{max\}=\; 1$ $+\; mathrm\; \{exp\}\; left(\; -\; pi\; frac\; \{\; ho\; \}\{\; mu\; \}\; ight)\; .$

The final value of the step response is 1, so the exponential is the actual overshoot itself. It is clear the overshoot is zero if μ = 0, which is the condition:

:

This quadratic is solved for the ratio of time constants by setting "x" = ( τ₁ / τ₂ )^{1 / 2} with the result

:

Because β "A"₀ >> 1, the 1 in the square root can be dropped, and the result is

:

In words, the first time constant must be much larger than the second. To be more adventurous than a design allowing for no overshoot we can introduce a factor α in the above relation:

:

and let α be set by the amount of overshoot that is acceptable.

Figure 4 illustrates the procedure. Comparing the top panel (α = 4) with the lower panel (α = 0.5) shows lower values for α increase the rate of response, but increase overshoot. The case α = 2 (center panel) is the "maximally flat" design that shows no peaking in the Bode gain vs. frequency plot. That design has the rule of thumb built-in safety margin to deal with non-ideal realities like multiple poles (or zeros), nonlinearity (signal amplitude dependence) and manufacturing variations, any of which can lead to too much overshoot. The adjustment of the pole separation (that is, setting α ) is the subject of frequency compensation, and one such method is pole splitting.

Control of settling time

The amplitude of ringing in the step response in Figure 3 is governed by the damping factor exp ( −ρ t ). That is, if we specify some acceptable step response deviation from final value, say Δ, that is:

:$S(t)\; le\; 1\; +\; Delta\; ,$

this condition is satisfied regardless of the value of β "A"_OL provided the time is longer than the settling time, say "t"_S, given by: [This estimate is a bit conservative (long) because the factor 1 /sin(φ) in the overshoot contribution to "S" ( "t" ) has been replaced by 1 /sin(φ) ≈ 1.]

:$Delta\; =\; e^\{-\; ho\; t\_S\; \}$   or   $t\_S\; =\; frac\; \{\; mathrm\{ln\}\; left(\; frac\{1\}\; \{\; Delta\}\; ight)\; \}\; \{\; ho\; \}\; =\; au\_2\; frac\; \{2\; mathrm\{ln\}\; left(\; frac\{1\}\; \{\; Delta\}\; ight)\; \}\; \{\; 1\; +\; frac\; \{\; au\_2\; \}\; \{\; au\_1\}\; \}\; approx\; 2\; au\_2\; mathrm\{ln\}\; left(\; frac\{1\}\; \{\; Delta\}\; ight)\; ,$

where the approximation τ₁ >> τ₂ is applicable because of the overshoot control condition, which makes τ₁ = α β"A"_OL τ₂. Often the settling time condition is referred to by saying the settling period is inversely proportional to the unity gain bandwidth, because 1/( 2π τ₂ ) is close to this bandwidth for an amplifier with typical dominant pole compensation. However, this result is more precise than this rule of thumb. As an example of this formula, if Δ = 1/e⁴ = 1.8 %, the settling time condition is "t"_S = 8 τ₂.

In general, control of overshoot sets the time constant ratio, and settling time "t"_S sets τ₂. cite book
author=David A. Johns & Martin K W
title=Analog integrated circuit design
year= 1997
pages=pp. 234-235
publisher=Wiley
location=New York
isbn=0-471-14448-7
url=http://worldcat.org/isbn/0-471-14448-7] cite book
author=Willy M C Sansen
title=Analog design essentials
page=§0528 p. 163
year= 2006
publisher=Springer
location=Dordrecht, The Netherlands
isbn=0-387-25746-2
url=http://worldcat.org/isbn/0-387-25746-2] [According to Johns and Martin, "op. cit.", settling time is significant in switched-capacitor circuits, for example, where an op amp settling time must be less than half a clock period for sufficiently rapid charge transfer.]

Phase margin

Next, the choice of pole ratio τ₁ / τ₂ is related to the phase margin of the feedback amplifier. [The gain margin of the amplifier cannot be found using a two-pole model, because gain margin requires determination of the frequency "f"₁₈₀ where the gain flips sign, and this never happens in a two-pole system. If we know "f"₁₈₀ for the amplifier at hand, the gain margin can be found approximately, but "f"₁₈₀ then depends on the third and higher pole positions, as does the gain margin, unlike the estimate of phase margin, which is a two-pole estimate.] The procedure outlined in the Bode plot article is followed. Figure 5 is the Bode gain plot for the two-pole amplifier in the range of frequencies up to the second pole position. The assumption behind Figure 5 is that the frequency "f"_0dB lies between the lowest pole at "f"₁ = 1 / ( 2π τ₁ ) and the second pole at "f"₂ = 1 / ( 2π τ₂ ). As indicated in Figure 5, this condition is satisfied for values of α ≥ 1.

Using Figure 5 the frequency (denoted by "f"_0dB ) is found where the loop gain β"A"₀ satisfies the unity gain or 0 dB condition , as defined by:

:

The slope of the downward leg of the gain plot is (20 dB/decade); for every factor of ten increase in frequency, the gain drops by the same factor: [For more detail on the use of logarithmic scales, see log scale.] See also Figure 3 in pole splitting.]

:

The phase margin is the departure of the phase at "f"_0dB from −180°. Thus, the margin is:

:$phi\_m\; =\; 180\; ^circ\; -\; mathrm\; \{atan\}\; (f\_\{0dB\}\; /f\_1)\; -\; mathrm\; \{atan\}\; (\; f\_\{0dB\}\; /f\_2)\; .$

Because "f"_0dB / "f"₁ = β"A"₀ >> 1, the term in "f"₁ is 90°. That makes the phase margin:

:$phi\_m\; =\; 90\; ^circ\; -\; mathrm\; \{atan\}\; (\; f\_\{0dB\}\; /f\_2)$::::$=\; 90\; ^circ\; -\; mathrm\; \{atan\}\; left(\; frac\; \{1\}\; \{alpha\; \}\; ight)$ $=\; mathrm\; \{atan\}\; left(\; alpha\; ight)\; .$

In particular, for case α = 1, φ_m = 45°, and for α = 2, φ_m = 63.4°. Sansencite book
author=Willy M C Sansen
title=§0526 p. 162
isbn=0-387-25746-2
url=http://worldcat.org/isbn/0-387-25746-2] recommends α = 3, φ_m = 71.6° as a "good safety position to start with".

If α is increased by shortening τ₂, the settling time "t"_S also is shortened. If α is increased by lengthening τ₁, the settling time "t"_S is little altered. More commonly, both τ₁ "and" τ₂ change, for example if the technique of pole splitting is used.

As an aside, for an amplifier with more than two poles, the diagram of Figure 5 still may be made to fit the Bode plots by making "f"₂ a fitting parameter, referred to as an "equivalent second pole" position. cite book
author=Gaetano Palumbo & Pennisi S
title=Feedback amplifiers: theory and design
year= 2002
pages=§ 4.4 pp. 97-98
publisher=Kluwer Academic Press
location=Boston/Dordrecht/London
isbn=0-7923-7643-9
url=http://books.google.com/books?id=Xb0W1VsQFe0C&pg=PA98&dq=%22equivalent+two-pole+amplifier%22&lr=&as_brr=0&sig=LVB6t0wlg7WL7U_PB4eavGTP7Aw ]

Formal mathematical description

This section provides a formal mathematical definition of step response in terms of the abstract mathematical concept of a dynamical system $scriptstylemathfrak\{S\}$: all notations and assumptions required for the following description are listed here.
*$scriptstyle\; tin\; T$ is the evolution parameter of the system, called "time" for the sake of simplicity,
* is the state of the system at time $t,$, called "output" for the sake of simplicity,
*$scriptstylePhi:T\; imes\; Mlongrightarrow\; M$ is the dynamical system evolution function,
* is the dynamical system initial state,
*$scriptstyle\; H(t),$ is the Heaviside step function

Nonlinear dynamical system

For a general dynamical system, the step response is defined as follows:

::

It is the evolution function when the control inputs (or source term, or forcing inputs) are Heaviside functions: the notation emphasizes this concept showing $H(t)$ as a subscript.

Linear dynamical system

For a linear time-invariant black box, let $scriptstylemathfrak$S equiv S for notational convenience: the step response can be obtained by convolution of the Heaviside step function control and the impulse response "h (t)" of the system itself

::$a(t)\; =\; \{h*H\}(t)\; =\; \{H*h\}(t)\; =\; intlimits\_\{-infty\; \}^\{+infty\}!!\{h(\; au\; )H(t\; -\; au\; )\}\; d\; au\; =\; intlimits\_\{-infty\}^t!!\{h(\; au)\}d\; au\; .$

See also

*impulse response
*overshoot
*rise time
*settling time
*pole splitting
* [http://bcs.wiley.com/he-bcs/Books?action=resource&bcsId=2357&itemId=0471134767&resourceId=5596 Kuo power point slides; Chapter 7 especially]

References and notes

Further reading

*Robert I. Demrow "Settling time of operational amplifiers" [http://www.analog.com/UploadedFiles/Application_Notes/466359863287538299597392756AN359.pdf]
*Cezmi Kayabasi "Settling time measurement techniques achieving high precision at high speeds" [http://www.wpi.edu/Pubs/ETD/Available/etd-050505-140358/unrestricted/ckayabasi.pdf]

*Vladimir Igorevic Arnol'd "Ordinary differential equations", various editions from MIT Press and from Springer Verlag, chapter 1 "Fundamental concepts"