# Pulse-Wave-Pattern Classification with a Convolutional Neural Network

Gaoyang Li, Kazuhiro Watanabe, Hitomi Anzai, Xiaorui Song, Aike Qiao, Makoto Ohta

► To cite this version:

Gaoyang Li, Kazuhiro Watanabe, Hitomi Anzai, Xiaorui Song, Aike Qiao, et al. Pulse-Wave-Pattern Classification with a Convolutional Neural Network. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-51334-2. hal-02357984

HAL Id: hal-02357984

https://hal.science/hal-02357984v1

Submitted on 14 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
There are amendments to this paper

OPEN

# Pulse-Wave-Pattern Classification with a Convolutional Neural Network

Gaoyang Li<fnref n="1" />,², Kazuhiro Watanabe<fnref n="1" />,², Hitomi Anzai², Xiaorui Song³, Aike Qiao⁴ & Makoto Ohta⁵*

Owing to the diversity of pulse-wave morphology, pulse-based diagnosis is difficult, especially pulse-wave-pattern classification (PWPC). A powerful method for PWPC is a convolutional neural network (CNN). It outperforms conventional methods in pattern classification due to extracting informative abstraction and features. For previous PWPC criteria, the relationship between pulse and disease types is not clear. In order to improve the clinical practicability, there is a need for a CNN model to find the one-to-one correspondence between pulse pattern and disease categories. In this study, five cardiovascular diseases (CVD) and complications were extracted from medical records as classification criteria to build pulse data set 1. Four physiological parameters closely related to the selected diseases were also extracted as classification criteria to build data set 2. An optimized CNN model with stronger feature extraction capability for pulse signals was proposed, which achieved PWPC with 95% accuracy in data set 1 and 89% accuracy in data set 2. It demonstrated that pulse waves are the result of multiple physiological parameters. There are limitations when using a single physiological parameter to characterise the overall pulse pattern. The proposed CNN model can achieve high accuracy of PWPC while using CVD and complication categories as classification criteria.

Pulse waves contain a large quantity of pathological and physiological information<fnref n="1" />,². Pulse-wave characteristics are closely related to diseases (hypertension, type 2 diabetes, atherosclerosis, etc.), especially cardiovascular diseases (CVD) and physiological parameters [pulse-wave velocity, cardio-ankle vascular index (CAVI), blood pressure, etc.]³,⁴. Therefore, pulse analysis is extensive used in cardiovascular function assessment and non-invasive early diagnosis of cardiovascular disease and related complications⁵. TCPD (Traditional Chinese Pulse Diagnosis) refers to the diagnosis of diseases via traditional Chinese medical practices by feeling the change in pulse at the patient's wrist, which is highly dependent on the doctor's skill and experience⁶. Computer-aided analysis has made some achievements in pulse diagnosis, especially in pulse-wave-pattern classification (PWPC). For example, Wang et al. divided 407 sets of pulse data into five pulse patterns by using a Bayesian network based on six pulse parameters: depth, width, length, frequency, rhythm and strength (84% successful classification rate)⁷. Moreover, Xu et al. divided 320 sets of pulse data into 16 pulse patterns by using a fuzzy neural network based on differences in pulse shapes, widths, positions and some specific local parameters (90% successful classification rate)⁸. However, the diverse morphology of pulse waves remains a difficulty for PWPC, which may lead to problems such as waveform local time shifting, as shown in Fig. 1⁹. In addition, the classification criteria of these studies are based on TCPD theory, which means that a pulse pattern may correspond to a variety of disease categories, as shown in Fig. 1¹⁰. It also leads to a decrease in the clinical practicality of pulse-based diagnosis. Thus, in this study, we selected new classification criteria—that is, the CVD and complication categories and the clinical physiological parameters—with the aim of developing a practical PWPC method with a high classification rate.

With the research and development of deep learning, various of neural network structures have been designed for signal processing. Recurrent neural network (RNN)¹¹, based on its internal memory, is used to process arbitrary time series input sequence such as non-segmented handwriting recognition, speech recognition, etc. Long

<footnote n="1">¹Institute of Fluid Science, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai, Miyagi, 980-8577, Japan. ²Graduate School of Biomedical Engineering, Tohoku University, 6-6 Aramaki-aza-aoba, Aoba-ku, Sendai, Miyagi, 980-8579, Japan. ³Department of Radiology, Taishan Medical University, No.619 Greatwall Road, Daiyue District, Taian, Shandong, 271000, China. ⁴College of Life Science and Bioengineering, Beijing University of Technology, No.100, Pingleyuan, Chaoyang District, Beijing, 100022, China. ⁵ELyTMAx UMI 3757, CNRS-Université de Lyon-Tohoku University, Sendai, Japan. *email: makoto.ohta@tohoku.ac.jp</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** According to previous studies' classification criteria, we show five pulse waves that exhibit a taut pulse pattern, which involves a pulse with a high second peak (local time shifting), as follows: (a) typical taut pulse, (b) taut pulse with high tidal wave and (c) taut pulse with tidal wave merged with percussion wave⁹. With the help of medical doctors, (d) and (e) were extracted from our database. Although (a-e) all feature a taut pulse pattern, there are still differences in some local waveform characteristics. In addition, the subject of (d) suffered from hyperlipidaemia, while the subject of (e) suffered from atherosclerosis. This shows that, under the previous classification criteria, a single pulse pattern might correspond to many disease categories.</figcaption></figure>

short term memory (LSTM)¹², as a variant of RNN, can effectively prevent the occurrence of gradient vanishing from processing time series signals. In recent years, remarkable achievements have been made in the field of pattern classification via the use of convolutional neural networks (CNNs) as deep learning structures¹³⁻¹⁶. CNNs provide an end-to-end learning model. The trained CNNs by the gradient descent method can learn the characteristics of input data and further complete the pattern classification. CNNs have strong ability of feature learning and pattern classification. The main reason is that the features of the lower layers are derived from the partial information and convolution kernel with sharing weights from the upper layer. CNNs have been applied in the classification of human physiological signal patterns. Based on a 34-layer CNN, Rajpurkar et al. classified the electrocardiogram (ECG) signals into 14 types¹⁷. Moreover, Rubin et al. performed heart-sound recordings based on deep CNN and Mel-frequency cepstral coefficients¹⁸. These studies used CNN to achieve pattern classification of relevant physiological signals and achieved higher accuracy than the diagnostic results of experienced physicians. Furthermore, Hu et al. used CNN to divide pulse waves into two types: health and subhealth¹⁹. In the present study, in view of the large amount of pathological and physiological information contained in pulse signals, we collected the required data under the guidance of the doctor and established two data sets based on either CVD/ complication categories or physiological parameters. We proposed an optimised CNN model for PWPC based on these two data sets. The purpose of this study was to identify a practical and efficient classification criterion for PWPC based on CNN, which contributed to non-invasive, practical and effective diagnosis of CVDs and related complications.

## Results

The average pulse waves of each pattern in the two data sets are shown in Fig. 2.

We showed the learning curves of data set 1 and data set 2 respectively to evaluate their PWPC performance on the proposed CNN model, as shown in Fig. 3. For cost-value curve, the decline rate of data set 1 was significantly higher than that of data set 2. For training error and test error, the minimum value of data set 1 was smaller than that of data set 2. Especially test error, data set 1 (When epoch was 90, the minimum test error was 0.08. Epoch was the number of iterations in CNN pattern classification) was much smaller than data set 2 (When epoch was 100, the minimum test error was 0.34). With the same proposed CNN, the six pulse patterns in data set 1 showed higher calculation efficiency and feature expression ability than those five patterns in data set 2.
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Average pulse-wave patterns in data set 1 (a) and data set 2 (b). Abbreviation Represents Pulse Wave Patterns: H1-Healthy Control Group in Data Set 1; Hn- hypertension; At-atherosclerosis; Ha-hyperlipidaemia; Td-type 2 diabetes; HCA-Hypertension complicated by atherosclerosis; H2-Healthy Control Group in Data Set 2; BP-blood pressure; baPWV- brachial-ankle pulse wave velocity; BV-blood viscosity.</figcaption></figure>

Table 1 shows the overall values of the evaluation parameters in the two data sets. The accuracy and other evaluation parameters of PWPC in data set 1 (overall accuracy = 0.95) were higher than those in data set 2 (overall accuracy = 0.89). Tables 2 and 3 show the details for each pattern in the two data sets separately. Pulse-wave patterns representing healthy subjects (H1 and H2) could be identified with high precision (precision H1 = 1, recall H1 = 0.99; precision H2 = 0.97, recall H2 = 0.97). The HCA, as the pulse pattern of complications, had the lowest classification rate in data set 1 (precision HCA = 0.89, recall HCA = 0.91). In addition, the classification performance of other pulse patterns in data set 1 was higher than that in data set 2.

To further assess the PWPC result of the proposed CNN model, the two data sets were put into different neural networks models for PWPC. Table 4 shows the accuracy of PWPC with those different models. It details network methods, classification criteria, number of subjects, and the accuracy. It demonstrated that compared with other neural networks or other CNN structures, our proposed CNN model achieved higher accuracy in PWPC under the new classification criteria, which also meant stronger feature extraction ability for pulse signals.

To further analyse the causes of errors in pattern classification, we determined the confusion matrix of the two data sets, as shown in Fig. 4. The cause of errors in data set 1 was mainly the erroneous classification of the four pulse patterns of Hn, At, HCA and Td. In data set 2, with the exception of the control group (H2), the remaining four pulse patterns (BP, CAVI, baPWV and BV) were found to interfere with each other and have higher error rates.
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>Figure 3. Learning curve in data set 1 (a) and data set 2 (b).</figcaption></figure>

<table><caption>Table 1. PWPC evaluation of per pulse patterns in two data sets.</caption><thead><tr><td>Data set</td><td>Overall accuracy</td><td>Overall precision</td><td>Overall recall</td><td>Overall F-measure</td></tr></thead><tbody><tr><td>Data set 1</td><td>0.95</td><td>0.95</td><td>0.95</td><td>0.95</td></tr><tr><td>Data set 2</td><td>0.89</td><td>0.89</td><td>0.89</td><td>0.89</td></tr></tbody></table>

<table><caption>Table 2. PWPC evaluation of per pulse patterns in data set 1.</caption><thead><tr><td>Pulse pattern</td><td>Precision</td><td>Recall</td><td>F-measure</td></tr></thead><tbody><tr><td>H1</td><td>1</td><td>0.99</td><td>0.99</td></tr><tr><td>Hn</td><td>0.94</td><td>0.93</td><td>0.94</td></tr><tr><td>At</td><td>0.90</td><td>0.94</td><td>0.92</td></tr><tr><td>Ha</td><td>1</td><td>0.99</td><td>0.99</td></tr><tr><td>Td</td><td>0.96</td><td>0.93</td><td>0.94</td></tr><tr><td>HCA</td><td>0.89</td><td>0.91</td><td>0.90</td></tr></tbody></table>

Discussion

In this study, CVD and associated complications as well as related physiological parameters were extracted, which were used as classification criteria. According to the new classification criteria, we screened the subjects' pulse waves and created data set 1 and data set 2, respectively. An optimised CNN model was proposed for PWPC. It
<table><caption>**Table 3.** PWPC evaluation of per pulse patterns in data set 2.</caption><thead><tr><td>Pulse pattern</td><td>Precision</td><td>Recall</td><td>F-measure</td></tr></thead><tbody><tr><td>H2</td><td>0.97</td><td>0.97</td><td>0.97</td></tr><tr><td>BP</td><td>0.89</td><td>0.89</td><td>0.89</td></tr><tr><td>CAVI</td><td>0.82</td><td>0.84</td><td>0.83</td></tr><tr><td>baPWV</td><td>0.84</td><td>0.87</td><td>0.85</td></tr><tr><td>BV</td><td>0.95</td><td>0.89</td><td>0.92</td></tr></tbody></table>

<table><caption>**Table 4.** PWPC accuracy of different methods.</caption><thead><tr><td>Network</td><td>Method</td><td>Classification criteria</td><td>Number of subjects</td><td>Accuracy</td></tr></thead><tbody><tr><td rowspan="2">The proposed CNN model</td><td>CNN</td><td>CVD and complications</td><td>412</td><td>0.95</td></tr><tr><td>CNN</td><td>Physiological parameters</td><td>412</td><td>0.89</td></tr><tr><td rowspan="2">LetNet<sup>38</sup></td><td>CNN</td><td>CVD and complications</td><td>412</td><td>0.69</td></tr><tr><td>CNN</td><td>Physiological parameters</td><td>412</td><td>0.63</td></tr><tr><td rowspan="2">AlexNet<sup>14</sup></td><td>CNN</td><td>CVD and complications</td><td>412</td><td>0.73</td></tr><tr><td>CNN</td><td>Physiological parameters</td><td>412</td><td>0.70</td></tr><tr><td rowspan="2">VGG-Net<sup>15</sup></td><td>CNN</td><td>CVD and complications</td><td>412</td><td>0.81</td></tr><tr><td>CNN</td><td>Physiological parameters</td><td>412</td><td>0.79</td></tr><tr><td>Wang's network<sup>7</sup></td><td>Bayesian Network</td><td>Based on TCPD</td><td>407</td><td>0.84</td></tr><tr><td>Xu's network<sup>8</sup></td><td>Fuzzy Neural Network</td><td>Based on TCPD</td><td>320</td><td>0.90</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** The confusion matrices of data set 1 (a) and data set 2 (b). The confusion matrix is an intuitive method for evaluating the results of pattern classification CNN models. The real categories (rows) and predicted categories (columns) of the classification results can be read directly. For example, in matrix (a), there were 70 (65 + 5) pulse waves which really belonged to the Hn pattern (the second row), while the CNN model predicted 69 (65 + 1 + 1 + 2) pulse waves in the Hn pattern (the second column).</figcaption></figure>

achieved the classification of six pulse patterns in data set 1 with an accuracy of 95% and the classification of six pulse patterns in data set 2 with an accuracy of 89%. The main contributions of this study are as follows:

1. Two pulse wave data sets were created, which contained a large amount of physiological and pathological information of subjects.

2. New classification criteria and optimized CNN model were proposed, which achieves higher accuracy than previous studies<sup>7,8,19-21</sup>.

This study demonstrates that CVD and complications are practical and efficient classification criteria, enabling the optimised CNN model to achieve high accuracy for PWPC.

We observed that the classification errors in data set 1 were mainly due to the erroneous classification of the Hn, At and HCA patterns. This was due to the simultaneous occurrence of hypertension and atherosclerosis on behalf of HCA. There must be some similar pulse characteristics between HCA and the other two diseases, which
<figure><img src="image_7.png" /><figcaption>**Figure 5.** The pulse-wave form baPWV pulse pattern with the different baPWV and different CAVI. Pulse waves from six subjects were selected.</figcaption></figure>

indicates that, in order to ensure that the characteristics of the different pulse patterns are typical, the selected data specimen must exclude the effect of complications at the same time. In addition, in data set 1, Td was also partially misclassified as Hn (n = 1), At (n = 3) and HCA (n = 1). Previous studies showed that type 2 diabetes could increase the risk and mortality of CVD, and they had similarities in the damage to the cardiovascular system<sup>22-24</sup>. Thus, there might have been similar pulse waveform characteristics between Td and Hn, At, HCA patterns, which led to classification errors.

In data set 2, four pulse patterns (BP, CAVI, baPWV and BV) were found to interfere with each other in pattern classification. Previous studies showed that the effect of a single physiological parameter on pulse waveform was mainly reflected in the change of some local characteristics<sup>25-27</sup>. The pulse waveform characteristics with the same value of one specific physiological parameter would change as a result of the differences of other physiological parameters, as shown in Fig. 5. It may have led to the errors of pattern classification in data set 2. Our study showed that the pulse-wave was the result of multiple physiological parameters. There are clearly limitations associated with using a single physiological parameter in characterising the overall pulse pattern. Disease was the result of multiple physiological parameters, which might explain the higher classification accuracy in data set 1.

This study had several limitations. The most important one was the relatively limited number of subjects. Limited by the number of subjects, the effects of some physiological information such as age, height and weight on pulse waveform were ignored, which inevitably led to errors in pattern classification<sup>1</sup>. However, in our study, the number of pulse waves in each pulse pattern was several times that in some previous studies<sup>21,28</sup>. To some extent, the findings indicated that each of our patterns could represent the typical pulse characteristics. In addition, this study focused on the classification criteria of pulse patterns. For this purpose, we used the same CNN model to classify two data sets. Regarding the low classification rate of data set 2, we did not explore whether it could be improved by optimising the architecture of the CNN model.

## Conclusions

In this study, we established pulse wave data set 1 and data set 2 based on the classification criteria: CVD categories and related physiological parameters. CNN was used to extract features from two data sets and to achieve PWPC with high accuracy. The main contribution of this study is to propose the new classification criteria for PWPC and construct a matching CNN model. The optimized CNN model achieved PWPC with 95% accuracy in data set 1 and 89% accuracy in data set 2. This study demonstrated that pulse waves are the result of multiple physiological parameters, so there are limitations when using a single physiological parameter to characterize the overall pulse pattern. The proposed CNN model can achieve high accuracy PWPC while using CVD and complication categories as classification criteria, which contributes to non-invasive, practical and effective diagnosis of CVD and associated complications.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** The process of screening the subjects in data set 1. a,b,c,d Screening criteria: The number of subjects for a selected disease should be more than 20. The disease or complications must be of the five types selected in this study. There is no serious abnormality in pulse waves caused by noise or incorrect data collection, among others. We show all cases and numbers of excluded subjects in c: type 2 diabetes complicated by hypertension (n = 4), type 2 diabetes complicated by atherosclerosis (n = 3), type 2 diabetes complicated by heart failure (n = 5) and diabetic foot disease (n = 8). Based on the screening criteria, we excluded these cases.</figcaption></figure>

**Method**

**Data collection.** The original pulse wave data were from the "Study on Evaluation Method of Cardiovascular System Based on Non-invasive Detection of Blood Pressure and Pulse-Wave of Limbs<sup>29</sup>", which recruited 412 subjects and determined their physiological parameters and more than 12,000 cycles of pulse waves. The pulse and blood pressure signal measuring device was Fukuda VS-1500A. In addition, the subjects' brachial ankle pulse-wave velocity (baPWV) and blood viscosity were collected. All subjects were registered at Beijing University of Technology Hospital, and information on their diseases was collected through the subjects' medical records.

The study with its experimental protocols and relevant details was approved by the Institutional Ethics Committee of Beijing University of Technology and Tohoku University. All experiments were performed in accordance with relevant guidelines and regulations. We explained the content of the study to the subjects in detail, and on this basis, the subjects signed the informed consent form.

**Pulse waveform denoising and normalisation.** In this study, we collected the pulse signals from the wrist of the subjects. The denoising and normalization of pulse signals were processed with the same method as the previous studies³⁰. Firstly, the noise was removed with wavelet transform decomposition method³¹. Then, in order to prevent the distortion of pulse signals, according to Nyquist theorem and actual sampling frequency⁸,¹⁹, the sampling points of single cycle of pulse wave were set at 200. Because the focus of this study was the change of pulse wave model, the amplitude of pulse wave was normalized to 0–200 in each cycle.

**Data sets.** Previous studies classified pulses into patterns based on the TCPD theory<sup>7–10,32</sup>. However, as mentioned previously, under this classification criterion, one pulse pattern may correspond to a variety of disease categories. Thus, in this study, based on subjects' clinical data, we directly selected five diseases as new classification criteria: hypertension, atherosclerosis, hyperlipidaemia, type 2 diabetes and hypertension complicated by atherosclerosis (HCA). Type 2 diabetes, as one of the common complications of CVD<sup>33</sup>, and HCA were used to study the effects of CVD complications on pulse waves. To ensure the typical characteristics of each pulse pattern, the pulse signals from subjects who only suffered from one of the five diseases and healthy subjects (a total of six types) were used as new pulse patterns to build data set 1, as shown in Fig. 6.

We simultaneously selected four physiological parameters closely related to the selected diseases as classification criteria: blood pressure, which can be used as an indicator for assessing hypertension³⁴; cardio-ankle vascular index (CAVI), which is one of the indicators for assessing atherosclerosis³⁵; and brachial ankle pulse-wave velocity (baPWV), which can be used as an indicator for evaluating cardiovascular function in type 2 diabetics³⁶; For patients with hyperlipidaemia, an increase of blood lipids often occurs simultaneously with increased blood viscosity³⁷. Based on the subjects in data set 2 and the medical reference range, we determined the range of each physiological parameter. The pulse waves of subjects in whom only one of the four parameters was beyond the range were selected. The pulse waves of subjects whose four parameters were all within the range were also selected as a healthy control group. Then the five types of pulse pattern were used to build data set 2, as shown in Fig. 7.

For the processing of pulse image, this study used the same method as previous studies³⁰. We extracted the pulse cycles from the selected subjects. To avoid data duplication affecting the accuracy of CNN prediction, all
<figure><img src="image_9.png" /><figcaption>**Figure 7.** The screening process of the subjects in data set 2. a,b,c,d Screening criteria: The number of subjects for selected parameters should be more than 20. Subjects' other parameters, such as stroke output and cardiac output, must be within the normal range of medical reference. There is no serious abnormality in pulse-wave caused by noise or incorrect data collection, among others. For a,b,c,d, most of the excluded subjects had three or even four parameter values outside of the range. To ensure that the characteristics of each pulse pattern were typical, we excluded these subjects.</figcaption></figure>

<table><caption>**Table 5.** The details of PWPC data sets.</caption><thead><tr><th></th><th>Pulse categories</th><th>Total pulse number</th><th>Training number</th><th>Test number</th><th>Type of disease or physiological parameters (range)</th></tr></thead><tbody><tr><td rowspan="6">Data set 1</td><td>H1</td><td>210</td><td>140</td><td>70</td><td>Healthy control group in data set 1</td></tr><tr><td>Hn</td><td>210</td><td>140</td><td>70</td><td>Only hypertension</td></tr><tr><td>At</td><td>210</td><td>140</td><td>70</td><td>Only atherosclerosis</td></tr><tr><td>Ha</td><td>210</td><td>140</td><td>70</td><td>Only hyperlipidaemia</td></tr><tr><td>Td</td><td>210</td><td>140</td><td>70</td><td>Only type 2 diabetes</td></tr><tr><td>HCA</td><td>210</td><td>140</td><td>70</td><td>Hypertension complicated by atherosclerosis</td></tr><tr><td rowspan="6">Data set 2</td><td>H2</td><td>210</td><td>140</td><td>70</td><td>Healthy control group in data set 2</td></tr><tr><td>BP</td><td>210</td><td>140</td><td>70</td><td>Only high blood pressure (&gt;90/140 mmHg)</td></tr><tr><td>CAVI</td><td>210</td><td>140</td><td>70</td><td>Only CAVI (&gt;9.0)</td></tr><tr><td>baPWV</td><td>210</td><td>140</td><td>70</td><td>Only baPWV (&gt;1400 cm/s)</td></tr><tr><td>BV</td><td>210</td><td>140</td><td>70</td><td>Only blood viscosity (&gt;5.0)</td></tr></tbody></table>

pulse waves in the two data sets were taken from different cycles. The total cycles of each pulse pattern were 210, which were divided into training set and test set, as shown in Table 5. As mentioned above, the number of sampling points in a single cycle of normalized pulse wave was 200, and the amplitude was 0–200. Therefore, the pulse wave signals were processed as input PNG pulse images with a size of 200 × 200 pixels.

**The proposed CNN.** In this study, an optimised CNN model (10-layer) was proposed based on DCNN¹⁹ and LeNet-5³⁸, which had been applied for PWPC, as shown in Fig. 8. Compared with the previous networks, we added dropout³⁹ between the third max pooling layer and the fully connected layer. When CVDs were used as classification criteria, each pulse pattern changed from local waveform difference under previous criteria to overall pulse waveform difference. This led to too many characteristic parameters of pulse wave extracted by CNN, which further led to over-fitting in the training process. Pre-experimental results showed that dropout layer could help reduce test errors and avoid over-fitting phenomenon in the training process (see Supplementary Fig. S1). In addition, the final Softmax activation produced a distribution over the output probability classes for each pulse pattern of two data sets. Besides the layers mentioned above, the CNN also included three convolution layers, three max pooling layers and two fully connected layers. The number of convolution layers was determined by the number of pulse wave characteristic. The insufficient layers led to the inadequate feature extraction ability of CNN, while the excessive layers increased the time cost and calculation cost. In this study, we determined the number of layers by pre-experimental results. The convolutional layers were used to extract complex parameters
<figure><img src="image_10.png" /><figcaption>**Figure 8.** An illustration of the CNN architecture. The size settings of convolution kernels and feature maps are shown in the figure.</figcaption></figure>

of the input feature maps by convolution with kernels. The max pooling layers achieved the down-sampling of the input signals by choosing the maximum value of the area as the value of the pooled area. The max pooling layers could retain the main features of the input signals while reducing the parameters and computation, which helped to avoid the occurrence of over-fitting and improve the generalization ability of the CNN model<sup>40</sup>. The final two fully connected layers combined all of the upper feature maps into a one-dimensional array, which was used to classify the output. In this study, we used the Adam optimiser, which is straightforward to implement, with high calculation efficiency and low memory requirements<sup>41</sup>. In accordance with previous studies and a preliminary experiment, the parameters of the Adam optimiser were as follows: learning rate = 0.001, ε = 0.001, ρ1 = 0.9, ρ2 = 0.999 and δ = 1E<sup>-8</sup>. During the optimisation process, we saved the best model configuration as evaluated on the test set. The CNN was trained by neural_network_console (Sony Company) on an Intel(R) HD Graphics 630 with batch size 64 for 100 epochs.

**Evaluation.** The proposed CNN was evaluated with the average of the operating parameters calculated over time. The overall accuracy, precision, recall and F-measure were determined to assess the classification performance of the network, as presented in the results section. To further evaluate the classification performance of each pulse pattern, we also present the evaluation parameters of each pattern and the confusion matrices for the two test sets. The evaluation parameters were calculated using the true positive (TP), true negative (TN), false positive (FP) and false negative (FN).

In order to further evaluate the PWPC capability of the CNN model proposed in this study, we selected three different neural networks (LetNet<sup>38</sup>, AlexNet<sup>14</sup>, VGG-Net<sup>15</sup>). Data set 1 and data set 2 were used as inputs of these three networks respectively. The PWPC results were compared with the CNN model proposed in this study.

**Data availability**

The datasets analysed during the current study are available from the corresponding author on reasonable request.

Received: 23 October 2018; Accepted: 24 September 2019;

Published online: 17 October 2019

**References**

1. O'Rourke, M. F., Pauca, A. & Jiang, X.-J. Pulse wave analysis. *Br. J. Clin. Pharmacol.* **51**, 507–522 (2001).

2. Korpas, D., Hálek, J. & Doležal, L. Parameters describing the pulse wave. *Physiol. Res.* **58**, 473–479 (2009).

3. Safar, M. E., Levy, B. I. & Struijker-Boudier, H. Current perspectives on arterial stiffness and pulse pressure in hypertension and cardiovascular diseases. *Circulation* **107**, 2864–2869 (2003).

4. Yamashina, A. et al. Brachial-ankle pulse wave velocity as a marker of atherosclerotic vascular damage and cardiovascular risk. *Hypertens. Res.* **26**, 615–622 (2003).

5. Cohn, J. N. et al. Noninvasive pulse wave analysis for the early detection of vascular disease. *Hypertension* **26**, 503–508 (1995).

6. Shu, J. J. & Sun, Y. Developing classification indices for Chinese pulse diagnosis. *Complement. Ther. Med.* **15**, 190–198 (2007).

7. Wang, H. & Cheng A quantitative system for pulse diagnosis in Traditional Chinese Medicine. *Conf. Proc. IEEE Eng. Med. Biol. Soc.* **6**, 5676–5679 (2005).

8. Xu, L., Meng, Q. H., Wang, K., Lu, W. & Li, N. Pulse images recognition using fuzzy neural network. *Expert Syst. Appl.* **36**, 3805–3811 (2009).

9. Zuo, W., Zhang, D., Zhang, D., Zhang, H. & Li, N. Classification of pulse waveforms using edit distance with real penalty. *EURASIP J. Adv. Signal Process.* **2010**, (2010).

10. Xu, L. S., Wang, K. Q. & Wang, L. Pulse Waveforms Classification Based on Wavelet Network. *IEEE EMBS.* 4596–4599 (2005).

11. Graves, A., Mohamed, A. & Hinton, G. Speech recognition with deep recurrent neural networks. *2013 IEEE International Conference on Acoustics, Speech and Signal Processing 6645–6649 (2013)*

12. Greff, K., Srivastava, R. K., Koutnik, J., Steunebrink, B. R. & Schmidhuber, J. LSTM: A Search Space Odyssey. *IEEE Trans. Neural Networks Learn. Syst.* **28**, 2222–2232 (2017).

13. Huang, G. B., Lee, H. & Learned-Miller, E. B. T.-I. C. on C. V. and P. R. Learning hierarchical representations for face verification with convolutional deep belief networks. In 2518–2525 (2012).
14. Krizhevsky, A., Sutskever, I. & Hinton, G. E. B. T.-I. C. on N. I. P. S. ImageNet classification with deep convolutional neural networks. NIPS. 1097–1105 (2012).

15. Simonyan, K. & Zisserman, A. Very Deep Convolutional Networks for Large-Scale Image Recognition. *Comput. Sci.* (2014).

16. zhang, J., liu, P., Zhang, F. & Song, Q. CloudNet: Ground-based Cloud Classification with Deep Convolutional Neural Network. *Geophys. Res. Lett.* 1–8, https://doi.org/10.1029/2018GL077787 (2018).

17. Rajpurkar, P., Hannun, A. Y., Haghpanahi, M., Bourn, C. & Ng, A. Y. Cardiologist-Level Arrhythmia Detection with Convolutional Neural Networks. (2017).

18. Rubin, J. et al. Classifying heart sound recordings using deep convolutional neural networks and mel-frequency cepstral coefficients. *Computing in cardiology conference (CinC)*. 813–816 (2016).

19. Hu, X., Zhu, H., Xu, J., Xu, D. & Dong. Wrist pulse signals analysis based on Deep Convolutional Neural Networks. 2014 IEEE Conference on, J. B. T.-C. I. in B. and C. B. In 1–7 (2014).

20. Bo, I. et al. Cardiovascular Morbidity and Mortality Associated With the Metabolic Syndrome. *Diabetes Care* **24**, 683–689 (2001).

21. Yang, W., Zhang, L. & Zhang, D. Wrist-Pulse Signal Diagnosis Using ICPulse. *Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (EMBS '05)*. 5676–5679, (2009).

22. Kannel, W. B. & Mcgee, D. L. Diabetes and Cardiovascular Disease: The Framingham Study. *Am. J. Cardiol.* **35**, 2035–2038 (1979).

23. Grøntved, A. & Hu, F. B. Television Viewing and Risk of Type 2 Diabetes, Cardiovascular Disease, and All-Cause Mortality A Meta-analysis. *Jama* **305**, 2448–2455 (2011).

24. Wilmot, E. G. et al. Sedentary time in adults and the association with diabetes, cardiovascular disease and death: systematic review and meta-analysis. *Diabetologia* **55**, 2895–2905 (2012).

25. Weber, T. et al. Wave reflections, assessed with a novel method for pulse wave separation, are associated with end-organ damage and clinical outcomes. *Hypertension* **60**, 534–541 (2012).

26. Hirata, K., Kawakami, M. & O'Rourke, M. F. Pulse Wave Analysis and Pulse Wave Velocity. *Circ. J.* **70**, 1231–1239 (2006).

27. Wang, A. et al. Quantification of radial arterial pulse characteristics change during exercise and recovery. *J. Physiol. Sci.* **68**, 1–8 (2018).

28. Zhang, S. R. & Sun Consumer and Control, Q. F. B. T.-I. S. on C. Human Pulse Recognition Based on Convolutional Neural Networks. In 366–369 (2016).

29. Song, X., Li, G., Qiao, A. & Chen, Z. Association of simultaneously measured four-limb blood pressures with cardiovascular function: a cross-sectional study. *Biomed. Eng. Online* **15**, 247–260 (2016).

30. Li, G., Anzai, H., Watanabe, K., Qiao, A. & Ohta, M. Pulse pattern classification of atherosclerotic patients based on convolutional neural network. *The Proceedings of the JSME Conference on Frontiers in Bioengineering 2018.29*, 2C11 (2018).

31. Chang, F., Hong, W., Zhang, T., Jing, J. & Liu, X. Research on Wavelet Denoising for Pulse Signal Based on Improved Wavelet Thresholding. In *2010 First International Conference on Pervasive Computing, Signal Processing and Applications* 564–567 https://doi.org/10.1109/PCSPA.2010.142 (2010).

32. Xu, L. et al. Objectifying Researches on Traditional Chinese Pulse Diagnosis. *Informatica Medica Slovenica*. **2**, 56–63 (2003).

33. Gaede, P. et al. Multifactorial Intervention and Cardiovascular Disease in Patients with Type 2 Diabetes. *N. Engl. J. Med.* **348**, 383–393 (2003).

34. Pessina, A. C. et al. Evaluation of hypertension and related target organ damage by average day-time blood pressure. *Clin. Exp. Hypertens.* **7**, 267 (1985).

35. Okura, T. et al. Relationship between Cardio-Ankle Vascular Index (CAVI) and Carotid Atherosclerosis in Patients with Essential Hypertension. *Hypertens. Res.* **30**, 335–340 (2007).

36. Arnett, D. K., Evans, G. W. & Riley, W. A. Arterial Stiffness: A New Cardiovascular Risk Factor? *Am. J. Epidemiol.* **140**, 669–682 (1994).

37. Irace, C. et al. Influence of blood lipids on plasma and blood viscosity. *Clin. Hemorheol. Microcirc.* **57**, 267–274 (2013).

38. Lécun, Y., Bottou, L., Bengio, Y. & Haffner, P. Gradient-based learning applied to document recognition. *Proc. IEEE* **86**, 2278–2324 (1998).

39. Srivastava, N., Hinton, G., Krizhevsky, A., Sutskever, I. & Salakhutdinov, R. Dropout: a simple way to prevent neural networks from overfitting. *J. Mach. Learn. Res.* **15**, 1929–1958 (2014).

40. Microbiana, B. et al. Lung Pattern Classification for Interstitial Lung Diseases Using a Deep Convolutional Neural Network. *IEEE Trans. Med. Imaging* **35**, 1207–1216 (2016).

41. Kingma, D. P. & Ba, J. Adam: A Method for Stochastic Optimization. *Comput. Sci.* https://doi.org/10.1063/1.4902458 (2014).

## Acknowledgements

This research is partially supported by the Creation of a development platform for implantable/wearable medical devices by a novel physiological data integration system of the Program on Open Innovation Platform with Enterprises, Research Institute and Academia (OPERA) from the Japan Science and Technology Agency (JST). This work was also supported by the Grant-in-Aid [A] (No16H01805), the ImPACT, (JST) "Bionic Humanoids Propelling New Industrial Revolution", and the AMED under Grant Number 18he1802004h0002 projects.

## Author contributions

G.L., X.S. and A.Q. acquired the data. G.L. and M.O. conceived and designed this research; G.L., K.W. and H.A. performed the experiments and analysed the data; All of the authors discussed and co-authored the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-51334-2.

Correspondence and requests for materials should be addressed to M.O.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019