Numerically Robust Implementations of Fast Recursive Least Squares Adaptive Filters Using Interval Arithmetic
Numerically Robust Implementations of Fast Recursive Least Squares Adaptive Filters using Interval Arithmetic
Christopher Peter Callender, B .Sc. AMIEE.
A thesis submitted for the degree of to the Faculty of Science at the University of Edinburgh.
1991
-1- University of Edinburgh
Abstract of Thesis
Name of Candidate Christopher Peter Callender Address Degree Ph.D. Date March 25, 1991 Title of Thesis Numerically Robust Implementations of Fast Recursive Least Squares Adaptive Filters using Interval Arithmetic Number of words in the main text of Thesis Approximately 28,000
Algorithms have been developed which perform least squares adaptive filtering with great computational efficiency. Unfortunately, the fast recursive least squares (RLS) algorithms all exhibit numerical instability due to finite precision computational errors, resulting in their failure to produce a useful solution after a short number of iterations.
In this thesis, a new solution to this instability problem is considered, making use of interval arithmetic. By modifying the algorithm so that upper and lower bounds are placed on all quantities calculated, it is possible to obtain a measure of confidence in the solution calculated by a fast RLS algorithm and if it is subject to a high degree of inaccuracy due to finite precision computational errors, then the algo- rithm may be rescued, using a reinitialisation procedure.
Simulation results show that the stabilised algorithms offer an accuracy of solution comparable with the standard recursive least squares algorithm. Both floating and fixed point implementations of the interval arithmetic method are simulated and long-term stability is demonstrated in both cases.
A hardware verification of the simulation results is also performed, using a digital signal processor(DSP). The results from this indicate that, the stabilised fast 'RLS algorithms are suitable for a number of applications requiring high speed, real time adaptive filtering.
A design study for a very large scale integration (VLSI) technology coprocessor, which provides hardware support for interval multiplication, is also considered. This device would enable the hardware realisation of a fast RLS algorithm to operate at far greater speed than that obtained by performing interval multiplication using a DSP.
Finally, the results presented in this thesis are summarised and the achievements and limitations of the work are identified. Areas for further research are suggested. Acknowledgements
There are many people who deserve thanks for the assistance and support which have made the work of this thesis possible. I would like to thank all of the members of the signal processing group, both past and present, for their helpful dis- cussions and comments on my work. They have helped to clarify many of my ideas.
Special thanks must go to Professor Cohn Cowan. In his role as my supervisor, he has contributed greatly to the project and his guidance has been very much appreci- ated. Similarly, I am most grateful to Dr. Bernie Muigrew for initially being my second supervisor and for taking over when Professor Cowan left the department. Thanks must also go to Professor Peter Grant for his encouragement and advice during my time in the group.
Finally, I would like to make it clear how much I appreciate the support of my parents during the time spent on this work. Contents
Chapter 1:Introduction 2
1.1. Adaptive Filters - Structures and Applications 2
1.2. Families of Adaptive Algorithms ...... 5
1.3. Applications of Adaptive Filters ...... 15
1.3.1. Prediction ...... 17
1.3.2. Noise Cancellation ...... 17
1.3.3. System Identification ...... 18
1.3.4. Inverse Modelling ...... 19
1.4. Organisation of Thesis ...... 19
Chapter 2:Least Squares Algorithms for Adaptive Filtering 22
2.1. Introduction ...... 22
2.2. The Least Squares Problem for Linear Transversal Adaptive Filtering ...... 23
2.3. The Conventional Recursive Least Squares Algorithm 28
2.4. Data Windows ...... 29
2.5. Computational Complexity ...... 31
2.6. The Fast Kalman Algorithm ...... 32
2.7. The Fast A Posteriori Error Sequential Technique 39
2.8. The Fast Transversal Filters Algorithm ...... 41
2.9. Comparison of the Least Squares Algorithms 44
-v - 2.10. Numerical Instability 45
2.10.1. Normalised Algorithms ...... 47
2.10.2. Lattice Algorithms ...... 48
2.10.3. Stabilisation by Regular Reinitialisation 49
2.10.4. Error Feedback ...... 50
2.1.1. Conclusions ...... 51
Chapter 3:Interval Arithmetic ...... 53
3.1. Introduction ...... 53
3.2. Interval Numbers