(57)【要約】【課題】人（被験者）の上体の姿勢の動きを適切に測定することができるシステムを提供する。【解決手段】上体運動計測システム１は、角速度センサ４及び加速度センサ５が各々搭載されてた複数の慣性センサユニット２を備える。複数の慣性センサユニット２は、被験者Ｐの上体の互いに異なる箇所に装着される。角速度センサ４及び加速度センサ５の検出出力を基に、各慣性センサユニット２の姿勢を推定し、さらに加速度を推定する。複数の慣性センサユニット２の推定された加速度を基に、被験者Ｐの上体の角加速度を推定する。【選択図】図１
2015186515.tif 000002
【特許請求の範囲】
【請求項１】
角速度を検出する角速度センサと加速度を検出する加速度センサとの組が各々搭載されており、被験者の上体の互いに異なる箇所に装着される複数の慣性センサユニットと、
前記被験者の上体に装着された各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力により示される角速度の検出値と加速度の検出値を基に、前記被験者の運動環境における各慣性センサユニットの姿勢を推定するセンサユニット姿勢推定手段と、
前記被験者の上体に装着された各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、前記センサユニット姿勢推定手段により推定された各慣性センサユニットの姿勢とを基に、前記被験者の運動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度を推定するセンサユニット加速度推定手段と、
前記複数の慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により推定された加速度の組を基に、前記グローバル座標系で見た前記被験者の上体の角加速度を推定する上体角加速度推定手段とを備えることを特徴とする上体運動計測システム。
【請求項２】
請求項１記載の上体運動計測システムにおいて、
前記複数の慣性センサユニットは、前記被験者の体幹軸に沿って該上体の上下方向に間隔を存するように該上体に装着される２つの慣性センサユニットを少なくとも含んでおり、
前記上体角加速度推定手段が推定する角加速度は、前記被験者のロール方向での上体の角加速度と前記被験者のピッチ方向での角加速度とのうちの少なくともいずれか一方を含むことを特徴とする上体運動計測システム。
【請求項３】
請求項２記載の上体運動計測システムにおいて、
前記２つの慣性センサユニットのうちの上側の慣性センサユニットは、前記被験者の上体のうち、該被験者の胸椎と腰椎との境目よりも上側の箇所に装着され、前記２つの慣性センサユニットのうちの下側の慣性センサユニットは、前記被験者の上体のうち、該被験者の胸椎と腰椎との境目よりも下側の箇所に装着されることを特徴とする上体運動計測システム。
【請求項４】
請求項１〜３のいずれか１項に記載の上体運動計測システムにおいて、
前記グローバル座標系で見た各慣性センサユニットの加速度、又は前記グローバル座標系で見た前記被験者の上体の角加速度を対象状態量と定義したとき、前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記複数の慣性センサユニットを上体に装着した前記被験者の歩行動作が行われた場合に、該歩行動作における前記対象状態量の推定値の時系列を生成する基本推定値生成処理と、該歩行動作における２歩分の歩容での前記対象状態量の平均波形データを生成する平均波形データ生成処理とを実行するように構成されており、
前記平均波形データ生成処理は、前記被験者の歩行動作に含まれるｎ個（ｎ：２以上の整数）の２歩分の歩容のそれぞれに対応して前記基本推定値生成処理により生成された前記対象状態量の推定値の時系列により示される波形データを時間軸方向にスケール変換することにより、２歩分の歩容の期間の時間幅を正規化してなる正規化波形データを生成し、さらに、前記ｎ個の２歩分の歩容のそれぞれに対応するｎ個の前記正規化波形データの平均の波形データを前記平均波形データとして生成するように構成されていることを特徴とする上体運動計測システム。
【請求項５】
請求項４記載の上体運動計測システムにおいて、
前記平均波形データ生成処理のうちの前記正規化波形データを生成する処理は、該正規化波形データにおける前記被験者の右脚側の１歩の期間の時間幅と左脚側の１歩の期間の時間幅との比率が、前記被験者の歩行動作に含まれるｎ個の２歩分の歩容における右脚側の１歩の期間の実際の時間幅の平均値と、該ｎ個の２歩分の歩容における左脚側の１歩の期間の実際の時間幅の平均値との比率に一致するように、前記正規化波形データを生成するように構成されていることを特徴とする上体運動計測システム。
【請求項６】
請求項４又は５記載の上体運動計測システムにおいて、
前記平均波形データ生成処理を実行する前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記平均波形データ生成処理により生成した前記平均波形データからオフセット成分を除去するオフセット成分除去処理をさらに実行するように構成されており、該オフセット成分除去処理は、前記オフセット成分を除去した後の平均波形データにより示される前記対象状態量の値の、前記２歩分の歩容の期間における平均値がゼロになるという条件を満たすように前記オフセット成分を除去するように構成されていることを特徴とする上体運動計測システム。
【請求項７】
請求項４〜６のいずれか１項に記載の上体運動計測システムにおいて、
前記平均波形データ生成処理を実行する前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記複数の慣性センサユニットのうちの少なくとも１つの慣性センサユニットの加速度センサの検出出力により示される加速度の検出値又は当該少なくとも１つの慣性センサユニットに関して前記センサユニット加速度推定手段により推定された加速度の変化に基づいて、前記被験者の歩行動作における１歩毎の切替わりのタイミングを認識するように構成されていることを特徴とする上体運動計測システム。
【請求項８】
請求項４〜７のいずれか１項に記載の上体運動計測システムにおいて、
前記センサユニット加速度推定手段は、前記各慣性センサユニットの加速度に関する前記平均波形データを生成するように構成されており、
前記上体角加速度推定手段は、前記複数の慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により生成された平均波形データに示される各慣性センサユニットの加速度の推定値を用いて前記２歩分の歩容の期間における前記被験者の上体の角加速度の推定値の時系列により構成される角加速度波形データを生成するように構成されており、
前記上体角加速度推定手段により生成された前記角加速度波形データにより示される角加速度の推定値を積分することにより前記被験者の上体の角速度の推定値を算出する上体角速度推定手段をさらに備えており、
前記上体角速度推定手段は、前記２歩分の歩容の期間における前記被験者の上体の角速度の推定値の平均値がゼロになるという条件を満たすように、前記被験者の上体の角速度の推定値を算出するように構成されていることを特徴とする上体運動計測システム。
【請求項９】
請求項４〜７のいずれか１項に記載の上体運動計測システムにおいて、
前記上体角加速度推定手段は、前記被験者の上体の角加速度に関する前記平均波形データを生成するように構成されており、
前記上体角加速度推定手段により生成された前記平均波形データにより示される角加速度の推定値を積分することにより前記被験者の上体の角速度の推定値を算出する上体角速度推定手段をさらに備えており、
前記上体角速度推定手段は、前記２歩分の歩容の期間における前記被験者の上体の角速度の推定値の平均値がゼロになるという条件を満たすように、前記被験者の上体の角速度の推定値を算出するように構成されていることを特徴とする上体運動計測システム。
【請求項１０】
角速度を検出する角速度センサと加速度を検出する加速度センサとの組が各々搭載された複数の慣性センサユニットを被験者の上体の互いに異なる箇所に装着した状態で各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力を取得する第１工程と、
前記第１工程で取得された各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力により示される角速度の検出値と加速度の検出値とを基に、前記被験者の運動環境における各慣性センサユニットの姿勢を推定する第２工程と、
前記第１工程で取得された各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、各慣性センサユニットに関して前記第２工程で推定され各慣性センサユニットの姿勢とを基に、前記被験者の移動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度を推定する第３工程と、
前記複数の慣性センサユニットのそれぞれに関して前記３工程で推定された加速度の組を基に、前記グローバル座標系で見た前記被験者の上体の角加速度を推定する第４工程とを備えることを特徴とする上体運動計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人（被験者）の上体の姿勢の動きを推定するシステム及び方法に関する。
【背景技術】
【０００２】
人の歩行状態を把握するために、歩行中の人の所定の部位の姿勢、あるいは、その姿勢の変化等をジャイロセンサ等の角速度センサ用いて計測することが従来より一般に行われている。
【０００３】
例えば特許文献１には、歩行者の腰部と各脚の大腿部、下腿部及び足部とにジャイロセンサを取り付けて、これらのジャイロセンサの検出値を用いて、股関節、膝関節、足首関節の関節角度を測定する技術が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１２−６５７２３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、本願発明者の種々様々な実験、検討によれば、歩行中の人の上体の姿勢の動き、特に、ピッチ方向もしくはロール方向での上体の姿勢の動きは、歩行中の人の脚の動きの影響が顕著に反映される。
【０００６】
また、特に、上体の角加速度は、歩行中の人に作用する床反力の状態にも密接に関連する。
【０００７】
従って、歩行中の人の上体の角加速度等、該上体の姿勢の動きを観測することは、人の歩行状態（例えば、左右の脚の動きのバランス状態、リハビリによる脚の動きの回復状態等）を把握し、もしくは評価する上で必要性の高い事項である。
【０００８】
この場合、人の上体にジャイロセンサ等の角速度センサを取り付けて、該角速度センサによる角速度の検出値に基づいて、上体の姿勢の動きを観測することが考えられる。そして、この場合、前記特許文献１に見られる腰部のジャイロセンサと同様に、角速度センサを、幅広のベルト等の固定板を介して人の上体に取り付けることが考えられる。
【０００９】
しかしながら、このように角速度センサを、幅広のベルト等の固定板を介して人の上体に取り付ける手法では、比較的大きな接触面で人の上体に固定板が装着されるため、人に不快感、違和感、異物の装着感等を及ぼしやすい。
【００１０】
さらに、上記不快感、違和感、異物の装着感等に起因して、あるいは、上記固定板の重量に起因して、人の歩行形態が、上記固定板が装着されていない通常時の歩行形態と異なるものとなってしまう場合もある。そして、このような場合には、人の歩行状態を高い信頼性で適切に評価することができなくなってしまう。
【００１１】
かかる不都合を解消するためには、角速度センサを、極力小さい接触面積で人の上体に装着することが考えられる。
【００１２】
しかるに、このようにした場合には、角速度センサの姿勢を人の上体の姿勢に精度よく安定に追従させることが難しく、角速度センサの姿勢の動きが、上体の実際の姿勢の動きに対してずれを生じやすい。その結果、人の上体の姿勢の動きを高い信頼性で観測することが困難となる。
【００１３】
また、例えば人の上体の角加速度を、該上体に装着した角速度センサを用いて計測しようとする場合には、該角速度センサによる角速度の検出値を微分する演算が必要となる。このため、上体の角加速度の計測値が、角速度の検出値に含まれるノイズ成分の影響を受けやすい。ひいては、上体の角加速度の測定精度が低下しやすいという不都合がある。
【００１４】
本発明はかかる背景に鑑みてなされたものであり、人（被験者）の上体の姿勢の動きを適切に測定することができるシステム及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１５】
本発明の上体運動計測システムは、かかる目的を達成するために、角速度を検出する角速度センサと加速度を検出する加速度センサとの組が各々搭載されており、被験者の上体の互いに異なる箇所に装着される複数の慣性センサユニットと、
前記被験者の上体に装着された各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力により示される角速度の検出値と加速度の検出値を基に、前記被験者の運動環境における各慣性センサユニットの姿勢を推定するセンサユニット姿勢推定手段と、
前記被験者の上体に装着された各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、前記センサユニット姿勢推定手段により推定された各慣性センサユニットの姿勢とを基に、前記被験者の運動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度を推定するセンサユニット加速度推定手段と、
前記複数の慣性センサユニットのうちの少なくとも２つの慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により推定された加速度を線形結合することにより、前記グローバル座標系で見た前記被験者の上体の角加速度を推定する上体角加速度推定手段とを備えることを特徴とする（第１発明）。
【００１６】
なお、本発明において、各慣性センサユニットの角速度センサが検出する角速度は、より詳しくは、空間的な（３次元空間の）角速度ベクトルである。同様に、各慣性センサユニットの加速度センサが検出する加速度は、より詳しくは、空間的な（３次元空間の）加速度ベクトルである。
【００１７】
上記第１発明によれば、前記複数の慣性センサユニットが装着された前記被験者の運動環境における各慣性センサユニットの姿勢（空間的な姿勢）が、前記センサユニット姿勢推定手段によって推定される。
【００１８】
ここで、前記各慣性センサユニットは、前記角速度センサ及び加速度センサの組が搭載されたものである。このため、各慣性センサユニットの姿勢を、該慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力により示される角速度の検出値と加速度の検出値とを基に推定できる。その推定の演算処理としては、例えばストラップダウン演算等の公知の演算処理を採用することができる。
【００１９】
そして、被験者の移動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度（並進加速度）が前記センサユニット加速度推定手段により推定される。
【００２０】
この場合、各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値、該加速度センサ、又は該慣性センサユニットに対して固定されたローカル座標系で表現され得る並進加速度である。そして、被験者の移動環境における各慣性センサユニットの姿勢は、前記センサユニット姿勢手段手段により推定される。
【００２１】
従って、各慣性センサユニット毎に、各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、前記センサユニット姿勢推定手段により推定された該慣性センサユニットの姿勢とを基に、グローバル座標系で見た各慣性センサユニットの加速度（並進加速度）を推定できる。
【００２２】
そして、前記上体角加速度推定手段により被験者の上体の角加速度が推定される。
【００２３】
ここで、上体の姿勢の動きによる角加速度が発生した場合、該角加速度が一定であっても、上体の互いに離れた各部の加速度（並進加速度）は一般には互いに異なる。例えば、上体の姿勢の回転中心付近からより遠い部分の加速度は、該回転中心により近い部分の加速度よりも大きくなる。
【００２４】
従って、被験者の歩行時等の上体の動きによって該上体の角加速度が変化すると、被験者の上体の互いに異なる箇所に装着された複数の慣性センサユニットのそれぞれの加速度の値の組が変化する。
【００２５】
この場合、各慣性センサユニットの前記被験者の上体への装着箇所が局所的なものであっても、各慣性センサユニットの加速度（並進加速度）が、該慣性センサユニットの装着箇所での上体の部位の加速度に対して高い追従性を有するように各慣性センサユニットを被験者の上体に装着し得る。
【００２６】
従って、前記複数の慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により推定された加速度の組を基に、前記上体角加速度推定手段によって、前記グローバル座標系で見た前記被験者の上体の角加速度を推定することで、実際の角加速度に対して追従性の高い推定値を得ることができる。
【００２７】
また、被験者の上体への各慣性センサユニットの装着箇所を局所的なものとすることができるので、被験者に違和感や不快感等を及ぼすことなく、各慣性センサユニットを被験者の上体に装着することができる。このため、複数の慣性センサユニットを被験者に装着した場合と装着していない場合とで、被験者の動きに違いが生じるのを防止できる。ひいては、複数の慣性センサユニットを装着していない場合の被験者の通常的な動きに対して、上体の角加速度を計測することができる。
【００２８】
よって、第１発明によれば、人（被験者）の上体の姿勢の動きを適切に測定することができる。
【００２９】
なお、本発明では、前記上体角加速度推定手段は、例えば、前記複数の慣性センサユニットの加速度の推定値を線形結合することよって被験者の上体の角加速度の推定値を算出することができる。この場合、当該線形結合の係数値としては、慣性センサユニット同士の間の距離又は位置関係に応じて設定した値、あるいは、実験的に設定した値等を使用できる。
【００３０】
また、前記上体角加速度推定手段が推定する角加速度は、特定の１方向、又は２方向の角加速度であってもよい。
【００３１】
上記第１発明では、前記複数の慣性センサユニットは、前記被験者の体幹軸に沿って該上体の上下方向に間隔を存するように該上体に装着される２つの慣性センサユニットを少なくとも含んでおり、前記上体角加速度推定手段が推定する角加速度は、前記被験者のロール方向での上体の角加速度と前記被験者のピッチ方向での角加速度とのうちの少なくともいずれか一方を含むことが好ましい（第２発明）。
【００３２】
この第２発明によれば、前記２つの慣性センサユニットが、前記被験者の上体の体幹軸に沿って該上体の上下方向に間隔を存するように該上体に装着されるので、該２つの慣性センサユニットの加速度（グローバル座標系見た加速度）の組が前記被験者のロール方向での上体の角加速度と前記被験者のピッチ方向での角加速度とに対して高い相関性を有するものとなる。
【００３３】
より詳しくは、２つの慣性センサユニットの加速度のうち、主に、被験者のピッチ軸方向（略左右方向）の加速度の値の組が被験者の上体のロール方向の角加速度に対して高い相関性を有するものとなる。
【００３４】
また、２つの慣性センサユニットの加速度のうち、主に、被験者のロール軸方向（略前後方向）の加速度の値の組が被験者の上体のピッチ方向の角加速度に対して高い相関性を有するものとなる。
【００３５】
従って、第２発明によれば、被験者のピッチ方向及びロール方向の一方又は両方の角加速度を高い信頼性で推定できる。ひいては、被験者の歩行動作等の運動状態を解析もしくは評価する上で有用性の高い計測データ（角加速度の推定値）を得ることができる。
【００３６】
上記第２発明では、前記２つの慣性センサユニットのうちの上側の慣性センサユニットは、前記被験者の上体のうち、該被験者の胸椎と腰椎との境目よりも上側の箇所に装着され、前記２つの慣性センサユニットのうちの下側の慣性センサユニットは、前記被験者の上体のうち、該被験者の胸椎と腰椎との境目よりも下側の箇所に装着されることことが好ましい（第３発明）。
【００３７】
ここで、前記被験者の上体のうち、該被験者の胸椎と腰椎との境目よりも上側の箇所は、内部に肋骨が配置されている箇所なので該上体において、比較的剛性が高い箇所である。
【００３８】
また、上記境目より下側の箇所は、被験者の上体のピッチ方向又はロール方向の姿勢の回転中心近辺の箇所に相当する。
【００３９】
従って、第３発明によれば、被験者のピッチ方向及びロール方向の一方又は両方の角加速度を、より一層、高い信頼性で推定できる。
【００４０】
また、前記第１〜第３発明では、前記グローバル座標系で見た各慣性センサユニットの加速度、又は前記グローバル座標系で見た前記被験者の上体の角加速度を対象状態量と定義したとき、前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記複数の慣性センサユニットを上体に装着した前記被験者の歩行動作が行われた場合に、該歩行動作における前記対象状態量の推定値の時系列を生成する基本推定値生成処理と、該歩行動作における２歩分の歩容での前記対象状態量の平均波形データを生成する平均波形データ生成処理とを実行するように構成されており、前記平均波形データ生成処理は、前記被験者の歩行動作に含まれるｎ個（ｎ：２以上の整数）の２歩分の歩容のそれぞれに対応して前記基本推定値生成処理により生成された前記対象状態量の推定値の時系列により示される波形データを時間軸方向にスケール変換することにより、２歩分の歩容の期間の時間幅を正規化してなる正規化波形データを生成し、さらに、前記ｎ個の２歩分の歩容のそれぞれに対応するｎ個の前記正規化波形データの平均の波形データを前記平均波形データとして生成するように構成されていることが好ましい（第４発明）。
【００４１】
なお、第４発明では、前記センサユニット加速度推定手段が、前記対象状態量としての各慣性センサユニットの加速度に関する前記基本推定値生成処理及び前記平均波形データ生成処理を実行する場合、前記基本推定値生成処理では、より詳しくは、各慣性センサユニット毎に、該慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、前記センサユニット姿勢推定手段により推定された該慣性センサユニットの姿勢とを基に、前記歩行動作における前記対象状態量（各慣性センサユニットの加速度）の推定値の時系列が算出される。
【００４２】
また、前記上体姿勢角加速度推定手段が、前記対象状態量としての被験者の上体の角か加速度に関する前記基本推定値生成処理及び前記平均波形データ生成処理を実行する場合、前記基本推定値生成処理では、より詳しくは、前記複数の慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により推定された加速度の組を基に、前記歩行動作における前記対象状態量（上体の角加速度）の推定値の時系列が算出される。
【００４３】
ここで、被験者の歩行状態での上体の動きを、解析あるいは評価する上では、被験者の複数歩の歩行動作における２歩分の歩容（換言すれば、歩行動作の繰り返しにおける１周期分の歩容）の期間での、平均的な上体の動き（加速度、各加速度等）の波形データを得ることが望ましいと考えられる。
【００４４】
そこで、第４発明では、前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記対象状態量の連続的な波形データを得るための前記基本推定値生成処理に加えて、前記平均波形データ生成処理を実行する。
【００４５】
この場合、前記被験者の歩行動作に含まれるｎ個の（複数の）２歩分の歩容のそれぞれの期間の時間幅は、定常的な歩行動作であっても、一般には、ｎ個の（複数の）２歩分の歩容のそれぞれに関して互いに異なる。
【００４６】
そこで、前記平均波形データ生成処理では、ｎ個の（複数の）２歩分の歩容のそれぞれについて、２歩分の歩容の期間の時間幅を正規化してなる前記正規化波形データを生成した上で、前記平均波形データを生成する。
【００４７】
なお、２歩分の歩容の期間の時間幅を正規化するということは、ｎ個の２歩分の歩容のそれぞれの時間幅を、互いに同一の所定値（一定値）にすることを意味する。
【００４８】
これにより、被験者の歩行動作における２歩分の歩容の期間における前記対象状態量の平均的な波形データとしての前記平均波形データを適切に生成することができる。
【００４９】
また、上記の如く生成される平均波形データを使用することで、被験者の歩行状態の解析あるいは評価を高い信頼性で適切に行うことが可能となる。
【００５０】
補足すると、本発明では、前記平均波形データにより示される前記対象状態量の値（瞬時値）は、該対象状態量の推定値であるとみなす。
【００５１】
従って、前記センサユニット加速度推定手段が、前記対象状態量としての各慣性センサユニットの加速度に関する前記基本推定値生成処理及び前記平均波形データ生成処理を実行する場合においては、各慣性センサユニット毎の前記平均波形データにより示される該慣性センサユニットの加速度の値は、前記センサユニット加速度生成手段が推定した加速度としての意味を持つ。
【００５２】
この場合、前記上体角加速度推定手段は、前記複数の慣性センサユニットのそれぞれに関する前記平均波形データにより示される加速度の値（推定値）を用いて、前記２歩分の歩容の期間における被験者の角加速度（瞬時値又はその時系列）を推定するようにしてもよい。
【００５３】
上記第４発明では、前記平均波形データ生成処理のうちの前記正規化波形データを生成する処理は、該正規化波形データにおける前記被験者の右脚側の１歩の期間の時間幅と左脚側の１歩の期間の時間幅との比率が、前記被験者の歩行動作に含まれるｎ個の２歩分の歩容における右脚側の１歩の期間の実際の時間幅の平均値と、該ｎ個の２歩分の歩容における左脚側の１歩の期間の実際の時間幅の平均値との比率に一致するように、前記正規化波形データを生成するように構成されていることが好ましい（第５発明）。
【００５４】
なお、前記被験者の右脚側の１歩の期間というのは、右脚が支持脚（又は遊脚）として動作する期間を意味し、前記被験者の左脚側の１歩の期間というのは、左脚が支持脚（又は遊脚）として動作する期間を意味する。この場合、両脚支持期では便宜上、例えば、該両脚支持期の直後に離床する脚を遊脚、他方の脚を支持脚と定義することができる。あるいは、逆に、両脚支持期の直後に離床する脚を支持脚、他方の脚を遊脚と定義してもよい。
【００５５】
上記第５発明よれば、前記ｎ個の２歩分の歩容のそれぞれに対応するｎ個の正規化波形データのそれぞれにおける右脚側の１歩の期間と左脚側の１歩の期間との切替わりのタイミングは、ｎ個の正規化波形データのいずれについても同じとなる。
【００５６】
従って、当該切替わりのタイミングの近辺での対象状態量の平均的な波形データを高い信頼性で適切に得ることができる。また、平均波形データによって、右脚側の１歩の期間での対象状態量の波形と左脚側の１歩の期間での対象状態量の波形との区別又はその比較を適切に行うことも可能となる。
【００５７】
前記第４発明又は第５発明では、前記平均波形データ生成処理を実行する前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記平均波形データ生成処理により生成した前記平均波形データからオフセット成分を除去するオフセット成分除去処理をさらに実行するように構成されており、該オフセット成分除去処理は、前記オフセット成分を除去した後の平均波形データにより示される前記対象状態量の値の、前記２歩分の歩容の期間における平均値がゼロになるという条件を満たすように前記オフセット成分を除去するように構成されていることが好ましい（第６発明）。
【００５８】
すなわち、各慣性センサユニットの角速度センサ又は加速度センサの検出出力のドリフト等に起因して、一般には、前記平均波形データ（前記ｎ個の２歩分の歩容のそれぞれに対応するｎ個の前記正規化波形データの平均の波形データ）は、オフセット成分を含む。
【００５９】
一方、被験者の定常的な歩行動作もしくはそれに近い歩行動作では、前記２歩分の歩容の期間における前記対象状態量の実際の値の平均値は、ゼロもしくはほぼゼロとなる。
【００６０】
従って、前記平均波形データからオフセット成分を適切に除去することができる。その結果、前記平均波形データにおける対象状態量の値の極性と、実際の極性との整合性を高めることができる。従って、平均波形データの信頼性をより一層高めることができる。
【００６１】
上記第４〜第６発明では、前記平均波形データを好適に生成するためには、２歩分の歩容の期間の時間幅、あるいは、１歩の期間の時間幅を特定するために、前記被験者の歩行動作における１歩毎の切替わりのタイミング（右脚及び左脚の一方側の１歩から他方側の一歩への切替わりのタイミング）を的確に認識する必要がある。
【００６２】
この場合、例えば、被験者の両脚のそれぞれの着地又は離床を検知する接地センサを両脚の足底部に備えることで、上記の認識を行うことができる。
【００６３】
また、本願発明者の各種実験、検討によれば、前記複数の慣性センサユニットのうちの少なくとも１つの慣性センサユニット（例えば前記第２発明又は第３発明における下側の慣性センサユニット）は、その加速度が、１歩毎の切替わりのタイミング付近で、比較的顕著な変化を呈するように配置できる。
【００６４】
そこで、第４〜第６発明では、前記平均波形データ生成処理を実行する前記センサユニット加速度推定手段又は前記上体角加速度推定手段は、前記複数の慣性センサユニットのうちの少なくとも１つの慣性センサユニットの加速度センサの検出出力により示される加速度の検出値又は当該少なくとも１つの慣性センサユニットに関して前記センサユニット加速度推定手段により推定された加速度の変化に基づいて、前記被験者の歩行動作における１歩毎の切替わりのタイミングを認識するように構成されていてもよい（第７発明）。
【００６５】
この第７発明によれば、上記接地センサ等を備えずとも、前記被験者の歩行動作における１歩毎の切替わりのタイミングを適切に認識できる。ひいては、信頼性の高い平均波形データを得ることができる。
【００６６】
また、前記第４〜第７発明では、被験者の上体の角加速度の推定値の時系列を用いて該上体の角速度を推定することもできる。この場合、次のような構成を採用することが望ましい。
【００６７】
すなわち、前記センサユニット加速度推定手段が、前記各慣性センサユニットに加速度に関する前記平均波形データを生成するように構成されている場合には、前記上体角加速度推定手段は、前記少なくとも２つの慣性センサユニットのそれぞれに関して前記センサユニット加速度推定手段により生成された平均波形データに示される各慣性センサユニットの加速度の推定値を用いて前記２歩分の歩容の期間における前記被験者の上体の角加速度の推定値の時系列により構成される角加速度波形データを生成するように構成される。
【００６８】
そして、この場合、前記上体角加速度推定手段により生成された前記角加速度波形データにより示される角加速度の推定値を積分することにより前記被験者の上体の角速度の推定値を算出する上体角速度推定手段をさらに備えており、前記上体角速度推定手段は、前記２歩分の歩容の期間における前記被験者の上体の角速度の推定値の平均値がゼロになるという条件を満たすように、前記被験者の上体の角速度の推定値を算出するように構成されていることが好ましい（第８発明）。
【００６９】
あるいは、前記上体角加速度推定手段が、前記被験者の上体の角加速度に関する前記平均波形データを生成するように構成されている場合には、前記上体角加速度推定手段により生成された前記平均波形データにより示される角加速度の推定値を積分することにより前記被験者の上体の角速度の推定値を算出する上体角速度推定手段をさらに備えており、前記上体角速度推定手段は、前記２歩分の歩容の期間における前記被験者の上体の角速度の推定値の平均値がゼロになるという条件を満たすように、前記被験者の上体の角速度の推定値を算出するように構成されていることが好ましい（第９発明）。
【００７０】
ここで、被験者の定常的な歩行動作もしくはそれに近い歩行動作では、前記２歩分の歩容の期間における被験者の上体の角速度の実際の値の平均値は、ゼロもしくはほぼゼロとなる。
【００７１】
そこで、第８発明又は第９発明では、前記上体角速度推定手段は、前記２歩分の歩容の期間における前記被験者の上体の角速度の推定値の平均値がゼロになるという条件を満たすように、前記被験者の上体の角速度の推定値を算出する。
【００７２】
これにより、前記２歩分の歩容の期間における被験者の平均的な上体の角速度の推定値（瞬時値又は時系列）を高い信頼性で適切に算出することができる。
【００７３】
なお、上記第８発明又は第９発明を前記第６発明と組み合わせた場合には、前記平均波形データは、前記オフセット成分を除去した平均波形データを意味する。
【００７４】
また、本発明の上体運動計測方法は、角速度を検出する角速度センサと加速度を検出する加速度センサとの組が各々搭載された複数の慣性センサユニットを被験者の上体の互いに異なる箇所に装着した状態で各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力を取得する第１工程と、
前記第１工程で取得された各慣性センサユニットの角速度センサ及び加速度センサのそれぞれの検出出力により示される角速度の検出値と加速度の検出値とを基に、前記被験者の運動環境における各慣性センサユニットの姿勢を推定する第２工程と、
前記第１工程で取得された各慣性センサユニットの加速度センサの検出出力により示される加速度の検出値と、各慣性センサユニットに関して前記第２工程で推定され各慣性センサユニットの姿勢とを基に、前記被験者の移動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度を推定する第３工程と、
前記複数の慣性センサユニットのそれぞれに関して前記３工程で推定された加速度の組を基に、前記グローバル座標系で見た前記被験者の上体の角加速度を推定する第４工程とを備えることを特徴とする（第１０発明）。
【００７５】
上記第１０発明によれば、前記複数の慣性センサユニットが装着された前記被験者の運動環境における各慣性センサユニットの姿勢（空間的な姿勢）が前記第２工程で推定される。この指定処理は、第１発明におけるセンサユニット姿勢推定手段の処理の場合と同様に行われる。
【００７６】
そして、被験者の移動環境に設定されたグローバル座標系で見た各慣性センサユニットの加速度（並進加速度）が第３工程で推定される。この推定処理は、第１発明におけるセンサユニット加速度推定手段の処理と同様に行わる。
【００７７】
そして、前記第４工程で、被験者の上体の角加速度が推定される。この推定処理は、第１発明における上体角加速度推定手段の処理と同様に行われる。
【００７８】
従って、前記複数の慣性センサユニットのそれぞれに関して推定された加速度の組を基に、前記グローバル座標系で見た前記被験者の上体の角加速度を推定することで、実際の角加速度に対して追従性の高い推定値を得ることができる。
【００７９】
また、第１発明と同様に、ｌ被験者の上体への各慣性センサユニットの装着箇所を局所的なものとすることができるので、複数の慣性センサユニットを被験者に装着した場合と装着していない場合とで、被験者の動きに違いが生じるのを防止できる。ひいては、複数の慣性センサユニットを装着していない場合の被験者の通常的な動きに対して、上体の角加速度を計測することができる。
【００８０】
よって、第１０発明によれば、人（被験者）の上体の姿勢の動きを適切に測定することができる。
【図面の簡単な説明】
【００８１】
【図１】（ａ）は第１実施形態の上体運動計測システムの全体構成を示す図、（ｂ）は被験者の上体に対する慣性センサユニットの配置形態を例示する図。
【図２】実施形態の各慣性センサユニットの構成を示す斜視図。
【図３】実施形態の上体運動計測システムの動作に関する構成を示すブロック図。
【図４】図１に示すデータ収集装置の処理を示すフローチャート。
【図５】（ａ）〜（ｃ）は、図４のＳＴＥＰ５，６の処理に関する説明図。
【図６】実施形態で得られたデータの表示例を示す図。
【図７】（ａ），（ｂ）は実施形態で得られるデータの評価処理の一例を説明するためのグラフ。
【図８】（ａ），（ｂ）は、第２実施形態において被験者の上体に装着する慣性センサユニットの配置形態の例を示す図。
【図９】被験者の上体に装着する慣性センサユニットの配置形態の他の例を示す図。
【図１０】本発明の実施形態に係る変形態様を説明するためのグラフ。
【図１１】本発明の実施形態に係る他の変形態様を説明するためのグラフ。
【発明を実施するための形態】
【００８２】
［第１実施形態］
本発明の第１実施形態を図１〜図８を参照して以下に説明する。
【００８３】
図１（ａ），（ｂ）を参照して、本実施形態の上体運動計測システム１は、被験者Ｐに装着される２つの慣性センサユニット２，２と、これらの慣性センサユニット２，２との通信を行いつつ測定データの収集等を行うデータ収集装置１０とを備える。
【００８４】
各慣性センサユニット２（以降、単にセンサユニット２という）は、図２に示すように、筐体３内に収容された角速度センサ４及び加速度センサ５と、通信機能等を有する回路ユニット６とを備える。角速度センサ４、加速度センサ５及び回路ユニット６は、筐体３に固定されている。
【００８５】
角速度センサ４は、空間的な（３次元空間の）角速度ベクトルを検出するセンサである。この角速度センサ４は、３つの検出軸方向のそれぞれの角速度（角速度ベクトルの各検出軸方向の成分）の検出出力を発生する３軸の角速度センサにより構成される。
【００８６】
加速度センサ５は、空間的な（３次元空間の）並進加速度ベクトルを検出するセンサである。この加速度センサ５は、３つの検出軸方向のそれぞれの加速度（並進加速度ベクトルの各検出軸方向の成分）の検出出力を発生する３軸の加速度センサにより構成される。
【００８７】
なお、３軸の角速度センサ４は、互いに検出軸方向が異なる各別の３つの１軸角度センサを組み合わせたもの、あるいは、３つの１軸角度センサを一体に構成したもののいずれであってもよい。このことは、３軸の加速度センサ５についても同様である。
【００８８】
本実施形態では、各センサユニット２毎の角速度センサ４の３つの検出軸方向と加速度センサ５の３つの検出軸方向とはいずれも互いに直交する方向である。また、角速度センサ４及び加速度センサ５は、それぞれの３つの検出軸方向が角速度センサ４と加速度センサ５とで互いに同じになるようにセンサユニット２に搭載されている。
【００８９】
例えば、図２に示すように、センサユニット２に対して固定的に設定されたｘｙｚ直交座標系Ｃｓ（以降、センサユニットローカル座標系Ｃｓという）の各軸方向（ｘ軸方向、ｙ軸方向及びｚ軸方向）が角速度センサ４及び加速度センサ５の共通の検出軸方向である。
【００９０】
ただし、各センサユニット２毎の角速度センサ４の３つの検出軸方向と加速度センサ５の３つの検出軸方向とは互いに異なっていてもよい。さらには、角速度センサ４の３つの検出軸方向、又は、加速度センサ５の３つの検出軸方向は、２つ以上の検出軸方向が互いに平行となる場合を除いて、互いに直交する方向でなくてもよい。
【００９１】
この場合、角速度センサ４及び加速度センサ５の相互の配置関係が一定に保たれるので、角速度センサ４による角速度ベクトルの検出値又は加速度センサ５による加速度ベクトルの検出値から、座標変換によって、角速度センサ４と加速度センサ５とで同一の検出軸方向の検出値を得ることができる。
【００９２】
各センサユニット２の回路ユニット６は、図３に示すように、演算処理部７及び通信機能部８を備える。
【００９３】
演算処理部７は、ＣＰＵ、メモリ等を含む回路部である。該演算処理部７は、実装されるプログラムを実行することで実現される機能として、該演算処理部７を含む回路ユニット６が搭載されたセンサユニット２の空間的な姿勢（空間的な向き）を推定するセンサユニット姿勢推定手段としての機能を有する。
【００９４】
この機能を有する演算処理部７は、角速度センサ４及び加速度センサ５のそれぞれの検出出力を増幅器、Ａ／Ｄ変換器等の入力回路（図示省略）を介して取得する。そして、演算処理部７は、角速度センサ４の検出出力により示される角速度ベクトルの検出値と、加速度センサ５の検出出力により示される加速度ベクトルの検出値とを用いて、公知のストラップダウン演算等、所定の演算処理を実行することで、被験者Ｐの運動環境におけるセンサユニット２の空間的な姿勢を推定する。
【００９５】
演算処理部７で推定される各センサユニット２の姿勢は、図１（ａ），（ｂ）に例示する如く、被験者Ｐの移動環境に設定されるグローバル座標系Ｃｇで見た姿勢として表現される。このグローバル座標系Ｃｇとしては、例えば、前記センサユニットローカル座標系Ｃｓのｘ軸を水平面に投影した軸の方向をＸ軸方向、鉛直方向（重力方向）をＺ軸方向、Ｘ軸及びＺ軸に直交する方向をＹ軸方向としたＸＹＺ直交座標系が使用される。
【００９６】
通信機能部８は、無線通信又は有線通信用の通信機器及びそれに付随する通信回路により構成される。この通信機能部８は、演算処理部７で推定されたセンサユニット２の姿勢を示す計測データと、加速度センサ５の検出出力により示されるセンサユニット２の加速度（並進加速度ベクトル）の検出値を示す計測データとを、無線又は有線によりデータ収集装置１０に送信する。
【００９７】
上記の如く構成された２つのセンサユニット２は、本実施形態では、図１（ａ），（ｂ）に示すように、被験者Ｐの上体の背中（背面）に、上下に間隔を有するように装着される。
【００９８】
具体的には、本実施形態では、２つのセンサユニット２，２は、被験者Ｐの上体の背中の脊椎沿いの線上で（体幹軸Ｌに沿って）で、第１２胸椎（胸椎のうちの最下段の椎骨）と第１腰椎（腰椎のうちの最上段の椎骨）との境界（以降、胸椎／腰椎境界という）よりも上側の箇所と、該胸椎／腰椎境界よりも下側の箇所（腰付近の箇所）に装着される。
【００９９】
より具体的には、上側のセンサユニット２の装着箇所は、例えば被験者Ｐの首の付け根付近の箇所（胸椎のうちの最上段の椎骨（第１胸椎）付近の箇所）である。また、下側のセンサユニット２の装着箇所は、例えば腰椎のうちの最下段の椎骨（第５腰椎）付近の箇所である。
【０１００】
ここで、人の上体のうち、腰椎が存在する部位は、該上体の姿勢をピッチ方向又はロール方向に変化させる場合における該上体の主たる折れ曲がり部分（回転中心部分）に相当する。従って、下側のセンサユニット２は、上体の姿勢の変化における回転中心付近の箇所で上体に装着されることとなる。
【０１０１】
また、胸椎には肋骨が接触するので、胸椎が存在する部位は、腰椎が存在する部位等に比して剛性が高い部位である。従って、上側のセンサユニット２は、上体の剛性が高い部位の上端寄りの箇所で上体に装着されることとなる。
【０１０２】
補足すると、センサユニット２，２の間の上下方向の距離が大きい方が、該距離が小さい場合よりも、上体の角速度の変化に対するセンサユニット２，２の加速度（Ｘ軸方向又はＹ軸方向の加速度）の差の変化の感度が高まる。このため、被験者Ｐの上体の角加速度（ピッチ方向の角加速度又はロール方向の角加速度）を後述する如く推定（計測）する場合、センサユニット２，２の上下方向の距離が大きい方が小さい場合よりも上体の角速度の推定の分解能を高めることができる。
【０１０３】
なお、センサユニット２，２の装着箇所は、上記の例に限られるものではない。例えば、センサユニット２，２の少なくともいずれか一方を被験者Ｐの前側（胸側あるいは腹側）に装着してもよい。また、センサユニット２，２を上記の例と異なる高さ位置で被験者Ｐに装着してもよい。
【０１０４】
被験者Ｐの上体の背中への各センサユニット２の装着は、適宜の装着用部材を介して行われる。例えば、被験者Ｐは、レオタード等、伸縮性を有すると共に被験者Ｐの上体への密着性が高い衣類を装着する。そして、各センサユニット２は、その筐体３が、例えば面ファスナー等の装着用部材を介して被験者Ｐの上体の衣類に固定される。これにより、各センサユニット２は、その装着箇所での上体の部位と一体又はほぼ一体に並進運動を行い得るように該上体に装着される。
【０１０５】
なお、各センサユニット２の被験者Ｐへの装着形態は、上記の形態に限らない。例えば各センサユニット２を、被験者Ｐの上体に装着するベルトに取り付けてもよい。あるいは、各センサユニット２を被験者Ｐの上体の皮膚にテープ等により直接的に貼り付けるようにしてもよい。
【０１０６】
この場合、各センサユニット２は、図２に示したセンサユニットローカル座標系Ｃｓのｘ軸方向が被験者Ｐの上体に対する該センサユニット２の接触面の法線方向（換言すれば、被験者Ｐの上体の前後方向）に一致もしくはほぼ一致し、ｙ軸方向が被験者Ｐの上体の左右方向（幅方向）に一致もしくはほぼ一致し、ｚ軸方向が被験者Ｐの上体の上下方向（体幹軸Ｌの方向）に一致もしくはほぼ一致するように被験者Ｐの上体に装着される。
【０１０７】
データ収集装置１０は、例えばディスプレイ１０ａを有するコンピュータにより構成される。このデータ収集装置１０は、各センサユニット２の回路ユニット６との通信によって、各センサユニット２の姿勢の推定値を示す信号を受信する。
【０１０８】
そして、データ収集装置１０は、実装されるプログラムにより実現される機能、あるいは、ハードウェア構成により実現される機能として、図３に示す如く、グローバル座標系Ｃｇで見た各慣性センサユニット２の加速度を推定するセンサユニット加速度推定部１１と、グローバル座標系で見た被験者Ｐの上体の角加速度を含む該上体の姿勢の動きを推定する上体姿勢運動推定部１２とを備える。
【０１０９】
上記センサユニット加速度推定部１１は、本発明におけるセンサユニット加速度推定手段としての機能を有するものである。また、上体姿勢運動推定部１２は、本発明における上体角加速度推定手段及び上体角速度推定手段としての機能を有するものである。
【０１１０】
本実施形態の上体運動計測システム１は、さらに、図１（ａ）に示す如く、被験者Ｐの各脚の着地及び離床を検知するための接地センサ２１，２１を備える。接地センサ２１は、感圧スイッチ等により構成され、被験者Ｐの各足平部に装着する靴の底部に装着されている。そして、各脚の接地センサ２１は、該脚の足平部が着地すると、オン信号を出力し、該脚の足平部が離床すると、オフ信号を出力する。この場合、各脚の接地センサ２１の出力信号は、図示を省略する通信装置を介して無線又は有線によりデータ収集装置１０に送信される。
【０１１１】
次に、データ収集装置１０の各機能部の処理の詳細と合わせて、本実施形態の上体運動計測システムによる被験者Ｐの上体の姿勢の動きの計測処理を説明する。
【０１１２】
まず、センサユニット２，２とデータ収集装置１０とが起動される。また、被験者Ｐの上体に前記した如く２つのセンサユニット２，２が装着されると共に、被験者Ｐがその各足平部に、前記接地センサ２１を有する靴を装着する。なお、このとき、センサユニット２，２とデータ収集装置１０との間でタイマー時刻の同期設定、あるいは、タイマー時刻の相互認識等の初期設定処理も行われる。
【０１１３】
この状態で、被験者Ｐの歩行（直進歩行）が行われる。その歩行は、２＊ｎ歩（ｎはあらかじめ定めた２以上の整数値）以上の歩数で行われる（“＊”は乗算記号を意味する）。なお、被験者Ｐは、適宜の歩行補助装具を装着していてもよい。
【０１１４】
この歩行中に、センサユニット姿勢推定手段としての機能を有する各センサユニット２の演算処理部７は、角速度センサ４及び加速度センサ５の検出出力を逐次取得しつつ、該センサユニット２の空間的な姿勢をストラップダウン演算等の手法によって逐次推定する。なお、センサユニット２の姿勢の推定値は、例えばオイラー角により表現される。ただし、センサユニット２の姿勢は、クオータニオン、あるいは、姿勢行列の形態で表現することもできる。
【０１１５】
具体的には、演算処理部７は、角速度センサ４の検出出力により示される角速度（３軸成分から成る角速度ベクトル）の検出値を積分する演算によって、センサユニット２の瞬時姿勢の基本値を逐次算出する。
【０１１６】
また、演算処理部７は、加速度センサ５の検出出力により示される加速度（３軸成分からなる並進加速度ベクトル）の検出値を基に、重力方向（鉛直方向）に対するセンサユニット２の瞬時姿勢（ピッチ方向及びロール方向の傾斜角）を逐次算出する。
【０１１７】
さらに、演算処理部７は、角速度センサ４の検出出力に基づく瞬時姿勢の基本値と、加速度センサ５の検出出力に基づく瞬時姿勢との偏差を算出する。
【０１１８】
ここで、加速度センサ５が検出する加速度は重力加速度ベクトルを含んでいる。また、加速度センサ５が検出する並進加速度ベクトルの水平方向成分の大きさは、一般に重力加速度ベクトルの大きさに比して十分に小さい。このため、加速度センサ５が検出する並進加速度ベクトルの長期的な平均値は、重力加速度ベクトルにほぼ一致する。
【０１１９】
そこで、演算処理部７は、上記偏差を長期的にゼロに収束させるように、角速度センサ４の検出出力に基づく瞬時姿勢の基本値を、該偏差に応じて補正することで、センサユニット２の姿勢の推定値を算出する。
【０１２０】
センサユニット２の姿勢の推定手法は、上記の如く、ストラップダウン演算等の公知の手法を採用できる。
【０１２１】
このように各センサユニット２の演算処理部７で逐次推定される該センサユニット２の姿勢を示す計測データが該センサユニット２の通信機能部８により、データ収集装置１０に送信される。さらに、該センサユニット２の加速度センサ５の検出出力により示される並進加速度ベクトルの検出値を示す計測データも、通信機能部８により、データ収集装置１０に送信される。
【０１２２】
なお、各センサユニット２の姿勢を推定する処理は、データ収集装置１０で実行するようにしてもよい。すなわち、データ収集装置１０にセンサユニット姿勢推定手段としての機能を持たせてもよい。この場合には、各センサユニット２の角速度センサ４及び加速度センサ５の検出出力を示す計測データを、各センサユニット２からデータ収集装置１０に送信する。
【０１２３】
補足すると、角速度センサ４及び加速度センサ５以外にセンサユニット２に搭載したカメラ等のセンサを使用して、センサユニット２の水平方向（床面とほぼ平行な方向）の加速度（並進加速度）を推定することも可能である。この場合には、加速度センサ５が検出する加速度ベクトルから、水平方向の加速度ベクトルの推定値を減算することで、センサユニットローカル座標系Ｃｓで見た重力加速度ベクトルを求めることができる。ひいては、この重力加速度ベクトルを基に、重力方向に対するセンサユニット２の瞬時姿勢を推定することもできる。
【０１２４】
また、各センサユニット２に地磁気センサを搭載した場合には、角速度センサ４の検出出力により示される角速度のうち、ヨー方向の成分のドリフト補正を行うこともできる。
【０１２５】
データ収集装置１０に各センサユニット２から送信される計測データ（詳しくは、被験者Ｐの歩行中の各時刻における各センサユニット２の姿勢の推定値及び並進加速度ベクトルの検出値を示す計測データ）は、該データ収集装置１０において、メモリ等の記憶媒体に時系列的に記憶保持される。
【０１２６】
そして、データ収集装置１０は、その計測データを用いて、前記被験者Ｐの上体の動きを解析又は評価するための処理（センサユニット加速度推定部１１及び上体姿勢運動推定部１２の処理を含む）を、図４のフローチャートに示す如く実行する。
【０１２７】
なお、以降の説明では、２つのセンサユニット２，２を区別するために、一方のセンサユニット２（例えば上側のセンサユニット２）を第１センサユニット２、他方のセンサユニット２（例えば下側のセンサユニット２）を第２センサユニット２ということがある。また、第１センサユニット２及び第２センサユニット２のうちの任意の１つを第ｉセンサユニット２（ｉ＝１，２）ということがある。
【０１２８】
また、第ｉセンサユニット２に対応するセンサユニットローカル座標系Ｃｓを第ｉローカル座標系Ｃs_iと称する。さらに、第ｉローカル座標系Ｃs_iのｘ軸を水平面に投影した軸の方向をＸ軸方向とするグローバル座標系Ｃｇを第ｉグローバル座標系Ｃg_iと称する。
【０１２９】
補足すると、各センサユニット２に対応する第ｉグローバル座標系Ｃg_iのＸ軸方向、Ｙ軸方向、Ｚ軸方向のそれぞれは、いずれのセンサユニット２についても互いに同一方向もしくはほぼ同一方向である。このため、第ｉグローバル座標系Ｃg_iを互いに区別せずに、単にグローバル座標系Ｃgと称する場合もある。
【０１３０】
なお、本実施形態の説明では、グローバル座標系ＣｇのＸ軸、Ｙ軸、Ｚ軸の正方向は、それぞれ、被験者Ｐから見て、前向き、左向き、上向きと定義する。また、センサユニットローカル座標系Ｃｓ（第ｉローカル座標系Ｃs_i）のｘ軸、ｙ軸、ｚ軸の正方向は、それぞれ、被験者Ｐの上体の前方向き、左向き、上向きと定義する。
【０１３１】
図４を参照して、データ収集装置１０は、まず、ＳＴＥＰ１において、記憶媒体に記憶保持された計測データから、２＊ｎ歩分の歩容の計測データを取得する。ここで、１歩分の歩容（以降、要素歩容という）の計測データは、被験者Ｐの支持脚（又は遊脚）が一方の脚から他方の脚に切替わる切替わりタイミングからその次の切替わりタイミングまでの期間の計測データである。
【０１３２】
例えば、本実施形態では、被験者Ｐの一方の脚の着地（又は離床）が、該一方の脚側の接地センサ２１により検出されるタイミングから、他方の脚の着地（又は離床）が、該他方の脚側の接地センサ２１により検出されるタイミングまでの期間が要素歩容の期間とされる。そして、データ収集装置１０は、連続した２＊ｎ個の要素歩容（以降、２ｎ歩容という）の期間の計測データを、２＊ｎ歩分の歩容の計測データとして取得する。
【０１３３】
なお、２ｎ歩容のうちの１歩目の要素歩容は、被験者Ｐが実際に歩行を開始した後の任意の歩数目の要素歩容でよい。例えば、被験者Ｐの歩行状態が定常的な状態になったとみなせる歩数目（実際の歩行における３歩目等）の歩容を、２ｎ歩容のうちの１歩目の要素歩容として計測データを取得するようにしてもよい。
【０１３４】
次いで、ＳＴＥＰ２において、データ収集装置１０は、センサユニット２毎に、２ｎ歩容の期間の各時刻ｔ（２ｎ歩容の開始時刻から所定の演算処理周期毎のサンプリング時刻）において第ｉグローバル座標系Ｃg_iで見た第ｉセンサユニット２の加速度ベクトル（瞬時値）である第ｉグローバル加速度ベクトルを算出する。
【０１３５】
以降、第ｉグローバル加速度ベクトルを表す参照符号として↑acc_i_global又は↑acc_i_global(t)を用いる。また、第ｉローカル座標系Ｃs_iで見た第ｉセンサユニット２の加速度ベクトル（瞬時値）を第ｉローカル加速度ベクトルと称し、その参照符号として↑acc_i_local又は↑acc_i_local(t)を用いる。
【０１３６】
なお、参照符号中の“↑”は、複数成分から成るベクトルを表す記号、末尾の(t)は、時刻ｔでの値を意味する記号である。
【０１３７】
第ｉグローバル加速度ベクトル↑acc_i_global(t)は、時刻ｔにおける第ｉセンサユニット２の加速度センサ５の検出出力により示される第ｉローカル加速度ベクトル↑acc_i_local(t)の検出値と、時刻ｔにおける第ｉセンサユニット２の姿勢（第ｉグローバル座標系Ｃg_iで見た姿勢）の推定値により規定される回転行列（姿勢変換行列）Ｍ_i(t)とから次式（１）により算出される。
【０１３８】
【数１】
2015186515.tif 000003
なお、上記式（１ａ），（１ｂ）の右辺の上付き添え字“Ｔ”は転置を意味する。また、式（１ｃ）により表される回転行列Ｍ_i(t)の第１列、第２列、第３列は、それぞれ、第ｉグローバル座標系で見た第ｉローカル座標系のｘ軸の基底ベクトル、ｙ軸の基底ベクトル、ｚ軸の基底ベクトルを意味する。
【０１３９】
なお、本実施形態では、第ｉグローバル加速度ベクトル↑acc_i_globalのＺ軸方向成分acc_i_z_global(t)は、以降の処理で使用しない。従って、Ｚ軸方向成分acc_i_z_global(t)を算出する処理を省略してもよい。
【０１４０】
次いで、ＳＴＥＰ３において、データ収集装置１０は、２ｎ歩容のうち、被験者Ｐの右脚が支持脚となるｎ個の要素歩容（以降、Ｒ支持歩容という）の期間のそれぞれの時間幅（各Ｒ支持歩容の開始時から終了時までの時間幅）であるＲ支持歩容周期period_R_k（ｋ＝１，２，…，ｎ）の平均値ave_period_Rと、被験者Ｐの左脚が支持脚となるｎ個の要素歩容（以降、Ｌ支持歩容という）の期間のそれぞれの時間幅であるＬ支持歩容周期period_L_k（ｋ＝１，２，…，ｎ）の平均値ave_period_Lとをそれぞれ、次式（２ａ），（２ｂ）により算出する。
【０１４１】
ave_period_R＝(period_R_1＋period_R_2＋……＋period_R_n)／ｎ  …（２ａ）
ave_period_L＝(period_L_1＋period_L_2＋……＋period_L_n)／ｎ  …（２ｂ）
次いで、ＳＴＥＰ４において、データ収集装置１０は、Ｒ支持歩容及びＬ支持歩容のそれぞれの時間軸を正規化するための正規化Ｒ支持歩容周期norm_period_Rと正規化Ｌ支持歩容周期norm_period_Lとを算出する。
【０１４２】
上記正規化Ｒ支持歩容周期norm_period_R及び正規化Ｌ支持歩容周期norm_period_Lは、それらの値の総和（＝norm_period_R＋norm_period_L）が任意の定数値（≠０）となり、且つ、norm_period_Rとnorm_period_Lとの比率が、ave_period_Rとave_period_Lとの比率に一致するように決定される。
【０１４３】
具体的には、データ収集装置１０は、次式（３ａ），（３ｂ）によりnorm_period_R及びnorm_period_Lを算出する。
【０１４４】
norm_period_R＝Ｔnorm＊ave_period_R／(ave_period_R＋ave_period_L)  …（３ａ）
norm_period_L＝Ｔnorm＊ave_period_L／(ave_period_R＋ave_period_L)  …（３ｂ）
式（３ａ），（３ｂ）のＴnorm（＝norm_period_R＋norm_period_L）は、２歩分の歩容の期間の時間幅を正規化して表す値として任意に設定した定数値Ｔnorm（≠０）であり、例えば“１”である。以降、Ｔnormを正規化基準時間幅という。
【０１４５】
なお、正規化基準時間幅Ｔnormは、“１”以外の定数値であってもよい。例えば、Ｔnorm＝ave_period_R＋ave_period_Lであってもよい。このようにした場合には、norm_period_R＝ave_period_R、norm_period_L＝ave_period_Lとなるため、ＳＴＥＰ４の処理は実質的に不要である。
【０１４６】
次いで、ＳＴＥＰ５において、データ収集装置１０は、前記２ｎ歩容を構成する各Ｒ支持歩容及び各Ｌ支持歩容のそれぞれの期間における第ｉグローバル加速度ベクトル↑acc_i_globalの各成分に波形を表す時間軸を、正規化基準時間幅Ｔnormを２歩分の歩容の期間の時間幅とする時間軸に正規化する。
【０１４７】
この正規化の処理は、２ｎ歩容のうちの各Ｒ支持歩容周期と各Ｌ支持歩容周期とをそれぞれ、norm_period_R、norm_period_Lに一致させるように、２ｎ歩容の各要素歩容毎に、第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形を、時間軸方向にスケール変換（拡大又は縮小）することでなされる。
【０１４８】
この場合、各Ｒ支持歩容における第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形の時間軸方向の拡大又は縮小の倍率は、該Ｒ支持歩容周期の実際の時間幅をＴ_Rとおくと、（norm_period_R／Ｔ_R）である。同様に、各Ｌ支持歩容における第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形の時間軸方向の拡大又は縮小の倍率は、該Ｌ支持歩容周期の実際の時間幅をＴ_Lとおくと、（norm_period_L／Ｔ_L）である。
【０１４９】
以上のＳＴＥＰ５の処理により、２ｎ歩容の各Ｒ支持歩容における第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形は、互いに同一の時間幅（＝norm_period_R）を有する波形に、時間軸方向でスケール変換される。同様に、２ｎ歩容の各Ｌ支持歩容における第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形は、互いに同一の時間幅（＝norm_period_L）を有する波形に、時間軸方向でスケール変換される。
【０１５０】
以降、上記の如く時間軸を正規化した第ｉグローバル加速度ベクトル↑acc_i_globalを、正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalと表記し、そのＸ軸方向成分、Ｙ軸方向成分、Ｚ軸方向成分をそれぞれnorm_acc_i_x_global、norm_acc_i_y_global、norm_acc_i_z_globalと表記する。また、正規化した時間軸を正規化時間軸という。
【０１５１】
第ｉグローバル加速度ベクトル↑acc_i_globalのいずれかの成分、例えばＸ軸方向成分acc_i_x_globalの波形（実際の時間軸上での波形）の例と、上記正規化処理により得られる正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalのＸ軸方向成分norm_acc_i_x_globalの波形（正規化時間軸上での波形）の例とをそれぞれ図５（ａ）、図５（ｂ）に示す。図５（ａ）は、２ｎ歩容のうちの１歩目、２歩目、３歩目、４歩目の各要素歩容におけるacc_i_x_globalの波形を例示し、図５（ｂ）は、１歩目、２歩目、３歩目、４歩目の各要素歩容におけるnorm_acc_i_x_globalの波形を例示している。
【０１５２】
なお、図示例では、奇数歩目（１歩目、３歩目）の要素歩容がＲ支持歩容、偶数歩目（２歩目、４歩目）の要素歩容がＬ支持歩容である。また、図５（ａ）では、参照符号ｔで示す横軸が時間軸を表し、図５（ｂ）では参照符号τで示す横軸が正規化時間軸を表している。
【０１５３】
図示の如く、Ｒ支持歩容及びＬ支持歩容のそれぞれについて、正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの波形は、正規化時間軸上で同一の時間幅を有する波形となる。
【０１５４】
次に、ＳＴＥＰ６において、データ収集装置１０は、２ｎ歩容のうちのｎ個のＲ支持歩容の正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの各成分の波形と、ｎ個のＬ支持歩容の正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの各成分の波形とを各々平均化してなる平均波形データを生成する。
【０１５５】
具体的には、Ｒ支持歩容に関する正規化時間軸の各時刻τ（正規化時間軸上での各Ｒ支持歩容の開始時刻から所定の刻み時間毎の各時刻）において、ｎ個のＲ支持歩容のそれぞれの正規化第ｉグローバル加速度ベクトル↑norm_acc_i_global(τ)の各成分毎の平均値を算出することで、ｎ個のＲ支持歩容の正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの各成分値の平均波形データが生成される。
【０１５６】
同様に、Ｌ支持歩容に関する正規化時間軸の各時刻τ（正規化時間軸上での各Ｌ支持歩容の開始時刻から所定の刻み時間毎の各時刻）において、ｎ個のＬ支持歩容のそれぞれの正規化第ｉグローバル加速度ベクトル↑norm_acc_i_global(τ)の各成分毎の平均値を算出することで、ｎ個のＬ支持歩容の正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの各成分の平均波形データが生成される。
【０１５７】
なお、参照符号の末尾の(τ)は、正規化時間軸上の時刻τでの値を意味する記号である。
【０１５８】
以降、Ｒ支持歩容に関する上記平均波形データにおける正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの値（瞬時値）をＲ側正規化第ｉ平均加速度ベクトル、Ｌ支持歩容に関する上記平均波形データにおける正規化第ｉグローバル加速度ベクトル↑norm_acc_i_globalの値（瞬時値）をＬ側正規化第ｉ平均加速度ベクトルと称する。そして、Ｒ側正規化第ｉ平均加速度ベクトルの参照符号として、↑norm_R_ave_acc_i_globalを使用し、Ｌ側正規化第ｉ平均加速度ベクトルの参照符号として、↑norm_L_ave_acc_i_globalを使用する。
【０１５９】
また、これらのＲ側正規化第ｉ平均加速度ベクトル及びＬ側正規化第ｉ平均加速度ベクトルを総称的に正規化第ｉ平均加速度ベクトルと称し、その参照符号として、↑norm_ave_acc_i_globalを使用する。
【０１６０】
以上のＳＴＥＰ６の処理により得られるＲ支持歩容の平均波形データにより示される波形（例えば↑norm_R_ave_acc_i_globalのＸ軸方向成分norm_R_ave_acc_i_x_globalの波形）の例と、Ｌ支持歩容の平均波形データにより示される波形（例えば↑norm_L_ave_acc_i_globalのＸ軸方向成分norm_L_ave_acc_i_x_globalの波形）の例とを、正規化時間軸上に連続的に並べた状態で図５（ｃ）に示す。
【０１６１】
本実施形態では、Ｒ側正規化第ｉ平均加速度ベクトル↑norm_R_ave_acc_i_globalの波形とＬ側正規化第ｉ平均加速度ベクトル↑norm_L_ave_acc_i_globalの波形とを各成分毎に正規化時間軸上で図５（ｃ）に例示した如く連続的に並べた波形（以降、ＲＬ正規化第ｉ平均加速度ベクトル波形という）が、被験者Ｐの２歩分の歩容の期間（正規化時間軸上でＴnormの時間幅の期間）における第ｉ加速度ベクトル↑acc_i_globalの平均波形データの波形として使用される。
【０１６２】
なお、ＲＬ正規化第ｉ平均加速度ベクトル波形におけるＲ側正規化第ｉ平均加速度ベクトル↑norm_R_ave_acc_i_globalの波形とＬ側正規化第ｉ平均加速度ベクトル↑norm_L_ave_acc_i_globalの波形との並び順はＲ側、Ｌ側の順番、あるいは、Ｌ側、Ｒ側の順番のいずれでもよい。
【０１６３】
Ｒ側正規化第ｉ平均加速度ベクトル↑norm_R_ave_acc_i_globalの波形は、被験者ＰのＲ支持歩容での平均的な第ｉ加速度ベクトル↑acc_i_globalの波形を正規化時間軸上で表現した波形に相当する。また、Ｌ側正規化第ｉ平均加速度ベクトル↑norm_L_ave_acc_i_globalの波形は、被験者ＰのＬ支持歩容での平均的な第ｉ加速度ベクトル↑acc_i_globalの波形を正規化時間軸上で表現した波形に相当する。
【０１６４】
以降の説明では、Ｒ側正規化第ｉ平均加速度ベクトル↑norm_R_ave_acc_i_globalの波形に対応するＲ支持歩容を正規化平均Ｒ支持歩容、Ｌ側正規化第ｉ平均加速度ベクトル↑norm_L_ave_acc_i_globalの波形に対応するＬ支持歩容を正規化平均Ｌ支持歩容という。また、これらの正規化平均Ｒ支持歩容と正規化平均Ｌ支持歩容とを連続させた２歩分の歩容（ＲＬ正規化第ｉ平均加速度ベクトル波形に対応する歩容）を２歩分正規化平均歩容という。
【０１６５】
補足すると、本実施形態では、以上説明したＳＴＥＰ２〜６の処理が、本発明におけるセンサユニット加速度推定手段に相当するセンサユニット加速度推定部１１により実行される処理である。この場合、各センサユニット２についての第ｉグローバル加速度ベクトル↑acc_i_globalの各成分が本発明における対象状態量に相当する。そして、ＳＴＥＰ２の処理が、本発明における基本推定値生成処理に相当し、ＳＴＥＰ３〜６の処理が、平均波形データ生成処理に相当する。
【０１６６】
次に、ＳＴＥＰ７において、データ収集装置１０は、上記ＲＬ正規化第ｉ平均加速度ベクトル波形のデータ（平均波形データ）を用いて、被験者Ｐの上体の角加速度（２歩分正規化平均歩容の期間の各時刻τでの瞬時値）を算出する。なお、このＳＴＥＰ７で算出する角加速度は、本実施形態では、被験者Ｐの上体のロール方向（Ｘ軸周り方向）の角加速度θx_dot2_with_offset、及び被験者Ｐの上体のピッチ方向（Ｙ軸周り方向）の角加速度θy_dot2_with_offsetである。そして、これらの角加速度θx_dot2_with_offset，θy_dot2_with_offsetは、後述するオフセット成分を除去する前の角加速度である。
【０１６７】
ここで、各センサユニット２は、被験者Ｐの上体に比して十分に小型なものであり、各センサユニット２の被験者Ｐの上体への装着面（接触面）の面積は比較的小さい。このため、各センサユニット２の瞬時的な姿勢は、被験者Ｐの上体の実際の瞬時的な姿勢に対してずれを生じ易い。ひいては、各センサユニット２自体の角加速度の瞬時値は、被験者Ｐの上体の実際の角加速度の瞬時値に対してずれを生じやすい。
【０１６８】
一方、各センサユニット２の並進加速度は、被験者Ｐの上体における該センサユニット２の装着箇所の並進加速度に対して比較的高い追従性を有する。このため、センサユニット２，２のそれぞれのＹ軸方向の並進加速度の差の変化は、被験者Ｐの上体のロール方向（Ｘ軸周り方向）の角加速度の変化に対して比較的高い追従性を有する。同様に、センサユニット２，２のそれぞれＸ軸方向の並進加速度の差の変化は、被験者Ｐの上体のピッチ方向（Ｙ軸周り方向）の角加速度の変化に対して高い追従性を有する。
【０１６９】
そこで、ＳＴＥＰ７では、データ収集装置１０は、センサユニット２のうちの第１センサユニット２（上側のセンサユニット２）についての正規化第１平均加速度ベクトル↑norm_ave_acc_1_globalのＹ軸方向成分norm_ave_acc_1_y_globalと、第２センサユニット２（下側のセンサユニット２）についての正規化第２平均加速度ベクトル↑norm_ave_acc_2_globalのＹ軸方向成分norm_ave_acc_2_y_globalとから、次式（４ａ）に示す線形結合によって被験者Ｐの上体のロール方向（Ｘ軸方向）の角加速度θx_dot2_with_offsetを算出する。
【０１７０】
また、データ収集装置１０は、第１センサユニット２についての正規化第１平均加速度ベクトル↑norm_ave_acc_1_globalのＸ軸方向成分norm_ave_acc_1_x_globalと、第２センサユニット２についての正規化第２平均加速度ベクトル↑norm_ave_acc_2_globalのＸ軸方向成分norm_ave_acc_2_x_globalとから、次式（４ｂ）に示す線形結合によって、被験者Ｐの上体のピッチ方向（Ｙ軸方向）の角加速度θy_dot2_with_offsetを算出する。
【０１７１】
θx_dot2_with_offset(τ)
＝−(norm_ave_acc_1_y_global(τ)−norm_ave_acc_2_y_global(τ))／(h1−h2)
…（４ａ）
θy_dot2_with_offset
＝(norm_ave_acc_1_x_global(τ)−norm_ave_acc_2_x_global(τ))／(h1−h2)
…（４ｂ）
ただし、
ｈ１：床面からの第１センサユニット２の高さ
ｈ２：床面からの第２センサユニット２の高さ（＜ｈ１）
上式（４ａ），（４ｂ）における右辺の各項の係数を規定するｈ１−ｈ２は、第１センサユニット２と第２センサユニット２との上下方向の距離に相当する。
【０１７２】
なお、ｈ１，ｈ２の値としては、被験者Ｐの直立姿勢状態での第１センサユニット２及び第２センサユニット２のそれぞれ高さの、あらかじめ計測された値（一定値）が用いられる。
【０１７３】
上記の如くＳＴＥＰ７で算出される角加速度の値は、一般に、各センサユニット２の角速度センサ４又は加速度センサ５の検出出力のドリフト等に起因して、オフセット成分を含む。
【０１７４】
そこで、データ収集装置１０は、次に、ＳＴＥＰ８において、上記オフセット成分を除去する処理を実行する。
【０１７５】
ここで、一般に、人の定常的な直進歩行の２歩分の歩容の期間における上体の実際の角加速度の平均値は、ゼロもしくはほぼゼロとなる。従って、前記２歩分正規化平均歩容におけるロール方向の角加速度θx_dot2_with_offset(τ)の平均値θx_dot2_aveと、ピッチ方向の角加速度θy_dot2_with_offset(τ)の平均値θy_dot2_aveとがそれぞれ、θx_dot2_with_offset(τ)のオフセット成分、θy_dot2_with_offset(τ)のオフセット成分に相当すると見なせる。
【０１７６】
そこで、ＳＴＥＰ８では、データ収集装置１０は、次式（５ａ），（５ｂ）により、前記２歩分正規化平均歩容における被験者Ｐの上体のロール方向の角加速度θx_dot2_with_offsetの平均値θx_dot2_aveとピッチ方向の角加速度θy_dot2_with_offsetの平均値θy_dot2_aveとを算出する。
【０１７７】
【数２】
2015186515.tif 000004
そして、データ収集装置１０は、次式（６ａ），（６ｂ）で示す如く、２歩分正規化平均歩容における各時刻τにおけるロール方向の角加速度θx_dot2_with_offset(τ)から上記の如く算出した平均値θx_dot2_ave（オフセット成分）を差し引いた値と、ピッチ方向の角加速度θy_dot2_with_offset(τ)から上記の如く算出した平均値θｙ_dot2_ave（オフセット成分）を差し引いた値とをそれぞれ、２歩分正規化平均歩容における被験者Ｐの上体の実際のロール方向の角加速度θx_dot2(τ)の推定値（オフセット成分の除去後の推定値）、ピッチ方向の角加速度θy_dot2(τ)の推定値（オフセット成分の除去後の推定値）として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体の角加速度θx_dot2(τ)，θy_dot2(τ)の推定値の波形データが得られることとなる。
【０１７８】
θx_dot2(τ)＝θx_dot2_with_offset(τ)−θx_dot2_ave    …（６ａ）
θy_dot2(τ)＝θy_dot2_with_offset(τ)−θy_dot2_ave    …（６ａ）
次に、ＳＴＥＰ９において、データ収集装置１０は、２歩分正規化平均歩容の各時刻τにおける被験者Ｐの上体の角速度（ロール方向の角速度θx_dot(τ)及びピッチ方向の角速度θy_dot(τ)）の推定値を算出する。
【０１７９】
この場合、θx_dot(τ)、θy_dot(τ)の推定値は次のように算出される。すなわち、まず、次式（７ａ），（７ｂ）で示す如く、ロール方向の角加速度θx_dot2(τ)の推定値、ピッチ方向の角加速度θy_dot2(τ)の推定値をそれぞれ積分することにより、２歩分正規化平均歩容の初期時刻（τ＝０）での角速度の値（初期値θx_dot(0)，θy_dot(0)）をゼロと仮定した場合における各時刻τでのロール方向の角速度の値θx_dot_with_offset(τ)及びピッチ方向の角速度の値θy_dot_with_offset(τ)が算出される。
【０１８０】
【数３】
2015186515.tif 000005
ここで、一般に、人の定常的な直進歩行の２歩分の歩容の期間における上体の角速度の平均値は、一般に、ゼロもしくはほぼゼロとなる。そこで、本実施形態では、データ収集装置１０は、２歩分正規化平均歩容におけるロール方向の角速度θx_dot(τ)及びピッチ方向の角速度θy_dot(τ)のそれぞれの平均値がゼロになるように、ロール方向の角速度θx_dot(τ)及びピッチ方向の角速度θy_dot(τ)のそれぞれの初期値θx_dot(0)，θy_dot(0)を決定する。
【０１８１】
具体的には、θx_dot(0)，θy_dot(0)は、それぞれ、次式（８ａ），（８ｂ）に算出される。
【０１８２】
【数４】
2015186515.tif 000006
すなわち、２歩分正規化平均歩容におけるθx_dot_with_offset(τ)の平均値θx_dot_aveの逆極性の値と、θy_dot_with_offset(τ)の平均値θy_dot_aveの逆極性の値とがそれぞれθx_dot(0)，θy_dot(0)として算出される。
【０１８３】
そして、データ収集装置１０は、次式（９ａ），（９ｂ）で示す如く、上記初期値θx_dot(0)，θy_dot(0)をそれぞれθx_dot_with_offset(τ)、θy_dot_with_offset(τ)に加算してなる値をそれぞれ、各時刻τでのロール方向の角速度θx_dot(τ)の推定値、ピッチ方向の角速度θy_dot(τ)の推定値として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体の角速度θx_dot(τ)，θy_dot(τ)の推定値の波形データが得られることとなる。
【０１８４】
θx_dot(τ)＝θx_dot_with_offset(τ)＋θx_dot(0)    …（９ａ）
θy_dot(τ)＝θy_dot_with_offset(τ)＋θy_dot(0)    …（９ｂ）
次に、ＳＴＥＰ１０において、データ収集装置１０は、２歩分正規化平均歩容の各時刻τにおける被験者Ｐの上体の姿勢（ロール方向の傾斜角θx(τ)及びピッチ方向の傾斜角θy(τ)）の推定値を算出する。
【０１８５】
この場合、θx(τ)、θy(τ)の推定値は次のように算出される。すなわち、まず、次式（１０ａ），（１０ｂ）で示す如く、ロール方向の角速度θx_dot(τ)の推定値、ピッチ方向の角速度θy_dot(τ)の推定値をそれぞれ積分することにより、２歩分正規化平均歩容の初期時刻（τ＝０）での傾斜角の値（初期値θx(0)，θy(0)）をゼロと仮定した場合における各時刻τでのロール方向の傾斜角の値θx_with_offset(τ)及びピッチ方向の傾斜角の値θy_with_offset(τ)が算出される。
【０１８６】
【数５】
2015186515.tif 000007
ここで、２歩分正規化平均歩容における被験者Ｐの上体のロール方向の傾斜角θx(τ)の平均値は、前記２ｎ歩容における第１センサユニット２（首の付け根付近のセンサユニット２）の姿勢のうちのロール方向の実際の傾斜角の平均値（以降、これをθx_neck_aveと表記する）と、被験者Ｐの直立姿勢における第１センサユニット２の姿勢のうちのロール方向の実際の傾斜角（以降、これをθx_neck_up_rightと表記する）との差（以降、これをθx_neck_offsetと表記する）に一致もしくはほぼ一致すると考えられる。
【０１８７】
同様に、２歩分正規化平均歩容における被験者Ｐの上体のピッチ方向の傾斜角θy(τ)の平均値は、前記２ｎ歩容における第１センサユニット２（首の付け根付近のセンサユニット２）の姿勢のうちのピッチ方向の実際の傾斜角の平均値（以降、これをθy_neck_aveと表記する）と、被験者Ｐの直立姿勢における第１センサユニット２のピッチ方向の実際の傾斜角の推定値（以降、これをθy_neck_up_rightと表記する）との差（以降、これをθｙ_neck_offsetと表記する）に一致もしくはほぼ一致すると考えられる。
【０１８８】
そこで、本実施形態では、データ収集装置１０は、２歩分正規化平均歩容における被験者Ｐの上体のロール方向の傾斜角θx_(τ)及びピッチ方向の傾斜角θy_dot(τ)のそれぞれの平均値が、上記θx_neck_offsetの推定値と上記θy_neck_offsetの推定値とに各々一致するように、ロール方向の傾斜角θx(τ)及びピッチ方向の傾斜角θy(τ)のそれぞれの初期値θx(0)，θy(0)を決定する。
【０１８９】
具体的には、θx_neck_offsetの推定値とθy_neck_offsetの推定値とが、それぞれ次式（１１ａ），（１１ｂ）により算出される。
【０１９０】
θx_neck_offset＝θx_neck_ave−θx_neck_up_right    …（１１ａ）
θy_neck_offset＝θy_neck_ave−θy_neck_up_right    …（１１ｂ）
ここで、θx_neck_ave，θy_neck_aveは、第１センサユニット２の演算処理部７で前記した如く逐次推定された第１センサユニット２の姿勢におけるロール方向の傾斜角の値と、ピッチ方向の傾斜角の値とからそれぞれ算出される。また、θx_neck_up_right，θy_neck_up_rightは、被験者Ｐの直立姿勢であらかじめ計測されてデータ収集装置１０に記憶保持された値である。
【０１９１】
また、２歩分正規化平均歩容について前記式（１０ａ），（１０ｂ）によりそれぞれ算出されたθx_with_offset(τ)，θy_with_offset(τ)のそれぞれの平均値θx_ave，θy_aveがそれぞれ次式（１２ａ），（１２ｂ）により算出される。
【０１９２】
【数６】
2015186515.tif 000008
さらに、次式（１３ａ）で示す如く、θx_aveの上記の算出値とθx_neck_offsetの上記の算出値との差の逆極性の値がθx(0)として算出される。また、次式（１３ｂ）で示す如く、θy_aveの上記の算出値とθy_neck_offsetの上記の算出値との差の逆極性の値がθy(0)として算出される。
【０１９３】
θx(0)＝−(θx_ave−θx_neck_offset)    …（１３ａ）
θy(0)＝−(θy_ave−θy_neck_offset)    …（１３ｂ）
そして、データ収集装置１０は、次式（１４ａ），（１４ｂ）で示す如く、上記初期値θx(0)，θy(0)をそれぞれθx_with_offset(τ)、θy_with_offset(τ)に加算してなる値を、それぞれ各時刻τでのロール方向の傾斜角θx(τ)の推定値、ピッチ方向の傾斜角θy(τ)の推定値として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体の傾斜角θx(τ)，θy(τ)の推定値の波形データが得られることとなる。
【０１９４】
θx(τ)＝θx_with_offset(τ)＋θx(0)    …（１４ａ）
θy(τ)＝θy_with_offset(τ)＋θy(0)    …（１４ｂ）
補足すると、本実施形態では、以上説明したＳＴＥＰ７〜１０の処理が、本発明における上体角加速度推定手段及び上体角速度推定手段としての機能を有する上体姿勢運動推定部１２により実行される処理である。この場合、ＳＴＥＰ７，８の処理が、上体角加速度推定手段として処理、ＳＴＥＰ９の処理が、上体角速度推定手段としての処理に相当する。
【０１９５】
本実施形態では、以上の如く、被験者Ｐの上体のロール方向の角加速度θx_dot2(τ)、角速度θx_do2(τ)、傾斜角θx(τ)の推定値の波形データと、ピッチ方向の角加速度θy_dot2(τ)、角速度θy_do2(τ)、傾斜角θy(τ)の推定値の波形データとが該上体の姿勢の実際の動きを代表的に示すデータとして得られる。
【０１９６】
そして、データ収集装置１０は、そのディスプレイ１０ａに上記波形データ、あるいは、正規化Ｒ支持歩容周期norm_period_R及び正規化Ｌ支持歩容周期norm_period_Lの比率（Ｒ支持歩容周期period_RとＬ支持歩容周期period_Lとの平均的な比率）等を視覚的に表示する。なお、波形データ等の表示は、印刷により行うようにしてもよい。
【０１９７】
ディスプレイ１０ａに表示される波形データ等の一例を図６に示す。図示例では、上体の角加速度の波形のグラフと、正規化Ｒ支持歩容周期norm_period_Rと正規化Ｌ支持歩容周期norm_period_Lとの大小を視覚的に示すグラフとが表示されている。
【０１９８】
波形のグラフは、角加速度以外の波形（角速度、傾斜角等の波形）に切替えることも可能である。
【０１９９】
なお、データ収集装置１０から出力する波形データは、上体の角加速度、角速度又は姿勢（傾斜角）の波形データに限らず、いずれかのセンサユニット２（例えば腰部付近で被験者Ｐの上体に装着されたセンサユニット２）のグローバル座標系で見た加速度の波形データを含んでいてもよい。
【０２００】
本実施形態の上体運動計測システム１によれば、以上の如く被験者Ｐの歩行動作における上体の角加速度等の波形データが得られる。そして、その波形データを利用して、被験者Ｐの歩行状態の解析又は評価を行うことができる。
【０２０１】
例えば、健常な被験者Ｐの歩行動作では、一般に、Ｒ支持歩容における上体の動きとＬ支持歩容における上体の動きとがバランスよく行われる。
【０２０２】
そして、本実施形態では、正規化Ｒ支持歩容周期norm_period_Rと正規化Ｌ支持歩容周期norm_period_Lとの比率（又は大小関係）、あるいは、正規化Ｒ支持歩容周期norm_period_R及び正規化Ｌ支持歩容周期norm_period_Lのそれぞれの期間での上体の角加速度又は角速度又は傾斜角度の波形から、被験者Ｐの歩行動作におけるＲ支持歩容とＬ支持歩容とのバランス状態を評価することができる。
【０２０３】
ひいては、被験者Ｐの片脚に不具合があるか否か、あるいは、不具合を生じた片脚の回復状況等を評価することができる。
【０２０４】
また、健常な被験者Ｐの歩行動作では、一般に、図７（ａ）に例示する如く、各脚の着地直前において、上体が後ろに傾く方向の角加速度（Ｙ軸周りで負方向の角加速度）が発生するように上体の姿勢が自然に（無意識的に）ピッチ方向に動く。この結果、当該上体の姿勢の動きが、着地衝撃に起因する上体のピッチ方向の角加速度（Ｙ軸周りで正方向の角加速度）とバランスする。ひいては、次の歩行サイクルへの移行が滑らかに行われる。
【０２０５】
一方、脚に障害を持つ被験者Ｐ等の歩行動作では、図７（ｂ）に例示する如く、着地衝撃に対する対応力を十分に発揮できず、各脚の着地直後において、上体のピッチ方向の角加速度がＹ軸周りの正方向にスパイク形状の波形となる場合がある。
【０２０６】
従って、被験者Ｐの各脚の着地近辺における上体のピッチ方向の角加速度の波形を観測することで、被験者Ｐの不具合状態を評価できる。
【０２０７】
［第２実施形態］
次に、本発明の第２実施形態を図８等を参照して説明する。なお、本実施形態の上体運動計測システムは、前記第１実施形態の上体運動計測システム１とセンサユニットの配置構成と、データ収集装置の一部の処理だけが第１実施形態と相違する。そこで、本実施形態の説明は、第１実施形態と相違する部分を中心に行い、第１実施形態と同一の事項については説明を省略する。
【０２０８】
図８（ａ）を参照して、本実施形態の測定システムは、４個のセンサユニット（慣性センサユニット）２，２，２，２を備える。各センサユニット２の構成は、第１実施形態と同じである。
【０２０９】
４個のセンサユニット２，２，２，２のうち、２つのセンサユニット２，２は、第１実施形態と同様に、被験者Ｐの上体の背中の脊椎沿いの線上（体幹軸沿いの線上）で、被験者Ｐの首の付け根の付近の箇所と腰付近の箇所とに装着されている。以降、第１実施形態と同様に、この２つのセンサユニット２，２のうちの上側の（首の付け根付近の）センサユニット２を第１センサユニット２、下側の（腰付近の）センサユニット２を第２センサユニット２ということがある。
【０２１０】
また、残りの２つのセンサユニット２，２は、人の上体の脊椎（体幹軸）の左側の箇所と右側の箇所とで該上体の背中（背面）に装着されている。それらの装着箇所は、人の上体の脊椎に対して左右対称になるようにして、上体の同じ高さ位置に配置されている。その高さ位置としては、例えば図８（ａ）に示すように、被験者Ｐの胸付近の高さ位置、あるいは、図８（ｂ）に示すように、被験者Ｐの肩と同程度の高さ位置等を採用できる。
【０２１１】
以降、上記の如く被験者Ｐの上体の左側の箇所と右側の箇所とに配置される２つのセンサユニット２，２の一方、例えば左側のセンサユニット２を第３センサユニット２、他方（右側）のセンサユニット２を第４センサユニット２ということがある。そして、第１〜第４センサユニット２の任意の１つを第ｉセンサユニット２（ｉ＝１，２，３，４）ということがある。
【０２１２】
本実施形態の上体運動計測システムの構成は、以上説明した事項以外は、第１実施形態と同じである。
【０２１３】
次に、本実施形態の上体運動計測システムよる計測処理を説明する。
【０２１４】
第１実施形態と同様に、データ収集装置１０に各センサユニット２から送信される計測データ（被験者Ｐの歩行中の各時刻における各センサユニット２の姿勢の推定値及び加速度ベクトルの検出値を示す計測データ）が、該データ収集装置１０の記憶媒体に時系列的に記憶保持される。
【０２１５】
そして、データ収集装置１０は、その計測データを用いて、前記被験者Ｐの上体の動きを解析するための処理を、第１実施形態と同様に、図４のフローチャートに示す如く実行する。
【０２１６】
ただし、本実施形態では、ＳＴＥＰ２，４の処理は、４個のセンサユニット２のそれぞれに関して実行される。
【０２１７】
また、本実施形態では、ＳＴＥＰ７からの処理において、第１実施形態で説明した処理に加えて、被験者Ｐの上体のヨー方向（Ｚ軸方向）の角加速度、角速度、及び回転角度を推定するための処理が追加される。
【０２１８】
具体的には、データ収集装置１０は、ＳＴＥＰ７において、被験者Ｐの上体のロール方向（Ｘ軸方向）の角加速度θx_dot2_with_offset、及びピッチ方向（Ｙ軸方向）の角加速度θy_dot2_with_offsetに加えて、ヨー方向（Ｚ軸方向）の角加速度θz_dot2_with_offsetを算出する。
【０２１９】
この場合、データ収集装置１０は、第３センサユニット２（右側のセンサユニット２）についての正規化第１平均加速度ベクトル↑norm_ave_acc_3_globalのＸ軸方向成分norm_ave_acc_1_x_globalと、第４センサユニット２（左側のセンサユニット２）についての正規化第２平均加速度ベクトル↑norm_ave_acc_4_globalのＸ軸方向成分norm_ave_acc_4_x_globalとから、次式（４ｃ）に示す線形結合によって被験者Ｐの上体のヨー方向（Ｚ軸方向）の角加速度θz_dot2_with_offsetを算出する。
【０２２０】
θz_dot2_with_offset(τ)
＝(norm_ave_acc_3_x_global(τ)−norm_ave_acc_4_x_global(τ))／Δｙ…（４ｃ）
ただし、
Δｙ：第３センサユニット２と第４センサユニット２間のＹ軸方向（左右方向）の距離
上式（４ｃ）のΔｙの値としては、被験者Ｐの直立姿勢状態での第３センサユニット２と第４センサユニット２との間のＹ軸方向の距離の、あらかじめ計測された値（一定値）が用いられる。
【０２２１】
また、次のＳＴＥＰ８では、データ収集装置１０は、ロール方向及びピッチ方向の角加速度θx_dot2_with_offset，θy_dot2_with_offsetの算出値からオフセット成分を除去することに加えて、ヨー方向の角加速度θz_dot2_with_offsetの算出値からオフセット成分を除去する処理も実行する。
【０２２２】
この処理では、データ収集装置１０は、次式（５ｃ）により、前記２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の角加速度θz_dot2_with_offsetの平均値θz_dot2_aveを算出する。
【０２２３】
【数７】
2015186515.tif 000009
そして、データ収集装置１０は、次式（６ｃ）で示す如く、２歩分正規化平均歩容における各時刻τにおけるヨー方向の角加速度θz_dot2_with_offset(τ)から上記の如く算出した平均値θz_dot2_ave（オフセット成分）を差し引いた値を、被験者Ｐの上体の実際のヨー方向の角加速度θz_dot2(τ)の推定値として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の角加速度θz_do2(τ)の推定値の波形データがさらに得られることとなる。
【０２２４】
θz_dot2(τ)＝θz_dot2_with_offset(τ)−θz_dot2_ave    …（６ｃ）
また、次のＳＴＥＰ９において、データ収集装置１０は、２歩分正規化平均歩容の各時刻τにおける被験者Ｐの上体のロール方向の角速度θx_dot(τ)及びピッチ方向の角速度θy_dot(τ)の推定値に加えて、ヨー方向の角速度θz_dot(τ)の推定値を算出する。
【０２２５】
この場合、θz_dot(τ)の推定値は次のように算出される。すなわち、まず、次式（７ｃ）で示す如く、ヨー方向の角加速度θz_dot2(τ)の推定値を積分することにより、２歩分正規化平均歩容の初期時刻（τ＝０）での角速度の値（初期値θz_dot(0)）をゼロと仮定した場合における各時刻τでのヨー方向の角速度の値θz_dot_with_offset(τ)が算出される。
【０２２６】
【数８】
2015186515.tif 000010
ここで、前記したように、人の定常的な直進歩行の２歩分の歩容の期間における上体の角速度の平均値は、一般に、ゼロもしくはほぼゼロとなる。そこで、データ収集装置１０は、２歩分正規化平均歩容におけるヨー方向の角速度θz_dot(τ)の平均値がゼロになるように、ヨー方向の角速度θz_dot(τ)の初期値θz_dot(0)を決定する。
【０２２７】
具体的には、θz_dot(0)は次式（８ｃ）に算出される。
【０２２８】
【数９】
2015186515.tif 000011
すなわち、２歩分正規化平均歩容におけるθz_dot_with_offset(τ)の平均値θz_dot_aveの逆極性の値がθz_dot(0)として算出される。
【０２２９】
そして、データ収集装置１０は、次式（９ｃ）で示す如く、上記初期値θz_dot(0)をθz_dot_with_offset(τ)に加算してなる値を、各時刻τでのヨー方向の角速度θz_dot(τ)の推定値として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の角速度θz_dot(τ)の推定値の波形データがさらに得られることとなる。
【０２３０】
θz_dot(τ)＝θz_dot_with_offset(τ)＋θz_dot(0)    …（９ｃ）
次のＳＴＥＰ１０において、データ収集装置１０は、２歩分正規化平均歩容の各時刻τにおける被験者Ｐの上体のロール方向の傾斜角θx(τ)及びピッチ方向の傾斜角θy(τ)の推定値に加えて、ヨー方向の回転角θz(τ)の推定値を算出する。
【０２３１】
この場合、θz(τ)の推定値は次のように算出される。すなわち、まず、次式（１０ｃ）で示す如く、ヨー方向の角速度θz_dot(τ)の推定値を積分することにより、２歩分正規化平均歩容の初期時刻（τ＝０）での回転角の値（初期値θz(0)）をゼロと仮定した場合における各時刻τでのヨー方向の回転角の値θz_with_offset(τ)が算出される。
【０２３２】
【数１０】
2015186515.tif 000012
ここで、２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の傾斜角θx(τ)の平均値は、一般にゼロもしくはほぼゼロになる。
【０２３３】
そこで、本実施形態では、データ収集装置１０は、２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の回転角θz_(τ)の平均値がゼロになるように、ヨー方向の回転角θz(τ)の初期値θz(0)を決定する。
【０２３４】
具体的には、次式（１３ｃ）で示す如く、２歩分正規化平均歩容について前記式（１０ｃ）により算出されたθz_with_offset(τ)の平均値θz_aveの逆極性の値がθx(0)として算出される。
【０２３５】
【数１１】
2015186515.tif 000013
そして、データ収集装置１０は、次式（１４ｃ）で示す如く、上記初期値θz(0)をθz_with_offset(τ)に加算してなる値を、各時刻τでのヨー方向の回転角θz(τ)の推定値として算出する。これにより、２歩分正規化平均歩容における被験者Ｐの上体のヨー方向の回転角θz(τ)の推定値の波形データがさらに得られることとなる。
【０２３６】
θz(τ)＝θz_with_offset(τ)＋θz(0)    …（１４ｃ）
本実施形態は、以上説明した事項以外は、第１実施形態と同じである。
【０２３７】
本実施形態によれば、被験者Ｐの歩行状態の解析あるいは評価を行うために、上体のロール方向及びピッチ方向の角加速度、角速度及び傾斜角度の他、さらに、ヨー方向の角加速度、角速度及び回転角度の推定値及びその波形データを活用することができる。
【０２３８】
［変形態様］
次に、以上説明した各実施形態に関連する変形態様をいくつか説明する。
【０２３９】
前記各実施形態では、前記式（４ａ），（４ｂ），（４ｃ）により、被験者Ｐの上体の角加速度の推定値を算出するようにした。ただし、例えば上体が直立姿勢状態からロール方向に傾いた状態では、上体のロール方向の角加速度の変化に応じて、各センサユニット２（特に、上体の上部側のセンサユニット２）の並進加速度ベクトルは、Ｙ軸方向成分に加えて、Ｚ軸方向成分も変化する。
【０２４０】
同様に、例えば上体が直立姿勢状態からピッチ方向に傾いた状態では、上体のピッチ方向の角加速度の変化に応じて、各センサユニット２（特に、上体の上部側のセンサユニット２）の並進加速度ベクトルは、Ｘ軸方向成分に加えて、Ｚ軸方向成分も変化する。
【０２４１】
これらのことを考慮すると、ロール方向、ピッチ方向、及びヨー方向の各方向の上体の角加速度の推定値は、より一般的には、上体に装着する複数のセンサユニット２の加速度を線形結合することによって算出するようにしてもよい。
【０２４２】
例えば、前記ＳＴＥＰ７の処理において、次式（１５ａ），（１５ｂ），（１５ｃ）により、ロール方向、ピッチ方向、及びヨー方向の各方向の上体の角加速度の推定値（オフセット成分を除去する前の推定値）を算出するようにしてもよい。
【０２４３】
【数１２】
2015186515.tif 000014
なお、Ｎは、上体に装着するセンサユニット２の個数、ｗ_xy_i，ｗxz_i，ｗ_yx_i，ｗyz_i，ｗ_zx_i、ｗ_zy_i（ｉ＝１，２，…，Ｎ）は重み係数である。これらの重み係数ｗ_xy_i，ｗxz_i，ｗ_yx_i，ｗyz_i，ｗ_zx_i、ｗ_zy_iの値は、あらかじめ実験等に基づいて設定すればよい。この場合、重み係数ｗ_xy_i，ｗxz_i，ｗ_yx_i，ｗyz_i，ｗ_zx_i、ｗ_zy_iのうちの一部の重み係数はゼロであってもよい。
【０２４４】
このように上体の角加速度の推定値を算出することで、該角加速度の推定値の信頼性をより高めることができる。
【０２４５】
また、前記各実施形態では、上体に装着する複数のセンサユニット２の配置構成として、図１（ｂ）あるいは、図８（ａ）もしくは図８（ｂ）に示した配置構成を例示したが、例えば、図９に示す配置構成を採用してもよい。この例では、３個のセンサユニット２のうちの一つが、被験者Ｐの腰付近の高さで上体に装着され、２個のセンサユニット２が、上体の上部（例えば肩付近の高さの部分）で、体幹軸Ｌに対して左右対称となるように上体に装着される。
【０２４６】
また、前記各実施形態では、各センサユニット２を上体の背中（背面）に装着したが、各センサユニット２を上体の前面側（腹側の面）に装着してもよい。
【０２４７】
また、前記各実施形態では、正規化第ｉ平均加速度ベクトル↑norm_ave_acc_i_globalの各成分の波形を算出した上で、２歩分正規化平均歩容における上体の角加速度の推定値を算出するようにした。ただし、正規化第ｉ平均加速度ベクトル↑norm_ave_acc_i_globalの各成分の代わりに、前記ＳＴＥＰ２で算出される第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の値をそのまま使用して、式（４ａ），（４ｂ），（４ｃ）又は式（１５ａ），（１５ｂ），（１５ｃ）と同様の演算を行うことで、２ｎ歩容の全体における上体の各方向の角加速度の推定値の時系列を算出するようにしてもよい。
【０２４８】
そして、この場合には、上記の如く算出した角加速度の推定値の波形を、第ｉグローバル加速度ベクトル↑acc_i_globalの各成分の波形の正規化及び平均化と同様の処理によって、正規化及び平均化することで、各方向の角加速度の平均波形データを生成するようにしてもよい。
【０２４９】
この場合、上体の角加速度の平均波形データを生成した後は、前記ＳＴＥＰ８〜１０と同様の処理を行うことで、オフセット成分を除去した角加速度の推定値の算出、角速度の算出、上体の姿勢（傾斜角あるいは回転角）の推定値を算出できる。
【０２５０】
また、前記各実施形態では、ＳＴＥＰ６で生成した正規化第ｉ平均加速度ベクトル↑norm_ave_acc_i_globalの各成分の波形データ（平均波形データ）に対して、前記ＳＴＥＰ８と同様の処理によって、オフセット成分を除去する処理を実行するようにしてもよい。
【０２５１】
すなわち、２歩分正規化平均歩容の期間における↑norm_ave_acc_i_globalの各成分の平均値をオフセット成分として求め、このオフセット成分を当該成分の値から差し引くことで、該オフセット成分を除去してなる平均波形データを求めるようにしてもよい。このようにした場合には、前記ＳＴＥＰ８の処理を省略してもよい。
【０２５２】
また、前記各実施形態では、Ｒ支持歩容及びＬ支持歩容の一方から他方への切替わりのタイミンングをデータ収集装置１０で認識するために、接地センサ２１，２１を用いた。ただし、被験者の上体に装着するいずれかのセンサユニット２の加速度センサ５の検出出力により示される所定方向の加速度の検出値（センサユニットローカル座標系Ｃｓで見た該センサユニット２の加速度）又は、いずれかのセンサユニット２についての第ｉグローバル加速度ベクトルの所定方向成分の変化に基づいて、上記切替わりのタイミングを認識することもできる。
【０２５３】
例えば、被験者Ｐの腰付近に装着されるセンサユニット２（第２センサユニット２）についての第２グローバル加速度ベクトル↑acc_2_globalのＸ軸方向成分acc_2__x_globalの推定値は、被験者Ｐの歩行動作において、図１０に例示する如く、各要素歩容における遊脚の着地タイミングで、減速方向（Ｘ軸の負方向）の加速度に急激に変化する。このことは、第２センサユニット２についての第２ローカル加速度ベクトル↑acc_i_localの検出値のＸ軸方向成分acc_2__x_localについても同様である。
【０２５４】
従って、↑acc_2_globalのＸ軸方向成分acc_2__x_globalの微分値（時間的変化率）、又は↑acc_2__x_localのＸ軸方向成分acc_2_x_localの微分値（時間的変化率）を逐次算出し、この微分値が負の値（第２センサユニット２がＸ軸方向で減速する方向の値）で、且つ、その絶対値又は二乗値が、所定の設定閾値以上になったタイミングを、Ｒ支持歩容及びＬ支持歩容の一方から他方への切替わりのタイミンングとして認識するようにしてもよい。
【０２５５】
さらに、このようにＲ支持歩容及びＬ支持歩容の一方から他方への切替わりのタイミンングとして認識する場合、例えば第２センサユニット２についての第２グローバル加速度ベクトル↑acc_2_globalのＹ軸方向成分acc_2_y_globalの変化パターンに基づいて、各要素歩容において支持脚となっている脚が右脚及び左脚のいずれであるかを認識することも可能である。
【０２５６】
すなわち、図１１に例示する如く、↑acc_2_globalのＹ軸方向成分acc_2_y_globalの波形は、Ｒ支持歩容とＬ支持歩容とで各別の特徴的な変化パータンを有するので、該変化パターンの特徴から、各要素歩容において支持脚となっている脚が右脚及び左脚のいずれであるかを認識することが可能である。
【０２５７】
この場合、acc_2_y_globalの波形の変化パターンの認識は、例えば、次のような手法で行うことができる。すなわち、図１１に示す如く、各要素歩容の期間を複数の期間に分割する。
【０２５８】
そして、その各分割期間ごとに、該分割期間の始端から終端までのacc_2_y_globalの変化量（終端における値と始端における値との差）と、該分割期間でのacc_2_y_globalの極性（詳しくは該分割期間内で正極性となっている時間と負極性となっている時間とのうちのより長い方の時間に対応する極性）とを指標として、各分割期間でのacc_2_y_globalの変化パターンを分類する。
【０２５９】
そして、各要素歩容での全ての分割期間での変化パターンと、Ｒ支持歩容とＬ支持歩容とに対応してあらかじめ設定した変化パターンとの一致度合を照合することで、各要素歩容で支持脚となっている脚が右脚及び左脚のいずれであるかを認識する。
【０２６０】
また、前記接地センサ２１の代わりに、例えば被験者Ｐの各足平部に加速度センサを装着してもよい。この場合、各足平部に装着した加速度センサにより検出される加速度の絶対値が、各脚の着地又は離床のタイミングで急激に大きなものとなる。従って、上記加速度の絶対値の変化に基づいて検出される各脚の着地又は離床のタイミングを、、Ｒ支持歩容及びＬ支持歩容の一方から他方への切替わりのタイミングとして認識することも可能である。
【符号の説明】
【０２６１】
１…上体運動計測システム、２…慣性センサユニット、４…角速度センサ、５…加速度センサ、７…演算処理部（センサユニット姿勢推定手段）、１１…センサユニット加速度推定部（センサユニット加速度推定手段）、１２…上体姿勢運動推定部（上体角加速度推定手段、上体角速度推定手段）。
【図６】
2015186515.tif 000015
【図１】
2015186515.tif 000016
【図２】
2015186515.tif 000017
【図３】
2015186515.tif 000018
【図４】
2015186515.tif 000019
【図５】
2015186515.tif 000020
【図７】
2015186515.tif 000021
【図８】
2015186515.tif 000022
【図９】
2015186515.tif 000023
【図１０】
2015186515.tif 000024
【図１１】
2015186515.tif 000025