Skip to main content

Tow-Thomas LPF derivation; realize −2000/(s²+500s+10⁶)

Short answer

For the Tow-Thomas low-pass biquad, choose equal capacitors C and equal resistors R2=R3=R. With ω0=1000 rad/s, Q=2, and DC gain magnitude 0.002, using C=1 μF: R=1/(ω0C)=1 kΩ, R1=QR=2 kΩ, R4=R2/0.002=500 kΩ, r=10 kΩ. Take the output from V3 for the negative sign.

Draw the circuit diagram of Tow Thomas low pass filter and derive its transfer function. Realize following low pass filter using Tow Thomas biquad circuit. T(s) = -2000/(s² + 500s + 1000000)

Tow-Thomas low-pass biquad

The Tow-Thomas biquad is a two-integrator-loop circuit using three op-amps: a lossy integrator, an ideal integrator and an inverter in a feedback loop. It gives low-pass and band-pass outputs at the same time, has low sensitivity, and lets ω0\omega_0, QQ and gain be tuned independently.

            R1 (sets Q)
          +--/\/\--+
          |   C1   |
          +---||---+
          |        |
Vin--R4---+--[OA1]-+--> V1 (band-pass)
          ^
          |  V1--R3-->[OA2 integrator, C2]--> V2 (LP)
          |  V2--r--->[OA3 inverter, r/r]---> V3=-V2
          |                                    |
          +------------------R2----------------+

All op-amp (+) inputs are grounded. OA1 is a lossy integrator (R1∥C1R_1 \parallel C_1 in feedback), OA2 an ideal integrator (C2C_2 in feedback), OA3 a unity-gain inverter.

Derivation (ideal op-amps, so each inverting input is a virtual ground):

KCL at the inverting input of OA1:

VinR4+V3R2+V1(1R1+sC1)=0\frac{V_{in}}{R_4} + \frac{V_3}{R_2} + V_1\left(\frac{1}{R_1} + sC_1\right) = 0

OA2 (integrator) and OA3 (inverter):

V2=−V1sR3C2,V3=−V2=V1sR3C2  ⇒  V1=sR3C2V3V_2 = -\frac{V_1}{sR_3C_2}, \qquad V_3 = -V_2 = \frac{V_1}{sR_3C_2} \;\Rightarrow\; V_1 = sR_3C_2V_3

Substituting V1V_1 into the KCL equation:

sR3C2V3(1R1+sC1)+V3R2=−VinR4V3(s2R3C1C2+sR3C2R1+1R2)=−VinR4\begin{aligned} sR_3C_2V_3\left(\frac{1}{R_1} + sC_1\right) + \frac{V_3}{R_2} &= -\frac{V_{in}}{R_4} \\ V_3\left(s^2R_3C_1C_2 + \frac{sR_3C_2}{R_1} + \frac{1}{R_2}\right) &= -\frac{V_{in}}{R_4} \end{aligned}

Dividing by R3C1C2R_3C_1C_2:

V3Vin=−1R3R4C1C2s2+sR1C1+1R2R3C1C2,V2Vin=−V3Vin\frac{V_3}{V_{in}} = \frac{-\dfrac{1}{R_3R_4C_1C_2}}{s^2 + \dfrac{s}{R_1C_1} + \dfrac{1}{R_2R_3C_1C_2}}, \qquad \frac{V_2}{V_{in}} = -\frac{V_3}{V_{in}}

Comparing with T(s)=±Kω02s2+(ω0/Q)s+ω02T(s) = \dfrac{\pm K\omega_0^2}{s^2 + (\omega_0/Q)s + \omega_0^2}:

ω0=1R2R3C1C2,ω0Q=1R1C1,K=R2R4 (dc gain)\omega_0 = \frac{1}{\sqrt{R_2R_3C_1C_2}}, \qquad \frac{\omega_0}{Q} = \frac{1}{R_1C_1}, \qquad K = \frac{R_2}{R_4}\ \text{(dc gain)}

V2V_2 is a non-inverting low-pass output and V3V_3 an inverting one; V1V_1 is a band-pass output.

Design equations (choose C1=C2=CC_1 = C_2 = C, R2=R3=RR_2 = R_3 = R):

R=1ω0C,R1=QR=Qω0C,R4=R2K,r=any convenient valueR = \frac{1}{\omega_0 C}, \qquad R_1 = QR = \frac{Q}{\omega_0 C}, \qquad R_4 = \frac{R_2}{K}, \qquad r = \text{any convenient value}

Realization of T(s)=−2000s2+500s+1,000,000T(s) = \dfrac{-2000}{s^2 + 500s + 1{,}000{,}000}

Step 1: Parameters.

ω02=106  ⇒  ω0=1000 rad/sω0Q=500  ⇒  Q=1000500=2dc gain=T(0)=−2000106=−0.002,so K=R2R4=0.002\begin{aligned} \omega_0^2 &= 10^6 \;\Rightarrow\; \omega_0 = 1000\ \text{rad/s} \\ \frac{\omega_0}{Q} &= 500 \;\Rightarrow\; Q = \frac{1000}{500} = 2 \\ \text{dc gain} &= T(0) = \frac{-2000}{10^6} = -0.002, \quad \text{so } K = \frac{R_2}{R_4} = 0.002 \end{aligned}

The gain is negative, so the output is taken at V3V_3 (the inverter output), which gives −1/(R3R4C1C2)s2+s/(R1C1)+1/(R2R3C1C2)-\dfrac{1/(R_3R_4C_1C_2)}{s^2 + s/(R_1C_1) + 1/(R_2R_3C_1C_2)}.

Step 2: Element values. Choose C1=C2=C=1 μC_1 = C_2 = C = 1\ \muF, R2=R3=RR_2 = R_3 = R (with C=0.1 μC = 0.1\ \muF, R4R_4 would be 5 MΩ\Omega, which is too large, so 1 μ1\ \muF is chosen).

R2=R3=R=1ω0C=11000×10−6=1 kΩR1=QR=2×1=2 kΩR4=R2K=1 kΩ0.002=500 kΩr=10 kΩ (inverter)\begin{aligned} R_2 = R_3 = R &= \frac{1}{\omega_0 C} = \frac{1}{1000 \times 10^{-6}} = 1\ \text{k}\Omega \\ R_1 &= QR = 2 \times 1 = 2\ \text{k}\Omega \\ R_4 &= \frac{R_2}{K} = \frac{1\ \text{k}\Omega}{0.002} = 500\ \text{k}\Omega \\ r &= 10\ \text{k}\Omega \ \text{(inverter)} \end{aligned}

Check: numerator 1R3R4C1C2=1103×5×105×10−12=2000\dfrac{1}{R_3R_4C_1C_2} = \dfrac{1}{10^3 \times 5\times10^5 \times 10^{-12}} = 2000; 1R1C1=12000×10−6=500\dfrac{1}{R_1C_1} = \dfrac{1}{2000\times10^{-6}} = 500; 1R2R3C2=106\dfrac{1}{R_2R_3C^2} = 10^6. All coefficients match.

Answer: C1=C2=1 μC_1 = C_2 = 1\ \muF, R1=2R_1 = 2 kΩ\Omega, R2=R3=1R_2 = R_3 = 1 kΩ\Omega, R4=500R_4 = 500 kΩ\Omega, r=10r = 10 kΩ\Omega, output at V3V_3 (OA3).

✓ Answer checked by Claude Sonnet 5.5

All Filter Design chapter 6 questions →