Next Article in Journal
A Context-Aware Edge-Based VANET Communication Scheme for ITS
Next Article in Special Issue
Correction: Screening and Biosensor-Based Approaches for Lung Cancer Detection. Sensors 2017, 17, 2420
Previous Article in Journal
Modeling and Characterization of Traffic Flows in Urban Environments
Previous Article in Special Issue
Towards an Online Seizure Advisory System—An Adaptive Seizure Prediction Framework Using Active Learning Heuristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast and Robust Non-Sparse Signal Recovery Algorithm for Wearable ECG Telemonitoring Using ADMM-Based Block Sparse Bayesian Learning

1
School of Computer Science and Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
2
School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(7), 2021; https://doi.org/10.3390/s18072021
Submission received: 28 April 2018 / Revised: 4 June 2018 / Accepted: 19 June 2018 / Published: 23 June 2018
(This article belongs to the Special Issue Sensors for Health Monitoring and Disease Diagnosis)

Abstract

:
Wearable telemonitoring of electrocardiogram (ECG) based on wireless body Area networks (WBAN) is a promising approach in next-generation patient-centric telecardiology solutions. In order to guarantee long-term effective operation of monitoring systems, the power consumption of the sensors must be strictly limited. Compressed sensing (CS) is an effective method to alleviate this problem. However, ECG signals in WBAN are usually non-sparse, and most traditional compressed sensing recovery algorithms have difficulty recovering non-sparse signals. In this paper, we proposed a fast and robust non-sparse signal recovery algorithm for wearable ECG telemonitoring. In the proposed algorithm, the alternating direction method of multipliers (ADMM) is used to accelerate the speed of block sparse Bayesian learning (BSBL) framework. We used the famous MIT-BIH Arrhythmia Database, MIT-BIH Long-Term ECG Database and ECG datasets collected in our practical wearable ECG telemonitoring system to verify the performance of the proposed algorithm. The experimental results show that the proposed algorithm can directly recover ECG signals with a satisfactory accuracy in a time domain without a dictionary matrix. Due to acceleration by ADMM, the proposed algorithm has a fast speed, and also it is robust for different ECG datasets. These results suggest that the proposed algorithm is very promising for wearable ECG telemonitoring.

1. Introduction

Wearable telemonitoring of electrocardiogram (ECG) via wireless body area networks (WBAN) is a very important topic in telemedicine, and many works focused on the devices, sensor networks and other subjects in ECG telemonitoring [1,2,3], but the most significant challenge of practical application of wearable remote ECG monitoring system is power consumption [4]. Some researchers have demonstrated that the main consumption of energy in WBAN is data acquisition and wireless data transmission [5,6]. In common WBAN-based ECG monitors, sensor nodes fall short of energy efficiency due to large data acquired from continuous monitoring and the energy wireless links [7,8], thus it is desirable to reduce the amount of data that need to be acquired and transmitted.
Compressed sensing (CS) [9,10,11] is an emerging signal sampling technique that requires a much lower sampling rate to recover the original signal compared with the traditional Nyquist sampling method. Figure 1 shows a typical CS-based wireless ECG telemonitoring system, where the ECG signals are firstly sampled via an analog-digital converter; then, these samples are compressed by a measurement matrix and transmitted to the remote terminal, and the ECG signals are recovered from the compressed data in the remote terminal. Many studies have focused on the application of compressed sensing in WBAN-based ECG monitoring [12,13,14,15], but it is still in its infancy [16].
Generally, the recorded ECG signals in wireless telemonitoring are non-sparse in the time domain and it is hard to find other sparse domains [17] for signal recovery, while compressed sensing algorithms usually assume that the signals are sparse in a time domain or some other transformed domains and the signals are recovered in the sparse domain. Most common compressed sensing recovery algorithms [18,19,20,21,22,23] usually can accurately recover sparse signals, but can not do much about non-sparse signals [24]. Therefore, accurate recovery of non-sparse signals is very important for wearable ECG telemonitoring.
In order to solve the problem of compressed sensing recovery on non-sparse signals, Zhang et al. [25,26] proposed a block sparse Bayesian learning (BSBL) framework-based BSBL-BO (bound optimization) algorithm to accurately recover non-sparse physiological signals. Exploiting the intra-block correlation of the signals, the BSBL framework has shown its outstanding performance on physiological signals. However, for the practical wearable wireless ECG monitoring system, BSBL-BO is too slow to recover the ECG signals in time. Zhang et al. [25] have proposed another BSBL-L1 algorithm that is very fast, but its recovery accuracy is much worse than BSBL-BO. Similarly, Liu et al. [27] proposed a fast algorithm based on the BSBL framework while the recovery accuracy is also unsatisfactory. Block sparse Bayesian learning algorithms usually can recover non-sparse directly by exploiting the intra-block correlation of signals such as BSBL-BO [25], but some other BSBL-based algorithms can not recover non-sparse signals accurately such as BSBL-L1 [25] and BSBL-FM [27]. In fact, most existing compressed sensing algorithms fail to meet the requirements of real-time ECG telemonitoring for speed and accuracy. Therefore, a fast algorithm that can accurately recover non-sparse signals is very necessary.
In this work, in order to recover non-sparse ECG signal in wireless ECG monitoring fast and accurately, a fast and robust non-sparse signal recovery algorithm termed BSBL-ADMM based on BSBL framework was proposed. The minimization of BSBL cost function was transformed to a group lasso optimization problem and alternating direction multiplier method (ADMM) [28] is used to solve the optimization problem because of that ADMM can make the proposed algorithm jumpily converged. The proposed BSBL-ADMM algorithm can directly recover the non-sparse ECG signals in the time domain quickly and satisfactory accuracy. We used the famous MIT-BIH Arrhythmia Database [29] and MIT-BIH Long-Term ECG Database [30] to demonstrate the outstanding performance of the proposed BSBl-ADMM algorithm. In addition, we also built a practical wearable ECG telemonitoring system based on digital compressed sensing and the proposed BSBL-ADMM algorithm is simulated on the ECG datasets collected in our practical wearable ECG telemonitoring system. The experimental results show that the proposed BSBL-ADMM algorithm not only has the ability to recover non-sparse ECG signals with fast speed and satisfactory recovery accuracy, but is also robust against different ECG datasets. Therefore, the peoposed BSBL-ADMM algorithm is promising for telemonitoring of ECG signals. The main contributions of this paper are as follows:
  • While the sparse domain of ECG signals and the corresponding dictionary matrix are hard to find, the proposed BSBL-ADMM algorithm has ability to directly recover the non-sparse ECG signal in a time domain without a dictionary matrix, instead of some CS algorithms that need to transform the signal to a sparse domain by a dictionary matrix. Therefore, the proposed BSBL-ADMM algorithm is suitable for non-sparse ECG signal recovery in wearable ECG telemonitoring.
  • The alternating direction multiplier method (ADMM) is introduced in BSBL framework to solve the optimization problem of cost function. ADMM can reach an approximate iteration result quickly in just a few iterations that make the proposed algorithm jumpily converged. Thus, the proposed BSBL-ADMM algorithm can meet the real-time requirement of the wearable ECG telemonitoring.
  • To ensure the generality of the proposed algorithm, a practical wearable ECG telemonitoring system based on digital compressed sensing is built to collect practical wearable ECG datasets. The proposed algorithm was simulated on three different datasets include our practical wearable ECG datasets and two MIT-BIH Databases. The experimental results demonstrate the outstanding performance and robustness of the proposed BSBL-ADMM algorithm. These advantages make the application of the proposed BSBL-ADMM algorithm in wearable ECG telemonitoring to be very promising.
The remainder of this paper is organized as follows: Section 2 states the block sparse Bayesian learning (BSBL) framework; Section 3 presents the proposed BSBL-ADMM algorithm; Section 4 describes datasets, experiment settings, and experimental results; conclusions are given in the last section.

2. BSBL Framework for ECG Recovery

To facilitate the study of algorithms, we firstly assume that signals have passed through the analog-digital converter (ADC), as in the “digital CS” paradigm [13].

2.1. Compressed Sensing Based on WBAN

Compressed sensing is a newly introduced data sampling and compression technique in wearable ECG telemonitoring based on wireless body area networks.
Most compressed sensing algorithms recover the signals in a sparse transformed domain. However, in WBAN, the signals are usually non-sparse and it is also hard to find a sparse transformed domain. The Block Sparse Bayesian Learning (BSBL) algorithms [25] can usually recover the non-sparse ECG signals in the time domain directly.
The original ECG signal is denoted by α ; then, the basic model of CS can be expressed as
s = Φ α ,
where s is the compressed data and Φ is the measurement matrix whose row number is smaller than column number; both s and Φ are known to CS algorithm for recovery. CS algorithms use the compressed data s and the measurement matrix Φ to recover the original ECG signal α .

2.2. Block Sparse Bayesian Learning

In applications, α generally has some structures, and a widely studied structure is the block structure. With this structure, α can be expressed as:
α = [ α 1 , , α d 1 α 1 T , , α d g 1 + 1 , , α d g α g T ] T .
It means that α has g blocks, and, among these blocks, only a few blocks are nonzero. Of course, if the signal is non-sparse, then the intra-block correlation will determine whether the signal can be recovered successfully. In BSBL [25], each block α i is modeled as a parameterized multivariate Gaussian distribution:
p ( α i ; γ i , B i ) = N ( α i ; 0 , γ i B i ) ,
where γ i is an unknown nonnegative parameter controlling the block-sparsity of α , and B i is an unknown positive definite matrix modeling the intra-block correlation structure of the block α i . The blocks are assumed to be mutually independent; then, the prior of α is p ( α ; { γ i , B i } ) = N ( 0 , Σ 0 ) , where Σ 0 = d i a g { γ 1 B 1 , , γ g B g } . Assume the measurement noise to be independent and satisfy a Gaussian distribution with zero mean and unknown variance λ , thus the Gaussian likelihood is
p ( s | α ; λ ) = N ( s ; Φ α , λ I ) .
In addition, the posterior of x is expressed as
p ( α | s ; { γ i , B i } i = 1 g , λ ) = N ( μ , Σ ) ,
where
μ = Σ 0 Φ T ( λ I + Φ Σ 0 Φ T ) 1 s ,
Σ = ( Σ 0 1 + 1 λ Φ T Φ ) 1 .
The parameters can be estimated by a Type II maximum likelihood procedure, which is equivalent to minimize the following cost function:
L ( { γ i , B i } i , λ ) = l o g | λ I + Φ Σ 0 Φ T | + s T ( λ I + Φ Σ 0 Φ T ) 1 s .
This cost function can be minimized by many optimization methods such as expectation maximization (EM), bound optimization (BO) [25] and fast marginal likelihood maximization (FM) [27]. The EM and BO algorithms can achieve a high accuracy on non-sparse signal recovery while the speed is slow. The FM algorithm is much faster than EM and BO, but the recovery accuracy of FM is not satisfactory.

3. Proposed ADMM-Based BSBL

In order to make a fast CS recovery algorithm based on the BSBL framework [25], the cost function will be transformed from a function of γ to a function of x. We introduce the alternating direction multiplier method (ADMM) to solve the minimization of the transformed cost function. The iteration of ADMM is the inner loop of the proposed algorithm, and it need not converge completely and an approximate iteration result can make the proposed algorithm jumpily converged. ADMM can get a approximate iteration result quickly in just a few iterations and the completely converged result has little influence to the final recovery accuracy, thus the main loop of the proposed algorithm can converge rapidly.

3.1. Transformation of the Cost Function

The BSBL cost function (8) can be transformed from γ -space to x-space using some identity and inequality. The parameters λ and B can be treated as fixed values first for convenience.
As presented in Zhang et al. [25], the minimization of cost function can be transformed to the following optimization
α k + 1 = arg min α s Φ α 2 2 + λ i σ i ( k ) α i T B i 1 α i ,
where σ i = 2 T r B i Φ i T λ I + Φ Σ 0 Φ T 1 Φ i 1 2 , and γ i = 2 α i T B i 1 α i σ i .
In addition, the learning rule of λ and B i can be derived using the EM method
λ = | | s Φ α | | 2 2 + T r ( Σ Φ T Φ ) M ,
B i = Σ i + μ i ( μ i ) T γ i ,
where M is the length of α , μ i is the corresponding i-th block in μ , and Σ i is the corresponding i-th principal diagonal block in Σ .
Then, in order to avoid overfitting, a strategy including parameter average and the first-order auto-regressive (AR) process is used to constrain each B i to be the same matrix B. The average of all B i is calculated as
B ¯ = 1 g i = 1 g B i ,
then the matrix B can be reconstructed by a Toeplitz matrix as
B = T o e p l i t z ( [ 1 , r , , r d 1 ] ) ,
where r is the AR coefficient empirically calculated by r = r 1 / r 0 , in which the r 0 and r 1 are the average of the elements along the main diagonal and main sub-diagonal of the matrix B ¯ .

3.2. Minimizing the Cost Function Using ADMM

As the minimization of BSBL cost function (8) was transformed to the optimization problem (9), it is easy to find that the problem (9) can be transformed into a L1/L2 regularization, which is called group lasso [31]. Let u i = σ i B i 1 / 2 α i , u = [ u 1 T , u 2 T , , u g T ] T , and H = Φ · d i a g { B 1 1 / 2 / σ 1 , , B g 1 / 2 / σ g } , then the problem (9) can be transformed to the following form:
u ( k + 1 ) = arg min u s H u 2 2 + λ i u i 2 , arg min u 1 2 s H u 2 2 + λ 0 i u i 2 ,
where λ 0 = 1 2 λ .
Furthermore, it can be seen that problem (14) has a similar form as ADMM framework [28], then let f ( u ) = 1 2 | | s H u | | 2 2 , g ( z ) = λ 0 i | | z i | | 2 , the problem (9) can be expressed as an ADMM form:
m i n i m i z e f ( u ) + g ( z ) , s u b j e c t t o u z = 0 .

3.3. Solving the ADMM

ADMM is an algorithm that solves convex optimization problems by breaking them into smaller pieces, each of which are then easier to handle. As in the method of multipliers, the form of the augmented Lagrangian is as the following:
L ρ ( u , z , b ) = f ( u ) + g ( z ) + b T ( u z ) + ( ρ / 2 ) | | u z | | 2 2 .
Then, the ADMM consists of the iterations
u j + 1 = arg min u L ρ ( u , z j , b j ) ,
z j + 1 = arg min z L ρ ( u j + 1 , z , b j ) ,
b j + 1 = b j + ρ ( u j + 1 z j + 1 ) ,
where ρ > 0 . In ADMM, the u and z are updated in an alternating fashion, which accounts for the term alternating direction.
The ADMM also can be written in a scaled form, which is more convenient. Defining the residual r = u z , then
b T ( u z ) + ( ρ / 2 ) | | u z | | 2 2 = b T r + ( ρ / 2 ) | | r | | 2 2 = ( ρ / 2 ) | | r + ( 1 / ρ ) b | | 2 2 ( 1 / 2 ρ ) | | b | | 2 2 = ( ρ / 2 ) | | r + y | | 2 2 ( ρ / 2 ) | | y | | 2 2 ,
where y = ( 1 / ρ ) b is the scaled dual variable. The ADMM can be expressed as:
u j + 1 = arg min u 1 2 | | s H u | | 2 2 + ( ρ / 2 ) | | u z j + y j | | 2 2 ,
z i j + 1 = arg min z i λ 0 i | | z i | | 2 + ( ρ / 2 ) | | u j + 1 z i + y j + 1 | | 2 2 ,
y j + 1 = y j + u j + 1 z j + 1 .
The step over u (21) is an unconstrained strongly convex quadratic minimization whose unique solution is given by a linear system as follows:
u j + 1 = ( H T H + ρ I ) 1 ( H T s + ρ ( z j y j ) ) .
The step over z i (22) has a simple closed-form solution as
z i j + 1 = S λ 0 / ρ ( u i j + 1 + y j ) ,
where the soft thresholding operator S is defined as
S k ( a ) a k , a > k , 0 , | a | < k , a + k , a < k .
The u, z and y are iteratively updated by (23)–(25) and until converging. Using ADMM to solve the problem (14), all the subproblems are easy to solve, and the algorithm does not need strict conditions to ensure the global convergence. In addition, in our problem, it need not converge completely because the result with complete convergence has little effect on the proposed algorithm. Once the u i is solved, the α i can be obtained by
α i = B i 1 / 2 σ i 1 u i .
Note that, through the approximate iteration result that was reached in just a few iterations, the ADMM makes the proposed algorithm jumpily converged, thus the proposed algorithm can get a fast recovery speed. The following Algorithm 1 gives the main procedure of the proposed algorithm, which we call the BSBL-ADMM algorithm.
Algorithm 1 Proposed BSBL-ADMM Algorithm
1:
Input: Φ , s
2:
Initialize: α = [ 1 , 1 , , 1 ] , B i = I , σ i = 1
3:
While not converge
4:
while not converge do
5:
 Calculate γ i = 2 α i T B i 1 α i σ i
6:
 Calculate λ = | | s Φ α | | 2 2 + T r ( Σ Φ T Φ ) M
7:
 Calculate B i = Σ i + μ i ( μ i ) T γ i
8:
 Calculate B = T o e p l i t z ( [ 1 , r , , r d 1 ] )
9:
 Calculate σ i = 2 T r B i Φ i T λ I + Φ Σ 0 Φ T 1 Φ i 1 2
10:
 Calculate H = Φ · d i a g { B 1 1 / 2 / σ 1 , , B g 1 / 2 / σ g }
11:
while not converge do
12:
  Calculate u j + 1 = ( H T H + ρ I ) 1 ( H T s + ρ ( z j y j ) )
13:
  Calculate z i j + 1 = S λ 0 / ρ ( u i j + 1 + y j )
14:
  Calculate y j + 1 = y j + u j + 1 z j + 1
15:
end while
16:
end while
17:
Calculate α i = B i 1 / 2 σ i 1 u i
18:
Output: α

4. Experiments and Results

4.1. Datasets

In order to validate the performance of the proposed BSBL-ADMM algorithm without loss of generality, we used three different datasets including the famous MIT-BIH Arrhythmia Database [29], MIT-BIH Long-Term ECG Database [30] and the practical wearable ECG datasets as the original signals to generate the compressed data.
(1)
MIT-BIH Arrhythmia Database: This is the most representative database for arrhythmia, and as such it has been used for most of the published research. It was also the first database available for analysis of the ECG signals and has been constantly refined over the years [29,32]. It includes 48 two-channel recordings at 360 samples per second with about 30 min.
(2)
MIT-BIH Long-Term ECG Database: This database contains seven long-term ECG recordings, and each recording lasts about 14 h to 22 h at the sampling rate of 128 Hz [30,33], and the long-term ECG recordings can verify the performance of the algorithms effectively.
(3)
Practical wearable ECG datasets: We built a practical wearable system for ECG collecting based on digital compressed sensing. Figure 2 and Figure 3 show the block diagram and the practical devices of the ECG collecting system, respectively. In this system, three-lead electrodes are used to obtain ECG signals via AD8232 module [34], the leads LA and RA are attached at the left and right chest respectively and the LL lead is attached at the left lower abdomen. After we obtain the analog ECG signals via AD8232, an analog-digital converter is used to sample the analog signals to digital signals at the sampling frequency 250 Hz [35], which is commonly used for ECG monitoring in body area networks. Then, a STM32 microcontroller is used to compress the digital ECG signals by a simple matrix-vector multiplication based on compressed sensing, these compressed data are transmitted to a computer via bluetooth, and the ECG signals are recovered from the compressed data on the computer. Finally, eight one-channel recordings that lasted about 30 min were collected using this system from eight different people.

4.2. Experimental Setup

In our experiments, the original signals were divided into several frames with the length N = 500 . To obtain the compressed measurements of the length M, a measurement matrix must be constructed. Here, we use a sparse binary matrix as the measurement matrix. Note that the sparse binary matrix has been widely used in CS-based telemonitoring for its efficiency in storage and matrix-vector multiplication. The compression ratio (CR) was defined as
C R = N M N .
In the experiments, we first constructed the measurement matrices with N = 500 and M = 100 , 150 , 200 , 250 , 300 , the corresponding C R = 80 % , 70 % , 60 % , 50 % , 40 % . Every measurement matrix was generated with M rows and N columns. Regardless of the size, its each column contained 12 entries of 1 s with random locations and other entries are all 0 s. Then, the prepared frame vector x was multiplied by the measurement matrix and the measurement vector was obtained. Then, the signal would be recovered by the recovery algorithm using the measurement matrix and measurement vector.

4.3. Performance Metrics

To measure recovery quality, we use the percentage root-mean squared distortion (PRD) to quantify the error percentage between the original signal α and the recovered signal α ^ :
P R D = | | α ^ α | | 2 | | α | | 2 × 100 .
The lower PRD means the better recovery performance. In addition, Pearson correlation between original signal and recovered signal is another performance index that reflects the reconstruction quality that higher Pearson correlation means better recovery performance.
Furthermore, the CPU time is used as the performance index, which reflects the speed of algorithms. The computer used in the experiments is a laptop equipped with an Intel Core-i7 6500U (Intel Corporation Co., Ltd., Santa Clara, CA, USA) at 2.50 GHz, 8-GB RAM, Windows 10 64-bit, and MATLAB R2016a (MathWorks, Inc., Natick, MA, USA).

4.4. The Compared Algorithms

In order to demonstrate the performance improvement of the proposed algorithm for BSBL framework, we compared it with some classical BSBL-based algorithms including BSBL-BO [25], BSBL-L1 [25] and BSBL-FM [27], in addition to another typical CS algorithm Accelerated Iterative Hard Thresholding (AIHT) [19], which is used for comparison in order to show the advantage of BSBL framework only for sparse signal recovery. All of these algorithms would be used to recover the ECG signals directly without a dictionary matrix.

4.5. Results on MIT-BIH Arrhythmia and Long-Term ECG Databases

In general, real-time wearable ECG telemonitoring has certain requirements of the speed and accuracy of signal recovery, thus the speed and accuracy of the signal recovery algorithm must achieve an acceptable range at the same time. Figure 4 shows the recovery results averaged over all the 48 recordings in the MIT-BIH Arrhythmia Database and seven recordings in the MIT-BIH Long-Term ECG Database. From Figure 4a–d, it can be seen that the proposed BSBL-ADMM algorithm has a very close recovery accuracy relative to BSBL-BO at different CRs, and, when the CR is 80%, the recovery accuracy of the proposed BSBL-ADMM algorithm is even slightly better than BSBL-BO [25]. The other three algorithms BSBL-L1 [25] and BSBL-FM [27] and AIHT [19] have much worse recovery accuracy compared with BSBL-ADMM and BSBL-BO.
In addition, Figure 4e–f shows the average CPU times that the algorithms cost for recovering a single frame of ECG signal, we can see that the speed of the proposed BSBL-ADMM algorithm is much faster than BSBL-BO and BSBL-FM. Although the speed of BSBL-L1 and AIHT is faster than BSBL-ADMM, the recovery accuracy of these two algorithms is very bad. In fact, higher CR usually brings about less CPU time, but the recovery accuracy will be worse, so it is necessary to choose a compromise between recovery accuracy and CPU time according to practical requirements. In our experiments, we choose CR = 60% as the result of a compromise that is acceptable considering both the accuracy and speed, and more detailed results at the CR of 60% will be presented.
To better compare the performance of the proposed BSBL-ADMM algorithm with the other four algorithms, Figure 5 plots the histogram of the difference between the PRD of BSBL-ADMM and the PRD of the other four algorithms at each frame over all the recordings in the MIT-BIH Arrhythmia Database and MIT-BIH Long-Term ECG Database at the CR (compression ratios) of 60%, i.e., the histogram of p ( i ) q ( i ) , where p ( i ) indicates the PRD of BSBL-ADMM at the ith frame, and q ( i ) indicates that of BSBL-BO or other compared algorithms at the ith frame. The histogram is a statistical comparison of the recovery accuracy of two algorithms. According to the distribution of the difference between PRDs of two algorithms in the histogram, if the part that is less than 0 is in the majority, it means that the former has the better recovery accuracy, and vice versa. In Figure 5, it can be seen that the PRD of BSBL-ADMM is significantly smaller than BSBL-L1, BSBL-FM and AIHT. In addition, the PRD of BSBL-ADMM is larger than BSBL-BO, but the difference is relatively small and the recovery accuracy of the proposed BSBL-ADMM and BSBL-BO is very closed.
Figure 6 demonstrates the robustness of the proposed BSBL-ADMM algorithm. It shows the specific results of eight randomly selected recordings in the MIT-BIH Arrhythmia Database and all seven of the recordings in the MIT-BIH Long-Term ECG Database at the CR of 60%. For each recording, the proposed BSBL-ADMM algorithm has similar recovery accuracy to the average of all recordings. In terms of CPU time, the proposed BSBL-ADMM algorithm has a very stable speed for signal recovery. Therefore, these results show the robustness of the proposed BSBL-ADMM algorithm having a satisfactory performance on different recordings. The robustness ensures that the proposed BSBL-ADMM algorithm is suitable for different ECG datasets in wearable ECG telemonitoring.
Moreover, it can be clearly seen that there is a significant difference in the recovery accuracy of the same algorithm between these two databases. This is mainly because the sampling rates of the two databases are different: the sampling rate of the MIT-BIH Arrhythmia Database is 360 Hz, and the sampling rate of the MIT-BIH Long-Term ECG Database is 128 Hz. Usually a higher sampling rate will bring better recovery accuracy, but the higher sampling rate means more energy consumption of the telemonitoring system. In general, the frequency range of an ECG signal varies between 0.05–100 Hz, thus 250 Hz is a commonly used sampling rate for ECG monitoring in body area networks [35]. We sampled the ECG signal at the sampling rate of 250 Hz in our wearable ECG telemonitoring system for the experiment in the next subsection.

4.6. Results on Practical Wearable ECG Datasets

In order to verify applicability of the proposed BSBL-ADMM algorithm in practical application, we built a wearable ECG telemonitoring system described as the above Figure 2 and Figure 3. The practical wearable ECG datasets were collected using our wearable ECG collecting system. The datasets contain eight one-channel recordings collected from eight different people, and every recording was sampled at the sampling frequency 250 Hz and lasted about 30 min. In our experiments, for convenience to calculate the error between the recovered ECG signals and original ECG signals, we transmitted the original ECG signals to the computer and then both the compression and recovery of ECG signals were done on the computer.
Similarly, the experimental results on practical wearable ECG datasets also demonstrate the outstanding performance of the proposed BSBL-ADMM algorithm. The average results of practical wearable ECG datasets are presented in Figure 7. From Figure 7a,b, we can see a very similar result to MIT-BIH Arrhythmia Database that the proposed BSBL-ADMM algorithm has a very close recovery accuracy relative to BSBL-BO, and the other three algorithms have much worse accuracy. Meanwhile, Figure 7c shows the fast speed of the proposed BSBL-ADMM algorithm. Although BSBL-L1 and AIHT have the fastest speeds among these algorithms, both of their recovery accuracies are very bad.
To further show the outstanding performance and robustness of the proposed algorithm, we presented the more detailed results of these algorithms on our practical wearable ECG datasets. In order to visually demonstrate the recovery quality of the algorithms, we selected a frame of ECG signal of 1000 sampling points in recording 3 and recording 7, respectively, and presented the original and recovered ECG signals at the CR of 60% in Figure 8. It can be seen that the recovered ECG signals by BSBL-ADMM and BSBL-BO are very similar and have no obvious distortion compared with the original signal, but the recovered ECG signals by BSBL-L1, BSBL-FM and AIHT are obviously distorted.
Moreover, Table 1 shows the details of average PRD and Pearson correlation of every recording in our practical wearable ECG datasets at the CR of 60%. From Table 1, we can see that the recovery accuracy of the proposed BSBL-ADMM algorithm for all eight of the recordings in the practical wearable ECG datasets has little difference, which reflects the robustness. The average PRD and Pearson correlation of the eight recordings are 6.92 and 0.9974, respectively, and the corresponding standard deviations are 0.43 and 0.0003, respectively. Thus, the proposed algorithm is robust to different ECG recordings. Furthermore, Table 2 shows the details of average CPU time of every recording in our practical wearable ECG datasets at the CR of 60%. In addition, it can be seen that the CPU time cost by the proposed algorithm on all eight of the recordings in our practical wearable ECG datasets is relatively stable, the average CPU time of the eight recordings is 0.0629 s, and the standard deviation is 0.0006 s. Therefore, these results demonstrate that the proposed BSBL-ADMM algorithm not only has a fast speed and satisfactory accuracy but also is robust for non-sparse ECG signals in practical wearable ECG telemonitoring systems.

5. Discussion

All of these experimental results on the three different datasets demonstrate that the proposed BSBL-ADMM algorithm can recover the non-sparse ECG signals with both satisfactory recovery accuracy and fast speed compared with the other four algorithms. At the same time, the performance of the proposed BSBL-ADMM algorithm is also robust for different datasets. These advantages make the proposed BSBL-ADMM algorithm suitable for real-time ECG telemonitoring.
However, these conclusions are only considered from the perspective of the algorithm itself, and the experiments are carried out on the simulation platform. When the algorithm is actually applied to a practical wearable ECG telemonitoring system, more aspects need to be considered, such as efficient implementation of algorithmic code, the possibility to operate in real time, power consumption and the wireless link. Despite this, the proposed BSBL-ADMM algorithm is still very promising for use in wearable ECG telemonitoring systems, and we will try to actually apply the proposed BSBL-ADMM algorithm to a practical wearable ECG telemonitoring system in our future work.

6. Conclusions

In this paper, we proposed a fast and robust non-sparse signal recovery algorithm termed BSBL-ADMM for wearable ECG telemonitoring. This algorithm can recover the non-sparse ECG signals accurately acquired in wearable ECG telemonitoring in real time by using ADMM to accelerate the speed of BSBL framework. Experimental results on the MIT-BIH Databases and our practical wearable ECG datasets showed that the proposed BSBL-ADMM algorithm can recover non-sparse ECG signals fast with a satisfactory accuracy, and it also has robust performance on different ECG datasets. All these advantages indicate that the proposed BSBL-ADMM algorithm will be suitable for wearable ECG telemonitoring.

Author Contributions

Conceptualization, Y.C.; Methodology, Y.C. and Y.Y.; Software, Y.C.; Validation, W.H. and X.D.; Formal Analysis, Y.L.; Investigation, X.D.; Resources, M.H.; Data Curation, Y.L.; Writing—Original Draft Preparation, Y.C.; Writing—Review & Editing, Y.Y.; Visualization, W.H.; Supervision, M.H.; Project Administration, Y.Y.; Funding Acquisition, M.H.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 61501096 and Grant No. 61472067), by the International Science and Technology Cooperation and Exchange Program of Sichuan Province, China (Grant No. 2016HH0020), by the Chengdu Research Institute of UESTC (Grant No. RWS-CYHKF-02-20150005), by the Sichuan Science and Technology Support Project (2015GZ0199 and 2016FZ0105), and by the Fundamental Research Funds for the Central Universities (Grant No. ZYGX2016J089).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martin, T.; Jovanov, E.; Raskovic, D. Issues in wearable computing for medical monitoring applications: A case study of a wearable ECG monitoring device. In Proceedings of the Fourth International Symposium on Wearable Computers, Atlanta, GA, USA, 16–17 October 2000; pp. 43–49. [Google Scholar]
  2. Nemati, E.; Deen, M.J.; Mondal, T. A wireless wearable ECG sensor for long–term applications. IEEE Commun. Mag. 2012, 50. [Google Scholar] [CrossRef]
  3. Dey, N.; Ashour, A.S.; Shi, F.; Fong, S.J.; Sherratt, R.S. Developing residential wireless sensor networks for ECG healthcare monitoring. IEEE Trans. Consum. Electron. 2017, 63, 442–449. [Google Scholar] [CrossRef]
  4. Craven, D.; McGinley, B.; Kilmartin, L.; Glavin, M.; Jones, E. Energy–efficient Compressed Sensing for ambulatory ECG monitoring. Comput. Biol. Med. 2016, 71, 1–3. [Google Scholar] [CrossRef] [PubMed]
  5. Sun, F.T.; Kuo, C.; Griss, M. PEAR: Power efficiency through activity recognition (for ECG–based sensing). In Proceedings of the 5th International Conference on Pervasive Computing Technologies for Healthcare (PervasiveHealth), Dublin, Ireland, 23–26 May 2011; pp. 115–122. [Google Scholar]
  6. Cao, H.; Leung, V.; Chow, C.; Chan, H. Enabling technologies for wireless body area networks: A survey and outlook. IEEE Commun. Mag. 2009, 47. [Google Scholar] [CrossRef]
  7. Ma, T.; Shrestha, P.L.; Hempel, M.; Peng, D.; Sharif, H.; Chen, H.H. Assurance of energy efficiency and data security for ECG transmission in BASNs. IEEE Trans. Biomed. Eng. 2012, 59, 1041–1048. [Google Scholar] [PubMed]
  8. Dixon, A.M.; Allstot, E.G.; Gangopadhyay, D.; Allstot, D.J. Compressed sensing system considerations for ECG and EMG wireless biosensors. IEEE Trans. Biomed. Circuits Syst. 2012, 6, 156–166. [Google Scholar] [CrossRef] [PubMed]
  9. Candès, E.J.; Romberg, J.; Tao, T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 2006, 52, 489–509. [Google Scholar] [CrossRef]
  10. Candes, E.J.; Tao, T. Near–optimal signal recovery from random projections: Universal encoding strategies? IEEE Trans. Biomed. Circuits Syst. 2006, 52, 5406–5425. [Google Scholar] [CrossRef]
  11. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  12. Craven, D.; McGinley, B.; Kilmartin, L.; Glavin, M.; Jones, E. Compressed sensing for bioelectric signals: A review. IEEE J. Biomed. Health Inf. 2015, 19, 529–540. [Google Scholar] [CrossRef] [PubMed]
  13. Mamaghanian, H.; Khaled, N.; Atienza, D.; Vandergheynst, P. Compressed sensing for real–time energy–efficient ECG compression on wireless body sensor nodes. IEEE Trans. Biomed. Eng. 2011, 58, 2456–2466. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, Y.; Doleschel, S.; Wunderlich, R.; Heinen, S. High energy efficient analog compressed sensing encoder for wireless ECG system. Microelectron. J. 2016, 56, 10–16. [Google Scholar] [CrossRef]
  15. Polania, L.F.; Carrillo, R.E.; Blanco–Velasco, M.; Barner, K.E. Compressed sensing based method for ECG compression. In Proceedings of the 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Prague, Czech Republic, 22–27 May 2011; pp. 761–764. [Google Scholar]
  16. Polania, L.F.; Carrillo, R.E.; Blanco–Velasco, M.; Barner, K.E. Exploiting prior knowledge in compressed sensing wireless ECG systems. IEEE J. Biomed. Health Inf. 2015, 19, 508–519. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, Z.; Rao, B.D.; Jung, T.P. Compressed sensing for energy–efficient wireless telemonitoring: Challenges and opportunities. In Proceedings of the 2013 Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 3–6 November 2013; pp. 80–85. [Google Scholar]
  18. Needell, D.; Tropp, J.A. CoSaMP: iterative signal recovery from incomplete and inaccurate samples. Commun. ACM 2010, 53, 93–100. [Google Scholar] [CrossRef]
  19. Blumensath, T. Accelerated iterative hard thresholding. Signal Process. 2012, 92, 752–756. [Google Scholar] [CrossRef] [Green Version]
  20. Mohimani, H.; Babaie–Zadeh, M.; Jutten, C. A Fast Approach for Overcomplete Sparse Decomposition Based on Smoothed l0 Norm. IEEE Trans. Signal Process. 2009, 57, 289–301. [Google Scholar] [CrossRef]
  21. Eldar, Y.C.; Kuppinger, P.; Bolcskei, H. Block–sparse signals: Uncertainty relations and efficient recovery. IEEE Trans. Signal Process. 2010, 58, 3042–3054. [Google Scholar] [CrossRef]
  22. Huang, J.; Zhang, T.; Metaxas, D. Learning with structured sparsity. J. Mach. Learn. Res. 2011, 12, 3371–3412. [Google Scholar]
  23. Peleg, T.; Eldar, Y.C.; Elad, M. Exploiting statistical dependencies in sparse representations for signal recovery. IEEE Trans. Signal Process. 2012, 60, 2286–2303. [Google Scholar] [CrossRef]
  24. Zhang, Z.; Jung, T.P.; Makeig, S.; Rao, B.D. Compressed sensing for energy–efficient wireless telemonitoring of noninvasive fetal ECG via block sparse Bayesian learning. IEEE Trans. Biomed. Eng. 2013, 60, 300–309. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, Z.; Rao, B.D. Extension of SBL algorithms for the recovery of block sparse signals with intra–block correlation. IEEE Trans. Signal Process. 2013, 61, 2009–2015. [Google Scholar] [CrossRef]
  26. Zhang, Z.; Jung, T.P.; Makeig, S.; Rao, B.D. Compressed sensing of EEG for wireless telemonitoring with low energy consumption and inexpensive hardware. IEEE Trans. Biomed. Eng. 2013, 60, 221–224. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, B.; Zhang, Z.; Xu, G.; Fan, H.; Fu, Q. Energy efficient telemonitoring of physiological signals via compressed sensing: A fast algorithm and power consumption evaluation. Biomed. Signal Process. Control 2014, 11, 80–88. [Google Scholar] [CrossRef] [Green Version]
  28. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed optimization and statistical learning via the alternating direction method of multipliers. Found. Trends Mach. Learn. 2011, 3, 1–122. [Google Scholar] [CrossRef]
  29. Moody, G.B.; Mark, R.G. The impact of the MIT–BIH arrhythmia database. IEEE Eng. Med. Biol. Mag. 2001, 20, 45–50. [Google Scholar] [CrossRef] [PubMed]
  30. Goldberger, A.L.; Amaral, L.A.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.K.; Stanley, H.E. Physiobank, physiotoolkit, and physionet. Circulation 2000, 13, e215–e220. [Google Scholar] [CrossRef]
  31. Yuan, M.; Lin, Y. Model selection and estimation in regression with grouped variables. J. R. Stat. Soc. Ser. B 2006, 68, 49–67. [Google Scholar] [CrossRef] [Green Version]
  32. Kiranyaz, S.; Ince, T.; Gabbouj, M. Real–time patient–specific ECG classification by 1–D convolutional neural networks. IEEE Trans. Biomed. Eng. 2016, 63, 664–675. [Google Scholar] [CrossRef] [PubMed]
  33. Bhanot, K.; Peddoju, S.K.; Bhardwaj, T. A model to find optimal percentage of training and testing data for efficient ECG analysis using neural network. Int. J. Syst. Assur. Eng. Manag. 2018, 9, 12–17. [Google Scholar] [CrossRef]
  34. Agung, M.A.; Basari. 3–lead acquisition using single channel ECG device developed on AD8232 analog front end for wireless ECG application. InAIP Conf. Proc. 2017. [Google Scholar] [CrossRef]
  35. Paradiso, R.; Loriga, G.; Taccini, N. A wearable health care system based on knitted integrated sensors. IEEE Trans. Inf. Technol. Biomed. 2005, 9, 337–344. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Typical compressed sensing (CS)-based ECG wireless telemonitoring system.
Figure 1. Typical compressed sensing (CS)-based ECG wireless telemonitoring system.
Sensors 18 02021 g001
Figure 2. Block diagram of a wearable ECG telemonitoring system.
Figure 2. Block diagram of a wearable ECG telemonitoring system.
Sensors 18 02021 g002
Figure 3. Devices of a wearable ECG collecting system.
Figure 3. Devices of a wearable ECG collecting system.
Sensors 18 02021 g003
Figure 4. Recovery results of all 48 of the recordings in MIT-BIH Arrhythmia Database and seven recordings in MIT-BIH Long-Term ECG Database at different CRs (compression ratios). (ad) show that the PRD (percentage root-mean squared distortion) and Pearson correlation of BSBL-ADMM are very close to BSBL-BO, the other three have much worse recovery accuracy; (e,f) show that the BSBL-ADMM is much faster compared with BSBL-BO. These results on MIT-BIH Databases demonstrate that the proposed BSBL-ADMM algorithm has both satisfactory speed and accuracy compared with the other algorithms.
Figure 4. Recovery results of all 48 of the recordings in MIT-BIH Arrhythmia Database and seven recordings in MIT-BIH Long-Term ECG Database at different CRs (compression ratios). (ad) show that the PRD (percentage root-mean squared distortion) and Pearson correlation of BSBL-ADMM are very close to BSBL-BO, the other three have much worse recovery accuracy; (e,f) show that the BSBL-ADMM is much faster compared with BSBL-BO. These results on MIT-BIH Databases demonstrate that the proposed BSBL-ADMM algorithm has both satisfactory speed and accuracy compared with the other algorithms.
Sensors 18 02021 g004
Figure 5. Histogram of the difference between the PRD of BSBL-ADMM and the PRD of the other four algorithms at each frame over all the recordings in MIT-BIH Arrhythmia Database and MIT-BIH Long-Term ECG Database at the CR of 60%. It can be seen that the PRD of BSBL-ADMM is significantly smaller than BSBL-L1, BSBL-FM and AIHT. In addition, the PRD of BSBL-ADMM is larger than BSBL-BO, but the difference is relatively small and the recovery accuracy of the proposed BSBL-ADMM and BSBL-BO is very close.
Figure 5. Histogram of the difference between the PRD of BSBL-ADMM and the PRD of the other four algorithms at each frame over all the recordings in MIT-BIH Arrhythmia Database and MIT-BIH Long-Term ECG Database at the CR of 60%. It can be seen that the PRD of BSBL-ADMM is significantly smaller than BSBL-L1, BSBL-FM and AIHT. In addition, the PRD of BSBL-ADMM is larger than BSBL-BO, but the difference is relatively small and the recovery accuracy of the proposed BSBL-ADMM and BSBL-BO is very close.
Sensors 18 02021 g005
Figure 6. Recovery results of eight randomly selected recordings in MIT-BIH Arrhythmia Database and all seven of the recordings in MIT-BIH Long-Term ECG Database. (a,b) show that the BSBL-ADMM has relatively stable recovery accuracy for different recordings; (c,d) show that the the BSBL-ADMM has relatively stable CPU times for different recordings. All these results on MIT-BIH Databases demonstrate the robustness of the proposed BSBL-ADMM algorithm.
Figure 6. Recovery results of eight randomly selected recordings in MIT-BIH Arrhythmia Database and all seven of the recordings in MIT-BIH Long-Term ECG Database. (a,b) show that the BSBL-ADMM has relatively stable recovery accuracy for different recordings; (c,d) show that the the BSBL-ADMM has relatively stable CPU times for different recordings. All these results on MIT-BIH Databases demonstrate the robustness of the proposed BSBL-ADMM algorithm.
Sensors 18 02021 g006
Figure 7. Recovery results of all the recordings in practical wearable ECG datasets at different CRs (compression ratios). (a,b) show that the PRD (percentage root-mean squared distortion) and Pearson correlation of BSBL-ADMM is very close to BSBL-BO, the other two have much worse recovery accuracy; (c) shows that the BSBL-ADMM has much faster speed compared with BSBL-BO. These results on practical wearable ECG datasets demonstrate that the proposed BSBL-ADMM algorithm has both satisfactory speed and accuracy compared with the other algorithms.
Figure 7. Recovery results of all the recordings in practical wearable ECG datasets at different CRs (compression ratios). (a,b) show that the PRD (percentage root-mean squared distortion) and Pearson correlation of BSBL-ADMM is very close to BSBL-BO, the other two have much worse recovery accuracy; (c) shows that the BSBL-ADMM has much faster speed compared with BSBL-BO. These results on practical wearable ECG datasets demonstrate that the proposed BSBL-ADMM algorithm has both satisfactory speed and accuracy compared with the other algorithms.
Sensors 18 02021 g007
Figure 8. Recovery results of two ECG signals selected from recording 3 and recording 7, respectively. It can be seen intuitively that the recovered ECG signals by BSBL-ADMM and BSBL-BO are very similar and have no obvious distortion compared with the original signal, but the recovered ECG signals by BSBL-L1, BSBL-FM and AIHT are obviously distorted.
Figure 8. Recovery results of two ECG signals selected from recording 3 and recording 7, respectively. It can be seen intuitively that the recovered ECG signals by BSBL-ADMM and BSBL-BO are very similar and have no obvious distortion compared with the original signal, but the recovered ECG signals by BSBL-L1, BSBL-FM and AIHT are obviously distorted.
Sensors 18 02021 g008
Table 1. The details of recovery accuracy of practical wearable ECG datasets (CR = 60%). These results show that the proposed BSBL-ADMM algorithm is robust for the recovery accuracy for different ECG recordings.
Table 1. The details of recovery accuracy of practical wearable ECG datasets (CR = 60%). These results show that the proposed BSBL-ADMM algorithm is robust for the recovery accuracy for different ECG recordings.
AlgorithmsRec1Rec2Rec3Rec4Rec5Rec6Rec7Rec8Mean (±Std)
PRDBSBL-ADMM7.346.636.697.636.356.686.827.226.92 (±0.43)
BSBL-BO [25]7.566.356.188.175.696.426.916.606.74 (±0.80)
BSBL-L1 [25]54.1059.7678.0071.6066.4960.3476.3775.4467.76 (±8.94)
BSBL-FM [27]35.7734.1835.1039.3335.7134.8835.9138.6336.19 (±1.82)
AIHT [19]101.98103.74105.54106.57104.40103.06103.04104.55104.11 (±1.48)
Pearson CorrelationBSBL-ADMM0.99720.99770.99770.99700.99800.99770.99760.99740.9976 (±0.0003)
BSBL-BO [25]0.99700.99790.99800.99660.99830.99790.99760.99780.9976 (±0.0006)
BSBL-L1 [25]0.84360.80800.66010.72360.76150.80400.68040.69370.7469 (±0.0675)
BSBL-FM [27]0.93290.93740.93210.91620.93270.93590.92790.91920.9293 (±0.0077)
AIHT [19]0.40900.39610.36890.35480.36690.39590.37440.36570.3790 (±0.0190)
Table 2. The details of recovery CPU Time of practical wearable ECG datasets (CR = 60%). These results show that the proposed BSBL-ADMM algorithm is robust for the recovery speed for different ECG recordings.
Table 2. The details of recovery CPU Time of practical wearable ECG datasets (CR = 60%). These results show that the proposed BSBL-ADMM algorithm is robust for the recovery speed for different ECG recordings.
Rec1Rec2Rec3Rec4Rec5Rec6Rec7Rec8Mean (±Std)
BSBL-ADMM0.06390.06350.06310.06250.06250.06240.06240.06280.0629 (±0.0006)
BSBL-BO [25]0.23590.23610.23530.23520.23550.23500.23500.23580.2355 (±0.0004)
BSBL-L1 [25]0.02940.02930.02910.02900.02910.02910.02900.02910.0292 (±0.0001)
BSBL-FM [27]0.16050.20850.28680.25410.21700.20380.28260.25900.2340 (±0.0044)
AIHT [19]0.04580.03770.03720.04440.04250.04590.04070.04140.0419 (±0.0034)

Share and Cite

MDPI and ACS Style

Cheng, Y.; Ye, Y.; Hou, M.; He, W.; Li, Y.; Deng, X. A Fast and Robust Non-Sparse Signal Recovery Algorithm for Wearable ECG Telemonitoring Using ADMM-Based Block Sparse Bayesian Learning. Sensors 2018, 18, 2021. https://doi.org/10.3390/s18072021

AMA Style

Cheng Y, Ye Y, Hou M, He W, Li Y, Deng X. A Fast and Robust Non-Sparse Signal Recovery Algorithm for Wearable ECG Telemonitoring Using ADMM-Based Block Sparse Bayesian Learning. Sensors. 2018; 18(7):2021. https://doi.org/10.3390/s18072021

Chicago/Turabian Style

Cheng, Yunfei, Yalan Ye, Mengshu Hou, Wenwen He, Yunxia Li, and Xuesong Deng. 2018. "A Fast and Robust Non-Sparse Signal Recovery Algorithm for Wearable ECG Telemonitoring Using ADMM-Based Block Sparse Bayesian Learning" Sensors 18, no. 7: 2021. https://doi.org/10.3390/s18072021

APA Style

Cheng, Y., Ye, Y., Hou, M., He, W., Li, Y., & Deng, X. (2018). A Fast and Robust Non-Sparse Signal Recovery Algorithm for Wearable ECG Telemonitoring Using ADMM-Based Block Sparse Bayesian Learning. Sensors, 18(7), 2021. https://doi.org/10.3390/s18072021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop