JP7428605B2 - Respiratory heart rate measuring device and respiratory heart rate measuring program - Google Patents
Respiratory heart rate measuring device and respiratory heart rate measuring program Download PDFInfo
- Publication number
- JP7428605B2 JP7428605B2 JP2020114641A JP2020114641A JP7428605B2 JP 7428605 B2 JP7428605 B2 JP 7428605B2 JP 2020114641 A JP2020114641 A JP 2020114641A JP 2020114641 A JP2020114641 A JP 2020114641A JP 7428605 B2 JP7428605 B2 JP 7428605B2
- Authority
- JP
- Japan
- Prior art keywords
- radar signal
- frequency
- heart rate
- subject
- maximum peak
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000000241 respiratory effect Effects 0.000 title claims description 57
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 102
- 238000006243 chemical reaction Methods 0.000 claims description 61
- 238000005259 measurement Methods 0.000 claims description 50
- 230000036387 respiratory rate Effects 0.000 claims description 23
- 238000000034 method Methods 0.000 description 21
- 238000009532 heart rate measurement Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 12
- 239000000284 extract Substances 0.000 description 6
- 238000012854 evaluation process Methods 0.000 description 3
- 238000009531 respiratory rate measurement Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本開示は、レーダを用いて被検体の呼吸数及び心拍数を計測する技術に関する。 The present disclosure relates to a technique for measuring a subject's breathing rate and heart rate using radar.
レーダを用いて被検体の呼吸数及び心拍数を計測する技術が、特許文献1等に開示されている。特許文献1等では、被検体の呼吸及び心拍に由来する被検体の体表面の位置の時間変化を、被検体から反射されたレーダ信号の振幅又は位相の時間変化として検出することにより、被検体と非接触で被検体の呼吸数及び心拍数を計測することができる。
A technique for measuring a subject's breathing rate and heart rate using a radar is disclosed in
特許文献1の呼吸心拍計測処理を図1に示す。まず、一般的な呼吸数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、周波数fB1を有する電力Pの大きなレーダ信号成分をレーダ信号から抽出することができ、被検体の呼吸数を計測することができる。次に、一般的な心拍数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、周波数fH1を有する電力Pの小さなレーダ信号成分をレーダ信号から抽出することができ、被検体の心拍数を計測することができる。
FIG. 1 shows the respiratory heart rate measurement process of
ここで、被検体の呼吸に由来するレーダ信号成分として、基本周波数fB1を有する電力Pの大きなレーダ信号成分のみならず、非正弦波信号又は非線形歪みに起因する高調波周波数fB2、fB3(更なる高調波周波数については省略している。)を有するレーダ信号成分が存在する。しかし、基本周波数fB1を有するレーダ信号成分は考慮されているが、高調波周波数fB2、fB3を有するレーダ信号成分は考慮されていない。そして、被検体の心拍に由来するレーダ信号成分として、基本周波数fH1を有する電力Pの小さなレーダ信号成分のみならず、非正弦波信号又は非線形歪みに起因する高調波周波数fH2、fH3を有するレーダ信号成分が存在する。しかし、基本周波数fH1を有するレーダ信号成分は考慮されているが、高調波周波数fH2、fH3を有するレーダ信号成分は考慮されていない。 Here, the radar signal component originating from the subject's respiration includes not only a radar signal component having a fundamental frequency f B1 and a large power P, but also harmonic frequencies f B2 and f B3 resulting from a non-sinusoidal signal or non-linear distortion. (further harmonic frequencies have been omitted). However, while the radar signal component with fundamental frequency f B1 is taken into account, the radar signal components with harmonic frequencies f B2 , f B3 are not taken into account. As the radar signal component derived from the heartbeat of the subject, not only the radar signal component with a small power P having the fundamental frequency f H1 but also harmonic frequencies f H2 and f H3 resulting from a non-sinusoidal signal or non-linear distortion are used. There is a radar signal component with However, while the radar signal component with the fundamental frequency f H1 is taken into account, the radar signal components with harmonic frequencies f H2 , f H3 are not taken into account.
図1の上段では、被検体の呼吸及び心拍は、通常と同等である。よって、一般的な呼吸数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、基本周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号から抽出することができる。そして、一般的な心拍数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、基本周波数fH1及び高調波周波数fH2、fH3を有するレーダ信号成分をレーダ信号から抽出することができる。つまり、被検体の呼吸に由来する基本・高調波周波数fB1、fB2、fB3を有するレーダ信号成分と、被検体の心拍に由来する基本・高調波周波数fH1、fH2、fH3を有するレーダ信号成分と、を正しく分離することができる。そして、被検体の呼吸数及び心拍数を正確に計測することができる。 In the upper part of FIG. 1, the subject's breathing and heartbeat are normal. Therefore, by applying a bandpass filter having a pass band corresponding to the general breathing rate to the radar signal, radar signal components having the fundamental frequency f B1 and harmonic frequencies f B2 and f B3 are extracted from the radar signal. be able to. Then, by applying a bandpass filter having a passband corresponding to a general heart rate to the radar signal, radar signal components having a fundamental frequency f H1 and harmonic frequencies f H2 and f H3 are extracted from the radar signal. be able to. In other words, radar signal components having fundamental/harmonic frequencies f B1 , f B2 , f B3 derived from the subject's breathing, and fundamental/harmonic frequencies f H1 , f H2 , f H3 originating from the subject's heartbeat are combined. It is possible to correctly separate the radar signal components that have Then, the breathing rate and heart rate of the subject can be accurately measured.
図1の中段では、被検体の呼吸及び心拍は、通常と比べて遅い。よって、一般的な呼吸数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、基本周波数fB1及び高調波周波数fB2を有するレーダ信号成分をレーダ信号から抽出することができるが、高調波周波数fB3を有するレーダ信号成分をレーダ信号から抽出することができない。そして、一般的な心拍数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、基本周波数fH1及び高調波周波数fH2、fH3を有するレーダ信号成分をレーダ信号から抽出することができるが、呼吸に由来する高調波周波数fB3を有するレーダ信号成分をレーダ信号から除去することができない。つまり、被検体の呼吸に由来する高調波周波数fB3を有するレーダ信号成分を、被検体の心拍に由来する基本周波数を有するレーダ信号成分として誤って抽出してしまう。よって、被検体の呼吸数を正確に計測することはできるが、被検体の心拍数を正確に計測することができない。 In the middle part of FIG. 1, the subject's breathing and heartbeat are slower than normal. Therefore, by applying a bandpass filter having a passband corresponding to the general breathing rate to the radar signal, the radar signal component having the fundamental frequency f B1 and the harmonic frequency f B2 can be extracted from the radar signal. However, the radar signal component with harmonic frequency f B3 cannot be extracted from the radar signal. Then, by applying a bandpass filter having a passband corresponding to a general heart rate to the radar signal, radar signal components having a fundamental frequency f H1 and harmonic frequencies f H2 and f H3 are extracted from the radar signal. However, the radar signal component having the harmonic frequency f B3 originating from breathing cannot be removed from the radar signal. In other words, the radar signal component having the harmonic frequency fB3 derived from the subject's breathing is erroneously extracted as the radar signal component having the fundamental frequency originating from the subject's heartbeat. Therefore, although it is possible to accurately measure the respiratory rate of the subject, it is not possible to accurately measure the heart rate of the subject.
図1の下段では、被検体の呼吸及び心拍は、通常と比べて早い。よって、一般的な呼吸数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、基本周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号から抽出することができるが、心拍に由来する基本周波数fH1を有するレーダ信号成分をレーダ信号から除去することができない。そして、一般的な心拍数に応じた通過帯域を有するバンドパスフィルタをレーダ信号に適用することにより、高調波周波数fH2、fH3を有するレーダ信号成分をレーダ信号から抽出することができるが、基本周波数fH1を有するレーダ信号成分をレーダ信号から抽出することができない。つまり、被検体の心拍に由来する高調波周波数fH2を有するレーダ信号成分を、被検体の心拍に由来する基本周波数を有するレーダ信号成分として誤って抽出してしまう。よって、被検体の呼吸数を正確に計測することはできるが、被検体の心拍数を正確に計測することができない。 In the lower part of FIG. 1, the subject's breathing and heartbeat are faster than normal. Therefore, by applying a bandpass filter having a pass band corresponding to the general breathing rate to the radar signal, radar signal components having the fundamental frequency f B1 and harmonic frequencies f B2 and f B3 are extracted from the radar signal. However, the radar signal component with the fundamental frequency f H1 originating from the heartbeat cannot be removed from the radar signal. Then, by applying a bandpass filter having a pass band corresponding to a general heart rate to the radar signal, radar signal components having harmonic frequencies f H2 and f H3 can be extracted from the radar signal. The radar signal component with fundamental frequency f H1 cannot be extracted from the radar signal. That is, a radar signal component having a harmonic frequency f H2 derived from the subject's heartbeat is erroneously extracted as a radar signal component having a fundamental frequency originating from the subject's heartbeat. Therefore, although it is possible to accurately measure the respiratory rate of the subject, it is not possible to accurately measure the heart rate of the subject.
そこで、前記課題を解決するために、本開示は、レーダを用いて被検体の呼吸数及び心拍数を計測するにあたり、被検体の呼吸及び心拍に由来するレーダ信号成分として、基本周波数を有するレーダ信号成分のみならず、高調波周波数を有するレーダ信号成分が存在するときでも、被検体の呼吸数及び心拍数を正確に計測することを目的とする。 Therefore, in order to solve the above problems, the present disclosure provides a radar having a fundamental frequency as a radar signal component derived from the breathing and heartbeat of the subject when measuring the breathing rate and heart rate of the subject using a radar. It is an object of the present invention to accurately measure the respiratory rate and heart rate of a subject even when not only signal components but also radar signal components having harmonic frequencies are present.
前記課題を解決するために、レーダ信号の周波数変換結果の最大ピーク周波数を有するレーダ信号成分を、被検体の呼吸に由来する基本周波数を有するレーダ信号成分として抽出することとした。すると、被検体の呼吸に由来する高調波周波数を有するレーダ信号成分を特定することができる。そして、被検体の呼吸に由来する基本周波数及び高調波周波数を有するレーダ信号成分を、レーダ信号から除去することができる。さらに、残りのレーダ信号の周波数変換結果の最大ピーク周波数を有するレーダ信号成分を、被検体の心拍に由来する基本周波数を有するレーダ信号成分として抽出することができる。 In order to solve the above problem, it was decided to extract the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal as the radar signal component having the fundamental frequency derived from the breathing of the subject. Then, a radar signal component having a harmonic frequency originating from the subject's breathing can be identified. Then, radar signal components having a fundamental frequency and harmonic frequencies derived from the subject's breathing can be removed from the radar signal. Furthermore, the radar signal component having the maximum peak frequency of the frequency conversion results of the remaining radar signals can be extracted as the radar signal component having the fundamental frequency derived from the heartbeat of the subject.
具体的には、本開示は、被検体から反射されたレーダ信号を周波数変換し、前記レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の呼吸数を計測する呼吸数計測部と、前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去部と、前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測部と、を備えることを特徴とする呼吸心拍計測装置である。 Specifically, the present disclosure provides a respiration rate measuring unit that frequency-converts a radar signal reflected from a subject and measures the respiratory rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal. and removing a radar signal component originating from the subject's respiration from the radar signal by removing from the radar signal a radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency. After that, a respiratory component removal unit generates a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject, and a respiratory component removal unit that converts the frequency of the heartbeat radar signal and detects the maximum peak of the frequency conversion result of the heartbeat radar signal. The respiratory heart rate measuring device is characterized by comprising: a heart rate measuring section that measures the heart rate of the subject based on frequency.
また、本開示は、被検体から反射されたレーダ信号を周波数変換し、前記レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の呼吸数を計測する呼吸数計測ステップと、前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去ステップと、前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測ステップと、を順にコンピュータに実行させるための呼吸心拍計測プログラムである。 The present disclosure also provides a respiration rate measuring step of converting the frequency of a radar signal reflected from the subject and measuring the respiration rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal; By removing the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency from the radar signal, the radar signal component originating from the respiration of the subject is removed from the radar signal. , a step of removing a respiratory component to generate a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject, and converting the frequency of the heartbeat radar signal based on the maximum peak frequency of the frequency conversion result of the heartbeat radar signal. and a heart rate measuring step of measuring the heart rate of the subject.
これらの構成によれば、被検体の呼吸及び心拍に由来するレーダ信号成分として、基本周波数を有するレーダ信号成分のみならず、高調波周波数を有するレーダ信号成分が存在するときでも、被検体の呼吸数及び心拍数を正確に計測することができる。 According to these configurations, even when there are not only radar signal components having a fundamental frequency but also radar signal components having harmonic frequencies as radar signal components derived from the subject's breathing and heartbeat, the subject's breathing and heartbeat can be detected. It is possible to accurately measure the number of heartbeats and heart rate.
また、本開示は、前記呼吸成分除去部は、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分の振幅及び位相を算出するときに、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分を前記レーダ信号から除去した誤差値の絶対値が最小値となるようにすることを特徴とする呼吸心拍計測装置である。 Further, in the present disclosure, when calculating the amplitude and phase of a radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal, the respiratory component removal unit converts the frequency of the radar signal. The respiratory heart rate measuring device is characterized in that the absolute value of the error value obtained by removing the resulting radar signal component having the maximum peak frequency and the harmonic frequency from the radar signal becomes a minimum value.
この構成によれば、被検体の呼吸に由来する基本周波数及び高調波周波数でのバンド除去フィルタをレーダ信号に適用することなく、被検体の呼吸に由来する基本周波数及び高調波周波数を有するレーダ信号成分を、容易にレーダ信号から除去することができる。 According to this configuration, the radar signal having the fundamental frequency and harmonic frequencies originating from the subject's breathing can be transmitted without applying a band removal filter at the fundamental frequency and harmonic frequencies originating from the subject's breathing to the radar signal. components can be easily removed from the radar signal.
また、本開示は、前記呼吸成分除去部は、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分の振幅及び位相を算出するときに、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数での振幅及び位相を参照することを特徴とする呼吸心拍計測装置である。 Further, in the present disclosure, when calculating the amplitude and phase of a radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal, the respiratory component removal unit converts the frequency of the radar signal. The respiratory heart rate measuring device is characterized in that it refers to the maximum peak frequency of the result and the amplitude and phase at the harmonic frequency.
この構成によっても、被検体の呼吸に由来する基本周波数及び高調波周波数でのバンド除去フィルタをレーダ信号に適用することなく、被検体の呼吸に由来する基本周波数及び高調波周波数を有するレーダ信号成分を、容易にレーダ信号から除去することができる。 With this configuration as well, the radar signal component having the fundamental frequency and harmonic frequencies originating from the subject's breathing can be removed without applying a band removal filter at the fundamental frequency and harmonic frequencies originating from the subject's breathing to the radar signal. can be easily removed from the radar signal.
また、本開示は、前記呼吸成分除去部は、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数でのバンド除去フィルタを前記レーダ信号に適用するときに、前記バンド除去フィルタの周波数特性に起因する線形歪み信号成分を補償することを特徴とする呼吸心拍計測装置である。 Further, in the present disclosure, when applying a band removal filter at the maximum peak frequency of the frequency conversion result of the radar signal and the harmonic frequency to the radar signal, This is a respiratory heart rate measuring device characterized by compensating for linear distortion signal components caused by characteristics.
この構成によれば、被検体の呼吸に由来する基本周波数及び高調波周波数でのバンド除去フィルタをレーダ信号に適用するときでも、被検体の呼吸に由来する基本周波数及び高調波周波数を有するレーダ信号成分を、確実にレーダ信号から除去することができる。 According to this configuration, even when applying a band removal filter at the fundamental frequency and harmonic frequencies originating from the breathing of the subject to the radar signal, the radar signal having the fundamental frequency and harmonic frequencies originating from the breathing of the subject components can be reliably removed from the radar signal.
また、本開示は、前記心拍数計測部は、前記心拍レーダ信号と、前記心拍レーダ信号を擬似的に模した基準心拍レーダ信号と、の間の相互相関結果を周波数変換し、前記相互相関結果の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測することを特徴とする呼吸心拍計測装置である。 Further, in the present disclosure, the heart rate measuring unit frequency-converts a cross-correlation result between the heartbeat radar signal and a reference heartbeat radar signal that pseudo-simulates the heartbeat radar signal, and The respiratory heart rate measuring device is characterized in that the heart rate of the subject is measured based on the maximum peak frequency of the frequency conversion result.
この構成によれば、呼吸心拍計測装置に由来する高い周波数を有するノイズ成分や、被検体の呼吸に由来する低い周波数を有する信号成分が、心拍レーダ信号に重畳又は残留しているときでも、被検体の心拍数を正確に計測することができる。 According to this configuration, even when a noise component having a high frequency originating from the respiratory heart rate measuring device or a signal component having a low frequency originating from the subject's breathing is superimposed on or remains on the heart rate radar signal, the The heart rate of the sample can be measured accurately.
また、本開示は、前記呼吸数計測部は、前記レーダ信号の周波数変換結果の最大ピーク周波数を有するレーダ信号成分が前記レーダ信号の周波数変換結果の全体中に占める重み付けの大小に基づいて、前記被検体の呼吸数の計測結果の信頼度の高低を評価し、前記心拍数計測部は、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数を有するレーダ信号成分が前記心拍レーダ信号の周波数変換結果の全体中に占める重み付けの大小に基づいて、前記被検体の心拍数の計測結果の信頼度の高低を評価することを特徴とする呼吸心拍計測装置である。 Further, in the present disclosure, the respiration rate measurement unit calculates the amount of weight that a radar signal component having a maximum peak frequency of the frequency conversion result of the radar signal occupies in the entire frequency conversion result of the radar signal. The heart rate measurement unit evaluates the degree of reliability of the measurement result of the respiratory rate of the subject, and determines whether the radar signal component having the maximum peak frequency of the frequency conversion result of the heart rate radar signal is the result of the frequency conversion of the heart rate radar signal. The respiration heart rate measuring device is characterized in that the degree of reliability of the measurement result of the heart rate of the subject is evaluated based on the magnitude of weighting in the whole.
この構成によれば、被検体の呼吸数及び心拍数の計測結果の信頼度を評価することができる。そして、計測結果の信頼度が閾値と比べて高い/低いときには、計測結果を採用/棄却することができる。さらに、計測結果の信頼度に対する閾値が高い/低いときには、正確性/リアルタイム性を重視して計測結果を出力することができる。 According to this configuration, the reliability of the measurement results of the breathing rate and heart rate of the subject can be evaluated. Then, when the reliability of the measurement result is higher/lower than the threshold value, the measurement result can be adopted/discarded. Furthermore, when the threshold value for the reliability of the measurement result is high/low, the measurement result can be output with emphasis on accuracy/real-time performance.
このように、本開示は、レーダを用いて被検体の呼吸数及び心拍数を計測するにあたり、被検体の呼吸及び心拍に由来するレーダ信号成分として、基本周波数を有するレーダ信号成分のみならず、高調波周波数を有するレーダ信号成分が存在するときでも、被検体の呼吸数及び心拍数を正確に計測することができる。 As described above, the present disclosure provides that when measuring the respiration rate and heart rate of a subject using radar, not only the radar signal component having a fundamental frequency but also the radar signal component derived from the respiration and heartbeat of the subject are used. Even when radar signal components having harmonic frequencies are present, the respiratory rate and heart rate of the subject can be accurately measured.
添付の図面を参照して本開示の実施形態を説明する。以下に説明する実施形態は本開示の実施の例であり、本開示は以下の実施形態に制限されるものではない。 Embodiments of the present disclosure will be described with reference to the accompanying drawings. The embodiments described below are examples of implementation of the present disclosure, and the present disclosure is not limited to the following embodiments.
本開示の呼吸心拍計測システムを図2に示す。本開示の呼吸心拍計測手順を図3に示す。呼吸心拍計測システムSは、レーダ送受信装置1、呼吸心拍計測装置2及び呼吸心拍表示装置3を備え、レーダ送受信装置1を用いて、被検体Tに照射されるレーダ信号を送信したうえで、被検体Tから反射されたレーダ信号を受信するとともに、呼吸心拍表示装置3を用いて、被検体Tの呼吸数及び心拍数を表示することができる。呼吸心拍計測装置2は、呼吸数計測部21、呼吸成分除去部22及び心拍数計測部23を備え、図3に示した呼吸心拍計測プログラムをコンピュータにインストールして実現することができる。
FIG. 2 shows a respiratory heart rate measurement system of the present disclosure. FIG. 3 shows the respiratory heart rate measurement procedure of the present disclosure. The respiratory heart rate measuring system S includes a radar transmitting/
本開示の第1の呼吸計測処理を図4に示す。図4では、In-phase成分及びQuadrature-phase成分を重畳する。呼吸数計測部21は、被検体Tから反射されたレーダ信号S(t)を周波数変換し、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1に基づいて、被検体Tの呼吸数=60fB1/minを計測する(ステップS1)。つまり、呼吸数計測部21は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1を有するレーダ信号成分を、被検体Tの呼吸に由来する基本周波数fB1を有するレーダ信号成分として抽出する。そして、呼吸数計測部21は、被検体Tの呼吸に由来する高調波周波数fB2、fB3(更なる高調波周波数については省略している。)を有するレーダ信号成分を特定する。
FIG. 4 shows the first respiration measurement process of the present disclosure. In FIG. 4, the In-phase component and the Quadrature-phase component are superimposed. The respiration
よって、被検体Tの呼吸及び心拍に由来するレーダ信号成分として、基本周波数fB1、fH1を有するレーダ信号成分のみならず、高調波周波数fB2、fB3、fH2、fH3を有するレーダ信号成分が存在するときでも、被検体Tの呼吸数=60fB1/minを正確に計測することができる。なお、被検体Tの呼吸に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を考慮してもよい。 Therefore, radar signal components derived from the breathing and heartbeat of the subject T include not only radar signal components having fundamental frequencies f B1 and f H1 but also radar signal components having harmonic frequencies f B2 , f B3 , f H2 , and f H3 . Even when a signal component exists, the respiratory rate of the subject T = 60 f B1 /min can be accurately measured. Note that as the radar signal component derived from the breathing of the subject T, a radar signal component having a subharmonic frequency such as half of the fundamental frequency may be considered.
呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及びその高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号S(t)から除去する(ステップS2)。つまり、呼吸成分除去部22は、被検体Tの呼吸に由来するレーダ信号成分をレーダ信号S(t)から除去したうえで、被検体Tの心拍に由来するレーダ信号成分を抽出した心拍レーダ信号SH(t)を生成する(ステップS3)。
The respiratory
具体的には、呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分の振幅及び位相を算出するときに、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号S(t)から除去した誤差値の絶対値が最小値となるようにする。すると、心拍レーダ信号SH(t)は、数1で表される。
よって、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3でのバンド除去フィルタをレーダ信号S(t)に適用することなく、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分を、容易にレーダ信号S(t)から除去することができる。なお、被検体Tの呼吸に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を除去してもよい。 Therefore, without applying a band removal filter at the fundamental frequency f B1 and harmonic frequencies f B2 and f B3 originating from the breathing of the subject T to the radar signal S(t), the fundamental frequency originating from the breathing of the subject T can be removed. Radar signal components having frequency f B1 and harmonic frequencies f B2 , f B3 can be easily removed from radar signal S(t). Note that as the radar signal component derived from the breathing of the subject T, a radar signal component having a subharmonic frequency such as half of the fundamental frequency may be removed.
本開示の第2の呼吸計測処理を図5に示す。図5では、In-phase成分及びQuadrature-phase成分を重畳する。ステップS1については、図5に示した第2の呼吸計測処理は、図4に示した第1の呼吸計測処理と同様である。 FIG. 5 shows the second respiration measurement process of the present disclosure. In FIG. 5, the In-phase component and the Quadrature-phase component are superimposed. Regarding step S1, the second respiration measurement process shown in FIG. 5 is similar to the first respiration measurement process shown in FIG. 4.
呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及びその高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号S(t)から除去する(ステップS2)。つまり、呼吸成分除去部22は、被検体Tの呼吸に由来するレーダ信号成分をレーダ信号S(t)から除去したうえで、被検体Tの心拍に由来するレーダ信号成分を抽出した心拍レーダ信号SH(t)を生成する(ステップS3)。
The respiratory
具体的には、呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分の振幅及び位相を算出するときに、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3での振幅及び位相を参照する。なお、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3について、振幅は図示されているが位相は図示されていない。心拍レーダ信号SH(t)は、数2で表される。
よって、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3でのバンド除去フィルタをレーダ信号S(t)に適用することなく、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分を、容易にレーダ信号S(t)から除去することができる。なお、被検体Tの呼吸に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を除去してもよい。 Therefore, without applying a band removal filter at the fundamental frequency f B1 and harmonic frequencies f B2 and f B3 originating from the breathing of the subject T to the radar signal S(t), the fundamental frequency originating from the breathing of the subject T can be removed. Radar signal components having frequency f B1 and harmonic frequencies f B2 , f B3 can be easily removed from radar signal S(t). Note that as the radar signal component derived from the breathing of the subject T, a radar signal component having a subharmonic frequency such as half of the fundamental frequency may be removed.
本開示の第3の呼吸計測処理を図6、7に示す。図6、7では、In-phase成分及びQuadrature-phase成分を重畳する。ステップS1については、図6、7に示した第3の呼吸計測処理は、図4に示した第1の呼吸計測処理と同様である。 The third respiration measurement process of the present disclosure is shown in FIGS. 6 and 7. In FIGS. 6 and 7, the In-phase component and the Quadrature-phase component are superimposed. Regarding step S1, the third respiration measurement process shown in FIGS. 6 and 7 is the same as the first respiration measurement process shown in FIG. 4.
呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及びその高調波周波数fB2、fB3を有するレーダ信号成分をレーダ信号S(t)から除去する(ステップS2)。つまり、呼吸成分除去部22は、被検体Tの呼吸に由来するレーダ信号成分をレーダ信号S(t)から除去したうえで、被検体Tの心拍に由来するレーダ信号成分を抽出した心拍レーダ信号SH(t)を生成する(ステップS3)。
The respiratory
具体的には、呼吸成分除去部22は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3でのバンド除去フィルタをレーダ信号S(t)に適用するときに、バンド除去フィルタの周波数特性に起因する線形歪み信号成分を補償する。なお、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1及び高調波周波数fB2、fB3は特定されており、バンド除去フィルタの周波数特性は既知であるため、バンド除去フィルタの周波数特性に起因する線形歪み信号成分を計算することができる。
Specifically, the respiratory
図6の第3段では、一括のバンド除去フィルタを適用しており、最大ピーク周波数fB1から高調波周波数fB3までの周波数を有するレーダ信号成分を一括に除去している。図7の第3段では、くし型のバンド除去フィルタを適用しており、最大ピーク周波数fB1、高調波周波数fB2及び高調波周波数fB3を有するレーダ信号成分を個別に除去している。 In the third stage of FIG. 6, a collective band removal filter is applied, and radar signal components having frequencies from the maximum peak frequency f B1 to the harmonic frequency f B3 are collectively removed. In the third stage of FIG. 7, a comb-shaped band rejection filter is applied, and radar signal components having maximum peak frequency f B1 , harmonic frequency f B2 and harmonic frequency f B3 are individually removed.
よって、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3でのバンド除去フィルタをレーダ信号S(t)に適用するときでも、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2、fB3を有するレーダ信号成分を、確実にレーダ信号S(t)から除去することができる。なお、被検体Tの呼吸に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を除去してもよい。 Therefore, even when applying a band removal filter at the fundamental frequency f B1 and harmonic frequencies f B2 and f B3 derived from the breathing of the subject T to the radar signal S(t), the fundamental frequency f B1 originating from the breathing of the subject T is applied to the radar signal S(t). Radar signal components having frequency f B1 and harmonic frequencies f B2 , f B3 can be reliably removed from radar signal S(t). Note that as the radar signal component derived from the breathing of the subject T, a radar signal component having a subharmonic frequency such as half of the fundamental frequency may be removed.
本開示の第1の心拍計測処理を図8に示す。図8では、In-phase成分及びQuadrature-phase成分を重畳する。心拍数計測部23は、心拍レーダ信号SH(t)を周波数変換し、心拍レーダ信号SH(t)の周波数変換結果の最大ピーク周波数fH1に基づいて、被検体Tの心拍数=60fH1/minを計測する(ステップS4)。つまり、心拍数計測部23は、心拍レーダ信号SH(t)の周波数変換結果の最大ピーク周波数fH1を有するレーダ信号成分を、被検体Tの心拍に由来する基本周波数fH1を有するレーダ信号成分として抽出する。そして、心拍数計測部23は、被検体Tの心拍に由来する高調波周波数fH2、fH3(更なる高調波周波数については省略している。)を有するレーダ信号成分を特定する。
FIG. 8 shows the first heart rate measurement process of the present disclosure. In FIG. 8, the In-phase component and the Quadrature-phase component are superimposed. The heart
ここで、心拍レーダ信号SH(t)の周波数変換結果では、被検体Tの呼吸に由来する基本周波数fB1及び高調波周波数fB2を有するレーダ信号成分が、若干残留しているが被検体Tの心拍に由来する基本周波数fH1を有するレーダ信号成分と比べて小さい。 Here, in the frequency conversion result of the heartbeat radar signal S H (t), some radar signal components having the fundamental frequency f B1 and harmonic frequency f B2 derived from the breathing of the subject T remain, but the is small compared to the radar signal component with fundamental frequency f H1 originating from the heartbeat of T.
よって、被検体Tの呼吸及び心拍に由来するレーダ信号成分として、基本周波数fB1、fH1を有するレーダ信号成分のみならず、高調波周波数fB2、fB3、fH2、fH3を有するレーダ信号成分が存在するときでも、被検体Tの心拍数=60fH1/minを正確に計測することができる。なお、被検体Tの心拍に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を考慮してもよい。 Therefore, radar signal components derived from the breathing and heartbeat of the subject T include not only radar signal components having fundamental frequencies f B1 and f H1 but also radar signal components having harmonic frequencies f B2 , f B3 , f H2 , and f H3 . Even when a signal component exists, the heart rate of the subject T = 60 f H1 /min can be accurately measured. Note that as the radar signal component derived from the heartbeat of the subject T, a radar signal component having a subharmonic frequency such as half of the fundamental frequency may be considered.
本開示の第2の心拍計測処理を図9に示す。図9では、In-phase成分及びQuadrature-phase成分を重畳する。心拍数計測部23は、心拍レーダ信号SH(t)と、心拍レーダ信号SH(t)を擬似的に模した基準心拍レーダ信号SS(t)と、の間の相互相関結果SC(t)を周波数変換し、相互相関結果SC(t)の周波数変換結果の最大ピーク周波数fH1に基づいて、被検体Tの心拍数=60fH1/minを計測する(ステップS4)。つまり、心拍数計測部23は、相互相関結果SC(t)の周波数変換結果の最大ピーク周波数fH1を有するレーダ信号成分を、被検体Tの心拍に由来する基本周波数fH1を有するレーダ信号成分として抽出する。そして、心拍数計測部23は、被検体Tの心拍に由来する高調波周波数fH2、fH3(更なる高調波周波数については省略している。)を有するレーダ信号成分を特定する。
FIG. 9 shows the second heart rate measurement process of the present disclosure. In FIG. 9, the In-phase component and the Quadrature-phase component are superimposed. The heart
ここで、基準心拍レーダ信号SS(t)として、cos基底である基準心拍レーダ信号SCOS(t)と、sin基底である基準心拍レーダ信号SSIN(t)と、を準備する。そして、相互相関結果SC(t)として、数3に示す相互相関を計算する。
よって、呼吸心拍計測装置2に由来する高い周波数を有するノイズ成分や、被検体Tの呼吸に由来する低い周波数fB1、fB2、fB3を有する信号成分が、心拍レーダ信号SH(t)に重畳又は残留しているときでも、被検体Tの心拍数=60fH1/minを正確に計測することができる。なお、被検体Tの心拍に由来するレーダ信号成分として、基本周波数の半分等の低調波周波数を有するレーダ信号成分を考慮してもよい。
Therefore, the noise component having a high frequency originating from the respiratory heart
本開示の呼吸数評価処理を図10に示す。図10では、In-phase成分及びQuadrature-phase成分を重畳する。呼吸数計測部21は、レーダ信号S(t)の周波数変換結果の最大ピーク周波数fB1を有するレーダ信号成分がレーダ信号S(t)の周波数変換結果の全体中に占める重み付けの大小に基づいて、被検体Tの呼吸数=60fB1/minの計測結果の信頼度の高低を評価する(ステップS5)。
FIG. 10 shows the respiration rate evaluation process of the present disclosure. In FIG. 10, the In-phase component and the Quadrature-phase component are superimposed. The respiration
具体的には、呼吸数計測部21は、被検体Tの呼吸数の計測結果の信頼度として、(最大ピーク周波数fB1の信号成分)/(全ピーク周波数又は上位ピーク周波数の信号成分和)を計算する。なお、被検体Tの呼吸数の計測結果の信頼度として、(ピーク周波数fB1、fB2、fB3の信号成分和)/(全ピーク周波数又は上位ピーク周波数の信号成分和)を用いてもよい。そして、呼吸数計測部21は、被検体Tの呼吸数の計測結果の信頼度が閾値と比べて高い/低いときには、被検体Tの呼吸数の計測結果を採用/棄却する。さらに、呼吸数計測部21は、被検体Tの呼吸数の計測結果に対する閾値が高い/低いときには、正確性/リアルタイム性を重視して被検体Tの呼吸数の計測結果を出力する。
Specifically, the respiration
本開示の心拍数評価処理を図11に示す。図11では、In-phase成分及びQuadrature-phase成分を重畳する。心拍数計測部23は、心拍レーダ信号SH(t)の周波数変換結果の最大ピーク周波数fH1を有するレーダ信号成分が心拍レーダ信号SH(t)の周波数変換結果の全体中に占める重み付けの大小に基づいて、被検体Tの心拍数=60fH1/minの計測結果の信頼度の高低を評価する(ステップS6)。
FIG. 11 shows the heart rate evaluation process of the present disclosure. In FIG. 11, the In-phase component and the Quadrature-phase component are superimposed. The heart
具体的には、心拍数計測部23は、被検体Tの心拍数の計測結果の信頼度として、(最大ピーク周波数fH1の信号成分)/(全ピーク周波数又は上位ピーク周波数の信号成分和)を計算する。なお、被検体Tの心拍数の計測結果の信頼度として、(ピーク周波数fH1、fH2、fH3の信号成分和)/(全ピーク周波数又は上位ピーク周波数の信号成分和)を用いてもよい。そして、心拍数計測部23は、被検体Tの心拍数の計測結果の信頼度が閾値と比べて高い/低いときには、被検体Tの心拍数の計測結果を採用/棄却する。さらに、心拍数計測部23は、被検体Tの心拍数の計測結果に対する閾値が高い/低いときには、正確性/リアルタイム性を重視して被検体Tの心拍数の計測結果を出力する。
Specifically, the heart
本開示の呼吸計測結果を図12に示す。図12の上部の縦線は、被検体Tの呼吸数の計測結果の信頼度が閾値と比べて低いため、被検体Tの呼吸数の計測結果が棄却されていることを示す。図12の下部の白丸は、被検体Tの呼吸数の計測結果の信頼度が閾値と比べて高いため、被検体Tの呼吸数の計測結果が採用されていることを示す。なお、被検体Tの呼吸数の計測結果の信頼度が閾値と比べて高いときでも、被検体Tの呼吸数の計測結果が急激に変化したときには、被検体Tの呼吸数の計測結果を棄却してもよい。 The respiration measurement results of the present disclosure are shown in FIG. 12. The vertical line at the top of FIG. 12 indicates that the measurement result of the respiratory rate of the subject T is rejected because the reliability of the measurement result of the respiratory rate of the subject T is lower than the threshold value. The white circle at the bottom of FIG. 12 indicates that the measurement result of the respiratory rate of the subject T is adopted because the reliability of the measurement result of the respiratory rate of the subject T is higher than the threshold value. Furthermore, even when the reliability of the measurement result of the respiration rate of the subject T is high compared to the threshold value, if the measurement result of the respiration rate of the subject T changes suddenly, the measurement result of the respiration rate of the subject T may be rejected. You may.
本開示の呼吸心拍計測装置及び呼吸心拍計測プログラムは、レーダを用いて被検体の呼吸数及び心拍数を計測するにあたり、被検体の呼吸及び心拍に由来するレーダ信号成分として、基本周波数を有するレーダ信号成分のみならず、高調波周波数を有するレーダ信号成分が存在するときでも、被検体の呼吸数及び心拍数を正確に計測することができる。 The respiration heart rate measurement device and the respiration heart rate measurement program of the present disclosure use a radar having a fundamental frequency as a radar signal component derived from the respiration and heartbeat of the object when measuring the respiration rate and heart rate of the object using radar. Even when not only signal components but also radar signal components having harmonic frequencies are present, the respiratory rate and heart rate of the subject can be accurately measured.
S:呼吸心拍計測システム
T:被検体
1:レーダ送受信装置
2:呼吸心拍計測装置
3:呼吸心拍表示装置
21:呼吸数計測部
22:呼吸成分除去部
23:心拍数計測部
S: Respiratory heart rate measuring system T: Subject 1: Radar transmitting/receiving device 2: Respiratory heart rate measuring device 3: Respiratory heart rate display device 21: Respiratory rate measuring section 22: Respiratory component removing section 23: Heart rate measuring section
Claims (8)
前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去部と、
前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測部と、を備え、
前記呼吸成分除去部は、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分の振幅及び位相を算出するときに、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分を前記レーダ信号から除去した誤差値の絶対値が最小値となるようにする
ことを特徴とする呼吸心拍計測装置。 a respiration rate measuring unit that converts the frequency of a radar signal reflected from the subject and measures the respiratory rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal;
By removing from the radar signal the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency, the radar signal component originating from the respiration of the subject is removed from the radar signal. a respiratory component removal unit that generates a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject;
a heart rate measuring unit that converts the frequency of the heart rate radar signal and measures the heart rate of the subject based on the maximum peak frequency of the frequency conversion result of the heart rate radar signal ,
The respiratory component removal unit calculates the maximum peak frequency and phase of the radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal. The absolute value of the error value obtained by removing the radar signal component having the harmonic frequency from the radar signal is set to a minimum value.
A respiratory heart rate measuring device characterized by:
前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去部と、
前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測部と、を備え、
前記心拍数計測部は、前記心拍レーダ信号と、前記心拍レーダ信号を擬似的に模したマザーウェーブレット関数を使った基準心拍レーダ信号と、の間の相互相関結果を周波数変換し、前記相互相関結果の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する
ことを特徴とする呼吸心拍計測装置。 a respiration rate measuring unit that converts the frequency of a radar signal reflected from the subject and measures the respiratory rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal;
By removing from the radar signal the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency, the radar signal component originating from the respiration of the subject is removed from the radar signal. a respiratory component removing unit that generates a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject;
a heart rate measuring unit that converts the frequency of the heart rate radar signal and measures the heart rate of the subject based on the maximum peak frequency of the frequency conversion result of the heart rate radar signal ,
The heart rate measuring unit frequency-converts the cross-correlation result between the heartbeat radar signal and a reference heartbeat radar signal using a mother wavelet function that simulates the heartbeat radar signal, and converts the cross-correlation result. Measure the heart rate of the subject based on the maximum peak frequency of the frequency conversion result.
A respiratory heart rate measuring device characterized by:
ことを特徴とする、請求項2に記載の呼吸心拍計測装置。 The respiratory component removal unit calculates the maximum peak frequency and phase of the frequency conversion result of the radar signal when calculating the amplitude and phase of the radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal. The respiratory heart rate measuring device according to claim 2 , wherein the absolute value of an error value obtained by removing the radar signal component having the harmonic frequency from the radar signal is a minimum value.
ことを特徴とする、請求項2に記載の呼吸心拍計測装置。 The respiratory component removal unit calculates the maximum peak frequency and phase of the radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal. The respiratory heart rate measuring device according to claim 2 , characterized in that the amplitude and phase at the harmonic frequency are referred to.
ことを特徴とする、請求項2に記載の呼吸心拍計測装置。 The respiratory component removing unit removes linear distortion due to frequency characteristics of the band removal filter when applying the band removal filter at the maximum peak frequency of the frequency conversion result of the radar signal and the harmonic frequency to the radar signal. The respiratory heart rate measuring device according to claim 2 , wherein the respiratory heart rate measuring device compensates for signal components.
前記心拍数計測部は、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数を有するレーダ信号成分が前記心拍レーダ信号の周波数変換結果の全体中に占める重み付けの大小に基づいて、前記被検体の心拍数の計測結果の信頼度の高低を評価する
ことを特徴とする、請求項1から5のいずれかに記載の呼吸心拍計測装置。 The respiration rate measurement unit measures the respiration rate of the subject based on the weighting that a radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal occupies in the entire frequency conversion result of the radar signal. Evaluate the reliability of measurement results,
The heart rate measuring unit calculates the heart rate of the subject based on the weight that a radar signal component having the maximum peak frequency of the frequency conversion result of the heart rate radar signal occupies in the entire frequency conversion result of the heart rate radar signal. The respiratory heart rate measuring device according to any one of claims 1 to 5, characterized in that the reliability level of the measurement result of the number of heartbeats is evaluated.
前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去ステップと、
前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測ステップと、
を順にコンピュータに実行させ、
前記呼吸成分除去ステップは、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分の振幅及び位相を算出するときに、前記レーダ信号の周波数変換結果の最大ピーク周波数及び前記高調波周波数を有するレーダ信号成分を前記レーダ信号から除去した誤差値の絶対値が最小値となるようにする
ことを特徴とする呼吸心拍計測プログラム。 a respiration rate measuring step of converting the frequency of a radar signal reflected from the subject and measuring the respiration rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal;
By removing from the radar signal the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency, the radar signal component originating from the respiration of the subject is removed from the radar signal. a respiratory component removal step of generating a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject;
a heart rate measuring step of converting the frequency of the heart rate radar signal and measuring the heart rate of the subject based on the maximum peak frequency of the frequency conversion result of the heart rate radar signal;
have the computer execute the following in order ,
The respiratory component removing step includes calculating the maximum peak frequency and phase of the radar signal component having the maximum peak frequency and the harmonic frequency of the frequency conversion result of the radar signal. The absolute value of the error value obtained by removing the radar signal component having the harmonic frequency from the radar signal is set to a minimum value.
A respiratory heart rate measurement program.
前記レーダ信号の周波数変換結果の最大ピーク周波数及びその高調波周波数を有するレーダ信号成分を前記レーダ信号から除去することにより、前記被検体の呼吸に由来するレーダ信号成分を前記レーダ信号から除去したうえで、前記被検体の心拍に由来するレーダ信号成分を抽出した心拍レーダ信号を生成する呼吸成分除去ステップと、
前記心拍レーダ信号を周波数変換し、前記心拍レーダ信号の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する心拍数計測ステップと、
を順にコンピュータに実行させ、
前記心拍数計測ステップは、前記心拍レーダ信号と、前記心拍レーダ信号を擬似的に模したマザーウェーブレット関数を使った基準心拍レーダ信号と、の間の相互相関結果を周波数変換し、前記相互相関結果の周波数変換結果の最大ピーク周波数に基づいて、前記被検体の心拍数を計測する
ことを特徴とする呼吸心拍計測プログラム。 a respiration rate measuring step of converting the frequency of a radar signal reflected from the subject and measuring the respiration rate of the subject based on the maximum peak frequency of the frequency conversion result of the radar signal;
By removing from the radar signal the radar signal component having the maximum peak frequency of the frequency conversion result of the radar signal and its harmonic frequency, the radar signal component originating from the respiration of the subject is removed from the radar signal. a respiratory component removal step of generating a heartbeat radar signal by extracting a radar signal component derived from the heartbeat of the subject;
a heart rate measuring step of converting the frequency of the heart rate radar signal and measuring the heart rate of the subject based on the maximum peak frequency of the frequency conversion result of the heart rate radar signal;
have the computer execute the following in order ,
The heart rate measuring step frequency-converts the cross-correlation result between the heartbeat radar signal and a reference heartbeat radar signal using a mother wavelet function that simulates the heartbeat radar signal, and converts the cross-correlation result. Measure the heart rate of the subject based on the maximum peak frequency of the frequency conversion result.
A respiratory heart rate measurement program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020114641A JP7428605B2 (en) | 2020-07-02 | 2020-07-02 | Respiratory heart rate measuring device and respiratory heart rate measuring program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020114641A JP7428605B2 (en) | 2020-07-02 | 2020-07-02 | Respiratory heart rate measuring device and respiratory heart rate measuring program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022012646A JP2022012646A (en) | 2022-01-17 |
JP7428605B2 true JP7428605B2 (en) | 2024-02-06 |
Family
ID=80148849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020114641A Active JP7428605B2 (en) | 2020-07-02 | 2020-07-02 | Respiratory heart rate measuring device and respiratory heart rate measuring program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7428605B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114983373B (en) * | 2022-06-02 | 2023-03-28 | 谢俊 | Method for detecting human heart rate |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009055997A (en) | 2007-08-30 | 2009-03-19 | Honda Motor Co Ltd | Biological vibration frequency detection device and vehicle |
JP2015521075A (en) | 2012-05-28 | 2015-07-27 | オービーエス メディカル リミテッド | Extraction of respiratory rate from cardiac signal |
JP2017169868A (en) | 2016-03-24 | 2017-09-28 | 新日本無線株式会社 | Cardiopulmonary function measuring device |
JP2018187129A (en) | 2017-05-09 | 2018-11-29 | シャープ株式会社 | Measurement device, measurement method and measurement program |
JP2019129996A (en) | 2018-01-31 | 2019-08-08 | 学校法人慶應義塾 | Heartbeat detection system and heartbeat detection method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02233032A (en) * | 1989-03-07 | 1990-09-14 | Kokusai Electric Co Ltd | Linear distortion compensation band stop filter |
-
2020
- 2020-07-02 JP JP2020114641A patent/JP7428605B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009055997A (en) | 2007-08-30 | 2009-03-19 | Honda Motor Co Ltd | Biological vibration frequency detection device and vehicle |
JP2015521075A (en) | 2012-05-28 | 2015-07-27 | オービーエス メディカル リミテッド | Extraction of respiratory rate from cardiac signal |
JP2017169868A (en) | 2016-03-24 | 2017-09-28 | 新日本無線株式会社 | Cardiopulmonary function measuring device |
JP2018187129A (en) | 2017-05-09 | 2018-11-29 | シャープ株式会社 | Measurement device, measurement method and measurement program |
JP2019129996A (en) | 2018-01-31 | 2019-08-08 | 学校法人慶應義塾 | Heartbeat detection system and heartbeat detection method |
Also Published As
Publication number | Publication date |
---|---|
JP2022012646A (en) | 2022-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Varghees et al. | A novel heart sound activity detection framework for automated heart sound analysis | |
US8655436B2 (en) | Heart rate meter and heart beat detecting method | |
CN104820786B (en) | A kind of instantaneous weighting is synchronous to squeeze small echo double-spectrum analysis method | |
JP7428605B2 (en) | Respiratory heart rate measuring device and respiratory heart rate measuring program | |
JPWO2020003910A1 (en) | Heart rate detector, heart rate detection method and program | |
JP2017012366A (en) | Heart disease diagnosis apparatus, heart disease diagnosis program and medium | |
Tan et al. | EMD-based electrocardiogram delineation for a wearable low-power ECG monitoring device | |
JP5975880B2 (en) | Speech recognition using multiple parallel complex filters for fast extraction of formants | |
CN108294736A (en) | Continuous BP measurement system and measurement method | |
US10980485B2 (en) | Measuring apparatus, measuring method and non-transitory computer readable medium | |
TWI504378B (en) | Denoising method and apparatus of pulse wave signal and pulse oximetry | |
TW201332512A (en) | Method and apparatus for heart rate measurement | |
US11911136B2 (en) | System and method for calculating cardiac pulse transit or arrival time information | |
US11721358B2 (en) | System and method for calculating cardiovascular heartbeat information from an electronic audio signal | |
Li et al. | Robust suppression of nonstationary power-line interference in electrocardiogram signals | |
CN109682958B (en) | Acceleration sensor signal compensation method for thrombelastogram instrument | |
KR20180067348A (en) | Apparatus and method for estimating heart beat using photo-plethysmogram | |
Shervegar | Sliding Window Adaptive Filter for PCG denoising | |
Xia et al. | An Effective Photoplethysmography Denosing Method Based on Diffusion Probabilistic Model | |
Chukka et al. | Analysis and comparison of ECG signal quality assessments methods | |
CN119097344B (en) | An electrocardiogram measurement method based on ultrasonic signals from smart speakers | |
JP2849711B2 (en) | Accelerometer | |
JP2020092783A (en) | Cardiac cycle measuring device and cardiac cycle measuring program | |
Sanz et al. | A novel approach to detecting pulse onset in photoplethysmographic signal using an automatic non assisted method | |
Zhang et al. | Detection of saturation and artifact |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230703 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231121 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240117 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240123 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240125 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7428605 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |