Mason's gain formula (MGF) is a method for finding the transfer function of a linear signal-flow graph (SFG). The formula was derived by Samuel Jefferson Mason, for whom it is named. MGF is an alternate method to finding the transfer function algebraically by labeling each signal, writing down the equation for how that signal depends on other signals, and then solving the multiple equations for the output signal in terms of the input signal. MGF provides a step by step method to obtain the transfer function from a SFG. Often, MGF can be determined by inspection of the SFG. The method can easily handle SFGs with many variables and loops including loops with inner loops. MGF comes up often in the context of control systems, microwave circuits and digital filters because these are often represented by SFGs.
Formula
The gain formula is as follows:
G = y out y in = ∑ k = 1 N G k Δ k Δ {\displaystyle G={\frac {y_{\text{out}}}{y_{\text{in}}}}={\frac {\sum _{k=1}^{N}{G_{k}\Delta _{k}}}{\Delta \ }}} Δ = 1 − ∑ L i + ∑ L i L j − ∑ L i L j L k + ⋯ + ( − 1 ) m ∑ ⋯ + ⋯ {\displaystyle \Delta =1-\sum L_{i}+\sum L_{i}L_{j}-\sum L_{i}L_{j}L_{k}+\cdots +(-1)^{m}\sum \cdots +\cdots }where:
- Δ = the determinant of the graph.
- yin = input-node variable
- yout = output-node variable
- G = complete gain between yin and yout
- N = total number of forward paths between yin and yout
- Gk = path gain of the kth forward path between yin and yout
- Li = loop gain of each closed loop in the system
- LiLj = product of the loop gains of any two non-touching loops (no common nodes)
- LiLjLk = product of the loop gains of any three pairwise nontouching loops
- Δk = the cofactor value of Δ for the kth forward path, with the loops touching the kth forward path removed.
Definitions
Source:2
- Path: a continuous set of branches traversed in the direction that they indicate.
- Forward path: A path from an input node to an output node in which no node is touched more than once.
- Loop: A path that originates and ends on the same node in which no node is touched more than once.
- Path gain: the product of the gains of all the branches in the path.
- Loop gain: the product of the gains of all the branches in the loop.
Procedure to find the solution
- Make a list of all forward paths, and their gains, and label these Gk.
- Make a list of all the loops and their gains, and label these Li (for i loops). Make a list of all pairs of non-touching loops, and the products of their gains (LiLj). Make a list of all pairwise non-touching loops taken three at a time (LiLjLk), then four at a time, and so forth, until there are no more.
- Compute the determinant Δ and cofactors Δk.
- Apply the formula.
Examples
Circuit containing two-port
The transfer function from Vin to V2 is desired.
There is only one forward path:
- Vin to V1 to I2 to V2 with gain G 1 = − y 21 R L {\displaystyle G_{1}=-y_{21}R_{L}\,}
There are three loops:
- V1 to I1 to V1 with gain L 1 = − R in y 11 {\displaystyle L_{1}=-R_{\text{in}}y_{11}\,}
- V2 to I2 to V2 with gain L 2 = − R L y 22 {\displaystyle L_{2}=-R_{L}y_{22}\,}
- V1 to I2 to V2 to I1 to V1 with gain L 3 = y 21 R L y 12 R in {\displaystyle L_{3}=y_{21}R_{L}y_{12}R_{\text{in}}\,}
Digital IIR biquad filter
Digital filters are often diagramed as signal flow graphs.
There are two loops- L 1 = − a 1 Z − 1 {\displaystyle L_{1}=-a_{1}Z^{-1}\,}
- L 2 = − a 2 Z − 2 {\displaystyle L_{2}=-a_{2}Z^{-2}\,}
- G 0 = b 0 {\displaystyle G_{0}=b_{0}\,}
- G 1 = b 1 Z − 1 {\displaystyle G_{1}=b_{1}Z^{-1}\,}
- G 2 = b 2 Z − 2 {\displaystyle G_{2}=b_{2}Z^{-2}\,}
Servo
The signal flow graph has six loops. They are:
- L 0 = − β s M {\displaystyle L_{0}=-{\frac {\beta }{sM}}\,}
- L 1 = − ( R M + R S ) s L M {\displaystyle L_{1}={\frac {-(R_{M}+R_{S})}{sL_{M}}}\,}
- L 2 = − G M K M s 2 L M M {\displaystyle L_{2}=\,{\frac {-G_{M}K_{M}}{s^{2}L_{M}M}}}
- L 3 = − K C R S s L M {\displaystyle L_{3}={\frac {-K_{C}R_{S}}{sL_{M}}}\,}
- L 4 = − K V K C K M G T s 2 L M M {\displaystyle L_{4}={\frac {-K_{V}K_{C}K_{M}G_{T}}{s^{2}L_{M}M}}\,}
- L 5 = − K P K V K C K M s 3 L M M {\displaystyle L_{5}={\frac {-K_{P}K_{V}K_{C}K_{M}}{s^{3}L_{M}M}}\,}
There is one forward path:
- g 0 = K P K V K C K M s 3 L M M {\displaystyle g_{0}={\frac {K_{P}K_{V}K_{C}K_{M}}{s^{3}L_{M}M}}\,}
The forward path touches all the loops therefore the co-factor Δ 0 = 1 {\displaystyle \Delta _{0}=1}
And the gain from input to output is θ L θ C = g 0 Δ 0 Δ {\displaystyle {\frac {\theta _{L}}{\theta _{C}}}={\frac {g_{0}\Delta _{0}}{\Delta }}\,}
Equivalent matrix form
Mason's rule can be stated in a simple matrix form. Assume T {\displaystyle \mathbf {T} } is the transient matrix of the graph where t n m = [ T ] n m {\displaystyle t_{nm}=\left[\mathbf {T} \right]_{nm}} is the sum transmittance of branches from node m toward node n. Then, the gain from node m to node n of the graph is equal to u n m = [ U ] n m {\displaystyle u_{nm}=\left[\mathbf {U} \right]_{nm}} , where
U = ( I − T ) − 1 {\displaystyle \mathbf {U} =\left(\mathbf {I} -\mathbf {T} \right)^{-1}} ,and I {\displaystyle \mathbf {I} } is the identity matrix.
Mason's Rule is also particularly useful for deriving the z-domain transfer function of discrete networks that have inner feedback loops embedded within outer feedback loops (nested loops). If the discrete network can be drawn as a signal flow graph, then the application of Mason's Rule will give that network's z-domain H(z) transfer function.
Complexity and computational applications
Mason's Rule can grow factorially, because the enumeration of paths in a directed graph grows dramatically. To see this consider the complete directed graph on n {\displaystyle n} vertices, having an edge between every pair of vertices. There is a path form y in {\displaystyle y_{\text{in}}} to y out {\displaystyle y_{\text{out}}} for each of the ( n − 2 ) ! {\displaystyle (n-2)!} permutations of the intermediate vertices. Thus Gaussian elimination is more efficient in the general case.
Yet Mason's rule characterizes the transfer functions of interconnected systems in a way which is simultaneously algebraic and combinatorial, allowing for general statements and other computations in algebraic systems theory. While numerous inverses occur during Gaussian elimination, Mason's rule naturally collects these into a single quasi-inverse. General form is
p 1 − q , {\displaystyle {\frac {p}{1-q}},}Where as described above, q {\displaystyle q} is a sum of cycle products, each of which typically falls into an ideal (for example, the strictly causal operators). Fractions of this form make a subring R ( 1 + ⟨ L i ⟩ ) − 1 {\displaystyle R(1+\langle L_{i}\rangle )^{-1}} of the rational function field. This observation carries over to the noncommutative case,3 even though Mason's rule itself must then be replaced by Riegle's rule.
See also
Notes
- Bolton, W. Newnes (1998). Control Engineering Pocketbook. Oxford: Newnes.
- Van Valkenburg, M. E. (1974). Network Analysis (3rd ed.). Englewood Cliffs, NJ: Prentice-Hall.
References
Mason, Samuel J. (July 1956). "Feedback Theory - Further Properties of Signal Flow Graphs" (PDF). Proceedings of the IRE. 44 (7): 920–926. doi:10.1109/jrproc.1956.275147. hdl:1721.1/4778. S2CID 18184015. http://dspace.mit.edu/bitstream/1721.1/4778/1/RLE-TR-303-15342712.pdf ↩
Kuo, Benjamin C. (1967). Automatic Control Systems (2nd ed.). Prentice-Hall. pp. 59–60. ↩
Pliam, J.O.; Lee, E.B. (1995). "On the global properties of interconnected systems". IEEE Trans. Circuits and Syst. I. 42 (12): 1013–1017. doi:10.1109/81.481196. /wiki/Doi_(identifier) ↩