(57)【要約】  筋シナジー解析装置（１）は、人体の運動所定部位の拮抗する筋ペア群における各筋に生じる筋電信号を検出する筋電計（１０）と、人体の所定の運動中に筋電計（１０）によって検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出する筋拮抗比算出部（２０２）と、時系列データの各時刻における筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第１の主成分解析部（２０３）とを備えた。これにより、検体の所定の運動毎における筋協調の関係である固有の筋シナジーを得、かつ多変量解析を行って、少ない種類の主成分による運動評価を可能にする。
2011030781.tif 000002
【特許請求の範囲】
【請求項１】
検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出するステップと、前記検体の所定の運動中に検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出するステップと、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出するステップとを備えたことを特徴とする筋シナジー解析方法。
【請求項２】
前記検体の筋拮抗比に基づいて得られた主成分を参照主成分とし、この参照主成分と、新たな検査によって算出された同検体の筋拮抗比に基づいて得られた主成分との差分を算出するステップと、算出された前記差分を第１の表示部に表示するステップとを備えたことを特徴とする請求項１に記載の筋シナジー解析方法。
【請求項３】
前記検体の運動中に検出された筋電信号から前記拮抗する筋ペアの筋電信号のレベルを加算して筋活性度の時系列データを算出するステップと、前記時系列データの各時刻における筋活性度を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出するステップとを備えたことを特徴とする請求項１又は２に記載の筋シナジー解析方法。
【請求項４】
前記検体の筋活性度に基づいて得られた主成分を参照主成分とし、この参照主成分と、新たな検査によって算出された同検体の筋活性度に基づいて得られた主成分との差分を算出するステップと、算出された前記差分を第２の表示部に表示するステップとを備えたことを特徴とする請求項３に記載の筋シナジー解析方法。
【請求項５】
検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出する検出器と、前記検体の所定の運動中に前記検出器によって検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出する筋拮抗比算出手段と、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第１の主成分解析手段とを備えたことを特徴とする筋シナジー解析装置。
【請求項６】
第１の検体と第２の検体における両主成分の差分を算出する第１の差分算出手段と、前記第１の差分算出手段によって算出された差分を第１の表示部に表示する第１の表示制御手段とを備えたことを特徴とする請求項５に記載の筋シナジー解析装置。
【請求項７】
前記第１の検体と前記第２の検体とは同一体であることを特徴とする請求項６に記載の筋シナジー解析装置。
【請求項８】
各主成分の寄与度の累積値が第１の閾値を超えていることを特徴とする請求項５〜７のいずれかに記載の筋シナジー解析装置。
【請求項９】
前記検体の運動中に前記検出器によって検出された筋電信号から前記拮抗する筋ペアの筋電信号のレベルを加算して筋活性度の時系列データを算出する筋活性度算出手段と、前記時系列データの各時刻における筋活性度を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第２の主成分解析手段とを備える請求項５〜８のいずれかに記載の筋シナジー解析装置。
【請求項１０】
第１の検体と第２の検体における両主成分の差分を算出する第２の差分算出手段と、前記第２の差分算出手段によって算出された差分を第２の表示部に表示する第２の表示制御手段とを備えたことを特徴とする請求項９に記載の筋シナジー解析装置。
【請求項１１】
各主成分の寄与度の累積値が第２の閾値を超えていることを特徴とする請求項９又は１０に記載の筋シナジー解析装置。
【請求項１２】
複数の筋に所与の運動を行わせる筋シナジーインターフェースにおいて、
前記複数の拮抗する筋対群に対する運動を規定する運動情報を入力する運動コマンド設定手段と、前記運動情報を主成分情報に変換する第１の変換手段と、前記主成分情報を前記各拮抗する筋対へ印加する電位信号に変換する第２の変換手段とを備えたことを特徴とする筋シナジーインターフェース。
【請求項１３】
前記第１の変換手段は、運動情報としての直交座標系の位置情報を極座標系に変換するものであることを特徴とする請求項１２記載の筋シナジーインターフェース。
【請求項１４】
前記第１の変換手段は、運動情報としての運動周期及び１周期での単位運動を極座標系に変換するものであることを特徴とする請求項１２記載の筋シナジーインターフェース。
【請求項１５】
前記第１の変換手段は、前記運動情報を２種類の主成分情報に変換し、前記第２の変換手段は、得られた主成分情報を前記各拮抗する筋対への電位信号に変換することを特徴とする請求項１２〜１４のいずれかに記載の筋シナジーインターフェース。
【請求項１６】
前記第１の変換手段は、少なくとも前記複数の拮抗する筋群から構成される運動基部と運動端部との間の距離と回転角度とを２種類の主成分情報に置換し、前記第２の変換手段は、得られた２種類の主成分情報を前記各拮抗する筋対への電位信号に変換することを特徴とする請求項１５に記載の筋シナジーインターフェース。
【請求項１７】
前記第１の変換手段は、前記運動実行時の制約に応じた固有値を算出することを特徴とする請求項１２〜１６のいずれかに記載の筋シナジーインターフェース。
【請求項１８】
前記第２の変換手段は、前記主成分情報を筋拮抗比及び活性度に変換するものであることを特徴とする請求項１２〜１７のいずれかに記載の筋シナジーインターフェース。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、筋電信号に基づいて運動する筋対群の筋シナジー関連技術に関する。
【背景技術】
【０００２】
例えばスポーツ運動において効率的にトレーニングを行うためには、対象となる身体運動を適切に評価する指導者が必要である。しかし、指導者はスポーツ人口に比して少なく、かつその評価は経験等に依存しており、必ずしも一定ではない。一方、対象運動を工学的に定量評価しようとする試みがあるが、かかる試みの多くは、運動全体のパフォーマンス又は身体表面に取り付けたセンサからの情報に基づく表面的評価に留まっている。【０００３】
近年、汎用的な運動中の人体内部の筋活動をリアルタイムに推定、可視化するシステムが記載されている（非特許文献１）。このシステムは、検体に取り付けたマーカを検知するモーションキャプチャーシステムと、検体の動きを検知するプレート状の感圧素子から構成される床反力計と、検知結果に対して所与の演算を実行する計算機とから構成され、身体運動のより詳細な解析、深部感覚推定への応用の可能性について多くの注目を集めている。このシステムでは、主に運動学データに基づき、汎用運動時の筋張力を可視化することを可能にしている。【０００４】
また、特許文献１には、筋電計で得られた複数の筋肉の筋電信号より各筋肉の張力を検出し、各筋肉の張力に応じて定まる複数の関節の状態と複数の筋肉の張力との履歴に基づいて、複数の関節の状態の組み合わせの時系列を推定する機械学習のモデルを生成する技術が記載されている。この技術によれば、筋肉の動きが制約されることなく、複数の筋肉の筋電信号に基づいて所望の複数の関節の状態やこれらの状態の組み合わせの時系列を精度よく高速に推定できる。【０００５】
また、非特許文献２には、運動のための筋シナジーの組合せに関して記載されている。より詳細には、非特許文献２には、１３箇所の筋からそれぞれ検出した時系列的な筋電信号を多変量解析によって５種類の主要成分で表される合成変数（得点）に変換する手法が開示されている。このようにして、制御対象モデルの変数の種類を減じることで、モデルに対する運動理解及び機械的な制御が容易となる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００７−１３６０４１号公報
【非特許文献】
【０００７】
【非特許文献１】村井昭彦、黒崎浩介、山根克、中村仁彦、“モーションキャプチャ、ＥＭＧ、筋の動特性モデルに基づく筋張力のリアルタイム推定及び可視化”第１４回ロボテシクスシンポジア、北海道、２００９年３月１６日
【非特許文献２】Ｅ．Ｂｉｚｚｉ，Ｖ．Ｃ．Ｋ，Ｃｈｅｕｎｇ，Ａ．ｄ’Ａｖｅｌｌａ，Ｐ．Ｓａｌｔｉｅｌ、Ｍ．Ｔｒｅｓｃｈ、“Ｃｏｍｂｉｎｉｎｇ  ｍｏｄｕｌｅｓ  ｆｏｒ  ｍｏｖｅｍｅｎｔ”,ｗｗｗ．ｅｌｓｅｖｉｅｒ．ｃｏｍ／ｌｏｃａｔｅ／ｂｒａｉｎｒｅｓｒｅｖ，受理２００７年８月４日、オンライン入手２００７年９月５日
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ところで、非特許文献１は、モーションキャプチャーシステムや床反力計を要するため、規模が大掛かりで、かつ高価であるという問題がある。また、汎用運動を解析対象としていることから膨大な自由度を扱うことになり、リアルタイム計算を実現する上で、計算コストが無視できない。また、運動者毎に身体モデルを同定する必要があるため、適切な身体の筋骨格モデルが必要となる。特許文献１は、装置への学習段階では、三次元位置計測装置が必要となり、装置が大掛かりである。これら非特許文献１及び特許文献１には、複数の変数に対して多変量解析（主成分解析）を施すことは記載されていない。【０００９】
一方、非特許文献２には、多数の筋電信号に対して多変量解析を施して変数を減少させる（すなわち運動パターンを圧縮する）点が記載されている。しかしながら、非特許文献２には、経時的な筋電信号に関して解析が施されているのみで、特定の運動に起因する、すなわちタスク志向型における筋協調に関する発想はない。従って、計算処理に要するコスト、時間は無視できないものとなる。【００１０】
また、本発明者らは、人体等の検体の上肢回転や歩行等の所定の運動における筋協調による固有の筋シナジーを得、かつ多変量解析を行って、少ない種類の主成分による運動評価を可能にする筋シナジー解析方法及びその装置を提案した（特願２００９−２１２１４９）。より具体的には、この筋シナジー解析方法は、検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出するステップと、前記検体の所定の運動中に検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出するステップと、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い例えば第１，第２主成分得点を算出するステップとを備えたものである。これにより、より少ない種類の主成分によって運動を評価可能にした。【００１１】
ところで、上記提案の通り、拮抗する筋ペアに生じる筋電信号と主成分得点との間に高い相関性が存在することが解明できたものの、大脳運動野からの運動コマンドと拮抗する筋ペアへの筋電位コマンドとの間のメカニズムは依然として不明のままであり、制御系としては不充分なものであった。【００１２】
本発明の目的は、検体の所定の運動毎における筋協調の関係である固有の筋シナジーを得、かつ多変量解析を行って、少ない種類の主成分による運動評価を可能にする筋シナジー解析方法及びその装置を提供することにある。【００１３】
また、本発明の他の目的は、多数の筋のそれぞれに対する活動コマンドを数少ない運動コマンドから作成することが可能な筋シナジーに基づく筋シナジーインターフェースを提供することにある。
【課題を解決するための手段】
【００１４】
本発明に係る筋シナジー解析方法は、検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出するステップと、前記検体の所定の運動中に検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出するステップと、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出するステップとを備えたことを特徴とするものである。【００１５】
また、本発明に係る筋シナジー解析装置は、検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出する検出器と、前記検体の所定の運動中に前記検出器によって検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出する筋拮抗比算出手段と、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第１の主成分解析手段とを備えたことを特徴とするものである。なお、拮抗する筋ペアの筋電信号レベルの正規化した比、すなわち正規化した筋拮抗比を用いて時系列データ化してもよい。【００１６】
これらの発明によれば、人間等の検体である、検査対象となる身体部位の拮抗する筋ペア群における各筋に生じる筋電信号が、所定の運動中に、検出器で経時方向に検出される。筋拮抗比算出手段によって、検出器によって検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比が時系列データとして算出される。次いで、第１の主成分解析手段によって、時系列データの各時刻における前記筋拮抗比を変数として多変量解析が行われて、相関の寄与度の高い少なくとも第１主成分を含む主成分が算出される。従って、従来に比して、装置がコンパクトになる。【００１７】
このように、各筋の筋電位行列に対して、単に多変量解析（ここでは主成分解析）を施すのではなく、運動における筋の協調を考慮して、筋拮抗比という概念を導入し、この概念に対して、主成分解析を施すことで、変数の種類を低減させて処理の容易、高速化を実現して、リアルタイム処理が可能となる。また、煩雑かつ困難な身体各部位のパラメータ同定が不要であるため、その分、処理が簡易となる。また、主成分解析後の変数の種類が、第３主成分までであると、表示に際して３次元座標系で表現できるという利点がある。しかも、筋拮抗比を採用することで、運動の種類毎に主成分解析の結果も異なり、運動に即した主成分の数、得点が得られる。【００１８】
また、この算出された主成分得点は、必要に応じて処理される。例えば、単に出力（典型的には、モニタに表示）される場合、他人との比較のために出力される態様、同一人の理想の運動（参照情報）との比較のために出力される態様、同一人の理想の運動との対比のために、自己の理想運動の情報と今回の測定した結果の差分を算出して、算出した差分が表示される態様において適用可能である。特に、自己との対比の態様に関しては、固有の運動に対し理想との対比という点で、タスク志向型を実現可能である。また、熟練運動における運動の巧みさを解明することも可能であり、スポーツ科学、運動生理学の分野に有用となり得る。【００１９】
また、本発明に係る筋シナジーインターフェースは、複数の筋に所与の運動を行わせるインターフェースにおいて、前記複数の拮抗する筋対群に対する運動を規定する運動情報を入力する運動コマンド設定手段と、前記運動情報を主成分情報に変換する第１の変換手段と、前記主成分情報を前記各拮抗する筋対へ印加する電位信号に変換する第２の変換手段とを備えたものである。【００２０】
この発明によれば、運動コマンド設定手段によって、前記複数の拮抗する筋対群に対する運動を規定する運動情報が入力されると、第１の変換手段によって、この運動情報が主成分情報に変換され、次いで、第２の変換手段によって、前記主成分情報が前記各拮抗する筋対へ印加するための電位信号に変換される。そして、各拮抗する筋対は、変換して得られた、対応する各電位信号によって筋活動させられ、これによって、規定された運動が実行される。従って、各拮抗する筋対への運動指示が、高い相関を有するより数少ない主成分情報で行うことができるため、従来に比して処理が高速となり、かつ構成がコンパクトになる。
【発明の効果】
【００２１】
本発明によれば、検体の所定の運動毎における筋協調の関係である固有の筋シナジーを得、かつ多変量解析を行って、少ない種類の主成分による運動評価ができる。【００２２】
また、本発明によれば、各拮抗する筋対への運動指示を、高い整合性を有するより数少ない主成分情報で行わせることができ、従来に比して処理を高速化でき、かつ構成をコンパクト化できる。
【図面の簡単な説明】
【００２３】
【図１】本発明に係る筋シナジー解析装置の一実施形態を示す全体構成図である。
【図２】検体の筋の一例として、歩行運動に関与する右片脚下肢８筋についての筋名と付着位置とを示す模式図である。
【図３】信号処理部によって実行される測定処理の一例を示すフローチャートである。
【図４】信号処理部によって実行される評価処理の一例を示すフローチャートである。
【図５】信号処理部によって実行される評価処理の一例を示すフローチャートである。
【図６】行列Ｒに対して、累積寄与率が９０％以上となる３つの主成分の主成分得点を第１主成分、第２主成分、第３主成分をｘ軸、ｙ軸、ｚ軸にとった図表である。
【図７】図６における各主成分の要素を筋拮抗比順に示した図表である。
【図８】筋活性度を示す図表である。
【図９】筋拮抗比の主成分得点のうち、ここでは第１、２主成分における差分を表示した図である。
【図１０】ＥＭＧ、手先位置および関節角の測定結果を示す図表である。
【図１１】各関節角の遷移を示す図表である。
【図１２】主成分得点の推移を示す図表である。
【図１３】第１主成分得点を横軸に、第２主成分得点を縦軸にプロットした様子を示す図表である。
【図１４】本発明に係る筋シナジーインターフェースの一例を示す構成図である。
【図１５】本発明に係る筋シナジーインターフェースによって実行される処理手順の概要を示すフロー図である。
【図１６】上肢の水平面内運動における運動コマンド生成の流れを示すフロー図である。
【図１７】被験者の頭部、肩から手先までの各筋の位置関係と、各拮抗する筋対の関係を示す筋の模式図である。
【図１８】手先を水平面内で時計回りの円運動をさせたとき、横軸に第１主成分ｗ１を、縦軸に第２主成分ｗ２を採って、（ｗ１，ｗ２）の点の運動位置（時間経過；図中の矢印参照）を表した図表である。
【図１９】手先を水平面内で反時計回りの円運動をさせたとき、横軸に第１主成分ｗ１を、縦軸に第２主成分ｗ２を採って、（ｗ１，ｗ２）の点の運動位置（時間経過；図中の矢印参照）を表した図表である。
【図２０】第１主成分得点ｗ１が肩と手先とを結ぶ直線上に沿った伸縮（ｒ）に対応することを示す重回帰分析結果の図表である。
【図２１】第２主成分得点ｗ２が肩あるいは頭部を中心とした手先の回転（θ）に対応することを示す重回帰分析結果の図表である。
【図２２】クランクタスク（時計回り）における実験値（ｒ、θ）から時計回りの（ｗ１，ｗ２）を推定した結果（上記重回帰モデルを利用）を示す図表である。
【図２３】クランクタスク（時計回り）における実験値（ｒ、θ）から反時計回りの（ｗ１，ｗ２）を推定した結果（上記重回帰モデルを利用）を示す図表である。
【図２４】筋拮抗比の変換の一例を示すフロー図である。
【図２５】人工筋アームの駆動を説明する図表で、（ａ）は実験１，２を示し、（ｂ）は実験３を示している。
【図２６】級数展開における近似結果を示す図表で、（ａ）は歩行（３km／h）を分解した図、（ｂ）（ｃ）はその時の主成分ｗ１，ｗ２を示している。
【図２７】ダフィング方程式による近似結果を示す図表で、（ａ）は歩行（３km／h）を分解した図、（ｂ）はその時の主成分ｗ１に関する情報を示し、（ｃ）はその時の主成分ｗ２に関する情報を示している。
【発明を実施するための形態】
【００２４】
図１は、本発明に係る筋シナジー解析装置の一実施形態を示す全体構成図である。本装置１は、検体の筋電位を検出する筋電計１０と、検出した筋電位に基づいて所定の信号処理を施すマイクロコンピュータ等で構成される信号処理部２０とを備える。本装置１には、必要に応じて、端子を介して有線又は無線でロボット３０が接続されている。ロボット３０は、機構体の所要部位を運動させるための複数のアクチュエータ３１を有する。アクチュエータ３１は、電気信号に応じて空圧を調整し、空気圧ゴムチューブを伸張、収縮可能にするMcKibben型人工筋肉等で構成されており、信号処理部２０からの信号を駆動信号として受信し、動作を行う。アクチュエータ３１はMcKibben型人工筋肉に限定されず、電気信号で、あるいは電気信号を変換した別の力によって機械的運動に変える電気機械変換素子、例えば電磁ソレノイド、圧電素子、モータ等であってもよい。
【００２５】
筋電計１０は、所要数の電極１１と、各電極１１に生じる電気信号を検出する筋電検出回路１２と、検出された電気信号に対して所定の前処理を施す前処理回路１３とを備える。電極１１は、導電性材料からなり、所定形状例えば円形を有し、直径は数ｍｍ〜十数ｍｍが好ましく、本実施形態では１０ｍｍ程度である。筋電検出回路１２は、所要数の接続端子を有し、それぞれの電極１１とリード線もしくは、無線を介して接続されている。筋電検出回路１２は、各端子に入力される筋電位を時分割で順番に繰り返し取り込むものである。あるいは、同期して取り込むものでもよい。筋電位信号は、レベルが数μＶ〜数ｍＶ、周波数が数Hz〜数10KHz程度の交流信号である。そこで、前処理回路１３は、筋電位を処理可能なレベル（数Ｖ）まで増幅するアンプ、筋電位の主要周波数帯である１０Ｈｚ〜１５０Ｈｚの信号のみを通過させるバンドパスフィルタ、及び全波整流回路を備えている。また、前処理回路１３は、筋電位信号をデジタルで処理するべく、出力側にＡＤ変換部を備えている。なお、前処理回路１３の一部、例えばフィルタは信号処理部２０内でソフトウエアとしてのデジタルフィルタで処理する態様としてもよい。
【００２６】
信号処理部２０は、ＣＰＵを有するマイクロコンピュータで構成されており、本発明に係る信号処理のための処理プログラムや所要の画像情報等を記憶するＲＯＭ２１、処理途中のデータを一時的に記憶するＲＡＭ２２、所要の指示を行うためのテンキーやマウス等で構成される操作部２３及びモニタ２４が接続されている。モニタ２４は、操作部２３からの入力情報の確認や処理結果を表示するものである。ＲＯＭ２１内の処理プログラムは、ＲＡＭ２２に読み出され、保持された状態で実行される。
【００２７】
信号処理部２０のＣＰＵは、ＲＡＭ２２上に保持された処理プログラムを実行することによって、筋電位の測定処理を実行する測定処理部２０１、取り込まれた筋電位から所定の拮抗する筋ペアに対する筋拮抗比を算出する筋拮抗比算出部２０２、算出された筋拮抗比に対して多変量解析（主成分解析）を実行して主成分の設定とそれらの得点を算出する第１の主成分解析部２０３、いずれも第１の主成分解析部２０３で算出された値であって、事前に取り込まれた基準となる参照主成分得点と新たな測定によって得られた主成分得点との差分を算出する第１の差分算出部２０４、及び主成分解析の結果をモニタ２４に表示する他、参照主成分得点との差分をモニタ２４に表示する表示制御部２０５として機能する。
【００２８】
また、信号処理部２０のＣＰＵは、ＲＡＭ２２上に保持された処理プログラムを実行することによって、取り込まれた筋電位から所定の拮抗する筋ペアに対する筋活性度を算出する筋活性度算出部２０６、算出された筋活性度に対して多変量解析（主成分解析）を実行して主成分の設定とそれらの得点を算出する第２の主成分解析部２０７、及びいずれも第２の主成分解析部２０７で算出された値であって、事前に取り込まれた基準となる参照主成分得点と新たな測定によって得られた主成分得点との差分を算出する第２の差分算出部２０８として機能する。
【００２９】
なお、表示制御部２０５は、第１の表示制御手段として、第１の差分算出部２０４によって算出された差分のモニタ２４、及び第２の表示制御手段として、第２の差分算出部２０８によって算出された差分のモニタ２４への表示を行う。
【００３０】
第１、第２の主成分解析部２０３，２０７では、さらに、それぞれの主成分の寄与度が算出され、その累積寄与率（累積値）が所定の閾値、例えば９０％を超えるまでの寄与度の高い順の、第１主成分、第２主成分、…が採用される。
【００３１】
測定に先立って、検体のうち所定の運動（固有のタスク、すなわちタスク志向型）を行わせる部位、すなわち測定対象となる筋の表面側の皮膚面に、例えば皮膚抵抗を低減する剤を塗り、その上に電極１１が貼着される。この状態で、検体が前記所定運動、例えば歩行運動を実行する。
【００３２】
図２は、検体の筋の一例として、歩行運動に関与する右片脚下肢８筋についての筋名と付着位置とを示す模式図である。図２には、股関節１関節筋、股関節膝関節２関節筋、膝関節１関節筋、足関節筋の各屈筋・伸筋の計８個の筋ｍ１〜ｍ８が示されている。また各屈筋・伸筋をペアとした８筋を表１に示す。
【００３３】
【表１】
2011030781.tif 000003
【００３４】
表１のｒ１〜ｒ８は、各筋拮抗比の番号を示している。ｒ１〜ｒ４は、一関節筋に関する拮抗する筋ペア群の各比であり、ｒ５〜ｒ８は、二関節筋を含む拮抗する筋ペア群の各比である。表１は、二関節構造における二関節筋と一関節筋との協調による動作に関する報告（熊本水頼編著：“ヒューマノイド工学、生物進化から学ぶ、２関節ロボット機構”、東京、2006）に基づいている。なお、例えば、腕及び手を用いての運動においては、表１とは異なる表が作成される（表２、表３参照。詳細は後述する）。
【００３５】
測定処理部２０１は、検体の運動に合わせて各電極１１からの筋電位を経時方向に検出するものである。例えば、一連の運動を複数に分割し、各区分内での筋電位を平均化する等して、各区分の筋電位を決定する。また筋電位は、予め測定された各筋の最大検出電位に対して百分率で表して（% MaximamVoluntaryContraction：％ＭＶＣ）、出力するものとしている。なお、所定運動を複数回連続して行わせ、各所定運動の検出データで平均するようにして検出精度をアップさせるようにしてもよい。
【００３６】
ここで、各区分の時間をｔ、各区分における各筋ｐの％ＭＶＣをｍｐｔと表すとき、各筋の時間毎の％ＭＶＣを行にとったものを筋順で並べた行列を筋活動データ行列Ｍとすると、行列Ｍは、式１で表される。
【００３７】
【数１】
2011030781.tif 000004
【００３８】
各区分の時間をｔとし、各区分における各筋拮抗比番号ｐに対して算出された筋拮抗比をｒｔｐと表し、各時間の８筋の筋拮抗比を行にとったものを、時系列順に縦に並べた行列をＲとすると、行列Ｒは式２で表される。筋拮抗比算出部２０２は、各行列の数値を算出する。
【００３９】
【数２】
2011030781.tif 000005
【００４０】
次いで、第１の主成分解析部２０３によって、行列Ｒに対して主成分解析が施されると、各筋拮抗比の平均値Ｅｒが、式３のように得られる。
【００４１】
【数３】
2011030781.tif 000006
【００４２】
行列Ｒは、主成分解析によって筋協調を求めるためのものである。また、筋活性度算出部２０６は、式１において、拮抗する筋ペア同士の％ＭＶＣの値を加算する。加算値はそれぞれ筋活性度という。上記と同様な考え方により、第２の主成分解析部２０７によって、筋活性度に対しても主成分解析が行われる。なお、筋活性度算出部２０６は、筋拮抗比算出部２０２における拮抗する筋ペアと同一のペアを対象とすることが好ましいが、限定される必要はない。また、第１、第２の主成分解析部２０３，２０７によって、筋数のパラメータが主成分の種類数に減少、すなわち圧縮されることで、その後の処理負担が減り、迅速となる。
【００４３】
主成分解析により算出された主成分のうち、寄与度の高いｎ個の主成分を決定する。これらの主成分の中で、第ｋ番目に高い寄与率を持つものを第ｋ主成分といい、第ｋ主成分に対する固有値λｋは、式４のように求められ、対応する寄与率ｃｋ、累積寄与率Ｃｋを、式５，式６のように表す。
【００４４】
【数４】
2011030781.tif 000007
【００４５】
【数５】
2011030781.tif 000008
【００４６】
【数６】
2011030781.tif 000009
【００４７】
また、第１、第２の差分算出部２０４，２０８は、事前に取り込まれた基準となる参照主成分得点と新たな測定によって得られた主成分得点との差分を算出するものである。
【００４８】
次に、図３は、信号処理部２０によって実行される測定処理の一例を示すフローチャートである。このフローチャーは、検体の基準となる参照情報を取得するものである。まず、処理に先立って（あるいは処理の最後であってもよい）、検体を識別する情報、例えば氏名の登録が行われる。
【００４９】
この後、各電極１１から前処理回路１３を介して周期的に筋電位が取得され（ステップＳ１）、取得された筋電位データは参照データとして例えばＲＡＭ２２に蓄積される（ステップＳ３）。次いで、筋拮抗比が算出され（ステップＳ５）、算出された筋拮抗比に対して第１の主成分解析が実行される（ステップＳ７）。
【００５０】
続いて、参照データから筋活性度が算出され（ステップＳ９）、算出された筋活性度に対して第２の主成分解析が実行される（ステップＳ１１）。次いで、各主成分解析の結果がモニタ２４に確認的に表示される（ステップＳ１３）。モニタ２４の画像を見て、異常との操作部２３からの入力がなければ、参照データを取得するための測定が終了する。なお、測定結果に異常があれば、同様にして測定をやり直せばよい。また、本実施形態では、筋電位データを格納する一方、第１、第２の主成分解析は確認表示のために行っているが、格納していない。一方、第１、第２の主成分解析の結果を、予め格納する態様としてもよい。予め格納する態様では、後の差分表示において、その都度計算する必要がない。
【００５１】
図４、図５は、信号処理部２０によって実行される評価処理の一例を示すフローチャートである。このフローチャートは、一般的には、図３の検体と同一の検体について行われる。すなわち、評価に先立って、評価対象となる検体の氏名が入力され、入力された氏名と既に登録されている氏名との照合が行われ、一致する氏名があれば、ステップＳ２１に移行する。一方、不一致であれば、未登録検体として、好ましくは図３のフローチャートが実行される。
【００５２】
次いで、各電極１１から前処理回路１３を介して周期的に筋電位が取得され（ステップＳ２１）、取得された筋電位データがＲＡＭ２２に保存される（ステップＳ２３）。次いで、筋拮抗比の評価指示の有無が判断される（ステップＳ２５）。この判断が肯定されると、保存データから筋拮抗比の算出が実行され（ステップＳ２７）、続いて、第１の主成分解析が実行されて、主成分得点が得られる（ステップＳ２９）。
【００５３】
続いて、同一の検体について格納されている参照筋電位データが読み出され、参照用の筋拮抗比が算出され（ステップＳ３１）、参照用の第１の主成分解析が実行される（ステップＳ３３）。
【００５４】
次いで、ステップＳ２９で得られた主成分得点と、ステップＳ３３で得られた主成分得点との差分が算出され（ステップＳ３５）、その結果がモニタ２４に表示される（ステップＳ３７）。主成分は第１主成分、第２主成分、…というように、複数種類であることから、それぞれを直交座標軸としてモニタ２４に表示する。両主成分得点の差分の計算は、各座標軸単位で行ってモニタ２４に表示する単元式でもよいし、複数元の差分をまとめて演算し、その結果を表示する方法でもよい。
【００５５】
差分のモニタ２４への表示中は、表示の終了、例えば終了キーが押下されたか否かが判断される（ステップＳ５１）。終了でなければ、ステップＳ５３に進む。
【００５６】
一方、ステップＳ２５で、判断が否定されると、保存データから筋活性度の算出が実行され（ステップＳ３９）、続いて、第２の主成分解析が実行されて、主成分得点が得られる（ステップＳ４１）。続いて、同一の検体について格納されている参照筋電位データが読み出され、参照用の筋活性度が算出され（ステップＳ４３）、参照用の第２の主成分解析が実行される（ステップＳ４５）。
【００５７】
次いで、ステップＳ４１で得られた主成分得点と、ステップＳ４５で得られた主成分得点との差分が算出され（ステップＳ４７）、その結果がモニタ２４に表示される（ステップＳ４９）。この表示は、ステップＳ３７と略同様の方法で行われる。
【００５８】
そして、ステップＳ５１で、終了か否かの判断が行われ、終了であれば、すなわち、終了キーなどが押下されたのであれば、本フローを抜ける。一方、終了でなければ、ステップＳ５３で、他方側の差分表示に移行する指示の有無が判断される。例えば、この指示は、例えば特定のキーを割り当てておけばよい。他方側の差分表示に移行する指示がなければ、ステップＳ５１に移行する。一方、他方側の差分表示に移行する指示があれば、現在が第１の主成分解析による得点の差分表示状態か否かが判断され（ステップＳ５５）、この判断が肯定されれば、第２の主成分解析による得点の算出、差分表示状態に切り換えられ（ステップＳ３９へ移行）、判断が否定されれば、第１の主成分解析による得点の算出、差分表示状態に切り換えられる（ステップＳ２７９へ移行）。なお、前記主成分解析の算出、及び差分表示が行われ、そのデータが一時的に保存されることで、かかる主成分解析の算出及び差分表示が行われていることが判断された場合に、ステップＳ５５の後、第２の主成分解析による得点の差分表示状態に切り換えられ（ステップＳ４９へ移行）、判断が否定されれば、第１の主成分解析による得点の差分表示状態に切り換えられる（ステップＳ３７へ移行）。
【００５９】
上記構成及び処理によれば、検電計１０とコンピュータからなる信号処理部２０とで装置が構成されるので、小型で安価な装置を提供することが可能となり、一般的な普及が容易となる。
【００６０】
また、筋拮抗比の主成分解析による主成分得点と参照主成分得点との差分、また筋活性度の主成分解析による主成分得点と参照主成分得点との差分を表示することで、定量評価が可能となり、理想運動に近づけるために必要な筋シナジーの修正法（身体の動かし方）を教示することができる。さらに、差分結果を評価情報とすることで、ユーザへの直接的運動教示が可能なロボットシステムや機能的電気刺激システムに応用することで、新たなコーチングやリハビリテーションへの適用が可能となる。
【００６１】
次に、実験例について説明する。
【００６２】
＜実験例１＞
検体として健常成人男性の歩行運動に対して、ＥＭＧ（筋活動を直接計測する筋電図：electromyogram）計測を行い、筋協調を重視し、筋協調の最小単位と考える拮抗する筋同士の比（筋拮抗比）の概念を導入し、それを用いて主成分解析を行う。これにより８つの筋拮抗比が、３つのパターンに圧縮できた。また、拮抗する筋同士の活動を足し合わせる、筋活性度の概念も導入し、今回用いた筋拮抗比という概念と合わせて、筋協調を考える上で拮抗する筋群の働きが非常に重要な役割を果たしていることを示す。
【００６３】
（実験方法）
本研究ではトレッドミル歩行中の筋活動に焦点を当て、健常成人男性（23歳、身長174.5[cm]、体重48[kg]、利き足右）に自然歩行時の右片脚下肢８筋を選択し、筋電位を計測した（図２参照）。
【００６４】
（実験機器）
筋電位の計測は、１０[mm]の表面電極を電極間距離２０[mm]にて使用し、筋電位取得装置および生体アンプWEB-9000（日本光電社製）で筋電位を増幅したのち、ＡＤ変換記録装置LX-110（TEAC社製）と画像記録装置AQVU（TEAC社製）とを用いて、カメラ画像を同期して記録した。カメラ画像データは歩行相を判断するために用いた。歩行相は右踵接地から次の右踵接地までを１歩行周期とし、両脚が接地している間を両脚支持期、右下肢のみが接地している間を単脚支持期、右下肢が接地していない間を遊脚期と呼び、１歩行周期は順に、両脚支持期・単脚支持期・両脚支持期、遊脚期と分別される。また、サンプリング周波数は１[KHz]とした。
【００６５】
（被験筋）
検体の被験筋は、二関節構造における二関節筋と一関節筋の協調による動作に関する前記報告（熊本水頼編著；前述）を元に、股関節１関節筋、股関節膝関節２関節筋、膝関節１関節筋、足関節筋の各屈筋・伸筋をペアとした８筋を対象とした。図２に筋の筋名と付着位置を示す。
【００６６】
（実験環境）
検体（被験者）に電極１１を貼り付ける際には、スキンピュア（日本光電社製）で皮膚処理を行い、皮膚抵抗を１０[kΩ]以下とした。歩行は、トレッドミルT650m（SportsArt Fitness社製）を用い、走行速度、3.0[km/h]、4.0[km/h]の設定で１分間ずつ行ってもらい、各速度での歩行運動時の筋電位を計測した。
【００６７】
（データ処理）
筋電位は、前述の通り、筋間での正規化をとるために最大等尺性収縮時の筋電位に対する百分率％ＭＶＣで表す。そのために予め各筋の最大等尺性収縮時の筋電位を計測しておく。実測値、正規化用データにそれぞれ整流化、フィルタリング(10-150[Hz] 帯域通過)、平滑化を施し、実測値を正規化用データで除すことで、データを百分率で表した後、１歩行周期毎でデータを平均化し、時間軸も百分率（1000ポイント）で表した。さらに、5[%]時間ごとにデータを平均し、時間を（t1 . . . t20）で表した。また、各筋のうち作用が相反するもの同士を拮抗する筋ペア群として、お互いの比率を算出し、筋拮抗比とした。今回はその中から、前記表１の組み合わせの筋拮抗比を用いた。
【００６８】
（解析）
得られたデータは８筋の情報が含まれており、これらの情報をより少ないパラメータで表現するために、これらデータに対し主成分分析を行った。これにより、８筋の自由度がさらに少ない要素に圧縮でき、この圧縮した情報が筋協調につながるのではないかと考える。筋拮抗比番号をｐ、各時間をt、各筋拮抗比をｒｔｐとして前記式２〜式６を得ている。なお、筋活動データ行列Ｒは、20×8の行列とした。
【００６９】
算出された主成分ベクトルは、拮抗する筋ペア群の組合せに対応すると考える。
【００７０】
行列Ｒに対して、累積寄与率が９０％以上となる３つの主成分の主成分得点を第１主成分、第２主成分、第３主成分をｘ軸、ｙ軸、ｚ軸にとったグラフが図６、各主成分の要素を筋拮抗比番号順に示したグラフが図７である。
【００７１】
他項目を以下に示す。
【００７２】
[走行速度4.0km/hのとき]
Ｅｒ＝{2.16, 1.80, 3.32, 1.56, 0.89, 2.85, 1.37, 0.96}
λ＝{4.62, 2.43, 0.91}
ｃ＝{0.56, 0.29, 0.11}
Ｃ＝{0.56, 0.85, 0.96}
[走行速度3.0km/hのとき]
Ｅｒ＝{2.10, 1.45, 3.20, 2.03, 0.99, 2.45, 0.98, 0.87}
λ＝{4.03, 2.76, 0.74}
ｃ＝{0.49, 0.33, 0.09}
Ｃ＝{0.49, 0.82, 0.91}
なお、図６において、第３主成分は、値−２．０００〜２．０００までを、０．２００ずつの２０区分に分けたもので、色彩で識別表示している。また、図６中、各プロットには、ｔ1，ｔ2，…ｔ20の経時情報が付されている。今、色彩の識別表示を、最下位から最上位に向けて、識別番号１（値−２．０００〜−１．８００）から識別番号２０（値１．８００〜２．０００）の番号を付すとき、
図６（ａ）のｔ1からｔ20までの識別番号は、｛１０，９，８，８，９，１２，１５，１９，１９，１６，１１，７，３，３，５，７，９，１２，１２，９｝であり、
図６（ｂ）のｔ1からｔ20までの識別番号は、｛１１，６，７，７，７，９，１３，１８，１９，１８，１３，９，８，７，６，７，９，１１，１３，１３｝である。但し、実際には、モニタ２４はカラーモニタが採用されるため、そのままで色識別が可能である。
【００７３】
また、図８は、筋活性度を示すグラフである。図９は、筋拮抗比の主成分得点のうち、ここでは第１、２主成分における差分を表示した図である。なお、図６〜図８において、（ａ）は４．０[km/h]、（ｂ）は３．０[km/h]である。
【００７４】
表示制御部２０５は、時間情報の第１主成分の差分を算出すると共に、算出した差分を矢印で表示している。あるいは、各時間情報中の同一時間となるプロット位置とプロット位置と間に線分を表示し、かつ参照側から新たな測定側に向けて矢印を付すように描画してもよい。
【００７５】
（考察）
G.Cappelliniら（“Motor Patterns in HumanWalking and Running”, J Neurophysiol Vol.95,p.3426-3437, 2006.）の報告では、３２筋の筋電位に対し５つのタイミングパターンで筋電位データを十分説明できるとしていたが、本処理は８筋に対する８種の筋拮抗比について３つの主成分で表した際、累積寄与率は９０％を超え、３つの要素基本パターンで標準化した筋電位データを説明できることが判った。
【００７６】
次に、筋協調について検討する。行列Ｒに対する主成分解析の結果の、累積寄与率が９０％を超えるまでの各主成分の要素と、筋拮抗比の主な作用を比較すると、第１主成分は股関節伸展、第２主成分は膝関節伸展、第３主成分は足関節伸展の役割をもっていると考えられる。第１主成分、第２主成分をそれぞれｘ軸、ｙ軸にとったグラフは楕円に近い形を描くのが特徴的である。歩行相と対応させてみると、両脚支持期を除いて、ほぼ一定間隔で主成分得点が変化しているのがわかる。また、第３主成分は両脚支持期付近で値が大きくなることがわかる。これら結果から歩行動作を３つの主成分の重みを変化させることで説明できることがわかる。さらに、これらの値の正負も重要な意味を持っており、符号の反転は、作用の反転を意味していると考えられる。すなわち、ｙ軸を例にとるならば、立脚期から遊脚期に移行する際に、ｙ軸の符号が反転している。これは歩行相と対応させて考えると、作用が反転していることがわかる。
【００７７】
筋活性度では、ほぼすべての拮抗する筋ペアで踵接地付近での活動が高くなっている。これは踵接地の瞬間、外乱に備えて関節剛性を高めていると考えられ、リハビリテーション分野などで一般的に言われている、踵接地時の衝撃に対する筋活動の特徴に関連付けられる（江原義弘: “人間の歩行、ロボットの歩行”, 計測と制御、Vol45,p.1018-1023,2006.）。またボールキャッチングでは、ボールが手に触れる前後だけ、手首や肘の拮抗する筋の活性度が上がり、関節剛性を高めて衝撃に備えるという報告にも類似性が予想される（Theodore E.Milner: “Adaptation to destabilizing dynamics by means ofmuscle cocontraction”, Exp BrainRes Vol43,p.406-416,2002.）。足関節拮抗筋の活性度だけが単脚支持期の間、高まっているのも特徴である。
【００７８】
なお、本実験では、歩行動作中の健常男性右片脚８筋の筋電位を計測し、筋拮抗比に基づく主成分解析を行い、筋協調を求めた。今回行った、筋拮抗比に基づいた解析結果では、先のG.Cappelliniらの重視する筋活動のタイミングに比して、運動中の筋協調が重視された結果となったといえる。これら筋拮抗比を用いた主成分解析は、歩行評価などに応用が考えられる。また、活性度によって外乱に対して安定性を高めている可能性も示された。以上のことから、関節運動において拮抗する筋同士の働きは非常に重要であると考えられる。
【００７９】
＜実検例２＞
（実験方法）
本実験では、手先の拘束を全く受けずに、平面内で時計回りの円軌道を描くという固有のタスクについて、健常成人男性（22歳，右利き）の協力のもと、筋電位、手先位置および関節角の測定を行った（図１０）。
【００８０】
筋電位測定を行った筋の名称を表２に示す。
【００８１】
【表２】
2011030781.tif 000010
【００８２】
また、各パラメータの定義及び筋の付着点を表３に示す。
【００８３】
【表３】
2011030781.tif 000011
【００８４】
手先位置及び関節角の初期値は、
（φ，θ１，θ２，θ３，）＝（−π／２，π／２，π／２，０）とする。
【００８５】
（実験タスク）
手先は全く拘束しないが、重力を補償するために十分に高い天井からひもを吊るし、上腕を支える。目標円軌道半径は１００[mm]とし、なるべく目標円軌道を描くよう被験者に指示した。全データのｘ座標およびｙ座標それぞれの平均値を回転中心とみなし、計測後に回転中心ｃ及び位相角φを算出する。また、φ＝０を１試行の開始地点/終了地点とみなす。そして、１試行を８秒〜１０秒で行うように指示し、６０秒間で６試行分の測定データを取得した。
【００８６】
（筋電位測定）
筋電位の測定は、１０[mm]の表面電極を電極間距離２０[mm]にて使用し、筋電位取得装置及び生体アンプWEB-9000（日本光電社製）で筋電位を増幅した後、PowerLab8/30 （ADInstruments 社製）を通じて行い、サンプリング周波数１[khz]で記録した。なお、被験者に電極を張り付ける際にはスキンピュア（日本光電社製）を用い、皮膚抵抗が１０[kΩ]以下となるように処理を行った。
【００８７】
（その他の計測機器）
本実験装置（システム）には、SenSable Technologies 社製（SensAble technologies: http://www.sensable.com/.）が開発した力覚提示デバイスであるPHANTOMを導入した。ただし、本運動（タスク）では力覚提示を行わず、手先位置計測のみを行った。また、関節角の計測にはゴニオメータ（Biometrics社製）を用い、ADボード（Contec 社製）を通じて取得した電圧値を角度値に変換した。
【００８８】
（解析）
各筋の間に存在する協調関係を明らかにするため、表４で定義する筋拮抗比ｒ１〜ｒ４を求め、さらに１試行の長さを正規化するなどした上で主成分分析を行った。解析手順を以下で述べる。
【００８９】
【表４】
2011030781.tif 000012
【００９０】
（％ＭＶＣおよび筋拮抗比の算出）
筋電位で取得したデータに対して、整流化、フィルタリング（50〜150[Hz]；通過帯域)、平滑化を順に施した上で、予め測定した最大等尺性収縮時の値で除する事により正規化し、その百分率を取って％ＭＶＣに変換する。更に、表４に従い、％ＭＶＣから筋拮抗比を求める。
【００９１】
（試行の正規化）
１試行にかかる時間が均一ではないため、各試行時間について正規化を行った。さらに各試行を時間的に２５分割し、各時間区間をｔ1,ｔ2,…,ｔ25と定義する。そして、各時間区間内のデータを平均し、さらに全試行の平均値を求め、その結果、得られたデータを代表データとして扱う。
【００９２】
（筋拮抗比の主成分解析）
代表データのうち、筋拮抗比を表わす４列を抜き出す。抜き出した２５行４列のデータについて主成分解析を行い、その結果を用いて、筋拮抗比に基づく人間の運動解析を行う。
【００９３】
（結果）
主成分解析の結果を、表４に示す。なお、第１主成分得点をＰ1、同様に第２主成分得点以下をＰ2,Ｐ3,Ｐ4で表す。また、各関節角の遷移について図１１に、主成分得点の推移を図１２に、及び第１主成分得点を横軸に、第２主成分得点を縦軸にプロットした様子を図１３に示す。なお、図１３においては、５ポイント毎に時間区間をラベリングしてある。
【００９４】
（考察）
第２主成分までの累積寄与率が９５％以上になるという結果が得られた。従って、第１主成分、第２主成分の主成分得点、及び主成分ベクトルに着目して考察を行う。
【００９５】
本実験で行った運動（タスク）は、平面内をゆっくりと移動するタスクであるため、各関節角パラメータ、すなわち手先位置を維持するための筋拮抗が生じていると考えることができる。θ1とＰ1、θ2とＰ2をそれぞれ見比べてみると、主成分得点のグラフと角度のグラフの位相がほぼ合致していることからも、上で述べた考えが正しいことが確認できた。
【００９６】
次に、第１主成分および第２主成分についての考察を行う。まず、第１主成分については、φ＝π／２およびφ＝−π／２付近で主成分得点が０に近づき、φ＝π付近で最大値、φ＝０付近で最小値に近づくことから、左右方向あるいは運動基部と運動端部とを結ぶ方向に垂直な方向（偏角方向）の運動を支配している成分であることが分かる。例えば、φ＝０は、円軌道上で最も右の位置であり、肩関節が最も伸展するため、結果は妥当であると考える。
【００９７】
次に、第２主成分については、第１主成分と同様に考えると、前後方向あるいは運動基部と運動端部とを結ぶ方向（動径方向）の運動を支配している成分であることが分かる。例えばφ＝π／２は円軌道上で最も手前の位置であることから、この結果は妥当であると考える。
【００９８】
また、φ＝０付近を除けば、ほぼ真円上に主成分得点Ｐ1 Ｐ2がプロットされている。手先位置と対応付けて見ると興味深い。
【００９９】
なお、本実験は、平面内で時計回りの円軌道を描くタスクを取り上げ、筋拮抗比を用いた動作解析の一例を示した。そして、筋拮抗比を用いて筋協調の説明を行った。本実験で用いた装置（システム）では、PHANTOMを導入することによって、極めて自由度の高いタスクの実験を行うことが可能となっている。PHANTOMPremium 3.0/6DOF は、作業空間内において最大２２[N]の力を発生することができ、３次元的なフォースフィールドを生成することが可能である。また、平面内でのタスクのみならず、３次元的な運動計測も可能であり、本装置（システム）を利用した上肢の幅広い運動解析を行うことができるもの思われる。
【０１００】
なお、本発明は、以下の態様を採用することが可能である。
【０１０１】
（１）本実施形態では、腕と脚の所要の運動（タスク）を対象としたが、本発明は、これに限定されず、拮抗する筋ペアを備える関節部位における所定の運動（タスク）に適用可能である。
【０１０２】
（２）本実施形態では、歩行相判断用にカメラを用いたが、カメラ必須ではなく、例えば時間情報を利用して歩行相を設定するようにしてもよい。
【０１０３】
（３）本実施形態では、検体から得られる基準及び差分信号を、図１に示すロボット３０のアクチュエータ３１に供給することで、差分状況をモニタ２４を介して観察する場合に比して、ロボット３０を用いて実際の動きの中で、所定の運動（タスク）のコーチングをより好適に行うことが可能となる。
【０１０４】
（４）本実施形態では、各主成分得点の差分をモニタ２４に表示するようにしたが、この差分に応じたレベルの電気刺激信号を検体の所与の部位（筋）に付与することで、検体の所定の運動（タスク）を理想（基準）に近づけるような機能的電気刺激システムにも適用可能である。
【０１０５】
図１４は、本発明に係る筋シナジーインターフェースの一例を示す構成図である。図１５は、本発明に係る筋シナジーインターフェースによって実行される処理手順の概要を示すフロー図である。
【０１０６】
筋シナジーインターフェース１００は、制御対象１６０に対する運動内容を指示する運動コマンド設定部１１０と、マイクロコンピュータ等で構成される制御部１２０とを備えている。制御部１２０は、運動コマンド設定部１１０からのコマンドを主成分情報に変換する第１の変換部１３０と、第１の変換部１３０で得られた主成分情報を制御対象１６０を運動させる信号に変換する第２の変換部１４０とを備えている。そして、筋シナジーインターフェース１００は、第２の変換部１４０で得られた信号を、運動指示のための電位信号として出力する出力部１５０とを備えている。
【０１０７】
制御対象１６０は、リハビリ対象者（ヒト）であってもよい。ヒトの場合、対象筋肉に対向した皮膚表面に電極１５１，１５１を貼付し、この電極１５１に出力部１５０からの電位信号を出力するようにすればよい。本実施形態では、電位信号は、後述するように拮抗する筋対毎に生成される。また、制御対象１６０は、ヒトの動きを支援する（リハビリ用の）装具（パワースーツ等）やロボットなども想定される。装具やロボットの場合、機構体の所要部位を運動させるための複数の人工筋等のアクチュエータを有する。アクチュエータは、電気信号に応じて空圧を調整し、空気圧ゴムチューブを伸張、収縮可能にするMcKibben型人工筋肉等で構成されており、出力部１５０からの信号を電位信号（駆動信号）として受信し、動作を行う。アクチュエータはMcKibben型人工筋肉に限定されず、電気信号で、あるいは電気信号を変換した別の力によって機械的運動に変える電気機械変換素子、例えば電磁ソレノイド、圧電素子、モータ等であってもよい。
【０１０８】
運動コマンド設定部１１０は、運動コマンドである情報が入力可能な構成、例えば、キーボードやマウス等からなる操作入力部を少なくとも備えている。また、運動コマンド設定部１１０は、入力された運動を規定する運動情報を、典型的には位置情報として出力するものである。位置情報は、例えば、軸上とか平面内での運動態様を、後述するように直交座標で規定している。
【０１０９】
制御部１２０は、更にモニタ１２１を備え、入力データや処理途中のデータが確認可能な態様としてもよい。また、制御部１２０は、制御プログラムや基本データ類を記憶するＲＯＭや処理途中の情報を一時的に記憶するＲＡＭで構成される所要の記憶部１２２を有する。
【０１１０】
第１の変換部１３０は、運動コマンド設定部１１０で設定される直交座標の位置情報を極座標に変換する座標変換部１３１、及び極座標と主成分情報との関係を規定する主成分得点変換部１３２から構成されている。第２の変換部１４０は、主成分得点情報を元に主成分に変換する主成分処理部１４１、主成分を各筋拮抗比及び活性度に変換するための出力変換部１４２、及び筋拮抗比及び活性度から、拮抗する筋対を構成する主動筋と拮抗筋とに印加されるレベルの電位信号を作成する筋活動コマンド生成部１４３から構成されている。
【０１１１】
図１５において、ステップＳ１０１は、運動コマンド設定部１１０によって、例えば手先位置情報（ｘ、ｙ）が目標運動情報として設定される。このステップＳ１０１で設定される情報は、位置情報の他、運動速度や加速度が、更には運動速度が速い場合における動的干渉や、手先の動きに規制（手先力）等をかける所定の制約がある場合における付属情報が含まれてもよい。なお、運動コマンド設定部１１０で入力される目標運動情報としては、手先の運動を規定する情報の他、歩行の運動を規定する情報であってもよく、また位置情報の他、運動に応じて該運動を規定するのに好ましい情報を採用する態様としてもよい。例えば、歩行運動の場合における情報としては、歩行速度情報と歩幅情報とが採用可能である。かかる情報であっても、足先の位置を規定する情報となり得る。
【０１１２】
ステップＳ１２１は、上述した手先に対する運動情報の設定の場合であり、ステップＳ１２２は、上述した歩行における運動情報の設定の場合である。なお、歩行の基底関数とは、歩行運動を記述する基本の関数をいう。ステップＳ１０３は、主成分得点を取得する処理を行う。ステップＳ１０４は、取得した主成分得点から主成分を生成する処理を行う。ステップＳ１０５は、主成分を筋拮抗比及び活性度に変換する処理を行う。ステップＳ１６０は、拮抗する筋対毎に主動筋と拮抗筋とに印加する電位信号を生成し、対応する電極１５１に出力する処理を行う。
【０１１３】
図１５に示す処理手順は、ステップＳ１０１〜ステップＳ１０３の処理と、ステップＳ１０４〜ステップＳ１０６の処理に大別される。ステップＳ１０１〜ステップＳ１０３までの処理を第１の変換処理とし、ステップＳ１０４〜ステップＳ１０６の処理を第２の変換処理とする。
【０１１４】
まず、ステップＳ１０４〜ステップＳ１０６の第２の変換処理について説明する。この第２の変換処理は、特願２００９−２１２１４９において本発明者らによって既に提案されている。
【０１１５】
すなわち、第２の変換処理は、多変量解析を利用して得られる主成分情報から、各拮抗する筋対を構成する主動筋と拮抗筋とにそれぞれ筋活動コマンドとしての電位信号を生成し、出力する処理である。これは、拮抗する筋対の各筋へ印加される電位信号のレベルと主成分情報との間に高い相関性が存在することに基づくものである。
【０１１６】
特願２００９−２１２１４９においては、検体に所定の運動動作、例えば上肢の水平面内での円運動、下肢による歩行運動を行わせると共に、これらの運動において活動する筋の電位を筋と対向する皮膚に貼付した電極を介して検出し、検出した各筋の電位を多変量解析を介して、主成分を抽出した。例えば、上肢では、図１７に示すように、肩から手先にかけての、三角筋前部繊維（ｍ１）、・・・、尺側手根伸筋（ｍ８）の筋電位を、筋の協調関係を表す各拮抗する筋対の筋拮抗比と筋活性度とに置き換え、これらの各拮抗する筋対の筋拮抗比と筋活性度とから、多変量解析を実行して、第１、第２、第３、・・・の主成分を算出し、寄与度が所定％、好ましくは略９０％を超える主成分を抽出するものである。この上肢の円運動動作では、第１、第２の主成分で累積寄与率９０％を超えている。従って、第１、第２の主成分の得点を求め、更に主成分ベクトルを求めて、運動に寄与する力として各拮抗する筋対の主動筋と拮抗筋とに対して、高い相関を有する電位信号が生成可能となる。
【０１１７】
また、下肢による歩行運動を行わせると共に、これらの運動において活動する筋の電位を筋と対向する皮膚に貼付した電極を介して検出し、検出した各筋の電位を多変量解析を介して、主成分を抽出した。例えば、足の付け根から足先にかけて、股関節１関節筋、股関節膝関節２関節筋、膝関節１関節筋、足関節筋の各屈筋・伸筋をペアとした８筋を対象にした。検出した筋電位を各拮抗する筋対の筋拮抗比と筋活性度とに置き換え、これらの各拮抗する筋対の筋拮抗比と筋活性度とから、多変量解析を実行して第１、第２、第３、・・・の主成分を算出し、寄与度が所定％、好ましくは略９０％を超える主成分を抽出するものである。この下肢の歩行動作では、第１、第２の主成分で少なくとも累積寄与率８０数％を超えている。従って、第１、第２の主成分の得点を求め、更に主成分ベクトルを求めて、運動に寄与する力として各拮抗する筋対の主動筋と拮抗筋に対して、高い相関を有する電位信号が生成可能となる。
【０１１８】
このように、予め運動を司る筋及び運動内容（上述の例では、上肢の円運動、下肢の歩行運動）が決まると、当該運動を規定する主成分情報が決まることになる。そこで、事前に決まる主成分情報を利用して、運動内容から主成分得点を得る処理を含む第１の変換処理（図１５のステップＳ１０１〜Ｓ１０３）が可能となる。以下、説明する。
【０１１９】
まず、運動内容が決まれば、対応する主成分の種類が決まる。例えば、上肢の水平面内運動であれば、第１、第２の主成分が所要の累積寄与率で予め決定される。そして、後述するように、運動コマンドと主成分との間に何等かの相関があれば、当該相関を利用することで、運動コマンドを解析して各筋に対する電位信号を生成するという複雑な演算を実行することなく、図１５のステップＳ１０１からステップＳ１０６まで、主成分情報を利用した一連の、すなわち簡易、迅速に電位信号を生成する処理が可能となる。
【０１２０】
そこで、まず、運動コマンドの内容と当該運動に対応した主成分との関係を説明する。運動コマンドは、大脳運動野で生成される、運動内容を記述する目標位置（ｘ、ｙ）として定義される。前述の上肢による水平面内での運動であれば、目標手先位置（ｘ、ｙ）である。
【０１２１】
ここで、上肢による上記水平面内運動を分析すると、肩と手先を結ぶ直線上に沿った伸縮（距離ｒ）と、肩あるいは頭部を中心とした手先の回転（角度θ）とで記述し得る。到達運動タスクやクランク回転タスクにおける検出筋電位の主成分解析の結果、多くの事例で運動コマンドはかかる２種類の意味を持つことが判っている。
【０１２２】
すなわち、手先の運動コマンドを極座標（ｒ、θ）で表現していると考えられる。大脳運動野では目標とする手先の運動が直交座標系で表現されるため、中枢・脊髄神経系で極座標への座標変換が行われていると考えられる。手先運動情報として目標手先位置の他、速度、加速度、またダイナミックな干渉や手先が環境に及ぼす力やモーメント等の付属的な情報も考慮すると、上肢の水平面内運動における運動コマンド生成の流れは、図１６のようになる。
【０１２３】
図１６は、上肢の水平面内運動における運動コマンド生成の流れを示すフロー図である。図１６において、まず、ステップＳ２０１で、運動コマンド設定部１１０によって目標手先位置（ｘ、ｙ）が設定される。目標位置の設定は、座標を直接入力する態様でもよいし、予め運動内容（座標情報）と識別情報とを関連付けて記憶しておき、識別情報を指定することで対応する運動内容が指定される態様としてもよい。この場合、運動内容は、制御部１２０内の記憶部１２２に記憶されていればよい。運動内容には、単に位置情報の他、運動速度等もあることから、必要な付属的な情報も指定可能な態様がよい。
【０１２４】
ステップＳ２０２、Ｓ２０３で、目標手先位置（ｘ、ｙ）は、公知の直交座標−極座標変換式を使用する等して、極座標（ｒ、θ）に変換される。ステップＳ２０４で、付属事項も含めて、目標手先位置（ｒ、θ）＋α（速度・加速度・動的干渉・手先力・姿勢）が得られる。
【０１２５】
そして、ステップＳ２０５で、主成分毎の得点（ｗ１，ｗ２）がそれぞれ算出される。具体的には、
ｗ１＝ｆ（ｒ、ｒ′、ｒ″）＋α１（動的干渉、手先力、姿勢）
ｗ２＝ｇ（θ、θ′、θ″）＋α２（動的干渉、手先力、姿勢）
と表される。なお、記号「′」は１次微分を示し、記号「″」は２次微分を示す。上記において、二次微分項は必ずしも存在するとは限らず、逆に、場合によっては３次微分項が必要となる態様も考えられる。ここでは、主成分Ｗ１が極座標の距離ｒに対応し、主成分ｗ２が極座標の角度θに対応している。また、運動コマンドが極座標（ｒ、θ）で表現されていることは、歩行運動中の下肢においても同様であり、腰に対する足先位置を極座標（ｒ、θ）で表現したものが、（ｗ１，ｗ２）に対応している。
【０１２６】
続いて、実験結果を説明する。表５は、本実験における筋電位計測の対象となる８つの筋の名前と、対応する電極貼付位置とを示す。
【０１２７】
【表５】
2011030781.tif 000013
【０１２８】
また、各筋拮抗比の定義を以下（表６）に示す。
【０１２９】
【表６】
2011030781.tif 000014
【０１３０】
図１７は、被験者７の頭部７０、肩から手先までの各筋の位置関係と、各拮抗する筋対の関係を示す筋の模式図である。
【０１３１】
＜実験タスク＞
平面内で手先位置を半径１００[mm]の目標円軌道に沿って移動させるタスクを行い、カメラ等による手先位置の計測と筋電位計測とを行った。なお、重力を補償するために十分に高い天井から紐を吊るし、上腕を支えた。
【０１３２】
＜実験結果１＞
時計回りの運動における主成分ベクトル（２つの主成分の累積寄与率９０％以上）と主成分得点との関係は、表７及び図１８の通りである。
【０１３３】
本実験では、筋電位の測定は、例えば１０[mm]の表面電極を電極間距離２０[mm]にて使用し、筋電位取得装置及び生体アンプWEB-9000（日本光電社製）で筋電位を増幅した後、PowerLab8/30 （ADInstruments 社製）を通じて行い、サンプリング周波数１[khz]で記録した。なお、被験者に電極１５１を張り付ける際にはスキンピュア（日本光電社製）を用い、皮膚抵抗が１０[kΩ]以下となるように処理を行った。また、本実験装置（システム）には、SenSable Technologies 社製（SensAble technologies: http://www.sensable.com/.）が開発した力覚提示デバイスであるPHANTOMを導入した。ただし、本運動（タスク）では力覚提示を行わず、手先位置計測のみを行った。また、関節角の計測にはゴニオメータ（Biometrics社製）を用い、ADボード（Contec 社製）を通じて取得した電圧値を角度値に変換した。各筋の間に存在する協調関係を明らかにするため、表６で定義する筋拮抗比ｒ１〜ｒ４を求め、さらに１試行の長さを正規化するなどした上で主成分分析を行った。さらに、筋電位で取得したデータに対して、整流化、フィルタリング（50〜150[Hz]；通過帯域)、平滑化を順に施した上で、予め測定した最大等尺性収縮時の値で除する事により正規化し、その百分率を取って％ＭＶＣ（% Maximum Voluntary Contraction）に変換した。そして、表６に従い、％ＭＶＣから筋拮抗比を求めた。
【０１３４】
【表７】
2011030781.tif 000015
【０１３５】
表７から、第１主成分ベクトルの第１，３要素（0.5807，-0.5635）は、肩伸展(＋)と肘屈曲(−)を表しており、第1主成分が肩と手先を結ぶ直線上に沿った伸縮（ｒ）に対応していることが判る。一方、第２主成分ベクトルの第２要素（0.9472）は、２関節筋による肩と肘の同時屈曲(＋)を表しており、第２主成分が肩あるいは頭部を中心とした手先の回転（θ）に対応していることが判る。図１８は、手先を水平面内で時計回りの円運動をさせたとき、横軸に第１主成分ｗ１を、縦軸に第２主成分ｗ２を採って、（ｗ１，ｗ２）の点の運動位置（時間経過；図中の矢印参照）を表した図表である。
【０１３６】
また、反時計回りの運動における主成分ベクトルについても同様であることを確認した。すなわち、反時計回りの運動における主成分ベクトル（２つの主成分の累積寄与率９０％以上）と主成分得点との関係は、表８及び図１９の通りである。
【０１３７】
【表８】
2011030781.tif 000016
【０１３８】
表８から、第１主成分ベクトルの第１，３要素（0.6216，-0.4704）は、肩伸展(＋)と肘屈曲(−)を表しており、第1主成分が肩と手先を結ぶ直線上に沿った伸縮（ｒ）に対応していることが判る。一方、第２主成分ベクトルの第２要素（0.8141）は、２関節筋による肩と肘の同時屈曲(＋)を表しており、第２主成分が肩あるいは頭部を中心とした手先の回転（θ）に対応していることが判る。図１９は、手先を水平面内で反時計回りの円運動をさせたとき、横軸に第１主成分ｗ１を、縦軸に第２主成分ｗ２を採って、（ｗ１，ｗ２）の点の運動位置（時間経過；図中の矢印参照）を表した図表である。
【０１３９】
＜実験結果２＞
図２０は、第１主成分得点ｗ１が肩と手先とを結ぶ直線上に沿った伸縮（ｒ）に対応することを示す重回帰分析結果の図表である。横軸は運動位置（時間）で、縦軸はｗ１であり、破線は実験に基づくもの、実線は重回帰分析したものである。クランクタスク（時計回り）におけるｗ１の（ｒ、ｒ′）による推定結果を、式７で示す。
【０１４０】
ｗ１＝−０．０２ｒ−０．０５ｒ′＋１２．６  ・・・（式７）
図２１は、第２主成分得点ｗ２が肩あるいは頭部を中心とした手先の回転（θ）に対応することを示す重回帰分析結果の図表である。横軸は運動位置（時間）で、縦軸はｗ２であり、破線は実験に基づくもの、実線は重回帰分析したものである。クランクタスク（時計回り）におけるｗ２の（θ、θ′）による推定結果を、式８で示す。
【０１４１】
ｗ２＝６．５θ＋７．７θ′−１５．３  ・・・（式８）
図２０、図２１に示すように、重回帰分析された曲線ｗ１、ｗ２はいずれも計測値から多変量解析された主成分得点と比べて近似し、高い相関関係にあることが判る。
【０１４２】
＜実験結果３＞
図２２は、クランクタスク（時計回り）における実験値（ｒ、θ）から時計回りの（ｗ１，ｗ２）を推定した結果（上記重回帰モデルを利用）を示す図表である。図２３は、クランクタスク（時計回り）における実験値（ｒ、θ）から反時計回りの（ｗ１，ｗ２）を推定した結果（上記重回帰モデルを利用）を示す図表である。図２２、図２３において、破線は実験に基づくもの、実線は重回帰分析したものである。
【０１４３】
図２３から判るように、重回帰分析された曲線ｗ１、ｗ２はいずれも計測値から多変量解析された主成分得点と比べて近似し、高い相関関係にあることが判る。従って、反時計回りの筋電位計測をすることなく、運動を生成するための筋電位信号を合成することができる。
【０１４４】
以上のように、種々の運動内容について予め主成分解析が施されて主成分が規定されている場合、座標変換部１３１によって直交座標系の運動コマンドが極座標に変換され、更に必要に応じて付属情報も含めて、主成分得点変換部１３２によって略一意的に（高い相関性を有して）、主成分得点に変換される（図１５のステップＳ１０１〜Ｓ１０３）。
【０１４５】
続いて、主成分処理部１４１によって主成分得点が主成分に変換される処理が行われ、出力変換部１４２によって、各拮抗する筋対に対して筋拮抗比及び筋活性度が算出される。そして、筋活動コマンド生成部１４３によって、得られた各拮抗する筋対の筋拮抗比及び筋活性度から、主動筋と拮抗筋とに印加される電位信号が生成され、制御対象１６０がヒトの場合、対応する電極１５１に出力される。筋協調として、筋拮抗比及び筋活性度を利用することで、制御信号をより数少なくすることができる。
【０１４６】
なお、本発明は、以下の態様を採用することが可能である。
【０１４７】
（１）制御対象１６０が、人工筋アームを備えた筋骨格ロボット、例えば上肢筋骨格ロボットである場合、人体の場合と同様に目標手先運動を人工筋活動コマンドに変換するメカニズムが必要となる。また、人体の動きを検知する筋電計及び上述のようにして筋活動コマンドを生成する制御部を備えたマスタ側とマスタ側からの出力信号に対応して運動を行うスレーブ機側であるロボットを制御する場合（いわゆるマスタースレーブシステム）には、ヒトとロボットとでは筋特性や構造が異なるため、ヒトの筋電位計測結果からリアルタイムで得られる筋拮抗比や筋活性度をロボットの筋拮抗比や筋活性度に変換する必要がある。
【０１４８】
図２４は、筋拮抗比の変換の一例を示すフロー図である。筋活性度についても同様である。なお、変換関数ｆ、ｇは予め例えば実験結果等から獲得しておけばよい。あるいは、筋拮抗比と手先位置との関係を実験的に求め、手先位置の肩に対する動径方向の変位が第１主成分得点、偏角方向の変位が第１主成分得点に対応するよう主成分ベクトルを決定してロボットの筋拮抗比に変換する方法もある。
【０１４９】
この図２４のフロー図は、例えば図１５のステップＳ１０５の段階で実行される。まず、ステップＳ３０１で、筋拮抗比ｒｈ（ヒト）（なお、ｒｈは前述のｒｔｐに相当）が生成される。次いで、ステップＳ３０２で、目標関節角＝ｆ（ｒｈ）の設定が行われる。さらに、ステップＳ３０３で、筋拮抗比ｒｈ（ロボット）＝ｇ（目標関節角）の変換が行われる。そして、図１５のステップＳ１０６に代えて、ステップＳ３０４で人工筋活動コマンドが出力される。人工筋活動コマンドは、ロボットの人工筋アーム（アクチュエータ）に対応する駆動信号として出力される。このようにすることで、前記マスタースレーブシステムへの適用が可能となる。
【０１５０】
図２５は、人工筋アームの駆動を説明する図表である。ヒトと同様の筋骨格構造を有するロボットでは、ヒトが用いている主成分ベクトルを人工筋アームに適用することで、ヒトの運動コマンド（ｗ１，ｗ２）をそのまま利用して人工筋アームを操作できる可能性がある。図２５、図２６は、ヒトと同様の筋骨格構造を有するロボットにヒトが用いている主成分ベクトル埋め込み、人工筋の制御コマンドを計算した後、人工筋アームを駆動した実験結果を示している。
【０１５１】
図２５は、主成分得点ｗ１，ｗ２を離散的に与え、人工筋の制御コマンドを計算した後、人工筋アームを駆動した場合のものである。図２５（ａ）の実験１（ｗ１）はそのときの手先位置をプロットしたもので、グラフの矢印は主成分得点ｗ１が離散的に大きくなる方向を表している。同様に、矢印に示すように、主成分得点ｗ２を離散的に大きくしていった場合の手先位置をプロットしたものが実験２（ｗ２）である。
【０１５２】
図２５（ｂ）において、実験３（ｗ１，ｗ２）は、ｗ１，ｗ２を同時に変えたときの結果を示している。図２５（ｂ）によれば、（ｗ１，ｗ２）を適当に与えることによって、水平面内の任意の位置（可動範囲内の）に手先を位置決めできることが判る。人工筋アームの側もヒトと同様の制御部を実装しておけば、ヒトの運動コマンド（ｗ１，ｗ２）をそのまま利用して人工筋アームを操作することが可能である。
【０１５３】
（２）本実施形態では、図１５のステップＳ１０１で設定する運動コマンドとして直交座標を採用したが、運動を規定する運動情報としては、これに限定されず、例えば歩行運動では、運動を規定できれば、他の情報、例えば歩幅と周期であってもよい。そして、この態様の場合、図１５のステップＳ１０１〜Ｓ１０３に代えて、以下の式９、式１０で表される方法が採用可能である。
【０１５４】
【数７】
2011030781.tif 000017
【０１５５】
式９は、歩行動作における主成分得点運動の特徴を、級数展開で表した数式である。式１０は、主成分得点ｗ１及びｗ２が周期Ｔ，Ｔ／２のみの級数展開式で記述できるとした場合、ｗ１，ｗ２が持っているある性質を利用するようにして表現したものである。ある性質とは、ｗ１，ｗ２を半周期だけ位相をずらし、それぞれの和と差を取ると、必ず正弦波様や余弦波様（正弦波や余弦波に近似した波形）が得られるというものである。そして、その周期は和を取る場合には歩行周期Tに対してＴ／２になり、一方、差をとるとＴになる。
【０１５６】
この性質を満たすための、ｗ１，ｗ２の式の形は、逆算を行うことで、式９のｗ１（ｔ）、ｗ２（ｔ）の級数展開の式を得ることができる。なお、数式中の変数は、w1：第1主成分得点、w2：第2主成分得点、Ｔ：１歩行周期、ａ１，ａ２，ａ３，ａ４：第１主成分得点ｗ１の級数展開の係数、ｂ１，ｂ２，ｂ３，ｂ４：第２主成分得点ｗ２の級数展開の係数である。
【０１５７】
図２６は、級数展開における近似結果を示す図表で、図２６（ａ）は歩行（３km／h）を分解した図、図２６（ｂ）（ｃ）はその時の主成分ｗ１，ｗ２を示している。図２６（ｂ）（ｃ）において、各点は計測結果から得られた主成分得点、実線は本級数展開の数式による近似結果であり、両者は極めて近似していることが認められる。
【０１５８】
また、式１０は、歩行動作における主成分得点運動の特徴を表すためのダフィング方程式である。
【０１５９】
【数８】
2011030781.tif 000018
【０１６０】
ヒトの歩行現象を説明する場合、生理学分野の長年の研究により、ＣＰＧ(CentralPattern Generator）と呼ばれる歩行パターンを発生する仕組みが中枢下位にあることが示唆されている。そして、かかるパターン生成器の挙動を非線形振動子を用いて説明する研究が行われている。また、振動子を定義することで運動指令の挙動を表す力学場が定義できる。一方、脳科学の分野においても脳内に力学場を仮定して現象を説明しようとする試みがあるが、同様に運動指令そのものである主成分得点運動を支配する力学場を見つけることが必要となり、そのためにダフィング方程式を採用した。この方程式自体は様々な工学分野に出てくる一般的なモデルである。この方程式で近似して得られる解を半周期分位相をずらして和や差を取ると、正弦波や余弦波に似た波形を示し、式９で示した級数展開法の動機となった現象を満足することとなる。外力が印加されたダフィング方程式は、パラメータの決定が難しい面があるものの、さまざまな歩行動作を、この式で説明することが可能である。すなわち、パラメータが決まれば、さまざまな歩行運動を記述可能できる。なお、この式におけるパラメータは、主成分得点運動を機械モデルと表した場合、それぞれの変数は、ｃ：粘性係数、ｋ：弾性係数、Ａ：周期外力（例えば重心運動）の振幅、Ｔ：歩行周期、ｚ：定常外力（例えば上位中枢からの速度指令）である。また、ダフィング式は、変数Ａ，ｚが他のパラメータに比べて大きく影響することから、かかる変数Ａ，ｚ、いわゆる歩幅と歩行周期とを適宜に調整設定することで、所要の運動が記述可能となる。
【０１６１】
図２７は、ダフィング方程式による近似結果を示す図表で、図２７（ａ）は歩行（３km／h）を分解した図、図２７（ｂ）はその時の主成分ｗ１に関する情報を示し、図２７（ｃ）はその時の主成分ｗ２に関する情報を示している。図２７（ｂ）（ｃ）において、破線は計測結果から得られた主成分得点、実線はダフィング方程式による近似結果であり、両者は極めて近似していることが認められる。
【０１６２】
また、ダフィング方程式においては、パラメータを変更することによって、図２７の３km／hとは異なる別の歩行速度、例えば４km／hでの主成分得点を算出することが可能である。この場合も、実測を行って、その結果から得られる主成分得点と、ダフィング方程式のうち、速度に寄与するパラメータｚの設定によって推定した４km／hでの主成分得点とが近似していることが認められている。
【０１６３】
（３）本実施形態では、２次元の運動について説明したが、本発明は、前述したように、軸上での１次元の運動の他、３次元の運動にも適用可能である。例えば３次元の運動の場合、直交座標から３次元極座標（ｒ、θ、φ）（φ：俯仰角）に変換され、各極座標が累積寄与度に依存する数の主成分に置換することで上述と同様に処理可能である。
【０１６４】
以上、説明したように、本発明に係る筋シナジー解析方法は、検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出するステップと、前記検体の所定の運動中に検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出するステップと、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出するステップとを備えることが好ましい。これによれば、各筋の筋電位行列に対して、単に多変量解析（ここでは主成分解析）を施すのではなく、運動における筋の協調を考慮して、筋拮抗比という概念を導入し、この概念に対して、主成分解析を施すことで、変数の種類を低減させて処理の容易、高速化を実現して、リアルタイム処理が可能となる。また、煩雑かつ困難な身体各部位のパラメータ同定が不要であるため、その分、処理が簡易となる。また、主成分解析後の変数の種類が、第３主成分までであると、表示に際して３次元座標系で表現できるという利点がある。しかも、筋拮抗比を採用することで、運動の種類毎に主成分解析の結果も異なり、運動に即した主成分の数、得点が得られる。
【０１６５】
また、本発明に係る筋シナジー解析方法において、前記検体の筋拮抗比に基づいて得られた主成分を参照主成分とし、この参照主成分と、新たな検査によって算出された同検体の筋拮抗比に基づいて得られた主成分との差分を算出するステップと、算出された前記差分を第１の表示部に表示するステップとを備えることが好ましい。この構成によれば、検体が同一人物である場合には、事前に取得しておいた自己の基準となる参照運動状態との対比が差分という形で示されることから、基準状態に近づけるアドバイスが的確に得られる。
【０１６６】
また、本発明に係る筋シナジー解析方法において、前記検体の運動中に検出された筋電信号から前記拮抗する筋ペアの筋電信号のレベルを加算して筋活性度の時系列データを算出するステップと、前記時系列データの各時刻における前記筋活性度を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出するステップとを備えることが好ましい。この構成によれば、筋拮抗比に対する多変量解析に加えて、好ましくは、拮抗する筋ペアについて筋活性度という新たな概念を導入し、この成分に対して、多変量解析を施したので、上記筋拮抗比と併用することで、運動の推定、評価がより正確となる。
【０１６７】
また、本発明に係る筋シナジー解析方法において、前記検体の筋活性度に基づいて得られた主成分を参照主成分とし、この参照主成分と、新たな検査によって算出された同検体の筋活性度に基づいて得られた主成分との差分を算出するステップと、算出された前記差分を第２の表示部に表示するステップとを備えることが好ましい。この構成によれば、事前に取得しておいた自己の基準となる参照運動状態との対比が差分という形で示されることから、基準状態に近づけるアドバイスが的確に得られる。
【０１６８】
また、本発明に係る筋シナジー解析装置は、検体の拮抗する筋ペア群における各筋に生じる筋電信号を検出する検出器と、前記検体の所定の運動中に前記検出器によって検出された筋電信号から拮抗する筋ペアの筋電信号のレベルの筋拮抗比を時系列データとして算出する筋拮抗比算出手段と、前記時系列データの各時刻における前記筋拮抗比を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第１の主成分解析手段とを備えることが好ましい。本発明によれば、筋シナジー解析方法と同様な作用を有する。
【０１６９】
また、本発明に係る筋シナジー解析装置において、第１の検体と第２の検体における両主成分の差分を算出する第１の差分算出手段と、前記第１の差分算出手段によって算出された差分を第１の表示部に表示する第１の表示制御手段とを備えることが好ましい。この構成によれば、第１の検体と第２の検体の差分が具体的に表示されるので、対比が的確、かつ容易となる。例えば、体格などが近似した者同士が同一の運動をした場合に、互いの差異を具体的に知ることが可能となる。
【０１７０】
また、本発明に係る筋シナジー解析装置において、前記第１の検体と前記第２の検体とは同一体であることが好ましい。この構成によれば、第１、第２の検体が同一人物である場合には、事前に取得しておいた自己の基準となる参照運動状態との対比が差分という形で示されることから、基準状態に近づけるアドバイスが的確に得られる。
【０１７１】
また、本発明に係る筋シナジー解析装置において、各主成分の寄与度の累積値が第１の閾値を超えていることが好ましい。この構成によれば、多変量解析で得られた主成分の種類は、寄与度の累積値が第１の閾値を超えたものであるので、分析結果を信頼性のあるものとすることができる。
【０１７２】
また、本発明に係る筋シナジー解析装置において、前記検体の運動中に前記検出器によって検出された筋電信号から前記拮抗する筋ペアの筋電信号のレベルを加算して筋活性度の時系列データを算出する筋活性度算出手段と、前記時系列データの各時刻における前記筋活性度を変数として多変量解析を行って相関の寄与度の高い少なくとも第１主成分を含む主成分を算出する第２の主成分解析手段とを備えることが好ましい。この構成によれば、筋拮抗比に対する多変量解析に加えて、好ましくは、拮抗する筋ペアについて筋活性度という新たな概念を導入し、この成分に対して、多変量解析を施したので、上記筋拮抗比と併用することで、運動の推定、評価がより正確となる。
【０１７３】
また、本発明に係る筋シナジー解析装置において、第１の検体と第２の検体における両主成分の差分を算出する第２の差分算出手段と、前記第２の差分算出手段によって算出された差分を第２の表示部に表示する第２の表示手段とを備えることが好ましい。この構成によれば、第１の検体と第２の検体の差分が具体的に表示されるので、対比が的確、かつ容易となる。例えば、体格などが近似した者同士が同一の運動をした場合に、互いの差異を具体的に知ることが可能となる。
【０１７４】
また、本発明に係る筋シナジー解析装置において、各主成分の寄与度の累積値が第２の閾値を超えていることが好ましい。この構成によれば、多変量解析で得られた主成分の種類は、寄与度の累積値が第１の閾値を超えているものであるので、分析結果を信頼性のあるものとすることができる。
【０１７５】
また、本発明に係る筋シナジーインターフェースは、複数の筋に所与の運動を行わせる筋シナジーインターフェースにおいて、前記複数の拮抗する筋対群に対する運動を規定する運動情報を入力する運動コマンド設定手段と、前記運動情報を主成分情報に変換する第１の変換手段と、前記主成分情報を前記各拮抗する筋対へ印加する電位信号に変換する第２の変換手段とを備えることが好ましい。この発明によれば、各拮抗する筋対への運動指示が、高い相関を有するより数少ない主成分情報で行うことができるため、従来に比して処理が高速となり、かつ構成がコンパクトになる。
【０１７６】
また、本発明に係る筋シナジーインターフェースにおいて、前記第１の変換手段は、運動情報としての直交座標系の位置情報を極座標系に変換するものであることが好ましい。
【０１７７】
また、本発明に係る筋シナジーインターフェースにおいて、前記第１の変換手段は、運動情報としての運動周期及び１周期での単位運動を極座標系に変換することが好ましい。この構成によれば、直交座標系の位置情報以外の情報も適用可能となるので、制御対象に合わせた高い汎用性が得られる。
【０１７８】
また、本発明に係る筋シナジーインターフェースにおいて、前記第１の変換手段は、前記運動情報を２種類の主成分情報に変換し、前記第２の変換手段は、得られた主成分情報を前記各拮抗する筋対への電位信号に変換することが好ましい。この構成によれば、主成分情報を寄与率の高い２種類とし、この２種類の主成分情報が少なくとも各拮抗する筋対への電位信号に変換されるので、相関性の高い運動が可能となる。
【０１７９】
また、本発明に係る筋シナジーインターフェースにおいて、前記第１の変換手段は、少なくとも前記複数の拮抗する筋群から構成される運動基部と運動端部との間の距離と回転角度とを２種類の主成分情報に置換し、前記第２の変換手段は、得られた２種類の主成分情報を前記各拮抗する筋対への電位信号に変換することが好ましい。この構成によれば、２種類の主成分情報が距離と回転角度とに置換される。
【０１８０】
また、本発明に係る筋シナジーインターフェースにおいて、前記第１の変換手段は、前記運動実行時の制約に応じた固有値を算出することが好ましい。この構成によれば、運動の速度、加速度、動的干渉、運動先端への運動への制約内容等の付属情報に応じて固有値が設定されるため、これらの付属情報に対応した主成分得点が設定される。
【０１８１】
また、本発明に係る筋シナジーインターフェースにおいて、前記第２の変換手段は、前記主成分情報を筋拮抗比及び活性度に変換するものであることが好ましい。この構成によれば、ヒトやヒトを模擬した人工筋を有するロボット等の制御対象への適用が可能となる。
【符号の説明】
【０１８２】
１  筋シナジー解析装置
１０  筋電計（検出器）
２０  信号処理部
２０１  測定処理部
２０２  筋拮抗比算出部（筋拮抗比算出手段）
２０３  第１の主成分解析部（第１の主成分解析手段）
２０４  第１の差分算出部（第１の差分算出手段）
２０５  表示制御部（第１の表示制御手段、第２の表示制御手段）
２０６  筋活性度算出部（筋活性度算出手段）
２０７  第２の主成分解析部（第２の主成分解析手段）
２０８  第２の差分算出部（第２１の差分算出手段）
２４  モニタ（第１の表示部、第２の表示部）
１００  筋シナジーインターフェース
１１０  運動コメント設定部
１２０  制御部
１３０  第１の変換部
１３１  座標変換部
１３２主成分得点変換部
１４０  第２の変換部
１４１  主成分処理部
１４２  出力変換部
１４３  筋活動コマンド生成部
１５０  出力部
１５１  電極
１６０  制御対象
【図１】
2011030781.tif 000019
【図２】
2011030781.tif 000020
【図３】
2011030781.tif 000021
【図４】
2011030781.tif 000022
【図５】
2011030781.tif 000023
【図６】
2011030781.tif 000024
【図７】
2011030781.tif 000025
【図８】
2011030781.tif 000026
【図９】
2011030781.tif 000027
【図１０】
2011030781.tif 000028
【図１１】
2011030781.tif 000029
【図１２】
2011030781.tif 000030
【図１３】
2011030781.tif 000031
【図１４】
2011030781.tif 000032
【図１５】
2011030781.tif 000033
【図１６】
2011030781.tif 000034
【図１７】
2011030781.tif 000035
【図１８】
2011030781.tif 000036
【図１９】
2011030781.tif 000037
【図２０】
2011030781.tif 000038
【図２１】
2011030781.tif 000039
【図２２】
2011030781.tif 000040
【図２３】
2011030781.tif 000041
【図２４】
2011030781.tif 000042
【図２５】
2011030781.tif 000043
【図２６】
2011030781.tif 000044
【図２７】
2011030781.tif 000045
【国際調査報告】
2011030781.tif 000046
2011030781.tif 000047
2011030781.tif 000048
2011030781.tif 000049