Finite wordlength properties of matrix inversion algorithms in fixed-point and logarithmic number systems

Matrix inversion is sensitive towards the number representation used. In this paper simulations of matrix inversion with numbers represented in the fixed-point and logarithmic number systems (LNS) are presented. A software framework has been implemented to allow extensive simulation of finite wordlength matrix inversion. Six different algorithms have been used and results on matrix condition number, wordlength, and to some extent matrix size are presented. The simulations among other things show that the wordlength requirements differ significantly between different algorithms in both fixed-point and LNS representations. The results can be used as a starting point for a matrix inversion hardware implementation.

[1]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[2]  Earl E. Swartzlander,et al.  The Sign/Logarithm Number System , 1975, IEEE Transactions on Computers.

[3]  George M. Papadourakis,et al.  High speed implementation of matrix inversion algorithms in orthogonal systolic architectures , 1988, Conference Proceedings '88., IEEE Southeastcon.

[4]  Nicholas J. Higham,et al.  INVERSE PROBLEMS NEWSLETTER , 1991 .

[5]  A. Happonen,et al.  Several approaches to fixed-point implementation of matrix inversion , 2005, International Symposium on Signals, Circuits and Systems, 2005. ISSCS 2005..

[6]  Johan Eilert,et al.  Efficient Complex Matrix Inversion for MIMO Software Defined Radio , 2007, 2007 IEEE International Symposium on Circuits and Systems.

[7]  Lei Ma,et al.  QR Decomposition-Based Matrix Inversion for High Performance Embedded MIMO Receivers , 2011, IEEE Transactions on Signal Processing.

[8]  Ryan Kastner,et al.  GUSTO: An automatic generation and optimization tool for matrix inversion architectures , 2010, TECS.