# A Cancellation-Free Symbolic Sensitivity Technique Based on Network Determinant Expansion

Authors: Vladimir Filaretov; Konstantin Gorshkov; Sergey Kurganov
Journal: Advances in Electrical Engineering
 (2015)

---

## Abstract

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.


---

## Body

## 1. Introduction

G1The influence of network components on transfer function is expressed by sensitivities [1]. 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 [6]. 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 [7] 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 Section2 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 ofT(s) with respect to change of the parameter X requires finding the corresponding derivative:
(1)S(H(s),X)=dH(s)H(s)XdX.The study of network response by means of formula1 is acceptable only when changes in its parameter X 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.If we consider the transfer function asT(s)=N(s)/D(s) formula 1 will transform into expression proposed by Hoang [5]:
(2)S(T(s),X)=XdN(s)N(s)-dN(s)D(s).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 [1] presented the sensitivity function T(s) in form of ratio of determinant and cofactors of nodal admittance or loop impedance matrix multiplication. For example, the sensitivity of T(s)=V3(s)/V1(s) of the circuit in Figure 1 is expressed as follows:
(3)SXT=-XΔ15Δ73ΔΔ13,
where X is a transfer impedance Z of current controlled voltage source (CCVS) or transfer admittance of voltage controlled current source (VCCS), Δ is a determinant of loop impedance matrix or Δ=Δ11 for nodal admittance matrix, and Δ11, Δ15, Δ13, and Δ73 are cofactors.Figure 1
Test circuit.Expression3 derived from 1 by means of Jacobi identity in form of
(4)Δ13,75Δ=Δ13Δ75-Δ15Δ73,
where Δ13,75 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 [5] 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. Thenth order cofactor is represented by a symbol Δr1,k1,r2,k2,…,rn,kn, where r1,r2,…,rn and k1,k2,…,kn are the numbers of deleted rows and columns, respectively. The summative cofactors are cofactors with at least one deletion in a form (a+d) or (a-d), which means the following: add or subtract the row (column) a to or from the row (column) d. For example the summative cofactors of the same matrix can be expressed as follows [15]:
(5)Δ(a+b)(c+d)=Δac-Δad-Δbc+Δbd,
where Δac=Δ(a+0)(c+0), Δad=Δ(a+0)(d+0), and so on.In accordance with GPEM basis [16–21] we consider the sums of row (column) indexes as a result of nullator (norator) connection to the network nodes corresponding to the indexes.Let us express the voltage gain function for circuit model in Figure2(a) as follows:
(6)T(s)=Vout(s)Vin(s)=Δ(1+2)(3+4)Δ(1+2)(1+2),
where HOSC in the numerator is the determinant of the network, in which the independent source and arbitrary response are replaced by oriented nullor [16] and HOSC in the denominator is the determinant of the network, in which the independent excitation and the arbitrary response are zero.Figure 2
The circuit for sensitivity computation of transfer functions: voltage gain (a), transfer admittance (b), transfer impedance (c), and current gain (d).
(a)
(b)
(c)
(d)The calculation procedure of the network determinants is based on the recursive usage of the parameter extraction formula [16]:
(7)Δ=χΔ(χ⟶∞)+Δ(χ=0),
where χ is a parameter of arbitrary linear circuit element and Δ(χ→∞) and Δ(χ=0) correspond to the determinants of the circuit matrix in which the parameter of extracted elements χ→∞ or χ=0, respectively.The parameter extraction formula7 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 [16]
(8)Δ=±Δn,
where χ is the determinant of the circuit matrix after the extraction procedure of the nullor number n. 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.Formula8 can be implemented not only for the nullor extraction from circuit model, but also for Coates flow-graph modification as shown in Section 4.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 fromhttp://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 [22].As can be seen from the comparison of symbolic expressions of notch filter transfer function [23] 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 [24] as example of compact size expression.


### 2.3. Derivation of Symbolic Sensitivity Computation Formulae

Let us express the transfer function6 by means of parameter extraction formula as
(9)T(s)=YΔ(1+2)(3+4),(7+8)(5+6)+Δ(1+2)(3+4)0YΔ(1+2)(1+2),(7+8)(5+6)+Δ(1+2)(1+2)0,
and then expressed the Hoang formula 2 as follows:
(10)SYT=Y·Δ(1+2)(1+2)∞Δ(1+2)(1+2)-Δ∞Δ.Note, that upper indexes〈〈0〉〉 and 〈〈∞〉〉 in formulas 9 and 10 mean that those HOSC are derived with Y=0 and Y→∞ correspondingly. In case of controlled source the operation χ→∞ leads to substitution of CS by the nullor.Formula10 is applicable for sensitivity computation of the arbitrary circuit elements. Just replace the 〈〈Y〉〉 by 〈〈χ〉〉.Now let us consider the transformation of Bode formula1 by substitution of expression 9 and cancellation by Δ(1+2)(1+2):
(11)SYT=YΔ1+23+4,7+85+6Δ1+21+2Δ1+21+2Δ1+23+4-YΔ(1+2)(3+4)Δ(1+2)(1+2),(7+8)(5+6)Δ(1+2)(1+2)Δ(1+2)(3+4).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Δ(1+2)(1+2) instead of determinant Δ:
(12)SYT=-YΔ(1+2)(5+6)Δ(1+2)(1+2),(7+8)(3+4)Δ(1+2)(1+2)Δ(1+2)(3+4).As can be seen,Δ(1+2)(1+2) is a numerator of transfer function from circuit input to controlling voltage of VCCS, Δ(1+2)(1+2),(7+8)(3+4) is a numerator of transfer function from controlled current source to circuit output, and Δ(1+2)(1+2)Δ(7+8)(3+4) 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 Table1 the generalized sensitivity formulae are proposed.Table 1
Sensitivity symbolic calculation formulae.
Number
Sensitivity type
Formula
1
S
χ
T
-
χ
Δ
(
1
+
2
)
(
5
+
6
)
Δ
(
1
+
2
)
(
1
+
2
)
,
(
7
+
8
)
(
3
+
4
)
Δ
(
1
+
2
)
(
1
+
2
)
Δ
(
1
+
2
)
(
3
+
4
)
2
S
χ
Y
3
S
χ
B
-
χ
Δ
(
1
+
2
)
(
5
+
6
)
Δ
(
7
+
8
)
(
3
+
4
)
Δ
Δ
(
1
+
2
)
(
3
+
4
)
4
S
χ
ZNote that the formula in first line of Table1 is appropriate not only for voltage gain, but also for transfer impedance function. Let us show it for circuit in Figure 2(b):
(13)Ys=IoutsVins=Δ1+23+4Δ1+21+2.Even if output nodes of transfer impedance function in Figure2(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:
(14)Z=VoutIin=Δ(1+2)(3+4)Δ,B=IoutIin=Δ(1+2)(3+4)Δ.The denominator of any sensitivity formulae in Table1 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:(15)SYT=-YΔ(1+2)(5+6)Δ(1+2)(1+2),(5+6)(3+4)Δ(1+2)(1+2)Δ(1+2)(3+4).The formula for transfer function sensitivity with respect to impedance parameter (χ=Z) 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 Table2. As it can be seen, the HOSC indexes 3 and 4 correspond to 7 and 8. Taking into account the identity χΔ(1+2)(5+6)=Δ(1+2)(7+8) we obtain the expressions consisting of just two HOSC.Table 2
The simplified sensitivity formulae in case of controlled source connecting with circuit output.
Number
Sensitivity type
Formula
Circuit
Circuit
1
S
μ
T
S
Z
K
Δ
(
1
+
2
)
(
1
+
2
)
,
(
3
+
4
)
(
3
+
4
)
Δ
(
1
+
2
)
(
1
+
2
)
2
S
μ
Z
S
Z
Z
Δ
(
3
+
4
)
(
3
+
4
)
Δ
3
S
Y
Y
S
β
Y
Δ
(
1
+
2
)
(
1
+
2
)
0
Δ
(
1
+
2
)
(
1
+
2
)
4
S
Y
B
S
β
B
Δ
0
Δ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 Table1 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:
(16)SβYZ=βYΔ1+25+6Δ7+89+10,1+21+2hhihhh×Δ11+123+4,1+21+2+Δ1+29+10hhhihh×Δ11+125+6,1+21+2Δ7+83+4,1+21+2×Δ2Δ1+23+4-1.Figure 3
The test circuit for second-order transfer impedance sensitivity computation.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 formula2. 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:
(17)
S
β
Y
T
=
β
Y
Δ
1
+
2
3
+
4
∞
Δ
2
-
Δ
Y
→
∞
Δ
1
+
2
3
+
4
β
→
∞
Δ
hhhhhh
-
Δ
1
+
2
3
+
4
Y
→
∞
Δ
β
→
∞
Δ
-
Δ
(
1
+
2
)
(
3
+
4
)
Δ
∞
Δ
hhhhhh
+
2
Δ
1
+
2
3
+
4
Δ
β
→
∞
Δ
Y
→
∞
×
Δ
2
Δ
(
1
+
2
)
(
3
+
4
)
-
1
,
where upper indexes 〈〈Y→∞〉〉 and 〈〈β→∞〉〉 mean that those HOSC are derived with Y→∞ and β→∞ correspondingly. The index 〈〈∞〉〉 means that this HOSC is derived with both parameters tending to infinity.


## 2.1. Symbolic Approach for Sensitivity Analysis

The sensitivity computation of the transfer function ofT(s) with respect to change of the parameter X requires finding the corresponding derivative:
(1)S(H(s),X)=dH(s)H(s)XdX.The study of network response by means of formula1 is acceptable only when changes in its parameter X 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.If we consider the transfer function asT(s)=N(s)/D(s) formula 1 will transform into expression proposed by Hoang [5]:
(2)S(T(s),X)=XdN(s)N(s)-dN(s)D(s).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 [1] presented the sensitivity function T(s) in form of ratio of determinant and cofactors of nodal admittance or loop impedance matrix multiplication. For example, the sensitivity of T(s)=V3(s)/V1(s) of the circuit in Figure 1 is expressed as follows:
(3)SXT=-XΔ15Δ73ΔΔ13,
where X is a transfer impedance Z of current controlled voltage source (CCVS) or transfer admittance of voltage controlled current source (VCCS), Δ is a determinant of loop impedance matrix or Δ=Δ11 for nodal admittance matrix, and Δ11, Δ15, Δ13, and Δ73 are cofactors.Figure 1
Test circuit.Expression3 derived from 1 by means of Jacobi identity in form of
(4)Δ13,75Δ=Δ13Δ75-Δ15Δ73,
where Δ13,75 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 [5] 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. Thenth order cofactor is represented by a symbol Δr1,k1,r2,k2,…,rn,kn, where r1,r2,…,rn and k1,k2,…,kn are the numbers of deleted rows and columns, respectively. The summative cofactors are cofactors with at least one deletion in a form (a+d) or (a-d), which means the following: add or subtract the row (column) a to or from the row (column) d. For example the summative cofactors of the same matrix can be expressed as follows [15]:
(5)Δ(a+b)(c+d)=Δac-Δad-Δbc+Δbd,
where Δac=Δ(a+0)(c+0), Δad=Δ(a+0)(d+0), and so on.In accordance with GPEM basis [16–21] we consider the sums of row (column) indexes as a result of nullator (norator) connection to the network nodes corresponding to the indexes.Let us express the voltage gain function for circuit model in Figure2(a) as follows:
(6)T(s)=Vout(s)Vin(s)=Δ(1+2)(3+4)Δ(1+2)(1+2),
where HOSC in the numerator is the determinant of the network, in which the independent source and arbitrary response are replaced by oriented nullor [16] and HOSC in the denominator is the determinant of the network, in which the independent excitation and the arbitrary response are zero.Figure 2
The circuit for sensitivity computation of transfer functions: voltage gain (a), transfer admittance (b), transfer impedance (c), and current gain (d).
(a)
(b)
(c)
(d)The calculation procedure of the network determinants is based on the recursive usage of the parameter extraction formula [16]:
(7)Δ=χΔ(χ⟶∞)+Δ(χ=0),
where χ is a parameter of arbitrary linear circuit element and Δ(χ→∞) and Δ(χ=0) correspond to the determinants of the circuit matrix in which the parameter of extracted elements χ→∞ or χ=0, respectively.The parameter extraction formula7 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 [16]
(8)Δ=±Δn,
where χ is the determinant of the circuit matrix after the extraction procedure of the nullor number n. 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.Formula8 can be implemented not only for the nullor extraction from circuit model, but also for Coates flow-graph modification as shown in Section 4.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 fromhttp://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 [22].As can be seen from the comparison of symbolic expressions of notch filter transfer function [23] 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 [24] as example of compact size expression.


## 2.3. Derivation of Symbolic Sensitivity Computation Formulae

Let us express the transfer function6 by means of parameter extraction formula as
(9)T(s)=YΔ(1+2)(3+4),(7+8)(5+6)+Δ(1+2)(3+4)0YΔ(1+2)(1+2),(7+8)(5+6)+Δ(1+2)(1+2)0,
and then expressed the Hoang formula 2 as follows:
(10)SYT=Y·Δ(1+2)(1+2)∞Δ(1+2)(1+2)-Δ∞Δ.Note, that upper indexes〈〈0〉〉 and 〈〈∞〉〉 in formulas 9 and 10 mean that those HOSC are derived with Y=0 and Y→∞ correspondingly. In case of controlled source the operation χ→∞ leads to substitution of CS by the nullor.Formula10 is applicable for sensitivity computation of the arbitrary circuit elements. Just replace the 〈〈Y〉〉 by 〈〈χ〉〉.Now let us consider the transformation of Bode formula1 by substitution of expression 9 and cancellation by Δ(1+2)(1+2):
(11)SYT=YΔ1+23+4,7+85+6Δ1+21+2Δ1+21+2Δ1+23+4-YΔ(1+2)(3+4)Δ(1+2)(1+2),(7+8)(5+6)Δ(1+2)(1+2)Δ(1+2)(3+4).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Δ(1+2)(1+2) instead of determinant Δ:
(12)SYT=-YΔ(1+2)(5+6)Δ(1+2)(1+2),(7+8)(3+4)Δ(1+2)(1+2)Δ(1+2)(3+4).As can be seen,Δ(1+2)(1+2) is a numerator of transfer function from circuit input to controlling voltage of VCCS, Δ(1+2)(1+2),(7+8)(3+4) is a numerator of transfer function from controlled current source to circuit output, and Δ(1+2)(1+2)Δ(7+8)(3+4) 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 Table1 the generalized sensitivity formulae are proposed.Table 1
Sensitivity symbolic calculation formulae.
Number
Sensitivity type
Formula
1
S
χ
T
-
χ
Δ
(
1
+
2
)
(
5
+
6
)
Δ
(
1
+
2
)
(
1
+
2
)
,
(
7
+
8
)
(
3
+
4
)
Δ
(
1
+
2
)
(
1
+
2
)
Δ
(
1
+
2
)
(
3
+
4
)
2
S
χ
Y
3
S
χ
B
-
χ
Δ
(
1
+
2
)
(
5
+
6
)
Δ
(
7
+
8
)
(
3
+
4
)
Δ
Δ
(
1
+
2
)
(
3
+
4
)
4
S
χ
ZNote that the formula in first line of Table1 is appropriate not only for voltage gain, but also for transfer impedance function. Let us show it for circuit in Figure 2(b):
(13)Ys=IoutsVins=Δ1+23+4Δ1+21+2.Even if output nodes of transfer impedance function in Figure2(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:
(14)Z=VoutIin=Δ(1+2)(3+4)Δ,B=IoutIin=Δ(1+2)(3+4)Δ.The denominator of any sensitivity formulae in Table1 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:(15)SYT=-YΔ(1+2)(5+6)Δ(1+2)(1+2),(5+6)(3+4)Δ(1+2)(1+2)Δ(1+2)(3+4).The formula for transfer function sensitivity with respect to impedance parameter (χ=Z) 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 Table2. As it can be seen, the HOSC indexes 3 and 4 correspond to 7 and 8. Taking into account the identity χΔ(1+2)(5+6)=Δ(1+2)(7+8) we obtain the expressions consisting of just two HOSC.Table 2
The simplified sensitivity formulae in case of controlled source connecting with circuit output.
Number
Sensitivity type
Formula
Circuit
Circuit
1
S
μ
T
S
Z
K
Δ
(
1
+
2
)
(
1
+
2
)
,
(
3
+
4
)
(
3
+
4
)
Δ
(
1
+
2
)
(
1
+
2
)
2
S
μ
Z
S
Z
Z
Δ
(
3
+
4
)
(
3
+
4
)
Δ
3
S
Y
Y
S
β
Y
Δ
(
1
+
2
)
(
1
+
2
)
0
Δ
(
1
+
2
)
(
1
+
2
)
4
S
Y
B
S
β
B
Δ
0
Δ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 Table1 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:
(16)SβYZ=βYΔ1+25+6Δ7+89+10,1+21+2hhihhh×Δ11+123+4,1+21+2+Δ1+29+10hhhihh×Δ11+125+6,1+21+2Δ7+83+4,1+21+2×Δ2Δ1+23+4-1.Figure 3
The test circuit for second-order transfer impedance sensitivity computation.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 formula2. 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:
(17)
S
β
Y
T
=
β
Y
Δ
1
+
2
3
+
4
∞
Δ
2
-
Δ
Y
→
∞
Δ
1
+
2
3
+
4
β
→
∞
Δ
hhhhhh
-
Δ
1
+
2
3
+
4
Y
→
∞
Δ
β
→
∞
Δ
-
Δ
(
1
+
2
)
(
3
+
4
)
Δ
∞
Δ
hhhhhh
+
2
Δ
1
+
2
3
+
4
Δ
β
→
∞
Δ
Y
→
∞
×
Δ
2
Δ
(
1
+
2
)
(
3
+
4
)
-
1
,
where upper indexes 〈〈Y→∞〉〉 and 〈〈β→∞〉〉 mean that those HOSC are derived with Y→∞ and β→∞ correspondingly. The index 〈〈∞〉〉 means that this HOSC is derived with both parameters tending to infinity.


## 3. Symbolic Sensitivity Analysis Examples

Let us consider the voltage transfer function symbolic AC sensitivity analysis of second-order low-pass filter shown in Figure4 [25]. The values of parameters elements are as follows: C1=18,7571 nF; C2=9,37857 nF; G1=G2=1 mS; μ=1.Figure 4
Second-order low-pass filter.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 Table2. Taking into account the nodes numeration of filter circuit in Figure 4 the sensitivity function can be expressed as follows:
(18)SμT=Δ1+01+0,2+02+0Δ1+01+0=Δ11,22Δ11.The HOSC expansion is performed by means of GPEM. In accordance with7 the VCVS parameter is extracted from circuit as shown in Figure 7.The second summand of expression (Figure7) is correspond to the HOSC Δ11,22, 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 (Figure8) is equal to unity because both admittances connected in short circuit. After expansion of element G2 from second circuit in (Figure 8) we obtain the determinant of short-circuited admittance G1 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 G2 and G1.The symbolic determinant expression for the first summand in (Figure7) is shown in Figure 9For computation of expression (Figure9) the following procedures were done: deletion of capacities admittance sC2 connecting in parallel with nullator, parameter extraction of admittance sC1 connecting in series with nullator, followed by deletion of the admittance G1 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 “-1.”Now we need to substitute the symbolic expression from (Figures8 and 9) in 18 to obtain the voltage gain sensitivity function:
(19)SμT=s2C1C2+s[C1G2+C2(G1+G2)]+G1G2s2C1C2+s[(1-μ)C1G2+C2(G1+G2)]+G1G2.The derivation process of sensitivity functions with respect to other parameters is similar. For example we derive the first- and second-order sensitivitiesSG1T and SG1G2T:
(20)SG1T=-s2C1C2+s[C1G2+C1G2(1-μ)]s2C1C2+s[(1-μ)C1G2+C2(G1+G2)]+G1G2,SG1G2T=G1G2Δ1+01+4Δ1+44+3,1+01+0hhihhhhh×Δ4+32+0,1+01+0+Δ(1+0)(4+3)hhhhhihh×Δ(4+3)(1+4),(1+0)(1+0)Δ(1+4)(2+0),(1+0)(1+0)×Δ2Δ(1+0)(2+0)-1,
where
(21)Δ(1+0)(2+0)=μG1G2,Δ=sC2(sC1+G2(1-μ))+sC1(G2+G1)+G1G2,Δ(1+0)(1+4)=s2(C1C2)+s(C1G2(1-μ)+C2G2),Δ1+44+3,1+01+0=-sC2,Δ4+32+0,1+01+0=s-µC1-µG1,Δ(1+0)(4+3)=sC2G1,Δ(4+3)(1+4),(1+0)(1+0)=s(C1µ-C2),Δ1+42+0,1+01+0=-μG20.Using Mathcad software symbolic sensitivity simulation results versus frequencyf are shown in Figure 5.Figure 5
First- and second-order symbolic sensitivities versus frequencyf.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 [6], the method based on the modified Coates flow-graph [7], and the proposed technique. The test circuit is shown in Figure 6 [7].Figure 6
The STAR network.Figure 7Figure 8Figure 9Let us suppose that voltage transfer functionT(s) and the first- and second-order symbolic sensitivities with respect to parameters G4 and C2 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 expressionsT(s), SG4T, SsC2T, and SG4sC2T 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(22)T(s)=Δ1+03+0Δ,
where
(23)y1=s*C1;y2=s*C2;Δ1+03+0=y1*G7+G2+G5+y2*G1hhhhhhG7+G2+G5+y2-G4*G6+G3+(G7+G2)*G1*(G5+y2+G4);Δ=y1*G7+G5+G4*-G6+G2+y2*G3+G1+G5+G4+y2*G3+G1*G2-G6*G7.The sensitivity functions in Hoang form(24)SG4T=G4·Δ(1+0)(3+0)∞Δ(1+0)(3+0)-Δ∞Δ,
where
(25)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)∞=y1*-G3+G6+G1*G7+G2;Δ∞=G3+G1*G2-G6*G7+y1.(26)SsC2T=sC2·Δ(1+0)(3+0)∞Δ(1+0)(3+0)-Δ∞Δ,
where
(27)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)∞=G1*(G7+G2+y1);Δ∞=G2+y1*G3+G1-G6*G7.(28)SG4sC2T=G4*y2*Δ1+03+0∞*Δ*Δ-ΔsC2→∞*Δ1+03+0G4→∞Δ1+03+0*Δ*Δ-1hhihhh*Δ-Δ1+03+0sC2→∞*ΔG4→∞hhhihh*Δ-Δ1+03+0*Δ∞*Δ+2hhihhh*Δ1+03+0*ΔG4→∞*ΔsC2→∞hhihhΔ1+03+0∞*Δ*Δ-ΔsC2→∞*Δ1+03+0G4→∞×Δ1+03+0*Δ*Δ-1,
where
(29)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)G4→∞=y1*-G6+G3+G1*(G7+G2);Δ(1+0)(3+0)sC2→∞=G1*(G7+G2+y1);Δ(1+0)(3+0)∞=0;ΔG4→∞=(G3+G1)*(G2)-G6*(G7+y1);ΔsC2→∞=G2+y1*G3+G1-G6*G7;Δ∞=0.The sensitivity functions in Bode form(30)SG4T=-G4Δ(1+0)(1+2)Δ(1+2)(3+0),(1+0)(1+0)ΔΔ(1+0)(3+0),
where
(31)y1=s*C1;y2=s*C2;Δ1+01+2=G1*y2*-G7+G2*G5+y1*G2+y2*G3-G6*G7+G5+((G5+y2))*G2*G3-G6*G7;Δ(1+2)(3+0),(1+0)(1+0)=y1*G6+G3+G1,(32)SsC2K=-y2Δ(1+0)(2+3)Δ(2+3)(3+0),(1+0)(1+0)ΔΔ(1+0)(3+0),
where
(33)y1=s*C1;y2=s*C2;Δ1+02+3=y1*G7+G5*-G1+G3*G4+G4*G1*-G7+G2*G3-G6*G7-G1*G5*G7+G2;Δ(2+3)(3+0),(1+0)(1+0)=(-y1*(G6+G3+G1)).(34)SG4sC2T=G4y2Δ1+01+2Δ1+22+3,1+01+0hhhhhhh×Δ2+33+0,1+01+0+Δ(1+0)(2+3)hhhhhhh×Δ(2+3)(1+2),(1+0)(1+0)Δ(1+2)(3+0),(1+0)(1+0)×Δ2Δ(1+0)(2+0)-1,
where
(35)y1=s*C1;y2=s*C2;Δ1+01+2=G1*y2*-G7+G2*G5+y1*G2+y2*G3-G6*G7+G5+G5+y2*G2*G3-G6*G7;Δ1+22+3,1+01+0=y1*-G3+G1-G3+G1*G2-G6*G7;Δ2+33+0,1+01+0=-y1*G6+G3+G1;Δ1+02+3=y1*G7+G5*-G1+G3*G4+G4*G1*-G7+G2*G3-G6*G7-G1*G5*(G7+G2);Δ2+31+2,1+01+0=y1*G6-((G3+G1)*(G2)-G6*G7);Δ(1+2)(3+0),(1+0)(1+0)=y1*(G6+G3+G1).


### 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 [7]. 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 [7]) 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 (-sC2 and -G2) into node 3 are moved to node 0;(2)
the self-loopY33 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-G6 into node 5 is moved to node 4;(2)
the self-loopY55 transformed into the incoming edge into node 4.Third, we combine the two supporting nodes 3 and 5:(1)
the edge-G6 between nodes 3 and 5 transforms into the self-loop in the new node named by 3;(2)
the parallel codirectional edges-G6 and Y55 are combined into the edge Y55.Consequently we obtained the modified Coates flow-graph shown in Figure  8(b) (see publication [7]). 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 [7]. It leads to complications of method usage and high error probability.Implementation of the oriented nullor concept [16] 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 formulaeT(s), SG4T, SsC2T, and SG4sC2T presented here were replicated from [7] with correction of typing error in subexpression Y55 and sign of determinants:
(36)Ts=-Y22*G1*Y44+y1*y1hhi*G1+y1*Y55*G4×Y22*G6*Y44-Y22*G2*Y55-y1hhhh*y1*G6-y1*y2*Y55-1,
where
(37)Y22=(G4+G5+y1+y2);Y44=(G2+G7+y1);Y55=(G1+G3+G6);y1=s*C1;y2=s*C2.Let us substitute the subexpressionsY22, Y44, and Y55 into the transfer function 36:
(38)Ts=-(G4+G5+y1+y2)*G1hhih*(G2+G7+y1)+y1*y1*G1hhih+y1*G1+G3+G6*G4×G4+G5+y1+y2*G6hhhhh*G2+G7+y1-G4+G5+y1+y2hhhhh*G2*(G1+G3+G6)-y1*y1*G6hhhhh-y1*y2*G1+G3+G6-1.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 function36 includes the 8 pairs of cancellation sum-of-product terms:(39)
S
G
4
T
=
G
4
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
*
1
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
Y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
*
y
1
*
Y
55
,(40)
S
s
C
2
T
=
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
h
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
×
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
-
1
*
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
+
G
2
*
G
4
*
Y
55
+
G
2
*
G
1
*
y
1
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
h
l
h
*
G
1
+
y
1
*
Y
55
*
G
4
*
y
1
*
Y
55
,(41)
S
G
4
s
C
2
T
=
-
G
4
*
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
h
-
G
4
*
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
×
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
^
2
*
G
6
-
y
1
*
y
2
*
Y
55
-
1
*
-
(
-
Y
44
*
G
6
+
Y
55
*
G
2
)
*
y
1
*
Y
55
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
h
h
*
-
Y
44
*
(
G
4
*
G
6
-
G
1
*
y
2
)
+
G
2
*
G
4
*
Y
55
+
G
2
*
G
1
*
y
1
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
(
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
)
h
h
-
(
-
Y
44
*
G
6
+
Y
55
*
G
2
+
y
1
*
Y
55
)
*
y
1
*
Y
55
h
h
*
1
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
.


### 4.4. Symbolic Expressions Obtained by means of Two-Port Transimpedance Method

The sequences of expressions, derived by symbolic analysis program STAINS [6], 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 G4 and C2:
(42)T(s)=ZoiZii,
where
(43)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=G2(y1);x4=G4*d1;x5=G2+G6+y2+y2*d1;x6=-y2-(G4+G5+y1+y2)*d1;d2=-(G2+G7+y1)(y1);x7=x1+G4*d2;x8=x2+y2*d2;x9=-y1-(G4+G5+y1+y2)*d2;d3=-G1(x3);x10=G1+G4-x7*d3;x11=-x8*d3;x12=-G4-x9*d3;d4=-G6(x3);x13=x4-x7*d4;x14=x5-x8*d4;x15=x6-x9*d4;Z21=1*x12;Zoi=Z21;Z11=-1*x15;Zii=Z11.After substitution of the subexpressions into42 we obtain the transfer function expression as follows:
(44)
T
s
=
-
1
*
G
4
*
G
1
+
G
3
+
G
6
-
G
1
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhh
*
-
s
*
C
1
-
G
4
+
G
5
+
s
*
C
1
+
C
2
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhhhh
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
×
-
1
*
G
6
*
(
-
s
*
C
1
)
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
*
G
2
s
*
C
1
hhihhhhhhhhhhhhh
-
G
4
+
G
5
+
s
*
C
1
+
C
2
hhihhhhhhhhhhhhh
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhhhh
+
-
s
*
C
2
-
G
4
+
G
5
+
s
*
C
1
+
C
2
*
G
2
s
*
C
1
hhhhhhhhhhhhh
*
G
2
s
*
C
1
hhhhhhhhhh
*
G
2
s
*
C
1
*
G
1
+
G
3
+
G
6
-
1
.The expansion of the transfer function numerator and denominator leads to obtaining the expressions with 16 and 28 terms correspondingly. Just as formula36 the transfer function 42 includes the 8 pairs of cancellation sum-of-product terms:
(45)SG4T=G4*Zki*(Zok*Zii-Zik*Zoi)(Zoi*D00*Zii),
where
(46)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=-G1(x3);x4=G1+G4-x1*d1;x5=-x2*d1;x6=y1*d1-G4;x7=-(G2+G7+y1)*d1;d2=-G6(x3);x8=-x1*d2;x9=G2+G6+y2-x2*d2;x10=y1*d2-y2;x11=-G2-(G2+G7+y1)*d2;Z21=1*(x6*y1+x7*(G4+G5+y1+y2));Zoi=Z21;Z11=-1*(x10*y1+x11*(G4+G5+y1+y2));Zii=Z11;Z13=1*(x8*y1-x11*G4);Zik=Z11-Z13;Z23=-1*(x4*y1-x7*G4);Zok=Z21-Z23;Z31=1*(x6*x11-x7*x10);Zki=Z11-Z31;D00=x4*Z11+x8*Z21-G4*Z31.(47)SsC2T=y2*Zki*(Zik*Zoi-Zok*Zii)(Zoi*D00*Zii),
where
(48)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=-G1(x3);x4=G1+G4-x1*d1;x5=-x2*d1;x6=y1*d1-G4;x7=-(G2+G7+y1)*d1;d2=-G6(x3);x8=-x1*d2;x9=G2+G6+y2-x2*d2;x10=y1*d2-y2;x11=-G2-(G2+G7+y1)*d2;Z21=1*(x6*y1+x7*(G4+G5+y1+y2));Zoi=Z21;Z11=-1*(x10*y1+x11*(G4+G5+y1+y2));Zii=Z11;Z12=-0;Z13=1*(x8*y1-x11*G4);Zik=Z12-Z13;Z22=0;Z23=-1*(x4*y1-x7*G4);Zok=Z22-Z23;Z31=1*(x6*x11-x7*x10);Zki=Z21-Z31;D00=x4*Z11+x8*Z21-G4*Z31.


### 4.5. Efficiency Comparison

In Table3 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.Table 3
Comparison of symbolic analysis methods efficiency.
Symbolic function
Symbolic analysis method
M/D
A/S
T
(
s
)
GPEM22
14
20
MCFG36
17
13
TI42
29
34
S
G
4
T
GPEM (Hoang form)24
20
28
GPEM (Bode form)30
27
32
MCFG39
28
17
TI45
37
42
S
s
C
2
T
GPEM (Hoang form)26
19
26
GPEM (Bode form)32
27
33
MCFG40
34
20
TI47
37
43
S
G
4
s
C
2
T
GPEM (Hoang form)28
50
72
GPEM (Bode form)34
54
63
MCFG41
68
38
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.


## 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 [6], the method based on the modified Coates flow-graph [7], and the proposed technique. The test circuit is shown in Figure 6 [7].Figure 6
The STAR network.Figure 7Figure 8Figure 9Let us suppose that voltage transfer functionT(s) and the first- and second-order symbolic sensitivities with respect to parameters G4 and C2 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 expressionsT(s), SG4T, SsC2T, and SG4sC2T 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(22)T(s)=Δ1+03+0Δ,
where
(23)y1=s*C1;y2=s*C2;Δ1+03+0=y1*G7+G2+G5+y2*G1hhhhhhG7+G2+G5+y2-G4*G6+G3+(G7+G2)*G1*(G5+y2+G4);Δ=y1*G7+G5+G4*-G6+G2+y2*G3+G1+G5+G4+y2*G3+G1*G2-G6*G7.The sensitivity functions in Hoang form(24)SG4T=G4·Δ(1+0)(3+0)∞Δ(1+0)(3+0)-Δ∞Δ,
where
(25)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)∞=y1*-G3+G6+G1*G7+G2;Δ∞=G3+G1*G2-G6*G7+y1.(26)SsC2T=sC2·Δ(1+0)(3+0)∞Δ(1+0)(3+0)-Δ∞Δ,
where
(27)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)∞=G1*(G7+G2+y1);Δ∞=G2+y1*G3+G1-G6*G7.(28)SG4sC2T=G4*y2*Δ1+03+0∞*Δ*Δ-ΔsC2→∞*Δ1+03+0G4→∞Δ1+03+0*Δ*Δ-1hhihhh*Δ-Δ1+03+0sC2→∞*ΔG4→∞hhhihh*Δ-Δ1+03+0*Δ∞*Δ+2hhihhh*Δ1+03+0*ΔG4→∞*ΔsC2→∞hhihhΔ1+03+0∞*Δ*Δ-ΔsC2→∞*Δ1+03+0G4→∞×Δ1+03+0*Δ*Δ-1,
where
(29)y1=s*C1;y2=s*C2;Δ(1+0)(3+0)G4→∞=y1*-G6+G3+G1*(G7+G2);Δ(1+0)(3+0)sC2→∞=G1*(G7+G2+y1);Δ(1+0)(3+0)∞=0;ΔG4→∞=(G3+G1)*(G2)-G6*(G7+y1);ΔsC2→∞=G2+y1*G3+G1-G6*G7;Δ∞=0.The sensitivity functions in Bode form(30)SG4T=-G4Δ(1+0)(1+2)Δ(1+2)(3+0),(1+0)(1+0)ΔΔ(1+0)(3+0),
where
(31)y1=s*C1;y2=s*C2;Δ1+01+2=G1*y2*-G7+G2*G5+y1*G2+y2*G3-G6*G7+G5+((G5+y2))*G2*G3-G6*G7;Δ(1+2)(3+0),(1+0)(1+0)=y1*G6+G3+G1,(32)SsC2K=-y2Δ(1+0)(2+3)Δ(2+3)(3+0),(1+0)(1+0)ΔΔ(1+0)(3+0),
where
(33)y1=s*C1;y2=s*C2;Δ1+02+3=y1*G7+G5*-G1+G3*G4+G4*G1*-G7+G2*G3-G6*G7-G1*G5*G7+G2;Δ(2+3)(3+0),(1+0)(1+0)=(-y1*(G6+G3+G1)).(34)SG4sC2T=G4y2Δ1+01+2Δ1+22+3,1+01+0hhhhhhh×Δ2+33+0,1+01+0+Δ(1+0)(2+3)hhhhhhh×Δ(2+3)(1+2),(1+0)(1+0)Δ(1+2)(3+0),(1+0)(1+0)×Δ2Δ(1+0)(2+0)-1,
where
(35)y1=s*C1;y2=s*C2;Δ1+01+2=G1*y2*-G7+G2*G5+y1*G2+y2*G3-G6*G7+G5+G5+y2*G2*G3-G6*G7;Δ1+22+3,1+01+0=y1*-G3+G1-G3+G1*G2-G6*G7;Δ2+33+0,1+01+0=-y1*G6+G3+G1;Δ1+02+3=y1*G7+G5*-G1+G3*G4+G4*G1*-G7+G2*G3-G6*G7-G1*G5*(G7+G2);Δ2+31+2,1+01+0=y1*G6-((G3+G1)*(G2)-G6*G7);Δ(1+2)(3+0),(1+0)(1+0)=y1*(G6+G3+G1).


## 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 [7]. 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 [7]) 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 (-sC2 and -G2) into node 3 are moved to node 0;(2)
the self-loopY33 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-G6 into node 5 is moved to node 4;(2)
the self-loopY55 transformed into the incoming edge into node 4.Third, we combine the two supporting nodes 3 and 5:(1)
the edge-G6 between nodes 3 and 5 transforms into the self-loop in the new node named by 3;(2)
the parallel codirectional edges-G6 and Y55 are combined into the edge Y55.Consequently we obtained the modified Coates flow-graph shown in Figure  8(b) (see publication [7]). 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 [7]. It leads to complications of method usage and high error probability.Implementation of the oriented nullor concept [16] 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 formulaeT(s), SG4T, SsC2T, and SG4sC2T presented here were replicated from [7] with correction of typing error in subexpression Y55 and sign of determinants:
(36)Ts=-Y22*G1*Y44+y1*y1hhi*G1+y1*Y55*G4×Y22*G6*Y44-Y22*G2*Y55-y1hhhh*y1*G6-y1*y2*Y55-1,
where
(37)Y22=(G4+G5+y1+y2);Y44=(G2+G7+y1);Y55=(G1+G3+G6);y1=s*C1;y2=s*C2.Let us substitute the subexpressionsY22, Y44, and Y55 into the transfer function 36:
(38)Ts=-(G4+G5+y1+y2)*G1hhih*(G2+G7+y1)+y1*y1*G1hhih+y1*G1+G3+G6*G4×G4+G5+y1+y2*G6hhhhh*G2+G7+y1-G4+G5+y1+y2hhhhh*G2*(G1+G3+G6)-y1*y1*G6hhhhh-y1*y2*G1+G3+G6-1.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 function36 includes the 8 pairs of cancellation sum-of-product terms:(39)
S
G
4
T
=
G
4
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
*
1
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
Y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
*
y
1
*
Y
55
,(40)
S
s
C
2
T
=
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
h
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
×
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
-
1
*
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
+
G
2
*
G
4
*
Y
55
+
G
2
*
G
1
*
y
1
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
h
l
h
*
G
1
+
y
1
*
Y
55
*
G
4
*
y
1
*
Y
55
,(41)
S
G
4
s
C
2
T
=
-
G
4
*
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
h
-
G
4
*
y
2
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
×
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
^
2
*
G
6
-
y
1
*
y
2
*
Y
55
-
1
*
-
(
-
Y
44
*
G
6
+
Y
55
*
G
2
)
*
y
1
*
Y
55
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
h
h
*
-
Y
44
*
(
G
4
*
G
6
-
G
1
*
y
2
)
+
G
2
*
G
4
*
Y
55
+
G
2
*
G
1
*
y
1
-
Y
22
*
G
1
*
Y
44
+
y
1
*
y
1
*
G
1
+
y
1
*
Y
55
*
G
4
(
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
)
h
h
-
(
-
Y
44
*
G
6
+
Y
55
*
G
2
+
y
1
*
Y
55
)
*
y
1
*
Y
55
h
h
*
1
-
Y
44
*
G
4
*
G
6
-
G
1
*
y
2
-
G
2
*
G
4
*
Y
55
-
G
2
*
G
1
*
y
1
Y
22
*
G
6
*
Y
44
-
Y
22
*
G
2
*
Y
55
-
y
1
*
y
1
*
G
6
-
y
1
*
y
2
*
Y
55
.


## 4.4. Symbolic Expressions Obtained by means of Two-Port Transimpedance Method

The sequences of expressions, derived by symbolic analysis program STAINS [6], 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 G4 and C2:
(42)T(s)=ZoiZii,
where
(43)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=G2(y1);x4=G4*d1;x5=G2+G6+y2+y2*d1;x6=-y2-(G4+G5+y1+y2)*d1;d2=-(G2+G7+y1)(y1);x7=x1+G4*d2;x8=x2+y2*d2;x9=-y1-(G4+G5+y1+y2)*d2;d3=-G1(x3);x10=G1+G4-x7*d3;x11=-x8*d3;x12=-G4-x9*d3;d4=-G6(x3);x13=x4-x7*d4;x14=x5-x8*d4;x15=x6-x9*d4;Z21=1*x12;Zoi=Z21;Z11=-1*x15;Zii=Z11.After substitution of the subexpressions into42 we obtain the transfer function expression as follows:
(44)
T
s
=
-
1
*
G
4
*
G
1
+
G
3
+
G
6
-
G
1
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhh
*
-
s
*
C
1
-
G
4
+
G
5
+
s
*
C
1
+
C
2
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhhhh
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
×
-
1
*
G
6
*
(
-
s
*
C
1
)
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
*
G
2
s
*
C
1
hhihhhhhhhhhhhhh
-
G
4
+
G
5
+
s
*
C
1
+
C
2
hhihhhhhhhhhhhhh
*
-
G
2
+
G
7
+
s
*
C
1
s
*
C
1
hhhhhhhhhh
+
-
s
*
C
2
-
G
4
+
G
5
+
s
*
C
1
+
C
2
*
G
2
s
*
C
1
hhhhhhhhhhhhh
*
G
2
s
*
C
1
hhhhhhhhhh
*
G
2
s
*
C
1
*
G
1
+
G
3
+
G
6
-
1
.The expansion of the transfer function numerator and denominator leads to obtaining the expressions with 16 and 28 terms correspondingly. Just as formula36 the transfer function 42 includes the 8 pairs of cancellation sum-of-product terms:
(45)SG4T=G4*Zki*(Zok*Zii-Zik*Zoi)(Zoi*D00*Zii),
where
(46)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=-G1(x3);x4=G1+G4-x1*d1;x5=-x2*d1;x6=y1*d1-G4;x7=-(G2+G7+y1)*d1;d2=-G6(x3);x8=-x1*d2;x9=G2+G6+y2-x2*d2;x10=y1*d2-y2;x11=-G2-(G2+G7+y1)*d2;Z21=1*(x6*y1+x7*(G4+G5+y1+y2));Zoi=Z21;Z11=-1*(x10*y1+x11*(G4+G5+y1+y2));Zii=Z11;Z13=1*(x8*y1-x11*G4);Zik=Z11-Z13;Z23=-1*(x4*y1-x7*G4);Zok=Z21-Z23;Z31=1*(x6*x11-x7*x10);Zki=Z11-Z31;D00=x4*Z11+x8*Z21-G4*Z31.(47)SsC2T=y2*Zki*(Zik*Zoi-Zok*Zii)(Zoi*D00*Zii),
where
(48)y1=s*C1;y2=s*C2;x1=-G1;x2=-G2-G6;x3=G1+G3+G6;d1=-G1(x3);x4=G1+G4-x1*d1;x5=-x2*d1;x6=y1*d1-G4;x7=-(G2+G7+y1)*d1;d2=-G6(x3);x8=-x1*d2;x9=G2+G6+y2-x2*d2;x10=y1*d2-y2;x11=-G2-(G2+G7+y1)*d2;Z21=1*(x6*y1+x7*(G4+G5+y1+y2));Zoi=Z21;Z11=-1*(x10*y1+x11*(G4+G5+y1+y2));Zii=Z11;Z12=-0;Z13=1*(x8*y1-x11*G4);Zik=Z12-Z13;Z22=0;Z23=-1*(x4*y1-x7*G4);Zok=Z22-Z23;Z31=1*(x6*x11-x7*x10);Zki=Z21-Z31;D00=x4*Z11+x8*Z21-G4*Z31.


## 4.5. Efficiency Comparison

In Table3 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.Table 3
Comparison of symbolic analysis methods efficiency.
Symbolic function
Symbolic analysis method
M/D
A/S
T
(
s
)
GPEM22
14
20
MCFG36
17
13
TI42
29
34
S
G
4
T
GPEM (Hoang form)24
20
28
GPEM (Bode form)30
27
32
MCFG39
28
17
TI45
37
42
S
s
C
2
T
GPEM (Hoang form)26
19
26
GPEM (Bode form)32
27
33
MCFG40
34
20
TI47
37
43
S
G
4
s
C
2
T
GPEM (Hoang form)28
50
72
GPEM (Bode form)34
54
63
MCFG41
68
38
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.


## 5. Conclusions

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.



---

*Source: 328517-2015-01-18.xml*