WO2024166481A1 - 生体情報計測装置 - Google Patents

生体情報計測装置 Download PDF

Info

Publication number
WO2024166481A1
WO2024166481A1 PCT/JP2023/041263 JP2023041263W WO2024166481A1 WO 2024166481 A1 WO2024166481 A1 WO 2024166481A1 JP 2023041263 W JP2023041263 W JP 2023041263W WO 2024166481 A1 WO2024166481 A1 WO 2024166481A1
Authority
WO
WIPO (PCT)
Prior art keywords
frequency
signal
bioinformation
measuring device
calculation unit
Prior art date
Application number
PCT/JP2023/041263
Other languages
English (en)
French (fr)
Japanese (ja)
Inventor
直明 松岡
卓 草薙
研人 藤木
愛 河原
Original Assignee
株式会社村田製作所
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 株式会社村田製作所 filed Critical 株式会社村田製作所
Priority to JP2024576115A priority Critical patent/JPWO2024166481A1/ja
Publication of WO2024166481A1 publication Critical patent/WO2024166481A1/ja

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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • 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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024—Measuring pulse rate or heart rate
    • A61B5/0245—Measuring 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/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb occurring during breathing

Definitions

  • the present invention relates to a bioinformation measuring device.
  • a bioinformation measuring device that measures bioinformation by analyzing biosignals such as pulse waves is known (Patent Document 1).
  • This bioinformation measuring device includes a phase-locked loop circuit to which the biosignal is input.
  • the phase-locked loop circuit includes a phase frequency comparator, a loop filter, and a voltage-controlled oscillator.
  • a variable low-pass filter blocks signals of a predetermined frequency band contained in the deviation signal that has passed through the loop filter. The bioinformation is obtained from the signal that has passed through the variable low-pass filter.
  • Biological signals include signals related to the biological information of the measurement target as well as signals related to biological information that is not the target of measurement.
  • a ballistocardiogram (BCG) obtained to measure heart rate includes low-frequency signals caused by breathing, etc. If signals generated by other biological phenomena are superimposed on the signals related to the biological information of the measurement target, the measurement accuracy of the biological information of the measurement target will decrease.
  • the obtained signals may also include noise caused by the environment, such as vibration.
  • the object of the present invention is to provide a bioinformation measuring device that can suppress a decrease in the measurement accuracy of bioinformation even if signals generated due to other biological phenomena or the environment are superimposed on a signal related to the bioinformation of the measurement subject.
  • a bandpass filter that receives a biological signal having a harmonic structure, passes a component of a frequency band including a fundamental frequency and one of a plurality of harmonics included in the biological signal, and attenuates components of other frequencies to output a first signal; a frequency calculation unit that receives the first signal and outputs a second signal including information about a frequency of the input signal; and a biological information acquiring unit that acquires biological information from the second signal.
  • FIG. 1A is a block diagram of a bioinformation measuring device according to a first embodiment and a diagram showing an example of a signal waveform
  • FIG. 1B is a graph showing a heart rate waveform as an example of a biosignal SigB.
  • FIG. 2 is a graph showing an example of a spectrum obtained by performing a spectral analysis on the biological signal SigB.
  • FIG. 3 is a block diagram of the bandpass filter and the frequency calculation unit.
  • FIG. 4 is a diagram showing an example of the spectrum and signal waveform of the biological signal SigB to be measured and the signal SigR caused by respiration or the like superimposed on the biological signal SigB.
  • FIG. 1A is a block diagram of a bioinformation measuring device according to a first embodiment and a diagram showing an example of a signal waveform
  • FIG. 1B is a graph showing a heart rate waveform as an example of a biosignal SigB.
  • FIG. 2 is a graph showing an example of a spectrum
  • FIG. 5 is a graph showing an example of waveforms of the biological signal SigB and the first signal Sig1, and a time change in the frequency of the fundamental wave of the biological signal SigB and the first signal Sig1.
  • FIG. 6 is a block diagram of a biological information measuring device according to the second embodiment.
  • FIG. 7 is a flowchart showing a procedure for the signal analysis unit to determine the target order n.
  • FIG. 8 is a flowchart showing a control procedure executed by the input control unit.
  • FIG. 9 is a graph for explaining the control of the input control unit.
  • FIG. 10A is a block diagram of a biological information measuring device according to the third embodiment
  • FIG. 10B is a block diagram of a frequency calculation unit.
  • FIG. 11 is a block diagram of a biological information measuring device according to the fourth embodiment.
  • FIG. 12 is a graph showing the relationship between the pass frequency bands of a plurality of filter units BPF i .
  • FIG. 13 is a block diagram of the output control unit.
  • FIG. 14 is a graph showing the relationship between the previous frequency ft p , the pass frequency band of the filter unit BPF i , and the selected phase locked loop PLL i .
  • FIG. 15 is a table showing an example of cutoff frequencies and threshold values for switching the phase locked loop.
  • FIG. 16 is a block diagram of a biological information measuring device according to a modification of the fourth embodiment.
  • FIG. 1A is a block diagram of the bioinformation measuring device according to the first embodiment, and a diagram showing an example of a signal waveform.
  • Fig. 1B is a graph showing a heart rate waveform as an example of a biosignal SigB.
  • the horizontal axis of Fig. 1B represents time, and the vertical axis represents sensor output.
  • the bioinformation measuring device includes a sensor 70, a bandpass filter 10, a frequency calculation unit 20, a bioinformation acquisition unit 30, and a display device 60.
  • the sensor 70 is, for example, an acceleration sensor that acquires a ballistocardiogram (BCG).
  • BCG ballistocardiogram
  • the sensor 70 is used by being placed around the human body, for example on a seat or bed, or by being in direct contact with the human body. Heartbeat vibrations are detected by the sensor 70.
  • the functions of the bandpass filter 10, frequency calculation unit 20, and bioinformation acquisition unit 30 are realized in software by a microcontrol unit (MCU).
  • the sensor 70 detects the biosignal SigB, and the biosignal SigB shown in FIG. 1B is input to the bandpass filter 10.
  • the signal output from the sensor 70 may be an analog signal or a digital signal. If the signal output from the sensor 70 is an analog signal, it is converted to a digital signal by an AD converter in the MCU.
  • a general formula expressing the biosignal SigB will be described.
  • a heartbeat signal is captured by a ballistocardiogram, an electrocardiogram, a pulse wave signal, or the like
  • these signals often mimic a structure composed of multiple sine waves.
  • a signal representing breathing can also be confirmed to have a similar structure if it is periodic.
  • the waveform y'(t) of these periodic biosignals SigB, such as heartbeat and breathing can often be described by the following equation:
  • t represents time
  • f r represents frequency
  • ⁇ r represents phase
  • a r represents amplitude.
  • the waveform y(t, f 0 ) of the biosignal SigB can be expressed by the following equation.
  • the first term on the right side of equation (2) represents the fundamental wave, and the second term represents the harmonic wave.
  • the fundamental frequency f0 In a biosignal such as a heartbeat signal, in addition to the change in intensity depending on time t, the fundamental frequency f0 also changes.
  • the change in fundamental frequency f0 is a factor that changes the heartbeat interval.
  • the argument of the function y includes the fundamental frequency f0 .
  • the amplitude components ( a0 , ar ), frequency components ( f0 , kr ), phase components ( ⁇ 0 , ⁇ r ), etc. become values specific to the vibration transmission path or the living body.
  • the biosignal SigB can be expressed as a sum of a fundamental wave and harmonics up to order N
  • the waveform y(t, f0 ) of the biosignal SigB can be expressed by the following formula.
  • the frequency components of a heartbeat signal (e.g., a BCG waveform) acquired by a specific sensor are f0 , 2f0 , 3f0 , 4f0 , and 5f0 , assuming that the maximum order N of the constituent signals is 5.
  • Each frequency component has its own amplitude A1 , A2 , A3 , A4 , and A5 .
  • the fundamental frequency f0 is called the heartbeat frequency
  • its reciprocal 1/ f0 is called the heartbeat (pulse wave) interval.
  • the biosignal SigB ( Figure 1B) has a harmonic structure expressed by equation (3). That is, the biosignal SigB is composed of a fundamental wave and its harmonics.
  • the period of the fundamental wave of the heartbeat waveform is marked as T.
  • the period T of the fundamental wave corresponds to the heartbeat interval, and its reciprocal corresponds to the heartbeat frequency.
  • the bandpass filter 10 passes signals in one of the fundamental frequency band and the frequency bands of the multiple harmonics of the inputted biosignal SigB, and attenuates signals in the other frequency bands.
  • the pass frequency band of the bandpass filter 10 corresponds to the frequency band of the nth harmonic of the biosignal SigB.
  • the bandpass filter 10 passes signals in the fundamental frequency band of the biosignal SigB.
  • the fundamental frequency of the biosignal SigB is f0
  • the waveform of the first signal Sig1 that has passed through the bandpass filter 10 becomes a shape close to a sine wave with a frequency of nf0 .
  • the frequency calculation unit 20 analyzes the first signal Sig1 that has passed through the band-pass filter 10, and outputs a second signal Sig2 that includes information about the frequency nf 0 of the first signal Sig1.
  • the second signal Sig2 has a value of frequency nf 0.
  • the second signal Sig2 has a voltage value corresponding to the frequency nf 0 .
  • the bioinformation acquiring unit 30 acquires the bioinformation infB from the second signal Sig2. For example, the bioinformation acquiring unit 30 obtains a fundamental frequency f0 from the second signal Sig2.
  • the biosignal infB is a heartbeat signal
  • the bioinformation infB is a heartbeat frequency
  • its value is given by the fundamental frequency f0
  • the heartbeat interval is given by its reciprocal 1/ f0 .
  • the display device 60 displays information related to the bioinformation infB acquired by the bioinformation acquisition unit 30. For example, it displays the heartbeat frequency or heartbeat interval as numbers or a graph.
  • FIG. 2 is a graph showing an example of a spectrum obtained by spectral analysis of the biosignal SigB.
  • the horizontal axis represents frequency, and the vertical axis represents signal strength for each frequency.
  • the biological signal SigB includes a fundamental wave with a fundamental frequency f0 and second to fifth harmonics with frequencies 2f0 , 3f0 , 4f0 , and 5f0 .
  • the bandpass filter 10 passes signals in a frequency band including the second harmonic frequency 2f0 , for example.
  • the lower cutoff frequency on the low frequency side of the bandpass filter 10 is higher than the fundamental frequency f0
  • the upper cutoff frequency on the high frequency side is lower than the third harmonic frequency 3f0 .
  • the bandpass filter 10 passes signals in one of the frequency bands of the fundamental wave and each of the multiple harmonics of the biological signal SigB, and attenuates signals in the other frequency bands.
  • the frequency calculation unit 20 uses the first signal Sig1 (second harmonic) that has passed through the bandpass filter 10 to calculate its frequency 2f0 , and outputs it as the second signal Sig2.
  • the order of the harmonic that is the subject of calculation by the frequency calculation unit 20 is referred to as the target order n.
  • the bioinformation acquisition unit 30 calculates the fundamental frequency f0 by dividing the value indicated by the second signal Sig2 by 2.
  • the fundamental frequency f0 is calculated without directly using the signal in the fundamental frequency band of the biological signal SigB. Note that when the bandpass filter 10 passes the signal in the fundamental frequency band and attenuates the signal in the harmonic frequency band, the second signal Sig2 output from the frequency calculation unit 20 represents the fundamental frequency f0 .
  • FIG. 3 is a block diagram of the bandpass filter 10 and the frequency calculation unit 20.
  • the functions of each block of the frequency calculation unit 20 are realized, for example, by software. Note that it is also possible to realize these functions by hardware circuits.
  • the bandpass filter 10 passes signals in a predetermined specific frequency band among the frequency bands of the fundamental wave and multiple harmonics of the biosignal SigB.
  • a fourth-order infinite impulse response (IIR) digital filter can be used as the bandpass filter 10 that passes signals in a frequency band including the second harmonic.
  • the lower cutoff frequency is set to, for example, 2 Hz, and the upper cutoff frequency is set to 4 Hz. If large output delay is not a problem or if a steep cutoff characteristic is not required, a finite impulse response (FIR) digital filter may be used as the bandpass filter 10.
  • the frequency calculation unit 20 includes a phase locked loop 21, a frequency conversion unit 26, and a low pass filter 27.
  • the phase locked loop 21 includes a phase comparison unit 22, a loop filter 23, and a numerically controlled oscillator 24.
  • the phase locked loop 21 is designed to be able to track signals in the frequency band that is passed by the band pass filter 10, among the frequency bands of the fundamental wave and multiple harmonics of the biosignal SigB.
  • the numerically controlled oscillator 24 changes the frequency and phase of the tracking signal Sigt that it outputs according to the output of the loop filter 23.
  • the initial frequency at the start of operation of the numerically controlled oscillator 24 and the range of frequencies that the tracking signal Sigt tracks (hereinafter sometimes referred to as tracking frequencies) can be initialized by an external control signal.
  • the operation of the phase locked loop circuit 21 can be stopped (tracking stopped) by an external control signal.
  • the phase locked loop circuit 21 may be a free-running phase locked loop circuit that can track a specific frequency band without parameter settings or external control signal input.
  • a voltage controlled oscillator is used instead of the numerically controlled oscillator 24.
  • the phase comparator 22 compares the input first signal Sig1 with the tracking signal Sigt output from the numerically controlled oscillator 24, and calculates the phase difference between them.
  • the loop filter 23 outputs an appropriate control signal for controlling the numerically controlled oscillator 24 based on the phase difference calculated by the phase comparator 22.
  • the frequency conversion unit 26 converts the control signal output from the loop filter 23 into frequency information. More specifically, it converts the control value of the control signal input to the numerically controlled oscillator 24 into the current tracking frequency of the phase-locked loop circuit 21. Depending on the configurations of the loop filter 23 and the numerically controlled oscillator 24, the output of the loop filter 23 may contain frequency information. In such cases, the frequency conversion unit 26 is not necessary.
  • the low-pass filter 27 smooths the time change in the control value of the control signal output from the loop filter 23. For example, depending on the design of the loop filter 23 and the numerically controlled oscillator 24, ripple noise of a magnitude that cannot be ignored may be superimposed on the output of the loop filter 23.
  • the low-pass filter 27 is provided for the purpose of removing this ripple noise. Note that if the design of the loop filter 23 and the numerically controlled oscillator 24 can suppress the ripple noise to a negligible level, or if the ripple noise is not a problem for downstream display control or applications, then the low-pass filter 27 does not need to be provided.
  • the initial frequency of the phase locked loop 21 is 2.5 Hz
  • the tracking frequency range of the numerically controlled oscillator 24 is 2 Hz to 4 Hz.
  • a fourth-order IIR digital filter is used as the low-pass filter 27, and the cutoff frequency of the low-pass filter 27 is 0.6 Hz. Note that if a large output delay is not a problem or if a steep cutoff characteristic is not required, an FIR digital filter may be used as the low-pass filter 27.
  • FIG. 4 is a diagram showing an example of the spectrum and signal waveform of the biological signal SigB to be measured and the signal SigR caused by respiration or the like superimposed on the biological signal SigB.
  • the spectrum Sph of the biological signal SigB shows peaks of the fundamental wave at frequency f 0 , the second harmonic at frequency 2f 0 , the third harmonic at frequency 3f 0 , the fourth harmonic at frequency 4f 0 , and the fifth harmonic at frequency 5f 0.
  • the spectrum Spr of the signal SigR caused by breathing, etc. appears.
  • the frequency band of the spectrum Spr overlaps with the frequency band of the fundamental wave of the biological signal SigB.
  • the heart rate is generally said to be between 60 and 85 bpm, which, when converted to frequency, gives a heart rate of between 1 and 1.4 Hz.
  • the respiratory rate is generally between 12 and 20 bpm, which, when converted to frequency, gives a respiratory rate of between 0.2 and 0.3 Hz. In terms of frequency, the two do not overlap.
  • movement caused by breathing is often several times greater than movement caused by the heart rate. When body surface movement is captured with an acceleration sensor, the signal due to movement caused by breathing appears relatively large.
  • the signal SigR caused by breathing has a gentle triangular wave shape, and as shown in FIG. 4, its harmonics can extend into the frequency band of the fundamental wave of the biosignal SigB while maintaining a large amplitude level. If a signal combining the signal SigR caused by breathing, etc. and the biosignal SigB caused by the heartbeat is directly input to the phase-locked loop 21 (FIG. 3), it will be influenced by the signal SigR caused by breathing, etc., and will make it difficult to accurately track the fundamental wave of the biosignal SigB caused by the heartbeat.
  • the second harmonic rather than the fundamental wave of the biological signal SigB, is input to the frequency calculation unit 20 (FIG. 1A), and the fundamental wave of the biological signal SigB and the signal SigR caused by breathing, etc. are attenuated by the bandpass filter 10 (FIG. 1A) and are not input to the frequency calculation unit 20 (FIG. 1A). Therefore, the frequency 2f 0 of the second harmonic can be calculated with high accuracy without being affected by the signal SigR caused by breathing, etc.
  • the biological information acquisition unit 30 calculates the frequency f 0 of the fundamental wave from the frequency 2f 0 calculated by the frequency calculation unit 20, so that the frequency f 0 of the fundamental wave can be calculated without being affected by the signal SigR caused by breathing, etc.
  • FIG. 5 is a graph showing an example of the waveforms of the biosignal SigB and the first signal Sig1, and the time change in frequency of the fundamental wave of the biosignal SigB and the first signal Sig1.
  • the horizontal axis represents time
  • the first graph represents the waveform of the biosignal SigB
  • the second graph represents the waveform of the first signal Sig1 that has passed through the bandpass filter 10. Note that the waveform shown in the upper part of FIG. 5 has signals other than BCG and noise excluded.
  • the target order n is set to 2.
  • the frequency of the first signal Sig1, which is the second harmonic is in the range of 2 Hz to 4 Hz.
  • the vertical axis of the third graph represents frequency.
  • the solid line in the third graph represents the change over time in frequency of the first signal Sig1 (i.e., the change over time in frequency of the second harmonic of the biosignal SigB), and the dashed line represents the change over time in the fundamental wave of the biosignal SigB (i.e., the change over time in heart rate frequency).
  • the frequency range of the spectrum Spr ( Figure 4) of the signal SigR resulting from breathing, etc. is hatched.
  • the heartbeat frequency may fall within the frequency range of the spectrum Spr of the signal SigR caused by breathing, etc.
  • the heartbeat frequency cannot be accurately determined based on the fundamental wave of the biosignal SigB. Since the frequency of the first signal Sig1, which is the second harmonic of the biosignal SigB, does not overlap with the frequency range of the spectrum of the signal caused by breathing, etc., the frequency of the first signal Sig1 can be accurately determined in any time period. As a result, it becomes possible to accurately determine the heartbeat frequency.
  • the bandpass filter 10 passes the second harmonic of the biosignal SigB, but may pass harmonics of other orders.
  • the frequency band of the third or higher harmonics does not overlap with the frequency band of the spectrum of the signal SigR caused by breathing or the like. Therefore, the bandpass filter 10 may pass any of the third or higher harmonics.
  • the bioinformation acquisition unit 30 may divide the value of the second signal Sig2 by the order of the harmonic passing through the bandpass filter 10.
  • the bandpass filter 10 may pass the fundamental wave and attenuate the second and higher harmonics.
  • an acceleration sensor is used as the sensor 70 (FIG. 1A), but other sensors capable of acquiring a BCG, such as a piezoelectric sensor, may also be used.
  • a signal representing a BCG waveform is adopted as the biosignal SigB input to the bandpass filter 10, but other signals may also be used.
  • a signal representing an electrocardiogram (ECG) waveform or a signal representing a pulse wave waveform may be used as the biosignal SigB.
  • ECG electrocardiogram
  • the signal representing an ECG waveform can be acquired, for example, using an electrocardiograph.
  • the signal representing a pulse wave waveform can be acquired, for example, using a photoplethysmograph sensor.
  • the biosignal SigB is a signal corresponding to the heartbeat of the living body, but the biosignal SigB may be a signal having another harmonic structure.
  • a signal that changes in response to respiration may be used as the biosignal SigB.
  • the signal to be calculated by the frequency calculation unit 20 is selected from a fundamental wave or a harmonic having a frequency that is an integer multiple of the frequency of the fundamental wave, but the frequency band of the signal to be calculated by the frequency calculation unit 20 does not have to be an integer multiple of the fundamental frequency.
  • FIG. 6 is a block diagram of a bioinformation measuring device according to the second embodiment.
  • the bioinformation measuring device according to the second embodiment includes a signal analysis unit 40 and an input control unit 50 in addition to the sensor 70, bandpass filter 10, frequency calculation unit 20, bioinformation acquisition unit 30, and display unit 60 of the bioinformation measuring device according to the first embodiment.
  • a biosignal SigB is input to the signal analysis unit 40.
  • the harmonic order n (target order n) that the bandpass filter 10 passes is determined in advance, but in the second embodiment, the signal analysis unit 40 analyzes the biosignal SigB to determine the target order n.
  • FIG. 7 is a flowchart showing the procedure by which the signal analysis unit 40 determines the target order n.
  • the signal analysis unit 40 analyzes the input signal (step SA1). Based on the analysis result, it is determined whether or not the input signal contains a biosignal SigB such as a heartbeat signal (step SA2). For example, in step SA1, the input signal is Fourier transformed, and if the result of the Fourier transform has a harmonic structure expressed by equation (3), it is determined in step SA2 that a biosignal SigB is present. Alternatively, in step SA1, the root mean square (RMS) of the intensity of the input signal is calculated, and if the calculation result is equal to or greater than a threshold value, it is determined in step SA2 that a biosignal SigB is present.
  • RMS root mean square
  • steps SA1 and SA2 are repeated until it is determined that the biosignal SigB is present. If it is determined that the input signal contains the biosignal SigB, the target order n is determined, and various parameters of the bandpass filter 10 and the frequency calculation unit 20 are determined and set (step SA3).
  • the target order n may be, for example, the order of a peak that appears in the frequency band with the smallest noise floor on the frequency axis in the spectrum obtained by performing a Fourier transform or the like on the biosignal SigB.
  • the cutoff frequency of the bandpass filter 10 the initial frequency of the numerically controlled oscillator 24 ( Figure 3), and the parameters of the loop filter 23 are determined and set.
  • the value of the second signal Sig2 output from the frequency calculation section 20 becomes the frequency nf0 of the nth harmonic.
  • the bioinformation acquisition unit 30 includes a divider 32 and an inverse calculator 33.
  • the signal analysis unit 40 notifies the divider 32 of the target order n.
  • the divider 32 divides the frequency nf 0 indicated by the second signal Sig2 input from the frequency calculation unit 20 by the target order n to generate bioinformation infB indicating the value of the fundamental frequency f 0 of the biosignal SigB.
  • the bioinformation infB represents the heartbeat frequency.
  • the inverse calculator 33 calculates the inverse of the bioinformation infB output from the divider 32 to generate information T representing the heartbeat interval.
  • the obtained heartbeat frequency and heartbeat interval information are displayed on the display device 60.
  • FIG. 8 is a flowchart showing the control procedure executed by the input control unit 50.
  • the input control unit 50 analyzes the first signal Sig1 output from the bandpass filter 10 (step SB1). Based on the analysis result, it is determined whether the bandpass filter 10 outputs a first signal Sig1 (hereinafter referred to as a significant first signal Sig1) whose magnitude can be calculated by the frequency calculation unit 20 (step SB2). If the bandpass filter 10 outputs a significant first signal Sig1, it turns on the signal input to the frequency calculation unit 20 (step SB3). This causes the frequency calculation unit 20 to calculate the frequency (step SB4). If the bandpass filter 10 does not output a significant first signal Sig1, it turns off the signal input to the frequency calculation unit 20 (step SB5). In other words, no signal is input to the frequency calculation unit 20. Furthermore, it initializes the frequency calculation unit 20 (step SB6).
  • FIG. 9 is a graph for explaining the control of the input control unit 50.
  • the upper graph in FIG. 9 represents the time change in the root mean square (RMS) of the intensity of the first signal Sig1 output from the bandpass filter 10, and the lower graph is a timing chart of the on/off of the signal input to the frequency calculation unit 20 (FIG. 6).
  • the horizontal axis represents time
  • the vertical axis of the upper graph represents RMS
  • the vertical axis of the lower graph represents the on/off of the signal input to the frequency calculation unit 20.
  • the input control unit 50 calculates the RMS of the first signal Sig1.
  • the input control unit 50 switches the signal input to the frequency calculation unit 20 from off to on.
  • the signal input to the frequency calculation unit 20 is currently on, if the RMS value becomes equal to or less than the off threshold THoff (time t2 ), the input control unit 50 switches the signal input to the frequency calculation unit 20 from on to off.
  • the on threshold THon is greater than the off threshold THoff.
  • the input of the signal to the frequency calculation unit 20 can be turned off to prevent the phase synchronization circuit 21 ( Figure 3) of the frequency calculation unit 20 from tracking a signal of an incorrect frequency.
  • the input control unit 50 may have functions to adjust the gain of the input signal to the frequency calculation unit 20 and adjust the sampling rate.
  • step SB1 the output signal of the bandpass filter 10 is Fourier transformed, and in the determination of step SB2, if a peak exists in the frequency band of the target order n, it may be determined that a significant first signal Sig1 is being output.
  • the target order n is determined by analyzing the signal input from the sensor 70 (FIG. 6). For this reason, among the fundamental wave and multiple harmonics of the biological signal SigB, the fundamental wave or the harmonic of the order that is least susceptible to noise is input to the frequency calculation unit 20. By determining a preferable target order n according to the noise situation, the calculation accuracy of the frequency calculation unit 20 can be improved.
  • FIG. 10A is a block diagram of a bioinformation measuring device according to the third embodiment
  • FIG. 10B is a block diagram of a frequency calculation unit 20.
  • the control value that the loop filter 23 (FIG. 3) uses to control the numerically controlled oscillator 24 is output from the frequency calculation unit 20 as a second signal Sig2 via the frequency conversion unit 26 and the low-pass filter 27.
  • the frequency calculation unit 20 includes a frequency counter 29.
  • the tracking signal Sigt generated by the numerically controlled oscillator 24 is input to the frequency counter 29.
  • the frequency counter 29 counts the frequency of the tracking signal Sigt and outputs the counting result as the second signal Sig2.
  • the frequency counter 29 has a function of sequentially resetting the counter value internally.
  • a function may be provided to initialize the counter value by a command from the input control unit 50 depending on the state of the first signal Sig1 that has passed through the bandpass filter 10.
  • the tracking signal Sigt generated by the numerically controlled oscillator 24 is divided by the frequency divider 25 and input to the phase comparison unit 22.
  • the frequency division ratio m is set so that the frequency of the signal divided by the frequency divider 25 and the frequency of the first signal Sig1 are included in the same frequency band.
  • a fractional PLL may be used as the frequency calculation unit 20.
  • the control value output by the loop filter 23 to the numerically controlled oscillator 24 becomes nf 0.
  • the numerically controlled oscillator 24 generates a signal having a frequency of mnf 0 as the tracking signal Sigt. Therefore, the value of the second signal Sig2 output from the frequency counter 29 becomes mnf 0 .
  • the input control unit 50 has the same functions as the input control unit 50 (FIG. 6) of the bioinformation measuring device according to the second embodiment.
  • the input control unit 50 may have a function of initializing the phase synchronization circuit 21 and the frequency counter 29 of the frequency calculation unit 20 in step SB6 of FIG. 8.
  • the signal analysis unit 40 has the same functions as the signal analysis unit 40 (Fig. 6) of the bioinformation measuring device according to the second embodiment. Furthermore, the signal analysis unit 40 determines the division ratio m to be set in the divider 25 based on the oscillation frequency of the numerically controlled oscillator 24 (Fig. 10B) and the target order n. The division ratio m can be determined by dividing the oscillation frequency of the numerically controlled oscillator 24 (Fig. 10B) when the heartbeat frequency is 1 Hz by the target order n.
  • the division ratio m is 512. Also, when the target order n is 4, the division ratio m is 256.
  • the oscillation frequency of the numerically controlled oscillator 24 is designed to be 1000 Hz when the heartbeat frequency is 1 Hz. If the heartbeat frequency range is assumed to be 0.5 Hz or more and 2.5 Hz or less, a numerically controlled oscillator 24 capable of oscillating in a frequency range of 500 Hz or more and 2500 Hz or less is used. When the target order n is 2, the division ratio m should be set to 500.
  • the signal analysis unit 40 notifies the phase synchronization circuit 21 of the frequency calculation unit 20 of the division ratio m, and notifies the divider 32 of the bioinformation acquisition unit 30 of the division ratio m and the target order n.
  • the divider 25 of the phase synchronization circuit 21 divides the tracking signal Sigt by the division ratio m, and provides the divided signal to the phase comparison unit 22.
  • the divider 32 obtains bioinformation infB by dividing the value mnf0 of the second signal Sig2 counted by the frequency counter 29 by the product of the target order n and the division ratio m.
  • the value of the bioinformation infB is equal to the fundamental frequency f0 of the bioinformation measuring device SigB.
  • the reciprocal calculator 33 has the same function as the reciprocal calculator 33 of the bioinformation measuring device according to the second embodiment.
  • the product of the target order n and the division ratio m may be stored in the divider 32 as the heartbeat calculation counter value.
  • the divider 32 can obtain the bioinformation infB by dividing the value of the second signal Sig2 by the heartbeat calculation counter value.
  • the bioinformation infB (heartbeat frequency) is 1.025.
  • the calculation accuracy of the frequency calculation unit 20 can be improved by determining a preferable target order n depending on the noise situation.
  • the loop filter 23 (FIG. 3) outputs the second signal Sig2 by counting the frequency of the tracking signal Sigt generated by the phase locked loop circuit 21 with the frequency counter 29 without referring to the control value that controls the numerically controlled oscillator 24. This makes it easier to realize the function of the frequency calculation unit 20 with hardware.
  • the counter value by the frequency counter 29 may be interpreted as a fixed-point number, and the divider 32 may be omitted.
  • the product of the target order n and the division ratio m is 1000, it may be interpreted that there is a decimal point between the hundreds and thousands digits of the counter value by the frequency counter 29.
  • the counter value is 1000, the bioinformation infB (heartbeat frequency) is 1.000, and when the counter value is 825, the bioinformation infB (heartbeat frequency) is 0.825.
  • the counter value may be displayed as it is on the display device 60.
  • Fig. 11 is a block diagram of a bioinformation measuring device according to the fourth embodiment.
  • the bandpass filter 10 includes a plurality of filter units BPF i having different pass frequency bands.
  • i represents a serial number assigned to the plurality of filter units.
  • a common biosignal SigB is input to the plurality of filter units BPF i .
  • Fig. 12 is a graph showing the relationship of the pass frequency band of a plurality of filter units BPF i .
  • the horizontal axis represents frequency, and the vertical axis represents pass rate.
  • Fig. 12 shows an example in which the number of filter units BPF i is five.
  • the plurality of filter units BPF i are labeled BPF 1 , BPF 2 , BPF 3 , BPF 4 , and BPF 5 in order from the lowest center frequency of the pass frequency band.
  • the pass frequency bands of two filter units BPF i whose center frequencies are adjacent to each other partially overlap. For example, a part of the high frequency side of the pass frequency band of filter unit BPF i and a part of the low frequency side of the pass frequency band of filter unit BPF i+1 overlap with each other. Also, the pass frequency bands of filter units BPF i whose center frequencies are not adjacent to each other do not overlap with each other. For example, the pass frequency band of filter unit BPF i and the pass frequency band of filter unit BPF i+2 do not overlap with each other.
  • the lower cutoff frequency of the filter unit BPF i is denoted as f iL
  • the upper cutoff frequency is denoted as f iH .
  • the cutoff frequencies of the multiple filter units BPF i are set so that the following inequality is satisfied. f i+1,L ⁇ f iH f iH ⁇ f i+2,L
  • the frequency calculation section 20 (FIG. 11) includes a plurality of phase synchronization circuits 21 and an output control section 28.
  • the number of phase synchronization circuits 21 is the same as the number of filter units BPF i , and there is a one-to-one correspondence between the phase synchronization circuits 21 and the filter units BPF i .
  • the phase synchronization circuits 21 corresponding to the filter units BPF 1 , BPF 2 , BPF 3 , BPF 4 , and BPF 5 are denoted as PLL 1 , PLL 2 , PLL 3 , PLL 4 , and PLL 5 , respectively.
  • the first signal Sig1 that has passed through the filter unit BPF i is input to the corresponding phase locked loop PLL i .
  • the phase locked loop PLL i outputs the frequency (tracking frequency) value nf 0 of the tracking signal Sigt.
  • the output control unit 28 selects one of the tracking frequency values output from each of the multiple phase locked loops PLL i and outputs it as the second signal Sig2.
  • the function of the biometric information acquisition unit 30 is the same as the function of the biometric information acquisition unit 30 of the biometric information measuring device according to the second embodiment.
  • the input control unit 50 inputs the first signal Sig1, which is the output of the filter unit BPF i , to the corresponding phase locked loop PLL i . Furthermore, the input control unit 50 performs the same control as the input control unit 50 of the bioinformation measuring device according to the second embodiment (FIG. 6) for the pair of the mutually corresponding filter unit BPF i and the phase locked loop PLL i .
  • the input control unit 50 also controls the output control unit 28.
  • the control of the output control unit 28 will be explained later.
  • the signal analysis unit 40 has the same functions as the signal analysis unit 40 of the bioinformation measuring device according to the second embodiment (FIG. 6). For example, the signal analysis unit 40 determines the target order n by analyzing the input biosignal SigB. Once the target order n is determined, the signal analysis unit 40 determines the cutoff frequency of each of the multiple filter units BPF i , and sets the cutoff frequency to each of the filter units BPF i . Furthermore, the signal analysis unit 40 individually controls each of the multiple phase-locked loops PLL i .
  • the signal analysis unit 40 also controls the output control unit 28.
  • the control of the output control unit 28 will be explained later.
  • Fig. 13 is a block diagram of the output control unit 28.
  • the output control unit 28 includes a selector 28A and a delay unit 28B.
  • the tracking frequency ft i output from each of the phase locked loops PLL i is input to the selector 28A.
  • the tracking frequency ft output from the selector 28A is delayed by the delay unit 28B and input to the selector 28A.
  • the tracking frequency delayed by the delay unit 28B and input to the selector 28A is referred to as the previous frequency ft p .
  • the selector 28A compares the previous frequency ftp with a threshold value set in the selector 28A, and determines the tracking frequency ft to be output this time from among a plurality of tracking frequencies fti based on the comparison result.
  • An initial value finit of the previous frequency ftp is set in the delay unit 28B.
  • the delay unit 28B outputs the initial value finit as the previous frequency ftp .
  • the initial value finit is set by a command from the signal analysis unit 40 or the input control unit 50.
  • Fig. 14 is a graph showing the relationship between the previous frequency ftp , the pass frequency band of the filter unit BPF i , and the selected phase locked loop PLL i .
  • the horizontal axis represents the previous tracking frequency ftp
  • the vertical axis represents the selected phase locked loop PLL i .
  • the lower cutoff frequency of the pass frequency band of the filter unit BPF i is marked as f iL
  • the upper cutoff frequency is marked as f iH .
  • An upper threshold value f i,i+1 and a lower threshold value f i+1,i are set in the range where the pass frequency bands of filter units BPF i and BPF i+1 overlap.
  • the upper threshold value f i,i+1 is higher than the lower threshold value f i+1,i .
  • a case where the output of phase locked loop circuit PLL i is currently selected will be described.
  • the output of phase locked loop circuit PLL i is also selected this time.
  • the output of the phase locked loop PLL i+1 is selected this time. If the previous tracking frequency ft p is less than the lower threshold value f i,i-1 , the output of the phase locked loop PLL i-1 is selected this time.
  • the selected phase locked loop is switched to the output of the phase locked loop PLL i+1 , which covers a higher frequency band.
  • the selected phase locked loop is switched to the output of the phase locked loop PLL i-1 , which covers a lower frequency band.
  • Hysteresis is provided to prevent frequent switching near the threshold for switching the output of the phase locked loop.
  • the filter unit BPF2 is 1.8 Hz
  • the upper cutoff frequency f2H is 2.8 Hz
  • the lower threshold f21 for switching the phase locked loop from PLL2 to PLL1 is 1.9 Hz
  • the upper threshold f23 for switching from PLL2 to PLL3 is 2.7 Hz.
  • the control of the output control unit 28 by the input control unit 50 (FIG. 11) will be described.
  • the output of the phase control circuit PLL i selected by the selector 28A is controlled to be fixed.
  • the previous tracking frequency ft p of the delay unit 28B (FIG. 13) is initialized. That is, the previous tracking frequency ft p is controlled to be set to an initial value f init .
  • the signal analysis unit 40 has a function of determining the initial value f init of the previous tracking frequency ftp of the delay unit 28B (FIG. 13) of the output control unit 28. For example, when the subject is at rest and there is no body movement, the signal analysis unit 40 extracts signals from the lower cutoff frequency f1L of the filter unit BPF1 on the lowest frequency side to the upper cutoff frequency f5H of the filter unit BPF5 on the highest frequency side through a bandpass filter. The signal that has passed through the bandpass filter is subjected to a spectrum analysis to detect a spectrum peak. The frequency at which the peak appears is adopted as the initial value f init of the previous tracking frequency ftp of the delay unit 28B.
  • a control may be performed in which a predetermined initial frequency is provided as the initial value f init of the delay section 28B.
  • a pass frequency band covered by a plurality of filter units BPF i as a whole i.e., a frequency band from the lower cutoff frequency f 1L of filter unit BPF 1 to the upper cutoff frequency f 5H of filter unit BPF 5 , is realized by one filter unit.
  • a specific harmonic of the target order n but also harmonics of other orders and noise are likely to be input to the frequency calculation unit 20.
  • the frequency calculation unit 20 is affected by signals other than the harmonic of the target order n, and the frequency calculation accuracy decreases.
  • the bandpass filter 10 ( Figure 11) is realized by a single filter unit with a narrow pass frequency bandwidth.
  • fluctuations in the frequency of the fundamental wave of the biological signal SigB make it easier for the frequency of the harmonic of the target order n to fall outside the pass frequency band of the bandpass filter 10. If the frequency of the harmonic of the target order n falls outside the pass frequency band of the bandpass filter 10, the frequency calculation unit 20 will no longer be able to track the harmonic of the target order n. This narrows the range of measurable heartbeat frequencies.
  • a wide pass frequency band can be covered as a whole by using a plurality of filter units BPF i , thereby widening the range of measurable frequencies of biosignals.
  • the number of filter units BPF i constituting the bandpass filter 10 (FIG. 11) is five, but it may be any other number.
  • FIG. 16 is a block diagram of a bioinformation measuring device according to this modified example.
  • the frequency calculation unit 20 is composed of a plurality of phase synchronization circuits PLL i , and an output control unit 28 is arranged in the subsequent stage of the plurality of phase synchronization circuits PLL i .
  • one phase synchronization circuit 21 is arranged in the subsequent stage of the output control unit 28.
  • the output control unit 28 selects one of the plurality of filter units BPF i , and inputs the first signal Sig1 output from the selected filter unit BPF i to the phase synchronization circuit 21.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Signal Processing (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
PCT/JP2023/041263 2023-02-10 2023-11-16 生体情報計測装置 WO2024166481A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2024576115A JPWO2024166481A1 (enrdf_load_stackoverflow) 2023-02-10 2023-11-16

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2023019368 2023-02-10
JP2023-019368 2023-02-10

Publications (1)

Publication Number Publication Date
WO2024166481A1 true WO2024166481A1 (ja) 2024-08-15

Family

ID=92262864

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2023/041263 WO2024166481A1 (ja) 2023-02-10 2023-11-16 生体情報計測装置

Country Status (2)

Country Link
JP (1) JPWO2024166481A1 (enrdf_load_stackoverflow)
WO (1) WO2024166481A1 (enrdf_load_stackoverflow)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003135434A (ja) * 2001-10-30 2003-05-13 Nippon Koden Corp 信号処理方法および脈波信号処理方法
US20080005838A1 (en) * 2006-07-05 2008-01-10 Stryker Corporation System for detecting and monitoring vital signs
US20090143692A1 (en) * 2007-11-30 2009-06-04 Transoma Medical, Inc. Physiologic Signal Processing To Determine A Cardiac Condition
JP2011194217A (ja) * 2010-03-23 2011-10-06 General Electric Co <Ge> 振動測定法におけるアーティファクトの周波数スペクトルの使用法
JP2012502671A (ja) * 2008-05-12 2012-02-02 アーリーセンス エルティディ 臨床症状のモニタリング、予測及び治療

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003135434A (ja) * 2001-10-30 2003-05-13 Nippon Koden Corp 信号処理方法および脈波信号処理方法
US20080005838A1 (en) * 2006-07-05 2008-01-10 Stryker Corporation System for detecting and monitoring vital signs
US20090143692A1 (en) * 2007-11-30 2009-06-04 Transoma Medical, Inc. Physiologic Signal Processing To Determine A Cardiac Condition
JP2012502671A (ja) * 2008-05-12 2012-02-02 アーリーセンス エルティディ 臨床症状のモニタリング、予測及び治療
JP2011194217A (ja) * 2010-03-23 2011-10-06 General Electric Co <Ge> 振動測定法におけるアーティファクトの周波数スペクトルの使用法

Also Published As

Publication number Publication date
JPWO2024166481A1 (enrdf_load_stackoverflow) 2024-08-15

Similar Documents

Publication Publication Date Title
Kohler et al. QRS detection using zero crossing counts
AU2004236400B2 (en) Operating method for a coriolis gyroscope and evaluation/adjustment electronic system and pulse modulator suitable therefor
JP2011115459A (ja) 生体情報検出装置及び生体情報検出方法
US11183092B2 (en) Capacitance-type touch panel device and operation method for same
US6802815B2 (en) Noninvasive continuous blood pressure measuring apparatus and a method of noninvasively measuring continuous blood pressure
KR20090096803A (ko) 맥박 신호에서 잡음을 제거하기 위한 방법 및 장치와기록매체
JP2016182165A (ja) 生体信号処理装置、生体信号処理プログラム、生体信号処理プログラムを記録したコンピュータ読み取り可能な記録媒体、生体信号処理方法
WO2024166481A1 (ja) 生体情報計測装置
WO2008029587A1 (fr) Dispositif de mesure médicale destiné à la production d&#39;informations utiles pour un diagnostic
JP2017169868A (ja) 心肺機能測定装置
WO2024166483A1 (ja) 生体情報計測装置及び可変フィルタ回路
US9717462B2 (en) Biological information measurement method and apparatus with variable loop filter
JP2003135412A (ja) 血圧計測装置
JP5800776B2 (ja) 生体動情報検出装置
JP2013183845A (ja) 拍動検出装置、電子機器及びプログラム
JP5741887B2 (ja) 周波数測定システム、周波数測定方法、及び周波数測定システムを備える電子機器
JP7596208B2 (ja) 周波数計測器
JP6649037B2 (ja) 検体情報処理装置、情報処理方法、情報処理プログラム、及び同プログラムを記録したコンピュータ読み取り可能な記録媒体
US9928419B2 (en) Periodicity analysis system
JP2005214932A5 (enrdf_load_stackoverflow)
JP2005337980A (ja) 交流信号測定装置
KR101107722B1 (ko) 광대역 디지털 주파수 합성기
JP6213569B2 (ja) 生体情報計測装置
US6425872B1 (en) System for measuring physical parameter utilizing vibration transmission
JPH06165764A (ja) 血管硬化度測定装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23921302

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2024576115

Country of ref document: JP