[go: up one dir, main page]

JP6667145B2 - Sensor and estimation method - Google Patents

Sensor and estimation method Download PDF

Info

Publication number
JP6667145B2
JP6667145B2 JP2017063507A JP2017063507A JP6667145B2 JP 6667145 B2 JP6667145 B2 JP 6667145B2 JP 2017063507 A JP2017063507 A JP 2017063507A JP 2017063507 A JP2017063507 A JP 2017063507A JP 6667145 B2 JP6667145 B2 JP 6667145B2
Authority
JP
Japan
Prior art keywords
living body
antenna elements
matrix
distance
sensor
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
JP2017063507A
Other languages
Japanese (ja)
Other versions
JP2018008021A (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 CN201710416994.9A priority Critical patent/CN107576957B/en
Priority to US15/631,367 priority patent/US10473756B2/en
Publication of JP2018008021A publication Critical patent/JP2018008021A/en
Priority to US16/598,628 priority patent/US10895631B2/en
Application granted granted Critical
Publication of JP6667145B2 publication Critical patent/JP6667145B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Description

本開示は、無線信号を利用することで生体の姿勢を推定するセンサーおよび推定方法に関する。   The present disclosure relates to a sensor and an estimation method for estimating a posture of a living body by using a wireless signal.

人物の位置などを知る方法として、無線信号を利用する方法が検討されている(例えば、特許文献1〜3参照)。特許文献1にはドップラーセンサを用いた生体検出の方法、特許文献2にはドップラーセンサとフィルタとを用いた人の動作や生体情報の検知方法が開示されている。特許文献3には、フーリエ変換を用いてドップラーシフトを含む成分を解析することで検出対象となる人物の位置や状態を知ることができることが開示されている。また、特許文献4には、複数アンテナのチャネル情報、各種センサ情報を元に機械学習により生体の位置や状態を推定する方法が開示されており、特許文献5には複数アンテナ、超音波レーダ、複数アンテナによる生体の状態推定方法が開示されている。   As a method of knowing the position of a person or the like, a method using a radio signal is being studied (for example, see Patent Documents 1 to 3). Patent Document 1 discloses a method of detecting a living body using a Doppler sensor, and Patent Document 2 discloses a method of detecting human motion and biological information using a Doppler sensor and a filter. Patent Document 3 discloses that the position and state of a person to be detected can be known by analyzing a component including a Doppler shift using Fourier transform. Patent Document 4 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, and Patent Document 5 discloses a method of estimating a plurality of antennas, an ultrasonic radar, A living body state estimation method using a plurality of antennas is disclosed.

特表2014−512526号公報JP 2014-512526 A 国際公開第2014/141519号International Publication No. 2014/141519 特開2015−117972号公報JP-A-2005-117792 特開2014−190724号公報JP 2014-190724 A 特開2005−292129号公報JP 2005-292129 A 特開2001−159678号公報JP 2001-159678 A

しかしながら、無線信号を利用することで、生体の姿勢を推定する精度を向上させるには、さらなる改善が求められている。   However, in order to improve the accuracy of estimating the posture of a living body by using a wireless signal, further improvement is required.

上記目的を達成するために、本発明の一形態に係るセンサーは、センサーであって、生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは2以上の自然数)の送信アンテナ素子を有する送信アンテナと、前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは2以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、メモリと、を備え、前記回路は、前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、前記第2行列を用いて、前記センサーに対する前記生体の存在する位置を推定し、前記推定した位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、前記第1距離および前記第2距離を用いて、前記生体に対するRCS(Rader cross-section)値を算出し、前記算出されたRCS値と、前記メモリに記憶されている、RCS値および前記生体の姿勢の対応関係を示す情報と、を用いて、前記生体の姿勢を推定する。   In order to achieve the above object, a sensor according to one embodiment of the present invention is a sensor, which transmits N transmission signals to a predetermined range where a living body can exist (N is a natural number of 2 or more). A transmission antenna having a plurality of transmission antenna elements, and N reception elements including a reflection signal obtained by reflecting a part of transmission signals of the N transmission signals transmitted by the N transmission antenna elements by the living body. A receiving antenna having M (M is a natural number of 2 or more) receiving antenna elements each receiving a signal; a circuit; and a memory, wherein the circuit is provided in each of the M receiving antenna elements. From each of the N received signals received during the period, a complex transmission signal indicating a propagation characteristic between each of the N transmitting antenna elements and each of the M receiving antenna elements. By calculating a first matrix of N × M having a function as a component and extracting a second matrix corresponding to a predetermined frequency range in the first matrix, at least one of respiration, heartbeat, and body motion of the living body Extracting the second matrix corresponding to the component affected by the vital activity, including: estimating a position of the living body with respect to the sensor using the second matrix; A first distance indicating a distance between the living body and the transmitting antenna and a second distance indicating a distance between the living body and the receiving antenna are calculated based on the position of the antenna and the position of the receiving antenna. Calculating an RCS (Rader cross-section) value for the living body using the first distance and the second distance, and calculating the calculated RCS value and the RCS value stored in the memory; With, and information indicating the correspondence between the orientation of the serial vivo, to estimate the posture of the living body.

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

図1は、実施の形態1におけるセンサーの構成の一例を示すブロック図である。FIG. 1 is a block diagram illustrating an example of a configuration of the sensor according to the first embodiment. 図2は、実施の形態1における回路およびメモリの機能的な構成を示すブロック図である。FIG. 2 is a block diagram showing a functional configuration of a circuit and a memory according to the first embodiment. 図3は、実施の形態1における対応関係を示す情報の一例を示す。FIG. 3 shows an example of information indicating the correspondence in the first embodiment. 図4は、実施の形態1におけるセンサーの動作の一例を示すフローチャートである。FIG. 4 is a flowchart illustrating an example of the operation of the sensor according to the first embodiment. 図5は、実施の形態2におけるセンサーの構成の一例を示すブロック図である。FIG. 5 is a block diagram illustrating an example of a configuration of a sensor according to the second embodiment. 図6は実施の形態2における回路およびメモリの機能的な構成を示すブロック図である。FIG. 6 is a block diagram showing a functional configuration of a circuit and a memory according to the second embodiment. 図7は、実施の形態2における対応関係を示す情報の一例を示す。FIG. 7 shows an example of information indicating the correspondence in the second embodiment. 図8は、実施の形態2のセンサーの効果を確認するために実施した実験の概要を示す図である。FIG. 8 is a diagram illustrating an outline of an experiment performed to confirm the effect of the sensor according to the second embodiment. 図9は、図8で示した実験系を用いた実験結果を示す図である。FIG. 9 is a diagram showing an experimental result using the experimental system shown in FIG. 図10は実験結果から得られた第1RCS範囲〜第4RCS範囲の具体例、および、第1高さ範囲〜第4高さ範囲の具体例を示す図である。FIG. 10 is a diagram showing a specific example of the first to fourth RCS ranges and a specific example of the first to fourth height ranges obtained from the experimental results.

(本発明の基礎となった知見)
無線信号を利用した生体の状態推定に関する従来技術について、発明者らは詳細な検討を行った。その結果、特許文献1および特許文献2の方法では、人物の在、不在は検知可能だが、人物の存在する方向、位置、大きさ、姿勢などを検出することはできない問題があることがわかった。
(Knowledge underlying the present invention)
The present inventors have conducted detailed studies on a conventional technique relating to estimation of a living body state using a radio signal. As a result, it was found that the methods of Patent Documents 1 and 2 can detect the presence or absence of a person, but cannot detect the direction, position, size, posture, and the like of the person. .

また、特許文献3の方法では、人物などの生体が存在する方向や生体が存在する位置を短時間かつ高精度に検出することは困難であるという問題があることがわかった。なぜなら、生体活動由来のドップラー効果による周波数変化は極めて小さく、フーリエ変換によってこの周波数変化を観測するためには、生体が静止した状態で長時間(例えば数十秒)の観測が必須であるからである。また、一般的に、生体は数十秒にわたって同じ姿勢や位置を継続することは少ないからである。   Further, it has been found that the method of Patent Document 3 has a problem that it is difficult to detect a direction in which a living body such as a person exists or a position where the living body exists in a short time and with high accuracy. This is because the frequency change due to the Doppler effect due to biological activity is extremely small, and in order to observe this frequency change by Fourier transform, it is necessary to observe for a long time (for example, several tens of seconds) while the living body is stationary. is there. Further, in general, the living body rarely keeps the same posture and position for several tens of seconds.

さらに、特許文献4では使用者ごとに機械学習をしないといけないという課題があり、また、特許文献5は、天井の広範囲に複数の超音波アンテナを設置する取り付け課題、コスト課題があることがわかった。   Furthermore, Patent Literature 4 has a problem that machine learning must be performed for each user. Patent Literature 5 has a mounting problem and a cost problem of installing a plurality of ultrasonic antennas over a wide area of a ceiling. Was.

発明者らは、以上の課題に対して研究を重ねた結果、異なる位置に置かれたアンテナ素子を含む送信アンテナから送信され、生体によって反射された反射信号の伝搬特性と散乱断面積とを用いることにより、当該生体が存在する方向、位置、大きさ、姿勢などの推定を短時間かつ高精度に行うことが可能であることを見出し、本開示に至った。   The inventors have conducted research on the above problems, and as a result, use the propagation characteristics and the scattering cross section of the reflected signal transmitted from the transmitting antenna including the antenna element placed at different positions and reflected by the living body As a result, it has been found that the direction, position, size, posture, and the like of the living body can be estimated in a short time and with high accuracy, and the present disclosure has been made.

すなわち、本発明の一態様に係るセンサーは、センサーであって、生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは2以上の自然数)の送信アンテナ素子を有する送信アンテナと、前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは2以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、メモリと、を備え、前記回路は、前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、前記第2行列を用いて、前記センサーに対する前記生体の存在する位置を推定し、前記推定した位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、前記第1距離および前記第2距離を用いて、前記生体に対するRCS(Rader cross-section)値を算出し、前記算出されたRCS値と、前記メモリに記憶されている、RCS値および前記生体の姿勢の対応関係を示す情報と、を用いて、前記生体の姿勢を推定する。   That is, the sensor according to one embodiment of the present invention is a sensor and includes N (N is a natural number of 2 or more) transmitting antenna elements each transmitting a transmission signal to a predetermined range where a living body can exist. A transmitting antenna and N received signals including a reflected signal of a part of the N transmitted signals transmitted by the N transmitting antenna elements, the reflected signals being reflected by the living body, respectively. A reception antenna having M reception antenna elements (M is a natural number of 2 or more), a circuit, and a memory, wherein the circuit receives the reception antenna elements in each of the M reception antenna elements for a predetermined period. From each of the N received signals, a complex transfer function indicating a propagation characteristic between each of the N transmit antenna elements and each of the M receive antenna elements is defined as a component. , N × M first matrix, and extracting a second matrix corresponding to a predetermined frequency range in the first matrix, the vital activity including at least one of respiration, heartbeat, and body movement of the living body Extracting the second matrix corresponding to the affected component, using the second matrix to estimate the position of the living body with respect to the sensor, the estimated position, the position of the transmitting antenna, Calculating a first distance indicating a distance between the living body and the transmitting antenna, and a second distance indicating a distance between the living body and the receiving antenna, based on the position of the receiving antenna; And calculating a RCS (Rader cross-section) value for the living body using the second distance, and calculating the calculated RCS value and the RCS value and the posture of the living body stored in the memory. Using the information indicating the response relationship, and to estimate the posture of the living body.

このため、生体が存在する位置および当該位置における生体の姿勢の推定を短時間かつ高精度に行うことができる。   Therefore, the position where the living body exists and the posture of the living body at the position can be estimated in a short time and with high accuracy.

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

このため、生体が存在する位置および当該位置における生体の姿勢の推定を効果的に行うことができる。   Therefore, it is possible to effectively estimate the position where the living body exists and the posture of the living body at the position.

また、前記回路は、前記生体が、前記送信アンテナおよび前記受信アンテナの並び方向に垂直な方向に対して正対している姿勢であるか否かを推定してもよい。   Further, the circuit may estimate whether or not the living body is in a posture facing the direction perpendicular to the direction in which the transmitting antenna and the receiving antenna are arranged.

また、前記Nは3以上の自然数であり、前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記Mは3以上の自然数であり、前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、前記対応関係を示す情報は、前記センサーに対する前記生体の存在する鉛直方向における位置である鉛直位置、RCS値、および、前記生体の姿勢の対応関係を示し、前記対応関係を示す情報において対応付けられている前記生体の姿勢は、直立、椅座、胡坐、および、仰臥を含み、前記回路は、前記第2行列を用いて、前記鉛直位置を含む3次元位置を推定し、前記推定した前記3次元位置と、前記算出されたRCS値と、前記メモリに記憶されている、前記対応関係を示す情報と、を用いて、前記生体が、直立、椅座、胡坐、および、仰臥のいずれの姿勢であるかを推定してもよい。   The N is a natural number of 3 or more, and at least three of the N transmission antenna elements are arranged at different positions in the vertical direction and the horizontal direction, respectively, and the M is 3 or more. A natural number, at least three receiving antenna elements of the M receiving antenna elements are respectively arranged at different positions in a vertical direction and a horizontal direction, and the information indicating the correspondence is the information of the living body with respect to the sensor. The vertical position, which is the position in the vertical direction, the RCS value, and the corresponding relationship between the posture of the living body are shown, and the posture of the living body associated with the information indicating the corresponding relationship is upright, chair, and crossed. , And supine, wherein the circuit estimates a three-dimensional position including the vertical position using the second matrix, and the estimated three-dimensional position , Using the calculated RCS value and the information indicating the correspondence stored in the memory to determine whether the living body is in an upright, chair, cross, or supine position. May be estimated.

このため、生体が存在する3次元位置および当該3次元位置における生体の姿勢の推定を短時間かつ高精度に行うことができる。   Therefore, the three-dimensional position where the living body exists and the posture of the living body at the three-dimensional position can be estimated in a short time and with high accuracy.

なお、本発明は、装置として実現するだけでなく、このような装置が備える処理手段を備える集積回路として実現したり、その装置を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを示す情報、データまたは信号として実現したりすることもできる。そして、それらプログラム、情報、データおよび信号は、CD−ROM等の記録媒体やインターネット等の通信媒体を介して配信してもよい。   It should be noted that the present invention can be realized not only as an apparatus, but also as an integrated circuit including processing means provided in such an apparatus, or as a method in which processing means constituting the apparatus are implemented as steps. Can be realized as a program to be executed by a computer, or can be realized as information, data or a signal indicating the program. These 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 each of the embodiments described below shows a preferred specific example of the present invention. Numerical values, shapes, materials, constituent elements, arrangement positions and connection forms of constituent elements, steps, order of steps, and the like shown in the following embodiments are merely examples, and do not limit the present invention. In addition, among the components in the following embodiments, components not described in the independent claims that represent the highest concept of the present invention are described as arbitrary components that constitute a more preferable embodiment. In the specification and the drawings, components having substantially the same functional configuration are denoted by the same reference numerals, and redundant description is omitted.

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

図1に示すように、センサー10は、送信アンテナ20、受信アンテナ30、回路40およびメモリ41を備える。センサー10は、ヒト等の生体50に対して送信アンテナ20よりマイクロ波を発射し、受信アンテナ30にて生体50で反射された反射波を受信する。ここで、送信アンテナ20に対して任意に設定された第1基準方向と、送信アンテナ20から生体50への方向である第1生体方向とのなす角をθとする。同様に、受信アンテナ30に対して任意に設定された第2基準方向と、受信アンテナ30から生体50への方向である第2生体方向とのなす角をθとする。なお、第1基準方向、第1生体方向、第2基準方向、および、第2生体方向は、水平面上の方向である。 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 a microwave from the transmitting antenna 20 to a living body 50 such as a human, and receives a reflected wave reflected by the living body 50 at the receiving antenna 30. Here, the first reference direction is arbitrarily set to the transmission antenna 20, and the angle theta T between the first biological direction is the direction from the transmitting antenna 20 to the living body 50. Similarly, a second reference direction which is arbitrarily set to the receiving antenna 30, and the angle theta R from the receiving antenna 30 and the second living direction is the direction of the living body 50. The first reference direction, the first living body direction, the second reference direction, and the second living body direction are directions on a horizontal plane.

送信アンテナ20は、N個(Nは2以上の自然数)の送信アンテナ素子21を有する。送信アンテナ20は、N個の送信アンテナ素子21が水平面上の第1所定方向に並んで配置されることで構成されるアレーアンテナを有する。N個の送信アンテナ素子21のそれぞれは、生体が存在しうる所定範囲に対して送信信号を送信する。つまり、送信アンテナ20は、異なるN箇所の位置からN個の送信信号を所定範囲に対して送信する。なお、生体が存在しうる所定範囲とは、センサー10が生体の存在を検知する検知範囲である。   The transmission antenna 20 has N (N is a natural number of 2 or more) transmission antenna elements 21. The transmission antenna 20 has an array antenna configured by arranging N transmission antenna elements 21 in a first predetermined direction on a horizontal plane. Each of the N transmission antenna elements 21 transmits a transmission signal to a predetermined range where a living body can exist. That is, the transmission antenna 20 transmits N transmission signals from different N positions to the predetermined range. Note that the predetermined range in which a living body can exist is a detection range in which the sensor 10 detects the presence of a 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 a microwave as a transmission signal to a living body 50 such as a human. The N transmission antenna elements 21 may transmit a signal on which modulation processing different for each transmission antenna element 21 has been performed as a transmission signal. Further, each of the N transmitting antenna elements 21 may sequentially switch a modulated signal or a non-modulated signal and transmit the signal. The modulation process may be performed by the transmission antenna 20. As described above, by making the transmission signals transmitted from the N transmission antenna elements 21 different transmission signals for each of the N transmission antenna elements 21, the transmission signal transmitted by the reception antenna 30 is transmitted. The antenna element 21 can be specified. As described above, the transmission antenna 20 may include a circuit for performing the modulation process.

受信アンテナ30は、M個(Mは2以上の自然数)の受信アンテナ素子31を有する。受信アンテナ30は、M個の受信アンテナ素子31が水平面上の第2所定方向に並んで配置されることで構成されるアレーアンテナを有する。M個の受信アンテナ素子31のそれぞれは、N個の送信信号のうち生体50により反射された信号である反射信号を含むN個の受信信号を受信する。受信アンテナ30は、マイクロ波からなる受信信号を周波数変換し、低周波数信号に変換する。受信アンテナ30は、低周波数信号に変換することにより得られた信号を回路40に出力する。つまり、受信アンテナ30は、受信信号を処理するための回路を含んでいてもよい。   The receiving antenna 30 has M (M is a natural number of 2 or more) receiving antenna elements 31. The receiving antenna 30 has an array antenna configured by arranging M receiving antenna elements 31 in a second predetermined direction on a horizontal plane. Each of the M reception antenna elements 31 receives N reception signals including a reflection signal that is a signal reflected by the living body 50 among the N transmission signals. The receiving antenna 30 converts the frequency of a received signal composed of microwaves and converts it into a low-frequency signal. The receiving antenna 30 outputs a signal obtained by converting the signal to 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 performs various processes for operating the sensor 10. The circuit 40 includes, for example, a processor that executes a control program, and a volatile storage area (main storage device) used as a work area used when executing the control program. The volatile storage area is, for example, a RAM (Random Access Memory). Note that 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), for example, a ROM (Read Only Memory), a flash memory, a HDD (Hard Disk Drive), or the like. The memory 41 stores, for example, information used for various processes for operating the sensor 10.

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

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

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

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

生体成分算出部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 a complex transfer function matrix component obtained from a received signal passing through the living body 50 and a complex transfer function matrix component obtained from a received signal not passing through the living body 50. The component that has passed through the living body 50 is a component that fluctuates over time due to life activity. Therefore, the components passing through the living body 50 are, for example, components other than the direct current from the components obtained by performing Fourier transform on the components of the complex transfer function matrix in the time direction when the components other than the living body 50 are stationary. Can be extracted by taking out Further, the components that have 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 is not within the predetermined range exceeds a predetermined threshold. As described above, the biological component calculating 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 Is extracted. The predetermined frequency range is, for example, a frequency derived from vital activity including at least one of the above-described respiration, heartbeat, and body motion of a living body. The predetermined frequency range is, for example, a frequency in a range from 0.1 Hz to 3 Hz. As a result, it is possible to extract the vital activity of the part of the living body 50 due to the heart, lung, diaphragm, and built-in movement, or the biological component affected by the vital activity of the hand, foot, or the like. The part of the living body 50 due to the heart, lung, diaphragm, and built-in movement is, for example, a person's soul.

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

生体成分算出部420で算出された生体成分は、位置推定処理部430に出力される。位置推定処理部430は、算出された生体成分を用いて生体の位置推定を行う。つまり、位置推定処理部430は、第2行列を用いて、センサー10に対する生体50の存在する位置を推定する。位置推定には、送信アンテナ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 a living body using the calculated biological component. That is, the position estimation processing unit 430 estimates the position of the living body 50 with respect to the sensor 10 using the second matrix. For the position estimation, both the angle of departure θ T from the transmitting antenna 20 and the angle of arrival θ R to the receiving antenna 30 are estimated, and the living body 50 is triangulated from the estimated departure angle θ T and angle of arrival θ R. Estimate the position of.

RCS算出部440は、生体成分と推定された位置とを用いて散乱断面積(RCS:Radar Cross Section)を算出する。RCS算出部440は、具体的には、散乱断面積を計算するために、推定された位置と、送信アンテナ20の位置と、受信アンテナ30の位置と、に基づいて、生体50と送信アンテナ20との距離を示す距離RT、および、生体50と受信アンテナ30との距離を示す距離RRを算出する。RCS算出部440は、算出した距離RTおよび距離RRから伝搬距離を算出し、算出した伝播距離と生体成分の強度とを用いてRCSを算出する。なお、送信アンテナ20の位置と、受信アンテナ30の位置とは、メモリ41に予め記憶されていてもよい。   The RCS calculation unit 440 calculates a scattering cross section (RCS: Radar Cross Section) using the biological component and the estimated position. Specifically, the RCS calculation unit 440 calculates the scattering cross section based on the estimated position, the position of the transmitting antenna 20, and the position of the receiving antenna 30, based on the position of the living body 50 and the transmitting antenna 20. Is calculated, and a distance RR indicating the distance between the living body 50 and the receiving antenna 30 is 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 strength of the biological component. Note that 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は、RCS算出部440により算出されたRCS値と、メモリ41に記憶されている、RCS値および生体50の姿勢の対応関係を示す情報42と、を用いて、生体50の姿勢を推定する。なお、メモリ41に記憶されているRCS値および生体50の姿勢の対応関係を示す情報42とは、図3に示すように、仰臥、胡坐、椅座および直立で示される各姿勢に予め対応付けられたRCS値の範囲を示す情報である。なお、仰臥は、仰向けの姿勢を示し、椅座は、椅子に座っている姿勢を示す。   The posture estimating unit 450 uses the RCS value calculated by the RCS calculating unit 440 and the information 42 indicating the correspondence between the RCS value and the posture of the living body 50 stored in the memory 41 to determine the posture of the living body 50. Is estimated. As shown in FIG. 3, the information 42 indicating the correspondence between the RCS value and the posture of the living body 50 stored in the memory 41 is previously associated with each of the postures of supine, crossed, chair, and upright. Information indicating the range of the assigned RCS value. In addition, the supine indicates a posture on the back, and the chair indicates a posture sitting on the chair.

例えば、仰臥は、第1RCS範囲と対応付けられており、胡坐は、第2RCS範囲と対応付けられており、椅座は、第3RCS範囲と対応付けられており、直立は、第4RCS範囲と対応付けられている。なお、第1RCS範囲〜第4RCS範囲は、それぞれ、異なるRCS値の範囲である。   For example, the supine is associated with the first RCS range, the cross is associated with the second RCS range, the chair is associated with the third RCS range, and the upright is associated with the fourth RCS range. It is attached. The first RCS range to the fourth RCS range are ranges of different RCS values.

次に実施の形態1のセンサー10の動作原理の詳細を、数式を用いて説明する。なおここでは、フーリエ変換を用いて生体成分を抽出する方法について示す。ここで説明する処理は、回路40により行われる。送信アンテナ20と受信アンテナ30との間の複素伝達関数行列を、

Figure 0006667145
と定義する。ここでtは時刻を表す。式1の各成分をフーリエ変換すると、
Figure 0006667145
のような周波数応答行列が得られる。ここでfは周波数を表しており、周波数応答行列の各成分は複素数である。この周波数応答行列には、生体50を経由する伝搬成分と、生体50以外を経由する伝搬成分との両方が含まれている。生体以外が静止していると考えられる場合、周波数応答行列の直流成分、すなわちG(0)は生体以外の伝搬成分を主として含んでいるものと考えられる。これは、生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動によってドップラーシフトが発生するため、生体を経由する成分はf=0以外に含まれると考えられるからである。さらに、生体の呼吸または心拍の周波数およびその高調波を考えると、f<3〔Hz〕の範囲に生体由来の成分が多く存在するものと考えられる。したがって、例えば0〔Hz〕<f<3〔Hz〕の所定周波数範囲のG(f)を取り出せば、生体成分を効果的に抽出することができる。 Next, the operation principle of the sensor 10 according to the first embodiment will be described in detail using mathematical expressions. Here, a method of extracting a biological component using Fourier transform will be described. The processing described here is performed by the circuit 40. The complex transfer function matrix between the transmitting antenna 20 and the receiving antenna 30 is
Figure 0006667145
Is defined. Here, t represents time. When each component of Equation 1 is Fourier transformed,
Figure 0006667145
The following frequency response matrix is obtained. Here, f represents a frequency, and each component of the frequency response matrix is a complex number. This frequency response matrix includes both a propagation component passing through the living body 50 and a propagation component passing through other than the living body 50. When it is considered that the body other than the living body is stationary, it is considered that the DC component of the frequency response matrix, that is, G (0) mainly includes a propagation component other than the living body. This is because a Doppler shift occurs due to a vital activity including at least one of respiration, heartbeat, and body motion of a living body, and thus components passing through the living body are considered to be included in other than f = 0. Furthermore, considering the frequency of the respiration or heartbeat of the living body and its harmonics, it is considered that there are many components derived from the living body in the range of f <3 [Hz]. Therefore, for example, if G (f) in a predetermined frequency range of 0 [Hz] <f <3 [Hz] is extracted, a biological component can be effectively extracted.

次に生体成分G(f)を用いた生体位置推定方法について説明する。生体成分行列G(f)を、

Figure 0006667145
のようにベクトル形式に並び替える。これを生体成分ベクトルと定義する。ここで{・}は転置を表す。生体成分ベクトルg(f)から相関行列を、
Figure 0006667145
によって計算する。ここで、{・}は複素共役転置を表す。さらにRは0〔Hz〕<f<3〔Hz〕の間で平均化されている。この平均化によって、後述の位置推定精度が向上することが知られている。次に、算出した相関行列Rを固有値分解することで、当該相関行列Rの固有ベクトルUとその複素共役転置ベクトルUを算出する。 Next, a living body position estimation method using the living body component G (f) will be described. The biological component matrix G (f) is
Figure 0006667145
Sort into vector format like This is defined as a biological component vector. Here, {·} T represents transposition. A correlation matrix is obtained from the biological component vector g (f),
Figure 0006667145
Calculate by Here, {·} H represents a complex conjugate transpose. Further, R is averaged between 0 [Hz] <f <3 [Hz]. It is known that this averaging improves the position estimation accuracy described later. Then, the calculated correlation matrix R by eigenvalue decomposition to calculate the eigenvector U and its complex conjugate transposed vector U H of the correlation matrix R.

Figure 0006667145
Figure 0006667145

なお、式5における、固有ベクトルは、以下に示す式6で表される。   Note that the eigenvector in Expression 5 is expressed by Expression 6 shown below.

Figure 0006667145
Figure 0006667145

ここで、uはi列目の固有ベクトルを表しており、要素数はNMである。Dは対角要素が固有値である対角行列であり、

Figure 0006667145
と表される。ここで、diag[・]とは[・]内の要素を対角項に持つ対角行列を表す。回路40は、以上の情報を用いて、検出対象となる生体50の位置推定を行う。ここでは一例としてMUSIC法に基づく位置推定方法について説明する。MUSIC法ではステアリングベクトルと呼ばれる方向ベクトルと式6で示した固有ベクトルとを用いることによって方向や位置を推定する方法である。式3で示した生体成分ベクトルは、本来のM×N行列を変形して得られる。生体50の位置推定を行うためには、これに対応してステアリングベクトルを定義する必要がある。送信アンテナ20の第1基準方向から出発角θの第1生体方向のステアリングベクトルと、受信アンテナ30の第2基準方向から到来角θの第2生体方向のステアリングベクトルとは、それぞれ、
Figure 0006667145
Figure 0006667145
と表される。ここで、kは波数を表し、dはアンテナアレーの各アンテナ素子の素子間隔を表す。なお、本実施の形態では素子間隔が一定のリニアアレーアンテナを想定している。dは、例えば、送信アンテナ20においては、複数の送信アンテナ素子21のうち互いに隣接する2つの送信アンテナ素子21の間隔を表す。また、dは、受信アンテナ30においては、複数の受信アンテナ素子31のうち互いに隣接する2つの受信アンテナ素子31の間隔を表す。これらのステアリングベクトルのクロネッカ積を求めると、
Figure 0006667145
となる。ここで
Figure 0006667145
はクロネッカ積を表す演算子である。a(θ,θ)は、MN×1の要素を持つベクトルであり、出発角θと到来角θとの2変数を持つ関数となる。以降、a(θ,θ)をステアリングベクトルと定義する。検出範囲内に存在する生体の数をLとすると、評価関数
Figure 0006667145
によって生体位置を特定する。ここで、式11の評価関数はMUSICスペクトラムと呼ばれており、送信アンテナ20および受信アンテナ30のそれぞれから検出対象へ向かう方向の組み合わせ(θ,θ)において極大を取る。極大に対応するθおよびθから、三角法を用いて検出対象である生体50の位置を特定できる。ここで、Lは検出対象の数を表す。つまり、固有ベクトルの数MNは検出対象の数Lより大きい必要がある。 Here, ui represents the eigenvector of the i- th column, and the number of elements is NM. D is a diagonal matrix whose diagonal elements are eigenvalues,
Figure 0006667145
It is expressed as Here, diag [•] represents a diagonal matrix having elements in [•] as diagonal terms. The circuit 40 estimates the position of the living body 50 to be detected using the above information. Here, a position estimation method based on the MUSIC method will be described as an example. The MUSIC method is a method of estimating a direction or a position by using a direction vector called a steering vector and an eigenvector shown in Expression 6. The biological component vector shown in Expression 3 is obtained by modifying the original M × N matrix. In order to estimate the position of the living body 50, it is necessary to define a steering vector corresponding to this. The steering vector in the first living body direction from the first reference direction of the transmitting antenna 20 to the departure angle θ T and the steering vector in the second living body direction to the arrival angle θ R from the second reference direction to the receiving antenna 30 are respectively
Figure 0006667145
Figure 0006667145
It is expressed as Here, k represents the wave number, and d represents the element interval of each antenna element of the antenna array. In this embodiment, a linear array antenna having a constant element spacing is assumed. For example, in the transmission antenna 20, d represents the interval between two adjacent transmission antenna elements 21 among the plurality of transmission antenna elements 21. In the receiving antenna 30, d represents the interval between two adjacent receiving antenna elements 31 among the plurality of receiving antenna elements 31. When calculating the Kronecker product of these steering vectors,
Figure 0006667145
Becomes here
Figure 0006667145
Is an operator representing the Kronecker product. a (θ T , θ R ) is a vector having an element of MN × 1, and is a function having two variables of a departure angle θ T and an arrival angle θ R. Hereinafter, a (θ T , θ R ) is defined as a steering vector. Assuming that the number of living organisms present in the detection range is L, the evaluation function
Figure 0006667145
The body position is specified by Here, the evaluation function of Expression 11 is called a MUSIC spectrum, and has a maximum in a combination (θ T , θ R ) of the direction from each of the transmitting antenna 20 and the receiving antenna 30 to the detection target. From θ T and θ R corresponding to the maximum, the position of the living body 50 to be detected can be specified using trigonometry. Here, L represents the number of detection targets. That is, the number MN of eigenvectors needs to be larger than the number L of detection targets.

更に式12よりRCS値を求め、生体50の位置およびRCS値により生体50の姿勢を推定する。前述の抽出する周波数範囲をf1〜f2(f1<f2)とすると、生体から反射され観測されるチャネル成分から電力の伝達係数を求めると、

Figure 0006667145
と計算することができる。ここでρijは行列
Figure 0006667145
の(i,j)番目の要素を表している。一方、j番目の送信アンテナ素子21から生体50を経由してi番目の受信アンテナ素子31に到達する電力は、
Figure 0006667145
と表される。ここで、Ptは、送信電力を表す。なお、全ての送信アンテナ素子21から等しい電力が送信されているものとする。Gtは送信アンテナ20の動作利得を表し、Grは受信アンテナ30の動作利得を表し、r1は送信アンテナ20から生体50までの距離を表し、r2は生体50から受信アンテナ30までの距離を表す。距離r1および距離r2は式11により推定した位置から容易に計算することができる。すると、式12で定義される電力伝達係数は、ρij=Prij/Ptで表されるので、散乱断面積は、
Figure 0006667145
により計算することができる。なお、ここで取り扱う散乱断面積では、生体50全体の散乱断面積ではなく、呼吸、心拍、体動などの生体50のバイタル活動の影響を受けることにより生じた変動成分に対応する散乱面積だけが考慮されている。さらに、全要素を平均して
Figure 0006667145
によって平均散乱断面積を求める。以降、
Figure 0006667145
を単に散乱断面積と呼称する。 Further, the RCS value is obtained from Expression 12, and the posture of the living body 50 is estimated from the position of the living body 50 and the RCS value. Assuming that the frequency range to be extracted is f1 to f2 (f1 <f2), a power transfer coefficient is obtained from a channel component reflected and observed from a living body.
Figure 0006667145
Can be calculated. Where ρij is a matrix
Figure 0006667145
(I, j) -th element. On the other hand, the power that reaches the i-th receiving antenna element 31 from the j-th transmitting antenna element 21 via the living body 50 is
Figure 0006667145
It is expressed as Here, Pt represents transmission power. It is assumed that equal power is transmitted from all transmission antenna elements 21. Gt represents the operating gain of the transmitting antenna 20, Gr represents the operating gain of the receiving antenna 30, r1 represents the distance from the transmitting antenna 20 to the living body 50, and r2 represents the distance from the living body 50 to the receiving antenna 30. The distance r1 and the distance r2 can be easily calculated from the position estimated by Expression 11. Then, since the power transfer coefficient defined by Equation 12 is expressed by ρij = Prij / Pt, the scattering cross section is
Figure 0006667145
Can be calculated by In the scattering cross section handled here, not the scattering cross section of the whole living body 50 but only the scattering area corresponding to the fluctuation component caused by the influence of the vital activity of the living body 50 such as respiration, heartbeat, and body movement. Is considered. Furthermore, on average, all elements
Figure 0006667145
To determine the average scattering cross section. Or later,
Figure 0006667145
Is simply referred to as a scattering cross section.

なお、例えば、生体のバイタル成分と体動成分とをさらに分離して測定する場合、式12の周波数範囲を調整し、その後の処理を行えば良い事は言うまでもない。   Note that, for example, when the vital component and the body motion component of the living body are further separated and measured, it is needless to say that the frequency range of Expression 12 is adjusted and the subsequent processing is performed.

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

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

センサー10では、送信アンテナ20のN個の送信アンテナ素子21が、生体50が存在しうる所定範囲に対してN個の送信アンテナ素子21を用いてN個の送信信号を送信する(S11)。   In the sensor 10, the N transmission antenna elements 21 of the transmission antenna 20 transmit the N transmission signals using the N transmission antenna elements 21 to a predetermined range where 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 received signals including the plurality of reflected signals of the N transmitted signals transmitted by the transmitting antenna 20 and reflected by the living body 50 (S12).

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

回路40は、第1行列における所定周波数範囲に対応する第2行列を抽出することで、生体50の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する第2行列を抽出する(S14)。   The circuit 40 extracts a second matrix corresponding to a predetermined frequency range in the first matrix, thereby obtaining a second matrix corresponding 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の存在する位置を推定する(S15)。   The circuit 40 estimates the position of the living body 50 with respect to the sensor 10 using the second matrix (S15).

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

回路40は、第1距離および第2距離を用いて、生体50に対するRCS値を算出する(S17)。   The circuit 40 calculates an 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 posture of the living body 50 using the calculated RCS value and the information 42 indicating the correspondence between the RCS value and the posture of the living body 50 stored in the memory 41 (S18).

本実施の形態に係るセンサー10によれば、生体50が存在する位置および当該位置における生体の姿勢の推定を短時間かつ高精度に行うことができる。   According to the sensor 10 according to the present embodiment, the position where the living body 50 exists and the posture of the living body at the position can be estimated 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, it is possible to estimate whether the human is alive and in any of the postures of upright, chair, cross, and supine. As a result, the survival of the human can be confirmed effectively. In addition, since the existence of a person can be confirmed without analyzing the image captured by the camera, the existence of the person can be confirmed while protecting the privacy of the person.

(実施の形態2)
図5は実施の形態2におけるセンサーの構成の一例を示すブロック図である。図6は実施の形態2における回路およびメモリの機能的な構成を示すブロック図である。
(Embodiment 2)
FIG. 5 is a block diagram illustrating an example of a configuration of a sensor according to the second embodiment. FIG. 6 is a block diagram showing a functional configuration of a circuit and a memory according to the second embodiment.

実施の形態2におけるセンサー10Aは、実施の形態1におけるセンサー10と比較して、送信アンテナ20Aが有するN個の送信アンテナ素子21A、および、受信アンテナ30Aが有するM個の受信アンテナ素子31Aの配置の仕方が異なる。また、NおよびMは、それぞれ、3以上の自然数である点も異なる。   The sensor 10A according to the second embodiment is different from the sensor 10 according to the first embodiment in the arrangement of N transmission antenna elements 21A included in the transmission antenna 20A and M reception antenna elements 31A included in the reception antenna 30A. In different ways. Also, the difference is that N and M are each a natural number of 3 or more.

送信アンテナ20Aは、N個の送信アンテナ素子21Aを有する。送信アンテナ20Aは、水平方向(x方向)にN個が並び、かつ、鉛直方向(z方向)にN個が並ぶように矩形配置された、N個(N=N×N)の送信アンテナ素子21Aで構成されるアレーアンテナを有する。つまり、N個の送信アンテナ素子21Aのうち少なくとも3個の送信アンテナ素子21Aは、鉛直方向および水平方向の異なる位置に配置されている。 The transmission antenna 20A has N transmission antenna elements 21A. Transmitting antenna 20A is 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 ) Has an array antenna composed of the transmission antenna elements 21A. That is, at least three transmission antenna elements 21A out of the N transmission antenna elements 21A are arranged at different positions in the vertical and horizontal directions.

受信アンテナ30Aは、M個の受信アンテナ素子31Aを有する。受信アンテナ30Aは、水平方向(x方向)にM個が並び、かつ、鉛直方向(z方向)にM個が並ぶように矩形配置された、M個(M=M×M)の受信アンテナ素子31Aで構成されるアレーアンテナを有する。つまり、M個の受信アンテナ素子31Aのうち少なくとも3個の受信アンテナ素子31Aは、鉛直方向および水平方向の異なる位置に配置されている。 The receiving antenna 30A has M receiving antenna elements 31A. Receiving antenna 30A is lined M x pieces in the horizontal direction (x-direction), and, M z pieces in the vertical direction (z-direction) is rectangular arranged side by side, M pieces (M = M x × M z ) Of the receiving antenna element 31A. That is, at least three of the M receiving antenna elements 31A are arranged at different positions in the vertical and horizontal directions.

ここで、送信アンテナ20Aに対して任意に設定された水平面上の方向である第1基準方向と、送信アンテナ20Aから生体50Aへの方向である第1生体方向とのなす角をφとする。また、鉛直方向と、第1生体方向との為す角である生体50Aの仰角をθとする。また、受信アンテナ30Aに対して任意に設定された水平面上の方向である第2基準方向と、受信アンテナ30Aから生体50Aへの方向である第2生体方向との為す角である生体50Aの仰角をφとする。また、鉛直方向と、第2生体方向とのなす角をθとする。生体50Aがバイタル活動を行っている部位の中心座標を(x,y,z)とすると、送信アンテナ20A、受信アンテナ30Aおよび生体50Aの位置関係によって、方向(θ,θ,φ,φ)と座標(x,y,z)は相互に変換可能である。 Here, the optionally the first reference direction is set direction on a horizontal plane, the angle phi T between the first biological direction is the direction from the transmitting antenna 20A to the living body 50A to the transmitting antenna 20A . Further, the vertical direction, the elevation of the living body 50A is to the corner for the first biological direction and theta T. In addition, the elevation angle of the living body 50A, which is an angle between a second reference direction that is a direction on the horizontal plane arbitrarily set with respect to the receiving antenna 30A and a second living body direction that is a direction from the receiving antenna 30A to the living body 50A. It is referred to as φ R. Further, the vertical direction, the angle between the second living direction and theta R. The center coordinates of the site where biological 50A is performing a vital activities (x b, y b, z b) When, by the positional relationship of the transmitting antenna 20A, the reception antenna 30A and biological 50A, the direction (θ T, θ R, φ T , φ R ) and coordinates (x b , y b , z b ) are mutually convertible.

また、センサー10Aは、センサー10と比較して、回路40Aが行う処理が異なる。生体50Aは、例えば、高さzにて周囲と比べて最も強くバイタル活動を行っている生体である。生体50Aは、例えば、呼吸による体表面変位が最も大きい腹部を有する。本実施の形態におけるセンサー10Aでは、実施の形態1におけるセンサー10で推定する水平面における生体50の位置に加えて、生体50Aの例えば腹部の鉛直方向における位置である鉛直位置zb(z軸方向の位置)である高さを含む3次元位置を推定する。つまり、本実施の形態のセンサー10Aの回路40Aは、実施の形態1の回路40の位置推定処理部430の代わりに、上記3次元位置を推定する処理を行う、3次元位置推定処理部430Aを有する。 The sensor 10A differs from the sensor 10 in the processing performed by the circuit 40A. Biological 50A is, for example, a living organism is carried out most strongly vital activities compared to the ambient at a height z b. The living body 50A has, for example, an abdomen where displacement of the body surface due to respiration is greatest. In the sensor 10A according to the present embodiment, in addition to the position of the living body 50 in the horizontal plane estimated by the sensor 10 in the first embodiment, a vertical position zb (a position in the z-axis direction) which is, for example, a vertical position of the abdomen of the living body 50A. ) Is estimated. That is, the circuit 40A of the sensor 10A according to the present embodiment includes a three-dimensional position estimation processing unit 430A that performs the process of estimating the three-dimensional position, instead of the position estimation processing unit 430 of the circuit 40 according to the first embodiment. Have.

また、センサー10Aは、センサー10と比較して、メモリ41に記憶されている対応関係を示す情報42Aが異なる。本実施の形態において、メモリ41に記憶されているRCS値および生体50の姿勢の対応関係を示す情報42Aとは、図7に示すように、仰臥、胡坐、椅座および直立で示される各姿勢に予め対応付けられたRCS値の範囲および高さの範囲を示す情報である。例えば、仰臥は、第1RCS範囲および第1高さ範囲と対応付けられており、胡坐は、第2RCS範囲および第2高さ範囲と対応付けられており、椅座は、第3RCS範囲および第3高さ範囲と対応付けられており、直立は、第4RCS範囲および第4高さ範囲と対応付けられている。なお、第1RCS範囲〜第4RCS範囲は、それぞれ異なるRCS値の範囲である。また、第1高さ範囲〜第4高さ範囲は、それぞれことなる高さの範囲である。   The sensor 10A differs from the sensor 10 in the information 42A indicating the correspondence stored in the memory 41. In the present embodiment, the information 42A indicating the correspondence between the RCS value and the posture of the living body 50 stored in the memory 41 is, as shown in FIG. Is information indicating the range of the RCS value and the range of the height which are associated in advance with. For example, the supine is associated with the first RCS range and the first height range, the cross is associated with the second RCS range and the second height range, and the chair is associated with the third RCS range and the third height range. The upright is associated with the fourth RCS range and the fourth height range. The first RCS range to the fourth RCS range are ranges of different RCS values. Further, the first to fourth height ranges are different height ranges.

センサー10Aのその他の構成は、センサー10と同じであるため、当該構成については、実施の形態1と同一の符号を付してその説明を省略する。   The other configuration of the sensor 10A is the same as that of the sensor 10, and therefore, the configuration is denoted by the same reference numeral as in the first embodiment, and the description is omitted.

本実施の形態のセンサー10Aは、実施の形態1におけるセンサー10で推定する、水平面における生体50の位置に加えて、鉛直方向における生体50の位置である鉛直位置を推定する。   The sensor 10A according to the present embodiment estimates a vertical position, which is the position of the living body 50 in the vertical direction, in addition to the position of the living body 50 in the horizontal plane, which is estimated by the sensor 10 according to the first embodiment.

次に実施の形態2のセンサー10Aの動作原理の詳細を、数式を用いて説明する。式1から式4までは、実施の形態1と同様な処理を行い、式4によって生体成分から相関行列

Figure 0006667145
を得る。Rは実施の形態1と同様に0〔Hz〕<f<3〔Hz〕の間で平均化されている。 Next, details of the operation principle of the sensor 10A according to the second embodiment will be described using mathematical expressions. Equations 1 to 4 perform the same processing as that of the first embodiment.
Figure 0006667145
Get. R is averaged in the range of 0 [Hz] <f <3 [Hz] as in the first embodiment.

次に、式5と同様に算出した相関行列Rを固有値分解することで

Figure 0006667145
および
Figure 0006667145
を計算する。以降では、MUSIC法を用いた高さ方向を含む位置推定方法について説明する。送信アンテナ20Aから生体50Aへ向かう(θT、φT)方向を示すステアリングベクトルと、受信アンテナ30Aから生体50Aへ向かう(θR、φR)方向を示すステアリングベクトルとは、それぞれ、
Figure 0006667145
Figure 0006667145
と表される。ここで、
Figure 0006667145
Figure 0006667145
と表される。ここで、kは波数を表し、dTxおよびdTzはそれぞれ送信アンテナ素子21Aのx方向およびz方向における素子間隔を表し、dRxおよびdRzはそれぞれ受信アンテナ素子31Aのx方向およびz方向における素子間隔を表す。なお、本実施の形態では素子間隔が同一の方向においては一定のリニアアレーアンテナを想定している。 Next, by eigenvalue decomposition of the correlation matrix R calculated in the same manner as in Equation 5,
Figure 0006667145
and
Figure 0006667145
Is calculated. Hereinafter, a position estimation method including the height direction using the MUSIC method will be described. A steering vector indicating a (θT, φT) direction from the transmitting antenna 20A toward the living body 50A and a steering vector indicating a (θR, φR) direction from the receiving antenna 30A toward the living body 50A,
Figure 0006667145
Figure 0006667145
It is expressed as here,
Figure 0006667145
Figure 0006667145
It is expressed as Here, k represents a wave number, d Tx and d Tz represent element intervals in the x direction and z direction of the transmitting antenna element 21A, respectively, and d Rx and d Rz represent elements in the x direction and z direction of the receiving antenna element 31A, respectively. Indicates the element spacing. In the present embodiment, it is assumed that the linear array antenna is constant in the direction in which the element spacing is the same.

Txは、例えば、複数の送信アンテナ素子21Aのうちx方向において互いに隣接する2つの送信アンテナ素子21Aの間隔を表す。dTzは、例えば、複数の送信アンテナ素子21Aのうちz方向において互いに隣接する2つの送信アンテナ素子21Aの間隔を表す。また、dRxは、例えば、複数の受信アンテナ素子31Aのうちx方向において互いに隣接する2つの受信アンテナ素子31Aの間隔を表す。また、dRzは、例えば、複数の受信アンテナ素子31Aのうちz方向において互いに隣接する2つの受信アンテナ素子31Aの間隔を表す。これらのステアリングベクトルのクロネッカ積を求めると、

Figure 0006667145
となる。ステアリングベクトルa(θ,φ,θ,φ)は、MN×1の要素を持つベクトルであり、出発角θT、φTと到来角θR、φRとの4変数を持つ関数となる。以降、a(θ,φ,θ,φ)をステアリングベクトルと定義する。検出範囲内に存在する生体の数をLとすると、評価関数
Figure 0006667145
によって生体位置を特定する。実施の形態1と同様に、式22のMUSICスペクトラムの極大点を探索することで、鉛直位置を含む、送信アンテナ20Aおよび受信アンテナ30Aから見た生体50Aの3次元位置を特定できる。 d Tx represents, for example, an interval between two transmission antenna elements 21A that are adjacent to each other in the x direction among the plurality of transmission antenna elements 21A. d Tz represents, for example, an interval between two transmission antenna elements 21A adjacent to each other in the z direction among the plurality of transmission antenna elements 21A. D Rx represents, for example, an interval between two reception antenna elements 31A that are adjacent to each other in the x direction among the plurality of reception antenna elements 31A. Further, d Rz represents, for example, an interval between two reception antenna elements 31A adjacent to each other in the z direction among the plurality of reception antenna elements 31A. When calculating the Kronecker product of these steering vectors,
Figure 0006667145
Becomes The steering vector a (θ T , φ T , θ R , φ R ) is a vector having MN × 1 elements, and is a function having four variables of departure angles θT, φT and arrival angles θR, φR. Hereinafter, a (θ T , φ T , θ R , φ R ) is defined as a steering vector. Assuming that the number of living organisms present in the detection range is L, the evaluation function
Figure 0006667145
The body position is specified by Similar to Embodiment 1, by searching for the maximum point of the MUSIC spectrum of Expression 22, the three-dimensional position of the living body 50A, including the vertical position, viewed from the transmitting antenna 20A and the receiving antenna 30A can be specified.

なお、本実施の形態におけるセンサー10Aの動作の説明は、図4で説明したフローチャートのステップS15の生体の位置の推定において、上記の3次元位置の推定を行うことで説明できるため省略する。   The description of the operation of the sensor 10A in the present embodiment will be omitted because the estimation of the three-dimensional position described above can be made in the estimation of the position of the living body in step S15 of the flowchart described with reference to FIG.

図8は、実施の形態2のセンサーの効果を確認するために実施した実験の概要を表す図である。   FIG. 8 is a diagram illustrating an outline of an experiment performed to confirm the effect of the sensor according to the second embodiment.

図8に示すように、送信アンテナ20Aは4×4=16の送信アンテナ素子21Aが配置されることで構成される方形アレーアンテナである。送信アンテナ素子21Aはパッチアンテナとする。受信アンテナ30Aは4×4=16の受信アンテナ素子31Aが配置されることで構成される方形アレーアンテナである。受信アンテナ素子31Aはパッチアンテナとする。送信アンテナ素子21Aおよび受信アンテナ素子31Aは、ともに、素子間隔が0.5波長で配置されている。   As shown in FIG. 8, the transmission antenna 20A is a rectangular array antenna configured by arranging 4 × 4 = 16 transmission antenna elements 21A. The transmitting antenna element 21A is a patch antenna. The receiving antenna 30A is a rectangular array antenna configured by arranging 4 × 4 = 16 receiving antenna elements 31A. The receiving antenna element 31A is a patch antenna. The transmitting antenna element 21A and the receiving antenna element 31A are both arranged at an element interval of 0.5 wavelength.

被検者は1名であり、直立、椅座(椅子に座った状態)、胡座(あぐら)、仰臥(仰向け)の4状態にて測定を行った。また、被験者の位置は(X,Y)=(5.0,1.0)mである。仰臥以外の姿勢では、被験者はアンテナ方向(Y方向)を向いている。仰臥では、被験者はアンテナ方向(Y方向)に足を向けている。実験1回につき約5秒間観測を行った。   The number of subjects was one, and the measurement was performed in four states: erect, chair (sitting on a chair), Goza (Agura), and supine (facing). The position of the subject is (X, Y) = (5.0, 1.0) m. In a posture other than the supine position, the subject faces the antenna direction (Y direction). In the supine position, the subject turns his or her feet in the antenna direction (Y direction). Observation was performed for about 5 seconds for each experiment.

図9は、図8で示した実験系を用いた実験結果を示す図である。図9において、横軸は推定された高さzbを示し、縦軸は算出した散乱断面積(RCS:Radar Cross Section)を示す。   FIG. 9 is a diagram showing an experimental result using the experimental system shown in FIG. In FIG. 9, the horizontal axis indicates the estimated height zb, and the vertical axis indicates the calculated scattering cross section (RCS: Radar Cross Section).

図9に示すように、姿勢によって高さおよびRCS値が分布している領域が異なることが分かる。つまり、推定された高さおよび算出したRCS値から被験者の姿勢を推定できることが確認できる。例えば、図10に示すように、対応関係を示す情報42Aを設定しておくことで、推定した3次元位置と、算出されたRCS値と、メモリ41に記憶されている、対応関係を示す情報42Aと、を用いて、生体50Aが、直立、椅座、胡坐、および、仰臥のいずれの姿勢であるかを推定できる。なお、図10は実験結果から得られた第1RCS範囲〜第4RCS範囲の具体例、および、第1高さ範囲〜第4高さ範囲の具体例を示す図である。なお、図10の第1RCS範囲〜第4RCS範囲は、実施の形態1の対応関係を示す情報42の第1RCS範囲〜第4RCS範囲にも適用可能である。   As shown in FIG. 9, it can be seen that the region where the height and the RCS value are distributed differs depending on the posture. That is, it can be confirmed that the posture of the subject can be estimated from the estimated height and the calculated RCS value. For example, as shown in FIG. 10, by setting information 42A indicating the correspondence, the estimated three-dimensional position, the calculated RCS value, and the information indicating the correspondence stored in the memory 41 are stored. Using 42A, it can be estimated whether the living body 50A is in an upright posture, a chair, a cross, or a supine position. FIG. 10 is a diagram showing a specific example of the first to fourth RCS ranges and a specific example of the first to fourth height ranges obtained from the experimental results. Note that the first to fourth RCS ranges in FIG. 10 are also applicable to the first to fourth RCS ranges in the information 42 indicating the correspondence in the first embodiment.

本実施の形態に係るセンサー10Aによれば、生体50が存在する3次元位置および当該位置における生体の姿勢の推定を短時間かつ高精度に行うことができる。   According to the sensor 10A according to the present embodiment, the three-dimensional position where the living body 50 exists and the posture of the living body at the position can be estimated in a short time and with high accuracy.

なお、上記各実施の形態において、各構成要素は、専用のハードウェアで構成されるか、各構成要素に適したソフトウェアプログラムを実行することによって実現されてもよい。各構成要素は、CPUまたはプロセッサなどのプログラム実行部が、ハードディスクまたは半導体メモリなどの記録媒体に記録されたソフトウェアプログラムを読み出して実行することによって実現されてもよい。ここで、上記各実施の形態のセンサー10、10Aなどを実現するソフトウェアは、次のようなプログラムである。   In each of the above embodiments, each component may be configured by dedicated hardware, or may be realized by executing a software program suitable for each component. Each component may be realized by a program execution unit such as a CPU or a processor reading and executing a software program recorded on a recording medium such as a hard disk or a semiconductor memory. Here, the software for realizing the sensors 10, 10A and the like in each of the above embodiments is the following program.

すなわち、このプログラムは、コンピュータに、N個(Nは2以上の自然数)の送信アンテナ素子を有する送信アンテナと、M個(Mは2以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、メモリとを備えるセンサーによる推定方法であって、生体が存在しうる所定範囲に対してN個の送信アンテナ素子を用いてN個の送信信号を送信し、送信した前記N個の送信信号の一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号を複数の受信アンテナ素子のそれぞれを用いて受信し、前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、前記第2行列を用いて、前記センサーに対する前記生体の存在する位置を推定し、前記推定した位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、前記第1距離および前記第2距離を用いて、前記生体に対するRCS(Rader cross-section)値を算出し、前記算出されたRCS値と、前記メモリに記憶されている、RCS値および前記生体の姿勢の対応関係を示す情報と、を用いて、前記生体の姿勢を推定する、推定方法を実行させる。   That is, the program includes a transmitting antenna having N (N is a natural number of 2 or more) transmitting antennas, a receiving antenna having M (M is a natural number of 2 or more) receiving antenna elements, a circuit, And a memory, wherein the N transmission signals are transmitted by using N transmission antenna elements to a predetermined range where a living body can exist, and the transmitted N transmission signals are transmitted. A part of the transmission signals is received N reception signals including a reflection signal reflected by the living body using each of the plurality of reception antenna elements, and is received by each of the M reception antenna elements for a predetermined period. From each of the N received signals obtained, each complex representing a propagation characteristic between each of the N transmit antenna elements and each of the M receive antenna elements. By calculating a first matrix of N × M having a transfer function as a component and extracting a second matrix corresponding to a predetermined frequency range in the first matrix, at least one of respiration, heartbeat, and body motion of the living body Extracting the second matrix corresponding to the component affected by the vital activity including the, using the second matrix to estimate the position of the living body with respect to the sensor, the estimated position, Based on a position of a transmitting antenna and a position of the receiving antenna, a first distance indicating a distance between the living body and the transmitting antenna and a second distance indicating a distance between the living body and the receiving antenna are calculated. Then, an RCS (Rader cross-section) value for the living body is calculated using the first distance and the second distance, and the calculated RCS value and the RCS value stored in the memory are calculated. With, and information indicating the correspondence between the orientation of the body, to estimate the posture of the living body to execute an estimation process.

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

上記実施の形態では、センサー10、10Aは、生体が直立、椅座、胡坐、および、仰臥のいずれの姿勢であるかを推定するとしたが、これに限らない。例えば、生体が、送信アンテナ20、20Aおよび受信アンテナ30、30Aの並び方向に垂直な方向に対して正対している姿勢であるか否かを推定してもよい。具体的には、図1または図5におけるy方向に対して正対している姿勢であるか、y方向に対して横向きの姿勢であるかを推定してもよい。例えば、ヒトがy方向に対して正対している場合、y方向に対して正対していない場合と比較してRCS値が小さくなる。なお、y方向に対して横向きの姿勢とは、x方向に対して正対している姿勢である。   In the above embodiment, the sensors 10 and 10A estimate whether the living body is in an upright posture, a chair, a cross, or a supine position. However, the present invention is not limited to this. For example, it may be estimated whether or not the living body is in a posture facing the direction perpendicular to the direction in which the transmitting antennas 20 and 20A and the receiving antennas 30 and 30A are arranged. Specifically, it may be estimated whether the posture is directly facing the y direction in FIG. 1 or FIG. 5 or the posture is lateral to the y direction. For example, when a person faces directly in the y direction, the RCS value becomes smaller than when the person does not face the y direction. In addition, the attitude that is lateral to the y direction is an attitude that faces directly to the x direction.

上記実施の形態1では、N個の送信アンテナ素子21は、水平面上の第1所定方向に並んで配置され、かつ、M個の受信アンテナ素子31は、水平面上の第2所定方向に並んで配置されるとしたが、これに限らない。つまり、第1所定方向および第2所定方向は、水平面上の方向に限らずに、例えば、鉛直方向を含む平面上の方向であってもよいし、水平面から傾いた平面上の方向であってもよい。この場合、回路40は、第1所定方向および第2所定方向を含む平面上における生体の位置を推定することとなる。   In the first embodiment, N transmitting antenna elements 21 are arranged in a first predetermined direction on a horizontal plane, and M receiving antenna elements 31 are arranged in a second predetermined direction on a horizontal plane. Although it is arranged, it is not limited to this. That is, the first predetermined direction and the second predetermined direction are not limited to the directions on the horizontal plane, and may be, for example, directions on a plane including a vertical direction, or directions on a plane inclined from the horizontal plane. Is also good. In this case, the circuit 40 estimates the position of the living body on a plane including the first predetermined direction and the second predetermined direction.

本発明は、無線信号を利用した動体の方向や位置を推定するセンサーおよび推定方法に利用できる。特に、生体と機械を含む動体の方向や位置を測定する測定器、動体の方向や位置に応じた制御を行う家電機器、動体の侵入を検知する監視装置などに搭載される方向推定方法、方向推定方法,位置推定方法および散乱断面積を用いた高々推定装置に利用できる。   INDUSTRIAL APPLICATION This invention can be utilized for the sensor and estimation method which estimate the direction and position of a moving body using a radio signal. In particular, direction estimation methods and directions 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, monitoring devices that detect intrusion of moving objects, etc. The present invention can be used for an estimation method, a position estimation method, and an estimating apparatus at most using a scattering cross section.

10、10A センサー
20、20A 送信アンテナ
21、21A 送信アンテナ素子
30、30A 受信アンテナ
31、31A 受信アンテナ素子
40、40A 回路
41 メモリ
42、42A 対応関係を示す情報
50、50A 生体
410 複素伝達関数算出部
420 生体成分算出部
430 位置推定処理部
430A 3次元位置推定処理部
440 RCS算出部
450 姿勢推定部
10, 10A sensor 20, 20A transmitting antenna 21, 21A transmitting antenna element 30, 30A receiving antenna 31, 31A receiving antenna element 40, 40A circuit 41 memory 42, 42A information indicating correspondence 50, 50A living body 410 complex transfer function calculating unit 420 biological component calculation section 430 position estimation processing section 430A three-dimensional position estimation processing section 440 RCS calculation section 450 attitude estimation section

Claims (5)

センサーであって、
生体が存在しうる所定範囲に対して送信信号をそれぞれが送信するN個(Nは2以上の自然数)の送信アンテナ素子を有する送信アンテナと、
前記N個の送信アンテナ素子により送信された前記N個の送信信号のうちの一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号をそれぞれが受信するM個(Mは2以上の自然数)の受信アンテナ素子を有する受信アンテナと、
回路と、
メモリと、を備え、
前記回路は、
前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、
前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、
前記第2行列を用いて、前記センサーに対する前記生体の存在する位置を推定し、
前記推定した位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、
前記第1距離および前記第2距離を用いて、前記生体に対するRCS(Rader cross-section)値を算出し、
前記算出されたRCS値と、前記メモリに記憶されている、RCS値および前記生体の姿勢の対応関係を示す情報と、を用いて、前記生体の姿勢を推定する、
センサー。
A sensor,
A transmission antenna having N (N is a natural number of 2 or more) transmission antenna elements each transmitting a transmission signal to a predetermined range where a living body can exist;
Some of the N transmission signals transmitted by the N transmission antenna elements are M (M) signals each of which receives N reception signals including a reflection signal reflected by the living body. Is a natural number of 2 or more) receiving antenna elements,
Circuit and
And a memory,
The circuit is
Propagation characteristics between each of the N transmission antenna elements and each of the M reception antenna elements from each of the N reception signals received in each of the M reception antenna elements for a predetermined period. A first matrix of N × M having each complex transfer function as a component is calculated,
By extracting a second matrix corresponding to a predetermined frequency range in the first matrix, the second matrix corresponding to a component affected by vital activity including at least one of respiration, heartbeat, and body movement of the living body Extract
Using the second matrix, estimate the position of the living body with respect to the sensor,
A first distance indicating a distance between the living body and the transmitting antenna based on the estimated position, a position of the transmitting antenna, and a position of the receiving antenna, and a distance between the living body and the receiving antenna Calculate the second distance indicating
Using the first distance and the second distance, calculate an RCS (Rader cross-section) value for the living body,
Estimating the posture of the living body using the calculated RCS value and information indicating the correspondence between the RCS value and the posture of the living body stored in the memory,
sensor.
前記所定期間は、前記生体の呼吸、心拍、および、体動の少なくとも1つの周期の略半分である、
請求項1に記載のセンサー。
The predetermined time period is substantially half of at least one cycle of respiration, heartbeat, and body movement of the living body.
The sensor according to claim 1.
前記回路は、
前記生体が、前記送信アンテナおよび前記受信アンテナの並び方向に垂直な方向に対して正対している姿勢であるか否かを推定する、
請求項1または2に記載のセンサー。
The circuit is
Estimating whether the living body is in a position facing directly with respect to a direction perpendicular to the arrangement direction of the transmitting antenna and the receiving antenna,
The sensor according to claim 1.
前記Nは3以上の自然数であり、
前記N個の送信アンテナ素子のうち少なくとも3個の送信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記Mは3以上の自然数であり、
前記M個の受信アンテナ素子のうち少なくとも3個の受信アンテナ素子は、それぞれ、鉛直方向および水平方向の異なる位置に配置され、
前記対応関係を示す情報は、前記センサーに対する前記生体の存在する鉛直方向における位置である鉛直位置、RCS値、および、前記生体の姿勢の対応関係を示し、
前記対応関係を示す情報において対応付けられている前記生体の姿勢は、直立、椅座、胡坐、および、仰臥を含み、
前記回路は、
前記第2行列を用いて、前記鉛直位置を含む3次元位置を推定し、
前記推定した前記3次元位置と、前記算出されたRCS値と、前記メモリに記憶されている、前記対応関係を示す情報と、を用いて、前記生体が、直立、椅座、胡坐、および、仰臥のいずれの姿勢であるかを推定する、
請求項1から3のいずれか1項に記載のセンサー。
N is a natural number of 3 or more;
At least three transmission antenna elements of the N transmission antenna elements are respectively arranged at different positions in a vertical direction and a horizontal direction,
M is a natural number of 3 or more;
At least three of the M receiving antenna elements are arranged at different positions in the vertical and horizontal directions, respectively.
The information indicating the correspondence relationship, a vertical position that is a position in the vertical direction in which the living body exists with respect to the sensor, an RCS value, and indicates the correspondence relationship between the posture of the living body,
The posture of the living body that is associated in the information indicating the correspondence relationship includes upright, chair, cross, and supine,
The circuit is
Using the second matrix, estimate a three-dimensional position including the vertical position,
Using the estimated three-dimensional position, the calculated RCS value, and the information indicating the correspondence stored in the memory, the living body is upright, a chair, a cross, and Estimate which posture is supine,
The sensor according to any one of claims 1 to 3.
N個(Nは2以上の自然数)の送信アンテナ素子を有する送信アンテナと、M個(Mは2以上の自然数)の受信アンテナ素子を有する受信アンテナと、回路と、メモリとを備えるセンサーによる推定方法であって、
生体が存在しうる所定範囲に対してN個の送信アンテナ素子を用いてN個の送信信号を送信し、
送信した前記N個の送信信号の一部の送信信号が前記生体により反射された反射信号を含むN個の受信信号を複数の受信アンテナ素子のそれぞれを用いて受信し、
前記M個の受信アンテナ素子のそれぞれにおいて所定期間で受信された前記N個の受信信号のそれぞれから、前記N個の送信アンテナ素子それぞれと、前記M個の受信アンテナ素子それぞれとの間の伝搬特性を示す各複素伝達関数を成分とする、N×Mの第1行列を算出し、
前記第1行列における所定周波数範囲に対応する第2行列を抽出することで、前記生体の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する前記第2行列を抽出し、
前記第2行列を用いて、前記センサーに対する前記生体の存在する位置を推定し、
前記推定した位置と、前記送信アンテナの位置と、前記受信アンテナの位置と、に基づいて、前記生体と前記送信アンテナとの距離を示す第1距離、および、前記生体と前記受信アンテナとの距離を示す第2距離を算出し、
前記第1距離および前記第2距離を用いて、前記生体に対するRCS(Rader cross-section)値を算出し、
前記算出されたRCS値と、前記メモリに記憶されている、RCS値および前記生体の姿勢の対応関係を示す情報と、を用いて、前記生体の姿勢を推定する、
推定方法。
Estimation by a sensor including N (N is a natural number of 2 or more) transmitting antenna elements, M receiving antennas (M is a natural number of 2 or more) receiving antenna elements, a circuit, and a memory The method
Transmit N transmission signals using N transmission antenna elements for a predetermined range where a living body can exist,
A part of the transmitted N transmission signals, N reception signals including a reflection signal reflected by the living body are received using each of the plurality of reception antenna elements,
Propagation characteristics between each of the N transmission antenna elements and each of the M reception antenna elements from each of the N reception signals received in each of the M reception antenna elements for a predetermined period. A first matrix of N × M having each complex transfer function as a component is calculated,
By extracting a second matrix corresponding to a predetermined frequency range in the first matrix, the second matrix corresponding to a component affected by vital activity including at least one of respiration, heartbeat, and body movement of the living body Extract
Using the second matrix, estimate the position of the living body with respect to the sensor,
A first distance indicating a distance between the living body and the transmitting antenna based on the estimated position, a position of the transmitting antenna, and a position of the receiving antenna, and a distance between the living body and the receiving antenna Calculate the second distance indicating
Using the first distance and the second distance, calculate an RCS (Rader cross-section) value for the living body,
Estimating the posture of the living body using the calculated RCS value and information indicating the correspondence between the RCS value and the posture of the living body stored in the memory,
Estimation method.
JP2017063507A 2016-07-05 2017-03-28 Sensor and estimation method Active JP6667145B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201710416994.9A CN107576957B (en) 2016-07-05 2017-06-06 Sensor and Estimation Method
US15/631,367 US10473756B2 (en) 2016-07-05 2017-06-23 Sensor and estimating method
US16/598,628 US10895631B2 (en) 2016-07-05 2019-10-10 Sensor and estimating method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016133283 2016-07-05
JP2016133283 2016-07-05

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2020019086A Division JP6839807B2 (en) 2016-07-05 2020-02-06 Estimator and program

Publications (2)

Publication Number Publication Date
JP2018008021A JP2018008021A (en) 2018-01-18
JP6667145B2 true JP6667145B2 (en) 2020-03-18

Family

ID=60994506

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2017063507A Active JP6667145B2 (en) 2016-07-05 2017-03-28 Sensor and estimation method
JP2020019086A Active JP6839807B2 (en) 2016-07-05 2020-02-06 Estimator and program

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2020019086A Active JP6839807B2 (en) 2016-07-05 2020-02-06 Estimator and program

Country Status (2)

Country Link
JP (2) JP6667145B2 (en)
CN (1) CN107576957B (en)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10729339B2 (en) * 2018-02-22 2020-08-04 Vayyar Imaging Ltd. Detecting and measuring correlated movement with MIMO radar
KR102427925B1 (en) 2018-03-05 2022-08-02 삼성전자주식회사 Electronic device and method for operating a predetermined function corresponding to a distance of object based on a signal output through an antenna
US11163057B2 (en) 2018-05-02 2021-11-02 Panasonic Intellectual Property Management Co., Ltd. Estimation device, living body count estimation device, estimation method, and recording medium
US10772511B2 (en) * 2018-05-16 2020-09-15 Qualcomm Incorporated Motion sensor using cross coupling
US11226408B2 (en) * 2018-07-03 2022-01-18 Panasonic Intellectual Property Management Co., Ltd. Sensor, estimating device, estimating method, and recording medium
JP6587199B1 (en) * 2018-07-03 2019-10-09 パナソニックIpマネジメント株式会社 Estimation apparatus and estimation method
JP7174921B2 (en) * 2018-07-03 2022-11-18 パナソニックIpマネジメント株式会社 Sensor, estimation device, estimation method, and program
IT201800008235A1 (en) * 2018-08-29 2020-02-29 De' Longhi Appliances Srl METHOD OF ACTIVATION AND CONTROL OF A CONDITIONING APPARATUS
JP7281784B2 (en) * 2018-12-28 2023-05-26 パナソニックIpマネジメント株式会社 Estimation method and estimation device
CN111381227B (en) * 2018-12-28 2024-04-26 松下知识产权经营株式会社 Inference method and inference device
CN111381228B (en) 2018-12-28 2024-10-01 松下知识产权经营株式会社 Inference method and inference device
US11378685B2 (en) * 2019-02-27 2022-07-05 Leolabs, Inc. Systems, devices, and methods for determining space object attitude stabilities from radar cross-section statistics
CN113260871B (en) * 2019-07-02 2025-06-17 松下知识产权经营株式会社 sensor
JP7617558B2 (en) * 2021-02-26 2025-01-20 パナソニックIpマネジメント株式会社 Discrimination system, discrimination method, and program
CN118414558A (en) 2021-12-27 2024-07-30 松下知识产权经营株式会社 Estimation device, estimation method and program
CN120360524B (en) * 2025-06-25 2025-10-03 北京理工雷科电子信息技术有限公司 Method, equipment and medium for detecting falling bed facing household health monitoring radar

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3750904B2 (en) * 1999-10-19 2006-03-01 三菱電機株式会社 Radar equipment
JP3783843B2 (en) * 2001-03-14 2006-06-07 三菱電機株式会社 Radar equipment
JP4016826B2 (en) * 2002-12-10 2007-12-05 株式会社デンソー Object labeling method and apparatus, program
JP2006081771A (en) * 2004-09-17 2006-03-30 Canon Inc Method and apparatus for estimating information about living organisms using electromagnetic waves
JP2007071605A (en) * 2005-09-05 2007-03-22 Optex Co Ltd Crime prevention sensor
US7567200B1 (en) * 2006-04-27 2009-07-28 Josef Osterweil Method and apparatus for body position monitor and fall detect ion using radar
US7864103B2 (en) * 2007-04-27 2011-01-04 Accipiter Radar Technologies, Inc. Device and method for 3D height-finding avian radar
BR112012024337A2 (en) * 2010-03-26 2017-10-03 Vayyar Imaging Ltd METHOD AND APPARATUS FOR ENHANCEMENT OF MICROWAVE IMAGE FORMATION OF AN OBJECT, MIMO RADAR SYSTEM, AND COMPUTER PROGRAM PRODUCT
US8970425B2 (en) * 2011-06-09 2015-03-03 Sony Corporation Radar apparatus and method
JP6273546B2 (en) * 2013-10-02 2018-02-07 三星電子株式会社Samsung Electronics Co.,Ltd. Position estimation apparatus and position estimation method
EP3055708A2 (en) * 2013-10-09 2016-08-17 Massachusetts Institute Of Technology Motion tracking via body radio reflections
JP6402398B2 (en) * 2013-12-17 2018-10-10 三星電子株式会社Samsung Electronics Co.,Ltd. Processing apparatus and processing method
JP6256681B2 (en) * 2013-12-17 2018-01-10 三星電子株式会社Samsung Electronics Co.,Ltd. Calibration apparatus, position estimation apparatus, calibration method, and position estimation method
JP5848469B1 (en) * 2015-01-23 2016-01-27 株式会社光波 Biological condition detection device

Also Published As

Publication number Publication date
JP2020075168A (en) 2020-05-21
JP6839807B2 (en) 2021-03-10
CN107576957A (en) 2018-01-12
CN107576957B (en) 2023-01-06
JP2018008021A (en) 2018-01-18

Similar Documents

Publication Publication Date Title
JP6667145B2 (en) Sensor and estimation method
US10895631B2 (en) Sensor and estimating method
JP6865394B2 (en) program
CN108279413B (en) Sensor and biological action estimation method
US11150336B2 (en) Positioning sensor and direction estimation method
CN110673135B (en) Sensor, estimation device, estimation method, and program recording medium
JP6868846B2 (en) Sensors and methods
EP3289971B1 (en) Biometric device and biometric method
JP2020109389A (en) Estimating method, estimating device, and program
JP2014228291A (en) Radio detection device and radio detection method
JP7174921B2 (en) Sensor, estimation device, estimation method, and program
US11047968B2 (en) Estimating method and estimating device
CN111381229B (en) Estimation method, estimation device, and recording medium
US12303236B2 (en) Contactless seismocardiography
CN111381228B (en) Inference method and inference device
JP2019197039A (en) Estimation device, living body number estimation device, estimation method, and program
JP7429877B2 (en) Estimation method and device
KR101735956B1 (en) Apparatus and method for estimating direction of vessel
CN116321425A (en) Method and device for indoor multi-target human trajectory tracking

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191120

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20191120

TRDD Decision of grant or rejection written
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20191202

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200114

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200206

R151 Written notification of patent or utility model registration

Ref document number: 6667145

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151