If you're seeing this message, it means we're having trouble loading external resources on our website.

If you're behind a web filter, please make sure that the domains *.kastatic.org and *.kasandbox.org are unblocked.

Main content

RLC natural response - derivation

A formal derivation of the natural response of the RLC circuit. Written by Willy McAllister.

Introduction

Let's take a deep look at the natural response of a resistor-inductor-capacitor circuit (RLC). This is the last circuit we'll analyze with the full differential equation treatment.
The RLC circuit is representative of real life circuits we can actually build, since every real circuit has some finite resistance. This circuit has a rich and complex behavior that finds application in many areas of electrical engineering.
The circuit for the RLC natural response.

What we're building to

We will model the RLC circuit with a 2nd-order linear differential equation with current, i, as the independent variable:
Ld2idt2+Rdidt+1Ci=0
The resulting characteristic equation is:
s2+RLs+1LC=0
We will solve for the roots of the characteristic equation using the quadratic formula:
s=R±R24L/C2L
By substituting variables α and ωo we can write s a little simpler as:
s=α±α2ωo2
where α=R2L, and ωo=1LC
α is called the damping factor and ωo is the resonant frequency.
We will solve an example RLC circuit with specific component values and discover what the current and voltages look like.

Strategy

We follow the same line of reasoning we used for solving the second-order LC circuit in an earlier article.
  1. Create a second-order differential equation based on the i-v equations for the R, L, and C components. We will use Kirchhoff's Voltage Law to build the equation.
  2. Make an informed guess at a solution. As usual, our guess will be an exponential function of the form Kest.
  3. Insert the proposed solution into the differential equation. The exponential terms will factor out and leave us with a characteristic equation in variable s.
  4. Find the roots of the characteristic equation. This time we will need to use the quadratic formula to solve for the roots.
  5. Find the constants by accounting for the initial conditions.
  6. Celebrate the solution.

Model the circuit with a differential equation

Circuit conditions just before the switch closes: The current is 0, and the capacitor is charged up to an initial voltage of V0 volts.
When the switch closes, the circuit looks like this (now with voltage labels on the inductor and resistor, vL and vR).
Circuit conditions just after the switch closes. We still have to find the current and voltage at t=0+. We work this out in the section titled Find the initial conditions.
For each individual element, we can write i-v equations.
vL=Ldidt
vR=iR
vC=1Cidt
We can write Kirchhoff's Voltage Law (KVL) starting in the lower left corner and summing voltages going around the loop clockwise. The inductor has a voltage rise, while the resistor and capacitor have voltage drops.
+vLvRvC=0
Replacing the v terms with the corresponding i terms gives us:
Ldidt+Ri+1Cidt=0
If we wanted to, we could attack this equation and try to solve it, but the integral term is awkward to deal with. We can retire the integral if we take the derivative of the entire equation.
ddt[Ldidt+Ri+1Cidt=0]
This gives us the following equation with a second derivative term, a first derivative term, and a plain i term, all still equal to 0.
Ld2idt2+Rdidt+1Ci=0
This is called a homogeneous second-order ordinary differential equation. It is homogeneous because every term is related to i and its derivatives. It is second order because the highest derivative is a second derivative. It is ordinary because there is only one independent variable (no partial derivatives). Now we set about solving our differential equation.

Propose a solution

Just like we did with previous natural response problems ( RC, RL, LC), we assume a solution with an exponential form. Exponential functions have the wondrous property that the derivatives look a lot like the original function. When you have multiple derivatives participating in a differential equation, it is really nice when they look alike. We assume a solution with this form:
i(t)=Kest
K is an adjustable parameter representing the amplitude of the current.
s is up in the exponent next to t, so it must represent some kind of frequency (s has to have units of 1/t). We call this the natural frequency.

Try the proposed solution

Next, substitute the proposed solution into the differential equation. If the equation turns out to be true, then our solution is a winner.
Ld2dt2Kest+RddtKest+1CKest=0
Now let's work on the terms with derivatives.
Middle term: The first derivative of the R term is
RddtKest=sRKest
Leading term: We take the derivative if the leading Kest term two times:
ddtKest=sKest
ddtsKest=s2Kest
so the leading term becomes:
Ld2dt2Kest=s2LKest
Plug these back into the differential equation:
s2LKest+sRKest+1CKest=0
Now we can factor out the common Kest term:
Kest(s2L+sR+1C)=0

Make the equation true

Now let's figure out how many ways we can make this equation true.
We could set K equal to 0. That means i=0 and we are putting nothing into the circuit and getting nothing out. Pretty boring.
The term est never becomes 0 unless we wait until t goes to . That's a long time from now. That leaves us with one interesting way to make the equation true: if the term with all the s's is zero.
s2L+sR+1C=0
This is called the characteristic equation of the LRC circuit.

Solve for the roots of the characteristic equation

Let's find values of s that make the characteristic equation true. (We want to find the roots of the characteristic equation.)
We have exactly the right tool for this, the quadratic formula:
For any quadratic equation: ax2+bx+c=0,
the quadratic formula gives us the roots (the zero-crossings):
x=b±b24ac2a
Looking back at the characteristic equation, we can plug in our circuit component values to get the roots. a=L, b=R, and c=1/C.
s=R±R24L/C2L
That's the answer for s, the natural frequency. We need to break this down a little more to get a feeling for the meaning of this solution.
We can make the notation a little more compact by replacing parts of the expression with two new variables, α and ωo.
α=R2L
ωo=1LC
Let me write the characteristic equation this way (dividing through by L):
s2+RLs+1LC=0
If we use α and ωo, the characteristic equation can be written:
s2+2αs+ωo2=0
We can revise the quadratic formula by mashing the 2L denominator up into each term of the numerator:
s=R2L±(R2L)2(4L/C4L2)
The second term under the square root reduces to:
(4L/C4L2)=(4L/C4L2)=1LC
And this lets us write s in terms of α and ωo as:
s=α±α2ωo2
We know s is some sort of frequency (it has to have units of 1/t). That means the two terms making up s are also some sort of frequency.
  • α is called the damping factor. It will determine how quickly the overall signal fades to zero.
  • ωo is called the resonant frequency. It will determine how fast the system swings back and forth. This is the same resonant frequency we found in the LC natural response.

Proposed solution, upgraded

The quadratic formula gave us two solutions for s, we'll call these s1 and s2. We need to include both of these in the proposed solution, so we update our proposed solution to be a linear combination (the superposition) of two separate exponential terms with four adjustable parameters:
i=K1es1t+K2es2t
s1 and s2 are natural frequencies,
K1 and K2 are amplitude terms.

Example circuit

At this point it's helpful to do a specific example with some actual component values, to see how one particular solution plays out. Here is our example circuit:
RLC natural response example. The capacitor has an initial voltage of 10 volts. There is no current flowing in the inductor at the moment the switch is closed.
The differential equation for the RLC circuit is
Ld2idt2+Rdidt+1Ci=0
With real component values it becomes:
1d2idt2+2didt+5i=0
As we always do, assume a solution of the form: i(t)=Kest
We perform the analysis we did above, which results in this characteristic equation:
s2+2s+5=0
Solving for the roots of the characteristic equation with the quadratic formula:
s=R±R24L/C2L
With real component values:
s=2±224152
s=2±4202
s=1±162
s=1±j2
(Electrical engineers use the letter j for the imaginary 1, since we use i as the symbol for current.)
We get a complex answer, just as we did with the LC natural response, only this time the complex answer includes both a real part and an imaginary part.
Solving for the roots of the characteristic equation gave us two possible answers for s, so the proposed solution for i is now written as the superposition of two different exponential terms:
i=K1e(1+j2)t+K2e(1j2)t
The terms up in the exponents are complex conjugates. Let's fuss around with the way this is written. We can tease apart the real and imaginary parts of the exponents:
i=K1e1te+j2t+K2e1tej2t,
and factor out the common e1t term:
i=et(K1e+j2t+K2ej2t)
Notice how the real part of s came through the factoring process to give us the leading term, a decaying exponential, et.
The terms in the parentheses are a sum of two imaginary exponentials where the exponents are complex conjugates. This looks just like what we saw in the LC natural response. As we did then, we call on Euler's formula to help us with these terms.

Euler's formula

Using Maclaurin series expansions for ejx, sinjx, and cosjx, it is possible to derive Euler's formula:
e+jx=cosx+jsinx
and
ejx=cosxjsinx
In the linked video, whenever Sal says i, we say j.
These formulas let us turn eimaginary into a normal complex number.

Use Euler's formula

We can use Euler's formula to transform the sum
K1e+j2t+K2ej2t
into
K1(cos2t+jsin2t)+K2(cos2tjsin2t).
Multiply through the constants K1 and K2:
K1cos2t+jK1sin2t+K2cos2tjK2sin2t,
and gather the cosine terms and sine terms:
(K1+K2)cos2t+j(K1K2)sin2t
Without messing up the equation, we can simplify how this appears if we replace the unknown K's with different unknown A's. Let A1=(K1+K2), and A2=j(K1K2).
The previous expression becomes:
A1cos2t+A2sin2t
And now put this back into our proposed solution:
i=et(A1cos2t+A2sin2t)
So far so good. Next, we need to figure out A1 and A2 using the initial conditions.

Find the initial conditions

For a second-order equation, you need two initial conditions to get a complete solution: one for the independent variable, i, and another for its first derivative, di/dt.
If we can figure out i and di/dt at a specific time, we can find A1 and A2.
Finding the initial conditions for RLC is pretty much the same as for the LC circuit. We just have account for the resistor.
Here's what we know about t=0 (the moment before the switch closes):
Circuit conditions just before the switch closes. At t=0
the current is 0 and the initial voltage on the capacitor is vC=10V.
  • The switch is open, so i(0)=0
  • The starting capacitor voltage is specified: vC(0)=V0
If t=0+ is the moment just after the switch closes, our goal is to find i(0+) and di/dt(0+). We know some properties of inductors and capacitors that tell us what happens when the switch closes, going from t=0 to t=0+:
  • Inductor current does not change instantly, so i(0+)=i(0)=0
  • Capacitor voltage does not change instantly, so vC(0+)=vC(0)=V0
Circuit conditions just after the switch closes, at t=0+. i(0+)=0 and vC(0+)=10V.
Now we know one initial condition, i(0+)=0, and we know something about the voltage, but we don't know di/dt(0+).
Let's go after the second initial condition, di/dt(0+). Any time I see di/dt it makes me think of the inductor i-v equation. If we can figure out the voltage across the inductor, we can figure out di/dt. Let's do that by process of elimination.
Kirchhoff's Voltage Law around the loop is:
+vLvRvC=0
Since i(0+)=0, that means the voltage across the resistor, vR, has to be 0. We also know the voltage across the capacitor is vC=V0. Let's fill in the KVL equation with these values:
vL0V0=0
And now we know the voltage across the inductor at t=0+:
vL=V0
We can use this to derive di/dt using the inductor i-v equation.
vL(0+)=Ldidt(0+)
10=1didt(0+)
didt(0+)=10A/sec
The moment just after the switch closes, the current in the inductor has an initial slope of 10 amperes per second.

Find constants A1 and A2 using the initial conditions

As a reminder, our proposed solution is:
i=et(A1cos2t+A2sin2t)
and the initial conditions are:
i(0+)=0
didt(0+)=10
If we evaluate i at t=0, we can find one of the A constants. Insert t=0 and i=0 into the proposed solution:
0=e0(A1cos20+A2sin20)
0=1(A1cos0+A2sin0)
0=(A11+A20)
A1=0
A1=0, so the cosine term drops out of the solution. One down, one to go. Our proposed solution now looks like:
i=A2etsin2t
Let's chase down A2 using the second initial condition:
We need an equation for the derivative of i. Where might we get such a thing? How about we take the derivative of the proposed solution?
didt=ddt(A2etsin2t)
The proposed solution is the product of two functions. To take its derivative we use the product rule:
(fg)=fg+fg
Identify the two parts of the product and their derivatives:
f=A2etg=sin2t
f=A2etg=2cos2t
Assemble the parts according to the product rule:
didt=A2etsin2t+A2et2cos2t
didt=A2et(2cos2tsin2t)
We can evaluate this expression at t=0+:
10=A2e0(2cos0sin0)
10=A21(20)=2A2
A2=5

Solution for current

And finally, after a bunch of hard work, the solution for current is:
i=5etsin2t
The graph of i as a function of time looks like this:
Natural response of an RLC circuit, R=2Ω, L=1H, and C=15F. The faint curves plot ±5et, the envelope of the decaying sine wave.
When the switch closes, the current takes a big surge upwards and takes on the shape of the first hump of a sine wave. The sine wave quickly fades away after a few swings because the energy in the system rapidly dissipates as heat as charge flows back and forth through the resistor.
The role of "friction" in this example, as played by the resistor value, represents a fairly high rate of energy dissipation. The current visibly changes sign only two times before settling to zero.
This is an example of an underdamped solution. We will introduce this descriptive term in the next section.

Solve the voltages

There is only one current in the circuit. Now that we know the natural response of the current, we can find the natural response of the three voltages.

Resistor voltage

We use Ohm's Law to find the resistor voltage: (there's a sign because i is backwards relative to vR)
vR=iR
vR=5etsin2t2Ω
vR=10etsin2t

Inductor voltage

The inductor voltage emerges from the inductor i-v equation:
vL=Ldidt
vL=1ddt(5etsin2t)
vL=5et(sin2t2cos2t)

Capacitor voltage

To find the capacitor voltage we can use the integral form of the capacitor i-v equation: (another extra sign because of the direction of i is reversed relative to vC)
vC=1Cidt
vC=11/55etsin2tdt
vC=5et(sin2t+2cos2t)
Here are all three voltages plotted together:

Summary

The RLC circuit is the electronic equivalent of a swinging pendulum with friction. The circuit can be modeled by this 2nd-order linear differential equation:
Ld2idt2+Rdidt+1Ci=0
The resulting characteristic equation is:
s2+RLs+1LC=0
We solved for the roots of the characteristic equation using the quadratic formula:
s=R±R24L/C2L
By substituting variables α and ωo we wrote s a little simpler as:
s=α±α2ωo2
where α=R2L and ωo=1LC
We finished up by solving an example circuit whose components produced a current (and voltages) that swing back and forth a few times.
The roots of the characteristic equation can take on both real and complex forms, depending on the relative size of α and ωo. In the next article, we will describe these three forms in additional detail:
  • overdamped, α>ω0, leads to the sum of two decaying exponentials
  • critically damped, α=ω0, leads to t decaying exponential
  • underdamped, α<ω0, leads to a decaying sine

Want to join the conversation?

  • blobby green style avatar for user Mark  Freiheit
    You guys should really make a video on this topic! Also, it is much easier to find v(0+), i(0+), dv(0+)/dt, di(0+)/dt, v(infinity), i(infinity) right off the bat so you don't have to pound through the problem the long way.

    Oh, and the alpha and omega for this only works if they are in series. You need to do a write up to show how to solve if they are in parallel.

    Thanks
    (25 votes)
    Default Khan Academy avatar avatar for user
  • blobby green style avatar for user ksanchezlacroix
    Why is the voltage across the capacitor (vc(t)) equals to the integral of the opposite of the current ? where does this minus sign come from ?
    (5 votes)
    Default Khan Academy avatar avatar for user
    • aqualine ultimate style avatar for user Ali Jaya Meilio
      Because, in this article, it is assumed that the positive sign of the capacitor is on the top, so, the direction of current that it used is different from the current that supposed to be in the circuit... so that's why it got the minus sign...

      If you flip the sign of the capacitor, it could be also worked though... so we don't get the minus sign, but in the KVL equation, it will be regarded as voltage rise, and it will get the plus sign... and in the end... it gets the same equation...
      (5 votes)
  • blobby green style avatar for user megaultrachickenrawr
    Hi,
    When finding A2 in derivative of the proposed solution, ​di​/dt​​=​​​d​​(e​^(−t​​)A​2​​sin2t)/dt
    Shouldn't the answer be e^(-t) (4 cos(2 t) - 2 sin(2 t)) since we are using the product rule? With the answer coming out as A2=10/4?
    (2 votes)
    Default Khan Academy avatar avatar for user
    • spunky sam orange style avatar for user Willy McAllister
      Thanks for pointing how where I skipped the product rule step when taking the derivative. I will ask to have that fixed. Here's how it should go:

      di/dt = d/dt ( e^-t A2 sin2t) = e^-t A2 (2cos2t - sin 2t)

      Evaluate this at t=0,

      10 = e^0 A2 (2cos 0 - sin 0)

      10 = 1 . A2 . (2 - 0) = 2 . A2

      A2 = 5
      (1 vote)
  • blobby green style avatar for user Emmen Bramasta
    How do you get the simpler s equation (the one with the variables α and ω​naught) from the quadratic formula?
    (1 vote)
    Default Khan Academy avatar avatar for user
    • spunky sam orange style avatar for user Willy McAllister
      Review the definitions of alpha and ωo. This happens right after it says, "We can make the notation a little more compact by replacing parts of the expression with two new variables, alpha and ωo:"

      Now look at the equation right after the sentence, "We can revise the quadratic formula by mashing the 2L denominator up into each term of the numerator:".

      The first part of the equation, before the +/- sign, looks like -alpha.
      The second part of the equation, under the square root symbol you find two terms, they look like alpha^2 and ωo^2.

      Please let me know if this helped you understand. This is a spot where the article needs to be really clear.
      (3 votes)
  • blobby green style avatar for user pushkarm27
    How can we include imaginary constant (j) in A2 and still get A2 as real (=5)?
    (2 votes)
    Default Khan Academy avatar avatar for user
    • spunky sam orange style avatar for user Willy McAllister
      Back in the section Use Euler's Formula, A2 was defined in terms of the K constants.

      A2 = j(K1 - K2)

      Then we pressed ahead and figured out A2 = 5.

      If you go back to the definition of A2 in terms of K's, that means the K's have to be imaginary numbers. They have to include a j term to get rid of the j in A2 = j(K1 - K2).
      (1 vote)
  • leaf green style avatar for user Frank Gu
    When you assume a solution in the exponential form in the "Propose a solution" section, that lets you factor out the exponential and set the other, quadratic factor to 0. Then you use the solutions for s in the "Proposed solution, upgraded" section to make a new form for the solution as a sum of 2 exponentials. But why were you able to do that when you assumed the form of the solution to be in the form of a single exponential to begin with? Wouldn't that be contradicting the assumption?
    (2 votes)
    Default Khan Academy avatar avatar for user
  • purple pi purple style avatar for user yssh
    Do we have to do all this calculus everytime we make a RLC circuit, And this will get more complicated with complicated RLC circuits. There is no general formula for the RLC like-LR or RC circuits?
    I wanted to make a function generator using a buck-boost converter, To predict the values of current and voltage at the given continuously changing PWM signal will be hell lot of work!
    (1 vote)
    Default Khan Academy avatar avatar for user
    • old spice man green style avatar for user Willy McAllister
      Most RLC circuits are analyzed with Laplace Transform theory. That topic is beyond the scope of this introductory KA class, but it is the standard method for working with frequency filters and control systems.

      The RLC circuit analyzed here is the parallel form. The solution to the natural response emerges from this long analysis. The answer is not a simple single exponential equation like we get for RC, but rather a choice of three different responses ("variations") depending on the value of R, L , and C. I did an improved writeup of this at spinningnumbers dot com. Look for the RLC articles in the Natural and Forced Response section.
      (2 votes)
  • mr pink red style avatar for user manavpandya31
    What role does s(natural frequency) play here ? What it determines, can you please explain about it ?
    (1 vote)
    Default Khan Academy avatar avatar for user
    • spunky sam orange style avatar for user Willy McAllister
      In the RLC derivation we found that s could have real or complex values. We described s in terms of alpha and omega, and found there were three general classes of behavior (under-damped, critically-damped, and over-damped). The solutions are rich and wonderful and varied. You get a drooping behavior from the real parts, and a wiggly behavior from the imaginary parts.

      If you want to think about the meaning of s, it is a parameter (a complex number) that captures both the natural droopy and wiggly aspects of current and voltage as they change with time. The units of s are 1/t, which is the unit of frequency. The real part of s specifies the rate of droop while the imaginary part of s tells you the rate of wiggle.
      (2 votes)
  • blobby green style avatar for user Rachel
    Why do we have 2L as the denominator in the quadratic formula?

    Rachel
    (1 vote)
    Default Khan Academy avatar avatar for user
    • old spice man green style avatar for user Willy McAllister
      Our differential equation for the RLC can be turned into a "characteristic equation". It turns out the characteristic equation follows the general form of a quadratic equation...

      ax^2 + bx + c = 0

      which can be solved with the quadratic formula...

      x = -b +/- sqrt(b^2 -4ac) / 2a

      In our case the RLC circuit assigns a, b, and c as...
      a = L, b = R, and c = 1/C

      The 2 comes from the quadratic formula.
      (1 vote)
  • blobby green style avatar for user Rachel
    Please , how do I find the value of R in terms of L and C that will result in repeated roots to be able to demonstrate an underdamped, overdamped and critically damped response in a RLC combination circuit.
    (1 vote)
    Default Khan Academy avatar avatar for user