(57)【要約】      （修正有）【課題】インパクト時間決定の精度を向上させることができる、ゴルフスイングにおけるインパクト時間の抽出方法を提供する。【解決手段】３軸回りの角速度を計測可能なセンサが取り付けられたゴルフクラブ１でゴルフボールＢを打撃して得られる時系列波形であって、ｘ軸をゴルフクラブヘッド１ａのトウ‐ヒール方向に沿う方向に配向し、ｙ軸を打球方向に沿うように配向し、ｚ軸をシャフトの軸方向に一致するように配向したときの、ｘ軸回りの角速度ωｘ、ｙ軸回りの角速度ωｙ及びｚ軸回りの角速度ωｚの時系列波形の少なくとも１つの時系列波形を取得する工程と、ゴルフクラブヘッド１ａがボールＢとインパクトする直前に設けられたレーザ光が当該ゴルフクラブヘッド１ａ又はシャフトで遮光されるタイミングを取得する工程と、前記少なくとも１つの時系列波形と、前記タイミングとを用いてインパクト時間を決定する決定工程とを含んでいる。【選択図】図１
2014147554.tif 000002
【特許請求の範囲】
【請求項１】
３軸回りの角速度を計測可能なセンサが取り付けられたゴルフクラブでゴルフボールを打撃して得られる時系列波形であって、ｘ軸をゴルフクラブヘッドのトウ‐ヒール方向に沿う方向に配向し、ｙ軸を打球方向に沿うように配向し、ｚ軸をシャフトの軸方向に一致するように配向したときの、ｘ軸回りの角速度ωｘ、ｙ軸回りの角速度ωｙ及びｚ軸回りの角速度ωｚの時系列波形の少なくとも１つの時系列波形を取得する工程と、
ゴルフクラブヘッドがボールとインパクトする直前に設けられたレーザ光が当該ゴルフクラブヘッド又はシャフトで遮光されるタイミングを取得する工程と、
前記少なくとも１つの時系列波形と、前記タイミングとを用いてインパクト時間を決定する決定工程と
を含むことを特徴とする、ゴルフスイングにおけるインパクト時間の抽出方法。
【請求項２】
前記決定工程は、
前記レーザ光を受光する受光器から送信されるレーザ遮断の信号が前記センサで受信されたタイミングをトリガタイミングとすると、このトリガタイミングよりも前の所定時間帯における前記ｘ軸回りの角速度変化量を算出する算出工程、及び
算出された角速度変化量が所定の閾値を超えるか否かを判断する判断工程
を含む、請求項１に記載のゴルフスイングにおけるインパクト時間の抽出方法。
【請求項３】
前記判断工程において、算出された角速度変化量が所定の閾値を超えないと判断された場合に、前記決定工程は、
所定の時間幅において角速度ωｚがレンジオーバーしているか否かを判定する判定工程と、
前記判定工程においてレンジオーバーしていると判定された場合に第１の仮インパクト時間を設定し、レンジオーバーしていないと判定された場合に第２の仮インパクト時間を設定する設定工程と、
設定された第１又は第２の仮インパクト時間を基準として、角速度ωｘ、角速度ωｙ及び角速度ωｚのそれぞれの時系列波形からインパクト時間の３つの候補を取得する候補取得工程と、
３つのインパクト時間の候補から、所定の判断基準にしたがいインパクト時間を決定する決定工程と
を更に含む、請求項２に記載のゴルフスイングにおけるインパクト時間の抽出方法。
【請求項４】
前記第１の仮インパクト時間は、角速度ωｚの時系列波形においてレンジオーバーが発生した時点の時間である請求項３に記載のゴルフスイングにおけるインパクト時間の抽出方法。
【請求項５】
前記第２の仮インパクト時間は、角速度ωｚの時系列波形において最大値が現れる時間から第１の所定時間前までに最小値が現れた場合に、当該最小値が表れた時間から第２の所定時間前までに最大値が現れたときの時間である、請求項３又は４に記載のゴルフスイングにおけるインパクト時間の抽出方法。
【請求項６】
ボールの飛行方向を基準として、前記レーザ光が、ボールの後端より１ｍｍ〜１０ｍｍ後方に設けられている、請求項１〜５のいずれか一項に記載のゴルフスイングにおけるインパクト時間の抽出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はゴルフスイングにおけるインパクト時間の抽出方法に関する。
【背景技術】
【０００２】
ゴルファーにとってボールの飛距離を延ばすとともにボールを狙った方向・角度に飛ばすことは永遠のテーマである。そのためには、自身のスイングに合ったゴルフクラブを使用することが重要である。ゴルファーに合ったゴルフクラブを選定することは一般にフィッティングと称されている。
【０００３】
このフィッティングに関連して本出願人は、３軸回りの加速度及び角速度を計測可能なセンサがグリップに取り付けられたゴルフクラブで実際にボールを打撃してもらい、得られた計測結果からゴルファーのスイングを分析する方法を提案している（特願２０１１−１９６７２７号。以下、「第１先願発明」ともいう）。この方法により適切にスイングを分析することで、ゴルファーにマッチしたゴルフクラブを選定することができる。
【０００４】
第１先願発明では、ゴルファーのスイングを分類するために、センサを用いて計測したシャフト軸回りの角速度の時系列波形からインパクト時間を抽出している。具体的には、インパクトの瞬間を含むゴルフスイングの角速度の時系列波形において、シャフト軸（ｚ軸）回りの角速度が最大となる時間が仮インパクト時間として決定される。そして、より精度良くインパクト時間を抽出するために、仮インパクト時間の周辺（仮インパクト時間から所定時間前の時刻から、仮インパクトから所定時間後の時刻までの時間幅）において、ｘ軸回りの角速度ωｘが最小となる時間及びｚ軸回りの加速度が最小となる時間のうち早い方の時間をインパクト時間としている。
【０００５】
ところで、フィッティングを希望するゴルファーがボールを試打するに際し、ボールとクラブフェースが当たった場所（打点）によっては、図２３に示されるように、インパクト時ではなくインパクト後にシャフト軸回りの角速度ωｚが最大となることがある。図２３及び後出する図２４は、シャフト軸回りの角速度ωｚの時系列波形の一部（インパクト周辺）を示しており、同図において、横軸はデータ取得開始からのデータ番号、縦軸は角速度（ｄｅｇ／ｓ）を表している。図示された例におけるデータのサンプリング周波数は１０００Ｈｚであるので、データひとつ分（１点）が１ｍｓｅｃになる。図２３の例では、Ａ点が本来抽出すべきインパクト時間であるが、角速度ωｚが最大であるＢ点をインパクト時間であると誤認識する惧れがある。
【０００６】
また、計測に用いるジャイロセンサの測定可能範囲を上回る角速度でスイングした場合などでは、図２４に示されるように、センサにレンジオーバーが発生することがある。この場合も、本来抽出すべきインパクト時間はＣ点であるが、角速度ωｚがレンジの最大となったＤ点をインパクト時間であると誤認識する惧れがある。
【０００７】
そして、正しいインパクト時間が決定されないと、当該インパクト時間を用いて分析されたゴルファーのスイングに基づくフィッティングの精度に影響が出ることも考えられる。
【０００８】
これに対し、本出願人は、インパクト時間決定の精度を向上させることができる、ゴルフスイングにおけるインパクト時間の抽出方法を提案している（特願２０１２−０８１７７３号。以下、「第２先願発明」ともいう）。この第２先願発明は、シャフト軸回りの角速度ωｚがレンジオーバーしているかを判定し、その判定結果に応じてインパクト時間の決定方法を変更している。
【０００９】
第２先願発明では、レンジオーバーの有無に応じてそれぞれ所定の仮インパクト時間を設定し、この仮インパクト時間を基準として３つの軸回りの角速度それぞれの時系列波形からインパクト時間の３つの候補を取得し、これら３つの候補から所定の判断基準にしたがってインパクト時間を決定するので、前述したレンジオーバーなどの影響を排除して精度よくインパクト時間を決定することができる。  
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかし、センサにより計測した角速度の時系列波形だけからインパクト時間を抽出する方法では、稀にではあるが正しいインパクト時間を決定することができない場合がある。
例えば、インパクト時に、ヘッドの先端（トウ側）で試打された場合は、インパクトの直前にシャフト軸回りの角速度ωｚが最大となり、このときを仮インパクト時間と誤認識することがある。また、図２５に示されるように、アドレス時にゴルファーが急なワッグル動作を行った場合、インパクト時以外で計測データ内において角速度ωｚが最大となることがあり、この場合も仮インパクト時間を誤認識する惧れがある。さらに、センサの角速度計測レンジが従来品よりも広くなったためレンジオーバーすることがなくなり、図２６に示されるように、従来であればフラットになっていた時間帯に２つ以上のピークが発生すると、インパクト時間を正しく決定することができなくなる。すなわち、本来抽出すべき仮インパクト時間ｔ１とは異なるｔ２（図２５）又はｔ３（図２６）を仮インパクト時間と誤認識する惧れがある。
【００１１】
また、第２先願発明では、抽出した最大のシャフト軸回りの角速度ωｚの時刻を起点として、予め定めた時間の範囲内でインパクト時間を抽出しているため、将来、現状とは特性が大きく異なるゴルフボールやゴルフクラブが販売された場合、シャフトの振動特性が大きく変化し、予め定めた時間内に最小値又は変曲点が現れないことも考えられる。
【００１２】
本発明は、このような事情に鑑みてなされたものであり、インパクト時間決定の精度を向上させることができる、ゴルフスイングにおけるインパクト時間の抽出方法を提供することを目的としている。
【課題を解決するための手段】
【００１３】
（１）本発明のゴルフスイングにおけるインパクト時間の抽出方法（以下、単に「抽出方法」ともいう）は、３軸回りの角速度を計測可能なセンサが取り付けられたゴルフクラブでゴルフボールを打撃して得られる時系列波形であって、ｘ軸をゴルフクラブヘッドのトウ‐ヒール方向に沿う方向に配向し、ｙ軸を打球方向に沿うように配向し、ｚ軸をシャフトの軸方向に一致するように配向したときの、ｘ軸回りの角速度ωｘ、ｙ軸回りの角速度ωｙ及びｚ軸回りの角速度ωｚの時系列波形の少なくとも１つの時系列波形を取得する工程と、
ゴルフクラブヘッドがボールとインパクトする直前に設けられたレーザ光が当該ゴルフクラブヘッド又はシャフトで遮光されるタイミングを取得する工程と、
前記少なくとも１つの時系列波形と、前記タイミングとを用いてインパクト時間を決定する決定工程と
を含むことを特徴としている。
【００１４】
本発明の抽出方法では、ゴルフクラブに取り付けられたセンサから得られるｘ軸回りの角速度ωｘ、ｙ軸回りの角速度ωｙ及びｚ軸回りの角速度ωｚの時系列波形の少なくとも１つの時系列波形と、ゴルフクラブヘッドがボールとインパクトする直前に設けられたレーザ光が遮光されるタイミングとを用いてインパクト時間を決定するので、インパクト時間決定の精度を向上させることができる。
一般的なヘッドスピードでボールを打球した場合、ボールとクラブとが接触してから５／１００００秒までの間にクラブからボールにエネルギーが伝達され、その反動でボールからクラブに振動が伝わる。本明細書では、かかる振動がクラブに伝達される瞬間を「インパクト時間」としている。
【００１５】
（２）前記（１）の抽出方法において、前記決定工程は、
前記レーザ光を受光する受光器から送信されるレーザ遮断の信号が前記センサで受信されたタイミングをトリガタイミングとすると、このトリガタイミングよりも前の所定時間帯における前記ｘ軸回りの角速度変化量を算出する算出工程、及び
算出された角速度変化量が所定の閾値を超えるか否かを判断する判断工程
を含むものとすることができる。
【００１６】
（３）前記（２）の抽出方法において、前記判断工程において、算出された角速度変化量が所定の閾値を超えないと判断された場合に、前記決定工程は、
所定の時間幅において角速度ωｚがレンジオーバーしているか否かを判定する判定工程と、
前記判定工程においてレンジオーバーしていると判定された場合に第１の仮インパクト時間を設定し、レンジオーバーしていないと判定された場合に第２の仮インパクト時間を設定する設定工程と、
設定された第１又は第２の仮インパクト時間を基準として、角速度ωｘ、角速度ωｙ及び角速度ωｚのそれぞれの時系列波形からインパクト時間の３つの候補を取得する候補取得工程と、
３つのインパクト時間の候補から、所定の判断基準にしたがいインパクト時間を決定する決定工程と
を更に含むものとすることができる。
【００１７】
（４）前記（３）の抽出方法において、前記第１の仮インパクト時間を、角速度ωｚの時系列波形においてレンジオーバーが発生した時点の時間とすることができる。
【００１８】
（５）前記（３）又は（４）の抽出方法において、前記第２の仮インパクト時間を、角速度ωｚの時系列波形において最大値が現れる時間から第１の所定時間前までに最小値が現れた場合に、当該最小値が表れた時間から第２の所定時間前までに最大値が現れたときの時間とすることができる。
【００１９】
また、前記（３）の抽出方法において、角速度ωｚの時系列波形において最大値が現れる時間から第１の所定時間前までに最小値が現れない場合に、当該最大値を第２の仮インパクト時間とすることができる。
また、前記第１の所定時間は１０〜３０ｍｓｅｃであり、第２の所定時間は２０〜４０ｍｓｅｃであってもよい。
【００２０】
また、角速度ωｘの時系列波形から取得されるインパクト時間の候補は、第１の仮インパクト時間前後の第１の所定時間幅において最初の極大値が現れる時間から第３の所定時間前までに最小値が現れる時間とし、
角速度ωｙの時系列波形から取得されるインパクト時間の候補は、第１の仮インパクト時間前後の第２の所定時間幅において最初に変極点が現れる時間とし、
角速度ωｚの時系列波形から取得されるインパクト時間の候補は、第１の仮インパクトから第４の所定時間前までに最初に変極点が現れる時間と
することができる。
【００２１】
前記第３の所定時間は２０〜４０ｍｓｅｃであり、第４の所定時間は３０〜５０ｍｓｅｃであってもよい。
また、角速度ωｘの時系列波形から取得されるインパクト時間の候補は、第２の仮インパクト時間前後の第１の所定時間幅において最初の極大値が現れる時間から第３の所定時間前までに最小値が現れる時間とし、
角速度ωｙの時系列波形から取得されるインパクト時間の候補は、第２の仮インパクト時間前後の第２の所定時間幅において最初に変極点が現れる時間とし、
角速度ωｚの時系列波形から取得されるインパクト時間の候補は、第２の仮インパクト前後の第３の所定時間幅において最初に変極点が現れる時間と
することができる。
【００２２】
また、前記第１の所定時間幅は、第１又は第２の仮インパクト時間より５０ｍｓｅｃ前から、第１又は第２の仮インパクト時間より２０ｍｓｅｃ後までの時間幅であり、
前記第２の所定時間幅は、第１又は第２の仮インパクト時間より３０ｍｓｅｃ前から、第１又は第２の仮インパクト時間より５０ｍｓｅｃ後までの時間幅であってもよい。
また、前記第３の所定時間幅は、第２の仮インパクト時間より４０ｍｓｅｃ前から、第２の仮インパクト時間より１５ｍｓｅｃ後までの時間幅であってもよい。
【００２３】
（６）前記（１）〜（５）の抽出方法において、ボールの飛行方向を基準として、前記レーザ光が、ボールの後端より１ｍｍ〜１０ｍｍ後方に設けられていてもよい。
【発明の効果】
【００２４】
本発明の抽出方法によれば、インパクト時間決定の精度を向上させることができる。
【図面の簡単な説明】
【００２５】
【図１】本発明におけるスイング特徴量を計測する方法を説明する図である。
【図２】センサが取り付けられたゴルフクラブの一部拡大斜視図である。
【図３】（ａ）は図２に示されるセンサの平面図であり、（ｂ）は同センサの側面図である。
【図４】レーザ発振器およびレーザ受光器の配置を説明する図である。
【図５】本発明の抽出方法の一実施の形態のフローチャートである。
【図６】角速度と角速度変化量の例を示す図である。
【図７】レンジオーバーしているｚ軸回りの角速度の時系列波形の例を示す図である。
【図８】ｘ軸回りの角速度の時系列波形の例を示す図である。
【図９】ｙ軸回りの角速度の時系列波形の例を示す図である。
【図１０】変極点の求め方を示す図である。
【図１１】レンジオーバーしていないｚ軸回りの角速度の時系列波形の例を示す図である。
【図１２】スイング中のシャフトの撓みの挙動を説明する図である。
【図１３】適用例に係るフィッティング方法において曲げ剛性が測定されるシャフトの４つの位置を説明する図である。
【図１４】曲げ剛性の測定方法を説明する図である。
【図１５】スイングにおけるアドレス及びテイクバックを示す図である。
【図１６】スイングにおけるトップ及びダウンスイングを示す図である。
【図１７】スイングにおけるダウンスイング及びインパクトを示す図である。
【図１８】スイングにおけるフォロースルー及びフィニッシュを示す図である。
【図１９】スイングにおけるコック方向の角速度の時間経過にしたがう変化を示す図である。
【図２０】振り易さの評価基準を説明する図である。
【図２１】「６インチ」の測定点におけるスイング特徴量（４）と曲げ剛性（ＥＩ値）との関係を示す図である。
【図２２】「６インチ」の測定点におけるＥＩ値算出のフローを示す図である。
【図２３】シャフト軸回りの角速度の時系列波形の一例である。
【図２４】シャフト軸回りの角速度の時系列波形の他の例である。
【図２５】インパクト時以外でωｚが最大となる時系列波形の一例である。
【図２６】計測レンジの拡大によりレンジオーバーがなくなった時系列波形の一例である。
【発明を実施するための形態】
【００２６】
以下、添付図面を参照しつつ、本発明の抽出方法の実施の形態を詳細に説明する。
本発明の抽出方法は、フィッティングを行うためにゴルファーのスイングの特徴を分析する際に用いられるインパクト時間を精度良く抽出する方法である。ゴルファーに特有のスイングの特徴は、後述するように定量的に把握することができ、本明細書では、かかる特徴を「スイング特徴量」と称する。
【００２７】
〔角速度の取得〕
スイング特徴量は、図１に示されるように、ゴルフクラブのフィッティングを希望するゴルファーに実際にスイングをしてもらい、そのスイングから計測することができる。ゴルフクラブ１のグリップエンドには、図２〜３に示されるように、３軸回りの角速度を計測することができるセンサ２がアダプター３を介して取り付けられている。本実施の形態におけるセンサ２は、平面視正方形状の箱体からなるケーシング２ａを備えており、このケーシング２ａは、両面テープ、接着剤、ねじ止めなどによりグリップエンドに固定することができる。図１に示される例では、ゴルファーＧは右利きであり、所定位置にセットされたボールＢを打撃するためのスイングを開始する直前のアドレス状態である。
【００２８】
センサ２は無線式であり、測定されたデータは、無線により、データ解析装置としてのコンピュータ１０に内蔵された無線受信装置（図示せず）に送信される。無線通信として、例えばＢｌｕｅｔｏｏｔｈ（ブルートゥース（Ｂｌｕｅｔｏｏｔｈは登録商標である。））の規格及び技術を用いることができる。
【００２９】
センサ２は、３軸方向（ｘ軸方向、ｙ軸方向、及びｚ軸方向）回りの角速度を計測することができる角速度センサ（図示せず）を内蔵している。センサ２は、更に、Ａ／Ｄ変換器、ＣＰＵ、無線インターフェース、無線アンテナ及び電源を備えている。電源としては、例えばボタン型のリチウムイオン電池などを用いることができる。電池は充電可能なものであってもよい。そして、センサ２は、電池を充電するための充電回路を備えていてもよい。使用可能なセンサ２の例としては、ワイヤレステクノロジー社製のＷＡＡ−０１０（商品名）をあげることができる。
【００３０】
なお、センサ２からの信号を受信する前記無線受信装置は、無線アンテナ、無線インターフェース、ＣＰＵ及びネットワークインターフェースを備えている。
【００３１】
データ解析装置としてのコンピュータ１０は、キーボード４及びマウス５からなる入力部６と、表示部７とを備えている。また、図示していないが、コンピュータ１０は、ハードディスク、メモリ、ＣＰＵ及びネットワークインターフェースを備えている。
【００３２】
センサ２は、ｘ軸、ｙ軸及びｚ軸の各軸回りの角速度を検知する。これらの角速度はアナログ信号として得られ、このアナログ信号は、センサ２に内蔵されているＡ／Ｄ変換器によってデジタル信号に変換される。Ａ／Ｄ変換器からの出力はＣＰＵに伝達されて１次フィルタリングなどの演算処理が実行される。こうしてセンサ２内で処理されたデータは、無線インターフェースを介して、無線アンテナからコンピュータ１０に内蔵された無線受信装置に送信される。
【００３３】
センサ２から送信されたデータは、無線受信装置側の無線アンテナを介して、無線インターフェースによって受信される。受信されたデータは、コンピュータ１０のＣＰＵで演算処理される。
【００３４】
コンピュータ１０に送られたデータは、ハードディスクなどのメモリ資源に記録される。ハードディスクには、データ処理などに必要なプログラム及びデータなどが記憶されている。このプログラムは、ＣＰＵに、必要なデータ処理を実行させる。ＣＰＵは、各種の演算処理を実行可能であり、演算結果は、表示部７又は図示しない印刷装置などによって出力される。
【００３５】
センサ２のグリップエンドへの取付に際しては、計測軸とゴルフクラブ１との関係が考慮される。本実施の形態では、センサ２のｚ軸がゴルフクラブ１のシャフト軸に一致している。センサ２のｘ軸は、ゴルフクラブ１のヘッド１ａのトウーヒール方向に沿うように配向される。また、センサ２のｙ軸は、ヘッド１ａのフェース面の法線方向に沿うように配向される。このようにセンサ２を取り付けることで、演算を簡略化することができる。
【００３６】
本実施の形態では、局座標系が考慮され、この局座標系のｘ軸、ｙ軸及びｚ軸は、三次元直交座標系である。本実施の形態では、ｚ軸が前記ゴルフクラブ１のシャフト軸とされ、ｘ軸は、ヘッド１ａのトウーヒール方向に沿うように配向される。また、ｙ軸は、打球方向に沿うように配向される。
【００３７】
すなわち、局座標系のｚ軸は、センサ２のｚ軸に一致しており、局座標系のｙ軸は、センサ２のｙ軸に一致している。また、局座標系のｘ軸は、センサ２のｘ軸に一致している。
【００３８】
センサ２によって、時系列的に連続した複数のデータを得ることができる。単位時間当たりのデータ数は、サンプリング周波数に依存する。サンプリング周波数は、例えば１０００Ｈｚとすることができる。
【００３９】
〔レーザ光遮光タイミングの取得〕
本実施の形態では、センサ２により取得することができる３軸回りの角速度の時系列波形と、ゴルフクラブヘッド１ａがボールＢとインパクトする直前に設けられたレーザ光が遮光されるタイミングとを用いてインパクト時間を決定している。このため、本実施の形態では、図４に示されるように、ボールＢの飛行方向（図４の矢印Ｆ参照）を基準として当該ボールＢの後端Ｂｅから後方に距離Ｌだけ離れた位置にレーザ光が設けられるように、通常のレーザ発振器３０およびレーザ受光器３１が配設されている。距離Ｌは、ゴルファーのヘッドスピードなどを考慮して選定することができ、例えば１〜１０ｍｍ程度とすることができる。なお、図４では、わかり易くするために前記距離Ｌの大きさを誇張して描いている。
【００４０】
〔インパクト時間の抽出〕
図５は、本発明の一実施の形態に係るインパクト時間の抽出方法のフローチャートである。
まず、ステップＳ１において、トリガタイミングを０として、−２５０〜−５０データの間でｘ軸回りの角速度から角速度変化量が算出される。ここで、「トリガタイミング」とは、クラブヘッドによってレーザ光が遮断ないし遮光されたタイミングではなく、レーザ受光器３１から送信されるレーザ遮断の信号がコンピュータ１０を経由してセンサ２にて受信されたタイミングのことである。本実施の形態では、実際のレーザ遮断から「トリガタイミング」までの間に、クラブヘッドがボールを打撃するインパクトタイミングが存在するという前提で前記データ数などを設定している。また、本実施の形態におけるデータのサンプリング周期は１／１０００秒である。したがって、トリガタイミングを０として、−２５０〜−５０データの間とは、トリガタイミングの２５０／１０００秒前から５０／１０００秒前までの間のことであり、特許請求の範囲における「トリガタイミングよりも前の所定時間帯」の一例である。なお、「角速度変化量」とは、１データ毎の角速度の差分を示している。
【００４１】
図６の（ａ）は、或るテスターがヘッドスピード５０〜５５ｍ／ｓで試打をしたときの角速度を示しており、図６の（ｂ）は、図６の（ａ）に示される角速度の角速度変化量を示している。図６の角速度は、前述したトリガタイミングを０として、−２５０〜−５０データの間でのｘ軸回りの角速度である。
【００４２】
ついで、ステップＳ２において、ステップＳ１で算出された角速度変化量において３０００ｄｅｇ／ｓを超える点があるか否かを判断する。本実施の形態では、角速度変化量が３０００ｄｅｇ／ｓを超える点を変極点と判断している。３０００ｄｅｇ／ｓを超える点があると判断される（Ｙｅｓ）と、ステップＳ３に進み、当該ステップＳ３において、前記データ範囲（−２５０〜−５０データ）において最初に３０００ｄｅｇ／ｓを超えた点がＳ点（変極点）とされる。
【００４３】
ついで、ステップＳ４において、ステップＳ３において設定したＳ点が−２４５〜−５０データの間に入っているか否かが判断される。そして、Ｓ点が−２４５〜−５０データの間に入っていると判断されると、ステップＳ５に進み、当該ステップＳ５において、Ｓ点から５データ前〜Ｓ点の間で角速度変化量が２０００ｄｅｇ／ｓを超える点をインパクトタイミングとする。
【００４４】
一方、ステップＳ２において、角速度変化量が３０００ｄｅｇ／ｓを超える点が存在しない（Ｎｏ）と判断された場合、およびステップＳ４において、Ｓ点が−２４５〜−５０データの間に入っていない（Ｎｏ）と判断された場合は、ステップＳ６に進み、当該ステップＳ６において、シャフト軸（ｚ軸）回りの角速度ωｚの時系列波形において最大値が抽出される。すなわち、前述した第２先願発明にしたがって、シャフト軸回りの角速度ωｚがレンジオーバーしているかを判定し、その判定結果に応じてインパクト時間の決定方法を変更している。ステップＳ２において変極点があると判断されても、当該変極点が−２４５〜−５０データの間に入っていない場合は、この領域以外にインパクトが存在する可能性があるので、変極点はないものとしている。
【００４５】
ついでステップＳ７において、ステップＳ６において抽出された最大値が所定時間（例えば、３ｍｓｅｃ以上）連続しているか否かが判断される。最大値が所定時間連続することは、図７に示されるように、手首の回転速度が速いなどの理由によりセンサがレンジオーバーしているものと考えることができる。なお、図７に示される例では、横軸を時間としているが、図２３〜２４に示される例のように時系列波形の横軸がデータ取得開始からのデータ番号である場合、所定周波数のサンプリングにより取得されるデータが何点連続するかでレンジオーバーしているか否かを判断することもできる。
【００４６】
ステップＳ７において、抽出された最大値が所定時間連続していると判断されると、ステップＳ８に進み、当該ステップＳ８において第１の仮インパクト時間が決定される。レンジオーバーしている場合は最大値が連続して続くことから、最初に最大値が現れるタイミング（時間）を第１の仮インパクト時間Ｅ１とする。
【００４７】
ついで、ステップＳ９において、ステップＳ８で決定された第１の仮インパクト時間Ｅ１を用いて、ｘ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ａ１を抽出する。このインパクト時間ａ１の抽出は、次のようにして行う。
【００４８】
図８は、ｘ軸回りの角速度ωｘの時系列波形の一例の一部（インパクト周辺）を示している。ステップＳ８において、第１の仮インパクト時間Ｅ１が決定されると、ｘ軸回りの角速度ωｘの時系列波形において当該第１の仮インパクト時間Ｅ１に対応する時間を基準として、インパクト時間ａ１が抽出される。具体的に、まず第１の仮インパクト時間Ｅ１前後の所定の時間幅（第１の所定時間幅）において最初の極大値が現れる時間を求める。所定の時間幅（第１の所定時間幅）としては、例えば、第１の仮インパクト時間Ｅ１より５０ｍｓｅｃ前から、第１の仮インパクト時間Ｅ１より２０ｍｓｅｃ後までの時間幅とすることができる。図８では、ｔ１が最初の極大値が現れる時間を示している。
【００４９】
ついで、最初の極大値が現れる時間ｔ１から、所定の時間（第３の所定時間）前までに最小値が現れる時間をインパクト時間ａ１とする。所定の時間（第３の所定時間）としては、例えば、３ｍｓｅｃとすることができる。
【００５０】
ついで、同じくステップＳ９において、ステップＳ８で決定された第１の仮インパクト時間Ｅ１を用いて、ｙ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ｂ１を抽出する。このインパクト時間ｂ１の抽出は、次のようにして行う。
【００５１】
図９は、ｙ軸回りの角速度ωｙの時系列波形の一例の一部（インパクト周辺）を示している。ステップＳ８において、第１の仮インパクト時間Ｅ１が決定されると、ｙ軸回りの角速度ωｙの時系列波形において当該第１の仮インパクト時間Ｅ１に対応する時間を基準として、インパクト時間ｂ１が抽出される。具体的に、第１の仮インパクト時間Ｅ１前後の所定の時間幅（第２の所定時間幅）において最初に変極点が現れる時間をインパクト時間ｂ１とする。所定の時間幅（第２の所定時間幅）としては、例えば、第１の仮インパクト時間Ｅ１より３０ｍｓｅｃ前から、第１の仮インパクト時間Ｅ１より５０ｍｓｅｃ後までの時間幅とすることができる。
【００５２】
変極点は、例えば次のようにして求めることができる。図１０は、変極点の求め方を示す図である。連続する３点で、真ん中の点を中心として、この中心点と前後の点とのなす角度θを求める。図１０における３つの連続する点、（ｉ−２）、（ｉ−１）、（ｉ）は、ほぼ直線上にあることから、点（ｉ−２）と点（ｉ−１）を結ぶ線分と、点（ｉ−１）と点（ｉ）とを結ぶ線分とのなす角度θ≒πである。そして、真ん中に来る点を１点ずつずらして角度θを算出し続ける。変極点以外の点では、θ≒πとなるが、変極点にきたときの角度θは小さくなる。θが最小になったときの点を「変極点」とする。図１０の例では、点（ｉ）が変極点である。
【００５３】
ついで、同じくステップＳ９において、ステップＳ８で決定された第１の仮インパクト時間Ｅ１を用いて、ｚ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ｃ１を抽出する。このインパクト時間ｃ１の抽出は、次のようにして行う。
【００５４】
この抽出の前提となっているｚ軸回りの角速度ωｚは、図７に示されるように、レンジオーバーしている。この場合、最大値が最初に現れる時間、すなわちレンジオーバーが発生した時点の時間を第１の仮インパクト時間Ｅ１としているが、この第１の仮インパクト時間Ｅ１から所定の時間（第４の所定時間）前までに最初に変極点が現れる時間をインパクトｃ１とする。所定の時間（第４の所定時間）としては、例えば、４０ｍｓｅｃとすることができる。なお、変極点は、ステップＳ９で説明した方法を用いて取得することができる。
【００５５】
一方、ステップＳ７において、最大値が所定時間連続していないと判断された場合は、ステップＳ１０に進み、第２の仮インパクト時間Ｅ２が決定される。この第２の仮インパクト時間Ｅ２は、次のようにして決定される。図１１は、レンジオーバーしていないｚ軸回りの角速度の時系列波形の例を示す図であり、同図において、ｍ１が最大値である。この最大値ｍ１が現れる時間ｔ３から、所定時間（第１の所定時間）前までに最小値が現れた場合、この最小値が現れた時間をｔ４とし、ついでこのｔ４から所定時間（第２の所定時間）前までに最大値が現れたとき、この最大値が現れた時間を第２の仮インパクト時間Ｅ２とする。所定時間（第１の所定時間）としては、例えば２０ｍｓｅｃとし、所定時間（第２の所定時間）としては、例えば３０ｍｓｅｃとすることができる。
【００５６】
ついで、ステップＳ１１において、ステップＳ１０で決定された第２の仮インパクト時間Ｅ２を用いて、ｘ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ａ２を抽出する。このインパクト時間ａ２の抽出は、前述したインパクト時間ａ１と同様にして行うことができる。
【００５７】
ついで、同じくステップＳ１１において、ステップＳ１０で決定された第２の仮インパクト時間Ｅ２を用いて、ｙ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ｂ２を抽出する。このインパクト時間ｂ２の抽出は、前述したインパクト時間ｂ１と同様にして行うことができる。
【００５８】
ついで、同じくステップＳ１１において、ステップＳ１０で決定された第２の仮インパクト時間Ｅ２を用いて、ｚ軸回りの時系列波形からインパクト時間の候補の１つであるインパクト時間ｃ２を抽出する。このインパクト時間ｃ２の抽出は、前後の所定の時間幅（第３の所定時間幅）において最初に変極点が現れる時間とする。所定の時間幅（第３の所定時間幅）としては、例えば、第２の仮インパクト時間Ｅ２より４０ｍｓｅｃ前から、第２の仮インパクト時間Ｅ２より１５ｍｓｅｃ後までの時間幅とすることができる。
【００５９】
ついで、ステップＳ１２において、ステップＳ９又はステップＳ１１によって抽出されたインパクト時間の候補ａ１、ｂ１及びｃ１、又は、インパクト時間の候補ａ２、ｂ２及びｃ２の各組について、抽出に失敗しているインパクト時間のリジェクトを行う。センサにより計測された各軸の時系列波形には、レンジオーバーしているものが含まれていることがあり、かかるレンジオーバーなどにより、１つの時系列波形において抽出ミスが発生した場合、その抽出ミスがインパクト時間の真値に影響を及ぼす可能性がある。そこで、本実施の形態では、このような影響を回避するために、抽出されたインパクト時間の候補の中から抽出に失敗していると考えられるものをリジェクトしている。
【００６０】
具体的に、ａ１、ｂ１及びｃ１（ａ２、ｂ２及びｃ２についても同様である）の３つの点の中で、最も速く発生している点、及び最も遅く発生している点を決定する。ついで、最も速く（遅く）発生している点が、他の２点に比べて、例えば５ｍｓｅｃ以上離れていれば、その点はインパクト時間の抽出に失敗していると判断し、リジェクトする。
【００６１】
ついで、ステップＳ１３において、ステップＳ１２で抽出に失敗していないとされた点の中で、最も早いタイミングで発生している点をインパクト時間とする。
【００６２】
〔本発明の抽出方法の適用例〕
本発明の抽出方法は、前述したように、フィッティングにおいてインパクト時間を決定する際に適用することができる。本発明の抽出方法は、センサで計測したシャフト軸回りの角速度からインパクト時間を抽出する方法であれば、特に制限されることなく適用が可能であるが、例えば、ゴルファーのスイング全体に着目する以下のフィッティング方法に好適に用いることができる。
【００６３】
このスイング全体に着目するフィッティング方法は、ゴルフクラブのシャフトの曲げが、トップからインパクトに向けてスイングが進行するにしたがい当該シャフトの手元側から先端側に伝わることに着目している。そして、或るゴルファーのトップからインパクトまでの時間に伴うスイング特徴（スイング特徴の内容については後述する）と、当該ゴルファーにマッチしたシャフトのインチ毎の硬さとの間には相関関係があるという仮定のもと、鋭意研究・検討を重ねた結果、完成されたものである。
【００６４】
すなわち、ボールを打撃するときのゴルファーのスイングは、アドレス→トップ→インパクトと推移するが、その際、ゴルフクラブのシャフトは、当該ゴルフクラブの先端に比較的重量が大きいヘッドが存在するため、その慣性により曲げが生じる。この曲げは、スイングの全過程において、シャフトの同一箇所に生じるのではなく、図１２に示されるように、トップからインパクトに向けてシャフトの手元側から先端側に伝わる。換言すれば、トップからインパクトに向けてスイングが進行するにしたがい、シャフトにおける曲げの位置が当該シャフトの手元側から先端側に移動する。
【００６５】
具体的に、アドレスからテイクバックを行い、トップに至った時点（図１２において（１）で示される時点）では、シャフトの手元付近に曲げが生じる。ついで、切り返しを行い、ダウンスイング初期（図１２において（２）で示される時点）に至ると、曲げはシャフトの先端側にやや移動する。さらに、ゴルファーの腕が水平になる時点（図１２において（３）で示される時点）では、曲げはシャフト中央よりも先端側に移動する。そして、インパクト直前（図１２において（４）で示される時点）では、曲げはシャフトの先端付近まで移動する。
【００６６】
このようにスイング時におけるシャフトの曲げが、トップからインパクトに向けてシャフトの手元側から先端側に伝わることに鑑み、前記時間帯（１）〜（４）におけるゴルファーのスイング特徴に着目することで、シャフトのインチ毎に最適な曲げ剛性を選定することができる。
【００６７】
具体的に、図１３に示されるように、シャフト２０を４つの領域に分割して、各領域中の１点の曲げ剛性を定義した。本適用例では、シャフト２０のチップ端２０ａから３６インチの箇所を測定点（１）とし、２６インチの箇所を測定点（２）とし、１６インチの箇所を測定点（３）とし、６インチの箇所を測定点（４）としている。そして、シャフト２０の４つの測定点における曲げ剛性を計測し、数値化している。なお、本明細書において、「インチ毎の」とは、１インチの、２インチの、・・・という意味ではなく、「シャフトの一端から所定インチの距離にある複数の箇所について、各箇所の」という意味であり、本適用例における「所定インチ」とは、前述したようにシャフト２０のチップ端２０ａから３６インチ、２６インチ、１６インチ、６インチである。
【００６８】
シャフトのインチ毎の曲げ剛性（ＥＩ値：Ｎ・ｍ２）は、例えばインテスコ社製の２０２０型計測機（最大荷重５００ｋｇｆ）を用いて図１４に示されるようにして測定することができる。
【００６９】
具体的に、２つの支持点１１、１２においてシャフト２０を下方から支えつつ、測定点Ｐに上方から荷重Ｆを加えたときのたわみ量αを測定する。測定点Ｐは、本適用例の場合、シャフト２０のチップ端２０ａから３６インチ、２６インチ、１６インチ及び６インチの４箇所である。支持点１１と支持点１２との間の距離（スパン）は２００ｍｍである。また、測定点Ｐは、支持点１１と支持点１２の中間点である。上方から荷重Ｆを加える圧子１３の先端は、シャフト２０を傷付けないように丸められている。圧子１３の先端の断面形状は、シャフト軸方向に平行な断面において、１０ｍｍの曲率半径を有している。シャフト軸方向に対して垂直な方向の断面において、圧子１３の先端の断面形状は直線であり、その長さは４５ｍｍである。
【００７０】
支持体１４は、支持点１１においてシャフト２０を下方から支持する。支持体１４の先端は、凸状の丸みを有している。支持体１４の先端の断面形状は、シャフト軸方向に対して平行な断面において１５ｍｍの曲率半径を有する。シャフト軸方向に対して垂直な断面において、支持体１４の先端の断面形状は直線であり、その長さは５０ｍｍである。支持体１５の形状は支持体１４と同一である。支持体１５は、支持点１２においてシャフト２０を下方から支持する。支持体１５の先端は、凸状の丸みを有している。支持体１５の先端の断面形状は、シャフト軸方向に対して平行な断面において１５ｍｍの曲率半径を有する。シャフト軸方向に対して垂直な断面において、支持体１５の先端の断面形状は直線であり、その長さは５０ｍｍである。
【００７１】
前述した支持体１４及び支持体１５を固定した状態で、圧子１３を５ｍｍ／ｍｉｎの速度で下方へ移動させる。そして、荷重Ｆが２０ｋｇに達した時点で圧子１３の移動を終了させる。圧子１３の移動を終了させた瞬間におけるシャフト２０のたわみ量α（ｍｍ）を測定し、以下の式（１）にしたがって曲げ剛性ＥＩ（Ｎ・ｍ２）を計算する。
曲げ剛性ＥＩ（Ｎ・ｍ２）＝３２．７／α    ・・・・・・（１）
【００７２】
そして、測定されたシャフトのインチ毎の曲げ剛性を指標としてフィッティングが行われる。インチ毎のシャフトの硬さ（曲げ剛性）は、トップからインパクトまでの時間に伴うスイング特徴と相関関係があり、個々のゴルファーにおける前記スイング特徴が分かれば、その特徴に合ったインチ毎のシャフトの硬さを決めることができる。スイング中のシャフトの曲げないし変形（たわみ量）は、前述したように、トップからダウンスイングにかけてシャフトの手元側から先端側に曲げが伝わっていく。本適用例では、この曲げの伝わりに着目し、トップ付近のシャフトの撓み量がトップ付近のグリップの挙動（角速度）の速い遅いに関係があることから、この速度が速いゴルファーには硬めのシャフトを提供し、遅いゴルファーには軟らかめのシャフトを提供する。
【００７３】
本適用例では、前述した方法によりインチ毎に測定されたシャフトの曲げ剛性に応じて複数段階のランク値のうちいずれかの値を付与している。具体的に、曲げ剛性に応じて１０段階のＩＦＣのうちいずれかの値を付与している。このＩＦＣは、Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｆｌｅｘ  Ｃｏｒｄ（インターナショナル・フレックス・コード）の略であり、シャフトの硬さを表すものとして本出願人により提案されている指標である。
【００７４】
以下の表１〜４は、それぞれ測定点（１）〜（４）におけるシャフトのＥＩ値からＩＦＣへの変換表である。シャフトの硬さを１０段階に分ける方法として、市販されている全てのシャフトを対象にして１０段階に分ける方法や、使用頻度などを考慮して、フィッターがユーザーに提供しようとして用意したシャフトの範囲内で１０段階に分ける方法など、いくつかの方法が考えられるが、本適用例では、後者の方法によりフィッティングを行っている。
【００７５】
【表１】
2014147554.tif 000003
【００７６】
【表２】
2014147554.tif 000004
【００７７】
【表３】
2014147554.tif 000005
【００７８】
【表４】
2014147554.tif 000006
【００７９】
<スイング特徴量>
図１５〜１８は、ゴルファーによるアドレスからフィニッシュまでのスイングを説明する図である。スイングには、インパクト後のフォールスルーが含まれるが、本適用例では、アドレスからインパクトまでのスイングの特徴に着目している。
【００８０】
図１５〜１８は、ゴルファーを正面から見た図である。スイングの始まりはアドレスであり、スイングの終わりはフィニッシュと呼ばれている。スイングは、（Ｓ１）、（Ｓ２）、（Ｓ３）、（Ｓ４）、（Ｓ５）、（Ｓ６）、（Ｓ７）、（Ｓ８）の順で進行する。図１５の（Ｓ１）がアドレスであり、（Ｓ２）がテイクバックである。図１６の（Ｓ３）がトップ（トップオブスイング）である。通常、トップでは、スイング中におけるヘッドの移動速度が最小である。図１６の（Ｓ４）はダウンスイングである。図１７の（Ｓ５）もダウンスイングであるが、図１６の（Ｓ４）よりもダウンスイングが進行した状態である。図１７の（Ｓ６）はインパクトであり、ゴルフクラブ１のヘッド１ａとボールＢとが衝突した瞬間である。図１８の（Ｓ７）はフォロースルーであり、（Ｓ８）はフィニッシュである。フィニッシュで、スイングが終了する。
【００８１】
本適用例では、前述したスイングにおける種々の段階のうちトップ付近からインパクトに至るダウンスイング中のコック方向の角速度ωｙに着目し、時間経過にしたがって当該角速度ωｙを細分化して定量化する。なお、本明細書において「トップ付近」とは、トップ直前の所定時間及びトップ直後の所定時間を含む時間帯を意味しており、具体的には、例えばトップ−５０ｍｓから、トップ＋５０ｍｓまでの１００ｍｓの時間帯を意味している。
【００８２】
図１９は、或るスイングにおけるアドレスからインパクトまでの時間（ｓ）とコック方向の角速度ωｙ（ｄｅｇ／ｓ）との関係を示している。本適用例では、図１９に示されるように、スイングの時間経過にしたがって４つのスイング特徴（１）〜（４）を設定し、各スイング特徴を定量化している。
【００８３】
スイング特徴量（１）は、トップ付近のコック方向の角速度ωｙの傾きであり、例えばトップから５０ｍｓ前の角速度ωｙと、トップから５０ｍｓ後の角速度ωｙとの和で表すことができる。このスイング特徴量（１）は、前述したシャフトのチップ端から３６インチの測定点における曲げ剛性と相関関係がある。
【００８４】
スイング特徴量（２）は、トップから、角速度ωｙが最大となる時点までの当該角速度ωｙの平均値である。トップからインパクトまでの角速度ωｙにおける最大値を求め、トップから、この最大値となる時点までの角速度ωｙの累積値を、トップから、前記最大値となる時点までの時間で除することにより前記平均値を求めることができる。このスイング特徴量（２）は、前述したシャフトのチップ端から２６インチの測定点における曲げ剛性と相関関係がある。
【００８５】
スイング特徴量（３）は、角速度ωｙが最大となる時点からインパクトまでの当該角速度ωｙの平均値であり、前記最大値となる時点からインパクトまでの角速度ωｙの累積値を、前記最大値となる時点からインパクトまでの時間で除することにより前記平均値を求めることができる。このスイング特徴量（３）は、前述したシャフトのチップ端から１６インチの測定点における曲げ剛性と相関関係がある。
【００８６】
スイング特徴量（４）は、トップからインパクトまでの角速度ωｙの平均値であり、トップからインパクトまでの角速度ωｙの累積値を、トップからインパクトまでの時間で除することにより前記平均値を求めることができる。このスイング特徴量（４）は、前述したシャフトのチップ端から６インチの測定点における曲げ剛性と相関関係がある。
【００８７】
以上のスイング特徴量（１）〜（４）は、フィッティングを希望するゴルファーに所定数のボール、例えば５球のボールを試打してもらい、各打球時に算出したスイング特徴量の平均を当該ゴルファーのスイング特徴量と設定することができる。
【００８８】
<インチ毎のシャフト剛性の算出>
ついで、算出されたスイング特徴量（１）〜（４）に基づいて、当該ゴルファーに適したインチ毎のシャフト剛性を計算する。この計算に先立って、予め各スイング特徴量について、スイング特徴量と好ましいシャフトの曲げ剛性（ＥＩ値）との関係を表す近似式を求めておく。この近似式は、複数のテスターに試打をしてもらうことでデータを収集し、このデータに基づいて作成する。近似式の信頼性を高める観点からは、テスターの数は多い方が好ましい。本適用例では、テスターを、ハンディが２０以下の中級者又は上級者とした。各テスターに対し、予め用意した複数のゴルフクラブ（ドライバー）の中から普段使用しているゴルフクラブの重量、長さ及びフレックス（又は調子）を基準にして数本のゴルフクラブを選定し、この数本のゴルフクラブについてそれぞれ数球、例えば６球試打（ミスショットを除く）してもらい、後述する基準にしたがって打ち易いゴルフクラブを選定してもらった。
【００８９】
打ち易いゴルフクラブは、「飛距離」、「方向性」及び「振り易さ」の３つの項目について、例えば以下の表５に示される評価基準にしたがい点数化し、得点合計が１．５点以上のクラブを打ち易いクラブと判断する。選定した数本のゴルフクラブは、シャフト部分が黒塗りであり、テスターに対してランダムに提供された。これにより、テスターはどのクラブで打っているのかを判断することができないようにした。試打は、３球×クラブ本数を２セット行い、振り易さのアンケートは、或るクラブについて３球試打したら、その都度点数を確認してもらった。
【００９０】
【表５】
2014147554.tif 000007
【００９１】
図２１は、前記のようにして予め作成された「６インチ」の測定点におけるスイング特徴量（４）と曲げ剛性（ＥＩ値）との関係を示す図であり、図２２は、「６インチ」の測定点におけるＥＩ値算出のフローを示す図である。なお、図２２において、「ωｙ４」は、コック方向の角速度のスイング特徴量（４）のことである。
【００９２】
本適用例では、ＥＩ値算出の精度を向上させるために、１つの近似式だけではなく３つの近似式を予め用意している。スイング特徴量（４）と「６インチ」の測定点におけるＥＩ値との関係を１つの近似式で表すことも可能であるが、トップ付近におけるｘ軸方向の角速度ωｘが極端に小さいゴルファーやヘッド速度がかなり小さいゴルファーなど平均的なスイングから外れたゴルファーに対しては、１つの近似式では信頼性の高いＥＩ値を算出することができない場合がある。そこで、平均的なゴルファーに対応する通常の近似式とともに、本適用例では、例外（１）近似式と、例外（２）近似式とを用意している。
【００９３】
まず、ステップＳ２１において、例外処理（１）が実行される。具体的に、トップ付近におけるｘ軸方向の角速度ωｘがωｘ＜−３００であり、且つ、トップからインパクトまでのｙ軸方向の角速度ωｙの平均値とトップからインパクトまでのｘ軸方向の角速度ωｘの平均値との比（ωｙ（全体）／ωｘ（全体））が０．５５以下であるか否かが判断され、Ｙｅｓの場合は、図２１において一点鎖線で示される例外（１）近似式を用いてＥＩ値が算出される。
【００９４】
一方、ステップＳ２１においてＮｏの場合は、ステップＳ２２に進んで例外処理（２）が実行される。この例外処理（２）では、以下の５つの条件（ａ）〜（ｅ）のうち少なくとも１つが満たされているか否かが判断される。
（ａ）トップからインパクトまでのｚ軸方向の角速度ωｚの累積値が極端に大きく、ωｚ＞２７０である。
（ｂ）ヘッド速度が４１．０ｍ／ｓ未満である。
（ｃ）トップ付近におけるｘ軸方向の角速度ωｘが、ωｘ＞０である。
（ｄ）トップからインパクトまでのｙ軸方向の角速度ωｙの平均値とトップからインパクトまでのｘ軸方向の角速度ωｘの平均値との比（ωｙ（全体）／ωｘ（全体））が２．０以上である。
（ｅ）トップからインパクトまでのｙ軸方向の角速度ωｙの平均値とトップからインパクトまでのｘ軸方向の角速度ωｘの平均値との比（ωｙ（全体）／ωｘ（全体））が１．５以上であり、且つ、スイング特徴量（４）が８００よりも大きい。
【００９５】
ステップＳ２２においてＹｅｓの場合は、図２１において破線で示される例外（２）近似式を用いてＥＩ値が算出される。一方、ステップＳ２２においてＮｏの場合は、図２１において実線で示される通常近似式を用いてＥＩ値が算出される。
【００９６】
通常近似式、例外（１）近似式及び例外（２）近似式は、いずれも最小二乗法による回帰直線を表す一次式とすることができる。各式の傾きと切片を次の表６に示す。
【００９７】
【表６】
2014147554.tif 000008
【００９８】
スイング特徴量（３）、スイング特徴量（２）及びスイング特徴量（１）についても、同様にして作成されたインチ毎の測定点におけるスイング特徴量と曲げ剛性（ＥＩ値）との関係を示す近似式を用いて曲げ剛性（ＥＩ値）を算出することができる。
【００９９】
<インチ毎のＩＦＣの算出>
前述した方法によりインチ毎に算出されたシャフトのＥＩ値について、例えば前記表１〜４に示される変換表を用いて、１０段階のＩＦＣのうちいずれかの値を算出する。本適用例では、前述したように、使用頻度などを考慮して、フィッターがユーザーに提供しようとして用意したシャフトの範囲内で１０段階に分ける方法を採用している。
【０１００】
<シャフトの選定>
以上のようにして、フィッティングを行うゴルファーについて、インチ毎のＩＦＣが算出される。算出されるＩＦＣの例として、３６インチ：５、２６インチ：４、１６インチ：４、６インチ：２をあげることができる。
【０１０１】
そして、この計算結果と最もマッチするシャフトをデータベースから選定する。データベースには、予め複数種類のシャフトについて計測されたインチ毎のＩＦＣ、重量などのデータが記憶されている。このデータベースを用いて、以下の式（３）に示される「一致度」をデータベースに記憶されている全てのクラブについて計算し、その値が最も小さいクラブをゴルファーに提案する。なお、本明細書における「一致度」は、一致している度合いを意味するのではなく、式（２）から明らかなように、その値が小さいほど計算結果に近い曲げ剛性を有するシャフトであるということを意味する指標である。一致度が０とは、計算結果と同じインチ毎の曲げ剛性を有するシャフトであるということである。
【０１０２】
【数１】
2014147554.tif 000009
【０１０３】
また、一致度が同程度に小さいシャフトが複数本存在する場合は、ユーザーに複数本提案してもよいし、ユーザーの要望を考慮して１本に絞り込んでもよい。絞込みの基準としては、振り易さを重視して、シャフト手元側の「３６インチ」又は「２６インチ」におけるＩＦＣの一致度を優先させる方法や、パフォーマンス（飛距離、方向性）を重視して、シャフト先端側の「１６インチ」又は「６インチ」におけるＩＦの一致度を優先させる方法がある。
【０１０４】
なお、今回開示された実施の形態はすべての点において単なる例示であって制限的なものではないと考えられるべきである。本発明の範囲は、前記した意味ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内のすべての変更が含まれることが意図される。
【符号の説明】
【０１０５】
１      ゴルフクラブ
２      センサ
２ａ    ケーシング
１０      コンピュータ（データ解析装置）
２０      シャフト
２０ａ    チップ端
２０ｂ    バット端
３０      レーザ発振器
３１      レーザ受光器
Ｇ      ゴルファー
Ｂ      ボール
Ｂｅ    ボール後端
【図１】
2014147554.tif 000010
【図２】
2014147554.tif 000011
【図３】
2014147554.tif 000012
【図４】
2014147554.tif 000013
【図５】
2014147554.tif 000014
【図６】
2014147554.tif 000015
【図１０】
2014147554.tif 000016
【図１２】
2014147554.tif 000017
【図１３】
2014147554.tif 000018
【図１４】
2014147554.tif 000019
【図１５】
2014147554.tif 000020
【図１６】
2014147554.tif 000021
【図１７】
2014147554.tif 000022
【図１８】
2014147554.tif 000023
【図１９】
2014147554.tif 000024
【図２０】
2014147554.tif 000025
【図２１】
2014147554.tif 000026
【図２２】
2014147554.tif 000027
【図２３】
2014147554.tif 000028
【図２４】
2014147554.tif 000029
【図２５】
2014147554.tif 000030
【図２６】
2014147554.tif 000031
【図７】
2014147554.tif 000032
【図８】
2014147554.tif 000033
【図９】
2014147554.tif 000034
【図１１】
2014147554.tif 000035