CN105476630A - Method for detecting traditional Chinese medicine visceral functions - Google Patents
Method for detecting traditional Chinese medicine visceral functions Download PDFInfo
- Publication number
- CN105476630A CN105476630A CN201510846258.8A CN201510846258A CN105476630A CN 105476630 A CN105476630 A CN 105476630A CN 201510846258 A CN201510846258 A CN 201510846258A CN 105476630 A CN105476630 A CN 105476630A
- Authority
- CN
- China
- Prior art keywords
- experimenter
- value range
- normal variation
- variation value
- heart rate
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention provides a method for detecting traditional Chinese medicine visceral functions, and relates to technologies for detecting indication of human bodies. The method includes acquiring electrocardiogram signals of subjects in preset time periods; extracting normal R-wave vertexes from the electrocardiogram signals and acquiring RR interval time sequences; acquiring data of long axes and short axes of Poincare scatter diagrams according to the RR interval time sequences; acquiring proportions of the long axes to the short axes and/or the areas of the Poincare scatter diagrams by means of computing; determining the visceral controlling blood and circulation functions of the subjects according to the data of the long axes, the data of the short axes, the proportions of the long axes to the short axes and/or the areas of the Poincare scatter diagrams. The method has the advantages that the method is beneficial to objectifying detection of the traditional Chinese medicine visceral functions, and the accuracy of detection of the traditional Chinese medicine visceral functions can be improved.
Description
Technical field
The present invention relates to human body indication detection technique, particularly relate to the detection method that the traditional Chinese medical science heart hides function.
Background technology
Be different from anatomical internal organs, the Tibetan (sound z à ng) of Chinese medicine as if refer to is hidden (sound c á ng) spleen, lung, kidney, liver, the heart etc. in body and is found expression in outer physiological and pathological sign, together with natural correlate and phenomenon, be also called spleen, lung, kidney, liver, the heart five hide system, be called for short five and hide.Wherein, Basic Theories of Chinese Medicine is thought, the heart is hidden has controlling blood and vessel and house mind's two functions.
The controlling blood and vessel function that the heart is hidden and the heart are hidden and are promoted the function run in vascular of blood, and cardiopalmus, uncomfortable in chest, the symptom such as to breathe hard can appear in this dysfunction, the main cardiac function corresponding to modern medicine.At present, cardiac function often carries out NYHA (New York Heart disease association) cardiac functional grading according to the physical limited degree of patient and determines; But the accuracy of this classification determination cardiac function is not high.Theory of Chinese medical science is thought: house mind's function that the heart is hidden refers to that the heart hides the ergasias such as main department consciousness, thinking, emotion, the exception of this function there will be the symptoms such as insomnia, dreaminess, lassitude, coma, wherein daytime the mental status and sleeping at night quality be important the measuring that the heart hides house mind function; At present, sleep quality is determined by Pittsburgh evaluation form, and the mental status is then determined by fatigue scale.But, the controlling blood and vessel function of in prior art, the heart being hidden and the detection subjectivity of house mind's function strong, its accuracy is very low.
Summary of the invention
For above-mentioned defect of the prior art, the invention provides the detection method that a kind of traditional Chinese medical science heart hides function, the detection that the method contributes to making the traditional Chinese medical science heart hide function objectifies, and can improve the accuracy traditional Chinese medical science heart being hidden to Function detection.
The invention provides the detection method that a kind of traditional Chinese medical science heart hides function, comprising: the step gathering the electrocardiosignal of experimenter in preset time period; From described electrocardiosignal, extract normal R wave crest point, and obtain the step of RR time interval sequence; The major axis of Poincare scatterplot and the step of minor axis data is obtained according to described RR time interval sequence; And calculate the acquisition ratio of described minor axis with major axis and/or the step of described Poincare scatterplot area.
Further, describedly from described electrocardiosignal, extract normal R wave crest point comprise: according to mobile minute window method and default window length, a minute window is carried out to described electrocardiosignal; Extract the summit of described electrocardiosignal in each window; Obtain slope characteristics and the R wave crest point slope characteristics on described summit; When the slope characteristics on described summit is mated with described R wave crest point slope characteristics, using described summit as R wave crest point.
Further, described method also comprises: the step obtaining the ratio of described major axis and/or minor axis and/or minor axis and major axis and/or the normal variation value range of described Poincare scatterplot area.
Further, described method also comprises: according to the step of the normal variation value range of area MSEQ and MSEQ of described RR time interval sequence determination multi-scale entropy; And/or according to the normal variation value range of described MSEQ and MSEQ, the heart obtaining described experimenter hides the step of controlling blood and vessel functional status information.
Further, described method also comprises: the step of described RR time interval sequence being carried out to moving average process; According to carrying out the RR time interval sequence after moving average, obtain the step of the changes in heart rate curve of described experimenter in preset time period; According to carrying out the RR time interval sequence after moving average, obtain tested recovery time and/or the step of the tested length of one's sleep of described experimenter; According to carrying out the RR time interval sequence after moving average, obtain the average heart rate of described experimenter within the described tested length of one's sleep and/or the step of the normal variation value range of average heart rate in the length of one's sleep; And/or according to carrying out the RR time interval sequence after moving average, obtain that heart rate criteria in the average heart rate of described experimenter in tested recovery time, tested recovery time is poor, the step of the normal variation value range of heart rate standard deviation in the normal variation value range of average heart rate and/or recovery time in recovery time.
Further, described method also comprises: according to described RR time interval sequence, is obtained the step of the normal variation value range of the normal variation value range of SDANN, SDNNID, rMSSD, PNN50, SDANN of described experimenter, the normal variation value range of SDNNID, the normal variation value range of rMSSD and/or PNN50 by time-domain analysis; Wherein, the standard deviation of the described SDANN every 5 minutes RR intervals that are described experimenter in preset time period; The standard deviation of every 5 minutes RR intervals that described SDNNID is described experimenter in preset time period; Described rMSSD is the described experimenter square root that adjacent R R interval difference is mean square in preset time period; Described PNN50 is the number that described experimenter differs >50 millisecond adjacent R R interval in preset time period.
Further, described method also comprises: according to described RR time interval sequence, obtains the HF of described experimenter and/or the step of normal HF excursion by frequency-domain analysis; Wherein, described HF is high frequency power value.
The traditional Chinese medical science heart provided by the invention hides the detection method of function, comprise: the electrocardiosignal gathering experimenter in preset time period, normal R wave crest point is extracted from described electrocardiosignal, and obtain RR time interval sequence, major axis and the minor axis data of Poincare scatterplot are obtained according to described RR time interval sequence, calculate the ratio and/or described Poincare scatterplot area that obtain described minor axis and major axis, with according to these major axis data, minor axis data, ratio and/or the Poincare scatterplot area of minor axis and major axis determine that the heart of experimenter hides controlling blood and vessel function, thus contribute to realizing objectifying of traditional Chinese medical science heart Tibetan Function detection, improve the accuracy traditional Chinese medical science heart being hidden to Function detection.
Accompanying drawing explanation
Fig. 1 is the detection method first pass schematic diagram that the traditional Chinese medical science heart of the present invention hides function;
Fig. 2 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides Poincare scatterplot in the detection method of function;
Fig. 3 is detection method second schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function;
Fig. 4 is detection method the 3rd schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function;
Fig. 5 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides multi-scale entropy in the detection method of function;
Fig. 6 is detection method the 4th schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function;
Fig. 7 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides heart rate change curve in the detection method of function;
Fig. 8 is detection method the 5th schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function.
Detailed description of the invention
Fig. 1 is the detection method first pass schematic diagram that the traditional Chinese medical science heart of the present invention hides function; Fig. 2 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides Poincare scatterplot in the detection method of function; Please refer to Fig. 1-2, the embodiment of the present invention provides a kind of traditional Chinese medical science heart to hide the detection method of function, comprising:
S11, in preset time period, gather the step of the electrocardiosignal of experimenter; Preferably, this preset time period is 24 hours, namely from before first day sleep, terminates before second night sleep.
In the embodiment that the present invention one class is concrete, by the electrocardiosignal of electrocardioelectrode Real-time Collection experimenter, and pretreatment can be carried out to electrocardiosignal; The analog electrocardiogram signal collected can also be converted to digital electrocardiosignal by analog-digital converter, and digital electrocardiosignal is input in data storage stores, detect to hide function according to the electrocardiosignal in preset time period to the heart.After the ecg signal acquiring in preset time period, the electrocardiosignal in preset time period can be saved to corresponding treatment facility, to be further processed electrocardiosignal by the mode of wireless transmit or the mode of data storage unloading storage.Preferably, data storage can adopt flash memory, safe digital card or USB flash disk; Sample frequency can be more than 150 hertz, and preferably more than 250 hertz to improve the quality of electrocardiosignal gathered.
S12, from electrocardiosignal, extract normal R wave crest point, and obtain the step of RR time interval sequence; Wherein, QRS wave group generally includes three ripples be closely connected, and first downward ripple is called Q ripple, and the upright ripple of a high point after Q ripple is called R ripple, and the ripple under R ripple is backward is called S ripple.
The step extracting normal R wave crest point from electrocardiosignal can adopt the recognition methods of existing R wave crest point; Such as: for the Morphological Features of ecg-r wave, such as: the speed rising and decline identifies R wave crest point, a distance i.e. corresponding cardiac cycle of two adjacent R wave crest points is also a RR interval.In the embodiment that a class of the present invention is concrete, the step extracting normal R wave crest point from electrocardiosignal can comprise such operation: carry out a minute window according to mobile minute window method and default window length to electrocardiosignal; Extract the summit of electrocardiosignal in each window; Obtain slope characteristics and the R wave crest point slope characteristics on summit; When the slope characteristics on summit is mated with R wave crest point slope characteristics, using summit as R wave crest point.
Such as: when the sample frequency of electrocardiosignal is 250Hz, also namely 250 points can be gathered in 1 second, can using the length corresponding to 18 points as the window length preset, 1-18 point, 2-19 point, 3-20 point is extracted successively according to a mobile point window method ... the summit of corresponding electrocardiosignal, the slope characteristics on the summit of extracting is mated with R wave crest point slope characteristics, if the slope characteristics on this summit is mated with R wave crest point slope characteristics, then using this summit as R wave crest point; Wherein, slope characteristics can be rising, decline, the slope value of upstroke and the slope value of decent, and R wave crest point slope characteristics can for determining through statistics and being stored in data storage in advance.
Be understandable that: the slope value of the upstroke of a normal R wave crest point is greater than 10, and the absolute value of the slope of decent is greater than 20.
S13, obtain the major axis of Poincare scatterplot and the step of minor axis data according to RR time interval sequence;
S14, calculating obtain the ratio of minor axis with major axis and/or the step of Poincare scatterplot area;
Above-mentioned step S12, S13, S14 can be performed by microprocessor, microprocessor can adopt ARM, and (English full name is AdvancedRISCMachine, Advanced Reduced Instruction Set machine), (English full name is digitalsingnalprocessor to DSP, digital signal processor), (English full name is Field-ProgrammableGateArray to FPGA, field programmable gate array) and MCU (English full name is MicrocontrollerUnit, micro-control unit) in any one.
As shown in Figure 2, because loose some D in Poincare scatterplot mostly concentrates on a certain region, this region is made to have an obvious boundary B S, so, can be obtained the boundary B S of Poincare scatterplot by image processing techniques, the part that boundary B S surrounds is the area of Poincare scatterplot; Become the straight line of miter angle to be the longitudinal axis L A of Poincare scatterplot with transverse axis; Become another straight line of 135 degree to be the minor axis SA of Poincare scatterplot with transverse axis.
The traditional Chinese medical science heart that the embodiment of the present invention provides hides the detection method of function, electrocardiosignal based on experimenter determines the RR time interval sequence of this experimenter, major axis and the minor axis data of Poincare scatterplot are obtained according to this time interval sequence, and calculate the ratio and/or Poincare scatterplot area that obtain minor axis and major axis, with according to these major axis data, minor axis data, ratio and/or the Poincare scatterplot area of minor axis and major axis determine that the heart of experimenter hides controlling blood and vessel function, thus contribute to realizing objectifying of traditional Chinese medical science heart Tibetan Function detection, improve the accuracy traditional Chinese medical science heart being hidden to Function detection.
Fig. 3 is detection method second schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function; Please refer to Fig. 3, the method can also comprise:
S21, obtain the step of major axis and/or the ratio of minor axis and/or minor axis and major axis and/or the normal variation value range of Poincare scatterplot area.Wherein, the normal variation value range of ratio, the normal variation value range of area of the normal variation value range of major axis, the normal variation value range of minor axis, minor axis major axis can be stored in advance in data storage, and the heart that the normal variation value range of the normal variation value range of major axis, minor axis, the normal variation value range of the ratio of minor axis major axis, the normal variation value range of area are Corpus--based Method to be determined hide function normal time major axis span, minor axis span, the ratio span of minor axis major axis and area span.
Above-mentioned step can be performed by microprocessor, and microprocessor can adopt any one in ARM, DSP, FPGA and MCU.In a class detailed description of the invention of the present invention, the ratio of the major axis data obtained, minor axis data, minor axis and major axis and/or Poincare scatterplot area can be contrasted with the normal span of corresponding major axis, the normal span of minor axis, the normal span of the ratio of minor axis major axis and/or the normal span of area respectively; If the ratio of major axis, minor axis, minor axis major axis and/or area are less than the normal variation value range of major axis, the normal variation value range of minor axis, the normal variation value range of the ratio of minor axis major axis and/or the normal variation value range of area respectively, then can determine: the heart of experimenter is hidden controlling blood and vessel function and reduced; Major axis, minor axis or above-mentioned Poincare scatterplot area are less, and the heart Tibetan controlling blood and vessel function of experimenter is lower; If the ratio of minor axis, major axis is greater than the normal variation scope of the ratio of minor axis major axis, then can determine that the heart Tibetan controlling blood and vessel function of experimenter is extremely low, disease enters whole late period; Because the present invention can provide the data of above-mentioned science and determine quantitative analysis conclusion, thus realize and ensure that the traditional Chinese medical science heart hides the quantification detection of controlling blood and vessel function, making Diagnostics of Chinese Medicine be provided with the analysis means of science, quantification.
Fig. 4 is detection method the 3rd schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function; Fig. 5 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides multi-scale entropy in the detection method of function; Please refer to Fig. 4-5, can also comprise after step s 12:
S31, step according to the normal variation value range of area MSEQ and MSEQ of RR time interval sequence determination multi-scale entropy; And/or
S32, normal variation value range according to MSEQ and MSEQ, the heart obtaining experimenter hides the step of controlling blood and vessel functional status information.Wherein, the normal variation value range of MSEQ be that Corpus--based Method is determined, the heart hide function normal time span, can be stored in advance in data storage.
In the embodiment that a class of the present invention is concrete, step S31, S32 can be performed by microprocessor, and microprocessor can adopt any one in ARM, DSP, FPGA and MCU.As shown in Figure 5, transverse axis represents 20 yardstick Scal, and the longitudinal axis represents entropy MSE; The area that forming curves L0, curve L0 and the longitudinal axis, transverse axis surround that coupled together by entropy MSE corresponding to each yardstick Scal is MSEQ.By the normal variation value range comparison of MSEQ and MSEQ, if MSEQ is less than the normal variation scope of MSEQ, then can illustrates and determine that the heart of experimenter is hidden controlling blood and vessel function and reduced; MSEQ is less, and the heart Tibetan controlling blood and vessel function of experimenter is lower.The mode of such employing data analysis, the quantification making the traditional Chinese medical science heart hide controlling blood and vessel function detects and is achieved.
Fig. 6 is detection method the 4th schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function; Fig. 7 is the schematic diagram that the traditional Chinese medical science heart of the present invention hides heart rate change curve in the detection method of function; Please refer to Fig. 6-7, can also comprise after step s 12:
S41, RR time interval sequence is carried out to the step of moving average process;
S42, basis carry out the RR time interval sequence after moving average, obtain the step of the changes in heart rate curve of experimenter in preset time period;
S43, basis carry out the RR time interval sequence after moving average, obtain tested recovery time and/or the step of the tested length of one's sleep of experimenter;
S44, basis carry out the RR time interval sequence after moving average, obtain the step of the normal variation value range of the average heart rate of experimenter within the tested length of one's sleep and/or the average heart rate in the length of one's sleep; Wherein, the normal variation value range of the average heart rate in the length of one's sleep be that Corpus--based Method is determined, the heart hide function normal time span, can be stored in advance in data storage.
In a class detailed description of the invention of the present invention, the normal variation value range of the average heart rate of experimenter within the tested length of one's sleep and the average heart rate in the length of one's sleep is contrasted, if the average heart rate of experimenter within the tested length of one's sleep is greater than the normal variation value range of the average heart rate in the length of one's sleep, then can determine that the sleep of experimenter is more shallow; Or, if experimenter declines comparatively slow at tested sleep initial stage heart rate, then can determine that experimenter has difficulty in going to sleep; Or, if experimenter is in the comparatively morning final heart rate of tested sleep rapid rise time, then can determine experimenter's early awakening; Or, if there are repeatedly the rapid rise and fall of heart rate in the tested length of one's sleep, then can determine that experimenter sleeps middle dreaminess or easily wake up, and achieves the detection by quantitative of the sleep quality to experimenter.
S45, according to carrying out the RR time interval sequence after moving average, obtain that heart rate criteria in the average heart rate of experimenter in tested recovery time, tested recovery time is poor, the step of the normal variation value range of heart rate criteria difference in the normal variation value range of average heart rate in recovery time and/or recovery time.Wherein, the normal variation value range of the heart rate criteria difference in the normal variation value range of the average heart rate in recovery time and recovery time be that Corpus--based Method is determined, the heart hide function normal time span, can be stored in advance in data storage.
In a class detailed description of the invention of the present invention, the average heart rate of experimenter in tested recovery time, heart rate criteria difference are contrasted with the normal variation value range of average heart rate in corresponding recovery time, the normal variation value range of heart rate criteria difference, if the average heart rate of experimenter in tested recovery time is greater than the normal variation value range of the average heart rate in recovery time, and heart rate criteria difference in tested recovery time is comparatively large, then can determine: experimenter is in and is overexcited or is in anxiety state in above-mentioned tested recovery time; Or, if the average heart rate of experimenter in tested recovery time is less than the normal variation scope of the average heart rate in recovery time, and heart rate criteria difference in tested recovery time is less, then determine experimenter's lassitude or be in depressive state, thus achieving the detection by quantitative of the mental status to experimenter.
Each step in the present embodiment can be performed by microprocessor, and microprocessor can adopt any one in ARM, DSP, FPGA and MCU.Such as: as shown in Figure 7, transverse axis represents the detection beginning and ending time of 24 hours, the longitudinal axis represents heart rate, the change of heart rate in region A between transverse axis with straight line L1 corresponding tested length of one's sleep, the change of heart rate in the corresponding tested recovery time of region B between straight line L1 with straight line L2, the change of heart rate during the region C of more than straight line L2 corresponding excitatory state; Heart rate major part corresponding in the T1-T2 time period is all positioned at region A, represent that experimenter is in sleep state, the T1-T2 time period is tested length of one's sleep of experimenter, in this time period, if when namely a certain moment heart rate heart rate be greater than corresponding to L1 is positioned at region B, represents in the sleeping of this moment experimenter and wake up or have a dream; Heart rate major part corresponding in the T2-T3 time period is all positioned at region B, represent that experimenter is in waking state, the T2-T3 time period is the tested recovery time of experimenter, in this time period, if namely a certain moment heart rate is positioned at region A lower than the heart rate corresponding to L1, then represent this moment experimenter's lassitude; If namely a certain moment heart rate is positioned at region C higher than the heart rate corresponding to L2, then represent that this moment experimenter is overexcited.It should be noted that: the moment in figure representated by TI is 22:44 moment of first day, the moment representated by T3 is 22:44 moment of second day.
Fig. 8 is detection method the 5th schematic flow sheet that the traditional Chinese medical science heart of the present invention hides function; Please refer to Fig. 8, also comprise after step s 12:
S51, according to RR time interval sequence, obtained the step of the normal variation value range of the normal variation value range of SDANN, SDNNID, rMSSD, PNN50, SDANN of experimenter, the normal variation value range of SDNNID, the normal variation value range of rMSSD and/or PNN50 by time-domain analysis; Wherein,
The standard deviation of every 5 minutes RR intervals that SDANN is experimenter in preset time period;
The standard deviation of every 5 minutes RR intervals that SDNNID is experimenter in preset time period;
RMSSD is experimenter's square root that adjacent R R interval difference is mean square in preset time period;
PNN50 is the number that experimenter differs >50 millisecond adjacent R R interval in preset time period.
The normal variation value range of SDANN, the normal variation value range of SDNNID, the normal variation value range of rMSSD and the normal variation value range of PNN50 be respectively the heart determined of Corpus--based Method hide function normal time span, can be stored in advance in data storage.The normal variation value range of the normal variation value range of SDANN, SDNNID, rMSSD, PNN50 of experimenter and corresponding SDANN, the normal variation value range of SDNNID, the normal variation value range of rMSSD and PNN50 is contrasted; If the normal variation scope that the SDANN of experimenter is greater than the normal variation scope of SDANN, the SDNNID of experimenter is less than SDNNID, then determine the sympathetic activation of experimenter; If the rMSSD of experimenter is greater than the normal variation scope of rMSSD, the PNN50 of experimenter is greater than the normal variation scope of PNN50, then determine that the parasympathetic nervous of experimenter is more excited, thus achieve the detection by quantitative of the mental status to experimenter.
Can also comprise after step s 12 or after step S51:
S52, according to RR time interval sequence, obtain the HF of experimenter and/or the step of normal HF excursion by frequency-domain analysis; Wherein, HF is high frequency power value.Normal HF excursion be that Corpus--based Method is determined, the heart hide function normal time span, can be stored in advance in data storage.Particularly, the HF of experimenter and normal HF excursion can be contrasted, if the HF of experimenter is greater than the normal variation scope of HF, then determines that the parasympathetic nervous of experimenter is more excited, thus achieve the detection by quantitative of the mental status to experimenter.
Each step in the present embodiment can be performed by microprocessor, and microprocessor can adopt any one in ARM, DSP, FPGA and MCU.
One of ordinary skill in the art will appreciate that: all or part of step realizing above-mentioned each embodiment of the method can have been come by the hardware that programmed instruction is relevant.Aforesaid program can be stored in a computer read/write memory medium.This program, when performing, performs the step comprising above-mentioned each embodiment of the method; And aforesaid storage medium comprises: ROM, RAM, magnetic disc or CD etc. various can be program code stored medium.
Last it is noted that above each embodiment is only in order to illustrate technical scheme of the present invention, be not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of various embodiments of the present invention technical scheme.
Claims (7)
1. the traditional Chinese medical science heart hides a detection method for function, it is characterized in that, comprising:
The step of the electrocardiosignal of experimenter is gathered in preset time period;
From described electrocardiosignal, extract normal R wave crest point, and obtain the step of RR time interval sequence;
The major axis of Poincare scatterplot and the step of minor axis data is obtained according to described RR time interval sequence; And
Calculate and obtain the ratio of described minor axis with major axis and/or the step of described Poincare scatterplot area.
2. method according to claim 1, is characterized in that, describedly from described electrocardiosignal, extracts normal R wave crest point comprise:
According to mobile point window method and default window length, a point window is carried out to described electrocardiosignal;
Extract the summit of described electrocardiosignal in each window;
Obtain slope characteristics and the R wave crest point slope characteristics on described summit;
When the slope characteristics on described summit is mated with described R wave crest point slope characteristics, using described summit as R wave crest point.
3. method according to claim 1 and 2, is characterized in that, also comprises:
Obtain the step of the ratio of described major axis and/or minor axis and/or minor axis and major axis and/or the normal variation value range of described Poincare scatterplot area.
4. method according to claim 1 and 2, is characterized in that, also comprises:
According to the step of the normal variation value range of area MSEQ and MSEQ of described RR time interval sequence determination multi-scale entropy; And/or
According to the normal variation value range of described MSEQ and MSEQ, the heart obtaining described experimenter hides the step of controlling blood and vessel functional status information.
5. method according to claim 1 and 2, is characterized in that, also comprises:
Described RR time interval sequence is carried out to the step of moving average process;
According to carrying out the RR time interval sequence after moving average, obtain the step of the changes in heart rate curve of described experimenter in preset time period;
According to carrying out the RR time interval sequence after moving average, obtain tested recovery time and/or the step of the tested length of one's sleep of described experimenter;
According to carrying out the RR time interval sequence after moving average, obtain the average heart rate of described experimenter within the described tested length of one's sleep and/or the step of the normal variation value range of average heart rate in the length of one's sleep; And/or
According to carrying out the RR time interval sequence after moving average, obtain that heart rate criteria in the average heart rate of described experimenter in tested recovery time, tested recovery time is poor, the step of the normal variation value range of heart rate standard deviation in the normal variation value range of average heart rate and/or recovery time in recovery time.
6. method according to claim 1 and 2, it is characterized in that, also comprise: according to described RR time interval sequence, obtained the step of the normal variation value range of the normal variation value range of SDANN, SDNNID, rMSSD, PNN50, SDANN of described experimenter, the normal variation value range of SDNNID, the normal variation value range of rMSSD and/or PNN50 by time-domain analysis; Wherein,
The standard deviation of every 5 minutes RR intervals that described SDANN is described experimenter in preset time period;
The standard deviation of every 5 minutes RR intervals that described SDNNID is described experimenter in preset time period;
Described rMSSD is the described experimenter square root that adjacent R R interval difference is mean square in preset time period;
Described PNN50 is the number that described experimenter differs >50 millisecond adjacent R R interval in preset time period.
7. method according to claim 1 and 2, is characterized in that, also comprises: according to described RR time interval sequence, obtains the HF of described experimenter and/or the step of normal HF excursion by frequency-domain analysis; Wherein, described HF is high frequency power value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510846258.8A CN105476630A (en) | 2015-11-26 | 2015-11-26 | Method for detecting traditional Chinese medicine visceral functions |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510846258.8A CN105476630A (en) | 2015-11-26 | 2015-11-26 | Method for detecting traditional Chinese medicine visceral functions |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105476630A true CN105476630A (en) | 2016-04-13 |
Family
ID=55664162
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510846258.8A Pending CN105476630A (en) | 2015-11-26 | 2015-11-26 | Method for detecting traditional Chinese medicine visceral functions |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105476630A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106343992A (en) * | 2016-08-23 | 2017-01-25 | 清华大学 | Heart rate variability analysis method and device and application |
| CN108185994A (en) * | 2018-02-02 | 2018-06-22 | 杭州拜伦医疗科技有限公司 | A kind of balanced capacity detection method and device based on nonlinear technology |
| CN109998527A (en) * | 2019-04-09 | 2019-07-12 | 湖北工业大学 | A kind of heart disease detection method based on multi-scale entropy |
| CN110916649A (en) * | 2019-12-25 | 2020-03-27 | 深圳市博英医疗仪器科技有限公司 | Processing device, processing method and detection device for long-range electrocardiogram scatter diagram |
| CN113397519A (en) * | 2021-08-05 | 2021-09-17 | 季华实验室 | Cardiovascular health state detection device |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020165460A1 (en) * | 2000-06-26 | 2002-11-07 | Starobin Joseph M. | Method and system for evaluating cardiac ischemia with RR-interval data sets and pulse or blood pressure monitoring |
| US20050038351A1 (en) * | 2003-07-23 | 2005-02-17 | Starobin Joseph M. | Method and system for evaluating cardiac ischemia based on heart rate fluctuations |
| CN102271576A (en) * | 2008-12-15 | 2011-12-07 | 艾瑞尔医学有限公司 | Method and system for measuring heart rate variability |
| CN103209637A (en) * | 2010-08-05 | 2013-07-17 | 莱福尔心脏病诊断有限公司 | Method for measuring heart rate variability |
| CN103717125A (en) * | 2011-05-18 | 2014-04-09 | V视股份有限公司 | System and method for determining sleep and sleep stages of a person |
| CN103815896A (en) * | 2014-01-22 | 2014-05-28 | 东软熙康健康科技有限公司 | Mental fatigue monitoring method, device and system and mobile processing terminal |
| CN103976717A (en) * | 2014-04-15 | 2014-08-13 | 德赛电子(惠州)有限公司 | Multidimensional sleeping quality monitoring method and system |
| CN104161509A (en) * | 2014-08-08 | 2014-11-26 | 申岱 | Heart rate variability analyzing method based on amplitude spectrum and instruments |
| CN104305958A (en) * | 2014-10-08 | 2015-01-28 | 西安交通大学 | Photo plethysmograph multi-parameter analyzing method of short autonomic nerve state |
-
2015
- 2015-11-26 CN CN201510846258.8A patent/CN105476630A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020165460A1 (en) * | 2000-06-26 | 2002-11-07 | Starobin Joseph M. | Method and system for evaluating cardiac ischemia with RR-interval data sets and pulse or blood pressure monitoring |
| US20050038351A1 (en) * | 2003-07-23 | 2005-02-17 | Starobin Joseph M. | Method and system for evaluating cardiac ischemia based on heart rate fluctuations |
| CN102271576A (en) * | 2008-12-15 | 2011-12-07 | 艾瑞尔医学有限公司 | Method and system for measuring heart rate variability |
| CN103209637A (en) * | 2010-08-05 | 2013-07-17 | 莱福尔心脏病诊断有限公司 | Method for measuring heart rate variability |
| CN103717125A (en) * | 2011-05-18 | 2014-04-09 | V视股份有限公司 | System and method for determining sleep and sleep stages of a person |
| CN103815896A (en) * | 2014-01-22 | 2014-05-28 | 东软熙康健康科技有限公司 | Mental fatigue monitoring method, device and system and mobile processing terminal |
| CN103976717A (en) * | 2014-04-15 | 2014-08-13 | 德赛电子(惠州)有限公司 | Multidimensional sleeping quality monitoring method and system |
| CN104161509A (en) * | 2014-08-08 | 2014-11-26 | 申岱 | Heart rate variability analyzing method based on amplitude spectrum and instruments |
| CN104305958A (en) * | 2014-10-08 | 2015-01-28 | 西安交通大学 | Photo plethysmograph multi-parameter analyzing method of short autonomic nerve state |
Non-Patent Citations (4)
| Title |
|---|
| YI-LWUN HO 等: "The Prognostic Value of Non-Linear Analysis of Heart Rate Variability in Patients with Congestive Heart Failure—A Pilot Study of Multiscale Entropy", 《PLOS ONE》 * |
| 伍金莉: "心率变异性分析方法的研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 刘丽晓等: "觉醒和睡眠对正常人24小时心率变异性的影响", 《中国心脏起搏与心电生理杂志》 * |
| 霍铖宇: "基于庞加莱散点图和递归图", 《中国博士学位论文全文数据库》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106343992A (en) * | 2016-08-23 | 2017-01-25 | 清华大学 | Heart rate variability analysis method and device and application |
| WO2018036368A1 (en) * | 2016-08-23 | 2018-03-01 | 清华大学 | Method for analyzing heart rate variability, apparatus and use thereof |
| US11154255B2 (en) | 2016-08-23 | 2021-10-26 | Tsinghua University | Method for analyzing heart rate variability, apparatus and use thereof |
| CN108185994A (en) * | 2018-02-02 | 2018-06-22 | 杭州拜伦医疗科技有限公司 | A kind of balanced capacity detection method and device based on nonlinear technology |
| CN108185994B (en) * | 2018-02-02 | 2020-06-09 | 杭州拜伦医疗科技有限公司 | Balance capability detection method and device based on nonlinear technology |
| CN109998527A (en) * | 2019-04-09 | 2019-07-12 | 湖北工业大学 | A kind of heart disease detection method based on multi-scale entropy |
| CN110916649A (en) * | 2019-12-25 | 2020-03-27 | 深圳市博英医疗仪器科技有限公司 | Processing device, processing method and detection device for long-range electrocardiogram scatter diagram |
| CN113397519A (en) * | 2021-08-05 | 2021-09-17 | 季华实验室 | Cardiovascular health state detection device |
| CN113397519B (en) * | 2021-08-05 | 2024-05-28 | 季华实验室 | Cardiovascular health status detection device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105476630A (en) | Method for detecting traditional Chinese medicine visceral functions | |
| US9314177B2 (en) | System and method of detecting abnormal movement of a physical object | |
| US10258250B2 (en) | GPU-based parallel electrocardiogram signal analysis method, computer readable storage medium and device | |
| Liu et al. | Multiscale permutation entropy analysis of electrocardiogram | |
| CN106473700B (en) | The detection system of obstructive sleep apnea syndrome based on HRV frequency-domain analysis | |
| CN106618560B (en) | Method and device for processing brain wave signals | |
| CN111067503A (en) | Sleep staging method based on heart rate variability | |
| CN1736324A (en) | Detection/disposal method and apparatus for obtaining respiratory disturbance index information | |
| US11504052B2 (en) | Modeling method for screening surgical patients | |
| CN110141249B (en) | Non-invasive blood glucose monitoring method, system, equipment and medium based on PPG signal | |
| CN107874750A (en) | Pulse frequency variability and the psychological pressure monitoring method and device of sleep quality fusion | |
| US11154255B2 (en) | Method for analyzing heart rate variability, apparatus and use thereof | |
| CN113303814A (en) | Single-channel ear electroencephalogram automatic sleep staging method based on deep transfer learning | |
| CN112294264A (en) | Sleep staging method based on BCG and blood oxygen saturation rate | |
| CN106295225B (en) | System for detecting sleep apnea syndrome based on mutual information network | |
| CN111317446B (en) | Sleep structure automatic analysis method based on human muscle surface electric signals | |
| CN103405225A (en) | Method, apparatus and device for obtaining pain feeling evaluation indexes | |
| CN204909917U (en) | Intelligent sheet of sleep guardianship and medical diagnosis based on cloth sensor | |
| EP3212075B1 (en) | System for the analysis of the daily heart rhythm autonomic nervous system balance | |
| CN112674779A (en) | Automatic sleep staging method based on ECG (electrocardiogram) signals | |
| Rooijakkers et al. | Low-complexity R-peak detection in ECG signals: A preliminary step towards ambulatory fetal monitoring | |
| CN103892822B (en) | A kind of ECG signal processing method and device | |
| CN113243902A (en) | Feature extraction method based on photoplethysmography | |
| CN110491504B (en) | Method for acquiring medical index data of heart sound signal | |
| Holmer et al. | Detection of ventricular premature beats based on the pressure signals of a hemodialysis machine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160413 |