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