JP6981735B2 - Ｇｕｉ表示装置 - Google Patents

Description

本発明は、ゴルフクラブのスイング動作に関するＧＵＩ（グラフィカルユーザーインターフェース）画面を表示するＧＵＩ表示装置に関する。

従来より、ゴルファーによるゴルフクラブのスイング動作を計測装置により計測し、計測データに基づいてスイング動作を解析する装置が公知である（特許文献１等）。このような解析の場面では、解析のアルゴリズム自体も勿論重要であるが、解析の結果をユーザに分かり易く説明することも重要である。特許文献１には、ゴルフスイングの解析結果をグラフ等を用いて分かり易く視認できるようにしたＧＵＩ画面が開示されている。なお、このようなゴルフスイングの解析結果は、ゴルファーに適したゴルフクラブのフィッティングや、ゴルファーのフォームの改善等、様々な用途で利用される。

特開２０１４−１２１４１２号公報

ところで、ゴルフスイング中には、コックを溜める動きと、コックを解放する動きとが行われ、このような動きは、ゴルファーのスイング動作を特徴付ける指標となる。従って、本発明者らは、ゴルファー固有のこれらの動きの強さを解析し、ゴルファーに提示できるようにすることが重要であると考えた。しかしながら、コックを溜める動き及びコックを解放する動きが何を意味するのか、多くの一般のゴルファーが理解しているわけではない。

本発明は、ゴルフスイング中のコックを溜める動き及びコックを解放する動きを視覚的に分かり易く説明することができるＧＵＩ表示装置を提供することを１つの目的とする。

第１観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、表示制御部を備える。前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ₁（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ₂（ｉ）とを同時に表示する。

第２観点に係るＧＵＩ表示装置は、第１観点に係るＧＵＩ表示装置であって、前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ₁（ｉ）及び前記第２図形Ｇ₂（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って動画の態様で表示する。

第３観点に係るＧＵＩ表示装置は、第１観点に係るＧＵＩ表示装置であって、前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ₁（ｉ）及び前記第２図形Ｇ₂（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って配列する態様で表示する。

第４観点に係るＧＵＩ表示装置は、第１観点から第３観点のいずれかに係るＧＵＩ表示装置であって、前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記第１図形Ｇ₁（ｉ）及び前記第２図形Ｇ₂（ｉ）に加え、第ｉ番目のタイミングにおける前記ゴルファーの腕を表す図形である腕図形Ｇ₃（ｉ）と、第ｉ番目のタイミングにおける前記ゴルフクラブを表す図形であるクラブ図形Ｇ₄（ｉ）とを同時に表示する。

第５観点に係るＧＵＩ表示装置は、第４観点に係るＧＵＩ表示装置であって、前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記腕図形Ｇ₃（ｉ）に重ねて又はその近傍に前記第１図形Ｇ₁（ｉ）を表示し、前記クラブ図形Ｇ₄（ｉ）に重ねて又はその近傍に前記第２図形Ｇ₂（ｉ）を表示する。

第６観点に係るＧＵＩ表示装置は、第１観点から第５観点のいずれかに係るＧＵＩ表示装置であって、取得部と、算出部とをさらに備える。前記取得部は、ゴルファーによる前記ゴルフクラブのスイング動作を計測した計測データを取得する。前記算出部は、前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する。前記表示制御部は、前記ＧＵＩ画面上において、前記第１指標及び前記第２指標を表示する。

第７観点に係るＧＵＩ表示装置は、第６観点に係るＧＵＩ表示装置であって、決定部をさらに備える。前記決定部は、前記第１指標及び前記第２指標に応じて、前記ゴルファーに適したシャフト重量である最適シャフト重量を決定する。前記表示制御部は、前記ＧＵＩ画面上において、前記最適シャフト重量を表示する。

第８観点に係るＧＵＩ表示装プログラムは、ＧＵＩ画面を表示するＧＵＩ表示プログラムであって、以下の（１）のステップをコンピュータに実行させる。
（１）ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ₁（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ₂（ｉ）とを同時に表示するステップ。

第９観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、取得部と、算出部と、表示制御部とを備える。前記取得部は、ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する。前記算出部は、前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する。前記表示制御部は、前記ＧＵＩ画面上において、前記コックを溜める強さを第１軸とし、前記コックを解放する強さを第２軸とするプロット平面内にプロットされる点として、前記第１指標及び前記第２指標を表示する。前記プロット平面は、前記ゴルファーに推奨すべきシャフト重量帯に分割されている。
ＧＵＩ表示装置。

第１０観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、取得部と、決定部と、表示制御部とを備える。前記取得部は、ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する。前記決定部は、前記計測データに基づいて、インパクト時の前記ゴルフクラブに含まれるヘッドの傾きを決定する。前記表示制御部は、前記ＧＵＩ画面上において、インパクト時の前記ヘッドを表す図形である第１ヘッド図形と、インパクトの直前のタイミングにおける前記ヘッドを表す図形である第２ヘッド図形とを同時に表示する。前記第１ヘッド図形は、前記ヘッドの傾きに応じて傾くように表示される。

第１観点によれば、ゴルフスイング中の複数のタイミングにおける、コックを溜める強さを表す第１図形及びコックを解放する強さを表す第２図形が、ＧＵＩ画面上に表示される。また、同じタイミングにおける第１図形及び第２図形は、同時に表示される。従って、このようなＧＵＩ画面を視たユーザは、ゴルフスイング中のコックを溜める動き及びコックを解放する動きの意味するところを直感的に理解することができる。すなわち、このようなＧＵＩ画面により、ゴルフスイング中のコックを溜める動き及びコックを解放する動きを視覚的に分かり易く説明することができる。

本発明の一実施形態に係るＧＵＩ表示装置（フィッティング装置）を備えるフィッティングシステムを示す図。 フィッティングシステムの機能ブロック図。 ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。 （Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。 フィッティング処理の流れを示すフローチャート。 スイング平面を説明する図。 二重振り子モデルを概念的に説明する図。 肩挙動導出工程を概念的に説明する図。 肩挙動導出工程を概念的に説明する別の図。 第１指標算出工程の流れを示すフローチャート。 二重振り子モデルを概念的に説明する別の図。 腕エネルギーを説明する図。 最適シャフト重量決定工程の流れを示すフローチャート。 複数の最適シャフト重量帯に対応する領域に分割された腕出力パワー−クラブ入力パワー平面を示す図。 インターナショナル・フレックス・コード（ＩＦＣ）を説明する図。 シャフトの曲げ剛性の測定方法を説明する図。 スイング中のシャフトの曲げを説明する図。 ヘッドをフェース面側から視た図。 解析の結果を表示するＧＵＩ画面を示す図。 解析の結果を表示するＧＵＩ画面を示す別の図。 解析の結果を表示するＧＵＩ画面を示すさらに別の図。 解析の結果を表示するＧＵＩ画面を示すさらに別の図。 スイング中に発生する力を説明する画像を示す図。 続くタイミングでのスイング中に発生する力を説明する画像を示す図。 さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。 さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。 さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。 さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。 解析の結果を表示するＧＵＩ画面を示すさらに別の図。

以下、図面を参照しつつ、本発明の一実施形態に係るＧＵＩ表示装置及びプログラムについて説明する。

＜１．フィッティングシステムの概略構成＞
図１及び図２に、本実施形態に係るＧＵＩ表示装置であるフィッティング装置２を備えるフィッティングシステム１００の全体構成を示す。フィッティング装置２は、ゴルファー７に適したゴルフクラブ４を選定するのを支援するための装置である。そのために、フィッティング装置２は、ゴルファー７によるゴルフクラブ４のスイング動作を計測した計測データを取得し、当該計測データに基づいて、スイング動作の解析を行う。解析の結果は、後述する画面Ｗ１〜Ｗ５（図１９〜図２４）等を介してユーザに提示される。なお、ここでいうユーザとは、ゴルファー７自身やそのインストラクター等の、解析の結果を必要とする者の総称である。本実施形態では、スイング動作の計測を行うセンサユニットは、慣性センサユニット１及びカメラシステム５であり、フィッティング装置２は、これらの慣性センサユニット１及びカメラシステム５とともに、フィッティングシステム１００を構成する。以下、慣性センサユニット１、カメラシステム５及びフィッティング装置２の構成について説明した後、フィッティング処理の流れについて説明する。

＜１−１．慣性センサユニットの構成＞
慣性センサユニット１は、図１及び図３に示すとおり、ゴルフクラブ４のグリップ４２におけるヘッド４１と反対側の端部に取り付けられており、グリップ４２の挙動を計測する。なお、ゴルフクラブ４は、一般的なゴルフクラブであり、シャフト４０と、シャフト４０の一端に設けられたヘッド４１と、シャフト４０の他端に設けられたグリップ４２とから構成される。本実施形態に係るシャフト４０は、カーボン製のシャフトである。慣性センサユニット１は、スイング動作の妨げとならないよう、小型且つ軽量に構成されている。慣性センサユニット１は、ゴルフクラブ４の外側に取り付けることができ、ゴルフクラブ４に対して着脱自在に構成することができる。

図２に示すように、本実施形態に係る慣性センサユニット１には、加速度センサ１１、角速度センサ１２及び地磁気センサ１３が搭載されている。また、慣性センサユニット１には、これらのセンサ１１〜１３から出力されるセンサデータ（計測データ）を外部のフィッティング装置２に送信するための通信装置１０も搭載されている。なお、本実施形態では、通信装置１０は、スイング動作の妨げにならないように無線式であるが、ケーブルを介して有線式にフィッティング装置２に接続するようにしてもよい。

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

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

＜１−２．カメラシステムの構成＞
次に、カメラシステム５の構成について説明する。カメラシステム５は、カメラ５１と、複数台のストロボ５３，５４とを備えており、ストロボ式の撮影を行う。カメラ５１は、ゴルフクラブ４のヘッド４１のクラウン部４１ｂ（図１８参照）を上方から撮影できるように、ゴルファー７の正面側において、支持台５７に固定されており、アドレス時のボールの位置の斜め上方に配置されている。ストロボ５３，５４も、支持台５７に固定されており、地面よりも上方で、カメラ５１の下方に配置されている。

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

また、カメラシステム５は、投光器５５Ａ，５５Ｂ及び受光器５６Ａ，５６Ｂを備えており、投光器５５Ａ及び受光器５６Ａが１つのタイミングセンサを構成しており、投光器５５Ｂ及び受光器５６Ｂがもう１つのタイミングセンサを構成している。これらのタイミングセンサにより生成される時刻は、後述するとおり、ストロボ５３，５４の発光及びそれに続くカメラ５１での撮影を行うタイミングを決定するのに使用される他、インパクト直前（インパクト時とみなし得る）のヘッド速度Ｖ_hiを算出するのに使用される。

さらに、カメラシステム５は、以上の装置５１，５３〜５６Ｂの動作を制御するための制御装置５０を備えている。制御装置５０は、ＣＰＵ，ＲＯＭ，ＲＡＭ等の他、通信部５０Ａ（図２参照）を備えており、通信部５０Ａは、以上の装置５１，５３〜５６Ｂに接続されている。また、通信部５０Ａは、フィッティング装置２の通信部２５（図２参照）にも接続されている。

投光器５５Ａ，５５Ｂは、ゴルファー７の正面側の地面付近において、カメラ５１の下方に配置されている。一方、受光器５６Ａ，５６Ｂは、ゴルファー７の足のつま先付近に配置されている。投光器５５Ａ及び受光器５６Ａは、Ｘ軸に概ね平行な直線上に配置されており、互いに対向している（図１参照）。投光器５５Ｂ及び受光器５６Ｂについても同様である。投光器５５Ａ，５５Ｂは、ゴルフスイング中、常時それぞれ受光器５６Ａ，５６Ｂに向けて光を照射しており、受光器５６Ａ，５６Ｂがこれを受光する。しかしながら、ゴルフクラブ４が投光器５５Ａ，５５Ｂと受光器５６Ａ，５６Ｂとの間を通過するタイミングでは、投光器５５Ａ，５５Ｂからの光がゴルフクラブ４により遮断されるため、受光器５６Ａ，５６Ｂはこれを受光することができない。受光器５６Ａ，５６Ｂはこのタイミングを検出し、これを受けて、制御装置５０がそれぞれタイミングｔ１，ｔ２を生成する。

制御装置５０は、このタイミングｔ１の後のタイミングｔ３において、ストロボ５３に発光を命令するとともに、カメラ５１に撮影を命令する。また、制御装置５０は、タイミングｔ２よりも後のタイミングｔ４において、ストロボ５４に発光を命令するとともに、カメラ５１に撮影を命令する。カメラ５１により撮影された画像データ（計測データ）は、制御装置５０に送信され、制御装置５０からさらにフィッティング装置２に送信される。また、制御装置５０は、タイミングｔ１，ｔ２の情報（計測データ）も、フィッティング装置２に送信する。

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

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

記憶部２３は、ハードディスク等の不揮発性の記憶装置により構成される。記憶部２３内には、プログラム３が格納されている他、センサユニット１，５から送られてくる計測データが保存される。また、記憶部２３内には、対応関係データ２８、ヘッドデータベース（ＤＢ）２７及びシャフトデータベース（ＤＢ）２９が格納されている。ヘッドＤＢ２７は、多数のヘッド４１のスペックを示す情報が、ヘッド４１の種類を特定する情報に関連付けて格納されたデータベースである。シャフトＤＢ２９は、多数のシャフト４０のスペックを示す情報が、シャフト４０の種類を特定する情報に関連付けて格納されたデータベースである。対応関係データ２８については、後述する。

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

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

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

＜２−１．計測工程＞
計測工程（Ｓ１）では、ゴルファー７により、上述の慣性センサユニット１付きゴルフクラブ４がスイングされる。以下、計測工程で使用されるゴルフクラブ４を、テストクラブと呼ぶことがある。そして、テストクラブのスイング動作中のグリップ加速度ａ_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の時系列データを積分することにより、アドレスからインパクトまでのＸＹＺ全体座標系でのグリップ速度ｖ_X，ｖ_Y，ｖ_Zを導出する。なお、局所座標系から全体座標系への値の変換方法については、様々知られている。従って、ここでは詳細な説明を省略するが、必要であれば、同出願人らによる特開２０１６−２４２９号公報や特開２０１６−２４３０号公報等に記載の方法に従うことができる。

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

また、グリップ挙動導出部２４Ｂは、グリップ速度（ｖ_pY，ｖ_pZ）を積分することにより、スイング平面Ｐ内でのグリップ４２の軌道を算出する。さらに、グリップ挙動導出部２４Ｂは、スイング平面Ｐに直交する軸周りのグリップ角速度ω_pXを算出する。なお、第２変換工程の具体的な計算方法は、適宜選択することができるが、必要であれば、同出願人らによる特開２０１６−２４２９号公報や特開２０１６−２４３０号公報等に記載の方法に従うことができる。

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

グリップの挙動から肩の挙動を特定するに当たり、本実施形態に係る二重振り子モデルは、以下の（１）〜（３）を前提とする。図７は、以下の前提条件を概念的に説明する図である。
（１）スイング平面Ｐ上において、グリップ４２（手首）は肩を中心として円運動する。
（２）スイング平面Ｐ上において、肩とグリップ４２との距離（半径）Ｒは、一定である。
（３）肩は、スイング動作中は動かない（ただし、回転する）。

以上の前提の下、肩挙動導出部２４Ｃは、第２変換工程で得られたスイング平面Ｐ内でのグリップ４２の軌道を、円弧（円）に近似する（図８参照）。そして、当該円弧（円）の中心を肩の位置Ｐ_s＝（Ｐ_sX，Ｐ_sY）とし、当該円弧（円）の中心からグリップ４２の軌道までの平均的な距離を、腕長さ（肩とグリップ４２との距離）Ｒとする。以下、グリップ４２の軌道上の点Ａ_iの座標を、（Ｘ_i，Ｙ_i）と表す（ｉ＝１，２，・・・）。なお、ここでいうｉは、時系列に沿ったデータ番号である。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

なお、以上の腕出力パワーＰ₁，Ｐ₁__AVEとは、スイング動作中にゴルファー７がコックを溜める強さを表す指標であると言い換えることができる。また、クラブ入力パワーＰ₂，Ｐ₂__AVEとは、スイング動作中にゴルファー７がコックを解放する強さを表す指標であると言い換えることができる。

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

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

＜２−６．最適シャフト重量決定工程＞
以下、図１３を参照しつつ、最適シャフト重量決定工程（Ｓ６）の流れについて説明する。最適シャフト重量決定工程は、腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEの大きさに応じて、最適シャフト重量の範囲（以下、最適シャフト重量帯）を決定するステップである。ここでは、Ｐ₁__AVE，Ｐ₂__AVEの値が大きい程、最適シャフト重量帯は段階的に大きな値に設定される。

まず、ステップＳ４１では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ₁__AVE，Ｐ₂__AVE）で表される点が、図１４に示すＰ₁__AVE−Ｐ₂__AVE平面において直線Ｌ１の上側にあるか、すなわち、図１４の領域Ａ１に属するか否か（以下、条件１）を判定する。そして、条件１が満たされる場合には、８０ｇ以上が最適シャフト重量帯であると判定する。一方、ステップＳ４１で条件１が満たされない場合には、ステップＳ４２に進む。ステップＳ４２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ₁__AVE，Ｐ₂__AVE）で表される点が、図１４に示すＰ₁__AVE−Ｐ₂__AVE平面において直線Ｌ２の上側かつ直線Ｌ１よりも下側にあるか、すなわち、図１４の領域Ａ２に属するか否か（以下、条件２）を判定する。そして、条件２が満たされる場合には、７０ｇ以上、８０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４２で条件２が満たされない場合には、ステップＳ４３に進む。ステップＳ４３では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ₁__AVE，Ｐ₂__AVE）で表される点が、図１４に示すＰ₁__AVE−Ｐ₂__AVE平面において直線Ｌ３の上側かつ直線Ｌ２よりも下側にあるか、すなわち、図１４の領域Ａ３に属するか否か（以下、条件３）を判定する。そして、条件３が満たされる場合には、６０ｇ以上、７０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４３で条件３が満たされない場合には、ステップＳ４４に進む。ステップＳ４４では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ₁__AVE，Ｐ₂__AVE）で表される点が、図１４に示すＰ₁__AVE−Ｐ₂__AVE平面において直線Ｌ４の上側かつ直線Ｌ３よりも下側にあるか、すなわち、図１４の領域Ａ４に属するか否か（以下、条件４）を判定する。そして、条件４が満たされる場合には、５０ｇ以上、６０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４４で条件４が満たされない場合、すなわち、第１指標算出工程で算出された（Ｐ₁__AVE，Ｐ₂__AVE）で表される点が、図１４に示すＰ₁__AVE−Ｐ₂__AVE平面において直線Ｌ４よりも下側にあり、図１４の領域Ａ５に属する場合には、５５ｇ以下が最適シャフト重量帯であると判定する。

以上の最適シャフト重量決定工程は、以下の知見に基づく。すなわち、本発明者らは、多数のゴルファーにテストクラブをスイングさせ、このときの腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEを算出した。また、同じゴルファーに、様々なシャフト重量のゴルフクラブをスイングさせ、最大の飛距離を与えるシャフト重量を算出し、これを最適シャフト重量とした。

上記の実験からは、最適シャフト重量は、同じゴルファーによるテストクラブのスイング時の腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEが共に大きいほど、大きくなる傾向にあることが確認された。その結果、本発明者らは、腕出力パワーＰ₁__AVE−クラブ入力パワーＰ₂__AVE空間を、図１４に示すように領域分割することで、最適シャフト重量帯に対応する領域Ａ１〜Ａ５を定義可能であることを発見した。すなわち、上述のステップＳ４１〜Ｓ４４は、腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEを示す点が、Ｐ₁__AVE−Ｐ₂__AVE空間におけるどの領域にプロットされるかに応じて、最適シャフト重量帯を判定するステップとなっている。ステップＳ４１〜Ｓ４４での判定に使用される閾値、言い換えると、図１４に示す分割領域Ａ１〜Ａ５の境界線Ｌ１〜Ｌ４を特定する情報等は、対応関係データ２８として記憶部２３内に格納されている。すなわち、対応関係データ２８とは、Ｐ₁__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だけでなく、最適フレックスも決定される。具体的には、決定部２４Ｅは、最適ランク値Ｋ_S1〜Ｋ_S4に基づいて、最適フレックスを決定する。フレックスとは、シャフト４０全体での硬さ（曲げ剛性）を評価する指標である。従って、シャフト４０の複数の位置におけるゴルファー７に適した曲げ剛性を表す最適ランク値Ｋ_S1〜Ｋ_S4が分かれば、これらの値に基づいて最適フレックスを算出することができる。例えば、特定の位置での最適ランク値を最適フレックスとすることもできるし、複数の位置での最適ランク値の平均値を最適フレックスをとすることもできる。フレックスには、例えば、「ＳＲ」「Ｓ」「Ｘ」などの様々な種類が存在するが、最適フレックスは、このうちの１つに特定される必要はなく、「ＳＲ」か「Ｓ」や、「Ｓ」か「Ｘ」等、幅を持つように設定することができる。

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

＜２−９．推奨シャフト選択工程＞
以上の工程Ｓ１〜Ｓ８により、最適シャフト重量帯及び最適剛性指標（最適ＥＩ値Ｊ_S1〜Ｊ_S4、最適ランク値Ｋ_S1〜Ｋ_S4及び最適フレックス）が決定されると、選択部２４Ｆは、推奨シャフト選択工程（Ｓ９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、ゴルファー７に適したシャフト４０(以下、推奨シャフト)が特定される。また、本実施形態では、推奨シャフトとして、最適シャフト重量帯を基準としたシャフトである第１推奨シャフトと、ゴルファー７の普段使用しているゴルフクラブの重量（以下、マイクラブ重量）を基準とした推奨シャフトである第２推奨シャフトとが決定される。

第１推奨シャフトの決定方法について説明する。まず、選択部２４Ｆは、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報には、製造メーカー、型番、４つの位置Ｈ１〜Ｈ４におけるＥＩ値Ｊ₁〜Ｊ₄及びランク値Ｋ₁〜Ｋ₄（ＩＦＣ）、シャフト４０の重量、フレックス、トルク、調子等が含まれる。そして、選択部２４Ｆは、これらの情報を参照して、シャフトＤＢ２９内に登録されている全てのシャフト４０の中から、第１絞り込みシャフトを特定する。第１絞り込みシャフトとは、シャフト４０の重量が最適シャフト重量帯に属し、かつ、フレックスが最適フレックスに一致するシャフト４０である。なお、第１絞り込みシャフトは、通常多数本存在する。

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

次に、第２推奨シャフトの決定方法について説明する。まず、選択部２４Ｆは、推奨ゴルフクラブに用いられるべきヘッド４１(以下、推奨ヘッド)の種類を決定する。推奨ヘッドの種類の決定は、本明細書では説明されないフィッティング処理により行うこともできるし、表示部２１及び入力部２２を介してユーザに質問する等してお好みのヘッド４１を選択させることにより行うこともできる。また、選択部２４Ｆは、表示部２１及び入力部２２を介してユーザに質問する等して、マイクラブ重量を特定する。

続いて、選択部２４Ｆは、推奨ヘッドのスペックを示す情報をヘッドＤＢ２７内から読み出すととともに、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。ヘッドＤＢ２７内に登録されているヘッド４１のスペックを示す情報には、製造メーカー、型番及び重量等が含まれる。そして、選択部２４Ｆは、これらの情報を参照して、シャフトＤＢ２９内に登録されている全てのシャフト４０の中から、第２絞り込みシャフトを特定する。第２絞り込みシャフトとは、最適フレックスに一致するフレックスを有し、推奨ヘッドと組み合わせたときの重量がマイクラブ重量±所定値（グリップの重量及びフェラルの重量等に誤差範囲を加えた値）の範囲内となる重量を有するシャフト４０である。なお、第２絞り込みシャフトは、通常多数本存在する。

続いて、選択部２４Ｆは、各第２絞り込みシャフトについて、当該シャフトのランク値Ｋ₁〜Ｋ₄と、最適剛性決定工程で決定された最適ランク値Ｋ_S1〜Ｋ_S4との一致度を算出し、第２推奨シャフトとして、一致度の高いシャフトを特定する。一致度は、数１３の式に従って算出することができる。なお、第１推奨シャフトとしては、１本のみ特定されてもよいし、複数本が特定されてもよい。第２推奨シャフトについても、同様である。

＜２−１０．ヘッド挙動決定工程＞
続いて、決定部２４Ｅは、計測工程で得られた画像データに基づいて、ヘッド４１の挙動を決定するヘッド挙動決定工程（Ｓ１０）を実行する。本実施形態では、ヘッド４１の挙動として、インパクト直前のヘッド速度Ｖ_hi、インパクト時のフェース角、ヘッド軌道の角度、打点を算出する。なお、ここでいうヘッド軌道とは、アウトサイドイン、インサイドアウト、インサイドイン等の軌道の種類を意味し、ヘッド軌道の角度として、定量的な角度が算出される。打点は、ヘッド４１のフェース面４１ａ上におけるボールの打点である。

インパクト直前のヘッド速度Ｖ_hiは、以下のように算出される。すなわち、投光器５５Ａ，５５Ｂと受光器５６Ａ，５６Ｂとの間隔は既知である。従って、タイミングｔ１，ｔ２が分かれば、インパクト直前のヘッド速度Ｖ_hiを算出することができるため、決定部２４Ｅは、タイミングｔ１，ｔ２の情報に基づいて、ヘッド速度Ｖ_hiを算出する。

また、インパクト時のフェース角及びヘッド軌道は、画像データを画像処理することにより決定される。図１８は、ヘッド４１をフェース面４１ａ側から視た図である。同図に示すように、本実施形態に係るヘッド４１のクラウン部４１ｂには、フェース面４１ａに沿うように帯状のマーカーＭ１が貼付されている。マーカーＭ１は、ストロボ５３，５４からの光を効率的に反射する素材で形成されている。従って、カメラ５１により撮影された画像上においては、マーカーＭ１の領域、言い換えると、ヘッド４１の平面視においてフェース面４１ａに沿った帯状の領域が鮮明に写り込むことになる。決定部２４Ｅは、記憶部２５内に記憶されているインパクト直前の２つのタイミングにおける画像（ストロボ５３，５４の発光のタイミングでの２枚の画像）から、マーカーＭ１の像を抽出する。そして、これらのマーカーＭ１の像に基づいて、インパクト時のフェース角を推定する。また、これらのマーカーＭ１の像の重心を結ぶベクトルと、飛球線方向のベクトルとの為す角度を算出し、この角度からアウトサイドイン等の軌道の種類を特定するとともに、この角度をヘッド軌道の角度として決定する。

また、決定部２４Ｅは、画像データを画像処理にすることにより、打点を決定する。より具体的には、インパクト時又はその付近の画像において、ボール及びフェース面４１ａの画像を抽出し、これらの位置関係から打点を算出する。本実施形態では、フェース面４１ａの幾何学中心であるフェースセンターから、打点がトゥ−ヒール方向にどれ位ずれているかが定量的に算出される。

＜２−１１．結果表示工程＞
次に、以上の解析の結果を示すＧＵＩ画面を作成し、これを表示する結果表示工程（Ｓ１１）について説明する。本実施形態では、ＧＵＩ画面として、図１９〜図２４に示すような画面Ｗ１〜Ｗ５が表示部２１上に表示される。

表示制御部２４Ｇは、入力部２２を介したユーザの操作を受けて、画面Ｗ１〜Ｗ５を適宜切り替えながら、表示部２１上に表示させる。本実施形態では、画面Ｗ１〜Ｗ４は、図１９、図２０、図２２及び図２４に示すようにタブ形式で表示され、現在選択されているタブに対応する画面が表示部２１上に表示される。

図１９に示す画面Ｗ１は、「スイングの傾向」というタブに対応する画面であり、文字どおり、主としてゴルファー７のスイングの傾向を表示する。まず、画面Ｗ１上の左半分の領域Ｃ１内では、「コックの動きの速さ度合い」、「手の甲の動き」、「手のローテーション」及び「コックの使い方」が、スライダー形式で定量的に表示される。従って、画面Ｗ１を見たユーザは、スライダーの位置といった図形的な観点から直感的に、各種指標の大小を正しく評価することができる。

「コックの動きの速さ度合い」としては、トップ付近、ダウンスイング前半、ダウンスイング後半及びダウンスイング全体の４つの期間における値が表示され、スライダーが左に配置されるほどこの度合いが小さく、右に配置されるほど大きいことを意味している。「コックの動きの速さ度合い」は、対応する期間の角速度ω_yに基づいて決定部２４Ｅにより決定される。「手の甲の動き」としては、トップ付近の値が表示され、スライダーが左に配置されるほど手の甲が開いており、右に配置されるほど閉じていることを意味している。「手の甲の動き」は、トップ付近の角速度ω_xに基づいて決定部２４Ｅにより決定される。「手のローテーション」としては、ダウンスイング全体及びインパクト直前の２つの期間における値が表示され、スライダーが左に配置されるほど手のローテーションが少なく、右に配置されるほど多いことを意味している。「手のローテーション」は、対応する期間における角速度ω_zに基づいて決定部２４Ｅにより決定される。

「コックの使い方」としては、「コックを溜める強さ」及び「コックを解放する強さ」が表示され、スライダーが左に配置されるほどコックの使い方が弱く、右に配置されるほど強いことを意味している。「コックを溜める強さ」は、上述したゴルファー７の腕出力パワーＰ₁__AVEに応じて決定部２４Ｅにより決定される。「コックを解放する強さ」は、上述したゴルファー７のクラブ入力パワーＰ₂__AVEに応じて決定部２４Ｅにより決定される。

画面Ｗ１上の右上には、ゴルファー７のスイングタイプが表示される。ゴルファー７のスイングタイプとは、最適ＥＩ値Ｊ_S1〜Ｊ_S4（又は、最適ランク値Ｋ_S1〜Ｋ_S4）、すなわち、ゴルファー７に適したシャフト４０のＩＦＣ（以下、理想ＩＦＣという）に基づいて、決定部２４Ｅにより決定される。より具体的には、スイングタイプは、理想ＩＦＣの４つのコードが、シャフト４０の手元から手先にかけて右上がりであるか、右下がりであるか、中央が突出しているか、一定であるか等を表す。ここでは、横軸をシャフト４０の長さ方向の位置とし、縦軸をシャフト４０の硬さをとしたグラフ領域内において、理想ＩＦＣをグラフ状に表す画像Ｃ２が表示される。ただし、スイングタイプは、ゴルファー７のタイプを知るための指標であるため、シャフトの位置とＩＦＣのコードとが厳密に描画されるのではなく、模式化されたグラフが表示される。

また、画面Ｗ１には、「平均ヘッドスピード」として、上述したヘッドスピードＶ_h（平均値）が表示され、「推奨シャフト硬度」として、最適フレックスが表示され、「推奨重量帯」として、最適シャフト重量帯が表示される。ヘッドスピードＶ_h（平均値）は、数値で表示される。一方、「推奨シャフト硬度」及び「推奨重量帯」の欄には、それぞれ、シャフトＤＢ２９内に登録されているシャフト４０のフレックス及びシャフト重量の範囲がスケールの形態で表示されており、当該スケールにおいてゴルファー７に対応する部分が強調表示される。これにより、ユーザは、スケールを参照して、ゴルファーに対応する最適フレックス及び最適シャフト重量帯を直感的に理解することができる。

図２０に示す画面Ｗ２は、「理想のＩＦＣ」というタブに対応する画面であり、主として理想ＩＦＣ及び推奨シャフトを表示する。画面Ｗ２には、理想ＩＦＣのコードが数値で表示され、「推奨シャフトタイプ」として、上述したスイングタイプが文字で表示される。また、画面Ｗ２には、「推奨シャフト１」及び「推奨シャフト２」として、上述した上位２本の第２推奨シャフトを特定する情報を表示する欄Ｃ３，Ｃ３が表示される。また、「推奨シャフトＷ」として、上述した１本の第１推奨シャフトを特定する情報を表示する欄Ｃ４が表示される。より具体的には、欄Ｃ３，Ｃ３，Ｃ４内には、製造メーカー、型番、質量、フレックス等が表示される。なお、これらの情報は、シャフトＤＢ２９内を参照して特定される。

また、画面Ｗ２には、以上の理想ＩＦＣに加え、上述した１本の第１推奨シャフト及び２本の第２推奨シャフトのＩＦＣをグラフ状に表示するためのグラフ領域Ｃ５が表示される。これらのグラフは、ＩＦＣのコードに対応する位置に配置され、当該コードが数値で記述されたプロット点と、これらのプロット点を連結する直線部分とから構成される。グラフ領域Ｃ５内では、理想ＩＦＣのグラフは、３本の推奨シャフトに対応するＩＦＣのグラフとは異なる態様で、好ましくは、ユーザの注意をより引きつけるような態様（図２０では、グラフの線が太く、プロット点が大きい）で表示される。また、どのグラフがどの欄Ｃ３，Ｃ３，Ｃ４により特定されるシャフトに対応するかが一目で分かるように、欄Ｃ３，Ｃ３，Ｃ４の装飾と、これらに対応するグラフの装飾とに統一感を持たせることが好ましい。例えば、欄Ｃ３，Ｃ３，Ｃ４内の文字の色と、グラフの色とを一致させることができる。

欄Ｃ３，Ｃ３，Ｃ４は、排他的に選択可能に構成されている。そして、領域Ｃ６には、現在選択されている欄に対応する推奨シャフトのデザイン画像が表示される。これにより、現在選択されている推奨シャフトが、どのシャフトであるのかを直感的に理解することができる。また、グラフ領域Ｃ５内では、現在選択されている推奨シャフトのＩＦＣのグラフは、その他の推奨シャフトのＩＦＣのグラフと異なる態様で、好ましくは、ユーザの注意をより引きつけるような態様（例えば、前者の方が色が濃く、後者の方が色が薄い）で表示される。

ところで、図２０の例では、画面Ｗ２の左側に、「コックの動きの速さ度合い」、「手の甲の動き」及び「手のローテーション」が、図１９と同じ形式で表示されている。そして、画面Ｗ１，Ｗ２上の「コックの動きの速さ度合い」を示す欄には、アイコンＣ７が表示されており、これが選択されると、図２１に示す画面Ｗ５がポップアップにより表示される。この画面Ｗ５上には、コックの動きの速さ度合いの意味するところを説明する画像が表示される。より具体的には、この画面Ｗ５上では、コックの動きとは手のどのような方向の動きであるのかが、手の絵、シャフトの絵、及び動きの方向を表す矢印等を用いて図形的に説明される。

次に、図２２に示す画面Ｗ３は、は、「推奨シャフト重量」というタブに対応する画面であり、文字どおり、主としてゴルファー７に推奨されるべきシャフト４０の重量帯を表示する。より具体的には、画面Ｗ３には、画面Ｗ１と同様に、ゴルファー７の最適シャフト重量帯を表示する「推奨重量帯」の欄Ｃ８が表示される。表示形態も画面Ｗ１と同様であり、欄Ｃ８には、シャフトＤＢ２９内に登録されているシャフト重量の範囲がスケールの形態で表示されており、当該スケールにおいて最適シャフト重量帯の部分が強調表示されている。これにより、ユーザは、スケールを参照して、最適シャフト重量帯を直感的に理解することができる。

また、画面Ｗ３には、最適シャフト重量帯についてのユーザの納得性を高めるべく、その算出の根拠を説明するための様々な欄Ｃ９〜Ｃ１０が表示される。欄Ｃ１０には、最適シャフト重量帯を決定するための指標であるゴルファー７の腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEが表示される。より具体的には、これらの指標Ｐ₁__AVE，Ｐ₂__AVEが、それぞれ「コックを溜める強さ」及び「コックを解放する強さ」として、画面Ｗ１と同様にスライダー形式で表示される。これにより、ユーザは、最適シャフト重量帯を決定する根拠となる、ゴルファー７固有のコックを溜める強さ及びコックを解放する強さの大小を、直感的に正しく理解することができる。

欄Ｃ９には、図１４に示したような、腕出力パワーＰ₁__AVEを横軸とし、クラブ入力パワーＰ₂__AVEを縦軸とするＰ₁__AVE−Ｐ₂__AVE平面であるプロット領域が表示され、このプロット領域内には、ゴルファー７の腕出力パワーＰ₁__AVE及びクラブ入力パワーＰ₂__AVEに対応する点がプロットされ、表示される。欄Ｃ９では、プロット領域の横軸及び縦軸が、それぞれ「コックを溜める強さ」「コックを解放する強さ」を表すことをが、軸に添えられた文字により明確に示されている。また、このプロット領域は、ゴルファー７に推奨すべきシャフト重量帯にそれぞれ対応する、上述した領域Ａ１〜Ａ５に分割されている。これらの領域は、プロット領域内において、互いに異なるシャフト重量帯に対応する領域であることが分かるように異なる態様で表示される。例えば、領域間に境界線が引かれたり、異なる領域に異なる色や模様が付される。また、各領域には、対応するシャフト重量帯に属するシャフト４０の仕様を表す記号が表示される。図２２の例では、「８Ｘ・Ｓ」「７Ｘ・Ｓ」等の記号が表示されているが、これらの８や７等の数値は、シャフト重量帯を意味しており、Ｘ、Ｓ等のアルファベットは、フレックスを示している。

ところで、コックを溜める強さ及びコックを解放する強さという概念は、全てのユーザが知っているものではない。そのため、欄Ｃ１１には、「スイング中に発生する力、腕とクラブの関係」を示す画像が表示される。本実施形態では、コックを溜める強さ及びコックを解放する強さを説明するための動画が表示される。この動画は、ゴルファー７の腕を模式的に示す図形（以下、腕図形）と、ゴルファー７により把持されたゴルフクラブ４を模式的に示す図形（以下、クラブ図形）とが、スイング動作に伴って二重振り子のように移動する様子を示す動画である。

今、スイング動作中のタイミングを、時系列に沿って序数詞ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）で表すこととする。このとき、欄Ｃ１１内で表示される動画は、時系列に沿ってコマ送りの形式で再生されるＮ枚の画像Ｉ（ｉ）を含む。画像Ｉ（ｉ）上には、第ｉ番目のタイミングにおいて腕が存在する位置に腕図形Ｇ₃（ｉ）が表示され、第ｉ番目のタイミングにおいてゴルフクラブ４が存在する位置にクラブ図形Ｇ₄（ｉ）が表示される。図２３Ａ〜図２３Ｆは、画像Ｉ（ｉ）の例を示しており、図２３Ａ〜図２３Ｆの画像は、この順に時系列に沿っている。

また、画像Ｉ（ｉ）上には、腕図形Ｇ₃（ｉ）及びクラブ図形Ｇ₄（ｉ）に加えて、第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ₁（ｉ）と、同じく第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ₂（ｉ）とが表示される。同じタイミングに対応するこれらの図形Ｇ₁（ｉ）〜Ｇ₄（ｉ）は、１枚の画像Ｉ（ｉ）上において同時に表示される。第ｉ番目のタイミングにおけるコックを溜める強さとは、第ｉ番目のタイミングでの腕出力パワーＰ₁であり、第ｉ番目のタイミングにおけるコックを解放する強さとは、第ｉ番目のタイミングでのクラブ入力パワーＰ₂である。画像Ｉ（ｉ）上の第１図形Ｇ₁（ｉ）は、第ｉ番目のタイミングでの腕出力パワーＰ₁の大きさに応じた形態をとり、第２図形Ｇ₂（ｉ）は、第ｉ番目のタイミングでのクラブ入力パワーＰ₂の大きさに応じた形態をとる。より具体的には、本実施形態では、腕出力パワーＰ₁が大きいタイミングである程、第１図形Ｇ₁（ｉ）が大きく表示され、クラブ入力パワーＰ₂が大きいタイミングである程、第２図形Ｇ₂（ｉ）が大きく表示される。なお、腕出力パワーＰ₁及びクラブ入力パワーＰ₂の大きさは、第１図形Ｇ₁（ｉ）及び第２図形Ｇ₂（ｉ）の色等を変化させることにより表現されてもよい。また、第１図形Ｇ₁（ｉ）及び第２図形Ｇ₂（ｉ）の形態を決定する腕出力パワーＰ₁及びクラブ入力パワーＰ₂は、計測データに基づくゴルファー７固有の値である必要はなく、一般的なスイング動作時の腕出力パワーＰ₁及びクラブ入力パワーＰ₂とすることができる。この場合、欄Ｃ１１内では、ゴルファー７に関わらず同じ動画が再生される。

ところで、観念上、コックを溜める力は、腕に関連付けられ、コックを解放する力は、ゴルフクラブ４に関連付けられる。従って、本実施形態では、第１図形Ｇ₁（ｉ）は、腕図形Ｇ₃（ｉ）に重ねて配置され、第２図形Ｇ₂（ｉ）は、クラブ図形Ｇ₄（ｉ）に重ねて配置される。ただし、第１図形Ｇ₁（ｉ）及び第２図形Ｇ₂（ｉ）の配置はこれに限られず、別の好ましい例としては、第１図形Ｇ₁（ｉ）を腕図形Ｇ₃（ｉ）の近傍に配置し、第２図形Ｇ₂（ｉ）をクラブ図形Ｇ₄（ｉ）の近傍に配置することもできる。また、第１図形Ｇ₁（ｉ）と腕図形Ｇ₃（ｉ）とを概ね一体の図形として表現することもできる。この場合、例えば、腕を模した図形の形態を、腕出力パワーＰ₁の大きさに応じて変化させるように構成することができる。第２図形Ｇ₂（ｉ）及びクラブ図形Ｇ₄（ｉ）についても同様であり、この場合、例えば、ゴルフクラブ４を模した図形の形態を、クラブ入力パワーＰ₂の大きさに応じて変化させるように構成することができる。また、第１図形Ｇ₁（ｉ）を腕図形の肩に対応する部分又はその近傍に配置し、かつ、第２図形Ｇ₂（ｉ）を腕図形の手首に対応する部分又はその近傍に配置することもできる。以上のような動画を見たユーザは、ゴルフスイング中のコックを溜める動き及びコックを解放する動きの意味するところを直感的に理解することができる。

図２４に示す画面Ｗ４は、「ヘッド軌道」というタブに対応する画面であり、主としてヘッド４１の挙動を表示する。画面Ｗ４には、ヘッド挙動決定工程で決定されたインパクト直前のヘッド速度Ｖ_hi、ヘッド軌道の角度、インパクト時のフェース角及び打点が定量的に数値で表示される。また、欄Ｃ１２には、インパクト付近の２つのタイミングにおけるヘッドを表す図形であるヘッド図形Ｇ₅，Ｇ₆が並べて同時に表示される。ヘッド図形Ｇ₅は、インパクト時のヘッドを表しており、ヘッド図形Ｇ₆は、その直前のヘッドを表している。また、欄Ｃ１２内のヘッド図形Ｇ₅は、上述したヘッド軌道の角度に応じて傾けて表示される。これにより、ユーザは、アウトサイドインかインサイドアウトかといったゴルファー７のヘッド軌道を直感的に正確に理解することができる。なお、ヘッド図形Ｇ₅の傾きは、インパクト時のフェース角度に応じて設定することもできる。

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

＜３−１＞
上記実施形態では、コックを溜める強さである腕出力パワーＰ₁，Ｐ₁__AVE、及び、コックを解放する強さであるクラブ入力パワーＰ₂，Ｐ₂__AVEが、慣性センサユニットによる計測データに基づいて計算されたが、これらの指標の少なくとも一方を、カメラユニットによる計測データに基づいて算出することもできる。

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

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

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

また、最適剛性指標として、ゴルファー７に適したＩＦＣやフレックスに代えて又は加えて、ゴルファー７に適した調子やトルクを決定するようにしてもよい。なお、トルクは、シャフト４０全体でのねじれ剛性を表す指標である。

＜３−３＞
上記実施形態では、ｉ＝１，２，・・・，Ｎについての画像Ｉ（ｉ）が、動画の形式で表示されたが、これらＮ枚の画像Ｉ（ｉ）を静止画の態様で表示することもできる。すなわち、図２３Ａ〜図２３Ｆのような画像を、同じ画面上に同時に表示することができる。このとき、時系列に沿って、画像Ｉ（ｉ）を配列することが好ましい。また、図２２に示すような１枚の画像のみを表示することもできる。図２２に示す画像では、様々なタイミングでの図形Ｇ₁（ｉ）〜Ｇ₄（ｉ）がストロボ画像のように１枚の画像上に同時に表示されている。従って、このような１枚の画像のみであっても、スイング動作中のコックを溜める強さ及びコックを解放する強さを直感的に理解することができる。

１ 慣性センサユニット
２ フィッティング装置（ＧＵＩ表示装置）
３ プログラム（ＧＵＩ表示プログラム）
４ ゴルフクラブ
５ カメラシステム
７ ゴルファー
２４Ａ 取得部
２４Ｂ グリップ挙動導出部
２４Ｃ 肩挙動導出部
２４Ｄ 算出部
２４Ｅ 決定部
２４Ｆ 選択部
２４Ｇ 表示制御部
４０ シャフト
４１ ヘッド
４２ グリップ
Ｗ１〜Ｗ５ ＧＵＩ画面
Ｇ₁（ｉ） 第１図形
Ｇ₂（ｉ） 第２図形
Ｇ₃（ｉ） 腕図形
Ｇ₄（ｉ） クラブ図形
Ｇ₅ ヘッド図形（第１ヘッド図形）
Ｇ₆ ヘッド図形（第２ヘッド図形）

Claims (7)

1. ＧＵＩ画面を表示するＧＵＩ表示装置であって、
   ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さとして、第ｉ番目のタイミングにおけるゴルファーの腕が出力するパワー（仕事率）である腕出力パワーの大きさを表す第１図形Ｇ₁（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さとして、第ｉ番目のタイミングにおける前記ゴルフクラブに入力されるパワー（仕事率）であるクラブ入力パワーの大きさを表す第２図形Ｇ₂（ｉ）と、第ｉ番目のタイミングにおける前記ゴルファーの腕を表す図形である腕図形Ｇ ₃ （ｉ）と、第ｉ番目のタイミングにおける前記ゴルフクラブを表す図形であるクラブ図形Ｇ ₄ （ｉ）とを同時に表示する表示制御部
   を備え、
   前記腕出力パワーは、第ｉ番目のタイミングにおいて前記腕で発揮される仕事率のうち、スイング平面内で前記ゴルファーの肩を中心として前記ゴルフクラブを回転させる動作での仕事率であり、前記クラブ入力パワーは、第ｉ番目のタイミングにおいて前記腕で発揮される仕事率のうち、前記スイング平面内で前記ゴルフクラブのグリップを中心として前記ゴルフクラブを回転させる動作での仕事率であり、
   前記表示制御部は、前記腕出力パワーの大きさに応じて、前記第１図形Ｇ₁（ｉ）の形態を変化させるとともに、前記クラブ入力パワーの大きさに応じて、前記第２図形Ｇ₂（ｉ）の形態を変化させ、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記腕図形Ｇ ₃ （ｉ）に重ねて又はその近傍に前記第１図形Ｇ ₁ （ｉ）を表示し、前記クラブ図形Ｇ ₄ （ｉ）に重ねて又はその近傍に前記第２図形Ｇ ₂ （ｉ）を表示する、
   ＧＵＩ表示装置。
2. 前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ₁（ｉ）及び前記第２図形Ｇ₂（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って動画の態様で表示する、
   請求項１に記載のＧＵＩ表示装置。
3. 前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ₁（ｉ）及び前記第２図形Ｇ₂（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って配列する態様で表示する、
   請求項１に記載のＧＵＩ表示装置。
4. 前記ゴルファーによる前記ゴルフクラブのスイング動作を計測した計測データを取得する取得部と、
   前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する算出部と
   をさらに備え、
   前記表示制御部は、前記ＧＵＩ画面上において、前記第１指標及び前記第２指標を表示する、
   請求項１から３のいずれかに記載のＧＵＩ表示装置。
5. 前記第１指標及び前記第２指標に応じて、前記ゴルファーに適したシャフト重量である最適シャフト重量を決定する決定部
   をさらに備え、
   前記表示制御部は、前記ＧＵＩ画面上において、前記最適シャフト重量を表示する、
   請求項４に記載のＧＵＩ表示装置。
6. ＧＵＩ画面を表示するＧＵＩ表示プログラムであって、
   ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さとして、第ｉ番目のタイミングにおけるゴルファーの腕が出力するパワー（仕事率）である腕出力パワーの大きさを表す第１図形Ｇ₁（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さとして、第ｉ番目のタイミングにおける前記ゴルフクラブに入力されるパワー（仕事率）であるクラブ入力パワーの大きさを表す第２図形Ｇ₂（ｉ）と、第ｉ番目のタイミングにおける前記ゴルファーの腕を表す図形である腕図形Ｇ ₃ （ｉ）と、第ｉ番目のタイミングにおける前記ゴルフクラブを表す図形であるクラブ図形Ｇ ₄ （ｉ）とを同時に表示するステップ
   をコンピュータに実行させ、
   前記腕出力パワーは、第ｉ番目のタイミングにおいて前記腕で発揮される仕事率のうち、スイング平面内で前記ゴルファーの肩を中心として前記ゴルフクラブを回転させる動作での仕事率であり、前記クラブ入力パワーは、第ｉ番目のタイミングにおいて前記腕で発揮される仕事率のうち、前記スイング平面内で前記ゴルフクラブのグリップを中心として前記ゴルフクラブを回転させる動作での仕事率であり、
   前記表示するステップは、前記腕出力パワーの大きさに応じて、前記第１図形Ｇ₁（ｉ）の形態を変化させるとともに、前記クラブ入力パワーの大きさに応じて、前記第２図形Ｇ₂（ｉ）の形態を変化させるステップと、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記腕図形Ｇ ₃ （ｉ）に重ねて又はその近傍に前記第１図形Ｇ ₁ （ｉ）を表示し、前記クラブ図形Ｇ ₄ （ｉ）に重ねて又はその近傍に前記第２図形Ｇ ₂ （ｉ）を表示するステップとを含む、
   ＧＵＩ表示プログラム。
7. ＧＵＩ画面を表示するＧＵＩ表示装置であって、
   ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する取得部と、
   前記計測データに基づいて、コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する算出部と、
   前記ＧＵＩ画面上において、前記コックを溜める強さを第１軸とし、前記コックを解放する強さを第２軸とするプロット平面内にプロットされる点として、前記第１指標及び前記第２指標を表示する表示制御部と
   を備え、
   前記算出部は、前記第１指標として、前記ゴルファーの腕が出力するパワー（仕事率）である腕出力パワーを算出するとともに、前記第２指標として、前記ゴルフクラブに入力されるパワー（仕事率）であるクラブ入力パワーを算出し、
   前記腕出力パワーは、前記腕で発揮される仕事率のうち、スイング平面内で前記ゴルファーの肩を中心として前記ゴルフクラブを回転させる動作での仕事率であり、前記クラブ入力パワーは、前記腕で発揮される仕事率のうち、前記スイング平面内で前記ゴルフクラブのグリップを中心として前記ゴルフクラブを回転させる動作での仕事率であり、
   前記プロット平面は、前記ゴルファーに推奨すべきシャフト重量帯に分割されている、ＧＵＩ表示装置。
   
   

Images