The branch method of circuit analysis is the basis for our proof. Let’s assume that inside of a box we have a circuit with two leads exiting the box (
Figure 2.6.3). The circuit in the box is comprised of a complicated network of resistors
\(R\text{,}\) \(\ell\) DC voltage sources
\(V\text{,}\) and
\(k\) ideal current sources
\(I\text{.}\) We will also assume
\(j\) unknown currents
\(i\) that are dependent on the size of the voltage difference
\(\Delta v_\text{ext}\) across the output terminals of the circuit.
Assume that our analysis results in \(N_1\) independent KCL equations. Each KCL equation can be expressed as
\begin{equation}
a_{r,\text{ext}} i_{\text{ext}} + \sum\limits_{n=1}^j a_{rn} i_n + \sum\limits_{n=1}^k b_{rn} I_n=0\tag{2.6.3}
\end{equation}
where \(a_{rn}\) and \(b_{rn}\) have values of +1, -1, or 0 depending on whether the currents \(i_n\) and \(I_n\) are entering, exiting, or not present at a junction and \(r\) is an index counting over KCL equations.
Analysis will also provide \(N_2\) independent KVL equations, analyzing loops that are absent known current sources \(I_n\text{.}\) This restriction will not cost us anything because this proof relies on knowledge of the currents through each circuit segment. Each KVL equation has the form
\begin{equation}
m_{r,\text{ext}} \Delta v_\text{ext} + \sum\limits_{n=1}^\ell m_{rn} V_n + c_{r,\text{ext}} i_\text{ext}
+ \sum\limits_{n=1}^j c_{rn} i_n = 0\text{.}\tag{2.6.4}
\end{equation}
The variables \(m_{rn}\) have values of either +1, or -1 depending on the battery orientation relative to the direction traveled around the loop, or a value of 0 for batteries not present in a given loop. Variable \(c_{r_n}\) has units of resistance because each nonzero term in these sums represents a voltage change described by Ohm’s Law. Here, the index \(r\) counts over KVL equations.
These
\(N_1\) instances of
(2.6.3) and
\(N_2\) instances of
(2.6.4) can be re-expressed as a matrix equation
\(\mathbf{A}\mathbf{x}=\mathbf{v}\text{.}\) First, let’s rearrange these equations so that all unknowns are on the left side of the equation and all of the source terms are on the right
\begin{align*}
a_{r,\text{ext}} i_{\text{ext}} + \sum\limits_{n=1}^j a_{rn} i_n \amp = - \sum\limits_{n=1}^k b_{rn} I_n\\
c_{r,\text{ext}} i_\text{ext} + \sum\limits_{n=1}^j c_{rn} i_n
\amp = - m_{r,\text{ext}} \Delta v_\text{ext} - \sum\limits_{n=1}^\ell m_{rn} V_n \text{.}
\end{align*}
We can express these equations in the form \(\mathbf{A}\mathbf{x}=\mathbf{v}\) if we let
\begin{equation*}
\mathbf{A} = \begin{pmatrix}
a_{1,\text{ext}} \amp a_{11} \amp \cdots \amp a_{1j}\\
\vdots \amp \vdots \amp \ddots \amp \vdots \\
a_{N_1,\text{ext}} \amp a_{N_1 1}\amp \cdots \amp a_{N_1 j} \\
c_{1,\text{ext}} \amp c_{11} \amp \cdots \amp c_{1j} \\
\vdots \amp \vdots \amp \ddots \amp \vdots \\
c_{N_2,\text{ext}} \amp c_{N_2 1} \amp \cdots \amp c_{N_2 j} \\
\end{pmatrix}\text{,}
\end{equation*}
\begin{equation*}
\mathbf{x} =
\begin{pmatrix}
i_\text{ext} \\ i_1 \\ \vdots \\ i_j \\
\end{pmatrix}
\quad , \quad
\mathbf{v} =
\begin{pmatrix}
-\sum\limits_{n=1}^{k} b_{1n} I_n \\
\vdots \\
-\sum\limits_{n=1}^{k} b_{N_1 n} I_n \\
-m_{1,\text{ext}}\Delta v_\text{ext} - \sum\limits_{n=1}^\ell m_{1n} V_n\\
\vdots \\
-m_{N_{2},\text{ext}}\Delta v_\text{ext} - \sum\limits_{n=1}^\ell m_{N_2 n} V_n
\end{pmatrix} \text{.}
\end{equation*}
Now, every element in matrix \(\mathbf{A}\) is either a unitless number or a function of resistances only. More to the point, there is no voltage or current dependence in \(\mathbf{A}\text{.}\) If we multiply both sides of our matrix equation by \(\mathbf{A}^{-1}\text{,}\) we get
\begin{equation*}
\mathbf{x}=\mathbf{A}^{-1}\mathbf{v}\text{.}
\end{equation*}
The solution that is most important to us is \(i_\text{ext}\) as we wish to examine the relationship between output current and voltage. If we actually perform the calculations necessary to invert \(\mathbf{A}\text{,}\) use matrix multiplication to find an expression for \(i_\text{ext}\text{,}\) and then group source terms and rename constants multiplying each source (\(v_\text{ext},V,\ I\)), the resulting solution is of the form
\begin{equation*}
i_\text{ext} = \alpha_0 v_\text{ext} + \sum\limits_{n=1}^\ell \alpha_n V_n + \sum\limits_{n=1}^k \beta_n I_n
\end{equation*}
where all \(\alpha_n\) and \(\beta_n\) terms are functions of resistances only and are thus constants. This result demonstrates that \(i_\text{ext}\) is a linear function of \(v_\text{ext}\) with slope \(\alpha_0\) and a y-intercept that equals the two summations on the right side of the equation. So, we have now shown that any circuit comprised of solely linear circuit elements is itself a linear circuit.