Skip to main content

Tow-Thomas LPF for ECG (gain 1.3)

Short answer

From poles –450 ± j893.03, ω0 = 1000 rad/s, Q = 1.111. With C = 0.1 µF, R2 = R3 = 10 kΩ, R1 = 11.11 kΩ, R4 = 7.692 kΩ, inverter resistors r = 10 kΩ.

Draw the circuit diagram of Tow-Thomas biquad circuit and derive its transfer function. For an ECG machine of a hospital, you are asked to design a low pass filter using Tow-Thomas biquad circuit with poles located at -450 ± j893.03 and DC gain of 1.3. Make sure the elements in the final circuit are of practically realizable values.

Tow-Thomas biquad circuit

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

Step 1: ω0\omega_0 and QQ from the poles. Poles s=−σ±jωds = -\sigma \pm j\omega_d with σ=450\sigma = 450, ωd=893.03\omega_d = 893.03:

ω02=σ2+ωd2=4502+893.032=1,000,002.6  ⇒  ω0≈1000 rad/sω0Q=2σ  ⇒  Q=ω02σ=1000900=1.111\begin{aligned} \omega_0^2 &= \sigma^2 + \omega_d^2 = 450^2 + 893.03^2 = 1{,}000{,}002.6 \;\Rightarrow\; \omega_0 \approx 1000\ \text{rad/s} \\ \frac{\omega_0}{Q} &= 2\sigma \;\Rightarrow\; Q = \frac{\omega_0}{2\sigma} = \frac{1000}{900} = 1.111 \end{aligned}

Required transfer function (dc gain K=1.3K = 1.3):

T(s)=1.3 ω02s2+2σs+ω02=1.3×106s2+900 s+106T(s) = \frac{1.3\,\omega_0^2}{s^2 + 2\sigma s + \omega_0^2} = \frac{1.3 \times 10^{6}}{s^2 + 900\,s + 10^{6}}

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}

Step 2: Element values. Choose C1=C2=C=0.1 μC_1 = C_2 = C = 0.1\ \muF (a standard value that gives resistors in the kΩ\Omega range).

R2=R3=R=1ω0C=11000×0.1×10−6=10 kΩR1=QR=1.111×10 kΩ=11.11 kΩR4=R2K=10 kΩ1.3=7.692 kΩr=10 kΩ (inverter, equal resistors)\begin{aligned} R_2 = R_3 = R &= \frac{1}{\omega_0 C} = \frac{1}{1000 \times 0.1\times10^{-6}} = 10\ \text{k}\Omega \\ R_1 &= QR = 1.111 \times 10\ \text{k}\Omega = 11.11\ \text{k}\Omega \\ R_4 &= \frac{R_2}{K} = \frac{10\ \text{k}\Omega}{1.3} = 7.692\ \text{k}\Omega \\ r &= 10\ \text{k}\Omega \ \text{(inverter, equal resistors)} \end{aligned}

A 1000 rad/s (about 159 Hz) cutoff with Q≈1.1Q \approx 1.1 passes the ECG band and removes higher-frequency noise and interference. Use 1% metal-film resistors (nearest standard values: 7.68 kΩ\Omega for R4R_4 and 11.0 kΩ\Omega for R1R_1) and low-noise op-amps.

Final circuit: Tow-Thomas biquad (figure above) with C1=C2=0.1 μFC_1 = C_2 = 0.1\ \mu\text{F}, R1=11.11 kΩR_1 = 11.11\ \text{k}\Omega, R2=R3=10 kΩR_2 = R_3 = 10\ \text{k}\Omega, R4=7.692 kΩR_4 = 7.692\ \text{k}\Omega, r=10 kΩr = 10\ \text{k}\Omega; take the output at V2V_2 (OA2 output) for a positive dc gain of 1.3.

Answer: C=0.1 μFC = 0.1\ \mu\text{F}, R1=11.11 kΩR_1 = 11.11\ \text{k}\Omega, R2=R3=10 kΩR_2 = R_3 = 10\ \text{k}\Omega, R4=7.692 kΩR_4 = 7.692\ \text{k}\Omega, r=10 kΩr = 10\ \text{k}\Omega (all practical values).

✓ Answer checked by Claude Sonnet 5.5

All Filter Design chapter 6 questions →