Nothing Special   »   [go: up one dir, main page]

JP2011045524A - Device for estimating exercise intensity and/or calorie consumption of carrier - Google Patents

Device for estimating exercise intensity and/or calorie consumption of carrier Download PDF

Info

Publication number
JP2011045524A
JP2011045524A JP2009196503A JP2009196503A JP2011045524A JP 2011045524 A JP2011045524 A JP 2011045524A JP 2009196503 A JP2009196503 A JP 2009196503A JP 2009196503 A JP2009196503 A JP 2009196503A JP 2011045524 A JP2011045524 A JP 2011045524A
Authority
JP
Japan
Prior art keywords
state
exercise intensity
acceleration
estimation
representative
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.)
Withdrawn
Application number
JP2009196503A
Other languages
Japanese (ja)
Inventor
Atsunobu Namikawa
敦宣 南川
Arei Kobayashi
亜令 小林
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.)
KDDI Corp
Original Assignee
KDDI Corp
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 KDDI Corp filed Critical KDDI Corp
Priority to JP2009196503A priority Critical patent/JP2011045524A/en
Publication of JP2011045524A publication Critical patent/JP2011045524A/en
Withdrawn legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a device for estimating exercise intensity with high accuracy, with no need to fix in a predetermined direction. <P>SOLUTION: The estimating device includes: a means for setting a plurality of transform periods so that their starting points differ from each other in a predetermined number and subjecting the accelerations included in the respective transform periods to Fourier transform to obtain a plurality of spectra, and substituting each of the plurality of obtained spectra with the representative spectrum with the least error among a plurality of predetermined representative spectra; a means for obtaining appearance probabilities in the plurality of representative spectra obtained by the substitution for each of a plurality of characteristic spectra selected from the plurality of representative spectra in advance, and generating data for estimation including the appearance probabilities of the obtained respective characteristic spectra; a means for determining states corresponding to the data for estimation based on the study data including the same data as the data for estimation measured at least one time in the respective states; and a means for deciding the exercise intensity based on the determined states. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、装置に搭載された加速度センサを利用して、当該装置を所持している者の運動強度及び消費カロリーを推定する技術に関する。   The present invention relates to a technique for estimating exercise intensity and calorie consumption of a person who possesses the device by using an acceleration sensor mounted on the device.

日々の生活における運動強度や消費カロリーを推定する携帯型の装置について各種の提案が行われている(例えば、特許文献1、2及び3、参照。)。   Various proposals have been made on portable devices for estimating exercise intensity and calorie consumption in daily life (see, for example, Patent Documents 1, 2, and 3).

特許文献1によると、推定装置が有する加速度センサの一定期間における出力値の標準偏差値に基づき運動強度を推定している。また、特許文献2によると、推定装置は3軸加速度センサを備えており、3軸加速度センサの各軸方向の出力値から、姿勢成分及び合成加速度の運動成分の平均値を演算し、これら演算した値から、利用者が現在、走行、歩行、階段上昇、階段下降、立位、座位、臥位のいずれの状態にあるかを判定し、さらに、判定した状態に基づきエネルギー消費量を算出している。さらに、特許文献3によると、推定装置は3軸加速度センサを備えており、3軸加速度センサの各軸方向の一定時間における出力値から、各軸方向の代表値を算出し、水平成分の代表値と垂直成分の代表値に基づいて運動強度を推定している。   According to Patent Document 1, the exercise intensity is estimated based on the standard deviation value of the output value of the acceleration sensor included in the estimation device over a certain period. According to Patent Document 2, the estimation device includes a triaxial acceleration sensor, and calculates an average value of the posture component and the motion component of the combined acceleration from the output value in each axis direction of the triaxial acceleration sensor, and calculates these values. From this value, it is determined whether the user is currently running, walking, climbing stairs, descending stairs, standing, sitting, or standing, and energy consumption is calculated based on the determined status. ing. Further, according to Patent Document 3, the estimation apparatus includes a triaxial acceleration sensor, calculates a representative value in each axial direction from an output value of each triaxial acceleration sensor in a certain time, and represents a representative horizontal component. The exercise intensity is estimated based on the value and the representative value of the vertical component.

特開2006−204446号公報JP 2006-204446 A 特開2007−160076号公報JP 2007-160076 A 特開2009−028312号公報JP 2009-028312 A

例えば、静止して体を動かしている状態と走っている状態では運動強度は異なることになるが、特許文献1に記載の方法は、利用者の状態を判定することなく、測定した加速度の一定期間における標準偏差値のみを利用して運動強度を推定するものであり、推定精度が低いという問題がある。また、特許文献2及び3に記載の推定装置は、3軸加速度センサの各軸が所定の方向となる様に利用者に固定する必要があり不便である。特に、推定装置を、他の機能を有する装置、例えば、移動通信端末と一体の装置として実装することを考えると、推定装置のその利用者への固定により、他の機能の利用ができなくなる、あるいは、利用しずらくなるという問題が生じる。つまり、この場合、特許文献2及び3に記載の構成を採用することはできない。   For example, the exercise intensity differs between a state where the body is stationary and a state where the body is running, but the method described in Patent Document 1 does not determine the state of the user, and the measured acceleration is constant. The exercise intensity is estimated by using only the standard deviation value in the period, and there is a problem that the estimation accuracy is low. In addition, the estimation devices described in Patent Documents 2 and 3 are inconvenient because it is necessary to fix the three-axis acceleration sensor to the user so that each axis is in a predetermined direction. In particular, when considering that the estimation device is implemented as a device having other functions, for example, a device integrated with a mobile communication terminal, other functions cannot be used due to fixing of the estimation device to the user. Or the problem that it becomes difficult to use arises. That is, in this case, the configurations described in Patent Documents 2 and 3 cannot be employed.

さらに、例えば、自転車、電車、自動車等に乗っている状態においては、利用者自身の動作以外の加速度も加速度センサは検出するが、上記従来技術は、利用者が何らかの乗り物を利用している状態を判別できないため、乗り物を利用している状態における推定精度が劣化するという問題がある。   Further, for example, in the state of riding a bicycle, a train, an automobile, etc., the acceleration sensor detects accelerations other than the user's own motion. However, the above-described conventional technology is a state where the user uses some vehicle. Therefore, there is a problem that the estimation accuracy in the state where the vehicle is used deteriorates.

したがって、本発明は、所定方向に固定して所持する必要がなく、従来技術より推定精度が高い、運動強度及び/又は消費カロリーの推定装置を提供することを目的とする。   Therefore, an object of the present invention is to provide an apparatus for estimating exercise intensity and / or calorie consumption, which does not need to be fixed and held in a predetermined direction and has higher estimation accuracy than the prior art.

本発明における運動強度の推定装置によれば、
加速度を測定する加速度測定手段と、所定長の推定期間内に、所定長の変換期間を互いにその開始位置をずらして所定数だけ設定し、加速度測定手段が前記推定期間内に測定した加速度のうち、各変換期間に含まれる加速度をそれぞれフーリエ変換して、複数のスペクトラムを求めるフーリエ変換手段と、フーリエ変換手段が求めた複数のスペクトラムのそれぞれを、あらかじめ定めた複数の代表スペクトラムの中で、誤差の最も少ない代表スペクトラムに置換する量子化手段と、前記複数の代表スペクトラムからあらかじめ選択した複数の特徴スペクトラムそれぞれについて、量子化手段が置換して得た複数の代表スペクトラムでの出現確率を求め、求めた各特徴スペクトラムの出現確率を含む推定用データを生成する推定用データ生成手段と、各状態において1回以上測定した、前記推定用データと同じデータを含む学習データを有し、前記学習データに基づき、前記推定用データに対応する状態を判定する分類手段と、判定した状態に基づき運動強度を決定する運動強度決定手段とを有することを特徴とする。
According to the exercise intensity estimation device of the present invention,
Among the accelerations measured by the acceleration measuring means and the acceleration measuring means, the predetermined number of conversion periods are set within the estimation period of the predetermined length while the start positions of the conversion periods of the predetermined length are shifted from each other. The Fourier transform means for obtaining a plurality of spectra by Fourier transforming the accelerations included in each conversion period, and each of the plurality of spectra obtained by the Fourier transform means within a plurality of predetermined representative spectra. For each of a plurality of representative spectrums selected in advance from the plurality of representative spectrums, a quantization means for substituting the representative spectrum with the smallest representative spectrum is obtained, An estimation data generator that generates estimation data including the appearance probability of each feature spectrum And classifying means for determining a state corresponding to the estimation data based on the learning data, having learning data including the same data as the estimation data, measured at least once in each state And exercise intensity determining means for determining exercise intensity based on the above.

本発明における推定装置の他の実施形態によれば、
前記所定数だけ設定した変換期間それぞれにおいて、加速度の分散値を求める手段を備えており、前記推定用データは、前記求めた所定数ある分散値の平均値及び分散値を含んでいることも好ましい。
According to another embodiment of the estimation device of the present invention,
It is also preferable that a means for obtaining a variance value of acceleration is provided for each of the predetermined number of conversion periods, and the estimation data includes an average value and a variance value of the obtained prescribed number of variance values. .

また、本発明における推定装置の他の実施形態によれば、
前記推定用データは、量子化手段が置換して得た複数の代表スペクトラムにおいて、最大の振幅を有する周波数を示す情報を含んでいることも好ましい。
Further, according to another embodiment of the estimation device of the present invention,
The estimation data preferably includes information indicating a frequency having the maximum amplitude in a plurality of representative spectra obtained by the replacement by the quantization means.

さらに、本発明における推定装置の他の実施形態によれば、
位置情報を取得する手段を備えており、分類手段は、位置情報を状態の判定に利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
It is also preferable that a means for acquiring position information is provided, and the classification means uses the position information for determining the state.

さらに、本発明における推定装置の他の実施形態によれば、
前記位置情報は、無線通信により通信している基地局を示す情報であり、分類手段は、前記基地局の変更頻度を状態の判定に利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
The position information is information indicating a base station that is communicating by wireless communication, and the classification unit preferably uses the change frequency of the base station to determine the state.

さらに、本発明における推定装置の他の実施形態によれば、
運動強度決定手段は、運動強度の決定に、前記分散値の平均値を利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
It is also preferable that the exercise intensity determining means uses an average value of the dispersion values for determining the exercise intensity.

さらに、本発明における推定装置の他の実施形態によれば、
前記装置に対する操作についての操作情報を収集する手段を備えており、運動強度決定手段は、前記装置に対する操作が行われた場合、前記加速度測定手段が測定した加速度のうち、操作が行われた前後一定期間の加速度に対してフィルタを適用し、フィルタ後の加速度から、前記分散値の平均値を計算して、運動強度の決定に利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
Means for collecting operation information about operations on the device, and the exercise intensity determining means, when an operation on the device is performed, out of the acceleration measured by the acceleration measuring means before and after the operation is performed; It is also preferable to apply a filter to acceleration for a certain period, calculate an average value of the variance values from the filtered acceleration, and use it for determination of exercise intensity.

さらに、本発明における推定装置の他の実施形態によれば、
分類手段が判定する状態は、電車又は自動車に乗っている状態を含み、運動強度決定手段は、状態が電車又は自動車に乗っている状態である場合、前記加速度測定手段が測定した加速度に対してフィルタを適用し、フィルタ後の加速度から、前記分散値の平均値を計算して、運動強度の決定に利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
The state determined by the classification unit includes a state where the vehicle is riding on a train or a car, and the exercise intensity determination unit is configured to detect the acceleration measured by the acceleration measuring unit when the state is a state where the vehicle is riding on a train or a car. It is also preferable to apply a filter and calculate an average value of the variance values from the filtered acceleration and use it for the determination of exercise intensity.

さらに、本発明における推定装置の他の実施形態によれば、
運動強度決定手段は、運動強度の決定に、量子化手段が置換して得た複数の代表スペクトラムの中に、最も多く出現する代表スペクトラムを示す情報を利用することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
It is also preferable that the exercise intensity determining means uses information indicating a representative spectrum that appears most frequently among a plurality of representative spectra obtained by the quantization means replacing the exercise intensity.

さらに、本発明における推定装置の他の実施形態によれば、
分類手段は、複数の分類器が多段接続されたものであり、各段の分類器は、その1つ前の段の分類器が判定した状態を、さらに分類した状態を判定することも好ましい。
Furthermore, according to another embodiment of the estimation device of the present invention,
It is also preferable that the classifying unit includes a plurality of classifiers connected in multiple stages, and the classifier at each stage determines a further classified state from the state determined by the classifier at the previous stage.

本発明における消費カロリーの推定装置によれば、
前記運動強度の推定装置が決定した運動強度に基づき、消費カロリーを計算することを特徴とする。
According to the calorie consumption estimation device of the present invention,
The calorie consumption is calculated based on the exercise intensity determined by the exercise intensity estimating device.

本発明においては、推定期間内の加速度から代表スペクトラム列を求め、あらかじめ測定した各状態における特徴スペクトラムの出現確率と、求めた代表スペクトラム列における特徴スペクトラムの出現確率により状態を判定する。したがって、ある瞬間に観測されるかもしれない特異なスペクトラムの影響を防ぐことができ、事前に測定した学習データにより正確に状態判定を行うことができる。また、本発明においては、加速度の方向についての情報を推定に利用しないため、装置を利用者に固定する必要もない。   In the present invention, the representative spectrum sequence is obtained from the acceleration within the estimation period, and the state is determined based on the appearance probability of the feature spectrum in each state measured in advance and the appearance probability of the feature spectrum in the obtained representative spectrum sequence. Therefore, it is possible to prevent the influence of a peculiar spectrum that may be observed at a certain moment, and to accurately determine the state based on the learning data measured in advance. In the present invention, since information about the direction of acceleration is not used for estimation, it is not necessary to fix the device to the user.

本発明による運動強度及び/又は消費カロリーを推定する装置の簡略化した構成図である。1 is a simplified configuration diagram of an apparatus for estimating exercise intensity and / or calorie consumption according to the present invention. FIG. データ変換部の構成図である。It is a block diagram of a data converter. データ収集方法を説明する図である。It is a figure explaining a data collection method. 測定したスペクトラム系列を示す図である。It is a figure which shows the measured spectrum series. 代表スペクトラム・テーブルを示す図である。It is a figure which shows a representative spectrum table. 事前データを示す図である。It is a figure which shows prior data. 特徴スペクトラムの決定を説明する図である。It is a figure explaining the determination of a characteristic spectrum. 特徴スペクトラム・テーブルを示す図である。It is a figure which shows the characteristic spectrum table. 学習データを示す図である。It is a figure which shows learning data. 運動強度の決定を説明する図である。It is a figure explaining determination of exercise intensity. 本発明による推定装置の他の実施形態を示す図である。It is a figure which shows other embodiment of the estimation apparatus by this invention.

本発明を実施するための形態について、以下では図面を用いて詳細に説明する。   EMBODIMENT OF THE INVENTION The form for implementing this invention is demonstrated in detail below using drawing.

本発明による推定装置は、例えば、当該推定装置を携帯している利用者の状態を、“停止状態”、“歩行状態”、“走行状態”、“自転車を運転している状態”、“電車又は自動車に乗っている状態”に分類して、運動強度及び消費カロリーを推定する。図1は、本発明による推定装置の簡略化した構成図であり、推定装置は、加速度センサ1と、合成加速度計算部2と、データ変換部3と、分類器4と、運動強度決定部5と、消費カロリー計算部6とを備えている。   The estimation device according to the present invention, for example, indicates the state of a user carrying the estimation device as “stop state”, “walking state”, “running state”, “state driving a bicycle”, “train” Or, the exercise intensity and calorie consumption are estimated. FIG. 1 is a simplified configuration diagram of an estimation apparatus according to the present invention. The estimation apparatus includes an acceleration sensor 1, a composite acceleration calculation unit 2, a data conversion unit 3, a classifier 4, and an exercise intensity determination unit 5. And a calorie consumption calculation unit 6.

加速度センサ1は、本実施形態においては3軸加速度センサであり、互いに直交する3方向それぞれの加速度を周期的に測定し、測定した値を合成加速度計算部2に出力する。合成加速度計算部2は、加速度センサ1が同時に収集した3方向の加速度から、各時刻の合成加速度、つまり、実際の加速度の大きさを、以下の式により算出する。
合成加速度=(x+y+z0.5 (1)
ここで、x、y、zは、それぞ、各方向の加速度である。
The acceleration sensor 1 is a triaxial acceleration sensor in the present embodiment, periodically measures accelerations in three directions orthogonal to each other, and outputs the measured values to the combined acceleration calculation unit 2. The combined acceleration calculation unit 2 calculates the combined acceleration at each time, that is, the actual magnitude of the acceleration, from the three directions of acceleration simultaneously collected by the acceleration sensor 1 by the following formula.
Composite acceleration = (x 2 + y 2 + z 2 ) 0.5 (1)
Here, x, y, and z are accelerations in the respective directions.

データ変換部3は、合成加速度計算部2が算出した各時刻の合成加速度から、状態を推定するための推定用データを生成して分類器4に出力する。図2は、データ変換部3の構成図であり、データ変換部3は、フーリエ変換部31と、分散計算部32と、量子化部33と、代表スペクトラム・テーブル34と、推定用データ生成部35と、特徴スペクトラム・テーブル36とを備えている。   The data conversion unit 3 generates estimation data for estimating the state from the combined acceleration at each time calculated by the combined acceleration calculation unit 2 and outputs the estimation data to the classifier 4. FIG. 2 is a configuration diagram of the data conversion unit 3. The data conversion unit 3 includes a Fourier transform unit 31, a variance calculation unit 32, a quantization unit 33, a representative spectrum table 34, and an estimation data generation unit. 35 and a characteristic spectrum table 36.

図3の上段に示す様に、各加速度測定時刻における合成加速度がデータ変換部3に入力されるが、フーリエ変換部31は、複数のFFT期間、図3においては、符号60−1から符号60−LのL個のFFT期間に含まれる合成加速度を、それぞれ、フーリエ変換して、合成加速度のスペクトラムを算出する。なお、各FFT期間の開始位置は、図3に示す様に、シフト期間Tsだけずらして設定する。また、総てのFFT期間を包含する期間を、図3に示す様に、推定期間と呼ぶものとする。フーリエ変換部31は、図4に示す様に、FFT期間60−1に含まれる合成加速度から求めたスペクトラム70−1、FFT期間60−2に含まれる合成加速度から求めたスペクトラム70−2、最後に、FFT期間60−Lに含まれる合成加速度から求めたスペクトラム70−Lを含む、スペクトラムの系列を量子化部33に出力する。   As shown in the upper part of FIG. 3, the combined acceleration at each acceleration measurement time is input to the data conversion unit 3, and the Fourier transform unit 31 includes a plurality of FFT periods. The resultant acceleration included in the L FFT periods of −L is Fourier transformed to calculate the spectrum of the resultant acceleration. Note that the start position of each FFT period is set by being shifted by the shift period Ts as shown in FIG. Further, a period including all FFT periods is referred to as an estimation period as shown in FIG. As shown in FIG. 4, the Fourier transform unit 31 includes a spectrum 70-1 obtained from the combined acceleration included in the FFT period 60-1, a spectrum 70-2 calculated from the combined acceleration included in the FFT period 60-2, and the last. In addition, a spectrum series including the spectrum 70 -L obtained from the combined acceleration included in the FFT period 60 -L is output to the quantization unit 33.

量子化部33は、入力されたスペクトラム70−1〜70−Lのそれぞれについて、代表スペクトラム・テーブル34に記載されている複数の代表スペクトラムから、所定の基準において一番近いスペクトラムを選択し、選択した代表スペクトラムでスペクトラム70−1〜70−Lを置換する。図5は、代表スペクトラム・テーブル34を示す図である。なお、代表スペクトラム・テーブル34の生成については後述する。図5によると、代表スペクトラム・テーブル34は、S〜Sの識別子が付与された合計N個の代表スペクトラムを有しており、例えば、識別子としてSが付与されている代表スペクトラムは、直流成分(周波数番号0)が3であり、第1番目、第2番目、第3番目の周波数成分が、それぞれ、6、1、9であることが示されている。 The quantization unit 33 selects and selects a spectrum closest to a predetermined reference from a plurality of representative spectra described in the representative spectrum table 34 for each of the input spectra 70-1 to 70-L. The spectrums 70-1 to 70-L are replaced with the representative spectrum. FIG. 5 is a diagram showing the representative spectrum table 34. The generation of the representative spectrum table 34 will be described later. According to FIG. 5, the representative spectrum table 34 has a total of N representative spectra to which identifiers S 1 to S N are assigned. For example, the representative spectrum to which S 1 is assigned as an identifier is It is shown that the DC component (frequency number 0) is 3, and the first, second, and third frequency components are 6, 1, and 9, respectively.

量子化部33は、入力された各スペクトラム70−1〜70−Lのそれぞれを、一番近い、言い換えると、誤差が最も少ない代表スペクトラムと置換して、置換後のスペクトラム系列を出力する。ここで誤差は、例えば、対応する周波数成分の振幅差の絶対値又は振幅差の2乗を、全周波数に渡り合計した値等により評価する。例えば、スペクトラム70−1に最も近い代表スペクトラムがS34であり、スペクトラム70−2に最も近い代表スペクトラムがS184であり、スペクトラム70−Lに最も近い代表スペクトラムがS10である場合、量子化部33は、(S34、S184、・・・、S10)の様に、置換した代表スペクトラムの識別子の系列を推定用データ生成部35に出力する。なお、この置換操作は、複数のスペクトラム状態を、1つの代表スペクトラムに集約するものであり、1種の量子化である。これにより、測定したスペクトラムの系列を、代表スペクトラムの識別子の系列で表現可能となる。 The quantization unit 33 replaces each of the input spectra 70-1 to 70-L with the representative spectrum that is closest, in other words, with the least error, and outputs a spectrum sequence after replacement. Here, the error is evaluated by, for example, a value obtained by summing the absolute value of the amplitude difference of the corresponding frequency component or the square of the amplitude difference over all frequencies. For example, the closest representative spectrum to the spectrum 70-1 is S 34, if the closest representative spectrum to the spectrum 70-2 is S 184, the closest representative spectrum to a spectrum 70-L is S 10, the quantization The unit 33 outputs the replacement representative spectrum identifier series to the estimation data generation unit 35 as in (S 34 , S 184 ,..., S 10 ). In addition, this replacement operation is to aggregate a plurality of spectrum states into one representative spectrum, which is one kind of quantization. As a result, the measured spectrum series can be expressed by the representative spectrum identifier series.

また、分散計算部32は、各FFT期間に含まれる合成加速度の分散値を計算し、推定用データ生成部35に出力する。つまり、合計L個の分散値を推定用データ生成部35に出力する。   In addition, the variance calculation unit 32 calculates the variance value of the composite acceleration included in each FFT period, and outputs it to the estimation data generation unit 35. That is, a total of L variance values are output to the estimation data generation unit 35.

推定用データ生成部35は、量子化部33からのスペクトラム系列と、分散計算部32からの分散値に基づき、特徴スペクトラム・テーブル36を利用して分類器4へ出力する推定用データを生成する。ここで、特徴スペクトラムとは、代表スペクトラムの内、各状態を判別するのに重要となる代表スペクトラムであり、事前に各状態それぞれで測定した多数の測定データより、各状態において、出現確率が所定値より高いと判定した代表スペクトラムを、特徴スペクトラムとしている。なお、特徴スペクトラム決定の詳細については後述する。図8は、特徴スペクトラム・テーブル36を示す図であり、事前の測定により合計7個の代表スペクトラムが、特徴スペクトラムとして選定されている。   The estimation data generation unit 35 generates estimation data to be output to the classifier 4 using the feature spectrum table 36 based on the spectrum series from the quantization unit 33 and the variance value from the variance calculation unit 32. . Here, the characteristic spectrum is a representative spectrum that is important for distinguishing each state in the representative spectrum, and the appearance probability is predetermined in each state from a large number of measurement data measured in each state in advance. The representative spectrum determined to be higher than the value is used as the feature spectrum. Details of the characteristic spectrum determination will be described later. FIG. 8 is a diagram showing the feature spectrum table 36, and a total of seven representative spectra are selected as feature spectra by prior measurement.

図8に示す様に、7個の特徴スペクトラムが設定されている場合、本実施形態において、推定用データ生成部35は、9つのデータを含む推定用データを生成する。推定用データの第1番目のデータは、入力されたL個の代表スペクトラムにおける、第1番目の特徴スペクトラムS21の出現確率である。推定用データの第2番目のデータは、入力されたL個の代表スペクトラムにおける、第2番目の特徴スペクトラムS31の出現確率である。以下、同様に、推定用データの第7番目までのデータが計算される。 As shown in FIG. 8, when seven feature spectra are set, in this embodiment, the estimation data generation unit 35 generates estimation data including nine data. The first data of the estimated data is in the L representative spectrum input, a probability of occurrence of the first aspect spectrum S 21. Second data of the estimated data is in the L representative spectrum input, a probability of occurrence of the second feature spectrum S 31. Hereinafter, similarly, data up to the seventh of the estimation data is calculated.

さらに、推定用データの第8番目のデータは、分散計算部32からのL個の分散値の平均値(以後、V平均値と呼ぶ。)であり、最後のデータは、分散計算部32からのL個の分散値の分散値(以後、V分散値と呼ぶ。)である。つまり、例えば、量子化部33からのL個の代表スペクトラムにおいて、S21の出現確率が0.3、S31の出現確率が0.2、S10の出現確率が0.1、S66の出現確率が0、S72の出現確率が0、S33の出現確率が0、Sの出現確率が0であり、分散計算部32からのL個の分散値の平均値及び分散値が、6及び4である場合、推定用データ生成35は、(0.3,0.2,0.1,0,0,0,0,6,4)を出力する。なお、各FFT期間内の合成加速度の平均ではなく、合成加速度の分散値を求め、この分散値の推定期間に渡る平均値及び分散値を使用するのは、重力加速度成分を除去するためであり、本願発明において、合成加速度の分散値は、合成加速度の強さを表す値として利用している。つまり、V平均値は、推定期間内の合成加速度の強さの平均を示す値である。 Further, the eighth data of the estimation data is an average value (hereinafter referred to as V average value) of L variance values from the variance calculation unit 32, and the last data is from the variance calculation unit 32. Of the L variance values (hereinafter referred to as V variance values). That is, for example, in the L representative spectrum from the quantizing unit 33, the probability of occurrence of S 21 is 0.3, the probability of occurrence of S 31 is 0.2, the appearance probability of 0.1 of S 10, the S 66 appearance probability is the appearance probability 0 probability of occurrence is 0, S 9 occurrence probability is 0, S 33 of 0, S 72, the mean value and variance values of L variance from the variance calculator 32, In the case of 6 and 4, the estimation data generation 35 outputs (0.3, 0.2, 0.1, 0, 0, 0, 0, 6, 4). The reason why the dispersion value of the composite acceleration is obtained instead of the average of the composite acceleration within each FFT period and the average value and the dispersion value over the estimation period of the dispersion value is used is to remove the gravitational acceleration component. In the present invention, the dispersion value of the composite acceleration is used as a value representing the strength of the composite acceleration. That is, the V average value is a value indicating the average of the strength of the combined acceleration within the estimation period.

分類器4は、例えば、単純ベイズ分類器や、ナイブ・ベイズ・ツリー分類器であり、あらかじめ作成され、分類器4に設定されている各状態に対する1つ以上の学習データに基づき、データ変換部3から入力される推定用データに対応する状態を判定する。図9は、分類器4に設定されている学習データを示す図である。なお、図9においては、簡単のため、状態を“A状態”、“B状態”、“C状態”の3つとしている。図9においては、第1番目の学習データは、A状態において測定されたものであり、図8に示す7つの特徴スペクトラムそれぞれの出現確率と、V平均値及びV分散値が示されている。なお、学習データの生成については後述する。   The classifier 4 is, for example, a naive Bayes classifier or a naive Bayes tree classifier, and is based on one or more learning data for each state created in advance and set in the classifier 4. A state corresponding to the estimation data input from 3 is determined. FIG. 9 is a diagram showing learning data set in the classifier 4. In FIG. 9, for simplicity, there are three states, “A state”, “B state”, and “C state”. In FIG. 9, the first learning data is measured in the A state, and the appearance probability, V average value, and V variance value of each of the seven feature spectra shown in FIG. 8 are shown. The generation of learning data will be described later.

分類器4は、入力された推定用データと各学習データを比較し、推定用データとの尤度が最も高い学習データを選択して、選択した学習データに対応する状態を運動強度決定部5に出力する。また、運動強度決定部5には、データ変換部3から推定用データの生成に使用したL個の代表スペクトラムとV平均値が入力され、合成加速度計算部2から、同じく推定用データの生成に使用した合成加速度が入力される。   The classifier 4 compares the input estimation data with each learning data, selects the learning data having the highest likelihood of the estimation data, and sets the state corresponding to the selected learning data to the exercise intensity determination unit 5. Output to. The exercise intensity determination unit 5 receives L representative spectra and V average values used for generating the estimation data from the data conversion unit 3, and the combined acceleration calculation unit 2 generates the estimation data. Used synthetic acceleration is input.

運動強度決定部5は、各状態とV平均値に基づきメッツ値を決定する。また、状態が歩行及び走行である場合には、単位時間当たりの歩数によりメッツ値が異なることから、運動強度決定部5は、合成加速度に基づき単位時間当たりの歩数を決定してメッツ値の決定に利用する。なお、合成加速度から歩数を決定する方法は、公知の種々の方法を利用できる。また、例えば、同じ状態であっても、出現するスペクトラムが異なることは状態の差異があることを意味し、当然メッツ値も異なるものとなる。このため、推定用データの元となったL個の代表スペクトラムの内、最も多く出現している代表スペクトラムを主要スペクトラムとし、主要スペクトラムもメッツ値の決定に使用することが好ましい。この場合、状態及びV平均値が同じであったとしても、主要スペクトラムが異なれば、異なるメッツ値が出力されることになる。図10に、運動強度決定部5がメッツ値を決定するために使用するテーブルの例を示す。なお、図10において、a0〜a2及びb0〜b2は、所定の係数である。なお、メッツ値決定に利用する各値及び/又は式と、メッツ値との関係は、事前に行う多数の測定により求めたものである。   The exercise intensity determination unit 5 determines a Mets value based on each state and the V average value. Further, when the state is walking and running, the Met's value varies depending on the number of steps per unit time. Therefore, the exercise intensity determination unit 5 determines the number of steps per unit time based on the resultant acceleration to determine the Met value. To use. Various known methods can be used for determining the number of steps from the combined acceleration. Further, for example, even in the same state, the appearance of different spectra means that there is a difference in state, and naturally the Mets value also differs. For this reason, it is preferable that the representative spectrum that appears most frequently among the L representative spectrums that are the sources of the estimation data is the main spectrum, and the main spectrum is also used to determine the Mets value. In this case, even if the state and the V average value are the same, if the main spectrum is different, different Mets values are output. FIG. 10 shows an example of a table used by the exercise intensity determination unit 5 to determine the Mets value. In FIG. 10, a0 to a2 and b0 to b2 are predetermined coefficients. It should be noted that the relationship between each value and / or formula used to determine the Met's value and the Met's value is obtained by a number of measurements performed in advance.

さらに、運動強度決定部5は、状態が、“電車又は自動車に乗っている状態”である場合、データ変換部3から受信するV平均値を使用するのではなく、合成加速度計算部2から入力される当該推定期間の合成加速度をフィルタ処理し、これにより、合成加速度に含まれている乗り物からの振動等により生じた加速度成分を除去し、フィルタ後の合成加速度からV平均値を求めて、これをメッツ値への変換に利用することが好ましい。   Furthermore, when the state is “a state of riding a train or a car”, the exercise intensity determination unit 5 does not use the V average value received from the data conversion unit 3 but inputs from the combined acceleration calculation unit 2 The combined acceleration of the estimated period is filtered, thereby removing an acceleration component caused by vibrations from the vehicle included in the combined acceleration, and obtaining a V average value from the combined acceleration after filtering, It is preferable to use this for conversion to a Mets value.

さらに、本発明による推定装置が、移動通信端末の様な、他の機能を有する装置に実装される場合、他の機能に対する操作、例えば、移動通信端末の開閉や、ボタンの押下等による加速度成分が生じ、メッツ値の判定に影響が出る。このため、利用者の操作情報を取得する機能を設けておき、利用者が何らかの操作を行った場合、運動強度決定部5は、その前後一定期間の合成加速度をフィルタ処理して、操作による加速度成分を除去し、フィルタ後の合成加速度からV平均値を求めて、これをメッツ値への変換に利用することが好ましい。   Furthermore, when the estimation device according to the present invention is mounted on a device having other functions, such as a mobile communication terminal, an acceleration component due to an operation on the other function, for example, opening / closing of the mobile communication terminal, button pressing, etc. This affects the determination of Mets values. For this reason, a function for acquiring user operation information is provided, and when the user performs some operation, the exercise intensity determination unit 5 filters the combined acceleration for a certain period of time before and after the operation to determine the acceleration due to the operation. It is preferable to remove the component, obtain the V average value from the combined acceleration after filtering, and use this for conversion to a Mets value.

最後に、消費カロリー計算部6は、運動強度決定部5が決定したメッツ値と、あらかじめ利用者が入力した利用者の体重と、そのメッツ値の継続時間を乗じて消費カロリーを求める。   Finally, the calorie consumption calculating unit 6 obtains calorie consumption by multiplying the Met's value determined by the exercise intensity determining unit 5, the weight of the user previously input by the user, and the duration of the Met value.

続いて、代表スペクトラム・テーブル34の作成について述べる。代表スペクトラム・テーブル34は、量子化ステップに相当するものであり、代表スペクトラムの数Nは、スペクトラムの周波数の数(図5におけるK)と、各周波数成分において取り得る振幅値から、適当な数を選択する。なお、各代表スペクトラムが似通ったものとならない様に、代表スペクトラム間の誤差の総和ができるだけ大きな値となる様に選択することが好ましい。また、以下に説明する学習データの作成の際には、各状態において多数の測定を行うが、このとき取得した各状態における多数のスペクトラムから遺伝的アルゴリズムにより代表スペクトラムを決定しても良い。   Next, creation of the representative spectrum table 34 will be described. The representative spectrum table 34 corresponds to a quantization step, and the number N of representative spectra is an appropriate number based on the number of spectrum frequencies (K in FIG. 5) and the amplitude values that can be taken for each frequency component. Select. It should be noted that the total sum of errors between representative spectra is preferably selected to be as large as possible so that the representative spectra are not similar. In creating learning data described below, a large number of measurements are performed in each state. A representative spectrum may be determined by a genetic algorithm from a large number of spectra in each state acquired at this time.

続いて、特徴スペクトラム・テーブル36及び分類器4に設定する学習データの作成について説明する。学習データを作成するために、まず、測定者は、本発明による推定装置と同様の大きさ、重さ及び形状であり、かつ、3軸加速度センサを含む装置を携帯し、所定の状態にて実際の推定の時と同様に各方向の加速度を収集し、収集した各方向の加速度から式(1)により合成加速度を求める。そして、合成加速度から、実際の推定の時と同様にL個のスペクトラムの系列を求め、これを、L個の代表スペクトラムの系列に置換する。さらに、実際の推定の時と同様に、V平均値及びV分散値を計算し、L個の代表スペクトラムの系列と、V平均値及びV分散値を、当該状態に対する事前データの1つとして記録する。なお上記事前データの取得については、各状態それぞれに対して複数回実施して、各状態について、それぞれ複数の事前データを取得する。なお、このとき、同一の測定者により測定するのではなく、複数の測定者が、各状態について、それぞれ、複数回実施することが好ましい。また、装置の所持方法も種々の形態にて測定し、測定の環境、例えば、歩行や走行の場合における道路状態、乗り物に乗っている場合における乗り物種別についても、種々の状態にて測定する。   Next, creation of learning data to be set in the feature spectrum table 36 and the classifier 4 will be described. To create the learning data, first, the measurer carries a device having the same size, weight, and shape as the estimation device according to the present invention and including a three-axis acceleration sensor, and in a predetermined state. The acceleration in each direction is collected in the same manner as in the actual estimation, and the resultant acceleration is obtained from Equation (1) from the collected acceleration in each direction. Then, the L spectrum series is obtained from the combined acceleration in the same manner as the actual estimation, and this is replaced with the L representative spectrum series. Further, as in the actual estimation, the V average value and the V variance value are calculated, and the L representative spectrum series, the V average value and the V variance value are recorded as one of the prior data for the state. To do. The acquisition of the prior data is performed a plurality of times for each state, and a plurality of prior data is acquired for each state. At this time, it is preferable that a plurality of measurers carry out the measurement several times for each state, instead of measuring by the same measurer. In addition, the method of possessing the apparatus is also measured in various forms, and the measurement environment, for example, the road state in the case of walking or running, and the vehicle type in the case of riding the vehicle are also measured in various states.

図6は、事前データを示す図である。なお、図6においては、簡単のため、状態を“A状態”、“B状態”、“C状態”の3つとしている。図6によると、1番目の学習データは、A状態において測定されたものであり、L個のスペクトラク系列が、代表スペクトラムS21、S10、S83、・・・、S91であることと、V平均値及びV分散値が、それぞれ、MとDであることが示されている。 FIG. 6 is a diagram showing prior data. In FIG. 6, for simplicity, there are three states, “A state”, “B state”, and “C state”. According to FIG. 6, the first learning data is measured in the A state, and the L spectral sequences are representative spectrums S 21 , S 10 , S 83 ,..., S 91. It is shown that the V average value and the V variance value are M 1 and D 1 , respectively.

続いて、取得した事前データから、特徴スペクトラムを決定する。図7は、特徴スペクトラムの決定を説明する図である。まず、図6の事前データの各状態それぞれについて、各スペクトラムの出現確率を計算する。具体的には、例えば、A状態の事前データが150個あるとすると、1つの事前データはL個のスペクトラムを含むため、A状態全体では、合計150L個の代表スペクトラムを含むことになる。例えば、この150L個のスペクトラムに含まれている代表スペクトラムS26の数が15個であるとすると、S26の出現確率は15/(150L)である。そして、図7に示す様に、各状態別に、出現確率の大きいものから順に並べ、出現確率が所定の閾値以上である代表スペクトラムを、特徴スペクトラムに決定する。例えば、図7において閾値を0.1とすると、図8に示す様に、S21、S31、S10、S66、S72、S33、Sの7つの代表スペクトラムが特徴スペクトラムとして選択される。なお、各状態においてほぼ同じ出現確率となる代表スペクトラムは、各状態の判定に役立つものではないため、閾値以上の出現確率であったとしても特徴スペクトラムから除外できる。 Subsequently, a characteristic spectrum is determined from the acquired prior data. FIG. 7 is a diagram for explaining the determination of the characteristic spectrum. First, the appearance probability of each spectrum is calculated for each state of the prior data in FIG. Specifically, for example, if there are 150 pieces of prior data in the A state, since one piece of prior data includes L spectrums, the entire A state includes a total of 150 L representative spectra. For example, if the number of representative spectra S 26 included in the 150 L spectrums is 15, the appearance probability of S 26 is 15 / (150 L). Then, as shown in FIG. 7, the representative spectrums having the appearance probabilities equal to or higher than a predetermined threshold are determined as the feature spectrums in order of increasing appearance probability for each state. For example, when the threshold value is 0.1 in FIG. 7, as shown in FIG. 8, seven representative spectra S 21 , S 31 , S 10 , S 66 , S 72 , S 33 , and S 9 are selected as feature spectra. Is done. Note that a representative spectrum having almost the same appearance probability in each state is not useful for determination of each state, so even if it has an appearance probability equal to or higher than a threshold, it can be excluded from the feature spectrum.

また、以下の式で定義される各代表スペクトラムの情報量利得を基準として特徴スペクトラムを選択することもできる。   Also, the characteristic spectrum can be selected based on the information amount gain of each representative spectrum defined by the following equation.

Figure 2011045524
ここで、cは分類すべき各状態、Cは全状態、Pr(S)はSの出現確率、Pr(S|c)は状態cのときにおけるSの出現確率、Pr(c|S)はSが出現したときに状態cである確率を示している。具体的には、Pr(S|C)>0.01を満たし、かつ、InfoGain(S)の上位所定数の代表スペクトラムを特徴スペクトラムとする。
Figure 2011045524
Here, c is each state to be classified, C is all states, Pr (S) is the appearance probability of S, Pr (S | c) is the appearance probability of S in state c, and Pr (c | S) is The probability of being in state c when S appears is shown. Specifically, Pr (S | C)> 0.01 is satisfied, and the upper-level predetermined number of representative spectrums of InfoGain (S) are set as the characteristic spectrum.

最後に、各事前データに含まれる代表スペクトラムでの、特徴スペクトラムの出現確率を求め、図9に示す学習データを作成する。つまり、図6の2番目の事前データにおいて、特徴スペクトラムS21、S31、S10、S66、S72、S33、Sそれぞれの出現確率が、0.2、0.19、0.12、0、0、0、0.01であり、この2番目の事前データから、図9に示す2番目の学習データが生成される。 Finally, the appearance probability of the characteristic spectrum in the representative spectrum included in each prior data is obtained, and learning data shown in FIG. 9 is created. That is, in the second prior data in FIG. 6, the appearance probabilities of the feature spectra S 21 , S 31 , S 10 , S 66 , S 72 , S 33 , S 9 are 0.2, 0.19, 0. The second learning data shown in FIG. 9 is generated from the second prior data.

以上、本発明においては、加速度の分散値に加えて、推定期間内の各スペクトラムの出現頻度を状態判定に利用する。したがって、ある瞬間に観測されるかもしれない特異なスペクトラムの影響を防ぐことができ、事前に測定した学習データにより正確に状態判定を行うことができる。また、本発明においては、3軸加速度センサの各軸の方向についての情報を推定に利用しないため、装置を利用者に固定する必要もない。   As described above, in the present invention, in addition to the acceleration dispersion value, the frequency of appearance of each spectrum within the estimation period is used for state determination. Therefore, it is possible to prevent the influence of a peculiar spectrum that may be observed at a certain moment, and to accurately determine the state based on the learning data measured in advance. In the present invention, since information about the direction of each axis of the triaxial acceleration sensor is not used for estimation, it is not necessary to fix the device to the user.

なお、上記実施形態においては、特徴スペクトラムの出現確率と、加速度の分散に基づく値を推定用データとしていたが、量子化部33が出力するL個の代表スペクトラムの内、振幅値が最大の周波数を示す情報を推定用データに追加しても良い。これは、最も強い加速度を有する周波数成分は、状態により異なるからであり、さらに、推定精度を高くすることができる。なお、このとき当然であるが、学習データにも振幅値が最大の周波数を追加することになる。   In the above embodiment, the value based on the appearance probability of the characteristic spectrum and the variance of the acceleration is used as the estimation data, but the frequency having the maximum amplitude value among the L representative spectra output from the quantization unit 33. May be added to the estimation data. This is because the frequency component having the strongest acceleration varies depending on the state, and the estimation accuracy can be further increased. Of course, at this time, the frequency having the maximum amplitude value is also added to the learning data.

また、本発明の推定装置が、移動通信端末に実装される場合、移動通信端末が周期的に取得する基地局を特定する情報を、分類器4に供給し、分類器4がこの情報を状態の判定に利用する形態であっても良い。具体的には、接続している基地局変更の頻度を、その頻度が多い程、大きくなる係数に変換し、この係数で“電車又は自動車に乗っている状態”の学習データに対する尤度が大きくなる様に調整する。“停止状態”と“電車又は自動車に乗っている状態”については、特徴スペクトラムの出現確率として互いに似た値となる場合があり、基地局変更の頻度を考慮することで、より正確に“停止状態”と“電車又は自動車に乗っている状態”を区別することができる。なお、移動通信端末がGPSを利用して位置情報を取得する機能を備えているのであれば、基地局変更の頻度に代えて、取得した位置情報から求めた速度を利用しても良い。   Further, when the estimation apparatus of the present invention is implemented in a mobile communication terminal, information specifying a base station periodically acquired by the mobile communication terminal is supplied to the classifier 4, and the classifier 4 It may be a form used for the determination. Specifically, the frequency of changing the connected base station is converted to a coefficient that increases as the frequency increases, and the likelihood of the learning data “in a train or a car” increases with this coefficient. Adjust as follows. The “stopped state” and “the state of riding in a train or car” may be similar to each other as the appearance probability of the characteristic spectrum, and the “stopped” can be more accurately considered by considering the frequency of base station changes. A “state” can be distinguished from a “state in a train or a car”. Note that if the mobile communication terminal has a function of acquiring position information using GPS, the speed obtained from the acquired position information may be used instead of the frequency of base station change.

さらに、分類器4を多段構成とし、より詳しい分類を行うことも可能である。具体的には、図11に示す様に、分類器4と運動強度決定部5の間に4つの分類器41〜44を追加する。本実施形態において、分類器41は、歩行状態の詳細状態を判定する分類器であり、分類器42は、走行状態の詳細状態を判定する分類器であり、分類器43は、自転車を運転している状態の詳細状態を判定する分類器であり、分類器44は、電車又は自動車に乗っている状態の詳細状態を判定する分類器である。分類器4は、推定用データから歩行状態と判定した場合には分類器41に推定用データを出力し、走行状態と判定した場合には分類器42に推定用データを出力し、自転車を運転している状態と判定した場合には分類器43に推定用データを出力し、電車又は自動車に乗っている状態と判定した場合には分類器44に推定用データを出力し、停止状態と判定した場合には、運動強度決定部5に停止状態との結果を出力する。   Further, the classifier 4 can have a multi-stage configuration to perform more detailed classification. Specifically, as shown in FIG. 11, four classifiers 41 to 44 are added between the classifier 4 and the exercise intensity determination unit 5. In this embodiment, the classifier 41 is a classifier that determines the detailed state of the walking state, the classifier 42 is a classifier that determines the detailed state of the running state, and the classifier 43 drives the bicycle. The classifier 44 is a classifier that determines the detailed state of a state in which the vehicle is on a train or a car. The classifier 4 outputs the estimation data to the classifier 41 when it is determined from the estimation data as the walking state, and outputs the estimation data to the classifier 42 when it is determined as the running state, and drives the bicycle. When it is determined that the vehicle is in a state of being in a train, the estimation data is output to the classifier 43. If it is, the result of the stop state is output to the exercise intensity determination unit 5.

分類器41は、歩行状態を、“通常”、“坂道上り”、“坂道下り”、“階段上り”、“階段下り”のいずれかに分類し、分類器42は、走行状態を、“通常”、“坂道上り”、“坂道下り”、“階段上り”、“階段下り”のいずれかに分類し、分類器43は、自転車を運転している状態を、“通常”、“坂道上り”、“坂道下り”のいずれかに分類し、分類器44は、電車又は自動車に乗っている状態を、“座位”、“立位”
のいずれかに分類し、運動強度決定部5に決定した詳細状態を出力する。運動強度決定部5は、詳細状態に基づきメッツ値を決定するため、例えば、歩行状態であり歩数が同じであっても、通常状態、つまり、平坦な道路を歩行している状態と、階段上り状態では、異なるメッツ値が算出される。
The classifier 41 classifies the walking state as one of “normal”, “uphill”, “downhill”, “upstairs”, and “downstairs”, and the classifier 42 sets the running state to “normal”. ”,“ Uphill ”,“ downhill ”,“ upstairs ”,“ downstairs ”, and the classifier 43 indicates that the state of driving a bicycle is“ normal ”,“ uphill ” , “Slope descent”, and the classifier 44 determines whether the train or car is in the “sitting position” or “standing position”.
And the detailed state determined by the exercise intensity determination unit 5 is output. Since the exercise intensity determination unit 5 determines the Met's value based on the detailed state, for example, even when the walking state is the same as the number of steps, the normal state, that is, the state walking on a flat road, In the state, different Mets values are calculated.

なお、各分類器41から44には、分類器4と同様に、事前に測定した各詳細状態の学習データがあらかじめ設定されている。分類器を多段接続構成にし、最初に大きく異なる状態を判定し、その後、各状態をさらに細かく分類した状態に対する学習データに基づき、これら詳細状態の判定を行うことで、詳細状態を精度高く推定することが可能になる。   As in the classifier 4, learning data for each detailed state measured in advance is set in advance in each of the classifiers 41 to 44. A classifier is configured in a multi-stage connection configuration, and a state that is greatly different first is determined, and then the detailed state is estimated based on the learning data for the state in which each state is further finely classified, thereby accurately estimating the detailed state. It becomes possible.

1 加速度センサ
2 合成加速度計算部
3 データ変換部
4、41、42、43、44 分類器
5 運動強度決定部
6 消費カロリー計算部
31 フーリエ変換部
32 分散計算部
33 量子化部
34 代表スペクトラム・テーブル
35 推定用データ生成部
36 特徴スペクトラム・テーブル
60−1〜60−L FFT期間
70−1〜70−L スペクトラム
DESCRIPTION OF SYMBOLS 1 Acceleration sensor 2 Synthetic acceleration calculation part 3 Data conversion part 4, 41, 42, 43, 44 Classifier 5 Exercise intensity determination part 6 Calorie consumption calculation part 31 Fourier transform part 32 Variance calculation part 33 Quantization part 34 Representative spectrum table 35 Estimation Data Generation Unit 36 Feature Spectrum Table 60-1 to 60-L FFT Period 70-1 to 70-L Spectrum

Claims (11)

加速度を測定する加速度測定手段と、
所定長の推定期間内に、所定長の変換期間を互いにその開始位置をずらして所定数だけ設定し、加速度測定手段が前記推定期間内に測定した加速度のうち、各変換期間に含まれる加速度をそれぞれフーリエ変換して、複数のスペクトラムを求めるフーリエ変換手段と、
フーリエ変換手段が求めた複数のスペクトラムのそれぞれを、あらかじめ定めた複数の代表スペクトラムの中で、誤差の最も少ない代表スペクトラムに置換する量子化手段と、
前記複数の代表スペクトラムからあらかじめ選択した複数の特徴スペクトラムそれぞれについて、量子化手段が置換して得た複数の代表スペクトラムでの出現確率を求め、求めた各特徴スペクトラムの出現確率を含む推定用データを生成する推定用データ生成手段と、
各状態において1回以上測定した、前記推定用データと同じデータを含む学習データを有し、前記学習データに基づき、前記推定用データに対応する状態を判定する分類手段と、
判定した状態に基づき運動強度を決定する運動強度決定手段と、
を備えている装置。
Acceleration measuring means for measuring acceleration;
Within a predetermined length of the estimation period, a predetermined number of conversion periods of a predetermined length are set with the start positions being shifted from each other, and the acceleration included in each conversion period among the accelerations measured by the acceleration measuring means within the estimation period is calculated. Fourier transform means for obtaining a plurality of spectra by performing Fourier transform respectively,
Quantization means for replacing each of the plurality of spectra obtained by the Fourier transform means with a representative spectrum with the least error among a plurality of predetermined representative spectra,
For each of a plurality of feature spectra selected in advance from the plurality of representative spectra, the appearance probability in the plurality of representative spectra obtained by the replacement by the quantizing means is obtained, and the estimation data including the appearance probability of each obtained feature spectrum is obtained. A data generation means for estimation to be generated;
Classification means for determining the state corresponding to the estimation data based on the learning data, having learning data including the same data as the estimation data, measured at least once in each state;
Exercise intensity determining means for determining exercise intensity based on the determined state;
A device equipped with.
前記所定数だけ設定した変換期間それぞれにおいて、加速度の分散値を求める手段を備えており、
前記推定用データは、前記求めた所定数ある分散値の平均値及び分散値を含んでいる、
請求項1に記載の装置。
In each of the conversion periods set by the predetermined number, a means for obtaining a dispersion value of acceleration is provided.
The estimation data includes an average value and a variance value of the determined predetermined number of variance values,
The apparatus of claim 1.
前記推定用データは、量子化手段が置換して得た複数の代表スペクトラムにおいて、最大の振幅を有する周波数を示す情報を含んでいる、
請求項1又は2に記載の装置。
The estimation data includes information indicating a frequency having the maximum amplitude in a plurality of representative spectra obtained by the replacement by the quantization means.
The apparatus according to claim 1 or 2.
位置情報を取得する手段を備えており、
分類手段は、位置情報を状態の判定に利用する、
請求項1から3のいずれか1項に記載の装置。
A means for acquiring location information;
The classifying means uses the position information to determine the state.
The apparatus according to any one of claims 1 to 3.
前記位置情報は、無線通信により通信している基地局を示す情報であり、
分類手段は、前記基地局の変更頻度を状態の判定に利用する、
請求項4に記載の装置。
The location information is information indicating a base station communicating by wireless communication,
The classifying means uses the change frequency of the base station to determine the state.
The apparatus according to claim 4.
運動強度決定手段は、運動強度の決定に、前記分散値の平均値を利用する、
請求項2から5のいずれか1項に記載の装置。
The exercise intensity determining means uses an average value of the variance values for determining the exercise intensity.
The device according to any one of claims 2 to 5.
前記装置に対する操作についての操作情報を収集する手段を備えており、
運動強度決定手段は、前記装置に対する操作が行われた場合、前記加速度測定手段が測定した加速度のうち、操作が行われた前後一定期間の加速度に対してフィルタを適用し、フィルタ後の加速度から、前記分散値の平均値を計算して、運動強度の決定に利用する、
請求項6に記載の装置。
Comprising means for collecting operation information on operations on the device;
The exercise intensity determining means applies a filter to the acceleration for a certain period before and after the operation is performed among the accelerations measured by the acceleration measuring means when an operation is performed on the device, and from the acceleration after the filter , Calculating an average value of the variance values and using it to determine exercise intensity;
The apparatus according to claim 6.
分類手段が判定する状態は、電車又は自動車に乗っている状態を含み、
運動強度決定手段は、状態が電車又は自動車に乗っている状態である場合、前記加速度測定手段が測定した加速度に対してフィルタを適用し、フィルタ後の加速度から、前記分散値の平均値を計算して、運動強度の決定に利用する、
請求項6又は7に記載の装置。
The state determined by the classification means includes a state of riding a train or a car,
The exercise intensity determining means applies a filter to the acceleration measured by the acceleration measuring means when the state is in a train or a car, and calculates the average value of the variance values from the filtered acceleration. And used to determine exercise intensity,
Apparatus according to claim 6 or 7.
運動強度決定手段は、運動強度の決定に、量子化手段が置換して得た複数の代表スペクトラムの中に、最も多く出現する代表スペクトラムを示す情報を利用する、
請求項1から8のいずれか1項に記載の装置。
The exercise intensity determination means uses information indicating the representative spectrum that appears most frequently among the plurality of representative spectra obtained by the quantization means for determining the exercise intensity.
Apparatus according to any one of claims 1 to 8.
分類手段は、複数の分類器が多段接続されたものであり、
各段の分類器は、その1つ前の段の分類器が判定した状態を、さらに分類した状態を判定する、
請求項1から9のいずれか1項に記載の装置。
The classification means is a multi-stage connection of multiple classifiers,
The classifier of each stage determines the state further classified from the state determined by the classifier of the previous stage.
Apparatus according to any one of claims 1 to 9.
請求項1から10のいずれか1項に記載の装置が決定した運動強度に基づき、消費カロリーを計算する装置。   The apparatus which calculates calorie consumption based on the exercise intensity which the apparatus of any one of Claim 1 to 10 determined.
JP2009196503A 2009-08-27 2009-08-27 Device for estimating exercise intensity and/or calorie consumption of carrier Withdrawn JP2011045524A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009196503A JP2011045524A (en) 2009-08-27 2009-08-27 Device for estimating exercise intensity and/or calorie consumption of carrier

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009196503A JP2011045524A (en) 2009-08-27 2009-08-27 Device for estimating exercise intensity and/or calorie consumption of carrier

Publications (1)

Publication Number Publication Date
JP2011045524A true JP2011045524A (en) 2011-03-10

Family

ID=43832456

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009196503A Withdrawn JP2011045524A (en) 2009-08-27 2009-08-27 Device for estimating exercise intensity and/or calorie consumption of carrier

Country Status (1)

Country Link
JP (1) JP2011045524A (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140117545A (en) * 2012-01-18 2014-10-07 나이키 이노베이트 씨.브이. Activity points
CN104169924A (en) * 2012-01-18 2014-11-26 耐克创新有限合伙公司 Activity identification
JP2016040072A (en) * 2014-08-11 2016-03-24 日立金属株式会社 Tool failure detection method
JP2016040071A (en) * 2014-08-11 2016-03-24 日立金属株式会社 Tool failure detection method
US9383220B2 (en) 2010-11-01 2016-07-05 Nike, Inc. Activity identification
US9415266B2 (en) 2011-11-01 2016-08-16 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
JP2016150102A (en) * 2015-02-17 2016-08-22 日本電信電話株式会社 Sequential attitude identification device, autonomic nerve function information acquisition device, method, and program
US9539486B2 (en) 2010-11-01 2017-01-10 Nike, Inc. Wearable device assembly having athletic functionality
JP2017023449A (en) * 2015-07-23 2017-02-02 沖電気工業株式会社 Activity amount estimation device, activity amount estimation method, program, model generation device, and model generation method
JP2017060899A (en) * 2017-01-17 2017-03-30 京セラ株式会社 Portable apparatus and activity amount calculation system
JP2017189267A (en) * 2016-04-12 2017-10-19 株式会社アコーズ Activity meter
WO2017199305A1 (en) * 2016-05-16 2017-11-23 富士通株式会社 Information processing device, information processing method, and information processing program
JP2021048965A (en) * 2019-09-24 2021-04-01 カシオ計算機株式会社 Feature amount extraction device, state estimation device, feature amount extraction method, state estimation method and program
CN112633093A (en) * 2020-12-10 2021-04-09 中国人民解放军32181部队 Equipment state detection method and system

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10456623B2 (en) 2010-11-01 2019-10-29 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US12125575B2 (en) 2010-11-01 2024-10-22 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US12062424B2 (en) 2010-11-01 2024-08-13 Nike, Inc. Wearable device assembly having athletic functionality
US11798673B2 (en) 2010-11-01 2023-10-24 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US11749395B2 (en) 2010-11-01 2023-09-05 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US9383220B2 (en) 2010-11-01 2016-07-05 Nike, Inc. Activity identification
US11735308B2 (en) 2010-11-01 2023-08-22 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US11495341B2 (en) 2010-11-01 2022-11-08 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US9616289B2 (en) 2010-11-01 2017-04-11 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
US9757640B2 (en) 2010-11-01 2017-09-12 Nike, Inc. Wearable device assembly having athletic functionality
US9539486B2 (en) 2010-11-01 2017-01-10 Nike, Inc. Wearable device assembly having athletic functionality
US9750976B2 (en) 2010-11-01 2017-09-05 Nike, Inc. Wearable device assembly having athletic functionality and trend tracking
US9415266B2 (en) 2011-11-01 2016-08-16 Nike, Inc. Wearable device assembly having athletic functionality and milestone tracking
KR101797421B1 (en) * 2012-01-18 2017-11-13 나이키 이노베이트 씨.브이. Method and apparatus for identifying activities
JP2015509755A (en) * 2012-01-18 2015-04-02 ナイキ イノベイト シーブイ Activity points
KR101669649B1 (en) * 2012-01-18 2016-10-27 나이키 이노베이트 씨.브이. Activity points
JP2016195788A (en) * 2012-01-18 2016-11-24 ナイキ イノベイト シーブイ Activity points
US10463278B2 (en) 2012-01-18 2019-11-05 Nike, Inc. Activity and inactivity monitoring
KR20140117545A (en) * 2012-01-18 2014-10-07 나이키 이노베이트 씨.브이. Activity points
CN104169924A (en) * 2012-01-18 2014-11-26 耐克创新有限合伙公司 Activity identification
KR101870906B1 (en) * 2012-01-18 2018-07-20 나이키 이노베이트 씨.브이. Activity points
US10034624B2 (en) 2012-01-18 2018-07-31 Nike, Inc. Activity points
KR20160127142A (en) * 2012-01-18 2016-11-02 나이키 이노베이트 씨.브이. Activity points
JP2016040072A (en) * 2014-08-11 2016-03-24 日立金属株式会社 Tool failure detection method
JP2016040071A (en) * 2014-08-11 2016-03-24 日立金属株式会社 Tool failure detection method
JP2016150102A (en) * 2015-02-17 2016-08-22 日本電信電話株式会社 Sequential attitude identification device, autonomic nerve function information acquisition device, method, and program
US11998360B2 (en) 2015-02-17 2024-06-04 Nippon Telegraph And Telephone Corporation Device and method for sequential posture identification and autonomic function information acquisition
US11350879B2 (en) 2015-02-17 2022-06-07 Nippon Telegraph And Telephone Corporation Device and method for sequential posture identification and autonomic function information acquisition
JP2017023449A (en) * 2015-07-23 2017-02-02 沖電気工業株式会社 Activity amount estimation device, activity amount estimation method, program, model generation device, and model generation method
JP2017189267A (en) * 2016-04-12 2017-10-19 株式会社アコーズ Activity meter
WO2017199305A1 (en) * 2016-05-16 2017-11-23 富士通株式会社 Information processing device, information processing method, and information processing program
JPWO2017199305A1 (en) * 2016-05-16 2019-01-17 富士通株式会社 Information processing apparatus, information processing method, and information processing program
JP2017060899A (en) * 2017-01-17 2017-03-30 京セラ株式会社 Portable apparatus and activity amount calculation system
JP7205433B2 (en) 2019-09-24 2023-01-17 カシオ計算機株式会社 State estimation device, state estimation method and program
JP2021048965A (en) * 2019-09-24 2021-04-01 カシオ計算機株式会社 Feature amount extraction device, state estimation device, feature amount extraction method, state estimation method and program
CN112633093A (en) * 2020-12-10 2021-04-09 中国人民解放军32181部队 Equipment state detection method and system
CN112633093B (en) * 2020-12-10 2024-04-16 中国人民解放军32181部队 Equipment state detection method and system

Similar Documents

Publication Publication Date Title
JP2011045524A (en) Device for estimating exercise intensity and/or calorie consumption of carrier
JP5356923B2 (en) Method and system for estimating movement state of portable terminal device
JP5923593B2 (en) System, method, and apparatus for classifying user activity using combined likelihood function values of mobile devices
US9700241B2 (en) Gait analysis system and method
ES2702312T3 (en) Method and apparatus to identify transitions between sitting and standing postures
CN105678222A (en) Human behavior identification method based on mobile equipment
JP2015506883A (en) Method and apparatus for elevator motion detection
JP5946364B2 (en) Time-series data processing system and time-series data processing method
JP2011152360A (en) Walking information extraction device, walking information extraction method, and computer readable storage medium
JP2017504440A (en) Improved gait detection in user movement measurements
CN104091053B (en) Method and apparatus for automatic detection behavior pattern
KR101418333B1 (en) Apparatus and method for recognizing user activity
WO2015066348A2 (en) Method and system for estimating multiple modes of motion
WO2012135726A1 (en) Devices, methods, and apparatuses for inferring a position of a mobile device
Altini et al. Self-calibration of walking speed estimations using smartphone sensors
EP2725975A1 (en) Sit-to-stand transfer detection
WO2016180653A1 (en) Pedometer and method for analyzing motion data
JP2014512045A (en) Method, device, and apparatus for activity classification using temporal scaling of time-based features
CN111189638A (en) HMM and QPSO optimization algorithm-based bearing fault degree identification method
JP5348941B2 (en) Method and system for estimating movement state of portable terminal device
CN104036283A (en) System and method for determining the activity of a mobile element
JP6794793B2 (en) Walking state estimation method, walking state estimation program and information processing device
CN114267430A (en) Route and environment-based tourist physical ability and consumption estimation method and system
KR20140037626A (en) Step counter detection apparatus and method using accelerometer
CN107239147A (en) A kind of human body context aware method based on wearable device, apparatus and system

Legal Events

Date Code Title Description
A300 Withdrawal of application because of no request for examination

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20121106