A Cancellation-Free Symbolic Sensitivity Technique Based on Network Determinant Expansion
The generalization of Bode’s sensitivity analysis technique for all types of the transfer functions and circuit elements is presented in the paper. The proposed formulae for first- and second-order symbolic sensitivity calculation provide the compact size of obtained expression and have the advantages of cancellation-free sum-of-product terms and matrix-free computation. This is achieved by means of the concept of high order summative cofactors and the generalized parameter extraction method. The proposed technique is implemented in symbolic circuit analysis program Cirsym. Illustrative example on symbolic sensitivity circuit analysis and comparison of the presented technique with the transimpedance method and the method based on the modified Coates flow-graph are given.
G1The influence of network components on transfer function is expressed by sensitivities . Sensitivity analysis is an important part of analog circuit design process. There are many methods of numerical calculation of circuit sensitivities, but only symbolic analysis of circuit sensitivity provides the way to get analytical sensitivity function of all interested circuit parameters. This is certainly advantageous in such applications as circuit parameters optimization and Monte Carlo simulation.
First researches on symbolic approach to sensitivity analysis of analog circuit were presented in [1–5], but those heuristic algorithms were not efficient enough and had limitations by the circuits scale.
Nowadays several more effective symbolic sensitivity analysis techniques have been developed. Some of them are based on sequence of expression computation such as two-port transimpedance method . But the generated expressions are not compact enough, do not fully exploit sharing of subexpressions, cannot be easily manipulated, and complicate the circuit insight.
The others methods are based on implementation of graph theory for symbolic circuit analysis [7–10]. In particular, the work in  shows the technique based on the modified Coates flow-graph for performing sensitivity analysis. However, the proposed concept requires the nodal admittance matrix expansion, which may consist of the equal summands with opposite sign. It leads to appearing of the cancellation sum-of-product terms. In addition, the usage of graph-based methods does not provide the optimal size of obtaining expression.
Another widespread approach is implementation of the modified nodal matrix [10, 11]. But obtained results are not compact and include cancellations so the expressions often present only in semisymbolic form.
In this paper we propose a cancellation-free symbolic sensitivity technique for computation of compact expression. The first- and second-order symbolic sensitivities formulas for all types of the transfer functions and circuit elements are derived by means of the concept of higher order summative cofactors (HOSC) [12–15] and generalized parameter extraction method (GPEM) [16–21] to avoid drawbacks of the previous methods mentioned above.
The rest of this paper is organized as follows. In Section 2 we explain the basic idea of sensitivity analysis by means of GPEM and derive the formulae for symbolic sensitivity of different transfer functions with respect to controlled sources (CS) and two ports. The process of obtaining symbolic expressions can be easily automated. The example of the AC sensitivity analysis is presented in Section 3. In Section 4 the comparison of the presented method with the transimpedance method and the method based on the modified Coates flow-graph is performed. Section 5 concludes this paper.
2. Theoretical Basis of the Proposed Method
2.1. Symbolic Approach for Sensitivity Analysis
The sensitivity computation of the transfer function of with respect to change of the parameter requires finding the corresponding derivative:
The study of network response by means of formula 1 is acceptable only when changes in its parameter are small. When some element’s parameters change significantly, this approximated result may be unsatisfied. Analytical sensitivity analysis would be much easier by using a symbolic approach than by a numerical approach where computation cost is a key concern.
Calculation of the transfer function derivatives leads to appearance of multiple algebraic cofactors with different signs and thus they may contain many cancellations. The well-known technique proposed by Bode  presented the sensitivity function in form of ratio of determinant and cofactors of nodal admittance or loop impedance matrix multiplication. For example, the sensitivity of of the circuit in Figure 1 is expressed as follows: where is a transfer impedance of current controlled voltage source (CCVS) or transfer admittance of voltage controlled current source (VCCS), is a determinant of loop impedance matrix or for nodal admittance matrix, and , , , and are cofactors.
Expression 3 derived from 1 by means of Jacobi identity in form of where is a second-order cofactor. The rows and columns with numbers 1 and 3, 5, and 7 correspond to input and output variables (controlled and controlling variables of controlled source).
Sensitivity analysis by means of Bode’s technique leads to appearance of cancellations because of multiple entries of elements parameters into the matrix. The sensitivity computation formulae with numbered ports of multiport circuit model presented in  are devoid of this drawback. But the tedious combinatorial search of the transfer function from one port to another is needed to use them.
We propose in this paper the generalization of Hoang formula and Bode’s technique for all types of transfer functions and circuit elements, based on implementation of the concept of higher order summative cofactors (HOSC) [12–15] instead of ordinary algebraic cofactors and computation of HOSC by means of GPEM [16–21].
2.2. Higher Order Summative Cofactors and Generalized Parameter Extraction Method
The higher order cofactor is a cofactor of a cofactor. The th order cofactor is represented by a symbol , where and are the numbers of deleted rows and columns, respectively. The summative cofactors are cofactors with at least one deletion in a form or , which means the following: add or subtract the row (column) to or from the row (column) . For example the summative cofactors of the same matrix can be expressed as follows : where , , and so on.
Let us express the voltage gain function for circuit model in Figure 2(a) as follows: where HOSC in the numerator is the determinant of the network, in which the independent source and arbitrary response are replaced by oriented nullor  and HOSC in the denominator 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 : where is a parameter of arbitrary linear circuit element and and correspond to the determinants of the circuit matrix in which the parameter of extracted elements or , respectively.
The parameter extraction formula 7 is applicable for determinant calculation of the network consisting of any kind of linear active circuit models including the pathological elements (nullors and pathological mirrors). In case of nullor extraction the special formula was proposed in  where is the determinant of the circuit matrix after the extraction procedure of the nullor number . Equation 8 means that the nullor extraction will change the sign of the initial determinant. The choice of the sign depends on the orientation of the nullator and norator. If these elements have got the same orientation with respect to the basic node that the sign will be positive, in the opposite case the sign will be negative.
The procedure of nullor extraction can be formalized by the following steps.(1)The choice of the supporting nodes: first of the supporting nodes should be connected to norator and the second to nullator. If there is a common node of a nullator and a norator it must be chosen as single supporting node. Note that supporting node may correspond to the ground node.(2)The terminals of non-extracted norators and current sources connected to supporting node are moved to the opposite node of extracted norator. In that case the non-extracted nullators and controlling voltages keep connections to supporting node. Then in the same way, the terminals of non-extracted nullators and controlling voltages connected to supported node are moved to the opposite node of extracted nullator. In that case non-extracted norators and current sources keep connections to supporting node.(3)A norator and a nullator of the extracted nullor are deleted from circuit. In case of two supporting nodes they must combine.(4)If extracted norator and nullator have got the same orientation with respect to the supporting node then the determinant sign will be positive. Otherwise the sign will be negative. In case of two supporting nodes the inverted rule is needed.
The generalized parameter extraction method does not have any limitation by circuit size. The GPEM has been realized by Filaretov in the program Cirsym (circuit simulator) as a part of the software tool SCAD. The program can be downloaded from http://intersyn.narod.ru/. The symbolic expressions obtained by means of Cirsym are comparable by the number of arithmetical operations with computation results by means of the professional mathematic programs as Maple and Matlab .
As can be seen from the comparison of symbolic expressions of notch filter transfer function  presented in “Symbolic circuit analysis: Library of benchmark circuits” (http://rodanski.net/ben/work/symbolic/index.htm) the result of program Cirsym named by Filaretov 2 is second ranked. Also this result was replicated in  as example of compact size expression.
2.3. Derivation of Symbolic Sensitivity Computation Formulae
Formula 10 is applicable for sensitivity computation of the arbitrary circuit elements. Just replace the by .
We obtain the transfer function sensitivity expression with respect to a VCCS parameter. It is possible to minimize the numerator by means of Jacobi identity. We need to use the HOSC instead of determinant :
As can be seen, is a numerator of transfer function from circuit input to controlling voltage of VCCS, is a numerator of transfer function from controlled current source to circuit output, and is a product of numerator and denominator of voltage gain.
The process of sensitivity computation with respect to other elements and/or for different transfer function types is the same. In Table 1 the generalized sensitivity formulae are proposed.
Even if output nodes of transfer impedance function in Figure 2(b) are connected in contrast to Figure 2(a), it is obvious that expression 13 is equal to 6. The formulae for transfer impedance (Figure 2(c)) and current gain (Figure 2(d)) are similar like expressions 6 and 13:
The denominator of any sensitivity formulae in Table 1 is a product of initial transfer function numerator and denominator and the numerator of the sensitivity formulae is a product of numerators of transfer functions from circuit input to controlled branch and from controlling source to circuit output.
The proposed formulae are easy to use for sensitivity computation of two ports. For example, if we change indexes 7 and 8 to 5 and 6 correspondingly, we will obtain the sensitivity function for admittance element:
The formula for transfer function sensitivity with respect to impedance parameter () can be derived in the same way. We just have to replace the number of index 8 by 6.
For sensitivity of the input impedance or input admittance functions derivation it is necessary to replace indexes 3 and 4 by 1 and 2 correspondingly.
2.4. The Special Topological Cases of Sensitivity Formulae Simplification
If the voltage source of CS is connected in parallel with a voltage output port or if the current source of CS is connected in series with a current output port, we can simplify the sensitivity formulae like that shown in Table 2. As it can be seen, the HOSC indexes 3 and 4 correspond to 7 and 8. Taking into account the identity we obtain the expressions consisting of just two HOSC.
In the same way, we can easily simplify the sensitivity formulae for other special cases if the input excitation transfers directly to the controlling branch of CS, if the controlled voltage source is connected to the input voltage port in series or the controlled current source is connected in parallel with input current source, and if the controlling voltage of CS has the same port as output circuit voltage or the controlling current of CS corresponds to output current port.
2.5. The Second-Order Sensitivity Computation
Let us consider the process of computation of second-order sensitivity function. The repeated differentiation of first-order formulae from Table 1 is needed in this case. So the second-order transfer impedance sensitivity of circuit in Figure 3 with respect to the VCCS and CCVS parameters can be expressed as follows:
We can see the product of squared determinant and the transfer function numerator in formula denominator and two terms in numerator. The first term here consists of the numerators of transfer functions, from circuit input to controlling voltage of VCVS, from voltage source of VCVS to controlling current of CCCS, and from current source of CCCS to circuit output. For the second term we take into account that succession of excitation transfer from input to output of circuit is different. In that case the signal transfers firstly to CCCS and then to VCVS.
For the purpose of obtaining the second-order sensitivity function in Hoang form we just need to differentiate and represent by means of HOSC formula 2. So for the circuit presented in Figure 3 the sensitivity expression for voltage transfer function with respect to parameters of CS is obtained as follows: where upper indexes and mean that those HOSC are derived with and correspondingly. The index means that this HOSC is derived with both parameters tending to infinity.
3. Symbolic Sensitivity Analysis Examples
As can be seen, the circuit topology makes it possible to simplify sensitivity computation with respect to VCVS parameter by means of formula from the first line of Table 2. Taking into account the nodes numeration of filter circuit in Figure 4 the sensitivity function can be expressed as follows:
The second summand of expression (Figure 7) is correspond to the HOSC , the numerator of formula 18. It is desirable to obtain the polynomial symbolic expression so the circuit expansion is performed by extracting capacitors admittances.
Note, that admittances extraction leads to appearance of two subcircuits with short circuit and deleting of the element correspondingly. The determinant of the first circuit in (Figure 8) is equal to unity because both admittances connected in short circuit. After expansion of element from second circuit in (Figure 8) we obtain the determinant of short-circuited admittance in the same way. In the third circuit with parallel admittances the determinant is equal to the elements sum. And the last circuit determinant is a product of and .
For computation of expression (Figure 9) the following procedures were done: deletion of capacities admittance connecting in parallel with nullator, parameter extraction of admittance connecting in series with nullator, followed by deletion of the admittance connecting in parallel with nullator. As a result we obtain norator-nullator pair connected in opposite direction and determinant of that kind of circuit is equal to “.”
The derivation process of sensitivity functions with respect to other parameters is similar. For example we derive the first- and second-order sensitivities and : where
Using Mathcad software symbolic sensitivity simulation results versus frequency are shown in Figure 5.
As it is seen, the symbolic sensitivity analysis technique based on the concept of high order summative cofactors and GPEM rules can be implemented in circuit design process and for circuit characterization, successfully.
4. Comparison Test
4.1. The Comparison Conditions
In this section we present the comparison of sensitivity expressions obtained by means of three different approaches: the two-port transimpedance method , the method based on the modified Coates flow-graph , and the proposed technique. The test circuit is shown in Figure 6 .
Let us suppose that voltage transfer function and the first- and second-order symbolic sensitivities with respect to parameters and are under consideration. We compare the formulae by following criteria in order to estimate the methods calculation efficiency: the number of arithmetical operations and occurrence of cancellations.
4.2. Symbolic Expressions Obtained by means of the GPEM-Based Technique
The symbolic expressions , , , and shown below are calculated by means of the symbolic circuit analysis program Cirsym. The sensitivity functions are presented here in two forms: Hoang form 10 and Bode form (Table 1).
The transfer function where
The sensitivity functions in Hoang form where where where
The sensitivity functions in Bode form where where where
4.3. Symbolic Expressions Obtained by means of Modified Coates Flow-Graph Concept
The Coates flow-graph is useful and often used in the network theory. The modification of Coates flow-graph provides the symbolic sensitivity analysis of nullor-based equivalent circuits . We will show here that modification of Coates flow-graph is corresponding to implementation of the nullor extraction formula 8.
Let us consider the flow-graph presented in Figure 8(a) (see publication ) for the passive part of the equivalent circuit of STAR network shown in Figure 6. The oriented norator connected between nodes 3 and 0. The oriented nullator connected between nodes 5 and 4.
First, we choose the node 3 as supporting and make the norator short-circuited:(1)the ends of all incoming edges ( and ) into node 3 are moved to node 0;(2)the self-loop transforms into the incoming edge into node 0;(3)the edges connected to node 0 are not shown on graph.
Second, we choose the node 5 as supporting and make the nullator short-circuited:(1)the end of incoming edge into node 5 is moved to node 4;(2)the self-loop transformed into the incoming edge into node 4.
Third, we combine the two supporting nodes 3 and 5:(1)the edge between nodes 3 and 5 transforms into the self-loop in the new node named by 3;(2)the parallel codirectional edges and are combined into the edge .
Consequently we obtained the modified Coates flow-graph shown in Figure 8(b) (see publication ). Note that in case of single supporting node the procedure of nodes combining is not needed.
The estimation process of the determinant sign of modified Coates flow-graph is not formalized in . It leads to complications of method usage and high error probability.
Implementation of the oriented nullor concept  provides the simple rule for estimation of the sign of Coates graph determinant: the sign is positive (negative), if norator and nullator have the same (different) orientations against the supporting nodes. In case of two supporting nodes the inverted rule is needed.
The symbolic formulae , , , and presented here were replicated from  with correction of typing error in subexpression and sign of determinants: where
Let us substitute the subexpressions , , and into the transfer function 36:
It is easy to verify that the expanded expressions of the numerator and denominator will consist of the 16 and 28 summands correspondingly. After simplification only 12 terms in numerator and 16 terms in denominator will remain. So, the transfer function 36 includes the 8 pairs of cancellation sum-of-product terms:
4.4. Symbolic Expressions Obtained by means of Two-Port Transimpedance Method
The sequences of expressions, derived by symbolic analysis program STAINS , are listed below. The calculation of second-order sensitivity is not support in the current version of the program. So we present here only the transfer function and two first-order symbolic sensitivities with respect to parameters and : where
After substitution of the subexpressions into 42 we obtain the transfer function expression as follows:
The expansion of the transfer function numerator and denominator leads to obtaining the expressions with 16 and 28 terms correspondingly. Just as formula 36 the transfer function 42 includes the 8 pairs of cancellation sum-of-product terms: where where
4.5. Efficiency Comparison
In Table 3 we present the comparison of symbolic analysis methods efficiency by the number of arithmetical operations. The multiplication and division are denoted as M/D and the addition and subtraction are denoted as A/S. The method based on modified Coates flow-graph is denoted as MCFG and the transimpedance method is denoted as TI.
As can be seen, the symbolic expressions obtained by transimpedance method consist of the largest arithmetical operations amount. The results of MCFG-computation are comparable to expressions presented in Bode form obtained by GPEM. The symbolic sensitivity functions presented in Hoang form are superior by number of arithmetical operations.
Both Bode and Hoang forms implement the four minor calculations. But the presentation of sensitivity function in Hoang form provides the optimal compact size of symbolic expressions. On the other hand the formulae in Bode form for special topological cases proposed in this paper provide the significant simplification of circuit analysis in many cases.
In this paper we have presented a method of computation the circuit functions sensitivities in symbolic form. The method requires neither matrix nor ordinary graph description of the circuit. The main advantages of the proposed technique are that it is cancellation-free and provided the compact size of obtaining sensitivity expression as distinct from other methods like the two-port transimpedance method and the method based on the modified Coates flow-graph. It is shown that presentation of sensitivity function in Hoang form provides the optimal compact size of symbolic expressions, and the formulae in Bode form provide the significant simplification of active circuit symbolic analysis. The process of calculation of first- and second-order and multiparameter symbolic sensitivity is automated by computer program Cirsym and allows obtaining all symbolic sensitivities simultaneously.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
H. W. Bode, Network Analysis and Feedback Amplifier Design, D. Van Nostrand, New York, NY, USA, 1945.
S. Hoang, “Direct topological method in network sensitivity analysis,” Archiwum Elektrotechniki, vol. 21, no. 4, pp. 767–784, 1975.View at: Google Scholar
H. Yang, M. Ranjan, W. Verhaegen, M. Ding, R. Vemuri, and G. Gielen, “Efficient symbolic sensitivity analysis of analog circuits using element coefficient diagrams,” in Proceedings of the Asia South-Pacific Design Automation Conference (ASPDAC '05), pp. 230–235, Yokohama, Japan, 2005.View at: Google Scholar
S. Rodriguez-Chavez, A. A. Palma-Rodriguez, E. Tlelo-Cuautle, and S. X.-D. Tan, “Graph-based symbolic and symbolic sensitivity analysis of analog integrated circuits,” in Analog/RF and Mixed-Signal Circuit Systematic Design, vol. 233 of Lecture Notes in Electrical Engineering, pp. 101–122, Springer, 2013.View at: Google Scholar
L. Mandache, M. Iordache, L. Dumitriu, I. Gabriela Sirbu, and D. Niculae, “Sensitivity analysis of analog circuits based on a modified nodal approach,” in Proceedings of the International Conference on Optimization of Electrical and Electronic Equipment (OPTIM '14), pp. 83–88, Bran, Romania, May 2014.View at: Google Scholar
P. M. Lin, Symbolic Network Analysis, Elsevier, Amsterdam, The Netherlands, 1991.
S. Lasota, “Parameter decision diagrams in the analysis and structural synthesis: part I, III,” in Proceedings of the 10th International Workshop on Symbolic and Numerical Methods, Modeling and Application to Circuit Design (SMACD '08), pp. 149–157, 180–187, Erfurt, Germany, 2008.View at: Google Scholar
S. Lasota, “Multilevel hierarchical always cancellation-free symbolic analysis method for large electric networks,” Elektronika: Konstrukcje, Technologie, Zastosowania, vol. 54, no. 2, pp. 51–57, 2013.View at: Google Scholar
V. V. Filaretov and A. S. Korotkov, “Generalized parameter extraction method in symbolic network analysis,” in Proceedings of the European Conference on Circuits Theory and Design (ECCTD '03), vol. 2, pp. 406–409, Kraków, Poland, September 2003.View at: Google Scholar
V. V. Filaretov and A. S. Korotkov, “Generalized parameter extraction method in case of multiple excitation,” in Proceedings of the 8th International Workshop on Symbolic Methods and Applications to Circuit Design (SMACD '04), pp. 8–11, Wrocław, Poland, September 2004.View at: Google Scholar
V. Filaretov, K. Gorshkov, and A. Mikheenko, “A circuit synthesis technique based on network determinant expansion,” in Proceedings of the International Conference on Synthesis, Modeling, Analysis and Simulation Methods and Applications to Circuit Design (SMACD '12), pp. 293–296, Seville Sevilla, Spain, September 2012.View at: Publisher Site | Google Scholar
V. V. Filaretov and K. S. Gorshkov, “Topological analysis of active networks containing pathological mirror elements,” in Proceedings of the IEEE 33rd International Scientific Conference Electronics and Nanotechnology (ELNANO '13), pp. 293–296, Kiev, Ukrain, April 2013.View at: Publisher Site | Google Scholar
R. Dmytryshyn and A. Kubaszek, “Sequence of expressions generation for the repetitive analysis acceleration,” in Proceedings of the International Workshop on Symbolic Methods, Modeling and Application in Circuit Design (SMACD ’98), pp. 154–159, Kaiserslautern, Germany, 1998.View at: Google Scholar
B. Rodanski, “Computational efficiency of symbolic sequential formulae,” in Proceedings of the 6th International Workshop on Symbolic Methods and Applications to Circuit Design (SMACD '00), pp. 45–50, Lisbon, Portugal, October 2000.View at: Google Scholar
J. Vlach and K. Singhal, Computer Methods for Circuit Analysis and Design, Van Nostrand Reinhold Company, New York, NY, USA, 1983.