JP2015119770A - Biological information measurement device and biological information measurement method - Google Patents
Biological information measurement device and biological information measurement method Download PDFInfo
- Publication number
- JP2015119770A JP2015119770A JP2013264172A JP2013264172A JP2015119770A JP 2015119770 A JP2015119770 A JP 2015119770A JP 2013264172 A JP2013264172 A JP 2013264172A JP 2013264172 A JP2013264172 A JP 2013264172A JP 2015119770 A JP2015119770 A JP 2015119770A
- Authority
- JP
- Japan
- Prior art keywords
- matrix
- person
- present
- measuring apparatus
- life activity
- 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.)
- Granted
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 28
- 238000000691 measurement method Methods 0.000 title 1
- 239000011159 matrix material Substances 0.000 claims abstract description 104
- 230000000694 effects Effects 0.000 claims abstract description 73
- 230000004044 response Effects 0.000 claims abstract description 27
- 230000004071 biological effect Effects 0.000 claims abstract description 17
- 230000005540 biological transmission Effects 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 abstract description 9
- 230000029058 respiratory gaseous exchange Effects 0.000 description 27
- 238000007781 pre-processing Methods 0.000 description 15
- 230000008859 change Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000036387 respiratory rate Effects 0.000 description 3
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000014155 detection of activity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
本発明は、生体情報測定装置及び生体情報測定方法に関する。 The present invention relates to a biological information measuring device and a biological information measuring method.
人間や動物等の被検体の心拍数や呼吸数等の身体情報を測定する生体情報測定装置が従来から知られている。このような生体情報測定装置に関する技術として、例えば特許文献1で開示されているような技術が開示されている。 2. Description of the Related Art Conventionally, a biological information measuring device that measures physical information such as heart rate and respiratory rate of a subject such as a human or an animal is known. As a technique related to such a biological information measuring device, for example, a technique disclosed in Patent Document 1 is disclosed.
特許文献1では、ベッドに内蔵されたマイクロ波送受信センサから、ベッドに横たわる被験者の胸部に向けて電磁波としてのマイクロ波を照射し、反射電磁波、すなわち、被験者から反射したマイクロ波を受信し、受信したマイクロ波に基づいて、被験者の心拍数や呼吸数等の身体状態を検出するものが提案されている。 In Patent Document 1, a microwave as an electromagnetic wave is irradiated from a microwave transmission / reception sensor built in a bed toward a chest of a subject lying on the bed, and a reflected electromagnetic wave, that is, a microwave reflected from the subject is received and received. Based on such microwaves, there has been proposed one that detects a physical state such as a heart rate and a respiration rate of a subject.
このような電波を利用した検出手法としては、従前からドップラーレーダが知られている。しかしドップラーレーザを利用した場合には、人との距離が遠くなると推定精度が悪くなるという問題があった。 As a detection method using such radio waves, Doppler radar has been known. However, when the Doppler laser is used, there is a problem in that the estimation accuracy deteriorates as the distance to the person increases.
そこで、MIMO(Multiple Input Multiple Output)技術を利用して、生体活動を検出する方法も従前から提案されている。 Therefore, a method for detecting a biological activity using a MIMO (Multiple Input Multiple Output) technique has been proposed.
しかし、MIMOを利用した場合には、人がいない場合の伝搬チャネルの周波数特性が測定するチャネルや伝搬環境の違いにより変化することがあり、伝搬チャネルの各要素での値が大きく異なることもある。基準にするチャネルをどれにするかで、人が入ることにより呼吸や心拍を推定するためのしきい値が変化する。この値の設定により、人がいない場合でも呼吸や心拍が誤って推定される場合が生じるという問題があった。 However, when MIMO is used, the frequency characteristics of the propagation channel when there is no person may change due to differences in the channel to be measured and the propagation environment, and the values at each element of the propagation channel may vary greatly. . Depending on which channel is used as a reference, the threshold for estimating respiration and heartbeat changes when a person enters. Due to the setting of this value, there is a problem that even when there is no person, breathing and heartbeat may be erroneously estimated.
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、MIMO技術を利用して、より正確に人物の活動の状態を推定することが可能な、新規かつ改良された生体情報測定装置及び生体情報測定方法を提供することにある。 Therefore, the present invention has been made in view of the above problems, and an object of the present invention is a novel technique that can estimate a person's activity state more accurately using MIMO technology. Another object of the present invention is to provide an improved biological information measuring apparatus and biological information measuring method.
上記課題を解決するために、本発明のある観点によれば、複数の送信アンテナと複数の受信アンテナとの間で形成される伝搬路を示す第1の行列を、前記複数の送信アンテナから送信され、前記複数の受信アンテナで受信される信号に基づいて測定する行列測定部と、前記伝搬路の間に対象物が存在する場合に前記第1の行列の時間応答を周波数応答に変換して得られる第2の行列及び前記伝搬路の間に対象物が存在しない場合に前記第1の行列の時間応答を周波数応答に変換して得られる第3の行列を算出する行列算出部と、前記第2の行列の要素毎に該第2の行列の各要素の複素共役を掛け合わせて得られる第1の結果及び前記第3の行列に該第2の行列の各要素の複素共役を掛け合わせて得られる第2の結果を用いて所定の条件を満たすか否かを判定し、該所定の条件を満たしていれば前記対象物が活動していると判定する判定部と、を備える、生体活動測定装置が提供される。 In order to solve the above-described problem, according to an aspect of the present invention, a first matrix indicating a propagation path formed between a plurality of transmission antennas and a plurality of reception antennas is transmitted from the plurality of transmission antennas. A matrix measuring unit for measuring based on signals received by the plurality of receiving antennas, and converting a time response of the first matrix into a frequency response when an object exists between the propagation paths. A matrix calculator that calculates a third matrix obtained by converting a time response of the first matrix into a frequency response when there is no object between the obtained second matrix and the propagation path; A first result obtained by multiplying each element of the second matrix by a complex conjugate of each element of the second matrix and a complex conjugate of each element of the second matrix multiplied by the third matrix Satisfy the predetermined condition using the second result obtained Whether determined, and a determination section determines that the object should satisfy the predetermined condition is active, the biological activity measuring device is provided.
かかる構成によれば、行列推定部は、複数の送信アンテナと複数の受信アンテナとの間で形成される伝搬路を示す第1の行列を、複数の送信アンテナから送信され、複数の受信アンテナで受信される信号に基づいて測定する。また行列算出部は、伝搬路の間に対象物が存在する場合に第1の行列の時間応答を周波数応答に変換して得られる第2の行列及び伝搬路の間に対象物が存在しない場合に第1の行列の時間応答を周波数応答に変換して得られる第3の行列を算出する。そして、判定部は、第2の行列の要素毎に該第2の行列の各要素の複素共役を掛け合わせて得られる第1の結果及び第3の行列に該第2の行列の各要素の複素共役を掛け合わせて得られる第2の結果を用いて所定の条件を満たすか否かを判定し、該所定の条件を満たしていれば対象物が活動していると判定する。 According to this configuration, the matrix estimation unit transmits the first matrix indicating the propagation path formed between the plurality of transmission antennas and the plurality of reception antennas from the plurality of transmission antennas, and uses the plurality of reception antennas. Measure based on the received signal. Further, the matrix calculation unit, when there is no object between the second matrix and the propagation path obtained by converting the time response of the first matrix into the frequency response when the object exists between the propagation paths A third matrix obtained by converting the time response of the first matrix into a frequency response is calculated. Then, the determination unit multiplies the first result obtained by multiplying the complex conjugate of each element of the second matrix for each element of the second matrix and the third matrix for each element of the second matrix. The second result obtained by multiplying the complex conjugate is used to determine whether or not a predetermined condition is satisfied. If the predetermined condition is satisfied, it is determined that the object is active.
本発明のある観点による生体活動測定装置は、かかる構成を有することによって、複数の送信アンテナと複数の受信アンテナとの間で電波を送受信するMIMO技術を利用して、より正確に人物の活動の状態を推定することが可能となる。 The life activity measuring apparatus according to a certain aspect of the present invention has such a configuration, and thus uses a MIMO technology for transmitting and receiving radio waves between a plurality of transmitting antennas and a plurality of receiving antennas, thereby more accurately detecting the activity of a person. The state can be estimated.
前記所定の条件は、前記第2の結果の絶対値より前記第1の結果の方が全ての要素で大きいことであってもよい。 The predetermined condition may be that the first result is larger in all elements than the absolute value of the second result.
前記所定の条件は、前記第2の結果を全ての要素で加算した値の絶対値より前記第1の結果を全ての要素で加算した値の方が大きいことであってもよい。 The predetermined condition may be that a value obtained by adding the first result for all elements is larger than an absolute value of a value obtained by adding the second result for all elements.
また、上記課題を解決するために、本発明の別の観点によれば、複数の送信アンテナと複数の受信アンテナとの間で形成される伝搬路を示す第1の行列を、前記複数の送信アンテナから送信され、前記複数の受信アンテナで受信される信号に基づいて測定する行列測定ステップと、前記伝搬路の間に対象物が存在する場合に前記第1の行列の時間応答を周波数応答に変換して得られる第2の行列及び前記伝搬路の間に対象物が存在しない場合に前記第1の行列の時間応答を周波数応答に変換して得られる第3の行列を算出する行列算出ステップと、前記第2の行列の各要素に該第2の行列の各要素の複素共役を掛け合わせて得られる結果及び前記第3の行列に該第2の行列の各要素の複素共役を掛け合わせて得られる結果を用いて所定の条件を満たすか否かを判定し、該所定の条件を満たしていれば前記対象物が活動していると判定する判定ステップと、を含む、生体活動測定方法が提供される。 In order to solve the above-described problem, according to another aspect of the present invention, a first matrix indicating a propagation path formed between a plurality of transmission antennas and a plurality of reception antennas is represented as the plurality of transmissions. A matrix measurement step of measuring based on signals transmitted from an antenna and received by the plurality of reception antennas, and a time response of the first matrix as a frequency response when an object exists between the propagation paths Matrix calculation step of calculating a third matrix obtained by converting the time response of the first matrix into a frequency response when there is no object between the second matrix obtained by conversion and the propagation path And the result obtained by multiplying each element of the second matrix by the complex conjugate of each element of the second matrix, and the complex conjugate of each element of the second matrix multiplied by the third matrix Using the results obtained It determines whether they meet, if they meet the predetermined condition including a determination step of determining said object is active, the biological activity measuring method is provided.
本発明の別の観点による生体活動測定方法は、かかるステップを含むことによって、複数の送信アンテナと複数の受信アンテナとの間で電波を送受信するMIMO技術を利用して、より正確に人物の活動の状態を推定することが可能となる。 The life activity measuring method according to another aspect of the present invention includes such steps, thereby more accurately using the MIMO technology that transmits and receives radio waves between a plurality of transmitting antennas and a plurality of receiving antennas. It is possible to estimate the state.
以上説明したように本発明によれば、MIMO技術を利用して、正確に状態を推定することが可能な、新規かつ改良された生体情報測定装置及び生体情報測定方法を提供することができる。 As described above, according to the present invention, it is possible to provide a new and improved biological information measuring apparatus and biological information measuring method capable of accurately estimating a state using MIMO technology.
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。 Exemplary embodiments of the present invention will be described below in detail with reference to the accompanying drawings. In addition, in this specification and drawing, about the component which has the substantially same function structure, duplication description is abbreviate | omitted by attaching | subjecting the same code | symbol.
<1.本発明の一実施形態>
[生体情報測定システム1の全体構成例]
まず、図面を参照しながら本発明の一実施形態に係る生体情報測定装置を含んだ生体情報測定システムの機能構成例について説明する。図1は、本発明の一実施形態に係る生体情報測定装置100を含んだ生体情報測定システム1の全体構成例を示す説明図である。以下、図1を用いて生体情報測定システム1の全体構成例について説明する。
<1. One Embodiment of the Present Invention>
[Example of overall configuration of biological information measuring system 1]
First, a functional configuration example of a biological information measuring system including a biological information measuring apparatus according to an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is an explanatory diagram showing an example of the overall configuration of a biological information measuring system 1 including a biological
図1に示した生体情報測定システム1は、MIMO(Multiple Input Multiple Output)技術を利用して、生体情報測定システム1が置かれた環境(例えば部屋の中)における人物2の生体活動を検出することを目的としたシステムである。図1に示したように、生体情報測定システム1は、送信機10と、生体活動測定装置100と、を含んで構成される。
The biological information measurement system 1 shown in FIG. 1 detects the biological activity of the
送信機10は、M本(M≧2)のアンテナ11a、11b、・・・、11mを備え、アンテナ11a、11b、・・・、11mからそれぞれ所定の周波数帯の電波を送信する。生体活動測定装置100は、N本(N≧2)のアンテナ101a、101b、・・・、101nを備え、送信機10のアンテナ11a、11b、・・・、11mから送信された所定の周波数帯の電波を受信する。
The
生体情報測定システム1が置かれた環境に人物2が存在すれば、送信機10のアンテナ11a、11b、・・・、11mから送信された電波は、人物2の体の表面で反射される場合がある。生体活動測定装置100は、送信機10から送信された電波を受信し、伝搬チャネル行列を解析することで、人物2の活動の有無や、人物2が活動する場合のその人物2の生体活動を測定する。
If the
なお、本発明において対象物は人物に限られない。例えば、本発明における対象物としては、固有の周波数成分を検出することが可能な任意の物体が挙げられる。以下では、本発明の実施形態に係る対象物が、人体である場合を例に挙げる。 In the present invention, the object is not limited to a person. For example, the object in the present invention includes any object capable of detecting a specific frequency component. Below, the case where the target object which concerns on embodiment of this invention is a human body is mentioned as an example.
本発明の一実施形態に係る生体活動測定装置100は、人物2の生体活動を測定する際に、人物2の生体活動をより正確に測定することを目的とするものである。MIMOを利用して人物の生体活動を測定しようとする場合には、人物がいない場合の伝搬チャネルの周波数特性が、測定するチャネルや伝搬環境の違いにより変化することがあり、また、伝搬チャネルの各要素での値が大きく異なることもある。基準にするチャネルをどれにするかで、人が入ることにより呼吸や心拍を推定するためのしきい値が変化する。この値の設定により、人がいない場合でも呼吸や心拍が誤って推定される場合が生じる。
The life
詳細な処理については後述するが、本発明の一実施形態に係る生体活動測定装置100は、送信機10との間の電波の伝搬路に、図1に示したように人物が存在している場合の伝搬チャネルの周波数特性と、人物が存在しない場合の伝搬チャネルの周波数特性と、をそれぞれ算出する。本発明の一実施形態に係る生体活動測定装置100は、人物が存在している場合の伝搬チャネルの周波数特性の各要素に対し、その各要素の複素共役を掛け合わせる。また本発明の一実施形態に係る生体活動測定装置100は、人物が存在していない場合の伝搬チャネルの周波数特性の各要素に対し、その人物が存在している場合の伝搬チャネルの周波数特性の各要素の複素共役を掛け合わせる。
Although detailed processing will be described later, in the life
本発明の一実施形態に係る生体活動測定装置100は、このように、人物が存在している場合の伝搬チャネルの周波数特性と、人物が存在しない場合の伝搬チャネルの周波数特性とに対する演算を行なうことで、人物がいないのに、ノイズの変化を人物の生体活動(例えば呼吸の実行)と誤検出してしまうことを防ぐことが出来る。また本発明の一実施形態に係る生体活動測定装置100は、人物がいる場合であっても、人物が呼吸を止めている期間なのに、ノイズの変化を人物の生体活動(例えば呼吸の実行)と誤検出してしまうことを防ぐことが出来る。
As described above, the life
以上、図1を用いて生体情報測定システム1の全体構成例について説明した。続いて、本発明の一実施形態に係る生体活動測定装置100の機能構成例について説明する。
The overall configuration example of the biological information measurement system 1 has been described above with reference to FIG. Subsequently, a functional configuration example of the life
[生体情報測定装置の機能構成例]
図2は、本発明の一実施形態に係る生体活動測定装置100の機能構成例を示す説明図である。以下、図2を用いて本発明の一実施形態に係る生体活動測定装置100の機能構成例について説明する。
[Functional configuration example of biological information measuring apparatus]
FIG. 2 is an explanatory diagram showing a functional configuration example of the life
図2に示した本発明の一実施形態に係る生体活動測定装置100は、送信機10のアンテナ11a、11b、・・・、11mから送信された電波を受信して、伝搬チャネル行列を解析することで、人物2の活動の有無や、人物2が活動する場合のその人物2の生体活動を測定することを目的とした装置である。
The life
図2に示したように、本発明の一実施形態に係る生体活動測定装置100は、前処理部110と、FFT(Fast Fourier Transform;高速フーリエ変換)処理部120と、判定部130と、を含んで構成される。
As shown in FIG. 2, the life
前処理部110は、アンテナ101a、101b、・・・、101nで受信した電波に対する前処理を実行する。前処理部110が実行する処理には、例えば、アンテナ101a、101b、・・・、101nが受信した信号のベースバンド信号への変換、N行M列の伝搬チャネル行列H(t)の測定が含まれ得る。j番目の送信アンテナとi番目の受信アンテナの間の伝搬路利得hij(t)を要素に持つN行M列の伝搬チャネル行列H(t)は、以下の数式(1)で表わされる。なお、数式1に示す“t”は、伝搬チャネル行列H(t)が取得された時間、すなわち、伝搬チャネル行列H(t)の算出元となる信号が受信された時間を示している。
The
前処理部110は、伝搬チャネル行列H(t)を測定すると、伝搬チャネル行列H(t)をFFT処理部120へ出力する。以上から、前処理部110は本発明の行列測定部の一例として機能し得る。
When the
FFT処理部120は、前処理部110が出力した時刻tの伝搬チャネル行列H(t)に対して高速フーリエ変換(FFT)処理を行う。FFT処理部120は、時間tにおける伝搬チャネル行列H(t)の時系列に沿った変化、すなわち、hij(t)の時間応答を高速フーリエ変換することによって、周波数ごとに各伝搬路の周波数応答Fij(f)を算出する。数式(1)の伝搬チャネル行列H(t)に対して高速フーリエ変換を実行すると、数式(2)で示す行列F(f)が得られる。
The
対象物が人体の場合には、数式2で示した行列F(f)には、呼吸や心拍などの人体の生体反応に固有の周波数が含まれる周波数成分が含まれ得る。FFT処理部120は、前処理部110が出力した時刻tの伝搬チャネル行列H(t)に対する高速フーリエ変換処理を、生体情報測定システム1が置かれた環境に人物2が存在する場合としない場合に対して行なう。生体情報測定システム1が置かれた環境に人物2が存在しない場合に、FFT処理部120が伝搬チャネル行列H(t)に対して高速フーリエ変換を行った結果をF(no)(f)、人物2が存在する場合で、FFT処理部120が伝搬チャネル行列H(t)に対して高速フーリエ変換を行った結果をF(f)とする。
When the object is a human body, the matrix F (f) expressed by
FFT処理部120は、伝搬チャネル行列H(t)に対して高速フーリエ変換を行って得られるF(f)及びF(no)(f)を判定部130に送る。以上から、FFT処理部120は本発明の行列算出部の一例として機能し得る。
The
判定部130は、FFT処理部120から送られるF(f)及びF(no)(f)を用いて、人物2の活動の有無の判定、人物2の生体情報の測定を実行する。
Using the F (f) and F (no) (f) sent from the
判定部130は、F(f)及びF(no)(f)に対して、それぞれ以下の数式(3)及び(4)の数式に基づいた関数の算出処理を実行する。なお、数式(3)及び(4)において、*は複素共役を表している。すなわち判定部130は、F(f)に対してはF(f)の複素共役を、F(no)(f)に対してはF(f)の複素共役を、それぞれ掛け合わせる処理を実行する。
The
続いて判定部130は、数式(3)及び(4)で得られるYij(f)及びYij (no)(f)を用いて、以下の数式(5)によってρijを算出する。
Subsequently, the
判定部130は、上記数式(5)で算出したρijが、iが1からNの間及びjが1からMの間で満たされれば、人物2が活動していると判定する。そして判定部130は、人物2が活動していると判定すると、対象周波数から心拍または呼吸の数を計数する。
The
上記数式(3)及び数式(4)が意味するところは、以下の通りである。単にF(no)(f)及びF(f)のそれぞれに対し、それぞれの複素共役が掛け合わされると、虚数部が無くなり値が最大化される。F(no)(f)及びF(f)のそれぞれに対し、それぞれの複素共役が掛け合わされる場合、人物2がいない場合であっても、ノイズの変化を人物2の生体活動(例えば呼吸の実行)と誤検出してしまうおそれがある。
The meanings of the above formulas (3) and (4) are as follows. Simply multiplying each of F (no) (f) and F (f) by their respective complex conjugates eliminates the imaginary part and maximizes the value. When F (no) (f) and F (f) are multiplied by their respective complex conjugates, even if no
しかし、上記数式(4)のように、Fij (no)(f)に、Fij(f)の複素共役Fij(f)*を掛け合わせることで、仮に呼吸成分がFij(f)に含まれていれば、ノイズ成分と掛け合わされることで、Yij (no)(f)の方が小さくなる。一方、呼吸成分がFij(f)に含まれていなければ、Yij (no)(f)はノイズと同等のレベルになるので、Yij (no)(f)の方が大きくなる。 However, as in the above equation (4), F ij (no ) to (f), F ij complex conjugate F ij of (f) (f) * By multiplying the, if respiration component F ij (f) If it is included, Y ij (no) (f) becomes smaller by being multiplied by the noise component. On the other hand, if the respiratory component is not included in F ij (f), Y ij (no) (f) is at the same level as noise, so Y ij (no) (f) is larger.
従って、上記数式(4)のように、Fij (no)(f)に、Fij(f)の複素共役Fij(f)*を掛け合わせることで、判定部130は、人物2がいない場合にノイズの変化を人物2の生体活動(例えば呼吸の実行)と誤検出してしまうことを防ぐことが出来る。また判定部130は、人物2がいる場合であっても、人物2が呼吸を止めている期間なのにノイズの変化を人物2の呼吸の実行と誤検出してしまうことを防ぐことが出来る。
Therefore, as in the above equation (4), the F ij (no) (f), by multiplying the complex conjugate F ij (f) * of F ij (f), the
なお判定部130は、以下の数式(6)によって算出されるρ’の条件を満たした場合に人物2が活動していると判定してもよい。チャネル行列H(t)の全ての要素で値が得られない場合があるので、そのような場合は、判定部130は、数式(5)ではなく、下記の数式(6)によって人物2の生体活動の有無やその生体活動の内容を判定してもよい。そして判定部130は、人物2が活動していると判定すると、対象周波数から人物2の心拍または呼吸の数を計数する。
Note that the
本発明の一実施形態に係る生体活動測定装置100は、図2に示したような構成を有することで、MIMO技術を利用して、より正確に人物の活動の状態を推定することが出来る。すなわち、本発明の一実施形態に係る生体活動測定装置100は、図2に示したような構成を有することで、人物が存在している場合の伝搬チャネルの周波数特性と、人物が存在しない場合の伝搬チャネルの周波数特性とに対する演算を行ない、演算結果が所定の条件を満たしているか否かで人物の生体活動の有無やその内容を判定することができる。
The life
従って、本発明の一実施形態に係る生体活動測定装置100は、人物がいないのに、ノイズの変化を人物の生体活動(例えば呼吸の実行)と誤検出してしまうことを防ぐことが出来る。また本発明の一実施形態に係る生体活動測定装置100は、人物がいる場合であっても、例えば人物が生体活動(例えば呼吸の実行)を止めている期間なのに、ノイズの変化を人物の生体活動(例えば呼吸の実行)と誤検出してしまうことを防ぐことが出来る。
Therefore, the life
以上、図2を用いて本発明の一実施形態に係る生体活動測定装置100の機能構成例について説明した。続いて、本発明の一実施形態に係る生体活動測定装置100の動作例について説明する。
The functional configuration example of the life
[生体情報測定装置の動作例]
図3は、本発明の一実施形態に係る生体活動測定装置100の動作例を示す流れ図である。図3に示した流れ図は、送信機10のアンテナ11a、11b、・・・、11mから送信された電波を生体活動測定装置100で受信して、伝搬チャネル行列を解析することで、人物2の活動の有無や、人物2が活動する場合のその人物2の生体活動を生体活動測定装置100で測定する際の動作例である。以下、図3を用いて本発明の一実施形態に係る生体活動測定装置100の動作例について説明する。
[Operation example of biological information measuring device]
FIG. 3 is a flowchart showing an operation example of the life
生体活動測定装置100は、送信機10のアンテナ11a、11b、・・・、11mから送信された電波をアンテナ101a、101b、・・・、101nで受信すると、受信した電波に対して前処理を実行する(ステップS101)。このステップS101の前処理は前処理部110が実行する。
The life
前処理部110が実行する処理は、例えば、アンテナ101a、101b、・・・、101nが受信した信号のベースバンド信号への変換、N行M列の伝搬チャネル行列H(t)の測定を含む。j番目の送信アンテナとi番目の受信アンテナの間の伝搬路利得hij(t)を要素に持つ、時刻tにおけるN行M列の伝搬チャネル行列H(t)は、上記数式(1)で表すことができる。
The processing executed by the
生体活動測定装置100は、受信した電波に対する前処理を上記ステップS101で実行すると、続いて、前処理によって得られた伝搬チャネル行列を高速フーリエ変換する(ステップS102)。この高速フーリエ変換は、FFT処理部120が実行する。
When the biological
FFT処理部120は、上述したように、時刻tにおける伝搬チャネル行列H(t)に対する上記数式(2)による高速フーリエ変換を、生体情報測定システム1が置かれた環境に人物2が存在する場合としない場合に対して行なう。生体情報測定システム1が置かれた環境に人物2が存在しない場合に、FFT処理部120が伝搬チャネル行列H(t)に対して高速フーリエ変換を行った結果得られる行列をF(no)(f)、人物2が存在する場合で、FFT処理部120が伝搬チャネル行列H(t)に対して高速フーリエ変換を行った結果得られる行列をF(f)とする。
As described above, the
生体活動測定装置100は、伝搬チャネル行列H(t)に対する高速フーリエ変換によって、行列F(no)(f)及び行列F(f)を上記ステップS102で算出すると、続いて、行列F(no)(f)及びF(f)を用いて人物の活動の判定処理を実行する(ステップS103)。ステップS103の判定処理は判定部130が実行する。
The life
判定部130は、上記の数式(3)及び(4)で得られる行列Yij(f)及びYij (no)(f)を用いて、上記の数式(5)によってρijを算出するか、または上記の数式(6)によってρ’を算出する。判定部130は、数式(5)または数式(6)のどちらの条件を用いて、人物2の生体活動の有無を判定してもよく、両方の条件を用いて人物2の生体活動の有無を判定してもよい。
Whether the
判定部130は、例えば数式(5)のρijの条件が、iが1からNの間及びjが1からMの間で満たされれば、人物2が存在した上で、人物2が生体活動を行っていると判定する。そして判定部130は、人物2が存在した上で、人物2が生体活動を行っていると判定すると、対象周波数から心拍または呼吸の数を計数する。また判定部130は、数式(6)のρ’の条件が満たされれば、人物2が活動していると判定し、対象周波数から心拍または呼吸の数を計数する。
For example, if the condition of ρ ij in Equation (5) satisfies the condition where i is between 1 and N and j is between 1 and M, then the
また判定部130は、例えば数式(6)の条件が満たされれば、上記数式(5)の条件の場合と同様に、人物2が存在した上で、人物2が生体活動を行っていると判定する。
Further, for example, if the condition of Expression (6) is satisfied, the
一方、判定部130は、例えば数式(5)のρijの条件が、iが1からNの間及びjが1からMの間で満たされなければ、人物2が存在しない、または人物2が存在していても、生体活動を行っていないと判定する。また判定部130は、例えば数式(6)の条件が満たされなければ、上記数式(5)の条件の場合と同様に、人物2が存在しない、または人物2が存在していても、生体活動を行っていないと判定する。
On the other hand, for example, if the condition of ρ ij in Expression (5) does not satisfy i between 1 and N and j between 1 and M, the
本発明の一実施形態に係る生体活動測定装置100は、図3に示したような動作を実行することで、MIMO技術を利用して、より正確に人物の活動の状態を推定することが出来る。すなわち、本発明の一実施形態に係る生体活動測定装置100は、図3に示したような動作を実行することで、人物が存在している場合の伝搬チャネルの周波数特性と、人物が存在しない場合の伝搬チャネルの周波数特性とに対する演算を行ない、演算結果が所定の条件を満たしているか否かで人物の生体活動の有無やその内容を判定することができる。
The life
ここで、従来のMIMO技術を利用した呼吸数の推定結果と、本発明の一実施形態に係る生体活動測定装置100による呼吸数の推定結果とを比較することで、本発明の一実施形態に係る生体活動測定装置100が、従来のMIMO技術を利用した場合と比較して、より正確に人物の活動の状態を推定できていることを示す。
Here, by comparing the respiration rate estimation result using the conventional MIMO technique with the respiration rate estimation result by the life
図4は、従来のMIMO技術を利用した呼吸数の推定結果をグラフで示す説明図であり、図5は、本発明の一実施形態に係る生体活動測定装置100による呼吸数の推定結果をグラフで示す説明図である。
FIG. 4 is an explanatory diagram illustrating a result of respiratory rate estimation using a conventional MIMO technique, and FIG. 5 is a graph illustrating a respiratory rate estimation result by the life
図4に示したグラフでは、人物が息を止めている期間であっても、誤って呼吸数が測定されてしまっている。これに対して本発明の一実施形態に係る生体活動測定装置100では、人物が息を止めている期間では呼吸数が測定されることはない。
In the graph shown in FIG. 4, the respiration rate is erroneously measured even during the period in which the person holds his / her breath. On the other hand, in the life
このように、本発明の一実施形態に係る生体活動測定装置100は、従来のMIMO技術を利用した場合と比較して、より正確に人物の活動の状態を推定できていることが分かる。
Thus, it can be seen that the life
<2.まとめ>
以上説明したように本発明の一実施形態によれば、送信機10のアンテナ11a、11b、・・・、11mから送信された電波を受信して、伝搬チャネル行列を解析することで、人物2の活動の有無や、人物2が活動する場合のその人物2の生体活動を測定する、生体活動測定装置100が提供される。
<2. Summary>
As described above, according to one embodiment of the present invention, the radio waves transmitted from the
本発明の一実施形態に係る生体活動測定装置100は、伝搬チャネル行列H(t)に対する上記数式(2)による高速フーリエ変換を、生体情報測定システム1が置かれた環境に人物2が存在する場合としない場合に対して実行する。そして本発明の一実施形態に係る生体活動測定装置100は、高速フーリエ変換の結果を用いて上述の数式(5)または数式(6)の条件を満たした場合に、人物2が活動していると判定するとともに、人物2が活動していると判定すると、対象周波数から心拍または呼吸の数を計数する。
A life
本発明の一実施形態に係る生体活動測定装置100は、伝搬チャネル行列H(t)に対する高速フーリエ変換を、生体情報測定システム1が置かれた環境に人物2が存在する場合としない場合に対して行なう。そして、上記数式(4)のように、F(no)(f)に対して、F(f)の複素共役を掛け合わせる。
The life
本発明の一実施形態に係る生体活動測定装置100は、生体情報測定システム1が置かれた環境に人物2が存在しない場合における、伝搬チャネル行列H(t)に対する高速フーリエ変換の結果に対して、人物2が存在しない場合における、伝搬チャネル行列H(t)に対する高速フーリエ変換を用いることで、ノイズの変化による影響を小さくし、人物2がいない場合にノイズの変化を人物2の生体活動(例えば呼吸の実行)と誤検出してしまったり、人物2がいる場合であっても人物2が呼吸を止めている期間なのにノイズの変化を人物2の呼吸の実行と誤検出してしまったりすることを防ぐことが出来る。
The life
本発明の一実施形態に係る生体活動測定装置100は、このような動作を実行することで、従来のMIMO技術を利用した場合と比較して、より正確に人物の活動の状態を推定することが可能となる。
The life
本明細書における装置が実行する処理における各ステップは、必ずしもフローチャートやシーケンス図として記載された順序に沿って時系列に処理される必要はない。例えば、各装置が実行する処理における各ステップは、フローチャートとして記載した順序と異なる順序で処理されても良く、並列的に処理されてもよい。 Each step in the processing executed by the apparatus in the present specification does not necessarily have to be processed in time series in the order described as a flowchart or a sequence diagram. For example, each step in the processing executed by each device may be processed in an order different from the order described as the flowchart, or may be processed in parallel.
また上述した各装置の構成と同等の機能を、本明細書における装置に内蔵されるCPU(Central Processing Unit)、ROM(Read Only Memory)およびRAM(Random Access Memory)などのハードウェア(図示せず)に発揮させるためのコンピュータプログラムも作成されることができる。また、該コンピュータプログラムを記憶させた記憶媒体も提供されることが可能である。また、機能ブロック図で示したそれぞれの機能ブロックをハードウェアで構成することで、一連の処理をハードウェアで実現することもできる。 In addition, functions (e.g., a CPU (Central Processing Unit), a ROM (Read Only Memory), and a RAM (Random Access Memory) incorporated in the device in this specification have functions equivalent to the configuration of each device described above. ) Can be created. A storage medium storing the computer program can also be provided. Moreover, a series of processes can also be realized by hardware by configuring each functional block shown in the functional block diagram with hardware.
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。 The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to such examples. It is obvious that a person having ordinary knowledge in the technical field to which the present invention pertains can come up with various changes or modifications within the scope of the technical idea described in the claims. Of course, it is understood that these also belong to the technical scope of the present invention.
例えば上記実施形態では、生体情報測定装置100がアンテナ101a、101b、・・・、101nを有している構成を例示したが、本発明は係る例に限定されるものではない。例えば、生体情報測定装置100にはアンテナが設けられておらず、生体情報測定装置100とは別の装置が送信機10からの電波を受信するような形態を採っても良い。
For example, in the above-described embodiment, the configuration in which the biological
1 生体情報測定システム
2 人物
10 送信機
11a、11b、・・・、11m アンテナ
100 生体情報測定装置
101a、101b、・・・、101n アンテナ
110 前処理部
120 FFT処理部
130 判定部
DESCRIPTION OF SYMBOLS 1 Biological
Claims (4)
前記伝搬路の間に対象物が存在する場合に前記第1の行列の時間応答を周波数応答に変換して得られる第2の行列及び前記伝搬路の間に対象物が存在しない場合に前記第1の行列の時間応答を周波数応答に変換して得られる第3の行列を算出する行列算出部と、
前記第2の行列の要素毎に該第2の行列の各要素の複素共役を掛け合わせて得られる第1の結果及び前記第3の行列に該第2の行列の各要素の複素共役を掛け合わせて得られる第2の結果を用いて所定の条件を満たすか否かを判定し、該所定の条件を満たしていれば前記対象物が活動していると判定する判定部と、
を備える、生体活動測定装置。 A first matrix indicating a propagation path formed between a plurality of transmission antennas and a plurality of reception antennas is measured based on signals transmitted from the plurality of transmission antennas and received by the plurality of reception antennas A matrix measurement unit;
The second matrix obtained by converting the time response of the first matrix into a frequency response when an object is present between the propagation paths and the object when the object is not present between the propagation paths. A matrix calculator that calculates a third matrix obtained by converting the time response of the matrix of 1 into a frequency response;
The first result obtained by multiplying each element of the second matrix by the complex conjugate of each element of the second matrix and the third matrix are multiplied by the complex conjugate of each element of the second matrix. A determination unit that determines whether or not a predetermined condition is satisfied using the second result obtained together, and determines that the object is active if the predetermined condition is satisfied;
A life activity measuring apparatus comprising:
前記伝搬路の間に対象物が存在する場合に前記第1の行列の時間応答を周波数応答に変換して得られる第2の行列及び前記伝搬路の間に対象物が存在しない場合に前記第1の行列の時間応答を周波数応答に変換して得られる第3の行列を算出する行列算出ステップと、
前記第2の行列の各要素に該第2の行列の各要素の複素共役を掛け合わせて得られる結果及び前記第3の行列に該第2の行列の各要素の複素共役を掛け合わせて得られる結果を用いて所定の条件を満たすか否かを判定し、該所定の条件を満たしていれば前記対象物が活動していると判定する判定ステップと、
を含む、生体活動測定方法。 A first matrix indicating a propagation path formed between a plurality of transmission antennas and a plurality of reception antennas is measured based on signals transmitted from the plurality of transmission antennas and received by the plurality of reception antennas A matrix measurement step;
The second matrix obtained by converting the time response of the first matrix into a frequency response when an object is present between the propagation paths and the object when the object is not present between the propagation paths. A matrix calculating step of calculating a third matrix obtained by converting the time response of the matrix of 1 into a frequency response;
The result obtained by multiplying each element of the second matrix by the complex conjugate of each element of the second matrix and the result obtained by multiplying the complex conjugate of each element of the second matrix by the third matrix. A determination step of determining whether or not a predetermined condition is satisfied using a result obtained, and determining that the object is active if the predetermined condition is satisfied;
A method for measuring life activity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013264172A JP6222454B2 (en) | 2013-12-20 | 2013-12-20 | Biological information measuring apparatus and biological information measuring method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013264172A JP6222454B2 (en) | 2013-12-20 | 2013-12-20 | Biological information measuring apparatus and biological information measuring method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2015119770A true JP2015119770A (en) | 2015-07-02 |
JP6222454B2 JP6222454B2 (en) | 2017-11-01 |
Family
ID=53532048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013264172A Expired - Fee Related JP6222454B2 (en) | 2013-12-20 | 2013-12-20 | Biological information measuring apparatus and biological information measuring method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6222454B2 (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3193189A1 (en) | 2016-01-15 | 2017-07-19 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor and direction estimation method |
JP2018036247A (en) * | 2016-08-29 | 2018-03-08 | パナソニックIpマネジメント株式会社 | Identification device and identification method |
KR101839686B1 (en) | 2016-06-27 | 2018-03-16 | 한국과학기술원 | Method, system and non-transitory computer-readable recording medium for monitoring hemodynamics |
JP2018112540A (en) * | 2017-01-06 | 2018-07-19 | パナソニックIpマネジメント株式会社 | Sensor and method |
US10241187B2 (en) | 2015-12-24 | 2019-03-26 | Panasonic Intellectual Property Management Co., Ltd. | Position sensor, direction estimation method, and system |
JP2020109389A (en) * | 2018-12-28 | 2020-07-16 | パナソニックIpマネジメント株式会社 | Estimation method, estimation device, and program |
US10914829B2 (en) | 2017-01-27 | 2021-02-09 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor, sensor, and method |
WO2022075467A1 (en) * | 2020-10-09 | 2022-04-14 | サクラテック株式会社 | Portable contactless biosignal detection device, driver monitor device, visitor screening system, and household healthcare system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH116874A (en) * | 1997-06-16 | 1999-01-12 | Tau Giken:Kk | Method and device for searching buried living body |
WO2005050786A1 (en) * | 2003-11-17 | 2005-06-02 | Jean-Marc Cortambert | Cross antenna comprising linear sub-antennas and associated processing |
JP2005156393A (en) * | 2003-11-27 | 2005-06-16 | Mitsubishi Electric Corp | Radar equipment |
US20090203972A1 (en) * | 2006-06-01 | 2009-08-13 | Biancamed Ltd. | Apparatus, system, and method for monitoring physiological signs |
WO2009124297A1 (en) * | 2008-04-03 | 2009-10-08 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
JP2012145332A (en) * | 2011-01-06 | 2012-08-02 | Mitsubishi Electric Corp | Mimo radar apparatus |
-
2013
- 2013-12-20 JP JP2013264172A patent/JP6222454B2/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH116874A (en) * | 1997-06-16 | 1999-01-12 | Tau Giken:Kk | Method and device for searching buried living body |
WO2005050786A1 (en) * | 2003-11-17 | 2005-06-02 | Jean-Marc Cortambert | Cross antenna comprising linear sub-antennas and associated processing |
JP2005156393A (en) * | 2003-11-27 | 2005-06-16 | Mitsubishi Electric Corp | Radar equipment |
US20090203972A1 (en) * | 2006-06-01 | 2009-08-13 | Biancamed Ltd. | Apparatus, system, and method for monitoring physiological signs |
WO2009124297A1 (en) * | 2008-04-03 | 2009-10-08 | Kai Sensors, Inc. | Non-contact physiologic motion sensors and methods for use |
JP2012145332A (en) * | 2011-01-06 | 2012-08-02 | Mitsubishi Electric Corp | Mimo radar apparatus |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10241187B2 (en) | 2015-12-24 | 2019-03-26 | Panasonic Intellectual Property Management Co., Ltd. | Position sensor, direction estimation method, and system |
US11150336B2 (en) | 2016-01-15 | 2021-10-19 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor and direction estimation method |
US10371808B2 (en) | 2016-01-15 | 2019-08-06 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor and direction estimation method |
US10663573B2 (en) | 2016-01-15 | 2020-05-26 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor and direction estimation method |
EP3193189A1 (en) | 2016-01-15 | 2017-07-19 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor and direction estimation method |
KR101839686B1 (en) | 2016-06-27 | 2018-03-16 | 한국과학기술원 | Method, system and non-transitory computer-readable recording medium for monitoring hemodynamics |
JP2018036247A (en) * | 2016-08-29 | 2018-03-08 | パナソニックIpマネジメント株式会社 | Identification device and identification method |
JP2018112540A (en) * | 2017-01-06 | 2018-07-19 | パナソニックIpマネジメント株式会社 | Sensor and method |
US10914829B2 (en) | 2017-01-27 | 2021-02-09 | Panasonic Intellectual Property Management Co., Ltd. | Positioning sensor, sensor, and method |
JP2020109389A (en) * | 2018-12-28 | 2020-07-16 | パナソニックIpマネジメント株式会社 | Estimation method, estimation device, and program |
JP7349661B2 (en) | 2018-12-28 | 2023-09-25 | パナソニックIpマネジメント株式会社 | Estimation method, estimation device and program |
WO2022075467A1 (en) * | 2020-10-09 | 2022-04-14 | サクラテック株式会社 | Portable contactless biosignal detection device, driver monitor device, visitor screening system, and household healthcare system |
JP7086441B1 (en) * | 2020-10-09 | 2022-06-20 | サクラテック株式会社 | Portable non-contact biosignal detectors, driver monitoring devices, attendee screening systems, and home healthcare systems |
Also Published As
Publication number | Publication date |
---|---|
JP6222454B2 (en) | 2017-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6222454B2 (en) | Biological information measuring apparatus and biological information measuring method | |
JP6865394B2 (en) | program | |
US11150336B2 (en) | Positioning sensor and direction estimation method | |
JP6290501B2 (en) | High frequency equipment | |
TWI476414B (en) | Distributed spectrum sensing | |
US20210255280A1 (en) | Radar vital signal tracking using a kalman filter | |
JP6975898B2 (en) | Estimating method, estimation device and program | |
US10413229B2 (en) | Biometric device and biometric method | |
US20210353156A1 (en) | Direct rf signal processing for heart-rate monitoring using uwb impulse radar | |
EP3342330B1 (en) | Phasor approach to signal to noise ratio measurement evaluation of physiological signals | |
JP7511189B2 (en) | Estimation device, estimation method, and program | |
Wang et al. | ResBeat: Resilient breathing beats monitoring with realtime bimodal CSI data | |
JP7349661B2 (en) | Estimation method, estimation device and program | |
CN111759305A (en) | Living body detection method and living body detection system | |
EP3384843A1 (en) | Determining breathing attributes by a detection device | |
JP2014228291A (en) | Radio detection device and radio detection method | |
EP3184031A1 (en) | Perspiration and heart rate estimation system, perspiration and heart rate estimation device, and perspiration and heart rate estimation method | |
JP7162192B2 (en) | Bionumber estimation device, bionumber estimation method, and program | |
JP7313018B2 (en) | sensor | |
US11047968B2 (en) | Estimating method and estimating device | |
JP7357217B2 (en) | Estimation method, estimation device, and program | |
EP3207864A1 (en) | Sensor information processing apparatus | |
CN110693454A (en) | Sleep characteristic event detection method and device based on radar and storage medium | |
US11561279B2 (en) | Radar estimating method, device and medium to extract living body vector information | |
Cho et al. | Novel design for heart rate detection using UWB impulse radar on Android platform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20161208 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20161208 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170224 |
|
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: 20170822 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170823 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20170919 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170919 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20170919 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6222454 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |