JP6299172B2 - Information processing apparatus, information processing method, and program - Google Patents

Information processing apparatus, information processing method, and program Download PDF

Info

Publication number
JP6299172B2
JP6299172B2 JP2013236921A JP2013236921A JP6299172B2 JP 6299172 B2 JP6299172 B2 JP 6299172B2 JP 2013236921 A JP2013236921 A JP 2013236921A JP 2013236921 A JP2013236921 A JP 2013236921A JP 6299172 B2 JP6299172 B2 JP 6299172B2
Authority
JP
Japan
Prior art keywords
time
period
unit
autoregressive
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
JP2013236921A
Other languages
Japanese (ja)
Other versions
JP2015096831A (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.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry 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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP2013236921A priority Critical patent/JP6299172B2/en
Publication of JP2015096831A publication Critical patent/JP2015096831A/en
Application granted granted Critical
Publication of JP6299172B2 publication Critical patent/JP6299172B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、情報処理装置、情報処理方法及びプログラムに関する。   The present invention relates to an information processing apparatus, an information processing method, and a program.

近年、部屋の見守りや対象物体の状態認識等のために、振動センサや加速度センサ、ドップラーセンサ等のセンサから得られる信号(センサ信号)の定常性を判定する技術が開発されている。定常性の判定は、信号の周期の時間的な変化の度合いを定量的に評価することにより行われる。センサ信号の定常性を判定することにより、信号の性質が明らかになり、センサが観測する対象物体や対象空間の状態や状態の変化を明らかにすることができる。   In recent years, techniques for determining the continuity of signals (sensor signals) obtained from sensors such as a vibration sensor, an acceleration sensor, and a Doppler sensor have been developed for watching a room and recognizing the state of a target object. The determination of continuity is performed by quantitatively evaluating the degree of temporal change in the signal cycle. By determining the continuity of the sensor signal, the nature of the signal is clarified, and the state of the target object and the target space observed by the sensor and the change in the state can be clarified.

このような技術としては例えば、下記特許文献1では、生体信号の超低周波帯域のゆらぎに対してゼロクロス値を適用して、時間窓毎に得られる周波数の平均値の時系列変化を周波数変動時系列波形として求める技術が開示されている。   As such a technique, for example, in Patent Document 1 below, a zero cross value is applied to fluctuations in a very low frequency band of a biological signal, and a time series change of an average value of frequencies obtained for each time window is changed in frequency. A technique for obtaining a time series waveform is disclosed.

また、下記特許文献2では、家庭内の生活音を音センサで収集し、短時間フーリエ変換で周波数ごとの成分を抽出して、利用者の生活パターンを把握する技術が開示されている。   Patent Document 2 below discloses a technique for collecting a living sound in a home with a sound sensor, extracting a component for each frequency by a short-time Fourier transform, and grasping a user's life pattern.

また、下記非特許文献1では、下記数式1で表される時変係数自己回帰モデルにおいて、ある時間ステップごとに得られるデータを用いて、現在から1個前からN個前(Nは自然数)の自己回帰係数を、カルマンフィルタを用いて求める技術が開示されている。   Further, in the following non-patent document 1, in the time-varying coefficient autoregressive model expressed by the following mathematical formula 1, using data obtained for each time step, the current one to the previous N (N is a natural number) The technique which calculates | requires the autoregressive coefficient of this using a Kalman filter is disclosed.

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

特開2013−111103号公報JP2013-111103A 特開2011−237865号公報JP2011-237865A

北川源四郎.時系列解析入門.岩波書店,2005.Genshiro Kitagawa. Introduction to time series analysis. Iwanami Shoten, 2005.

しかし、上記文献に開示された技術では、センサから出力された信号の定常性をリアルタイムに判定する場合に、精度や計算コストの面で問題があった。   However, the technique disclosed in the above document has a problem in terms of accuracy and calculation cost when determining the continuity of the signal output from the sensor in real time.

例えば、上記特許文献1に開示された技術で用いられるゼロクロス値は、大きな振幅の周波数が支配的になってしまうため、信号対雑音比(SN比)が低く、雑音の振幅が相対的に大きい場合に所望信号の正確な周波数を検出することが困難であった。   For example, the zero-cross value used in the technique disclosed in Patent Document 1 has a large amplitude frequency, so the signal-to-noise ratio (S / N ratio) is low and the noise amplitude is relatively large. In some cases, it was difficult to detect the exact frequency of the desired signal.

上記特許文献2に開示された技術では、フーリエ変換を用いているため、時間窓ごとのデータを使って演算する際の計算量が多い。このため、計算能力の低い端末では、逐次的に周波数成分を抽出するといったリアルタイムな動作を行うことが困難である。また、データ数が少ない場合は十分な周波数分解能が得られないという問題がある。さらに、周波数領域で所望信号とノイズの周波数とが重なる場合、ノイズと所望信号の分離が困難である。   Since the technique disclosed in Patent Document 2 uses Fourier transform, the amount of calculation when performing calculations using data for each time window is large. For this reason, it is difficult for a terminal with low calculation capability to perform a real-time operation of sequentially extracting frequency components. In addition, when the number of data is small, there is a problem that sufficient frequency resolution cannot be obtained. Furthermore, when the desired signal and the noise frequency overlap in the frequency domain, it is difficult to separate the noise and the desired signal.

上記非特許文献1に開示された技術では、短期間の時間ステップの自己回帰係数を用いるときは計算コストが小さいものの、局所的な信号の特徴のみが自己回帰係数に反映されるため、一周期あたりのサンプル数が多い定常信号を精度よく予測することが困難である。また、多サンプル数の観測値を用いて自己回帰係数を求めることは計算コストが大きく、計算能力が低い端末での動作は困難である。   In the technique disclosed in Non-Patent Document 1, when the autoregressive coefficient of a short time step is used, the calculation cost is small, but only local signal features are reflected in the autoregressive coefficient. It is difficult to accurately predict a stationary signal having a large number of samples per round. In addition, obtaining an autoregressive coefficient using observation values of a large number of samples is expensive in calculation and difficult to operate on a terminal having low calculation capability.

そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、センサから出力される定常信号のリアルタイムな予測をより少ない計算量で精度良く行うことが可能な、新規かつ改良された情報処理装置、情報処理方法及びプログラムを提供することにある。   Accordingly, the present invention has been made in view of the above problems, and an object of the present invention is to perform real-time prediction of a steady signal output from a sensor with a small amount of calculation with high accuracy. Another object of the present invention is to provide a new and improved information processing apparatus, information processing method and program.

上記課題を解決するために、本発明のある観点によれば、センサから出力された時系列信号の周期を抽出する抽出部と、前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、を備える情報処理装置が提供される。   In order to solve the above problems, according to an aspect of the present invention, an extraction unit that extracts a period of a time-series signal output from a sensor, and an autoregressive coefficient corresponding to the period extracted by the extraction unit are provided. Provided is an information processing apparatus comprising: a coefficient calculation unit that sequentially calculates; and a prediction unit that predicts the time-series signal using an autoregressive model equation using the autoregressive coefficient calculated by the coefficient calculating unit. Is done.

前記抽出部は、逐次的に前記周期を抽出してもよい。   The extraction unit may sequentially extract the period.

前記抽出部は、逐次自己共分散関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出してもよい。   The extraction unit may extract a time step that takes at least one of a maximum value and a maximum value of a sequential autocovariance function as the period.

前記抽出部は、逐次自己相関関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出してもよい。   The extraction unit may extract a time step taking at least one of a maximum value and a maximum value of a sequential autocorrelation function as the period.

前記抽出部は、所定区間の前記時系列信号において前記最大値または極大値の少なくともいずれをとる時間ステップのうち、閾値を超える時間ステップを、前記周期として抽出してもよい。   The extraction unit may extract, as the period, a time step that exceeds a threshold among time steps that take at least either the maximum value or the maximum value in the time-series signal in a predetermined section.

前記係数算出部は、単位数量の時間ステップ分の前記時系列信号が出力される都度、前記自己回帰係数を算出してもよい。   The coefficient calculation unit may calculate the autoregressive coefficient each time the time series signal corresponding to a unit quantity time step is output.

前記係数算出部は、カルマンフィルタにより前記自己回帰係数を算出してもよい。   The coefficient calculation unit may calculate the autoregressive coefficient using a Kalman filter.

前記係数算出部は、前記周期に対応する自己回帰係数に加えて、任意の時間ステップの自己回帰係数を逐次的に算出してもよい。   The coefficient calculation unit may sequentially calculate an autoregressive coefficient at an arbitrary time step in addition to the autoregressive coefficient corresponding to the period.

前記情報処理装置は、前記予測部により予測された前記時系列信号の予測値と前記センサから出力された観測値との比較結果に基づいて前記時系列信号の定常性を判定する判定部をさらに備えてもよい。   The information processing apparatus further includes a determination unit that determines the continuity of the time series signal based on a comparison result between the predicted value of the time series signal predicted by the prediction unit and the observation value output from the sensor. You may prepare.

前記判定部は、前記予測値と前記観測値との差分が閾値以上である場合に非定常であると判定し、前記閾値未満である場合に定常であると判定してもよい。   The determination unit may determine that the difference between the predicted value and the observed value is greater than or equal to a threshold value, and may determine that the difference is less than the threshold value.

前記センサは、ドップラーセンサ、振動センサ、加速度センサの少なくともいずれかであってもよい。   The sensor may be at least one of a Doppler sensor, a vibration sensor, and an acceleration sensor.

また、上記課題を解決するために、本発明の別の観点によれば、センサから出力された時系列信号の周期を抽出するステップと、抽出された前記周期に対応する自己回帰係数を逐次的に算出するステップと、算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行うステップと、を備える情報処理方法が提供される。   In order to solve the above problem, according to another aspect of the present invention, a step of extracting a period of a time-series signal output from a sensor and an autoregressive coefficient corresponding to the extracted period are sequentially obtained. There is provided an information processing method comprising the steps of: calculating the time series signal by using an equation of an autoregressive model using the calculated autoregressive coefficient.

また、上記課題を解決するために、本発明の別の観点によれば、コンピュータを、センサから出力された時系列信号の周期を抽出する抽出部と、前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、として機能させるためのプログラムが提供される。   In order to solve the above problem, according to another aspect of the present invention, a computer is provided with an extraction unit that extracts a period of a time-series signal output from a sensor, and the period extracted by the extraction unit. A coefficient calculation unit that sequentially calculates a corresponding autoregressive coefficient, and a prediction unit that performs prediction of the time series signal using an equation of an autoregressive model using the autoregressive coefficient calculated by the coefficient calculation unit, and A program for functioning is provided.

以上説明したように本発明によれば、センサから出力される定常信号のリアルタイムな予測をより少ない計算量で精度良く行うことが可能である。   As described above, according to the present invention, real-time prediction of a steady signal output from a sensor can be accurately performed with a smaller amount of calculation.

本発明の一実施形態に係る定常性判定装置の構成の一例を示すブロック図である。It is a block diagram which shows an example of a structure of the stationarity determination apparatus which concerns on one Embodiment of this invention. 本発明の一実施形態に係る定常性判定装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the stationarity determination apparatus which concerns on one Embodiment of this invention. 比較例に係る定常性判定装置による予測誤差を示す図である。It is a figure which shows the prediction error by the stationarity determination apparatus which concerns on a comparative example. 本発明の一実施形態に係る定常性判定装置による予測誤差を示す図である。It is a figure which shows the prediction error by the stationarity determination apparatus which concerns on one Embodiment of this invention.

以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。   Exemplary embodiments of the present invention will be described below in detail with reference to the accompanying drawings. In addition, in this specification and drawing, about the component which has the substantially same function structure, duplication description is abbreviate | omitted by attaching | subjecting the same code | symbol.

<1.定常性判定装置の概要>
まず、本発明の一実施形態に係る定常性判定装置の概要を説明する。本実施形態に係る定常性判定装置は、センサ信号の周期の定常性を、リアルタイムに判定する。本実施形態に係る定常性判定装置は、振動センサや加速度センサ、ドップラーセンサ等から得られるセンサ信号の定常性を解析することで、対象とする物体の動作や空間の状況についての変化点を発見することができる。上述したように、上記文献に記載された、ゼロクロス値やフーリエ変換、時変係数自己回帰モデルを用いた技術では、精度や計算コストの面で問題があった。
<1. Overview of stationarity determination device>
First, an outline of a continuity determination device according to an embodiment of the present invention will be described. The stationarity determination device according to the present embodiment determines the stationarity of the period of the sensor signal in real time. The stationarity determination device according to this embodiment finds a change point about the motion of the target object and the state of space by analyzing the stationarity of sensor signals obtained from vibration sensors, acceleration sensors, Doppler sensors, and the like. can do. As described above, the technique using the zero cross value, the Fourier transform, and the time-varying coefficient autoregressive model described in the above document has problems in terms of accuracy and calculation cost.

そこで、本実施形態に係る定常性判定装置は、上記非特許文献1に開示された技術を改良して、精度よく且つ計算量を削減しつつリアルタイムでセンサ信号の定常性を判定することを可能にした。具体的には、定常性判定装置は、逐次自己共分散関数によってセンサ信号の周期を求め、その周期に対応する時間ステップの項により構成される自己回帰モデルにより予測値を算出する。定常性判定装置は、周期に対応する時間ステップに着目することにより、一周期あたりのサンプル数が多い定常的な信号の自己回帰モデルによる予測値を、高精度かつ低計算コストに求めることができる。ここで、一周期あたりのサンプル数が多いとは、上記文献に記載された技術では計算コストの問題によりリアルタイムな動作が困難となるサンプル数を意味し、サンプル間隔や周期の長さに依存する値である。サンプル間隔が短い場合は短時間でもサンプル数が多く、サンプル間隔が長い場合は長時間でもサンプル数が少なくなり得る。定常性判定装置は、この予測値を用いてセンサ信号の定常性を判定する。   Therefore, the stationarity determination device according to the present embodiment can improve the technique disclosed in Non-Patent Document 1 and determine the stationarity of a sensor signal in real time while accurately reducing the amount of calculation. I made it. Specifically, the stationarity determination device obtains a cycle of the sensor signal by a sequential autocovariance function, and calculates a predicted value by an autoregressive model configured by terms of time steps corresponding to the cycle. The stationarity determination device can obtain a predicted value based on an autoregressive model of a stationary signal with a large number of samples per cycle by focusing on the time step corresponding to the cycle with high accuracy and low calculation cost. . Here, a large number of samples per cycle means the number of samples that is difficult to operate in real time due to the problem of calculation cost with the technique described in the above document, and depends on the sample interval and the length of the cycle. Value. If the sample interval is short, the number of samples is large even in a short time, and if the sample interval is long, the number of samples can be small even in a long time. The stationarity determination device determines the stationarity of the sensor signal using the predicted value.

定常性判定装置による処理について、より詳しく説明する。定常性判定装置は、センサから出力された時系列信号から、計算量が比較的小さい逐次自己共分散関数を利用して、信号の周期を抽出する。次いで、定常性判定装置は、抽出した周期に対応する時間ステップの自己回帰係数を用いて、逐次的に自己回帰係数を算出する。そして、定常性判定装置は、算出した自己回帰係数により構成される自己回帰モデルによる予測値と実際の値との差を予測誤差とし、予測誤差を信号の揺らぎとして定量的に定常性を判定する。   The process by the continuity determination device will be described in more detail. The stationarity determination device extracts a signal period from a time-series signal output from a sensor using a sequential autocovariance function with a relatively small amount of calculation. Next, the stationarity determination apparatus sequentially calculates autoregressive coefficients using autoregressive coefficients of time steps corresponding to the extracted periods. The stationarity determination device quantitatively determines stationarity using the difference between the predicted value and the actual value of the autoregressive model configured by the calculated autoregressive coefficient as a prediction error, and using the prediction error as a signal fluctuation. .

本実施形態では、所望信号と同じような周波数のノイズが存在する場合であっても、周期の揺らぎを評価することにより、所望信号の定常性を判定することができる。また、抽出した周期に対応する時間ステップの自己回帰係数を用いることによって、一周期あたりのサンプル数が多い定常信号を精度よくかつ計算量が少なく予測することが可能となる。これらによって、本実施形態では、定常的な信号と非定常な信号を区別する精度が向上する。   In the present embodiment, even when noise having the same frequency as that of the desired signal exists, the continuity of the desired signal can be determined by evaluating the fluctuation of the period. In addition, by using the autoregressive coefficient of the time step corresponding to the extracted period, it is possible to predict a steady signal having a large number of samples per period with high accuracy and a small amount of calculation. As a result, in this embodiment, the accuracy of distinguishing between stationary signals and non-stationary signals is improved.

以上、定常性判定装置の概要を説明した。続いて、図1を参照して、定常性判定装置の構成例を説明する。   The overview of the stationarity determination device has been described above. Next, a configuration example of the continuity determination device will be described with reference to FIG.

<2.定常性判定装置の構成例>
図1は、本発明の一実施形態に係る定常性判定装置の構成の一例を示すブロック図である。図1に示すように、定常性判定装置1は、センサ10、AD(Analog−to−Digital)コンバータ20、信号周期抽出部30、定常性判定部40、および出力部50を有する。
<2. Configuration example of stationarity determination device>
FIG. 1 is a block diagram illustrating an example of a configuration of a continuity determination device according to an embodiment of the present invention. As shown in FIG. 1, the stationarity determination device 1 includes a sensor 10, an AD (Analog-to-Digital) converter 20, a signal period extraction unit 30, a stationarity determination unit 40, and an output unit 50.

(センサ10)
センサ10は、任意の対象空間または対象物体を観測して、観測された時系列信号(センサ信号)を出力する装置である。センサ10は、例えば振動センサや加速度センサ、ドップラーセンサとして実現される。本実施形態では、センサ信号はアナログ信号であるものとする。
(Sensor 10)
The sensor 10 is an apparatus that observes an arbitrary target space or target object and outputs an observed time series signal (sensor signal). The sensor 10 is realized as a vibration sensor, an acceleration sensor, or a Doppler sensor, for example. In the present embodiment, it is assumed that the sensor signal is an analog signal.

(ADコンバータ20)
ADコンバータ20は、アナログ電気信号をデジタル電気信号に変換する電子回路である。ADコンバータ20は、センサ10から出力されたセンサ信号を、アナログ信号からデジタル信号に変換して出力する。
(AD converter 20)
The AD converter 20 is an electronic circuit that converts an analog electric signal into a digital electric signal. The AD converter 20 converts the sensor signal output from the sensor 10 from an analog signal to a digital signal and outputs the signal.

(信号周期抽出部30)
信号周期抽出部30は、センサ10から出力された時系列信号の周期を抽出する抽出部としての機能を有する。信号周期抽出部30は、定常性判定部40によるリアルタイムな定常性の判定精度を担保するため、逐次的に時系列信号の周期を抽出する。信号周期抽出部30が、逐次的に周期を抽出する手段は多様に考えられる。例えば、信号周期抽出部30は、逐次自己共分散関数または逐次自己相関関数のピークをとる時間ステップを、周期として抽出する。ここでいうピークとは、最大値または極大値の少なくともいずれかである。最大値は、信号のメインの周期を表す。極大値は、異なる周期を有する複数の信号が混在している場合に、その個々の信号の周期を表す。なお、センサ信号の各時間ステップと現時刻との時間間隔を、ラグタイムとも称する。
(Signal period extraction unit 30)
The signal period extraction unit 30 has a function as an extraction unit that extracts the period of the time-series signal output from the sensor 10. The signal period extraction unit 30 sequentially extracts the period of the time-series signal in order to ensure real-time continuity determination accuracy by the continuity determination unit 40. Various means can be conceived for the signal period extraction unit 30 to sequentially extract the periods. For example, the signal period extraction unit 30 extracts a time step that takes a peak of a sequential autocovariance function or a sequential autocorrelation function as a period. The peak here is at least one of a maximum value and a maximum value. The maximum value represents the main period of the signal. The maximum value represents the period of each signal when a plurality of signals having different periods are mixed. The time interval between each time step of the sensor signal and the current time is also referred to as lag time.

逐次共分散関数および逐次自己相関関数は、計算量を少なく、かつ周波数の分解能を細かくとることができるため、リアルタイムに周期を抽出する手段として適している。逐次自己共分散関数は、相関に分散を乗算したものであるため、信号周期抽出部30は、逐次自己共分散関数を用いた場合、振幅が小さいノイズの影響を除去して周期を抽出することができる。一方で、信号周期抽出部30は、逐次自己相関関数を用いた場合、微小な信号であっても周期を抽出することができる。信号周期抽出部30は、センサ10の種別やノイズの大きさ等に基づいて、逐次自己共分散関数または逐次自己相関関数のいずれを用いるかを選択する。   The sequential covariance function and the sequential autocorrelation function are suitable as means for extracting a period in real time because the amount of calculation is small and the frequency resolution can be finely taken. Since the sequential autocovariance function is obtained by multiplying correlation by variance, the signal period extraction unit 30 extracts the period by removing the influence of noise having a small amplitude when the sequential autocovariance function is used. Can do. On the other hand, when the sequential autocorrelation function is used, the signal period extraction unit 30 can extract the period even for a minute signal. The signal period extraction unit 30 selects whether to use a sequential autocovariance function or a sequential autocorrelation function based on the type of sensor 10, the magnitude of noise, and the like.

信号周期抽出部30は、単数または複数個のピークを抽出する。ここで、予測部43による予測精度を考慮すると、信号周期抽出部30は、時系列信号のすべてのピークを抽出して、予測部43がそのピークに対応する時間ステップの項を有する自己回帰モデルで予測することが望ましい。しかし、計算コストを削減する観点から、信号周期抽出部30は、最大値または極大値の少なくともいずれをとる時間ステップのうち、信号値が閾値を超える時間ステップを周期として抽出する。また、計算コストを削減する観点から、信号周期抽出部30は、期間を区切り、所定区間の時系列信号からピークを抽出してもよい。   The signal period extraction unit 30 extracts one or a plurality of peaks. Here, in consideration of the prediction accuracy by the prediction unit 43, the signal period extraction unit 30 extracts all the peaks of the time-series signal, and the prediction unit 43 has a time step term corresponding to the peak. It is desirable to predict with. However, from the viewpoint of reducing the calculation cost, the signal period extraction unit 30 extracts, as a period, a time step in which the signal value exceeds the threshold among time steps that take at least one of the maximum value and the maximum value. Further, from the viewpoint of reducing the calculation cost, the signal period extraction unit 30 may extract a peak from a time-series signal in a predetermined section by dividing the period.

逐次自己共分散関数や逐次自己相関関数を用いた方法の他にも、信号周期抽出部30は、時系列信号についての既知の周期性に対応する時間ステップを、周期として抽出してもよい。既知の周期性とは、例えばフーリエ変換やゼロクロス値等により抽出される周期性である。   In addition to the method using the sequential autocovariance function or the sequential autocorrelation function, the signal period extraction unit 30 may extract a time step corresponding to a known periodicity for the time series signal as a period. The known periodicity is a periodicity extracted by, for example, Fourier transform or zero cross value.

(定常性判定部40)
定常性判定部40は、信号周期抽出部30により抽出された周期を用いて、センサ10から出力された時系列信号の定常性を判定する機能を有する。図1に示すように、定常性判定部40は、係数算出部41、予測部43、および判定部45として機能する。
(Stability determination unit 40)
The stationarity determination unit 40 has a function of determining the stationarity of the time-series signal output from the sensor 10 using the period extracted by the signal period extraction unit 30. As illustrated in FIG. 1, the continuity determination unit 40 functions as a coefficient calculation unit 41, a prediction unit 43, and a determination unit 45.

・係数算出部41
係数算出部41は、信号周期抽出部30により抽出された周期に対応する時間ステップの自己回帰係数を逐次的に算出する機能を有する。詳しくは、係数算出部41は、単位数量の時間ステップ分の時系列信号がセンサ10から出力される都度、自己回帰係数を算出する。単位数量の時間ステップとしては、例えば1ステップでもよいし、所定秒数分の時間ステップであってもよい。
Coefficient calculation unit 41
The coefficient calculation unit 41 has a function of sequentially calculating autoregressive coefficients of time steps corresponding to the period extracted by the signal period extraction unit 30. Specifically, the coefficient calculation unit 41 calculates an autoregressive coefficient each time a time-series signal for a unit quantity of time steps is output from the sensor 10. The unit quantity time step may be, for example, one step or a time step of a predetermined number of seconds.

係数算出部41が、逐次的に自己回帰係数を抽出する手段は多様に考えられる。例えば、係数算出部41は、カルマンフィルタにより自己回帰係数を算出する。カルマンフィルタ以外に逐次的に自己回帰係数を算出する手段として、例えばSDAR(sequentially discounting AR model estimating)アルゴリズムが考えられる。しかし、SDARアルゴリズムでは、すべての時間ステップの自己回帰係数を算出する必要があるため、信号周期抽出部30により抽出された周期に対応する時間ステップの自己回帰係数のみを選択的に算出することができない。よって、係数算出部41は、カルマンフィルタを用いた方が、SDARアルゴリズムを用いる場合と比較して、低い計算コストで自己回帰係数を算出することができる。他にも、係数算出部41は、例えばパーティクルフィルタにより、自己回帰係数を算出してもよい。   There are various means by which the coefficient calculation unit 41 sequentially extracts autoregressive coefficients. For example, the coefficient calculation unit 41 calculates an autoregressive coefficient using a Kalman filter. As a means for sequentially calculating autoregressive coefficients other than the Kalman filter, for example, a SDAR (Sequentially Distinguishing AR model Estimating) algorithm can be considered. However, in the SDAR algorithm, since it is necessary to calculate the autoregressive coefficients of all time steps, it is possible to selectively calculate only the autoregressive coefficients of the time steps corresponding to the periods extracted by the signal period extracting unit 30. Can not. Therefore, the coefficient calculation unit 41 can calculate the autoregressive coefficient at a lower calculation cost when the Kalman filter is used than when the SDAR algorithm is used. In addition, the coefficient calculation unit 41 may calculate the autoregressive coefficient using, for example, a particle filter.

係数算出部41が算出した自己回帰係数は、信号周期抽出部30により抽出された周期に対応する時間ステップの自己回帰係数であり、後述の予測部43による予測に用いられる自己下記モデルを構成する。このため、自己回帰モデルは時系列信号の周期性を良く表現するものとなるので、予測部43は、一周期あたりのサンプル数が多い定常信号を精度よく予測することができる。   The autoregressive coefficient calculated by the coefficient calculating unit 41 is an autoregressive coefficient of a time step corresponding to the period extracted by the signal period extracting unit 30, and constitutes the following self model used for prediction by the prediction unit 43 described later. . For this reason, since the autoregressive model expresses the periodicity of the time-series signal well, the prediction unit 43 can accurately predict a steady signal having a large number of samples per cycle.

また、係数算出部41は、信号周期抽出部30により抽出された周期に対応する時間ステップの自己回帰係数に加えて、任意の時間ステップの自己回帰係数を逐次的に算出してもよい。例えば、係数算出部41は、現時刻の直近の数ステップ分の自己回帰係数を算出する。これにより、予測部43は、直近の数ステップについての観測値を加味した予測を行うことができる。予測部43による予測精度を考慮すると、信号周期抽出部30が時系列信号のすべてのピークを抽出して、予測部43がそのピークに対応する時間ステップの項を有する自己回帰モデルで予測することが望ましい。しかし、計算コストを削減する観点から、信号周期抽出部30は期間を区切ってピークを抽出することが考えられる。このような場合に、予測部43が、直近の数ステップについての観測値を加味した予測を行うことにより、期間を区切ったことによる予測精度の低下を補償することができる。   The coefficient calculation unit 41 may sequentially calculate an autoregressive coefficient at an arbitrary time step in addition to the autoregressive coefficient at a time step corresponding to the period extracted by the signal period extraction unit 30. For example, the coefficient calculation unit 41 calculates autoregressive coefficients for several steps nearest to the current time. Thereby, the prediction part 43 can perform the prediction which considered the observation value about the last several steps. In consideration of the prediction accuracy by the prediction unit 43, the signal period extraction unit 30 extracts all the peaks of the time-series signal, and the prediction unit 43 predicts with an autoregressive model having a time step term corresponding to the peak. Is desirable. However, from the viewpoint of reducing the calculation cost, it is conceivable that the signal period extraction unit 30 extracts peaks by dividing the period. In such a case, the prediction unit 43 can compensate for a decrease in prediction accuracy due to the division of the period by performing the prediction in consideration of the observation values for the most recent steps.

・予測部43
予測部43は、係数算出部41により算出された自己回帰係数を用いた自己回帰モデルにより、時系列信号の予測を行う機能を有する。予測部43は、後述の数式13により、時系列信号の予測を行う。数式1および数式13を参照すると、予測部43は、上記非特許文献1に記載された技術と同様に、自己回帰モデル(時変係数自己回帰モデル)により時系列信号の予測を行っている。ただし、上記非特許文献1では、上記数式1に示すように、対象期間のすべての時間ステップについて計算する必要がある。これに対し、予測部43は、数式13に示すように、信号周期抽出部30により抽出された周期の時間ステップについてのみ計算すればよい。このため、本実施形態では、上記非特許文献1に記載された技術と比較して、計算量を削減することができる。予測部43は、予測結果を判定部45に出力する。
Prediction unit 43
The prediction unit 43 has a function of predicting a time series signal using an autoregressive model using the autoregressive coefficient calculated by the coefficient calculating unit 41. The prediction unit 43 predicts a time-series signal using Equation 13 described later. Referring to Equation 1 and Equation 13, the prediction unit 43 predicts a time-series signal using an autoregressive model (time-varying coefficient autoregressive model) as in the technique described in Non-Patent Document 1. However, in the said nonpatent literature 1, as shown to the said Numerical formula 1, it is necessary to calculate about all the time steps of an object period. On the other hand, the prediction unit 43 only has to calculate only the time step of the period extracted by the signal period extraction unit 30 as shown in Equation 13. For this reason, in this embodiment, compared with the technique described in the said nonpatent literature 1, calculation amount can be reduced. The prediction unit 43 outputs the prediction result to the determination unit 45.

・判定部45
判定部45は、予測部43により予測された時系列信号の予測値と、センサ10から出力された実際の値(観測値)との比較結果に基づいて、時系列信号の定常性を判定する機能を有する。より詳しくは、判定部45は、予測値と実際の値との差分が閾値以上である場合に非定常であると判定し、閾値未満である場合に定常であると判定する。予測値と実際の値との差分を、以下では予測誤差とも称する。判定部45は、判定結果を出力部50に出力する。
・ Determining unit 45
The determination unit 45 determines the continuity of the time-series signal based on the comparison result between the predicted value of the time-series signal predicted by the prediction unit 43 and the actual value (observed value) output from the sensor 10. It has a function. More specifically, the determination unit 45 determines that the difference between the predicted value and the actual value is not less than the threshold value, and determines that the difference is less than the threshold value. The difference between the predicted value and the actual value is hereinafter also referred to as a prediction error. The determination unit 45 outputs the determination result to the output unit 50.

(出力部50)
出力部50は、映像、画像、音声などによって、定常性判定部40による判定結果を出力する。出力部50は、例えばCRT(Cathode Ray Tube)ディスプレイ装置、液晶ディスプレイ(Liquid Crystal Display)装置、スピーカー等により実現される。
(Output unit 50)
The output unit 50 outputs the determination result by the continuity determination unit 40 by video, image, audio, or the like. The output unit 50 is realized by, for example, a CRT (Cathode Ray Tube) display device, a liquid crystal display (Liquid Crystal Display) device, a speaker, or the like.

以上、本実施形態に係る定常性判定装置の構成例を説明した。   The configuration example of the continuity determination device according to the present embodiment has been described above.

<3.定常性判定装置の動作>
本節では、扇風機のような機械的で周期的な動きをする物体が存在する環境で、人の有無を判定する場合の動作について説明する。人が存在する場合のセンサ信号は非定常信号となり、人が存在しない場合のセンサ信号は定常信号となる。定常性判定装置1の動作は、第1にセンサ10からの信号取得およびその変換、第2に信号周期の抽出、第3に自己回帰係数を用いた予測および予測誤差による定常性の判定、の3つの段階に分類される。以下、図2〜図4を参照して、各段階の動作についてそれぞれ説明する。なお、図2は、本発明の一実施形態に係る定常性判定装置1の動作を示すフローチャートである。
<3. Operation of Stationaryness Determination Device>
In this section, the operation for determining the presence or absence of a person in an environment where there is an object that moves mechanically and periodically, such as a fan, will be described. The sensor signal when a person is present is an unsteady signal, and the sensor signal when no person is present is a steady signal. The operation of the stationarity determination device 1 includes first acquisition of a signal from the sensor 10 and conversion thereof, second extraction of a signal period, third prediction using an autoregressive coefficient and determination of stationarity based on a prediction error. There are three stages. Hereinafter, the operation of each stage will be described with reference to FIGS. FIG. 2 is a flowchart showing the operation of the stationarity determination device 1 according to an embodiment of the present invention.

[3−1.センサからの信号取得およびその変換]
図2に示すように、まず、ステップS102で、定常性判定装置1は、センサ信号を取得する。詳しくは、センサ10は、センシング結果をアナログの時系列信号として出力する。
[3-1. Signal acquisition from sensor and its conversion]
As shown in FIG. 2, first, in step S102, the stationarity determination device 1 acquires a sensor signal. Specifically, the sensor 10 outputs the sensing result as an analog time series signal.

次いで、ステップS104で、定常性判定装置1は、センサ信号をデジタル信号に変換する。詳しくは、ADコンバータ20は、センサ10から出力されたアナログの時系列信号をデジタル信号に変換する。なお、センサ10により取得されるアナログ信号は一般に微小信号となる。このため、定常性判定装置1は、雑音の影響を抑制しSN比を改善するために、増幅器によって増幅したり、アナログフィルタやデジタルフィルタによって周波数フィルタリングを行ったりしてもよい。   Next, in step S104, the stationarity determination device 1 converts the sensor signal into a digital signal. Specifically, the AD converter 20 converts an analog time series signal output from the sensor 10 into a digital signal. The analog signal acquired by the sensor 10 is generally a minute signal. For this reason, in order to suppress the influence of noise and improve the S / N ratio, the stationarity determination apparatus 1 may amplify with an amplifier or perform frequency filtering with an analog filter or a digital filter.

[3−2.信号周期の抽出]
次に、ステップS106で、定常性判定装置1は、各ラグタイムにおける逐次自己共分散関数を求める。詳しくは、信号周期抽出部30は、ADコンバータ20から出力されたセンサ信号に基づいて、逐次自己共分散関数を求める。ある時刻nのラグタイムτにおける逐次自己共分散関数は、以下の数式3で表される。より具体的には、信号周期抽出部30は、観測値y(n=1,2,・・・N)について、ある時刻nにおける逐次平均u、ある時刻nのあるラグタイムτにおける逐次自己共分散関数C(n,τ)を、それぞれ以下の数式2、数式3の漸化式で計算する。
[3-2. Extraction of signal period]
Next, in step S106, the stationarity determination apparatus 1 obtains a sequential autocovariance function at each lag time. Specifically, the signal period extraction unit 30 sequentially obtains an autocovariance function based on the sensor signal output from the AD converter 20. The sequential autocovariance function at a lag time τ at a certain time n is expressed by the following Equation 3. More specifically, the signal period extraction unit 30 sequentially determines the observed value y n (n = 1, 2,... N) at the sequential average u n at a certain time n and the lag time τ at a certain time n. The autocovariance function C (n, τ) is calculated by the recurrence formulas of the following formulas 2 and 3, respectively.

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

ここで、rは忘却係数であり、0<r<1の範囲の値をとる。数式2と数式3で異なる忘却係数をとっても良い。信号周期抽出部30は、信号の周期が存在する期間を0≦τ≦τmaxとして、その時間範囲における逐次自己共分散関数を求める。 Here, r is a forgetting factor and takes a value in the range of 0 <r <1. Different forgetting factors may be used in Equation 2 and Equation 3. The signal period extraction unit 30 obtains a sequential autocovariance function in the time range, where 0 ≦ τ ≦ τ max is a period in which the signal period exists.

次いで、ステップS108で、定常性判定装置1は、逐次自己共分散関数が最大値または極大値となるラグライムを単数または複数抽出する。詳しくは、信号周期抽出部30は、上記数式3で示した逐次自己共分散関数が、最大値となるラグタイムτまたは極大値となるラグタイムτを、単数または複数抽出する。   Next, in step S108, the stationarity determination apparatus 1 extracts one or a plurality of raglimes whose sequential autocovariance function has a maximum value or a maximum value. Specifically, the signal period extraction unit 30 extracts one or a plurality of lag times τ at which the sequential autocovariance function expressed by Equation 3 has a maximum value or a lag time τ at which the maximum value is reached.

[3−3.自己回帰係数を用いた予測および予測誤差による定常性の判定]
次に、ステップS110で、定常性判定装置1は、抽出されたラグタイムに対応する時変係数自己回帰モデルを構成する。詳しくは、まず、係数算出部41は、上記ステップS108において信号周期抽出部30により抽出された周期に対応する時間ステップの自己回帰係数を、逐次的に算出する。係数算出部41は、以下に説明するように、カルマンフィルタにより自己回帰係数を算出するものとする。
[3-3. Prediction using autoregressive coefficient and determination of stationarity by prediction error]
Next, in step S110, the stationarity determination device 1 configures a time-varying coefficient autoregressive model corresponding to the extracted lag time. Specifically, first, the coefficient calculation unit 41 sequentially calculates the autoregressive coefficient of the time step corresponding to the period extracted by the signal period extraction unit 30 in step S108. The coefficient calculation unit 41 calculates autoregressive coefficients using a Kalman filter as described below.

カルマンフィルタは、m次元の未知のパラメータから構成される状態ベクトルxとスカラ観測値yとの関係を与える、数式4で示す状態方程式および数式5で示す観測方程式でモデル化される。   The Kalman filter is modeled by a state equation represented by Equation 4 and an observation equation represented by Equation 5, which give a relationship between a state vector x composed of m-dimensional unknown parameters and a scalar observation value y.

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

ここで、uはm次元のシステム雑音、vは1次元の観測ノイズであり、以下の数式6、数式7のように、平均0、分散共分散行列がそれぞれQ、Rの正規分布に従うと仮定する。 Here, u n is m-dimensional system noise, v n is the one-dimensional observation noise, the following equation 6, as in Equation 7, the mean 0, variance-covariance matrix follows Q, the normal distribution of the R, respectively Assume that

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

ここで、Iは単位行列である。Kはカルマンゲインと呼ばれ、以下の数式12のように表される。 Here, I is a unit matrix. K n is called the Kalman gain is expressed by the following equation 12.

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

Figure 0006299172
Figure 0006299172

ユールウォーカ方程式で自己回帰係数を求める場合はラグステップが1から順番に求めていく必要があるが、上記数式13に示すように、カルマンフィルタを用いた場合は必要な序数の係数のみを求めることができる。即ち、係数算出部41は、信号周期抽出部30により抽出されたラグタイムに対応する自己回帰係数のみを計算対象とすることにより、より低い計算コストで自己回帰係数を算出することができる。以上説明したように、ステップS110において、係数算出部41は、信号周期抽出部30により抽出されたラグタイムに対応する時変係数自己回帰モデルとして、上記数式13に示した自己回帰モデルを構成する。   When obtaining the autoregressive coefficient using the Yule-Walker equation, it is necessary to obtain the lag step in order from 1. However, as shown in the above equation 13, if the Kalman filter is used, only the necessary ordinal coefficient can be obtained. it can. That is, the coefficient calculating unit 41 can calculate the autoregressive coefficient at a lower calculation cost by using only the autoregressive coefficient corresponding to the lag time extracted by the signal period extracting unit 30 as a calculation target. As described above, in step S110, the coefficient calculation unit 41 configures the autoregressive model shown in Equation 13 as a time-varying coefficient autoregressive model corresponding to the lag time extracted by the signal period extraction unit 30. .

そして、ステップS112で、定常性判定装置1は、予測誤差を求める。詳しくは、まず、予測部43は、上記数式13に示した自己回帰モデルに基づき予測値を求める。そして、判定部45は、予測部43により予測された予測値とセンサ10により出力された実際の値との差を、以下の数式14で表される予測誤差eとして求める。 In step S112, the stationarity determination device 1 obtains a prediction error. Specifically, first, the prediction unit 43 obtains a predicted value based on the autoregressive model shown in the mathematical formula 13. Then, the determination unit 45, the difference between the actual value output by the prediction value and the sensor 10 predicted by the prediction unit 43 is obtained as the prediction error e n represented by the following equation 14.

Figure 0006299172
Figure 0006299172

ここで、図3、図4を参照して、現時刻から1ステップ前からNステップ前までの自己回帰係数により構成された自己回帰モデル(比較例)と、信号の周期に対応する自己回帰係数により構成された自己回帰モデル(本実施形態)との、予測誤差の違いを説明する。図3は、比較例に係る定常性判定装置による予測誤差を示す図である。図4は、本発明の一実施形態に係る定常性判定装置1による予測誤差を示す図である。   Here, referring to FIG. 3 and FIG. 4, an autoregressive model (comparative example) composed of autoregressive coefficients from one step before to N steps before the current time, and an autoregressive coefficient corresponding to the period of the signal The difference in prediction error from the autoregressive model (this embodiment) configured by FIG. 3 is a diagram illustrating a prediction error by the continuity determination device according to the comparative example. FIG. 4 is a diagram showing a prediction error by the stationarity determination device 1 according to an embodiment of the present invention.

図3では、比較例として、次数を3として、直近の3つの時間ステップd1、d2、d3についての自己回帰モデルによる予測誤差を示している。図3に示すように、過去の直近の観測値を用いて予測されるため、自己回帰モデルは対象とした3つの時間ステップを超える周期を表現することができず、予測値が直近の観測値に依存してしまうため、予測誤差が大きくなってしまう場合がある。特に、対象とした3つの時間ステップ内またはその前後において信号が局所的に不規則である場合、予測誤差が大きくなってしまう。   In FIG. 3, as a comparative example, the order is 3, and prediction errors based on the autoregressive model for the three most recent time steps d1, d2, and d3 are shown. As shown in FIG. 3, since the most recent observation value in the past is used for prediction, the autoregressive model cannot express the period exceeding the three time steps of interest, and the prediction value is the latest observation value. The prediction error may increase. In particular, if the signal is locally irregular within or around the three time steps of interest, the prediction error will be large.

図4では、次数を3として、直近の2つの時間ステップd2およびd3、逐次自己共分散関数により抽出した周期に対応する時間ステップd1についての自己回帰モデルによる予測誤差を示している。図4に示すように、逐次自己共分散関数により抽出された周期に対応する項を自己回帰モデルに採用することにより、自己回帰モデルは信号の周期性を表現することが可能となり、図3に示した比較例と比較して予測誤差を小さくすることができる。この自己回帰モデルは、直近の観測値に過度に依存しないため、対象とした3つの時間ステップ内またはその前後において信号が局所的に不規則である場合であっても、予測誤差が大きくなることを抑制することができる。   FIG. 4 shows the prediction error by the autoregressive model for the most recent two time steps d2 and d3 and the time step d1 corresponding to the period extracted by the sequential autocovariance function, where the order is 3. As shown in FIG. 4, by adopting a term corresponding to the period extracted by the sequential autocovariance function in the autoregressive model, the autoregressive model can express the periodicity of the signal. The prediction error can be reduced as compared with the comparative example shown. This autoregressive model does not depend excessively on the most recent observations, so the prediction error increases even if the signal is locally irregular within or around the three time steps of interest. Can be suppressed.

以上、ステップS112における処理について説明した。次いで、ステップS114で、定常性判定装置1は、予測誤差が閾値以上であるか否かを判定する。詳しくは、判定部45は、数式14で示した予測誤差が閾値以上であるか否かを判定する。   The processing in step S112 has been described above. Next, in step S114, the stationarity determination device 1 determines whether the prediction error is greater than or equal to a threshold value. Specifically, the determination unit 45 determines whether or not the prediction error expressed by Equation 14 is greater than or equal to a threshold value.

予測誤差が以上であると判定された場合(S114/YES)、ステップS116で、定常性判定部40は、センサ信号は非定常信号であると判定する。例えば、センサ10のセンシング対象空間内に人が存在する場合、人の呼吸や動作に起因する信号が観測されるため、定常性判定部40は、センサ信号は非定常信号であると判定する。   If it is determined that the prediction error is greater than or equal to (S114 / YES), in step S116, the stationarity determination unit 40 determines that the sensor signal is an unsteady signal. For example, when a person is present in the sensing target space of the sensor 10, a signal resulting from the person's breathing or movement is observed, so the continuity determination unit 40 determines that the sensor signal is an unsteady signal.

一方で、予測誤差が未満であると判定された場合(S114/NO)、ステップS118で、定常性判定部40は、センサ信号は定常信号であると判定する。例えば、センサ10のセンシング対象空間内に人が存在しない場合、扇風機の周期的な動きに起因する定常信号のみが観測されるため、定常性判定部40は、センサ信号は定常信号であると判定する。   On the other hand, when it is determined that the prediction error is less (S114 / NO), in step S118, the stationarity determination unit 40 determines that the sensor signal is a stationar signal. For example, when there is no person in the sensing target space of the sensor 10, only the stationary signal due to the periodic movement of the electric fan is observed, so the continuity determination unit 40 determines that the sensor signal is a stationary signal. To do.

そして、ステップS120で、定常性判定装置1は、判定結果を出力する。詳しくは、出力部50は、上記ステップS116またはS118において判定部45により判定された判定結果を、出力する。   In step S120, the stationarity determination device 1 outputs a determination result. Specifically, the output unit 50 outputs the determination result determined by the determination unit 45 in step S116 or S118.

次に、ステップS122で、定常性判定装置1は、処理を終了するか否かを判定する。
終了しないと判定された場合(S122/NO)、処理は再度ステップS102に戻り、定常性判定装置1は定常性の判定を繰り返す。これにより、定常性判定装置1は、リアルタイムにセンサ信号の定常性を判定することができる。一方で、終了すると判定された場合(S122/YES)、定常性判定装置1は処理を終了する。
Next, in step S122, the stationarity determination device 1 determines whether or not to end the process.
When it is determined not to end (S122 / NO), the process returns to step S102 again, and the stationarity determination device 1 repeats the stationarity determination. Thereby, the stationarity determination apparatus 1 can determine the stationarity of the sensor signal in real time. On the other hand, when it determines with complete | finishing (S122 / YES), the stationarity determination apparatus 1 complete | finishes a process.

以上、本実施形態に係る定常性判定装置1の動作を説明した。   The operation of the continuity determination device 1 according to this embodiment has been described above.

(補足)
上述したように、信号周期抽出部30は、逐次自己共分散関数の他にも、逐次自己相関関数を用いて周期を抽出することができる。この場合、信号周期抽出部30は、上記数式3に代えて、数式3で示した時刻nにおける逐次自己共分散関数C(n,τ)を逐次分散v(n)で除した、下記数式15で示す逐次自己相関関数Rxx(n,τ)を用いる。
(Supplement)
As described above, the signal period extraction unit 30 can extract the period using the sequential autocorrelation function in addition to the sequential autocovariance function. In this case, the signal period extraction unit 30 replaces the equation 3 above with the sequential autocovariance function C (n, τ) at time n shown in the equation 3 divided by the sequential variance v (n). The sequential autocorrelation function Rxx (n, τ) shown in FIG.

Figure 0006299172
Figure 0006299172

ここで、逐次分散v(n)は、下記の数式16で示す漸化式で与えられる。   Here, the sequential variance v (n) is given by a recurrence formula shown by the following formula 16.

Figure 0006299172
Figure 0006299172

ここで、rは忘却係数であり、0<r<1の範囲の値をとる。   Here, r is a forgetting factor and takes a value in the range of 0 <r <1.

<4.まとめ>
以上説明したように、本実施形態によれば、センサ信号の定常性をより少ない計算量で精度良くリアルタイムに判定することが可能である。より詳しくは、定常性判定装置1は、逐次自己共分散関数または逐次自己相関関数を用いることにより、フーリエ変換等を用いた場合と比較してより少ない計算量でセンサ信号の周期を抽出することができる。また、定常性判定装置1は、抽出した周期に対応する時間ステップについての自己回帰モデルを構築することにより、周期的な信号の予測誤差を低減することができるため、より精度よくセンサ信号の定常性を判定することができる。
<4. Summary>
As described above, according to the present embodiment, it is possible to determine the continuity of a sensor signal with high accuracy and in real time with a smaller amount of calculation. More specifically, the stationarity determination device 1 uses the sequential autocovariance function or the sequential autocorrelation function to extract the period of the sensor signal with a smaller amount of calculation compared to the case where Fourier transform or the like is used. Can do. In addition, the stationarity determination apparatus 1 can reduce the prediction error of a periodic signal by constructing an autoregressive model with respect to a time step corresponding to the extracted period, so that the steady state of the sensor signal can be improved more accurately. Gender can be determined.

定常性判定装置1は、センサ信号の定常性を判定することにより、例えばドップラーセンサからの出力結果に基づいて、扇風機のように一定の周期で機械的な動作をする機器の動きと、人のような周期の揺らぎが大きい不規則動きとを区別することが可能である。また、定常性判定装置1は、人の有無を高精度で検知する人感センサにも応用することができる。   The stationarity determination device 1 determines the stationarity of the sensor signal, for example, based on the output result from the Doppler sensor, the movement of a device that performs a mechanical operation at a constant period, such as an electric fan, It is possible to distinguish from such irregular movement with a large fluctuation of the period. The stationarity determination device 1 can also be applied to a human sensor that detects the presence or absence of a person with high accuracy.

以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。   The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to such examples. It is obvious that a person having ordinary knowledge in the technical field to which the present invention pertains can come up with various changes or modifications within the scope of the technical idea described in the claims. Of course, it is understood that these also belong to the technical scope of the present invention.

例えば、上記実施形態では、センサ10、ADコンバータ20、信号周期抽出部30、定常性判定部40、および出力部50が一体化しているものとしたが、本発明はかかる例に限定されない。例えば、各構成要素が別々に構成されていてもよい。具体的には、PC(Personal Conputer)やモバイル端末、ネットワーク上のサーバ等が信号周期抽出部30および定常性判定部40として機能し、センサ10、ADコンバータ20からの出力に基づいて判定結果を出力部50に出力してもよい。   For example, in the above embodiment, the sensor 10, the AD converter 20, the signal period extraction unit 30, the continuity determination unit 40, and the output unit 50 are integrated, but the present invention is not limited to such an example. For example, each component may be configured separately. Specifically, a PC (Personal Computer), a mobile terminal, a server on the network, and the like function as the signal period extraction unit 30 and the continuity determination unit 40, and the determination result is obtained based on the output from the sensor 10 and the AD converter 20. You may output to the output part 50.

また、情報処理装置に内蔵されるCPU、ROM及びRAM等のハードウェアに、上記定常性判定装置1の各構成と同等の機能を発揮させるためのコンピュータプログラムも作成可能である。また、当該コンピュータプログラムを記録した記録媒体も提供される。   Further, it is possible to create a computer program for causing hardware such as a CPU, a ROM, and a RAM built in the information processing apparatus to perform the same functions as the components of the continuity determination apparatus 1. A recording medium recording the computer program is also provided.

1 定常性判定装置
10 センサ
20 ADコンバータ
30 信号周期抽出部
40 定常性判定部
41 係数算出部
43 予測部
45 判定部
50 出力部
DESCRIPTION OF SYMBOLS 1 Steadyness determination apparatus 10 Sensor 20 AD converter 30 Signal period extraction part 40 Steadyness determination part 41 Coefficient calculation part 43 Prediction part 45 Determination part 50 Output part

Claims (13)

センサから出力された時系列信号の周期を逐次的に抽出する抽出部と、
前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、
前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、
を備え
前記抽出部は、逐次自己共分散関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、情報処理装置。
An extraction unit that sequentially extracts the period of the time-series signal output from the sensor;
A coefficient calculating unit that sequentially calculates autoregressive coefficients corresponding to the period extracted by the extracting unit;
A prediction unit that predicts the time-series signal by an equation of an autoregressive model using the autoregressive coefficient calculated by the coefficient calculating unit;
Equipped with a,
The information processing apparatus , wherein the extraction unit extracts a time step taking at least one of a maximum value and a maximum value of a sequential autocovariance function as the period .
センサから出力された時系列信号の周期を逐次的に抽出する抽出部と、
前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、
前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、
を備え
前記抽出部は、逐次自己相関関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、情報処理装置。
An extraction unit that sequentially extracts the period of the time-series signal output from the sensor;
A coefficient calculating unit that sequentially calculates autoregressive coefficients corresponding to the period extracted by the extracting unit;
A prediction unit that predicts the time-series signal by an equation of an autoregressive model using the autoregressive coefficient calculated by the coefficient calculating unit;
Equipped with a,
The information processing apparatus , wherein the extraction unit extracts, as the period, a time step that takes at least one of a maximum value and a maximum value of a sequential autocorrelation function .
前記抽出部は、所定区間の前記時系列信号において前記最大値または極大値の少なくともいずれをとる時間ステップのうち、閾値を超える時間ステップを、前記周期として抽出する、請求項1または2に記載の情報処理装置。 The extraction section, among the time step taking any at least of the maximum value or the maximum in the time-series signal of a predetermined period, the time step that exceeds the threshold value is extracted as the period, according to claim 1 or 2 Information processing device. 前記係数算出部は、単位数量の時間ステップ分の前記時系列信号が出力される都度、前記自己回帰係数を算出する、請求項1〜のいずれか一項に記載の情報処理装置。 The coefficient calculation unit, each time the time-series signal time step portion of the unit quantity is output, calculates the autoregressive coefficients, the information processing apparatus according to any one of claims 1-3. 前記係数算出部は、カルマンフィルタにより前記自己回帰係数を算出する、請求項に記載の情報処理装置。 The information processing apparatus according to claim 4 , wherein the coefficient calculation unit calculates the autoregressive coefficient using a Kalman filter. 前記係数算出部は、前記周期に対応する自己回帰係数に加えて、任意の時間ステップの自己回帰係数を逐次的に算出する、請求項1〜のいずれか一項に記載の情報処理装置。 The coefficient calculation unit, in addition to the autoregressive coefficients corresponding to the period to sequentially calculate the autoregressive coefficients of the arbitrary time step, the information processing apparatus according to any one of claims 1-5. 前記情報処理装置は、前記予測部により予測された前記時系列信号の予測値と前記センサから出力された観測値との比較結果に基づいて前記時系列信号の定常性を判定する判定部をさらに備える、請求項1〜のいずれか一項に記載の情報処理装置。 The information processing apparatus further includes a determination unit that determines the continuity of the time series signal based on a comparison result between the predicted value of the time series signal predicted by the prediction unit and the observation value output from the sensor. The information processing apparatus according to any one of claims 1 to 6 , further comprising: 前記判定部は、前記予測値と前記観測値との差分が閾値以上である場合に非定常であると判定し、前記閾値未満である場合に定常であると判定する、請求項に記載の情報処理装置。 The determination unit, the difference between the predicted value and the observed value is determined to be non-stationary if it is more than the threshold value, it determines that the constant in the case is less than the threshold value, according to claim 7 Information processing device. 前記センサは、ドップラーセンサ、振動センサ、加速度センサの少なくともいずれかである、請求項1〜のいずれか一項に記載の情報処理装置。 The sensor Doppler sensor, the vibration sensor is at least one of an acceleration sensor, the information processing apparatus according to any one of claims 1-8. センサから出力された時系列信号の周期を逐次的に抽出するステップと、
抽出された前記周期に対応する自己回帰係数を逐次的に算出するステップと、
算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行うステップと、
を備え
前記抽出するステップは、逐次自己共分散関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、情報処理方法。
Sequentially extracting the period of the time-series signal output from the sensor;
Sequentially calculating autoregressive coefficients corresponding to the extracted periods;
Predicting the time-series signal according to an equation of an autoregressive model using the calculated autoregressive coefficient;
Equipped with a,
The information extracting method, wherein the extracting step extracts a time step taking at least one of a maximum value and a maximum value of a sequential autocovariance function as the period .
センサから出力された時系列信号の周期を逐次的に抽出するステップと、
抽出された前記周期に対応する自己回帰係数を逐次的に算出するステップと、
算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行うステップと、
を備え
前記抽出するステップは、逐次自己相関関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、情報処理方法。
Sequentially extracting the period of the time-series signal output from the sensor;
Sequentially calculating autoregressive coefficients corresponding to the extracted periods;
Predicting the time-series signal according to an equation of an autoregressive model using the calculated autoregressive coefficient;
Equipped with a,
The information extracting method wherein the extracting step extracts a time step taking at least one of a maximum value and a maximum value of a sequential autocorrelation function as the period .
コンピュータを、
センサから出力された時系列信号の周期を逐次的に抽出する抽出部と、
前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、
前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、
として機能させ
前記抽出部は、逐次自己共分散関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、プログラム。
Computer
An extraction unit that sequentially extracts the period of the time-series signal output from the sensor;
A coefficient calculating unit that sequentially calculates autoregressive coefficients corresponding to the period extracted by the extracting unit;
A prediction unit that predicts the time-series signal by an equation of an autoregressive model using the autoregressive coefficient calculated by the coefficient calculating unit;
To function as,
The said extraction part is a program which extracts the time step which takes at least one of the maximum value or maximum value of a sequential autocovariance function as said period .
コンピュータを、
センサから出力された時系列信号の周期を逐次的に抽出する抽出部と、
前記抽出部により抽出された前記周期に対応する自己回帰係数を逐次的に算出する係数算出部と、
前記係数算出部により算出された前記自己回帰係数を用いた自己回帰モデルの式により前記時系列信号の予測を行う予測部と、
として機能させ
前記抽出部は、逐次自己相関関数の最大値または極大値の少なくともいずれかをとる時間ステップを、前記周期として抽出する、プログラム。
Computer
An extraction unit that sequentially extracts the period of the time-series signal output from the sensor;
A coefficient calculating unit that sequentially calculates autoregressive coefficients corresponding to the period extracted by the extracting unit;
A prediction unit that predicts the time-series signal by an equation of an autoregressive model using the autoregressive coefficient calculated by the coefficient calculating unit;
To function as,
The said extraction part is a program which extracts the time step which takes at least one of the maximum value or maximum value of a sequential autocorrelation function as said period .
JP2013236921A 2013-11-15 2013-11-15 Information processing apparatus, information processing method, and program Active JP6299172B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013236921A JP6299172B2 (en) 2013-11-15 2013-11-15 Information processing apparatus, information processing method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013236921A JP6299172B2 (en) 2013-11-15 2013-11-15 Information processing apparatus, information processing method, and program

Publications (2)

Publication Number Publication Date
JP2015096831A JP2015096831A (en) 2015-05-21
JP6299172B2 true JP6299172B2 (en) 2018-03-28

Family

ID=53374191

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013236921A Active JP6299172B2 (en) 2013-11-15 2013-11-15 Information processing apparatus, information processing method, and program

Country Status (1)

Country Link
JP (1) JP6299172B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017022312A1 (en) * 2015-08-06 2017-02-09 新日鉄住金ソリューションズ株式会社 Information processing device, information processing system, information processing method, and program
WO2019026193A1 (en) * 2017-08-02 2019-02-07 日本電気株式会社 Information processing device, information processing system, information processing method, and recording medium
WO2021153188A1 (en) * 2020-01-31 2021-08-05 パナソニックIpマネジメント株式会社 Signal processing system, sensor system, signal processing method, and program

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57127820A (en) * 1981-01-30 1982-08-09 Sumitomo Metal Ind Ltd Method of monitoring cyclically moving body
JPS61212715A (en) * 1985-03-18 1986-09-20 Nippon Kishiyou Kyokai Analysis of wave characteristic
JPH04244997A (en) * 1991-01-31 1992-09-01 Toshiba Corp Monitoring method for stability of core of boiling water reactor
JPH06264704A (en) * 1993-03-10 1994-09-20 Toshiba Corp Vibration diagnostic device for rotational machine
JPH07123589A (en) * 1993-10-19 1995-05-12 Toshiba Corp Demand estimation system
JP2959376B2 (en) * 1994-01-26 1999-10-06 松下電器産業株式会社 Monitoring device
JP2001074616A (en) * 1999-09-06 2001-03-23 Mitsubishi Electric Corp Device for diagnozing abnormality of rotating machine
JP4598809B2 (en) * 2007-08-07 2010-12-15 アプライドリサーチ株式会社 Soundness diagnosis method by sound analysis
JP2009257862A (en) * 2008-04-15 2009-11-05 Original Engineering Consultants Co Ltd Method for diagnosing integrity in equipment such as rotary machine by sound signal
JP5327699B2 (en) * 2008-09-29 2013-10-30 Toto株式会社 Human body detection device and urinal equipped with the same

Also Published As

Publication number Publication date
JP2015096831A (en) 2015-05-21

Similar Documents

Publication Publication Date Title
JP6515670B2 (en) Sleep depth estimation device, sleep depth estimation method, and program
US20160025847A1 (en) State recognizing device, state recognizing method, and recording medium
JP6135280B2 (en) Feature amount calculation apparatus, feature amount calculation method, and program
JP2018173285A (en) Activity recognition system and activity recognition method
JP6056389B2 (en) Heart rate estimation device, heart rate estimation method and program
CN110673135A (en) Sensor, estimation device, estimation method, and program recording medium
JP6299172B2 (en) Information processing apparatus, information processing method, and program
KR20120113530A (en) Method and apparatus for detecting peak from biological signal
JP2016206033A (en) Oscillation state estimating device, oscillation state estimating method, and program
JP2015052491A (en) Signal processing apparatus, signal processing method and computer program
JP6871718B2 (en) Sound source search device, sound source search method and its program
JP6019659B2 (en) Apnea condition determination device, apnea condition determination method, and apnea condition determination program
JP4122003B2 (en) Extraction method and device for biological signals such as heartbeat and respiration
CN107405086B (en) Measurement device, measurement method, and program
JP6536038B2 (en) Period estimation apparatus, period estimation method and program
US20130211273A1 (en) Method and apparatus for heart rate measurement
JP6225079B2 (en) Air conditioner operation detection method and system
CN110801214A (en) Heart rate real-time detection method and system
JP2016182165A (en) Biological signal processing device, biological signal processing program, computer readable recording medium recording biological signal processing program and biological signal processing method
JP2022106705A (en) Autonomous full spectrum biometric monitoring
CN115137329A (en) Heart rate detection method and wearable device
JP6379709B2 (en) Signal processing apparatus, signal processing method, and program
CN113693582A (en) Vital sign information monitoring method and device, storage medium and processor
JP2011139409A (en) Audio signal processor, audio signal processing method, and computer program
CN114533010A (en) Heart rate detection method and device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160816

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170720

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170725

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170907

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180212

R150 Certificate of patent or registration of utility model

Ref document number: 6299172

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150