CN106510737B - A method of Psychic stress state during real-time detection is spoken is estimated with heart rate - Google Patents

A method of Psychic stress state during real-time detection is spoken is estimated with heart rate Download PDF

Info

Publication number
CN106510737B
CN106510737B CN201611122898.5A CN201611122898A CN106510737B CN 106510737 B CN106510737 B CN 106510737B CN 201611122898 A CN201611122898 A CN 201611122898A CN 106510737 B CN106510737 B CN 106510737B
Authority
CN
China
Prior art keywords
heart rate
interphase
stress state
spoken
psychic
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.)
Active
Application number
CN201611122898.5A
Other languages
Chinese (zh)
Other versions
CN106510737A (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.)
Southwest University
Original Assignee
Southwest University
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 Southwest University filed Critical Southwest University
Publication of CN106510737A publication Critical patent/CN106510737A/en
Application granted granted Critical
Publication of CN106510737B publication Critical patent/CN106510737B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
    • 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
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405—Determining heart rate variability
    • 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
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • 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
    • 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]
    • 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
    • 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/7271—Specific aspects of physiological measurement analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • Educational Technology (AREA)
  • Hospice & Palliative Care (AREA)
  • Psychology (AREA)
  • Social Psychology (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

A method of Psychic stress state during real-time detection is spoken, the acquisition including electrocardiosignal, the positioning of R wave, the calculating of heart rate, the calculating of respiratory sinus arrhythmia index, the calculating of Heart Rate Variability and the calculating of heart rate mean value are estimated with heart rate.Under non-strenuous exercise's state, there is short RR interphase mean value and small heart rate variability rate, and (or) caused RR interphase repeatedly continuous elongated in adjacent periods by shallow and slow breathing or continuously shortened, prompt speaker to be in Psychic stress state.Particularly, only repeatedly there is the continuous elongated of the adjacent interphase of RR during speaking or continuously shorten, speaker has usually been prompted to be in weak Psychic stress state；Occur above-mentioned the first two feature simultaneously or all three features occur, then speaker has been prompted to be in strong Psychic stress state；No above three feature occurs, then speaker has been prompted to occur without psychological stress.By detecting small RR interphase mean value, sinus arrhythmia index caused by small heart rate variability rate and shallow and slow breathing, using empirical value method or classical taxonomy device, the method that obtains Psychic stress state during a kind of real-time detection is spoken.

Description

A method of Psychic stress state during real-time detection is spoken is estimated with heart rate
Technical field
The present invention relates to a kind of methods that Psychic stress state during real-time detection is spoken is estimated with heart rate.
Background technique
The long-term existence of psychological stress may cause a variety of physical impairments, for example immunity reduces, and induces cardiovascular disease Deng.In document by statistics and hypothesis testing method, it has already been proven that psychological stress will lead to autonomic nervous system as Lower reaction: sympathetic activation and parasympathetic nerve inhibit.Above-mentioned autonomic nerve stress reaction being presented as in RR interphase, makes RR interval reduction, and the aberration rate of RR interval series is made to become smaller.Although having above-mentioned determining conclusion, psychology is detected at present and is answered Electrocardiosignal feature used in swashing is also in the laboratory research stage, it usually needs baseline characteristic is used as reference to identify psychology Stress situation lacks the method without baseline information for being suitble to real-time psychological stress detection.
Summary of the invention
The problem to be solved in the present invention is: providing a kind of psychological during real-time detection is spoken without being estimated with heart rate for baseline The method of stress situation accurately judges whether have strong or weak psychological stress to occur during speaking.
To solve the above problems, estimating Psychic stress state during real-time detection is spoken with heart rate the present invention provides a kind of Method, which is characterized in that detection method includes the following steps for the psychological stress:
(S1), the electrocardiosignal during acquisition is spoken；
(S2), it detects R wave and RR time interval sequence is calculated；
(S3), respiratory sinus arrhythmia index is calculated, and judges whether to cause RR interphase more by shallow and slow breathing It is secondary continuous elongated in adjacent periods or continuously shorten；
(S4), RR interphase aberration rate index is calculated, and judges whether weak heart rate variability feature occur；
(S5), RR interphase mean value is calculated；
(S6), detect whether Psychic stress state occur, which includes:
Judge whether RR interphase mean value and RR aberration rate index reach the condition of strong psychological stress；
If the condition of not up to strong psychological stress, judges whether respiratory sinus arrhythmia index reaches weak psychological stress Condition；
If being both not up to strong psychological stress condition, also the condition of not up to weak psychological stress, then be judged as no psychological stress.
The beneficial effects of the present invention are: firstly, three hearts rate used in the present invention are estimated all by laboratory data Verifying, the Psychic stress state during can effectively indicating to speak.Secondly, novel RR interphase variation of the present invention Rate index also not by other researchs and technology for detecting Psychic stress state, but inventor pass through laboratory data it was demonstrated that The appearance of strong psychological stress during the index can indicate to speak under conditions of being not necessarily to baseline.Finally, of the present invention Novel respiratory sinus arrhythmia index is yet not by other researchs and technology for detecting Psychic stress state, and the index The weak Psychic stress state in laboratory data can be accurately detected under conditions of being not necessarily to baseline.
Detailed description of the invention
Fig. 1 is flow chart of the invention；
Fig. 2 is the calculating of R wave positioning and RR interval series in embodiment；
Fig. 3 is the RR interval series schematic diagram of weak heart rate variability feature in embodiment；
Fig. 4 is the corresponding adjacent R R interphase difference sequence schematic diagram of RR interphase of weak heart rate variability feature in embodiment；
Fig. 5 is the RR interval series schematic diagram of sinus arrhythmia feature caused by breathing shallow and slow in embodiment；
Fig. 6 is the corresponding adjacent R R interphase difference sequence of sinus arrhythmia feature caused by breathing shallow and slow in embodiment List intention.
Specific embodiment
The present invention is described further in the following with reference to the drawings and specific embodiments:
As shown in Figure 1, the method for the invention that Psychic stress state during real-time detection is spoken is estimated with heart rate, including the heart The acquisition of electric signal, the positioning of R wave, the calculating of heart rate, the calculating of heart rate mean value, the calculating and respiratory sinus property of Heart Rate Variability The calculating of cardiac arrhythmia index.
Wherein heart rate is calculated from the positioning of the R wave of electrocardiosignal and the interval RR, and electrocardiosignal is directly to pass through life Reason signal collecting device obtains.The positioning of R wave is to detect R wave crest point with threshold method.The case where losing R wave is acquired when there are data, Extremely long RR interphase, abnormal short RR interphase and adjacent R R interphase abnormal big caused by being lost with threshold method detection R wave Difference.If R wave loses the detection for seriously stopping psychological stress and data reporting acquisition is not normal；If R wave only has a small number of loss, RR time interval sequence is calculated according to R wave vertex position and sample frequency, as shown in Figure 2.
The mean value of heart rate is added using RR gap lengths all in the data buffering time again divided by the total interval RR number Mode is calculated, and exception RR interphase is first removed before calculating.In the analysis of experimental data carried out, the data buffering time is set It is 120 seconds.Normal cardiac rate range is widely considered to be 60-100 beats/min.It is tachycardia when heart rate is greater than 100 beats/min. One of the characteristics of psychological stress is exactly to cause tachycardia.Therefore, can be with 100 beats/min, i.e. the interval 0.6 second/RR, as strong Psychological stress empirical value Tm.When RR gap length mean value is less than the empirical value, in conjunction with heart rate variability rate index Judge whether strong psychological stress jointly.
The empirical value T of weak heart rate variability is arranged in the calculating of Heart Rate Variability firstw.It is 400 in data sampling rate When point/second, the empirical value of weak heart rate variability is Tw=0.0075 second.When adjacent R R interphase difference is less than the empirical value, sentence Break as weak heart rate variability.And judging whether strong psychological stress also needs to calculate as follows:
(1) within the data buffering time, between the adjacent R R that interception was approximately equal in a slow breathing cycle time without overlapping Phase difference data.The respiratory rate of normal adult is about 12 beats/min to 22 beats/min, i.e., the cycle length of one time slow breathing is about 5 Second, it can be in this, as the data length empirical value T for investigating weak heart rate variability frequency of occurrenceb。
It (2) is T to intercepted lengthbAdjacent R R interphase difference data, if RR interphase difference be less than TwThe number that situation occurs is greater than TsIt is secondary, then it is assumed that heart rate variability feature occurs primary.Wherein, TsIt is the empirical value for judging weak heart rate variability feature and whether occurring, To existing laboratory data, which is taken as 7.
(3) it to the data cutout each time in the entire data buffering time, counts total time that weak heart rate variability feature occurs Number.When the total degree that this feature occurs is greater than empirical value TnWhen, in conjunction with the tachycardia embodied with RR interphase mean value Judge that strong psychological stress occurs.
Steps are as follows for the calculating of respiratory sinus arrhythmia index:
(1) it is similar to Heart Rate Variability to calculate, non-overlapping T is intercepted within the data buffering timebThe adjacent R R of length Whether interphase difference sequence corresponds to shallow breathing with threshold decision adjacent R R interphase difference sequence.The experience threshold that existing experimental data uses Value is TrCorrespond to shallow breathing within < 0.03 second.
(2) to the adjacent R R interphase difference sequence for meeting shallow breathing condition, the component that value is positive in the sequence is calculated separately Number, the number for the component that the number and value for the component that value is negative are zero, and by value be zero component number close And into positive/negative component number biggish one.
(3) whether it is greater than threshold value T with both positive/negative component numbers of RR interphase difference sequence biggish onelTo determine whether There is shallow and slow breathing.To existing experimental data, TlEmpirical value be 6.
(4) to the data cutout each time in the entire data buffering time, Dou Xingxin caused by the shallow and slow breathing of statistics Restrain the total degree that uneven feature occurs.When the total degree that this feature occurs is greater than empirical value ToWhen, that is, it can determine whether weak psychology Stress occur.
It not will lead to more inspections of weak heart rate variability feature with abnormal small RR interphase greatly in view of exception and shallowly exhale slowly More inspections of feature are inhaled, in order to retain the time sequencing of adjacent R R interval series, are calculating Heart Rate Variability and sinus rhythm When uneven index, removal processing is not done to the abnormal RR interphase occurred once in a while.Fig. 3 shows that is in a strong Psychic stress state Speaker 120 seconds length RR interval series, corresponding adjacent R R interphase difference sequence (Fig. 4 short line) and weak heart rate variability Feature (Fig. 4 long lines section).Know that wherein weak heart rate variability feature occurs 10 times.Fig. 5 shows that is in a weak psychology 120 seconds length RR interval series of the speaker of stress situation, corresponding adjacent R R interphase difference sequence (Fig. 6 short line) and are exhaled The shallow and slow respiratory characteristic (Fig. 6 long lines section) that absorption sinus arrhythmia index reflects, it is known that shallow and slow breathing Feature occurs 7 times.
Above-mentioned threshold value Tm, TnAnd ToJudge whether Psychic stress state occur suitable for empirical value method.To being computed RR interphase mean value, RR interphase aberration rate index value and the respiratory sinus arrhythmia index value arrived, can also first use given data Any one classical taxonomy device of training, such as support vector machine classifier, then detect new data with trained listening group In whether there is Psychic stress state.
The above description of this invention is illustrative and not restrictive, those skilled in the art understand that wanting in right Ask it can be carried out within the spirit and scope of restriction it is many modification, variation or it is equivalent, but they fall within it is of the invention In protection scope.

Claims (2)

1. a kind of method for estimating Psychic stress state during real-time detection is spoken with heart rate, which is characterized in that psychological stress inspection Survey method the following steps are included:
(S1), the electrocardiosignal during acquisition is spoken；
(S2), it detects R wave and RR time interval sequence is calculated；
(S3), respiratory sinus arrhythmia index is calculated, and judges whether to cause RR interphase repeatedly to exist by shallow and slow breathing It is continuous elongated or continuously shorten in adjacent periods；
(S4), RR interphase aberration rate index is calculated, and judges whether weak heart rate variability feature occur；
(S5), RR interphase mean value is calculated；
(S6), detect whether Psychic stress state occur, which includes:
Judge whether RR interphase mean value and weak heart rate variability characteristic index reach the condition of strong psychological stress；
If the condition of not up to strong psychological stress, judges whether respiratory sinus arrhythmia index reaches the item of weak psychological stress Part；
If being both not up to strong psychological stress condition, also the condition of not up to weak psychological stress, then be judged as no psychological stress.
2. estimating the method for Psychic stress state during real-time detection is spoken with heart rate as described in claim 1, feature exists In the primary data buffer time of real-time detection process is generally no less than 120 seconds.
CN201611122898.5A 2015-12-16 2016-12-08 A method of Psychic stress state during real-time detection is spoken is estimated with heart rate Active CN106510737B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2015109455793 2015-12-16
CN201510945579 2015-12-16

Publications (2)

Publication Number Publication Date
CN106510737A CN106510737A (en) 2017-03-22
CN106510737B true CN106510737B (en) 2019-07-26

Family

ID=58341864

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611122898.5A Active CN106510737B (en) 2015-12-16 2016-12-08 A method of Psychic stress state during real-time detection is spoken is estimated with heart rate

Country Status (1)

Country Link
CN (1) CN106510737B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109044335B (en) * 2018-07-17 2020-11-10 西安交通大学 Heart function evaluation method based on instantaneous sound stimulation
CN108937916A (en) * 2018-08-03 2018-12-07 西南大学 A kind of electrocardiograph signal detection method, device and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014061079A (en) * 2012-09-20 2014-04-10 Crosswell:Kk Autonomic nervous function evaluation apparatus and program
CN103815896A (en) * 2014-01-22 2014-05-28 东软熙康健康科技有限公司 Mental fatigue monitoring method, device and system and mobile processing terminal
CN104173064A (en) * 2014-09-04 2014-12-03 西双版纳生物医学研究院 Heart rate variability analysis based lie detection method and lie detection device
WO2015060267A1 (en) * 2013-10-21 2015-04-30 テイ・エス テック株式会社 Alertness device, seat, and method for determining alertness
CN104755020A (en) * 2012-10-23 2015-07-01 皇家飞利浦有限公司 Stress-measuring system
WO2015121503A1 (en) * 2014-02-17 2015-08-20 Vasile Zoicas Method and system for monitoring the autonomic nervous system of a subject

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8155733B2 (en) * 2009-05-15 2012-04-10 Nanjing University Lie detection method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014061079A (en) * 2012-09-20 2014-04-10 Crosswell:Kk Autonomic nervous function evaluation apparatus and program
CN104755020A (en) * 2012-10-23 2015-07-01 皇家飞利浦有限公司 Stress-measuring system
WO2015060267A1 (en) * 2013-10-21 2015-04-30 テイ・エス テック株式会社 Alertness device, seat, and method for determining alertness
CN103815896A (en) * 2014-01-22 2014-05-28 东软熙康健康科技有限公司 Mental fatigue monitoring method, device and system and mobile processing terminal
WO2015121503A1 (en) * 2014-02-17 2015-08-20 Vasile Zoicas Method and system for monitoring the autonomic nervous system of a subject
CN104173064A (en) * 2014-09-04 2014-12-03 西双版纳生物医学研究院 Heart rate variability analysis based lie detection method and lie detection device

Also Published As

Publication number Publication date
CN106510737A (en) 2017-03-22

Similar Documents

Publication Publication Date Title
JP4743534B2 (en) Heart rate detector
RU2015128271A (en) SYSTEM AND METHOD FOR DETECTING SIGNIFICANT ARRHYTHM EVENTS BY PHOTOPLETISMOGRAM (PPG) AND ACCELEROMETER
US9254095B2 (en) Electrocardiogram signal detection
CN104905785B (en) ECG signal detection method and system
WO2015025187A8 (en) Method and device for improving prediction and detection of change in a physiological condition
WO2015101698A3 (en) A method and system for measuring stress and recovery
CN106510737B (en) A method of Psychic stress state during real-time detection is spoken is estimated with heart rate
JP2011528239A5 (en)
JP2017042269A (en) Blinking detection device
RU2012101554A (en) ECG CONTROL WITH REDUCED FALSE ALARMS IN CONNECTION WITH ASYSTOLIA
WO2008007236A3 (en) Atrial fibrillation detection
JP2007512043A5 (en)
JP2012502670A5 (en)
WO2010123542A3 (en) Methods and systems for detecting epileptic events using nonlinear analysis parameters
CN104939810B (en) A kind of method and device controlled the emotion
CN106037671A (en) Method and system for apnea event detection based on BCG signal
WO2010068588A3 (en) Integrated instrumentation system and method for assessing feeding readiness and competence in preterm infants
CN104720794A (en) Method for detecting electrocardiogram wave form R points through superposing
KR101245445B1 (en) R-peak detection method in ECG Signal and apparatus thereof
WO2012087515A4 (en) Rate initialization and overdrive pacing for capture threshold testing
Kang et al. Heart sound segmentation toward automated heart murmur classification in pediatric patents
EP1774907A1 (en) Method and device for measuring heart rate
CN110507299A (en) Heart rate signal detection device and method
Reyes et al. Analysis of discontinuous adventitious lung sounds by Hilbert-Huang spectrum
WO2011061606A3 (en) Methods and systems for atrial fibrillation and st segment deviation detection

Legal Events

Date Code Title Description
C06 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