Skip to main content

Leapfrog simulation of 4th order Butterworth LPF

Short answer

Leapfrog simulation realizes the ladder equations using op-amp integrators. For the 4th-order Butterworth prototype (L1=0.7654, C2=1.848, L3=1.848, C4=0.7654) scaled to ω0=6283.2 rad/s (1 kHz) with R=10 kΩ, capacitors are A1/A5: 12.18 nF, A2/A4: 29.41 nF; two inverters complete the structure.

Design (simulate) the 4th order Butterworth LPF in doubly-terminated (resistively-terminated lossless) network using Leapfrog simulation. The necessary information is listed in the given table: Order (n) = 4 and LPF; R1 = 1; L1 = 0.7654; C2 = 1.848; L3 = 1.848; C4 = 0.7654; R2 = 1

Leapfrog (operational) simulation realises the doubly terminated LC ladder by simulating its equations rather than its elements: each reactive element is replaced by an op-amp integrator, and the integrators are coupled exactly as the ladder voltages and currents are coupled. The result keeps the low sensitivity of the passive ladder.

Given prototype (normalised to 1 rad/s, 1 Ω terminations): R1=1R_1 = 1, L1=0.7654L_1 = 0.7654 H, C2=1.848C_2 = 1.848 F, L3=1.848L_3 = 1.848 H, C4=0.7654C_4 = 0.7654 F, R2=1 ΩR_2 = 1\ \Omega.

 Rs=R1=1      L1=0.7654        L3=1.848
 o----[R1]-----(LLL)-----+-----(LLL)-----+-----+-----o Vo
                         |               |     |
                     C2=1.848     C4=0.7654  [R2=1]
                         |               |     |
 GND --------------------+---------------+-----+-----o

Write the state equations of the doubly terminated ladder (normalised, Rs=RL=1 ΩR_s = R_L = 1\ \Omega). Every current is turned into a voltage by multiplying by a 1 Ω resistance, so each equation becomes an integrator:

I1=Vin−V2sL1+1V2=I1−I3sC2I3=V2−V4sL3V4=I3sC4+1Vo=V4\begin{aligned} I_1 &= \frac{V_{in} - V_2}{sL_1 + 1} \\ V_2 &= \frac{I_1 - I_3}{sC_2} \\ I_3 &= \frac{V_2 - V_4}{sL_3} \\ V_4 &= \frac{I_3}{sC_4 + 1} \\ V_o &= V_4 \end{aligned}

The first and last equations are lossy integrators (the terminations give the "+1"); the middle ones are lossless integrators. Neighbouring integrators feed each other in alternate directions, which gives the "leapfrog" structure.

Using inverting op-amp integrators, Vo=−1sRC∑ViV_o = -\frac{1}{sRC}\sum V_i (lossless) or Vo=−∑Vi/RsC+1/RV_o = -\frac{\sum V_i/R}{sC + 1/R} (lossy, resistor RR across CC), the signs are arranged as follows. Inverters supply the signals of opposite sign:

Op-ampTypeInputs (each through RR)OutputNormalised CC (R = 1 Ω)
A1lossy integratorVinV_{in}, −V2-V_2−I1-I_10.7654 F
A2integrator−I1-I_1, I3I_3V2V_21.848 F
A3inverterV2V_2−V2-V_2R / R
A4integrator−V2-V_2, V4V_4I3I_31.848 F
A5lossy integratorI3I_3−V4-V_40.7654 F
A6inverter−V4-V_4V4=VoV_4 = V_oR / R

Check of one stage: A2 gives −(−I1)+I3sCA2=I1−I3s(1.848)-\frac{(-I_1) + I_3}{sC_{A2}} = \frac{I_1 - I_3}{s(1.848)}, which is the ladder equation for V2V_2.

Denormalisation. Frequency scale by ω0=2π(1000)=6283.2 rad/s\omega_0 = 2\pi(1000) = 6283.2\ \text{rad/s} and choose all resistors RR = 10 kΩ (impedance scale kz=104k_z = 10^4). Each capacitor becomes C=g/(ω0R)C = g/(\omega_0 R):

CA1=0.76546283.19×104=12.18 nFCA2=1.8486283.19×104=29.41 nFCA4=1.8486283.19×104=29.41 nFCA5=0.76546283.19×104=12.18 nF\begin{aligned} C_{A1} &= \frac{0.7654}{6283.19 \times 10^4} = \text{12.18 nF} \\ C_{A2} &= \frac{1.848}{6283.19 \times 10^4} = \text{29.41 nF} \\ C_{A4} &= \frac{1.848}{6283.19 \times 10^4} = \text{29.41 nF} \\ C_{A5} &= \frac{0.7654}{6283.19 \times 10^4} = \text{12.18 nF} \end{aligned}

Final element values:

Op-ampResistorsCapacitor
A1input(s) 10 kΩ, feedback 10 kΩ across C12.18 nF
A2input(s) 10 kΩ29.41 nF
A4input(s) 10 kΩ29.41 nF
A5input(s) 10 kΩ, feedback 10 kΩ across C12.18 nF
inverters10 kΩ in, 10 kΩ feedbacknone
 Lossy integrator (A1 and A5)
              +-----[ C ]-----+
              |               |
              +-----[ R ]-----+
              |               |
 Va --[R]-----+---|-\         |
 Vb --[R]-----+   |  >--------+---- Vo
              GND-|+/
 Lossless integrator (A2 and A4): the same circuit
 without the feedback resistor (only C from output to (-)).

The passband gain is 1/2 (−6 dB), the same as the doubly terminated prototype; the response is that of the passive ladder but with no inductors, and the low sensitivity of the ladder is kept.

Since no frequency is given, the normalised design (all R = 1 Ω, C = 0.7654, 1.848, 1.848, 0.7654 F) is the answer; the 1 kHz values above show how it is made practical.

Answer (1 kHz example): R = 10 kΩ everywhere, capacitors 12.18 nF, 29.41 nF, 29.41 nF, 12.18 nF.

✓ Answer checked by Claude Sonnet 5.5

All Filter Design chapter 8 questions →