Journal of Healthcare Engineering
Volume 2018 (2018), Article ID 5396030, 10 pages
https://doi.org/10.1155/2018/5396030

Research Article

Analysis for the Influence of ABR Sensitivity on PTT-Based Cuff-Less Blood Pressure Estimation before and after Exercise

Yang Xu, Peng Ping, Dong Wang, and Weigong Zhang

School of Instrument Science and Engineering at Southeast University, Nanjing, China

Correspondence should be addressed to Dong Wang; kingeast16@seu.edu.cn and Weigong Zhang; weig_zhang@163.com

Received 18 June 2018; Revised 23 August 2018; Accepted 6 September 2018; Published 8 October 2018

Academic Editor: Vincenzo Positano

Copyright © 2018 Yang Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An accurate and continuous measurement of blood pressure (BP) is of great importance for the prognosis of some cardiovascular diseases in out-of-hospital settings. Pulse transit time (PTT) is a well-known cardiovascular parameter which is highly correlated with BP and has been widely applied in the estimation of continuous BP. However, due to the complexity of cardiovascular system, the accuracy of PTT-based BP estimation is still unsatisfactory. Recent studies indicate that, for the subjects before and after exercise, PTT can track the high-frequency BP oscillation (HF-BP) well, but is inadequate to follow the low-frequency BP variance (LF-BP). Unfortunately, the cause for this failure of PTT in LF-BP estimation is still unclear. Based on these previous researches, we investigated the cause behind this failure of PTT in LF-BP estimation. The heart rate- (HR-) related arterial baroreflex (ABR) model was introduced to analyze the failure of PTT in LF-BP estimation. Data from 42 healthy volunteers before and after exercise were collected to evaluate the correlation between the ABR sensitivity and the estimation error of PTT-based BP in LF and HF components. In the correlation plot, an obvious difference was observed between the LF and HF groups. The correlation coefficient r for the ABR sensitivity with the estimation error of systolic BP (SBP) and diastolic BP (DBP) in LF was 0.817 ± 0.038 and 0.757 ± 0.069, respectively. However, those correlation coefficient r for the ABR sensitivity with the estimation error of SBP and DBP in HF was only 0.403 ± 0.145 and 0.274 ± 0.154, respectively. These results indicated that there is an ABR-related complex LF autonomic regulation mechanism on BP, PTT, and HR, which influences the effect of PTT in LF-BP estimation.

1. Introduction

Blood pressure (BP) is a vital sign which is defined as the pressure of circulating blood on the walls of blood vessels. In morphology, BP is usually expressed in terms of the systolic BP (maximum during one heartbeat) over diastolic BP (minimum in between two heartbeats). Continuous blood pressure monitoring is of great clinical significance which could provide a long-term character of the cardiac system. It is important to the diagnosis of cardiac diseases. However, the traditional BP monitoring involves the inflating and deflating of a cuff which causes unavoidable intermittent monitoring of the blood pressure. Because of this shortcoming of cuff-based BP monitoring and the profound need for continuous blood pressure monitoring, the scientific community has paid a lot of attention to the cuff-less continuous BP estimation methods in recent decades [1].

Among these cuff-less BP measurement methods, the PTT-based method has been extensively investigated in clinical settings. A mountain of research studies has demonstrated that there is a high correlation between BP and PTT [2–5]. PTT is a cardiovascular parameter which is defined as the time delay for the pressure wave to travel from a proximal point to a distal point in the arterial within the same cardiac cycle. It could be easily calculated from the feature points of ECG and PPG. The principal fundamental behind this method is that the blood flow in the arteries is physically modeled as the propagation of pressure waves inside elastic tubes. In this way, the physical characteristics of blood flow wave in arteries could be described by the Moens–Korteweg (MK) equation as follows [6]:

where , , , , and are the length of the vessel, the elastic modulus of the tube wall, the vessel wall thickness, the vessel radius, and the constant blood density, respectively. PWV is the velocity of the blood flow wave propagating in the entire arterial tree, which is inversely related with PTT. is exponentially correlated with BP through

where is the elastic modulus at zero pressure, is a coefficient depending on a particular vessel, and is BP. The MK equation and its variants could also be simplified into different mathematical models, in which the relevant hemodynamic parameters (blood density, vessel radius, etc.,) were converted into a number of individual-specific model parameters (see Section 2.1 for details). Therefore, with an initial calibration of these individual-specific model parameters, the PTT-BP relationship is obtained. Then, measuring PTT could offer a continuous and cuff-less BP monitoring [2].

Over the past 20 years, various calibration models of the BP-PTT relationship have been proposed to achieve the continuous BP measurement with PTT [7–21]. Chen et al. [7] established a mathematical model involving the patients’ ages and genders. It indicated that this model could capture the PTT-DBP relationship for subjects on a wider range. Gesche et al. [11] developed a one-point calibration model for the PTT-BP relationship, which needs only one measurement of BP using a cuff-based reference. The results indicated that the SBP calculated from PTT correlates significantly with the cuff-based SBP (). Esmaili et al. [18] proposed a PTT-BP nonlinear model for the accurate estimation of both SBP and DBP. It attained a high accuracy with evaluated error and variance of 0.12 ± 6.15 mmHg for SBP and 1.31 ± 5.36 mmHg for DBP.

Although the PTT-based method has been considered as the most promising cuff-less continuous BP monitoring technique, there are still several problems that need to be solved before its widespread application [8]. The major and most important challenge is that the accuracy of PTT-based BP estimation is still not satisfactory enough. The possible influence factors include: arterial compliance, cardiac output, peripheral resistance, and blood volume [10]. However, due to the complexity of the cardiovascular system, the influences of these physiological factors on the PTT-based BP are hard to be evaluated. It is impracticable to employ these factors directly to improve the accuracy of PTT-based BP estimation. The possible solution is, based on the experiment phenomenon, to introduce the applicable regulation models or parameters into the PTT-BP calibration model to weaken the influences from other physiological factors.

From the differences of influence factors, there are two special frequency bands of the BP changes [22–24]: (1) the high-frequency (HF, 0.2–0.35 Hz) BP oscillation influenced by the physical factors (e.g., respiration) and (2) the low-frequency (LF, 0.1–0.15 Hz) BP variance under the control of the autonomic nervous system (ANS). Recently, various researches [25–27] have investigated the time-frequency correlation between the variability of PTT, HR, and BP before and after exercise. The results indicate that, relatively speaking, PTT can track the HF BP oscillation well, but inadequate to follow the LF BP variance (see Section 2.3 for details). From the experimental results, this poor performance of PTT in LF-BP estimation is considered as a reason for the nonaccuracy of PTT-based BP estimation [10]. However, the cause for this failure of PTT in LF-BP estimation is still unknown.

Based on these previous researches, we try to introduce the ABR regulation model [28] (arterial baroreflex, a major BP-related autonomic nervous regulation mechanism, see section 2.2 for details) to analyze this failure of PTT in LF-BP estimation in this research. From the view of the ABR model and the interaction between the cardiovascular parameters, it was found that the ABR sensitivity may influence the estimation accuracy of PTT-based LF-BP. Then, synchronous physiological data (BP, PPG, and ECG) from 42 subjects before and after exercises were collected to quantitatively analyze the influence of ABR sensitivity on the PTT-based HF- and LF-BP with the help of VMD (variation mode decomposition, a time-frequency analysis technique). The correlation between the ABR sensitivity and the estimation error of PTT-based HF- and LF-BP was evaluated. The aim was to analyze the cause behind failure of PTT in LF-BP estimation in the experimental condition. It is expected that the results of this work may help to construct a better PTT-BP calibration model in practice.

The rest of this paper is organized as follows. Section 2 introduces the backgrounds of this paper including the current PTT-based calibration model, the ABR model, the time-frequency correlation between BP, HR, and PTT before and after exercise, and the qualitative analysis of the possible cause for the failure of PTT in LF-BP estimation. Section 3 explains the quantitative analysis methodology including the data collection and preprocessing, the calculation method of ABR sensitivity, and the correlation analyses between the ABR sensitivity and the PTT-based BP estimation error. Section 4 presents and discusses the results of the experiments. The conclusion and further works are given in Section 5.

2. Backgrounds

2.1. The Calibration Models of PTT-BP Relationship

In the PTT-based BP estimation method, the arterial vessels are physically modeled as an elastic tube. The pressure wave propagating on the vessels follows a function of the pulse propagating position and time ( and , respectively) as [4]

where is a constant that represents the arterial inertance per unit length. The vessel compliance is defined as the rate of tube cross section changes in terms of blood pressure as follows:

where , , and are subject-specific parameters [5]. Accordingly, the PWV is equal to , and PTT is the time interval for the pressure wave traversing a tube of length , which is expressed as:

Then, combining equations (3)–(5) with the aforementioned equations (1) and (2), the two most popular calibration models in the literature could be obtained. By substituting (1) into (2), a popular calibration model is given as follows [7]:

The other popular physical calibration model is derived by substituting (4) into (5), given as follows [8]:

where and are unknown subject-specific parameters which are obtained by fitting reference BP with PTT in the regression model.

Experimental studies have shown that 1/PTT, rather than PTT, is linearly related to BP [11, 12]. In other words, model (7) performs better in PTT-BP fitting than model (6). However, no matter model (6) or model (7) are the MK equation-based models, which could only indicate the rhythmic HF-BP oscillation that was caused by the physical activity [6]. The BP changes are not only HF-oscillating caused by the exogenous driving but also LF-fluctuating due to the autonomic nervous regulation [23].

2.2. Arterial Baroreflex Regulation Model

Arterial baroreflex (ABR) is a major autonomic nervous regulation mechanism which is responsible to stabilize BP [28]. The baroreceptors, locating along all major arteries of the human body, sense the BP changes and deliver a signal to the autonomous nervous system (ANS). Then, the ANS stabilizes BP through a feedback regulation mechanism by HR like this: decrease (increase) in BP leads to the consequent reduction (enhancement) in the HR. Meanwhile, the increase (decrease) in HR causes the increase (decrease) in BP as a direct feedforward effect. Therefore, the ABR mechanism could be described by a simplified double-loop feedback diagram (Figure 1).

Figure 1: The simplified diagram of ABR.

The overall ABR regulation nonlinearly drives the heart rate on the basis of the arterial pressure, which follows a logistic model [28]:

where, and are the lower and upper levels of heart rate, respectively, denotes the arterial pressure at the midpoint of the heart rate range, and determines the slope of the linear region (or sensitive region) in the overall ABR regulation curve (Figure 2).

Figure 2: The schematic diagram of the overall ABR regulation model.

2.3. The Time-Frequency Correlation Between BP, HR, and PTT before and after Exercise

In the past studies, numerous studies [25–27] have investigated the time-frequency correlation between the variability of PTT, HR, and BP before and after exercise. The variability of these physiological data is defined as the variation between two consecutive data points. Drinnan et al. [25] analyzed the cross-correlation function between PTT and HR from 15 normal healthy subjects during paced respiration. The results suggest that the high HRV would significantly influence the PTTV. This relationship between HRV and PTTV is negatively correlated, i.e.,

Recently, using the recursive autoregressive model, Ma and Zhang [26] and Liu et al. [27] investigated the time-frequency correlation between BP, HR, and PTT in the LF and HF components for the subjects before and after exercise. The results reveal that PTT is highly correlated with HF-BP changes, but insignificantly correlated with LF-BP changes. However, the cause behind this phenomenon is still unclear. The possible reason is the mediation mechanisms of ANS on the cardiovascular system and the inherent correlations between the cardiovascular parameters.

2.4. Qualitative Analysis for the Failure of PTT in LF-BP Estimation

Based on these research studies and the ABR model, a qualitative analysis for the ineffective of PTT in LF-BP estimation before and after exercise was proposed here.

As shown in Figure 2, under the sensitive ABR regulation (the red box region), HR is increasing rapidly (high HRV) while BP is slow varying (LF-BP changes). As aforementioned in (9), PTTV is highly coupled with HRV. In this scenario, the high HRV significantly influences the fluctuation of PTT. The PTT would change unpredictably that may no longer follow the PTT-BP relationship model in (6) or (7), and leads to the failure of PTT in the estimation of LF-BP changes. On the other hand, under the nonsensitive ABR (region outside the red box), HR is slow varying (low HRV) as BP is increasing rapidly (HF-BP). Therefore, PTT is less affected by the HRV during the nonsensitive ABR regulation and predicts the HF-BP well. In brief, the high HRV may have unexpected effects on PTT under the ABR regulation. These ABR-model-based qualitative analysis results could explain the poor performance of PTT in LF-BP estimation well. It inspired us that ABR regulation is possibly the main cause for the failure of PTT in LF-BP estimation before and after exercise. However, further quantitative analysis should be implemented to confirm that. In the following sections, we investigated the correlation between the ABR sensitivity and the estimation error of PTT-based BP in HF-band and LF-band separately.

3. Methodology

3.1. Data Collection and Protocol

In our previous research [29], synchronous ECG, PPG, and reference BP data were collected on 42 healthy adults (21 males) with a mean age of 25.6 ± 2.1 years (range 21–31 years) from 8:00 to 11:00 a.m. in a quiet environment and at a constant room temperature of 22–25°C. These volunteers were nonsmokers with no history of cardiovascular diseases (CVDs) and no caffeine ingestion 6 h prior to the examination. Reference BP including systolic blood pressure (SBP) and diastolic blood pressure (DBP) was measured by Finapres (Finapres Medical System), a noninvasive BP measurement system, with the finger cuff on the right thumb and brachial cuff on the right upper arm. ECG and PPG were acquired with lead-II ECG electrode and PPG sensor (Biopac Systems) on left middle finger, respectively. All the data collection was performed with subjects in sitting position before and after a treadmill running exercise. Treadmill running exercise is a commonly used drugless method to increase the subject’s BP variation range in trials [30]. The signals were recorded at the sampling rate of 1000 Hz. The ECG signal was filtered with a 0.5 Hz Butterworth high-pass filter and a 35 Hz Butterworth low-pass filter. The PPG signal was filtered with a 0.05 Hz Butterworth high-pass filter and a 10 Hz Butterworth low-pass filter. These filter parameters are the default setting that is recommended by the Biopac physiological data acquisition system. Table 1 summarized the whole experimental procedure for each volunteer. In total, 3 ∗ 42 ∗ 20-min ECG, PPG, and reference BP signal were collected.

Table 1: The experiment procedure.

Physiological conditions
Trial no.
Time (min)
Recording Length (min)

Acclimatization
1
10
0

Sitting rest
2
10
10

Treadmill running (8 km/h)
3
5
0

Sitting recovery
4
10
10

3.2. Parameter Extraction and Data Processing

PTT is usually calculated as the time interval between the ECG R-waveform and (1) the peak point of PPG (PTT1), (2) the maximal first derivate point of PPG (PTT2), or (3) the minimum point of PPG (PTT3) among RR intervals (RRi) in the same cardiac cycle (see Figure 3). However, due to the motion artifact, the peak point or the minimum point of PPG is easy to deform, which influences the measurement of PTT1 and PTT3 in the experiment. Therefore, we chose the relatively more robust PTT2 [6] as the measurement of PTT to ensure the test results, here. Then, the obtained PTT was applied to match the reference SBP and DBP with the calibration model (7), which is recommended by Mukkamala et al. [6]. The PTT-BP relationship model was been adaptively trained with the ordinary least squares (OLS) algorithm which can be mathematically expressed as an optimization problem as follows:

where , as given in (7) and is the referenced SBP and DBP.

Figure 3: The definition of PTT and RR interval.

In order to analyze the estimation accuracy of PTT-based BP in LF and HF range separately, VMD technology was utilized to decompose the estimated BP and reference BP (including SBP and DBP) signal. VMD is a newly developed time-frequency analysis technique proposed by Dragomiretskiy and Zosso [31]. Using VMD, the signal could be adaptively decomposed into an ensemble of N () band-limited modes (BLMs) without the need of setting frequency range artificially. Each BLM compacts around a center frequency determined by the signal itself. The bandwidth for each BLM is calculated with the L2-norm of the gradient of its Hilbert transformed signal. Thus, these BLM could capture the intrinsic characters of the signal in the time-frequency domain without artificial influence. To benefit from this unique property of BLM, VMD technology has been widely used in the biosignal processing [32–34]. Here, the reference and estimated BP signals were decomposed into LF and HF modes. This process was implemented with the help of VMD Matlab toolbox software programmed by Dragomiretskiy and Zosso [31].

Figure 4 shows one segment for the time-frequency decomposition results of the reference and the corresponding estimated BP signal. It could be observed that the spectral decomposition results of the estimated BP are obviously different from those of reference BP, especially in the LF regions. Similarly, in the time domain, the differences between the estimated BP and the reference BP are more obvious in the LF region than those in the HF region. These results are consistent with the research results of Ma and Zhang [26] and Liu et al. [27], i.e., the worse performance of PTT in LF-BP estimation.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 4: One exemplary reference BP signal decomposition results in time and frequency domain. (a) The reference DBP spectral decomposition. (b) The reference DBP. (c) The estimated DBP spectral decomposition. (d) The estimated DBP. (e) The reference SBP spectral decomposition. (f) The reference SBP. (g) The estimated SBP spectral decomposition and (h) the estimated SBP.

Then, the estimation errors of the PTT-based SBP/DBP were calculated in LF, HF, and overall as follows:

where is the estimation error, is the estimation value, the rightmost subtrahend in the equations is the corresponding reference BP value, represents the subject, and is the data point.

The ABR sensitivity is quantified by the clinical parameter baroreflex sensitivity (BRS) which is defined as the RRi variations as a reaction to the change of BP. The most widely applied BRS measuring method is the time-domain approach, in which BRS (measured in ms/mmHg) is calculated as the averaging regression coefficient between the first-order difference of RRi (measured in ms) and SBP (measured in mmHg) changes in the same direction (dRRi > 0, dSBP > 0 or dRRi < 0, dSBP < 0, i.e., 1st and 3rd quadrants) during a sliding window [35]. In this study, the BRS is dynamically calculated with a one-point-overlapping sliding window in size of 16 data points [35]. Figure 5 shows the process of dynamic BRS calculation. The details of the algorithm were described as the pseudocode in Table 2.

Figure 5: Dynamic BRS calculation. dRRi and dSBP are the first-order forward difference of RRi and reference SBP, respectively.

Table 2: The pseudocode of BRS calculation.

Initialize

The sliding-window start data-point number

Searching the sliding window:

Set inside-window start data-point number

    Repeat

        ,

    if (dRRi() < 0 and dSBP() < 0)

        Register in the first-quadrant array

    else if (dRRi() > 0 and dSBP() > 0)

        Register in the third-quadrant array

    Until

Calculate the regression coefficient between dRRi (i) and dSBP (i) in the first-quadrant and third-quadrant separately: R₁ and R₂

    BRS in the sliding-window is calculated as:

    BRS  = (R₁ + R₂)/2, .

go to Searching the sliding window;

until the last data-point is reached

3.3. Correlation Evaluation

To confirm and pinpoint the influence of ABR sensitivity on the PTT-based BP estimation, the Pearson’s correlation coefficients (represented in terms of and ) was utilized to calculate the correlation between the estimation error of PTT-based BP and BRS for each subject. Also, we applied the minimum absolute value (MAV) to test the significance of the correlation coefficients. The specific formulas of these criteria are given as follows:

where is the subject, is the number of subjects, is the number of data points, and is the minimum absolute value.

The obtained data were contrastively analyzed stage by stage. At first, the overall correlation coefficients were compared between two groups: (1) versus BRS and (2) versus BRS. It is to confirm if there is an overall correlation between BRS and the estimation error of PTT-based BP. Then, the correlations are analyzed in LF and HF sections comparatively: (1) versus BRS, (2) versus BRS, (3) versus BRS, and (4) versus BRS. It is designed to further test if there is relativity between the ABR sensitivity and the failure of PTT in LF-BP estimation. In each stage, the significance test of the difference between two groups was implemented with a T-test. It was expected to verify and confirm the influence of ABR on PTT-based BP estimation with these contrastive analyses.

4. Experiment Results and Discussion

4.1. Overall Correlation Analysis

Figure 6 shows the boxplot of the overall correlation coefficients for two groups: (1) versus BRS and (2) versus BRS. The red line represents the median value. The bottom and top of the blue box are the first and third quartiles of data distributions, respectively. And, the black line represents 1.5 times the interquartile range of upper and lower quartiles. It is observed that the mean value of overall correlation coefficients in the group versus BRS is higher than that in the group versus BRS, indicating a closer relationship between and BRS. However, it is worth noting that the correlation coefficients of each group are changing in both positive and negative regions. It reveals that the overall correlation for -BRS and -BRS is uncertain.

Figure 6: Boxplot of the overall correlation coefficients between the PTT-based SBP and DBP estimated error and BRS.

Table 3 lists the evaluation criteria for the correlation coefficients in detail. Especially, the MAV of the correlation coefficients in the group of versus BRS and versus BRS are 0.031 and 0.043, respectively. This near-zero value of MAV indicates the lack of correlation between ABR sensitivity and the overall PTT-based BP estimation error. It indicated that the ABR mechanism could not fully explain the estimation errors of PTT-based BP.

Table 3: The overall correlation for BRS and BP estimation error.

Groups
Pearson’s correlation coefficients () ()

vs BRS
0.533
0.109
0.031

vs BRS
0.411
0.145
0.043

4.2. Correlation Analysis of the ABR Sensitivity and the Estimation Error for PTT-Based LF- and HF-BP

In addition, we compared the correlation coefficients between the PTT-based LF- and HF-BP estimation error and BRS, separately. As shown in Figure 7, it is observed that the correlation coefficients in the groups of versus BRS and versus BRS are more central around their mean values than those in the groups of versus BRS and versus BRS. Also, the mean values of the correlation coefficients in LF sections ( versus BRS, versus BRS) are apparently higher than those in HF sections ( versus BRS, versus BRS). In other words, there is a higher correlation between the ABR sensitivity and PTT-based LF-BP estimation error.

Figure 7: Boxplots of the correlation coefficients in LF ( versus BRS, versus BRS) and HF ( versus BRS, versus BRS) sections.

The detailed evaluation criteria for the correlation coefficients listed in Table 4 also indicate the higher correlation between the ABR sensitivity and PTT-based LF-BP estimation error. The mean value ± standard variation of the correlation coefficients in the group of versus BRS and versus BRS is 0.817 ± 0.038 and 0.411 ± 0.145, respectively. In contrast, the correlation coefficients for the group of versus BRS and versus BRS are only 0.403 ± 0.145 and 0.274 ± 0.154, respectively. The MAV value is 0.522 and 0.355 in the LF group of versus BRS and versus BRS, respectively. However, the MAV is only 0.071 and 0.080 in the HF group of versus BRS and versus BRS, respectively. In the t-test of the difference between the LF and HF class, significant differences were found ().

Table 4: The correlation for BRS and BP estimation error in LF and HF sections.

Classes
Groups
Pearson’s correlation coefficients ()

LF
vs BRS
0.817
0.038
0.522

vs BRS
0.757
0.069
0.355

HF
vs BRS
0.403
0.145
0.071

vs BRS
0.274
0.154
0.080

These results clearly suggested that there is a high correlation between the PTT-based LF-BP estimation error and ABR sensitivity. Nevertheless, no correlation has been found between the PTT-based HF-BP estimation error and ABR sensitivity. It is consistent with the results of qualitative analyses in Section 2.4.

5. Conclusion and Further Works

In this paper, we analyzed the poor performance of PTT in LF-BP estimation for the subjects before and after exercise. At first, based on the ABR model, the influence of ABR sensitivity on PTT-based BP estimation was analyzed. It was found that HR may have the unexpected effects on PTT in the LF section under the ABR regulation, which could help to explain the cause behind the failure of PTT in LF-BP estimation. Then, the physiological data from 42 volunteers were collected to verify this hypothesis through correlation analysis between the ABR sensitivity and the estimation error of PTT-based BP in different frequencies. For the overall correlation analysis (Table 3), no obvious correlation was found between the ABR sensitivity and the PTT-based SBP and DBP estimated error. The value is 0.533 ± 0.109 and 0.411 ± 0.145 for the SBP and DBP group, respectively. However, in the comparative correlation analysis (Table 4), remarkable differences () on the correlations of ABR sensitivity and PTT-based BP estimation error have been observed among LF and HF groups. Especially, in the correlation analysis of ABR sensitivity and the PTT-based LF-BP estimation error, very high values of 0.817 ± 0.038 and 0.757 ± 0.069 are observed in the SBP and DBP group, respectively. In contrast, for the correlation of ABR sensitivity and PTT-based HF-BP estimation error, only values of 0.403 ± 0.145 and 0.274 ± 0.154 are viewed in the SBP and DBP group, respectively. Moreover, it is worth noting that the MAV values in the LF group are all higher than those in the HF group.

These experiment results indicated that there is a definitely high correlation between the error of PTT-based LF-BP estimation and the ABR sensitivity for the subjects before and after exercise. To the best of our knowledge, this phenomenon has not been mentioned in the previous studies. This is also the most important finding in this study. It suggested that there is an ABR-related complex LF autonomic regulation mechanism on BP, PTT, and HR, which influences the estimation accuracy of PTT in LF-BP for the subjects before and after exercise. This finding also gave some possible inspirations to help constructing a better PTT-based BP estimation model that: (1) to weaken the influence of ABR, HR is an important parameter that should be considered in improving the accuracy of PTT-based LF-BP estimation; and (2) a weighted frequency-dependent model is required to better estimate BP from PTT. However, it should be noted that there are some limitations in this study. Due to the complexity of cardiovascular system, ABR could not be the only reason for the failure of PTT in LF-BP estimation. The details should be investigated with more different cardiovascular parameters from the subjects under different conditions in further work.

Data Availability

The physiological data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank the volunteers from the School of Instrument Science. This study was supported by the “Fundamental Research Funds for the Central Universities” and “Postgraduate Research & Practice Innovation Program of Jiangsu Province” (KYCX17_0066)

References

L. Peter, N. Noury, and M. Cerny, “A review of methods for non-invasive and continuous blood pressure monitoring: pulse transit time method is promising?” IRBM, vol. 35, no. 5, pp. 271–282, 2014.

G. V. Marie, C. R. Lo, J. Van Jones, and D. W. Johnston, “The relationship between arterial blood pressure and pulse transit time during static and dynamic exercise,” Psychophysiology, vol. 21, no. 5, pp. 521–527, 1984.

R. A. Payne, C. N. Symeonides, D. J. Webb, and S. R. J. Maxwell, “Pulse transit time measured from the ECG: an unreliable marker of beat-to-beat blood pressure,” Journal of Applied Physiology, vol. 100, no. 1, pp. 136–141, 2006.

C. Douniama, C. U. Sauter, and R. Couronne, “Blood pressure tracking capabilities of pulse transit times in different arterial segments: a clinical evaluation,” Computers in Cardiology, vol. 2009, pp. 201–204, 2009.

M. Masè, W. Mattei, R. Cucino, L. Faes, and G. Nollo, “Feasibility of cuff-free measurement of systolic and diastolic arterial blood pressure,” Journal of Electrocardiology, vol. 44, no. 2, pp. 201–207, 2011.

R. Mukkamala, J. Hahn, O. T. Inan, L. K. Mestha, C. Kim, and T. Hakan, “Toward ubiquitous blood pressure monitoring via pulse transit time : theory and practice,” IEEE Transactions on Biomedical Engineering, vol. 62, no. 8, pp. 1879–1901, 2015.

Y. Chen, C. Wen, G. Tao, M. Bi, and G. Li, “Continuous and noninvasive blood pressure measurement: a novel modeling methodology of the relationship between blood pressure and pulse wave velocity,” Annals of Biomedical Engineering, vol. 37, no. 11, pp. 2222–2233, 2009.

Y. Yoon, J. H. Cho, and G. Yoon, “Non-constrained blood pressure monitoring using ECG and PPG for personal healthcare,” Journal of Medical Systems, vol. 33, no. 4, pp. 261–266, 2009.

C. Ahlstrom, A. Johansson, F. Uhlin, T. Länne, and P. Ask, “Noninvasive investigation of blood pressure changes using the pulse wave transit time: a novel approach in the monitoring of hemodialysis patients,” Journal of Artificial Organs, vol. 8, no. 3, pp. 192–197, 2005.

X. R. Ding, Y. T. Zhang, J. Liu, W. X. Dai, and H. K. Tsang, “Continuous cuffless blood pressure estimation using pulse transit time and photoplethysmogram intensity ratio,” IEEE Transactions on Biomedical Engineering, vol. 63, no. 5, pp. 964–972, 2015.

H. Gesche, D. Grosskurth, G. Küchler, and A. Patzak, “Continuous blood pressure measurement by using the pulse transit time: comparison to a cuff-based method,” European Journal of Applied Physiology, vol. 112, no. 1, pp. 309–315, 2012.

E. Chung, G. Chen, B. Alexander, and M. Cannesson, “Non-invasive continuous blood pressure monitoring: a review of current applications,” Frontiers of Medicine, vol. 7, no. 1, pp. 91–101, 2013.

E. Monte-Moreno, “Non-invasive estimate of blood glucose and blood pressure from a photoplethysmograph by means of machine learning techniques,” Artificial Intelligence in Medicine, vol. 53, no. 2, pp. 127–138, 2011.

M. Kachuee, M. M. Kiani, H. Mohammadzade, and M. Shabany, “Cuff-less blood pressure estimation algorithms for continuous health-care monitoring,” IEEE Transactions on Biomedical Engineering, vol. 64, no. 4, pp. 859–869, 2016.

M. Gao, N. B. Olivier, and R. Mukkamala, “Comparison of noninvasive pulse transit time estimates as markers of blood pressure using invasive pulse transit time measurements as a reference,” Physiological Reports, vol. 4, no. 10, Article ID e12768, 2016.

X. R. Ding, N. Zhao, G. Z. Yang et al., “Continuous blood pressure measurement from invasive to unobtrusive: celebration of 200th birth anniversary of carl ludwig,” IEEE Journal of Biomedical and Health Informatics, vol. 20, no. 6, pp. 1455–1465, 2016.

S. L. O. Martin, A. M. Carek, C. S. Kim et al., “Weighing scale-based pulse transit time is a superior marker of blood pressure than conventional pulse arrival time,” Scientific Reports, vol. 6, no. 1, pp. 1–8, 2016.

A. Esmaili, M. Kachuee, and M. Shabany, “Nonlinear cuffless blood pressure estimation of healthy subjects using pulse transit time and arrival time,” IEEE Transactions on Instrumentation and Measurement, vol. 66, no. 12, pp. 3299–3308, 2017.

X. Ding, B. P. Yan, Y. T. Zhang, J. Liu, N. Zhao, and H. K. Tsang, “Pulse transit time based continuous cuffless blood pressure estimation: a new extension and a comprehensive evaluation,” Scientific Reports, vol. 7, no. 1, p. 11554, 2017.

F. Miao, N. Fu, Y. T. Zhang et al., “A novel continuous blood pressure estimation approach based on data mining techniques,” IEEE Journal of Biomedical and Health Informatics, vol. 21, no. 6, pp. 1730–1740, 2017.

A. Rapalis, A. Janušauskas, V. Marozas, and A. Lukoševičius, “Estimation of blood pressure variability during orthostatic test using instantaneous photoplethysmogram frequency and pulse arrival time,” Biomedical Signal Processing and Control, vol. 32, pp. 82–89, 2017.

L. Faes, G. Nollo, and A. Porta, “Mechanisms of causal interaction between short-term RR interval and systolic arterial pressure oscillations during orthostatic challenge,” Journal of Applied Physiology, vol. 114, no. 12, pp. 1657–1667, 2013.

G. Mancia, “Short- and long-term blood pressure variability: present and future,” Hypertension, vol. 60, no. 2, pp. 512–517, 2012.

G. Parati, J. P. Saul, M. Di Rienzo, and G. Mancia, “Spectral analysis of blood pressure and heart rate variability in evaluating cardiovascular regulation a critical appraisal,” Hypertension, vol. 25, no. 6, pp. 1276–1286, 1995.

M. J. Drinnan, J. Allen, and A. Murray, “Relation between heart rate and pulse transit time during paced respiration,” Physiological Measurement, vol. 22, no. 3, pp. 425–432, 2001.

H. T. Ma and Y. T. Zhang, “Spectral analysis of pulse transit time variability and its coherence with other cardiovascular variabilities,” Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 1, pp. 6442–6445, 2006.

Q. Liu, C. C. Y. Poon, and Y. T. Zhang, “Time-frequency analysis of variabilities of heart rate, systolic blood pressure and pulse transit time before and after exercise using the recursive autoregressive model,” Biomedical Signal Processing and Control, vol. 6, no. 4, pp. 364–369, 2011.

S. Cavalcanti, “Arterial baroreflex influence on heart rate variability: a mathematical model-based analysis,” Medical & Biological Engineering & Computing, vol. 38, no. 2, pp. 189–197, 2000.

X. Yang, W. Zhang, D. Wang, and P. Ping, “A calibration method for cuffless continue blood pressure measurement using gaussian normalized pulse transit time,” in Proceedings of 2018 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), IEEE, China, May 2018.

P. Castiglioni, G. Parati, A. Civijian, L. Quintin, and M. Di Rienzo, “Local scale exponents of blood pressure and heart rate variability by detrended fluctuation analysis: Effects of posture, exercise, and aging,” IEEE Transactions on Biomedical Engineering, vol. 56, no. 3, pp. 675–684, 2009.

K. Dragomiretskiy and D. Zosso, “Variational mode decomposition,” IEEE Transactions on Signal Processing, vol. 62, no. 3, pp. 531–544, 2014.

Y. Xu, D. Wang, W. Zhang, P. Ping, and L. Feng, “Detection of ventricular tachycardia and fibrillation using adaptive variational mode decomposition and boosted-CART classifier,” Biomedical Signal Processing and Control, vol. 39, pp. 219–229, 2018.

Ş. Yücelbaş, C. Yücelbaş, G. Tezel, S. Özşen, and Ş. Yosunkaya, “Automatic sleep staging based on SVD, VMD, HHT and morphological features of single-lead ECG signal,” Expert Systems with Applications, vol. 102, pp. 193–206, 2018.

R. K. Tripathy, M. R. A. Paternina, J. G. Arrieta, and P. Pattanaik, “Automated detection of atrial fibrillation ECG signals using two stage VMD and atrial fibrillation diagnosis index,” Journal of Mechanics in Medicine and Biology, vol. 17, no. 7, Article ID 1740044, 2017.

M. I. Bogachev, O. V. Mamontov, A. O. Konradi, Y. D. Uljanitski, J. W. Kantelhardt, and E. V. Schlyakhto, “Analysis of blood pressure-heart rate feedback regulation under non-stationary conditions: beyond baroreflex sensitivity,” Physiological Measurement, vol. 30, no. 7, pp. 631–645, 2009.