Abstract

A new method of the symbolic circuit tolerance analysis has been proposed. The approach is based on Middlebrook's extra element theorem and the generalized parameter extraction method. It does not need the matrix network description or presentation as algebraic sets or topological graphs. The proposed techniques have been realized in the computer program (Toleralize). To verify the theoretical analysis, computer simulation results are included.

1. Introduction

The transfer function of the linear time-invariant network is generally presented as a rational fraction. But in many cases, especially in network synthesis, it is more appropriate to write the transfer function as follows [1]:𝑆=𝑆0(1+𝛾),(1) where 𝑆0 is a nominal transfer coefficient, and Ξ³ is a multiplicative relative error. Both the variables 𝑆0 and Ξ³ can be presented as a rational fraction.

The transfer function in the form (1) can be obtained by the usage of the Middlebrook’s extra element theorem [1–3] 𝑆=𝑆0ξ‚΅1+𝑍𝑛1/πœ’1+𝑍𝑑1ξ‚Ά,/πœ’π‘†=π‘†ξ…ž0ξ‚΅1+πœ’/𝑍𝑛11+πœ’/𝑍𝑑1ξ‚Ά,(2) where πœ’ is a parameter of an arbitrary circuit element (impedance, conductance, and controlled source), 𝑆0 and π‘†ξ…ž0 are nominal transfer functions of the initial circuit in which πœ’β†’βˆž or πœ’β†’0 correspondingly, 𝑍𝑛1 is a driving-point impedance of the network from the port where πœ’βˆ’element is connected, with the input port open, and 𝑍𝑑1 is the driving-point impedance of the network at the port where πœ’-element is connected, with the input port short.

In case of πœ’β†’βˆž, the selected element must be omitted from the circuit if πœ’ is an impedance and replaced by a short or by a nullor [4–6] if πœ’ is a conductance or a controlled source correspondingly. In case of πœ’β†’0, the selected element must be replaced by a short if πœ’ is an impedance, and deleted if πœ’ is a conductance or a controlled source.

So, the usage of formulae (2) required computation of three different circuit functions 𝑍𝑛1, 𝑍𝑑1 and 𝑆0 or π‘†ξ…ž0. This requirement makes problems for the symbolic circuit analysis of complicated network.

In Section 2 of this paper, the technique for improvement of the estimation efficiency of the extra element theorem is proposed. The obtained results provide the development of new symbolic circuit tolerance analysis method presented in Section 3. A practical example is discussed in Section 4. The conclusion is finally drawn in Section 5.

2. Extra Element Theorem and Generalized Parameter Extraction Method

In accordance with [7], the transfer function of electronic network can be expressed as 𝑁𝑆=𝐷,(3) where N is the determinant of the network, in which the independent source and arbitrary response are replaced by nullor, and D is the determinant of the network, in which the independent excitation and the arbitrary response are zero.

The calculation procedure of the network determinants is based on the recursive usage of the parameter extraction formula which generalizes the known Feussner's equations [8, 9]Ξ”=πœ’Ξ”(πœ’βŸΆβˆž)+Ξ”(πœ’=0),(4) where Ο‡ is a parameter of arbitrary circuit element, Ξ”(πœ’β†’βˆž) and Ξ”(πœ’β†’0) correspond to determinants of the circuit matrix in which the parameter of extracted elements πœ’β†’βˆž or πœ’β†’0, respectively.

The described techniques develop the singular elements approach discussed by Braun [7], Parten and Seacat [10], and Hashemian [11] and have been realized in the generalized parameter extraction method, an effective tool for network symbolic analysis [12–14] and synthesis [15]. The method can be used for the analysis of the circuit with all types of controlled sources. It requires neither matrix, nor ordinary graph description of the circuit. The main advantage of the method is that there are no cancellations among the generated terms.

The generalized parameter extraction method has been realized in the program CIRSYM as a part of the software tool SCAD (http://intersyn.narod.ru/scad.htm).

So, it is possible to present the formulae (2) with the help of network determinants as 𝑁𝑆=πœ’π·πœ’ξ‚΅πœ’+π‘πœ’/π‘πœ’πœ’+π·πœ’/π·πœ’ξ‚Ά,𝑁(5)𝑆=πœ’π·πœ’ξƒ©ξ€·π‘1+πœ’πœ’/π‘πœ’ξ€Έξ€·π·1+πœ’πœ’/π·πœ’ξ€Έξƒͺ,(6) where π‘πœ’ and π·πœ’ are determinants of networks, in which the extracted parameter πœ’β†’βˆž; NΟ‡ and DΟ‡ are determinants of networks, in which the extracted parameter πœ’β†’0.

As evident from the expressions (5) and (6), we have two pairs of identical network determinants. So, for estimation of the transfer function in the form (1), we need to obtain only four network determinants: π‘πœ’, π·πœ’, π‘πœ’, and π·πœ’. It is simpler than a separate computation of circuit functions from (2): 𝑆0=π‘πœ’/π·πœ’ or π‘†ξ…ž0=π‘πœ’/π·πœ’, 𝑍𝑛1=π‘πœ’/π‘πœ’, and 𝑍𝑑1=π·πœ’/π·πœ’. So, the proposed technique provides more than 30% increase of the computation efficiency.

In the case of a multiple extraction, that is, a simultaneous extraction of several parameters πœ’1,πœ’2,…,πœ’π‘›, it is possible to modify the formulae (5) and (6) as follows:𝑆=𝑆0πœ’1ξ‚€π‘πœ’1πœ’2,…,𝑛/π‘πœ’1,2,…,𝑛+β‹―+πœ’π‘›ξ‚€π‘πœ’π‘›πœ’1,2…,π‘›βˆ’1/π‘πœ’1,2,…,𝑛+πœ’1πœ’2ξ‚€π‘πœ’1,2πœ’3,…,𝑛/π‘πœ’1,2,…,𝑛+β‹―+πœ’1πœ’2β‹―πœ’π‘›+π‘πœ’(1,2,…,𝑛)πœ’(1,2,…,𝑛)/π‘πœ’1,2,…,π‘›πœ’1ξ‚€π·πœ’1πœ’2,…,𝑛/π·πœ’1,2,…,𝑛+β‹―+πœ’π‘›ξ‚€π·πœ’π‘›πœ’1,2,…,π‘›βˆ’1/π·πœ’1,2,…,𝑛+πœ’1πœ’2ξ‚€π·πœ’1,2πœ’2,…,𝑛/π·πœ’1,2,…,𝑛+β‹―+πœ’1πœ’2β‹―πœ’π‘›+π·πœ’(1,2,…,𝑛)πœ’(1,2,…,𝑛)/π·πœ’1,2,…,𝑛,𝑆=𝑆′01+πœ’1ξ‚€π‘πœ’1πœ’2,…,𝑛/π‘πœ’1,2,…,𝑛+β‹―+πœ’π‘›ξ‚€π‘πœ’π‘›πœ’1,2,…,π‘›βˆ’1/π‘πœ’1,2,…,𝑛+πœ’1πœ’2ξ‚€π‘πœ’1,2πœ’3,…,𝑛/π‘πœ’1,2,…,𝑛+β‹―+πœ’1πœ’2β‹―πœ’π‘›ξ‚€π‘πœ’(1,2,…,𝑛)πœ’(1,2,…,𝑛)/π‘πœ’1,2,…,𝑛1+πœ’1ξ‚€π·πœ’1πœ’2,…,𝑛/π·πœ’1,2,…,𝑛+β‹―+πœ’π‘›ξ‚€π·πœ’π‘›πœ’1,2,…,π‘›βˆ’1/π·πœ’1,2,…,𝑛+πœ’1πœ’2ξ‚€π·πœ’1,2πœ’2,…,𝑛/π·πœ’1,2,…,𝑛+β‹―+πœ’1πœ’2β‹―πœ’π‘›ξ‚€π·πœ’(1,2,…,𝑛)πœ’(1,2,…,𝑛)/π·πœ’1,2,…,𝑛,(7)where 𝑆0=π‘πœ’1,2,…,𝑛/π·πœ’1,2,…,𝑛 and π‘†ξ…ž0=π‘πœ’1,2,…,𝑛/π·πœ’1,2,…,𝑛.

3. Symbolic Circuit Tolerance Analysis

The network multiplicative relative error from (1) has a direct connection with the parameters spread of the circuit elements. In accordance with [16], the influence of the elements tolerances to the network characteristics can be simulated by including additional elements into the network. If we consider the additional elements as β€œextra” elements, it is possible to use the extra element theorem for the network tolerance analysis.

The tolerance 𝛿(πœ’π‘–) of the two-port circuit element πœ’i can be simulated by series resistances or by shunt admittances as shown in Figures 1(a) and 1(b) correspondingly. The circuit simulation of the controlled sources tolerances is presented in Figures 1(c)–1(f).

The symbolic expression of the network multiplicative relative error corresponding to the notations of circuit functions (1) and (4) can be defined as𝑁𝛾=𝛿(πœ’)𝛿(πœ’)𝐷𝛿(πœ’)βˆ’π‘π›Ώ(πœ’)𝐷𝛿(πœ’)𝛿(πœ’)𝑁𝛿(πœ’)𝐷𝛿(πœ’)+𝑁𝛿(πœ’)𝐷𝛿(πœ’).(8)

The fractional formulae for obtaining the tolerance of the arbitrary circuit elements Ξ΄(πœ’i) can be expressed in the same way𝛿(πœ’)=𝛾𝑁𝛿(πœ’)𝐷𝛿(πœ’)𝑁𝛿(πœ’)𝐷𝛿(πœ’)βˆ’π‘π›Ώ(πœ’)𝐷𝛿(πœ’)(𝛾+1).(9)

It is possible to generalize the expression (8) in case when multiple tolerance-simulated elements are added into a circuit simultaneously,𝛾(1,2,…,𝑛)=ξ‚Έπ›Ώξ€·πœ’1ξ€Έξ‚΅π‘π›Ώξ€·πœ’1ξ€Έπ›Ώξ€·πœ’2,…,π‘›ξ€Έπ·π›Ώξ€·πœ’1,2,…,π‘›ξ€Έβˆ’π‘π›Ώξ€·πœ’1,2,…,π‘›ξ€Έπ·π›Ώξ€·πœ’1ξ€Έπ›Ώξ€·πœ’2,…,π‘›ξ€Έξ‚Άξ€·πœ’+β‹―+π›Ώπ‘›ξ€Έξ‚΅π‘π›Ώξ€·πœ’π‘›ξ€Έπ›Ώξ€·πœ’1,2,…,π‘›βˆ’1ξ€Έπ·π›Ώξ€·πœ’1,2,…,π‘›βˆ’1ξ€Έβˆ’π‘π›Ώξ€·πœ’1,2,…,π‘›ξ€Έπ·π›Ώξ€·πœ’1ξ€Έπ›Ώξ€·πœ’1,2,…,π‘›βˆ’1ξ€Έξ‚Άξ€·πœ’+𝛿1πœ’2ξ€Έξ‚΅π‘π›Ώξ€·πœ’1,2ξ€Έπ›Ώξ€·πœ’3,…,π‘›ξ€Έπ·π›Ώξ€·πœ’1,2,…,π‘›ξ€Έβˆ’π‘π›Ώξ€·πœ’1,2,…,π‘›ξ€Έπ·π›Ώξ€·πœ’1,2ξ€Έπ›Ώξ€·πœ’3,…,π‘›ξ€Έξ‚Άξ€·πœ’+β‹―+𝛿1πœ’2β‹―πœ’π‘›ξ€Έξ‚΅π‘π›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έπ›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έπ·π›Ώξ€·πœ’1,2,…,π‘›ξ€Έβˆ’π‘π›Ώξ€·πœ’1,2,…,π‘›ξ€Έπ·π›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έπ›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έ/ξ‚Έπ‘ξ‚Άξ‚Ήπ›Ώξ€·πœ’1,2,…,π‘›ξ€Έξ‚΅π›Ώξ€·πœ’1ξ€Έπ·π›Ώξ€·πœ’1ξ€Έπ›Ώξ€·πœ’2,…,π‘›ξ€Έξ€·πœ’+β‹―+π›Ώπ‘›ξ€Έπ·π›Ώξ€·πœ’π‘›ξ€Έπ›Ώξ€·πœ’1,2,…,π‘›βˆ’1ξ€Έξ‚Άξ€·πœ’+𝛿1πœ’2ξ€Έπ·π›Ώξ€·πœ’1,2ξ€Έπ›Ώξ€·πœ’3,…,π‘›ξ€Έξ€·πœ’+β‹―+𝛿1πœ’2β‹―πœ’π‘›ξ€Έπ·π›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έπ›Ώξ€·πœ’(1,2,…,𝑛)ξ€Έ+π·π›Ώξ€·πœ’1,2,…,𝑛(10)

Note that the value of network determinants π‘πœ’ and π·πœ’ in formulae (8) and (9) as well as determinants π‘πœ’1,2,…,𝑛и π·πœ’1,2,…,𝑛 in (10) is invariable and can be found in the stage of nominal transfer function estimation.

The proposed techniques of tolerance analysis have been realized in the computer program Toleralize as a part of the software tool SCAD. Toleralize provides the estimation of symbolic expressions of a network transfer function in forms (5) or (6), a network multiplicative relative error in forms (8) or (10), and an element tolerance in form (9), as well as their numerical values.

4. Example

As an example, let us consider the procedure of the tolerance analysis of the transistor amplifier shown in Figure 2 [16]. The elements parameters are presented in Table 1 (first column).

It is possible to use the original formulae (2) for the tolerance analysis, but the technique on the base of network determinants approach considered in Section 3 makes the calculation process much simpler and faster.

First, we compute the nominal transfer function of the network by means of a generalized parameter extraction method𝑁𝛿(πœ’π‘–)=𝑅6𝑅3+𝑅4𝑅2𝐡1+ξ€·+1ξ€Έξ€Έβˆ’π΅2𝑅3ξ€Έξ€·βˆ’π΅1𝑅1+𝑅2𝐷=101020200000000,(11)𝛿(πœ’π‘–)=𝑅3+𝑅4ξ€Έξ€·βˆ’π΅1𝑅1𝑅2+𝑅1𝑅2+𝑅5ξ€Έ+𝑅2𝑅5ξ€Έ+𝑅3+𝑅4𝑅6ξ€Έξ€·βˆ’π΅1𝑅2+𝑅2+𝑅5ξ€Έ+𝑅3𝑅6𝐡2𝑅2𝐡1ξ€Έ=1244400000000.(12)

So, the nominal transfer function of the amplifier is 𝑆0=81,18.

Then, we use formula (8) for the estimation of network multiplicative relative errors provided that the elements parameters spread is 𝛿(πœ’π‘–)=Β±10%. Note that we already know the value of network determinants 𝑁𝛿(πœ’π‘–) (11) and 𝐷𝛿(πœ’π‘–) (12), so we need to obtain only two network determinants 𝑁𝛿(πœ’π‘–) and 𝐷𝛿(πœ’π‘–) for each πœ’π‘– circuit elements by means of a tolerances simulation approach. The results of the computation of the determinants symbolic expressions and their numerical values are presented in Table 1 (second and third columns, correspondingly).

Now, we estimate the elements tolerances provided that the multiplicative relative error of the amplifier transfer function is 𝛾=Β±15%. The numeric values of network determinants included in expression (9) are already obtained, so the computation process becomes really simple. The calculated results of the parameters spread of circuit elements are presented in Table 1 (fourth column).

The results obtained have been verificated by the computer program Toleralize.

5. Conclusions

In this paper, a network determinant approach implemented for symbolic circuit analysis by means of the extra element theorem has been considered. The formulae for analytical evaluation of a network multiplicative relative error and the tolerance of the arbitrary circuit elements are presented. The proposed techniques provide more than 30% increase of estimation efficiency as compared to the original Middlebrook’s method. The techniques proposed have been realized in the computer program Toleralize.