(57)【要約】【課題】ゴルフスイング中のコックを溜める動き及びコックを解放する動きを視覚的に分かり易く説明することができるＧＵＩ表示装置を提供する。【解決手段】ＧＵＩ画面を表示するＧＵＩ表示装置が提供される。前記ＧＵＩ表示装置は、表示制御部を備える。前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とを同時に表示する。【選択図】図２２
2017217423.tif 000002
【特許請求の範囲】
【請求項１】
ＧＵＩ画面を表示するＧＵＩ表示装置であって、
ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とを同時に表示する表示制御部
を備える、
ＧＵＩ表示装置。
【請求項２】
前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って動画の態様で表示する、
請求項１に記載のＧＵＩ表示装置。
【請求項３】
前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って配列する態様で表示する、
請求項１に記載のＧＵＩ表示装置。
【請求項４】
前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）に加え、第ｉ番目のタイミングにおける前記ゴルファーの腕を表す図形である腕図形Ｇ3（ｉ）と、第ｉ番目のタイミングにおける前記ゴルフクラブを表す図形であるクラブ図形Ｇ4（ｉ）とを同時に表示する、
請求項１から３のいずれかに記載のＧＵＩ表示装置。
【請求項５】
前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記腕図形Ｇ3（ｉ）に重ねて又はその近傍に前記第１図形Ｇ1（ｉ）を表示し、前記クラブ図形Ｇ4（ｉ）に重ねて又はその近傍に前記第２図形Ｇ2（ｉ）を表示する、
請求項４に記載のＧＵＩ表示装置。
【請求項６】
ゴルファーによる前記ゴルフクラブのスイング動作を計測した計測データを取得する取得部と、
前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する算出部と
をさらに備え、
前記表示制御部は、前記ＧＵＩ画面上において、前記第１指標及び前記第２指標を表示する、
請求項１から５のいずれかに記載のＧＵＩ表示装置。
【請求項７】
前記第１指標及び前記第２指標に応じて、前記ゴルファーに適したシャフト重量である最適シャフト重量を決定する決定部
をさらに備え、
前記表示制御部は、前記ＧＵＩ画面上において、前記最適シャフト重量を表示する、
請求項６に記載のＧＵＩ表示装置。
【請求項８】
ＧＵＩ画面を表示するＧＵＩ表示プログラムであって、
ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とを同時に表示するステップ
をコンピュータに実行させる、
ＧＵＩ表示プログラム。
【請求項９】
ＧＵＩ画面を表示するＧＵＩ表示装置であって、
ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する取得部と、
前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する算出部と、
前記ＧＵＩ画面上において、前記コックを溜める強さを第１軸とし、前記コックを解放する強さを第２軸とするプロット平面内にプロットされる点として、前記第１指標及び前記第２指標を表示する表示制御部と
を備え、
前記プロット平面は、前記ゴルファーに推奨すべきシャフト重量帯に分割されている、
ＧＵＩ表示装置。
【請求項１０】
ＧＵＩ画面を表示するＧＵＩ表示装置であって、
ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する取得部と、
前記計測データに基づいて、インパクト時の前記ゴルフクラブに含まれるヘッドの傾きを決定する決定部と、
前記ＧＵＩ画面上において、インパクト時の前記ヘッドを表す図形である第１ヘッド図形と、インパクトの直前のタイミングにおける前記ヘッドを表す図形である第２ヘッド図形とを同時に表示する表示制御部と
を備え、
前記第１ヘッド図形は、前記ヘッドの傾きに応じて傾くように表示される、
ＧＵＩ表示装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴルフクラブのスイング動作に関するＧＵＩ（グラフィカルユーザーインターフェース）画面を表示するＧＵＩ表示装置に関する。
【背景技術】
【０００２】
従来より、ゴルファーによるゴルフクラブのスイング動作を計測装置により計測し、計測データに基づいてスイング動作を解析する装置が公知である（特許文献１等）。このような解析の場面では、解析のアルゴリズム自体も勿論重要であるが、解析の結果をユーザに分かり易く説明することも重要である。特許文献１には、ゴルフスイングの解析結果をグラフ等を用いて分かり易く視認できるようにしたＧＵＩ画面が開示されている。なお、このようなゴルフスイングの解析結果は、ゴルファーに適したゴルフクラブのフィッティングや、ゴルファーのフォームの改善等、様々な用途で利用される。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−１２１４１２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、ゴルフスイング中には、コックを溜める動きと、コックを解放する動きとが行われ、このような動きは、ゴルファーのスイング動作を特徴付ける指標となる。従って、本発明者らは、ゴルファー固有のこれらの動きの強さを解析し、ゴルファーに提示できるようにすることが重要であると考えた。しかしながら、コックを溜める動き及びコックを解放する動きが何を意味するのか、多くの一般のゴルファーが理解しているわけではない。
【０００５】
本発明は、ゴルフスイング中のコックを溜める動き及びコックを解放する動きを視覚的に分かり易く説明することができるＧＵＩ表示装置を提供することを１つの目的とする。
【課題を解決するための手段】
【０００６】
第１観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、表示制御部を備える。前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とを同時に表示する。
【０００７】
第２観点に係るＧＵＩ表示装置は、第１観点に係るＧＵＩ表示装置であって、前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って動画の態様で表示する。
【０００８】
第３観点に係るＧＵＩ表示装置は、第１観点に係るＧＵＩ表示装置であって、前記表示制御部は、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）を同時に表示する画像Ｉ（ｉ）を、ｉ＝１，２，・・・，Ｎの時系列に沿って配列する態様で表示する。
【０００９】
第４観点に係るＧＵＩ表示装置は、第１観点から第３観点のいずれかに係るＧＵＩ表示装置であって、前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記第１図形Ｇ1（ｉ）及び前記第２図形Ｇ2（ｉ）に加え、第ｉ番目のタイミングにおける前記ゴルファーの腕を表す図形である腕図形Ｇ3（ｉ）と、第ｉ番目のタイミングにおける前記ゴルフクラブを表す図形であるクラブ図形Ｇ4（ｉ）とを同時に表示する。
【００１０】
第５観点に係るＧＵＩ表示装置は、第４観点に係るＧＵＩ表示装置であって、前記表示制御部は、ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記腕図形Ｇ3（ｉ）に重ねて又はその近傍に前記第１図形Ｇ1（ｉ）を表示し、前記クラブ図形Ｇ4（ｉ）に重ねて又はその近傍に前記第２図形Ｇ2（ｉ）を表示する。
【００１１】
第６観点に係るＧＵＩ表示装置は、第１観点から第５観点のいずれかに係るＧＵＩ表示装置であって、取得部と、算出部とをさらに備える。前記取得部は、ゴルファーによる前記ゴルフクラブのスイング動作を計測した計測データを取得する。前記算出部は、前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する。前記表示制御部は、前記ＧＵＩ画面上において、前記第１指標及び前記第２指標を表示する。
【００１２】
第７観点に係るＧＵＩ表示装置は、第６観点に係るＧＵＩ表示装置であって、決定部をさらに備える。前記決定部は、前記第１指標及び前記第２指標に応じて、前記ゴルファーに適したシャフト重量である最適シャフト重量を決定する。前記表示制御部は、前記ＧＵＩ画面上において、前記最適シャフト重量を表示する。
【００１３】
第８観点に係るＧＵＩ表示装プログラムは、ＧＵＩ画面を表示するＧＵＩ表示プログラムであって、以下の（１）のステップをコンピュータに実行させる。
（１）ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）に対し、前記ＧＵＩ画面上において、前記ゴルフクラブのスイング動作中の第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とを同時に表示するステップ。
【００１４】
第９観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、取得部と、算出部と、表示制御部とを備える。前記取得部は、ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する。前記算出部は、前記計測データに基づいて、前記コックを溜める強さを表す第１指標と、前記コックを解放する強さを表す第２指標とを算出する。前記表示制御部は、前記ＧＵＩ画面上において、前記コックを溜める強さを第１軸とし、前記コックを解放する強さを第２軸とするプロット平面内にプロットされる点として、前記第１指標及び前記第２指標を表示する。前記プロット平面は、前記ゴルファーに推奨すべきシャフト重量帯に分割されている。
ＧＵＩ表示装置。
【００１５】
第１０観点に係るＧＵＩ表示装置は、ＧＵＩ画面を表示するＧＵＩ表示装置であって、取得部と、決定部と、表示制御部とを備える。前記取得部は、ゴルファーによるゴルフクラブのスイング動作を計測した計測データを取得する。前記決定部は、前記計測データに基づいて、インパクト時の前記ゴルフクラブに含まれるヘッドの傾きを決定する。前記表示制御部は、前記ＧＵＩ画面上において、インパクト時の前記ヘッドを表す図形である第１ヘッド図形と、インパクトの直前のタイミングにおける前記ヘッドを表す図形である第２ヘッド図形とを同時に表示する。前記第１ヘッド図形は、前記ヘッドの傾きに応じて傾くように表示される。
【発明の効果】
【００１６】
第１観点によれば、ゴルフスイング中の複数のタイミングにおける、コックを溜める強さを表す第１図形及びコックを解放する強さを表す第２図形が、ＧＵＩ画面上に表示される。また、同じタイミングにおける第１図形及び第２図形は、同時に表示される。従って、このようなＧＵＩ画面を視たユーザは、ゴルフスイング中のコックを溜める動き及びコックを解放する動きの意味するところを直感的に理解することができる。すなわち、このようなＧＵＩ画面により、ゴルフスイング中のコックを溜める動き及びコックを解放する動きを視覚的に分かり易く説明することができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の一実施形態に係るＧＵＩ表示装置（フィッティング装置）を備えるフィッティングシステムを示す図。
【図２】フィッティングシステムの機能ブロック図。
【図３】ゴルフクラブのグリップを基準とするｘｙｚ局所座標系を説明する図。
【図４】（Ａ）アドレス状態を示す図。（Ｂ）トップ状態を示す図。（Ｃ）インパクト状態を示す図。（Ｄ）フィニッシュ状態を示す図。
【図５】フィッティング処理の流れを示すフローチャート。
【図６】スイング平面を説明する図。
【図７】二重振り子モデルを概念的に説明する図。
【図８】肩挙動導出工程を概念的に説明する図。
【図９】肩挙動導出工程を概念的に説明する別の図。
【図１０】第１指標算出工程の流れを示すフローチャート。
【図１１】二重振り子モデルを概念的に説明する別の図。
【図１２】腕エネルギーを説明する図。
【図１３】最適シャフト重量決定工程の流れを示すフローチャート。
【図１４】複数の最適シャフト重量帯に対応する領域に分割された腕出力パワー−クラブ入力パワー平面を示す図。
【図１５】インターナショナル・フレックス・コード（ＩＦＣ）を説明する図。
【図１６】シャフトの曲げ剛性の測定方法を説明する図。
【図１７】スイング中のシャフトの曲げを説明する図。
【図１８】ヘッドをフェース面側から視た図。
【図１９】解析の結果を表示するＧＵＩ画面を示す図。
【図２０】解析の結果を表示するＧＵＩ画面を示す別の図。
【図２１】解析の結果を表示するＧＵＩ画面を示すさらに別の図。
【図２２】解析の結果を表示するＧＵＩ画面を示すさらに別の図。
【図２３Ａ】スイング中に発生する力を説明する画像を示す図。
【図２３Ｂ】続くタイミングでのスイング中に発生する力を説明する画像を示す図。
【図２３Ｃ】さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。
【図２３Ｄ】さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。
【図２３Ｅ】さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。
【図２３Ｆ】さらに続くタイミングでのスイング中に発生する力を説明する画像を示す図。
【図２４】解析の結果を表示するＧＵＩ画面を示すさらに別の図。
【発明を実施するための形態】
【００１８】
以下、図面を参照しつつ、本発明の一実施形態に係るＧＵＩ表示装置及びプログラムについて説明する。
【００１９】
＜１．フィッティングシステムの概略構成＞
図１及び図２に、本実施形態に係るＧＵＩ表示装置であるフィッティング装置２を備えるフィッティングシステム１００の全体構成を示す。フィッティング装置２は、ゴルファー７に適したゴルフクラブ４を選定するのを支援するための装置である。そのために、フィッティング装置２は、ゴルファー７によるゴルフクラブ４のスイング動作を計測した計測データを取得し、当該計測データに基づいて、スイング動作の解析を行う。解析の結果は、後述する画面Ｗ１〜Ｗ５（図１９〜図２４）等を介してユーザに提示される。なお、ここでいうユーザとは、ゴルファー７自身やそのインストラクター等の、解析の結果を必要とする者の総称である。本実施形態では、スイング動作の計測を行うセンサユニットは、慣性センサユニット１及びカメラシステム５であり、フィッティング装置２は、これらの慣性センサユニット１及びカメラシステム５とともに、フィッティングシステム１００を構成する。以下、慣性センサユニット１、カメラシステム５及びフィッティング装置２の構成について説明した後、フィッティング処理の流れについて説明する。
【００２０】
＜１−１．慣性センサユニットの構成＞
慣性センサユニット１は、図１及び図３に示すとおり、ゴルフクラブ４のグリップ４２におけるヘッド４１と反対側の端部に取り付けられており、グリップ４２の挙動を計測する。なお、ゴルフクラブ４は、一般的なゴルフクラブであり、シャフト４０と、シャフト４０の一端に設けられたヘッド４１と、シャフト４０の他端に設けられたグリップ４２とから構成される。本実施形態に係るシャフト４０は、カーボン製のシャフトである。慣性センサユニット１は、スイング動作の妨げとならないよう、小型且つ軽量に構成されている。慣性センサユニット１は、ゴルフクラブ４の外側に取り付けることができ、ゴルフクラブ４に対して着脱自在に構成することができる。
【００２１】
図２に示すように、本実施形態に係る慣性センサユニット１には、加速度センサ１１、角速度センサ１２及び地磁気センサ１３が搭載されている。また、慣性センサユニット１には、これらのセンサ１１〜１３から出力されるセンサデータ（計測データ）を外部のフィッティング装置２に送信するための通信装置１０も搭載されている。なお、本実施形態では、通信装置１０は、スイング動作の妨げにならないように無線式であるが、ケーブルを介して有線式にフィッティング装置２に接続するようにしてもよい。
【００２２】
加速度センサ１１、角速度センサ１２及び地磁気センサ１３はそれぞれ、グリップ４２を基準としたｘｙｚ局所座標系におけるグリップ加速度、グリップ角速度及びグリップ地磁気を計測する。より具体的には、加速度センサ１１は、ｘ軸、ｙ軸及びｚ軸方向のグリップ加速度ａ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に変換する第１変換工程
（Ｓ３）ＸＹＺ全体座標系でのグリップ４２の挙動を、スイング平面Ｐ（後述する）内でのグリップ４２の挙動へと変換する第２変換工程
（Ｓ４）スイング平面Ｐ内でのグリップ４２の挙動に基づいて、スイング平面Ｐ内でのゴルファー７の疑似的な肩の挙動を導出する肩挙動導出工程
（Ｓ５）スイング平面Ｐ内でのグリップ４２の挙動及び疑似的な肩の挙動に基づいて、後述する腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVE（以下、これらをまとめて、第１スイング指標と呼ぶことがある。）を算出する第１指標算出工程
（Ｓ６）第１スイング指標に基づいて、ゴルファー７に適したシャフト４０の重量である最適シャフト重量を決定する最適シャフト重量決定工程
（Ｓ７）センサデータに基づいて、第２スイング指標（本実施形態では、後述する第１〜第４特徴量Ｆ1〜Ｆ4）を算出する第２指標算出工程
（Ｓ８）第２スイング指標に基づいて、ゴルファー７に適したシャフト４０の剛性（本実施形態では、後述するＩＦＣ及びフレックス）を示す最適剛性指標を決定する最適剛性決定工程
（Ｓ９）最適シャフト重量及び最適剛性指標に合致するシャフト４０を選択する推奨シャフト選択工程
（Ｓ１０）計測工程で得られた画像データ（計測データ）に基づいて、ヘッド４１の挙動を決定するヘッド挙動決定工程
（Ｓ１１）解析の結果をＧＵＩ画面上に表示する結果表示工程
【００３５】
以下、これらの工程を順に説明する。なお、ＸＹＺ全体座標系は、図１に示すとおりに定義される３軸直交座標系である。すなわち、Ｚ軸は、鉛直下方から上方に向かう方向であり、Ｘ軸は、ゴルファー７の背から腹に向かう方向であり、Ｙ軸は、地平面に平行でボールの打球地点から目標地点に向かう方向である。
【００３６】
＜２−１．計測工程＞
計測工程（Ｓ１）では、ゴルファー７により、上述の慣性センサユニット１付きゴルフクラブ４がスイングされる。以下、計測工程で使用されるゴルフクラブ４を、テストクラブと呼ぶことがある。そして、テストクラブのスイング動作中のグリップ加速度ａ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（スカラー）を算出する。
【数１】
2017217423.tif 000003
【００４１】
また、グリップ挙動導出部２４Ｂは、グリップ速度（ｖpY，ｖpZ）を積分することにより、スイング平面Ｐ内でのグリップ４２の軌道を算出する。さらに、グリップ挙動導出部２４Ｂは、スイング平面Ｐに直交する軸周りのグリップ角速度ωpXを算出する。なお、第２変換工程の具体的な計算方法は、適宜選択することができるが、必要であれば、同出願人らによる特開２０１６−２４２９号公報や特開２０１６−２４３０号公報等に記載の方法に従うことができる。
【００４２】
＜２−４．肩挙動導出工程＞
以下、スイング平面Ｐ内でのグリップの挙動に基づいて、スイング平面Ｐ内の疑似的な肩の挙動を導出する肩挙動導出工程（Ｓ４）について説明する。本実施形態では、ゴルフクラブ４の挙動は、ゴルファー７の肩及びグリップ４２（或いは、これを握るゴルファーの手首）を節点とし、ゴルファー７の腕及びゴルフクラブ４をリンクとする二重振り子モデル（図７参照）に基づいて解析される。ただし、肩の挙動は直接的に実測されるのではなく、実測されたグリップの挙動に基づいて、疑似的な肩の挙動として導出される。以下では、特に断らない限り、単に「肩」という場合も、このような疑似的な肩を意味し得るものとする。疑似的な肩とグリップ４２（手首）との間を直線的に延びるものとして定義される疑似的な「腕」についても同様である。
【００４３】
グリップの挙動から肩の挙動を特定するに当たり、本実施形態に係る二重振り子モデルは、以下の（１）〜（３）を前提とする。図７は、以下の前提条件を概念的に説明する図である。
（１）スイング平面Ｐ上において、グリップ４２（手首）は肩を中心として円運動する。
（２）スイング平面Ｐ上において、肩とグリップ４２との距離（半径）Ｒは、一定である。
（３）肩は、スイング動作中は動かない（ただし、回転する）。
【００４４】
以上の前提の下、肩挙動導出部２４Ｃは、第２変換工程で得られたスイング平面Ｐ内でのグリップ４２の軌道を、円弧（円）に近似する（図８参照）。そして、当該円弧（円）の中心を肩の位置Ｐs＝（ＰsX，ＰsY）とし、当該円弧（円）の中心からグリップ４２の軌道までの平均的な距離を、腕長さ（肩とグリップ４２との距離）Ｒとする。以下、グリップ４２の軌道上の点Ａiの座標を、（Ｘi，Ｙi）と表す（ｉ＝１，２，・・・）。なお、ここでいうｉは、時系列に沿ったデータ番号である。
【００４５】
以下に、グリップ４２の軌道の近似円（円弧）を導出する方法の一例を示す。まず、グリップ４２の軌道上の任意の３点、例えば、Ａi，Ａi+30，Ａi+60を考える。このとき、任意の三角形の外接円の中心は、当該三角形の三辺の垂直二等分線の交点となることから、Ａi，Ａi+30，Ａi+60を頂点とする三角形を考えたとき（図９参照）、以下の数２の式、ひいては数３の式が成り立つ。
【数２】
2017217423.tif 000004
【数３】
2017217423.tif 000005
【００４６】
そして、様々なｉに対する複数の数３の式から以下の数４の式を作成し、疑似逆行列を導出する。これにより、グリップ４２の軌道の近似円（円弧）の中心Ｐs＝（ＰsX，ＰsY）を導出することができる。
【数４】
2017217423.tif 000006
【００４７】
続いて、肩挙動導出部２４Ｃは、近似円（円弧）の中心Ｐs＝（ＰsX，ＰsY）からグリップ４２の軌道上の各点Ａi＝（Ｘi，Ｙi）までの距離の平均値を算出し、腕長さＲとする。そして、肩挙動導出部２４Ｃは、この腕長さＲに基づいて、スイング平面Ｐ内におけるトップからインパクトまでの肩周りの角速度（腕の角速度）ω1＝ＶGE／Ｒを算出する。
すなわち、腕の角速度ω1は、計測によるグリップ速度ＶGEが反映された値となる。
【００４８】
＜２−５．第１指標算出工程＞
以下、図１０を参照しつつ、グリップ４２の挙動及び肩の挙動に基づいて、第１スイング指標を算出する第１指標算出工程（Ｓ５）について説明する。第１スイング指標とは、最適シャフト重量を決定するための指標であり、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態の第１スイング指標は、後述される腕出力パワーＰ1_AVE、及びクラブ入力パワーＰ2_AVEである。
【００４９】
具体的には、まず、ステップＳ３１では、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω1を積分し、トップからインパクトまでの腕の回転角度θ1を算出する。なお、回転角度θ1は、図１１のように定義され、図１１の紙面は、スイング平面Ｐに等しい。以下では、図１１に示されるスイング平面Ｐ内での新たなＸＹ座標系に基づいて、解析が進められる。スイング平面Ｐ内での新たなＸＹ座標系のＸ軸は、上述したＸＹＺ全体座標系のＹ軸に等しく、新たなＸＹ座標系のＹ軸は、ＸＹＺ全体座標系のＺ軸をスイング平面Ｐ内に投影した軸である。
【００５０】
また、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の角速度ω1を微分し、トップからインパクトまでの角加速度ω1'を算出する。次に、肩挙動導出部２４Ｃは、トップからインパクトまでの腕の重心の位置（Ｘ1，Ｙ1）、速度（ＶX1，ＶY1）及び加速度（ＡX1，ＡY1）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。
【数５】
2017217423.tif 000007
【００５１】
ただし、ｒは、肩から腕の重心までの距離である。本実施形態では、腕の重心は、腕の中心にあるものと仮定される。従って、Ｒ＝２ｒである。
【００５２】
次に、ステップＳ３２では、グリップ挙動導出部２４Ｂは、ステップＳ３１と同様の演算をグリップ４２周りについても行う。すなわち、トップからインパクトまでのグリップ角速度ωpX＝グリップ４２周りのゴルフクラブ４の角速度ω2を積分し、トップからインパクトまでのグリップ４２周りのゴルフクラブ４（シャフト４０）の回転角度θ2を算出する。回転角度θ2は、図１１のように定義される。
【００５３】
続いて、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の角速度ω2を微分し、トップからインパクトまでの角加速度ω2'を算出する。次に、グリップ挙動導出部２４Ｂは、トップからインパクトまでのゴルフクラブ４の重心の位置（Ｘ2，Ｙ2）、速度（ＶX2，ＶY2）及び加速度（ＡX2，ＡY2）を算出する。これらの値は、上述した計算結果を以下の式に代入することにより算出される。
【数６】
2017217423.tif 000008
【００５４】
ただし、Ｌは、グリップ４２からゴルフクラブ４の重心までの距離である。Ｌの値は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【００５５】
次に、ステップＳ３３では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの肩に発生する拘束力Ｒ1＝（ＲX1，ＲY1）を算出するとともに、トップからインパクトまでのグリップ４２に発生する拘束力Ｒ2＝（ＲX2，ＲY2）を算出する。以下の式は、並進方向の力の釣り合いに基づくものである。ただし、ｍ1は、腕の質量であり、本実施形態では、腕の質量ｍ1は、適宜予め定められているものとする。例えば、解析を開始する前に、ゴルファー７の体重を入力しておき、入力された体重に所定の係数を掛ける等して、自動的に腕の質量が算出される。ｍ2は、ゴルフクラブ４の質量であり、ｇは、重力加速度である。また、ｍ2は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【数７】
2017217423.tif 000009
【００５６】
続くステップＳ３４では、算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の重心周りのトルクＴg1及びゴルフクラブ４の重心周りのトルクＴg2を算出する。
【数８】
2017217423.tif 000010
ただし、Ｉ1は、腕の重心周りの慣性モーメントであり、Ｉ2は、ゴルフクラブ４の重心周りの慣性モーメントである。本実施形態では、腕の重心周りの慣性モーメントＩ1は、腕の重心が腕の中心にあるとの仮定の下、Ｉ1＝ｍ1ｒ2／３として算出される。また、I2は、ゴルフクラブ４のスペックであり、予め定められているものとする。
【００５７】
続くステップＳ３５では、算出部２４Ｄは、上述した計算結果に基づいて、トップからインパクトまでの腕の仕事率（パワー）Ｅ1'を算出する。具体的には、Ｅ1'は、肩の速度ベクトルをｖsとし、グリップ４２の速度ベクトルをｖgとして、以下の式に従って表される。また、ｖs，ｖgはそれぞれ、肩の位置ベクトルｄs、グリップ４２の位置ベクトルｄg＝ｄs＋（２Ｘ1，２Ｙ1）を一階微分することにより算出可能である。
【数９】
2017217423.tif 000011
【００５８】
また、本実施形態では、肩は動かないため、ｖs＝（０，０）となり、腕の仕事率Ｅ1'は、以下の式に従って算出される。算出部２４Ｄは、上述した計算結果を以下の式に代入することにより、トップからインパクトまでの腕の仕事率Ｅ1'を算出する。
【数１０】
2017217423.tif 000012
【００５９】
ところで、ゴルフスイングにおいて、ゴルフクラブ４の先端（ヘッド４１）を最も加速させるためには、まず腕を十分に加速させて、その後、腕の動きを止めてゴルフクラブ４に勢いを与えることが求められると考えられる。ここでいう腕の加速具合とは、腕が出力するパワー（腕出力パワー）Ｐ1という物理指標に置き換えることができ、ゴルフクラブ４に与える勢いとは、ゴルフクラブ４に入力されるパワー（クラブ入力パワー）Ｐ2という物理指標に置き換えることができる。そして、腕出力パワーＰ1とは、腕の仕事率Ｅ1'を表す数１０の右辺の第２項及び第３項部分に相当する。また、クラブ入力パワーＰ2とは、数１０の式中の右辺の第１項部分に相当する。すなわち、腕出力パワーＰ1及びクラブ入力パワーＰ2は、以下のとおり表すことができる。ステップＳ３５では、算出部２４Ｄは、腕の仕事率Ｅ1'に加え、トップからインパクトまでの腕出力パワーＰ1及びクラブ入力パワーＰ2を算出する。
【数１１】
2017217423.tif 000013
【００６０】
なお、スイング動作中にゴルフクラブ４で発揮される仕事率Ｅ2'は、下式のように表すことができる。すなわち、クラブ入力パワーＰ2＝Ｒ2ｖgTが橋渡しになって、腕からゴルフクラブ４へとエネルギーが伝達される。
【数１２】
2017217423.tif 000014
【００６１】
続くステップＳ３６では、算出部２４Ｄは、トップの時刻ｔtから仕事率Ｅ1'が最大となる時刻ｔmaxまでの腕の仕事量Ｅ1を算出する。腕の仕事量Ｅ1は、時刻ｔt〜ｔmaxの区間で腕の仕事率Ｅ1'を積分することにより、算出される（図１２参照）。なお、仕事量Ｅ1は、時刻ｔt〜ｔmaxの間に腕で発揮される仕事量（エネルギー）を表す指標と考えることができるから、この意味で、スイング動作中の腕エネルギーと呼ぶことができる。また、算出部２４Ｄは、時刻ｔt〜ｔmaxの間に腕で発揮される平均仕事率ＥAVE＝Ｅ1／ｔmax―ｔtを算出する。平均仕事率ＥAVEは、スイング動作中に単位時間当たりに平均的に発揮ないし消費される腕エネルギーである。
【００６２】
また、算出部２４Ｄは、トップの時刻ｔtから腕出力パワーＰ1が最大値をとる時刻ｔmまでの区間で腕出力パワーＰ1を積分し、この積分値Ｄ1をこの積分区間で除算することで、スイング動作中の平均的な腕出力パワーＰ1_AVEを算出する。なお、この積分値Ｄ1は、スイング動作中にゴルファーの腕が行う仕事量であり、腕出力パワーを示す指標となり得る。同様に、算出部２４Ｄは、トップの時刻ｔtからクラブ入力パワーＰ2が最大値をとる時刻ｔnまでの区間でクラブ入力パワーＰ2を積分し、この積分値Ｄ2をこの積分区間で除算することで、スイング動作中の平均的なクラブ入力パワーＰ2_AVEを算出する。なお、この積分値Ｄ2は、スイング動作中のゴルフクラブ４に対してされる仕事量であり、クラブ入力パワーを示す指標となり得る。なお、ここで示した積分区間は例示であり、例えば、時刻ｔtからインパクトの時刻ｔiまでの区間等、適宜設定し得る。
【００６３】
なお、以上の腕出力パワーＰ1，Ｐ1_AVEとは、スイング動作中にゴルファー７がコックを溜める強さを表す指標であると言い換えることができる。また、クラブ入力パワーＰ2，Ｐ2_AVEとは、スイング動作中にゴルファー７がコックを解放する強さを表す指標であると言い換えることができる。
【００６４】
続くステップＳ３７では、算出部２４Ｄは、スイング動作中のコック解放タイミングｔrを算出する。なお、本発明者らは、実験を通して、インパクト時のヘッド速度Ｖhが、スイング動作中のコック解放タイミングｔr、及び、腕エネルギーＥ1又は平均仕事率ＥAVEと相関があることを発見した。そこで、ここでは、インパクト時のヘッド速度Ｖhを算出すべく、コック解放タイミングｔrが算出される。本実施形態では、コック解放タイミングｔrは、時刻ｔt〜ｔiまでの区間で腕の仕事率Ｅ1'が最大となる時刻が、コック解放タイミングｔrとして特定される。
【００６５】
続くステップＳ３８では、算出部２４Ｄは、コック解放タイミングｔr及び腕エネルギーＥAVEに基づいて、インパクト時のヘッド速度Ｖhを算出する。具体的には、インパクト時のヘッド速度Ｖhは、下の式に従って算出される。なお、ｋ1，ｋ2，ｋ3は、予め行われた多数の実験結果から重回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。以上より、指標算出工程が終了する。
Ｖh＝ｋ1・ＥAVE＋ｋ2・ｔr＋ｋ3
【００６６】
＜２−６．最適シャフト重量決定工程＞
以下、図１３を参照しつつ、最適シャフト重量決定工程（Ｓ６）の流れについて説明する。最適シャフト重量決定工程は、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEの大きさに応じて、最適シャフト重量の範囲（以下、最適シャフト重量帯）を決定するステップである。ここでは、Ｐ1_AVE，Ｐ2_AVEの値が大きい程、最適シャフト重量帯は段階的に大きな値に設定される。
【００６７】
まず、ステップＳ４１では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ１の上側にあるか、すなわち、図１４の領域Ａ１に属するか否か（以下、条件１）を判定する。そして、条件１が満たされる場合には、８０ｇ以上が最適シャフト重量帯であると判定する。一方、ステップＳ４１で条件１が満たされない場合には、ステップＳ４２に進む。ステップＳ４２では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ２の上側かつ直線Ｌ１よりも下側にあるか、すなわち、図１４の領域Ａ２に属するか否か（以下、条件２）を判定する。そして、条件２が満たされる場合には、７０ｇ以上、８０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４２で条件２が満たされない場合には、ステップＳ４３に進む。ステップＳ４３では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ３の上側かつ直線Ｌ２よりも下側にあるか、すなわち、図１４の領域Ａ３に属するか否か（以下、条件３）を判定する。そして、条件３が満たされる場合には、６０ｇ以上、７０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４３で条件３が満たされない場合には、ステップＳ４４に進む。ステップＳ４４では、決定部２４Ｅは、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ４の上側かつ直線Ｌ３よりも下側にあるか、すなわち、図１４の領域Ａ４に属するか否か（以下、条件４）を判定する。そして、条件４が満たされる場合には、５０ｇ以上、６０ｇ未満が最適シャフト重量帯であると判定する。一方、ステップＳ４４で条件４が満たされない場合、すなわち、第１指標算出工程で算出された（Ｐ1_AVE，Ｐ2_AVE）で表される点が、図１４に示すＰ1_AVE−Ｐ2_AVE平面において直線Ｌ４よりも下側にあり、図１４の領域Ａ５に属する場合には、５５ｇ以下が最適シャフト重量帯であると判定する。
【００６８】
以上の最適シャフト重量決定工程は、以下の知見に基づく。すなわち、本発明者らは、多数のゴルファーにテストクラブをスイングさせ、このときの腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEを算出した。また、同じゴルファーに、様々なシャフト重量のゴルフクラブをスイングさせ、最大の飛距離を与えるシャフト重量を算出し、これを最適シャフト重量とした。
【００６９】
上記の実験からは、最適シャフト重量は、同じゴルファーによるテストクラブのスイング時の腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEが共に大きいほど、大きくなる傾向にあることが確認された。その結果、本発明者らは、腕出力パワーＰ1_AVE−クラブ入力パワーＰ2_AVE空間を、図１４に示すように領域分割することで、最適シャフト重量帯に対応する領域Ａ１〜Ａ５を定義可能であることを発見した。すなわち、上述のステップＳ４１〜Ｓ４４は、腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEを示す点が、Ｐ1_AVE−Ｐ2_AVE空間におけるどの領域にプロットされるかに応じて、最適シャフト重量帯を判定するステップとなっている。ステップＳ４１〜Ｓ４４での判定に使用される閾値、言い換えると、図１４に示す分割領域Ａ１〜Ａ５の境界線Ｌ１〜Ｌ４を特定する情報等は、対応関係データ２８として記憶部２３内に格納されている。すなわち、対応関係データ２８とは、Ｐ1_AVE,Ｐ2_AVEの大きさと最適シャフト重量帯との対応関係を定めるデータである。なお、境界線Ｌ１〜Ｌ４は互いに概ね平行であり、いずれもＰ1_AVE−Ｐ2_AVE平面において負の傾きを持つ直線である。ステップＳ４１〜Ｓ４４では、この記憶部２３内の対応関係データ２８が参照され、上記の判定が行われる。なお、図２では、対応関係データ２８は、プログラム３とは別のデータとして示されているが、プログラム３内に組み込まれていてもよい。
【００７０】
＜２−７．第２指標算出工程＞
以下、計測工程で得られたセンサデータに基づいて、第２スイング指標を算出する第２指標算出工程（Ｓ７）について説明する。第２スイング指標とは、最適剛性指標を決定するための指標であり、ゴルファー７によるスイング動作を特徴付ける特徴量である。本実施形態では、第２スイング指標として、後述される第１〜第４特徴量Ｆ1〜Ｆ4が算出される。
【００７１】
第１〜第４特徴量Ｆ1〜Ｆ4について理解するためには、まず、最適剛性指標について理解することが重要である。最適剛性指標とは、ゴルファー７に適したシャフト４０の剛性を示す指標のことであり、本実施形態では、シャフト４０の剛性は、シャフト４０の複数の位置における曲げ剛性の分布（以下、ＥＩ分布）として評価される。本実施形態に係るＥＩ分布は、定量的に数値を用いて表現され、より具体的には、インターナショナル・フレックス・コード（ＩＦＣ）を用いて算出される。そのため、まず、このＩＦＣについて説明する。なお、ＩＦＣは、本出願人により広く提案されているシャフトの特性を示す公知の指標であり、例えば、特許文献１をはじめとして、既に様々な文献で詳しく説明されている。従って、ここで改めて説明する必要は必ずしもないが、参考のため、ここでも説明を行う。
【００７２】
ＩＦＣは、図１５に示すとおり、シャフト４０の延びる方向に沿った４つの位置Ｈ１〜Ｈ４におけるシャフト４０の曲げ剛性をそれぞれ０〜９の１桁の数値で表し、この４つの数値をシャフト４０の延びる方向に沿って配列したコードである。より具体的には、シャフト４０のバット端からチップ端に向かってこの順に概ね一定間隔で、４つの測定点Ｈ１〜Ｈ４が定義される。例えば、シャフト４０のチップ端から３６インチの箇所を測定点Ｈ１とし、２６インチの箇所を測定点Ｈ２とし、１６インチの箇所を測定点Ｈ３とし、６インチの箇所を測定点Ｈ４とすることができる。そして、これらの４つの測定点Ｈ１〜Ｈ４のそれぞれにおける曲げ剛性の値（以下、ＥＩ値）Ｊ1〜Ｊ4が計測される。
【００７３】
シャフト４０の各測定点Ｈ（Ｈ１〜Ｈ４）におけるＥＩ値（Ｎ・ｍ2）は、様々な方法で測定することができ、例えば、インテスコ社製の２０２０型計測機（最大荷重５００ｋｇｆ）を用いて図１６に示すようにして測定することができる。この測定方法では、２つの支持点１１１，１１２でシャフト４０を下方から支持しつつ、測定点Ｈに上方から荷重Ｆを加えたときのたわみ量を測定する。支持点１１１と支持点１１２との間の距離（スパン）は、例えば、２００ｍｍとすることができ、測定点Ｈは、支持点１１１と支持点１１２の中間点とすることができる。より具体的には、支持点１１１，１１２を支える支持体１１４，１１５を固定した状態で、測定点Ｈにおいて圧子１１３を一定速度（例えば、５ｍｍ／分）で下方へ移動させる。そして、荷重Ｆが２０ｋｇｆに達した時点で圧子１１３の移動を終了させ、この瞬間のシャフト４０のたわみ量（ｍｍ）を測定し、このたわみ量をＥＩ値（Ｎ・ｍ2）に換算する。
【００７４】
次に、以上の４つの測定点Ｈ１〜Ｈ４におけるＥＩ値Ｊ1〜Ｊ4を、それぞれ１０段階のランク値Ｋ1〜Ｋ4に変換する。具体的には、ランク値Ｋ1〜Ｋ4は、それぞれ測定点Ｈ１〜Ｈ４用の以下の変換表（表１〜表４）に従って、ＥＩ値Ｊ1〜Ｊ4から算出することができる（表１〜４中、変換後のランク値をＩＦＣの欄に示している）。そして、このようにして測定点Ｈ１〜Ｈ４にそれぞれ付与された４つのランク値Ｋ1〜Ｋ4を、よりバット側に対応する値がより左に、よりチップ側に対応する値がより右にくるように配列する。こうして得られた４桁のコードが、ＩＦＣである。ＩＦＣでは、各桁の数値が大きい程、対応する位置での剛性が高いことを意味する。
【表１】
2017217423.tif 000015
【表２】
2017217423.tif 000016
【表３】
2017217423.tif 000017
【表４】
2017217423.tif 000018
【００７５】
第２指標算出工程では、算出部２４Ｄにより、第１〜第４特徴量Ｆ1〜Ｆ4が算出される。本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4は、それぞれゴルファー７に適したＥＩ値Ｊ1〜Ｊ4である最適ＥＩ値ＪS1〜ＪS4、ひいてはゴルファー７に適したランク値Ｋ1〜Ｋ4である最適ランク値ＫS1〜ＫS4を決定するための指標である。そのため、本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4としては、それぞれ最適ＥＩ値ＪS1〜ＪS4と相関を有する特徴量が選択される。また、本実施形態では、第１〜第４特徴量Ｆ1〜Ｆ4として以下の指標が用いられるが、第２スイング指標としては、最適剛性指標との相関が認められる限り、その他の任意の特徴量を用いることができる。
【００７６】
第１特徴量Ｆ1は、トップ付近のコック方向の角速度ωyの傾きであり、例えばトップから５０ｍｓ前の角速度ωyの絶対値と、トップから５０ｍｓ後の角速度ωyの絶対値との和で表すことができる。
【００７７】
第２特徴量Ｆ2は、トップから、角速度ωyが最大となる時点までの当該角速度ωyの平均値である。第２特徴量Ｆ2は、まず、トップからインパクトまでの間で角速度ωyが最大となる時点を求め、トップからこの時点までの角速度ωyの累積値を、トップからこの時点までの時間で除することにより算出される。
【００７８】
第３特徴量Ｆ3は、角速度ωyが最大となる時点からインパクトまでの当該角速度ωyの平均値である。第３特徴量Ｆ3は、角速度ωyが最大となる時点からインパクトまでの角速度ωyの累積値を、角速度ωyが最大となる時点からインパクトまでの時間で除することにより算出される。
【００７９】
第４特徴量Ｆ4は、トップからインパクトまでの角速度ωyの平均値であり、トップからインパクトまでの角速度ωyの累積値を、トップからインパクトまでの時間で除することにより算出される。
【００８０】
ところで、スイング動作中、ゴルフクラブ４のシャフト４０は、当該ゴルフクラブの先端に比較的重量が大きいヘッド４１が存在するため、その慣性により曲げが生じる。この曲げは、スイングの全過程において、シャフト４０の同一箇所に生じるのではなく、図１７に示されるように、トップからインパクトに向けてシャフト４０の手元側から先端側に伝わる。換言すれば、トップからインパクトに向けてスイングが進行するにしたがい、シャフト４０における曲げの位置が当該シャフト４０の手元側から先端側に移動する。
【００８１】
より具体的には、アドレスからテイクバックを行い、トップに至った時点（図１７において（１）で示される時点）では、シャフト４０の手元付近に曲げが生じる。ついで、切り返しを行い、ダウンスイング初期（図１７において（２）で示される時点）に至ると、曲げはシャフト４０の先端側にやや移動する。さらに、ゴルファー７の腕が水平になる時点（図１７において（３）で示される時点）では、曲げはシャフト４０中央よりも先端側に移動する。そして、インパクト直前（図１７において（４）で示される時点）では、曲げはシャフト４０の先端付近まで移動する。
【００８２】
従って、第１〜第４特徴量Ｆ1〜Ｆ4は、それぞれスイング動作中のトップ付近からインパクト付近までの間の第１〜第４区間において算出することができる。また、ここでの第１〜第３区間は、この順に時間経過に沿った区間であり、互いに一部重複する又は重複することのない区間となっている。
【００８３】
＜２−８．最適剛性決定工程＞
次に、最適剛性決定工程（Ｓ８）の流れについて説明する。本工程では、決定部２４Ｅが、第２スイング指標（第１〜第４特徴量Ｆ1〜Ｆ4）と最適剛性指標（最適ＥＩ値ＪS1〜ＪS4）との相関関係を表す予め定められた近似式に従って、最適剛性指標（最適ＥＩ値ＪS1〜ＪS4）を決定する。本実施形態に係る近似式は、線形近似式であり、以下のように表される。
ＪS1＝ａ1・Ｆ1＋ｂ1
ＪS2＝ａ2・Ｆ2＋ｂ2
ＪS3＝ａ3・Ｆ3＋ｂ3
ＪS4＝ａ4・Ｆ4＋ｂ4
【００８４】
決定部２４Ｅは、第２指標算出工程で算出された第１〜第４特徴量Ｆ1〜Ｆ4をかかる近似式に代入することにより、最適ＥＩ値ＪS1〜ＪS4を算出する。また、決定部２４Ｅは、上述の表１〜表４の換算表に従って、最適ＥＩ値ＪS1〜ＪS4をそれぞれ最適ランク値ＫS1〜ＫS4に変換する。
【００８５】
なお、上式中、ａ１〜ａ４及びｂ１〜ｂ４は、予め行われた多数の実験結果から回帰分析により得られた定数であり、記憶部２３内に予め保持されている値である。ここでいう実験とは、例えば、特許文献１と同様に、以下のように行うことができる。すなわち、まず、多数のゴルファーの各々に複数のゴルフクラブを振ってもらい、そのときの飛距離、打球の方向性（左右ずれ）及び官能試験による振り易さを数値化する。そして、その数値から各ゴルファーに適しているゴルフクラブを決定し、当該ゴルフクラブのＥＩ値を当該ゴルファーの最適ＥＩ値とする。また、上記と同様の方法で、各ゴルファーの第１〜第４特徴量Ｆ1〜Ｆ4を算出する。そして、かかる実験の後、多数のゴルファー分の第１〜第４特徴量Ｆ1〜Ｆ4及び最適ＥＩ値のデータを回帰分析することで、ａ１〜ａ４及びｂ１〜ｂ４が算出される。
【００８６】
また、より信頼性の高い近似式とするために、ａ１〜ａ４及びｂ１〜ｂ４の値を条件に応じて変更することができる。例えば、ヘッド速度Ｖhに応じて、近似式を用意することができる。一例としては、上記実験データを、ヘッド速度帯に応じて分類し（例えば、４５ｍ／ｓ以上、４１〜４５ｍ／ｓ、４１ｍ／ｓ以下）、同じ分類に属するデータのみを対象に上記近似式を作成し、ａ１〜ａ４及びｂ１〜ｂ４を決定しておくことができる。そして、最適剛性決定工程では、ゴルファー７のヘッド速度Ｖhがどのヘッド速度帯に属するかを判断し、当該ヘッド速度帯に対応する近似式を用いて、最適剛性指標である最適ランク値ＫS1〜ＫS4を算出する。
【００８７】
また、本実施形態では、最適剛性指標として、最適ランク値ＫS1〜ＫS4だけでなく、最適フレックスも決定される。具体的には、決定部２４Ｅは、最適ランク値ＫS1〜ＫS4に基づいて、最適フレックスを決定する。フレックスとは、シャフト４０全体での硬さ（曲げ剛性）を評価する指標である。従って、シャフト４０の複数の位置におけるゴルファー７に適した曲げ剛性を表す最適ランク値ＫS1〜ＫS4が分かれば、これらの値に基づいて最適フレックスを算出することができる。例えば、特定の位置での最適ランク値を最適フレックスとすることもできるし、複数の位置での最適ランク値の平均値を最適フレックスをとすることもできる。フレックスには、例えば、「ＳＲ」「Ｓ」「Ｘ」などの様々な種類が存在するが、最適フレックスは、このうちの１つに特定される必要はなく、「ＳＲ」か「Ｓ」や、「Ｓ」か「Ｘ」等、幅を持つように設定することができる。
【００８８】
なお、他の実施形態では、最適フレックスを最適ランク値ＫS1〜ＫS4ではなく、最適ＥＩ値ＪS1〜ＪS4から直接算出することもできるし、これらの値ＪS1〜ＪS4，ＫS1〜ＫS4に基づかずに算出することもきる。後者の場合、最適フレックスを特定可能な適切な特徴量を第２スイング指標決定工程で算出しておけばよい。
【００８９】
＜２−９．推奨シャフト選択工程＞
以上の工程Ｓ１〜Ｓ８により、最適シャフト重量帯及び最適剛性指標（最適ＥＩ値ＪS1〜ＪS4、最適ランク値ＫS1〜ＫS4及び最適フレックス）が決定されると、選択部２４Ｆは、推奨シャフト選択工程（Ｓ９）を実行する。本工程では、シャフトＤＢ２９内に登録されている多数のシャフトの中から、ゴルファー７に適したシャフト４０(以下、推奨シャフト)が特定される。また、本実施形態では、推奨シャフトとして、最適シャフト重量帯を基準としたシャフトである第１推奨シャフトと、ゴルファー７の普段使用しているゴルフクラブの重量（以下、マイクラブ重量）を基準とした推奨シャフトである第２推奨シャフトとが決定される。
【００９０】
第１推奨シャフトの決定方法について説明する。まず、選択部２４Ｆは、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。シャフトＤＢ２９内に登録されているシャフト４０のスペックを示す情報には、製造メーカー、型番、４つの位置Ｈ１〜Ｈ４におけるＥＩ値Ｊ1〜Ｊ4及びランク値Ｋ1〜Ｋ4（ＩＦＣ）、シャフト４０の重量、フレックス、トルク、調子等が含まれる。そして、選択部２４Ｆは、これらの情報を参照して、シャフトＤＢ２９内に登録されている全てのシャフト４０の中から、第１絞り込みシャフトを特定する。第１絞り込みシャフトとは、シャフト４０の重量が最適シャフト重量帯に属し、かつ、フレックスが最適フレックスに一致するシャフト４０である。なお、第１絞り込みシャフトは、通常多数本存在する。
【００９１】
続いて、選択部２４Ｆは、各第１絞り込みシャフトについて、当該シャフトのランク値Ｋ1〜Ｋ4と、最適剛性決定工程で決定された最適ランク値ＫS1〜ＫS4との一致度を算出し、第１推奨シャフトとして、一致度の高いシャフトを特定する。一致度は、例えば、以下の数２５の式に従って算出することができ、値が小さいほど一致度が高い。
【数１３】
2017217423.tif 000019
【００９２】
次に、第２推奨シャフトの決定方法について説明する。まず、選択部２４Ｆは、推奨ゴルフクラブに用いられるべきヘッド４１(以下、推奨ヘッド)の種類を決定する。推奨ヘッドの種類の決定は、本明細書では説明されないフィッティング処理により行うこともできるし、表示部２１及び入力部２２を介してユーザに質問する等してお好みのヘッド４１を選択させることにより行うこともできる。また、選択部２４Ｆは、表示部２１及び入力部２２を介してユーザに質問する等して、マイクラブ重量を特定する。
【００９３】
続いて、選択部２４Ｆは、推奨ヘッドのスペックを示す情報をヘッドＤＢ２７内から読み出すととともに、シャフトＤＢ２９内に登録されている全てのシャフト４０のスペックを示す情報を読み出す。ヘッドＤＢ２７内に登録されているヘッド４１のスペックを示す情報には、製造メーカー、型番及び重量等が含まれる。そして、選択部２４Ｆは、これらの情報を参照して、シャフトＤＢ２９内に登録されている全てのシャフト４０の中から、第２絞り込みシャフトを特定する。第２絞り込みシャフトとは、最適フレックスに一致するフレックスを有し、推奨ヘッドと組み合わせたときの重量がマイクラブ重量±所定値（グリップの重量及びフェラルの重量等に誤差範囲を加えた値）の範囲内となる重量を有するシャフト４０である。なお、第２絞り込みシャフトは、通常多数本存在する。
【００９４】
続いて、選択部２４Ｆは、各第２絞り込みシャフトについて、当該シャフトのランク値Ｋ1〜Ｋ4と、最適剛性決定工程で決定された最適ランク値ＫS1〜ＫS4との一致度を算出し、第２推奨シャフトとして、一致度の高いシャフトを特定する。一致度は、数１３の式に従って算出することができる。なお、第１推奨シャフトとしては、１本のみ特定されてもよいし、複数本が特定されてもよい。第２推奨シャフトについても、同様である。
【００９５】
＜２−１０．ヘッド挙動決定工程＞
続いて、決定部２４Ｅは、計測工程で得られた画像データに基づいて、ヘッド４１の挙動を決定するヘッド挙動決定工程（Ｓ１０）を実行する。本実施形態では、ヘッド４１の挙動として、インパクト直前のヘッド速度Ｖhi、インパクト時のフェース角、ヘッド軌道の角度、打点を算出する。なお、ここでいうヘッド軌道とは、アウトサイドイン、インサイドアウト、インサイドイン等の軌道の種類を意味し、ヘッド軌道の角度として、定量的な角度が算出される。打点は、ヘッド４１のフェース面４１ａ上におけるボールの打点である。
【００９６】
インパクト直前のヘッド速度Ｖhiは、以下のように算出される。すなわち、投光器５５Ａ，５５Ｂと受光器５６Ａ，５６Ｂとの間隔は既知である。従って、タイミングｔ１，ｔ２が分かれば、インパクト直前のヘッド速度Ｖhiを算出することができるため、決定部２４Ｅは、タイミングｔ１，ｔ２の情報に基づいて、ヘッド速度Ｖhiを算出する。
【００９７】
また、インパクト時のフェース角及びヘッド軌道は、画像データを画像処理することにより決定される。図１８は、ヘッド４１をフェース面４１ａ側から視た図である。同図に示すように、本実施形態に係るヘッド４１のクラウン部４１ｂには、フェース面４１ａに沿うように帯状のマーカーＭ１が貼付されている。マーカーＭ１は、ストロボ５３，５４からの光を効率的に反射する素材で形成されている。従って、カメラ５１により撮影された画像上においては、マーカーＭ１の領域、言い換えると、ヘッド４１の平面視においてフェース面４１ａに沿った帯状の領域が鮮明に写り込むことになる。決定部２４Ｅは、記憶部２５内に記憶されているインパクト直前の２つのタイミングにおける画像（ストロボ５３，５４の発光のタイミングでの２枚の画像）から、マーカーＭ１の像を抽出する。そして、これらのマーカーＭ１の像に基づいて、インパクト時のフェース角を推定する。また、これらのマーカーＭ１の像の重心を結ぶベクトルと、飛球線方向のベクトルとの為す角度を算出し、この角度からアウトサイドイン等の軌道の種類を特定するとともに、この角度をヘッド軌道の角度として決定する。
【００９８】
また、決定部２４Ｅは、画像データを画像処理にすることにより、打点を決定する。より具体的には、インパクト時又はその付近の画像において、ボール及びフェース面４１ａの画像を抽出し、これらの位置関係から打点を算出する。本実施形態では、フェース面４１ａの幾何学中心であるフェースセンターから、打点がトゥ−ヒール方向にどれ位ずれているかが定量的に算出される。
【００９９】
＜２−１１．結果表示工程＞
次に、以上の解析の結果を示すＧＵＩ画面を作成し、これを表示する結果表示工程（Ｓ１１）について説明する。本実施形態では、ＧＵＩ画面として、図１９〜図２４に示すような画面Ｗ１〜Ｗ５が表示部２１上に表示される。
【０１００】
表示制御部２４Ｇは、入力部２２を介したユーザの操作を受けて、画面Ｗ１〜Ｗ５を適宜切り替えながら、表示部２１上に表示させる。本実施形態では、画面Ｗ１〜Ｗ４は、図１９、図２０、図２２及び図２４に示すようにタブ形式で表示され、現在選択されているタブに対応する画面が表示部２１上に表示される。
【０１０１】
図１９に示す画面Ｗ１は、「スイングの傾向」というタブに対応する画面であり、文字どおり、主としてゴルファー７のスイングの傾向を表示する。まず、画面Ｗ１上の左半分の領域Ｃ１内では、「コックの動きの速さ度合い」、「手の甲の動き」、「手のローテーション」及び「コックの使い方」が、スライダー形式で定量的に表示される。従って、画面Ｗ１を見たユーザは、スライダーの位置といった図形的な観点から直感的に、各種指標の大小を正しく評価することができる。
【０１０２】
「コックの動きの速さ度合い」としては、トップ付近、ダウンスイング前半、ダウンスイング後半及びダウンスイング全体の４つの期間における値が表示され、スライダーが左に配置されるほどこの度合いが小さく、右に配置されるほど大きいことを意味している。「コックの動きの速さ度合い」は、対応する期間の角速度ωyに基づいて決定部２４Ｅにより決定される。「手の甲の動き」としては、トップ付近の値が表示され、スライダーが左に配置されるほど手の甲が開いており、右に配置されるほど閉じていることを意味している。「手の甲の動き」は、トップ付近の角速度ωxに基づいて決定部２４Ｅにより決定される。「手のローテーション」としては、ダウンスイング全体及びインパクト直前の２つの期間における値が表示され、スライダーが左に配置されるほど手のローテーションが少なく、右に配置されるほど多いことを意味している。「手のローテーション」は、対応する期間における角速度ωzに基づいて決定部２４Ｅにより決定される。
【０１０３】
「コックの使い方」としては、「コックを溜める強さ」及び「コックを解放する強さ」が表示され、スライダーが左に配置されるほどコックの使い方が弱く、右に配置されるほど強いことを意味している。「コックを溜める強さ」は、上述したゴルファー７の腕出力パワーＰ1_AVEに応じて決定部２４Ｅにより決定される。「コックを解放する強さ」は、上述したゴルファー７のクラブ入力パワーＰ2_AVEに応じて決定部２４Ｅにより決定される。
【０１０４】
画面Ｗ１上の右上には、ゴルファー７のスイングタイプが表示される。ゴルファー７のスイングタイプとは、最適ＥＩ値ＪS1〜ＪS4（又は、最適ランク値ＫS1〜ＫS4）、すなわち、ゴルファー７に適したシャフト４０のＩＦＣ（以下、理想ＩＦＣという）に基づいて、決定部２４Ｅにより決定される。より具体的には、スイングタイプは、理想ＩＦＣの４つのコードが、シャフト４０の手元から手先にかけて右上がりであるか、右下がりであるか、中央が突出しているか、一定であるか等を表す。ここでは、横軸をシャフト４０の長さ方向の位置とし、縦軸をシャフト４０の硬さをとしたグラフ領域内において、理想ＩＦＣをグラフ状に表す画像Ｃ２が表示される。ただし、スイングタイプは、ゴルファー７のタイプを知るための指標であるため、シャフトの位置とＩＦＣのコードとが厳密に描画されるのではなく、模式化されたグラフが表示される。
【０１０５】
また、画面Ｗ１には、「平均ヘッドスピード」として、上述したヘッドスピードＶh（平均値）が表示され、「推奨シャフト硬度」として、最適フレックスが表示され、「推奨重量帯」として、最適シャフト重量帯が表示される。ヘッドスピードＶh（平均値）は、数値で表示される。一方、「推奨シャフト硬度」及び「推奨重量帯」の欄には、それぞれ、シャフトＤＢ２９内に登録されているシャフト４０のフレックス及びシャフト重量の範囲がスケールの形態で表示されており、当該スケールにおいてゴルファー７に対応する部分が強調表示される。これにより、ユーザは、スケールを参照して、ゴルファーに対応する最適フレックス及び最適シャフト重量帯を直感的に理解することができる。
【０１０６】
図２０に示す画面Ｗ２は、「理想のＩＦＣ」というタブに対応する画面であり、主として理想ＩＦＣ及び推奨シャフトを表示する。画面Ｗ２には、理想ＩＦＣのコードが数値で表示され、「推奨シャフトタイプ」として、上述したスイングタイプが文字で表示される。また、画面Ｗ２には、「推奨シャフト１」及び「推奨シャフト２」として、上述した上位２本の第２推奨シャフトを特定する情報を表示する欄Ｃ３，Ｃ３が表示される。また、「推奨シャフトＷ」として、上述した１本の第１推奨シャフトを特定する情報を表示する欄Ｃ４が表示される。より具体的には、欄Ｃ３，Ｃ３，Ｃ４内には、製造メーカー、型番、質量、フレックス等が表示される。なお、これらの情報は、シャフトＤＢ２９内を参照して特定される。
【０１０７】
また、画面Ｗ２には、以上の理想ＩＦＣに加え、上述した１本の第１推奨シャフト及び２本の第２推奨シャフトのＩＦＣをグラフ状に表示するためのグラフ領域Ｃ５が表示される。これらのグラフは、ＩＦＣのコードに対応する位置に配置され、当該コードが数値で記述されたプロット点と、これらのプロット点を連結する直線部分とから構成される。グラフ領域Ｃ５内では、理想ＩＦＣのグラフは、３本の推奨シャフトに対応するＩＦＣのグラフとは異なる態様で、好ましくは、ユーザの注意をより引きつけるような態様（図２０では、グラフの線が太く、プロット点が大きい）で表示される。また、どのグラフがどの欄Ｃ３，Ｃ３，Ｃ４により特定されるシャフトに対応するかが一目で分かるように、欄Ｃ３，Ｃ３，Ｃ４の装飾と、これらに対応するグラフの装飾とに統一感を持たせることが好ましい。例えば、欄Ｃ３，Ｃ３，Ｃ４内の文字の色と、グラフの色とを一致させることができる。
【０１０８】
欄Ｃ３，Ｃ３，Ｃ４は、排他的に選択可能に構成されている。そして、領域Ｃ６には、現在選択されている欄に対応する推奨シャフトのデザイン画像が表示される。これにより、現在選択されている推奨シャフトが、どのシャフトであるのかを直感的に理解することができる。また、グラフ領域Ｃ５内では、現在選択されている推奨シャフトのＩＦＣのグラフは、その他の推奨シャフトのＩＦＣのグラフと異なる態様で、好ましくは、ユーザの注意をより引きつけるような態様（例えば、前者の方が色が濃く、後者の方が色が薄い）で表示される。
【０１０９】
ところで、図２０の例では、画面Ｗ２の左側に、「コックの動きの速さ度合い」、「手の甲の動き」及び「手のローテーション」が、図１９と同じ形式で表示されている。そして、画面Ｗ１，Ｗ２上の「コックの動きの速さ度合い」を示す欄には、アイコンＣ７が表示されており、これが選択されると、図２１に示す画面Ｗ５がポップアップにより表示される。この画面Ｗ５上には、コックの動きの速さ度合いの意味するところを説明する画像が表示される。より具体的には、この画面Ｗ５上では、コックの動きとは手のどのような方向の動きであるのかが、手の絵、シャフトの絵、及び動きの方向を表す矢印等を用いて図形的に説明される。
【０１１０】
次に、図２２に示す画面Ｗ３は、は、「推奨シャフト重量」というタブに対応する画面であり、文字どおり、主としてゴルファー７に推奨されるべきシャフト４０の重量帯を表示する。より具体的には、画面Ｗ３には、画面Ｗ１と同様に、ゴルファー７の最適シャフト重量帯を表示する「推奨重量帯」の欄Ｃ８が表示される。表示形態も画面Ｗ１と同様であり、欄Ｃ８には、シャフトＤＢ２９内に登録されているシャフト重量の範囲がスケールの形態で表示されており、当該スケールにおいて最適シャフト重量帯の部分が強調表示されている。これにより、ユーザは、スケールを参照して、最適シャフト重量帯を直感的に理解することができる。
【０１１１】
また、画面Ｗ３には、最適シャフト重量帯についてのユーザの納得性を高めるべく、その算出の根拠を説明するための様々な欄Ｃ９〜Ｃ１０が表示される。欄Ｃ１０には、最適シャフト重量帯を決定するための指標であるゴルファー７の腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEが表示される。より具体的には、これらの指標Ｐ1_AVE，Ｐ2_AVEが、それぞれ「コックを溜める強さ」及び「コックを解放する強さ」として、画面Ｗ１と同様にスライダー形式で表示される。これにより、ユーザは、最適シャフト重量帯を決定する根拠となる、ゴルファー７固有のコックを溜める強さ及びコックを解放する強さの大小を、直感的に正しく理解することができる。
【０１１２】
欄Ｃ９には、図１４に示したような、腕出力パワーＰ1_AVEを横軸とし、クラブ入力パワーＰ2_AVEを縦軸とするＰ1_AVE−Ｐ2_AVE平面であるプロット領域が表示され、このプロット領域内には、ゴルファー７の腕出力パワーＰ1_AVE及びクラブ入力パワーＰ2_AVEに対応する点がプロットされ、表示される。欄Ｃ９では、プロット領域の横軸及び縦軸が、それぞれ「コックを溜める強さ」「コックを解放する強さ」を表すことをが、軸に添えられた文字により明確に示されている。また、このプロット領域は、ゴルファー７に推奨すべきシャフト重量帯にそれぞれ対応する、上述した領域Ａ１〜Ａ５に分割されている。これらの領域は、プロット領域内において、互いに異なるシャフト重量帯に対応する領域であることが分かるように異なる態様で表示される。例えば、領域間に境界線が引かれたり、異なる領域に異なる色や模様が付される。また、各領域には、対応するシャフト重量帯に属するシャフト４０の仕様を表す記号が表示される。図２２の例では、「８Ｘ・Ｓ」「７Ｘ・Ｓ」等の記号が表示されているが、これらの８や７等の数値は、シャフト重量帯を意味しており、Ｘ、Ｓ等のアルファベットは、フレックスを示している。
【０１１３】
ところで、コックを溜める強さ及びコックを解放する強さという概念は、全てのユーザが知っているものではない。そのため、欄Ｃ１１には、「スイング中に発生する力、腕とクラブの関係」を示す画像が表示される。本実施形態では、コックを溜める強さ及びコックを解放する強さを説明するための動画が表示される。この動画は、ゴルファー７の腕を模式的に示す図形（以下、腕図形）と、ゴルファー７により把持されたゴルフクラブ４を模式的に示す図形（以下、クラブ図形）とが、スイング動作に伴って二重振り子のように移動する様子を示す動画である。
【０１１４】
今、スイング動作中のタイミングを、時系列に沿って序数詞ｉ＝１，２，・・・，Ｎ（Ｎは、２以上の整数）で表すこととする。このとき、欄Ｃ１１内で表示される動画は、時系列に沿ってコマ送りの形式で再生されるＮ枚の画像Ｉ（ｉ）を含む。画像Ｉ（ｉ）上には、第ｉ番目のタイミングにおいて腕が存在する位置に腕図形Ｇ3（ｉ）が表示され、第ｉ番目のタイミングにおいてゴルフクラブ４が存在する位置にクラブ図形Ｇ4（ｉ）が表示される。図２３Ａ〜図２３Ｆは、画像Ｉ（ｉ）の例を示しており、図２３Ａ〜図２３Ｆの画像は、この順に時系列に沿っている。
【０１１５】
また、画像Ｉ（ｉ）上には、腕図形Ｇ3（ｉ）及びクラブ図形Ｇ4（ｉ）に加えて、第ｉ番目のタイミングにおけるコックを溜める強さを表す第１図形Ｇ1（ｉ）と、同じく第ｉ番目のタイミングにおけるコックを解放する強さを表す第２図形Ｇ2（ｉ）とが表示される。同じタイミングに対応するこれらの図形Ｇ1（ｉ）〜Ｇ4（ｉ）は、１枚の画像Ｉ（ｉ）上において同時に表示される。第ｉ番目のタイミングにおけるコックを溜める強さとは、第ｉ番目のタイミングでの腕出力パワーＰ1であり、第ｉ番目のタイミングにおけるコックを解放する強さとは、第ｉ番目のタイミングでのクラブ入力パワーＰ2である。画像Ｉ（ｉ）上の第１図形Ｇ1（ｉ）は、第ｉ番目のタイミングでの腕出力パワーＰ1の大きさに応じた形態をとり、第２図形Ｇ2（ｉ）は、第ｉ番目のタイミングでのクラブ入力パワーＰ2の大きさに応じた形態をとる。より具体的には、本実施形態では、腕出力パワーＰ1が大きいタイミングである程、第１図形Ｇ1（ｉ）が大きく表示され、クラブ入力パワーＰ2が大きいタイミングである程、第２図形Ｇ2（ｉ）が大きく表示される。なお、腕出力パワーＰ1及びクラブ入力パワーＰ2の大きさは、第１図形Ｇ1（ｉ）及び第２図形Ｇ2（ｉ）の色等を変化させることにより表現されてもよい。また、第１図形Ｇ1（ｉ）及び第２図形Ｇ2（ｉ）の形態を決定する腕出力パワーＰ1及びクラブ入力パワーＰ2は、計測データに基づくゴルファー７固有の値である必要はなく、一般的なスイング動作時の腕出力パワーＰ1及びクラブ入力パワーＰ2とすることができる。この場合、欄Ｃ１１内では、ゴルファー７に関わらず同じ動画が再生される。
【０１１６】
ところで、観念上、コックを溜める力は、腕に関連付けられ、コックを解放する力は、ゴルフクラブ４に関連付けられる。従って、本実施形態では、第１図形Ｇ1（ｉ）は、腕図形Ｇ3（ｉ）に重ねて配置され、第２図形Ｇ2（ｉ）は、クラブ図形Ｇ4（ｉ）に重ねて配置される。ただし、第１図形Ｇ1（ｉ）及び第２図形Ｇ2（ｉ）の配置はこれに限られず、別の好ましい例としては、第１図形Ｇ1（ｉ）を腕図形Ｇ3（ｉ）の近傍に配置し、第２図形Ｇ2（ｉ）をクラブ図形Ｇ4（ｉ）の近傍に配置することもできる。また、第１図形Ｇ1（ｉ）と腕図形Ｇ3（ｉ）とを概ね一体の図形として表現することもできる。この場合、例えば、腕を模した図形の形態を、腕出力パワーＰ1の大きさに応じて変化させるように構成することができる。第２図形Ｇ2（ｉ）及びクラブ図形Ｇ4（ｉ）についても同様であり、この場合、例えば、ゴルフクラブ４を模した図形の形態を、クラブ入力パワーＰ2の大きさに応じて変化させるように構成することができる。また、第１図形Ｇ1（ｉ）を腕図形の肩に対応する部分又はその近傍に配置し、かつ、第２図形Ｇ2（ｉ）を腕図形の手首に対応する部分又はその近傍に配置することもできる。以上のような動画を見たユーザは、ゴルフスイング中のコックを溜める動き及びコックを解放する動きの意味するところを直感的に理解することができる。
【０１１７】
図２４に示す画面Ｗ４は、「ヘッド軌道」というタブに対応する画面であり、主としてヘッド４１の挙動を表示する。画面Ｗ４には、ヘッド挙動決定工程で決定されたインパクト直前のヘッド速度Ｖhi、ヘッド軌道の角度、インパクト時のフェース角及び打点が定量的に数値で表示される。また、欄Ｃ１２には、インパクト付近の２つのタイミングにおけるヘッドを表す図形であるヘッド図形Ｇ5，Ｇ6が並べて同時に表示される。ヘッド図形Ｇ5は、インパクト時のヘッドを表しており、ヘッド図形Ｇ6は、その直前のヘッドを表している。また、欄Ｃ１２内のヘッド図形Ｇ5は、上述したヘッド軌道の角度に応じて傾けて表示される。これにより、ユーザは、アウトサイドインかインサイドアウトかといったゴルファー７のヘッド軌道を直感的に正確に理解することができる。なお、ヘッド図形Ｇ5の傾きは、インパクト時のフェース角度に応じて設定することもできる。
【０１１８】
＜３．変形例＞
以上、本発明の幾つかの実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、その趣旨を逸脱しない限りにおいて、種々の変更が可能である。例えば、以下の変更が可能である。また、以下の変形例の要旨は、適宜組み合わせることができる。
【０１１９】
＜３−１＞
上記実施形態では、コックを溜める強さである腕出力パワーＰ1，Ｐ1_AVE、及び、コックを解放する強さであるクラブ入力パワーＰ2，Ｐ2_AVEが、慣性センサユニットによる計測データに基づいて計算されたが、これらの指標の少なくとも一方を、カメラユニットによる計測データに基づいて算出することもできる。
【０１２０】
＜３−２＞
上記実施形態では、シャフトの剛性として、曲げ剛性が評価されたが、これに代えて、ねじれ剛性を評価してもよい。ねじれ剛性の値（以下、ＧＪ値）も、シャフト４０の延びる方向に沿った複数の位置において測定又は算出することができる。すなわち、シャフト４０の延びる方向に沿った複数の位置におけるねじれ剛性の分布を、シャフトの剛性としてもよい。この場合、最適剛性指標としては、ゴルファー７に適したＧＪ値（最適ＧＪ値）が決定されることになるが、最適ＧＪ値を決定するための第２スイング指標としては、最適ＧＪ値との相関が認められる任意の指標を用いることができる。このような第２スイング指標としては、例えば、特開２０１４−２１２８６２号公報に記載されているような、以下の指標を用いることができる。
【０１２１】
（１）グリップ角速度ωyが最大となるときからインパクトまでの単位時間あたりのグリップ角速度ωxの変化量の大きさ
（２）トップ付近でのグリップ角速度ωzの変化量
（３）トップからダウンスイング途中であってグリップ角速度ωyが最大となるときまでのグリップ角速度ωzの変化量の大きさ
【０１２２】
本変形例でも、第２スイング指標と最適ＧＪ値との関係を表す近似式を予め実験により算出し、記憶部２３内に格納しておくことで、計測工程で得られる計測データに基づく第２スイング指標から、最適ＧＪ値を決定することができる。
【０１２３】
また、最適剛性指標として、ゴルファー７に適したＩＦＣやフレックスに代えて又は加えて、ゴルファー７に適した調子やトルクを決定するようにしてもよい。なお、トルクは、シャフト４０全体でのねじれ剛性を表す指標である。
【０１２４】
＜３−３＞
上記実施形態では、ｉ＝１，２，・・・，Ｎについての画像Ｉ（ｉ）が、動画の形式で表示されたが、これらＮ枚の画像Ｉ（ｉ）を静止画の態様で表示することもできる。すなわち、図２３Ａ〜図２３Ｆのような画像を、同じ画面上に同時に表示することができる。このとき、時系列に沿って、画像Ｉ（ｉ）を配列することが好ましい。また、図２２に示すような１枚の画像のみを表示することもできる。図２２に示す画像では、様々なタイミングでの図形Ｇ1（ｉ）〜Ｇ4（ｉ）がストロボ画像のように１枚の画像上に同時に表示されている。従って、このような１枚の画像のみであっても、スイング動作中のコックを溜める強さ及びコックを解放する強さを直感的に理解することができる。
【符号の説明】
【０１２５】
１  慣性センサユニット
２  フィッティング装置（ＧＵＩ表示装置）
３  プログラム（ＧＵＩ表示プログラム）
４  ゴルフクラブ
５  カメラシステム
７  ゴルファー
２４Ａ  取得部
２４Ｂ  グリップ挙動導出部
２４Ｃ  肩挙動導出部
２４Ｄ  算出部
２４Ｅ  決定部
２４Ｆ  選択部
２４Ｇ  表示制御部
４０  シャフト
４１  ヘッド
４２  グリップ
Ｗ１〜Ｗ５  ＧＵＩ画面
Ｇ1（ｉ）  第１図形
Ｇ2（ｉ）  第２図形
Ｇ3（ｉ）  腕図形
Ｇ4（ｉ）  クラブ図形
Ｇ5  ヘッド図形（第１ヘッド図形）
Ｇ6  ヘッド図形（第２ヘッド図形）
【図１】
2017217423.tif 000020
【図２】
2017217423.tif 000021
【図３】
2017217423.tif 000022
【図４】
2017217423.tif 000023
【図５】
2017217423.tif 000024
【図６】
2017217423.tif 000025
【図７】
2017217423.tif 000026
【図８】
2017217423.tif 000027
【図９】
2017217423.tif 000028
【図１０】
2017217423.tif 000029
【図１１】
2017217423.tif 000030
【図１２】
2017217423.tif 000031
【図１３】
2017217423.tif 000032
【図１４】
2017217423.tif 000033
【図１５】
2017217423.tif 000034
【図１６】
2017217423.tif 000035
【図１７】
2017217423.tif 000036
【図１８】
2017217423.tif 000037
【図１９】
2017217423.tif 000038
【図２０】
2017217423.tif 000039
【図２１】
2017217423.tif 000040
【図２２】
2017217423.tif 000041
【図２３Ａ】
2017217423.tif 000042
【図２３Ｂ】
2017217423.tif 000043
【図２３Ｃ】
2017217423.tif 000044
【図２３Ｄ】
2017217423.tif 000045
【図２３Ｅ】
2017217423.tif 000046
【図２３Ｆ】
2017217423.tif 000047
【図２４】
2017217423.tif 000048