# Circuit poles from the Cochrun-Grabel method

In the first of this series of related articles on circuit dynamics, “Circuit Dynamics for Design”, the groundwork was laid in polynomial algebra for the development of some simple methods for finding circuit poles and zeros. One of them, the Cochrun-Grabel method, is presented here. It is a tabular method that in manual use is often simple enough to not even require a table. Though more powerful methods have emerged, this one is so easy to use it is worth knowing. It is quite useful for multistage amplifiers with interacting stages.

**Cochrun-Grabel Method**

An improvement in accuracy over the OCTC method was introduced by Cochrun and Grabel in 1973 (“A Method for the Determination of the Transfer Function of Electronic Circuits”, Basil L. Cochrun, Arvin Grabel, *IEEE Transactions on Circuit Theory*, Vol. CT-20, No. 1, JAN73). It proceduralizes determination of the polynomial coefficients of *D*(*s*) for circuits with *n* capacitors but does not include zeros. It is based on an expansion of *s*-domain pole polynomial *D*(*s*) as in the OCTC method, where the linear term is the sum of the OCTCs. However, the method does not stop there and includes the coefficients of all the powers of *s*. In normalized form,

*a*_{0} = 1

Designate the OCTCs as *τ _{i}*,

*i*= 1,...,

*n*. Then the coefficient of the linear term is

The coefficient of *s*^{2} is

where *τ _{m}*

_{;k}is

*τ*of the

*C*port with the port of

_{m}*C*shorted:

_{k}*τ*

_{m}_{;shorted}with the remaining ports open. For the

*τ*

_{1}factor, find the time constants of the other ports with port 1 shorted. Then short port 2 and do the rest of the ports. For each successive

*τ*, the number of product terms decreases by one. For

_{i}*n*= 2, it is the product of two time constants, an OCTC and the other with the OCTC port shorted.

For *n* = 3,

By the third degree, the general procedure is quite cumbersome and non-intuitive:

For *n* = 3, it is not nearly so daunting;

For a large (> 3) value of *n*, keeping track of the combinations of terms in each coefficient can be arduous, and Sol Rosenstark streamlined the method by tabulating it in his book, *Feedback Amplifier Principles* (Macmillan, 1986). (I also describe it in *Designing High-Performance Amplifiers *(D. Feucht, SciTech, 2010); www.scitechpub.com). The use of *Rosenstark tables* is central to a refined Cochrun-Grabel procedure. The tables are triangular in form; the general second- and third-degree tables are shown below.

The columns are added to give the *a _{m}* after multiplying each column entry by the time constants in the columns to its left. The first (leftmost) column is

*a*

_{1}, the sum of the OCTCs. The second collects the combinations of

*a*

_{2}. Each column entry is then multiplied by the time constants in each of the columns to its left before adding it to the other (similarly-multiplied) column entries. Rosenstark truncates the tables at

*a*

_{2}or

*a*

_{3}for tractable approximations to

*D*(

*s*). While the Cochrun-Grabel method is not difficult to learn using Rosenstark tables, it is easy to retain intuitively in the case of quadratic or cubic polynomials. Most of the work is in finding the OCTCs. Because the order of choosing capacitors is arbitrary there are two possible ways of computing the quadratic coefficient and they are equivalent. It is demonstrated here for the textbook CE stage (with

*R*= 0 Ω), and using

_{E}*r*

_{π}= (

*β*

_{0}+ 1)·

*r*. First, find the OCTCs for both

_{e}*C*and

_{e}*C*:

_{c}*R _{be}* =

*R*||

_{B}*r*

_{π};

*R*=

_{bc}*R*+

_{L}*R*·(1 +

_{be}*R*/

_{L}*r*)

_{m}Now short *C _{c}* and find

*R _{e}*

_{;c}=

*R*||(

_{be}*R*||

_{L}*r*)

_{m}A shorted *C _{c}* places the dependent current source across the voltage upon which it is dependent. By the substitution theorem, if a current source of

*v*/

_{be}*r*is placed across

_{m}*v*, then it is equivalent to a resistance of

_{be}*v*/(

_{be}*v*/

_{be}*r*) =

_{m}*r*.

_{m}Next, short *C _{e}* instead and find the resistance across the

*b-c*port with

*b*-

*e*port shorted:

*R _{c}*

_{;e}=

*R*

_{L}We only need to short one of the two capacitances for the procedure, but because the order of the capacitors does not matter, there are two ways of working the method and both are given here.

The linear coefficient, *b*, is the sum of the OCTCs:

The quadratic coefficient can be found in two ways:

or

After much algebra, the two expressions for *a* are found to be* *equivalent*.* In both cases, the OCTC of the first C is multiplied by the SCTC of the second capacitor with the first C shorted. For either alternative, the resulting pole polynomial is

*s*^{2}·(*τ _{e}*·

*τ*

_{c}_{;e}) +

*s*·(

*τ*+

_{e}*t*) + 1 =

_{c}*s*

^{2}·

*a*+

*s*·

*b*+ 1

**Next: Cascaded CE Stages**

“Proof” that the Laplace Transform is Nothing

Transistor and PWM-Switch Analogs

Is there a high-frequency μ transform?

Advanced Technology: Flying Saucers & the Underworld

Advanced Technology: Beginning with Roswell

Is Technology More Advanced than it Appears?

Accidental engineering: 10 mistakes turned into innovation

That 60W-equivalent LED: What you don’t know, and what no one will tell you…

6 famous people you may not know are engineers

DC distribution in your house and 42-V cars

10 tips for a successful engineering resume

The 5 greatest engineers of all time

Higgs Pt. 9: What makes King Carl XVI Gustaf think it’s the Higgs Boson?

10 things you may not know about Tesla

Analog Fundamentals: Instrumentation for impedance measurement

What are the most useful sensors?

Welcome back to The EMC Blog, now with Magic!

Browser suffers from JavaScript-creep disease

New approaches to dominate in embedded development

Do your LEDs have this defect?

National Instruments announces world’s first mmWave software-defined 5G transceiver – all in PXI

Free PCB ECAD: The Ultimate list

Megavolt/kiloamp tests reveal extreme engineering challenges