Abstract
Introduction
The photoplethysmogram (PPG) is a low-cost and non-invasive method capable of analyzing blood pressure, heart rate, and other physiological information. It is widely used for diagnosing physical health and has important medical research value. 1 A PPG is produced mainly by irradiating human vascular tissue using light with a wavelength of 900 nm and collecting changes in light intensity through a receiver. Due to the continuous expansion and contraction of blood vessels in the human body, the internal blood volume changes periodically, which leads to a difference in the amount of light absorbed by the blood and a change in the intensity of received light. A PPG signal is formed from photoelectric conversion, which reflects the change in blood vessel volume and the health of the cardiovascular system.2,3 PPG signals express many biological information. 4 After the decomposition and reconstruction of PPG signals, a potential relationship can be studied between different PPG signal characteristics and cardiovascular health in humans. The effects of the sympathetic and parasympathetic nervous system on cardiovascular regulation can also be analyzed.5–8 Considering that the effect of moxibustion treatment is difficult to define directly using biochemical indicators, PPG signal analysis is a promising method for quantitatively defining the relationship between moxibustion treatment and human cardiovascular health.
Moxibustion is a type of traditional Chinese medicine (TCM) treatment that can benefit health by adjusting far-infrared (FIR) radiation. It is practiced by burning moxa cones on the acupoints of the body, thereby influencing physiological and biochemical functions. Moxibustion has far-reaching influence in Asia and around the world, and it represents a non-pharmaceutical intervention that has long-term applications in clinical practice. 7
Moxibustion can promote cardiovascular function and improve the regulation of the nervous system by stimulating the acupoints. Fang and Zhou 9 found that after electrical stimulation of acupoints such as st 36 (zusanli), sp 6 (sanyinjiao), and li 11 (quchi), the mean arterial pressure (MAP) of animals recovered rapidly, and myocardial ischemia and cardiac function recovery were improved, thereby protecting the myocardium. Liu 10 studied the literature on the treatment of cardiovascular diseases and found that the stimulation of st 36 (zusanli), sp 6 (sanyinjiao), li 11 (quchi), and other acupoints played an important role in the recovery of cardiovascular function. Peng et al. studied the regulation of endocrine function by vascular endothelial cells and found that the generation and clearance of oxygen free radicals when moxibustion treatment is used lowers the blood pressure of patients with hypertension. Therefore, moxibustion stimulation on acupoints such as st 36 (zusanli), sp 6 (sanyinjiao), and li 11 (quchi) can promote the recovery of brain cell function; improve the working ability of cerebral cortical cells; improve heart function; regulate heart rhythm; and increase red blood cell count, white blood cell count, blood sugar, blood oxygen saturation, vasodilation, blood rheology, and neuromodulation, all of which have significant effects. 11 Coyle et al. 12 applied moxibustion for patients with chronic obstructive pulmonary disease (COPD) and achieved positive therapeutic effects. Ni et al. 13 conducted studies on the effects of acupuncture and moxibustion on patients with thoracic arteritis, and research results showed that acupuncture combined with moxibustion could regulate humoral immunity and cellular immunity. In addition, moxibustion stimulation on acupoints such as st 36 (zusanli), sp 6 (sanyinjiao), and li 11 (quchi) can promote cerebral function; improve the working ability of the cerebral cortex cells; benefit heart function; adjust the rhythm of the heart; increase the counts of red blood cells and white blood cells; and bring about significant changes in the neuromodulation of blood sugar, blood oxygen saturation, blood vessels, diastolic blood pressure, 14 and blood rheology.
The above studies have expounded the relationship between the effect of moxibustion treatment and the change in cardiovascular function in the human body from the perspective of medical science. However, this research can only give a rough description using the change of pulse signal speed, and it is difficult to use real-time human physiological information data for characterization because these measurements lack objectivity and convenience, which limits the clinical application of moxibustion. A study confirmed that PPG signals reflect changes in blood pressure.
15
In order to further study the relationship between moxibustion treatment and human health, it is necessary to use PPG signals to analyze the influence of moxibustion on the amplitude and frequency of pulse signals. A research method for analyzing the frequency corresponding to the maximum energy (
Basic theory and methods
HHT and HT power spectrum analysis of PPG signals
Since the pulse wave signal is a non-linear and non-stationary physiological signal, it is difficult to extract the correct signal characteristics and analyze the real-time physiological signal through the traditional Fourier transform and wavelet transform analysis. 14
In recent years, with the development of HHT, the application of HHT for pulse waves has been analyzed.16,17 The HHT is a time–frequency signal analysis method proposed by Huang et al. in 1998, and it mainly studies the signal frequency characteristics through the signal’s intrinsic mode function (IMF), which is obtained by empirical mode decomposition (EMD). Since the IMF is able to provide a quantitative analysis of energy changes at different frequencies, its results are more accurate than the wavelet or Fourier transform, highly efficient and adaptive, and have important applications in the study of energy–frequency–time distributions.18,19 Therefore, the HHT can overcome the problems of non-adaptive functions and allow for non-linear and non-stationary signal analysis.
In the HHT process, EMD is first used to obtain the IMF and find the instantaneous component of the signal associated with the pulse wave transform. Then, the analytic function is constructed through the Hilbert transform (HT) to calculate the instantaneous frequency and amplitude of the signal and then analyze the information contained in the signal.
For example, Huang et al. 20 applied this method in their study on the HHT and Poincaré plot analysis (PPA) heartbeat segmentation algorithm, and Yimin et al. 21 studied the relationship between pulse signal and cardiovascular status. Haidong. 22 researched the frequency behaviors of the Hilbert spectrum and Hilbert marginal spectrum (HMS) to estimate the spectral traits of heart rate variability (HRV) signals.
The HHT is based on the instantaneous frequency obtained using EEMD (Ensemble Empirical Mode Decomposition), and the HT of the IMF component signal
It can be expressed in polar coordinates as
which includes
where
Therefore, the real part of the signal
where
Accordingly, the marginal spectrum
Definition of frequency extreme point
According to the Sadrawi et al.,
17
In Figure 1, the IMF component of the PPG signal obtained by EEMD before moxibustion (Figure 1(a)) and the marginal power spectrum of IMF components after HT (Figure 1(b)) are shown. The frequency corresponding to the maximum energy in the figure is the value of

The intrinsic mode function of PPG and its corresponding power spectrum
Experimental results
Experimental test
Research subjects
Six healthy young volunteers (five males and one female) were recruited in the study. People with known systemic diseases as well as those with smoking and drinking habits were excluded from the study. All participants were required not to drink alcohol or coffee for the 24 h before the test. In addition, all the volunteers’ basic personal information (including age, gender, weight, height, body mass index, etc.) was recorded. All of the volunteers had not received moxibustion treatment in the month prior to the test. Test participants were informed of moxibustion therapy, including experimental devices, stimulation methods, treatment intensity and frequency, and possible adverse events. Each participant signed a written consent form. The agreement and procedure of the study are in line with the principles of Helsinki Declaration and approved by the medical ethics review committee of Ningxia Medical University.
Equipment and data acquisition
In the study, the signal of decomposed short-time PPG was acquired using a self-developed device designed for human fingers. The volunteers’ index fingers were selected as the data acquisition object. In order to get better accuracy of data, the signal frequency was chosen at 500 Hz, and each data record lasted for 20 min. The collected data were uploaded to the computer through the analog-digital acquisition card-USB-6009 DAQ (National Instruments, Austin, TX, USA). LabVIEW 2014 and MATLAB 2015 are used for data acquisition and processing operations.
The acquired PPG signal was decomposed by EEMD and then used to calculate the IMF because pulse parameters could be measured accurately and reliably in the PPG signal. In order to ensure that the data were not affected by the environment, all PPG signals collected within 1 min before and after the collection are discarded.
26
The remaining data were decomposed by the HHT at 6-s intervals. The data involved in HHT decomposition were 3000 points. In total, 180 time decomposition processes were implemented, and 180 instantaneous frequencies were obtained. On this basis, 15 maximum values and 15 minimum values of all instantaneous frequencies were discarded, and the remaining 150 maximum instantaneous frequencies were averaged as the
Experimental test
The temperature in the room was kept at 26°C, and a quiet environment was ensured. The testers were in a relaxed supine position and not affected by the external environment during the entire test. The PPG infrared sensor was used for data collection using the index finger of the volunteers’ left hand. 26
Before data collection, the volunteers rested for 5 min. Then, the data were collected for 20 min before moxibustion, recording approximately 600,000 numbers of signals. During the moxibustion process, professional doctors were responsible for sticking the moxa cone to the li 11 (quchi) point, st 36 (zusanli) point, and sp 6 (sanyinjiao) point and keeping it burning at the point for 5 min, and then instructing the volunteers to rest for 10 min. After the end of moxibustion treatment, the data were collected again after the moxibustion, with an acquisition time of 20 min.
During processing, both the amount of data collected before and after moxibustion was 600,000 signals. The obtained data were then filtered. The filtering method was mainly adopted from the method of subsection comparison: firstly, all data were statistically counted to find the data cycle benchmark, and then the data benchmark was compared with each group of data to remove the abnormally large local data and avoid data pollution. The filtered data segment length was 6 s of PPG signals with a sampling point number of 3000. After all the data was decomposed by EEMD, IMFs and residuals could be generated, as shown in Figure 2(a). By Hilbert–Huang spectral analysis, the instantaneous frequency of different energies corresponding to each IMF could be obtained, as shown in Figure 2(b).

IMF and HT power spectrum decomposed from the standard PPG segment: (a) IMF intrinsic function of the volunteers’ pulse signals measured before moxibustion and (b) HT power spectrum of the volunteers’ pulse signals measured before moxibustion.
Results and analysis
Based on the above research method, the moxibustion effect tracking study was carried out with the help of six volunteers. The study lasted for 6 days and mainly involved recording
In order to track the continuous change of cardiovascular status during the moxibustion treatment cycle, a volunteer was selected and observed to analyze the changing trend of

Effects of moxibustion treatment on
In order to further study the effect of moxibustion on

The change of mean
Figure 4(b) shows the increasing trend of the
In order to further study the correlation between

Changes of the mean
Discussion
According to the theory of TCM, the points of li 11 (quchi), st 36 (zusanli), and sp 6 (sanyinjiao) are important acupoints of human body. In this study, all participating volunteers recorded their pulse signals before and after moxibustion, and HHT was used to obtain
The frequency information of PPG signal can reflect the state of blood vessels—distinguish between healthy people and diabetics.
26
In the field of TCM treatment, relevant scholars have studied that the heat-sensitive moxibustion had neuroprotective effects against focal cerebral ischemia/reperfusion injury.
27
In this study,
During the 6-day experiment, the changes in PPG signals of the volunteers reflected the effects of moxibustion on the human body, both individually and as a group. The value of
In addition, the amplitude of
Conclusion
Through using the HHT, the effect of moxibustion treatment on
After moxibustion, the amplitude of human PPG signal (except IMF9) decreases, and the energy is concentrated from the high-frequency part of the human pulse signal to the low-frequency part, which is consistent with the idea that moxibustion can reduce the effect of organs on blood supply according to the theory of TCM. In the follow-up study, the plan for further research is to focus on the efficacy of moxibustion, improve the comparative test, study the comprehensive effect of different points of moxibustion on human and pulse signals, and get further confirmation on the actual efficacy of moxibustion.
