Network synthesis

This is a good article. Click here for more information.
Source: Wikipedia, the free encyclopedia.

Network synthesis is a design technique for linear electrical circuits. Synthesis starts from a prescribed impedance function of frequency or frequency response and then determines the possible networks that will produce the required response. The technique is to be compared to network analysis in which the response (or other behaviour) of a given circuit is calculated. Prior to network synthesis, only network analysis was available, but this requires that one already knows what form of circuit is to be analysed. There is no guarantee that the chosen circuit will be the closest possible match to the desired response, nor that the circuit is the simplest possible. Network synthesis directly addresses both these issues. Network synthesis has historically been concerned with synthesising passive networks, but is not limited to such circuits.

The field was founded by Wilhelm Cauer after reading Ronald M. Foster's 1924 paper A reactance theorem. Foster's theorem provided a method of synthesising LC circuits with arbitrary number of elements by a partial fraction expansion of the impedance function. Cauer extended Foster's method to RC and RL circuits, found new synthesis methods, and methods that could synthesise a general RLC circuit. Other important advances before World War II are due to Otto Brune and Sidney Darlington. In the 1940s Raoul Bott and Richard Duffin published a synthesis technique that did not require transformers in the general case (the elimination of which had been troubling researchers for some time). In the 1950s, a great deal of effort was put into the question of minimising the number of elements required in a synthesis, but with only limited success. Little was done in the field until the 2000s when the issue of minimisation again became an active area of research, but as of 2023, is still an unsolved problem.

A primary application of network synthesis is the design of

directional couplers, and equalisation. In the 2000s, network synthesis began to be applied to mechanical systems as well as electrical, notably in Formula One
racing.

Overview

Network synthesis is all about designing an electrical network that behaves in a prescribed way without any preconception of the network form. Typically, an

(C). Such networks always have an impedance, denoted , in the form of a
complex frequency variable s. That is, the impedance is the ratio of two polynomials in s.[1]

There are three broad areas of study in network synthesis; approximating a requirement with a rational function, synthesising that function into a network, and determining equivalents of the synthesised network.[2]

Approximation

The idealised prescribed function will rarely be capable of being exactly described by polynomials. It is therefore not possible to synthesise a network to exactly reproduce it.

piecewise continuous response is impossible to represent with polynomials because of the discontinuities. To overcome this difficulty, a rational function is found that closely approximates the prescribed function using approximation theory.[4] In general, the closer the approximation is required to be, the higher the degree of the polynomial and the more elements will be required in the network.[5]

There are many polynomials and functions used in network synthesis for this purpose. The choice depends on which parameters of the prescribed function the designer wishes to optimise.

Chebyshev approximation in which the designer specifies how much the passband response can deviate from the ideal in exchange for improvements in other parameters.[8] Other approximations are available for optimising time delay, impedance matching, roll-off, and many other requirements.[9]

Realisation

Given a rational function, it is usually necessary to determine whether the function is realisable as a discrete passive network. All such networks are described by a rational function, but not all rational functions are realisable as a discrete passive network.[10] Historically, network synthesis was concerned exclusively with such networks. Modern active components have made this limitation less relevant in many applications,[11] but at the higher radio frequencies passive networks are still the technology of choice.[12] There is a simple property of rational functions that predicts whether the function is realisable as a passive network. Once it is determined that a function is realisable, there a number of algorithms available that will synthesise a network from it.[13]

Equivalence

A network realisation from a rational function is not unique. The same function may realise many equivalent networks. It is known that affine transformations of the impedance matrix formed in mesh analysis of a network are all impedance matrices of equivalent networks (further information at Analogue filter § Realisability and equivalence).[14] Other impedance transformations are known, but whether there are further equivalence classes that remain to be discovered is an open question.[15]

A major area of research in network synthesis has been to find the realisation which uses the minimum number of elements. This question has not been fully solved for the general case,[16] but solutions are available for many networks with practical applications.[17]

History

Wilhelm Cauer

The field of network synthesis was founded by German mathematician and scientist Wilhelm Cauer (1900–1945). The first hint towards a theory came from American mathematician Ronald M. Foster (1896–1998) when he published A reactance theorem in 1924. Cauer immediately recognised the importance of this work and set about generalising and extending it. His thesis in 1926 was on "The realisation of impedances of prescibed frequency dependence" and is the beginning of the field. Cauer's most detailed work was done during World War II, but he was killed shortly before the end of the war. His work could not be widely published during the war, and it was not until 1958 that his family collected his papers and published them for the wider world. Meanwhile, progress had been made in the United States based on Cauer's pre-war publications and material captured during the war.[18]

English self-taught mathematician and scientist

doctoral thesis.[21]

Raoul Bott

Foster's realisation was limited to LC networks and was in one of two forms; either a number of series LC circuits in parallel, or a number of parallel LC circuits in series. Foster's method was to expand into

partial fractions
. Cauer showed that Foster's method could be extended to RL and RC networks. Cauer also found another method; expanding as a
ideal transformers in their realisations of RLC networks. Having to include transformers is undesirable in a practical implementation of a circuit.[23]

A method of realisation that did not require transformers was provided in 1949 by Hungarian-American mathematician Raoul Bott (1923–2005) and American physicist Richard Duffin (1909–1996).[24] The Bott and Duffin method provides an expansion by repeated application of Richards' theorem, a 1947 result due to American physicist and applied mathematician Paul I. Richards (1923–1978).[25] The resulting Bott-Duffin networks have limited practical use (at least for rational functionals of high degree) because the number of components required grows exponentially with the degree.[26] A number of variations of the original Bott-Duffin method all reduce the number of elements in each section from six to five, but still with exponentially growing overall numbers.[27] Papers achieving this include Pantell (1954), Reza (1954), Storer (1954) and Fialkow & Gest (1955).[28] As of 2010, there has been no further significant advance in synthesising rational functions.[29]

In 1939, American electrical engineer

Boolean algebra of digital circuits.[34]

Cauer used

Butterworth polynomials.[37] Butterworth's work was entirely independent of Cauer, but it was later found that the Butterworth polynomials were a limiting case of the Chebyshev polynomials.[38] Even earlier (1929) and again independently, American engineer and scientist Edward Lawry Norton (1898–1983) designed a maximally-flat mechanical filter with a response entirely analogous to Butterworth's electrical filter.[39]

In the 2000s, interest in further developing network synthesis theory was given a boost when the theory started to be applied to large mechanical systems.[40] The unsolved problem of minimisation is much more important in the mechanical domain than the electrical due to the size and cost of components.[41] In 2017, researchers at the University of Cambridge, limiting themselves to considering biquadratic rational functions, determined that Bott-Duffin realisations of such functions for all series-parallel networks and most arbitrary networks had the minimum number of reactances (Hughes, 2017). They found this result surprising as it showed that the Bott-Duffin method was not quite so non-minimal as previously thought.[42] This research partly centred on revisiting the Ladenheim Catalogue. This is an enumeration of all distinct RLC networks with no more than two reactances and three resistances. Edward Ladenheim carried out this work in 1948 while a student of Foster. The relevance of the catalogue is that all these networks are realised by biquadratic functions.[43]

Applications

The single most widely used application of network synthesis is in the design of

network synthesis filter.[44]

Hendrik Bode

Another application is the design of

resistance then it is only possible to achieve a perfect match at a number of discrete frequencies; the match over the band as a whole must be approximated.[45] The designer first prescribes the frequency band over which the matching network is to operate, and then designs a band-pass filter for that band. The only essential difference between a standard filter and a matching network is that the source and load impedances are not equal.[46]

There are differences between filters and matching networks in which parameters are important. Unless the network has a dual function, the designer is not too concerned over the behaviour of the impedance matching network outside the passband. It does not matter if the transition band is not very narrow, or that the stopband has poor attenuation. In fact, trying to improve the bandwidth beyond what is strictly necessary will detract from the accuracy of the impedance match. With a given number of elements in the network, narrowing the design bandwidth improves the matching and vice versa. The limitations of impedance matching networks were first investigated by American engineer and scientist Hendrik Wade Bode in 1945, and the principle that they must necessarily be filter-like was established by Italian-American computer scientist Robert Fano in 1950.[47] One parameter in the passband that is usually set for filters is the maximum insertion loss. For impedance matching networks, a better match can be obtained by also setting a minimum loss. That is, the gain never rises to unity at any point.[48]

Time-delay networks can be designed by network synthesis with filter-like structures. It is not possible to design a delay network that has a constant delay at all frequencies in a band. An approximation to this behaviour must be used limited to a prescribed bandwidth. The prescribed delay will occur at most at a finite number of spot frequencies. The Bessel filter has maximally-flat time-delay.[49]

The application of network synthesis is not limited to the electrical domain. It can be applied to systems in any energy domain that can be represented as a network of linear components. In particular, network synthesis has found applications in mechanical networks in the mechanical domain. Consideration of mechanical network synthesis led Malcolm C. Smith to propose a new mechanical network element, the inerter, which is analogous to the electrical capacitor.[50] Mechanical components with the inertance property have found an application in the suspensions of Formula One racing cars.[51]

Synthesis techniques

Synthesis begins by choosing an approximation technique that delivers a rational function approximating the required function of the network. If the function is to be implemented with passive components, the function must also meet the conditions of a positive-real function (PRF).[52] The synthesis technique used depends in part on what form of network is desired, and in part how many kinds of elements are needed in the network. A one-element-kind network is a trivial case, reducing to an impedance of a single element. A two-element-kind network (LC, RC, or RL) can be synthesised with Foster or Cauer synthesis. A three-element-kind network (an RLC network) requires more advanced treatment such as Brune or Bott-Duffin synthesis.[53]

Which, and how many kinds of, elements are required can be determined by examining the

poles and zeroes (collectively called critical frequencies) of the function.[54]
The requirement on the critical frequencies is given for each kind of network in the relevant sections below.

Foster synthesis

Foster's synthesis, in its original form, can be applied only to LC networks. A PRF represents a two-element-kind LC network if the critical frequencies of all exist on the axis of the complex plane of (the

conjugate pairs
. must be the ratio of an even and odd polynomial and their degrees must differ by exactly one. These requirements are a consequence of
Foster's reactance theorem.[55]

Foster I form

Example of a Foster I form realisation

Foster's first form (Foster I form) synthesises as a set of parallel LC circuits in series. For example,

can be expanded into partial fractions as,

The first term represents a series inductor, a consequence of having a pole at infinity. If it had had a pole at the origin, that would represent a series capacitor. The remaining two terms each represent conjugate pairs of poles on the axis. Each of these terms can be synthesised as a parallel LC circuit by comparison with the impedance expression for such a circuit,[56]

The resulting circuit is shown in the figure.

Foster II form

Example of a Foster II form realisation

Foster II form synthesises as a set of series LC circuits in parallel. The same method of expanding into partial fractions is used as for Foster I form, but applied to the admittance, , instead of . Using the same example PRF as before,

Expanded in partial fractions,

The first term represents a shunt inductor, a consequence of having a pole at the origin (or, equivalently, has a zero at the origin). If it had had a pole at infinity, that would represent a shunt capacitor. The remaining two terms each represent conjugate pairs of poles on the axis. Each of these terms can be synthesised as a series LC circuit by comparison with the admittance expression for such a circuit,[57]

The resulting circuit is shown in the figure.

Extension to RC or RL networks

Foster synthesis can be extended to any two-element-kind network. For instance, the partial fraction terms of an RC network in Foster I form will each represent an R and C element in parallel. In this case, the partial fractions will be of the form,[58]

Other forms and element kinds follow by analogy. As with an LC network, The PRF can be tested to see if it is an RC or RL network by examining the critical frequencies. The critical frequencies must all be on the negative real axis and alternate between poles and zeroes, and there must be an equal number of each. If the critical frequency nearest, or at, the origin is a pole, then the PRF is an RC network if it represents a , or it is an RL network if it represents a . Vice versa if the critical frequency nearest, or at, the origin is a zero. These extensions of the theory also apply to the Cauer forms described below.[59]

Immittance

In the Foster synthesis above, the expansion of the function is the same procedure in both the Foster I form and Foster II form. It is convenient, especially in theoretical works, to treat them together as an immittance rather than separately as either an impedance or an admittance. It is only necessary to declare whether the function represents an impedance or an admittance at the point that an actual circuit needs to be realised. Immittance can also be used in the same way with the Cauer I and Cauer II forms and other procedures.[60]

Cauer synthesis

Cauer synthesis is an alternative synthesis to Foster synthesis and the conditions that a PRF must meet are exactly the same as Foster synthesis. Like Foster synthesis, there are two forms of Cauer synthesis, and both can be extended to RC and RL networks.

Cauer I form

Example of a Cauer I form realisation

The Cauer I form expands into a continued fraction. Using the same example as used for the Foster I form,

or, in more compact notation,

The terms of this expansion can be directly implemented as the component values of a ladder network as shown in the figure.[61] The given PRF may have a denominator that has a greater degree than the numerator. In such cases, the multiplicative inverse of the function is expanded instead. That is, if the function represents , then is expanded instead and vice versa.[62]

Cauer II form

Example of a Cauer II form realisation

Cauer II form expands in exactly the same way as Cauer I form except that lowest degree term is extracted first in the continued fraction expansion rather than the highest degree term as is done in Cauer I form.

mutual inductances.[65]

The essential reason that the example cannot be realised as a Cauer II form is that this form has a high-pass topology. The first element extracted in the continued fraction is a series capacitor. This makes it impossible for the zero of at the origin to be realised. The Cauer I form, on the other hand, has a low-pass topology and naturally has a zero at the origin.[66] However, the of this function can be realised as a Cauer II form since the first element extracted is a shunt inductor. This gives a pole at the origin for , but that translates to the necessary zero at the origin for . The continued fraction expansion is,

and the realised network is shown in the figure.

Brune synthesis

The Brune synthesis can synthesise any arbitrary PRF, so in general will result in a 3-element-kind (i.e. RLC) network. The poles and zeroes can lie anywhere in the left-hand half of the complex plane.

cascade of Brune sections.[69]

Removal of critical frequencies on the imaginary axis

Poles and zeroes on the axis represent L and C elements that can be extracted from the PRF. Specifically,

  • a pole at the origin represents a series capacitor
  • a pole at infinity represents a series inductance
  • a zero at the origin represents a shunt inductor
  • a zero at infinity represents a shunt capacitor
  • a pair of poles at represents a parallel LC circuit of resonant frequency in series
  • a pair of zeroes at represents a series LC circuit of resonant frequency in shunt

After these extractions, the remainder PRF has no critical frequencies on the imaginary axis and is known as a minimum reactance, minimum susceptance function. Brune synthesis proper begins with such a function.[70]

Broad outline of method

The essence of the Brune method is to create a conjugate pair of zeroes on the axis by extracting the real and imaginary parts of the function at that frequency, and then extract the pair of zeroes as a resonant circuit. This is the first Brune section of the synthesised network. The resulting remainder is another minimum reactance function that is two degrees lower. The cycle is then repeated, each cycle producing one more Brune section of the final network until just a constant value (a resistance) remains.[71] The Brune synthesis is canonical, that is, the number of elements in the final synthesised network is equal to the number of arbitrary coefficients in the impedance function. The number of elements in the synthesised circuit cannot therefore be reduced any further.[72]

Removal of minimum resistance

Extracting minimum resistance

A minimum reactance function will have a minimum real part, , at some frequency . This resistance can be extracted from the function leaving a remainder of another PRF called a minimum positive-real function, or just minimum function.[73] For example, the minimum reactance function

has and . The minimum function, , is therefore,

Removal of a negative inductance or capacitance

Extraction of negative inductance

Since has no real part, we can write,

For the example function,

In this case, is negative, and we interpret it as the reactance of a negative-valued inductor, . Thus,

and

after substituting in the values of and . This inductance is then extracted from , leaving another PRF, ,

[74]

The reason for extracting a negative value is because is a PRF, which it would not be if were positive. This guarantees that will also be PRF (because the sum of two PRFs is also PRF).[75] For cases where is a positive value, the admittance function is used instead and a negative capacitance is extracted.[76] How these negative values are implemented is explained in a later section.

Removal of a conjugate pair of zeroes

Extraction of a pair of zeroes

Both the real and imaginary parts of have been removed in previous steps. This leaves a pair of zeroes in at as shown by factorising the example function;[77]

Since such a pair of zeroes represents a shunt resonant circuit, we extract it as a pair of poles from the admittance function,

The rightmost term is the extracted resonant circuit with and .[78] The network synthesised so far is shown in the figure.

Removal of a pole at infinity

Extraction of a pole at infinity

must have a pole at infinity, since one was created there by the extraction of a negative inductance. This pole can now be extracted as a positive inductance.[79]

Thus as shown in the figure.

Replacing negative inductance with a transformer

Eliminate negative inductance by using a transformer

The negative inductance cannot be implemented directly with passive components. However, the "tee" of inductors can be converted into mutually coupled inductors which absorbs the negative inductance.[80] With a coupling coefficient of unity (tightly coupled) the mutual inductance, , in the example case is 2.0.

Rinse and repeat

In general, will be another minimum reactance function and the Brune cycle is then repeated to extract another Brune section[81] In the example case, the original PRF was of degree 2, so after reducing it by two degrees, only a constant term is left which, trivially, synthesises as a resistance.

Positive X

In step two of the cycle it was mentioned that a negative element value must be extracted in order to guarantee a PRF remainder. If is positive, the element extracted must be a shunt capacitor instead of a series inductor if the element is to be negative. It is extracted from the admittance instead of the impedance . The circuit topology arrived at in step four of the cycle is a

Π (pi)
of capacitors plus an inductor instead of a tee of inductors plus a capacitor. It can be shown that this Π of capacitors plus inductor is an equivalent circuit of the tee of inductors plus capacitor. Thus, it is permissible to extract a positive inductance and then proceed as though were PRF, even though it is not. The correct result will still be arrived at and the remainder function will be PRF so can be fed into the next cycle.[82]

Bott-Duffin synthesis

Example of Bott-Duffin synthesis steps 1 and 2
Example of Bott-Duffin synthesis steps 3 and 4

The Bott-Duffin synthesis begins as with the Brune synthesis by removing all poles and zeroes on the axis. Then Richards' theorem is invoked, which states for,

if is a PRF then is a PRF for all real, positive values of .[83]

Making the subject of the expression results in,[84]

An example of one cycle of Bott-Duffin synthesis is shown in the figures. The four terms in this expression are, respectively, a PRF ( in the diagram), an inductance, , in parallel with it, another PRF ( in the diagram), and a capacitance, , in parallel with it. A pair of critical frequencies on the axis is then extracted from each of the two new PRFs (details not given here) each realised as a resonant circuit. The two residual PRFs ( and in the diagram) are each two degrees lower than .[85] The same procedure is then repeatedly applied to the new PRFs generated until just a single element remains.[86] Since the number of PRFs generated doubles with each cycle, the number of elements synthesised will grow exponentially. Although the Bott-Duffin method avoids the use of transformers and can be applied to any expression capable of realisation as a passive network, it has limited practical use due to the high component count required.[87]

Bayard synthesis

Bayard synthesis is a

state-space synthesis method based on the Gauss factorisation procedure. This method returns a synthesis using the minimum number of resistors and contains no gyrators. However, the method is non-canonical and will, in general, return a non-minimal number of reactance elements.[88]

Darlington synthesis

Darlington synthesis starts from a different perspective to the techniques discussed so far, which all start from a prescribed rational function and realise it as a

one-port impedance. Darlington synthesis starts with a prescribed rational function that is the desired transfer function of a two-port network. Darlington showed that any PRF can be realised as a two-port network using only L and C elements with a single resistor terminating the output port.[89] The Darlington and related methods are called the insertion loss method.[90] The method can be extended to multi-port networks with each port terminated with a single resistor.[91]

The Darlington method, in general, will require transformers or coupled inductors. However, most common filter types can be constructed by the Darlington method without these undesirable features.[92]

Active and digital realisations

An example of a Sallen-Key cell low-pass filter

If the requirement to use only passive elements is lifted, then the realisation can be greatly simplified. Amplifiers can be used to

biquadratic filter.[95] Like the Darlington approach, Butterworth and Sallen-Key start with a prescribed transfer function rather than an impedance. A major practical advantage of active implementation is that it can avoid the use of wound components (transformers and inductors) altogether.[96] These are undesirable for manufacturing reasons.[97] Another feature of active designs is that they are not limited to PRFs.[98]

Digital realisations, like active circuits, are not limited to PRFs and can implement any rational function simply by programming it in. However, the function may not be stable. That is, it may lead to oscillation. PRFs are guaranteed to be stable, but other functions may not be. The stability of a rational function can be determined by examining the poles and zeroes of the function and applying the Nyquist stability criterion.[99]

References

  1. ^ Aatre, p. 259
  2. ^ E. Cauer et al., p. 4
  3. ^ Storer, p.3
  4. ^ Robertson et al., p. 107
  5. ^ Robertson et al., p. 108
  6. ^ Paarman, p. 15
  7. ^ Robertson et al., p. 107
  8. ^ Shenoi, pp. 4, 18
  9. ^ Paarman, pp. 15–16
  10. ^ Bakshi & Chitode, p. 6-1
  11. ^ Anderson & Vongpanitlerd, p. 509
  12. ^ Wanhammar, p. 10
  13. ^ Bakshi & Chitode, p. 6-1
  14. ^ E. Cauer et al., p.4
  15. ^ Hughes et al., p. 288
  16. ^ Hughes et al., p. 288
  17. ^ Bakshi & Chitode, p. 6-1
  18. ^ E. Cauer et al., pp. 3-4
  19. ^ Kalman, p. 4
  20. ^ E. Cauer et al., pp. 6-7
  21. ^ Kalman, p. 6
  22. ^ Kalman, p. 4
  23. ^ Hubbard, p. 3
  24. ^ Hubbard, p. 3
  25. ^ Wing, p. 122
  26. ^ Kalman, p. 7
  27. ^ Chen & Smith, p. 38
  28. ^ Chen & Smith, pp. 38, 50
  29. ^
    • Wing, p. 128
    • Kalman, p. 10
  30. ^ E. Cauer et al., p. 7
  31. ^ Wing, p. 128
  32. ^ Belevitch, p. 854
  33. ^ Lee, pp. 756-757
  34. ^ Kalman, p. 10
  35. ^ E. Cauer et al., p. 5
  36. ^ Swanson, p. 58
  37. ^ McCarthy, p. 51
  38. ^ Swanson, p. 58
  39. ^ Darlington, p. 7
  40. ^ Chen & Hu, p. 8
  41. ^ Chen & Smith, p. 35
  42. ^ Hughes et al., p. 286
  43. ^ Hughes et al., pp. 287–288
  44. ^ Awang, p. 227
  45. ^ Matthaei et al., pp. 3-5
  46. ^ Matthaei et al., p. 681
  47. ^ Matthaei et al., p. 5
  48. ^ Matthaei et al., pp. 121
  49. ^ Wanhammar, p. 58
  50. ^ Chen & Smith, p. 35
  51. ^ Chen & Hu, p. 8
  52. ^ Bakshi & Bakshi, p. 3-13–3-14
  53. ^ Wing, p. 115
  54. ^
    • Bakshi & Bakshi, p. 3-23, 3-30–3-31, 3-37
    • Wing, p. 115
  55. ^ Bakshi & Bakshi, pp. 3-23–3-24
  56. ^ Wing, pp. 115-116
  57. ^ Wing, pp. 115-116
  58. ^ Bakshi & Bakshi, pp. 3-31, 3-33
  59. ^ Bakshi & Bakshi, pp. 3-30–3-31, 3-37
  60. ^ Ghosh & Chakroborty, p. 771
  61. ^ Bakshi & Bakshi, pp. 3-28–3-29
  62. ^ Bakshi & Bakshi, p. 3-21
  63. ^ Bakshi & Bakshi, p. 3-29
  64. ^ Bakshi & Bakshi, p. 3-21
  65. ^ Houpis & Lubelfeld, p. 183
  66. ^ Bakshi & Bakshi, p. 3-30
  67. ^ Wing, p. 115
  68. ^ Hughes et al., p. 283
  69. ^ Carlin & Civalleri, p. 254
  70. ^
    • Wing, pp. 115-116
    • Ghosh & Chakroborty, pp. 792-793
  71. ^ Wing, pp. 117-118
  72. ^ Wing, p. 122
  73. ^
    • Wing, p. 116
    • Ghosh & Chakroborty, p. 793
  74. ^ Wing, pp. 116-117
  75. ^ Wing, p. 117
  76. ^ Wing, p. 119
  77. ^ Wing, p. 117
  78. ^ Wing, p. 117
  79. ^ Wing, p. 117
  80. ^ Wing, pp. 117-118
  81. ^ Wing, pp. 117, 120
  82. ^ Wing, p. 121
  83. ^ Wing, p. 122
  84. ^ Hughes et al., p. 284
  85. ^ Hughes et al., pp. 284–285
  86. ^
    • Wing, pp. 122–126
    • Youla, pp. 65–66
  87. ^ Kalman, p. 7
  88. ^ Anderson & Vongpanitlerd, p. 427
  89. ^ E. Cauer et al., p. 7
  90. ^ Sisodia, p. 5.13
  91. ^ Belevitch, p. 853
  92. ^ Wing, p. 164
  93. ^ Belevitch, p. 852
  94. ^ Belevitch, p. 850
  95. ^ Glisson, p. 727
  96. ^ Vaisband et al., p. 280
  97. ^ Comer & Comer, p. 435
  98. ^ Wing, p. 91
  99. ^ Chao & Athans, p. 524

Bibliography

Sources

Primary documents