JP3146691B2 - 音色評価方法及びその装置 - Google Patents
音色評価方法及びその装置Info
- Publication number
- JP3146691B2 JP3146691B2 JP28938892A JP28938892A JP3146691B2 JP 3146691 B2 JP3146691 B2 JP 3146691B2 JP 28938892 A JP28938892 A JP 28938892A JP 28938892 A JP28938892 A JP 28938892A JP 3146691 B2 JP3146691 B2 JP 3146691B2
- Authority
- JP
- Japan
- Prior art keywords
- amplitude
- frequency
- modulated wave
- sound
- fluctuation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、振幅変動音の人間の聴
感に与える変動感あるいは心地よさの感覚を適切に評価
する音色評価方法及びその装置に関する。
感に与える変動感あるいは心地よさの感覚を適切に評価
する音色評価方法及びその装置に関する。
【０００２】
【従来の技術】日常、人間が接する音のほとんどは、そ
の音圧振幅が時間的に変動している。こういった音は、
その変動の仕方によって、不快感を催したり心地よい音
になったりする。従来、このような振幅変動音を評価す
る方法として、音の各周波数成分の変動量を測定し、変
動の音圧レベル（音圧振幅）の大きさに応じた聴感補正
を行った後、評価値を算出する方法が提案されている。
（例えば特開平３−８２９２２）しかし、上記の方法は
人間の聴感の特性を考慮しているとはいえない。たとえ
ば１kHz の音の 音圧振幅が３０Hzで変動している場合
と、１kHzの音の音圧振幅が３００Hzで変動している場
合を比べると、たとえ変動の音圧レベル（音圧振幅）の
大きさが同じであっても、前者のほうが後者より変動感
は強く、後者の音はほとんど変動感を感じない。それに
もかかわらず上記の従来技術では、両者は同じ補正が行
なわれるため聴感にあった評価を行うことができない。
の音圧振幅が時間的に変動している。こういった音は、
その変動の仕方によって、不快感を催したり心地よい音
になったりする。従来、このような振幅変動音を評価す
る方法として、音の各周波数成分の変動量を測定し、変
動の音圧レベル（音圧振幅）の大きさに応じた聴感補正
を行った後、評価値を算出する方法が提案されている。
（例えば特開平３−８２９２２）しかし、上記の方法は
人間の聴感の特性を考慮しているとはいえない。たとえ
ば１kHz の音の 音圧振幅が３０Hzで変動している場合
と、１kHzの音の音圧振幅が３００Hzで変動している場
合を比べると、たとえ変動の音圧レベル（音圧振幅）の
大きさが同じであっても、前者のほうが後者より変動感
は強く、後者の音はほとんど変動感を感じない。それに
もかかわらず上記の従来技術では、両者は同じ補正が行
なわれるため聴感にあった評価を行うことができない。
【０００３】
［着眼点］振幅変動音に対する人間の聴感は、振幅変調
波の特徴量に対して極値もつような特性を示す。例えば
図３に破線で示す音圧波形の包絡線によって定義される
振幅変調波の周波数に関して例を挙げれば、１kHzの音
の振幅変調波の周波数が３Hzである場合と、３０Hzであ
る 場合と、３００Hzである場合とを比べると、たとえ
３者の変動の大きさが同じであっても、人間の聴感上は
３０Hzの音の変動感が最も強く、３Hzおよび３００Hzの
音はほとんど変動感を感じない。また、前記振幅変調波
の周波数ｆの時間的ゆらぎに関して考えた場合、図１０
に示す周波数ｆの時間的ゆらぎの周波数分布（周波数毎
の頻度）がゆらぎの周波数ｆに対して、１／ｆ0.5の 音
と、１／ｆの音と、１／ｆ2の音とを比べると、ゆらぎ
の周波数分布が１／ｆの音が聴感上最も心地よい。本発
明は、上記のような人間の聴感上の特性に着眼し、音の
振幅変調波の特徴量に対して極値を持つような聴感特性
係数を用いることにより、人間の聴感にあった評価を行
なうことを目的とする。
波の特徴量に対して極値もつような特性を示す。例えば
図３に破線で示す音圧波形の包絡線によって定義される
振幅変調波の周波数に関して例を挙げれば、１kHzの音
の振幅変調波の周波数が３Hzである場合と、３０Hzであ
る 場合と、３００Hzである場合とを比べると、たとえ
３者の変動の大きさが同じであっても、人間の聴感上は
３０Hzの音の変動感が最も強く、３Hzおよび３００Hzの
音はほとんど変動感を感じない。また、前記振幅変調波
の周波数ｆの時間的ゆらぎに関して考えた場合、図１０
に示す周波数ｆの時間的ゆらぎの周波数分布（周波数毎
の頻度）がゆらぎの周波数ｆに対して、１／ｆ0.5の 音
と、１／ｆの音と、１／ｆ2の音とを比べると、ゆらぎ
の周波数分布が１／ｆの音が聴感上最も心地よい。本発
明は、上記のような人間の聴感上の特性に着眼し、音の
振幅変調波の特徴量に対して極値を持つような聴感特性
係数を用いることにより、人間の聴感にあった評価を行
なうことを目的とする。
【０００４】
【課題を解決するための手段】請求項１に記載の本発明
の第１発明は、評価対象とする音を電気信号に変換する
工程と、該周波数帯域成分の信号波形の包絡線である振
幅変調波を検出する工程と、前記振幅変調波を分析し振
幅変調波の周波数毎の振幅変動量を演算する工程と、該
振幅変動量と、前記振幅変調波の特徴量に対して極値を
持つ聴感特性係数とに基づき評価値を演算する工程と、
該評価値に基づき音色を評価する工程とから成ることを
特徴とする音色評価方法である。
の第１発明は、評価対象とする音を電気信号に変換する
工程と、該周波数帯域成分の信号波形の包絡線である振
幅変調波を検出する工程と、前記振幅変調波を分析し振
幅変調波の周波数毎の振幅変動量を演算する工程と、該
振幅変動量と、前記振幅変調波の特徴量に対して極値を
持つ聴感特性係数とに基づき評価値を演算する工程と、
該評価値に基づき音色を評価する工程とから成ることを
特徴とする音色評価方法である。
【０００５】
【作用】本第１発明は、以下に述べるような作用を奏す
る。まず、評価対象とする音を電気信号に変換し、該電
気信号波形の包絡線である振幅変調波を検出する。次
に、前記振幅変調波を分析し振幅変調波の周波数毎の振
幅変動量を演算する。さらに、該振幅変動量と、前記振
幅変調波の特徴量に対して極値を持つ聴感特性係数とに
基づき評価値を演算し、該評価値に基づき音色を評価す
る。
る。まず、評価対象とする音を電気信号に変換し、該電
気信号波形の包絡線である振幅変調波を検出する。次
に、前記振幅変調波を分析し振幅変調波の周波数毎の振
幅変動量を演算する。さらに、該振幅変動量と、前記振
幅変調波の特徴量に対して極値を持つ聴感特性係数とに
基づき評価値を演算し、該評価値に基づき音色を評価す
る。
【０００６】
【発明の効果】本第１発明は、前記振幅変調波の特徴量
に対して極値をもつような補正を行なうことによって、
音の振幅変調波の特徴量が聴感上の音色に強く影響する
という人間の聴感上の特性に着眼し、音の振幅変調波の
特徴量に対して極値を持つような聴感特性係数を用いる
ことにより、人間の聴感にあった音色評価を行なうこと
が可能になる。
に対して極値をもつような補正を行なうことによって、
音の振幅変調波の特徴量が聴感上の音色に強く影響する
という人間の聴感上の特性に着眼し、音の振幅変調波の
特徴量に対して極値を持つような聴感特性係数を用いる
ことにより、人間の聴感にあった音色評価を行なうこと
が可能になる。
【０００７】（その他の発明）請求項２に記載の本発明
の第２発明は、評価対象とする音を収音し電気信号に変
換するマイクロホンと、該電気信号波形の包絡線である
振幅変調波を検出する振幅変調波検出部と、前記振幅変
調波を分析し振幅変調波の周波数毎の振幅変動量を演算
する振幅変動量演算部と、該振幅変動量と該振幅変調波
の特徴量に対して極値を持つ聴感特性係数とに基づき評
価値を演算する評価値演算部と、該評価値算出部の出力
に基づき音色を評価する音色評価部とから成ることを特
徴とする音色評価装置である。本第２発明は、前記第１
発明と同様の作用・効果を奏する。
の第２発明は、評価対象とする音を収音し電気信号に変
換するマイクロホンと、該電気信号波形の包絡線である
振幅変調波を検出する振幅変調波検出部と、前記振幅変
調波を分析し振幅変調波の周波数毎の振幅変動量を演算
する振幅変動量演算部と、該振幅変動量と該振幅変調波
の特徴量に対して極値を持つ聴感特性係数とに基づき評
価値を演算する評価値演算部と、該評価値算出部の出力
に基づき音色を評価する音色評価部とから成ることを特
徴とする音色評価装置である。本第２発明は、前記第１
発明と同様の作用・効果を奏する。
【０００８】
（第１実施例）本発明の第１実施例として、自動車など
の機械から発生し、周期的に振幅が変動する音の変動感
を評価する音色評価装置に関して説明する。図1は、本
発明の第１実施例にかかわる音色評価装置の全体構成を
示す。マイクロホン１は、評価対象とする音を収録し電
気信号に変換する。Ａ／Ｄ変換器８は、前記電気信号を
デジタル信号に変換する。Ａ／Ｄ変換器８及び以下に説
明するデータ記憶部９、帯域フィルタ２、振幅変調波検
出部３、振幅変調波周波数分析部４、評価値演算部５、
音色評価部６は、パーソナルコンピュータで構成され、
表示部７は該パーソナルコンピュータに接続されたディ
スプレイ装置で構成される。
の機械から発生し、周期的に振幅が変動する音の変動感
を評価する音色評価装置に関して説明する。図1は、本
発明の第１実施例にかかわる音色評価装置の全体構成を
示す。マイクロホン１は、評価対象とする音を収録し電
気信号に変換する。Ａ／Ｄ変換器８は、前記電気信号を
デジタル信号に変換する。Ａ／Ｄ変換器８及び以下に説
明するデータ記憶部９、帯域フィルタ２、振幅変調波検
出部３、振幅変調波周波数分析部４、評価値演算部５、
音色評価部６は、パーソナルコンピュータで構成され、
表示部７は該パーソナルコンピュータに接続されたディ
スプレイ装置で構成される。
【０００９】データ記憶部９はＡ／Ｄ変換器８から出力
されるデジタル信号を記憶する。帯域フィルタ２は、デ
ータ記憶部９に保存されたデジタル信号から各周波数帯
域成分の信号を抽出するものである。ここで、帯域通過
フィルタ２の中心周波数および帯域幅は人間の聴覚系の
特性である臨界帯域（E.Zwicker:Subdivision of theau
dible frequency range into critical bands",J. Acou
st.Soc.Am. 33,248 (1961)）に示されるよう に、中心
周波数５０、１５０、２５０、・・・、１３．５kHz、
帯域幅８０、１００、１０００、・・・・、３５００Hz
の２４の帯域フィルタである。振幅変調波検出部３は、
帯域フィルタ２からの出力波形（図４に実線で示す）の
包絡線（図４に破線で示す）である振幅変調波を検出す
る。振幅変調波周波数分析部４は振幅変調波をＦＦＴ分
析し、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）を算
出する。評価値演算部５は、振幅変調波の周波数毎の振
幅変動量Ｄ（ｆ）と振幅変調波の周波数に応じた聴感特
性係数Ｃ（ｆ）より次式のように評価値Ｅを求める。 Ｅ＝∫Ｄ（ｆ）・Ｃ（ｆ）ｄｆ ・・・ （１）
されるデジタル信号を記憶する。帯域フィルタ２は、デ
ータ記憶部９に保存されたデジタル信号から各周波数帯
域成分の信号を抽出するものである。ここで、帯域通過
フィルタ２の中心周波数および帯域幅は人間の聴覚系の
特性である臨界帯域（E.Zwicker:Subdivision of theau
dible frequency range into critical bands",J. Acou
st.Soc.Am. 33,248 (1961)）に示されるよう に、中心
周波数５０、１５０、２５０、・・・、１３．５kHz、
帯域幅８０、１００、１０００、・・・・、３５００Hz
の２４の帯域フィルタである。振幅変調波検出部３は、
帯域フィルタ２からの出力波形（図４に実線で示す）の
包絡線（図４に破線で示す）である振幅変調波を検出す
る。振幅変調波周波数分析部４は振幅変調波をＦＦＴ分
析し、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）を算
出する。評価値演算部５は、振幅変調波の周波数毎の振
幅変動量Ｄ（ｆ）と振幅変調波の周波数に応じた聴感特
性係数Ｃ（ｆ）より次式のように評価値Ｅを求める。 Ｅ＝∫Ｄ（ｆ）・Ｃ（ｆ）ｄｆ ・・・ （１）
【００１０】上記聴感特性係数Ｃ（ｆ）はあらかじめ官
能評価実験などにより求められるものであり、振幅変調
波の周波数に対して極値を持つ特性となる。例えば帯域
フィルタの中心周波数が２５０Hzの周波数帯域での聴感
特性係数Ｃ（ｆ）を図２に、４kHzの周波数帯域での聴
感特性係数Ｃ（ｆ）を図３に示す。ここで特徴的なの
は、振幅変動の特徴量として振幅変調波の周波数に着目
し、振幅変調波の周波数に対して極値をもつような特性
を用いることである。これによって、「発明が解決しよ
うとする課題」の「着眼点」で説明しているように、１
kHzの音の振幅変調波の周波数が３Hzである場合と、３
０Hzである場合と、３００Hzで ある場合とを比べる
と、たとえ３者の変動の大きさが同じであっても、人間
の聴感上は３０Hzの音の変動感が最も強く、３Hzおよび
３００Hzの音はほとんど変動感を感じないという人間の
聴感上の特性に着眼し、音の振幅変動の特徴量に対して
極値を持つような聴感補正を行うことにより、人間の聴
感にあった評価を行なうことが可能になる。さらに、評
価値演算部の出力は、音色評価部６で加算して集計され
音色評価値となり、その結果が表示部７で表示される。
能評価実験などにより求められるものであり、振幅変調
波の周波数に対して極値を持つ特性となる。例えば帯域
フィルタの中心周波数が２５０Hzの周波数帯域での聴感
特性係数Ｃ（ｆ）を図２に、４kHzの周波数帯域での聴
感特性係数Ｃ（ｆ）を図３に示す。ここで特徴的なの
は、振幅変動の特徴量として振幅変調波の周波数に着目
し、振幅変調波の周波数に対して極値をもつような特性
を用いることである。これによって、「発明が解決しよ
うとする課題」の「着眼点」で説明しているように、１
kHzの音の振幅変調波の周波数が３Hzである場合と、３
０Hzである場合と、３００Hzで ある場合とを比べる
と、たとえ３者の変動の大きさが同じであっても、人間
の聴感上は３０Hzの音の変動感が最も強く、３Hzおよび
３００Hzの音はほとんど変動感を感じないという人間の
聴感上の特性に着眼し、音の振幅変動の特徴量に対して
極値を持つような聴感補正を行うことにより、人間の聴
感にあった評価を行なうことが可能になる。さらに、評
価値演算部の出力は、音色評価部６で加算して集計され
音色評価値となり、その結果が表示部７で表示される。
【００１１】次に上記構成における作用を説明する。本
第１実施例のフローチャートを図１１に示す。評価対象
音はマイクロホン１により、電気信号に変換される。該
電気信号はＡ／Ｄ変換器８によってデジタル信号に変換
され、その信号のうち２００msec程度の長さの信号がデ
ータ記憶部９において保存され、さらに帯域通過フィル
タ２により各帯域に分割される。この出力信号の例を図
４に示す。次にその出力信号が振幅変調波検出部３にお
いて振幅変調波に変換される。この振幅変調波の例を図
５に示す。次に、振幅変調波周波数分析部４において振
幅変調波の周波数毎の振幅変動量Ｄ（ｆ）が算出され
る。Ｄ（ｆ）の例を図６に示す。さらに、評価値算出部
５において、振幅変調波の周波数毎の振幅変動量Ｄ
（ｆ）と振幅変調波の周波数に応じた聴感特性係数Ｃ
（ｆ）より前記（１）式により評価値Ｅを求める。
第１実施例のフローチャートを図１１に示す。評価対象
音はマイクロホン１により、電気信号に変換される。該
電気信号はＡ／Ｄ変換器８によってデジタル信号に変換
され、その信号のうち２００msec程度の長さの信号がデ
ータ記憶部９において保存され、さらに帯域通過フィル
タ２により各帯域に分割される。この出力信号の例を図
４に示す。次にその出力信号が振幅変調波検出部３にお
いて振幅変調波に変換される。この振幅変調波の例を図
５に示す。次に、振幅変調波周波数分析部４において振
幅変調波の周波数毎の振幅変動量Ｄ（ｆ）が算出され
る。Ｄ（ｆ）の例を図６に示す。さらに、評価値算出部
５において、振幅変調波の周波数毎の振幅変動量Ｄ
（ｆ）と振幅変調波の周波数に応じた聴感特性係数Ｃ
（ｆ）より前記（１）式により評価値Ｅを求める。
【００１２】ここで特徴的なのは、振幅変調波の特徴量
として振幅変調波の周波数に着目し、振幅変調波の周波
数に対して極値をもつような特性を用いることである。
これによって、「発明が解決しようとする課題」の「着
眼点」で説明しているように、１kHzの音の振幅変調波
の周波数が３Hzである場合と、３０Hzである場合と、３
００Hzで ある場合とを比べると、たとえ３者の変動の
大きさが同じであっても、人間の聴感上は３０Hzの音の
変動感が最も強く、３Hzおよび３００Hzの音はほとんど
変動感を感じないという人間の聴感上の特性に着眼し、
音の振幅変動の特徴量に対して極値を持つような特性係
数を用いることにより、人間の聴感にあった評価を行な
うことが可能になる。そして、評価値算出部５の出力
は、集計部６で加算して集計され、その結果が表示部７
で表示される。
として振幅変調波の周波数に着目し、振幅変調波の周波
数に対して極値をもつような特性を用いることである。
これによって、「発明が解決しようとする課題」の「着
眼点」で説明しているように、１kHzの音の振幅変調波
の周波数が３Hzである場合と、３０Hzである場合と、３
００Hzで ある場合とを比べると、たとえ３者の変動の
大きさが同じであっても、人間の聴感上は３０Hzの音の
変動感が最も強く、３Hzおよび３００Hzの音はほとんど
変動感を感じないという人間の聴感上の特性に着眼し、
音の振幅変動の特徴量に対して極値を持つような特性係
数を用いることにより、人間の聴感にあった評価を行な
うことが可能になる。そして、評価値算出部５の出力
は、集計部６で加算して集計され、その結果が表示部７
で表示される。
【００１３】本第１実施例の音色評価装置を用いて自動
車の排気音の評価を行なった。車種の異なる１０種類の
排気音の変動感の評価実験をした結果を図７に示す。縦
軸は本装置による変動感の評価値の時間的平均値で、横
軸は専門家による変動感の官能評価値である。図７より
本装置による変動感の評価値の時間的平均値と、専門家
による変動感の官能評価値との相関が強いことが分か
る。以上のように、本第１実施例の音色評価装置は、振
幅変動の特徴量として振幅変調波の周波数に着目し、振
幅変調波の周波数に対して極値をもつような特性係数を
用いることによって人間の聴感と合致した評価値を得る
ことができる。
車の排気音の評価を行なった。車種の異なる１０種類の
排気音の変動感の評価実験をした結果を図７に示す。縦
軸は本装置による変動感の評価値の時間的平均値で、横
軸は専門家による変動感の官能評価値である。図７より
本装置による変動感の評価値の時間的平均値と、専門家
による変動感の官能評価値との相関が強いことが分か
る。以上のように、本第１実施例の音色評価装置は、振
幅変動の特徴量として振幅変調波の周波数に着目し、振
幅変調波の周波数に対して極値をもつような特性係数を
用いることによって人間の聴感と合致した評価値を得る
ことができる。
【００１４】（第２実施例）本発明の第２実施例とし
て、音のゆらぎの心地よさを評価する音色評価装置に関
して説明する。第８図は本発明の第２実施例にかかわる
音色評価装置の全体構成を示す。マイクロホン１１は評
価対象とする音を収録し、電気信号に変換する。Ａ／Ｄ
変換器１９は、前記電気信号をデジタル信号に変換す
る。Ａ／Ｄ変換器１９及び以下に説明するデータ記憶部
２０、帯域フィルタ１２、振幅変調波検出部１３、振幅
変調波周波数分析部１４、ゆらぎ分布指標演算部１５、
評価値演算部１６、音色評価部１７は、パーソナルコン
ピュータで構成され、表示部１８は該パーソナルコンピ
ュータに接続されたディスプレイ装置で構成される。
て、音のゆらぎの心地よさを評価する音色評価装置に関
して説明する。第８図は本発明の第２実施例にかかわる
音色評価装置の全体構成を示す。マイクロホン１１は評
価対象とする音を収録し、電気信号に変換する。Ａ／Ｄ
変換器１９は、前記電気信号をデジタル信号に変換す
る。Ａ／Ｄ変換器１９及び以下に説明するデータ記憶部
２０、帯域フィルタ１２、振幅変調波検出部１３、振幅
変調波周波数分析部１４、ゆらぎ分布指標演算部１５、
評価値演算部１６、音色評価部１７は、パーソナルコン
ピュータで構成され、表示部１８は該パーソナルコンピ
ュータに接続されたディスプレイ装置で構成される。
【００１５】データ記憶部２０は、Ａ／Ｄ変換器１９か
ら出力されるデジタル信号を記憶する。帯域フィルタ１
２は、データ記憶部９に保存されたデジタル信号から各
周波数帯域成分の信号を抽出するものである。ここで、
帯域フィルタ１２の中心周波数および帯域幅は前記第１
実施例と同様である。振幅変調波検出部１３は、帯域フ
ィルタ１２からの出力波形（図４に実線で示す）の包絡
線（図４に破線で示す）である振幅変調波を検出する。
振幅変調波周波数分析部１４は振幅変調波をＦＦＴ分析
し、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）を算出
する。ゆらぎ分布指数演算部１５は、振幅変調波の周波
数ｆ毎の振幅変動量Ｄ（ｆ）よりａ及びｂを（２）式で
求める。 Ｄ（ｆ）＝ｂ／（ｆa） ・・・ （２） 具体的には、（２）式の両辺の対数をとり、最小二乗法
によりａ、ｂを定める。前記ａは、ゆらぎ分布指標であ
る。評価値算出部１６は、ゆらぎ分布指標ａに応じた係
数Ｐ（ａ）より次式のように評価値Ｅを求める。 Ｅ＝Ｐ（ａ） ・・・ （３） 上記係数Ｐ（ａ）はあらかじめ官能評価実験などにより
求められるものであり、ａ＝１のとき極値を持つ特性と
なる。図９に例を示す。さらに、音色評価部１７は、評
価値演算部１６の出力を加算し集計して音色評価値と
し、表示部１８でその結果が表示される。
ら出力されるデジタル信号を記憶する。帯域フィルタ１
２は、データ記憶部９に保存されたデジタル信号から各
周波数帯域成分の信号を抽出するものである。ここで、
帯域フィルタ１２の中心周波数および帯域幅は前記第１
実施例と同様である。振幅変調波検出部１３は、帯域フ
ィルタ１２からの出力波形（図４に実線で示す）の包絡
線（図４に破線で示す）である振幅変調波を検出する。
振幅変調波周波数分析部１４は振幅変調波をＦＦＴ分析
し、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）を算出
する。ゆらぎ分布指数演算部１５は、振幅変調波の周波
数ｆ毎の振幅変動量Ｄ（ｆ）よりａ及びｂを（２）式で
求める。 Ｄ（ｆ）＝ｂ／（ｆa） ・・・ （２） 具体的には、（２）式の両辺の対数をとり、最小二乗法
によりａ、ｂを定める。前記ａは、ゆらぎ分布指標であ
る。評価値算出部１６は、ゆらぎ分布指標ａに応じた係
数Ｐ（ａ）より次式のように評価値Ｅを求める。 Ｅ＝Ｐ（ａ） ・・・ （３） 上記係数Ｐ（ａ）はあらかじめ官能評価実験などにより
求められるものであり、ａ＝１のとき極値を持つ特性と
なる。図９に例を示す。さらに、音色評価部１７は、評
価値演算部１６の出力を加算し集計して音色評価値と
し、表示部１８でその結果が表示される。
【００１６】ここで上記構成における作用を説明する。
本第２実施例のフローチャートを図１２に示す。評価対
象音はマイクロホン１１により、電気信号に変換され
る。その出力信号は、Ａ／Ｄ変換器１９でデジタル信号
に変換され、その信号のうち５分程度の長さの信号がデ
ータ記憶部２０に記憶される。データ記憶部２０に記憶
された信号は、帯域フィルタ１２によって各帯域に分割
され、さらに振幅変調波検出部１３において振幅変調波
に変換される。その振幅変調波より振幅変動量演算部１
４において振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）
が算出される。さらに、ゆらぎ分布指標演算部１５にお
いて、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）より
ゆらぎ分布指標ａが算出される。次に、評価値演算部１
６でゆらぎ分布指標ａより評価値Ｅが前記（２）式によ
り求まる。
本第２実施例のフローチャートを図１２に示す。評価対
象音はマイクロホン１１により、電気信号に変換され
る。その出力信号は、Ａ／Ｄ変換器１９でデジタル信号
に変換され、その信号のうち５分程度の長さの信号がデ
ータ記憶部２０に記憶される。データ記憶部２０に記憶
された信号は、帯域フィルタ１２によって各帯域に分割
され、さらに振幅変調波検出部１３において振幅変調波
に変換される。その振幅変調波より振幅変動量演算部１
４において振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）
が算出される。さらに、ゆらぎ分布指標演算部１５にお
いて、振幅変調波の周波数毎の振幅変動量Ｄ（ｆ）より
ゆらぎ分布指標ａが算出される。次に、評価値演算部１
６でゆらぎ分布指標ａより評価値Ｅが前記（２）式によ
り求まる。
【００１７】ここで特徴的なのは、振幅変調波の特徴量
として振幅変調波の周波数のゆらぎ分布指標に着目し、
振幅変調波の周波数のゆらぎ分布指標に対して極値をも
つような特性を用いることである。図９にゆらぎ分布指
標ａに対する補正係数Ｐ（ａ）の関係を示す。これによ
って、「発明が解決しようとする課題」の「着眼点」で
説明しているように、振幅変調波の周波数ｆの時間的ゆ
らぎに関して考えた場合、該周波数ｆの時間的ゆらぎの
周波数分布（周波数毎の頻度）がゆらぎの周波数ｆに対
して、１／ｆ0.5の音と、１／ｆ の音と、１／ｆ2の音
とを比べると、たとえ３者の変動の大きさが同じであっ
ても、ゆらぎの周波数分布が１／ｆの音が聴感上最も心
地よいという人間の聴感上の特性を考慮した評価が可能
となる。さらに、評価値演算部１６の出力は、音色評価
部１７で加算して集計され音色評価値となり、その結果
が表示部１８で表示される。以上のように、本第２実施
例によれば、人間の聴感上の心地よさと良く対応がとれ
た適切な音色評価を行なうことができる。
として振幅変調波の周波数のゆらぎ分布指標に着目し、
振幅変調波の周波数のゆらぎ分布指標に対して極値をも
つような特性を用いることである。図９にゆらぎ分布指
標ａに対する補正係数Ｐ（ａ）の関係を示す。これによ
って、「発明が解決しようとする課題」の「着眼点」で
説明しているように、振幅変調波の周波数ｆの時間的ゆ
らぎに関して考えた場合、該周波数ｆの時間的ゆらぎの
周波数分布（周波数毎の頻度）がゆらぎの周波数ｆに対
して、１／ｆ0.5の音と、１／ｆ の音と、１／ｆ2の音
とを比べると、たとえ３者の変動の大きさが同じであっ
ても、ゆらぎの周波数分布が１／ｆの音が聴感上最も心
地よいという人間の聴感上の特性を考慮した評価が可能
となる。さらに、評価値演算部１６の出力は、音色評価
部１７で加算して集計され音色評価値となり、その結果
が表示部１８で表示される。以上のように、本第２実施
例によれば、人間の聴感上の心地よさと良く対応がとれ
た適切な音色評価を行なうことができる。
【図１】本発明の第１実施例の全体構成を示す図
【図２】本発明の第１実施例の振幅変調波の周波数ｆと
補正係数Ｃ（ｆ）との関係を表す線図（帯域フィルタの
中心周波数２５０Hz）
補正係数Ｃ（ｆ）との関係を表す線図（帯域フィルタの
中心周波数２５０Hz）
【図３】本発明の第１実施例の振幅変調波の周波数ｆと
補正係数Ｃ（ｆ）との関係を表す線図（帯域フィルタの
中心周波数４kHz）
補正係数Ｃ（ｆ）との関係を表す線図（帯域フィルタの
中心周波数４kHz）
【図４】本発明の第１および第２実施例の評価対象とす
る音を電気信号に変換した波形を表す線図
る音を電気信号に変換した波形を表す線図
【図５】本発明の第１および第２実施例の振幅変調波の
波形を表す線図
波形を表す線図
【図６】本発明の第１および第２実施例の振幅変調波の
周波数ｆと振幅変動量Ｄ（ｆ）との関係を表す線図
周波数ｆと振幅変動量Ｄ（ｆ）との関係を表す線図
【図７】本発明の第１実施例の音色評価値と官能評価値
との相関を表す線図
との相関を表す線図
【図８】本発明の第２実施例の全体構成を示す図
【図９】本発明の第２実施例のゆらぎ分布指標ａと補正
係数Ｐ（ａ）との関係を表す線図
係数Ｐ（ａ）との関係を表す線図
【図１０】本発明の第２実施例のゆらぎ分布指標ａの定
義を表す線図
義を表す線図
【図１１】本発明の第１実施例のフローチャート
【図１２】本発明の第２実施例のフローチャート
１ マイクロホン ２ 帯域フィルタ ３ 振幅変調波検出部 ４ 振幅変動量演算部 ５ 評価値演算部 ６ 音色評価部 ７ 表示部 ８ Ａ／Ｄ変換器
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7，ＤＢ名) G01H 3/00 G01H 17/00
Claims (2)
- 【請求項１】 評価対象とする音を電気信号に変換する
工程と、 該電気信号波形の包絡線である振幅変調波を検出する工
程と、 前記振幅変調波を分析し振幅変調波の周波数毎の振幅変
動量を演算する工程と、 該振幅変動量と、前記振幅変調波の特徴量に対して極値
を持つ聴感特性係数とに基づき評価値を演算する工程
と、 該評価値に基づき音色を評価する工程と、から成ること
を特徴とする音色評価方法。 - 【請求項２】 評価対象とする音を収音し電気信号に変
換するマイクロホンと、 該電気信号波形の包絡線である振幅変調波を検出する振
幅変調波検出部と、 前記振幅変調波を分析し振幅変調波の周波数毎の振幅変
動量を演算する振幅変動量演算部と、 該振幅変動量と、該振幅変調波の特徴量に対して極値を
持つ聴感特性係数とに基づき評価値を演算する評価値演
算部と、 該評価値算出部の出力に基づき音色を評価する音色評価
部と、から成ることを特徴とする音色評価装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28938892A JP3146691B2 (ja) | 1992-10-02 | 1992-10-02 | 音色評価方法及びその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28938892A JP3146691B2 (ja) | 1992-10-02 | 1992-10-02 | 音色評価方法及びその装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06117912A JPH06117912A (ja) | 1994-04-28 |
| JP3146691B2 true JP3146691B2 (ja) | 2001-03-19 |
Family
ID=17742579
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28938892A Expired - Fee Related JP3146691B2 (ja) | 1992-10-02 | 1992-10-02 | 音色評価方法及びその装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3146691B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3388556B2 (ja) * | 1994-05-31 | 2003-03-24 | 株式会社アドバンテスト | 定常音波の評価方法及びその装置 |
| EP1772712A1 (en) | 2004-07-26 | 2007-04-11 | Sango Co., Ltd. | Method for analyzing sound quality of exhaust sound |
| JP2007271485A (ja) * | 2006-03-31 | 2007-10-18 | Ono Sokki Co Ltd | 音色評価支援方法および音色評価支援装置 |
| JP4902364B2 (ja) * | 2007-01-11 | 2012-03-21 | 帝人ファイバー株式会社 | 布帛摩擦時の摩擦音測定方法 |
| JP6073185B2 (ja) * | 2013-05-16 | 2017-02-01 | 株式会社小野測器 | 波形変換装置および波形変換方法 |
-
1992
- 1992-10-02 JP JP28938892A patent/JP3146691B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06117912A (ja) | 1994-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fastl | Fluctuation strength and temporal masking patterns of amplitude-modulated broadband noise | |
| Moore | Cochlear hearing loss: physiological, psychological and technical issues | |
| JP5191750B2 (ja) | 音響装置 | |
| CN101208991B (zh) | 具有加强的高频再现功能的助听器以及处理声频信号的方法 | |
| Carlyon et al. | Excitation produced by Schroeder-phase complexes: Evidence for fast-acting compression in the auditory system | |
| US6675114B2 (en) | Method for evaluating sound and system for carrying out the same | |
| CN105163262A (zh) | 一种扬声器异音检测方法及检测系统 | |
| US8147422B2 (en) | Method for detecting hearing impairment and device thereof | |
| JPH10267742A (ja) | 音質評価装置および音質評価方法 | |
| JP3146691B2 (ja) | 音色評価方法及びその装置 | |
| JP3922061B2 (ja) | 音質評価装置及び音質評価方法 | |
| Turner et al. | The relation between vowel recognition and measures of frequency resolution | |
| US5805457A (en) | System for analyzing sound quality in automobiles using musical intervals | |
| Wakita et al. | Objective rating of rumble in vehicle passenger compartment during acceleration | |
| CN117690452A (zh) | 电机信号处理方法、装置、设备及介质 | |
| US5535131A (en) | System for analyzing sound quality in automobile using musical intervals | |
| Genuit | The use of psychoacoustic parameters combined with A-weighted SPL in noise description | |
| US5944672A (en) | Digital hearing impairment simulation method and hearing aid evaluation method using the same | |
| Szabo et al. | Methods to study pre-school teachers' voice at work: simultaneous recordings with a voice accumulator and a DAT recorder | |
| JPH06105828A (ja) | 聴覚検査・訓練装置およびそれに用いるヘッドホン | |
| JPH11142231A (ja) | 騒音分析装置 | |
| Blommer et al. | A new method for calculating Fluctuation Strength in electric motors | |
| JP2002062186A (ja) | 地域環境騒音の計測・心理評価の方法及び装置 | |
| JP3420831B2 (ja) | 骨伝導音声のノイズ除去装置 | |
| JP2837640B2 (ja) | 補聴器 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080112 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090112 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |