(57)【特許請求の範囲】
【請求項１】
物理量のｍ階時間積分値（ｍは１以上の整数）の真値が既知である複数の第１期間と運動解析の対象である第２期間とを含む期間のセンサー部の出力データを取得するデータ取得部と、
前記出力データの前記複数の第１期間のｍ階時間積分値と前記真値との差を用いて、前記複数の第１期間と前記第２期間とを含む前記期間における、前記センサー部が検出する前記物理量の値の真値に対する誤差の時間関数を推定する誤差時間関数推定部と、
前記出力データをｍ階時間積分し、前記誤差時間関数推定部が推定した前記誤差の時間関数を用いて、前記出力データの前記第２期間のｍ階時間積分値を補正するデータ補正部と、
前記データ補正部が補正した前記第２期間のｍ階時間積分値を用いて、前記対象物の運動解析情報を生成する運動解析情報生成部と、を含み、
前記複数の第１期間の少なくとも１つは、前記第２期間の開始前の期間であり、
前記複数の第１期間の少なくとも１つは、前記第２期間の終了後の期間であり、
前記複数の第１期間は、前記対象物が静止している期間であり、
前記第２期間は、前記対象物がスイング運動をしている期間である、運動解析装置。
【請求項２】
請求項１において、
前記運動解析は、テニスの運動解析である、運動解析装置。
【請求項３】
請求項１又は２において、
前記対象物は、テニスラケットである、運動解析装置。
【請求項４】
請求項１乃至３のいずれか一項において、
前記誤差時間関数推定部は、
前記誤差の時間関数を多項式で近似して当該多項式の係数を算出することにより、前記誤差の時間関数を推定する、運動解析装置。
【請求項５】
請求項１乃至３のいずれか一項において、
前記誤差時間関数推定部は、
前記誤差の時間関数ε（ｔ）を式（８）で表されるｎ次多項式ｇ（ｔ）で近似し、
【数８】
0005948011.tif 000019
前記出力データのｍ階時間積分値Ｘｍ（ｔ）と前記物理量のｍ階時間積分値の真値Ｆｍ（ｔ）との差Ｅｍ（ｔ）が、積分定数ｃｋを考慮して前記ｎ次多項式をｍ階時間積分した式（１１）で表される多項式Ｇｍ（ｔ）と近似するものとして、
【数１１】
0005948011.tif 000020
前記第１期間における複数の時刻ｔｒｉにおける、前記出力データのｍ階時間積分値Ｘｍ（ｔｒｉ）と前記物理量のｍ階時間積分値の真値Ｆｍ（ｔｒｉ）を用いて、式（１３）で表される前記ｎ次多項式ｇ（ｔ）の係数ａｋ及び前記積分定数ｃｋについての過決定連立方程式を立て、
【数１３】
0005948011.tif 000021
前記過決定連立方程式（１３）から前記係数ａｋ及び前記積分定数ｃｋを求めて式（１１）に代入することにより、多項式Ｇｍ（ｔ）を求め、
前記出力データの前記第２期間のｍ階時間積分値Ｘｍ（ｔ）から前記多項式Ｇｍ（ｔ）を減算することにより、前記出力データの前記第２期間のｍ階時間積分値を補正する、運動解析装置。
【請求項６】
請求項１乃至５のいずれか一項において、
前記センサー部は、
前記物理量として加速度および角速度の少なくとも一方を検出する、運動解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、運動解析装置に関する。
【背景技術】
【０００２】
様々な分野において人や物体の運動を解析する装置が必要とされている。例えば、テニスラケットやゴルフクラブのスイング軌道、野球のピッチングやバッティングのフォーム等を解析し、解析結果から改善点を明らかにすることで競技力の向上につなげることができる。
【０００３】
現在、実用的な運動解析装置としては、マークがつけられた被測定物を赤外線カメラ等で連続撮影し、撮影された連続画像を用いてマークの移動軌跡を算出することで、運動を解析するものが一般的である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−２４４８８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、このような運動解析装置では、画像を撮影するための赤外線カメラが必要であるため装置が大がかりなものになってしまい、取り扱いにくいという問題がある。例えば、テニスの練習中の画像を複数の角度から撮影したい場合、撮影したい角度に合わせて赤外線カメラの位置を移動させるかプレイヤーの向きを変える必要がある。
【０００６】
これに対して、近年、被測定物に小型の慣性センサーを取り付け、センサーの出力データから被測定物の運動を解析する装置が提案がされており、赤外線カメラが不要であるため取り扱いが容易であるという利点がある。例えば、加速度センサーが検出する加速度値ａ（ｔ）に対して、それぞれ次式（１）及び（２）に示すような時間積分処理を行うことで、被測定物の速度ｖ（ｔ）及び位置ｐ（ｔ）を算出することができる。
【０００７】
【数１】
0005948011.tif 000002
【０００８】
【数２】
0005948011.tif 000003
【０００９】
ところが、一般に慣性センサーの出力値には観測しようとする値の他に誤差が含まれている。従って、例えば、加速度センサーの出力データｘ（ｔ）は、加速度値ａ（ｔ）と誤差ε（ｔ）を用いて次式（３）のように表すことができる。
【００１０】
【数３】
0005948011.tif 000004
【００１１】
そのため、加速度センサーの出力データｘ（ｔ）から、それぞれ次式（４）及び（５）に示すような時間積分処理を行って被測定物の速度ｖ（ｔ）及び位置ｐ（ｔ）を算出すると、誤差ε（ｔ）も時間積分されるため、時間ｔの経過に伴って速度ｖ（ｔ）及び位置ｐ（ｔ）の誤差が急速に増大する。
【００１２】
【数４】
0005948011.tif 000005
【００１３】
【数５】
0005948011.tif 000006
【００１４】
要するに、慣性センサーを用いた運動解析装置では、実際にはセンサー自体の特性が十分ではなく、慣性センサーの出力データを積分処理して姿勢、速度、位置等を算出する場合、センサーの出力に含まれる誤差が積分処理によって急速に増大するため、十分な解析（測定）性能が得られていないのが現状である。
【００１５】
本発明は、以上のような問題点に鑑みてなされたものであり、本発明のいくつかの態様によれば、取り扱いが容易であるとともに十分な精度の解析情報が得られる運動解析装置を提供することができる。
【課題を解決するための手段】
【００１６】
（１）本発明は、対象物に取り付けられ、物理量を検出するセンサー部と、前記物理量のｍ階時間積分値（ｍは１以上の整数）の真値が既知である第１期間と運動解析の対象である第２期間とを含む期間の前記センサー部の出力データを取得するデータ取得部と、前記出力データをｍ階時間積分し、前記出力データの前記第１期間のｍ階時間積分値と前記真値との差に基づいて、前記センサー部が検出する前記物理量の値の前記真値に対する誤差の時間関数を推定する誤差時間関数推定部と、前記誤差時間関数推定部の推定結果に基づいて、前記出力データの前記第２期間のｍ階時間積分値を補正するデータ補正部と、前記データ補正部が補正した前記第２期間のｍ階時間積分値に基づいて、前記対象物の運動解析情報を生成する運動解析情報生成部と、を含む、運動解析装置である。
【００１７】
解析の対象となる対象物は、人であってもよいし、人以外の物（例えば運動器具や車両等）であってもよい。
【００１８】
対象物の運動を解析するための情報は、例えば、対象物の軌跡の情報や対象物の速度変化の情報等であってもよい。
【００１９】
ｍ階時間積分とは、連続時間系でのｍ階時間積分であってもよいし、離散時間系でのｍ階時間積分（ｍ階時間差分）であってもよい。
【００２０】
本発明によれば、センサー部の検出誤差を時間関数として推定し、推定した誤差の時間関数を用いて、検出対象の物理量のｍ階時間積分値を補正することにより、十分な精度の解析情報を生成することができる。また、本発明によれば、赤外線カメラの代わりにセンサーを用いるので簡便な構成にすることができ、取り扱いも容易である。
【００２１】
（２）この運動解析装置において、前記誤差時間関数推定部は、前記誤差の時間関数を多項式で近似して当該多項式の係数を算出することにより、前記誤差の時間関数を推定するようにしてもよい。
【００２２】
このようにすれば、比較的簡単な演算により、検出誤差の時間関数を十分な精度で推定することができる。なお、多項式の次数は、運動解析に要求される精度に応じて決定すればよい。
【００２３】
また、例えば、前記誤差時間関数推定部は、前記データ取得部が取得したデータの前記第１期間におけるｍ階時間積分値の真値に対する誤差が前記多項式の前記第１期間におけるｍ階時間積分値と近似するものとする過決定連立方程式を解いて、前記多項式の係数を算出するようにしてもよい。
【００２４】
このように、第１期間でより多くのデータを取得して過決定連立方程式を立てることで、検出誤差の時間関数の推定精度を高めることができる。なお、過決定連立方程式を、例えば、最小二乗法で解いてもよい。
【００２５】
（３）この運動解析装置において、前記第１期間は複数設定され、前記誤差時間関数推定部は、前記データ取得部が複数の前記第１期間の各々において取得したデータに基づいて、前記誤差の時間関数を推定するようにしてもよい。
【００２６】
このように複数の第１期間を設けることで、検出誤差の時間関数の推定精度をより高めることができる。
【００２７】
（４）この運動解析装置において、複数の前記第１期間の少なくとも１つは、前記第２期間の開始前の期間であり、複数の前記第１期間の少なくとも１つは、前記第２期間の終了後の期間であるようにしてもよい。
【００２８】
このようにすれば、運動解析の対象である第２期間における検出誤差の時間関数の推定精度をより高めることができるので、より精度の高い運転解析情報を生成することができる。
【００２９】
（５）この運動解析装置において、前記第１期間は、前記対象物が静止している期間であるようにしてもよい。
【００３０】
このようにすれば、例えば、第１期間において対象物の速度、姿勢、位置を既知にすることができる。
【００３１】
（６）この運動解析装置において、前記センサー部は、前記物理量として加速度および角速度の少なくとも一方を検出するようにしてもよい。
【図面の簡単な説明】
【００３２】
【図１】本実施形態の運動解析装置の構成を示す図。
【図２】処理部による運動解析情報の生成処理の一例を示すフローチャート図。
【図３】データ取得期間、第１期間及び第２期間の例を示す図。
【図４】誤差時間関数の推定処理及びデータ補正処理のフローチャート図。
【図５】本実験例におけるセンサー部の概略構成図。
【図６】本実験例におけるセンサー部の取り付け例を示す図。
【図７】本実験例における被験者の動作シーケンスについて説明するための図。
【図８】本実験例における座標系の定義について説明するための図。
【図９】本実験例における処理部による処理のフローチャート図。
【図１０】本実験例における軌跡データを示す図。
【図１１】本実施形態の手法と従来手法による２つの軌跡データを比較する図。
【発明を実施するための形態】
【００３３】
以下、本発明の好適な実施形態について図面を用いて詳細に説明する。なお、以下に説明する実施の形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。
【００３４】
図１は、本実施形態の運動解析装置の構成を示す図である。
【００３５】
本実施形態の運動解析装置１は、１又は複数のセンサー部１０とホスト端末２０を含んで構成され、対象物の運動を解析する。センサー部１０とホスト端末２０は無線接続されていてもよいし、有線接続されていてもよい。
【００３６】
センサー部１０は、運動解析の対象物に取り付けられ、所与の物理量を検出する処理を行う。本実施形態では、センサー部１０は、１又は複数のセンサー１００、データ処理部１１０、通信部１２０を含んで構成されている。
【００３７】
センサー１００は、所与の物理量を検出し、検出した物理量（例えば、加速度、角速度、速度、角加速度など）の大きさに応じた信号（データ）を出力するセンサーである。センサー１００は、例えば慣性センサーであってもよい。
【００３８】
データ処理部１１０は、各センサー１００の出力データの同期を取り、当該データを時刻情報などと組合せたパケットにして通信部１２０に出力する処理を行う。さらに、データ処理部１１０は、センサー１００のバイアス補正や温度補正の処理を行うようにしてもよい。なお、バイアス補正や温度補正の機能をセンサー１００に組み込んでもよい。
【００３９】
通信部１２０は、データ処理部１１０から受け取ったパケットデータをホスト端末２０に送信する処理を行う。
【００４０】
ホスト端末２０は、処理部（ＣＰＵ）２００、通信部２１０、操作部２２０、ＲＯＭ２３０、ＲＡＭ２４０、不揮発性メモリー２５０、表示部２６０を含んで構成されている。
【００４１】
通信部２１０は、センサー部１０から送信されたデータを受信し、処理部２００に送る処理を行う。
【００４２】
操作部２２０は、ユーザーからの操作データを取得し、処理部２００に送る処理を行う。操作部２２０は、例えば、タッチパネル型ディスプレイ、ボタン、キー、マイクなどである。
【００４３】
ＲＯＭ２３０は、処理部２００が各種の計算処理や制御処理を行うためのプログラムや、アプリケーション機能を実現するための各種プログラムやデータ等を記憶している。
【００４４】
ＲＡＭ２４０は、処理部２００の作業領域として用いられ、ＲＯＭ２３０から読み出されたプログラムやデータ、操作部２２０から入力されたデータ、処理部２００が各種プログラムに従って実行した演算結果等を一時的に記憶する記憶部である。
【００４５】
不揮発性メモリー２５０は、処理部２００の処理により生成されたデータのうち、長期的な保存が必要なデータを記録する記録部である。
【００４６】
表示部２６０は、処理部２００の処理結果を文字やグラフ、その他の画像として表示するものである。表示部２６０は、例えば、ＣＲＴ、ＬＣＤ、タッチパネル型ディスプレイ、ＨＭＤ（ヘッドマウントディスプレイ）などである。なお、１つのタッチパネル型ディスプレイで操作部２２０と表示部２６０の機能を実現するようにしてもよい。
【００４７】
処理部２００は、ＲＯＭ２４０に記憶されているプログラムに従って、センサー部１０から通信部２１０を介して受信したデータに対する各種の計算処理や、各種の制御処理（表示部２６０に対する表示制御等）を行う。
【００４８】
特に、本実施形態では、処理部２００は、以下に説明するデータ取得部２０２、誤差時間関数推定部２０４、データ補正部２０６、運動解析情報生成部２０８として機能する。
【００４９】
データ取得部２０２は、センサー１００の検出対象である物理量のｍ階時間積分値の真値が既知である第１期間と運動解析の対象である第２期間とを含む期間において、センサー部１０の出力データを取得する処理を行う。取得したデータは例えばＲＡＭ２４０に記憶される。
【００５０】
誤差時間関数推定部２０４は、センサー部１０の出力データをｍ階時間積分し、当該出力データの第１期間のｍ階時間積分値と真値との差に基づいて、センサー部１０が検出する物理量の値の真値に対する誤差の時間関数（以下、「誤差時間関数」という）を推定する処理を行う。
【００５１】
データ補正部２０６は、誤差時間関数推定部２０４の推定結果に基づいて、センサー部１０の出力データの第２期間のｍ階時間積分値を補正する処理を行う。
【００５２】
運動解析情報生成部２０８は、データ補正部２０６が補正した第２期間のｍ階時間積分値に基づいて、対象物の運動を解析するための情報（以下、「運動解析情報」という）を生成する処理を行う。生成した運動解析情報は、文字、グラフ、図などで表示部２６０に表示させてもよいし、ホスト端末２０の外部に出力してもよい。
【００５３】
図２は、処理部２００による運動解析情報の生成処理の一例を示すフローチャート図である。
【００５４】
まず、処理部２００は、データ取得部２０２として、データ取得期間が終了するまで（ステップＳ２０のＮ）センサー部１０から新たなデータを周期的に取得する（ステップＳ１０）。
【００５５】
次に、データ取得期間が終了すると（ステップＳ２０のＹ）、処理部２００は、データをｍ階積分し（ステップＳ２１）、誤差時間関数推定部２０４として、ステップＳ１０で取得したデータの第１期間におけるｍ階時間積分値と真値との差に基づいて、誤差時間関数を推定する（ステップＳ３０）。
【００５６】
次に、処理部２００は、データ補正部２０６として、ステップＳ３０で推定した時間関数に基づいて、ステップＳ１０で取得したデータの第２期間におけるｍ階時間積分値を補正する（ステップＳ４０）。
【００５７】
最後に、処理部２００は、運動解析情報生成部２０８として、ステップＳ４０で補正した第２期間におけるｍ階時間積分値に基づいて、運動解析情報を生成する（ステップＳ５０）。
【００５８】
図３（Ａ）及び図３（Ｂ）は、データ取得期間、第１期間及び第２期間の例を示す図である。
【００５９】
図３（Ａ）の例では、時刻ｔ２〜ｔ３に解析対象物が運動を行う第２期間が設けられ、第２期間の前後に、時刻ｔ０〜ｔ１とｔ４〜ｔ５の時間的に離れた２つの第１期間が設けられている。そして、時刻ｔ０〜ｔ５がデータ取得期間になっており、このデータ取得期間では、ホスト端末２０により、センサー部１０の出力データが一定の周期でサンプリング（取得）される。２つの第１期間では、センサー部１０の検出対象である物理量のｍ階積分値の真値が既知であるので、センサー部１０の出力データのｍ階時間積分値と真値との差がわかる。この差の情報から、センサー部１０の出力データに対する誤差時間関数をデータ取得期間全体にわたって推定することができる。なお、１番目の第１期間（時刻ｔ０〜ｔ１）と２番目の第１期間（時刻ｔ４〜ｔ５）のいずれか一方はなくてもよいが、誤差時間関数の推定精度を高めるためには、第２期間の前後に第１期間を設けた方がよい。要するに、誤差時間関数の推定精度を高めるためには、電源変動や温度変動などに起因するランダムな誤差の変動を反映させて誤差時間関数を推定することが有効であるため、時間的に離れた複数の第１期間を設けることが望ましい。特に、第２期間の前後に第１期間を設けることで、第２期間における推定誤差の精度が高まるので、第２期間におけるデータの補正精度を向上させることができる。
【００６０】
また、図３（Ｂ）の例では、時刻ｔ２〜ｔ３と時刻ｔ４〜ｔ５に、解析対象物が運動を行う２つの第２期間が設けられ、１番目の第２期間（時刻ｔ２〜ｔ３）の前に１番目の第１期間（時刻ｔ３〜ｔ４）、２つの第２期間の間に２番目の第１期間（時刻ｔ３〜ｔ４）、２番目の第２期間の後に３番目の第１期間（時刻ｔ６〜ｔ７）が設けられている。そして、時刻ｔ０〜ｔ７がデータ取得期間になっている。３つの第１期間では、センサー部１０の検出対象である物理量のｍ階積分値の真値が既知であるので、センサー部１０の出力データのｍ階時間積分値と真値との差がわかる。この差の情報から、センサー部１０の出力データに対する誤差時間関数をデータ取得期間全体にわたって推定することができる。なお、図３（Ｂ）の例では、運動解析の対象となる第２期間が２つ設けられているので、２つの第２期間を挟むように、時間的に離れた位置に３つの第１期間を設けることで、２つの第２期間における誤差時間関数の推定精度を高めることができる。すなわち、運動解析の対象となるそれぞれの第２期間の前後に第１期間を設けることで、解析対象の運動を時間をおいて繰り返し行うような場合でも、それぞれの第２期間におけるデータの補正精度を向上させることができる。
【００６１】
［誤差時間関数の推定及びデータ補正］
次に、誤差時間関数の推定及びデータ補正の手法の一例について詳しく説明する。
【００６２】
まず、処理部２００での算出対象となる物理量の時刻ｔにおける値をＦｍ（ｔ）とし、センサー部１０がそのｍ次導関数の値ｆ（ｔ）を測定する場合、次式（６）が成立する。
【００６３】
【数６】
0005948011.tif 000007
【００６４】
ここで、センサー部１０の出力データｘ（ｔ）は誤差ε（ｔ）を含むものとすると、ｘ（ｔ）は次式（７）のように表される。
【００６５】
【数７】
0005948011.tif 000008
【００６６】
次式（８）のように、誤差時間関数ε（ｔ）をｎ次の多項式ｇ（ｔ）で近似することを考える。
【００６７】
【数８】
0005948011.tif 000009
【００６８】
センサー部１０の出力データｘ（ｔ）をｍ階時間積分した結果であるＸｍ（ｔ）にも算出対象となる物理量Ｆｍ（ｔ）の他にε（ｔ）と初期状態誤差（積分定数）による誤差成分Ｅｍ（ｔ）が含まれる。従って、Ｘｍ（ｔ）は次式（９）のように表される。
【００６９】
【数９】
0005948011.tif 000010
【００７０】
この誤差成分Ｅｍ（ｔ）もｇ（ｔ）のｍ階時間積分に対して積分定数（初期状態誤差）ｃｋを考慮した多項式Ｇｍ（ｔ）で近似できると考えると、次式（１０）、（１１）が成立する。
【００７１】
【数１０】
0005948011.tif 000011
【００７２】
【数１１】
0005948011.tif 000012
【００７３】
従って、ある時刻ｔｒにおける物理量Ｆｍ（ｔｒ）がわかっている場合、次式（１２）の関係が成立する。
【００７４】
【数１２】
0005948011.tif 000013
【００７５】
この関係式（１２）を算出対象となる物理量の値がわかっている時刻の数だけ準備することで、近似多項式（１１）の係数ａｋおよびｃｋについて以下のような過決定連立方程式（１３）を立てることができる。
【００７６】
【数１３】
0005948011.tif 000014
【００７７】
この過決定連立方程式（１３）から例えば最小二乗法によって近似多項式（１１）の係数ａｋおよびｃｋを求める。
【００７８】
【数１４】
0005948011.tif 000015
【００７９】
【数１５】
0005948011.tif 000016
【００８０】
この係数ａｋおよびｃｋを使って近似多項式ｇ（ｔ）およびＧｍ（ｔ）が決定するので、物理量Ｆｍ（ｔ）とそのｍ次導関数の値ｆ（ｔ）をそれぞれ次式（１６）、（１７）のように推定できる。
【００８１】
【数１６】
0005948011.tif 000017
【００８２】
【数１７】
0005948011.tif 000018
【００８３】
以上に説明した手法に基づく誤差時間関数の推定処理及びデータ補正処理のフローチャートを図４に示す。
【００８４】
まず、取得したデータｘ（ｔ）をｍ階時間積分してＸｍ（ｔ）を計算する（ステップＳ３２）。
【００８５】
次に、誤差時間関数ε（ｔ）を多項式ｇ（ｔ）で近似し、第１期間の各時刻ｔｒにおけるｍ階時間積分値Ｘｍ（ｔｒ）と真値Ｆｍ（ｔｒ）を用いて過決定連立方程式（１３）を生成する（ステップＳ３４）。
【００８６】
次に、ステップＳ３４で生成した過決定連立方程式（１３）を解いてｇ（ｔ）の係数値ａｋ，ｃｋを算出する（ステップＳ３６）。
【００８７】
次に、ステップＳ３６で算出した係数値ａｋ，ｃｋを用いて、式（１１）よりＧｍ（ｔ）を計算する（ステップＳ３８）。
【００８８】
最後に、ステップＳ３２で計算したＸｍ（ｔ）とステップＳ３６で計算したＧｍ（ｔ）を用いて、式（１６）より、Ｆｍ（ｔ）を計算する（ステップＳ４２）。
【００８９】
なお、ステップＳ３２〜Ｓ３８の処理は、図２のフローチャートのステップＳ３０の処理に対応し、ステップＳ４２の処理は、図２のフローチャートのステップＳ４０の処理に対応する。
【００９０】
以上に説明したように、本実施形態の運動解析装置によれば、センサー部１０の出力データの誤差時間関数を推定してセンサー部１０の出力データのｍ階時間積分値を補正することにより、十分な精度の運動解析情報を生成することができる。また、本実施形によれば、赤外線カメラの代わりにセンサーを用いるので簡便な構成かつ取り扱い易い運動解析装置を実現することができる。
【００９１】
また、本実施形態によれば、誤差時間関数を多項式で近似することで、例えば式（１５）のような比較的簡単な演算により、誤差時間関数を十分な精度で推定することができる。さらに、第１期間でより多くのデータを取得して過決定連立方程式（１３）を立てることで、誤差時間関数の推定精度を高めることができる。
【００９２】
［運動解析の実験例］
次に、本実施形態の運動解析手法を適用した実験例について説明する。本実験例では、図５に示すように構成されたセンサー部１０を、図６に示すように、解析対象物であるテニスラケットのグリップエンドに取り付け、被験者がこのテニスラケットを持ってテニスボールを打った時のテニスラケットのトップ３０２とグリップエンド３０４の軌跡（運動解析情報の一例）を表示する。
【００９３】
図５に示すように、本実験例で使用したセンサー部１０は、図１に示したセンサー１００として、ｘ軸、ｙ軸、ｚ軸方向の加速度を検出する３軸加速度センサー１０２ｘ，１０２ｙ，１０２ｚ（慣性センサーの一例）と、ｘ軸、ｙ軸、ｚ軸方向の角速度を検出する３軸ジャイロセンサー（角速度センサー）１０４ｘ，１０４ｙ，１０４ｚ（慣性センサーの一例）からなる６軸モーションセンサーを備えている。ｘ軸、ｙ軸、ｚ軸は、例えば右手系に従って決められる。
【００９４】
データ処理部１１０は、６軸モーションセンサーの出力データの同期を取って通信部１２０に出力する処理の他、６軸モーションセンサーの取付け角のずれによる検出誤差を補正する処理なども行う。
【００９５】
通信部１２０は、データ処理部１１０から受け取ったデータをホスト端末２０に送信する処理を行う。
【００９６】
このセンサー部１０を、例えば、図６に示すようにテニスラケット３００のグリップエンド３０４にｘ軸がフェイス（打球面）と垂直になるように取り付ける。センサー部１０を取り付ける向きは任意の向きでよい。例えば、図６のように、ｘ軸方向が紙面奥から手前に向かう垂線方向、ｙ軸方向が水平右方向、ｚ軸方向が垂直上方向になるように、センサー部１０を取り付ける。
【００９７】
この実験例では、被験者に一定の動作シーケンスを行わせる。この動作シーケンスについて図７を用いて説明する。まず、テニスラケット３００をあらかじめ決められた第１の位置に置いて少なくとも１秒程度静止させる（時刻ｔ０〜ｔ１）。次に、被験者がテニスラケット３００を持って第２の位置に移動し、スイングの準備をする（時刻ｔ１〜ｔ２）。次に、被験者に向けてテニスボールを送り、被験者がテニスラケット３００でテニスボールを打つ（時刻ｔ２〜ｔ３）。次に、スイング終了後、被験者がテニスラケット３００を持って第１の位置に移動し、テニスラケットを第１の位置に置く（時刻ｔ３〜ｔ４）。最後に、テニスラケット３００を少なくとも１秒程度静止させる（時刻ｔ４〜ｔ５）。時刻ｔ０〜ｔ５の期間がデータ取得期間に相当し、例えば毎秒５００サンプルのサンプリングレート（０．５ｋＨｚ）でセンサー部１０の出力データをサンプリングする。また、時刻ｔ０〜ｔ１の期間と時刻ｔ４〜ｔ５の期間は、センサー部１０の位置が既知であり、第１期間に相当する。また、時刻ｔ２〜ｔ３の期間が運動解析の対象となる第２期間に相当する。
【００９８】
また、この実験例では、図８に示すように、テニスラケット３００のトップ３０２が最大速度の時（テニスラケット３００のフェイスにテニスボール４００が当たる直前）のセンサー部１０の位置を原点、トップ３０２の最大速度の方向をＸ軸として右手系でＹ軸、Ｚ軸を決める。そして、第２期間（時刻ｔ２〜ｔ３の期間）におけるＸＹＺ座標系でのテニスラケット３００のトップ３０２とグリップエンド３０４の軌跡をグラフ化して表示する。
【００９９】
図９は、処理部２００がセンサー部１０の出力データの取得を開始してから、第２期間におけるＸＹＺ座標系でのテニスラケット３００のトップ３０２とグリップエンド３０４の軌跡をグラフ表示するまでの処理のフローチャート図である。
【０１００】
まず、データ取得期間が終了するまで（ステップＳ１２０のＮ）センサー部１０から新たな３軸加速度データと３軸角速度データを周期的に取得する（ステップＳ１１０）。
【０１０１】
次に、データ取得期間が終了すると（ステップＳ１２０のＹ）、ステップＳ１１０で２つの第１期間（時刻ｔ０〜ｔ１の期間と時刻ｔ４〜ｔ５の期間）に取得した３軸角速度データの真値（０）に対する誤差を算出し、３軸ジャイロセンサーの出力誤差（角速度誤差）の時間関数を推定する（ステップＳ１３０）。例えば、角速度誤差の時間関数を多項式で近似して推定してもよい。
【０１０２】
次に、ステップＳ１３０で推定した時間関数を用いて、ステップ１１０で取得した３軸角速度データの誤差を除去して積分し、センサー部１０のＸＹＺ座標系での姿勢を算出する（ステップＳ１４０）。
【０１０３】
次に、ステップＳ１４０で算出したセンサー部１０のＸＹＺ座標系での姿勢を用いて、ステップＳ１１０で取得した３軸加速度データ（ｘｙｚ座標系の加速度ベクトル）をＸＹＺ座標系の加速度ベクトルに座標変換する（ステップＳ１５０）。
【０１０４】
次に、ステップＳ１５０の座標変換で得られたＸＹＺ座標系の加速度ベクトルを２階積分し、データ取得期間（時刻ｔ０〜ｔ５の期間）におけるＸＹＺ座標系でのセンサー部１０の位置を算出する（ステップＳ１６０）。
【０１０５】
次に、２つの第１期間（時刻ｔ０〜ｔ１の期間と時刻ｔ４〜ｔ５の期間）におけるＸＹＺ座標系でのセンサー部１０の位置の真値（第１の位置）に対する誤差を算出し、ＸＹＺ座標系の加速度ベクトルのＸ軸、Ｙ軸、Ｚ軸の各方向の加速度誤差の時間関数を推定する（ステップＳ１７０）。
【０１０６】
次に、ステップＳ１７０で推定した加速度誤差の時間関数を用いて、ＸＹＺ座標系の加速度ベクトルの誤差を除去して２階積分し、ＸＹＺ座標系でのセンサー部１０の位置（テニスラケット３００のグリップエンド３０４の位置）を算出する（ステップＳ１８０）。
【０１０７】
次に、センサー部１０からトップまでの距離と方向があらかじめ計測されていて既知であるとして、ステップＳ１６０で算出したＸＹＺ座標系でのセンサー部１０の位置とステップＳ１４０で算出したＸＹＺ座標系でのセンサー部１０の姿勢からＸＹＺ座標系でのテニスラケット３００のトップ３０２の位置を算出する（ステップＳ１９０）。
【０１０８】
最後に、運動解析の対象となる第２期間（時刻ｔ２〜ｔ３の期間）におけるＸＹＺ座標系でのテニスラケット３００のトップ３０２とグリップエンド３０４の位置の座標を抽出してグラフ表示する（ステップＳ２００）。
【０１０９】
図１０（Ａ）及び図１０（Ｂ）は、第２期間（時刻ｔ２〜ｔ３の期間）におけるテニスラケット３００のトップ３０２とグリップエンド３０４の軌跡の一例を示す図である。図１０（Ａ）はＸ−Ｙ平面上の軌跡であり、図１０（Ｂ）はＸ−Ｚ平面上の軌跡である。図１０（Ａ）では、Ｌ１で示す曲線がトップ３０２の軌跡であり、Ｌ２で示す曲線がグリップエンド３０４の軌跡である。また、図１０（Ｂ）では、Ｌ３で示す曲線がトップ３０２の軌跡であり、Ｌ４で示す曲線がグリップエンド３０４の軌跡である。図１０（Ａ）及び図１０（Ｂ）の軌跡は、実際のスイングの軌跡に適合するものになっている。
【０１１０】
比較のため、図１１（Ａ）及び図１１（Ｂ）のそれぞれには、図１０（Ａ）及び図１０（Ｂ）の軌跡に３軸加速度データの誤差を補正せずに積分する従来手法を適用した場合の軌跡を重ねて表示した図を示す。図１１（Ａ）において、軌跡グラフＧ１は図１０（Ａ）の軌跡グラフ（本実施形態の手法を適用した場合のＸ−Ｙ平面上の軌跡グラフ）であり、軌跡グラフＧ２は従来手法を適用した場合のＸ−Ｙ平面上の軌跡グラフである。また、図１１（Ｂ）において、軌跡グラフＧ３は図１０（Ｂ）の軌跡グラフ（本実施形態の手法を適用した場合のＸ−Ｚ平面上の軌跡グラフ）であり、軌跡グラフＧ４は従来手法を適用した場合のＸ−Ｚ平面上の軌跡グラフである。図１１（Ａ）及び図１１（Ｂ）から、従来手法を適用した場合の軌跡グラフＧ２とＧ４ではＸ軸方向の変位が約４ｍもあり、実際のスイングの軌跡と合っていないことが明らかである。この結果より、本実施形態の手法を適用することで、スイングの軌跡の精度が大幅に向上することがわかる。
【０１１１】
本発明は本実施形態に限定されず、本発明の要旨の範囲内で種々の変形実施が可能である。
【０１１２】
例えば、本実施形態では、加速度データを２階時間積分した位置のデータを補正する場合を例に挙げて説明したが、他の例としては、加速度データを１階時間積分した速度のデータを補正することもできる。その場合、例えば、第１期間を対象物の静止期間とすれば、第１期間では速度が０であるとして、加速度誤差の時間関数を推定すればよい。このように速度を補正することで、例えば、テニスラケット、ゴルフグラブ、バット等のスイング速度を精度よく計測することができる。また、他の例としては、ジャイロセンサーが出力する角速度データを１階時間積分した１軸回りの角度（回転角度）のデータを補正することもできる。その場合、例えば、第１期間を対象物の静止期間とすれば、第１期間では回転角度が０であるとして、角速度誤差の時間関数を推定すればよい。このように回転角度を補正することで、例えば、テニスラケットやゴルフグラブ等にボールが当たった直後の（インパクトの直後）の打球面の回転角度を精度よく計測することができる。
【０１１３】
本発明は、実施の形態で説明した構成と実質的に同一の構成（例えば、機能、方法及び結果が同一の構成、あるいは目的及び効果が同一の構成）を含む。また、本発明は、実施の形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施の形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施の形態で説明した構成に公知技術を付加した構成を含む。
【符号の説明】
【０１１４】
１  運動解析装置、１０  センサー部、２０  ホスト端末、１００  センサー、１０２ｘ，１０２ｙ，１０２ｚ  加速度センサー、１０４ｘ，１０４ｙ，１０４ｚ  ジャイロセンサー、１１０  データ処理部、１２０  通信部、２００  処理部（ＣＰＵ）、２０２  データ取得部、２０４  誤差時間関数推定部、２０６  データ補正部、２０８  運動解析情報生成部、２１０  通信部、２２０  操作部、２３０  ＲＯＭ、２４０  ＲＡＭ、２５０  不揮発性メモリー、２６０  表示部、３００  テニスラケット、３０２  トップ、３０４  グリップエンド、４００  テニスボール
【図１】
0005948011.tif 000022
【図２】
0005948011.tif 000023
【図３】
0005948011.tif 000024
【図４】
0005948011.tif 000025
【図５】
0005948011.tif 000026
【図６】
0005948011.tif 000027
【図７】
0005948011.tif 000028
【図８】
0005948011.tif 000029
【図９】
0005948011.tif 000030
【図１０】
0005948011.tif 000031
【図１１】
0005948011.tif 000032