"3軸方向の加速度及び3軸回りの角速度あるいは角度を計測可能なセンサが取り付けられたゴルフクラブでのスイングを計測するステップと、上記センサの計測結果に基づいて、スイング分類のための指標を得るステップと、上記指標に基づいてスイング分類するステップとを含み、上記指標が、以下のaまたはbを含むスイング分析方法。aダウンスイング中の少なくとも1時刻におけるグリップ角速度bダウンスイング中の少なくとも1時刻におけるグリップ速度請求項1に記載のスイング分析方法において、アドレス、トップ及びインパクトの時刻が決定される分析方法。上記指標が、以下のc、d及びeから選択される少なくとも2つを含む請求項1または2に記載のスイング分析方法。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。eゴルファーへ向かう方向のグリップ速度。上記指標が、以下のc、d及びfから選択される少なくとも2つである請求項1または2に記載のスイング分析方法。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。fヘッドの動的ロフト角。上記指標が、以下のc、d及びgから選択される少なくとも2つである請求項1または2に記載のスイング分析方法。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。gインパクトでのグリップ姿勢角度。請求項1から5のいずれかに記載のスイング分析方法を用いたゴルフクラブのフィッティング方法であって、上記スイング分類毎に推奨シャフトが定められているフィッティング方法。上記文献の計測装置は、大型且つ複雑である。このような計測装置は、店頭でフィッティングを行うためのハードウェアには適さない。モーションキャプチャシステムも実用化されているが、高価である。一方、複雑な動作であるスイングを適切に分類できれば、フィッティングが容易且つ正確となる。本発明の目的は、簡易な装置によってスイングを分類しうるスイング分析方法の提供にある。本発明に係るスイング分析方法は、3軸方向の加速度及び3軸回りの角速度あるいは角度を計測可能なセンサが取り付けられたゴルフクラブでのスイングを計測するステップと、上記センサの計測結果に基づいて、スイング分類のための指標を得るステップと、上記指標に基づいてスイング分類するステップとを含む。上記指標が、以下のaまたはbを含む。aダウンスイング中の少なくとも1時刻におけるグリップ角速度bダウンスイング中の少なくとも1時刻におけるグリップ速度好ましくは、上記分析方法において、アドレス、トップ及びインパクトの時刻が決定される。好ましくは、上記指標が、以下のc、d及びeから選択される少なくとも2つを含む。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。eゴルファーへ向かう方向のグリップ速度。他の好ましい態様は、上記指標が、以下のc、d及びfから選択される少なくとも2つを含む。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。fヘッドの動的ロフト角。他の好ましい態様は、上記指標が、以下のc、d及びgから選択される少なくとも2つを含む。cトウダウン方向のグリップ角速度。dシャフト軸回りのグリップ角速度。gインパクトでのグリップ姿勢角度。本発明にかかるフィッティング方法は、上述のいずれかに記載のスイング分析方法を用いたゴルフクラブのフィッティング方法である。このフィッティング方法では、上記スイング分類毎に推奨シャフトが定められている。このスイング分析方法では、スイング分類が簡易な装置によってなされうる。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。図1は、本発明の一実施形態を示す図である。本実施形態では、スイング分析システム2が用いられる。スイング分析システム2は、センサ4と、コンピュータに内蔵された無線受信装置6と、データ解析装置8とを備えている。図1では、センサ4がゴルファーg1の腕によって隠されている。図2は、ゴルフクラブc1のグリップエンドに取り付けられたセンサ4を示す斜視図である。図3aは、センサ4の平面図である。図3bは、センサ4の正面図である。センサ4は、グリップ7の後端に取り付けられている。センサ4は、アダプター9を介して、グリップ7に取り付けられている。センサ4の取り付けは、両面テープ、接着剤、ねじ止め等によって容易になされうる。センサ4は、小型且つ軽量である。センサ4は、無線式である。センサ4は、無線により、測定データを発信しうる。この無線通信の詳細については後述される。無線通信として、例えばBluetoothブルートゥースの規格及び技術が好適に用いられ得る。センサ4は、3軸方向x1軸、y1軸及びz1軸のそれぞれの加速度を測定できる加速度センサーを内蔵している。またセンサ4は、3軸方向x1軸、y1軸及びz1軸周りの角速度を計測できる角速度センサを内蔵している。更にセンサ4は、A/D変換器、CPU、無線インターフェース、無線アンテナ及び電源を備えている。電源として、電池が用いられている。電池は、例えばリチウムイオン電池等の小型電池が好適に用いられる。いわゆるボタン型電池が好適に用いられ得る。電池は、充電可能なものであってもよい。センサ4は、電池を充電するための充電回路を備えていてもよい。センサ4の一例として、ワイヤレステクノロジー社製の商品名WAA-010が挙げられる。なお、上記角速度センサに代えて、角度センサが用いられてもよい。即ち、3軸回りの角度を計測可能なセンサが用いられても良い。図示しないが、無線受信装置6は、無線アンテナ、無線インターフェース、CPU及びネットワークインターフェースを備えている。データ解析装置8として、例えば、コンピュータが用いられる。データ解析装置8は、入力部12及び表示部14を備えている。図示しないが、データ解析装置8は、ハードディスク、メモリ、CPU及びネットワークインターフェースを備えている。入力部12は、キーボード16とマウス18とを備えている。図1には、スイング分析システム2の他、ゴルファーg1と、ゴルフクラブc1と、ゴルフボールb1とが示されている。図1に描かれたゴルファーg1は、アドレス状態である。ゴルファーg1は、右利きである。センサ4は、x1軸、y1軸及びz1軸のそれぞれの方向の加速度を検知する。図2が示すように、x1軸、y1軸及びz1軸は、3軸直交座標系を構成している。これらの加速度は、アナログ信号として得られる。このアナログ信号は、A/D変換器によって、デジタル信号に変換される。A/D変換器からの出力は、例えばCPUに伝達されて1次フィルタリング等の演算処理が実行される。このようにセンサ4内で処理されたデータは、無線インタフェースを介して、無線アンテナから送信される。センサ4は、x1軸、y1軸及びz1軸の各軸周りの角速度を検知する。これらの角速度は、アナログ信号として得られる。このアナログ信号は、A/D変換器によって、デジタル信号に変換される。A/D変換器からの出力は、例えばCPUに伝達されて1次フィルタリング等の演算処理が実行される。このようにセンサ4内で処理されたデータは、無線インタフェースを介して、無線アンテナから送信される。センサ4の無線アンテナから送信されたデータは、無線受信装置6側の無線アンテナを介して、無線インタフェースによって受信される。この受信されたデータは、例えばCPUで演算処理され、データ解析装置8に送られる。データ解析装置8に送られたデータは、ハードディスク等のメモリ資源に記録される。ハードディスクは、データ処理等に必要なプログラム及びデータ等を記憶している。このプログラムは、CPUに、必要なデータ処理を実行させる。CPUは、各種の演算処理を実行しうる。演算処理の例については、後述される。演算結果は、表示部14または図示しない印刷装置等によって出力される。センサ4の取り付けにおいては、計測軸とゴルフクラブとの関係が考慮される。本実施形態では、センサ4のz1軸が、ゴルフクラブc1のシャフト軸Zsに一致している。センサ4のx1軸は、ヘッドのトウ-ヒール方向にできるだけ沿うように配向される。センサ4のy1軸は、フェース面の法線方向にできるだけ沿うように配向される。このようにセンサ4を取り付けることで、演算が簡略化される。[全体座標系と局座標系]本願では、全体座標系及び局座標系が考慮される。説明の便宜のため、全体座標系の各軸を、大文字のX軸、Y軸及びZ軸で表すこととし、局座標系の各軸を、小文字のx軸、y軸及びz軸で表すこととする。全体座標系のX軸、Y軸及びZ軸は、三次元直交座標系である。解析を容易とするため、本実施形態では、Z軸が鉛直方向とされ、X軸がゴルファーg1の前後方向とされる。Y軸は目標方向とされる。より詳細には、Y軸は、打球地点と目標地点とを結び且つ地面と平行な方向とされる。局座標系のx軸、y軸及びz軸は、三次元直交座標系である。解析を容易とするため、本実施形態では、z軸が上記シャフト軸Zsとされる。x軸は、ヘッドのトウ-ヒール方向にできるだけ沿うように設定される。y軸は、フェース面の法線方向にできるだけ沿うように設定される。局座標系のz軸は、センサ4のz1軸に一致している。局座標系のy軸は、センサ4のy1軸に一致している。局座標系のx軸は、センサ4のx1軸に一致している。この設定により、下記の6つのデータが、測定値として直接得られる。センサ4の測定値として直接得られるデータは、次の6つである。x軸回りの角速度ωxグリップ角速度ωxy軸回りの角速度ωyグリップ角速度ωyz軸回りの角速度ωzグリップ角速度ωzx軸方向の加速度Axy軸回りの加速度Ayz軸回りの加速度Azセンサ4により、時系列的に連続した複数のデータが得られる。単位時間当たりのデータ数は、サンプリング周波数に依存する。アドレス状態のゴルファーg1においては、全体座標系のY軸と、局座標系のy軸とは、ほぼ平行である。全体座標系に対する局座標系の向きは、スイング中、時々刻々と変化する。局座標系の向きは、ゴルフクラブc1の動きに連動して変化する。図4から図7は、スイングを説明するための図である。図4から図7は、ゴルファーg1の正面前側から見た図である。スイングの始まりは、アドレスである。スイングの終わりは、フィニッシュと称される。スイングは、S1、S2、S3、S4、S5、S6、S7、S8の順で進行する。S1及びS2が、図4に示されている。S3及びS4が、図5に示されている。S5及びS6が、図6に示されている。S7及びS8が、図7に示されている。図4のS1は、アドレスである。図4のS2は、テークバックである。図5のS3は、トップトップオブスイングである。通常、トップオブスイングでは、スイング中におけるヘッドの移動速度が最小である。図5のS4は、ダウンスイングである。図6のS5も、ダウンスイングである。S5は、S4よりもダウンスイングが進行した状態である。図6のS6は、インパクトである。インパクトは、ゴルフクラブc1のヘッドとゴルフボールb1とが衝突した瞬間である。図7のS7は、フォロースルーである。図7のS8は、フィニッシュである。フィニッシュで、スイングは終了する。スイング分析システム2を用いるための準備として、スイング分類のための指標を探索するステップがなされた。この指標の探索では、モーションキャプチャシステムが用いられた。この指標の探索では、20名のゴルファーのスイングが分析され、43種類のスイング特徴が定量化された。また別途、これら20名のゴルファーのそれぞれに5本のクラブを試打させ、飛距離、打球の曲がりにくさ及び振りやすさに関する評価データを取得した。そして、上記43種類のスイング特徴と上記評価データとに基づき、統計的な解析がなされた。この統計的な解析では、判別分析と称される手法が用いられた。この結果、上記43種類のスイング特徴のうちのいくつかが、スイング分類の指標として有効であることが判明した。特に、以下の3つの指標は、統計的にみて70以上のゴルファーに有効であるという知見を得た。指標1:ダウンスイング中における、トウダウン方向の角速度指標2:ダウンスイング中における、シャフト軸Zs回りの角速度指標3:インパクト近傍好ましくはインパクト時における、ゴルファーへ向かう方向のグリップ速度指標1の角速度は、上記y軸回りの角速度に対応する。よってこの角速度が、本願においてωyと表記される。指標2の角速度は、上記z軸回りの角速度に対応する。よってこの角速度が、本願においてωzと表記される。これらの有効な指標を簡便に計測するためのシステムが検討され、上記スイング分析システム2を得た。このスイング分析方法では、センサ4からのデータのみで、スイング分析が可能である。スイング分析システム2を用いたスイング分析方法は、次のステップst1からst3を含む。3軸方向の加速度及び3軸回りの角速度を計測可能なセンサ4が取り付けられたゴルフクラブc1がスイングされるステップst1上記センサの計測結果に基づいて、スイング分類のための指標を得るステップst2上記指標に基づいてスイング分類するステップst3上記指標は、以下のaまたはbを含む。aダウンスイング中の少なくとも1時刻におけるグリップ角速度bダウンスイング中の少なくとも1時刻におけるグリップ速度ダウンスイング中とは、アドレスからインパクトまでの間の時間である。このましくは、このスイング分析方法では、上記センサの計測結果に基づいて、アドレス、トップ及びインパクトの時刻が決定される。これらの時刻が把握されることで、データを取得すべき時刻が容易に把握されうる。好ましくは、上記指標は、以下のc、d及びeから選択される少なくとも2つを含む。より好ましくは、以下のc、d及びeの全てを含む。cトウダウン方向のグリップ角速度ωydシャフト軸回りのグリップ角速度ωzeゴルファーへ向かう方向のグリップ速度Vx上記cは、少なくとも1時刻における値である。スイング分類の精度を高める観点から、好ましくは、上記cは平均値とされる。平均値は、時刻t1からt2までの間での計測値を平均することで得られる。時刻t1及びt2は、ダウンスイング中の時刻である。時刻t1がトップとされ、時刻t2がインパクトとされてもよい。上記dは、少なくとも1時刻における値である。スイング分類の精度を高める観点から、好ましくは、上記dは平均値とされる。平均値は、時刻t1からt2までの間での計測値を平均することで得られる。時刻t1及びt2は、ダウンスイング中の時刻である。時刻t1がトップとされ、時刻t2がインパクトとされてもよい。スイング分類の精度を高める観点から、好ましくは、上記eは、インパクト近傍における値である。インパクト近傍とは、インパクトの0.01秒前からインパクトまでを意味する。より好ましくは、上記eは、インパクト時における値である。ゴルファー方向のグリップ速度Vxとは、x軸方向の速度である。このグリップ速度Vxは、グリップとゴルファーとの距離変化の指標となりうる。ゴルファーによっては、インパクト近傍において、グリップを自分の体の方向に引き込むような動作をする。この動作の度合いにより、グリップ速度Vxが変化しうる。なお、グリップ速度Vxに代えて、グリップとゴルファーとの距離変化を示しうる他の速度が採用されてもよい。図8は、上記c、d及びeを用いたスイング分類ステップst3を示すフローチャートである。このスイング分類では、上記cのグリップ角速度ωyが判断されるステップst10。このグリップ角速度ωyが500deg/s未満である場合、スイングパターンAに分類される。このグリップ角速度ωyが500deg/s以上600deg/s未満である場合、スイングパターンBに分類される。このグリップ角速度ωyが600deg/s以上である場合、次のステップst20へと進む。500deg/s及び600deg/s以外の値が用いられてもよい。即ち閾値は適宜設定されうる。なお、この実施形態では、グリップ角速度ωyとして、ダウンスイング全体での平均値が採用されている。次に、上記dのグリップ角速度ωzが判断されるステップst20。このグリップ角速度ωzが300deg/s以上である場合、スイングパターンEに分類される。このグリップ角速度ωzが300deg/s未満である場合、次のステップst30へと進む。300deg/s以外の値が用いられてもよい。即ち閾値は適宜設定されうる。なお、この実施形態では、グリップ角速度ωzとして、ダウンスイング全体での平均値が採用されている。次に、上記eのグリップ速度Vxが判断されるステップst30。このグリップ速度Vxが2.0m/s未満である場合、スイングパターンCに分類される。このグリップ速度Vxが2.0m/s以上である場合、スイングパターンDに分類される。2.0m/s以外の値が用いられてもよい。即ち閾値は適宜設定されうる。なお、この実施形態では、グリップ速度Vxとして、インパクト時の値が採用されている。以上のようにして、スイングが、5つのパターンAからEに分類されうる。他の好ましい態様では、上記指標は、以下のc、d及びfから選択される少なくとも2つを含む。より好ましくは、以下のc、d及びfの全てを含む。cトウダウン方向のグリップ角速度ωydシャフト軸回りのグリップ角速度ωzfヘッドの動的ロフト角DL。動的ロフト角DLとは、鉛直方向に対するフェース面の傾斜角度である。スイング分類の精度を高める観点から、上記fは、好ましくはインパクト近傍における値であり、より好ましくはインパクト時における値である。動的ロフト角DLは、例えば、インパクト近傍におけるヘッド画像に基づいて測定される。図9は、上記c、d及びfを用いたスイング分類ステップst3を示すフローチャートである。ステップst10及びステップst20は、前述した図8の実施形態と同じである。この実施形態では、上記ステップst30に代えて、ステップst40がなされる。このステップst40では、インパクト時の動的ロフト角DLが判断される。この動的ロフト角DLが15degree未満である場合、スイングパターンCに分類される。この動的ロフト角DLが15degree以上である場合、スイングパターンDに分類される。15degree以外の値が用いられてもよい。即ち閾値は適宜設定されうる。本発明に係るフィッティング方法では、前述したスイング分析方法を用いる。このフィッティング方法では、上記スイング分類毎に推奨シャフトが定められている。このフィッティング方法では、スイング分類することによって、推奨シャフトが決定される。一般に、フィッティングにおいては、シャフトの選定が難しい。このフィッティング方法では、シャフトの選定が容易になされうる。次に、各指標の算出方法について説明する。[グリップ角速度]グリップ角速度は、センサ4により得られた角速度から直接的に得られうる。上記実施形態では、センサ4のy1軸回りの角速度が、グリップ角速度ωyに等しい。また、センサ4のz1軸回りの角速度が、グリップ角速度ωzに等しい。センサ4のx1軸回りの角速度が、x軸回りの角速度グリップ角速度ωxに等しい。[グリップ傾斜角度α]後述される姿勢行列を求めることによって、グリップ傾斜角度αが決定される。姿勢行列から角度αを計算する方法は、後述される。角度αは、正面から見たときの、Y軸とシャフト軸線の成す角度である。[グリップ速度]通常、加速度を積分することによって、速度が求められる。センサ4では、3軸方向の加速度が計測されるので、これらの加速度を積分することが考えられる。しかしこの単純な積分では、グリップ速度を算出することはできない。センサ4で計測された加速度には、次の3つの加速度が含まれている。加速度1並進成分の加速度加速度2遠心成分の加速度加速度3重力加速度グリップ速度並進速度を求めるには、並進成分の加速度のみを積分することが必要である。そのためには、遠心成分の加速度及び重力加速度を排除する必要がある。グリップ速度の算出では、姿勢行列が用いられる。図10は、この姿勢行列を示す。この姿勢行列はaからiまでの9成分を有している。この姿勢行列には、グリップの姿勢が記述されている。この姿勢行列における各成分の意味は、次の通りである。成分a:全体座標系のX軸と、局座標系のx軸とのなす角度のコサインcos成分b:全体座標系のY軸と、局座標系のx軸とのなす角度のコサインcos成分c:全体座標系のZ軸と、局座標系のx軸とのなす角度のコサインcos成分d:全体座標系のX軸と、局座標系のy軸とのなす角度のコサインcos成分e:全体座標系のY軸と、局座標系のy軸とのなす角度のコサインcos成分f:全体座標系のZ軸と、局座標系のy軸とのなす角度のコサインcos成分g:全体座標系のX軸と、局座標系のz軸とのなす角度のコサインcos成分h:全体座標系のY軸と、局座標系のz軸とのなす角度のコサインcos成分i:全体座標系のZ軸と、局座標系のz軸とのなす角度のコサインcos3つの成分a,b,cは、局座標系のx軸ベクトルの向きを示している。3つの成分d,e,fは、局座標系のy軸ベクトルの向きを示している。3つの成分g,h,iは、局座標系のz軸ベクトルの向きを示している。スイング中、グリップの姿勢は時々刻々と変化する。よって、上記姿勢行列も、スイング中、時々刻々と変化する。クリップ速度の算出方法では、遠心成分の加速度及び重力加速度を除いた上で、加速度の積分がなされる。この算出方法では、先ず、アドレスでの姿勢行列が算出されるステップst100。次に、各時刻における姿勢行列が算出されるステップst200。次に、計測された加速度データが、全体座標系に変換されるステップst300。この変換により、遠心成分の加速度が除かれる。次に、アドレス時の加速度がオフセットされるステップst400。このオフセットにより、重力加速度が除かれる。最後に、遠心成分の加速度及び重力加速度が除かれた加速度を積分するステップst500。この積分により、グリップ速度が得られる。ステップst100では、重力が鉛直方向に作用していることを利用して、アドレス時の姿勢行列が算出されうる。このステップst100では、先ず、姿勢行列の成分c、成分f及び成分iが求められるステップst110。重力は、鉛直方向、即ち、全体座標系のZ軸方向に作用しているので、成分c、成分f及び成分iは、センサーによる加速度データに基づき直接的に算出されうる。次に、求められた成分c、成分f及び成分iを用いて、他の6成分が算出されるステップst120。この6成分も、上記成分c、f及びiと同様に、幾何学的演算によって算出されうる。また、この6成分の算出では、重回帰分析が用いられうる。また、この6成分の算出に、オイラー変換行列が用いられてもよい。図11は、オイラー変換行列を示す。本実施形態では、上記6成分のうち、成分a及び成分gは、重回帰分析によって高い精度で推定することができた。この推定の精度を示す重相関係数は0.98であった。一方、他の4成分b、d、e及びhは、重回帰分析での推定では高い精度が得られなかった。成分bについては、オイラー変換行列を用いて推定した。即ち、オイラー変換行列と、既知の値である成分a,c,f,g及びiとを対比して、オイラー角度であるθ、φ及びψを求め、これらオイラー角度に基づいて成分bが求められた。モーションキャプチャにより求めた成分bと比較したところ、相関係数は0.86であった。なお、上記オイラー角度のうち、ψのプラスマイナスは、角度情報のみからは決定できない。このプラスマイナスは、例えば、アドレス時におけるグリップ角度の計算結果に基づいて決定される。グリップ角度の計算方法については、後述される。その他の成分のうち、成分hは、x軸とz軸との成す角度が90度であることを利用して計算した。また、成分d及びeは、他の成分の値を用いて、次の式により算出した。d=bi-che=ai-cg以上に基づいて算出された姿勢行列の精度を検証した。この検証は、上記姿勢行列のグリップ姿勢と、モーションキャプチャによって求めたグリップ姿勢との比較により行った。5名のテスターテスターA、B、C、D及びEの計測データに基づいて検証がなされた。結果が図12から14のグラフで示される。図12は、アドレスでのグリップ角度αに関するグラフである。図13は、アドレスでのグリップ角度βに関するグラフである。図14は、アドレスでのグリップ角度γに関するグラフである。図12から14において、横軸は上記姿勢行列に基づき算出された角度であり、縦軸はモーションキャプチャに基づき算出された角度である。いずれのグラフにおいても、縦軸の値と横軸の値とはおおよそ一致している。ほとんどのデータにおいて、誤差は±3degree以内に収まっている。よって、上記姿勢行列の精度が高いと考えられる。なお、グリップ角度α、β及びγについては、後述される。ステップst200では、アドレス時の姿勢行列N1と回転行列T1との積が求められるステップst210。この積により、姿勢行列N2が得られる。即ち、N2=N1T1である。図15は、この回転行列T1を示す。回転行列T1は、アドレス時から、最初のサンプリング時刻T1までの角度変化量により定まる。この角度変化量は、サンプリング時刻T1におけるセンサ4の角速度ωx1、ωy1及びωz1のそれぞれと、時間Δtとの積である。時間Δtは、計測サンプリングの単位時間に相当する。このΔtは、サンプリング時間とも称される。回転行列T1において、θx1=ωx1Δtであり、θy1=ωy1Δtであり、θz1=ωz1Δtである。次に、上記姿勢行列N2と回転行列T2との積が求められるステップst220。この積により、姿勢行列N3が得られる。即ち、即ち、N3=N2T2=N1T1T2である。図16は、この回転行列T2を示す。回転行列T2は、前サンプリング時刻T1から、このサンプリング時刻T2までの角度変化量により定まる。この角度変化量は、サンプリング時刻T2におけるセンサ4の角速度ωx2、ωy2及びωz2のそれぞれと、時間Δtとの積である。回転行列T2において、θx2=ωx2Δtであり、θy2=ωy2Δtであり、θz2=ωz2Δtである。以後、同様にして、各時刻の姿勢行列Niと、回転行列Tiとの積が繰り返されるステップst230。この繰り返しにより、全てのサンプリング時刻における姿勢行列Niが求められる。ステップst300では、各サンプリング時刻の加速度Aiと、その時刻の姿勢行列Niとの積が求められるステップst310。加速度Aiは、センサ4から得られるデータそのものである。これらの積により、各サンプリング時刻における加速度Aiが全体座標系の加速度に変換される。この変換により、遠心成分の加速度が除かれたことになる。ステップst400では、アドレス時の加速度として、データ初期の100データの平均値が採用される。この100のデータの選択では、例えば、後述の方法で決定されるアドレス時刻が用いられる。このアドレス時刻の直後の50データと、このアドレス時刻の直前の50データとが用いられ得る。平均されるデータの数は、適宜設定される。ステップst500では、、遠心成分の加速度及び重力加速度が除かれた加速度の波形が積分される。積分時間は、例えば、2.0秒とされる。センサ4には、ドリフト誤差と呼ばれる誤差が生じうる。ドリフト誤差は、校正終了後に発生する測定器または測定システムの性能変化に起因する。ドリフト誤差の主な原因として、測定器内部の接続ケーブルの熱膨張、及び、周波数変換器の温度ドリフトが挙げられる。周囲温度の変化に応じて頻繁に校正を実行したり、周囲温度をできるだけ一定に保ったりすることで、ドリフト誤差は低減されうる。しかし、これらの対策によっては、ドリフト誤差を完全に排除することは難しい。ドリフト誤差は、基本的に、低周波成分である。このため、ステップst500の積分において完全積分が用いられた場合、低周波での振幅ゲインが増大する。即ちこの場合、積分によりドリフト誤差が増大しやすい。ドリフト誤差を抑制する観点から、ステップst500の積分では、不完全積分が用いられる。ここで、完全積分及び不完全積分について説明する。例えば、変位Xが[X=Asinωt]で表されるとき、速度Vは、このXを微分して、次の式で表される。V=X’=Aωcosωt即ち、微分すると、振幅ゲインはω倍となる。よって、逆に積分すると、振幅ゲインが1/ω倍になる。ω=2πfであるから、積分すると振幅は1/2πf倍となる。これが通常の積分であり、完全積分とも称される。ゴルフスイングにおいて、テイクバックからダウンスイングまでの運動に要する時間は、概ね2秒程度である。よってこの運動の周波数は、概ね0.5Hzということになる。そこで、0.5Hz以上の高周波成分については完全積分を行い、一方、0.5Hz未満の低周波成分については、積分による振幅の増大が抑制されるのが好ましい。この場合、低周波成分であるドリフト誤差が効果的に抑制されうる。この積分による振幅の増大を抑制するために、不完全積分が用いられる。例えば、0.5Hz未満の低周波成分について、不完全積分が用いられる。加速度がakとされ、速度がvkとされ、1ステップ後の速度がvk1とされ、Δtがサンプリング時間とされ、kが積分の刻みステップとされ、Tiが積分時間とされるとき、この不完全積分では、vk1が次の式1で求められる。vk1=pvkqak1ただし、a0=v0である。式1において、pはフィルタ係数と称され、次の式により算出される。p=exp-Δt/Ti=exp-2πfiΔt式1において、qはフィルタ入力係数と称され、次の式により算出される。q=1-pTi=1-p/2πfi次に、アドレス、トップ及びインパクトの時刻を決定する方法について説明する。この方法は、3つの段階に分けられる。第1段階では、インパクト時刻が決定される。第2段階では、インパクト時刻に基づいてトップ時刻が決定される。第3段階では、トップ時刻に基づいて、アドレス時刻が決定される。上記第1段階では、先ず、仮インパクト時刻Ttが決定されるステップst1000。このステップst1000では、計測されたデータの中で、グリップ角速度ωzが最大値となっている時刻が、仮インパクト時刻とされる。次に、上記仮インパクト時刻Ttの10msec前の時刻Tpと、上記仮インパクト時刻の10msec後の時刻Tfとが決定されるステップst1010。次に、上記時刻Tpから上記時刻Tfの間で、グリップ角速度ωxが最小となる時刻Taが決定されるステップst1020。また、上記時刻Tpから上記時刻Tfの間で、z軸回りの加速度Azが最小となる時刻Tbが決定されるステップst1030。次に、上記時刻Taと上記時刻Tbとを比較して、早いほうの時刻が、インパクト時刻とされるステップst1040。図17は、センサ4による測定結果の一例を示すグラフである。図17の上側のグラフは、角速度ωx、角速度ωy及び角速度ωzの測定結果である。丸いプロットは、角速度ωxを示す。四角のプロットは、角速度ωyを示す。三角のプロットは、角速度ωzを示す。図17の下側のグラフは、加速度Ax、加速度Ay及び加速度Azの測定結果である。丸いプロットは、加速度Axを示す。四角のプロットは、加速度Ayを示す。三角のプロットは、加速度Azを示す。図17において、縦に延びる太い実線は、仮インパクト時刻Ttを示す。図17において、縦に延びる太い破線は、上記時刻Tpを示す。図17において、上記時刻Tfは示されていない。図17の実施例では、時刻Tbが、時刻Taよりも早い。よってこの実施例では、時刻Tbがインパクト時刻とされる。第2段階の、トップ時刻の決定では、グリップ角速度ωyが用いられる。第2段階では、インパクト時刻の0.5秒前からインパクト時刻までの間で、グリップ角速度ωyがゼロとなっている時刻が、トップ時刻とされるステップst1050。第3段階では、アドレス時刻が決定される。通常、ゴルファーは、アドレスに、ある程度の時間をかける。これに対して、ここでアドレス時刻を決定する意義は、スイング分析に最適なアドレス時刻を決定することにある。また、アドレス時刻の決定により、スイング分析のための適切な分析開始時刻を決定することができる。この第3段階では、グリップ角速度ωxが用いられる。図18は、この第3段階アドレス時刻決定方法を示すフローチャートである。好ましいアドレス時刻決定方法は、以下のステップst1100から1210から選ばれる少なくとも1つを含むこのアドレス時刻決定方法では、計測開始からトップ時刻までの間で角速度ωxが最大となる時刻Tmaxが決定されるステップst1100。次に、この時刻Tmaxを初期値とする時刻ttが設定されるステップst1110。また、固定値としての時間Tsが設定されるステップst1110。図18の実施形態では、時間Tsは50msである。50ms以外の値が採用されてもよい。時刻ttにおける角速度ωxがωxttである。この時刻tt-Tsにおける角速度ωxがωxtt-Tsである。[ωxtt-ωxtt-Ts]が1.0deg/s以下であるか否かが判断されるステップst1120。即ち、時刻ttからさかのぼり、時間Tsにおける角速度ωxが所定値以下であるか否かが判断される。1.0deg/s以外の値が採用されてもよい。ステップst1120がいいえNoである場合、時刻ttが、1ms前の時刻に更新されるst1130。1ms以外の値が採用されてもよい。更新された時刻ttが、0より小さいか否かが判断されるステップst1140。時刻0とは測定開始時刻である。時刻ttが0より小さい場合、アドレス時刻における計測がなされていないとみなされる。この場合、再計測がなされるステップst1150。ステップst1140において、時刻ttが0以上である場合、判断ステップst1120に戻る。ステップst1120がはいYESである場合、角速度ωxは小さい。角速度ωxが小さい状態は、アドレス状態に近い。よって、この場合の時刻ttが、仮アドレス時刻Ttに決定されるステップst1160。一方、本実施形態では、角速度ではなく加速度のデータが用いられる。x軸方向の加速度Ax、y軸回りの加速度Ay及びz軸回りの加速度Azを用いて、合成加速度Amを求めるステップst1170。合成加速度Amの計算式は次の通りである。Am=Ax2Ay2Az21/2次に、上記時刻Ttからさかのぼり、合成加速度Amが重力加速度近傍で安定している最初の時間帯Tzが決定されるステップst1180。具体的には、上記時刻Ttから時間をさかのぼり、合成加速度Amについて、連続した30個のデータが9.7以上10.5以下である最初の時間帯が決定される。連続したデータの数は、30以外であってもよい。9.7以外の下限値が採用されてもよい。10.5以外の上限値が採用されてもよい。次に、この時間帯Tzに基づいて、時刻t1が決定されるステップst1190。時刻t1は、例えば、上記時間帯Tzにおける最初の時刻でもよいし、上記時間帯Tzにおける最後の時刻でもよいし、上記時刻帯Tzの中心時刻であってもよい。9.7以上10.5以下の範囲は、重力加速度を除くと、ほぼゼロに近い。この数値範囲では、重力加速度以外の加速度がゼロに近い。これは、アドレスに近い状態である。次に、y軸回りの加速度Ayについて、移動平均波形Wsが作成されるステップst1200。具体的には、加速度Ayの波形について、測定開始時から上記時刻t1まで50点移動平均し、これを更に50点移動平均して、滑らかな波形Wsを得る。移動平均の処理数は2回以外であってもよい。移動平均の方法は、50点以外であってもよい。次に、上記波形Wsについて上記時刻t1から時間をさかのぼり、区間差が0.01以下となる最初の時刻tiが決定されるステップst1210。具体的な区間差は、時刻tiと時刻ti-20との差である。時刻ti-20とは、時刻tiから、20データ分だけ時間をさかのぼった時刻である。この区間差が0.01以下となる最初の時刻tiが、アドレス時刻t2とされる。0.01以外の閾値が採用されてもよい。20以外のデータ数が採用されてもよい。このようにして、アドレス時刻t2が決定される。アドレス時のデータを選択する場合、好ましくは、このアドレス時刻t2が基準とされる。好ましくは、アドレス時のデータを選択する場合、このアドレス時刻t2近傍の複数のデータが平均される。好ましい具体例では、アドレス時刻t2から、時刻t2-50の範囲の50個のデータが平均される。時刻t2-50とは、時刻t2から、50データ分だけ時間をさかのぼった時刻である。50以外の個数が採用されてもよい。次に、グリップ角度α、β及びγについて説明する。。図19は、グリップ角度αを示す。グリップ角度αは、Y軸とグリップ軸シャフト軸との成す角度である。グリップ角度αは、正面から見たときの2次元画像における角度である。このグリップ角度αは、前述したグリップ傾斜角度に等しい。図20は、グリップ角度βを示す。グリップ角度βは、X軸とグリップ軸シャフト軸との成す角度である。グリップ角度βは、目標方向の後方から見たときの2次元画像における角度である。図21は、グリップ角度γを示す。グリップ角度γは、Y軸とグリップ軸シャフト軸との成す角度である。グリップ角度γは、上方から見たときの2次元画像における角度である。グリップ角度α、β及びγは、姿勢行列に基づいて計算されうる。図10で示される姿勢行列の成分f及びiを用いて、グリップ角度αは、次のように計算されうる。α=arctanf/-i図10で示される姿勢行列の成分c及びiを用いて、グリップ角度βは、次のように計算されうる。β=arctanc/-i図10で示される姿勢行列の成分c及びfを用いて、グリップ角度γは、次のように計算されうる。γ=arctanc/f図22は、図8及び図9で示したスイング分類ステップst3の変形例である。この例では、以下のc、d及びgが用いられている。c上記角速度ωyd上記角速度ωzgインパクトでのグリップ姿勢角度上記gの、グリップ姿勢角度として、上記グリップ角度α、β及びγが挙げられる。好ましくは、グリップ角度αが用いられる。図22のフローチャートにおいて、ステップst10及びステップst20は、前述した図8の実施形態と同じである。この実施形態では、上記ステップst30に代えて、ステップst50がなされる。このステップst50では、インパクト時のグリップ角度αが判断される。この角度αが2degree未満である場合、スイングパターンCに分類される。この角度αが2degree以上である場合、スイングパターンDに分類される。2degree以外の値が用いられてもよい。即ち閾値は適宜設定されうる。センサ4の大きさ及び重量は限定されない。スイングを妨げない観点から、センサ4は、小型でかつ軽量であるのが好ましい。この観点から、センサ4の重量は、10g以下であるのが好ましく、6g以下であるのがより好ましい。以下、実施例によって本発明の効果が明らかにされるが、この実施例の記載に基づいて本発明が限定的に解釈されるべきではない。16名のゴルファーが、テストを行った。ゴルファーは、全員上級者とされた。全てのゴルファーが、ドライバーでのヘッドスピードが42m/sであり、普段使用しているドライバーの重量が300g以上であった。5本のテストクラブと基準クラブとが用意された。全てのクラブで、ヘッド及びグリップは共通とされた。いずれのクラブも、番手はドライバーとされ、スイングバランスはD2とされ、ロフト角は9.5度とされた。クラブの総重量は、307g以上323g以下であった。テストクラブに装着されているシャフトは、シャフトA、シャフトB、シャフトC、シャフトD及びシャフトEの5種類である。これらのシャフトAからEの硬さフレックスは、Sである。モーションキャプチャシステムを用いた事前検証により、シャフトとスイングパターンとの適合性が確認されていた。シャフトAは、スイングパターンAのゴルファーに適している。シャフトBは、スイングパターンBのゴルファーに適している。シャフトCは、スイングパターンCのゴルファーに適している。シャフトDは、スイングパターンDのゴルファーに適している。シャフトEは、スイングパターンEのゴルファーに適している。図2と同様にして、基準クラブにセンサが取り付けられた。この基準クラブを用いて、図8のフローチャートに基づき、スイング分類がなされた。このスイング分類により、適合するシャフトが決定された。次に、この適合シャフトと、その適合シャフトにスペックが比較的近い2本のシャフトとが選定され、これら3本のシャフトの比較がなされた。たとえば、あるテスターXの場合、スイング分類の結果が、スイングパターンBであり、適合シャフトはシャフトBであった。このシャフトBにスペックが比較的近いのは、シャフトA及びCであった。そこで、テスターXは、3本のシャフトA、B及びCを比較し、最も良好なシャフトを選定した。この選定の基準は、飛距離、打球方向性及び振りやすさとされた。最も良好なシャフトがシャフトBである場合、正解とされた。最も良好なシャフトがシャフトAまたはCを選んだ場合、不正解とされた。このテストの結果は次の通りであった。飛距離に関して正解であったのは、16名中、11名であった。即ち、飛距離に関する正解率は69であった。一方、不正解は2名であり、有意差が認められない人が3名であった。なお有意差の有無は、有意水準10で判断した。打球方向性に関して正解であったのは、16名中、11名であった。即ち、打球方向性に関する正解率は69であった。一方、不正解は3名であり、有意差が認められない人が2名であった。なお有意差の有無は、有意水準10で判断した。振りやすさに関して正解であったのは、16名中、15名であった。即ち、打球方向性に関する正解率は94であった。一方、不正解は1名であり、有意差が認められない人が0名であった。なお有意差の有無は、有意水準10で判断した。これらの正解率は高い。本発明の優位性は明らかである。"