"実測又はシミュレーションによってゴルフクラブの打球音を取得する取得ステップと、上記打球音を解析する解析ステップと、上記解析結果に基づいて、打球音高さを評価する評価ステップとを含み、上記解析ステップが、音圧特性を取得するステップを含み、上記評価ステップが、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップを含む打球音高さの評価方法。上記音圧系評価ステップでは、上記音圧特性が大きいほど打球音高さが高いと評価される請求項1に記載の評価方法。上記解析ステップが、音圧特性を取得するステップと、周波数特性を取得するステップとを含み、上記評価ステップが、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップと、上記周波数特性に基づいて打球音高さを評価する周波数系評価ステップとを含む請求項1又は2に記載の評価方法。周波数特性よりも音圧特性に依存して上記打球音の高さを評価する音圧派被験者と、音圧特性よりも周波数特性に依存して上記打球音の高さを評価する周波数派被験者とを分類する分類ステップと、上記音圧派被験者による評価結果を得るステップと、上記周波数派被験者による評価結果を得るステップとを更に含み、上記音圧系評価ステップでは、上記音圧派被験者による評価結果が用いられ、上記周波数系評価ステップでは、上記周波数派被験者による評価結果が用いられる請求項3に記載の評価方法。上記周波数特性が、1次ピーク周波数又は1次固有振動数である請求項3又は4に記載の評価方法。上記音圧特性が、オーバーオール音圧レベルである請求項1から5のいずれかに記載の評価方法。音圧特性が大きいほど打球音高さが高いことを考慮した設計を含むゴルフクラブヘッドの設計方法。上記従来技術とは異なる技術思想により、より精度よく打球音の高さが評価できることが判明した。本発明の目的は、ゴルファーによる評価に近く高精度な打球音高さの評価方法を提供することにある。本発明に係る打球音高さの評価方法は、実測又はシミュレーションによってゴルフクラブの打球音を取得する取得ステップと、上記打球音を解析する解析ステップと、上記解析結果に基づいて、打球音高さを評価する評価ステップとを含む。上記解析ステップが、音圧特性を取得するステップを含む。上記評価ステップが、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップを含む。好ましくは、上記音圧系評価ステップでは、上記音圧特性が大きいほど打球音高さが高いと評価される。好ましくは、上記解析ステップが、音圧特性を取得するステップと、周波数特性を取得するステップとを含む。好ましくは、上記評価ステップが、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップと、上記周波数特性に基づいて打球音高さを評価する周波数系評価ステップとを含む。好ましくは、上記評価方法は、周波数特性よりも音圧特性に依存して上記打球音の高さを評価する音圧派被験者と、音圧特性よりも周波数特性に依存して上記打球音の高さを評価する周波数派被験者とを分類する分類ステップを更に含む。好ましくは、上記評価方法は、上記音圧派被験者による評価結果を得るステップを更に含む。好ましくは、上記評価方法は、上記周波数派被験者による評価結果を得るステップを更に含む。好ましくは、上記音圧系評価ステップでは、上記音圧派被験者による評価結果が用いられる。好ましくは、上記周波数系評価ステップでは、上記周波数派被験者による評価結果が用いられる。好ましくは、上記周波数特性が、1次ピーク周波数又は1次固有振動数である。好ましくは、上記音圧特性が、オーバーオール音圧レベルである。好ましいゴルフクラブヘッドの設計方法は、音圧特性が大きいほど打球音高さが高いことを考慮した設計を含む。ゴルファーによる評価に近く高精度な打球音高さの評価が達成されうる。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。従来、打球音高さの評価は、周波数特性に基づいてなされていた。例えば前述の特開2008-132262号公報では、最大ピーク周波数又はエネルギー加重平均周波数が用いられていた。音の高さを評価する以上、周波数特性を用いるのが当然であると考えられてきた。ところが、本発明者は、音圧特性が音の高さに相関しうることを見いだした。音の高さは、人間の感覚である。人間の聴覚特性により、何らかの理由で、音圧特性が音の高さに相関しうるものと考えられるが、この理由は不明である。更に、打球音の高さ評価に関して、ゴルファー被験者が2つのグループに分類されうることが判明した。このグループは、次のグループA及びグループBである。[グループA]:音圧特性よりも周波数特性に依存して上記打球音の高さを評価する周波数派被験者[グループB]:周波数特性よりも音圧特性に依存して上記打球音の高さを評価する音圧派被験者グループBの音圧派被験者は、音圧特性が大きいほど打球音が高いと判断する傾向にあることが分かった。一方、グループAの周波数派被験者は、周波数特性が高いほど打球音が高いと判断する。したがって、周波数特性のみを基準に打球音高さを評価する場合、実際にゴルファーが感じる打球音との間に差異が生じうることが分かった。音圧特性をも加味して打球音高さを評価することで、より汎用性の高い評価が可能となる。音圧特性が、音の大きさを示す指標であることは知られている。この音圧特性が、打球音高さと相関しうる理由は、上述の通り不明である。ただし、本願にいう打球音はゴルフクラブヘッドでの打球音である。よって、音圧派被験者を生じさせる前提には、ヘッドの打球音における何らかの特徴が関わっていると推測できる。即ち、ヘッド打球音は一般の音とは異なる周波数特性等を有しているので、このヘッド打球音特有の事情が、音圧派被験者を生じさせていると推測される。例えば、一般に、ヘッド打球音では、3kHz以上の複数のピークが観測される。このような周波数特性が、音圧派被験者を生じさせる原因である可能性がある。打球音が比較的大きくなるのは、体積が大きく且つ中空部を有するヘッドである。このヘッドの典型例は、いわゆるドライバーヘッドである。打球音が大きいほど、音の高さが問題となりやすい。この観点から、ヘッド体積は、200cc以上が好ましく、300cc以上がより好ましく、400cc以上がより好ましく、420cc以上がより好ましく、450cc以上が特に好ましい。ルールの観点から、ヘッド体積は、470cc以下が好ましい。ヘッドのタイプとしては、ウッド型、ユーティリティ型及びハイブリッド型が好ましく、ウッド型がより好ましい。本実施形態に係る打球音高さの評価方法は、以下のステップを含む。[St1]:実測又はシミュレーションによってゴルフクラブの打球音を取得する取得ステップSt1[St2]:上記打球音を解析する解析ステップSt2[St3]:上記解析結果に基づいて、打球音高さを評価する評価ステップSt3取得ステップSt1では、実測又はシミュレーションが用いられる。実測では、例えば、所定の位置にマイクロフォン等を配置し、打球音が録音される。シミュレーションでは、例えば、公知の有限要素法が用いられる。この有限要素法では、ヘッドの3次元データ計算モデルが、メッシュ分割される。このメッシュ分割は、例えば、市販のプリプロセッサHyperMesh等によってなされうる。また、ヘッド材質の物性値として、縦弾性率、密度、ポアソン比等が用いられる。シミュレーションによって打球音を解析する場合、ボールが衝突した際のヘッドの表面振動速度又は加速度を境界条件として、音-振動の連成問題を解く。このヘッドの表面振動によって空気が振動し、この空気の振動が音を生じさせる。ヘッド表面の自由振動は、複数の固有モードを重ね合わせることによって表現されうる。それぞれの固有モードは、固有振動数及び振動の形を有している。振動の形は、固有モード形とも称される。ある一つのモードにおいて、表面振動速度を決定する要素は、固有モード形、振幅、固有振動数及びモード減衰比である。表面振動速度は、全てのモードの表面振動速度を重ね合わせることによって表現される。固有モード形及び固有振動数は、固有値解析モード解析によって計算されうる。モード減衰比は、実験等によって求められうる。モーダル解析ソフトウェアとして、例えば、バイブラントテクノロジーVibrantTechnology社製の商品名ME’scopeVESが挙げられる。衝撃応答解析には、CAEソフトウェアが用いられる。CAEソフトウェアとして、リバモアソフトウェアテクノロジー社LivermoreSoftwareTechnologyCorporation製の商品名LS-DYNAが例示される。好ましくは、CAEソフトウェアにより、有限要素法の解析が実行される。衝撃応答解析における拘束条件境界条件は限定されない。打球音はヘッドの自由振動によって発生していることから、拘束条件は自由フリーが好ましい。このましくは、解析ステップSt2は、音圧特性を取得するステップSt21を含む。この音圧特性は、音圧に係る評価指標である。音圧は、空気の圧力変動である。通常、音圧特性が大きいほど、音は大きく感じられる。音圧特性の一例は、音圧レベルである。音圧レベルは、音圧の大きさを、基準値との比の常用対数によって表現した値である。音圧レベルの単位は、デジベルdBである。音圧特性の他の例としてラウドネスが挙げられる。ラウドネスは、心理的な音の大きさを評価する指標として知られている。ISO532Bには、定常ラウドネスレベルの算出方法が定められている。この定常ラウドネスレベルは、ラウドネスの一例である。好ましい音圧特性として、オーバーオール音圧レベルが挙げられる。オーバーオール音圧レベルは、時々刻々と正負の値をとって変化する音圧の実効値に対する音圧レベルであり、単に音圧レベルとも称される。より詳細には、音圧pの実効値peが下記の下記の式1で表され、オーバーオール音圧レベルLpが下記の式2で表される。但し、式1において、Tは音圧の周期であり、pe0は基準音圧である。オーバーオール音圧レベルは、ヘッドとマイクとの位置関係等の条件によって変化しうる。後述される実施形態では、オーバーオール音圧レベルの算出条件は次の通りとされた。即ち、マイク位置がゴルファーに対向する位置に設定された。マイクとボールとの水平方向距離が1mに設定され、マイクとボールとの鉛直方向距離は1.5mに設定された。オーバーオール音圧レベルは、聴覚により感じられる音の大きさとの相関が大きい。このオーバーオール音圧レベルが、音圧派被験者にとっては、打球音高さとの相関が高いことが判明した。好ましくは、評価ステップSt3は、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップSt31を含む。このましくは、この音圧系評価ステップSt31では、上記音圧特性が大きいほど打球音高さが高いと評価される。好ましくは、解析ステップSt2は、音圧特性を取得するステップSt21と、周波数特性を取得するステップSt22とを含む。周波数特性は、周波数に係る評価指標である。この周波数特性として、最大ピーク周波数、1次ピーク周波数、1次固有振動数、エネルギー加重平均周波数及びシャープネスが例示される。1次固有振動数は、前述した固有値解析によって求められうる。1次ピーク周波数と1次固有振動数とは、近似値となることが多い。1次ピーク周波数とは、ピーク周波数のうち最低の周波数である。ただし、特に実測の場合、ノイズ等に起因する微小なピークが発生しやすい。この場合、ノイズ等に起因するピークを除去する観点から、1次ピーク周波数として、所定の閾値Th以上の音圧を有するピークから選択されてもよい。好ましくは、この閾値Thは、最大ピークの音圧を基準として決定される。例えば、この閾値Thは、最大ピークの音圧-20dB以上であり且つ最大ピークの音圧-10dB以下である範囲から選択されうる。また、実測の場合、ボールから発生する音も観測される。ボールから発生する音は、通常、1500Hz近傍である。よって、1500Hz近傍のピーク、特に、1000Hz以上2000Hz以下のピークは、無視するのが好ましい。好ましくは、評価ステップSt3は、上記音圧特性に基づいて打球音高さを評価する音圧系評価ステップSt31と、上記周波数特性に基づいて打球音高さを評価する周波数系評価ステップSt32とを含む。音圧特性と周波数特性との両者が考慮されることで、より汎用性に優れた打球音高さの評価が可能となる。好ましくは、上記評価方法は、以下のステップを更に含む。[St4]周波数特性よりも音圧特性に依存して上記打球音の高さを評価する音圧派被験者と、音圧特性よりも周波数特性に依存して上記打球音の高さを評価する周波数派被験者とを分類する分類ステップSt4[St5]上記音圧派被験者による評価結果を得るステップSt5[St6]上記周波数派被験者による評価結果を得るステップSt6ただし、これらのステップSt4、ステップSt5及びステップSt6は、予め実施しておくことができる。これらの実施により、予め、上記音圧派被験者による評価結果と上記周波数派被験者による評価結果を得ておくことができる。一旦、これらのこれらの評価結果を得ておけば、その後は、これらの評価結果を利用して打球音高さを評価することができる。よって、打球音高さの評価において、上記ステップSt4、ステップSt5及びステップSt6は不要とされうる。好ましくは、音圧系評価ステップSt31では、上記音圧派被験者による評価結果が用いられる。この評価結果の一例は、官能評価の評点である。好ましくは、音圧系評価ステップSt31では、音圧系高さ評価基準が用いられる。この音圧系高さ評価基準として次が例示される。基準b1音圧特性が大きいほど打球音高さが高いと判断される基準b2音圧特性が所定値以上である場合、打球音が高いと判断される。好ましくは、周波数系評価ステップSt32では、上記周波数派被験者による評価結果が用いられる。この評価結果の一例は、官能評価の評点である。好ましくは、周波数系評価ステップSt32では、上記周波数系高さ評価基準が用いられる。この周波数系高さ評価基準として次が例示される。基準a1周波数特性が大きい高いほど打球音高さが高いと判断される基準a2周波数特性が所定値以上である場合、打球音が高いと判断される。好ましくは、上記周波数特性は、1次ピーク周波数又は1次固有振動数である。周波数派被験者において、1次ピーク周波数及び1次固有振動数は、打球音高さとの相関が高い。好ましくは、上記音圧特性が、オーバーオール音圧レベルである。音圧派被験者において、オーバーオール音圧レベルは、打球音高さとの相関が高い。本発明に係るゴルフクラブヘッドの設計方法の一例は、音圧特性が大きいほど打球音高さが高いことを考慮した設計を含む。上述の通り、音圧派被験者の存在が明らかとなったた。よって、この音圧派被験者による評価を設計に取り入れることによって、打球音高さの設計が精度よくなされうる。更に好ましい設計方法の一例は、音圧特性が大きいほど打球音が高いことを考慮した設計と、周波数特性が大きいほど打球音が高いことを考慮した設計とを含む。音圧派被験者と周波数派被験者との両者の評価を取り入れることによって、評価の汎用性がより一層高まる。よって、打球音高さをより精度良く設計することができる。音圧特性を大きくする設計手法として、例えば、ヘッド体積の増加、ヘッド中空部の体積の増加、ヘッド形状の変更、ヘッド厚み例えばクラウン厚み及び/又はソール厚みの調整、材質の変更等が挙げられる。周波数特性を大きくする設計手法として、例えばヘッド厚み例えばクラウン厚み及び/又はソール厚みの調整、ヘッド形状の変更、ヘッド内面へのリブの設置、材質の変更等が挙げられる。図1は、本発明の実施形態の一例を示すフローチャートである。このフローチャートでは、4つの系統のフローが示されている。なお、図1のフローチャートでは、上記4つの系統が1つに纏められているが、上記4つの系統はそれぞれ独立して実施されうることは当然である。[第1系統]この第1系統のフローでは、比較基準とされる打球音が録音されるステップSt10。このステップSt10は、上記取得ステップSt1の一例である。次に、被験者グループによる官能評価がなされるステップSt20。このステップSt20では、被験者は、周波数派被験者と音圧派被験者とに分類されていない。ステップSt20では、ステップSt10で録音された音を被験者グループに聞かせて、官能評価がなされる。この官能評価の一例は、後述される実施例に記載されている。次に、ステップSt20での結果に基づいて、被験者グループが分類されるステップSt70。このステップSt70では、打球音の大きさと打球音の高さとの相関係数が考慮される。この相関係数が低い被験者は、グループA周波数派に分類される。この相関係数が高い被験者は、グループB音圧派に分類される。このステップSt70は、上記分類ステップSt4の一例である。次に、打球音の官能評価値がグループごとに算出されるステップSt80。このステップSt80により、グループA周波数派及びグループB音圧派のそれぞれにおいて、官能評価値が算出される。グループA周波数派の官能評価値を算出することは、上記ステップSt6の一例である。グループB音圧派の官能評価値を算出することは、上記ステップSt5の一例である。このステップSt80により、周波数派被験者による評価値と音圧派被験者による評価値の両方が得られる。これらの評価値を総合して、打球音の高さを評価することができる。即ち、周波数派被験者及び音圧派被験者の両方の要素を考慮して、打球音高さを評価することができる。[第2系統]この第2系統のフローでは、比較基準とされる打球音が録音されるステップSt10。このステップSt10は、上記取得ステップSt1の一例である。次に、打球音の物理指標が算出されるステップSt30。この物理指標は、周波数特性及び音圧特性を含む。このステップSt30は、上記解析ステップSt2の一例である。次に、打球音の官能評価推定値がグループごとに算出されるステップSt90。これらの官能評価推定値は、ステップSt80に基づいてなされうる。即ち、ステップSt80は、物理指標と官能評価との相関が得られうるので、この相関を利用して、ステップSt90では、官能評価値を推定することができる。この官能評価推定値は、グループA周波数派とグループB音圧派とのそれぞれについて得られる。よって、これらの官能評価推定値により、周波数派被験者及び音圧派被験者の両方の要素を考慮して、打球音高さを評価することができる。この第2系統では、分類ステップSt4を行うことなく、官能評価推定値を得ることができる。この第2系統では、被験者による官能評価を行うことなく、官能評価推定値を得ることができる。この第2系統のフローでは、被験者は不要である。ステップSt80において取得されたデータを用いることで、被験者が不要とされている。即ち、周波数派被験者による評価結果と音圧派被験者による評価結果とを予め取得しておけば、それらの評価結果に基づいて、ステップSt90において官能評価値を算出することができる。よって、打球音を取得するだけで打球音高さが評価可能であり、被験者が不要とされうる。[第3系統]この第3系統のフローでは、音の高さを評価したいクラブの打球音が録音されるステップSt40。このステップSt40は、上記取得ステップSt1の一例である。次に、被験者グループによる官能評価がなされるステップSt50。このステップSt50では、被験者は、周波数派被験者と音圧派被験者とに分類されていない。ステップSt50では、ステップSt40で録音された音を被験者グループに聞かせて、官能評価がなされる。この官能評価の一例は、後述される実施例に記載されている。次に、ステップSt50での結果に基づいて、被験者グループが分類されるステップSt70。このステップSt70では、打球音の大きさと打球音の高さとの相関係数が考慮される。この相関係数が低い被験者は、グループA周波数派に分類される。この相関係数が高い被験者は、グループB音圧派に分類される。このステップSt70は、上記分類ステップSt4の一例である。次に、打球音の官能評価値がグループごとに算出されるステップSt80。このステップSt80により、グループA周波数派及びグループB音圧派のそれぞれにおいて、官能評価値が算出される。グループA周波数派の官能評価値を算出することは、上記ステップSt6の一例である。グループB音圧派の官能評価値を算出することは、上記ステップSt5の一例である。このステップSt80により、周波数派被験者による評価値と音圧派被験者による評価値の両方が得られる。これらの評価値を総合して、打球音の高さを評価することができる。即ち、周波数派被験者及び音圧派被験者の両方の要素を考慮して、打球音高さを評価することができる。[第4系統]この第4系統のフローでは、打球音の高さを評価したいクラブの打球音が録音されるステップSt40。このステップSt40は、上記取得ステップSt1の一例である。次に、打球音の物理指標が算出されるステップSt60。この物理指標は、周波数特性及び音圧特性を含む。このステップSt60は、上記解析ステップSt2の一例である。次に、打球音の官能評価推定値がグループごとに算出されるステップSt90。これらの官能評価推定値は、ステップSt80に基づいてなされうる。即ち、ステップSt80では、物理指標と官能評価との相関が得られうるので、この相関を利用して、ステップSt90では、官能評価値を推定することができる。この官能評価推定値は、グループA周波数派とグループB音圧派とのそれぞれについて得られる。よって、これらの官能評価推定値により、周波数派被験者及び音圧派被験者の両方の要素を考慮して、打球音高さを評価することができる。この第4系統では、分類ステップSt4を行うことなく、官能評価推定値を得ることができる。この第4系統では、被験者による官能評価を行うことなく、官能評価推定値を得ることができる。この第4系統のフローでは、被験者は不要である。即ち、周波数派被験者による評価結果と音圧派被験者による評価結果とを予め取得しておけば、それらの評価結果に基づいて、ステップSt90において官能評価値を算出することができる。図2は、他の実施形態を示すフローチャートである。この実施形態は、ヘッドの設計に適用されうる。この実施形態では、音の高さを評価したいクラブの打球音が、実打の録音又はシミュレーションによって得られるステップSt100。このステップSt100は、上記取得ステップSt1の一例である。次に、ステップSt100で得られたデータに基づいて、周波数特性が算出されるステップSt200。このステップSt200は、解析ステップSt2の一例である。このステップSt200は、上記ステップSt22の一例である。また、ステップSt100で得られたデータに基づいて、音圧特性が算出されるステップSt300。このステップSt300は、解析ステップSt2の一例である。また、ステップSt300は、上記ステップSt21の一例である。次に、打球音の高さの評価値がグループごとに算出されるステップSt400。前述したステップSt70で分類されたグループごとに、評価値が算出される。このステップSt400は、音圧派被験者による評価結果を得るステップSt5と、周波数派被験者による評価結果を得るステップSt6とで得られた評価結果が用いられる。即ち、周波数派被験者による評価結果と音圧派被験者による評価結果とを予め取得しておけば、それらの評価結果に基づいて、官能評価値を算出することができる。したがって、ステップSt400では、被験者は不要とされている。次に、打球音高さが判断されるステップSt500。好ましくは、このステップSt500では、グループA周波数派による評価値と、グループB音圧派による評価値とを総合して、打球音高さが判断される。よって、打球音高さの判断精度が向上しうる。打球音高さが適切でないと判断された場合、打球音が低い構造要因が例示されるステップSt600。ソフトウェアには、音圧特性及び周波数特性を変更しうる構造要因が予め記憶されている。このソフトウェアにより、打球音高さを改善するための構造が例示される。ステップSt600において打球音高さが適切であると判断された場合、設計は終了する。この実施形態では、分類ステップSt4が不要とされうる。上記ステップSt400においては、被験者による官能評価に変えて、官能評価推定値が算出されてもよい。即ち、このステップSt400での算出値は、前述のステップSt90と同様に、推定値とされてもよい。この場合、被験者による官能評価が省略されうる。以下、実施例によって本発明の効果が明らかにされるが、この実施例の記載に基づいて本発明が限定的に解釈されるべきではない。肉厚分布及び体積が相違する7種類のヘッドH1からH7が用意された。全てウッド型の中空ヘッドとされた。全てがドライバーヘッドであった。ヘッド体積は、ヘッドH1が305ccであり、ヘッドH2が350ccであり、ヘッドH3が407ccであり、ヘッドH4が432ccであり、ヘッドH5が455ccであり、ヘッドH6が460ccであり、ヘッドH7が460ccであった。全てのヘッドにおいて、材質はチタン合金であった。ヘッドH1にシャフト及びグリップを装着して、クラブ1を得た。ヘッドH2にシャフト及びグリップを装着して、クラブ2を得た。ヘッドH3にシャフト及びグリップを装着して、クラブ3を得た。ヘッドH4にシャフト及びグリップを装着して、クラブ4を得た。ヘッドH5にシャフト及びグリップを装着して、クラブ5を得た。ヘッドH6にシャフト及びグリップを装着して、クラブ6を得た。ヘッドH7にシャフト及びグリップを装着して、クラブ7を得た。シミュレーションにより、各ヘッドの打球音を得た。前述した方法によりヘッド表面の自由振動を求め、更に、空気中の伝搬を計算して、ヘッドから所定距離隔てた位置に到達する音のスペクトル波形を得た。なお、このスペクトル波形の算出は以下の通りに行った。衝突解析により、インパクト中にフェース面が受ける力Ffを得た。ヘッド速度を40m/sとし、静止したボールに、所定のロフト角となる姿勢でヘッドを衝突させた。衝突解析ソフトウェアとしてLS-DYNAを用いた。ボールの弾性率は50MPa~250MPaとし、粘性は無視した。ボールのポアソン比は0.463とした。ボールの密度は1.24×10-9t/mm3とした。ヘッドについては、弾性率は120MPaとされ、ポアソン比は0.3とされ、密度は4.2×10-9t/mm3とした。次に、固有値解析を行い、ヘッドのモード形及び固有振動数を得た。固有値解析ソフトウェアとして、、MSCソフトウェア社のNastranを用いた。ヘッドの物性値は上述の衝突解析と同じとし、自由支持条件における固有モードを算出した。次に、音響解析を行った。音響解析ソフトとして日本イーエスアイ社のVA-oneを用い、打球音のスペクトル波形を得た。前述の力Ffと固有モードとを入力し、ヘッドと観測点との位置関係を、前述の通り、水平距離1m且つ鉛直距離1.5mに設定して、放射音を計算した。入力された固有モードとは、具体的には、固有モード形、固有振動数及びモード質量である。放射音の計算では、空気の密度を1.18kg/m3とし、音速を343m/sとした。ヘッドH1クラブ1のスペクトル波形が、図3に示される。ヘッドH2クラブ2のスペクトル波形が、図4に示される。ヘッドH3クラブ3のスペクトル波形が、図5に示される。ヘッドH4クラブ4のスペクトル波形が、図6に示される。ヘッドH5クラブ5のスペクトル波形が、図7に示される。ヘッドH6クラブ6のスペクトル波形が、図8に示される。ヘッドH7クラブ7のスペクトル波形が、図9に示される。なお、図3から図9には、ピーク周波数の位置が縦線で示されている。例えば図3のヘッドH1では、1次ピーク周波数は、5225Hzである。全てのヘッドH1からH7において、1次ピーク周波数は、3500Hz以上であった。全てのヘッドH1からH7において、最大ピーク周波数は、3500Hz以上であった。図3から図9には、スペクトル波形のグラフとともに、各ピーク周波数と音圧との対応表が示されている。これらのスペクトル波形では、ピークの位置が示されている。ピークの判別方法は次の通りである。スペクトルのピーク周波数はヘッドの固有振動数に相当することから、本実施例では、計算で得られた固有振動数をピーク周波数とした。実測された音の場合は、公知の実験モーダル解析手法を用いて、ピーク周波数を判別できる。スペクトル波形から判別しにくい場合、振動実験を行うことでピーク周波数を判別することができる。これらのスペクトル波形を逆フーリエ変換して、時間を横軸とする波形を得た。アンプ及びスピーカーを有する再生機を用いて、この波形の音を再生した。11名のゴルファーが、再生された音を聞いて、官能評価を行った。この評価に用いられた評価用紙を図10に示す。この官能評価では、一対比較法が用いられた。本実施例では、この一対比較法として、シェッフェの方法の浦の変法が採用された。この一対比較法により、7種類の音に対して、順序効果を考慮し、42通りの比較を行った。比較される2つの音を順次再生し、先の音に対して後の音がどのように感じられたかが評価された。評価は、図10に示される5段階評価によってなされた。音の大きさ及び音の高さのそれぞれについて、評価がなされた。得られた評価結果をクラスタ分析した。この結果、11名のゴルファーのうちの7名が周波数派被験者グループAに分類された。残り4名は、音圧派被験者グループBに分類された。このクラスタ分析の具体的方法は次の通りである。距離類似度としてユークリッド距離を用い、クラスター合併のアルゴリズムとしてウォード法を用いた。ソフトウェアとして統合解析ソフトRを用いた。図11は、音圧派被験者グループBにおける、大きさと高さとの相関を示す散布図である。横軸の大きさは、一対対比法で得られた打球音の大きさの評点42データの平均値であり、縦軸の高さは、一対対比法で得られた打球音の高さの評点42データの平均値である。この図11のグラフが示すように、音圧派被験者では、打球音の大きさと打球音の高さとの相関が高い。この図11のグラフにおいて、相関係数Rは0.97であった。本願にいう相関係数Rは、ピアソンの積率相関係数である。この図11により、音圧派被験者の存在が示されている。図12は、周波数派被験者グループAにおける、大きさと高さとの相関を示す散布図である。横軸及び縦軸の意味は図11のグラフと同じである。この図12のグラフが示すように、音圧派被験者では、打球音の大きさと打球音の高さとの相関が低い。この図12のグラフにおいて、相関係数Rは-0.65であった。このように、周波数派被験者と音圧派被験者とで、相関係数Rに大きな差が見られた。図13は、オーバーオール音圧レベルと大きさの評価点との相関である。この図13は、11名全員のデータに基づく。この図13が示すように、オーバーオール音圧レベルと、聴覚上の打球音の大きさとは、相関が高い。この図13の結果と、図11との結果とを考慮すると、音圧派被験者グループBでは、オーバーオール音圧レベルと打球音の高さとの相関が高いことが理解できる。これらの結果は、音圧特性に基づいて打球音の高さを評価することの有効性を示している。次に、周波数派被験者グループAにおいて、好ましい周波数特性の指標を探索した。上記7種のヘッドH1からH7のシミュレーション結果及びそれぞれの官能評価を纏めたのが、下記の表1である。この表1の官能評価の欄には、上記7名の周波数派被験者による、打球音高さの官能評価点が示されている。これらの官能評価点をP2とし、上記一対比較法で得られた上記評価点をP1とするとき、P2は次のように計算された。P2=P1×3/24表1のうち、1次固有振動数以外は、図3から図9のシミュレーション波形により求めた。1次固有振動数は、モーダル解析の結果から求めた。図14から図19は、打球音高さ評価と各周波数特性との相関を示す散布図である。また、上記表1には、これら散布図で示されていないものも含めて、音の高さ評価と各周波数特性との相関係数Rが示されている。表1が示すように、相関係数Rが特に高い周波数特性は、1次ピーク周波数及び1次固有振動数である。表1には、エネルギー加重平均周波数が示されている。このエネルギー加重平均周波数の計算方法は、前述した特開2008-132262号公報に記載された方法と同じである。周波数範囲を限定して、複数のエネルギー加重平均周波数を算出した。これらのうち、特開2008-132262号公報に記載されているエネルギー加重平均周波数に相当するのは、10kHz以下の全ピークである。この場合の相関係数Rは0.77であり、1次ピーク周波数の相関係数R0.90よりも低く、1次固有振動数の相関係数R0.86よりも低い。上記クラブ1から7を用いて,実打での打球音を用いた官能評価テストも行った。この官能評価においても、図10の評価用紙が用いられた。テスターは、上記11名であった。この実打音でのテストの結果、図11、図12及び13と同様の傾向が見られた。"