Adaptive Predictors for Extracting Physiological Signals in Two
ADAPTIVE PREDICTORS FOR EXTRACTING PHYSIOLOGICAL SIGNALS IN TWO
MODERN BIOINSTRUMENTS
by
Brently W. Robinson
APPROVED BY SUPERVISORY COMMITTEE:
Dr. Mohammad Saquib
Dr. Issa Panahi
Dr. Naofal Al-Dhahir
Dr. P.K. Rajasekaran Copyright 2019
Brently W. Robinson
All rights reserved To All My Family. ADAPTIVE PREDICTORS FOR EXTRACTING PHYSIOLOGICAL SIGNALS IN TWO
MODERN BIOINSTRUMENTS
by
BRENTLY W. ROBINSON, BS, MS
DISSERTATION
Presented to the Faculty of
The University of Texas at Dallas
in Partial Fulfillment
of the Requirements
for the Degree of
DOCTOR OF PHILOSOPHY IN
ELECTRICAL ENGINEERING
THE UNIVERSITY OF TEXAS AT DALLAS
August 2019 ACKNOWLEDGMENTS
I would like to express deep gratitude and appreciation to Dr. Mohammad Saquib, my PhD advisor, who consistently and patiently encouraged me through periods of personal challenges and academic development. Apart from his support and guidance this research would not have come to fruition.
I would also like to acknowledge and thank, my supervisory committee members, Dr. Issa Panahi, Dr. Naofal Al-Dhahir, and Dr. P.K. Rajasekaran. They provided feedback and direction through their comments and suggestions that helped to develop my research into a more comprehensive work.
Finally, I would like to acknowledge and remember Dr. Phillip Loizou, an original supervisory committee member, whose guidance brought me to UTD and who was always wiling to share both personal and academic advice.
My research was, in part, funded by the Raytheon Educational Assistance and Advanced Studies Program.
July 2019
v ADAPTIVE PREDICTORS FOR EXTRACTING PHYSIOLOGICAL SIGNALS IN TWO
MODERN BIOINSTRUMENTS
Brently W. Robinson, PhD The University of Texas at Dallas, 2019
Supervising Professor: Dr. Mohammad Saquib
Physiological signals are at the core of understanding, diagnosing, and treating the human body. Those provide valuable insight into the internal function and state of systems within the anatomy. Depending on the system being observed, a physiological signal can either be a source of information or a source of interference.
This dissertation first examines physiological hand tremor, the body’s response to stress, tiredness, or hunger, as a source of interference during microsurgery. It then examines the electrocardiogram (ECG), a source of vital information about the condition of the heart, when corrupted by broadband interference. Examination of these two physiological signals, obtained by bioinstruments, leads us to develop novel real-time adaptive predictors. Based on Kalman adaptation principle, an adaptive predictor is developed for removing physio- logical hand tremor and a scalable, cascaded predictor is designed for removing broadband interference from the ECG.
Due to the real-time requirement of bioinstruments, this dissertation addresses the issues with implementing adaptive algorithms in fixed-point representation. A proposed modified binary floating-point format is presented and is shown to overcome the prior known issues associated with fixed-point implementations and demonstrated for removal of physiological hand tremor.
vi TABLE OF CONTENTS
ACKNOWLEDGMENTS ...... v ABSTRACT ...... vi LIST OF FIGURES ...... ix LIST OF TABLES ...... xi CHAPTER 1 INTRODUCTION ...... 1 1.1 Overview ...... 1 1.2 Background ...... 2 1.2.1 Cancelling Broadband Interference ...... 3 1.2.2 Cancelling Narrowband Interference ...... 6 1.3 Motivation ...... 7 1.4 Outline ...... 9 1.5 Contributions ...... 9 CHAPTER 2 ADAPTIVE PREDICTORS ...... 11 2.1 Introduction ...... 11 2.2 Wiener Filter ...... 11 2.3 Adaptive Filter ...... 14 2.3.1 Cost Function ...... 16 2.3.2 Minimization Algorithm ...... 17 2.4 Adaptive Predictor ...... 19 2.5 Conclusion ...... 25 CHAPTER 3 PREDICTING PHYSIOLOGICAL HAND TREMOR DURING MICRO- SURGERY ...... 26 3.1 Introduction ...... 26 3.2 System Model ...... 28 3.3 Proposed Algorithm ...... 29 3.4 Performance ...... 33 3.5 Hardware Implementation ...... 33 3.5.1 Fixed-Point Representation ...... 34
vii 3.5.2 Floating-Point Representation ...... 40 3.5.3 Comparison ...... 48 3.6 Conclusion ...... 50 CHAPTER 4 ENHANCEMENT OF THE ELECTROCARDIOGRAM USING ADAP- TIVE COHERENT AVERAGING ...... 51 4.1 Introduction ...... 51 4.2 Proposed Structure ...... 55 4.3 Cascaded Architecture ...... 58 4.4 Analysis ...... 61 4.4.1 Derivation of Optimum Filter at Stage 1 ...... 61 4.4.2 Derivation of Optimum Filter at Stage n ...... 62 4.4.3 Derivation of MSE at Stage n ...... 64 4.5 Performance ...... 66 4.5.1 Least Mean Squares ...... 67 4.5.2 Recursive Least Squares ...... 67 4.5.3 Numerical Results ...... 68 4.5.4 Synthetic ECG Results ...... 69 4.6 Conclusion ...... 70 CHAPTER 5 CONCLUSIONS ...... 80 5.1 Summary ...... 80 5.2 Future Works ...... 81 REFERENCES ...... 82 BIOGRAPHICAL SKETCH ...... 87 CURRICULUM VITAE
viii LIST OF FIGURES
1.1 The adaptive noise canceller [55]...... 2 1.2 The implementation of adaptive cancellation for electrosurgical interference [57]. 3 1.3 The filter model for the cancellation of broadband interference [55]...... 4 1.4 The time-sequenced adaptive filter [13]...... 5 1.5 The adaptive recurrent filter [47]...... 5 1.6 The filter model for the cancellation of narrowband interference [55]...... 6 1.7 The Fourier linear combiner [51]...... 7 1.8 The weighted-frequency Fourier linear combiner [35]...... 8 2.1 The FIR transversal filter structure...... 13 2.2 The Wiener filter structure...... 13 2.3 The adaptive filter structure...... 14 2.4 The adaptive predictor structure...... 20 2.5 The adaptive predictor for cancelling broadband interference...... 22 2.6 The adaptive predictor for cancelling narrowband interference...... 23 2.7 The adaptive predictor with signal leakage...... 24 3.1 The hand tremor measurement and cancellation system model...... 28 3.2 The adaptive predictor used for hand tremor cancellation...... 29 3.3 Example of explosive divergence and the stalling effect...... 35 3.4 The mean square error of 64 bit fixed-point RLS...... 37 3.5 The mean square error of 64 bit fixed-point Kalman...... 38 3.6 The mean square error of 32 bit fixed-point RLS...... 38 3.7 The mean square error of 32 bit fixed-point Kalman...... 39 3.8 The mean square error of 16 bit fixed-point RLS...... 39 3.9 The mean square error of 16 bit fixed-point Kalman...... 40 3.10 The 16 bit half precision binary floating-point format...... 41 3.11 The binary floating-point unpack process...... 41 3.12 The mean square error of 64 bit floating-point RLS...... 44 3.13 The mean square error of 64 bit floating-point Kalman...... 44
ix 3.14 The mean square error of 32 bit floating-point RLS...... 45 3.15 The mean square error of 32 bit floating-point Kalman...... 45 3.16 The mean square error of 16 bit floating-point RLS...... 46 3.17 The mean square error of 16 bit floating-point Kalman...... 46 3.18 The 24 bit extended mantissa half precision binary floating-point format. . . . . 47 3.19 The mean square error of 24 bit floating-point RLS...... 47 3.20 The mean square error of 24 bit fixed-point Kalman...... 48 4.1 The synthetic ECG signal model...... 52 4.2 The frequency spectrum of the synthetic ECG signal...... 53 4.3 The reoccurence interval of the synthetic ECG signal model...... 54 4.4 The adaptive predictor structure...... 55 4.5 The adaptive coherent averaging structure...... 56 4.6 The cascade form of the ACA structure...... 58 4.7 The mean square error for the ACA Stage 1 (SNR = −5dB)...... 72 4.8 The mean square error for the ACA Stage 2 (SNR = −5dB)...... 72 4.9 The mean square error for the ACA Stage 3 (SNR = −5dB)...... 73 4.10 The mean square error for the ACA Stage 4 (SNR = −5dB)...... 73 4.11 The mean square error for the ACA Stage 5 (SNR = −5dB)...... 74 4.12 The mean square error for the ACA Stage 6 (SNR = −5dB)...... 74 4.13 The mean square error for the ACA Stage 7 (SNR = −5dB)...... 75 4.14 The mean square error for the ACA Stage 8 (SNR = −5dB)...... 75 4.15 The mean square error for the ACA Stage 1 (SNR = 5dB)...... 76 4.16 The mean square error for the ACA Stage 2 (SNR = 5dB)...... 76 4.17 The mean square error for the ACA Stage 3 (SNR = 5dB)...... 77 4.18 The mean square error for the ACA Stage 4 (SNR = 5dB)...... 77 4.19 The mean square error for the ACA Stage 5 (SNR = 5dB)...... 78 4.20 The mean square error for the ACA Stage 6 (SNR = 5dB)...... 78 4.21 The mean square error for the ACA Stage 7 (SNR = 5dB)...... 79 4.22 The mean square error for the ACA Stage 8 (SNR = 5dB)...... 79
x LIST OF TABLES
3.1 The comparison between RLS and derived Kalman algorithms...... 32 3.2 The RLS and Kalman simulation results for hand tremor cancellation...... 33 3.3 The 3-stage binary floating-point multiplication...... 42 3.4 The 3-stage binary floating-point divide...... 42 3.5 The 3-stage binary floating-point addition/subtraction...... 43 3.6 The fixed-point RLS results for hand tremor cancellation...... 48 3.7 The fixed-point Kalman results for hand tremor cancellation...... 49 3.8 The floating-point RLS results for hand tremor cancellation...... 49 3.9 The floating-point Kalman results for hand tremor cancellation...... 50 4.1 The output SNDR for different input SNR of synthetic ECG and AWGN. . . . . 70 4.2 The output SNDR for different input SNR of synthetic ECG and AWGN. . . . . 70
xi CHAPTER 1
INTRODUCTION
1.1 Overview
Physiological signals are generated within the human body and represent the internal func- tion and state of the eleven major organ systems. These signals are transmitted and measured in various forms, including: electrical, acoustic, magnetic, optical, mechanical, and chemi- cal. A physiological signal can represent a single independent response or a complex set of time dependent responses. Since these signals originate within the body but are typically measured outside of the body, they are subject to both internal interference and external measurement noise. Because each of the major organ systems generates physiological signals within the body, each signal can be considered either a source of information or a source of interference with respect to determining the state of a particular system.
Biointruments are used to measure a physiological signal from the body. The measure- ment is either processed in real-time by the bioinstrument or offline by a general purpose processor in order to extract the desired signal. Adaptive signal processing has made sig- nificant contributions toward offline processing of physiological signals, but contributions toward real-time processing has been slow. This is primarily due to the barriers present when implementing adaptive algorithms in real-time specialized processors within a bioin- strument.
In this chapter we first look at the origin of the adaptive filter and limitations of its implementation. We then discuss how the adaptive filter is used to address two different classes of noise cancellation. The first class is the cancellation of broadband interference in the presence of a narrowband signal. The second class is the cancellation of narrowband interference in the presence of a broadband signal. In each of these classes we show exam- ples of how variants of the adaptive filter have been proposed to extend and specialize the
1 general adaptive filter. We then consider cancellation of broadband noise in the presence of a quasi-periodic signal. This quasi-periodicity motivates a novel variant of the adaptive
filter proposed in this dissertation.
1.2 Background
Research on the enhancement of physiological signals gained significant momentum with the introduction of the adaptive filter, shown in Fig. 1.1, by Widrow et al. in 1975. In this article they demonstrated multiple applications of the adaptive filter, but emphasized its use toward the cancellation of noise for two physiological signals, the electrocardiogram
(ECG) and speech. The adaptive filter is shown to obtain the equivalent performance as the Wiener filter for stationary signals and close to optimal performance for nonstationary signals. The adaptive filter has two significant advantages over the Wiener filter. The first advantage is the adaptive filter does not require the same amount of a priori knowledge about the statistical characteristics of the noise. The second advantage is the adaptive filter significantly reduced the computational complexity by the introduction of recursion. The trade-off incurred for these two advantages is that the adaptive filter did not perform as well as the Wiener filter on its first iteration, but required a number of iterations (time) in order to converge to the optimum weights obtained by the Wiener filter [55].
Figure 1.1. The adaptive noise canceller [55].
2 Over the next eight years, significant research was invested by Treichler [50] and Zeidler et al. [2, 3, 33, 58] into the analysis of the adaptive filter parameters required to minimize the time to convergence to the optimum weights obtained by the Wiener filter.
Following the introduction and analysis of the adaptive filter was the implementation into a bioinstrument in 1983 by Yelderman et al. [57]. In this article they presented a software implementation of the adaptive filter using a minicomputer programmed in assembly language. This implementation, shown in Fig. 1.2, was used for the adaptive cancellation of electrosurgical interference. Due to the computational complexity and the limited dynamic range of fixed-point, adaptive filters remain implemented in software using floating-point and do not advance into real-time hardware for many decades.
Figure 1.2. The implementation of adaptive cancellation for electrosurgical interference [57].
1.2.1 Cancelling Broadband Interference
The first class of adaptive filters we consider are focused on cancelling broadband interference in the presence of a narrowband signal. The filter model for this class of adaptive filters is shown in Fig. 1.3. In this model, the output is obtained from the filter output. The filter will
3 iteratively adaptive to form a bandpass filter around the narrowband signal. Unfortunately, broadband interference within the narrowband spectrum will not be eliminated. Within this class of adaptive filters multiple variants have been proposed. We first consider a variant that uses a delayed reference input and then consider a variant using a synthetic reference input.
Figure 1.3. The filter model for the cancellation of broadband interference [55].
In 1981, the time-sequenced adaptive filter, shown in Fig. 1.4, was proposed by Ferrara and Widrow [13], which focused on optimum estimation of nonstationary signals with recur- ring statistical character. This adaptive filter variant specifically considered applications of the ECG because of the quasi-periodic nature of the signal. The time-sequenced adaptive fil- ter increased performance by introducing an array of transversal filters, which also increased computational complexity. The introduction of the array required an input sequence num- ber, which was expected to be obtained by a reference signal. The trade-off for performance was both computational complexity and additional knowledge about the signal of interest in order to provide an accurate sequence number. In 1991, the adaptive recurrent filter, shown in Fig. 1.5, was proposed by Thakor and Zhu [47], which reduced the computational complexity of the time-sequence adaptive filter by significantly reducing the number of transversal filters. The adaptive recurrent filter still required a precise reference input to denote the beginning of a sequence. Both the
4 Figure 1.4. The time-sequenced adaptive filter [13]. time-sequenced adaptive filter and the recurrent filter suffered from the same limitation and require an reference input similar to that obtained by a matched filter in order to improve upon the adaptive filter.
Figure 1.5. The adaptive recurrent filter [47].
5 1.2.2 Cancelling Narrowband Interference
The second class of adaptive filters we consider focused on cancelling narrowband interfer- ence in the presence of a broadband signal. The filter model for this class of adaptive filters is shown in Fig. 1.6. In this model, the output is obtained by subtracting the filter output from the primary input. The filter will iteratively adapt to form a bandpass filter around the narrowband interference, which is then used to subtract the narrowband interference from the broadband signal. Unlike the first class where residual broadband interference per- sisted, the second class is subject to distortion of the broadband signal when the narrowband interference has overlapping signal spectrum. Similar to the prior class, multiple adaptive filters variants have been proposed, each slightly improving on the ability to estimate the narrowband interference.
Figure 1.6. The filter model for the cancellation of narrowband interference [55].
In 1994, the adaptive Fourier linear combiner (FLC), shown in Fig. 1.7, was proposed by Vaz and Thakor [51], which adaptively adjusted the amplitudes of sinusoidal compo- nents to maximally represent the signal to be removed. The FLC increased performance by synthetically estimating the narrowband interference through the combining of multiple si- nusoids, however, it was limited by predetermined frequencies of the sinusoidal components. The trade-off for performance was both computational complexity and additional knowledge about the frequency components of the narrowband interference.
6 Figure 1.7. The Fourier linear combiner [51].
In 1996, the weighted-frequency Fourier linear combiner (WFLC), shown in Fig. 1.8, was proposed by Riviere and Thakor [35], which improved on the prior FLC by adaptively ad- justing both the amplitudes and phases of the sinusoidal components to maximally represent the signal to be removed. The WFLC is not limited like the FLC requiring knowledge of the frequency components of the narrowband interference, however, this is a trade-off obtained at the cost of increased computational complexity and latency.
1.3 Motivation
From the background review of adaptive filters, we observe that the adaptive filter output is either subject to residual interference within the signal spectrum or signal distortion. We also note that adaptive filters are constrained to software implementations in floating- point due to the dynamic range required for the adaptive parameters and the computational complexity of each variant.
7 Figure 1.8. The weighted-frequency Fourier linear combiner [35].
This dissertation first examines physiological hand tremor, as a source of interference during microsurgery. Using the adaptive filter for cancellation of narrowband interference in the presence of a broadband signal, we develop a novel real-time adaptive predictor based on Kalman adaptation principle suitable for implementation into hardware.
We then demonstrate the issues with implementing the adaptive filter in fixed-point representation, which is a common hardware implementation. A proposed modified binary
floating-point format is presented and is shown to overcome the prior known issues associated with fixed-point implementations. This allows the adaptive filter to achieve the real-time requirement of modern bioinstruments.
Finally we examine the ECG, a source of vital information about the condition of the heart, when corrupted by broadband noise. Using the adaptive filter for cancellation of broadband interference in the presence of a periodic signal, we develop a scalable, cascaded structure based on the adaptive predictor capable of removing broadband interference from the ECG while preserving the morphological features of the signal.
8 1.4 Outline
The remainder of this dissertation is organized in the following manner.
In Chapter 2, we give an overview of the adaptive predictor, providing a foundation for this dissertation.
In Chapter 3, we investigate physiological hand tremor as a source of interference limiting precision during microsurgery. Using an adaptive linear prediction system model, we develop a modified Kalman adaptation algorithm for the prediction and removal of physiological hand tremor. The derived Kalman algorithm is compared with the existing recursive least squares
(RLS) adaptation algorithm. Additionally, we examine the barriers and tradeoffs in real-time implementations of adaptive predictors. A modified binary floating-point format is proposed that efficiently uses native hardware multipliers and overcome the issues. Analysis of both
fixed and floating-point precision is performed.
In Chapter 4, we examine the ECG when corrupted by broadband interference. Ex- panding on the adaptive predictor structure presented in Chapter 3, we develop an adaptive coherent averaging (ACA) structure. The ACA structure is capable of enhancing a quasi- periodic signal in the presence of broadband interference while preserving the morphological features of the signal. Analysis is performed to derive the optimum filter and the mean square error (MSE) for each stage of the cascaded ACA structure. This analysis demon- strates the cascading benefit over the existing adaptive predictor structure. We then present the performance of the ACA structure using two common adaptation algorithms: the least mean squares (LMS) and the RLS.
Finally, in Chapter 5 we present the conclusion and future work.
1.5 Contributions
Some of the results of this dissertation have been presented previously.
9 • Chapter 3 was presented at the 2010 IEEE International Conference on Acoustics Speech and Signal Processing (ICASSP) and also published in IEEE Journal of Selected Topics in Signal Processing (JSTSP) in June 2010.
• Chapter 4 was presented at the 2014 International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) and will be submitted to the IEEE Journal of Engineering in Medicine and Biology Society (EMBS) in Fall of 2019.
10 CHAPTER 2
ADAPTIVE PREDICTORS
2.1 Introduction
In this chapter we review the foundations of signal processing leading up to the adaptive
predictor. We begin by introducing the Wiener filter, which preceded the adaptive filter. The
Wiener filter obtains the optimum solution in the mean square sense, but is computationally
complex and requires noise characteristics that are not always available. This motivated the
development of the adaptive filter, which overcame many of the prior issues by introducing
a recursive adaptation algorithm. Finally we review the adaptive predictor, a specialized
version of the adaptive filter that only requires a single input.
2.2 Wiener Filter
One focus of signal processing is the use of a priori information about the statistics of a input signal in order to provide an estimated or enhanced output signal. In 1931, Wiener and Hopf developed the set of equations for optimum estimation, now commonly referred to as the Wiener-Hopf equations. These equations describes the relationship between a primary input vector d(k), of length p, given by
> d(k) = d(k), d(k − 1), . . . , d(k − p), d(k − p + 1) , (2.1)
and a reference input vector x(k), of length p, given by
> x(k) = x(k), x(k − 1), . . . , x(k − p), x(k − p + 1) . (2.2)
Pxd, is the cross-correlation vector between the reference input vector x(k), and the primary input d(k) denoted by
11 Pxd = E[x(k)d(k)]. (2.3)
Rxx, is the autocorrelation matrix of the reference input vector x(k) denoted by
rxx(0) rxx(1) rxx(2) . . . rxx(p − 1) rxx(1) rxx(2) rxx(3) . . . rxx(p − 2) . . . . Rxx = r (2) ...... , (2.4) xx ...... . . . . . rxx(p − 1) rxx(p − 2) ...... rxx(0) where
rxx(j) = E[x(k)x(k − j)]. (2.5)
In 1947, Levinson described the Wiener-Hopf equation for discrete time in its more com- mon matrix form [17] is written as
Rxxwopt = Pxd. (2.6)
The Wiener-Hopf equation can provide the optimum coefficients wopt, in the mean square
−1 sense by multiplying both sides of (2.6) by Rxx . This form of the equation, commonly referred to as the Wiener filter solution is expressed as
−1 wopt = Rxx Pxd, (2.7) where the optimum coefficients wopt, are given by
> wopt = w1, w2, . . . , wp−1, wp . (2.8)
12 x(k) x(k-1) x(k-p) x(k-p+1) z-1 z-1 … z-1
w1 x w2 x wp-1 x wp x y(k) + … + +
Figure 2.1. The FIR transversal filter structure.
Using the optimum coefficients wopt, in a Finite Impulse Response (FIR) transversal filter, shown in Fig. 2.1, to process the reference input vector x(k), we can describe the Wiener filter output y(k) as
p X y(k) = wix(k − i), (2.9) i=1 which can also be written in vector form as
> y(k) = woptx(k). (2.10)
The Wiener filter structure, shown in Fig. 2.2, shows the relationship between the primary input d(k), the reference input x(k), and the Wiener filter output y(k).