JP4130779B2 - 音場制御システム及び音場制御方法 - Google Patents
音場制御システム及び音場制御方法 Download PDFInfo
- Publication number
- JP4130779B2 JP4130779B2 JP2003067814A JP2003067814A JP4130779B2 JP 4130779 B2 JP4130779 B2 JP 4130779B2 JP 2003067814 A JP2003067814 A JP 2003067814A JP 2003067814 A JP2003067814 A JP 2003067814A JP 4130779 B2 JP4130779 B2 JP 4130779B2
- Authority
- JP
- Japan
- Prior art keywords
- characteristic
- sound field
- sound
- dividing
- space
- 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
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
- H04S7/30—Control circuits for electronic adaptation of the sound field
- H04S7/307—Frequency adjustment, e.g. tone control
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
- H04S7/30—Control circuits for electronic adaptation of the sound field
- H04S7/305—Electronic adaptation of stereophonic audio signals to reverberation of the listening space
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Reverberation, Karaoke And Other Acoustics (AREA)
- Stereophonic System (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、音場制御システム及び音場制御方法の技術分野に属する。
【０００２】
【従来の技術】
近年新しいオーディオメディアの出現により、オーディオ視聴環境が多様化してきている。このため、どのよう視聴環境においても、臨場感に富んだ音響を得るためのシステムが望まれている。
【０００３】
このような臨場感に富んだ音響を得るシステムとして、いわゆる残響付加装置が数多く提案されている。
【０００４】
従来の残響付加装置は、システムの前段に有限長の初期反射音生成部を有し、その出力信号を後段の残響生成部において再生音源に付加する。この残響生成部は、ＩＩＲ（Infinite Impulse Response）タイプのフィルタが並列に配置される形式が一般的に用いられている。更には、これらの他に実音場のインパルス応答を再生音源に畳み込む形式のものも提案されている。
【０００５】
また同様に臨場感に富んだ音響を得るシステムとして、例えばトランスオーラルシステムのような逆フィルタを使用するシステムも提案されている（例えば特許文献１参照）。このトランスオーラルシステムとは、目標とする音空間における受聴者の位置で収録した音を再生音場において受聴することで、目標とする音空間で受聴しているのと同様の臨場感を得ようとするシステムである。
【特許文献１】
特開平２０００−３１４８７号公報
【０００６】
【発明が解決しようとする課題】
しかしながら、従来の残響付加装置は、空間の広さ等の再生音場の特性を考慮することなく一義的に設計者の勘や経験に基づいて残響特性等が決定され提供されていた。そのため各再生音場の音響特性に適した音響再生が行われていなかった。
【０００７】
また、従来のトランスオーラルシステムにおいては、最適位置で受聴した場合には原音場での臨場感が得られるものの、最適位置を少しでも離れると原音場とは異なった臨場感となるといった制御エリアの狭さが問題であった。またとりわけ狭い空間での再生に際しては、厳密な原音場の定位、残響等の制御が要求されるため、精緻な逆フィルタの設計が困難であった。
【０００８】
本発明は、上記の各問題点に鑑みて為されたものであり、その課題の一例としては、それぞれの再生音場の音響特性に応じて音響情報を付加し、原音場と同様の臨場感のある再生音を得る音場制御システム及び音場制御方法等を提供することにある。
【０００９】
【課題を解決するための手段】
上記の課題を解決するために、請求項１に記載の音場制御システムの発明は、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段と、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段と、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段と、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段は、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１０】
また、請求項６に記載の音場制御方法の発明は、特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、差異検出手段が、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出工程と、音源分割手段が、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割工程と、補正手段が、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正工程と、音源合成手段が、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成工程と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定工程においては、前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１１】
また、請求項７に記載の音場制御用プログラムの発明は、コンピュータによって、音場制御を行う音場制御用プログラムであって、前記コンピュータを、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段、として機能させ、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段としての前記コンピュータを、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴として構成する。
【００１２】
また、請求項８に記載の記録媒体の発明は、請求項７に記載の音場制御用プログラムをコンピュータに読み取り可能に記録したことを特徴として構成する。
【００１３】
また、請求項９に記載の音場空間特性判定システムの発明は、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段は、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１４】
また請求項１１に記載の音場空間特性判定方法の発明は、特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定工程においては、前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１５】
また、請求項１２に記載の音場空間特性判定用プログラムの発明は、コンピュータによって、音場空間の特性判定を行う音場空間特性判定用プログラムであって、前記コンピュータを、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、として機能させ、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段としての前記コンピュータを、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴として構成する。
【００１６】
また、請求項１３に記載の記録媒体の発明は、請求項１２に記載の音場空間特性判定用プログラムをコンピュータに読み取り可能に記録したことを特徴として構成する。
【００１７】
【発明の実施の形態】
続いて、上述の人間が体感する音場の広さ感判定方法を利用した音場空間特性判定システム及び音場制御システムについての実施の形態について説明する。
【００１８】
まず、本願に係る音場空間特性判定システムの概要構成について図１及び図２を用いて説明する。
【００１９】
なお、図１は本願に係る音場空間特性判定システム及び音場制御システムの概要構成の一例を示す機能ブロック図であり、図２は、本願に係る音場空間特性判定システム２の細部構成の一例を示す図である。
【００２０】
図１に示すように、本実施形態の音場制御システムＳＦは、特性測定部２と、特性分割部３と、特性判定部４と、差異検出部５と、音源分割部６と、音源補正部７と、音源合成部８と、を含んで構成される。
【００２１】
なお、特性測定部２と、特性分割部３と、特性判定部４とで、音場空間特性判定システムＳＩを構成している。
【００２２】
また、例えば本実施形態の特性測定部２は本発明の特性測定手段を構成し、特性分割部３は本発明の特性分割手段を構成し、特性判定部４は本発明の空間特性判定手段を構成する。さらに、差異検出部５は本発明の差異検出手段を構成し、音源分割部６は本発明の音源分割手段を構成し、音源補正部７は本発明の補正手段、付加情報算出手段及び情報付加手段を構成する。また音源合成部８は本発明の音源合成手段を構成する。
【００２３】
まず、特性測定部２は、音場再生の対象となる音場としての対象音場のインパルス応答を測定するようになっている。特性測定部２の細部の構成については後述する。
【００２４】
特性分割部３は、測定されたインパルス応答を予め設定された周波数帯域毎に複数に分割するとともに、予め設定された経過時間単位毎に複数に分割することによって、時間−周波数軸上でマトリクス状に複数のブロックデータに分割するようになっている。
【００２５】
なお、予め設定された周波数帯域による分割とは、一定の帯域幅に限らず、それぞれ異なる帯域幅で分割する場合を含む。
【００２６】
また同様に、予め設定された経過時間単位とは、一定の時間単位に限らず、それぞれ異なる時間幅で分割する場合を含む。ここで、経過時間とは、直接音の到来時刻を０ｍｓｅｃとした場合の経過時間をいい、それぞれの分割方法については後述する。
【００２７】
特性判定部４は、分割されたブロックデータに基づいて、対象音場の空間の音響特性としての広さ感を判定し一定の評価値として数値化するようになっている。
【００２８】
なお、対象音場とは、受聴者が受聴する再生音場をいう。再生音場の一例としては、受聴者の自宅の一室や自動車内の座席空間が挙げられる。
【００２９】
また、広さ感とは、音場の広狭に対する感覚として人間が体感する音場の空間の音響特性の一つをいう。その数値特性とは、音場の容積等によって異なる残響を聴くことで、人間がその音場の広さを感覚的に認識する度合いを評価値として数値化されたものをいう。この広さ感を表す評価値の算出方法については後述する。
【００３０】
差異検出部５は、判定された対象音場の広さ感を表す評価値と、所望する音場の空間特性として予め設定された所望空間特性としての所望する音場の広さ感を表す評価値と、を比較してその差異を検出するようになっている。
【００３１】
なお、所望する音場とは、例えばコンサートホールやライブハウス等、受聴者が音楽等を受聴するにあたって、そこで聴いているかのような臨場感を得ることを希望する音場をいう。一般的に、対象音場は所望する音場よりより狭い空間である。
【００３２】
音源分割部６は、再生対象としての音源１を所定の周波数帯域毎に複数の音源成分に分割するようになっている。なお分割される周波数帯域は、上記の音場空間特性判定システムＳＩにおける分割と同じ帯域である。
【００３３】
音源補正部７は、上記周波数成分の分割数と同じ数の反射音追加フィルタを備え、差異検出部５により検出された対象音場の広さ感を示す評価値と所望する音場の広さ感を示す評価値との差異に基づいて、評価値が等しくなるように音源成分を補正するようになっている。
【００３４】
音源合成部８は、反射音パターンが音源成分に付加された周波数帯域毎の信号を再合成して、出力装置ＯＡに対して合成信号を出力するようになっている。
【００３５】
次に、特性測定部２の細部構成について図２を用いて説明する。
【００３６】
特性測定部２は、対象音場のインパルス応答を測定するため、例えば図２に示すように対象音場としての音場１０の中に配設される。
【００３７】
図２に示すように、特性測定部２は、測定用信号音を拡声するスピーカ２１と、スピーカ２１から一定距離離れた位置に設けられたマイクロフォン２２と、スピーカ２１を駆動するテスト信号を発生するテスト信号発生器２３と、スピーカ２１を駆動する増幅器２４と、マイクロフォン２２で受信した音のレベルを調整する増幅器２５と、テスト信号発生器２３のテスト信号とマイクロフォン２２で収音した音波信号からインパルス応答を算出するアナライザ２６と、算出されたインパルス応答を記録する記録装置２７と、を含んで構成されている。
【００３８】
なお、例えば、本実施形態のスピーカ２１は、本発明の拡声手段を構成し、テスト信号発生器２３は本発明の信号発生手段を構成する。さらに、マイクロフォン２２は本発明の収音手段を構成し、アナライザ２６は本発明の応答算出手段を構成する。
【００３９】
続いて、特性測定部２におけるインパルス応答の測定方法について説明する。
【００４０】
まず、テスト信号発生器２３はインパルス信号等のテスト信号（試験信号）を発生し、テスト信号発生器２３からのテスト信号を増幅器２４で所定のレベルにまで増幅調整し、その出力をスピーカ２１に入力させて音を音場１０内に発生させる。
【００４１】
マイクロフォン２２は音場１０内に発せられた音を収音する。マイクロフォン２２は、スピーカ２１に対して予め設定された位置に配設される。
【００４２】
スピーカ２１からの音は図２に示すように、直接にマイクロフォン２２へ到達する直接音と、音場１０の壁、天井、床等に反射してから到達する反射音（間接音）とがあり、これらがインパルス応答の特徴を示すことになる。
【００４３】
増幅器２５は、マイクロフォン２２によって収音された音の信号を所定のレベルに増幅し、アナライザ２６は増幅された信号からインパルス応答を算出し記録装置２７へ記録する。
【００４４】
なお、本実施形態では、アナライザ２６や記録装置２７を用いているが、アナライザ２６や記録装置２７に代えて測定プログラムが備わったパーソナルコンピュータ等を用いることも可能である。
【００４５】
また、本実施形態ではテスト信号（試験信号）としてインパルス信号を用いているが、入力する信号はインパルス信号に限るものではなく、例えばＭ系列雑音やタイムストレッチドパルス信号等を用いることも可能である。
【００４６】
以上のようにして音場のインパルス応答の測定が行われ、測定されたインパルス応答は音場の広さ感を示す数値特性としての評価値の算出に用いられることとなる。
【００４７】
続いて、音場空間特性判定システムＳＩにおいて行われる、音場の広さ感を示す数値特性としての評価値の算出原理について説明する。
【００４８】
まず、例えばさまざまな部屋のインパルス応答を用いて主観評価実験を行うことによりこの評価値を算出する。
【００４９】
具体的には、容積等の異なる複数の部屋で、図２で示した構成で特性測定部２により測定したインパルス応答を音源に畳み込むことにより複数の刺激音を作成する。なお、再生時における各刺激音の音色の変化を最小限に抑えるため、全てのインパルス応答は可聴帯域で周波数特性がほぼフラットになるように補正しておく。
【００５０】
この主観評価実験では、例えばこれら複数の刺激音のうち２つの音を一対にして、ヘッドフォンを通してランダムに再生するようになっており、各音に対する複数の被験者の聴感上の音の広さに対するイメージを得るようになっている。
【００５１】
なお、音源による結果の相違を考慮するため、例えば性質がまったく異なる音声（以下、「ソース１」という）とドラム（以下、「ソース２」という）の２種類のソースを用いている。この実験結果に対して、サーストンの比較判断の原則（ケースＶ）を適用して評価値を算出する。
【００５２】
ここで、サーストンの比較判断の原則（ケースＶ）とは、間隔尺度を用いて評価値を算出する原則である。この間隔尺度とは、心理学上感覚的に等しい距離を数値的に等しい距離で表す値である。したがって、原点が設けられておらず、その値の各差に対して意義を有するものである。
【００５３】
この算出された評価値の一例を図３に示す。
【００５４】
なお、図３は上記実験により得られた評価値の一例を示す表である。
【００５５】
図３に示すように、ソース１、ソース２における各ルーム毎の評価値の傾向は近似している。
【００５６】
次に、以上のようにして得られた評価値データに基づいて、対象音場の広さ感を示す評価値を算出する方法について説明する。
【００５７】
まず、特性分割部３では、特性測定部２で測定された音場のインパルス応答を、時間−周波数軸上で任意の周波数帯域及び経過時間単位で複数のブロックデータにマトリクス状に分割する。
【００５８】
ここで、インパルス応答を時間−周波数軸でそれぞれ分割する理由について図４を用いて説明する。
【００５９】
なお、図４はインパルス応答から算出される残響特性と周波数帯域との関係の一例を示した図である。このうち図４（Ａ）は低い周波数帯域における残響特性の一例を示した図であり、図４（Ｂ）は高い周波数帯域における残響特性の一例を示した図である。
【００６０】
また、音には、前述のように直接音と反射音とが含まれているが、インパルス応答から算出される残響特性とは、反射音の減衰に係わる特性である。
【００６１】
図４（Ａ）（Ｂ）に示すように、一定時刻Ｔを過ぎると時間が経過するほど減衰する残響特性を有する。
【００６２】
また、図４（Ａ）に示す高い周波数帯域では、図４（Ｂ）に示す低い周波数帯域に比べて減衰が早い残響特性を有することとなる。
【００６３】
このように、インパルス応答から算出される残響特性は、経過時間及び周波数帯域により異なる。このため、本実施形態における特性分割部３は、インパルス応答を時間−周波数軸で複数のブロックデータに分割する。
【００６４】
ここで、特性分割部３で行われる時間周波数分割の一例について図５を用いて説明する。
【００６５】
なお、図５は、時間周波数分割の一例を示した図である。
【００６６】
特性分割部３は、図５に示すように、例えば時間軸方向に４分割、周波数軸方向に３分割の計１２ブロックにマトリクス状に分割する。
【００６７】
まず、周波数軸方向については、特性分割部３は、例えば人間の音声帯域（５００Ｈｚ〜２ｋＨｚ）を基準として３分割する。
【００６８】
また、時間軸方向については、特性分割部３は、例えば人間が聴いたときに音が分離して聞こえるか否かの経過時刻として一般的に知られている８０ｍｓｅｃを基準として４分割する。
【００６９】
なお、図５において、時間軸の０ｍｓｅｃは、直接音が到来する時刻を示している。
【００７０】
次に、特性判定部４は、特性分割部３で分割された各ブロックデータのエネルギーを算出する。各ブロックデータに含まれるインパルス応答のエネルギーをＥijとして、特性判定部４で行われるＥijの算出方法の一例としてＥ22の算出手順を以下に示す。
【００７１】
まず特性判定部４は、解析対象となるインパルス応答ｐ（ｔ）を、Ｅ22は５００Ｈｚから２ｋＨｚまでの成分に対応するバンドパスフィルタｈ2（ｔ）でフィルタリングする。これを帯域幅で正規化し、エネルギー系列ｅ（ｔ）を算出する。
【００７２】
なお、バンドパスフィルタｈ2（ｔ）としては、例えば直線位相ＦＩＲ（Finite Impulse Response）を設定する。
【００７３】
【数１】
特性判定部４は、ｅ（ｔ）を時間軸上で積分し、かつ時間幅で正規化しＥ22を算出する。ここで、ｔ0は直接音が到来する時刻とする。
【００７４】
【数２】
特性判定部４は、Ｅ22以外の他のブロックデータのエネルギーについても同様に算出し、算出された後に、Ｅijを正規化する。
【００７５】
【数３】
このように、ある部屋のインパルス応答ｐ（ｔ）が得られているとき、エネルギーＥijが得られる。広さ感を示す評価値Ｒが重み係数（回帰係数）Ｗijを用いてＥijの重み付きの線形和（加重線形和）で表すとすると、
【００７６】
【数４】
となる。
【００７７】
上記の実験等で得られた評価値Ｒを目的変数とし、対応するインパルス応答のエネルギーＥijを説明変数として重回帰分析が行われ、重み係数Ｗijが算出される。
【００７８】
このように重み係数Ｗijは、評価値Ｒの算出に対する寄与度を示す係数としてブロックデータ毎に予め定められる。この重み係数Ｗijがプラスの場合、そのブロックデータのエネルギー値が増加すると評価値Ｒは大きくなり、反対に重み係数Ｗijがマイナスのブロックデータのエネルギー値が増加すると、評価値Ｒは小さくなる。
【００７９】
ここで、この重み係数Ｗijの一例について図６及び図７を用いて説明する。
【００８０】
なお、図６は実験結果に基づくソース１についての重み係数の一例を示したグラフであり、図７は実験結果に基づくソース２についての重み係数の一例を示したグラフである。
【００８１】
図６及び図７に示すように、いずれのソースに対する重み係数も同様の形状を示すグラフとなっており、各ブロックデータのエネルギー値の増減による全体の広さ感Ｒの変化の度合いは、各ブロックデータで夫々異なっている。その中でも周波数２ｋＨｚ超の帯域かつ経過時間８０ｍｓｅｃ−１６０ｍｓｅｃのブロックにおけるＥ43の重み係数の値が大きくなっていることが示される。
【００８２】
これは、重み係数の高いＥ43の成分値の増減により全体としての広さ感を示す評価値Ｒが大きく左右されることを意味している。
【００８３】
以上のように、広さ感を示す評価値は、時間周波数平面での離散的エネルギー分布の加重線形和により（式４）で計算される。また、反射音の時間周波数構造を任意に変更し、この評価値を制御することによって人間が体感する広さ感を変更することが可能となる。
【００８４】
そこで次に、この原理をふまえて本実施形態における音源補正部７において行われる具体的な音源成分の補正方法について説明する。
【００８５】
上述のように、重み係数Ｗijがプラスの場合、そのブロックデータのエネルギー値が増加すると評価値Ｒは大きくなり、反対に重み係数Ｗijがマイナスのブロックデータのエネルギー値が増加すると、評価値Ｒは小さくなる。
【００８６】
したがって、音源補正部７は、音場１０の広さ感を示す評価値Ｒsmallを大きく、すなわちより広い音場の残響特性に変更する場合には、重み係数Ｗijがプラスのブロックのエネルギー値を引き上げるように付加情報を算出し、対応する帯域の反射音追加フィルタの修正を行い、音源成分の補正を行う。
【００８７】
一方、例えば修正により音場１０の広さ感を示す評価値評価値Ｒsmallが大きくなりすぎて、より狭い音場の残響特性に変更する場合には、音源補正部７は、重み係数Ｗijがマイナスのブロックのエネルギー値を引き上げるように付加情報を算出し、対応する帯域の反射音追加フィルタの修正を行い、音源成分の補正を行う。
【００８８】
また、具体的な付加情報の決定方法としては、音源補正部７は、例えば目標とする音場のエネルギー分布と対象音場のエネルギー分布とをブロック単位で比較し、差分の大きいブロックから優先的に補正するように付加情報を決定する。
【００８９】
なお、音源補正部７は、補正後の音がより自然な音となるように付加情報の算出及びフィルタの修正を行う。
【００９０】
すなわち、自然音は時間が経過するにつれてエネルギーが減衰する性質を有しており、経験的には指数的に減衰していくのが人間の耳にとって自然な音となる。また周波数が高いほど、減衰曲線のカーブは鋭くなる。そこで、音源補正部７は、これらのパターンに沿う形で音源成分を補正する。
【００９１】
なお、音源補正部７は、付加情報の算出及びフィルタの修正の際には、上記ブロック全てを対象とせず、上記ブロックの中でも重み係数が大きく広さ感への影響が大きいＥ43を対象とし、Ｅ43のエネルギー値を上げるように音源成分を補正することも可能である。
【００９２】
また、同様の理由により、音源補正部７は、例えば経過時間が所定の時間、本実施形態においては８０ｍｓｅｃを超える経過時間に属するＥ43及びＥ42のエネルギー値を上げるように音源成分を補正することも可能である。
【００９３】
続いて、上述の音場空間特性判定システムによる広さ感判定及び音源成分の補正方法を用いた音場制御処理の動作について図８を用いて説明する。
【００９４】
なお、図８は、実施形態に係る音場制御システムＳＦによって行われる音場制御処理の動作を示すフローチャートである。
【００９５】
まず、図示しない制御部から音場制御指示を受けた場合、音源補正部７は、再生を行う音場１０のインパルス応答ｈsmallを設定する（ステップＳ１１）。
【００９６】
次いで、差異検出部５は、例えばコンサートホール等の所望する音場の広さ感を示す評価値Ｒlargeを設定する（ステップＳ１２）。
【００９７】
次いで、音源補正部７は、分割する周波数帯域数分、ここでは３つの帯域の反射音付加フィルタを設定する（ステップＳ１３）。
【００９８】
次いで、特性判定部４は、再生を行う音場１０のインパルス応答ｈsmallとこれらを合成したフィルタの畳み込み計算を行い、値Ｈを計算する（ステップＳ１４）。
【００９９】
【数５】
次いで、特性判定部４は、前述の広さ感を示す評価値の計算式（式４）を用いて畳み込み演算後の音場１０の広さ感を示す評価値Ｒsmallを算出し、算出結果を差異検出部５へ出力する（ステップＳ１５）。
【０１００】
次いで、差異検出部５は、特性判定部４で算出された音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとの差異を検出し、値が一致するか否かの判断を行う（ステップＳ１６）。
【０１０１】
差異検出部５が音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとが一致しないと判断した場合には（ステップＳ１６、Ｎｏ）、音源補正部７は、差異検出部５において検出された差異に基づいて、反射音追加フィルタのいずれかを微修正すべく音源成分に付加する情報を算出し、音源補正部７の反射音追加フィルタの修正を行う（ステップＳ１７）。この音源補正部７における付加情報の算出及び反射音追加フィルタの修正方法については前述したとおりである。
【０１０２】
なお、フィルタの修正を行った後は、ステップＳ１４の処理へと戻り、差異検出部５は、Ｒlargeと修正後のＲsmallとの差異の検出を行う。
【０１０３】
一方、ステップＳ１６の判断において、差異検出部５が音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとが一致すると判断した場合には（ステップＳ１６、Ｙｅｓ）、一連の処理を終了する（ステップＳ１８）。
【０１０４】
これにより、一連の音場制御処理の終了後、修正された反射音追加フィルタを用いて例えば音楽の再生が開始されることによって、受聴者は所望する広さ感を体感することが可能となる。
【０１０５】
なお、インパルス応答ｈsmallのデータは、例えば前述した測定方法によりあらかじめ自動車内の座席空間等の対象音場で測定しておき、記録装置２７に記録したデータを用いることが可能である。また所望する音場の広さ感を示す評価値Ｒlargeのデータについても同様に、あらかじめ記録装置２７に記録しておき、これを用いることが可能である。
【０１０６】
以上説明したように、本実施形態によれば、音場再生の対象となる音場としての音場１０のインパルス応答を測定する特性測定部２と、測定されたインパルス応答を、予め設定された周波数帯域毎に複数に分割するとともに、予め設定された経過時間単位毎に複数に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割部３と、分割されたブロックデータ夫々に基づいて、音場１０における空間の音響特性としての対象空間特性を判定する特性判定部４と、判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出部５と、音場１０で受聴に供する音源１を予め設定された周波数帯域毎に音源成分を分割する音源分割部６と、検出された空間特性の差異に基づいて、音源成分毎に少なくとも１の音源成分を補正する音源補正部７と、補正された各音源成分に基づいて音源を再合成する音源合成部８と、を備えた構成を有している。
【０１０７】
したがって、対象音場のインパルス応答を測定し、予め設定された周波数及び経過時間で複数のブロックに分割することによって対象音場の空間特性を判定し、所望する空間特性との差異を検出するとともに、その検出結果に応じて、対象音場で受聴に供する音に対して少なくともいずれか１の前記ブロックの音源を補正する。
【０１０８】
このように、周波数及び経過時間によって異なる反射音パターンを複数のブロックに分けることにより対象音場の空間特性、例えば人間が体感する音場の広さ感について判定することができるため、音場制御や音響ホール設計に活用することが可能となる音場空間特性判定システムを提供することができる。
【０１０９】
また、所望する音場との空間特性の差異に応じて、音源に付加すべき情報を算出し、所定のブロック単位で再生すべき音源の補正を行うので、例えば狭い受聴空間においても広いコンサートホールなど所望する音場で受聴した場合と同様の臨場感を容易に得ることが可能となる音場制御システムを提供することができる。
【０１１０】
また本実施形態によれば、空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、特性判定部４は、夫々のブロックデータのエネルギー値Ｅijと各ブロックデータ毎に定められた重み係数Ｗijとの加重線形和に基づいて空間特性Ｒを判定することを特徴として構成されている。
【０１１１】
したがって、周波数及び時間によりブロックごとに異なる広さ感への影響度を重み付けを行うことによって計算し判定するので、人間が体感する音場の広さ感を客観的な評価値として表すことが可能となる音場空間特性判定システムを提供することができる。
【０１１２】
また、この評価値の差異に基づいてブロック単位で音場制御を行えばよいので、従来必要とされたような厳密なフィルタ制御を行う必要がなく、短時間にフィルタを設計することが可能となり、また分割する帯域数を減らす効果もあり実時間処理を大幅に減らすことが可能となる。
【０１１３】
さらに評価値の差異に基づいて音場制御を行うので、設計者の勘や経験に頼ることなく、客観的な音場制御を行うことが可能な音場制御システムを提供することができる。
【０１１４】
また本実施形態によれば、特性測定部２は、試験信号に基づいて音を拡声するスピーカ２１と、試験信号を発生させるテスト信号発生器２３と、スピーカ２１から発せられる音を収音するマイクロフォン２２と、試験信号と、収音された音とに基づいて、スピーカ２１とマイクロフォン２２との間のインパルス応答を算出するアナライザ２６と、を含む構成を有している。
【０１１５】
したがって、音場空間特性判定システムＳＩにおいてインパルス応答を算出することにより、対象音場の空間特性を容易に測定することが可能な音場空間特性判定システム及びこれを利用した音場制御システムを提供することができる。
【０１１６】
また本実施形態によれば、予め設定された周波数帯域に、予め定められた低中域を含む帯域と、当該低中域を超過する帯域とを少なくとも含む場合に、音源補正部７は、当該低中域を超過する帯域に含まれる音源成分を補正する構成を有している。
【０１１７】
したがって、人間の聴覚のメカニズム上、音場の広さに対する感覚に対して重要な要素である中高域での反射音成分を補正するので、所望する広さ感を効果的に得ることが可能な音場制御システムを提供することができる。
【０１１８】
また本実施形態によれば、周波数軸において分割するために予め設定された周波数としての低中域の値は、５００Ｈｚ〜２ｋＨｚとする構成を有している。
【０１１９】
したがって、低中域の上限として一般的に認識されている値である２ｋＨｚを超える周波数帯域で反射音成分を補正するので、所望する広さ感をより効果的に得ることが可能な音場制御システムを提供することができる。
【０１２０】
また本実施形態によれば、音源補正部７は、直接音が到来する時刻を基準として、当該時刻から予め定められた時間を経過した経過時間に含まれる反射音成分を補正することを特徴とする構成を有している。
【０１２１】
したがって、人間の聴覚のメカニズム上、音場の広さに対する感覚に対して重要な要素である所定時刻（例えば８０ｍｓｅｃ）を超える経過時間で音源に反射音を付加するので、所望する広さ感を効果的に得ることが可能な音場制御システムを提供することができる。
【０１２２】
なお、この所定時刻は例えばクラシックやロック等再生する音楽のジャンルにより任意に変更することができるので、ジャンル等に応じてさらに効果的な音場制御を行うことが可能となる。
【０１２３】
また、本実施形態によれば、音源補正部７は、検出された空間特性の差異に基づいて、少なくともいずれか１の音源成分に付加する付加情報を算出し、当該音源成分に対して算出された付加情報を付加することによって補正を行うことを特徴とする構成を有している。
【０１２４】
したがって、補正方法としては付加情報を算出し音源成分に付加するだけでよいので、従来必要とされたような厳密なフィルタ制御を行う必要がなく、短時間にフィルタを設計することが可能となり、実時間処理を更に大幅に減らすことが可能な音場制御システムを提供することができる。
【０１２５】
また、本実施形態では、音場空間特性判定システムＳＩを対象音場において所望する音場の広さ感を得ることができるように音源に情報を付加し調整する音場制御システムＳＦに適用する構成をとっているが、本判定システムの適用はこれに限るものではない。すなわち、例えば音響ホールやリスニングルーム等再生空間の設計の分野に対しても当該音場空間特性判定システムは適用することができる。またその他空間印象の尺度として本広さ感判定を利用して各種設計及び制御を行うあらゆる分野に対して本判定システムは適用可能である。
【０１２６】
また、本実施形態では、上述のシステム構成によって音場空間特性の判定及び音場の制御を行うようになっているが、システム中にコンピュータおよび記録媒体を備え、この記録媒体に上述の広さ感判定動作及び音場制御動作を行うプログラムを格納し、コンピュータによって当該広さ感判定プログラム及び音場制御プログラムを読み込むことによって上述と同様の動作を行うようにすることも可能である。
【０１２７】
また、この音場空間特性判定プログラム及び音場制御プログラムを実行するシステムにおいて、記録媒体をＤＶＤやＣＤなどの記録媒体により構成するようにしてもよい。
【０１２８】
この場合、当該システムは、記録媒体からプログラムを読み出す読出装置を備えることとなる。
【０１２９】
また、本実施形態では、一つのスピーカに対して一箇所の受音点がある場合を想定しているが、これに限るものではなく、２スピーカを備えたステレオシステムやマルチチャンネルシステムに対しても適用可能である。例えば複数のスピーカを所定の位置において再生する場合には、それぞれのスピーカについて個別にインパルス応答の測定を行い、図８におけるステップＳ１３乃至１７で説明した処理を音源の数だけ繰り返し行えばよい。
【０１３０】
なお、本実施の形態においては、時間軸方向に４分割、周波数軸方向に３分割の計１２ブロックに分割しているが、これに限るものではなく、時間−周波数軸における分割方法及び分割数は任意である。したがって、例えば再生する音楽のジャンル、使用するフィルタ等さまざまな条件を考慮して適宜時間軸方向の分割数及び周波数軸方向の分割数や分割の基準とする時間・周波数を決定することが可能である。
【図面の簡単な説明】
【図１】本願に係る音場空間特性判定システム及び音場制御システムの概要構成の一例を示す機能ブロック図である。
【図２】本願に係る音場空間特性判定システム２の細部構成の一例を示す図である。
【図３】広さ感を示す評価値の一例を示す表である。
【図４】インパルス応答から算出される残響特性と周波数帯域との関係の一例を示した図であり、このうち（Ａ）は低い周波数帯域における残響特性の一例を示した図であり、（Ｂ）は高い周波数帯域における残響特性の一例を示した図である。
【図５】時間周波数分割の一例を示した図である。
【図６】実験結果に基づくソース１についての重み係数の一例を示したグラフである。
【図７】実験結果に基づくソース２についての重み係数の一例を示したグラフである。
【図８】実施形態に係る音場制御システムＳＦによって行われる音場制御処理の動作を示すフローチャートである。
【符号の説明】
ＳＦ … 音場制御システム
ＳＩ … 音場空間特性判定システム
１ … 音源
２ … 特性測定部
３ … 特性分割部
４ … 特性判定部
５ … 差異検出部
６ … 音源分割部
７ … 音源補正部
８ … 音源合成部
１０ … 音場
２１ … スピーカ
２２ … マイクロフォン
２３ … テスト信号発生器
２４ … 増幅器
２５ … 増幅器
２６ … アナライザ
２７ … 記録装置
【発明の属する技術分野】
本発明は、音場制御システム及び音場制御方法の技術分野に属する。
【０００２】
【従来の技術】
近年新しいオーディオメディアの出現により、オーディオ視聴環境が多様化してきている。このため、どのよう視聴環境においても、臨場感に富んだ音響を得るためのシステムが望まれている。
【０００３】
このような臨場感に富んだ音響を得るシステムとして、いわゆる残響付加装置が数多く提案されている。
【０００４】
従来の残響付加装置は、システムの前段に有限長の初期反射音生成部を有し、その出力信号を後段の残響生成部において再生音源に付加する。この残響生成部は、ＩＩＲ（Infinite Impulse Response）タイプのフィルタが並列に配置される形式が一般的に用いられている。更には、これらの他に実音場のインパルス応答を再生音源に畳み込む形式のものも提案されている。
【０００５】
また同様に臨場感に富んだ音響を得るシステムとして、例えばトランスオーラルシステムのような逆フィルタを使用するシステムも提案されている（例えば特許文献１参照）。このトランスオーラルシステムとは、目標とする音空間における受聴者の位置で収録した音を再生音場において受聴することで、目標とする音空間で受聴しているのと同様の臨場感を得ようとするシステムである。
【特許文献１】
特開平２０００−３１４８７号公報
【０００６】
【発明が解決しようとする課題】
しかしながら、従来の残響付加装置は、空間の広さ等の再生音場の特性を考慮することなく一義的に設計者の勘や経験に基づいて残響特性等が決定され提供されていた。そのため各再生音場の音響特性に適した音響再生が行われていなかった。
【０００７】
また、従来のトランスオーラルシステムにおいては、最適位置で受聴した場合には原音場での臨場感が得られるものの、最適位置を少しでも離れると原音場とは異なった臨場感となるといった制御エリアの狭さが問題であった。またとりわけ狭い空間での再生に際しては、厳密な原音場の定位、残響等の制御が要求されるため、精緻な逆フィルタの設計が困難であった。
【０００８】
本発明は、上記の各問題点に鑑みて為されたものであり、その課題の一例としては、それぞれの再生音場の音響特性に応じて音響情報を付加し、原音場と同様の臨場感のある再生音を得る音場制御システム及び音場制御方法等を提供することにある。
【０００９】
【課題を解決するための手段】
上記の課題を解決するために、請求項１に記載の音場制御システムの発明は、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段と、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段と、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段と、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段は、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１０】
また、請求項６に記載の音場制御方法の発明は、特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、差異検出手段が、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出工程と、音源分割手段が、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割工程と、補正手段が、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正工程と、音源合成手段が、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成工程と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定工程においては、前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１１】
また、請求項７に記載の音場制御用プログラムの発明は、コンピュータによって、音場制御を行う音場制御用プログラムであって、前記コンピュータを、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段、として機能させ、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段としての前記コンピュータを、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴として構成する。
【００１２】
また、請求項８に記載の記録媒体の発明は、請求項７に記載の音場制御用プログラムをコンピュータに読み取り可能に記録したことを特徴として構成する。
【００１３】
また、請求項９に記載の音場空間特性判定システムの発明は、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段は、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１４】
また請求項１１に記載の音場空間特性判定方法の発明は、特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、を備え、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定工程においては、前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴として構成する。
【００１５】
また、請求項１２に記載の音場空間特性判定用プログラムの発明は、コンピュータによって、音場空間の特性判定を行う音場空間特性判定用プログラムであって、前記コンピュータを、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、として機能させ、前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、前記空間特性判定手段としての前記コンピュータを、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴として構成する。
【００１６】
また、請求項１３に記載の記録媒体の発明は、請求項１２に記載の音場空間特性判定用プログラムをコンピュータに読み取り可能に記録したことを特徴として構成する。
【００１７】
【発明の実施の形態】
続いて、上述の人間が体感する音場の広さ感判定方法を利用した音場空間特性判定システム及び音場制御システムについての実施の形態について説明する。
【００１８】
まず、本願に係る音場空間特性判定システムの概要構成について図１及び図２を用いて説明する。
【００１９】
なお、図１は本願に係る音場空間特性判定システム及び音場制御システムの概要構成の一例を示す機能ブロック図であり、図２は、本願に係る音場空間特性判定システム２の細部構成の一例を示す図である。
【００２０】
図１に示すように、本実施形態の音場制御システムＳＦは、特性測定部２と、特性分割部３と、特性判定部４と、差異検出部５と、音源分割部６と、音源補正部７と、音源合成部８と、を含んで構成される。
【００２１】
なお、特性測定部２と、特性分割部３と、特性判定部４とで、音場空間特性判定システムＳＩを構成している。
【００２２】
また、例えば本実施形態の特性測定部２は本発明の特性測定手段を構成し、特性分割部３は本発明の特性分割手段を構成し、特性判定部４は本発明の空間特性判定手段を構成する。さらに、差異検出部５は本発明の差異検出手段を構成し、音源分割部６は本発明の音源分割手段を構成し、音源補正部７は本発明の補正手段、付加情報算出手段及び情報付加手段を構成する。また音源合成部８は本発明の音源合成手段を構成する。
【００２３】
まず、特性測定部２は、音場再生の対象となる音場としての対象音場のインパルス応答を測定するようになっている。特性測定部２の細部の構成については後述する。
【００２４】
特性分割部３は、測定されたインパルス応答を予め設定された周波数帯域毎に複数に分割するとともに、予め設定された経過時間単位毎に複数に分割することによって、時間−周波数軸上でマトリクス状に複数のブロックデータに分割するようになっている。
【００２５】
なお、予め設定された周波数帯域による分割とは、一定の帯域幅に限らず、それぞれ異なる帯域幅で分割する場合を含む。
【００２６】
また同様に、予め設定された経過時間単位とは、一定の時間単位に限らず、それぞれ異なる時間幅で分割する場合を含む。ここで、経過時間とは、直接音の到来時刻を０ｍｓｅｃとした場合の経過時間をいい、それぞれの分割方法については後述する。
【００２７】
特性判定部４は、分割されたブロックデータに基づいて、対象音場の空間の音響特性としての広さ感を判定し一定の評価値として数値化するようになっている。
【００２８】
なお、対象音場とは、受聴者が受聴する再生音場をいう。再生音場の一例としては、受聴者の自宅の一室や自動車内の座席空間が挙げられる。
【００２９】
また、広さ感とは、音場の広狭に対する感覚として人間が体感する音場の空間の音響特性の一つをいう。その数値特性とは、音場の容積等によって異なる残響を聴くことで、人間がその音場の広さを感覚的に認識する度合いを評価値として数値化されたものをいう。この広さ感を表す評価値の算出方法については後述する。
【００３０】
差異検出部５は、判定された対象音場の広さ感を表す評価値と、所望する音場の空間特性として予め設定された所望空間特性としての所望する音場の広さ感を表す評価値と、を比較してその差異を検出するようになっている。
【００３１】
なお、所望する音場とは、例えばコンサートホールやライブハウス等、受聴者が音楽等を受聴するにあたって、そこで聴いているかのような臨場感を得ることを希望する音場をいう。一般的に、対象音場は所望する音場よりより狭い空間である。
【００３２】
音源分割部６は、再生対象としての音源１を所定の周波数帯域毎に複数の音源成分に分割するようになっている。なお分割される周波数帯域は、上記の音場空間特性判定システムＳＩにおける分割と同じ帯域である。
【００３３】
音源補正部７は、上記周波数成分の分割数と同じ数の反射音追加フィルタを備え、差異検出部５により検出された対象音場の広さ感を示す評価値と所望する音場の広さ感を示す評価値との差異に基づいて、評価値が等しくなるように音源成分を補正するようになっている。
【００３４】
音源合成部８は、反射音パターンが音源成分に付加された周波数帯域毎の信号を再合成して、出力装置ＯＡに対して合成信号を出力するようになっている。
【００３５】
次に、特性測定部２の細部構成について図２を用いて説明する。
【００３６】
特性測定部２は、対象音場のインパルス応答を測定するため、例えば図２に示すように対象音場としての音場１０の中に配設される。
【００３７】
図２に示すように、特性測定部２は、測定用信号音を拡声するスピーカ２１と、スピーカ２１から一定距離離れた位置に設けられたマイクロフォン２２と、スピーカ２１を駆動するテスト信号を発生するテスト信号発生器２３と、スピーカ２１を駆動する増幅器２４と、マイクロフォン２２で受信した音のレベルを調整する増幅器２５と、テスト信号発生器２３のテスト信号とマイクロフォン２２で収音した音波信号からインパルス応答を算出するアナライザ２６と、算出されたインパルス応答を記録する記録装置２７と、を含んで構成されている。
【００３８】
なお、例えば、本実施形態のスピーカ２１は、本発明の拡声手段を構成し、テスト信号発生器２３は本発明の信号発生手段を構成する。さらに、マイクロフォン２２は本発明の収音手段を構成し、アナライザ２６は本発明の応答算出手段を構成する。
【００３９】
続いて、特性測定部２におけるインパルス応答の測定方法について説明する。
【００４０】
まず、テスト信号発生器２３はインパルス信号等のテスト信号（試験信号）を発生し、テスト信号発生器２３からのテスト信号を増幅器２４で所定のレベルにまで増幅調整し、その出力をスピーカ２１に入力させて音を音場１０内に発生させる。
【００４１】
マイクロフォン２２は音場１０内に発せられた音を収音する。マイクロフォン２２は、スピーカ２１に対して予め設定された位置に配設される。
【００４２】
スピーカ２１からの音は図２に示すように、直接にマイクロフォン２２へ到達する直接音と、音場１０の壁、天井、床等に反射してから到達する反射音（間接音）とがあり、これらがインパルス応答の特徴を示すことになる。
【００４３】
増幅器２５は、マイクロフォン２２によって収音された音の信号を所定のレベルに増幅し、アナライザ２６は増幅された信号からインパルス応答を算出し記録装置２７へ記録する。
【００４４】
なお、本実施形態では、アナライザ２６や記録装置２７を用いているが、アナライザ２６や記録装置２７に代えて測定プログラムが備わったパーソナルコンピュータ等を用いることも可能である。
【００４５】
また、本実施形態ではテスト信号（試験信号）としてインパルス信号を用いているが、入力する信号はインパルス信号に限るものではなく、例えばＭ系列雑音やタイムストレッチドパルス信号等を用いることも可能である。
【００４６】
以上のようにして音場のインパルス応答の測定が行われ、測定されたインパルス応答は音場の広さ感を示す数値特性としての評価値の算出に用いられることとなる。
【００４７】
続いて、音場空間特性判定システムＳＩにおいて行われる、音場の広さ感を示す数値特性としての評価値の算出原理について説明する。
【００４８】
まず、例えばさまざまな部屋のインパルス応答を用いて主観評価実験を行うことによりこの評価値を算出する。
【００４９】
具体的には、容積等の異なる複数の部屋で、図２で示した構成で特性測定部２により測定したインパルス応答を音源に畳み込むことにより複数の刺激音を作成する。なお、再生時における各刺激音の音色の変化を最小限に抑えるため、全てのインパルス応答は可聴帯域で周波数特性がほぼフラットになるように補正しておく。
【００５０】
この主観評価実験では、例えばこれら複数の刺激音のうち２つの音を一対にして、ヘッドフォンを通してランダムに再生するようになっており、各音に対する複数の被験者の聴感上の音の広さに対するイメージを得るようになっている。
【００５１】
なお、音源による結果の相違を考慮するため、例えば性質がまったく異なる音声（以下、「ソース１」という）とドラム（以下、「ソース２」という）の２種類のソースを用いている。この実験結果に対して、サーストンの比較判断の原則（ケースＶ）を適用して評価値を算出する。
【００５２】
ここで、サーストンの比較判断の原則（ケースＶ）とは、間隔尺度を用いて評価値を算出する原則である。この間隔尺度とは、心理学上感覚的に等しい距離を数値的に等しい距離で表す値である。したがって、原点が設けられておらず、その値の各差に対して意義を有するものである。
【００５３】
この算出された評価値の一例を図３に示す。
【００５４】
なお、図３は上記実験により得られた評価値の一例を示す表である。
【００５５】
図３に示すように、ソース１、ソース２における各ルーム毎の評価値の傾向は近似している。
【００５６】
次に、以上のようにして得られた評価値データに基づいて、対象音場の広さ感を示す評価値を算出する方法について説明する。
【００５７】
まず、特性分割部３では、特性測定部２で測定された音場のインパルス応答を、時間−周波数軸上で任意の周波数帯域及び経過時間単位で複数のブロックデータにマトリクス状に分割する。
【００５８】
ここで、インパルス応答を時間−周波数軸でそれぞれ分割する理由について図４を用いて説明する。
【００５９】
なお、図４はインパルス応答から算出される残響特性と周波数帯域との関係の一例を示した図である。このうち図４（Ａ）は低い周波数帯域における残響特性の一例を示した図であり、図４（Ｂ）は高い周波数帯域における残響特性の一例を示した図である。
【００６０】
また、音には、前述のように直接音と反射音とが含まれているが、インパルス応答から算出される残響特性とは、反射音の減衰に係わる特性である。
【００６１】
図４（Ａ）（Ｂ）に示すように、一定時刻Ｔを過ぎると時間が経過するほど減衰する残響特性を有する。
【００６２】
また、図４（Ａ）に示す高い周波数帯域では、図４（Ｂ）に示す低い周波数帯域に比べて減衰が早い残響特性を有することとなる。
【００６３】
このように、インパルス応答から算出される残響特性は、経過時間及び周波数帯域により異なる。このため、本実施形態における特性分割部３は、インパルス応答を時間−周波数軸で複数のブロックデータに分割する。
【００６４】
ここで、特性分割部３で行われる時間周波数分割の一例について図５を用いて説明する。
【００６５】
なお、図５は、時間周波数分割の一例を示した図である。
【００６６】
特性分割部３は、図５に示すように、例えば時間軸方向に４分割、周波数軸方向に３分割の計１２ブロックにマトリクス状に分割する。
【００６７】
まず、周波数軸方向については、特性分割部３は、例えば人間の音声帯域（５００Ｈｚ〜２ｋＨｚ）を基準として３分割する。
【００６８】
また、時間軸方向については、特性分割部３は、例えば人間が聴いたときに音が分離して聞こえるか否かの経過時刻として一般的に知られている８０ｍｓｅｃを基準として４分割する。
【００６９】
なお、図５において、時間軸の０ｍｓｅｃは、直接音が到来する時刻を示している。
【００７０】
次に、特性判定部４は、特性分割部３で分割された各ブロックデータのエネルギーを算出する。各ブロックデータに含まれるインパルス応答のエネルギーをＥijとして、特性判定部４で行われるＥijの算出方法の一例としてＥ22の算出手順を以下に示す。
【００７１】
まず特性判定部４は、解析対象となるインパルス応答ｐ（ｔ）を、Ｅ22は５００Ｈｚから２ｋＨｚまでの成分に対応するバンドパスフィルタｈ2（ｔ）でフィルタリングする。これを帯域幅で正規化し、エネルギー系列ｅ（ｔ）を算出する。
【００７２】
なお、バンドパスフィルタｈ2（ｔ）としては、例えば直線位相ＦＩＲ（Finite Impulse Response）を設定する。
【００７３】
【数１】
【００７４】
【数２】
【００７５】
【数３】
【００７６】
【数４】
【００７７】
上記の実験等で得られた評価値Ｒを目的変数とし、対応するインパルス応答のエネルギーＥijを説明変数として重回帰分析が行われ、重み係数Ｗijが算出される。
【００７８】
このように重み係数Ｗijは、評価値Ｒの算出に対する寄与度を示す係数としてブロックデータ毎に予め定められる。この重み係数Ｗijがプラスの場合、そのブロックデータのエネルギー値が増加すると評価値Ｒは大きくなり、反対に重み係数Ｗijがマイナスのブロックデータのエネルギー値が増加すると、評価値Ｒは小さくなる。
【００７９】
ここで、この重み係数Ｗijの一例について図６及び図７を用いて説明する。
【００８０】
なお、図６は実験結果に基づくソース１についての重み係数の一例を示したグラフであり、図７は実験結果に基づくソース２についての重み係数の一例を示したグラフである。
【００８１】
図６及び図７に示すように、いずれのソースに対する重み係数も同様の形状を示すグラフとなっており、各ブロックデータのエネルギー値の増減による全体の広さ感Ｒの変化の度合いは、各ブロックデータで夫々異なっている。その中でも周波数２ｋＨｚ超の帯域かつ経過時間８０ｍｓｅｃ−１６０ｍｓｅｃのブロックにおけるＥ43の重み係数の値が大きくなっていることが示される。
【００８２】
これは、重み係数の高いＥ43の成分値の増減により全体としての広さ感を示す評価値Ｒが大きく左右されることを意味している。
【００８３】
以上のように、広さ感を示す評価値は、時間周波数平面での離散的エネルギー分布の加重線形和により（式４）で計算される。また、反射音の時間周波数構造を任意に変更し、この評価値を制御することによって人間が体感する広さ感を変更することが可能となる。
【００８４】
そこで次に、この原理をふまえて本実施形態における音源補正部７において行われる具体的な音源成分の補正方法について説明する。
【００８５】
上述のように、重み係数Ｗijがプラスの場合、そのブロックデータのエネルギー値が増加すると評価値Ｒは大きくなり、反対に重み係数Ｗijがマイナスのブロックデータのエネルギー値が増加すると、評価値Ｒは小さくなる。
【００８６】
したがって、音源補正部７は、音場１０の広さ感を示す評価値Ｒsmallを大きく、すなわちより広い音場の残響特性に変更する場合には、重み係数Ｗijがプラスのブロックのエネルギー値を引き上げるように付加情報を算出し、対応する帯域の反射音追加フィルタの修正を行い、音源成分の補正を行う。
【００８７】
一方、例えば修正により音場１０の広さ感を示す評価値評価値Ｒsmallが大きくなりすぎて、より狭い音場の残響特性に変更する場合には、音源補正部７は、重み係数Ｗijがマイナスのブロックのエネルギー値を引き上げるように付加情報を算出し、対応する帯域の反射音追加フィルタの修正を行い、音源成分の補正を行う。
【００８８】
また、具体的な付加情報の決定方法としては、音源補正部７は、例えば目標とする音場のエネルギー分布と対象音場のエネルギー分布とをブロック単位で比較し、差分の大きいブロックから優先的に補正するように付加情報を決定する。
【００８９】
なお、音源補正部７は、補正後の音がより自然な音となるように付加情報の算出及びフィルタの修正を行う。
【００９０】
すなわち、自然音は時間が経過するにつれてエネルギーが減衰する性質を有しており、経験的には指数的に減衰していくのが人間の耳にとって自然な音となる。また周波数が高いほど、減衰曲線のカーブは鋭くなる。そこで、音源補正部７は、これらのパターンに沿う形で音源成分を補正する。
【００９１】
なお、音源補正部７は、付加情報の算出及びフィルタの修正の際には、上記ブロック全てを対象とせず、上記ブロックの中でも重み係数が大きく広さ感への影響が大きいＥ43を対象とし、Ｅ43のエネルギー値を上げるように音源成分を補正することも可能である。
【００９２】
また、同様の理由により、音源補正部７は、例えば経過時間が所定の時間、本実施形態においては８０ｍｓｅｃを超える経過時間に属するＥ43及びＥ42のエネルギー値を上げるように音源成分を補正することも可能である。
【００９３】
続いて、上述の音場空間特性判定システムによる広さ感判定及び音源成分の補正方法を用いた音場制御処理の動作について図８を用いて説明する。
【００９４】
なお、図８は、実施形態に係る音場制御システムＳＦによって行われる音場制御処理の動作を示すフローチャートである。
【００９５】
まず、図示しない制御部から音場制御指示を受けた場合、音源補正部７は、再生を行う音場１０のインパルス応答ｈsmallを設定する（ステップＳ１１）。
【００９６】
次いで、差異検出部５は、例えばコンサートホール等の所望する音場の広さ感を示す評価値Ｒlargeを設定する（ステップＳ１２）。
【００９７】
次いで、音源補正部７は、分割する周波数帯域数分、ここでは３つの帯域の反射音付加フィルタを設定する（ステップＳ１３）。
【００９８】
次いで、特性判定部４は、再生を行う音場１０のインパルス応答ｈsmallとこれらを合成したフィルタの畳み込み計算を行い、値Ｈを計算する（ステップＳ１４）。
【００９９】
【数５】
【０１００】
次いで、差異検出部５は、特性判定部４で算出された音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとの差異を検出し、値が一致するか否かの判断を行う（ステップＳ１６）。
【０１０１】
差異検出部５が音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとが一致しないと判断した場合には（ステップＳ１６、Ｎｏ）、音源補正部７は、差異検出部５において検出された差異に基づいて、反射音追加フィルタのいずれかを微修正すべく音源成分に付加する情報を算出し、音源補正部７の反射音追加フィルタの修正を行う（ステップＳ１７）。この音源補正部７における付加情報の算出及び反射音追加フィルタの修正方法については前述したとおりである。
【０１０２】
なお、フィルタの修正を行った後は、ステップＳ１４の処理へと戻り、差異検出部５は、Ｒlargeと修正後のＲsmallとの差異の検出を行う。
【０１０３】
一方、ステップＳ１６の判断において、差異検出部５が音場１０の広さ感を示す評価値Ｒsmallと目標とする広さ感を示す評価値Ｒlargeとが一致すると判断した場合には（ステップＳ１６、Ｙｅｓ）、一連の処理を終了する（ステップＳ１８）。
【０１０４】
これにより、一連の音場制御処理の終了後、修正された反射音追加フィルタを用いて例えば音楽の再生が開始されることによって、受聴者は所望する広さ感を体感することが可能となる。
【０１０５】
なお、インパルス応答ｈsmallのデータは、例えば前述した測定方法によりあらかじめ自動車内の座席空間等の対象音場で測定しておき、記録装置２７に記録したデータを用いることが可能である。また所望する音場の広さ感を示す評価値Ｒlargeのデータについても同様に、あらかじめ記録装置２７に記録しておき、これを用いることが可能である。
【０１０６】
以上説明したように、本実施形態によれば、音場再生の対象となる音場としての音場１０のインパルス応答を測定する特性測定部２と、測定されたインパルス応答を、予め設定された周波数帯域毎に複数に分割するとともに、予め設定された経過時間単位毎に複数に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割部３と、分割されたブロックデータ夫々に基づいて、音場１０における空間の音響特性としての対象空間特性を判定する特性判定部４と、判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出部５と、音場１０で受聴に供する音源１を予め設定された周波数帯域毎に音源成分を分割する音源分割部６と、検出された空間特性の差異に基づいて、音源成分毎に少なくとも１の音源成分を補正する音源補正部７と、補正された各音源成分に基づいて音源を再合成する音源合成部８と、を備えた構成を有している。
【０１０７】
したがって、対象音場のインパルス応答を測定し、予め設定された周波数及び経過時間で複数のブロックに分割することによって対象音場の空間特性を判定し、所望する空間特性との差異を検出するとともに、その検出結果に応じて、対象音場で受聴に供する音に対して少なくともいずれか１の前記ブロックの音源を補正する。
【０１０８】
このように、周波数及び経過時間によって異なる反射音パターンを複数のブロックに分けることにより対象音場の空間特性、例えば人間が体感する音場の広さ感について判定することができるため、音場制御や音響ホール設計に活用することが可能となる音場空間特性判定システムを提供することができる。
【０１０９】
また、所望する音場との空間特性の差異に応じて、音源に付加すべき情報を算出し、所定のブロック単位で再生すべき音源の補正を行うので、例えば狭い受聴空間においても広いコンサートホールなど所望する音場で受聴した場合と同様の臨場感を容易に得ることが可能となる音場制御システムを提供することができる。
【０１１０】
また本実施形態によれば、空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、特性判定部４は、夫々のブロックデータのエネルギー値Ｅijと各ブロックデータ毎に定められた重み係数Ｗijとの加重線形和に基づいて空間特性Ｒを判定することを特徴として構成されている。
【０１１１】
したがって、周波数及び時間によりブロックごとに異なる広さ感への影響度を重み付けを行うことによって計算し判定するので、人間が体感する音場の広さ感を客観的な評価値として表すことが可能となる音場空間特性判定システムを提供することができる。
【０１１２】
また、この評価値の差異に基づいてブロック単位で音場制御を行えばよいので、従来必要とされたような厳密なフィルタ制御を行う必要がなく、短時間にフィルタを設計することが可能となり、また分割する帯域数を減らす効果もあり実時間処理を大幅に減らすことが可能となる。
【０１１３】
さらに評価値の差異に基づいて音場制御を行うので、設計者の勘や経験に頼ることなく、客観的な音場制御を行うことが可能な音場制御システムを提供することができる。
【０１１４】
また本実施形態によれば、特性測定部２は、試験信号に基づいて音を拡声するスピーカ２１と、試験信号を発生させるテスト信号発生器２３と、スピーカ２１から発せられる音を収音するマイクロフォン２２と、試験信号と、収音された音とに基づいて、スピーカ２１とマイクロフォン２２との間のインパルス応答を算出するアナライザ２６と、を含む構成を有している。
【０１１５】
したがって、音場空間特性判定システムＳＩにおいてインパルス応答を算出することにより、対象音場の空間特性を容易に測定することが可能な音場空間特性判定システム及びこれを利用した音場制御システムを提供することができる。
【０１１６】
また本実施形態によれば、予め設定された周波数帯域に、予め定められた低中域を含む帯域と、当該低中域を超過する帯域とを少なくとも含む場合に、音源補正部７は、当該低中域を超過する帯域に含まれる音源成分を補正する構成を有している。
【０１１７】
したがって、人間の聴覚のメカニズム上、音場の広さに対する感覚に対して重要な要素である中高域での反射音成分を補正するので、所望する広さ感を効果的に得ることが可能な音場制御システムを提供することができる。
【０１１８】
また本実施形態によれば、周波数軸において分割するために予め設定された周波数としての低中域の値は、５００Ｈｚ〜２ｋＨｚとする構成を有している。
【０１１９】
したがって、低中域の上限として一般的に認識されている値である２ｋＨｚを超える周波数帯域で反射音成分を補正するので、所望する広さ感をより効果的に得ることが可能な音場制御システムを提供することができる。
【０１２０】
また本実施形態によれば、音源補正部７は、直接音が到来する時刻を基準として、当該時刻から予め定められた時間を経過した経過時間に含まれる反射音成分を補正することを特徴とする構成を有している。
【０１２１】
したがって、人間の聴覚のメカニズム上、音場の広さに対する感覚に対して重要な要素である所定時刻（例えば８０ｍｓｅｃ）を超える経過時間で音源に反射音を付加するので、所望する広さ感を効果的に得ることが可能な音場制御システムを提供することができる。
【０１２２】
なお、この所定時刻は例えばクラシックやロック等再生する音楽のジャンルにより任意に変更することができるので、ジャンル等に応じてさらに効果的な音場制御を行うことが可能となる。
【０１２３】
また、本実施形態によれば、音源補正部７は、検出された空間特性の差異に基づいて、少なくともいずれか１の音源成分に付加する付加情報を算出し、当該音源成分に対して算出された付加情報を付加することによって補正を行うことを特徴とする構成を有している。
【０１２４】
したがって、補正方法としては付加情報を算出し音源成分に付加するだけでよいので、従来必要とされたような厳密なフィルタ制御を行う必要がなく、短時間にフィルタを設計することが可能となり、実時間処理を更に大幅に減らすことが可能な音場制御システムを提供することができる。
【０１２５】
また、本実施形態では、音場空間特性判定システムＳＩを対象音場において所望する音場の広さ感を得ることができるように音源に情報を付加し調整する音場制御システムＳＦに適用する構成をとっているが、本判定システムの適用はこれに限るものではない。すなわち、例えば音響ホールやリスニングルーム等再生空間の設計の分野に対しても当該音場空間特性判定システムは適用することができる。またその他空間印象の尺度として本広さ感判定を利用して各種設計及び制御を行うあらゆる分野に対して本判定システムは適用可能である。
【０１２６】
また、本実施形態では、上述のシステム構成によって音場空間特性の判定及び音場の制御を行うようになっているが、システム中にコンピュータおよび記録媒体を備え、この記録媒体に上述の広さ感判定動作及び音場制御動作を行うプログラムを格納し、コンピュータによって当該広さ感判定プログラム及び音場制御プログラムを読み込むことによって上述と同様の動作を行うようにすることも可能である。
【０１２７】
また、この音場空間特性判定プログラム及び音場制御プログラムを実行するシステムにおいて、記録媒体をＤＶＤやＣＤなどの記録媒体により構成するようにしてもよい。
【０１２８】
この場合、当該システムは、記録媒体からプログラムを読み出す読出装置を備えることとなる。
【０１２９】
また、本実施形態では、一つのスピーカに対して一箇所の受音点がある場合を想定しているが、これに限るものではなく、２スピーカを備えたステレオシステムやマルチチャンネルシステムに対しても適用可能である。例えば複数のスピーカを所定の位置において再生する場合には、それぞれのスピーカについて個別にインパルス応答の測定を行い、図８におけるステップＳ１３乃至１７で説明した処理を音源の数だけ繰り返し行えばよい。
【０１３０】
なお、本実施の形態においては、時間軸方向に４分割、周波数軸方向に３分割の計１２ブロックに分割しているが、これに限るものではなく、時間−周波数軸における分割方法及び分割数は任意である。したがって、例えば再生する音楽のジャンル、使用するフィルタ等さまざまな条件を考慮して適宜時間軸方向の分割数及び周波数軸方向の分割数や分割の基準とする時間・周波数を決定することが可能である。
【図面の簡単な説明】
【図１】本願に係る音場空間特性判定システム及び音場制御システムの概要構成の一例を示す機能ブロック図である。
【図２】本願に係る音場空間特性判定システム２の細部構成の一例を示す図である。
【図３】広さ感を示す評価値の一例を示す表である。
【図４】インパルス応答から算出される残響特性と周波数帯域との関係の一例を示した図であり、このうち（Ａ）は低い周波数帯域における残響特性の一例を示した図であり、（Ｂ）は高い周波数帯域における残響特性の一例を示した図である。
【図５】時間周波数分割の一例を示した図である。
【図６】実験結果に基づくソース１についての重み係数の一例を示したグラフである。
【図７】実験結果に基づくソース２についての重み係数の一例を示したグラフである。
【図８】実施形態に係る音場制御システムＳＦによって行われる音場制御処理の動作を示すフローチャートである。
【符号の説明】
ＳＦ … 音場制御システム
ＳＩ … 音場空間特性判定システム
１ … 音源
２ … 特性測定部
３ … 特性分割部
４ … 特性判定部
５ … 差異検出部
６ … 音源分割部
７ … 音源補正部
８ … 音源合成部
１０ … 音場
２１ … スピーカ
２２ … マイクロフォン
２３ … テスト信号発生器
２４ … 増幅器
２５ … 増幅器
２６ … アナライザ
２７ … 記録装置
Claims (13)
- 音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、
前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、
前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、
前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段と、
前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段と、
前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段と、
前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段と、
を備え、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定手段は、
夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴とする音場制御システム。 - 請求項１に記載の音場制御システムにおいて、
前記特性測定手段は、
試験信号に基づいて音を拡声する拡声手段と、
前記試験信号を発生させる信号発生手段と、
前記拡声手段から発せられる前記音を収音する収音手段と、
前記試験信号と、前記収音された前記音とに基づいて、前記拡声手段と前記収音手段との間のインパルス応答を算出する応答算出手段と、
を有することを特徴とする音場制御システム。 - 請求項１又は２に記載の音場制御システムにおいて、
前記補正手段は、前記予め設定された周波数帯域のうち、予め定められた周波数を超過する帯域に含まれる前記音源成分を補正することを特徴とする音場制御システム。 - 請求項１乃至３いずれか１項に記載の音場制御システムにおいて、
前記補正手段は、前記予め設定された前記経過時間のうち、予め定められた経過時刻を超過する前記経過時間に含まれる前記音源成分を補正することを特徴とする音場制御システム。 - 請求項１乃至４いずれか１項に記載の音場制御システムにおいて、
前記補正手段は、
前記検出された前記空間特性の差異に基づいて、少なくともいずれか１の前記音源成分に付加する付加情報を算出する付加情報算出手段と、
前記音源成分に対して当該算出された付加情報を付加することによって補正を行う情報付加手段と、
を有することを特徴とする音場制御システム。 - 特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、
特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、
空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、
差異検出手段が、前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出工程と、
音源分割手段が、前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割工程と、
補正手段が、前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正工程と、
音源合成手段が、前記補正された前記音源成分に基づいて前記音源を再合成する音源合成工程と、
を備え、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定工程においては、
前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴とする音場制御方法。 - コンピュータによって、音場制御を行う音場制御用プログラムであって、
前記コンピュータを、
音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、
前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、
前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、
前記判定された対象空間特性と、所望する音場の空間の音響特性として予め設定された所望空間特性と、の差異を検出する差異検出手段、
前記対象音場で受聴に供する音源を前記予め設定された周波数帯域毎に音源成分を分割する音源分割手段、
前記検出された前記空間特性の差異に基づいて、前記音源成分毎に少なくとも１の前記音源成分を補正する補正手段、
前記補正された前記音源成分に基づいて前記音源を再合成する音源合成手段、
として機能させ、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定手段としての前記コンピュータを、
夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴とする音場制御用プログラム。 - 請求項７に記載の音場制御用プログラムをコンピュータに読み取り可能に記録したことを特徴とする記録媒体。
- 音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段と、
前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段と、
前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段と、
を備え、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定手段は、
夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴とする音場空間特性判定システム。 - 請求項９に記載の音場空間特性判定システムにおいて、
前記特性測定手段は、
試験信号に基づいて音を拡声する拡声手段と、
前記試験信号を発生させる信号発生手段と、
前記拡声手段から発せられる前記音を収音する収音手段と、
前記試験信号と、前記収音された前記音とに基づいて、前記拡声手段と前記収音手段との間のインパルス応答を算出する応答算出手段と、
を有することを特徴とする音場空間特性判定システム。 - 特性測定手段が、音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定工程と、
特性分割手段が、前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割工程と、
空間特性判定手段が、前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定工程と、
を備え、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定工程においては、
前記空間特性判定手段が、夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定することを特徴とする音場空間特性判定方法。 - コンピュータによって、音場空間の特性判定を行う音場空間特性判定用プログラムであって、
前記コンピュータを、
音場再生の対象となる音場としての対象音場のインパルス応答を測定する特性測定手段、
前記測定されたインパルス応答を、予め設定された同一または異なる帯域幅で複数の周波数帯域に分割するとともに、予め設定された同一または異なる時間幅で複数の経過時間単位に分割することによって、時間及び周波数軸上で複数のブロックデータに分割する特性分割手段、
前記分割されたブロックデータ夫々に基づいて、前記対象音場における空間の音響特性としての対象空間特性を判定する空間特性判定手段、
として機能させ、
前記空間の音響特性が、人間が体感する音場の広狭に対する感覚としての広さ感を示す数値特性である場合に、
前記空間特性判定手段としての前記コンピュータを、
夫々の前記ブロックデータのエネルギー値と各ブロックデータ毎に定められた重み係数との加重線形和に基づいて前記数値特性を判定するように機能させることを特徴とする音場空間特性判定用プログラム。 - 請求項１２に記載の音場空間特性判定用プログラムをコンピュータに読み取り可能に記録したことを特徴とする記録媒体。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003067814A JP4130779B2 (ja) | 2003-03-13 | 2003-03-13 | 音場制御システム及び音場制御方法 |
| US10/798,944 US20040179696A1 (en) | 2003-03-13 | 2004-03-12 | 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 |
| EP04006014A EP1458218A3 (en) | 2003-03-13 | 2004-03-12 | Sound field control system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003067814A JP4130779B2 (ja) | 2003-03-13 | 2003-03-13 | 音場制御システム及び音場制御方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004279525A JP2004279525A (ja) | 2004-10-07 |
| JP4130779B2 true JP4130779B2 (ja) | 2008-08-06 |
Family
ID=32767955
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003067814A Expired - Fee Related JP4130779B2 (ja) | 2003-03-13 | 2003-03-13 | 音場制御システム及び音場制御方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20040179696A1 (ja) |
| EP (1) | EP1458218A3 (ja) |
| JP (1) | JP4130779B2 (ja) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7864631B2 (en) * | 2005-06-09 | 2011-01-04 | Koninklijke Philips Electronics N.V. | Method of and system for determining distances between loudspeakers |
| EP1740016B1 (en) * | 2005-06-28 | 2010-02-24 | AKG Acoustics GmbH | Method for the simulation of a room impression and/or sound impression |
| JP4396646B2 (ja) | 2006-02-07 | 2010-01-13 | ヤマハ株式会社 | 応答波形合成方法、応答波形合成装置、音響設計支援装置および音響設計支援プログラム |
| JP4668118B2 (ja) | 2006-04-28 | 2011-04-13 | ヤマハ株式会社 | 音場制御装置 |
| US9015051B2 (en) * | 2007-03-21 | 2015-04-21 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Reconstruction of audio channels with direction parameters indicating direction of origin |
| JP5286739B2 (ja) * | 2007-10-18 | 2013-09-11 | ヤマハ株式会社 | 音像定位パラメータ算出装置、音像定位制御装置、音像定位装置及びプログラム |
| TWI475896B (zh) * | 2008-09-25 | 2015-03-01 | Dolby Lab Licensing Corp | 單音相容性及揚聲器相容性之立體聲濾波器 |
| WO2011015932A1 (en) | 2009-08-03 | 2011-02-10 | Imax Corporation | Systems and method for monitoring cinema loudspeakers and compensating for quality problems |
| JP4893789B2 (ja) * | 2009-08-10 | 2012-03-07 | ヤマハ株式会社 | 音場制御装置 |
| JP5658588B2 (ja) * | 2011-02-07 | 2015-01-28 | 日本放送協会 | 聴覚臨場感評価装置及び聴覚臨場感評価プログラム |
| KR20190106902A (ko) * | 2019-08-29 | 2019-09-18 | 엘지전자 주식회사 | 사운드 분석 방법 및 장치 |
| WO2022220036A1 (ja) * | 2021-04-12 | 2022-10-20 | パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカ | 音響特徴量推定方法、音響特徴量推定システム、プログラム、及び、レンダリング方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL8800745A (nl) * | 1988-03-24 | 1989-10-16 | Augustinus Johannes Berkhout | Werkwijze en inrichting voor het creeren van een variabele akoestiek in een ruimte. |
| US5146507A (en) * | 1989-02-23 | 1992-09-08 | Yamaha Corporation | Audio reproduction characteristics control device |
| US5572443A (en) * | 1993-05-11 | 1996-11-05 | Yamaha Corporation | Acoustic characteristic correction device |
| JP2758846B2 (ja) * | 1995-02-27 | 1998-05-28 | 埼玉日本電気株式会社 | ノイズキャンセラ装置 |
| JP2956642B2 (ja) * | 1996-06-17 | 1999-10-04 | ヤマハ株式会社 | 音場制御ユニットおよび音場制御装置 |
| CA2382786A1 (en) * | 2001-04-30 | 2002-10-30 | The Government Of The United States Of America, As Represented By The Se Cretary, Department Of Health And Human Services, Centers For Disease Co | Auscultatory training system |
-
2003
- 2003-03-13 JP JP2003067814A patent/JP4130779B2/ja not_active Expired - Fee Related
-
2004
- 2004-03-12 US US10/798,944 patent/US20040179696A1/en not_active Abandoned
- 2004-03-12 EP EP04006014A patent/EP1458218A3/en not_active Withdrawn
Also Published As
| Publication number | Publication date |
|---|---|
| US20040179696A1 (en) | 2004-09-16 |
| JP2004279525A (ja) | 2004-10-07 |
| EP1458218A3 (en) | 2008-05-07 |
| EP1458218A2 (en) | 2004-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10104485B2 (en) | Headphone response measurement and equalization | |
| JP5666023B2 (ja) | 残響知覚レベルの大きさを決定する装置及び方法、オーディオプロセッサ並びに信号処理方法 | |
| US9264834B2 (en) | System for modifying an acoustic space with audio source content | |
| JP4059478B2 (ja) | 音場制御方法及び音場制御システム | |
| Larsen et al. | On the minimum audible difference in direct-to-reverberant energy ratio | |
| EP1571884A2 (en) | Sound reproducing method and apparatus | |
| KR101989062B1 (ko) | 오디오 신호를 향상시키기 위한 장치 및 방법 및 음향 향상 시스템 | |
| US20070025559A1 (en) | Audio tuning system | |
| JP4130779B2 (ja) | 音場制御システム及び音場制御方法 | |
| EP2949133B1 (en) | Automatic loudspeaker polarity detection | |
| JP4904461B2 (ja) | 音声周波数応答処理システム | |
| Lee et al. | Comparison of psychoacoustic-based reverberance parameters | |
| Lavandier et al. | Identification of some perceptual dimensions underlying loudspeaker dissimilarities | |
| JP7010649B2 (ja) | オーディオ信号処理装置及びオーディオ信号処理方法 | |
| JP6998823B2 (ja) | マルチチャンネル客観評価装置及びプログラム | |
| US20230199419A1 (en) | System, apparatus, and method for multi-dimensional adaptive microphone-loudspeaker array sets for room correction and equalization | |
| JP2009027331A (ja) | 音場再現システム | |
| JP2023012347A (ja) | 音響装置および音響制御方法 | |
| JP2009027332A (ja) | 音場再現システム | |
| US20180286427A1 (en) | System and a method for determining an interference or distraction | |
| Zhu et al. | Relative Contribution of Frequency and Parameter Values to Selectivity for Interaural Correlation | |
| KR20090065788A (ko) | 청력 측정 장치 및 그 제어방법 | |
| Soeta et al. | Application to Sound Design | |
| Ovans | Equalization in Sound Reinforcement: Psychoacoustics, Methods, and Issues | |
| JP2006279242A (ja) | ステレオ音声用ラウドネスレベルメータ |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060223 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080212 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080403 |
|
| 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: 20080520 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080523 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110530 Year of fee payment: 3 |
|
| LAPS | Cancellation because of no payment of annual fees |