Pulse wave data processing method and related equipment

Abstract

The embodiment of the present invention provides a pulse wave data processing method, including: obtaining a pulse wave signal of a designated target in a target drug state; extracting harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal; calculating the variables of the harmonic frequency characteristic data to obtain variable distribution data of the pulse wave signal. After obtaining the pulse wave signal of the designated target in a target drug state, the variable distribution of the pulse wave frequency domain characteristics in the target drug state on the arterial system can be quantified, aiming to solve the problem that the pulse wave frequency domain characteristics provided by the prior art cannot be quantitatively analyzed because the complete characteristics of a Chinese herbal medicine include many complex volatile components, hydrophilic components, insoluble components, etc., resulting in different medicinal properties of some complex components obtained by short cooking and long cooking, high fire and low fire, processing, essential oil purification, etc., thereby converting the pulse waves of Chinese herbal medicines in different drug states into measurable, reproducible, and intuitive quantitative analysis data.

Description

Technical Field

[0001] The present invention relates to the field of signal processing, and in particular to a pulse wave data processing method and related equipment. Background Art

[0002] Pulse wave frequency domain feature analysis is also called pulse wave harmonic feature analysis. It is an important part of hemodynamic research and analysis. The amplitude and phase of the obtained pulse wave harmonics are important physiological parameters in the human blood circulation system, reflecting the state of the arterial system and conveying important information about the pathological changes of cardiovascular diseases. Therefore, the amplitude and phase characteristics of pulse wave harmonics are also used as a quantitative test method for the impact of traditional Chinese medicine on the circulatory system. They are applied in pharmacokinetics and are called pulse pharmacodynamics research. Therefore, it is very important to accurately obtain stable pulse wave frequency domain features. In the classic Chinese medicine book "Huangdi Neijing", the "Suwen Sanbu Jiuhou Lun Pian" is a summary of the clinical pulse analysis by Chinese medical predecessors, which shows that the change of pulse wave waveform characteristics between different arteries is an important pulse representation of the human body state. However, the complete characteristics of a Chinese herbal medicine include many complex volatile components, hydrophilic components, insoluble components, etc. Therefore, due to the different medicinal properties of some complex components obtained by short and long cooking, high and low fire, processing, essential oil purification, etc., the characteristics of Chinese herbal medicine cannot be quantitatively analyzed. Summary of the invention

[0003] The embodiment of the present invention aims to provide a pulse wave data processing method, which can quantify the variable distribution of the pulse wave frequency domain characteristics in the target drug state on the arterial system after obtaining the pulse wave signal of a specified target in the target drug state, and aims to solve the problem that the pulse wave frequency domain characteristics provided by the prior art cannot be quantitatively analyzed because the complete characteristics of a Chinese herbal medicine include many complex volatile components, hydrophilic components, insoluble components, etc., resulting in short cooking and long cooking, high fire and low fire, processing, essential oil purification, etc., which results in different medicinal properties of some complex components. Therefore, the pulse waves of the Chinese herbal medicine in different drug states are converted into measurable, reproducible and intuitive quantitative analysis data.

[0004] In a first aspect, an embodiment of the present invention provides a pulse wave data processing method, the pulse wave data processing method comprising:

[0005] Acquire a pulse wave signal of a designated target in a target drug state;

[0006] Extracting harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal;

[0007] The variables of the harmonic frequency characteristic data are calculated to obtain the variable distribution data of the pulse wave signal.

[0008] Optionally, the target drug state includes three drug states corresponding to pure dew, decoction, and solid matter, and obtaining the pulse wave signal of the specified target in the target drug state includes:

[0009] The pulse wave signals corresponding to the designated target in the pure dew drug state, the decoction drug state, and the solid drug state are obtained.

[0010] Optionally, performing harmonic frequency feature extraction on the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal includes:

[0011] Dividing the pulse wave signal according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence including segmentation signals having the same number of periods as the pulse wave signal;

[0012] Compound series conversion is performed on each segmentation signal in the segmentation sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0013] Optionally, performing compound series conversion processing on each of the segmented signals in the segmented sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal includes:

[0014] Performing compound series conversion processing on each of the segmentation signals in the segmentation sequence to obtain harmonic frequency data of the segmentation signal;

[0015] The harmonic frequency data of the segmented signal is normalized to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0016] Optionally, performing compound series conversion processing on each of the segmented signals in the segmented sequence to obtain harmonic frequency data of the segmented signals includes:

[0017] Calculating the average amplitude and average phase of the pulse wave signal;

[0018] The harmonic frequency data of the segmented signal is normalized according to the average amplitude and average phase of the pulse wave signal to obtain the harmonic frequency characteristic data corresponding to the pulse wave signal.

[0019] Optionally, calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal includes:

[0020] The percentage variable of the harmonic frequency characteristic data is calculated to obtain the percentage variable distribution data of the pulse wave signal.

[0021] Optionally, calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal includes:

[0022] The phase difference variable of the harmonic frequency characteristic data is calculated to obtain the phase difference variable distribution data of the pulse wave signal.

[0023] In a second aspect, an embodiment of the present invention provides a pulse wave data processing device, the device comprising:

[0024] An acquisition module, used to acquire a pulse wave signal of a designated target in a target drug state;

[0025] An extraction module, used for extracting harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal;

[0026] The calculation module is used to calculate the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal.

[0027] In a third aspect, an embodiment of the present invention provides an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the pulse wave data processing method provided in the embodiment of the present invention when executing the computer program.

[0028] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps in the pulse wave data processing method provided in the embodiment of the invention are implemented.

[0029] In an embodiment of the present invention, a pulse wave signal of a designated target in a target drug state is obtained; harmonic frequency feature extraction is performed on the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal; and variables of the harmonic frequency feature data are calculated to obtain variable distribution data of the pulse wave signal. After obtaining the pulse wave signal of a designated target in a target drug state, the variable distribution of the pulse wave frequency domain feature in the target drug state on the arterial system can be quantified, aiming to solve the problem that the pulse wave frequency domain feature provided by the prior art cannot be quantitatively analyzed because the complete characteristics of a Chinese herbal medicine include many complex volatile components, hydrophilic components, insoluble components, etc., resulting in different medicinal properties of some complex components obtained by short cooking and long cooking, high fire and low fire, processing, essential oil purification, etc., thereby converting the pulse wave of the Chinese herbal medicine in different drug states into measurable, reproducible, and intuitive quantitative analysis data. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0031] Figure 1 is an architecture diagram of a pulse wave data processing system provided by an embodiment of the present invention;

[0032] Figure 2 It is a structural schematic diagram of a medicinal material distillation device provided by an embodiment of the present invention;

[0033] Figure 3 It is a flow chart of a pulse wave data processing method provided by an embodiment of the present invention;

[0034] Figure 4 This is a pulse wave diagram before and after taking Chinese herbal medicine provided by an embodiment of the present invention;

[0035] Figure 5 A schematic diagram of the change of the average value of harmonic frequency amplitude is provided in an embodiment of the present invention;

[0036] Figure 6 A schematic diagram of the variation of harmonic frequency amplitude variation coefficient provided by an embodiment of the present invention;

[0037] Figure 7 A schematic diagram of harmonic phase change is provided in an embodiment of the present invention;

[0038] Figure 8 A schematic diagram of a standard phase deviation change is provided in an embodiment of the present invention;

[0039] Fig. 9 It is a structural schematic diagram of a pulse wave data processing device provided by an embodiment of the present invention;

[0040] Fig.10 It is a structural schematic diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0041] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0042] See also Figure 1 , Figure 1is an architecture diagram of a pulse wave data processing system provided by an embodiment of the present invention, such as Figure 1 As shown, the pulse wave data processing system includes: a medicinal material distillation device, a pressure signal acquisition module, a pulse harmonic frequency analysis unit and a harmonic frequency characteristic variation analysis unit.

[0043] Among them, the medicinal material distillation device includes a constant temperature timed heating boiler, a fan, a condenser, and a pure dew collection bottle; the pressure signal acquisition module includes a pressure sensor and a DC and AC pulse signal acquisition circuit; the pressure signal acquisition module is connected to the pulse harmonic frequency analysis unit signal, and the pressure signal acquisition module is used to transmit the collected pulse signal to the pulse wave harmonic frequency analysis unit; the pulse wave harmonic frequency analysis unit is connected to the harmonic frequency characteristic variation analysis unit signal.

[0044] For details, please refer to Figure 2 , Figure 2 Schematic diagram of a medicinal material distillation device provided by an embodiment of the present invention. Figure 2 As shown, the condenser is arranged above the constant temperature timing heating boiler, the fan is arranged above the condenser, the outlet of the condenser is provided with a liquid collecting nozzle, and a pure dew collecting bottle is arranged below the liquid collecting nozzle. More specifically, the Chinese herbal medicine distillation can be carried out through the following steps:

[0045] D1. Prepare raw materials: wash, slice, chop and flatten the target Chinese herbal medicines to facilitate cooking and purification. The target Chinese herbal medicines are fresh Chinese herbal medicines.

[0046] D2. Put the target Chinese herbal medicine into a constant temperature and timed heating boiler, fully soak and cover the target Chinese herbal medicine with water or liquid solvent, heat the constant temperature and timed heating boiler to a specific constant temperature for distillation, and collect the condensed liquid hydrosol through the fan, condenser and hydrosol collection bottle above the constant temperature and timed heating boiler. The above hydrosol may contain essential oils, water, and low-boiling volatiles. After distillation for a fixed time, open the lid to cool, and then separate the distilled soup and solid matter in the boiler by filtering.

[0047] In a possible embodiment, the temperature of the constant temperature heating in the above step D2 can be maintained at 80-120 degrees Celsius, the distillation time in the above step D2 can be 30-420 minutes, the liquid solvent in the above step D2 can be a hydrophilic solvent, the volume of the liquid solvent in the above step D2 can be 1 to 10 times that of the medicinal material, and the low-pressure and low-temperature distillation method can be used in the above step D2, wherein the reduced-pressure distillation temperature is 35-80 degrees Celsius.

[0048] In an embodiment of the present invention, a pulse wave signal of a designated target in a target drug state is obtained, and the designated target may be a mammal, such as a white mouse, a white rabbit, or other mammal, wherein the target drug state includes three drug states, namely, before and after taking a pure dew, before and after taking a decoction, and before and after taking a solid. By distilling Chinese herbal medicines, pure dew, decoction, and solids of Chinese herbal medicines are obtained, and pulse wave signals of mammals before taking, after taking the pure dew, decoction, and solids are recorded, and the collected pulse wave signals are analyzed as data samples.

[0049] In an embodiment of the present invention, three groups of mammals can be selected. The first group of mammals records the pulse wave signals of the mammals before and after taking the pure dew, the second group of mammals records the pulse wave signals of the mammals before and after taking the soup, and the third group of mammals records the pulse wave signals of the mammals before and after taking the solid matter.

[0050] Specifically, the pulse wave signal of a mammal before taking food, or after taking pure water, soup, or solid matter can be collected through the pressure signal collection module, the collected pulse wave signal can be harmonically processed through the pulse harmonic frequency analysis unit, and the pulse wave signal after harmonically processed can be subjected to variable analysis through the harmonic frequency characteristic variation analysis unit.

[0051] The above-mentioned signal connection can be a wired connection or a wireless connection. The above-mentioned wireless connection method can include but is not limited to 3G / 4G / 5G connection, WiFi connection, Bluetooth connection, WiMAX connection, Zigbee connection, UWB (ultrawideband) connection, and other wireless connection methods currently known or to be developed in the future.

[0052] See also Figure 3 , Figure 3 is a flow chart of a pulse wave data processing method provided by an embodiment of the present invention, such as Figure 3 As shown, the pulse wave data processing method specifically includes the following steps:

[0053] 301. Obtain a pulse wave signal of a designated target in a target drug state.

[0054] In an embodiment of the present invention, a pulse wave signal of a designated target in a target drug state is obtained, and the designated target may be a mammal, such as a white mouse, a white rabbit, or other mammal, wherein the target drug state includes three drug states, namely, before and after taking a pure dew, before and after taking a decoction, and before and after taking a solid. By distilling Chinese herbal medicines, pure dew, decoction, and solids of Chinese herbal medicines are obtained, and pulse wave signals of mammals before taking, after taking the pure dew, decoction, and solids are recorded, and the collected pulse wave signals are analyzed as data samples.

[0055] In an embodiment of the present invention, three groups of mammals can be selected. The first group of mammals records the pulse wave signals of the mammals before and after taking the pure dew, the second group of mammals records the pulse wave signals of the mammals before and after taking the soup, and the third group of mammals records the pulse wave signals of the mammals before and after taking the solid matter.

[0056] The pulse wave signal may be a pulse wave signal obtained by collecting the AC pressure signal from the pressure signal collection module, collecting the signal through an ADC, and removing high-frequency noise through digital filtering.

[0057] Optionally, in the step of acquiring the pulse wave signal of the specified target in the target drug state, the pulse wave signals of multiple specified human body parts can be collected under the condition that the difference in the average static pressure readings of any two pulse wave signals is kept less than a preset value to obtain the pulse wave signals of multiple specified human body parts.

[0058] Specifically, during the pulse wave signal acquisition process, the pressure control module can be used to maintain a constant pressure reading of the cuff airbag in the pressure signal acquisition module for 10-60 seconds, with the pressure value ranging from 60 to 100 mmHg. The difference in the average static pressure readings of any two cuff airbags is <3 mmHg. When extracting stable features from the frequency domain of a series of collected single-cycle signals, when the pulse rate of any single pulse wave is greater than 110% of the average pulse rate or the pulse rate of a single pulse wave is less than 90% of the average pulse rate, the single pulse wave will not be included in the subsequent calculation.

[0059] Optionally, the target drug state includes three drug states corresponding to pure dew, decoction, and solid matter. Specifically, in the step of obtaining the pulse wave signal of the designated target under the target drug state, the target drug state may include three drug states before and after taking pure dew, before and after taking decoction, and before and after taking solid matter. The pulse wave signal includes six types of pulse wave signals before and after taking pure dew, before and after taking decoction, and before and after taking solid matter.

[0060] 302. Extract harmonic frequency features of the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal.

[0061] In the embodiment of the present invention, it should be noted that the pulse wave frequency domain feature is also called the pulse wave harmonic feature. The initial pulse wave signal belongs to the time domain signal. The harmonic feature of the pulse wave signal can be extracted to obtain the harmonic feature data belonging to the frequency domain signal.

[0062] Optionally, the pulse wave signal can be segmented according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence includes segmentation signals with the same number of pulse wave signal periods; each segmentation signal in the segmentation sequence is subjected to compound series conversion processing to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0063] Specifically, the continuous pulse wave can be divided into individual pulse waves according to the lowest point, and the pulse wave of each part can be converted sequentially by compound series to obtain the harmonic frequency characteristic data Cn.

[0064] Optionally, in the step of performing compound series conversion processing on each segmented signal in the segmented sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal, each segmented signal in the segmented sequence can also be subjected to compound series conversion processing to obtain the harmonic frequency data of the segmented signal; and the harmonic frequency data of the segmented signal can be normalized to obtain the harmonic frequency characteristic data corresponding to the pulse wave signal.

[0065] Specifically, the continuous pulse wave can be divided into individual pulse waves according to the lowest point, and the harmonic frequency characteristic data Cn can be calculated sequentially for the pulse wave of each part by means of compound series conversion and normalization.

[0066] Optionally, in the step of performing compound series conversion processing on each segmentation signal in the segmentation sequence to obtain harmonic frequency data of the segmentation signal, the average amplitude and average phase of the pulse wave signal can be calculated; based on the average amplitude and average phase of the pulse wave signal, the harmonic frequency data of the segmentation signal is normalized to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0067] In an embodiment of the present invention, the harmonic frequency data includes amplitude data and phase data. The amplitude data of the segmented signal can be normalized to obtain harmonic frequency characteristic data of the pulse wave signal corresponding to the amplitude data. Specifically, the average amplitude of the pulse wave signal can be calculated; according to the average amplitude of the pulse wave signal, the amplitude data of the segmented signal can be normalized to obtain harmonic frequency characteristic data of the pulse wave signal corresponding to the amplitude data.

[0068] Furthermore, the harmonic frequency characteristic data corresponding to the amplitude data also includes a normalized amplitude standard deviation CnSD and a normalized harmonic frequency amplitude variation coefficient CnCV.

[0069] The phase data of the segmented signal can be normalized to obtain the harmonic frequency characteristic data of the pulse wave signal corresponding to the phase data. Specifically, the average phase of the pulse wave signal can be calculated; according to the average phase of the pulse wave signal, the phase data of the segmented signal can be normalized to obtain the harmonic frequency characteristic data of the pulse wave signal corresponding to the phase data, that is, the harmonic frequency phase average value.

[0070] Furthermore, the harmonic frequency characteristic data corresponding to the phase data also includes a standard deviation PnSD of the normalized phase.

[0071] For further explanation, the embodiment of the present invention takes a continuous pulse wave signal of more than 12 cycles as an example, and the continuous pulse wave signal is divided into cycles with the lowest point of each pulse wave as the dividing point; each pulse wave waveform contains the amplitude An and phase θn characteristics of each harmonic frequency, and the amplitude An and phase θn of each pulse wave can be obtained through the characteristic calculation of the compound series; wherein, the general expression of each pulse wave waveform is shown in the following formula 1:

[0072]

[0073] Where Y(K) is the discrete digitized value of a single pulse wave signal, An is the amplitude of the nth compound series, θn is the phase of the nth compound series, and n represents the nth harmonic frequency. N represents the total number of data points. f0 is the main frequency (fundamental frequency) of the pulse wave signal. The calculation expression of the above compound series can be shown as the following formula 2:

[0074]

[0075] The amplitude An is divided by the average value A0 of the pulse wave waveform to obtain the normalized harmonic frequency characteristic data Cn (also called harmonic frequency amplitude). After obtaining the amplitude An of each pulse wave, the normalized harmonic frequency amplitude average value CnAVG of a continuous pulse wave can be calculated by the following formulas 3 and 4. Formulas 3 and 4 are as follows:

[0076]

[0077]

[0078] Wherein, M is the total number of pulse waves in a continuous period of time, An,q is the amplitude of the nth compound series after the qth pulse wave is calculated, A0,q is the average value of the qth pulse wave, and CnAVG is the normalized average value of the harmonic amplitude of the nth harmonic frequency in a continuous period of time.

[0079] After obtaining the normalized harmonic frequency amplitude average value CnAVG of a continuous pulse wave, the normalized standard amplitude deviation CnSD of the nth harmonic frequency in a continuous period of time and the normalized harmonic frequency amplitude variation coefficient CnCV of the nth harmonic frequency in a continuous period of time can be calculated by the following formulas 5 and 6. Formulas 5 and 6 are as follows:

[0080]

[0081]

[0082] Where CnCV is the coefficient of variation of the normalized amplitude of the nth harmonic frequency in a continuous period of time.

[0083] After obtaining the phase θn of each pulse wave, the normalized harmonic phase average value Pn and the normalized standard phase deviation PnSD of a continuous pulse wave can be calculated by the following equations 7 and 8. The equations 7 and 8 are as follows:

[0084]

[0085]

[0086] Where M is the total number of pulse waves in a continuous period of time, θn,q is the phase of the nth compound series after the qth pulse wave is calculated, and n represents the nth harmonic frequency. Pn is the average phase value of the nth harmonic frequency in a continuous period of time. PnSD is the normalized standard phase deviation of the nth harmonic frequency in a continuous period of time.

[0087] 303. Calculate the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal.

[0088] In an embodiment of the present invention, the variables of the harmonic frequency characteristic data may include percentage variables of the harmonic frequency characteristic data and phase difference variables of the harmonic frequency characteristic data. The variables of the harmonic frequency characteristic data may be variables before and after taking Chinese herbal medicine.

[0089] Through the variables of the harmonic frequency characteristic data, the variable distribution data of the pulse wave signal before and after taking the Chinese herbal medicine can be obtained, and then the influence of the Chinese herbal medicine on the pulse wave signal in different drug states can be obtained. Specifically, through the percentage variables of the harmonic frequency characteristic data and the phase difference variables of the harmonic frequency characteristic data, the influence of the Chinese herbal medicine on the pulse wave signal in different drug states can be obtained.

[0090] Optionally, in the step of calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal, the percentage variables of the harmonic frequency characteristic data may be calculated to obtain the percentage variable distribution data of the pulse wave signal.

[0091] Specifically, the percentage variables of the harmonic frequency characteristic data include the harmonic frequency amplitude average value CnAVG percentage variable under the target drug state, the standard phase deviation PnSD percentage variable, and the harmonic frequency amplitude variation coefficient CnCV percentage variable.

[0092] More specifically, the percentage variables of the harmonic frequency characteristic data include the harmonic frequency amplitude average value CnAVG percentage variable, standard phase deviation PnSD percentage variable and harmonic frequency amplitude variation coefficient CnCV percentage variable before and after taking the pure dew; the harmonic frequency amplitude average value CnAVG percentage variable, standard phase deviation PnSD percentage variable and harmonic frequency amplitude variation coefficient CnCV percentage variable before and after taking the decoction; the harmonic frequency amplitude average value CnAVG percentage variable, standard phase deviation PnSD percentage variable and harmonic frequency amplitude variation coefficient CnCV percentage variable before and after taking the solid substance.

[0093] Specifically, the percentage variable of the above harmonic frequency characteristic data can be calculated by the following equations 9 to 11:

[0094]

[0095]

[0096]

[0097] CnAVG 0 Represents the average value of the harmonic frequency amplitude of the pulse wave signal before taking Chinese herbal medicine, CnCV 0 Table 2. PnSD coefficient of variation of harmonic frequency amplitude of pulse wave signal before taking Chinese herbal medicine 0 Represents the standard phase deviation of the pulse wave signal before taking Chinese herbal medicine; CnAVG 1 Represents the average value of harmonic frequency amplitude of pulse wave signal after taking Chinese herbal medicine, CnCV 1 Represents the harmonic frequency amplitude variation coefficient, PnSD, of the pulse wave signal after taking Chinese herbal medicine 1 Represents the standard phase deviation of the pulse wave signal after taking Chinese herbal medicine. △CnAVG (%) represents the harmonic frequency amplitude average value CnAVG percentage variable under the target drug state, △CnCV (%) represents the harmonic frequency amplitude variation coefficient CnCV percentage variable under the target drug state, and △PnSD (%) represents the standard phase deviation PnSD percentage variable under the target drug state.

[0098] Optionally, in the step of calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal, the phase difference variables of the harmonic frequency characteristic data may be calculated to obtain the phase difference variable distribution data of the pulse wave signal.

[0099] Specifically, the percentage variable of the harmonic frequency characteristic data includes the harmonic frequency phase average value Pn variable under the target drug state. The above phase difference variable distribution data is the harmonic frequency phase average value Pn variable.

[0100] More specifically, the percentage variables of the harmonic frequency characteristic data include the harmonic frequency phase average value Pn variable before and after taking the pure dew; the harmonic frequency phase average value Pn variable before and after taking the decoction; and the harmonic frequency phase average value Pn variable before and after taking the solid substance.

[0101] The above harmonic phase average value Pn variable can be calculated by the following formula 12:

[0102] ΔPn(%)=Pn 1 -Pn 0 (12)

[0103] The above Pn 0 Represents the harmonic phase average of the pulse wave signal before taking Chinese herbal medicine, Pn 1 Represents the harmonic frequency phase average value of the pulse wave signal after taking Chinese herbal medicine. △Pn (%) represents the phase deviation value of the harmonic frequency phase average value under the target drug state. The phase deviation value can be used to present the changes in the harmonic frequency characteristics of the pulse before and after taking Chinese herbal medicine.

[0104] In an embodiment of the present invention, a pulse wave signal of a designated target in a target drug state is obtained; harmonic frequency feature extraction is performed on the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal; and variables of the harmonic frequency feature data are calculated to obtain variable distribution data of the pulse wave signal. After obtaining the pulse wave signal of a designated target in a target drug state, the variable distribution of the pulse wave frequency domain feature in the target drug state on the arterial system can be quantified, aiming to solve the problem that the pulse wave frequency domain feature provided by the prior art cannot be quantitatively analyzed because the complete characteristics of a Chinese herbal medicine include many complex volatile components, hydrophilic components, insoluble components, etc., resulting in different medicinal properties of some complex components obtained by short cooking and long cooking, high fire and low fire, processing, essential oil purification, etc., thereby converting the pulse wave of the Chinese herbal medicine in different drug states into measurable, reproducible, and intuitive quantitative analysis data.

[0105] The embodiment of the present invention can reflect the changes in the state of the circulatory system caused by taking traditional Chinese medicine through the variable distribution data of the pulse wave signal, and can be applied to the research related to the identification and analysis of hemodynamics and the pharmacodynamics of Chinese herbal medicines, or as the cornerstone of the quantitative data of the three parts and nine signs of traditional Chinese medicine. It can see more detailed and multi-dimensional change characteristics in the pulse or hemodynamic state, which has a great positive promoting effect on the scientific quantitative research of traditional Chinese medicine pulse diagnosis. It can intuitively view the medicinal effect data of Chinese herbal medicines, which has more iconic features in the identification and analysis of hemodynamic state and pharmacodynamics of Chinese herbal medicines, and can fully grasp the influence of a single Chinese herbal medicine on the pulse, as well as the main physicochemical properties of its source components. At the same time, it can also observe the metabolic characteristics of these separated composite components over time, which is a very important cornerstone of data quantification for the study of the pharmacodynamics of Chinese herbal medicines.

[0106] Please refer to Figure 4 , Figure 4 1 is a pulse wave diagram before and after taking Chinese herbal medicine provided by an embodiment of the present invention, wherein part (A) is the pulse wave before taking the Chinese herbal medicine, and part (B) is the pulse wave after taking the Chinese herbal medicine. It can be seen that there is a difference in waveform between the pulse wave before taking the Chinese herbal medicine and the pulse wave after taking the Chinese herbal medicine. By implementing the embodiment of the present invention, the pulse waves of Chinese herbal medicine in different drug states can be converted into measurable, reproducible and intuitive quantitative analysis data.

[0107] Specifically, taking the quantitative full spectrum analysis of pulse and medicinal properties of 15-year-old Changbai Mountain ginseng as an example, ginseng pure water, ginseng decoction and ginseng solid are separated by the Chinese herbal medicine distillation method provided by the embodiment of the present invention, and the pulse wave signals of healthy people or healthy mammals are collected before taking the Chinese herbal medicine. 30 minutes after the healthy people or healthy mammals take the ginseng pure water, ginseng decoction and ginseng solid, respectively, the pulse wave signals after taking the Chinese herbal medicine are collected.

[0108] First, the harmonic frequency amplitude and change can be obtained according to the pulse wave data processing method provided in the embodiment of the present invention, such as Figure 5 As shown, Figure 5 Schematic diagram of the change of the average value of the harmonic amplitude provided by an embodiment of the present invention. Figure 5 As shown in the figure, it can be seen that the changes in the amplitude of the harmonic frequency of the pulse wave after the subjects took the "pure dew, soup, and solid matter" of 15-year-old white mountain ginseng for 30 minutes. n represents the nth harmonic frequency of the compound leaf series. Specifically, it can be seen that drinking the soup separated from 15-year-old white mountain ginseng has the greatest impact on the amplitude of the harmonic frequency, among which the amplitude of the third harmonic frequency increases the most, so it can be concluded that "ginseng soup" has a greater impact on blood dynamics. Relatively speaking, "pure dew and solid matter" have less impact on the characteristics of the harmonic frequency amplitude.

[0109] Secondly, the change in the average value of the harmonic frequency amplitude can be obtained according to the pulse wave data processing method provided in the embodiment of the present invention, such as Figure 6 As shown, Figure 6 Schematic diagram of the variation of harmonic frequency amplitude variation coefficient provided by an embodiment of the present invention. Figure 6 As shown in the figure, it can be seen that the change in the pulse wave harmonic amplitude variation coefficient after the subjects took the 15-year-old Changbai Mountain ginseng "pure dew, soup, solid" for 30 minutes. n represents the nth harmonic frequency of the compound leaf series. Specifically, it can be seen that the change in the pulse wave harmonic amplitude variation coefficient CnCV is still most directly affected by the Changbai Mountain ginseng "soup", while the "pure dew, solid" has a smaller impact.

[0110] Then, the harmonic frequency phase change can be obtained according to the pulse wave data processing method provided by the embodiment of the present invention. Figure 7 As shown, Figure 7 Schematic diagram of harmonic phase change provided by an embodiment of the present invention. Figure 7 As shown in the figure, it can be seen that the change in the pulse wave harmonic amplitude variation coefficient CnCV after the subjects took the 15-year-old Changbai Mountain ginseng "pure dew, soup, solid" for 30 minutes. n represents the nth harmonic frequency of the compound leaf series. Specifically, it can be seen that the change in the pulse wave harmonic phase average value Pn is still in the same direction as the influence of Changbai Mountain ginseng "decoction and pure dew". The change intensity of the influence of "decoction" is higher, while the direction of "solid" affecting Pn is opposite to the former two, showing an antagonistic effect, and has a greater impact on the 9-11 high-frequency harmonics.

[0111] Finally, the standard phase deviation can be obtained according to the pulse wave data processing method provided by the embodiment of the present invention. Figure 8 As shown, Figure 8 Schematic diagram of a standard phase deviation change provided by an embodiment of the present invention. Figure 8 As shown in the figure, it can be seen that the change in the standard phase deviation of the pulse wave after the subjects took the "pure dew, soup, and solid" of 15-year-old Changbai Mountain ginseng for 30 minutes. n represents the nth harmonic frequency of the compound leaf series. Specifically, it can be seen that the "decoction" of Changbai Mountain ginseng still has the greatest impact on the change of the harmonic phase standard deviation PnSD of the pulse wave, and the increase in the 3rd and 4th harmonic frequencies exceeds 300%. It means that the "ginseng soup" has a greater impact on the phase standard deviation (blood) of the spleen meridian and lung meridian, while the effects of ginseng pure dew and solid are far less than those of the soup.

[0112] In summary, the embodiment of the present invention can obtain the complete quantization process of the harmonic amplitude average value CnAVG (the harmonic amplitude average value can also be called the harmonic amplitude), the harmonic amplitude variation coefficient CnCV, the harmonic phase average value Pn (the harmonic phase average value can also be called the harmonic phase), and the phase standard deviation PnSD. The quantified data can be used to more comprehensively understand each harmonic characteristic aspect of the volatile components, water-soluble components, and solid components of the 15-year Changbai Mountain ginseng that affect the circulatory system, which explains that the main influencing component of ginseng on the change of pulse is the water-soluble decoction, followed by the volatile pure dew, while the solid matter has less impact on the circulatory system.

[0113] It should be noted that the pulse wave data processing method provided in the embodiment of the present invention can be applied to devices such as smart phones, computers, servers, etc. that can process pulse wave data.

[0114] Optional, see Fig. 9 , Fig. 9 is a schematic diagram of the structure of a pulse wave data processing device provided by an embodiment of the present invention, such as Fig. 9 As shown, the device comprises:

[0115] An acquisition module 901 is used to acquire a pulse wave signal of a designated target in a target drug state;

[0116] An extraction module 902 is used to extract harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal;

[0117] The calculation module 903 is used to calculate the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal.

[0118] Optionally, the target drug state includes three drug states corresponding to hydrosol, decoction, and solid matter, and the acquisition module 901 includes:

[0119] The acquisition submodule is used to obtain the pulse wave signals corresponding to the designated target in the pure dew drug state, the decoction drug state, and the solid drug state.

[0120] Optionally, the extraction module 902 includes:

[0121] A segmentation submodule, used for segmenting the pulse wave signal according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence including segmentation signals having the same number of periods as the pulse wave signal;

[0122] The processing submodule is used to perform compound series conversion processing on each segmentation signal in the segmentation sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0123] Optionally, the processing submodule includes:

[0124] A first processing unit, configured to perform compound series conversion processing on each of the segmented signals in the segmented sequence to obtain harmonic frequency data of the segmented signals;

[0125] The second processing unit is used to normalize the harmonic frequency data of the segmented signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0126] Optionally, the first processing unit includes:

[0127] A calculation subunit, used for calculating the average amplitude and average phase of the pulse wave signal;

[0128] The processing subunit is used to normalize the harmonic frequency data of the segmented signal according to the average amplitude and average phase of the pulse wave signal to obtain the harmonic frequency characteristic data corresponding to the pulse wave signal.

[0129] Optionally, the calculation module 903 includes:

[0130] The first calculation submodule is used to calculate the percentage variable of the harmonic frequency characteristic data to obtain the percentage variable distribution data of the pulse wave signal.

[0131] Optionally, the calculation module 903 includes:

[0132] The second calculation submodule is used to calculate the phase difference variable of the harmonic frequency characteristic data to obtain the phase difference variable distribution data of the pulse wave signal.

[0133] It should be noted that the data center device provided in the embodiment of the present invention can be applied to devices such as smart phones, computers, servers, etc. that can process pulse wave data.

[0134] The data center device provided in the embodiment of the present invention can implement each process implemented by the pulse wave data processing method in the above method embodiment and can achieve the same beneficial effects. To avoid repetition, it will not be described here.

[0135] See also Fig.10 , Fig.10 is a schematic diagram of the structure of an electronic device provided by an embodiment of the present invention, such as Fig.10 As shown, it includes: a memory 1002, a processor 1001, and a computer program of a pulse wave data processing method stored in the memory 1002 and executable on the processor 1001, wherein:

[0136] The electronic device is applied to a data center, and the processor 1001 is used to call the computer program stored in the memory 1002 to perform the following steps:

[0137] Acquire a pulse wave signal of a designated target in a target drug state;

[0138] Extracting harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal;

[0139] The variables of the harmonic frequency characteristic data are calculated to obtain the variable distribution data of the pulse wave signal.

[0140] Optionally, the target drug state executed by the processor 1001 includes three drug states corresponding to pure dew, decoction, and solid matter, and obtaining the pulse wave signal of the specified target in the target drug state includes:

[0141] The pulse wave signals corresponding to the designated target in the pure dew drug state, the decoction drug state, and the solid drug state are obtained.

[0142] Optionally, the processor 1001 performs harmonic frequency feature extraction on the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal, including:

[0143] Dividing the pulse wave signal according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence including segmentation signals having the same number of periods as the pulse wave signal;

[0144] Compound series conversion is performed on each segmentation signal in the segmentation sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0145] Optionally, the processor 1001 performs compound series conversion processing on each of the segmented signals in the segmented sequence to obtain harmonic frequency characteristic data corresponding to the pulse wave signal, including:

[0146] Performing compound series conversion processing on each of the segmentation signals in the segmentation sequence to obtain harmonic frequency data of the segmentation signal;

[0147] The harmonic frequency data of the segmented signal is normalized to obtain harmonic frequency characteristic data corresponding to the pulse wave signal.

[0148] Optionally, the processor 1001 performs compound series conversion processing on each of the segmented signals in the segmented sequence to obtain harmonic frequency data of the segmented signals, including:

[0149] Calculating the average amplitude and average phase of the pulse wave signal;

[0150] The harmonic frequency data of the segmented signal is normalized according to the average amplitude and average phase of the pulse wave signal to obtain the harmonic frequency characteristic data corresponding to the pulse wave signal.

[0151] Optionally, the processor 1001 performs the step of calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal, including:

[0152] The percentage variable of the harmonic frequency characteristic data is calculated to obtain the percentage variable distribution data of the pulse wave signal.

[0153] Optionally, the processor 1001 performs the step of calculating the variables of the harmonic frequency characteristic data to obtain the variable distribution data of the pulse wave signal, including:

[0154] The phase difference variable of the harmonic frequency characteristic data is calculated to obtain the phase difference variable distribution data of the pulse wave signal.

[0155] It should be noted that the electronic device provided by the embodiment of the present invention can be applied to devices such as smart phones, computers, servers, etc. that can process pulse wave data.

[0156] The electronic device provided by the embodiment of the present invention can implement each process implemented by the pulse wave data processing method in the above method embodiment and can achieve the same beneficial effects. To avoid repetition, it will not be described here.

[0157] An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the various processes of the pulse wave data processing method or the application-end pulse wave data processing method provided in the embodiment of the present invention are implemented, and the same technical effect can be achieved. To avoid repetition, it will not be repeated here.

[0158] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM).

[0159] The above disclosure is only the preferred embodiment of the present invention, which certainly cannot be used to limit the scope of the present invention. Therefore, equivalent changes made according to the claims of the present invention are still within the scope of the present invention.

Claims

1. A pulse wave data processing method, It is characterized in that The following steps are involved: Acquire the pulse wave signal of the designated target in the target drug state; the target drug state includes three drug states corresponding to pure dew, decoction, and solid matter; Extracting harmonic frequency characteristics of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal; specifically dividing the pulse wave signal according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence including segmentation signals having the same number of periods as the pulse wave signal; performing compound series conversion processing on each segmentation signal in the segmentation sequence to obtain harmonic frequency data of the segmentation signal; calculating the average amplitude and average phase of the pulse wave signal; normalizing the harmonic frequency data of the segmentation signal according to the average amplitude and average phase of the pulse wave signal to obtain harmonic frequency characteristic data corresponding to the pulse wave signal; The percentage variables of the harmonic frequency characteristic data are calculated to obtain the percentage variable distribution data of the pulse wave signal. The percentage variables of the harmonic frequency characteristic data include the harmonic frequency phase average value variable before and after taking the pure dew; the harmonic frequency phase average value variable before and after taking the decoction; and the harmonic frequency phase average value variable before and after taking the solid substance.

2. The pulse wave data processing method according to claim 1, It is characterized in that Obtain the pulse wave signal of the specified target in the target drug state, including: The pulse wave signals corresponding to the designated target in the pure dew drug state, the decoction drug state, and the solid drug state are obtained.

3. The pulse wave data processing method according to claim 1, It is characterized in that The phase difference variable of the harmonic frequency characteristic data is calculated to obtain the phase difference variable distribution data of the pulse wave signal.

4. A pulse wave data processing device, It is characterized in that The pulse wave data processing device comprises: An acquisition module is used to acquire the pulse wave signal of a designated target in a target drug state; the target drug state includes three drug states corresponding to pure dew, decoction, and solid matter; An extraction module is used to extract harmonic frequency features from the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal; specifically, the pulse wave signal is segmented according to the signal period to obtain a segmentation sequence of the pulse wave signal, each segmentation sequence includes segmentation signals having the same number of periods as the pulse wave signal; each segmentation signal in the segmentation sequence is subjected to compound series conversion processing to obtain harmonic frequency data of the segmentation signal; an average amplitude and an average phase of the pulse wave signal are calculated; and the harmonic frequency data of the segmentation signal is normalized according to the average amplitude and the average phase of the pulse wave signal to obtain harmonic frequency feature data corresponding to the pulse wave signal; The calculation module is used to calculate the percentage variable of the harmonic frequency characteristic data to obtain the percentage variable distribution data of the pulse wave signal. The percentage variable of the harmonic frequency characteristic data includes the harmonic frequency phase average value variable before and after taking the pure dew; the harmonic frequency phase average value variable before and after taking the soup; and the harmonic frequency phase average value variable before and after taking the solid matter.

5. An electronic device, It is characterized in that include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps in the pulse wave data processing method according to any one of claims 1 to 3 when executing the computer program.

6. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the pulse wave data processing method according to any one of claims 1 to 3 are implemented.

Images