CN107616786A - Pulse wave recognition methods and pulse wave identification device - Google Patents
Pulse wave recognition methods and pulse wave identification device Download PDFInfo
- Publication number
- CN107616786A CN107616786A CN201711001882.3A CN201711001882A CN107616786A CN 107616786 A CN107616786 A CN 107616786A CN 201711001882 A CN201711001882 A CN 201711001882A CN 107616786 A CN107616786 A CN 107616786A
- Authority
- CN
- China
- Prior art keywords
- wave
- pulse wave
- dicrotic
- identification
- main ripple
- 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.)
- Withdrawn
Links
Abstract
The invention provides a kind of pulse wave recognition methods and pulse wave identification device.Wherein, this method includes：Original pulse wave is filtered, to obtain dicrotic wave identification curve；The crest position of curve and the main ripple identified in advance identification dicrotic wave is identified according to dicrotic wave, to identify dicrotic wave.In the present invention, on the basis of the identification of main ripple, on the one hand strengthen amplitude of the signal in dicrotic wave position by high-pass filtering, so as to effectively prevent the changeable interference to dicrotic wave waveform recognition of pulse wave form, the accuracy rate and identification crest, the precision of wave trough position of dicrotic wave identification can be effectively improved；On the other hand, by the way that the decent of pulse wave to be resolved into the decent of main ripple and being superimposed for tidal wave ascending branch, using difference method, the position where tidal wave crest is identified, and then has efficiently identified out tidal wave, providing accurate, comprehensive data for clinical diagnosis supports.
Description
Technical field
The present invention relates to pulse wave identification technology field, in particular to a kind of pulse wave recognition methods and pulse wave identification dress
Put.
Background technology
Digitalized pulse diagnosis is to collect pulse wave signal by the pressure sensor being attached on radial artery.This technology makes
Quantification, standardization are progressively moved towards in diagnosis by feeling the pulse, and more objective, accurate clinical diagnosis, auxiliary diagnosis, curative effect are provided for pulse wave spectrum
The technical support of judgement, and then promoted the process of Chinese medicine modernization.The mainstream technology framework of Contemporary Digital diagnosis by feeling the pulse is from arteries and veins
The key positions such as main ripple, tidal wave, dicrotic notch, dicrotic wave are extracted on ripple, then pulse wave shape parameter are calculated by these positions, based on this
A little pulse wave parameters provide for diagnosis arteries and veins type quantifies basis for estimation.
The pulse measurement of the people of different constitutions, even same person different times, can all there is difference physiologically, it is this
Difference can cause the change of pulse wave signal, pulse wave signal is showed obvious high complexity.The pulse wave identification of current main-stream
Algorithm is directly to be identified on original pulse wave signal based on morphological method, and a small amplitude change is likely to lead
Cause recognition result significantly deviation.Current algorithm tackles the height of pulse wave signal using excessive parameter and frequently parameter adjustment
Complicated and form is changeable, generally existing bad adaptability, identifies the problems such as unstable, parameter system is numerous and jumbled, computation complexity is too high.
The content of the invention
In consideration of it, the present invention proposes a kind of pulse wave recognition methods and pulse wave identification device, it is intended to solves current pulse wave and knows
The problem of accuracy that method for distinguishing identifies to pulse wave is low.
On one side, the present invention proposes a kind of pulse wave recognition methods, and this method comprises the following steps：Original pulse wave is entered
Row filtering, to obtain dicrotic wave identification curve；The ripple of curve and the main ripple identified in advance identification dicrotic wave is identified according to dicrotic wave
Peak position, to identify dicrotic wave.
Further, in above-mentioned pulse wave recognition methods, the recognition methods of the crest position of dicrotic wave is：By dicrotic wave
The crest of identification curve is denoted as k, and the starting point of next main ripple of the crest is denoted as into v；On k-v sections, the pole of minimum is identified
Small value x1And minimum x1First maximum x afterwards2, maximum x2The crest position of corresponding dicrotic wave.
Further, in above-mentioned pulse wave recognition methods, the scope of filtering is less than 4Hz and is more than 10Hz.
Further, in above-mentioned pulse wave recognition methods, also include after dicrotic wave is identified：Tidal wave is identified according to main ripple.
Further, in above-mentioned pulse wave recognition methods, the recognition methods of tidal wave is：The ascending branch of main ripple is carried out on master
The mirror image of the crest of ripple, to obtain conversion pulse wave；According to original pulse wave signal and conversion pulse wave signal identification tidal wave.
Further, it is above-mentioned according to original pulse wave signal and conversion pulse wave signal identification tide in above-mentioned pulse wave recognition methods
The method of ripple is：Original pulse wave signal and conversion pulse wave signal are subtracted each other, to obtain difference；Difference is carried out to difference, to obtain
Tidal wave identifies curve；On k-v sections, maximum is identified, position corresponding to the maximum is the crest position of tidal wave.
Further, in above-mentioned pulse wave recognition methods, the recognition methods of main ripple is：Original pulse wave is pre-processed, with
Obtain main ripple identification curve；The crest position of the main ripple of Curves Recognition is identified according to main ripple, to identify main ripple.
Further, in above-mentioned pulse wave recognition methods, the pretreatment to original pulse wave is：The signal of original pulse wave is carried out
Timing Difference, to obtain difference sequence；Each value progress square to difference sequence, to obtain difference square；Difference is put down
Fang Jinhang convolution, to obtain main ripple identification curve.
Further, in above-mentioned pulse wave recognition methods, the recognition methods of the crest position of main ripple is：Identify that main ripple is known
The maximum of other curve, the maximum point correspond to the steepest rising point of the ascending branch of main ripple；Respectively in the both sides of steepest rising point
Preset time range in identification minimum value, the waveform between two minimum values is main ripple；Identified most in preset time range
Big value, maximum are the crest of main ripple.
Further, in above-mentioned pulse wave recognition methods, preset time range 0.02s-0.5s.
On the other hand, the invention also provides a kind of pulse wave identification device, the device include：Filtration module, for original
Beginning pulse wave is filtered, to obtain dicrotic wave identification curve；Identification module, for identifying curve and in advance identification according to dicrotic wave
Main ripple identification dicrotic notch position and dicrotic wave crest position, to identify dicrotic wave.
In the present invention, on the basis of the identification of main ripple, signal shaking in dicrotic wave position is on the one hand made by high-pass filtering
Width strengthens, and so as to effectively prevent the changeable interference to dicrotic wave waveform recognition of pulse wave form, and is filtering enhanced signal
It is upper carry out dicrotic wave identification, can effectively improve dicrotic wave identification accuracy rate and identification crest, wave trough position it is accurate
Degree；On the other hand, caused by the stacking pattern that the morphological differences of the decent of pulse wave and ascending branch is regarded as to different ripples, pass through
The decent of pulse wave is resolved into the decent of main ripple and being superimposed for tidal wave ascending branch, using difference method, identifies tidal wave crest
The position at place, and then tidal wave has been efficiently identified out, and tidal wave recognition methods provided by the invention has obvious crest in tidal wave
Or protrusion it is unconspicuous in the case of can keep very high accuracy rate, provide accurate, comprehensive data branch for clinical diagnosis
Hold.In addition, compared with prior art, the parameter knot required for dicrotic wave recognition methods and tidal wave recognition methods provided by the invention
Structure is simple, and then reduces computation complexity.
Brief description of the drawings
By reading the detailed description of hereafter preferred embodiment, it is various other the advantages of and benefit it is common for this area
Technical staff will be clear understanding.Accompanying drawing is only used for showing the purpose of preferred embodiment, and is not considered as to the present invention
Limitation.And in whole accompanying drawing, identical part is denoted by the same reference numerals.In the accompanying drawings：
Fig. 1 is the flow chart of pulse wave recognition methods provided in an embodiment of the present invention；
Fig. 2 is the pulse wave schematic diagram of a cycle provided in an embodiment of the present invention；
Fig. 3 is original waveform provided in an embodiment of the present invention and the Contrast on effect schematic diagram of waveform after filtering；
Fig. 4 is the flow chart of the recognition methods of the crest position of dicrotic wave provided in an embodiment of the present invention；
Fig. 5 is that dicrotic wave provided in an embodiment of the present invention identifies curve synoptic diagram；
Fig. 6 is the unconspicuous taut pulse of one group of dicrotic wave provided in an embodiment of the present invention and dicrotic notch-dicrotic wave recognition result
Schematic diagram；
Fig. 7 is the another flow chart of pulse wave recognition methods provided in an embodiment of the present invention；
Fig. 8 is the flow chart of the recognition methods of tidal wave provided in an embodiment of the present invention；
Fig. 9 is the signal for the mirror image that the ascending branch of main ripple provided in an embodiment of the present invention carries out the crest on main ripple
Figure；
Figure 10 be it is provided in an embodiment of the present invention according to original pulse wave signal and it is described conversion signal identification tidal wave method
Flow chart；
Figure 11 is that tidal wave provided in an embodiment of the present invention identifies curve synoptic diagram；
Figure 12 is the flow chart of main ripple recognition methods provided in an embodiment of the present invention；
Figure 13 is the flow chart provided in an embodiment of the present invention pre-processed to original pulse wave；
Figure 14 is the schematic diagram that main ripple provided in an embodiment of the present invention identifies curve；
Figure 15 is the flow chart of the recognition methods of the crest position of main ripple provided in an embodiment of the present invention；
Figure 16 is the structured flowchart of pulse wave identification device provided in an embodiment of the present invention.
Embodiment
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
Limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure
Completely it is communicated to those skilled in the art.It should be noted that in the case where not conflicting, embodiment in the present invention and
Feature in embodiment can be mutually combined.Describe the present invention in detail below with reference to the accompanying drawings and in conjunction with the embodiments.
Embodiment of the method：
Referring to Fig. 1, Fig. 1 is the flow chart for the pulse wave recognition methods that the present embodiment provides.As illustrated, this method is included such as
Lower step：
Step S110, original pulse wave is filtered, to obtain dicrotic wave identification curve.
Specifically, main ripple 1, tidal wave 2, dicrotic notch 3 and dicrotic wave 4, arteries and veins are included successively referring to Fig. 2, a Pulse period
The starting point of ripple to main ripple 1 crest to rise to, it to the terminal of pulse wave is descending branch that the crest of main ripple 1, which is,.By observing pulse wave signal
Understand, because the form of dicrotic wave 4 has diversity, so if directly knowing on original pulse wave signal to dicrotic wave 4
Not, higher recognition accuracy is extremely difficult to, also, the pulse wave signal of actual acquisition and the structure of theory make a big difference,
Some, or even most of waveform None- identifieds are might have, and because dicrotic wave 4 is blood flow ripple caused by acral bounce-back
Dynamic, amplitude is smaller, and in the pulse wave signal of many types, dicrotic wave 4 may show as the small projection in decent, so that knowing
Do not become extremely difficult.It is therefore desirable to carry out high-pass filtering to original pulse wave signal, i.e., original pulse wave signal is protruded
Radio-frequency component suppresses the processing of low-frequency component, and high-pass filtering can obtain amplitude of the original pulse wave signal in the position of dicrotic wave 4
To enhancing.Waveform A' after being filtered is filtered to original waveform A and may refer to Fig. 3, as illustrated, waveform A' after filtering
On, the fluctuation of the correspondence position of dicrotic wave 4 becomes apparent from.Dicrotic wave identification curve can be obtained after filtering, and the curve can
Using the feature as identification dicrotic wave 4.When it is implemented, original pulse wave can be filtered again after main ripple 1 is identified.
Step S120, identify that curve and the main ripple identified in advance identify the crest position of dicrotic wave according to dicrotic wave,
To identify dicrotic wave.
Specifically, identify that curve and the main ripple 1 identified in advance identify dicrotic wave 4 according to the dicrotic wave of above-mentioned acquisition
Crest position, that is, dicrotic wave 4 is identified, it is also possible to identify the position of dicrotic notch 3.
As can be seen that in the present embodiment, first by being filtered to original pulse wave, original pulse wave signal can be made in weight
The amplitude of position of ripple 4 of fighting is strengthened, and then can accurately obtain dicrotic wave identification curve；Then identified according to dicrotic wave
Curve and the main ripple 1 identified in advance can identify the crest position of dicrotic wave 4, and then realize and accurately identify dicrotic wave 4
Purpose.
Referring to Fig. 4, Fig. 4 is the recognition methods of the crest position for the dicrotic wave that the present embodiment provides.As illustrated, should
Method can include：
Step S410, the crest that dicrotic wave identifies curve is denoted as k, the starting point of next main ripple of the crest is denoted as v.
Specifically, the crest for the dicrotic wave identified in above-described embodiment being identified into curve is denoted as k, and by the crest
The starting point of next main ripple 1 be denoted as v.
Step S420, on k-v sections, identify the minimum x of minimum1And minimum x1First maximum afterwards
x2, maximum x2The crest position of corresponding dicrotic wave.
Specifically, as shown in figure 5, on k-v sections, the minimum x of minimum is identified1, and identify minimum x1Afterwards
First maximum x2, minimum x1The corresponding position of dicrotic notch 3, maximum x2The crest institute for then corresponding to dicrotic wave 4 is in place
Put.Compare and understand with Fig. 6, the dicrotic wave recognition methods that the present embodiment is provided is on the pulse wave signal of a variety of pulse condition types
All achieve good adaptability.
In above-described embodiment, energy ingredient of the dicrotic wave 4 on frequency spectrum is mainly distributed on 4Hz-10Hz, main ripple 1, tidal wave 2
And the energy ingredient of other positions is below this scope, so, filter range when being filtered to original pulse wave can be with
Less than 4Hz and it is more than 10Hz.
Referring to Fig. 7, Fig. 7 is the another flow chart for the pulse wave recognition methods that the present embodiment provides.As illustrated, this method can
With including：
Step S710, original pulse wave is filtered, to obtain dicrotic wave identification curve.
Step S720, identify that curve and the main ripple identified in advance identify the crest position of dicrotic wave according to dicrotic wave,
To identify dicrotic wave.
Step S730, tidal wave is identified according to main ripple.
Specifically, during pulse wave is identified, main ripple 1 can be first identified in advance, then identifies dicrotic wave 4 again, most
Tidal wave 2 is identified afterwards, and then can ensure the accuracy that tidal wave 2 identifies.And tidal wave 2 can be carried out according to the main ripple 1 identified in advance
Identification.
Specifically, Fig. 8 is may refer to, Fig. 8 is the flow chart of the recognition methods for the tidal wave that the present embodiment provides.As schemed
Show, this method can include：
Step S810, the mirror image of the crest on main ripple is carried out to the ascending branch of main ripple, to obtain conversion signal.
Specifically, in the signal of actual acquisition, the pulse wave signal of many types be present, wherein, tidal wave 2 and main ripple 1
Positional distance it is closer, or even fusion occurs.In this case, accurately identifying for tidal wave 2 can become extremely difficult.By ripple
The superposition theory of shape understands that the superposition of ripple can cause the deformation of curve.Theoretically, if main ripple 1 is not disturbed by other waveforms,
Its decent will be with ascending branch in the form of symmetrical or near symmetrical.On the decent of main ripple 1, due to the rising of tidal wave 2
The decent superposition of Zhi Yuzhu ripples 1, the symmetry of waveform can be by bigger interference.Become based on waveform caused by this interference
Change feature, when identifying tidal wave 2, referring to Fig. 9, the ascending branch B of the main ripple 1 to identifying in advance first makees the crest mirror on main ripple 1
Face maps, and obtains the mapping branch B' after mirror image, and other positions keep constant, and then obtains converting pulse wave signal.
Step S820, according to original pulse wave signal and conversion pulse wave signal identification tidal wave.
Specifically, tidal wave can be identified according to original pulse wave signal and conversion signal, specifically may refer to Figure 10,
As illustrated, it can be included according to the method for original pulse wave signal and conversion pulse wave signal identification tidal wave：
Step S1010, original pulse wave signal and conversion pulse wave signal are subtracted each other, to obtain difference.
Specifically, referring to Figure 11, original pulse wave signal is denoted as y0, by the conversion pulse wave after the conversion of crest minute surface
Signal is denoted as y1.By y0And y1Make the difference, the calculation formula for obtaining converting difference d, d is：D=y0–y1。
Step S1020, difference is carried out to difference, to obtain tidal wave identification curve.
Specifically, d is subjected to difference, to be identified the feature f used in tidal wave 2d。fdCalculation formula be：fd(i)
=d (i)-d (i-1) (i=1,2 ...), that is, obtain tidal wave identification curve C.
Step S1030, on k-v sections, identify maximum, position corresponding to the maximum for tidal wave crest institute in place
Put.
Specifically, due to having been realized in main ripple crest k, dicrotic notch v identification in above-described embodiment, therefore, to tide
The scope of the crest position identification of ripple 2 can be narrowed down between k and v.Identified on tidal wave identification curve C k-v sections
Maximum, position corresponding to the maximum are the crest position of tidal wave 2.
As can be seen that in the present embodiment, converted first by doing crest minute surface to original pulse wave signal, then by original pulse wave
Signal makes the difference with pulse wave signal after conversion, then carries out difference to difference, to obtain tidal wave identification curve, is finally identified in tidal wave
Tidal wave 2 is identified on curve C, tidal wave 2 is identified using the method, other kinds of arteries and veins signal wave can be avoided to tidal wave 2
The interference of identification, it ensure that the accuracy that tidal wave 2 identifies.
Referring to Figure 12, Figure 12 is the flow chart for the main ripple recognition methods that the present embodiment provides.As illustrated, this method can be with
Including：
Step S1210, original pulse wave is pre-processed, to obtain main ripple identification curve.
Specifically, because pulse wave signal has form polytropy, if directly detecting the ripple of main ripple 1 on pulse wave signal
Paddy-crest, can factor be influenceed by baseline drift, dicrotic pulse wave-amplitude be excessive etc..And the energy of baseline drift and dicrotic wave 4 exists
Low frequency part is focused primarily upon on frequency spectrum, therefore first has to pre-process original pulse wave, i.e., original pulse wave signal is carried out
Prominent radio-frequency component simultaneously suppresses the processing of low-frequency component, and then obtains main ripple identification curve.
Step S1220, the crest position of the main ripple of Curves Recognition is identified according to main ripple, to identify main ripple.
Specifically, the crest position of main ripple 1 is identified on the main ripple identification curve of above-mentioned acquisition, and then identifies master
Ripple 1.
Referring to Figure 13, Figure 13 is the flow chart pre-processed to original pulse wave that the present embodiment provides.It is as illustrated, pre-
Processing includes：
Step S1310, Timing Difference is carried out to the signal of original pulse wave, to obtain difference sequence.
Specifically, original pulse wave signal is denoted as y, y is sampling time t function, i.e. y=f (t).Sequential is carried out to y
Difference, obtain difference sequence dy.The calculation of its i-th of value is：
Dy (i)=y (i)-y (i-1) (i=1,2 ..., t, t+1, t+2 ...).
Step S1320, each value progress square to difference sequence, to obtain difference square.
Specifically, in order to protrude the radio-frequency component of signal, difference sequence dy each value is carried out square, and then
Obtain difference square dy2.Calculation formula is：
dy2(i)=dy (i) × dy (i)
Step S1330, convolution is carried out to difference square, to obtain main ripple identification curve.
Specifically, the signal after difference square still can have shake in high frequency position, and then crest can be caused to detect
Stability is poor.Therefore, the rectangular window of seclected time width is introduced on the signal after difference square, to enter line slip convolution,
And then smoother curve can be obtained.When it is implemented, convolution time width is adjustable, preferably 0.05s-0.25s.Referring to figure
14, using the smoothed curve that convolution obtains as the feature for identifying main ripple 1, i.e., main ripple identification curve E.
As can be seen that in the present embodiment, by being pre-processed to original pulse wave, i.e., by successively to original pulse wave signal
Carry out difference, carry out convolution to difference sequence progress square and to difference square, it is possible to achieve original pulse wave signal is dashed forward
Go out radio-frequency component and suppress low-frequency component, to obtain main ripple identification curve E, and then the accuracy of main ripple identification can be ensured, also entered
One step ensure that the identification accuracy of follow-up dicrotic wave 4 and tidal wave 2.
Referring to Figure 15, Figure 15 is the flow chart of the recognition methods of the crest position for the main ripple that the present embodiment provides.Such as
Shown in figure, this method comprises the following steps：
Step S1510, identifies the maximum of main ripple identification curve, and the maximum point is corresponded on the steepest of the ascending branch of main ripple
Rise point.
Specifically, on the main ripple identification curve obtained in the above-described embodiments, the maximum of main ripple identification curve is identified,
The maximum point corresponds to the steepest rising point of the ascending branch of main ripple 1, and the steepest rising point is that signal intensity increase is most in ascending branch
Fast position.
Step S1520, minimum value, two minimum values are identified in the preset time range of the both sides of steepest rising point respectively
Between waveform be main ripple.
Specifically, a minimum value is identified in the preset time range in the left side of steepest rising point, equally, in steepest
One minimum value of identification, i.e., in preset time range, can recognize pulse wave in the preset time range on the right side of rising point
The beginning and end in cycle, therefore it is main ripple 1 that the waveform between two minimum values, which is,.When it is implemented, preset time range can be with
For 0.02s-0.5s.
Step S1530, identifies maximum in preset time range, and maximum is the crest of main ripple.
Specifically, maximum is identified in preset time range, the maximum is the crest for main ripple 1.
As can be seen that in the present embodiment, steepest rising point is identified on curve E by being identified in main ripple, master can be identified
Ripple 1, and then the crest of main ripple 1 is identified, it is simple and easy.
To sum up, in the present embodiment, on the basis of the identification of main ripple, signal is on the one hand made in dicrotic wave position by high-pass filtering
The amplitude enhancing put, so as to effectively prevent the changeable interference to dicrotic wave waveform recognition of pulse wave form, and after filtering strengthens
Signal on carry out dicrotic wave identification, can effectively improve dicrotic wave identification accuracy rate and identification crest, wave trough position
Precision；On the other hand, the stacking pattern for regarding the morphological differences of the decent of pulse wave and ascending branch as different ripples causes
, by the way that the decent of pulse wave to be resolved into the decent of main ripple and being superimposed for tidal wave ascending branch, utilize difference method, identification tide
Position where ripple crest, and then tidal wave has been efficiently identified out, and the tidal wave recognition methods that the present embodiment provides has in tidal wave
Very high accuracy rate can be kept in the case of obvious crest or protrusion are unconspicuous, is provided for clinical diagnosis accurate, comprehensive
Data support.In addition, compared with prior art, needed for the dicrotic wave recognition methods and tidal wave recognition methods that the present embodiment provides
The argument structure wanted is simple, and then reduces computation complexity.
Device embodiment：
Referring to Figure 16, Figure 16 is the structured flowchart for the pulse wave identification device that the present embodiment provides.As illustrated, the device can
With including：Filtration module 100 and identification module 200.Wherein, filtration module 100 is used to be filtered original pulse wave, to obtain
Dicrotic wave identifies curve.Identification module 200 is used to identify that curve and the main ripple identified in advance identify dicrotic notch institute according to dicrotic wave
In position and the crest position of dicrotic wave, to identify dicrotic wave.Wherein, the specific implementation process side of being referred to of the device
Implementation process in method embodiment, here is omitted.
As can be seen that in the present embodiment, first by being filtered to original pulse wave, original pulse wave signal can be made in weight
The amplitude of position of ripple of fighting is strengthened, and then can accurately obtain dicrotic wave identification curve；Then identified according to dicrotic wave
Curve and the main ripple identified in advance can identify the crest position of dicrotic wave, and then realize the mesh for accurately identifying dicrotic wave
's.
It should be noted that the present embodiment provide pulse wave recognition methods it is identical with pulse wave identification device principle, it is similar it
Place can be with cross-referenced.
To sum up, in the present embodiment, first by being filtered to original pulse wave, original pulse wave signal can be made in dicrotic wave
The amplitude of position strengthened, and then can accurately obtain dicrotic wave identification curve；Then curve is identified according to dicrotic wave
The main ripple identified in advance can identify the crest position of dicrotic wave, and then realize the purpose for accurately identifying dicrotic wave.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention
God and scope.So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to comprising including these changes and modification.
Claims (11)
1. a kind of pulse wave recognition methods, it is characterised in that comprise the following steps：
Original pulse wave is filtered, to obtain dicrotic wave identification curve；
The crest position of curve and the main ripple identified in advance identification dicrotic wave is identified according to the dicrotic wave, with described in identification
Dicrotic wave.
2. pulse wave recognition methods according to claim 1, it is characterised in that the knowledge of the crest position of the dicrotic wave
Other method is：
The crest of dicrotic wave identification curve is denoted as k, the starting point of next main ripple of the crest is denoted as v；
On k-v sections, the minimum x of minimum is identified1And minimum x1First maximum x afterwards2, the maximum x2
The crest position of the corresponding dicrotic wave.
3. pulse wave recognition methods according to claim 1, it is characterised in that
The scope of the filtering is less than 4Hz and is more than 10Hz.
4. pulse wave recognition methods according to claim 2, it is characterised in that also include after the dicrotic wave is identified：
According to the main ripple identification tidal wave.
5. pulse wave recognition methods according to claim 4, it is characterised in that the recognition methods of the tidal wave is：
The mirror image of crest on the main ripple is carried out to the ascending branch of the main ripple, to obtain conversion pulse wave；
According to original pulse wave signal and the conversion pulse wave signal identification tidal wave.
6. pulse wave recognition methods according to claim 5, it is characterised in that above-mentioned according to original pulse wave signal and the change
Change pulse wave signal and identify that the method for tidal wave is：
The original pulse wave signal and the conversion pulse wave signal are subtracted each other, to obtain difference；
Difference is carried out to the difference, to obtain tidal wave identification curve；
On the k-v sections, maximum is identified, position corresponding to the maximum is the crest position of the tidal wave.
7. pulse wave recognition methods according to any one of claim 1 to 6, it is characterised in that the identification side of the main ripple
Method is：
The original pulse wave is pre-processed, to obtain main ripple identification curve；
According to the crest position of main ripple described in the main ripple identification Curves Recognition, to identify the main ripple.
8. pulse wave recognition methods according to claim 7, it is characterised in that the pretreatment to the original pulse wave is：
Timing Difference is carried out to the signal of the original pulse wave, to obtain difference sequence；
Each value progress square to the difference sequence, to obtain difference square；
Convolution is carried out to the difference square, to obtain the main ripple identification curve.
9. the pulse wave recognition methods according to claim 7 or 8, it is characterised in that the crest position of the main ripple
Recognition methods is：
The maximum of the main ripple identification curve is identified, the maximum point corresponds to the steepest rising point of the ascending branch of the main ripple；
Minimum value is identified in the preset time range of the both sides of the steepest rising point respectively, between two minimum values
Waveform is the main ripple；
Maximum is identified in the preset time range, the maximum is the crest of the main ripple.
10. pulse wave recognition methods according to claim 9, it is characterised in that
The preset time range is 0.02s-0.5s.
11. a kind of pulse wave identification device of pulse wave recognition methods using any one of claim 1-10, its feature exist
In, including：
Filtration module, for being filtered to original pulse wave, to obtain dicrotic wave identification curve；
Identification module, for identifying that curve and the main ripple identified in advance identify dicrotic notch position and institute according to the dicrotic wave
The crest position of dicrotic wave is stated, to identify the dicrotic wave.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711001882.3A CN107616786A (en) | 2017-10-24 | 2017-10-24 | Pulse wave recognition methods and pulse wave identification device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711001882.3A CN107616786A (en) | 2017-10-24 | 2017-10-24 | Pulse wave recognition methods and pulse wave identification device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107616786A true CN107616786A (en) | 2018-01-23 |
Family
ID=61091925
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711001882.3A Withdrawn CN107616786A (en) | 2017-10-24 | 2017-10-24 | Pulse wave recognition methods and pulse wave identification device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107616786A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108403094A (en) * | 2018-03-23 | 2018-08-17 | 中国科学院微电子研究所 | A method of identification pulse wave wave crest |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050010116A1 (en) * | 2003-07-07 | 2005-01-13 | Ilkka Korhonen | Method and apparatus based on combination of physiological parameters for assessment of analgesia during anesthesia or sedation |
| CN101732033A (en) * | 2008-11-07 | 2010-06-16 | 中国科学院计算技术研究所 | Method and device for extracting characteristic parameter in human body waveform |
| CN102293639A (en) * | 2011-06-30 | 2011-12-28 | 芜湖圣美孚科技有限公司 | Pulse condition signal time domain feature extraction method |
| CN105286815A (en) * | 2015-11-02 | 2016-02-03 | 重庆大学 | Pulse wave signal feature point detection method based on waveform time domain features |
| CN107432736A (en) * | 2017-06-06 | 2017-12-05 | 新绎健康科技有限公司 | A kind of method for identifying pulse wave signal |
| GB2552455A (en) * | 2016-06-16 | 2018-01-31 | Digital & Future Tech Limited | Blood monitoring |
-
2017
- 2017-10-24 CN CN201711001882.3A patent/CN107616786A/en not_active Withdrawn
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050010116A1 (en) * | 2003-07-07 | 2005-01-13 | Ilkka Korhonen | Method and apparatus based on combination of physiological parameters for assessment of analgesia during anesthesia or sedation |
| CN101732033A (en) * | 2008-11-07 | 2010-06-16 | 中国科学院计算技术研究所 | Method and device for extracting characteristic parameter in human body waveform |
| CN102293639A (en) * | 2011-06-30 | 2011-12-28 | 芜湖圣美孚科技有限公司 | Pulse condition signal time domain feature extraction method |
| CN105286815A (en) * | 2015-11-02 | 2016-02-03 | 重庆大学 | Pulse wave signal feature point detection method based on waveform time domain features |
| GB2552455A (en) * | 2016-06-16 | 2018-01-31 | Digital & Future Tech Limited | Blood monitoring |
| CN107432736A (en) * | 2017-06-06 | 2017-12-05 | 新绎健康科技有限公司 | A kind of method for identifying pulse wave signal |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108403094A (en) * | 2018-03-23 | 2018-08-17 | 中国科学院微电子研究所 | A method of identification pulse wave wave crest |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3087915B1 (en) | Method and system for noise cleaning of photoplethysmogram signals for estimating blood pressure | |
| CN103027667B (en) | Characteristic parameter extraction of pulse wave | |
| CN105030228B (en) | The method and device of its P wave position is determined in electrocardiosignal | |
| CN105030233B (en) | A kind of electrocardiosignal ST sections of recognition methods | |
| CN109982633B (en) | Pulse detection device, image analysis device, and biological information generation system | |
| CN104138254A (en) | Non-contact type automatic heart rate measurement system and measurement method | |
| CN105266772A (en) | Measuring method for physiological parameter | |
| KR102358325B1 (en) | Method and apparatus for measuring robust pulse rate and respiratory rate using face images | |
| CN101912256A (en) | Method for processing dynamic spectral data based on single-edge extraction | |
| CN106137175A (en) | Physiological parameter is estimated | |
| CN112494001B (en) | PPG signal quality evaluation method and device, wearable device and storage medium | |
| CN107616786A (en) | Pulse wave recognition methods and pulse wave identification device | |
| Fukunishi et al. | Improvements in remote video based estimation of heart rate variability using the Welch FFT method | |
| KR101787828B1 (en) | Heartrate measuring system using skin color filter | |
| CN109009034A (en) | blood pressure measuring method, terminal and storage medium | |
| CN105046209B (en) | A kind of contactless method for measuring heart rate based on canonical correlation analysis | |
| KR101531994B1 (en) | Apparatus and method for selectively collecting electroencephalogram data through motion recognition | |
| CN104000581B (en) | ECG's data compression method and device | |
| CN104739384B (en) | Arterial wave analysis method and system | |
| CN103767690B (en) | computer analysis method and device of physiological information | |
| CN111166325B (en) | Electrocardiosignal QRS complex wave detection method and system based on IPCMM algorithm | |
| CN114073534B (en) | Heart function analysis algorithm | |
| CN110197172A (en) | One kind carrying out identity authentication method and device based on photoelectricity capacity of blood vessel information | |
| CN110037667A (en) | The human body physiological parameter monitoring method based on recognition of face applied to work station | |
| CN104688199A (en) | Non-contact type pulse measurement method based on skin pigment concentration difference |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180123 |