CN109044335A - A kind of cardiac function evaluation method based on distance speech stimulation - Google Patents
A kind of cardiac function evaluation method based on distance speech stimulation Download PDFInfo
- Publication number
- CN109044335A CN109044335A CN201810783946.8A CN201810783946A CN109044335A CN 109044335 A CN109044335 A CN 109044335A CN 201810783946 A CN201810783946 A CN 201810783946A CN 109044335 A CN109044335 A CN 109044335A
- Authority
- CN
- China
- Prior art keywords
- cardiac function
- stimulation
- interphase
- time point
- testee
- 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
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]
-
- 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
-
- 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/7253—Details of waveform analysis characterised by using transforms
Abstract
The invention discloses a kind of cardiac function evaluation methods based on distance speech stimulation, pass through the multiple auditory stimulation experiment of progress and records the electrocardiosignal for stimulating testee in preceding and stimulating course, using testee by auditory stimulation after, it will appear the temporary decline of heart rate in a short time, that is the temporary increase of RR interphase is pre-processed to obtain RR interphase RR (n) by shooting to the electrocardiosignal of acquisition；Then by method for normalizing, the average RR interphase variation after stimulating in a period of time is calculated, obtains testee interior RR interphase average increase Δ RRI after stimulationin, to obtain heart and brain coupled relation according to the cardiac function of tester, this method is convenient and easy, it is widely applicable, it is a kind of useful mode for evaluating heart responding ability, and it is widely applicable.
Description
Technical field
The invention belongs to processing of biomedical signals technical fields, it is intended to propose a kind of heart based on distance speech stimulation
Function evaluation methods.
Background technique
Heart is organ particularly important in human body, it pushes blood to press in the blood vessel by Rythmic contractions characteristic and diastole
Ceaselessly circulated according to certain direction, thus ensure that organismic internal environment it is relative constant and metabolic it is normal into
Row, malfunction will lead to atrial fibrillation, ventricular arrhythmia, heart failure pernicious cardiovascular event, and then cause patient its
The anoxic of his internal organs is impaired or even dead." global non-communicable diseases status report in 2010 " display, had 1730 in 2008
Ten thousand people die of cardiovascular disease, account for the 30% of global dead sum.Expect the year two thousand thirty, death toll caused by cardiovascular disease
It will be increased to 23,300,000 people.
Research confirms that heart is in addition to the basilic rhythm bounce according to sinoatrial node, while also by body fluid and nerve modulation
Effect, exactly have ignored heart to adapt to body and external environment variation, the above method and correctly receive nerve modulation and adaptation
The ability of environment.
And in numerous mechanism for adjusting heart, brain is undoubtedly most important and is also most direct one.
In fact, between brain and heart by numerous nerves, hormone and biophysics mechanism carry out in real time exchange with
Ensure the homeostasis of heart, while complicated angiocarpy and system are made to various external stimulus and responded.Based on close coupling
Conjunction relationship, brain (or heart) structure and/or function lesion can cause heart (or brain) dysmotility.It is handed over as typical brain-heart
Logical circulation road has furtherd investigate the regulation of cardiomotility using brain stem oblongata as the basic autonomic nerves system for adjusting maincenter.
However, cerebral cortex can regulate and control cardiovascular activity by autonomic nerves system, the pass between high-level center and heart
System's (brain-heart coupling) receives more and more attention recently.Research brain-heart coupling facilitates we have appreciated that body most important two
Physiology connection between big organ, and preferably recognize cardiac-related diseases mechanism such as cardiac arrhythmia, hypertension, heart failure and its interior
It is contacting, is more effectively preventing and treating strategy to make.
Studies have shown that environmental stimuli play in the pathogenesis of hypertension and other cardiovascular diseases it is extremely important
Role.Stress situation will lead to the change dramatically of autonomic nervous function, causes cardiovascular danger event and even dies suddenly, especially right
There is the people at highest risk of cardiovascular disease history and history of operation.Research finds emotional stress and stress is coronary heart disease (CHD) sudden death
One big induced factor.It is very wide to Cardiovascular diseases affect covering surface about psychological stress, including it is based on particular disease, it is based on group
The major issue (such as earthquake) of body and the research about Chronic Pressure in epidemiology field.It is depressed to cause angiocarpy
Patient's acute myocardial infarction, and angry rear this level of significance can be higher.There is research to confirm that the indignation of height can improve implanted heart rate
The probability of happening of defibrillator patient's ventricular arrhythmia.In addition to the research for the patient for having heart disease, also there is minority to be directed to
Without coronary heart disease, the research of the normal patient of heart (such as congenital ventricular fibrillation patient) finds psychological stress and sudden cardiac arrest
Sudden death have important relationship.
Scaring stimulation can lead to arrhythmia cordis and sudden death by increasing Sympathetic Nerve impulsion.Under Psychic stress state, heart rate
The corresponding brain area activity of the not normal unstable and autonomic nerve along with myocardium state is extremely related.This shows under stress situation,
There is abnormal autonomic nerve to get excited and is passed to heart.And this exception has occurred often in cardiovascular disease and hypertension history or family
In the crowd of race's history, in these individuals, abnormal ventricular wall motion and myocardial function degeneration can be put under stress situation
Greatly, lead to the disturbance of heart afferent impulse.The incoming feedback of heart can further influence autonomic nerve and make to cardiovascular control
With.The incoming feedback of i.e. abnormal heart, can increase the risk of cardiac arrhythmia.And for most of absent cardiovasculars and hypertension history
And the healthy population of family history, although the scaring stimulation of general level and psychological pressure can cause heart and autonomic nerve certain
Response in degree, but this response can't induce serious cardiovascular event and fatal arrhythmia.Therefore, heart function
The evaluation and test of energy, establishes cardiac function and heart and brain coupled relation, for preventing and reducing pernicious cardiovascular event, examines related drugs
Curative effect has important clinical value, is the research hotspot of field of biomedicine all the time.Existing cardiac function is commented
Valence method is based primarily upon the bioelectrical signals and video signal of tranquillization state, passes through the place to quiescent condition biomedicine signals
Reason, obtains heart rate, often a series of cardiac functions such as the amount of fighting, cardiac output, ventricle thickness, ejection fraction.But the existing heart
Still there are limitations for dirty functional parameter, have only evaluated and tested the intrinsic rhythm and pace of moving things bounce of heart and the ability of pump blood, have not studied
To the coupled relation between heart and brain.
Summary of the invention
It is existing to overcome the purpose of the present invention is to provide a kind of cardiac function evaluation method based on distance speech stimulation
The deficiency of technology.
In order to achieve the above objectives, the present invention adopts the following technical scheme:
A kind of cardiac function evaluation method based on distance speech stimulation, comprising the following steps:
Step 1) carries out multiple auditory stimulation experiment and records the electrocardio letter for stimulating testee in preceding and stimulating course
Number；
Step 2) is pre-processed to obtain RR interphase RR (n) by shooting to the electrocardiosignal of acquisition；
Step 3) carries out the calculating of RR interphase response sequence to the RR interphase RR (n) that step 2) obtains, and obtains normalized RR
Interval series RR ' (t ', i)；
The normalized RR interval series RR ' (t ', i) that step 4), basis obtain calculates testee interior RR after stimulation
Interphase average increase Δ RRIin, calculation method is as follows:
ΔRRIin=(RR'(0.5)+RR'(1))/2；
Step 5), by Δ RRIinAs the supplement index of basic cardiac function, cardiac function evaluation is carried out, if Δ
RRIin≤ 0, then it is assumed that Neuroprotective Mechanisms are not perfect.
Further, to collected original electro-cardiologic signals ecg (t), first with wavelet decomposition remove baseline drift and
High-frequency noise, the electrocardiosignal x (t) after being denoised；Then the method for utilizing Wavelet Modulus Maxima, the small echo of selection are to support
Compactly support and there is single order vanishing moment db2 small echo, R wave crest is detected, then inspect modification, acquisition one to testing result
The set R (N) of the continuous R peak time point of consecutive；To the set R (N) of R peak time point, the latter peak time point R (n is utilized
+ 1) the time point R (n) of previous wave crest is subtracted, RR interphase RR (n) by shooting is obtained, wherein (1≤n < N).
Further, it to collected original electro-cardiologic signals ecg (t), selects db2 small echo to carry out 9 layers of decomposition to signal, goes
Electrocardiosignal x (t) except the baseline drift of the last layer and the high-frequency noise of three first layers, after being denoised.
Further, R wave crest is detected using the method for Wavelet Modulus Maxima, using support compactly support and has one
The equivalence filter coefficient of the db2 small echo of rank vanishing moment, db2 small echo is as follows:
H1=0.3750, h2=0.1250, h3=0.0000
G1=0.5798.g2=0.0869, g3=0.0061
hk=h1-k,gk=-g1-k
If k > 3, hk=gk=0
Wavelet decomposition is made to the collected electrocardiosignal Mallat algorithm of institute.
Further, the time point t (i) of sonic stimulation is recorded, wherein i=1,2 ..., 20, it will be every obtained in step 2)
The RR interval series of a stimulation time point t (i) nearby carry out adopting again from the 2Hz in the times s such as t (i) time point front and back respectively
Sample obtains the unified RR interval series RR (t, i) of sample rate, and when the stimulation appearance point for newly obtaining RR interval series is defined as
Between put 0, data start time point is-s, and end time point is s, i.e. RR (t ', i).
Further, it specifically, carrying out the resampling from t (i) -2.5s to the 2Hz of t (i)+2.5s respectively, is sampled
The unified RR interval series RR (t, i) of rate, and the stimulation appearance point for newly obtaining RR interval series is defined as time point 0, data
Start time point is -2.5, and end time point is 2.5, i.e. RR (t ', i).
Further, in step 5), specifically, carrying out cardiac function diagnosis first；
If testee's items cardiac function base values is normal, and index Δ RRIin> 0, then heart regulating power is strong；
If testee's items cardiac function base values is normal, and index Δ RRIin≤ 0, then heart regulating power is weak；
If testee's cardiac function base values has abnormal index, and index Δ RRIin≤ 0, then heart regulating power
Difference.
Further, obtained RR interval series RR (t ', i) is normalized, obtains normalized RR interval series
RR ' (t ', i):
Wherein, RR_baseline is average value of the RR interval series in -2.5 < t < 0 section.
Further, RR'(t', i) is averaging along stimulation number, the RR interphase sonic stimulation for obtaining testee is rung
Answer sequence RR'(t').
Compared with prior art, the invention has the following beneficial technical effects:
A kind of cardiac function evaluation method based on distance speech stimulation of the present invention, by carrying out multiple auditory stimulation experiment
And record stimulation before and stimulating course in testee electrocardiosignal, using testee by auditory stimulation after, in short-term
The interior temporary decline that will appear heart rate, the i.e. temporary increase of RR interphase pre-process to the electrocardiosignal of acquisition
To RR interphase RR (n) by shooting；Then by method for normalizing, the average RR interphase variation after stimulating in a period of time is calculated, is obtained
To testee after stimulation in RR interphase average increase Δ RRIin, to obtain heart and brain coupling according to the cardiac function of tester
Conjunction relationship, this method is convenient and easy, widely applicable, is a kind of useful mode for evaluating heart responding ability, widely applicable.
Detailed description of the invention
Fig. 1 is a kind of system block diagram of cardiac function evaluating method based on distance speech stimulation；
Fig. 2 is auditory stimulation experiment flow schematic diagram；
Fig. 3 is R wave crest point detection algorithm flow chart；
Fig. 4 is R wave crest point testing result figure；
Fig. 5 is the RR interphase response sequence calculated result figure of family history of hypertension and Healthy People；
Fig. 6 is the Δ RRI of family history of hypertension and Healthy PeopleinThe statistical results chart of index.
Specific embodiment
The invention will be described in further detail with reference to the accompanying drawing:
As shown in Figure 1, a kind of cardiac function evaluation method based on distance speech stimulation, comprising the following steps:
Step 1) carries out multiple auditory stimulation experiment and records the electrocardio letter for stimulating testee in preceding and stimulating course
Number；
Step 2) is pre-processed to obtain RR interphase RR (n) by shooting to the electrocardiosignal of acquisition；
Step 3) carries out the calculating of RR interphase response sequence to each mode of heart rate that step 2) obtains, and obtains normalized RR
Interval series RR ' (t ', i)；
Step 4) seeks its average increase Δ RRI to the normalized RR interval series RR ' (t ', i) that step 3) obtainsin,
Calculation method is as follows:
ΔRRIin=(RRt(0.5)+RRt(1))/2；
The implementation steps of modules in attached drawing are specifically described in the following with reference to the drawings and specific embodiments:
1, auditory stimulation process described in step (1) of the present invention, as shown in Fig. 2, including following procedure:
Before on-test, test process is introduced to testee, testee enters temperature and illumination is suitable for, and sound insulation is good
Laboratory after seated resting 30 minutes；
Connect electrode apparatus, the electrocardiosignal ecg (t) of start recording testee in a particular state；
After laboratory keeps peace and quiet in 5 minutes, start that 20 100 milliseconds 100 decibels of white noise Sound stimulats are presented, between stimulation
Every being randomly dispersed between 40 seconds to 60 seconds；Record the electrocardiosignal of testee in stimulating course, record end；
2, ECG signal processing:
To collected original electro-cardiologic signals ecg (t), baseline drift and high-frequency noise are removed first with wavelet decomposition,
Electrocardiosignal x (t) after being denoised；Then the method for utilizing Wavelet Modulus Maxima, the small echo of selection be support compactly support and
With single order vanishing moment db2 small echo, R wave crest is detected, then testing result is carried out to inspect modification, is obtained a succession of continuous
The set R (N) of R peak time point；To the set R (N) of R peak time point, subtracted using the latter peak time point R (n+1)
The time point R (n) of previous wave crest obtains RR interphase RR (n) by shooting, wherein (1≤n < N)；It is obtained by shooting using HR=60/RR
Heart rate HR (n)；
To collected original electro-cardiologic signals ecg (t), selects db2 small echo to carry out 9 layers of decomposition to signal, remove last
The baseline drift of layer and the high-frequency noise of three first layers, the electrocardiosignal x (t) after being denoised；
R wave crest is detected using the method for Wavelet Modulus Maxima, using support compactly support and there is single order vanishing moment
Db2 small echo, the equivalence filter coefficient of db2 small echo is as follows:
H1=0.3750, h2=0.1250, h3=0.0000
G1=0.5798.g2=0.0869, g3=0.0061
hk=h1-k,gk=-g1-k
If k > 3, hk=gk=0
Wavelet decomposition is made to the collected electrocardiosignal Mallat algorithm of institute；
1) selection of characteristic dimension:
Signal decomposition is the ingredient of different frequency range by wavelet transformation, and high fdrequency component and noise are mainly fallen on small scale, low
Frequency component and noise are mainly fallen in large scale.For different people, the frequency spectrum of QRS wave is slightly different in ECG signal, but energy
It is concentrated mainly on scale 23With 24On, and 23Upper energy is maximum.The more QRS wave of high fdrequency component, 22Upper energy is maximum, and low
The more QRS wave of frequency component, 24Upper energy is maximum.Therefore, 2 are selected in this paper1-24Four scales.Using db2 small echo, the heart
Electric signal carries out four layers of decomposition, obtains wavelet transformation ingredientWith smooth signal
2) determination of R wave modulus maximum column:
R wave can generate a pair of of wavelet modulus maxima, the i.e. negative minimum of positive maximum-on each characteristic dimension
It is right.High-frequency noise can be in scale 21With 22Upper generation modulus maximum, and the high P wave of low frequency or high T wave can be in scales 23With 24Upper life
At modulus maximum, therefore, the modulus maximum column on this 4 scales is detected, are examined with reducing noise and high P wave and high T wave to R wave
The influence of survey.
The corresponding modulus maximum point of R wave is determined by following steps:
(1) from scale 24Wavelet transform result in find out greater than threshold epsilon4Modulus maximum point, obtain these point position
Set
(2) from scale 23Wavelet transform result in,It finds out in neighborhood greater than threshold epsilon3And withLocate modulus maximum
Point with the modulus maximum point of symbol, be selected as by " neighborhood " hereinEach 10 points in left and right, its position is set toIfNear
There are several modulus maximum points, then selects maximum one.If the width of amplitude several modulus maximum points other less than 1.2 times of this point
Value, then choosing nearPoint.IfIn neighborhood not withLocate modulus maximum point with the modulus maximum point of symbol, then enables
It is 0；Obtain the set of all these positions
(3) step (2) are repeated, finds scale 22、21On modulus maximum point position setWith
In this algorithm, the different amplitude threshold { ε of different scale uses4, ε3, ε2, ε1According to being most recently detected
Wavelet modulus maxima refreshes formula come what is refreshed are as follows:
IfThen,
Otherwise
And
WhereinRepresent the wavelet modulus maxima detected.This method can guarantee, the QRS of detection
It will not influence later amplitude features when wave amplitude increases suddenly.
Use it is this from large scale to small dimension search wavelet modulus maxima point method can remove it is small
The modulus maximum point generated on scale by noise accurately detects R wave crest point, while can save operation time.
3) calculating of singular point Singularity Degree:
It enablesAssuming that a is the upper limit of Lipschitz index.Enable a ≈ log2 aj+1(nk)-log2
aj(nk).Pass through four layers of wavelet decomposition, available a1、a2And a3.There must be a at R wave crest point1> 0, and a under normal conditions2
> 0, and even if a2When < 0, a1+a2Inherently it is greater than 0.For most of R waves, usually there is a3< 0 and a1+a2+a3≤ 0,
And acutely for high-frequency noise and interference, a1≤ 0, a2≤ 0, a3≤ 0, and a1+a2+a3≤0.Therefore, from a1+a2+a3Value cannot
Resolution R wave, high-frequency noise and interference, and a1+a2There is good resolving effect.So having selected a when detecting R wave1、a2, and
Enable a'=(a1+a2)/2.If a' > 0, corresponding modulus maximum point is corresponding to R crest value point；If a' reduces become suddenly
Negative value, then corresponding modulus maximum point must be that noise or interference are corresponding, should delete corresponding modulus maximum column.
4) the isolated and extra modulus maximum column of removal:
Band overlapping of the frequency band of motion artifacts and myoelectricity noise usually with QRS wave.It therefore, should be from the obtained mould in front
In maximum column set, rejects and arranged by the modulus maximum that artefact or myoelectricity noise introduce.
(1) isolated modulus maximum column are deleted.
R wave both corresponds to a pair of of modulus maximum column, the i.e. negative minimum pair of positive maximum-on each characteristic dimension.This two
The spacing of a modulus maximum point is in scale 21On it is smaller than the width of R wave.IfFor scale 21OnA positive maximum
Point,For on same scaleNegative minimum point, ifWithSpacing is greater than given threshold value, if
Spacing threshold is 120ms, then claimsTo isolate modulus maximum point, should be deleted.
(2) extra modulus maximum column are deleted.
One R wave only generates a pair of of modulus maximum point on each scale.But some R waves with noise, in a pair of of anode
In the neighborhood (threshold value 120ms) of the negative minimum pair of big value-, two pairs or more of modulus maximum column can be generated, and are wherein only had
A pair is useful.Following criterion can be used to delete in extra modulus maximum column:
Here the scale 2 for selecting the energy of QRS wave mainly to concentrate3On modulus maximum differentiate.If two negative minimums
Point is respectively Min1 and Min2, and amplitude is respectively A1 and A2, and they are respectively L1 at a distance from corresponding positive maximum point
And L2.
Criterion 1: ifThen Min2 is extra modulus minimum point；
Criterion 2: ifThen Min1 is extra modulus minimum point；
Criterion 3: otherwise, if Min1, Min2 in the ipsilateral of the positive maximum, then remote from the positive maximum is extra mould
Maximum point；If Min1, Min2 are in the two sides of the positive maximum point, then that point after the positive maximum point is extra point.
5) R wave crest point detects:
Front algorithm has eliminated noise and interference and the corresponding negative minimum pair of positive maximum-of P wave T wave, and
Isolated and extra modulus maximum column are deleted, scale 2 is obtained1The corresponding positive maximum-of upper ECG signal R wave crest point
Negative minimum pair；The position for detecting these positive maximum-negative maximum pair zero crossings, has just obtained the position of R wave crest point；
As shown in Figure 4；
In order to improve verification and measurement ratio, two strategies have also been used:
1) refractory period:
The heart rate of common people will not generate another secondary heart less than 300 beats/min within a period of time after a heartbeat
It jumps, that is, has one section of refractory period；Therefore, after detecting a R wave, the extreme value in 200ms thereafter is all ignored, it can
To avoid erroneous detection caused by due to noise jamming.
2) reverse search:
In the case that arrhythmia cordis or other, R wave amplitude and frequency can become smaller suddenly, and the amplitude of modulus maximum point is caused to reach
Less than threshold value, lead to missing inspection；It in our algorithm, is first averaged, is obtained recently to RR interphase detected in first 30 seconds
The average heart cycle T of a period of time, if the RR interphase of this detection is greater than 1.5T, in this interphase intrinsic time theory 23Upper use
0.5ε3Detect modulus maximum.If the negative minimum of a pair of of positive maximum-in this section is less than 140ms to the interval between point,
Think there is missing inspection, detects the zero cross point between them, and corrected with 3 points of time shift, the R wave examined again.Using this side
Method can reduce missing inspection in most cases.
Detection algorithm process as shown in figure 3,
It to above method testing result, is finally inspected and is modified again, after guarantee is errorless, obtain a succession of continuous R wave
The set R (N) of peak time point.
3, RR interphase response sequences calculate:
The time point t (i) of sonic stimulation, wherein i=1,2 ..., 20 are recorded, it will be obtained in step 2) when each stimulation
Between the RR interval series of point t (i) nearby carry out the resampling from t (i) -2.5s to the 2Hz of t (i)+2.5s respectively, sampled
The unified RR interval series RR (t, i) of rate, and the stimulation appearance point for newly obtaining RR interval series is defined as time point 0, and will
The stimulation appearance point for newly obtaining RR interval series is defined as time point 0, and data start time point is-s, and end time point is s, i.e.,
RR(t',i)；S value range is 1-5s.
In order to exclude excessive interference of the single test to whole result, obtained RR interval series RR (t ', i) is returned
One changes, and obtains normalized RR interval series RR ' (t ', i), the method is as follows:
Wherein, RR_baseline is average value of the RR interval series in -2.5 < t < 0 section.
RR'(t', i) is averaging along stimulation number, obtains the RR interphase sonic stimulation response sequence RR' of testee
(t')。
Wherein RR interphase response sequence such as Fig. 5 of normal health subject is shown in solid, the RR interphase response of hypertension subject
Sequence is as shown in Fig. 5 dotted line；
4, index Δ RRIinIt calculates:
Calculate testee interior RR interphase average increase Δ RRI after stimulationin, calculation method is as follows:
ΔRRIin=(RR'(0.5)+RR'(1))/2
According to the calculation method, statistical result of the index in family history of hypertension testee and healthy testee
As shown in fig. 6, the index is in two groups of testees, there were significant differences, the Δ RRI of family history of hypertension testeeinIndex is aobvious
It writes and weakens, show that parasympathetic protection mechanism of the heart in auditory stimulation response process is significantly reduced compared to normal testee；
I.e. by Δ RRIinAs the supplement index of basic cardiac function, if Δ RRIin≤ 0, then it is assumed that neuroprotection
Mechanism is not perfect.
Specifically, in conjunction with Δ RRIinCardiac function is evaluated, progress cardiac function diagnosis first,
If testee's items base values is normal, index Δ RRIin> 0, then heart regulating power is strong, i.e., heart is basic
Function is good, and the neuroprotection response mechanism that heart faces burst stimulation is perfect, and cardiovascular event risk is lower；
If testee's items base values is normal, index Δ RRIin≤ 0, then heart regulating power is weak, i.e., heart is basic
Function is good, but heart face burst stimulation neuroprotection response mechanism it is not perfect, cardiovascular event risk is higher；
If testee has abnormal base values, and index Δ RRIin≤ 0, then heart regulating power is poor, i.e. the heart
Dirty basic function has been damaged, and heart faces not perfect, the cardiovascular event risk of Neuroprotective Mechanisms response of burst stimulation
It is very high.
This has especially selected the subject data of Healthy subjects and family history of hypertension in illustrating.With high blood
The subject of pressure family history is widely believed that the bigger risk with hypertension and other cardiovascular diseases, but in this hair
In the experiment that bright people carries out, if with traditional evaluating method, the subject quiescent condition of Healthy subjects and family history of hypertension
Under heart basic function index not there is significant difference.This explanation is in existing cardiac function evaluating method, hypertension
Family history subject and the cardiac function of Healthy subjects are almost the same, it is clear that and this does not simultaneously meet the known fact, therefore, existing index
The health of heart degree of individual can not be evaluated well.But method of the present invention is applied, stimuli responsive state is conceived to
Under heart responding ability, using the RR interphase variable quantity in 1s after sonic stimulation as evaluation index, it can be seen that hypertension man
Race's history is tested that there are significantly different with the functional parameters of Healthy subjects.Adopting said method, can be from the new side of heart responding ability
Method is set about, and is evaluated cardiac function, is improved existing health of heart degree detecting, make up the limitation of existing index.
By research, auditory stimulation can induce the cardiovascular typical scaring response of all subjects, i.e., sympathetic nerve is emerging
It puts forth energy and heart rate increases.But for the subject for having family history of hypertension, their cardiovascular response has significantly with normal subjects
Difference.Normal subjects will appear the temporary decline of heart rate after by auditory stimulation in 1s, i.e., RR interphase is temporary
Increase, but in the subject for having family history of hypertension, this responsiveness will appear decrease and even disappear.Brain electricity Evoked ptential
Processing find that the decrease of this responsiveness is related to the couple variations between heart and brain, the results showed that, brainstem auditory evoked
The α wave shock range of prefrontal cortex significantly less than normal subjects, and with RR interphase, elevated-levels are positively correlated in 1s,
The raising of this RR interphase instantaneous after auditory stimulation present in normal subjects is in the response of the normal brain activity heart and coupling process
A kind of common phenomenon, which can improve the parasympathetic excitement degree of heart, shield to heart, in order to avoid heart
Due to sympathetic nerve paradoxical discharge and induce arrhythmia cordis and other cardiovascular times.
Claims (9)
1. a kind of cardiac function evaluation method based on distance speech stimulation, which comprises the following steps:
Step 1) carries out multiple auditory stimulation experiment and records the electrocardiosignal for stimulating testee in preceding and stimulating course；
Step 2) is pre-processed to obtain RR interphase RR (n) by shooting to the electrocardiosignal of acquisition；
Step 3) carries out the calculating of RR interphase response sequence to the RR interphase RR (n) that step 2) obtains, and obtains normalized RR interphase
Sequence RR ' (t ', i)；
The normalized RR interval series RR ' (t ', i) that step 4), basis obtain calculates testee's interior RR interphase after stimulation
Average increase Δ RRIin, calculation method is as follows:
ΔRRIin=(RR'(0.5)+RR'(1))/2；
Step 5), by Δ RRIinAs the supplement index of basic cardiac function, cardiac function evaluation is carried out, if Δ RRIin≤
0, then it is assumed that Neuroprotective Mechanisms are not perfect.
2. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 1, which is characterized in that right
Collected original electro-cardiologic signals ecg (t) removes baseline drift and high-frequency noise first with wavelet decomposition, after obtaining denoising
Electrocardiosignal x (t)；Then the method for utilizing Wavelet Modulus Maxima, the small echo of selection are to support compactly support and disappear with single order
Square db2 small echo is lost, R wave crest is detected, then testing result is carried out to inspect modification, obtains a succession of continuous R peak time
The set R (N) of point；To the set R (N) of R peak time point, previous wave crest is subtracted using the latter peak time point R (n+1)
Time point R (n), RR interphase RR (n) by shooting is obtained, wherein (1≤n < N).
3. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 2, which is characterized in that right
Collected original electro-cardiologic signals ecg (t) selects db2 small echo to carry out 9 layers of decomposition to signal, removes the baseline drift of the last layer
The high-frequency noise of shifting and three first layers, the electrocardiosignal x (t) after being denoised.
4. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 2, which is characterized in that benefit
R wave crest is detected with the method for Wavelet Modulus Maxima, using support compactly support and with single order vanishing moment db2 small echo,
The equivalence filter coefficient of db2 small echo is as follows:
H1=0.3750, h2=0.1250, h3=0.0000
G1=0.5798.g2=0.0869, g3=0.0061
hk=h1-k,gk=-g1-k
If k > 3, hk=gk=0
Wavelet decomposition is made to the collected electrocardiosignal Mallat algorithm of institute.
5. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 1, which is characterized in that note
The time point t (i) of sonic stimulation is recorded, wherein i=1,2 ..., 20, each stimulation time point t (i) obtained in step 2) is attached
Close RR interval series carry out the resampling from the 2Hz in the times s such as t (i) time point front and back respectively, obtain sample rate unification
RR interval series RR (t, i), and the stimulation appearance point for newly obtaining RR interval series is defined as time point 0, data initial time
Point for-s, end time point is s, i.e. RR (t ', i).
6. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 5, which is characterized in that tool
Body, the resampling from t (i) -2.5s to the 2Hz of t (i)+2.5s is carried out respectively, obtains the unified RR interval series RR of sample rate
(t, i), and the stimulation appearance point for newly obtaining RR interval series is defined as time point 0, data start time point is -2.5, is terminated
Time point is 2.5, i.e. RR (t ', i).
7. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 1, which is characterized in that step
It is rapid 5) in, specifically, first carry out cardiac function diagnosis；
If testee's items cardiac function base values is normal, and index Δ RRIin> 0, then heart regulating power is strong；
If testee's items cardiac function base values is normal, and index Δ RRIin≤ 0, then heart regulating power is weak；
If testee's cardiac function base values has abnormal index, and index Δ RRIin≤ 0, then heart regulating power is poor.
8. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 1, which is characterized in that right
Obtained RR interval series RR (t ', i) is normalized, and obtains normalized RR interval series RR ' (t ', i):
Wherein, RR_baseline is average value of the RR interval series in -2.5 < t < 0 section.
9. a kind of cardiac function evaluation method based on distance speech stimulation according to claim 8, which is characterized in that will
RR'(t', i) it is averaging along stimulation number, obtain the RR interphase sonic stimulation response sequence RR'(t' of testee).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810783946.8A CN109044335B (en) | 2018-07-17 | 2018-07-17 | Heart function evaluation method based on instantaneous sound stimulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810783946.8A CN109044335B (en) | 2018-07-17 | 2018-07-17 | Heart function evaluation method based on instantaneous sound stimulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109044335A true CN109044335A (en) | 2018-12-21 |
| CN109044335B CN109044335B (en) | 2020-11-10 |
Family
ID=64816841
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810783946.8A Active CN109044335B (en) | 2018-07-17 | 2018-07-17 | Heart function evaluation method based on instantaneous sound stimulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109044335B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113712566A (en) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | Method and device for generating an interval of heart beat difference data sequence |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014071126A1 (en) * | 2012-11-01 | 2014-05-08 | The Johns Hopkins University | Method and system for determining strain relaxation of left ventricular diastolic function |
| CN103892812A (en) * | 2014-02-24 | 2014-07-02 | 南京丰生永康软件科技有限责任公司 | Matching degree analysis technology |
| US8862214B2 (en) * | 2008-04-04 | 2014-10-14 | Draeger Medical Systems, Inc. | Cardiac condition detection system |
| WO2016067119A1 (en) * | 2014-10-27 | 2016-05-06 | Uab Metapro Holding | System for the analysis of the daily heart rhythm autonomic nervous system balance |
| US20160256063A1 (en) * | 2013-09-27 | 2016-09-08 | Mayo Foundation For Medical Education And Research | Analyte assessment and arrhythmia risk prediction using physiological electrical data |
| CN105982661A (en) * | 2015-02-02 | 2016-10-05 | 四川理工学院 | Portable electrocardiogram monitoring system and data processing method thereof |
| US20160374580A1 (en) * | 2015-06-29 | 2016-12-29 | Kaunas University Of Technology | Method and System for Predicting of Acute Hypotensive Episodes |
| CN106419898A (en) * | 2016-08-12 | 2017-02-22 | 武汉中旗生物医疗电子有限公司 | Method removing electrocardiosignal baseline drift |
| CN106510737A (en) * | 2015-12-16 | 2017-03-22 | 西南大学 | Method for real-time detection of mental stress state during speech through heart rate measurement |
| CN106535742A (en) * | 2014-07-01 | 2017-03-22 | 美敦力公司 | Apparatus for verifying discriminating of tachycardia events in a medical device having dual sensing vectors |
| CN107595305A (en) * | 2017-09-18 | 2018-01-19 | 西南大学 | Anxiety state detection method and device |
| CN107688553A (en) * | 2017-08-16 | 2018-02-13 | 安徽心之声医疗科技有限公司 | Method based on wavelet transformation and logistic regression algorithm detection ecg wave form feature |
-
2018
- 2018-07-17 CN CN201810783946.8A patent/CN109044335B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8862214B2 (en) * | 2008-04-04 | 2014-10-14 | Draeger Medical Systems, Inc. | Cardiac condition detection system |
| WO2014071126A1 (en) * | 2012-11-01 | 2014-05-08 | The Johns Hopkins University | Method and system for determining strain relaxation of left ventricular diastolic function |
| US20160256063A1 (en) * | 2013-09-27 | 2016-09-08 | Mayo Foundation For Medical Education And Research | Analyte assessment and arrhythmia risk prediction using physiological electrical data |
| CN103892812A (en) * | 2014-02-24 | 2014-07-02 | 南京丰生永康软件科技有限责任公司 | Matching degree analysis technology |
| CN106535742A (en) * | 2014-07-01 | 2017-03-22 | 美敦力公司 | Apparatus for verifying discriminating of tachycardia events in a medical device having dual sensing vectors |
| US20170203114A1 (en) * | 2014-07-01 | 2017-07-20 | Medtronic, Inc. | Method and apparatus for verifying discriminating of tachycardia events in a medical device having dual sensing vectors |
| WO2016067119A1 (en) * | 2014-10-27 | 2016-05-06 | Uab Metapro Holding | System for the analysis of the daily heart rhythm autonomic nervous system balance |
| CN105982661A (en) * | 2015-02-02 | 2016-10-05 | 四川理工学院 | Portable electrocardiogram monitoring system and data processing method thereof |
| US20160374580A1 (en) * | 2015-06-29 | 2016-12-29 | Kaunas University Of Technology | Method and System for Predicting of Acute Hypotensive Episodes |
| CN106510737A (en) * | 2015-12-16 | 2017-03-22 | 西南大学 | Method for real-time detection of mental stress state during speech through heart rate measurement |
| CN106419898A (en) * | 2016-08-12 | 2017-02-22 | 武汉中旗生物医疗电子有限公司 | Method removing electrocardiosignal baseline drift |
| CN107688553A (en) * | 2017-08-16 | 2018-02-13 | 安徽心之声医疗科技有限公司 | Method based on wavelet transformation and logistic regression algorithm detection ecg wave form feature |
| CN107595305A (en) * | 2017-09-18 | 2018-01-19 | 西南大学 | Anxiety state detection method and device |
Non-Patent Citations (3)
| Title |
|---|
| P. SUN, Q. H. WU*,ET AL: "An Improved Morphological Approach to Background Normalization of ECG Signals", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
| 于晓琳: "意识控制心率过程中头皮脑电活动和心率变异性分析", 《中国医疗器械杂志》 * |
| 刘澄玉等: "心电间期序列归一化直方图的构造方法及其在评价心力衰竭中的应用", 《生物物理学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113712566A (en) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | Method and device for generating an interval of heart beat difference data sequence |
| CN113712566B (en) * | 2020-05-12 | 2024-02-06 | 深圳市科瑞康实业有限公司 | Method and device for generating heart beat interval difference value data sequence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109044335B (en) | 2020-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5351687A (en) | Method and apparatus for non-invasively determining a patient's susceptibility to ventricular arrhythmias | |
| Malberg et al. | Advanced analysis of spontaneous baroreflex sensitivity, blood pressure and heart rate variability in patients with dilated cardiomyopathy | |
| Schels et al. | Frequency analysis of the electrocardiogram with maximum entropy method for identification of patients with sustained ventricular tachycardia | |
| US20150223746A1 (en) | System with kinesthetic stimulation medical device for the non-invasive assessment of the sympathovagal balance of a patient | |
| COLOMBO et al. | Comparison between spectral analysis and the phenylephrine method for the assessment of baroreflex sensitivity in chronic heart failure | |
| Chung et al. | Noninvasive heart rate variability analysis using loadcell-installed bed during sleep | |
| CN109044335A (en) | A kind of cardiac function evaluation method based on distance speech stimulation | |
| US8255051B2 (en) | Skin response monitoring for neural and cardiac therapies | |
| AU2005261908B2 (en) | Defibrillator with cardiac blood flow determination | |
| Orini et al. | Detection of transient, regional cardiac repolarization alternans by time-frequency analysis of synthetic electrograms | |
| US10820820B2 (en) | Physiologic signal analysis using multiple frequency bands | |
| Mangin et al. | Simultaneous analysis of heart rate variability and myocardial contractility during head‐up tilt in patients with vasovagal syncope | |
| Akşahin et al. | Obstructive sleep apnea classification with artificial neural network based on two synchronic hrv series | |
| Partin et al. | Monitoring driver physiological parameters for improved safety | |
| Pyko et al. | Systolic blood pressure and pulse intervals synchronization | |
| CN113242716A (en) | Method, equipment and system for monitoring arrhythmia event | |
| Kuo et al. | Quantification of respiratory sinus arrhythmia using Hilbert–Huang transform | |
| Chen | Signal processing of electrocardiogram for arrhythmia prediction and classification | |
| Dambal | Statistical Analysis On Heart Rate Variability (Hrv) Features To Detect Congestive Heart Failure | |
| Tsukahara et al. | The recruitment pattern of single vasoconstrictor neurons in human | |
| Kenny et al. | Autonomic reflexes in patients with cardioinhibitory carotid sinus syncope | |
| Mejía-Rodríguez et al. | Time varying heart rate variability analysis of active orthostatic and cold face tests applied both independently and simultaneously | |
| Eftestøl et al. | Analysis of intracardiac electrogram changes | |
| Zhang et al. | Research of HRV and PRV under Electronic Synchronous Acquisition System | |
| Sema | The impact of continuous ECG-monitoring in patients with unexplained syncope |
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 |