JP5825430B2 - ゴルフ用具フィッティングシステム、及びゴルフ用具フィッティングプログラム - Google Patents

Description

本発明は、ゴルフ用具フィッティングシステム、及びゴルフ用具フィッティングプログラムに関する。

２００８年、ゴルフ規則によりヘッドの反発規制が施行され、ヘッド単体で飛距離を追及することが難しくなった。これを受けて各メーカーはシャフトに着目し、シャフトのしなりを活かして飛距離アップを試みるようになった。そのため、２０１２年時点ではシャフトのバリエーションが増大し、ヘッドとの組み合わせによってクラブのバリエーションはさらに増加することとなった。この弊害として、プレイヤーがクラブ購入の際、そのプレイヤーにとってベストな用具、とりわけベストなシャフトを選び出すことが困難な状態にある。

これを解決するものがフィッティング技術である。ゴルフクラブ全体に着目したフィッティング技術として、特許文献１に記載の技術が挙げられる。これは、高速度カメラによりインパクトの瞬間のヘッド挙動を撮影し、ＤＬＴ法によって３次元座標に置き換えてヘッド姿勢を定量化したものである。これによって、クラブごとにインパクト時のヘッド姿勢が特定できることになり、適切なクラブを選択することができる。しかしながら、本手法ではプレイヤーが実際に多くのクラブでゴルフボールを打つ必要があり、ベターな選択はできてもベストな選択が困難であった。

この問題点を解決し得る技術として、特許文献２に記載の技術がある。こちらは、まず始めにスイングを計測し、そのスイングデータに基づいてヘッド挙動をシミュレートするものである。この技術を用いれば、プレイヤーは１回のスイング計測を行うだけで何百本もの試打に相当する結果が得られる。しかしながら、実際にはゴルファーは硬さや重量が異なるクラブを使用した場合、スイング自体を変化させてボールを打つ。そのため、１つのスイングデータに関して、極端に性能の異なるクラブやシャフトをシミュレートさせても、シミュレート結果が現実の解と合致しないという問題点があった。

この問題点を解決し得る技術として、特許文献３の技術が挙げられる。これは、シャフトの性能（フレックス・調子・トルクの３つの性能）ごとに変化するスイングを、応答曲面法を用いて算出し、変化後のスイングに基づいてシミュレートすることで、より現実の現象に近いシミュレートができる。

特開２００５―３１２７３４号公報 特許４８７１２１８号公報 特開２０１１−４２５号公報

しかしながら、特許文献３に記載の技術においては、例えばゴルフクラブを構成する部材の「重量」や「長さ」を変化させてシミュレートすることは難しいという問題がある。この理由は、重量や長さを変化させた場合、同一のプレイヤーであってもアドレス（ゴルフクラブを地面に付けた状態）からトップ（ゴルフクラブを振り上げた状態）、トップからインパクト（ヘッドがボールに当たった状態）まで等のスイング時間が大きく異なり、十分な計算精度が得られないためである。そのため、特許文献３に記載の技術では、変化させることができるゴルフクラブのスペックが重量や長さ等のスイング時間に大きく影響するものを除くスペックに限定されていた。

また、同文献に記載の技術においても膨大な計算結果の中から即座に適切なスペックを見つけ出すことは困難であった。そのため、即座にフィードバックが求められるフィッティングの場面に適用することは難しかった。

本発明は、このような事情に鑑みてなされたもので、スイング時間に影響するスペックを変化させた場合にも正確な計算結果を算出し、短時間でシミュレーション結果を出力することができるゴルフ用具フィッティングシステム、及びゴルフ用具フィッティングプログラムを提供する。

以下では、本発明の理解を容易にするために添付図面の参照符号を便宜的に括弧書きにて付記するが、それにより本発明が図示の態様に限定されるものではない。

上述した課題を解決するために、本発明の一態様は、スペックの異なる複数のゴルフクラブに取り付けられたセンサからスイングデータを取得するスイングデータ取得部（２１）と、前記スイングデータ取得部により取得されたスイングデータに基づき、応答曲面法によってスイング応答曲面を算出するスイング応答曲面算出部（２２１）と、前記スイングデータ取得部により取得されたスイングデータに基づき、応答曲面法によって時間応答曲面を算出する時間応答曲面算出部（２２２）と、前記スイング応答曲面と前記時間応答曲面とから計測していないゴルフクラブのスイングデータを算出し、算出したスイングデータに基づき、ゴルフスイングのシミュレーションを実行するシミュレーション実行部（２３）と、前記シミュレーション実行部によるシミュレーション結果を出力する結果出力部（２４）と、を有するゴルフ用具フィッティングシステム（２）。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくとも重量を含むことを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくともクラブ長さを含むことを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくとも２つのスペックが異なり、前記少なくとも２つの異なるスペックは互いの影響度の差が所定以内であることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、異なる３つのスペックと該各々のスペックが３水準を有する２７のスペックのうち、実験計画法におけるＬ９型直交表に基づいて選定した９つのスペックであることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記異なる３つのスペックが、重量、調子、及びフレックスであることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記異なる３つのスペックが、クラブ硬さ、クラブ長さ、クラブ総重量、及びクラブバランスのうちから選ばれることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記異なる３つのスペックが、ヘッドの重心高さ、重心深度、重心距離、及び重心角のうちから選ばれることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記異なる３つのスペックが、ヘッドのロフト角、ライ角、及びフェイス角であることを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記異なる３つのスペックのうち、第１のスペックの分割水準数をｉ、その水準をｌ、第２のスペックの分割水準数をｊ、その水準をｍ、第３のスペックの分割水準数をｋ、その水準をｎとした場合、前記結果出力部は、前記シミュレーション結果に対して（ｌ−１）ｉ＋（ｍ−１）ｊ＋（ｎ−１）ｋで定義される結果番号を割り振り、前記結果番号から（ｊ−１）×ｋ×（ｌ−１）を減じた値をプロット番号として前記結果番号と置換え、該プロット番号を横軸、シミュレーション結果を縦軸としてプロットしたデータを前記フィッティング結果として出力することを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記フィッティング結果の出力において、前記結果出力部は、下記（ｉ）〜（ｉｉｉ）の条件に当てはまる前記プロットしたデータ同士を線分で結んで出力することを特徴とする。
（ｉ）第１のスペックが第２、第３のスペックと異なり、第２、第３のスペックが同一であるデータ。
（ｉｉ）第２のスペックが第１、第３のスペックと異なり、第１、第３のスペックが同一であるデータ。
（ｉｉｉ）第３のスペックが第１、第２のスペックと異なり、第１、第２のスペックが同一であるデータ。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記結果出力部は、前記フィッティング結果を、各々の絶対値と傾きとの関連に応じて自然言語に変換して出力することを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記結果出力部は、前記シミュレーション結果の最大値を示すスペックを、製品名に変換して出力することを特徴とする。

また、本発明の一態様は、上記ゴルフ用具フィッティングシステムにおいて、前記出力部は、前記シミュレーション実行部によるシミュレーション結果に基づき、ゴルフクラブを選定して表示することを特徴とする。

また、上述した課題を解決するために、本発明の一態様は、スペックの異なる複数のゴルフクラブに取り付けられたセンサからスイングデータを取得するスイングデータ取得手順と、前記スイングデータ取得手順で取得されたスイングデータに基づき、応答曲面法によってスイング応答曲面を算出するスイング応答曲面算出手順と、前記スイングデータ取得手順で取得されたスイングデータに基づき、応答曲面法によって時間応答曲面を算出する時間応答曲面算出手順と、前記スイング応答曲面と前記時間応答曲面とから計測していないゴルフクラブのスイングデータを算出し、算出したスイングデータに基づき、ゴルフスイングのシミュレーションを実行するシミュレーション実行手順と、前記シミュレーション実行手順によるシミュレーション結果を出力する結果出力手順と、をコンピュータに実行させるためのゴルフ用具フィッティングプログラム。

本発明によれば、スイング時間に影響するスペックを変化させた場合にも正確な計算結果を算出し、短時間でシミュレーション結果を出力することができる。

本発明の第１の実施形態の構成を示すブロック図である。 ゴルフシャフトのスペックの一例を示す説明図である。 ゴルフシャフト及びゴルフクラブヘッドのスペックの一例を示す説明図である。 スイングデータを得るための試打に用いる９本のゴルフクラブのスペックを示す図である。 図１に示すゴルフ用具フィッティングシステム２の動作を示すフローチャートである。 シャフト重量が同じ場合のスイングデータを示す図である。 シャフト重量が異なる場合のスイングデータを示す図である。 スイング応答曲面のみで作成したスイングデータである。 スイング応答曲面と時間応答曲面で作成したスイングデータである。 特別な処理をしない場合のシミュレーション結果を示す図である。 シミュレーション番号とゴルフクラブのスペックとの関係を示す図である。 結果番号の割り振り後の、結果番号とゴルフクラブのスペックとの関係を示す図である。 結果番号を横軸に選んだときのシミュレーション結果を示す図である。 プロット番号を横軸に選んだときのシミュレーション結果を示す図である。 異なる９つの弾道を示す模式図である。 縦軸に複合的な条件としてＦ値を選んだシミュレーション結果を示す図である。 スペックごとの影響度を示す図である。 ３つの異なるスペックとしてヘッド重量、シャフト捩り剛性分布、ヘッド重心高を選んでシミュレーションした結果を示す図である。

以下、本発明の一実施形態によるゴルフ用具フィッティングシステムについて、図面を参照して説明する。ここではゴルフ用具のスペックとしてゴルフシャフトを例示して説明するが、これに限定されない。図１は同実施形態の構成を示すブロック図である。この図において、符号１は、予めスペックが既知であるゴルフシャフトが装着されたゴルフクラブであり、グリップ部分のシャフト内部にセンサ１１と、センサ１１の出力を無線通信によって外部へ送信する送信部１２を備えている。センサ１１はシャフト外部に取り付ける形態であってもよい。また、センサ１１は、３軸の加速度と、３軸の角速度を検出して出力する６軸センサである。さらに、センサ１１は、３軸の加速度と、３軸の角速度に加えて、地磁気を用いて３軸の方位を計測する９軸センサ等であってもよい。あるいは、センサ１１として、６軸センサと３軸の地磁気センサの２つが用られても良い。

符号２は、コンピュータ装置で構成するゴルフ用具フィッティングシステムである。符号２０は、送信部１２が送信するセンサ出力データ（スイングデータ）を受信する受信部である。符号２１は、受信部２０を介して、スイングデータを取得するスイングデータ取得部である。符号２２１は、取得したスイングデータから、ゴルフクラブ１を試打した時のスイング応答曲面を求めるスイング応答曲面算出部である。符号２２２は、取得したスイングデータから、ゴルフクラブ１を試打した時のスイング時間の応答曲面を求める時間応答曲面算出部である。符号２３は、スイング応答曲面と時間応答曲面を用いて、ＦＥＭ（Ｆｉｎｉｔｅ Ｅｌｅｍｅｎｔ Ｍｅｔｈｏｄ）によるシミュレーションを行うシミュレーション実行部である。符号２４は、シミュレーション実行部２３において実行したシミュレーションの結果を出力する結果出力部である。

次に、スイング計測データを得るための試打に用いるゴルフクラブ１について説明する。ゴルファーは、スイングを行う際に使用したシャフトの重量、曲げ剛性、捩り剛性などを考慮に入れ、そのクラブ特性にあわせたスイングを行う。このため、所定の特性を有する１本のゴルフクラブを使用してスイングの計測データを取得して、スイングの解析を行っても妥当な解析結果を得ることが難しい場合がある。すなわち、クラブのスペックが異なればスイングが異なるのである。そこで、一例として、実験計画法に基づき、数多くあるゴルフクラブのスペックのうち、３つの異なるスペックが選定され、それぞれのスペックに３水準を設けた合計２７本がクラブを用意された。その中から、Ｌ９型直交表に基づき９本のゴルフクラブが選出された。実験計画法を用いずに３要因３水準のデータを得ようとすると３^３＝２７本のクラブについて実験しなければならず、現実に実験することは難しいためである。ここでは例として、シャフトの重量、シャフトの曲げ剛性（シャフトスペックのフレックスに相当、以下フレックスという）、シャフトの曲げ剛性分布（シャフトスペックの調子に相当、以下調子という）のスペックが用いられる。用いられるクラブのスペックは、３スペックが３水準で異なるものを例示しているが、異なるスペックを複数含むゴルフクラブであればこれに限定されない。本実施形態は、特に重量や長さのスペックが異なるゴルフクラブに好適に使用できる。

図２に示すように、今回例示したシャフトの３つのスペックは、シャフトのみの重量と、たわみで規定されるフレックス（図２（ａ））と、シャフトのしなりが最も大きくなる点を示す関数Ｐの係数Ｃで規定される調子（図２（ｂ））とからなる。

図２（ａ）は、一例として、シャフトの細径側端部から９２０ｍｍの位置が下側から支持され、そこからさらに１５０ｍｍ太径側方向の位置（細径側端部から１０７０ｍｍ）が上側から支持され、細径側から１０ｍｍの位置に３．０ｋｇｆの荷重が加えられた状態を示す。このときの細径側端部の変位量がフレックス（曲げ剛性）を表す。

図２（ｂ）は、調子（曲げ剛性分布Ｌ（ｘ））の一例を示す。この図において、ＥＩ（ｉ）は曲げ剛性を示し、ｘはシャフトの先端を基準としたシャフト上の位置を示す。曲げ剛性分布Ｌ（ｘ）は、Ｌ（ｘ）＝Ｃ０×Ｐ０（ｘ）＋Ｃ１×Ｐ１（ｘ）＋Ｃ２×Ｐ２（ｘ）＋Ｃ３×Ｐ３（ｘ）…の式で表現される。

調子は、シャフトのしなりが最も大きくなる点がグリップ寄りにある元調子、ヘッド寄りにある先調子、及び中間にある中調子に分類される。

また、調子は関数Ｐの係数Ｃで規定する値によって、元調子、先調子、中調子のいずれかに割り当てられている。スイングデータの取得は、前述した９本のゴルフクラブ１を使用して行われる。

なお、他のスペックとして、シャフトの捩り剛性（シャフトスペックのトルクに相当、以下トルクという）、シャフトの捩り剛性分布、シャフト重量分布、ゴルフクラブ長さ、ヘッド重量、クラブバランス、ヘッド重心深度、ヘッド重心高、ヘッド重心距離、グリップ重量、ロフト角、ライ角及びフェイス角を適用することが可能である。

図３に示すように、トルクは、シャフトの捩り角で規定される。図３（ａ）は、一例として、シャフトの細径側端部から１０３５ｍｍの位置が固定され、シャフトの先端から４５ｍｍの位置に捩り荷重が与えられた状態を示す。なお、Ａはしなり方向が固定された状態で回転可能であり、Ｂは完全に固定されている。また、シャフト軸線上から１２０ｍｍ離れた位置に１．１５２ｋｇｆの捩り荷重が加えられている。このときのシャフト細径側端部の捩れ角度が、トルクとして定義される。

また、ヘッド重心深度は、フェイス面からヘッドの重心までの深さ（距離）で規定される（図３（ｂ））。また、ヘッド重心高は、リーディングエッジからフェイス面上の重心までの長さで規定され、重心距離は、シャフト軸線からフェイス面上の重心に向かって垂線を延ばした際の長さで規定される（図３（ｃ））。また、シャフトの捩り剛性分布とシャフトの重量分布についても、曲げ剛性分布と調子の関係と全く同一に表現することができる。

クラブバランス（スイングウェイトともいう）は、ゴルフクラブのヘッドの利き具合（スイングやワッグルの際にヘッドの重量感をどれだけ感じやすいか）を示す。クラブバランスはＫｅｎｎｅｔｈ Ｓｍｉｔｈ社製クラブバランス計「Ｇｏｌｆ Ｃｌｕｂ Ｓｃａｌｅ」を用いて測定される。

クラブ硬さは、グリップとヘッドを取り付けた状態での振動数で表される。この振動数は、藤倉ゴム工業社製「Ｇｏｌｆ Ｃｌｕｂ Ｔｉｍｉｎｇ Ｈａｒｍｏｎｉｚｅｒ」を用いて測定される。例えば、グリップ端部から１８０ｍｍが固定された状態で、クラブヘッドを振動させたときの、グリップ端部から７６０ｍｍの位置の１分間あたりの振動数がクラブ硬さとして規定される。

図４は、実験計画法Ｌ９型直交表に基づき、スイングデータを得るための試打に用いる９本のゴルフクラブのスペックを示す。シャフト重量は、正規化した値を示しており、値が小さいほど重量が重いことを意味する。また、フレックスは、正規化した値を示しており、値が小さいほど剛性が高いことを意味する。シャフト重量、フレックスは共に、規定の長さと径を有する場合にゴルフシャフトとして機能する範囲を採用した。ここでは、シャフト重量に関しては、０の場合はシャフト重量が８０ｇ、０．５の場合はシャフト重量が７０ｇ、１の場合はシャフト重量が６０ｇとした。また、フレックスに関しては、０の場合はフレックス（曲げ剛性）の値が１３０ｍｍ、０．５の場合はフレックス（曲げ剛性）の値が１８０ｍｍ、フレックスが１の場合は、フレックス（曲げ剛性）の値が２２０ｍｍとした。また、調子に関しては、０の場合は先調子、０．５の場合は中調子、１の場合は元調子とした。用いられるクラブの本数は９本に限らず、現実的に実験が行える数量であればよい。

なお、図４に示すゴルフクラブに取り付けられるヘッドは同一のものを使用しており、シャフト長さとゴルフクラブの重量も同一である。これは、フィッティングしたいスペック以外は同一のスペックのゴルフクラブを使用することで、その他の要因による変化をなくすための措置である。また、９本のゴルフクラブ１のそれぞれについて、センサ１１、送信部１２等をシャフト内に挿入することによりゴルフクラブ１の総重量が重くならないように、センサ１１、送信部１２及びこれらを動作させるために必要な機器全体の重量が２０ｇに抑えられている。これにより、市販の軽量グリップを用いることで総重量の増加を抑えることができるため、クラブの重量増加によるスイングへの悪影響が抑制される。

次に、図５を参照して、図１に示すゴルフ用具フィッティングシステム２の動作を説明する。まず、試打を行う者（フィッティング対象のユーザ）は、９本のゴルフクラブ１のうち、クラブ番号が「１」であるゴルフクラブ１を使用して試打を行う。センサ１１は、この試打動作中の検出結果（スイングデータ）を送信部１２に対して出力する。送信部１２は、センサ１１からスイングデータの出力が行われると、無線通信を使用してスイングデータを外部へ送信する。このスイングデータは、受信部２０によって受信され、スイングデータ取得部２１に出力される。スイングデータ取得部２１は、このスイングデータを内部に保持する（ステップＳ１）。

この動作を他のゴルフクラブ１（クラブ番号が「２」〜「９」）について繰り返し行うことによって、スイングデータ取得部２１の内部には、少なくとも９本のゴルフクラブ１のそれぞれによって１回試打した場合の時系列のスイングデータが保持されることになる。スイングデータ取得部２１は、取得したスイングデータをゴルフクラブ１のグリップ部分の移動速度データと、シャフトの軸回転データとに変換する。この変換は、ゴルフクラブ１内に挿入されているセンサ１１の取付位置と、ゴルフクラブ１のグリップ部分の予め決められた２点の位置関係とから幾何学的に変換することにより行われる。

本実施形態では、スイングデータ取得部２１は、スイングデータを参照して、スイングデータの中からスイングのトップ（ゴルフクラブを振り上げた状態）からインパクト（ヘッドがボールに当たった状態）までのスイングデータのみを移動速度データと軸回転データに変換する対象とした。しかし、アドレス（ゴルフクラブを地面に付けた状態）からインパクトまでのスイングデータを解析対象としてもよい。

次に、スイング応答曲面算出部２２１は、スイングデータ取得部２１に保持されている９本のクラブのスイングデータを読み出して、試打者の技量と癖を１次関数化したスイング応答曲面を算出する（ステップＳ２）。スイング応答曲面とは、図４に示す９種類のゴルフクラブで試打した場合に得られたゴルフクラブのグリップ部分の移動速度データと軸回転データと、ゴルフクラブの３つのスペック（ここではシャフト重量、フレックス、調子）の関係式である。

ここで数式１に関して説明する。複数のゴルフクラブを試打することにより計測されたスイングデータをｆ_１〜ｆ_９で表現し、時間ｔはｔ＝｛ｔ_１，…，ｔ_ｎ｝に離散化する。なお、本実施形態では９本のゴルフクラブで試打したため、スイングデータはｆ_１〜ｆ_９までだが、その数値はゴルフクラブの試打本数によって変わる。さらに、ｆ_ｊ（ｔ_ｉ）はｊ番目のゴルフクラブによる計測量で、具体的には３方向加速度｛ａ_ｘ，ａ_ｙ，ａ_ｚ｝、３方向角速度｛ω_ｘ，ω_ｙ，ω_ｚ｝のそれぞれの量を示す。

９本のゴルフクラブのそれぞれの３つのスペック（設計変数）を｛ｘ_ｊ，ｙ_ｊ，ｚ_ｊ｝（ｊ＝１…９）とすると、数式１の関係が得られる。そして、各ｔ_ｉに対して、数式１を解いていく。

ここで、ｘ，ｙ，ｚはそれぞれ設計変数であり、ｘが第１のスペック（シャフト重量）、ｙが第２のスペック（フレックス）、ｚが第３のスペック（調子）である。ｘ_１〜ｘ_ｎ、ｙ_１〜ｙ_ｎ、ｚ_１〜ｚ_ｎにおける１〜ｎの番号は各クラブの番号に対応する。ｘバー，ｙバー，ｚバーは各解析に都合の良い任意値が与えられる（例えば、設計変数の中間値）。

数式１を解くことにより、数式２のように応答曲面の係数ａ_１〜ａ_４を求めることができる。数式１はｎ＝５以上の場合、一般には過条件にあり厳密な解は存在しない。そこで一般化逆行列Ａ＋（ムーア・ペンローズ逆行列、擬似逆行列とも言う）が用いられる。これは、厳密な解が存在しない場合に近似解を求める手法である。すなわち、誤差｜Ａｘ−ｂ｜^２を最小とする解を求める手法である。一般的な数学的手法であるため、詳細は省略する。なお、ここでは数式１を解くために、ＭａｔｈＷｏｒｋｓ社製数値計算ソフト「ＭＡＴＬＡＢ」が用いられた。

数式２で求められた係数ａ_１〜ａ_４が試打者の技量とスイングの癖に相当する値である。つまり、この処理によって、ｘ、ｙ、ｚを実際には計測していないスペックに変更したとしても、数式３に示したように関数ｆで表されるスイングデータが求められることになる。換言すると、数式３によって任意の｛ｘ，ｙ，ｚ｝に対する３方向加速度、３方向角速度の近似値として下記が得られる。数式３から、計測していない任意のｍ番目のシャフトによるスイングデータは数式４で示すように表わされる。数式４で表されるｆ_ｍ（ｔ）が、スイング応答曲面である。

ゴルフクラブによるスイングの時間が変化する場合は、各ゴルフクラブによるスイングの時間をｔ_ｋ（ｋ＝１，…，９）とすると、各クラブｋに対して数式５のようにサンプリング時間が正規化される。

スイング時間に大きな差がない場合は、正規化の手間を省くために、９本のクラブのうち最もスイング時間の短かったｔ_ｍｉｎで打ち切ってもよい。

次に、時間応答曲面算出部２２２は、スイングデータ取得部２１に保持されている９本のクラブのスイングデータを読み出して、試打者のスイング時間を定式化した時間応答曲面を算出する（ステップＳ３）。時間応答曲面とは、図４に示す９種類のゴルフクラブで試打した場合に得られたスイング時間と、ゴルフクラブの３つのスペック（ここではシャフト重量、フレックス、調子）の関係式である。時間応答曲面は、数式６〜数式８によって算出することができる。

数式６におけるｇ_１〜ｇ_９は、各ゴルフクラブのスイング時間である。そして、数式７で求められた係数ｂ_１〜ｂ_４が試打者のスイング時間に相当する値である。この処理によって、重量の違いによるスイング時間の違いを計算することができる。

次に、数式８に基づき、それぞれのｔがｇ_ｍに変換される。数式８で表されるｇ_ｍが、時間応答曲面である。また、数式９から導きだされるｆ_ｍ’（ｔ）が、スイング応答曲面と時間応答曲面に基づき新たに算出したスイングデータである。

シミュレーション実行部２３は、スイングデータ取得部２１に保存されている計測済みのグリップ部分の移動速度データ及び軸回転データと、スイング応答曲面算出部２２１が算出したスイング応答曲面と時間応答曲面算出部２２２が算出した時間応答曲面とから、計測していないゴルフクラブのスイングデータを算出する。そして、シミュレーション実行部２３は、算出されたスイングデータに基づきゴルフクラブヘッドの運動を動的有限要素法によりシミュレーションを行う（ステップＳ４）。

この解析によって得られるゴルフクラブの運動を示すシミュレーション結果は、インパクト時のクラブスピード、フェイスアングル（クラブフェイスの飛球線に対する水平方向角度）、インパクトロフト（インパクト時における地面に対してのロフト角）である。クラブスピードは、球の飛距離に対応する。フェイスアングルは、球の方向性に対応する。インパクトロフトは、弾道高さに対応する。また、上記以外にも得られるシミュレーション結果はいくつも挙げられるが、一例を示すと、クラブパス、アタックアングル、ボールスピード、飛距離等あらゆる値をとり得る。

この動的有限要素法によるゴルフクラブヘッドの運動の解析方法は、公知の方法（例えば、市販の有限要素法ソフトウェア）を用いるため、詳細な処理動作の説明を省略する。シミュレーション実行部２３は、考えられるスペック全てのゴルフクラブを使用して、試打者の技量、癖を１次関数化したスイング応答曲面と、試打者のスイング時間を定式化した時間応答曲面とを入力としてゴルフクラブヘッドの運動を解析する。このため、試打者の技量、癖及びスイング時間（スペックが異なる場合のスイングの性質とスイング時間の変化）を考慮したゴルフクラブの運動をシミュレートすることができる。

前述の通り、重量やクラブ長さはスイング時間に大きな影響を及ぼす。図６、図７にシャフト重量が同じ場合（図６）と異なる場合（図７）のスイングデータ（一例としてω_ｘ）を示す。シャフト重量が同じ場合は、図６に示すように、スイング時間のばらつきは少ない。一方、シャフト重量がそれぞれ、４０ｇ、６０ｇ、８０ｇと異なる場合は、図７に示すようにスイング時間が大きくばらつく。

この場合、前述の通り数式２を一般化逆行列を用いて解くには、ｔを一定にする必要がある。図６のようにスイング時間が大きく変わらない場合はどのような方法でスイング時間を一定にしても大きな誤差は生まれない。しかし、図７のようにスイング時間が大きく異なる場合に無理やり一定にすると大きな誤差が生まれる。そこで、まずはそれぞれのｔが標準化された後、時間応答曲面算出部２２２は時間応答曲面を求める。求められた時間応答曲面に応じて再度ｔを変換することによって、スイング時間を反映したスイングデータを求めることができる。

このように、重量や長さを変更した場合にも適切なスイング時間を反映したスイングデータを作成することができる。そのため、重量や長さ等のスイング時間に大きく影響するスペックにおいても正確なシミュレーション結果を得ることができる。また、スイング時間に影響が少ないスペックにおいても厳密に言えばスイング時間は異なる。スイング応答曲面と時間応答曲面を導入することで、スイング時間に影響が少ないスペックにおいても計算精度を向上させることができる。図８がスイング応答曲面のみを用いて算出したデータ、図９がスイング応答曲面と時間応答曲面を用いて算出したデータである。

次に、結果出力部２４が出力する結果データの生成処理動作について説明する。ここでは、算出すべきゴルフシャフトのスペックであるシャフト重量の水準（第１のスペックの分割水準数ｉ）が３、フレックスの水準（第２のスペックの分割水準数ｊ）が５、調子の水準（第３のスペックの分割水準数ｋ）が５であるものとして説明する。分割水準数は、実際に試打した９本のスペックを任意の数に分割したものである。例えば、フレックスを例に挙げると、実際に試打した９本のうち、フレックスにはＸ，Ｒ，Ｌが存在するが、フレックスの分割水準数ｊを５とすることで、前記３つのフレックスの間のフレックス、すなわちＳ，Ａについても結果を算出することができる（フレックスは硬い方からＸ，Ｓ，Ｒ，Ａ，Ｌで表される）。

この場合、シミュレーション実行部２３においては、３×５×５＝７５本のシャフトを装着したゴルフクラブの運動を模擬したシミュレーション結果が得られることになる。この結果データを、例えばクラブスピードを縦軸、シミュレーション番号を横軸として画面表示の出力が行われた場合、図１０に示すような表示となる。これでは、どの値がどのようなスペック値を有するシャフトを示しているかを把握することが難しい。なお、図１１はシミュレーション番号とゴルフクラブのスペックとの関係を示している。

そこで、次のような処理を行うことによって、ユーザがゴルフシャフトのスペック値を読み取ることができる。

まず、シミュレーション番号に下記の定義に基づき結果番号を割り振る。
結果番号：第１のスペックの分割水準数をｉ、その水準をｌ、第２のスペックの分割水準数をｊ、その水準をｍ、第３のスペックの分割水準数をｋ、その水準をｎとし、シミュレーション結果に対して（ｌ−１）ｉ＋（ｍ−１）ｊ＋（ｎ−１）ｋで定義される。換言すると下記の通りである。

［１］第１のスペックと第２のスペックの水準を最小値とし、第３のスペックの水準を最小値〜最大値まで変化させる。
［２］第１のスペックの水準を最小値とし、第２のスペックの水準を１水準大きくし、第３のスペックの水準を最小値〜最大値まで変化させる。
［３］第２スペックの水準が最大値となるまで［２］を繰り返す。
［４］第１スペックの水準を１水準大きくし、［２］、［３］を実行する。
［５］第１スペックの水準が最大値となるまで［４］を繰り返す。
上記の手順をより詳細に説明すると下記の通りである。
（Ａ）シャフト重量を１水準目で固定、フレックスを１水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｂ）シャフト重量を１水準目で固定、フレックスを２水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｃ）シャフト重量を１水準目で固定、フレックスを３水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｄ）シャフト重量を１水準目で固定、フレックスを４水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｅ）シャフト重量を１水準目で固定、フレックスを５水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｆ）シャフト重量を２水準目で固定、フレックスを１水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｇ）シャフト重量を２水準目で固定、フレックスを２水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｈ）シャフト重量を２水準目で固定、フレックスを３水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｉ）シャフト重量を２水準目で固定、フレックスを４水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｊ）シャフト重量を２水準目で固定、フレックスを５水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｋ）シャフト重量を３水準目で固定、フレックスを１水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｌ）シャフト重量を３水準目で固定、フレックスを２水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｍ）シャフト重量を３水準目で固定、フレックスを３水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｎ）シャフト重量を３水準目で固定、フレックスを４水準目で固定、調子を先調子から元調子に向かって５水準で変える。
（Ｏ）シャフト重量を３水準目で固定、フレックスを５水準目で固定、調子を先調子から元調子に向かって５水準で変える。

前述の（Ａ）〜（Ｏ）の順序で結果番号を定義すると、結果番号は図１２に示すように定義される。図１２において、シャフト重量０は、シャフトが重い（例えば、８０ｇ）ことを示し、シャフト重量０．５は、シャフトの重さが中程度である（例えば、７０ｇ）ことを示し、シャフト重量１は、シャフトが軽い（例えば、６０ｇ）ことを示している。また、フレックス０は、最も硬い（例えば、フレックスＸ）ことを示し、フレックス０．２５（例えば、フレックスＳ）、０．５（例えば、フレックスＲ）、０．７５（例えば、フレックスＡ）へと値が増加するにつれて徐々に柔らかくなる。フレックス１は最も柔らかい（例えば、フレックスＬ）ことを示している。また、調子０は、先調子を示し、調子０．２５は、先中調子を示し、調子０．５は、中調子を示し、調子０．７５は、中元調子を示し、調子１は、元調子を示す。

これに基づき、縦軸をクラブスピード、横軸を結果番号としてプロットした場合、図１３に示すような結果が得られる。シミュレーション実行部２３が実行順序を決めてシミュレーションを行い、その実行順序に基づき結果出力部２４が結果出力を行うことにより、図１３に示すようなクラブスピードの傾向を把握することができるようになる。しかし、図１３に示す出力形式ではまだ不十分である。

結果出力部２４は、結果番号から（ｊ−１）×ｋ×（ｌ−１）を減じた値をプロット番号として、再度プロットする。この方法でデータをプロットして出力したフィッティング結果が図１４となる。

さらに詳細に説明すると、結果出力部２４は、次の（ａ）〜（ｅ）ように結果番号を変換して新たにプロット番号を付与し、そのプロット番号とシミュレーション結果を再度プロットして出力する。これによって、ゴルフ用具のスペック値とシミュレーション結果との関係の傾向が把握できるようになる。

（ａ）シミュレーション順序にある（Ａ）から（Ｅ）によって求めた結果をプロットする。
（ｂ）シミュレーション順序にある（Ｆ）から（Ｊ）によって求めた結果を左にずらしてプロットする。ずらす量はｘ軸２０個分（＝７５／３−５＝２０）である。
（ｃ）シミュレーション順序にある（Ｋ）から（Ｏ）によって求めた結果を（ｂ）でプロットしたものに対して、上記と同様の量だけさらにずらしてプロットする。
（ｄ）調子だけが異なるもの同士、フレックスだけが異なるもの同士、シャフト重量だけが異なるもの同士を直線で結ぶ。

なお、煩雑さを防ぐため、調子は全て結び、フレックスとシャフト重量は端のスペックのみ結んだ。また、ここでは実際のシミュレーション結果値を直線で結んだが、直線で結ばれている点を４次関数で近似し、その近似値をシミュレーション結果値として用いることもできる。このように４次関数近似を用いることでゴルフ初心者にありがちなスイングの乱れを吸収し、より意義のある結果を表示することができる。

図１４において、縦軸は、クラブスピードであり、値が大きいほどスピードが速い（飛距離が出る）、値が小さいほどスピードが遅い（飛距離が出ない）ことを示す。横軸は、新たに設定したプロット番号である。○印がスペックと関連付けられたゴルフクラブを示す。各値にはゴルフシャフトのスペックである「フレックス」、「調子」、「シャフト重量」の値が関係づけられている。

図１４に示す結果の出力例を見ると、クラブスピードを最大にするためには、シャフト重量は「６０ｇ」、フレックスは一番硬い「Ｘ」、調子は「中元調子」を選択すると良いことが分かる。また、ここでは一例としてクラブスピードを示したが、フェイスアングルやインパクトロフト等も同様の手法で表現することができる。

さらに、結果出力部２４は、スペックの傾向を自然言語によって表現して表示するようにしてもよい。すなわち、結果出力部２４は、フィッティング結果を、各々の絶対値と傾きとの関連に応じて自然言語に変換して表示する。例えば、フィッティングの目的がクラブスピードを上げることであれば、結果出力部２４は、「フレックスが硬いほどクラブスピードが上がる」、「シャフトが軽いほどクラブスピードが上がる」、「調子が元調子であるほどクラブスピードが上がる」、「シャフト重量の影響が最も大きい」等の表示を行う。フィッティングの目的がインパクトロフトを大きくして球を上げることであれば、結果出力部２４は、「フレックスが硬いほど球が上がる」、「シャフトが軽いほど球が上がる」、「調子が先調子であるほど球が上がる」、「フレックスの影響が最も大きい」等の表示を行う。これにより、ユーザは、スペックの傾向を自然言語で把握することができるようになるため、容易にゴルフ用具の選択ができるようになる。

このように、フィッティングの目的である「クラブスピード」とスペックとを関係つけて結果を出力するようにしたため、一目でスペックとの関連付けができる。また、結果出力部２４は、クラブスピードが最大のときのゴルフシャフトのスペックを１点だけ出力するのではなく、ありうるスペックすべてを出力するようにした。これによって、ほどよい飛距離、ほどよい弾道高さ、ほどよい弾道の曲がりを達成できるスペックを把握することができるようになり、ユーザの好みに合わせてゴルフシャフトのスペックを選択することができるようになる。これは特に、フェイスアングルを算出した際に有効であり、ほどよい弾道の曲がりが選べることはユーザにとって好ましい。弾道を曲げたいがために、最大に曲がるシャフトを選んでしまうと、弾道の曲がりすぎの問題が出てくるためである。

なお、図１４において、結果出力部２４が○印に対してスペックを表す文字情報を表示する例を示したが、色分けをして関係づけたスペックを表示するようにしてもよい。また、さらに見やすくするために、結果出力部２４は、ベストな値を含むものを目立つように、線の色の濃さ等を変更して出力するようにしてもよい。また、結果出力部２４は、図１４をさらに簡略化したデータとして自然言語に変わるユーザへの伝達手段を出力してもよい。

また、前述した説明では、第１のスペックをシャフト重量（３水準）、第２のスペックをフレックス（５水準）、第３のスペックを調子（５水準）として説明したが、他のスペックを用いてもよい。

例えば、第１のスペックをトルク（５水準）、第２のスペックをフレックス（５水準）、第３のスペックを調子（５水準）としてもよい。また、第１のスペックをフレックス（５水準）、第２のスペックをトルク（５水準）、第３のスペックを重量（３水準）としてもよい。

次に、フィッティング結果として、クラブスピードと弾道高さ、クラブスピードと弾道の曲がりなど複合的な条件を出力する動作について説明する。一般に、ゴルファーはゴルフクラブに１つの性能だけを求めるのではなく、複合的な性能を求める。例えば、クラブスピードを最大にしつつスライスを防ぐ、同時に弾道を高くするなどである。複合的な条件を出力する場合、ユーザの要望に基づいて図１５に示される９つの弾道のうち１つが選択される。球の弾道は、弾道の高さと弾道の方向によって決まる。弾道高さを「Ｈｉｇｈ（高弾道）」、「Ｍｉｄ（中弾道）」、「Ｌｏｗ（低弾道）」の３つに分け、弾道方向を「Ｆａｄｅ（フェード）」、「Ｓｔｒａｉｇｈｔ（ストレート）」、「Ｄｒａｗ（ドロー）」の３つに分けると、図１５に示すように９つの弾道に分けることができる。ここで、フェードとは、右打ちの場合に右へ曲がる弾道であり、ドローとは、右打ちの場合に左へ曲がる弾道のことである。そして、例えば、クラブスピードを最大にしつつ、選ばれた弾道を満たすような条件が出力される。複合的な条件を出力するため、目的関数Ｆを最大にするスペックが求められる。目的関数Ｆは、Ｆ＝α×ｆ_１＋β×ｆ_２＋γ×ｆ_３で表される。

上記の式において、ｆ_１は、第１の結果データ（例えばクラブスピード）であり、ｆ_２は第２の結果データ（例えば、フェイスアングル）であり、ｆ_３は、第３の結果データ（例えば、インパクトロフト）であり、α、β、γは重み係数である。α、β、γは、各ゴルファーの要望に応じて適宜選択される。概ね、αは、（β＋γ）の値の１〜３倍であることが好ましい。これはゴルファーにとって、クラブスピードが最も重要であるためである。

図１６は、複合的な条件として、弾道を低弾道・ドロー（Ｌｏｗ・Ｄｒａｗ）、クラブスピードを最大とした場合のフィッティング結果の出力例を示す図である。ここでは縦軸がクラブスピードを高めつつ、ＬｏｗＤｒａｗの弾道を打つために適しているかを示している。この図から、このプレイヤーにとっての最適なスペックは、シャフト重量が６０ｇ、フレックスがＳ、調子が中調子であることがわかる。

なお、結果出力部２４は、シミュレーション結果（クラブスピード、インパクトロフト、フェイスアングル等）の最大値を示すスペックを、ゴルフ用具の製品名に変換して表示するようにしてもよい。

このように、弾道を選択し、その選択した弾道を実現するために適しているか否かを表示することができるため、ユーザは、所望の弾道を実現するゴルフ用具を容易に選ぶことができるようになる。すなわち、前述したゴルフ用具フィッティングシステムによってゴルフ用具のフィッティングが可能となる。また、短時間で計測・解析・結果表示を行うことができ、最適なゴルフシャフトのスペックを視覚的に判断できる。

前述の例では、シャフト重量・フレックス・調子を選んだがこれに限定されない。また、３スペックに限らずこれより多くても少なくても良い。上述の通り、スイング時間が変わりやすいスペック（例えば重量やクラブ長さ）が異なるものが好ましい。他にはシャフトの曲げ剛性、シャフトの捩り剛性、シャフト重量、シャフトの曲げ剛性分布、シャフトの捩り剛性分布、シャフト重量分布、ゴルフクラブ長さ、ヘッド重量、クラブバランス、ヘッド重心深度、ヘッド重心高、ヘッド重心距離、グリップ重量、ロフト角、ライ角及びフェイス角のようなスペックが考えられる。

また、時間応答曲面としては、全てのスイング時間を同等とするように、各スイングデータを最も長いものに引き伸ばしたものを用いてもよい。この場合、ある程度の精度と計算時間の短縮を図ることができる。

また、上述同様複数の異なるスペックを選ぶ場合、なるべく影響度が近いものを組み合わせる方がよい。影響度とは、そのスペックを変更した場合に、最終的なヘッド挙動にどの程度影響するかを示したものである。影響度を図１７に示す。この理由としては次のとおりである。

図１８は３つの異なるスペックとしてヘッド重量（影響度：５）、シャフト捩り剛性分布（影響度：１）、ヘッド重心高（影響度：１）を選んでシミュレーションした結果を示す図である。このように、影響度が大きく異なるものを選んだ場合、影響度の大きいものにデータが引っ張られ、残りの影響度が小さいスペックの結果がほとんどわからない状況となる。そのため、異なるスペックとしては影響度が所定以内（本実施形態では２以内）のスペックである方が好ましい。また、３つの異なるスペックの影響度は、それぞれ同じであることがより好ましい。図１８の例では影響度の大きいヘッド重量の違いによってヘッドスピードが変わることがわかるが、他の２つのスペックでは差が見られない。

３つのスペックを選択して組み合わせる場合、シャフト重量、曲げ剛性（フレックス）、曲げ剛性分布（調子）を用いるのがよい。この場合、シャフトのフィッティング及び設計として好適に使用できる。

また、曲げ剛性（フレックス）、クラブ長さ、ヘッド重量を組み合わせてもよい。この場合、クラブのフィッティング及び設計として好適に使用できる。

また、ヘッドの重心高さ、重心深度、重心距離、及び重心角を組み合わせてもよい。この場合、ヘッドのフィッティング及び設計として好適に使用できる。

さらに、ヘッドのロフト角、ライ角、フェイス角を組み合わせてもよい。この場合、ヘッドのフィッティング及び設計として好適に使用できる。また、クラブ硬さ、クラブ長さ、クラブ総重量、及びクラブバランスを組み合わせてもよい。

ここでは上記の通り、好適な組み合わせを例示したが、これに限定されるものではない。

なお、ゴルフ用具フィッティングシステムの各機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより上述の各部の処理を行ってもよい。ここで、「記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行する」とは、コンピュータシステムにプログラムをインストールすることを含む。ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータシステム」は、インターネットやＷＡＮ、ＬＡＮ、専用回線等の通信回線を含むネットワークを介して接続された複数のコンピュータ装置を含んでもよい。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。このように、プログラムを記憶した記録媒体は、ＣＤ−ＲＯＭ等の非一過性の記録媒体であってもよい。また、記録媒体には、当該プログラムを配信するために配信サーバからアクセス可能な内部または外部に設けられた記録媒体も含まれる。配信サーバの記録媒体に記憶されるプログラムのコードは、端末装置で実行可能な形式のプログラムのコードと異なるものでもよい。すなわち、配信サーバからダウンロードされて端末装置で実行可能な形でインストールができるものであれば、配信サーバで記憶される形式は問わない。なお、プログラムを複数に分割し、それぞれ異なるタイミングでダウンロードした後に端末装置で合体される構成や、分割されたプログラムのそれぞれを配信する配信サーバが異なっていてもよい。さらに「コンピュータ読み取り可能な記録媒体」とは、ネットワークを介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。また、上記プログラムは、上述した機能の一部を実現するためのものであってもよい。さらに、上述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。

また、上述した機能の一部または全部を、ＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｉｏｎ）等の集積回路として実現してもよい。上述した各機能は個別にプロセッサ化してもよいし、一部、または全部を集積してプロセッサ化してもよい。また、集積回路化の手法はＬＳＩに限らず専用回路、または汎用プロセッサで実現してもよい。また、半導体技術の進歩によりＬＳＩに代替する集積回路化の技術が出現した場合、当該技術による集積回路を用いてもよい。

以上のように、本発明にかかるゴルフ用具フィッティングシステム、及びゴルフ用具フィッティングプログラムは、ゴルファーの要望や癖、技量に応じたゴルフ用具のフィッティングに適用できる。

１・・・ゴルフクラブ、１１・・・センサ、１２・・・送信部、２・・・ゴルフシャフトフィッティングシステム、２０・・・受信部、２１・・・スイングデータ取得部、２２１・・・スイング応答曲面算出部、２２２・・・時間応答曲面算出部、２３・・・シミュレーション実行部、２４・・・結果出力部

Claims (15)

1. スペックの異なる複数のゴルフクラブに取り付けられたセンサからスイングデータを取得するスイングデータ取得部と、
   前記スイングデータ取得部により取得されたスイングデータに基づき、応答曲面法によってスイング応答曲面を算出するスイング応答曲面算出部と、
   前記スイングデータ取得部により取得されたスイングデータに基づき、応答曲面法によって時間応答曲面を算出する時間応答曲面算出部と、
   前記スイング応答曲面と前記時間応答曲面とから計測していないゴルフクラブのスイングデータを算出し、算出したスイングデータに基づき、ゴルフスイングのシミュレーションを実行するシミュレーション実行部と、
   前記シミュレーション実行部によるシミュレーション結果を出力する結果出力部と、
   を有するゴルフ用具フィッティングシステム。
2. 前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくともシャフト重量を含む請求項１記載のゴルフ用具フィッティングシステム。
3. 前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくともクラブ長さを含む請求項１または２記載のゴルフ用具フィッティングシステム。
4. 前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、少なくとも２つのスペックが異なり、前記少なくとも２つの異なるスペックは互いの影響度の差が所定以内である請求項１から３のいずれか１項に記載のゴルフ用具フィッティングシステム。
5. 前記シミュレーションに用いられる前記複数のゴルフクラブのスペックは、異なる３つのスペックと該各々のスペックが３水準を有する２７のスペックのうち、実験計画法におけるＬ９型直交表に基づいて選定した９つのスペックである請求項１から４のいずれか１項に記載のゴルフ用具フィッティングシステム。
6. 前記異なる３つのスペックが、シャフト重量、調子、及びフレックスである請求項５に記載のゴルフ用具フィッティングシステム。
7. 前記異なる３つのスペックが、クラブ硬さ、クラブ長さ、クラブ総重量、及びクラブバランスのうちから選ばれる請求項１から４のいずれか１項に記載のゴルフ用具フィッティングシステム。
8. 前記異なる３つのスペックが、ヘッドの重心高さ、重心深度、重心距離、及び重心角のうちから選ばれる請求項１から４のいずれか１項に記載のゴルフ用具フィッティングシステム。
9. 前記異なる３つのスペックが、ヘッドのロフト角、ライ角、及びフェイス角である請求項１から４のいずれか１項に記載のゴルフ用具フィッティングシステム。
10. 前記異なる３つのスペックのうち、第１のスペックの分割水準数をｉ、その水準をｌ、第２のスペックの分割水準数をｊ、その水準をｍ、第３のスペックの分割水準数をｋ、その水準をｎとした場合、
    前記結果出力部は、
    前記シミュレーション結果に対して（ｌ−１）ｉ＋（ｍ−１）ｊ＋（ｎ−１）ｋで定義される結果番号を割り振り、
    前記結果番号から（ｊ−１）×ｋ×（ｌ−１）を減じた値をプロット番号として前記結果番号と置換え、該プロット番号を横軸、シミュレーション結果を縦軸としてプロットしたデータを前記フィッティング結果として出力する請求項５〜９のいずれか１項に記載のゴルフ用具フィッティングシステム。
11. 前記フィッティング結果の出力において、前記結果出力部は、下記（i）〜（iii）の条件に当てはまる前記プロットしたデータ同士を線分で結んで出力する請求項１０に記載のゴルフ用具フィッティングシステム。
    （i）第１のスペックが第２、第３のスペックと異なり、第２、第３のスペックが同一であるデータ。
    （ii）第２のスペックが第１、第３のスペックと異なり、第１、第３のスペックが同一であるデータ。
    （iii）第３のスペックが第１、第２のスペックと異なり、第１、第２のスペックが同一であるデータ。
12. 前記結果出力部は、前記フィッティング結果を、各々の絶対値と傾きとの関連に応じて自然言語に変換して出力する請求項１０または１１に記載のゴルフ用具フィッティングシステム。
13. 前記結果出力部は、前記シミュレーション結果の最大値を示すスペックを、製品名に変換して出力する請求項１から４のいずれか１項に記載のゴルフ用具フィッティングシステム。
14. 前記出力部は、前記シミュレーション実行部によるシミュレーション結果に基づき、ゴルフクラブを選定して表示する請求項１に記載のゴルフ用具フィッティングシステム。
15. スペックの異なる複数のゴルフクラブに取り付けられたセンサからスイングデータを取得するスイングデータ取得手順と、
    前記スイングデータ取得手順で取得されたスイングデータに基づき、応答曲面法によってスイング応答曲面を算出するスイング応答曲面算出手順と、
    前記スイングデータ取得手順で取得されたスイングデータに基づき、応答曲面法によって時間応答曲面を算出する時間応答曲面算出手順と、
    前記スイング応答曲面と前記時間応答曲面とから計測していないゴルフクラブのスイングデータを算出し、算出したスイングデータに基づき、ゴルフスイングのシミュレーションを実行するシミュレーション実行手順と、
    前記シミュレーション実行手順によるシミュレーション結果を出力する結果出力手順と、
    をコンピュータに実行させるためのゴルフ用具フィッティングプログラム。

Images