View source for Transfer function matrix

{{Short description|Matrix relating system inputs and outputs}}

In [[Control system|control system]] theory, and various branches of engineering, a '''transfer function matrix''', or just '''transfer matrix''' is a generalisation of the [[Transfer function|transfer function]]s of [[Single-input single-output system|single-input single-output]] (SISO) systems to multiple-input and multiple-output (MIMO)<!-- Do not link this to [[MIMO]], that is the wrong article --> systems.<ref>Chen, p. 1038</ref>  The [[Matrix (mathematics)|matrix]] relates the outputs of the system to its inputs.  It is a particularly useful construction for linear time-invariant (LTI) systems because it can be expressed in terms of the [[S-plane|s-plane]].

In some systems, especially ones consisting entirely of [[Engineering:Passivity (engineering)|passive]] components, it can be ambiguous which variables are inputs and which are outputs.  In electrical engineering, a common scheme is to gather all the voltage variables on one side and all the current variables on the other regardless of which are inputs or outputs.  This results in all the elements of the transfer matrix being in units of [[Physics:Electrical impedance|impedance]].  The concept of impedance (and hence impedance matrices) has been borrowed into other energy domains by analogy, especially mechanics and acoustics.

Many control systems span several different energy domains.  This requires transfer matrices with elements in mixed units.  This is needed both to describe [[Engineering:Transducer|transducer]]s that make connections between domains and to describe the system as a whole.  If the matrix is to properly model energy flows in the system, compatible variables must be chosen to allow this.

== General ==

A MIMO system with {{mvar|m}} outputs and {{mvar|n}} inputs is represented by a {{math|''m'' × ''n''}} matrix.  Each entry in the matrix is in the form of a transfer function relating an output to an input.  For example, for a three-input, two-output system, one might write,

:<math>
\begin{bmatrix} y_1 \\ y_2 \end{bmatrix} = 
\begin{bmatrix} 
g_{11} & g_{12} & g_{13} \\
g_{21} & g_{22} & g_{23}
\end{bmatrix}
\begin{bmatrix} u_1 \\ u_2 \\ u_3 \end{bmatrix}
</math>

where the {{mvar|u<sub>n</sub>}} are the inputs, the {{mvar|y<sub>m</sub>}} are the outputs, and the {{mvar|g<sub>mn</sub>}} are the transfer functions.  This may be written more succinctly in matrix operator notation as,

:<math> \mathbf Y = \mathbf G \mathbf U </math>

where {{math|'''Y'''}} is a column vector of the outputs, {{math|'''G'''}} is a matrix of the transfer functions, and {{math|'''U'''}} is a column vector of the inputs.

In many cases, the system under consideration is a linear time-invariant (LTI) system.  In such cases, it is convenient to express the transfer matrix in terms of the [[Laplace transform]] (in the case of continuous time variables) or the [[Z-transform|z-transform]] (in the case of discrete time variables) of the variables.  This may be indicated by writing, for instance,

:<math> \mathbf Y (s) = \mathbf G (s) \mathbf U (s) </math>

which indicates that the variables and matrix are in terms of {{mvar|s}}, the complex frequency variable of the [[S-plane|s-plane]] arising from Laplace transforms, rather than time.  The examples in this article are all assumed to be in this form, although that is not explicitly indicated for brevity.  For discrete time systems {{mvar|s}} is replaced by {{mvar|z}} from the z-transform, but this makes no difference to subsequent analysis.  The matrix is particularly useful when it is a '''proper rational matrix''', that is, all its elements are proper rational functions.  In this case the [[State-space representation|state-space representation]] can be applied.<ref>{{multiref|Levine, p. 481|Chen, pp. 1037–1038}}</ref>

In systems engineering, the overall system transfer matrix {{math|'''G''' (''s'')}} is decomposed into two parts: {{math|'''H''' (''s'')}} representing the system being controlled, and {{math|'''C'''(''s'')}} representing the control system.  {{math|'''C''' (''s'')}} takes as its inputs the inputs of {{math|'''G''' (''s'')}} and the outputs of {{math|'''H''' (''s'')}}.  The outputs of {{math|'''C''' (''s'')}} form the inputs for {{math|'''H''' (''s'')}}.<ref>Kavanagh, p. 350</ref>

== Electrical systems ==
In electrical systems it is often the case that the distinction between input and output variables is ambiguous.  They can be either, depending on circumstance and point of view.  In such cases the concept of [[Engineering:Port (circuit theory)|port]] (a place where energy is transferred from one system to another) can be more useful than input and output.  It is customary to define two variables for each port ({{mvar|p}}): the [[Physics:Voltage|voltage]] across it ({{mvar|V<sub>p</sub>}}) and the [[Physics:Electric current|current]] entering it ({{mvar|I<sub>p</sub>}}).  For instance, the transfer matrix of a [[Two-port network|two-port network]] can be defined as follows,

:<math>
\begin{bmatrix} V_1 \\ V_2 \end{bmatrix} = 
\begin{bmatrix}
z_{11} & z_{12} \\
z_{21} & z_{22} \\
\end{bmatrix}
\begin{bmatrix} I_1 \\ I_2 \end{bmatrix}
</math>

where the {{mvar|z<sub>mn</sub>}} are called the [[Impedance parameters|impedance parameters]], or ''z''-parameters.  They are so called because they are in units of [[Physics:Electrical impedance|impedance]] and relate port currents to a port voltage.  The z-parameters are not the only way that transfer matrices are defined for two-port networks.  There are six basic matrices that relate voltages and currents each with advantages for particular system network topologies.<ref>{{multiref|Chen, pp. 54–55|Iyer, p. 240|Bakshi & Bakshi, p. 420}}</ref>  However, only two of these can be extended beyond two ports to an arbitrary number of ports.  These two are the ''z''-parameters and their inverse, the [[Admittance parameters|admittance parameters]] or ''y''-parameters.<ref>Choma, p. 197</ref>
thumb|upright=0.7|Voltage divider circuit
To understand the relationship between port voltages and currents and inputs and outputs, consider the simple voltage divider circuit.  If we only wish to consider the output voltage ({{math|''V''<sub>2</sub>}}) resulting from applying the input voltage ({{math|''V''<sub>1</sub>}}) then the transfer function can be expressed as,

: <math>
\begin{bmatrix} V_2 \end{bmatrix} = 
\begin{bmatrix} 
\dfrac{R_2}{R_1 + R_2}
\end{bmatrix}
\begin{bmatrix} V_1 \end{bmatrix}
</math>

which can be considered the trivial case of a 1×1 transfer matrix.  The expression correctly predicts the output voltage if there is no current leaving port 2, but is increasingly inaccurate as the load increases.  If, however, we attempt to use the circuit in reverse, driving it with a voltage at port 2 and calculate the resulting voltage at port 1 the expression gives completely the wrong result even with no load on port 1.  It predicts a greater voltage at port 1 than was applied at port 2, an impossibility with a purely resistive circuit like this one.  To correctly predict the behaviour of the circuit, the currents entering or leaving the ports must also be taken into account, which is what the transfer matrix does.<ref>Yang & Lee, pp. 37–38</ref>  The impedance matrix for the voltage divider circuit is,

:<math>
\begin{bmatrix} V_1 \\ V_2 \end{bmatrix} = 
\begin{bmatrix}
R_1 + R_2 & R_2 \\
R_2 & R_2
\end{bmatrix}
\begin{bmatrix} I_1 \\ I_2 \end{bmatrix}
</math>

which fully describes its behaviour under all input and output conditions.<ref>Bessai, pp. 4–5</ref>

At [[Physics:Microwave|microwave]] frequencies, none of the transfer matrices based on port voltages and currents are convenient to use in practice.  Voltage is difficult to measure directly, current next to impossible, and the open circuits and short circuits required by the measurement technique cannot be achieved with any accuracy.  For [[Physics:Waveguide|waveguide]] implementations, circuit voltage and current are entirely meaningless.  Transfer matrices using different sorts of variables are used instead.  These are the [[Physics:Power|powers]] transmitted into, and reflected from a port which are readily measured in the [[Engineering:Transmission line|transmission line]] technology used in [[Engineering:Distributed-element circuit|distributed-element circuit]]s in the microwave band.  The most well known and widely used of these sorts of parameters is the [[Scattering parameters|scattering parameters]], or s-parameters.<ref>{{multiref|Nguyen, p. 271|Bessai, p. 1}}</ref>

==Mechanical and other systems==
[[File:VIEW TO SOUTHWEST. DETAIL, GEAR TRAIN SYSTEM, CONTROL CABIN INTERIOR. - Gianella Bridge, Spanning Sacramento River at State Highway 32, Hamilton City, Glenn County, CA HAER CAL,11-HAMCI.V,1-22.tif|thumb|A gear train in the control cabin of the former Gianella Bridge which operated this swing bridge.  Gear trains are two-ports.]]
The concept of [[Physics:Mechanical impedance|impedance]] can be extended into the mechanical, and other domains through a mechanical-electrical analogy, hence the impedance parameters, and other forms of 2-port network parameters, can be extended to the mechanical domain also.  To do this an effort variable and a flow variable are made analogues of voltage and current respectively.  For mechanical systems under translation these variables are [[Physics:Force|force]] and [[Physics:Velocity|velocity]] respectively.<ref>Busch-Vishniac, pp. 19–20</ref>

Expressing the behaviour of a mechanical component as a two-port or multi-port with a transfer matrix is a useful thing to do because, like electrical circuits, the component can often be operated in reverse and its behaviour is dependent on the loads at the inputs and outputs.  For instance, a [[Engineering:Gear train|gear train]] is often characterised simply by its gear ratio, a SISO transfer function.  However, the gearbox output [[Physics:Shaft (mechanical engineering)|shaft]] can be driven round to turn the input shaft requiring a MIMO analysis.  In this example the effort and flow variables are [[Physics:Torque|torque]] {{mvar|T}} and [[Physics:Angular velocity|angular velocity]] {{mvar|ω}} respectively.  The transfer matrix in terms of z-parameters will look like,

:<math>
\begin{bmatrix} T_1 \\ T_2 \end{bmatrix} = 
\begin{bmatrix}
z_{11} & z_{12} \\
z_{21} & z_{22}
\end{bmatrix}
\begin{bmatrix} \omega_1 \\ \omega_2 \end{bmatrix}
</math>

However, the z-parameters are not necessarily the most convenient for characterising gear trains.  A gear train is the analogue of an electrical [[Engineering:Transformer|transformer]] and the h-parameters (''hybrid'' parameters) better describe transformers because they directly include the turns ratios (the analogue of gear ratios).<ref>Olsen, pp. 239–240</ref>  The gearbox transfer matrix in h-parameter format is,

:<math>
\begin{bmatrix} T_1 \\ \omega_2 \end{bmatrix} = 
\begin{bmatrix}
h_{11} & h_{12} \\
h_{21} & h_{22}
\end{bmatrix}
\begin{bmatrix} \omega_1 \\ T_2 \end{bmatrix}
</math>

:where
:{{math|''h''<sub>21</sub>}} is the velocity ratio of the gear train with no load on the output,
:{{math|''h''<sub>12</sub>}} is the reverse direction torque ratio of the gear train with input shaft clamped, equal to the forward velocity ratio for an ideal gearbox,
:{{math|''h''<sub>11</sub>}} is the input rotational mechanical impedance with no load on the output shaft, zero for an ideal gearbox, and,
:{{math|''h''<sub>22</sub>}} is the output rotational mechanical [[Physics:Admittance|admittance]] with the input shaft clamped.

For an ideal gear train with no losses (friction, distortion etc), this simplifies to,

:<math>
\begin{bmatrix} T_1 \\ \omega_2 \end{bmatrix} =
\begin{bmatrix} 
0 & N \\
N & 0
\end{bmatrix}
\begin{bmatrix} \omega_1 \\ T_2 \end{bmatrix}
</math>

where {{mvar|N}} is the gear ratio.<ref>{{multiref|Busch-Vishniac, p. 20|Koenig & Blackwell, p. 170}}</ref>

== Transducers and actuators ==
[[File:Filterbaustein MF200 1 800 293.jpg|thumb|A mechanical filter opened to show the mechanical-electrical transducers at either end]]
In a system that consists of multiple energy domains, transfer matrices are required that can handle components with ports in different domains.  In robotics and [[Mechatronics|mechatronics]], [[Engineering:Actuator|actuator]]s are required.  These usually consist of a [[Engineering:Transducer|transducer]] converting, for instance, signals from the control system in the electrical domain into motion in the mechanical domain.  The control system also requires [[Engineering:Sensor|sensor]]s that detect the motion and convert it back into the electrical domain through another transducer so that the motion can be properly controlled through a feedback loop. Other sensors in the system may be transducers converting yet other energy domains into electrical signals, such as optical, audio, thermal, fluid flow and chemical.  Another application is the field of [[Physics:Mechanical filter|mechanical filter]]s which require transducers between the electrical and mechanical domains in both directions.

A simple example is an electromagnetic electromechanical actuator driven by an electronic controller.  This requires a transducer with an input port in the electrical domain and an output port in the mechanical domain.  This might be represented simplistically by a SISO transfer function, but for similar reasons to those already stated, a more accurate representation is achieved with a two-input, two-output MIMO transfer matrix.  In the z-parameters, this takes the form,

:<math>
\begin{bmatrix} V \\ F \end{bmatrix} = 
\begin{bmatrix}
z_{11} & z_{12} \\
z_{21} & z_{22}
\end{bmatrix}
\begin{bmatrix} I \\ v \end{bmatrix}
</math>

where {{mvar|F}} is the force applied to the actuator and {{mvar|v}} is the resulting velocity of the actuator.  The impedance parameters here are a mixture of units; {{math|''z''<sub>11</sub>}} is an electrical impedance, {{math|''z''<sub>22</sub>}} is a mechanical impedance and the other two are transimpedances in a hybrid mix of units.<ref>Pierce, p. 200</ref>

== Acoustic systems ==
[[Physics:Acoustics|Acoustic systems]] are a subset of [[Physics:Fluid dynamics|fluid dynamics]], and in both fields the primary input and output variables are [[Physics:Pressure|pressure]], {{mvar|P}}, and [[Physics:Volumetric flow rate|volumetric flow rate]], {{mvar|Q}}, except in the case of sound travelling through solid components.  In the latter case, the primary variables of mechanics, force and velocity, are more appropriate.  An example of a two-port acoustic component is a [[Filter (signal processing)|filter]] such as a [[Engineering:Muffler|muffler]] on an [[Engineering:Exhaust system|exhaust system]].  A transfer matrix representation of it may look like,

:<math>
\begin{bmatrix} P_2 \\ Q_2 \end{bmatrix} = 
\begin{bmatrix} 
T_{11} & T_{12} \\
T_{21} & T_{22}
\end{bmatrix}
\begin{bmatrix} P_1 \\ Q_1 \end{bmatrix}
</math>

Here, the {{math|''T<sub>mn</sub>''}} are the transmission parameters, also known as ABCD-parameters.  The component can be just as easily described by the z-parameters, but transmission parameters have a mathematical advantage when dealing with a system of two-ports that are connected in a cascade of the output of one into the input port of another.  In such cases the overall transmission parameters are found simply by the matrix multiplication of the transmission parameter matrices of the constituent components.<ref>Munjal, p. 81</ref>

== Compatible variables ==
[[File:Pneumatic Rack and Pinion Actuators.JPG|thumb|A pneumatic rack and pinion actuator controlling a valve in a water pipe.  The actuator is a two-port device that converts from the pneumatic domain to the mechanical domain.  Together with the valve itself it comprises a three-port system; the pneumatic control port and the fluid flow input and output water pipe ports of the valve.]]  
When working with mixed variables from different energy domains consideration needs to be given on which variables to consider analogous.  The choice depends on what the analysis is intended to achieve.  If it is desired to correctly model energy flows throughout the entire system then a pair of variables whose product is power (power conjugate variables) in one energy domain must map to power conjugate variables in other domains.  Power conjugate variables are not unique so care needs to be taken to use the same mapping of variables throughout the system.<ref>Busch-Vishniac, p. 18</ref>

A common mapping (used in some of the examples in this article) maps the effort variables (ones that initiate an action) from each domain together and maps the flow variables (ones that are a property of an action) from each domain together.  Each pair of effort and flow variables is power conjugate.  This system is known as the [[Engineering:Impedance analogy|impedance analogy]] because a ratio of the effort to the flow variable in each domain is analogous to electrical impedance.<ref>Busch-Vishniac, p. 20</ref>

There are two other power conjugate systems on the same variables that are in use.  The [[Engineering:Mobility analogy|mobility analogy]] maps mechanical force to electric current instead of voltage.  This analogy is widely used by mechanical filter designers and frequently in audio electronics also.  The mapping has the advantage of preserving network topologies across domains but does not maintain the mapping of impedances.  The Trent analogy classes the power conjugate variables as either ''across'' variables, or ''through'' variables depending on whether they act across an element of a system or through it.  This largely ends up the same as the mobility analogy except in the case of the fluid flow domain (including the acoustics domain).  Here pressure is made analogous to voltage (as in the impedance analogy) instead of current (as in the mobility analogy).  However, force in the mechanical domain ''is'' analogous to current because force acts ''through'' an object.<ref>Busch-Vishniac, pp. 19–20</ref>

There are some commonly used analogies that do not use power conjugate pairs.  For sensors, correctly modelling energy flows may not be so important.  Sensors often extract only tiny amounts of energy into the system.  Choosing variables that are convenient to measure, particularly ones that the sensor is sensing, may be more useful.  For instance, in the [[Physics:Thermal resistance|thermal resistance]] analogy, thermal resistance is considered analogous to electrical resistance, resulting in temperature difference and thermal power mapping to voltage and current respectively.  The power conjugate of temperature difference is not thermal power, but rather [[Physics:Entropy|entropy]] flow rate, something that cannot be directly measured.  Another analogy of the same sort occurs in the magnetic domain.  This maps [[Physics:Magnetic reluctance|magnetic reluctance]] to electrical resistance, resulting in [[Physics:Magnetic flux|magnetic flux]] mapping to current instead of magnetic flux rate of change as required for compatible variables.<ref>Busch-Vishniac, pp. 18, 20</ref>

== History ==
The matrix representation of [[Linear algebra|linear algebra]]ic equations has been known for some time.  [[Biography:Henri Poincaré|Poincaré]] in 1907 was the first to describe a transducer as a pair of such equations relating electrical variables (voltage and current) to mechanical variables (force and velocity).  Wegel, in 1921, was the first to express these equations in terms of mechanical impedance as well as electrical impedance.<ref>Pierce, p. 200</ref>

The first use of transfer matrices to represent a MIMO control system was by Boksenbom and Hood in 1950, but only for the particular case of the gas turbine engines they were studying for the [[Organization:National Advisory Committee for Aeronautics|National Advisory Committee for Aeronautics]].<ref>{{multiref|Kavanagh, p. 350|Bokenham & Hood, p. 581}}</ref>  Cruickshank provided a firmer basis in 1955 but without complete generality.  Kavanagh in 1956 gave the first completely general treatment, establishing the matrix relationship between system and control and providing criteria for realisability of a control system that could deliver a prescribed behaviour of the system under control.<ref>Kavanagh, pp. 349–350</ref>

== See also ==
* [[Physics:Transfer-matrix method (optics)|Transfer-matrix method (optics)]]

==References==
{{reflist|23em}}

==Bibliography==
* Bessai, Horst, ''MIMO Signals and Systems'', Springer, 2006 {{isbn|038727457X}}.
* Bakshi, A.V.; Bakshi, U.A., ''Network Theory'', Technical Publications, 2008 {{isbn|8184314027}}.
* Boksenbom, Aaron S.; Hood, Richard, [https://ntrs.nasa.gov/search.jsp?R=19930092041 "General algebraic method applied to control analysis of complex engine types"], [[Organization:National Advisory Committee for Aeronautics|NACA]] Report 980, 1950.
* Busch-Vishniac, Ilene J., ''Electromechanical Sensors and Actuators'', Springer, 1999 {{isbn|038798495X}}.
* Chen, Wai Kai, ''The Electrical Engineering Handbook'', Academic Press, 2004 {{isbn|0080477488}}.
* Choma, John, ''Electrical Networks: Theory and Analysis'', Wiley, 1985 {{isbn|0471085286}}.
* Cruickshank, A. J. O., "Matrix formulation of control system equations", ''The Matrix and Tensor Quarterly'', vol. 5, no. 3, p. 76, 1955.
* Iyer, T. S. K. V., ''Circuit Theory'', Tata McGraw-Hill Education, 1985 {{isbn|0074516817}}.
* Kavanagh, R. J., [https://web.archive.org/web/20161226061512/http://booksc.org/book/2062381 "The application of matrix methods to multi-variable control systems"], ''Journal of the Franklin Institute'', vol. 262, iss. 5, pp. 349–367, November 1956.
* Koenig, Herman Edward; Blackwell, William A., ''Electromechanical System Theory'', McGraw-Hill, 1961 {{OCLC|564134}}
* Levine, William S., '' The Control Handbook'', CRC Press, 1996 {{isbn|0849385709}}.
* Nguyen, Cam, ''Radio-Frequency Integrated-Circuit Engineering'', Wiley, 2015 {{isbn|1118936485}}.
* Olsen A., [https://ieeexplore.ieee.org/document/1082571/ "Characterization of Transformers by h-Paraameters"], ''IEEE Transactions on Circuit Theory'', vol. 13, iss. 2, pp. 239–240, June 1966.
* Pierce, Allan D. ''Acoustics: an Introduction to its Physical Principles and Applications'', Acoustical Society of America, 1989 {{isbn|0883186128}}.
* Poincaré, H., [http://henripoincarepapers.univ-lorraine.fr/bibliohp/index.php?a=on&art=%C3%89tude+du+r%C3%A9cepteur+t%C3%A9l%C3%A9phonique&action=Chercher "Etude du récepteur téléphonique"], ''Eclairage Electrique'', vol. 50, pp. 221–372, 1907.
* Wegel, R. L., [https://ieeexplore.ieee.org/document/6594447/ "Theory of magneto-mechanical systems as applied to telephone receivers and similar structures"], ''Journal of the American Institute of Electrical Engineers'', vol. 40, pp. 791–802, 1921.
* Yang, Won Y.; Lee, Seung C., ''Circuit Systems with MATLAB and PSpice'', Wiley 2008, {{isbn|0470822406}}.

[[Category:Control theory]]
[[Category:Control engineering]]
[[Category:Systems engineering]]
[[Category:Systems theory]]
[[Category:Automation]]
[[Category:Signal processing]]
[[Category:Frequency-domain analysis]]
[[Category:Types of functions]]
[[Category:Matrices]]

{{Sourceattribution|Transfer function matrix}}