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

Description

本発明は、ゴルファーに適したゴルフクラブの重量（以下、最適クラブ重量）を決定するフィッティング装置、方法及びプログラムに関する。

従来より、ゴルファーにテストクラブを試打させてその動作を計測機器により計測し、当該計測値に基づいてゴルファーに適したゴルフクラブを決定する様々なフィッティング方法が提案されている。特許文献１では、ゴルファーに適したゴルフクラブを選択するために、ゴルフクラブの重量やゴルフクラブの長さ等を考慮することが重要であることが指摘されている。

特開２０１３−２２６３７５号公報

ところで、最適クラブ重量は、そのゴルファーが現在使っているゴルフクラブの重量や、様々なゴルフクラブの試打結果や見た目のスイングのテンポ等に基づいて、熟練のクラブフィッターにより経験的に決定されることが多い。しかしながら、このような経験や勘に依存したフィッティング方法は客観的ではなく、フィッティングの結果に個人差が生じるという問題がある。

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

本発明の第１観点に係るフィッティング装置は、ゴルファーに適したゴルフクラブの重量である最適クラブ重量を決定するフィッティング装置であって、取得部と、算出部と、決定部とを備える。前記取得部は、前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する。前記算出部は、前記計測値に基づいて、前記スイング動作時に前記ゴルファーにより発揮されるエネルギー又はトルクを表す指標、或いはこれらと相関のある指標を算出する。前記決定部は、前記指標の大きさに応じて、前記最適クラブ重量を決定する。

本発明の第２観点に係るフィッティング装置は、第１観点に係るフィッティング装置であって、前記算出部は、複数の種類の前記指標を算出する。前記決定部は、前記複数の種類の指標の大きさに応じて、前記最適クラブ重量を決定する。

本発明の第３観点に係るフィッティング装置は、第１観点又は第２観点に係るフィッティング装置であって、前記指標には、前記スイング動作時の前記ゴルファーの腕エネルギー及び前記スイング動作時の前記ゴルファーの肩周りのトルクの少なくとも一方が含まれる。

本発明の第４観点に係るフィッティング装置は、第３観点に係るフィッティング装置であって、前記算出部は、前記スイング動作時のヘッドスピードを算出する。前記決定部は、前記腕エネルギー及び前記肩周りのトルクの少なくとも一方の大きさに加え、前記ヘッドスピードの大きさに応じて、前記最適クラブ重量を決定する。

本発明の第５観点に係るフィッティング装置は、第１観点から第４観点のいずれかに係るフィッティング装置であって、前記決定部は、前記指標が大きい又は小さい程、前記最適クラブ重量を大きな値に決定する。

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

本発明の第７観点に係るフィッティング装置は、第１観点から第６観点のいずれかに係るフィッティング装置であって、前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定する最適クラブ特定部をさらに備える。

本発明の第８観点に係るフィッティング方法は、ゴルファーに適したゴルフクラブの重量である最適クラブ重量を決定するフィッティング方法であって、以下の（１）〜（３）のステップを備える。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測するステップ。
（２）前記計測機器による計測値に基づいて、前記スイング動作時に前記ゴルファーにより発揮されるエネルギー又はトルクを表す指標、或いはこれらと相関のある指標を算出するステップ。
（３）前記指標の大きさに応じて、前記最適クラブ重量を決定するステップ。

本発明の第９観点に係るフィッティング方法は、第８観点に係るフィッティング方法であって、以下の（４）のステップをさらに備える。
（４）前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定するステップ。

本発明の第１０観点に係るフィッティングプログラムは、ゴルファーに適したゴルフクラブの重量である最適クラブ重量を決定するフィッティングプログラムであって、以下の（１）〜（３）のステップをコンピュータに実行させる。
（１）前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップ。
（２）前記計測機器による計測値に基づいて、前記スイング動作時に前記ゴルファーにより発揮されるエネルギー又はトルクを表す指標、或いはこれらと相関のある指標を算出するステップ。
（３）前記指標の大きさに応じて、前記最適クラブ重量を決定するステップ。

本発明の第１１観点に係るフィッティングプログラムは、第１０観点に係るフィッティングプログラムであって、以下の（４）のステップをさらにコンピュータに実行させる。
（４）前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定するステップ。

本発明者らは、最適クラブ重量と、スイング動作時にゴルファーにより発揮されるエネルギー又はトルクとの間には、相関関係があることを発見した。そこで、本発明によれば、ゴルファーによるテストクラブのスイング動作時に当該ゴルファーにより発揮されるエネルギー又はトルクを表す指標、或いはこれらと相関のある指標の大きさに応じて、最適クラブ重量が決定される。これにより、最適クラブ重量を客観的に精度よく決定することができる。

本発明の一実施形態に係るフィッティング装置を備えるゴルフスイングシステムを示す図。 フィッティングシステムの機能ブロック図。 ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。 （Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。 第１変換処理の流れを示すフローチャート。 アドレスの時刻を導出する処理の流れを示すフローチャート。 スイング平面を説明する図。 肩挙動導出工程の流れを示すフローチャート。 二重振り子モデルを概念的に説明する図。 指標算出工程の流れを示すフローチャート。 二重振り子モデルを概念的に説明する別の図。 重量指標の１つである腕エネルギーを説明する図。 最適総重量決定工程の流れを示すフローチャート。 プロモデルユーザ１０名及びアベレージモデルユーザ１０名によるスイング動作時の腕エネルギー及び平均肩トルクの値をプロットした図。 プロモデル領域及びアベレージモデル領域を示す図。 ゴルフクラブの重量とヘッド速度との関係を示す図。 実験により求めた、プロモデルユーザによるスイング動作時のヘッド速度及び最適クラブ重量を示す図。 実験により求めた、アベレージモデルユーザによるスイング動作時のヘッド速度及び最適クラブ重量を示す図。 プロモデルユーザに最適クラブ重量を割り当てるための分割領域を示す図。 アベレージモデルユーザに最適クラブ重量を割り当てるための分割領域を示す図。 シミュレーションによるヘッド速度の等高線図。 別のシミュレーションによるヘッド速度の等高線図。 さらに別のシミュレーションによるヘッド速度の等高線図。 シミュレーションによるコック角の等高線図。 最適クラブ特定工程を説明する図。 変形例に係る重量指標である各種角度を説明する図。

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

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

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

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

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

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

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

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

記憶部２３は、ハードディスク等の不揮発性の記憶装置により構成される。記憶部２３内には、フィッティングプログラム３が格納されている他、センサユニット１から送られてくる計測データが保存される。また、記憶部２３には、対応関係データ２８が格納されている。対応関係データ２８とは、詳細は後述するが、様々なモデルのゴルフクラブ４毎に規定されており、最適クラブ重量を決定するための条件を示すデータである。通信部２５は、フィッティング装置２と外部装置との通信を可能にする通信インターフェースであり、センサユニット１からデータを受信する。

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

＜２．スイング動作の解析処理＞
続いて、フィッティングシステム１００による、ゴルフクラブ４のフィッティングのためのスイング動作の解析処理について説明する。本実施形態に係る解析処理は、以下の７つの工程から構成されている。
（１）ｘｙｚ局所座標系でのグリップ加速度ａ_x，ａ_y，ａ_z、グリップ角速度ω_x，ω_y，ω_z及びグリップ地磁気ｍ_x，ｍ_y，ｍ_zの計測データを計測する計測工程
（２）計測工程で得られたｘｙｚ局所座標系での計測データを、ＸＹＺ全体座標系でのグリップ加速度ａ_X，ａ_Y，ａ_Z及びグリップ角速度ω_X，ω_Y，ω_Zに変換する第１変換工程（第１変換工程では、ＸＹＺ全体座標系でのグリップ速度ｖ_X，ｖ_Y，ｖ_Zも導出される。）
（３）ＸＹＺ全体座標系でのグリップ４２の挙動（グリップ角速度ω_X，ω_Y，ω_Z及びグリップ速度ｖ_X，ｖ_Y，ｖ_Z）を、スイング平面Ｐ（後述する）内でのグリップ４２の挙動へと変換する第２変換工程
（４）スイング平面Ｐ内でのグリップ４２の挙動に基づいて、スイング平面Ｐ内でのゴルファー７の疑似的な肩の挙動を導出する肩挙動導出工程
（５）グリップ４２の挙動及び疑似的な肩の挙動に基づいて、最適クラブ重量を決定するための仕様となる重量指標（本実施形態では、ゴルファー７により発揮される平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hに関する３つの指標）を算出する指標算出工程
（６）重量指標に基づいて、最適クラブ重量を決定する最適総重量決定工程
（７）ゴルファー７に適したシャフト４０の重量（以下、最適シャフト重量）及び曲げ剛性（以下、最適曲げ剛性）を算出する最適シャフト決定工程
（８）最適総重量決定工程で特定された最適クラブ重量に合致するとともに、最適シャフト決定工程で算出された最適シャフト重量及び最適曲げ剛性に合致する複数本のゴルフクラブの中から、特にヘッド速度を高めることができるゴルフクラブを特定する最適クラブ特定工程
以下、これらの工程を順に説明する。

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

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

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

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

また、計測工程では、以上のゴルフクラブ４が複数回、好ましくは５回以上試打されることが好ましい。この場合、計測データの平均値を算出し、以降の演算に使用することができる。また、ミスショットや計測ミス等による異常値を取り除くため、計測データの標準偏差σを算出し、試打全ての計測データが好ましくは平均値±１．６５σ以内に、より好ましくは平均値±１．２８σ以内に収まるような計測データを得ることが好ましい。そして、当該チェックを行うべく、制御部２４により計測データの標準偏差σを算出するようにし、このσの値が以上の条件を満たしていない場合には、計測の追加又はやり直しを求めるメッセージを表示部２１上に表示させるようにしてもよい。なお、計測データ自体の平均値ではなく、計測データに基づいて算出される加工値（例えば、後述する平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_h）の平均値を算出するようにしてもよい。加工値の平均値を算出する場合も、同じく標準偏差σに基づくデータの信頼性のチェックを行うことができる。

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

次に、ステップＳ１で読み出されたｘｙｚ局所座標系での時系列の計測データに基づいて、グリップ挙動導出部２４Ｂが、インパクト、トップ及びアドレスの時刻ｔ_i，ｔ_t，ｔ_aを導出する（ステップＳ２）。本実施形態では、まずインパクトの時刻ｔ_iが導出され、インパクトの時刻ｔ_iに基づいてトップの時刻ｔ_tが導出され、トップの時刻ｔ_tに基づい
てアドレスの時刻ｔ_aが導出される。

具体的には、グリップ角速度ω_xのサンプリング周期Δｔ当たりの増分が閾値である３００ｄｅｇ／ｓを最初に超えた時刻が、仮のインパクトの時刻として設定される。そして、この仮のインパクトの時刻から所定の時間を溯った時刻から、仮のインパクトの時刻までで、グリップ角速度ω_xのサンプリング周期Δｔ当たりの増分が２００ｄｅｇ／ｓを超えた時刻が検出され、インパクトの時刻ｔ_iとして設定される。

次に、インパクトの時刻ｔ_iよりも前の時刻であって、グリップ角速度ω_yが負から正へ切り替わった時刻が、トップの時刻ｔ_tとして特定される。また、アドレスの時刻ｔ_aは、図６のフローチャートに従って算出される。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

＜２−５．指標算出工程＞
以下、図１０を参照しつつ、グリップ４２の挙動及び肩の挙動に基づいて、最適クラブ重量を決定するための重量指標を算出する指標算出工程について説明する。本実施形態では、重量指標として、後述される平均肩トルクＴ_AVE、腕エネルギーＥ_AVE及びヘッド速度Ｖ_hが算出される。

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

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

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

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

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

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

次に、ステップＳ３３では、指標算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでのグリップ４２に発生する拘束力Ｒ₂＝（Ｒ_X2，Ｒ_Y2）を算出する。以下の式は、並進方向の力の釣り合いに基づくものである。ただし、ｍ₂は、ゴルフクラブの質量であり、ｇは、重力加速度である。また、ｍ₂は、ゴルフクラブ４のスペックであり、予め定められているものとする。

続くステップＳ３４では、指標算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの肩周りのトルクＴ₁及びグリップ４２周りのトルクＴ₂を算出する。

ただし、Ｉ₁は、腕の重心周りの慣性モーメントであり、Ｉ₂は、ゴルフクラブ４の重心周りの慣性モーメントである。本実施形態では、腕の重心周りの慣性モーメントＩ₁は、
腕の重心が腕の中心にあるとの仮定の下、Ｉ₁＝ｍ₁ｒ²／３として算出される。ｍ₁は、腕の質量であり、本実施形態では、腕の質量ｍ₁は、適宜予め定められているものとする。
例えば、解析を開始する前に、ゴルファー７の体重を入力しておき、入力された体重に所定の係数を掛ける等して、自動的に腕の質量が算出される。また、Ｉ₂は、ゴルフクラブ４のスペックであり、予め定められているものとする。

また、本実施形態では、指標算出部２４Ｄは、肩周りのトルクＴ₁をトップからインパクトまでの区間で積分した値Ｔ_tiを算出する。Ｔ_tiは、トップからインパクトまでの間にゴルファー７の肩周りで発揮されるトルクの総量（以下、総肩トルク）であり、この意味でスイング動作中の肩周りのトルクを表す指標である。また、指標算出部２４Ｄは、トルクＴ_tiをトップからインパクトまでの時間で除算したスイング動作中の平均的な肩周りのトルク（以下、平均肩トルク）Ｔ_AVE＝Ｔ_ti／ｔ_i−ｔ_cを算出する。平均肩トルクＴ_AVEも、スイング動作中の肩周りのトルクを表す指標となり、重量指標の１つとなる。なお、総肩トルクＴ_tiを算出するに当たり、正のトルクＴ₁だけを積分してもよいし、トルクＴ₁の平均値を積分してもよい。

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

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

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

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

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

＜２−６．最適総重量決定工程＞
以下、図１３を参照しつつ、最適クラブ重量を決定する最適総重量決定工程の流れについて説明する。

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

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

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

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

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

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

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

なお、最適クラブ重量は、以下のアルゴリズムに従って、シミュレーション上の最適クラブ重量として算出することも可能である。すなわち、本発明者らは、ヘッド速度Ｖ_hとゴルフクラブ４の重量との間には、概ね図１６に示すような関係があることを発見した。具体的に説明すると、スイングされるゴルフクラブ４が重くなる程、ゴルファーの力ではゴルフクラブ４を振り切ることができなくなり、ヘッド速度Ｖ_hが小さくなる。かといって、ある一定以上軽くしても、ヘッド速度Ｖ_hは頭打ちとなる（図１６参照）。フルスイング以上の力では振ることができないからである。言い換えると、ある一定以上軽く振り易くしても、ゴルファーの限界に達すれば、ヘッド速度Ｖ_hが伸びなくなる。従って、ヘッド速度Ｖ_hは、ある点（図１６の変化点）を境にして、ヘッド速度Ｖ_hがゴルフクラブ４の重量に比例する比例領域と、ヘッド速度Ｖ_hがゴルフクラブ４の重量に関わらず概ね一定となる一定領域とに分かれる。また、飛距離を伸ばすためには、ヘッド速度Ｖ_h及びゴルフクラブ４の重量は共に大きい方がよいため、以上の変化点に対応するゴルフクラブ４の重量が、飛距離を最大にする最適クラブ重量であると言うことができる。従って、上記実験に参加するゴルファーに様々な重量のゴルフクラブ４をスイングさせて、ヘッド速度Ｖ_hとゴルフクラブ４の重量との関係をヘッド速度Ｖ_h−ゴルフクラブの重量平面内にプロットする。そして、比例領域における回帰直線及び一定領域における回帰直線を算出し、その交点である変化点を求め、当該変化点に対応するゴルフクラブ重量を、最適ゴルフクラブ重量とすることもできる。

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

ところで、プロモデル領域では、Ｔ_AVE≧８５Ｎ・ｍ、かつ、Ｅ_AVE≧３５０Ｎ・ｍ／ｓの場合には、ヘッド速度Ｖ_hが４５ｍ／ｓか否かにより、最適重量帯が３２０ｇ以上か３１０ｇ〜３２０ｇかのいずれかに振り分けられるようになっている。これは、図１７に示す実験結果に示されるとおり、Ｔ_AVE≧８５Ｎ・ｍ、かつ、Ｅ_AVE≧３５０Ｎ・ｍ／ｓの領域の境界線付近においては、Ｖ_h≦４５ｍ／ｓの場合の最適クラブ重量が３１６ｇで、Ｖ_h≧４５ｍ／ｓの場合の最適クラブ重量が３２４ｇ以上になっているためである。

＜２−７．最適シャフト決定工程＞
以下、ゴルファー７に適した最適シャフト重量及び最適曲げ剛性を算出する最適シャフト決定工程について説明する。

まず、決定部２４Ｅは、最適重量帯に基づいて、最適シャフト重量の範囲（以下、最適シャフト重量帯）を決定する。具体的には、ゴルフクラブ４のモデル毎の各最適重量帯には、以下の表１及び２とおり、予め最適シャフト重量帯が定められており、これに基づいて最適シャフト重量帯が決定される。

続いて、決定部２４Ｅは、シャフト４０の最適曲げ剛性の範囲（以下、最適曲げ剛性帯）を決定する。最適曲げ剛性帯の決定方法は、公知であるため（要すれば、特開２０１３−２０８３６６号公報参照）、ここでは詳細な説明を省略する。

以上のステップにより、最適重量帯、最適シャフト重量帯及び最適曲げ剛性帯が決定されると、決定部２４Ｅは、ゴルフクラブ４全体での重量、シャフト４０の重量及び曲げ剛性が既知の様々なゴルフクラブ（以下、対象クラブ）の中から、最適重量帯、最適シャフト重量帯及び最適曲げ剛性帯に属するゴルフクラブ（以下、候補クラブ）を特定する。表示制御部２４Ｆは、特定された候補クラブの種類を示す情報とともに、最適重量帯、最適シャフト重量帯及び最適曲げ剛性帯を表示部２１上に表示させる。これにより、ユーザは、ゴルファー７に適したゴルフクラブの種類を知ることができるとともに、最適重量帯、最適シャフト重量帯及び最適曲げ剛性帯を知ることができる。

＜２−８．最適クラブ特定工程＞
以下、上記の工程で特定された候補クラブの中から、特にヘッド速度を高めることができるゴルフクラブ（以下、最適クラブ）を特定する最適クラブ特定工程について説明する。最適クラブは、フィッティングを行っているゴルファー７が各候補クラブを用いた場合のグリップエンド慣性モーメントＩ_G及びスイング慣性モーメントＩ_Sの値に基づいて特定される。

なお、グリップエンド慣性モーメントＩ_Gとは、グリップエンド周りの慣性モーメントであり、Ｉ_G＝Ｉ₂＋ｍ₂Ｌ²として算出される。一方、スイング慣性モーメントＩ_Sは、スイング中の肩周りの慣性モーメントであり、以下の式に従って算出可能である。
Ｉ_S＝Ｉ₂＋ｍ₂（２ｒ＋Ｌ）²＋Ｉ₁＋ｍ₁ｒ²

ただし、各ゴルファー７については、ゴルフクラブが変わっても腕の重量は同じである。従って、簡単のため、本実施形態におけるスイング慣性モーメントＩ_Sは、腕の回転分の慣性モーメントを省略し、以下の式に従って定義される。
Ｉ_S＝Ｉ₂＋ｍ₂（２ｒ＋Ｌ）²

以下に説明する最適クラブ特定工程のアルゴリズムは、ヘッド速度Ｖ_h、グリップエンド慣性モーメントＩ_G及びスイング慣性モーメントＩ_Sが、図２１〜図２３に示すような関係を有することに基づく。図２１〜図２３は、本発明者らが行ったシミュレーションの結果である。図２１中のＩ_G−Ｉ_S平面内の１３個の点（丸）は、被験者がスペックの異なる１３本のゴルフクラブをスイングしたときの（Ｉ_G，Ｉ_S）を示す点である。ただし、本データは、実際に１３本のゴルフクラブを被験者に試打させて得たデータではなく、１本のゴルフクラブ（以下、基準クラブ）を試打させ、このときのデータからシミュレーションにより得たデータである。

基準クラブのスペック（ゴルフクラブ重量ｍ₂[ｇ]，重心位置Ｌ[ｍｍ]，慣性モーメントＩ₂[ｋｇ・ｃｍ²]）は、（２７２，９３６，４５３）であった。また、基準クラブには、上記計測工程で用いられたようなセンサユニット１が取り付けられた。そして、当該基準クラブを被験者に試打させ、上述したアルゴリズムと同様のアルゴリズムに従って逆動力学解析を行い、肩周りのトルクＴ₁、グリップ４２周りのトルクＴ₂及び腕長さＲを算出した。続いて、これらのパラメータＴ₁，Ｔ₂，Ｒを一定と仮定した上で、１３本のゴルフクラブのスペック（ゴルフクラブ重量ｍ₂、重心距離Ｌ、慣性モーメントＩ₂等）の値を用いて、順運動力学解析を行った。なお、図２１のシミュレーションでは、ゴルフクラブ重量ｍ₂＝２７２ｇ（一定）のゴルフクラブが用いられ、基準クラブを除く残りの１２本のゴルフクラブのスペック（Ｌ，Ｉ₂）は、（９４１，１９２）（９４１，１３８）（９３８．５，３２４）（９３６，５０７）（９４１，８３）（９３８．５，２７０）（９３３．５，６３３）（９３１，８０９）（９３６，３９９）（９３３．５，５７８）（９３１，７５５）（９３１，７００）とされた。以上の順運動力学解析により、各ゴルフクラブによるヘッド速度Ｖ_h、グリップエンド慣性モーメントＩ_G及びスイング慣性モーメントＩ_Sが算出された。図２１中に引かれている１０本の斜線は、１３個のヘッド速度Ｖ_hの値に基づいて導出されたヘッド速度Ｖ_hの等高線である。

同様のシミュレーションを異なる１３本のゴルフクラブに対し行った結果が、図２２及び図２３である。具体的には、図２２のシミュレーションでは、重心位置Ｌ＝９３６ｃｍ（一定）のゴルフクラブが用いられ、基準クラブを除く残りの１２本のゴルフクラブのスペック（ｍ₂，Ｉ₂）は、（２８４，３４９）（２８４，２９４）（２７８，４００）（２７２，５０７）（２８４，２４０）（２７８，３４６）（２６６，５６０）（２６０，６６６）（２７２，３９９）（２６６，５０６）（２６０，６１２）（２６０，５５７）とされた。また、図２３のシミュレーションでは、慣性モーメントＩ₂＝４５３ｋｇ・ｃｍ²（一定）のゴルフクラブが用いられ、基準クラブを除く残りの１２本のゴルフクラブのスペック（Ｌ，ｍ₂）は、（９３２．６，２９３）（９３０．９，２９７）（９３４．３，２８２）（９３７．８，２６７．６）（９２９．２，３０２）（９３２．５，２８７）（９３３．５，２５８）（９４３．５，２４４）（９３４．２，２７６．４）（９３７．８，２６２）（９４１．５，２４８）（９３９．７，２５２）とされた。

以上の図２１〜図２３の等高線図は、０．１ｍ／ｓ刻みで引かれている。これらの等高線に添えられている数値は、基準クラブでのヘッド速度Ｖ_hを０．０ｍ／ｓとしたときの相対的なヘッド速度Ｖ_hの値である。これらの等高線図からは、Ｉ_G−Ｉ_S平面内では、右下にいくにつれてヘッド速度Ｖ_hが向上し、左上にいくほど低下することが分かる。言い換えると、スイング慣性モーメントＩ_Sの値が小さくなり、グリップエンド慣性モーメントＩ_Gが大きくなる程、ヘッド速度Ｖ_hが向上することが分かる。この傾向は、ゴルフクラブ重量ｍ₂が一定であるか、重心位置Ｌが一定であるか、慣性モーメントＩ₂が一定であるか等の条件には依存しない。

ところで、一般に、グリップエンド慣性モーメントＩ_Gやスイング慣性モーメントＩ_Sが増大すると、重心の位置がヘッド側に近づくためゴルフクラブが振りにくくなり、ヘッド速度Ｖ_hが低下すると考えられている。しかしながら、本発明者らは、図２１〜図２３のシミュレーション結果から、たとえグリップエンド慣性モーメントＩ_Gが増大したとしても、当該増分に対するスイング慣性モーメントＩ_Sの増分を一定値以下とすれば、むしろヘッド速度Ｖ_hを向上させることができることに気が付いた。発明者らはさらに検討を進め、このことは、図２４に示す別のシミュレーションの結果から説明できることが分かった。

図２４は、図２１のシミュレーションにおいて、被験者が１３本のゴルフクラブをスイングしたときの１３個のコック角の値に基づいて導出されたコック角の等高線図である。なお、ここでいうコック角とは、コック解放タイミングｔ_rにおける腕とゴルフクラブとの為す角度（図１１に一定鎖線で示される角度）である。

図２４の等高線は、０．５°刻みで引かれている。これらの等高線に添えられている数値は、基準クラブでのコック角を０．０°としたときの相対的な角度である。これらの等高線図からは、Ｉ_G−Ｉ_S平面内では、右へいくにつれてコック角が小さくなり、左にいくほど大きくなることが分かる。言い換えると、グリップエンド慣性モーメントＩ_Gが大きくなる程、コック角が小さくなることが分かる。一方で、等高線が概ね上下方向に延びていることから、コック角はスイング慣性モーメントＩ_Sには影響を受けない。なお、簡単のため、図２１に対応するシミュレーションの結果のみを示したが、図２２及び図２３に対応するシミュレーションにおいても、コック角に関し同様の傾向が確認された。

また、コック角が小さいとは、リストコックが溜まっており、スイング中にゴルフクラブがゴルファーの胴体の近くを通ることを意味する。従って、コック角が小さい場合には、実効的なスイング慣性モーメントＩ_Sが小さくなり、ヘッド速度Ｖ_hの上昇が期待される。

以上より、グリップエンド慣性モーメントＩ_Gが増大したとしても、スイング慣性モーメントＩ_Sの増分を一定値以下とすれば、ゴルフクラブが振りにくくなることによるデメリットよりも、コック角が小さくなることによるメリットが勝り、ヘッド速度Ｖ_hが向上することが分かる。すなわち、グリップエンド慣性モーメントＩ_Gを大きくし、スイング慣性モーメントＩ_Sを小さくすることができれば、ヘッド速度Ｖ_hを向上させることができる。

以上の知見に基づいて、最適クラブ特定工程が実行される。まず、最適クラブ特定部２４Ｇは、フィッティングの対象となる複数本の対象クラブの中から、候補クラブを絞り込む。具体的には、記憶部２３内には、各対象クラブのスペックを示す情報（以下、スペック情報）が予め格納されている。本実施形態では、ここでいうスペックとして、各対象クラブに対し、ゴルフクラブ重量ｍ₂、ゴルフクラブ４の重心周りの慣性モーメントＩ₂、グリップ４２からゴルフクラブ４の重心までの距離（重心距離）Ｌ等の値が格納されている。従って、最適クラブ特定部２４Ｇは、記憶部２３内の当該スペック情報を参照することにより、対象クラブの中から、候補クラブとして、最適重量帯、最適シャフト重量帯及び最適曲げ剛性帯に属するゴルフクラブを全て特定する。なお、候補クラブが１本しかない場合には、以下の処理は省略され、当該一本の候補クラブが、最適クラブとして特定される。

一方、候補クラブが複数本存在する場合には、これらの候補クラブの中から特にヘッド速度を高めることが可能なゴルフクラブである、最適クラブが特定される。具体的には、最適クラブ特定部２４Ｇは、上述した定義式に従って、フィッティングを行っているゴルファー７が各候補クラブをスイングしたときの、グリップエンド慣性モーメントＩ_G及びスイング慣性モーメントＩ_Sを導出する。このとき、計測工程から指標算出工程までの工程で既に算出されているパラメータＲと、候補クラブのスペック（ゴルフクラブ重量ｍ₂、重心距離Ｌ、慣性モーメントＩ₂等）の値を用いて、順運動力学解析が実行される。この場合、再度ゴルファー７に、候補クラブをスイングしてもらう必要はない。

各候補クラブに対応する慣性モーメントＩ_G，Ｉ_Sが判明すると、最適クラブ特定部２４Ｇは、複数本の候補クラブの中から、スイング慣性モーメントＩ_Sがより小さく、グリップエンド慣性モーメントＩ_Gのより大きいゴルフクラブを特定し、これを最適クラブとする。具体的には、最適クラブ特定部２４Ｇは、各候補クラブに対応する（Ｉ_G，Ｉ_S）の点をＩ_G−Ｉ_S平面内にプロットし、最も右下にくる点に対応する候補クラブを最適クラブと判定する。例えば、図２５（Ａ）に示すような３本の候補クラブＡ，Ｂ，Ｃが存在する場合には、候補クラブＣが最適クラブと判定される。一方、図２５（Ｂ）のような３本の候補クラブＡ，Ｂ，Ｃが存在する場合には、Ｂ，Ｃの候補クラブに関しては、どちらが右下に存在しているのか判定することができない。この場合には、Ｂ，Ｃの両方が最適クラブとして判定される。

なお、図２５（Ｂ）のＢ，Ｃに対応する複数本の候補クラブも、Ｉ_G−Ｉ_S平面内に図２１〜図２３に示したような等高線を引くことができれば、優劣を判断することができる。従って、別の実施形態では、最適クラブ特定部２４Ｇが、既に算出済みのパラメータＴ₁，Ｔ₂，Ｒと、様々な所定のゴルフクラブのスペック（ゴルフクラブ重量ｍ₂、重心距離Ｌ、慣性モーメントＩ₂等）の値を用いて、順動力学解析により、Ｉ_G−Ｉ_S平面内におけるヘッド速度Ｖ_hの等高線を導出するようにしてもよい。なお、この等高線は、ゴルファー７に依存するものである。そして、この等高線の傾きから、Ｂ，Ｃのような複数本の候補クラブに対応するヘッド速度Ｖ_hの高低を判断し（図２５（Ｃ）参照）、最も高いものを最適クラブと判定することができる。

以上の処理が終わると、表示制御部２４Ｆは、最適クラブを特定する情報を表示部２１上に表示させる。これにより、ゴルファー７は、自身に最適なゴルフクラブを把握することができる。また、以上の出力値の説得性を向上させるべく、表示制御部２４Ｆは、図２５（Ａ）〜（Ｃ）に示すグラフも併せて表示部２１上に表示させることもできる。
＜３．変形例＞
以上、本発明の一実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。

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

＜３−２＞
ヘッド速度Ｖ_hは、上述した重回帰式に基づく統計的手法ではなく、他の方法で算出することもできる。例えば、以下の式に従って幾何学的に算出してもよい。ただし、Ｌ_clubは、ゴルフクラブのスペックであるゴルフクラブの長さである。
シャフト４０の先端のヘッドの位置ベクトル（ｄ_hX，ｄ_hY）
ｄ_hX＝２Ｘ₁＋Ｌ_clubｃｏｓθ₂
ｄ_hY＝２Ｙ₁＋Ｌ_clubｓｉｎθ₂
シャフト４０の先端のヘッドの速度ベクトル（Ｖ_hX，Ｖ_hY）
Ｖ_hX＝２Ｖ_X1−Ｌ_clubω₂ｓｉｎθ₂
Ｖ_hY＝２Ｖ_Y1＋Ｌ_clubω₂ｃｏｓθ₂
Ｖ_h＝ｓｑｒｔ（Ｖ_hX ²＋Ｖ_hY ²）

＜３−３＞
最適総重量決定工程は、フィッティング装置２ではなく、ユーザ自身が実行することができる。すなわち、指標算出工程の終了後に、算出された重量指標を表示部２１上に表示させる。そして、図１９及び図２０に示すような対応関係データ２８に基づいて、ユーザが自ら最適重量帯を判断してもよい。このとき、対応関係データ２８を表示部２１上に表示させたり、或いは対応関係データ２８を印刷した紙媒体を用意したりすると、ユーザの判断が容易となる。同様に、最適シャフト決定工程及び／又は最適クラブ特定工程における処理の一部又は全部も、フィッティング装置２に代えて、ユーザ自身が実行することができる。

＜３−４＞
上記実施形態で示された重量指標は、例示であり、最適クラブ重量と相関があるその他の様々な指標を用いることができる。例えば、腕エネルギーＥ_AVE及び平均肩トルクＴ_AVEに限らず、腕エネルギーＥ₁及び肩トルクＴ_tiを重量指標とすることができる。また、スイング動作時にゴルファー７により発揮される任意の部位のエネルギー又は任意の部位のトルクを表す指標を、重量指標とすることができる。また、ここでいう「スイング動作時」とは、上記実施形態の中で説明したトップからインパクトないしコック解放タイミングまでの時間に限られず、重量指標と最適クラブ重量とに相関がある限り、スイング動作中の任意の時刻又は期間を意味するものということができる。また、ゴルファー７により発揮されるエネルギー又はトルクと相関があると考えられる以下の指標も、重量指標として用いることができる。なお、以下の（７）〜（９）以外の指標は、値が大きい程、最適クラブ重量が大きくなる傾向があり、（７）〜（９）の指標は、値が小さい程、最適クラブ重量が大きくなる傾向がある。
（１）トップの時刻でシャフト４０と全体座標系Ｚ軸（グリップよりも下方）とが為す角度θ３（図２６参照）
（２）スイング動作中の角速度ω₂の平均値
（３）トップからインパクトまでの間での角速度ω₂の最大値
（４）トップからインパクトまでのグリップ速度Ｖ_GEの平均値
（５）トップからインパクトまでのグリップ速度Ｖ_GEの最大値
（６）トップからインパクトまでの間でのグリップ４２の移動距離Ｄ
（７）コック解放タイミングｔ_rとインパクトの時間との差分（ここでいうコック解放タイミングｔ_rは、腕とシャフト４０が為すコック角θ４の解放スピードが速まり、腕のエネルギーがシャフト４０のエネルギーへと変わり始めるタイミングと定義することができる。）
（８）コック解放タイミングｔ_rでの腕とシャフト４０が為すコック角θ４（図２６参照）
（９）ダウンスイング時間、すなわち、トップからインパクトまでの時間
（１０）トップの時刻からグリップ４２周りのトルクＴ₂が正負逆転する時刻までのトルクＴ₂の積分値

＜３−５＞
上記実施形態では、３つの重量指標の大きさに応じて、最適クラブ重量を決定したが、１つ、２つ、又は４つ以上の重量指標の大きさに応じて、最適クラブ重量を決定してもよい。例えば、腕エネルギーＥ_AVE、平均肩トルクＴ_AVE及びヘッド速度Ｖ_hの３つから選択される任意の１つ又は２つの指標の大きさのみに応じて、最適クラブ重量を決定することができる。なお、図１７及び図１８に示す実験の結果からは、腕エネルギーＥ_AVE及び平均肩トルクＴ_AVEは、ヘッド速度Ｖ_hよりも最適クラブ重量との相関が高いことが分かる。従って、ヘッド速度Ｖ_hを重量指標として用いる場合には、腕エネルギーＥ_AVE及び平均肩トルクＴ_AVEの少なくとも一方と組み合わせて用いることが好ましい。また、Ｅ_AVE，Ｔ_AVE，Ｖ_hに加え、変形例３−４で例示した重量指標の中から選択される任意の数の指標の大きさに応じて、最適クラブ重量を決定することもできる。

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

プロモデルユーザ及びアベレージモデルユーザ５名ずつに、それぞれプロモデルクラブ及びアベレージモデルクラブを試打させて平均肩トルクＴ_AVE、腕エネルギーＥ_AVE、ヘッド速度Ｖ_hを算出したところ、以下の表３，４に示す結果を得た。なお、Ｔ_AVE，Ｅ_AVEの値は、上記した計測工程から指標算出工程と同様の工程に従って算出された。また、Ｖ_hの値は、ＢＣＳシステムで実測した。ここで使用されたＢＣＳシステムとは、出願人らによる特開２０１２−１７０５４７及び特開２０１２−１７０５３２に開示される計測システムである。また、上記実施形態に係る方法で、最適重量帯を導出するとともに、上記実施形態の中で説明した実験と同様の方法で、飛距離を最大にする最適クラブ重量を導出した。

以上の結果より、プロモデルユーザに関しては、５名中Ｐ３の１名を除き、本実施例に係る最適重量帯と、実験による最適クラブ重量とが一致した。また、不一致の１名に関しても、両値に大きなずれはなかった。また、アベレージモデルユーザに関しては、５名中Ａ４の１名を除き、本実施例に係る方法での最適重量帯と、実験による最適クラブ重量とが一致した。また、不一致の１名に関しても、両値に大きなずれはなかった。これにより、上記実施形態に係るフィッティング方法の精度の高さが確認された。

また、以上のゴルファーＰ１，Ｐ２，Ｐ４，Ｐ５及びＡ１，Ａ２，Ａ５が最適クラブ重量の導出のための実験に用いた複数本のゴルフクラブのうち、最適クラブ重量のもの（以下、実験クラブ）のスペック（ゴルフクラブの長さ、ヘッドの重量、シャフトの重量、グリップの重量、ゴルフクラブの総重量ｍ₂）は、以下の表５のとおりであった。また、上記の定義式に従って、各ゴルファーが実験クラブを用いたときのスイング慣性モーメントＩ_S及びグリップエンド慣性モーメントＩ_G）も算出された。そして、これらのゴルファーＰ１，Ｐ２，Ｐ４，Ｐ５及びＡ１，Ａ２，Ａ５に、最適重量帯に属するゴルフクラブであって、実験クラブよりもグリップエンド慣性モーメントＩ_Gが大きく、スイング慣性モーメントＩ_Sの小さい別のゴルフクラブ（以下、検証クラブ）をスイングさせた。検証クラブのスペックは、表５のとおりであった。また、上記と同様、ＢＣＳシステムを用いて、検証クラブのスイング時のヘッド速度Ｖ_hを計測したところ、以下の表５に示すとおりの結果となった。なお、ゴルファーＰ３，Ａ４については、本実施例に係る方法での最適重量帯と、実験による最適クラブ重量とが一致しなかったため、検証の対象としなかった。また、Ａ３については、適切な検証クラブを用意できなかったため、検証の対象外とした。

以上の結果より、実験クラブよりも、検証クラブの方が、ヘッド速度Ｖ_hが速くなることが確認された。すなわち、上記実施形態に係る最適クラブ特定工程によれば、さらにヘッド速度Ｖ_hを高めることができるゴルフクラブを絞り込むことが可能であることが確認された。

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

Claims (8)

1. ゴルファーに最大の飛距離を与えるゴルフクラブの重量である最適クラブ重量を決定するフィッティング装置であって、
   前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得する取得部と、
   前記計測値に基づいて、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される腕の仕事量、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される肩周りのトルク、並びにこれらの少なくとも一方と相関のある指標からなる群から選択される少なくとも２つの要素を算出する算出部と、
   前記少なくとも２つの要素の大きさの組み合わせと、最大の飛距離を与える前記最適クラブ重量との対応関係を定める対応関係データを記憶する記憶部と、
   前記対応関係データを参照することにより、前記少なくとも２つの要素の大きさに応じて、前記最適クラブ重量を決定する決定部と
   を備え、
   前記対応関係データでは、前記少なくとも２つの要素の各々について、当該要素が大きい又は小さい程、前記最適クラブ重量のより大きな値が対応付けられている、
   フィッティング装置。
2. 前記算出部は、前記相関のある指標として、前記スイング動作時のヘッドスピードを算出し、
   前記決定部は、前記腕の仕事量及び前記肩周りのトルクの少なくとも一方の大きさに加え、前記ヘッドスピードの大きさに応じて、前記最適クラブ重量を決定する、
   請求項１に記載のフィッティング装置。
3. 前記記憶部は、前記テストクラブの種類毎に、前記対応関係データを記憶し、
   前記決定部は、前記テストクラブの種類に応じて、前記記憶部内の前記対応関係データを参照することにより、前記最適クラブ重量を決定する、
   請求項１又は２に記載のフィッティング装置。
4. 前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定する最適クラブ特定部
   をさらに備える、請求項１から３のいずれかに記載のフィッティング装置。
5. ゴルファーに最大の飛距離を与えるゴルフクラブの重量である最適クラブ重量を決定するフィッティング方法であって、
   前記ゴルファーによるテストクラブのスイング動作を計測機器により計測するステップと、
   前記計測機器による計測値に基づいて、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される腕の仕事量、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される肩周りのトルク、並びにこれらの少なくとも一方と相関のある指標からなる群から選択される少なくとも２つの要素を算出するステップと、
   前記少なくとも２つの要素の大きさの組み合わせと、最大の飛距離を与える前記最適クラブ重量との対応関係を定める対応関係データを参照することにより、前記少なくとも２つの要素の大きさに応じて、前記最適クラブ重量を決定するステップとを備え、
   前記対応関係データでは、前記少なくとも２つの要素の各々について、当該要素が大きい又は小さい程、前記最適クラブ重量のより大きな値が対応付けられている、
   フィッティング方法。
6. 前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定するステップ
   をさらに備える、請求項５に記載の方法。
7. ゴルファーに最大の飛距離を与えるゴルフクラブの重量である最適クラブ重量を決定するフィッティングプログラムであって、
   前記ゴルファーによるテストクラブのスイング動作を計測機器により計測した計測値を取得するステップと、
   前記計測機器による計測値に基づいて、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される腕の仕事量、ダウンスイング中の期間を含む前記スイング動作時に前記ゴルファーにより発揮される肩周りのトルク、並びにこれらの少なくとも一方と相関のある指標からなる群から選択される少なくとも２つの要素を算出するステップと、
   前記少なくとも２つの要素の大きさの組み合わせと、最大の飛距離を与える前記最適クラブ重量との対応関係を定める対応関係データを参照することにより、前記少なくとも２つの要素の大きさに応じて、前記最適クラブ重量を決定するステップと
   をコンピュータに実行させ、
   前記対応関係データでは、前記少なくとも２つの要素の各々について、当該要素が大きい又は小さい程、前記最適クラブ重量のより大きな値が対応付けられている、
   フィッティングプログラム。
8. 前記最適クラブ重量に合致する複数本のゴルフクラブの中から、スイング慣性モーメントが小さく、グリップエンド慣性モーメントの大きいゴルフクラブを特定するステップ
   をさらにコンピュータに実行させる、請求項７に記載のプログラム。