JP2022181564A - ゴルフクラブのフィッティング装置、方法及びプログラム - Google Patents

Abstract

【課題】ゴルファーに適した長さのゴルフクラブを精度よく選定するためのフィッティング装置を提供する。【解決手段】フィッティング装置は、取得部と、算出部と、決定部とを備える。取得部は、ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得する。算出部は、前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標を算出する。決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの長さである最適クラブ長さを決定する。【選択図】図６

Description

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

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

特開２０１７－１７０１０５号公報

特許文献１の方法によれば、最適振り易さ指標及び最適剛性指標にできる限り合致するシャフト、ひいてはゴルフクラブを絞り込むことができるが、ゴルフクラブとしてのクラブ長さを調整するため、異なるシャフト長さのオプションが用意されていることがある。このため、特許文献１に開示する方法に限らず、従来の方法により絞り込まれたゴルフクラブの中から、さらにゴルファーに適した長さのゴルフクラブを選択することができる方法が望まれていた。なお、このことは、何らかの方法によってゴルファーに適するシャフト及びゴルフクラブが絞り込まれていない場合に、ゴルファーに適したシャフトまたはゴルフクラブを選択しようとする場合についても同様である。

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

第１観点に係るフィッティング装置は、取得部と、算出部と、決定部とを備える。取得部は、ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得する。算出部は、前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標を算出する。決定部は、前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの長さである最適クラブ長さを決定する。

第２観点に係るフィッティング装置は、第１観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標を所定の第１閾値と比較し、前記第１スイング指標が前記第１閾値以上である場合には、前記最適クラブ長さが第１クラブ長さであると決定する。

第３観点に係るフィッティング装置は、第２観点に係るフィッティング装置であって、前記決定部は、前記第１スイング指標を前記第１閾値よりも小さい所定の第２閾値と比較し、前記第１スイング指標が前記第２閾値以下である場合には、前記最適クラブ長さが前記第１クラブ長さよりも長い第２クラブ長さであると決定する。

第４観点に係るフィッティング装置は、第１観点から第３観点のいずれかに係るフィッティング装置であって、前記計測データには、前記ゴルフクラブのグリップエンドにおける角速度が含まれ、前記第１スイング指標は、前記ゴルフクラブのシャフトに平行な軸周りの角速度に基づいて算出される。

第５観点に係るフィッティング装置は、第４観点に係るフィッティング装置であって、前記第１スイング指標は、ダウンスイング時における前記角速度の平均値または積分値である。

第６観点に係るフィッティング装置は、第１観点から第５観点のいずれかに係るフィッティング装置であって、前記最適クラブ長さに最も合致するゴルフクラブ及び当該ゴルフクラブに含まれるべきシャフトの少なくとも一方を選択する選択部をさらに備える。

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

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

第９観点に係るフィッティング装置は、第７観点または第８観点に係るフィッティング装置であって、前記計測データには、前記ゴルフクラブのグリップエンドにおける角速度が含まれ、前記第１スイング指標は、前記ゴルフクラブのシャフトに平行な軸周りの角速度に基づいて算出される。

第１０観点に係るフィッティング方法は、以下の（１）～（３）のステップを含む。
（１）ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測デー タを取得するステップ
（２）前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標をコンピュータを用いて算出するステップ
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定するステップ

第１１観点に係るフィッティング方法は、以下の（１）～（３）のステップをコンピュータに実行させる。
（１）ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得するステップ
（２）前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標をコンピュータを用いて算出するステップ
（３）前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定するステップ

本発明の第１観点によれば、ゴルフクラブによる計測値に基づいて、ゴルフクラブの捩り度合いを表す第１スイング指標が算出される。そして、第１スイング指標に応じて、ゴルファーに適したゴルフクラブの長さである最適クラブ長さが決定される。すなわち、ゴルファーのスイング動作の特徴に基づき、これに適した最適クラブ長さが決定されるため、ゴルファーに適したクラブ長さのゴルフクラブを精度よく選定することができる。

本発明の第７観点によれば、ゴルフクラブによる計測値に基づいて、第１スイング指標、第２スイング指標及び第３スイング指標が算出される。そして、第１スイング指標に応じてゴルファーに適したゴルフクラブの長さである最適クラブ長さが決定され、第２スイング指標に応じてゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標が決定され、第３スイング指標に応じてゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適振り易さ指標、最適剛性指標及び最適クラブ長さに最も合致するゴルフクラブ及びこれに含まれるべきシャフトのうち少なくとも一方が選定される。第１スイング指標は、ゴルフクラブの捩り度合いを表す。これにより、様々な観点からゴルファーに適したゴルフクラブが絞り込まれるので、ゴルフクラブを精度よく選定することができる。

本発明の第８観点によれば、ゴルフクラブによる計測値に基づいて、第１スイング指標、第２スイング指標及び第３スイング指標が算出される。そして、第１スイング指標に応じてゴルファーに適したゴルフクラブの長さである最適クラブ長さが決定され、第２スイング指標に応じてゴルファーに適したゴルフクラブ特定の部位の特性を表す最適特性指標が決定され、第３スイング指標に応じてゴルファーに適したシャフトの剛性を示す最適剛性指標が決定される。そして、最適特性指標、最適剛性指標及び最適クラブ長さに最も合致するゴルフクラブ及びこれに含まれるべきシャフトのうち少なくとも一方が選定される。第１スイング指標は、ゴルフクラブの捩り度合いを表す。これにより、様々な観点からゴルファーに適したゴルフクラブが絞り込まれるので、ゴルフクラブを精度よく選定することができる。

本発明の一実施形態に係るフィッティング装置を備えるフィッティングシステムを示す図。 フィッティングシステムの機能ブロック図。 ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。 フィッティング処理の流れを示すフローチャート。 （Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。 第１指標算出工程の流れを示すフローチャート。 トップの時刻の導出方法を説明する図。 最適クラブ長さ決定工程の流れを示すフローチャート。 スイング挙動の違いによる角速度ω_zの波形の違いを示すグラフ。 スイング平面を説明する図。 二重振り子モデルを概念的に説明する図。 二重振り子モデルを概念的に説明する別の図。 プロモデル領域を示す図。 アベレージモデル領域を示す図。 インターナショナル・フレックス・コード（ＩＦＣ）を説明する図。 シャフトの曲げ剛性の測定方法を説明する図。 スイング中のシャフトの曲げを説明する図。 変形例に係る第１スイング指標となる角度を説明する図。 最適シャフト重量帯に対応する分割領域に分割された第２スイング指標を示す空間を示す図。 特定のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。 別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。 さらに別のフレックスに対する、最適シャフト重量帯に対応する分割領域に分割された第１スイング指標を示す空間を示す図。 捩りのローテーションの動きを示す概念図。 プッシュのローテーションの動きを示す概念図。 例外処理の流れを示す図。 ゴルファーによる実施例及び比較例による試打結果を比較する図。 別のゴルファーによる実施例及び比較例による試打結果を比較する図。 さらに別のゴルファーによる実施例及び比較例による試打結果を比較する図。 さらに別のゴルファーによる実施例及び比較例による試打結果を比較する図。

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

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

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

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

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

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

＜１－２．フィッティング装置の構成＞
フィッティング装置２は、ハードウェアとしては汎用のコンピュータであり、例えば、デスクトップ型コンピュータ、ノート型コンピュータ、タブレットコンピュータ、スマートフォン等として実現される。図２に示す通り、フィッティング装置２は、本実施形態に係るフィッティングプログラム３（以下、単に「プログラム３」とも称する）を、汎用のパーソナルコンピュータにインストールすることにより製造される。プログラム３は、コンピュータで読み取り可能なＣＤ－ＲＯＭ等の記録媒体２０から、或いは通信部２５に接続されるローカルエリアネットワーク（ＬＡＮ）やインターネット等の通信ネットワークを介して、フィッティング装置２に取得される。プログラム３は、センサユニット１から送られてくる計測データに基づいてスイング動作を解析し、ゴルファーＧに適したゴルフクラブ４を選択するのを支援する情報を出力するためのソフトウェアである。プログラム３は、フィッティング装置２に後述する動作を実行させる。

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

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

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

＜２．フィッティング処理＞
続いて、システム１００により実行されるフィッティング処理について説明する。本実施形態に係るフィッティング処理は、図４に示すとおり、以下の１１の工程（Ｓ１～Ｓ１１）から構成されている。
（Ｓ１）ｘｙｚ局所座標系でのグリップ加速度ａ_x，ａ_y，ａ_z、グリップ角速度ω_x，ω_y，ω_z及びグリップ地磁気ｍ_x，ｍ_y，ｍ_zの計測データを計測する計測工程
（Ｓ２）グリップ角速度ω_zから第１スイング指標を算出する第１指標算出工程
（Ｓ３）第１スイング指標に基づいて、最適クラブ長さを決定する最適クラブ長さ決定工程
（Ｓ４）計測工程で得られたｘｙｚ局所座標系での計測データを、ＸＹＺ全体座標系でのグリップ加速度ａ_X，ａ_Y，ａ_Z及びグリップ角速度ω_X，ω_Y，ω_Zに変換する第１変換工程（第１変換工程では、ＸＹＺ全体座標系でのグリップ速度ｖ_X，ｖ_Y，ｖ_Zも導出される。）
（Ｓ５）ＸＹＺ全体座標系でのグリップ４２の挙動（グリップ角速度ω_X，ω_Y，ω_Z及びグリップ速度ｖ_X，ｖ_Y，ｖ_Z）を、スイング平面Ｐ（後述する）内でのグリップ４２の挙動へと変換する第２変換工程
（Ｓ６）スイング平面Ｐ内でのグリップ４２の挙動に基づいて、スイング平面Ｐ内でのゴルファーＧの疑似的な肩の挙動を導出する肩挙動導出工程
（Ｓ７）スイング平面Ｐ内でのグリップ４２の挙動及び疑似的な肩の挙動に基づいて、第２スイング指標（本実施形態では、後述する腕出力パワーＰ₁__AVE、クラブ入力パワーＰ₂__AVE及びヘッド速度Ｖ_h）を算出する第２指標算出工程
（Ｓ８）第２スイング指標に基づいて、ゴルファーＧに適した振り易さ指標（本実施形態では、後述するスイング慣性モーメントＩ_S）である最適振り易さ指標（本実施形態では、最適スイングＭＩ）を決定する最適振り易さ決定工程
（Ｓ９）計測データに基づいて、第３スイング指標（本実施形態では、後述する第１～第４特徴量Ｆ₁～Ｆ₄）を算出する第３指標算出工程
（Ｓ１０）第３スイング指標に基づいて、ゴルファーＧに適したシャフト４０の剛性を示す最適剛性指標（本実施形態では、後述するＥＩ分布）を決定する最適剛性決定工程
（Ｓ１１）最適振り易さ指標、最適剛性指標及び最適クラブ長さに合致するゴルフクラブ４、特にシャフト４０を選択する最適クラブ選択工程
以下、これらの工程を順に説明する。

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

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

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

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

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

＜２－２．第１指標算出工程＞
以下、図６を参照しつつ、第１指標算出工程（Ｓ２）について説明する。第１指標算出工程では、グリップ角速度ω_x，ω_y，ω_zの計測データに基づいて、第１スイング指標ＳＷ１が算出される。第１スイング指標ＳＷ１は、スイング動作中のゴルフクラブ４の捩り度合い、言い換えると、ゴルフクラブ４のシャフト軸周りの捻りのローテーションの動き（図２１Ａ参照）を表す指標である。本実施形態では、第１スイング指標ＳＷ１として、トップからインパクトまでのダウンスイング時のグリップ角速度ω_zの平均値が算出される。第１スイング指標ＳＷ１は、後述する工程で最適クラブ長さを決定するための指標である。最適クラブ長さを決定するための指標として、このような捩り度合いを表す指標が算出される理由については、後述する。

具体的には、まず、取得部２４Ａが、記憶部２３内に格納されているｘｙｚ局所座標系でのグリップ角速度ω_x，ω_y，ω_zの時系列の計測データを読み出す（ステップＳ２１）。

次に、グリップ挙動導出部２４Ｂが、グリップ角速度ω_x，ω_y，ω_zの計測データに基づいてトップの時刻ｔ_i及びインパクトの時刻ｔ_tを導出する（ステップＳ２２）。トップの時刻ｔ_iの導出方法としては、様々な方法が公知であるが、ステップＳ２２では図７に示すようなグリップ角速度ω_yの波形において、グリップ角速度ω_yの符号が負から正へと転換するタイミング、すなわちω_y＝０となるタイミングをトップのタイミングと定義する。従って、グリップ挙動導出部２４Ｂは、ω_yの計測データから、ω_y＝０となる時刻をトップの時刻ｔ_iとして特定する。インパクトの時刻ｔ_tの導出方法としては、様々な方法が公知であるため、詳細な説明を省略する。なお、本実施形態では、特許第6059878号で説明されている方法に従い、グリップ角速度ω_x，ω_y，ω_zの計測データに基づいてインパクトの時刻ｔ_tが導出される。

ステップＳ２３では、算出部２４Ｄが、第１スイング指標ＳＷ１として、トップの時刻ｔ_iからインパクトの時刻ｔ_tまでのグリップ角速度ω_zの平均値を算出する。本実施形態では、下式に従ってグリップ角速度ω_zの平均値が算出される。

すなわち、本実施形態の第１スイング指標ＳＷ１は、トップの時刻ｔ_iからインパクトの時刻ｔ_tまでの区間でのグリップ角速度ω_zの絶対値の積分値を、トップからインパクトまでの期間で除することで算出される。しかしながら、第１スイング指標ＳＷ１は、これに限られず、トップの時刻ｔ_iからインパクトの時刻ｔ_tまでの区間でのグリップ角速度ω_zの絶対値の積分値であってもよい。算出された第１スイング指標ＳＷ１がＲＡＭまたは記憶部２３に保存されると、第１指標算出工程は終了する。

＜２－３．最適クラブ長さ決定工程＞
以下、図８を参照しつつ、第１指標算出工程で算出された第１スイング指標ＳＷ１に基づいて、最適クラブ長さを決定する最適クラブ長さ決定工程（Ｓ３）について説明する。本実施形態では、第１スイング指標ＳＷ１の大きさに応じて、基本クラブ長さ、基本クラブ長さよりも短い第１クラブ長さ、及び基本クラブ長さよりも長い第２クラブ長さのいずれかが最適クラブ長さとして決定される。各クラブ長さは、予め定められており、基本クラブ長さは、これに限定されないが、例えば４５．２５インチである。また、第１クラブ長さは、これに限定されないが、例えば４４．７５インチであり、第２クラブ長さは、これに限定されないが、例えば４５．７５インチである。なお、クラブ長さの計測方法は、シャフトＤＢ２９に登録されている多数のシャフト４０間で統一されていれば、６０度法に基づいて定義されてもよいし、ヒールエンド法に基づいて定義されてもよい。本実施形態のクラブ長さの計測方法は、ヒールエンド法に基づいている。

ステップＳ３１では、決定部２４Ｅが、第１スイング指標ＳＷ１が所定の第１閾値ＴＨ１以上（ＳＷ１≧ＴＨ１）であるか否かを判定する。第１閾値ＴＨ１は、多数の実験データ等に基づいて予め定められた値であり、記憶部２３内に格納されていてもよいし、プログラム３に組み込まれていてもよい。本実施形態では、ＴＨ１＝２５０と定められている。決定部２４Ｅが、ＳＷ１≧ＴＨ１であると判定すると、処理はステップＳ３２に進む。一方、決定部２４Ｅが、第１スイング指標ＳＷ１が所定の第１閾値ＴＨ１未満、つまりＳＷ１＜ＴＨ１であると判定すると、処理はステップＳ３３に進む。

ステップＳ３２では、決定部２４Ｅが、最適クラブ長さを、基本クラブ長さよりも短い第１クラブ長さに決定する。その後、最適クラブ長さ決定工程は終了する。

ステップＳ３３では、決定部２４Ｅが、第１スイング指標ＳＷ１が所定の第２閾値ＴＨ２以下（ＳＷ１≦ＴＨ２）であるか否かを判定する。第２閾値ＴＨ２は、第１閾値ＴＨ１と同様に、多数の実験データ等に基づいて予め定められた値であり、記憶部２３内に格納されていてもよいし、プログラム３に組み込まれていてもよい。本実施形態では、ＴＨ２＝７５と定められている。決定部２４Ｅが、ＳＷ１≦ＴＨ２であると判定すると、処理はステップＳ３４に進む。一方、決定部２４Ｅが、第１スイング指標ＳＷ１が所定の第２閾値ＴＨ２を超えていると判定すると、処理はステップＳ３５に進む。

ステップＳ３４では、決定部２４Ｅが、最適クラブ長さを、基本クラブ長さよりも長い第２クラブ長さに決定する。その後、最適クラブ長さ決定工程は終了する。

第１スイング指標ＳＷ１が第２閾値ＴＨ２を超え、第１閾値ＴＨ１未満（ＴＨ２＜ＳＷ１＜ＴＨ１）である場合の処理は、ステップＳ３５の処理となる。ステップＳ３５では、決定部２４Ｅが、最適クラブ長さを、基本クラブ長さに決定する。その後、最適クラブ長さ決定工程は終了する。

以上のように、第１スイング指標ＳＷ１に基づいて最適クラブ長さを決定可能であることの裏付けは、発明者らの行った実験により確認された。発明者らは、多数のゴルファーにセンサユニット１付きゴルフクラブ４を試打させ、グリップ角速度ω_zの時系列のデータを取得した。すると、図９に示すように、ダウンスイング中に捻りのローテーションの動きを多く使用するゴルファーと、捻りのローテーションの動きを殆ど使用しないゴルファーとではグリップ角速度ω_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をスイング平面Ｐ内へ射影したグリップ速度ｖ_pX，ｖ_pY，ｖ_pZ及びグリップ角速度ω_pX，ω_pY，ω_pZを算出する。具体的な算出方法については、特許文献１で詳細に説明されているため、ここでは説明を省略する。

＜２－６．肩挙動導出工程＞
肩挙動導出工程（Ｓ６）では、スイング平面Ｐ内でのグリップの挙動（グリップ速度Ｖ_GE及びグリップ角速度ω_pX）に基づいて、スイング平面Ｐ内の疑似的な肩の挙動を導出する処理が行われる。本実施形態では、ゴルフクラブ４の挙動は、ゴルファーＧの肩及びグリップ４２（或いは、これを握るゴルファーの手首）を節点とし、ゴルファーＧの腕及びゴルフクラブ４をリンクとする、図１１に示す二重振り子モデルに基づいて解析される。肩挙動導出部２４Ｃは、肩の挙動として、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω₁を算出する。具体的な算出方法については、特許文献１で詳細に説明されているため、ここでは説明を省略する。

＜２－７．第２指標算出工程＞
第１指標算出工程（Ｓ７）では、グリップ４２の挙動及び肩の挙動に基づいて、第２スイング指標を算出する処理が行われる。第２スイング指標とは、最適振り易さ指標を決定するための、ゴルファーＧによるスイング動作を特徴付ける特徴量である。本実施形態では、算出部２４Ｄが、第２スイング指標として、腕出力パワーＰ₁__AVE、クラブ入力パワーＰ₂__AVE及びヘッド速度Ｖ_hを算出する。ここで、第２スイング指標は、図１２に示されるスイング平面Ｐ内での新たなＸＹ座標系に準拠して算出される。図１２の紙面は、スイング平面Ｐに等しい。スイング平面Ｐ内での新たなＸＹ座標系のＸ軸は、上述したＸＹＺ全体座標系のＹ軸に等しく、新たなＸＹ座標系のＹ軸は、ＸＹＺ全体座標系のＺ軸をスイング平面Ｐ内に投影した軸である。腕出力パワーＰ₁__AVE、クラブ入力パワーＰ₂__AVE及びヘッド速度Ｖ_hの具体的な算出方法については、特許文献１や本出願人による先行技術文献である特開２０１７－２１７３２４号公報で詳細に説明されているため、ここでは説明を省略する。

特開２０１７－２１７３２４号公報でも説明されている通り、腕出力パワーＰ₁__AVEとは、スイング動作中にゴルファーＧがコックを溜める強さを表す指標であると言い換えることができる。また、クラブ入力パワーＰ₂__AVEとは、スイング動作中にゴルファーＧがコックを解放する強さを表す指標であると言い換えることができる。

＜２－８．最適振り易さ決定工程＞
最適振り易さ決定工程（Ｓ８）では、第２指標算出工程で算出された第２スイング指標に基づいて、最適振り易さ指標が決定される。具体的な最適振り易さ指標の決定手順については、特許文献１で詳細に説明されているため、以下、参考のために簡単に説明する。ここで、最適振り易さ指標とは、最適スイングＭＩであり、最適スイングＭＩとは、ゴルファーに適したゴルフクラブ４のスイング慣性モーメントＩ_Sのことを言う。スイング慣性モーメントＩ_sとは、スイング中の肩回りの慣性モーメントであり、本実施形態では、以下の式に従って定義される。
Ｉ_s＝Ｉ₂＋ｍ₂（Ｒ＋Ｌ）²＋Ｉ₁＋ｍ₁（Ｒ／２）²
ただし、Ｉ₂はゴルフクラブ４の重心周りの慣性モーメントであり、ゴルフクラブ４のスペックとして予め定められているものとする。ｍ₂はゴルフクラブ４の質量である。ＲはゴルファーＧの腕長さであり、本実施形態ではＲ＝６０ｃｍ（一定）である。しかし、肩挙動導出工程（Ｓ６）における過程で算出される、スイング平面Ｐ内におけるトップからインパクトまでの肩とグリップ４２との距離Ｒ（図１１参照）としてもよい。また、Ｌはグリップ４２からゴルフクラブ４の重心までの距離であり、ゴルフクラブ４のスペックとして予め定められているものとする。さらに、ｍ₁はゴルファーＧの腕の質量であり、例えば、解析を開始する前にゴルファーＧの体重を入力しておき、入力された体重に所定の係数を掛ける等して自動的に算出されるものとする。
なお、各ゴルファーＧについては、ゴルフクラブ４が変わっても腕の重量は同じである。従って、本実施形態では簡単のため、スイング慣性モーメントＩ_sは、腕の回転分の慣性モーメントを省略し、以下の式に従って算出される。
Ｉ_s＝Ｉ₂＋ｍ₂（Ｒ＋Ｌ）²
ところで、Ｉ_Sを決定するパラメータであるｍ₂，Ｉ₂，Ｌは、ゴルフクラブ４の諸元である。このため、本実施形態のスイング慣性モーメントＩ_Sも、ゴルフクラブ４の諸元である。

決定部２４Ｅは、計測工程で試打されたテストクラブがプロモデルクラブであった場合は、腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEを示す点が、図１３に示すＰ₁__AVE－Ｐ₂__AVE空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定する。Ｐ₁__AVE－Ｐ₂__AVE空間を分割する境界線Ｌ１～Ｌ４を特定する情報、及び境界線Ｌ１～Ｌ４で分割される分割領域Ａ１～Ａ５と最適スイングＭＩ帯との対応関係を定めるデータは、対応関係データ２８として記憶部２３内に格納されている。

また、決定部２４Ｅは、計測工程で試打されたテストクラブがアベレージモデルクラブであった場合は、腕出力パワーＰ₁__AVE、クラブ入力パワーＰ₂__AVE及びヘッド速度Ｖ_hを示す点が、図１４に示す腕出力パワーＰ₁__AVE－クラブ入力パワーＰ₂__AVE－ヘッド速度Ｖ_h空間におけるどの領域にプロットされるかに応じて、最適スイングＭＩ帯を判定する。ただし、簡単のため、図１４では、ヘッド速度Ｖ_hを表す軸は省略され、腕出力パワーＰ₁__AVE－クラブ入力パワーＰ₂__AVE平面が示されている。Ｐ₁__AVE－Ｐ₂__AVE空間を分割する境界線Ｌ５～Ｌ６を特定する情報、及び境界線Ｌ１～Ｌ６で分割される分割領域Ｂ１～Ｂ３と最適スイングＭＩ帯との対応関係を定めるデータは、対応関係データ２８として記憶部２３内に格納されている。なお、図２では、対応関係データ２８は、フィッティングプログラム３とは別のデータとして示されているが、プログラム３内に組み込まれていてもよい。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

第２絞り込みシャフトとしては、１本のみ特定されてもよいし、複数本が特定されてもよい。また、第２絞り込みシャフトとしては、第１絞り込みシャフトの中で相対的に一致度が高い所定数のシャフトが特定されてもよいし、一定以上の一致度を有する全てのシャフトが特定されてもよい。

続いて、選択部２４Ｆは、ヘッドＤＢ２７及びシャフトＤＢ２９内に含まれる情報から、各第２絞り込みシャフトと推奨ヘッドとを組み合わせた場合のゴルフクラブ４のクラブ長さを算出し、その値が最適クラブ長さ決定工程（Ｓ３）で決定された最適クラブ長さに合致するか、最も近くなるようなシャフトを特定する。なお、クラブ長さは、通常ヘッド側の基準点からグリップのグリップキャップラインまでの長さを含むため、選択部２４Ｆは、各第２絞り込みシャフトと推奨ヘッドとを組み合わせた場合のゴルフクラブ４のクラブ長さとして、グリップ部分による長さを考慮したクラブ長さを算出する。選択部２４Ｆは、特定されたシャフトを推奨シャフトとして選択する。また、選択部２４Ｆは、推奨シャフトと推奨ヘッドとを組み合わせた場合のゴルフクラブ４を推奨ゴルフクラブとして選択する。

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

＜３．特徴＞
本実施形態によるシステム１００によれば、共通のセンサユニット１による計測データに基づいて、例えばゴルファーＧの嗜好するヘッド４１を含む、ゴルファーＧに適した推奨ゴルフクラブが選択される。推奨ゴルフクラブは、最適クラブ長さ、最適振り易さ指標及び最適剛性指標に合致するか、一致度が高いゴルフクラブである。すなわち、クラブ長さ、スイングＭＩ帯及びＩＦＣの観点からゴルファーＧに適したゴルフクラブが選択されるため、精度の高いフィッティングを行うことができる。

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

＜４－１＞
上記実施形態では、最適クラブ選択工程（Ｓ１１）において、まずスイング慣性モーメントＩ_Sの値が最適スイングＭＩ帯に属するゴルフクラブまたはこれに含まれるシャフトが絞られ、次に最適剛性指標に合致するシャフトが絞られ、最終的に最適クラブ長さに合致するゴルフクラブまたはこれに含まれるシャフトが特定された。しかし、ゴルフクラブまたはシャフトの絞り込みは、この順序で行われなくてもよい。例えば、最初に最適クラブ長さに合致するゴルフクラブまたはこれに含まれるシャフトが絞り込まれ、その後に最適剛性指標に合致するシャフトが絞られ、あるいはスイング慣性モーメントＩ_Sの値が最適スイングＭＩ帯に属するゴルフクラブまたはこれに含まれるシャフトが絞られてもよい。

＜４－２＞
上記実施形態のシステム１００では、クラブ長さのフィッティングと振り易さのフィッティング及び剛性のフィッティングが合わせて行われた。しかし、システム１００は、振り易さのフィッティング及び剛性のフィッティングを省略し、単にクラブ長さのみをフィッティングするシステムとしても構成することができる。つまり、システム１００は、工程Ｓ１～Ｓ３で第１スイング指標ＳＷ１に基づいて最適クラブ長さを決定した後、最適クラブ選択工程（Ｓ１１）で最適クラブ長さに合致するゴルフクラブを選択するように構成されてもよい。この場合のゴルフクラブの選択は、例えば上記実施形態と同様にまず推奨ヘッドを決定し、推奨ヘッドと組み合わせた場合のゴルフクラブ４のクラブ長さが最適クラブ長さに合致するか、最も近くなるようなシャフト４０を推奨シャフトとして特定することにより行うことができる。また例えば、記憶部２３に予め固定されたヘッド４１とシャフト４０との組み合わせのデータを格納しておき、最適クラブ長さに合致するか、最も近くなるようなヘッド４１とシャフト４０とが組み合わせられたゴルフクラブを工程Ｓ１１で選択するようにしてもよい。さらに、第１スイング指標ＳＷ１に基づくクラブ長さのフィッティングは、振り易さのフィッティング及び剛性のフィッティングのうち一方と、あるいは別の方法によるフィッティングと組み合わされてもよい。

＜４－３＞
上記実施形態では、最適振り易さ決定工程（Ｓ８）において、最終的に最適スイングＭＩ帯が決定された。しかしながら、最適振り易さ決定工程（Ｓ８）では、最適スイングＭＩから、ゴルファーＧに適したゴルフクラブ４の特定の部位の特性を表す最適特性指標が決定されてもよい。例えば、最終的にゴルファーＧに適したゴルフクラブ４のシャフト４０の重量（最適シャフト重量）が最適特性指標として決定され、最適クラブ選択工程（Ｓ１１）において、最適シャフト重量に基づいて最適シャフトが絞り込まれてもよい。この場合、第２スイング指標（Ｐ₁__AVE、Ｐ₂__AVE、及びＶ_h）を説明変数とし、最適スイングＭＩを目的変数とする重回帰式を予め算出し、記憶しておくことで、第２スイング指標から最適スイングＭＩを直接算出してもよい。最適スイングＭＩと、最適シャフト重量帯との対応関係を対応関係データとして予め定めておくことにより、対応関係データを照合し、算出した最適スイングＭＩから最適シャフト重量帯を決定することができる。ここで、ゴルフクラブ４のスイング慣性モーメントＩ_Sとは、ヘッド４１のスイング慣性モーメントと、シャフト４０のスイング慣性モーメントと、その他の部位（グリップ４２、フェラル等）のスイング慣性モーメントとの合計値である。従って、その他の部位のスイング慣性モーメントが記憶部２３に記憶されている場合には、これらのデータを参照し、ゴルフクラブ４全体としてのスイング慣性モーメントＩ_Sが最適化されるよう、推奨シャフト及び推奨ヘッドを選択することができる。その他、シャフト重量帯のみならず、ゴルファーＧに適したシャフト４０及びグリップ４２を合わせた重量帯、ヘッド４１分を除くゴルフクラブ４の重量帯等が最適スイングＭＩとの対応から最適特性指標として決定されてもよい。すなわち、ゴルフクラブ４の特定の部位はシャフト４０、シャフト４０及びグリップ４２、並びにヘッド４１を除くゴルフクラブ４の部位であり得る。また、対応関係データは、ヘッド４１の種類ごとに用意されてもよいし、ヘッド４１の種類に依存しないデータであってもよい。

また、上記実施形態では、最適剛性決定工程（Ｓ１０）において、最適ランク値Ｋ_S1～Ｋ_S4（ＩＦＣ）が決定されたが、これに代えて、ゴルファーＧに適したゴルフクラブ４のシャフト４０のフレックス（最適フレックス）が決定されてもよい。最適シャフト重量及び最適フレックスの決定方法、並びに最適シャフト重量及び最適フレックスに基づいたシャフト４０の選択方法については、特許文献１で詳細に説明されているため、ここでは説明を省略する。最終的には、最適シャフト重量及び最適フレックスに合致するシャフト４０の中から、最適クラブ長さが実現されるようなシャフト４０が推奨シャフトとして選択されてもよい。あるいは、最適クラブ長さに合致するゴルフクラブ４の中から、最適シャフト重量及び最適フレックスに合致するシャフト４０を含むものが推奨ゴルフクラブとして選択されてもよい。

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

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

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

＜４－７＞
また、最適振り易さ指標としては、上述したクラブ重量ｍ₂、グリップエンド周りの慣性モーメントＩ_G、及び重心周りの慣性モーメントＩ₂等、最適スイングＭＩに限らず様々な指標を設定することができるが、これらの最適振り易さ指標との間に一定の関係（相関）が認められる限り、任意の指標を第２スイング指標とすることができる。例えば、スイング動作時にゴルファーにより発揮される腕エネルギーＥ_AVE、Ｅ₁、総肩トルクＴ_ti及び平均肩トルクＴ_AVEを第２スイング指標として用いることができる。Ｔ_ti及びＴ_AVEは、スイング動作時にゴルファーにより発揮される肩周りのトルクを表す指標である。Ｅ_AVE、Ｅ₁、Ｔ_ti及びＴ_AVEの算出方法は、特許文献１で詳細に説明されている。また、第２スイング指標として腕エネルギーＥ_AVE、平均肩トルクＴ_AVE及びヘッド速度Ｖ_hが算出され、これらの指標Ｅ_AVE、Ｔ_AVE、Ｖ_hに応じて最適クラブ重量の範囲（以下、最適重量帯）または最適スイングＭＩ帯が算出される場合の手順については、特許文献１で詳細に説明されているため、ここでは説明を省略する。

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

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

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

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

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

＜４－９＞
上記実施形態では、最適振り易さ指標決定工程（Ｓ８）の後に最適剛性指標決定工程（Ｓ１０）が行われたが、シャフト４０及びゴルフクラブ４の絞り込み方法はこれに限定されない。例えば、第２スイング指標の大きさに応じて、最適シャフト重量帯を直接算出することもできる。例えば、図１９に示すように、（好ましくは、ヘッド４１の種類毎に）Ｐ₁__AVE－Ｐ₂__AVE（－Ｖ_h）空間を領域分割して、最適シャフト重量帯に対応する領域を定義しておき、当該定義情報を図２の対応関係データ２８に代わる対応関係データとして予め記憶しておくことができる。この場合、第２指標算出工程で導出される第２スイング指標がどの領域に属するかに応じて、最適特性指標を導出することができる。

また、対応関係データとして、最適フレックスごとに、第２スイング指標空間（例えば、Ｐ₁__AVE－Ｐ₂__AVE平面）を各最適シャフト重量帯に分割する境界線を特定するデータが予め算出され、記憶されていてもよい。そして、最適クラブ長さによるシャフト４０またはゴルフクラブ４の絞り込みは、最適フレックスに合致し、最適シャフト重量帯に属するシャフト４０の絞り込み前に行われてもよいし、後に行われてもよい。

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

＜４－１０＞
上記実施形態や変形例において、フィッティングの精度をより向上させるべく、第２スイング指標が、第２スイング指標を表す空間を分割することにより定義される複数の領域の境界付近に存在する場合には、既出の最適スイングＭＩ帯や最適シャフト重量を修正する境界処理が行われてもよい。最適シャフト重量は、例えば最適シャフト重量帯に存在するとともに、最適フレックスに合致するか、最も近いフレックスを有するシャフトの重量であり、このようなシャフトが１本存在する場合は、そのシャフトの重量を最適シャフト重量とすることができる。また、このようなシャフトが複数本存在する場合は、例えばヘッド４１の種類を固定した場合に、ゴルフクラブ４全体としてのスイング慣性モーメントＩ_Sが最適化されるようなシャフトの重量とすることもできる。例えば、第２スイング指標が、Ｐ₁__AVE－Ｐ₂__AVE（－Ｖ_h）空間内のいずれかの分割領域の上限境界線と閾値以内の距離に存在する場合、決定部２４Ｅは、既出の最適スイングＭＩ帯を一段階上の領域の最適スイングＭＩ帯に修正し、または既出の最適シャフト重量を、所定量（例えば、２．５ｇ）だけ増量するように調整してもよい。また、第２スイング指標が、Ｐ₁__AVE－Ｐ₂__AVE（－Ｖ_h）空間内のいずれかの分割領域の下限境界線と閾値以内の距離に存在する場合、決定部２４Ｅは、既出の最適スイングＭＩ帯を一段階下の領域の最適スイングＭＩ帯に修正し、または既出の最適シャフト重量を所定量、所定量（例えば、２．５ｇ）だけ減量するように調整してもよい。

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

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

以下では、一例として、図１９のように最適シャフト重量帯が決定される場合に適用される例外処理について説明する。具体的な処理の流れは、図２２に示すとおりである。図２２に示す例外処理は、最適シャフト重量が、図１９に示す判定基準（分割領域Ｃ１～Ｃ５の境界線Ｌ７～Ｌ１０）に従って決定された後に実行される。

ステップＳ７１では、算出部２４Ｄが、ゴルフクラブ４のシャフト軸周りの捻りのローテーションの動きを表す指標（捩り指標）を導出する。捻り指標は、既に第１指標算出工程（Ｓ２）が実行されていれば、そこで算出された第１スイング指標ＳＷ１を使用することができ、第１指標算出工程が未だ実行されていない場合、上述した第１スイング指標ＳＷ１の定義に従って算出することができる。続いて、決定部２４Ｅが、ヘッド速度Ｖ_hが所定範囲内であり、かつ、捩り指標が所定値以上であるかを判定する。そして、かかる条件が満たされる場合には、既出の最適シャフト重量を所定量だけ増量するように調整する（ステップＳ７２）。一方、かかる条件が満たされない場合には、処理はステップＳ７３に進む。

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

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

以上のステップが完了すると、例外処理は終了する。本例外処理と、変形例４－１０の境界処理とは、いずれも最適振り易さ指標及び最適シャフト重量帯をより精度よく決定するための処理であり、両処理は、組み合わせることもできる。この場合、境界処理の後、本例外処理を実行することが好ましい。

＜４－１２＞
上記実施形態では、第１スイング指標算出工程（Ｓ２）で算出された第１スイング指標ＳＷ１と第１閾値ＴＨ１及び第２閾値ＴＨ２とが比較され、比較結果に応じて最適クラブ長さが基本クラブ長さ、第１クラブ長さ、及び第２クラブ長さのいずれかに決定された。しかし、決定される最適クラブ長さは２種類であってもよく、４種類以上であってもよく、これに応じて第１スイング指標ＳＷ１の閾値の数も増減してよい。また、上述した最適クラブ長さ決定工程（Ｓ３）において、ステップＳ３１及びステップＳ３２、あるいはステップＳ３３及びステップＳ３４のいずれかが省略されてもよい。具体的には、決定部２４Ｅは、第１スイング指標ＳＷ１が第１閾値ＴＨ１以上である場合、最適クラブ長さを第１クラブ長さと決定し、それ以外の場合、最適クラブ長さを基本クラブ長さと決定するように構成されてもよい。また、決定部２４Ｅは、第１スイング指標ＳＷ１が第２閾値ＴＨ２以下である場合、最適クラブ長さを第２クラブ長さと決定し、それ以外の場合、最適クラブ長さを基本クラブ長さと決定するように構成されてもよい。

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

＜実施例＞
複数のゴルファーに対し、特許文献１のフィッティング方法（比較例）を適用するとともに、上記実施形態に係るフィッティング方法（実施例）を適用した。

比較例の方法では、スイング慣性モーメントＩ_Sが各ゴルファーの最適スイングＭＩ帯に属するとともに、最適剛性指標に最も合致したＩＦＣのシャフトを有するゴルフクラブのうち１本が、推奨ゴルフクラブとして選択された。推奨ゴルフクラブのクラブ長さは、いずれも基本クラブ長さ（４５．２５インチ）であった。そして、推奨ゴルフクラブを用いて、各ゴルファーに何球かの試打を行わせ、打球の飛距離及び左右ズレを計測した。

実施例の方法では、最適スイングＭＩ帯及び最適剛性指標に加え、上記実施形態の通り、各ゴルファーに適した最適クラブ長さも決定された。そして、最適クラブ長さが第１クラブ長さ（４４．７５インチ）または第２クラブ長さ（４５．７５インチ）であると決定されたゴルファーについては、スイング慣性モーメントＩ_Sが各ゴルファーの最適スイングＭＩ帯に属するとともに、最適剛性指標に最も合致したＩＦＣのシャフトを有し、さらにクラブ長さが最適クラブ長さにカスタマイズされたゴルフクラブを推奨ゴルフクラブとした。そして、推奨ゴルフクラブを用いて、ゴルファーに何球かの試打を行わせ、打球の飛距離及び左右ズレを計測した。

図２３Ａ～Ｄは、実施例で基本クラブ長さ以外のクラブ長さが推奨された各ゴルファーについて、実施例及び比較例による推奨ゴルフクラブの飛距離及び左右ズレの計測結果を比較する図である。図２３Ａは、実施例において、第１スイング指標ＳＷ１が第１閾値である２５０ｄｅｇ/ｓ以上であったため、第１クラブ長さのゴルフクラブが推奨されたゴルファーの例である。図２３Ａから分かるように、比較例による推奨ゴルフクラブでも飛距離は充分であったところ、基本クラブ長さの推奨ゴルフクラブでは、ボールの方向が左にずれる傾向があり、ボールが捕まりすぎる傾向があることを示唆していた。この点、実施例による推奨ゴルフクラブでは、飛距離は一定水準を保ちながらボールの方向性が安定し、左右ズレが改善される結果となった。図２３Ｂは、実施例において、第１クラブ長さのゴルフクラブが推奨された別のゴルファーの例である。図２３Ｂから分かるように、この例では、実施例によるゴルフクラブを用いた試打で、最大飛距離が大幅に増大した。

図２４Ｃは、実施例において、第１スイング指標ＳＷ１が第２閾値である７５ｄｅｇ/ｓ以下であったため、第２クラブ長さのゴルフクラブが推奨されたゴルファーの例である。図２４Ｃから分かるように、比較例による推奨ゴルフクラブと実施例による推奨ゴルフクラブとでは、左右ズレの程度に殆ど差はなかったものの、実施例によるゴルフクラブを用いた試打で、最大飛距離が大幅に増大した。図２４Ｄは、実施例において、第２クラブ長さのゴルフクラブが推奨された別のゴルファーの例である。図２４Ｄから分かるように、実施例によるゴルフクラブを用いた試打では、左右ズレが比較的安定し、最大飛距離が大幅に増大した。

以上の結果から、実施例の方法により選択された推奨ゴルフクラブでは、比較例の方法により選択された推奨ゴルフクラブと比較して、より好ましいスイングが可能となることが分かった。つまり、実施例の方法によれば、ゴルファーにより適したゴルフクラブを推奨し得ることが確認された。

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

Claims (11)

1. ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得する取得部と、
   前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標を算出する算出部と、
   前記第１スイング指標の大きさに応じて、前記ゴルファーに適した前記ゴルフクラブの長さである最適クラブ長さを決定する決定部と
   を備える、
   フィッティング装置。
2. 前記決定部は、前記第１スイング指標を所定の第１閾値と比較し、前記第１スイング指標が前記第１閾値以上である場合には、前記最適クラブ長さが第１クラブ長さであると決定する、
   請求項１に記載のフィッティング装置。
3. 前記決定部は、前記第１スイング指標を前記第１閾値よりも小さい所定の第２閾値と比較し、前記第１スイング指標が前記第２閾値以下である場合には、前記最適クラブ長さが前記第１クラブ長さよりも長い第２クラブ長さであると決定する、
   請求項２に記載のフィッティング装置。
4. 前記計測データには、前記ゴルフクラブのグリップエンドにおける角速度が含まれ、前記第１スイング指標は、前記ゴルフクラブのシャフトに平行な軸周りの角速度に基づいて算出される、
   請求項１から３のいずれかに記載のフィッティング装置。
5. 前記第１スイング指標は、ダウンスイング時における前記角速度の平均値または積分値である、
   請求項４に記載のフィッティング装置。
6. 前記最適クラブ長さに最も合致するゴルフクラブ及び当該ゴルフクラブに含まれるべきシャフトの少なくとも一方を選択する選択部、
   をさらに備える、
   請求項１から５のいずれかに記載のフィッティング装置。
7. ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得する取得部と、
   前記計測データに基づいて、前記スイング動作に関する第１スイング指標、第２スイング指標及び第３スイング指標を算出する算出部と、
   前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定し、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの振り易さ指標である最適振り易さ指標を決定し、前記第３スイング指標の大きさに応じて、前記ゴルファーに適した前記シャフトの剛性を示す最適剛性指標を決定する決定部と、
   前記最適振り易さ指標、前記最適剛性指標及び前記最適クラブ長さに最も合致するゴルフクラブ及び当該ゴルフクラブに含まれるべきシャフトの少なくとも一方を選択する選択部と
   を備え、
   前記第１スイング指標は、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する指標である、
   フィッティング装置。
8. ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得する取得部と、
   前記計測データに基づいて、前記スイング動作に関する第１スイング指標、第２スイング指標及び第３スイング指標を算出する算出部と、
   前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定し、前記第２スイング指標の大きさに応じて、前記ゴルファーに適したゴルフクラブの特定の部位の特性を表す最適特性指標を決定し、前記第３スイング指標の大きさに応じて、前記ゴルファーに適した前記シャフトの剛性を示す最適剛性指標を決定する決定部と、
   前記最適特性指標、前記最適剛性指標、及び前記最適クラブ長さに最も合致するゴルフクラブ及び当該ゴルフクラブに含まれるべきシャフトの少なくとも一方を選択する選択部と
   を備え、
   前記第１スイング指標は、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する指標である、
   フィッティング装置。
9. 前記計測データには、前記ゴルフクラブのグリップエンドにおける角速度が含まれ、前記第１スイング指標は、前記ゴルフクラブのシャフトに平行な軸周りの角速度に基づいて算出される、
   請求項７または８に記載のフィッティング装置。
10. ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得するステップと、
    前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標をコンピュータを用いて算出するステップと、
    前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定するステップと、
    を含む、
    フィッティング方法。
11. ゴルファーによるゴルフクラブのスイング動作を計測機器により計測した計測データを取得するステップと、
    前記計測データに基づいて、前記スイング動作中の前記ゴルフクラブの捩り度合いに関する第１スイング指標を算出するステップと、
    前記第１スイング指標の大きさに応じて、前記ゴルファーに適したクラブ長さである最適クラブ長さを決定するステップと、
    をコンピュータに実行させる、
    フィッティングプログラム。

Images