(57)【要約】【課題】人間等の動物の生理状態に応じて、その運動を適当なスケールおよびリズムで誘導しうるように当該動物の運動を管理するシステム等を提供する。【解決手段】本発明の運動管理システム１によれば、歩行運動誘導装置２０の動作によりユーザの歩行運動が誘導されている過程で、歩行比ｋおよび消費エネルギー等の生理変数が測定される。そして、このユーザの身体機能の活性化または身体負荷の軽減に鑑みて生理変数が適当な値となった状態における歩行比ｋが推奨歩行比ｋ0として設定される。これにより、推奨歩行比ｋ0がこのユーザの身体機能の活性化等に鑑みて適当に設定されうる。そして、歩行運動誘導装置２０によりこの推奨歩行比ｋ0に基づいて歩行運動が誘導されることで、歩幅ｑおよび歩行率ｐがこのユーザの身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。【選択図】  図１
2007275283.tif 000002
【特許請求の範囲】
【請求項１】
動物の運動を管理するシステムであって、
動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら該動物の運動を誘導する運動誘導装置（以下、適宜「第２運動誘導装置」という。）の動作を制御する制御部と、
動物の生理状態に応じた生理変数を測定する生理変数測定部と、
動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定部と、
制御部により運動変数が変化するように制御されている第２運動誘導装置の動作によって動物の運動が誘導されている過程で、生理変数測定部により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、運動変数測定部により測定された運動変数を該運動変数の推奨値として設定する推奨値設定部とを備えていることを特徴とする運動管理システム。
【請求項２】
前記制御部が、動物の運動振動子を測定する運動振動子測定部と、
運動振動子測定部により測定された運動振動子と、固有角速度が反映された形で相互に引き込み合う１次振動子を生成する１次振動子生成部と、
運動振動子測定部により測定された運動振動子と１次振動子生成部により生成された１次振動子との位相差に基づき、新たな固有角速度を設定する固有角速度設定部と、
固有角速度設定部により設定された新たな固有角速度が反映されたリズムで振動する２次振動子を生成する２次振動子生成部と、
２次振動子生成部により生成された２次振動子に基づいて第２運動誘導装置の動作のスケールおよびリズムを特定する誘導振動子を生成する誘導振動子生成部とを備えていることを特徴とする請求項１記載の運動管理システム。
【請求項３】
生理変数測定部が、動物の酸素摂取量、呼吸周波数、心拍数、脈拍数、血圧、血中飽和酸素濃度、乳酸値、筋電位および消費エネルギーのうち一部または全部を生理変数として測定することを特徴とする請求項１記載の運動管理システム。
【請求項４】
運動変数測定部が、
動物の運動の速さを調節しながら第２運動誘導装置とは別に該動物の運動を誘導する運動誘導装置（以下、適宜「第１運動誘導装置」という。）の動作により運動が誘導されている動物の運動の速さを表す第１運動変数を第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定部と、
動物の運動のリズムまたはスケールを表す第２運動変数を第１運動誘導装置と該動物との相互作用状態または第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定部とにより構成され、かつ、
第１および第２運動変数の関数を前記運動変数として測定することを特徴とする請求項１記載の運動管理システム。
【請求項５】
第１運動変数測定部が動物の歩行または走行速度を第１運動変数として測定し、第２運動変数測定部が動物の歩幅または歩行率（単位時間当たりの歩数を意味する。）を第２運動変数として測定することを特徴とする請求項４記載の運動管理システム。
【請求項６】
運動変数測定部が第１運動変数としての歩行または走行速度に対する第２運動変数としての歩幅の二乗の比率、または第２運動変数としての歩行率の二乗に対する第１運動変数としての歩行または走行速度の比率である歩行比を前記運動変数として測定することを特徴とする請求項５記載の運動管理システム。
【請求項７】
第１運動変数測定部が、第１運動誘導装置が有する回動体の回動方向に逆らう方向に運動が誘導されている動物の第１運動変数を、当該回動体の回動速度に基づいて測定することを特徴とする請求項４記載の運動管理システム。
【請求項８】
第１運動変数測定部が、第１運動誘導装置としてのトレッドミルが有する複数のローラに掛け渡された回動体としての無端ベルトの動きに逆らう方向に誘導されている動物の歩行または走行運動の速さを表す第１運動変数としての歩行または走行速度を、当該無端ベルトの駆動速度に基づいて測定することを特徴とする請求項４記載の運動管理システム。
【請求項９】
推奨値設定部が、複数の動物のそれぞれの身体的特徴を識別するための識別子を認識した上で各動物の複数の推奨値を予備的に設定し、認識した識別子ごとに予備的に設定済みの推奨値に基づいて新たな推奨値を設定することにより、識別子および推奨値が対応付けられているデータベースを構築することを特徴とする請求項１記載の運動管理システム。
【請求項１０】
動物の運動を管理する方法であって、
動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら、動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数が変化するように該動物の運動を誘導し、
動物の生理状態に応じた生理変数を測定し、
動物の運動変数を測定し、
測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で測定された運動変数を該運動変数の推奨値として設定することを特徴とする運動管理方法。
【請求項１１】
動物の運動を管理する機能をコンピュータに付与するプログラムであって、
動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら該動物の運動を誘導する運動誘導装置の動作を制御する制御機能と、
動物の生理状態に応じた生理変数を測定する生理変数測定機能と、
動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定機能と、
制御機能により運動変数が変化するように制御されている第２運動誘導装置の動作によって動物の運動が誘導されている過程で、生理変数測定機能により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、運動変数測定機能により測定された運動変数を該運動変数の推奨値として設定する推奨値設定機能とをコンピュータに付与することを特徴とする運動管理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動物の運動を管理するシステム、方法および当該管理機能をコンピュータに付与するプログラムに関する。
【背景技術】
【０００２】
ユーザの歩行訓練のレベルを設定する等の目的のため、このユーザの心電信号等の生理指標に鑑みて適当な歩行速度を決定する手法が提案されている（たとえば、特許文献１参照）。
【特許文献１】特開２００３−１５４０２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、ユーザの身体に力を作用させることで脚等の動きを誘導（または補助）する装置がユーザに装着されている場合、このユーザにとって適当な誘導スケールおよびリズムを設定することは考慮されていなかった。
【０００４】
そこで、本発明は、人間等の動物の生理状態に応じて、その運動を適当なスケールおよびリズムで誘導しうるように当該動物の運動を管理するシステム、方法、および当該管理機能をコンピュータに付与するプログラムを提供することを解決課題とする。
【課題を解決するための手段】
【０００５】
前記課題を解決するための本発明の運動管理システムは、動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら該動物の運動を誘導する運動誘導装置（以下、適宜「第２運動誘導装置」という。）の動作を制御する制御部と、動物の生理状態に応じた生理変数を測定する生理変数測定部と、動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定部と、制御部により運動変数が変化するように制御されている第２運動誘導装置の動作によって動物の運動が誘導されている過程で、生理変数測定部により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、運動変数測定部により測定された運動変数を該運動変数の推奨値として設定する推奨値設定部とを備えていることを特徴とする。
【０００６】
本発明の運動管理システムによれば、第２運動誘導装置の動作により動物の運動が誘導されている過程で「生理変数」および「運動変数」が測定される。動物の酸素摂取量、呼吸周波数、心拍数、脈拍数、血圧、血中飽和酸素濃度、乳酸値、筋電位、および消費エネルギーのうち一部または全部が「生理変数」として測定されうる。そして、この動物の「身体機能の活性化」または「身体負荷の軽減」に鑑みて生理変数が適当な値となった状態における運動変数の値が当該運動変数の推奨値として設定される。これにより、動物の運動変数の推奨値が、この動物の身体機能の活性化等に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の運動スケールおよびリズムがこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。すなわち、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうる。
【０００７】
また、本発明の運動管理システムは、前記制御部が、動物の運動振動子を測定する運動振動子測定部と、運動振動子測定部により測定された運動振動子と、固有角速度が反映された形で相互に引き込み合う１次振動子を生成する１次振動子生成部と、運動振動子測定部により測定された運動振動子と１次振動子生成部により生成された１次振動子との位相差に基づき、新たな固有角速度を設定する固有角速度設定部と、固有角速度設定部により設定された新たな固有角速度が反映されたリズムで振動する２次振動子を生成する２次振動子生成部と、２次振動子生成部により生成された２次振動子に基づいて第２運動誘導装置の動作のスケールおよびリズムを特定する誘導振動子を生成する誘導振動子生成部とを備えていることを特徴とする。
【０００８】
本発明の運動管理システムによれば、動物の運動リズムと、第２運動誘導装置の動作リズム（動物の運動を誘導するリズム）との相互の歩み寄りまたは調和が図られる。これにより、動物の運動および第２運動誘導装置の動作の調和を図りながら、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうる。なお、ここで「振動」とは、現実的または仮想的なもの（物体）がほぼ一定の周期をもって揺れ動くことのほか、広く時間変化することを意味する。また「振動子」とは、値が時間変化する電気信号や、ソフトウェアにおいて値が時間変化するものとして定義される関数等を意味する。
【０００９】
また、本発明の運動管理システムは、運動変数測定部が、動物の運動の速さを調節しながら第２運動誘導装置とは別に該動物の運動を誘導する運動誘導装置（以下、適宜「第１運動誘導装置」という。）の動作により運動が誘導されている動物の運動の速さを表す第１運動変数を第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定部と、動物の運動のリズムまたはスケールを表す第２運動変数を第１運動誘導装置と該動物との相互作用状態または第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定部とにより構成され、かつ、第１および第２運動変数の関数を前記運動変数として測定することを特徴とする。
【００１０】
本発明の運動管理システムによれば「第１運動変数」が第１運動誘導装置の動作の速さに基づいて測定される。また「第２運動変数」が第１運動誘導装置と動物との相互作用または第２運動誘導装置による誘導スケールもしくはリズムに基づいて測定される。これにより、第１および第２運動変数、ひいてはこれらの関数である運動変数の測定精度の向上が図られる。したがって、動物の身体機能の活性化等に鑑みて、より適当な運動変数の推奨値が設定されうる。
【００１１】
さらに、本発明の運動管理システムは、第１運動変数測定部が動物の歩行または走行速度を第１運動変数として測定し、第２運動変数測定部が動物の歩幅または歩行率（単位時間当たりの歩数を意味する。）を第２運動変数として測定することを特徴とする。
【００１２】
本発明の運動管理システムによれば、動物の歩行または走行の速さ（速度）および歩幅（スケール）または歩行率（リズム）の関数である運動変数の推奨値が、この動物の身体機能の活性化や身体負荷の軽減に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の歩幅および歩行率がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００１３】
また、本発明の運動管理システムは、運動変数測定部が第１運動変数としての歩行または走行速度に対する第２運動変数としての歩幅の二乗の比率、または第２運動変数としての歩行率の二乗に対する第１運動変数としての歩行または走行速度の比率である歩行比を前記運動変数として測定することを特徴とする。
【００１４】
本発明の運動管理システムによれば、動物の歩行比の推奨値が、この動物の身体機能の活性化や身体負荷の軽減に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の歩行比がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００１５】
さらに、本発明の運動管理システムは、第１運動変数測定部が、第１運動誘導装置が有する回動体の回動方向に逆らう方向に運動が誘導されている動物の第１運動変数を当該回動体の回動速度に基づいて測定することを特徴とする。
【００１６】
本発明の運動管理システムによれば、第１運動誘導装置が有する回動体を回動させることにより、回動体と身体部分が接触する動物の、回動体の回動方向とは逆向きの運動が誘導される。「回動体」には、複数のローラに掛け渡された無端ベルト、中心または中心から外れた点を通る軸回りに回動される球体、楕円球体、中心軸または中心軸から離れた中心軸に平行な軸回りに回動する円筒、四角筒等の筒体、その他任意の軸回りに回動される一塊の物体が含まれる。これにより、当該運動に伴う動物またはその身体部分の移動が回動体の回動によって相殺されるので、第１運動誘導装置の設置スペースさえあれば動物の運動変数の推奨値設定が可能であり、当該設定スペースが比較的自由に選択されうる。
【００１７】
また、本発明の運動管理システムは、第１運動変数測定部が、第１運動誘導装置としてのトレッドミルが有する複数のローラに掛け渡された回動体としての無端ベルトの動きに逆らう方向に誘導されている動物の歩行または走行運動の速さを表す第１運動変数としての歩行または走行速度を、当該無端ベルトの駆動速度に基づいて測定することを特徴とする。
【００１８】
本発明の運動管理システムによれば、トレッドミル（第１運動誘導装置）の無端ベルト（回動体）を回動させることにより、この無端ベルトの動きとは逆向きの歩行または走行が誘導される。これにより、当該歩行または走行に伴う動物の移動が無端ベルトの回動によって相殺されるので、第１運動誘導装置の設置スペースさえあれば動物の歩行または走行運動に伴う運動変数の推奨値設定が可能であり、当該設定スペースが比較的自由に選択されうる。
【００１９】
さらに、本発明の運動管理システムは、推奨値設定部が複数の動物のそれぞれの身体的特徴を識別するための識別子を認識した上で各動物の複数の推奨値を予備的に設定し、認識した識別子ごとに予備的に設定済みの推奨値に基づいて新たな推奨値を設定することにより、識別子および推奨値が対応付けられているデータベースを構築することを特徴とする。
【００２０】
本発明の運動管理システムによれば、個々の動物の身体的特徴（身体のサイズ、ウェイト、年齢、性別等が含まれる。）に応じて標準的な推奨値が設定される。そして、２運動誘導装置によりこの標準的な推奨値に基づいて動物の運動が誘導されることで、動物の運動変数がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００２１】
前記課題を解決するための本発明の運動管理方法は、動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら、動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数が変化するように該動物の運動を誘導し、動物の生理状態に応じた生理変数を測定し、動物の運動変数を測定し、測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で測定された運動変数を該運動変数の推奨値として設定することを特徴とする。
【００２２】
本発明の運動管理方法によれば、前記運動管理システムと同様に、これにより、動物の運動変数の推奨値が、この動物の身体機能の活性化等に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の運動スケールおよびリズムがこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００２３】
前記課題を解決するための本発明の運動管理プログラムは、動物の運動リズムと調和したリズムで該動物の運動スケールを調節しながら該動物の運動を誘導する運動誘導装置の動作を制御する制御機能と、動物の生理状態に応じた生理変数を測定する生理変数測定機能と、動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定機能と、制御機能により運動変数が変化するように制御されている第２運動誘導装置の動作によって動物の運動が誘導されている過程で、生理変数測定機能により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、運動変数測定機能により測定された運動変数を該運動変数の推奨値として設定する推奨値設定機能とをコンピュータに付与することを特徴とする。
【００２４】
本発明の運動管理プログラムによれば、これにより、動物の運動変数の推奨値が、この動物の身体機能の活性化等に鑑みて適当に設定されうるように、動物の運動を管理する機能がコンピュータに付与される。なお、当該複数の機能のうち全てが一のコンピュータに付与されてもよく、当該複数の機能が複数のコンピュータに分散して付与されてもよい。
【発明を実施するための最良の形態】
【００２５】
本発明の運動管理システム、運動管理方法および運動管理プログラムの実施形態について図面を用いて説明する。
図１は本発明の運動管理システムの構成例示図であり、図２〜図３は本発明の運動管理方法の説明図であり、図４は歩行運動を誘導する仮想的なバネおよびダンパの説明図であり、図５は推奨歩行比の設定方法の説明図であり、図６は本発明の運動管理方法の説明図であり、図７〜図９は本発明の運動管理システムの効果に関する実験結果の説明図である。
【００２６】
以下、歩行者の脚体等について左右を区別するためにパラメータに添字Ｌ、Ｒを添付するが、表記の簡単のため左右を区別する必要が特にない場合には添字Ｌ、Ｒを省略する。
【００２７】
本発明の運動管理システムの構成について図１を用いて説明する。
【００２８】
図１に示されている運動管理システム１は、トレッドミル（第１運動誘導装置）１０および歩行運動誘導装置（第２運動誘導装置）を用いて、ユーザ（人間（動物））の身体機能の活性化等に鑑みて適当な推奨歩行比を設定するためのものである。また。運動管理システム１は、ユーザの歩行（または走行）運動と、トレッドミル１０の動作と、歩行運動誘導装置２０の動作との調和を図りながら、推奨歩行比に応じたスケールおよびリズムで歩行するようにこのユーザを訓練するためのものである。
【００２９】
トレッドミル１０は、人間の標準的な横幅よりも若干広い程度の幅の駆動ローラ１１と、駆動ローラ１１とほぼ同じ幅の従動ローラ１２と、駆動ローラ１１および従動ローラ１２に掛け渡されている無端ベルト（回動体）１３と、無端ベルト１３のうちユーザが乗る部分を下から支える踏板１４とを備えている。駆動ローラ１１は、モータ、変速機等より構成される駆動機構１０１によって駆動される。駆動ローラ１１が図中時計回りに駆動されることにより、従動ローラ１２の同方格への従動を伴ってベルト１３も時計回りに回る。これによりベルト１３に乗っているユーザの図中左への歩行（または走行）が誘導される。また、ベルト１３の移動速度に応じた信号を出力する速度センサ１０２と、踏板１４が受ける圧力に応じた信号を出力する圧力センサ１０４とが設けられている。ユーザの身体には、その心拍数等の生理変数ｂに応じた信号を出力する生理変数センサ１０６が取り付けられている。生理変数センサ１０６は生理変数ｂの内容によって適当な方法でユーザの身体に取り付けられる。なお、トレッドミル１０としては、一般に市販されているもの等、公知のあらゆる構成のトレッドミルが採用されてもよい。
【００３０】
歩行運動誘導装置２０は、腰部装具２１と、大腿部装具２２と、力伝達部材２３と、バッテリ２４と、アクチュエータ（電動モータ）２５と、股関節角度センサ２６とを備えている。
【００３１】
腰部装具２１は剛性のある素材と柔軟性のある素材とが組み合わせられて作られており、ユーザの腰部に装着される。大腿部装具２２も剛性のある素材と柔軟性のある素材とが組み合わせられて作られており、ユーザの大腿部の前後それぞれに装着される。力伝達部材２３は、軽量の硬質プラスチック等の定形性のある素材より作られており、ユーザの大腿部に沿って、ユーザの腰部の横から下方に延びた後で大腿部の前後に向けて二股に分かれた形状であり、アクチュエータ２５および前後の大腿部装具２２のそれぞれに連結されている。バッテリ２４は腰部装具２１に収納されており（たとえば、腰部装具２１を構成する複数枚の素材の間に固定されており）、アクチュエータ２５等に対して電力を供給する。アクチュエータ２５は腰部装具２１に取り付けられており、力伝達部材２３および大腿部装具２２を介してユーザの大腿部に力を作用させる。股関節角度センサ２６はユーザの腰部の横に設けられたロータリエンコーダ等により構成され、股関節角度に応じた信号を出力する。
【００３２】
運動管理システム１は、第１制御部１００と、第２制御部２００とを備えている。
【００３３】
第１制御部１００は、トレッドミル１０に付属するマイクロコンピュータ等のコンピュータにより構成されており、駆動機構１０２による駆動ローラ１１の駆動速度等を制御する。第１制御部１００は、運動変数測定部１１０と、生理変数測定部１２０と、推奨歩行比設定部１３０とを備えている。運動変数測定部１１０、生理変数測定部１２０および推奨歩行比設定部１３０はそれぞれハードウェアとしてのコンピュータ（ＣＰＵ、ＲＯＭ、ＲＡＭ、Ｉ／Ｏ（入出力装置）等により構成されている。）と、このコンピュータに機能を付与するソフトウェアとしての本発明の「運動管理プログラム」の一部とにより構成されている。
【００３４】
運動変数測定部１１０は第１運動変数測定部１１１と、第２運動変数測定部１１２とを備えている。
【００３５】
第１運動変数測定部１１１は、速度センサ１０２の出力に基づいてユーザの歩行速度ｖを当該ユーザの運動の速さを表す第１運動変数として測定する。第２運動変数測定部１１２は、圧力センサ１０４の出力に基づいてユーザの歩行率（単位時間当たりの歩数）ｐを当該ユーザの運動のリズムを表す第２運動変数として測定する。そして、運動変数測定部１１０は、歩行速度ｖおよび歩行率ｐの関数である運動変数として歩行比ｋ（＝ｖ／ｐ2）を測定する。
【００３６】
生理変数測定部１２０は、ユーザの身体に取り付けられた生理変数センサ１０６の出力に基づいてこのユーザの生理状態に応じた生理変数ｂを測定する。
【００３７】
推奨歩行比設定部１３０は、運動変数測定部１１０により測定されたユーザの歩行比ｋと、生理変数測定部１２０により測定された生理変数ｂとに基づき、推奨歩行比ｋ0を設定する。
【００３８】
第２制御部（本発明の「制御部」に該当する。）２００は、歩行運動誘導装置２０の腰部装具２１に収納されたコンピュータと、このコンピュータに対して歩行運動誘導装置２０の制御機能等を付与するソフトウェアとしての本発明の「運動管理プログラム」とにより構成されている。
【００３９】
第２制御部２００は、運動振動子測定部２１０と、１次振動子生成部２２０と、固有角速度設定部２３０と、２次振動子生成部２４０と、誘導振動子生成部２５０とを備えている。
【００４０】
運動振動子測定部２１０は、股関節角度センサ２６の出力に基づき、ユーザの運動スケールに応じた運動振動子として股関節角度φHを測定する。また、運動振動子測定部２１０は、股関節角度センサ２６の出力に基づき、ユーザの運動スケールに応じた運動振動子の時間微分値として股関節角速度ｄφH／ｄｔを測定する。
【００４１】
１次振動子生成部２２０は運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔと、固有角速度ωMとに基づき、「第１モデル」にしたがって１次振動子ｘを生成する。
【００４２】
固有角速度設定部２３０は、第１位相差設定部２３１と、第２位相差設定部２３２と、相関係数設定部２３３と、第１角速度設定部２３４と、第２角速度設定部２３５とを備えている。
【００４３】
第１位相差設定部２３１は股関節角速度ｄφH／ｄｔの角速度ωHと、ファン・デル・ポル方程式に含まれる固有角速度ωMが反映された振動子ｘとの位相差を第１位相差δθ1として設定する。
【００４４】
第２位相差設定部２３２は仮想運動振動子θhと仮想誘導振動子（仮想誘導振動子）θmとの関係を表す「仮想モデル」にしたがって、仮想運動振動子θhと仮想誘導振動子θmとの位相差を第２位相差δθ2（＝θh−θm）として設定する。
【００４５】
相関係数設定部２３３は、第２位相差設定部２３２により設定された第２位相差δθ2が、第１位相差設定部２３１により設定された第１位相差δθ1に近づくように仮想運動振動子θhと仮想誘導振動子θmとの相関係数εを設定する。
【００４６】
第１角速度設定部２３４は、相関係数設定部２３３により設定された相関係数εに基づき、仮想運動振動子θhの角速度ωhを設定する。
【００４７】
第２角速度設定部２３５は、第１角速度設定部２３４により設定された仮想運動振動子θhの角速度ωhに基づき、第２位相差設定部２３２により設定された第２位相差δθ2が、目標位相差設定部２１２により設定された目標位相差δθｄに近づくように仮想誘導振動子θmの角速度ωmを新たな固有角速度ωMとして設定する。
【００４８】
２次振動子生成部２４０は、運動振動子測定部２１０により測定された股関節角度φHと、固有角速度設定部２３０により設定された新たな固有角速度ωMとに基づき「第２モデル」にしたがって２次振動子ｙを生成する。
【００４９】
誘導振動子生成部２５０は、第１誘導振動子生成部２５１と、第２誘導振動子生成部２５２とを備えている。
【００５０】
第１誘導振動子生成部２５１は２次振動子生成部２４０により生成された２次振動子ｙと、固有角速度設定部２３０により設定された固有角速度ωMと、運動振動子測定部２１０により測定された股関節角度φHをその目標値φ0に近付ける第１ポテンシャル（後述）とに基づく第１誘導振動子ｚ1を生成する。第２誘導振動子生成部２５２は２次振動子生成部２４０により生成された２次振動子ｙと、固有角速度設定部２３０により設定された固有角速度ωMと、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔに応じて股関節角度φHの絶対値の増大を抑制する第２ポテンシャル（後述）に基づく第２誘導振動子ｚ2を生成する。そして、誘導振動子生成部２５０は、第１誘導振動子ｚ1および第２誘導振動子ｚ2に基づき、歩行運動誘導装置２０によりユーザに作用させられる股関節回りのトルクＴを特定する誘導振動子ｚを生成する。
【００５１】
第１制御部１００および第２制御部２００は無線通信が可能とされている。なお、第１制御部１００に代えて第２制御部２００が推奨歩行比設定部１３０を備えている等、運動管理システムを構成する複数の処理部１１０（１１１，１１２）１，１２０，１３０，２１０，２２０，‥が、第１制御部１００および第２制御部２００に任意のパターンで配設されていてもよい。また、第１制御部１００および第２制御部２００は同一のコンピュータにより構成されていてもよい。第１制御部１００および第２制御部２００の通信は有線通信であってもよい。
【００５２】
前記構成の運動管理システム１の機能、特にユーザの推奨歩行比（運動変数の推奨値）の設定方法について図２〜図５を用いて説明する。
【００５３】
第１制御部１００により、ベルト１３の速度が一定になるようにトレッドミル１０の動作が制御されることで、ベルト１３の動きとは逆方向へのほぼ一定速度でのユーザの歩行運動が誘導される。
【００５４】
この状態で、運動振動子測定部２１０は股関節角度センサ２６の出力に基づき、ユーザの左右の股関節角度（ユーザの運動スケールに応じた運動振動子）φH＝（φHL，φHR）を測定する（図２／ｓ０１１）。また、運動振動子測定部２１０は股関節角度センサ２６の出力に基づき、ユーザの左右の股関節角速度（ユーザの運動スケールに応じた運動振動子の時間微分値）ｄφH／ｄｔ＝（ｄφHL／ｄｔ，ｄφHR／ｄｔ）を測定する（図２／ｓ０１２）。
【００５５】
さらに、１次振動子生成部２２０が、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔおよびメモリに記憶されている最新の固有角速度ωM＝（ωML，ωMR）に基づき「第１モデル」にしたがって１次振動子ｘ＝（ｘL，ｘR）を設定する（図２／ｓ０２０）。第１モデルは、股関節角速度ｄφH／ｄｔ等の運動振動子に応じて出力が変動する仮想的な左右の脚体等、仮想的な複数の第１要素の関係を、次式（１）で表されるファン・デル・ポル（ｖａｎ  ｄｅｒ  Ｐｏｌ）方程式によって表現するものである。
【００５６】
（ｄ2ｘL／ｄｔ2）＝ξ（１−ｘL2）（ｄｘL／ｄｔ）−ωML2ｘL
＋ｇ（ｘL−ｘR）＋Ｋ（ｄφHL／ｄｔ），
（ｄ2ｘR／ｄｔ2）＝ξ（１−ｘR2）（ｄｘR／ｄｔ）−ωMR2ｘR
＋ｇ（ｘR−ｘL）＋Ｋ（ｄφHR／ｄｔ）  ‥（１）
ここで「ξ」は１次振動子ｘおよびその１回時間微分（ｄｘ／ｄｔ）がｘ−（ｄｘ／ｄｔ）平面で安定なリミットサイクルを描くように設定される係数（＞０）である。ｇは第１モデルにおける仮想的な左右の脚体（第１要素）の相関関係を表す第１相関係数である。Ｋはフィードバック係数である。なお、固有角速度ωMは、歩行運動誘導装置２０による実際の歩行補助リズム（歩行誘導リズム）から大きく外れない範囲で任意に設定されてよい。
【００５７】
１次振動子ｘ＝（ｘL，ｘR）はルンゲ・クッタ法により算定（設定）される。１次振動子ｘの成分ｘLおよびｘRはそれぞれ左右の脚体に関する仮想的な歩行補助リズムを表す。また、振動子ｘはファン・デル・ポル方程式の１つの性質である「相互引き込み」により、実際の歩行運動リズムとほぼ同じリズム（角速度）で時間変化する股関節角速度ｄφH／ｄｔのリズムと調和しながらも「固有角速度」ωMが反映された自律的なリズムをもって時間変化に応じて振動するという性質がある。
【００５８】
なお、股関節角速度ｄφH／ｄｔに代えてまたは加えて、股関節角度φHや、膝関節、足関節、肩関節、肘関節の角度や角速度、さらには歩行者の着地音、呼吸音、断続的な発声音等、ユーザの歩行運動リズム（運動リズム）が反映されたリズムで変動する種々の振動子に基づき、１次振動子ｘが生成されてもよい。
【００５９】
また、式（１）で表現されるファン・デル・ポル方程式とは異なる形のファン・デル・ポル方程式によって第１モデルが表現されてもよく、股関節角速度ｄφH／ｄｔ等の運動振動子と相互引き込み効果をもって振動子が生成され得るあらゆる方程式によって第１モデルが表現されてもよい。
【００６０】
前記のように股関節角速度ｄφH／ｄｔ等の運動振動子に応じて出力が変動する仮想的な複数の第１要素（左右の脚体）の関係を表現する第１モデルにしたがって、第１要素の出力として１次振動子ｘが生成される（式（１），図２／ｓ０２０）。これにより、ユーザの実際の運動に関係する現実の複数の要素の関係を第１モデルにおける仮想的な複数の第１要素の係数ｇ等により表される相関関係に反映させることで、当該現実の複数の要素間の関係に鑑みて適当な１次振動子ｘが生成され得る。たとえば、現実の複数の要素として左右の脚体や同一脚体の複数の関節が想定された場合、交互に前後に動く等の左右の脚体の定性的関係や股関節回りの脚体運動と膝関節回りの脚体運動との周期や位相差等の同一脚体の関節間の定性的関係等が反映された形で１次振動子ｘが生成される。したがって、ユーザの運動を誘導する誘導振動子のリズムおよびスケールを当該関係に鑑みて適当なものとすることができる。
【００６１】
続いて、固有角速度設定部２３０がメモリに記憶されている目標位相差δθdと、１次振動子設定部２１０により生成された１次振動子ｘとに基づき、仮想的な２つの振動子が含まれている仮想モデルにしたがって固有角速度ωMを新たに設定する（図２／ｓ０３０）。
【００６２】
具体的には、まず、第１位相差測定部２３１が、左右各成分について、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔの位相θHと、１次振動子設定部２１０により生成された１次振動子ｘの位相θMとの位相差θH−θMを、第１位相差δθ1として設定する（図２／ｓ０３１）。
【００６３】
次に、第２位相差設定部２３２が、過去３歩行周期にわたって第１位相差設定部２３１により設定された第１位相差δθ1が一定であったことを要件として、左右各成分について、次式（２．１）および（２．２）によって表される「仮想モデル」にしたがって、次式（２．３）によって表される仮想運動振動子θhと仮想誘導振動子θmとの位相差θh−θmを第２位相差δθ2として設定する（図２／ｓ０３２）。
【００６４】
ｄθh／ｄｔ＝ωh＋ε・ｓｉｎ（θmL−θhL）  ‥（２．１）
ｄθm／ｄｔ＝ωm＋ε・ｓｉｎ（θhL−θmL）  ‥（２．２）
δθ2＝ａｒｃｓｉｎ［（ωh−ωm）／２ε］  ‥（２．３）
ここで、ε＝（εL，εR）は仮想モデルにおける仮想運動振動子θh＝（θhL，θhR）および仮想誘導振動子θm＝（θmL，θmR）の左右成分ごとの相関係数である。また、ωhは仮想運動振動子θhの角速度であり、ωmは仮想誘導振動子θmの角速度である。
【００６５】
続いて、相関係数設定部２３３が、第１位相差設定部２３１により設定された第１位相差δθ1と、第２位相差設定部２３２により設定された第２位相差δθ2との差δθ1−δθ2が最小になるように、相関係数εを設定する（図２／ｓ０３３）。
【００６６】
具体的には次式（２．４）にしたがって、左右各成分について、股関節角速度（運動振動子）ｄφH／ｄｔが０となる離散的な時間ｔid（ｄ＝１，２，・・）における相関係数εが逐次設定される。
【００６７】
ε（ｔid+1）
＝ε（ｔid）−η｛Ｖ1（ｔid+1）−Ｖ1（ｔid）｝
／｛ε（ｔid）−ε（ｔid-1）｝，
Ｖ（ｔid+1）≡（１／２）｛δθ1（ｔid+1）−δθ2（ｔid）｝2  ‥（２．４）
ここで、η＝（ηL，ηR）の各成分は、第１位相差δθ1の左右各成分と、第２位相差δθ2の左右各成分とを近づけるポテンシャルＶ＝（ＶL，ＶR）の安定性を表す係数である。
【００６８】
次に、第１角速度設定部２３４が、相関係数設定部２３３により設定された相関係数εに基づき、仮想誘導振動子θmの固有角速度ωmが一定であるという条件下で、左右各成分について、第１および第２位相差の差δθ1−δθ2の各成分が最小となるように仮想運動振動子θhの角速度ωhを次式（２．５）にしたがって設定する（図２／ｓ０３４）。
【００６９】
ωh（ｔid）
＝−α∫ｄｔ・（［４ε（ｔid）2−｛ωh（ｔ）−ωm（ｔid）｝2］1/2
×ｓｉｎ［ｓｉｎ-1｛（ωh（ｔ）−ωm（ｔid-1））／２ε（ｔid）｝
−δθ1（ｔid）］）  ‥（２．５）
ここで、α＝（αL，αR）の各成分は系の安定性を表す係数である。
【００７０】
続いて、第２角速度設定部２３５が、左右各成分について、第１角速度設定部２３４により設定された仮想運動振動子θhの角速度ωhに基づき、仮想誘導振動子θmの角速度ωmを新たな固有角速度ωMとして設定する（図２／ｓ０３５）。具体的には、第２角速度設定部２３５が、左右各成分について、第２位相差δθ2が目標位相差δθdに近づくように、次式（２．６）にしたがって仮想誘導振動子θmの角速度ωm＝（ωmL，ωmR）を設定する。
【００７１】
ωm（ｔid）
＝β∫ｄｔ・（［４ε（ｔid）2−｛ωh（ｔid）−ωm（ｔ）｝2）
×ｓｉｎ［ｓｉｎ-1｛（ωh（ｔid）−ωm（ｔ））／２ε（ｔid）｝−δθd］）
‥（２．６）
ここで、β＝（βL，βR）の各成分は系の安定性を表す係数である。
【００７２】
続いて、２次振動子生成部２４０が、運動振動子測定部２１０により測定された股関節角度φHと、固有角速度設定部２３０により設定された新たな固有角速度ωMとに基づき、「第２モデル」にしたがって、２次振動子ｙ＝（ｙL+，ｙL-，ｙR+，ｙR-）を生成する（図２／ｓ０４０）。第２モデルは、股関節角度φH等の運動振動子に応じて出力が変動する仮想的な複数の神経要素等、仮想的な複数の第２要素の関係を表現するものである。より具体的には、第２モデルは、左大腿部の屈曲方向（前方）および伸展方向（後方）のそれぞれへの運動を支配する神経要素L+およびL-、並びに右大腿部の屈曲方向および伸展方向のそれぞれへの運動を支配する神経要素R+およびR-の膜電位の変動に対応する状態変数ｕi（i=L+,L-,R+,R-）と、神経要素ｉの順応効果が反映される自己抑制因子ｖiとを含む、次の連立微分方程式（３）によって表現される。
【００７３】
τ1L+・ｄｕL+／ｄｔ＝−ｕL+＋ｗL+/L-ｙL-＋ｗL+/R+ｙR+−λLｖL+
＋ｆ1（ωML）＋ｆ2（ωML）Ｋ（φL），
τ1L-・ｄｕL-／ｄｔ＝−ｕL-＋ｗL-/L+ｙL+＋ｗL-/R-ｙR-−λLｖL-
＋ｆ1（ωML）＋ｆ2（ωML）Ｋ（φL），
τ1R+・ｄｕR+／ｄｔ＝−ｕR+＋ｗR+/L+ｙL+＋ｗR+/R-ｙR-−λRｖR+
＋ｆ1（ωMR）＋ｆ2（ωMR）Ｋ（φR），
τ1R-・ｄｕR-／ｄｔ＝−ｕR-＋ｗR-/L-ｙL-＋ｗR-/R+ｙR+−λRｖR-
＋ｆ1（ωMR）＋ｆ2（ωMR）Ｋ（φR），
τi・ｄｖi／ｄｔ＝−ｖi＋ｙi，
ｙi＝Ｈ（ｕi−ｕth），Ｈ≡０（ｕi＜ｕth），１（ｕi≧ｕth）または
ｙi＝ｆs（ｕi）≡１／｛１＋ｅｘｐ（−ｕi／ξ）｝  ‥（３）
ここで、τ1iは状態変数ｕiの変化特性を規定する時定数であり、左右各成分について、次式（３．１）によって表されているように、新たな固有角速度ωMへの依存性を有する。
【００７４】
τ1i≡ｔ（ωML）／ωML−γL（i=L+,L-），
ｔ（ωMR）／ωMR−γR（i=R+,R-）‥（３．１）
ｔ（ω）はω依存性を有する係数である。γ＝（γL，γR）は定数である。
【００７５】
また、τ2iは自己抑制因子ｖiの変化特性を規定する時定数である。さらに、ｗi/j（＜０）は仮想的な複数の第２要素（神経要素）ｉおよびｊの関係を表す第２相関係数（定数）である。「λL」および「λR」は慣れ係数である。Ｋは股関節角度φHに応じたフィードバック係数である。
【００７６】
「ｆ1」および「ｆ2」はそれぞれ次式（３．２）および（３．３）により定義される関数である。
【００７７】
ｆ1（ω）≡ｃ・ω（ｃ＞０）  ‥（３．２）
ｆ2（ω）≡ｃ0＋ｃ1ω＋ｃ2ω2  ‥（３．３）
新たな固有角速度ωMの関数であるｆ1（ωM）およびｆ2（ωM）の係数ｃ，ｃ0，ｃ1，ｃ2は、目標運動設定部２１１によって設定された目標となる運動リズムに応じた係数として設定され得る。
【００７８】
なお、股関節角度φHに代えてまたは加えて、股関節角速度ｄφH／ｄｔや、膝関節、足関節、肩関節、肘関節の角度や角速度、さらには歩行者の着地音、呼吸音、意図的な発声音等、歩行運動リズムと連関したリズムで変動する種々の振動子に基づき、２次振動子ｙiが生成されてもよい。
【００７９】
２次振動子ｙiは、状態変数ｕiの値が閾値ｕth未満である場合は０、状態変数ｕiの値が閾値ｕth以上である場合はこのｕiの値をとる。或いは、２次振動子ｙiは、シグモイド関数ｆsによって定義されている（式（３）参照）。これにより、大腿部の屈曲方向（前方）への動きについてはこの動きを支配する第２要素（神経要素）L+、R+の出力である２次振動子ｙL+およびｙR+が、他の第２要素の出力よりも大きくなる。また、大腿部の伸展方向（後方）への動きについてはこの動きを支配する第２素子L-、R-の出力である２次振動子ｙL-およびｙR-が、他の第２要素の出力よりも大きくなる。脚体（大腿部）の前方または後方への動きは、たとえば、股関節角速度ｄφH／ｄｔの極性によって識別される。
【００８０】
前記のように股関節角速度ｄφH／ｄｔ等の運動振動子に応じて出力が変動する仮想的な複数の第２要素の関係を表現する第２モデルにしたがって、第２要素の出力として２次振動子ｙiが生成される（式（３），図２／ｓ０４０）。これにより、ユーザの実際の運動に関係する現実の複数の要素の関係を第２モデルにおける仮想的な複数の第２要素の係数ｗi/jにより表される相関関係に反映させることで、当該現実の複数の要素間の関係に鑑みて適当な２次振動子ｙiが生成され得る。たとえば、現実の複数の要素としてユーザの複数の神経（ニューロン）が想定された場合、左右の脚体による歩行を支配するニューロン間の定性的関係等が反映された形で２次振動子ｙiが生成される。したがって、ユーザの運動を誘導する誘導振動子のリズムおよびスケールを当該関係に鑑みて適当なものとすることができる。
【００８１】
次に、誘導振動子生成部２５０が運動振動子測定部２１０によって測定された股関節角度φHおよび股関節角速度ｄφH／ｄｔと、２次振動子生成部２４０により生成された２次振動子ｙiと、固有角速度設定部２３０によって設定された新たな固有角速度ωMとに基づき、誘導振動子ｚを設定する（図３／ｓ０５０）。
【００８２】
具体的には、次式（４）にしたがって第１誘導振動子ｚ1が生成される（図３／ｓ０５１）。
【００８３】
ｚ1L＝ｇ1+（ωmL）ｇ+（φHL）ｙL+−ｇ1-（ωmL）ｇ-（φHL）ｙL-，
ｚ1R＝ｇ1+（ωmR）ｇ+（φHR）ｙR+−ｇ1-（ωmR）ｇ-（φHR）ｙＲR-  ‥（４）
ここで「ｇ1+」「ｇ1-」「ｇ+」および「ｇ-」は次式（４．１）〜（４．４）のそれぞれによって定義される関数である。
【００８４】
ｇ1+（ｐ，ω）≡Σkａk+ωk  （ａk+：係数，ｋ＝０〜３）‥（４．１）
ｇ1-（ｐ，ω）≡Σkａk-ωk  （ａk-：係数，ｋ＝０〜３）‥（４．２）
ｇ+（φ）≡ｃ1+（φ−φ0+）＋ｃ2+（φ−φ0+）3
（ｃ1+，ｃ2;：係数，φ0+：屈曲方向の股関節角度φHの目標値）‥（４．３）
ｇ-（φ）≡ｃ1-（φ−φ0-）＋ｃ2-（φ−φ0-）3
（ｃ1-，ｃ2-：係数，φ0-：伸展方向の股関節角度φHの目標値）‥（４．４）
【００８５】
第１誘導振動子ｚ1は、第１係数ｇ1+およびｇ1-をそれぞれバネ係数（弾性係数）とする、図４に示されている２つの仮想的なバネＧ1+およびＧ1-の弾性力として認識される。第１係数ｇ1+およびｇ1-は、固有角速度ωMに応じて、股関節角度（ユーザの運動スケールに応じた運動振動子）φHを目標運動スケールに応じた目標角度φ0+（＞０）およびφ0-（＜０）に近付ける第１ポテンシャル（仮想的なバネ（弾性要素）のポテンシャル）のグラディエントを表す（式（４．１）（４．２）参照）。すなわち、第１誘導振動子ｚ1は、第１ポテンシャルに応じた第１係数ｇ1+，ｇ1-を弾性係数（バネ係数）とし、かつ、股関節角度φHの値を目標角度φ0+およびφ0-に復元させる仮想的なバネ等の弾性要素による弾性力として表現される。これにより、筋肉の収縮状態から伸展状態への移行時の弾性力等、ユーザの身体の弾性要素が反映されたリズムおよびスケールをもってユーザの運動が誘導され得る。
【００８６】
仮想的なバネＧ1+による弾性力は、そのバネ係数ｇ1+に応じて股関節角度φHをその目標角度φ0+に近付けるようにユーザの大腿部に作用する（式（４）参照）。すなわち、股関節角度φHが目標角度φ0+未満である場合、バネＧ1+による弾性力が、股関節角度φHを増加させる方向（前方）に大腿部を動かすようにこの大腿部に作用する。また、股関節角度φHが目標角度φ0+を超えた場合、バネＧ1+による弾性力が、股関節角度φHを減少させる方向（後方）に大腿部を動かすようにこの大腿部に作用する。
【００８７】
また、他の仮想的なバネＧ1-による弾性力は、そのバネ係数ｇ1-に応じて、股関節角度φHをこの目標角度φ0-に近付けるようにユーザの大腿部に作用する（式（４）参照）。すなわち、股関節角度φHが目標角度φ0-を超えている場合、バネＧ1-による弾性力が、股関節角度φHを減少させる方向（後方）に大腿部を動かすようにこの大腿部に作用する。また、股関節角度φHが目標角度φ0-を下回った場合、バネＧ1-による弾性力が、股関節角度φHを増加させる方向（前方）に大腿部を動かすようにこの大腿部に作用する。
【００８８】
また、前記のように大腿部の前方への動きおよび後方への動きの別に応じて、複数の第２要素ｉ（=L+,L-,R+,R-）のうち一部から偏重的に出力があるので、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。
【００８９】
すなわち、左の大腿部が前方に動いているとき、この動きを支配する第２要素L+に応じた２次振動子ｙL+の値が他の第２要素L-に応じた２次振動子ｙL-の値より大きくなって、式（４）により表されている第１誘導振動ｚ1Lが次式（４ａ）のように表される（または近似される）。
【００９０】
ｚ1L＝ｇ1+（ωmL）ｇ+（φHL）ｙL+  ‥（４ａ）
【００９１】
このため、左の大腿部が前方に動いているとき、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力のうち、股関節角度φHを前側の目標角度φ0+に近付けるようにユーザの大腿部に作用するバネＧ1+の弾性力が偏重的に反映されることとなる。これにより、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。
【００９２】
また、左の大腿部が後方に動いているとき、この動きを支配する第２要素L-の出力が他の第２要素L+の出力より大きくなり、これによって第２要素L-に応じた２次振動子ｙL-の値が第２要素L+に応じた２次振動子ｙL+の値より大きくなって、式（４）によって表されている第１誘導振動子ｚ1Lは次式（４ｂ）のように表される（または近似される）。
【００９３】
ｚ1L＝−ｇ1-（ωmL）ｇ-（φHL）ｙL-  ‥（４ｂ）
【００９４】
このため、左の大腿部が後方に動いているとき、２つの仮想のバネＧ1+およびＧ１−のそれぞれの弾性力のうち、股関節角度φHを後側の目標角度φ0-に近付けるようにユーザの大腿部に作用する仮想のバネＧ1-の弾性力が偏重的に反映されることとなる。これにより、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。これは、右の脚体（大腿部）の動きについても同様である。
【００９５】
なお、股関節角速度ｄφH／ｄｔを変数とするシグモイド関数ｆs（式（３）参照）が第１係数ｇ1+，ｇ1-に組み込まれ、これにより股関節角速度ｄφH／ｄｔの極性により特定される大腿部の前後への動きの別に応じて、複数の第２要素ｉの出力としての２次振動子ｙiのうち一部が偏重的に反映された形で第１トルクＴ1が生成されてもよい。これによっても、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避され得る。
【００９６】
さらに、次式（５）にしたがって第２誘導振動子ｚ2が設定される（図２／ｓ０５２）。
【００９７】
ｚ2L＝−ｇ2+（ωmL）（ｄφHL／ｄｔ）Ｈ+（φHL）ｙL+
＋ｇ2-（ωmL）（ｄφHL／ｄｔ）Ｈ-（φHL）ｙL-，
ｚ2R＝−ｇ2+（ωmR）（ｄφHR／ｄｔ）Ｈ+（φHR）ｙR+
＋ｇ2-（ωmR）（ｄφHR／ｄｔ）Ｈ-（φHR）ｙR-‥（５）
ここで「ｇ2+」「ｇ2-」「Ｈ+」および「Ｈ-」は次式（５．１）〜（５．４）のそれぞれによって定義される関数である。
【００９８】
ｇ2+（ω）≡Σkｂk+ωk  （ｂK+：係数，ｋ＝０〜３）‥（５．１）
ｇ2-（ω）≡Σkｂk-ωk  （ｂk-：係数，ｋ＝０〜３）‥（５．２）
Ｈ+（φ）≡０（φ≦０），１（φ＞０）  ‥（５．３）
Ｈ-（φ）≡０（φ＞０），１（φ≦０）  ‥（５．４）
【００９９】
第２誘導振動子ｚ2は、第２係数ｇ2+およびｇ2-をそれぞれダンパ係数（減衰係数）とする、図４に示されている２つの仮想的なダンパＧ2;およびＧ2-の減衰力として把握される。第２係数ｇ2+およびｇ2-は、固有角速度ωMに応じて股関節角度φHの絶対値の増大を抑制する第２ポテンシャル（仮想的なダンパ（減衰要素）のポテンシャル）のグラディエントを表す（式（５．１）（５．２）参照）。すなわち、第２誘導振動子ｚ2は、第２ポテンシャルに応じた第２係数ｇ2+，ｇ2-を減数係数（ダンパ係数）とし、かつ、股関節角速度ｄφH／ｄｔに応じて股関節角度φHの絶対値の増大を抑制する仮想的なダンパ等の減衰要素による減衰力として表現される。これにより、筋肉の伸展状態から屈曲状態への移行時の粘性力等、ユーザの身体の減衰要素が反映されたリズムおよびスケールをもってユーザの運動が誘導され得る。
【０１００】
一方の仮想的なダンパＧ2+による減衰力は、そのダンパ係数ｇ2+および股関節角速度ｄφH／ｄｔに応じて、前側（屈曲側）への股関節角度φＨの絶対値の増大を抑制するようにユーザの大腿部に作用する（式（５）参照）。すなわち、仮想的なダンパＧ2+による減衰力は、大腿部の前方への過剰な動きを抑制するようにこの大腿部に作用する。
【０１０１】
また、他方の仮想的なダンパＧ2-による弾性力は、そのダンパ係数ｇ2-および股関節角速度ｄφH／ｄｔに応じて、後側（伸展側）への股関節角度φHの絶対値の増大を抑制するようにユーザの大腿部に作用する（式（５）参照）。すなわち、仮想的なダンパＧ2-による減衰力は、大腿部の後方への過剰な動きを抑制するようにこの大腿部に作用する。
【０１０２】
また、第２誘導振動子ｚ2には、股関節角度φHの関数としての階段関数Ｈ+，Ｈ-が含まれている。したがって、２つの仮想的なダンパＧ2+およびＧ2-のそれぞれの減衰力が相殺される事態が回避される。
【０１０３】
そして、誘導振動子生成部２５０により生成された第１誘導振動子ｚ1＝（ｚ1L，ｚ1R）と、第２誘導振動子ｚ2＝（ｚ2L，ｚ2R） とを含む誘導振動子ｚ（＝ｚ1＋ｚ2）に応じた電流Ｉ＝（ＩL，ＩR）が電池２０６から左右のアクチュエータ２１０にそれぞれ供給され、ユーザの大腿部に力（股関節回りのトルク）Ｔが作用する。
【０１０４】
以後、前記処理（図２／ｓ０１１，ｓ０１２，・・，ｓ０４０，図３／ｓ０５０）が繰り返されることで、ユーザは歩行運動誘導装置２０の動作によって股関節回りのトルクＴが作用している状態で歩行運動する。
【０１０５】
そして、ユーザが前記のようにトレッドミル１０および歩行運動誘導装置２０の動作により歩行運動が誘導されている状態で、運動変数測定部１１０が歩行比ｋを測定する（図３／ｓ１１０）。具体的には、第１運動変数測定部１０１がトレッドミル１０のベルト１３の速度に応じた速度センサ１０２の出力に基づき、ユーザの歩行速度ｖを測定する（図３／ｓ１１１）。また、第２運動変数測定部１１２がトレッドミル１０の踏板１４が受ける圧力に応じた圧力センサ１０４の出力がピークを示す単位時間当たりの回数に基づき、ユーザの歩行率（単位時間当たりの歩数）ｐを測定する（図２／Ｓ１１２）。なお、歩行運動誘導装置２０によりユーザに作用させられる股関節回りのトルクＴの時間変化に基づいて歩行率ｐが測定されてもよい。また、ユーザの身体に加速度センサが取り付けられ、ユーザの鉛直方向の加速度に応じた加速度センサの出力に基づいて歩行率ｐが測定されてもよい。そして、運動変数測定部１００が歩行速度ｖおよび歩行率ｐの関数としての歩行比ｋ＝ｖ／ｐ2を測定する。
【０１０６】
また、生理変数測定部１２０が生理変数センサ１０６の出力に基づいて整理変数ｂを測定する（図３／ｓ１２０）。生理変数ｂには、ユーザの酸素摂取量、呼吸周波数、心拍数、脈拍数、血圧、血中飽和酸素濃度、乳酸値、筋電位、および消費エネルギーのうち一部または全部が含まれる。
【０１０７】
生理変数ｂが測定されるたび、股関節角度φHの目標角度φ0+およびφ0-、第１係数ｇ1+（ωM）およびｇ1-（ωM）のそれぞれに含まれる係数ａk+およびａk-、ならびに第２係数ｇ2+（ωM）およびｇ2-（ωM）のそれぞれに含まれる係数ｂk+およびｂk-のうち一部または全部が変更設定される。目標角度φ0+およびφ0-が変更設定されることで歩行運動誘導装置２０による誘導スケールが変化し、ユーザの歩幅ｑが変更されうる。また、係数ａk+，ａk-，ｂk+およびｂk-の一部または全部が変更設定されることで歩行運動誘導装置２０による誘導リズムが変化し、ユーザの歩行率（目標運動リズム）ｐが変更されうる。これにより、歩行運動誘導装置２０によって歩行比ｋが逐次変化するようにユーザの歩行運動が誘導されていく過程で、生理変数ｂが逐次測定される。そして、歩行比ｋおよび生理変数ｂの関係が求められ、推奨歩行比設定部１３０がこの関係に基づいて推奨歩行比ｋ0を設定する（図３／ｓ１３０）。
【０１０８】
たとえば、生理変数ｂとしてユーザの消費エネルギーが測定され、図５に丸（○）で示されているような歩行比ｋおよび消費エネルギー（生理変数）の関係が得られた場合を考える。この場合、当該関係は図５に示されている歩行比ｋの２次曲線によって近似表現され、この２次曲線が最小値をとる歩行比ｋ＝０．０７５がこのユーザにとっての推奨歩行比ｋ0として設定される。ユーザの消費エネルギーが最小となるように設定された推奨歩行比ｋ0は、ユーザの身体負荷の軽減に鑑みて適当な歩行比ｋである。なお、心拍数、乳酸値（身体の疲労度を表す。）、筋群活性度（身体機能の活性化の程度を表す。）等、消費エネルギーとは異なる生理変数ｂが測定され、当該測定変数ｂに基づいて推奨歩行比ｋ0が設定されてもよい。
【０１０９】
次に、歩行運動誘導装置２０による推奨歩行比ｋ0に応じたユーザの歩行運動の誘導方法について図６を用いて説明する。推奨歩行比ｋ0に応じたユーザの歩行運動の誘導方法は、推奨歩行比ｋ0の設定時におけるユーザの歩行運動の誘導方法と、誘導振動子ｚの生成方法（図３／ｓ０５０）において相違するのみである。そこで、運動振動子の測定（図２／ｓ０１１，ｓ０１２）、１次振動子ｘの生成（図２／ｓ０２０）、固有角速度ωMの設定（図２／ｓ０３０）および２次振動子ｙの生成（図２／ｓ０４０）については説明を省略する。
【０１１０】
第１運動変数測定部１１１が、トレッドミル１０の動作および歩行運動誘導装置２０の動作により歩行運動が誘導されているユーザの歩行速度ｖを測定する（図６／ｓ１１１）。
【０１１１】
また、誘導振動子生成部２５０が、運動振動子測定部２１０により測定された股関節角度φHおよび股関節角速度ｄφH／ｄｔと、２次振動子生成部２４０により生成された２次振動子ｙiと、固有角速度設定部２３０によって設定された固有角速度ωMとに加え、さらに、第１運動変数測定部１１１により測定されたユーザの歩行速度ｖと、推奨歩行比設定部１３０により測定された推奨歩行比ｋ0とに基づき、誘導振動子ｚを設定する（図２／ｓ２５０）。
【０１１２】
具体的には、次式（６）にしたがって第１誘導振動子ｚ1が生成される（図６／ｓ２５１）。
【０１１３】
ｚ1L＝ｇ1+（ｐ0，ωmL）ｇ+（ｑ0，φHL）ｙL+
−ｇ1-（ｐ0，ωmL）ｇ-（ｑ0，φHL）ｙL-，
ｚ1R＝ｇ1+（ｐ0，ωmR）ｇ+（ｑ0，φHR）ｙR+
−ｇ1-（ｐ0，ωmR）ｇ-（ｑ0，φHR）ｙＲR-，
ｐ0≡（ｖ／ｋ0）1/2，ｑ0≡（ｖｋ0）1/2‥（６）
ここで「ｇ1+」「ｇ1-」「ｇ+」および「ｇ-」は次式（６．１）〜（６．４）のそれぞれによって定義される関数である。
【０１１４】
ｇ1+（ｐ0，ω）≡Σkａk+（ｐ0）ωk
（ａk+（ｐ0）：係数，ｋ＝０〜３）‥（６．１）
ｇ1-（ｐ0，ω）≡Σkａk-（ｐ0）ωk
（ａk-（ｐ0）：係数，ｋ＝０〜３）‥（６．２）
ｇ+（ｑ0，φ）≡ｃ1+（φ−φ0+（ｑ0））＋ｃ2+（φ−φ0+（ｑ0））3
（ｃ1+，ｃ2;：係数，φ0+：屈曲方向の股関節角度φHの目標値）‥（６．３）
ｇ-（ｑ0，φ）≡ｃ1-（φ−φ0-（ｑ0））＋ｃ2-（φ−φ0-（ｑ0））3
（ｃ1-，ｃ2-：係数，φ0-：伸展方向の股関節角度φHの目標値）‥（６．４）
【０１１５】
第１誘導振動子ｚ1により、股関節角度φHが推奨歩行比ｋ0および歩行速度ｖに応じた目標角度φ0+，φ0-に一致するようにユーザの歩行運動が誘導されうる。
【０１１６】
さらに、次式（７）にしたがって第２誘導振動子ｚ2が設定される（図６／ｓ２５２）。
【０１１７】
ｚ2L＝−ｇ2+（ｐ0，ωmL）（ｄφHL／ｄｔ）Ｈ+（φHL）ｙL+
＋ｇ2-（ｐ0，ωmL）（ｄφHL／ｄｔ）Ｈ-（φHL）ｙL-，
ｚ2R＝−ｇ2+（ｐ0，ωmR）（ｄφHR／ｄｔ）Ｈ+（φHR）ｙR+
＋ｇ2-（ｐ0，ωmR）（ｄφHR／ｄｔ）Ｈ-（φHR）ｙR-‥（７）
ここで「ｇ2+」「ｇ2-」は次式（７．１）〜（７．２）のそれぞれによって定義される関数である。
【０１１８】
ｇ2+（ｐ0，ω）≡Σkｂk+（ｐ0）ωk
（ｂK+（ｐ0）：係数，ｋ＝０〜３）‥（５．１）
ｇ2-（ｐ0，ω）≡Σkｂk-（ｐ0）ωk
（ｂk-（ｐ0）：係数，ｋ＝０〜３）‥（５．２）
【０１１９】
第２誘導振動子ｚ2により、股関節角速度ｄφH／ｄｔに応じて股関節角度φHの絶対値の増大が抑制されるようにユーザの歩行運動が誘導されうる。
【０１２０】
股関節角度φHの目標角度φ0+（ｑ0）およびφ0-（ｑ0）は推奨歩行比ｋ0および歩行速度ｖに応じた推奨歩幅ｑ0（＝（ｖｋ0）1/2）関数であり、ユーザの歩幅ｑ（＝ｖ／ｐ）と推奨歩幅ｑ0との偏差δｑに基づいて補正されうる。
【０１２１】
また、第１係数ｇ1+（ｐ0，ωM），ｇ1-（ｐ0，ωM）に含まれる係数ａk+（ｐ0）およびａk-（ｐ0）、第２係数ｇ2+（ｐ0，ωM），ｇ2-（ｐ0，ωM）に含まれる係数ｂk+（ｐ0）およびｂk-（ｐ0）は推奨歩行比ｋ0および歩行速度ｖに応じた推奨歩行率ｐ0（＝（ｖ／ｋ0）1/2）の関数であり、ユーザの歩行率ｐ（＝ｖ／ｑ）と推奨補効率ｐ0との偏差δｐに基づいて補正されうる。
【０１２２】
そして、誘導振動子生成部２５０により生成された第１誘導振動子ｚ1と、第２誘導振動子ｚ2とを含む誘導振動子ｚに応じた電流Ｉが電池２０６から左右のアクチュエータ２１０にそれぞれ供給され、ユーザの大腿部に力（股関節回りのトルク）Ｔが作用する。
【０１２３】
以後、前記処理（図２／ｓ０１１，ｓ０１２，・・，ｓ０４０，図６／ｓ２５０）が繰り返されることで、ユーザは歩行運動誘導装置２０の動作によって股関節回りのトルクＴが作用している状態で歩行する。
【０１２４】
前記機能を発揮する運動管理システム１によれば、歩行運動誘導装置（第２運動誘導装置）２０の動作によりユーザの運動が誘導されている過程で歩行比（運動変数）ｋおよび生理変数ｂが測定される（図３／ｓ１１０，ｓ１２０）。そして、この動物の「身体機能の活性化」または「身体負荷の軽減」に鑑みて生理変数ｂが適当な値となった状態における歩行比ｋの値が推奨目標比ｋ0として設定される。これにより、推奨歩行比ｋ0が、このユーザの身体機能の活性化等に鑑みて適当に設定されうる。そして、歩行運動誘導装置２０によりこの推奨歩行比ｋ0に基づいてユーザの歩行運動が誘導されることで、ユーザの歩幅（運動スケール）ｑおよび歩行率（運動リズム）ｐがこのユーザの身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【０１２５】
運動管理システム１の効果に関する実験結果について図７〜図９を用いて説明する。
【０１２６】
図７に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の歩行比ｋが異なる。すなわち、ユーザが歩行補助装置２０を装着せずに歩行した場合の歩行比ｋは約０．００６５に制御されている。これに対して、ユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合の歩行比ｋは、前記のように設定された推奨歩行比ｋ（＝０．００７５）に制御されている。このユーザは、歩行運動誘導装置２０の動作によって歩行比ｋ、ひいては歩幅ｑが大きくなるように歩行運動が誘導されている。
【０１２７】
また、図８に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の心拍数（生理変数）が異なる。すなわち、ユーザが歩行補助装置２０を装着せずに歩行した場合の心拍数は、ユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合の心拍数よりも高い。これは、図７に示されているように、ユーザの歩幅が大きくなるように歩行運動が誘導されている状態で、その心拍数が低く抑えられており、身体への負荷が軽減されていることを意味する。
【０１２８】
また、図９に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の筋群活性度（生理変数）が異なる。図９（ｂ）に示されているユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合のその股関節周囲の筋群活性度は、図９（ａ）に示されているユーザが歩行補助装置２０を装着せずに歩行した場合の同じ箇所の筋群活性度よりも高い。これは、図７に示されているように、ユーザは歩幅が大きくなるように歩行運動が誘導されることで、ユーザの身体機能が活性化されていることを意味する。
【０１２９】
図７〜図９に示されている実験結果から、たとえば、運動管理システム１が高齢者等、歩行運動機能が低下しているユーザの運動管理に利用された場合の意義が大きいことがわかる。すなわち、本発明の運動管理システム１が利用されることで、高齢者等の身体負荷を軽減しながら、その身体機能の低下抑制、さらには活性化を促すことができる。
【０１３０】
また、ユーザの歩行速度（第１運動変数）ｖがトレッドミル１０のベルト１３の速さに基づいて測定される（図３／ｓ１１１）。またユーザの歩行率（第２運動変数）ｐが、ユーザの着床および離床の繰り返しにより変化するトレッドミル１０への圧力に基づいて測定される（図３／ｓ１１２）。これにより、歩行速度ｖおよび歩行率ｐの関数である歩行比ｋの測定精度の向上が図られる（図３／ｓ１１０）。したがって、ユーザの身体機能の活性化や身体負荷の軽減に鑑みて適当な運動変数の推奨値がより正確に設定されうる。
【０１３１】
また、トレッドミル１０の無端ベルト（回動体）１３を回動させることにより、この無端ベルト１３の動きとは逆向きの歩行または走行運動が誘導される。これにより、当該歩行または走行に伴うユーザの移動が無端ベルト１３の回動によって相殺されるので、トレッドミル１０の設置スペースさえあればユーザの歩行または走行運動に伴う推奨歩行比ｋ0の設定が可能である。トレッドミル１０が省略され、歩行運動誘導装置２０のみによってユーザの歩行運動が誘導されている状態で推奨歩行比ｋ0が設定されてもよい。
【０１３２】
なお、運動管理システム１は、ユーザの歩行運動以外のあらゆる運動に伴って変化する運動変数の推奨値の設定に利用されてもよい。たとえば、車椅子の左右の車輪に手で力を加える運動について、両腕が車輪を押すリズム等、他の運動についての運動変数の推奨値設定に利用されてもよい。運動管理システム１は馬の走行運動等、人間以外の動物の運動についての運動変数の推奨値設定に利用されてもよい。また、歩行比ｋのほか、歩行速度ｖおよび歩行率ｐの任意の関数としての運動変数の推奨値が設定されてもよい。
【０１３３】
また、トレッドミル１０のほか、複数のローラに掛け渡された無端ベルト、中心または中心から外れた点を通る軸回りに回動される球体、楕円球体、中心軸または中心軸から離れた中心軸に平行な軸回りに回動する円筒、四角筒等の筒体、その他任意の軸回りに回動される一塊の物体等の回動体を回動させることでこの回動体に接触するユーザの腕や脚などの身体部分の運動を誘導する装置が第１運動誘導装置として採用されてもよい。
【０１３４】
前記実施形態では誘導振動子ｚに応じた左右の股関節回りのトルクＴ＝（ＴL，ＴR）がユーザの身体に作用させられたが、他の実施形態として膝関節、足関節、肩関節、肘関節、手根関節等、種々の関節回りのトルクがユーザの身体に作用させされてもよい。トルク作用対象となる関節の組合せは、ユーザに応じてさまざまに変更されてもよい。
【０１３５】
また、第２制御部２００の２次振動子生成部２４０が、運動振動子測定部２１０により測定された股関節角度φHまたは股関節角速度ｄφH／ｄｔ（運動振動子）の周期変化の大きさが閾値を超えた場合、固有角速度設定部２３０により設定された新たな固有角速度ωMに代えて、運動振動子測定部２１０により測定された股関節角度φH等の角速度（位相の時間変化）と、１次振動子生成部２２０により生成された１次振動子ｘの角速度とのうち一方または両方が反映されたリズムで振動する２次振動子ｙを生成してもよい。
【０１３６】
当該構成によれば、ユーザの運動リズムが急に変化した場合でも、この変化後の運動リズムに応じた適当なリズムをもってユーザの運動が誘導されうる。
【０１３７】
また、特開２００４−７３６４９号公報に開示されている方法にしたがって誘導振動子ｚが生成されてもよい。
【０１３８】
さらに、次のような形で誘導振動子ｚが生成されてもよい。
【０１３９】
すなわち、運動振動子測定部２１０がユーザの２つの異なる身体部分の動きに応じた運動振動子（たとえば肩関節角速度および股関節角度）のそれぞれを「第１運動振動子」および「第２運動振動子」として測定する。また、１次振動子生成部２２０が第１運動振動子と、固有角速度ωMが反映された形で相互に引き込み合う１次振動子ｘを生成する。さらに、固有角速度設定部２３０が第１運動振動子と１次振動子ｘとの位相差に基づき、新たな固有角速度ωMを設定する。また、２次振動子生成部２４０が第２運動振動子に基づき、新たな固有角速度ωMが反映されたリズムで振動する２次振動子ｙを生成する。そして、誘導振動子生成部２５０が２次振動子ｙに加え、第１運動変数測定部１０１により測定された歩行速度ｖ（第１運動変数）と、第２運動変数測定部１０２により測定された歩行率ｐ（第２運動変数）とのうち一方または両方に基づいて誘導振動子ｚを生成する。
【０１４０】
当該構成によれば、ユーザの異なる身体部分のそれぞれの運動リズムと、この運動を誘導するリズムとの調和を図りながら、ユーザの運動リズムが目標運動リズムに近づくようにこの運動が誘導されうる。
【０１４１】
さらに、本発明の運動管理システム１は、推奨歩行比設定部１３０が複数のユーザのそれぞれの身体的特徴（身体のサイズ、ウェイト、年齢、性別等が含まれる。）を識別するための識別子を認識した上で各ユーザの複数の推奨歩行比ｋ0を予備的に設定し、認識した識別子ごとに予備的に設定済みの推奨歩行比ｋ0に基づいて新たな推奨歩行比ｋ0を設定することにより、識別子および推奨値が対応付けられているデータベースを構築してもよい。この場合、ユーザの識別子、ひいては身体的特徴に応じてデータベースから検索された標準的な推奨値に基づいて、ユーザの歩行運動が誘導される。これにより、ユーザの歩行比ｋがこのユーザの身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【図面の簡単な説明】
【０１４２】
【図１】本発明の運動管理システムの構成例示図
【図２】本発明の運動管理方法の説明図
【図３】本発明の運動管理方法の説明図
【図４】歩行運動を誘導する仮想的なバネおよびダンパの説明図
【図５】推奨歩行比の設定方法の説明図
【図６】本発明の運動管理方法説明図
【図７】本発明の運動管理システムの効果に関する実験結果の説明図
【図８】本発明の運動管理システムの効果に関する実験結果の説明図
【図９】本発明の運動管理システムの効果に関する実験結果の説明図
【符号の説明】
【０１４３】
１‥運動管理システム、１０‥トレッドミル（第１運動誘導装置）、１３‥回動体、２０‥歩行運動誘導装置（第２運動誘導装置）、２５‥アクチュエータ、２６‥股関節角度センサ、１００‥第１制御部、１１０‥運動変数測定部、１１１‥第１運動変数測定部、１１２‥第２運動変数測定部、１２０‥生理変数測定部、１３０‥推奨歩行比設定部（推奨値設定部）、２００‥第２制御部（制御部）、２１０‥運動振動子測定部、２２０‥１次振動子生成部、２３０‥固有角速度設定部、２３１‥第１位相差設定部、２３２‥第２位相差設定部、２３３‥相関係数設定部、２３４‥第１角速度設定部、２３５‥第２角速度設定部、２４０‥２次振動子生成部、２５０‥誘導振動子生成部、２５１‥第１誘導振動子生成部、２５２‥第２誘導振動子生成部
【図１】
2007275283.tif 000003
【図２】
2007275283.tif 000004
【図３】
2007275283.tif 000005
【図４】
2007275283.tif 000006
【図５】
2007275283.tif 000007
【図６】
2007275283.tif 000008
【図７】
2007275283.tif 000009
【図８】
2007275283.tif 000010
【図９】
2007275283.tif 000011
【手続補正書】
【提出日】平成19年6月22日(2007.6.22)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
動物の運動を管理するシステムであって、
前記動物の運動リズムと調和したリズムで該動物に力を作用させることによって該動物の運動スケールを調節しながら該動物の運動を誘導する第２運動誘導装置の動作を制御する制御部と、
前記動物の生理状態に応じた生理変数を測定する生理変数測定部と、
前記動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定部と、
前記制御部により前記運動変数が変化するように制御されている前記第２運動誘導装置の動作によって前記動物の運動が誘導されている過程で、前記生理変数測定部により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、前記運動変数測定部により測定された前記運動変数を該運動変数の推奨値として設定する推奨値設定部とを備え、
前記運動変数測定部が、
前記動物の運動の速さを調節しながら前記第２運動誘導装置とは別に該動物の運動を誘導する第１運動誘導装置の動作により運動が誘導されている前記動物の運動の速さを表す前記運動変数としての第１運動変数を該第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定部と、
前記動物の運動のリズムまたはスケールを表す前記運動振動子としての第２運動変数を該第１運動誘導装置と該動物との相互作用状態または該第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定部とにより構成され、かつ、
前記第１運動変数および前記第２運動変数の関数を前記運動変数として測定することを特徴とする運動管理システム。
【請求項２】
請求項１記載の運動管理システムにおいて、
前記制御部が、
前記動物の身体運動に応じて時間変化するパラメータとして、該動物の第１運動振動子および第２運動振動子を測定する運動振動子測定部と、
入力振動信号と相互に引き込み合うことで固有角速度に基づいて定まる角速度で時間変化する出力振動信号を生成する第１モデルに、前記運動振動子測定部により測定された前記第１運動振動子を該入力振動信号として入力することにより、該出力振動信号として１次振動子を生成する１次振動子生成部と、
該運動振動子測定部により測定された該第１運動振動子と該１次振動子生成部により生成された該１次振動子との位相差である第１位相差を目標位相差に近づけるように該固有角速度を新たに設定する固有角速度設定部と、
入力振動信号に基づき、該固有角速度設定部により設定された該固有角速度に基づいて定まる角速度で時間変化する出力振動信号を生成する第２モデルに、該運動振動子測定部により測定された該第２運動振動子を該入力振動信号として入力することにより、該出力振動信号として２次振動子を生成する２次振動子生成部と、
該２次振動子生成部により生成された該２次振動子に基づいて前記第２運動誘導装置の動作のスケールおよびリズムを特定する誘導振動子を生成する誘導振動子生成部とを備えていることを特徴とする運動管理システム。
【請求項３】
請求項１記載の運動管理システムにおいて、
前記第１運動変数測定部が前記動物の歩行または走行速度を第１運動変数として測定し、第２運動変数測定部が動物の歩幅または歩行率（単位時間当たりの歩数を意味する。）を第２運動変数として測定することを特徴とする運動管理システム。
【請求項４】
請求項３記載の運動管理システムにおいて、
前記運動変数測定部が前記第１運動変数としての歩行または走行速度に対する第２運動変数としての歩幅の二乗の比率、または第前記２運動変数としての歩行率の二乗に対する前記第１運動変数としての歩行または走行速度の比率である歩行比を前記運動変数として測定することを特徴とする運動管理システム。
【請求項５】
請求項１記載の運動管理システムにおいて、
前記第１運動変数測定部が、前記第１運動誘導装置が有する回動体の回動方向に逆らう方向に運動が誘導されている前記動物の前記第１運動変数を、当該回動体の回動速度に基づいて測定することを特徴とする運動管理システム。
【請求項６】
請求項１記載の運動管理システムにおいて、
前記第１運動変数測定部が、前記第１運動誘導装置としてのトレッドミルが有する複数のローラに掛け渡された回動体としての無端ベルトの動きに逆らう方向に誘導されている前記動物の歩行または走行運動の速さを表す前記第１運動変数としての歩行または走行速度を当該無端ベルトの駆動速度に基づいて測定することを特徴とする運動管理システム。
【請求項７】
請求項１記載の運動管理システムにおいて、
前記推奨値設定部が、複数の前記動物のそれぞれの身体的特徴を識別するための識別子を認識した上で各動物の複数の前記推奨値を予備的に設定し、認識した該識別子ごとに予備的に設定済みの該推奨値に基づいて新たな該推奨値を設定することにより、該識別子および該推奨値が対応付けられているデータベースを構築することを特徴とする運動管理システム。
【請求項８】
動物の運動を管理する方法であって、
前記動物の運動リズムと調和したリズムで該動物に力を作用させることによって該動物の運動スケールを調節しながら該動物の運動を誘導する第２運動誘導装置の動作を制御する制御処理と、
前記動物の生理状態に応じた生理変数を測定する生理変数測定処理と、
前記動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定処理と、
前記制御処理において前記運動変数が変化するように制御されている前記第２運動誘導装置の動作によって前記動物の運動が誘導されている過程で、前記生理変数測定処理において測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、前記運動変数測定処理において測定された前記運動変数を該運動変数の推奨値として設定する推奨値設定処理とを実行し、
前記運動変数測定処理において、
前記動物の運動の速さを調節しながら前記第２運動誘導装置とは別に該動物の運動を誘導する第１運動誘導装置の動作により運動が誘導されている前記動物の運動の速さを表す前記運動変数としての第１運動変数を該第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定処理と、
前記動物の運動のリズムまたはスケールを表す前記運動振動子としての第２運動変数を該第１運動誘導装置と該動物との相互作用状態または該第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定処理とを実行し、かつ、
前記第１運動変数および前記第２運動変数の関数を前記運動変数として測定することを特徴とする運動管理方法。
【請求項９】
動物の運動を管理する装置としてコンピュータを機能させるプログラムであって、
前記動物の運動リズムと調和したリズムで該動物に力を作用させることによって該動物の運動スケールを調節しながら該動物の運動を誘導する第２運動誘導装置の動作を制御する制御部と、
前記動物の生理状態に応じた生理変数を測定する生理変数測定部と、
前記動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定部と、
前記制御部により前記運動変数が変化するように制御されている前記第２運動誘導装置の動作によって前記動物の運動が誘導されている過程で、前記生理変数測定部により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、前記運動変数測定部により測定された前記運動変数を該運動変数の推奨値として設定する推奨値設定部とを備え、
前記運動変数測定部が、
前記動物の運動の速さを調節しながら前記第２運動誘導装置とは別に該動物の運動を誘導する第１運動誘導装置の動作により運動が誘導されている前記動物の運動の速さを表す前記運動変数としての第１運動変数を該第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定部と、
前記動物の運動のリズムまたはスケールを表す前記運動振動子としての第２運動変数を該第１運動誘導装置と該動物との相互作用状態または該第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定部とにより構成され、かつ、
前記第１運動変数および前記第２運動変数の関数を前記運動変数として測定する前記装置として前記コンピュータを機能させることを特徴とする運動管理プログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動物の運動を管理するシステム、方法および当該管理機能をコンピュータに付与するプログラムに関する。
【背景技術】
【０００２】
ユーザの歩行訓練のレベルを設定する等の目的のため、このユーザの心電信号等の生理指標に鑑みて適当な歩行速度を決定する手法が提案されている（たとえば、特許文献１参照）。
【特許文献１】特開２００３−１５４０２９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、ユーザの身体に力を作用させることで脚等の動きを誘導（または補助）する装置がユーザに装着されている場合、このユーザにとって適当な誘導スケールおよびリズムを設定することは考慮されていなかった。
【０００４】
そこで、本発明は、人間等の動物の生理状態に応じて、その運動を適当なスケールおよびリズムで誘導しうるように当該動物の運動を管理するシステム、方法、および当該管理機能をコンピュータに付与するプログラムを提供することを解決課題とする。
【課題を解決するための手段】
【０００５】
第１発明の運動管理システムは、動物の運動を管理するシステムであって、前記動物の運動リズムと調和したリズムで該動物に力を作用させることによって該動物の運動スケールを調節しながら該動物の運動を誘導する第２運動誘導装置の動作を制御する制御部と、前記動物の生理状態に応じた生理変数を測定する生理変数測定部と、前記動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定部と、前記制御部により前記運動変数が変化するように制御されている前記第２運動誘導装置の動作によって前記動物の運動が誘導されている過程で、前記生理変数測定部により測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、前記運動変数測定部により測定された前記運動変数を該運動変数の推奨値として設定する推奨値設定部とを備え、前記運動変数測定部が、前記動物の運動の速さを調節しながら前記第２運動誘導装置とは別に該動物の運動を誘導する第１運動誘導装置の動作により運動が誘導されている前記動物の運動の速さを表す前記運動変数としての第１運動変数を該第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定部と、前記動物の運動のリズムまたはスケールを表す前記運動振動子としての第２運動変数を該第１運動誘導装置と該動物との相互作用状態または該第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定部とにより構成され、かつ、前記第１運動変数および前記第２運動変数の関数を前記運動変数として測定することを特徴とする。
【０００６】
第１発明の運動管理システムによれば、第２運動誘導装置から力を受けながら動物の運動が誘導されている過程で「生理変数」が測定される。動物の酸素摂取量、呼吸周波数、心拍数、脈拍数、血圧、血中飽和酸素濃度、乳酸値、筋電位、および消費エネルギーのうち一部または全部が「生理変数」として測定されうる。また、「第１運動変数」が第１運動誘導装置の動作の速さに基づいて測定され、「第２運動変数」が第１運動誘導装置と動物との相互作用または第２運動誘導装置による誘導スケールもしくはリズムに基づいて測定される。これにより、第１および第２運動変数、ひいてはこれらの関数である運動変数の測定精度の向上が図られる。さらに、この動物の「身体機能の活性化」または「身体負荷の軽減」に鑑みて生理変数が適当な値となった状態における運動変数の値が当該運動変数の推奨値として設定される。これにより、動物の運動変数の推奨値が、この動物の身体機能の活性化等に鑑みて適当に設定されうる。したがって、第２運動誘導装置によりこの推奨値に基づく力が加えられながら動物の運動が誘導されることで、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうる。
【０００７】
第２発明の運動管理システムは、第１発明の運動管理システムにおいて、前記制御部が、前記動物の身体運動に応じて時間変化するパラメータとして、該動物の第１運動振動子および第２運動振動子を測定する運動振動子測定部と、入力振動信号と相互に引き込み合うことで固有角速度に基づいて定まる角速度で時間変化する出力振動信号を生成する第１モデルに、前記運動振動子測定部により測定された前記第１運動振動子を該入力振動信号として入力することにより、該出力振動信号として１次振動子を生成する１次振動子生成部と、該運動振動子測定部により測定された該第１運動振動子と該１次振動子生成部により生成された該１次振動子との位相差である第１位相差を目標位相差に近づけるように該固有角速度を新たに設定する固有角速度設定部と、入力振動信号に基づき、該固有角速度設定部により設定された該固有角速度に基づいて定まる角速度で時間変化する出力振動信号を生成する第２モデルに、該運動振動子測定部により測定された該第２運動振動子を該入力振動信号として入力することにより、該出力振動信号として２次振動子を生成する２次振動子生成部と、該２次振動子生成部により生成された該２次振動子に基づいて前記第２運動誘導装置の動作のスケールおよびリズムを特定する誘導振動子を生成する誘導振動子生成部とを備えていることを特徴とする。
【０００８】
第２発明の運動管理システムによれば、動物の運動リズムと、第２運動誘導装置の動作リズム（動物の運動を誘導するリズム）との相互の歩み寄りまたは調和が図られながら、当該動物の運動が当該第２運動誘導装置により誘導される。また、このように動物の運動が誘導されている過程で「生理変数」および「運動変数」が測定される。したがって、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の運動および第２運動誘導装置の動作の調和を図りながら、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうる。なお、ここで「振動」とは、現実的または仮想的なもの（物体）がほぼ一定の周期をもって揺れ動くことのほか、広く時間変化することを意味する。また「振動子」とは、値が時間変化する電気信号や、ソフトウェアにおいて値が時間変化するものとして定義される関数等を意味する。
【０００９】
第３発明の運動管理システムは、第１発明の運動管理システムにおいて、前記第１運動変数測定部が前記動物の歩行または走行速度を第１運動変数として測定し、第２運動変数測定部が動物の歩幅または歩行率（単位時間当たりの歩数を意味する。）を第２運動変数として測定することを特徴とする運動管理。
【００１０】
第３発明の運動管理システムによれば、動物の歩行または走行の速さ（速度）および歩幅（スケール）または歩行率（リズム）の関数である運動変数の推奨値が、この動物の身体機能の活性化や身体負荷の軽減に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の歩幅および歩行率がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００１１】
第４発明の運動管理システムは、第３発明の運動管理システムにおいて、前記運動変数測定部が前記第１運動変数としての歩行または走行速度に対する第２運動変数としての歩幅の二乗の比率、または第前記２運動変数としての歩行率の二乗に対する前記第１運動変数としての歩行または走行速度の比率である歩行比を前記運動変数として測定することを特徴とする。
【００１２】
第４発明の運動管理システムによれば、動物の歩行比の推奨値が、この動物の身体機能の活性化や身体負荷の軽減に鑑みて適当に設定されうる。そして、第２運動誘導装置によりこの推奨値に基づいて動物の運動が誘導されることで、動物の歩行比がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００１３】
第５発明の運動管理システムは、第１発明の運動管理システムにおいて、前記第１運動変数測定部が、前記第１運動誘導装置が有する回動体の回動方向に逆らう方向に運動が誘導されている前記動物の前記第１運動変数を、当該回動体の回動速度に基づいて測定することを特徴とする。
【００１４】
第５発明の運動管理システムによれば、第１運動誘導装置が有する回動体を回動させることにより、回動体と身体部分が接触する動物の、回動体の回動方向とは逆向きの運動が誘導される。「回動体」には、複数のローラに掛け渡された無端ベルト、中心または中心から外れた点を通る軸回りに回動される球体、楕円球体、中心軸または中心軸から離れた中心軸に平行な軸回りに回動する円筒、四角筒等の筒体、その他任意の軸回りに回動される一塊の物体が含まれる。これにより、当該運動に伴う動物またはその身体部分の移動が回動体の回動によって相殺されるので、第１運動誘導装置の設置スペースさえあれば動物の運動変数の推奨値設定が可能であり、当該設定スペースが比較的自由に選択されうる。
【００１５】
第６発明の運動管理システムは、第１発明の運動管理システムにおいて、前記第１運動変数測定部が、前記第１運動誘導装置としてのトレッドミルが有する複数のローラに掛け渡された回動体としての無端ベルトの動きに逆らう方向に誘導されている前記動物の歩行または走行運動の速さを表す前記第１運動変数としての歩行または走行速度を当該無端ベルトの駆動速度に基づいて測定することを特徴とする。
【００１６】
第６発明の運動管理システムによれば、トレッドミル（第１運動誘導装置）の無端ベルト（回動体）を回動させることにより、この無端ベルトの動きとは逆向きの歩行または走行が誘導される。これにより、当該歩行または走行に伴う動物の移動が無端ベルトの回動によって相殺されるので、第１運動誘導装置の設置スペースさえあれば動物の歩行または走行運動に伴う運動変数の推奨値設定が可能であり、当該設定スペースが比較的自由に選択されうる。
【００１７】
第７発明の運動管理システムは、第１発明の運動管理システムにおいて、前記推奨値設定部が、複数の前記動物のそれぞれの身体的特徴を識別するための識別子を認識した上で各動物の複数の前記推奨値を予備的に設定し、認識した該識別子ごとに予備的に設定済みの該推奨値に基づいて新たな該推奨値を設定することにより、該識別子および該推奨値が対応付けられているデータベースを構築することを特徴とする。
【００１８】
第７発明の運動管理システムによれば、個々の動物の身体的特徴（身体のサイズ、ウェイト、年齢、性別等が含まれる。）に応じて標準的な推奨値が設定される。そして、２運動誘導装置によりこの標準的な推奨値に基づいて動物の運動が誘導されることで、動物の運動変数がこの動物の身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【００１９】
第８発明の運動管理方法は、動物の運動を管理する方法であって、前記動物の運動リズムと調和したリズムで該動物に力を作用させることによって該動物の運動スケールを調節しながら該動物の運動を誘導する第２運動誘導装置の動作を制御する制御処理と、前記動物の生理状態に応じた生理変数を測定する生理変数測定処理と、前記動物の運動スケールおよびリズムのうち一方または両方に応じた運動変数を測定する運動変数測定処理と、前記制御処理において前記運動変数が変化するように制御されている前記第２運動誘導装置の動作によって前記動物の運動が誘導されている過程で、前記生理変数測定処理において測定された生理変数が該動物の身体機能の活性化または身体負荷の軽減に鑑みて適当な値となった状態で、前記運動変数測定処理において測定された前記運動変数を該運動変数の推奨値として設定する推奨値設定処理とを実行し、前記運動変数測定処理において、前記動物の運動の速さを調節しながら前記第２運動誘導装置とは別に該動物の運動を誘導する第１運動誘導装置の動作により運動が誘導されている前記動物の運動の速さを表す前記運動変数としての第１運動変数を該第１運動誘導装置の動作の速さに基づいて測定する第１運動変数測定処理と、前記動物の運動のリズムまたはスケールを表す前記運動振動子としての第２運動変数を該第１運動誘導装置と該動物との相互作用状態または該第２運動誘導装置の動作状態に基づいて測定する第２運動変数測定処理とを実行し、かつ、前記第１運動変数および前記第２運動変数の関数を前記運動変数として測定することを特徴とする。
【００２０】
第８発明の運動管理方法によれば、第２運動誘導装置から力を受けながら動物の運動が誘導されている過程でこの動物の「身体機能の活性化」または「身体負荷の軽減」に鑑みて生理変数が適当な値となった状態における運動変数の値が当該運動変数の推奨値として設定される。これにより、動物の運動変数の推奨値が、この動物の身体機能の活性化等に鑑みて適当に設定されうる。したがって、第２運動誘導装置によりこの推奨値に基づく力が加えられながら動物の運動が誘導されることで、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうる。
【００２１】
第９発明の運動管理プログラムは、第１発明の運動管理装置としてコンピュータを機能させることを特徴とする。
【００２２】
第９発明の運動管理プログラムによれば、第２運動誘導装置によりこの推奨値に基づく力が加えられながら動物の運動が誘導されることで、動物の生理状態に応じてこの動物の運動がその身体機能の活性化等に鑑みて適当なスケールおよびリズムで誘導されうるように動物の運動を管理する装置としてコンピュータを機能させることができる。
【発明を実施するための最良の形態】
【００２３】
本発明の運動管理システム、運動管理方法および運動管理プログラムの実施形態について図面を用いて説明する。以下、歩行者の脚体等について左右を区別するためにパラメータに添字Ｌ、Ｒを添付するが、表記の簡単のため左右を区別する必要が特にない場合には添字Ｌ、Ｒを省略する。
【００２４】
本発明の運動管理システムの構成について図１を用いて説明する。
【００２５】
図１に示されている運動管理システム１は、トレッドミル（第１運動誘導装置）１０および歩行運動誘導装置（第２運動誘導装置）を用いて、ユーザ（人間（動物））の身体機能の活性化等に鑑みて適当な推奨歩行比を設定するためのものである。また。運動管理システム１は、ユーザの歩行（または走行）運動と、トレッドミル１０の動作と、歩行運動誘導装置２０の動作との調和を図りながら、推奨歩行比に応じたスケールおよびリズムで歩行するようにこのユーザを訓練するためのものである。
【００２６】
トレッドミル１０は、人間の標準的な横幅よりも若干広い程度の幅の駆動ローラ１１と、駆動ローラ１１とほぼ同じ幅の従動ローラ１２と、駆動ローラ１１および従動ローラ１２に掛け渡されている無端ベルト（回動体）１３と、無端ベルト１３のうちユーザが乗る部分を下から支える踏板１４とを備えている。駆動ローラ１１は、モータ、変速機等より構成される駆動機構１０１によって駆動される。駆動ローラ１１が図中時計回りに駆動されることにより、従動ローラ１２の同方格への従動を伴ってベルト１３も時計回りに回る。これによりベルト１３に乗っているユーザの図中左への歩行（または走行）が誘導される。また、ベルト１３の移動速度に応じた信号を出力する速度センサ１０２と、踏板１４が受ける圧力に応じた信号を出力する圧力センサ１０４とが設けられている。ユーザの身体には、その心拍数等の生理変数ｂに応じた信号を出力する生理変数センサ１０６が取り付けられている。生理変数センサ１０６は生理変数ｂの内容によって適当な方法でユーザの身体に取り付けられる。なお、トレッドミル１０としては、一般に市販されているもの等、公知のあらゆる構成のトレッドミルが採用されてもよい。
【００２７】
歩行運動誘導装置２０は、腰部装具２１と、大腿部装具２２と、力伝達部材２３と、バッテリ２４と、アクチュエータ（電動モータ）２５と、股関節角度センサ２６とを備えている。
【００２８】
腰部装具２１は剛性のある素材と柔軟性のある素材とが組み合わせられて作られており、ユーザの腰部に装着される。大腿部装具２２も剛性のある素材と柔軟性のある素材とが組み合わせられて作られており、ユーザの大腿部の前後それぞれに装着される。力伝達部材２３は、軽量の硬質プラスチック等の定形性のある素材より作られており、ユーザの大腿部に沿って、ユーザの腰部の横から下方に延びた後で大腿部の前後に向けて二股に分かれた形状であり、アクチュエータ２５および前後の大腿部装具２２のそれぞれに連結されている。バッテリ２４は腰部装具２１に収納されており（たとえば、腰部装具２１を構成する複数枚の素材の間に固定されており）、アクチュエータ２５等に対して電力を供給する。アクチュエータ２５は腰部装具２１に取り付けられており、力伝達部材２３および大腿部装具２２を介してユーザの大腿部に力を作用させる。股関節角度センサ２６はユーザの腰部の横に設けられたロータリエンコーダ等により構成され、股関節角度に応じた信号を出力する。
【００２９】
運動管理システム１は、第１制御部１００と、第２制御部２００とを備えている。
【００３０】
第１制御部１００は、トレッドミル１０に付属するマイクロコンピュータ等のコンピュータにより構成されており、駆動機構１０２による駆動ローラ１１の駆動速度等を制御する。第１制御部１００は、運動変数測定部１１０と、生理変数測定部１２０と、推奨歩行比設定部１３０とを備えている。運動変数測定部１１０、生理変数測定部１２０および推奨歩行比設定部１３０はそれぞれハードウェアとしてのコンピュータ（ＣＰＵ、ＲＯＭ、ＲＡＭ、Ｉ／Ｏ（入出力装置）等により構成されている。）と、このコンピュータに機能を付与するソフトウェアとしての本発明の「運動管理プログラム」の一部とにより構成されている。
【００３１】
運動変数測定部１１０は第１運動変数測定部１１１と、第２運動変数測定部１１２とを備えている。
【００３２】
第１運動変数測定部１１１は、速度センサ１０２の出力に基づいてユーザの歩行速度ｖを当該ユーザの運動の速さを表す第１運動変数として測定する。第２運動変数測定部１１２は、圧力センサ１０４の出力に基づいてユーザの歩行率（単位時間当たりの歩数）ｐを当該ユーザの運動のリズムを表す第２運動変数として測定する。そして、運動変数測定部１１０は、歩行速度ｖおよび歩行率ｐの関数である運動変数として歩行比ｋ（＝ｖ／ｐ2）を測定する。
【００３３】
生理変数測定部１２０は、ユーザの身体に取り付けられた生理変数センサ１０６の出力に基づいてこのユーザの生理状態に応じた生理変数ｂを測定する。
【００３４】
推奨歩行比設定部１３０は、運動変数測定部１１０により測定されたユーザの歩行比ｋと、生理変数測定部１２０により測定された生理変数ｂとに基づき、推奨歩行比ｋ0を設定する。
【００３５】
第２制御部（本発明の「制御部」に該当する。）２００は、歩行運動誘導装置２０の腰部装具２１に収納されたコンピュータと、このコンピュータを歩行運動誘導装置２０の当該制御部として機能させるソフトウェアとしての本発明の「運動管理プログラム」とにより構成されている。
【００３６】
第２制御部２００は、運動振動子測定部２１０と、１次振動子生成部２２０と、固有角速度設定部２３０と、２次振動子生成部２４０と、誘導振動子生成部２５０とを備えている。
【００３７】
運動振動子測定部２１０は、股関節角度センサ２６の出力に基づき、第２運動振動子として股関節角度φH（ユーザの運動スケールに応じた運動振動子である「第３運動振動子」にも該当する。）を測定する。また、運動振動子測定部２１０は、股関節角度センサ２６の出力に基づき、第１運動振動子として股関節角速度ｄφH／ｄｔ（第３運動振動子の時間微分値である「第４運動振動子」にも該当する。）を測定する。
【００３８】
１次振動子生成部２２０は運動振動子測定部２１０により測定された第１運動振動子としての股関節角速度ｄφH／ｄｔと、固有角速度ωMとに基づき、「第１モデル」にしたがって１次振動子ｘを生成する。
【００３９】
固有角速度設定部２３０は、第１位相差設定部２３１と、第２位相差設定部２３２と、相関係数設定部２３３と、第１角速度設定部２３４と、第２角速度設定部２３５とを備えている。
【００４０】
第１位相差設定部２３１は股関節角速度ｄφH／ｄｔの角速度ωHと、ファン・デル・ポル方程式に含まれる固有角速度ωMに基づいて定まる角速度で時間変化する１次振動子ｘとの位相差を第１位相差δθ1として設定する。第２位相差設定部２３２は仮想運動振動子θhと仮想誘導振動子（仮想誘導振動子）θmとの関係を表す「仮想モデル」にしたがって、仮想運動振動子θhと仮想誘導振動子θmとの位相差を第２位相差δθ2（＝θh−θm）として設定する。
【００４１】
相関係数設定部２３３は、第２位相差設定部２３２により設定された第２位相差δθ2が、第１位相差設定部２３１により設定された第１位相差δθ1に近づくように仮想運動振動子θhと仮想誘導振動子θmとの相関係数εを設定する。第１角速度設定部２３４は、相関係数設定部２３３により設定された相関係数εに基づき、仮想運動振動子θhの角速度ωhを設定する。第２角速度設定部２３５は、第１角速度設定部２３４により設定された仮想運動振動子θhの角速度ωhに基づき、第２位相差設定部２３２により設定された第２位相差δθ2が、目標位相差設定部２１２により設定された目標位相差δθｄに近づくように仮想誘導振動子θmの角速度ωmを新たな固有角速度ωMとして設定する。
【００４２】
２次振動子生成部２４０は、運動振動子測定部２１０により測定された第２運動振動子としての股関節角度φHと、固有角速度設定部２３０により設定された新たな固有角速度ωMとに基づき「第２モデル」にしたがって２次振動子ｙを生成する。
【００４３】
誘導振動子生成部２５０は、第１誘導振動子生成部２５１と、第２誘導振動子生成部２５２とを備えている。
【００４４】
第１誘導振動子生成部２５１は２次振動子生成部２４０により生成された２次振動子ｙと、固有角速度設定部２３０により設定された固有角速度ωMと、運動振動子測定部２１０により測定された股関節角度φH（第３運動振動子）をその目標値φ0に近付けるための仮想的男性要素の弾性力を表す第１誘導振動子ｚ1を生成する。第２誘導振動子生成部２５２は２次振動子生成部２４０により生成された２次振動子ｙと、固有角速度設定部２３０により設定された固有角速度ωMと、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔ（第４運動振動子）に応じて股関節角度φH（第３運動振動子）の絶対値の増大を抑制するための仮想的な減衰要素の減衰力を表す第２誘導振動子ｚ2を生成する。そして、誘導振動子生成部２５０は、第１誘導振動子ｚ1および第２誘導振動子ｚ2に基づき、歩行運動誘導装置２０によりユーザに作用させられる股関節回りのトルクＴを特定する誘導振動子ｚを生成する。
【００４５】
第１制御部１００および第２制御部２００は無線通信が可能とされている。なお、第１制御部１００に代えて第２制御部２００が推奨歩行比設定部１３０を備えている等、運動管理システムを構成する複数の処理部１１０（１１１，１１２）１，１２０，１３０，２１０，２２０，‥が、第１制御部１００および第２制御部２００に任意のパターンで配設されていてもよい。また、第１制御部１００および第２制御部２００は同一のコンピュータにより構成されていてもよい。第１制御部１００および第２制御部２００の通信は有線通信であってもよい。
【００４６】
前記構成の運動管理システム１の機能、特にユーザの推奨歩行比（運動変数の推奨値）の設定方法について図２〜図５を用いて説明する。
【００４７】
第１制御部１００により、ベルト１３の速度が一定になるようにトレッドミル１０の動作が制御されることで、ベルト１３の動きとは逆方向へのほぼ一定速度でのユーザの歩行運動が誘導される。
【００４８】
この状態で、運動振動子測定部２１０は股関節角度センサ２６の出力に基づき、ユーザの左右の股関節角度φH＝（φHL，φHR）（第２運動振動子および第３運動振動子に該当する。）を測定する（図２／ｓ０１１）。また、運動振動子測定部２１０は股関節角度センサ２６の出力に基づき、ユーザの左右の股関節角速度ｄφH／ｄｔ＝（ｄφHL／ｄｔ，ｄφHR／ｄｔ）（第１運動振動子および第４運動振動子に相当する。）を測定する（図２／ｓ０１２）。
【００４９】
さらに、１次振動子生成部２２０が、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔおよびメモリに記憶されている最新の固有角速度ωM＝（ωML，ωMR）に基づき「第１モデル」にしたがって１次振動子ｘ＝（ｘL，ｘR）を設定する（図２／ｓ０２０）。第１モデルは、股関節角速度ｄφH／ｄｔ等の第１運動振動子に応じて出力が変動する、左右の脚体等、複数の第１要素の相関関係を表現するモデルである。具体的には「第１モデル」は入力振動信号と相互に引き込み合うことで、固有角速度に基づいて定まる角速度で時間変化する出力振動信号（すなわち、当該第１要素のそれぞれの出力振動信号）を生成するモデルである。第１モデルは、たとえば次式（１）で表されるファン・デル・ポル（ｖａｎ  ｄｅｒ  Ｐｏｌ）方程式によって表現される。１次振動子生成部２２０によって股関節角速度ｄφH／ｄｔ（第１運動振動子）が当該入力振動信号として第１モデルに入力されることにより、当該出力振動信号として１次振動子ｘ＝（ｘL，ｘR）が生成される。
【００５０】
（ｄ2ｘL／ｄｔ2）＝ξ（１−ｘL2）（ｄｘL／ｄｔ）−ωML2ｘL
＋ｇ（ｘL−ｘR）＋Ｋ（ｄφHL／ｄｔ），
（ｄ2ｘR／ｄｔ2）＝ξ（１−ｘR2）（ｄｘR／ｄｔ）−ωMR2ｘR
＋ｇ（ｘR−ｘL）＋Ｋ（ｄφHR／ｄｔ）  ‥（１）
【００５１】
ここで「ξ」は１次振動子ｘおよびその１回時間微分（ｄｘ／ｄｔ）がｘ−（ｄｘ／ｄｔ）平面で安定なリミットサイクルを描くように設定される係数（＞０）である。ｇは第１モデルにおいて左右の脚体（第１要素）の相関関係を表す第１相関係数である。Ｋはフィードバック係数である。なお、固有角速度ωMは、歩行運動誘導装置２０による実際の歩行補助リズム（歩行誘導リズム）から大きく外れない範囲で任意に設定されてよい。
【００５２】
１次振動子ｘ＝（ｘL，ｘR）はルンゲ・クッタ法により算定（設定）される。１次振動子ｘの成分ｘLおよびｘRはそれぞれ左右の脚体に関する仮想的な歩行補助リズムを表す。また、振動子ｘはファン・デル・ポル方程式の１つの性質である「相互引き込み」により、実際の歩行運動リズムとほぼ同じリズム（角速度）で時間変化する股関節角速度ｄφH／ｄｔのリズムと調和しながらも「固有角速度」ωMに基づいて定まる自律的なリズムまたは角速度をもって時間変化または振動するという性質がある。
【００５３】
なお、股関節角速度ｄφH／ｄｔに代えてまたは加えて、股関節角度φHや、膝関節、足関節、肩関節、肘関節の角度や角速度、さらには歩行者の着地音、呼吸音、断続的な発声音等、ユーザの歩行運動リズム（運動リズム）が反映されたリズムで変動する種々の第１運動振動子に基づき、１次振動子ｘが生成されてもよい。
【００５４】
また、式（１）で表現されるファン・デル・ポル方程式とは異なる形のファン・デル・ポル方程式によって第１モデルが表現されてもよく、股関節角速度ｄφH／ｄｔ等の入力振動信号と相互引き込み効果をもって出力振動信号が生成され得るあらゆる方程式によって第１モデルが表現されてもよい。
【００５５】
前記のように股関節角速度ｄφH／ｄｔ等の第１運動振動子に応じて出力が変動する仮想的な複数の第１要素（左右の脚体）の関係を表現する第１モデルにしたがって、第１要素の出力として１次振動子ｘが生成される（式（１），図２／ｓ０２０）。これにより、ユーザの実際の運動に関係する複数の第１要素の関係を、第１相関係数ｇ等により表される第１モデルにおける複数の第１要素の相関関係に反映させることで、当該複数の第１要素間の関係に鑑みて適当な１次振動子ｘが生成され得る。たとえば、複数の第１要素として左右の脚体や同一脚体の複数の関節が想定された場合、交互に前後に動く等の左右の脚体の定性的関係や股関節回りの脚体運動と膝関節回りの脚体運動との周期や位相差等の同一脚体の関節間の定性的関係等が反映された形で１次振動子ｘが生成される。したがって、ユーザの運動を誘導する誘導振動子のリズムおよびスケールを当該関係に鑑みて適当なものとすることができる。
【００５６】
続いて、固有角速度設定部２３０がメモリに記憶されている目標位相差δθdと、１次振動子設定部２１０により生成された１次振動子ｘとに基づき、仮想的な２つの振動子が含まれている仮想モデルにしたがって固有角速度ωMを新たに設定する（図２／ｓ０３０）。
【００５７】
具体的には、まず、第１位相差測定部２３１が、左右各成分について、運動振動子測定部２１０により測定された股関節角速度ｄφH／ｄｔ（第１運動振動子）の位相θHと、１次振動子設定部２１０により生成された１次振動子ｘの位相θMとの位相差θH−θMを、第１位相差δθ1として設定する（図２／ｓ０３１）。
【００５８】
次に、第２位相差設定部２３２が、過去３歩行周期にわたって第１位相差設定部２３１により設定された第１位相差δθ1が一定であったことを要件として、左右各成分について、次式（２．１）および（２．２）によって表される「仮想モデル」にしたがって、次式（２．３）によって表される仮想運動振動子θhと仮想誘導振動子θmとの位相差θh−θmを第２位相差δθ2として設定する（図２／ｓ０３２）。
【００５９】
ｄθh／ｄｔ＝ωh＋ε・ｓｉｎ（θmL−θhL）  ‥（２．１）
ｄθm／ｄｔ＝ωm＋ε・ｓｉｎ（θhL−θmL）  ‥（２．２）
δθ2＝ａｒｃｓｉｎ［（ωh−ωm）／２ε］  ‥（２．３）
【００６０】
ここで、ε＝（εL，εR）は仮想モデルにおける仮想運動振動子θh＝（θhL，θhR）および仮想誘導振動子θm＝（θmL，θmR）の左右成分ごとの相関係数である。また、ωhは仮想運動振動子θhの角速度であり、ωmは仮想誘導振動子θmの角速度である。
【００６１】
続いて、相関係数設定部２３３が、第１位相差設定部２３１により設定された第１位相差δθ1と、第２位相差設定部２３２により設定された第２位相差δθ2との差δθ1−δθ2が最小になるように、相関係数εを設定する（図２／ｓ０３３）。
【００６２】
具体的には次式（２．４）にしたがって、左右各成分について、股関節角速度ｄφH／ｄｔ（第１運動振動子）が０となる離散的な時間ｔid（ｄ＝１，２，・・）における相関係数εが逐次設定される。
【００６３】
ε（ｔid+1）
＝ε（ｔid）−η｛Ｖ1（ｔid+1）−Ｖ1（ｔid）｝／｛ε（ｔid）−ε（ｔid-1）｝，
Ｖ（ｔid+1）≡（１／２）｛δθ1（ｔid+1）−δθ2（ｔid）｝2  ‥（２．４）
【００６４】
ここで、η＝（ηL，ηR）の各成分は、第１位相差δθ1の左右各成分と、第２位相差δθ2の左右各成分とを近づけるポテンシャルＶ＝（ＶL，ＶR）の安定性を表す係数である。
【００６５】
次に、第１角速度設定部２３４が、相関係数設定部２３３により設定された相関係数εに基づき、仮想誘導振動子θmの固有角速度ωmが一定であるという条件下で、左右各成分について第１および第２位相差の差δθ1−δθ2の各成分が最小となるように仮想運動振動子θhの角速度ωhを次式（２．５）にしたがって設定する（図２／ｓ０３４）。
【００６６】
ωh（ｔid）＝−α∫ｄｔ・（［４ε（ｔid）2−｛ωh（ｔ）−ωm（ｔid）｝2］1/2
×ｓｉｎ［ｓｉｎ-1｛（ωh（ｔ）−ωm（ｔid-1））／２ε（ｔid）｝−δθ1（ｔid）］）  ‥（２．５）
【００６７】
ここで、α＝（αL，αR）の各成分は系の安定性を表す係数である。
【００６８】
続いて、第２角速度設定部２３５が、左右各成分について、第１角速度設定部２３４により設定された仮想運動振動子θhの角速度ωhに基づき、仮想誘導振動子θmの角速度ωmを新たな固有角速度ωMとして設定する（図２／ｓ０３５）。具体的には、第２角速度設定部２３５が、左右各成分について、第２位相差δθ2が目標位相差δθdに近づくように、次式（２．６）にしたがって仮想誘導振動子θmの角速度ωm＝（ωmL，ωmR）を設定する。
【００６９】
ωm（ｔid）＝β∫ｄｔ・（［４ε（ｔid）2−｛ωh（ｔid）−ωm（ｔ）｝2）
×ｓｉｎ［ｓｉｎ-1｛（ωh（ｔid）−ωm（ｔ））／２ε（ｔid）｝−δθd］）
‥（２．６）
【００７０】
ここで、β＝（βL，βR）の各成分は系の安定性を表す係数である。
【００７１】
続いて、２次振動子生成部２４０が、運動振動子測定部２１０により測定された股関節角度φH（第２運動振動子）と、固有角速度設定部２３０により設定された新たな固有角速度ωMとに基づき、第２モデルにしたがって、２次振動子ｙ＝（ｙL+，ｙL-，ｙR+，ｙR-）を生成する（図２／ｓ０４０）。「第２モデル」は、股関節角度φH等の第２運動振動子に応じて出力が変動する複数の神経要素等、複数の第２要素の相関関係を表現するモデルである。具体的には「第２モデル」は、入力振動信号に応じて当該新たな固有角速度に基づいて定まる角速度で時間変化する出力振動信号（すなわち、当該第２要素のそれぞれの出力振動信号）を生成するモデルである。第２モデルは次式（３）の連立微分方程式により表現される。当該連立微分方程式には、左大腿部の屈曲方向（前方）および伸展方向（後方）のそれぞれへの運動を支配する神経要素L+およびL-、並びに右大腿部の屈曲方向および伸展方向のそれぞれへの運動を支配する神経要素R+およびR-の膜電位の変動に対応する状態変数ｕi（i=L+,L-,R+,R-）と、神経要素ｉの順応効果が反映される自己抑制因子ｖiとが含まれている。２次振動子生成部２４０により、股関節角度φH（第２運動振動子）が当該入力振動信号として第２モデルに入力されることにより、当該出力振動信号として２次振動子ｙが生成される。
【００７２】
τ1L+・ｄｕL+／ｄｔ＝−ｕL+＋ｗL+/L-ｙL-＋ｗL+/R+ｙR+−λLｖL+
＋ｆ1（ωML）＋ｆ2（ωML）Ｋ（φL），
τ1L-・ｄｕL-／ｄｔ＝−ｕL-＋ｗL-/L+ｙL+＋ｗL-/R-ｙR-−λLｖL-
＋ｆ1（ωML）＋ｆ2（ωML）Ｋ（φL），
τ1R+・ｄｕR+／ｄｔ＝−ｕR+＋ｗR+/L+ｙL+＋ｗR+/R-ｙR-−λRｖR+
＋ｆ1（ωMR）＋ｆ2（ωMR）Ｋ（φR），
τ1R-・ｄｕR-／ｄｔ＝−ｕR-＋ｗR-/L-ｙL-＋ｗR-/R+ｙR+−λRｖR-
＋ｆ1（ωMR）＋ｆ2（ωＭMR）Ｋ（φR），
τ2i・ｄｖi／ｄｔ＝−ｖi＋ｙi，
ｙi＝Ｈ（ｕi−ｕth），Ｈ≡０（ｕi＜ｕth），１（ｕi≧ｕth）または
ｙi＝ｆs（ｕi）≡１／｛１＋ｅｘｐ（−ｕi／ξ）｝  ‥（３）
【００７３】
ここで、τ1iは状態変数ｕiの変化特性を規定する時定数であり、左右各成分について、次式（３．１）によって表されているように、新たな固有角速度ωMへの依存性を有する。
【００７４】
τ1i≡ｔ（ωML）／ωML−γL（i=L+,L-），
ｔ（ωMR）／ωMR−γR（i=R+,R-）‥（３．１）
【００７５】
ｔ（ω）はω依存性を有する係数である。γ＝（γL，γR）は定数である。
また、τ2iは自己抑制因子ｖiの変化特性を規定する時定数である。さらに、ｗi/j（＜０）は複数の第２要素（神経要素）ｉおよびｊの相関関係を第２モデルにしたがって生成される複数の（第２要素の）出力振動信号の相関関係に反映させるための第２相関係数（定数）である。「λL」および「λR」は慣れ係数である。Ｋは股関節角度φHに応じたフィードバック係数である。
【００７６】
「ｆ1」および「ｆ2」はそれぞれ次式（３．２）および（３．３）により定義される関数である。
【００７７】
ｆ1（ω）≡ｃ・ω（ｃ＞０）  ‥（３．２）
ｆ2（ω）≡ｃ0＋ｃ1ω＋ｃ2ω2  ‥（３．３）
【００７８】
新たな固有角速度ωMの関数であるｆ1（ωM）およびｆ2（ωM）の係数ｃ，ｃ0，ｃ1，ｃ2は、目標運動設定部２１１によって設定された目標となる運動リズムに応じた係数として設定され得る。
【００７９】
なお、股関節角度φHに代えてまたは加えて、股関節角速度ｄφH／ｄｔや、膝関節、足関節、肩関節、肘関節の角度や角速度、さらには歩行者の着地音、呼吸音、意図的な発声音等、歩行運動リズムと連関したリズムで変動する種々の振動子に基づき、２次振動子ｙiが生成されてもよい。
【００８０】
２次振動子ｙiは、状態変数ｕiの値が閾値ｕth未満である場合は０、状態変数ｕiの値が閾値ｕth以上である場合はこのｕiの値をとる。或いは、２次振動子ｙiは、シグモイド関数ｆsによって定義されている（式（３）参照）。これにより、大腿部の屈曲方向（前方）への動きについてはこの動きを支配する第２要素（神経要素）L+、R+の出力である２次振動子ｙL+およびｙR+が、他の第２要素の出力よりも大きくなる。また、大腿部の伸展方向（後方）への動きについてはこの動きを支配する第２素子L-、R-の出力である２次振動子ｙL-およびｙR-が、他の第２要素の出力よりも大きくなる。脚体（大腿部）の前方または後方への動きは、たとえば、股関節角速度ｄφH／ｄｔの極性によって識別される。
【００８１】
前記のように股関節角度φH等の第２運動振動子に応じて出力が変動する複数の第２要素の関係を表現する第２モデルにしたがって、第２要素の出力として２次振動子ｙiが生成される（式（３），図２／ｓ０４０参照）。これにより、ユーザの実際の運動に関係する複数の第２要素の相関関係を、第２相関係数ｗi/jにより表される第２モデルにおける複数の第２要素の相関関係に反映させることで、当該複数の第２要素間の関係に鑑みて適当な２次振動子ｙiが生成され得る。たとえば、複数の第２要素としてユーザの複数の神経（ニューロン）が想定された場合、左右の脚体による歩行を支配するニューロン間の定性的関係等が反映された形で２次振動子ｙiが生成される。したがって、ユーザの運動を誘導する誘導振動子のリズムおよびスケールを当該関係に鑑みて適当なものとすることができる。
【００８２】
次に、誘導振動子生成部２５０が運動振動子測定部２１０によって測定された股関節角度φH（第３運動振動子）および股関節角速度ｄφH／ｄｔ（第４運動振動子）と、２次振動子生成部２４０により生成された２次振動子ｙiと、固有角速度設定部２３０によって設定された新たな固有角速度ωMとに基づき、誘導振動子ｚを設定する（図３／ｓ０５０）。
【００８３】
具体的には、次式（４）にしたがって第１誘導振動子ｚ1が生成される（図３／ｓ０５１）。
【００８４】
ｚ1L＝ｇ1+（ωmL）ｇ+（φHL）ｙL+−ｇ1-（ωmL）ｇ-（φHL）ｙL-，
ｚ1R＝ｇ1+（ωmR）ｇ+（φHR）ｙR+−ｇ1-（ωmR）ｇ-（φHR）ｙＲR-  ‥（４）
【００８５】
ここで「ｇ1+」「ｇ1-」「ｇ+」および「ｇ-」は次式（４．１）〜（４．４）のそれぞれによって定義される関数である。
【００８６】
ｇ1+（ｐ，ω）≡Σkａk+ωk  （ａk+：係数，ｋ＝０〜３）‥（４．１）
ｇ1-（ｐ，ω）≡Σkａk-ωk  （ａk-：係数，ｋ＝０〜３）‥（４．２）
ｇ+（φ）≡ｃ1+（φ−φ0+）＋ｃ2+（φ−φ0+）3
（ｃ1+，ｃ2;：係数，φ0+：屈曲方向の股関節角度φHの目標値）‥（４．３）
ｇ-（φ）≡ｃ1-（φ−φ0-）＋ｃ2-（φ−φ0-）3
（ｃ1-，ｃ2-：係数，φ0-：伸展方向の股関節角度φHの目標値）‥（４．４）
【００８７】
第１誘導振動子ｚ1は、第１係数ｇ1+およびｇ1-をそれぞれバネ係数（弾性係数）とする、図４に示されている２つの仮想的なバネＧ1+およびＧ1-の弾性力として認識される。第１係数ｇ1+およびｇ1-は、固有角速度ωMに応じて、股関節角度φH（ユーザの運動スケールに応じた運動変数としての第３運動振動子）を目標運動スケールに応じた目標角度φ0+（＞０）およびφ0-（＜０）に近付ける第１ポテンシャル（仮想的なバネ（弾性要素）のポテンシャル）のグラディエントを表す（式（４．１）（４．２）参照）。すなわち、第１誘導振動子ｚ1は、第１ポテンシャルに応じた第１係数ｇ1+，ｇ1-を弾性係数（バネ係数）とし、かつ、股関節角度φH（第３運動振動子）の値を目標角度φ0+およびφ0-に復元させる仮想的なバネ等の弾性要素による弾性力として表現される。これにより、筋肉の収縮状態から伸展状態への移行時の弾性力等、ユーザの身体の弾性要素が反映されたリズムおよびスケールをもってユーザの運動が誘導され得る。
【００８８】
仮想的なバネＧ1+による弾性力は、そのバネ係数ｇ1+に応じて股関節角度φH（第３運動振動子）をその目標角度φ0+に近付けるようにユーザの大腿部に作用する（式（４）参照）。すなわち、股関節角度φHが目標角度φ0+未満である場合、バネＧ1+による弾性力が、股関節角度φHを増加させる方向（前方）に大腿部を動かすようにこの大腿部に作用する。また、股関節角度φHが目標角度φ0+を超えた場合、バネＧ1+による弾性力が、股関節角度φHを減少させる方向（後方）に大腿部を動かすようにこの大腿部に作用する。
【００８９】
また、他の仮想的なバネＧ1-による弾性力は、そのバネ係数ｇ1-に応じて、股関節角度φH（第３運動振動子）をこの目標角度φ0-に近付けるようにユーザの大腿部に作用する（式（４）参照）。すなわち、股関節角度φHが目標角度φ0-を超えている場合、バネＧ1-による弾性力が、股関節角度φHを減少させる方向（後方）に大腿部を動かすようにこの大腿部に作用する。また、股関節角度φHが目標角度φ0-を下回った場合、バネＧ1-による弾性力が、股関節角度φHを増加させる方向（前方）に大腿部を動かすようにこの大腿部に作用する。
【００９０】
また、前記のように大腿部の前方への動きおよび後方への動きの別に応じて、複数の第２要素ｉ（=L+,L-,R+,R-）のうち一部から偏重的に出力があるので、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。
【００９１】
すなわち、左の大腿部が前方に動いているとき、この動きを支配する第２要素L+に応じた２次振動子ｙL+の値が他の第２要素L-に応じた２次振動子ｙL-の値より大きくなって、式（４）により表されている第１誘導振動ｚ1Lが次式（４ａ）のように表される（または近似される）。
【００９２】
ｚ1L＝ｇ1+（ωmL）ｇ+（φHL）ｙL+  ‥（４ａ）
【００９３】
このため、左の大腿部が前方に動いているとき、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力のうち、股関節角度φHを前側の目標角度φ0+に近付けるようにユーザの大腿部に作用するバネＧ1+の弾性力が偏重的に反映されることとなる。これにより、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。
【００９４】
また、左の大腿部が後方に動いているとき、この動きを支配する第２要素L-の出力が他の第２要素L+の出力より大きくなり、これによって第２要素L-に応じた２次振動子ｙL-の値が第２要素L+に応じた２次振動子ｙL+の値より大きくなって式（４）によって表されている第１誘導振動子ｚ1Lは次式（４ｂ）のように表される（または近似される）。
【００９５】
ｚ1L＝−ｇ1-（ωmL）ｇ-（φHL）ｙL-  ‥（４ｂ）
【００９６】
このため、左の大腿部が後方に動いているとき、２つの仮想のバネＧ1+およびＧ１−のそれぞれの弾性力のうち、股関節角度φHを後側の目標角度φ0-に近付けるようにユーザの大腿部に作用する仮想のバネＧ1-の弾性力が偏重的に反映されることとなる。これにより、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避される。これは、右の脚体（大腿部）の動きについても同様である。
【００９７】
なお、股関節角速度ｄφH／ｄｔを変数とするシグモイド関数ｆs（式（３）参照）が第１係数ｇ1+，ｇ1-に組み込まれ、これにより股関節角速度ｄφH／ｄｔの極性により特定される大腿部の前後への動きの別に応じて、複数の第２要素ｉの出力としての２次振動子ｙiのうち一部が偏重的に反映された形で第１トルクＴ1が生成されてもよい。これによっても、２つの仮想的なバネＧ1+およびＧ1-のそれぞれの弾性力が相殺される事態が回避され得る。
【００９８】
さらに、次式（５）にしたがって第２誘導振動子ｚ2が設定される（図２／ｓ０５２）。
【００９９】
ｚ2L＝−ｇ2+（ωmL）（ｄφHL／ｄｔ）Ｈ+（φHL）ｙL+
＋ｇ2-（ωmL）（ｄφHL／ｄｔ）Ｈ-（φHL）ｙL-，
ｚ2R＝−ｇ2+（ωmR）（ｄφHR／ｄｔ）Ｈ+（φHR）ｙR+
＋ｇ2-（ωmR）（ｄφHR／ｄｔ）Ｈ-（φHR）ｙR-‥（５）
【０１００】
ここで「ｇ2+」「ｇ2-」「Ｈ+」および「Ｈ-」は次式（５．１）〜（５．４）のそれぞれによって定義される関数である。
ｇ2+（ω）≡Σkｂk+ωk  （ｂK+：係数，ｋ＝０〜３）‥（５．１）
ｇ2-（ω）≡Σkｂk-ωk  （ｂk-：係数，ｋ＝０〜３）‥（５．２）
Ｈ+（φ）≡０（φ≦０），１（φ＞０）  ‥（５．３）
Ｈ-（φ）≡０（φ＞０），１（φ≦０）  ‥（５．４）
【０１０１】
第２誘導振動子ｚ2は、第２係数ｇ2+およびｇ2-をそれぞれダンパ係数（減衰係数）とする、図４に示されている２つの仮想的なダンパＧ2;およびＧ2-の減衰力として把握される。第２係数ｇ2+およびｇ2-は、固有角速度ωMに応じて股関節角度φHの絶対値の増大を抑制する第２ポテンシャル（仮想的なダンパ（減衰要素）のポテンシャル）のグラディエントを表す（式（５．１）（５．２）参照）。すなわち、第２誘導振動子ｚ2は、第２ポテンシャルに応じた第２係数ｇ2+，ｇ2-を減数係数（ダンパ係数）とし、かつ、股関節角速度ｄφH／ｄｔ（第４運動振動子）に応じて股関節角度φH（第３運動振動子）の絶対値の増大を抑制する仮想的なダンパ等の減衰要素による減衰力として表現される。これにより、筋肉の伸展状態から屈曲状態への移行時の粘性力等、ユーザの身体の減衰要素が反映されたリズムおよびスケールをもってユーザの運動が誘導され得る。
【０１０２】
一方の仮想的なダンパＧ2+による減衰力は、そのダンパ係数ｇ2+および股関節角速度ｄφH／ｄｔ（第４運動振動子）に応じて、前側（屈曲側）への股関節角度φＨの絶対値の増大を抑制するようにユーザの大腿部に作用する（式（５）参照）。すなわち、仮想的なダンパＧ2+による減衰力は、大腿部の前方への過剰な動きを抑制するようにこの大腿部に作用する。
【０１０３】
また、他方の仮想的なダンパＧ2-による弾性力は、そのダンパ係数ｇ2-および股関節角速度ｄφH／ｄｔ（第４運動振動子）に応じて、後側（伸展側）への股関節角度φHの絶対値の増大を抑制するようにユーザの大腿部に作用する（式（５）参照）。すなわち、仮想的なダンパＧ2-による減衰力は、大腿部の後方への過剰な動きを抑制するようにこの大腿部に作用する。
【０１０４】
また、第２誘導振動子ｚ2には、股関節角度φHの関数としての階段関数Ｈ+，Ｈ-が含まれている。したがって、２つの仮想的なダンパＧ2+およびＧ2-のそれぞれの減衰力が相殺される事態が回避される。
【０１０５】
そして、誘導振動子生成部２５０により生成された第１誘導振動子ｚ1＝（ｚ1L，ｚ1R）と、第２誘導振動子ｚ2＝（ｚ2L，ｚ2R） とを含む誘導振動子ｚ（＝ｚ1＋ｚ2）に応じた電流Ｉ＝（ＩL，ＩR）が電池２０６から左右のアクチュエータ２１０にそれぞれ供給され、ユーザの大腿部に力（股関節回りのトルク）Ｔが作用する。
【０１０６】
以後、前記処理（図２／ｓ０１１，ｓ０１２，・・，ｓ０４０，図３／ｓ０５０）が繰り返されることで、ユーザは歩行運動誘導装置２０の動作によって股関節回りのトルクＴが作用している状態で歩行運動する。
【０１０７】
そして、ユーザが前記のようにトレッドミル１０および歩行運動誘導装置２０の動作により歩行運動が誘導されている状態で、運動変数測定部１１０が歩行比ｋを測定する（図３／ｓ１１０）。具体的には、第１運動変数測定部１０１がトレッドミル１０のベルト１３の速度に応じた速度センサ１０２の出力に基づき、ユーザの歩行速度ｖを測定する（図３／ｓ１１１）。また、第２運動変数測定部１１２がトレッドミル１０の踏板１４が受ける圧力に応じた圧力センサ１０４の出力がピークを示す単位時間当たりの回数に基づき、ユーザの歩行率（単位時間当たりの歩数）ｐを測定する（図２／Ｓ１１２）。なお、歩行運動誘導装置２０によりユーザに作用させられる股関節回りのトルクＴの時間変化に基づいて歩行率ｐが測定されてもよい。また、ユーザの身体に加速度センサが取り付けられ、ユーザの鉛直方向の加速度に応じた加速度センサの出力に基づいて歩行率ｐが測定されてもよい。そして、運動変数測定部１００が歩行速度ｖおよび歩行率ｐの関数としての歩行比ｋ＝ｖ／ｐ2を測定する。
【０１０８】
また、生理変数測定部１２０が生理変数センサ１０６の出力に基づいて整理変数ｂを測定する（図３／ｓ１２０）。生理変数ｂには、ユーザの酸素摂取量、呼吸周波数、心拍数、脈拍数、血圧、血中飽和酸素濃度、乳酸値、筋電位、および消費エネルギーのうち一部または全部が含まれる。
【０１０９】
生理変数ｂが測定されるたび、股関節角度φHの目標角度φ0+およびφ0-、第１係数ｇ1+（ωM）およびｇ1-（ωM）のそれぞれに含まれる係数ａk+およびａk-、ならびに第２係数ｇ2+（ωM）およびｇ2-（ωM）のそれぞれに含まれる係数ｂk+およびｂk-のうち一部または全部が変更設定される。目標角度φ0+およびφ0-が変更設定されることで歩行運動誘導装置２０による誘導スケールが変化し、ユーザの歩幅ｑが変更されうる。また、係数ａk+，ａk-，ｂk+およびｂk-の一部または全部が変更設定されることで歩行運動誘導装置２０による誘導リズムが変化し、ユーザの歩行率（目標運動リズム）ｐが変更されうる。これにより、歩行運動誘導装置２０によって歩行比ｋが逐次変化するようにユーザの歩行運動が誘導されていく過程で、生理変数ｂが逐次測定される。そして、歩行比ｋおよび生理変数ｂの関係が求められ、推奨歩行比設定部１３０がこの関係に基づいて推奨歩行比ｋ0を設定する（図３／ｓ１３０）。
【０１１０】
たとえば、生理変数ｂとしてユーザの消費エネルギーが測定され、図５に丸（○）で示されているような歩行比ｋおよび消費エネルギー（生理変数）の関係が得られた場合を考える。この場合、当該関係は図５に示されている歩行比ｋの２次曲線によって近似表現され、この２次曲線が最小値をとる歩行比ｋ＝０．０７５がこのユーザにとっての推奨歩行比ｋ0として設定される。ユーザの消費エネルギーが最小となるように設定された推奨歩行比ｋ0は、ユーザの身体負荷の軽減に鑑みて適当な歩行比ｋである。なお、心拍数、乳酸値（身体の疲労度を表す。）、筋群活性度（身体機能の活性化の程度を表す。）等、消費エネルギーとは異なる生理変数ｂが測定され、当該測定変数ｂに基づいて推奨歩行比ｋ0が設定されてもよい。
【０１１１】
次に、歩行運動誘導装置２０による推奨歩行比ｋ0に応じたユーザの歩行運動の誘導方法について図６を用いて説明する。推奨歩行比ｋ0に応じたユーザの歩行運動の誘導方法は、推奨歩行比ｋ0の設定時におけるユーザの歩行運動の誘導方法と、誘導振動子ｚの生成方法（図３／ｓ０５０）において相違するのみである。そこで、運動振動子の測定（図２／ｓ０１１，ｓ０１２）、１次振動子ｘの生成（図２／ｓ０２０）、固有角速度ωMの設定（図２／ｓ０３０）および２次振動子ｙの生成（図２／ｓ０４０）については説明を省略する。
【０１１２】
第１運動変数測定部１１１が、トレッドミル１０の動作および歩行運動誘導装置２０の動作により歩行運動が誘導されているユーザの歩行速度ｖを測定する（図６／ｓ１１１）。
【０１１３】
また、誘導振動子生成部２５０が、運動振動子測定部２１０により測定された股関節角度φHおよび股関節角速度ｄφH／ｄｔと、２次振動子生成部２４０により生成された２次振動子ｙiと、固有角速度設定部２３０によって設定された固有角速度ωMとに加え、さらに、第１運動変数測定部１１１により測定されたユーザの歩行速度ｖと、推奨歩行比設定部１３０により測定された推奨歩行比ｋ0とに基づき、誘導振動子ｚを設定する（図２／ｓ２５０）。
【０１１４】
具体的には、次式（６）にしたがって第１誘導振動子ｚ1が生成される（図６／ｓ２５１）。
【０１１５】
ｚ1L＝ｇ1+（ｐ0，ωmL）ｇ+（ｑ0，φHL）ｙL+−ｇ1-（ｐ0，ωmL）ｇ-（ｑ0，φHL）ｙL-，
ｚ1R＝ｇ1+（ｐ0，ωmR）ｇ+（ｑ0，φHR）ｙR+−ｇ1-（ｐ0，ωmR）ｇ-（ｑ0，φHR）ｙＲR-，
ｐ0≡（ｖ／ｋ0）1/2，ｑ0≡（ｖｋ0）1/2‥（６）
【０１１６】
ここで「ｇ1+」「ｇ1-」「ｇ+」および「ｇ-」は次式（６．１）〜（６．４）のそれぞれによって定義される関数である。
【０１１７】
ｇ1+（ｐ0，ω）≡Σkａk+（ｐ0）ωk（ａk+（ｐ0）：係数，ｋ＝０〜３）‥（６．１）
ｇ1-（ｐ0，ω）≡Σkａk-（ｐ0）ωk（ａk-（ｐ0）：係数，ｋ＝０〜３）‥（６．２）
ｇ+（ｑ0，φ）≡ｃ1+（φ−φ0+（ｑ0））＋ｃ2+（φ−φ0+（ｑ0））3
（ｃ1+，ｃ2;：係数，φ0+：屈曲方向の股関節角度φHの目標値）‥（６．３）
ｇ-（ｑ0，φ）≡ｃ1-（φ−φ0-（ｑ0））＋ｃ2-（φ−φ0-（ｑ0））3
（ｃ1-，ｃ2-：係数，φ0-：伸展方向の股関節角度φHの目標値）‥（６．４）
【０１１８】
第１誘導振動子ｚ1により、股関節角度φH（第３運動振動子）が推奨歩行比ｋ0および歩行速度ｖに応じた目標角度φ0+，φ0-に一致するようにユーザの歩行運動が誘導されうる。
【０１１９】
さらに、次式（７）にしたがって第２誘導振動子ｚ2が設定される（図６／ｓ２５２）。
【０１２０】
ｚ2L＝−ｇ2+（ｐ0，ωmL）（ｄφHL／ｄｔ）Ｈ+（φHL）ｙL+
＋ｇ2-（ｐ0，ωmL）（ｄφHL／ｄｔ）Ｈ-（φHL）ｙL-，
ｚ2R＝−ｇ2+（ｐ0，ωmR）（ｄφHR／ｄｔ）Ｈ+（φHR）ｙR+
＋ｇ2-（ｐ0，ωmR）（ｄφHR／ｄｔ）Ｈ-（φHR）ｙR-‥（７）
【０１２１】
ここで「ｇ2+」「ｇ2-」は次式（７．１）〜（７．２）のそれぞれによって定義される関数である。
【０１２２】
ｇ2+（ｐ0，ω）≡Σkｂk+（ｐ0）ωk
（ｂK+（ｐ0）：係数，ｋ＝０〜３）‥（５．１）
ｇ2-（ｐ0，ω）≡Σkｂk-（ｐ0）ωk
（ｂk-（ｐ0）：係数，ｋ＝０〜３）‥（５．２）
【０１２３】
第２誘導振動子ｚ2により、股関節角速度ｄφH／ｄｔ（第４運動振動子）に応じて股関節角度φH（第３運動振動子）の絶対値の増大が抑制されるようにユーザの歩行運動が誘導されうる。
【０１２４】
股関節角度φHの目標角度φ0+（ｑ0）およびφ0-（ｑ0）は推奨歩行比ｋ0および歩行速度ｖに応じた推奨歩幅ｑ0（＝（ｖｋ0）1/2）関数であり、ユーザの歩幅ｑ（＝ｖ／ｐ）と推奨歩幅ｑ0との偏差δｑに基づいて補正されうる。
【０１２５】
また、第１係数ｇ1+（ｐ0，ωM），ｇ1-（ｐ0，ωM）に含まれる係数ａk+（ｐ0）およびａk-（ｐ0）、第２係数ｇ2+（ｐ0，ωM），ｇ2-（ｐ0，ωM）に含まれる係数ｂk+（ｐ0）およびｂk-（ｐ0）は推奨歩行比ｋ0および歩行速度ｖに応じた推奨歩行率ｐ0（＝（ｖ／ｋ0）1/2）の関数であり、ユーザの歩行率ｐ（＝ｖ／ｑ）と推奨補効率ｐ0との偏差δｐに基づいて補正されうる。
【０１２６】
そして、誘導振動子生成部２５０により生成された第１誘導振動子ｚ1と、第２誘導振動子ｚ2とを含む誘導振動子ｚに応じた電流Ｉが電池２０６から左右のアクチュエータ２１０にそれぞれ供給され、ユーザの大腿部に力（股関節回りのトルク）Ｔが作用する。
【０１２７】
以後、前記処理（図２／ｓ０１１，ｓ０１２，・・，ｓ０４０，図６／ｓ２５０）が繰り返されることで、ユーザは歩行運動誘導装置２０の動作によって股関節回りのトルクＴが作用している状態で歩行する。
【０１２８】
前記機能を発揮する運動管理システム１によれば、歩行運動誘導装置（第２運動誘導装置）２０の動作によりユーザの運動が誘導されている過程で歩行比（運動変数）ｋおよび生理変数ｂが測定される（図３／ｓ１１０，ｓ１２０）。そして、この動物の「身体機能の活性化」または「身体負荷の軽減」に鑑みて生理変数ｂが適当な値となった状態における歩行比ｋの値が推奨目標比ｋ0として設定される。これにより、推奨歩行比ｋ0が、このユーザの身体機能の活性化等に鑑みて適当に設定されうる。そして、歩行運動誘導装置２０によりこの推奨歩行比ｋ0に基づいてユーザの歩行運動が誘導されることで、ユーザの歩幅（運動スケール）ｑおよび歩行率（運動リズム）ｐがこのユーザの身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【０１２９】
運動管理システム１の効果に関する実験結果について図７〜図９を用いて説明する。
【０１３０】
図７に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の歩行比ｋが異なる。すなわち、ユーザが歩行補助装置２０を装着せずに歩行した場合の歩行比ｋは約０．００６５に制御されている。これに対して、ユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合の歩行比ｋは、前記のように設定された推奨歩行比ｋ（＝０．００７５）に制御されている。このユーザは、歩行運動誘導装置２０の動作によって歩行比ｋ、ひいては歩幅ｑが大きくなるように歩行運動が誘導されている。
【０１３１】
また、図８に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の心拍数（生理変数）が異なる。すなわち、ユーザが歩行補助装置２０を装着せずに歩行した場合の心拍数は、ユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合の心拍数よりも高い。これは、図７に示されているように、ユーザの歩幅が大きくなるように歩行運動が誘導されている状態で、その心拍数が低く抑えられており、身体への負荷が軽減されていることを意味する。
【０１３２】
また、図９に示されているように、歩行運動誘導装置２０による歩行運動の誘導の有無に応じて同一のユーザが同一の速度ｖで歩行した場合の筋群活性度（生理変数）が異なる。図９（ｂ）に示されているユーザが歩行運動誘導装置２０を装着して歩行運動が誘導された状態で歩行した場合のその股関節周囲の筋群活性度は、図９（ａ）に示されているユーザが歩行補助装置２０を装着せずに歩行した場合の同じ箇所の筋群活性度よりも高い。これは、図７に示されているように、ユーザは歩幅が大きくなるように歩行運動が誘導されることで、ユーザの身体機能が活性化されていることを意味する。
【０１３３】
図７〜図９に示されている実験結果から、たとえば、運動管理システム１が高齢者等、歩行運動機能が低下しているユーザの運動管理に利用された場合の意義が大きいことがわかる。すなわち、本発明の運動管理システム１が利用されることで、高齢者等の身体負荷を軽減しながら、その身体機能の低下抑制、さらには活性化を促すことができる。
【０１３４】
さらに、ユーザの歩行速度ｖ（第１運動変数）がトレッドミル１０のベルト１３の速さに基づいて測定される（図３／ｓ１１１）。またユーザの歩行率ｐ（第２運動変数）が、ユーザの着床および離床の繰り返しにより変化するトレッドミル１０への圧力に基づいて測定される（図３／ｓ１１２）。これにより、歩行速度ｖおよび歩行率ｐの関数である歩行比ｋの測定精度の向上が図られる（図３／ｓ１１０）。したがって、ユーザの身体機能の活性化や身体負荷の軽減に鑑みて適当な運動変数の推奨値がより正確に設定されうる。
【０１３５】
また、トレッドミル１０の無端ベルト（回動体）１３を回動させることにより、この無端ベルト１３の動きとは逆向きの歩行または走行運動が誘導される。これにより、当該歩行または走行に伴うユーザの移動が無端ベルト１３の回動によって相殺されるので、トレッドミル１０の設置スペースさえあればユーザの歩行または走行運動に伴う推奨歩行比ｋ0の設定が可能である。トレッドミル１０が省略され、歩行運動誘導装置２０のみによってユーザの歩行運動が誘導されている状態で推奨歩行比ｋ0が設定されてもよい。
【０１３６】
なお、運動管理システム１はユーザの歩行運動以外のあらゆる運動に伴って変化する運動変数の推奨値の設定に利用されてもよい。たとえば、車椅子の左右の車輪に手で力を加える運動について、両腕が車輪を押すリズム等、他の運動についての運動変数の推奨値設定に利用されてもよい。運動管理システム１は馬の走行運動等、人間以外の動物の運動についての運動変数の推奨値設定に利用されてもよい。また、歩行比ｋのほか、歩行速度ｖおよび歩行率ｐの任意の関数としての運動変数の推奨値が設定されてもよい。
【０１３７】
また、トレッドミル１０のほか、複数のローラに掛け渡された無端ベルト、中心または中心から外れた点を通る軸回りに回動される球体、楕円球体、中心軸または中心軸から離れた中心軸に平行な軸回りに回動する円筒、四角筒等の筒体、その他任意の軸回りに回動される一塊の物体等の回動体を回動させることでこの回動体に接触するユーザの腕や脚などの身体部分の運動を誘導する装置が第１運動誘導装置として採用されてもよい。
【０１３８】
前記実施形態では誘導振動子ｚに応じた左右の股関節回りのトルクＴ＝（ＴL，ＴR）がユーザの身体に作用させられたが、他の実施形態として膝関節、足関節、肩関節、肘関節、手根関節等、種々の関節回りのトルクがユーザの身体に作用させられてもよい。トルク作用対象となる関節の組合せは、ユーザに応じてさまざまに変更されてもよい。
【０１３９】
また、第２制御部２００の２次振動子生成部２４０が、運動振動子測定部２１０により測定された股関節角度φH（第２運動振動子）または股関節角速度ｄφH／ｄｔ（第１運動振動子）の周期変化の大きさが閾値を超えた場合、固有角速度設定部２３０により設定された新たな固有角速度ωMに代えて、運動振動子測定部２１０により測定された股関節角度φH等の角速度（位相の時間変化）と、１次振動子生成部２２０により生成された１次振動子ｘの角速度とのうち一方または両方が反映されたリズムで振動する２次振動子ｙを生成してもよい。
【０１４０】
当該構成によれば、ユーザの運動リズムが急に変化した場合でも、この変化後の運動リズムに応じた適当なリズムをもってユーザの運動が誘導されうる。
【０１４１】
また、特開２００４−７３６４９号公報に開示されている方法にしたがって誘導振動子ｚが生成されてもよい。
【０１４２】
さらに、次のような形で誘導振動子ｚが生成されてもよい。すなわち、運動振動子測定部２１０がユーザの２つの異なる身体部分の動きに応じて時間変化するパラメータとしての運動振動子（たとえば肩関節角速度および股関節角度）のそれぞれを「第１運動振動子」および「第２運動振動子」として測定する。また、１次振動子生成部２２０が第１運動振動子と相互に引き込み合い、固有角速度ωMに基づいて定まる角速度で時間変化する１次振動子ｘを生成する。さらに、固有角速度設定部２３０が第１運動振動子と１次振動子ｘとの位相差に基づき、新たな固有角速度ωMを設定する。また、２次振動子生成部２４０が第２運動振動子に基づき、新たな固有角速度ωMに基づいて定まる角速度で時間変化する２次振動子ｙを生成する。そして、誘導振動子生成部２５０が２次振動子ｙに加え、第１運動変数測定部１０１により測定された歩行速度ｖ（第１運動変数）と、第２運動変数測定部１０２により測定された歩行率ｐ（第２運動変数）とのうち一方または両方に基づいて誘導振動子ｚを生成する。
【０１４３】
当該構成によれば、ユーザの異なる身体部分のそれぞれの運動リズムと、この運動を誘導するリズムとの調和を図りながら、ユーザの運動リズムが目標運動リズムに近づくようにこの運動が誘導されうる。
【０１４４】
さらに、本発明の運動管理システム１は、推奨歩行比設定部１３０が複数のユーザのそれぞれの身体的特徴（身体のサイズ、ウェイト、年齢、性別等が含まれる。）を識別するための識別子を認識した上で各ユーザの複数の推奨歩行比ｋ0を予備的に設定し、認識した識別子ごとに予備的に設定済みの推奨歩行比ｋ0に基づいて新たな推奨歩行比ｋ0を設定することにより、識別子および推奨値が対応付けられているデータベースを構築してもよい。この場合、ユーザの識別子、ひいては身体的特徴に応じてデータベースから検索された標準的な推奨値に基づいて、ユーザの歩行運動が誘導される。これにより、ユーザの歩行比ｋがこのユーザの身体機能の活性化等に鑑みて適当なものとなるように、当該誘導スケールおよびリズムが適当に制御されうる。
【図面の簡単な説明】
【０１４５】
【図１】本発明の運動管理システムの構成例示図
【図２】本発明の運動管理方法の説明図
【図３】本発明の運動管理方法の説明図
【図４】歩行運動を誘導する仮想的なバネおよびダンパの説明図
【図５】推奨歩行比の設定方法の説明図
【図６】本発明の運動管理方法説明図
【図７】本発明の運動管理システムの効果に関する実験結果の説明図
【図８】本発明の運動管理システムの効果に関する実験結果の説明図
【図９】本発明の運動管理システムの効果に関する実験結果の説明図
【符号の説明】
【０１４６】
１‥運動管理システム、１０‥トレッドミル（第１運動誘導装置）、１３‥回動体、２０‥歩行運動誘導装置（第２運動誘導装置）、２５‥アクチュエータ、２６‥股関節角度センサ、１００‥第１制御部、１１０‥運動変数測定部、１１１‥第１運動変数測定部、１１２‥第２運動変数測定部、１２０‥生理変数測定部、１３０‥推奨歩行比設定部（推奨値設定部）、２００‥第２制御部（制御部）、２１０‥運動振動子測定部、２２０‥１次振動子生成部、２３０‥固有角速度設定部、２３１‥第１位相差設定部、２３２‥第２位相差設定部、２３３‥相関係数設定部、２３４‥第１角速度設定部、２３５‥第２角速度設定部、２４０‥２次振動子生成部、２５０‥誘導振動子生成部、２５１‥第１誘導振動子生成部、２５２‥第２誘導振動子生成部