CN108670231A - Blood pressure measuring method, terminal and computer readable storage medium - Google Patents
Blood pressure measuring method, terminal and computer readable storage medium Download PDFInfo
- Publication number
- CN108670231A CN108670231A CN201810216746.4A CN201810216746A CN108670231A CN 108670231 A CN108670231 A CN 108670231A CN 201810216746 A CN201810216746 A CN 201810216746A CN 108670231 A CN108670231 A CN 108670231A
- Authority
- CN
- China
- Prior art keywords
- blood pressure
- matrix
- signal data
- coefficient
- value
- 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/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
Abstract
The invention discloses a kind of blood pressure measuring method, terminal and computer readable storage medium, the method includes：Acquire blood pressure data, with from blood pressure data obtain ecg signal data to pulse signal data each time difference；The calibration blood pressure of preset times is measured, and obtains the heart rate of blood pressure data, to calculate coefficient matrix according to heart rate, each time difference and each calibration blood pressure；Recursive calculation is carried out to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；Multiple coefficient matrixes are input in preset pressure value calculation formula, multiple pressure values are obtained.The present invention improves the accuracy and convenience of blood pressure measurement.
Description
Technical field
The present invention relates to a kind of medical data processing technology field more particularly to blood pressure measuring method, terminal and computers
Readable storage medium storing program for executing.
Background technology
With the development of science and technology, the technology of medical data processing also gradually develops, wherein blood pressure is a kind of important
Physiological parameter accurately measures blood pressure and has great importance for the diagnosis of doctor.
Existing frequently-used blood pressure measurement mode has cuff sphygmomanometer, and it not is very convenient to use, for this purpose, passing through list
Ppg (photoplethysmography, the arteries and veins led ECG (electrocardiogram, electrocardiosignal) signals and synchronized
Fight signal) come to measure blood pressure be very promising, the measurement for singly leading ECG and ppg is relatively easy, is easy to survey for relationship between signal
Amount.
But the existing measurement for singly leading ECG and ppg, it generally requires and measures a large amount of numbers, can just obtain accurately
As a result, losing the convenient original intention for measuring blood pressure, if pendulous frequency is few, then measurement result can be caused inaccurate.
Invention content
The main purpose of the present invention is to provide a kind of blood pressure measuring method, terminal and computer readable storage medium, purports
Existing blood pressure measurement mode is being solved, is needing pendulous frequency more, it is not convenient enough to the measurement of blood pressure, if pendulous frequency is not
It is enough, the technical problem of measurement result inaccuracy can be caused.
To achieve the above object, the present invention provides a kind of blood pressure measuring method, and the blood pressure measuring method includes：
Acquire blood pressure data, with from blood pressure data obtain ecg signal data to pulse signal data each time
Difference；
The calibration blood pressure for measuring preset times, and obtains the heart rate of blood pressure data, with according to heart rate, each time difference and each
A calibration blood pressure calculates coefficient matrix；
Recursive calculation is carried out to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；
Multiple coefficient matrixes are input in preset pressure value calculation formula, multiple pressure values are obtained.
Optionally, the step that each time difference of the ecg signal data to pulse signal data is obtained from blood pressure data
Suddenly include：
The wave crest point of ecg signal data is obtained, and obtains starting point, wave crest point and the starting point of pulse signal data
Midpoint between wave crest point；
The first time of the wave crest point of calculating ecg signal data to the starting point of pulse signal data is poor；
The wave crest point of ecg signal data is calculated to second time difference at the midpoint of pulse signal data；
Third time difference of the wave crest point of calculating ecg signal data to the wave crest point of pulse signal data.
Optionally, described the step of calculating coefficient matrix according to heart rate, each time difference and each calibration blood pressure, includes：
Constant matrices is obtained based on heart rate and each time difference；
The transposed matrix of constant matrices is obtained, and symmetrical matrix is obtained according to constant matrices and transposed matrix；
Calculate the Matrix condition number of symmetrical matrix；
According to the value of Matrix condition number, matrix parameter is determined；
According to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor matrix.
Optionally, the formula of the Matrix condition number for calculating symmetrical matrix includes：
γ(X*XT)=| | X*XT||*||X*XT-1||
Wherein, γ (X*XT) be symmetrical matrix conditional number.
Optionally, described according to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor matrix
Formula include：
Wherein, C is coefficient matrix, X*XTIt is symmetrical matrix,It is matrix parameter, I indicates unit matrix, SBPiIndicate ith
The calibration blood pressure of measurement.
Optionally, described that recursive calculation is carried out to coefficient matrix using preset filtering equations group, obtain multiple coefficient squares
Battle array the step of include：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, actual blood pressure value is obtained；
Calculate actual blood pressure value variance corresponding with the calibration error of blood pressure；
Variance, coefficient matrix, the covariance of coefficient matrix and the pressure value of calibration are input to preset filtering equations
In group, to carry out recursive calculation to coefficient matrix, multiple coefficient matrixes are obtained, wherein preset filtering equations group includes：
Wherein,
Indicate actual blood pressure value variance corresponding with the calibration error of blood pressure, SBPrIt is the pressure value of the r times calibration, C0
It is the covariance of coefficient matrix C for the C in the calculation formula of coefficient matrix, ∑ _ C.
Optionally, described to substitute into the coefficient matrix calculated in preset pressure value calculation formula, obtain practical blood
The step of pressure value includes：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, corresponding blood of each time difference is obtained
Pressure value；
In corresponding pressure value of each time difference, determines each pressure value and demarcate the error between blood pressure, and obtain
Error vector, to calculate each desired value according to error vector；
Minimization expected value is determined in each desired value, according to minimization expected value, is obtained the minimization expected value and is corresponded to
Pressure value as actual blood pressure value.
Optionally, described that multiple coefficient matrixes are input in preset pressure value calculation formula, obtain multiple pressure values
The step of include：
Multiple coefficient matrix C are input to SBP=XTIn C, multiple pressure values are obtained, wherein XTIt is the transposition of constant matrices
Matrix.
In addition, to achieve the above object, the present invention also provides a kind of terminal, the terminal include memory, processor and
It is stored in the blood pressure measurement program that can be run on the memory and on the processor, the blood pressure measurement program is described
The step of processor realizes blood pressure measuring method as described above when executing.
In addition, to achieve the above object, it is described computer-readable the present invention also provides a kind of computer readable storage medium
Blood pressure measurement program is stored on storage medium, the blood pressure measurement program realizes blood as described above when being executed by processor
The step of pressing measurement method.
Blood pressure measuring method proposed by the present invention, the measurement to blood pressure are the calibration blood pressures for first obtaining preset times, with knot
Heart rate and each time difference design factor matrix are closed, actual blood pressure value is calculated further according to the coefficient matrix, to reduce the time difference
To the error of blood pressure measurement, in addition, carrying out recursive calculation to coefficient matrix to expand out multiple coefficient squares by filtering equations group
Multiple coefficient matrixes, are subsequently directly input in preset pressure value calculation formula by battle array, you can multiple pressure values are obtained, even if
The pendulous frequency for demarcating blood pressure is less, and accurate measurement result also can be obtained, achieved the effect that convenient measurement blood pressure.
Description of the drawings
Fig. 1 is the terminal structure schematic diagram for the hardware running environment that the embodiment of the present invention is related to；
Fig. 2 is the flow diagram of blood pressure measuring method first embodiment of the present invention；
Fig. 3 is the refinement flow diagram of step S20 in Fig. 2；
Fig. 4 is the refinement flow diagram of step S30 in Fig. 2.
The embodiments will be further described with reference to the accompanying drawings for the realization, the function and the advantages of the object of the present invention.
Specific implementation mode
It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not intended to limit the present invention.
The solution of the embodiment of the present invention is mainly：Blood pressure data is acquired, to obtain electrocardiosignal from blood pressure data
Each time difference of the data to pulse signal data；The calibration blood pressure of preset times is measured, and obtains the heart rate of blood pressure data, with
Coefficient matrix is calculated according to heart rate, each time difference and each calibration blood pressure；Using preset filtering equations group to coefficient matrix
Recursive calculation is carried out, multiple coefficient matrixes are obtained；Multiple coefficient matrixes are input in preset pressure value calculation formula, are obtained
Multiple pressure values.To solve in a manner of existing blood pressure measurement, convenience is relatively low, the problem of measurement result inaccuracy.
As shown in Figure 1, the terminal structure schematic diagram for the hardware running environment that Fig. 1, which is the embodiment of the present invention, to be related to.
Terminal of the embodiment of the present invention is the medical measurement equipment for blood pressure measurement.
As shown in Figure 1, the terminal may include：Processor 1001, such as CPU, communication bus 1002, user interface
1003, network interface 1004, memory 1005.Wherein, communication bus 1002 is for realizing the connection communication between these components.
User interface 1003 may include display screen (Display), input unit such as keyboard (Keyboard), optional user interface
1003 can also include the wireline interface (such as connecting wired keyboard, wire mouse etc.) and/or wireless interface (example of standard
It is such as used to connect Wireless Keyboard, wireless mouse).Network interface 1004 may include optionally the wireline interface of standard (for connecting
Connect cable network), (such as WI-FI interfaces, blue tooth interface, infrared interface, probe interface, 3G/4G/5G networkings are logical for wireless interface
Interface etc. is believed, for connecting wireless network).Memory 1005 can be high-speed RAM memory, can also be stable memory
(non-volatile memory), such as magnetic disk storage.Memory 1005 optionally can also be independently of aforementioned processor
1001 storage device.
Optionally, terminal can also include camera, RF (Radio Frequency, radio frequency) circuit, sensor, audio
Circuit, WiFi module etc..
It will be understood by those skilled in the art that the restriction of the not structure paired terminal of terminal structure shown in Fig. 1, can wrap
It includes than illustrating more or fewer components, either combines certain components or different components arrangement.
As shown in Figure 1, as may include operating system, net in a kind of memory 1005 of computer readable storage medium
Network communication module, Subscriber Interface Module SIM and blood pressure measurement program.Wherein, operating system is that management and control terminal are provided with software
The program in source supports the operation of network communication module, Subscriber Interface Module SIM, blood pressure measurement program and other programs or software；
Network communication module is for managing and controlling network interface 1002；Subscriber Interface Module SIM is for managing and controlling user interface
1003。
In terminal shown in Fig. 1, network interface 1004 is mainly used for Connection Service device or external equipment, with server or
External equipment is into row data communication；User interface 1003 is mainly used for connecting display interface；The terminal passes through processor 1001
The blood pressure measurement program stored in memory 1005 is called, to realize following steps：
Acquire blood pressure data, with from blood pressure data obtain ecg signal data to pulse signal data each time
Difference；
The calibration blood pressure for measuring preset times, and obtains the heart rate of blood pressure data, with according to heart rate, each time difference and each
A calibration blood pressure calculates coefficient matrix；
Recursive calculation is carried out to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；
Multiple coefficient matrixes are input in preset pressure value calculation formula, multiple pressure values are obtained.
Further, the terminal calls the blood pressure measurement program stored in memory 1005 by processor 1001, with
Realize the step of each time difference of the ecg signal data to pulse signal data is obtained from blood pressure data：
The wave crest point of ecg signal data is obtained, and obtains starting point, wave crest point and the starting point of pulse signal data
Midpoint between wave crest point；
The first time of the wave crest point of calculating ecg signal data to the starting point of pulse signal data is poor；
The wave crest point of ecg signal data is calculated to second time difference at the midpoint of pulse signal data；
Third time difference of the wave crest point of calculating ecg signal data to the wave crest point of pulse signal data.
Further, the terminal calls the blood pressure measurement program stored in memory 1005 by processor 1001, with
Realize the step of coefficient matrix is calculated according to heart rate, each time difference and each calibration blood pressure：
Constant matrices is obtained based on heart rate and each time difference；
The transposed matrix of constant matrices is obtained, and symmetrical matrix is obtained according to constant matrices and transposed matrix；
Calculate the Matrix condition number of symmetrical matrix；
According to the value of Matrix condition number, matrix parameter is determined；
According to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor matrix.
Further, the formula for calculating the Matrix condition number of symmetrical matrix includes：
γ(X*XT)=| | X*XT||*||X*XT-1||
Wherein, γ (X*XT) be symmetrical matrix conditional number.
Further, described according to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor square
Battle array formula include：
Wherein, C is coefficient matrix, X*XTIt is symmetrical matrix,It is matrix parameter, I indicates unit matrix, SBPiIndicate ith
The calibration blood pressure of measurement.
Further, the terminal calls the blood pressure measurement program stored in memory 1005 by processor 1001, with
Realize the step of carrying out recursive calculation to coefficient matrix using preset filtering equations group, obtain multiple coefficient matrixes：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, actual blood pressure value is obtained；
Calculate actual blood pressure value variance corresponding with the calibration error of blood pressure；
Variance, coefficient matrix, the covariance of coefficient matrix and the pressure value of calibration are input to preset filtering equations
In group, to carry out recursive calculation to coefficient matrix, multiple coefficient matrixes are obtained, wherein preset filtering equations group includes：
Wherein,
Indicate actual blood pressure value variance corresponding with the calibration error of blood pressure, SBPrIt is the pressure value of the r times calibration, C0
It is the covariance of coefficient matrix C for the C in the calculation formula of coefficient matrix, ∑ _ C.
Further, the terminal calls the blood pressure measurement program stored in memory 1005 by processor 1001, with
The step of realizing and substitute into the coefficient matrix calculated in preset pressure value calculation formula, obtaining actual blood pressure value：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, corresponding blood of each time difference is obtained
Pressure value；
In corresponding pressure value of each time difference, determines each pressure value and demarcate the error between blood pressure, and obtain
Error vector, to calculate each desired value according to error vector；
Minimization expected value is determined in each desired value, according to minimization expected value, is obtained the minimization expected value and is corresponded to
Pressure value as actual blood pressure value.
Further, the terminal calls the blood pressure measurement program stored in memory 1005 by processor 1001, with
The step of multiple coefficient matrixes are input in preset pressure value calculation formula, obtain multiple pressure values by realization：
Multiple coefficient matrix C are input to SBP=XTIn C, multiple pressure values are obtained, wherein XTIt is the transposition of constant matrices
Matrix.
Based on above-mentioned terminal hardware structure, each embodiment of blood pressure measuring method of the present invention is proposed.
It is the flow diagram of blood pressure measuring method first embodiment of the present invention with reference to Fig. 2, Fig. 2.
In the present embodiment, the blood pressure measuring method includes：
Step S10 acquires blood pressure data, to obtain ecg signal data from blood pressure data to each of pulse signal data
A time difference；
Step S20, measures the calibration blood pressure of preset times, and obtains the heart rate of blood pressure data, with according to heart rate, Ge Geshi
Between difference and each calibration blood pressure calculating coefficient matrix；
Step S30 carries out recursive calculation to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；
Multiple coefficient matrixes are input in preset pressure value calculation formula, obtain multiple pressure values by step S40.
In the present embodiment, the blood pressure measuring method is optional is applied to terminal, and the terminal is chosen as the end described in Fig. 1
End.The terminal is measuring apparatus, and the measuring apparatus is connected with cuff, and cuff is for monitoring blood pressure, and measuring apparatus is for obtaining
The blood pressure data for taking cuff to measure.
The specific steps of blood pressure measurement are realized in this programme described in detail below：
Step S10 acquires blood pressure data, to obtain ecg signal data from blood pressure data to each of pulse signal data
A time difference；
In the present embodiment, the mode of terminal acquisition blood pressure data is chosen as：Terminal acquires blood pressure by the cuff of connection
Data, it should be appreciated that in the case that terminal connects cuff, cuff deflation curve is obtained by cuff, is getting cuff
In the case of deflation curve, cuff deflation curve removal linear components are extracted with obtaining Oscillographic Signal in Oscillographic Signal
Ecg signal data and pulse signal data, the ecg signal data and pulse signal data are denoted as blood pressure data.
Specifically, described " each time difference of the acquisition ecg signal data to pulse signal data from blood pressure data "
Step includes：
Step a, obtains the wave crest point of ecg signal data, and obtain the starting points of pulse signal data, wave crest point and
Midpoint between starting point and wave crest point；
Step b, the first time for calculating the wave crest point of ecg signal data to the starting point of pulse signal data are poor；
Step c calculates the wave crest point of ecg signal data to second time difference at the midpoint of pulse signal data；
Step d, calculate the wave crest point of ecg signal data to the wave crest point of pulse signal data the third time difference.
In the present embodiment, the wave crest point (being indicated with R points) of ecg signal data is obtained, pulse signal data are then obtained
Starting point (being indicated with Ps points), the midpoint between wave crest point (being indicated with Pp points) and starting point and wave crest point is (with Pm point tables
Show), later, the first time for calculating the wave crest point of ecg signal data to the starting point of pulse signal data is poor, that is, calculates the heart
The R points of electrical signal data and calculate the wave crest of ecg signal data to the first time poor PTT1 of the Ps points of pulse signal data
Point arrives the second time difference of the starting point of pulse signal data, i.e. the R points of calculating ecg signal data to pulse signal data
Second time difference PTT2 of Pm points, in addition, calculate the wave crest point of ecg signal data to pulse signal data wave crest point the
Three time differences, i.e. the third time difference PTT3 of the Pp points of the R points of calculating ecg signal data to pulse signal data.
Step S20, measures the calibration blood pressure of preset times, and obtains the heart rate of blood pressure data, with according to heart rate, Ge Geshi
Between difference and each calibration blood pressure calculating coefficient matrix；
After calculating each time difference of the ecg signal data to pulse signal data, the calibration of preset times is measured
Blood pressure, to initialize coefficient matrix by the calibration blood pressure of preset times, wherein the concrete numerical value of the preset times does not limit
It is fixed, it is arranged according to actual needs, it should be noted that the concrete numerical value of the preset times is less than the first definite value, to indicate to mark
Determine blood pressure pendulous frequency it is less when, can also realize that subsequent pressure value calculating process, first definite value are chosen as 10 times.
In the present embodiment, the corresponding calculation formula of heart rate (being indicated with HR) of blood pressure data includes：HR=(number
Of peaks in OMW)/T, wherein OMW indicates Oscillographic Signal.
In the calibration blood pressure for measuring preset times, and after obtaining the heart rate of blood pressure data, according to heart rate, each time difference
Coefficient matrix is calculated with each calibration blood pressure, specifically, includes with reference to Fig. 3, the step S20：
Step S21 obtains constant matrices based on heart rate and each time difference；
Step S22 obtains the transposed matrix of constant matrices, and obtains symmetrical matrix according to constant matrices and transposed matrix；
Step S23 calculates the Matrix condition number of symmetrical matrix；
Step S24 determines matrix parameter according to the value of Matrix condition number；
Step S25, according to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor matrix.
In the present embodiment, pressure value calculation formula SBP=c1*PTT+c2*HR+c3 is first introduced, definition coefficient matrix is C
=[c1, c2, c3]T, constant matrices X=[PTT, HR, 1]T, then pressure value calculation formula can be reduced to：SBP=XTC, wherein
SBP is the pressure value calculated, and X is constant matrices, XTFor the transposed matrix of constant matrices, C is coefficient matrix, needs to illustrate
It is that pressure value SBP and coefficient matrix C are values to be calculated, i.e. it is unknown that pressure value SBP and coefficient matrix C, which are values to be calculated,
Amount.
Since heart rate HR and each time difference PTT have been calculated above, constant matrices X=[PTT, HR,
1]TValue can directly obtain.After obtaining constant matrices based on heart rate and each time difference, the transposition square of constant matrices is obtained
Battle array, transposed matrix, that is, XT, after obtaining transposed matrix, symmetrical matrix is obtained according to constant matrices and transposed matrix, wherein described
Symmetrical matrix=[X*XT]。
Later, it calculates the Matrix condition number of symmetrical matrix, in the present embodiment, calculates the formula packet of the Matrix condition number of symmetrical matrix
It includes：
γ(X*XT)=| | X*XT||*||X*XT-1| |, wherein γ (X*XT) be symmetrical matrix conditional number.
After calculating the Matrix condition number of symmetrical matrix, according to the value of Matrix condition number, matrix parameter, the square are determined
Battle array parameter can be selectedIt indicates, in the present embodiment, Matrix condition number can be selected inWhen, thenγ(X*
XT) ≈ 1,It should be understood thatValue be only exemplary, can also be set as other according to actual needs
Value, does not do and repeats one by one herein.
After determining matrix parameter according to the Matrix condition number of symmetrical matrix, according to symmetrical matrix, matrix parameter, unit matrix
And each calibration blood pressure, design factor matrix, specifically, the calculation formula includes:
Wherein, C is coefficient matrix, X*XTIt is symmetrical matrix,It is matrix parameter, I indicates unit matrix, SBPiIndicate ith
The calibration blood pressure of measurement.It should be appreciated that due to symmetrical matrix X*XT, matrix parameterUnit matrix I and calibration blood pressure SBPiIt is logical
Crossing process above can all get, and therefore, coefficient matrix can be calculated according to these data.
In the present embodiment, since coefficient matrix C includes c1, c2 and c3, it is therefore desirable to could most by three equations
Coefficient matrix C is determined eventually, then three equations may respectively be C1, C2 and C3 expression, for C1, the PTT in constant matrices X
Component is PTT1, and for C2, the PTT components in constant matrices X are PTT2, and for C3, the PTT components in constant matrices X are
PTT3。
That is, when measuring calibration blood pressure for the first time, then SBPiFor SBP1, then PPT is PTT1；It measures for the first time
When to calibration blood pressure, SBPiFor SBP2, PPT PTT2；When third time measures calibration blood pressure, SBPiFor SBP3, PPT PTT3,
Three coefficient matrixes C1, C2 and C3 so can be obtained, pass through above three coefficient matrix, you can calculate final coefficient matrix
The value of C.It should be noted that the calculated C C of the calculation formula of above-mentioned coefficient matrix0It indicates.
Step S30 carries out recursive calculation to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；
In the present embodiment, after calculating coefficient matrix, coefficient matrix is carried out using preset filtering equations group
Recursive calculation, to obtain multiple coefficient matrixes, specifically, reference Fig. 4, the step S30 include：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, obtains actual blood pressure by step S31
Value；
Step S32 calculates actual blood pressure value variance corresponding with the calibration error of blood pressure；
Variance, coefficient matrix, the covariance of coefficient matrix and the pressure value of calibration are input to preset by step S33
In filtering equations group, to carry out recursive calculation to coefficient matrix, multiple coefficient matrixes are obtained, wherein preset filtering equations group
Including：
Wherein,
Indicate actual blood pressure value variance corresponding with the calibration error of blood pressure, SBPrIt is the pressure value of the r times calibration, C0
It is the covariance of coefficient matrix C for the C in the calculation formula of coefficient matrix, ∑ _ C.
In the present embodiment, first by above-mentioned calculated coefficient matrix C (C0) it is updated to pressure value calculation formula SBP=XTC
In, obtain actual blood pressure value.
Then, it calculates actual blood pressure value and demarcates the corresponding variance of the error of blood pressure, then by variance, coefficient matrix, coefficient
The covariance of matrix and the pressure value of calibration are input in preset filtering equations group, to be carried out in terms of recurrence to coefficient matrix
It calculates, obtains multiple coefficient matrixes, specific calculation formula is such as：
It should be noted that by the calculation formula, new coefficient matrix C1, C2 etc. can be calculated, it is more to obtain
A coefficient matrix, this process are the equal of being adjusted to coefficient matrix C, so that the value of coefficient matrix C is by adjusting more
Accurately, follow-up calculated pressure value also can be more accurate.
Multiple coefficient matrixes are input in preset pressure value calculation formula, you can obtain multiple blood pressures by step S40
Value.
Later, multiple coefficient matrixes are input in preset pressure value calculation formula, obtain multiple pressure values, i.e., it will be more
A coefficient matrix is input to pressure value calculation formula SBP=XTIn C, to obtain multiple pressure values, such as SBP1、SBP2、SBP3, answer
Work as understanding, XTIt is the transposed matrix of constant matrices.
It should be appreciated that by above-mentioned processing mode, as long as calculating coefficient square most starting to measure calibration blood pressure several times
After the value of battle array C, new pressure value is subsequently can be obtained, correct result can be also obtained without calibration blood pressure is measured again so that
The measurement of blood pressure is very convenient,
Blood pressure measurement mode in the prior art does not consider the heart rate of blood pressure data, the time difference is only considered, to blood pressure
Measurement result there are large errors, and when measuring blood pressure, need to measure a large amount of calibration blood pressure, can just calculate blood pressure
Value so that the measurement process of blood pressure is relatively complicated.
The blood pressure measuring method that the present embodiment proposes, the measurement to blood pressure are the calibration blood pressures for first obtaining preset times, with
In conjunction with heart rate and each time difference design factor matrix, actual blood pressure value is calculated further according to the coefficient matrix, to reduce the time
Difference is to the error of blood pressure measurement, in addition, carrying out recursive calculation to coefficient matrix to expand out multiple coefficients by filtering equations group
Multiple coefficient matrixes are subsequently directly input in preset pressure value calculation formula, you can obtain multiple pressure values, i.e., by matrix
Keep the pendulous frequency of calibration blood pressure less, accurate measurement result also can be obtained, achievees the effect that convenient measurement blood pressure.
Further, in order to improve the accuracy of pressure value calculating, blood pressure measurement of the present invention is proposed based on first embodiment
The second embodiment of method.
Difference lies in the steps for the second embodiment of blood pressure measuring method and the first embodiment of blood pressure measuring method
S31 includes：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, obtains each time difference pair by step e
The pressure value answered；
Step f determines the error between each pressure value and calibration blood pressure in corresponding pressure value of each time difference,
And error vector is obtained, to calculate each desired value according to error vector；
Step g determines minimization expected value in each desired value, according to minimization expected value, obtains the minimum expectation
It is worth corresponding pressure value as actual blood pressure value.
In the present embodiment, it is when calculating three when calculating time difference of the ecg signal data to pulse signal data
Between it is poor, therefore, the pressure value finally obtained can also obtain three pressure values, respectively SBP1, SBP2, SBP3.Obtaining three
After pressure value, first determines each pressure value and demarcate the error between blood pressure, and obtain error vector, with according to error vector
Calculate each desired value, to determine minimization expected value in each desired value, finally according to minimization expected value, obtain described in most
The corresponding pressure value of small desired value is as actual blood pressure value.I.e.：
Wherein, Ei=E (e*eT), e indicates actual blood pressure value and demarcates the error vector of blood pressure, and E indicates it is expected.
In the present embodiment, indefinite problem is defined for PTT, in the case of calculating 3 kinds of PTT, to using Kalman
3 kinds of pressure values that filtering calculates are merged, and optimal blood pressure estimated value are obtained, to improve the accuracy of blood pressure measurement.
In addition, the embodiment of the present invention also proposes a kind of computer readable storage medium, the computer readable storage medium
On be stored with blood pressure measurement program, the blood pressure measurement program realizes blood pressure measurement side as described above when being executed by processor
The step of method.
The specific implementation mode of computer readable storage medium and each embodiment of above-mentioned blood pressure measuring method are basic
Unanimously, it is not repeated herein.
It should be noted that herein, the terms "include", "comprise" or its any other variant are intended to non-row
His property includes, so that process, method, article or device including a series of elements include not only those elements, and
And further include the other elements being not explicitly listed, or further include for this process, method, article or device institute it is intrinsic
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including this
There is also other identical elements in the process of element, method, article or device.
The embodiments of the present invention are for illustration only, can not represent the quality of embodiment.
Through the above description of the embodiments, those skilled in the art can be understood that above-described embodiment side
Method can add the mode of required general hardware platform to realize by software, naturally it is also possible to by hardware, but in many cases
The former is more preferably embodiment.Based on this understanding, technical scheme of the present invention substantially in other words does the prior art
Going out the part of contribution can be expressed in the form of software products, which is stored in a storage medium
In (such as ROM/RAM, magnetic disc, CD), including some instructions are used so that a station terminal equipment (can be fixed terminal, such as object
The smart homes such as networked intelligent device, including intelligent air condition, intelligent electric lamp, intelligent power, intelligent router；It can also be movement
Many networked devices such as terminal, including smart mobile phone, wearable networking AR/VR devices, intelligent sound box, autonomous driving vehicle)
Execute the method described in each embodiment of the present invention.
It these are only the preferred embodiment of the present invention, be not intended to limit the scope of the invention, it is every to utilize this hair
Equivalent structure or equivalent flow shift made by bright specification and accompanying drawing content is applied directly or indirectly in other relevant skills
Art field, is included within the scope of the present invention.
Claims (10)
1. a kind of blood pressure measuring method, which is characterized in that the blood pressure measuring method includes：
Acquire blood pressure data, with from blood pressure data obtain ecg signal data to pulse signal data each time difference；
The calibration blood pressure of preset times is measured, and obtains the heart rate of blood pressure data, with according to heart rate, each time difference and each mark
Determine blood pressure and calculates coefficient matrix；
Recursive calculation is carried out to coefficient matrix using preset filtering equations group, obtains multiple coefficient matrixes；
Multiple coefficient matrixes are input in preset pressure value calculation formula, multiple pressure values are obtained.
2. blood pressure measuring method as described in claim 1, which is characterized in that described to obtain electrocardiosignal number from blood pressure data
Include according to the step of each time difference to pulse signal data：
The wave crest point of ecg signal data is obtained, and obtains starting point, wave crest point and the starting point and wave of pulse signal data
Midpoint between peak dot；
The first time of the wave crest point of calculating ecg signal data to the starting point of pulse signal data is poor；
The wave crest point of ecg signal data is calculated to second time difference at the midpoint of pulse signal data；
Third time difference of the wave crest point of calculating ecg signal data to the wave crest point of pulse signal data.
3. blood pressure measuring method as described in claim 1, which is characterized in that described according to heart rate, each time difference and each
Demarcating the step of blood pressure calculates coefficient matrix includes：
Constant matrices is obtained based on heart rate and each time difference；
The transposed matrix of constant matrices is obtained, and symmetrical matrix is obtained according to constant matrices and transposed matrix；
Calculate the Matrix condition number of symmetrical matrix；
According to the value of Matrix condition number, matrix parameter is determined；
According to symmetrical matrix, matrix parameter, unit matrix and each calibration blood pressure, design factor matrix.
4. blood pressure measuring method as claimed in claim 3, which is characterized in that the public affairs of the Matrix condition number for calculating symmetrical matrix
Formula includes：
γ(X*XT)=| | X*XT||*||X*XT-1||
Wherein, γ (X*XT) be symmetrical matrix conditional number.
5. blood pressure measuring method as claimed in claim 3, which is characterized in that described according to symmetrical matrix, matrix parameter, unit square
Battle array and each calibration blood pressure, the formula of design factor matrix include：
Wherein, C is coefficient matrix, X*XTIt is symmetrical matrix,It is matrix parameter, I indicates unit matrix, SBPiIndicate that ith measures
Calibration blood pressure.
6. blood pressure measuring method as described in claim 1, which is characterized in that described to use preset filtering equations group to coefficient
Matrix carries out recursive calculation, and the step of obtaining multiple coefficient matrixes includes：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, actual blood pressure value is obtained；
Calculate actual blood pressure value variance corresponding with the calibration error of blood pressure；
Variance, coefficient matrix, the covariance of coefficient matrix and the pressure value of calibration are input in preset filtering equations group,
To carry out recursive calculation to coefficient matrix, multiple coefficient matrixes are obtained, wherein preset filtering equations group includes：
Wherein,
Indicate actual blood pressure value variance corresponding with the calibration error of blood pressure, SBPrIt is the pressure value of the r times calibration, C0To be
C in the calculation formula of matrix number, ∑ _ C are the covariance of coefficient matrix C.
7. blood pressure measuring method as claimed in claim 6, which is characterized in that described to substitute into the coefficient matrix calculated in advance
If pressure value calculation formula in, the step of obtaining actual blood pressure value includes：
The coefficient matrix calculated is substituted into preset pressure value calculation formula, corresponding blood pressure of each time difference is obtained
Value；
In corresponding pressure value of each time difference, determines each pressure value and demarcate the error between blood pressure, and obtain error
Vector, to calculate each desired value according to error vector；
Minimization expected value is determined in each desired value, according to minimization expected value, obtains the corresponding blood of the minimization expected value
Pressure value is as actual blood pressure value.
8. the blood pressure measuring method as described in any one of claim 1-7, which is characterized in that described that multiple coefficient matrixes are defeated
The step of entering into preset pressure value calculation formula, obtaining multiple pressure values include：
Multiple coefficient matrix C are input to SBP=XTIn C, multiple pressure values are obtained, wherein XTIt is the transposition square of constant matrices
Battle array.
9. a kind of terminal, which is characterized in that the terminal includes memory, processor and is stored on the memory and can be
The blood pressure measurement program run on the processor is realized when the blood pressure measurement program is executed by the processor as right is wanted
The step of seeking the blood pressure measuring method described in any one of 1 to 8.
10. a kind of computer readable storage medium, which is characterized in that be stored with blood pressure survey on the computer readable storage medium
Range sequence realizes such as blood pressure measurement described in any item of the claim 1 to 8 when the blood pressure measurement program is executed by processor
The step of method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810216746.4A CN108670231B (en) | 2018-03-14 | 2018-03-14 | Blood pressure measuring method, terminal and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810216746.4A CN108670231B (en) | 2018-03-14 | 2018-03-14 | Blood pressure measuring method, terminal and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108670231A true CN108670231A (en) | 2018-10-19 |
| CN108670231B CN108670231B (en) | 2022-02-15 |
Family
ID=63799259
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810216746.4A Active CN108670231B (en) | 2018-03-14 | 2018-03-14 | Blood pressure measuring method, terminal and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108670231B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114041765A (en) * | 2021-10-20 | 2022-02-15 | 广东宝莱特医用科技股份有限公司 | Self-calibration noninvasive continuous blood pressure measuring method and device |
| CN114340483A (en) * | 2019-09-25 | 2022-04-12 | 长桑医疗（海南）有限公司 | Blood pressure calibration selection method and modeling method thereof |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5857975A (en) * | 1996-10-11 | 1999-01-12 | Dxtek, Inc. | Method and apparatus for non-invasive, cuffless continuous blood pressure determination |
| TW201021763A (en) * | 2008-09-29 | 2010-06-16 | Qualcomm Inc | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
| CN105748051A (en) * | 2016-02-18 | 2016-07-13 | 京东方科技集团股份有限公司 | Blood pressure measuring method and device |
| CN105943005A (en) * | 2016-06-01 | 2016-09-21 | 合肥芯福传感器技术有限公司 | Non-invasive blood pressure detection method based on mixing of photoelectric green-light pulses and electrocardiogram |
| CN106413534A (en) * | 2015-08-08 | 2017-02-15 | 深圳先进技术研究院 | Blood-pressure continuous-measurement device, measurement model establishment method, and system |
| CN106889979A (en) * | 2016-12-30 | 2017-06-27 | 中国科学院电子学研究所 | A kind of continuity blood pressure measuring method based on electrocardiosignal and blood oxygen volume ripple |
-
2018
- 2018-03-14 CN CN201810216746.4A patent/CN108670231B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5857975A (en) * | 1996-10-11 | 1999-01-12 | Dxtek, Inc. | Method and apparatus for non-invasive, cuffless continuous blood pressure determination |
| TW201021763A (en) * | 2008-09-29 | 2010-06-16 | Qualcomm Inc | Method and apparatus for non-invasive cuff-less blood pressure estimation using pulse arrival time and heart rate with adaptive calibration |
| CN106413534A (en) * | 2015-08-08 | 2017-02-15 | 深圳先进技术研究院 | Blood-pressure continuous-measurement device, measurement model establishment method, and system |
| CN105748051A (en) * | 2016-02-18 | 2016-07-13 | 京东方科技集团股份有限公司 | Blood pressure measuring method and device |
| CN105943005A (en) * | 2016-06-01 | 2016-09-21 | 合肥芯福传感器技术有限公司 | Non-invasive blood pressure detection method based on mixing of photoelectric green-light pulses and electrocardiogram |
| CN106889979A (en) * | 2016-12-30 | 2017-06-27 | 中国科学院电子学研究所 | A kind of continuity blood pressure measuring method based on electrocardiosignal and blood oxygen volume ripple |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114340483A (en) * | 2019-09-25 | 2022-04-12 | 长桑医疗（海南）有限公司 | Blood pressure calibration selection method and modeling method thereof |
| CN114041765A (en) * | 2021-10-20 | 2022-02-15 | 广东宝莱特医用科技股份有限公司 | Self-calibration noninvasive continuous blood pressure measuring method and device |
| CN114041765B (en) * | 2021-10-20 | 2023-11-21 | 广东宝莱特医用科技股份有限公司 | Self-calibration noninvasive continuous blood pressure measurement method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108670231B (en) | 2022-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10296835B2 (en) | Automated quality assessment of physiological signals | |
| US20180000424A1 (en) | On-demand heart rate estimation based on optical measurements | |
| US20190036545A1 (en) | Apparatus and method for data compression in a wearable device | |
| CN103168425A (en) | Method and apparatus for low complexity compression of signals | |
| EP2255583A2 (en) | Distributed spectrum sensing | |
| CN108261190B (en) | Continuous blood pressure estimation method, device and equipment | |
| WO2020098264A1 (en) | Body composition measurement method, apparatus and system, and storage medium | |
| US10959641B2 (en) | Apparatus and method for measuring bioelectrical impedance and apparatus and method for measuring biometric information | |
| CN108670231A (en) | Blood pressure measuring method, terminal and computer readable storage medium | |
| CN113647921A (en) | Blood pressure measuring method, system, device and storage medium | |
| WO2022068589A1 (en) | Dynamic heart rate prediction method and apparatus, and device and readable storage medium | |
| CN111387936A (en) | Sleep stage identification method, device and equipment | |
| WO2021034960A1 (en) | Systems and methods for imputing real-time physiological signals | |
| US20220192607A1 (en) | Apparatus and method for estimating blood pressure | |
| CN108124419A (en) | A kind of blood pressure measuring method, Intelligent heart rate earphone and system | |
| CN108567420A (en) | Blood pressure measuring method, terminal and computer readable storage medium | |
| JP2022520497A (en) | Assessment of muscle fatigue | |
| US20220233079A1 (en) | Apparatus and method for estimating blood pressure | |
| DE102019101353A1 (en) | INTELLIGENT PERSONAL PORTABLE BLOOD PRESSURE MEASUREMENT SYSTEM AND METHOD FOR CALIBRATING THE BLOOD PRESSURE MEASUREMENT IN USE | |
| CN110432882B (en) | Maximum heart rate prediction method and device based on metabolic equivalent and physiological parameters | |
| CN105326482A (en) | Device and method for recording physiological signal | |
| Ochi et al. | Non-contact Blood Pressure Estimation By Microwave Reflection Employing Machine Learning | |
| Kosmadakis et al. | Employing savitzky-golay smoothing in a low cost eHealth platform | |
| Guo et al. | Research on home healthcare management system based on the improved internet of things architecture | |
| Augustyniak | Autoadaptivity and optimization in distributed ECG interpretation |
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 |