(57)【要約】【課題】ゴルファーに適したゴルフクラブを精度よく選定するためのフィッティング装置を提供する。【解決手段】ゴルファーに適したゴルフクラブを選定するためのフィッティング装置が提供される。前記フィッティング装置は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部とを備える。【選択図】図４
2017170105.tif 000002
【特許請求の範囲】
【請求項１】
ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備える、
フィッティング装置。
【請求項２】
前記算出部は、複数の種類の前記第１スイング指標を算出し、
前記決定部は、前記複数の種類の第１スイング指標の大きさに応じて、前記最適振り易さ指標を決定する、
請求項１に記載のフィッティング装置。
【請求項３】
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる、
請求項１又は２に記載のフィッティング装置。
【請求項４】
前記第１スイング指標には、前記スイング動作時のヘッドスピードが含まれる、
請求項１から３のいずれかに記載のフィッティング装置。
【請求項５】
前記決定部は、前記第１スイング指標が大きい又は小さい程、前記最適振り易さ指標を大きな値に決定する、
請求項１から４のいずれかに記載のフィッティング装置。
【請求項６】
前記テストクラブの種類毎に、前記第１スイング指標の大きさと前記最適振り易さ指標の大きさとの対応関係を定める対応関係データを記憶する記憶部
をさらに備え、
前記決定部は、前記テストクラブの種類に応じて、前記記憶部内の前記対応関係データを参照することにより、前記最適振り易さ指標を決定する、
請求項１から５のいずれかに記載のフィッティング装置。
【請求項７】
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
請求項１から６のいずれかに記載のフィッティング装置。
【請求項８】
前記計測値には、前記テストクラブのグリップエンドにおける加速度、角速度及び地磁気の少なくとも１つが含まれる、
請求項１から７のいずれかに記載のフィッティング装置。
【請求項９】
前記決定部は、前記最適剛性指標として、前記ゴルファーに適した、シャフトの複数の位置における剛性分布を決定する、
請求項１から８のいずれかに記載のフィッティング装置。
【請求項１０】
前記算出部は、前記第２スイング指標として、第１〜第３特徴量を算出し、
前記決定部は、前記第１〜第３特徴量のそれぞれの大きさに応じて、前記最適剛性指標として、前記シャフトのバット端からチップ端に向かってこの順に並ぶ第１〜第３位置における、前記ゴルファーに適したシャフトの剛性を示す第１〜第３剛性値を決定する、
請求項９に記載のフィッティング装置。
【請求項１１】
前記決定部は、前記第２スイング指標と前記最適剛性指標との関係を表す予め定められた式に従って、前記最適剛性指標を決定する、
請求項９又は１０に記載のフィッティング装置。
【請求項１２】
ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップと、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含む、
フィッティング方法。
【請求項１３】
ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップと、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させる、
フィッティングプログラム。
【請求項１４】
ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備える、
フィッティング装置。
【請求項１５】
前記決定部は、前記第１スイング指標の大きさに加え、前記ゴルフクラブに含まれるべきヘッドの種類に応じて、前記最適特性指標を決定し、
前記特定の部位は、前記ヘッド以外の部位である、
請求項１４に記載のフィッティング装置。
【請求項１６】
前記特定の部位は、シャフト又はグリップである、
請求項１４又は１５に記載のフィッティング装置。
【請求項１７】
前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさに応じて、前記最適特性指標を決定する、
請求項１４〜１６のいずれかに記載のフィッティング装置。
【請求項１８】
前記決定部は、前記第１スイング指標を説明変数とし、前記最適振り易さ指標を目的変数する所定の回帰式に、前記算出部により算出された前記第１スイング指標を代入することにより、前記最適振り易さ指標を算出する、
請求項１７に記載のフィッティング装置。
【請求項１９】
前記決定部は、前記第１スイング指標を説明変数とし、前記最適特性指標を目的変数する所定の回帰式に、前記算出部により算出された前記スイング指標を代入することにより、前記最適特性指標を算出する、
請求項１４〜１６のいずれかに記載のフィッティング装置。
【請求項２０】
前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である、
請求項１４〜１９のいずれかに記載のフィッティング装置。
【請求項２１】
ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
ユーザから前記ゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付ける入力受付部と、
前記ゴルファーのスイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量を決定する決定部と
を備える、
フィッティング装置。
【請求項２２】
前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記最適シャフト重量を決定する、
請求項２１に記載のフィッティング装置。
【請求項２３】
ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正する決定部と、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備える、
フィッティング装置。
【請求項２４】
ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定する決定部と、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備え、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる、
フィッティング装置。
【請求項２５】
ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップと、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含む、
フィッティング方法。
【請求項２６】
ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップと、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させる、
フィッティングプログラム。
【請求項２７】
ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含む、
フィッティング方法。
【請求項２８】
ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させる、
フィッティングプログラム。
【請求項２９】
ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含み、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる、
フィッティング方法。
【請求項３０】
ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をンピュータに実行させ、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる、
フィッティングプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴルファーに適したゴルフクラブ、特にシャフトを選定するためのフィッティング装置、方法及びプログラムに関する。
【背景技術】
【０００２】
従来より、ゴルファーにテストクラブを試打させてその動作を計測機器により計測し、当該計測値に基づいて当該ゴルファーに適したゴルフクラブを選定する様々なフィッティング方法が提案されている。その１つとして、特許文献１には、ゴルファーに適したゴルフクラブのシャフトを選定するためのフィッティング方法が開示されている。具体的には、特許文献１では、まずテストクラブによる計測値に基づいて、ゴルファーに適したシャフトの曲げ剛性が決定される。そして、データベースに登録されている多数のシャフトの中から、当該曲げ剛性の値に合致するシャフトが抽出される。かかる方法は、シャフトのフィッティングの精度を向上させる技術として、大いに期待されるものである。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−２２６３７５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、特許文献１の方法において、データベース内には、同じような曲げ剛性を有するシャフトが多数登録されていることも十分に想定される。そうすると、ゴルファーに適した曲げ剛性を決定するだけでは、ゴルファーに適したシャフトを十分に絞り込むことができない場合がある。特に、同じような曲げ剛性を有するが、重量帯の異なるシャフトが多数存在する場合が多々あり、どの重量帯のシャフトを選定してよいのか判断に困ることがある。
【０００５】
この点、特許文献１では、上述した曲げ剛性に加え、ゴルファーが普段使用しているゴルフクラブの重量（以下、マイクラブ重量）を基準に、シャフトの選定を行っている。具体的には、新たに提案されるゴルフクラブの重量がマイクラブ重量±５ｇとなるように、シャフトの選定を行っている。すなわち、新たに使用するゴルフクラブの重量がマイクラブ重量から大きく変化すると、タイミングがとりにくくなったり、振りにくく感じたりし、ゴルファーのパフォーマンスが低下する虞がある。この点、マイクラブ重量からの重量の変化が制限される特許文献１の方法は優れていると言える。
【０００６】
しかしながら、マイクラブ重量がそもそも当該ゴルファーに全く適していない、すなわち、オーバースペックだったり、アンダースペックだったりすることも考えられる。このような場合には、特許文献１の方法では、フィーリングが改善されても弾道が改善されなかったり、フィーリングが合わなかったりといった事態が生じ得る。
【０００７】
本発明は、ゴルファーに適したゴルフクラブを精度よく選定するためのフィッティング装置、方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
第１観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部とを備える。
【０００９】
第２観点に係るフィッティング装置は、第１観点に係るフィッティング装置であって、前記算出部は、複数の種類の前記第１スイング指標を算出する。前記決定部は、前記複数の種類の第１スイング指標の大きさに応じて、前記最適振り易さ指標を決定する。
【００１０】
第３観点に係るフィッティング装置は、第１観点又は第２観点に係るフィッティング装置であって、前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる。
【００１１】
第４観点に係るフィッティング装置は、第１観点から第３観点のいずれかに係るフィッティング装置であって、前記第１スイング指標には、前記スイング動作時のヘッドスピードが含まれる。
【００１２】
第５観点に係るフィッティング装置は、第１観点から第４観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標が大きい又は小さい程、前記最適振り易さ指標を大きな値に決定する。
【００１３】
第６観点に係るフィッティング装置は、第１観点から第５観点のいずれかに係るフィッティング装置であって、前記テストクラブの種類毎に、前記第１スイング指標の大きさと前記最適振り易さ指標の大きさとの対応関係を定める対応関係データを記憶する記憶部をさらに備える。前記決定部は、前記テストクラブの種類に応じて、前記記憶部内の前記対応関係データを参照することにより、前記最適振り易さ指標を決定する。
【００１４】
第７観点に係るフィッティング装置は、第１観点から第６観点のいずれかに係るフィッティング装置であって、前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる。
【００１５】
第８観点に係るフィッティング装置は、第１観点から第７観点のいずれかに係るフィッティング装置であって、前記計測値には、前記テストクラブのグリップエンドにおける加速度、角速度及び地磁気の少なくとも１つが含まれる。
【００１６】
第９観点に係るフィッティング装置は、第１観点から第８観点のいずれかに係るフィッティング装置であって、前記決定部は、前記最適剛性指標として、前記ゴルファーに適した、シャフトの複数の位置における剛性分布を決定する。
【００１７】
第１０観点に係るフィッティング装置は、第９観点に係るフィッティング装置であって、前記算出部は、前記第２スイング指標として、第１〜第３特徴量を算出する。前記決定部は、前記第１〜第３特徴量のそれぞれの大きさに応じて、前記最適剛性指標として、前記シャフトのバット端からチップ端に向かってこの順に並ぶ第１〜第３位置における、前記ゴルファーに適したシャフトの剛性を示す第１〜第３剛性値を決定する。
【００１８】
第１１観点に係るフィッティング装置は、第９観点又は第１０観点に係るフィッティング装置であって、前記決定部は、前記第２スイング指標と前記最適剛性指標との関係を表す予め定められた式に従って、前記最適剛性指標を決定する。
【００１９】
第１２観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００２０】
第１３観点に係るフィッティングプログラムは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００２１】
第１４観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部とを備える。
【００２２】
第１５観点に係るフィッティング装置は、第１４観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに加え、前記ゴルフクラブに含まれるべきヘッドの種類に応じて、前記最適特性指標を決定する。前記特定の部位は、前記ヘッド以外の部位である。
【００２３】
第１６観点に係るフィッティング装置は、第１４観点又は第１５観点に係るフィッティング装置であって、前記特定の部位は、シャフト又はグリップである。
【００２４】
第１７観点に係るフィッティング装置は、第１４観点から第１６観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさに応じて、前記最適特性指標を決定する。
【００２５】
第１８観点に係るフィッティング装置は、第１７観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標を説明変数とし、前記最適振り易さ指標を目的変数する所定の回帰式に、前記算出部により算出された前記第１スイング指標を代入することにより、前記最適振り易さ指標を算出する。
【００２６】
第１９観点に係るフィッティング装置は、第１４観点から第１６観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標を説明変数とし、前記最適特性指標を目的変数する所定の回帰式に、前記算出部により算出された前記スイング指標を代入することにより、前記最適特性指標を算出する。
【００２７】
第２０観点に係るフィッティング装置は、第１４観点から第１９観点のいずれかに係るフィッティング装置であって、前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である。
【００２８】
第２１観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、ユーザから前記ゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付ける入力受付部と、前記ゴルファーのスイング動作に関する第１スイング指標を算出する算出部と、前記第１スイング指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量を決定する決定部とを備える。
【００２９】
第２２観点に係るフィッティング装置は、第２１観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記最適シャフト重量を決定する。
【００３０】
第２３観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、取得部と、算出部と、決定部と、選択部とを備える。前記取得部は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する。前記算出部は、前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する。前記決定部は、前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正する。前記選択部は、前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する。
【００３１】
第２４観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、取得部と、算出部と、決定部と、選択部とを備える。前記取得部は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する。前記算出部は、前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する。前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定する。前記選択部は、前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。
【００３２】
第２５観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００３３】
第２６観点に係るフィッティングプログラムでは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００３４】
第２７観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。
（３）前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００３５】
第２８観点に係るフィッティングプログラムでは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。
（３）前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００３６】
第２９観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【００３７】
第３０観点に係るフィッティングプログラムは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。
【発明の効果】
【００３８】
本発明の第１観点によれば、テストクラブによる計測値に基づいて、第１スイング指標及び第２スイング指標が算出される。そして、第１スイング指標に応じて、ゴルファーに適したゴルフクラブの振り易さ指標（例えば、ゴルフクラブのスイング慣性モーメント、グリップエンド周りの慣性モーメント、重量等）である最適振り易さ指標が決定されるとともに、第２スイング指標に応じて、ゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適振り易さ指標及び最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方が選択される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適した振り易さ指標を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。
【００３９】
本発明の第１４観点によれば、テストクラブによる計測値に基づいて、第１スイング指標及び第２スイング指標が算出される。そして、第１スイング指標に応じて、ゴルファーに適したゴルフクラブの特定の部位の特性（例えば、シャフトの重量等）を表す最適特性指標が決定されるとともに、第２スイング指標に応じて、ゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適特性指標及び最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方が選択される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適したゴルフクラブの特定の部位の特性を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。
【００４０】
本発明の第２１観点によれば、ユーザからゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付けるとともに、ゴルファーのスイング動作に関する第１スイング指標が算出される。そして、ユーザにより指定されたヘッドの種類に加え、第１スイング指標に応じて、ゴルファーに適したシャフトの重量（最適シャフト重量）が決定される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適したシャフトの重量を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。
【図面の簡単な説明】
【００４１】
【図１】本発明の一実施形態に係るフィッティング装置を備えるフィッティングシステムを示す図。
【図２】フィッティングシステムの機能ブロック図。
【図３】ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。
【図４】フィッティング処理の流れを示すフローチャート。
【図５】（Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。
【図６】第１変換処理の流れを示すフローチャート。
【図７】スイング平面を説明する図。
【図８】肩挙動導出工程の流れを示すフローチャート。
【図９】二重振り子モデルを概念的に説明する図。
【図１０】第１指標算出工程の流れを示すフローチャート。
【図１１】二重振り子モデルを概念的に説明する別の図。
【図１２】腕エネルギーを説明する図。
【図１３】最適振り易さ決定工程の流れを示すフローチャート。
【図１４】プロモデル領域を示す図。
【図１５】アベレージモデル領域を示す図。
【図１６】インターナショナル・フレックス・コード（ＩＦＣ）を説明する図。
【図１７】シャフトの曲げ剛性の測定方法を説明する図。
【図１８】スイング中のシャフトの曲げを説明する図。
【図１９】第２実施形態に係るフィッティングシステムの機能ブロック図。
【図２０】第２実施形態に係るフィッティング処理の流れを示すフローチャート。
【図２１】最適特性決定工程の流れを示すフローチャート。
【図２２】様々なゴルフクラブのＩS，ＩGの関係をプロットしたグラフ。
【図２３】様々なゴルフクラブのＩS，ｍ2の関係をプロットしたグラフ。
【図２４】変形例に係る最適振り易さ決定工程の流れを示すフローチャート。
【図２５】別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。
【図２６】変形例に係る第１スイング指標となる角度を説明する図。
【図２７】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。
【図２８】図２７の変形例に係るプロモデル領域を示す図。
【図２９】図２７の変形例に係るアベレージモデル領域を示す図。
【図３０】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。
【図３１】図３０の変形例に係る平均肩トルク−腕エネルギー平面を示す図。
【図３２】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。
【図３３】変形例に係る肩挙動導出工程を概念的に説明する図。
【図３４】変形例に係る肩挙動導出工程を概念的に説明する別の図。
【図３５Ａ】第１実施形態に係る二重振り子モデルのシミュレーション結果を示す図。
【図３５Ｂ】変形例に係る二重振り子モデルのシミュレーション結果を示す図。
【図３６】最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。
【図３７Ａ】特定のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。
【図３７Ｂ】別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。
【図３７Ｃ】さらに別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。
【図３８】境界処理の流れを示す図。
【図３９Ａ】捩じりのローテーションの動きを示す概念図。
【図３９Ｂ】プッシュのローテーションの動きを示す概念図。
【図４０】例外処理の流れを示す図。
【図４１Ａ】実施例１及び比較例による飛距離を示すグラフ。
【図４１Ｂ】実施例１及び比較例による左右ぶれを示すグラフ。
【図４２Ａ】実施例２の試打結果を示すグラフ。
【図４２Ｂ】参考例の試打結果を示すグラフ。
【発明を実施するための形態】
【００４２】
以下、図面を参照しつつ、本発明の幾つかの実施形態に係るゴルフクラブのフィッティング装置、方法及びプログラムについて説明する。
【００４３】
＜１．第１実施形態＞
＜１−１．フィッティングシステムの概略構成＞
図１及び図２に、本実施形態に係るフィッティング装置２を備えるフィッティングシステム１００の全体構成を示す。フィッティング装置２は、ゴルファー７によるゴルフクラブ４のスイング動作を計測した計測データに基づいて、当該ゴルファー７に適したゴルフクラブ４を選定するための装置である。本実施形態では、スイング動作の計測は、ゴルフクラブ４のグリップ４２に取り付けられたセンサユニット１により行われ、フィッティング装置２は、このセンサユニット１とともに、フィッティングシステム１００を構成する。
【００４４】
以下、センサユニット１及びフィッティング装置２の構成について説明した後、フィッティング処理の流れについて説明する。
【００４５】
＜１−１−１．センサユニットの構成＞
センサユニット１は、図１及び図３に示すとおり、ゴルフクラブ４のグリップ４２におけるヘッド４１と反対側の端部に取り付けられており、グリップ４２の挙動を計測する。なお、ゴルフクラブ４は、一般的なゴルフクラブであり、シャフト４０と、シャフト４０の一端に設けられたヘッド４１と、シャフト４０の他端に設けられたグリップ４２とから構成される。本実施形態に係るシャフト４０は、カーボン製のシャフトである。センサユニット１は、スイング動作の妨げとならないよう、小型且つ軽量に構成されている。図２に示すように、本実施形態に係るセンサユニット１には、加速度センサ１１、角速度センサ１２及び地磁気センサ１３が搭載されている。また、センサユニット１には、これらのセンサ１１〜１３による計測データを外部のフィッティング装置２に送信するための通信装置１０も搭載されている。なお、本実施形態では、通信装置１０は、スイング動作の妨げにならないように無線式であるが、ケーブルを介して有線式にフィッティング装置２に接続するようにしてもよい。
【００４６】
加速度センサ１１、角速度センサ１２及び地磁気センサ１３はそれぞれ、グリップ４２を基準としたｘｙｚ局所座標系におけるグリップ加速度、グリップ角速度及びグリップ地磁気を計測する。より具体的には、加速度センサ１１は、ｘ軸、ｙ軸及びｚ軸方向のグリップ加速度ａx，ａy，ａzを計測する。角速度センサ１２は、ｘ軸、ｙ軸及びｚ軸周りの
グリップ角速度ωx，ωy，ωzを計測する。地磁気センサ１３は、ｘ軸、ｙ軸及びｚ軸方
向のグリップ地磁気ｍx，ｍy，ｍzを計測する。これらの計測データは、所定のサンプリ
ング周期Δｔの時系列データとして取得される。なお、ｘｙｚ局所座標系は、図３に示すとおりに定義される３軸直交座標系である。すなわち、ｚ軸は、シャフト４０の延びる方向に一致し、ヘッド４１からグリップ４２に向かう方向が、ｚ軸正方向である。ｘ軸は、ヘッド４１のトゥ−ヒール方向にできる限り沿うように配向され、ｙ軸は、ヘッド４１のフェース面の法線方向にできる限り沿うように配向される。
【００４７】
本実施形態では、加速度センサ１１、角速度センサ１２及び地磁気センサ１３による計測データは、通信装置１０を介してリアルタイムにフィッティング装置２に送信される。しかしながら、例えば、センサユニット１内の記憶装置に計測データを格納しておき、スイング動作の終了後に当該記憶装置から計測データを取り出して、フィッティング装置２に受け渡すようにしてもよい。
【００４８】
＜１−１−２．フィッティング装置の構成＞
図２を参照しつつ、フィッティング装置２の構成について説明する。フィッティング装置２は、ＣＤ−ＲＯＭ、ＵＳＢメモリ等のコンピュータで読み取り可能な記録媒体２０に格納された本実施形態に係るフィッティングプログラム３を、当該記録媒体２０から汎用のパーソナルコンピュータにインストールすることにより製造される。フィッティングプログラム３は、センサユニット１から送られてくる計測データに基づいてスイング動作を解析し、ゴルファー７に適したゴルフクラブ４を選択するのを支援する情報を出力するためのソフトウェアである。フィッティングプログラム３は、フィッティング装置２に後述する動作を実行させる。
【００４９】
フィッティング装置２は、表示部２１、入力部２２、記憶部２３、制御部２４及び通信部２５を備える。そして、これらの部２１〜２５は、バス線２６を介して接続されており、相互に通信可能である。本実施形態では、表示部２１は、液晶ディスプレイ等で構成され、後述する情報をユーザに対し表示する。なお、ここでいうユーザとは、ゴルファー７自身やそのインストラクター等の、フィッティングの結果を必要とする者の総称である。また、入力部２２は、マウス、キーボード、タッチパネル等で構成することができ、フィッティング装置２に対するユーザからの操作を受け付ける。通信部２５は、フィッティング装置２と外部装置との通信を可能にする通信インターフェースであり、センサユニット１からデータを受信する。
【００５０】
記憶部２３は、ハードディスク等の不揮発性の記憶装置により構成される。記憶部２３内には、フィッティングプログラム３が格納されている他、センサユニット１から送られてくる計測データが保存される。また、記憶部２３内には、対応関係データ２８、ヘッドデータベース（ＤＢ）２７及びシャフトデータベース（ＤＢ）２９が格納されている。対応関係データ２８とは、詳細は後述するが、ゴルフクラブ４の様々なモデル（シリーズ）毎に規定されており、最適振り易さ指標を決定するための条件を示すデータである。同様に詳細は後述するが、ヘッドＤＢ２７は、多数のヘッド４１のスペックを示す情報が、ヘッド４１の種類を特定する情報に関連付けて格納されたデータベースである。シャフトＤＢ２９は、多数のシャフト４０のスペックを示す情報が、シャフト４０の種類を特定する情報に関連付けて格納されたデータベースである。シャフトＤＢ２９には、例えば、１００を超える種類のシャフト４０のデータが格納されている。
【００５１】
制御部２４は、ＣＰＵ、ＲＯＭおよびＲＡＭ等から構成することができる。制御部２４は、記憶部２３内のフィッティングプログラム３を読み出して実行することにより、仮想的に取得部２４Ａ、グリップ挙動導出部２４Ｂ、肩挙動導出部２４Ｃ、算出部２４Ｄ、決定部２４Ｅ、選択部２４Ｆ及び表示制御部２４Ｇとして動作する。各部２４Ａ〜２４Ｇの動作の詳細については、後述する。
【００５２】
＜１−２．フィッティング処理＞
続いて、フィッティングシステム１００により実行されるフィッティング処理について説明する。本実施形態に係るフィッティング処理は、図４に示すとおり、以下の９つの工程（Ｓ１〜Ｓ９）から構成されている。
（Ｓ１）ｘｙｚ局所座標系でのグリップ加速度ａx，ａy，ａz、グリップ角速度ωx，ωy，ωz及びグリップ地磁気ｍx，ｍy，ｍzの計測データを計測する計測工程
（Ｓ２）計測工程で得られたｘｙｚ局所座標系での計測データを、ＸＹＺ全体座標系でのグリップ加速度ａX，ａY，ａZ及びグリップ角速度ωX，ωY，ωZに変換する第１変換工程（第１変換工程では、ＸＹＺ全体座標系でのグリップ速度ｖX，ｖY，ｖZも導出される。）
（Ｓ３）ＸＹＺ全体座標系でのグリップ４２の挙動（グリップ角速度ωX，ωY，ωZ及びグリップ速度ｖX，ｖY，ｖZ）を、スイング平面Ｐ（後述する）内でのグリップ４２の挙動へと変換する第２変換工程
（Ｓ４）スイング平面Ｐ内でのグリップ４２の挙動に基づいて、スイング平面Ｐ内でのゴルファー７の疑似的な肩の挙動を導出する肩挙動導出工程
（Ｓ５）スイング平面Ｐ内でのグリップ４２の挙動及び疑似的な肩の挙動に基づいて、第１スイング指標（本実施形態では、後述する腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖh）を算出する第１指標算出工程
（Ｓ６）第１スイング指標に基づいて、ゴルファー７に適した振り易さ指標（本実施形態では、後述するスイング慣性モーメントＩS）である最適振り易さ指標（本実施形態では、最適スイングＭＩ）を決定する最適振り易さ決定工程
（Ｓ７）計測データに基づいて、第２スイング指標（本実施形態では、後述する第１〜第４特徴量Ｆ1〜Ｆ4）を算出する第２指標算出工程
（Ｓ８）第２スイング指標に基づいて、ゴルファー７に適したシャフト４０の剛性を示す最適剛性指標（本実施形態では、後述するＥＩ分布）を決定する最適剛性決定工程
（Ｓ９）最適振り易さ指標及び最適剛性指標に合致するゴルフクラブ４、特にシャフト４０を選択する最適クラブ選択工程
以下、これらの工程を順に説明する。
【００５３】
なお、ＸＹＺ全体座標系は、図１に示すとおりに定義される３軸直交座標系である。すなわち、Ｚ軸は、鉛直下方から上方に向かう方向であり、Ｘ軸は、ゴルファー７の背から腹に向かう方向であり、Ｙ軸は、地平面に平行でボールの打球地点から目標地点に向かう方向である。
【００５４】
＜１−２−１．計測工程＞
計測工程（Ｓ１）では、ゴルファー７により、上述のセンサユニット１付きゴルフクラブ４がスイングされる。計測工程でスイングされるゴルフクラブ４は、２本のテストクラブのうちの１本である。これらのテストクラブは、異なる種類のゴルフクラブであり、本実施形態では、１本はプロ仕様のゴルフクラブ（以下、プロモデルクラブ）であり、もう１本はアベレージユーザーに適したゴルフクラブ（以下、アベレージモデルクラブ）である。また、本実施形態では、プロモデルクラブは、アベレージモデルクラブよりも重量が大きい。計測工程でいずれのテストクラブがスイングされるかは、ゴルファー７の好みや経験等に基づいて、決定される。
【００５５】
なお、フィッティングの精度を向上させるために、ユーザのマイクラブの長さが、シャフトＤＢ２９に記憶されているシャフトに基づくゴルフクラブの長さと異なる場合は、テストクラブの総重量をマイクラブの長さ相当のクラブ総重量に変更する。これにより、当該ユーザにマッチしたシャフトを選定することができる。例えば、シャフトＤＢ２９に記憶されているクラブ長さが全て４５インチ（＝１１４３ｍｍ）であり、ユーザのマイクラブの長さＡ（ｍｍ）が４５インチと異なる場合、テストクラブの総重量を次の式により算出したもの（４５インチ相当の総重量）に変更して、フィッティングを行う。
（テストクラブ総重量）＝（Ａ−１１４３）×０．３７７＋（マイクラブのクラブ総重量）
【００５６】
続いて、以上のようなゴルフクラブ４のスイング動作中のグリップ加速度ａx，ａy，ａz、グリップ角速度ωx，ωy，ωz及びグリップ地磁気ｍx，ｍy，ｍzの計測データが、センサユニット１により計測される。この計測データは、センサユニット１の通信装置１０を介してフィッティング装置２に送信される。一方、フィッティング装置２側では、取得部２４Ａが通信部２５を介してこれを受信し、記憶部２３内に格納する。本実施形態では、少なくともアドレスからインパクトまでの時系列の計測データが計測される。
【００５７】
なお、ゴルフクラブのスイング動作は、一般に、アドレス、トップ、インパクト、フィニッシュの順に進む。アドレスとは、図５（Ａ）に示すとおり、ゴルフクラブ４のヘッド４１をボール近くに配置した初期の状態を意味し、トップとは、図５（Ｂ）に示すとおり、アドレスからゴルフクラブ４をテイクバックし、最もヘッド４１が振り上げられた状態を意味する。インパクトとは、図５（Ｃ）に示すとおり、トップからゴルフクラブ４が振り下ろされ（ダウンスイング）、ヘッド４１がボールと衝突した瞬間の状態を意味し、フィニッシュとは、図５（Ｄ）に示すとおり、インパクト後、ゴルフクラブ４を前方へ振り抜いた状態を意味する。
【００５８】
計測工程では、以上のゴルフクラブ４が複数回、好ましくは５回以上試打されることが好ましい。この場合、計測データの平均値を算出し、以降の演算に使用することができる。また、ミスショットや計測ミス等による異常値を取り除くため、計測データの標準偏差σを算出するようにし、全ての計測データが平均値±ｋ・σ（ｋは、定数）以内に収まっていない場合には、計測の追加又はやり直しを求めるメッセージを表示部２１上に表示させるようにしてもよい。なお、計測データ自体の平均値ではなく、計測データに基づいて算出される加工値（例えば、後述する腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖh）の平均値を算出するようにしてもよい。加工値の平均値を算出する場合も、同じく標準偏差σに基づくデータの信頼性のチェックを行うことができる。
【００５９】
＜１−２−２．第１変換工程＞
以下、図６を参照しつつ、ｘｙｚ局所座標系の計測データをＸＹＺ全体座標系の値へと変換する第１変換工程（Ｓ２）について説明する。具体的には、まず、取得部２４Ａが、記憶部２３内に格納されているｘｙｚ局所座標系でのグリップ加速度ａx，ａy，ａz、グリップ角速度ωx，ωy，ωz及びグリップ地磁気ｍx，ｍy，ｍzの時系列の計測データを読み出す（ステップＳ１１）。
【００６０】
次に、ステップＳ１１で読み出されたｘｙｚ局所座標系での時系列の計測データに基づいて、グリップ挙動導出部２４Ｂが、インパクト、トップ及びアドレスの時刻ｔi，ｔt，ｔaを導出する（ステップＳ１２）。これらの時刻ｔi，ｔt，ｔaの導出方法としては、様々なものが公知であるため、ここでは詳細な説明を省略する。
【００６１】
続くステップＳ１３では、グリップ挙動導出部２４Ｂが、アドレスからインパクトまでの時刻ｔにおける姿勢行列Ｎ（ｔ）を算出する。今、姿勢行列を以下の式で表すとする。姿勢行列Ｎ（ｔ）は、時刻ｔにおけるＸＹＺ全体座標系をｘｙｚ局所座標系に変換するための行列である。
【数１】
2017170105.tif 000003
【００６２】
姿勢行列Ｎ（ｔ）の９つの成分の意味は、以下のとおりである。
成分ａ：全体座標系のＸ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｂ：全体座標系のＹ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｃ：全体座標系のＺ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｄ：全体座標系のＸ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｅ：全体座標系のＹ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｆ：全体座標系のＺ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｇ：全体座標系のＸ軸と、局所座標系のｚ軸とのなす角度の余弦
成分ｈ：全体座標系のＹ軸と、局所座標系のｚ軸とのなす角度の余弦
成分ｉ：全体座標系のＺ軸と、局所座標系のｚ軸とのなす角度の余弦
ここで、ベクトル（ａ，ｂ，ｃ）は、ｘ軸方向の単位ベクトルを表し、ベクトル（ｄ，ｅ，ｆ）は、ｙ軸方向の単位ベクトルを表し、ベクトル（ｇ，ｈ，ｉ）は、ｚ軸方向の単位ベクトルを表している。
【００６３】
また、姿勢行列Ｎ（ｔ）は、Ｚ−Ｙ−Ｚ系のオイラー角の考え方に従うと、以下の式で表すことができる。ただし、φ，θ，ψは、それぞれＺ軸、Ｙ軸、Ｚ軸周りの回転角度とする。
【数２】
2017170105.tif 000004
【００６４】
アドレスからインパクトまでの姿勢行列Ｎ（ｔ）を算出するに当たり、まず、アドレスの時刻ｔaにおける姿勢行列Ｎ（ｔa）が算出される。具体的には、以下の式に従って、アドレス時のφ，θが算出される。なお、以下の式は、アドレス時にはゴルフクラブ４は静止しており、加速度センサ１１によって鉛直方向の重力のみが検出されることを利用している。以下の式中のグリップ加速度ａx，ａy，ａzは、アドレス時の値である。
【数３】
2017170105.tif 000005
【数４】
2017170105.tif 000006
【００６５】
続いて、以下の式に従って、アドレス時のψが算出される。
【数５】
2017170105.tif 000007
ただし、上式中のｍxi，ｍyiの値は、以下の式に従って算出される。また、以下の式中のグリップ地磁気ｍx，ｍy，ｍzは、アドレス時の値である。
【数６】
2017170105.tif 000008
【００６６】
以上より、アドレス時のφ，θ，ψが、ｘｙｚ局所座標系でのグリップ加速度ａx，ａy，ａz及びグリップ地磁気ｍx，ｍy，ｍzに基づいて算出される。そして、これらのφ，θ，ψの値を数２の式に代入することにより、アドレス時の姿勢行列Ｎ（ｔa）が算出される。
【００６７】
続いて、アドレス時の姿勢行列Ｎ（ｔa）をサンプリング周期Δｔ間隔で時々刻々更新してゆくことにより、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）が算出される。具体的に説明すると、まず、姿勢行列Ｎ（ｔ）は、クォータニオンの４変数ｑ1，ｑ2，ｑ3，ｑ4（ｑ4がスカラー部）を用いて、以下の式で表される。
【数７】
2017170105.tif 000009
【００６８】
従って、数１及び数７より、クォータニオンの４変数ｑ1，ｑ2，ｑ3，ｑ4は、以下の式に従って、算出することができる。
【数８】
2017170105.tif 000010
【００６９】
今、アドレス時の姿勢行列Ｎ（ｔa）を規定するａ〜ｉの値は既知である。よって、以上の式に従って、まず、アドレス時のクォータニオンの４変数ｑ1，ｑ2，ｑ3，ｑ4が算出される。
【００７０】
そして、時刻ｔから微小時刻経過後のクォータニオンｑ'は、時刻ｔにおけるクォータニオンｑを用いて以下の式で表される。
【数９】
2017170105.tif 000011
【００７１】
また、クォータニオンの４変数ｑ1，ｑ2，ｑ3，ｑ4の時間変化を表す１階微分方程式は、以下の式で表される。
【数１０】
2017170105.tif 000012
【００７２】
数９及び数１０の式を用いれば、時刻ｔのクォータニオンを順次、次の時刻ｔ＋Δｔのクォータニオンへと更新することができる。ここでは、アドレスからインパクトまでのクォータニオンが算出される。そして、アドレスからインパクトまでのクォータニオンの４変数ｑ1，ｑ2，ｑ3，ｑ4を数７の式に順次代入してゆくことにより、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）が算出される。
【００７３】
次に、ステップＳ１４では、グリップ挙動導出部２４Ｂが、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）に基づいて、アドレスからインパクトまでのｘｙｚ局所座標系でのグリップ加速度ａx，ａy，ａz及びグリップ角速度ωx，ωy，ωzの時系列データを、ＸＹＺ全体座標系での時系列データに変換する。変換後のグリップ加速度ａX，ａY，ａZ及びグリップ角速度ωX，ωY，ωZは、以下の式に従って算出される。
【数１１】
2017170105.tif 000013
【００７４】
続くステップＳ１５では、グリップ挙動導出部２４Ｂは、グリップ加速度ａX，ａY，ａZの時系列データを積分することにより、アドレスからインパクトまでのＸＹＺ全体座標系でのグリップ速度ｖX，ｖY，ｖZを導出する。このとき、アドレスからインパクトまでのグリップ速度ｖX，ｖY，ｖZを、トップにおいて０ｍ／ｓとなるように、オフセットを行うことが好ましい。
【００７５】
＜１−２−３．第２変換工程＞
以下、第１変換工程で算出されたＸＹＺ全体座標系でのグリップ４２の挙動を、スイング平面Ｐ内でのグリップ４２の挙動へと変換する第２変換工程（Ｓ３）について説明する。本実施形態では、スイング平面Ｐは、ＸＹＺ全体座標系の原点を含み、Ｙ軸及びインパクト時のシャフト４０と平行な面として定義される（図７参照）。第２変換工程では、グリップ挙動導出部２４Ｂは、ＸＹＺ全体座標系でのグリップ速度ｖX，ｖY，ｖZ及びグリップ角速度ωX，ωY，ωZをスイング平面Ｐ内へ射影したグリップ速度ｖpX，ｖpY，ｖpZ及びグリップ角速度ωpX，ωpY，ωpZを算出する。
【００７６】
具体的には、シャフト４０の延びる方向を表す、姿勢行列Ｎ（ｔ）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）に基づいて、Ｘ軸正方向から見た（ゴルファー７を正面から見た）シャフト４０の傾きの時系列データを算出する。そして、この時系列データに基づいて、Ｘ軸正方向から見てシャフト４０がＺ軸と平行になる時刻を特定し、これをインパクトの時刻ｔiとする。なお、ここでのインパクトの時刻ｔiは、既出のインパクトの時刻ｔiと一致するとは限らない。続いて、このインパクトの時刻ｔiにおける姿勢行列Ｎ（ｔi）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）に基づいて、Ｙ軸負方向から見たシャフト４０の傾きを算出する。すなわち、インパクト時にＹ軸負方向から見たシャフト４０とＸ軸との為す角度α'を算出し、これをスイング平面角度とする。
【００７７】
スイング平面角度α'が求まると、これを用いてＸＹＺ全体座標系における任意の点をスイング平面Ｐに射影するための射影変換行列Ａを、以下のとおり算出することができる。ただし、α＝９０°−α'である。
【数１２】
2017170105.tif 000014
【００７８】
ここでは、以上の射影変換行列Ａに基づいて、以下の式に従って、アドレスからインパクトまでの射影変換後のグリップ速度ｖpX，ｖpY，ｖpZ及びグリップ角速度ωpX，ωpY，ωpZの時系列データが算出される。
【数１３】
2017170105.tif 000015
【００７９】
なお、以上の演算により得られるグリップ速度（ｖpY，ｖpZ）は、スイング平面Ｐ内でのグリップ速度（ベクトル）を表しており、グリップ角速度ωpXは、スイング平面Ｐに対して垂直な軸周りの角速度を表している。ここでは、以下の式に従って、アドレスからインパクトまでのスイング平面Ｐ内でのグリップ速度（スカラー）が算出される。
【数１４】
2017170105.tif 000016
【００８０】
また、ここでは、後の計算に必要となる、スイング平面Ｐ内におけるトップでのシャフト４０の傾きも算出される。具体的には、まず、トップでの姿勢行列Ｎ（ｔt）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）を、射影変換行列Ａを用いて、以下の式に従ってスイング平面Ｐ内に射影する。ただし、射影後のベクトルを（ｇ'，ｈ'，ｉ'）とする。
【数１５】
2017170105.tif 000017
【００８１】
以上の式により特定されるベクトル（ｈ'，ｉ'）は、スイング平面Ｐ内におけるトップでのシャフト４０の傾きを表すベクトルである。従って、以上の計算結果を以下の式に代入することにより、スイング平面Ｐ内におけるトップでのシャフト４０の傾きβが算出される。
【数１６】
2017170105.tif 000018
【００８２】
＜１−２−４．肩挙動導出工程＞
以下、図８を算出しつつ、スイング平面Ｐ内でのグリップの挙動（グリップ速度ＶGE及びグリップ角速度ωpX）に基づいて、スイング平面Ｐ内の疑似的な肩の挙動を導出する肩挙動導出工程（Ｓ４）について説明する。本実施形態では、ゴルフクラブ４の挙動は、ゴルファー７の肩及びグリップ４２（或いは、これを握るゴルファーの手首）を節点とし、ゴルファー７の腕及びゴルフクラブ４をリンクとする二重振り子モデルに基づいて解析される。ただし、肩の挙動は直接的に実測されるのではなく、実測されたグリップの挙動に基づいて、疑似的な肩の挙動として導出される。以下では、特に断らない限り、単に「肩」という場合も、このような疑似的な肩を意味するものとする。疑似的な肩とグリップ４２（手首）との間を直線的に延びるものとして定義される疑似的な「腕」についても同様である。
【００８３】
グリップの挙動から肩の挙動を特定するに当たり、本実施形態に係る二重振り子モデルは、以下の（１）〜（５）を前提とする。図９は、以下の前提条件を概念的に説明する図である。
（１）スイング平面Ｐ上において、グリップ４２（手首）は肩を中心として円運動する。（２）スイング平面Ｐ上において、肩とグリップ４２との距離（半径）Ｒは、一定である。
（３）肩は、スイング動作中は動かない（ただし、回転する）。
（４）スイング平面Ｐ上において、トップでの腕とゴルフクラブ４との為す角度は９０°である。
（５）インパクト時の腕は、Ｘ軸正方向から見てＺ軸下方を向く。
【００８４】
以上の前提の下、肩挙動導出部２４Ｃは、スイング平面Ｐ内におけるトップからインパクトまでのグリップ４２の移動距離Ｄを算出する（ステップＳ２１）。移動距離Ｄは、トップからインパクトまでのグリップ速度ＶGEを積分することにより導出される。
【００８５】
続いて、肩挙動導出部２４Ｃは、スイング平面Ｐ内におけるトップからインパクトまでの腕の回転角度γを算出する（ステップＳ２２）。回転角度γは、第２変換工程で算出されたトップでのシャフト４０の傾きβに基づいて算出される。次に、肩挙動導出部２４Ｃは、半径Ｒ＝Ｄ／γを算出する（ステップＳ２３）。
【００８６】
そして、肩挙動導出部２４Ｃは、以下の式に従って、肩の挙動として、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω1を算出する（ステップＳ２４）。すなわち、腕の角速度ω1は、計測によるグリップ速度ＶGEが反映さ
れた値となる。
ω1＝ＶGE／Ｒ
【００８７】
＜１−２−５．第１指標算出工程＞
以下、図１０を参照しつつ、グリップ４２の挙動及び肩の挙動に基づいて、第１スイング指標を算出する第１指標算出工程（Ｓ５）について説明する。第１スイング指標とは、最適振り易さ指標を決定するための、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態では、第１スイング指標として、後述される腕出力パワーＰ1_AVE
、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhが算出される。
【００８８】
具体的には、まず、ステップＳ３１では、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω1を積分し、トップからインパクトまでの腕の回転角度θ1を算出する。このとき、台形積分を用いることが好ましい。なお、回転角度θ1は、図１１のように定義され、図１１の紙面は、スイング平面Ｐに等しい。以下では、図１１に示されるスイング平面Ｐ内での新たなＸＹ座標系に基づいて、解析が進められる。スイング平面Ｐ内での新たなＸＹ座標系のＸ軸は、上述したＸＹＺ全体座標系のＹ軸に等しく、新たなＸＹ座標系のＹ軸は、ＸＹＺ全体座標系のＺ軸をスイング平面Ｐ内に投影した軸である。
【００８９】
また、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω1を微分し、トップからインパクトまでの角加速度ω1'を算出する。次に、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の重心の位置（Ｘ1，Ｙ1）、速度（ＶX1，ＶY1）及び加速度（ＡX1，ＡY1）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。
【数１７】
2017170105.tif 000019
【００９０】
ただし、ｒは、肩から腕の重心までの距離である。本実施形態では、腕の重心は、腕の中心にあるものと仮定される。従って、Ｒ＝２ｒである。
【００９１】
次に、ステップＳ３２では、グリップ挙動導出部２４Ｂは、ステップＳ３１と同様の演算をグリップ４２周りについても行う。すなわち、トップからインパクトまでのグリップ角速度ωpX＝グリップ４２周りのゴルフクラブ４の角速度ω2を積分し、トップからインパクトまでのグリップ４２周りのゴルフクラブ４（シャフト４０）の回転角度θ2を算出する。このときも、台形積分を用いることが好ましく、回転角度θ2は、図１１のように定義される。
【００９２】
続いて、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の角速度ω2を微分し、トップからインパクトまでの角加速度ω2'を算出する。次に、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の重心の位置（Ｘ2，Ｙ2）、速度（ＶX2，ＶY2）及び加速度（ＡX2，ＡY2）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。
【数１８】
2017170105.tif 000020
【００９３】
ただし、Ｌは、グリップ４２からゴルフクラブ４の重心までの距離である。Ｌの値は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【００９４】
次に、ステップＳ３３では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの肩に発生する拘束力Ｒ1＝（ＲX1，ＲY1）を算出するとともに、トップからインパクトまでのグリップ４２に発生する拘束力Ｒ2＝（ＲX2，ＲY2）を算出する。以下の式は、並進方向の力の釣り合いに基づくものである。ただし、ｍ1は、腕の質量であり、本実施形態では、腕の質量ｍ1は、適宜予め定められているものとする。例えば、解析を開始する前に、ゴルファー７の体重を入力しておき、入力された体重に所定の係数を掛ける等して、自動的に腕の質量が算出される。ｍ2は、ゴルフクラブ４の質量であり、ｇは、重力加速度である。また、ｍ2は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【数１９】
2017170105.tif 000021
【００９５】
続くステップＳ３４では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の重心周りのトルクＴg1及びゴルフクラブ４の重心周りのトルクＴg2を算出する。
【数２０】
2017170105.tif 000022
ただし、Ｉ1は、腕の重心周りの慣性モーメントであり、Ｉ2は、ゴルフクラブ４の重心周りの慣性モーメントである。本実施形態では、腕の重心周りの慣性モーメントＩ1は、腕の重心が腕の中心にあるとの仮定の下、Ｉ1＝ｍ1ｒ2／３として算出される。また、I2は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【００９６】
続くステップＳ３５では、算出部２４Ｄは、上述した計算結果に基づいて、トップからインパクトまでの腕の仕事率（パワー）Ｅ1'を算出する。具体的には、Ｅ1'は、肩の速度ベクトルをｖsとし、グリップ４２の速度ベクトルをｖgとして、以下の式に従って表される。また、ｖs，ｖgはそれぞれ、肩の位置ベクトルｄs、グリップ４２の位置ベクトルｄg＝ｄs＋（２Ｘ1，２Ｙ1）を一階微分することにより算出可能である。
【数２１】
2017170105.tif 000023
【００９７】
また、本実施形態では、肩は動かないため、ｖs＝（０，０）となり、腕の仕事率Ｅ1'は、以下の式に従って算出される。算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の仕事率Ｅ1'を算出する。
【数２２】
2017170105.tif 000024
【００９８】
ところで、ゴルフスイングにおいて、ゴルフクラブ４の先端（ヘッド４１）を最も加速させるためには、まず腕を十分に加速させて、その後、腕の動きを止めてゴルフクラブ４に勢いを与えることが求められると考えられる。ここでいう腕の加速具合とは、腕が出力するパワー（腕出力パワー）Ｐ1という物理指標に置き換えることができ、ゴルフクラブ４に与える勢いとは、ゴルフクラブ４に入力されるパワー（クラブ入力パワー）Ｐ2という物理指標に置き換えることができる。そして、腕出力パワーＰ1とは、腕の仕事率Ｅ1'を表す数２２の右辺の第２項及び第３項部分に相当する。また、クラブ入力パワーＰ2とは、数２２の式中の右辺の第１項部分に相当する。すなわち、腕出力パワーＰ1及びクラブ入力パワーＰ2は、以下のとおり表すことができる。ステップＳ３５では、算出部２４Ｄは、腕の仕事率Ｅ1'に加え、トップからインパクトまでの腕出力パワーＰ1及びクラブ入力パワーＰ2を算出する。
【数２３】
2017170105.tif 000025
【００９９】
なお、スイング動作中にゴルフクラブ４で発揮される仕事率Ｅ2'は、下式のように表すことができる。すなわち、クラブ入力パワーＰ2＝Ｒ2ｖgTが橋渡しになって、腕からゴルフクラブ４へとエネルギーが伝達される。
【数２４】
2017170105.tif 000026
【０１００】
続くステップＳ３６では、算出部２４Ｄは、トップ以降で腕の仕事率Ｅ1'が正から負へ転じる時刻ｔcを特定し、トップの時刻ｔtから時刻ｔcまでの腕の仕事量Ｅ1を算出する。腕の仕事量Ｅ1は、時刻ｔt〜ｔcの区間で腕の仕事率Ｅ1'を積分することにより、算出される（図１２参照）。なお、仕事量Ｅ1は、時刻ｔt〜ｔcの間に腕で発揮される仕事量（エネルギー）を表す指標と考えることができるから、この意味で、スイング動作中の腕エネルギーと呼ぶことができる。また、算出部２４Ｄは、時刻ｔt〜ｔcの間に腕で発揮される平均仕事率ＥAVE＝Ｅ1／ｔc―ｔtを算出する。平均仕事率ＥAVEは、スイング動作中に単位時間当たりに平均的に発揮ないし消費される腕エネルギーである。
【０１０１】
また、算出部２４Ｄは、トップの時刻ｔtから腕出力パワーＰ1が最大値をとる時刻ｔmまでの区間で腕出力パワーＰ1を積分し、この積分値Ｄ1を積分区間（ｔc―ｔm）で除算することで、スイング動作中の平均的な腕出力パワーＰ1_AVEを算出する。なお、この積分値Ｄ1は、スイング動作中にゴルファーの腕が行う仕事量であり、腕出力パワーを示す指標となり得る。同様に、算出部２４Ｄは、トップの時刻ｔtからクラブ入力パワーＰ2が最大値をとる時刻ｔnまでの区間でクラブ入力パワーＰ2を積分し、この積分値Ｄ2を積分区間（ｔc―ｔn）で除算することで、スイング動作中の平均的なクラブ入力パワーＰ2_AVEを算出する。なお、この積分値Ｄ2は、スイング動作中のゴルフクラブ４に対してされる仕事量であり、クラブ入力パワーを示す指標となり得る。なお、ここで示した積分区間は例示であり、例えば、時刻ｔt〜ｔcまでの区間等、適宜設定し得る。
【０１０２】
続くステップＳ３７では、算出部２４Ｄは、スイング動作中のコック解放タイミングｔrを算出する。なお、本発明者らは、実験を通して、インパクト時のヘッド速度Ｖhが、スイング動作中のコック解放タイミングｔr、及び、腕エネルギーＥ1又は平均仕事率ＥAVEと相関があることを発見した。そこで、ここでは、インパクト時のヘッド速度Ｖhを算出すべく、コック解放タイミングｔrが算出される。本実施形態では、コック解放タイミングｔrは、時刻ｔt〜ｔcまでの区間で腕の仕事率Ｅ1'が最大となる時刻が、コック解放タイミングｔrとして特定される（図１２参照）。
【０１０３】
続くステップＳ３８では、算出部２４Ｄは、コック解放タイミングｔr及び腕エネルギーＥAVEに基づいて、インパクト時のヘッド速度Ｖhを算出する。具体的には、インパクト時のヘッド速度Ｖhは、下の式に従って算出される。なお、ｋ1，ｋ2，ｋ3は、予め行われた多数の実験結果から重回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。以上より、指標算出工程が終了する。
Ｖh＝ｋ1・ＥAVE＋ｋ2・ｔr＋ｋ3
【０１０４】
以上の処理により、最適振り易さ指標を決定するための第１スイング指標として、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhが算出される。
【０１０５】
＜１−２−６．最適振り易さ決定工程＞
以下、図１３を参照しつつ、最適振り易さ決定工程（Ｓ６）の流れについて説明する。本実施形態に係る最適振り易さ指標とは、最適スイングＭＩであり、最適スイングＭＩとは、ゴルファーに適したゴルフクラブ４のスイング慣性モーメントＩSのことを言う。
【０１０６】
また、スイング慣性モーメントＩSとは、スイング中の肩周りの慣性モーメントであり、例えば、以下の式に従って定義することができる。
ＩS＝Ｉ2＋ｍ2（Ｒ＋Ｌ）2＋Ｉ1＋ｍ1（Ｒ／２）2
なお、各ゴルファー７については、ゴルフクラブ４が変わっても腕の重量は同じである。従って、本実施形態では、簡単のため、スイング慣性モーメントＩSは、腕の回転分の慣性モーメントを省略し、以下の式に従って算出される。
ＩS＝Ｉ2＋ｍ2（Ｒ＋Ｌ）2
さらに、本実施形態では、腕長さＲ＝６０ｃｍ（一定）として、スイング慣性モーメントＩSが計算される。ただし、上式中のＲには、ステップＳ２３で算出される腕長さＲの値を代入することもできる。ところで、ＩSを決定するパラメータであるｍ2，Ｉ2，Ｌは、ゴルフクラブ４の諸元である。従って、本実施形態でのスイング慣性モーメントも、ゴルフクラブ４の諸元となる。
【０１０７】
まず、ステップＳ４０において、決定部２４Ｅは、計測工程で試打されたテストクラブの種類を判定する。プロモデルクラブが試打されていた場合には、ステップＳ４１に進み、アベレージモデルクラブが試打されていた場合には、ステップＳ５１に進む。いずれのテストクラブが試打されたかは、入力部２２を介してユーザにより入力される情報に基づいて判定されるものとする。
【０１０８】
なお、ステップＳ４０においてテストクラブの種類が判定されるのは、テストクラブの種類に応じて第１スイング指標の分布する領域が異なるからである。具体的に説明すると、本発明者らは、プロモデルのゴルフクラブを通常使用しているゴルファー（以下、プロモデルユーザ）２１名及びアベレージモデルのゴルフクラブを通常使用しているゴルファー（以下、アベレージモデルユーザ）２２名に、それぞれプロモデルクラブ及びアベレージモデルクラブを試打させて第１スイング指標を算出したところ、図１４及び図１５に示す結果を得た。なお、この実験で算出された第１スイング指標は、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhであり、具体的な値は、上記と同様の工程に従って算出された。そして、本発明者らは、図１４及び図１５に示す結果から、第１スイング指標を示す空間は、プロモデル領域とアベレージモデル領域とに概ね分割されることを発見した。なお、プロモデル領域とは、プロモデルユーザによるスイング動作時の第１スイング指標が分布する領域であり、アベレージモデル領域とは、アベレージモデルユーザによるスイング動作時の第１スイング指標が分布する領域である。図１４及び図１５を比較すれば分かるように、プロモデル領域は、アベレージモデル領域よりも、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhがいずれも大きい領域である。しかしながら、これらの数値は、テストクラブの種類等のスイング条件によって異なるものとなり得る。本実験では、プロモデルクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ  Ｋｏｓｕｍａ  Ｂｌｕｅ６  Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１４ｇ、バランスは、Ｄ３）のドライバーが用いられ、アベレージモデルクラブとして、ダンロップスポーツ株式会社製のＸＸＩＯ（登録商標）８（シャフトは、ＭＰ−８００  Ｒ−Ｆｌｅｘ、ゴルフクラブの重量は、２７２ｇ、バランスは、Ｄ３）のドライバーが用いられた。
【０１０９】
続いて、ステップＳ４１及びこれに続くステップＳ４２〜Ｓ４４について説明する。ステップＳ４１〜Ｓ４４は、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEの大きさに応じて、最適スイングＭＩの範囲（以下、最適スイングＭＩ帯）を決定するステップである。ここでは、Ｐ1_AVE，Ｐ2_AVEの値が大きい程、最適スイングＭＩ帯は段階的に大きな値に設定される。
【０１１０】
具体的には、ステップＳ４１では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ１の上側にあるか、すなわち、図１４の領域Ａ１に属するか否か（以下、プロ条件１）を判定する。そして、プロ条件１が満たされる場合には、５６２３ｋｇ・ｃｍ2以上が最適スイングＭＩ帯であると判定する。一方、ステップＳ４１でプロ条件１が満たされない場合には、ステップＳ４２に進む。ステップＳ４２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ２の上側かつ直線Ｌ１よりも下側にあるか、すなわち、図１４の領域Ａ２に属するか否か（以下、プロ条件２）を判定する。そして、プロ条件２が満たされる場合には、５５４０〜５６２３ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ４２でプロ条件２が満たされない場合には、ステップＳ４３に進む。ステップＳ４３では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ３の上側かつ直線Ｌ２よりも下側にあるか、すなわち、図１４の領域Ａ３に属するか否か（以下、プロ条件３）を判定する。そして、プロ条件３が満たされる場合には、５４８０〜５５４０ｋｇ・ｃｍ2が最適スイングＭ帯Ｉであると判定する。一方、ステップＳ４３でプロ条件３が満たされない場合には、ステップＳ４４に進む。ステップＳ４４では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ４の上側かつ直線Ｌ３よりも下側にあるか、すなわち、図１４の領域Ａ４に属するか否か（以下、プロ条件４）を判定する。そして、プロ条件４が満たされる場合には、５４００〜５４８０ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ４４でプロ条件４が満たされない場合、すなわち、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ４よりも下側にあり、図１４の領域Ａ５に属する場合には、５４００ｋｇ・ｃｍ2以下が最適スイングＭＩ帯であると判定する。
【０１１１】
一方、ステップＳ５１及びこれに続くステップＳ５２，Ｓ５３は、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhの大きさに応じて、最適スイングＭＩ帯を決定するステップである。ここでは、Ｐ1_AVE，Ｐ2_AVE及びＶhの値が大きい程、最適スイングＭＩ帯は段階的に大きな値に設定される。
【０１１２】
具体的には、決定部２４Ｅは、ステップＳ５１において、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１５に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ５の上側にあるか、すなわち、図１５の領域Ｂ１に属するか否か（以下、アベレージ条件１）を判定する。そして、アベレージ条件１が満たされる場合には、５３００ｋｇ・ｃｍ2以上が最適スイングＭＩ帯であると判定する。一方、ステップＳ５１でアベレージ条件１が満たされない場合には、ステップＳ５２に進む。ステップＳ５２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１５に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ６の上側かつ直線Ｌ５よりも下側にあるか、すなわち、図１５の領域Ｂ２に属するか否か（以下、アベレージ条件２）を判定する。そして、アベレージ条件２が満たされる場合には、５２００〜５３００ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ５２でアベレージ条件２が満たされない場合には、ステップＳ５３に進む。ステップＳ５３では、決定部２４Ｅは、第１指標算出工程で算出されたヘッド速度Ｖhが３５．５ｍ／ｓより大きいか否か（アベレージ条件３）を判定する。そして、アベレージ条件３が満たされる場合には、５２００〜５３００ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ５３でアベレージ条件３が満たされない場合には、すなわち、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１５に示す直線Ｌ６よりも下側にあり、図１４の領域Ｂ３に属する場合には、５２００ｋｇ・ｃｍ2以下が最適スイングＭＩ帯であると判定する。
【０１１３】
以上のステップＳ４１〜Ｓ４４及びステップＳ５１〜Ｓ５３は、以下の知見に基づく。すなわち、図１４及び図１５を参照して説明した上記実験では、ヘッド速度Ｖh及び、飛距離を最大にする最適スイングＭＩも算出された。ヘッド速度Ｖhは、上記と同様の工程に従って算出された。一方、最適スイングＭＩについては、ゴルファーに様々なスイング慣性モーメントＩsのゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブのスイング慣性モーメントＩsを特定し、これを最適スイングＭＩとした。より具体的には、プロモデルユーザには、スイング慣性モーメントＩsが５６５０ｋｇ・ｃｍ2、５６１０ｋｇ・ｃｍ2、５５２０ｋｇ・ｃｍ2、５４４０ｋｇ・ｃｍ2、５３６０ｋｇ・ｃｍ2の５種類のゴルフクラブをスイングさせた。一方、アベレージモデルユーザには、スイング慣性モーメントＩsが５３３０ｋｇ・ｃｍ2、５２３０ｋｇ・ｃｍ2、５１３０ｋｇ・ｃｍ2の３種類のゴルフクラブをスイングさせた。図１４及び図１５には、本実験により得られたヘッド速度Ｖh[ｍ／ｓ]及び最適スイングＭＩの値も示されている。
【０１１４】
図１４及び図１５に示すとおり、上記の実験からは、最適スイングＭＩは、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEが共に大きいほど、大きくなることが分かる。その結果、本発明者らは、腕出力パワーＰ1_AVE−クラブ入力パワーＰ2_AVE空間（プロモデル領域）を、図１４に示すように領域分割することで、最適スイングＭＩ帯に対応する領域Ａ１〜Ａ５を定義可能であることを発見した。また、腕出力パワーＰ1_AVE−クラブ入力パワーＰ2_AVE−ヘッド速度Ｖh空間（アベレージモデル領域）を、図１５に示すように領域分割することで、最適スイングＭＩ帯に対応する領域Ｂ１〜Ｂ３を定義し得ることを発見した。ただし、簡単のため、図１５では、ヘッド速度Ｖhを表す軸は省略され、腕出力パワーＰ1_AVE−クラブ入力パワーＰ2_AVE平面が示されている。すなわち、上述のステップＳ４１〜Ｓ４４は、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEを示す点が、Ｐ1_AVE−Ｐ2_AVE空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定するステップとなっている。一方、ステップＳ５１〜Ｓ５３は、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhを示す点が、Ｐ1_AVE−Ｐ2_AVE−Ｖh空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定するステップとなっている。ステップＳ４１〜Ｓ４４及びステップＳ５１〜Ｓ５３での判定に使用される閾値、言い換えると、図１４及び図１５に示す分割領域Ａ１〜Ａ５，Ｂ１〜Ｂ３の境界線Ｌ１〜Ｌ６を特定する情報等は、テストクラブのモデル毎に整理され、対応関係データ２８として記憶部２３内に格納されている。すなわち、対応関係データ２８とは、Ｐ1_AVE,Ｐ2_AVE及びＶhの大きさと最適スイングＭＩ帯との対応関係を定めるデータである。なお、境界線Ｌ１〜Ｌ４は互いに概ね平行であり、境界線Ｌ５，Ｌ６も互いに概ね平行である。また、境界線Ｌ１〜Ｌ６は、いずれもＰ1_AVE−Ｐ2_AVE平面において負の傾きを持つ直線である。ステップＳ４１〜Ｓ４４及びステップＳ５１〜Ｓ５３では、この記憶部２３内の対応関係データ２８が参照され、上記の判定が行われる。なお、図２では、対応関係データ２８は、フィッティングプログラム３とは別のデータとして示されているが、プログラム３内に組み込まれていてもよい。
【０１１５】
ところで、アベレージモデル領域では、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEを示す点が領域Ｂ３に属する場合には、ヘッド速度Ｖhが３５．５ｍ／ｓより大きいか否かにより、異なる最適スイングＭＩ帯に振り分けられるようになっている。これは、図１５に示す実験結果に示されるとおり、これらの領域においては、Ｖh＞３５．５ｍ／ｓか否かで最適スイングＭＩが異なる値になっているためである。
【０１１６】
ところで、図１４及び図１５には、上述したとおり、実際にプロモデルユーザ２１名及びアベレージモデルユーザ２２名にテストクラブを試打させた実験の結果がプロットされている。そして、プロモデルユーザに関しては、２１個のデータうち、Ｆのマークを付している１名を除く２０名において、最適スイングＭＩが上述の最適振り易さ決定工程で決定される最適スイングＭＩ帯に属する結果となった。すなわち、正答率９５％以上でフィッティングができていることが確認された。また、アベレージモデルユーザに関しては、２２個のデータうち、Ｆのマークを付している１名を除く２１名において、最適スイングＭＩが上述の最適振り易さ決定工程で決定される最適スイングＭＩ帯に属する結果となった。すなわち、正答率９５％以上でフィッティングができていることが確認された。
【０１１７】
＜１−２−７．第２指標算出工程＞
以下、計測工程で得られた計測データに基づいて、第２スイング指標を算出する第２指標算出工程（Ｓ７）について説明する。第２スイング指標とは、最適剛性指標を決定するための、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態では、第２スイング指標として、後述される第１〜第４特徴量Ｆ1〜Ｆ4が算出される。
【０１１８】
第１〜第４特徴量Ｆ1〜Ｆ4について理解するためには、まず、最適剛性指標について理解することが重要である。最適剛性指標とは、ゴルファー７に適したシャフト４０の剛性を示す指標のことであり、本実施形態では、シャフト４０の剛性は、シャフト４０の複数の位置における曲げ剛性の分布（以下、ＥＩ分布）として評価される。本実施形態に係るＥＩ分布は、定量的に数値を用いて表現され、より具体的には、インターナショナル・フレックス・コード（ＩＦＣ）を用いて算出される。そのため、まず、このＩＦＣについて説明する。なお、ＩＦＣは、本出願人により広く提案されているシャフトの特性を示す公知の指標であり、例えば、特許文献１をはじめとして、既に様々な文献で詳しく説明されている。従って、ここで改めて説明する必要は必ずしもないが、参考のため、ここでも説明を行う。
【０１１９】
ＩＦＣは、図１６に示すとおり、シャフト４０の延びる方向に沿った４つの位置Ｈ１〜Ｈ４におけるシャフト４０の曲げ剛性をそれぞれ０〜９の１桁の数値で表し、この４つの数値をシャフト４０の延びる方向に沿って配列したコードである。より具体的には、シャフト４０のバット端からチップ端に向かってこの順に概ね一定間隔で、４つの測定点Ｈ１〜Ｈ４が定義される。例えば、シャフト４０のチップ端から３６インチの箇所を測定点Ｈ１とし、２６インチの箇所を測定点Ｈ２とし、１６インチの箇所を測定点Ｈ３とし、６インチの箇所を測定点Ｈ４とすることができる。そして、これらの４つの測定点Ｈ１〜Ｈ４のそれぞれにおける曲げ剛性の値（以下、ＥＩ値）Ｊ1〜Ｊ4が計測される。
【０１２０】
シャフト４０の各測定点Ｈ（Ｈ１〜Ｈ４）におけるＥＩ値（Ｎ・ｍ2）は、様々な方法で測定することができ、例えば、インテスコ社製の２０２０型計測機（最大荷重５００ｋｇｆ）を用いて図１７に示すようにして測定することができる。この測定方法では、２つの支持点１１１，１１２でシャフト４０を下方から支持しつつ、測定点Ｈに上方から荷重Ｆを加えたときのたわみ量を測定する。支持点１１１と支持点１１２との間の距離（スパン）は、例えば、２００ｍｍとすることができ、測定点Ｈは、支持点１１１と支持点１１２の中間点とすることができる。より具体的には、支持点１１１，１１２を支える支持体１１４，１１５を固定した状態で、測定点Ｈにおいて圧子１１３を一定速度（例えば、５ｍｍ／分）で下方へ移動させる。そして、荷重Ｆが２０ｋｇｆに達した時点で圧子１１３の移動を終了させ、この瞬間のシャフト４０のたわみ量（ｍｍ）を測定し、このたわみ量をＥＩ値（Ｎ・ｍ2）に換算する。
【０１２１】
次に、以上の４つの測定点Ｈ１〜Ｈ４におけるＥＩ値Ｊ1〜Ｊ4を、それぞれ１０段階のランク値Ｋ1〜Ｋ4に変換する。具体的には、ランク値Ｋ1〜Ｋ4は、それぞれ測定点Ｈ１〜Ｈ４用の以下の変換表（表１〜表４）に従って、ＥＩ値Ｊ1〜Ｊ4から算出することができる（表１〜４中、変換後のランク値をＩＦＣの欄に示している）。そして、このようにして測定点Ｈ１〜Ｈ４にそれぞれ付与された４つのランク値Ｋ1〜Ｋ4を、よりバット側に対応する値がより左に、よりチップ側に対応する値がより右にくるように配列する。こうして得られた４桁のコードが、ＩＦＣである。ＩＦＣでは、各桁の数値が大きい程、対応する位置での剛性が高いことを意味する。
【表１】
2017170105.tif 000027
【表２】
2017170105.tif 000028
【表３】
2017170105.tif 000029
【表４】
2017170105.tif 000030
【０１２２】
第２指標算出工程では、算出部２４Ｄにより、第１〜第４特徴量Ｆ1〜Ｆ4が算出される。本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4は、それぞれゴルファー７に適したＥＩ値Ｊ1〜Ｊ4である最適ＥＩ値ＪS1〜ＪS4、ひいてはゴルファー７に適したランク値Ｋ1〜Ｋ4である最適ランク値ＫS1〜ＫS4を決定するための指標である。そのため、本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4としては、それぞれ最適ＥＩ値ＪS1〜ＪS4と相関を有する特徴量が選択される。また、本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4として以下の指標が用いられるが、第２スイング指標としては、最適剛性指標との相関が認められる限り、その他の任意の特徴量を用いることができる。
【０１２３】
第１特徴量Ｆ1は、トップ付近のコック方向の角速度ωyの傾きであり、例えばトップから５０ｍｓ前の角速度ωyと、トップから５０ｍｓ後の角速度ωyとの和で表すことができる。
【０１２４】
第２特徴量Ｆ2は、トップから、角速度ωyが最大となる時点までの当該角速度ωyの平均値である。第２特徴量Ｆ2は、まず、トップからインパクトまでの間で角速度ωyが最大となる時点を求め、トップからこの時点までの角速度ωyの累積値を、トップからこの時点までの時間で除することにより算出される。
【０１２５】
第３特徴量Ｆ3は、角速度ωyが最大となる時点からインパクトまでの当該角速度ωyの平均値である。第３特徴量Ｆ3は、角速度ωyが最大となる時点からインパクトまでの角速度ωyの累積値を、角速度ωyが最大となる時点からインパクトまでの時間で除することにより算出される。
【０１２６】
第４特徴量Ｆ4は、トップからインパクトまでの角速度ωyの平均値であり、トップからインパクトまでの角速度ωyの累積値を、トップからインパクトまでの時間で除することにより算出される。
【０１２７】
ところで、スイング動作中、ゴルフクラブ４のシャフト４０は、当該ゴルフクラブの先端に比較的重量が大きいヘッド４１が存在するため、その慣性により曲げが生じる。この曲げは、スイングの全過程において、シャフト４０の同一箇所に生じるのではなく、図１８に示されるように、トップからインパクトに向けてシャフト４０の手元側から先端側に伝わる。換言すれば、トップからインパクトに向けてスイングが進行するにしたがい、シャフト４０における曲げの位置が当該シャフト４０の手元側から先端側に移動する。
【０１２８】
より具体的には、アドレスからテイクバックを行い、トップに至った時点（図１８において（１）で示される時点）では、シャフト４０の手元付近に曲げが生じる。ついで、切り返しを行い、ダウンスイング初期（図１８において（２）で示される時点）に至ると、曲げはシャフト４０の先端側にやや移動する。さらに、ゴルファー７の腕が水平になる時点（図１８において（３）で示される時点）では、曲げはシャフト４０中央よりも先端側に移動する。そして、インパクト直前（図１８において（４）で示される時点）では、曲げはシャフト４０の先端付近まで移動する。
【０１２９】
従って、第１〜第４特徴量Ｆ1〜Ｆ4は、それぞれスイング動作中のトップ付近からインパクト付近までの間の第１〜第４区間において算出することができる。また、ここでの第１〜第３区間は、この順に時間経過に沿った区間であり、互いに一部重複する又は重複することのない区間となっている。
【０１３０】
＜１−２−８．最適剛性決定工程＞
次に、最適剛性決定工程（Ｓ８）の流れについて説明する。本工程では、決定部２４Ｅが、第２スイング指標（第１〜第４特徴量Ｆ1〜Ｆ4）と最適剛性指標（最適ＥＩ値ＪS1〜ＪS4）との相関関係を表す予め定められた近似式に従って、最適剛性指標（最適ＥＩ値ＪS1〜ＪS4）を決定する。本実施形態に係る近似式は、線形近似式であり、以下のように表される。
ＪS1＝ａ1・Ｆ1＋ｂ1
ＪS2＝ａ2・Ｆ2＋ｂ2
ＪS3＝ａ3・Ｆ3＋ｂ3
ＪS4＝ａ4・Ｆ4＋ｂ4
【０１３１】
決定部２４Ｅは、第２指標算出工程で算出された第１〜第４特徴量Ｆ1〜Ｆ4をかかる近似式に代入することにより、最適ＥＩ値ＪS1〜ＪS4を算出する。また、決定部２４Ｅは、上述の表１〜表４の換算表に従って、最適ＥＩ値ＪS1〜ＪS4をそれぞれ最適ランク値ＫS1〜ＫS4に変換する。
【０１３２】
なお、上式中、ａ１〜ａ４及びｂ１〜ｂ４は、予め行われた多数の実験結果から回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。ここでいう実験とは、例えば、特許文献１と同様に、以下のように行うことができる。すなわち、まず、多数のゴルファーの各々に複数のゴルフクラブを振ってもらい、そのときの飛距離、打球の方向性（左右ずれ）及び官能試験による振り易さを数値化する。そして、その数値から各ゴルファーに適しているゴルフクラブを決定し、当該ゴルフクラブのＥＩ値を当該ゴルファーの最適ＥＩ値とする。また、上記と同様の方法で、各ゴルファーの第１〜第４特徴量Ｆ1〜Ｆ4を算出する。そして、かかる実験の後、多数のゴルファー分の第１〜第４特徴量Ｆ1〜Ｆ4及び最適ＥＩ値のデータを回帰分析することで、ａ１〜ａ４及びｂ１〜ｂ４が算出される。
【０１３３】
また、より信頼性の高い近似式とするために、ａ１〜ａ４及びｂ１〜ｂ４の値を条件に応じて変更することができる。例えば、ヘッド速度Ｖhに応じて、近似式を用意することができる。一例としては、上記実験データを、ヘッド速度帯に応じて分類し（例えば、４５ｍ／ｓ以上、４１〜４５ｍ／ｓ、４１ｍ／ｓ以下）、同じ分類に属するデータのみを対象に上記近似式を作成し、ａ１〜ａ４及びｂ１〜ｂ４を決定しておくことができる。そして、最適剛性決定工程では、ゴルファー７のヘッド速度Ｖhがどのヘッド速度帯に属するかを判断し、当該ヘッド速度帯に対応する近似式を用いて、最適剛性指標を算出する。
【０１３４】
＜１−２−９．最適クラブ選択工程＞
以上の工程Ｓ１〜Ｓ８により、最適振り易さ指標（最適スイングＭＩ帯）及び最適剛性指標（最適ＥＩ値ＪS1〜ＪS4、最適ランク値ＫS1〜ＫS4）が決定されると、選択部２４Ｆは、最適クラブ選択工程（Ｓ９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、ゴルファー７に適したシャフト４０(以下、推奨シャフト)が特定される。また、フィッティングの結果として、ゴルファー７に適したゴルフクラブ４（以下、推奨ゴルフクラブ）も特定される。
【０１３５】
本実施形態では、まず、選択部２４Ｆは、推奨ゴルフクラブに用いられるべきヘッド４１(以下、推奨ヘッド)の種類を決定する。推奨ヘッドの種類の決定は、本明細書では説明されないフィッティング処理により行うこともできるし、表示部２１及び入力部２２を介してユーザに質問する等してお好みのヘッド４１を選択させることにより行うこともできる。そして、選択部２４Ｆは、推奨ヘッドのスペックを示す情報をヘッドＤＢ２７内から読み出すととともに、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。ヘッドＤＢ２７内に登録されているヘッド４１のスペックを示す情報には、製造メーカー、型番及び重量等が含まれる。一方、シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報には、製造メーカー、型番、４つの位置Ｈ１〜Ｈ４におけるＥＩ値Ｊ1〜Ｊ4及びランク値Ｋ1〜Ｋ4（ＩＦＣ）、シャフト４０の重量、フレックス、トルク、調子、長さ及び重心位置等が含まれる。選択部２４Ｆは、これらの情報から、各シャフト４０と推奨ヘッドとを組み合わせた場合のゴルフクラブ４のスイング慣性モーメントＩSを算出し、その値が最適スイングＭＩ帯に属することとなるゴルフクラブ４（以下、第１絞り込みゴルフクラブ）及びそれに含まれるシャフト４０（以下、第１絞り込みシャフト）を特定する。なお、第１絞り込みゴルフクラブ及びシャフトは、通常多数本存在する。
【０１３６】
続いて、選択部２４Ｆは、各第１絞り込みシャフトについて、当該シャフトのランク値Ｋ1〜Ｋ4と、最適剛性決定工程で決定された最適ランク値ＫS1〜ＫS4との一致度を算出し、推奨シャフトとして、一致度の高いシャフトを特定する。一致度は、例えば、以下の数２５の式に従って算出することができ、値が小さいほど一致度が高い。
【数２５】
2017170105.tif 000031
【０１３７】
なお、推奨シャフトとしては、１本のみ特定されてもよいし、複数本が特定されてもよい。また、推奨シャフトとしては、第１絞り込みシャフトの中で相対的に一致度が高い所定数のシャフトが特定されてもよいし、一定以上の一致度を有する全てのシャフトが特定されてもよい。表示制御部２４Ｇは、シャフトＤＢ２９を参照して、推奨シャフトの種類を示す情報とともに、推奨シャフトのスペックを示す情報（ＩＦＣの値を含む）を表示部２１上に表示させる。さらに、表示制御部２４Ｇは、ゴルファーの最適スイングＭＩ帯及び最適ランク値ＫS1〜ＫS4を表示部２１上に表示させる。また、表示制御部２４Ｇは、推奨ゴルフクラブが推奨シャフトと推奨ヘッドとを組み合わせたゴルフクラブであることを示す情報を、表示部２１上に表示させる。これにより、ユーザは、ゴルファー７に適したゴルフクラブ４及びシャフト４０の種類を知ることができるとともに、当該ゴルフクラブのＩＦＣの値、最適スイングＭＩ帯及び最適ランク値ＫS1〜ＫS4を知ることができる。
【０１３８】
＜２．第２実施形態＞
図１９に、第２実施形態に係るフィッティングシステム２００の全体構成を示す。フィッティングシステム２００は、第１実施形態に係るフィッティングシステム１００と多くの点で共通する。従って、以下では、簡単のため、第１実施形態と第２実施形態の差異点を中心に説明し、第１実施形態と同様の構成には同様の符号を付して、説明を省略する。
【０１３９】
図２０は、第２実施形態に係るフィッティング処理の流れを示すフローチャートである。図４と図２０とを比較すれば明らかなとおり、第１実施形態と第２実施形態の解析処理の相違点は、主として最適振り易さ決定工程（Ｓ６）に代えて、最適特性決定工程（Ｓ２０６）が実行され、最適剛性決定工程（Ｓ８）に代えて、最適剛性決定工程（Ｓ２０８）が実行され、最適クラブ選択工程（Ｓ９）に代えて、最適クラブ選択工程（Ｓ２０９）が実行される点にある。より具体的には、第１実施形態に係る最適振り易さ決定工程（Ｓ６）では、最終的に最適スイングＭＩ帯が決定されたが、第２実施形態に係る最適特性決定工程（Ｓ２０６）では、最終的にゴルファーに適したゴルフクラブ４のシャフト４０の重量（最適シャフト重量）が決定される。すなわち、第１実施形態では、最適振り易さ指標が決定されたが、第２実施形態では、ゴルファーに適したゴルフクラブ４の特定の部位の特性を表す最適特性指標が決定される。また、第１実施形態に係る最適剛性決定工程（Ｓ８）では、最終的に最適剛性指標として、４桁のＩＦＣのコード（最適ランク値ＫS1〜ＫS4）が決定されたが、第２実施形態に係る最適剛性決定工程（Ｓ２０８）では、最終的に最適剛性指標として、ゴルファーに適したゴルフクラブ４のシャフト４０のフレックス（最適フレックス）が決定される。また、第１実施形態に係る最適クラブ選択工程（Ｓ９）では、ヘッド４１を固定した上で、ゴルフクラブ４全体での最適振り易さ指標（最適スイングＭＩ帯）を達成することのできるシャフト４０を絞り込んだ後、最適剛性指標（最適ランク値ＫS1〜ＫS4）にできる限り合致するシャフト４０が決定されたが、第２実施形態に係る最適クラブ選択工程（Ｓ２０９）では、ヘッド４１を固定した上で、最適特性指標（最適シャフト重量）及び最適剛性指標（最適フレックス）にできる限り合致するシャフト４０が決定される。
【０１４０】
図１９に示すとおり、フィッティングシステム２００は、フィッティング装置２に代えて、フィッティング装置２０２を備える。フィッティング装置２０２は、フィッティング装置２と同様のハードウェア構成を有するが、フィッティング装置２０２では、フィッティングプログラム３に代えて、フィッティングプログラム２０３がインストールされている。そのため、制御部２４は、記憶部２３内のフィッティングプログラム２０３を読み出して実行することにより、仮想的に取得部２４Ａ、グリップ挙動導出部２４Ｂ、肩挙動導出部２４Ｃ、算出部２４Ｄとして動作する他、決定部２２４Ｅ、選択部２２４Ｆ及び表示制御部２２４Ｇとして動作することができる。決定部２２４Ｅは、第１実施形態との相違点である最適特性決定工程（Ｓ２０６）及び最適剛性決定工程（Ｓ２０８）を実行する仮想ユニットであり、選択部２２４Ｆ及び表示制御部２２５Ｇは、第１実施形態との相違点である最適クラブ選択工程（Ｓ２０９）を実行する仮想ユニットである。また、最適特性決定工程（Ｓ２０６）を実行できるよう、フィッティング装置２０２の記憶部２３には、対応関係データ２８に代えて、対応関係データ２２８が格納されている。対応関係データ２２８は、最適シャフト重量を決定するための条件を示すデータである。
【０１４１】
第２実施形態でも、第１実施形態と同じく、計測工程（Ｓ１）、第１変換工程（Ｓ２）、第２変換工程（Ｓ３）、肩挙動導出工程（Ｓ４）及び第１指標算出工程（Ｓ５）が順次実行され、続いて、最適特性決定工程（Ｓ２０６）が実行される。そしてその後、第２指標算出工程（Ｓ７）が実行され、続いて、最適剛性決定工程（Ｓ２０８）及び最適クラブ選択工程（Ｓ２０９）が順次実行される。以下では、第１実施形態との相違点である最適特性決定工程（Ｓ２０６）、最適剛性決定工程（Ｓ２０８）及び最適クラブ選択工程（Ｓ２０９）について、順に説明する。
【０１４２】
＜２−１．最適特性決定工程＞
最適特性決定工程（Ｓ２０６）では、決定部２２４Ｅは、第１スイング指標（腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖh）の大きさに応じて、最適シャフト重量の範囲（以下、最適シャフト重量帯）を決定する。このとき、Ｐ1_AVE，Ｐ2_AVE及びＶhの値が大きい程、最適シャフト重量帯は段階的に大きな値に設定される。具体的には、図２１に示すフローチャートに従って、最適特性決定工程は進行する。
【０１４３】
ステップＳ１０１では、決定部２２４Ｅは、ユーザの入力に基づいて、ユーザのお好みのヘッド４１（以下、嗜好ヘッド）の種類を決定する。嗜好ヘッドとは、推奨ゴルフクラブに含まれるべきヘッド４１である。嗜好ヘッドの種類の決定は、表示部２１及び入力部２２を介して、ユーザから嗜好ヘッドの種類を直接的に指定させることにより行うこともできるし、ユーザに各種質問を行うとともにその回答を受け付け、当該回答から所定のアルゴリズムに従って嗜好ヘッドを絞り込むことにより行うこともできる。
【０１４４】
続くステップＳ１０２では、決定部２２４Ｅは、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhの大きさに応じて、最適振り易さ指標である最適スイングＭＩを決定する。最適スイングＭＩは、第１実施形態と同様に、図１３に示すフローチャートに従って算出することも可能であるが、本実施形態では、以下の重回帰式に従って算出される。
（最適スイングＭＩ）＝ｅ1・Ｐ1_AVE＋ｅ2・Ｐ2_AVE＋ｅ3・Ｖh＋ｅ4
【０１４５】
上記重回帰式は、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhを説明変数とし、最適スイングＭＩを目的変数とする重回帰式である。係数ｅ1〜ｅ4は、予め実験により変数Ｐ1_AVE，Ｐ2_AVE,Ｖh及び最適スイングＭＩの多数のデータセットを取得し、これらのデータセットに基づいて重回帰分析を行うことにより算出され、対応関係データ２２８として記憶部２３内に格納されている。すなわち、対応関係データ２２８は、Ｐ1_AVE,Ｐ2v_AVE及びＶhの大きさと、最適スイングＭＩとの対応関係を定めるデータを含む。なお、上記重回帰式は、テストクラブのモデル毎に用意しておくこともできる。この場合、第１実施形態と同様に、解析の精度が向上し得る。また、図１９では、対応関係データ２２８は、フィッティングプログラム２０３とは別のデータとして示されているが、プログラム２０３内に組み込まれていてもよい。
【０１４６】
続くステップＳ１０３では、決定部２２４Ｅは、最適スイングＭＩの大きさ及び嗜好ヘッドの種類に応じて、最適シャフト重量帯を決定する。具体的には、以下の表５及び表６に従って、最適スイングＭＩを最適シャフト重量帯に変換する。表５は、ヘッドＡ用のデータであり、表６は、ヘッドＡとは種類の異なるヘッドＢ用のデータである。
【表５】
2017170105.tif 000032
【表６】
2017170105.tif 000033
【０１４７】
すなわち、決定部２２４Ｅは、嗜好ヘッドがヘッドＡである場合には、表５に最適スイングＭＩを照合し、最適シャフト重量帯を決定する。一方、嗜好ヘッドがヘッドＢである場合には、表６に最適スイングＭＩを照合し、最適シャフト重量帯を決定する。なお、本実施形態では、ステップＳ１０１において、嗜好ヘッドとしてヘッドＡ又はヘッドＢのいずれかが選択されるものとする。
【０１４８】
表５及び表６のデータは、最適スイングＭＩと最適シャフト重量帯との対応関係を定めるデータであり、対応関係データ２２８に含まれる。つまり、対応関係データ２２８とは、第１スイング指標（Ｐ1_AVE,Ｐ2_AVE及びＶh）と最適特性指標（最適シャフト重量帯）との対応関係を定めるデータであると言うこともできる。最適スイングＭＩと最適シャフト重量帯との対応関係を定める対応関係データ２２８は、ヘッド４１の種類毎に用意することができ、本実施形態では、同じモデルで重量の異なる２つのヘッド（ヘッドＡ及びヘッドＢ）に対しそれぞれ用意されている。ただし、最適スイングＭＩと最適シャフト重量帯との対応関係を定める対応関係データ２２８としては、ヘッド４１の種類に依存しないデータを用意することもできる。
【０１４９】
ところで、ゴルフクラブ４のスイング慣性モーメントＩSとは、ヘッド４１のスイング慣性モーメントと、シャフト４０のスイング慣性モーメントと、その他の部位（グリップ４２、フェラル等）のスイング慣性モーメントとの合計値である。従って、嗜好ヘッドとしてスイング慣性モーメントの小さい（又は大きい）ものが選択される場合には、その分だけ、推奨ゴルフクラブ４に含まれるべき推奨シャフトはスイング慣性モーメントの大きい（又は小さい）ものが選択されることが適切である。この点、本実施形態のように、ヘッド４１の種類毎に最適シャフト重量帯を決定するための対応関係データ２２８が用意される場合には、嗜好ヘッドに応じてシャフト４０が選択されることになり、ゴルフクラブ４全体でのスイング慣性モーメントＩSを最適化することができる。以上のとおり、本実施形態では、第１実施形態に引き続き、ヘッド４１の種類を決定した上で、ゴルフクラブ４全体での最適振り易さ指標を達成できるような最適シャフト重量帯が決定される。
【０１５０】
また、本発明者らは、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−７４５に様々な重量のシャフトを取り付け、これらのゴルフクラブを３７名のテスターに試打させた。そして、各テスターについて、飛距離を最大にするシャフトの重量を決定し、実験による最適シャフト重量とした。また、３７名のテスターの各々に対し、上述した最適特性決定工程と同じ方法で最適シャフト重量帯を導出した。そして、これらの結果を比較したところ、３７名中３３名において、実験による最適シャフト重量が最適振り易さ決定工程で決定された最適シャフト重量帯に属する結果となった。すなわち、約９０％の正答率でフィッティングができていることが確認された。
【０１５１】
＜２−２．最適剛性決定工程＞
次に、最適剛性決定工程（Ｓ２０８）の流れについて説明する。本工程では、決定部２２４Ｅが、第２スイング指標（第１〜第４特徴量Ｆ1〜Ｆ4）の大きさに応じて、最適剛性指標（最適フレックス）を決定する。
【０１５２】
具体的には、まず、決定部２２４Ｅは、第１〜第４特徴量Ｆ1〜Ｆ4に基づいて、最適ＥＩ値ＪS1〜ＪS4及び最適ランク値ＫS1〜ＫS4を決定する。これらの値ＪS1〜ＪS4，ＫS1〜ＫS4の決定方法は、第１実施形態と同様とすることができる。
【０１５３】
続いて、決定部２２４Ｅは、最適ランク値ＫS1〜ＫS4に基づいて、最適フレックスを決定する。フレックスとは、シャフト４０全体での硬さ（曲げ剛性）を評価する指標である。従って、シャフト４０の複数の位置におけるゴルファー７に適した曲げ剛性を表す最適ランク値ＫS1〜ＫS4が分かれば、これらの値に基づいて最適フレックスを算出することができる。例えば、特定の位置での最適ランク値を最適フレックスとすることもできるし、複数の位置での最適ランク値の平均値を最適フレックスをとすることもできる。
【０１５４】
なお、他の実施形態では、最適フレックスを最適ランク値ＫS1〜ＫS4ではなく、最適ＥＩ値ＪS1〜ＪS4から直接算出することもできるし、これらの値ＪS1〜ＪS4，ＫS1〜ＫS4に基づかずに算出することもきる。後者の場合、最適フレックスを特定可能な適切な特徴量を第２スイング指標決定工程で算出しておけばよい。
【０１５５】
＜２−３．最適クラブ選択工程＞
以上の工程により、最適特性指標（最適シャフト重量帯）及び最適剛性指標（最適フレックス）が決定されると、選択部２２４Ｆは、最適クラブ選択工程（Ｓ２０９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、推奨シャフトが特定される。また、推奨ゴルフクラブも特定される。
【０１５６】
本実施形態では、まず、選択部２２４Ｆは、シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報を参照して、最適剛性決定工程で決定された最適フレックスと同じフレックスを有するシャフト４０（以下、第１絞り込みシャフト）を特定する。第１絞り込みシャフトは、通常多数本存在する。続いて、選択部２２４Ｆは、第１絞り込みシャフトの中から、最適特性決定工程で決定された最適シャフト重量帯に属するシャフト４０を特定し、推奨シャフトとして選択する。このとき、第１絞り込みシャフトの中に最適シャフト重量帯に属するシャフト４０が存在しない場合には、第１絞り込みシャフトのうち最適シャフト重量帯に最も近いシャフト４０が、推奨シャフトとして選択される。推奨シャフトは、１本のみ特定されてもよいし、複数本が特定されてもよい。さらに、選択部２２４Ｆは、推奨シャフトと嗜好ヘッドを備えるゴルフクラブ４を、推奨ゴルフクラブとして選択する。なお、本実施形態では、最適シャフト重量帯よりも最適フレックスを優先して推奨シャフトが決定されるが、逆に、最適フレックスよりも最適シャフト重量帯を優先して推奨シャフトが決定されてもよい。
【０１５７】
表示制御部２２４Ｇは、シャフトＤＢ２９を参照して、推奨シャフトの種類を示す情報とともに、推奨シャフトのスペックを示す情報（ＩＦＣ、フレックス、トルク、調子及び重量の値を含む）を表示部２１上に表示させる。さらに、表示制御部２２４Ｇは、ゴルファーの最適スイングＭＩ帯及び最適ランク値ＫS1〜ＫS4を表示部２１上に表示させる。また、表示制御部２２４Ｇは、推奨ゴルフクラブが推奨シャフトと嗜好ヘッドとを組み合わせたゴルフクラブであることを示す情報を、表示部２１上に表示させる。これにより、ユーザは、推奨ゴルフクラブ及び推奨シャフトの詳細を知ることができる。
【０１５８】
＜３．変形例＞
以上、本発明の幾つかの実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。
【０１５９】
＜３−１＞
上記実施形態では、ゴルファー７のスイング動作を計測する計測機器として、加速度センサ、角速度センサ及び地磁気センサの３つを有するセンサユニット１が使用されたが、計測機器を他の構成とすることもできる。例えば、地磁気センサを省略することもできる。この場合には、統計的手法により、ｘｙｚ局所座標系からＸＹＺ全体座標系へと計測データを変換することが可能である。なお、このような手法については、公知技術であるため（要すれば、特開２０１３−５６０７４号公報参照）、ここでは詳細な説明を省略する。或いは、計測機器として、三次元計測カメラを使用することもできる。三次元計測カメラにより、ゴルファーやゴルフクラブ、ゴルフボールの挙動を計測する手法についても、公知であるため、ここでは詳細な説明を省略する。なお、三次元計測カメラを用いた場合には、計測データのｘｙｚ局所座標系からＸＹＺ全体座標系への変換工程を省略することもでき、直接的にＸＹＺ全体座標系でのグリップの挙動を計測することができる。
【０１６０】
＜３−２＞
上記実施形態では、最適振り易さ指標として、最適スイングＭＩを算出することが例示された。しかしながら、ゴルフクラブ４の振り易さを表すその他の様々な指標について、最適振り易さ指標を算出してもよい。例えば、ゴルフクラブ４の重量ｍ2（厳密には、質量をｍ2とするならば重量はｍ2ｇであるが、簡単のためここではいずれもｍ2で表すこととする）、グリップエンド周りの慣性モーメントＩG、及び重心周りの慣性モーメントＩ2も、ゴルフクラブ４の振り易さを表す指標であり、これらの指標に対して最適振り易さ指標を算出してもよい。また、複数の最適振り易さ指標を算出して、これら全ての最適振り易さ指標に基づいてフィッティングを行ってもよい。例えば、最適スイングＭＩ、ゴルファーに適したグリップエンド周りの慣性モーメントＩG（以下、最適グリップエンドＭＩ）、ゴルファーに適した重心周りの慣性モーメントＩ2（以下、最適重心ＭＩ）及びゴル
ファーに適したゴルフクラブ４の重量ｍ2（以下、最適クラブ重量）の全てを算出し、これら４つの条件に合致するゴルフクラブ４のシャフト４０をシャフトＤＢ２９内から検索してもよい。また、これら４つの条件に合致する最適特性指標を決定してもよい。
【０１６１】
なお、グリップエンド周りの慣性モーメントＩGとは、例えば、以下の式に従って算出
される。ＩGを決定するパラメータであるｍ2，Ｉ2，Ｌは、ゴルフクラブ４の諸元である
ため、ＩGもゴルフクラブ４の諸元となる。
ＩG＝Ｉ2＋ｍ2Ｌ2
【０１６２】
また、ゴルフクラブ４の重量ｍ2、グリップエンド周りの慣性モーメントＩG、及び重心周りの慣性モーメントＩ2は、スイング慣性モーメントＩSと概ね比例関係にある。例えば、図２２は、様々なゴルフクラブのＩS，ＩGの関係をプロットしたグラフであり、図２３は、様々なゴルフクラブのＩS，ｍ2の関係をプロットしたグラフである。従って、上記実施形態に係る最適振り易さ決定工程において、最適スイングＭＩを最適グリップエンドＭＩ、最適重心ＭＩ又は最適クラブ重量に変換すれば、ＩG、Ｉ2、ｍ2についての最適振り易さを決定することができる。例えば、図２４及び図２５は、図１３のフローチャートを改変したものであり、図２５は、最適クラブ重量を決定するためのフローチャートであり、図２４は、最適グリップエンドＭＩを決定するためのフローチャートである。
【０１６３】
＜３−３＞
上記実施形態では、最適振り易さ指標を決定するための第１スイング指標として、腕出力パワーＰ1_AVE、クラブ入力パワーＰ2_AVE及びヘッド速度Ｖhの３つが組み合わせて用いられた。しかしながら、この例に限られず、例えば、スイング指標としてＰ1_AVEだけ、Ｐ2_AVEだけを用いてもよいし、（Ｐ1_AVE，Ｖh）、（Ｐ2_AVE，Ｖh）、（Ｐ1_AVE，Ｐ2_AVE）のような２つの指標の組み合わせを用いてもよい。
【０１６４】
また、変形例３−２に記載のとおり、最適振り易さ指標としては、最適スイングＭＩに限らず様々な指標を設定することができるが、最適振り易さ指標との間に一定の関係（相関）が認められる限り、任意の指標を第１スイング指標とすることができる。例えば、スイング動作時にゴルファーにより発揮される腕エネルギーＥAVE、Ｅ1は、最適振り易さ指標と相関があり、第１スイング指標として用いることができる。また、以下の総肩トルクＴti及び平均肩トルクＴAVEも、最適振り易さ指標と相関があり、第１スイング指標とし
て用いることができる。Ｔti及びＴAVEは、スイング動作時にゴルファーにより発揮され
る肩周りのトルクを表す指標である。
【０１６５】
総肩トルクＴti及び平均肩トルクＴAVEは、以下のように算出することができる。まず、以下の式に従って、トップからインパクトまでの肩周りのトルクＴ1及びグリップ４２周りのトルクＴ2を算出する。
【数２６】
2017170105.tif 000034
【数２７】
2017170105.tif 000035
【０１６６】
総肩トルクＴtiは、上述の肩周りのトルクＴ1をトップからインパクトまでの区間で積分した値Ｔtiである。Ｔtiは、トップからインパクトまでの間にゴルファー７の肩周りで発揮されるトルクの総量である。一方、平均肩トルクＴAVEは、総肩トルクＴtiをトップからインパクトまでの時間で除算したスイング動作中の平均的な肩周りのトルクであり、ＴAVE＝Ｔti／（ｔi−ｔt）の式に従って算出される。なお、総肩トルクＴtiを算出するに当たり、正のトルクＴ1だけを積分してもよい。
【０１６７】
また、以下の指標も、第１スイング指標として用いることができる。なお、以下の（７）〜（９）以外の指標は、値が大きい程、最適振り易さ指標が大きくなる傾向があり、（７）〜（９）の指標は、値が小さい程、最適振り易さ指標が大きくなる傾向がある。
（１）トップの時刻でシャフト４０と全体座標系Ｚ軸（グリップよりも下方）とが為す角度θ３（図２６参照）
（２）スイング動作中の角速度ω2の平均値
（３）トップからインパクトまでの間での角速度ω2の最大値
（４）トップからインパクトまでのグリップ速度ＶGEの平均値
（５）トップからインパクトまでのグリップ速度ＶGEの最大値
（６）トップからインパクトまでの間でのグリップ４２の移動距離Ｄ
（７）コック解放タイミングｔrとインパクトの時間との差分（ここでいうコック解放タイミングｔrは、腕とシャフト４０が為すコック角θ４の解放スピードが速まり、腕のエネルギーがシャフト４０のエネルギーへと変わり始めるタイミングと定義することができる。）
（８）コック解放タイミングｔrでの腕とシャフト４０が為すコック角θ４（図２６参照）
（９）ダウンスイング時間、すなわち、トップからインパクトまでの時間
（１０）トップの時刻からグリップ４２周りのトルクＴ2が正負逆転する時刻までのトルクＴ2の積分値
【０１６８】
＜３−３−１＞
以下に、図２７〜図２９を参照しつつ、第１スイング指標として腕エネルギーＥAVE、平均肩トルクＴAVE及びヘッド速度Ｖhが算出され、これらの指標ＥAVE、ＴAVE、Ｖhに応じて最適クラブ重量の範囲（以下、最適重量帯）が算出される場合の、最適振り易さ決定工程について具体的に説明する。
【０１６９】
まず、ステップＳ１４０において、決定部２４Ｅは、計測工程で試打されたテストクラブの種類を判定する。プロモデルクラブが試打されていた場合には、ステップＳ１４１に進み、アベレージモデルクラブが試打されていた場合には、ステップＳ１５１に進む。いずれのテストクラブが試打されたかは、入力部２２を介してユーザにより入力される情報に基づいて判定されるものとする。
【０１７０】
続くステップＳ１４１及びこれに続くステップＳ１４２〜Ｓ１４５は、平均肩トルクＴAVE、腕エネルギーＥAVE及びヘッド速度Ｖhの大きさに応じて、最適重量帯を決定するステップである。ここでは、ＴAVE，ＥAVE及びＶhの値が大きい程、最適重量帯は段階的に大きな値に設定される。
【０１７１】
具体的には、ステップＳ１４１では、決定部２４Ｅは、ヘッド速度Ｖhが４５ｍ／ｓ以上であり、平均肩トルクＴAVEが８５Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが３５０Ｎ・ｍ／ｓ以上（以下、プロ条件１）であるか否かを判定する。そして、プロ条件１が満たされる場合には、３２０ｇ以上が最適重量帯であると判定する。一方、ステップＳ１４１でプロ条件１が満たされない場合には、ステップＳ１４２に進む。ステップＳ１４２では、決定部２４Ｅは、平均肩トルクＴAVEが７５Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが２５０Ｎ・ｍ／ｓ以上（以下、プロ条件２）であるか否かを判定する。そして、プロ条件２が満たされる場合には、３１０ｇ〜３２０ｇが最適重量帯であると判定する。一方、ステップＳ１４２でプロ条件２が満たされない場合には、ステップＳ１４３に進む。ステップＳ１４３では、決定部２４Ｅは、平均肩トルクＴAVEが７０Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１７５Ｎ・ｍ／ｓ以上（以下、プロ条件３）であるか否かを判定する。そして、プロ条件３が満たされる場合には、３０５ｇ〜３１５ｇが最適重量帯であると判定する。一方、ステップＳ１４３でプロ条件３が満たされない場合には、ステップＳ１４４に進む。ステップＳ１４４では、決定部２４Ｅは、平均肩トルクＴAVEが５５Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１２０Ｎ・ｍ／ｓ以上（以下、プロ条件４）であるか否かを判定する。そして、プロ条件４が満たされる場合には、３００ｇ〜３１０ｇが最適重量帯であると判定する。一方、ステップＳ１４４でプロ条件４が満たされない場合には、ステップＳ１４５に進む。ステップＳ１４５では、決定部２４Ｅは、ヘッド速度Ｖhが３８ｍ／ｓより大きいか否か（以下、プロ条件５）を判定する。そして、プロ条件５が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１４５でプロ条件５が満たされない場合には、プロモデルクラブよりも、アベレージモデルクラブの方が適切であると判定する。これを受けて、表示制御部２４Ｇは、アベレージモデルクラブを用いて計測工程からやり直す旨を示すメッセージを表示部２１上に表示させ、ステップＳ１５１に進む。
【０１７２】
一方、ステップＳ１５１及びこれに続くステップＳ１５２，Ｓ１５３は、平均肩トルクＴAVE、腕エネルギーＥAVE及びヘッド速度Ｖhの大きさに応じて、最適重量帯を決定するステップである。ここでは、ＴAVE，ＥAVE及びＶhの値が大きい程、最適重量帯は段階的に大きな値に設定される。
【０１７３】
具体的には、決定部２４Ｅは、ステップＳ１５１において、ヘッド速度Ｖhが４０ｍ／ｓ以上であり、平均肩トルクＴAVEが６５Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが２００Ｎ・ｍ／ｓ以上（以下、アベレージ条件１）であるか否かを判定する。そして、アベレージ条件１が満たされる場合には、アベレージモデルクラブよりも、プロモデルクラブの方が適切であると判定する。これを受けて、表示制御部２４Ｇは、プロモデルクラブを用いて計測工程からやり直す旨を示すメッセージを表示部２１上に表示させ、ステップＳ１４１に進む。一方、ステップＳ１５１でアベレージ条件１が満たされない場合には、ステップＳ１５２に進む。ステップＳ１５２では、決定部２４Ｅは、平均肩トルクＴAVEが５５Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１５０Ｎ・ｍ／ｓ以上（以下、アベレージ条件２）であるか否かを判定する。そして、アベレージ条件２が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１５２でアベレージ条件２が満たされない場合には、ステップＳ１５３に進む。ステップＳ１５３では、決定部２４Ｅは、ヘッド速度Ｖhが３８ｍ／ｓより大きく、平均肩トルクＴAVEが４８Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１００Ｎ・ｍ／ｓ以上（以下、アベレージ条件３）であるか否かを判定する。そして、アベレージ条件３が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１５３でアベレージ条件３が満たされない場合には、２９０ｇ以下が最適重量帯であると判定する。
【０１７４】
以上のステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３は、以下の知見に基づく。すなわち、本発明者らは、プロモデルのゴルフクラブを通常使用しているゴルファー（以下、プロモデルユーザ）１０名及びアベレージモデルのゴルフクラブを通常使用しているゴルファー（以下、アベレージモデルユーザ）１０名に、それぞれプロモデルクラブ及びアベレージモデルクラブを試打させたところ、図２８及び図２９に示す結果を得た。そして、本発明者らは、図２８及び図２９に示す結果から、第１スイング指標を示す空間は、プロモデル領域とアベレージモデル領域とに分割されることを発見した。図２８及び図２９の例では、プロモデル領域とは、平均肩トルクＴAVE≧５５Ｎ・ｍかつ腕エネ
ルギーＥAVE≧１５０Ｎ・ｍ／ｓの領域であり、アベレージモデル領域とは、平均肩トル
クＴAVE＜５５Ｎ・ｍ又は腕エネルギーＥAVE＜１５０Ｎ・ｍ／ｓの領域である。しかしながら、これらの数値は、テストクラブの種類等のスイング条件によって異なるものとなり得る。本実験では、プロモデルクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ  ＫＥＮＡ  Ｂｌｕｅ６  Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１５ｇ、バランスは、Ｄ２）のドライバーが用いられ、アベレージモデルクラブとして、ダンロップスポーツ株式会社製のＸＸＩＯ（登録商標）７（シャフトは、ＭＰ−７００  Ｒ−Ｆｌｅｘ、ゴルフクラブの重量は、２８５ｇ、バランスは、Ｄ１）のドライバーが用いられた。
【０１７５】
図２８及び図２９を参照して説明した上記実験では、ヘッド速度Ｖh及び、飛距離を最大にする最適クラブ重量も算出された。ヘッド速度Ｖhは、上記と同様の工程に従って算出された。一方、最適クラブ重量については、ゴルファーに様々な重量のゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブの重量を特定し、これを最適クラブ重量とした。より具体的には、プロモデルユーザには、３２４ｇ，３１６ｇ，３１２ｇ，３０８ｇ，３００ｇの５種類のゴルフクラブをスイングさせ、アベレージモデルユーザには、２９５ｇ，２９０ｇ，２８５ｇのゴルフクラブをスイングさせた。図２８及び図２９に、本実験により得られたヘッド速度Ｖh及び最適クラブ重量の値を示す。
【０１７６】
図２８及び図２９に示すとおり、上記の実験からは、最適クラブ重量は、腕エネルギーＥAVE及び平均肩トルクＴAVEが共に大きいほど、大きくなることが分かる。その結果、本発明者らは、平均肩トルクＴAVE−腕エネルギーＥAVE−ヘッド速度Ｖh空間を、プロモデルユーザについては図２８に示すように、アベレージモデルユーザについては図２９に示すように領域分割することで、最適重量帯を示す領域を定義可能であることを発見した。ただし、簡単のため、図２８及び図２９では、ヘッド速度Ｖhを表す軸は省略され、平均肩トルクＴAVE−腕エネルギーＥAVE平面が示されている。すなわち、上述のステップＳ１４１〜Ｓ１４５及びステップＳ５１〜Ｓ５３は、平均肩トルクＴAVE、腕エネルギーＥAVE及びヘッド速度Ｖhを示す点が、ＴAVE−ＥAVE−Ｖh空間におけるどの領域にプロットされるかに応じて、最適重量帯を判定するステップとなっている。なお、ステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３での判定に使用されるＴAVE，ＥAVE及びＶhの閾値、言い換えると、図２８及び図２９に示す分割領域の境界を示す値は、テストクラブのモデル毎に整理され、対応関係データ２８として記憶部２３内に格納されている。すなわち、本変形例に係る対応関係データ２８とは、ＴAVE，ＥAVE及びＶhの大きさと最適重量帯との対応関係を定めるデータである。ステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３では、この対応関係データ２８が参照され、上記の判定が行われる。
【０１７７】
＜３−３−２＞
以下に、図３０及び図３１を参照しつつ、第１スイング指標として腕エネルギーＥAVE、平均肩トルクＴAVE及びヘッド速度Ｖhが算出され、これらの指標ＥAVE、ＴAVE、Ｖhに応じて最適スイングＭＩ帯が算出される場合の、最適振り易さ決定工程について具体的に説明する。
【０１７８】
本変形例に先立つ計測工程では、プロモデルクラブ及びアベレージモデルクラブの２本のテストクラブに代えて、一本のセンサユニット１付きテストクラブがゴルファー７により試打される。ただし、本変形例でも、２本のテストクラブを打ち分けるようにし、フィッティングの精度を向上させることも可能である。
【０１７９】
具体的には、ステップＳ１６０では、決定部２４Ｅは、平均肩トルクＴAVEが７７Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１８０Ｎ・ｍ／ｓ以上（以下、条件１）であるか否かを判定する。条件１が満たされる場合には、決定部２４Ｅは、ヘッド速度Ｖhが４５ｍ／ｓ以上であるか否か（以下、条件２）を判定する（ステップＳ１６１）。そして、条件２が満たされる場合には、５６００ｋｇ・ｃｍ2以上が最適スイングＭＩ帯であると判定し、そうでなければ、５５９０〜５６３０ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ１６０で条件１が満たされない場合には、ステップＳ１６２に進む。ステップＳ１６２では、決定部２４Ｅは、平均肩トルクＴAVEが６０Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１７０Ｎ・ｍ／ｓ以上（以下、条件３）であるか否かを判定する。条件３が満たされる場合には、決定部２４Ｅは、ヘッド速度Ｖhが４５ｍ／ｓ以上であるか否か（以下、条件４）を判定する（ステップＳ１６３）。そして、条件４が満たされる場合には、５５９０〜５６３０ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定し、そうでなければ、５５１０〜５５９０ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定する。一方、ステップＳ１６２で条件３が満たされない場合には、ステップＳ１６４に進む。ステップＳ１６４では、決定部２４Ｅは、平均肩トルクＴAVEが５０Ｎ・ｍ以上であり、かつ、腕エネルギーＥAVEが１３０Ｎ・ｍ／ｓ以上（以下、条件５）であるか否かを判定する。そして、条件５が満たされる場合には、５４６０〜５５１０ｋｇ・ｃｍ2が最適スイングＭＩ帯であると判定し、そうでなければ、５４８０ｋｇ・ｃｍ2以下が最適スイングＭＩ帯であると判定する。
【０１８０】
以上のステップＳ１６０〜Ｓ１６４は、以下の知見に基づく。すなわち、本発明者らは、本実施形態に係るテストクラブをゴルファー２１名に試打させてスイング指標を算出したところ、図３１に示す結果を得た。なお、この実験で算出されたスイング指標も、平均肩トルクＴAVE及び腕エネルギーＥAVEであり、具体的な値は、上記と同様の工程に従って算出された。本実験では、テストクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ  Ｋｏｓｕｍａ  Ｂｌｕｅ６  Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１４ｇ、バランスは、Ｄ３）のドライバーが用いられた。また、本実験では、ヘッド速度Ｖh及び、飛距離を最大にする最適スイングＭＩも算出された。ヘッド速度Ｖhは、上記と同様の工程に従って算出された。一方、最適スイングＭＩについては、ゴルファーに様々なスイング慣性モーメントのゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブのスイング慣性モーメントを特定し、これを最適スイングＭＩとした。より具体的には、スイング慣性モーメントが５６５０ｋｋｇ・ｃｍ2、５６１０ｋｇ・ｃｍ2、５５５０ｋｇ・ｃｍ2、５４８５ｋｇ・ｃｍ2、５４００ｋｇ・ｃｍ2の５種類のゴルフクラブをスイングさせた。図３１には、本実験により得られたヘッド速度Ｖh及び最適スイングＭＩの値が示されている。
【０１８１】
図３１に示すとおり、上記の実験からは、最適スイングＭＩは、腕エネルギーＥAVE及び平均肩トルクＴAVEが共に大きいほど、大きくなることが分かる。その結果、本発明者らは、平均肩トルクＴAVE−腕エネルギーＥAVE−ヘッド速度Ｖh空間を、図３１に示すように領域分割することで、最適スイングＭＩ帯を示す領域を定義可能であることを発見した。ただし、簡単のため、図３１では、ヘッド速度Ｖhを表す軸は省略され、平均肩トルクＴAVE−腕エネルギーＥAVE平面が示されている。すなわち、上述のステップＳ１６０〜Ｓ１６４は、平均肩トルクＴAVE、腕エネルギーＥAVE及びヘッド速度Ｖhを示す点が、ＴAVE−ＥAVE−Ｖh空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定するステップとなっている。なお、ステップＳ１６０〜Ｓ１６４での判定に使用されるＴAVE，ＥAVE及びＶhの閾値、言い換えると、図３１に示す分割領域の境界を示す値は、対応関係データ２８として記憶部２３内に格納されている。すなわち、本変形例に係る対応関係データ２８とは、ＴAVE，ＥAVE及びＶhの大きさと最適スイングＭＩ帯との対応関係を定めるデータである。ステップＳ１６０〜Ｓ１６４では、この記憶部２３内の対応関係データ２８が参照され、上記の判定が行われる。
【０１８２】
また、上記のとおり、グリップエンド周りの慣性モーメントＩGは、スイング慣性モーメントＩSと概ね比例関係にある。従って、図３０の最適振り易さ決定工程で決定される最適スイングＭＩ帯を、ゴルファーに適した慣性モーメントＩGの範囲に変換することで、ＩGに基づく最適振り易さ決定工程を実行することができる。図３２は、図３０のフローチャートを改変したものであり、最適グリップエンドＭＩを決定するためのフローチャートである。なお、図３１の右側の凡例には、最適スイングＭＩを最適グリップエンドＭＩに換算した値も示されている。
【０１８３】
＜３−４＞
第１及び第２実施形態では、シャフトの剛性として、曲げ剛性が評価されたが、これに代えて、ねじれ剛性を評価してもよい。ねじれ剛性の値（以下、ＧＪ値）も、シャフト４０の延びる方向に沿った複数の位置において測定又は算出することができる。すなわち、シャフト４０の延びる方向に沿った複数の位置におけるねじれ剛性の分布を、シャフトの剛性としてもよい。この場合、最適剛性指標としては、ゴルファー７に適したＧＪ値（最適ＧＪ値）が決定されることになるが、最適ＧＪ値を決定するための第２スイング指標としては、最適ＧＪ値との相関が認められる任意の指標を用いることができる。このような第２スイング指標としては、例えば、特開２０１４−２１２８６２号公報に記載されているような、以下の指標を用いることができる。
【０１８４】
（１）グリップ角速度ωyが最大となるときからインパクトまでの単位時間あたりのグリップ角速度ωxの変化量の大きさ
（２）トップ付近でのグリップ角速度ωzの変化量
（３）トップからダウンスイング途中であってグリップ角速度ωyが最大となるときまでのグリップ角速度ωzの変化量の大きさ
【０１８５】
本変形例でも、第２スイング指標と最適ＧＪ値との関係を表す近似式を予め実験により算出し、記憶部２３内に格納しておくことで、計測工程で得られる計測データに基づく第２スイング指標から、最適ＧＪ値を決定することができる。
【０１８６】
また、第１実施形態において、最適剛性指標として、ゴルファー７に適したシャフト４０の複数の位置における剛性分布ではなく、ゴルファー７に適したシャフト４０のフレックス、調子又はトルクを決定するようにしてもよい。同様に、第２実施形態においても、最適剛性指標として、ゴルファー７に適したシャフト４０のフレックスではなく、ゴルファー７に適したシャフト４０の調子、トルク、又は複数の位置における剛性分布を決定するようにしてもよい。なお、トルクは、シャフト４０全体でのねじれ剛性を表す指標である。
【０１８７】
＜３−５＞
上記実施形態では、上述した（１）〜（５）を前提とした肩挙動導出工程の一例について説明したが、肩挙動導出工程のアルゴリズムはこれに限定されず、例えば、以下の態様とすることができる。
【０１８８】
本変形例に係る肩挙動導出工程では、上述の（１）〜（５）の前提のうち、（４）及び（５）が省略され、（１）〜（３）のみが前提とされる。具体的には、肩挙動導出部２４Ｃは、第２変換工程で得られるスイング平面Ｐ内でのグリップ４２の軌道を、円弧（円）に近似する（図３０参照）。そして、当該円弧（円）の中心を肩の位置Ｐs＝（ＰsX，ＰsY）とし、当該円弧（円）の中心からグリップ４２の軌道までの平均的な距離を、腕長さ（肩とグリップ４２との距離）Ｒとする。この場合、肩挙動導出処理に先行する第２変換工程においては、グリップ挙動導出部２４Ｂにより、スイング平面Ｐ内でのグリップ４２の軌道が算出されているものとする。例えば、グリップ４２の軌道は、スイング平面Ｐ内でのグリップ速度（ｖpY，ｖpZ）を積分し、グリップ４２の軌道上の点Ａi＝（Ｘi，Ｙi）（ｉ＝１，２，・・・）を算出することにより導出可能である。なお、ここでいうｉは、時系列に沿ったデータ番号である。
【０１８９】
以下に、グリップ４２の軌道の近似円（円弧）を導出する方法の一例を示す。まず、グリップ４２の軌道上の任意の３点、例えば、Ａi，Ａi+30，Ａi+60を考える。このとき、任意の三角形の外接円の中心は、当該三角形の三辺の垂直二等分線の交点となることから、Ａi，Ａi+30，Ａi+60を頂点とする三角形を考えたとき（図３１参照）、以下の数２４の式、ひいては数２５の式が成り立つ。
【数２８】
2017170105.tif 000036
【数２９】
2017170105.tif 000037
【０１９０】
そして、様々なｉに対する複数の数２５の式から以下の数２６の式を作成し、疑似逆行列を導出する。これにより、グリップ４２の軌道の近似円（円弧）の中心Ｐs＝（ＰsX，ＰsY）を導出することができる。
【数３０】
2017170105.tif 000038
【０１９１】
続いて、肩挙動導出部２４Ｃは、近似円（円弧）の中心Ｐs＝（ＰsX，ＰsY）からグリップ４２の軌道上の各点Ａi＝（Ｘi，Ｙi）までの距離の平均値を算出し、腕長さＲとする。そして、肩挙動導出部２４Ｃは、この腕長さＲに基づいて、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω1＝ＶGE／Ｒを算出する。
【０１９２】
図３２Ａは、上記実施形態に係る方法で再現したトップからインパクトまでの二重振り子モデルを示しており、図３２Ｂは、本変形例に係る方法で再現したトップからインパクトまでの二重振り子モデルを示している。これらの図を比較すると明らかであるが、本変形例では、前提（４）が省略されることにより、トップでの腕とゴルフクラブ４との為す角度が評価されるため、コック角の浅い／深いを評価することができる。また、前提（５）が省略されることにより、インパクト時の腕の角度が評価されるため、ハンドファスト／ハンドレイトを評価することができる。
【０１９３】
＜３−６＞
第２実施形態では、第１スイング指標に基づいて最適振り易さ指標（最適スイングＭＩ）を決定した後、この最適スイングＭＩに基づいて最適特性指標（最適シャフト重量）が決定された。しかしながら、第１スイング指標の大きさに応じて、最適特性指標を直接算出することもできる。具体的には、例えば、第１スイング指標を説明変数とし、最適特性指標を目的変数とする重回帰式を、対応関係データ２２８として予め算出し記憶しておくことができる。この場合、第１指標算出工程で導出される第１スイング指標を当該重回帰式に代入することにより、最適特性指標を導出することができる。あるいは、第１実施形態の図１４及び図１５と同様に、（好ましくは、ヘッド４１の種類毎に）Ｐ1_AVE−Ｐ2_AVE（−Ｖh）空間を領域分割して、最適特性指標に対応する領域を定義しておき（図３６参照）、当該定義情報を対応関係データ２２８として予め記憶しておくことができる。この場合、第１指標算出工程で導出される第１スイング指標がどの領域に属するかに応じて、最適特性指標を導出することができる。
【０１９４】
また、図３６の対応関係データ２２８に代えて、図３７Ａ〜図３７Ｃに示すような対応関係データ２２８を用いることができる。図３７Ａは、最適フレックスが「Ｘ」の場合に最適シャフト重量帯を決定するためのデータであり、図３７Ｂは、最適フレックスが「Ｓ」の場合に最適シャフト重量帯を決定するためのデータであり、図３７Ｃは、最適フレックスが「ＳＲ」の場合に最適シャフト重量帯を決定するためのデータである。すなわち、まず最適剛性指標を決定し、これに対応する対応関係データ２２８を選択する。そして、選択された対応関係データ２２８に基づき、第１スイング指標の大きさに応じて、最適特性指標を決定する。この方法によれば、最適振り易さ指標を決定することなく、第１スイング指標から最適特性指標を直接算出することができる他、最適特性指標よりも最適剛性指標を優先して、推奨シャフトを決定することができる。
【０１９５】
＜３−７＞
第２実施形態では、最適特性指標として最適シャフト重量が導出されたが、これに代えて又は加えて、ゴルファー７に適したシャフト４０以外の部位の特性を表す指標、例えば、ゴルファーに適したグリップ４２の重量を導出することもできる。また、ゴルフクラブ４の特定の部位の特性として、重量以外の特性、例えば、慣性モーメントを導出することもできる。
【０１９６】
＜３−８＞
上述した最適振り易さ指標及び最適特性指標は、ゴルファー７に適したゴルフクラブ４の全体又は特定の部位の特性を表す指標（以下、最適クラブ指標という）である。第１実施形態及び変形例３−６では、第１スイング指標が所定の領域Ａ１〜Ａ５，Ｂ１〜Ｂ３，Ｃ１〜Ｃ５のいずれに属するかに応じて最適クラブ指標が決定されたが、より精度の高いフィッティングのため、最適クラブ指標の決定時に以下の境界処理を実行することができる。本変形例に係る境界処理は、計測データに基づいて算出された第１スイング指標が、第１スイング指標を表す空間を分割することにより定義される複数の領域の境界付近に存在する場合には、既出の最適クラブ指標を修正する処理である。
【０１９７】
ここでは、一例として、変形例３−６と組み合わせる場合の境界処理について説明する。具体的な処理の流れは、図３８に示すとおりである。図３８に示す境界処理は、最適クラブ指標である最適シャフト重量が、図３６に示す判定基準（分割領域Ｃ１〜Ｃ５の境界線Ｌ６〜Ｌ１０）に従って決定された後に実行される。
【０１９８】
まず、ステップＳ６５において、決定部２２４Ｅは、判定指標が判定領域の上限境界付近に存在するか否かを判定する。なお、判定指標とは、計測データに基づいて算出されたＰ1_AVE，Ｐ2_AVE（第１スイング指標）を意味し、判定領域とは、Ｐ1_AVE−Ｐ2_AVE（−Ｖh）空間内に規定される領域Ｃ１〜Ｃ５のうち、判定指標が属する領域である。また、判定指標が上限境界付近に存在するか否かは、例えば、Ｐ1_AVE−Ｐ2_AVE（−Ｖh）空間内において、判定指標を表す点と判定領域の上限境界線との距離が閾値以内であるかによって判定することができる。判定指標が領域Ｃ２に属する場合には、境界線Ｌ７が上限境界線である。
【０１９９】
ステップＳ６５において判定指標が判定領域の上限境界付近に存在すると判定された場合には、決定部２２４Ｅは、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ６６）。
【０２００】
一方、ステップＳ６５において判定指標が判定領域の上限境界付近に存在しないと判定された場合には、処理はステップＳ６７に進む。ステップＳ６７では、判定指標が判定領域の下限境界付近に存在するか否かを判定する。なお、判定指標が下限境界付近に存在するか否かは、例えば、Ｐ1_AVE−Ｐ2_AVE（−Ｖh）空間内において、判定指標を表す点と判定領域の下限境界線との距離が閾値以内であるかによって判定することができる。判定指標が領域Ｃ２に属する場合には境界線Ｌ８が下限境界線である。
【０２０１】
ステップＳ６７において判定指標が判定領域の下限境界付近に存在すると判定された場合には、決定部２２４Ｅは、既出の最適シャフト重量を所定量だけ減量するように調整する（ステップＳ６８）。そして、以上のステップが終了すると、境界処理は終了する。
【０２０２】
＜３−９＞
変形例３−８では、最適クラブ指標をより精度よく決定するための境界処理について説明したが、本変形例に係る例外処理も、最適クラブ指標をより精度よく決定するための処理である。上記実施形態では、第１スイング指標を算出するに当たり、スイング平面Ｐ内での解析が行われたが、上述した第１スイング指標であるＰ1_AVE，Ｐ2_AVEや変形例３−３で示した指標は、スイング平面Ｐ内に投影された二次元的なゴルファー７のスイング動作を表す指標である。しかしながら、実際のスイング動作は、三次元的に行われる。本変形例に係る例外処理は、二次元的な解析による誤差を軽減し、最適クラブ指標の精度を向上させるための処理である。
【０２０３】
本変形例に係る例外処理では、第１スイング指標として、上述したものに加えて、スイング平面Ｐに現れないローテーションの動きを示す指標、及びコックの動きを表す指標（以下、コック指標という）が導出される。スイング平面Ｐに現れないローテーションの動きとは、スイング平面Ｐに投影されたゴルフクラブ４のシャフト軸周りの捻りのローテーションの動き（図３９Ａ参照）を表す指標（以下、捻り指標という）や、スイング平面Ｐから飛び出す方向（プッシュ方向）のゴルフクラブ４のローテーションの動き（図３９Ｂ参照）を表す指標（以下、プッシュ指標という）により評価することができる。ゴルファー７によっては、このようなローテーションの動きによりヘッド速度を獲得している。また、ゴルファー７によっては、ダウンスイング後半のコックの動きによりヘッド速度を獲得することもあるが、そのようなコックの動きの全体像は、必ずしもスイング平面Ｐ上に投影されない。ここで説明する例外処理では、このようなゴルファー７についても精度よく、最適クラブ指標を導出することができる。
【０２０４】
以下では、一例として、変形例３−６と組み合わせる場合の例外処理について説明する。具体的な処理の流れは、図４０に示すとおりである。図４０に示す例外処理は、最適クラブ指標である最適シャフト重量が、図３６に示す判定基準（分割領域Ｃ１〜Ｃ５の境界線Ｌ６〜Ｌ１０）に従って決定された後に実行される。
【０２０５】
まず、ステップ７１では、算出部２４Ｄが、計測データに基づいて、捻り指標を算出する。捻り指標は、例えば、シャフト軸周りの角速度ωzのトップからインパクトまでの平均値として算出することができる。続いて、決定部２２４Ｅが、ヘッド速度Ｖhが所定範囲内であり、かつ、捩り指標が所定値以上であるかを判定する。そして、かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７２）。一方、かかる条件が満たされない場合には、処理はステップＳ７３に進む。
【０２０６】
ステップＳ７３では、算出部２４Ｄが、計測データに基づいて、プッシュ指標を算出する。続いて、決定部２２４Ｅが、ヘッド速度Ｖhが所定範囲内であり、プッシュ指標が所定値以上であり、かつ、既出の最適シャフト重量が所定範囲内であるか否かを判定する。かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７４）。一方、かかる条件が満たされない場合には、処理はステップＳ７５に進む。なお、プッシュ指標が所定値以上の場合とは、例えば、トップのωx≧０、かつ、トップからインパクトまでのωyの平均値／トップからインパクトまでのωxの平均値≧１．５が満たされる場合とすることができる。
【０２０７】
ステップＳ７５では、算出部２４Ｄが、計測データに基づいて、コック指標を算出する。コック指標は、例えば、トップ付近のωyの差と、ωyが最大となる時刻からインパクトまでのωyの平均値とを比較した値として算出することができる。続いて、決定部２２４Ｅが、ヘッド速度Ｖhが所定範囲内であり、コック指標が所定値以上であり、かつ、既出の最適シャフト重量が所定範囲内であるか否かを判定する。かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７６）。
【０２０８】
以上のステップが完了すると、例外処理は終了する。本例外処理と、変形例３−８の境界処理とは、いずれも最適クラブ指標をより精度よく決定するための処理であり、両処理は、組み合わせることもできる。この場合、境界処理の後、本例外処理を実行することが好ましい。
【実施例】
【０２０９】
以下、本発明の実施例について説明する。但し、本発明は、以下の実施例に限定されない。
【０２１０】
＜実施例１＞
スイング慣性モーメントＩSが５６１５ｋｇ・ｃｍ2のゴルフクラブを普段使用しているゴルファーに対し、特許文献１のフィッティング方法（比較例）を適用するとともに、第１実施形態に係るフィッティング方法（実施例１）を適用した。
【０２１１】
比較例の方法では、マイクラブ重量±５ｇの範囲に収まるようなゴルフクラブしか選択されないため、Miyazaki Indio 7S(ＩS＝約５６１０ｋｇ・ｃｍ2)のシャフトが選択された。そして、当該シャフトを備えるゴルフクラブを用いて、当該ゴルファーに５球の試打を行わせ、打球の飛距離及び左右ぶれの平均値を計測した。
【０２１２】
一方、実施例１の方法では、Miyazaki Indio 6S(ＩS＝約５５４０ｋｇ・ｃｍ2)のシャフトが選択された。そして、当該シャフトを備えるゴルフクラブであって、比較例の場合とはシャフトのみが異なる（すなわち、ヘッドは同一の）ゴルフクラブを用いて、当該ゴルファーに５球の試打を行わせ、打球の飛距離及び左右ぶれの平均値を計測した。
【０２１３】
図４１Ａ及び図４１Ｂは、実施例１及び比較例による飛距離及び左右ぶれの平均値及びばらつきを示すグラフである。同図に示すように、飛距離及び左右ぶれの両方について、比較例の方法で選定されたゴルフクラブよりも実施例１の方法で選定されたゴルフクラブを用いた方が、好ましい結果が得られた。従って、実施例１の優位性が確認された。
【０２１４】
＜実施例２＞
第２実施形態に係るフィッティング方法（実施例２）を適用して、あるゴルファーに適切なゴルフクラブを選択したところ、ヘッドＡ（重量：２００ｇ）にシャフトＡ（フレックス：ＳＲ、重量：６５．５ｇ）を組み合わせたゴルフクラブＡが選択された。当該ゴルファーは、ヘッドＡと同等の重量のゴルフクラブを普段使用していたため、ヘッドＡが嗜好クラブされた。また、ゴルフクラブＡと重量が同等であり、振り易さ指標が同等であるゴルフクラブとして、ヘッドＢ（重量：２０２ｇ）にシャフトＢ（フレックス：ＳＲ、重量：５６．５ｇ）を組み合わせたゴルフクラブＢを用意した。そして、当該ゴルファーにゴルフクラブＡ及びＢの両方を５回ずつ試打させたところ、図４２Ａ及び図４２Ｂに示す結果が得られた。図４２Ａは、ゴルフクラブＡの試打結果であり、図４２Ｂは、ゴルフクラブＢの試打結果である。
【０２１５】
図４２Ａ及び図４２Ｂを比較すれば明らかなとおり、普段使用しているヘッドと同等の重量のヘッドＡを含むゴルフクラブＡをスイングした方が、左右ずれ及び飛距離の両方において優れていることが分かる。すなわち、ヘッドの重量を調整して最適振り易さ指標を実現するよりも、ヘッドをゴルファーの嗜好に合わせた上でシャフトの重量を調整して最適振り易さ指標を実現する方が適切であることが分かる。従って、ヘッドを固定した上で、最適振り易さ指標を達成できるようなシャフトを選択する上記実施形態に係るフィッティング方法の優位性が確認された。
【符号の説明】
【０２１６】
１  センサユニット（計測装置）
２  フィッティング装置
３  フィッティングプログラム
４  ゴルフクラブ
７  ゴルファー
２４Ａ  取得部
２４Ｂ  グリップ挙動導出部
２４Ｃ  肩挙動導出部
２４Ｄ  算出部
２４Ｅ  決定部
２４Ｆ  選択部
４０  シャフト
４１  ヘッド
４２  グリップ
【図１】
2017170105.tif 000039
【図２】
2017170105.tif 000040
【図３】
2017170105.tif 000041
【図４】
2017170105.tif 000042
【図５】
2017170105.tif 000043
【図６】
2017170105.tif 000044
【図７】
2017170105.tif 000045
【図８】
2017170105.tif 000046
【図９】
2017170105.tif 000047
【図１０】
2017170105.tif 000048
【図１１】
2017170105.tif 000049
【図１２】
2017170105.tif 000050
【図１３】
2017170105.tif 000051
【図１４】
2017170105.tif 000052
【図１５】
2017170105.tif 000053
【図１６】
2017170105.tif 000054
【図１７】
2017170105.tif 000055
【図１８】
2017170105.tif 000056
【図１９】
2017170105.tif 000057
【図２０】
2017170105.tif 000058
【図２１】
2017170105.tif 000059
【図２２】
2017170105.tif 000060
【図２３】
2017170105.tif 000061
【図２４】
2017170105.tif 000062
【図２５】
2017170105.tif 000063
【図２６】
2017170105.tif 000064
【図２７】
2017170105.tif 000065
【図２８】
2017170105.tif 000066
【図２９】
2017170105.tif 000067
【図３０】
2017170105.tif 000068
【図３１】
2017170105.tif 000069
【図３２】
2017170105.tif 000070
【図３３】
2017170105.tif 000071
【図３４】
2017170105.tif 000072
【図３５Ａ】
2017170105.tif 000073
【図３５Ｂ】
2017170105.tif 000074
【図３６】
2017170105.tif 000075
【図３７Ａ】
2017170105.tif 000076
【図３７Ｂ】
2017170105.tif 000077
【図３７Ｃ】
2017170105.tif 000078
【図３８】
2017170105.tif 000079
【図３９Ａ】
2017170105.tif 000080
【図３９Ｂ】
2017170105.tif 000081
【図４０】
2017170105.tif 000082
【図４１Ａ】
2017170105.tif 000083
【図４１Ｂ】
2017170105.tif 000084
【図４２Ａ】
2017170105.tif 000085
【図４２Ｂ】
2017170105.tif 000086