(57)【要約】【課題】  低周波数領域におけるパワーが相対的に大となる周波数特性を有したノイズが存在する場合であっても、高精度なインパルス応答測定を可能とする。【解決手段】  高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる評価音信号を生成し、スピーカ２に入力して音響空間Ｒのインパルス応答を測定するようにした。【選択図】      図３
2005338521.tif 000002
【特許請求の範囲】
【請求項１】
音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定装置において、
高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる前記評価音信号を生成し、前記発音点に入力する
ことを特徴とするインパルス応答測定装置。
【請求項２】
音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定装置において、
高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、時間領域における信号波形のエンベロープが略一定値となる前記評価音信号を生成し、前記発音点に入力する
ことを特徴とするインパルス応答測定装置。
【請求項３】
前記評価音信号は、位相の周波数特性を２階微分した値が、振幅の周波数特性の絶対値を２乗した値に比例する信号であることを特徴とする請求項１または２に記載のインパルス応答測定装置。
【請求項４】
前記評価音信号は、前記振幅の値が、周波数の逆数の平方根の値に比例する周波数特性を有することを特徴とする請求項３に記載のインパルス応答測定装置。
【請求項５】
音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定方法において、
高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる前記評価音信号を生成し、前記発音点に入力する
ことを特徴とするインパルス応答測定方法。
【請求項６】
音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定方法において、
高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、時間領域における信号波形のエンベロープが略一定値となる前記評価音信号を生成し、前記発音点に入力する
ことを特徴とするインパルス応答測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コンサートホールや劇場などの音響空間におけるインパルス応答を測定するための技術に関する。
【背景技術】
【０００２】
近年、ユーザが家庭などの部屋において音楽を聴く場合に、あたかもコンサートホールや劇場などで生演奏を聴いているかのような臨場感を味わえるようにする技術として、コンサートホールなどの音響空間の音場を別の音響空間に再現するための音場再現技術が知られている。この音場再現技術におていは、評価音を音響室内に放ったときのインパルス応答を音響特性として測定するといったことが行われている。また、上記評価音としては、例えば図５及び図６に示すＭ系列信号や図７及び図８に示す時間伸張パルス（ＴＳＰ：Time Stretched Pulse）信号、図９及び図１０に示すピンクノイズ等が用いられている（例えば、特許文献１参照）。
【特許文献１】特開２０００−６９５９７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、Ｍ系列信号及び時間伸張パルス信号は、図６及び図８に示すように、周波数領域における信号レベルが略一定であるため、これらＭ系列信号或いは時間伸張パルス信号を評価音として用いた場合、低周波数領域のパワーが相対的に大きくなるといった周波数特性を有するノイズが音響空間に存在したときに、低周波数領域におけるインパルス応答の測定精度（Ｓ／Ｎ比）が劣化するといった問題がある。
【０００４】
また、ピンクノイズは、図１０に示すよういに、低周波数領域で比較的大きな信号パワーを有するため、上記のようなノイズが音響空間に存在したとしても、低周波数領域におけるインパルス応答測定精度（Ｓ／Ｎ比）を良好なものに維持できるが、次のような問題がある。すなわち、ピンクノイズは、図９に示すように、信号の平均パワーに対する最大振幅の比が大きい信号であるため、ピンクノイズの平均パワーを大にしてインパルス応答を測定すると、過度な最大振幅が音響空間に与えられてしまい、線形性が崩れ、インパルス応答の波形が歪んでしまうといった問題がある。従って、ピンクノイズの平均パワーを小にしてインパルス応答を測定する必要があるものの、小さな平均パワーでの測定しかできないため、インパルス応答測定装置のダイナミックレンジを有効に活用しての測定ができないといった問題がある。
【０００５】
本発明は、上述した事情に鑑みてなされたものであり、低周波数領域におけるパワーが相対的に大となる周波数特性を有したノイズが存在する場合であっても、高精度にインパルス応答を測定することのできるインパルス応答生成装置及びインパルス応答生成方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記目的を達成するために、本発明は、音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定装置において、高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる前記評価音信号を生成し、前記発音点に入力することを特徴とする。
【０００７】
また上記目的を達成するために、本発明は、音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定装置において、高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、時間領域における信号波形のエンベロープが略一定値となる前記評価音信号を生成し、前記発音点に入力することを特徴とする。
【０００８】
また本発明は、上記発明において、前記評価音信号は、位相の周波数特性を２階微分した値が、振幅の周波数特性の絶対値を２乗した値に比例する信号であることを特徴とする。
【０００９】
また本発明は、上記発明において、前記評価音信号は、前記振幅の値が、周波数の逆数の平方根の値に比例する周波数特性を有することを特徴とする。
【００１０】
また上記目的を達成するために、本発明は、音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定方法において、高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる前記評価音信号を生成し、前記発音点に入力することを特徴とする。
【００１１】
また上記目的を達成するために、本発明は、音響空間に配置された発音点に評価音信号を入力して評価音を発せさせ、受音点にて測定された測定信号と前記評価音信号とに基づいて前記発音点と前記受音点間のインパルス応答を測定するインパルス応答測定方法において、高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、時間領域における信号波形のエンベロープが略一定値となる前記評価音信号を生成し、前記発音点に入力することを特徴とする。
【発明の効果】
【００１２】
本発明によれば、低周波数領域におけるパワーが相対的に大となる周波数特性を有したノイズが存在する場合であっても、高精度にインパルス応答を測定することができる。
【発明を実施するための最良の形態】
【００１３】
以下、図面を参照して本発明の実施の形態について説明する。図１は、本実施の形態に係るインパルス応答測定システムＳの構成を模式的に示す図である。この図に示すように、インパルス応答測定システムＳは、インパルス応答測定装置１と、音響特性の測定対象となる音響空間Ｒに配置された放音装置としてのスピーカ２及び収音装置としてのマイク３とを備えている。
【００１４】
インパルス応答測定装置１は、インパルス応答の測定用信号としての評価信号４を生成し、上記スピーカ２に出力して評価音を発せさせると共に、当該評価音をマイク３が収音した際の収音信号５を受け取り、評価信号４及び収音信号５に基づいて、スピーカ２とマイク３間、すなわち、発音点と受音点間のインパルス応答を算出するものであり、上記評価信号４を生成する音源装置や、インパルス応答算出演算を実行する演算装置などを備えて構成されている。
【００１５】
なお、スピーカ２（発音点）やマイク３（受音点）の数は適宜に変更可能である。また、マイク３の代わりに、人体の頭部を模したダミーヘッドを用いても良い。このように発音点と受音点との数を増やした場合には、インパルス応答測定装置１は、発音点の各々と受音点の各々との間のインパルス応答をそれぞれ算出する。
【００１６】
インパルス応答の生成について詳細には、評価信号４のフーリエ変換をＸ（ω）、音響空間Ｒにおけるスピーカ２とマイク３間のインパルス応答のフーリエ変換をＨ（ω）、マイク３にて収音された収音信号５のフーリエ変換をＹ（ω）とすると、次式（１）の関係が成り立つ。
【００１７】
【数１】
2005338521.tif 000003
【００１８】
そこで、インパルス応答測定装置１は、評価信号４及び収音信号５に基づいて、次式（１）に示すＨ（ω）（いわゆる伝達周波数特性）を求め、このＨ（ω）を逆フーリエ変換することで上記インパルス応答を算出することとしている。
【００１９】
ところで、図２に示すように、音響空間Ｒにノイズが存在する場合、そのノイズ信号のフーリエ変換をＮ（ω）とすると、マイク３で測定される信号は、Ｙ（ω）＋Ｎ（ω）となり、そのときのインパルス応答の測定値Ｈ’（ω）は次式（２）のようになる。
【００２０】
【数２】
2005338521.tif 000004
【００２１】
したがって、音響空間Ｒにノイズが存在する場合には、そのノイズ信号により、測定されたインパルス応答に誤差が含まれる。特に、ノイズ信号が周波数領域において一定の信号レベルではなく、低周波数領域において比較的信号レベルが大きくなるといった周波数特性を有している場合には、低周波数領域において、上記（２）式の第２項の値が大きくなるため、当該低周波数領域における誤差が増大する。
【００２２】
そこで、本実施の形態では、低周波数領域における誤差を抑えるために、上記評価信号４として次のような信号を用いることとしている。すなわち、高周波数領域のパワーよりも低周波数領域のパワーが大きく、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内（すなわち、突出した最大振幅を有しない）となる信号である。
【００２３】
詳述すると、評価信号４として、高周波数領域のパワーよりも低周波数領域のパワーが大きい信号を用いることで、低周波数領域における上記（２）式第２項の値を小さくし、当該低周波数領域での誤差を抑え、結果として、インパルス応答の測定精度を向上させることができるのである。
【００２４】
さらに、評価信号４として、評価信号４の平均パワーに対する評価信号４の最大振幅の比が所定範囲内となる信号を用いることで、インパルス応答測定装置１が有する測定可能範囲（ダイナミックレンジ）を最大限に生かした測定が可能となる。すなわち、平均パワーに対する評価信号４の最大振幅の比が所定範囲内であるため、評価信号４の平均パワーを大としても、評価信号４の最大振幅が突出することがなく、これにより、当該評価信号４をスピーカ２に入力して評価音を発せさせた場合であっても、音響空間Ｒの応答が線形性を維持する範囲内でインパルス応答を測定できる。この結果、インパルス応答測定装置１が有する測定可能範囲（ダイナミックレンジ）を有効に活用してのインパルス応答測定が可能となる。
【００２５】
以下、上記のような特性を有する評価信号４について詳述する。
【００２６】
周波数領域における振幅特性をＡ（ω）、位相特性をＰ（ω）とした場合、上記評価信号４は次式（３）を満たせば良い。
【００２７】
【数３】
2005338521.tif 000005
【００２８】
すなわち、この評価信号４に対して逆フーリエ変換処理を施し、時間波形を求めると、振幅特性がＡ（ω）であり、その信号波形が、時間領域において略一定となるエンベロープを有する信号となる。
【００２９】
次いで、このような信号波形を有するＮ個のデジタル信号列ｘｉ（ｉ＝０・・・Ｎ−１）を生成する場合について説明する。サンプリング周波数をｆｓとした場合、周波数ωは式（４）のように離散化され、また、上記振幅特性Ａ（ω）及び位相特性Ｐ（ω）はそれぞれ式（５）、式（６）にて表される。
【００３０】
【数４】
2005338521.tif 000006
【００３１】
【数５】
2005338521.tif 000007
【００３２】
【数６】
2005338521.tif 000008
【００３３】
また、デジタル信号列ｘｉの逆フーリエ変換であるＸｋは、次式（７）、式（８）のように示され、そして、式（９）に示すように、これらのＸｋを離散フーリエ逆変換することでデジタル信号列ｘｉを得ることができる。したがって、振幅特性Ａｋ及び位相特性Ｐｋが上記式（３）を満足するように与えられることで、上記のような信号波形を有するデジタル信号ｘｉを得ることができる。
【００３４】
【数７】
2005338521.tif 000009
【００３５】
【数８】
2005338521.tif 000010
【００３６】
【数９】
2005338521.tif 000011
【００３７】
ここで、振幅特性Ａ（ω）が周波数によらず一定である場合には、Ｐ（ω）∝ω2となり、従来のＴＳＰ信号と同じになる。これに対して、本実施の形態では、振幅特性Ａ（ω）に周波数特性を持たせた信号を生成することとしている。例えば、データ数がＮであり、信号パワーが１／ｆ特性（周波数に反比例する特性）を有する信号は、周波数領域での振幅Ａｋと位相特性Ｐｋとを次式（１０）乃至（１４）ように設定し、式（１５）を用いて離散フーリエ逆変換することで生成することができる。但し、ｍは任意の整数である。
【００３８】
【数１０】
2005338521.tif 000012
【００３９】
【数１１】
2005338521.tif 000013
【００４０】
【数１２】
2005338521.tif 000014
【００４１】
【数１３】
2005338521.tif 000015
【００４２】
【数１４】
2005338521.tif 000016
【００４３】
【数１５】
2005338521.tif 000017
【００４４】
【数１６】
2005338521.tif 000018
【００４５】
【数１７】
2005338521.tif 000019
【００４６】
例えば、上記式（１０）乃至（１５）において、Ｎ＝２16、ｍ＝２15とすれば、図３に示すように、時間領域においてエンベロープが略一定値となる信号波形であり、なおかつ、図４に示すように、低周波数領域の信号パワーが高周波数領域よりも高いデジタル信号列ｘｉを得ることができ、上記評価信号４として用いるに好適なデジタル信号列ｘｉが得られる。
【００４７】
以上説明したように、本実施の形態によれば、インパルス応答測定装置１が、高周波数領域のパワーよりも低周波数領域のパワーが大きくなる周波数特性を有し、なおかつ、信号の平均パワーに対する信号の最大振幅の比が所定範囲内となる評価音信号を生成し、スピーカ２に入力して音響空間Ｒのインパルス応答を測定するようにしたため、低周波数領域におけるパワーが相対的に大となる周波数特性を有したノイズが音響空間Ｒに存在する場合であっても、当該低周波数領域におけるインパルス応答の測定精度（Ｓ／Ｎ比）を向上させることができ、以って、高精度にインパルス応答を測定することができる。
【００４８】
また、評価信号４として、評価信号４の平均パワーに対する評価信号４の最大振幅の比が所定範囲内となる信号を用いるため、評価信号４の平均パワーを大としても、評価信号４の最大振幅が突出することがなく、結果として、当該評価信号４をスピーカ２に入力して評価音を発せさせた場合であっても、音響空間Ｒの応答が線形性を維持する範囲内でインパルス応答を測定できる。これにより、インパルス応答測定装置１が有する測定可能範囲を有効に活用してのインパルス応答測定が可能となる。
【００４９】
上述した実施の形態は、あくまでも本発明の一態様を示すものであり、本発明の範囲内で任意に変形可能である。例えば、上述した実施の形態では、インパルス応答測定対象の音響空間Ｒにスピーカ２とマイク３とを実際に配置してインパルス応答を測定するようにしたが、これに限らず、例えば、上記音響空間Ｒの構造や壁材等の音響特性を予めデータ化し、インパルス応答測定装置１が当該データ化された音響空間Ｒに対してコンピュータシミュレーションによりインパルス応答を算出する構成としても良い。
【図面の簡単な説明】
【００５０】
【図１】本発明の実施の形態に係るインパルス応答測定システムの構成を示す図である。
【図２】インパルス応答の算出方法を説明するための図である。
【図３】評価音信号の時間領域における振幅特性を示す図である。
【図４】評価音信号の周波数特性を示す図である。
【図５】Ｍ系列信号の時間領域における振幅特性を示す図である。
【図６】Ｍ系列信号の周波数特性を示す図である。
【図７】時間伸張パルス信号の時間領域における振幅特性を示す図である。
【図８】時間伸張パルス信号の周波数特性を示す図である。
【図９】ピンクノイズの時間領域における振幅特性を示す図である。
【図１０】ピンクノイズの周波数特性を示す図である。
【符号の説明】
【００５１】
１  インパルス応答測定装置
２  スピーカ（発音点）
３  マイク（受音点）
４  評価信号
５  収音信号
Ｒ  音響空間
Ｓ  インパルス応答測定システム
【図１】
2005338521.tif 000020
【図２】
2005338521.tif 000021
【図３】
2005338521.tif 000022
【図４】
2005338521.tif 000023
【図５】
2005338521.tif 000024
【図６】
2005338521.tif 000025
【図７】
2005338521.tif 000026
【図８】
2005338521.tif 000027
【図９】
2005338521.tif 000028
【図１０】
2005338521.tif 000029