Intelligent Wearable Wrist Pulse Detection System Based on Piezoelectric Sensor Array

The human radial artery pulse carries a rich array of biomedical information. Accurate detection of pulse signal waveform and the identification of the corresponding pulse condition are helpful in understanding the health status of the human body. In the process of pulse detection, there are some problems, such as inaccurate location of radial artery key points, poor signal noise reduction effect and low accuracy of pulse recognition. In this system, the pulse signal waveform is collected by the main control circuit and the new piezoelectric sensor array combined with the wearable wristband, creating the hardware circuit. The key points of radial artery are located by an adaptive pulse finding algorithm. The pulse signal is denoised by wavelet transform, iterative sliding window and prediction reconstruction algorithm. The slippery pulse and the normal pulse are recognized by feature extraction and classification algorithm, so as to analyze the health status of the human body. The system has accurate pulse positioning, good noise reduction effect, and the accuracy of intelligent analysis is up to 98.4%, which can meet the needs of family health care.


Introduction
The heart shoots blood into the artery and transmits it to the whole body through periodic and regular contraction and relaxation. After the blood is exchanged with tissues, it returns to the heart from the veins, forming a complete blood circulation cycle; in medicine, we call this the circulatory system. The pulse waveform formed by blood impact on the radial artery wall has periodic characteristics and carries biomedical information [1]. Because the radial artery of the wrist of the upper limb has the advantages of being close to the body surface, convenient to measure, and less affected by external factors, traditional Chinese medicine (TCM) looks there for Cun, Guan, and Chi; according to the academic language of pulse, "Take three parts of pulse, and each part is one cun" [2]. Namely, there are three key points of the radial artery: the radial protruding point is Guan, one cun above the Guan pulse is Cun, and one cun below the Guan pulse is Chi. Different forces are applied at key points to make the pulse in three states of floating, middle and sinking [3], and the amplitude, frequency and shape of the pulse signal are sensed, so as to discover, dig and extract the pulse condition information hidden behind the pulse signal [4]. According to the internal relationship between human health state and pulse condition, the human health state is analyzed to achieve the purpose of disease screening [5][6][7]. Although this method can quickly, safely and non-invasively realize the patient's health detection and disease diagnosis, it is greatly affected by subjective factors and can find it difficult to accurately identify the pulse condition, which affects the analysis of human health status.
With the development of sensor technology and the rise of artificial intelligence, modern medicine tends to analyze human health status by collecting pulse signal waveforms

System Constitutes
The hardware of this system is composed of the main control circuit, the new piezoelectric sensor array, and the wrist pressure device, as shown in Figure 1. The main control circuit of the system takes STM32F103C8T6 chip as the core, including four 12-bit AD acquisition circuits, XGZP6847A air pressure sensor, SC300IPN micro air pump and EVS-HD-3003 solenoid valve composed of air pressure detection and pressure circuit pressure detection and pressure circuits, laser positioning module composed of an ADNS-9800 chip, power management circuit composed of FM3209 chip, indicator circuit, HC-05 Bluetooth communication circuit, which realizes analog to digital signal conversion, wrist band pressure control, sensor displacement measurement, and Bluetooth data communication. The wrist pressure device is composed of an inflatable wristband, an air bag, radial stem bone ferrules, air filling and discharging ports, and corresponding hoses, which has the characteristics of small volume, fitting wrist skin and good air tightness. The new piezoelectric sensor array is composed of three 11 mm diameter cylindrical silica gel packaged new piezoelectric sensors arranged at an interval of 15 mm and integrated on a flexible circuit board, it can simulate the person's index finger, middle finger and ring finger in the close state, which can effectively cover the corresponding position of Cun, Guan and Chi, so as to realize the acquisition and amplification of weak pulse signal.

New Piezoelectricity Sensor
The new piezoelectric sensor is composed of copper sheet, piezoelectric ceramics (PZT), silver chloride sheet combined with differential charge amplifier circuit, voltage amplifier circuit and filter circuit. When the inner metal surface of the new piezoelectric sensor is under pressure, the piezoelectric ceramics (PZT) between the copper sheet and the silver chloride sheet is deformed. Due to the piezoelectric effect, the surface of the

New Piezoelectricity Sensor
The new piezoelectric sensor is composed of copper sheet, piezoelectric ceramics (PZT), silver chloride sheet combined with differential charge amplifier circuit, voltage amplifier circuit and filter circuit. When the inner metal surface of the new piezoelectric sensor is under pressure, the piezoelectric ceramics (PZT) between the copper sheet and the silver chloride sheet is deformed. Due to the piezoelectric effect, the surface of the material forms a large number of ordered and isotropic charges, and the two ends of the metal are equivalent to a charged capacitor. After the U1 differential charge amplifier circuit, the U2 and U3 voltage amplifier circuits, and the U4 filter circuit, the composite material can identify 0 pC to 1 pC level charge transformation (corresponding to pressures in the range of 0 g to 2 g) and output with 0 V to 3.3 V working voltage. The experimental setup consists of MW MPS-3003C DC source, Tektronix MSO 2012 oscilloscope, LUILEC DC electric pushrod motor, motor operation handle, pressure probe, gasket, KFS-A electronic scale, and the new piezoelectric sensor. The constant speed motor push rod is controlled by the handle, and the pressure is applied to the sensor under the gasket by the pressure probe, and then the pressure is read by the electronic scale, and the oscilloscope shows the waveform. After that, the handle controls the constant speed motor push rod reset. The experiment is repeated to obtain the relationship between the pressure and the sensor output voltage. After the pressure test in the range of 0 g to 2 g, the output voltage of the sensor is positively correlated with the force of the piezoelectric material. When the external force disappears, the material recovers its original shape, the isotropic charge on the material surface disappears, and the sensor output voltage is zero. The piezoelectric effect of the new piezoelectric sensor is shown as follows: where K is a constant, Q is the amount of charge, C 1 is the feedback capacitance, u is the output voltage, and m represents the external pressure. The frequency response range of the sensor is 0.5 Hz to 100 Hz, that is, the noise below 0.5 Hz and the high-frequency noise above 100 Hz are filtered through the circuit design, and the effective information of the pulse signal is retained to the maximum extent. The pressure response range is less than 0.206 kPa, that is, the sensor force range is 0 g to 2 g, the sensor surface area is 0.3025π × 10 −4 m 2 , the gravity acceleration g = 9.8 m/s 2 , according to the pressure calculation formula P = F/S = mg/S ≤ 0.206 kPa. The sensitivity is 16.019 mV/Pa, that is, when the sensor range is 3.3 V, the maximum pressure is 0.206 kPa, so the sensor sensitivity is 3300 mV/206 Pa = 16.019 mV/Pa. The error is less than 5%. The sensor has the characteristics of strong anti-interference ability, high sensitivity, wide pressure range, good repeatability and strong consistency, which can be used for accurate detection of pulse signal waveform. The working principle of the new piezoelectric sensor is shown in Figure 2: is less than 0.206 kPa, that is, the sensor force range is 0 g to 2 g, the sensor surface area is 0.3025 × 10 −4 m 2 , the gravity acceleration g = 9.8 m/s 2 , according to the pressure calculation formula P = F/S = mg/S ≤ 0.206 kPa. The sensitivity is 16.019 mV/Pa, that is, when the sensor range is 3.3 V, the maximum pressure is 0.206 kPa, so the sensor sensitivity is 3300 mV/206Pa = 16.019 mV/Pa. The error is less than 5%. The sensor has the characteristics of strong anti-interference ability, high sensitivity, wide pressure range, good repeatability and strong consistency, which can be used for accurate detection of pulse signal waveform. The working principle of the new piezoelectric sensor is shown in Figure 2:

System Acquisition Process
Firstly, the system scanned the fixed position of the wrist through the new piezoelectric sensor array and the laser positioning module, collected three pulse signal data through three ADC at the sampling frequency of 125 Hz, obtained the waveform amplitude of the pulse signal and uploaded it through Bluetooth, combined with the adaptive pulse finding algorithm for data processing, positioned the three points with the largest average amplitude of the pulse waveform as the key points of the radial artery, and the new piezoelectric sensor array was deployed to the corresponding position by moving the wrist band. Secondly, the system applied a continuous decreasing dynamic pressure of 4

System Acquisition Process
Firstly, the system scanned the fixed position of the wrist through the new piezoelectric sensor array and the laser positioning module, collected three pulse signal data through three ADC at the sampling frequency of 125 Hz, obtained the waveform amplitude of the pulse signal and uploaded it through Bluetooth, combined with the adaptive pulse finding algorithm for data processing, positioned the three points with the largest average amplitude of the pulse waveform as the key points of the radial artery, and the new piezoelectric sensor array was deployed to the corresponding position by moving the wrist band. Secondly, the system applied a continuous decreasing dynamic pressure of 4 mmHg/s to the wristband in the range of 60-160 mmHg, three pulse signals and one wristband pressure values were collected by four-channel ADC at a sampling frequency of 125 Hz; the data was uploaded to the computer for noise reduction by Bluetooth, and the waveform amplitude of pulse signal at different key points of radial artery with dynamic pressure was obtained so as to obtain the optimal pulse pressure P and the waveform amplitude of the corresponding pulse signal y = u max 4096/3.3. Finally, the system pressured the wrist band to 60 mmHg, 80 mmHg, 100 mmHg through the air pump, three pulse signal data and one wristband pressure value were collected by four-channel ADC with a sampling frequency of 125 Hz similarly, the data were uploaded to the computer by Bluetooth for noise reduction, feature extraction and analysis. Thus, the pulse signal waveform under three static air pressure, the recognition results of the slippery pulse and the normal pulse, and the analysis of human health state were obtained. The system acquisition process is shown in Figure 3: were uploaded to the computer by Bluetooth for noise reduction, feature extraction and analysis. Thus, the pulse signal waveform under three static air pressure, the recognition results of the slippery pulse and the normal pulse, and the analysis of human health state were obtained. The system acquisition process is shown in Figure 3:

Adaptive Pulse Finding
The key points of the radial artery are different in different people. Direct positioning of the key points of the radial artery by mechanical devices leads to large positioning errors, which affects the subsequent acquisition of pulse signal waveform and is not conducive to the analysis of human health status. In this system, the protrusion of the radial stem bone is first located as the origin, and the coordinates are established with the junction between the bottom of the palm and the arm as the x-axis, the arm direction as the yaxis, and the waveform amplitude of the pulse signal as the z-axis. The laser module is placed on the outside of the wrist band, and its laser head and lens are close to the arm to locate the new piezoelectric sensor and obtain the displacement of the new piezoelectric sensor. Then, the new piezoelectric sensor array, combined with the laser positioning module, is used to scan the fixed position of the wrist from top to bottom and from left to right along the x axis and y axis directions, as shown in Formula (3), a ternary data set ( , , ) consisting of the position of the scanning point and the waveform amplitude of the pulse signal can be obtained. Then, the ternary data set is used for surface fitting and regression analysis to obtain Formula (4), and the relationship between pulse signal waveform amplitude and scanning point position is obtained. Then, conditional extreme value extraction is carried out (the distance between the three sensors is fixed at 1.5 cm), and the three points with the maximum average amplitude of pulse waveform ( 0 , 0 ), ( 1 , 1 ), ( 2 , 2 ) are obtained as the key points of the radial artery, as shown in Formula (5). Finally, by moving the wristband, the middle sensor is deployed to the position of the Guan pulse, and then the sensors at both ends are deployed to the position of the Cun pulse and

Adaptive Pulse Finding
The key points of the radial artery are different in different people. Direct positioning of the key points of the radial artery by mechanical devices leads to large positioning errors, which affects the subsequent acquisition of pulse signal waveform and is not conducive to the analysis of human health status. In this system, the protrusion of the radial stem bone is first located as the origin, and the coordinates are established with the junction between the bottom of the palm and the arm as the x-axis, the arm direction as the y-axis, and the waveform amplitude of the pulse signal as the z-axis. The laser module is placed on the outside of the wrist band, and its laser head and lens are close to the arm to locate the new piezoelectric sensor and obtain the displacement of the new piezoelectric sensor. Then, the new piezoelectric sensor array, combined with the laser positioning module, is used to scan the fixed position of the wrist from top to bottom and from left to right along the x axis and y axis directions, as shown in Formula (3), a ternary data set (x, y, z) consisting of the position of the scanning point and the waveform amplitude of the pulse signal can be obtained. Then, the ternary data set is used for surface fitting and regression analysis to obtain Formula (4), and the relationship between pulse signal waveform amplitude and scanning point position is obtained. Then, conditional extreme value extraction is carried out (the distance between the three sensors is fixed at 1.5 cm), and the three points with the maximum average amplitude of pulse waveform (x 0 , y 0 ), (x 1 , y 1 ), (x 2 , y 2 ) are obtained as the key points of the radial artery, as shown in Formula (5). Finally, by moving the wristband, the middle sensor is deployed to the position of the Guan pulse, and then the sensors at both ends are deployed to the position of the Cun pulse and the Chi pulse by adjusting the angle θ between the sensor array and the Y-axis, as shown in Formula (6), so where (x,y) is the coordinates of the scanning point, a[i] is the parameter, z is the waveform amplitude of the pulse signal, f (x i ,y i ) is the maximum average amplitude, (x i , y i ) are the coordinates of Cun, Guan and Chi, z i is the waveform amplitude of pulse signal at Cun, Guan and Chi, and θ is the angle between the sensor array and the Y-axis. The principle of adaptive pulse finding is shown in Figure 4:

Denoise Processing
The acquisition of pulse signals is usually affected by noise disturbances from hardware circuits and the baseline drift caused by self-respiration and body displacement. The short-time Fourier transform (STFT) of the signal can reflect the frequency characteristics of the signal well, the pulse signal with fixed frequency is selected through the window function to filter the noise, as shown in Formula (7), but its time resolution is fixed, and the details of the non-stationary signal can't be extracted. The mean filtering method can

Denoise Processing
The acquisition of pulse signals is usually affected by noise disturbances from hardware circuits and the baseline drift caused by self-respiration and body displacement. The short-time Fourier transform (STFT) of the signal can reflect the frequency characteristics of the signal well, the pulse signal with fixed frequency is selected through the window function to filter the noise, as shown in Formula (7), but its time resolution is fixed, and the details of the non-stationary signal can't be extracted. The mean filtering method can filter out high-frequency noise, but there are extreme values and the phase angle offset. The baseline can be removed by cubic spline interpolation of the signal, but the stationarity of the signal edge becomes worse. In this system, the collected pulse data are decomposed by wavelet to obtain the scaling function coefficient CA5[n] and the wavelet function coefficient CD5[n], CD4[n], CD3[n], CD2[n], CD1[n], as shown in Formula (8). Then the wavelet function coefficient is denoised by the soft and hard threshold compromise method, as shown in Formula (9). Then the scaling function coefficient and the wavelet function coefficient are reconstructed by wavelet to obtain the smooth pulse signal waveform. Finally, the iterative sliding window algorithm is used to remove the phase angle offset by setting the sliding window width as the pulse signal period and using the coordinate rotation formula as the window function recursion, as shown in Formula (10), remove the baseline, as shown in Formula (11), and the complete pulse signal waveform with smooth and steady edge is obtained. [n, x(n)] = [n, f (n)] * 1 −tanθ 0 1 , (−π/2 ≤ θ ≤ π/2) (10) In Formula (7) x ∑ n CDx[n]ψ x,n (t) is the detail component. In Formula (9), N is the threshold and λ is a constant. In Formula (10), θ is the angle between the connecting line and the horizontal line of the starting point and the ending point of the single period pulse, f (n) is the pulse signal processed by discrete wavelet transform (DWT). In Formula (11), T is the period of pulse signal, x(n) is the singly periodic pulse signal with the phase angle offset removed, and y(n) is the singly periodic pulse signal with the baseline drift removed. The wavelet transform and iterative sliding window processing of pulse signal are shown in Figure 5: mula (10), is the angle between the connecting line and the horizontal line of the starting point and the ending point of the single period pulse, ( ) is the pulse signal processed by discrete wavelet transform (DWT). In Formula (11), is the period of pulse signal, ( ) is the singly periodic pulse signal with the phase angle offset removed, and ( ) is the singly periodic pulse signal with the baseline drift removed. The wavelet transform and iterative sliding window processing of pulse signal are shown in Figure 5: (e) (f) After wavelet transform and iterative sliding window processing, the mixed waveform which contains external interference in the pulse signal and some waveform amplitude and period mutations is predicted and reconstructed according to the pulse waveform of the last period. The waveform reconstruction by curve fitting requires a large number of signal sampling points, and the algorithm complexity is high. The pulse waveform is a composite waveform formed by superposition of main wave, tide wave and dicrotic wave according to frequency and modal reconstruction. Firstly, the main wave is analyzed by regression analysis according to the peak and the period , as shown in Formula (12), combined with the corresponding curvature of the main wave, as shown in Formula (13), and the waveform of the main wave is obtained, as shown in Formula (14). After wavelet transform and iterative sliding window processing, the mixed waveform which contains external interference in the pulse signal and some waveform amplitude and period mutations is predicted and reconstructed according to the pulse waveform of the last period. The waveform reconstruction by curve fitting requires a large number of signal sampling points, and the algorithm complexity is high. The pulse waveform is a composite waveform formed by superposition of main wave, tide wave and dicrotic wave according to frequency and modal reconstruction. Firstly, the main wave is analyzed by regression analysis according to the peak Y m and the period T, as shown in Formula (12), combined with the corresponding curvature of the main wave, as shown in Formula (13), and the waveform of the main wave is obtained, as shown in Formula (14). Then, the waveform of the tidal wave can be obtained by regression analysis based on peak Y m1 and period (T 1 − b), as shown in Formula (15). Then, the waveform of the dicrotic wave can be obtained by regression analysis based on peak Y m2 and period (T − T 1 ), as shown in Formula (16). Then, the spatial and temporal domain relationships of the predicted main wave, the predicted tidal wave, and the predicted dicrotic wave are obtained in turn by the height of the peak and the trough of the main wave, tidal wave and dicrotic wave of the pulse signal of the last period h1, h2, h3, h4, h5, h6, and the occupied period t1, t2, t3, t4, t5, t6, as shown in Formula (17). Finally, the predicted main wave, the predicted tidal wave and the predicted dicrotic wave are superimposed in the period T and deployed to the corresponding position, as shown in Formula (18), and the predicted and reconstructed pulse signal waveform is obtained.
where T 1 is the time from starting point to the tidal wave valley, b is the time from starting point to the main wave valley, t s is the time occupied by the ascending branch, t s 1 is the time from starting point to the tidal wave peak, t s 2 is the time from starting point to the dicrotic wave peak, k 1 (t) is the corresponding curvature of the signal, f 1 (t) is the fundamental signal. c is a constant, y 1 (t) is the predicted main wave waveform, y 2 (t) is the predicted tidal wave waveform, y 3 (t) is the predicted dicrotic wave waveform, and Y(t) is the predicted reconstructed pulse signal waveform. The neural network is used to adjust b to change the width of the main wave, adjust t s to change the proportion of the rising branch, and adjust c to change the curvature of the waveform, so that the predicted reconstructed pulse waveform has high fitting degree, obvious characteristics, and meets the requirements of amplitude and period. The prediction and reconstruction of the pulse signal is shown in Figure 6: the dicrotic wave peak, 1 ( ) is the corresponding curvature of the signal, 1 ( ) is the fundamental signal. is a constant, 1 ( ) is the predicted main wave waveform, 2 ( ) is the predicted tidal wave waveform, 3 ( ) is the predicted dicrotic wave waveform, and ( ) is the predicted reconstructed pulse signal waveform. The neural network is used to adjust to change the width of the main wave, adjust to change the proportion of the rising branch, and adjust to change the curvature of the waveform, so that the predicted reconstructed pulse waveform has high fitting degree, obvious characteristics, and meets the requirements of amplitude and period. The prediction and reconstruction of the pulse signal is shown in Figure 6:

Feature Extraction and Classification
The difference in waveform between the sliding pulse and the normal pulse is the absence or presence of the dicrotic wave. Due to the influence of human factors, there is a condition in which the dicrotic wave are inconspicuous in the collected pulse signal waveform. The extreme value extraction method is used to determine whether the pulse signal waveform contains dicrotic wave, which leads to the inaccurate identification of the sliding pulse and the normal pulse, thus affecting the analysis of human health status. We commissioned Shanghai Bayes Health Technology Co., Ltd., Shanghai, China, to collect the pulse of 140 subjects from various TCM clinics through this system, including 100 healthy men and women aged 20-30 years, 30 men and women aged 20-30 years who were in the state of cold or fever, and 10 women aged 20-30 years who were in pregnancy. A total of 1260 sets of corresponding pulse signals were divided into the normal pulse and the slippery pulse labels. Then the pulse dataset was randomly shuffled and rearranged, of which 60% was used as the training set, 20% was used as the validation set, and 20% was used as the test set. Through the feature extraction and classification algorithm, this system first extracts the spatial-temporal features and the modal energy features of pulse signal, and then uses the machine learning method to train the classification model to identify the slip-pulse and normal pulse, so as to realize the accurate analysis of human health status. The principle of feature extraction and classification is shown in Figure 7:

Feature Extraction and Classification
The difference in waveform between the sliding pulse and the normal pulse is the absence or presence of the dicrotic wave. Due to the influence of human factors, there is a condition in which the dicrotic wave are inconspicuous in the collected pulse signal waveform. The extreme value extraction method is used to determine whether the pulse signal waveform contains dicrotic wave, which leads to the inaccurate identification of the sliding pulse and the normal pulse, thus affecting the analysis of human health status. We commissioned Shanghai Bayes Health Technology Co., Ltd., Shanghai, China, to collect the pulse of 140 subjects from various TCM clinics through this system, including 100 healthy men and women aged 20-30 years, 30 men and women aged 20-30 years who were in the state of cold or fever, and 10 women aged 20-30 years who were in pregnancy. A total of 1260 sets of corresponding pulse signals were divided into the normal pulse and the slippery pulse labels. Then the pulse dataset was randomly shuffled and rearranged, of which 60% was used as the training set, 20% was used as the validation set, and 20% was used as the test set. Through the feature extraction and classification algorithm, this system first extracts the spatial-temporal features F t and the modal energy features F y of pulse signal, and then uses the machine learning method to train the classification model to identify the slip-pulse and normal pulse, so as to realize the accurate analysis of human health status. The principle of feature extraction and classification is shown in Figure 7:   Firstly, the system processes the signal by wavelet transform, iterative sliding window and prediction reconstruction algorithm. Then, the ratio of ascending and descending branches t s /t f , the height ratio of main wave, tidal wave and dicrotic wave h 1 : h 2 : h 3 , and the occupied period ratio t 1 : t 2 : t 3 can be obtained by the period segmentation and the extreme value extraction, so as to obtain the spatial-temporal features F t t s /t f , t 1 : t 2 : t 3 , h 1 : h 2 : h 3 , as shown in Table 1: Then, by calculating the correlation coefficient r = cov(Y, Y i )/σσ i between the processed pulse signal Y and the components after wavelet decomposition Y i , the residual components Y 1 , Y 2 , Y 3 with r ≤ 0.2, and the detail components Y 4 , Y 5 and the scale component Y 6 with r > 0.2 are obtained. Extract wavelet components of r > 0.2 to calculate the corresponding mode energy y i = ||Y i ||/||Y||, so as to obtain the modal energy features F y {y 6 , y 5 , y 4 }, as shown in Table 2: Finally, in the case of a small data set and three independent repeated experiments, the classification models of the slippery pulse and the normal pulse are trained, respectively, by neural network (NN) and support vector machine (SVM), linear kernel function (Linear) and radial basis function (Rbf), combined with spatial time domain characteristic F t and mixed characteristic F F t , F y , and the recognition effects of different classification models are compared to select the best classification model, as shown in Table 3: Since the slippery pulse corresponds to the pregnant women or the human body is in the state of anemia, rheumatism, acute infection and fever, while the normal pulse corresponds to the human body is in the state of good nutrition and health [5][6][7], the system uses support vector machine (SVM), radial basis function (Rbf), combined with mixed features to train the slippery pulse and the normal pulse classification model, identify the slippery pulse and the normal pulse, realize the intelligent analysis of human health status, and its accuracy is up to 98.4%.

Discussion
In this paper, an intelligent wearable wrist pulse detection system based on a new piezoelectric sensor array is designed. The system adopts the main control circuit, new piezoelectric sensor array, wearable wristband hardware circuit design, combined with the algorithm of adaptive pulse finding, wavelet transformation, iterative sliding window, prediction and reconstruction, feature extraction and classification to locate the key points of the radial artery, reproduce the pulse signal waveform under the air pressure of different wristbands, eliminate the hardware circuit noise and external interference, extract the pulse signal features and identify the slippery pulse and the normal pulse, so as to realize the intelligent analysis of human health status. The sensitivity of the new piezoelectric sensor in the system is 16.019 mV/Pa, the pressure response range is 0.206 kPa, and the error is less than 5%. The pulse finding algorithm can realize millimeter-level location of radial artery key points, and the pulse signal waveform with smooth waveform, stable edge and obvious waveform characteristics can be obtained by noise reduction processing. The accuracy rate of the slippery pulse and the normal pulse identified by feature extraction and classification algorithm is 98.4%, which can accurately identify human health status.
The proposed system and other pulse measurement devices presented in recent studies are summarized in Table 4. Clemente F [22] proposed a piezo-film-based measurement system. Jessica [9] proposed a wearable pulse-taking device, and Sun Y [16] proposed a wearable pulse wave monitoring system. In conclusion, the wrist pulse acquisition system is feasible, and has the characteristics of small size, accurate pulse positioning, good noise reduction effect, and accurate intelligent analysis. It is suitable for different groups of people and meets the needs of family health care. In the future, we would integrate the accelerometer and gyroscope into the wearable system to effectively detect the arm's motion state and eliminate the inhibitory motion artifacts. At the same time, more experiments will be carried out on the basis of traditional Chinese medicine to perfect the pulse database, and the relationship between pulse condition and human health status will be established through machine learning, and the patients with different diseases and healthy personnel can be classified and identified precisely, which is of great significance for the objective research of traditional Chinese medicine pulse diagnosis.