WO2016084473A1 - 高周波装置 - Google Patents

高周波装置 Download PDF

Info

Publication number
WO2016084473A1
WO2016084473A1 PCT/JP2015/077834 JP2015077834W WO2016084473A1 WO 2016084473 A1 WO2016084473 A1 WO 2016084473A1 JP 2015077834 W JP2015077834 W JP 2015077834W WO 2016084473 A1 WO2016084473 A1 WO 2016084473A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
heart rate
digital signal
biological
respiration
Prior art date
Application number
PCT/JP2015/077834
Other languages
English (en)
French (fr)
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 US15/516,713 priority Critical patent/US10582877B2/en
Priority to JP2016561437A priority patent/JP6290448B2/ja
Priority to CN201580065123.6A priority patent/CN106999105B/zh
Publication of WO2016084473A1 publication Critical patent/WO2016084473A1/ja

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0507Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/0205Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405Determining heart rate variability
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02438Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02444Details of sensor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816Measuring devices for examining respiratory frequency
    • 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
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • 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/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Definitions

  • the present invention relates to a high-frequency device for sensing biological information such as heart rate.
  • the biological signal detection apparatus includes a sensor unit 101, a biological signal extraction unit 102, a distance calculation unit 103, and a biological signal output determination unit 104.
  • FIG. 13 shows a schematic configuration of the sensor unit 101.
  • the signal transmitted from the local oscillator 301 is divided into two by the distributor 302.
  • One signal is transmitted to the transmission antenna 303, the other signal is further divided into two by the distributor 308, and the signal is input to the mixers 306 and 307.
  • the two signals divided by the distributor 308 are 90 degrees out of phase.
  • the signal radiated from the transmission antenna 303 is directed toward the occupant and reflected mainly on the occupant's body surface.
  • the reflected signal reflected from the body surface is input to the receiving antenna 304.
  • the transmission signal to the occupant undergoes a Doppler shift. Therefore, the reflected signal is input to the receiving antenna 304 as a received signal in a state where the frequency is phase-modulated.
  • the received signal input to the receiving antenna 304 is divided into two by the distributor 305 and input to the mixers 306 and 307, respectively.
  • the signal input from the distributor 308 is also input to the mixers 306 and 307.
  • the mixers 306 and 307 perform multiplication processing, and pass the low-pass filters 309 and 310 to output baseband signals subjected to Doppler shift.
  • the baseband signal is analog-digital converted by AD converters 311 and 312 and output as a Bi (t) signal and a Bq (t) signal.
  • the Bi (t) signal and the Bq (t) signal are signals that are instantaneously 90 degrees out of phase.
  • the Bi (t) signal and Bq (t) signal are input to the biological signal extraction unit 102 as shown in FIG.
  • the biological signal extraction unit 102 extracts a biological signal
  • the distance calculation unit 103 calculates an estimated distance between the sensor unit 101 and the occupant.
  • the biological signal output determination unit 104 determines whether or not to output a biological signal based on the estimated distance by setting a certain threshold value. Specifically, when the reliability determination unit 408 in the biological signal output determination unit 104 determines that the distance to the occupant is constant or within a threshold, a biological signal is output.
  • the digitized Bi (t) signal and Bq (t) signal are signals that have undergone Doppler shift, that is, signals that have undergone phase modulation at the Doppler frequency, and can be represented by signals of sine waves and cosine waves, respectively.
  • the phase signal calculation unit 401 calculates the phase change with time, that is, the time waveform of each heartbeat and respiration, from the two signals, and the heartbeat signal extraction unit 402 and the respiratory signal extraction. A process of extraction by the unit 403 is performed.
  • JP 2010-120493 A (published on June 3, 2010)
  • a conventional microwave sensor device is typically the one disclosed in Patent Document 1 described above, and has a mechanism for detecting minute movements of a subject using microwaves and obtaining biological information such as a heartbeat. Specifically, it is intended to detect fine movement of the body surface such as a subject's heartbeat by detecting a phase change (phase difference) of a reflected wave with respect to an irradiation wave.
  • heart rate and respiratory information can be acquired, but the subject is sleeping while lying down (prone,
  • the heart rate / respiration component could not be measured and monitored at all times.
  • the heartbeat signal is a weaker signal than the respiratory signal, it is particularly difficult to accurately extract the heartbeat component.
  • the present invention has been made to solve the above-described problems, and an object of the present invention is to provide a high-frequency device that can wirelessly and accurately detect biological information representing the heartbeat and respiration of the living body. It is in.
  • a high-frequency device includes: A biological signal extraction unit that extracts a biological signal indicating a specific frequency component of the digital signal from a digital signal representing a biological phenomenon that is a movement appearing in the living body; The biological signal extracted by the biological signal extraction unit is sampled, an autocorrelation function is calculated based on the predetermined number of the biological signals sampled, and the periodicity of the autocorrelation function is calculated from the peak value of the correlation coefficient. And an autocorrelation function processing unit for calculating biological information representing the biological phenomenon.
  • a high frequency signal that accurately detects biological information representing heartbeat and respiration using an autocorrelation function.
  • (A) is a side view showing an example of a method of using the high-frequency device according to Embodiments 1 to 4 of the present invention.
  • (B) is a top view showing an example of a method of using the high-frequency device according to Embodiments 1 to 4 of the present invention.
  • It is a block diagram which shows schematic structure of the radio wave radar part with which the high frequency apparatus which concerns on Embodiment 1 and 2 of this invention is provided.
  • It is a block diagram which shows schematic structure of the signal processing circuit with which the high frequency device which concerns on Embodiment 1 of this invention is provided.
  • (A) shows the relationship between the amplitude value of the heartbeat signal and the number of samples regarding the heartbeat signal sampled in the first heartbeat autocorrelation function processing unit and the second heartbeat autocorrelation function processing unit provided in the high-frequency device. It is a graph.
  • (B) is a graph which shows the relationship between the correlation coefficient of an autocorrelation function, and the number of samples regarding the said heart rate signal.
  • It is a block diagram which shows schematic structure of the signal processing circuit with which the high frequency device which concerns on Embodiment 4 of this invention is provided. It is a block diagram which shows schematic structure of the digital signal processing part with which the said high frequency device is equipped, and a display. It is the block diagram which showed schematic structure of the sensor part with which the conventional microwave sensor apparatus is provided. It is the block diagram which showed schematic structure of the conventional microwave sensor apparatus.
  • FIG. 1 is a side view which shows an example of the usage method of the high frequency device 1 which concerns on this embodiment.
  • FIG. 1B is a top view showing an example of how to use the high-frequency device 1 according to the present embodiment.
  • the high-frequency device 1 is a device for obtaining biological information representing a biological phenomenon of the living body 10 based on a reflected signal 12 (see FIG. 2) in which the transmission signal 11 is reflected by the living body 10.
  • the biological phenomenon is a movement that appears mainly in the upper body of the living body 10, and indicates a heartbeat component, a respiratory component, a body motion component, and the like.
  • the body movement refers to movements that appear on the upper body of the living body 10 other than heartbeat and respiration.
  • examples of the biological information include a heart rate and a respiration rate for a predetermined time, a heart rate waveform, a respiration waveform, and a body motion waveform.
  • the high-frequency device 1 includes a radio wave radar unit 5a and a signal processing circuit 400 (see FIG. 2).
  • the transmission signal 11 (see FIG. 2) radiated from the high-frequency device 1 is desirably radiated to the upper half of the living body 10, so that the transmission antenna 25 (
  • the directivity pattern 93 (see FIG. 2) is adjusted. As an example of the adjustment, it can be adjusted by the number of patches of the microstrip patch antenna.
  • the reflected signal 12 becomes a received signal of the high frequency device 1.
  • the high-frequency device 1 is attached to the handrail 91 b of the bed 91, so that the transmission signal 11 is directed substantially parallel to the contact surface of the mattress 91 a of the bed 91 with the living body 10. Radiate. Therefore, it is possible to always measure and monitor the heart rate component, the respiratory component, and the body motion component with respect to the sleeping posture of the living body 10 (prone, lying on its back, sideways, moving its body while sleeping, etc.).
  • the transmission signal 11 is irradiated in the substantially horizontal direction, so that the transmission signal 11 and the reflection signal 12 are irradiated. This is because there are fewer. That is, unnecessary reflection signals received by the high-frequency device 1 are reduced.
  • the high-frequency device 1 processes a heartbeat signal representing a heartbeat component and a respiratory signal representing a respiratory component using an autocorrelation function described later. Therefore, the high-frequency device 1 can measure and monitor a heartbeat component and the like from an indirect reflection signal described later even in a bedroom where there are many unnecessary reflection signals.
  • the attachment position of the high frequency apparatus 1 it is not limited to the handrail 91b of the bed 91 like this embodiment.
  • the reflection component of the transmission signal 11 with respect to the mattress 91a is small, and therefore, substantially the same effect as when the high-frequency device 1 is mounted as in the present embodiment can be obtained.
  • the method of using the high frequency device 1 in the present embodiment is similarly adopted for the high frequency devices 2 to 4 described later.
  • FIG. 2 is a block diagram showing a schematic configuration of the radio wave radar unit 5a.
  • the radio wave radar unit 5 a radiates the transmission signal 11 toward the living body 10, receives the reflected signal 12 reflected by the body surface of the living body 10, and generates an analog signal 33 s based on the reflected signal 12.
  • the radio wave radar unit 5 a includes an oscillator 21, an amplifier 22, a transmission antenna 25, a reception antenna 30, a low noise amplifier 31, a mixer 32, and a filter 33.
  • a microwave sine wave (not shown) output from the oscillator 21 is amplified by the amplifier 22 and then radiated from the transmission antenna 25 as the transmission signal 11.
  • the transmitted transmission signal 11 hits the body surface of the living body 10, for example, the chest surface of the living body 10, and is reflected.
  • the living body 10 has chest movement due to heartbeat, respiration, and body movement, and this movement causes a Doppler shift in the transmission signal 11 that hits the chest. That is, by adding the frequency of the chest motion including heartbeat and respiration to the frequency of the transmission signal 11, the transmission signal 11 is frequency / phase-modulated to become the reflected signal 12. Then, the reflected signal 12 is input to the receiving antenna 30.
  • the reflected signal 12 is received by the receiving antenna 30 through direct reflection on the living body 10 and directly received by the receiving antenna 30 after reflection and reflection from an object other than the living body 10. Both indirect reflected signals are included.
  • the reflected signal 12 input to the receiving antenna 30 is amplified by the low noise amplifier 31, and then input to the mixer 32 as an analog signal 31s.
  • the mixer 32 is also supplied with one analog signal 22 s which is transmitted from the oscillator 21 and amplified by the amplifier 22 and then divided into two.
  • the analog signal 31s is input to the filter 33 in a state where the frequency is down-converted by the analog signal 22s.
  • the analog signal 33 s output from the filter 33 is a baseband signal, and is output as a signal that has undergone Doppler shift due to chest movement. That is, the analog signal 33s includes the heartbeat component, the respiratory component, and the body motion component of the living body 10.
  • the analog signal 33s is input to a signal processing circuit 400 described later and subjected to signal processing, and then displayed on a display 90 described later as heartbeat information, respiratory information, and body movement information.
  • FIG. 3 is a block diagram showing a schematic configuration of the signal processing circuit 400.
  • FIG. 4 is a block diagram showing a schematic configuration of a digital signal processing unit (hereinafter referred to as “DSP”) 500 and a display 90 described later.
  • DSP digital signal processing unit
  • the signal processing circuit 400 Based on the input analog signal 33s, the signal processing circuit 400, the heart rate information (heart rate and heart rate waveform per minute) of the living body 10, respiratory information (respiration rate and respiratory waveform per minute), and body motion information. Calculate and determine (body motion waveform).
  • the signal processing circuit 400 includes a digital signal generation unit 40 a and a DSP 500.
  • the digital signal generator 40a generates a digital signal 46s including a heartbeat component, a respiratory component, and a body motion component of the living body 10 based on the input analog signal 33s.
  • the digital signal generation unit 40 a includes an input terminal 41, a high-pass filter 43, a low-pass filter 44, an amplifier 45, and an AD converter 46.
  • the DSP 500 calculates and determines the heartbeat information, the respiratory information, and the body motion information based on the input digital signal 46s.
  • the DSP 500 includes an input terminal 51, a heartbeat signal extraction unit (biological signal extraction unit) 53, a heartbeat autocorrelation function processing unit (autocorrelation function processing unit) 54, a heart rate determination unit 55, and a respiratory signal extraction.
  • the analog signal 33 s (see FIG. 2) is input to the input terminal 41 and then band-limited by the high-pass filter 43 (for the purpose of suppressing the dc component) and the low-pass filter 44.
  • the high-pass filter 43 has a pass-band cutoff of 0.1 Hz
  • the low-pass filter 44 has a pass-band cutoff of 10 Hz, which passes analog signals in frequency bands related to heartbeat, breathing, and body movement. Then, the signal is amplified by the amplifier 45.
  • the amplifier 45 prevents a high-speed analog signal that has passed through the high-pass filter 43 and the low-pass filter 44 from performing a saturation operation, and consequently a non-linear operation.
  • the low-pass filter 44 also functions as an anti-alias filter for the AD converter 46.
  • An analog signal (not shown) output from the amplifier 45 is input to the AD converter 46 and subjected to analog-digital conversion.
  • analog-to-digital conversion is performed at a sampling rate of 12 bits and 50 Hz
  • the cutoff of the pass band in the low-pass filter 44 is 25 Hz or less (in this embodiment, 10 Hz).
  • the digital signal 46 s generated by analog-digital conversion by the AD converter 46 is input to a digital signal processing unit (hereinafter referred to as “DSP”) 500. Then, after the heart rate information, respiration information, and body movement information are calculated by the DSP 500, these information are displayed on the display 90.
  • DSP digital signal processing unit
  • the digital signal 46s (see FIG. 3) input to the input terminal 51 is divided into a first digital signal 58a, a second digital signal 58b, and a third digital signal 58c. Two digital signals are processed in parallel.
  • the first digital signal 58 a is input to the heartbeat signal extraction unit 53.
  • the input first digital signal 58a is band-limited by the high-pass filter 53a and the low-pass filter 53b included in the heartbeat signal extraction unit 53, so that a heartbeat signal indicating a heartbeat of a specific frequency is extracted.
  • the high-pass filter 53a has a passband cutoff frequency of 0.8 Hz
  • the low-pass filter 53b has a passband cutoff frequency of 3.0 Hz.
  • the respiratory component is 0.2 Hz to 0.33 Hz (about 12 to 20 times per minute) at the time of calm, and therefore the pass band of the high-pass filter 53a is cut to suppress the respiratory component as much as possible.
  • An off frequency of 0.8 Hz is selected.
  • the low-pass filter 53b has a heartbeat component of about 50 to 90 (frequency is 0.8 Hz to 1.5 Hz) at the time of calm, so a filter with a passband cutoff frequency of 3.0 Hz is selected.
  • the heartbeat signal extracted by the heartbeat signal extraction unit 53 is divided into two, and one of the signals is input to the heartbeat autocorrelation function processing unit 54. Then, the heartbeat autocorrelation function processing unit 54 performs an arithmetic process using the autocorrelation function, thereby calculating a heart rate per minute (a value of a biological phenomenon; hereinafter, simply referred to as “heart rate”).
  • the heartbeat autocorrelation function processor 54 includes a sampling processor 54a, a heartbeat autocorrelation function calculator 54b, and a peak detector 54c.
  • the sampling processing unit 54a samples the heartbeat signal 256 times in one sample 20 ms from the first digital signal 58a. .
  • the heartbeat autocorrelation function calculator 54b calculates an autocorrelation function by inputting the sampled heartbeat signal.
  • the peak detector 54c detects the peak value of the correlation coefficient from the autocorrelation function, and calculates the period of the autocorrelation function. In other words, the periodicity of the autocorrelation function is determined by the peak detector 54c.
  • the peak detection unit 54c converts the calculated cycle into a heart rate. In the autocorrelation function, since the data after 0 second has the largest correlation coefficient, the peak value of the autocorrelation coefficient is detected except for the data at 0 second.
  • the sampling processing by the sampling processing unit 54a is 5 seconds (256 samplings per 20 ms)
  • the autocorrelation function is calculated by the heartbeat autocorrelation function calculating unit 54b
  • the period is calculated by the peak detecting unit 54c. It takes about 5 seconds. Therefore, one heart rate is calculated in about 10 seconds.
  • the heart rate calculated by the heart rate autocorrelation function processing unit 54 is input to the heart rate determining unit 55 (not shown) as a signal indicating the heart rate.
  • the heart rate determination unit 55 determines a heart rate to be displayed on the display device 90 (specifically, a heart rate display unit 90a described later) based on the signal, and the heart rate indicating the determined heart rate
  • the signal 55s is output.
  • the heart rate signal 55s becomes an input signal to the display 90.
  • the determination of the heart rate in the heart rate determination unit 55 is not limited to one method.
  • one signal is output from the heartbeat autocorrelation function processing unit 54 every about 10 seconds, for example, calculation for obtaining an average value, a moving average value, a median value, or the like for every three signals. May be made.
  • any method that outputs a stable numerical value as the heart rate displayed on the display 90 may be used.
  • the other divided into two is a heartbeat waveform signal 71a representing a change in heartbeat amplitude with the passage of time, and is an input signal to the display 90.
  • the second digital signal 58b is a body motion waveform signal 71b representing a change in body motion amplitude with the passage of time, and is an input signal to the display 90.
  • the third digital signal 58 c is input to the respiratory signal extraction unit 63. Then, the input third digital signal 58c is band-limited by the low-pass filter 63a included in the respiration signal extraction unit 63, so that a respiration signal indicating respiration at a specific frequency is extracted.
  • the low-pass filter 63a has a passband cutoff frequency of 0.75 Hz. This is because the breathing component is 0.2 Hz to 0.33 Hz (about 12 to 20 times per minute) at the time of calm, and thus the breathing signal extraction unit 63 suppresses the heart rate component as much as possible.
  • the respiratory signal extracted by the respiratory signal extraction unit 63 is divided into two, and one of the signals is input to the respiratory autocorrelation function processing unit 64. Then, the respiratory autocorrelation function processing unit 64 performs an arithmetic process using the autocorrelation function, thereby calculating a respiratory rate per minute (a value of a biological phenomenon; hereinafter, simply referred to as “respiration rate”).
  • the respiratory autocorrelation function processor 64 includes a sampling processor 64a, a respiratory autocorrelation function calculator 64b, and a peak detector 64c.
  • the sampling processing unit 64a samples the respiratory signal 512 times in 20 ms per sample from the third digital signal 58c. .
  • the respiratory autocorrelation function calculation unit 64b calculates an autocorrelation function by inputting the sampled respiratory signal.
  • the peak detector 64c detects the peak value of the correlation coefficient from the autocorrelation function, and calculates the period of the autocorrelation function. In other words, the periodicity of the autocorrelation function is determined by the peak detector 64c.
  • the peak detection unit 64c converts the calculated cycle into a respiration rate.
  • the sampling processing by the sampling processing unit 64a is 10 seconds (512 samplings per 20 ms)
  • the autocorrelation function is calculated by the respiratory autocorrelation function calculating unit 64b
  • the period is calculated by the peak detecting unit 64c. It takes about 5 seconds. Therefore, one respiration rate is calculated in about 15 seconds.
  • Embodiment 3 The details of the method for determining the periodicity related to the respiratory component using the autocorrelation function will be described in Embodiment 3.
  • the respiration rate calculated by the respiration autocorrelation function processing unit 64 is input to the respiration rate determination unit 65 (not shown) as a signal indicating the respiration rate. Based on the signal, the respiration rate determination unit 65 determines a respiration rate to be displayed on the display device 90 (specifically, a respiration rate display unit 90e described later), and the respiration rate indicating the determined respiration rate The signal 65s is output.
  • the respiration rate signal 65s is an input signal to the display 90.
  • the determination of the respiration rate in the respiration rate determination unit 65 is not limited to one method.
  • one signal is output from the respiratory autocorrelation function processing unit 64 about every 15 seconds, for example, calculation for obtaining an average value, a moving average value, a median value, or the like for every three signals. May be made.
  • any method that outputs a stable numerical value as the respiratory rate displayed on the display 90 may be used.
  • the other of the two distributed signals becomes a respiratory waveform signal 71c representing a change in respiratory amplitude with the passage of time, and becomes an input signal to the display 90.
  • the display device 90 includes a heart rate display unit 90a, a heart rate waveform display unit 90b, a body motion waveform display unit 90c, a respiration waveform display unit 90d, and a respiration rate display unit 90e.
  • Each signal output from the DSP 500 is input to the display 90.
  • the heart rate signal 55s is displayed on the heart rate display unit 90a
  • the heart rate waveform signal 71a is displayed on the heart rate waveform display unit 90b
  • the body motion waveform signal 71b is displayed on the body motion waveform display unit 90c
  • the respiration waveform signal 71c is displayed on the respiration waveform.
  • the respiration rate signal 65s is input to the respiration rate display unit 90e.
  • the heart rate display unit 90a has a heart rate (numerical value)
  • the heart rate waveform display unit 90b has a heart rate waveform (graph)
  • the body motion waveform display unit 90c has a body motion waveform (graph)
  • a respiratory waveform display unit is displayed in 90d
  • a respiratory rate (numerical value) is displayed in the respiratory rate display section 90e.
  • Examples of the display device 90 include a display device and a data storage function such as a personal computer, a mobile phone, or a smartphone.
  • the DSP 500 and the display device 90 may be integrated.
  • the digital signal generation unit 40a wirelessly transmits a signal that has been converted from analog to digital, and a personal computer or the like equipped with the DSP 500 and the display device 90 receives the signal, thereby performing signal processing and display by the personal computer or the like. May be.
  • the high-frequency device 1 extracts only a band necessary for the analog filter (the high-pass filter 43 and the low-pass filter 44) and the digital filter (the high-pass filter 53a, the low-pass filters 53b and 63a). But pass only the necessary band.
  • the high-frequency device 1 increases the dynamic range in the heartbeat or breathing band by extracting only a necessary band with an analog filter, removing unnecessary signals, amplifying, and then performing analog-digital conversion. Then, the SN ratio of the digital signal is increased in the digital filter. Therefore, the detection accuracy of the peak value of the correlation coefficient detected from the autocorrelation function can be increased.
  • the high-frequency device 1 detects the peak value of the correlation coefficient using the autocorrelation function obtained by sampling the heartbeat signal and the respiratory signal, and determines the period of the autocorrelation function, thereby determining the heart rate and Calculate respiratory rate. Therefore, even when the living body 10 is moving or there are many noise components such as unnecessary reflection signals, the heart rate and the respiration rate can be calculated with high accuracy.
  • the heartbeat signal and the respiration signal can be extracted from the direct reflection signal and the indirect reflection signal using the heartbeat signal extraction unit 53 and the respiration signal extraction unit 63. Therefore, for the indirect reflection signal, the heart rate, the respiration rate, and the like can be calculated using the autocorrelation function of the heart rate signal and the respiration signal.
  • the high-frequency device 1 can calculate the heart rate and respiration rate of the living body 10 existing in the space using the indirect reflection signal even when the direct reflection signal cannot be captured.
  • the transmission signal 11 when the living body 10 exists in the space, for example, even when it is difficult to directly irradiate the desired part of the living body 10 with the transmission signal 11 such as the living body 10 lying down.
  • the transmission signal 11 By reflecting the transmission signal 11 with a wall or the like, the transmission signal 11 can be indirectly irradiated onto the part. Therefore, by using the indirect reflection signal, the heart rate and the respiration rate can be calculated regardless of the posture of the living body 10.
  • the high frequency device 2 according to the present embodiment is different from the high frequency device 1 according to the first embodiment in that a signal processing circuit 401 is provided instead of the signal processing circuit 400.
  • FIG. 5 is a block diagram showing a schematic configuration of the signal processing circuit 401.
  • FIG. 6 is a block diagram showing a schematic configuration of the DSP 501 and the display device 90.
  • the signal processing circuit 401 includes a digital signal generation unit 40 b and a DSP 501.
  • the digital signal generation unit 40b generates a digital signal 46sa including a heartbeat component and a body motion component of the living body 10 and a digital signal 46sb including a respiratory component of the living body 10 based on the input analog signal 33s (see FIG. 2). To do.
  • the digital signal generator 40b includes an input terminal 41, high-pass filters 43a and 43b, low-pass filters 44a and 44b, amplifiers 45a and 45b, and AD converters 46a and 46b.
  • the terminal 48 distributes the analog signal 33s into the first analog signal 48a and the second analog signal 48b.
  • the first analog signal 48a is limited and amplified in the band for heart rate and body movement.
  • the cut-off of the pass band in the high-pass filter 43a is 0.75 Hz
  • the cut-off of the pass band in the low-pass filter 44a is 10 Hz.
  • Amplification is performed by the amplifier 45a. The reason why the band limitation by the high-pass filter 43a is as described above is to suppress the respiratory component as much as possible, as in the first embodiment.
  • An analog signal (not shown) output from the amplifier 45a is input to the AD converter 46a and subjected to analog-digital conversion.
  • the first biological digital signal 46sa relating to the heartbeat component and the body motion component generated by the analog-digital conversion by the AD converter 46a is input to the DSP 501.
  • the movement of the chest surface during breathing is usually 4 mm to 12 mm, and the movement of the chest surface during heartbeat is about 0.5 mm, so the analog signal related to heartbeat is smaller than the analog signal related to breathing. It is.
  • the amplifier 45a is prevented from being saturated by movement of the chest surface during breathing, and only the heartbeat component is analog-digital converted in the AD converter 46a. be able to. Therefore, the dynamic range of the amplitude of the heartbeat waveform over time can be increased, and the sensitivity and accuracy of digital signal processing by the DSP 501 can be improved.
  • only analog signals of a heartbeat component and a body motion component having a small amplitude can be converted into an analog-digital signal by 16 bits, and the effect is particularly exerted at rest or sleep when the body motion component is small.
  • the second analog signal 48b is limited to the breathing band and amplified.
  • the cutoff of the pass band in the high pass filter 43b is 0.1 Hz
  • the cutoff of the pass band in the low pass filter 44b is 0.75 Hz
  • the amplifier Amplify at 45b after passing the signal in the frequency band related to the respiratory component, the amplifier Amplify at 45b.
  • An analog signal (not shown) output from the amplifier 45b is input to the AD converter 46b and subjected to analog-digital conversion.
  • the second biological digital signal 46 sb relating to the respiratory component generated by analog-digital conversion by the AD converter 46 b is input to the DSP 501.
  • the high-pass filter 43b suppresses the dc component
  • the low-pass filter 44b suppresses the heart rate component and the body motion component.
  • the sensitivity and accuracy of the digital signal processing by the DSP 501 can be improved by increasing the dynamic range of the respiratory component signal obtained from the second analog signal 48b.
  • only analog signals of respiratory components can be converted from analog to digital in 16 bits, which is effective when the amplitude of body motion components is large or when there is body motion before going to sleep even during sleep. .
  • the first biological digital signal 46sa (see FIG. 5) is inputted to the input terminal 51, and the second biological digital signal 46sb (see FIG. 5) is inputted to the input terminal 61, respectively.
  • the first biological digital signal 46sa input to the input terminal 51 is distributed to the heartbeat digital signal 59a related to the heartbeat component and the body motion digital signal 59b related to the body motion component, and these two signals are processed in parallel.
  • the signal processing for obtaining the heart rate signal 55s and the heart rate waveform signal 71a based on the heart rate digital signal 59a is the same as that of the first embodiment, the description thereof is omitted.
  • the second biological digital signal 46sb input to the input terminal 62 is directly processed as a respiratory digital signal 59c.
  • the signal processing for obtaining the respiration rate signal 65s and the respiration waveform signal 72c2 based on the respiration digital signal 59c is the same as that in the first embodiment, and thus the description thereof is omitted.
  • the body motion digital signal 59b is added with the branch signal 71c1 of the respiratory signal extracted by the respiratory signal extraction unit 63 in order to add the chest movement of the living body 10 related to the respiratory component as the body motion component. . Then, it is output from the DSP 501 as a body motion waveform signal 71b that expresses a wide range from slow body motion to fast body motion.
  • the digital signal generation unit 40b can increase the dynamic range for each of the first biological digital signal 46sa and the second biological digital signal 46sb. Therefore, the heart rate signal extraction unit 53 and the respiration signal extraction unit 63 can extract the heart rate signal and the respiration signal with high sensitivity and accuracy.
  • the digital signal generation unit 40b suppresses the respiratory component of the first analog signal 48a even during sleep of a living body having a small amplitude of the body motion component (that is, when the influence of the respiratory component is large). By doing so, it is possible to generate the first biological digital signal 46sa having an increased dynamic range. Further, even when the living body has a large amplitude of the body motion component (that is, when the influence of the body motion component is large), by suppressing the body motion component of the second analog signal 48b in the digital signal generation unit 40b, The second biological digital signal 46sb with an increased dynamic range can be generated.
  • the digital signal generation unit 40b generates the first biological digital signal 46sa related to the heartbeat component and the body motion component, and the second biological digital signal 46sb related to the respiratory component.
  • the signal generation in the digital signal generation unit 40b is not limited to the above case.
  • the digital signal generation unit 40b receives three digital signals: a digital signal related to a heartbeat component and a body motion component having a speed equivalent to the heartbeat component, a digital signal related to a body motion component moving faster than the heartbeat component, and a respiratory component signal. It may be generated. Also, four or more digital signals may be generated.
  • the DSP 501 and the display device 90 may be integrated.
  • the digital signal generation unit 40b wirelessly transmits a signal that has been converted from analog to digital, and a personal computer or the like that includes the DSP 501 and the display 90 receives the signal, thereby performing signal processing and display by the personal computer or the like. May be.
  • the high-frequency device 3 according to the present embodiment is different from the first and second embodiments in that a radio wave radar unit 5b is provided instead of the radio wave radar unit 5a, and a signal processing circuit 402 is provided instead of the signal processing circuits 400 and 401. Different from the high-frequency devices 1 and 2.
  • FIG. 7 is a block diagram showing a schematic configuration of the radio wave radar unit 5b.
  • the radio wave radar unit 5b includes an oscillator 21, an amplifier 22, a transmission antenna 25, a reception antenna 30, a low noise amplifier 31, an I mixer 32i, a Q mixer 32q, a filter 33i, a filter 33q, and a phase shifter 38. I have.
  • the signal processing up to the point where the transmission signal 11 is radiated from the transmission antenna 25 and the reflected signal 12 is input to the reception antenna 30 is the same as in the first and second embodiments.
  • the point that the reflected signal 12 includes both a direct reflected signal and an indirect reflected signal is the same as in the first and second embodiments.
  • the signal input to the receiving antenna 25 is amplified by the low noise amplifier 31 and then divided into two, an I analog signal 31is and a Q analog signal 31iq.
  • the I analog signal 31is is input to the phase shifter 38, so that the phase is shifted by 90 degrees with respect to the Q analog signal 31iq.
  • the I analog signal 31is output from the phase shifter 38 is input to the I mixer 32i, and the Q analog signal 31iq is input to the Q mixer 32q.
  • one of the two analog signals 22 s divided into two is further divided into the I-side local oscillation signal 22 is and the Q-side local oscillation signal 22 qs, and the I mixers 32 i and Q Each is input to the mixer 32q.
  • the phases of the I analog signal 31is and the Q analog signal 31iq are shifted by 90 degrees, but the phases of the I-side local oscillation signal 22is and the Q-side local oscillation signal 22qs are shifted by 90 degrees. It doesn't matter.
  • the I analog signal 31is output from the phase shifter 38 is frequency down-converted by the I mixer 32i and the Q analog signal 31iq by the Q mixer 32q, respectively, and then input to the filter 33i and the filter 33q. Then, an I baseband signal 33is is output from the filter 33i, and a Q baseband signal 33qs is output from the filter 33q.
  • the I baseband signal 33is is a cosine wave, and the Q baseband signal 33qs is a sine wave.
  • Both the I baseband signal 33is and the Q baseband signal 33qs are analog signals that have undergone Doppler shift due to chest movements including heartbeat and respiration.
  • the difference between the I baseband signal 33is and the Q baseband signal 33qs is whether or not the signal has passed through the phase shifter 38.
  • both signals are 90 degrees out of phase. For this reason, the phase relationship between the I baseband signal 33is and the Q baseband signal 33qs as the integrated signal output from the radio wave radar unit 5b constantly varies in accordance with the magnitude and direction of the signal speed.
  • the autocorrelation function is calculated from the two channels for the I series and the Q series, respectively, when the living body 10 has a delicate body motion, the I baseband signal 33is and the Q baseband.
  • the amplitude value of the signal 33qs with respect to time is likely to change. Therefore, for example, the Q baseband signal 33qs may be a weak signal compared to the I baseband signal 33is. Therefore, only the detection of the peak value by the autocorrelation function of one of the channels will deteriorate the detection sensitivity and detection accuracy of the peak value.
  • FIG. 8 is a block diagram showing a schematic configuration of the signal processing circuit 402.
  • FIG. 9 is a block diagram showing a schematic configuration of the DSP 502 and the display device 90.
  • the signal processing circuit 402 includes a first digital signal generation unit (digital signal generation unit) 40ci, a second digital signal generation unit (digital signal generation unit) 40cq, and a DSP 502.
  • the I baseband signal 33is (see FIG. 7) is input to the input terminal 41i, and the Q baseband signal 33qs (see FIG. 7) is input to the input terminal 41q.
  • the input I baseband signal 33is is processed by the I-sequence first digital signal generator 40ci, and the input Q baseband signal 33qs is processed by the Q-sequence second digital signal generator 40cq. Since the first and second digital signal generation units 40ci and 40cq have the same configuration as the digital signal generation unit 40a (see FIG. 3) for the I series and the Q series, the first and second digital signal generation units 40ci and 40cq A description of signal processing in the digital signal generation units 40ci and 40cq is omitted.
  • the I digital signal 46si output from the first digital signal generation unit 40ci and the Q digital signal 46sq output from the second digital signal generation unit 40cq are respectively body motion components (including a heartbeat component and a respiratory component). After being limited to the frequency band and output as an amplified digital signal, the DSP 502 processes the signal.
  • the I digital signal 46si is a cosine wave
  • the Q digital signal 46sq is a sine wave.
  • the DSP 502 includes input terminals 52i and 52q, a first heartbeat signal extraction unit (heartbeat signal extraction unit) 53i, a second heartbeat signal extraction unit (heartbeat signal extraction unit) 53q, and a first respiratory signal extraction unit.
  • Respiration signal extraction unit 63i second respiration signal extraction unit (respiration signal extraction unit) 63q, first heart rate autocorrelation function processing unit (heart rate autocorrelation function processing unit) 54i, second heart rate autocorrelation function processing unit (heart rate) Autocorrelation function processor) 54q, first respiratory autocorrelation function processor (respiratory autocorrelation function processor) 64i, second respiratory autocorrelation function processor (respiratory autocorrelation function processor) 64q, heart rate determination unit 55a, A respiration rate determining unit 65a is provided.
  • the I digital signal 46si (see FIG. 8) input to the input terminal 52i is 3 into the first I digital signal 58ai, the second I digital signal 58bi, and the third I digital signal 58ci. Distributed.
  • the second I digital signal 58bi is output as it is as a body motion waveform signal 71b representing a change in body motion amplitude over time, and becomes an input signal to the display 90.
  • the Q digital signal 46sq (see FIG. 8) input to the input terminal 52q is divided into two parts, the first Q digital signal 58aq and the second Q digital signal 58bq.
  • the first heart rate that is the heart rate calculated by the first heart rate autocorrelation function processing unit 54i and the second heart rate that is the heart rate calculated by the second heart rate autocorrelation function processing unit 54q are: Each is input to the heart rate determination unit 55a (not shown) as a signal indicating the heart rate. Specifically, the first heart rate is input to the first heart rate determining unit 55i provided in the heart rate determining unit 55a. Also, the second heart rate is input to the second heart rate determination unit 55q provided in the heart rate determination unit 55a.
  • the method for determining the first heart rate and the second heart rate is the same as that for the heart rate determination unit 55 described in the first embodiment, and thus the description thereof is omitted.
  • the first heart rate and the second heart rate are the same in the state where the living body 10 is calm, but when the balance between the I-sequence channel and the Q-sequence channel is poor, that is, when the body movements, etc. In some cases, one channel may show abnormal values. Therefore, as a method for determining the heart rate for display in the display heart rate determination unit 55b provided in the heart rate determination unit 55a, for example, the numerical values of the first heart rate and the second heart rate are within an error of several percent. If so, these average values may be used as the heart rate. Alternatively, the first heart rate and the second heart rate may be compared, and the higher value may be used as the heart rate.
  • the first respiration rate which is the respiration rate calculated by the first respiration autocorrelation function processing unit 64i
  • the second respiration rate which is the respiration rate calculated by the second respiration autocorrelation function processing unit 64q
  • the respiration rate determination unit 65a (not shown) as a signal indicating the respiration rate.
  • the first respiratory rate is input to the first respiratory rate determining unit 65i provided in the respiratory rate determining unit 65a.
  • the second respiratory rate is also input to the second respiratory rate determining unit 65q provided in the respiratory rate determining unit 65a.
  • the description is abbreviate
  • the method of determining the display respiratory rate in the display respiratory rate determining unit 65b provided in the respiratory rate determining unit 65a is the same as the method of determining the display heart rate, and the description thereof is omitted. To do.
  • FIG. 10A is a graph showing the relationship between the amplitude value of the heartbeat signal and the number of samples regarding the heartbeat signal sampled by the first heartbeat autocorrelation function processing unit 54i and the second heartbeat autocorrelation function processing unit 54q. It is.
  • FIG. 10B shows the relationship between the correlation coefficient of the autocorrelation function and the number of samples for the heartbeat signals sampled by the first heartbeat autocorrelation function processing unit 54i and the second heartbeat autocorrelation function processing unit 54q. It is a graph.
  • a biological signal such as a heartbeat signal and a respiration signal can be regarded as a signal having a finite average power, and thus the signal is expressed as x [n].
  • n is a sample number by sampling and is an integer.
  • the autocorrelation function r xx [l] is calculated by the equation (1).
  • Equation (1) N1 is an integer sufficiently large with respect to the assumed period, and Equation (1) is
  • the vertical axis is the amplitude value
  • the vertical axis is the correlation coefficient
  • the horizontal axis represents the number of samples of the heartbeat signal. A larger value of both the amplitude value and the correlation coefficient indicates that the correlation with the past signal is higher.
  • the sampling frequency is 50 Hz as in the first and second embodiments, the sampling frequency is 20 milliseconds per sample.
  • the periodicity is calculated from the peak values of the two correlation coefficients using two autocorrelation functions for I series and Q series. Determine. In general, since breathing is performed about once every 4 seconds, the breathing signal needs to be sampled about 512 times, which is different from the case of the heartbeat signal.
  • the high-frequency device 3 finally determines the heart rate to be displayed on the display 90 using the first heart rate and the second heart rate. For example, if the error between the numerical values of the first heart rate and the second heart rate is within a few percent, it is assumed that the living body 10 is stationary (that is, the living body 10 is in a state suitable for calculating the heart rate). Then you can determine your heart rate. A similar method can be used for determining the respiratory rate. Therefore, the high frequency device 3 can finally display the heart rate and the respiration rate with high accuracy on the display 90.
  • the high-frequency device 3 determines the first heart rate and the second heart rate based on two digital signals (I digital signal 46si and Q digital signal 46sq) that are 90 degrees out of phase.
  • the I digital signal 46si which is a cosine wave
  • the Q digital signal 46sq is a signal having a large amplitude because it is a sine wave
  • the two signals compensate for the phase change. It becomes. Therefore, even when there is a relatively large body movement in the living body 10, at least one of the heart rates is kept highly accurate. The same applies to the relationship between the first respiratory rate and the second respiratory rate. Therefore, the high-frequency device 3 can calculate a heart rate and a respiratory rate with high accuracy even when the living body 10 has a relatively large body movement.
  • the levels of the sampled I-sequence heartbeat signal and the Q-sequence heartbeat signal are different (the Q-sequence heartbeat signal is used for the I-sequence.
  • the signal level is lower than the heart rate signal). Therefore, for example, in the discrimination of the periodicity of the heart rate component and the respiration component based only on the autocorrelation function of the Q series heart rate signal, the accuracy of the discrimination is lowered.
  • an autocorrelation function is used.
  • the high-frequency device 4 according to the present embodiment is different from the high-frequency devices 1 to 3 according to the first to third embodiments in that a signal processing circuit 403 is provided instead of the signal processing circuits 400 to 402.
  • FIG. 11 is a block diagram illustrating a schematic configuration of the signal processing circuit 403.
  • FIG. 12 is a block diagram showing a schematic configuration of the DSP 503 and the display device 90.
  • the signal processing circuit 403 includes an I heartbeat digital signal generation unit 49a, an I breathing digital signal generation unit 49b, a Q heartbeat digital signal generation unit 49c, a Q breathing digital signal generation unit 49d, and a DSP 503. .
  • the I baseband signal 33is (see FIG. 7) is input to the input terminal 41b, and the Q baseband signal 33qs (see FIG. 7) is input to the input terminal 41c.
  • the input I baseband signal 33is is divided into two I-line I heart rate analog signals 48ia and I breath analog signals 48ib.
  • the input Q baseband signal 33iq is divided into two Q-sequence Q heart rate analog signals 48qa and Q respiration analog signals 48qb.
  • the I heart rate analog signal 48ia and the Q heart rate analog signal 48qa are analog signals including a heart rate component and a body motion component
  • the I breath analog signal 48ib and the Q breath analog signal 48qb are analog signals including a breath component.
  • the I heart rate analog signal 48ia is processed by the I heart rate digital signal generator 49a, and the I breath analog signal 48ib is processed by the I breath digital signal generator 49b.
  • the Q heart rate analog signal 48qa is signal-processed by the Q heart rate digital signal generation unit 49c, and the Q respiration analog signal 48qb is signal-processed by the Q respiration digital signal generation unit 49d.
  • Each of the digital signal generation units 49a to 49d has the same configuration as that of the digital signal generation unit 40a (see FIG. 3) for the I series and the Q series. The description of the signal processing is omitted.
  • the I heart rate digital signal 46 sia output from the I heart rate digital signal generation unit 49 a and the I respiration digital signal 46 sib output from the I respiration digital signal generation unit 49 b are each processed by the DSP 503.
  • the Q heart rate digital signal 46 sqa output from the Q heart rate digital signal generation unit 49 c and the Q respiration digital signal 46 sqb output from the Q respiration digital signal generation unit 49 d are also processed by the DSP 503.
  • the digital signals 46sia and 46sib for I series are cosine waves, and the digital signals 46sqa and 46sqb for Q series are sine waves.
  • the I heartbeat digital signal 46sia is input to the input terminal 52i. Then, after being output from the input terminal 52i, it is divided into two parts, the first I digital signal 58ai and the second I digital signal 58bi. Thereafter, the first I digital signal 58ai is input to the first heartbeat signal extraction unit 53i. The second I digital signal 58bi is output as it is as a body motion waveform signal 71b representing a change in body motion amplitude over time, and becomes an input signal to the display 90.
  • the I respiration digital signal 46sib is input to the input terminal 62i, and then output to the first respiration signal extraction unit 63i as the third I digital signal 58ci.
  • the Q heartbeat digital signal 46sqa is input to the input terminal 52q and then output to the second heartbeat signal extraction unit 53q as the first Q digital signal 58aq.
  • the Q respiration digital signal 46sqb is input to the input terminal 62q and then output to the second respiration signal extraction unit 63q as the second Q digital signal 58bq.
  • the high-frequency device 4 includes the biological signal including the heartbeat component and the body motion component (hereinafter referred to as “heartbeat component signal”) and the biological signal including the respiratory component (hereinafter referred to as “respiration”). It is divided into the component signal "), and analog signal processing and digital signal processing are performed. Therefore, the dynamic ranges of the heartbeat component signal and the respiratory component signal can be increased at each stage of analog signal processing and digital signal processing. Therefore, each heart rate autocorrelation function processing unit and each respiratory autocorrelation function processing unit can calculate each heart rate and each respiration rate with high accuracy.
  • the heart rate component signal is effective at rest and sleep with little body movement, and the respiratory component signal is when the body motion is non-stationary and when there is still body movement before going to sleep even during sleep To be effective.
  • the high frequency device 4 calculates the heart rate and the respiration rate for display based on the heart rate (first heart rate and second heart rate) and the respiration rate (first respiration rate and second respiration rate). decide. Therefore, the high-frequency device 4 can finally display the heart rate and the respiration rate with higher accuracy on the display 90, and the heart rate with high accuracy even when the living body 10 has relatively large body movement. And the respiratory rate can be calculated.
  • the control blocks (particularly, the radio wave radar units 5a and 5b and the signal processing circuits 400 to 403) of the high frequency devices 1 to 4 may be realized by a logic circuit (hardware) formed in an integrated circuit (IC chip) or the like. Alternatively, it may be realized by software using a CPU (Central Processing Unit).
  • a logic circuit hardware
  • IC chip integrated circuit
  • CPU Central Processing Unit
  • the high-frequency devices 1 to 4 include a CPU that executes instructions of a program that is software that realizes each function, and a ROM (Read Only Memory) in which the program and various data are recorded so as to be readable by the computer (or CPU). ) Or a storage device (these are referred to as “recording media”), a RAM (Random Access Memory) for expanding the program, and the like.
  • recording media a “non-temporary tangible medium” such as a tape, a disk, a card, a semiconductor memory, a programmable logic circuit, or the like can be used.
  • the program may be supplied to the computer via an arbitrary transmission medium (such as a communication network or a broadcast wave) that can transmit the program.
  • a transmission medium such as a communication network or a broadcast wave
  • the present invention can also be realized in the form of a data signal embedded in a carrier wave in which the program is embodied by electronic transmission.
  • the high-frequency device (1, 2, 3, 4) is A biological signal extraction unit (heart rate signal extraction unit 53, respiratory signal extraction unit 63) that extracts a biological signal indicating a specific frequency component of the digital signal from a digital signal representing a biological phenomenon that is a movement appearing in the biological body;
  • the biological signal extracted by the biological signal extraction unit is sampled, an autocorrelation function is calculated based on the predetermined number of the biological signals sampled, and the periodicity of the autocorrelation function is calculated from the peak value of the correlation coefficient.
  • an autocorrelation function processing unit (heart rate autocorrelation function processing unit 54, respiratory autocorrelation function processing unit 63) that calculates biological information representing the biological phenomenon.
  • the high-frequency device extracts only the necessary band with the analog filter and passes only the necessary band with the digital filter. Specifically, the high-frequency device increases the dynamic range in the heartbeat or breathing band by extracting only the necessary band with an analog filter, removing unnecessary signals, amplifying, and then performing analog-digital conversion. Then, the SN ratio of the digital signal is increased in the digital filter. Therefore, the detection accuracy of the peak value of the correlation coefficient detected from the autocorrelation function can be increased.
  • the high-frequency device calculates biological information by detecting a peak value of a correlation coefficient using an autocorrelation function obtained by sampling a biological signal and discriminating the period of the autocorrelation function. Therefore, even when the living body is moving or when there are many noise components such as unnecessary reflection signals, the biological information can be calculated with high accuracy.
  • the high-frequency device is the above aspect 1,
  • the biological signal extraction unit includes a direct reflection signal (reflected signal 12) in which a transmission signal transmitted from the own device toward the living body to extract the biological signal is directly reflected by the living body, and the transmission signal is an object.
  • the biological signal may be extracted from an indirect reflection signal (reflected signal 12) reflected by the living body through reflection at the surface.
  • the biological signal can be extracted from the direct reflection signal and the indirect reflection signal using the biological signal extraction unit. Therefore, biological information can also be calculated for an indirect reflection signal using the autocorrelation function of the biological signal.
  • the high-frequency device can calculate the biological information of the living body existing in the space using the indirect reflection signal.
  • a high-frequency device is the above aspect 1 or 2, wherein
  • the biological signal extraction unit includes a heartbeat signal extraction unit (53) that extracts a heartbeat signal indicating a heartbeat of a specific frequency, and a respiration signal extraction unit (63) that extracts a respiration signal indicating respiration of a specific frequency.
  • the autocorrelation function processing unit includes a heartbeat autocorrelation function processing unit (54) that calculates a heart rate at a predetermined time and a respiratory autocorrelation function processing unit (64) that calculates a respiration rate at a predetermined time. May be.
  • the high-frequency device can determine the predetermined frequency by determining the periodicity of the autocorrelation function in the heartbeat signal and the respiratory signal even when the living body is moving or there are many noise components such as unnecessary reflection signals. Heart rate and respiration rate over time can be calculated. Therefore, it is possible to realize a high-frequency device with improved detection sensitivity regarding the calculation of the heart rate and the respiration rate.
  • the high-frequency device (2) according to aspect 4 of the present invention is the above aspect 3, By inputting analog signals (first analog signal 48a and second analog signal 48b) related to the biological phenomenon, the first biological digital signal (46sa) related to the heart rate component and the body motion component and the second biological digital signal related to the respiratory component.
  • the heartbeat signal extraction unit extracts the heartbeat signal based on the first biological digital signal output from the digital signal generation unit,
  • the respiration signal processing unit may extract the respiration signal based on the second biological digital signal output from the digital signal generation unit.
  • the digital signal generation unit can increase the dynamic range for each of the first biological digital signal and the second biological digital signal. Therefore, the heartbeat signal extraction unit and the respiration signal extraction unit can extract the heartbeat signal and the respiration signal with high sensitivity and accuracy.
  • the high-frequency device (3) according to aspect 5 of the present invention is the above aspect 3, A radio wave radar unit (5b) that outputs an I analog signal (31is) related to the biological phenomenon and a Q analog signal (31qs) that is 90 degrees out of phase with the I analog signal; A heart rate determining unit (55a) for determining the heart rate; A respiratory rate determining unit (55b) for determining the respiratory rate, Based on the I analog signal output from the radio wave radar unit, the digital signal generation unit (first digital signal generation unit 40ci, second digital signal generation unit 40cq) Generating an I digital signal (46si) and a Q digital signal (46sq) relating to a heart rate component, a respiratory component, and a body motion component based on the Q analog signal output from the radio wave radar unit, The heartbeat autocorrelation function processing unit (first heartbeat autocorrelation function processing unit 54i, second heartbeat autocorrelation function processing unit 54q) calculates a first heart rate at a predetermined time based on the I digital signal
  • the high frequency device determines the final heart rate using the first heart rate and the second heart rate. For example, if the error in the numerical values of the first heart rate and the second heart rate is within a few percent, it is considered that the living body is stationary (that is, the living body is in a state suitable for calculating the heart rate). Heart rate can be determined. A similar method can be used for determining the respiratory rate. Therefore, the high frequency device can finally determine the heart rate and the respiratory rate with high accuracy.
  • the high-frequency device determines the first heart rate and the second heart rate based on two digital signals (I digital signal and Q digital signal) that are 90 degrees out of phase.
  • I digital signal which is a cosine wave
  • Q digital signal is a sine wave, so that the signal has a large amplitude, and both signals compensate for the phase change. . Therefore, even when there is a relatively large body movement in the living body, at least one of the heart rates is kept highly accurate. The same applies to the relationship between the first respiratory rate and the second respiratory rate. Therefore, the high-frequency device can calculate a heart rate and a respiratory rate with high accuracy even when there is a relatively large body motion in the living body.
  • the high frequency device (4) includes:
  • the digital signal generation units (I heartbeat digital signal generation unit 49a, I breathing digital signal generation unit 49b, Q heartbeat digital signal generation unit 49c, Q breathing digital signal generation unit 49d) are configured to output the I signal output from the radio wave radar unit.
  • an I heart rate digital signal (46 sia) related to the heart rate component and the body motion component and an I respiration digital signal (46 sib) related to the respiratory component are generated,
  • a Q heart rate digital signal (46 sqa) relating to a heart rate component and a body motion component, and a respiratory component Q digital breathing signal (46 sqb) The heart rate autocorrelation function processing unit calculates the first heart rate based on the I heart rate digital signal, and calculates the second heart rate based on the Q heart rate digital signal
  • the respiration autocorrelation function processing unit may calculate the first respiration rate based on the I respiration digital signal and calculate the second respiration rate based on the Q respiration digital signal.
  • the high-frequency device performs analog signal processing and digital signal processing separately for the heartbeat component signal and the respiratory component signal. Therefore, the dynamic ranges of the heartbeat component signal and the respiratory component signal can be increased at each stage of analog signal processing and digital signal processing. Therefore, each heart rate autocorrelation function processing unit and each respiratory autocorrelation function processing unit can calculate each heart rate and each respiration rate with high accuracy.
  • the high-frequency device determines the heart rate and respiration rate for display based on the first heart rate and the second heart rate, the first respiration rate and the second respiration rate. Therefore, the high-frequency device can finally calculate heart rate and respiration rate with higher accuracy, and can calculate heart rate and respiration rate with high accuracy even when there is relatively large body movement in the living body. be able to.
  • the high-frequency device (1, 2, 3, 4) according to each aspect of the present invention may be realized by a computer.
  • the computer is operated as each unit (limited to software elements) of the high-frequency device.
  • a high-frequency device control program for realizing the above-described high-frequency device by a computer and a computer-readable recording medium on which the control program is recorded fall within the scope of the present invention.
  • the high-frequency device of the present invention is useful for detecting a biological signal in a room, detecting a biological signal during driving of a car, and the like, and can be applied to welfare care and medical care.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Pulmonology (AREA)
  • Artificial Intelligence (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Power Engineering (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

 心拍および呼吸等に関する生体情報を精度よく検出する高周波装置を実現する。高周波装置(1)は、特定の周波数成分を示す生体信号を抽出する生体信号抽出部(心拍信号抽出部(53)、呼吸信号抽出部(63))と、自己相関関数の周期性を判別することにより、生体情報を算出する自己相関関数処理部(心拍自己相関関数処理部(54)、呼吸自己相関関数処理部(64))と、を備える。

Description

高周波装置
 本発明は、心拍数等の生体情報をセンシングするための高周波装置に関する。
 従来から、電磁波を測定対象物に照射し、当該測定対象物で反射する反射波のドップラーシフトを利用して、測定対象物の振動状態や変位を求める手法が広く知られている。また、マイクロ波-ミリ波帯の電磁波は、誘電体等の媒質を透過する性質も有していることから、この電磁波を利用して、人体(被験者)において振動として現れる心臓の拍動や呼吸を、被験者へのマイクロ波照射で検出しようとする試みが近年提案されている。マイクロ波を用いることで、人体に非接触、かつ、被験者が衣服を着たままでの測定が可能となり、センシング時の被験者の負担を軽減することができる。このようなマイクロ波を用いてセンシングする装置の一例としては、特許文献1に開示された生体信号検知装置がある。
 上記の生体信号検知装置について、図13および図14を用いて説明する。生体信号検知装置は、センサ部101、生体信号抽出部102、距離算出部103および生体信号出力判定部104を備える。
 図13はセンサ部101の概略構成を示したものである。図13に示すように、局部発振器301から発信された信号は、分配器302で2分配される。一方の信号は送信アンテナ303へ送信され、もう一方の信号は分配器308でさらに2分配され、混合器306および307へ信号が入力される。分配器308で2分配された2つの信号は、90度位相が異なっている。
 一方、送信アンテナ303から放射された信号は、乗員に向かい、主として乗員の体表面で反射する。そして、体表面で反射した反射信号は受信アンテナ304に入力される。このとき、体表面は、乗員の体動(呼吸、心拍を含む)による振動があるため、乗員への送信信号はドップラーシフトを受ける。そのため、反射信号は、周波数が位相変調した状態で、受信信号として受信アンテナ304に入力される。
 受信アンテナ304に入力された受信信号は、分配器305で2分配され、混合器306および307にそれぞれ入力される。分配器308で入力された信号も、同様に、混合器306および307に入力される。混合器306および307では乗算処理がなされ、ローパスフィルタ309および310を通して、ドップラーシフトを受けたベースバンド信号が出力される。さらに、前記ベースバンド信号は、ADコンバータ311および312によりアナログデジタル変換され、Bi(t)信号およびBq(t)信号として出力される。Bi(t)信号およびBq(t)信号は、瞬時的には90度位相がずれた信号である。
 このBi(t)信号およびBq(t)信号は、図14に示すように、生体信号抽出部102に入力される。生体信号抽出部102では生体信号が抽出され、距離算出部103ではセンサ部101と乗員との推定距離が算出される。そして、生体信号出力判定部104では、前記推定距離に基づいて、生体信号を出力するか否かが、ある閾値を設定して判定される。具体的には、生体信号出力判定部104内の信頼度判定部408において乗員までの距離が一定または閾値以内と判定された場合、生体信号が出力される。
 ここで、デジタル化されたBi(t)信号およびBq(t)信号は、ドップラーシフトを受けた信号、つまりドップラ周波数で位相変調された信号であり、それぞれサイン波、コサイン波の信号で表現できる。生体信号抽出部102内では、これら2つの信号から、位相信号算出部401での演算を経て、位相の時間変化、すなわち、心拍・呼吸それぞれの時間波形を、心拍信号抽出部402および呼吸信号抽出部403で抽出する処理がなされる。
日本国公開特許公報「特開2010-120493号公報(2010年6月3日公開)」
 従来のマイクロ波センサ装置は、上記特許文献1に示されるものが代表的であり、マイクロ波を用いて被験者の微小な動きを検出して心拍等の生体情報を得る仕組みとなっている。具体的には、照射波に対する反射波の位相変化(位相差)を検出することで、被験者の心拍などの体表の微動を検出しようとするものである。
 ここで、直接Bi(t)信号およびBq(t)信号から位相の時間変化を検出しようとすると、体動(心拍および呼吸を除く)が優勢である(すなわち、体が動いている)場合には、心拍信号および呼吸信号が微弱な信号であることから、これらの信号が体動に埋もれやすかった。また、無線環境下で上記検出を行うがゆえに、不要反射信号等のノイズ成分(クラッター)に埋もれやすく、直接波(被験者に直接照射され、且つ、反射波が直接受信されるもの)のみでしか、心拍・呼吸情報を抽出できないという課題もあった。
 例えば、被験者が椅子に座っている状態で、電波が直接胸(または背中)に照射される場合、心拍・呼吸情報は取得できるものの、被験者が寝ている状況で、被験者の寝姿(うつ伏せ、仰向け、横向き、寝ている間で変化する場合等)に対して、常時、心拍・呼吸成分を測定・モニターすることができなかった。
 加えて、呼吸信号に対して、心拍信号はさらに微弱な信号であるために、特に心拍成分を精度よく抽出することが困難であった。
 本発明は、上記の問題点を解決するためになされたものであり、その目的は、生体の心拍および呼吸等を表す生体情報を、無線で精度よく検出することができる高周波装置を提供することにある。
 上記の課題を解決するために、本発明の一態様に係る高周波装置は、
 生体に現れる動きである生体現象を表すデジタル信号から、上記デジタル信号の特定の周波数成分を示す生体信号を抽出する生体信号抽出部と、
 上記生体信号抽出部により抽出された上記生体信号をサンプリングし、上記サンプリングされた所定数の上記生体信号に基づいて自己相関関数を算出し、相関係数のピーク値から上記自己相関関数の周期性を判別することにより、上記生体現象を表す生体情報を算出する自己相関関数処理部と、を備える。
 本発明の一態様によれば、生体が動いている場合または不要反射信号等のノイズ成分が多い場合においても、自己相関関数を用いて、心拍および呼吸等を表す生体情報を精度よく検出する高周波装置を実現することができる。
(a)は、本発明の実施形態1~4に係る高周波装置の使用方法の一例を示す側面図である。(b)は、本発明の実施形態1~4に係る高周波装置の使用方法の一例を示す上面図である。 本発明の実施形態1および2に係る高周波装置が備える電波レーダ部の概略構成を示すブロック図である。 本発明の実施形態1に係る高周波装置が備える信号処理回路の概略構成を示すブロック図である。 上記高周波装置が備えるデジタル信号処理部および表示器の概略構成を示すブロック図である。 本発明の実施形態2に係る高周波装置が備える信号処理回路の概略構成を示すブロック図である。 上記高周波装置が備えるデジタル信号処理部および表示器の概略構成を示すブロック図である。 本発明の実施形態3および4に係る高周波装置が備える電波レーダ部の概略構成を示すブロック図である。 本発明の実施形態3に係る高周波装置が備える信号処理回路の概略構成を示すブロック図である。 上記高周波装置が備えるデジタル信号処理部および表示器の概略構成を示すブロック図である。 (a)は、上記高周波装置が備える、第1心拍自己相関関数処理部および第2心拍自己相関関数処理部においてサンプリングされた心拍信号に関して、当該心拍信号の振幅値とサンプル数との関係を示すグラフである。(b)は、上記心拍信号に関して、自己相関関数の相関係数とサンプル数との関係を示すグラフである。 本発明の実施形態4に係る高周波装置が備える信号処理回路の概略構成を示すブロック図である。 上記高周波装置が備えるデジタル信号処理部および表示器の概略構成を示すブロック図である。 従来のマイクロ波センサ装置が備えるセンサ部の概略構成を示したブロック図である。 従来のマイクロ波センサ装置の概略構成を示したブロック図である。
 〔実施形態1〕
 本発明に係る一実施形態について、図1~4に基づいて説明すれば、以下のとおりである。
 <高周波装置1の使用方法>
 以下に、図1を用いて、高周波装置1の使用方法について、ベッド91の手すり91bに高周波装置1を取り付けた場合を例に挙げて説明する。図1の(a)は、本実施形態に係る高周波装置1の使用方法の一例を示す側面図である。図1の(b)は、本実施形態に係る高周波装置1の使用方法の一例を示す上面図である。
 高周波装置1は、送信信号11が生体10で反射された反射信号12(図2参照)に基づいて、生体10の生体現象を表す生体情報を得るための装置である。ここで、生体現象とは、主として生体10の上半身に現れる動きのことであり、心拍成分、呼吸成分および体動成分等を指す。また、体動とは、生体10の上半身に現われる、心拍および呼吸以外の動きをいう。さらに、生体情報としては、例えば、所定時間の心拍数および呼吸数、心拍波形、呼吸波形、体動波形等がある。高周波装置1は、電波レーダ部5aおよび信号処理回路400を備える(図2参照)。
 図1に示すように、高周波装置1から放射される送信信号11(図2参照)は、望ましくは、生体10の上半身に放射されるように、電波レーダ部5aに備えられた送信アンテナ25(図2参照)の指向性パターン93が調整される。前記調整の一例としては、マイクロストリップパッチアンテナのパッチ数により調整することができる。反射信号12は高周波装置1の受信信号となる。上記指向性パターン93の調整によって、生体10が睡眠中であっても、心拍成分、呼吸成分および体動成分をモニターすることができる。
 とりわけ、高周波装置1は、ベッド91の手すり91bに取り付けられることによって、生体10の上半身に向けて、送信信号11を、ベッド91のマットレス91aにおける生体10との接触面に対して略平行方向に放射する。そのため、生体10の寝姿(うつ伏せ、仰向け、横向き、寝ている間で体が動くなど)に対して、常に、心拍成分、呼吸成分および体動成分を測定・モニターすることができる。
 これは、生体10に比較的近い位置から送信信号11を放射することに加え、送信信号11を上記略水平方向に照射することにより、送信信号11および反射信号12が照射される反射物の存在が少なくなるからである。すなわち、高周波装置1が受信する不要反射信号が少なくなるからである。また、高周波装置1は、後述する自己相関関数を用いて心拍成分を表す心拍信号および呼吸成分を表す呼吸信号を処理している。そのため、高周波装置1は、不要反射信号が多数存在する寝室においても、後述する間接反射信号から心拍成分等を計測・モニターすることができるためである。
 なお、高周波装置1の取り付け位置については、本実施形態のようなベッド91の手すり91bに限定されない。例えば、ベッド91の床下部(望ましくは、ベッド91のマットレス91aにおける、生体10との接触面と対向する部分)に取り付けてもよい。高周波装置1をこのように取り付けても、マットレス91aに対する送信信号11の反射成分は少ないため、本実施形態のように高周波装置1が取り付けられた場合と、ほぼ同様な効果が得られる。
 また、本実施形態における高周波装置1の使用方法は、後述する高周波装置2~4についても、同様に採用される。
 <電波レーダ部5aによる信号処理>
 以下に、図2を用いて、電波レーダ部5aによる信号処理を説明する。図2は、電波レーダ部5aの概略構成を示すブロック図である。電波レーダ部5aは、生体10に向けて送信信号11を放射するとともに、生体10の体表で反射された反射信号12を受信し、当該反射信号12に基づいてアナログ信号33sを生成する。図2に示すように、電波レーダ部5aは、発振器21、増幅器22、送信アンテナ25、受信アンテナ30、低雑音アンプ31、ミキサ32およびフィルタ33を備えている。
 図2に示すように、発振器21から出力されるマイクロ波正弦波(図示しない)は、増幅器22により増幅された後、送信アンテナ25から送信信号11として放射される。放射された送信信号11は、生体10の体表、例えば、生体10の胸部表面に当たって反射される。このとき、生体10には、心拍、呼吸および体動による胸部の動きがあり、この動きによって胸部に当たった送信信号11にドップラーシフトが生じる。すなわち、送信信号11の周波数に、心拍および呼吸を含む胸の動きの周波数が付加されることで、送信信号11が周波数・位相変調され、反射信号12となる。そして、反射信号12は、受信アンテナ30に入力される。
 なお、反射信号12には、生体10に直接照射され、且つ、反射後、受信アンテナ30によって直接受信される直接反射信号と、生体10以外の物体での反射を介して受信アンテナ30に受信される間接反射信号との、両方の信号が含まれる。
 受信アンテナ30に入力された反射信号12は、低雑音アンプ31によって増幅された後、アナログ信号31sとしてミキサ32に入力される。ミキサ32には、発振器21から発信され、増幅器22で増幅された後、2分配されたうちの一方のアナログ信号22sも入力される。アナログ信号31sは、アナログ信号22sによって周波数がダウンコンバージョンされた状態でフィルタ33に入力される。
 フィルタ33から出力されるアナログ信号33sは、ベースバンド信号となり、胸の動きによってドップラーシフトを受けた信号として出力される。すなわち、アナログ信号33sは、生体10の心拍成分、呼吸成分および体動成分を含んでいる。アナログ信号33sは、後述する信号処理回路400に入力され、信号処理された後、心拍情報、呼吸情報および体動情報として、後述する表示器90に表示される。
 <信号処理回路400による信号処理>
 以下に、図3および図4を用いて、信号処理回路400による信号処理について説明する。図3は、信号処理回路400の概略構成を示すブロック図である。図4は、デジタル信号処理部(以下、「DSP」とする)500および後述する表示器90の概略構成を示すブロック図である。信号処理回路400は、入力されたアナログ信号33sに基づいて、生体10の心拍情報(1分間当たりの心拍数および心拍波形)、呼吸情報(1分間当たりの呼吸数および呼吸波形)並びに体動情報(体動波形)を算出・決定する。図3に示すように、信号処理回路400は、デジタル信号生成部40aおよびDSP500を備えている。
 デジタル信号生成部40aは、入力されたアナログ信号33sに基づいて、生体10の心拍成分、呼吸成分および体動成分を含むデジタル信号46sを生成する。図3に示すように、デジタル信号生成部40aは、入力端子41、ハイパスフィルタ43、ローパスフィルタ44、増幅器45およびAD変換器46を備えている。
 DSP500は、入力されたデジタル信号46sに基づいて、上記心拍情報、呼吸情報および体動情報を算出・決定する。図4に示すように、DSP500は、入力端子51、心拍信号抽出部(生体信号抽出部)53、心拍自己相関関数処理部(自己相関関数処理部)54、心拍数決定部55、呼吸信号抽出部(生体信号抽出部)63、呼吸自己相関関数処理部(自己相関関数処理部)64および呼吸数決定部65を備えている。
 (デジタル信号生成部40aによる信号処理)
 図3に示すように、アナログ信号33s(図2参照)は、入力端子41に入力された後、ハイパスフィルタ43(dc成分の抑圧を目的とする)およびローパスフィルタ44によって帯域制限される。一例として、ハイパスフィルタ43での通過帯域のカットオフは0.1Hz、ローパスフィルタ44での通過帯域のカットオフは10Hzであり、心拍、呼吸および体動にそれぞれ関連した周波数帯域のアナログ信号を通過させた後、増幅器45で増幅する。
 アナログ信号33sをハイパスフィルタ43およびローパスフィルタ44に通過させることにより、例えば、扇風機等の高速動作に基づくドップラーシフトを受けたアナログ信号33sから、心拍、呼吸および体動以外のアナログ信号を抑圧することができる。また、増幅器45によって、ハイパスフィルタ43およびローパスフィルタ44を通過した高速なアナログ信号が飽和動作し、ひいては非線形動作となることを防いでいる。加えて、ローパスフィルタ44は、AD変換器46のアンチエリアスフィルタとしても機能している。
 増幅器45から出力されたアナログ信号(図示しない)は、AD変換器46に入力されて、アナログデジタル変換される。一例として、12ビット50Hzのサンプリングレートでアナログデジタル変換する場合、ローパスフィルタ44での通過帯域のカットオフは、25Hz以下となる(本実施形態では10Hzとしている)。
 AD変換器46によるアナログデジタル変換によって生成されたデジタル信号46sは、デジタル信号処理部(以下、「DSP」とする)500に入力される。そして、DSP500によって上記心拍情報、呼吸情報および体動情報が算出された後、これらの情報は表示器90に表示される。
 (DSP500による信号処理)
 図4に示すように、入力端子51に入力されたデジタル信号46s(図3参照)は、第1のデジタル信号58a、第2のデジタル信号58bおよび第3のデジタル信号58cに分配され、これら3つのデジタル信号の処理が並列してなされる。
 まず、第1のデジタル信号58aは、心拍信号抽出部53に入力される。そして、入力された第1のデジタル信号58aが、心拍信号抽出部53が有するハイパスフィルタ53aとローパスフィルタ53bとで帯域制限されることで、特定の周波数の心拍を示す心拍信号が抽出される。
 本実施形態では、ハイパスフィルタ53aは、通過帯域のカットオフ周波数が0.8Hzのものを、ローパスフィルタ53bは、通過帯域のカットオフ周波数が3.0Hzのものを用いた。これは、呼吸成分が、平静時0.2Hz~0.33Hz(1分間に12回~20回程度)であることから、極力、呼吸成分を抑圧するために、ハイパスフィルタ53aの通過帯域のカットオフ周波数として0.8Hzのものを選んでいる。一方、ローパスフィルタ53bは、心拍成分が、平静時50~90程度(周波数としては、0.8Hz~1.5Hz)であるため、通過帯域のカットオフ周波数として3.0Hzものを選んでいる。
 心拍信号抽出部53で抽出された心拍信号は2分配され、その一方の信号は、心拍自己相関関数処理部54に入力される。そして、心拍自己相関関数処理部54において自己相関関数を用いた演算処理がなされることにより、1分間の心拍数(生体現象の値;以下、単に「心拍数」とする)が算出される。心拍自己相関関数処理部54は、サンプリング処理部54a、心拍自己相関関数演算部54bおよびピーク検出部54cを有している。
 本実施形態では、AD変換器46が12ビット50Hzのサンプリングレートでアナログデジタル変換するため、サンプリング処理部54aは、第1のデジタル信号58aから、1サンプル20m秒で256回、心拍信号をサンプリングする。心拍自己相関関数演算部54bは、サンプリングされた心拍信号が入力されることにより、自己相関関数を算出する。そして、ピーク検出部54cにおいて、自己相関関数から相関係数のピーク値が検出され、当該自己相関関数の周期が算出される。換言すれば、ピーク検出部54cによって自己相関関数の周期性が判別される。ピーク検出部54cは、算出された周期を心拍数に換算する。自己相関関数においては、0秒後のデータが最も相関係数が大きいため、0秒のデータを除いて、自己相関係数のピーク値が検出される。
 本実施形態では、サンプリング処理部54aによるサンプリング処理で5秒(1サンプル20m秒で256回のサンプリング)、心拍自己相関関数演算部54bによる自己相関関数の算出およびピーク検出部54cによる周期の算出で約5秒かかる。そのため、約10秒で一つの心拍数が算出される。
 なお、自己相関関数を用いた心拍成分に関する周期性の判別方法の詳細については、実施形態3にて説明する。
 心拍自己相関関数処理部54で算出された心拍数は、当該心拍数を示す信号として(図示しない)心拍数決定部55に入力される。心拍数決定部55は、前記信号に基づいて、表示器90(具体的には、後述する心拍数表示部90a)に表示するための心拍数を決定し、当該決定した心拍数を示す心拍数信号55sを出力する。心拍数信号55sは表示器90への入力信号となる。
 心拍数決定部55における心拍数の決定は、一つの方法に限定されない。本実施形態では、心拍自己相関関数処理部54から約10秒ごとに一つの信号が出力されるため、例えば、3つの信号毎に、平均値、移動平均値または中央値等を得るための演算がなされてもよい。換言すれば、表示器90に表示される心拍数として、安定した数値が出力される方法であればよい。
 一方、心拍信号抽出部53で抽出された心拍信号のうち、2分配された他方は、時間の経過に伴う心拍振幅の変化を表す心拍波形信号71aとなり、表示器90への入力信号となる。また、第2のデジタル信号58bは、そのまま、時間の経過に伴う体動振幅の変化を表す体動波形信号71bとなり、表示器90への入力信号となる。
 次に、第3のデジタル信号58cは、呼吸信号抽出部63に入力される。そして、入力された第3のデジタル信号58cが、呼吸信号抽出部63が有するローパスフィルタ63aで帯域制限されることで、特定の周波数の呼吸を示す呼吸信号が抽出される。
 本実施形態では、ローパスフィルタ63aは、通過帯域のカットオフ周波数が0.75Hzのものを用いた。これは、呼吸成分が、平静時0.2Hz~0.33Hz(1分間に12回~20回程度)であることから、呼吸信号抽出部63において、極力、心拍成分を抑圧するためである。
 呼吸信号抽出部63で抽出された呼吸信号は2分配され、その一方の信号は、呼吸自己相関関数処理部64に入力される。そして、呼吸自己相関関数処理部64において自己相関関数を用いた演算処理がなされることにより、1分間の呼吸数(生体現象の値;以下、単に「呼吸数」とする)が算出される。呼吸自己相関関数処理部64は、サンプリング処理部64a、呼吸自己相関関数演算部64bおよびピーク検出部64cを有している。
 本実施形態では、AD変換器46が12ビット50Hzのサンプリングレートでアナログデジタル変換するため、サンプリング処理部64aは、第3のデジタル信号58cから、1サンプル20m秒で512回、呼吸信号をサンプリングする。呼吸自己相関関数演算部64bは、サンプリングされた呼吸信号が入力されることにより、自己相関関数を算出する。そして、ピーク検出部64cにおいて、自己相関関数から相関係数のピーク値が検出され、当該自己相関関数の周期が算出される。換言すれば、ピーク検出部64cによって自己相関関数の周期性が判別される。ピーク検出部64cは、算出された周期を呼吸数に換算する。
 本実施形態では、サンプリング処理部64aによるサンプリング処理で10秒(1サンプル20m秒で512回のサンプリング)、呼吸自己相関関数演算部64bによる自己相関関数の算出およびピーク検出部64cによる周期の算出で約5秒かかる。そのため、約15秒で一つの呼吸数が算出される。
 なお、自己相関関数を用いた呼吸成分に関する周期性の判別方法の詳細については、実施形態3にて説明する。
 呼吸自己相関関数処理部64で算出された呼吸数は、当該呼吸数を示す信号として(図示しない)呼吸数決定部65に入力される。呼吸数決定部65は、前記信号に基づいて、表示器90(具体的には、後述する呼吸数表示部90e)に表示するための呼吸数を決定し、当該決定した呼吸数を示す呼吸数信号65sを出力する。呼吸数信号65sは表示器90への入力信号となる。
 呼吸数決定部65における呼吸数の決定は、一つの方法に限定されない。本実施形態では、呼吸自己相関関数処理部64から約15秒ごとに一つの信号が出力されるため、例えば、3つの信号毎に、平均値、移動平均値または中央値等を得るための演算がなされてもよい。換言すれば、表示器90に表示される呼吸数として、安定した数値が出力される方法であればよい。
 一方、呼吸信号抽出部63で抽出された呼吸信号のうち、2分配された他方は、時間の経過に伴う呼吸振幅の変化を表す呼吸波形信号71cとなり、表示器90への入力信号となる。
 <表示器90による心拍・呼吸・体動情報の表示>
 以下に、図4を用いて、表示器90による生体10の心拍情報、呼吸情報および体動情報の表示について説明する。図4に示すように、表示器90は、心拍数表示部90a、心拍波形表示部90b、体動波形表示部90c、呼吸波形表示部90dおよび呼吸数表示部90eを備えている。
 DSP500から出力されたそれぞれの信号は、表示器90に入力される。具体的には、心拍数信号55sは心拍数表示部90aに、心拍波形信号71aは心拍波形表示部90bに、体動波形信号71bは体動波形表示部90cに、呼吸波形信号71cは呼吸波形表示部90dに、呼吸数信号65sは呼吸数表示部90eに、それぞれ入力される。そして、心拍数表示部90aには心拍数(数値)が、心拍波形表示部90bには心拍波形(グラフ)が、体動波形表示部90cには体動波形(グラフ)が、呼吸波形表示部90dには呼吸波形(グラフ)が、呼吸数表示部90eには呼吸数(数値)が、それぞれ表示される。
 表示器90の例としては、パソコン、携帯電話またはスマートフォン等の表示機能・データ保存機能を有するものが挙げられる。なお、DSP500と表示器90とが、一体となっていても構わない。例えば、デジタル信号生成部40aにおいてアナログデジタル変換された信号を無線で送信し、DSP500と表示器90とを備えたパソコン等が当該信号を受信することで、当該パソコン等による信号処理および表示がなされてもよい。
 <効果>
 以上のように、本実施形態によれば、高周波装置1は、アナログフィルタ(ハイパスフィルタ43、ローパスフィルタ44)で必要な帯域のみ抽出し、かつデジタルフィルタ(ハイパスフィルタ53a、ローパスフィルタ53bおよび63a)でも必要な帯域のみを通過させる。具体的には、高周波装置1は、アナログフィルタで必要な帯域のみ抽出して不要信号を除き、且つ増幅した後、アナログデジタル変換することによって、心拍または呼吸帯域でのダイナミックレンジを大きくする。そして、デジタルフィルタにおいてデジタル信号のSN比を大きくする。そのため、自己相関関数から検出される相関係数のピーク値の検知精度を高くすることができる。
 また、高周波装置1は、心拍信号および呼吸信号をサンプリングして得られた自己相関関数を用いて相関係数のピーク値の検出し、当該自己相関関数の周期を判別することで、心拍数および呼吸数を算出する。そのため、生体10が動いている場合または不要反射信号等のノイズ成分が多い場合においても、心拍数および呼吸数を精度高く算出することができる。
 また、一般に、間接反射信号については、波形が変形したり、ノイズが付加したりするため、当該間接反射信号から生体情報を算出することは困難とされている。
 その点、本実施形態によれば、直接反射信号および間接反射信号から、心拍信号抽出部53および呼吸信号抽出部63を用いて、心拍信号および呼吸信号を抽出することができる。そのため、間接反射信号についても、心拍信号および呼吸信号の自己相関関数を用いて、心拍数および呼吸数等を算出することができる。
 また、寝室等の比較的狭い空間に生体が存在する場合には、当該空間内に間接反射信号が多数存在する。そのため、高周波装置1は、直接反射信号を捕捉できない場合においても、上記間接反射信号を用いて当該空間内に存在する生体10の心拍数および呼吸数等を算出することができる。
 さらに、上記空間に生体10が存在する場合においては、例えば、生体10がうつ伏せで横たわっている等、生体10の所望の部位に直接送信信号11を照射することが困難な状況であっても、壁等で送信信号11を反射させることによって、当該部位に間接的に送信信号11を照射することができる。そのため、間接反射信号を用いることで、生体10の姿勢に拘らず心拍数および呼吸数等を算出することができる。
 〔実施形態2〕
 本発明の他の実施形態について、図5および図6に基づいて説明すれば、以下のとおりである。なお、説明の便宜上、前記実施形態にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を省略する。
 本実施形態に係る高周波装置2は、信号処理回路400に代えて、信号処理回路401を備えている点で、実施形態1に係る高周波装置1と異なる。
 <信号処理回路401による信号処理>
 以下に、図5および図6を用いて、信号処理回路401による信号処理について説明する。図5は、信号処理回路401の概略構成を示すブロック図である。図6は、DSP501および表示器90の概略構成を示すブロック図である。図5に示すように、信号処理回路401は、デジタル信号生成部40bおよびDSP501を備えている。
 デジタル信号生成部40bは、入力されたアナログ信号33s(図2参照)に基づいて、生体10の心拍成分および体動成分を含むデジタル信号46sa、並びに生体10の呼吸成分を含むデジタル信号46sbを生成する。
 (デジタル信号生成部40bによる信号処理)
 図5に示すように、デジタル信号生成部40bは、入力端子41、ハイパスフィルタ43aおよび43b、ローパスフィルタ44aおよび44b、増幅器45aおよび45b並びにAD変換器46aおよび46bを備えている。
 図5に示すように、アナログ信号33sは、入力端子41に入力された後、端子48で第1アナログ信号48aと第2アナログ信号48bに2分配される。
 第1アナログ信号48aは、この後、心拍および体動用の帯域に制限され増幅される。一例として、ハイパスフィルタ43aでの通過帯域のカットオフは0.75Hz、ローパスフィルタ44aでの通過帯域のカットオフは10Hzであり、心拍および体動に関連した周波数帯のアナログ信号を通過させたあと、増幅器45aで増幅する。ハイパスフィルタ43aによる帯域制限を上記のようにしたのは、実施形態1と同様、極力、呼吸成分を抑圧するためである。
 増幅器45aから出力されたアナログ信号(図示しない)は、AD変換器46aに入力されて、アナログデジタル変換される。AD変換器46aによるアナログデジタル変換によって生成された、心拍成分および体動成分に関する第1生体デジタル信号46saは、DSP501に入力される。
 成人の場合、通常、呼吸時の胸部表面の動きは4mm~12mmであり、心拍時の胸部表面の動きは0.5mm程度であるため、心拍に関するアナログ信号は、呼吸に関するアナログ信号に比べて微小である。その点、ハイパスフィルタ43aおよびローパスフィルタ44aによる帯域制限を行うことにより、呼吸時の胸部表面の動きによって増幅器45aが飽和するのを防ぎ、且つ、AD変換器46aにおいて心拍成分のみをアナログデジタル変換することができる。そのため、時間経過に対する心拍波形の振幅のダイナミックレンジを大きくし、DSP501によるデジタル信号処理の感度および精度を向上させることができる。一例として、振幅の小さい心拍成分および体動成分のアナログ信号のみを、16ビットでアナログデジタル変換することができ、特に、体動成分が小さい安静時や睡眠時に効果を発揮する。
 一方、第2アナログ信号48bは、呼吸用の帯域に制限され増幅される。一例として、ハイパスフィルタ43bでの通過帯域のカットオフは0.1Hz、ローパスフィルタ44bでの通過帯域のカットオフは0.75Hzとし、呼吸成分に関連した周波数帯域の信号を通過させたあと、増幅器45bで増幅する。
 増幅器45bから出力されたアナログ信号(図示しない)は、AD変換器46bに入力されて、アナログデジタル変換される。AD変換器46bによるアナログデジタル変換によって生成された、呼吸成分に関する第2生体デジタル信号46sbは、DSP501に入力される。
 ハイパスフィルタ43bではdc成分を抑圧し、ローパスフィルタ44bでは心拍成分および体動成分を抑圧する。そして、増幅器45bおよびAD変換器46bにおいて、第2アナログ信号48bから得られた呼吸成分信号のダイナミックレンジを大きくすることで、DSP501によるデジタル信号処理の感度および精度を向上させることができる。一例として、呼吸成分のアナログ信号のみを、16ビットでアナログデジタル変換することができ、体動成分の振幅が大きいとき、または睡眠時でも眠りに入る前の体動があるときに効果を発揮する。
 (DSP501による信号処理)
 図6に示すように、第1生体デジタル信号46sa(図5参照)は入力端子51に、第2生体デジタル信号46sb(図5参照)は入力端子61に、それぞれ入力される。
 入力端子51に入力された第1生体デジタル信号46saは、心拍成分に関する心拍デジタル信号59aおよび体動成分に関する体動デジタル信号59bに分配され、これらの2つの信号が並列して処理される。このうち、心拍デジタル信号59aに基づいて心拍数信号55sおよび心拍波形信号71aを得る信号処理については、実施形態1と同様であるため、その説明を省略する。
 また、入力端子62に入力された第2生体デジタル信号46sbは、呼吸デジタル信号59cとしてそのまま信号処理される。呼吸デジタル信号59cに基づいて呼吸数信号65sおよび呼吸波形信号72c2を得る信号処理については、実施形態1と同様であるため、その説明を省略する。
 一方、体動デジタル信号59bについては、体動成分として、呼吸成分に関連した生体10の胸部の動きを付加するため、呼吸信号抽出部63によって抽出された呼吸信号の分岐信号71c1が加算される。そして、ゆっくりとした体動から速い体動まで幅広く表現する体動波形信号71bとして、DSP501から出力される。
 <効果>
 以上のように、本実施形態によれば、デジタル信号生成部40bは、第1生体デジタル信号46saおよび第2生体デジタル信号46sbのそれぞれについてダイナミックレンジを大きくすることができる。そのため、心拍信号抽出部53および呼吸信号抽出部63は、高い感度・精度で心拍信号および呼吸信号を抽出することができる。
 具体的には、例えば、体動成分の振幅が小さい生体の睡眠時等(すなわち、呼吸成分の影響が大きい場合)でも、デジタル信号生成部40bにおいて第1アナログ信号48aのうちの呼吸成分を抑圧することで、ダイナミックレンジを大きくした第1生体デジタル信号46saを生成することができる。また、体動成分の振幅が大きい生体の動作時(すなわち、体動成分の影響が大きい場合)でも、デジタル信号生成部40bにおいて第2アナログ信号48bのうちの体動成分を抑圧することで、ダイナミックレンジを大きくした第2生体デジタル信号46sbを生成することができる。
 <変形例>
 なお、本実施形態では、デジタル信号生成部40bにおいて、心拍成分および体動成分に関する第1生体デジタル信号46saと、呼吸成分に関する第2生体デジタル信号46sbとを生成している。しかし、デジタル信号生成部40bにおける信号生成は前記場合に限定されない。例えば、デジタル信号生成部40bは、心拍成分およびこれと同等の速さの体動成分に関するデジタル信号、心拍成分より高速な動きの体動成分に関するデジタル信号、呼吸成分信号の、3つのデジタル信号を生成してもよい。また、4つ以上のデジタル信号を生成してもよい。
 また、実施形態1と同様に、DSP501と表示器90とが、一体となっていても構わない。例えば、デジタル信号生成部40bにおいてアナログデジタル変換された信号を無線で送信し、DSP501と表示器90とを備えたパソコン等が当該信号を受信することで、当該パソコン等による信号処理および表示がなされてもよい。
 〔実施形態3〕
 本発明の他の実施形態について、図7~図10に基づいて説明すれば、以下のとおりである。なお、説明の便宜上、前記実施形態にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を省略する。
 本実施形態に係る高周波装置3は、電波レーダ部5aに代えて電波レーダ部5bを、信号処理回路400および401に代えて、信号処理回路402を備えている点で、実施形態1および2に係る高周波装置1および2と異なる。
 <電波レーダ部5bによる信号処理>
 以下に、図7を用いて、電波レーダ部5bによる信号処理を説明する。図7は、電波レーダ部5bの概略構成を示すブロック図である。図7に示すように、電波レーダ部5bは、発振器21、増幅器22、送信アンテナ25、受信アンテナ30、低雑音アンプ31、Iミキサ32i、Qミキサ32q、フィルタ33i、フィルタ33qおよび位相器38を備えている。
 図7に示すように、送信アンテナ25から送信信号11が放射され、反射信号12が受信アンテナ30に入力されるところまでの信号処理は、実施形態1および2と同様である。反射信号12が、直接反射信号および間接反射信号の両方を含む点も、実施形態1および2と同様である。
 受信アンテナ25に入力された信号は、低雑音アンプ31によって増幅された後、Iアナログ信号31isとQアナログ信号31iqとに2分配される。Iアナログ信号31isは、位相器38に入力されることで、Qアナログ信号31iqに対して位相が90度ずれる。そして、位相器38から出力されたIアナログ信号31isはIミキサ32iに入力され、Qアナログ信号31iqはQミキサ32qに入力される。
 一方、増幅器22で増幅された後、2分配された方の一方のアナログ信号22sも、さらに2分配されてI側の局部発振信号22isおよびQ側の局部発振信号22qsとなり、Iミキサ32iおよびQミキサ32qにそれぞれ入力される。なお、本実施形態では、Iアナログ信号31isとQアナログ信号31iqとの位相を90度ずらしているが、I側の局部発振信号22isとQ側の局部発振信号22qsとの位相を90度ずらしても構わない。
 位相器38から出力されたIアナログ信号31isはIミキサ32iによって、Qアナログ信号31iqはQミキサ32qによって、それぞれ周波数ダウンコンバージョンされた上で、フィルタ33iおよびフィルタ33qに入力される。そして、フィルタ33iからIベースバンド信号33isが、フィルタ33qからQベースバンド信号33qsが、それぞれ出力される。なお、Iベースバンド信号33isはコサイン波となっており、Qベースバンド信号33qsはサイン波となっている。
 Iベースバンド信号33isおよびQベースバンド信号33qsは、ともに、心拍および呼吸を含む胸の動きによって、ドップラーシフトを受けたアナログ信号である。Iベースバンド信号33isとQベースバンド信号33qsとの差異は、位相器38を通過した信号か否かであるが、受信アンテナ25に入力される信号速度が時間とともに変化するため、瞬時的には、両信号は90度位相がずれる。そのため、信号速度の大きさや方向に応じて、電波レーダ部5bから出力される積分信号としてのIベースバンド信号33isとQベースバンド信号33qsとの位相関係は、常に時間変動している。
 このように、I系列用およびQ系列用の2つのチャンネルから、それぞれ自己相関関数を算出するようにしたのは、生体10に微妙な体動がある場合、Iベースバンド信号33isおよびQベースバンド信号33qsの時間に対する振幅値は変わりやすくなる。そのため、例えば、Iベースバンド信号33isに比べてQベースバンド信号33qsが微弱な信号となる可能性がある。したがって、片方のチャンネルの自己相関関数によるピーク値の検出だけでは、当該ピーク値の検知感度および検知精度が悪くなってしまうからである。
 <信号処理回路402による信号処理>
 以下に、図8および図9を用いて、信号処理回路402による信号処理について説明する。図8は、信号処理回路402の概略構成を示すブロック図である。図9は、DSP502および表示器90の概略構成を示すブロック図である。図8に示すように、信号処理回路402は、第1デジタル信号生成部(デジタル信号生成部)40ci、第2デジタル信号生成部(デジタル信号生成部)40cqおよびDSP502を備えている。
 (第1デジタル信号生成部40ciおよび第2デジタル信号生成部40cqによる信号処理)
 図8に示すように、Iベースバンド信号33is(図7参照)は入力端子41iに入力され、Qベースバンド信号33qs(図7参照)は入力端子41qに入力される。そして、入力されたIベースバンド信号33isはI系列の第1デジタル信号生成部40ciで信号処理され、入力されたQベースバンド信号33qsはQ系列の第2デジタル信号生成部40cqで処理される。第1および第2デジタル信号生成部40ciおよび40cqは、デジタル信号生成部40a(図3参照)と同一の構成ものが、I系列用およびQ系列用として備わっているだけなので、第1および第2デジタル信号生成部40ciおよび40cqでの信号処理についての説明は省略する。
 第1デジタル信号生成部40ciから出力されたIデジタル信号46si、および、第2デジタル信号生成部40cqから出力されたQデジタル信号46sqは、それぞれ、体動成分(心拍成分および呼吸成分を含む)の周波数帯に制限され、増幅されたデジタル信号として出力された後、DSP502で信号処理される。なお、Iデジタル信号46siはコサイン波となっており、Qデジタル信号46sqはサイン波となっている。
 (DSP502による信号処理)
 図9に示すように、DSP502は、入力端子52iおよび52q、第1心拍信号抽出部(心拍信号抽出部)53i、第2心拍信号抽出部(心拍信号抽出部)53q、第1呼吸信号抽出部(呼吸信号抽出部)63i、第2呼吸信号抽出部(呼吸信号抽出部)63q、第1心拍自己相関関数処理部(心拍自己相関関数処理部)54i、第2心拍自己相関関数処理部(心拍自己相関関数処理部)54q、第1呼吸自己相関関数処理部(呼吸自己相関関数処理部)64i、第2呼吸自己相関関数処理部(呼吸自己相関関数処理部)64q、心拍数決定部55aおよび呼吸数決定部65aを備えている。
 図9に示すように、入力端子52iに入力されたIデジタル信号46si(図8参照)は、第1のIデジタル信号58ai、第2のIデジタル信号58biおよび第3のIデジタル信号58ciに3分配される。
 第1のIデジタル信号58aiおよび第3のIデジタル信号58ciに基づいて、自己相関関数の周期を算出するまでの信号処理の流れについては、実施形態1および2と同様であるため、その説明を省略する。また、第2のIデジタル信号58biは、そのまま、時間の経過に伴う体動振幅の変化を表す体動波形信号71bとして出力され、表示器90への入力信号となる。
 一方、入力端子52qに入力されたQデジタル信号46sq(図8参照)は、第1のQデジタル信号58aqおよび第2のQデジタル信号58bqに2分配される。
 第1のQデジタル信号58aqおよび第2のQデジタル信号58bqに基づいて、自己相関関数の周期を算出するまでの信号処理の流れについては、実施形態1および2と同様であるため、その説明を省略する。
 次に、第1心拍自己相関関数処理部54iで算出された心拍数である第1心拍数、および、第2心拍自己相関関数処理部54qで算出された心拍数である第2心拍数は、それぞれ当該心拍数を示す信号として(図示しない)心拍数決定部55aに入力される。具体的には、第1心拍数は、心拍数決定部55aに備えられた第1心拍数決定部55iに入力される。また、第2心拍数は、同じく心拍数決定部55aに備えられた第2心拍数決定部55qに入力される。なお、第1心拍数および第2心拍数の決定方法については、実施形態1で説明した心拍数決定部55と同様であるため、その説明を省略する。
 第1心拍数および第2心拍数は、生体10が平静にしている状態では、同じような数値となるが、I系列用チャンネルおよびQ系列用チャンネルのバランスが悪い場合、すなわち身体の動作等がある場合、片方のチャンネルが異常な数値を示すことがある。したがって、心拍数決定部55aに備えられた表示心拍数決定部55bにおける、表示用の心拍数を決定する方法としては、例えば、第1心拍数および第2心拍数の数値が数パーセント以内の誤差であれば、これらの平均値を心拍数としてもよい。または、第1心拍数と第2心拍数とを比較して、数値の高い方を心拍数としてもよい。
 次に、第1呼吸自己相関関数処理部64iで算出された呼吸数である第1呼吸数、および、第2呼吸自己相関関数処理部64qで算出された呼吸数である第2呼吸数は、それぞれ当該呼吸数を示す信号として(図示しない)呼吸数決定部65aに入力される。具体的には、第1呼吸数は、呼吸数決定部65aに備えられた第1呼吸数決定部65iに入力される。また、第2呼吸数は、同じく呼吸数決定部65aに備えられた第2呼吸数決定部65qに入力される。なお、第1呼吸数および第2呼吸数の決定方法については、実施形態1で説明した呼吸数決定部65と同様であるため、その説明を省略する。
 また、呼吸数決定部65aに備えられた表示呼吸数決定部65bにおける、表示用の呼吸数を決定する方法は、上記表示用の心拍数を決定する方法と同様であるため、その説明を省略する。
 <自己相関関数を用いた心拍成分に関する周期性の判別方法>
 次に、図10の(a)および(b)を用いて、自己相関関数を用いた心拍成分に関する周期性の判別方法について説明する。図10の(a)は、第1心拍自己相関関数処理部54iおよび第2心拍自己相関関数処理部54qにおいてサンプリングされた心拍信号に関して、当該心拍信号の振幅値とサンプル数との関係を示すグラフである。図10の(b)は、第1心拍自己相関関数処理部54iおよび第2心拍自己相関関数処理部54qにおいてサンプリングされた心拍信号に関して、自己相関関数の相関係数とサンプル数との関係を示すグラフである。
 一般に、心拍信号および呼吸信号等の生体信号は、有限の平均パワーをもつ信号と見做せることから、当該信号をx〔n〕と表現する。nはサンプリングによるサンプル番号で整数である。このとき自己相関関数rxx〔l〕は式(1)によって算出される。
Figure JPOXMLDOC01-appb-M000001
 l:サンプルの遅れ(整数)
 N1:整数
 また、周期がLである周期信号の自己相関関数rxx〔l〕は式(2)によって算出される。
Figure JPOXMLDOC01-appb-M000002
 このとき、周期Lが未知な場合には、式(2)の計算をそのまま行うことができない。したがって、式(1)において、N1を、想定される周期に対して十分大きな整数にとり、式(1)を、
Figure JPOXMLDOC01-appb-M000003
に置き換えて算出し、算出結果を、式(2)で定義される自己相関関数rxx〔l〕の近似として用いる。x〔n〕が周期関数であれば、式(3)で算出した自己相関関数の相関係数は、l=0以外に、等間隔で、L、2L、3L、4L…においてピーク値をとる。しかもl=0におけるピーク値とほぼ同じ正のピーク値を持てば、x〔n〕は周期信号であると判定でき、この場合の整数Lに基づいて、周期信号の周期が推定される。
 DSP502において、上記方法によって算出された自己相関関数を用いて、心拍成分に関する周期性を判別する場合、以下のようになる。
 ここで、図10の(a)のグラフでは、縦軸は振幅値であり、図10の(b)のグラフでは、縦軸は相関係数である。両グラフとも、横軸は心拍信号のサンプル数である。振幅値および相関係数ともその値が大きいほど、過去の信号との相関性が高いことを示す。また、式(3)における、n=l=0以外の最初のピーク値(振幅値および相関係数)に対応するサンプル数に基づいて、自己相関関数の周期が算出される。
 図10の(b)のグラフでは、I系列用およびQ系列用の2つの自己相関関数が示されている。前記2つの自己相関関数ともに、横軸のサンプル数がL=46のときに相関係数が最初のピーク値となっている。ここで、サンプリング周波数は、実施形態1および2と同様に50Hzであることから、1サンプル当たり20m秒となる。また、本実施形態では、式(3)に示されたN1=128回のサンプリングを行っていることから、トータルで、20m秒×128=2.56秒のサンプリング時間となる。
 したがって、上記2つの自己相関関数の周期Tは、ともに46×20m秒=0.92秒となる。周期Tを心拍信号の周波数に換算すれば1.087Hz、1分間の心拍数に換算すれば65bpmとなる。
 なお、I系列用およびQ系列用の2つの呼吸信号についても、同様に、I系列用およびQ系列用の2つの自己相関関数を用いて、2つの相関係数のピーク値からその周期性を判別する。なお、通常、呼吸は4秒に1回程度であることから、呼吸信号のサンプリング回数が512回程度必要になるところが、心拍信号の場合と異なる。
 <効果>
 以上のように、本実施形態によれば、高周波装置3は、第1心拍数および第2心拍数を用いて、最終的に表示器90に表示する心拍数を決定する。例えば、第1心拍数および第2心拍数の数値の誤差が数パーセント以内であれば、生体10が静止している(すなわち、生体10が心拍数の算出に適した状態にある)と見做して心拍数を決定することができる。呼吸数の決定についても同様の方法を採ることができる。そのため、高周波装置3は、最終的に精度の高い心拍数および呼吸数を表示器90に表示させることができる。
 また、高周波装置3は、第1心拍数および第2心拍数を、それぞれ位相が90度異なる2つのデジタル信号(Iデジタル信号46si、Qデジタル信号46sq)に基づいて決定する。ここで、コサイン波であるIデジタル信号46siが振幅の小さい信号であっても、Qデジタル信号46sqはサイン波であるため振幅の大きな信号となり、両信号は位相の変化に対して補償し合う関係となる。したがって、生体10に比較的大きな体動があった場合でも、少なくともいずれか一方の心拍数は精度を高く保っている。第1呼吸数と第2呼吸数との関係についても同様である。そのため、高周波装置3は、生体10に比較的大きな体動があった場合においても、精度の高い心拍数および呼吸数を算出することができる。
 具体的には、図10の(a)に示すように、サンプリングされたI系列用の心拍信号とQ系列用の心拍信号のレベルは異なっている(Q系列用の心拍信号がI系列用の心拍信号に比して信号レベルが小さくなっている)。したがって、例えば、Q系列用の心拍信号の自己相関関数のみによる心拍成分および呼吸成分の周期性の判別では、判別の精度が低下してしまう。しかし、本実施形態のように、I系列用の心拍信号およびQ系列用の心拍信号の双方を用いることにより、レベルの高い心拍信号でも抽出することができ、それに伴って、自己相関関数を用いてI系列用とQ系列用の両方で、周期性判別ができる。その結果、高周波装置3における周期性の判別精度を高くすることができる。一例として、図14の(b)では、自己相関関数から求めた周期性は、I系列用の心拍信号およびQ系列用の心拍信号の双方でピーク値L=46となり、1分間の心拍数65bpmが双方で得られる。
 〔実施形態4〕
 本発明の他の実施形態について、図11および図12に基づいて説明すれば、以下のとおりである。なお、説明の便宜上、前記実施形態にて説明した部材と同じ機能を有する部材については、同じ符号を付記し、その説明を省略する。
 本実施形態に係る高周波装置4は、信号処理回路400~402に代えて、信号処理回路403を備えている点で、実施形態1~3に係る高周波装置1~3と異なる。
 <信号処理回路403による信号処理>
 以下に、図11および図12を用いて、信号処理回路403による信号処理について説明する。図11は、信号処理回路403の概略構成を示すブロック図である。図12は、DSP503および表示器90の概略構成を示すブロック図である。図11に示すように、信号処理回路403は、I心拍デジタル信号生成部49a、I呼吸デジタル信号生成部49b、Q心拍デジタル信号生成部49c、Q呼吸デジタル信号生成部49dおよびDSP503を備えている。
 (各デジタル信号生成部49a~49dによる信号処理)
 図11に示すように、Iベースバンド信号33is(図7参照)は入力端子41bに入力され、Qベースバンド信号33qs(図7参照)は入力端子41cに入力される。そして、入力されたIベースバンド信号33isは、I系列のI心拍アナログ信号48iaとI呼吸アナログ信号48ibとに2分配される。また、入力されたQベースバンド信号33iqは、Q系列のQ心拍アナログ信号48qaとQ呼吸アナログ信号48qbとに2分配される。
 ここで、I心拍アナログ信号48iaおよびQ心拍アナログ信号48qaは、心拍成分および体動成分を含むアナログ信号であり、I呼吸アナログ信号48ibおよびQ呼吸アナログ信号48qbは、呼吸成分を含むアナログ信号である。
 そして、I心拍アナログ信号48iaはI心拍デジタル信号生成部49aで信号処理され、I呼吸アナログ信号48ibはI呼吸デジタル信号生成部49bで信号処理される。また、Q心拍アナログ信号48qaはQ心拍デジタル信号生成部49cで信号処理され、Q呼吸アナログ信号48qbはQ呼吸デジタル信号生成部49dで信号処理される。各デジタル信号生成部49a~49dは、デジタル信号生成部40a(図3参照)と同一の構成ものが、I系列用およびQ系列用として備わっているだけなので、各デジタル信号生成部49a~49dでの信号処理についての説明は省略する。
 I心拍デジタル信号生成部49aから出力されたI心拍デジタル信号46sia、および、I呼吸デジタル信号生成部49bから出力されたI呼吸デジタル信号46sibは、それぞれDSP503で信号処理される。また、Q心拍デジタル信号生成部49cから出力されたQ心拍デジタル信号46sqa、および、Q呼吸デジタル信号生成部49dから出力されたQ呼吸デジタル信号46sqbも、それぞれDSP503で信号処理される。なお、I系列用のデジタル信号46siaおよび46sibはコサイン波となっており、Q系列用のデジタル信号46sqaおよび46sqbはサイン波となっている。
 (DSP503による信号処理)
 図12に示すように、I心拍デジタル信号46siaは入力端子52iに入力される。そして、入力端子52iから出力された後、第1のIデジタル信号58aiと第2のIデジタル信号58biとに2分配される。その後、第1のIデジタル信号58aiは、第1心拍信号抽出部53iに入力される。また、第2のIデジタル信号58biは、そのまま、時間の経過に伴う体動振幅の変化を表す体動波形信号71bとして出力され、表示器90への入力信号となる。
 また、I呼吸デジタル信号46sibは入力端子62iに入力された後、第3のIデジタル信号58ciとして第1呼吸信号抽出部63iに出力される。Q心拍デジタル信号46sqaは入力端子52qに入力された後、第1のQデジタル信号58aqとして第2心拍信号抽出部53qに出力される。Q呼吸デジタル信号46sqbは入力端子62qに入力された後、第2のQデジタル信号58bqとして第2呼吸信号抽出部63qに出力される。
 第1のIデジタル信号58ai、第3のIデジタル信号58ci、第1のQデジタル信号58aqおよび第2のQデジタル信号58bqに基づいて、各自己相関関数の周期を算出するまでの信号処理の流れについては、実施形態1~3と同様であるため、その説明を省略する。また、第1のIデジタル信号58aiに基づいて心拍波形信号71aを得る信号処理、第3のIデジタル信号58ciに基づいて呼吸波形信号71cを得る信号処理についても、実施形態1~3と同様であるため、その説明を省略する。
 次に、心拍数決定部55aにおいて表示用の心拍数を決定する方法、および、呼吸数決定部65aにおいて表示用の呼吸数を決定する方法とも、実施形態3と同様であるため、その説明を省略する。
 <効果>
 以上のように、本実施形態によれば、高周波装置4は、心拍成分および体動成分を含む生体信号(以下、「心拍成分信号」とする)と呼吸成分を含む生体信号(以下、「呼吸成分信号」)とに分けて、アナログ信号処理およびデジタル信号処理する。したがって、心拍成分信号および呼吸成分信号のダイナミックレンジを、アナログ信号処理およびデジタル信号処理の各段階で大きくすることができる。そのため、各心拍自己相関関数処理部および各呼吸自己相関関数処理部は、それぞれ高い精度で、各心拍数および各呼吸数を算出することができる。
 具体的には、心拍成分信号については、呼吸成分の影響を低減し、呼吸成分信号については、体動成分の影響を低減することが可能である。心拍成分信号については、体動の少ない安静時および睡眠時に効果を発揮し、呼吸成分信号については、体動が大きい非静止時、および、睡眠時でも眠りに入る前の未だ体動があるときに効果を発揮する。
 また、高周波装置4は、上記各心拍数(第1心拍数、第2心拍数)および各呼吸数(第1呼吸数、第2呼吸数)に基づいて、表示用の心拍数および呼吸数を決定する。そのため、高周波装置4は、最終的により精度の高い心拍数および呼吸数を表示器90に表示することができるとともに、生体10に比較的大きな体動があった場合においても、精度の高い心拍数および呼吸数を算出することができる。
 〔ソフトウェアによる実現例〕
 高周波装置1~4の制御ブロック(特に電波レーダ部5aおよび5b、信号処理回路400~403)は、集積回路(ICチップ)等に形成された論理回路(ハードウェア)によって実現してもよいし、CPU(Central Processing Unit)を用いてソフトウェアによって実現してもよい。
 後者の場合、高周波装置1~4は、各機能を実現するソフトウェアであるプログラムの命令を実行するCPU、上記プログラムおよび各種データがコンピュータ(またはCPU)で読み取り可能に記録されたROM(Read Only Memory)または記憶装置(これらを「記録媒体」と称する)、上記プログラムを展開するRAM(Random Access Memory)などを有している。そして、コンピュータ(またはCPU)が上記プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記プログラムは、該プログラムを伝送可能な任意の伝送媒体(通信ネットワークや放送波等)を介して上記コンピュータに供給されてもよい。なお、本発明は、上記プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。
 〔まとめ〕
 本発明の態様1に係る高周波装置(1、2、3、4)は、
 生体に現れる動きである生体現象を表すデジタル信号から、上記デジタル信号の特定の周波数成分を示す生体信号を抽出する生体信号抽出部(心拍信号抽出部53、呼吸信号抽出部63)と、
 上記生体信号抽出部により抽出された上記生体信号をサンプリングし、上記サンプリングされた所定数の上記生体信号に基づいて自己相関関数を算出し、相関係数のピーク値から上記自己相関関数の周期性を判別することにより、上記生体現象を表す生体情報を算出する自己相関関数処理部(心拍自己相関関数処理部54、呼吸自己相関関数処理部63)と、を備えている。
 上記の構成によれば、高周波装置は、アナログフィルタで必要な帯域のみ抽出し、かつデジタルフィルタでも必要な帯域のみを通過させる。具体的には、高周波装置は、アナログフィルタで必要な帯域のみ抽出して不要信号を除き、且つ増幅した後、アナログデジタル変換することによって、心拍または呼吸帯域でのダイナミックレンジを大きくする。そして、デジタルフィルタにおいてデジタル信号のSN比を大きくする。そのため、自己相関関数から検出される相関係数のピーク値の検知精度を高くすることができる。
 また、高周波装置は、生体信号をサンプリングして得られた自己相関関数を用いて相関係数のピーク値の検出し、当該自己相関関数の周期を判別することで、生体情報を算出する。そのため、生体が動いている場合または不要反射信号等のノイズ成分が多い場合においても、生体情報を精度高く算出することができる。
 本発明の態様2に係る高周波装置は、上記態様1において、
 上記生体信号抽出部は、上記生体信号を抽出するために自装置から上記生体に向けて送信される送信信号が上記生体で直接反射した直接反射信号(反射信号12)および、上記送信信号が物体での反射を介して上記生体で反射した間接反射信号(反射信号12)から、上記生体信号を抽出してもよい。
 一般に、間接反射信号については、波形が変形したり、ノイズが付加したりするため、当該間接反射信号から生体情報を算出することは困難とされている。
 その点、上記の構成によれば、直接反射信号および間接反射信号から、生体信号抽出部を用いて生体信号を抽出することができる。そのため、間接反射信号についても、生体信号の自己相関関数を用いて生体情報を算出することができる。
 また、寝室等の比較的狭い空間に生体が存在する場合には、当該空間内に間接反射信号が多数存在する。そのため、高周波装置は、直接反射信号を捕捉できない場合においても、上記間接反射信号を用いて当該空間内に存在する生体の生体情報を算出することができる。
 さらに、上記空間に生体が存在する場合においては、例えば、生体がうつ伏せで横たわっている等、生体の所望の部位に直接送信信号を照射することが困難な状況であっても、壁等で送信信号を反射させることによって、当該部位に間接的に送信信号を照射することができる。そのため、間接反射信号を用いることで、生体の姿勢に拘らず生体情報を算出することができる。
 本発明の態様3に係る高周波装置は、上記態様1または2において、
 上記生体信号抽出部は、特定の周波数の心拍を示す心拍信号を抽出する心拍信号抽出部(53)と、特定の周波数の呼吸を示す呼吸信号を抽出する呼吸信号抽出部(63)と、を有し、
 上記自己相関関数処理部は、所定時間における心拍数を算出する心拍自己相関関数処理部(54)と、所定時間における呼吸数を算出する呼吸自己相関関数処理部(64)と、を有していてもよい。
 上記の構成によれば、高周波装置は、生体が動いている場合または不要反射信号等のノイズ成分が多い場合においても、心拍信号および呼吸信号における自己相関関数の周期性を判別することにより、所定時間における心拍数および呼吸数を算出することができる。そのため、前記心拍数および呼吸数の算出に関して、検知感度の向上した高周波装置を実現することができる。
 本発明の態様4に係る高周波装置(2)は、上記態様3において、
 上記生体現象に関するアナログ信号(第1アナログ信号48a、第2アナログ信号48b)が入力されることにより、心拍成分および体動成分に関する第1生体デジタル信号(46sa)および呼吸成分に関する第2生体デジタル信号(46sb)を生成するデジタル信号生成部(40b)をさらに備え、
 上記心拍信号抽出部は、上記デジタル信号生成部から出力された上記第1生体デジタル信号に基づいて上記心拍信号を抽出し、
 上記呼吸信号処理部は、上記デジタル信号生成部から出力された上記第2生体デジタル信号に基づいて上記呼吸信号を抽出してもよい。
 上記の構成によれば、デジタル信号生成部は、第1生体デジタル信号および第2生体デジタル信号のそれぞれについてダイナミックレンジを大きくすることができる。そのため、心拍信号抽出部および呼吸信号抽出部は、高い感度・精度で心拍信号および呼吸信号を抽出することができる。
 本発明の態様5に係る高周波装置(3)は、上記態様3において、
 上記生体現象に関するIアナログ信号(31is)と、上記Iアナログ信号と位相が90度異なるQアナログ信号(31qs)とを出力する電波レーダ部(5b)と、
 上記心拍数を決定する心拍数決定部(55a)と、
 上記呼吸数を決定する呼吸数決定部(55b)と、をさらに備え、
 上記デジタル信号生成部(第1デジタル信号生成部40ci、第2デジタル信号生成部40cq)は、上記電波レーダ部から出力された上記Iアナログ信号に基づいて、心拍成分と呼吸成分と体動成分とに関するIデジタル信号(46si)を生成し、かつ、上記電波レーダ部から出力された上記Qアナログ信号に基づいて、心拍成分と呼吸成分と体動成分とに関するQデジタル信号(46sq)を生成し、
 上記心拍自己相関関数処理部(第1心拍自己相関関数処理部54i、第2心拍自己相関関数処理部54q)は、上記Iデジタル信号に基づいて所定時間における第1心拍数を算出し、かつ、上記Qデジタル信号に基づいて所定時間における第2心拍数を算出し、
 上記呼吸自己相関関数処理部(第1呼吸自己相関関数処理部64i、第2呼吸自己相関関数処理部64q)は、上記Iデジタル信号に基づいて所定時間における第1呼吸数を算出し、かつ、上記Qデジタル信号に基づいて所定時間における第2呼吸数を算出し、
 上記心拍数決定部は、上記心拍自己相関関数処理部から出力された上記第1心拍数および上記第2心拍数に基づいて上記心拍数を決定し、
 上記呼吸数決定部は、上記呼吸自己相関関数処理部から出力された上記第1呼吸数および上記第2呼吸数に基づいて上記呼吸数を決定してもよい。
 上記の構成によれば、高周波装置は、第1心拍数および第2心拍数を用いて、最終的な心拍数を決定する。例えば、第1心拍数および第2心拍数の数値の誤差が数パーセント以内であれば、生体が静止している(すなわち、生体が心拍数の算出に適した状態にある)と見做して心拍数を決定することができる。呼吸数の決定についても同様の方法を採ることができる。そのため、高周波装置は、最終的に精度の高い心拍数および呼吸数を決定することができる。
 また、高周波装置は、第1心拍数および第2心拍数を、それぞれ位相が90度異なる2つのデジタル信号(Iデジタル信号、Qデジタル信号)に基づいて決定する。ここで、コサイン波であるIデジタル信号が振幅の小さい信号であっても、Qデジタル信号はサイン波であるため振幅の大きな信号となり、両信号は位相の変化に対して補償し合う関係となる。したがって、生体に比較的大きな体動があった場合でも、少なくともいずれか一方の心拍数は精度を高く保っている。第1呼吸数と第2呼吸数との関係についても同様である。そのため、高周波装置は、生体に比較的大きな体動があった場合においても、精度の高い心拍数および呼吸数を算出することができる。
 なお、上記態様5に係る高周波装置(4)は、
 上記デジタル信号生成部(I心拍デジタル信号生成部49a、I呼吸デジタル信号生成部49b、Q心拍デジタル信号生成部49c、Q呼吸デジタル信号生成部49d)は、上記電波レーダ部から出力された上記Iアナログ信号(I心拍アナログ信号48ia、I呼吸アナログ信号48ib)に基づいて、心拍成分および体動成分に関するI心拍デジタル信号(46sia)と、呼吸成分に関するI呼吸デジタル信号(46sib)とを生成し、かつ、上記電波レーダ部から出力された上記Qアナログ信号(Q心拍アナログ信号48qa、Q呼吸アナログ信号48qb)に基づいて、心拍成分および体動成分に関するQ心拍デジタル信号(46sqa)と、呼吸成分に関するQ呼吸デジタル信号(46sqb)とを生成し、
 上記心拍自己相関関数処理部は、上記I心拍デジタル信号に基づいて上記第1心拍数を算出し、かつ、上記Q心拍デジタル信号に基づいて上記第2心拍数を算出し、
 上記呼吸自己相関関数処理部は、上記I呼吸デジタル信号に基づいて上記第1呼吸数を算出し、かつ、上記Q呼吸デジタル信号に基づいて上記第2呼吸数を算出してもよい。
 上記の構成によれば、高周波装置は、心拍成分信号と呼吸成分信号とに分けて、アナログ信号処理およびデジタル信号処理する。したがって、心拍成分信号および呼吸成分信号のダイナミックレンジを、アナログ信号処理およびデジタル信号処理の各段階で大きくすることができる。そのため、各心拍自己相関関数処理部および各呼吸自己相関関数処理部は、それぞれ高い精度で、各心拍数および各呼吸数を算出することができる。
 また、高周波装置は、第1心拍数および第2心拍数、第1呼吸数および第2呼吸数に基づいて、表示用の心拍数および呼吸数を決定する。そのため、高周波装置は、最終的により精度の高い心拍数および呼吸数を算出することができるとともに、生体に比較的大きな体動があった場合においても、精度の高い心拍数および呼吸数を算出することができる。
 本発明の各態様に係る高周波装置(1、2、3、4)は、コンピュータによって実現してもよく、この場合には、コンピュータを上記高周波装置が有する各部(ソフトウェア要素に限る)として動作させることにより上記高周波装置をコンピュータにて実現させる高周波装置の制御プログラム、およびそれを記録したコンピュータ読み取り可能な記録媒体も、本発明の範疇に入る。
 本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。さらに、各実施形態にそれぞれ開示された技術的手段を組み合わせることにより、新しい技術的特徴を形成することができる。
 本発明の高周波装置は、部屋での生体信号検知、自動車の運転時における生体信号検知等に有用であるとともに、福祉介護および医療等に応用することができる。
 1、2、3、4 高周波装置
 5a、5b   電波レーダ部
 10      生体
 11      送信信号
 12      反射信号(直接反射信号、間接反射信号)
 31is    Iアナログ信号
 31qs    Qアナログ信号
 40a、40b デジタル信号生成部
 40ci    第1デジタル信号生成部(デジタル信号生成部)
 40cq    第2デジタル信号生成部(デジタル信号生成部)
 46sa    第1生体デジタル信号
 46sb    第2生体デジタル信号
 46si    Iデジタル信号
 46sq    Qデジタル信号
 46sia   I心拍デジタル信号
 46sib   I呼吸デジタル信号
 46sqa   Q心拍デジタル信号
 46sqb   Q呼吸デジタル信号
 48a     第1アナログ信号(アナログ信号)
 48b     第2アナログ信号(アナログ信号)
 48ia    I心拍アナログ信号(Iアナログ信号)
 48ib    I呼吸アナログ信号(Iアナログ信号)
 48qa    Q心拍アナログ信号(Qアナログ信号)
 48qb    Q呼吸アナログ信号(Qアナログ信号)
 49a     I心拍デジタル信号生成部(デジタル信号生成部)
 49b     I呼吸デジタル信号生成部(デジタル信号生成部)
 49c     Q心拍デジタル信号生成部(デジタル信号生成部)
 49d     Q呼吸デジタル信号生成部(デジタル信号生成部)
 53      心拍信号抽出部
 53i     第1心拍信号抽出部(心拍信号抽出部)
 53q     第2心拍信号抽出部(心拍信号抽出部)
 54      心拍自己相関関数処理部
 54i     第1心拍自己相関関数処理部(心拍自己相関関数処理部)
 54q     第2心拍自己相関関数処理部(心拍自己相関関数処理部)
 55、55a  心拍数決定部
 63      呼吸信号抽出部
 63i     第1呼吸信号抽出部(呼吸信号抽出部)
 63q     第2呼吸信号抽出部(呼吸信号抽出部)
 64      呼吸自己相関関数処理部
 64i     第1呼吸自己相関関数処理部(呼吸自己相関関数処理部)
 64q     第2呼吸自己相関関数処理部(呼吸自己相関関数処理部)
 65、65a  呼吸数決定部

Claims (5)

  1.  生体に現れる動きである生体現象を表すデジタル信号から、上記デジタル信号の特定の周波数成分を示す生体信号を抽出する生体信号抽出部と、
     上記生体信号抽出部により抽出された上記生体信号をサンプリングし、上記サンプリングされた所定数の上記生体信号に基づいて自己相関関数を算出し、相関係数のピーク値から上記自己相関関数の周期性を判別することにより、上記生体現象に関する生体情報を算出する自己相関関数処理部と、を備えていることを特徴とする高周波装置。
  2.  上記生体信号抽出部は、上記生体信号を抽出するために自装置から上記生体に向けて送信される送信信号が上記生体で直接反射した直接反射信号および、上記送信信号が物体での反射を介して上記生体で反射した間接反射信号から、上記生体信号を抽出することを特徴とする請求項1に記載の高周波装置。
  3.  上記生体信号抽出部は、特定の周波数の心拍を示す心拍信号を抽出する心拍信号抽出部と、特定の周波数の呼吸を示す呼吸信号を抽出する呼吸信号抽出部と、を有し、
     上記自己相関関数処理部は、所定時間における心拍数を算出する心拍自己相関関数処理部と、所定時間における呼吸数を算出する呼吸自己相関関数処理部と、を有することを特徴とする請求項1または2に記載の高周波装置。
  4.  上記生体現象に関するアナログ信号が入力されることにより、心拍成分および体動成分に関する第1生体デジタル信号および呼吸成分に関する第2生体デジタル信号を生成するデジタル信号生成部をさらに備え、
     上記心拍信号抽出部は、上記デジタル信号生成部から出力された上記第1生体デジタル信号に基づいて上記心拍信号を抽出し、
     上記呼吸信号処理部は、上記デジタル信号生成部から出力された上記第2生体デジタル信号に基づいて上記呼吸信号を抽出することを特徴とする請求項3に記載の高周波装置。
  5.  上記生体現象に関するIアナログ信号と、上記Iアナログ信号と位相が90度異なるQアナログ信号とを出力する電波レーダ部と、
     上記心拍数を決定する心拍数決定部と、
     上記呼吸数を決定する呼吸数決定部と、をさらに備え、
     上記デジタル信号生成部は、上記電波レーダ部から出力された上記Iアナログ信号に基づいて、心拍成分と呼吸成分と体動成分とに関するIデジタル信号を生成し、かつ、上記電波レーダ部から出力された上記Qアナログ信号に基づいて、心拍成分と呼吸成分と体動成分とに関するQデジタル信号を生成し、
     上記心拍自己相関関数処理部は、上記Iデジタル信号に基づいて所定時間における第1心拍数を算出し、かつ、上記Qデジタル信号に基づいて所定時間における第2心拍数を算出し、
     上記呼吸自己相関関数処理部は、上記Iデジタル信号に基づいて所定時間における第1呼吸数を算出し、かつ、上記Qデジタル信号に基づいて所定時間における第2呼吸数を算出し、
     上記心拍数決定部は、上記心拍自己相関関数処理部から出力された上記第1心拍数および上記第2心拍数に基づいて上記心拍数を決定し、
     上記呼吸数決定部は、上記呼吸自己相関関数処理部から出力された上記第1呼吸数および上記第2呼吸数に基づいて上記呼吸数を決定することを特徴とする請求項3に記載の高周波装置。
PCT/JP2015/077834 2014-11-28 2015-09-30 高周波装置 WO2016084473A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US15/516,713 US10582877B2 (en) 2014-11-28 2015-09-30 High-frequency device
JP2016561437A JP6290448B2 (ja) 2014-11-28 2015-09-30 高周波装置
CN201580065123.6A CN106999105B (zh) 2014-11-28 2015-09-30 高频装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014-242161 2014-11-28
JP2014242161 2014-11-28

Publications (1)

Publication Number Publication Date
WO2016084473A1 true WO2016084473A1 (ja) 2016-06-02

Family

ID=56074052

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/077834 WO2016084473A1 (ja) 2014-11-28 2015-09-30 高周波装置

Country Status (4)

Country Link
US (1) US10582877B2 (ja)
JP (2) JP6290448B2 (ja)
CN (1) CN106999105B (ja)
WO (1) WO2016084473A1 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017225834A (ja) * 2014-11-28 2017-12-28 シャープ株式会社 高周波装置
CN108697352A (zh) * 2017-06-29 2018-10-23 深圳和而泰智能控制股份有限公司 生理信息测量方法及生理信息监测装置、设备
JP2018187129A (ja) * 2017-05-09 2018-11-29 シャープ株式会社 計測装置、計測方法、および計測プログラム
JP2020513929A (ja) * 2017-01-09 2020-05-21 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 聴診器及び被験体により生成される呼吸音を処理する方法
KR102381262B1 (ko) * 2021-12-17 2022-04-01 (주)디지탈엣지 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법
JP7146142B1 (ja) * 2021-12-17 2022-10-03 三菱電機株式会社 バイタル測定装置、バイタル測定方法及びバイタル測定システム
JP7166687B1 (ja) 2022-02-07 2022-11-08 トミタテクノロジー・ジャパン株式会社 生体信号処理装置

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107928250A (zh) * 2017-11-22 2018-04-20 宁波德葳智能科技有限公司 一种无接触毫米波智能婴儿摇篮
KR102117748B1 (ko) * 2018-03-20 2020-06-01 한양대학교 산학협력단 생체 정보 측정용 신호의 품질 평가 장치 및 방법
CN108272447A (zh) * 2018-03-27 2018-07-13 厦门金网科技有限公司 一种基于uwb技术的呼吸和心跳速率检测系统
CN108888271B (zh) * 2018-07-25 2020-05-22 佛山市丈量科技有限公司 一种生理参数测量系统及配备该测量系统的智能座椅
WO2020042444A1 (zh) * 2018-08-31 2020-03-05 深圳迈睿智能科技有限公司 人体存在探测器及其人体存在探测方法
DE102018218215A1 (de) * 2018-10-24 2020-04-30 Robert Bosch Gmbh Insassenüberwachungssystem für ein Fahrzeug
WO2020105836A1 (ko) * 2018-11-20 2020-05-28 한국과학기술원 연속파 도플러 레이더 기반 고정확도 실시간 심장 박동 수 측정 방법 및 그 레이더 시스템
JP7224163B2 (ja) * 2018-12-11 2023-02-17 日本無線株式会社 心拍周期測定装置及び心拍周期測定プログラム
JP7099353B2 (ja) * 2019-02-14 2022-07-12 コニカミノルタ株式会社 生体情報測定システム
CN110215214A (zh) * 2019-05-29 2019-09-10 深圳和而泰家居在线网络科技有限公司 对床的监测方法、监测装置及监测设备
WO2021095893A1 (ko) * 2019-11-11 2021-05-20 주식회사 제이씨에프테크놀러지 생체신호 측정용 레이더에서 레이더 수집신호의 랜덤노이즈 제거방법 및 그 장치
DE112020005588B9 (de) 2020-01-08 2024-02-15 Mitsubishi Electric Corporation Biometrische Authentifizierungsvorrichtung, biometrisches Authentifizierungsverfahren und biometrisches Authentifizierungsprogramm
CN113925481A (zh) * 2021-10-20 2022-01-14 友达光电股份有限公司 测量生理信号方法及生理信号测量装置
CN114287885B (zh) * 2021-12-28 2023-12-08 深圳数联天下智能科技有限公司 一种人体体征监测方法、装置、系统以及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010120493A (ja) * 2008-11-19 2010-06-03 Panasonic Corp 生体信号検知装置
JP2014176427A (ja) * 2013-03-13 2014-09-25 Kansai Univ データ解析装置、及びデータ解析プログラム
JP2014215200A (ja) * 2013-04-26 2014-11-17 沖電気工業株式会社 特徴量算出装置、特徴量算出方法およびプログラム

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2654095C (en) * 2006-06-01 2015-12-22 Biancamed Ltd. Apparatus, system, and method for monitoring physiological signs
CN100488453C (zh) * 2007-05-10 2009-05-20 浙江大学 基于射频识别的非接触式生物特征识别系统
CN101324666A (zh) * 2007-06-16 2008-12-17 电子科技大学 隐蔽目标生命迹象探测方法及隐蔽目标探测装置
WO2009009722A2 (en) * 2007-07-12 2009-01-15 University Of Florida Research Foundation, Inc. Random body movement cancellation for non-contact vital sign detection
US8814805B2 (en) * 2007-12-07 2014-08-26 University Of Florida Research Foundation, Inc. Complex signal demodulation and angular demodulation for non-contact vital sign detection
US8454528B2 (en) * 2008-04-03 2013-06-04 Kai Medical, Inc. Non-contact physiologic motion sensors and methods for use
US20100152600A1 (en) * 2008-04-03 2010-06-17 Kai Sensors, Inc. Non-contact physiologic motion sensors and methods for use
CN102215746B (zh) * 2008-09-15 2017-03-01 新加坡南洋理工大学 用于检测心跳和/或呼吸的方法
JP5469574B2 (ja) * 2010-09-22 2014-04-16 学校法人北里研究所 電気的インピーダンス断層像測定装置。
US8740793B2 (en) * 2011-08-29 2014-06-03 General Electric Company Radar based systems and methods for monitoring a subject
CN102715920B (zh) * 2012-07-06 2013-10-16 电子科技大学 人体目标生命特征检测方法
CN102988051B (zh) * 2012-12-13 2014-07-02 中国人民解放军第四军医大学 用于计算机操作者健康的监测装置
CN103110422B (zh) * 2012-12-18 2014-10-15 中国人民解放军第四军医大学 基于生物雷达检测的呼吸和心跳实时分离方法
WO2016084473A1 (ja) * 2014-11-28 2016-06-02 シャープ株式会社 高周波装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010120493A (ja) * 2008-11-19 2010-06-03 Panasonic Corp 生体信号検知装置
JP2014176427A (ja) * 2013-03-13 2014-09-25 Kansai Univ データ解析装置、及びデータ解析プログラム
JP2014215200A (ja) * 2013-04-26 2014-11-17 沖電気工業株式会社 特徴量算出装置、特徴量算出方法およびプログラム

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017225834A (ja) * 2014-11-28 2017-12-28 シャープ株式会社 高周波装置
JP2020513929A (ja) * 2017-01-09 2020-05-21 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 聴診器及び被験体により生成される呼吸音を処理する方法
US10682117B2 (en) 2017-01-09 2020-06-16 Koninklijke Philips N.V. Stethoscope apparatus and a method of processing breathing sounds produced by a subject
JP2018187129A (ja) * 2017-05-09 2018-11-29 シャープ株式会社 計測装置、計測方法、および計測プログラム
CN108697352A (zh) * 2017-06-29 2018-10-23 深圳和而泰智能控制股份有限公司 生理信息测量方法及生理信息监测装置、设备
WO2019000338A1 (zh) * 2017-06-29 2019-01-03 深圳和而泰智能控制股份有限公司 生理信息测量方法及生理信息监测装置、设备
KR102381262B1 (ko) * 2021-12-17 2022-04-01 (주)디지탈엣지 레이더 시스템 및 그에서 수행되는 생체 신호 탐지방법
JP7146142B1 (ja) * 2021-12-17 2022-10-03 三菱電機株式会社 バイタル測定装置、バイタル測定方法及びバイタル測定システム
WO2023112295A1 (ja) * 2021-12-17 2023-06-22 三菱電機株式会社 バイタル測定装置、バイタル測定方法及びバイタル測定システム
JP7166687B1 (ja) 2022-02-07 2022-11-08 トミタテクノロジー・ジャパン株式会社 生体信号処理装置
JP2023114575A (ja) * 2022-02-07 2023-08-18 トミタテクノロジー・ジャパン株式会社 生体信号処理装置

Also Published As

Publication number Publication date
CN106999105B (zh) 2019-10-18
JP2017225834A (ja) 2017-12-28
JPWO2016084473A1 (ja) 2017-06-08
US20170311840A1 (en) 2017-11-02
JP6290448B2 (ja) 2018-03-07
CN106999105A (zh) 2017-08-01
JP6290501B2 (ja) 2018-03-07
US10582877B2 (en) 2020-03-10

Similar Documents

Publication Publication Date Title
JP6290501B2 (ja) 高周波装置
US9993175B2 (en) Pulsed ultra-wideband sensor and the method thereof
US10791965B2 (en) Method and apparatus for wirelessly monitoring repetitive bodily movements
Wang et al. SonarBeat: Sonar phase for breathing beat monitoring with smartphones
US8597196B2 (en) Cardiac signal processing apparatus and cardiac signal processing method
JP5848469B1 (ja) 生体状態検出装置
US8698636B2 (en) Wireless detection apparatus and method
JP2009213881A6 (ja) パルス型超広帯域センサおよびその方法
WO2021035197A1 (en) Wireless communication with enhanced maximum permissible exposure (mpe) compliance based on vital signs detection
CN104394773B (zh) 使用二维阵列数据的移动超声诊断系统以及用于该系统的移动超声诊断探测装置和超声诊断设备
JP5578683B2 (ja) 身体情報測定装置及び身体情報測定方法
JP5672604B2 (ja) 身体情報測定装置及び身体情報測定プログラム
JP2016202516A (ja) 心容積及び心拍出量の推定装置
JP6222454B2 (ja) 生体情報測定装置及び生体情報測定方法
TW201225912A (en) Method for measuring physiological parameters
KR102206238B1 (ko) 생체신호를 검출하는 방법 및 장치
JP2018187129A (ja) 計測装置、計測方法、および計測プログラム
KR20170019785A (ko) 원형편파 헬릭스 안테나를 이용한 고감도 비침습 생체신호 검출 장치 및 방법
EP3207864A1 (en) Sensor information processing apparatus
US8679017B2 (en) Biological information measurement apparatus and method
CN105796051B (zh) 三维生理检测系统及其操作方法
RU2392852C2 (ru) Импульсный сверхширокополосный датчик дистанционного мониторинга дыхания и сердцебиения
US20220330920A1 (en) Ultrasonic diagnostic apparatus and medical image processing apparatus
JP5505981B2 (ja) 磁気共鳴イメージング装置
Song et al. Adaptive frequency compounding for speckle reduction

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: 15862749

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016561437

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15516713

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15862749

Country of ref document: EP

Kind code of ref document: A1