特許 6911298 ゴルフクラブのフィッティング装置、方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6911298

(24)【登録日】2021年7月12日

(45)【発行日】2021年7月28日

(54)【発明の名称】ゴルフクラブのフィッティング装置、方法及びプログラム

(51)【国際特許分類】

   A63B 69/36 20060101AFI20210715BHJP

【ＦＩ】

   A63B69/36 541S

   A63B69/36 541P

   A63B69/36 541W

【請求項の数】25

【全頁数】62

(21)【出願番号】特願2016-158230(P2016-158230)

(22)【出願日】2016年8月10日

(65)【公開番号】特開2017-170105(P2017-170105A)

(43)【公開日】2017年9月28日

【審査請求日】2019年6月20日

(31)【優先権主張番号】特願2015-174207(P2015-174207)

(32)【優先日】2015年9月3日

(33)【優先権主張国】JP

(31)【優先権主張番号】特願2016-57502(P2016-57502)

(32)【優先日】2016年3月22日

(33)【優先権主張国】JP

(73)【特許権者】

【識別番号】000183233

【氏名又は名称】住友ゴム工業株式会社

(74)【代理人】

【識別番号】100124039

【弁理士】

【氏名又は名称】立花 顕治

(74)【代理人】

【識別番号】100156845

【弁理士】

【氏名又は名称】山田 威一郎

(74)【代理人】

【識別番号】100179213

【弁理士】

【氏名又は名称】山下 未知子

(72)【発明者】

【氏名】岡崎 弘祐

(72)【発明者】

【氏名】植田 勝彦

【審査官】 佐々木 崇

(56)【参考文献】

【文献】 国際公開第２０１４／１３２８８５（ＷＯ，Ａ１）

【文献】 特開２０１１−０３０７６１（ＪＰ，Ａ）

【文献】 特許第４１４５６１８（ＪＰ，Ｂ２）

【文献】 特開２００５−０２１３２９（ＪＰ，Ａ）

【文献】 特開２０１４−０７３３１３（ＪＰ，Ａ）

【文献】 特開２０１５−１３４１０３（ＪＰ，Ａ）

【文献】 特開２０１３−２０８３６６（ＪＰ，Ａ）

【文献】 特開２０１３−１６５８０８（ＪＰ，Ａ）

【文献】 特開２０１２−０１６５８２（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６３Ｂ６９／００−６９／４０

(57)【特許請求の範囲】

【請求項1】

ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定する決定部と、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備え、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
フィッティング装置。

【請求項2】

前記算出部は、複数の種類の前記第１スイング指標を算出し、
前記決定部は、前記複数の種類の第１スイング指標の大きさに応じて、前記最適振り易さ指標を決定する、
請求項１に記載のフィッティング装置。

【請求項3】

前記第１スイング指標には、前記スイング動作時のヘッドスピードがさらに含まれる、
請求項１又は２に記載のフィッティング装置。

【請求項4】

前記決定部は、前記第１スイング指標が大きい又は小さい程、前記最適振り易さ指標を大きな値に決定する、
請求項１から３のいずれかに記載のフィッティング装置。

【請求項5】

前記テストクラブの種類毎に、前記第１スイング指標の大きさと前記最適振り易さ指標の大きさとの対応関係を定める対応関係データを記憶する記憶部
をさらに備え、
前記決定部は、前記テストクラブの種類に応じて、前記記憶部内の前記対応関係データを参照することにより、前記最適振り易さ指標を決定する、
請求項１から４のいずれかに記載のフィッティング装置。

【請求項6】

前記計測値には、前記テストクラブのグリップエンドにおける加速度、角速度及び地磁気の少なくとも１つが含まれる、
請求項１から５のいずれかに記載のフィッティング装置。

【請求項7】

前記決定部は、前記最適剛性指標として、前記ゴルファーに適した、シャフトの複数の位置における剛性分布を決定する、
請求項１から６のいずれかに記載のフィッティング装置。

【請求項8】

前記算出部は、前記第２スイング指標として、第１〜第３特徴量を算出し、
前記決定部は、前記第１〜第３特徴量のそれぞれの大きさに応じて、前記最適剛性指標として、前記シャフトのバット端からチップ端に向かってこの順に並ぶ第１〜第３位置における、前記ゴルファーに適したシャフトの剛性を示す第１〜第３剛性値を決定する、
請求項７に記載のフィッティング装置。

【請求項9】

前記決定部は、前記第２スイング指標と前記最適剛性指標との関係を表す予め定められた式に従って、前記最適剛性指標を決定する、
請求項７又は８に記載のフィッティング装置。

【請求項10】

ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定するステップと、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含み、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
フィッティング方法。

【請求項11】

ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定するステップと、
前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させ、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
フィッティングプログラム。

【請求項12】

ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定する決定部と、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備え、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である、
フィッティング装置。

【請求項13】

前記決定部は、前記第１スイング指標の大きさに加え、前記ゴルフクラブに含まれるべきヘッドの種類に応じて、前記最適特性指標を決定する、
請求項１２に記載のフィッティング装置。

【請求項14】

前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさに応じて、前記最適特性指標を決定し、
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
請求項１２又は１３に記載のフィッティング装置。

【請求項15】

前記決定部は、前記第１スイング指標を説明変数とし、前記最適振り易さ指標を目的変数する所定の回帰式に、前記算出部により算出された前記第１スイング指標を代入することにより、前記最適振り易さ指標を算出する、
請求項１４に記載のフィッティング装置。

【請求項16】

前記決定部は、前記第１スイング指標を説明変数とし、前記最適特性指標を目的変数する所定の回帰式に、前記算出部により算出された前記第１スイング指標を代入することにより、前記最適特性指標を算出する、
請求項１２又は１３に記載のフィッティング装置。

【請求項17】

ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
ユーザから前記ゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付ける入力受付部と、
前記ゴルファーのスイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量を決定する決定部と
を備え、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時にテストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記最適シャフト重量を決定し、
前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる、
フィッティング装置。

【請求項18】

ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正する決定部と、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備え、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる、
フィッティング装置。

【請求項19】

ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する算出部と、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定する決定部と、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部と
を備え、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標であって、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれるとともに、前記スイング動作時のヘッドスピード及びコックの動きを表す指標がさらに含まれ、
前記決定部は、前記ゴルファーの腕が出力するパワーを示す前記指標、前記テストクラブに入力されるパワーを示す前記指標、前記ゴルファーにより発揮されるエネルギーを示す前記指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す前記指標の少なくとも１つの大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、決定された当該最適クラブ指標が所定の閾値以下である場合に、前記ヘッドスピードが所定の閾値以上であり、かつ、前記コックの動きを表す指標が所定の閾値以上であるときに、前記決定された最適クラブ指標を修正する、
フィッティング装置。

【請求項20】

ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定するステップと、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含み、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である、
フィッティング方法。

【請求項21】

ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、シャフトの剛性を示す剛性指標であって、前記ゴルファーに適した前記剛性指標である最適剛性指標を決定するステップと、
前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させ、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれ、
前記第２スイング指標は、前記最適剛性指標との相関が認められる、前記ゴルファーによる前記スイング動作を特徴付ける特徴量であり、
前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である、
フィッティングプログラム。

【請求項22】

ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含み、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる、
フィッティング方法。

【請求項23】

ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させ、
前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる、
フィッティングプログラム。

【請求項24】

ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
を含み、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標であって、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれるとともに、前記スイング動作時のヘッドスピード及びコックの動きを表す指標がさらに含まれ、
前記最適クラブ指標を決定するステップは、前記ゴルファーの腕が出力するパワーを示す前記指標、前記テストクラブに入力されるパワーを示す前記指標、前記ゴルファーにより発揮されるエネルギーを示す前記指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す前記指標の少なくとも１つの大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、決定された当該最適クラブ指標が所定の閾値以下である場合に、前記ヘッドスピードが所定の閾値以上であり、かつ、前記コックの動きを表す指標が所定の閾値以上であるときに、前記決定された最適クラブ指標を修正するステップを含む、
フィッティング方法。

【請求項25】

ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、
前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップと、
前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップと、
前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップと
をコンピュータに実行させ、
前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標であって、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれるとともに、前記スイング動作時のヘッドスピード及びコックの動きを表す指標がさらに含まれ、
前記最適クラブ指標を決定するステップは、前記ゴルファーの腕が出力するパワーを示す前記指標、前記テストクラブに入力されるパワーを示す前記指標、前記ゴルファーにより発揮されるエネルギーを示す前記指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す前記指標の少なくとも１つの大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、決定された当該最適クラブ指標が所定の閾値以下である場合に、前記ヘッドスピードが所定の閾値以上であり、かつ、前記コックの動きを表す指標が所定の閾値以上であるときに、前記決定された最適クラブ指標を修正するステップを含む、
フィッティングプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、ゴルファーに適したゴルフクラブ、特にシャフトを選定するためのフィッティング装置、方法及びプログラムに関する。

【背景技術】

【0002】

従来より、ゴルファーにテストクラブを試打させてその動作を計測機器により計測し、当該計測値に基づいて当該ゴルファーに適したゴルフクラブを選定する様々なフィッティング方法が提案されている。その１つとして、特許文献１には、ゴルファーに適したゴルフクラブのシャフトを選定するためのフィッティング方法が開示されている。具体的には、特許文献１では、まずテストクラブによる計測値に基づいて、ゴルファーに適したシャフトの曲げ剛性が決定される。そして、データベースに登録されている多数のシャフトの中から、当該曲げ剛性の値に合致するシャフトが抽出される。かかる方法は、シャフトのフィッティングの精度を向上させる技術として、大いに期待されるものである。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１３−２２６３７５号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

ところで、特許文献１の方法において、データベース内には、同じような曲げ剛性を有するシャフトが多数登録されていることも十分に想定される。そうすると、ゴルファーに適した曲げ剛性を決定するだけでは、ゴルファーに適したシャフトを十分に絞り込むことができない場合がある。特に、同じような曲げ剛性を有するが、重量帯の異なるシャフトが多数存在する場合が多々あり、どの重量帯のシャフトを選定してよいのか判断に困ることがある。

【0005】

この点、特許文献１では、上述した曲げ剛性に加え、ゴルファーが普段使用しているゴルフクラブの重量（以下、マイクラブ重量）を基準に、シャフトの選定を行っている。具体的には、新たに提案されるゴルフクラブの重量がマイクラブ重量±５ｇとなるように、シャフトの選定を行っている。すなわち、新たに使用するゴルフクラブの重量がマイクラブ重量から大きく変化すると、タイミングがとりにくくなったり、振りにくく感じたりし、ゴルファーのパフォーマンスが低下する虞がある。この点、マイクラブ重量からの重量の変化が制限される特許文献１の方法は優れていると言える。

【0006】

しかしながら、マイクラブ重量がそもそも当該ゴルファーに全く適していない、すなわち、オーバースペックだったり、アンダースペックだったりすることも考えられる。このような場合には、特許文献１の方法では、フィーリングが改善されても弾道が改善されなかったり、フィーリングが合わなかったりといった事態が生じ得る。

【0007】

本発明は、ゴルファーに適したゴルフクラブを精度よく選定するためのフィッティング装置、方法及びプログラムを提供することを目的とする。

【課題を解決するための手段】

【0008】

第１観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部とを備える。

【0009】

第２観点に係るフィッティング装置は、第１観点に係るフィッティング装置であって、前記算出部は、複数の種類の前記第１スイング指標を算出する。前記決定部は、前記複数の種類の第１スイング指標の大きさに応じて、前記最適振り易さ指標を決定する。

【0010】

第３観点に係るフィッティング装置は、第１観点又は第２観点に係るフィッティング装置であって、前記第１スイング指標には、前記スイング動作時に前記ゴルファーの腕が出力するパワーを示す指標、前記スイング動作時に前記テストクラブに入力されるパワーを示す指標、前記スイング動作時に前記ゴルファーにより発揮されるエネルギーを示す指標、及び、前記スイング動作時に前記ゴルファーにより発揮されるトルクを示す指標の少なくとも１つが含まれる。

【0011】

第４観点に係るフィッティング装置は、第１観点から第３観点のいずれかに係るフィッティング装置であって、前記第１スイング指標には、前記スイング動作時のヘッドスピードが含まれる。

【0012】

第５観点に係るフィッティング装置は、第１観点から第４観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標が大きい又は小さい程、前記最適振り易さ指標を大きな値に決定する。

【0013】

第６観点に係るフィッティング装置は、第１観点から第５観点のいずれかに係るフィッティング装置であって、前記テストクラブの種類毎に、前記第１スイング指標の大きさと前記最適振り易さ指標の大きさとの対応関係を定める対応関係データを記憶する記憶部をさらに備える。前記決定部は、前記テストクラブの種類に応じて、前記記憶部内の前記対応関係データを参照することにより、前記最適振り易さ指標を決定する。

【0014】

第７観点に係るフィッティング装置は、第１観点から第６観点のいずれかに係るフィッティング装置であって、前記振り易さ指標には、前記ゴルファーの肩周りの慣性モーメント、前記ゴルフクラブの慣性モーメント、及び前記ゴルフクラブの重量の少なくとも１つが含まれる。

【0015】

第８観点に係るフィッティング装置は、第１観点から第７観点のいずれかに係るフィッティング装置であって、前記計測値には、前記テストクラブのグリップエンドにおける加速度、角速度及び地磁気の少なくとも１つが含まれる。

【0016】

第９観点に係るフィッティング装置は、第１観点から第８観点のいずれかに係るフィッティング装置であって、前記決定部は、前記最適剛性指標として、前記ゴルファーに適した、シャフトの複数の位置における剛性分布を決定する。

【0017】

第１０観点に係るフィッティング装置は、第９観点に係るフィッティング装置であって、前記算出部は、前記第２スイング指標として、第１〜第３特徴量を算出する。前記決定部は、前記第１〜第３特徴量のそれぞれの大きさに応じて、前記最適剛性指標として、前記シャフトのバット端からチップ端に向かってこの順に並ぶ第１〜第３位置における、前記ゴルファーに適したシャフトの剛性を示す第１〜第３剛性値を決定する。

【0018】

第１１観点に係るフィッティング装置は、第９観点又は第１０観点に係るフィッティング装置であって、前記決定部は、前記第２スイング指標と前記最適剛性指標との関係を表す予め定められた式に従って、前記最適剛性指標を決定する。

【0019】

第１２観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0020】

第１３観点に係るフィッティングプログラムは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適振り易さ指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0021】

第１４観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出する算出部と、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定する決定部と、前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する選択部とを備える。

【0022】

第１５観点に係るフィッティング装置は、第１４観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに加え、前記ゴルフクラブに含まれるべきヘッドの種類に応じて、前記最適特性指標を決定する。前記特定の部位は、前記ヘッド以外の部位である。

【0023】

第１６観点に係るフィッティング装置は、第１４観点又は第１５観点に係るフィッティング装置であって、前記特定の部位は、シャフト又はグリップである。

【0024】

第１７観点に係るフィッティング装置は、第１４観点から第１６観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさに応じて、前記最適特性指標を決定する。

【0025】

第１８観点に係るフィッティング装置は、第１７観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標を説明変数とし、前記最適振り易さ指標を目的変数する所定の回帰式に、前記算出部により算出された前記第１スイング指標を代入することにより、前記最適振り易さ指標を算出する。

【0026】

第１９観点に係るフィッティング装置は、第１４観点から第１６観点のいずれかに係るフィッティング装置であって、前記決定部は、前記第１スイング指標を説明変数とし、前記最適特性指標を目的変数する所定の回帰式に、前記算出部により算出された前記スイング指標を代入することにより、前記最適特性指標を算出する。

【0027】

第２０観点に係るフィッティング装置は、第１４観点から第１９観点のいずれかに係るフィッティング装置であって、前記最適特性指標は、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量である。

【0028】

第２１観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、ユーザから前記ゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付ける入力受付部と、前記ゴルファーのスイング動作に関する第１スイング指標を算出する算出部と、前記第１スイング指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記ゴルファーに適した前記ゴルフクラブのシャフトの重量である最適シャフト重量を決定する決定部とを備える。

【0029】

第２２観点に係るフィッティング装置は、第２１観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記最適振り易さ指標の大きさ及び前記ユーザにより指定された前記ヘッドの種類に応じて、前記最適シャフト重量を決定する。

【0030】

第２３観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、取得部と、算出部と、決定部と、選択部とを備える。前記取得部は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する。前記算出部は、前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する。前記決定部は、前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正する。前記選択部は、前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する。

【0031】

第２４観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブを選定するためのフィッティング装置であって、取得部と、算出部と、決定部と、選択部とを備える。前記取得部は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する。前記算出部は、前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出する。前記決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定する。前記選択部は、前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択する。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。

【0032】

第２５観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0033】

第２６観点に係るフィッティングプログラムでは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標及び第２スイング指標を算出するステップ。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの特定の部位の特性を表す最適特性指標を決定するとともに、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したシャフトの剛性を示す最適剛性指標を決定するステップ。
（４）前記最適特性指標及び前記最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0034】

第２７観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。
（３）前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0035】

第２８観点に係るフィッティングプログラムでは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。
（３）前記第１スイング指標を表す空間を分割することにより定義される複数の領域のいずれに、前記計測値に基づいて算出された前記第１スイング指標が属するのか応じて、前記ゴルファーに適したゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定した後、前記計測値に基づいて算出された前記第１スイング指標が前記複数の領域の境界付近に存在する場合には、前記最適クラブ指標を修正するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0036】

第２９観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブを選定するためのフィッティング方法であって、以下の（１）〜（４）のステップを含む。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【0037】

第３０観点に係るフィッティングプログラムは、ゴルファーに適したゴルフクラブを選定するためのフィッティングプログラムであって、以下の（１）〜（４）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測値に基づいて、前記スイング動作に関する第１スイング指標を算出するステップ。前記第１スイング指標には、スイング平面上での前記スイング動作を表す指標が含まれるとともに、前記スイング平面に現れないローテーションの動き、及びコックの動きの少なくとも一方を表す指標が含まれる。
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブ全体又は特定の部位の特性を表す最適クラブ指標を決定するステップ。
（４）前記最適クラブ指標に合致するゴルフクラブ及びシャフトの少なくとも一方を選択するステップ。

【発明の効果】

【0038】

本発明の第１観点によれば、テストクラブによる計測値に基づいて、第１スイング指標及び第２スイング指標が算出される。そして、第１スイング指標に応じて、ゴルファーに適したゴルフクラブの振り易さ指標（例えば、ゴルフクラブのスイング慣性モーメント、グリップエンド周りの慣性モーメント、重量等）である最適振り易さ指標が決定されるとともに、第２スイング指標に応じて、ゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適振り易さ指標及び最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方が選択される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適した振り易さ指標を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。

【0039】

本発明の第１４観点によれば、テストクラブによる計測値に基づいて、第１スイング指標及び第２スイング指標が算出される。そして、第１スイング指標に応じて、ゴルファーに適したゴルフクラブの特定の部位の特性（例えば、シャフトの重量等）を表す最適特性指標が決定されるとともに、第２スイング指標に応じて、ゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適特性指標及び最適剛性指標に合致するゴルフクラブ及びシャフトの少なくとも一方が選択される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適したゴルフクラブの特定の部位の特性を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。

【0040】

本発明の第２１観点によれば、ユーザからゴルフクラブに含まれるべきヘッドの種類を決定するための情報の入力を受け付けるとともに、ゴルファーのスイング動作に関する第１スイング指標が算出される。そして、ユーザにより指定されたヘッドの種類に加え、第１スイング指標に応じて、ゴルファーに適したシャフトの重量（最適シャフト重量）が決定される。すなわち、ゴルファーに適しているかどうか定かではないマイクラブ重量のような指標ではなく、ゴルファーに適したシャフトの重量を基準にゴルフクラブ及びシャフトの少なくとも一方の絞り込みが行われる。従って、ゴルファーに適したゴルフクラブを精度よく選定することができる。

【図面の簡単な説明】

【0041】

【図1】本発明の一実施形態に係るフィッティング装置を備えるフィッティングシステムを示す図。

【図2】フィッティングシステムの機能ブロック図。

【図3】ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。

【図4】フィッティング処理の流れを示すフローチャート。

【図5】（Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。

【図6】第１変換処理の流れを示すフローチャート。

【図7】スイング平面を説明する図。

【図8】肩挙動導出工程の流れを示すフローチャート。

【図9】二重振り子モデルを概念的に説明する図。

【図10】第１指標算出工程の流れを示すフローチャート。

【図11】二重振り子モデルを概念的に説明する別の図。

【図12】腕エネルギーを説明する図。

【図13】最適振り易さ決定工程の流れを示すフローチャート。

【図14】プロモデル領域を示す図。

【図15】アベレージモデル領域を示す図。

【図16】インターナショナル・フレックス・コード（ＩＦＣ）を説明する図。

【図17】シャフトの曲げ剛性の測定方法を説明する図。

【図18】スイング中のシャフトの曲げを説明する図。

【図19】第２実施形態に係るフィッティングシステムの機能ブロック図。

【図20】第２実施形態に係るフィッティング処理の流れを示すフローチャート。

【図21】最適特性決定工程の流れを示すフローチャート。

【図22】様々なゴルフクラブのＩ_S，Ｉ_Gの関係をプロットしたグラフ。

【図23】様々なゴルフクラブのＩ_S，ｍ₂の関係をプロットしたグラフ。

【図24】変形例に係る最適振り易さ決定工程の流れを示すフローチャート。

【図25】別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。

【図26】変形例に係る第１スイング指標となる角度を説明する図。

【図27】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。

【図28】図２７の変形例に係るプロモデル領域を示す図。

【図29】図２７の変形例に係るアベレージモデル領域を示す図。

【図30】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。

【図31】図３０の変形例に係る平均肩トルク−腕エネルギー平面を示す図。

【図32】さらに別の変形例に係る最適振り易さ決定工程の流れを示すフローチャート。

【図33】変形例に係る肩挙動導出工程を概念的に説明する図。

【図34】変形例に係る肩挙動導出工程を概念的に説明する別の図。

【図35A】第１実施形態に係る二重振り子モデルのシミュレーション結果を示す図。

【図35B】変形例に係る二重振り子モデルのシミュレーション結果を示す図。

【図36】最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。

【図37A】特定のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。

【図37B】別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。

【図37C】さらに別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。

【図38】境界処理の流れを示す図。

【図39A】捩じりのローテーションの動きを示す概念図。

【図39B】プッシュのローテーションの動きを示す概念図。

【図40】例外処理の流れを示す図。

【図41A】実施例１及び比較例による飛距離を示すグラフ。

【図41B】実施例１及び比較例による左右ぶれを示すグラフ。

【図42A】実施例２の試打結果を示すグラフ。

【図42B】参考例の試打結果を示すグラフ。

【発明を実施するための形態】

【0042】

以下、図面を参照しつつ、本発明の幾つかの実施形態に係るゴルフクラブのフィッティング装置、方法及びプログラムについて説明する。

【0043】

＜１．第１実施形態＞
＜１−１．フィッティングシステムの概略構成＞
図１及び図２に、本実施形態に係るフィッティング装置２を備えるフィッティングシステム１００の全体構成を示す。フィッティング装置２は、ゴルファー７によるゴルフクラブ４のスイング動作を計測した計測データに基づいて、当該ゴルファー７に適したゴルフクラブ４を選定するための装置である。本実施形態では、スイング動作の計測は、ゴルフクラブ４のグリップ４２に取り付けられたセンサユニット１により行われ、フィッティング装置２は、このセンサユニット１とともに、フィッティングシステム１００を構成する。

【0044】

以下、センサユニット１及びフィッティング装置２の構成について説明した後、フィッティング処理の流れについて説明する。

【0045】

＜１−１−１．センサユニットの構成＞
センサユニット１は、図１及び図３に示すとおり、ゴルフクラブ４のグリップ４２におけるヘッド４１と反対側の端部に取り付けられており、グリップ４２の挙動を計測する。なお、ゴルフクラブ４は、一般的なゴルフクラブであり、シャフト４０と、シャフト４０の一端に設けられたヘッド４１と、シャフト４０の他端に設けられたグリップ４２とから構成される。本実施形態に係るシャフト４０は、カーボン製のシャフトである。センサユニット１は、スイング動作の妨げとならないよう、小型且つ軽量に構成されている。図２に示すように、本実施形態に係るセンサユニット１には、加速度センサ１１、角速度センサ１２及び地磁気センサ１３が搭載されている。また、センサユニット１には、これらのセンサ１１〜１３による計測データを外部のフィッティング装置２に送信するための通信装置１０も搭載されている。なお、本実施形態では、通信装置１０は、スイング動作の妨げにならないように無線式であるが、ケーブルを介して有線式にフィッティング装置２に接続するようにしてもよい。

【0046】

加速度センサ１１、角速度センサ１２及び地磁気センサ１３はそれぞれ、グリップ４２を基準としたｘｙｚ局所座標系におけるグリップ加速度、グリップ角速度及びグリップ地磁気を計測する。より具体的には、加速度センサ１１は、ｘ軸、ｙ軸及びｚ軸方向のグリップ加速度ａ_x，ａ_y，ａ_zを計測する。角速度センサ１２は、ｘ軸、ｙ軸及びｚ軸周りの
グリップ角速度ω_x，ω_y，ω_zを計測する。地磁気センサ１３は、ｘ軸、ｙ軸及びｚ軸方
向のグリップ地磁気ｍ_x，ｍ_y，ｍ_zを計測する。これらの計測データは、所定のサンプリ
ング周期Δｔの時系列データとして取得される。なお、ｘｙｚ局所座標系は、図３に示すとおりに定義される３軸直交座標系である。すなわち、ｚ軸は、シャフト４０の延びる方向に一致し、ヘッド４１からグリップ４２に向かう方向が、ｚ軸正方向である。ｘ軸は、ヘッド４１のトゥ−ヒール方向にできる限り沿うように配向され、ｙ軸は、ヘッド４１のフェース面の法線方向にできる限り沿うように配向される。

【0047】

本実施形態では、加速度センサ１１、角速度センサ１２及び地磁気センサ１３による計測データは、通信装置１０を介してリアルタイムにフィッティング装置２に送信される。しかしながら、例えば、センサユニット１内の記憶装置に計測データを格納しておき、スイング動作の終了後に当該記憶装置から計測データを取り出して、フィッティング装置２に受け渡すようにしてもよい。

【0048】

＜１−１−２．フィッティング装置の構成＞
図２を参照しつつ、フィッティング装置２の構成について説明する。フィッティング装置２は、ＣＤ−ＲＯＭ、ＵＳＢメモリ等のコンピュータで読み取り可能な記録媒体２０に格納された本実施形態に係るフィッティングプログラム３を、当該記録媒体２０から汎用のパーソナルコンピュータにインストールすることにより製造される。フィッティングプログラム３は、センサユニット１から送られてくる計測データに基づいてスイング動作を解析し、ゴルファー７に適したゴルフクラブ４を選択するのを支援する情報を出力するためのソフトウェアである。フィッティングプログラム３は、フィッティング装置２に後述する動作を実行させる。

【0049】

フィッティング装置２は、表示部２１、入力部２２、記憶部２３、制御部２４及び通信部２５を備える。そして、これらの部２１〜２５は、バス線２６を介して接続されており、相互に通信可能である。本実施形態では、表示部２１は、液晶ディスプレイ等で構成され、後述する情報をユーザに対し表示する。なお、ここでいうユーザとは、ゴルファー７自身やそのインストラクター等の、フィッティングの結果を必要とする者の総称である。また、入力部２２は、マウス、キーボード、タッチパネル等で構成することができ、フィッティング装置２に対するユーザからの操作を受け付ける。通信部２５は、フィッティング装置２と外部装置との通信を可能にする通信インターフェースであり、センサユニット１からデータを受信する。

【0050】

記憶部２３は、ハードディスク等の不揮発性の記憶装置により構成される。記憶部２３内には、フィッティングプログラム３が格納されている他、センサユニット１から送られてくる計測データが保存される。また、記憶部２３内には、対応関係データ２８、ヘッドデータベース（ＤＢ）２７及びシャフトデータベース（ＤＢ）２９が格納されている。対応関係データ２８とは、詳細は後述するが、ゴルフクラブ４の様々なモデル（シリーズ）毎に規定されており、最適振り易さ指標を決定するための条件を示すデータである。同様に詳細は後述するが、ヘッドＤＢ２７は、多数のヘッド４１のスペックを示す情報が、ヘッド４１の種類を特定する情報に関連付けて格納されたデータベースである。シャフトＤＢ２９は、多数のシャフト４０のスペックを示す情報が、シャフト４０の種類を特定する情報に関連付けて格納されたデータベースである。シャフトＤＢ２９には、例えば、１００を超える種類のシャフト４０のデータが格納されている。

【0051】

制御部２４は、ＣＰＵ、ＲＯＭおよびＲＡＭ等から構成することができる。制御部２４は、記憶部２３内のフィッティングプログラム３を読み出して実行することにより、仮想的に取得部２４Ａ、グリップ挙動導出部２４Ｂ、肩挙動導出部２４Ｃ、算出部２４Ｄ、決定部２４Ｅ、選択部２４Ｆ及び表示制御部２４Ｇとして動作する。各部２４Ａ〜２４Ｇの動作の詳細については、後述する。

【0052】

＜１−２．フィッティング処理＞
続いて、フィッティングシステム１００により実行されるフィッティング処理について説明する。本実施形態に係るフィッティング処理は、図４に示すとおり、以下の９つの工程（Ｓ１〜Ｓ９）から構成されている。
（Ｓ１）ｘｙｚ局所座標系でのグリップ加速度ａ_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）である最適振り易さ指標（本実施形態では、最適スイングＭＩ）を決定する最適振り易さ決定工程
（Ｓ７）計測データに基づいて、第２スイング指標（本実施形態では、後述する第１〜第４特徴量Ｆ₁〜Ｆ₄）を算出する第２指標算出工程
（Ｓ８）第２スイング指標に基づいて、ゴルファー７に適したシャフト４０の剛性を示す最適剛性指標（本実施形態では、後述するＥＩ分布）を決定する最適剛性決定工程
（Ｓ９）最適振り易さ指標及び最適剛性指標に合致するゴルフクラブ４、特にシャフト４０を選択する最適クラブ選択工程
以下、これらの工程を順に説明する。

【0053】

なお、ＸＹＺ全体座標系は、図１に示すとおりに定義される３軸直交座標系である。すなわち、Ｚ軸は、鉛直下方から上方に向かう方向であり、Ｘ軸は、ゴルファー７の背から腹に向かう方向であり、Ｙ軸は、地平面に平行でボールの打球地点から目標地点に向かう方向である。

【0054】

＜１−２−１．計測工程＞
計測工程（Ｓ１）では、ゴルファー７により、上述のセンサユニット１付きゴルフクラブ４がスイングされる。計測工程でスイングされるゴルフクラブ４は、２本のテストクラブのうちの１本である。これらのテストクラブは、異なる種類のゴルフクラブであり、本実施形態では、１本はプロ仕様のゴルフクラブ（以下、プロモデルクラブ）であり、もう１本はアベレージユーザーに適したゴルフクラブ（以下、アベレージモデルクラブ）である。また、本実施形態では、プロモデルクラブは、アベレージモデルクラブよりも重量が大きい。計測工程でいずれのテストクラブがスイングされるかは、ゴルファー７の好みや経験等に基づいて、決定される。

【0055】

なお、フィッティングの精度を向上させるために、ユーザのマイクラブの長さが、シャフトＤＢ２９に記憶されているシャフトに基づくゴルフクラブの長さと異なる場合は、テストクラブの総重量をマイクラブの長さ相当のクラブ総重量に変更する。これにより、当該ユーザにマッチしたシャフトを選定することができる。例えば、シャフトＤＢ２９に記憶されているクラブ長さが全て４５インチ（＝１１４３ｍｍ）であり、ユーザのマイクラブの長さＡ（ｍｍ）が４５インチと異なる場合、テストクラブの総重量を次の式により算出したもの（４５インチ相当の総重量）に変更して、フィッティングを行う。
（テストクラブ総重量）＝（Ａ−１１４３）×０．３７７＋（マイクラブのクラブ総重量）

【0056】

続いて、以上のようなゴルフクラブ４のスイング動作中のグリップ加速度ａ_x，ａ_y，ａ_z、グリップ角速度ω_x，ω_y，ω_z及びグリップ地磁気ｍ_x，ｍ_y，ｍ_zの計測データが、センサユニット１により計測される。この計測データは、センサユニット１の通信装置１０を介してフィッティング装置２に送信される。一方、フィッティング装置２側では、取得部２４Ａが通信部２５を介してこれを受信し、記憶部２３内に格納する。本実施形態では、少なくともアドレスからインパクトまでの時系列の計測データが計測される。

【0057】

なお、ゴルフクラブのスイング動作は、一般に、アドレス、トップ、インパクト、フィニッシュの順に進む。アドレスとは、図５（Ａ）に示すとおり、ゴルフクラブ４のヘッド４１をボール近くに配置した初期の状態を意味し、トップとは、図５（Ｂ）に示すとおり、アドレスからゴルフクラブ４をテイクバックし、最もヘッド４１が振り上げられた状態を意味する。インパクトとは、図５（Ｃ）に示すとおり、トップからゴルフクラブ４が振り下ろされ（ダウンスイング）、ヘッド４１がボールと衝突した瞬間の状態を意味し、フィニッシュとは、図５（Ｄ）に示すとおり、インパクト後、ゴルフクラブ４を前方へ振り抜いた状態を意味する。

【0058】

計測工程では、以上のゴルフクラブ４が複数回、好ましくは５回以上試打されることが好ましい。この場合、計測データの平均値を算出し、以降の演算に使用することができる。また、ミスショットや計測ミス等による異常値を取り除くため、計測データの標準偏差σを算出するようにし、全ての計測データが平均値±ｋ・σ（ｋは、定数）以内に収まっていない場合には、計測の追加又はやり直しを求めるメッセージを表示部２１上に表示させるようにしてもよい。なお、計測データ自体の平均値ではなく、計測データに基づいて算出される加工値（例えば、後述する腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_h）の平均値を算出するようにしてもよい。加工値の平均値を算出する場合も、同じく標準偏差σに基づくデータの信頼性のチェックを行うことができる。

【0059】

＜１−２−２．第１変換工程＞
以下、図６を参照しつつ、ｘｙｚ局所座標系の計測データをＸＹＺ全体座標系の値へと変換する第１変換工程（Ｓ２）について説明する。具体的には、まず、取得部２４Ａが、記憶部２３内に格納されているｘｙｚ局所座標系でのグリップ加速度ａ_x，ａ_y，ａ_z、グリップ角速度ω_x，ω_y，ω_z及びグリップ地磁気ｍ_x，ｍ_y，ｍ_zの時系列の計測データを読み出す（ステップＳ１１）。

【0060】

次に、ステップＳ１１で読み出されたｘｙｚ局所座標系での時系列の計測データに基づいて、グリップ挙動導出部２４Ｂが、インパクト、トップ及びアドレスの時刻ｔ_i，ｔ_t，ｔ_aを導出する（ステップＳ１２）。これらの時刻ｔ_i，ｔ_t，ｔ_aの導出方法としては、様々なものが公知であるため、ここでは詳細な説明を省略する。

【0061】

続くステップＳ１３では、グリップ挙動導出部２４Ｂが、アドレスからインパクトまでの時刻ｔにおける姿勢行列Ｎ（ｔ）を算出する。今、姿勢行列を以下の式で表すとする。姿勢行列Ｎ（ｔ）は、時刻ｔにおけるＸＹＺ全体座標系をｘｙｚ局所座標系に変換するための行列である。

【数1】

【0062】

姿勢行列Ｎ（ｔ）の９つの成分の意味は、以下のとおりである。
成分ａ：全体座標系のＸ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｂ：全体座標系のＹ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｃ：全体座標系のＺ軸と、局所座標系のｘ軸とのなす角度の余弦
成分ｄ：全体座標系のＸ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｅ：全体座標系のＹ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｆ：全体座標系のＺ軸と、局所座標系のｙ軸とのなす角度の余弦
成分ｇ：全体座標系のＸ軸と、局所座標系のｚ軸とのなす角度の余弦
成分ｈ：全体座標系のＹ軸と、局所座標系のｚ軸とのなす角度の余弦
成分ｉ：全体座標系のＺ軸と、局所座標系のｚ軸とのなす角度の余弦
ここで、ベクトル（ａ，ｂ，ｃ）は、ｘ軸方向の単位ベクトルを表し、ベクトル（ｄ，ｅ，ｆ）は、ｙ軸方向の単位ベクトルを表し、ベクトル（ｇ，ｈ，ｉ）は、ｚ軸方向の単位ベクトルを表している。

【0063】

また、姿勢行列Ｎ（ｔ）は、Ｚ−Ｙ−Ｚ系のオイラー角の考え方に従うと、以下の式で表すことができる。ただし、φ，θ，ψは、それぞれＺ軸、Ｙ軸、Ｚ軸周りの回転角度とする。

【数2】

【0064】

アドレスからインパクトまでの姿勢行列Ｎ（ｔ）を算出するに当たり、まず、アドレスの時刻ｔ_aにおける姿勢行列Ｎ（ｔ_a）が算出される。具体的には、以下の式に従って、アドレス時のφ，θが算出される。なお、以下の式は、アドレス時にはゴルフクラブ４は静止しており、加速度センサ１１によって鉛直方向の重力のみが検出されることを利用している。以下の式中のグリップ加速度ａ_x，ａ_y，ａ_zは、アドレス時の値である。

【数3】

【数4】

【0065】

続いて、以下の式に従って、アドレス時のψが算出される。

【数5】

ただし、上式中のｍ_xi，ｍ_yiの値は、以下の式に従って算出される。また、以下の式中のグリップ地磁気ｍ_x，ｍ_y，ｍ_zは、アドレス時の値である。

【数6】

【0066】

以上より、アドレス時のφ，θ，ψが、ｘｙｚ局所座標系でのグリップ加速度ａ_x，ａ_y，ａ_z及びグリップ地磁気ｍ_x，ｍ_y，ｍ_zに基づいて算出される。そして、これらのφ，θ，ψの値を数２の式に代入することにより、アドレス時の姿勢行列Ｎ（ｔ_a）が算出される。

【0067】

続いて、アドレス時の姿勢行列Ｎ（ｔ_a）をサンプリング周期Δｔ間隔で時々刻々更新してゆくことにより、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）が算出される。具体的に説明すると、まず、姿勢行列Ｎ（ｔ）は、クォータニオンの４変数ｑ₁，ｑ₂，ｑ₃，ｑ₄（ｑ₄がスカラー部）を用いて、以下の式で表される。

【数7】

【0068】

従って、数１及び数７より、クォータニオンの４変数ｑ₁，ｑ₂，ｑ₃，ｑ₄は、以下の式に従って、算出することができる。

【数8】

【0069】

今、アドレス時の姿勢行列Ｎ（ｔ_a）を規定するａ〜ｉの値は既知である。よって、以上の式に従って、まず、アドレス時のクォータニオンの４変数ｑ₁，ｑ₂，ｑ₃，ｑ₄が算出される。

【0070】

そして、時刻ｔから微小時刻経過後のクォータニオンｑ'は、時刻ｔにおけるクォータニオンｑを用いて以下の式で表される。

【数9】

【0071】

また、クォータニオンの４変数ｑ₁，ｑ₂，ｑ₃，ｑ₄の時間変化を表す１階微分方程式は、以下の式で表される。

【数10】

【0072】

数９及び数１０の式を用いれば、時刻ｔのクォータニオンを順次、次の時刻ｔ＋Δｔのクォータニオンへと更新することができる。ここでは、アドレスからインパクトまでのクォータニオンが算出される。そして、アドレスからインパクトまでのクォータニオンの４変数ｑ₁，ｑ₂，ｑ₃，ｑ₄を数７の式に順次代入してゆくことにより、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）が算出される。

【0073】

次に、ステップＳ１４では、グリップ挙動導出部２４Ｂが、アドレスからインパクトまでの姿勢行列Ｎ（ｔ）に基づいて、アドレスからインパクトまでのｘｙｚ局所座標系でのグリップ加速度ａ_x，ａ_y，ａ_z及びグリップ角速度ω_x，ω_y，ω_zの時系列データを、ＸＹＺ全体座標系での時系列データに変換する。変換後のグリップ加速度ａ_X，ａ_Y，ａ_Z及びグリップ角速度ω_X，ω_Y，ω_Zは、以下の式に従って算出される。

【数11】

【0074】

続くステップＳ１５では、グリップ挙動導出部２４Ｂは、グリップ加速度ａ_X，ａ_Y，ａ_Zの時系列データを積分することにより、アドレスからインパクトまでのＸＹＺ全体座標系でのグリップ速度ｖ_X，ｖ_Y，ｖ_Zを導出する。このとき、アドレスからインパクトまでのグリップ速度ｖ_X，ｖ_Y，ｖ_Zを、トップにおいて０ｍ／ｓとなるように、オフセットを行うことが好ましい。

【0075】

＜１−２−３．第２変換工程＞
以下、第１変換工程で算出されたＸＹＺ全体座標系でのグリップ４２の挙動を、スイング平面Ｐ内でのグリップ４２の挙動へと変換する第２変換工程（Ｓ３）について説明する。本実施形態では、スイング平面Ｐは、ＸＹＺ全体座標系の原点を含み、Ｙ軸及びインパクト時のシャフト４０と平行な面として定義される（図７参照）。第２変換工程では、グリップ挙動導出部２４Ｂは、ＸＹＺ全体座標系でのグリップ速度ｖ_X，ｖ_Y，ｖ_Z及びグリップ角速度ω_X，ω_Y，ω_Zをスイング平面Ｐ内へ射影したグリップ速度ｖ_pX，ｖ_pY，ｖ_pZ及びグリップ角速度ω_pX，ω_pY，ω_pZを算出する。

【0076】

具体的には、シャフト４０の延びる方向を表す、姿勢行列Ｎ（ｔ）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）に基づいて、Ｘ軸正方向から見た（ゴルファー７を正面から見た）シャフト４０の傾きの時系列データを算出する。そして、この時系列データに基づいて、Ｘ軸正方向から見てシャフト４０がＺ軸と平行になる時刻を特定し、これをインパクトの時刻ｔ_iとする。なお、ここでのインパクトの時刻ｔ_iは、既出のインパクトの時刻ｔ_iと一致するとは限らない。続いて、このインパクトの時刻ｔ_iにおける姿勢行列Ｎ（ｔ_i）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）に基づいて、Ｙ軸負方向から見たシャフト４０の傾きを算出する。すなわち、インパクト時にＹ軸負方向から見たシャフト４０とＸ軸との為す角度α'を算出し、これをスイング平面角度とする。

【0077】

スイング平面角度α'が求まると、これを用いてＸＹＺ全体座標系における任意の点をスイング平面Ｐに射影するための射影変換行列Ａを、以下のとおり算出することができる。ただし、α＝９０°−α'である。

【数12】

【0078】

ここでは、以上の射影変換行列Ａに基づいて、以下の式に従って、アドレスからインパクトまでの射影変換後のグリップ速度ｖ_pX，ｖ_pY，ｖ_pZ及びグリップ角速度ω_pX，ω_pY，ω_pZの時系列データが算出される。

【数13】

【0079】

なお、以上の演算により得られるグリップ速度（ｖ_pY，ｖ_pZ）は、スイング平面Ｐ内でのグリップ速度（ベクトル）を表しており、グリップ角速度ω_pXは、スイング平面Ｐに対して垂直な軸周りの角速度を表している。ここでは、以下の式に従って、アドレスからインパクトまでのスイング平面Ｐ内でのグリップ速度（スカラー）が算出される。

【数14】

【0080】

また、ここでは、後の計算に必要となる、スイング平面Ｐ内におけるトップでのシャフト４０の傾きも算出される。具体的には、まず、トップでの姿勢行列Ｎ（ｔ_t）に含まれるｚ軸ベクトル（ｇ，ｈ，ｉ）を、射影変換行列Ａを用いて、以下の式に従ってスイング平面Ｐ内に射影する。ただし、射影後のベクトルを（ｇ'，ｈ'，ｉ'）とする。

【数15】

【0081】

以上の式により特定されるベクトル（ｈ'，ｉ'）は、スイング平面Ｐ内におけるトップでのシャフト４０の傾きを表すベクトルである。従って、以上の計算結果を以下の式に代入することにより、スイング平面Ｐ内におけるトップでのシャフト４０の傾きβが算出される。

【数16】

【0082】

＜１−２−４．肩挙動導出工程＞
以下、図８を算出しつつ、スイング平面Ｐ内でのグリップの挙動（グリップ速度Ｖ_GE及びグリップ角速度ω_pX）に基づいて、スイング平面Ｐ内の疑似的な肩の挙動を導出する肩挙動導出工程（Ｓ４）について説明する。本実施形態では、ゴルフクラブ４の挙動は、ゴルファー７の肩及びグリップ４２（或いは、これを握るゴルファーの手首）を節点とし、ゴルファー７の腕及びゴルフクラブ４をリンクとする二重振り子モデルに基づいて解析される。ただし、肩の挙動は直接的に実測されるのではなく、実測されたグリップの挙動に基づいて、疑似的な肩の挙動として導出される。以下では、特に断らない限り、単に「肩」という場合も、このような疑似的な肩を意味するものとする。疑似的な肩とグリップ４２（手首）との間を直線的に延びるものとして定義される疑似的な「腕」についても同様である。

【0083】

グリップの挙動から肩の挙動を特定するに当たり、本実施形態に係る二重振り子モデルは、以下の（１）〜（５）を前提とする。図９は、以下の前提条件を概念的に説明する図である。
（１）スイング平面Ｐ上において、グリップ４２（手首）は肩を中心として円運動する。（２）スイング平面Ｐ上において、肩とグリップ４２との距離（半径）Ｒは、一定である。
（３）肩は、スイング動作中は動かない（ただし、回転する）。
（４）スイング平面Ｐ上において、トップでの腕とゴルフクラブ４との為す角度は９０°である。
（５）インパクト時の腕は、Ｘ軸正方向から見てＺ軸下方を向く。

【0084】

以上の前提の下、肩挙動導出部２４Ｃは、スイング平面Ｐ内におけるトップからインパクトまでのグリップ４２の移動距離Ｄを算出する（ステップＳ２１）。移動距離Ｄは、トップからインパクトまでのグリップ速度Ｖ_GEを積分することにより導出される。

【0085】

続いて、肩挙動導出部２４Ｃは、スイング平面Ｐ内におけるトップからインパクトまでの腕の回転角度γを算出する（ステップＳ２２）。回転角度γは、第２変換工程で算出されたトップでのシャフト４０の傾きβに基づいて算出される。次に、肩挙動導出部２４Ｃは、半径Ｒ＝Ｄ／γを算出する（ステップＳ２３）。

【0086】

そして、肩挙動導出部２４Ｃは、以下の式に従って、肩の挙動として、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω₁を算出する（ステップＳ２４）。すなわち、腕の角速度ω₁は、計測によるグリップ速度Ｖ_GEが反映さ
れた値となる。
ω₁＝Ｖ_GE／Ｒ

【0087】

＜１−２−５．第１指標算出工程＞
以下、図１０を参照しつつ、グリップ４２の挙動及び肩の挙動に基づいて、第１スイング指標を算出する第１指標算出工程（Ｓ５）について説明する。第１スイング指標とは、最適振り易さ指標を決定するための、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態では、第１スイング指標として、後述される腕出力パワーＰ_{1_AVE}
、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hが算出される。

【0088】

具体的には、まず、ステップＳ３１では、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω₁を積分し、トップからインパクトまでの腕の回転角度θ₁を算出する。このとき、台形積分を用いることが好ましい。なお、回転角度θ₁は、図１１のように定義され、図１１の紙面は、スイング平面Ｐに等しい。以下では、図１１に示されるスイング平面Ｐ内での新たなＸＹ座標系に基づいて、解析が進められる。スイング平面Ｐ内での新たなＸＹ座標系のＸ軸は、上述したＸＹＺ全体座標系のＹ軸に等しく、新たなＸＹ座標系のＹ軸は、ＸＹＺ全体座標系のＺ軸をスイング平面Ｐ内に投影した軸である。

【0089】

また、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω₁を微分し、トップからインパクトまでの角加速度ω₁'を算出する。次に、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の重心の位置（Ｘ₁，Ｙ₁）、速度（Ｖ_X1，Ｖ_Y1）及び加速度（Ａ_X1，Ａ_Y1）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。

【数17】

【0090】

ただし、ｒは、肩から腕の重心までの距離である。本実施形態では、腕の重心は、腕の中心にあるものと仮定される。従って、Ｒ＝２ｒである。

【0091】

次に、ステップＳ３２では、グリップ挙動導出部２４Ｂは、ステップＳ３１と同様の演算をグリップ４２周りについても行う。すなわち、トップからインパクトまでのグリップ角速度ω_pX＝グリップ４２周りのゴルフクラブ４の角速度ω₂を積分し、トップからインパクトまでのグリップ４２周りのゴルフクラブ４（シャフト４０）の回転角度θ₂を算出する。このときも、台形積分を用いることが好ましく、回転角度θ₂は、図１１のように定義される。

【0092】

続いて、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の角速度ω₂を微分し、トップからインパクトまでの角加速度ω₂'を算出する。次に、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の重心の位置（Ｘ₂，Ｙ₂）、速度（Ｖ_X2，Ｖ_Y2）及び加速度（Ａ_X2，Ａ_Y2）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。

【数18】

【0093】

ただし、Ｌは、グリップ４２からゴルフクラブ４の重心までの距離である。Ｌの値は、ゴルフクラブ４のスペックであり、予め定められているものとする。

【0094】

次に、ステップＳ３３では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの肩に発生する拘束力Ｒ₁＝（Ｒ_X1，Ｒ_Y1）を算出するとともに、トップからインパクトまでのグリップ４２に発生する拘束力Ｒ₂＝（Ｒ_X2，Ｒ_Y2）を算出する。以下の式は、並進方向の力の釣り合いに基づくものである。ただし、ｍ₁は、腕の質量であり、本実施形態では、腕の質量ｍ₁は、適宜予め定められているものとする。例えば、解析を開始する前に、ゴルファー７の体重を入力しておき、入力された体重に所定の係数を掛ける等して、自動的に腕の質量が算出される。ｍ₂は、ゴルフクラブ４の質量であり、ｇは、重力加速度である。また、ｍ₂は、ゴルフクラブ４のスペックであり、予め定められているものとする。

【数19】

【0095】

続くステップＳ３４では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の重心周りのトルクＴ_g1及びゴルフクラブ４の重心周りのトルクＴ_g2を算出する。

【数20】

ただし、Ｉ₁は、腕の重心周りの慣性モーメントであり、Ｉ₂は、ゴルフクラブ４の重心周りの慣性モーメントである。本実施形態では、腕の重心周りの慣性モーメントＩ₁は、腕の重心が腕の中心にあるとの仮定の下、Ｉ₁＝ｍ₁ｒ²／３として算出される。また、I₂は、ゴルフクラブ４のスペックであり、予め定められているものとする。

【0096】

続くステップＳ３５では、算出部２４Ｄは、上述した計算結果に基づいて、トップからインパクトまでの腕の仕事率（パワー）Ｅ₁'を算出する。具体的には、Ｅ₁'は、肩の速度ベクトルをｖ_sとし、グリップ４２の速度ベクトルをｖ_gとして、以下の式に従って表される。また、ｖ_s，ｖ_gはそれぞれ、肩の位置ベクトルｄ_s、グリップ４２の位置ベクトルｄ_g＝ｄ_s＋（２Ｘ₁，２Ｙ₁）を一階微分することにより算出可能である。

【数21】

【0097】

また、本実施形態では、肩は動かないため、ｖ_s＝（０，０）となり、腕の仕事率Ｅ₁'は、以下の式に従って算出される。算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の仕事率Ｅ₁'を算出する。

【数22】

【0098】

ところで、ゴルフスイングにおいて、ゴルフクラブ４の先端（ヘッド４１）を最も加速させるためには、まず腕を十分に加速させて、その後、腕の動きを止めてゴルフクラブ４に勢いを与えることが求められると考えられる。ここでいう腕の加速具合とは、腕が出力するパワー（腕出力パワー）Ｐ₁という物理指標に置き換えることができ、ゴルフクラブ４に与える勢いとは、ゴルフクラブ４に入力されるパワー（クラブ入力パワー）Ｐ₂という物理指標に置き換えることができる。そして、腕出力パワーＰ₁とは、腕の仕事率Ｅ₁'を表す数２２の右辺の第２項及び第３項部分に相当する。また、クラブ入力パワーＰ₂とは、数２２の式中の右辺の第１項部分に相当する。すなわち、腕出力パワーＰ₁及びクラブ入力パワーＰ₂は、以下のとおり表すことができる。ステップＳ３５では、算出部２４Ｄは、腕の仕事率Ｅ₁'に加え、トップからインパクトまでの腕出力パワーＰ₁及びクラブ入力パワーＰ₂を算出する。

【数23】

【0099】

なお、スイング動作中にゴルフクラブ４で発揮される仕事率Ｅ₂'は、下式のように表すことができる。すなわち、クラブ入力パワーＰ₂＝Ｒ₂ｖ_g^Tが橋渡しになって、腕からゴルフクラブ４へとエネルギーが伝達される。

【数24】

【0100】

続くステップＳ３６では、算出部２４Ｄは、トップ以降で腕の仕事率Ｅ₁'が正から負へ転じる時刻ｔ_cを特定し、トップの時刻ｔ_tから時刻ｔ_cまでの腕の仕事量Ｅ₁を算出する。腕の仕事量Ｅ₁は、時刻ｔ_t〜ｔ_cの区間で腕の仕事率Ｅ₁'を積分することにより、算出される（図１２参照）。なお、仕事量Ｅ₁は、時刻ｔ_t〜ｔ_cの間に腕で発揮される仕事量（エネルギー）を表す指標と考えることができるから、この意味で、スイング動作中の腕エネルギーと呼ぶことができる。また、算出部２４Ｄは、時刻ｔ_t〜ｔ_cの間に腕で発揮される平均仕事率Ｅ_AVE＝Ｅ₁／ｔ_c―ｔ_tを算出する。平均仕事率Ｅ_AVEは、スイング動作中に単位時間当たりに平均的に発揮ないし消費される腕エネルギーである。

【0101】

また、算出部２４Ｄは、トップの時刻ｔ_tから腕出力パワーＰ₁が最大値をとる時刻ｔ_mまでの区間で腕出力パワーＰ₁を積分し、この積分値Ｄ₁を積分区間（ｔ_c―ｔ_m）で除算することで、スイング動作中の平均的な腕出力パワーＰ_{1_AVE}を算出する。なお、この積分値Ｄ₁は、スイング動作中にゴルファーの腕が行う仕事量であり、腕出力パワーを示す指標となり得る。同様に、算出部２４Ｄは、トップの時刻ｔ_tからクラブ入力パワーＰ₂が最大値をとる時刻ｔ_nまでの区間でクラブ入力パワーＰ₂を積分し、この積分値Ｄ₂を積分区間（ｔ_c―ｔ_n）で除算することで、スイング動作中の平均的なクラブ入力パワーＰ_{2_AVE}を算出する。なお、この積分値Ｄ₂は、スイング動作中のゴルフクラブ４に対してされる仕事量であり、クラブ入力パワーを示す指標となり得る。なお、ここで示した積分区間は例示であり、例えば、時刻ｔ_t〜ｔ_cまでの区間等、適宜設定し得る。

【0102】

続くステップＳ３７では、算出部２４Ｄは、スイング動作中のコック解放タイミングｔ_rを算出する。なお、本発明者らは、実験を通して、インパクト時のヘッド速度Ｖ_hが、スイング動作中のコック解放タイミングｔ_r、及び、腕エネルギーＥ₁又は平均仕事率Ｅ_AVEと相関があることを発見した。そこで、ここでは、インパクト時のヘッド速度Ｖ_hを算出すべく、コック解放タイミングｔ_rが算出される。本実施形態では、コック解放タイミングｔ_rは、時刻ｔ_t〜ｔ_cまでの区間で腕の仕事率Ｅ₁'が最大となる時刻が、コック解放タイミングｔ_rとして特定される（図１２参照）。

【0103】

続くステップＳ３８では、算出部２４Ｄは、コック解放タイミングｔ_r及び腕エネルギーＥ_AVEに基づいて、インパクト時のヘッド速度Ｖ_hを算出する。具体的には、インパクト時のヘッド速度Ｖ_hは、下の式に従って算出される。なお、ｋ₁，ｋ₂，ｋ₃は、予め行われた多数の実験結果から重回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。以上より、指標算出工程が終了する。
Ｖ_h＝ｋ₁・Ｅ_AVE＋ｋ₂・ｔ_r＋ｋ₃

【0104】

以上の処理により、最適振り易さ指標を決定するための第１スイング指標として、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hが算出される。

【0105】

＜１−２−６．最適振り易さ決定工程＞
以下、図１３を参照しつつ、最適振り易さ決定工程（Ｓ６）の流れについて説明する。本実施形態に係る最適振り易さ指標とは、最適スイングＭＩであり、最適スイングＭＩとは、ゴルファーに適したゴルフクラブ４のスイング慣性モーメントＩ_Sのことを言う。

【0106】

また、スイング慣性モーメントＩ_Sとは、スイング中の肩周りの慣性モーメントであり、例えば、以下の式に従って定義することができる。
Ｉ_S＝Ｉ₂＋ｍ₂（Ｒ＋Ｌ）²＋Ｉ₁＋ｍ₁（Ｒ／２）²
なお、各ゴルファー７については、ゴルフクラブ４が変わっても腕の重量は同じである。従って、本実施形態では、簡単のため、スイング慣性モーメントＩ_Sは、腕の回転分の慣性モーメントを省略し、以下の式に従って算出される。
Ｉ_S＝Ｉ₂＋ｍ₂（Ｒ＋Ｌ）²
さらに、本実施形態では、腕長さＲ＝６０ｃｍ（一定）として、スイング慣性モーメントＩ_Sが計算される。ただし、上式中のＲには、ステップＳ２３で算出される腕長さＲの値を代入することもできる。ところで、Ｉ_Sを決定するパラメータであるｍ₂，Ｉ₂，Ｌは、ゴルフクラブ４の諸元である。従って、本実施形態でのスイング慣性モーメントも、ゴルフクラブ４の諸元となる。

【0107】

まず、ステップＳ４０において、決定部２４Ｅは、計測工程で試打されたテストクラブの種類を判定する。プロモデルクラブが試打されていた場合には、ステップＳ４１に進み、アベレージモデルクラブが試打されていた場合には、ステップＳ５１に進む。いずれのテストクラブが試打されたかは、入力部２２を介してユーザにより入力される情報に基づいて判定されるものとする。

【0108】

なお、ステップＳ４０においてテストクラブの種類が判定されるのは、テストクラブの種類に応じて第１スイング指標の分布する領域が異なるからである。具体的に説明すると、本発明者らは、プロモデルのゴルフクラブを通常使用しているゴルファー（以下、プロモデルユーザ）２１名及びアベレージモデルのゴルフクラブを通常使用しているゴルファー（以下、アベレージモデルユーザ）２２名に、それぞれプロモデルクラブ及びアベレージモデルクラブを試打させて第１スイング指標を算出したところ、図１４及び図１５に示す結果を得た。なお、この実験で算出された第１スイング指標は、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hであり、具体的な値は、上記と同様の工程に従って算出された。そして、本発明者らは、図１４及び図１５に示す結果から、第１スイング指標を示す空間は、プロモデル領域とアベレージモデル領域とに概ね分割されることを発見した。なお、プロモデル領域とは、プロモデルユーザによるスイング動作時の第１スイング指標が分布する領域であり、アベレージモデル領域とは、アベレージモデルユーザによるスイング動作時の第１スイング指標が分布する領域である。図１４及び図１５を比較すれば分かるように、プロモデル領域は、アベレージモデル領域よりも、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hがいずれも大きい領域である。しかしながら、これらの数値は、テストクラブの種類等のスイング条件によって異なるものとなり得る。本実験では、プロモデルクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ Ｋｏｓｕｍａ Ｂｌｕｅ６ Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１４ｇ、バランスは、Ｄ３）のドライバーが用いられ、アベレージモデルクラブとして、ダンロップスポーツ株式会社製のＸＸＩＯ（登録商標）８（シャフトは、ＭＰ−８００ Ｒ−Ｆｌｅｘ、ゴルフクラブの重量は、２７２ｇ、バランスは、Ｄ３）のドライバーが用いられた。

【0109】

続いて、ステップＳ４１及びこれに続くステップＳ４２〜Ｓ４４について説明する。ステップＳ４１〜Ｓ４４は、腕出力パワーＰ_{1_AVE}及びクラブ入力パワーＰ_{2_AVE}の大きさに応じて、最適スイングＭＩの範囲（以下、最適スイングＭＩ帯）を決定するステップである。ここでは、Ｐ_{1_AVE}，Ｐ_{2_AVE}の値が大きい程、最適スイングＭＩ帯は段階的に大きな値に設定される。

【0110】

具体的には、ステップＳ４１では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１４に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ１の上側にあるか、すなわち、図１４の領域Ａ１に属するか否か（以下、プロ条件１）を判定する。そして、プロ条件１が満たされる場合には、５６２３ｋｇ・ｃｍ²以上が最適スイングＭＩ帯であると判定する。一方、ステップＳ４１でプロ条件１が満たされない場合には、ステップＳ４２に進む。ステップＳ４２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１４に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ２の上側かつ直線Ｌ１よりも下側にあるか、すなわち、図１４の領域Ａ２に属するか否か（以下、プロ条件２）を判定する。そして、プロ条件２が満たされる場合には、５５４０〜５６２３ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ４２でプロ条件２が満たされない場合には、ステップＳ４３に進む。ステップＳ４３では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１４に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ３の上側かつ直線Ｌ２よりも下側にあるか、すなわち、図１４の領域Ａ３に属するか否か（以下、プロ条件３）を判定する。そして、プロ条件３が満たされる場合には、５４８０〜５５４０ｋｇ・ｃｍ²が最適スイングＭ帯Ｉであると判定する。一方、ステップＳ４３でプロ条件３が満たされない場合には、ステップＳ４４に進む。ステップＳ４４では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１４に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ４の上側かつ直線Ｌ３よりも下側にあるか、すなわち、図１４の領域Ａ４に属するか否か（以下、プロ条件４）を判定する。そして、プロ条件４が満たされる場合には、５４００〜５４８０ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ４４でプロ条件４が満たされない場合、すなわち、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１４に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ４よりも下側にあり、図１４の領域Ａ５に属する場合には、５４００ｋｇ・ｃｍ²以下が最適スイングＭＩ帯であると判定する。

【0111】

一方、ステップＳ５１及びこれに続くステップＳ５２，Ｓ５３は、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hの大きさに応じて、最適スイングＭＩ帯を決定するステップである。ここでは、Ｐ_{1_AVE}，Ｐ_{2_AVE}及びＶ_hの値が大きい程、最適スイングＭＩ帯は段階的に大きな値に設定される。

【0112】

具体的には、決定部２４Ｅは、ステップＳ５１において、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１５に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ５の上側にあるか、すなわち、図１５の領域Ｂ１に属するか否か（以下、アベレージ条件１）を判定する。そして、アベレージ条件１が満たされる場合には、５３００ｋｇ・ｃｍ²以上が最適スイングＭＩ帯であると判定する。一方、ステップＳ５１でアベレージ条件１が満たされない場合には、ステップＳ５２に進む。ステップＳ５２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１５に示すＰ_{1_AVE}−Ｐ_{2_AVE}平面において直線Ｌ６の上側かつ直線Ｌ５よりも下側にあるか、すなわち、図１５の領域Ｂ２に属するか否か（以下、アベレージ条件２）を判定する。そして、アベレージ条件２が満たされる場合には、５２００〜５３００ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ５２でアベレージ条件２が満たされない場合には、ステップＳ５３に進む。ステップＳ５３では、決定部２４Ｅは、第１指標算出工程で算出されたヘッド速度Ｖ_hが３５．５ｍ／ｓより大きいか否か（アベレージ条件３）を判定する。そして、アベレージ条件３が満たされる場合には、５２００〜５３００ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ５３でアベレージ条件３が満たされない場合には、すなわち、第１指標算出工程で算出された（Ｐ_{1_AVE}，Ｐ_{2_AVE}）で表される点が、図１５に示す直線Ｌ６よりも下側にあり、図１４の領域Ｂ３に属する場合には、５２００ｋｇ・ｃｍ²以下が最適スイングＭＩ帯であると判定する。

【0113】

以上のステップＳ４１〜Ｓ４４及びステップＳ５１〜Ｓ５３は、以下の知見に基づく。すなわち、図１４及び図１５を参照して説明した上記実験では、ヘッド速度Ｖ_h及び、飛距離を最大にする最適スイングＭＩも算出された。ヘッド速度Ｖ_hは、上記と同様の工程に従って算出された。一方、最適スイングＭＩについては、ゴルファーに様々なスイング慣性モーメントＩ_sのゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブのスイング慣性モーメントＩ_sを特定し、これを最適スイングＭＩとした。より具体的には、プロモデルユーザには、スイング慣性モーメントＩ_sが５６５０ｋｇ・ｃｍ²、５６１０ｋｇ・ｃｍ²、５５２０ｋｇ・ｃｍ²、５４４０ｋｇ・ｃｍ²、５３６０ｋｇ・ｃｍ²の５種類のゴルフクラブをスイングさせた。一方、アベレージモデルユーザには、スイング慣性モーメントＩ_sが５３３０ｋｇ・ｃｍ²、５２３０ｋｇ・ｃｍ²、５１３０ｋｇ・ｃｍ²の３種類のゴルフクラブをスイングさせた。図１４及び図１５には、本実験により得られたヘッド速度Ｖ_h[ｍ／ｓ]及び最適スイングＭＩの値も示されている。

【0114】

図１４及び図１５に示すとおり、上記の実験からは、最適スイングＭＩは、腕出力パワーＰ_{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}平面において負の傾きを持つ直線である。ステップＳ４１〜Ｓ４４及びステップＳ５１〜Ｓ５３では、この記憶部２３内の対応関係データ２８が参照され、上記の判定が行われる。なお、図２では、対応関係データ２８は、フィッティングプログラム３とは別のデータとして示されているが、プログラム３内に組み込まれていてもよい。

【0115】

ところで、アベレージモデル領域では、腕出力パワーＰ_{1_AVE}及びクラブ入力パワーＰ_{2_AVE}を示す点が領域Ｂ３に属する場合には、ヘッド速度Ｖ_hが３５．５ｍ／ｓより大きいか否かにより、異なる最適スイングＭＩ帯に振り分けられるようになっている。これは、図１５に示す実験結果に示されるとおり、これらの領域においては、Ｖ_h＞３５．５ｍ／ｓか否かで最適スイングＭＩが異なる値になっているためである。

【0116】

ところで、図１４及び図１５には、上述したとおり、実際にプロモデルユーザ２１名及びアベレージモデルユーザ２２名にテストクラブを試打させた実験の結果がプロットされている。そして、プロモデルユーザに関しては、２１個のデータうち、Ｆのマークを付している１名を除く２０名において、最適スイングＭＩが上述の最適振り易さ決定工程で決定される最適スイングＭＩ帯に属する結果となった。すなわち、正答率９５％以上でフィッティングができていることが確認された。また、アベレージモデルユーザに関しては、２２個のデータうち、Ｆのマークを付している１名を除く２１名において、最適スイングＭＩが上述の最適振り易さ決定工程で決定される最適スイングＭＩ帯に属する結果となった。すなわち、正答率９５％以上でフィッティングができていることが確認された。

【0117】

＜１−２−７．第２指標算出工程＞
以下、計測工程で得られた計測データに基づいて、第２スイング指標を算出する第２指標算出工程（Ｓ７）について説明する。第２スイング指標とは、最適剛性指標を決定するための、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態では、第２スイング指標として、後述される第１〜第４特徴量Ｆ₁〜Ｆ₄が算出される。

【0118】

第１〜第４特徴量Ｆ₁〜Ｆ₄について理解するためには、まず、最適剛性指標について理解することが重要である。最適剛性指標とは、ゴルファー７に適したシャフト４０の剛性を示す指標のことであり、本実施形態では、シャフト４０の剛性は、シャフト４０の複数の位置における曲げ剛性の分布（以下、ＥＩ分布）として評価される。本実施形態に係るＥＩ分布は、定量的に数値を用いて表現され、より具体的には、インターナショナル・フレックス・コード（ＩＦＣ）を用いて算出される。そのため、まず、このＩＦＣについて説明する。なお、ＩＦＣは、本出願人により広く提案されているシャフトの特性を示す公知の指標であり、例えば、特許文献１をはじめとして、既に様々な文献で詳しく説明されている。従って、ここで改めて説明する必要は必ずしもないが、参考のため、ここでも説明を行う。

【0119】

ＩＦＣは、図１６に示すとおり、シャフト４０の延びる方向に沿った４つの位置Ｈ１〜Ｈ４におけるシャフト４０の曲げ剛性をそれぞれ０〜９の１桁の数値で表し、この４つの数値をシャフト４０の延びる方向に沿って配列したコードである。より具体的には、シャフト４０のバット端からチップ端に向かってこの順に概ね一定間隔で、４つの測定点Ｈ１〜Ｈ４が定義される。例えば、シャフト４０のチップ端から３６インチの箇所を測定点Ｈ１とし、２６インチの箇所を測定点Ｈ２とし、１６インチの箇所を測定点Ｈ３とし、６インチの箇所を測定点Ｈ４とすることができる。そして、これらの４つの測定点Ｈ１〜Ｈ４のそれぞれにおける曲げ剛性の値（以下、ＥＩ値）Ｊ₁〜Ｊ₄が計測される。

【0120】

シャフト４０の各測定点Ｈ（Ｈ１〜Ｈ４）におけるＥＩ値（Ｎ・ｍ²）は、様々な方法で測定することができ、例えば、インテスコ社製の２０２０型計測機（最大荷重５００ｋｇｆ）を用いて図１７に示すようにして測定することができる。この測定方法では、２つの支持点１１１，１１２でシャフト４０を下方から支持しつつ、測定点Ｈに上方から荷重Ｆを加えたときのたわみ量を測定する。支持点１１１と支持点１１２との間の距離（スパン）は、例えば、２００ｍｍとすることができ、測定点Ｈは、支持点１１１と支持点１１２の中間点とすることができる。より具体的には、支持点１１１，１１２を支える支持体１１４，１１５を固定した状態で、測定点Ｈにおいて圧子１１３を一定速度（例えば、５ｍｍ／分）で下方へ移動させる。そして、荷重Ｆが２０ｋｇｆに達した時点で圧子１１３の移動を終了させ、この瞬間のシャフト４０のたわみ量（ｍｍ）を測定し、このたわみ量をＥＩ値（Ｎ・ｍ²）に換算する。

【0121】

次に、以上の４つの測定点Ｈ１〜Ｈ４におけるＥＩ値Ｊ₁〜Ｊ₄を、それぞれ１０段階のランク値Ｋ₁〜Ｋ₄に変換する。具体的には、ランク値Ｋ₁〜Ｋ₄は、それぞれ測定点Ｈ１〜Ｈ４用の以下の変換表（表１〜表４）に従って、ＥＩ値Ｊ₁〜Ｊ₄から算出することができる（表１〜４中、変換後のランク値をＩＦＣの欄に示している）。そして、このようにして測定点Ｈ１〜Ｈ４にそれぞれ付与された４つのランク値Ｋ₁〜Ｋ₄を、よりバット側に対応する値がより左に、よりチップ側に対応する値がより右にくるように配列する。こうして得られた４桁のコードが、ＩＦＣである。ＩＦＣでは、各桁の数値が大きい程、対応する位置での剛性が高いことを意味する。

【表1】

【表2】

【表3】

【表4】

【0122】

第２指標算出工程では、算出部２４Ｄにより、第１〜第４特徴量Ｆ₁〜Ｆ₄が算出される。本実施形態では、第１〜第４特徴量Ｆ₁〜Ｆ₄は、それぞれゴルファー７に適したＥＩ値Ｊ₁〜Ｊ₄である最適ＥＩ値Ｊ_S1〜Ｊ_S4、ひいてはゴルファー７に適したランク値Ｋ₁〜Ｋ₄である最適ランク値Ｋ_S1〜Ｋ_S4を決定するための指標である。そのため、本実施形態では、第１〜第４特徴量Ｆ₁〜Ｆ₄としては、それぞれ最適ＥＩ値Ｊ_S1〜Ｊ_S4と相関を有する特徴量が選択される。また、本実施形態では、第１〜第４特徴量Ｆ₁〜Ｆ₄として以下の指標が用いられるが、第２スイング指標としては、最適剛性指標との相関が認められる限り、その他の任意の特徴量を用いることができる。

【0123】

第１特徴量Ｆ₁は、トップ付近のコック方向の角速度ω_yの傾きであり、例えばトップから５０ｍｓ前の角速度ω_yと、トップから５０ｍｓ後の角速度ω_yとの和で表すことができる。

【0124】

第２特徴量Ｆ₂は、トップから、角速度ω_yが最大となる時点までの当該角速度ω_yの平均値である。第２特徴量Ｆ₂は、まず、トップからインパクトまでの間で角速度ω_yが最大となる時点を求め、トップからこの時点までの角速度ω_yの累積値を、トップからこの時点までの時間で除することにより算出される。

【0125】

第３特徴量Ｆ₃は、角速度ω_yが最大となる時点からインパクトまでの当該角速度ω_yの平均値である。第３特徴量Ｆ₃は、角速度ω_yが最大となる時点からインパクトまでの角速度ω_yの累積値を、角速度ω_yが最大となる時点からインパクトまでの時間で除することにより算出される。

【0126】

第４特徴量Ｆ₄は、トップからインパクトまでの角速度ω_yの平均値であり、トップからインパクトまでの角速度ω_yの累積値を、トップからインパクトまでの時間で除することにより算出される。

【0127】

ところで、スイング動作中、ゴルフクラブ４のシャフト４０は、当該ゴルフクラブの先端に比較的重量が大きいヘッド４１が存在するため、その慣性により曲げが生じる。この曲げは、スイングの全過程において、シャフト４０の同一箇所に生じるのではなく、図１８に示されるように、トップからインパクトに向けてシャフト４０の手元側から先端側に伝わる。換言すれば、トップからインパクトに向けてスイングが進行するにしたがい、シャフト４０における曲げの位置が当該シャフト４０の手元側から先端側に移動する。

【0128】

より具体的には、アドレスからテイクバックを行い、トップに至った時点（図１８において（１）で示される時点）では、シャフト４０の手元付近に曲げが生じる。ついで、切り返しを行い、ダウンスイング初期（図１８において（２）で示される時点）に至ると、曲げはシャフト４０の先端側にやや移動する。さらに、ゴルファー７の腕が水平になる時点（図１８において（３）で示される時点）では、曲げはシャフト４０中央よりも先端側に移動する。そして、インパクト直前（図１８において（４）で示される時点）では、曲げはシャフト４０の先端付近まで移動する。

【0129】

従って、第１〜第４特徴量Ｆ₁〜Ｆ₄は、それぞれスイング動作中のトップ付近からインパクト付近までの間の第１〜第４区間において算出することができる。また、ここでの第１〜第３区間は、この順に時間経過に沿った区間であり、互いに一部重複する又は重複することのない区間となっている。

【0130】

＜１−２−８．最適剛性決定工程＞
次に、最適剛性決定工程（Ｓ８）の流れについて説明する。本工程では、決定部２４Ｅが、第２スイング指標（第１〜第４特徴量Ｆ₁〜Ｆ₄）と最適剛性指標（最適ＥＩ値Ｊ_S1〜Ｊ_S4）との相関関係を表す予め定められた近似式に従って、最適剛性指標（最適ＥＩ値Ｊ_S1〜Ｊ_S4）を決定する。本実施形態に係る近似式は、線形近似式であり、以下のように表される。
Ｊ_S1＝ａ₁・Ｆ₁＋ｂ₁
Ｊ_S2＝ａ₂・Ｆ₂＋ｂ₂
Ｊ_S3＝ａ₃・Ｆ₃＋ｂ₃
Ｊ_S4＝ａ₄・Ｆ₄＋ｂ₄

【0131】

決定部２４Ｅは、第２指標算出工程で算出された第１〜第４特徴量Ｆ₁〜Ｆ₄をかかる近似式に代入することにより、最適ＥＩ値Ｊ_S1〜Ｊ_S4を算出する。また、決定部２４Ｅは、上述の表１〜表４の換算表に従って、最適ＥＩ値Ｊ_S1〜Ｊ_S4をそれぞれ最適ランク値Ｋ_S1〜Ｋ_S4に変換する。

【0132】

なお、上式中、ａ１〜ａ４及びｂ１〜ｂ４は、予め行われた多数の実験結果から回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。ここでいう実験とは、例えば、特許文献１と同様に、以下のように行うことができる。すなわち、まず、多数のゴルファーの各々に複数のゴルフクラブを振ってもらい、そのときの飛距離、打球の方向性（左右ずれ）及び官能試験による振り易さを数値化する。そして、その数値から各ゴルファーに適しているゴルフクラブを決定し、当該ゴルフクラブのＥＩ値を当該ゴルファーの最適ＥＩ値とする。また、上記と同様の方法で、各ゴルファーの第１〜第４特徴量Ｆ₁〜Ｆ₄を算出する。そして、かかる実験の後、多数のゴルファー分の第１〜第４特徴量Ｆ₁〜Ｆ₄及び最適ＥＩ値のデータを回帰分析することで、ａ１〜ａ４及びｂ１〜ｂ４が算出される。

【0133】

また、より信頼性の高い近似式とするために、ａ１〜ａ４及びｂ１〜ｂ４の値を条件に応じて変更することができる。例えば、ヘッド速度Ｖ_hに応じて、近似式を用意することができる。一例としては、上記実験データを、ヘッド速度帯に応じて分類し（例えば、４５ｍ／ｓ以上、４１〜４５ｍ／ｓ、４１ｍ／ｓ以下）、同じ分類に属するデータのみを対象に上記近似式を作成し、ａ１〜ａ４及びｂ１〜ｂ４を決定しておくことができる。そして、最適剛性決定工程では、ゴルファー７のヘッド速度Ｖ_hがどのヘッド速度帯に属するかを判断し、当該ヘッド速度帯に対応する近似式を用いて、最適剛性指標を算出する。

【0134】

＜１−２−９．最適クラブ選択工程＞
以上の工程Ｓ１〜Ｓ８により、最適振り易さ指標（最適スイングＭＩ帯）及び最適剛性指標（最適ＥＩ値Ｊ_S1〜Ｊ_S4、最適ランク値Ｋ_S1〜Ｋ_S4）が決定されると、選択部２４Ｆは、最適クラブ選択工程（Ｓ９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、ゴルファー７に適したシャフト４０(以下、推奨シャフト)が特定される。また、フィッティングの結果として、ゴルファー７に適したゴルフクラブ４（以下、推奨ゴルフクラブ）も特定される。

【0135】

本実施形態では、まず、選択部２４Ｆは、推奨ゴルフクラブに用いられるべきヘッド４１(以下、推奨ヘッド)の種類を決定する。推奨ヘッドの種類の決定は、本明細書では説明されないフィッティング処理により行うこともできるし、表示部２１及び入力部２２を介してユーザに質問する等してお好みのヘッド４１を選択させることにより行うこともできる。そして、選択部２４Ｆは、推奨ヘッドのスペックを示す情報をヘッドＤＢ２７内から読み出すととともに、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。ヘッドＤＢ２７内に登録されているヘッド４１のスペックを示す情報には、製造メーカー、型番及び重量等が含まれる。一方、シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報には、製造メーカー、型番、４つの位置Ｈ１〜Ｈ４におけるＥＩ値Ｊ₁〜Ｊ₄及びランク値Ｋ₁〜Ｋ₄（ＩＦＣ）、シャフト４０の重量、フレックス、トルク、調子、長さ及び重心位置等が含まれる。選択部２４Ｆは、これらの情報から、各シャフト４０と推奨ヘッドとを組み合わせた場合のゴルフクラブ４のスイング慣性モーメントＩ_Sを算出し、その値が最適スイングＭＩ帯に属することとなるゴルフクラブ４（以下、第１絞り込みゴルフクラブ）及びそれに含まれるシャフト４０（以下、第１絞り込みシャフト）を特定する。なお、第１絞り込みゴルフクラブ及びシャフトは、通常多数本存在する。

【0136】

続いて、選択部２４Ｆは、各第１絞り込みシャフトについて、当該シャフトのランク値Ｋ₁〜Ｋ₄と、最適剛性決定工程で決定された最適ランク値Ｋ_S1〜Ｋ_S4との一致度を算出し、推奨シャフトとして、一致度の高いシャフトを特定する。一致度は、例えば、以下の数２５の式に従って算出することができ、値が小さいほど一致度が高い。

【数25】

【0137】

なお、推奨シャフトとしては、１本のみ特定されてもよいし、複数本が特定されてもよい。また、推奨シャフトとしては、第１絞り込みシャフトの中で相対的に一致度が高い所定数のシャフトが特定されてもよいし、一定以上の一致度を有する全てのシャフトが特定されてもよい。表示制御部２４Ｇは、シャフトＤＢ２９を参照して、推奨シャフトの種類を示す情報とともに、推奨シャフトのスペックを示す情報（ＩＦＣの値を含む）を表示部２１上に表示させる。さらに、表示制御部２４Ｇは、ゴルファーの最適スイングＭＩ帯及び最適ランク値Ｋ_S1〜Ｋ_S4を表示部２１上に表示させる。また、表示制御部２４Ｇは、推奨ゴルフクラブが推奨シャフトと推奨ヘッドとを組み合わせたゴルフクラブであることを示す情報を、表示部２１上に表示させる。これにより、ユーザは、ゴルファー７に適したゴルフクラブ４及びシャフト４０の種類を知ることができるとともに、当該ゴルフクラブのＩＦＣの値、最適スイングＭＩ帯及び最適ランク値Ｋ_S1〜Ｋ_S4を知ることができる。

【0138】

＜２．第２実施形態＞
図１９に、第２実施形態に係るフィッティングシステム２００の全体構成を示す。フィッティングシステム２００は、第１実施形態に係るフィッティングシステム１００と多くの点で共通する。従って、以下では、簡単のため、第１実施形態と第２実施形態の差異点を中心に説明し、第１実施形態と同様の構成には同様の符号を付して、説明を省略する。

【0139】

図２０は、第２実施形態に係るフィッティング処理の流れを示すフローチャートである。図４と図２０とを比較すれば明らかなとおり、第１実施形態と第２実施形態の解析処理の相違点は、主として最適振り易さ決定工程（Ｓ６）に代えて、最適特性決定工程（Ｓ２０６）が実行され、最適剛性決定工程（Ｓ８）に代えて、最適剛性決定工程（Ｓ２０８）が実行され、最適クラブ選択工程（Ｓ９）に代えて、最適クラブ選択工程（Ｓ２０９）が実行される点にある。より具体的には、第１実施形態に係る最適振り易さ決定工程（Ｓ６）では、最終的に最適スイングＭＩ帯が決定されたが、第２実施形態に係る最適特性決定工程（Ｓ２０６）では、最終的にゴルファーに適したゴルフクラブ４のシャフト４０の重量（最適シャフト重量）が決定される。すなわち、第１実施形態では、最適振り易さ指標が決定されたが、第２実施形態では、ゴルファーに適したゴルフクラブ４の特定の部位の特性を表す最適特性指標が決定される。また、第１実施形態に係る最適剛性決定工程（Ｓ８）では、最終的に最適剛性指標として、４桁のＩＦＣのコード（最適ランク値Ｋ_S1〜Ｋ_S4）が決定されたが、第２実施形態に係る最適剛性決定工程（Ｓ２０８）では、最終的に最適剛性指標として、ゴルファーに適したゴルフクラブ４のシャフト４０のフレックス（最適フレックス）が決定される。また、第１実施形態に係る最適クラブ選択工程（Ｓ９）では、ヘッド４１を固定した上で、ゴルフクラブ４全体での最適振り易さ指標（最適スイングＭＩ帯）を達成することのできるシャフト４０を絞り込んだ後、最適剛性指標（最適ランク値Ｋ_S1〜Ｋ_S4）にできる限り合致するシャフト４０が決定されたが、第２実施形態に係る最適クラブ選択工程（Ｓ２０９）では、ヘッド４１を固定した上で、最適特性指標（最適シャフト重量）及び最適剛性指標（最適フレックス）にできる限り合致するシャフト４０が決定される。

【0140】

図１９に示すとおり、フィッティングシステム２００は、フィッティング装置２に代えて、フィッティング装置２０２を備える。フィッティング装置２０２は、フィッティング装置２と同様のハードウェア構成を有するが、フィッティング装置２０２では、フィッティングプログラム３に代えて、フィッティングプログラム２０３がインストールされている。そのため、制御部２４は、記憶部２３内のフィッティングプログラム２０３を読み出して実行することにより、仮想的に取得部２４Ａ、グリップ挙動導出部２４Ｂ、肩挙動導出部２４Ｃ、算出部２４Ｄとして動作する他、決定部２２４Ｅ、選択部２２４Ｆ及び表示制御部２２４Ｇとして動作することができる。決定部２２４Ｅは、第１実施形態との相違点である最適特性決定工程（Ｓ２０６）及び最適剛性決定工程（Ｓ２０８）を実行する仮想ユニットであり、選択部２２４Ｆ及び表示制御部２２５Ｇは、第１実施形態との相違点である最適クラブ選択工程（Ｓ２０９）を実行する仮想ユニットである。また、最適特性決定工程（Ｓ２０６）を実行できるよう、フィッティング装置２０２の記憶部２３には、対応関係データ２８に代えて、対応関係データ２２８が格納されている。対応関係データ２２８は、最適シャフト重量を決定するための条件を示すデータである。

【0141】

第２実施形態でも、第１実施形態と同じく、計測工程（Ｓ１）、第１変換工程（Ｓ２）、第２変換工程（Ｓ３）、肩挙動導出工程（Ｓ４）及び第１指標算出工程（Ｓ５）が順次実行され、続いて、最適特性決定工程（Ｓ２０６）が実行される。そしてその後、第２指標算出工程（Ｓ７）が実行され、続いて、最適剛性決定工程（Ｓ２０８）及び最適クラブ選択工程（Ｓ２０９）が順次実行される。以下では、第１実施形態との相違点である最適特性決定工程（Ｓ２０６）、最適剛性決定工程（Ｓ２０８）及び最適クラブ選択工程（Ｓ２０９）について、順に説明する。

【0142】

＜２−１．最適特性決定工程＞
最適特性決定工程（Ｓ２０６）では、決定部２２４Ｅは、第１スイング指標（腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_h）の大きさに応じて、最適シャフト重量の範囲（以下、最適シャフト重量帯）を決定する。このとき、Ｐ_{1_AVE}，Ｐ_{2_AVE}及びＶ_hの値が大きい程、最適シャフト重量帯は段階的に大きな値に設定される。具体的には、図２１に示すフローチャートに従って、最適特性決定工程は進行する。

【0143】

ステップＳ１０１では、決定部２２４Ｅは、ユーザの入力に基づいて、ユーザのお好みのヘッド４１（以下、嗜好ヘッド）の種類を決定する。嗜好ヘッドとは、推奨ゴルフクラブに含まれるべきヘッド４１である。嗜好ヘッドの種類の決定は、表示部２１及び入力部２２を介して、ユーザから嗜好ヘッドの種類を直接的に指定させることにより行うこともできるし、ユーザに各種質問を行うとともにその回答を受け付け、当該回答から所定のアルゴリズムに従って嗜好ヘッドを絞り込むことにより行うこともできる。

【0144】

続くステップＳ１０２では、決定部２２４Ｅは、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hの大きさに応じて、最適振り易さ指標である最適スイングＭＩを決定する。最適スイングＭＩは、第１実施形態と同様に、図１３に示すフローチャートに従って算出することも可能であるが、本実施形態では、以下の重回帰式に従って算出される。
（最適スイングＭＩ）＝ｅ₁・Ｐ_{1_AVE}＋ｅ₂・Ｐ_{2_AVE}＋ｅ₃・Ｖ_h＋ｅ₄

【0145】

上記重回帰式は、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hを説明変数とし、最適スイングＭＩを目的変数とする重回帰式である。係数ｅ₁〜ｅ₄は、予め実験により変数Ｐ_{1_AVE}，Ｐ_{2_AVE,}Ｖ_h及び最適スイングＭＩの多数のデータセットを取得し、これらのデータセットに基づいて重回帰分析を行うことにより算出され、対応関係データ２２８として記憶部２３内に格納されている。すなわち、対応関係データ２２８は、Ｐ_{1_AVE},Ｐ₂v__AVE及びＶ_hの大きさと、最適スイングＭＩとの対応関係を定めるデータを含む。なお、上記重回帰式は、テストクラブのモデル毎に用意しておくこともできる。この場合、第１実施形態と同様に、解析の精度が向上し得る。また、図１９では、対応関係データ２２８は、フィッティングプログラム２０３とは別のデータとして示されているが、プログラム２０３内に組み込まれていてもよい。

【0146】

続くステップＳ１０３では、決定部２２４Ｅは、最適スイングＭＩの大きさ及び嗜好ヘッドの種類に応じて、最適シャフト重量帯を決定する。具体的には、以下の表５及び表６に従って、最適スイングＭＩを最適シャフト重量帯に変換する。表５は、ヘッドＡ用のデータであり、表６は、ヘッドＡとは種類の異なるヘッドＢ用のデータである。

【表5】

【表6】

【0147】

すなわち、決定部２２４Ｅは、嗜好ヘッドがヘッドＡである場合には、表５に最適スイングＭＩを照合し、最適シャフト重量帯を決定する。一方、嗜好ヘッドがヘッドＢである場合には、表６に最適スイングＭＩを照合し、最適シャフト重量帯を決定する。なお、本実施形態では、ステップＳ１０１において、嗜好ヘッドとしてヘッドＡ又はヘッドＢのいずれかが選択されるものとする。

【0148】

表５及び表６のデータは、最適スイングＭＩと最適シャフト重量帯との対応関係を定めるデータであり、対応関係データ２２８に含まれる。つまり、対応関係データ２２８とは、第１スイング指標（Ｐ_{1_AVE},Ｐ_{2_AVE}及びＶ_h）と最適特性指標（最適シャフト重量帯）との対応関係を定めるデータであると言うこともできる。最適スイングＭＩと最適シャフト重量帯との対応関係を定める対応関係データ２２８は、ヘッド４１の種類毎に用意することができ、本実施形態では、同じモデルで重量の異なる２つのヘッド（ヘッドＡ及びヘッドＢ）に対しそれぞれ用意されている。ただし、最適スイングＭＩと最適シャフト重量帯との対応関係を定める対応関係データ２２８としては、ヘッド４１の種類に依存しないデータを用意することもできる。

【0149】

ところで、ゴルフクラブ４のスイング慣性モーメントＩ_Sとは、ヘッド４１のスイング慣性モーメントと、シャフト４０のスイング慣性モーメントと、その他の部位（グリップ４２、フェラル等）のスイング慣性モーメントとの合計値である。従って、嗜好ヘッドとしてスイング慣性モーメントの小さい（又は大きい）ものが選択される場合には、その分だけ、推奨ゴルフクラブ４に含まれるべき推奨シャフトはスイング慣性モーメントの大きい（又は小さい）ものが選択されることが適切である。この点、本実施形態のように、ヘッド４１の種類毎に最適シャフト重量帯を決定するための対応関係データ２２８が用意される場合には、嗜好ヘッドに応じてシャフト４０が選択されることになり、ゴルフクラブ４全体でのスイング慣性モーメントＩ_Sを最適化することができる。以上のとおり、本実施形態では、第１実施形態に引き続き、ヘッド４１の種類を決定した上で、ゴルフクラブ４全体での最適振り易さ指標を達成できるような最適シャフト重量帯が決定される。

【0150】

また、本発明者らは、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−７４５に様々な重量のシャフトを取り付け、これらのゴルフクラブを３７名のテスターに試打させた。そして、各テスターについて、飛距離を最大にするシャフトの重量を決定し、実験による最適シャフト重量とした。また、３７名のテスターの各々に対し、上述した最適特性決定工程と同じ方法で最適シャフト重量帯を導出した。そして、これらの結果を比較したところ、３７名中３３名において、実験による最適シャフト重量が最適振り易さ決定工程で決定された最適シャフト重量帯に属する結果となった。すなわち、約９０％の正答率でフィッティングができていることが確認された。

【0151】

＜２−２．最適剛性決定工程＞
次に、最適剛性決定工程（Ｓ２０８）の流れについて説明する。本工程では、決定部２２４Ｅが、第２スイング指標（第１〜第４特徴量Ｆ₁〜Ｆ₄）の大きさに応じて、最適剛性指標（最適フレックス）を決定する。

【0152】

具体的には、まず、決定部２２４Ｅは、第１〜第４特徴量Ｆ₁〜Ｆ₄に基づいて、最適ＥＩ値Ｊ_S1〜Ｊ_S4及び最適ランク値Ｋ_S1〜Ｋ_S4を決定する。これらの値Ｊ_S1〜Ｊ_S4，Ｋ_S1〜Ｋ_S4の決定方法は、第１実施形態と同様とすることができる。

【0153】

続いて、決定部２２４Ｅは、最適ランク値Ｋ_S1〜Ｋ_S4に基づいて、最適フレックスを決定する。フレックスとは、シャフト４０全体での硬さ（曲げ剛性）を評価する指標である。従って、シャフト４０の複数の位置におけるゴルファー７に適した曲げ剛性を表す最適ランク値Ｋ_S1〜Ｋ_S4が分かれば、これらの値に基づいて最適フレックスを算出することができる。例えば、特定の位置での最適ランク値を最適フレックスとすることもできるし、複数の位置での最適ランク値の平均値を最適フレックスをとすることもできる。

【0154】

なお、他の実施形態では、最適フレックスを最適ランク値Ｋ_S1〜Ｋ_S4ではなく、最適ＥＩ値Ｊ_S1〜Ｊ_S4から直接算出することもできるし、これらの値Ｊ_S1〜Ｊ_S4，Ｋ_S1〜Ｋ_S4に基づかずに算出することもきる。後者の場合、最適フレックスを特定可能な適切な特徴量を第２スイング指標決定工程で算出しておけばよい。

【0155】

＜２−３．最適クラブ選択工程＞
以上の工程により、最適特性指標（最適シャフト重量帯）及び最適剛性指標（最適フレックス）が決定されると、選択部２２４Ｆは、最適クラブ選択工程（Ｓ２０９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、推奨シャフトが特定される。また、推奨ゴルフクラブも特定される。

【0156】

本実施形態では、まず、選択部２２４Ｆは、シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報を参照して、最適剛性決定工程で決定された最適フレックスと同じフレックスを有するシャフト４０（以下、第１絞り込みシャフト）を特定する。第１絞り込みシャフトは、通常多数本存在する。続いて、選択部２２４Ｆは、第１絞り込みシャフトの中から、最適特性決定工程で決定された最適シャフト重量帯に属するシャフト４０を特定し、推奨シャフトとして選択する。このとき、第１絞り込みシャフトの中に最適シャフト重量帯に属するシャフト４０が存在しない場合には、第１絞り込みシャフトのうち最適シャフト重量帯に最も近いシャフト４０が、推奨シャフトとして選択される。推奨シャフトは、１本のみ特定されてもよいし、複数本が特定されてもよい。さらに、選択部２２４Ｆは、推奨シャフトと嗜好ヘッドを備えるゴルフクラブ４を、推奨ゴルフクラブとして選択する。なお、本実施形態では、最適シャフト重量帯よりも最適フレックスを優先して推奨シャフトが決定されるが、逆に、最適フレックスよりも最適シャフト重量帯を優先して推奨シャフトが決定されてもよい。

【0157】

表示制御部２２４Ｇは、シャフトＤＢ２９を参照して、推奨シャフトの種類を示す情報とともに、推奨シャフトのスペックを示す情報（ＩＦＣ、フレックス、トルク、調子及び重量の値を含む）を表示部２１上に表示させる。さらに、表示制御部２２４Ｇは、ゴルファーの最適スイングＭＩ帯及び最適ランク値Ｋ_S1〜Ｋ_S4を表示部２１上に表示させる。また、表示制御部２２４Ｇは、推奨ゴルフクラブが推奨シャフトと嗜好ヘッドとを組み合わせたゴルフクラブであることを示す情報を、表示部２１上に表示させる。これにより、ユーザは、推奨ゴルフクラブ及び推奨シャフトの詳細を知ることができる。

【0158】

＜３．変形例＞
以上、本発明の幾つかの実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。

【0159】

＜３−１＞
上記実施形態では、ゴルファー７のスイング動作を計測する計測機器として、加速度センサ、角速度センサ及び地磁気センサの３つを有するセンサユニット１が使用されたが、計測機器を他の構成とすることもできる。例えば、地磁気センサを省略することもできる。この場合には、統計的手法により、ｘｙｚ局所座標系からＸＹＺ全体座標系へと計測データを変換することが可能である。なお、このような手法については、公知技術であるため（要すれば、特開２０１３−５６０７４号公報参照）、ここでは詳細な説明を省略する。或いは、計測機器として、三次元計測カメラを使用することもできる。三次元計測カメラにより、ゴルファーやゴルフクラブ、ゴルフボールの挙動を計測する手法についても、公知であるため、ここでは詳細な説明を省略する。なお、三次元計測カメラを用いた場合には、計測データのｘｙｚ局所座標系からＸＹＺ全体座標系への変換工程を省略することもでき、直接的にＸＹＺ全体座標系でのグリップの挙動を計測することができる。

【0160】

＜３−２＞
上記実施形態では、最適振り易さ指標として、最適スイングＭＩを算出することが例示された。しかしながら、ゴルフクラブ４の振り易さを表すその他の様々な指標について、最適振り易さ指標を算出してもよい。例えば、ゴルフクラブ４の重量ｍ₂（厳密には、質量をｍ₂とするならば重量はｍ₂ｇであるが、簡単のためここではいずれもｍ₂で表すこととする）、グリップエンド周りの慣性モーメントＩ_G、及び重心周りの慣性モーメントＩ₂も、ゴルフクラブ４の振り易さを表す指標であり、これらの指標に対して最適振り易さ指標を算出してもよい。また、複数の最適振り易さ指標を算出して、これら全ての最適振り易さ指標に基づいてフィッティングを行ってもよい。例えば、最適スイングＭＩ、ゴルファーに適したグリップエンド周りの慣性モーメントＩ_G（以下、最適グリップエンドＭＩ）、ゴルファーに適した重心周りの慣性モーメントＩ₂（以下、最適重心ＭＩ）及びゴル
ファーに適したゴルフクラブ４の重量ｍ₂（以下、最適クラブ重量）の全てを算出し、これら４つの条件に合致するゴルフクラブ４のシャフト４０をシャフトＤＢ２９内から検索してもよい。また、これら４つの条件に合致する最適特性指標を決定してもよい。

【0161】

なお、グリップエンド周りの慣性モーメントＩ_Gとは、例えば、以下の式に従って算出
される。Ｉ_Gを決定するパラメータであるｍ₂，Ｉ₂，Ｌは、ゴルフクラブ４の諸元である
ため、Ｉ_Gもゴルフクラブ４の諸元となる。
Ｉ_G＝Ｉ₂＋ｍ₂Ｌ²

【0162】

また、ゴルフクラブ４の重量ｍ₂、グリップエンド周りの慣性モーメントＩ_G、及び重心周りの慣性モーメントＩ₂は、スイング慣性モーメントＩ_Sと概ね比例関係にある。例えば、図２２は、様々なゴルフクラブのＩ_S，Ｉ_Gの関係をプロットしたグラフであり、図２３は、様々なゴルフクラブのＩ_S，ｍ₂の関係をプロットしたグラフである。従って、上記実施形態に係る最適振り易さ決定工程において、最適スイングＭＩを最適グリップエンドＭＩ、最適重心ＭＩ又は最適クラブ重量に変換すれば、Ｉ_G、Ｉ₂、ｍ₂についての最適振り易さを決定することができる。例えば、図２４及び図２５は、図１３のフローチャートを改変したものであり、図２５は、最適クラブ重量を決定するためのフローチャートであり、図２４は、最適グリップエンドＭＩを決定するためのフローチャートである。

【0163】

＜３−３＞
上記実施形態では、最適振り易さ指標を決定するための第１スイング指標として、腕出力パワーＰ_{1_AVE}、クラブ入力パワーＰ_{2_AVE}及びヘッド速度Ｖ_hの３つが組み合わせて用いられた。しかしながら、この例に限られず、例えば、スイング指標としてＰ_{1_AVE}だけ、Ｐ_{2_AVE}だけを用いてもよいし、（Ｐ_{1_AVE}，Ｖ_h）、（Ｐ_{2_AVE}，Ｖ_h）、（Ｐ_{1_AVE}，Ｐ_{2_AVE}）のような２つの指標の組み合わせを用いてもよい。

【0164】

また、変形例３−２に記載のとおり、最適振り易さ指標としては、最適スイングＭＩに限らず様々な指標を設定することができるが、最適振り易さ指標との間に一定の関係（相関）が認められる限り、任意の指標を第１スイング指標とすることができる。例えば、スイング動作時にゴルファーにより発揮される腕エネルギーＥ_AVE、Ｅ₁は、最適振り易さ指標と相関があり、第１スイング指標として用いることができる。また、以下の総肩トルクＴ_ti及び平均肩トルクＴ_AVEも、最適振り易さ指標と相関があり、第１スイング指標とし
て用いることができる。Ｔ_ti及びＴ_AVEは、スイング動作時にゴルファーにより発揮され
る肩周りのトルクを表す指標である。

【0165】

総肩トルクＴ_ti及び平均肩トルクＴ_AVEは、以下のように算出することができる。まず、以下の式に従って、トップからインパクトまでの肩周りのトルクＴ₁及びグリップ４２周りのトルクＴ₂を算出する。

【数26】

【数27】

【0166】

総肩トルクＴ_tiは、上述の肩周りのトルクＴ₁をトップからインパクトまでの区間で積分した値Ｔ_tiである。Ｔ_tiは、トップからインパクトまでの間にゴルファー７の肩周りで発揮されるトルクの総量である。一方、平均肩トルクＴ_AVEは、総肩トルクＴ_tiをトップからインパクトまでの時間で除算したスイング動作中の平均的な肩周りのトルクであり、Ｔ_AVE＝Ｔ_ti／（ｔ_i−ｔ_t）の式に従って算出される。なお、総肩トルクＴ_tiを算出するに当たり、正のトルクＴ₁だけを積分してもよい。

【0167】

また、以下の指標も、第１スイング指標として用いることができる。なお、以下の（７）〜（９）以外の指標は、値が大きい程、最適振り易さ指標が大きくなる傾向があり、（７）〜（９）の指標は、値が小さい程、最適振り易さ指標が大きくなる傾向がある。
（１）トップの時刻でシャフト４０と全体座標系Ｚ軸（グリップよりも下方）とが為す角度θ３（図２６参照）
（２）スイング動作中の角速度ω₂の平均値
（３）トップからインパクトまでの間での角速度ω₂の最大値
（４）トップからインパクトまでのグリップ速度Ｖ_GEの平均値
（５）トップからインパクトまでのグリップ速度Ｖ_GEの最大値
（６）トップからインパクトまでの間でのグリップ４２の移動距離Ｄ
（７）コック解放タイミングｔ_rとインパクトの時間との差分（ここでいうコック解放タイミングｔ_rは、腕とシャフト４０が為すコック角θ４の解放スピードが速まり、腕のエネルギーがシャフト４０のエネルギーへと変わり始めるタイミングと定義することができる。）
（８）コック解放タイミングｔ_rでの腕とシャフト４０が為すコック角θ４（図２６参照）
（９）ダウンスイング時間、すなわち、トップからインパクトまでの時間
（１０）トップの時刻からグリップ４２周りのトルクＴ₂が正負逆転する時刻までのトルクＴ₂の積分値

【0168】

＜３−３−１＞
以下に、図２７〜図２９を参照しつつ、第１スイング指標として腕エネルギーＥ_AVE、平均肩トルクＴ_AVE及びヘッド速度Ｖ_hが算出され、これらの指標Ｅ_AVE、Ｔ_AVE、Ｖ_hに応じて最適クラブ重量の範囲（以下、最適重量帯）が算出される場合の、最適振り易さ決定工程について具体的に説明する。

【0169】

まず、ステップＳ１４０において、決定部２４Ｅは、計測工程で試打されたテストクラブの種類を判定する。プロモデルクラブが試打されていた場合には、ステップＳ１４１に進み、アベレージモデルクラブが試打されていた場合には、ステップＳ１５１に進む。いずれのテストクラブが試打されたかは、入力部２２を介してユーザにより入力される情報に基づいて判定されるものとする。

【0170】

続くステップＳ１４１及びこれに続くステップＳ１４２〜Ｓ１４５は、平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hの大きさに応じて、最適重量帯を決定するステップである。ここでは、Ｔ_AVE，Ｅ_AVE及びＶ_hの値が大きい程、最適重量帯は段階的に大きな値に設定される。

【0171】

具体的には、ステップＳ１４１では、決定部２４Ｅは、ヘッド速度Ｖ_hが４５ｍ／ｓ以上であり、平均肩トルクＴ_AVEが８５Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが３５０Ｎ・ｍ／ｓ以上（以下、プロ条件１）であるか否かを判定する。そして、プロ条件１が満たされる場合には、３２０ｇ以上が最適重量帯であると判定する。一方、ステップＳ１４１でプロ条件１が満たされない場合には、ステップＳ１４２に進む。ステップＳ１４２では、決定部２４Ｅは、平均肩トルクＴ_AVEが７５Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが２５０Ｎ・ｍ／ｓ以上（以下、プロ条件２）であるか否かを判定する。そして、プロ条件２が満たされる場合には、３１０ｇ〜３２０ｇが最適重量帯であると判定する。一方、ステップＳ１４２でプロ条件２が満たされない場合には、ステップＳ１４３に進む。ステップＳ１４３では、決定部２４Ｅは、平均肩トルクＴ_AVEが７０Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１７５Ｎ・ｍ／ｓ以上（以下、プロ条件３）であるか否かを判定する。そして、プロ条件３が満たされる場合には、３０５ｇ〜３１５ｇが最適重量帯であると判定する。一方、ステップＳ１４３でプロ条件３が満たされない場合には、ステップＳ１４４に進む。ステップＳ１４４では、決定部２４Ｅは、平均肩トルクＴ_AVEが５５Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１２０Ｎ・ｍ／ｓ以上（以下、プロ条件４）であるか否かを判定する。そして、プロ条件４が満たされる場合には、３００ｇ〜３１０ｇが最適重量帯であると判定する。一方、ステップＳ１４４でプロ条件４が満たされない場合には、ステップＳ１４５に進む。ステップＳ１４５では、決定部２４Ｅは、ヘッド速度Ｖ_hが３８ｍ／ｓより大きいか否か（以下、プロ条件５）を判定する。そして、プロ条件５が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１４５でプロ条件５が満たされない場合には、プロモデルクラブよりも、アベレージモデルクラブの方が適切であると判定する。これを受けて、表示制御部２４Ｇは、アベレージモデルクラブを用いて計測工程からやり直す旨を示すメッセージを表示部２１上に表示させ、ステップＳ１５１に進む。

【0172】

一方、ステップＳ１５１及びこれに続くステップＳ１５２，Ｓ１５３は、平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hの大きさに応じて、最適重量帯を決定するステップである。ここでは、Ｔ_AVE，Ｅ_AVE及びＶ_hの値が大きい程、最適重量帯は段階的に大きな値に設定される。

【0173】

具体的には、決定部２４Ｅは、ステップＳ１５１において、ヘッド速度Ｖ_hが４０ｍ／ｓ以上であり、平均肩トルクＴ_AVEが６５Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが２００Ｎ・ｍ／ｓ以上（以下、アベレージ条件１）であるか否かを判定する。そして、アベレージ条件１が満たされる場合には、アベレージモデルクラブよりも、プロモデルクラブの方が適切であると判定する。これを受けて、表示制御部２４Ｇは、プロモデルクラブを用いて計測工程からやり直す旨を示すメッセージを表示部２１上に表示させ、ステップＳ１４１に進む。一方、ステップＳ１５１でアベレージ条件１が満たされない場合には、ステップＳ１５２に進む。ステップＳ１５２では、決定部２４Ｅは、平均肩トルクＴ_AVEが５５Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１５０Ｎ・ｍ／ｓ以上（以下、アベレージ条件２）であるか否かを判定する。そして、アベレージ条件２が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１５２でアベレージ条件２が満たされない場合には、ステップＳ１５３に進む。ステップＳ１５３では、決定部２４Ｅは、ヘッド速度Ｖ_hが３８ｍ／ｓより大きく、平均肩トルクＴ_AVEが４８Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１００Ｎ・ｍ／ｓ以上（以下、アベレージ条件３）であるか否かを判定する。そして、アベレージ条件３が満たされる場合には、２９０ｇ〜３００ｇが最適重量帯であると判定する。一方、ステップＳ１５３でアベレージ条件３が満たされない場合には、２９０ｇ以下が最適重量帯であると判定する。

【0174】

以上のステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３は、以下の知見に基づく。すなわち、本発明者らは、プロモデルのゴルフクラブを通常使用しているゴルファー（以下、プロモデルユーザ）１０名及びアベレージモデルのゴルフクラブを通常使用しているゴルファー（以下、アベレージモデルユーザ）１０名に、それぞれプロモデルクラブ及びアベレージモデルクラブを試打させたところ、図２８及び図２９に示す結果を得た。そして、本発明者らは、図２８及び図２９に示す結果から、第１スイング指標を示す空間は、プロモデル領域とアベレージモデル領域とに分割されることを発見した。図２８及び図２９の例では、プロモデル領域とは、平均肩トルクＴ_AVE≧５５Ｎ・ｍかつ腕エネ
ルギーＥ_AVE≧１５０Ｎ・ｍ／ｓの領域であり、アベレージモデル領域とは、平均肩トル
クＴ_AVE＜５５Ｎ・ｍ又は腕エネルギーＥ_AVE＜１５０Ｎ・ｍ／ｓの領域である。しかしながら、これらの数値は、テストクラブの種類等のスイング条件によって異なるものとなり得る。本実験では、プロモデルクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ ＫＥＮＡ Ｂｌｕｅ６ Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１５ｇ、バランスは、Ｄ２）のドライバーが用いられ、アベレージモデルクラブとして、ダンロップスポーツ株式会社製のＸＸＩＯ（登録商標）７（シャフトは、ＭＰ−７００ Ｒ−Ｆｌｅｘ、ゴルフクラブの重量は、２８５ｇ、バランスは、Ｄ１）のドライバーが用いられた。

【0175】

図２８及び図２９を参照して説明した上記実験では、ヘッド速度Ｖ_h及び、飛距離を最大にする最適クラブ重量も算出された。ヘッド速度Ｖ_hは、上記と同様の工程に従って算出された。一方、最適クラブ重量については、ゴルファーに様々な重量のゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブの重量を特定し、これを最適クラブ重量とした。より具体的には、プロモデルユーザには、３２４ｇ，３１６ｇ，３１２ｇ，３０８ｇ，３００ｇの５種類のゴルフクラブをスイングさせ、アベレージモデルユーザには、２９５ｇ，２９０ｇ，２８５ｇのゴルフクラブをスイングさせた。図２８及び図２９に、本実験により得られたヘッド速度Ｖ_h及び最適クラブ重量の値を示す。

【0176】

図２８及び図２９に示すとおり、上記の実験からは、最適クラブ重量は、腕エネルギーＥ_AVE及び平均肩トルクＴ_AVEが共に大きいほど、大きくなることが分かる。その結果、本発明者らは、平均肩トルクＴ_AVE−腕エネルギーＥ_AVE−ヘッド速度Ｖ_h空間を、プロモデルユーザについては図２８に示すように、アベレージモデルユーザについては図２９に示すように領域分割することで、最適重量帯を示す領域を定義可能であることを発見した。ただし、簡単のため、図２８及び図２９では、ヘッド速度Ｖ_hを表す軸は省略され、平均肩トルクＴ_AVE−腕エネルギーＥ_AVE平面が示されている。すなわち、上述のステップＳ１４１〜Ｓ１４５及びステップＳ５１〜Ｓ５３は、平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hを示す点が、Ｔ_AVE−Ｅ_AVE−Ｖ_h空間におけるどの領域にプロットされるかに応じて、最適重量帯を判定するステップとなっている。なお、ステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３での判定に使用されるＴ_AVE，Ｅ_AVE及びＶ_hの閾値、言い換えると、図２８及び図２９に示す分割領域の境界を示す値は、テストクラブのモデル毎に整理され、対応関係データ２８として記憶部２３内に格納されている。すなわち、本変形例に係る対応関係データ２８とは、Ｔ_AVE，Ｅ_AVE及びＶ_hの大きさと最適重量帯との対応関係を定めるデータである。ステップＳ１４１〜Ｓ１４５及びステップＳ１５１〜Ｓ１５３では、この対応関係データ２８が参照され、上記の判定が行われる。

【0177】

＜３−３−２＞
以下に、図３０及び図３１を参照しつつ、第１スイング指標として腕エネルギーＥ_AVE、平均肩トルクＴ_AVE及びヘッド速度Ｖ_hが算出され、これらの指標Ｅ_AVE、Ｔ_AVE、Ｖ_hに応じて最適スイングＭＩ帯が算出される場合の、最適振り易さ決定工程について具体的に説明する。

【0178】

本変形例に先立つ計測工程では、プロモデルクラブ及びアベレージモデルクラブの２本のテストクラブに代えて、一本のセンサユニット１付きテストクラブがゴルファー７により試打される。ただし、本変形例でも、２本のテストクラブを打ち分けるようにし、フィッティングの精度を向上させることも可能である。

【0179】

具体的には、ステップＳ１６０では、決定部２４Ｅは、平均肩トルクＴ_AVEが７７Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１８０Ｎ・ｍ／ｓ以上（以下、条件１）であるか否かを判定する。条件１が満たされる場合には、決定部２４Ｅは、ヘッド速度Ｖ_hが４５ｍ／ｓ以上であるか否か（以下、条件２）を判定する（ステップＳ１６１）。そして、条件２が満たされる場合には、５６００ｋｇ・ｃｍ²以上が最適スイングＭＩ帯であると判定し、そうでなければ、５５９０〜５６３０ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ１６０で条件１が満たされない場合には、ステップＳ１６２に進む。ステップＳ１６２では、決定部２４Ｅは、平均肩トルクＴ_AVEが６０Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１７０Ｎ・ｍ／ｓ以上（以下、条件３）であるか否かを判定する。条件３が満たされる場合には、決定部２４Ｅは、ヘッド速度Ｖ_hが４５ｍ／ｓ以上であるか否か（以下、条件４）を判定する（ステップＳ１６３）。そして、条件４が満たされる場合には、５５９０〜５６３０ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定し、そうでなければ、５５１０〜５５９０ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定する。一方、ステップＳ１６２で条件３が満たされない場合には、ステップＳ１６４に進む。ステップＳ１６４では、決定部２４Ｅは、平均肩トルクＴ_AVEが５０Ｎ・ｍ以上であり、かつ、腕エネルギーＥ_AVEが１３０Ｎ・ｍ／ｓ以上（以下、条件５）であるか否かを判定する。そして、条件５が満たされる場合には、５４６０〜５５１０ｋｇ・ｃｍ²が最適スイングＭＩ帯であると判定し、そうでなければ、５４８０ｋｇ・ｃｍ²以下が最適スイングＭＩ帯であると判定する。

【0180】

以上のステップＳ１６０〜Ｓ１６４は、以下の知見に基づく。すなわち、本発明者らは、本実施形態に係るテストクラブをゴルファー２１名に試打させてスイング指標を算出したところ、図３１に示す結果を得た。なお、この実験で算出されたスイング指標も、平均肩トルクＴ_AVE及び腕エネルギーＥ_AVEであり、具体的な値は、上記と同様の工程に従って算出された。本実験では、テストクラブとして、ダンロップスポーツ株式会社製のＳＲＩＸＯＮ（登録商標）Ｚ−５２５（シャフトは、Ｍｉｙａｚａｋｉ Ｋｏｓｕｍａ Ｂｌｕｅ６ Ｓ−Ｆｌｅｘ、ゴルフクラブの重量は、３１４ｇ、バランスは、Ｄ３）のドライバーが用いられた。また、本実験では、ヘッド速度Ｖ_h及び、飛距離を最大にする最適スイングＭＩも算出された。ヘッド速度Ｖ_hは、上記と同様の工程に従って算出された。一方、最適スイングＭＩについては、ゴルファーに様々なスイング慣性モーメントのゴルフクラブをスイングさせ、最大の飛距離を与えるゴルフクラブのスイング慣性モーメントを特定し、これを最適スイングＭＩとした。より具体的には、スイング慣性モーメントが５６５０ｋｋｇ・ｃｍ²、５６１０ｋｇ・ｃｍ²、５５５０ｋｇ・ｃｍ²、５４８５ｋｇ・ｃｍ²、５４００ｋｇ・ｃｍ²の５種類のゴルフクラブをスイングさせた。図３１には、本実験により得られたヘッド速度Ｖ_h及び最適スイングＭＩの値が示されている。

【0181】

図３１に示すとおり、上記の実験からは、最適スイングＭＩは、腕エネルギーＥ_AVE及び平均肩トルクＴ_AVEが共に大きいほど、大きくなることが分かる。その結果、本発明者らは、平均肩トルクＴ_AVE−腕エネルギーＥ_AVE−ヘッド速度Ｖ_h空間を、図３１に示すように領域分割することで、最適スイングＭＩ帯を示す領域を定義可能であることを発見した。ただし、簡単のため、図３１では、ヘッド速度Ｖ_hを表す軸は省略され、平均肩トルクＴ_AVE−腕エネルギーＥ_AVE平面が示されている。すなわち、上述のステップＳ１６０〜Ｓ１６４は、平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hを示す点が、Ｔ_AVE−Ｅ_AVE−Ｖ_h空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定するステップとなっている。なお、ステップＳ１６０〜Ｓ１６４での判定に使用されるＴ_AVE，Ｅ_AVE及びＶ_hの閾値、言い換えると、図３１に示す分割領域の境界を示す値は、対応関係データ２８として記憶部２３内に格納されている。すなわち、本変形例に係る対応関係データ２８とは、Ｔ_AVE，Ｅ_AVE及びＶ_hの大きさと最適スイングＭＩ帯との対応関係を定めるデータである。ステップＳ１６０〜Ｓ１６４では、この記憶部２３内の対応関係データ２８が参照され、上記の判定が行われる。

【0182】

また、上記のとおり、グリップエンド周りの慣性モーメントＩ_Gは、スイング慣性モーメントＩ_Sと概ね比例関係にある。従って、図３０の最適振り易さ決定工程で決定される最適スイングＭＩ帯を、ゴルファーに適した慣性モーメントＩ_Gの範囲に変換することで、Ｉ_Gに基づく最適振り易さ決定工程を実行することができる。図３２は、図３０のフローチャートを改変したものであり、最適グリップエンドＭＩを決定するためのフローチャートである。なお、図３１の右側の凡例には、最適スイングＭＩを最適グリップエンドＭＩに換算した値も示されている。

【0183】

＜３−４＞
第１及び第２実施形態では、シャフトの剛性として、曲げ剛性が評価されたが、これに代えて、ねじれ剛性を評価してもよい。ねじれ剛性の値（以下、ＧＪ値）も、シャフト４０の延びる方向に沿った複数の位置において測定又は算出することができる。すなわち、シャフト４０の延びる方向に沿った複数の位置におけるねじれ剛性の分布を、シャフトの剛性としてもよい。この場合、最適剛性指標としては、ゴルファー７に適したＧＪ値（最適ＧＪ値）が決定されることになるが、最適ＧＪ値を決定するための第２スイング指標としては、最適ＧＪ値との相関が認められる任意の指標を用いることができる。このような第２スイング指標としては、例えば、特開２０１４−２１２８６２号公報に記載されているような、以下の指標を用いることができる。

【0184】

（１）グリップ角速度ω_yが最大となるときからインパクトまでの単位時間あたりのグリップ角速度ω_xの変化量の大きさ
（２）トップ付近でのグリップ角速度ω_zの変化量
（３）トップからダウンスイング途中であってグリップ角速度ω_yが最大となるときまでのグリップ角速度ω_zの変化量の大きさ

【0185】

本変形例でも、第２スイング指標と最適ＧＪ値との関係を表す近似式を予め実験により算出し、記憶部２３内に格納しておくことで、計測工程で得られる計測データに基づく第２スイング指標から、最適ＧＪ値を決定することができる。

【0186】

また、第１実施形態において、最適剛性指標として、ゴルファー７に適したシャフト４０の複数の位置における剛性分布ではなく、ゴルファー７に適したシャフト４０のフレックス、調子又はトルクを決定するようにしてもよい。同様に、第２実施形態においても、最適剛性指標として、ゴルファー７に適したシャフト４０のフレックスではなく、ゴルファー７に適したシャフト４０の調子、トルク、又は複数の位置における剛性分布を決定するようにしてもよい。なお、トルクは、シャフト４０全体でのねじれ剛性を表す指標である。

【0187】

＜３−５＞
上記実施形態では、上述した（１）〜（５）を前提とした肩挙動導出工程の一例について説明したが、肩挙動導出工程のアルゴリズムはこれに限定されず、例えば、以下の態様とすることができる。

【0188】

本変形例に係る肩挙動導出工程では、上述の（１）〜（５）の前提のうち、（４）及び（５）が省略され、（１）〜（３）のみが前提とされる。具体的には、肩挙動導出部２４Ｃは、第２変換工程で得られるスイング平面Ｐ内でのグリップ４２の軌道を、円弧（円）に近似する（図３０参照）。そして、当該円弧（円）の中心を肩の位置Ｐ_s＝（Ｐ_sX，Ｐ_sY）とし、当該円弧（円）の中心からグリップ４２の軌道までの平均的な距離を、腕長さ（肩とグリップ４２との距離）Ｒとする。この場合、肩挙動導出処理に先行する第２変換工程においては、グリップ挙動導出部２４Ｂにより、スイング平面Ｐ内でのグリップ４２の軌道が算出されているものとする。例えば、グリップ４２の軌道は、スイング平面Ｐ内でのグリップ速度（ｖ_pY，ｖ_pZ）を積分し、グリップ４２の軌道上の点Ａ_i＝（Ｘ_i，Ｙ_i）（ｉ＝１，２，・・・）を算出することにより導出可能である。なお、ここでいうｉは、時系列に沿ったデータ番号である。

【0189】

以下に、グリップ４２の軌道の近似円（円弧）を導出する方法の一例を示す。まず、グリップ４２の軌道上の任意の３点、例えば、Ａ_i，Ａ_i+30，Ａ_i+60を考える。このとき、任意の三角形の外接円の中心は、当該三角形の三辺の垂直二等分線の交点となることから、Ａ_i，Ａ_i+30，Ａ_i+60を頂点とする三角形を考えたとき（図３１参照）、以下の数２４の式、ひいては数２５の式が成り立つ。

【数28】

【数29】

【0190】

そして、様々なｉに対する複数の数２５の式から以下の数２６の式を作成し、疑似逆行列を導出する。これにより、グリップ４２の軌道の近似円（円弧）の中心Ｐ_s＝（Ｐ_sX，Ｐ_sY）を導出することができる。

【数30】

【0191】

続いて、肩挙動導出部２４Ｃは、近似円（円弧）の中心Ｐ_s＝（Ｐ_sX，Ｐ_sY）からグリップ４２の軌道上の各点Ａ_i＝（Ｘ_i，Ｙ_i）までの距離の平均値を算出し、腕長さＲとする。そして、肩挙動導出部２４Ｃは、この腕長さＲに基づいて、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω₁＝Ｖ_GE／Ｒを算出する。

【0192】

図３２Ａは、上記実施形態に係る方法で再現したトップからインパクトまでの二重振り子モデルを示しており、図３２Ｂは、本変形例に係る方法で再現したトップからインパクトまでの二重振り子モデルを示している。これらの図を比較すると明らかであるが、本変形例では、前提（４）が省略されることにより、トップでの腕とゴルフクラブ４との為す角度が評価されるため、コック角の浅い／深いを評価することができる。また、前提（５）が省略されることにより、インパクト時の腕の角度が評価されるため、ハンドファスト／ハンドレイトを評価することができる。

【0193】

＜３−６＞
第２実施形態では、第１スイング指標に基づいて最適振り易さ指標（最適スイングＭＩ）を決定した後、この最適スイングＭＩに基づいて最適特性指標（最適シャフト重量）が決定された。しかしながら、第１スイング指標の大きさに応じて、最適特性指標を直接算出することもできる。具体的には、例えば、第１スイング指標を説明変数とし、最適特性指標を目的変数とする重回帰式を、対応関係データ２２８として予め算出し記憶しておくことができる。この場合、第１指標算出工程で導出される第１スイング指標を当該重回帰式に代入することにより、最適特性指標を導出することができる。あるいは、第１実施形態の図１４及び図１５と同様に、（好ましくは、ヘッド４１の種類毎に）Ｐ_{1_AVE}−Ｐ_{2_AVE}（−Ｖ_h）空間を領域分割して、最適特性指標に対応する領域を定義しておき（図３６参照）、当該定義情報を対応関係データ２２８として予め記憶しておくことができる。この場合、第１指標算出工程で導出される第１スイング指標がどの領域に属するかに応じて、最適特性指標を導出することができる。

【0194】

また、図３６の対応関係データ２２８に代えて、図３７Ａ〜図３７Ｃに示すような対応関係データ２２８を用いることができる。図３７Ａは、最適フレックスが「Ｘ」の場合に最適シャフト重量帯を決定するためのデータであり、図３７Ｂは、最適フレックスが「Ｓ」の場合に最適シャフト重量帯を決定するためのデータであり、図３７Ｃは、最適フレックスが「ＳＲ」の場合に最適シャフト重量帯を決定するためのデータである。すなわち、まず最適剛性指標を決定し、これに対応する対応関係データ２２８を選択する。そして、選択された対応関係データ２２８に基づき、第１スイング指標の大きさに応じて、最適特性指標を決定する。この方法によれば、最適振り易さ指標を決定することなく、第１スイング指標から最適特性指標を直接算出することができる他、最適特性指標よりも最適剛性指標を優先して、推奨シャフトを決定することができる。

【0195】

＜３−７＞
第２実施形態では、最適特性指標として最適シャフト重量が導出されたが、これに代えて又は加えて、ゴルファー７に適したシャフト４０以外の部位の特性を表す指標、例えば、ゴルファーに適したグリップ４２の重量を導出することもできる。また、ゴルフクラブ４の特定の部位の特性として、重量以外の特性、例えば、慣性モーメントを導出することもできる。

【0196】

＜３−８＞
上述した最適振り易さ指標及び最適特性指標は、ゴルファー７に適したゴルフクラブ４の全体又は特定の部位の特性を表す指標（以下、最適クラブ指標という）である。第１実施形態及び変形例３−６では、第１スイング指標が所定の領域Ａ１〜Ａ５，Ｂ１〜Ｂ３，Ｃ１〜Ｃ５のいずれに属するかに応じて最適クラブ指標が決定されたが、より精度の高いフィッティングのため、最適クラブ指標の決定時に以下の境界処理を実行することができる。本変形例に係る境界処理は、計測データに基づいて算出された第１スイング指標が、第１スイング指標を表す空間を分割することにより定義される複数の領域の境界付近に存在する場合には、既出の最適クラブ指標を修正する処理である。

【0197】

ここでは、一例として、変形例３−６と組み合わせる場合の境界処理について説明する。具体的な処理の流れは、図３８に示すとおりである。図３８に示す境界処理は、最適クラブ指標である最適シャフト重量が、図３６に示す判定基準（分割領域Ｃ１〜Ｃ５の境界線Ｌ６〜Ｌ１０）に従って決定された後に実行される。

【0198】

まず、ステップＳ６５において、決定部２２４Ｅは、判定指標が判定領域の上限境界付近に存在するか否かを判定する。なお、判定指標とは、計測データに基づいて算出されたＰ_{1_AVE}，Ｐ_{2_AVE}（第１スイング指標）を意味し、判定領域とは、Ｐ_{1_AVE}−Ｐ_{2_AVE}（−Ｖ_h）空間内に規定される領域Ｃ１〜Ｃ５のうち、判定指標が属する領域である。また、判定指標が上限境界付近に存在するか否かは、例えば、Ｐ_{1_AVE}−Ｐ_{2_AVE}（−Ｖ_h）空間内において、判定指標を表す点と判定領域の上限境界線との距離が閾値以内であるかによって判定することができる。判定指標が領域Ｃ２に属する場合には、境界線Ｌ７が上限境界線である。

【0199】

ステップＳ６５において判定指標が判定領域の上限境界付近に存在すると判定された場合には、決定部２２４Ｅは、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ６６）。

【0200】

一方、ステップＳ６５において判定指標が判定領域の上限境界付近に存在しないと判定された場合には、処理はステップＳ６７に進む。ステップＳ６７では、判定指標が判定領域の下限境界付近に存在するか否かを判定する。なお、判定指標が下限境界付近に存在するか否かは、例えば、Ｐ_{1_AVE}−Ｐ_{2_AVE}（−Ｖ_h）空間内において、判定指標を表す点と判定領域の下限境界線との距離が閾値以内であるかによって判定することができる。判定指標が領域Ｃ２に属する場合には境界線Ｌ８が下限境界線である。

【0201】

ステップＳ６７において判定指標が判定領域の下限境界付近に存在すると判定された場合には、決定部２２４Ｅは、既出の最適シャフト重量を所定量だけ減量するように調整する（ステップＳ６８）。そして、以上のステップが終了すると、境界処理は終了する。

【0202】

＜３−９＞
変形例３−８では、最適クラブ指標をより精度よく決定するための境界処理について説明したが、本変形例に係る例外処理も、最適クラブ指標をより精度よく決定するための処理である。上記実施形態では、第１スイング指標を算出するに当たり、スイング平面Ｐ内での解析が行われたが、上述した第１スイング指標であるＰ_{1_AVE}，Ｐ_{2_AVE}や変形例３−３で示した指標は、スイング平面Ｐ内に投影された二次元的なゴルファー７のスイング動作を表す指標である。しかしながら、実際のスイング動作は、三次元的に行われる。本変形例に係る例外処理は、二次元的な解析による誤差を軽減し、最適クラブ指標の精度を向上させるための処理である。

【0203】

本変形例に係る例外処理では、第１スイング指標として、上述したものに加えて、スイング平面Ｐに現れないローテーションの動きを示す指標、及びコックの動きを表す指標（以下、コック指標という）が導出される。スイング平面Ｐに現れないローテーションの動きとは、スイング平面Ｐに投影されたゴルフクラブ４のシャフト軸周りの捻りのローテーションの動き（図３９Ａ参照）を表す指標（以下、捻り指標という）や、スイング平面Ｐから飛び出す方向（プッシュ方向）のゴルフクラブ４のローテーションの動き（図３９Ｂ参照）を表す指標（以下、プッシュ指標という）により評価することができる。ゴルファー７によっては、このようなローテーションの動きによりヘッド速度を獲得している。また、ゴルファー７によっては、ダウンスイング後半のコックの動きによりヘッド速度を獲得することもあるが、そのようなコックの動きの全体像は、必ずしもスイング平面Ｐ上に投影されない。ここで説明する例外処理では、このようなゴルファー７についても精度よく、最適クラブ指標を導出することができる。

【0204】

以下では、一例として、変形例３−６と組み合わせる場合の例外処理について説明する。具体的な処理の流れは、図４０に示すとおりである。図４０に示す例外処理は、最適クラブ指標である最適シャフト重量が、図３６に示す判定基準（分割領域Ｃ１〜Ｃ５の境界線Ｌ６〜Ｌ１０）に従って決定された後に実行される。

【0205】

まず、ステップ７１では、算出部２４Ｄが、計測データに基づいて、捻り指標を算出する。捻り指標は、例えば、シャフト軸周りの角速度ω_zのトップからインパクトまでの平均値として算出することができる。続いて、決定部２２４Ｅが、ヘッド速度Ｖ_hが所定範囲内であり、かつ、捩り指標が所定値以上であるかを判定する。そして、かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７２）。一方、かかる条件が満たされない場合には、処理はステップＳ７３に進む。

【0206】

ステップＳ７３では、算出部２４Ｄが、計測データに基づいて、プッシュ指標を算出する。続いて、決定部２２４Ｅが、ヘッド速度Ｖ_hが所定範囲内であり、プッシュ指標が所定値以上であり、かつ、既出の最適シャフト重量が所定範囲内であるか否かを判定する。かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７４）。一方、かかる条件が満たされない場合には、処理はステップＳ７５に進む。なお、プッシュ指標が所定値以上の場合とは、例えば、トップのω_x≧０、かつ、トップからインパクトまでのω_yの平均値／トップからインパクトまでのω_xの平均値≧１．５が満たされる場合とすることができる。

【0207】

ステップＳ７５では、算出部２４Ｄが、計測データに基づいて、コック指標を算出する。コック指標は、例えば、トップ付近のω_yの差と、ω_yが最大となる時刻からインパクトまでのω_yの平均値とを比較した値として算出することができる。続いて、決定部２２４Ｅが、ヘッド速度Ｖ_hが所定範囲内であり、コック指標が所定値以上であり、かつ、既出の最適シャフト重量が所定範囲内であるか否かを判定する。かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７６）。

【0208】

以上のステップが完了すると、例外処理は終了する。本例外処理と、変形例３−８の境界処理とは、いずれも最適クラブ指標をより精度よく決定するための処理であり、両処理は、組み合わせることもできる。この場合、境界処理の後、本例外処理を実行することが好ましい。

【実施例】

【0209】

以下、本発明の実施例について説明する。但し、本発明は、以下の実施例に限定されない。

【0210】

＜実施例１＞
スイング慣性モーメントＩ_Sが５６１５ｋｇ・ｃｍ²のゴルフクラブを普段使用しているゴルファーに対し、特許文献１のフィッティング方法（比較例）を適用するとともに、第１実施形態に係るフィッティング方法（実施例１）を適用した。

【0211】

比較例の方法では、マイクラブ重量±５ｇの範囲に収まるようなゴルフクラブしか選択されないため、Miyazaki Indio 7S(Ｉ_S＝約５６１０ｋｇ・ｃｍ²)のシャフトが選択された。そして、当該シャフトを備えるゴルフクラブを用いて、当該ゴルファーに５球の試打を行わせ、打球の飛距離及び左右ぶれの平均値を計測した。

【0212】

一方、実施例１の方法では、Miyazaki Indio 6S(Ｉ_S＝約５５４０ｋｇ・ｃｍ²)のシャフトが選択された。そして、当該シャフトを備えるゴルフクラブであって、比較例の場合とはシャフトのみが異なる（すなわち、ヘッドは同一の）ゴルフクラブを用いて、当該ゴルファーに５球の試打を行わせ、打球の飛距離及び左右ぶれの平均値を計測した。

【0213】

図４１Ａ及び図４１Ｂは、実施例１及び比較例による飛距離及び左右ぶれの平均値及びばらつきを示すグラフである。同図に示すように、飛距離及び左右ぶれの両方について、比較例の方法で選定されたゴルフクラブよりも実施例１の方法で選定されたゴルフクラブを用いた方が、好ましい結果が得られた。従って、実施例１の優位性が確認された。

【0214】

＜実施例２＞
第２実施形態に係るフィッティング方法（実施例２）を適用して、あるゴルファーに適切なゴルフクラブを選択したところ、ヘッドＡ（重量：２００ｇ）にシャフトＡ（フレックス：ＳＲ、重量：６５．５ｇ）を組み合わせたゴルフクラブＡが選択された。当該ゴルファーは、ヘッドＡと同等の重量のゴルフクラブを普段使用していたため、ヘッドＡが嗜好クラブされた。また、ゴルフクラブＡと重量が同等であり、振り易さ指標が同等であるゴルフクラブとして、ヘッドＢ（重量：２０２ｇ）にシャフトＢ（フレックス：ＳＲ、重量：５６．５ｇ）を組み合わせたゴルフクラブＢを用意した。そして、当該ゴルファーにゴルフクラブＡ及びＢの両方を５回ずつ試打させたところ、図４２Ａ及び図４２Ｂに示す結果が得られた。図４２Ａは、ゴルフクラブＡの試打結果であり、図４２Ｂは、ゴルフクラブＢの試打結果である。

【0215】

図４２Ａ及び図４２Ｂを比較すれば明らかなとおり、普段使用しているヘッドと同等の重量のヘッドＡを含むゴルフクラブＡをスイングした方が、左右ずれ及び飛距離の両方において優れていることが分かる。すなわち、ヘッドの重量を調整して最適振り易さ指標を実現するよりも、ヘッドをゴルファーの嗜好に合わせた上でシャフトの重量を調整して最適振り易さ指標を実現する方が適切であることが分かる。従って、ヘッドを固定した上で、最適振り易さ指標を達成できるようなシャフトを選択する上記実施形態に係るフィッティング方法の優位性が確認された。

【符号の説明】

【0216】

１ センサユニット（計測装置）
２ フィッティング装置
３ フィッティングプログラム
４ ゴルフクラブ
７ ゴルファー
２４Ａ 取得部
２４Ｂ グリップ挙動導出部
２４Ｃ 肩挙動導出部
２４Ｄ 算出部
２４Ｅ 決定部
２４Ｆ 選択部
４０ シャフト
４１ ヘッド
４２ グリップ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】

【図25】

【図26】

【図27】

【図28】

【図29】

【図30】

【図31】

【図32】

【図33】

【図34】

【図35A】

【図35B】

【図36】

【図37A】

【図37B】

【図37C】

【図38】

【図39A】

【図39B】

【図40】

【図41A】

【図41B】

【図42A】

【図42B】