JP2019211458A - Estimation device and estimation method - Google Patents

Estimation device and estimation method Download PDF

Info

Publication number
JP2019211458A
JP2019211458A JP2018160851A JP2018160851A JP2019211458A JP 2019211458 A JP2019211458 A JP 2019211458A JP 2018160851 A JP2018160851 A JP 2018160851A JP 2018160851 A JP2018160851 A JP 2018160851A JP 2019211458 A JP2019211458 A JP 2019211458A
Authority
JP
Japan
Prior art keywords
transmission
living body
antenna element
antenna elements
predetermined range
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
Application number
JP2018160851A
Other languages
Japanese (ja)
Other versions
JP7129648B2 (en
Inventor
翔一 飯塚
Shoichi Iizuka
翔一 飯塚
武司 中山
Takeshi Nakayama
武司 中山
尚樹 本間
Naoki Honma
尚樹 本間
大 笹川
Dai Sasagawa
大 笹川
之雄 袁
Zhixiong Chen
之雄 袁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Intellectual Property Management Co Ltd
Original Assignee
Panasonic Intellectual Property Management Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Panasonic Intellectual Property Management Co Ltd filed Critical Panasonic Intellectual Property Management Co Ltd
Priority to CN201910011041.3A priority Critical patent/CN110018522B/en
Priority to US16/240,919 priority patent/US10969479B2/en
Publication of JP2019211458A publication Critical patent/JP2019211458A/en
Application granted granted Critical
Publication of JP7129648B2 publication Critical patent/JP7129648B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

To provide an estimation device capable of improving resolution and accuracy of a biological body direction estimated.SOLUTION: An estimation device 10 for estimating a biological body direction comprises: transmission/reception parts 20A-20H for transmitting a transmitter signal using M transmission antenna elements which are arranged to surround a prescribed range A1 including a biological body 100, and receiving a reception signal using N reception antenna elements; and a circuit 40 for calculating a feature amount which indicates a large amplitude and more regular waveform as N reception signals have a larger value, for each of M set of N reception signals respectively corresponding to M transmission signals, comparing M feature amounts acquired by calculation for identifying a first transmission antenna element corresponding to the first feature amount which is maximum out of M feature amounts, and then, estimating the prescribed direction on the basis of the identified first transmission antenna element as the biological body direction.SELECTED DRAWING: Figure 1

Description

本開示は、生体に無線信号を照射し、その反射信号を受信して生体の正面が向いている方向である生体方向を推定する推定装置及び推定方法に関する。   The present disclosure relates to an estimation apparatus and an estimation method that irradiate a living body with a radio signal, receive a reflection signal thereof, and estimate a living body direction that is a direction in which the front of the living body is facing.

特許文献1には、ベッド上の生体に対して電磁波を照射し、その反射波に含まれる心拍及び呼吸の振動を検出することで、生体が仰向けであるかうつ伏せであるかを推定する装置が開示されている。   Patent Document 1 discloses an apparatus that estimates whether a living body is lying on its back or lying down by irradiating a living body on a bed with electromagnetic waves and detecting vibrations of heartbeat and respiration contained in the reflected wave. It is disclosed.

また、特許文献2には、災害救助の現場において瓦礫の山に向けて電磁波を照射し、その反射波から生存者の有無を検出する装置が開示されている。   Patent Document 2 discloses an apparatus that irradiates electromagnetic waves toward a pile of rubble at a disaster rescue site and detects the presence or absence of a survivor from the reflected waves.

特開2014−207935号公報JP 2014-207935 A 特開平11−6874号公報Japanese Patent Laid-Open No. 11-6874

しかしながら、従来技術では、生体の方向をより細かい分解能で推定することができない。   However, in the prior art, the direction of the living body cannot be estimated with a finer resolution.

上記目的を達成するために、本開示に係る推定装置は、生体の正面が向いている方向である生体方向を推定する推定装置であって、前記生体を含む所定範囲を囲んで配置されるM個(Mは3以上)の送信アンテナ素子を用いて前記所定範囲にM個の送信信号を送信する送信部と、前記所定範囲の周囲を囲んで配置されるN個(Nは2以上)の受信部であって、前記生体によって前記M個の送信信号が反射、透過または散乱された信号を含む受信信号を、前記N個の受信部のそれぞれが有する受信アンテナ素子を用いて、それぞれが所定期間受信するN個の受信部と、回路と、を備え、前記回路は、前記M個の送信アンテナ素子が送信したM個の前記送信信号にそれぞれが対応するMセットあるN個の前記受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づいて、大きい値ほど、振幅が大きく、かつ、より規則的な波形であることを示す特徴量を算出し、前記算出により得られたM個の前記特徴量を互いに比較することで、前記M個の特徴量のうちで最大となる第1特徴量に対応する第1送信アンテナ素子を特定し、特定した前記第1送信アンテナ素子に基づく所定の方向を前記生体方向として推定する。   In order to achieve the above object, an estimation apparatus according to the present disclosure is an estimation apparatus that estimates a biological direction that is a direction in which the front of a living body is facing, and is arranged to surround a predetermined range including the living body. Transmitters that transmit M transmission signals to the predetermined range using M (M is 3 or more) transmission antenna elements, and N (N is 2 or more) arranged around the predetermined range Each of the reception units includes a reception signal including a signal obtained by reflecting, transmitting, or scattering the M transmission signals by the living body using a reception antenna element included in each of the N reception units. N reception units for receiving a period, and a circuit, and the circuit includes N reception signals of M sets each corresponding to the M transmission signals transmitted by the M transmission antenna elements. For each set of Based on the N received signals, a feature value indicating that the larger the value is, the larger the amplitude is and the more regular waveform is calculated, and the M feature values obtained by the calculation are calculated. By comparing with each other, a first transmission antenna element corresponding to the maximum first feature quantity among the M feature quantities is identified, and a predetermined direction based on the identified first transmission antenna element is determined as the living body. Estimated as direction.

なお、これらの全般的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なCD−ROMなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。   These general or specific aspects may be realized by a system, a method, an integrated circuit, a computer program, or a recording medium such as a computer-readable CD-ROM. The system, method, integrated circuit, computer program And any combination of recording media.

本開示に係る推定装置によれば、生体の正面が向いている方向である生体方向の推定において、推定した生体方向の分解能及び精度を向上させることができる。   According to the estimation apparatus according to the present disclosure, it is possible to improve the resolution and accuracy of the estimated biological direction in the estimation of the biological direction, which is the direction in which the front of the biological body is facing.

図1は、実施の形態における推定装置の一例を示す構成図である。FIG. 1 is a configuration diagram illustrating an example of an estimation apparatus according to an embodiment. 図2は、実施の形態における回路の機能構成の一例を示すブロック図である。FIG. 2 is a block diagram illustrating an example of a functional configuration of a circuit in the embodiment. 図3は、実施の形態における受信波形の極大値高さの総和の算出方法について説明するための図である。FIG. 3 is a diagram for describing a method of calculating the sum of the maximum heights of the received waveforms in the embodiment. 図4は、図5に示す間取りの部屋で生体の方向推定を行った結果を示す一例である。FIG. 4 is an example showing the result of estimating the direction of the living body in the room layout shown in FIG. 図5は、生体の方向推定を行った部屋の間取りを示す図である。FIG. 5 is a diagram showing a floor plan of a room in which the direction of the living body is estimated. 図6は、図5の状況において生体が東側を向いているときに、テレビに注意が向いているとの判断結果を示す表示の一例である。FIG. 6 is an example of a display indicating a determination result that attention is directed to the television when the living body is facing the east side in the situation of FIG. 図7は、実施の形態における推定装置の動作の一例を示すフローチャートである。FIG. 7 is a flowchart illustrating an example of the operation of the estimation device according to the embodiment. 図8は、実施の形態における推定装置による試験で観測された、正面方向の送信アンテナのチャネル応答の絶対値を示す図である。FIG. 8 is a diagram illustrating the absolute value of the channel response of the transmitting antenna in the front direction, which is observed in the test by the estimation apparatus according to the embodiment. 図9は、実施の形態における推定装置による試験で観測された、側面方向の送信アンテナのチャネル応答の絶対値を示す図である。FIG. 9 is a diagram illustrating the absolute value of the channel response of the transmitting antenna in the lateral direction, which is observed in the test by the estimation apparatus according to the embodiment. 図10は、実施の形態における推定装置による試験で観測された、背面方向の送信アンテナのチャネル応答の絶対値を示す図である。FIG. 10 is a diagram illustrating the absolute value of the channel response of the transmitting antenna in the rear direction, which is observed in the test by the estimation apparatus according to the embodiment. 図11は、8つのチャネル応答の絶対値から特徴量を計算した結果の一例を示す図である。FIG. 11 is a diagram illustrating an example of a result of calculating feature amounts from absolute values of eight channel responses. 図12は、図11の結果を二次関数で近似した例を示す図である。FIG. 12 is a diagram illustrating an example in which the result of FIG. 11 is approximated by a quadratic function.

(本発明の基礎となった知見)
電磁波を利用した生体の方向推定に関する従来技術について、発明者らは詳細な検討を行った。特許文献1の方法では、生体に電磁波を照射してその反射波を受信したとき、受信信号の変化量が閾値以上ならば仰向け、閾値未満ならばうつ伏せという判定を行うことが開示されている。しかしながら、この方法では、生体の正面が向いているおおよその方向は推定可能であるが、介護で必要とされるような離床動作に伴う体の正面の方向の変化といった、細かな体の正面の方向を推定することができない。
(Knowledge that became the basis of the present invention)
The inventors conducted a detailed study on the prior art related to the direction estimation of a living body using electromagnetic waves. In the method of Patent Document 1, it is disclosed that when a reflected wave is received by irradiating a living body with electromagnetic waves, a determination is made that if the amount of change in the received signal is greater than or equal to a threshold value, it is supine, and if it is less than the threshold value, prone. However, in this method, it is possible to estimate the approximate direction in which the front of the living body is facing, but there is a small amount of the front of the body, such as a change in the direction of the front of the body that accompanies getting-off movement as required for nursing care. The direction cannot be estimated.

また、特許文献2には、生体に電磁波を照射したとき、生体の正面での反射波の波形には、生体の他の面での反射波形と比較して振幅が大きく、かつ、規則的であるという知見が開示されている。しかしながら、特許文献2には、振幅が大きく、かつ、規則的であるということを定量的に測定する方法は開示されておらず、生体の正面が向いている詳細な方向を知るためには、ユーザが得られた反射波形を目視し経験を踏まえて生体の正面が向いている詳細な方向を判断する必要がある。   Further, in Patent Document 2, when a living body is irradiated with electromagnetic waves, the waveform of the reflected wave on the front surface of the living body has a larger amplitude and regularity than the reflected waveform on the other surface of the living body. The knowledge that there is is disclosed. However, Patent Document 2 does not disclose a method for quantitatively measuring that the amplitude is large and regular, and in order to know the detailed direction in which the front of the living body is facing, It is necessary to visually determine the reflected waveform obtained by the user and determine the detailed direction in which the front of the living body is facing based on experience.

発明者らは、以上の課題に対して研究を重ねた結果、推定装置により推定された生体の正面が向いている生体方向の分解能及び精度を向上させるには、特許文献1のように1方向から送信信号を送信し、その反射波を受信するだけでなく、生体の周囲に複数のアンテナを設置し様々な方向から送信波を送信し、かつ、様々な方向にて反射波、透過波または散乱波を受信することで生体の特徴をより多く補足する事が重要であると考えた。また、特許文献2に開示されている生体正面による反射波形の特徴、つまり、振幅が大きく、かつ、規則的であるということを、定量的に評価できる特徴量を用いることが重要であると考えた。   As a result of repeated research on the above problems, the inventors have improved the resolution and accuracy in the direction of the living body in which the front of the living body estimated by the estimation device is facing, as in Patent Document 1, one direction In addition to transmitting transmission signals from and receiving the reflected waves, multiple antennas are installed around the living body to transmit transmission waves from various directions, and reflected waves, transmitted waves or We thought that it was important to capture more features of the living body by receiving scattered waves. In addition, it is considered important to use features that can quantitatively evaluate the characteristics of the reflected waveform from the front of the living body disclosed in Patent Document 2, that is, that the amplitude is large and regular. It was.

この結果、発明者らは、以下のような推定装置及び推定方法を発明するに至った。   As a result, the inventors have invented the following estimation apparatus and estimation method.

すなわち、本開示に係る推定装置は、生体の正面が向いている方向である生体方向を推定する推定装置であって、前記生体を含む所定範囲を囲んで配置されるM個(Mは3以上)の送信アンテナ素子を用いて前記所定範囲にM個の送信信号を送信する送信部と、前記所定範囲の周囲を囲んで配置されるN個(Nは2以上)の受信部であって、前記生体によって前記M個の送信信号が反射、透過または散乱された信号を含む受信信号を、前記N個の受信部のそれぞれが有する受信アンテナ素子を用いて、それぞれが所定期間受信するN個の受信部と、回路と、を備え、前記回路は、前記M個の送信アンテナ素子が送信したM個の前記送信信号にそれぞれが対応するMセットあるN個の前記受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づいて、大きい値ほど、振幅が大きく、かつ、より規則的な波形であることを示す特徴量を算出し、前記算出により得られたM個の前記特徴量を互いに比較することで、前記M個の特徴量のうちで最大となる第1特徴量に対応する第1送信アンテナ素子を特定し、特定した前記第1送信アンテナ素子に基づく所定の方向を前記生体方向として推定する。   That is, the estimation apparatus according to the present disclosure is an estimation apparatus that estimates the direction of the living body, which is the direction in which the front of the living body is facing, and is arranged around a predetermined range including the living body (M is 3 or more). ) Transmitting unit that transmits M transmission signals to the predetermined range using the transmitting antenna element of N), and N (N is 2 or more) receiving units arranged around the predetermined range, N reception signals each including a signal obtained by reflecting, transmitting or scattering the M transmission signals by the living body are received for a predetermined period using reception antenna elements included in the N reception units. A reception unit and a circuit, and the circuit includes, for each set of N reception signals, M sets, each of which corresponds to the M transmission signals transmitted by the M transmission antenna elements. N pieces included in the set Based on the signal, the larger the value, the larger the amplitude and the feature amount indicating that the waveform is more regular, and the M feature amounts obtained by the calculation are compared with each other, A first transmission antenna element corresponding to the maximum first feature quantity among the M feature quantities is identified, and a predetermined direction based on the identified first transmission antenna element is estimated as the biological direction.

これによれば、所定範囲の周囲に配置されたM個の送信アンテナ素子から送信されたM個の送信信号に対応するM個の特徴量のうち最大の第1特徴量に対応する送信アンテナ素子に基づく所定の方向を生体方向として判定する。このように、送信アンテナ素子および受信アンテナ素子を所定範囲の周囲に配置することで推定した生体方向の分解能及び精度を向上させることができる。   According to this, the transmission antenna element corresponding to the maximum first feature quantity among the M feature quantities corresponding to the M transmission signals transmitted from the M transmission antenna elements arranged around the predetermined range. The predetermined direction based on is determined as the biological direction. As described above, the resolution and accuracy of the biological direction estimated by arranging the transmitting antenna element and the receiving antenna element around a predetermined range can be improved.

例えば、前記回路は、前記MセットあるN個の受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づく第1受信波形のDC(Direct Current)成分を、前記第1受信波形から除去することにより得られる第2受信波形から前記特徴量を算出してもよい。   For example, for each set of the M received signals of the M sets, the circuit converts a DC (Direct Current) component of a first received waveform based on the N received signals included in the set into the first received waveform. The feature amount may be calculated from a second received waveform obtained by removing from the second received waveform.

これにより、受信信号から、生体の識別には不要な成分であるDC成分を抑制することができ、生体識別を効率よく行うことができる。   Thereby, it is possible to suppress a DC component, which is an unnecessary component for identification of a living body, from the received signal, and to perform living body identification efficiently.

例えば、前記回路は、前記第2受信波形の極大値高さの総和を前記特徴量として算出してもよい。   For example, the circuit may calculate the sum of the maximum heights of the second received waveform as the feature amount.

例えば、前記回路は、前記第2受信波形の時間方向の分散を前記特徴量として算出してもよい。   For example, the circuit may calculate a variance in the time direction of the second received waveform as the feature amount.

例えば、前記回路は、前記第2受信波形を周波数領域変換することで得られる関数の所定の周波数範囲での積分値を前記特徴量として算出してもよい。   For example, the circuit may calculate an integral value in a predetermined frequency range of a function obtained by performing frequency domain conversion on the second received waveform as the feature amount.

例えば、前記回路は、特定された前記第1送信アンテナ素子の前記所定範囲に対する第1方向と、前記第1送信アンテナ素子に対応する前記第1特徴量と、前記M個の送信アンテナ素子のうちで前記第1送信アンテナ素子から所定距離以内の範囲に配置されている2以上の第2送信アンテナ素子の前記所定範囲に対する2以上の第2方向と、前記2以上の第2送信アンテナ素子にそれぞれが対応する2以上の第2特徴量とに対して、上に凸の曲線を有する近似関数であって、前記所定範囲に対する方向と前記特徴量との関係を示す近似関数による近似を行うことで得られた近似曲線が最大値を取る方向を前記生体方向として推定してもよい。   For example, the circuit includes a first direction of the identified first transmission antenna element with respect to the predetermined range, the first feature amount corresponding to the first transmission antenna element, and the M transmission antenna elements. And two or more second directions of the two or more second transmission antenna elements arranged within a predetermined distance from the first transmission antenna element with respect to the predetermined range and the two or more second transmission antenna elements, respectively. An approximation function having an upwardly convex curve with respect to two or more second feature quantities corresponding to is performed by approximation using an approximation function indicating a relationship between the direction with respect to the predetermined range and the feature quantity. The direction in which the obtained approximate curve takes the maximum value may be estimated as the biological direction.

これによれば、最大となる第1特徴量と、第1特徴量に対応する第1送信アンテナ素子の近隣の2以上の第2送信アンテナ素子に対応する2以上の第2特徴量と、第1送信アンテナ素子及び2以上の第2送信アンテナ素子の所定範囲に対する第1方向及び2以上の第2方向を用いて、近似関数による近似を行うことで、さらに細かな分解能での方向推定ができる。   According to this, the maximum first feature value, two or more second feature values corresponding to two or more second transmit antenna elements in the vicinity of the first transmit antenna element corresponding to the first feature value, By performing approximation using an approximation function using a first direction and two or more second directions with respect to a predetermined range of one transmission antenna element and two or more second transmission antenna elements, direction estimation with a finer resolution can be performed. .

なお、これらの全般的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータで読み取り可能なCD−ROMなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。   These general or specific aspects may be realized by a system, a method, an integrated circuit, a computer program, or a computer-readable recording medium such as a CD-ROM. The system, method, integrated circuit, computer You may implement | achieve with arbitrary combinations of a program and a recording medium.

以下、本開示の実施の形態について、図面を用いて詳細に説明する。なお、以下で説明する実施の形態は、いずれも本開示の好ましい一具体例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本開示を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本開示の最上位概念を示す独立請求項に記載されていない構成要素については、より好ましい形態を構成する任意の構成要素として説明される。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。   Hereinafter, embodiments of the present disclosure will be described in detail with reference to the drawings. Note that each of the embodiments described below shows a preferred specific example of the present disclosure. Numerical values, shapes, materials, components, arrangement positions and connection forms of components, steps, order of steps, and the like shown in the following embodiments are merely examples, and are not intended to limit the present disclosure. In addition, among the constituent elements in the following embodiments, constituent elements that are not described in the independent claims indicating the highest concept of the present disclosure are described as optional constituent elements that constitute a more preferable embodiment. 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は、実施の形態における推定装置の構成の一例を示す構成図である。
(Embodiment)
FIG. 1 is a configuration diagram illustrating an example of a configuration of an estimation device according to an embodiment.

図1に示すように、推定装置10は、8個の送受信部20A〜20Hと、回路40とを備える。推定装置10は、ヒト等の生体100を含む所定範囲A1に対して送受信部20A〜20Hより送信信号を送信し、送受信部20A〜20Hにて生体100で反射、透過または散乱された信号を含む受信信号を受信する。推定装置10は、送受信部20A〜20Hにより受信された受信信号を、回路40において処理することで、生体100の正面が向いている方向である生体方向を推定する。   As illustrated in FIG. 1, the estimation device 10 includes eight transmission / reception units 20 </ b> A to 20 </ b> H and a circuit 40. The estimation apparatus 10 transmits transmission signals from the transmission / reception units 20A to 20H to a predetermined range A1 including the living body 100 such as a human, and includes signals reflected, transmitted, or scattered by the living body 100 by the transmission / reception units 20A to 20H. Receive a received signal. The estimation device 10 processes a reception signal received by the transmission / reception units 20 </ b> A to 20 </ b> H in the circuit 40, thereby estimating a living body direction that is a direction in which the front surface of the living body 100 is facing.

8個の送受信部20A〜20Hのそれぞれは、1個のアンテナ素子30A〜30Hを有する。送受信部20A〜20Hは、それぞれ、アンテナ素子30A〜30Hを用いて所定範囲A1に送信信号を送信する。これにより、アンテナ素子30A〜30Hは、具体的には、ヒトなどの生体100に対して、マイクロ波を送信信号として発射する。送受信部20A〜20Hは、アンテナ素子30A〜30Hから、無変調の送信信号を送信してもよいし、変調処理が行われた送信信号を送信してもよい。変調処理が行われる場合、送受信部20A〜20Hは、変調処理を行うための回路を含んでいてもよい。アンテナ素子30A〜30Hは、所定範囲A1の周囲の位置、つまり、所定範囲A1の外側の位置に所定範囲A1を囲んで配置される。例えば、アンテナ素子30A〜30Hは、それぞれ、所定範囲A1に対して、異なる8方向の位置に配置される。異なる8方向とは、例えば、所定範囲A1の中心の位置を中心にして360°を8等分することで定められる方向である。8個の送受信部20A〜20Hは、それぞれ、アンテナ素子30A〜30Hとともに所定範囲A1を囲んで配置されていてもよい。また、8個の送受信部20A〜20Hは、アンテナ素子30A〜30Hとは異なる位置に配置されていてもよい。   Each of the eight transmission / reception units 20A to 20H includes one antenna element 30A to 30H. The transmission / reception units 20A to 20H transmit transmission signals to the predetermined range A1 using the antenna elements 30A to 30H, respectively. Thereby, the antenna elements 30A to 30H specifically emit microwaves as transmission signals to the living body 100 such as a human. The transmission / reception units 20A to 20H may transmit unmodulated transmission signals from the antenna elements 30A to 30H, or may transmit transmission signals subjected to modulation processing. When the modulation process is performed, the transmission / reception units 20A to 20H may include a circuit for performing the modulation process. The antenna elements 30A to 30H are arranged so as to surround the predetermined range A1 at positions around the predetermined range A1, that is, positions outside the predetermined range A1. For example, the antenna elements 30A to 30H are respectively arranged at different eight positions with respect to the predetermined range A1. The different eight directions are, for example, directions determined by dividing 360 ° into eight equal parts around the center position of the predetermined range A1. The eight transmission / reception units 20A to 20H may be arranged so as to surround the predetermined range A1 together with the antenna elements 30A to 30H. Further, the eight transmission / reception units 20A to 20H may be arranged at positions different from the antenna elements 30A to 30H.

なお、所定範囲A1とは、推定装置10が生体100の方向を推定するために予め定められた範囲の空間である。   The predetermined range A1 is a space in a predetermined range for the estimation apparatus 10 to estimate the direction of the living body 100.

送受信部20A〜20Hのそれぞれは、生体100によって送受信部20A〜20Hからの送信信号が反射、透過または散乱された信号を含む受信信号を、それぞれが有するアンテナ素子30A〜30Hを用いて、所定期間受信する。送受信部20A〜20Hのそれぞれは、受信信号を周波数変換し、低周波信号に変換してもよい。また、送受信部20A〜20Hのそれぞれは、受信信号に復調処理を行ってもよい。送受信部20A〜20Hのそれぞれは、周波数変換および/または復調処理することにより得られた信号を回路40に出力する。送受信部20A〜20Hのそれぞれは、受信信号を処理するための回路を含んでいてもよい。   Each of the transmission / reception units 20A to 20H uses the antenna elements 30A to 30H each having a reception signal including a signal in which the transmission signal from the transmission / reception units 20A to 20H is reflected, transmitted, or scattered by the living body 100, for a predetermined period. Receive. Each of the transmission / reception units 20A to 20H may convert the frequency of the received signal into a low frequency signal. In addition, each of the transmission / reception units 20A to 20H may perform demodulation processing on the reception signal. Each of the transmission / reception units 20 </ b> A to 20 </ b> H outputs a signal obtained by frequency conversion and / or demodulation processing to the circuit 40. Each of the transmission / reception units 20A to 20H may include a circuit for processing a reception signal.

また、アンテナ素子30A〜30Hは、送信アンテナ素子及び受信アンテナ素子を兼用する構成であるが、これに限らずに、送信アンテナ素子及び受信アンテナ素子は別体の構成であってもよい。また、送信アンテナ素子及び受信アンテナ素子が別体の構成である場合、送信アンテナ素子の数と受信アンテナ素子の数とは互いに異なっていてもよい。送信アンテナ素子は、8個に限らずに、例えば、M個(Mは3以上)であればよい。受信アンテナ素子は、8個に限らずに、例えば、N個(Nは2以上)であればよい。また、送受信部20A〜20Hは、送信信号を送信する送信部と、受信信号を受信する受信部とが一体の構成であるが、送信部及び受信部が別体の構成であってもよい。   In addition, the antenna elements 30A to 30H are configured to serve as both the transmitting antenna element and the receiving antenna element. However, the configuration is not limited thereto, and the transmitting antenna element and the receiving antenna element may be configured separately. In addition, when the transmission antenna element and the reception antenna element have separate structures, the number of transmission antenna elements and the number of reception antenna elements may be different from each other. The number of transmitting antenna elements is not limited to eight, and may be, for example, M (M is 3 or more). The number of receiving antenna elements is not limited to eight, and may be N (N is 2 or more), for example. In addition, although the transmission / reception units 20A to 20H have a configuration in which a transmission unit that transmits a transmission signal and a reception unit that receives a reception signal are integrated, the transmission unit and the reception unit may be configured separately.

次に、回路40の機能的な構成について説明する。   Next, a functional configuration of the circuit 40 will be described.

図2は、実施の形態における回路の機能的な構成の一例を示すブロック図である。   FIG. 2 is a block diagram illustrating an example of a functional configuration of a circuit in the embodiment.

回路40は、算出部41と、比較部42と、方向推定部43とを有し、推定装置10を動作させる各種処理を実行する。回路40は、例えば、制御プログラムを実行するプロセッサと、当該制御プログラムを実行するときに使用するワークエリアとして用いられる揮発性の記憶領域(主記憶装置)とにより構成される。揮発性の記憶領域は、例えば、RAM(Randdom Access Memory)である。また、回路40は、推定装置10を動作させる各種処理を実行する専用回路により実現されてもよい。つまり、回路40は、ソフトウェア処理を実行する回路として実現されてもよいし、ハードウェア処理を実行する回路として実現されてもよい。   The circuit 40 includes a calculation unit 41, a comparison unit 42, and a direction estimation unit 43, and executes various processes for operating the estimation device 10. The circuit 40 includes, for example, a processor that executes a control program and a volatile storage area (main storage device) that is used as a work area used when the control program is executed. The volatile storage area is, for example, a RAM (Random Access Memory). The circuit 40 may be realized by a dedicated circuit that executes various processes for operating the estimation device 10. That is, the circuit 40 may be realized as a circuit that executes software processing, or may be realized as a circuit that executes hardware processing.

回路40は、送受信部20A〜20Hのそれぞれから取得した信号を揮発性の記憶領域に所定期間、一時的に記憶する。回路40は、当該信号の位相および振幅を揮発性の記憶領域所定期間、一時的に記憶してもよい。なお、回路40は、不揮発性の記憶領域を有していてもよく、当該信号を不揮発性の記憶領域に所定期間、一時的に記憶してもよい。   The circuit 40 temporarily stores a signal acquired from each of the transmission / reception units 20A to 20H in a volatile storage area for a predetermined period. The circuit 40 may temporarily store the phase and amplitude of the signal for a volatile storage area for a predetermined period. Note that the circuit 40 may have a nonvolatile storage area, and the signal may be temporarily stored in the nonvolatile storage area for a predetermined period.

以降、回路40を構成する算出部41、比較部42および方向推定部43による各処理について順に説明する。   Hereinafter, each process performed by the calculation unit 41, the comparison unit 42, and the direction estimation unit 43 included in the circuit 40 will be described in order.

まず、算出部41で行う特徴量の算出について説明する。   First, calculation of the feature amount performed by the calculation unit 41 will be described.

算出部41は、回路40の記憶領域に記憶している受信信号を用いて、受信信号の伝搬チャネルH(t)を算出する。   The calculation unit 41 uses the received signal stored in the storage area of the circuit 40 to calculate the propagation channel H (t) of the received signal.

ここで、M個の送信アンテナ素子とN個の受信アンテナ素子とで構成されるMIMO(Multiple−Input and Multiple−Output)アレーアンテナを生体100の周囲に配置した場合に得られる伝搬チャネルH(t)は、以下の式1で表される。   Here, a propagation channel H (t obtained by arranging a MIMO (Multiple-Input and Multiple-Output) array antenna composed of M transmitting antenna elements and N receiving antenna elements around the living body 100. ) Is represented by Equation 1 below.

Figure 2019211458
Figure 2019211458

式1において、hijはj番目の送信アンテナからi番目の受信アンテナの複素チャネル応答を示し、tは観測時間を示す。 In Equation 1, h ij represents the complex channel response from the j-th transmitting antenna to the i-th receiving antenna, and t represents the observation time.

次に、算出部41は、以下の式2で表される、受信信号の伝搬チャネルから、生体100の識別に不要な成分であるDC成分を除去したDC除去チャネルを算出する。算出部41は、MセットあるN個の受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づく第1受信波形のDC成分を、第1受信波形から除去することにより第2受信波形、つまり、DC除去チャネルを得る。算出部41は、算出したDC除去チャネルを回路40の記憶領域に記憶してもよい。   Next, the calculation unit 41 calculates a DC removal channel obtained by removing a DC component, which is a component unnecessary for identification of the living body 100, from the propagation channel of the received signal represented by the following Expression 2. For each set of M received signals of M sets, the calculation unit 41 removes the DC component of the first received waveform based on the N received signals included in the set from the first received waveform, thereby removing the second received signal. A received waveform, ie, a DC rejection channel is obtained. The calculation unit 41 may store the calculated DC removal channel in the storage area of the circuit 40.

ここで、受信信号のDC除去チャネルc(t)は、以下の式3及び式4に示されるように、伝搬チャネルの各成分から、各成分の測定時間平均で算出されるDC成分が減算されることで算出される。   Here, the DC removal channel c (t) of the received signal is obtained by subtracting the DC component calculated by the measurement time average of each component from each component of the propagation channel, as shown in Equation 3 and Equation 4 below. Is calculated.

Figure 2019211458
Figure 2019211458
Figure 2019211458
Figure 2019211458
Figure 2019211458
Figure 2019211458

ここで、Nはスナップショット数を示し、Fはサンプリング周波数を示し、Tは測定時間を示す。 Here, N indicates the number of snapshots, F s represents the sampling frequency, T is shows the measurement time.

なお、DC成分の除去方法は、式3の右辺で表される方法に限らない。例えば、生体100の居ない無人の所定範囲A1において得た伝搬チャネルを減算することでDC成分を除去してもよい。   The DC component removal method is not limited to the method represented by the right side of Equation 3. For example, the DC component may be removed by subtracting a propagation channel obtained in an unattended predetermined range A1 where the living body 100 is not present.

算出部41は、次にアンテナ素子30A〜30H毎のN個のDC除去チャネルをもとに、送信アンテナ素子毎に、当該送信アンテナに対応する第2受信波形について、大きい値ほど、振幅が大きく、かつ、より規則的な波形であることを示す特徴量を算出する。算出部41が算出する特徴量は、例えば、受信波形の極大値高さの総和、受信波形の絶対値の和、受信波形の分散、及び、受信波形を周波数領域変換して得られる関数の所定の周波数範囲での積分値の少なくともいずれか1つである。   The calculation unit 41 then increases the amplitude of the second reception waveform corresponding to the transmission antenna for each transmission antenna element based on the N DC removal channels for each of the antenna elements 30A to 30H as the value increases. In addition, a feature amount indicating a more regular waveform is calculated. The feature amount calculated by the calculation unit 41 is, for example, a predetermined function that is obtained by summing the maximum values of the received waveform, the sum of the absolute values of the received waveform, the variance of the received waveform, and the frequency domain conversion of the received waveform. At least one of the integral values in the frequency range.

まず、算出部41が特徴量として受信波形の極大値高さの総和を算出する場合について説明する。   First, the case where the calculation unit 41 calculates the sum of the maximum heights of the received waveforms as the feature amount will be described.

図3は、受信波形の極大値高さの総和の算出方法について説明するための図である。   FIG. 3 is a diagram for explaining a method for calculating the sum of the maximum heights of the received waveforms.

はじめに、算出部41は、対象となる波形、ここでは、第2受信波形であるの極大値を検出する。図3では、極大値は、三角の点P2で示されている。次に、点P2から水平に線を引き、波形と交わるか、信号の左端または右端に到達するまで左右に伸ばすことで、左右それぞれについて、2つの区間R1、R2が特定される。そして、2つの区間R1、R2のそれぞれの区間内での信号の最小値である点P1、P2が特定される。そして、二つの間隔の最小値である点P1、P3のうち高い方の点P3を基準レベルとし、基準レベルから極大値までの垂直距離を極大値高さとして決定する。この極大値高さを、送信アンテナ素子に対応するDC除去チャネルの受信波形の全ての極大値について算出し、全極大値の極大値高さの総和を計算することで、j番目の送信アンテナ素子に対応する特徴量s(j)が得られる。ここで、算出部41は、雑音による受信波形の乱れによる影響を緩和するため、極大値高さが閾値以下の極大値を極大値高さの総和に含めずに特徴量を算出してもよい。   First, the calculation unit 41 detects the maximum value of the target waveform, here, the second received waveform. In FIG. 3, the maximum value is indicated by a triangular point P2. Next, two sections R1 and R2 are specified for each of the left and right by drawing a line horizontally from the point P2 and extending it to the left and right until it intersects with the waveform or reaches the left or right end of the signal. Then, the points P1 and P2, which are the minimum values of the signals in each of the two sections R1 and R2, are specified. Then, the higher point P3 of the points P1 and P3 which are the minimum values of the two intervals is set as the reference level, and the vertical distance from the reference level to the maximum value is determined as the maximum value height. The maximum value height is calculated for all the maximum values of the received waveform of the DC removal channel corresponding to the transmission antenna element, and the sum of the maximum value heights of all the maximum values is calculated. A feature quantity s (j) corresponding to is obtained. Here, the calculation unit 41 may calculate the feature amount without including the maximum value whose maximum value height is equal to or smaller than the threshold value in the sum of the maximum value heights in order to reduce the influence of the disturbance of the reception waveform due to noise. .

また、算出部41は、N個のDC除去チャネルの絶対値の和を特徴量s(j)として算出してもよい。この場合、算出部41は、以下の式5を用いて特徴量s(j)を算出する。   Further, the calculation unit 41 may calculate the sum of absolute values of the N DC removal channels as the feature amount s (j). In this case, the calculation unit 41 calculates the feature quantity s (j) using the following Expression 5.

Figure 2019211458
Figure 2019211458

また、算出部41は、M個の送信アンテナ素子ごとの伝搬チャネルのDC除去チャネルの受信波形の分散を特徴量s(j)として算出してもよい。この場合、算出部41は、以下の式6を用いて特徴量s(j)を算出する。   Further, the calculation unit 41 may calculate the variance of the received waveform of the DC removal channel of the propagation channel for each of the M transmission antenna elements as the feature quantity s (j). In this case, the calculation unit 41 calculates the feature quantity s (j) using the following Expression 6.

Figure 2019211458
Figure 2019211458

ここで、var_kは、変数kに対する分散を示す。   Here, var_k indicates the variance for the variable k.

また、算出部41は、M個の送信アンテナ素子毎の伝搬チャネルのDC除去チャネルの受信波形を周波数領域変換して得られる関数Fの所定の周波数範囲での積分値を特徴量s(j)として算出してもよい。この場合、算出部41は、以下の式7を用いて特徴量s(j)を算出する。   Further, the calculation unit 41 calculates an integral value in a predetermined frequency range of the function F obtained by performing frequency domain conversion on the reception waveform of the DC removal channel of the propagation channel for each of the M transmission antenna elements as a feature quantity s (j). May be calculated as In this case, the calculation unit 41 calculates the feature quantity s (j) using the following Expression 7.

Figure 2019211458
Figure 2019211458

算出部41はこの3種類の特徴量のうち、少なくとも1つの特徴量を算出し、算出した少なくとも1つの特徴量を比較部42へ出力する。   The calculation unit 41 calculates at least one feature amount among the three types of feature amounts, and outputs the calculated at least one feature amount to the comparison unit 42.

比較部42は、算出部41での算出により得られたM個の送信アンテナ素子に対応するM個の特徴量を互いに比較することで、M個の特徴量のうち最大となる第1特徴量に対応する第1送信アンテナ素子を特定する。   The comparison unit 42 compares the M feature amounts corresponding to the M transmission antenna elements obtained by the calculation in the calculation unit 41 with each other, thereby making the first feature amount that is the maximum among the M feature amounts. The first transmitting antenna element corresponding to is specified.

方向推定部43は、比較部42で特定された第1送信アンテナ素子に基づく所定の方向を生体100の正面が向いている方向である生体方向として推定する。最も単純な例では、方向推定部43は、M個の特徴量のうちで最大である第1特徴量に対応する第1送信アンテナ素子が位置する所定範囲A1に対する方向を生体方向として推定し、推定した方向を推定結果として出力する。   The direction estimation unit 43 estimates a predetermined direction based on the first transmission antenna element specified by the comparison unit 42 as a biological direction that is a direction in which the front of the living body 100 is facing. In the simplest example, the direction estimation unit 43 estimates the direction with respect to the predetermined range A1 where the first transmission antenna element corresponding to the first feature amount that is the largest among the M feature amounts is the biological direction, The estimated direction is output as an estimation result.

また、方向推定部43は、次の方法を用いることで、送信アンテナ素子が位置する所定範囲A1に対する方向以外の方向、例えば、2つの送信アンテナ素子の間の位置への方向を生体方向として推定してもよい。具体的には、方向推定部43は、比較部42において特定された第1送信アンテナ素子の所定範囲A1に対する第1方向と、第1送信アンテナ素子に対応する第1特徴量と、M個の送信アンテナ素子のうちで第1送信アンテナ素子から所定距離以内の範囲に配置されている2以上の第2送信アンテナ素子の所定範囲A1に対する2以上の第2方向と、2以上の第2送信アンテナ素子にそれぞれが対応する2以上の第2特徴量とに対して、近似関数を適用することで得られた近似曲線が最大値を取る方向を生体方向として推定してもよい。なお、方向推定部43は、例えば、2以上の第2送信アンテナ素子として、第1送信アンテナ素子に隣接して配置される、つまり、両隣に配置される2つの送信アンテナ素子を特定してもよい。また、近似関数は、上に凸の曲線を有する近似関数であって、所定範囲A1に対する方向と特徴量との関係を示す近似関数である。   In addition, the direction estimation unit 43 estimates the direction other than the direction with respect to the predetermined range A1 where the transmission antenna element is located, for example, the direction to the position between the two transmission antenna elements as the biological direction by using the following method. May be. Specifically, the direction estimation unit 43 includes a first direction with respect to the predetermined range A1 of the first transmission antenna element specified by the comparison unit 42, a first feature amount corresponding to the first transmission antenna element, and M number of pieces. Among the transmission antenna elements, two or more second directions with respect to a predetermined range A1 of two or more second transmission antenna elements arranged within a predetermined distance from the first transmission antenna element, and two or more second transmission antennas The direction in which the approximate curve obtained by applying the approximate function with respect to two or more second feature quantities each corresponding to the element takes the maximum value may be estimated as the biological direction. Note that the direction estimation unit 43 may identify, for example, two transmission antenna elements that are arranged adjacent to the first transmission antenna element as two or more second transmission antenna elements. Good. The approximate function is an approximate function having an upwardly convex curve, and is an approximate function indicating the relationship between the direction and the feature amount with respect to the predetermined range A1.

このように、方向推定部43では、近似関数が最大値をとる方向を生体方向として推定するため、送信アンテナ素子の数を増やすことなく生体方向の分解能を向上させることができる。また、方向推定部43は、所定の回数分の推定結果を記憶しておき、その平均値を最終的な生体方向とすることで、受信波形に含まれる雑音の影響を軽減することができる。   Thus, since the direction estimation unit 43 estimates the direction in which the approximate function takes the maximum value as the biological direction, the biological direction resolution can be improved without increasing the number of transmission antenna elements. In addition, the direction estimation unit 43 stores the estimation results for a predetermined number of times and sets the average value as the final living body direction, thereby reducing the influence of noise included in the received waveform.

なお、推定装置10は、さらに、推定された生体の方向を表示する表示部50を備えていてもよい。図4は図5に示す間取りの部屋60で生体の方向推定を行った結果を示す一例である。部屋60は、所定範囲A1の一例である。この例に示すように生体の向きを表示する方法には例えば生体の正面の向いている方角を表示する方法や、生体を示すアイコンを回転させる方法などがある。   The estimation device 10 may further include a display unit 50 that displays the estimated direction of the living body. FIG. 4 is an example showing the result of estimating the direction of the living body in the room 60 shown in FIG. The room 60 is an example of the predetermined range A1. As shown in this example, there are a method of displaying the direction of the living body, for example, a method of displaying the direction of the front of the living body, a method of rotating an icon indicating the living body, and the like.

また、生体の正面が向いている方向はその生体の注意が向いている、または生体の意図の向いている方向である可能性が高い。そのため、前記表示部50は、前述した生体の方向を表示する代わりに、注意を向けているものや意図の向いている対象を表示してもよい。図6は図5の状況において生体が東側を向いているときに、テレビに注意が向いているとの判断結果を示す表示の一例である。この場合、推定装置10では、部屋60において、方向と、部屋60の当該方向側に配置されている物体の種類とが関連付けられた情報を記憶しており、当該情報を用いて、推定した生体が推定した方向から生体が向いている方向に対応する物体の種類を、生体の注意の向いている対象として表示部50に表示させる。   In addition, the direction in which the front of the living body is facing is likely to be the direction in which the living body is paying attention or the direction of the living body's intention. For this reason, the display unit 50 may display a target with attention or a target with an intention instead of displaying the direction of the living body described above. FIG. 6 is an example of a display showing a determination result that attention is directed to the television when the living body is facing the east side in the situation of FIG. In this case, the estimation apparatus 10 stores information in which the direction in the room 60 is associated with the type of the object arranged on the direction side of the room 60, and the estimated living body is stored using the information. The type of the object corresponding to the direction in which the living body is facing from the estimated direction is displayed on the display unit 50 as a target to which the living body is paying attention.

また、表示部50が生体の注意が向いている対象を表示する代わりに、機器または照明の遠隔制御を行ってもよい。例えば、生体の注意の向いている方向の照明を点灯して良好な視界を提供することができたり、音声入力による操作が可能な機器が複数存在する場合に、ユーザがどの機器を操作しようとしているかを判断し、適切な機器のみを操作したりすることができる。   Moreover, you may perform remote control of an apparatus or an illumination instead of displaying the object to which the display part 50 is paying attention of a biological body. For example, if there are multiple devices that can provide good visibility by turning on the illumination in the direction in which the biological attention is directed, or that can be operated by voice input, which device the user tries to operate It is possible to operate only appropriate devices.

図7は、実施の形態における推定装置10の動作の一例を示すフローチャートである。   FIG. 7 is a flowchart illustrating an example of the operation of the estimation apparatus 10 according to the embodiment.

推定装置10では、生体100を所定範囲A1内に配置した状態で、アンテナ素子30A〜30Hが所定範囲に送信信号を送信する(S11)。つまり、推定装置10では、8個の送受信部20A〜20Hにより各アンテナ素子30A〜30Hを用いて8個の送信信号を、生体100を含む所定範囲A1に送信する。   In the estimation device 10, the antenna elements 30A to 30H transmit transmission signals to the predetermined range in a state where the living body 100 is disposed within the predetermined range A1 (S11). That is, in the estimation apparatus 10, the eight transmission signals are transmitted to the predetermined range A <b> 1 including the living body 100 using the antenna elements 30 </ b> A to 30 </ b> H by the eight transmission / reception units 20 </ b> A to 20 </ b> H.

次に、送受信部20A〜20Hは、8個のアンテナ素子30A〜30Hを用いて、生体100によって送信信号が反射、透過または散乱された信号を含む受信信号を所定期間受信する(S12)。   Next, the transmission / reception units 20A to 20H receive the reception signal including the signal in which the transmission signal is reflected, transmitted, or scattered by the living body 100 using the eight antenna elements 30A to 30H for a predetermined period (S12).

次に、算出部41は、送受信部20A〜20Hにより受信された受信信号に基づく第1受信波形のDC成分を、第1受信波形から除去する(S13)。これにより、8個の送信アンテナ素子が送信した8個の送信信号にそれぞれが対応する8セットある8個の受信信号のセット毎の8個の第1受信波形からDC成分が除去されるため、8個のセット毎に8個の第2受信波形が得られる。   Next, the calculation unit 41 removes the DC component of the first reception waveform based on the reception signals received by the transmission / reception units 20A to 20H from the first reception waveform (S13). As a result, the DC component is removed from the eight first received waveforms for each set of eight received signals, each of which has eight sets corresponding to the eight transmitted signals transmitted by the eight transmitting antenna elements. Eight second received waveforms are obtained for every eight sets.

算出部41は、8個のセット毎に得られた8個の第2受信波形から、特徴量を算出する(S14)。算出部41は、受信波形の極大値高さの総和、受信波形の絶対値の和、受信波形の分散、及び、受信波形を周波数領域変換して得られる関数の所定の周波数範囲での積分値の少なくとも1つを特徴量として算出する。   The calculating unit 41 calculates a feature amount from the eight second received waveforms obtained for each of the eight sets (S14). The calculation unit 41 calculates the sum of the maximum heights of the received waveform, the sum of the absolute values of the received waveform, the variance of the received waveform, and the integral value in a predetermined frequency range of a function obtained by frequency domain conversion of the received waveform. Is calculated as a feature amount.

比較部42は、算出部41での算出により得られた8個の特徴量を互いに比較することで、8個の特徴量のうちで最大となる第1特徴量に対応する第1送信アンテナ素子を特定する(S15)。   The comparison unit 42 compares the eight feature amounts obtained by the calculation in the calculation unit 41 with each other, thereby the first transmission antenna element corresponding to the first feature amount that is the largest among the eight feature amounts. Is specified (S15).

方向推定部43は、比較部で特定された第1送信アンテナ素子に基づく所定の方向を生体方向として推定する(S16)。   The direction estimation unit 43 estimates a predetermined direction based on the first transmission antenna element specified by the comparison unit as a biological direction (S16).

発明者らは推定装置10の検証のために試験を行った。試験では、送受信部20A〜20Hに相当する送受信機を8台用いた。この8台の送受信機は、生体100を中心として0.5mの半径、かつ、45度間隔で円形に並べて配置されている。ここで8個の送受信機のそれぞれが有する送信アンテナ素子は、1素子から構成され、方形パッチアンテナである。また、8台の送受信機のそれぞれが有する受信アンテナ素子は、1素子から構成され、方形パッチアンテナである。また、床面から受信アンテナ素子が設置される位置までの高さは0.86mである。送信アンテナ素子は、受信アンテナ素子のマイクロ波の1波長真上に配置されている。   The inventors conducted a test to verify the estimation device 10. In the test, eight transceivers corresponding to the transceivers 20A to 20H were used. The eight transceivers are arranged in a circle with a radius of 0.5 m centering on the living body 100 and at intervals of 45 degrees. Here, each of the eight transceivers has a transmitting antenna element that is composed of one element and is a rectangular patch antenna. Each of the eight transceivers has a receiving antenna element that is composed of one element and is a rectangular patch antenna. The height from the floor surface to the position where the receiving antenna element is installed is 0.86 m. The transmitting antenna element is arranged immediately above one wavelength of the microwave of the receiving antenna element.

ここで、図8は、生体100の正面方向に位置する送信アンテナ素子が送信信号を送信した場合のチャネル応答の絶対値|H(t)|を示す。図9は、生体100の側面方向に位置する送信アンテナ素子が送信信号を送信した場合のチャネル応答の絶対値|H(t)|を示す。図10は、生体100の背面方向に位置する送信アンテナ素子が送信信号を送信した場合のチャネル応答の絶対値|H(t)|を示す。 Here, FIG. 8 shows the absolute value | H 1 (t) | of the channel response when the transmission antenna element located in the front direction of the living body 100 transmits a transmission signal. FIG. 9 shows the absolute value | H 3 (t) | of the channel response when the transmission antenna element located in the lateral direction of the living body 100 transmits a transmission signal. FIG. 10 shows the absolute value | H 5 (t) | of the channel response when the transmission antenna element located in the back direction of the living body 100 transmits a transmission signal.

図8〜図10を比較すると、正面から送信している|H(t)|のみ、振幅が大きく、かつ、規則的な波形が観測されていることが分かる。これらの8つのチャネル応答の絶対値から特徴量を計算した結果の一例を図11に示す。この場合、特徴量が最大となる第1送信アンテナ素子は2番であり、生体の正面は、45°の方向を向いていると判定される。つまり、比較部42は、第1送信アンテナ素子として2番の送信アンテナ素子を特定する。また、特徴量の最大値をとる送信アンテナ及びその近隣の送信アンテナの特徴量を用いて関数による近似を行うことでさらに細かな分解能での方向推定ができる。 8 to 10, it can be seen that only | H 1 (t) | transmitted from the front has a large amplitude and a regular waveform is observed. An example of the result of calculating the feature amount from the absolute values of these eight channel responses is shown in FIG. In this case, it is determined that the first transmitting antenna element having the maximum feature amount is No. 2, and the front of the living body is oriented in a 45 ° direction. That is, the comparison unit 42 specifies the second transmission antenna element as the first transmission antenna element. In addition, direction approximation with a finer resolution can be performed by performing approximation using a function using the feature values of the transmission antenna having the maximum feature value and the neighboring transmission antennas.

図12は、図11の結果を二次関数で近似した例を示す図である。   FIG. 12 is a diagram illustrating an example in which the result of FIG. 11 is approximated by a quadratic function.

この場合、方向推定部43は、第1送信アンテナ素子として特定された2番の送信アンテナ素子に対応する方向である45°と、その特徴量0.88とを抽出する。また、方向推定部43は、2番の送信アンテナ素子に隣接する1番の送信アンテナ素子に対応する方向である0°と、その特徴量0.79とを抽出する。また、方向推定部は、2番の送信アンテナ素子に隣接する3番の送信アンテナ素子に対応する方向である90°と、その特徴量0.15とを抽出する。つまり、方向推定部43は、3つの(「送信アンテナ素子に対する方向」、「特徴量」)の組み合わせである、P11(0°、0.79)、P12(45°、0.88)、及び、P13(90°、0.15)の3点を二次関数で近似することにより得られる曲線の極大値を示すP14に対応する方向、例えば30°を生体方向として推定することができる。   In this case, the direction estimation unit 43 extracts 45 ° that is the direction corresponding to the second transmission antenna element specified as the first transmission antenna element and the feature amount 0.88. In addition, the direction estimation unit 43 extracts 0 ° that is the direction corresponding to the first transmission antenna element adjacent to the second transmission antenna element and the feature amount 0.79. In addition, the direction estimation unit extracts 90 ° that is the direction corresponding to the third transmission antenna element adjacent to the second transmission antenna element and the feature amount 0.15. That is, the direction estimation unit 43 is a combination of three (“direction with respect to the transmission antenna element” and “feature”), P11 (0 °, 0.79), P12 (45 °, 0.88), and , P13 (90 °, 0.15), a direction corresponding to P14 indicating the maximum value of the curve obtained by approximating the three points with a quadratic function, for example, 30 ° can be estimated as the biological direction.

(効果など)
本実施の形態に係る推定装置10によれば、M個の送信アンテナ素子から送信されたM個の送信信号に対応するM個の特徴量のうち最大の第1特徴量に対応する送信アンテナ素子に基づく所定の方向を生体方向として判定する。このように、送信アンテナ素子および受信アンテナ素子を所定範囲の周囲に配置することで推定した生体方向の分解能及び精度を向上させることができる。
(Effect etc.)
According to the estimation apparatus 10 according to the present embodiment, the transmission antenna element corresponding to the maximum first feature quantity among the M feature quantities corresponding to the M transmission signals transmitted from the M transmission antenna elements. The predetermined direction based on is determined as the biological direction. As described above, the resolution and accuracy of the biological direction estimated by arranging the transmitting antenna element and the receiving antenna element around a predetermined range can be improved.

また、本実施の形態に係る推定装置10によれば、最大となる第1特徴量と、第1特徴量に対応する第1送信アンテナ素子の近隣の2以上の第2送信アンテナ素子に対応する2以上の第2特徴量と、第1送信アンテナ素子及び2以上の第2送信アンテナ素子の所定範囲に対する第1方向及び2以上の第2方向を用いて、近似関数による近似を行うことで得られる所定の方向を生体方向として推定するため、さらに細かな分解能での方向推定ができる。   Moreover, according to the estimation apparatus 10 according to the present embodiment, it corresponds to the maximum first feature quantity and two or more second transmission antenna elements in the vicinity of the first transmission antenna element corresponding to the first feature quantity. It is obtained by performing approximation using an approximation function using two or more second feature values, a first direction and two or more second directions with respect to a predetermined range of the first transmission antenna element and the two or more second transmission antenna elements. Since the predetermined direction is estimated as the biological direction, the direction can be estimated with a finer resolution.

また、本実施の形態に係る推定装置10は、マイクロ波などの無線信号を用いて、ヒト等の生体100の正面が向いている方向を推定することができる。つまり、本実施の形態に係る推定装置10は、カメラ等で撮像した画像に対して画像解析することなくヒト等の生体100の生体方向を推定できるため、ヒトのプライバシーを保護した状態で、ヒトの生体方向の推定を行うことができる。   Moreover, the estimation apparatus 10 according to the present embodiment can estimate the direction in which the front surface of the living body 100 such as a human is facing using a radio signal such as a microwave. That is, since the estimation apparatus 10 according to the present embodiment can estimate the living body direction of the living body 100 such as a human without performing image analysis on an image captured by a camera or the like, Can be estimated.

上記実施の形態に係る推定装置10の、送信部の機能と、受信部の機能とを入れ替えた装置を推定装置として採用してもよい。   You may employ | adopt the apparatus which replaced the function of the transmission part and the function of the receiving part of the estimation apparatus 10 which concerns on the said embodiment as an estimation apparatus.

本発明は、無線信号を利用した生体の方向を推定する推定装置に利用でき、特に、生体の方向に応じた制御を行う家電機器、高齢者の起き上がりや寝返りを検知して見守りを行うシステムなどに利用できる。   INDUSTRIAL APPLICABILITY The present invention can be used in an estimation device that estimates the direction of a living body using a radio signal, and in particular, home appliances that perform control according to the direction of the living body, a system that detects and monitors the rising or turning of an elderly person, and the like Available to:

10 推定装置
20A〜20H 送受信部
30A〜30H アンテナ素子
40 回路
41 算出部
42 比較部
43 方向推定部
50 表示部
100 生体
DESCRIPTION OF SYMBOLS 10 Estimation apparatus 20A-20H Transmission / reception part 30A-30H Antenna element 40 Circuit 41 Calculation part 42 Comparison part 43 Direction estimation part 50 Display part 100 Living body

Claims (7)

生体の正面が向いている方向である生体方向を推定する推定装置であって、
前記生体を含む所定範囲を囲んで配置されるM個(Mは3以上)の送信アンテナ素子を用いて前記所定範囲にM個の送信信号を送信する送信部と、
前記所定範囲の周囲を囲んで配置されるN個(Nは2以上)の受信部であって、前記生体によって前記M個の送信信号が反射、透過または散乱された信号を含む受信信号を、前記N個の受信部のそれぞれが有する受信アンテナ素子を用いて、それぞれが所定期間受信するN個の受信部と、
回路と、を備え、
前記回路は、
前記M個の送信アンテナ素子が送信したM個の前記送信信号にそれぞれが対応するMセットあるN個の前記受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づいて、大きい値ほど、振幅が大きく、かつ、より規則的な波形であることを示す特徴量を算出し、
前記算出により得られたM個の前記特徴量を互いに比較することで、前記M個の特徴量のうちで最大となる第1特徴量に対応する第1送信アンテナ素子を特定し、
特定した前記第1送信アンテナ素子に基づく所定の方向を前記生体方向として推定する
推定装置。
An estimation device for estimating a direction of a living body, which is a direction in which the front of the living body is facing,
A transmission unit that transmits M transmission signals to the predetermined range using M (M is 3 or more) transmission antenna elements arranged around the predetermined range including the living body;
N reception units (N is 2 or more) arranged around the predetermined range, and a reception signal including a signal in which the M transmission signals are reflected, transmitted, or scattered by the living body, N receiving units each receiving for a predetermined period using a receiving antenna element included in each of the N receiving units;
A circuit,
The circuit is
For each set of N received signals that are M sets corresponding to the M transmit signals transmitted by the M transmit antenna elements, the set is large based on the N received signals included in the set. Calculate the feature value that indicates that the value is larger and the waveform is more regular,
By comparing the M feature quantities obtained by the calculation with each other, the first transmission antenna element corresponding to the first feature quantity that is the maximum among the M feature quantities is specified,
An estimation device that estimates a predetermined direction based on the identified first transmission antenna element as the biological direction.
前記回路は、前記MセットあるN個の受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づく第1受信波形のDC(Direct Current)成分を、前記第1受信波形から除去することにより得られる第2受信波形から前記特徴量を算出する
請求項1に記載の推定装置。
The circuit removes, from the first received waveform, a DC (Direct Current) component of the first received waveform based on the N received signals included in the set for each set of the M received signals of the M sets. The estimation device according to claim 1, wherein the feature amount is calculated from a second received waveform obtained by performing the processing.
前記回路は、前記第2受信波形の極大値高さの総和を前記特徴量として算出する
請求項2に記載の推定装置。
The estimation device according to claim 2, wherein the circuit calculates a sum of maximum heights of the second received waveform as the feature amount.
前記回路は、前記第2受信波形の時間方向の分散を前記特徴量として算出する
請求項2に記載の推定装置。
The estimation device according to claim 2, wherein the circuit calculates a variance in the time direction of the second received waveform as the feature amount.
前記回路は、前記第2受信波形を周波数領域変換することで得られる関数の所定の周波数範囲での積分値を前記特徴量として算出する
請求項2に記載の推定装置。
The estimation device according to claim 2, wherein the circuit calculates, as the feature amount, an integral value in a predetermined frequency range of a function obtained by performing frequency domain conversion on the second received waveform.
前記回路は、特定された前記第1送信アンテナ素子の前記所定範囲に対する第1方向と、前記第1送信アンテナ素子に対応する前記第1特徴量と、前記M個の送信アンテナ素子のうちで前記第1送信アンテナ素子から所定距離以内の範囲に配置されている2以上の第2送信アンテナ素子の前記所定範囲に対する2以上の第2方向と、前記2以上の第2送信アンテナ素子にそれぞれが対応する2以上の第2特徴量とに対して、上に凸の曲線を有する近似関数であって、前記所定範囲に対する方向と前記特徴量との関係を示す近似関数による近似を行うことで得られた近似曲線が最大値を取る方向を前記生体方向として推定する
請求項1から5のいずれか1項に記載の推定装置。
The circuit includes a first direction of the identified first transmission antenna element with respect to the predetermined range, the first feature amount corresponding to the first transmission antenna element, and the M transmission antenna elements. Two or more second transmitting antenna elements arranged within a predetermined distance from the first transmitting antenna element correspond to two or more second directions with respect to the predetermined range, and each of the two or more second transmitting antenna elements. An approximation function having an upwardly convex curve with respect to two or more second feature quantities to be obtained, and obtained by performing approximation using an approximation function indicating the relationship between the direction with respect to the predetermined range and the feature quantity. The estimation apparatus according to any one of claims 1 to 5, wherein a direction in which the approximate curve takes a maximum value is estimated as the biological direction.
生体の正面が向いている方向である生体方向を推定する推定装置が実行する推定方法であって、
前記推定装置は、
前記生体を含む所定範囲を囲んで配置されるM個(Mは3以上)の送信アンテナ素子を用いて前記所定範囲にM個の送信信号を送信する送信部と、前記所定範囲の周囲を囲んで配置されるN個(Nは2以上)の受信部と、回路と、を備え、
前記推定方法では、
前記M個の送信アンテナ素子を用いて、前記生体を含む前記所定範囲に前記送信信号を送信し、
前記N個の受信部が有する前記N個の受信アンテナ素子を用いて、前記生体によって前記送信信号が反射、透過または散乱された信号を含む受信信号を所定期間受信し、
前記M個の送信アンテナ素子が送信したM個の前記送信信号にそれぞれが対応するMセットあるN個の前記受信信号のセット毎に、当該セットに含まれるN個の受信信号に基づいて、大きい値ほど、振幅が大きく、かつ、より規則的な波形であることを示す特徴量を算出し、
前記算出により得られたM個の前記特徴量を互いに比較することで、前記M個の特徴量のうちで最大となる第1特徴量に対応する第1送信アンテナ素子を特定し、
特定した前記第1送信アンテナ素子に基づく所定の方向を前記生体方向として推定する
推定方法。
An estimation method executed by an estimation device that estimates a biological direction that is a direction in which the front of the biological body is facing,
The estimation device includes:
A transmission unit that transmits M transmission signals to the predetermined range using M (M is 3 or more) transmission antenna elements arranged around the predetermined range including the living body, and surrounds the periphery of the predetermined range N receivers (N is 2 or more) arranged in the circuit and a circuit,
In the estimation method,
Using the M transmitting antenna elements, the transmission signal is transmitted to the predetermined range including the living body,
Using the N reception antenna elements of the N reception units, a reception signal including a signal in which the transmission signal is reflected, transmitted, or scattered by the living body is received for a predetermined period.
For each set of N received signals that are M sets corresponding to the M transmit signals transmitted by the M transmit antenna elements, the set is large based on the N received signals included in the set. Calculate the feature value that indicates that the value is larger and the waveform is more regular,
By comparing the M feature quantities obtained by the calculation with each other, the first transmission antenna element corresponding to the first feature quantity that is the maximum among the M feature quantities is identified,
An estimation method for estimating a predetermined direction based on the identified first transmission antenna element as the biological direction.
JP2018160851A 2018-01-09 2018-08-29 Estimation device and estimation method Active JP7129648B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910011041.3A CN110018522B (en) 2018-01-09 2019-01-07 Estimation device and estimation method
US16/240,919 US10969479B2 (en) 2018-01-09 2019-01-07 Estimation device and estimation method

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2018001045 2018-01-09
JP2018001045 2018-01-09
JP2018104092 2018-05-30
JP2018104092 2018-05-30

Publications (2)

Publication Number Publication Date
JP2019211458A true JP2019211458A (en) 2019-12-12
JP7129648B2 JP7129648B2 (en) 2022-09-02

Family

ID=68846606

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018160851A Active JP7129648B2 (en) 2018-01-09 2018-08-29 Estimation device and estimation method

Country Status (1)

Country Link
JP (1) JP7129648B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021140988A1 (en) 2020-01-07 2021-07-15 パナソニックIpマネジメント株式会社 Estimation device and estimation method
WO2023276945A1 (en) 2021-06-29 2023-01-05 パナソニックIpマネジメント株式会社 Estimation device and estimation method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005314847A (en) * 2004-04-30 2005-11-10 Shima Seiki Mfg Ltd Method and apparatus for measuring body shape
JP2013235012A (en) * 2013-07-22 2013-11-21 Oki Electric Ind Co Ltd Data processing apparatus, operation recognition system, operation determination method, and program
JP2014228291A (en) * 2013-05-20 2014-12-08 三星電子株式会社Samsung Electronics Co.,Ltd. Radio detection device and radio detection method
US20160320852A1 (en) * 2015-04-30 2016-11-03 Google Inc. Wide-Field Radar-Based Gesture Recognition
JP2018036247A (en) * 2016-08-29 2018-03-08 パナソニックIpマネジメント株式会社 Identification device and identification method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005314847A (en) * 2004-04-30 2005-11-10 Shima Seiki Mfg Ltd Method and apparatus for measuring body shape
JP2014228291A (en) * 2013-05-20 2014-12-08 三星電子株式会社Samsung Electronics Co.,Ltd. Radio detection device and radio detection method
JP2013235012A (en) * 2013-07-22 2013-11-21 Oki Electric Ind Co Ltd Data processing apparatus, operation recognition system, operation determination method, and program
US20160320852A1 (en) * 2015-04-30 2016-11-03 Google Inc. Wide-Field Radar-Based Gesture Recognition
JP2018036247A (en) * 2016-08-29 2018-03-08 パナソニックIpマネジメント株式会社 Identification device and identification method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
袁 之雄 他: "アレーアンテナを用いた生体向き推定法", 電子情報通信学会2018年総合大会講演論文集 通信1, JPN6022030319, 6 March 2018 (2018-03-06), JP, pages 132, ISSN: 0004833611 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021140988A1 (en) 2020-01-07 2021-07-15 パナソニックIpマネジメント株式会社 Estimation device and estimation method
WO2023276945A1 (en) 2021-06-29 2023-01-05 パナソニックIpマネジメント株式会社 Estimation device and estimation method

Also Published As

Publication number Publication date
JP7129648B2 (en) 2022-09-02

Similar Documents

Publication Publication Date Title
US10880301B2 (en) Identification device and identification method
JP6751923B2 (en) program
US10895631B2 (en) Sensor and estimating method
US20140378809A1 (en) Systems and methods for extracting physiological characteristics using frequency harmonics
US20190008459A1 (en) Method and a system for detecting a vital sign of a subject
US8142361B2 (en) Method and apparatus for automatic optimization of doppler imaging parameters
JP2017513656A (en) Monitoring vital signs by radio reflection
US10561358B2 (en) Biometric device and biometric method
US10456115B2 (en) Ultrasound system and clutter filtering method thereof
US11638562B2 (en) Brain activity analysis method and apparatus thereof
CN111759305A (en) Living body detection method and living body detection system
JP2019211458A (en) Estimation device and estimation method
US20170296152A1 (en) Ultrasound coupling medium detection
JP6887089B2 (en) Identification device and identification method
KR20180088019A (en) Method and Device for Measuring Biometric Data using Radar
US20140336510A1 (en) Enhancement in Diagnostic Ultrasound Spectral Doppler Imaging
JP2019197039A (en) Estimation device, living body number estimation device, estimation method, and program
RU2384860C2 (en) Method of detecting people and moving objects behind barrier and device for realising said method
JP7065404B2 (en) Identification device and identification method
CN110693454A (en) Sleep characteristic event detection method and device based on radar and storage medium
US10969479B2 (en) Estimation device and estimation method
US20210321938A1 (en) System and method for wireless detection and measurement of a subject rising from rest
WO2021140988A1 (en) Estimation device and estimation method
WO2023276945A1 (en) Estimation device and estimation method
US20210219958A1 (en) Ultrasound diagnostic apparatus and method of controlling the same, and computer program product

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210721

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220706

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220810

R151 Written notification of patent or utility model registration

Ref document number: 7129648

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151