CN108652611B - Heart rate calculation method, device, equipment and computer readable storage medium - Google Patents

Heart rate calculation method, device, equipment and computer readable storage medium Download PDF

Info

Publication number
CN108652611B
CN108652611B CN201810473882.1A CN201810473882A CN108652611B CN 108652611 B CN108652611 B CN 108652611B CN 201810473882 A CN201810473882 A CN 201810473882A CN 108652611 B CN108652611 B CN 108652611B
Authority
CN
China
Prior art keywords
heart rate
sampling
calculating
waveform
current round
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810473882.1A
Other languages
Chinese (zh)
Other versions
CN108652611A (en
Inventor
谢卓延
杨其宇
黄昭岚
蔡熊科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201810473882.1A priority Critical patent/CN108652611B/en
Publication of CN108652611A publication Critical patent/CN108652611A/en
Application granted granted Critical
Publication of CN108652611B publication Critical patent/CN108652611B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G06F2218/04—Denoising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks

Abstract

The application discloses a heart rate calculation method, which comprises the steps of obtaining an electrocardiosignal, sampling the electrocardiosignal according to a preset sampling frequency, and obtaining a sampling sequence; then, calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round; acquiring the waveform count of a preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time of the current round; and calculating the heart rate of the current round by utilizing the waveform counting and the sampling number. Therefore, the waveform counting and the sampling number are obtained according to the sampling frequency and the heart rate of the previous round, the self-adaptive time is adjusted according to the change of the heart rate, the calculation error caused by using fixed window time in the prior art is avoided, and the situation that the calculation time is too long when the time window is large is avoided, so that the accuracy and the speed of calculating the heart rate are improved. The application also discloses a heart rate calculating device, heart rate calculating equipment and a computer readable storage medium, which have the beneficial effects.

Description

Heart rate calculation method, device, equipment and computer readable storage medium
Technical Field
The present invention relates to the field of electronic information, and in particular, to a method, an apparatus, a device, and a computer-readable storage medium for calculating a heart rate.
Background
Along with the gradual improvement of health consciousness of people, the cognitive degree of each function of a body is higher and higher, the heart rate is used as one of important parameters of human physiological information, and the requirement of people on the accuracy of the heart rate is higher and higher.
Currently, the best method for measuring heart rate is to use Electrocardiogram (ECG) to measure, by detecting the R-wave position in the ECG, and then by acquiring the number of R-waves within a fixed time window, and dividing the window time by the number of R-waves, thereby obtaining the heart rate. However, since the variation range of the heart rate is large, if the selected time window is too small, the error of the calculated heart rate is large; if the selected time window is too large, the response speed is slow when the heart rate is calculated, and the heart rate cannot be calculated in time.
Therefore, how to improve the efficiency of calculating the heart rate is a technical problem that needs to be solved by those skilled in the art at present.
Disclosure of Invention
In view of this, the present invention provides a method for calculating a heart rate, which can improve efficiency of calculating a heart rate; another object of the present invention is to provide a heart rate calculating device, a heart rate calculating apparatus and a computer readable storage medium, all of which have the above advantages.
In order to solve the technical problem, the invention provides a heart rate calculation method, which comprises the following steps:
acquiring an electrocardiosignal and sampling the electrocardiosignal according to a preset sampling frequency to obtain a sampling sequence;
calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round;
acquiring a waveform count of a preset waveform and a sampling number corresponding to the waveform count according to the self-adaptive time of the current round;
and calculating the heart rate of the current round by using the waveform count and the sampling number.
Preferably, the calculating the adaptive time of the current round according to the sampling frequency and the heart rate of the previous round specifically includes:
according to
Figure BDA0001663943160000021
Calculating the self-adaptive time of the previous round;
by using
Figure BDA0001663943160000022
Calculating the self-adaptive time of the current wheel
Figure BDA0001663943160000023
Wherein m is an inequality
Figure BDA0001663943160000024
The smallest positive integer of (d); y isnRepresenting the heart rate of the current round; y isn-1Representing the heart rate of the previous round; the constant c is the maximum value of the difference between the measurement interval value and the ideal interval value of the adjacent preset waveforms; f is the sampling frequency; the constant epsilon is the allowable error in calculating the heart rate for the current round.
Preferably, the obtaining of the waveform count of the preset waveform according to the self-adaptive time of the current round and the number of samples corresponding to the waveform count specifically include:
judging whether sampling points in the sampling sequence are the preset waveforms or not, and increasing the number of the samples by 1;
if yes, further increasing the waveform count by 1;
and stopping counting after the sampling time corresponding to the sampling point is longer than the self-adaptive time of the sampling point and the sampling point is judged to be the preset waveform again, and obtaining the sampling number and the waveform count.
Preferably, after the obtaining the electrocardiographic signal and sampling the electrocardiographic signal according to a preset sampling frequency to obtain a sampling sequence, the method further includes:
and marking the preset waveform in the sampling sequence.
Preferably, further comprising:
calculating a heart rate sequence formed by combining a plurality of heart rates;
and carrying out filtering processing on the heart rate sequence.
Preferably, the preset waveform is specifically an R-wave.
In order to solve the above technical problem, the present invention further provides a heart rate calculating device, including:
the sampling module is used for acquiring the electrocardiosignals and sampling the electrocardiosignals according to a preset sampling frequency to obtain a sampling sequence;
the first calculation module is used for calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round;
the acquisition module is used for acquiring the waveform count of a preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time of the current round;
and the second calculation module is used for calculating the heart rate of the current round by utilizing the waveform count and the sampling number.
In order to solve the above technical problem, the present invention further provides a heart rate calculating device, including:
a memory for storing a computer program;
a processor for implementing the steps of any of the above methods of calculating heart rate when executing the computer program.
In order to solve the above technical problem, the present invention further provides a computer-readable storage medium, on which a computer program is stored, and the computer program, when executed by a processor, implements the steps of any one of the above methods for calculating a heart rate.
The heart rate calculation method provided by the invention comprises the following steps: acquiring an electrocardiosignal and sampling the electrocardiosignal according to a preset sampling frequency to obtain a sampling sequence; then, calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round; acquiring the waveform count of the preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time of the current round; and calculating the heart rate of the current round by using the waveform counting and the sampling number.
Therefore, the method calculates the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round, and obtains the waveform count of the preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time, namely, obtains the waveform count and the sampling number according to the sampling frequency and the heart rate of the previous round, and adjusts the self-adaptive time according to the change of the heart rate, thereby avoiding the calculation error caused by using fixed window time and the situation that the calculation time is too long when the time window is large in the prior art, and further improving the accuracy and the speed of calculating the heart rate.
In order to solve the technical problem, the invention also provides a heart rate calculating device, heart rate calculating equipment and a computer readable storage medium, which have the beneficial effects.
Drawings
In order to more clearly illustrate the embodiments or technical solutions of the present invention, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
Fig. 1 is a flowchart of a method for calculating a heart rate according to an embodiment of the present invention;
FIG. 2 is a flow chart of another method for calculating heart rate according to an embodiment of the present invention;
FIG. 3 is a block diagram of a heart rate calculating device according to an embodiment of the present invention;
fig. 4 is a block diagram of a heart rate computing device according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The heart rate calculating method has the advantages that the heart rate calculating method can improve the heart rate calculating efficiency; another core of the present invention is to provide a heart rate calculating device, a heart rate calculating apparatus, and a computer readable storage medium, all having the above beneficial effects.
In order that those skilled in the art will better understand the disclosure, the invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
Fig. 1 is a flowchart of a method for calculating a heart rate according to an embodiment of the present invention, and as shown in fig. 1, the method for calculating a heart rate includes:
s10: acquiring the electrocardiosignals, and sampling the electrocardiosignals according to a preset sampling frequency to obtain a sampling sequence.
In a specific implementation, the electrocardiograph is generally used to obtain an electrocardiograph signal, which is a device capable of automatically recording an electrocardiograph signal, which is a biological electrical signal generated by exciting the heart muscle during heart activity. Specifically, a sampling frequency is set according to actual requirements, after the electrocardiosignals are obtained, the electrocardiosignals are sampled according to the preset sampling frequency, and a sampling sequence is obtained, that is, continuous electrocardiosignals in a time domain are converted into discrete quantities to obtain the sampling sequence.
It should be noted that the acquired electrocardiographic signals may cause interference to the electrocardiographic signals due to myoelectric interference, respiratory interference, and the like in the living body, crosstalk, and the like caused by power frequency interference, and the like outside the living body, so as to be an optimal implementation manner, after the electrocardiographic signals are acquired, the acquired electrocardiographic signals may be further filtered, and then the electrocardiographic signals are sampled according to a preset sampling frequency to obtain a sampling sequence, so that the influence of interference waves on the electrocardiographic signals can be eliminated, and the accuracy of a subsequent calculation process is improved.
S20: and calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round.
Specifically, the heart rate calculated in the previous round is obtained, and the self-adaptive time of the round is calculated by using the sampling frequency and the heart rate of the previous round. It should be noted that the adaptive time is a window time for acquiring a waveform count of the preset waveform and a corresponding number of samples. Since this uncertainty of the heart rate of the previous round is used, the calculated window time will be adaptively changed according to the heart rate of the previous round.
It should be noted that, when the adaptive time of the current round is calculated, the heart rate of the current round is calculated for the first time, that is, the heart rate of the previous round is not calculated when the adaptive time is currently calculated, so that the initial heart rate can be set as the heart rate of the previous round according to actual experience or operation requirements, so as to calculate the adaptive time of the current round.
S30: and acquiring the waveform count of the preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time of the current round.
Specifically, the calculated self-adaptive time of the current round is used for obtaining the waveform count of the preset waveform and obtaining the sampling number corresponding to the waveform count. More specifically, a start time is preset, the start time may be a directly preset time point, or an end time obtained by calculating the adaptive time for multiple times at the preset time point may be used as the start time, and the waveform count and the number of samples of the current round are obtained by using the start time and the adaptive time of the current round.
S40: and calculating the heart rate of the current round by using the waveform counting and the sampling number.
Specifically, after the waveform count and the sampling number are obtained, the heart rate of the current round is calculated by using the obtained waveform count and the sampling number. More specifically, the heart rate calculation formula of the current round is as follows:
Figure BDA0001663943160000051
namely, it is
Figure BDA0001663943160000052
Wherein, ynRepresenting heart rate of the current round, m1For waveform counting of preset waveform obtained by using self-adaptive time of the present round, m2F is a preset sampling frequency in order to count the number of samples obtained by utilizing the self-adaptive time and the waveform of the current round. Note that, in the calculation process, m is1、m2And f are calculated in seconds, but in general, ynIs generally bpm, i.e. heart beat times/min, and therefore the heart rate calculation formula of the current round needs to be adjusted, i.e.:
Figure BDA0001663943160000053
according to the heart rate calculation method provided by the embodiment of the invention, the self-adaptive time of the current round is calculated according to the sampling frequency and the heart rate of the previous round, and the waveform count of the preset waveform and the sampling number corresponding to the waveform count are obtained according to the self-adaptive time, namely, the waveform count and the sampling number are obtained according to the sampling frequency and the heart rate of the previous round, and the self-adaptive time is adjusted according to the change of the heart rate, so that the calculation error caused by fixed window time in the prior art and the situation that the calculation time is too long when the time window is larger are avoided, and the accuracy and the speed of calculating the heart rate are improved.
On the basis of the above embodiment, the embodiment further explains and optimizes the technical scheme, and specifically, calculating the adaptive time of the current round according to the sampling frequency and the heart rate of the previous round specifically includes:
according to
Figure BDA0001663943160000061
Calculating the self-adaptive time of the previous round;
by using
Figure BDA0001663943160000062
Calculating the self-adaptive time of the current wheel
Figure BDA0001663943160000063
Wherein m is an inequality
Figure BDA0001663943160000064
The smallest positive integer of (d); y isnRepresenting the heart rate of the current round; y isn-1Representing the heart rate of the previous round; the constant c is the maximum value of the difference between the measurement interval value and the ideal interval value of the adjacent preset waveforms; the constant epsilon is the allowable error in calculating the heart rate for the current round.
In addition, y isnThe heart rate of the current round is a value to be calculated, and the heart rate y of the current round is calculated by calculating the self-adaptive time of the current round and acquiring the corresponding waveform counting and sampling number according to the self-adaptive time of the current roundnSpecific values of (a). The constant c is the maximum value of the difference between the measured interval value and the ideal interval value of the adjacent preset waveforms. Specifically, the ideal interval value refers to an ideal electrocardiosignal which is generated by the electrocardio simulator and meets the standard, and the ideal interval value of the preset waveform is calculated according to the preset waveform corresponding to the ideal electrocardiosignal. According to the electrocardiosignals acquired by the embodiment, the measurement interval values between the actual adjacent preset waveforms are calculated, the difference values between the ideal interval value and the measurement interval value are respectively calculated, the maximum value in the difference values is found, and then the maximum value is set as a constant c. The constant ε is the allowable error in calculating the heart rate for the current round, in particular, since the calculation is correctThe calculated heart rate has different purposes of use, so the accuracy of the calculated heart rate is required to be different, namely the allowable error of the calculated heart rate is different.
It should be noted that, because the sampling frequency is a preset fixed value, the adaptive time of the previous round is calculated by using the heart rate and the sampling frequency of the previous round, and the adaptive time satisfying the inequality is calculated first
Figure BDA0001663943160000071
The self-adaptive time of the round is calculated by using the minimum positive integer m, the sampling frequency f, the constant c and the constant epsilon, and the calculation process specifically comprises the following steps:
Figure BDA0001663943160000072
as a preferred embodiment, the present embodiment further comprises:
calculating a heart rate sequence formed by combining a plurality of heart rates;
and carrying out filtering processing on the heart rate sequence.
It should be noted that, when the next round of calculation needs to be performed, the calculated heart rate of the current round is set as the heart rate of the previous round, so that the adaptive time of the next round is calculated by using the heart rate of the previous round and the sampling frequency, and the next round of calculation can be performed, that is, the calculated y is usednIs set to yn-1(ii) a And then combining the plurality of calculated heart rates into a heart rate sequence, and carrying out filtering processing on the heart rate sequence by using a preset filtering method.
In a preferred embodiment, after the plurality of heart rates are calculated and combined into a heart rate sequence, the heart rate sequence is subjected to median filtering, so that the calculated heart rate value is more stable and accurate.
Fig. 2 is a flowchart of another heart rate calculation method according to an embodiment of the present invention. As shown in fig. 2, the waveform count of the preset waveform and the sampling number corresponding to the waveform count obtained according to the self-adaptive time of the current round are specifically:
s31: judging whether sampling points in the sampling sequence are preset waveforms or not, and increasing the number of the samples by 1;
s32: if yes, further increasing the waveform count by 1;
s33: and stopping counting when the sampling time corresponding to the sampling point is longer than the self-adaptive time of the round and the sampling point is judged to be the preset waveform again, so as to obtain the sampling number and the waveform count.
Generally, whether a waveform corresponding to a sampling point in a sampling sequence is a preset waveform is sequentially judged, and a count value of the number of samples is increased by 1 after each judgment; and if the waveform corresponding to the sampling point in the sampling sequence is judged to be the preset waveform, increasing the waveform count by 1. It can be understood that when the sampling point in the sampling sequence is determined not to be the preset waveform, the number of samples is increased by 1, but S32 does not need to be performed, that is, the waveform count does not need to be increased by 1, and then it is continuously determined whether the next sampling point is the preset waveform. And when the sampling time corresponding to the sampling point is longer than the self-adaptive time of the current round and the waveform corresponding to the current sampling point is a preset waveform, stopping counting, and thus obtaining the sampling number and the waveform count. In addition, when the next round of calculation is performed, the waveform counting and the sampling number of the preset waveform need to be restarted to be counted; if the count is the accumulated count, when the heart rate of the next round is calculated, the count of the current round needs to be subtracted from the accumulated count, that is, the count of the next round is obtained.
As a preferred embodiment, after acquiring the cardiac electrical signal and sampling the cardiac electrical signal according to a preset sampling frequency to obtain a sampling sequence, the method further includes:
and marking a preset waveform in the sampling sequence.
Specifically, after obtaining the electrocardiographic signal and sampling the electrocardiographic signal according to a preset sampling frequency to obtain a sampling sequence, a preset waveform in the sampling sequence may be marked, specifically, the preset waveform may be marked by a specific identifier, or a sampling point in the sampling sequence corresponding to the preset waveform may also be marked, which is not limited in this embodiment. It can be understood that whether the waveform corresponding to the sampling point is the preset waveform can be judged more conveniently by marking the preset waveform in the sampling sequence.
As a preferred embodiment, the preset waveform is specifically an R-wave. It will be appreciated that a typical wave set of an electrocardiographic signal includes a P-wave, a QRS complex, a T-wave, and a U-wave. Wherein, the QRS complex comprises three closely linked waves: the first downward wave is called Q wave, a high-tip vertical wave following the Q wave is called R wave, and a downward wave following the R wave is called S wave, that is, the waveform of the R wave has a characteristic remarkably different from other waveforms, is not easily interfered, and is not easily filtered relative to other waves when filtering is performed. That is to say, use R wave as presetting the waveform and calculate, it is more convenient accurate.
The above detailed description is given for the embodiment of the method for calculating a heart rate according to the present invention, and the present invention further provides a device, an apparatus, and a computer-readable storage medium for calculating a heart rate corresponding to the method.
Fig. 3 is a block diagram of a heart rate calculating apparatus according to an embodiment of the present invention, as shown in fig. 3, including:
the sampling module 31 is configured to acquire an electrocardiographic signal and sample the electrocardiographic signal according to a preset sampling frequency to obtain a sampling sequence;
the first calculation module 32 is used for calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round;
the obtaining module 33 is configured to obtain a waveform count of a preset waveform and a sampling number corresponding to the waveform count according to the self-adaptive time of the current round;
and the second calculating module 34 is used for calculating the heart rate of the current round by utilizing the waveform counting and the sampling number.
As a preferred embodiment, the first calculation module 32 specifically includes:
first computation submoduleA block for according to
Figure BDA0001663943160000091
Calculating the self-adaptive time of the previous round;
a second calculation submodule for utilizing
Figure BDA0001663943160000092
Calculating the self-adaptive time of the current wheel
Figure BDA0001663943160000093
Wherein m is an inequality
Figure BDA0001663943160000094
The smallest positive integer of (d); y isnRepresenting the heart rate of the current round; y isn-1Representing the heart rate of the previous round; the constant c is the maximum value of the difference between the measurement interval value and the ideal interval value of the adjacent preset waveforms; the constant epsilon is the allowable error in calculating the heart rate for the current round.
As a preferred embodiment, the obtaining module 33 specifically includes:
the judgment submodule is used for judging whether the sampling points in the sampling sequence are preset waveforms or not and increasing the number of the samples by 1;
the first counting submodule is used for increasing the waveform counting by 1 when the sampling point in the sampling sequence is judged to be a preset waveform;
and the second counting submodule is used for stopping counting after the sampling time corresponding to the sampling point is longer than the self-adaptive time of the round and the sampling point is judged to be the preset waveform again, so that the sampling number and the waveform count are obtained.
The heart rate calculating device provided by the invention has the beneficial effects of the heart rate calculating method.
Fig. 4 is a block diagram of a heart rate computing device according to an embodiment of the present invention, as shown in fig. 4, including:
a memory 41 for storing a computer program;
a processor 42 for implementing the steps of the above-described method of calculating a heart rate when executing a computer program.
The heart rate calculating device has the beneficial effects of the heart rate calculating method.
In order to solve the above technical problem, the present invention further provides a computer readable storage medium, on which a computer program is stored, and the computer program, when executed by a processor, implements the steps of the method for calculating a heart rate.
The computer-readable storage medium provided by the embodiment has the beneficial effects of the method for calculating the heart rate.
The heart rate calculating method, the heart rate calculating device, the heart rate calculating equipment and the heart rate calculating computer-readable storage medium are described in detail above. The principles and embodiments of the present invention are explained herein using specific examples, which are set forth only to help understand the method and its core ideas of the present invention. It should be noted that, for those skilled in the art, it is possible to make various improvements and modifications to the present invention without departing from the principle of the present invention, and those improvements and modifications also fall within the scope of the claims of the present invention.
The embodiments are described in a progressive manner in the specification, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative components and steps have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may reside in Random Access Memory (RAM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

Claims (8)

1. A method of calculating heart rate, comprising:
acquiring an electrocardiosignal and sampling the electrocardiosignal according to a preset sampling frequency to obtain a sampling sequence;
calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round;
acquiring a waveform count of a preset waveform and a sampling number corresponding to the waveform count according to the self-adaptive time of the current round;
calculating the heart rate of the current round by using the waveform count and the sampling number;
wherein, the calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round specifically comprises:
according to
Figure FDA0002787948770000011
Calculating the self-adaptive time of the previous round;
by using
Figure FDA0002787948770000012
Calculating the self-adaptive time of the current wheel
Figure FDA0002787948770000013
Wherein m is an inequality
Figure FDA0002787948770000014
The smallest positive integer of (d); y isnRepresenting the heart rate of the current round; y isn-1Representing the heart rate of the previous round; the constant c is the maximum value of the difference between the measurement interval value and the ideal interval value of the adjacent preset waveforms; f is the sampling frequency; the constant epsilon is an allowable error for calculating the heart rate of the current round;
wherein the calculating the heart rate of the current round using the waveform count and the number of samples comprises:
by using
Figure FDA0002787948770000015
Calculating the heart rate of the current round;
wherein, ynRepresenting the heart rate of the current round, m1The waveform count, m, of the preset waveform obtained using the self-adaptation time of the current round2F is the preset sampling frequency, and the sampling number is acquired by utilizing the self-adaptive time of the current round and the waveform counting.
2. The calculation method according to claim 1, wherein the obtaining of the waveform count of the preset waveform according to the self-adaptive time of the current round and the number of samples corresponding to the waveform count specifically include:
judging whether sampling points in the sampling sequence are the preset waveforms or not, and increasing the number of the samples by 1;
if yes, further increasing the waveform count by 1;
and stopping counting after the sampling time corresponding to the sampling point is longer than the self-adaptive time of the sampling point and the sampling point is judged to be the preset waveform again, and obtaining the sampling number and the waveform count.
3. The computing method according to claim 1, wherein after obtaining the cardiac electrical signal and sampling the cardiac electrical signal according to a preset sampling frequency to obtain a sampling sequence, the method further comprises:
and marking the preset waveform in the sampling sequence.
4. The computing method of claim 1, further comprising:
calculating a heart rate sequence formed by combining a plurality of heart rates;
and carrying out filtering processing on the heart rate sequence.
5. Calculation method according to any one of claims 1 to 4, characterized in that said preset waveform is in particular an R-wave.
6. A heart rate calculation device, comprising:
the sampling module is used for acquiring the electrocardiosignals and sampling the electrocardiosignals according to a preset sampling frequency to obtain a sampling sequence;
the first calculation module is used for calculating the self-adaptive time of the current round according to the sampling frequency and the heart rate of the previous round;
the acquisition module is used for acquiring the waveform count of a preset waveform and the sampling number corresponding to the waveform count according to the self-adaptive time of the current round;
the second calculation module is used for calculating the heart rate of the current round by utilizing the waveform count and the sampling number;
wherein, the first calculating module specifically comprises:
a first calculation submodule for calculating based on
Figure FDA0002787948770000021
Calculating the self-adaptive time of the previous round;
a second calculation submodule for utilizing
Figure FDA0002787948770000022
Calculating the self-adaptive time of the current wheel
Figure FDA0002787948770000023
Wherein m is an inequality
Figure FDA0002787948770000024
The smallest positive integer of (d); y isnRepresenting the heart rate of the current round; y isn-1Representing the heart rate of the previous round; the constant c is the maximum value of the difference between the measurement interval value and the ideal interval value of the adjacent preset waveforms; constant epsilon is the allowable error for calculating the heart rate of the current round
Wherein the calculating the heart rate of the current round using the waveform count and the number of samples comprises:
by using
Figure FDA0002787948770000031
Calculating the heart rate of the current round;
wherein, ynRepresenting the heart rate of the current round, m1The waveform count, m, of the preset waveform obtained using the self-adaptation time of the current round2F is the preset sampling frequency, and the sampling number is acquired by utilizing the self-adaptive time of the current round and the waveform counting.
7. A heart rate computing device, comprising:
a memory for storing a computer program;
a processor for implementing the steps of the method of calculating a heart rate according to any one of claims 1 to 5 when executing the computer program.
8. A computer-readable storage medium, characterized in that a computer program is stored on the computer-readable storage medium, which computer program, when being executed by a processor, carries out the steps of the method of calculating a heart rate according to any one of claims 1 to 5.
CN201810473882.1A 2018-05-17 2018-05-17 Heart rate calculation method, device, equipment and computer readable storage medium Active CN108652611B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810473882.1A CN108652611B (en) 2018-05-17 2018-05-17 Heart rate calculation method, device, equipment and computer readable storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810473882.1A CN108652611B (en) 2018-05-17 2018-05-17 Heart rate calculation method, device, equipment and computer readable storage medium

Publications (2)

Publication Number Publication Date
CN108652611A CN108652611A (en) 2018-10-16
CN108652611B true CN108652611B (en) 2021-03-16

Family

ID=63776126

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810473882.1A Active CN108652611B (en) 2018-05-17 2018-05-17 Heart rate calculation method, device, equipment and computer readable storage medium

Country Status (1)

Country Link
CN (1) CN108652611B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109635786B (en) * 2019-01-22 2022-05-17 佛山市百步梯医疗科技有限公司 Heart rate calculation method based on heart sound autocorrelation function
CN111657905B (en) * 2020-06-23 2022-06-07 中国医学科学院生物医学工程研究所 Feature point detection method, device, equipment and storage medium
CN112057067A (en) * 2020-09-03 2020-12-11 青岛歌尔智能传感器有限公司 Heart rate detection method, wearable device and readable storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5357970A (en) * 1993-04-08 1994-10-25 Critikon, Inc. Method for determining dominant heart rates
CN103153175A (en) * 2010-05-14 2013-06-12 先进计算发展中心 Diagnostic classifications of pulse signal waveform data
CN105326477A (en) * 2015-10-30 2016-02-17 复旦大学 Synchronous breath and heartbeat signal extraction method and synchronous breath and heartbeat signal extraction device for collecting multiple medical signals
CN106974637A (en) * 2017-04-11 2017-07-25 西南大学 Heart rate real-time computing technique and device
CN107569227A (en) * 2017-10-26 2018-01-12 太仓市同维电子有限公司 The processing method and monitoring device of heart rate under a kind of motion state

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5357970A (en) * 1993-04-08 1994-10-25 Critikon, Inc. Method for determining dominant heart rates
CN103153175A (en) * 2010-05-14 2013-06-12 先进计算发展中心 Diagnostic classifications of pulse signal waveform data
CN105326477A (en) * 2015-10-30 2016-02-17 复旦大学 Synchronous breath and heartbeat signal extraction method and synchronous breath and heartbeat signal extraction device for collecting multiple medical signals
CN106974637A (en) * 2017-04-11 2017-07-25 西南大学 Heart rate real-time computing technique and device
CN107569227A (en) * 2017-10-26 2018-01-12 太仓市同维电子有限公司 The processing method and monitoring device of heart rate under a kind of motion state

Also Published As

Publication number Publication date
CN108652611A (en) 2018-10-16

Similar Documents

Publication Publication Date Title
Mahdiani et al. Is 50 Hz high enough ECG sampling frequency for accurate HRV analysis?
JP6219942B2 (en) Real-time QRS period measurement in ECG
US10973465B2 (en) Pacing signal processing method, system and electrocardiogram monitor
CN108652611B (en) Heart rate calculation method, device, equipment and computer readable storage medium
US11253204B2 (en) Method for assessing electrocardiogram signal quality
US5511554A (en) Real-time artifact removal from waveforms using a dynamic filter having a fixed delay
KR101910982B1 (en) Method and apparatus for eliminating motion artifact of biosignal using personalized biosignal pattern
EP0434856A1 (en) Method of deriving a respiration signal and/or a cardiac artifact signal from a physiological signal
CN110505839B (en) Method and system for processing EMG signals
US20190365268A1 (en) Method and apparatus for processing ecg signals, and ecg machine
CN110507317B (en) Self-adaptive CA-CFAR (Carrier frequency-constant false alarm) positioning method for electrocardiosignal R wave
US9706945B2 (en) Respiration rate determination in impedance pneumography
McSharry et al. Open-source software for generating electrocardiogram signals
Berset et al. Robust heart rhythm calculation and respiration rate estimation in ambulatory ECG monitoring
CN111419219A (en) PPG heart beat signal preprocessing method and device and atrial fibrillation detection equipment
Vega-Martínez et al. Wavelet packet based algorithm for QRS region detection and R/S wave identification
WO1984003032A1 (en) System and method for predicting ventricular tachycardia
US9826937B2 (en) Method and apparatus for reducing motion artifacts in ECG signals
JPH05212006A (en) Heartbeat interval measuring instrument
KR20090040721A (en) Method for detecting heart-rate
WO2018167854A1 (en) Analysis apparatus and analysis program
Brown et al. Real-time estimation of the ECG-derived respiration (EDR) signal
Yadav et al. Time and frequency exploration of ECG signal
Ahmad et al. An improved QRS complex detection for online medical diagnosis
Rashkovska et al. Signal processing methods for ST variability assessment in ECG

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