Neumann series explained

T . This has the generator form\sum ^{inf ty} _{k =0} T ^{k} where T ^{k} is the k-times repeated application of T ; T ^{0} is the identity operator I and T ^{k} : = {}T ^{k-1} \circ{T} for k > 0 . This is a special case of the generalization of a geometric series of real or complex numbers to a geometric series of operators. The generalized initial term of the series is the identity operator T ^{0} = I and the generalized common ratio of the series is the operator T.

The series is named after the mathematician Carl Neumann, who used it in 1877 in the context of potential theory. The Neumann series is used in functional analysis. It is closely connected to the resolvent formalism for studying the spectrum of bounded operators and, applied from the left to a function, it forms the Liouville-Neumann series that formally solves Fredholm integral equations.

Properties

Suppose that

T is a bounded linear operator on the normed vector space X . If the Neumann series converges in the operator norm, then I - T is invertible and its inverse is the series:(I - T) ^{-1} = \sum ^{inf ty} _{k =0} T ^{k} ,where I is the identity operator in X . To see why, consider the partial sumsS _{n} : = \sum ^{n} _{k =0} T ^{k} .Then we have\lim _{n} (I - T)S _{n} = \lim _{n} \left(\sum ^{n} _{k =0} T ^{k} - \sum ^{n} _{k =0} T ^{k+1} \right) = \lim _{n} \left(I - T ^{n+1} \right) = I This result on operators is analogous to geometric series in R .

One case in which convergence is guaranteed is when

X is a Banach space and |T| < 1 in the operator norm; another compatible case is that $\backslash sum\_^\backslash infty\; |T^k|$ converges. However, there are also results which give weaker conditions under which the series converges.

Example

Let

C\in \R ^{3} be given by:\begin{pmatrix} 0 & (1)/(2) & (1)/(\\ 4) (5)/(7) & 0 & (1)/(\\ 7) (3)/(10) & (3)/(5) & 0 \end{pmatrix}. For the Neumann series $\backslash sum\_^n\; C^k$ to converge to (I - C) ^{-1} as n goes to infinity, the matrix norm of C must be smaller than unity. This norm is\begin{align} ||C|| _{infty} & =max _{i} \sum _{j} |c _{ij} | = max \left\lbrace (3)/(4) , (6)/(7) , (9)/(10) \right\rbrace = (9)/(10) < 1 , \end{align} confirming that the Neumann series converges.

Approximate matrix inversion

A truncated Neumann series can be used for approximate matrix inversion. To approximate the inverse of an invertible matrix

A , consider that\begin{align} A ^{-1} & = (I - I + A) ^{-1} \\ & = (I - (I - A)) ^{-1} \\ & = (I - T) ^{-1} \end{align}

for

T = (I - A). Then, using the Neumann series identity that $\backslash sum\_^\backslash infty\; T^k\; =\; (1\; -\; T)^$ if the appropriate norm condition on T = (I - A) is satisfied, $A^\; =\; (I\; -\; (I\; -\; A))^\; =\; \backslash sum\_^\backslash infty\; (I\; -\; A)^k.$ Since these terms shrink with increasing k, given the conditions on the norm, then truncating the series at some finite n may give a practical approximation to the inverse matrix:A ^{-1} ≈ \sum ^{n} _{k =0} (I - A) ^{k.}

The set of invertible operators is open

A corollary is that the set of invertible operators between two Banach spaces

B and B' is open in the topology induced by the operator norm. Indeed, let S : B \to B' be an invertible operator and let T: B \to B' be another operator. If |S - T| < |S ^{-1} | ^{-1} , then T is also invertible. Since |I - S ^{-1} T| < 1 , the Neumann series $\backslash sum\_^\backslash infty\; (I-S^T)^k$ is convergent. Therefore, we haveT ^{-1} S = (I - (I - S ^{-1} T)) ^{-1} = \sum ^{inf ty(I} _{k =0} - S ^{-1} T) ^{k}

Taking the norms, we get

|T ^{-1} S| \le (1)/(1- | I - (S ^{-1} T)|)

The norm of

T ^{-1} can be bounded by|T ^{-1} | \le \tfrac{ 1 }{ 1- q} |S ^{-1} | where q = |S - T| |S ^{-1} |.

Applications

The Neumann series has been used for linear data detection in massive multiuser multiple-input multiple-output (MIMO) wireless systems. Using a truncated Neumann series avoids computation of an explicit matrix inverse, which reduces the complexity of linear data detection from cubic to square.^[1]

Another application is the theory of propagation graphs which takes advantage of Neumann series to derive closed form expressions for transfer functions.

References

• Book: Werner, Dirk. 2005. Funktionalanalysis . de. Springer Verlag. 3-540-43586-7.

Notes and References

1. Book: Wu. M.. Yin. B.. Vosoughi. A.. Studer. C.. Cavallaro. J. R.. Dick. C.. 2013 IEEE International Symposium on Circuits and Systems (ISCAS2013) . Approximate matrix inversion for high-throughput data detection in the large-scale MIMO uplink . May 2013. https://ieeexplore.ieee.org/document/6572301. 2155–2158. 10.1109/ISCAS.2013.6572301. 1911/75011. 978-1-4673-5762-3 . 389966 . free.