(57)【特許請求の範囲】
【請求項１】
ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mから、n個の単位縦ベクトルW(1)，W(2)，...，W(n)を横に並べた筋シナジ行列Wと、n個の横ベクトルC(1)，C(2)，...，C(n)を縦に並べた制御行列Cと、誤差行列Eと、を、
M = WC + E
を満たすように計算するシナジ計算部、
前記計算された筋シナジ行列Wに含まれる単位縦ベクトルW(1)，W(2)，...，W(n)のまとまりを表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えることを特徴とする訓練装置。
【請求項２】
請求項１に記載の訓練装置であって、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、縦ベクトルx,yの内積演算p(x,y)により、
SCI = 2/〔n(n+2)〕×Σi=1n Σj=1n,j≠i p(W(i),W(j))
と計算されるSCIである
ことを特徴とする訓練装置。
【請求項３】
請求項２に記載の訓練装置であって、
前記SCIが所定の閾値以上であることが、前記被験者に出力される
ことを特徴とする訓練装置。
【請求項４】
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算部、
前記計算された筋シナジ行列W1，W2，...，Wpに含まれる単位縦ベクトルW1(1)，W2(1)，...，Wp(1)のまとまり、単位縦ベクトルW1(2)，W2(2)，...，Wp(2)のまとまり、...、単位縦ベクトルW1(n)，W2(n)，...，Wp(n)のまとまりを表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えることを特徴とする訓練装置。
【請求項５】
請求項４に記載の訓練装置であって、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、ベクトルx,yの相関係数演算r(x,y)により、
SSI = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Wk(i),Wh(i))
と計算されるSSIである
ことを特徴とする訓練装置。
【請求項６】
請求項５に記載の訓練装置であって、
前記SSIが所定の閾値以上であること、もしくは、前記pの増加に対して前記SSIが増加傾向にあることが、前記被験者に出力される
ことを特徴とする訓練装置。
【請求項７】
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算部、
前記計算された制御行列C1，C2，...，Cpに含まれる横ベクトルC1(1)，C2(1)，...，Cp(1)のまとまり、横ベクトルC1(2)，C2(2)，...，Cp(2)のまとまり、横ベクトルC1(n)，C2(n)，...，Cp(n)のまとまり、を表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えることを特徴とする訓練装置。
【請求項８】
請求項７に記載の訓練装置であって、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、ベクトルx,yの相関係数演算r(x,y)により、
SSIC = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Ck(i),Ch(i))
と計算されるSSICである
ことを特徴とする訓練装置。
【請求項９】
請求項８に記載の訓練装置であって、
前記SSICが所定の閾値以下であること、もしくは、前記pの増加に対して前記SSICが減少傾向にあることが、前記被験者に出力される
ことを特徴とする訓練装置。
【請求項１０】
請求項１、２、４、５、７、８のいずれか１項に記載の訓練装置において、
前記被験者の右半身の筋肉に対する右特徴量と、前記被験者の左半身の筋肉に対する左特徴量と、が計算され、
前記右特徴量および前記左特徴量の両者、もしくは、前記右特徴量および前記左特徴量の類似度が、前記被験者に出力される
ことを特徴とする訓練装置。
【請求項１１】
ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mから、n個の単位縦ベクトルW(1)，W(2)，...，W(n)を横に並べた筋シナジ行列Wと、n個の横ベクトルC(1)，C(2)，...，C(n)を縦に並べた制御行列Cと、誤差行列Eと、を、
M = WC + E
を満たすように計算するシナジ計算工程、
前記計算された筋シナジ行列Wに含まれる単位縦ベクトルW(1)，W(2)，...，W(n)のまとまりを表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えることを特徴とする訓練方法。
【請求項１２】
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算工程、
前記計算された筋シナジ行列W1，W2，...，Wpに含まれる単位縦ベクトルW1(1)，W2(1)，...，Wp(1)のまとまり、単位縦ベクトルW1(2)，W2(2)，...，Wp(2)のまとまり、...、単位縦ベクトルW1(n)，W2(n)，...，Wp(n)のまとまりを表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えることを特徴とする訓練方法。
【請求項１３】
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算工程、
前記計算された制御行列C1，C2，...，Cpに含まれる横ベクトルC1(1)，C2(1)，...，Cp(1)のまとまり、横ベクトルC1(2)，C2(2)，...，Cp(2)のまとまり、横ベクトルC1(n)，C2(n)，...，Cp(n)のまとまり、を表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えることを特徴とする訓練方法。
【請求項１４】
コンピュータを、請求項１から１０のいずれか１項に記載の訓練装置の各部として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被験者が１つの振舞いをしたり、１つ以上の振舞いからなるタスクをしたりする間の筋電位の時系列から、当該被験者の当該振舞いや当該タスクに対する熟練度を判定するのに好適な訓練装置、訓練方法、ならびに、当該訓練装置をコンピュータにより実行するためのプログラムに関する。
【背景技術】
【０００２】
従来から、複数の筋肉を制御する中枢神経系の働きについての研究が進められている。たとえば、非特許文献１に開示されるように、人間が１つの振舞いをしたり、１つ以上の振舞いからなるタスクをしたりする際には、複数の筋肉が冗長性を持ちながら協働して働く筋シナジという現象が生じる。
【先行技術文献】
【特許文献】
【０００３】
【非特許文献１】Vincent C.K.Cheung, Andrea Turolla, Michela Agostini, Stefano Silvoni, Caoimhe Bennis, Patrick Kasi, Sabrina Paganoni, Paolo Bonato，and Emilio Bizzi, Muscle synergy patterns as physilogical markers of motor cortical damage, PNAS, vol.109，no.36，pp.14652-14656, www.pnas.org/cgi/doi/10.1073/pnas.1212056109,２０１２年９月４日
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ここで、怪我や病気によりリハビリが必要となった患者が、リハビリにおいて課せられた振舞いやタスクをどの程度できるようになったのか、を客観的に判断する特徴量が求められている。この特徴量を求めるにあたっては、人間の中枢神経系における筋シナジという現象を考慮する必要がある。
【０００５】
本発明は、以上のような課題を解決するためのもので、被験者が１つの振舞いをしたり、１つ以上の振舞いからなるタスクをしたりする間の筋電位の時系列から、当該被験者の当該振舞いや当該タスクに対する熟練度を判定するのに好適な訓練装置、訓練方法、ならびに、当該訓練装置をコンピュータにより実行するためのプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の第１の観点に係る訓練装置は、
ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mから、n個の単位縦ベクトルW(1)，W(2)，...，W(n)を横に並べた筋シナジ行列Wと、n個の横ベクトルC(1)，C(2)，...，C(n)を縦に並べた制御行列Cと、誤差行列Eと、を、
M = WC + E
を満たすように計算するシナジ計算部、
前記計算された筋シナジ行列Wに含まれる単位縦ベクトルW(1)，W(2)，...，W(n)のまとまりを表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えるように構成する。
【０００７】
また、本発明の訓練装置において、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、縦ベクトルx,yの内積演算p(x,y)により、
SCI = 2/〔n(n+2)〕×Σi=1n Σj=1n,j≠i p(W(i),W(j))
と計算されるSCIである
ように構成することができる。
【０００８】
また、本発明の訓練装置において、
前記SCIが所定の閾値以上であることが、前記被験者に出力される
ように構成することができる。
【０００９】
本発明の第２の観点に係る訓練装置は、
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算部、
前記計算された筋シナジ行列W1，W2，...，Wpに含まれる単位縦ベクトルW1(1)，W2(1)，...，Wp(1)のまとまり、単位縦ベクトルW1(2)，W2(2)，...，Wp(2)のまとまり、...、単位縦ベクトルW1(n)，W2(n)，...，Wp(n)のまとまりを表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えるように構成する。
【００１０】
また、本発明の訓練装置において、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、ベクトルx,yの相関係数演算r(x,y)により、
SSI = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Wk(i),Wh(i))
と計算されるSSIである
ように構成することができる。
【００１１】
また、本発明の訓練装置において
前記SSIが所定の閾値以上であること、もしくは、前記pの増加に対して前記SSIが増加傾向にあることが、前記被験者に出力される
ように構成することができる。
【００１２】
本発明の第３の観点に係る訓練装置は、
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得部、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算部、
前記計算された制御行列C1，C2，...，Cpに含まれる横ベクトルC1(1)，C2(1)，...，Cp(1)のまとまり、横ベクトルC1(2)，C2(2)，...，Cp(2)のまとまり、横ベクトルC1(n)，C2(n)，...，Cp(n)のまとまり、を表す特徴量を計算する特徴量計算部、
前記計算された特徴量を出力する出力部
を備えるように構成する。
【００１３】
また、本発明の訓練装置において、
前記シナジ計算部では、non-negative matrix factorizationが適用され、
前記特徴量は、ベクトルx,yの相関係数演算r(x,y)により、
SSIC = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Ck(i),Ch(i))
と計算されるSSICである
ように構成することができる。
【００１４】
また、本発明の訓練装置において、
前記SSICが所定の閾値以下であること、もしくは、前記pの増加に対して前記SSICが減少傾向にあることが、前記被験者に出力される
ように構成することができる。
【００１５】
また、本発明の訓練装置において、
前記被験者の右半身の筋肉に対する右特徴量と、前記被験者の左半身の筋肉に対する左特徴量と、が計算され、
前記右特徴量および前記左特徴量の両者、もしくは、前記右特徴量および前記左特徴量の類似度が、前記被験者に出力される
ように構成することができる。
【００１６】
本発明の第４の観点に係る訓練方法は、
ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mから、n個の単位縦ベクトルW(1)，W(2)，...，W(n)を横に並べた筋シナジ行列Wと、n個の横ベクトルC(1)，C(2)，...，C(n)を縦に並べた制御行列Cと、誤差行列Eと、を、
M = WC + E
を満たすように計算するシナジ計算工程、
前記計算された筋シナジ行列Wに含まれる単位縦ベクトルW(1)，W(2)，...，W(n)のまとまりを表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えるように構成する。
【００１７】
本発明の第５の観点に係る訓練方法は、
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算工程、
前記計算された筋シナジ行列W1，W2，...，Wpに含まれる単位縦ベクトルW1(1)，W2(1)，...，Wp(1)のまとまり、単位縦ベクトルW1(2)，W2(2)，...，Wp(2)のまとまり、...、単位縦ベクトルW1(n)，W2(n)，...，Wp(n)のまとまりを表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えるように構成する。
【００１８】
本発明の第６の観点に係る訓練方法は、
振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する取得工程、
前記振舞いk = 1，2，...，pのそれぞれにおいて、前記取得された前記筋肉1，2，...，mのそれぞれの筋電位の時系列を表す横ベクトルを縦に並べた筋電位行列Mkから、n個の単位縦ベクトルWk(1)，Wk(2)，...，Wk(n)を横に並べた筋シナジ行列Wkと、n個の横ベクトルCk(1)，Ck(2)，...，Ck(n)を縦に並べた制御行列Ckと、誤差行列Ekと、を、
Mk = WkCk + Ek
を満たすように計算するシナジ計算工程、
前記計算された制御行列C1，C2，...，Cpに含まれる横ベクトルC1(1)，C2(1)，...，Cp(1)のまとまり、横ベクトルC1(2)，C2(2)，...，Cp(2)のまとまり、横ベクトルC1(n)，C2(n)，...，Cp(n)のまとまり、を表す特徴量を計算する特徴量計算工程、
前記計算された特徴量を出力する出力工程
を備えるように構成する。
【００１９】
本発明の第７の観点に係るプログラムは、コンピュータを、上記の訓練装置の各部として機能させるように構成する。
【００２０】
また、本発明のプログラムは、コンパクトディスク、フレキシブルディスク、ハードディスク、光磁気ディスク、デジタルビデオディスク、磁気テープ、半導体メモリ等のコンピュータ読取可能な情報記憶媒体に記録することができる。
【００２１】
上記プログラムは、プログラムが実行されるコンピュータとは独立して、コンピュータ通信網を介して配布・販売することができる。また、上記情報記憶媒体は、コンピュータやデジタル信号プロセッサとは独立して配布・販売することができる。
【発明の効果】
【００２２】
本発明によれば、被験者が１つの振舞いをしたり、１つ以上の振舞いからなるタスクをしたりする間の筋電位の時系列から、当該被験者の当該振舞いや当該タスクに対する熟練度を判定するのに好適な訓練装置、訓練方法、ならびに、当該訓練装置をコンピュータにより実行するためのプログラムを提供することができる。
【図面の簡単な説明】
【００２３】
【図１】訓練装置の概要構成を示す説明図である。
【図２】訓練装置が実行する訓練処理の制御の流れを示す説明図である。
【図３】ある振舞いについてシナジ数nで計算された類似度Lを示すグラフである。
【図４Ａ】被験者があるタスクを一日おきに行った際の、そのタスクに対するスコアの変化を示すグラフである。
【図４Ｂ】被験者があるタスクを一日おきに行った際の、SSI、SSICの変化を示すグラフである。
【図４Ｃ】被験者があるタスクを一日おきに行った際の、SCIの変化を示すグラフである。
【発明を実施するための形態】
【００２４】
以下に本発明の実施形態を説明する。なお、以下に説明する実施形態は説明のためのものであり、本発明の範囲を制限するものではない。したがって、当業者であればこれらの各要素もしくは全要素をこれと均等なものに置換した実施形態を採用することが可能であるが、これらの実施形態も本発明の範囲に含まれる。
【実施例１】
【００２５】
図１は、本実施形態に係る訓練装置の概要構成を示す説明図である。以下、本図を参照して説明する。
【００２６】
本実施形態に係る訓練装置１０１は、取得部１０２、シナジ計算部１０３、特徴量計算部１０４、出力部１０５を備える。訓練装置１０１は、所定のプログラムをコンピュータ上で実行することにより実現することが可能であるが、専用の電子回路により実装したり、ＦＰＧＡ（Field Programmable Gate Array）等のように、ソフトウェアからハードウェアを構成することによって実現することもである。
【００２７】
本実施形態では、ある振舞いを被験者がした間、あるいは、振舞いk = 1，2，...，pからなるタスクを被験者がした間の筋電位に基づいて、当該被験者の当該振舞いあるいは当該タスクに対する熟練度を計算して出力する。
【００２８】
たとえば、右手のリハビリを行っている被験者の場合、右手首を曲げて伸ばすことを繰り返す訓練を行うことがある。この訓練において、一回の右手首の曲げ伸ばしが、１つの振舞いに相当する。また、繰り返しをp回行うことで、1回目の曲げ伸ばし、2回目の曲げ伸ばし、...、p回目の曲げ伸ばしのそれぞれが、振舞い1、振舞い2、...、振舞いpに相当する。なお、タスクにおける各振舞いは、同じものとするのが典型的であるが、１つのタスクにおいて、異なる振舞いをすることとしても良い。
【００２９】
以下では、まず、ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列に基づく場合について説明する。
【００３０】
図２は、本実施形態に係る訓練装置が実行する訓練処理の制御の流れを示す説明図である。以下、本図を参照して説明する。
【００３１】
まず、取得部１０２は、ある振舞いを被験者がした間の筋肉1，2，...，mの筋電位の時系列を取得する（ステップＳ２０１）。本実施形態では、被験者の体のm箇所の筋電位を測定する。筋電位を測定すべき箇所は、訓練を行う部位（たとえば右手）の振舞い（たとえば右手の曲げ伸ばし）に関わる筋肉である。
【００３２】
測定は、一定の時間間隔で、１つの振舞いが開始されてから終了するまで行われる。i番目の筋肉のj番目の時点における筋電位の値は、筋電位行列Mのi行j列の要素M[i,j]に格納される。すなわち、
筋肉1の筋電位の時系列からなる横ベクトルM(1)と、
筋肉2の筋電位の時系列からなる横ベクトルM(2)と、…、
筋肉mの筋電位の時系列からなる横ベクトルM(m)と、
を縦に並べたものが、筋電位行列Mである。
【００３３】
したがって、筋電位行列Mの行数は、mである。また、筋電位行列Mの列数は測定の時間長、すなわち、振舞いの時間長と、当該振舞いの間の測定頻度、測定間隔によって変わる。
【００３４】
このようにして、筋電位行列Mが取得されると、シナジ計算部１０３は、
M = WC+E
が成立するように、筋シナジ行列W、制御行列C、誤差行列Eを計算する（ステップＳ２０２）。この際には、non-negative matrix factorizationが使用される。
【００３５】
以下、理解を容易にするため、添字kを適宜省略して説明する。
【００３６】
さて、non-negative matrix factorizationでは、誤差度を最小化、あるいは、類似度Lを最大化する。
【００３７】
ここで、筋電位行列Mの列数、制御行列Cの列数、誤差行列Eの列数がいずれもtであるとし、筋電位行列Mの行数、筋シナジ行列Wの行数、誤差行列Eの行数がいずれもmであるとし、筋シナジ行列Wの列数、制御行列Cの行数がいずれもnであるとすると、類似度Lは、以下のように定義できる。
L = 1 - 1/m × Σi=1m √〔Σj=1t E[i,j]2〕/√〔Σj=1t (WC)[i,j]2〕
【００３８】
ここで、nは、シナジ数を表す数値である。一般に、nを大きくすれば、Lも大きくなるが、適切なnの大きさも、non-negative matrix factorizationを適宜用いることで、以下のように定めることができる。
【００３９】
一般に、non-negative matrix factorizationにおいては、類似度Lが７０％以上となるようなシナジ数nを選択することが望ましいとされている。一方で、シナジ数nが大きすぎると、計算負荷が高まるほか、過適応が生じて、かえって適切な処理ができなくなる。
【００４０】
そこで、以下のような手法を用いる。
【００４１】
すなわち、n = 1，2，3，4，...のそれぞれについて、上記の類似度Lを計算する。
【００４２】
図３は、ある振舞いについてシナジ数nで計算された類似度Lを示すグラフである。以下、本図を参照して説明する。
【００４３】
本図の横軸Number of synergies nは、シナジ数nを表し、縦軸Similarity L(%)は、類似度Lを表す。本図に示すように、シナジ数nが大きくなると、類似度Lも大きくなっていくが、シナジ数nが５程度で、増加の度合が飽和しており、しかも７０％以上となっていることがわかる。したがって、飽和し始めの前後の数値、たとえば、４、５あるいは６を、以降の計算を行うためのシナジ数nとして採用することができる。
【００４４】
シナジ数nは、各被験者ごとに異なる値としても良いし、人間がある振舞いを行う際のシナジ数に大きな差はないと考えられるので、全被験者に共通する値としても良い。後者の場合には、あらかじめ実験的に何人かの被験者に振舞いをさせ、適当なnの値をnon-negative matrix factorizationにより定めてから、以降は、これで定めたnの値を他の被験者に対しても、そのまま用いる。
【００４５】
このモデルでは、被験者の中枢神経がn個の制御用信号C(1)，C(2)，...，C(n)をm個の筋肉に与えると、筋肉1は、筋電位W C(1)を満たすように動こうとし、筋肉2は、筋電位W C(2)を満たすように動こうとし、...、筋肉mは、筋電位W C(m)を満たすように動こうとする、と想定している。
【００４６】
発明者らの研究によれば、後に実験結果を示す通り、制御用信号C(1)，C(2)，...，C(n)が、互いにできるだけ独立して、非冗長なものとなっているほど、すなわち、筋シナジ行列Wにおける単位縦ベクトルW(1)，W(2)，...,W(n)が、互いにばらつかず、まとまっているほど、その振舞いに対する中枢神経の働きが良好であり、熟練度が高い、との結果が得られている。
【００４７】
そこで、特徴量計算部１０４は、制御行列Cや筋シナジ行列Wから、中枢神経の働きや熟練度と関連する特徴量を計算する（ステップＳ２０３）。
【００４８】
単位縦ベクトルW(1)，W(2)，...,W(n)が、互いにばらつかず、まとまっている、ということは単位縦ベクトルW(1)，W(2)，...,W(n)によって現わされる空間の広がりが狭いことを意味する。すなわち、各辺の長さが１で、各辺の方向がW(1)，W(2)，...,W(n)のいずれかであるn次元平行多面体のn次元体積が小さいことになる。
【００４９】
また、単位縦ベクトルW(1)，W(2)，...,W(n)同士のなす角がそれぞれ直角から遠ければ遠いほど、単位縦ベクトルW(1)，W(2)，...,W(n)は互いにばらつかず、まとまっていることになる。
【００５０】
１つの振舞いに対する筋シナジ行列Wに含まれる単位縦ベクトルW(1)，W(2)，...,W(n)がばらつかず、まとまっていることを表す特徴量としては、種々のものを採用することが可能である。
【００５１】
たとえば、ベクトルの内積を利用することが考えられる。２つのベクトルx，yに対する内積p(x,y)は、ベクトルxの各要素をx[1]，x[2]，...，x[u]と、ベクトルyの各要素をy[1]，y[2]，...，y[u]と、それぞれ表記することにより以下のように定義できる。
p(x,y) = Σi=1u(x[i]×y[i])
【００５２】
すると、筋シナジ行列Wに対する特徴量SCIを、以下のように定めることができる。
SCI = 2/〔n(n+2)〕×Σi=1n Σj=1n,j≠i p(W(i),W(j))
【００５３】
特徴量SCIは、振舞いに対する熟練度が高ければ高いほど、大きくなる数値である。
【００５４】
出力部１０５は、このようにして計算された特徴量を出力して（ステップＳ２０４）、本処理を終了する。
【００５５】
計算されたSCI等の特徴量は、リハビリ等を行っている被験者に対して出力しても良いし、その医師や監督者、トレーナーなどに対して出力することとしても良い。
【００５６】
上記のように、特徴量SCIは、振舞いに対する熟練度が高ければ高いほど、大きくなる数値である。発明者の研究によれば、特徴量SCIは、健常な人においてはある閾値以上の数値をとることがわかっている。
【００５７】
一方で、リハビリにおいては、振舞いがうまくできていることを被験者にフィードバックすることが良い、と考えられている。
【００５８】
そこで、リハビリ中の被験者がある振舞いをしたときに、その振舞いに対して計算された特徴量SCIが、上記の閾値以上であれば、その旨を被験者にフィードバックすることによって、訓練がうまく進んでいることを被験者に知得させることができる。
【００５９】
フィードバックの手法としては、たとえば、画面に特徴量SCIのグラフを表示してこれが閾値以上となればグラフの色を変化させる手法、特徴量SCIが閾値以上となれば電気刺激をリハビリ中の体の部位に与える手法、特徴量SCIが閾値以上となれば音声でその旨を知らせる手法などを採用することができる。さらに、特徴量SCIが閾値未満の場合には、リハビリに対して被験者が前向きになるような音声や画像などを出力することしても良い。
【００６０】
さて、同じ振舞いの繰り返しからなるタスクによってリハビリをしている場合、タスクに含まれる振舞いk = 1，2，...，pのそれぞれに対して、SCIkを以下のように計算しても良い。
SCIk = 2/〔n(n+2)〕×Σi=1n Σj=1n,j≠i r(Wk(i),Wk(j))
【００６１】
すると、リハビリが順調に進んでいるときは、各振舞いに対するSCI1，SCI2，...,SCIpは増加の傾向にあると考えられ、ある程度の成果が得られた後は、飽和するものと考えられる。そこで、各振舞いについて順次SCIを計算して出力することにより、当該振舞いにおける成果を被験者に確認させてから、次の振舞いを開始させることによって、リハビリの効果を高めることができる。
【００６２】
なお、上記の例では、右手首の１回の曲げ伸ばしを１つの振舞いと考えて説明したが、右手首の曲げ伸ばしの繰り返しを１つの振舞いと考えても良い。この場合には、複数の振舞いからなるタスク全体や複数のタスクからなる訓練メニュー全体を、まとめて１つの振舞いと解釈することで、SCIを計算することが可能である。たとえば、「ある日におこなったリハビリ全体のSCI」などを計算することとしても良い。
【００６３】
この場合、一人の被験者について日毎にSCIの変化を蓄積し、SCIが増加の傾向にあるようであれば、現在被験者が行っているタスクや訓練メニューは適切であると考えることができる。SCIの増加が飽和した段階で、このタスクや訓練メニューは終了し、新たなタスクや訓練メニューを開始することで、リハビリの進行を客観的に調整することができる。
【００６４】
さて、上記実施形態では、１回の振舞いに対する特徴量としてSCIを計算したが、以下の実施形態では、振舞いk = 1，2，...，pからなる１つのタスクに対する特徴量であるSSIもしくはSSICを計算する。典型的には、１つのタスクに含まれる振舞いk = 1，2，...，pは、同じ動作とすべきであるが、異なる動作としても良い。
【００６５】
取得部１０２は、上記と同様に、p個の振舞いk = 1，2，...，pからなるタスクに対して、筋肉1，2，...，mの振舞い1を被験者がした間の筋電位の時系列からなる筋電位行列M1，振舞い2を被験者がした間の筋電位の時系列からなる筋電位行列M2，...,振舞いpを被験者がした間の筋電位の時系列からなる筋電位行列Mpを取得する（ステップＳ２０１）。
【００６６】
ここで取得される筋電位行列M1，M2，...,Mpの行数はいずれもmであるが、筋電位行列M1，M2，...,Mpの列数は、それぞれの振舞いをする時間長によって異なっても良い。
【００６７】
ついで、シナジ計算部１０３は、上記と同様に、
Mk = WkCk+Ek    (k = 1，2，...，p)
が成立するように、筋シナジ行列Wk、制御行列Ck、誤差行列Ekを計算する（ステップＳ２０２）。
【００６８】
さて、振舞いを繰り返すことでその振舞いに対する熟練度が増すと、筋シナジ行列W1，W2，...,Wpは、次第に変化しなくなっていくと考えられる。
【００６９】
このため、熟練度が増すと、筋シナジ行列W1，W2，...,Wpの同じ位置に含まれる単位縦ベクトルは、次第にまとまり、ばらつきがなくなっていく。
【００７０】
SSIは、筋シナジ行列W1，W2，...,Wpの同じ位置に含まれる単位縦ベクトルのまとまりの度合を示す特徴量である。
【００７１】
一方で、振舞いを繰り返すとしても、被験者の状態や環境は刻一刻と変化する。したがって、各シナジに与える制御信号は、この変化に呼応するため、次第にまとまるのではなく、大きくばらつく可能性がある。
【００７２】
このため、筋シナジ行列W1，W2，...,Wpの同じ位置に含まれる縦ベクトルが次第にまとまるのとは裏腹に、制御行列C1，C2，...,Cpの同じ位置に含まれる横ベクトルはばらつく。
【００７３】
SSICは、制御行列C1，C2，...,Cpの同じ位置に含まれる横ベクトルのまとまりの度合を示す特徴量である。
【００７４】
以下では、ベクトルxに含まれる要素の平均e(x)、ベクトルxに含まれる要素の分散v(x)、ベクトルxに含まれる要素の標準偏差s(x)、２つのベクトルx,yに対する相関係数演算r(x,y)は、ベクトルxの各要素をx[1]，x[2]，...，x[u]と、ベクトルyの各要素をy[1]，y[2]，...，y[u]と、それぞれ表記すると、以下のように定義できる。
e(x) = (1/u)×Σi=1u x[i],
v(x) = (1/u)×Σi=1u (x[i]-e(x))2,
s(x) = v(x)1/2,
r(x,y) = Σi=1u(x[i]-e(x))×(y[i]-e(y))/〔m×s(x)×s(y)〕
【００７５】
特徴量計算部１０３は、特徴量SSIあるいはSSICを、以下のように計算する。
SSI = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Wk(i),Wh(i)),
SSIC = 2/〔n×p(p-1)〕×Σi=1n Σk=1p Σh=1p,h≠k r(Ck(i),Ch(i))
【００７６】
上記のように、SSIは、大きければ大きいほど
単位縦ベクトルW1(1)，W2(1)，...，Wp(1)のまとまり、
単位縦ベクトルW1(2)，W2(2)，...，Wp(2)のまとまり、...、
単位縦ベクトルW1(n)，W2(n)，...，Wp(n)のまとまり
が大きく、熟練度が高いものと考えられる。
【００７７】
一方、SSICは、小さければ小さいほど
横ベクトルC1(1)，C2(1)，...，Cp(1)のばらつき、
横ベクトルC1(2)，C2(2)，...，Cp(2)のばらつき、
横ベクトルC1(n)，C2(n)，...，Cp(n)のばらつき
が大きく、すなわち、それぞれのまとまりが小さく、熟練度が高いものと考えられる。
【００７８】
また、タスクが被験者にとって有用なものであれば、pが増加していくにつれ、SSIは増加の傾向になりSSICは減少の傾向になる、と考えられる。
【００７９】
そこで、出力部１０４は、上記と同様に計算された特徴量を出力して（ステップＳ２０４）、本処理を終了する。
【００８０】
本実施形態では、p回目の振舞いが終了した段階で、振舞い1，2，...，pに対するSSIやSSICを計算して、これらをグラフ等で出力しても良いし、SSIが、ある閾値以上であるか増加傾向であれば、および／またはSCIkが、ある閾値以下であるか減少傾向であれば、その旨を電気刺激や音声などで被験者に知得させることとしても良い。
【００８１】
また、これと合わせて、SCIkをグラフ等で出力しても良いし、SCIkがある閾値以上であるか増加傾向であれば、その旨を電気刺激や音声などで被験者に知得させることとしても良い。
【００８２】
図４Ａ、４Ｂ、４Ｃは、被験者があるタスクを一日おきに行った際の、そのタスクに対するスコア、SSI、SSIC、SCIの変化を示すグラフである。以下、これらの図を参照して説明する。
【００８３】
図４Ａは、縦軸がスコア、横軸が経過日数を表すグラフである。スコアは、被験者に与えられたタスクに対してあらかじめ定められた動作を行った場合には加点がされ、動作に失敗した場合には減点がされるように定められており、リハビリが進むにつれて、スコアが向上の傾向にあることがわかる。
【００８４】
図４Ｂは、縦軸がSSIならびにSSIC、横軸が経過日数を表すグラフである。リハビリが進むにしたがって、SSI（W）が増加の傾向にあり、SSIC（C）が減少の傾向にあることがわかる。
【００８５】
図４Ｃは、縦軸がSCI、横軸が経過日数を表すグラフである。リハビリが進むにしたがって、SCIが増加の傾向にあることがわかる。
【００８６】
このように、リハビリ等で被験者にあるタスクをさせるときを考えると、最初は、被験者は熟練度が低い。
【００８７】
１つ以上の振舞いからなるタスクを被験者が繰り返して訓練すると、SSIの増加、SSICの減少が見られる。また、タスク全体のSCIの増加が見られるほか、タスクに含まれる振舞いの多くについて、SCI1，SCI2，SCI3，...は次第に増加する。
【００８８】
すなわち、これらの現象が見られる間は、リハビリ用の運動として当該タスクが有効に機能していることを意味する。
【００８９】
この後、被験者が十分にタスクに熟練すると、SSI、SSIC、SCI、SCIkがあまり変化しなくなる。
【００９０】
この段階に至れば、他のタスクを課したり、タスクにおける振舞いの回数を変化させたり、などによって、被験者のリハビリをさらに進めることができる。また、リハビリを終了しても良い、と判定することも可能である。
【００９１】
このように、本実施形態によれば、リハビリにおいて被験者に与えられるタスクや振舞いに対して被験者がどの程度熟練したかを客観的に把握することができるようになる。
【実施例２】
【００９２】
上記実施例では、ある振舞いやあるタスクを行った際に計算された特徴量を、あらかじめ求められた閾値と比較した結果や、その特徴量の増減の傾向によって、電気刺激や音声の出力をすることとしていたが、本実施形態は、人間の体が左右対称であることを利用して、出力を行う。
【００９３】
たとえば、被験者の右半身と左半身の一方に障害があってリハビリを行う場合を考える。
【００９４】
本実施形態では、被験者には、左半身と右半身で同時に同じ振舞いあるいはタスクを試行させる。そして、左右対称に設けた測定点から、左半身の筋電位に対する左特徴量と、右半身の筋電位に対する右特徴量と、を、上記実施例の特徴量の計算手法により計算する。
【００９５】
たとえば、右手に障害があるときには、被験者は、右手首と左手首を同時に曲げ伸ばしさせる。そして、右手の筋肉と、左手の筋肉と、に、左右対称に測定点を設け、それぞれの筋電位を測定して、右手に対する特徴量と、左手に対する特徴量と、を、計算するのである。
【００９６】
両半身が健常である場合には、右特徴量と、左特徴量と、は、ほぼ等しい値をとるものと思われる。
【００９７】
そこで、本実施形態では、右特徴量と、左特徴量と、を、同時にグラフ表示したり、右特徴量と左特徴量の類似度が一定以上（右特徴量と左特徴量の差が一定以下）であれば、振舞いやタスクに成功したものとして、電気刺激や画像、音声で被験者にフィードバックをする。右特徴量と左特徴量がほぼ一致するようになれば、リハビリが十分にできたことになる。
【００９８】
本実施形態によれば、左右対称な運動をリハビリとして行うことで、あらかじめ閾値を定めることなく、客観的にリハビリの進行度を判断することができるようになる。
【産業上の利用可能性】
【００９９】
本発明によれば、被験者が１つの振舞いをしたり、複数の振舞いからなるタスクをしたりする間の筋電位の時系列から、当該被験者の当該振舞いや当該タスクに対する熟練度を判定するのに好適な訓練装置、訓練方法、ならびに、当該訓練装置をコンピュータにより実行するためのプログラムを提供することができる。
【符号の説明】
【０１００】
１０１  訓練装置
１０２  取得部
１０３  シナジ計算部
１０４  特徴量計算部
１０５  出力部
【図１】
0006011807.tif 000002
【図２】
0006011807.tif 000003
【図３】
0006011807.tif 000004
【図４Ａ】
0006011807.tif 000005
【図４Ｂ】
0006011807.tif 000006
【図４Ｃ】
0006011807.tif 000007