[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN112155523A - Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification - Google Patents

Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification Download PDF

Info

Publication number
CN112155523A
CN112155523A CN202011035168.8A CN202011035168A CN112155523A CN 112155523 A CN112155523 A CN 112155523A CN 202011035168 A CN202011035168 A CN 202011035168A CN 112155523 A CN112155523 A CN 112155523A
Authority
CN
China
Prior art keywords
modal
signal
component
signals
follows
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011035168.8A
Other languages
Chinese (zh)
Other versions
CN112155523B (en
Inventor
吕玉祥
张琦
李广
朱中艳
胡智君
崔程
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanxi Hi Tan Ke Technology Co ltd
Original Assignee
Taiyuan University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Taiyuan University of Technology filed Critical Taiyuan University of Technology
Priority to CN202011035168.8A priority Critical patent/CN112155523B/en
Publication of CN112155523A publication Critical patent/CN112155523A/en
Application granted granted Critical
Publication of CN112155523B publication Critical patent/CN112155523B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4854—Diagnosis based on concepts of traditional oriental medicine
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Physiology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Cardiology (AREA)
  • Alternative & Traditional Medicine (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The invention relates to a pulse signal feature extraction and classification method based on modal energy principal component ratio quantification, belonging to the technical field of pulse signal feature extraction and classification; the technical problem to be solved is as follows: the improvement of a pulse signal characteristic extraction and classification method based on modal energy principal component ratio quantification is provided; the technical scheme for solving the technical problem is as follows: the method comprises the following steps: acquiring a section of normal human pulse signals through a pulse sensor, carrying out CEEMDAN self-adaptive full integration empirical mode decomposition on the acquired signals, decomposing original complex signals into mode components with different time scales, further decomposing the decomposed signals based on a Mode Energy Principal Component Ratio (MEPCR), and quantizing the decomposed signals for describing the energy principal component distribution of pulses; the invention is applied to pulse signal feature extraction and classification.

Description

Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification
Technical Field
The invention discloses a pulse signal feature extraction and classification method based on modal energy principal component ratio quantization, and belongs to the technical field of pulse signal feature extraction and classification.
Background
Traditional Chinese medicine and modern medicine show that the pulse signals of human bodies contain rich health condition information of the human bodies, in recent years, many researchers are dedicated to detection and analysis of the pulse signals of the human bodies, and diagnosis methods for extracting and classifying various pulse signal characteristics are proposed and applied; the existing pulse signal processing mainly comprises the steps of carrying out characteristic quantization on signals and then classifying the signals according to characteristic vectors, the existing extraction and classification method can basically achieve the classification effect, but the classification efficiency and the effectiveness are still insufficient, and the method process needs to be improved and optimized; because the pulse signal is a typical nonlinear and non-stationary signal, the oscillation mode of the pulse signal is complex, and the problems of low accuracy and low robustness can occur when the signal is classified by adopting a single characteristic.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention aims to solve the technical problems that: the improvement of the pulse signal feature extraction and classification method based on modal energy principal component ratio quantification is provided.
In order to solve the technical problems, the invention adopts the technical scheme that: a pulse signal feature extraction and classification method based on modal energy principal component ratio quantification comprises the following steps:
the method comprises the following steps: acquiring a section of normal human body pulse signals through a pulse sensor, wherein the length of the signals is required to be not less than 3 periods; step two: performing CEEMDAN self-adaptive full-integration empirical mode decomposition on the acquired signals, and decomposing original complex signals into modal components with different time scales;
before decomposition, an operator E is definedk(. to solve the k-th modal component IMF of the EMD decompositionk；
Definition of ωj(t) white noise satisfying the standard normal distribution added in the ith experiment;
definition ofkThe amplitude coefficient of the white noise added for the Kth time;
defining X (t) as an original signal sequence;
step 2.1: to signal X (t) +0ωi(t) performing I times of tests, and decomposing by an EMD method to obtain a first modal component, wherein the calculation formula of the first modal component is as follows:
Figure BDA0002704837180000011
step 2.2: calculating a unique residual signal of a first stage, wherein the calculation formula of the unique residual signal is as follows:
r1(t)＝X(t)-IMF1；
step 2.3: constructing a signal r1(t)+1E1(ωi(t)), then performing EMD, calculating a second modal component, the second modal component being calculated by the formula:
Figure BDA0002704837180000021
step 2.4: in each subsequent stage, the Kth residual component r is calculatedk(t)＝rk-1(t)-IMFkConstructing a new signal, executing EMD, and obtaining a K +1 mode component, wherein a calculation formula of the K +1 mode component is as follows:
Figure BDA0002704837180000022
step 2.5: repeating the step 2.4 until the value of the residual component is less than the two extreme values, and stopping decomposition;
finally obtaining the residual component
Figure BDA0002704837180000023
The original signal is finally decomposed into K modal components, which are expressed as:
Figure BDA0002704837180000024
step three: calculating a Modal Energy Principal Component Ratio (MEPCR), further processing the decomposed signals, and quantizing the signals for describing the energy principal component distribution of the pulse;
step 3.1: decomposing an original signal into K IMF components through a CEEMDAN algorithm, performing correlation calculation between each component and the original signal, and screening out IMF components close to the original signal, wherein the specific calculation steps are as follows:
correlation coefficient
Figure BDA0002704837180000025
The calculation formula of (2) is as follows:
Figure BDA0002704837180000026
wherein Cov (x, IMF)i) IMF for original sequence and modal componentsiCovariance of (a)xFor the original sequence variance, σ IMFiIs the modal component variance;
the fractions were screened based on the above formula, with the screening rules defined as follows:
{IMFj：pxIMFi≧ 0.5}, where j ∈ i, i ═ 1, 2,. k;
it is known that, in signal correlation, a correlation coefficient of 0.5 or more is significant correlation;
step 3.2: the calculation formula of the modal energy principal component ratio MEPCR is as follows:
Figure BDA0002704837180000027
in the formula:
Figure BDA0002704837180000031
compared with the prior art, the invention has the beneficial effects that: the invention provides a new pulse parameter to quantify the signal characteristics, so that the pulse parameter is applied to the extraction and classification of pulse signals; the newly proposed pulse parameter MEPCR has good signal characteristic description effect in the aspect of signal energy distribution, and especially plays an important role in extracting and analyzing nonlinear and non-stable pulse signal characteristics; corresponding experiments are carried out to verify that the new parameter MEPCR has obvious effect on characterization of the pulse signal features, and in practical application, the information contained in the pulse signal can be further comprehensively and comprehensively quantized by using a combination mode of the traditional features and the MEPCR so as to improve the consistency and reliability of the pulse signal diagnosis result.
Drawings
The invention is further described below with reference to the accompanying drawings:
FIG. 1 is a flowchart illustrating a process of extracting pulse signal features according to the present invention;
FIG. 2 is a diagram illustrating the effect of decomposing a pulse signal by an adaptive fully integrated empirical mode according to the present invention;
FIG. 3 is a diagram illustrating the correlation between components and their respective original signals according to an embodiment of the present invention;
FIG. 4 is a graph comparing pre-meal post-meal MEPCR values according to embodiments of the present invention;
FIG. 5 is a scatter plot of MEPCR parameters after a meal for all pulse signals of a data set according to an embodiment of the present invention;
FIG. 6 is an analysis diagram of a classifier model constructed according to an embodiment of the present invention.
Detailed Description
As shown in fig. 1, the present invention relates to the field of signal feature extraction and classification, and proposes a new pulse parameter to quantify a signal feature for accurate classification of a pulse signal, where the pulse parameter is specifically a modal energy principal component ratio MEPCR, and can be applied to pulse signal feature extraction, and the specific steps are as follows:
(1) signal acquisition:
firstly, a section of normal human pulse signals are collected, the length of the signals is required to be 3-5 periods, and in a specific experiment, the pulse signals are collected through an HK-2000C pulse sensor.
(2) Performing CEEMDAN self-adaptive fully-integrated empirical mode decomposition on the acquired signals:
the main function of adopting CEEMDAN decomposition is to decompose original complex signals into modal components with different time scales, and the specific steps of performing CEEMDAN decomposition on the signals are as follows:
first, an operator E is definedk(. the effect of which is to solve the k-th modal component IMF of the EMD decompositionk；
Let omegai(t) white noise satisfying the standard normal distribution added in the ith experiment,kfor the amplitude coefficient of white noise added at the kth time, X (t) is the original signal sequence:
step 2.1: to signal X (t) +0ωi(t) carrying out I times of tests, decomposing by an EMD method to obtain a first modal component, wherein the calculation formula is as follows:
Figure BDA0002704837180000041
step 2.2: and calculating the unique residual signal of the first stage by the following formula:
r1(t)＝X(t)-IMF1；
step 2.3: constructing a signal r1(t)+1E1(ωi(t)), then performing EMD, calculating a second modal component, the calculation formula being:
Figure BDA0002704837180000042
step 2.4: at each stage, the Kth residual component r is calculatedk(t)＝rk-1(t)-IMFkConstructing a new signal, and executing EMD, wherein the calculation formula of the K +1 mode component is as follows:
Figure BDA0002704837180000043
step 2.5: repeating the step 2.4 until the value of the residual component is less than the two extreme values, and stopping decomposition; finally obtaining the residual component
Figure BDA0002704837180000044
The original signal is finally decomposed into K modal components, expressed as:
Figure BDA0002704837180000045
the signal decomposition diagram is shown in fig. 2;
(3) further processing the decomposed signal and quantizing;
when the physiological condition of a human body changes, the energy of certain frequency bands in the human body pulse signal changes, so that the energy distribution is influenced, and therefore the energy characteristics can sensitively detect the tiny change of the pulse signal, and the pulse signal analysis is facilitated.
The invention further provides a novel pulse parameter, namely a Principal component ratio of Modal Energy (MEPCR) based on the CEEMDAN theory, which is used for describing the energy Principal component distribution of the pulse.
The specific calculation method is as follows:
step 3.1: screening main components: decomposing K IMF components of the original signal through a CEEMDAN algorithm, performing correlation calculation between each component and the original signal, and screening out IMF components close to the original signal.
The specific implementation steps are as follows: calculating a correlation coefficient
Figure BDA0002704837180000046
ComputingThe formula is as follows:
Figure BDA0002704837180000047
wherein Cov (x, IMF)i) IMF for original sequence and modal componentsiCovariance of (a)xIn the form of the variance of the original sequence,
Figure BDA0002704837180000051
is the modal component variance. Then, the components are screened according to the following screening rules:
{IMFj：pxIMFi≧ 0.5}, where j ∈ i, i ═ 1, 2,. k;
it is known that a correlation number of 0.5 or more is significant in signal correlation.
Step 3.2: the modal energy principal component ratio MEPCR is:
Figure BDA0002704837180000052
wherein
Figure BDA0002704837180000053
(4) And (3) experimental verification:
in order to evaluate the effectiveness of the new characteristics, the invention collects and analyzes the pulse signals of the same person after meals under the same measuring condition, and respectively calculates new parameters of the two groups of signals according to the steps, and the specific process is as follows:
firstly, two groups of signals are subjected to CEEMDAN decomposition, the standard deviation of the added noise is 0.2, the noise adding times are 500, the maximum iteration times are 240, and the two groups of signals are decomposed into 8 IMF components and a residual component through CEEMDAN.
Then, MEPCR calculations were performed:
component screening is performed first, and correlations between the two sets of signals and respective modal components are calculated according to the aforementioned method, and the correlations between the two sets of components and respective original signals are shown in fig. 3.
And (4) performing relevance screening on the components according to the screening rule, wherein the IMF components corresponding to the circular labeling part in the figure 3 are screened IMF components.
The table below shows the results of screening two sets of modal components of the signal.
Figure BDA0002704837180000054
After component screening is completed, MEPCR calculation is carried out according to steps, a comparison graph of the pre-meal and post-meal MEPCR values is shown in FIG. 4, as can be seen from FIG. 4, the MEPCR parameters of the pre-meal and post-meal pulse signals of the same person are obviously different, and the MEPCR of the post-meal pulse signals is obviously higher than the MEPCR parameters of the pre-meal pulse signals.
In order to verify the effectiveness and universality of the feature in the pulse signal classification process, the embodiment constructs a data set containing 200 pulses, the data source is 200 volunteers aged 23-26 years and healthy, and pulse signals of experimenters are respectively collected before and after breakfast, and the specific data set composition information is shown in the following table:
Figure BDA0002704837180000061
after data acquisition is finished, performing MEPCR parameter calculation on all pulse signals according to steps, and then analyzing results.
FIG. 5 is a scatter plot of the MEPCR parameters of all the pulse signals in the data set, and it can be seen from FIG. 5 that the MEPCR parameters of the human body pre-meal post-meal pulse signals are obviously different and regularly distributed, the MEPCR parameters of the pre-meal pulse signals (star icon) are obviously smaller than the MEPCR values of the post-meal pulse signals (diamond icon), the MEPCR values of the pre-meal pulse signals are distributed between 0.5 and 0.75, and the MEPCR values of the post-meal pulse signals are basically distributed above 0.78. The quantitative effect of the integral characteristics of the MEPCR parameters is obvious.
To further verify the pulse signal validity, the present invention will use the above data set for classifier construction:
firstly, sample marking and label setting are carried out:
wherein the MEPCRqMEPCR value, MEPCR for Pre-meal SignalhMEPCR value as postprandial signal;
mix MEPCRqSet as tag 1, MEPCRhSet as tag 2;
when the classifier is constructed, a Cubic SVM is selected as a classification model, five-layer cross inspection is set, the model construction and analysis results are shown in fig. 6, a triangular icon is shown in the graph as a label 2, a circular icon is shown in the graph as a label 1, the recognition accuracy of the trained model signal is known to be 95.0%, and the result shows that the new parameter feature MEPCR feature is effective.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (1)

1. A pulse signal feature extraction and classification method based on modal energy principal component ratio quantification is characterized in that: the method comprises the following steps:
the method comprises the following steps: acquiring a section of normal human body pulse signals through a pulse sensor, wherein the length of the signals is required to be not less than 3 periods;
step two: performing CEEMDAN self-adaptive full-integration empirical mode decomposition on the acquired signals, and decomposing original complex signals into modal components with different time scales;
before decomposition, an operator E is definedk(. to solve the k-th modal component IMF of the EMD decompositionk；
Definition of ωi(t) white noise satisfying the standard normal distribution added in the ith experiment;
definition ofkThe amplitude coefficient of the white noise added for the Kth time;
defining X (t) as an original signal sequence;
step 2.1: to signal X (t) +0ωi(t) performing I times of tests, and decomposing by an EMD method to obtain a first modal component, wherein the calculation formula of the first modal component is as follows:
Figure FDA0002704837170000011
step 2.2: calculating a unique residual signal of a first stage, wherein the calculation formula of the unique residual signal is as follows:
r1(t)＝X(t)-IMF1；
step 2.3: constructing a signal r1(t)+1E1(ωi(t)), then performing EMD, calculating a second modal component, the second modal component being calculated by the formula:
Figure FDA0002704837170000012
step 2.4: in each subsequent stage, the Kth residual component r is calculatedk(t)＝rk-1(t)-IMFkConstructing a new signal, executing EMD, and obtaining a K +1 mode component, wherein a calculation formula of the K +1 mode component is as follows:
Figure FDA0002704837170000013
step 2.5: repeating the step 2.4 until the value of the residual component is less than the two extreme values, and stopping decomposition;
finally obtaining the residual component
Figure FDA0002704837170000014
The original signal is finally decomposed into K modal components, which are expressed as:
Figure FDA0002704837170000015
step three: calculating a Modal Energy Principal Component Ratio (MEPCR), further processing the decomposed signals, and quantizing the signals for describing the energy principal component distribution of the pulse;
step 3.1: decomposing an original signal into K IMF components through a CEEMDAN algorithm, performing correlation calculation between each component and the original signal, and screening out IMF components close to the original signal, wherein the specific calculation steps are as follows:
correlation coefficient
Figure FDA0002704837170000026
The calculation formula of (2) is as follows:
Figure FDA0002704837170000021
wherein Cov (x, IMF)i) IMF for original sequence and modal componentsiCovariance of (a)xIn the form of the variance of the original sequence,
Figure FDA0002704837170000022
is the modal component variance;
the fractions were screened based on the above formula, with the screening rules defined as follows:
Figure FDA0002704837170000023
wherein j ∈ i, i ═ 1, 2,. k;
it is known that, in signal correlation, a correlation coefficient of 0.5 or more is significant correlation;
step 3.2: the calculation formula of the modal energy principal component ratio MEPCR is as follows:
Figure FDA0002704837170000024
in the formula:
Figure FDA0002704837170000025
CN202011035168.8A 2020-09-27 2020-09-27 Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification Active CN112155523B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011035168.8A CN112155523B (en) 2020-09-27 2020-09-27 Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011035168.8A CN112155523B (en) 2020-09-27 2020-09-27 Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification

Publications (2)

Publication Number Publication Date
CN112155523A true CN112155523A (en) 2021-01-01
CN112155523B CN112155523B (en) 2022-09-16

Family

ID=73861361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011035168.8A Active CN112155523B (en) 2020-09-27 2020-09-27 Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification

Country Status (1)

Country Link
CN (1) CN112155523B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114587291A (en) * 2022-03-14 2022-06-07 常熟理工学院 Pulse rate variability signal extraction method and system

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1076586A1 (en) * 1998-05-06 2001-02-21 Exogen, Inc. Ultrasound bandages
EP1482142A2 (en) * 2001-08-06 2004-12-01 Southwest Research Institute Method and apparatus for testing catalytic converter durability
US20110004110A1 (en) * 2000-05-30 2011-01-06 Vladimir Shusterman Personalized Monitoring and Healthcare Information Management Using Physiological Basis Functions
US20110057818A1 (en) * 2006-01-18 2011-03-10 Lg Electronics, Inc. Apparatus and Method for Encoding and Decoding Signal
US20110158202A1 (en) * 1995-06-30 2011-06-30 Interdigital Technology Corporation Code division multiple access (cdma) communication system
US20110245628A1 (en) * 2010-03-31 2011-10-06 Nellcor Puritan Bennett Llc Photoplethysmograph Filtering Using Empirical Mode Decomposition
US20110270579A1 (en) * 2010-04-28 2011-11-03 Nellcor Puritan Bennett Ireland Systems and methods for signal monitoring using lissajous figures
US20120138805A1 (en) * 2010-12-02 2012-06-07 Xtreme Technologies Gmbh Method for the Spatially Resolved Measurement of Parameters in a Cross Section of a Beam Bundle of High-Energy Radiation of High Intensity
CN103496625A (en) * 2013-10-17 2014-01-08 太原理工大学 Multi-rope friction lifter load identification method based on vibration analysis
CN104414688A (en) * 2013-08-23 2015-03-18 北京大学 Ensemble empirical mode decomposition-based vasovagal syncope precursor detection method
CN104921715A (en) * 2015-06-09 2015-09-23 上海华旌科技有限公司 Multi-parameter vital sign measurement device
US20170135644A1 (en) * 2015-11-18 2017-05-18 National Chiao Tung University Physiological signal measuring system and method thereof
CN108013548A (en) * 2018-01-24 2018-05-11 太原理工大学 Can pinpoint human body physical sign monitoring Intelligent bracelet
US20190021615A1 (en) * 2017-07-18 2019-01-24 STMicroelecteronics s.r.l. Processing of Electrophysiological Signals
CN109498041A (en) * 2019-01-15 2019-03-22 吉林大学 Driver road anger state identification method based on brain electricity and pulse information
US10357163B1 (en) * 2012-06-01 2019-07-23 Vital Connect, Inc. Respiratory rate detection using decomposition of ECG
CN110090010A (en) * 2019-06-17 2019-08-06 北京心数矩阵科技有限公司 A kind of contactless blood pressure measuring method and system
CN110309817A (en) * 2019-07-19 2019-10-08 北京理工大学 A kind of pulse wave motion artifacts minimizing technology of parameter adaptive optimization VMD
CN111243739A (en) * 2020-01-07 2020-06-05 四川大学 Anti-interference physiological parameter telemetering method and system
EP3669757A1 (en) * 2018-12-18 2020-06-24 Koninklijke Philips N.V. System and method for detecting fluid accumulation

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110158202A1 (en) * 1995-06-30 2011-06-30 Interdigital Technology Corporation Code division multiple access (cdma) communication system
EP1076586A1 (en) * 1998-05-06 2001-02-21 Exogen, Inc. Ultrasound bandages
US20110004110A1 (en) * 2000-05-30 2011-01-06 Vladimir Shusterman Personalized Monitoring and Healthcare Information Management Using Physiological Basis Functions
EP1482142A2 (en) * 2001-08-06 2004-12-01 Southwest Research Institute Method and apparatus for testing catalytic converter durability
US20110057818A1 (en) * 2006-01-18 2011-03-10 Lg Electronics, Inc. Apparatus and Method for Encoding and Decoding Signal
US20110245628A1 (en) * 2010-03-31 2011-10-06 Nellcor Puritan Bennett Llc Photoplethysmograph Filtering Using Empirical Mode Decomposition
US20110270579A1 (en) * 2010-04-28 2011-11-03 Nellcor Puritan Bennett Ireland Systems and methods for signal monitoring using lissajous figures
US20120138805A1 (en) * 2010-12-02 2012-06-07 Xtreme Technologies Gmbh Method for the Spatially Resolved Measurement of Parameters in a Cross Section of a Beam Bundle of High-Energy Radiation of High Intensity
US10357163B1 (en) * 2012-06-01 2019-07-23 Vital Connect, Inc. Respiratory rate detection using decomposition of ECG
CN104414688A (en) * 2013-08-23 2015-03-18 北京大学 Ensemble empirical mode decomposition-based vasovagal syncope precursor detection method
CN103496625A (en) * 2013-10-17 2014-01-08 太原理工大学 Multi-rope friction lifter load identification method based on vibration analysis
CN104921715A (en) * 2015-06-09 2015-09-23 上海华旌科技有限公司 Multi-parameter vital sign measurement device
US20170135644A1 (en) * 2015-11-18 2017-05-18 National Chiao Tung University Physiological signal measuring system and method thereof
US20190021615A1 (en) * 2017-07-18 2019-01-24 STMicroelecteronics s.r.l. Processing of Electrophysiological Signals
CN108013548A (en) * 2018-01-24 2018-05-11 太原理工大学 Can pinpoint human body physical sign monitoring Intelligent bracelet
EP3669757A1 (en) * 2018-12-18 2020-06-24 Koninklijke Philips N.V. System and method for detecting fluid accumulation
CN109498041A (en) * 2019-01-15 2019-03-22 吉林大学 Driver road anger state identification method based on brain electricity and pulse information
CN110090010A (en) * 2019-06-17 2019-08-06 北京心数矩阵科技有限公司 A kind of contactless blood pressure measuring method and system
CN110309817A (en) * 2019-07-19 2019-10-08 北京理工大学 A kind of pulse wave motion artifacts minimizing technology of parameter adaptive optimization VMD
CN111243739A (en) * 2020-01-07 2020-06-05 四川大学 Anti-interference physiological parameter telemetering method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
DANHUA ZHOU: "State of Health Monitoring and Remaining Useful Life Prediction of Lithium-Ion Batteries Based on Temporal Convolutional Network", 《IEEE ACCESS》 *
KANG SUN: "Extraction of partial discharge pulses from the complex noisy signals of power cables based on CEEMDAN and wavelet packetpower", 《ENERGIES》 *
张兴: "光电容积脉搏信号特征提取与分类识别研究", 《中国优秀硕士学位论文全文数据库》 *
杨杰: "新型三部脉象采集系统与脉象信息分析研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114587291A (en) * 2022-03-14 2022-06-07 常熟理工学院 Pulse rate variability signal extraction method and system

Also Published As

Publication number Publication date
CN112155523B (en) 2022-09-16

Similar Documents

Publication Publication Date Title
US8315970B2 (en) Method and rhythm extractor for detecting and isolating rhythmic signal features from an input signal using the wavelet packet transform
Khan et al. Automated classification of lung sound signals based on empirical mode decomposition
CN103310113B (en) A kind ofly to be separated and the general blood glucose prediction method of data modeling based on frequency band
CN109602417A (en) Sleep stage method and system based on random forest
Güler et al. A modified mixture of experts network structure for ECG beats classification with diverse features
CN113069117A (en) Electroencephalogram emotion recognition method and system based on time convolution neural network
CN108567418A (en) A kind of pulse signal inferior health detection method and detecting system based on PCANet
CN103690160A (en) Electroencephalogram feature extraction method based on non-Gaussian time sequence model
CN117349741A (en) Raman spectrum classification method, species blood semen and species classification method
CN112155523B (en) Pulse signal feature extraction and classification method based on modal energy principal component ratio quantification
CN114113471A (en) Method and system for detecting food freshness of artificial nose refrigerator based on machine learning
CN107045624B (en) Electroencephalogram signal preprocessing and classifying method based on maximum weighted cluster
CN113796873B (en) Wearable dynamic electrocardiosignal classification method and system
CN114947850A (en) Mental load grade objective detection method based on pulse Bouss model characteristics
Liu et al. Evaluating and improving automatic sleep spindle detection by using multi-objective evolutionary algorithms
Nabilah et al. Human bone age estimation of carpal bone X-ray using residual network with batch normalization classification
Hu et al. A novel nonlinear bispectrum analysis for dynamical complex oscillations
CN112378942B (en) White spirit grade classification and identification method based on nuclear magnetic resonance fingerprint
CN115269679A (en) Multidimensional time series overall complexity evaluation method
CN110327043B (en) Event-related potential waveform map solving method based on sparse modeling
CN112200228A (en) Epileptic seizure state identification method based on two-dimensional convolutional neural network
CN111631711B (en) Method for analyzing and processing brain electrical data of schizophrenia and auditory hallucination symptoms
Huang et al. Bandpass empirical mode decomposition using a rolling ball algorithm
Navarro et al. Performance analysis of hurst's exponent estimators in higly immature breathing patterns of preterm infants
CN116712089B (en) Epileptiform discharge enriching epileptiform interval and method for predicting focus

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20230517

Address after: AW03, AW08, 4th Floor, Building 4, No. 11 Kangshou Street, Tanghuai Industrial Park, Shanxi Transformation and Comprehensive Reform Demonstration Zone, Taiyuan City, Shanxi Province, 030032 (China Shanxi Overseas Students Entrepreneurship Park)

Patentee after: Shanxi Hi Tan Ke Technology Co.,Ltd.

Address before: 030024 No. 79 West Main Street, Taiyuan, Shanxi, Yingze

Patentee before: Taiyuan University of Technology

TR01 Transfer of patent right