JP6868846B2 - Sensors and methods - Google Patents

Sensors and methods Download PDF

Info

Publication number
JP6868846B2
JP6868846B2 JP2017180438A JP2017180438A JP6868846B2 JP 6868846 B2 JP6868846 B2 JP 6868846B2 JP 2017180438 A JP2017180438 A JP 2017180438A JP 2017180438 A JP2017180438 A JP 2017180438A JP 6868846 B2 JP6868846 B2 JP 6868846B2
Authority
JP
Japan
Prior art keywords
living body
estimated
antenna elements
transmitting antenna
receiving antenna
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017180438A
Other languages
Japanese (ja)
Other versions
JP2018112540A (en
Inventor
武司 中山
武司 中山
翔一 飯塚
翔一 飯塚
尚樹 本間
尚樹 本間
大 笹川
大 笹川
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 US15/851,858 priority Critical patent/US10912493B2/en
Priority to CN201711454518.2A priority patent/CN108279413B/en
Publication of JP2018112540A publication Critical patent/JP2018112540A/en
Application granted granted Critical
Publication of JP6868846B2 publication Critical patent/JP6868846B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、無線信号を利用して生体の行動を推定するセンサーおよび方法に関する。 The present invention relates to a sensor and a method for estimating the behavior of a living body using a radio signal.

人物の位置、行動などを知る方法として、無線信号を利用する方法が検討されている(例えば、特許文献1〜7参照)。具体的には、特許文献1には、受信電波の変化量により人の動作を監視し、人の在または不在を判断する方法が開示されている。特許文献2には、THz波を用い、生体の頭部および四肢を把握する方法が開示されている。特許文献3には、電波レーダにより対象物の大きさを推定する方法が開示されている。特許文献4には、ミリ波レーダにより対象物の位置の軌跡を測定する方法が開示されている。特許文献5には、ドップラーレーダのRCS測定により対象物が人か否かを判断する方法が開示されている。特許文献6には、複数アンテナのチャネル情報および各種センサ情報をもとに機械学習により生体の位置や状態を推定する方法が開示されている。特許文献7には、FMCWレーダの測定結果より人の臥位または座位を判定する方法が開示されている。 As a method of knowing the position, behavior, etc. of a person, a method of using a wireless signal has been studied (see, for example, Patent Documents 1 to 7). Specifically, Patent Document 1 discloses a method of monitoring the movement of a person based on the amount of change in received radio waves and determining the presence or absence of the person. Patent Document 2 discloses a method of grasping the head and limbs of a living body using THz waves. Patent Document 3 discloses a method of estimating the size of an object by a radio wave radar. Patent Document 4 discloses a method of measuring a trajectory of a position of an object by a millimeter-wave radar. Patent Document 5 discloses a method of determining whether or not an object is a human by RCS measurement of a Doppler radar. Patent Document 6 discloses a method of estimating the position and state of a living body by machine learning based on channel information of a plurality of antennas and various sensor information. Patent Document 7 discloses a method of determining a person's lying position or sitting position based on the measurement result of the FMCW radar.

特開2015−184149号公報Japanese Unexamined Patent Publication No. 2015-184149 特開2006−81771号公報Japanese Unexamined Patent Publication No. 2006-81771 特開2001−159678号公報Japanese Unexamined Patent Publication No. 2001-159678 特開2013−160730号公報Japanese Unexamined Patent Publication No. 2013-160730 特開2004−340729号公報Japanese Unexamined Patent Publication No. 2004-340729 特開2014−190724号公報Japanese Unexamined Patent Publication No. 2014-190724 特開2016−135233号公報Japanese Unexamined Patent Publication No. 2016-135233

しかしながら、無線信号を利用することで、生体の動作を推定する制度を向上させるには、さらなる改善が求められている。 However, further improvement is required to improve the system for estimating the movement of living organisms by using wireless signals.

上記目的を達成するために、本発明の一形態に係るセンサーは、センサーであって、生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは3以上の自然数)の送信アンテナ素子を有する送信アンテナと、前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは3以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、前記センサーに対する前記生体の存在する鉛直方向における位置である鉛直位置、および、RCS(Radar cross-section)値の時間的変化と、前記生体の動作との対応関係を示す情報を記憶しているメモリと、を備え、前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記回路は、前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、前記第2行列を用いて、前記センサーに対する前記生体の存在する3次元位置であって、前記鉛直位置を含む3次元位置を推定し、前記推定した前記3次元位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、前記第1距離および前記第2距離を用いて、前記生体に対するRCS値を算出し、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定する。 In order to achieve the above object, the antenna according to one embodiment of the present invention is a sensor, each of which transmits N transmission signals to a predetermined range in which a living body can exist (N is a natural number of 3 or more). A transmitting antenna having the transmitting antenna element of the above, and N receiving signals including a reflected signal in which a part of the N transmitting signals transmitted by the N transmitting antenna elements is reflected by the living body. A receiving antenna having M receiving antenna elements, each of which receives a signal (M is a natural number of 3 or more), a circuit, a vertical position which is a position in the vertical direction in which the living body exists with respect to the sensor, and RCS ( Radar cross-section) A memory that stores information indicating the correspondence between the temporal change of the value and the operation of the living body is provided, and at least three of the N transmitting antenna elements are transmitted antenna elements. Are arranged at different positions in the vertical and horizontal directions, respectively, and at least three of the M receiving antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively. Is between each of the N transmitting antenna elements and each of the M receiving antenna elements from each of the N receiving signals received in each of the M receiving antenna elements in a predetermined period. By calculating the first matrix of N × M having each complex transmission function showing the propagation characteristics as a component and extracting the second matrix corresponding to the predetermined frequency range in the first matrix, the respiration and heartbeat of the living body. And the second matrix corresponding to the component affected by the vital activity including at least one of the body movements is extracted, and the second matrix is used to be the three-dimensional position where the living body exists with respect to the sensor. , The three-dimensional position including the vertical position is estimated, and the distance between the living body and the transmitting antenna is determined based on the estimated three-dimensional position, the position of the transmitting antenna, and the position of the receiving antenna. The first distance shown and the second distance indicating the distance between the living body and the receiving antenna were calculated, and the RCS value for the living body was calculated and estimated using the first distance and the second distance. The movement of the living body is estimated using the three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory.

本発明によれば、無線信号を利用することで、生体の動作の推定を、短時間かつ高精度に行うことができる。 According to the present invention, by using a wireless signal, it is possible to estimate the movement of a living body in a short time and with high accuracy.

図1は、実施の形態におけるセンサーの構成の一例を示すブロック図である。FIG. 1 is a block diagram showing an example of a sensor configuration according to an embodiment. 図2は、実施の形態におけるセンサーの設置例を示す図である。FIG. 2 is a diagram showing an installation example of the sensor in the embodiment. 図3は実施の形態における回路およびメモリの機能的な構成を示すブロック図である。FIG. 3 is a block diagram showing a functional configuration of a circuit and a memory according to an embodiment. 図4Aは、異なる測定期間で測定するセンサーの位置分解能の一例について説明するための図である。FIG. 4A is a diagram for explaining an example of the position resolution of the sensor measuring in different measurement periods. 図4Bは、異なる測定期間で測定するセンサーの位置分解能の他の一例について説明するための図である。FIG. 4B is a diagram for explaining another example of the position resolution of the sensor measuring in different measurement periods. 図5は、高さ位置(Height)またはRCE値の時系列データから動作期間を抽出する例について説明するための図である。FIG. 5 is a diagram for explaining an example of extracting an operation period from time series data of a height position (Height) or an RCE value. 図6は、方向ベクトルに変換する処理について説明するための図である。FIG. 6 is a diagram for explaining a process of converting into a direction vector. 図7は、方向コード表の一例を示す図である。FIG. 7 is a diagram showing an example of a direction code table. 図8は、算出された方向コードおよび距離の時系列データの一例を示す図である。FIG. 8 is a diagram showing an example of time-series data of the calculated direction code and distance. 図9は、対応関係を示す情報の一例を示す図である。FIG. 9 is a diagram showing an example of information showing a correspondence relationship. 図10は、測定により得られたテストデータと、モデルコードであるモデルデータとを示す図である。FIG. 10 is a diagram showing test data obtained by measurement and model data which is a model code. 図11は、実施の形態におけるセンサーの動作の一例を示すフローチャートである。FIG. 11 is a flowchart showing an example of the operation of the sensor according to the embodiment. 図12は、推定処理の詳細の一例を示すフローチャートである。FIG. 12 is a flowchart showing an example of details of the estimation process. 図13は、事前学習処理におけるセンサーの動作の一例を示すフローチャートである。FIG. 13 is a flowchart showing an example of the operation of the sensor in the pre-learning process.

(本発明の基礎となった知見)
無線信号を利用した生体の状態推定に関する従来技術について、発明者らは詳細な検討を行った。その結果、特許文献1の方法では、人物の在、不在は検知可能だが、人物の存在する方向や位置、状態、動作などを検出することは困難であるという問題があることがわかった。
(Knowledge that became the basis of the present invention)
The inventors have conducted a detailed study on the prior art for estimating the state of a living body using a wireless signal. As a result, it was found that the method of Patent Document 1 can detect the presence or absence of a person, but has a problem that it is difficult to detect the direction, position, state, movement, etc. of the person.

また、特許文献2の方法では、人物の頭部および四肢を検知し、人物の存在する方向、位置、状態などを推定可能だが、テラヘルツ帯デバイスを用いており、コストが高い問題があることがわかった。 Further, in the method of Patent Document 2, it is possible to detect the head and limbs of a person and estimate the direction, position, state, etc. of the person, but there is a problem that the cost is high because a terahertz band device is used. all right.

また、特許文献3の方法では、対象物の大きさを推定することは可能だが、人物などの生体の状態を推定することは困難であるという問題があることがわかった。 Further, it has been found that the method of Patent Document 3 has a problem that it is possible to estimate the size of an object, but it is difficult to estimate the state of a living body such as a person.

また、特許文献4の方法では、対象生体の位置の軌跡を推定することは可能だが、生体の状態を推定することは困難であるという問題があることがわかった。 Further, it has been found that the method of Patent Document 4 has a problem that it is possible to estimate the locus of the position of the target living body, but it is difficult to estimate the state of the living body.

また、特許文献5の方法では、RCSにより対象物が人物か否か推定することは可能だが、人物などの生体の状態を推定することは困難であるという問題があることがわかった。 Further, it has been found that the method of Patent Document 5 has a problem that it is possible to estimate whether or not the object is a person by RCS, but it is difficult to estimate the state of a living body such as a person.

また、特許文献6では、使用者ごとに機械学習をする必要があるという課題があることがわかった。 Further, in Patent Document 6, it was found that there is a problem that machine learning needs to be performed for each user.

また、特許文献7では、人物の動作を推定することは困難であるという問題があることが分かった。 Further, in Patent Document 7, it has been found that there is a problem that it is difficult to estimate the movement of a person.

発明者らは、以上の課題に対して研究を重ねた結果、垂直方向および水平方向に異なる位置に置かれた複数のアンテナ素子を含む送信アンテナから送信され、生体によって反射された反射信号の伝搬特性と散乱断面積とを用いることにより、生体が存在する方向、位置、大きさ、姿勢、動作などの推定を短時間かつ高精度に行うことが可能であることを見出し、本発明に至った。 As a result of repeated research on the above problems, the inventors have propagated a reflected signal transmitted from a transmitting antenna including a plurality of antenna elements placed at different positions in the vertical direction and the horizontal direction and reflected by a living body. We have found that it is possible to estimate the direction, position, size, posture, motion, etc. of a living body in a short time and with high accuracy by using the characteristics and the scattering cross section, and have reached the present invention. ..

更に静止中の生体の位置を推定するには、主に生体の呼吸、心拍成分をもとに位置推定を行っていることから、生体の動作の推定には数秒程度のデータが必要となる。このため、例えば、転倒など動きの早い動作の推定には、1秒未満の短いデータを用いて推定が必要であることがわかり、生体の静止中、生体の動作中、または静止中および動作中の両方における生体の姿勢の推定には、1度の測定データについて、異なる時間データ長を用いて、生体の位置および姿勢の推定を行う必要があることがわかった。 Furthermore, in order to estimate the position of a living body at rest, the position is estimated mainly based on the respiratory and heartbeat components of the living body, so data of about several seconds is required to estimate the movement of the living body. Therefore, for example, it was found that estimation of a fast-moving motion such as a fall requires estimation using short data of less than 1 second, and the living body is stationary, the living body is moving, or is stationary and moving. It was found that it is necessary to estimate the position and posture of the living body by using different time data lengths for one measurement data in order to estimate the posture of the living body in both cases.

すなわち、本発明の一態様に係るセンサーは、センサーであって、生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは3以上の自然数)の送信アンテナ素子を有する送信アンテナと、前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは3以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、前記センサーに対する前記生体の存在する鉛直方向における位置である鉛直位置、および、RCS(Radar cross-section)値の時間的変化と、前記生体の動作との対応関係を示す情報を記憶しているメモリと、を備え、前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記回路は、前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、前記第2行列を用いて、前記センサーに対する前記生体の存在する3次元位置であって、前記鉛直位置を含む3次元位置を推定し、前記推定した前記3次元位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、前記第1距離および前記第2距離を用いて、前記生体に対するRCS値を算出し、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定する。 That is, the sensor according to one aspect of the present invention is a sensor and has N transmitting antenna elements (N is a natural number of 3 or more) in which transmission signals are transmitted to a predetermined range in which a living body can exist. Each of the transmitting antenna and the N transmitting signals transmitted by the N transmitting antenna elements receive N receiving signals including a reflected signal reflected by the living body. A receiving antenna having M (M is a natural number of 3 or more), a circuit, a vertical position which is a position in the vertical direction in which the living body exists with respect to the sensor, and an RCS (Radar cross-section) value. A memory that stores information indicating a correspondence relationship between the temporal change of the living body and the operation of the living body is provided, and at least three of the N transmitting antenna elements are in the vertical direction, respectively. And at different positions in the horizontal direction, at least three of the M receiving antenna elements are arranged at different positions in the vertical and horizontal directions, respectively, and the circuit is composed of the M receiving antenna elements. Each complex showing propagation characteristics between each of the N transmitting antenna elements and each of the M receiving antenna elements from each of the N received signals received in each of the receiving antenna elements in a predetermined period. By calculating the first matrix of N × M having the transmission function as a component and extracting the second matrix corresponding to the predetermined frequency range in the first matrix, at least one of the breathing, heartbeat, and body movement of the living body. The second matrix corresponding to the component affected by the vital activity including the antenna is extracted, and the second matrix is used to be the three-dimensional position where the living body exists with respect to the sensor and includes the vertical position. A first distance indicating the distance between the living body and the transmitting antenna based on the estimated three-dimensional position, the position of the transmitting antenna, and the position of the receiving antenna, and , The second distance indicating the distance between the living body and the receiving antenna is calculated, the RCS value with respect to the living body is calculated using the first distance and the second distance, and the estimated three-dimensional position and the above-mentioned three-dimensional position and , The operation of the living body is estimated by using the calculated RCS value with time and the information indicating the correspondence relationship stored in the memory.

このため、無線信号を利用することで、生体の動作の推定を、短時間かつ高精度に行うことができる。 Therefore, by using the wireless signal, it is possible to estimate the movement of the living body in a short time and with high accuracy.

また、前記対応関係を示す情報において対応付けられている前記生体の動作は、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換を含み、前記回路は、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体が、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換のいずれの動作を行ったかを推定してもよい。 In addition, the movement of the living body associated with the information indicating the correspondence relationship includes falling, sitting on a chair, sitting on the floor, standing up from the chair, standing up from the floor, jumping, and changing direction. In the circuit, the living body falls over using the estimated three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory. , Sitting on a chair, sitting on a floor, standing up from a chair, standing up from the floor, jumping, and turning.

このため、生体の動作の推定を、より短時間に行うことができる。 Therefore, the movement of the living body can be estimated in a shorter time.

また、前記回路は、前記生体の動作の推定において、前記推定した鉛直位置、または、前記算出したRCS値の時間的変化が所定値よりも大きい期間を、前記生体が動作している動作期間として抽出し、抽出した前記動作期間における時間的変化であって、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定してもよい。 Further, in the estimation of the operation of the living body, the circuit sets a period in which the estimated vertical position or the time change of the calculated RCS value is larger than a predetermined value as an operating period in which the living body is operating. Information indicating the temporal change in the extracted and extracted operating period, the temporal change of the estimated three-dimensional position and the calculated RCS value, and the correspondence relationship stored in the memory. And, may be used to estimate the movement of the living body.

このため、生体の動作の推定にかかる処理負荷を低減することができる。 Therefore, the processing load required for estimating the operation of the living body can be reduced.

また、前記回路は、時系列で得られた複数の前記鉛直位置、または、複数の前記RCS値から、所定のフィルタを用いて瞬時ノイズ成分を除去することで得られる時系列データを用いて、前記動作期間を抽出してもよい。 Further, the circuit uses time-series data obtained by removing an instantaneous noise component from a plurality of the vertical positions or a plurality of the RCS values obtained in a time series using a predetermined filter. The operating period may be extracted.

このため、生体の動作の推定を、より高精度に行うことができる。 Therefore, it is possible to estimate the movement of the living body with higher accuracy.

また、前記対応関係を示す情報において前記生体の動作に対応付けられている、前記鉛直位置および前記RCS値は、事前に、前記生体が前記生体の動作のうちの一の動作を前記所定範囲で行ったときに、前記回路により推定された鉛直位置、および、前記回路により算出されたRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換し、変換した前記方向ベクトルを正規化することで得られた方向コードで示されており、前記回路は、前記生体の動作の推定において、抽出した前記動作期間における時間的変化であって、前記推定した前記3次元位置から得られる前記鉛直位置、および、前記算出したRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換し、変換した前記方向ベクトルを正規化することで方向コードを算出し、算出した前記方向コードと、前記対応関係を示す情報と、を用いて、前記生体の動作を推定してもよい。 Further, the vertical position and the RCS value, which are associated with the movement of the living body in the information indicating the correspondence, are set in advance by the living body in the predetermined range of one of the movements of the living body. When this is done, the vertical position estimated by the circuit and the temporal change of the RCS value calculated by the circuit are converted into a direction vector using a predetermined method, and the converted direction vector is normalized. It is shown by the direction code obtained by the above, and the circuit is a temporal change in the extracted operation period in the estimation of the operation of the living body, and the circuit is obtained from the estimated three-dimensional position. The direction code is calculated by converting the vertical position and the temporal change of the calculated RCS value into a direction vector using a predetermined method, and normalizing the converted direction vector, and the calculated direction code. And the information indicating the correspondence relationship may be used to estimate the movement of the living body.

このため、生体の動作の推定を、より高精度に行うことができる。 Therefore, it is possible to estimate the movement of the living body with higher accuracy.

また、前記回路は、第1の動作期間の次の第2の動作期間において、前記第1の動作期間の終了時の前記生体の姿勢を用いて、前記生体の動作を推定してもよい。 Further, in the second operation period following the first operation period, the circuit may estimate the operation of the living body by using the posture of the living body at the end of the first operation period.

このため、推定した生体の動作を利用して、次の生体の動作を推定できる。よって、生体の動作の推定を、より効果的に行うことができる。 Therefore, the next movement of the living body can be estimated by using the estimated movement of the living body. Therefore, the movement of the living body can be estimated more effectively.

また、前記回路は、さらに、前記推定した前記3次元位置における水平方向への変動が所定の距離以上の変動である場合、前記生体が水平方向に移動していると推定してもよい。 Further, the circuit may further estimate that the living body is moving in the horizontal direction when the horizontal fluctuation at the estimated three-dimensional position is a fluctuation of a predetermined distance or more.

このため、生体の水平方向への移動を、短時間かつ高精度に推定することができる。 Therefore, the horizontal movement of the living body can be estimated in a short time and with high accuracy.

また、前記回路は、さらに、前記生体が水平方向に移動していると推定した場合の前記3次元位置に含まれる前記鉛直位置を用いて、前記生体の身長を推定してもよい。 Further, the circuit may further estimate the height of the living body by using the vertical position included in the three-dimensional position when the living body is estimated to be moving in the horizontal direction.

このため、生体の身長を短時間かつ高精度に推定することができる。これにより、例えば、予め所定範囲に存在しうる複数の生体が分かっており、複数の生体の身長がそれぞれ異なっていれば、複数の生体のうちのどの生体が存在しているか、動作しているかを特定することに利用することができる。 Therefore, the height of the living body can be estimated in a short time and with high accuracy. As a result, for example, if a plurality of organisms that can exist in a predetermined range are known in advance and the heights of the plurality of organisms are different, which of the plurality of organisms exists or is operating. Can be used to identify.

また、前記回路は、さらに、前記生体が水平方向に移動していると推定した場合の前記RCS値を用いて、前記生体の身体のサイズを推定してもよい。 Further, the circuit may further estimate the size of the body of the living body by using the RCS value when it is estimated that the living body is moving in the horizontal direction.

このため、生体の身体のサイズを短時間かつ高精度に推定することができる。これにより、例えば、予め所定範囲に存在しうる複数の生体が分かっており、複数の生体の身体のサイズがそれぞれ異なっていれば、複数の生体のうちのどの生体が存在しているか、動作しているかを特定することに利用することができる。 Therefore, the size of the body of the living body can be estimated in a short time and with high accuracy. As a result, for example, if a plurality of organisms that can exist in a predetermined range are known in advance and the body sizes of the plurality of organisms are different, which organism among the plurality of organisms exists will operate. It can be used to identify if it is.

また、前記所定期間は、前記生体の呼吸、心拍、および、体動の少なくとも1つの周期の略半分であってもよい。 In addition, the predetermined period may be approximately half of at least one cycle of respiration, heartbeat, and body movement of the living body.

このため、生体の動作の推定を効果的に行うことができる。 Therefore, it is possible to effectively estimate the movement of the living body.

なお、本発明は、センサーとして実現するだけでなく、このようなセンサーが備える処理手段を備える集積回路として実現したり、その装置を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを示す情報、データまたは信号として実現したりすることもできる。そして、それらプログラム、情報、データおよび信号は、CD−ROM等の記録媒体やインターネット等の通信媒体を介して配信してもよい。 The present invention is not only realized as a sensor, but also realized as an integrated circuit provided with processing means provided by such a sensor, or as a method in which processing means constituting the device is used as a step. Can be realized as a program that causes a computer to execute, or can be realized as information, data, or a signal indicating the program. Then, the programs, information, data and signals may be distributed via a recording medium such as a CD-ROM or a communication medium such as the Internet.

以下、本発明の実施の形態について、図面を用いて詳細に説明する。なお、以下で説明する実施の形態は、いずれも本発明の好ましい一具体例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本発明の最上位概念を示す独立請求項に記載されていない構成要素については、より好ましい形態を構成する任意の構成要素として説明される。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. It should be noted that all of the embodiments described below show a preferred specific example of the present invention. Numerical values, shapes, materials, components, arrangement positions and connection forms of components, steps, order of steps, etc. shown in the following embodiments are examples, and are not intended to limit the present invention. Further, among the components in the following embodiments, the components not described in the independent claims indicating the highest level concept of the present invention will be described as arbitrary components constituting the more preferable form. In the present specification and the drawings, components having substantially the same functional configuration are designated by the same reference numerals, so that duplicate description will be omitted.

(実施の形態)
図1は実施の形態におけるセンサーの構成の一例を示すブロック図である。図2は、実施の形態におけるセンサーの設置例を示す図である。
(Embodiment)
FIG. 1 is a block diagram showing an example of a sensor configuration according to an embodiment. FIG. 2 is a diagram showing an installation example of the sensor in the embodiment.

図1に示すように、センサー10は、送信アンテナ20、受信アンテナ30、回路40およびメモリ41を備える。センサー10は、ヒト等の生体50に対して送信アンテナ20よりマイクロ波を発射し、受信アンテナ30にて生体50で反射された反射波を受信する。ここで、送信アンテナ20に対して任意に設定された水平面上の方向である第1基準方向と、送信アンテナ20から生体50への方向である第1生体方向とのなす角をφとする。また、鉛直方向と、第1生体方向との為す角である生体50の仰角をθとする。また、受信アンテナ30に対して任意に設定された水平面上の方向である第2基準方向と、受信アンテナ30から生体50への方向である第2生体方向との為す角である生体50の仰角をφとする。また、鉛直方向と、第2生体方向とのなす角をθとする。生体50がバイタル活動を行っている部位の中心座標を(x,y,z)とすると、送信アンテナ20、受信アンテナ30および生体50の位置関係によって、方向(θ,θ,φ,φ)と座標(x,y,z)は相互に変換可能である。 As shown in FIG. 1, the sensor 10 includes a transmitting antenna 20, a receiving antenna 30, a circuit 40, and a memory 41. The sensor 10 emits microwaves from the transmitting antenna 20 to the living body 50 such as a human, and receives the reflected waves reflected by the living body 50 by the receiving antenna 30. Here, let φ T be the angle formed by the first reference direction, which is a direction on the horizontal plane arbitrarily set with respect to the transmitting antenna 20, and the first living body direction, which is the direction from the transmitting antenna 20 to the living body 50. .. Further, the elevation angle of the living body 50, which is the angle formed by the vertical direction and the first living body direction, is set to θ T. Further, the elevation angle of the living body 50, which is the angle formed by the second reference direction, which is a direction on the horizontal plane arbitrarily set with respect to the receiving antenna 30, and the second living body direction, which is the direction from the receiving antenna 30 to the living body 50. Let φ R. Further, the angle formed by the vertical direction and the second living body direction is θ R. Assuming that the central coordinates of the part where the living body 50 is performing vital activity are (x b , y b , z b ), the directions (θ T , θ R , depending on the positional relationship between the transmitting antenna 20, the receiving antenna 30, and the living body 50) φ T , φ R ) and coordinates (x b , y b , z b ) can be converted to each other.

送信アンテナ20は、N個の送信アンテナ素子21を有する。送信アンテナ20は、水平方向(x方向)にN個が並び、かつ、鉛直方向(z方向)にN個が並ぶように矩形配置された、N個(N=N×N)の送信アンテナ素子21で構成されるアレーアンテナを有する。つまり、N個の送信アンテナ素子21のうち少なくとも3個の送信アンテナ素子21は、鉛直方向および水平方向の異なる位置に配置されている。N個の送信アンテナ素子21のそれぞれは、生体が存在しうる所定範囲に対して送信信号を送信する。つまり、送信アンテナ20は、異なるN箇所の位置からN個の送信信号を所定範囲に対して送信する。なお、生体が存在しうる所定範囲とは、センサー10が生体の存在を検知する検知範囲である。 The transmitting antenna 20 has N transmitting antenna elements 21. Transmitting antenna 20, lined N x pieces in the horizontal direction (x-direction), and vertical direction N z number in (z direction) is rectangular arranged side by side, N number (N = N x × N z ) It has an array antenna composed of the transmitting antenna element 21 of the above. That is, at least three of the N transmitting antenna elements 21 are arranged at different positions in the vertical direction and the horizontal direction. Each of the N transmitting antenna elements 21 transmits a transmission signal to a predetermined range in which a living body can exist. That is, the transmitting antenna 20 transmits N transmission signals to a predetermined range from different N positions. The predetermined range in which the living body can exist is a detection range in which the sensor 10 detects the presence of the living body.

N個の送信アンテナ素子21のそれぞれは、具体的には、ヒトなどの生体50に対して、マイクロ波を送信信号として発射する。N個の送信アンテナ素子21は、送信アンテナ素子21毎に異なる変調処理が行われた信号を送信信号として送信してもよい。また、N個の送信アンテナ素子21のそれぞれは、変調信号または無変調の信号を逐次的に切り替えて送信してもよい。変調処理は、送信アンテナ20により行われてもよい。このように、N個の送信アンテナ素子21毎に、N個の送信アンテナ素子21から送信される送信信号をそれぞれ異なる送信信号とすることで、受信アンテナ30により受信された送信信号を送信した送信アンテナ素子21を特定できる。このように、送信アンテナ20は、変調処理を行うための回路を含んでいてもよい。 Specifically, each of the N transmitting antenna elements 21 emits microwaves as transmission signals to a living body 50 such as a human being. The N transmitting antenna elements 21 may transmit a signal subjected to different modulation processing for each transmitting antenna element 21 as a transmitting signal. Further, each of the N transmitting antenna elements 21 may sequentially switch and transmit a modulated signal or an unmodulated signal. The modulation process may be performed by the transmitting antenna 20. In this way, by setting the transmission signals transmitted from the N transmission antenna elements 21 to different transmission signals for each of the N transmission antenna elements 21, the transmission signal received by the reception antenna 30 is transmitted. The antenna element 21 can be specified. As described above, the transmitting antenna 20 may include a circuit for performing the modulation process.

受信アンテナ30は、M個の受信アンテナ素子31を有する。受信アンテナ30は、水平方向(x方向)にM個が並び、かつ、鉛直方向(z方向)にM個が並ぶように矩形配置された、M個(M=M×M)の受信アンテナ素子31で構成されるアレーアンテナを有する。つまり、M個の受信アンテナ素子31のうち少なくとも3個の受信アンテナ素子31は、鉛直方向および水平方向の異なる位置に配置されている。M個の受信アンテナ素子31のそれぞれは、N個の送信信号のうち生体50により反射された信号である反射信号を含むN個の受信信号を受信する。受信アンテナ30は、マイクロ波からなる受信信号を周波数変換し、低周波数信号に変換する。受信アンテナ30は、低周波数信号に変換することにより得られた信号を回路40に出力する。つまり、受信アンテナ30は、受信信号を処理するための回路を含んでいてもよい。 The receiving antenna 30 has M receiving antenna elements 31. Receiving antenna 30, lined M x pieces in the horizontal direction (x-direction), and vertical M z number in (z direction) is rectangular arranged side by side, M pieces (M = M x × M z ) It has an array antenna composed of the receiving antenna element 31 of the above. That is, at least three of the M receiving antenna elements 31 are arranged at different positions in the vertical direction and the horizontal direction. Each of the M receiving antenna elements 31 receives N reception signals including a reflection signal which is a signal reflected by the living body 50 among the N transmission signals. The receiving antenna 30 frequency-converts a received signal composed of microwaves and converts it into a low-frequency signal. The receiving antenna 30 outputs the signal obtained by converting it into a low frequency signal to the circuit 40. That is, the receiving antenna 30 may include a circuit for processing the received signal.

回路40は、センサー10を動作させる各種処理を実行する。回路40は、例えば、制御プログラムを実行するプロセッサと、当該制御プログラムを実行するときに使用するワークエリアとして用いられる揮発性の記憶領域(主記憶装置)とにより構成される。揮発
性の記憶領域は、例えば、RAM(Randdom Access Memory)である。なお、回路40は、センサー10を動作させる各種処理を行うための専用回路により構成されていてもよい。つまり、回路40は、ソフトウェア処理を行う回路であってもよいし、ハードウェア処理を行う回路であってもよい。
The circuit 40 executes various processes for operating the sensor 10. The circuit 40 is composed of, 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, RAM (Random Access Memory). The circuit 40 may be configured by a dedicated circuit for performing various processes for operating the sensor 10. That is, the circuit 40 may be a circuit that performs software processing or may be a circuit that performs hardware processing.

メモリ41は、不揮発性の記憶領域(補助記憶装置)であり、例えば、ROM(Read Only Memory)、フラッシュメモリ、HDD(Hard Disk Drive)などである。メモリ41は、例えば、センサー10を動作させる各種処理に利用される情報を記憶している。 The memory 41 is a non-volatile storage area (auxiliary storage device), and is, for example, a ROM (Read Only Memory), a flash memory, an HDD (Hard Disk Drive), or the like. The memory 41 stores, for example, information used for various processes for operating the sensor 10.

次に、回路40の機能的な構成について図3を用いて説明する。 Next, the functional configuration of the circuit 40 will be described with reference to FIG.

図3は実施の形態における回路およびメモリの機能的な構成を示すブロック図である。 FIG. 3 is a block diagram showing a functional configuration of a circuit and a memory according to an embodiment.

回路40は、複素伝達関数算出部410と、生体成分算出部420と、位置推定処理部430と、RCS算出部440と、動作推定部450とを有する。 The circuit 40 includes a complex transfer function calculation unit 410, a biological component calculation unit 420, a position estimation processing unit 430, an RCS calculation unit 440, and an operation estimation unit 450.

複素伝達関数算出部410は、低周波信号に変換された受信信号から複素伝達関数を算出する。複素伝達関数とは、各送信アンテナ素子21と各受信アンテナ素子31との間の伝搬損失および位相回転を表すものである。複素伝達関数は、送信アンテナ素子数がN個であり、受信アンテナ素子数がM個の場合、M×Nの成分を持つ複素行列となる。以降、この複素行列を複素伝達関数行列と呼ぶ。推定した複素伝達関数行列は、生体成分算出部420に出力される。つまり、複素伝達関数算出部410は、M個の受信アンテナ素子31のそれぞれにおいて所定期間で受信された複数の受信信号のそれぞれから、N個の送信アンテナ素子21それぞれと、M個の受信アンテナ素子31それぞれとの間の伝播特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出する。 The complex transfer function calculation unit 410 calculates the complex transfer function from the received signal converted into the low frequency signal. The complex transfer function represents the propagation loss and phase rotation between each transmitting antenna element 21 and each receiving antenna element 31. The complex transfer function is a complex matrix having an M × N component when the number of transmitting antenna elements is N and the number of receiving antenna elements is M. Hereinafter, this complex matrix will be referred to as a complex transfer function matrix. The estimated complex transfer function matrix is output to the biological component calculation unit 420. That is, the complex transfer function calculation unit 410 has N transmission antenna elements 21 and M reception antenna elements from each of the plurality of reception signals received in each of the M reception antenna elements 31 in a predetermined period. 31 A first matrix of N × M is calculated, which is composed of each complex transfer function showing the propagation characteristics with each of them.

生体成分算出部420は、生体50を経由した受信信号から得られた複素伝達関数行列成分と、生体50を経由していない受信信号から得られた複素伝達関数行列成分とに分離する。生体50を経由した成分とは、生体活動により時変動する成分である。よって、生体50を経由した成分は、例えば、生体50以外は静止しているものとした場合、複素伝達関数行列の成分を時間方向にフーリエ変換することで得られた成分から、直流以外の成分を取り出すことによって抽出することが可能である。また、生体50を経由した成分は、例えば、生体50が所定範囲に存在しないときに観測された結果との差分が所定の閾値を超えている成分を取り出すことによって抽出することも可能である。このように、生体成分算出部420は、生体50を経由した反射信号を含む受信信号から得られた複素伝達関数行列成分を抽出することで、抽出した複素伝達関数行列成分を生体成分として算出する。つまり、生体成分算出部420は、第1行列における所定周波数範囲に対応する第2行列を抽出することで、生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する第2行列を抽出する。所定周波数範囲は、例えば、上述した生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動に由来する周波数である。所定周波数範囲は、例えば、0.1Hz以上3Hz以下の範囲の周波数である。これにより、心臓、肺、横隔膜、内蔵の動きによる生体50の部位のバイタル活動、または、手、足などによるバイタル活動の影響を受けた生体成分を抽出できる。なお、心臓、肺、横隔膜、内蔵の動きによる生体50の部位とは、例えば、人のみぞおちである。 The biological component calculation unit 420 separates the complex transfer function matrix component obtained from the received signal passing through the living body 50 and the complex transfer function matrix component obtained from the received signal not passing through the living body 50. The component that has passed through the living body 50 is a component that fluctuates with time due to biological activity. Therefore, the component that has passed through the living body 50 is, for example, a component other than the direct current from the component obtained by Fourier transforming the component of the complex transfer function matrix in the time direction, assuming that the components other than the living body 50 are stationary. Can be extracted by taking out. Further, the component that has passed through the living body 50 can be extracted, for example, by extracting a component whose difference from the result observed when the living body 50 does not exist in a predetermined range exceeds a predetermined threshold value. In this way, the biological component calculation unit 420 calculates the extracted complex transfer function matrix component as a biological component by extracting the complex transfer function matrix component obtained from the received signal including the reflected signal passing through the living body 50. .. That is, the biological component calculation unit 420 extracts the second matrix corresponding to the predetermined frequency range in the first matrix, and thereby, the component affected by the vital activity including at least one of the respiration, heartbeat, and body movement of the living body. The second matrix corresponding to is extracted. The predetermined frequency range is, for example, a frequency derived from vital activity including at least one of the above-mentioned respiration, heartbeat, and body movement of a living body. The predetermined frequency range is, for example, a frequency in the range of 0.1 Hz or more and 3 Hz or less. This makes it possible to extract biological components affected by vital activity of 50 parts of the living body due to the movement of the heart, lungs, diaphragm, and internal organs, or vital activity of hands, feet, and the like. The parts of the living body 50 due to the movement of the heart, lungs, diaphragm, and internal organs are, for example, the epigastrium of a human being.

ここで生体成分は、M×Nの成分を持つ行列であり、所定期間に受信アンテナ30において観測された受信信号から得られる複素伝達関数から抽出される。このため、生体成分
は、周波数応答あるいは時間応答情報を持っているものとする。なお、所定期間は、生体の呼吸、心拍、および、体動の少なくとも1つの周期の略半分の期間である。
Here, the biological component is a matrix having an M × N component, and is extracted from the complex transfer function obtained from the received signal observed by the receiving antenna 30 in a predetermined period. Therefore, it is assumed that the biological component has frequency response or time response information. The predetermined period is approximately half the period of at least one cycle of respiration, heartbeat, and body movement of the living body.

生体成分算出部420で算出された生体成分は、位置推定処理部430に出力される。位置推定処理部430は、算出された生体成分を用いて生体の位置推定を行う。つまり、位置推定処理部430は、第2行列を用いて、センサー10に対する生体50の存在する3次元位置であって、センサー10に対する生体の存在する鉛直位置を含む3次元位置を推定する。位置推定には、送信アンテナ20からの出発角θと受信アンテナ30への到来角θとの両方の角度を推定し、推定した出発角θおよび到来角θから三角法によって生体50の位置を推定する。 The biological component calculated by the biological component calculation unit 420 is output to the position estimation processing unit 430. The position estimation processing unit 430 estimates the position of the living body using the calculated biological components. That is, the position estimation processing unit 430 estimates the three-dimensional position where the living body 50 exists with respect to the sensor 10 and includes the vertical position where the living body exists with respect to the sensor 10 by using the second matrix. The position estimation, both angle and departure angle theta T and the arrival angle theta R to the receiving antenna 30 from the transmitting antenna 20 is estimated, biological 50 by trigonometry from the departure angle theta T were estimated and the arrival angle theta R Estimate the position of.

なお、位置推定には、観測対象である生体50と送信アンテナ20との第1距離、または、生体50と受信アンテナ30との第2距離が、送信アンテナ20または受信アンテナ30を構成するアレーアンテナの開口と比べて同程度となるほど近接する場合は、球面波モードベクトルと呼ばれるものを用いて生体50の位置推定を行ってもよい。これは、送信アンテナ20または受信アンテナ30を構成するアレーアンテナ素子毎に出発角または到来角が異なるためである。この場合、平面波のように出発角または到来角を定義できないため、送信アンテナ20または受信アンテナ30を構成するアレーアンテナ素子と対象の生体50との位置関係を表す2次元または3次元の座標を用いて後述のステアリングベクトルを定義する。 For position estimation, the first distance between the living body 50 and the transmitting antenna 20 to be observed, or the second distance between the living body 50 and the receiving antenna 30, constitutes the transmitting antenna 20 or the receiving antenna 30. When the antennae are close to each other to the same extent as the opening of the living body 50, the position of the living body 50 may be estimated using what is called a spherical wave mode vector. This is because the departure angle or the arrival angle is different for each of the array antenna elements constituting the transmitting antenna 20 or the receiving antenna 30. In this case, since the departure angle or the arrival angle cannot be defined like a plane wave, two-dimensional or three-dimensional coordinates representing the positional relationship between the array antenna element constituting the transmitting antenna 20 or the receiving antenna 30 and the target living body 50 are used. The steering vector described later is defined.

なお、位置推定処理部430は、用途に応じて、例えば、3秒以上の第1測定期間において測定することにより得られた測定データにより静止中の生体50の位置を推定すると共に、例えば、1秒、0.5秒などの第1測定期間よりも短い第2測定期間において測定することにより得られた測定データを用いて、動作中の生体50の位置を推定しても良い。この時の概念を図4Aおよび図4Bに示す。 The position estimation processing unit 430 estimates the position of the stationary living body 50 from the measurement data obtained by measuring in the first measurement period of, for example, 3 seconds or more, depending on the application, and for example, 1 The position of the living body 50 in operation may be estimated using the measurement data obtained by measuring in the second measurement period shorter than the first measurement period such as seconds and 0.5 seconds. The concept at this time is shown in FIGS. 4A and 4B.

図4Aは、異なる測定期間で測定するセンサーの位置分解能の一例について説明するための図である。 FIG. 4A is a diagram for explaining an example of the position resolution of the sensor measuring in different measurement periods.

図4Aに示すように、例えば、格子301による複数の四角形の領域は、例えば、格子301の全体の範囲がセンサー10の測定範囲である所定範囲を示しており、かつ、8m×8mの範囲の空間である場合、100cm四方の分解能で第1測定期間において測定される場合の複数の領域を示す。また、格子302による複数の四角形の領域は、格子302の全体の範囲がセンサー10の測定範囲である所定範囲を示しており、かつ、格子301と同じ8m×8mの範囲の空間である場合、5cm四方の分解能で第2測定期間において測定される場合の複数の領域を示す。上述したように、第1測定期間は、例えば、3秒であり、第2測定期間は、例えば、0.5秒である。また、格子301で区切られた複数の領域のそれぞれを第1測定期間で測定するのに要する第1期間と、格子302で区切られた複数の領域のそれぞれを第2測定期間で測定するのに要する第2期間とは、互いに重複する期間である。つまり、格子301で区切られた複数の領域のそれぞれを測定する処理と、格子302で区切られた複数の領域のそれぞれを測定する処理とは、並列処理により、実行される。 As shown in FIG. 4A, for example, the plurality of quadrangular regions formed by the grid 301 indicate a predetermined range in which the entire range of the grid 301 is the measurement range of the sensor 10, and the range is 8 m × 8 m. In the case of space, a plurality of regions when measured in the first measurement period with a resolution of 100 cm square are shown. Further, when the plurality of quadrangular regions formed by the grid 302 indicate a predetermined range in which the entire range of the grid 302 is the measurement range of the sensor 10 and the space is the same as the grid 301 in the range of 8 m × 8 m. Multiple regions when measured in the second measurement period with a resolution of 5 cm square are shown. As described above, the first measurement period is, for example, 3 seconds, and the second measurement period is, for example, 0.5 seconds. Further, in order to measure each of the plurality of regions separated by the grid 301 in the first measurement period and each of the plurality of regions separated by the grid 302 in the second measurement period. The required second period is a period that overlaps with each other. That is, the process of measuring each of the plurality of regions separated by the grid 301 and the process of measuring each of the plurality of regions separated by the grid 302 are executed by parallel processing.

異なる測定期間で生体50を測定する処理は、下記のように行ってもよい。 The process of measuring the living body 50 in different measurement periods may be performed as follows.

図4Bは、異なる測定期間で測定するセンサーの位置分解能の他の一例について説明するための図である。 FIG. 4B is a diagram for explaining another example of the position resolution of the sensor measuring in different measurement periods.

図4Bに示すように、例えば、格子301による複数の四角形の領域は、図4Aの説明と同様に、格子301の全体の範囲がセンサー10の測定範囲である所定範囲を示しており、かつ、8m×8mの範囲の空間である場合、100cm四方の分解能で第1測定期間において測定される場合の複数の領域を示す。また、格子303による複数の四角形の領域は、格子303の全体の範囲が格子301による複数の領域の測定により検出された生体50の位置を含む領域である4m×4mの範囲の空間において、30cm四方の分解能で第2測定期間において測定される場合の複数の領域を示す。また、格子304による複数の四角形の領域は、格子304の全体の範囲が格子303による複数の領域の測定により検出された生体50の位置を含む領域である2m×2mの範囲の空間において、10cm四方の分解能で第2測定期間よりも短い第3測定期間において測定される場合の複数の領域を示す。第1測定期間は、例えば、3秒であり、第2測定期間は、例えば、1秒であり、第3測定期間は、例えば、0.5秒である。なお、静止中の生体50の位置の検出精度をより向上させる場合、第1測定期間を10〜20秒など3秒よりも長い時間に設定してもよい。 As shown in FIG. 4B, for example, the plurality of quadrangular regions formed by the grid 301 indicate a predetermined range in which the entire range of the grid 301 is the measurement range of the sensor 10, as in the description of FIG. 4A. In the case of a space in the range of 8 m × 8 m, a plurality of regions when measured in the first measurement period with a resolution of 100 cm square are shown. Further, the plurality of quadrangular regions formed by the grid 303 are 30 cm in a space of 4 m × 4 m, which is a region in which the entire range of the grid 303 includes the position of the living body 50 detected by the measurement of the plurality of regions by the grid 301. A plurality of regions when measured in the second measurement period with four-sided resolution are shown. Further, the plurality of quadrangular regions formed by the lattice 304 are 10 cm in a space having a range of 2 m × 2 m, in which the entire range of the lattice 304 includes the position of the living body 50 detected by the measurement of the plurality of areas by the lattice 303. A plurality of regions when measured in a third measurement period shorter than the second measurement period with four-sided resolution are shown. The first measurement period is, for example, 3 seconds, the second measurement period is, for example, 1 second, and the third measurement period is, for example, 0.5 seconds. In order to further improve the detection accuracy of the position of the living body 50 at rest, the first measurement period may be set to a time longer than 3 seconds, such as 10 to 20 seconds.

RCS算出部440は、生体成分と推定された位置とを用いて散乱断面積(RCS:Radar Cross Section)を算出する。RCS算出部440は、具体的には、散乱断面積を計算するために、推定された3次元位置と、送信アンテナ20の位置と、受信アンテナ30の位置と、に基づいて、生体50と送信アンテナ20との第1距離を示す距離RT、および、生体50と受信アンテナ30との第2距離を示す距離RRを算出する。RCS算出部440は、算出した距離RTおよび距離RRから伝搬距離を算出し、算出した伝播距離と生体成分の強度とを用いてRCSを算出する。なお、送信アンテナ20の位置と、受信アンテナ30の位置とは、メモリ41に予め記憶されていてもよい。 The RCS calculation unit 440 calculates the scattering cross section (RCS: Radar Cross Section) using the biological component and the estimated position. Specifically, the RCS calculation unit 440 transmits the living body 50 and the living body 50 based on the estimated three-dimensional position, the position of the transmitting antenna 20, and the position of the receiving antenna 30 in order to calculate the scattering cross section. The distance RT indicating the first distance to the antenna 20 and the distance RR indicating the second distance between the living body 50 and the receiving antenna 30 are calculated. The RCS calculation unit 440 calculates the propagation distance from the calculated distance RT and the distance RR, and calculates the RCS using the calculated propagation distance and the intensity of the biological component. The position of the transmitting antenna 20 and the position of the receiving antenna 30 may be stored in the memory 41 in advance.

動作推定部450は、位置推定処理部430により推定された3次元位置、および、RCS算出部440により算出されたRCS値の時間的変化を示す時系列データと、予めメモリ41に記憶されている対応関係を示す情報42と、を用いて、生体50の動作を推定する。動作推定部450は、動作期間抽出部451と、方向コード算出部452と、動作比較部453とを有する。 The motion estimation unit 450 stores in advance the three-dimensional position estimated by the position estimation processing unit 430 and the time-series data indicating the temporal change of the RCS value calculated by the RCS calculation unit 440 in the memory 41. The movement of the living body 50 is estimated using the information 42 indicating the correspondence relationship. The motion estimation unit 450 includes an operation period extraction unit 451, a direction code calculation unit 452, and an operation comparison unit 453.

動作期間抽出部451は、図5に示すように位置推定処理部430により推定された生体50の3次元位置、または、RCS算出部440により算出されたRCS値の時間的変化の変化幅が所定値よりも大きい期間を動作期間として抽出する。なお、図5は、高さ位置(Height)またはRCE値の時系列データから動作期間を抽出する例について説明するための図である。 As shown in FIG. 5, the operation period extraction unit 451 determines the three-dimensional position of the living body 50 estimated by the position estimation processing unit 430 or the change width of the time change of the RCS value calculated by the RCS calculation unit 440. The period larger than the value is extracted as the operation period. Note that FIG. 5 is a diagram for explaining an example of extracting the operation period from the time series data of the height position (Height) or the RCE value.

動作期間抽出部451は、例えば、鉛直位置である高さ位置またはRCS値を用いて動作期間を抽出する場合、瞬時ノイズの影響を避けるため、得られた3次元位置またはRCS値の時系列のデータに対し、例えば、メディアンフィルタ、FIRフィルタ、平均値などを用い、高さ位置、および、RCS値のノイズ成分を除去し、フィルタ処理後の高さ情報の変化区間やRCS変化区間を当該生体の動作期間として抽出してもよい。つまり、動作期間抽出部451は、時系列で得られた複数の鉛直位置、または、複数のRCS値から、所定のフィルタを用いて瞬時ノイズ成分を除去することで得られる時系列データを用いて、動作期間を抽出してもよい。図4では、一例として測定期間が0.6秒程度の測定データを用いて高さ位置およびRCS値を測定し、所定のフィルタ処理をした後の様子、および動作期間として抽出した様子を示す。なお、動作期間抽出部451は、計算量削減目的などで、推定する対象となる期間を限定したい時に有効であるが必ずしも設けられる必要はない。つまり、全区間について状態推定を行う場合、動作期間抽出部45を省略し、全区間について推定しても良いことは言うまでもない。 When the operation period extraction unit 451 extracts the operation period using, for example, a height position that is a vertical position or an RCS value, in order to avoid the influence of instantaneous noise, the operation period extraction unit 451 is a time series of the obtained three-dimensional position or RCS value. For the data, for example, a median filter, an FIR filter, an average value, etc. are used to remove the noise component of the height position and the RCS value, and the change section of the height information and the RCS change section after the filtering process are set to the living body. It may be extracted as the operation period of. That is, the operation period extraction unit 451 uses time-series data obtained by removing an instantaneous noise component from a plurality of vertical positions or a plurality of RCS values obtained in a time series using a predetermined filter. , The operating period may be extracted. FIG. 4 shows, as an example, a state in which the height position and the RCS value are measured using measurement data having a measurement period of about 0.6 seconds, a state after performing a predetermined filter process, and a state of extraction as an operation period. The operation period extraction unit 451 is effective when it is desired to limit the period to be estimated for the purpose of reducing the amount of calculation, but it is not always necessary to provide the operation period extraction unit 451. That is, when performing state estimation for all sections, it goes without saying that the operation period extraction unit 45 may be omitted and estimation may be performed for all sections.

方向コード算出部452は、動作期間抽出部451により抽出された動作期間における時間的変化であって、推定された3次元位置から得られる鉛直位置(高さ位置)、および、算出されたRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換する。具体的には、方向コード算出部452は、図6示すように高さ位置とRCS値とを2次元プロットし、その時間変化における軌跡において、当該軌跡の距離ΔPと、当該軌跡の方向θとを算出する。方向コード算出部452は、例えば、第1のタイミングにおける高さ位置H1およびRCS値R1で示される第1の座標p1(H1、R1)から、第1のタイミングの次の第2のタイミングにおける高さ位置H2およびRCS値R2で示される第2の座標p2(H2、R2)への軌跡において、第1の座標p1(H1、R1)と第2の座標p2(H2、R2)との間の距離ΔPと、第1の座標p1(H1、R1)から第2の座標p2(H2、R2)を見たときの方向θとを算出することで方向ベクトルに変換する。図6は、方向ベクトルに変換する処理について説明するための図である。 The direction code calculation unit 452 is a temporal change in the operation period extracted by the operation period extraction unit 451, and is a vertical position (height position) obtained from the estimated three-dimensional position and a calculated RCS value. The temporal change of is converted into a direction vector using a predetermined method. Specifically, the direction code calculation unit 452 plots the height position and the RCS value two-dimensionally as shown in FIG. 6, and in the locus in the time change, the distance ΔP of the locus and the direction θ of the locus Is calculated. The direction code calculation unit 452 is, for example, from the height position H1 at the first timing and the height at the second timing following the first timing from the first coordinates p1 (H1, R1) indicated by the RCS value R1. In the locus to the second coordinate p2 (H2, R2) indicated by the position H2 and the RCS value R2, between the first coordinate p1 (H1, R1) and the second coordinate p2 (H2, R2). It is converted into a direction vector by calculating the distance ΔP and the direction θ when the second coordinates p2 (H2, R2) are viewed from the first coordinates p1 (H1, R1). FIG. 6 is a diagram for explaining a process of converting into a direction vector.

次に、方向コード算出部452は、変換した方向ベクトルを正規化することで方向コードを算出する。具体的には、方向コード算出部452は、図7に示す方向コード表を参照することで、方向コードを算出する。例えば、方向コード算出部452は、1〜8で示される方向コードのうち方向θが最も近い方向コードを特定する。図7は、方向コード表の一例を示す図である。 Next, the direction code calculation unit 452 calculates the direction code by normalizing the converted direction vector. Specifically, the direction code calculation unit 452 calculates the direction code by referring to the direction code table shown in FIG. 7. For example, the direction code calculation unit 452 specifies the direction code closest to the direction θ among the direction codes shown in 1 to 8. FIG. 7 is a diagram showing an example of a direction code table.

方向コード算出部452は、上述のように方向コードと、距離ΔPとを算出することで図8に示すような方向コードの時系列データを得る。図8は、算出された方向コードおよび距離の時系列データの一例を示す図である。なお、この時、方向コード算出部452は、個人差による影響を避けるため、方向コードの正規化を行っても良い。 The direction code calculation unit 452 obtains time-series data of the direction code as shown in FIG. 8 by calculating the direction code and the distance ΔP as described above. FIG. 8 is a diagram showing an example of time-series data of the calculated direction code and distance. At this time, the direction code calculation unit 452 may normalize the direction code in order to avoid the influence of individual differences.

動作比較部453は、方向コード算出部452で算出された方向コードの時系列データと、メモリに記憶されている対応関係を示す情報42とを比較することにより、対応関係を示す情報42において当該時系列データに対応付けられている動作を特定することで、生体50の動作を推定する。 The operation comparison unit 453 compares the time-series data of the direction code calculated by the direction code calculation unit 452 with the information 42 indicating the correspondence relationship stored in the memory, thereby indicating the correspondence relationship in the information 42. The motion of the living body 50 is estimated by specifying the motion associated with the time series data.

なお、メモリ41に記憶されている対応関係を示す情報42とは、センサー10に対する生体50の存在する鉛直方向における位置である鉛直位置およびRCS値の時間的変化を示す複数のモデルコードと、生体50の動作との対応関係を示す情報である。また、対応関係を示す情報42において対応付けられている生体50の動作は、図9に示すように、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換を含む。つまり、動作推定部450は、位置推定処理部430により推定された3次元位置、および、RCS算出部440により算出されたRCS値の時間的変化と、予めメモリ41に記憶されている対応関係を示す情報42と、を用いて、生体50が、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換のいずれの動作を行ったかを推定する。なお、モデルコードは、図10に示されるような時系列データとして表される。 The information 42 indicating the correspondence relationship stored in the memory 41 includes a plurality of model codes indicating the vertical position of the living body 50 with respect to the sensor 10 in the vertical direction and the temporal change of the RCS value, and the living body. It is information which shows the correspondence relation with the operation of 50. Further, as shown in FIG. 9, the movements of the living body 50 associated with the information 42 indicating the correspondence relationship include falling, sitting on a chair, sitting on the floor, standing up from the chair, standing up from the floor, and jumping. , And change direction. That is, the motion estimation unit 450 determines the three-dimensional position estimated by the position estimation processing unit 430, the temporal change of the RCS value calculated by the RCS calculation unit 440, and the correspondence relationship stored in the memory 41 in advance. Using the information 42 shown, it is estimated whether the living body 50 has performed a fall, sitting on a chair, sitting on the floor, standing up from the chair, standing up from the floor, jumping, or changing direction. To do. The model code is represented as time-series data as shown in FIG.

なお、回路40では、複数の異なるタイミングで、繰り返し上記で説明した各部410〜450における処理を行うことにより、時系列データを得る。例えば、回路40は、図4Aまたは図4Bを用いて説明したような、所定のサンプリング周期で処理を繰り返し行うことで時系列の複数の3次元位置および時系列の複数のRCS値からなる時系列データを得る。 In the circuit 40, time-series data is obtained by repeatedly performing the processes in the respective units 41 to 450 described above at a plurality of different timings. For example, the circuit 40 is a time series consisting of a plurality of three-dimensional positions in the time series and a plurality of RCS values in the time series by repeating the process at a predetermined sampling cycle as described with reference to FIG. 4A or FIG. 4B. Get the data.

次に実施の形態のセンサー10の動作原理の詳細を、数式を用いて説明する。なおここでは、フーリエ変換を用いて生体成分を抽出する方法について示す。ここで説明する処理は、回路40により行われる。L人が存在する屋内環境で、送信アンテナ20にM素子、受信アンテナ30にM素子の平面アレーアンテナを使用したとき測定されるM×M時変動MIMOチャネルH(t)は、次の式1と表される。 Next, the details of the operating principle of the sensor 10 of the embodiment will be described using a mathematical formula. Here, a method of extracting a biological component using a Fourier transform will be described. The process described here is performed by the circuit 40. In indoor environments where L person is present, M T elements to the transmission antenna 20 is measured when using a planar array antenna of M R element to the receiving antenna 30 M R × M T when variation MIMO channel H (t) is It is expressed as the following equation 1.

Figure 0006868846
Figure 0006868846

ここで、tは観測時間を表し、第(i,j)要素であるhijはj番目の送信アンテナ素子21からi番目の受信アンテナ素子31へのチャネル応答を表す。 Here, t represents the observation time, and h ij, which is the (i, j) th element, represents the channel response from the j-th transmitting antenna element 21 to the i-th receiving antenna element 31.

×MMIMOアレーは、M×1SIMO(Single−Input Multiple−Output)構成で表現されるMIMO仮想アレーへの変換が可能である。このとき,M×MMIMOチャネルh(t)は、次の式2で表されるM×1仮想SIMOチャネルへと変換される。 M R × M T MIMO array is capable of converting into M R M T × 1SIMO (Single -Input Multiple-Output) MIMO virtual array represented in the configuration. In this case, M R × M T MIMO channel h (t) is converted into M R M T × 1 virtual SIMO channel represented by the following formula 2.

Figure 0006868846
Figure 0006868846

ここで、{・}は転置を表す。ここで差分時間Tを用い、差分チャネルhsb(t,T)を次の式3と定義する。 Here, {・} T represents transpose. Here, the difference time T is used, and the difference channel h sb (t, T) is defined as the following equation 3.

Figure 0006868846
Figure 0006868846

実際の複素チャネルには、直接波、固定物由来の反射波など、生体を経由しない反射波が含まれているが、差分チャネル行列には、差分演算によって生体を経由しない全ての反射波が消去される。このため、複素チャネルには、生体由来の反射波のみが含まれる。 The actual complex channel includes reflected waves that do not pass through the living body, such as direct waves and reflected waves derived from fixed objects, but the difference channel matrix eliminates all reflected waves that do not pass through the living body by the difference calculation. Will be done. Therefore, the complex channel contains only the reflected wave derived from the living body.

ここで差分チャネルhsb(t,T)を用い、ある観測時間tにおける差分時間Tの瞬時相関行列R(t,T)を次の式4と定義する。 Here, using the difference channel h sb (t, T), the instantaneous correlation matrix R (t, T) of the difference time T at a certain observation time t is defined as the following equation 4.

Figure 0006868846
Figure 0006868846

ここで、{・}は複素共役転置を表す。この瞬時相関行列のランクは1であるが、平均演算によって相関行列のランクを回復することが可能であり、これによって複数の到来波の同時推定が可能になる。 Here, {・} H represents a complex conjugate transpose. Although the rank of this instantaneous correlation matrix is 1, it is possible to recover the rank of the correlation matrix by averaging, which enables simultaneous estimation of a plurality of incoming waves.

次に、差分チャネル行列から求めた相関行列を用いて生体の3次元方向推定を行う方法について説明する。ここではMUSICアルゴリズムに基づく推定法について説明する。先の相関行列Rを固有値分解すると、次の式5〜7と書ける。 Next, a method of estimating the three-dimensional direction of the living body using the correlation matrix obtained from the difference channel matrix will be described. Here, an estimation method based on the MUSIC algorithm will be described. When the above correlation matrix R is decomposed into eigenvalues, it can be written as the following equations 5 to 7.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

ここでUは固有ベクトル、Λは固有ベクトルに対応する固有値であり、次の式8で示す順であるものとする。 Here, U is an eigenvector and Λ is an eigenvalue corresponding to the eigenvector, and it is assumed that the order is shown by the following equation 8.

Figure 0006868846
Figure 0006868846

Lは、到来波の数つまり検出対象の生体の数である。以降の説明では、送信アンテナ20または受信アンテナ30を構成するアレーアンテナの開口と比べて検出対象である生体との距離が比較的近い状態を想定して説明を行う。よって当該アレーアンテナでは球面波が観測されるものと仮定する。なお、検出対象である生体との距離が十分に遠方である場合であっても下記の数式は成立するため、検出に支障はない。十分に遠方であることが分かっている場合は、出発角θ,到来角θを使って対象の位置を推定しても良く、計算が比較的簡易になるという利点が得られる。送信アンテナ20側のアレーアンテナのステアリングベクトルは、次式9〜11と定義される。 L is the number of incoming waves, that is, the number of living organisms to be detected. In the following description, it is assumed that the distance to the living body to be detected is relatively short compared to the opening of the array antenna constituting the transmitting antenna 20 or the receiving antenna 30. Therefore, it is assumed that a spherical wave is observed in the array antenna. Even when the distance from the living body to be detected is sufficiently long, the following mathematical formula holds, so that there is no problem in detection. If it is known that the distance is sufficiently large, the position of the target may be estimated using the starting angle θ T and the arrival angle θ R, which has the advantage of making the calculation relatively simple. The steering vector of the array antenna on the transmitting antenna 20 side is defined by the following equations 9 to 11.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

受信アンテナ30側のアレーアンテナのステアリングベクトルも同様に、次式12〜14と定義される。 Similarly, the steering vector of the array antenna on the receiving antenna 30 side is also defined by the following equations 12 to 14.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

ここで、dmxmy,dnxnyは、波源、つまり生体とm,n番のアレー素子との距離を示し、dmxcmyc,dnxcnycは、波源と基準素子との距離を示し、Θnxny,Φnxnyは位相遅れを示し、λは波長を示す。 Here, d mxmy, d nxny is the wave source, i.e. represents the distance between the living body and the m x m y, n x n y th array element, d mxcmyc, d nxcnyc represents the distance between the wave source and the reference element , Θ nxny , Φ nxny indicate the phase lag, and λ indicates the wavelength.

さらに、送受信のステアリングベクトルを乗算し、次の式15に示すように、送受信双方の角度情報を考慮したステアリングベクトルを定義する。 Further, the transmission / reception steering vector is multiplied, and as shown in the following equation 15, the steering vector is defined in consideration of the angle information of both transmission / reception.

Figure 0006868846
Figure 0006868846

これにMUSIC法を適用すると、このステアリングベクトルを用いて、次の式16で示す評価関数の極大値を探索することで到来波の方向を推定する。 When the MUSIC method is applied to this, the direction of the incoming wave is estimated by searching for the maximum value of the evaluation function represented by the following equation 16 using this steering vector.

Figure 0006868846
Figure 0006868846

ここで,この探索は空間の座標(x,y,z)について行うため,3次元の探索処理を実施する。 Here, since this search is performed on the coordinates of space (x, y, z), a three-dimensional search process is performed.

さらに探索によって得た位置情報(x,y)を用いて生体からの散乱断面積(RCS)を計算する。 Furthermore, the scattering cross section (RCS) from the living body is calculated using the position information (x, y) obtained by the search.

観測した伝搬チャネル行列H(t)をフーリエ変換しベクトル化した周波数応答行列F(ω)を用いると、受信電力Pγ(ω)は、次式17と表される。 Using the frequency response matrix F (ω) obtained by Fourier transforming the observed propagation channel matrix H (t) and vectorizing it, the received power P γ (ω) is expressed by the following equation 17.

Figure 0006868846
Figure 0006868846

このとき、生体の散乱断面積σは、次式18〜20と表すことができる。 At this time, the scattering cross section σ of the living body can be expressed by the following equations 18 to 20.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

このとき,Rは受信アンテナ30から推定された生体50の位置までの距離を示し、Rは送信アンテナ20から推定された生体50の位置までの距離を示す。また、Gは受信アンテナ30の利得を示し、Gは受信、送信アンテナ20の利得を示す。また、ωは、生体活動の最小周波数を示し、ωは生体活動の最大周波数を示す。 In this case, R R represents a distance to the position of the biometric 50 estimated from the received antenna 30, the R T shows the distance to the position of the biological 50 estimated from the transmitting antenna 20. Also, G R represents a gain of the receiving antenna 30, G T is received, indicating the gain of the transmission antenna 20. Further, ω 1 indicates the minimum frequency of biological activity, and ω 2 indicates the maximum frequency of biological activity.

このように、生体活動に対応する周波数に対応する成分のみを抽出することにより生体からの反射電力のみを抽出することができる。生体の姿勢によってアンテナから見た体表面積が異なって見えることから、生体の状態を推定することが可能であると共に、生体が行動をした際に、体動による周波数成分が増えるためRCSが変動し、推定された生体の高さzと散乱断面積σとの軌道をモデル化することによって、生体の動作の推定が可能となる。 In this way, only the reflected power from the living body can be extracted by extracting only the components corresponding to the frequencies corresponding to the biological activity. Since the body surface area seen from the antenna looks different depending on the posture of the living body, it is possible to estimate the state of the living body, and when the living body behaves, the frequency component due to body movement increases, so the RCS fluctuates. By modeling the orbits of the estimated height z of the living body and the scattering cross section σ, it is possible to estimate the movement of the living body.

以下では動作推定の方法について説明する。 The method of motion estimation will be described below.

前述の方法により,生体の散乱断面積σおよび高さzを異なる複数のタイミングで連続して推定すると、σ−z特性の軌跡を観測することができる。この時、生体が静止状態にある場合は、RCS値の変化がほとんどないため、生体が行動しているとき、すなわちRCS値の変動の大きな軌跡点の流れを抽出する。この軌跡点群に対してi−1番目の軌跡点とi番目の軌跡点との間の距離である軌跡点移動量をΔPiとし、i−1番目の軌跡点とi番目の軌跡点との間の成す角を角度パラメータαiとし、それぞれ次の式21及び式22と定義する。 By continuously estimating the scattering cross section σ and height z of the living body at a plurality of different timings by the above-mentioned method, the locus of the σ−z characteristic can be observed. At this time, when the living body is in a resting state, there is almost no change in the RCS value. Let ΔPi be the amount of movement of the locus point, which is the distance between the i-1st locus point and the i-th locus point with respect to this locus point group, and the i-1st locus point and the i-th locus point. The angle formed between them is defined as the angle parameter αi, and is defined as the following equations 21 and 22, respectively.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

次に、角度パラメータの値を用いて軌跡点の移動方向に方向コードを割り当てる。角度に対する方向コードの割り当ては360°を8分割し、各方向に1〜8の数字を割り当てている。この際、垂直、水平方向への移動の際にコードが頻繁に変化しないように、その方向にコードの境を設けないようにしてもよい。動作の速度には個人差があるため、同じ動作であっても動作の速さによって軌跡点の数が異なることによる誤認識を避けるために、動作の速さの違いを考慮した方向コードの正規化を行ってもよい。例えば、方向推定により得られた原方向コード列c(j=1〜jmax)から、軌跡点移動量の総和に対する軌跡点移動量の割合を考慮して、K項からなる正規化コード列を生成する。ここでjmaxは軌跡点数であり、動作時間により異なってくる。正規化コード列C(k=1,2,・・・,K)はi番目の軌跡移動量ΔPi、軌跡点移動量の総和ΔPsumから以下のようにして作成する。 Next, the direction code is assigned to the moving direction of the locus point using the value of the angle parameter. The direction code is assigned to the angle by dividing 360 ° into eight parts and assigning numbers 1 to 8 in each direction. At this time, in order to prevent the cord from changing frequently when moving in the vertical and horizontal directions, the boundary of the cord may not be provided in that direction. Since there are individual differences in the speed of movement, normalization of the direction code that takes into account the difference in speed of movement in order to avoid misrecognition due to the difference in the number of locus points depending on the speed of movement even for the same movement It may be converted. For example, from the original direction code sequence c j (j = 1 to j max ) obtained by direction estimation, a normalized code sequence consisting of K terms is taken into consideration in consideration of the ratio of the locus point movement amount to the total locus point movement amount. To generate. Here, j max is the number of locus points and varies depending on the operating time. The normalized code sequence C k (k = 1, 2, ..., K) is created from the i-th locus movement amount ΔPi and the sum total of the locus point movement amounts ΔP sum as follows.

1)j=1では、次の式23を満たすkの範囲の正規化コード列は、原方向コード列におけるj=1のコードとする。 1) In j = 1, the normalized code string in the range of k satisfying the following equation 23 is the code of j = 1 in the original direction code string.

Figure 0006868846
Figure 0006868846

2)jが2〜jmaxの範囲では、次式24を満たすkの範囲の正規化コード列はそれぞれのjでのコードとする。ただし、ΔPsumは次式25を満たす。 2) In the range where j is 2 to j max , the normalized code string in the range of k that satisfies the following equation 24 is the code for each j. However, ΔP sum satisfies the following equation 25.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

動作の推定は、軌跡点データから作成したK項からなる正規化数列(テストデータ)を複数の動作のそれぞれに対応するK項からなるモデルデータ数列と比較することで行う。モデルデータ数列は、あらかじめ複数の動作のそれぞれを行っている間での動作の測定を複数回行い、複数回の動作の測定によって得られた正規化コード列の各項で最も多い方向コードを、その項でのモデルデータとする。方向コードは環状になっているので、最大差は4である。方向コード番号と実際の差を算出するため、次の式26とする。ただし、δC>4のとき、次の式27とする。 The motion is estimated by comparing the normalized sequence (test data) consisting of K terms created from the locus point data with the model data sequence consisting of K terms corresponding to each of the plurality of motions. In the model data sequence, the motion is measured multiple times while each of the plurality of motions is performed in advance, and the most directional code in each term of the normalized code sequence obtained by the measurement of the multiple motions is used. Use the model data in that section. Since the direction code is circular, the maximum difference is 4. In order to calculate the difference between the direction code number and the actual difference, the following equation 26 is used. However, when δC i > 4, the following equation 27 is used.

Figure 0006868846
Figure 0006868846

Figure 0006868846
Figure 0006868846

ただし、Ctest,iはテストデータ列の第項の要素を示し、Cmodel,iはモデルデータ列の第i項の要素を示す。次に、テストデータとモデルデータとの方向コード差の2乗総和を、次の式28で示す偏差として算出する。例えば、図10に示されるように、テストデータと、モデルデータとを比較することで、偏差を算出する。図10は、測定により得られたテストデータと、モデルコードであるモデルデータとを示す図である。 However, C test, i denotes an element of the i-th term of the test data sequence, C model, i denotes an element of the i-th term of the model data strings. Next, the sum of squares of the direction code differences between the test data and the model data is calculated as the deviation shown by the following equation 28. For example, as shown in FIG. 10, the deviation is calculated by comparing the test data with the model data. FIG. 10 is a diagram showing test data obtained by measurement and model data which is a model code.

Figure 0006868846
Figure 0006868846

そして、偏差が最小となるモデルデータ列に対応する行動を認識結果として出力する。 Then, the action corresponding to the model data string having the minimum deviation is output as the recognition result.

次に、実施の形態におけるセンサー10の動作についてフローチャートを用いて説明する。 Next, the operation of the sensor 10 in the embodiment will be described with reference to a flowchart.

図11は、実施の形態におけるセンサーの動作の一例を示すフローチャートである。 FIG. 11 is a flowchart showing an example of the operation of the sensor according to the embodiment.

センサー10では、送信アンテナ20のN個の送信アンテナ素子21が、生体50が存在しうる所定範囲に対してN個の送信アンテナ素子21を用いてN個の送信信号を送信する(S11)。 In the sensor 10, the N transmitting antenna elements 21 of the transmitting antenna 20 transmit N transmitting signals using the N transmitting antenna elements 21 to a predetermined range in which the living body 50 can exist (S11).

受信アンテナ30のM個の受信アンテナ素子31が送信アンテナ20により送信されたN個の送信信号が生体50により反射された複数の反射信号を含むN個の受信信号を受信する(S12)。 The M receiving antenna elements 31 of the receiving antenna 30 receive the N receiving signals including the plurality of reflected signals reflected by the living body 50 from the N transmitting signals transmitted by the transmitting antenna 20 (S12).

回路40は、M個の受信アンテナ素子31のそれぞれにおいて所定期間で受信されたN個の受信信号のそれぞれから、N個の送信アンテナ素子21それぞれと、M個の受信アンテナ素子31それぞれとの間の伝播特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出する(S13)。 In the circuit 40, from each of the N received signals received in each of the M receiving antenna elements 31 in a predetermined period, between each of the N transmitting antenna elements 21 and each of the M receiving antenna elements 31. The first matrix of N × M is calculated with each complex transfer function showing the propagation characteristics of (S13) as a component (S13).

回路40は、第1行列における所定周波数範囲に対応する第2行列を抽出することで、生体50の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する第2行列を抽出する(S14)。 The circuit 40 extracts a second matrix corresponding to a predetermined frequency range in the first matrix, and thereby corresponds to a component affected by vital activity including at least one of respiration, heartbeat, and body movement of the living body 50. Two matrices are extracted (S14).

回路40は、第2行列を用いて、センサー10に対する生体50の存在する3次元位置を推定する(S15)。 The circuit 40 estimates the three-dimensional position of the living body 50 with respect to the sensor 10 using the second matrix (S15).

回路40は、推定した3次元位置と、送信アンテナ20の位置と、受信アンテナ30の位置と、に基づいて、生体50と送信アンテナ20との距離を示す距離r1、および、生体50と受信アンテナ30との距離を示す距離r2を算出する(S16)。 The circuit 40 has a distance r1 indicating the distance between the living body 50 and the transmitting antenna 20 based on the estimated three-dimensional position, the position of the transmitting antenna 20, and the position of the receiving antenna 30, and the living body 50 and the receiving antenna. The distance r2 indicating the distance from 30 is calculated (S16).

回路40は、第1距離および第2距離を用いて、生体50に対するRCS値を算出する(S17)。 The circuit 40 calculates the RCS value for the living body 50 using the first distance and the second distance (S17).

回路40は、算出されたRCS値と、メモリ41に記憶されている、RCS値および生体50の動作の対応関係を示す情報42と、を用いて、生体50の動作を推定する(S18)。 The circuit 40 estimates the operation of the living body 50 by using the calculated RCS value and the information 42 stored in the memory 41 indicating the correspondence between the RCS value and the operation of the living body 50 (S18).

次に、生体50の動作を推定する推定処理の詳細を説明する。 Next, the details of the estimation process for estimating the operation of the living body 50 will be described.

図12は、推定処理の詳細の一例を示すフローチャートである。 FIG. 12 is a flowchart showing an example of details of the estimation process.

回路40は、推定した3次元位置のうちの鉛直位置、または、算出したRCS値の時間的変化が所定値よりも大きい期間を、生体50が動作している動作期間として抽出する(S21)。 The circuit 40 extracts a period in which the vertical position of the estimated three-dimensional positions or the time change of the calculated RCS value is larger than a predetermined value as an operating period in which the living body 50 is operating (S21).

回路40は、抽出した動作期間における時間的変化であって、推定した3次元位置から得られる鉛直位置、および、算出したRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換し、変換した方向ベクトルを正規化することで方向コードを算出する(S22)。 The circuit 40 converts the time change in the extracted operation period, the vertical position obtained from the estimated three-dimensional position, and the time change of the calculated RCS value into a direction vector using a predetermined method. , The direction code is calculated by normalizing the converted direction vector (S22).

回路40は、算出した方向コードの時系列データと、メモリに記憶されている対応関係を示す情報42とを比較することにより、対応関係を示す情報42において当該時系列データに対応付けられている動作を特定することで、生体50の動作を推定する(S23)。 The circuit 40 is associated with the time-series data in the information 42 indicating the correspondence by comparing the calculated time-series data of the direction code with the information 42 indicating the correspondence stored in the memory. By specifying the movement, the movement of the living body 50 is estimated (S23).

次に、対応関係を示す情報42を取得するためのセンサー10の事前学習の動作について説明する。 Next, the operation of pre-learning of the sensor 10 for acquiring the information 42 indicating the correspondence relationship will be described.

図13は、事前学習処理におけるセンサーの動作の一例を示すフローチャートである。 FIG. 13 is a flowchart showing an example of the operation of the sensor in the pre-learning process.

回路40は、図示しない入力手段により所定の動作を指定するための入力を受け付ける(S31)。これにより、回路40は、所定の期間の間に行われる動作が、受け付けた入力により示される動作であると認識する。 The circuit 40 receives an input for designating a predetermined operation by an input means (not shown) (S31). As a result, the circuit 40 recognizes that the operation performed during the predetermined period is the operation indicated by the received input.

次に、上記で説明したセンサー10の動作のステップS11〜ステップS17およびステップS21、S22と同様の処理が順に行われる。 Next, the same processes as in steps S11 to S17 and steps S21 and S22 of the operation of the sensor 10 described above are sequentially performed.

次に、回路40は、ステップS31で受け付けた入力が示す動作と、算出した方向コードの時系列データとを対応付けることにより得られた対応関係を示す情報を教師データとしてメモリ41に記憶する(S32)。 Next, the circuit 40 stores in the memory 41 as teacher data information indicating the correspondence relationship obtained by associating the operation indicated by the input received in step S31 with the time series data of the calculated direction code (S32). ).

本実施の形態に係るセンサー10によれば、生体50が存在する位置および当該位置における生体の動作の推定を短時間かつ高精度に行うことができる。 According to the sensor 10 according to the present embodiment, it is possible to estimate the position where the living body 50 exists and the movement of the living body at the position in a short time and with high accuracy.

センサー10は、動いている部位を検出することで、生体50の存在を検出する。このため、例えば、これを利用することにより、ヒトが生きている状態で、かつ、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換のどの動作を行ったかを推定できる。これにより、ヒトの生存確認を効果的に行うことができる。また、カメラで撮像した画像を画像解析することなくヒトの生存確認を行うことができるため、ヒトのプライバシーを保護した状態で、ヒトの生存確認を行うことができる。 The sensor 10 detects the presence of the living body 50 by detecting a moving part. Therefore, for example, by using this, a human being is alive and falls, sits on a chair, sits on a floor, stands up from a chair, stands up from the floor, jumps, and changes direction. It is possible to estimate which operation of the throat was performed. This makes it possible to effectively confirm the survival of humans. In addition, since it is possible to confirm the survival of a human without analyzing the image captured by the camera, it is possible to confirm the survival of a human while protecting the privacy of the human.

以上、本発明の一つまたは複数の態様に係るセンサー10について、実施の形態に基づいて説明したが、本発明は、この実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の一つまたは複数の態様の範囲内に含まれてもよい。 Although the sensor 10 according to one or more aspects of the present invention has been described above based on the embodiment, the present invention is not limited to this embodiment. As long as it does not deviate from the gist of the present invention, one or more of the present embodiments may be modified by those skilled in the art, or may be constructed by combining components in different embodiments. It may be included within the scope of the embodiment.

上記実施の形態によれば、センサー10は、対応関係を示す情報42は、鉛直位置およびRCS値が方向コード化されたモデルコードと、動作とが対応付けられた情報であるとしたが、これに限らない。例えば、対応関係を示す情報として、鉛直位置およびRCS値の時間的変化そのものと、動作とが対応付けられた情報を採用してもよい。この場合、回路40の動作推定部450は、方向コード算出部452を有していなくてもよい。 According to the above embodiment, the sensor 10 states that the information 42 indicating the correspondence relationship is the information in which the vertical position and the RCS value are direction-coded and the operation is associated with the model code. Not limited to. For example, as the information indicating the correspondence relationship, the information in which the vertical position and the temporal change of the RCS value itself and the operation are associated with each other may be adopted. In this case, the motion estimation unit 450 of the circuit 40 does not have to have the direction code calculation unit 452.

上記実施の形態によれば、センサー10の回路40は、推定した動作を用いて生体50により次に行われた動作を推定してもよい。つまり、回路40は、第1の動作期間の次の第2の動作期間において、第1の動作期間の終了時の生体50の姿勢を用いて、生体50の動作を推定してもよい。例えば、回路40は、生体50が起立したことを推定した場合、既に生体50が立位の状態であることはあきらかであるので、次に、起立することはあり得ないと判断し、対応関係を示す情報42と比較するときに、対応関係を示す情報42の比較の対象となる動作から起立する動作を除いてもよい。これにより、生体50の動作の推定をより高精度に行うことができる。 According to the above embodiment, the circuit 40 of the sensor 10 may estimate the next operation performed by the living body 50 using the estimated operation. That is, in the second operation period following the first operation period, the circuit 40 may estimate the operation of the living body 50 by using the posture of the living body 50 at the end of the first operating period. For example, when the circuit 40 estimates that the living body 50 has stood up, it is clear that the living body 50 is already in a standing position. Therefore, the circuit 40 determines that it cannot stand up next, and has a correspondence relationship. When comparing with the information 42 indicating the correspondence, the operation of standing up may be excluded from the operation to be compared with the information 42 indicating the correspondence. As a result, the motion of the living body 50 can be estimated with higher accuracy.

上記実施の形態によれば、センサー10は、無線信号を用いて得られたデータから、生体50の動作を推定するとしたが、動作を推定することに限らない。 According to the above embodiment, the sensor 10 estimates the operation of the living body 50 from the data obtained by using the wireless signal, but the operation is not limited to the estimation.

センサー10の回路40は、例えば、生体50の3次元位置の時系列データから水平方向への変動が所定の距離以上の変動であるか否かを判定し、当該変動が所定の距離以上である場合、生体50が水平方向に移動していると推定してもよい。この場合、回路40は、水平方向への変動が所定の閾値より大きい場合、生体50が走っていると推定してもよいし、所定の閾値より小さい場合、生体50が歩行していると推定してもよい。このため、生体50の水平方向への移動を、短時間かつ高精度に推定することができる。 The circuit 40 of the sensor 10 determines, for example, whether or not the fluctuation in the horizontal direction is a fluctuation of a predetermined distance or more from the time series data of the three-dimensional position of the living body 50, and the fluctuation is a predetermined distance or more. In this case, it may be estimated that the living body 50 is moving in the horizontal direction. In this case, the circuit 40 may presume that the living body 50 is running when the horizontal fluctuation is larger than a predetermined threshold value, or presume that the living body 50 is walking when it is smaller than the predetermined threshold value. You may. Therefore, the horizontal movement of the living body 50 can be estimated in a short time and with high accuracy.

また、センサー10の回路40は、例えば、さらに、生体50が水平方向に移動していると推定した場合の3次元位置に含まれる鉛直位置を用いて、生体50の身長を推定してもよい。この場合、回路40は、具体的には、得られた鉛直位置に所定の係数を乗じることで生体50の身長を推定してもよい。上述したように、生体50の鉛直位置としては、例えば、腹部の鉛直位置が得られるため、例えば、所定の係数として1.5〜2.0の範囲の係数を乗じることで生体50の身長を推定することができる。なお、水平方向へ移動している場合、生体50である人は、立っている状態であることを前提として上記のように所定の係数を乗じることで身長を推定している。 Further, the circuit 40 of the sensor 10 may further estimate the height of the living body 50 by using the vertical position included in the three-dimensional position when the living body 50 is estimated to be moving in the horizontal direction. .. In this case, the circuit 40 may specifically estimate the height of the living body 50 by multiplying the obtained vertical position by a predetermined coefficient. As described above, as the vertical position of the living body 50, for example, the vertical position of the abdomen can be obtained. Therefore, for example, the height of the living body 50 can be obtained by multiplying a predetermined coefficient by a coefficient in the range of 1.5 to 2.0. Can be estimated. When moving in the horizontal direction, the person who is the living body 50 estimates his / her height by multiplying it by a predetermined coefficient as described above on the premise that he / she is in a standing state.

これにより、例えば、予め所定範囲に存在しうる複数の生体が分かっており、複数の生体の身長がそれぞれ異なっていれば、複数の生体のうちのどの生体が存在しているか、動作しているかを特定することに利用することができる。例えば、所定範囲が家の中の部屋であり、当該家に住んでいる人が決まっているのであれば、得られた鉛直位置から当該家に住んでいる人の身長の大きさを参照することで、得られた生体50が当該家に住んでいる人のうちの誰であるかを特定することに利用することができる。 As a result, for example, if a plurality of organisms that can exist in a predetermined range are known in advance and the heights of the plurality of organisms are different, which of the plurality of organisms exists or is operating. Can be used to identify. For example, if a predetermined range is a room in a house and the person who lives in the house is decided, refer to the height size of the person who lives in the house from the obtained vertical position. Therefore, it can be used to identify who the obtained living body 50 is among the people living in the house.

また、センサー10の回路40は、例えば、さらに、生体50が水平方向に移動していると推定した場合のRCS値を用いて、生体50の身体のサイズを推定してもよい。これにより、例えば、予め所定範囲に存在しうる複数の生体が分かっており、複数の生体の身体のサイズがそれぞれ異なっていれば、複数の生体のうちのどの生体が存在しているか、動作しているかを特定することに利用することができる。つまり、身長の場合と同様に、誰が存在しているのかを特定することに利用することができる。 Further, the circuit 40 of the sensor 10 may further estimate the size of the body of the living body 50 by using the RCS value when the living body 50 is estimated to be moving in the horizontal direction. As a result, for example, if a plurality of organisms that can exist in a predetermined range are known in advance and the body sizes of the plurality of organisms are different, which organism among the plurality of organisms exists will operate. It can be used to identify if it is. In other words, it can be used to identify who exists, as in the case of height.

本発明は、無線信号を利用した動体の方向や位置を推定するセンサーおよび方法に利用できる。特に、生体と機械を含む動体の方向や位置を測定する測定器、動体の方向や位置に応じた制御を行う家電機器、動体の侵入を検知する監視装置などに搭載される方向推定方法、方向推定方法,位置推定方法および散乱断面積を用いた高々推定装置に利用できる。 The present invention can be used as a sensor and a method for estimating the direction and position of a moving body using a radio signal. In particular, the direction estimation method and direction mounted on measuring instruments that measure the direction and position of moving objects including living bodies and machines, home appliances that control according to the direction and position of moving objects, and monitoring devices that detect the intrusion of moving objects. It can be used for estimation methods, position estimation methods, and at most estimation devices using scattered cross sections.

10 センサー
20 送信アンテナ
21 送信アンテナ素子
30 受信アンテナ
31 受信アンテナ素子
40 回路
41 メモリ
42 対応関係を示す情報
50 生体
301〜304 格子
410 複素伝達関数算出部
420 生体成分算出部
430 位置推定処理部
440 RCS算出部
450 動作推定部
451 動作期間抽出部
452 方向コード算出部
453 動作比較部
10 Sensor 20 Transmitting antenna 21 Transmitting antenna element 30 Receiving antenna 31 Receiving antenna element 40 Circuit 41 Memory 42 Information indicating correspondence 50 Living body 301-304 Grid 410 Complex transfer function calculation unit 420 Biological component calculation unit 430 Position estimation processing unit 440 RCS Calculation unit 450 Operation estimation unit 451 Operation period extraction unit 452 Direction code calculation unit 453 Operation comparison unit

Claims (11)

センサーであって、
生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは3以上の自然数)の送信アンテナ素子を有する送信アンテナと、
前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは3以上の自然数)の受信アンテナ素子を有する受信アンテナと、
回路と、
前記センサーに対する前記生体の存在する鉛直方向における位置である鉛直位置、および、RCS(Radar cross-section)値の時間的変化と、前記生体の動作との対応関係を示す情報を記憶しているメモリと、を備え、
前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記回路は、
前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、
前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、
前記第2行列を用いて、前記センサーに対する前記生体の存在する3次元位置であって、前記鉛直位置を含む3次元位置を推定し、
前記推定した前記3次元位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、
前記第1距離および前記第2距離を用いて、前記生体に対するRCS値を算出し、
前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定する、
センサー。
It ’s a sensor,
A transmitting antenna having N (N is a natural number of 3 or more) transmitting antenna elements, each of which transmits a transmission signal to a predetermined range in which a living body can exist, and a transmitting antenna.
A part of the N transmission signals transmitted by the N transmission antenna elements receives N reception signals including a reflection signal reflected by the living body (M). Is a receiving antenna having a receiving antenna element (3 or more natural numbers) and
Circuit and
A memory that stores information indicating the correspondence between the vertical position, which is the position of the living body in the vertical direction with respect to the sensor, and the temporal change of the RCS (Radar cross-section) value, and the movement of the living body. And with
Of the N transmitting antenna elements, at least three transmitting antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively.
At least three of the M receiving antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively.
The circuit
Propagation characteristics between each of the N transmitting antenna elements and each of the M receiving antenna elements from each of the N received signals received in each of the M receiving antenna elements in a predetermined period. Calculate the first matrix of N × M with each complex transfer function indicating
By extracting the second matrix corresponding to the predetermined frequency range in the first matrix, the second matrix corresponding to the component affected by the vital activity including at least one of the respiration, heartbeat and body movement of the living body. Extract and
Using the second matrix, the three-dimensional position where the living body exists with respect to the sensor, including the vertical position, is estimated.
Based on the estimated three-dimensional position, the position of the transmitting antenna, and the position of the receiving antenna, a first distance indicating the distance between the living body and the transmitting antenna, and the living body and the receiving antenna. Calculate the second distance indicating the distance to and
Using the first distance and the second distance, the RCS value for the living body was calculated.
The movement of the living body is estimated by using the estimated three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory.
sensor.
前記対応関係を示す情報において対応付けられている前記生体の動作は、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換を含み、
前記回路は、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体が、転倒、椅子への着座、床への着座、椅子からの起立、床からの起立、ジャンプ、および、方向転換のいずれの動作を行ったかを推定する
請求項1に記載のセンサー。
The movements of the living body associated with the information indicating the correspondence include falling, sitting on a chair, sitting on the floor, standing up from the chair, standing up from the floor, jumping, and changing direction.
In the circuit, the living body falls over using the estimated three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory. The sensor according to claim 1, which estimates whether the movement of sitting on a chair, sitting on the floor, standing up from a chair, standing up from the floor, jumping, or changing direction is performed.
前記回路は、前記生体の動作の推定において、
前記推定した鉛直位置、または、前記算出したRCS値の時間的変化が所定値よりも大きい期間を、前記生体が動作している動作期間として抽出し、
抽出した前記動作期間における時間的変化であって、前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定する、
請求項1または2に記載のセンサー。
The circuit is used in estimating the movement of the living body.
The period in which the estimated vertical position or the time change of the calculated RCS value is larger than the predetermined value is extracted as the operating period in which the living body is operating.
The extracted temporal change in the operation period, the estimated three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory. To estimate the behavior of the living body,
The sensor according to claim 1 or 2.
前記回路は、時系列で得られた複数の前記鉛直位置、または、複数の前記RCS値から、所定のフィルタを用いて瞬時ノイズ成分を除去することで得られる時系列データを用いて、前記動作期間を抽出する
請求項3に記載のセンサー。
The circuit uses time-series data obtained by removing an instantaneous noise component from a plurality of the vertical positions or a plurality of the RCS values obtained in a time series using a predetermined filter to perform the operation. The sensor according to claim 3, wherein the period is extracted.
前記対応関係を示す情報において前記生体の動作に対応付けられている、前記鉛直位置および前記RCS値は、事前に、前記生体が前記生体の動作のうちの一の動作を前記所定範囲で行ったときに、前記回路により推定された鉛直位置、および、前記回路により算出されたRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換し、変換した前記方向ベクトルを正規化することで得られた方向コードで示されており、
前記回路は、前記生体の動作の推定において、
抽出した前記動作期間における時間的変化であって、前記推定した前記3次元位置から得られる前記鉛直位置、および、前記算出したRCS値の時間的変化を、所定の方法を用いて方向ベクトルに変換し、
変換した前記方向ベクトルを正規化することで方向コードを算出し、
算出した前記方向コードと、前記対応関係を示す情報と、を用いて、前記生体の動作を推定する
請求項3または4に記載のセンサー。
The vertical position and the RCS value, which are associated with the movement of the living body in the information indicating the correspondence, are obtained by the living body performing one of the movements of the living body in the predetermined range in advance. Occasionally, the vertical position estimated by the circuit and the temporal change of the RCS value calculated by the circuit are converted into a direction vector by using a predetermined method, and the converted direction vector is normalized. Indicated by the direction code obtained in
The circuit is used in estimating the movement of the living body.
The extracted temporal change in the operation period, the vertical position obtained from the estimated three-dimensional position, and the temporal change of the calculated RCS value are converted into a direction vector by using a predetermined method. And
The direction code is calculated by normalizing the converted direction vector.
The sensor according to claim 3 or 4, which estimates the movement of the living body by using the calculated direction code and the information indicating the correspondence relationship.
前記回路は、第1の動作期間の次の第2の動作期間において、前記第1の動作期間の終了時の前記生体の姿勢を用いて、前記生体の動作を推定する
請求項1から5のいずれか1項に記載のセンサー。
The circuit of claims 1 to 5 estimates the movement of the living body by using the posture of the living body at the end of the first operating period in the second operating period following the first operating period. The sensor according to any one item.
前記回路は、さらに、前記推定した前記3次元位置における水平方向への変動が所定の距離以上の変動である場合、前記生体が水平方向に移動していると推定する
請求項1から6のいずれか1項に記載のセンサー。
The circuit further presumes that the living body is moving in the horizontal direction when the horizontal fluctuation at the estimated three-dimensional position is a fluctuation of a predetermined distance or more. Or the sensor described in item 1.
前記回路は、さらに、前記生体が水平方向に移動していると推定した場合の前記3次元位置に含まれる前記鉛直位置を用いて、前記生体の身長を推定する
請求項7に記載のセンサー。
The sensor according to claim 7, wherein the circuit further estimates the height of the living body by using the vertical position included in the three-dimensional position when the living body is estimated to be moving in the horizontal direction.
前記回路は、さらに、前記生体が水平方向に移動していると推定した場合の前記RCS値を用いて、前記生体の身体のサイズを推定する
請求項7または8に記載のセンサー。
The sensor according to claim 7 or 8, wherein the circuit further estimates the size of the body of the living body by using the RCS value when it is estimated that the living body is moving in the horizontal direction.
前記所定期間は、前記生体の呼吸、心拍、および、体動の少なくとも1つの周期の略半分である
請求項1から9のいずれか1項に記載のセンサー。
The sensor according to any one of claims 1 to 9, wherein the predetermined period is approximately half of at least one cycle of respiration, heartbeat, and body movement of the living body.
センサーによって生体の動作を推定する方法であって、
前記センサーは、N個(Nは3以上の自然数)の送信アンテナ素子を有する送信アンテナと、M個(Mは3以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、前記センサーに対する生体の存在する鉛直方向における位置である鉛直位置、RCS(Radar cross-section)値、および、前記生体の動作の対応関係を示す情報を記憶しているメモリとを備え、
前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記方法は、
生体が存在しうる所定範囲に対して前記N個の送信アンテナを用いてN個の送信信号を送信し、
送信した前記N個の送信信号の一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号を前記M個の受信アンテナ素子のそれぞれを用いて受信し、
前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、
前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、
前記第2行列を用いて、前記センサーに対する前記生体の存在する3次元位置であって、前記鉛直位置を含む3次元位置を推定し、
前記推定した前記3次元位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、
前記第1距離および前記第2距離を用いて、前記生体に対するRCS値を算出し、
前記推定した前記3次元位置、および、前記算出したRCS値の時間的変化と、前記メモリに記憶されている前記対応関係を示す情報と、を用いて、前記生体の動作を推定する、
方法。
It is a method of estimating the movement of a living body by a sensor.
The sensor is a transmitting antenna having N (N is a natural number of 3 or more) transmitting antenna elements, a receiving antenna having M receiving antenna elements (M is a natural number of 3 or more), a circuit, and the sensor. It is provided with a memory that stores information indicating a vertical position, which is a position in the vertical direction in which a living body exists, an RCS (Radar cross-section) value, and information indicating a correspondence relationship between the movements of the living body.
Of the N transmitting antenna elements, at least three transmitting antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively.
At least three of the M receiving antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively.
The method is
N transmission signals are transmitted using the N transmission antennas in a predetermined range in which a living body can exist.
A part of the transmitted signals of the N transmitted signals received the N received signals including the reflected signal reflected by the living body by using each of the M receiving antenna elements.
Propagation characteristics between each of the N transmitting antenna elements and each of the M receiving antenna elements from each of the N received signals received in each of the M receiving antenna elements in a predetermined period. Calculate the first matrix of N × M with each complex transfer function indicating
By extracting the second matrix corresponding to the predetermined frequency range in the first matrix, the second matrix corresponding to the component affected by the vital activity including at least one of the respiration, heartbeat and body movement of the living body. Extract and
Using the second matrix, the three-dimensional position where the living body exists with respect to the sensor, including the vertical position, is estimated.
Based on the estimated three-dimensional position, the position of the transmitting antenna, and the position of the receiving antenna, a first distance indicating the distance between the living body and the transmitting antenna, and the living body and the receiving antenna. Calculate the second distance indicating the distance to and
Using the first distance and the second distance, the RCS value for the living body was calculated.
The movement of the living body is estimated by using the estimated three-dimensional position, the temporal change of the calculated RCS value, and the information indicating the correspondence relationship stored in the memory.
Method.
JP2017180438A 2017-01-06 2017-09-20 Sensors and methods Active JP6868846B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US15/851,858 US10912493B2 (en) 2017-01-06 2017-12-22 Sensor and method
CN201711454518.2A CN108279413B (en) 2017-01-06 2017-12-28 Sensor and biological action estimation method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017001317 2017-01-06
JP2017001317 2017-01-06

Publications (2)

Publication Number Publication Date
JP2018112540A JP2018112540A (en) 2018-07-19
JP6868846B2 true JP6868846B2 (en) 2021-05-12

Family

ID=62911145

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017180438A Active JP6868846B2 (en) 2017-01-06 2017-09-20 Sensors and methods

Country Status (1)

Country Link
JP (1) JP6868846B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018235755A1 (en) 2017-06-20 2018-12-27 三ツ星ベルト株式会社 V-ribbed belt and method for manufacturing same
JP7054583B2 (en) * 2018-09-18 2022-04-14 株式会社東芝 Wireless power transmission device, wireless power transmission system and wireless power transmission method
JP7429877B2 (en) 2018-12-28 2024-02-09 パナソニックIpマネジメント株式会社 Estimation method and device
JP7357217B2 (en) 2018-12-28 2023-10-06 パナソニックIpマネジメント株式会社 Estimation method, estimation device, and program
CN111381228A (en) * 2018-12-28 2020-07-07 松下知识产权经营株式会社 Estimation method and estimation device
JP7313018B2 (en) * 2019-07-02 2023-07-24 パナソニックIpマネジメント株式会社 sensor

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9019148B1 (en) * 2012-04-24 2015-04-28 Sandia Corporation Remote sensing using MIMO systems
JP6222454B2 (en) * 2013-12-20 2017-11-01 三星電子株式会社Samsung Electronics Co.,Ltd. Biological information measuring apparatus and biological information measuring method

Also Published As

Publication number Publication date
JP2018112540A (en) 2018-07-19

Similar Documents

Publication Publication Date Title
CN108279413B (en) Sensor and biological action estimation method
JP6868846B2 (en) Sensors and methods
JP6839807B2 (en) Estimator and program
US10895631B2 (en) Sensor and estimating method
US11226408B2 (en) Sensor, estimating device, estimating method, and recording medium
Ding et al. WiFi CSI-based human activity recognition using deep recurrent neural network
JP6751923B2 (en) program
JP6900390B2 (en) Methods, devices, servers and systems for detecting and monitoring vital signs
US9568594B2 (en) Human posture feature extraction in personal emergency response systems and methods
US10561358B2 (en) Biometric device and biometric method
CN111142102B (en) Respiratory data calculation method and related equipment
US20210321222A1 (en) Systems and methods for mapping a given environment
Liu et al. An automatic in-home fall detection system using Doppler radar signatures
US20200341115A1 (en) Subject identification in behavioral sensing systems
JP6887089B2 (en) Identification device and identification method
CN111381227B (en) Speculation method and speculation device
CN111381228A (en) Estimation method and estimation device
US11796304B2 (en) System and method for tracking a shape
JP7174921B2 (en) Sensor, estimation device, estimation method, and program
JP7429877B2 (en) Estimation method and device
CN113033407B (en) Non-contact type fitness monitoring method using intelligent sound box
Hino et al. Activity recognition using array antenna
WO2023276945A1 (en) Estimation device and estimation method
US20230371820A1 (en) Contactless Seismocardiography
Civan Through the wall human movement detection with software defined radio

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200403

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210310

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210323

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210329

R151 Written notification of patent or utility model registration

Ref document number: 6868846

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151