CN102389302B - Analysis method of dynamic characteristics of electrocardiosignal - Google Patents
Analysis method of dynamic characteristics of electrocardiosignal Download PDFInfo
- Publication number
- CN102389302B CN102389302B CN2011102038673A CN201110203867A CN102389302B CN 102389302 B CN102389302 B CN 102389302B CN 2011102038673 A CN2011102038673 A CN 2011102038673A CN 201110203867 A CN201110203867 A CN 201110203867A CN 102389302 B CN102389302 B CN 102389302B
- Authority
- CN
- China
- Prior art keywords
- electrocardiosignal
- dynamic characteristic
- network
- max
- data
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000009826 distribution Methods 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims description 3
- 206010003658 Atrial Fibrillation Diseases 0.000 abstract description 4
- 238000000034 method Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 3
- 210000004903 cardiac system Anatomy 0.000 description 2
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005183 dynamical system Methods 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention relates to an analysis method of dynamic characteristics of electrocardiosignal, which comprises the following steps of: firstly, utilizing a surrogate data algorithm to carry out the dynamic characteristic recognition of collected signals; secondly, converting the recognized electrocardiosignal into a weighted network, and further capturing the difference between the dynamic characteristics of different types of electrocardiosignal in a framework of a complex network by a point intensity distribution map; and finally, defining a statistic Rs, and successfully distinguishing normal electrocardiosignal from atrial fibrillation electrocardiosignal based on the statistic. In the technical scheme provided by the invention, the classification accuracy of the electrocardisignal is improved, and the information contained in the electrocardisignal is deeply discovered.
Description
Technical field
The present invention relates to a kind of analytical method of electrocardiosignal, relate in particular to a kind of analytical method of electrocardiosignal dynamic characteristic.
Background technology
Heart disease has become one of main killer who threatens human health, therefore, develops that the automatic analytical tool of electrocardiosignal is significant fast and effectively.Just be based on this needs, nearly decades, people proposed the method for many analysis electrocardiosignaies: the sequence hypothesis check algorithm, the complexity algorithm, analysis of spectrum, method of wavelet analysis, these methods all are to extract feature from the time domain of electrocardiosignal or frequency domain, according to these features the electrocardiosignal that collects are carried out discriminator.Wherein, certain methods has been used to automatic defibrillator, and has obtained good effect.The non-linear dynamic feature of the cardiac system that but these methods all do not relate in the electrocardiosignal to be implied.In this simultaneously, some research worker begin to utilize Nonlinear Dynamics to analyze electrocardiosignal, and phase space reconfiguration has been widely applied in the ECG Signal Analysis as the bridge between the fluctuation of observation signal system time and the dynamical system space characteristics.But unlike signal is different to the requirement of phase space reconfiguration parameter (reconstruct dimension and time delay), the way that all adopts identical parameters to be reconstructed to all signals is difficult to the geometry of cardiac system is opened fully at present, thereby has limited the raising of electrocardiosignal analysis precision in phase space.
Summary of the invention
The technical problem that the present invention solves is: a kind of analytical method of electrocardiosignal dynamic characteristic is provided, overcome prior art analyze the nicety of grading of electrocardiosignal not high, fail deep layer ground and excavate the technical problem of the information that electrocardiosignal contains.
Technical scheme of the present invention is: a kind of analytical method of electrocardiosignal dynamic characteristic is provided, comprises the steps:
Gather electrocardiosignal: gather electrocardiosignal with electrocardiogram acquisition equipment, be designated as:
{ x
i, i=1,2,3 ..., n; Wherein: i represents sampling number;
Electrocardiosignal is carried out the dynamic characteristic analysis: the network that electrocardiosignal is mapped as weighting, each electrocardiosignal of gathering is a data segment, the network of described weighting is with a node in the corresponding described network of each data segment, between a pair of node the weights on limit be this to the distance between the corresponding data segment of node, be set to cut-point with R wave-wave peak position target electrocardiosignal section be cut into some data segment { s
1, s
2... s
m, a node in the corresponding described complex network of each data segment; Then, define the simple distance between a pair of data segment, and this distance is corresponded in the network weights on limit between corresponding node; Obtain the node intensity distribution of gained network;
Step 300: judge electrocardiosignal, that is: the node intensity distribution that obtains is carried out Gauss curve fitting, define statistic R on this basis
s
Further technical scheme of the present invention is: in gathering the electrocardiosignal step, also comprise and adopt the data alternate algorithm that the electrocardiosignal dynamic characteristic is identified.
Further technical scheme of the present invention is: adopt the data alternate algorithm to carry out the identification of electrocardiosignal dynamic characteristic and comprise the steps:
Given null hypothesis: suppose { x
iSequence is stochastic signal;
Generate alternate data: with signal { x
i, i=1,2,3 ..., n generates alternate data at random
J=1,2 ..., T;
Select for use correlation dimension as statistic: to select to embed dimension d from small to large, calculate the different { x under the dimensions that embed respectively
iCorrelation dimension D
0,
J=1,2 ..., the correlation dimension D of T
1, D
2..., D
T, D
1, D
2..., D
TAverage
With standard deviation D
VarAnd maximum D
MaxMinima D
Min
Dynamic characteristic identification: if D
0Not in the interval
In, then refuse null hypothesis; If D
0In the interval
In, then assert the hypothesis establishment.
Further technical scheme of the present invention is: electrocardiosignal is being carried out in the dynamic characteristic analytical procedure s between two data segments
iWith s
jBetween the determining of distance, adopt following formula:
Wherein: l
iBe s
iLength, l
jBe s
jLength, establish l
i≤ l
j, s
I, k, s
J, k+lBe respectively s
i, s
jK on the data segment, k+1 point.
Electrocardiosignal is carried out come resulting some intensity distribution of match with the Gaussian function model in the dynamic characteristic analysis in step, define statistic R on this basis
s:
Wherein: ω is corresponding some intensity level of Gauss curve fitting peak of curve, and pro (st) is that gained network mid point intensity size is the statistical probability of st,
The maximum of expression gained network node mid point intensity, m is the node number of gained network.
Technique effect of the present invention is: the present invention relates to a kind of analytical method of electrocardiosignal dynamic characteristic, at first utilize the alternate data algorithm to carry out dynamic characteristic identification to gathering the signal that comes; Then, the electrocardiosignal that identifies is converted into the complex network of weighting, and then under the framework of complex network, utilizes the some intensity distribution to catch the difference of dynamic characteristic between dissimilar electrocardiosignaies; At last, define statistic Rs, and successfully normal electrocardiosignal and atrial fibrillation electrocardiosignal are distinguished according to this statistic.Technical solution of the present invention has improved the nicety of grading of electrocardiosignal, and more deep layer ground excavates the information that electrocardiosignal contains.
Description of drawings
Fig. 1 is flow chart of the present invention.
Fig. 2 is electrocardiosignal dynamic characteristic identification process figure of the present invention.
The specific embodiment
Below in conjunction with specific embodiment, technical solution of the present invention is further specified.
As shown in Figure 1, the specific embodiment of the present invention is: a kind of analytical method of electrocardiosignal dynamic characteristic comprises the steps:
{ x
i, i=1,2,3 ..., n; Wherein: i represents sampling number.
Specific implementation process is as follows: gather electrocardiosignal with electrocardiogram acquisition equipment, in the embodiment of the invention, adopt the rhythm and pace of moving things of the electrocardiogram acquisition equipment records heart that singly leads.
Step 200: electrocardiosignal is carried out the dynamic characteristic analysis: the network that electrocardiosignal is mapped as weighting, each electrocardiosignal of gathering is a data segment, the network of described weighting is with a node in the corresponding described network of each data segment, between a pair of node the weights on limit be this to the distance between the corresponding data segment of node, be set to cut-point with R wave-wave peak position target electrocardiosignal section be cut into some data segment { s
1, s
2... s
m, a node in the corresponding described complex network of each data segment; Then, define the simple distance between a pair of data segment, and this distance is corresponded in the network weights on limit between corresponding node; Obtain the node intensity distribution of gained network;
Specific implementation process is as follows: be set to cut-point with R wave-wave peak position target electrocardiosignal section is cut into some data segments, might as well be designated as { s
1, s
2... s
m, each data segment corresponds to a node in the complex network; The weights on limit are that this is to the distance between the corresponding data segment of node between a pair of node.The weights on limit, that is: in the network between two nodes distance reflected the distance of corresponding data section corresponding track in phase space.And s between two data segments here
i, s
jThe determining of distance, adopt following formula:
Wherein: l
iBe s
iLength, l
jBe s
jLength, establish l
i≤ l
j, s
I, k, s
J, k+lBe respectively s
i, s
jK on the data segment, k+1 point.
Above conversion method is mapped as the nearer node of network middle distance with near the track same unstable cycle rail in the phase space, thereby the dynamic characteristic that contains in the electrocardiosignal has been embedded in the topological structure of complex network.
Utilize the some intensity distributions to investigate the different of atrial fibrillation signal and normal electrocardiosignal then.The point strength formula of i node is:
M and d
IjGiven identical in represented meaning and front.The statistical probability that corresponding some intensity level is st is:
Here,
At last, come resulting some intensity distribution of match with the Gaussian function model, namely use shape as G
j(x)=a
jExp (((x-μ
j)/σ
j)
2) Gaussian function stack obtain fitting function
Wherein: the number of function in the M representation model, parameter a
j, μ
j, σ
jEstimated to obtain by the non-linear least square of classics.Determine corresponding some intensity level ω of matched curve peak value, define statistic R on this basis
s:
Wherein: ω is corresponding some intensity level of Gauss curve fitting peak of curve, and pro (st) is that gained network mid point intensity size is the statistical probability of st,
The maximum of expression gained network node mid point intensity, m is the node number of gained network.
Wherein: ω is corresponding some intensity level of Gauss curve fitting peak of curve.
In the specific embodiment, training and test sample book all adopt the data among the MIT-BIH data base that Massachusetts Institute Technology provides.
Preferred implementation of the present invention is: in gathering the electrocardiosignal step, also comprise and adopt the data alternate algorithm that the electrocardiosignal dynamic characteristic is identified, as shown in Figure 2, adopt data alternate algorithm identification electrocardiosignal dynamic characteristic to comprise the steps:
Step 101: given null hypothesis: suppose { x
iSequence is stochastic signal.
Step 102: generate alternate data, that is: with signal { x
i, i=1,2,3 ..., n generates alternate data at random
J=1,2 ..., T.Detailed process is as follows: this signal is upset order at random generate T (generally selecting T 〉=30) group alternate data, the steps include: at first, generate one group and { x
iIsometric random number
Secondly, will
Ascending order is arranged, and the footnote sequence after the note ordering is (k
1, k
2..., k
n); Then, order again
At last, repeat above process T time, obtain T group alternate data:
J=1,2 ..., T.
Step 103: select for use correlation dimension as statistic, that is: select to embed dimension d from small to large, calculate under the different embedding dimensions { x respectively
iCorrelation dimension D
0,
J=1,2 ..., the correlation dimension D of T
1, D
2..., D
T, D
1, D
2..., D
TAverage
With standard deviation D
VarAnd maximum D
MaxMinima D
Min
Step 104: dynamic characteristic identification, that is: if D
0Not in the interval
In, then refuse null hypothesis, that is: this signal has more rich dynamic characteristic information; If D
0In the interval
In, then assert the hypothesis establishment, that is: this signal may be subjected to noise jamming bigger, if it is carried out the dynamic characteristic analysis, can cause erroneous judgement largely.
Technique effect of the present invention is: the present invention relates to a kind of analytical method of electrocardiosignal dynamic characteristic, at first utilize the alternate data algorithm to carry out dynamic characteristic identification to gathering the signal that comes; Then, the electrocardiosignal that identifies is converted into the complex network of weighting, and then under the framework of complex network, utilizes the some intensity distribution to catch the difference of dynamic characteristic between dissimilar electrocardiosignaies; At last, define statistic Rs, and successfully normal electrocardiosignal and atrial fibrillation electrocardiosignal are distinguished according to this statistic.Technical solution of the present invention has improved the nicety of grading of electrocardiosignal, and more deep layer ground excavates the information that electrocardiosignal contains.
Above content be in conjunction with concrete preferred implementation to further describing that the present invention does, can not assert that concrete enforcement of the present invention is confined to these explanations.For the general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, can also make some simple deduction or replace, all should be considered as belonging to protection scope of the present invention.
Claims (3)
1. the analytical method of an electrocardiosignal dynamic characteristic comprises the steps:
Gather electrocardiosignal: gather electrocardiosignal with electrocardiogram acquisition equipment, be designated as:
{ x
i, i=1,2,3 ..., n; Wherein: i represents sampling number;
Electrocardiosignal is carried out the dynamic characteristic analysis: electrocardiosignal is mapped as the network of weighting, to the electrocardiosignal of each collection, is set to cut-point with R wave-wave peak position target electrocardiosignal section is cut into some data segment { s
1, s
2... s
m, a node in the corresponding described complex network of each data segment; Then, define the simple distance between a pair of data segment, this simple distance is s between two data segments
iWith s
jBetween distance, adopt following formula:
Wherein: l
iBe s
iLength, l
jBe s
jLength, establish l
i≤ l
j, s
I, k, s
J, k+lBe respectively s
i, s
jK on the data segment, k+l point;
And this distance corresponded in the network weights on limit between corresponding node; Obtain the some intensity distribution of gained network; The point intensity distribution that obtains is carried out Gauss curve fitting, define statistic R on this basis
s,
2. according to the analytical method of the described electrocardiosignal dynamic characteristic of claim 1, it is characterized in that: in gathering the electrocardiosignal step, also comprise and adopt the data alternate algorithm that the electrocardiosignal dynamic characteristic is identified.
3. according to the analytical method of the described electrocardiosignal dynamic characteristic of claim 2, it is characterized in that: adopt the data alternate algorithm to carry out the identification of electrocardiosignal dynamic characteristic and comprise the steps:
Given null hypothesis: suppose { x
iSequence is stochastic signal;
Generate alternate data: with signal { x
i, i=1,2,3 ..., n generates alternate data at random
Select for use correlation dimension as statistic: to select to embed dimension d from small to large, calculate the different { x under the dimensions that embed respectively
iCorrelation dimension D
0,
Correlation dimension D
1, D
2..., D
T, D
1, D
2..., D
TAverage
With standard deviation D
VarAnd maximum D
MaxMinima D
Min
Dynamic characteristic identification: if D
0Not in the interval
In, then refuse null hypothesis; If D
0In the interval
In, then assert the hypothesis establishment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102038673A CN102389302B (en) | 2011-07-20 | 2011-07-20 | Analysis method of dynamic characteristics of electrocardiosignal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102038673A CN102389302B (en) | 2011-07-20 | 2011-07-20 | Analysis method of dynamic characteristics of electrocardiosignal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102389302A CN102389302A (en) | 2012-03-28 |
CN102389302B true CN102389302B (en) | 2013-07-10 |
Family
ID=45856791
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102038673A Expired - Fee Related CN102389302B (en) | 2011-07-20 | 2011-07-20 | Analysis method of dynamic characteristics of electrocardiosignal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102389302B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103071236B (en) * | 2012-12-31 | 2015-03-11 | 深圳邦健生物医疗设备股份有限公司 | Method and device for detecting electric shock signal |
CN103007431B (en) * | 2012-12-31 | 2015-05-20 | 深圳邦健生物医疗设备股份有限公司 | Device and method for shockable signal detection |
CN103136465A (en) * | 2013-03-06 | 2013-06-05 | 天津理工大学 | Method of using electrocardiosignals for identity recognition |
CN104523266B (en) * | 2015-01-07 | 2017-04-05 | 河北大学 | A kind of electrocardiosignal automatic classification method |
CN106562781A (en) * | 2016-05-20 | 2017-04-19 | 杨燕 | Novel intelligent therapeutic apparatus for neurology department |
CN106562782A (en) * | 2016-05-20 | 2017-04-19 | 彭慧敏 | Dedicated ECG monitor for pediatric nursing |
CN106038143A (en) * | 2016-07-08 | 2016-10-26 | 王宗梅 | Multifunctional nursing frame for department of cardiology |
CN106264597A (en) * | 2016-10-08 | 2017-01-04 | 张秀芸 | A kind of child uses CT protective clothing |
CN106510619B (en) * | 2016-10-11 | 2019-04-05 | 天津大学 | ECG Signal Analysis method based on complex network and in the application being intelligently worn by |
CN107260161B (en) * | 2017-07-18 | 2019-02-19 | 上海图灵医疗科技有限公司 | A kind of electrocardio dynamics data quantitative analysis method |
-
2011
- 2011-07-20 CN CN2011102038673A patent/CN102389302B/en not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
一种灵活的加权复杂网络演化模型及其仿真;苏凯;《系统仿真学报》;20090131;第21卷(第1期);266-271 * |
王鼐.非线性动力学方法在时间序列分析中的应用.《中国优秀博硕士学位论文全文数据库 (博士) 基础科学辑》.2005,A002-18. |
苏凯.一种灵活的加权复杂网络演化模型及其仿真.《系统仿真学报》.2009,第21卷(第1期),266-271. |
非线性动力学方法在时间序列分析中的应用;王鼐;《中国优秀博硕士学位论文全文数据库 (博士) 基础科学辑》;20051115;A002-18 * |
Also Published As
Publication number | Publication date |
---|---|
CN102389302A (en) | 2012-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102389302B (en) | Analysis method of dynamic characteristics of electrocardiosignal | |
CN109171707A (en) | A kind of intelligent cardiac figure classification method | |
CN106108889B (en) | Electrocardiogram classification method based on deep learning algorithm | |
CN106658590A (en) | Design and implementation of multi-person indoor environment state monitoring system based on WiFi channel state information | |
CN104063645B (en) | A kind of personal identification method based on the dynamic self refresh sample of electrocardiosignal | |
CN107358250B (en) | Body gait recognition methods and system based on the fusion of two waveband radar micro-doppler | |
CN107290741B (en) | Indoor human body posture identification method based on weighted joint distance time-frequency transformation | |
CN102512158B (en) | Electrocardiosignal quality evaluation method and device based on high-dimensional fuzzy recognition | |
CN109171708B (en) | Defibrillation heart rhythm recognition device | |
CN105997055A (en) | Automatic classification method, system and device of electrocardiosignal ST band | |
CN102085095A (en) | Method, system and electrocardioscanner for detecting ST segment in electrocardiogram | |
CN107945817A (en) | Heart and lung sounds signal sorting technique, detection method, device, medium and computer equipment | |
CN103699873A (en) | Lower-limb flat ground walking gait recognition method based on GA-BP (Genetic Algorithm-Back Propagation) neural network | |
CN103049741A (en) | Foot-to-ground acting force-based gait feature extraction method and gait identification system | |
CN103093133B (en) | Biological identity authentication method facing institute of electrical and electronic engineers (IEEE) 802.15.6 | |
Sidek et al. | ECG biometric recognition in different physiological conditions using robust normalized QRS complexes | |
CN102961129A (en) | Method for analyzing abnormal electrocardiogram tension for remote medical care | |
CN104887222A (en) | Reversible electroencephalogram analysis method | |
CN105139029A (en) | Activity recognition method and activity recognition device for persons serving sentences | |
CN111597991A (en) | Rehabilitation detection method based on channel state information and BilSTM-Attention | |
CN112364296B (en) | P-wave arrival time automatic pickup method based on deep learning | |
CN112957056B (en) | Method and system for extracting muscle fatigue grade features by utilizing cooperative network | |
CN103886184A (en) | Construction method for heart pathology recognition model | |
CN110059384B (en) | Method for simulating crowd jump load based on generating type confrontation network | |
CN114159079B (en) | Multi-type muscle fatigue detection method based on feature extraction and GRU deep learning model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |