JP7451040B2 - Biological state estimation device - Google Patents

Biological state estimation device Download PDF

Info

Publication number
JP7451040B2
JP7451040B2 JP2020083450A JP2020083450A JP7451040B2 JP 7451040 B2 JP7451040 B2 JP 7451040B2 JP 2020083450 A JP2020083450 A JP 2020083450A JP 2020083450 A JP2020083450 A JP 2020083450A JP 7451040 B2 JP7451040 B2 JP 7451040B2
Authority
JP
Japan
Prior art keywords
living body
state
signal
task
unit
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
JP2020083450A
Other languages
Japanese (ja)
Other versions
JP2021177830A (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.)
Japan Radio Co Ltd
Original Assignee
Japan Radio 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 Japan Radio Co Ltd filed Critical Japan Radio Co Ltd
Priority to JP2020083450A priority Critical patent/JP7451040B2/en
Publication of JP2021177830A publication Critical patent/JP2021177830A/en
Application granted granted Critical
Publication of JP7451040B2 publication Critical patent/JP7451040B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

この発明は、生体状態推定装置に関し、具体的には、レーダーセンサを用いて生体の移動状態、安静状態、および不在状態を推定する技術に関する。 The present invention relates to a biological state estimating device, and specifically relates to a technique for estimating the moving state, resting state, and absent state of a living body using a radar sensor.

従来の人感センサとして、焦電センサや赤外線センサ、カメラが主に用いられており、生体の動作を検出して見守りなどが行われる。具体的には例えば、介護施設などにおいて利用者の見守りを行う技術として、各部屋に設置された利用者用センサと、各部屋に設置された利用者用ビデオカメラと、利用者の情報を表示する介護者用端末とを備える見守りシステムが知られている(特許文献1参照)。 Pyroelectric sensors, infrared sensors, and cameras are mainly used as conventional human sensors to detect movements of living organisms and monitor them. Specifically, for example, as a technology for monitoring users in nursing care facilities, etc., we use sensors for users installed in each room, video cameras for users installed in each room, and display of user information. A monitoring system is known that includes a terminal for a caregiver (see Patent Document 1).

特開2018-132795号公報Japanese Patent Application Publication No. 2018-132795

ところで、焦電センサや赤外線センサでは、ドアの開閉や生体の歩行などによって受信レベルが著しく変化することを利用して存在の有無を判別している。この際、受信レベル(別言すると、振幅)に関する閾値を予め設定してこの閾値と実際に観測される受信レベル(振幅)とを比較することによって判別するため、受信レベル(振幅)に関する閾値の設定が困難であり、結果として、動作の小さい安静状態のときに不在状態と誤認識するケースが多々ある、という問題がある。さらに、センサを使用する環境ごとに手動で受信レベル(振幅)に関する閾値を設定する必要があり、手間がかかるとともに、使用環境が変化するたびに閾値の設定を変更しないと誤認識してしまう、という問題がある。また、特許文献1に記載の技術のようにカメラを用いる場合には、プライバシーが侵害される、という問題があり、さらに、映像では睡眠時と生体の危機的状況とを判別することができない、という問題がある。 By the way, in pyroelectric sensors and infrared sensors, presence or absence is determined by utilizing the fact that the reception level changes significantly due to the opening/closing of a door or the walking of a living body. At this time, the threshold regarding the reception level (amplitude) is set in advance and the determination is made by comparing this threshold with the actually observed reception level (amplitude). The setting is difficult, and as a result, there are many cases where a resting state with little movement is mistakenly recognized as an absent state. Furthermore, it is necessary to manually set the reception level (amplitude) threshold for each environment in which the sensor is used, which is time-consuming and can lead to erroneous recognition unless the threshold setting is changed each time the usage environment changes. There is a problem. Furthermore, when a camera is used as in the technique described in Patent Document 1, there is a problem that privacy is violated, and furthermore, it is not possible to distinguish between sleeping and a critical state of the living body using images. There is a problem.

そこでこの発明は、プライバシーを侵害することなく、また、使用環境ごとに閾値を手動で設定することなく、生体の移動状態、安静状態、および不在状態を高精度に判別することが可能な、生体状態推定装置を提供することを目的とする。 Therefore, the present invention aims to provide a biological system that can accurately determine the moving state, resting state, and absent state of a living body without violating privacy or manually setting thresholds for each usage environment. The purpose of the present invention is to provide a state estimation device.

上記課題を解決するために、請求項1に記載の発明は、送信電波を放射して前記送信電波の反射波を受信するとともに前記送信電波と前記反射波とを混合してビート信号を生成する送受信部と、前記送受信部から出力される前記ビート信号の位相を計算する位相計算部と、前記送受信部から出力される前記ビート信号のデータと前記位相計算部から出力される前記位相のデータとのそれぞれから静止物からの反射波に相当する成分を除去する静止物反射除去部と、前記送受信部から出力される前記ビート信号のデータと前記位相計算部から出力される前記位相のデータとのそれぞれから体動に相当する成分を除去する体動成分除去部と、前記2つの成分の除去後の前記位相のデータを用いて前記送受信部と物標としての生体との相対速度を算出する速度算出部と、前記相対速度が前記生体の呼吸の速度よりも大きい場合に前記生体が移動状態である判断する移動状態推定部と、前記相対速度が前記生体の前記呼吸の速度以下であり、且つ、前記2つの成分の除去後の前記ビート信号のデータを高速フーリエ変換して得られる周波数スペクトルの分布のうちの特定の周波数範囲のスペクトル強度の合計前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計よりも大きい場合に前記生体が安静状態である判断する安静状態推定部と、前記特定の周波数範囲のスペクトル強度の合計が前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計以下である場合に前記2つの成分の除去後の前記ビート信号の振幅の大きさ別の度数分布形が正規分布であるか否かに基づいて前記生体が不在状態であるか否かを判断する不在状態推定部と、を有前記特定の周波数範囲が、0Hz~5Hz未満に設定される、ことを特徴とする生体状態推定装置である。 In order to solve the above problem, the invention according to claim 1 emits a transmitted radio wave, receives a reflected wave of the transmitted radio wave, and mixes the transmitted radio wave and the reflected wave to generate a beat signal. a transmitting/receiving section; a phase calculating section that calculates the phase of the beat signal output from the transmitting/receiving section; data of the beat signal output from the transmitting/receiving section; and data of the phase output from the phase calculating section; a stationary object reflection removing section that removes a component corresponding to a reflected wave from a stationary object from each of the beat signal data output from the transmitting/receiving section and the phase data output from the phase calculating section. a body movement component removing unit that removes a component corresponding to body movement from each; and a speed that calculates a relative velocity between the transmitting and receiving unit and a living body as a target using the phase data after removing the two components. a calculation unit; a moving state estimation unit that determines that the living body is in a moving state when the relative velocity is greater than the breathing rate of the living body; and the relative velocity is less than or equal to the breathing rate of the living body; and the sum of the spectral intensities in a specific frequency range of the frequency spectrum distribution obtained by fast Fourier transforming the data of the beat signal after the removal of the two components is within a predetermined region other than the specific frequency range. a resting state estimation unit that determines that the living body is in a resting state when the sum of spectral intensities in the specific frequency range is greater than the sum of spectral intensities in the specific frequency range; and a spectral intensity in a predetermined region other than the specific frequency range. If the sum is less than or equal to the sum of and an absent state estimating unit that makes a determination , and the specific frequency range is set from 0 Hz to less than 5 Hz .

請求項2に記載の発明は、請求項1に記載の生体状態推定装置において、前記生体の前記呼吸の速度が、1cm~10cm/秒に設定される、ことを特徴とする。 The invention according to claim 2 is characterized in that, in the biological state estimating device according to claim 1 , the respiration rate of the living body is set to 1 cm to 10 cm/sec.

請求項に記載の発明は、請求項1または2に記載の生体状態推定装置において、前記生体の状態の推定結果が、前記安静状態から、前記移動状態を経ることなく、前記不在状態へと推移した場合に、前記生体が危険状態にあると判断する危険状態検知部をさらに有する、ことを特徴とする。 According to a third aspect of the invention, in the biological state estimation device according to the first or second aspect, the estimation result of the biological state is changed from the resting state to the absent state without passing through the moving state. The present invention is characterized in that it further includes a dangerous state detection unit that determines that the living body is in a dangerous state when the living body changes.

請求項1及び2に記載の発明によれば、例えば壁や床などの静止物からの反射波に相当する成分を自律的に除去するようにしているので、使用環境ごとに閾値を設定することなく、言い換えると、使用環境に影響を受けることなく、また、使用環境に依存することなく、生体の不在状態、安静状態、および移動状態を高精度に判別することが可能となる。請求項1に記載の発明によれば、また、移動状態、安静状態、そして不在状態というように動作の大きい順に段階的に状態を区分して特定するようにしているので、無駄なく効率的に且つ精度高く、生体の状態を推定することが可能となる。 According to the inventions described in claims 1 and 2 , components corresponding to reflected waves from stationary objects such as walls and floors are autonomously removed, so a threshold value can be set for each usage environment. In other words, it is possible to determine with high precision whether the living body is absent, at rest, or moving, without being affected by or depending on the usage environment. According to the invention set forth in claim 1, the states are classified and identified step by step in descending order of movement, such as a moving state, a resting state, and an absent state. Moreover, it becomes possible to estimate the state of the living body with high accuracy.

また、請求項1及び請求項2に記載の発明によれば、観測対象の生体が人である場合に特に、移動状態であるか否かを適切に判断することが可能となる。 Moreover, according to the invention described in claims 1 and 2, it becomes possible to appropriately determine whether or not the living body to be observed is in a moving state, especially when the living body to be observed is a human.

また、請求項1及び請求項に記載の発明によれば、観測対象の生体が人である場合に特に、安静状態であるか否かを適切に判断することが可能となる。 Moreover, according to the inventions described in claims 1 and 2 , it becomes possible to appropriately determine whether or not the living body to be observed is in a resting state, especially when the living body to be observed is a human being.

請求項に記載の発明によれば、生体が移動状態、安静状態、および不在状態のうちのいずれであるかに加えて生体が危険状態にあることを検知することが可能となる。 According to the third aspect of the invention, it is possible to detect whether the living body is in a moving state, a resting state, or an absent state, as well as whether the living body is in a dangerous state.

この発明の実施の形態に係る生体状態推定装置の概略構成を示す機能ブロック図である。1 is a functional block diagram showing a schematic configuration of a biological state estimation device according to an embodiment of the present invention. 図1の生体状態推定装置における処理手順を示すフロー図である。FIG. 2 is a flow diagram showing a processing procedure in the biological state estimation device of FIG. 1. FIG. 位相変位を速度に変換する処理、および、生体が移動状態であるか否かを判断する処理を説明する図である。FIG. 3 is a diagram illustrating a process of converting phase displacement into velocity and a process of determining whether a living body is in a moving state. 生体が安静状態であるか否かを判断する処理を説明する図である。(A)は生体が不在状態であるときの周波数スペクトルの分布を示す図である。(B)は生体が安静状態であるときの周波数スペクトルの分布を示す図である。It is a figure explaining the process of determining whether a living body is in a resting state. (A) is a diagram showing a frequency spectrum distribution when a living body is absent. (B) is a diagram showing the frequency spectrum distribution when the living body is in a resting state. 生体が不在状態であるか否かを判断する処理を説明する図である。It is a figure explaining the process which judges whether a living body is absent.

以下、この発明を図示の実施の形態に基づいて説明する。図1は、この発明の実施の形態に係る生体状態推定装置1の概略構成を示す機能ブロック図である。また、図2は、実施の形態に係る生体状態推定装置1における処理手順を示すフロー図である。生体状態推定装置1は、レーダーセンサを用いて生体の移動状態、安静状態、および不在状態を推定するための仕組みであり、送受信部2と、推定装置3と、を有する。なお、生体としては、例えば、人や、飼育されたり保護されたりしている種々の動物が挙げられる。 The present invention will be described below based on the illustrated embodiments. FIG. 1 is a functional block diagram showing a schematic configuration of a biological condition estimation device 1 according to an embodiment of the present invention. Further, FIG. 2 is a flow diagram showing a processing procedure in the biological state estimation device 1 according to the embodiment. The biological state estimating device 1 is a mechanism for estimating the moving state, resting state, and absent state of a living body using a radar sensor, and includes a transmitting/receiving section 2 and an estimating device 3. Note that living organisms include, for example, humans and various animals that are kept or protected.

そして、この実施の形態に係る生体状態推定装置1は、送信電波を放射して送信電波の反射波を受信するとともに送信電波と反射波とを混合してビート信号を生成する送受信部2と、送受信部2から出力されるビート信号の位相を計算する位相計算部(位相計算タスク33a)と、送受信部2から出力されるビート信号のデータと位相計算部(位相計算タスク33a)から出力される位相のデータとのそれぞれから静止物からの反射波に相当する成分を除去する静止物反射除去部(静止物反射除去タスク33b)と、送受信部2から出力されるビート信号のデータと位相計算部(位相計算タスク33a)から出力される位相のデータとのそれぞれから体動に相当する成分を除去する体動成分除去部(体動成分除去タスク33c)と、前記2つの成分の除去後の位相のデータを用いて送受信部2と物標としての生体との相対速度を算出する速度算出部(速度算出タスク33d)と、相対速度に基づいて生体が移動状態であるか否かを判断する移動状態推定部(移動状態推定タスク33e)と、前記2つの成分の除去後のビート信号のデータを高速フーリエ変換して得られる周波数スペクトルの分布のうちの特定の周波数範囲のスペクトル強度の合計と前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計との比較に基づいて生体が安静状態であるか否かを判断する安静状態推定部(安静状態推定タスク33g)と、前記2つの成分の除去後のビート信号の振幅の大きさ別の度数分布形が正規分布であるか否かに基づいて生体が不在状態であるか否かを判断する不在状態推定部(不在状態推定タスク33h)と、を有する、ようにしている。 The biological condition estimation device 1 according to this embodiment includes a transmitting/receiving section 2 that emits a transmitted radio wave, receives a reflected wave of the transmitted radio wave, and mixes the transmitted radio wave and the reflected wave to generate a beat signal. A phase calculation unit (phase calculation task 33a) that calculates the phase of the beat signal output from the transmission and reception unit 2, and data of the beat signal output from the transmission and reception unit 2 and output from the phase calculation unit (phase calculation task 33a) a stationary object reflection removal unit (stationary object reflection removal task 33b) that removes a component corresponding to a reflected wave from a stationary object from the phase data, and a phase calculation unit and the beat signal data output from the transmitting/receiving unit 2. A body motion component removal unit (body motion component removal task 33c) that removes a component corresponding to body motion from each of the phase data output from the phase calculation task 33a, and the phase after the removal of the two components. a speed calculation unit (speed calculation task 33d) that calculates the relative speed between the transmitting/receiving unit 2 and the living body as a target using data; and a movement unit that determines whether or not the living body is in a moving state based on the relative speed. A state estimating unit (moving state estimation task 33e) and a sum of spectrum intensities in a specific frequency range of a frequency spectrum distribution obtained by performing fast Fourier transform on the data of the beat signal after removal of the two components, and the above. a resting state estimation unit (resting state estimation task 33g) that determines whether the living body is in a resting state based on a comparison with the sum of spectral intensities in a predetermined region other than a specific frequency range; an absent state estimation unit (absent state estimation task 33h) that determines whether or not the living body is in the absent state based on whether the frequency distribution shape of each amplitude of the removed beat signal is a normal distribution; , have, like.

この実施の形態に係る生体状態推定装置1は、さらに、生体の状態の推定結果が、安静状態から、移動状態を経ることなく、不在状態へと推移した場合に、生体が危険状態にあると判断する危険状態検知部(危険状態検知タスク33i)を有する、ようにしている。 The living body state estimation device 1 according to this embodiment further determines that the living body is in a dangerous state when the estimation result of the living body's state changes from a resting state to an absent state without passing through a moving state. It has a dangerous state detection unit (dangerous state detection task 33i) that makes a judgment.

送受信部2は、所定の領域についてレーダースキャンを行うとともに送信波と反射波とを混合してビート信号を生成するための機序であり、送信波としてマイクロ波を照射する送信部21と、前記マイクロ波の反射波を検出する受信部22とを備える。送信部21から照射されるマイクロ波の周波数帯は、例えば24GHz帯が挙げられるものの、24GHz帯には限定されない。また、前記における所定の領域は、観測対象の生体が居る空間であり、例えば居住空間や飼育スペース、保護区画などの生体の生活空間/活動空間が挙げられる。 The transmitting/receiving unit 2 is a mechanism for performing a radar scan on a predetermined area and mixing a transmitted wave and a reflected wave to generate a beat signal. It includes a receiving section 22 that detects reflected waves of microwaves. The frequency band of the microwave irradiated from the transmitter 21 is, for example, a 24 GHz band, but is not limited to the 24 GHz band. Further, the predetermined area mentioned above is a space where a living organism to be observed is present, and examples thereof include a living space/activity space of a living organism such as a living space, a rearing space, and a protection area.

送信部21は、FMCW(Frequency Modulated Continuous Wave の略;周波数変調連続波)方式のレーダーを送信用アンテナ211を介して送信する。FMCW方式のレーダーは、時間の経過に応じて周波数変調を行う連続発振レーダーであり、複数のチャープを含むバースト波を生成して照射する。バースト波に含まれる各チャープは、周波数が時間的に掃引されて(具体的には、時間経過に伴って周波数が線形的に変化するように)生成される。 The transmitting unit 21 transmits an FMCW (Frequency Modulated Continuous Wave) radar via a transmitting antenna 211 . The FMCW radar is a continuous wave radar that performs frequency modulation over time, and generates and irradiates burst waves containing multiple chirps. Each chirp included in the burst wave is generated with the frequency swept over time (specifically, the frequency changes linearly over time).

送信部21は、例えば、時間軸上で周波数が次第に高くなる変調区間と周波数が次第に低くなる変調区間とが交互に連続的に設けられて三角波形状に変化する送信信号を生成するための制御電圧値のデジタルデータをデジタル-アナログ変換して制御電圧信号を生成するとともに、前記制御電圧信号に応じて送信信号を生成して送信用アンテナ211へと伝送する。また、送信信号の一部が、所定の分配比で、ローカル信号として受信部22へも伝送される。 The transmitting unit 21 uses a control voltage for generating a transmission signal that changes in the shape of a triangular wave by alternately and continuously providing modulation sections in which the frequency gradually increases and modulation sections in which the frequency gradually decreases on the time axis, for example. A control voltage signal is generated by digital-to-analog conversion of the digital data of the value, and a transmission signal is generated according to the control voltage signal and transmitted to the transmission antenna 211. Further, a part of the transmission signal is also transmitted to the receiving section 22 as a local signal at a predetermined distribution ratio.

そして、送信部21の送信用アンテナ211は、伝送される送信信号に基づく送信電波を放射する。 The transmitting antenna 211 of the transmitter 21 emits a transmitting radio wave based on the transmitted transmitting signal.

FMCW方式のレーダーでは、送信波と受信波(即ち、反射波)との差分を解析することにより、送受信部2と物標との間の距離や送受信部2と物標との相対速度などが計算される。なお、チャープの周波数の変調幅や変調の周期(即ち、チャープの繰返し周期)は適宜調節されるようにしてよく、また、複数のチャープをすべて同一波形とするか一部のチャープを異なる波形とするかが適宜選択されるようにしてよい。 In the FMCW radar, by analyzing the difference between the transmitted wave and the received wave (i.e., the reflected wave), the distance between the transmitting/receiving unit 2 and the target object, the relative speed between the transmitting/receiving unit 2 and the target object, etc. can be determined. Calculated. Note that the modulation width and modulation period of the chirp frequency (that is, the repetition period of the chirp) may be adjusted as appropriate, and it is also possible to make all the chirps have the same waveform or to make some of the chirps have different waveforms. It may be possible to select as appropriate whether to do so or not.

受信部22は、送信部21の送信用アンテナ211から放射されるFMCW方式のレーダーの反射波(具体的には、反射電波)を受信用アンテナ221を介して受信するとともに、反射波のToF(Time of Flight の略)に比例する中間周波数へのダウンコンバート処理を施した信号を生成して出力する。 The receiving unit 22 receives reflected waves (specifically, reflected radio waves) of the FMCW radar radiated from the transmitting antenna 211 of the transmitting unit 21 via the receiving antenna 221, and also receives the ToF (ToF) of the reflected waves. It generates and outputs a signal that has been down-converted to an intermediate frequency proportional to Time of Flight (abbreviation for Time of Flight).

例えば、受信部22の受信用アンテナ221は、送信用アンテナ211から放射される送信電波が反射して戻ってくる反射電波を受信して受信信号を生成して出力する。受信信号は、同レベル比で分割されて、後の処理へと供給される。 For example, the reception antenna 221 of the reception unit 22 receives reflected radio waves that are reflected back from the transmission radio waves radiated from the transmission antenna 211, generates and outputs a reception signal. The received signal is divided at the same level ratio and provided for subsequent processing.

ここで、送信部21から受信部22へと伝送されるローカル信号は、同レベル比で分割され、分割されたうちの一方はそのまま後の処理へと供給され、分割されたうちの他方はπ/2の位相回転処理が施されたうえで後の処理へと供給される。 Here, the local signal transmitted from the transmitter 21 to the receiver 22 is divided at the same level ratio, one of the divided signals is supplied as is to subsequent processing, and the other of the divided signals is π The signal is subjected to phase rotation processing of /2 and then supplied to subsequent processing.

そして、受信部22は、受信用アンテナ221から出力される受信信号のうちの分割された一方の信号と、送信部21から伝送されるローカル信号のうちの分割された一方の信号(即ち、π/2位相回転処理が施されていない方の信号)とを混合して、中間周波数信号(IF信号:Intermediate Frequency signal)であるビート信号を、具体的には、送信信号と同相(In Phase)成分のI信号を生成する。 Then, the receiving section 22 receives one of the divided signals of the received signal output from the receiving antenna 221 and one of the divided signals of the local signal transmitted from the transmitting section 21 (i.e., π /2 signal that has not been subjected to phase rotation processing) to produce a beat signal that is an intermediate frequency signal (IF signal), specifically, in phase with the transmission signal. Generate component I signals.

受信部22は、また、受信用アンテナ221から出力される受信信号のうちの分割された他方の信号と、送信部21から伝送されるローカル信号のうちの分割されたうえでπ/2位相回転処理が施された方の信号とを混合して、中間周波数信号(IF信号)であるビート信号を、具体的には、送信信号に対して直交位相(Quadrature Phase)成分のQ信号を生成する。 The receiving unit 22 also performs π/2 phase rotation on the other divided signal of the received signal output from the receiving antenna 221 and the divided local signal transmitted from the transmitting unit 21. The processed signal is mixed to generate a beat signal that is an intermediate frequency signal (IF signal), specifically, a Q signal that is a quadrature phase component with respect to the transmitted signal. .

受信部22は、続いて、I信号に対して、必要に応じて不要な周波数成分の除去処理や信号レベルの増幅処理を施したうえで、アナログ-デジタル変換処理を施して、デジタル信号化したI信号を生成して出力する。 The receiving unit 22 then performs unnecessary frequency component removal processing and signal level amplification processing on the I signal as necessary, and then performs analog-to-digital conversion processing to convert the I signal into a digital signal. Generate and output an I signal.

受信部22は、また、Q信号に対して、必要に応じて不要な周波数成分の除去処理や信号レベルの増幅処理を施したうえで、アナログ-デジタル変換処理を施して、デジタル信号化したQ信号を生成して出力する。 The receiving unit 22 also performs unnecessary frequency component removal processing and signal level amplification processing on the Q signal as necessary, and then performs analog-to-digital conversion processing to convert the Q signal into a digital signal. Generate and output a signal.

なお、ToFは、送信されたレーダーが反射波として受信されるまでの時間であり、つまり送受信部2とレーダーを反射させた物との間の距離に比例する値である。ToFを用いることにより、送受信部2とレーダーを反射させた物(物標)との間の距離やレーダーを反射させた物(物標)に纏わる速度が計算される。 Note that ToF is the time it takes for the transmitted radar to be received as a reflected wave, that is, it is a value proportional to the distance between the transmitter/receiver 2 and the object that reflected the radar. By using ToF, the distance between the transmitting/receiving unit 2 and the object (target) that reflected the radar and the speed of the object (target) that reflected the radar are calculated.

送受信部2によるレーダースキャンは、所定の時間間隔で連続的に繰り返し行われる。そして、受信部22から、デジタル信号化されたビート信号(具体的には、I信号とQ信号との組み合わせ)が、前記所定の時間間隔で連続して継続的に出力される(ステップS1)。 Radar scanning by the transmitting/receiving unit 2 is continuously and repeatedly performed at predetermined time intervals. Then, a beat signal converted into a digital signal (specifically, a combination of an I signal and a Q signal) is continuously output from the receiving unit 22 at the predetermined time interval (step S1). .

推定装置3は、送受信部2から出力されるビート信号に基づいて生体の状態を推定するための機序であり、主として、入力受部31と、記憶部32と、メインタスク33と、中央処理部34と、を備える。 The estimation device 3 is a mechanism for estimating the state of the living body based on the beat signal output from the transmitting/receiving section 2, and mainly includes an input receiving section 31, a storage section 32, a main task 33, and a central processing unit. A section 34 is provided.

中央処理部34は、推定装置3を構成する各部を統制して制御などする機能を備え、例えば、CPU(Central Processing Unit の略)などを用いたプロセッサによって構成されたり、記憶部32に格納されたプログラムに従って各機能を実現するものとして構成されたりする。 The central processing unit 34 has a function of regulating and controlling each unit constituting the estimation device 3, and is configured by, for example, a processor using a CPU (abbreviation for Central Processing Unit), or is stored in the storage unit 32. It may be configured to realize each function according to a set program.

入力受部31は、受信部22から出力されるビート信号(具体的には、I信号とQ信号との組み合わせ)の入力を受け、入力された前記ビート信号(I信号とQ信号との組み合わせ;図2においては「IQデータ」)を記憶部32に記憶させる(ステップS2)。 The input receiving section 31 receives the input of the beat signal (specifically, the combination of the I signal and the Q signal) output from the receiving section 22, and receives the input beat signal (the combination of the I signal and the Q signal). ; "IQ data" in FIG. 2) is stored in the storage unit 32 (step S2).

記憶部32は、例えばハードディスクやメモリであり、各種の情報、プログラム、およびデータなどを記憶する記憶領域/記憶装置として機能する。記憶部32として、生体の状態の推定に纏わる処理を実行する際の作業領域として機能するRAM(Random Access Memory の略)が含まれるようにしてもよい。 The storage unit 32 is, for example, a hard disk or a memory, and functions as a storage area/storage device that stores various information, programs, data, and the like. The storage unit 32 may include a RAM (abbreviation for Random Access Memory) that functions as a work area when executing processes related to estimation of the state of the living body.

メインタスク33は、記憶部32に記憶されているビート信号(I信号とQ信号との組み合わせ)を用いて、生体の状態の推定に纏わる処理を実行するためのタスク・プログラム群である。このメインタスク33は、主として、位相計算タスク33a、静止物反射除去タスク33b、体動成分除去タスク33c、速度算出タスク33d、移動状態推定タスク33e、FFT処理タスク33f、安静状態推定タスク33g、不在状態推定タスク33h、および危険状態検知タスク33iを備える。 The main task 33 is a group of tasks and programs for executing processing related to estimating the state of the living body using the beat signal (combination of I signal and Q signal) stored in the storage unit 32. The main tasks 33 mainly include a phase calculation task 33a, a stationary object reflection removal task 33b, a body motion component removal task 33c, a speed calculation task 33d, a moving state estimation task 33e, an FFT processing task 33f, a resting state estimation task 33g, and an absent task. It includes a state estimation task 33h and a dangerous state detection task 33i.

位相計算タスク33aは、ビート信号の位相を計算するためのタスク・プログラムである。位相計算タスク33aは、具体的には、記憶部32に記憶されているビート信号(I信号とQ信号との組み合わせ)を読み込み、相互に対応するI信号とQ信号との組み合わせに基づいて直交検波を行って前記組み合わせそれぞれの位相を計算する(ステップS3)。 The phase calculation task 33a is a task program for calculating the phase of a beat signal. Specifically, the phase calculation task 33a reads the beat signal (combination of I signal and Q signal) stored in the storage unit 32, and performs orthogonal calculation based on the combination of the mutually corresponding I signal and Q signal. Detection is performed to calculate the phase of each of the combinations (step S3).

位相計算タスク33aは、計算した位相を記憶部32に記憶させる。 The phase calculation task 33a causes the storage unit 32 to store the calculated phase.

静止物反射除去タスク33bは、受信信号に対して信号処理を施して得られる信号から、例えば壁や床などの静止物からの反射波に相当する成分を除去するためのタスク・プログラムである。 The stationary object reflection removal task 33b is a task program for removing components corresponding to reflected waves from stationary objects such as walls and floors, for example, from a signal obtained by performing signal processing on a received signal.

ここで、受信用アンテナ221から出力される、反射波に対応する信号(即ち、受信信号)には、静止物からの反射波、生体からの反射波、およびノイズが含まれる。静止物反射除去タスク33bは、生体からの反射波を良好に検出して生体の状態を精度よく推定するために、前記のうちの静止物からの反射波に相当する成分を除去する。前記における静止物は、観測対象の生体が居る空間に存在する固定物であり、例えば前記空間を囲む壁や床などが挙げられる。 Here, the signal corresponding to the reflected wave (ie, the received signal) output from the reception antenna 221 includes a reflected wave from a stationary object, a reflected wave from a living body, and noise. The stationary object reflection removal task 33b removes the component corresponding to the reflected wave from the stationary object, in order to detect the reflected wave from the living body well and estimate the state of the living body with high accuracy. The stationary object mentioned above is a fixed object that exists in the space where the living organism to be observed resides, and includes, for example, a wall or a floor surrounding the space.

静止物からの反射波は、直流波である(もしくは、殆ど直流波である)と仮定すると、ステップS3の処理によって計算される位相の時間平均として把握される。そこで、静止物反射除去タスク33bは、記憶部32に記憶されている位相のデータを読み込んで所定の分だけ蓄積するとともに蓄積した位相の時間平均を計算し、計算した時間平均に相当する信号を、ステップS3の処理によって計算される位相のそれぞれから差し引くことによって静止物からの反射波に相当する成分を除去する(ステップS4)。 Assuming that the reflected wave from a stationary object is a DC wave (or is almost a DC wave), it is understood as a time average of the phase calculated by the process in step S3. Therefore, the stationary object reflection removal task 33b reads the phase data stored in the storage unit 32, accumulates it by a predetermined amount, calculates the time average of the accumulated phase, and generates a signal corresponding to the calculated time average. , components corresponding to reflected waves from stationary objects are removed by subtracting them from each of the phases calculated by the processing in step S3 (step S4).

静止物からの反射波に相当する成分を除去するための、位相のデータの蓄積の程度は、特定の時間長さあるいはデータ数に限定されるものではなく、静止物からの反射波に相当する成分を良好に除去することができるような適切な平均が計算され得る適当な時間長さやデータ数に適宜設定される。 The degree of accumulation of phase data to remove components corresponding to reflected waves from stationary objects is not limited to a specific time length or number of data, but is equivalent to reflected waves from stationary objects. The time length and number of data are appropriately set so that an appropriate average can be calculated so that the components can be removed satisfactorily.

静止物反射除去タスク33bは、静止物からの反射波に相当する成分の除去処理後の位相のデータ(「第1の除去処理後位相データ」と呼ぶ)を記憶部32に記憶させる。 The stationary object reflection removal task 33b causes the storage unit 32 to store phase data after removal processing of a component corresponding to a reflected wave from a stationary object (referred to as "first removed phase data").

体動成分除去タスク33cは、受信信号に対して信号処理を施して得られる信号から、観測対象の生体の意識とは関係なく起きる低周期の体全体の揺れ(「体動」とも呼ばれる)に相当する成分を除去するためのタスク・プログラムである。 The body motion component removal task 33c extracts low-period shaking of the entire body (also called "body motion") that occurs regardless of the consciousness of the living organism being observed, from a signal obtained by performing signal processing on the received signal. This is a task program for removing corresponding components.

体動に相当する成分を除去する仕法は、特定の手法や手順に限定されるものではなく、位相データから体動に相当する成分を除去(別言すると、抑圧)し得る手法や手順の中から適当なものが適宜選択される。具体的には例えば、位相データを所定の時間長さで区切り、区切られた区間それぞれのデータのN次関数近似曲線(但し、Nは自然数)を推定することによって体動に相当する成分を表す回帰曲線を求め、求めた回帰曲線に相当する信号を、位相データから差し引く手法が用いられ得る。 The method of removing the component corresponding to body movement is not limited to a specific method or procedure, but is any method or procedure that can remove (in other words, suppress) the component corresponding to body movement from phase data. An appropriate one is selected from among them. Specifically, for example, phase data is divided into predetermined time lengths, and components corresponding to body movements are expressed by estimating an N-order function approximation curve (where N is a natural number) of the data in each divided section. A method may be used in which a regression curve is determined and a signal corresponding to the determined regression curve is subtracted from the phase data.

体動成分除去タスク33cは、記憶部32に記憶されている第1の除去処理後位相データを読み込み、読み込んだ前記第1の除去処理後位相データから体動に相当する成分を除去する(ステップS5)。 The body motion component removal task 33c reads the first removed phase data stored in the storage unit 32, and removes a component corresponding to body movement from the read first removed phase data (step S5).

体動成分除去タスク33cは、体動に相当する成分の除去処理後の位相データ(「第2の除去処理後位相データ」と呼ぶ)を記憶部32に記憶させる。 The body motion component removal task 33c causes the storage unit 32 to store phase data after the removal process of the component corresponding to the body movement (referred to as "second post-removal process phase data").

速度算出タスク33dは、生体に纏わる速度を算出するためのタスク・プログラムである。 The speed calculation task 33d is a task program for calculating the speed related to a living body.

ここで、生体からの反射波は、生体の腕や脚の動作、呼吸、および拍動などの生体の動きに応じて位相や振幅が変化し、生体の腕や脚の動作に対応する反射波、生体の呼吸の動きに対応する反射波、さらに、生体の拍動に対応する反射波を含む。そして、生体からの反射波の位相や振幅の変化に伴って受信信号の位相や振幅が変化する。 Here, the reflected waves from the living body change in phase and amplitude depending on the movement of the living body, such as the movements of the arms and legs of the living body, breathing, and heartbeat, and the reflected waves corresponding to the movements of the arms and legs of the living body change. , a reflected wave corresponding to the respiratory movement of the living body, and a reflected wave corresponding to the pulsation of the living body. Then, the phase and amplitude of the received signal change with the change in the phase and amplitude of the reflected wave from the living body.

速度算出タスク33dは、記憶部32に記憶されている第2の除去処理後位相データを読み込み、位相変位を速度に変換して(図3参照;左の図の2本の線は、I信号とQ信号との各々の振幅の時間変化を表す)、送受信部2と物標(ここでは、特に生体)との相対速度を算出する(ステップS6)。 The speed calculation task 33d reads the phase data after the second removal process stored in the storage unit 32, converts the phase displacement into speed (see FIG. 3; the two lines in the left diagram indicate the I signal and the Q signal), and calculate the relative velocity between the transmitter/receiver 2 and the target object (here, particularly a living body) (step S6).

なお、送受信部2と物標としての生体との相対速度の算出の仕法は、特定の手法や手順に限定されるものではなく、生体の腕や脚の動作、呼吸、および拍動などの生体の動きに纏わる速度を計算することができる手法や手順の中から適当なものが適宜選択される。 Note that the method for calculating the relative velocity between the transmitting/receiving unit 2 and a living body as a target is not limited to a specific method or procedure, and is based on the movements of the arms and legs of the living body, breathing, and heartbeat. An appropriate method or procedure is selected as appropriate from among the methods and procedures that can calculate the speed associated with the movement of.

速度算出タスク33dは、算出した相対速度を移動状態推定タスク33eへと転送する。 The speed calculation task 33d transfers the calculated relative speed to the movement state estimation task 33e.

移動状態推定タスク33eは、速度算出タスク33dによる相対速度の算出結果に基づいて、物標としての生体が移動しているか否かを判断するためのタスク・プログラムである。移動状態推定タスク33eは、速度算出タスク33dから転送される相対速度に基づいて、物標としての生体が移動しているか否かを判断する(ステップS7)。 The movement state estimation task 33e is a task program for determining whether or not a living body serving as a target is moving, based on the relative speed calculation result by the speed calculation task 33d. The movement state estimation task 33e determines whether or not the living body serving as the target is moving based on the relative speed transferred from the speed calculation task 33d (step S7).

移動状態推定タスク33eは、具体的には、速度算出タスク33dから転送される相対速度の値を確認し、生体の呼吸の速度を上回る速度が検出された場合に、物標としての生体が移動していると判断する(図3参照)。生体の呼吸であるのか移動であるのかを推定/判別するための呼吸の速度は、特定の値に限定されるものではなく、公知の知見などに基づいて適当な値に適宜設定される。呼吸の速度は、具体的には例えば、人の呼吸について、苦しいときには腹部が大凡1秒に1回、最大で5cm程度の幅で拡縮するので、1cm~10cm/秒程度に設定されることが考えられ、5cm/秒程度に設定されることが好ましい。 Specifically, the movement state estimation task 33e checks the value of the relative velocity transferred from the speed calculation task 33d, and if a speed exceeding the respiration rate of the living body is detected, the movement state estimation task 33e determines that the living body as a target is moving. (See Figure 3). The respiration rate for estimating/discriminating whether it is respiration of a living body or movement is not limited to a specific value, but is appropriately set to an appropriate value based on known knowledge. Specifically, the breathing rate can be set to about 1 cm to 10 cm/sec, for example, because when a person's breathing is difficult, the abdomen expands and contracts approximately once every second, with a maximum width of about 5 cm. It is preferable to set the speed to about 5 cm/sec.

なお、呼吸の速度は、後のステップS13の処理とも関係して、移動状態と安静状態とを区別するための速度でもある。呼吸の速度は、このため、安静状態における呼吸の速度と、例えば寝返りの速度とを区別できるように設定される。具体的には例えば、呼吸の速度が5cm/秒程度に設定されると、寝返りのような生体の動きを移動と判断し、一方で、安静状態における呼吸に伴う生体の動きを移動とは判断しないようにすることができる。 Note that the breathing rate is also related to the process in step S13 later, and is also a rate for distinguishing between a moving state and a resting state. The rate of respiration is therefore set in such a way that the rate of respiration in a resting state can be distinguished from, for example, the rate of turning over. Specifically, for example, when the breathing rate is set to about 5 cm/sec, movements of the living body such as turning over are judged as movement, while movements of the living body accompanying breathing in a resting state are judged as movement. You can prevent it from happening.

速度算出タスク33dから転送される相対速度が呼吸の速度よりも大きい場合には(ステップS7:Yes)、移動状態推定タスク33eは、生体が移動している(移動状態である)と判断する(ステップS8)。そして、メインタスク33は、処理手順をステップS19へと進める。 If the relative velocity transferred from the velocity calculation task 33d is greater than the breathing velocity (step S7: Yes), the movement state estimation task 33e determines that the living body is moving (in a moving state). Step S8). Then, the main task 33 advances the processing procedure to step S19.

一方、速度算出タスク33dから転送される相対速度が呼吸の速度以下である場合には(ステップS7:No)、移動状態推定タスク33eは、生体が移動していないと判断し、処理手順をステップS9へと進める。 On the other hand, if the relative velocity transferred from the velocity calculation task 33d is less than or equal to the breathing velocity (step S7: No), the movement state estimation task 33e determines that the living body is not moving and moves the processing procedure to step S7. Proceed to S9.

ステップS9以降の処理では、ステップS2の処理において記憶部32に記憶されるビート信号(I信号とQ信号との組み合わせ:IQデータ)が用いられる。 In the processing after step S9, the beat signal (combination of I signal and Q signal: IQ data) stored in the storage unit 32 in the processing in step S2 is used.

静止物反射除去タスク33bは、ビート信号(I信号とQ信号との組み合わせ:IQデータ)を記憶部32から読み込み(ステップS9)、I信号およびQ信号のそれぞれから、例えば壁や床などの静止物からの反射波に相当する成分を除去する。 The stationary object reflection removal task 33b reads the beat signal (combination of I signal and Q signal: IQ data) from the storage unit 32 (step S9), and extracts the beat signal (combination of I signal and Q signal: IQ data) from each of the I signal and Q signal. Removes components corresponding to reflected waves from objects.

静止物からの反射波は、直流波である(もしくは、殆ど直流波である)と仮定すると、I信号とQ信号との各々の時間平均として把握される。そこで、静止物反射除去タスク33bは、記憶部32に記憶されているI信号およびQ信号のデータを読み込んで所定の分だけ蓄積するとともに蓄積したI信号とQ信号との各々の時間平均を計算し、計算した時間平均に相当する信号を、I信号とQ信号とのそれぞれから差し引くことによって静止物からの反射波に相当する成分を除去する(ステップS10)。 Assuming that the reflected wave from a stationary object is a DC wave (or is almost a DC wave), it is understood as the time average of each of the I signal and Q signal. Therefore, the stationary object reflection removal task 33b reads the data of the I signal and Q signal stored in the storage unit 32, accumulates a predetermined amount, and calculates the time average of each of the accumulated I signal and Q signal. Then, by subtracting a signal corresponding to the calculated time average from each of the I signal and Q signal, components corresponding to reflected waves from stationary objects are removed (step S10).

静止物からの反射波に相当する成分を除去するための、I信号およびQ信号のデータの蓄積の程度は、特定の時間長さあるいはデータ数に限定されるものではなく、静止物からの反射波に相当する成分を良好に除去することができるような適切な平均が計算され得る適当な時間長さやデータ数に適宜設定される。 The degree of accumulation of I signal and Q signal data in order to remove components corresponding to reflected waves from stationary objects is not limited to a specific time length or number of data; An appropriate length of time and number of data are appropriately set so that an appropriate average can be calculated so that components corresponding to waves can be satisfactorily removed.

静止物反射除去タスク33bは、静止物からの反射波に相当する成分の除去処理後のI信号とQ信号との組み合わせのデータ(「第1の除去処理後IQデータ」と呼ぶ)を記憶部32に記憶させる。 The stationary object reflection removal task 33b stores data of a combination of an I signal and a Q signal (referred to as "first post-removal IQ data") after processing to remove a component corresponding to a reflected wave from a stationary object in a storage unit. 32 to be memorized.

体動成分除去タスク33cは、第1の除去処理後IQデータを記憶部32から読み込み、I信号およびQ信号のそれぞれから、観測対象の生体の意識とは関係なく起きる低周期の体全体の揺れ(「体動」とも呼ばれる)に相当する成分を除去する。 The body motion component removal task 33c reads the IQ data after the first removal process from the storage unit 32, and extracts from each of the I signal and Q signal the low-period shaking of the entire body that occurs regardless of the consciousness of the living organism being observed. (also called "body movement").

体動に相当する成分を除去する仕法は、上述したとおり特定の手法や手順に限定されないものの、具体的には例えば、I信号のデータとQ信号のデータとのそれぞれを所定の時間長さで区切り、区切られた区間それぞれのデータのN次関数近似曲線(但し、Nは自然数)を推定することによって体動に相当する成分を表す回帰曲線を求め、求めた回帰曲線に相当する信号を、I信号のデータとQ信号のデータとのそれぞれから差し引く手法が用いられ得る。 The method for removing components corresponding to body movements is not limited to a specific method or procedure as described above, but specifically, for example, removing I signal data and Q signal data for a predetermined length of time, A regression curve representing a component corresponding to body movement is obtained by estimating an N-order function approximation curve (where N is a natural number) of data for each section and section, and a signal corresponding to the obtained regression curve is A method of subtracting from each of the data of the I signal and the data of the Q signal may be used.

体動成分除去タスク33cは、記憶部32に記憶されている第1の除去処理後IQデータを読み込み、読み込んだ第1の除去処理後IQデータから体動に相当する成分を除去する(ステップS11)。 The body movement component removal task 33c reads the IQ data after the first removal process stored in the storage unit 32, and removes the component corresponding to the body movement from the read IQ data after the first removal process (step S11 ).

体動成分除去タスク33cは、体動に相当する成分の除去処理後のI信号とQ信号との組み合わせのデータ(「第2の除去処理後IQデータ」と呼ぶ)を記憶部32に記憶させる。 The body motion component removal task 33c causes the storage unit 32 to store data of the combination of the I signal and the Q signal after the component corresponding to the body motion has been removed (referred to as "second post-removal IQ data"). .

FFT処理タスク33fは、高速フーリエ変換(FFT:Fast Fourier Transform の略)処理を用いてビート信号の周波数スペクトルを生成するためのタスク・プログラムである。 The FFT processing task 33f is a task program for generating a frequency spectrum of a beat signal using fast Fourier transform (FFT) processing.

FFT処理タスク33fは、第2の除去処理後IQデータについて、I信号を実部とするとともにQ信号を虚部とする複素数のビート信号として高速フーリエ変換処理を用いてビート信号の周波数スペクトルを生成する(ステップS12)。 The FFT processing task 33f generates a frequency spectrum of the beat signal using fast Fourier transform processing for the IQ data after the second removal process as a complex number beat signal with the I signal as the real part and the Q signal as the imaginary part. (Step S12).

そして、FFT処理タスク33fは、生成した周波数スペクトルを記憶部32に記憶させる。 Then, the FFT processing task 33f causes the storage unit 32 to store the generated frequency spectrum.

安静状態推定タスク33gは、生体が安静状態であるか否かを判断するためのタスク・プログラムである。安静状態推定タスク33gは、周波数領域の信号分布に注目して、生体が安静状態であるか否かを判断する。 The resting state estimation task 33g is a task program for determining whether the living body is in a resting state. The resting state estimation task 33g determines whether the living body is in a resting state by paying attention to the signal distribution in the frequency domain.

生体が不在のときの周波数スペクトルの分布ではすべての周波数領域(即ち、0Hz~∞Hz)においてスペクトル強度が概ね同じ強度で(言い換えると、大きく変動することなく)分布するのに対して(図4(A)参照)、生体が存在して安静状態にあるときの周波数スペクトルの分布では特定の周波数範囲においてスペクトル強度が他の周波数と比べて大きく変動する(同図(B)参照)。 In the frequency spectrum distribution when no living body is present, the spectral intensity is distributed at approximately the same intensity (in other words, without large fluctuations) in all frequency regions (i.e., 0 Hz to ∞ Hz) (Figure 4). In the frequency spectrum distribution when a living body is present and in a resting state, the spectral intensity fluctuates significantly in a specific frequency range compared to other frequencies (see (B) in the same figure).

そこで、安静状態推定タスク33gは、特定の周波数範囲のスペクトル強度の合計(「特定範囲強度合計値」と呼ぶ)と、前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計(「その他領域強度合計値」と呼ぶ)とを比較することにより、生体が安静状態であるか否かを判断する。 Therefore, the resting state estimation task 33g calculates the sum of the spectral intensities in a specific frequency range (referred to as "specific range intensity total value") and the sum of the spectral intensities in a predetermined area other than the specific frequency range (referred to as "other area intensity total value"). It is determined whether the living body is in a resting state by comparing the total strength value (called "total strength value").

生体が安静状態であるときにスペクトル強度が変動する周波数の範囲(「特定範囲」と呼ぶ)は、大凡0.25Hz~3.25Hzの範囲である。具体的には、呼吸は概ね3秒に1回とすると呼吸の動きに対応するスペクトル強度は0.3Hz付近に分布し、心拍は概ね1秒に1回とすると拍動に対応するスペクトル強度は1Hz付近に分布する。また、人の場合、起こり得る生体活動の周波数として、心拍が大凡200回/分(尚、痙攣を起こすような状態になる)が上限であると言われている。 The frequency range in which the spectral intensity varies when the living body is at rest (referred to as the "specific range") is approximately 0.25 Hz to 3.25 Hz. Specifically, if breathing is approximately once every 3 seconds, the spectral intensity corresponding to the movement of breathing will be distributed around 0.3Hz, and if the heartbeat is approximately once every 1 second, the spectral intensity corresponding to the pulsation will be distributed around 0.3Hz. Distributed around 1Hz. Furthermore, in the case of humans, it is said that the upper limit of the frequency of biological activity that can occur is approximately 200 heartbeats/minute (which can cause convulsions).

上記を踏まえ、例えば、特定範囲が0Hz~5Hz未満程度に設定されるとともに特定範囲以外の所定の領域(「その他領域」と呼ぶ)が5Hz~20Hz程度に設定されたり、特定範囲が0.25Hz~3.25Hz未満程度に設定されるとともにその他領域が3.25Hz~10Hz程度に設定されたりする。 Based on the above, for example, the specific range may be set to about 0Hz to less than 5Hz, and a predetermined area other than the specific range (referred to as "other area") may be set to about 5Hz to 20Hz, or the specific range may be set to about 0.25Hz. 3.25 Hz to less than 3.25 Hz, and other regions are set to about 3.25 Hz to 10 Hz.

安静状態推定タスク33gは、記憶部32に記憶されている周波数スペクトルを読み込み、特定範囲強度合計値とその他領域強度合計値とをそれぞれ計算する。安静状態推定タスク33gは、続いて、特定範囲強度合計値がその他領域強度合計値よりも大きいか否かを判断する(ステップS13)。 The resting state estimation task 33g reads the frequency spectrum stored in the storage unit 32, and calculates the specific range intensity total value and the other area intensity total value, respectively. The resting state estimation task 33g then determines whether the specific range total strength value is larger than the other area strength total value (step S13).

特定範囲強度合計値がその他領域強度合計値よりも大きい場合には(ステップS13:Yes)、安静状態推定タスク33gは、生体が安静状態であると判断する(ステップS14)。そして、メインタスク33は、処理手順をステップS19へと進める。 If the specific range intensity total value is larger than the other area intensity total value (step S13: Yes), the resting state estimation task 33g determines that the living body is in a resting state (step S14). Then, the main task 33 advances the processing procedure to step S19.

一方、特定範囲強度合計値がその他領域強度合計値以下である場合には(ステップS13:No)、安静状態推定タスク33gは、処理手順をステップS15へと進める。 On the other hand, if the specific range intensity total value is less than or equal to the other area intensity total value (step S13: No), the resting state estimation task 33g advances the processing procedure to step S15.

なお、安静状態推定タスク33gは、特定範囲のスペクトル強度の平均値とその他領域のスペクトル強度の平均値とを比較して、特定範囲のスペクトル強度の平均値がその他領域のスペクトル強度の平均値よりも大きい場合に生体が安静状態であると判断するようにしてもよい。 Note that the resting state estimation task 33g compares the average value of the spectral intensity in the specific range with the average value of the spectral intensity in other areas, and determines whether the average value of the spectral intensity in the specific range is higher than the average value of the spectral intensity in the other areas. The living body may be determined to be in a resting state when the amount is also large.

ステップS15以降の処理では、ステップS11の処理において記憶部32に記憶される第2の除去処理後IQデータ(図2においては「減算信号」)が用いられる。 In the processes after step S15, the second removed IQ data ("subtraction signal" in FIG. 2) stored in the storage unit 32 in the process in step S11 is used.

不在状態推定タスク33hは、生体が不在であるか否かを判断するためのタスク・プログラムである。不在状態推定タスク33hは、ホワイトノイズの振幅の大きさ別の度数分布が正規分布に従うことに注目して、生体が不在であるか否かを判断する(ステップS16)。 The absent state estimation task 33h is a task program for determining whether or not a living body is absent. The absent state estimation task 33h determines whether or not the living body is absent, paying attention to the fact that the frequency distribution according to the amplitude of the white noise follows a normal distribution (step S16).

不在状態推定タスク33hは、第2の除去処理後IQデータとしてのビート信号(具体的には、I信号,Q信号)の振幅の大きさ別の度数分布について正規性の検定を行う。不在状態推定タスク33hは、具体的には、記憶部32に記憶されている第2の除去処理後IQデータとしてのI信号およびQ信号を読み込んで所定の分だけ蓄積するとともに、蓄積したI信号とQ信号とのそれぞれについて振幅の大きさ別の度数を計算する(図5参照)。図5について、上の図の1本の線はI信号の振幅の時間変化を表し、また、下の図はI信号についての振幅の大きさ別の度数分布を表し、特に、下の左の図は度数分布形が正規分布である場合の例であり、下の右の図は度数分布形が正規分布でない場合の例である。 The absent state estimation task 33h tests the normality of the frequency distribution according to the amplitude of the beat signal (specifically, the I signal and the Q signal) as the IQ data after the second removal process. Specifically, the absent state estimation task 33h reads the I signal and Q signal as the second post-removal IQ data stored in the storage unit 32, accumulates a predetermined amount, and stores the accumulated I signal. The frequency of each amplitude is calculated for each of the signal and the Q signal (see FIG. 5). Regarding Figure 5, one line in the upper figure represents the temporal change in the amplitude of the I signal, and the lower figure represents the frequency distribution of the I signal by amplitude size. The figure shows an example where the frequency distribution type is a normal distribution, and the lower right figure is an example where the frequency distribution type is not a normal distribution.

振幅の大きさ別の度数を計算するための、I信号やQ信号のデータの蓄積の程度は、特定の時間長さあるいはデータ数に限定されるものではなく、振幅の大きさ別の度数分布が正規分布に従うか否かを良好に判定することができるような適切な振幅の大きさ別の度数が計算され得る適当な時間長さやデータ数に適宜設定される。また、振幅の大きさについて所定のピッチで区切られてランクが設定され、振幅の大きさのランク別に度数が計算されるようにしてもよい。 The degree of accumulation of I signal and Q signal data for calculating the frequency for each amplitude size is not limited to a specific time length or number of data, but is based on the frequency distribution for each amplitude size. The appropriate time length and number of data are appropriately set so that an appropriate frequency for each amplitude size can be calculated so that it can be determined satisfactorily whether or not the amplitude follows a normal distribution. Alternatively, ranks may be set for the magnitude of the amplitude divided by a predetermined pitch, and the frequency may be calculated for each rank of the magnitude of the amplitude.

I信号やQ信号の振幅の大きさ別の度数分布形が正規分布である場合には、生体は不在であると判断される。なお、I信号やQ信号の振幅の大きさ別の度数分布形が正規分布であるか否かの検定は、度数分布形の正規分布への適合度が所定の適合度以上であるか否かによって行われるようにしてもよい。また、振幅の大きさ別の度数分布形が正規分布であるか否かの検定は、I信号とQ信号とのうちのどちらか一方のみについて行われるようにしてもよく、或いは、I信号とQ信号との両方について行われるようにしてもよい。 If the frequency distribution shape of each amplitude of the I signal or Q signal is a normal distribution, it is determined that a living body is absent. Note that the test of whether the frequency distribution shape of each amplitude of the I signal or Q signal is a normal distribution is to check whether the degree of conformity of the frequency distribution shape to the normal distribution is greater than or equal to a predetermined degree of conformity. It may also be done by Further, the test as to whether or not the frequency distribution shape of each amplitude is a normal distribution may be performed on only one of the I signal and the Q signal, or may be performed on only one of the I signal and the Q signal. The process may be performed for both the Q signal and the Q signal.

ビート信号の振幅の大きさ別の度数分布形が正規分布であると判定される場合には(ステップS16:Yes)、不在状態推定タスク33hは、生体は不在である(不在状態である)と判断する(ステップS17)。そして、メインタスク33は、処理手順をステップS19へと進める。 If it is determined that the frequency distribution shape according to the magnitude of the amplitude of the beat signal is a normal distribution (step S16: Yes), the absent state estimation task 33h determines that the living body is absent (in an absent state). A judgment is made (step S17). Then, the main task 33 advances the processing procedure to step S19.

以上のように、生体状態推定装置1は、はじめに、移動状態に対応する大きな動作を生体がしているか否かを判断し(ステップS7)、次に、大きな動作をしていないと判断された場合に安静状態に対応する小さな(或いは、微小な)動作を生体がしているか否かを判断し(ステップS13)、さらに、小さな動作をしていないと判断された場合に不在状態であるか否かを判断する(ステップS16)ようにしている。このように、移動状態、安静状態、そして不在状態というように動作の大きい順に段階的に状態を区分して特定することにより、無駄なく効率的に且つ精度高く、生体の状態が推定される。 As described above, the living body state estimation device 1 first determines whether the living body is making a large movement corresponding to the moving state (step S7), and then it is determined that the living body is not making a large movement. In this case, it is determined whether the living body is making small (or minute) movements corresponding to a resting state (step S13), and further, if it is determined that the living body is not making small movements, it is determined whether the living body is in an absent state or not. It is determined whether or not (step S16). In this way, the state of the living body can be estimated efficiently and accurately without waste by classifying and specifying the states step by step in descending order of movement, such as the moving state, the resting state, and the absent state.

一方、ビート信号の振幅の大きさ別の度数分布形が正規分布でないと判定される場合には(ステップS16:No)、不在状態推定タスク33hは、混信の状態であると判断する(ステップS18)。そして、メインタスク33は、生体状態の推定結果を保留して、処理手順をステップS19へと進める。 On the other hand, if it is determined that the frequency distribution shape according to the magnitude of the amplitude of the beat signal is not a normal distribution (step S16: No), the absent state estimation task 33h determines that there is interference (step S18). ). Then, the main task 33 suspends the estimation result of the biological state and advances the processing procedure to step S19.

危険状態検知タスク33iは、生体が危険状態にあることを検知するためのタスク・プログラムである。危険状態検知タスク33iは、時系列における直前の生体の状態と新しく推定された生体の状態との組み合わせに基づいて、言い換えると、時系列における生体の状態の変化のパターンに基づいて、生体が危険状態にあるか否かを判定する(ステップS19)。 The dangerous state detection task 33i is a task program for detecting that a living body is in a dangerous state. The dangerous state detection task 33i detects whether the living body is in danger based on the combination of the immediately previous state of the living body in the time series and the newly estimated state of the living body, in other words, based on the pattern of changes in the state of the living body in the time series. It is determined whether it is in the state (step S19).

危険状態検知タスク33iは、例えば、生体の状態の推定結果の時系列での推移において、安静状態から、移動状態を経ることなく、不在状態へと推移した場合には、バイタルがある状態からバイタルが無い状態へと直接推移したので、言い換えると、呼吸をしている状態から呼吸をしていない状態へと直接推移したので、生体が危険状態にあると判断する。この場合、危険状態検知タスク33iは、例えば、警報音を発したり、予め登録されているアドレスへとメールを送信したりして、生体が危険状態にあることを知らせるようにしてもよい。 For example, in the time-series transition of the estimated results of the state of the living body, if the state of the living body changes from a resting state to an absent state without passing through a moving state, the dangerous state detection task 33i changes the vital state from a state with vitals to an absent state without passing through a moving state. In other words, there was a direct transition from a breathing state to a non-breathing state, so it is determined that the living body is in a dangerous state. In this case, the dangerous state detection task 33i may notify that the living body is in a dangerous state by emitting an alarm sound or sending an email to a pre-registered address, for example.

なお、生体の状態の推定結果の時系列での推移において、不在状態から、移動状態を経ることなく、安静状態へと推移した場合には、観測対象の所定の領域において例えば遠隔操作や自律制御によって電気機器等が作動したり、地震が発生したりするなどの外乱が発生したと判別して、生体の状態についてではないその他の異常事象が発生したと判断するようにしてもよい。 In addition, in the time-series transition of the estimation results of the state of the living organism, if the state changes from the absent state to the resting state without passing through the moving state, for example, remote control or autonomous control may be performed in the predetermined area of the observation target. It may be determined that a disturbance such as the activation of electrical equipment or the occurrence of an earthquake has occurred, and it may be determined that another abnormal event that is not related to the state of the living body has occurred.

そして、メインタスク33は、処理手順をステップS2へと戻す。なお、上述の通り、送受信部2によるレーダースキャンは所定の時間間隔で連続的に繰り返し行われ、受信部22から、デジタル信号化されたI信号とQ信号との組み合わせが、前記所定の時間間隔で連続して継続的に出力される(ステップS1)。 The main task 33 then returns the processing procedure to step S2. As mentioned above, the radar scan by the transmitting/receiving section 2 is continuously and repeatedly performed at predetermined time intervals, and the combination of the digitalized I signal and Q signal is sent from the receiving section 22 at the predetermined time interval. It is output continuously (step S1).

上記のような生体状態推定装置1によれば、例えば壁や床などの静止物からの反射波に相当する成分を自律的に除去するようにしているので、使用環境ごとに閾値を設定することなく、言い換えると、使用環境に影響を受けることなく、また、使用環境に依存することなく、生体の不在状態、安静状態、および移動状態を高精度に判別することが可能となる。上記のような生体状態推定装置1によれば、また、移動状態、安静状態、そして不在状態というように動作の大きい順に段階的に状態を区分して特定するようにしているので、無駄なく効率的に且つ精度高く、生体の状態を推定することが可能となる。 According to the biological state estimation device 1 described above, components corresponding to reflected waves from stationary objects such as walls and floors are autonomously removed, so a threshold value can be set for each usage environment. In other words, it is possible to determine with high precision whether the living body is absent, at rest, or moving, without being affected by or depending on the usage environment. According to the biological state estimating device 1 described above, the state is classified and specified in stages in descending order of movement, such as a moving state, a resting state, and an absent state, so that efficiency can be achieved without waste. It becomes possible to estimate the condition of the living body in a highly accurate and accurate manner.

以上、この発明の実施の形態について説明したが、具体的な構成は、上記の実施の形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計の変更等があっても、この発明に含まれる。例えば、上記の実施の形態では呼吸の速度や特定範囲を確定する閾値の設定について人を観測対象とする場合を中心に説明したが、この発明における観測対象の生体は人に限定されるものではなく、種々の動物を観測対象としてこの発明が適用されるようにしてもよい。 Although the embodiments of the present invention have been described above, the specific configuration is not limited to the above embodiments, and even if there are changes in the design within the scope of the gist of the present invention. Included in invention. For example, in the above embodiment, the setting of the threshold value for determining the breathing rate and specific range was mainly explained in the case where the observation target is a human, but the living organism to be observed in this invention is not limited to humans. However, the present invention may be applied to various animals as observation objects.

また、上記の実施の形態ではメインタスク33が危険状態検知タスク33iを含むようにしているが、危険状態検知タスク33iはこの発明において必須の構成ではない。すなわち、移動状態、安静状態、および不在状態のうちのいずれであるかを推定する構成が、この発明における基本の構成である。 Further, in the above embodiment, the main task 33 includes the dangerous state detection task 33i, but the dangerous state detection task 33i is not an essential component in the present invention. That is, the basic configuration of the present invention is a configuration for estimating which of the moving state, the resting state, and the absent state.

なお、所定の領域(即ち、観測対象の生体が居る空間)に対して複数の送受信部を設置し、得られる受信データを周知の技術を用いて物標ごとに分離して、複数の生体の状態を同時に推定するようにしてもよい。 In addition, multiple transmitter/receivers are installed in a predetermined area (i.e., the space where the living organisms to be observed are present), and the received data obtained is separated for each target object using well-known technology. The states may be estimated simultaneously.

この発明は、人や種々の動物のバイタル情報に基づいて生体の状態を推定することができるので、例えば、在宅の高齢者の遠隔での見守り、療養中の患者の容体観察、および飼育動物や保護下の動物の状態監視などの分野に適用され得る。また、住宅や施設への侵入監視の分野にも適用され得る。 This invention can estimate the condition of a living body based on the vital information of people and various animals, so it can be used, for example, to remotely monitor elderly people at home, observe the condition of patients undergoing medical treatment, and use domestic animals and other animals. It can be applied in fields such as monitoring the condition of animals under protection. It can also be applied to the field of intrusion monitoring into homes and facilities.

1 生体状態推定装置
2 送受信部
21 送信部
211 送信用アンテナ
22 受信部
221 受信用アンテナ
3 推定装置
31 入力受部
32 記憶部
33 メインタスク
33a 位相計算タスク
33b 静止物反射除去タスク
33c 体動成分除去タスク
33d 速度算出タスク
33e 移動状態推定タスク
33f FFT処理タスク
33g 安静状態推定タスク
33h 不在状態推定タスク
33i 危険状態検知タスク
34 中央処理部
1 biological state estimation device 2 transmitting/receiving section 21 transmitting section 211 transmitting antenna 22 receiving section 221 receiving antenna 3 estimating device 31 input receiving section 32 storage section 33 main task 33a phase calculation task 33b stationary object reflection removal task 33c body motion component removal Tasks 33d Speed calculation task 33e Moving state estimation task 33f FFT processing task 33g Resting state estimation task 33h Absence state estimation task 33i Dangerous state detection task 34 Central processing unit

Claims (3)

送信電波を放射して前記送信電波の反射波を受信するとともに前記送信電波と前記反射波とを混合してビート信号を生成する送受信部と、
前記送受信部から出力される前記ビート信号の位相を計算する位相計算部と、
前記送受信部から出力される前記ビート信号のデータと前記位相計算部から出力される前記位相のデータとのそれぞれから静止物からの反射波に相当する成分を除去する静止物反射除去部と、
前記送受信部から出力される前記ビート信号のデータと前記位相計算部から出力される前記位相のデータとのそれぞれから体動に相当する成分を除去する体動成分除去部と、
前記2つの成分の除去後の前記位相のデータを用いて前記送受信部と物標としての生体との相対速度を算出する速度算出部と、
前記相対速度が前記生体の呼吸の速度よりも大きい場合に前記生体が移動状態である判断する移動状態推定部と、
前記相対速度が前記生体の前記呼吸の速度以下であり、且つ、前記2つの成分の除去後の前記ビート信号のデータを高速フーリエ変換して得られる周波数スペクトルの分布のうちの特定の周波数範囲のスペクトル強度の合計前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計よりも大きい場合に前記生体が安静状態である判断する安静状態推定部と、
前記特定の周波数範囲のスペクトル強度の合計が前記特定の周波数範囲以外の所定の領域のスペクトル強度の合計以下である場合に前記2つの成分の除去後の前記ビート信号の振幅の大きさ別の度数分布形が正規分布であるか否かに基づいて前記生体が不在状態であるか否かを判断する不在状態推定部と、を有
前記特定の周波数範囲が、0Hz~5Hz未満に設定される、
ことを特徴とする生体状態推定装置。
a transmitting/receiving unit that emits a transmitted radio wave, receives a reflected wave of the transmitted radio wave, and mixes the transmitted radio wave and the reflected wave to generate a beat signal;
a phase calculation unit that calculates the phase of the beat signal output from the transmission and reception unit;
a stationary object reflection removing section that removes a component corresponding to a reflected wave from a stationary object from each of the beat signal data output from the transmitting and receiving section and the phase data output from the phase calculation section;
a body movement component removing unit that removes a component corresponding to body movement from each of the beat signal data output from the transmitting and receiving unit and the phase data output from the phase calculation unit;
a speed calculating unit that calculates a relative speed between the transmitting/receiving unit and a living body as a target using the phase data after removing the two components;
a moving state estimator that determines that the living body is in a moving state when the relative velocity is greater than the respiration rate of the living body ;
The relative velocity is equal to or lower than the respiration rate of the living body, and the frequency range is within a specific frequency range of the frequency spectrum distribution obtained by fast Fourier transforming the beat signal data after the two components have been removed. a resting state estimation unit that determines that the living body is in a resting state when the sum of spectral intensities is greater than the sum of spectral intensities in a predetermined region other than the specific frequency range;
If the sum of the spectral intensities in the specific frequency range is less than or equal to the sum of the spectral intensities in a predetermined area other than the specific frequency range, the frequency of each amplitude of the beat signal after removal of the two components. an absent state estimation unit that determines whether the living body is in an absent state based on whether the distribution shape is a normal distribution;
the specific frequency range is set from 0Hz to less than 5Hz;
A biological state estimation device characterized by:
記生体の前記呼吸の速度が、1cm~10cm/秒に設定される、
ことを特徴とする請求項1に記載の生体状態推定装置。
The respiration rate of the living body is set to 1 cm to 10 cm/sec.
The biological condition estimation device according to claim 1, characterized in that:
前記生体の状態の推定結果が、前記安静状態から、前記移動状態を経ることなく、前記不在状態へと推移した場合に、前記生体が危険状態にあると判断する危険状態検知部をさらに有する、
ことを特徴とする請求項1または2に記載の生体状態推定装置。
further comprising a dangerous state detection unit that determines that the living body is in a dangerous state when the estimation result of the state of the living body transitions from the resting state to the absent state without passing through the moving state;
The biological condition estimation device according to claim 1 or 2, characterized in that:
JP2020083450A 2020-05-11 2020-05-11 Biological state estimation device Active JP7451040B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020083450A JP7451040B2 (en) 2020-05-11 2020-05-11 Biological state estimation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020083450A JP7451040B2 (en) 2020-05-11 2020-05-11 Biological state estimation device

Publications (2)

Publication Number Publication Date
JP2021177830A JP2021177830A (en) 2021-11-18
JP7451040B2 true JP7451040B2 (en) 2024-03-18

Family

ID=78509908

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020083450A Active JP7451040B2 (en) 2020-05-11 2020-05-11 Biological state estimation device

Country Status (1)

Country Link
JP (1) JP7451040B2 (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012005745A (en) 2010-06-28 2012-01-12 Tateyama System Laboratory Ltd Bed nursing system
JP2017000484A (en) 2015-06-11 2017-01-05 富士通株式会社 Non-contact activity amount sensor and air conditioner
WO2017183528A1 (en) 2016-04-19 2017-10-26 コニカミノルタ株式会社 Device, method, and system for monitoring of monitor subject
WO2018216363A1 (en) 2017-05-25 2018-11-29 コニカミノルタ株式会社 Care support system and radio wave control method
JP2019061643A (en) 2017-09-27 2019-04-18 シャープ株式会社 Monitoring device
JP2019155071A (en) 2018-03-08 2019-09-19 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカPanasonic Intellectual Property Corporation of America Event prediction system, sensor signal processing system, event prediction method, and program

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012005745A (en) 2010-06-28 2012-01-12 Tateyama System Laboratory Ltd Bed nursing system
JP2017000484A (en) 2015-06-11 2017-01-05 富士通株式会社 Non-contact activity amount sensor and air conditioner
WO2017183528A1 (en) 2016-04-19 2017-10-26 コニカミノルタ株式会社 Device, method, and system for monitoring of monitor subject
WO2018216363A1 (en) 2017-05-25 2018-11-29 コニカミノルタ株式会社 Care support system and radio wave control method
JP2019061643A (en) 2017-09-27 2019-04-18 シャープ株式会社 Monitoring device
JP2019155071A (en) 2018-03-08 2019-09-19 パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカPanasonic Intellectual Property Corporation of America Event prediction system, sensor signal processing system, event prediction method, and program

Also Published As

Publication number Publication date
JP2021177830A (en) 2021-11-18

Similar Documents

Publication Publication Date Title
CA2891839C (en) Remote sensing of human breathing at a distance
US8742935B2 (en) Radar based systems and methods for detecting a fallen person
Forouzanfar et al. Event recognition for contactless activity monitoring using phase-modulated continuous wave radar
JP6874679B2 (en) Monitoring device
JP2017153850A (en) Living body monitoring system and living body monitoring method
WO2016136400A1 (en) Living-body monitoring device, living-body monitoring method, and living-body monitoring system
JP2010066877A (en) Resident monitoring system and resident monitoring method
KR102323343B1 (en) Patient monitering system using rada
CN110461230B (en) Information processing apparatus and detection apparatus
Takabatake et al. FMCW radar-based anomaly detection in toilet by supervised machine learning classifier
JP6856071B2 (en) Respiratory rate display device and respiratory rate display method
JP7451040B2 (en) Biological state estimation device
JP7072403B2 (en) Vital sensor
JP2009060989A (en) Respiration monitoring method and apparatus
US11399723B2 (en) Non-contact method of physiological characteristic detection
WO2020071374A1 (en) Condition monitoring device and condition monitoring method
KR102103260B1 (en) A accupancy sensing apparatus and method using Ultra-Wideband Radar
JP2019092726A (en) Biological information detection device and control method of biological information detection device
KR102060102B1 (en) Method for managing aging and vulnerability behaviors using ultra wide band
Kocur et al. Estimation of breathing frequency and heart rate by biometric UWB radar
JP6766585B2 (en) Observed person monitoring device and method
JP7294845B2 (en) Biological information detector
WO2021187308A1 (en) Living body detection device, living body detection method, and program
KR20160139651A (en) Monitering System for the old and the weak
JP5845372B1 (en) Health condition monitoring system

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230428

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231018

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231114

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240112

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240305

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240305

R150 Certificate of patent or registration of utility model

Ref document number: 7451040

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150