(57)【要約】  体動リズムの同調性を連続的に評価する。  所定の時刻における信号の位相から位相が所定角度変化した時刻、および、所定の時刻における信号の位相から位相が前記所定角度の２倍変化した時刻を算出する第１算出過程と、これらの時刻に基づいて、評価係数を算出する第２算出過程と、をそなえる。
2012036135.tif 000002
【特許請求の範囲】
【請求項１】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をそなえたことを特徴とする、情報処理方法。
【請求項２】
前記評価係数算出過程は、前記所定の時刻から前記第１時刻までの時間（以下、第１時間という）、および、前記第１時刻から前記第２時刻までの時間（以下、第２時間という）に基づいて、前記評価係数を算出する過程であることを特徴とする請求項１記載の情報処理方法。
【請求項３】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の−１倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をそなえたことを特徴とする、情報処理方法。
【請求項４】
前記評価係数算出過程は、前記所定の時刻から前記第１時刻までの時間（以下、第１時間という）、および、前記第２時刻から前記所定の時刻までの時間（以下、第２時間という）に基づいて、前記評価係数を算出する過程であることを特徴とする請求項３記載の情報処理方法。
【請求項５】
前記時刻算出過程は、ヒルベルト変換法またはパターンマッチング法を用いることで、前記第１時刻および前記第２時刻を算出する過程であることを特徴とする、請求項１〜４のいずれか１項に記載の情報処理方法。
【請求項６】
前記時刻算出過程は、前記所定の時刻における前記信号の位相から前記信号の位相が３６０度変化した前記第１時刻、および、前記所定の時刻における前記信号の位相から前記信号の位相が７２０度変化した前記第２時刻を算出する
ことを特徴とする請求項１〜５のいずれか１項に記載の情報処理方法。
【請求項７】
前記時刻算出過程は、前記所定の時刻における前記信号の位相から前記信号の位相が１８０度変化した前記第１時刻、および、前記所定の時刻における前記信号の位相から前記信号の位相が３６０度変化した前記第２時刻を算出する
ことを特徴とする請求項１〜５のいずれか１項に記載の情報処理方法。
【請求項８】
前記時刻算出過程は、一つの検出部により検出された、生体の繰り返し随意運動を表す信号から、前記第１時刻、および、前記第２時刻を算出する
ことを特徴とする請求項１〜７のいずれか１項に記載の情報処理方法。
【請求項９】
前記評価係数算出過程により算出された前記評価係数が、所定時間にわたって所定範囲内であるかを判定する判定過程をさらにそなえたことを特徴とする、請求項１〜８のいずれか１項に情報処理方法。
【請求項１０】
前記繰り返し随意運動が、歩行、ジョギング、ランニング、自転車走行、水泳、又は体操のいずれかであることを特徴とする、請求項１〜９のいずれか１項に記載の情報処理方法。
【請求項１１】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、
前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、をそなえたことを特徴とする、情報処理装置。
【請求項１２】
前記評価係数算出部は、前記所定の時刻から前記第１時刻までの時間（以下、第１時間という）、および、前記第１時刻から前記第２時刻までの時間（以下、第２時間という）に基づいて、前記評価係数を算出することを特徴とする請求項１１記載の情報処理装置。
【請求項１３】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の−１倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、
前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、をそなえたことを特徴とする、情報処理装置。
【請求項１４】
前記評価係数算出部は、前記所定の時刻から前記第１時刻までの時間（以下、第１時間という）、および、前記第２時刻から前記所定の時刻までの時間（以下、第２時間という）に基づいて、前記評価係数を算出することを特徴とする請求項１３記載の情報処理装置。
【請求項１５】
前記時刻算出部は、ヒルベルト変換法またはパターンマッチング法を用いることで、前記第１時刻および前記第２時刻を算出する過程であることを特徴とする、請求項１１〜１４のいずれか１項に記載の情報処理装置。
【請求項１６】
前記第１角度は、３６０度である
ことを特徴とする請求項１１〜１５のいずれか１項に記載の情報処理装置。
【請求項１７】
前記第１角度は、１８０度である
ことを特徴とする請求項１１〜１５のいずれか１項に記載の情報処理装置。
【請求項１８】
前記繰り返し随意運動を表す信号は、一つの検出部により検出された信号である
ことを特徴とする請求項１１〜１７のいずれか１項に記載の情報処理装置。
【請求項１９】
前記評価係数算出部により算出された前記評価係数が、所定時間にわたって所定範囲内であるかを判定する判定部をさらにそなえたことを特徴とする、請求項１１〜１８のいずれか１項に情報処理装置。
【請求項２０】
前記繰り返しリズム運動が、歩行、ジョギング、ランニング、自転車走行、水泳、又は体操のいずれかであることを特徴とする、請求項１１〜１９のいずれか１項に記載の情報処理装置。
【請求項２１】
所定の時刻における生体の繰り返し随意運動を表す信号の位相から前記信号の位相が所定角度変化した時刻、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻に基づく、前記繰り返し随意運動の同調性を評価するための評価係数を、前記評価係数を識別可能な形態にて出力する出力部をそなえたことを特徴とする出力装置。
【請求項２２】
所定の時刻における生体の繰り返し随意運動を表す信号の位相から前記信号の位相が所定角度変化した時刻、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の−１倍変化した時刻に基づく、前記繰り返し随意運動の同調性を評価するための評価係数を、前記評価係数を識別可能な形態にて出力する出力部をそなえたことを特徴とする出力装置。
【請求項２３】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、
前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、
前記評価係数算出部により算出された評価係数を、前記評価係数を識別可能な形態にて出力する出力部と、をそなえた
ことを特徴とする、情報処理システム。
【請求項２４】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の−１倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、
前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、
前記評価係数算出部により算出された評価係数を、前記評価係数を識別可能な形態にて出力する出力部と、をそなえた
ことを特徴とする、情報処理システム。
【請求項２５】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をコンピュータに実行させることを特徴とする、情報処理用プログラム。
【請求項２６】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の−１倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をコンピュータに実行させることを特徴とする、情報処理用プログラム。
【請求項２７】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記第１角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する第２算出過程と、をコンピュータに実行させることを特徴とする、情報処理用プログラムを記録したコンピュータ読取可能な記録媒体。
【請求項２８】
生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記第１角度の−１倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、
前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する第２算出過程と、をコンピュータに実行させることを特徴とする、情報処理用プログラムを記録したコンピュータ読取可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体に関する。
【背景技術】
【０００２】
日常生活における体の動きを連続的にモニタすることは、健康管理の観点から極めて重要である。また、外傷または脳神経系の疾患により、患者の歩行または走行等の規則的な体動リズムがどのような影響を受けているかを知ることは、病態を正確に把握する上で非常に役立つ。例えば、パーキンソン病の患者には１日のうちで症状が軽減する時（以下、単にオン時という場合がある）と症状が悪化する時（以下、単にオフ時という場合がある）がある。オフ時には突進歩行またはすくみ足が多くみられる。【０００３】
突進歩行またはすくみ足などの異常歩行時には左右の足運びのバランスが崩れることが経験的に知られている。このため、左右の足運び（歩行）の同調性をリアルタイムでモニタできれば、異常歩行の兆候をとらえて転倒防止等の適切な処置を講じることができる。また、片脚に障害を持つ患者の歩行では歩行のステップ間隔に大きな左右差がある。このためステップ間隔の左右差を定量化できれば、リハビリまたは治療の効果の適切な理解につながる。ここで、同調性とは、歩行中において、右足着地から左足着地までの時間と、左足着地から右足着地までの時間との比が一定値（例えば１）を保っているかどうかを指す指標である。すなわち、右足着地から左足着地までの時間と、左足着地から右足着地までの時間との比が一定値（１）を保っていれば同調性は高いと判断され、右足着地から左足着地までの時間と、左足着地から右足着地までの時間との比が一定値（１）を保っていなければ同調性は低いと判断される。【０００４】
パーキンソン病患者の左右の足運び（ステップ）の同調性の評価については、複数のフットスイッチ型センサを用いて、パーキンソン病患者の左右の足運びの同調性の評価することが知られている。
例えば、歩行を表すリズムが１サイクル変化するのに要する時間と、歩行を表すリズムが次の１サイクル変化するのに要する時間との比を求めることで、歩行の同調性を判断する。従来、この比は、踵離地または爪先離地の特徴点の間隔から求められていた。
【先行技術文献】
【特許文献】
【０００５】
【非特許文献１】Ronny Bartsch et al. “Fluctuation and synchronization of gait intervals and gait force profiles distinguish stages of Parkinson's disease” Physica A: Statistical Mechanics and its Applications 383(2):455-465 Sep 2007.
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来におけるパーキンソン病患者の左右の足運びの同調性の評価では、蓄積されたデータについて統計的な計算処理を行っているため、連続的な処理またはリアルタイム処理には不適である。
また、従来におけるパーキンソン病患者の左右の足運びの同調性の評価では、特徴点の間隔を用いて同調性の評価を行なっている。このため、特徴点として波形のピークを用いた場合において、ピークが分裂するなど明瞭なピークを持たない波形から上記比を算出すると大きな誤差が生じる。すなわち、特徴点の間隔を用いた同調性の評価は、明瞭なピークを持たない波形には不適である。また、特徴点の間隔を用いた場合、波形のピーク位置に対応する時間毎にしか比を計算できない。すなわち、特徴点の間隔を用いた同調性の評価は、連続的な処理またはリアルタイム処理には不適である。【０００７】
本発明は、このような課題に鑑み創案されたもので、体動リズムの同調性を連続的に評価することを目的とする。
なお、前記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本発明の他の目的の1つとして位置付けることができる。
【課題を解決するための手段】
【０００８】
本発明者らの鋭意研究により、信号の位相に着目することで、信号の測定時間間隔を下限とする任意の時間間隔で、連続的に運動リズムの同調性を評価できることを見出した。
すなわち、本発明の要旨は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をそなえた情報処理方法に存する。【０００９】
また、本発明の要旨は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、をそなえた情報処理装置に存する。【００１０】
さらに、本発明の要旨は、所定の時刻における生体の繰り返し随意運動を表す信号の位相から前記信号の位相が所定角度変化した時刻、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻に基づく、前記繰り返し随意運動の同調性を評価するための評価係数を、前記評価係数を識別可能な形態にて出力する出力部をそなえた出力装置に存する。
また、本発明の要旨は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出部と、前記時刻算出部により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出部と、前記評価係数算出部により算出された評価係数を、前記評価係数を識別可能な形態にて出力する出力部と、をそなえた情報処理システムに存する。【００１１】
さらに、本発明の要旨は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をコンピュータに実行させる情報処理用プログラムに存する。【００１２】
また、本発明の要旨は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度変化した時刻（以下、第１時刻という）、および、所定の時刻における前記信号の位相から前記信号の位相が前記所定角度の２倍変化した時刻（以下、第２時刻という）を算出する時刻算出過程と、前記時刻算出過程により算出された、前記第１時刻および前記第２時刻に基づいて、前記繰り返し随意運動の同調性を評価するための評価係数を算出する評価係数算出過程と、をコンピュータに実行させる情報処理用プログラムを記録したコンピュータ読取可能な記録媒体に存する。
【発明の効果】
【００１３】
運動リズムの同調性を連続的かつリアルタイムで評価することが可能となる。
【図面の簡単な説明】
【００１４】
【図１】図１は、実施形態の一例としてのシステムの構成を示す図である。
【図２】図２は、実施形態の一例としての運動リズムを示す図である。
【図３】図３は、実施形態の一例としての判定部の動作を説明するための図である。
【図４】図４は、実施形態の一例としてのシステムの動作を説明するためのフローチャートである。
【図５】図５（Ａ）及び図５（Ｂ）は、実施形態の一例としての体動信号および運動リズムを示す図であり、図５（Ａ）は、健常な被験者の腹部中央に体動信号検出装置を装着した状態で、体動信号検出装置が１００Hzサンプリングにて測定した１５分間のウォーキング中の加速度信号の一部であり、図５（Ｂ）は、体動信号から抽出された運動リズムの一例を示す図である。
【図６】図６は、実施形態の一例としての時定数決定部の動作を説明するためのフローチャートである。
【図７】図７は、実施形態の一例としての時定数の変化に対する分散の変化を示す図である。
【図８】図８は、実施形態の一例としてのリズム抽出部の動作を説明するためのフローチャートである。
【図９】図９は、実施形態の一例としての評価係数決定部の動作を説明するためのフローチャートである。
【図１０】図１０（Ａ）〜図１０（Ｃ）は、実施形態の一例としてのパターンマッチング法を説明するための図であり、図１０（Ａ）は、図５（Ｂ）のデータから一部の波形を抜き出したものであり、図１０（Ｂ）は、図１０（Ａ）中＊印で示した時間を中心として、幅０．４秒の基準波を選び自己相関係数を計算した結果であり、図１０（Ｃ）は、図１０（Ａ）中○印で示した時間を中心として、幅０．４秒の基準波を選び自己相関係数を計算した結果である。
【図１１】図１１は、実施形態の一例としての評価係数を示す図である。
【図１２】図１２は、実施形態の一例としての判定部の動作を説明するためのフローチャートである。
【図１３】図１３は、実施形態の一例としての判定部の動作を説明するためのフローチャートである。
【図１４】図１４は、実施形態の一例としての評価係数を示す図である。
【図１５】図１５は、実施形態の一例としての運動リズムを示す図である。
【図１６】図１６（Ａ）〜図１６（Ｃ）は、実施形態の一例としてのシステムにより得られる評価係数を説明するための図であり、図１６（Ａ）は、体動信号検出装置１０によって検出された加速度を示す図であり、図１６（Ｂ）は、２階積分後の加速度信号すなわち運動軌道を示す図であり、図１６（Ｃ）は、２階積分後の加速度信号から求められた評価係数を示す図である。
【図１７】図１７（Ａ）〜図１７（Ｃ）は、実施形態の一例としてのシステムにより得られる評価係数を説明するための図であり、図１７（Ａ）は、体動信号検出装置１０によって検出された加速度を示す図であり、図１７（Ｂ）は、２階積分後の加速度信号すなわち運動軌道を示す図であり、図１７（Ｃ）は、２階積分後の加速度信号から求められた評価係数を示す図である。
【図１８】図１８（Ａ）〜図１８（Ｃ）は、実施形態の一例としてのシステムにより得られる評価係数を説明するための図であり、図１８（Ａ）は、体動信号検出装置１０によって検出された加速度および角速度を示す図であり、図１８（Ｂ）は、２階積分後の加速度信号および角速度信号を示す図であり、図１８（Ｃ）は、２階積分後の加速度信号および角速度からそれぞれ求められた評価係数を示す図である。
【図１９】図１９（Ａ）は、平滑化処理を施した後の歩行指数の時間変化を示す図であり、図１９（Ｂ）は、平滑化処理を施した後の体動指数の時間変化を示す図である。
【発明を実施するための形態】
【００１５】
以下、図面を参照して本情報処理方法に係る実施の形態を説明する。
〔Ａ〕実施形態の説明
図１は、実施形態の一例にかかるシステムの構成を示す図である。システム１は、体動信号検出装置１０と情報処理装置２０とをそなえる。体動信号検出装置１０は、情報処理装置２０と、例えば、インターネット等の有線または無線LAN（Local Area Network）またはBluetooth（登録商標）等の無線を介して通信可能に接続されている。
【００１６】
体動信号検出装置１０は、例えば、被験者（生体）の繰り返し随意運動（以下、単に随意運動という場合がある）に伴う繰り返しリズム運動を非侵襲的かつ連続的に検出（測定）するものである。以下、随意運動に伴う繰り返しリズム運動を、単にリズム運動という場合がある。
ここで、随意運動とは、例えば、歩行、ジョギング、ランニング、自転車、水泳、体操、ウェートトレーニング、体力測定（踏み台昇降、反復横とび）、ジャグリング（お手玉、サッカーボールのリフティング）等である。また、随意運動に伴う繰り返しリズム運動には、例えば、随意運動が歩行である場合には、歩行自体のリズム運動が含まれる。
【００１７】
また、随意運動が規則的に繰り返される運動である場合に、より正確に随意運動自体のリズムを抽出することができる。なお、規則的に繰り返される随意運動とは、完全に同じ運動が繰り返されることのみではなく、略同じ運動が繰り返されることを含む。
また、非侵襲的とは、例えば、被験者の体に傷をつけないこと、または、被験者に対して負担を与えないことを意味する。
【００１８】
体動信号検出装置１０は、例えば、携帯可能に構成される。なお、体動信号検出装置１０の被験者への取り付け位置は、体の動きを検知できる部位であれば、特に制限はない。
体動信号検出装置１０は、例えば、体動信号検出部１１，記憶部１２およびインターフェース部１３をそなえる。体動信号検出部１１，記憶部１２およびインターフェース部１３は、相互に通信可能に接続されている。
【００１９】
体動信号検出部１１は、例えば、随意運動に伴う繰り返しリズム運動を体動信号（生体の随意運動に伴う繰り返しリズム運動に基づく信号）として検出（測定）する。すなわち、体動信号検出部１１は、生体の随意運動に伴う繰り返しリズム運動に基づく信号を体動信号として検出する。異なる観点から見れば、一つの体動信号検出部１１は、例えば、生体の随意運動に伴う繰り返しリズム運動を体動信号として検出する。具体的には、被験者のリズム運動による、例えば、力の変化、空間的な身体の位置の変化、身体から発する音、電磁波等の波または微細エネルギーの変化または身体の周りにおける場の変化等を体動信号として検出（測定）する。すなわち、体動信号検出部１１は、随意運動に伴う繰り返しリズム運動に基づく信号を測定する。
【００２０】
ここで、体動信号検出部１１は、例えば、加速度センサ，速度センサ，ジャイロセンサ等の慣性センサにより実現される。上記体動信号を検出する慣性センサについては、例えば、検出する信号の種類に応じて適宜選択される。通常、歩行リズムを検出する場合には、体の動きの加速度を測定する加速度センサが好ましく用いられるが、加速度センサに限定されるものではない。
【００２１】
また、加速度センサとしては、一軸〜三軸のものを任意に用いることができるが、歩行時における鉛直方向、水平前後方向、及び水平左右方向の三方向へ作用する加速度を検出するための三軸加速度センサを用いることが好ましいが、三軸加速度センサに限定されるものではない。
なお、体動信号検出部１１は、例えば、所定のサンプリング周波数（例えば、１００Hz）で体動信号を測定する。
【００２２】
記憶部１２は、例えば、RAM（Random Access Memory），HDD（Hard Disk Drive），SSD（Solid State Drive），フラッシュメモリ等の各種情報を記憶可能な記憶装置である。記憶部１２は、具体的には、例えば、体動信号検出部１１によって得られた体動信号を記憶する。また、記憶部１２は、例えば、体動信号検出装置１０に対して着脱自在に設けられてもよいし、体動信号検出装置１０に固定されてもよい。
【００２３】
インターフェース部１３は、例えば、情報処理装置２０と通信可能に接続されるインターフェースである。体動信号検出装置１０が情報処理装置２０と無線を介して接続されている場合には、例えば、インターフェース部１３はアンテナを含むものである。また、体動信号検出装置１０が情報処理装置２０と有線を介して接続されている場合には、例えば、インターフェース部１３は、有線に接続可能な接続端子である。なお、記憶部１２が、例えば、体動信号検出装置１０に設けられたスロットに対して着脱自在に設けられており、記憶部１２を体動信号検出装置１０から取り外して情報処理装置２０に接続する場合には、インターフェース部１３は、体動信号検出装置１０に設けられなくてもよい。
情報処理装置２０は、例えば、PC（Personal Computer）であり、体動信号検出装置１０により得られた体動信号から後述する評価係数を算出する。
【００２４】
情報処理装置２０は、例えば、中央演算部２１，記憶部２２，出力部２３およびインターフェース部２４をそなえる。ここで、中央演算部２１，記憶部２２，出力部２３およびインターフェース部２４は、相互に通信可能に接続されている。
記憶部２２は、例えば、RAM，HDD，SSD等の、アプリケーションプログラム及びデータを格納可能な記憶装置であり、各種情報を記憶する。
【００２５】
中央演算部２１は、例えば、記憶部２２に記憶された各種アプリケーションプログラムを実行することにより種々の演算または制御を行ない、これにより、各種機能を実現する処理装置である。
例えば、中央演算部２１は、記憶部２２に記憶された情報処理用プログラムを実行することにより、時定数決定部２１１，リズム抽出部２１２，評価係数決定部２１３，判定部２１４および出力制御部２１５として機能する。すなわち、情報処理用プログラムは、中央演算部２１を、時定数決定部２１１，リズム抽出部２１２，評価係数決定部２１３，判定部２１４および出力制御部２１５として機能させるプログラムである。
【００２６】
時定数決定部２１１は、例えば、体動信号検出部１１により得られた体動信号を記憶部１２から取得し、取得した体動信号に基づいて、随意運動自体のリズムを表す信号（以下、単に運動リズムという場合がある）を抽出するために用いる時定数を算出する。すなわち、運動リズムは、生体の繰り返し随意運動を表す信号に相当する。また、運動リズムが含まれる体動信号は一つの体動信号検出部１１により検出された信号であるため、生体の繰り返し随意運動を表す信号は一つの体動信号検出部１１（検出部）により検出された信号である。
【００２７】
具体的には、時定数決定部２１１は、以下の処理を行なうことで、時定数を決定する。
（１）時定数決定部２１１は、体動信号Xを、ある時定数Aで特徴付けられるハイパスフィルタ、あるいはバンドパスフィルタにかける。ここで、例えば、ハイパスフィルタとしては、体動信号Xを時間幅（時定数）Aのゼロ位相移動平均フィルタで平滑化する処理をF(X,A)と記載するとき、出力波形（信号）Y=X-F(X,A)で表現される処理を行なう。また、例えば、バンドパスフィルタとしては、Y=F(X-F(X,A),A/2.5)で表現される処理を行なう。ここで、ゼロ位相移動平均フィルタとは、位相ずれが０である移動平均フィルタを指す。なお、ゼロ位相移動平均フィルタは既知の種々の手法を用いて実現可能であり、詳細な説明は省略する。なお、バンドパスフィルタの時定数の一例として、A/2.5としているが、これに限定されるものではない。
【００２８】
（２）次に、時定数決定部２１１は、出力波形Yの規則性を数値化する。時定数決定部２１１は、例えば、出力波形Yの極大及び極小のピークにおける値の絶対値のCV（＝標準偏差／平均値：Coefficient of Variation）を求める。より具体的には、時定数決定部２１１は、出力波形Yから、極大および極小のピークを抽出し、ピーク位置における波形の振幅の絶対値を求め、求められた波形の振幅から標準偏差および平均値を求める。そして、時定数決定部２１１は、標準偏差を平均値で除することで、分散すなわちCVを算出する。
【００２９】
（３）さらに、時定数決定部２１１は、時定数Aを変化させた場合の出力波形Yの規則性を求める。時定数決定部２１１は、例えば、時定数Tを変化させた場合の出力波形Yの規則性をグラフ化する。具体的には、時定数決定部２１１は、例えば、時定数Aを変化させた場合の、各時定数AにおけるCVを算出することで、時定数Aに対する出力波形Yの規則性の変化を求める。
なお、時定数決定部２１１における出力波形Yの規則性の数値化として、出力波形Yのピークにおける絶対値のCVを利用する例を述べたが、この例に限定されるものではない。
（４）時定数決定部２１１は、上記（３）の処理で求められた時定数Aに対する出力波形Yの規則性の変化（CVの変化）から、極小点を求める。そして、時定数決定部２１１は、例えば、求められた極小点（例えば２つの極小点）のうち時定数Aの小さい方の極小点に対応する時定数Aを、運動リズムを体動信号から分離して抽出するための時定数（以下、第１時定数という場合がある）として決定する。すなわち、時定数決定部２１１は、なるべく規則性のあるような波形を分離して抽出することが可能となる時定数を選択する。言い換えれば、時定数決定部２１１は、複数の極小点（例えば２点）を選択し時定数Aの小さい方の極小点に対応する時定数を第１時定数として決定する。すなわち、時定数決定部２１１は、体動信号に対してバンドパスフィルタを適用した場合の出力に基づいて第１時定数を決定する。
【００３０】
なお、上記（１）〜（４）の過程で、体動信号をN階積分（N=1,2,3.…）する過程を含めた方がノイズ除去の観点から好ましい。N階積分する場合、上記（１）の処理（フィルタ処理）をN回以上施すことが好ましい。積分操作とフィルタ処理の順番は、例えば、以下の条件の下であればどのような順番でも良い。
・フィルタ処理は２回以上連続して行わない。
・フィルタ処理のうち１回は積分処理がすべて終了してから行う。
・フィルタ処理を２回以上行う場合、最後の１回を除いて、そのフィルタはハイパスフィルタとする。
【００３１】
リズム抽出部２１２は、例えば、体動信号検出部１１により得られた体動信号を記憶部１２から取得し、時定数決定部２１１により決定された第１時定数（A１）を用いて、フィルタリング処理を行なうことで体動信号から運動リズムを抽出する。
具体的には、リズム抽出部２１２は、例えば、以下の処理（フィルタリング処理）を行なうことで体動信号から運動リズムを抽出する。なお、下記の処理は、例示であり、この例に限定されるものではない。
【００３２】
（５）例えば、リズム抽出部２１２は、体動信号に対して、第１時定数のゼロ位相移動平均フィルタを施すことで、体動信号の低周波成分を求める。すなわち、リズム抽出部２１２は、体動信号に対して、F(X,A1)で表される処理を行なう。次に、リズム抽出部２１２は、体動信号から、F(X,A1)で表される処理を行なうことで求められた低周波成分を引く。すなわち、リズム抽出部２１２は、Y=X-F(X,A1)を求める。言い換えれば、リズム抽出部２１２は、第１時定数を用いて、体動信号に対してハイパスフィルタリングを行なう。
【００３３】
なお、このハイパスフィルタリングによって、例えば、後述する積分処理において発散防止が可能となる。
（６）次に、リズム抽出部２１２は、ハイパスフィルタリング後の信号に対して積分を行なう。積分の階数は任意であり、例えば、リズム抽出部２１２は、ハイパスフィルタリング後の信号に対して２階積分を行なう。
【００３４】
（７）リズム抽出部２１２は、処理（１１）における積分後の信号X2に対して、（１０）と同様の処理を行なう。すなわち、リズム抽出部２１２は、積分後の体動信号に対して、F(X2,A1)で表される処理を行なう。言い換えれば、リズム抽出部２１２は、第１時定数を用いて、積分後の信号である積分信号に対してハイパスフィルタリングを行なう。
（８）さらに、リズム抽出部２１２は、上記（１２）の処理で得られた信号の極大値を結ぶ包絡線と、極小値を結ぶ包絡線とを作成し、これらの２本の包絡線の平均からなる信号を、上記（１２）の処理で得られた信号から引く処理を行なう。この処理により得られた信号の振幅は、振幅の値０を挟んで変化することとなるため、後述するヒルベルト変換法を用いる場合において、位相を正確に求めることが可能となる。なお、この処理は省略してもよい。
【００３５】
上記（５）〜（８）の過程により、リズム抽出部２１２は、体動信号から運動リズムを抽出する。すなわち、リズム抽出部２１２は、体動信号にフィルタリング処理を施すことにより、体動信号から、随意運動のリズム（運動リズム）を抽出するリズム抽出部として機能する。従って、体動信号から、例えば歩行リズムの同調性を評価する場合には、あらかじめ歩行領域を抽出してもよいし、そのような前処理を施さなくてもよい。
なお、上記（６）の過程で、体動信号を２階積分しているが、積分処理の実施に限定されるものではなく、積分処理を行なわなくてもよい。また、２階積分に限定されるものではなく、N階積分を行なってもよい（N=1,2,3,…）。なお、N階積分する場合、上記（５）［または（７）］の処理（フィルタ処理）をN回以上施すことが好ましい。積分操作とフィルタ処理の順番は、例えば、以下の条件の下であればどのような順番でも良い。
・フィルタ処理は２回以上連続して行わない。
・フィルタ処理のうち１回は積分処理がすべて終了してから行う。
・フィルタ処理を２回以上行う場合、最後の１回を除いて、そのフィルタはハイパスフィルタとする。
【００３６】
評価係数決定部２１３は、例えば、運動リズムの位相に基づいて、運動リズムの同調性を評価するための評価係数を決定する。ここで、位相とは、時間的に繰り返されるリズムを、円周上を回転する運動と見立て、リズム波形上の各点が何度の回転角度に相当するかを表す指標である。例えば、正弦波において隣り合う二つのピーク点の位相には３６０度の差がある。図２は、運動リズムの一例を示す図であり、図中のピークは、例えば、被験者の右足の着地、または、左足の着地に対応する。ここで、評価係数の一例は、例えば、図２に示す運動リズムの一例において、一方の足（例えば、右足）に対応するピーク間における一方の足に対応するピークと他方の足に対応するピークとの時間間隔をそれぞれ第１時間T１，第２時間T２として、下記（１）式により求められる。なお、運動リズムの同調性を評価するための評価係数は、この評価係数に限定されるものではない。例えば、下記（１）式では、第１時間T１と第２時間T２とが等しいときの評価係数の値は０となるが、０に限定されるものではない。例えば、第１時間T１と第２時間T２とが等しいときの評価係数の値を１等０以外の値となるような評価係数の算出式を用いてもよい。すなわち、評価係数の算出式は下記（１）式に限定されるものではない。なお、下記（１）式は、一方の足の着地から他方の足の着地までの時間と他方の足の着地から一方の足の着地までの時間との差の半分と、一方の足の着地から他方の足の着地までの時間と他方の足の着地から一方の足の着地までの時間との和との比を表している。
【００３７】
【数１】
2012036135.tif 000003
【００３８】
なお、図２においては、一方の足（例えば、右足）に対応するピーク間における一方の足に対応するピークと他方の足に対応するピークとの時間間隔をそれぞれ第１時間T１，第２時間T２としたが、ピーク間の時間間隔に限定されるものではない。例えば、第１時間T１は、運動リズムのピークに相当しない所定の位置から位相が３６０度変化するまでの時間であってもよい。さらに、例えば、第２時間T２は、運動リズムのピークに相当しない所定の位置から位相が３６０度変化した位置から、さらに位相が３６０度変化するまでの時間であってもよい。
【００３９】
なお、第１時間T１は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が３６０度変化するまでの時間であり、第２時間T２は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が３６０度変化した位置から、さらに位相が３６０度変化するまでの時間であるとしたが、これに限定されない。例えば、第１時間T１は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が−３６０度変化するまでの時間（所定の位置から３６０度遡った時間）であり、第２時間T２は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が−３６０度変化した位置から、さらに位相が−３６０度変化するまでの時間（所定の位置から７２０度遡った時間）であってもよい。また、例えば、第１時間T１は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が３６０度変化するまでの時間であり、第２時間T２は、運動リズムのピーク、またはピークに相当しない所定の位置から位相が−３６０度変化するまでの時間（所定の位置から３６０度遡った時間）であってもよい。
【００４０】
ところで、第１時間T１を所定の位置から３６０度遡った時間として設定する場合、または、第２時間T２を所定の位置から３６０度または７２０度遡った時間として設定する場合、第１時間T１及び第２時間T２の値が負となる場合がある。このような場合には、数１に代入する第１時間T１及び第２時間T２に留意する必要がある。
【００４１】
評価係数決定部２１３は、時刻算出部２２３および評価係数算出部２３３をそなえる。
時刻算出部２２３は、例えば、運動リズムの位相を算出することで、所定の時刻（任意の時刻）における運動リズムの位相から位相が所定角度（以下、第１角度という場合がある）変化した時刻（第１時刻）および、所定の時刻における運動リズムの位相から位相が第１角度の２倍変化した時刻（第２時刻）を算出する。
【００４２】
また、時刻算出部２２３は、例えば、所定の時刻から第１時刻までの時間（第１時間： T１）を算出することができる。また、時刻算出部２２３は、例えば、第１時刻から第２時刻までの時間（第２時間：T２）を算出することができる。
ここで、第１角度とは、例えば、１８０度または３６０度の正の値である。なお、本願において、３６０度とは、厳密に３６０度である場合および略３６０度である場合を含み、同様に、１８０度とは、厳密に１８０度である場合および略１８０度である場合を含む。
【００４３】
なお、時刻算出部２２３は、例えば、運動リズムの位相を算出することで、所定の時刻（任意の時刻）における運動リズムの位相から位相が第１角度遡った時刻（所定角度遡った時刻も、第１時刻に含めることとする。）および、所定の時刻における運動リズムの位相から位相が第１角度の２倍変化した時刻（第１角度の倍遡った時刻も、第２時刻に含めることとする。）を算出するようにしてもよい。つまり第１角度は、例えば、−１８０度または−３６０度の負の値であっても構わない。この場合、時刻算出部２２３は、第１時刻から所定の時刻までの時間（第１時間： T１）を算出することができる。また、時刻算出部２２３は、第２時刻から所定の時刻までの時間（第２時間：T２）を算出することができる。
【００４４】
さらに、時刻算出部２２３は、所定の時刻（任意の時刻）における運動リズムの位相から位相が第１角度変化した時刻（第１時刻）および、所定の時刻における運動リズムの位相から位相が第１角度の−１倍変化した時刻を算出するようにしてもよい。この場合、時刻算出部２２３は、所定の時刻から第１時刻までの時間（第１時間： T１）を算出することができる。また、時刻算出部２２３は、第２時刻から所定の時刻までの時間（第２時間：T２）を算出することができる。
以下では、評価係数算出部２３３が第１時刻（第１時刻は所定の時刻よりも大きい）及び第２時刻（第２時刻は第１時刻よりも大きい）を算出する場合について説明することとする。第１時刻（第１時刻は所定の時刻よりも小さい）及び第２時刻（第２時刻は第１時刻よりも小さい）を算出する場合、及び第１時刻（第１時刻は所定の時刻よりも大きい）及び第２時刻（第２時刻は所定の時刻よりも小さい）を算出する場合についても、その演算手法は同様である。
【００４５】
時刻算出部２２３は、例えば、ヒルベルト変換法またはパターンマッチング法を用いて、運動リズムの位相を算出する。なお、位相を算出する方法は、ヒルベルト変換法またはパターンマッチング法に限定されるものではない。
ここで、ヒルベルト変換法は、波形の所定の位置における位相を具体的に算出する手法であり、パターンマッチング法は直接的に位相を特定しないが、互いに位相が３６０度ずれた点を見つける手法である。すなわち、パターンマッチング法は、波形の所定の位置から位相が３６０度の整数倍ずれた点を見つける手法である。言い換えれば、パターンマッチング法は、波形の所定の位置に対する相対的な位相（例えば３６０度）を算出する。
【００４６】
評価係数算出部２３３は、例えば、時刻算出部２２３によって算出された、第１時刻および第２時刻に基づいて、評価係数を算出する。例えば、評価係数算出部２３３は、時刻算出部２２３によって算出された第１時刻および第２時刻に基づいて第１時間および第２時間を算出し、第１時間および第２時間を上記（１）式に代入することで、評価係数を算出する。すなわち、評価係数算出部２３３は、第１時刻および第２時刻に基づいて、繰り返し随意運動の同調性を評価するための評価係数を算出する。
【００４７】
まず、時刻算出部２２３が、ヒルベルト変換法を用いて位相変化量の比を算出する場合について説明する。
（Ａ）ヒルベルト変換法を用いる場合
ヒルベルト変換は、任意の実数時系列信号X(t)から対応する虚数部分Y(t)を導く数学的手法である。これにより、時刻算出部２２３は、下記（２）式からX(t)の位相θ(t)を直接的に求めることができる。すなわち、位相算出部２２４は、下記（２）式から所定時刻（所定位置）における運動リズムの位相を求めることができる。
【００４８】
【数２】
2012036135.tif 000004
【００４９】
また、任意の時刻tにおける運動リズムの位相をθG(t)度とすると、時刻算出部２２３は、上記（２）式から、運動リズムの位相を求めることで、例えば、θG (t1)=θG (t)+360となるような時刻t１を求める。さらに、時刻算出部２２３は、例えば、θG (t2)=θG (t)+720となるような時刻t２を求める。
【００５０】
そして、評価係数算出部２３３は、例えば、時刻算出部２２３によって算出された時刻t１および時刻t２を用いて、例えば、下記（３）式から評価係数を算出する。
【００５１】
【数３】
2012036135.tif 000005
【００５２】
なお、（３）式は、（１）式と等価な式である。なお、評価係数算出部２３３は、時刻t１および時刻t２に基づいて算出された第１時間T１および第２時間T２を用いて、上記（１）式から評価係数を算出してもよい。
次に、時刻算出部２２３が、パターンマッチング法を用いて位相変化量の比を算出する場合について説明する。
（Ｂ）パターンマッチング法を用いる場合
パターンマッチング法は二つの信号の類似性を定量化する手法である。類似度の定義または計算法は多数あるが、具体的には例えば、「画像処理工学（末松良一・山田宏尚著、コロナ社）」等に記載の方法が用いられる。最も代表的なのが自己相関係数である。例えば3次元の体動信号については以下のような計算を行う。
【００５３】
まず体動データから適当な基準波を選び出す。基準波の座標を下記（４）式とする。
【００５４】
【数４】
2012036135.tif 000006
【００５５】
ここでp個のx、y、zはそれぞれの平均値がゼロになっているとする。下記（５）式で自己相関係数が計算される。
【００５６】
【数５】
2012036135.tif 000007
【００５７】
これはいわゆるスカラー量であり、座標系のとり方に依存しない。すなわち、体動測定中に体動信号検出装置１０が装着部位で回転のずれを起こしても同じ値となる。なお、1次元の信号についても同様にして計算する。
上記パターンマッチング法を用いると、任意の時刻tにおけるリズム波形の位相か３６０度だけ位相がずれた時刻t１または７２０度だけ位相がずれた時刻t２等を容易に求めることができる。言い換えれば、時刻算出部２２３は、パターンマッチングを用いることで、時刻tにおける運動リズムの位相から、それぞれ３６０度の整数倍だけ位相がずれた時刻を求めることができる。
【００５８】
そして、評価係数算出部２３３は、例えば、時刻算出部２２３によって算出された時刻t１および時刻t２を用いて、例えば、上記（３）式から評価係数を算出する。なお、評価係数算出部２３３は、時刻t１および時刻t２に基づいて算出された第１時間T１および第２時間T２を用いて上記（１）式から評価係数を算出してもよい。
【００５９】
判定部２１４は、例えば、評価係数決定部２１３により求められた評価係数に基づいて、運動リズムの同調性の評価（判定）を行なう。すなわち、判定部２１４は、評価係数に基づいて、運動リズムの同期またはバランスを評価する。具体的には、例えば、随意運動が歩行等である場合には、判定部２１４は、左右の足運びの同調性、すなわち左右の足運びの同期またはバランスの評価を行なう。ここで、左右の足運びの同期の良否は、評価係数が時間によらず一定値を保っているか、または、評価係数が周期的に変化しているか否かにより判断する。例えば、判定部２１４は、評価係数が時間によらず一定値を保っている場合または評価係数が周期的に変化している場合、には、左右の足運びの同期は良いと判定する。また、左右の足運びのバランス良否は、上記（１）式を用いて評価係数を決定する場合、評価係数が０に近いか否かにより判断する。例えば、判定部２１４は、評価係数が０に近いほど左右の足運びのバランスが良いと判定する。
【００６０】
具体的には、判定部２１４は、所定時間にわたる評価係数の変化から、例えば、以下のR，SおよびCの３つの指標を算出する。所定時間とは、例えば、随意運動が歩行である場合には、被験者が１０歩程度の歩行に要する時間（例えば５秒）であるが、この時間に限定されるものではない。
（i）評価係数の範囲（０から最も離れた値）：０±R
（ii）評価係数の標準偏差（ゆらぎの大きさ）：S
（iii）平均歩行周期（図２または後述する図１５におけるT１＋T２）だけ時間をずらしたときの評価係数の自己相関係数（評価係数の周期性）：C
判定部２１４は、これらのR，SおよびCの３つの指標に基づいて、同調性の評価を行なう。図３は、判定部２１４の動作の一例を説明するための図である。
【００６１】
図３に示すように、判定部２１４は、例えば、Rが０．０２以下であれば、左右の足運び（歩行リズム）の同期および左右の足運びのバランスは非常に良い（図３中◎参照）と判断する。すなわち、判定部２１４は、評価係数が、所定時間（例えば５秒）にわたって、所定範囲（例えば０±０．０２）内であるかを判定する判定部として機能する。
また、判定部２１４は、例えば、Rが０．０２以下ではないが、Sが、０．０１以下である場合には、左右の足運びの同期および左右の足運びのバランスは良い（図３中○参照）と判断する。
【００６２】
さらに、判定部２１４は、例えば、Rが０．０２以下ではなく、かつ、Sが０．０１以下ではない、場合には、Cに基づいて、左右の足運びの同期および左右の足運びのバランスを評価する。
判定部２１４は、例えば、Cが０．５以上の場合には、左右の足運びの同期は良いが、左右の足運びのバランスは悪い（図３中×参照）と判断し、Cが０．５より小さく０．２以上の場合には、左右の足運びの同期は普通（図３中△参照）だが、左右の足運びのバランスは悪い（図３中×参照）と判断する。さらに、判定部２１４は、例えば、Cが０．２未満の場合には、左右の足運びの同期は悪く（図３中×参照）、左右の足運びのバランスは判定不能であると判定する。
【００６３】
図３における数値（０．０２，０．０１，０．５および０．２）は例示であり、これらの数値に限定されるものではない。また、図３では、R，SおよびCの３つの指標を用いて５パターンの判定を行なっているが、この判定手法に限定されるものではなく、より細かい判定を行なってもよいし、より大まかな評価を行なうこととしてもよい。例えば、Rの値のみに基づいて２パターンの判定を行なうこととしてもよい。
以上の他にも、判定部２１４では、所定時間にわたる評価係数の変化から、左右の足運びの同期および左右の足運びのバランスを、例えば以下のように直接指標化することもできる。所定時間とは、例えば、随意運動が歩行である場合には、被験者が１０歩程度の歩行に要する時間（例えば５秒）であるが、この時間に限定されるものではない。
【００６４】
ここで、運動リズムの位相から算出された評価係数を時刻tの関数として、H(t)と記す。任意の時刻t1を起点として運動リズムの位相が720度変化したときの時刻をt2とする。
（ｉ）左右の足運びのバランス：S1
時刻t1≦t≦時刻t2の間におけるH(t)の標準偏差を求める。所定時間にわたり時刻t1を逐次変化させて標準偏差の時間変化を算出する。こうして得られた標準偏差の時系列データについて、所定時間内における平均値を求め、左右の足運びのバランスの指標S1とする。評価係数が（３）式で計算される場合、左右の足運びのバランスが完璧であればこの指標は０となり、バランスが崩れてくると０．５に近づく。歩行に問題のない成人であれば、この値は通常0.02以下となる。
（ｉｉ）左右の足運びの同期：S2
H(t2)‐H(t1)の絶対値を求める。所定時間にわたり時刻t1を逐次変化させて絶対値の時間変化を算出する。次にこうして得られた絶対値の時系列データの包絡線を求める。具体的には例えば、ヒルベルト変換の式を用いる。すなわち、（２）式において絶対値をX(t)とした場合のA(t)が包絡線に相当する。所定時間内における包絡線の平均値を求め、左右の足運びの同期の指標S2とする。評価係数が（３）式で計算される場合、左右の足運びの同期が完璧であればこの指標は０となり、同期が崩れてくると１に近づく。歩行に問題のない成人であれば、この値は通常0.02以下となる。
【００６５】
また、前述した時定数決定部２１１における出力波形Yの規則性の数値化としては、時刻算出部２２３で得られる位相や、判定部２１４で得られる指標を利用することもできる。例えば、位相の周期のCVあるいは左右の足運びの同期の指標S2を出力波形Yの規則性を反映しているとみなす。そして、フィルタの時定数Aを変化させて（例えば、歩行リズムの場合は、0から1秒の間で）、位相の周期のCVあるいは指標S2を求め、その値が最小値をとるような時定数Aを最適な時定数として決定する。
【００６６】
出力制御部２１５は、出力部２３を制御する。例えば、出力部２３がディスプレイの表示部をそなえる場合には、出力制御部２１５は、出力部２３の表示状態を制御することで、各種の情報を出力部２３に表示させる。例えば、出力制御部２１５は、リズム抽出部２１２によって抽出された運動リズムまたは評価係数決定部２１３によって求められ評価係数を出力部２３に表示させる。評価係数の表示は、グラフとして出力部２３に表示させてもよいし、「評価係数は＊＊です」のようにメッセージとして表示させてもよい。また、出力制御部２１５は、例えば、判定部２１４による判定結果を、評価係数の表示に加え、または、評価係数の表示に替えて、出力部２３に表示させてもよい。
【００６７】
また、出力部２３が、例えば、アラームまたは振動等により警告を行なうものであれば、判定部２１４により、例えば、評価係数が所定の範囲内の値ではないと判断された場合に、出力制御部２１５は、出力部２３を制御することで出力部２３に警告を行なわせる。
また、出力部２３が、情報処理装置２０外にそなえられている場合には、出力制御部２１５は、例えば、リズム抽出部２１２によって抽出された運動リズムまたは評価係数決定部２１３によって求められ評価係数等の各種情報を、無線または有線を介して出力部２３に送信することで、出力部２３の表示状態を制御する。
【００６８】
出力部２３は、出力制御部２１５の制御の下、各種の情報を出力する。例えば、出力部２３は、ディスプレイであり、評価係数決定部２１３によって求められた評価係数、または、判定部２１４により、例えば、評価係数が所定の範囲内の値ではないと判断された場合には、警告を表示する。
また、出力部２３は、例えばアラームまたは振動等によって、状態の変化または突然の異常について報知するもの等であってもよく、判定部２１４により、例えば、評価係数が所定の範囲内の値ではないと判断された場合には、アラームまたは振動等により警告を行なう。
【００６９】
なお、出力部２３は、情報処理装置２０にそなえられなくともよく、情報処理装置２０外にそなえられてもよい。ここで、出力部２３が、情報処理装置２０外にそなえられている場合とは、例えば、出力部２３が体動信号検出装置１０に含まれる場合、または、体動信号検出装置１０および情報処理装置２０にどちらにも含まれない場合である。情報処理装置２０外にそなえられている場合において、例えば、情報処理装置２０は、無線または有線を介して出力部２３と接続されている。すなわち、出力部２３は、情報処理装置により求められた結果（例えば、評価係数）を出力する出力部として機能する。また、出力部２３は、情報処理装置により求められた結果（例えば、評価係数）を出力する出力部をそなえた出力装置として機能する。
【００７０】
インターフェース部２４は、例えば、体動信号検出装置１０と通信可能に接続されるインターフェースである。情報処理装置２０が体動信号検出装置１０と無線を介して接続されている場合には、例えば、インターフェース部２４はアンテナを含むものである。また、情報処理装置２０が体動信号検出装置１０と有線を介して接続されている場合には、例えば、インターフェース部２４は、有線の接続端子である。なお、記憶部１２が、例えば、体動信号検出装置１０に対して着脱自在に設けられており、記憶部１２を体動信号検出装置１０から取り外して情報処理装置２０に接続する場合には、インターフェース部２４は、例えば、記憶部１２が接続されるスロットとしても機能する。
【００７１】
上述の如く構成された、実施形態の一例としてのシステム１の動作を、図４に示すフローチャート（ステップＡ１〜Ａ６）を参照しながら説明する。以下、積分信号をもとに運動リズムの同調性の評価を行う具体的な方法の一例を、歩行中の加速度信号を用いた場合について詳細に説明するが、他の慣性センサまたは他の運動中の信号の場合でも同様の処理が適用できる。
【００７２】
まず、体動信号検出部１１が、随意運動に伴う繰り返しリズム運動を体動信号として検出する（ステップＡ１）。図５（Ａ）は、健常な被験者の腹部中央に体動信号検出装置１０（３軸の加速度センサ）を装着した状態で、体動信号検出装置１０（体動信号検出部１１）が１００Hzサンプリングにて測定した１５分間のウォーキング中の加速度信号の一部である。なお、図５（Ａ）においては、３軸のうち、上下方向の加速度変化のみを示している。次に、体動信号検出部１１により検出された体動信号に基づいて、時定数決定部２１１は、運動リズムを抽出するための時定数を決定する（ステップＡ２）。そして、ステップＡ２にて決定された時定数を用いて、リズム抽出部２１２は、体動信号から運動リズムを抽出する（ステップＡ３）。図５（Ｂ）は、体動信号から抽出された運動リズムの一例を示す図である。
【００７３】
そして、評価係数決定部２１３は、抽出された運動リズムから評価係数を決定する（ステップＡ４）。判定部２１４は、この評価係数に基づいて、例えば、運動リズムの同調性の判定を行なう（ステップＡ５：判定過程）。次に、判定部２１４の判定結果に基づいて、出力制御部２１５は出力部２３に判定結果を出力させる（ステップＡ６）。
次に、時定数決定部２１１の詳細、すなわち図２におけるステップＡ２の詳細な動作を、図６示すフローチャート（ステップＡ２１〜Ａ２４）を参照しながら説明する。
【００７４】
まず、時定数決定部２１１は、所定の時定数Aを用いて、体動信号検出部１１により検出された体動信号に対してフィルタリングを行なう（ステップＡ２１）。例えば、時定数決定部２１１は、体動信号に対してY=X-F(X,A)で表現されるハイパスフィルタリングを行なった後に、２階積分を行ない、積分後の信号に再度Y=X-F(X,A)で表現されるハイパスフィルタリングを行なう。さらに、時定数決定部２１１は、例えば、ハイパスフィルタリング後に、Y=F(X,A/2.5)で表現されるローパスフィルタリングを行なう。すなわち、時定数決定部２１１は、体動信号に対して、バンドバスフィルタリングを行なう。
【００７５】
次に、時定数決定部２１１は、例えば、出力波形Yの極大及び極小のピークにおける値の絶対値のCVを求める（ステップＡ２２）。そして、時定数決定部２１１は、例えば、時定数Aを変化させ、上記ステップＡ２１，Ａ２２の処理を行なうことで、時定数Aの変化に対するCVの変化を求める（ステップＡ２３）。図７は、ステップＡ２３により求められた時定数（フィルタ時間）Aの変化に対するCV（分散）の変化の一例を示す図である。なお、この図７は、体動信号に呼吸リズムが含まれている場合の、時定数（フィルタ時間）Aの変化に対するCV（分散）の変化の一例を示す図である。体動信号に呼吸リズムが含まれていない場合には、図７における◆印で示される極小点は現れず、極小点は●印で示される極小点のみになる。次に、時定数決定部２１１は、ステップＡ２３により求められた時定数Tの変化に対するCVの変化から、極小点を求める。時定数決定部２１１は、例えば、図７に示すように、CVの値が小さい２つの極小点（図７中、●印および◆印）を求める。そして、時定数決定部２１１は、例えば、時定数Aの小さい方の極小点（図７中、●印）に対応する時定数A近傍の時定数の値（例えば０．４）を、第１時定数として決定する（ステップＡ２４）。このように、時定数決定部２１１は、極小点そのものに対応する時定数を第１時定数とするのではなく、極小点の近傍の時定数を第１時定数として決定してもよく、いずれの場合においても得られる結果に大差は生じない。
【００７６】
次に、リズム抽出部２１２の詳細、すなわち図４におけるステップＡ３の詳細な動作を、図８示すフローチャート（ステップＡ３７〜３９）を参照しながら説明する。
まず、リズム抽出部２１２は、体動信号検出部１１によって検出された体動信号に対して、第１時定数を用いてハイパスフィルタリングを行なう（ステップＡ３７）。リズム抽出部２１２は、N（例えば２）階積分を行なう（ステップＡ３８）。そして、積分後の信号に対して、再び、第１時定数を用いてハイパスフィルタリングを行なう（ステップＡ３９）。図５（Ｂ）は、ステップＡ３９にて得られる信号の一例を示す図である。この波形は、例えば、歩行である被験者の随意運動に相当する。より具体的には、ステップ３９にて得られた波形は、歩行による一歩一歩の相対運動軌道に相当する。すなわち、図５（Ｂ）に示された波形は、運動リズムを示す波形に相当する。
【００７７】
次に、評価係数決定部２１３の詳細、すなわち図４におけるステップＡ４の詳細な動作を、図９に示すフローチャート（ステップＡ４１〜Ａ４３）を参照しながら説明する。
まず時刻算出部２２３は、運動リズムの任意の時刻tにおける任意の点を選択する（ステップＡ４１）。時刻算出部２２３は、例えば、ヒルベルト変換またはパターンマッチング等を用いて位相を算出することで、任意の点から位相が３６０度ずれた点の時刻t１および、任意の点から位相が７２０度ずれた点の時刻t２を算出する（ステップＡ４２：時刻算出過程）。次に、評価係数算出部２３３は、時刻t１，t２を用いて、例えば上記（１）式から評価係数を算出する（ステップＡ４３：評価係数算出過程）。すなわち、ステップＡ４２は、生体の繰り返し随意運動を表す信号から、所定の時刻における前記信号の位相から前記信号の位相が所定角度（第１角度）変化した時刻（第１時刻）、および、所定の時刻における信号の位相から信号の位相が第１角度の２倍変化した時刻（第２時刻）を算出する時刻算出過程の一例である。また、ステップＡ４３は、第１時刻および第２時刻に基づいて、繰り返し随意運動の同調性を評価するための評価係数を算出する第２算出過程の一例である。
【００７８】
ここで、図１０（Ａ）〜図１０（Ｃ）は、パターンマッチング法を用いた場合に得られる信号の一例を説明するための図である。図１０（Ａ）〜図１０（Ｃ）を用いてパターンマッチングを用いた場合の評価係数決定部２１３の詳細な動作の一例を説明する。
図１０（Ａ）は、図５（Ｂ）のデータから一部の波形を抜き出したものである。ここで時刻算出部２２３が、＊印で示した時間を中心として、幅０．４秒の基準波を選び自己相関係数を計算した結果が図１０（Ｂ）である。同様に、図１０（Ａ）の○印を中心として、幅０．４秒の基準波から得られた自己相関係数が図１０（Ｃ）である。図１０（Ｂ）および図１０（Ｃ）において、位相が３６０度ずつずれた３点の位置を点線で示している。
【００７９】
評価係数算出部２３３が、例えば、図１０（Ｂ）または図１０（Ｃ）に示された波形から、第１時間T１および第２時間T２を算出し、時間T１，T２を上記（１）式に代入することで評価係数を算出する。
次に、図５（Ｂ）で例示した積分信号について、ヒルベルト変換法を用いて運動リズムの位相が３．６度増加するごとに評価係数を求めた例の結果を図１１の実線で示す。これは時間間隔に換算すると平均０．０５秒ごとに評価係数を求めていることに相当する。図１１から、評価係数が０近傍で周期的にゆらいでいることがわかる。これは、どのような生体でも体が完全に左右対称であることは稀であり、特に疾病または障害を持たなくても左右のステップ間隔に差が生じ、跛行気味に歩くからである。すなわち図２においてT１/T２≠1となっているからである。ただし、評価係数のゆらぎ自体は極めて規則的に変化しているので、左右の足運びの同期は高いといえる。
【００８０】
また、図５（Ｂ）で例示した積分信号について、パターンマッチング法を用いて時間間隔０．０５秒ごとに評価係数を求めた結果を図１１の破線で示す。図１１からヒルベルト変換法で求めた結果とほぼ同じ傾向の結果が得られているのがわかる。
次に、判定部２１４の詳細、すなわち図４におけるステップＡ５の詳細な動作を、図１２に示すフローチャート（ステップＡ５１〜Ａ５３）を参照しながら説明する。
【００８１】
まず、評価係数決定部２１３が、所定時間にわたって評価係数を算出する（ステップＡ５１）すなわち、所定時間にわたって上記ステップＡ４１〜Ａ４４が繰り返される。次に、判定部２１４は、所定時間にわたって算出された評価係数を用いて、R，SおよびCの３つの指標を算出する（ステップＡ５２）。そして、判定部２１４は、この３つの指標に基づいて運動リズムの同調性の判定を行なう（ステップＡ５３）。
【００８２】
さらに、判定部２１４による同調性の判定の詳細な処理、すなわち図１２におけるステップＡ５３の詳細な動作を、図１３に示すフローチャート（ステップＡ５３１〜Ａ５３９）を参照しながら説明する。
まず、判定部２１４は、Rが０．０２以下か否か判定する（ステップＡ５３１）。Rが０．０２以下の場合（ステップＡ５３１のＹｅｓルート参照）、所定時間にわたり、評価係数は略０であるとして、判定部２１４は、運動リズムの同期およびバランスは非常に良いと判定する（ステップＡ５３２）。すなわち、ステップＡ５３１は、評価係数が、所定時間にわたって、所定範囲（例えば、０±０．０２）内であるかを判定する判定過程の一例である。一方、Rが０．０２より大きい場合（ステップＡ５３１のＮｏルート参照）、判定部２１４は、Sが０．０１以下か否か判定する（ステップＡ５３３）。Sが０．０１以下の場合（ステップＡ５３３のＹｅｓルート参照）、評価係数は略０ではないが一定値を保っているものとして、判定部２１４は、運動リズムの同期およびバランスは良いと判定する（ステップＡ５３４）。一方、Sが０．０１より大きい場合（ステップＡ５３３のＮｏルート参照）、判定部２１４は、Cが０．５以上か否か判定する（ステップＡ５３５）。Cが０．５以上の場合（ステップＡ５３５のＹｅｓルート参照）、評価係数は周期的にゆらいでいるものとして、判定部２１４は、運動リズムの同期は良いがバランスは悪いと判定する（ステップＡ５３６）。一方、Cが０．５未満の場合（ステップＡ５３５のＮｏルート参照）、判定部２１４は、Cが０．２以上か否か判定する（ステップＡ５３７）。Cが０．２以上の場合（ステップＡ５３７のＹｅｓルート参照）、判定部２１４は、運動リズムの同期は普通でありバランスは悪いと判定する（ステップＡ５３８）。一方、Cが０．２未満の場合（ステップＡ５３７のＮｏルート参照）、判定部２１４は、運動リズムの同期は悪くバランスは判定不能であると判定する（ステップＡ５３９）。
【００８３】
〔Ｂ〕リアルタイム処理の説明
位相ずれの算出には体動リズムがちょうど２サイクル分変化するだけのデータが必要である。図５（Ａ）及び図５（Ｂ）で例示した歩行データの場合、これは1秒程度なので、データの計測から結果の出力までには少なくとも1秒の時間遅れが生じるが、実用的なリアルタイム処理には十分な許容範囲である。
【００８４】
また、リアルタイム処理のためには、現在の時刻までに測定および蓄積されたデータを用いて解析を行う必要がある。その際問題となるのは、データの境界の効果である。すなわち、スペクトル解析またはヒルベルト変換等の信号処理を施すと、データの端点（測定開始点および終了点）近傍で誤差が大きくなる現象である。つまり、現在の時刻における評価係数は、信号処理に基づく誤差を含んでいる可能性がある。これを見積もるために、図５（Ａ）の加速度データについて時間１２．１秒を現在時刻とし、この時間までのデータのみを用いてヒルベルト変換法による位相変化から評価係数を算出した結果を図１４の破線で示す。１５分間のデータ全体を使用した正確な計算法による結果（図１４中の実線）と比較すると、現在時刻から１．３秒程度さかのぼった時間でほとんど誤差の無い値が得られている。これは実用的なリアルタイム処理としては充分な性能といえる。
【００８５】
このように、本実施形態の一例におけるシステム１によれば、位相に着目し、位相の変化に基づいて、運動リズムの同調性を評価するための評価係数を算出しているので、明瞭なピークを持たないデータからも、評価係数を算出し、運動リズムの同調性を評価することができる。すなわち、本実施形態の一例におけるシステム１によれば、位相に着目しているため、運動リズムのピークに着目しなくても、評価係数を求め、運動リズムの同調性を評価することができる。
【００８６】
また、本実施形態の一例におけるシステム１によれば、運動リズムのピークに着目しなくてもよいため、任意の時間においてほぼ連続的に評価係数を求めることができ、さらには、リアルタイムに評価係数を求めることができる。すなわち、本実施形態の一例におけるシステム１によれば、連続的に運動リズムの同調性を評価することができ、さらには、リアルタイムに運動リズムの同調性を評価することができる。
【００８７】
また、本実施形態の一例におけるシステム１によれば、パターンマッチングを用いることで、自己相関をとった後の波形は元の体動信号に比べてノイズが少なくなっているのでピークが明瞭であり、ピーク位置を正確に特定できるため、評価係数が精度よく求まる。
〔Ｃ〕その他
なお、開示の技術は上述した実施形態に限定されるものではなく、本実施形態の趣旨を逸脱しない範囲で種々変形して実施することができる。
【００８８】
例えば、本実施形態の一例においては、図５（Ａ）において、３軸のうち、上下方向の加速度変化のみを示し、この信号に基づいて、運動リズムを抽出し、評価係数を算出しているが上下方向の加速度変化のみからの信号に限定されるものではない。例えば、前後方向または左右方向の加速度信号から運動リズムを抽出し、評価係数を算出してもよい。
また、本実施形態の一例においては、図２に示したように、所定位置から位相が７２０度変化する範囲において、所定位置から位相が３６０度ずつずれた点を探し、評価係数を算出しているが、３６０度に限定されるものではない。例えば、左右方向の加速度信号から評価係数を算出する場合、左右方向の加速度信号は、図１５に示すように隣り合うピークが同一足の着地に相当するので、所定位置（ピークの位置でなくてもよい）から位相が３６０度変化する範囲において、所定位置から位相が１８０度ずつずれた点を探して、それらの間の時間間隔をそれぞれ第１時間T１，第２時間T２として評価係数を算出する。同様に、体動信号検出部１１としてジャイロセンサにより測定された上下軸周りの角速度信号を用いた場合も、左右方向の加速度信号から評価係数を算出する場合と同様の処理を行なう。
【００８９】
なお、パターンマッチングを用いる場合には、相関値が極小となる点を、所定位置から位相が１８０度ずれた点とする。
さらに、本実施形態の一例では、体動信号検出装置１０は、中央演算部２１をそなえていないが、この構成に限定されるものではなく、体動信号検出装置１０が中央演算部２１をそなえてもよい。この場合、中央演算部２１は、体動信号検出装置１０内の記憶装置（例えば、記憶部１２）または体動信号検出装置１０外の図示しない記憶部に格納された情報処理用プログラムを実行することで、上述の機能を発揮する。すなわち、この場合、体動信号検出装置１０は、情報処理装置としても機能する。
【００９０】
また、本実施形態の一例では、体動信号検出装置１０は、体動信号検出部１１と記憶部１２とをそなえて構成されているが、この構成に限定されるものではない。例えば、体動信号検出装置１０は体動信号検出部１１をそなえるが、記憶部１２をそなえない構成としてもよい。この場合、例えば、体動信号検出装置１０と記憶部１２とは、有線または無線により接続され、体動信号検出部１１により検出された体動信号は、有線または無線を介して記憶部１２に格納される。また、この場合、記憶部１２と情報処理装置２０とは有線または無線により接続され、情報処理装置２０は、記憶部１２から体動信号を取得する。
【００９１】
さらに、体動信号検出装置１０と、情報処理装置２０（もしくは、記憶部１２）とが無線を介して接続される場合、体動信号検出装置１０は、体動信号検出部１１が検出した体動信号をアンテナであるインターフェース部１３を介して情報処理装置２０（もしくは、記憶部１２）に送信する機能を有する。この送信機能は、体動信号検出装置１０がそなえる図示しないたとえば中央演算部である処理部が、図示しない記憶部に格納されたプログラムを実行することで実現される。
【００９２】
また、本実施形態の一例では、情報処理装置２０は、体動信号から第１時定数を決定し、第１時定数に基づいて運動リズムを抽出しているが、この抽出手法に限定されるものではない。例えば、FFT（Fast Fourier Transform）またはウェイブレット解析等のスペクトル解析により体動信号の周波数特性を求め、スペクトル解析の結果から運動リズムの周波数域を特定する。そして、例えば、特定された周波数域に対応するバンドパスフィルタを体動信号にかけて運動リズムを抽出することとしてもよい。
【００９３】
また、体動信号をEMD（Empirical Mode Decomposition）またはEEMD（Ensemble Empirical Mode Decomposition）により各モード波形（Intrinsic Mode Function）に分解し、この結果から、運動リズムに対応するモード波形（例えば、強度の強い波形）を選択することとしてもよい。
さらに、本実施形態の一例では、時定数決定部２１１が時定数を決定した後に、リズム抽出部２１２が、体動信号から体動リズムをフィルタリング処理により抽出しているが、これに限定されるものではない。例えば、上記（３）の処理において、時定数Tを変化させた場合の出力波形Yを全て記憶部２２に記憶しておく。そして、リズム抽出部２１２は、全ての出力波形Yの中から、上記（４）の処理において決定された第１時定数に対応する波形を選択することとしてもよい。
【００９４】
また、本実施形態の一例は、人に限らず、ペット、家畜、馬等の動物に適応可能である。
さらに、本実施形態の一例では、主に随意運動が歩行である場合について述べているが、歩行に限定されるものではなく、本実施形態の一例は他の随意運動にも適用可能である。例えば、随意運動が、ジャグリング（お手玉）であれば、足ではなく、手の運動リズムに着目すればよい。
【００９５】
なお、情報処理装置２０にそなえられる中央演算部２１の各機能を実現するための種々のアプリケーションプログラムは、例えばフレキシブルディスク，ＣＤ（ＣＤ−ＲＯＭ，ＣＤ−Ｒ，ＣＤ−ＲＷ等），ＤＶＤ（ＤＶＤ−ＲＯＭ，ＤＶＤ−ＲＡＭ，ＤＶＤ−Ｒ，ＤＶＤ＋Ｒ，ＤＶＤ−ＲＷ，ＤＶＤ＋ＲＷ，ＨＤ  ＤＶＤ等），ブルーレイディスク，磁気ディスク，光ディスク，光磁気ディスク等の、コンピュータ読取可能な記録媒体に記録された形態で提供される。そして、コンピュータはその記録媒体からプログラムを読み取って内部記憶装置または外部記憶装置に転送し格納して用いる。又、そのプログラムを、例えば磁気ディスク，光ディスク，光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信経路を介してコンピュータに提供するようにしてもよい。
【００９６】
また、体動信号検出装置１０にそなえられる図示しない中央演算部の各機能を実現するための種々のアプリケーションプログラムは、例えばフレキシブルディスク，ＣＤ（ＣＤ−ＲＯＭ，ＣＤ−Ｒ，ＣＤ−ＲＷ等），ＤＶＤ（ＤＶＤ−ＲＯＭ，ＤＶＤ−ＲＡＭ，ＤＶＤ−Ｒ，ＤＶＤ＋Ｒ，ＤＶＤ−ＲＷ，ＤＶＤ＋ＲＷ，ＨＤ  ＤＶＤ等），ブルーレイディスク，磁気ディスク，光ディスク，光磁気ディスク等の、コンピュータ読取可能な記録媒体に記録された形態で提供される。そして、コンピュータはその記録媒体からプログラムを読み取って内部記憶装置または外部記憶装置に転送し格納して用いる。又、そのプログラムを、例えば磁気ディスク，光ディスク，光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信経路を介してコンピュータに提供するようにしてもよい。
【実施例】
【００９７】
以下、本発明を実施例によりさらに詳細に説明するが、本発明は、その要旨を逸脱しない限り、以下の実施例に限定されるものではない。
【実施例１】
【００９８】
体動信号検出装置１０として三菱化学社製加速度レコーダー「見守りゲイト」を専用ベルトに入れてパーキンソン病患者の腹部に巻き、１日の体動信号をサンプリング周波数１００Hzでサンプリングした。レコーダーは腹部中央に位置するようにした。データからオン時とオフ時の歩行時の上下方向加速度信号を取り出し、ヒルベルト変換法とパターンマッチング法により評価係数を求めた。図１６（Ａ）〜図１６（Ｃ）はオン時の歩行結果を示す図である。図１６（Ａ）は、体動信号検出装置１０によって検出された加速度を示す図である。また、図１６（Ｂ）は、２階積分後の加速度信号すなわち運動軌道を示す図である。図１６（Ｃ）は、２階積分後の加速度信号から求められた評価係数を示す図である。図１６（Ｃ）において、評価係数は一貫して規則的なゆらぎを持つので、歩行ステップに左右差はあるが左右の足運びの同期は高いといえる。【００９９】
一方、図１７（Ａ）〜図１７（Ｃ）は、時間が３〜１７秒の間は比較的安定な歩行であるが、その前後で突進歩行が認められたオフ時の歩行結果を示す図である。図１７（Ａ）は、体動信号検出装置１０によって検出された加速度を示す図である。図１７（Ｂ）は、２階積分後の加速度信号すなわち運動軌道を示す図である。図１７（Ｃ）は、２階積分後の加速度信号から求められた評価係数を示す図である。図１７（Ｃ）において、実線がヒルベルト変換法を用いて算出された評価係数を示し、破線がパターンマッチング法を用いて算出された評価係数を示す。図１６（Ａ）〜図１６（Ｃ）と図１７（Ａ）〜図１７（Ｃ）とを比較すると、図１７（Ｃ）においては、安定歩行時での評価係数のゆらぎが不規則になっており、左右のステップ間の同期が崩れていることがわかる。また、図１７（Ｃ）においては、突進歩行時では位相ずれのゆらぎが非常に大きくなっており、左右のステップのバランスが悪化していることがわかる。
【実施例２】
【０１００】
体動信号検出装置１０としてワイヤレステクノロジー社製小型無線ハイブリッドセンサWAA-006を健常被験者の腹部中央に装着し、重いカバンを肩にかけて街中を歩いた際の３軸加速度信号と３軸角速度信号をサンプリング周波数２００Hzで同時計測した。上下方向の加速度信号および上下軸回りの角速度信号についてそれぞれ
・  時定数１を用いてハイパスフィルタリング
・  ２階積分
・  時定数１を用いてハイパスフィルタリング
を施すことで運動リズムを抽出した後、ヒルベルト変換法により位相を求めた。加速度信号については位相差が３６０度となる２点から、角速度信号については位相差が１８０度となる２点から、評価係数を算出した。図１８（Ａ）は、体動信号検出装置１０によって検出された加速度および角速度を示す図である。図１８（Ｂ）は、２階積分後の加速度信号および角速度信号を示す図である。図１８（Ｃ）は、２階積分後の加速度信号および角速度からそれぞれ求められた評価係数を示す図である。図１８（Ａ）〜図１８（Ｃ）において、実線が加速度信号を示し、破線が角速度信号を示す。評価係数は両者でほぼ同じ強度のゆらぎを示している。また、重いカバンを肩にかけていることを反映して、左右のバランスが悪くなっていることがわかる。【０１０１】
このように慣性センサの種類を問わずに体動リズムの同調性を評価することができる。
【実施例３】
【０１０２】
体動信号検出装置１０として三菱化学社製加速度レコーダー「見守りゲイト」を専用ベルトに入れてパーキンソン病患者の腹部に巻き、体動信号をサンプリング周波数１００Hzで38時間連続サンプリングした。被験者には同時に、体の動きやすさの程度と転倒した場合はその時刻を日誌に記入してもらった。【０１０３】
まず、体動信号の時系列データから10秒以上連続して歩いている領域を抽出した。この領域を10秒間隔で時分割し、部分時系列を作成した。それぞれの部分時系列についてフィルタ時定数を最適化し、左右の足運びの同期S2を求めた。時定数の最適化はS2を最小にする方法で行った。得られたS2からexp(-10×S2)を計算し、これを歩行指数とした。歩行指数は左右の足運びの同期が優れているほど大きな値をとり、同期が完璧な場合は１となるような指標である。図１９（Ａ）に平滑化処理を施した後の歩行指数の時間変化を示す。図１９（Ａ）中には、被験者の自己申告による動きやすさの５段階評価（動きやすい、動きにくい、動けないの３段階に、動きやすい−動きにくいの中間評価、及び動きにくい−動けないの中間評価の２段階を加えた５段階評価。）と、実際に転倒した時刻ta，tb，tc，tdが記されている。図１９（Ａ）は、歩行指数の推移を点線で記し、被験者の自己申告の推移を実線で記している。【０１０４】
図１９（Ａ）を参照すると、１日目の９時〜１０時の間は被験者の自己申告は「動きにくい」または「動きやすい」−「動きにくい」の中間評価であり、時刻ta，tbにて転倒が起こっている。これに符合して歩行指数も低い値となっている。１日目の１７時〜１８時、２日目の１０時〜１１時の間は、被験者の自己申告は「動きやすい」であるが、歩行指数は前後の時間帯に比べ低下しており、実際に転倒が時刻tc，tdにて起こっている。しかも転倒時刻よりかなり前から歩行指数が低下傾向にあることがわかる。【０１０５】
このように、歩行リズムの同調性は、転倒危険性を予測する指標となる。【０１０６】
次に、同じ時系列データを10秒間隔で時分割し、部分時系列を作成した。ここでは歩行領域を取り出すことをしない。それぞれの部分時系列についてフィルタ時定数を最適化し、左右の足運びの同期S2および位相の周期のCVを求めた。時定数の最適化はS2を最小にする方法で行った。得られたS2からexp(-10×S2)を計算し、これを体動指数とした。このうち、位相の周期のCVが0.12以上になるものを除去した。これは周期性の悪い部分時系列を除外し、規則的なリズム（特に歩行リズム）の性質に着目することを意味する。図１９（Ｂ）に平滑化処理を施した後の体動指数の時間変化を示す。図１９（Ｂ）中には、被験者の自己申告による動きやすさの５段階評価（動きやすい、動きにくい、動けないの３段階に、動きやすい−動きにくいの中間評価、及び動きにくい−動けないの中間評価の２段階を加えた５段階評価。）と、実際に転倒した時刻がta，tb，tc，tdとが記されている。図１９（Ｂ）は、体動指数の推移を点線で記し、被験者の自己申告の推移を実線で記している。【０１０７】
図１９（Ｂ）を参照すると、転倒のあった領域近辺では、図１９（Ａ）と同じような振る舞いを示している。すなわち、１日目の９時〜１０時の間は被験者の自己申告は「動きにくい」または「動きやすい」−「動きにくい」の中間評価であり、時刻ta，tbにて転倒が起こっている。これに符合して体動指数も低い値となっている。１日目の１７時〜１８時、２日目の１０時〜１１時の間は、被験者の自己申告は「動きやすい」であるが、体動指数は前後の時間帯に比べ低下しており、実際に転倒が時刻tc，tdにて起こっている。しかも転倒時刻よりかなり前から体動指数が低下傾向にあることがわかる。【０１０８】
このように、体動リズムの同調性は転倒危険性を予測する指標となる。【０１０９】
本発明を詳細にまた特定の実施態様を参照して説明したが、本発明の精神と範囲を逸脱することなく様々な変更または修正を加えることができることは当業者にとって明らかである。【０１１０】
本出願は、2010年9月17日出願の日本特許出願（特願2010−209112）に基づくものであり、その内容はここに参照として取り込まれる。
【産業上の利用可能性】
【０１１１】
本発明の情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体によれば、運動リズムの同調性を連続的かつリアルタイムで評価することができる、という効果を奏する。
【符号の説明】
【０１１２】
１  システム
１０  体動信号検出装置
１１  体動信号検出部
１２  記憶部
１３  インターフェース部
２０  情報処理装置
２１  中央演算部
２２  記憶部
２３  出力部
２４  インターフェース部
２１１  時定数決定部
２１２  リズム抽出部
２１３  評価係数決定部
２１４  判定部
２１５  出力制御部
２２３  時刻算出部
２３３  評価係数算出部
【図１】
2012036135.tif 000008
【図２】
2012036135.tif 000009
【図３】
2012036135.tif 000010
【図４】
2012036135.tif 000011
【図５】
2012036135.tif 000012
【図６】
2012036135.tif 000013
【図７】
2012036135.tif 000014
【図８】
2012036135.tif 000015
【図９】
2012036135.tif 000016
【図１０】
2012036135.tif 000017
【図１１】
2012036135.tif 000018
【図１２】
2012036135.tif 000019
【図１３】
2012036135.tif 000020
【図１４】
2012036135.tif 000021
【図１５】
2012036135.tif 000022
【図１６】
2012036135.tif 000023
【図１７】
2012036135.tif 000024
【図１８】
2012036135.tif 000025
【図１９】
2012036135.tif 000026
【国際調査報告】
2012036135.tif 000027
2012036135.tif 000028
2012036135.tif 000029
2012036135.tif 000030