JP3546236B2 - Noise psychological evaluation method, apparatus and medium - Google Patents
Noise psychological evaluation method, apparatus and medium Download PDFInfo
- Publication number
- JP3546236B2 JP3546236B2 JP2002205577A JP2002205577A JP3546236B2 JP 3546236 B2 JP3546236 B2 JP 3546236B2 JP 2002205577 A JP2002205577 A JP 2002205577A JP 2002205577 A JP2002205577 A JP 2002205577A JP 3546236 B2 JP3546236 B2 JP 3546236B2
- Authority
- JP
- Japan
- Prior art keywords
- acf
- iacf
- noise
- factor
- iacc
- 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 - Lifetime
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、航空機騒音や自動車騒音などの地域環境騒音の計測・心理評価の方法及び装置に関するものである。特にバイノーラル方式による騒音の計測・心理評価の方法及び装置に関するものである。
【0002】
【従来の技術】
従来、航空機騒音や自動車騒音などの地域環境騒音は、モノオーラル方式による騒音計を用いて測定した音圧レベルやその周波数特性に関して議論されてきた。しかし、上述したモノオーラル方式により測定された物理的ファクターのみでは人間の主観的応答を表わすには不十分かつ不適切であることがわかってきた。また、コンサートホール音響学では、バイノーラル方式により、ホールの物理的なデータと心理的(主観的)な関連性が明らかとなってきているが、騒音の分野においてはモノオーラル方式に関するものが殆どである。
【0003】
【発明が解決しようとする課題】
長年の間、環境騒音は、音圧レベル(SPL;Sound Pressure Level)の統計値を用いて評価されてきた。このSPLは、LxまたはLeqで表わされ、これのパワースペクトルは、モノオーラル騒音計で測定する。しかしながら、このSPL及びパワースペクトルだけでは環境騒音の主観的な評価には適さない。
【0004】
即ち、本発明の目的は、人間の聴覚−大脳機能システムにもとづき、時間領域において時々刻々変化する自己相関関数及び相互相関関数から導出される物理ファクターを用いて、騒音源の種類を特定する方法、装置及び媒体を提供することである。また本発明の他の目的は、人間の聴覚−大脳機能システムにもとづき、時間領域において時々刻々変化する自己相関関数及び相互相関関数から導出される物理ファクターを用いて、より的確にラウドネス、ピッチ、音色、心理的時間感覚をはじめ、主観的拡がり感、騒音場の見かけの音源の幅などの心理評価を行う方法、装置及び媒体を提供することである。
【0005】
【課題を解決するための手段】
上述した目的を達成するために、音声採取手段を用いて環境騒音の音響信号を採取・記録する音響信号記録ステップと、この記録された音響信号からフーリエ変換を用いて演算手段により自己相関関数(ACF)を算出するACF演算ステップと、この算出されたACFから演算手段により各ACFファクターを求めるACFファクター演算ステップと、この求めた各ACFファクターを用いて演算手段により騒音源の種類を判定する判定ステップと、を含むことを特徴とする騒音源の種類を特定する方法を提供する。
【0006】
また、好適には、上述した騒音源の種類を特定する方法において、前記ACFファクター演算ステップが、前記計算されたACFからACFファクターである遅れ時間が0で表わされるエネルギー(Φ(0))、有効継続遅延時間(τe)、ACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)を計算する演算ステップを含み、 前記騒音源の種類を判定する判定ステップが、これらの計算されたACFファクターである遅れ時間が0で表わされるエネルギー(Φ(0))、有効継続遅延時間(τe)、ACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)からその対数と、予め作成してある騒音源の各ACFファクター毎の対応するテンプレートの対数との差の絶対値である距離をそれぞれ求めるステップと、予めACFファクターの各々の算術平均の標準偏差であるS2を、ACFファクターの全カテゴリーに対する標準偏差の算術平均であるS1で除算し、この除算したものの平方根である重み係数を各ACFファクター毎に求めるステップと、求めたそれぞれの距離に、予め求めておいた対応する各ACFファクターの重み係数を乗算し、合計の距離を求める合計距離演算ステップと、この求めた合計距離と、格納されているテンプレートの距離とを比較し、最も近いテンプレートの1つを選択する比較・選択ステップと、を含むことを、特徴とする騒音源の種類を特定する方法を提供する。
【0007】
本発明の他の目的を達成するためには、音声採取手段を用いて環境騒音の音響信号をバイノーラル方式で記録する音響信号記録ステップと、このバイノーラル方式で記録された音響信号から演算手段を用いて自己相関関数(ACF)及び左右の各チャンネル間の相互相関関数(IACF)を計算するACF及びIACF演算ステップと、この計算されたACFから前記演算手段を用いて各ACFファクターを計算し、及び/またはこの計算されたIACFから各IACFファクターを計算するACF・IACFファクター演算ステップと、この計算されたACF及び/またはIACFファクターの各々に基づき演算手段を用いて心理評価を行う心理評価ステップと、を含むことを特徴とする騒音源について心理評価を行う方法を提供する。
【発明の実施の形態】
【0008】
ラウドネス、ピッチ、音色などの基本的な知覚データと同様に、嗜好や拡散性などの多くの主観的なデータの記述は、人間の聴覚−大脳システムの音場に対する応答モデルに基づいている。この応答モデルは予測されてきたが、それは経験的に得られた結果と一致することが知られている。例えば最近、周波数帯域幅を制限したノイズのラウドネスは、SPLによって影響をうけるのと同様に、自己相関関数(ACF)における有効継続時間(τe)によって影響を受けることが知られている。また、複合音の基本周波数が約1200Hzよりも低い場合、ピッチ及びその強さは、それぞれACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)によって影響を受ける。特に、ある時間内におおいて求められたτeの最小値(τe)minで得られるACFファクターは、騒音源及び騒音場の主観的評価の差異を良く表わすものである。
【0009】
このモデルは、2つのそれぞれの経路における音響信号同士の自己相関と、これらの音響信号の間における相互相関とから構成され、人間の大脳半球の処理特性も考慮するものである。即ち、両耳に入ってくる音響信号を用いて、自己相関関数(ACF)及び相互相関関数(IACF)を計算する。直交ファクターである遅れ時間が0で表わされるエネルギー(Φ(0))、有効継続遅延時間(τe)、ACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)はACFから導出される。また、IACFファクターである聴取音圧レベル(LL)、最大振幅(IACC)、最大振幅までの遅延時間(τIACC)、最大振幅における幅(WIACC)は、IACFから導出される。
【0010】
図1は、本発明による装置の具体的な構成を示す装置概略図である。図1に示すように本発明による装置の具体例は、聴者の頭部の模型1に装着された騒音源からの音響信号を採取するバイノーラル方式の音声採取手段2(マイクロフォン)を、LPF3(ローパスフィルタ)、A/Dコンバータ4、コンピュータ5から構成される。この頭部としては、人体の頭部が最も望ましいがそれでは不便であるため、人体の頭部を模したダミーヘッドを用いることもできる。しかし、このダミーヘッドは高価であり、ダミーヘッド以外の頭部の模型1(発砲スチロールなどの材料を用いた球体(直径を20cm)としたもの)でも本発明で測定するACF、IACFでは、有意差がないため、発砲スチロール製の頭部の模型を用いた。このコンピュータ5は、採取された音響信号を格納する音響信号記憶手段6と、この格納された音響信号(左右2チャンネル)を読み出し、これらの音響信号に基づきACFを計算するACF演算手段7、とこれらの音響信号に基づきにIACFを計算するIACF演算手段8、この計算されたACFに基づきACFファクターを計算するACFファクター演算手段9、この計算されたIACFに基づきIACFファクターを計算するIACFファクター演算手段10、この計算されたACFファクターに基づき騒音源の種類を特定する騒音源の種類を特定する手段11、この計算されたACFファクター及び/またはIACFファクターに基づき心理評価を行う手段12、騒音源の種類の特定及び心理評価に用いるデータに関するデータベース13を具える。
【0011】
聴者の頭部の模型1の両端に取り付けた左右2チャンネルのコンデンサマイクロフォン(マイクアンプ付き)を、ローパスフィルタを介して可搬型パーソナルコンピュータ5のサウンド入出力端子(A/D変換部4)と接続する。このマイクロフォン(音響信号採取手段2)から周りの騒音の取り込みを行う。コンピュータ上のプログラムの管理下、計測、各物理ファクタの算出、騒音源の種類の特定、心理評価、などを行う。また、騒音源の種類の特定及び心理評価に用いるデータに関するデータベースを構築する。
【0012】
図2は、本発明による騒音源の種類の特定、心理評価を行う方法のフローチャートである。図2に示すように、ステップS1では、騒音源からの音響信号を音源採取手段2により採取する。この採取された音響信号はLPF3を介してA/Dコンバータ4によりデジタル信号に変換する。ステップS2では、ステップS1で採取された音響信号を音響信号記憶手段に格納する。ステップS3では、ステップS2で格納された音響信号を読み出す。ステップS4では、ステップS3で読み出された音響信号に基づきACF及びIACFをACF演算手段7及びIACF演算手段8により計算する。ステップS5では、ステップS4で計算されたACF及びIACFに基づきACFファクター演算手段9及びIACF演算手段10によりACFファクター及びIACFファクターを計算する。ステップS6では、ステップS5で計算されたACFファクター及びIACFファクターに基づき、騒音源種類特定手段11、心理評価手段12により騒音源の種類の特定、心理評価を行う。その特定、評価の際には、テンプレートを格納するデータベース13からデータを読み出し比較・検討を行う。
【0013】
まず初めに、ピーク検知プロセスにより、採取した音響信号から複数の測定セッションを抽出する。連続的な騒音から自動的に環境騒音や目的の騒音を抽出するために、左右それぞれの耳の入り口部位におけるエネルギーであるモノオーラルのエネルギーΦll(0)、Φrr(0)を連続的に分析する。図3は、ピーク検知処理手順を説明する図であって、縦軸にノイズレベル、横軸に時間をとったグラフであって、その下段に積分間隔を示す図である。騒音が航空機騒音や列車騒音などの連続騒音の場合、Φ(0)の計算のための間隔を、かなり長く(例えば1秒など)設定することができるが、騒音が短時間や断続的である場合は、より短い間隔を用いる必要がある。しかしながら、後述する式(1)で連続計算する場合、積分間隔よりも長い間隔を選ぶ必要がある。従って、この間隔は、騒音源の種類に応じて決定する必要がある。
【0014】
これによって、長い時間の間隔で普通の騒音計を用いてΦ(0)を決定するより、より正確にΦ(0)を決定することができる。ピークを検出するためには、前もってトリガーレベルLtrigを適切に設定しておく必要がある。適当なLtrig値は、目標とする騒音の種類、目標とする騒音と観察者との距離、大気の条件などに応じて変化するものである。従って、この値を予備測定によって決定する必要がある。目的騒音と観察者との距離が近くて、かつ、観察者の近くに干渉する騒音源がない場合、Ltrig値を決定することは容易である。
【0015】
最大値Φ(0)を中心とする騒音を、システムを用いて単一のセッションで記録する。各々の目的とする騒音に対する1つのセッションの継続時間すなわちtsは、Ltrig値を超えた後にΦ(0)のピークを含むように選択する。航空機騒音や列車騒音などの普通の環境騒音の場合は、ts値は約10秒である。これは、継続時間が長い定常状態の騒音と短い継続時間の断続的な騒音とでは異なる。このシステムは、干渉する騒音がある場合には使えないことに留意されたい。図3に示すように、一連のセッション(S1(t),S2(t),S3(t),…SN(t)、N:セッションの数、0<t<ts)をシステム上に自動的に格納する。
【0016】
図3に示すように、継続時間tsでの各セッションSN(t)に対するランニングACF及びランニングIACFを分析する。ここでは、「ランニング」のプロセスを説明するために単一のセッションのみを考えることとする。計算の前に、適切な積分間隔2T及び連続ステップtstepの値を決定する。前述したように、推奨される積分間隔は約30×(τe)min[ms]であり、この(τe)minは一連の値τeの最小値であり、予備測定で容易に発見し得るものである。これは、違う種類の環境騒音のデータを用いて見つけるものである。大抵の場合、隣接する積分間隔をお互いに重ね合わせる。
【0017】
ACFとIACFを、2Tの範囲での1セッションごとの各ステップ(n=1,2,…,M)につき計算する。各ステップは、
【数1】
のようにtstepずつシフトする。物理ファクターは、ACF及びIACFの各ステップから導出する。2Tは予測されるτeの値よりも十分長くする必要がある。また、これは、各ステップに対する知覚の「聴覚の時間窓」に大きく関連する。環境騒音に対する2Tとしては、概ね0.1〜0.5秒が適している。2Tがこの範囲よりも小さい場合、(τe)minがある値に収束する。一般的に、tstepは0.1秒が好適である。変動が細かい場合は、より短いtstepを選択する。よく知られているように、バイノーラル信号をFFT(高速フーリエ変換)と、その後逆FFTの処理を行うことにより、ACF及びIACFを得ることができる。A特性フィルター及び、マイクロフォンの周波数特性は、FFT処理の後で考慮する。
【0018】
左右の耳の部位におけるACFを、それぞれ、Φll(τ)、Φrr(τ)で表わす。特定の数字の場合は、Φll (i)、Φrr (i)で表わす(1<i<Tf、 f:サンプリング周波数(Hz)、i:整数)。左右のΦ(0)を計算するためには、Φll (i)とΦrr (i)を下記のように平均する。
【数2】
SPLの正確な値は、次式で得られる。
【数3】
Φref(0)は、基準音圧値20μPにおけるΦ(0)である。
【0019】
バイノーラルの聴取音圧レベルは、Φll(0)及びΦrr(0)の相乗平均である。
【数4】
このΦ(0)は、IACFを正規化する際の分母となるものであるため、IACFファクターの一方のもの、或いは右半球の空間ファクターに分類されるものと考える。正規化したACFの振幅が0.1(10%の遅延)になる時の遅延時間によって、有効継続時間τeを定義する。正規化した左右の耳におけるACF、φll,rr(τ)は、次式で得られる。
【数5】
【0020】
図4は、縦軸にACFの対数の絶対値、横軸に遅延時間をとったグラフである。
図4に示すように初期のACFが線形に減少するのが一般的に観察できるため、縦軸をデシベル(対数)に変換するとτeを容易に得ることができる。線形回帰の場合は、ある一定の短い時間Δτにおいて得られるACFのピークに対して最小平均自乗法(LMS)を使用する。このΔτは、ACFのピークを検知するために使用され、計算前に慎重に決定しておく必要がある。τeを計算する際、原点が回帰線上にない場合、ACFの原点(ACF=0、τ=0)を、考慮に入れなくても良い場合も多い。極端な例では、目的とする騒音が純音とホワイトノイズとを含む場合、原点において急激な減衰が観察される。その後の減衰は、純音成分のため一定に保たれる。この場合、ACF関数の解は求まらない。
【0021】
図5は、縦軸に正規化したACF、横軸に遅延時間をとったグラフである。
図5に示すように、τ1は正規化したACFの第1のピークまでの遅延時間、φ1はその第1ピークでの振幅である。第1ピークは、局所的な小さなピークは無視して、主要なピークに基づき決定する。ファクターτnとφn(N≧2)とは考慮に入れない。なぜなら、τnとφnは、一般的にτ1とφ1とに相関関係があるからである。
【0022】
図6は、縦軸に正規化したIACF、横軸に左右の信号の遅延時間をとったグラフである。左右の耳の音響信号の間のIACFは、φlr(τ)(-1<τ<+1[ms])で表わされる。デジタル形式では、Φlr (i)(-f/103≦i≦f/103、iは整数であり、これが負の場合は左のチャンネルに遅れがあるIACFであることを示す)。両耳の間の最大遅延としては−1から+1msを考慮すれば十分である。最大振幅IACCは主観的拡散に関連するファクターである。図6に示すように、正規化されたIACFΦlr (i)の最大振幅は遅延範囲内で得られる。即ち
IACC={φlr (i)}max (5)
正規化されたIACFは次式で得られる。
【数6】
【0023】
τIACCの値は、最大振幅の遅延時間において容易に求まる。例えば、τIACCが正の場合、音源は聴者の右側に位置する、或いは音源が右側にあるかのように知覚する。図6に示すように、最大振幅における幅WIACCを、最大値から0.1(IACC)下の部分のピーク幅で得ることができる。この係数0.1はIACC=1.0におけるJNDとして概算的に用いられるものである。聴取音圧レベルLLは、式(2)でSPLをLLと置き換えることによって得られる。このようにして、各物理ファクターを、ACF及びIACFから求めることができる。
【0024】
次に、ACFファクターに基づき騒音源の種類の特定する方法について説明する。
騒音源の種類は、4つのACFファクター遅れ時間が0で表わされるエネルギー(Φ(0))、有効継続遅延時間(τe)、ACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)を用いて特定する。Φ(0)は騒音源と聴者との距離に応じて変化するため、距離が不明の場合は、計算の条件には特別に注意を払う必要がある。たとえファクターΦ(0)が有効でない場合であっても、その他の3つのファクターを用いて騒音源の種類を特定することができる。空間情報が変化する場合、残りのIACFファクターを考慮に入れることもできる。音響信号の最も大きく変動する部分である最小τe:(τe)minを用いる理由の1つは、この部分が主観的な応答に最も深く関与するものであるということである。
【0025】
未知の対象データ(下記の式(7)~(10)では記号aで示す)用の(τe)minにおける各ファクターの値とデータベースに格納されたテンプレート用(記号bで示す)の値との差、即ち「距離」を計算する。ここで「対象」とは、システムによって特定されるオブジェクトとしての環境騒音のことを意味する。テンプレート値は、ある特定の環境騒音に対する典型的なACFファクターのセットであり、これらの複数のテンプレートを未知の騒音と比較する。
距離D(x)(x:Φ(0)、τe、τ1、φ1)を次式により計算する。
【数7】
【0026】
目的とする騒音源の合計距離Dは、次式で表わされる。
【数8】
W(x)(x;Φ(0)、(τe)min、τ1、φ1)は、重み係数である。この算出された距離Dに最も近いDを有するテンプレートを、求める騒音源であると判断する。これにより、未知の騒音源が、何であるのか、例えば鉄道、自動車、航空機、工場騒音であるのか、更にその車種、機種などを特定することが可能となる。
【0027】
図7は重み係数の計算方法を説明するブロック図である。式(11)の重み係数W(x)(x;Φ(0)、τe、τ1、φ1)は、統計値S1 (i)とS2 (i)とを用いて得ることができる。図7に示すように、S1 (i)は、ACFファクターの全カテゴリーに対する標準偏差(SD)の算術平均である。ここでカテゴリーとは、同じ種類の騒音に対するデータのセットを意味する。S2 (i)は、各カテゴリの算術平均の標準偏差である。W(x)は、ファクター{(S2/S1)1/2}maxの中の最大値で正規化した後、(S2/S1)1/2で得られる。この平方根の処理は経験的に得られたものである。騒音源の間におけるより大きなSDと、ある騒音の間におけるより小さなSDとのファクターとは他の種類の騒音とは区別できるため、このようなファクターの重みはその他のファクターのものよりも大きくなる。テンプレートを改善する学習機能がある場合、システム上においてテンプレートは、システム内でACFの各ファクターについての最新の値と、元の値との平均によって上書きすることもできる。
【0028】
図8は、聴覚−大脳機能システムのモデルを説明するブロック図である。聴覚−大脳機能システムのモデルは、自己相関(ACF)メカニズム、両耳間相互相関(IACF)メカニズム、左右大脳の機能分化を含んでいる。信号のパワースペクトルに含まれる情報は、音響信号のACFにも含まれていることは注目すべきことである。また騒音場の空間的感覚を示すため、IACFより抽出される空間的ファクターを考慮する。音色は音の基本的感覚と空間的感覚を含む総合的な感覚として定義される
【0029】
聴覚−大脳機能モデル(図8)を使って、自由空間内に存在する聴者の正面にある与えられた音響信号p(t)の基本的な感覚を考える。ここで長時間ACFを次式で得ることができる。
【数9】
p’(t)=p(t)*s(t)で、s(t)は耳の感度である。便宜上s(t)はA特性のインパルス応答が用いられる。パワースペクトルも次式のようにACFから得ることができる。
【数10】
このように、ACFとパワースペクトルは数学的には同じ情報を含んでいる。
【0030】
ACFの解析において3つの重要な事項として、遅れ時間が0で表わされるエネルギーΦp(0)と、正規化したACFのエンベロープから抽出される有効継続時間τeと、ピークやディップやその遅れ時間とを含む微細構造とがある。図4に示すように、この有効継続時間τeは、10パーセント遅れ時間として定義でき、騒音響信号それ自身に含まれる繰り返し成分、または残響成分として表わされる。前述したように正規化したACFはΦp(τ)=Φp(τ)/Φp(0)で得ることができる。
【0031】
ラウドネスSLは次式で表わされる。
SL=fL(Φ(0),τ1,φ1,τe) (15)
即ち、ACFファクターである、遅れ時間が0で表わされるエネルギー(Φ(0))、有効継続遅延時間(τe)、ACFの第1ピークまでの遅延時間(τ1)、正規化したACFの第1ピークの振幅(φ1)からラウドネスを求めることができる。
ここでτ1は騒音のピッチまたは後述するミッシングファンダメンタル現象に関係するものである。また、p’(t)が音圧レベルL(t)を与えるための圧力20μPaを基準として測定されるなら、等価騒音レベルLeqは次式で求めることができる。
【数11】
このLeqは10logΦp(0)に相当するものである。また、サンプリング周波数は、最大可聴周波数域の2倍以上としなければならないので、通常の騒音計で測定されたLeqよりも極めて精度良く測定できる。
【0032】
図9は、縦軸にラウドネス尺度値、横軸にバンド幅をとったグラフである。このグラフは、Φp(0)を一定とした条件下での一対比較テスト(1080dB/octaveのスロープを持つフィルタを使用)で得られた臨界帯域内のラウドネス尺度値を示したものである。明らかに純音のような騒音が同じ繰り返し成分を持つとき、τeは大きな値となり、ラウドネスが大きくなる。このように、ラウドネス対バンド幅の関係は、臨界帯域内でも平坦にならないことがわかる。なお、この結果は中心周波数1kHzの周波数帯域で得られたものである。
【0033】
騒音のピッチまたはミッシングファンダメンタルは次式で表わされる。
Sp=fp(τ1,φ1) (17)
ここで、ミッシングファンダメンタル現象とは、いくつかの倍音構造が存在するとき、実際にはない高さの音が聞こえるという現象である。
【0034】
最も複雑な知覚である音色は、次式で表わされる。音色には、ラウドネスやピッチも含まれるものである。
ST=fT[Φ(0),τe,(τ1,φ1),…,(τn,φn)] (18)
τn,φn(n=1,2,…)の中でτ1,φ1が最も顕著な直交ファクターであるため式(18)は以下のように書き直すことができる。
ST=fT[Φ(0),τe,τ1,φ1] (19)
【0035】
信号の時間的長さの知覚に関する感覚は、次式で表わされる。
SD=fD[Φ(0),τe,τ1,φ1] (20)
【0036】
長時間IACFは次式で求めることができる。
【数12】
ここでp’l,r(t)=p(t)l,r*s(t)、はp(t)l,r は左右外耳道入り口の音圧である。
【0037】
騒音源の水平面の方向の知覚を含む空間情報の知覚は次式で表わされる。
S=f(LL,IACC,τIACC,WIACC) (22)
ここで聴取音圧レベルLLは{Φll(0),Φrr(0)}である。記号{}は、左右の耳の入り口に到来する信号のτ=0のときのACFであるΦll(0)、Φrr(0)の組を表わす。数学的にはLLは、両耳に到来する音響信号のエネルギーの算術平均で次式のように表わされる。
【数13】
式(22)で示す4つのIACFファクター(直交ファクター)の中で、−1〜+1msの範囲内のτIACCは、水平方向の音源の水平方向の知覚に関する重要なファクターである。正規化したIACFが1つの鋭いピークを持ち、IACCが大きく、高周波数成分によってWIACCが小さい値であるとき、明確な方向感が得られる。逆に主観的拡がり感やあいまいな方向感はIACCが小さい値(<0.15)の時に起こる。
【0038】
正中面に位置する騒音源の知覚については、耳の入り口に到来する音響信号の長時間ACFから抽出される時間的ファクターを式(22)に加えるべきであろう。
図8に示すように、注目すべきはIACCに相当する下丘付近に存在する神経活動の存在である。また、室内音場においては、LLとIACCとは右大脳半球に支配的に関連があり、時間的ファクターであるΔt1やTsubは左大脳半球と関わっていることを発見した。
【0039】
主観的拡がり感の尺度値を得るため、2つの対称な反射音の水平入射角度を変更し、ホワイトノイズを用いて一対比較テストを行った。被験者は、LL、τIACC、WIACCが一定の条件下で、提示された2つの音場のうち、どちらの音場がより広がって聞こえるかを判断した。図10は、左縦軸に拡がり感の尺度値、右縦軸に最大振幅IACC、横軸に反射音の水平入射角度をとったグラフである。図10に示すように、250Hz~4kHz(図10(a):250Hz、(b):500Hz、(c):1kHz、(d):2kHz、(e):4kHz)の周波数帯域の結果において、尺度値と最大振幅IACCとは強い負の相関関係を示した。従って、上述した実験結果により、主観的尺度値を、IACCの3/2乗で次式のように求めることができる。
Sdiffuseness=−α(IACC)β (24)
実験により求めた係数αは2.9、乗数βは3/2である。
【0040】
騒音場の見かけの音源の幅(ASW)を求める方法について説明する。低域の周波数成分が大きい騒音場では、長時間IACFは遅れ時間τが−1〜+1msの範囲内に明確なピークを持たず、WIACCは大きくなる。このWIACCは次式で求めることができる。
【数14】
ここで、Δωcは2π(f1+f2)、f1とf2とは、それぞれ理想的なバンドパスフィルターの下限値と上限値である。便宜上、δは0.1(IACC)と定義する。
【0041】
注目すべきことは、大きなASWは低周波数帯域が多く、IACCが小さいときに知覚されるということである。すなわち、LLが一定でτIACC=0の条件下では、ASWはIACCとWIACCのIACFファクターに基づき求めることができる。ASWの尺度値を10名の被験者を用いて一対比較テストで求めた。WIACCの値を制御するため、1/3オクターブバンドパスノイズの中心周波数を250Hz〜2kHzで変化させた。IACCは直接音に対する反射音のレベルの比を制御して調整した。聴取音圧レベルLLは、ASWに影響するので、全ての音場の耳の入り口でのトータル音圧レベルはピーク値が75dBAで一定とした。被験者は提示された2つの音場のうちどちらかが広がって聞こえるかを判断した。尺度値SASWの分散分析の結果、IACC、WIACCの両方のIACFファクター共に有意であり(p<0.01)、以下のようにSASWに対して独立に寄与している。従って、SASWを次式で求めることができる。
SASW=a(IACC)3/2+b(WIACC)1/2 (26)
ここで係数a=−1.64、b=2.44であり、これらの係数は、図11(a)(b)に示す10名の被験者の尺度値の回帰曲線から得られたものである。図11(a)は縦軸にASW、横軸にIACCをとったグラフであり、図11(b)は縦軸にASW、横軸にWIACCをとったグラフである。また、図11は、縦軸に実際に測定したASWの尺度値、横軸に計算されたASWの尺度値をとったグラフである。図12に示すように、この式から求めたSASWの尺度値と、SASWの測定値はよく対応することを確かめた(r=0.97、p<0.01)。
【0042】
時間的に変動する環境騒音を評価するため、短時間ランニングACF及び短時間ランニングIACFを用いる。前述と同様の方法で抽出された短時間ランニングの空間的・時間的ファクターは、時変動する騒音場の基本的感覚を示すのに用いられている。短時間ACFは次式で求めることができる。
【数15】
ここで2Tは解析される信号の長さである。この長さ2Tは、ランニングACFの有効継続時間の最小値(τe)minを少なくとも含む範囲で決定すべきである。(τe)minを示す騒音は信号が最も急速に変動することを表わしており、この部分が最も主観的応答に影響を及ぼしている。
【0043】
各騒音の部分におけるラウドネスSLに関して、式(15)は次式のように書き換えることができる。
SL=fL(LL,τ1,φ1,τe) (28)
ここで各ファクターは各騒音の部分について得られ、式(15)のΦ(0)はLLに置き換えられる。ACFから抽出された時間的ファクターが、室内の反射音群(Δt1,Δt2,...)と後続残響時間Tsubに影響を及ぼしているはずだということに注目すべきである。
【0044】
環境騒音場のピッチの記述で、騒音場において有意な時間的ファクターはτ1とφ1とであり、従って式(17)はそのまま保たれる。
【0045】
環境騒音場の音色は時間的・空間的ファクター全てで、次式のように表わされる。
ST=fT(τe,τ1,φ1;LL,IACC,τIACC,WIACC) (29)
ここで、人間の大脳半球が時間的ファクターが左大脳半球に関連し、空間的ファクターが右大脳半球に関連しているという専門化を考えると、式(29)は以下のように置き換えることができる。
ST=fT(τe,τ1,φ1)left+fT(LL、IACC,τIACC,WIACC)right(29)
弱い反射音の閾値をΔt1の関数として図13に示す。式(29)に含まれる、聴者に対する反射音の空間的方向(IACCとτIACC)と反射音の遅れ時間Δt1はこの閾値を示している。
【0046】
耳の感度は外耳と中耳とを含む物理システムにより特徴づけられる。音響信号を解析する前に、便宜上、A特性をかけておく。
単一反射音の遅れ時間を関数とした単音節の明瞭度は、母音と子音との間の部分の短時間ACFから抽出された4つの直交ファクターを解析することにより予測できる。最近の調査では、音色や比類似度の判断は、コンサートホール内の音場の主観的プリファレンスと同じく、総合的な主観的応答であることを明確に示している。音色と同様に、主観的プリファレンスは、τeの最小値を用いて表わされる。短時間積分時間は次式で表わされる。
(2T)=30(τe)min (30)
精神作業に関する騒音の影響は、作業能率と大脳の専門化との間の妨害現象として解釈することができる。ACFから抽出された時間的ファクターは、左大脳半球に関連しており、IACFから抽出されるファクターは右大脳半球に主に関わっている。
【図面の簡単な説明】
【図1】本発明による装置の具体的な構成を示す装置概略図である。
【図2】本発明による騒音源の種類の特定、心理評価を行う方法のフローチャートである。
【図3】図3は、ピーク検知処理手順を説明する図であって、縦軸にノイズレベル、横軸に時間をとったグラフであって、その下段に積分間隔を示す図である。
【図4】縦軸にACFの絶対値の対数、横軸に遅延時間をとったグラフである。
【図5】縦軸に正規化したACF、横軸に遅延時間をとったグラフである。
【図6】縦軸に正規化したIACF、横軸に左右の信号の遅延時間をとったグラフである。
【図7】重み係数の計算方法を説明するブロック図である。
【図8】聴覚−大脳機能システムのモデルを説明するブロック図である。
【図9】縦軸にラウドネス尺度値、横軸にバンド幅をとったグラフである。
【図10】左縦軸に拡がり感の尺度値、右縦軸に最大振幅IACC、横軸に反射音の水平入射角度をとったグラフである。
【図11】(a)は縦軸にASW、横軸にIACCをとったグラフであり、(b)は縦軸にASW、横軸にWIACCをとったグラフである。
【図12】縦軸に実際に測定したASWの尺度値、横軸に計算されたASWの尺度値をとったグラフである。
【図13】縦軸に信号の閾値、横軸に遅延時間をとったグラフである。
【符号の説明】
1 頭部の模型
2 バイノーラル方式の音声採取手段
3 LPF(ローパスフィルタ)、
4 A/Dコンバータ
5 コンピュータ
6 音響信号記憶手段
7 ACF演算手段
8 IACF演算手段
9 ACFファクター演算手段9
10 IACFファクター演算手段10
11 騒音源種類特定手段
12 心理評価手段
13 データベース[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a method and apparatus for measuring and psychological evaluation of local environmental noise such as aircraft noise and automobile noise. In particular, the present invention relates to a method and apparatus for measuring and psychological evaluation of noise by a binaural method.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, local environmental noises such as aircraft noise and automobile noise have been discussed with respect to sound pressure levels measured using a mono-aural sound level meter and their frequency characteristics. However, it has been found that the physical factors measured by the above-described mono-oral method alone are insufficient and inadequate to express human subjective responses. Also, in concert hall acoustics, the binaural method has revealed the psychological (subjective) relationship with the physical data of the hall, but in the noise field, most are related to the mono-aural method. is there.
[0003]
[Problems to be solved by the invention]
For many years, environmental noise has been evaluated using the statistics of Sound Pressure Level (SPL). This SPL is LxOr LeqAnd its power spectrum is measured with a mono-aural sound level meter. However, this SPL and power spectrum alone are not suitable for subjective evaluation of environmental noise.
[0004]
That is, an object of the present invention is to provide a method for specifying the type of a noise source using physical factors derived from an autocorrelation function and a cross-correlation function that change every moment in a time domain based on a human auditory-cerebral function system. , Devices and media. Another object of the present invention is to provide a more accurate loudness, a pitch, and a physical factor derived from an autocorrelation function and a cross-correlation function that change every moment in a time domain based on a human auditory-cerebral function system. An object of the present invention is to provide a method, an apparatus, and a medium for performing psychological evaluation such as a tone, a sense of psychological time, a subjective feeling of spreading, and an apparent sound source width of a noise field.
[0005]
[Means for Solving the Problems]
In order to achieve the above-mentioned object, an audio signal recording step of collecting and recording an acoustic signal of environmental noise using a voice sampling means, and an autocorrelation function (Fourier transform) from the recorded acoustic signal by a calculation means by an arithmetic means. ACF calculation step for calculating (ACF), ACF factor calculation step for calculating each ACF factor by the calculation means from the calculated ACF, and determination of the type of the noise source by the calculation means using the obtained ACF factors And a method for identifying the type of the noise source.
[0006]
Preferably, in the above-mentioned method for specifying the type of a noise source, the step of calculating the ACF factor includes an energy (Φ (0)) in which a delay time, which is an ACF factor, is represented by 0 from the calculated ACF. Effective duration delay time (τe), The delay time until the first peak of the ACF (τ1), The amplitude of the first peak of the normalized ACF (φ1), And the determining step of determining the type of the noise source includes: an energy (Φ (0)) having a delay time of 0, which is the calculated ACF factor; τe), The delay time until the first peak of the ACF (τ1), The amplitude of the first peak of the normalized ACF (φ1) To obtain distances which are absolute values of the difference between the logarithm of the logarithm and the logarithm of the corresponding template for each ACF factor of the noise source prepared in advance, and the standard deviation of the arithmetic mean of each of the ACF factors in advance Is S2Is the arithmetic mean of the standard deviations for all categories of the ACF factor, S1And a weighting factor that is a square root of the divided value is obtained for each ACF factor. Each obtained distance is multiplied by a weighting factor of the corresponding ACF factor obtained in advance to obtain a total distance. , And a comparing / selecting step of comparing the calculated total distance with the distance of the stored template and selecting one of the closest templates. A method is provided for identifying the type of noise source.
[0007]
In order to achieve another object of the present invention, an acoustic signal recording step of recording an acoustic signal of environmental noise in a binaural manner using a voice sampling unit, and using an arithmetic unit from the acoustic signal recorded in the binaural manner. ACF and IACF calculation steps for calculating an autocorrelation function (ACF) and a cross-correlation function (IACF) between the left and right channels, and calculating each ACF factor from the calculated ACF using the calculation means; And / or an ACF / IACF factor calculation step of calculating each IACF factor from the calculated IACF, and a psychological evaluation step of performing a psychological evaluation by using calculation means based on each of the calculated ACF and / or IACF factors, The present invention provides a method for performing a psychological evaluation on a noise source characterized by including:
BEST MODE FOR CARRYING OUT THE INVENTION
[0008]
As well as basic perceptual data such as loudness, pitch, and timbre, the description of many subjective data such as preferences and diffusivity is based on a human auditory-cerebral system response model to the sound field. Although this response model has been predicted, it is known to be consistent with empirical results. For example, recently, the loudness of frequency band limited noise is affected by the effective duration (τ) in the autocorrelation function (ACF), as well as by SPL.e) Are known to be affected. Further, when the fundamental frequency of the complex sound is lower than about 1200 Hz, the pitch and its intensity are respectively determined by the delay time (τ) until the first peak of the ACF.1), The amplitude of the first peak of the normalized ACF (φ1). In particular, the τ determined within a certain timee(Τe)minThe ACF factor obtained in (1) well represents the difference between the subjective evaluation of the noise source and the noise field.
[0009]
This model is composed of the autocorrelation between acoustic signals in two respective paths and the cross-correlation between these acoustic signals, and also takes into account the processing characteristics of the human cerebral hemisphere. That is, an auto-correlation function (ACF) and a cross-correlation function (IACF) are calculated using acoustic signals entering both ears. Energy (Φ (0)) in which the delay time, which is an orthogonal factor, is represented by 0, the effective continuation delay time (τe), The delay time until the first peak of the ACF (τ1), The amplitude of the first peak of the normalized ACF (φ1) Is derived from the ACF. Also, the listening sound pressure level (LL), the maximum amplitude (IAC), and the delay time (τ)IACC), Width at maximum amplitude (WIACC) Is derived from the IACF.
[0010]
FIG. 1 is an apparatus schematic diagram showing a specific configuration of the apparatus according to the present invention. As shown in FIG. 1, a specific example of the apparatus according to the present invention includes a binaural-type sound sampling means 2 (microphone) for sampling an acoustic signal from a noise source mounted on a
[0011]
The left and right two-channel condenser microphones (with microphone amplifier) attached to both ends of the listener's
[0012]
FIG. 2 is a flowchart of a method for specifying the type of noise source and performing psychological evaluation according to the present invention. As shown in FIG. 2, in step S1, an acoustic signal from a noise source is collected by the sound
[0013]
First, a plurality of measurement sessions are extracted from the collected acoustic signal by a peak detection process. In order to automatically extract environmental noise and target noise from continuous noise, mono-oral energy Φ, which is the energy at the entrance of the left and right earsll(0), Φrr(0) is analyzed continuously. FIG. 3 is a diagram for explaining the procedure of the peak detection processing, in which the vertical axis represents the noise level, the horizontal axis represents the time, and the lower part thereof represents the integration interval. When the noise is continuous noise such as aircraft noise or train noise, the interval for calculating Φ (0) can be set to be considerably long (for example, 1 second), but the noise is short or intermittent. If so, shorter intervals need to be used. However, in the case of performing the continuous calculation using the expression (1) described later, it is necessary to select an interval longer than the integration interval. Therefore, this interval needs to be determined according to the type of the noise source.
[0014]
As a result, Φ (0) can be determined more accurately than when Φ (0) is determined using an ordinary sound level meter at long time intervals. To detect a peak, the trigger level L must be set in advance.trigMust be set appropriately. Suitable LtrigThe value changes according to the type of target noise, the distance between the target noise and the observer, the atmospheric conditions, and the like. Therefore, this value needs to be determined by preliminary measurement. If the distance between the target noise and the observer is short and there is no interfering noise source near the observer, LtrigIt is easy to determine the value.
[0015]
The noise centered on the maximum value Φ (0) is recorded in a single session using the system. The duration of one session for each desired noise, ie tsIs LtrigSelect to include the peak of Φ (0) after exceeding the value. For normal environmental noise such as aircraft noise and train noise, tsThe value is about 10 seconds. This is different for long duration steady state noise and short duration intermittent noise. Note that this system cannot be used in the presence of interfering noise. As shown in FIG. 3, a series of sessions (S1(t), STwo(t), SThree(t),… SN(t), N: number of sessions, 0 <t <ts) Is automatically stored on the system.
[0016]
As shown in FIG.sEach session S atNAnalyze the running ACF and running IACF for (t). Here, only a single session will be considered to explain the "running" process. Before the calculation, a suitable integration interval 2T and a continuous step tstepDetermine the value of. As mentioned above, the recommended integration interval is about 30 × (τe)min[ms] and this (τe)minIs a series of values τe, Which can be easily found by preliminary measurement. This is found using data on different types of environmental noise. In most cases, adjacent integration intervals overlap each other.
[0017]
ACF and IACF are calculated for each step (n = 1, 2,..., M) for each session in the range of 2T. Each step is
(Equation 1)
Like tstepShift by one. Physical factors are derived from each step of ACF and IACF. 2T is the expected τeMust be much longer than the value of. This is also largely related to the perceptual "auditory time window" for each step. As the 2T for the environmental noise, 0.1 to 0.5 seconds is generally suitable. If 2T is smaller than this range, (τe)minConverges to a value. In general, tstepIs preferably 0.1 second. For small fluctuations, a shorter tstepSelect As is well known, ACF and IACF can be obtained by performing an FFT (Fast Fourier Transform) and then an inverse FFT on a binaural signal. The A characteristic filter and the frequency characteristics of the microphone are considered after the FFT processing.
[0018]
The ACF at the left and right ear parts isll(τ), ΦrrExpressed by (τ). For specific numbers, Φll (i), Φrr (i)(1 <i <Tf, f: sampling frequency (Hz), i: integer). To calculate left and right Φ (0), Φll (i)And Φrr (i)Are averaged as follows:
(Equation 2)
The exact value of SPL is given by:
(Equation 3)
Φref(0) is Φ (0) at the reference sound pressure value of 20 μP.
[0019]
The binaural listening sound pressure level is Φll(0) and ΦrrThe geometric mean of (0).
(Equation 4)
Since Φ (0) is a denominator when normalizing the IACF, it is considered that the Φ (0) is classified into one of the IACF factors or the space factor of the right hemisphere. The effective duration τ is determined by the delay time when the amplitude of the normalized ACF becomes 0.1 (delay of 10%).eIs defined. ACF, φ at the normalized left and right earsll, rr(Τ) is obtained by the following equation.
(Equation 5)
[0020]
FIG. 4 is a graph in which the vertical axis represents the absolute value of the logarithm of the ACF and the horizontal axis represents the delay time.
Since it is generally observed that the initial ACF decreases linearly as shown in FIG. 4, when the vertical axis is converted to decibel (logarithmic), τeCan be easily obtained. In the case of linear regression, the least mean square method (LMS) is used for the ACF peak obtained at a certain short time Δτ. This Δτ is used to detect the ACF peak, and needs to be carefully determined before calculation. τeIn the calculation of, if the origin is not on the regression line, the origin of the ACF (ACF = 0, τ = 0) often does not need to be taken into account. In an extreme example, when the target noise includes a pure sound and white noise, a sharp attenuation is observed at the origin. Subsequent decay is kept constant for pure tone components. In this case, no solution of the ACF function is obtained.
[0021]
FIG. 5 is a graph with the normalized ACF on the vertical axis and the delay time on the horizontal axis.
As shown in FIG.1Is the delay time to the first peak of the normalized ACF, φ1Is the amplitude at the first peak. The first peak is determined based on the main peak, ignoring local small peaks. Factor τnAnd φn(N ≧ 2) is not taken into account. Because τnAnd φnIs generally τ1And φ1This is because there is a correlation with
[0022]
FIG. 6 is a graph in which the ordinate represents the normalized IACF and the abscissa represents the delay time of the left and right signals. The IACF between the acoustic signals of the left and right ears is φlr(τ) (-1 <τ <+1 [ms]). In digital form, Φlr (i)(-F / 10Three≤i≤f / 10Three, I is an integer, and a negative value indicates an IACF with a delay in the left channel). It is sufficient to consider the maximum delay between both ears from -1 to +1 ms. The maximum amplitude IACC is a factor related to subjective spreading. As shown in FIG. 6, the normalized IACFΦlr (i)Are obtained within the delay range. That is
IACC = {φlr (i)}max (5)
The normalized IACF is obtained by the following equation.
(Equation 6)
[0023]
τIACCIs easily obtained at the maximum amplitude delay time. For example, τIACCIs positive, the sound source is located on the right side of the listener, or the sound source is perceived as if it were on the right side. As shown in FIG. 6, the width W at the maximum amplitudeIACCCan be obtained with a peak width of 0.1 (IACC) below the maximum value. This coefficient of 0.1 is roughly used as JND at IACC = 1.0. The listening sound pressure level LL is obtained by replacing SPL with LL in equation (2). In this way, each physical factor can be obtained from the ACF and the IACF.
[0024]
Next, a method of specifying the type of the noise source based on the ACF factor will be described.
The type of noise source is energy (Φ (0)) in which four ACF factor delay times are represented by 0, and effective continuation delay time (τ).e), The delay time until the first peak of the ACF (τ1), The amplitude of the first peak of the normalized ACF (φ1). Since Φ (0) changes according to the distance between the noise source and the listener, special attention must be paid to the calculation conditions when the distance is unknown. Even if the factor Φ (0) is not valid, the type of the noise source can be specified using the other three factors. If the spatial information changes, the remaining IACF factors can also be taken into account. The minimum τ which is the largest part of the acoustic signale: (Τe)minOne reason for using is that this part is most closely involved in the subjective response.
[0025]
(Τ) for unknown target data (indicated by symbol a in equations (7) to (10) below)e)minAnd the difference between the value of each factor and the value for the template (indicated by the symbol b) stored in the database, that is, the “distance” is calculated. Here, “target” means environmental noise as an object specified by the system. The template value is a set of typical ACF factors for a particular environmental noise, comparing these multiple templates with the unknown noise.
Distance D (x) (x: Φ (0), τe, Τ1, Φ1) Is calculated by the following equation.
(Equation 7)
[0026]
The total distance D of the target noise source is expressed by the following equation.
(Equation 8)
W(x)(X; Φ (0), (τe)min, Τ1, Φ1) Is a weight coefficient. The template having D closest to the calculated distance D is determined to be the noise source to be obtained. As a result, it is possible to specify what the unknown noise source is, for example, a railway, automobile, aircraft, or factory noise, as well as its model and model.
[0027]
FIG. 7 is a block diagram illustrating a method of calculating a weight coefficient. Weighting coefficient W in equation (11)(x)(X; Φ (0), τe, Τ1, Φ1) Is the statistic S1 (i)And STwo (i)And can be obtained by using As shown in FIG.1 (i)Is the arithmetic mean of the standard deviations (SD) for all categories of ACF factor. Here, the category means a set of data for the same type of noise. STwo (i)Is the standard deviation of the arithmetic mean of each category. W(x)Is the factor {(STwo/ S1)1/2}maxAfter normalizing with the maximum value in, (STwo/ S1)1/2Is obtained. The processing of this square root has been obtained empirically. The weight of such factors is greater than that of other factors, because the factor of a larger SD between noise sources and a smaller SD between certain noises can be distinguished from other types of noise. . If there is a learning function to improve the template, the template on the system can also be overwritten by the average of the latest value for each factor of the ACF and the original value in the system.
[0028]
FIG. 8 is a block diagram illustrating a model of the auditory-cerebral function system. Models of the auditory-cerebral function system include an autocorrelation (ACF) mechanism, a binaural cross-correlation (IACF) mechanism, and functional differentiation of the left and right cerebrum. It should be noted that the information contained in the power spectrum of the signal is also contained in the ACF of the audio signal. In order to indicate the spatial sensation of the noise field, a spatial factor extracted from the IACF is considered. Tone is defined as an overall sensation that includes the fundamental and spatial senses of sound
[0029]
Using the auditory-cerebral function model (FIG. 8), consider the basic sensation of a given acoustic signal p (t) in front of a listener in free space. Here, the long-time ACF can be obtained by the following equation.
(Equation 9)
p '(t) = p (t) * s (t), where s (t) is ear sensitivity. For convenience, s (t) uses an impulse response of the A characteristic. The power spectrum can also be obtained from the ACF as in the following equation.
(Equation 10)
Thus, the ACF and the power spectrum mathematically contain the same information.
[0030]
In the analysis of the ACF, three important matters are the energy Φ with a delay time of 0.p(0) and the effective duration τ extracted from the normalized ACF envelopeeAnd a fine structure including a peak, a dip, and a delay time thereof. As shown in FIG. 4, this effective duration τeCan be defined as a 10% delay time, and is represented as a repetitive component or a reverberant component included in the noise signal itself. The ACF normalized as described above is Φp(τ) = Φp(τ) / Φp(0).
[0031]
Loudness SLIs represented by the following equation.
SL= fL(Φ (0), τ1, φ1, τe) (15)
That is, the energy (Φ (0)), which is the delay time represented by 0, which is the ACF factor, and the effective continuation delay time (τe), The delay time until the first peak of the ACF (τ1), The amplitude of the first peak of the normalized ACF (φ1) Can be used to determine loudness.
Where τ1Is related to the noise pitch or the missing fundamental phenomenon described later. Also, if p ′ (t) is measured based on a pressure of 20 μPa for providing the sound pressure level L (t), the equivalent noise level LeqCan be obtained by the following equation.
(Equation 11)
This LeqIs 10logΦpIt is equivalent to (0). Since the sampling frequency must be at least twice the maximum audible frequency range, LeqMeasurement can be performed with much higher accuracy than the above.
[0032]
FIG. 9 is a graph with the loudness scale value on the vertical axis and the bandwidth on the horizontal axis. This graph is ΦpFIG. 9 shows loudness scale values within a critical band obtained by a pairwise comparison test (using a filter having a slope of 1080 dB / octave) under the condition that (0) is constant. Obviously, when noise like pure tone has the same repetitive component, τeHas a large value, and the loudness increases. Thus, it can be seen that the relationship between loudness and bandwidth is not flat even within the critical band. Note that this result was obtained in a frequency band with a center frequency of 1 kHz.
[0033]
The pitch or missing fundamental of noise is expressed by the following equation.
Sp= fp(τ1, φ1) (17)
Here, the missing fundamental phenomenon is a phenomenon in which, when some overtone structures are present, a sound having a height which is not actually present is heard.
[0034]
The tone that is the most complex perception is expressed by the following equation. Tones include loudness and pitch.
ST= fT[Φ (0), τe, (τ1, φ1),…, (Τn, φn)] (18)
τn, φnΤ in (n = 1,2,…)1, φ1Is the most prominent orthogonal factor, equation (18) can be rewritten as:
ST= fT[Φ (0), τe, τ1, φ1] (19)
[0035]
The sensation related to the perception of the temporal length of the signal is expressed by the following equation.
SD= fD[Φ (0), τe, τ1, φ1] (20)
[0036]
The long-time IACF can be obtained by the following equation.
(Equation 12)
Where p 'l, r(t) = p (t)l, r* s (t) is p (t)l, r Is the sound pressure at the entrance of the left and right ear canal.
[0037]
The perception of spatial information including the perception of the direction of the horizontal plane of the noise source is represented by the following equation.
S = f (LL, IACC, τIACC, WIACC) (22)
Here, the listening sound pressure level LL is {Φll(0), Φrr(0)}. The symbol {} is the ACF of the signal arriving at the entrance of the left and right ears when τ = 0, Φll(0), ΦrrRepresents the set (0). Mathematically, LL is the arithmetic mean of the energy of the acoustic signal arriving at both ears and is expressed as:
(Equation 13)
Of the four IACF factors (orthogonal factors) shown in equation (22), τ within the range of −1 to +1 msIACCIs an important factor for the horizontal perception of a horizontal sound source. The normalized IACF has one sharp peak, the IACC is large, and WIACCIs small, a clear sense of direction is obtained. Conversely, a sense of subjective spread and an ambiguous sense of direction occur when IACC is a small value (<0.15).
[0038]
For the perception of a noise source located in the median plane, the time factor extracted from the long-term ACF of the acoustic signal arriving at the ear entrance should be added to equation (22).
As shown in FIG. 8, what should be noted is the presence of neural activity near the inferior colliculus corresponding to IACC. In a room sound field, LL and IACC are predominantly related to the right cerebral hemisphere, and a time factor Δt1And TsubFound that he was involved in the left cerebral hemisphere.
[0039]
In order to obtain a scale value of the subjective spread, a pairwise comparison test was performed using white noise while changing the horizontal incident angle of two symmetric reflected sounds. Subjects are LL, τIACC, WIACCUnder certain conditions, it was determined which of the two presented sound fields could be heard more widely. FIG. 10 is a graph in which the left vertical axis represents the scale value of the feeling of spreading, the right vertical axis represents the maximum amplitude IACC, and the horizontal axis represents the horizontal incident angle of the reflected sound. As shown in FIG. 10, in the result of the frequency band of 250 Hz to 4 kHz (FIG. 10 (a): 250 Hz, (b): 500 Hz, (c): 1 kHz, (d): 2 kHz, (e): 4 kHz) The scale value and the maximum amplitude IACC showed a strong negative correlation. Therefore, based on the above-described experimental results, the subjective scale value can be obtained by the following equation using the 3/2 power of IACC.
Sdiffuseness= −α (IACC)β (24)
The coefficient α obtained by the experiment is 2.9, and the multiplier β is 3/2.
[0040]
A method of obtaining the apparent sound source width (ASW) of the noise field will be described. In a noise field having a large low-frequency component, the long-time IACF does not have a clear peak in the range of the delay time τ in the range of −1 to +1 ms.IACCBecomes larger. This WIACCCan be obtained by the following equation.
[Equation 14]
Where ΔωcIs 2π (f1+ fTwo), F1And fTwoAre the lower limit and upper limit of the ideal bandpass filter, respectively. For convenience, δ is defined as 0.1 (IACC).
[0041]
It should be noted that a large ASW is perceived when the low frequency band is high and the IACC is low. That is, if LL is constant and τIACC= 0, ASW is equal to IACC and WIACCCan be determined based on the IACF factor. ASW scale values were determined in a pairwise comparison test using 10 subjects. WIACCWas controlled, the center frequency of the 1/3 octave bandpass noise was changed from 250 Hz to 2 kHz. The IACC controlled and adjusted the ratio of the reflected sound level to the direct sound. Since the listening sound pressure level LL affects the ASW, the total sound pressure level at the entrances of the ears in all sound fields was set to a constant peak value of 75 dBA. The subject determined which of the two sound fields presented was widespread and audible. Scale value SASWAnalysis of variance, IACC, WIACCAre significant (p <0.01) for both IACF factors, and SASWIndependent contribution to Therefore, SASWCan be obtained by the following equation.
SASW= A (IACC)3/2+ B (WIACC)1/2 (26)
Here, the coefficients a = -1.64 and b = 2.44, and these coefficients are obtained from the regression curves of the scale values of the ten subjects shown in FIGS. 11 (a) and 11 (b). . FIG. 11A is a graph with ASW on the vertical axis and IACC on the horizontal axis, and FIG. 11B is ASW on the vertical axis and W on the horizontal axis.IACCIt is a graph that takes FIG. 11 is a graph in which the vertical axis indicates the scale value of the actually measured ASW, and the horizontal axis indicates the calculated scale value of the ASW. As shown in FIG. 12, SASWScale value and SASWIt was confirmed that the measured values corresponded well (r = 0.97, p <0.01).
[0042]
In order to evaluate temporally fluctuating environmental noise, a short-time running ACF and a short-time running IACF are used. The spatial and temporal factors of short running extracted in the same manner as described above are used to indicate the basic sensation of a time-varying noise field. The short-time ACF can be obtained by the following equation.
(Equation 15)
Here, 2T is the length of the signal to be analyzed. This length 2T is the minimum value (τ) of the effective duration of the running ACF.e)minShould be determined in a range including at least. (Τe)minIndicates that the signal fluctuates most rapidly, and this part has the most influence on the subjective response.
[0043]
Loudness S at each noise partLEquation (15) can be rewritten as:
SL= FL(LL, τ1, φ1, τe) (28)
Here, each factor is obtained for each noise portion, and Φ (0) in Expression (15) is replaced by LL. The temporal factor extracted from the ACF is the reflected sound group (Δt1, Δt2,. . . ) And subsequent reverberation time TsubIt should be noted that
[0044]
In the description of the pitch of the environmental noise field, the significant temporal factor in the noise field is τ1And φ1Therefore, equation (17) is kept as it is.
[0045]
The timbre of the environmental noise field is expressed by the following equation with all temporal and spatial factors.
ST= FT(Τe, τ1, φ1LL, IACC, τIACC, WIACC) (29)
Here, considering the specialization that the human cerebral hemisphere has a temporal factor related to the left cerebral hemisphere and a spatial factor related to the right cerebral hemisphere, equation (29) can be replaced as follows: it can.
ST= FT(Τe, τ1, φ1)left+ FT(LL, IACC, τIACC, WIACC)right(29)
Δt is the threshold value for weak reflected sound1FIG. The spatial directions (IACC and τ) of the reflected sound for the listener included in equation (29)IACC) And the delay time Δt1 of the reflected sound indicate this threshold value.
[0046]
Ear sensitivity is characterized by a physical system that includes the outer and middle ears. Before analyzing the acoustic signal, the A characteristic is applied for convenience.
The clarity of a single syllable as a function of the delay time of a single reflected sound can be predicted by analyzing four orthogonal factors extracted from the short-time ACF of the portion between the vowel and the consonant. Recent research has clearly shown that the judgment of timbre and relative similarity is an overall subjective response, as is the subjective preference of the sound field in a concert hall. Like the timbre, the subjective preference is τeUsing the minimum value of The short integration time is represented by the following equation.
(2T) = 30 (τe)min (30)
The effect of noise on mental work can be interpreted as a disturbing phenomenon between work efficiency and cerebral specialization. Temporal factors extracted from the ACF are related to the left cerebral hemisphere, and factors extracted from the IACF are mainly related to the right cerebral hemisphere.
[Brief description of the drawings]
FIG. 1 is an apparatus schematic diagram showing a specific configuration of an apparatus according to the present invention.
FIG. 2 is a flowchart of a method for specifying a type of noise source and performing psychological evaluation according to the present invention.
FIG. 3 is a diagram illustrating a procedure of a peak detection process, in which a vertical axis represents a noise level, a horizontal axis represents time, and a lower part thereof represents an integration interval.
FIG. 4 is a graph in which the vertical axis represents the logarithm of the absolute value of the ACF and the horizontal axis represents the delay time.
FIG. 5 is a graph with normalized ACF on the vertical axis and delay time on the horizontal axis.
FIG. 6 is a graph in which the vertical axis represents the normalized IACF and the horizontal axis represents the delay time of the left and right signals.
FIG. 7 is a block diagram illustrating a method of calculating a weight coefficient.
FIG. 8 is a block diagram illustrating a model of the auditory-cerebral function system.
FIG. 9 is a graph with the loudness scale value on the vertical axis and the bandwidth on the horizontal axis.
FIG. 10 is a graph in which a left vertical axis represents a scale value of a feeling of spreading, a right vertical axis represents a maximum amplitude IACC, and a horizontal axis represents a horizontal incident angle of a reflected sound.
11A is a graph with ASW on the vertical axis and IACC on the horizontal axis, and FIG. 11B is a graph with ASW on the vertical axis and W on the horizontal axis.IACCIt is a graph that takes
FIG. 12 is a graph in which the vertical axis represents the actually measured ASW scale value, and the horizontal axis represents the calculated ASW scale value.
FIG. 13 is a graph in which a vertical axis indicates a signal threshold and a horizontal axis indicates a delay time.
[Explanation of symbols]
1 Model of the head
2 Binaural sound sampling means
3 LPF (low pass filter),
4 A / D converter
5 Computer
6 Acoustic signal storage means
7 ACF calculation means
8 IACF calculation means
9 ACF factor calculation means 9
10 IACF factor calculation means 10
11 Noise source type identification means
12 psychological evaluation means
13 Database
Claims (9)
このバイノーラル方式で記録された音響信号から演算手段を用いて自己相関関数ACF及び左右の各チャンネル間の相互相関関数IACFを計算するACF及びIACF演算ステップと、
前記演算手段を用いて、この計算されたACFからACFファクターを計算し、この計算されたIACFからIACFファクターを計算するACF・IACFファクター演算ステップと、
前記演算手段を用いて、この計算されたACFファクター、及び、IACFファクターに基づき心理評価を行う心理評価ステップと、
を含むことを特徴とする騒音心理評価方法。An acoustic signal recording step of recording the acoustic signal of the environmental noise in a binaural manner using a voice sampling means,
ACF and IACF calculation steps of calculating an autocorrelation function ACF and a cross-correlation function IACF between each of the left and right channels from the acoustic signal recorded in the binaural method by using calculation means;
An ACF / IACF factor calculating step of calculating an ACF factor from the calculated ACF using the calculating means and calculating an IACF factor from the calculated IACF;
A psychological evaluation step of performing a psychological evaluation based on the calculated ACF factor and the IACF factor using the arithmetic means;
A noise psychological evaluation method comprising:
前記心理評価を行うステップが、
前記計算されたIACFファクターである最大振幅IACCと係数αとに基づき、主観的拡がり感Sdiffusenessを、
Sdiffuseness=−α(IACC)3/2
により求める演算ステップ、
を含むことを特徴とする騒音心理評価方法。In the noise psychological evaluation method according to claim 1,
Performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC and the coefficient α, the subjective spread feeling S diffuseness is calculated as follows :
S diffuseness = -α (IACC) 3/2
Calculation step obtained by
A noise psychological evaluation method comprising:
前記心理評価を行うステップが、
前記計算されたIACFファクターである最大振幅IACCと、最大振幅における幅WIACCと、係数a及びbとに基づき、見かけの音源の幅SASWを、
SASW=−a(IACC)β+b(WIACC)1 /2
により求める演算ステップ、
を含むことを特徴とする騒音心理評価方法。In the noise psychological evaluation method according to claim 1,
Performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC, the width W IACC at the maximum amplitude, and the coefficients a and b, the apparent sound source width S ASW is calculated as follows:
S ASW = -a (IACC) β + b (W IACC) 1/2
Calculation step obtained by
A noise psychological evaluation method comprising:
このバイノーラル方式で記録された音響信号から演算手段を用いて自己相関関数ACF及び左右の各チャンネル間の相互相関関数IACFを計算するACF及びIACF演算手段と、
前記演算手段を用いて、この計算されたACFからACFファクターを計算し、この計算されたIACFからIACFファクターを計算するACF・IACFファクター演算手段と、
前記演算手段を用いて、この計算されたACFファクター、及び、IACFファクターに基づき心理評価を行う心理評価手段と、
を含むことを特徴とする騒音心理評価装置。Acoustic signal recording means for recording an acoustic signal of environmental noise in a binaural manner using a voice sampling means,
ACF and IACF calculating means for calculating an autocorrelation function ACF and a cross-correlation function IACF between each of the left and right channels from the acoustic signal recorded in the binaural method by using calculating means;
ACF / IACF factor calculating means for calculating an ACF factor from the calculated ACF using the calculating means and calculating an IACF factor from the calculated IACF;
A psychological evaluation means for performing a psychological evaluation based on the calculated ACF factor and the IACF factor using the arithmetic means;
A noise psychological evaluation device characterized by including:
前記心理評価を行う手段が、
前記計算されたIACFファクターである最大振幅IACCと係数αとに基づき、主観的拡がり感Sdiffusenessを、
Sdiffuseness=−α(IACC)3/2
により求める演算手段、
を含むことを特徴とする騒音心理評価装置。The noise psychological evaluation device according to claim 4,
The means for performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC and the coefficient α, the subjective spread feeling S diffuseness is calculated as follows :
S diffuseness = -α (IACC) 3/2
Calculation means obtained by
A noise psychological evaluation device characterized by including:
前記心理評価を行う手段が、
前記計算されたIACFファクターである最大振幅IACCと、最大振幅における幅WIACCと、係数a及びbとに基づき、見かけの音源の幅SASWを、
SASW=−a(IACC)β+b(WIACC)1 /2
により求める演算手段、
を含むことを特徴とする騒音心理評価装置。The noise psychological evaluation device according to claim 4,
The means for performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC, the width W IACC at the maximum amplitude, and the coefficients a and b, the apparent sound source width S ASW is calculated as follows:
S ASW = -a (IACC) β + b (W IACC) 1/2
Calculation means obtained by
A noise psychological evaluation device characterized by including:
このバイノーラル方式で記録された音響信号から演算手段を用いて自己相関関数ACF及び左右の各チャンネル間の相互相関関数IACFを計算するACF及びIACF演算ステップと、
前記演算手段を用いてこの計算されたACFからACFファクターを計算し、この計算されたIACFからIACFファクターを計算するACF・IACFファクター演算ステップと、
前記演算手段を用いてこの計算されたACFファクター、及び、IACFファクターに基づき心理評価を行う心理評価ステップと、を含むことを特徴とする、
騒音心理評価を行うプログラムを記録したコンピュータ読み取り可能な媒体。An acoustic signal recording step of recording the acoustic signal of the environmental noise in a binaural manner using a voice sampling means,
ACF and IACF calculation steps of calculating an autocorrelation function ACF and a cross-correlation function IACF between each of the left and right channels from the acoustic signal recorded in the binaural method by using calculation means;
An ACF / IACF factor calculation step of calculating an ACF factor from the calculated ACF using the calculation means and calculating an IACF factor from the calculated IACF;
A psychological evaluation step of performing a psychological evaluation based on the ACF factor calculated using the arithmetic means and the IACF factor.
A computer-readable medium on which a program for evaluating noise psychology is recorded.
前記心理評価を行うステップが、
前記計算されたIACFファクターである最大振幅IACCと係数αとに基づき、主観的拡がり感Sdiffusenessを、
Sdiffuseness=−α(IACC)3/2
により求める演算ステップ、を含むことを特徴とする、
騒音心理評価を行うプログラムを記録したコンピュータ読み取り可能な媒体。A computer-readable medium storing the program according to claim 7,
Performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC and the coefficient α, the subjective spread feeling S diffuseness is calculated as follows :
S diffuseness = -α (IACC) 3/2
A calculation step of:
A computer-readable medium on which a program for evaluating noise psychology is recorded.
前記心理評価を行うステップが、
前記計算されたIACFファクターである最大振幅IACCと、最大振幅における幅WIACCと、係数a及びbとに基づき、見かけの音源の幅SASWを、
SASW=−a(IACC)β+b(WIACC)1 /2
により求める演算ステップ、を含むことを特徴とする、
騒音心理評価を行うプログラムを記録したコンピュータ読み取り可能な媒体。A computer-readable medium storing the program according to claim 7,
Performing the psychological evaluation,
Based on the calculated IACF factor, the maximum amplitude IACC, the width W IACC at the maximum amplitude, and the coefficients a and b, the apparent sound source width S ASW is calculated as follows:
S ASW = -a (IACC) β + b (W IACC) 1/2
A calculation step of:
A computer-readable medium on which a program for evaluating noise psychology is recorded.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002205577A JP3546236B2 (en) | 2002-07-15 | 2002-07-15 | Noise psychological evaluation method, apparatus and medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002205577A JP3546236B2 (en) | 2002-07-15 | 2002-07-15 | Noise psychological evaluation method, apparatus and medium |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000246418A Division JP3350713B2 (en) | 2000-08-15 | 2000-08-15 | Method, apparatus and medium for identifying type of noise source |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003121253A JP2003121253A (en) | 2003-04-23 |
JP3546236B2 true JP3546236B2 (en) | 2004-07-21 |
Family
ID=19195765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002205577A Expired - Lifetime JP3546236B2 (en) | 2002-07-15 | 2002-07-15 | Noise psychological evaluation method, apparatus and medium |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3546236B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2021025800A (en) * | 2019-07-31 | 2021-02-22 | 日東電工株式会社 | Evaluation system, evaluation device, method for evaluation, program, and noise insulation device |
-
2002
- 2002-07-15 JP JP2002205577A patent/JP3546236B2/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2003121253A (en) | 2003-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6675114B2 (en) | Method for evaluating sound and system for carrying out the same | |
Ando | A theory of primary sensations and spatial sensations measuring environmental noise | |
Plomp | The role of modulation in hearing | |
Šrámková et al. | The softest sound levels of the human voice in normal subjects | |
EP2881948A1 (en) | Spectral comb voice activity detection | |
Fujii et al. | Temporal and spatial factors of traffic noise and its annoyance | |
Soeta et al. | Three dimensional psychological evaluation of aircraft noise and prediction by physical parameters | |
JP3350713B2 (en) | Method, apparatus and medium for identifying type of noise source | |
US20040179696A1 (en) | Sound field control system and sound field controlling method, as well as sound field space characteristic decision system and sound field space characteristic deciding method | |
JP2007534990A (en) | Judgment method of subjective characteristics of binaural sound signal | |
JP3584287B2 (en) | Sound evaluation method and system | |
JP4590545B2 (en) | Acoustic evaluation method and system | |
JP3546236B2 (en) | Noise psychological evaluation method, apparatus and medium | |
Rubak | Coloration in room impulse responses | |
Lavandier et al. | Identification of some perceptual dimensions underlying loudspeaker dissimilarities | |
Reuter et al. | Distortion product otoacoustic emission of symphony orchestra musicians before and after rehearsal | |
Sakurai et al. | A computational software for noise measurement and toward its identification | |
US20050004792A1 (en) | Speech characteristic extraction method speech charateristic extraction device speech recognition method and speech recognition device | |
CN106710602A (en) | Acoustic reverberation time estimation method and device | |
Sottek et al. | Sound quality evaluation of fan noise based on hearing-related parameters | |
Hotehama et al. | Dissimilarity judgments in relation to temporal and spatial factors for the sound field in an existing hall | |
US11832936B2 (en) | Methods and systems for evaluating hearing using cross frequency simultaneous masking | |
Sottek | Sound quality evaluation of noises with spectro-temporal patterns | |
Mendonça | Noise perception, psychoacoustic indicators and traffic noise | |
Baasch et al. | Parkinson-Speech Analysis: Methods and Aims |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040113 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20040316 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3546236 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |