CN113476039B - Ballistocardiogram signal acquisition method and device, storage medium and computer equipment - Google Patents

Ballistocardiogram signal acquisition method and device, storage medium and computer equipment Download PDF

Info

Publication number
CN113476039B
CN113476039B CN202110873353.2A CN202110873353A CN113476039B CN 113476039 B CN113476039 B CN 113476039B CN 202110873353 A CN202110873353 A CN 202110873353A CN 113476039 B CN113476039 B CN 113476039B
Authority
CN
China
Prior art keywords
signal
ballistocardiogram
signals
time window
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110873353.2A
Other languages
Chinese (zh)
Other versions
CN113476039A (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.)
Shenzhen Shuliantianxia Intelligent Technology Co Ltd
Original Assignee
Shenzhen Shuliantianxia Intelligent Technology Co Ltd
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 Shenzhen Shuliantianxia Intelligent Technology Co Ltd filed Critical Shenzhen Shuliantianxia Intelligent Technology Co Ltd
Priority to CN202110873353.2A priority Critical patent/CN113476039B/en
Publication of CN113476039A publication Critical patent/CN113476039A/en
Application granted granted Critical
Publication of CN113476039B publication Critical patent/CN113476039B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1102Ballistocardiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6887Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
    • A61B5/6892Mats
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0247Pressure sensors

Abstract

The invention discloses a method, a device, a storage medium and computer equipment for acquiring a ballistocardiogram signal, wherein the method comprises the following steps: at least three piezoelectric sensors are arranged at intervals in the length direction of the mattress to collect multi-channel ballistocardiogram signals, and a preset time window is slid on the multi-channel ballistocardiogram signals to screen out optimal signal data in the multi-channel ballistocardiogram signals, and the output ballistocardiogram signals are updated based on the ballistocardiogram signals of a data channel corresponding to the optimal signal data, so that even if the sleeping posture of a human body is changed continuously, the ballistocardiogram signals collected based on the optimal data channel can be used for updating the output ballistocardiogram signals, the strength of the collected ballistocardiogram signals is improved, and the physiological information of the sleeping human body is analyzed conveniently.

Description

Ballistocardiogram signal acquisition method and device, storage medium and computer equipment
Technical Field
The present invention relates to the field of signal processing technologies, and in particular, to a method and apparatus for acquiring a ballistocardiogram signal, a storage medium, and a computer device.
Background
Ballistocardiogram (BCG) is derived from the fact that heart pumping causes blood to flow in large blood vessels, creating an impact force on a supporting object in close contact with the human body. Physiological information related to sleep, such as heart rate/respiration, etc., can be calculated from the ballistocardiogram. Heart rate variability data may then be obtained to assess sleep quality, stress conditions, or cardiac function, etc.
In the process of acquiring the ballistocardiogram signal, the existing piezoelectric sensor can cause the intensity of the ballistocardiogram signal acquired by the piezoelectric sensor to be weaker and even can not acquire the ballistocardiogram signal because the sleeping posture of a human body in sleeping is not controlled.
Disclosure of Invention
In view of the foregoing, it is necessary to provide a method, an apparatus, a storage medium, and a computer device for acquiring a ballistocardiogram signal.
The first aspect of the invention provides a method for acquiring a ballistocardiogram signal, which comprises the following steps:
acquiring multipath heart attack signals acquired by all piezoelectric sensors on a mattress, wherein the mattress at least comprises three piezoelectric sensors arranged along the length direction of the mattress; respectively sliding a preset time window on the plurality of paths of ballistocardiogram signals by taking a preset first time length as a step length to obtain signal data of the plurality of paths of ballistocardiogram signals in the preset time window; and determining a data channel according to the optimal signal data in the plurality of paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
Wherein, the determining a data channel according to the optimal signal data in the multiple paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel comprises: acquiring optimal signal data in multiple paths of signal data, wherein the optimal signal data are used for indicating that the energy values of the ballistocardiogram signals acquired by all piezoelectric sensors are the largest in the same time period; and acquiring a data channel corresponding to the optimal signal data, and using the ballistocardiogram signal of the last step length on the data channel to update the output ballistocardiogram signal.
The sliding a preset time window on the multiple ballistocardiogram signals respectively by taking a preset first time length as a step length to obtain signal data of the multiple ballistocardiogram signals in the preset time window comprises the following steps: taking a preset first time length as a step length, and taking out the heart attack signal in the preset time window backwards; and carrying out iterative computation on the ballistocardiogram signals in the preset time window to obtain signal data in the preset time window.
Wherein, after iteratively calculating the ballistocardiogram signal in the preset time window to obtain the signal data in the preset time window, the method further comprises: judging whether the signal data in each preset time window are smaller than a first preset threshold value or not, wherein the first preset threshold value is the lower limit value of the heart attack signal of a normal user; if yes, the ballistocardiogram signal of the data channel determined by the last sliding time window is used for updating the output ballistocardiogram signal.
Wherein, after iteratively calculating the ballistocardiogram signal in the preset time window to obtain the signal data in the preset time window, the method further comprises: judging whether the signal data in each preset time window are larger than a second preset threshold value or not, wherein the second preset threshold value is the upper limit value of the ballistocardiogram signal of a normal user; if yes, the ballistocardiogram signal of the data channel determined by the last sliding time window is used for updating the output ballistocardiogram signal.
After the multi-channel ballistocardiographic signals acquired by all piezoelectric sensors on the mattress are acquired, the method further comprises the following steps: respectively carrying out power frequency interference filtering treatment on the multipath ballistocardiogram signals; respectively denoising the multipath heart attack signals subjected to the power frequency interference filtering treatment; by utilizing a first-order difference method, J points corresponding to central beats of the de-noised multipath heart attack signals are highlighted, and the obtained multipath heart attack signals are favorable for calculating heart rate signals of a human body.
After the multi-channel ballistocardiographic signals acquired by all piezoelectric sensors on the mattress are acquired, the method further comprises the following steps: respectively carrying out power frequency interference filtering treatment on the multipath ballistocardiogram signals; denoising the multipath heart attack signals subjected to the power frequency interference filtering treatment; and a sliding tie filtering method of 5 points is utilized to carry out smooth denoising on the multipath ballistocardiogram signals after denoising treatment, and the obtained multipath ballistocardiogram signals are beneficial to calculating human respiratory signals.
A second aspect of the present invention provides a ballistocardiographic signal acquisition device, the device comprising:
a mattress; at least three piezoelectric sensors arranged on the mattress in a arrayed manner along the length direction of the mattress; the host unit is electrically connected with all the piezoelectric sensors and is used for converting voltage signals of the piezoelectric sensors into ballistocardiographic signals; a signal acquisition subunit, a data acquisition subunit and a signal update subunit; the signal acquisition subunit is used for acquiring multipath heart attack signals acquired by all piezoelectric sensors on the mattress; the data acquisition subunit is used for respectively sliding a preset time window on the plurality of paths of ballistocardiographic signals by taking a preset first time length as a step length so as to obtain signal data of the plurality of paths of ballistocardiographic signals in the preset time window; the signal updating subunit is used for determining a data channel according to the optimal signal data in the plurality of paths of signal data and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
A third aspect of the invention provides a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of:
acquiring multipath heart attack signals acquired by all piezoelectric sensors on a mattress, wherein the mattress at least comprises three piezoelectric sensors arranged along the length direction of the mattress; respectively sliding a preset time window on the plurality of paths of ballistocardiogram signals by taking a preset first time length as a step length to obtain signal data of the plurality of paths of ballistocardiogram signals in the preset time window; and determining a data channel according to the optimal signal data in the plurality of paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
A fourth aspect of the invention provides a computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the method of any one of the preceding claims.
The embodiment of the invention has the following beneficial effects: at least three piezoelectric sensors are arranged at intervals in the length direction of the mattress to collect multi-channel ballistocardiogram signals, and a preset time window is slid on the multi-channel ballistocardiogram signals to screen out optimal signal data in the multi-channel ballistocardiogram signals, and the output ballistocardiogram signals are updated based on the ballistocardiogram signals of a data channel corresponding to the optimal signal data, so that even if the sleeping posture of a human body is changed continuously, the ballistocardiogram signals collected based on the optimal data channel can be used for updating the output ballistocardiogram signals, the strength of the collected ballistocardiogram signals is improved, and the physiological information of the sleeping human body is analyzed conveniently.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Wherein:
FIG. 1 is a flow chart of a method of acquiring ballistocardiogram signals according to one embodiment;
FIG. 2 is a flow chart of a method of acquiring ballistocardiographic signals that facilitates calculation of heart rate according to one embodiment;
FIG. 3 is a schematic diagram of voltage signals and ballistocardiogram signals collected by three piezoelectric sensors according to an embodiment of a method for collecting ballistocardiogram signals that is useful for calculating heart rate;
FIG. 4 is a flow chart of a method of acquisition of ballistocardiographic signals useful for computing respiration in one embodiment;
FIG. 5 is a schematic diagram of voltage signals and ballistocardiographic signals acquired by three piezoelectric sensors according to an embodiment of an acquisition method for calculating ballistocardiographic signals for respiration;
FIG. 6 is a block diagram of a center impact signaling device according to one embodiment;
FIG. 7 is a block diagram of a computer device in one embodiment.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1, the present embodiment provides a method for acquiring a ballistocardiogram signal, which includes:
s101, acquiring multipath ballistocardiogram signals acquired by all piezoelectric sensors on a mattress;
s102, respectively sliding a preset time window on the multipath ballistocardiogram signals by taking a preset first time length as a step length to obtain signal data of the multipath ballistocardiogram signals in the preset time window;
s103, determining a data channel according to the optimal signal data in the multipath signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
In step S101, the mattress includes at least three piezoelectric sensors arranged along the length direction of the mattress, and the collecting area formed by all the piezoelectric sensors completely covers the heart center position and the abdomen center position of the human body when the human body is laid on the mattress. The heart impact signals acquired from the heart center position and the abdomen center position are closest to real data, so that the distances among all piezoelectric sensors completely cover the heart center position and the abdomen center position of a human body, and the heart impact signals can be acquired more accurately.
In step S102, after the preset time window is slid on the multiple cardiac impulse signals, signal data of the multiple cardiac impulse signals within the preset time window may be obtained, and specifically, the signal data may be data for calculating a heart rate or data for calculating respiration.
According to the method for acquiring the ballistocardiogram signals, at least three piezoelectric sensors are arranged at intervals in the length direction of the mattress to acquire the multipath ballistocardiogram signals, the optimal signal data in the multipath ballistocardiogram signals are screened out by sliding a preset time window on the multipath ballistocardiogram signals, and the output ballistocardiogram signals are updated based on the ballistocardiogram signals of the data channels corresponding to the optimal signal data, so that even if the sleeping posture of a human body is changed continuously, the ballistocardiogram signals acquired based on the optimal data channels can be still used for updating the output ballistocardiogram signals, the strength of the acquired ballistocardiogram signals is improved, and the physiological information of the human body during sleeping is analyzed conveniently.
In one embodiment, step S102, with a preset first time length as a step, slides a preset time window on the multiple ballistocardiogram signals respectively, so as to obtain signal data of the multiple ballistocardiogram signals in the preset time window, including:
s1021, taking a preset first time length as a step length, and taking out a heart attack signal in a preset time window backwards;
and S1022, iteratively calculating the ballistocardiogram signals in the preset time window to obtain the signal data in the preset time window.
In step S1021, the first duration is a preset duration, which may be 1 second or may be 2 seconds or less than the duration of the time window, and in this embodiment, the first duration is 1S, and the time window is slid by using a step size of 1S, so as to take out the ballistocardiogram in the time window.
In step S1022, iterative calculation is to calculate the ballistocardiogram signal in the time window after each preset time window slides, for example, the duration of the preset time window is 5S, the ballistocardiogram signal of 0-5 seconds is framed for the first time, 1S is used as step length, the ballistocardiogram signal of 1-6S is framed for the second time, and at this time, the ballistocardiogram signals of the first and second time frames of the time window are iteratively calculated to obtain signal data of two frame selections; in addition, in this embodiment, taking 3 piezoelectric sensors as an example, when the time window is framed for the first time and the second time, the 3-path ballistocardiogram signals are actually framed on the 3-path ballistocardiogram signals, and the 3-path ballistocardiogram signals framed for the first time and the 3-path ballistocardiogram signals framed for the second time are obtained respectively.
In one embodiment, in step S1022, after iteratively calculating the ballistocardiogram signal in the preset time window to obtain the signal data in the preset time window, the ballistocardiogram signal acquisition method further includes:
s1023, judging whether signal data in each preset time window are smaller than a first preset threshold value, wherein the first preset threshold value is the lower limit value of a heart attack signal of a normal user;
and S1024, if so, using the ballistocardiogram signal of the data channel determined by the last sliding time window for updating the output ballistocardiogram signal.
In step S1023, since the first preset threshold is the lower limit value of the ballistocardiogram signal of the normal user, if the signal data in each preset time window is smaller than the first preset value, it indicates that no one is on the mattress at this time, and the ballistocardiogram signal acquired by the multi-channel ballistocardiogram signal is wrong, so that the ballistocardiogram signal is closest to the real data, step S1024 is used to update the ballistocardiogram signal of the data channel determined by the last sliding time window.
In one embodiment, in step S1022, after iteratively calculating the ballistocardiogram signal in the preset time window to obtain the signal data in the preset time window, the ballistocardiogram signal acquisition method further includes:
s1025, judging whether the signal data in each preset time window are larger than a second preset threshold value, wherein the second preset threshold value is the upper limit value of the heart attack signal of a normal user;
and S1026, if yes, using the ballistocardiogram signal of the data channel determined by the last sliding time window for updating the output ballistocardiogram signal.
In step S1025, since the second preset threshold is the upper limit value of the ballistocardiogram signal of the normal user, if the signal data in each preset time window is greater than the second preset value, it indicates that the pressure on the mattress is not brought by a person at this time, and therefore the ballistocardiogram signal acquired by the multi-channel ballistocardiogram signal is erroneous, and in order to make the ballistocardiogram signal closest to the real data, step S1026 is used to update the ballistocardiogram signal of the data channel determined by the last sliding time window.
In one embodiment, step S103, determining a data channel according to the optimal signal data in the multiple signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel, includes:
s1031, obtaining optimal signal data in the multipath signal data;
s1032, acquiring a data channel corresponding to the optimal signal data, and using the ballistocardiogram signal of the last step length on the data channel to update the output ballistocardiogram signal.
In step S1031, since multiple ballistocardiographic signals are collected, after the sliding time window, a signal data is generated in the data channel where each ballistocardiographic signal is located, and the optimal signal data is used for indicating the signal data with the largest energy value of the ballistocardiographic signals collected in all piezoelectric sensors in the same time period.
In addition, the same period in the above refers to a period within the same time window.
The energy value may be a heart rate energy value or a respiratory energy value, and in particular, different energy values may be obtained according to different settings.
In step S1032, the data channel where the optimal energy data is located is obtained, and the last step of the ballistocardiogram signal is updated to output the ballistocardiogram signal, so that the most recently generated ballistocardiogram signal with one step can be ensured to be the signal closest to the real situation, and thus the error between the most recently generated ballistocardiogram signal and the actual data is reduced.
In one embodiment, step S101, after acquiring the multiple ballistocardiogram signals acquired by all the piezoelectric sensors on the mattress, the method for acquiring the ballistocardiogram signals further includes:
s111, respectively performing power frequency interference filtering treatment on the multipath ballistocardiogram signals;
s121, respectively denoising the multipath heart attack signals subjected to the power frequency interference filtering treatment;
s131, using a first-order difference method to highlight J points corresponding to central beats of the de-noised multipath heart attack signals, and the multipath heart attack signals are beneficial to calculating heart rate signals of a human body.
In this embodiment, in step S111, the frequency interference filtering process may be performed by using a trap, where the parameter of the trap is 50Hz, in step S121, the denoising process may be performed by using a butterworth band-pass filter, where the parameter of the butterworth band-pass filter is 0.67-10Hz, and after steps S111 and S121, the frequency interference and the noise of the ballistocardiogram signal are removed to a certain extent, so that the voltage signal is converted into the ballistocardiogram signal, thereby improving the accuracy of the voltage signal, and the J point corresponding to the ballistocardiogram is the data at the ballistocardiogram signal of the ballistocardiogram, so that using step S131, the voltage signal can be converted into the ballistocardiogram signal for calculating the heart rate data.
In this embodiment, after step S131 is performed, heart beat data of the human body can be obtained by calculating the heart attack signal; and in the corresponding subsequent step S1023, the first preset threshold is a heart rate lower limit threshold, and in the corresponding step S1025, the second preset threshold is a heart rate upper limit threshold.
Specifically, referring to fig. 2, in combination with the above embodiment, when acquiring the ballistocardiogram signal that is beneficial to calculating the heart rate, the method for acquiring the ballistocardiogram signal is as follows:
step S051, step 1 second, iteratively calculating heart beat energy values in a 5 second time window.
Step S052, judging whether the current data time length is more than or equal to 5 seconds. If yes, step S054 is entered, and if no, step S053 is entered.
Step S053, step length is 1 second, heart beat energy in a time window of 1 second is calculated iteratively, and the channel with the largest energy is selected as the optimal channel. When the calculation window does not meet 5 seconds, the calculation window is temporarily changed into 1 second in order to output the screened data in real time.
In step S054, it is determined whether all three channels are smaller than the LOW ENERGY threshold LOW_HB_ENERGY_THR. The threshold is mainly a threshold for distinguishing whether or not a person is in bed. If so, no one is present and the process advances to step S055. If not, then someone goes to step S056.
In step S055, the optimal channel is the last selection result. Mainly, no one is in the state, and subsequent judgment is needed to determine the optimal channel.
In step S056, it is determined whether at least one of the three channels is greater than the HIGH ENERGY threshold HIGH_HB_ENERGY_THR. The threshold is primarily to distinguish whether there is body movement while the person is in the bed. If so, S057 is entered. If not, S058 is entered.
In step S057, the optimal channel is the last selection result. Mainly in this state there is a person in bed and there is body movement, and further subsequent judgment is required to determine the optimal channel.
In step S058, the channel with the largest energy is selected as the optimal channel. In this state, a person is in the bed and does not physically move, so that the judgment result is reliable.
Step S059 outputs the ballistocardiogram signal of the screened optimal channel. As shown in fig. 3, the three channels of raw data and the ballistocardiogram signal of the screened optimal channel are shown.
In one embodiment, step S101, after acquiring the multiple ballistocardiogram signals acquired by all the piezoelectric sensors on the mattress, the method for acquiring the ballistocardiogram signals further includes:
s141, respectively performing power frequency interference filtering treatment on the multipath ballistocardiogram signals;
s151, denoising the multipath ballistocardiogram signals subjected to the power frequency interference filtering treatment;
s161, a sliding tie filtering method of 5 points is utilized to carry out smooth denoising on the multipath ballistocardiogram signals after denoising treatment, and the obtained multipath ballistocardiogram signals are beneficial to calculating human respiratory signals.
In this embodiment, in step S141, the frequency interference filtering is performed by using a trap, the parameter of the trap is 50Hz, in step S1032, the denoising is performed by using a butterworth band-pass filter, the parameter of the butterworth band-pass filter is 0.2-0.4 Hz, after step S151 and step S141, the frequency interference and the noise of the ballistocardiogram signal are removed to a certain extent, thereby improving the accuracy of the voltage signal, and then, in step S61, the voltage signal can be converted into the ballistocardiogram signal for calculating the respiratory data.
In this embodiment, after step S161 is performed, respiratory data of the human body can be obtained by calculating the heart attack signal; and in the corresponding subsequent step S1023, the first preset threshold is a lower respiratory threshold, and in the corresponding step S1025, the second preset threshold is an upper heart rate threshold.
In one embodiment, in step S102, the preset time window includes a first time window and a second time window, where the duration of the first time window is a second duration, and the duration of the second time window is a third duration, and the third duration is less than the second duration and greater than or equal to the first duration.
Specifically, referring to fig. 4, in the above embodiment, when acquiring the ballistocardiogram signal beneficial to calculation of respiration, the method for acquiring the ballistocardiogram signal is as follows:
step S071, step 1S, iteratively calculating the respiratory energy value in the 5 second time window.
Step S072, judging whether the current data time length is more than or equal to 5 seconds. If yes, step S074 is entered, and if no, step S073 is entered.
Step S073, step 1S, iteratively calculating the breathing energy within 1S of the time window, and selecting the channel with the largest energy as the optimal channel. When the calculation window does not meet 5 seconds, the calculation window is temporarily changed into 1 second in order to output the screened data in real time.
Step S074, determining whether all three channels are less than the LOW ENERGY threshold LOW_RESP_ENERGY_THR. The threshold is mainly a threshold for distinguishing whether or not a person is in bed. If so, then no one is present and step S075 is entered. If not, then there is a person, proceeding to step S076.
In step S075, the optimal channel is the last selection result. Mainly, no one is in the state, and subsequent judgment is needed to determine the optimal channel.
Step S076, judging whether at least one channel in the three channels is larger than a HIGH ENERGY threshold value HIGH_RESP_ENERGY_THR. The threshold is primarily to distinguish whether there is body movement while the person is in the bed. If so, S077 is entered. If not, S078 is entered.
Step S077, the optimal channel is the last selection result. Mainly in this state there is a person in bed and there is body movement, and further subsequent judgment is required to determine the optimal channel.
In step S078, the channel with the largest energy is selected as the optimal channel. In this state, a person is in the bed and does not physically move, so that the judgment result is reliable.
Step S079, outputting the heart attack signal of the screened optimal channel. As shown in fig. 5, the three channels of raw data and the output ballistocardiogram signal of the screened optimal channel.
In step S101, after acquiring the multiple ballistocardiogram signals acquired by all the piezoelectric sensors on the mattress, the method for acquiring the ballistocardiogram signals further includes:
s201, judging whether the duration of the ballistocardiogram signal is greater than or equal to a second duration;
s301, if the time duration is less than the second time duration, executing step S102 by using the second time window as a preset time window;
s401, if the time duration is longer than the second time duration, executing step S102 by using the first time window as a preset time window.
In this embodiment, the first time period is 1 second, the second time period is 5 seconds, and the third time period is 1 second.
In step S101, for some reasons, such as sleeping posture of the human body, the ballistocardiogram signal acquired by a certain piezoelectric sensor may not be less than 5 seconds, that is, the second duration requirement of the second time window cannot be met, so it is required to determine whether the duration of the ballistocardiogram signal is less than 5 seconds on the basis of step S101, if so, step S301 is executed, and if so, step S401 is executed, if so, the second duration is greater than or equal to 5 seconds.
In step S301, since the time period of the ballistocardiogram signal is less than 5 seconds and exceeds the duration of the first time window, in order to output data in real time, it is necessary to slide the second time window by a step of 1 second to perform iterative calculation.
In step S1026, since the time period of the ballistocardiogram exceeds 5 seconds, the duration of the second time window is too short, and if the second time window is used, the calculation amount and the resource occupancy rate need to be increased, and the calculation speed is reduced, so that the calculation amount and the resource occupancy rate can be reduced by using the first time window, and the calculation speed is improved.
In one embodiment, the method for acquiring a ballistocardiogram signal further comprises: and displaying the voltage signal and the optimal ballistocardiogram signal.
In this embodiment, the ballistocardiogram signal and the optimal ballistocardiogram signal are displayed, so that each signal can be observed and researched more intuitively.
In one embodiment, before acquiring the multi-channel cardiac impulse signals acquired by all the piezoelectric sensors on the mattress in step S01, a hardware device with three piezoelectric sensors needs to be configured, and the positions of the three piezoelectric sensors have great influence on the accuracy of data acquisition, so that the distances between the three piezoelectric sensors need to be calculated, so that the three piezoelectric sensors are arranged at equal intervals, and the distance between the piezoelectric sensors at two ends is greater than the distance between the heart and the abdomen of a human body.
The distance from the heart to the abdomen of the human body is obtained from pre-acquired human body data, and the method for pre-acquiring the human body data comprises the following steps: marking the central positions of the heart and the abdomen in sleeping postures such as lying forward/recumbent/lying sideways/lying forward/lying sideways and the like of at least 100 adult experimenters (age 36.3+ -4.5, height 168.6 cm+ -11.2 cm); the positions of the three piezoelectric sensors are determined based on the center positions of the marked heart and abdomen so that the three piezoelectric sensors are equally spaced and fully cover the center positions of the marked heart and abdomen.
The center positions of the heart and the abdomen are marked, and the positions of the three piezoelectric sensors are determined by adopting manual marking and manual calculation methods.
In other embodiments, marking the center position of the heart and abdomen, determining the position of the three piezoelectric sensors may also employ automatic marking and automatic calculation methods. In this embodiment, the method for acquiring a ballistocardiogram signal further includes: the distance between the different piezoelectric sensors is calculated in advance so that the distance between the piezoelectric sensors at the two ends is larger than the distance between the heart and the abdomen of the human body. Wherein pre-calculating the distance between the different piezoelectric sensors includes: marking the heart center position and the abdomen center position of the human body in advance under various sleeping postures of the human body; in case all piezoelectric sensors are arranged at equal intervals so that all piezoelectric sensors can cover the heart center position and the abdomen center position, the distance between the different piezoelectric sensors is calculated.
Referring to fig. 6, the embodiment of the invention further provides a device for acquiring a ballistocardiographic signal, which comprises a mattress 1, at least three piezoelectric sensors 2 and a host unit 3.
At least three piezoelectric sensors 2 are arranged on the mattress 1 along the length direction of the mattress 1, and the collection areas formed by all the piezoelectric sensors 2 completely cover the heart center position and the abdomen center position of a human body when the human body is laid on the mattress; the host unit 3 is electrically connected to all piezoelectric sensors 2, and the host unit 3 is used for converting the voltage signals of the piezoelectric sensors 2 into ballistocardiographic signals.
Wherein the host unit 3 includes: a signal acquisition subunit, a data acquisition subunit and a signal update subunit; the signal acquisition subunit is used for acquiring multipath ballistocardiographic signals acquired by all piezoelectric sensors on the mattress; the data acquisition subunit is used for respectively sliding a preset time window on the multipath ballistocardiogram signals by taking a preset first time length as a step length so as to obtain signal data of the multipath ballistocardiogram signals in the preset time window; the signal updating subunit is used for determining a data channel according to the optimal signal data in the multipath signal data and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
According to the heart attack signal acquisition device provided by the embodiment, if a human body is sleeping, the sleeping gesture changes, at least three piezoelectric sensors can be used for acquiring heart attack signals in different areas, and as the areas acquired by the piezoelectric sensors are different, in the acquired signal data, the optimal signal data is closest to real data, and the heart attack signals are updated by using the optimal signal data, so that the heart attack signals are accurate as a whole, the errors between the heart attack signals and the real data are small, and the related data of the human body during sleeping can be calculated accurately according to the heart attack signals.
In one embodiment, the collecting device of ballistocardiogram signals further comprises: the device comprises a first trap, a first Butterworth band-pass filter and a differential unit, wherein the frequency of the first trap is 50Hz, and the first trap is used for performing work-frequency interference filtering treatment on a heart attack signal; the frequency of the first Butterworth band-pass filter is 0.67-10Hz, and the first Butterworth band-pass filter is used for denoising the heart attack signal; the differential unit is used for highlighting J points corresponding to heart beats of the heart shock signals by using a first-order differential method to obtain heart shock signals corresponding to voltage signals acquired by each piezoelectric sensor, and the heart shock signals are used for calculating heart rate data.
In one embodiment, the collecting device of ballistocardiogram signals further comprises: the frequency of the second trap is 50Hz, and the second trap is used for carrying out work-frequency interference filtering treatment on the heart attack signal; the frequency of the second Butterworth band-pass filter is 0.2-0.4 Hz, and the second Butterworth band-pass filter is used for denoising the heart attack signal; the filtering unit is used for smoothing and denoising the ballistocardiogram signals by using a sliding tie filtering method of 5 points to obtain ballistocardiogram signals corresponding to the voltage signals acquired by the piezoelectric sensor, wherein the ballistocardiogram signals are used for calculating respiratory data.
In one embodiment, the signal update subunit comprises: and the optimal signal data acquisition module and the updating module.
The optimal signal data acquisition module is used for acquiring optimal signal data in the multipath signal data, wherein the optimal signal data are used for indicating that the energy values of the ballistocardiogram signals acquired by all the piezoelectric sensors are the largest in the same time period;
the updating module is used for acquiring a data channel corresponding to the optimal signal data, and using the ballistocardiogram signal of the last step length on the data channel to update the output ballistocardiogram signal.
In one embodiment, the data acquisition subunit comprises: the heart attack signal acquisition module and the iterative calculation module; the heart attack signal acquisition module is used for taking a preset first time length as a step length and taking back heart attack signals in a preset time window; the iterative computation module is used for iteratively computing the ballistocardiogram signal in a preset time window so as to obtain signal data in the preset time window.
In one embodiment, the collecting device of ballistocardiogram signals further comprises: a first judging unit and a first output unit; the first judging unit is used for judging whether the signal data in each preset time window are smaller than a first preset threshold value or not, wherein the first preset threshold value is the lower limit value of the heart attack signal of a normal user; and the first output unit is used for updating the output ballistocardiogram signal of the data channel determined by the last sliding time window if the judgment of the first judgment unit is yes.
In one embodiment, the collecting device of ballistocardiogram signals further comprises: a second judging unit and a second output unit; the second judging unit is used for judging whether the signal data in each preset time window are all larger than a second preset threshold value, and the second preset threshold value is the upper limit value of the heart attack signal of the normal user; and the second output unit is used for updating the ballistocardiogram signal of the data channel determined by the last sliding time window if the judgment of the second judgment unit is yes.
In one embodiment, the collecting device of ballistocardiogram signals further comprises: the display module is used for displaying and processing the ballistocardiogram signal, the first optimal ballistocardiogram signal and the second optimal ballistocardiogram signal.
FIG. 7 illustrates an internal block diagram of a computer device in one embodiment. The computer device may specifically be a terminal or a server. As shown in fig. 5, the computer apparatus includes a processor, a memory, and a network interface connected by a device bus. The memory includes a nonvolatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores the operating means and may also store a computer program which, when executed by the processor, causes the processor to implement a method of acquiring ballistocardiogram signals. The internal memory may also store a computer program which, when executed by the processor, causes the processor to perform a method of acquiring ballistocardiogram signals. It will be appreciated by those skilled in the art that the structure shown in fig. 5 is merely a block diagram of some of the structures associated with the present application and is not limiting of the computer device to which the present application may be applied, and that a particular computer device may include more or fewer components than shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided comprising a memory and a processor, the memory storing a computer program that, when executed by the processor, causes the processor to perform the steps of:
the first aspect of the invention provides a method for acquiring a ballistocardiogram signal, which comprises the following steps: acquiring multipath heart attack signals acquired by all piezoelectric sensors on a mattress, wherein the mattress at least comprises three piezoelectric sensors arranged along the length direction of the mattress; respectively sliding a preset time window on the plurality of paths of ballistocardiogram signals by taking a preset first time length as a step length to obtain signal data of the plurality of paths of ballistocardiogram signals in the preset time window; and determining a data channel according to the optimal signal data in the plurality of paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel.
In one embodiment, a computer readable storage medium is provided, storing a computer program which, when executed by a processor, causes the processor to perform the steps of: acquiring voltage signals acquired by all piezoelectric sensors on a mattress, wherein the mattress at least comprises three piezoelectric sensors arranged along the length direction of the mattress; preprocessing the voltage signals acquired by each piezoelectric sensor respectively to obtain ballistocardiogram signals; dividing the ballistocardiogram signals according to preset time periods, and extracting the optimal result of the ballistocardiogram signals in each time period; and combining the optimal results of the ballistocardiographic signals of each period to generate an optimal ballistocardiographic signal.
Those skilled in the art will appreciate that all or part of the processes in the methods of the above embodiments may be implemented by a computer program for instructing relevant hardware, where the program may be stored in a non-volatile computer readable storage medium, and where the program, when executed, may include processes in the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the various embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile memory can include Read Only Memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples only represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the present application. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.

Claims (8)

1. A method of acquiring a ballistocardiogram signal, the method comprising:
acquiring multipath heart attack signals acquired by all piezoelectric sensors on a mattress, wherein the mattress at least comprises three piezoelectric sensors arranged along the length direction of the mattress;
respectively sliding a preset time window on the plurality of paths of ballistocardiogram signals by taking a preset first time length as a step length to obtain signal data of the plurality of paths of ballistocardiogram signals in the preset time window;
determining a data channel according to the optimal signal data in the multiple paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel;
the determining a data channel according to the optimal signal data in the multiple paths of signal data, and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel comprises:
acquiring optimal signal data in multiple paths of signal data, wherein the optimal signal data are used for indicating that the energy values of the ballistocardiogram signals acquired by all piezoelectric sensors are the largest in the same time period; wherein the same time period refers to a time period within the same time window;
acquiring a data channel corresponding to the optimal signal data, and using the ballistocardiogram signal of the last step length on the data channel to update an output ballistocardiogram signal;
the step of sliding a preset time window on the multiple ballistocardiograph signals with a preset first time length as a step length to obtain signal data of the multiple ballistocardiograph signals in the preset time window includes:
taking a preset first time length as a step length, and taking out the heart attack signal in the preset time window backwards;
iteratively calculating the ballistocardiogram signal in the preset time window to obtain signal data in the preset time window; and the iterative calculation is to frame and select the ballistocardiogram signal in the time window after each preset time window slides.
2. The method of claim 1, wherein after iteratively computing the ballistocardiogram signal within the preset time window to obtain signal data within the preset time window, the method further comprises:
judging whether the signal data in each preset time window are smaller than a first preset threshold value or not, wherein the first preset threshold value is the lower limit value of the heart attack signal of a normal user;
if yes, the ballistocardiogram signal of the data channel determined by the last sliding time window is used for updating the output ballistocardiogram signal.
3. The method of claim 1, wherein after iteratively computing the ballistocardiogram signal within the preset time window to obtain signal data within the preset time window, the method further comprises:
judging whether the signal data in each preset time window are larger than a second preset threshold value or not, wherein the second preset threshold value is the upper limit value of the ballistocardiogram signal of a normal user;
if yes, the ballistocardiogram signal of the data channel determined by the last sliding time window is used for updating the output ballistocardiogram signal.
4. A method according to any one of claims 1-3, wherein after the acquisition of the multiple ballistocardiographic signals acquired by all the piezoelectric sensors on the mattress, the method further comprises:
respectively carrying out power frequency interference filtering treatment on the multipath ballistocardiogram signals;
respectively denoising the multipath heart attack signals subjected to the power frequency interference filtering treatment;
by utilizing a first-order difference method, J points corresponding to central beats of the de-noised multipath heart attack signals are highlighted, and the obtained multipath heart attack signals are favorable for calculating heart rate signals of a human body.
5. A method according to any one of claims 1-3, wherein after the acquisition of the multiple ballistocardiographic signals acquired by all the piezoelectric sensors on the mattress, the method further comprises:
respectively carrying out power frequency interference filtering treatment on the multipath ballistocardiogram signals;
denoising the multipath heart attack signals subjected to the power frequency interference filtering treatment;
and a sliding tie filtering method of 5 points is utilized to carry out smooth denoising on the multipath ballistocardiogram signals after denoising treatment, and the obtained multipath ballistocardiogram signals are beneficial to calculating human respiratory signals.
6. A ballistocardiographic signal acquisition device, the device comprising:
a mattress;
at least three piezoelectric sensors arranged on the mattress in a arrayed manner along the length direction of the mattress;
the host unit is electrically connected with all the piezoelectric sensors and is used for converting voltage signals of the piezoelectric sensors into ballistocardiographic signals;
the host unit includes: a signal acquisition subunit, a data acquisition subunit and a signal update subunit;
the signal acquisition subunit is used for acquiring multipath heart attack signals acquired by all piezoelectric sensors on the mattress;
the data acquisition subunit is used for respectively sliding a preset time window on the plurality of paths of ballistocardiographic signals by taking a preset first time length as a step length so as to obtain signal data of the plurality of paths of ballistocardiographic signals in the preset time window;
the signal updating subunit is used for determining a data channel according to the optimal signal data in the plurality of paths of signal data and updating the output ballistocardiogram signal based on the ballistocardiogram signal of the data channel; acquiring optimal signal data in multiple paths of signal data, wherein the optimal signal data are used for indicating that the energy values of the ballistocardiogram signals acquired by all piezoelectric sensors are the largest in the same time period; wherein the same time period refers to a time period within the same time window; acquiring a data channel corresponding to the optimal signal data, and using the ballistocardiogram signal of the last step length on the data channel to update an output ballistocardiogram signal; taking a preset first time length as a step length, and taking out the heart attack signal in the preset time window backwards; iteratively calculating the ballistocardiogram signal in the preset time window to obtain signal data in the preset time window; and the iterative calculation is to frame and select the ballistocardiogram signal in the time window after each preset time window slides.
7. A computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of the method of any one of claims 1 to 5.
8. A computer device comprising a memory and a processor, the memory storing a computer program that, when executed by the processor, causes the processor to perform the steps of the method of any of claims 1 to 5.
CN202110873353.2A 2021-07-30 2021-07-30 Ballistocardiogram signal acquisition method and device, storage medium and computer equipment Active CN113476039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110873353.2A CN113476039B (en) 2021-07-30 2021-07-30 Ballistocardiogram signal acquisition method and device, storage medium and computer equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110873353.2A CN113476039B (en) 2021-07-30 2021-07-30 Ballistocardiogram signal acquisition method and device, storage medium and computer equipment

Publications (2)

Publication Number Publication Date
CN113476039A CN113476039A (en) 2021-10-08
CN113476039B true CN113476039B (en) 2024-03-19

Family

ID=77943633

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110873353.2A Active CN113476039B (en) 2021-07-30 2021-07-30 Ballistocardiogram signal acquisition method and device, storage medium and computer equipment

Country Status (1)

Country Link
CN (1) CN113476039B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8812091B1 (en) * 2013-03-15 2014-08-19 Apn Health, Llc Multi-channel cardiac measurements
CN104545870A (en) * 2015-01-23 2015-04-29 苏州本草芙源医疗设备有限公司 Heart rate detection mattress and method
FI20155577A (en) * 2015-08-10 2017-02-11 Murata Manufacturing Co Detection of sleep phenomenon using ballistocardiography
CN107028426A (en) * 2017-04-12 2017-08-11 佛山市量脑科技有限公司 A kind of intelligent mattress system
CN108042108A (en) * 2017-12-06 2018-05-18 中国科学院苏州生物医学工程技术研究所 A kind of sleep quality monitoring method and system based on body shake signal
CN109381174A (en) * 2017-08-09 2019-02-26 三星电子株式会社 Method and electronic device for electronic device
CN110115574A (en) * 2018-02-07 2019-08-13 普天信息技术有限公司 The method and apparatus of rhythm of the heart
US10426399B1 (en) * 2007-06-08 2019-10-01 Cleveland Medial Devices Inc. Method and device for in-home sleep and signal analysis
CN112869733A (en) * 2021-01-08 2021-06-01 广州中科新知科技有限公司 Real-time heart beat interval measuring and calculating method for ballistocardiogram

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8262582B2 (en) * 2008-08-22 2012-09-11 Valtion Teknillinen Tutkimuskeskus Extraction of heart inter beat interval from multichannel measurements
US9375152B2 (en) * 2012-03-07 2016-06-28 Cardiac Pacemakers, Inc. Heart sound detection systems and methods using updated heart sound expectation window functions
CN106999072B (en) * 2014-11-05 2021-03-12 新加坡科技研究局 Multi-channel ballistocardiograph with cepstral smoothing and quality-based dynamic channel selection
CA3090948A1 (en) * 2018-02-15 2019-08-22 Biosency Monitoring device for monitoring a physiological parameter and methods thereof

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10426399B1 (en) * 2007-06-08 2019-10-01 Cleveland Medial Devices Inc. Method and device for in-home sleep and signal analysis
US8812091B1 (en) * 2013-03-15 2014-08-19 Apn Health, Llc Multi-channel cardiac measurements
CN104545870A (en) * 2015-01-23 2015-04-29 苏州本草芙源医疗设备有限公司 Heart rate detection mattress and method
FI20155577A (en) * 2015-08-10 2017-02-11 Murata Manufacturing Co Detection of sleep phenomenon using ballistocardiography
CN107028426A (en) * 2017-04-12 2017-08-11 佛山市量脑科技有限公司 A kind of intelligent mattress system
CN109381174A (en) * 2017-08-09 2019-02-26 三星电子株式会社 Method and electronic device for electronic device
CN108042108A (en) * 2017-12-06 2018-05-18 中国科学院苏州生物医学工程技术研究所 A kind of sleep quality monitoring method and system based on body shake signal
CN110115574A (en) * 2018-02-07 2019-08-13 普天信息技术有限公司 The method and apparatus of rhythm of the heart
CN112869733A (en) * 2021-01-08 2021-06-01 广州中科新知科技有限公司 Real-time heart beat interval measuring and calculating method for ballistocardiogram

Also Published As

Publication number Publication date
CN113476039A (en) 2021-10-08

Similar Documents

Publication Publication Date Title
US20230131876A1 (en) Systems and methods of identity analysis of electrocardiograms
WO2019100562A1 (en) Automatic recognition and classification method for electrocardiogram heartbeat based on artificial intelligence
US9380954B2 (en) Method for physiological signal analysis and its system and computer program product storing physiological signal analysis program
US9198634B2 (en) Medical decision support system
US20160242665A1 (en) Devices and methods for real-time denoising of electrocardiograms
KR101910982B1 (en) Method and apparatus for eliminating motion artifact of biosignal using personalized biosignal pattern
JP2015510810A5 (en)
JP6689183B2 (en) State estimation device, method and program
CN107233093B (en) R wave detection method and device and electronic equipment
CN109411042A (en) Ecg information processing method and electro cardio signal workstation
KR101366101B1 (en) System and method for classificating normal signal of personalized ecg
CN110693488A (en) Electrocardiosignal processing system, electrocardiosignal processing method, electronic equipment and readable storage medium
CN110522443B (en) Atrioventricular conduction block detection method and device based on electrocardiosignals and electronic equipment
JP6243254B2 (en) Heart rate detection method and heart rate detection device
US20220125328A1 (en) Time features calculation apparatus, calculation method and its program
CN115315214A (en) Fusion signal processing for maternal uterine activity detection
US9549681B2 (en) Matrix-based patient signal analysis
CN113476039B (en) Ballistocardiogram signal acquisition method and device, storage medium and computer equipment
CN110179451B (en) Electrocardiosignal quality detection method and device, computer equipment and storage medium
EP3265937B1 (en) Systems for displaying medical data
Yadav et al. ECG signal characterization using Lagrange-Chebyshev polynomials
CN109044338B (en) Atrial fibrillation detection apparatus and storage medium
JP2018094156A (en) Instantaneous heart rate reliability evaluation device, method and program
JP6857582B2 (en) Instantaneous heart rate time series data complement device, complement method and its program
JP6468986B2 (en) Noise determination apparatus, method, and program

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