(57)【特許請求の範囲】
【請求項１】
物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップと、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップと、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納するステップと、
を含み、コンピュータにより実行される情報処理方法。
【請求項２】
前記物体は人体であり、
前記記憶装置に格納されている前記座標値を用いて、前記人体の状態を特定する状態特定ステップ
をさらに含む、請求項１記載の情報処理方法。
【請求項３】
前記状態特定ステップが、
前記記憶装置に格納されている前記座標値により表される、前記水平面上の点列についての中心点を通る直線によって前記水平面を所定角度毎の領域に分割するステップと、
前記領域毎に、前記中心点からの距離が最大となり且つ前記座標値により表される点である第１の点を特定し、当該第１の点の座標値を前記中心点を原点とする場合の座標値に変換し、前記記憶装置に格納するステップと、
前記人体の状態の候補毎に当該人体の状態の候補に対応する、各前記領域内の点である第２の点の座標値を格納する状態データベースを用いて、前記記憶装置に格納されている前記第１の点の座標値と前記第２の点の座標値とから算出され且つ全ての前記領域に関する距離が最小となる前記人体の状態の候補を特定するステップと、
を含む、請求項２記載の情報処理方法。
【請求項４】
前記物体は人であり、
前記補正ステップが、
前記人の歩行時の標準加速度を表すテンプレートを格納するテンプレート格納部に格納されている前記テンプレートを、歩行していたとみなされる期間において計測された前記加速度の値から特定される歩行周期及び加速度の振幅に合うように変形するステップと、
変形後の前記テンプレートを用いて、歩行していたとみなされる期間において計測された前記加速度の値から前記歩行時の標準加速度を取り除くステップと、
を含む、請求項１乃至３記載の情報処理方法。
【請求項５】
前記初速度算出ステップが、
前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値ａと、前記所定期間と前記所定のサンプリング間隔Δｔとから算出されるサンプル数Ｎとを用いて、
【数１】
0005477238.tif 000034
（ｉはサンプル番号を示す）
に従って、変数Ｋの値を算出するステップと、
前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値ａと、前記変数Ｋの値と、前記所定のサンプリング間隔Δｔと、前記サンプル数Ｎとを用いて、
【数２】
0005477238.tif 000035
（ｊ及びｌはサンプル番号を表す変数である）
に従って、前記物体の初速度ｖ0の値を算出するステップと、
を含む、請求項１乃至４記載の情報処理方法。
【請求項６】
物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップと、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップと、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納するステップと、
をコンピュータに実行させるためのプログラム。
【請求項７】
所定のサンプリング間隔で計測された物体の加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正部と、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出部と、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納する重心位置算出部と、
を有する情報処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、重心動揺を分析するための技術に関する。
【背景技術】
【０００２】
人は、立位時、常に細かい揺れ（重心動揺と呼ばれる）を行うことでバランスを保っている。従来から、重心動揺計を用いて重心動揺を測定し、測定結果から体のふらつき度合いを調べるような検査が行われている。重心動揺計では、被測定者が重心動揺計の上に乗ることで、被測定者の重心動揺を測定する。
【０００３】
また、従来から、加速度センサを用いて重心動揺を推定する技術が存在している。この技術では、加速度センサによって計測された加速度値を２階積分することで、各計測時刻の時点における位置を求め、重心動揺を推定する。なお、加速度値を２階積分する際には、初期位置及び初速度といった積分定数を与えることになるが、この技術では、初期位置及び初速度を０として計算を行うようになっている。すなわち、この技術では、完全に静止した状態で計測を開始することを想定している。そのため、完全に静止していない状態で計測が開始されると、重心動揺の推定において誤差が生じてしまう。また、例えば人が動いているような場合には、重心動揺を推定することはできない。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】ＷＯ２００４／０１４２３０号公報
【特許文献２】特開２００８−２２９２６６号公報
【特許文献３】特開２００８−０７３２６７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従って、本技術の目的は、一側面において、計測開始時に物体が静止しているか否かにかかわらず、加速度センサによって計測された加速度値から物体の重心動揺を高精度に推定するための技術を提供することである。
【課題を解決するための手段】
【０００６】
本情報処理方法は、（Ａ）物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップと、（Ｂ）物体の位置を加速度及び初速度から求めるための第１の関係式並びに物体の重心動揺で物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、データ格納部に格納されている、所定期間内の加速度の補正後の値を代入して、物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップと、（Ｃ）第１の関係式に、データ格納部に格納されている加速度の補正後の値と記憶装置に格納されている初速度の値とを代入して、物体の重心位置を表す、水平面における座標値を算出し、記憶装置に格納するステップとを含む。
【発明の効果】
【０００７】
計測開始時に物体が静止しているか否かにかかわらず、加速度センサによって計測された加速度値から物体の重心動揺を高精度に推定できる。
【図面の簡単な説明】
【０００８】
【図１】図１は、第１の実施の形態に係る携帯端末の機能ブロック図である。
【図２】図２は、加速度データに対して所定の補正を行う処理（前処理）の処理フローを示す図である。
【図３】図３は、計測結果格納部に格納されるデータの一例を示す図である。
【図４】図４は、パラメータ格納部に格納されるデータの一例を示す図である。
【図５】図５は、オフセット補正前の加速度とオフセット補正後の加速度とを示す図である。
【図６】図６は、単位ベクトルとｙｚ平面のなす角度を説明するための図である。
【図７】図７は、区間データ格納部に格納されるデータの一例を示す図である。
【図８】図８は、直前の４歩分の歩行周期を取得する処理を説明するための図である。
【図９】図９は、加速度データに対して所定の補正を行う処理（前処理）の処理フローを示す図である。
【図１０】図１０は、回転行列算出処理の処理フローを示す図である。
【図１１】図１１は、進行方向（ｙ軸）とベクトル（ａy，ａz）のなす角φを説明するための図である。
【図１２】図１２は、標準加速度差引処理の処理フローを示す図である。
【図１３】図１３は、歩行時における標準加速度を表すテンプレートの一例を示す図である。
【図１４】図１４は、図１３に示したテンプレートを時間軸方向に縮めた場合の例を示す図である。
【図１５】図１５は、図１４に示したテンプレートを加速度方向に伸ばした場合の例を示す図である。
【図１６】図１６は、重心動揺及び体調を推定する処理の処理フローを示す図である。
【図１７】図１７は、データ格納部からデータを読み出す処理を模式的に示す図である。
【図１８】図１８は、変数値格納部に格納されるデータの一例を示す図である。
【図１９】図１９は、初速度格納部に格納されるデータの一例を示す図である。
【図２０】図２０は、重心動揺データ格納部に格納されるデータの一例を示す図である。
【図２１】図２１は、クラスタ分析に用いる代表データを模式的に示す図である。
【図２２】図２２は、体調推定処理の処理フローを示す図である。
【図２３】図２３は、中心変位Ｄを原点とする相対位置に変換後のデータを示す図である。
【図２４】図２４は、第１の点を特定する処理を模式的に示す図である。
【図２５】図２５は、第１の点に関するデータを示す図である。
【図２６】図２６は、体調ＤＢに格納されるデータの一例を示す図である。
【図２７】図２７は、体調推定処理の処理フローを示す図である。
【図２８】図２８は、体調履歴データ格納部に格納される体調履歴データの一例を示す図である。
【図２９】図２９は、体調履歴データ格納部に格納される体調管理データの一例を示す図である。
【図３０】図３０（ａ）乃至（ｄ）は、変調変化の例を示す図である。
【図３１】図３１は、第２の実施の形態に係るシステムの構成図である。
【図３２】図３２は、第２の実施の形態に係るサーバの機能ブロック図である。
【図３３】図３３は、本実施の形態の第１の態様に係る情報処理方法の処理フローである。
【図３４】図３４は、本実施の形態の第２の態様に係る情報処理装置の機能ブロック図である。
【図３５】図３５は、コンピュータの機能ブロック図である。
【発明を実施するための形態】
【０００９】
本技術の実施の形態では、重心動揺では真の中心位置に戻ろうとする力が働くという特性に着目して、バネモデルを適用することで、高精度に重心動揺を推定する。具体的には、
・重心動揺で身体に働く力（Ｆ）は、真の中心点（Ｃ）に向かう
・重心動揺で身体に働く力（Ｆ）は、真の中心点（Ｃ）からの距離に比例する
という特性を用いて、重心動揺で身体に働く力（Ｆ）をバネモデルの運動方程式で表すと、以下のような関係式で表される。
Ｆ＝ｍ＊ａt＝−ｋ＊（ｘt−Ｃ）          （１）
ここで、ｍは質量を表し、ａtは加速度を表し、ｋはバネモデルにおけるバネ定数を表し、ｘtは重心位置を表す。なお、ｔは０からＮ−１までの整数を表し、Ｎは加速度のサンプル数を表す。
【００１０】
また、加速度を２階積分することで位置が計算できる。ここで、加速度ａt及び位置ｘtは離散値を取るため、加速度と位置の関係は、以下のような式で表される。
【００１１】
【数１】
0005477238.tif 000002
ここで、Δｔはサンプリング間隔を表し、ｖ0は初速度を表し、ｘ0は初期位置を表す。
【００１２】
また、（２）式の添え字をｔ−１に変更すると、以下のような式になる。
【００１３】
【数２】
0005477238.tif 000003
【００１４】
そして、（２）及び（３）式を用いて、ｘt−ｘt-1を表すと、以下のような式で表される。
【００１５】
【数３】
0005477238.tif 000004
【００１６】
一方、（１）式に自己回帰モデル（ＡＲモデル：autoregression model）を適用すると、以下のように変形される。
ａt＝−ｋ／ｍ＊ｘt-1＋ｗt      （５）
ここで、ｗtは誤差項を表す。
【００１７】
また、（５）式の添え字をｔ−１に変更すると、以下のような式になる。
ａt-1＝−ｋ／ｍ＊ｘt-2＋ｗt-1      （６）
【００１８】
そして、（５）及び（６）式を用いて、ａt−ａt-1を表すと、以下のような式で表される。
ａt−ａt-1＝−ｋ／ｍ＊（ｘt-1−ｘt-2）＋ｗt−ｗt-1      （７）
【００１９】
そして、（４）式を（７）式に代入すると、以下のような式になる。
【００２０】
【数４】
0005477238.tif 000005
【００２１】
また、（８）式は、以下のように変形される。
【００２２】
【数５】
0005477238.tif 000006
【００２３】
【数６】
0005477238.tif 000007
【００２４】
また、（９）式の添え字をｔ−１に変更すると、以下のような式になる。
【００２５】
【数７】
0005477238.tif 000008
【００２６】
そして、（９）及び（１０）式を用いて、ａt−ａt-1を表すと、以下のような式で表される。
【００２７】
【数８】
0005477238.tif 000009
【００２８】
また、（１１）式は、以下のように変形される。
ａt＝（２−ｋ／ｍ＊Δｔ2）−ａt-2＋ｗt−2ｗt-1＋ｗt-2      （１２）
【００２９】
ここで、変数Ｋ＝ｋ／ｍ＊Δｔ2とし、誤差項をＷt＝ｗt−2ｗt-1＋ｗt-2とすると、（１２）式は、以下のように変形される。
ａt＝（２−Ｋ）＊ａt-1−ａt-2＋Ｗt      （１３）
【００３０】
そして、（１３）式に最小二乗法を適用すると、以下のような式が導出される。
【００３１】
【数９】
0005477238.tif 000010
【００３２】
また、（１４）式は、以下のように変形される。
【００３３】
【数１０】
0005477238.tif 000011
【００３４】
そして、変数Ａ0、Ａ1、Ａ2及びγを以下のように定義する。
【００３５】
【数１１】
0005477238.tif 000012
【００３６】
そして、（１５）式を変数Ａ0、Ａ1、Ａ2及びγを用いて表すと、以下のような式で表される。
Ａ0γ2−４Ａ1γ＋２Ａ0＋２Ａ2      （１６）
【００３７】
そして、（１６）式をγで偏微分すると、２Ａ0γ−４Ａ1となり、ここでは、２Ａ0γ−４Ａ1＝０とする。そうすると、γ＝２＊Ａ1／Ａ0となり、変数Ｋは、以下のように定義される。
【００３８】
【数１２】
0005477238.tif 000013
【００３９】
一方で、（５）式に（２）式を代入すると、以下のような式になる。
【００４０】
【数１３】
0005477238.tif 000014
【００４１】
また、（１８）式の添え字をｔ−１に変更すると、以下のような式になる。
【００４２】
【数１４】
0005477238.tif 000015
【００４３】
そして、（１８）及び（１９）式を用いて、ａt−ａt-1を表すと、以下のような式で表される。
【００４４】
【数１５】
0005477238.tif 000016
【００４５】
また、（２０）式は、以下のように変形される。
【００４６】
【数１６】
0005477238.tif 000017
【００４７】
ここで、変数Ｋ＝ｋ／ｍ＊Δｔ2とし、誤差項をＹt＝ｙt−ｙt-1とすると、（２１）式は、以下のように変形される。
【００４８】
【数１７】
0005477238.tif 000018
【００４９】
そして、（２２）式に最小二乗法を適用すると、以下のような式が導出される。
【００５０】
【数１８】
0005477238.tif 000019
【００５１】
そして、変数Ｂjを以下のように定義する。
【００５２】
【数１９】
0005477238.tif 000020
【００５３】
そして、（２３）式を変数Ｂjを用いて表すと、以下のような式で表される。
【００５４】
【数２０】
0005477238.tif 000021
【００５５】
また、（２４）式は、以下のように変形される。
【００５６】
【数２１】
0005477238.tif 000022
【００５７】
そして、（２５）式をｖ0で偏微分したものを０とすると、以下のような式になる。
【００５８】
【数２２】
0005477238.tif 000023
【００５９】
そして、（２６）式を、ｖ0を求めるように変形すると、以下のような式になる。
【００６０】
【数２３】
0005477238.tif 000024
【００６１】
そして、（２７）式において、変数Ｂjを元の式に戻すと、以下のような式になる。
【００６２】
【数２４】
0005477238.tif 000025
【００６３】
本技術の実施の形態では、水平面におけるＮ個の加速度値ａtを（１７）式に代入することで変数Ｋの値を算出し、Ｎ個の加速度値ａt及び変数Ｋの値を（２８）式に代入することで、初速度値ｖ0を算出する。そして、Ｎ個の加速度値ａt及び初速度値ｖ0を（２）式に代入することで、重心位置ｘtを算出する。なお、重心動揺は相対的な位置であるため、初期位置ｘ0＝０を（２）式に代入する。これにより、計測開始時に静止していなくても、加速度センサによって計測された加速度値から身体の重心動揺を高精度に推定することができるようになる。
【００６４】
以下、このような処理を実現するための具体的な構成について説明する。
【００６５】
［実施の形態１］
図１に、本技術の第１の実施の形態に係る携帯端末１の機能ブロック図を示す。本実施の形態に係る携帯端末１は、加速度センサ１０１と、計測結果格納部１０３と、補正部１０５と、データ格納部１０７と、データ読出部１０９と、定数格納部１１１と、変数値算出部１１３と、変数値格納部１１５と、初速度算出部１１７と、初速度格納部１１９と、重心動揺計算部１２１と、重心動揺データ格納部１２３と、体調データベース（ＤＢ）１２５と、体調推定部１２７と、体調履歴データ格納部１２９と、出力部１３１とを有する。なお、携帯端末１は、被計測者の衣類等に取り付けられているものとする。
【００６６】
加速度センサ１０１は、所定のサンプリング間隔で加速度を計測し、計測結果である加速度データを計測結果格納部１０３に格納する。補正部１０５は、計測結果格納部１０３からデータを読み出し、所定の補正を行い、補正後の加速度データをデータ格納部１０７に格納する。なお、補正部１０５は、オフセット補正部１０５１と、重力方向補正部１０５２と、進行方向補正部１０５３と、歩行データ処理部１０５４と、パラメータ格納部１０５５と、区間データ格納部１０５６と、テンプレート格納部１０５７とを含む。また、データ読出部１０９は、所定期間毎に、データ格納部１０７から補正後の加速度データを読み出し、変数値算出部１１３と初速度算出部１１７と重心動揺計算部１２１とに出力する。変数値算出部１１３は、データ読出部１０９からのデータと定数格納部１１１に格納されているデータとを用いて、上記（１７）式に従って、変数Ｋの値を算出し、算出結果を変数値格納部１１５に格納する。初速度算出部１１７は、データ読出部１０９からのデータと定数格納部１１１及び変数値格納部１１５に格納されているデータとを用いて、上記（２８）式に従って初速度値ｖ0を算出し、算出結果を初速度格納部１１９に格納する。重心動揺計算部１２１は、データ読出部１０９からのデータと定数格納部１１１及び初速度格納部１１９に格納されているデータとを用いて被測定者の重心位置を計算し、計算結果を重心動揺データ格納部１２３に格納する。体調推定部１２７は、重心動揺データ格納部１２３と体調ＤＢ１２５と体調履歴データ格納部１２９とに格納されているデータを用いて、後で説明する体調推定処理を実施し、処理結果に基づき体調履歴データ格納部１２９を更新したり、体調の変化があった場合には体調変化通知を出力部１３１に出力したりする。出力部１３１は、体調推定部１２７から体調変化通知を受け取ると、体調履歴データ格納部１２９に格納されているデータに基づき、変化後の体調を出力する。
【００６７】
次に、図２乃至図３０を用いて、携帯端末１の処理について説明する。なお、本実施の形態では、重心動揺を推定するための前処理として、加速度センサ１０１によって計測された加速度データに対して所定の補正を行う。最初に、図２乃至図１５を用いて、この前処理について説明する。
【００６８】
まず、加速度センサ１０１は、所定のサンプリング間隔で加速度を計測し、計測結果格納部１０３に順次格納する（図２：ステップＳ１）。なお、サンプリング間隔は、サンプリング周波数（例えば２０Ｈｚ）の逆数で求められる。また、計測結果格納部１０３には、例えば図３に示すようなデータが格納される。図３の例では、計測結果格納部１０３には、時刻と、ｘ軸加速度と、ｙ軸加速度と、ｚ軸加速度とが格納されるようになっている。
【００６９】
そして、補正部１０５は、処理タイミングになったか判断する（ステップＳ３）。ここでは、所定期間（例えば１５秒）ごとに処理を実施するようになっているものとし、前回の処理から所定期間経過したか判断することによって、処理タイミングになったかどうか判断する。処理タイミングではないと判断された場合（ステップＳ３：Ｎｏルート）、ステップＳ１の処理に戻る。
【００７０】
一方、処理タイミングになったと判断された場合（ステップＳ３：Ｙｅｓルート）、補正部１０５は、処理対象となる加速度データを計測結果格納部１０３から読み出す（ステップＳ５）。すなわち、前回の処理から今回の処理までの間に計測された所定期間分の加速度データを読み出す。そして、オフセット補正部１０５１が、パラメータ格納部１０５５に格納されているデータを用いて、読み出された加速度データに対してオフセット補正を実施し（ステップＳ７）、補正後の加速度データを記憶装置に格納する。例えば加速度センサのｘ軸が重力方向、ｙ軸が進行方向、ｚ軸が左右方向となっている場合に、静止状態であれば、ｘ軸方向の加速度ａx＝１０００ｍＧ、ｙ軸方向の加速度ａy＝０ｍＧ、ｚ軸方向の加速度ａz＝０ｍＧといった値が本来ならば計測されるはずである。しかしながら、加速度センサなどのセンサには個体差が存在するため、静止状態で、例えばａx＝１００５ｍＧ、ａy＝１０ｍＧ、ａz＝２０ｍＧのように、誤差を含む値が計測されることがある。そこで、本実施の形態では、加速度が０Ｇの時には、加速度センサの測定値も０Ｇとなるようにし、加速度が１Ｇの時には、加速度センサの測定値も１Ｇとなるようにするためにオフセット補正を行う。
【００７１】
なお、パラメータ格納部１０５５には、例えば図４に示すようなデータが予め格納されている。図４の例では、パラメータ格納部１０５５には、ＸＦｕｌｌＳｃａｌｅ、ＸＯｆｆｓｅｔ、ＹＦｕｌｌＳｃａｌｅ、ＹＯｆｆｓｅｔ、ＺＦｕｌｌＳｃａｌｅ及びＺＯｆｆｓｅｔのパラメータ値が格納されている。例えばｘ軸方向の加速度については、ＸＦｕｌｌＳｃａｌｅ及びＸＯｆｆｓｅｔを用いて、ａｘ／ＸＦｕｌｌＳｃａｌｅ−ＸＯｆｆｓｅｔによって補正後の加速度を算出する。例えば図５に、オフセット補正前後の加速度データの一例を示す。なお、図５は、ＸＦｕｌｌＳｃａｌｅ＝１．０２６、ＸＯｆｆｓｅｔ＝０．０１５である場合の例を示す。例えば、時刻が「10:30:46.002」のオフセット補正前のｘ軸加速度についてオフセット補正を行うと、−９０６／１．０２６−０．０１５≒−８８３．０５６がオフセット補正後の加速度となる。
【００７２】
なお、これらのパラメータ値は、事前に、以下のような計測を行うことによって算出される。ここでは、ＸＦｕｌｌＳｃａｌｅ及びＸＯｆｆｓｅｔの算出の仕方を一例として説明する。具体的には、
［１］ｘ軸を重力方向に合わせて、ｘ軸の加速度を計測する
［２］ｘ軸の上下を逆にして、ｘ軸の加速度を計測する
［３］ｚ軸を重力方向に合わせて、ｘ軸の加速度を計測する
［４］ｚ軸の上下を逆にして、ｘ軸の加速度を計測する
そして、［１］で計測された値をＸＦｕｌｌＳｃａｌｅ１とし、［２］で計測された値をＸＦｕｌｌＳｃａｌｅ２とすると、ＸＦｕｌｌＳｃａｌｅ＝（ＸＦｕｌｌＳｃａｌｅ１−ＸＦｕｌｌＳｃａｌｅ２）／２によって算出される。また、［３］で計測された値をＸＯｆｆｓｅｔ１とし、［４］で計測された値をＸＯｆｆｓｅｔ２とすると、ＸＯｆｆｓｅｔ＝（ＸＯｆｆｓｅｔ１＋ＸＯｆｆｓｅｔ２）／２によって算出される。なお、［１］及び［２］において、ｙ軸を重力方向に合わせて、ｙ軸の加速度（すなわち、ＹＦｕｌｌＳｃａｌｅ１及び２）を計測することで、ＹＦｕｌｌＳｃａｌｅを算出できる。また、［３］及び［４］において、ｘ軸を重力方向に合わせて、ｙ軸の加速度（すなわち、ＹＯｆｆｓｅｔ１及び２）を計測することで、ＹＯｆｆｓｅｔを算出できる。さらに、［１］及び［２］において、ｚ軸を重力方向に合わせて、ｚ軸の加速度（すなわち、ＺＦｕｌｌＳｃａｌｅ１及び２）を計測することで、ＺＦｕｌｌＳｃａｌｅを算出できる。また、［３］及び［４］において、ｙ軸を重力方向に合わせて、ｚ軸の加速度（すなわち、ＺＯｆｆｓｅｔ１及び２）を計測することで、ＺＯｆｆｓｅｔを算出できる。
【００７３】
その後、重力方向補正部１０５２が、オフセット補正後の加速度データに対して重力方向補正を実施し（ステップＳ９）、補正後の加速度データを記憶装置に格納する。例えばフレームマトリックス法による座標変換を行うことで、加速度データを補正する。例えば図６に示すように単位ベクトルｉ、ｊ、ｋと、ｙｚ平面とのなす角度をそれぞれθｙ、θｚ、θｘとする。ここで、重力加速度ｇのｉ、ｊ、ｋ方向の成分をａy、ａz、ａxとすると、θy、θz、θxは以下のように定義される。
θy＝ｓｉｎ-1ａy／ｇ
θz＝ｓｉｎ-1ａz／ｇ
θx＝ｓｉｎ-1ａx／ｇ
【００７４】
また、これらの式からフレームマトリクスＥ0は以下のように定義される。
【数２５】
0005477238.tif 000026
【００７５】
そして、行列Ｅ0を用いて、Ｅ0＊ａによって、計測された加速度ａを、ｘ軸を重力方向とする座標系における加速度に変換することができる。
【００７６】
その後、進行方向補正部１０５３が、重力方向補正後の加速度データを用いて、重力方向の加速度データにおいて極大値の現れる時刻ＴＭを特定する（ステップＳ１１）。そして、進行方向補正部１０５３は、極大値の位置で区切られる各区間について、ステップＳ１１において特定されたＴＭを用いて歩行周期Ｔｗを算出し、区間データ格納部１０５６に格納する（ステップＳ１３）。具体的には、歩行周期をＴｗ＝ＴＭi−ＴＭi-1によって算出する。なお、区間データ格納部１０５６には、例えば図７に示すようなデータが格納される。図７の例では、区間データ格納部１０５６には、区間の列と、開始時刻の列と、終了時刻の列と、歩行周期Ｔｗの列と、歩行区間であるか否かを表す歩行フラグ（１：歩行区間、０：静止区間）の列と、後で説明する標準加速度差引処理が完了しているか否かを表す完了フラグ（１：完了、０：未完了）の列とが含まれる。例えば区間１について、ＴＭ1とＴＭ2とからＴｗ（＝ＴＭ2−ＴＭ1）を算出すると、ＴＭ1を開始時刻の欄に設定し、ＴＭ2を終了時刻の欄に設定し、（ＴＭ2−ＴＭ1）を歩行周期Ｔｗの欄に設定する。また、歩行フラグ及び完了フラグの欄には、０を設定しておく。
【００７７】
そして、進行方向補正部１０５３は、区間データ格納部１０５６から直前の４歩分の歩行周期Ｔｗを取得する（ステップＳ１５）。例えば図８に示すように、Ｔnの時点では、直前の４歩に該当する期間Ａにおける歩行周期Ｔｗを取得する。なお、図８の例では、左から右に時間が進んでいるものとし、黒塗りの丸印は、極大値の出現位置を表している。その後、端子Ａを介して図９の処理に移行する。
【００７８】
図９の説明に移行して、端子Ａの後、進行方向補正部１０５３は、直前の４歩分の歩行周期Ｔｗを用いて、歩行中であるかどうか判断する（図９：ステップＳ１７）。具体的には、８０〜１５０［ステップ／分］であり且つ４歩分の歩行周期Ｔｗのばらつきが±１０％以内である場合には、歩行中であると判断する。
【００７９】
歩行中であると判断された場合（ステップＳ１７：Ｙｅｓルート）、進行方向補正部１０５３は、区間データ格納部１０５６において今回処理対象となっている各区間についての歩行フラグを１に設定する。そして、進行方向補正部１０５３は、回転行列算出処理を実施する（ステップＳ１９）。この回転行列算出処理については図１０を用いて説明する。
【００８０】
まず、進行方向補正部１０５３は、区間データ格納部１０５６に格納されているデータに従って、重力方向補正後の加速度データが格納されている記憶装置から直前の４歩分の加速度データを読み出す（図１０：ステップＳ３１）。なお、直前の４歩に該当する期間の開始時刻及び終了時刻は、区間データ格納部１０５６に格納されているデータから取得できる。
【００８１】
そして、進行方向補正部１０５３は、４歩分の加速度データのうち、ｙｚ平面における加速度の大きさが最大となる加速度データ（ａy，ａz）を特定する（ステップＳ３３）。例えば４歩に該当する期間におけるデータ数がｎであった場合、（ａy1，ａz1）、（ａy2，ａz2）、（ａy3，ａz3）、・・・、（ａyn，ａyn）の中から、ベクトルの大きさが最大となる（ａy，ａz）を特定する。なお、歩行時には、進行方向に最も大きな加速度が生じる。
【００８２】
そして、進行方向補正部１０５３は、特定された加速度データ（ａy，ａz）を用いて、軸のずれを補正するための角度φを算出する（ステップＳ３５）。なお、角度φは、図１１に示すように、進行方向（ｙ軸）と、ベクトル（ａy，ａz）のなす角である。具体的には、以下の式によって角度φを算出する。
φ＝ｔａｎ-1（ａz／ａy）
【００８３】
そして、進行方向補正部１０５３は、算出された角度φを用いて回転行列Ｆ0を算出し、パラメータ格納部１０５５に格納する（ステップＳ３７）。具体的には、以下の式によって回転行列Ｆ0を算出する。
【００８４】
【数２６】
0005477238.tif 000027
その後、処理を終了し、元の処理に戻る。
【００８５】
図９の説明に戻って、ステップＳ１９の後、進行方向補正部１０５３は、ステップＳ１９において算出された回転行列Ｆ0を用いて処理対象の加速度データを補正し（ステップＳ２１）、補正後の加速度データを記憶装置に格納する。例えば補正前の加速度データをａy[i]及びａz[i]とし、補正後の加速度データをａ’y[i]及びａ’zi[i]とすると、以下の式によって補正後の加速度データを求める。
【数２７】
0005477238.tif 000028
【００８６】
なお、ここまでの処理で、軸の向きを合わせるための変換が完了したことになる。
【００８７】
その後、歩行データ処理部１０５４が、区間データ格納部１０５６及びテンプレート格納部１０５７に格納されているデータを用いて標準加速度差引処理を実施する（ステップＳ２３）。この標準加速度差引処理については図１２を用いて説明する。
【００８８】
まず、歩行データ処理部１０５４は、テンプレート格納部１０５７から、歩行時における進行方向の標準加速度を表すテンプレートｐyと、歩行時における左右方向の標準加速度を表すテンプレートｐzを取得する（ステップＳ４１）。
【００８９】
なお、図１３に、テンプレート格納部１０５７に予め格納されるテンプレートの一例を示す。なお、図１３は、歩行時における進行方向の標準加速度を表すテンプレートの一例を示している。図１３では、横軸は歩行周期を表しており、縦軸は加速度を表している。テンプレートは、例えば歩行加速度の特徴的パターンから生成されるデータである。なお、歩行周期については、６０［ステップ／分］となるように正規化されており、加速度については、極大値が１Ｇ且つ極小値が０Ｇとなるように正規化されているものとする。
【００９０】
なお、非特許文献１（小椋，大渕，小島，古名，潮見，「通常歩行時の骨盤加速度に注目した歩行分析−若年者における歩行加速度の特徴的パターンの同定に関する検討−」，理学療法科学，第２０巻第２号，ｐ１７１−ｐ１７７，２００５）で示されているように、左右方向の標準加速度は、左足の１歩と右足の１歩とで加速度の特徴的パターンが異なってくる。従って、本実施の形態では、左足の１歩についての標準加速度を表すテンプレートと、右足の１歩についての標準加速度を表すテンプレートとの２種類を用意しておく。
【００９１】
そして、歩行データ処理部１０５４は、区間データ格納部１０５６に格納されているデータに従って、歩行区間のうち未処理の区間を１つ特定する。具体的には、区間データ格納部１０５６において、歩行フラグが１であり且つ完了フラグが０である区間の中から１つ特定する。そして、歩行データ処理部１０５４は、特定した区間に対応する歩行周期Ｔｗを区間データ格納部１０５６から取得し、特定した区間内の加速度データａy及びａzを記憶装置から取得する（ステップＳ４３）。ここでは、進行方向の加速度データをａyと記し、左右方向の加速度データをａzと記す。そして、歩行データ処理部１０５４は、取得した加速度データａyから極大値Ｍy及び極小値ｍyを特定する。同様に、歩行データ処理部１０５４は、取得した加速度データａzから極大値Ｍz及び極小値ｍzを特定する（ステップＳ４５）。
【００９２】
そして、歩行データ処理部１０５４は、テンプレートｐyをＴｗ倍することで時間軸方向に伸縮する（ステップＳ４７）。例えば、図１３に示したテンプレートをＴｗ倍して時間軸方向に伸縮させると、図１４に示すようなテンプレートに変形される。図１４は、図１３に示したテンプレートを時間軸方向（横軸）に縮めた場合の一例を示している。なお、テンプレートを時間軸方向に伸ばす際には、補間処理などを行う。
【００９３】
そして、歩行データ処理部１０５４は、時間軸方向に伸縮されたテンプレートｐyを、極大値Ｍy及び極小値ｍyを用いて加速度方向に伸縮する（ステップＳ４９）。上で述べたように、テンプレートは、極大値が１Ｇ且つ極小値が０Ｇとなるように正規化されている。従って、本ステップでは、テンプレートを、計測された加速度の振幅に合うように変形する。具体的には、以下の式を用いて変形する。
（Ｍy−ｍy）＊ｐy＋ｍy
【００９４】
例えば、図１４に示したテンプレートを加速度方向に伸縮させると、図１５に示すようなテンプレートに変形される。図１５は、図１４に示したテンプレートを加速度方向（縦軸）に伸ばした場合の一例を示している。
【００９５】
そして、歩行データ処理部１０５４は、伸縮後のテンプレートｐyを用いて、特定された区間内のａyから標準加速度を取り除き、処理結果をデータ格納部１０７に格納する（ステップＳ５１）。具体的には、伸縮後のテンプレートｐyから得られる標準加速度をａpyとすると、ａy−ａpyによって、標準加速度を取り除いた加速度を算出する。
【００９６】
そして、歩行データ処理部１０５４は、特定された区間内のａzから、特定された区間が、左足の１歩の区間であるか又は右足の１歩の区間であるか判別する。なお、上で説明した処理では、進行方向が正の方向になるように回転行列Ｆ0を適用するため、例えば右手座標系の場合、回転行列Ｆ0を適用することで、左方向が正の方向となる。この場合、特定された区間内のａzの平均値を算出し、平均値が正の値であれば、右足の一歩と判断でき、平均値が負の値であれば、左足の一歩と判断できる。なお、以下の処理では、左足の１歩についての標加速度を表すテンプレートと、右足の１歩についての標準加速度を表すテンプレートとのうち、該当するテンプレートが用いられる。
【００９７】
そして、歩行データ処理部１０５４は、テンプレートｐzをＴｗ倍することで時間軸方向に伸縮する（ステップＳ５３）。そして、歩行データ処理部１０５４は、時間軸方向に伸縮されたテンプレートｐzを、極大値Ｍz及び極小値ｍzを用いて加速度方向に伸縮する（ステップＳ５５）。具体的には、以下の式を用いて変形する。
（Ｍz−ｍz）＊ｐｚ＋ｍz
【００９８】
そして、歩行データ処理部１０５４は、伸縮後のテンプレートｐzを用いて、特定された区間内のａzから標準加速度を取り除き、処理結果をデータ格納部１０７に格納する（ステップＳ５７）。具体的には、伸縮後のテンプレートｐzから得られる標準加速度をａpzとすると、ａz−ａpzによって、標準加速度を取り除いた加速度を算出する。
【００９９】
そして、歩行データ処理部１０５４は、区間データ格納部１０５６において、特定された区間についての完了フラグを１（完了）に設定する。
【０１００】
そして、歩行データ処理部１０５４は、歩行区間の中に未処理の区間があるか判断する（ステップＳ５９）。歩行区間の中に未処理の区間がある場合（ステップＳ５９：Ｙｅｓルート）、ステップＳ４３に戻り、上で述べた処理を繰り返す。一方、歩行区間の中に未処理の区間がなければ（ステップＳ５９：Ｎｏルート）、処理を終了し、元の処理（図９）に戻る。以上のような処理を実施することで、計測された加速度から、歩行時の標準加速度を取り除くことができる。
【０１０１】
図９の説明に戻って、ステップＳ２３の後、処理を終了する。
【０１０２】
一方、図９のステップＳ１７において、歩行中ではないと判断された場合（ステップＳ１７：Ｎｏルート）、ステップＳ２５の処理に移行する。そして、進行方向補正部１０５３は、過去に算出された回転行列Ｆ0があるかどうか判断する（ステップＳ２５）。上で説明したように、回転行列算出処理によって算出された回転行列Ｆ0は、パラメータ格納部１０５５に格納されるので、パラメータ格納部１０５５に回転行列Ｆ0が格納されているか否かで判断する。過去に算出された回転行列Ｆ0がある場合には（ステップＳ２５：Ｙｅｓルート）、進行方向補正部１０５３は、過去に算出された回転行列Ｆ0を用いて処理対象の加速度データを補正し（ステップＳ２１）、補正後の加速度データをデータ格納部１０７に格納する。なお、補正処理自体は、ステップＳ２１の処理と変わらないので、ここではこれ以上述べない。その後、処理を終了する。
【０１０３】
一方、過去に算出された回転行列Ｆ0がない場合には（ステップＳ２５：Ｎｏルート）、ステップＳ２７の処理を行わず、処理を終了する。
【０１０４】
以上のような処理を実施することにより、加速度センサによって取得された加速度データを補正することができる。なお、補正後の加速度データは、データ格納部１０７に格納される。また、取得した加速度データに対して、ローパスフィルタをかけるなどして、ノイズを除去するようにしてもよい。
【０１０５】
次に、図１６乃至図３０を用いて、重心動揺及び体調を推定する処理について説明する。本実施の形態では、図１６に示すような処理フローを、所定期間（例えば１５秒）ごとに実施する。まず、処理タイミングになると、データ読出部１０９は、処理対象となる加速度データをデータ格納部１０７から読み出し（図１６：ステップＳ７１）、変数値算出部１１３と初速度算出部１１７と重心動揺計算部１２１とに出力する。図１７に模式的に示すように、例えば所定期間におけるサンプリング回数をＮとすると、データ読出部１０９は、データ格納部１０７からＮ個の補正後の加速度値ａiを読み出す。
【０１０６】
そして、変数値算出部１１３は、データ読出部１０９からのＮ個の加速度値ａiと定数格納部１１１に格納されている所定の定数とを上記（１７）式に代入して、変数Ｋの値を算出し、算出結果を変数値格納部１１５に格納する（ステップＳ７３）。本ステップでは、進行方向のＮ個の加速度値ａiを用いて進行方向についての変数Ｋの値を算出し、左右方向のＮ個の加速度値ａiを用いて左右方向についての変数Ｋの値を算出する。ここでは、ｙ軸を進行方向とし、ｚ軸を左右方向として説明する。なお、定数格納部１１１には、サンプリング間隔Δｔ、サンプリング回数Ｎ、所定期間などを表すデータが予め格納されている。なお、変数値格納部１１５には、例えば図１８に示すようなデータが格納される。図１８の例では、変数値格納部１１５には、期間の列と、Ｋの値（進行方向）の列と、Ｋの値（左右方向）の列とが含まれる。
【０１０７】
その後、初速度算出部１１７は、データ読出部１０９からのＮ個の加速度値と、変数値格納部１１５に格納されている変数Ｋの値と、定数格納部１１１に格納されている所定の定数とを上記（２８）式に代入して、初速度値ｖ0を算出し、算出結果を初速度格納部１１９に格納する（ステップＳ７５）。本ステップでは、進行方向のＮ個の加速度値ａi及び進行方向についての変数Ｋの値を用いて、進行方向についての初速度値ｖ0を算出する。同様に、左右方向のＮ個の加速度値ａi及び左右方向についての変数Ｋの値を用いて、左右方向についての初速度値ｖ0を算出する。なお、初速度格納部１１９には、例えば図１９に示すようなデータが格納される。図１９の例では、初速度格納部１１９には、期間の列と、初速度（進行方向）の列と、初速度（左右方向）の列とが含まれる。
【０１０８】
そして、重心動揺計算部１２１は、データ読出部１０９からのＮ個の加速度値と、初速度格納部１１９に格納されている初速度値ｖ0と、定数格納部１１１に格納されている所定の定数とを上記（２）式に代入して、重心位置を計算し、計算結果を重心動揺データ格納部１２３に格納する（ステップＳ７７）。本ステップでは、進行方向のＮ個の加速度値ａi及び進行方向についての初速度値ｖ0を用いて、ｙ軸の座標値を算出する。同様に、左右方向のＮ個の加速度値ａi及び左右方向についての初速度値ｖ0を用いて、ｚ軸の座標値を算出する。なお、重心動揺データ格納部１２３には、例えば図２０に示すようなデータが格納される。図２０の例では、重心動揺データ格納部１２３には、時刻の列と、位置（ｙ軸）の列と、位置（ｚ軸）の列とが含まれる。
【０１０９】
その後、体調推定部１２７が、重心動揺データ格納部１２３と体調ＤＢ１２５と体調履歴データ格納部１２９とに格納されているデータを用いて体調推定処理を実施し（ステップＳ７９）、その後、処理を終了する。なお、体調推定処理では、複数の体調候補（例えば、疲労、睡眠不足、風邪、熱中症、酔い、緊張等）についての代表データを用意しておき、クラスタ分析を行うことによって体調を推定する。例えば図２１にクラスタ分析に用いる代表データを模式的に示す。図２１では、風邪状態についての代表データと、睡眠不足についての代表データと、熱中症についての代表データとが示されている。以下、図２２を用いて体調推定処理について説明する。
【０１１０】
まず、体調推定部１２７は、重心動揺データ格納部１２３から位置データを読み出し、重心動揺軌跡の中心変位Ｄを算出する（図２２：ステップＳ８１）。本ステップでは、ｙ軸及びｚ軸の座標値の平均値を求め、中心変位Ｄとする。そして、体調推定部１２７は、重心動揺データ格納部１２３に格納されている座標値を、中心変位Ｄを原点とする場合の座標値に変換し、記憶装置に格納する（ステップＳ８３）。例えば図２０に示したデータを、中心変位Ｄを原点とする相対位置に変換すると、図２３に示すようなデータとなる。そして、体調推定部１２７は、中心変位Ｄを通る３度間隔の直線により区切られる領域毎に、重心動揺データ格納部１２３に格納されている座標値で表される点のうち、中心変位Ｄからの距離が最大となる点を特定する（ステップＳ８５）。本ステップの処理を図２４に模式的に示す。図２４に示すように、中心変位Ｄを通る直線によってｙｚ平面を３度毎の領域に分割する。そして、領域毎に、中心変位Ｄからの距離が最大となる点（以下、第１の点と呼ぶ）を特定する。例えば図２５に示すように、領域毎に、領域の範囲と領域番号と第１の点の座標値と中心変位Ｄから当該第１の点までの距離とを記憶装置に格納しておく。
【０１１１】
そして、体調推定部１２７は、体調ＤＢ１２５から未処理の体調候補を１つ特定し、その体調候補に対応する、各領域内の代表点（以下、第２の点と呼ぶ）を抽出する（ステップＳ８７）。体調ＤＢ１２５には、例えば図２６に示すようなデータが予め格納されている。図２６の例では、体調ＤＢ１２５には、領域の範囲の列と、領域番号の列と、正常状態の列と、疲労状態の列と、・・・とが含まれる。体調ＤＢ１２５には、体調候補毎に、当該体調候補に対応する、各領域内の代表点のデータが格納されている。なお、代表点の位置は、座標値やベクトル値などを用いて表される。
【０１１２】
そして、体調推定部１２７は、第１の点と第２の点との間のユークリッド距離を算出し、記憶装置に格納する（ステップＳ８９）。本実施の形態では、第１の点のデータ（座標値又は中心変位Ｄから第１の点までの距離）をｒi（ｉは領域番号を表す変数である）とし、ｋ番目の体調候補についての代表点のデータ値（座標値又はベクトル値）をＲｋiとして、以下の式によってユークリッド距離を算出する。
【数２８】
0005477238.tif 000029
【０１１３】
そして、体調推定部１２７は、体調ＤＢ１２５において未処理の体調候補があるか判断する（ステップＳ９１）。未処理の体調候補がある場合（ステップＳ９１：Ｙｅｓルート）、ステップＳ８７に戻り、上で述べた処理を繰り返す。
【０１１４】
一方、未処理の体調候補がなければ（ステップＳ９１：Ｎｏルート）、体調推定部１２７は、記憶装置に格納されているデータに基づき、算出された距離が最小となる体調候補を特定する（ステップＳ９３）。その後、端子Ｂを介して図２７の処理に移行する。
【０１１５】
図２７の説明に移行して、端子Ｂの後、体調推定部１２７は、算出された距離のうち最小距離が所定の閾値以下であるか判断する（図２７：ステップＳ９５）。最小距離が所定の閾値以下ではない場合（ステップＳ９５：Ｎｏルート）、すなわち最小距離が所定の閾値より大きい場合には、体調推定部１２７は、体調不明と判断する。そして、体調推定部１２７は、体調不明を表すデータを体調履歴データ格納部１２９に登録する（ステップＳ９７）。体調履歴データ格納部１２９には、図２８及び図２９に示すようなデータが格納される。図２８は、体調履歴データ格納部１２９に格納される体調履歴データの一例を示す。また、図２９は、体調履歴データ格納部１２９に格納される体調管理データの一例を示す。ステップ９７では、体調不明を表すデータを含むレコードを体調履歴データ（図２８）の一番下に追加する。
【０１１６】
そして、体調推定部１２７は、体調履歴データに従って、体調不明が４回連続したか判断する（ステップＳ９９）。体調不明が４回連続した場合（ステップＳ９９：Ｙｅｓルート）、体調推定部１２７は、体調管理データに含まれるリセットフラグをＯＮに設定する（ステップＳ１０１）。その後、処理を終了し、元の処理に戻る。一方、体調不明が４回連続していなければ（ステップＳ９９：Ｎｏルート）、処理を終了し、元の処理に戻る。
【０１１７】
一方、ステップＳ９５において最小距離が閾値以下であると判断された場合（ステップＳ９５：Ｙｅｓルート）、体調推定部１２７は、特定された体調候補を表すデータを体調履歴データ格納部１２９に登録する（ステップＳ１０３）。具体的には、特定された体調候補を表すデータを含むレコードを体調履歴データの一番下に追加登録する。そして、体調推定部１２７は、体調履歴データに従って、同じ体調候補が４回連続したか判断する（ステップＳ１０５）。同じ体調候補が４回連続していなければ（ステップＳ１０５：Ｎｏルート）、以下で説明する処理を行わず、処理を終了し、元の処理に戻る。
【０１１８】
一方、同じ体調候補が４回連続した場合（ステップＳ１０５：Ｙｅｓルート）、体調推定部１２７は、体調管理データに含まれるリセットフラグがＯＮに設定されているか判断する（ステップＳ１０７）。リセットフラグがＯＦＦに設定されている場合（ステップＳ１０７：Ｎｏルート）、体調推定部１２７は、体調が変化したか判断する（ステップＳ１０９）。具体的には、４回連続した体調候補と、体調管理データに含まれる現在の体調とを比較する。そして、両者が一致する場合には変化なしと判断し、一致しない場合には変化ありと判断する。体調が変化していなければ（ステップＳ１０９：Ｎｏルート）、以下で説明する処理を行わず、処理を終了し、元の処理に戻る。
【０１１９】
一方、体調が変化したと判断された場合（ステップＳ１０９：Ｙｅｓルート）、体調推定部１２７は、体調変化通知を出力するよう出力部１３１に指示する。そして、出力部１３１は、体調推定部１２７からの指示に応じて、体調履歴データ格納部１２９に格納されている体調履歴データに従って、変化後の体調を表すデータを含む体調変化通知を出力する（ステップＳ１１１）。例えば、携帯端末１が表示装置を備えている場合には、変化後の体調を表示したり、携帯端末１がスピーカーなどを備えている場合には、音声データを出力したりする。その後、ステップＳ１１５に移行する。
【０１２０】
一方、リセットフラグがＯＮに設定されている場合（ステップＳ１０７：Ｙｅｓルート）、体調推定部１２７は、リセットフラグをＯＦＦに設定する（ステップＳ１０９）。その後、ステップＳ１１５に移行する。
【０１２１】
そして、ステップＳ１１１又はＳ１１３の後、体調推定部１２７は、特定された体調候補を現在の体調として特定し、体調履歴データ格納部１２９に格納する（ステップＳ１１５）。具体的には、特定された体調候補を表すデータが、体調履歴データ格納部１２９内の体調管理データに設定される。その後、処理を終了し、元の処理に戻る。
【０１２２】
以下、体調推定処理をより分かりやすくするために、図３０を用いて具体例を説明する。例えば図３０（ａ）に体調変化の第１の例を示す。図３０（ａ）では、体調Ａの後、体調Ｂが連続して検出されている。図３０（ａ）に示すように、体調Ｂが４回連続した時点で、体調Ｂに変化したと判断し、体調変化を通知する。
【０１２３】
また、図３０（ｂ）に体調変化の第２の例を示す。図３０（ｂ）では、体調Ａの後、体調Ｂが２回検出され、その後に体調Ｃが１回検出されている。そして、体調Ｃの後、体調Ｂが連続して検出されている。図３０（ｂ）に示すように、体調Ｃが検出された後、体調Ｂが４回連続した時点で、体調Ｂに変化したと判断し、体調変化を通知する。
【０１２４】
さらに、図３０（ｃ）に体調変化の第３の例を示す。図３０（ｃ）では、体調Ａの後、体調不明が２回検出され、その後に体調Ｂが連続して検出されている。図３０（ｃ）に示すように、体調不明が検出された後、体調Ｂが４回連続した時点で、体調Ｂに変化したと判断し、体調変化を通知する。
【０１２５】
また、図３０（ｄ）に体調変化の第４の例を示す。図３０（ｄ）では、体調Ａの後、体調不明が４回連続して検出され、その後に体調Ｂが連続して検出されている。この場合、体調不明が４回連続した時点でリセットフラグがＯＮに設定されるため、体調Ｂが４回連続しても、体調変化を通知しない。
【０１２６】
以上のような処理を実施することによって、計測開始時に静止していなくても、加速度センサ１０１によって計測された加速度値から重心位置を精度良く算出できるようになる。さらに、算出された重心位置から重心動揺を高精度に推定でき、体調を推定できるようになる。
【０１２７】
［第２の実施の形態］
第１の実施の形態では、スタンドアロン型の実施の形態を説明したが、クライアント・サーバ型に変形することも可能である。図３１に、クライアント・サーバ型システムで実現する場合のシステム構成の一例を示す。図３１に示すように、例えばインターネット等であるネットワーク５０には、サーバ１０と、センサ付き端末２０と、車載装置３０とが無線又は有線にて接続されている。センサ付き端末２０は、例えばユーザの衣類等に取り付けられているものとする。
【０１２８】
センサ付き端末２０は、加速度センサ２０１と近接検知部２０３と送信部２０５とを有する。加速度センサ２０１は、所定のサンプリング間隔で加速度を計測し、送信部２０５に順次出力する。送信部２０５は、加速度センサ２０１により計測された加速度の値をネットワーク５０を介してサーバ１０に順次送信する。また、近接検知部２０３は、車載装置３０に近づいたことを検出すると、センサ付き端末２０及び車載装置３０の識別情報（例えばＩＰアドレスなど）を送信部２０５に出力する。送信部２０５は、近接検知部２０３から車載装置３０の識別情報を受け取ると、車載装置３０の識別情報を含む体調情報送信要求をネットワーク５０を介してサーバ１０に送信する。
【０１２９】
また、図３２に、サーバ１０の機能ブロック図を示す。サーバ１０は、受信部１４１と、計測結果格納部１０３と、補正部１０５と、データ格納部１０７と、データ読出部１０９と、定数格納部１１１と、変数値算出部１１３と、変数値格納部１１５と、初速度算出部１１７と、初速度格納部１１９と、重心動揺計算部１２１と、重心動揺データ格納部１２３と、体調ＤＢ１２５と、体調推定部１２７と、体調履歴データ格納部１２９と、送信部１４３とを有する。なお、計測結果格納部１０３と、補正部１０５と、データ格納部１０７と、データ読出部１０９と、定数格納部１１１と、変数値算出部１１３と、変数値格納部１１５と、初速度算出部１１７と、初速度格納部１１９と、重心動揺計算部１２１と、重心動揺データ格納部１２３と、体調ＤＢ１２５と、体調推定部１２７と、体調履歴データ格納部１２９とは、第１の実施の形態と同じものである。
【０１３０】
受信部１４１は、ネットワーク５０を介してセンサ付き端末２０から送信される加速度の値を順次受信し、計測結果格納部１０３に格納する。また、受信部１４１は、体調情報送信要求を受信すると、送信部１４３に出力する。送信部１４３は、受信部１４１から体調情報送信要求を受信すると、要求元のセンサ付き端末２０を装着しているユーザの体調を体調履歴データ格納部１２９から読み出す。なお、体調履歴データ格納部１２９には、上で説明した処理により、体調履歴データや体調管理データなどが格納される。そして、送信部１４３は、体調情報送信要求に含まれるＩＰアドレスなどを用いて、ユーザの体調を表すデータを車載装置３０に送信する。
【０１３１】
また、車載装置３０は、受信部３０１と近接検知部３０３と出力部３０５とを有する。受信部３０１は、サーバ１０からユーザの体調を表すデータを受信すると、出力部３０５に出力する。出力部３０５は、受信部３０１からデータを受信すると、表示装置に表示したり、音声にて通知したりする。
【０１３２】
例えばユーザが車に乗ろうとした際に、ユーザの体調が「体調不良」であった場合には、車載装置３０は、ユーザに運転を控える旨を通知するようにしてもよい。また、飲酒状態であった場合には、ユーザに警告したりやエンジンを始動させないようにしてもよい。
【０１３３】
なお、ここでは、車載装置３０を有するシステム構成を示したが、車載装置３０以外の装置を有するシステム構成であってもよい。
【０１３４】
以上本技術の実施の形態を説明したが、本技術はこれに限定されるものではない。例えば、上で説明した携帯端末１、サーバ１０、センサ付き端末２０及び車載装置３０の機能ブロック図は必ずしも実際のプログラムモジュール構成に対応するものではない。データ格納部の構成も同様に一例にすぎない。また、処理フローにおいても、処理結果が変わらなければ処理の順番を入れ替えることも可能である。さらに、並列に実行させるようにしてもよい。
【０１３５】
なお、上で述べた携帯端末１、サーバ１０、センサ付き端末２０及び車載装置３０は、コンピュータ装置であって、図３５に示すように、メモリ２５０１とＣＰＵ２５０３とハードディスク・ドライブ（ＨＤＤ）２５０５と表示装置２５０９に接続される表示制御部２５０７とリムーバブル・ディスク２５１１用のドライブ装置２５１３と入力装置２５１５とネットワークに接続するための通信制御部２５１７とがバス２５１９で接続されている。オペレーティング・システム（ＯＳ：Operating System）及び本実施例における処理を実施するためのアプリケーション・プログラムは、ＨＤＤ２５０５に格納されており、ＣＰＵ２５０３により実行される際にはＨＤＤ２５０５からメモリ２５０１に読み出される。ＣＰＵ２５０３は、アプリケーション・プログラムの処理内容に応じて表示制御部２５０７、通信制御部２５１７、ドライブ装置２５１３を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ２５０１に格納されるが、ＨＤＤ２５０５に格納されるようにしてもよい。本技術の実施例では、上で述べた処理を実施するためのアプリケーション・プログラムはコンピュータ読み取り可能なリムーバブル・ディスク２５１１に格納されて頒布され、ドライブ装置２５１３からＨＤＤ２５０５にインストールされる。インターネットなどのネットワーク及び通信制御部２５１７を経由して、ＨＤＤ２５０５にインストールされる場合もある。このようなコンピュータ装置は、上で述べたＣＰＵ２５０３、メモリ２５０１などのハードウエアとＯＳ及びアプリケーション・プログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。
【０１３６】
なお、携帯端末１及びセンサ付き端末２０が、携帯電話機の場合には、ＨＤＤ２５０５やドライブ装置２５１３は有していない。その場合、代わりにフラッシュメモリ型のメモリや記録媒体インタフェースなどを有する。また、携帯端末１及びセンサ付き端末２０については、表示制御部２５０７や表示装置２５０９を有しない場合もある。
【０１３７】
以上本実施の形態をまとめると以下のようになる。
【０１３８】
本実施の形態の第１の態様に係る情報処理方法は、（Ａ）物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップ（図３３：Ｓ１００１）と、（Ｂ）物体の位置を加速度及び初速度から求めるための第１の関係式並びに物体の重心動揺で物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、データ格納部に格納されている、所定期間内の加速度の補正後の値を代入して、物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップ（図３３：Ｓ１００３）と、（Ｃ）第１の関係式に、データ格納部に格納されている加速度の補正後の値と記憶装置に格納されている初速度の値とを代入して、物体の重心位置を表す、水平面における座標値を算出し、記憶装置に格納するステップ（図３３：Ｓ１００５）とを含む。
【０１３９】
このように、新たに第３の関係式を導入することで、計測開始時に物体が静止していなくても、加速度センサによって計測された加速度の値から物体の重心位置を表す、水平面の座標値を精度良く算出できるようになる。さらに、時系列の座標値から物体の重心動揺を推定することができる。
【０１４０】
また、上で述べた物体は人体である場合もある。そして、第１の態様において、記憶装置に格納されている座標値を用いて、人体の状態を特定する状態特定ステップをさらに含むようにしてもよい。これにより、人体がどのような状態であるかを把握できるようになる。
【０１４１】
さらに、上で述べた状態特定ステップが、記憶装置に格納されている座標値により表される、水平面上の点列についての中心点を通る直線によって水平面を所定角度毎の領域に分割するステップと、領域毎に、中心点からの距離が最大となり且つ座標値により表される点である第１の点を特定し、当該第１の点の座標値を中心点を原点とする場合の座標値に変換し、記憶装置に格納するステップと、人体の状態の候補毎に当該人体の状態の候補に対応する、各領域内の点である第２の点の座標値を格納する状態データベースを用いて、記憶装置に格納されている第１の点の座標値と第２の点の座標値とから算出され且つ全ての領域に関する距離が最小となる人体の状態の候補を特定するステップとを含むようにしてもよい。このようにすれば、予め用意される人体の状態の候補の中から、現在の状態に近いものを精度良く特定することができるようになる。
【０１４２】
また、上で述べた物体は人である場合もある。そして、上で述べた補正ステップが、人の歩行時の標準加速度を表すテンプレートを格納するテンプレート格納部に格納されているテンプレートを、歩行していたとみなされる期間において計測された加速度の値から特定される歩行周期及び加速度の振幅に合うように変形するステップと、変形後のテンプレートを用いて、歩行していたとみなされる期間において計測された加速度の値から歩行時の標準加速度を取り除くステップとを含むようにしてもよい。このように、加速度センサによって計測された加速度の値から歩行時の標準加速度を取り除くことで、歩行時の重心動揺を推定することができるようになる。
【０１４３】
さらに、上で述べた初速度算出ステップが、データ格納部に格納されている、所定期間内の加速度の補正後の値ａと、所定期間と所定のサンプリング間隔Δｔとから算出されるサンプル数Ｎとを用いて、上に示した（１７）式に従って、変数Ｋの値を算出するステップと、データ格納部に格納されている、所定期間内の加速度の補正後の値ａと、変数Ｋの値と、所定のサンプリング間隔Δｔと、サンプル数Ｎとを用いて、上に示した（２８）式に従って、物体の初速度ｖ0の値を算出するステップとを含むようにしてもよい。これにより、物体の加速度を高精度に算出することができる。
【０１４４】
また、第１の態様において、加速度センサにより所定のサンプリング間隔Δｔで加速度の値を計測し、計測結果格納部に順次格納するステップをさらに含むようにしてもよい。そして、所定期間におけるサンプリング回数をＮとし、上で述べた初速度算出ステップが、所定期間ごとに、データ格納部から、処理対象の所定期間に生成されたＮ個の加速度の補正後の値ａi（ｉは、０からＮ−１までの整数）を順次読み出すステップと、読み出されたＮ個の加速度の補正後の値ａiを、上に示した（１７）式に代入することで、処理対象の所定期間における変数Ｋの値を算出し、記憶装置に格納するステップと、読み出されたＮ個の加速度の補正後の値ａi及び記憶装置に格納された変数Ｋの値を、上に示した（２８）式に代入することで、処理対象の所定期間における物体の初速度ｖ0の値を算出し、記憶装置に格納するステップとを含むようにしてもよい。
【０１４５】
本実施の形態の第２の態様に係る情報処理装置（図３４）は、（Ａ）所定のサンプリング間隔で計測された物体の加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正部（図３４：１５０１）と、（Ｂ）物体の位置を加速度及び初速度から求めるための第１の関係式並びに物体の重心動揺で物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、データ格納部に格納されている、所定期間内の加速度の補正後の値を代入して、物体の初速度の値を算出し、記憶装置に格納する初速度算出部（図３４：１５０３）と、（Ｃ）第１の関係式に、データ格納部に格納されている加速度の補正後の値と記憶装置に格納されている初速度の値とを代入して、物体の重心位置を表す、水平面における座標値を算出し、記憶装置に格納する重心位置算出部（図３４：１５０５）とを有する。
【０１４６】
また、第２の態様において、所定のサンプリング間隔で物体の加速度の値を計測し、計測結果格納部に順次格納する加速度センサをさらに有する場合もある。そして、上で述べた補正部は、計測結果格納部から所定期間において計測された物体の加速度の値を読み出して、所定の補正を行うようにしてもよい。
【０１４７】
さらに、上で述べた物体は人体である場合もある。そして、第２の態様において、記憶装置に格納されている座標値により表される、水平面上の点列についての中心点を通る直線によって水平面を所定角度毎の領域に分割し、領域毎に、中心点からの距離が最大となり且つ座標値により表される点である第１の点を特定し、当該第１の点の座標値を中心点を原点とする場合の座標値に変換し、記憶装置に格納し、人体の状態の候補毎に当該人体の状態の候補に対応する、各領域内の点である第２の点の座標値を格納する状態データベースを用いて、記憶装置に格納されている第１の点の座標値と第２の点の座標値とから算出され且つ全ての領域に関する距離が最小となる人体の状態の候補を特定する状態特定部をさらに有するようにしてもよい。
【０１４８】
なお、上で述べたような処理をコンピュータに実施させるためのプログラムを作成することができ、当該プログラムは、例えばフレキシブル・ディスク、ＣＤ−ＲＯＭ、光磁気ディスク、半導体メモリ（例えばＲＯＭ）、ハードディスク等のコンピュータ読み取り可能な記憶媒体又は記憶装置に格納される。なお、処理途中のデータについては、ＲＡＭ等の記憶装置に一時保管される。
【０１４９】
以上の実施例を含む実施形態に関し、さらに以下の付記を開示する。
【０１５０】
（付記１）
物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップと、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップと、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納するステップと、
を含み、コンピュータにより実行される情報処理方法。
【０１５１】
（付記２）
前記物体は人体であり、
前記記憶装置に格納されている前記座標値を用いて、前記人体の状態を特定する状態特定ステップ
をさらに含む、付記１記載の情報処理方法。
【０１５２】
（付記３）
前記状態特定ステップが、
前記記憶装置に格納されている前記座標値により表される、前記水平面上の点列についての中心点を通る直線によって前記水平面を所定角度毎の領域に分割するステップと、
前記領域毎に、前記中心点からの距離が最大となり且つ前記座標値により表される点である第１の点を特定し、当該第１の点の座標値を前記中心点を原点とする場合の座標値に変換し、前記記憶装置に格納するステップと、
前記人体の状態の候補毎に当該人体の状態の候補に対応する、各前記領域内の点である第２の点の座標値を格納する状態データベースを用いて、前記記憶装置に格納されている前記第１の点の座標値と前記第２の点の座標値とから算出され且つ全ての前記領域に関する距離が最小となる前記人体の状態の候補を特定するステップと、
を含む、付記２記載の情報処理方法。
【０１５３】
（付記４）
前記物体は人であり、
前記補正ステップが、
前記人の歩行時の標準加速度を表すテンプレートを格納するテンプレート格納部に格納されている前記テンプレートを、歩行していたとみなされる期間において計測された前記加速度の値から特定される歩行周期及び加速度の振幅に合うように変形するステップと、
変形後の前記テンプレートを用いて、歩行していたとみなされる期間において計測された前記加速度の値から前記歩行時の標準加速度を取り除くステップと、
を含む、付記１乃至３記載の情報処理方法。
【０１５４】
（付記５）
前記初速度算出ステップが、
前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値ａと、前記所定期間と前記所定のサンプリング間隔Δｔとから算出されるサンプル数Ｎとを用いて、
【数２９】
0005477238.tif 000030
（ｉはサンプル番号を示す）
に従って、変数Ｋの値を算出するステップと、
前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値ａと、前記変数Ｋの値と、前記所定のサンプリング間隔Δｔと、前記サンプル数Ｎとを用いて、
【数３０】
0005477238.tif 000031
（ｊ及びｌはサンプル番号を表す変数である）
に従って、前記物体の初速度ｖ0の値を算出するステップと、
を含む、付記１乃至４記載の情報処理方法。
【０１５５】
（付記６）
前記加速度センサにより前記所定のサンプリング間隔Δｔで前記加速度の値を計測し、計測結果格納部に順次格納するステップ
をさらに含み、
前記所定期間におけるサンプリング回数をＮとし、
前記初速度算出ステップが、
前記所定期間ごとに、
前記データ格納部から、処理対象の前記所定期間に生成されたＮ個の前記加速度の補正後の値ａi（ｉは、０からＮ−１までの整数）を順次読み出すステップと、
読み出されたＮ個の前記加速度の補正後の値ａiを、
【数３１】
0005477238.tif 000032
に代入することで、処理対象の前記所定期間における変数Ｋの値を算出し、前記記憶装置に格納するステップと、
読み出されたＮ個の前記加速度の補正後の値ａi及び前記記憶装置に格納された前記変数Ｋの値を、
【数３２】
0005477238.tif 000033
（ｊ及びｌは、サンプル番号を表す変数である）
に代入することで、処理対象の前記所定期間における前記物体の初速度ｖ0の値を算出し、前記記憶装置に格納するステップと、
を含む、付記１乃至４記載の情報処理方法。
【０１５６】
（付記７）
物体に備えられた加速度センサによって所定のサンプリング間隔で計測された加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正ステップと、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出ステップと、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納するステップと、
をコンピュータに実行させるためのプログラム。
【０１５７】
（付記８）
所定のサンプリング間隔で計測された物体の加速度の値に対して所定の補正を行い、加速度の補正後の値をデータ格納部に格納する補正部と、
前記物体の位置を加速度及び初速度から求めるための第１の関係式並びに前記物体の重心動揺で前記物体にかかる力をバネモデルの運動方程式で表すことにより得られる第２の関係式を所定期間内の加速度の離散値と前記所定のサンプリング間隔とから初速度を得られるように変形した第３の関係式に、前記データ格納部に格納されている、前記所定期間内の前記加速度の補正後の値を代入して、前記物体の初速度の値を算出し、記憶装置に格納する初速度算出部と、
前記第１の関係式に、前記データ格納部に格納されている前記加速度の補正後の値と前記記憶装置に格納されている前記初速度の値とを代入して、前記物体の重心位置を表す、水平面における座標値を算出し、前記記憶装置に格納する重心位置算出部と、
を有する情報処理装置。
【０１５８】
（付記９）
前記所定のサンプリング間隔で前記物体の前記加速度の値を計測し、計測結果格納部に順次格納する加速度センサ
をさらに有し、
前記補正部は、前記計測結果格納部から前記所定期間において計測された前記物体の前記加速度の値を読み出して、前記所定の補正を行う
付記８記載の情報処理装置。
【０１５９】
（付記１０）
前記物体は人体であり、
前記記憶装置に格納されている前記座標値により表される、前記水平面上の点列についての中心点を通る直線によって前記水平面を所定角度毎の領域に分割し、前記領域毎に、前記中心点からの距離が最大となり且つ前記座標値により表される点である第１の点を特定し、当該第１の点の座標値を前記中心点を原点とする場合の座標値に変換し、前記記憶装置に格納し、前記人体の状態の候補毎に当該人体の状態の候補に対応する、各前記領域内の点である第２の点の座標値を格納する状態データベースを用いて、前記記憶装置に格納されている前記第１の点の座標値と前記第２の点の座標値とから算出され且つ全ての前記領域に関する距離が最小となる前記人体の状態の候補を特定する状態特定部
をさらに有する付記８又は９記載の情報処理装置。
【符号の説明】
【０１６０】
１  携帯端末
１０  サーバ    ２０  センサ付き端末
３０  車載装置    ５０  ネットワーク
１０１  加速度センサ    １０３  計測結果格納部
１０５  補正部    １０７  データ格納部
１０９  データ読出部    １１１  定数格納部
１１３  変数値算出部    １１５  変数値格納部
１１７  初速度算出部    １１９  初速度格納部
１２１  重心動揺計算部    １２３  重心動揺データ格納部
１２５  体調ＤＢ    １２７  体調推定部
１２９  体調履歴データ格納部    １３１  出力部
１４１  受信部    １４３  送信部
２０１  加速度センサ    ２０３  近接検知部    ２０５  送信部
３０１  受信部    ３０３  近接検知部    ３０５  出力部
１０５１  オフセット補正部    １０５２  重力方向補正部
１０５３  進行方向補正部    １０５４  歩行データ処理部
１０５５  パラメータ格納部    １０５６  区間データ格納部
１０５７  テンプレート格納部
【図１】
0005477238.tif 000036
【図２】
0005477238.tif 000037
【図３】
0005477238.tif 000038
【図４】
0005477238.tif 000039
【図５】
0005477238.tif 000040
【図６】
0005477238.tif 000041
【図７】
0005477238.tif 000042
【図８】
0005477238.tif 000043
【図９】
0005477238.tif 000044
【図１０】
0005477238.tif 000045
【図１１】
0005477238.tif 000046
【図１２】
0005477238.tif 000047
【図１３】
0005477238.tif 000048
【図１４】
0005477238.tif 000049
【図１５】
0005477238.tif 000050
【図１６】
0005477238.tif 000051
【図１７】
0005477238.tif 000052
【図１８】
0005477238.tif 000053
【図１９】
0005477238.tif 000054
【図２０】
0005477238.tif 000055
【図２１】
0005477238.tif 000056
【図２２】
0005477238.tif 000057
【図２３】
0005477238.tif 000058
【図２４】
0005477238.tif 000059
【図２５】
0005477238.tif 000060
【図２６】
0005477238.tif 000061
【図２７】
0005477238.tif 000062
【図２８】
0005477238.tif 000063
【図２９】
0005477238.tif 000064
【図３０】
0005477238.tif 000065
【図３１】
0005477238.tif 000066
【図３２】
0005477238.tif 000067
【図３３】
0005477238.tif 000068
【図３４】
0005477238.tif 000069
【図３５】
0005477238.tif 000070