(57)【特許請求の範囲】
【請求項１】
被験者の身体の状態を判定する身体能力判定装置であって、
前記被験者の身体の動きを計測する身体状態計測センサから供給されるセンシング情報に基づいて、前記被験者の身体の動きの特徴を示す運動特徴量を算出する運動特徴量観測部と、
前記運動特徴量と身体能力の評価値との対応を特定するための基準情報を格納する基準情報格納部と、
前記基準情報格納部に格納されている前記基準情報に基づいて、前記運動特徴量観測部で算出された前記被験者の運動特徴量から、前記被験者の身体能力を示す評価値を特定する身体能力判定部とを、
有し、
前記身体状態計測センサが、４点の荷重センサを備え、前記被験者の足踏み運動時の前後左右の重心移動を計測し、
前記運動特徴量から特定されるリャプノフ指数と身体能力の評価値との対応が前記基準情報に記載されており、
前記身体能力判定部が、前記被験者の運動特徴量からカオス力学系の計算技術を用いて発散性を示すリャプノフ指数を算出し、前記基準情報に基づいて、前記被験者のリャプノフ指数から前記被験者の身体能力を示す評価値を特定し、
前記運動特徴量が、時系列カオス解析によって得られるアトラクタであり、再構成された前記アトラクタから、Sano-Sawada法を用いて前記リャプノフ指数を推定し、前記リャプノフ指数を推定する際に、超球サイズ、近傍点数、発展時間係数をパラメータとする身体能力判定装置。
【請求項２】
前記基準情報格納部に複数の基準情報が格納されており、
前記センシング情報に基づいて、トレンド除去変動解析法を利用した解析によりスケーリング指数を算出するトレンド除去変動解析部を有し、
前記身体能力判定部が、前記トレンド除去変動解析部で算出された前記スケーリング指数に応じた基準情報を選択し、前記選択された基準情報に基づいて、前記被験者の身体能力を示す評価値を特定するよう構成されている請求項１に記載の身体能力判定装置。
【請求項３】
ターケンス埋め込み定理を用いて前記アトラクタを再構成するためのパラメータを探索する請求項１に記載の身体能力判定装置。
【請求項４】
前記運動特徴量として用いられる時系列カオス解析によって得られるアトラクタと、トレンド除去変動解析法を利用した解析によって得られるスケーリング指数と、前記被験者の運動特徴量から特定される前記被験者の身体能力を示す前記評価値とを、リアルタイムに１つの表示画面内に表示する請求項１から３のいずれか１つに記載の身体能力判定装置。
【請求項５】
被験者の身体能力を示す評価値を出力する装置におけるデータ処理方法であって、
前記装置は、運動特徴量観測部、基準情報格納部、評価値特定部を有し、
前記運動特徴量観測部が、前記被験者の身体の動きを計測する身体状態計測センサから供給されるセンシング情報に基づいて、前記被験者の身体の動きの特徴を示す運動特徴量を算出するステップと、
前記評価値特定部が、前記基準情報格納部に格納されている前記運動特徴量と身体能力の評価値との対応を特定するための基準情報に基づいて、前記運動特徴量観測部で算出された前記被験者の運動特徴量から、前記被験者の身体能力を示す評価値を特定するステップとを、
有し、
前記身体状態計測センサが、４点の荷重センサを備え、前記被験者の足踏み運動時の前後左右の重心移動を計測し、
前記運動特徴量から特定されるリャプノフ指数と身体能力の評価値との対応が前記基準情報に記載されており、
前記評価値特定部が、前記被験者の運動特徴量からカオス力学系の計算技術を用いて発散性を示すリャプノフ指数を算出し、前記基準情報に基づいて、前記被験者のリャプノフ指数から前記被験者の身体能力を示す評価値を特定し、
前記運動特徴量が、時系列カオス解析によって得られるアトラクタであり、再構成された前記アトラクタから、Sano-Sawada法を用いて前記リャプノフ指数を推定し、前記リャプノフ指数を推定する際に、超球サイズ、近傍点数、発展時間係数をパラメータとするデータ処理方法。
【請求項６】
前記装置が、前記センシング情報に基づいて、トレンド除去変動解析法を利用した解析によりスケーリング指数を算出するトレンド除去変動解析部を有し、
前記評価値特定部が、前記トレンド除去変動解析部で算出された前記スケーリング指数に応じた基準情報を複数の基準情報の中から選択し、前記選択された基準情報に基づいて、前記被験者の身体能力を示す評価値を特定する請求項５に記載のデータ処理方法。
【請求項７】
ターケンス埋め込み定理を用いて前記アトラクタを再構成するためのパラメータを探索する請求項５に記載のデータ処理方法。
【請求項８】
前記運動特徴量として用いられる時系列カオス解析によって得られるアトラクタと、トレンド除去変動解析法を利用した解析によって得られるスケーリング指数と、前記被験者の運動特徴量から特定される前記被験者の身体能力を示す前記評価値とを、リアルタイムに１つの表示画面内に表示する請求項５から７のいずれか１つに記載のデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被験者の身体能力を評価する身体能力判定装置及びデータ処理方法に関する。
【背景技術】
【０００２】
ロコモティブシンドローム（運動器症候群）は、メタボリックシンドローム（代謝症候群）や認知症と並び、健康寿命の短縮や寝たきり状態又は要介護状態を誘発する３大要因の１つとして注目されており、この施策が急務とされている。
【０００３】
身体能力の低下には、筋力の低下や骨格の劣化と共に、脳機能の劣化という要因も存在する。筋力や骨格の強化に関しては整形外科の知識を必要とする一方、脳機能の劣化に関しては脳科学の知識が必要となり、整形外科及び脳科学の両方の観点から横断的に評価（診断）及び施術指導を行う必要がある。
【０００４】
身体の状態を測定する技術として、身体重心動揺や歩行運動のゆらぎを計測する技術などが存在する。例えば、非特許文献１では、下顎位の変化がヒトの直立姿勢時の身体重心動揺のゆらぎや筋肉の活動にどのような影響を及ぼすのかについて検討がなされている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】「下顎位の変化による身体重心動揺のゆらぎ解析」岐阜歯科学会３５巻３号１６７〜１８１、２００９年２月
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の技術においては、整形外科及び脳科学の両方の観点から、人間の身体能力を効率的に判定するシステム、例えば、脳機能劣化と整形劣化を同時計測でき、かつ低コストで簡便であるシステムは確立されていない。
【０００７】
上記の問題を解決するため、本発明は、身体能力を効率的に判定するための身体能力判定装置及びデータ処理方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記の目的を達成するため、本発明の身体能力判定装置は、被験者の身体の状態を判定する身体能力判定装置であって、
前記被験者の身体の動きを計測する身体状態計測センサから供給されるセンシング情報に基づいて、前記被験者の身体の動きの特徴を示す運動特徴量を算出する運動特徴量観測部と、
前記運動特徴量と身体能力の評価値との対応を特定するための基準情報を格納する基準情報格納部と、
前記基準情報格納部に格納されている前記基準情報に基づいて、前記運動特徴量観測部で算出された前記被験者の運動特徴量から、前記被験者の身体能力を示す評価値を特定する身体能力判定部とを、
有し、
前記身体状態計測センサが、４点の荷重センサを備え、前記被験者の足踏み運動時の前後左右の重心移動を計測し、
前記運動特徴量から特定されるリャプノフ指数と身体能力の評価値との対応が前記基準情報に記載されており、
前記身体能力判定部が、前記被験者の運動特徴量からカオス力学系の計算技術を用いて発散性を示すリャプノフ指数を算出し、前記基準情報に基づいて、前記被験者のリャプノフ指数から前記被験者の身体能力を示す評価値を特定し、
前記運動特徴量が、時系列カオス解析によって得られるアトラクタであり、再構成された前記アトラクタから、Sano-Sawada法を用いて前記リャプノフ指数を推定し、前記リャプノフ指数を推定する際に、超球サイズ、近傍点数、発展時間係数をパラメータとするものである。
【０００９】
また、上記の目的を達成するため、本発明のデータ処理方法は、被験者の身体能力を示す評価値を出力する装置におけるデータ処理方法であって、
前記装置は、運動特徴量観測部、基準情報格納部、評価値特定部を有し、
前記運動特徴量観測部が、前記被験者の身体の動きを計測する身体状態計測センサから供給されるセンシング情報に基づいて、前記被験者の身体の動きの特徴を示す運動特徴量を算出するステップと、
前記評価値特定部が、前記基準情報格納部に格納されている前記運動特徴量と身体能力の評価値との対応を特定するための基準情報に基づいて、前記運動特徴量観測部で算出された前記被験者の運動特徴量から、前記被験者の身体能力を示す評価値を特定するステップとを、
有し、
前記身体状態計測センサが、４点の荷重センサを備え、前記被験者の足踏み運動時の前後左右の重心移動を計測し、
前記運動特徴量から特定されるリャプノフ指数と身体能力の評価値との対応が前記基準情報に記載されており、
前記評価値特定部が、前記被験者の運動特徴量からカオス力学系の計算技術を用いて発散性を示すリャプノフ指数を算出し、前記基準情報に基づいて、前記被験者のリャプノフ指数から前記被験者の身体能力を示す評価値を特定し、
前記運動特徴量が、時系列カオス解析によって得られるアトラクタであり、再構成された前記アトラクタから、Sano-Sawada法を用いて前記リャプノフ指数を推定し、前記リャプノフ指数を推定する際に、超球サイズ、近傍点数、発展時間係数をパラメータとするものである。
【発明の効果】
【００１０】
本発明は、上記の構成又は処理を有しており、被験者の身体能力を効率的に判定するために脳機能劣化と整形劣化を同時計測（記録）可能とし、かつ、低コストで簡便なシステムを構築できるという効果を有している。
【図面の簡単な説明】
【００１１】
【図１Ａ】本発明の実施の形態における身体能力判定装置の構成の一例を示すブロック図である。
【図１Ｂ】本発明の実施の形態における身体能力判定装置の身体能力判定部の構成の一例を示すブロック図である。
【図２】本発明の実施の形態における運動特徴量の計測処理の一例を示すフローチャートである。
【図３Ａ】本発明の実施の形態における運動能力の判定処理の一例を示すフローチャートである。
【図３Ｂ】本発明の実施の形態におけるゆらぎ解析処理の一例を示すフローチャートである。
【図４】本発明の実施の形態における身体能力判定処理の一例を示すフローチャートである。
【図５】本発明の実施の形態におけるゆらぎ解析処理の具体的な一例を示す図である。
【図６】本発明の実施の形態において、被験者の身体の状態を計測するための身体状態計測センサ及び計測結果の一例を示す図である。
【図７】本発明の実施の形態において、被験者の身体能力を評価するための基準情報の一例を示す図である。
【図８】本発明の別の実施の形態における身体能力判定装置の構成の一例を示すブロック図である。
【図９】本発明の別の実施の形態において、被験者の身体能力を評価するための基準情報の一例を示す図である。
【図１０】本発明の別の実施の形態において、被験者の身体能力を評価するための基準情報の別の一例を示す図である。
【図１１】本発明に係る実験において、時系列データを空間ベクトルＰとして、良好なアトラクタを再構成するためのパラメータを探索し、空間ベクトルＰをプロットすることでアトラクタを描画した一例を示すイメージ図である。
【図１２】本発明に係る実験において、埋め込み次元数を４、埋め込み遅延時間を３０（ｓ）とした場合に、再構成されたアトラクタをある平面上に投影した写像（六角形状）の一例を示す図である。
【図１３】本発明に係る実験において、６人の被験者に係る重心位置データに基づいて、埋め込み次元数を４、埋め込み遅延時間を３０（ｓ）とした場合に、再構成されたアトラクタをある平面上に投影した写像（六角形状）の一例を示す図である。
【図１４】本発明に係る実験において、アトラクタ重畳時の写像空間における時系列の点軌跡を、ｄ１からｄ３の超球サイズ範囲内で軌道を描く運動のみを対象に軌道離散を計算する様子を示すイメージ図である。
【図１５】本発明に係る実験において、超球サイズを０．１、近傍点数を２０、発展時間係数を３０としてSano-Sawada法を用いて計算されたリャプノフ指数の値を縦軸とし、被験者の年齢を横軸としたグラフの一例を示す図である。
【図１６】本発明に係る実験において、各年代の５人の被験者に係る重心位置データに基づいて、埋め込み次元数を４、埋め込み遅延時間を８（ｓ）とした場合に、再構成されたアトラクタをある平面上に投影した写像（蝶型）の一例を示す図である。
【図１７】本発明に係る実験において、超球サイズを０．０８、近傍点数を２０、発展時間係数を３０としてSano-Sawada法を用いて計算されたリャプノフ指数の値を縦軸とし、被験者の年齢を横軸としたグラフの一例を示す図である。
【図１８】本発明に係る実験において、再構成されたアトラクタの安定性の指標としてＤＦＡにより算出したスケーリング指数を縦軸とし、被験者の年齢を横軸としたグラフの一例を示す図である。
【図１９】本発明において、リャプノフ指数、ＤＦＡのスケーリング指数、アトラクタの形状やリャプノフ指数などから得られる評価値（脳機能年齢）をリアルタイムに１画面に表示した場合の表示画面の一例を示す図である。
【発明を実施するための形態】
【００１２】
以下、図面を参照しながら、本発明の実施の形態について説明する。
【００１３】
本発明は、被験者の身体能力を判定することで、被験者の現在の身体能力の評価を行うことが可能な身体能力判定装置及びデータ処理方法を提供する。なお、本発明は、簡素な構成によって実現可能であるよう構成されており、また、被験者にとって負担にならない簡単な方法で被験者の身体能力の脳機能劣化と整形劣化を同時計測（記録）可能としし、被験者の身体能力を判定することが可能である。
【００１４】
まず、本発明の実施の形態における身体能力判定装置の構成について説明する。図１Ａは、本発明の実施の形態における身体能力判定装置の構成の一例を示すブロック図であり、図１Ｂは、本発明の実施の形態における身体能力判定装置の身体能力判定部の構成の一例を示すブロック図る。図１Ａに図示されている身体能力判定装置１００は、運動特徴量観測部１１０、身体能力判定部１２０、基準情報格納部１３０を有している。なお、図１Ａ及び図１Ｂには、各機能がブロックによって図示されているが、これらの各機能は、ハードウェア及び／又は（プログラムをコンピュータによって実行可能なプログラム）によって実現可能である。
【００１５】
運動特徴量観測部１１０は、被験者の身体の状態（被験者の身体の動き）に係る情報から、被験者の運動特徴量を計測する機能を有している。運動特徴量観測部１１０は、例えば、身体状態計測センサ１９１から供給される被験者の身体の状態を検出したセンシング情報に基づいて、被験者の運動特徴量を算出することが可能である。また、運動特徴量観測部１１０は、力学系の技術を用いた場合には、被験者の運動特徴量をアトラクタ化する（例えば、後述する図２の処理を行う）ことも可能である。
【００１６】
なお、身体状態計測センサ１９１として、被験者の身体における筋肉運動を数値として計測する任意のセンサが利用可能である。例えば、荷重・加速度センサを使って被験者の直立姿勢時の重心運動（重心動揺）を検出する方法や、被験者が歩行運動を行う際の身体の動きや歩調を検出する方法を利用することが可能である。また、ビデオカメラなどによるモーションキャプチャから体の特徴点の軌道をセンシングして被験者の直立姿勢時の重心運動（重心動揺）を検出する方法や、被験者が歩行運動を行う際の身体の動きや歩調を検出する方法を利用することが可能である。図６は、本発明の実施の形態において、被験者の歩行運動に伴って検出されるセンシング情報の一例を示す図である。図６（Ａ）に図示されているように、筋活動電位、３軸加速度変化、温度変化を計測することが可能なワイアレス生体センサ６００を被験者に装着してもらい（図６（Ａ）では、頭部、腰部、両膝部に装着）、被験者に歩行運動を行ってもらうことによって歩行運動に伴った被験者の動きを計測する。この計測によって、例えば図６（Ｂ）に図示されているように、歩行時の被験者の頭部、腰部、両膝部の位置の変動、被験者の重心運動、被験者の歩調などが計測可能である。なお、ここでは、主に、被験者に歩行運動を行ってもらう場合を一例として説明を行うが、例えば、図６（Ｃ）に図示されているように、被験者を直立させて被験者の重心運動（垂直重心動揺）を荷重センサ７００によって計測することも可能である。
【００１７】
また、身体能力判定部１２０は、運動特徴量観測部１１０で計測された被験者の運動特徴量に基づいて、被験者の身体能力を判定する機能を有している。身体能力判定部１２０は、例えば、被験者の運動特徴量から得られる数値（例えば、後述のリャプノフ指数）と、基準情報格納部１３０に格納されている基準情報とを比較することで、被験者の運動特徴量がどのレベルに位置するかを特定することが可能である。なお、身体能力判定部１２０による判定結果は、例えば評価値としてモニタ／スピーカ１９２などの出力装置から出力可能である。
【００１８】
身体能力判定部１２０は、図１Ｂに図示されているように、脳機能判定部１２１及び整形判定部１２２を有している。脳機能判定部１２１は、本発明の実施の形態において説明する動作（運動特徴量から身体能力の評価を行う動作）を実行する機能を有しており、一方、整形判定部１２２は、骨格劣化・筋力劣化を整形外科の観点で観測する機能を有している。例えば骨格劣化に関しては、右足・左足の足踏み運動の劣化を観測することで、背骨の歪みを観測することが可能である。
【００１９】
また、基準情報格納部１３０は、身体能力判定部１２０による判定の際に参照される基準情報を事前に格納する機能を有している。基準情報としては、例えば、被験者の運動特徴量から得られる数値（例えば、後述のリャプノフ指数）に対応する評価値（例えば、後述の脳機能年齢）を規定した対応表などを用いることが可能である。また、例えばアトラクタ群（アトラクタ行列）の各パターンを保持し、当該各パターンに対応する評価値を規定しておいてもよい。
【００２０】
次に、図１Ａ及び図１Ｂに図示されている身体能力判定装置１００の構成に基づき、本発明に係る処理について説明する。
【００２１】
図２は、本発明の実施の形態における運動特徴量の計測処理の一例を示すフローチャートである。なお、図２に図示されている運動特徴量の計測処理は、図１Ａに図示されている運動特徴量観測部１１０によって、ミリ秒オーダ（ミリ秒単位）で実行される処理である。
【００２２】
図２において、運動特徴量観測部１１０は、重心運動の計測値を取得し（ステップＳ２０１）、歩調の計測値を取得し（ステップＳ２０２）、３軸加速度変化の計測値（３次元空間における加速度の値）を取得する（ステップＳ２０３）。これらのステップＳ２０１〜Ｓ２０３の各センシング情報は、身体状態計測センサ１９１から、例えばミリ秒単位で運動特徴量観測部１１０へ供給される。なお、ステップＳ２０１〜Ｓ２０３の各処理は並列処理である。また、本発明の実施の形態では、重心運動、歩調、３軸加速度変化を計測する場合について説明するが、これらの計測値のうちのいずれか１つ又は複数（あるいはすべて）を利用してもよく、また、その他のセンシング情報を利用してもよい。
【００２３】
ここで、運動特徴量観測部１１０は、身体状態計測センサ１９１から供給される計測値をアトラクタ化する。運動特徴量観測部１１０は、例えば、重心運動の計測値をアトラクタ化し（ステップＳ２０４）、歩調の計測値をアトラクタ化し（ステップＳ２０５）、３軸加速度変化の計測値をアトラクタ化する（ステップＳ２０６）。
【００２４】
以上の処理のように、運動特徴量観測部１１０は、ミリ秒オーダで、被験者の運動特徴量を計測してアトラクタ化することが可能である。
【００２５】
また、図３Ａは、本発明の実施の形態における身体能力の判定処理の一例を示すフローチャートである。なお、図３Ａに図示されている身体能力の判定処理は、図１Ａ及び図１Ｂに図示されている身体能力判定部１２０によって、数秒〜１分オーダ（数秒〜１分単位）で実行される処理である。
【００２６】
図３Ａにおいて、身体能力判定部１２０は、直近の被験者の運動特徴量を運動特徴量観測部１１０から取得する（ステップＳ３１０）。なお、ステップＳ３１０において取得する被験者の運動特徴量は、図２に図示されているステップＳ２０４〜Ｓ２０６で算出された結果（各運動特徴量がアトラクタされたもの）である。次に、身体能力判定部１２０は、取得した被験者の運動特徴量のゆらぎ解析処理を行う（ステップＳ３３０）。なお、ステップＳ３３０におけるゆらぎ解析処理において、身体能力判定部１２０は、例えば、各運動特徴量からリャプノフ指数を算出する処理を行う。ステップＳ３３０におけるゆらぎ解析処理の詳細については、図３Ｂを参照しながら後述する。そして、身体能力判定部１２０は、ステップＳ３３０におけるゆらぎ解析処理の処理結果（具体的には、算出されたリャプノフ指数）の遷移をログとして蓄積（ロギング）する（ステップＳ３５０）。
【００２７】
また、図３Ｂは、本発明の実施の形態におけるゆらぎ解析処理の一例を示すフローチャートである。なお、図３Ｂに図示されているゆらぎ解析処理は、上述のように、図３ＡのステップＳ３３０におけるゆらぎ解析処理の詳細を説明するものである。
【００２８】
図３Ｂにおいて、身体能力判定部１２０は、重心運動のアトラクタ群（アトラクタ行列）を取得して（ステップＳ３３１１）、重心運動のリャプノフ指数を算定し（ステップＳ３３１２）、過去のアトラクタ群からリャプノフ指数の安定傾向を解析する（ステップＳ３３１３）。
【００２９】
また、身体能力判定部１２０は、歩調のアトラクタ群（アトラクタ行列）を取得して（ステップＳ３３２１）、歩調のリャプノフ指数を算定し（ステップＳ３３２２）、過去のアトラクタ群からリャプノフ指数の安定傾向を解析する（ステップＳ３３２３）。
【００３０】
また、身体能力判定部１２０は、３軸加速度変化の計測値のアトラクタ群（アトラクタ行列）を取得して（ステップＳ３３３１）、３軸加速度変化の計測値のリャプノフ指数を算定し（ステップＳ３３３２）、過去のアトラクタ群からリャプノフ指数の安定傾向を解析する（ステップＳ３３３３）。
【００３１】
なお、本発明の実施の形態におけるゆらぎ解析処理の具体的な計算については、例えば、図５に図示されているように実行することが可能である。図５は、本発明の実施の形態におけるゆらぎ解析処理の具体的な計算方法の一例を示す図である。なお、ゆらぎ解析処理の計算方法は、図５に図示されているものに限定されるものではない。複雑系の解析方法については、様々な研究が行われており、現在及び今後確立される任意の解析方法を本発明に適用することが可能である（例えば、此処まで来た複雑系解析ツール、http://www.ieice.org/cs/csbn/program/papers/04_1_miao.pdfを参照）。
【００３２】
本発明の実施の形態では、例えば、重心運動、歩調、３軸加速度変化などの各計測値を取得する（図２のステップＳ２０１〜Ｓ２０３）。こうした計測値（センシングデータ）は、例えば図５（Ａ）に図示されているように、時間軸に対して変動する曲線として取得され、被験者の身体の状態の特徴を表す運動特徴量として取り扱うことが可能である。
【００３３】
ここで、例えば、所定の微小ウィンドウ（時間幅τ）の範囲をずらしながら、微小ウィンドウ内における運動特徴量の局所的な傾向（ローカルトレンド）のアトラクタ（軌道）を取得することで、運動特徴量のアトラクタ群（アトラクタ行列）を取得し（図３ＢのＳ３３１１、Ｓ３３２１、Ｓ３３３１、図５（Ｂ））、このようにして得られたアトラクタ群（アトラクタ行列）からリャプノフ指数を算定することが可能である（図３ＢのＳ３３１２、Ｓ３３２２、Ｓ３３３２、図５（Ｃ））。なお、このとき、算定されるリャプノフ指数は１つのみであるが、所定の微小ウィンドウ（時間幅τ）を複数個含むウィンドウ（メタウィンドウ）をずらすことによって、リャプノフ指数の時系列データ（時間的に異なる複数個のリャプノフ指数）を得ることができる。
【００３４】
なお、身体能力判定部１２０は、判定に用いる運動特徴量として、上述のようにリャプノフ指数（あるいは、リャプノフ指数の時系列データ）を出力してもよく、また、アトラクタ群（アトラクタ行列）そのものを出力してもよい。例えば、基準情報として、被験者の運動特徴量から得られる数値（例えば、リャプノフ指数）に対応する評価値を規定した対応表を用いる場合には、判定に用いる運動特徴量としてリャプノフ指数を出力し、また、評価値に対応するアトラクタ群（アトラクタ行列）の各パターンを保持している場合には、アトラクタ群（アトラクタ行列）そのものを出力してもよい。
【００３５】
また、図４は、本発明の実施の形態における身体能力判定処理の一例を示すフローチャートである。なお、図４に図示されている身体能力判定処理は、図１Ａ及び図１Ｂに図示されている身体能力判定部１２０によって、ゆらぎ解析処理の解析結果（図３ＡのステップＳ３３０の処理結果）が取得されるたびに実行されてもよく、あるいは、数十秒から数分のオーダで被験者の運動特徴量を計測した後に実行されてもよい。
【００３６】
図４において、身体能力判定部１２０は、基準情報格納部１３０に格納されている基準情報を取得するとともに（ステップＳ４０１）、身体能力判定部１２０によって算出された判定に用いる運動特徴量（例えば、運動特徴量に基づくリャプノフ指数によって表される）を取得する（ステップＳ４０２）。なお、ステップＳ４０１において取得する基準情報は、基準情報格納部１３０に格納されているリャプノフ指数と評価値との対応表や、評価値に対応するアトラクタ群（アトラクタ行列）の各パターンなどである。また、ステップＳ４０２において取得する運動特徴量は、図３ＡのステップＳ３３０のゆらぎ解析処理の解析結果（図３ＡのステップＳ３３０の処理結果）であり、具体的には、各運動特徴量から取得されるアトラクタ群や各運動特徴量から算出されたリャプノフ指数である。
【００３７】
身体能力判定部１２０は、身体能力判定部１２０によって算出された判定に用いる運動特徴量と基準情報とに基づいて、被験者の身体能力の評価値を特定する（ステップＳ４０３）。具体的には、身体能力判定部１２０は、身体能力判定部１２０によって算出された判定に用いる運動特徴量が、基準情報においてどのレベルに対応するものであるかを判定する。この判定では、被験者の運動特徴量から得られる数値（例えば、リャプノフ指数）に対応する評価値（例えば、後述の脳機能年齢）を規定した対応表を用いる場合には、被験者の運動特徴量から得られるリャプノフ指数に対応する評価値を特定し、また、評価値に対応するアトラクタ群（アトラクタ行列）の各パターンを用いる場合には、被験者の運動特徴量から得られるアトラクタ群のパターンを解析して、どの評価値に対応するアトラクタ群（アトラクタ行列）のパターンに類似しているかを特定する。また、パターンの類似を調べる場合には、アトラクタの特定の局所領域の広がり幅や曲面形状などを、基準となる評価値に対応するアトラクタの形状と比較してもよい。
【００３８】
図７には、リャプノフ指数に対応する評価値（ここでは、脳機能年齢と記載）を規定した対応表（基準情報の一例）が図示されている。なお、一般的に、リャプノフ指数の値が大きいほど、身体能力を制御する脳の機能は活性化されており、リャプノフ指数の値が小さいほど、身体能力を制御する脳の機能は低減している傾向にあるが、後述のように（別の実施の形態として説明）、被験者の年齢（年代）に応じて利用する対応表を適宜選択することで、より精度の高い評価値を出力することが可能である。
【００３９】
身体能力判定部１２０は、例えば図７に図示されている対応表を用いて、リャプノフ指数に対応する評価値（脳機能年齢）を特定し、この評価値（脳機能年齢）を、モニタ／スピーカ１９２などから出力したり、通信部（図１Ａには不図示）を介して配信先（サーバなど）へ送信したりすることが可能である。また、被験者の登録（ユーザ登録）を行うことで被験者の評価値を逐次蓄積し、被験者の脳機能の状態をトレースしたり、被験者による整形外科施術（リハビリテーションなど）の進行状況などの把握に活用したりすることも可能である。
【００４０】
次に、本発明の別の実施の形態について説明する。この実施の形態によれば、ＤＦＡ（Detrended Fluctuation Analysis：トレンド除去変動解析法）解析を用いて推定した被験者の年齢（年代）を反映させることで、被験者の身体能力をより詳細に評価できるようになる。なお、上述のように、リャプノフ指数の値が大きいほど身体能力を制御する脳の機能は活性化されており、リャプノフ指数の値が小さいほど身体能力を制御する脳の機能は低減している傾向にあるが、被験者の年齢（年代）によって、リャプノフ指数の値に対応する評価値は異なっている。従って、被験者の年齢（年代）に応じた基準情報を選択することが有用である。
【００４１】
図８は、本発明の別の実施の形態における身体能力判定装置の構成の一例を示すブロック図である。図７に図示されている身体能力判定装置１００は、運動特徴量観測部１１０、身体能力判定部１２０、基準情報格納部１３０、ＤＦＡ解析部１４０を有している。なお、図８には、各機能がブロックによって図示されているが、これらの各機能は、ハードウェア及び／又は（プログラムをコンピュータによって実行可能なプログラム）によって実現可能である。
【００４２】
図８に図示されている身体能力判定装置１００は、図１Ａに図示されている構成と比較して、ＤＦＡ解析部１４０を更に有している。また、身体能力判定部１２０に機能が追加され、さらに、基準情報格納部１３０には、より詳細な評価を可能とする基準情報が格納される。以下、主に、本発明の別の実施の形態に特徴的な機能及び動作について説明する。
【００４３】
ＤＦＡ解析部１４０は、被験者の身体の状態（被験者の身体の動き）に係る情報（例えば、身体状態計測センサ１９１からのセンシング情報）に基づいてＤＦＡ解析を行い、そのＤＦＡ解析結果（スケーリング指数α）を出力する機能を有している。なお、ＤＦＡ解析は、非定常的な時系列データの長期相関特性を解析する手法であり、広く知られている技術なので、ここでは詳細な説明は省略する。ＤＦＡ解析は、入力データを積分することで作成した新たな時系列データを長さｎのボックスに分割し、各ボックスにおいて最小自乗法による近似直線を計算してローカルトレンドを求め、信号からトレンドを除去してから分散を平均した値（Ｆ（ｎ））を算出するものである。このＤＦＡ解析によって求められるＦ（ｎ）は、Ｆ（ｎ）∝ｎαによって表すことが可能であり、スケーリング指数αによって特徴付けられる。
【００４４】
本発明の発明者は、ＤＦＡ解析によって得られるスケーリング指数が被験者の年齢（年代）の推測に有用であり、推測された被験者の年齢（年代）に基づいて被験者の身体能力の判定を行うことが有用であることを発見した。具体的には、例えばα＝０．８を閾値として、α≧０．８とα＜０．８とで年代を分けることが可能である。以下では、α≧０．８と算出された場合を若年層、α＜０．８と算出された場合を老年層と呼ぶ。なお、ここでは一例として、１つの閾値（例えば、α＝０．８）によって年代を２つに分けているが（若年層と老年層の２つのカテゴリに分類）、閾値として設定される値や閾値の個数及び分類するカテゴリ数は、適宜設定することが可能である。
【００４５】
ＤＦＡ解析部１４０によるＤＦＡ解析結果は、身体能力判定部１２０へ供給される。身体能力判定部１２０は、運動特徴量観測部１１０で計測された被験者の運動特徴量に基づいて、被験者の身体能力を示す評価値を特定する機能を有しており、例えば、被験者の運動特徴量から得られる数値と、基準情報格納部１３０に格納されている基準情報とを比較することで、被験者の運動特徴量がどのレベルに位置するかを評価することが可能である
【００４６】
本実施の形態では、身体能力判定部１２０は、ＤＦＡ解析部１４０のＤＦＡ解析結果に基づいて、基準情報格納部１３０に格納されている基準情報を選択する機能を更に有している。身体能力判定部１２０は、例えばＤＦＡ解析部１４０から供給されたスケーリング指数αの値に基づいて、基準情報格納部１３０に格納されている複数の基準情報の中から、スケーリング指数αの値に対応する基準情報を選択し、被験者の運動特徴量と、選択された基準情報とを比較することで、被験者の身体能力の評価値を特定する。
【００４７】
例えば、基準情報格納部１３０には、若年層用及び老年層用の２つの基準情報が格納されている。例えば、被験者の身体能力の評価にリャプノフ指数が用いられる場合には、図９に図示されているような若年層用の対応表と、図１０に図示されているような老年層用の対応表とが格納されている。また、例えば、被験者の身体能力の評価にアトラクタ群のパターンの類似性の判断が行われる場合には、若年層用のアトラクタ群と老年層用のアトラクタ群とが格納されている。
【００４８】
上述の例に示されているように、ＤＦＡ解析部１４０から供給されたスケーリング指数αが０．８以上の場合には、身体能力判定部１２０は、若年層用の対応表（図９）を基準情報として選択し、選択された基準情報と比較することで、被験者の身体能力の評価値を特定する。同様に、ＤＦＡ解析部１４０から供給されたスケーリング指数αが０．８未満の場合には、身体能力判定部１２０は、老年層用の対応表（図１０）を基準情報として選択し、選択された基準情報と比較することで、被験者の身体能力の評価値を特定する。
【００４９】
以上のように、本実施の形態によれば、ＤＦＡ解析部１４０によるＤＦＡ解析結果に基づいて被験者の年齢（年代）を推定し、推定された被験者の年齢（年代）に応じた基準情報を選択し、選択された基準情報に基づいて、被験者の身体能力を評価することが可能である。
【００５０】
次に、上述の本発明の実施の形態に関連する実施例として、試作システムを用いて行われた実験の詳細について説明する。
【００５１】
本発明に係る実験では、図１に図示されている身体状態計測センサ１９１として、図６（Ｃ）に図示されているような４点の荷重センサ（例えば、ストレインゲージ式フォースセンサ）を備えた板状コントローラを用いて被験者の重心動揺を計測した。なお、リハビリテーション医学やスポーツ医学などの分野で一般的に使用されている重心動揺計は、重心動揺の正常／異常を判定するために微小なバランスのぶれを抽出する必要があることから、高精度・高感度の荷重センサを具備するとともに床を伝搬する振動などのノイズを除去するためのハードウェア機構を持ち、非常に高価な装置となっている。一方、本発明に係る実験では、例えばゲーム用コントローラとして広く普及しているものなどのような極めて簡易かつ安価な板状ボードコントローラ（以下、民生用ボードコントローラと記載する）を用いた。
【００５２】
本発明に係る実験では、被験者に民生用ボードコントローラ上で足踏み運動をさせ、民生用ボードコントローラから得られた重心動揺データから、床平面上の重心移動ベクトルを算出した。次に、得られた重心移動ベクトルに関し、民生用ボードコントローラの前後左右２軸の成分の時系列データを使ってアトラクタの再構成を試みた。その後、再構成されたアトラクタの重畳の自己相関の不安定の度合いを示すリャプノフ指数を推定し、この数値を脳活性度を示す生体ゆらぎの定量的な指標として用いた。最後に、足踏み運動における運動モーメントの安定性を定量的に確認するために、非定常的な信号の解析に対して有効であるスペクトル解析手法（ＤＦＡ）を用いてスケーリング指数を算出した。なお、ＤＦＡで算出されたスケーリング指数は、運動モーメントに起因する歩行の安定性を示す指標として取り扱うことができ、安定性が高いほどスケーリング指数は低い値となる。
【００５３】
足踏み動作時の足踏み位置やその荷重は時間と共に複雑に変化し、そのデータは非線形の時系列データとして計測可能である。しかしながら、計測環境の影響によるノイズや被験者の年齢、体調など生態的な影響も加わることから、その複雑な挙動を詳細に解析することは容易ではない。時系列カオス解析手法は、このような現実世界の複雑で不規則な時系列信号から、決定論的性質を単純な非線形信号モデルとして導き出す解析手法であり、具体的には、非線形、非定常な波形信号から自己相関性（傾向）を持った力学運動の軌跡を可視化（アトラクタ再構成）し、その傾向特性を定量化する手法である。
【００５４】
本発明に係る実験で用いられる民生用ボードコントローラから得られたデータは、計測時に接地面から伝搬する様々なノイズを含んでおり、かつ低解像度である。このようなノイズが多く含まれている低解像度の重心動揺データに対しても、傾向特性の分析を行う時系列カオス解析手法は有効であると考えられるが、本発明に係る実験に先立って、まず、このようなノイズが多く含まれている重心動揺データからのアトラクタ再構成の可否について調べる予備実験を行った。具体的には、民生用ボードコントローラを使用して３０（Ｈｚ）のサンプリングレートで計測した足踏み動作の重心動揺データに対してターケンスの埋め込み定理を使用し、図１１に図示されているように、時系列データを空間ベクトルＰ（埋め込みベクトルＰ）として、良好なアトラクタを再構成するためのパラメータ（埋め込み次元数ｎ、埋め込み遅延時間τ）を探索した。なお、図１１のイメージ図では、３次元の遅延空間をＸＹＺの観測点とし、Ｐ１、Ｐ２、Ｐ３の順にプロット移動したときの軌跡がアトラクタとなる。
【００５５】
その結果、埋め込み次元数を４、埋め込み遅延時間を３０（ｓ）とした場合に、再構成されたアトラクタをある平面上に投影した写像が、図１２に図示されているような六角形状となった。この結果から、広い帯域において様々なノイズが多く含まれているデータに対してアトラクタの再構成が可能であることが確認できた。なお、見つけ出されたアトラクタは、足踏み運動時の運動モーメントの傾向特性を視覚的に表現していることとなる。
【００５６】
本発明に係る実験では、３０代から７０代の被験者１３名に民生用ボードコントローラ上で３分間の足踏み運動を行わせ、足踏み姿勢時の被験者の左右方向をＸ軸、前後方向をＹ軸とした重心位置データを記録して、この重心位置データから重心移動ベクトルを算出した。そして、各被験者の重心移動ベクトルのＸ軸成分及びＹ軸成分それぞれに対し、上記の予備実験で得られたパラメータ（埋め込み次元数４、埋め込み遅延時間３０（ｓ））を用いてアトラクタの再構成を行った。Ｘ軸成分についてはアトラクタの再構成が可能であり、図１３には、得られたアトラクタの代表的な形状が図示されている。なお、図１３の右下に示すアトラクタは、脳梗塞半身不随者に係る重心移動ベクトルから再構成されたものであり、他者の形状とは異なり著しく歪んでいた。その他のアトラクタは歪みながらも六角形状に近いアトラクタ形状となり、傾向特性が抽出できた。
【００５７】
また、図１４に模式的に図示されているように、力学系において接近した軌道が離れていく度合いを表す量としてリャプノフ指数が用いられる。図１４のイメージ図には、アトラクタ重畳時の写像空間における時系列の点軌跡を、ｄ１からｄ３の超球サイズ範囲内で軌道を描く運動のみを対象に軌道離散を計算する様子が図示されている。
【００５８】
本発明に係る実験では、アトラクタとして表現された時系列データの傾向特性の軌道不安定性の度合いをリャプノフ指数で示している。ここでは、リャプノフ指数の計算にリアルタイム性を考慮したSano-Sawada法を用い、超球サイズ、近傍点数、発展時間係数をパラメータとしてリャプノフ指数を調整した。具体的には、超球サイズを０．１、近傍点数を２０、発展時間係数を３０とし、Sano-Sawada法を用いてリャプノフ指数の値の計算を行った。その結果推定されたリャプノフ指数の値を縦軸、被験者の年齢を横軸としたグラフを図１５に示す。この結果からリャプノフ指数と年齢との相関はなく、加齢と共に発散する傾向が見られることから、このアトラクタ形状（六角形状）は生体ゆらぎの脳活性傾向を示すものではないと考えられる。しかしながら、一方で、上述のように脳梗塞半身不随者のアトラクタは著しく歪んでいることから、このアトラクタ形状は運動モーメントの安定傾向を示すものであると考えられる。したがって、上記の運動特徴量解析で使用したパラメータを用いた解析は、足踏み運動に関わる整形外科的な異常をアトラクタ形状によって視覚的に、あるいは、リャプノフ指数によって定量的に視覚的に把握するのに有効であると考えられる。
【００５９】
また、上記の運動特徴量解析結果を踏まえ、ターケンスの埋め込み定理を利用して、運動特徴量として有効なアトラクタを再構成するための別のパラメータの探索を新たに行った。その結果、埋め込み次元数４、埋め込み遅延時間８（ｓ）のときに、重心移動ベクトルのＸ軸成分について、すべての被験者のデータから一様に蝶型のアトラクタが確認された。図１６は、年代毎の代表的なアトラクタ形状を示す図であり、（ａ）３０代、（ｂ）４０代、（ｃ）５０代、（ｄ）６０代以上健常者、（ｅ）６０代以上非健常者の被験者の重心移動ベクトルから得られたアトラクタ形状の一例が図示されている。図１６に図示されているように、この蝶型のアトラクタ形状では、加齢と共に歪み破れが大きくなっている傾向が確認できた。また、歩行障害を持つ高齢者の足踏み運動から抽出されたアトラクタ（図１６の（ｅ））には、著しい歪みが観察できた。
【００６０】
次に、Ｘ軸成分の再構成されたアトラクタから、超球サイズを０．０８、近傍点数を２０、発展時間係数を３０としてSano-Sawada法を用いて推定したリャプノフ指数を縦軸とし、被験者の年齢を横軸としたグラフを図１７に示す。被験者の母数が少なかったことから、この結果からは強い相関は確認できなかったが、加齢と共にリャプノフ指数が小さくなる傾向は確認できた。また、再構成されたアトラクタの安定性の指標としてＤＦＡにより算出したスケーリング指数を縦軸とし、被験者の年齢を横軸としたグラフを図１８に示す。この結果によれば、加齢と共にスケーリング指数が高くなる傾向が確認できた。
【００６１】
このように、上記の運動特徴量解析で使用したパラメータ（埋め込み次元数４、埋め込み遅延時間８（ｓ））を用いた場合には、重心移動ベクトルのＸ軸成分のアトラクタの写像が、すべての被験者について一様な蝶型となったことから、足踏み運動に関する有用な傾向特性を抽出することができたと考えられる。また、再構成されたアトラクタから推定されたリャプノフ指数は、加齢に伴って低下傾向にあることが確認できた。従って、各被験者の重心移動ベクトルから抽出されたリャプノフ指数は脳活性傾向を示す定量的な指標として有効であり、上述のように、例えば脳活性傾向を示す脳機能年齢としてマッピングすることで、各被験者の脳機能に係る特徴量を評価値として特定することができると考えられる。また、脳活性度と認知能力は相関があることから、本手法で求められたリャプノフ指数は、認知症を診断するための定量的な指標となり得ると思われる。さらに、ＤＦＡにより算出されたスケーリング指数については加齢に伴って増加する傾向が見られ、年齢が高いほど整形外科的な運動モーメントのぶれが大きいことが分かった。また、歩行障害を持つ被験者のアトラクタは、その写像が他の被験者と比較して視覚的に著しい歪みがあることも確認できた頃から、本手法は整形外科的な検査の定量的指標としても用いることができると考えられる。
【００６２】
以上のように、本発明に係る実験によって、非常に簡素なセンサ（例えば、民生用ボードコントローラ）から得られる様々なノイズが含まれたデータ（例えば、重心動揺データ）であっても、カオス力学系の計算技術を用いて解析を行うことによって、アトラクタに基づいて推定されたリャプノフ指数から生体ゆらぎの不安定性（脳活性度）を把握できることが確認された。
【００６３】
また、本発明に係る実験によって、再構成されたアトラクタ（アトラクタの歪みの度合い）から運動モーメントの安定性を把握できる可能性があることも分かった。また、本発明の時系列カオス解析の手法として、特に、ターケンスの埋め込み定理を用いたアトラクタ再構成のパラメータ（埋め込み次元数、埋め込み遅延時間）の探索、Sano-Sawada法を用いたリャプノフ指数の算定（超球サイズ、近傍点数、発展時間点数をパラメータとして調整）が適用できることが確認された。また、本発明に係る実験から、ＤＦＡにより算出されたスケーリング指数が被験者の年齢（年代）の推測に有用であることが確認された。
【００６４】
なお、上記本発明の実験では、民生用ボードコントローラによるデータ収集と同時に、３次元センサ（例えば、ゲーム用の３次元モーションセンサ）からのデータを収集して被験者の骨格動作（身体姿勢）を推定し、骨格動作の推定結果と重心動揺計から得られた情報（アトラクタ、ＤＦＡ、評価値（脳機能年齢）など）をタブレットコンピュータ上に表示する試みも行われた。
【００６５】
また、民生用ボードコントローラから得られた重心動揺データ、あるいは、いったん計測が完了して記録された推進動揺データから、アトラクタを再構成する際に、図１９に図示されているように、リャプノフ指数、ＤＦＡのスケーリング指数、アトラクタの形状やリャプノフ指数などから得られる評価値（脳機能年齢）をリアルタイムに１画面に表示することで、ユーザの視認性を向上させるとともに利便性の高い表示画面を実現してもよい。さらに、３次元センサによって収集されたデータに基づく被験者の骨格動作やその他の生体情報（例えば、脈拍情報など）を同期して表示させてもよい。
【産業上の利用可能性】
【００６６】
本発明は、被験者の身体能力を効率的に判定することができるという効果、脳機能劣化と整形劣化を同時計測でき、かつ、低コストで簡便なシステムを構築できるという効果を有しており、被験者の身体能力の評価を行う技術、及び、被験者の身体能力の評価結果に基づいて機能回復の施術を行う医療技術に適用可能である。
【符号の説明】
【００６７】
１００  身体能力判定装置
１１０  運動特徴量観測部
１２０  身体能力判定部
１２１  脳機能判定部
１２２  整形判定部
１３０  基準情報格納部
１４０  ＤＦＡ解析部
１９１  身体状態計測センサ
１９２  モニタ／スピーカ
７００  荷重センサ
【図１Ａ】
0005612627.tif 000002
【図１Ｂ】
0005612627.tif 000003
【図２】
0005612627.tif 000004
【図３Ａ】
0005612627.tif 000005
【図３Ｂ】
0005612627.tif 000006
【図４】
0005612627.tif 000007
【図５】
0005612627.tif 000008
【図６】
0005612627.tif 000009
【図７】
0005612627.tif 000010
【図８】
0005612627.tif 000011
【図９】
0005612627.tif 000012
【図１０】
0005612627.tif 000013
【図１１】
0005612627.tif 000014
【図１２】
0005612627.tif 000015
【図１３】
0005612627.tif 000016
【図１４】
0005612627.tif 000017
【図１５】
0005612627.tif 000018
【図１６】
0005612627.tif 000019
【図１７】
0005612627.tif 000020
【図１８】
0005612627.tif 000021
【図１９】
0005612627.tif 000022