(57)【要約】  【課題】上記問題を解決するために、評価対象とするゴルフクラブヘッドの打音が、任意のゴルファからどの程度好まれるか定量的に評価する、ゴルフクラブヘッドの打音評価方法を提供する。  【解決手段】評価対象のゴルフクラブヘッドで所定のゴルフボールを打撃した際の振動波形をサンプリングし、サンプリングした振動波形から、打音の大きさを表す音圧パラメータと、打音の減衰の度合を表す残響パラメータと、打音の周波数特性を表す周波数パラメータと、打音の澄み感を表す純音度パラメータとを、評価パラメータとして算出し、各評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値を算出し、この評価値に応じて、評価対象のゴルフクラブの打音を評価する。【選択図】図２
2006239132.tif 000002
【特許請求の範囲】
【請求項１】
ゴルフクラブヘッドでゴルフボールを打撃した際の打音を評価するゴルフクラブヘッドの打音評価方法であって、
評価対象のゴルフクラブヘッドでゴルフボールを打撃した際の振動波形のサンプリングを行なうサンプリングステップと、
サンプリングされた振動波形から、前記打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを、評価パラメータとして算出するパラメータ算出ステップと、
前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を表す好感評価値を算出する評価値算出ステップと、
前記好感評価値に応じて、前記評価対象のゴルフクラブの打音を評価する評価ステップとを有することを特徴とするゴルフクラブヘッドの打音評価方法。
【請求項２】
前記評価値算出ステップは、前記好感評価値に加え、前記周波数パラメータおよび前記純音度パラメータそれぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際、この打音に対して感じる音程の高さの程度を表す高低評価値を算出することを特徴とする請求項１に記載のゴルフクラブヘッドの打音評価方法。
【請求項３】
前記評価値算出ステップは、前記好感評価値に加え、前記周波数パラメータおよび前記純音度パラメータそれぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際、この打音に対して感じる音の澄み具合の程度を表す澄度評価値を算出することを特徴とする請求項１または２に記載のゴルフクラブヘッドの打音評価方法。
【請求項４】
前記パラメータ算出ステップは、サンプリングした振動波形全体から、振動波形の絶対値が予め定められたトリガー値以上となる、時間軸上で最初のサンプル番号を抽出し、このサンプル番号から所定のサンプル数だけ時間軸上で遡ったサンプル番号を始点番号とし、この始点番号から所定の数のサンプルで構成された第１サンプル群を用いて、前記評価パラメータを算出することを特徴とする請求項１〜３のいずれかに記載のゴルフクラブヘッドの打音評価方法。
【請求項５】
前記パラメータ算出ステップは、サンプリングした振動波形を周波数解析し、得られたスペクトル波形を用い、前記音圧パラメータとして前記打音のラウドネスレベルの値を算出することを特徴とする請求項１〜４のいずれかに記載のゴルフクラブヘッドの打音評価方法。
【請求項６】
前記パラメータ算出ステップは、前記第１サンプル群において、前記始点番号から所定のサンプル数で振動波形が順次区切られた複数のブロックからなるブロック群を作成し、このブロック群の各ブロック毎に、振動波形のレベルの二乗和の対数値を各ブロックのブロック値として求め、このブロック値を用いて前記残響パラメータを算出する請求項４または５に記載のゴルフクラブヘッドの打音評価方法。
【請求項７】
前記パラメータ算出ステップは、前記ブロック群におけるブロック値の最大値を用い、この最大値から所定の値だけ低い終点ブロック値を設定し、ブロック値が前記最大値となる起点ブロック、およびこの起点ブロック以降の複数のブロックにおいて、ブロック値が最初に前記起点ブロック値以下となる終点ブロックを抽出し、前記起点ブロックを構成するサンプルのうち最初のサンプルと、前記終点ブロックを構成するサンプルのうち最初のサンプルとの、取得時間の間隔を前記残響パラメータとして算出することを特徴とする請求項６記載のゴルフクラブヘッドの打音評価方法。
【請求項８】
前記パラメータ算出ステップは、前記第１サンプル群から前記始点番号以降の所定の数のサンプルを除いて得られる第２サンプル群を用いて、前記周波数パラメータおよび前記純音度パラメータをそれぞれ求めることを特徴とする請求項４〜７のいずれか１項に記載のゴルフクラブヘッドの打音評価方法。
【請求項９】
前記パラメータ算出ステップは、前記第２サンプル群について周波数解析を行なってスペクトル波形を求め、このスペクトル波形の周波数に対するスペクトル値を重み付け係数として周波数の加重平均値を求め、この加重平均値を用いて周波数パラメータを求める請求項８記載のゴルフクラブヘッドの打音評価方法。
【請求項１０】
前記パラメータ算出ステップは、前記第２サンプル群について周波数解析を行なってスペクトル波形を求め、このスペクトル波形の周波数解析の対象とする周波数範囲における平均スペクトル値より高いスペクトル波形の第１の合計値と、この周波数範囲全体におけるスペクトル波形の第２の合計値を求め、この第２の合計値に対する第１の合計値の比率を求め、この比率を純音度パラメータとする請求項８または９記載のゴルフクラブヘッドの打音評価方法。
【請求項１１】
前記サンプリングステップにおいてサンプリングする前記打音の振動波形は、人間の頭部を模したダミーヘッドの耳に対応する部分に設置された、音波情報取得手段によって取得された波形であることを特徴とする請求項１〜１０のいずれかに記載ゴルフクラブヘッドの打音評価方法。
【請求項１２】
ゴルフクラブヘッドの打撃面を打撃した際のゴルフクラブヘッドの打音をコンピュータに評価させるコンピュータが実行可能なプログラムであって、
評価対象のゴルフクラブヘッドでゴルフボールを打撃した際の振動波形から、前記打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを評価パラメータとして、コンピュータの演算手段に算出させ、コンピュータの記憶手段に記憶させるパラメータ算出手段と、
前記記憶手段に予め記憶されている重み付け係数を、前記評価パラメータとともに呼びださせ、前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値を算出させ、この評価値に応じた前記打音の評価を前記演算手段に行なわせる評価手順とを有するプログラム。
【請求項１３】
ゴルフクラブヘッドの打撃面を打撃した際のゴルフクラブヘッドの打音を評価するゴルフクラブヘッドの打音評価システムであって、
評価対象のゴルフクラブヘッドでゴルフボールを打撃した際の振動波形を計測する計測装置と、
振動波形から、打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを、評価パラメータとして算出するパラメータ算出部と、
前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値を算出する評価値算出部と、
前記評価値に応じて前記打音を評価する評価部を備える処理装置とを有することを特徴とするゴルフクラブヘッドの打音評価システム。
【請求項１４】
前記計測装置は、人間の頭部を模したダミーヘッドと、このダミーヘッドの耳に対応する部分に設置された音波情報取得手段とを備えることを特徴とする請求項１３に記載のゴルフクラブヘッドの打音評価システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴルフクラブによってゴルフボールを打撃した際の打音の評価を、定量的に精度良く行なうゴルフクラブヘッドの打音評価方法、この方法をコンピュータに実行させるプログラム、およびゴルフクラブの打音評価システムに関する。
【背景技術】
【０００２】
近年、ゴルフボールをより遠くに飛ばすことのできる反発係数の高い材料を打撃面に用いたゴルフクラブヘッドや、スイートスポットが大きく安定したゴルフボールの飛距離を確保することのできる容積の大きなゴルフクラブヘッドが種々提供されている。このようなゴルフクラブヘッドは、例えば、打撃面に特殊な材料を用いた中空ヘッドであるため、打撃面への打撃の際に発生する打音がそれぞれ大きく異なる。このため、実際にはゴルフボールを遠くに飛ばすことができるにもかかわらず、打音自体が鈍い音を発して、飛距離の出ないゴルフクラブとしての印象をゴルファに与えたり、不快なイメージを与える場合も多い。このような背景から、ゴルフクラブによるゴルフボールの打撃時の打音を調整することが望まれている。
【０００３】
下記特許文献１では、ドライバー等のウッド系のゴルフクラブヘッド又はアイアンのゴルフクラブヘッドの打撃面を打撃した際のゴルフクラブヘッドの打音を評価する際、設計変更による変更部分が僅かなゴルフクラブヘッドであっても精度高く打音を評価できるとともに、ゴルファの官能評価に対応した評価パラメータを算出して評価を行なうことができ、しかも、打音の評価結果の主要因となるスペクトル波形のピーク値またはピーク周波数を容易に特定することができるゴルフクラブヘッドの打音評価方法が提案されている。この特許文献１では、評価対象のゴルフクラブヘッドの打撃面を打撃した際の振動波形のサンプリングを行ない、サンプリングされた振動波形から、打音の大きさを表す音圧パラメータと、打音の減衰の度合を表す残響パラメータと、打音の高低を表す周波数パラメータと、打音の澄み具合を表す音色パラメータを評価パラメータとして算出する。そして、算出された評価パラメータを用い、この評価パラメータが既知であるゴルフクラブヘッドに対する評価対象のゴルフクラブヘッドの相対評価を行なっている。
【特許文献１】特開２００３−３２５７０７号公報
【０００４】
上記特許文献１では、既知の複数のゴルフクラブヘッドについて、上述の音圧パラメータ、残響パラメータ、周波数パラメータ、および音色パラメータのそれぞれを算出しておき、評価対象とするゴルフクラブヘッドの各パラメータを対応するパラメータに加えて母集団として平均値と標準偏差とを求め、この母集団の平均値を５０、標準偏差を１０とする変換を行なって、評価対象とするゴルフクラブヘッドの各パラメータの偏差値をそれぞれ求めている。そして、各ゴルフクラブ毎の各パラメータの偏差値の値をレーダチャートに表すことで、偏差値を基準にして各ゴルフクラブの打音の特性を評価している。特許文献１では、このように、複数のゴルフクラブヘッドについて評価パラメータを算出し、評価対象とするゴルフクラブの各評価パラメータの偏差値を求めている。特許文献１では、このようにして、ゴルフクラブヘッドの設計開発段階において効果的に活用することがゴルフクラブヘッドの打音評価結果を求めている。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ゴルフクラブの設計開発段階において、多くのゴルファから好まれる打音のゴルフクラブヘッドを設計開発するには、どのような打音が任意のゴルファから好まれやすいかの指標を得ることが重要である。このような指標を得るためには、任意のゴルファがゴルフクラブヘッドの打音を聞いた際に、この任意のゴルファがもつ打音に対する好感の程度を、ゴルフクラブヘッド毎に定量的に評価することが重要である。このように、ゴルフクラブヘッドの打音を任意のゴルファが聞いた際に、任意のゴルファがもつ好感の程度を、ゴルフクラブヘッド毎に定量的に評価できれば、複数のゴルフクラブヘッドの打音について、どのゴルフクラブヘッドの打音が、任意のゴルファにどの程度好まれやすいかが定量的に判断できる。このような情報を、ゴルフクラブヘッドの仕様の選択や、設計の指針に反映させることで、任意のゴルファからより好まれやすいゴルフクラブヘッドを設計することが可能となる。特許文献１記載のゴルフクラブヘッドの打音評価方法では、評価対象とするゴルフクラブヘッドの打音が、任意のゴルファからどの程度好まれるかの度合（打音の好感の程度）に関しては、評価することができない。また、このような特許文献１の評価方法では、母集団自体が異なると（すなわち、母集団を構成する複数のゴルフクラブヘッドが異なると）、評価対象とするゴルフクラブヘッドの打音の評価パラメータの偏差値も異なってしまう。このような評価方法では、評価対象とするゴルフクラブヘッドの打音が、母集団において比較的小さい、または比較的大きいなど、打音の特徴を定性的に評価することは可能であるが、各評価パラメータを定量的に評価することはできないといった問題もあった。
【０００６】
そこで、本発明は、上記問題を解決するために、評価対象とするゴルフクラブヘッドの打音が、任意のゴルファからどの程度好まれるかの度合（打音の好感の程度）を定量的に評価することを可能とする、ゴルフクラブヘッドの打音評価方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明は、ゴルフクラブヘッドでゴルフボールを打撃した際の打音を評価するゴルフクラブヘッドの打音評価方法であって、評価対象のゴルフクラブヘッドで所定のゴルフボールを打撃した際の振動波形のサンプリングを行なうサンプリングステップと、サンプリングされた振動波形から、前記打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを、評価パラメータとして算出するパラメータ算出ステップと、前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値を算出する評価値算出ステップとを有することを特徴とするゴルフクラブヘッドの打音評価方法を提供する。
【０００８】
また、前記評価値算出ステップは、前記好感評価値に加え、前記周波数パラメータおよび前記純音度パラメータそれぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際、この打音に対して感じる音程の高さの程度を表す高低評価値を算出することが好ましい。
【０００９】
また、前記評価値算出ステップは、前記好感評価値に加え、前記周波数パラメータおよび前記純音度パラメータそれぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際、この打音に対して感じる音の澄み具合の程度を表す澄度評価値を算出することが好ましい。
【００１０】
なお、前記パラメータ算出ステップは、サンプリングした振動波形を周波数解析し、得られたスペクトル波形を用い、前記音圧パラメータとして前記打音のラウドネスレベルの値を算出することが好ましい。
【００１１】
また、前記パラメータ算出ステップは、サンプリングした振動波形全体から、振動波形の絶対値が予め定められたトリガー値以上となる、時間軸上で最初のサンプル番号を抽出し、このサンプル番号から所定のサンプル数だけ時間軸上で遡ったサンプル番号を始点番号とし、この始点番号から所定の数のサンプルで構成された第１サンプル群を用いて、前記残響パラメータ、前記周波数パラメータ、および前記純音度パラメータを算出するとよい。
【００１２】
また、前記パラメータ算出ステップは、前記第１サンプル群において、前記始点番号から所定のサンプル数で振動波形が順次区切られた複数のブロックからなるブロック群を作成し、このブロック群の各ブロック毎に、振動波形のレベルの二乗和の対数値を各ブロックのブロック値として求め、このブロック値を用いて前記残響パラメータを算出することが好ましい。
【００１３】
また、前記パラメータ算出ステップは、前記ブロック群におけるブロック値の最大値を用い、この最大値から所定の値だけ低い終点ブロック値を設定し、ブロック値が前記最大値となる起点ブロック、およびこの起点ブロック以降の複数のブロックにおいて、ブロック値が最初に前記起点ブロック値以下となる終点ブロックを抽出し、前記起点ブロックを構成するサンプルのうち最初のサンプルと、前記終点ブロックを構成するサンプルのうち最初のサンプルとの、取得時間の間隔を前記残響パラメータとして算出することが好ましい。
【００１４】
また、前記パラメータ算出ステップは、前記第１サンプル群から前記始点番号以降の所定の数のサンプルを除いて得られる第２サンプル群を用いて、前記周波数パラメータおよび前記純音度パラメータをそれぞれ求めることが好ましい。
【００１５】
また、前記パラメータ算出ステップは、前記第２サンプル群について周波数解析を行なってスペクトル波形を求め、このスペクトル波形の周波数に対するスペクトル値を重み付け係数として周波数の加重平均値を求め、この加重平均値を用いて周波数パラメータを求めることが好ましい。
【００１６】
また、前記パラメータ算出ステップは、前記第２サンプル群について周波数解析を行なってスペクトル波形を求め、このスペクトル波形の周波数解析の対象とする周波数範囲における平均スペクトル値より高いスペクトル波形の第１の合計値と、この周波数範囲全体におけるスペクトル波形の第２の合計値を求め、この第２の合計値に対する第１の合計値の比率を求め、この比率を純音度パラメータとすることが好ましい。
【００１７】
なお、前記サンプリングステップにおいてサンプリングする前記打音の振動波形は、人間の頭部を模したダミーヘッドの耳に対応する部分に設置された、音情報取得手段によって取得された波形であることが好ましい。
【００１８】
また、本発明は、ゴルフクラブヘッドの打撃面を打撃した際のゴルフクラブヘッドの打音をコンピュータに評価させるコンピュータが実行可能なプログラムであって、評価対象のゴルフクラブヘッドでゴルフボールを打撃した際の振動波形から、前記打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを評価パラメータとして、コンピュータの演算手段に算出させ、コンピュータの記憶手段に記憶させるパラメータ算出手段と、前記記憶手段に予め記憶されている重み付け係数を、前記評価パラメータとともに呼びださせ、前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値の算出を前記演算手段に行なわせる評価手順とを有するプログラムも併せて提供する。
【００１９】
本発明は、また、ゴルフクラブヘッドの打撃面を打撃した際のゴルフクラブヘッドの打音を評価するゴルフクラブヘッドの打音評価システムであって、評価対象のゴルフクラブヘッドでゴルフボールを打撃した際の振動波形を計測する計測装置と、振動波形から、打音の大きさを表す音圧パラメータと、前記打音の減衰の度合を表す残響パラメータと、前記打音の周波数特性を表す周波数パラメータと、前記打音の澄み感を表す純音度パラメータとを、評価パラメータとして算出するパラメータ算出部、および、前記評価パラメータの値それぞれに対して重み付けを行なって、任意の人物が前記打音を聞いた際に感じる、前記打音に対する好感の程度を定量的に表す好感評価値を算出する評価部を備える処理装置とを有することを特徴とするゴルフクラブヘッドの打音評価システムも併せて提供する。
【００２０】
また、前記計測装置は、人間の頭部を模したダミーヘッドと、このダミーヘッドの耳に対応する部分に設置された音情報取得手段とを備えることが好ましい。
【発明の効果】
【００２１】
本発明のゴルフクラブヘッドの打音評価方法によれば、評価対象とするゴルフクラブヘッドの打音を任意のゴルファが聴いた際に、この任意のゴルファがもつ打音に対する好感の程度を、高い精度で定量的に表すことができる。このような好感の程度の情報は、ゴルフクラブの設計開発段階において効果的に活用することができる。
【発明を実施するための最良の形態】
【００２２】
本発明のゴルフクラブヘッドの打音評価方法、この方法を実行させるプログラムおよびゴルフクラブヘッドの打音評価システムについて、以下詳細に説明する。
図１（ａ）は、本発明のゴルフクラブヘッドの打音評価システムの一例である、ゴルフクラブヘッドの打音評価システム１０（以降、評価システム１０とする）の概略の構成を示す概略斜視図である。また、図１（ｂ）は、ゴルフクラブヘッドの打音評価システム１０の概略の構成を示す概略平面図である。
【００２３】
評価システム１０は、スウィングロボット１２と、ゴルフボール１３と、バイノーラル録音システム２０と、マイクロフォンアンプ１９と、処理装置１８とを有して構成される。
スウィングロボット１２は、評価対象とするゴルフクラブヘッド３２を有するゴルフクラブ３０を把持した状態で、所定のスウィングパターンでゴルフスウィングを行い、ゴルフクラブヘッド３２によって、ゴルフボール１３を打撃する公知のゴルフスウィングロボットである。
【００２４】
バイノーラル録音システム２０は、人間の頭部を模したダミーヘッド１６の、外耳道入口部分に対応する箇所に騒音マイク１４が備えられて構成されている。騒音マイク１４は、スウィングロボット１２によって、ゴルフクラブヘッド３２がゴルフボール１３を打撃した際に発する打音の音波情報（音圧の変化）を受けて電気信号に変える公知の騒音マイクである。バイノーラル録音システム２０では、このように、人間の頭部を模したダミーヘッドの外耳道入口部分にマイクをセットすることにより、任意の人物の頭部がこのダミーヘッド１６と同じ位置に位置した際に、この人物の外耳道入口部分に到達する打音の音波情報と同様の音波情報を、騒音マイク１４が受けることを可能としている。これにより、騒音マイク１４によって、任意の人物の頭部がこのダミーヘッド１６と同じ位置に位置した際に、この人物が聞く打音の音波情報（音圧の変化）とほぼ同じ音波情報を表す電気信号を取得することができる。
【００２５】
本実施形態において、ダミーヘッド１６は、ゴルフボールの設置位置を挟んでスウィングロボット１２と対向する側に、ゴルフクラブヘッド３２によるゴルフボール１３の打撃地点（すなわち、ゴルフボールの設置位置）を基準位置とし、地面と水平方向に１．００（ｍ）、地面からの高さ約１．７（ｍ）の位置に位置するよう設定した。スウィングロボット１２は、所定のゴルファがゴルフスウィングを行なった際のゴルフスウィングを再現するものであるが、このダミーヘッド１６の位置は、ゴルフボールの設置位置を通る地面に垂直な仮想線を軸とし、このスウィングロボットが再現するゴルフスウィングにおけるゴルファの頭部位置に対応する位置（図中Ｈで示す）と、対称な位置である。このような位置に設置されたダミーヘッドに届く打音は、スウィングロボット１２が再現するゴルフスウィングを実際にゴルファが行なった際、このゴルファの頭部（図１中のＨの位置）に届く打音（ゴルファが聞く打音）とほぼ同じになる。本実施形態では、このような位置にダミーヘッド１６を配置し、ゴルフスウィングを行なうゴルファが実際に聞く打音の音波情報と同様の音波情報を取得できる。バイノーラル録音システム１０では、ダミーヘッド１６の左右の耳に対応する部分それぞれに、騒音マイク１４が設けられており、各マイクそれぞれにおいて音波情報を取得する。
【００２６】
騒音マイク１４で取得された音波情報の電気信号は、マイクロフォンアンプ１９に送られる。マイクロフォンアンプ１９は、騒音マイク１４で取得された音波情報の電気信号を増幅して処理装置１８に送る。処理装置１８は、バイノーラル録音システム２０の騒音マイク１４で得られた電気信号をサンプリング（サンプリング周波数３０ｋＨｚ〜１００ｋＨｚ）して取り込み、得られた振動波形のデータを打音ファイルとして図示しない記憶手段に記憶する。そして、図示しない記憶手段２４に記憶されたこの打音ファイル（振動波形のデータ）から、音圧パラメータ、残響パラメータ、周波数パラメータおよび純音度パラメータをそれぞれ算出する。そして、これらのパラメータを重み付けして、任意の人物が上記打音を聞いた際に感じる、打音に対する好感の程度を表す好感評価値を算出して出力する。音圧パラメータ、残響パラメータ、周波数パラメータおよび純音度パラメータの算出、およびこれらのパラメータを用いて算出する好感評価値については後述する。
【００２７】
処理装置１８は、図示しないＣＰＵおよびメモリを備え、プログラムを実行することによって各部位が機能するコンピュータによって構成されている。処理装置１８では、ゴルフクラブヘッドによってゴルフボールを打撃した際の打音の振動波形を用いて、ゴルフクラブヘッドの打音をコンピュータに評価させる、コンピュータが実行可能なプログラムが用いられている。この場合、プログラムは、振動波形から、打音の大きさを表す音圧パラメータと、打音の減衰の度合を表す残響パラメータと、打音の周波数特性を表す周波数パラメータと、打音の硬さ感を表す純音度パラメータを評価パラメータとして、コンピュータのＣＰＵに算出させ、コンピュータのメモリに記憶させる。そして、予めメモリに記憶されている各パラメータ毎の重み付け係数を、記憶された各評価パラメータそれぞれとともに呼びださせ、評価対象のゴルフクラブヘッドの打音に対する好感の程度を表す好感評価値を算出する。プログラムがコンピュータに実行させる各手順については、後に詳述する。なお処理装置１８には、ディスプレイなどの表示手段２６、およびキーボードやマウスなどの入力手段２８が備えられており、プログラムによる各手順の実行は、オペレータが表示手段２６に表示される画面を見ながら入力手段２６を操作して行う、いわゆるＧＵＩ（ｇｒａｐｈｉｃａｌ  ｕｓｅｒ  ｉｎｔｅｒｆａｃｅ）による指示入力に応じて行われる。処理装置１８は、プログラムを実行することによって機能を発揮するコンピュータによって構成されることに限定されず、各手順をそれぞれ実行する複数の専用回路で構成されたものであってもよい。
【００２８】
図２は、評価システム１０で行われるゴルフクラブヘッドの打音評価方法の一実施形態の流れを示すフローチャートである。以降、評価システム１０で実施する、本発明のゴルフクラブヘッドの打音評価方法の一実施形態について説明する。
【００２９】
まず、スウィングロボット１２において、ゴルフクラブヘッド３２を把持した状態でゴルフスウィングが実施され、ゴルフクラブヘッド３２でゴルフボール１３を打撃した打音が、バイノーラル録音システム２０の騒音マイク１４で取得されて、打音の音波情報に応じた電気信号が処理装置１８に送られる。処理装置１８では、この電気信号をサンプリングして取り込み、得られた振動波形のデータを図示しない記憶手段に記憶する（ステップＳ１０２）。本実施形態では、電気信号を、サンプリング周期４４．１（ｋＨｚ）でサンプリングして取り込む。
【００３０】
処理装置１８では、この記憶したゴルフクラブヘッド３２の振動波形のデータを用い、打音評価を行なう。振動波形のデータが記憶されると、処理装置１８において音質評価プログラムが実行され、各種パラメータおよび評価値の算出が開始される。図３は、本発明の打音評価プログラムの一例である音質評価プログラムを実行した際、評価装置１８のディスプレイ２６に表示される画面の一例を示す図である。音質評価プログラムが実行されると、ディスプレイ２２には、この図３のようなプログラム実行制御画面３４が表示される。プログラム実行制御画面３４には、入力ボックス３６ａ〜３６ｃが設けられており、オペレータが入力手段２８（キーボードおよびマウスなど）を操作することで、所望の数値やファイル名などが入力可能となっている。また、プログラム実行制御画面には、計算実行指示ボタン３８ａ〜３８ｅが設けられており、オペレータが入力手段２８を操作することで、処理装置１８による各評価パラメータの算出や、好感評価値の算出の実行を指示する。
【００３１】
処理装置１８では、まず、サンプリングして得られた振動波形のデータ全体から、打音の評価に用いる振動波形データを既定する条件を設定する（ステップＳ１０４）。具体的には、評価対象のゴルフクラブヘッド３２でゴルフボール１３を打撃した瞬間を基準に、この打撃の前後一定時間の振動波形のデータを、打音の評価に用いる振動波形データとして既定するための条件を設定する。この設定は、オペレータが、ディスプレイ２６を見ながら入力装置２８を操作して、プログラム実行制御画面の、評価データ個数入力ボックス３６ａ、トリガーレベル入力ボックス３６ｂ、プリトリガー値入力ボックス３６ｃ、の各入力ボックスに所望の設定値を入力することで行なわれる。本実施形態では、打音の評価に用いる振動波形データデータ個数を１６３８４（個）、トリガーレベルを０．１（Ｖ）、プリトリガー値を２００（個）と設定する。以降のステップにおいては、振動波形が取得した振動波形データのうち、振動波形の電圧値がこのトリガーレベル（０．１（Ｖ））を超えたタイミングのサンプリングデータを基準に、プリトリガー数値として設定された所定のデータ数（２００個）だけ、時間軸上で遡ったサンプリングデータを始点データ（始点番号のサンプル）とし、この始点データから所定数のデータ（１６３８４個）を抽出する。そして、抽出したこの所定数のデータについてハイパスフィルタをかけることで、例えば５００Ｈｚ未満の低周波数成分を除去し、この５００Ｈｚ未満の周波数成分を除去したデータを、ゴルフクラブヘッド３２の打音の評価に用いる振動波形データ（第１サンプル群）として用いる。このように、ハイパスフィルタをかけて低周波成分を除去することで、振動波形に含まれる、ゴルフスウィングにおけるゴルフクラブの風切り音などの余分な音の成分を除去することができる。
【００３２】
図４（ａ）および（ｂ）は、バイノーラル録音システム２０の騒音マイク１４で取得して、処理装置１８に記憶した、評価対象のゴルフクラブヘッド３２の打音の振動波形データを示している。図４（ａ）および（ｂ）は、トリガーレベルを０．１（Ｖ）、プリトリガー数値を２００（データ数）として規定された、打音の評価に用いる振動波形データを示し、時間軸上で最前部分の４０００個のデータを示している（打音の評価に用いるデータは、全体では１６３８４個である）。図４（ａ）は、ダミーヘッド１６の左耳に対応する側に設けられたマイク１４において取得された、打音の評価に用いる振動波形データ（左側振動波形データ）であり、（ｂ）は、ダミーヘッド１６の右耳に対応する側に設けられたマイク１４において取得された、打音の評価に用いる振動波形データ（右側振動波形データ）である。このように、トリガーレベルおよびプリトリガー値を、騒音マイクの特性などに応じて適切な値に設定することで、各ゴルフクラブによるゴルフボールの打撃のタイミングを確実に含んだ振動波形のデータを、評価データとして取得することができる。
【００３３】
次に、処理装置１８において、この振動波形データを用い、音圧パラメータを算出する（ステップＳ１０６）。この音圧パラメータ算出ステップは、オペレータが、ディスプレイ２６に表示されたプログラム実行制御画面３４を見ながら入力装置２８を操作して、音圧パラメータ計算実行指示ボックス３８ａを選択することで、実行が指示される。音圧パラメータとしては、ＩＳＯ５３２Ｂ規格に規定されるラウドネス（単位；ｓｏｎｅ)を対数表示した、公知のラウドネスレベル（単位；ｐｈｏｎ)の値を用いる。具体的には、打音の評価に用いる振動波形データ全体（１６３８４個）でＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍａｔｉｏｎ)による周波数解析を行なってスペクトル波形を求め、ＩＳＯ５３２Ｂ規格におけるラウドネス曲線を適用することで、振動波形のラウドネス（単位；ｓｏｎｅ）を求め、このラウドネスを対数表示したラウドネスレベル（単位；ｐｈｏｎ）の値を求める。処理装置１８では、このラウドネスレベルの値を、右側振動波形データおよび左側振動波形データそれぞれについて求める。そして、右側振動波形および左側振動波形それぞれのラウドネスの平均値を求め、この平均値を音圧パラメータの値とする。求められた音圧パラメータの値は、評価装置１８の図示しない記憶手段に記憶される。
【００３４】
次に、残響パラメータを算出する（ステップＳ１０８）。この残響パラメータ算出ステップは、オペレータが、ディスプレイ２６に表示されたプログラム実行制御画面３４を見ながら入力装置２８を操作して、音圧パラメータ計算実行指示ボックス３８ｂを選択することで実行が指示される。残響パラメータは、打音の減衰の度合いを表すパラメータである。残響パラメータの算出では、まず、振動波形データを、時間軸上で一番最初のサンプルデータ（始点データ）を基準番号とし、この番号から所定のサンプル数、例えば、サンプル数１００で音圧振動波形が順次区切られたブロック群αを作成する。そして、ブロック群αのブロックＡｋ毎に振動波形のレベルの二乗和を求め、各ブロック毎にこの二乗和の対数値（以降、ブロック値とする）を求める。そして、ブロック群αの各ブロックから、ブロック値が最大となる最大音圧ブロックを、起点ブロックとして抽出する。そして、この最大音圧ブロックのブロック値から所定値（本実施形態では３０ｄＢ）だけ低い値を終点ブロック値として設定し、起点ブロック以降の複数のブロックにおいて、ブロック値が最初に終点ブロック値以下となる終点ブロックを抽出する。そして、これら起点ブロックと終点ブロックそれぞれの、時間軸上で最前のサンプルデータ間での時間幅を残響パラメータとする。
【００３５】
詳述すると、４４，１ｋＨｚのサンプリング周波数でサンプリングされた音圧振動波形の信号をｙ（ｉ）（ｉ＝１〜１６３８４の自然数）とすると、下記式（１）に従ってサンプル数１００毎に区分けしたブロック数Ａｋ毎に、相対音圧レベルの値Ｅ（ｋ）（Ｋは１〜１６２８５）を計算し、例えば１６２８５個のＥ（ｋ）を求め、この中で最大となるＥ（ｋ）を最大音圧ブロックとして抽出する。図５（ｂ）は、図５（ａ）に示す時間軸について一部を拡大して示す時間軸において、各ブロックＡｋが占める時間範囲を示している。各ブロックＡｋに含まれるサンプルは、隣接するブロックと一部分がオーバーラップしている。図５（ａ）は、ブロック番号ｋに対するＥ（ｋ）の変化の一例を示すグラフである。図５(ａ）において、ブロック番号ＭでＥ（ｋ）の値が最大となっている。そして、ブロック番号Ｅで、Ｅ（ｋ）の値が、最大音圧ブロックＭのブロック値（Ｅ（ｋ）の値）から３０ｄＢだけ低い終点ブロック値を最初に下回っている。このブロック番号Ｍとブロック番号Ｅそれぞれのサンプル群における、最初のサンプル番号のサンプル間の時間幅の値を残響パラメータとして算出する。処理装置１８では、この残響パラメータの値を、右側振動波形データおよび左側振動波形データそれぞれについて求める。そして、右側振動波形および左側振動波形それぞれの残響パラメータの値の平均値を求め、この平均値を残響パラメータの値とする。求められた残響パラメータの値は、評価装置１８の図示しない記憶手段に記憶される。
【数１】
2006239132.tif 000003
【００３６】
次に、周波数パラメータを算出する（ステップＳ１１０）。この周波数パラメータ算出ステップは、オペレータが、ディスプレイ２６に表示されたプログラム実行制御画面３４を見ながら入力装置２８を操作して、周波数パラメータ計算実行指示ボックス３８ｃを選択することで実行が指示される。具体的には、得られた振動波形のデータについて、時間軸上で最前部分の所定個数のサンプル（本実施形態では、１５００個）を除いた振動波形データのうち、最前部分の所定数のデータ（本実施形態では、８１９２データ）について周波数解析を行なって、スペクトル波形を求める。そして、このスペクトル波形のスペクトル値を重み付け係数として用いた、低周波成分を除いた所定の周波数範囲（本実施形態では、１ｋＨｚ〜２０ｋＨｚ）における、振動波形のエネルギーの加重平均周波数を周波数パラメータとする。図６に示すグラフは、それぞれ図４（ａ）および（ｂ）に示す振動波形データのうち、時間軸上で最前部分の１５００個のサンプルを除いた振動波形のうち、時間軸上で最前部分の８１９２データについてＦＦＴ解析を行なって得られた、スペクトル波形データのうち０ｋＨｚ〜２０ｋＨｚの範囲のスペクトル波形である。周波数パラメータ算出ステップでは、これらのスペクトル波形のうち、１ｋＨｚ〜２０ｋＨｚの範囲（図６（ａ）および（ｂ）における周波数範囲αを除いた範囲）で、上記エネルギー加重平均周波数を求める。以下、周波数パラメータの算出について、より詳細に説明する。
【００３７】
エネルギー加重平均周波数は下記式（２）に従って求める。式（２）において、ｆは周波数、Ｌはスペクトル波形における周波数ｆにおけるスペクトル値から求めた、周波数ｆにおける振動波形の音圧レベル（ｄＢ値）であり、式（２）に従って、低周波成分を除いた所定の周波数範囲（１ｋＨｚ〜２０ｋＨｚ）においてエネルギー加重平均周波数を求める。このようなエネルギー加重平均周波数は、打撃後一定の時間を経た振動波形データを用いて求めるので、パラメータとしての精度、さらには再現性が悪化するといった問題は解消する。また、例えば周波数が１ｋＨｚ未満の低周波成分を除いた周波数範囲でエネルギー加重平均周波数を求めるので、振動波形について周波数解析を行なった際に発生するＤＣ成分の影響を除くことができ、パラメータとしての精度、さらには再現性が悪化するといった問題をさらに解消している。処理装置１８では、この周波数パラメータの値を、右側振動波形データおよび左側振動波形データそれぞれについて求める。そして、右側振動波形および左側振動波形それぞれの周波数パラメータの値の平均値を求め、この周波数パラメータの平均値を残響パラメータの値とする。求められた残響パラメータの値は、評価装置１８の図示しない記憶手段に記憶される。
【数２】
2006239132.tif 000004
【００３８】
次に、純音度パラメータを算出する（ステップＳ１１２）。この純音度パラメータ算出ステップは、オペレータが、ディスプレイ２６に表示されたプログラム実行制御画面３４を見ながら入力装置２８を操作して、純音度パラメータ計算実行指示ボックス３８ｄを選択することで実行が指示される。純音度パラメータの算出は、具体的には、得られた振動波形のデータについて、時間軸上で最前部分の所定個数のサンプル（本実施形態では、２５００個）を除いた振動波形データのうち、最前部分の所定数のデータ（本実施形態では、８１９２データ）について周波数解析を行なって、スペクトル波形を求める。そして、この求めたスペクトル波形のうち、低周波成分を除いた所定の周波数範囲（本実施形態では、１ｋＨｚ〜２０ｋＨｚ）における平均スペクトル値より高いスペクトル波形の第１の合計値と、この周波数範囲全体におけるスペクトル波形の第２の合計値とを求め、この第２の合計値に対する第１の合計値の比率（単位％）を純音度パラメータとして算出する。このような純温度パラメータは、打撃後一定の時間を経た振動波形データを用いて求めるので、パラメータとしての精度、さらには再現性が悪化するといった問題は解消する。また、例えば周波数が１ｋＨｚ未満の低周波成分を除いた周波数範囲でエネルギー加重平均周波数を求めるので、振動波形について周波数解析を行なった際に発生するＤＣ成分の影響を除くことができ、パラメータとしての精度、さらには再現性が悪化するといった問題をさらに解消している。処理装置１８では、この純音度パラメータの値を、右側振動波形データおよび左側振動波形データそれぞれについて求める。そして、右側振動波形および左側振動波形それぞれの純音度パラメータの値の平均値を求め、この純音度パラメータの平均値をゴルフクラブヘッド３２の打音の純音度パラメータの値とする。求められた純音度パラメータの値は、評価装置１８の図示しない記憶手段に記憶される。
【００３９】
次に、算出した４つの評価パラメータを用いて、任意の人物が、ゴルフクラブヘッド３２の打音を聞いた際に感じる、前記打音に対する好感の度合いの大きさを表す好感評価値を算出する（ステップＳ１１４）。この好感評価値算出ステップは、オペレータが、ディスプレイ２６に表示されたプログラム実行制御画面３４を見ながら入力装置２８を操作して、好感評価値計算実行指示ボックス３８ｅを選択することで開始される。好感評価値算出ステップが開始されると、評価装置１８の図示しない記憶手段に記憶された、上述の各評価パラメータと、予め記憶された各評価パラメータ毎の重み付け係数が読み出される。評価装置１８は、読み出された各評価パラメータ毎の重み付け係数で、各評価パラメータを重み付けして１つにまとめて好感評価値を算出する。この各評価パラメータ毎の重み付け係数は、予め１つのファイルにまとめられて評価装置１８に記憶されている。評価装置１８には、例えば騒音マイク１４の特性などに応じて、それぞれ異なる重み付け係数の値が記録された、複数種類のファイルが記録されていてもよい。複数種類のファイル名が記憶されている場合、プログラム実行制御画面３４の重み付け係数ファイル名入力ボックスに、所望の重み付け係数が記録されたファイル名を入力することで、好感評価値の算出に用いるための任意の重み付け係数を選択指示する。
【００４０】
具体的には、音圧パラメータをＸ１、残響パラメータをＸ２、周波数パラメータをＸ３、純音度パラメータをＸ４とすると、好感評価値Ｙを下記式（３）に従い算出する。
Ｙ１＝ａ０＋ａ１・Ｘ１＋ａ２・Ｘ２＋ａ３・Ｘ３＋ａ４・Ｘ４    （３）
ここで、ａ０、ａ１、ａ２、ａ３、ａ４は、上述のように評価装置１８に予め記憶されており、本実施形態においては、ａ０＝−１３５１．６４、ａ１＝８．７７１、ａ２＝−０．８８３、ａ３＝０．０１７、ａ４＝４．９２３の各数値が記憶されている。
【００４１】
図７（ａ）は、ゴルフクラブヘッドＡ〜Ｉの９種類のゴルフクラブヘッドの打音について、本実施形態において算出された好感評価値と、好感心理量Ｚ１との対応を示す相関図である。この好感心理量Ｚ１は、ゴルフクラブヘッドの打音を聞いた被験者が、打音に対して好ましいと感じる程度を表す心理量であり、打音が好ましいと感じるほど高い量である。官能評価の方法は、一対比較法を用いている。具体的には、ゴルフクラブヘッドＡ〜Ｉの９本のゴルフクラブヘッドの打音から２本のゴルフクラブヘッドの打音を選択し、２本のゴルフクラブヘッドの打音のうち、どちらのゴルフクラブヘッドの打音が好ましいと感じるか、被験者に選定評価してもらうといった評価を繰り返した。被験者は５４人であり、被験者１人１人に対し、９種類のゴルフクラブヘッドから２種類のゴルフクラブヘッドを選択した際の全ての組み合わせ（３６通り）について、上述の選定評価をしてもらい、回答結果全てを集計した。そして、この集計結果に基づき、各ゴルフクラブの打音それぞれについて、各ゴルフクラブそれぞれがより好ましいと選択される割合を求め、この割合から正規分布のＺ値を求めた。そして、９種類のゴルフクラブヘッドの打音のＺ値を母集団とし、各ゴルフクラブヘッドのＺ値の偏差値を求め、この偏差値を好感心理量Ｚ１とした。このような一対比較法では、被験者は、２種類のゴルフクラブヘッドの打音のうち、どちらのゴルフクラブヘッドの打音が好ましいと感じるかのみを選択すればよく、被験者が打音を聴いた際の直感的な感覚に基づいた回答が得られる。
【００４２】
なお、各ゴルフクラブヘッドの打音を被験者に聞かせる際は、ヘッドフォンを用い、被験者の左右の耳それぞれに、ダミーヘッドの対応する側の耳部分に備えたマイク１４で取得した打音情報を聞かせている。このようにして被験者に打音情報を聞かせることで、被験者は、被験者自身が実際にゴルフクラブを把持してゴルフスウィングを行なった際に、被験者自身に聞こえる打音と同様の音を聞くことができる。被験者は、このようにして各ゴルフクラブの打音を聴き、上述の官能評価を行なっている。このようにして求められた好感心理量Ｚ１は、任意のゴルファが、ゴルフクラブヘッドＡ〜Ｉを備えた９本のゴルフクラブそれぞれを把持してゴルフスウィングを行なって、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対する好感の程度を高精度に表していると言える。
【００４３】
図７（ａ）に示すように、打音に対する好感の程度を表す好感心理量Ｚ１と、上述の４つの評価パラメータを用いて算出した好感評価値との関係は直線によって回帰され、相関係数は０．８２９と高い数値となっている。４つの評価パラメータを重み付けして表した上述の好感評価値は、任意のゴルファが、所定のゴルフクラブヘッドを備えたゴルフクラブを把持してゴルフスウィングを行なって、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対する好感の程度を高精度に表していると言える。
【００４４】
本願発明は、上述の感応評価によって求められた好感心理量に基づき、この好感心理量を高精度に再現する、ゴルフクラブの打音の振動波形の物理的特徴量（評価パラメータ）の組み合わせ、およびこれら評価パラメータに対応する重み付け係数を、本願発明者が技術的考察に基づいて導出することで得られた発明であり、従来は定量的に評価することができなかった打音に対する好感の度合を、各評価パラメータを重み付けして得られる評価値に基づいて、定量的に評価するといった大きな効果を有する。
【００４５】
なお、この好感評価値算出ステップでは、上述の周波数パラメータＸ３および純音度パラメータＸ４のみを用い、ゴルフクラブヘッドＡ〜Ｉを備えた９種類のゴルフクラブそれぞれを把持してゴルフスウィングを行なって、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対して感じる音程の高さの程度を表す、高低評価値を併せて算出する。
【００４６】
具体的には、周波数パラメータをＸ３、純音度パラメータをＸ４を用い、高低評価値Ｙ２を下記式（４）を用いて算出する。
Ｙ２＝ｂ０＋ｂ１・Ｘ３＋ｂ２・Ｘ４    （４）
本実施形態においては、ｂ０＝−１３４．０３、ｂ１＝０．００１８５、ｂ２＝１．８９９５、の各数値が設定されている。
【００４７】
図７（ｂ）は、９種類のゴルフクラブヘッドＡ〜Ｉを備えるゴルフクラブそれぞれの打音について、本実施形態において算出された高低評価値Ｙ２と、高低心理量Ｚ２との対応を示す相関図である。この高低心理量Ｚ２は、ゴルフクラブヘッドの打音を聞いた被験者が、打音に対して感じる音程の高さの程度を表す心理量であり、打音の音程が高いと感じるほど高い量である。この高低心理量Ｚ２は、好感心理量Ｚ１と同様、一対評価法によって得られた、各ゴルフクラブヘッドの打音それぞれについて、各ゴルフクラブヘッドの打音それぞれがより好ましいと選択される割合に基づく、偏差値である。
図７（ｂ）に示すように、高低心理量Ｚ２と、周波数パラメータＸ３と純音度パラメータＸ４とを用いて算出した高低評価値Ｙ２との関係は直線によって回帰され、相関係数は０．７８８と高い数値となっている。２つの物理パラメータを重み付けして表した上述の高低評価値Ｙ２は、任意のゴルファが、所定のゴルフクラブヘッドを把持してゴルフスウィングを行なった際、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対して感じる音程の高さの程度を高精度に表していると言える。
【００４８】
さらに、この好感評価値算出ステップでは、上述の周波数パラメータＸ３と純音度パラメータＸ４のみを用い、ゴルフクラブヘッドＡ〜Ｉを備えた９本のゴルフクラブそれぞれを把持してゴルフスウィングを行なって、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対して感じる音の澄みの程度を表す、澄度評価値Ｙ３を併せて算出する。
【００４９】
具体的には、周波数パラメータＸ３、純音度パラメータＸ４を用い、澄度評価値Ｙ３を下記式（５）を用いて算出する。
Ｙ３＝ｃ０＋ｃ１・Ｘ３＋ｃ２・Ｘ４    （５）
ここで、ｃ０、ｃ１、ｃ２は、評価装置１８に予め記憶されている。本実施形態においては、ｃ０＝−１６６．６７３、ｃ１＝０．００５９６、ｃ２＝２．０３３６７の各数値が設定されている。
【００５０】
図７（ｃ）は、９種類のゴルフクラブヘッドＡ〜Ｉを備えるゴルフクラブそれぞれの打音について、本実施形態において算出された澄度評価値Ｙ３と、澄度心理量Ｚ３との対応を示す相関図である。この澄度心理量Ｚ３は、ゴルフクラブヘッドの打音を聞いた被験者が、打音に対して感じる音の澄みの程度を表す心理量であり、打音が澄んでいると感じるほど高い量である。この澄度心理量Ｚ３は、好感心理量Ｚ１と同様、一対評価法によって得られた、各ゴルフクラブヘッドの打音それぞれについて、各ゴルフクラブヘッドの打音それぞれがより澄んでいると選択される割合に基づく、偏差値である。
【００５１】
図７（ｃ）に示すように、澄度心理量Ｚ３と、周波数パラメータＸ３と純音度パラメータＸ４とを用いて算出した澄度評価値Ｙ３との関係は直線によって回帰され、相関係数は０．７６７と高い数値となっている。２つの物理パラメータを重み付けして表した上述の澄度評価値は、任意のゴルファが、所定のゴルフクラブヘッドを把持してゴルフスウィングを行なった際、ゴルファ自身がゴルフクラブヘッドの打音を聴いた際の、この打音に対して感じる音の澄みの程度を高精度に表していると言える。
【００５２】
本発明のゴルフクラブヘッドの音質評価方法では、このように、４つの物理パラメータから、任意の被験者が所定の打音を聴いた場合、被験者がこの打音に対して感じる好感の程度を表す好感評価値を求める。また、周波数パラメータＸ３と純音度パラメータＸ４との２つの物理パラメータから、任意の被験者が所定の打音を聴いた場合、被験者がこの打音に対して感じる、打音の音程の高さの度合を表す高低評価値を求める。また、同様に、周波数パラメータＸ３と純音度パラメータＸ４との２つの物理パラメータから、任意の被験者が所定の打音を聴いた場合、被験者がこの打音に対して感じる、この打音の澄みの程度を表す澄度評価値を求める。求めた各評価値は、例えば評価装置１８のディスプレイ２６に表示出力される。
【００５３】
そして、次に、算出した各評価値を用いて、ゴルフクラブヘッドの打音を評価する（ステップＳ１１６）。この評価ステップでは、例えば、予め設定された基準評価値範囲と、求められた好感評価値とを比較し、算出した好感評価値に対応する基準評価値範囲に応じて打音の好感の程度を判定して評価すればよい。例えば、標準的な好感の打音の範囲を表す基準評価範囲として、評価値４５〜５５の範囲を設定しておき、好感評価値が４５未満の打音を、好感が低い打音であると判定し、好感評価値が４５〜５５の打音を、標準的な好感の打音であると判定し、好感評価値が５５より大きい打音を、好感が高い打音であると判定すればよい。また、複数のゴルフクラブヘッドの打音について好感評価値をそれぞれ求め、それぞれの好感評価値に基づいて、各ゴルフクラブヘッド打音を相対評価してもよい。このような評価は、高低評価値および澄度評価値についても同様に行なえばよい。これらの判定結果は、例えば評価装置１８のディスプレイ２６に表示出力される。
【００５４】
なお、本発明のゴルフクラブヘッドの音質評価方法において求められる音圧パラメータの値は、任意の被験者が所定の打音を聴いた場合、この打音に対して感じる音の大きさの程度を表す音量心理量Ｚ４を高精度に表している。また、本発明のゴルフクラブヘッドの音質評価方法において求められる残響パラメータの値は、任意の被験者が所定の打音を聴いた場合、この打音に対して感じる、この音の響きの程度を表す、響き心理量Ｚ５を高精度に表している。
図８（ａ）および（ｂ）は、それぞれ、本発明のゴルフクラブヘッドの音質評価方法において求められた評価パラメータと、任意の被験者が所定の打音を聴いた場合に、この打音に対して感じる心理量の大きさとの関係を表すグラフである。図８（ａ）は、音圧パラメータＸ１と音量心理量Ｚ４との関係を表すグラフであり、図８（ｂ）は、残響パラメータＸ２と響き心理量Ｚ５との関係を表すグラフである。なお、ここで、各心理量は、上述の好感心理Ｚ１と同様、一対比較法で得られた集計結果に基づく偏差値である。
【００５５】
図８（ａ）に示すように、音量心理量Ｚ４と、音圧パラメータＸ１の関係は直線によって回帰され、相関係数は０．８３６と絶対値が１に近い値となっている。また、図８（ｂ）に示すように、響き心理量Ｚ５と、残響パラメータＸ２の関係は直線によって回帰され、相関係数は０．８０２と絶対値が１に近い値となっている。このように、音圧パラメータＸ１は、音量心理量Ｚ４を高い精度で表し、残響パラメータＸ２は、響き心理量Ｚ５を高い精度で表す。このように、音圧パラメータＸ１、残響パラメータＸ２、周波数パラメータＸ３、純音度パラメータＸ４は、それぞれ単独に、または所定の評価パラメータを組み合わせ、所定の重み付け係数で重み付けして１つにまとめることで、好感心理量Ｚ１、高低心理量Ｚ２、澄度心理量Ｚ３、音量心理量Ｚ４、響き心理量Ｚ５の各心理量を表す評価値を、それぞれ高精度に（各心理量と高い相関で）求めることができる。
【００５６】
本発明のゴルフクラブヘッドの音質評価方法によれば、ゴルフクラブヘッドの打音について４つの評価パラメータを算出し、この評価パラメータを用いて、評価対象のゴルフクラブヘッドの打音を被験者が聞いた場合の上記各心理量を表す評価値を算出することができる。このような評価値により、評価対象のゴルフクラブヘッドの打音が好まれる程度を、面倒な官能評価を実施することなく定量的に評価することができる。
これに加え、より好まれる打音のゴルフクラブヘッドを作製するには、どのような心理量を高くするようにすればよいかの指針が得られ、さらに、どのような物理パラメータが高くなるようにすればよいかの指針を得ることができる。本発明のゴルフクラブヘッドの音質評価方法は、ゴルフクラブの設計開発段階において効果的に活用することができる。
【００５７】
以上、本発明のゴルフクラブヘッドの打音評価方法、この方法をコンピュータに実行させるプログラムおよびゴルフクラブヘッドの打音評価システムについて詳細に説明したが、本発明は上記実施例には限定はされず、本発明の要旨を逸脱しない範囲において、各種の改良および変更を行ってもよいのはもちろんである。 
【図面の簡単な説明】
【００５８】
【図１】（ａ）および（ｂ）は、本発明のゴルフクラブヘッドの打音評価システムの一例の構成を説明する図であり、（ａ）は概略斜視図、（ｂ）は概略平面図である。
【図２】本発明のゴルフクラブヘッドの打音評価方法の一例のフローチャート図である。
【図３】本発明の打音評価プログラムの一例である音質評価プログラムを実行した際、評価装置のディスプレイに表示される画面の一例を示す図である。
【図４】（ａ）および（ｂ）は、本発明のゴルフクラブヘッドの打音評価システムで取得される振動波形データの一部をそれぞれ示している。
【図５】（ａ）は、図４（ａ）に示す振動波形データをブロック化した際のブロック値の変化の一例を示すグラフであり、（ｂ）は、（ａ）に示す時間軸において、各ブロックが占める時間範囲を説明する図である。
【図６】（ａ）および（ｂ）は、本発明のゴルフクラブの打音評価方法で求められるスペクトル波形の一例を示す図である。
【図７】（ａ）〜（ｃ）は、それぞれ、本発明のゴルフクラブの打音評価方法で算出された評価値と心理量との対応を示す相関図である。
【図８】（ａ）および（ｂ）は、それぞれ、本発明のゴルフクラブの打音評価方法において算出された評価パラメータの値と心理量との対応を示す相関図である。
【符号の説明】
【００５９】
１０  ゴルフクラブヘッドの打音評価システム
１２  スウィングロボット
１３  ゴルフボール
１４  騒音マイク
１６  ダミーヘッド
１８  処理装置
２０  バイノーラル録音システム
２６  表示手段
２８  入力手段
３０  ゴルフクラブ
３２  ゴルフクラブヘッド
３４  プログラム実行制御画面
３６ａ〜３６ｃ  入力ボックス
３８ａ〜３８ｅ  計算実行指示ボタン
【図１】
2006239132.tif 000005
【図２】
2006239132.tif 000006
【図３】
2006239132.tif 000007
【図４】
2006239132.tif 000008
【図５】
2006239132.tif 000009
【図６】
2006239132.tif 000010
【図７】
2006239132.tif 000011
【図８】
2006239132.tif 000012