(57)【要約】【課題】符号長が２ｎ−１の周期関数からなる測定信号を用いて、精度良く音場環境の周波数特性を測定すること。【解決手段】音場測定装置１は、符号長が２ｎ−１（ｎは自然数）の周期関数からなる測定信号をスピーカ９へ出力する外部出力手段６と、スピーカ９より出力された測定信号を収音するマイクロフォン７と、収音された測定音を２ｍ（ｍは自然数）のサンプル長でフーリエ変換して周波数特性を求めるフーリエ変換手段１２と、求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外（ｋ＝０，１,２,・・・）の線スペクトルを間引き処理する間引き手段１３と、間引き処理された周波数特性に基づいて、信号レベルが平均化された音場の周波数特性を求める平均化処理手段１４とを有する。【選択図】図２
2015075345.tif 000002
【特許請求の範囲】
【請求項１】
符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、当該スピーカに対して前記測定信号を出力する外部出力手段と、
該スピーカより出力された前記測定信号を収音するマイクロフォンと、
該マイクロフォンで収音された測定音を２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換手段と、
該フーリエ変換手段によって求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、前記周波数特性におけるノイズを除去する間引き手段と、
該間引き手段により間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された音場の周波数特性を求める平均化処理手段と
を有し、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定装置。
【請求項２】
符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、当該スピーカに対して前記測定信号を出力する外部出力手段と、
該スピーカより出力された前記測定信号を収音するマイクロフォンと、
該マイクロフォンで収音された測定音を２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換手段と、
該フーリエ変換手段により求められた周波数特性を帯域分割することにより、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成する帯域分割手段と、
該帯域分割手段によって生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、前記第２周波数特性におけるノイズを除去する間引き手段と、
前記帯域分割手段により生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された第１周波数特性を生成する第１平均化処理手段と、
前記間引き手段により間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された第２周波数特性を生成する第２平均化処理手段と、
前記第１平均化処理手段により平均化処理された前記第１周波数特性と、前記第２平均化処理手段により平均化処理された前記第２周波数特性とを合成することにより、全帯域の信号成分を備えた音場の周波数特性を求める合成手段と
を有し、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定装置。
【請求項３】
符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、外部出力手段が、当該スピーカに対して前記測定信号を出力する外部出力ステップと、
該外部出力ステップにおいて前記スピーカより出力された前記測定信号を、マイクロフォンで収音する収音ステップと、
該収音ステップにおいて前記マイクロフォンにより収音された測定音を、フーリエ変換手段が、２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換ステップと、
該フーリエ変換ステップにおいて求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、間引き手段が、前記周波数特性におけるノイズを除去する間引きステップと、
該間引きステップにおいて間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理手段が、平均化処理された音場の周波数特性を求める平均化処理ステップと
を有し、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定装置の音場測定方法。
【請求項４】
符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、外部出力手段が、当該スピーカに対して前記測定信号を出力する外部出力ステップと、
該外部出力ステップにおいて前記スピーカより出力された前記測定信号を、マイクロフォンで収音する収音ステップと、
該収音ステップにおいて前記マイクロフォンで収音された測定音を、フーリエ変換手段が、２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換ステップと、
該フーリエ変換ステップにおいて求められた周波数特性を帯域分割することにより、帯域分割手段が、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成する帯域分割ステップと、
該帯域分割ステップにおいて生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、間引き手段が、前記第２周波数特性におけるノイズを除去する間引きステップと、
前記帯域分割ステップにおいて生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、第１平均化処理手段が、平均化処理された第１周波数特性を生成する第１平均化処理ステップと、
前記間引きステップにおいて間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、第２平均化処理手段が、平均化処理された第２周波数特性を生成する第２平均化処理ステップと、
前記第１平均化処理ステップにおいて平均化処理された前記第１周波数特性と、前記第２平均化処理ステップにおいて平均化処理された前記第２周波数特性とを合成することにより、合成手段が、全帯域の信号成分を備えた音場の周波数特性を求める合成ステップと
を有し、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定装置の音場測定方法。
【請求項５】
符号長が２ｎ−１の周期関数からなる測定信号を用いて音場の周波数特性を測定する音場測定装置の音場測定プログラムであって、
該音場測定装置のコンピュータに、
スピーカより前記測定信号を出力させるために当該スピーカに対して前記測定信号を出力させる外部出力機能と、
該外部出力機能により前記スピーカより出力された前記測定信号を、マイクロフォンを用いて収音させる収音機能と、
該収音機能によって収音された測定音を、２ｍのサンプル長でフーリエ変換させることにより、周波数特性を求めさせるフーリエ変換機能と、
該フーリエ変換機能によって求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理させることにより、前記周波数特性におけるノイズを除去させる間引き機能と、
該間引き機能により間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された前記音場の周波数特性を求めさせる平均化処理機能と
を実現させるためのプログラムであって、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定プログラム。
【請求項６】
符号長が２ｎ−１の周期関数からなる測定信号を用いて音場の周波数特性を測定する音場測定装置の音場測定プログラムであって、
該音場測定装置のコンピュータに、
スピーカより前記測定信号を出力させるために当該スピーカに対して前記測定信号を出力させる外部出力機能と、
該外部出力機能により前記スピーカより出力された前記測定信号を、マイクロフォンを用いて収音させる収音機能と、
該収音機能によって収音された測定音を、２ｍのサンプル長でフーリエ変換させることにより、周波数特性を求めさせるフーリエ変換機能と、
該フーリエ変換機能によって求められた周波数特性を帯域分割することにより、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成させる帯域分割機能と、
該帯域分割機能によって生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理させることにより、前記第２周波数特性におけるノイズを除去させる間引き機能と、
前記帯域分割機能により生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された第１周波数特性を生成させる第１平均化処理機能と、
前記間引き機能により間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された第２周波数特性を生成させる第２平均化処理機能と、
前記第１平均化処理機能により平均化処理された前記第１周波数特性と、前記第２平均化処理機能により平均化処理された前記第２周波数特性とを合成することにより、全帯域の信号成分を備えた前記音場の周波数特性を求めさせる合成機能と
を実現させるためのプログラムであって、
前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であること
を特徴とする音場測定装置の音場測定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は音場測定装置、音場測定方法および音場測定プログラムに関し、より詳細には、符号長が２ｎ−１（ｎは自然数）の周期関数からなる測定信号を用いて、音場環境の周波数特性を精度良く測定することが可能な音場測定装置、音場測定方法および音場測定プログラムに関する。
【背景技術】
【０００２】
従来より、スピーカ等が設置される音場環境において周波数特性を測定し、測定された周波数特性に基づいてオーディオ機器のイコライザを調節したり、出力音を音場に合わせて予め補正することにより、音場環境に最適な音質の音楽を提供する方法が知られている。
【０００３】
周波数特性を測定するための測定信号として、ＰＮ（Pseudorandom Noise：疑似ランダム雑音）符号やＴＳＰ信号（Time Stretched Pulse）が知られている。ＰＮ符号は一般的にランダムな雑音から成る人工的な測定信号であり、Ｍ系列符号（Maximum Length Sequence）やＧｏｌｄ系列符号などが該当する。
【０００４】
Ｍ系列符号やＧｏｌｄ系列符号は、所定の長さのシフトレジスタと排他的論理和を用いてフィードバックすることで生成することができる。シフトレジスタの長さ（段数）がｎの場合（ｎは自然数）、符号の周期（符号長）は２ｎ−１となり、シフトレジスタのフィードバック位置は生成多項式で求められる。Ｍ系列符号の場合には、出力信号が０と１の２値で直流成分が多く含まれた信号になるため、０を−１に変換して出力される。このように音場の周波数特性の測定に用いられる測定信号として、符号長が２ｎ−１の周期関数が用いられる。
【０００５】
これらの測定信号を用いて周波数特性を測定する場合には、例えば、スピーカより出力された測定信号を、聴取位置に設置されるマイクロフォンで収音したのちに、収音された信号をフーリエ変換することにより、音場環境の周波数特性を求める方法を用いる（例えば、特許文献１、特許文献２参照）。また、出力される測定信号をリファレンスとして用い、マイクロフォンで収音された測定信号との相互相関特性を求めることによって、インパルス応答を求めることも可能である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平７−７５１９０号公報
【特許文献２】特開２００７−２３２４９２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述したように、符号長が２ｎ−１の周期関数からなる測定信号を用いて音場の周波数特性を求める場合には、マイクロフォンで収音した測定音をフーリエ変換処理する必要がある。フーリエ変換処理を行う場合には、フーリエ変換のサンプル長を測定信号の符号長に対して２倍以上に設定することが多く、２倍以上に設定することにより、フーリエ変換毎の振幅スペクトルの変動が抑えられ、ほぼ一定の周波数特性を得ることが可能になる。
【０００８】
しかしながら、フーリエ変換のサンプル長は、一般的に２ｍ（ｍは自然数、但し、ｍ＞ｎ）となる一方で、測定信号の符号長は、２ｎ−１となる。このため、測定信号を用いて周波数特性を求める場合には、測定信号の符号長に対して、フーリエ変換のサンプル長が整数倍の関係にならず、非同期の関係になってしまう傾向がある。このようにフーリエ変換のサンプル長と、測定信号の符号長とが非同期の関係になると、求められた一様な線スペクトルの間に変動するレベルの小さなスペクトルが発生してしまい、ノイズとして検出されてしまうという問題があった。
【０００９】
しかしながら、符号長が２ｎ−１の周期関数からなる測定信号を用いる場合であっても、長い符号長の測定信号を用いることにより、振幅スペクトルの変動を低減させることができる。例えば、図１２（ａ）は、３２，７６７の符号長を持つＭ系列符号を用いた場合の周波数特性を示しており、図１２（ｂ）は、１／３オクターブバンド幅で対数的な平均化処理を行った場合の周波数特性を示している。長い符号長のＭ系列符号を用いた場合、振幅スペクトル毎にレベルの差異が発生してしまうが、対数的な平均化処理を行うことによって、周波数特性をほぼ一様な状態にすることができる。
【００１０】
このように、測定信号の符号長が長い場合には、フーリエ変換によって振幅スペクトル数が多くなり、周波数間隔が密になるため、平均化処理によりノイズを低減させることが可能となる。しかしながら、測定信号の符号長が長い場合には、フーリエ変換等に必要とされるメモリ量等が増加し、必要とされる処理時間や処理負担も増大してしまうという問題がある。
【００１１】
一方で、短い符号長の測定信号を用いることにより、フーリエ変換処理に必要とされるメモリ量を減少させ、また処理時間および処理負担を低減させることができる。図１３は、短い符号長の測定信号を用いて求められた周波数特性を示している。図１３（ａ）は、４,０９６の符号長を持つＭ系列符号を用いた場合の周波数特性を示しており、図１３（ｂ）は、１／３オクターブバンド幅で対数的な平均化処理を行った場合の周波数特性を示している。符号長を短くすると、測定時間や測定負担等を低減させ、使用するメモリ量を少なくすることができるが、図１３（ａ）に示すように、振幅スペクトルの周波数が粗くなってしまうという問題があった。
【００１２】
また、対数的な平均化処理を行った場合であっても、周波数特性に対して信号レベルの変動が生じてしまうという問題があった。図１３（ｂ）は、人間の聴覚特性に合わせて、１／３オクターブバンド幅で対数的な平均化処理を行った場合を示しているが、低中域において信号レベルが大きく変動する様子が示されている。
【００１３】
このように、符号長が２ｎ−１の周期関数からなる測定信号を用いて周波数特性を測定する場合には、フーリエ変換処理によって変動する小さな変動（ノイズ）が発生してしまい、精度良く音場の周波数特性を測定することが容易ではないという問題があった。
【００１４】
本発明は上記問題に鑑みてなされたものであり、符号長が２ｎ−１の周期関数からなる測定信号を用いて、精度良く音場環境の周波数特性を測定することが可能な音場測定装置、音場測定方法および音場測定プログラムを提供することを課題とする。
【課題を解決するための手段】
【００１５】
上記課題を解決するために、本発明に係る音場測定装置は、符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、当該スピーカに対して前記測定信号を出力する外部出力手段と、該スピーカより出力された前記測定信号を収音するマイクロフォンと、該マイクロフォンで収音された測定音を２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換手段と、該フーリエ変換手段によって求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、前記周波数特性におけるノイズを除去する間引き手段と、該間引き手段により間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された音場の周波数特性を求める平均化処理手段とを有し、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００１６】
また、本発明に係る音場測定装置の音場測定方法は、符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、外部出力手段が、当該スピーカに対して前記測定信号を出力する外部出力ステップと、該外部出力ステップにおいて前記スピーカより出力された前記測定信号を、マイクロフォンで収音する収音ステップと、該収音ステップにおいて前記マイクロフォンにより収音された測定音を、フーリエ変換手段が、２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換ステップと、該フーリエ変換ステップにおいて求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、間引き手段が、前記周波数特性におけるノイズを除去する間引きステップと、該間引きステップにおいて間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理手段が、平均化処理された音場の周波数特性を求める平均化処理ステップとを有し、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００１７】
さらに、本発明に係る音場測定装置の音場測定プログラムは、符号長が２ｎ−１の周期関数からなる測定信号を用いて音場の周波数特性を測定する音場測定装置の音場測定プログラムであって、該音場測定装置のコンピュータに、スピーカより前記測定信号を出力させるために当該スピーカに対して前記測定信号を出力させる外部出力機能と、該外部出力機能により前記スピーカより出力された前記測定信号を、マイクロフォンを用いて収音させる収音機能と、該収音機能によって収音された測定音を、２ｍのサンプル長でフーリエ変換させることにより、周波数特性を求めさせるフーリエ変換機能と、該フーリエ変換機能によって求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理させることにより、前記周波数特性におけるノイズを除去させる間引き機能と、該間引き機能により間引き処理された周波数特性に基づいて、所定の周波数間隔における信号レベルの平均値を、当該周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された前記音場の周波数特性を求めさせる平均化処理機能とを実現させるための音場測定プログラムであって、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００１８】
符号長が２ｎ−１の周期関数からなる測定信号をスピーカ等から出力させて、収音された信号をサンプル数２ｍでフーリエ変換させた場合には、フーリエ変換長（サンプル長）が測定信号の符号長の整数倍の関係にはならない。このようにフーリエ変換長と測定信号の符号長とが整数倍の関係とはならず、非同期となってしまう場合には、フーリエ変換毎に、変動するレベルの小さな線スペクトルが、一様な線スペクトルの間に発生してしまうおそれがある。この変動するレベルの小さな線スペクトルは、検出された周波数特性におけるノイズになってしまう。
【００１９】
このため、本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、フーリエ変換処理により求められた周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、周波数特性に生ずるノイズを効果的に除去することが可能となる。このように、フーリエ変換長と測定信号の符号長とが非同期になって、変動するレベルの小さな線スペクトルが周波数特性にノイズとして発生する場合であっても、間引き処理によってノイズを除去することができ、音場の周波数特性における測定精度を向上させることが可能となる。
【００２０】
さらに、測定信号の符号長が短い場合には、周波数特性の測定処理における処理負担および処理時間を低減でき、また処理に必要とされるメモリ量を減らすことができたが、低中域において線スペクトルの周波数間隔が粗くなると共に、検出される線スペクトルが変動するという問題があった。このため測定信号の符号長が短い場合には、十分な測定精度で周波数特性を測定することが容易でないという問題があった。
【００２１】
本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、線スペクトルの周波数間隔が粗い場合であっても、間引き処理により変動するレベルの小さな線スペクトルを取り除くことができるため、周波数特性における低中域の測定精度を十分に確保することが可能になる。
【００２２】
また、本発明に係る音場測定装置は、符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、当該スピーカに対して前記測定信号を出力する外部出力手段と、該スピーカより出力された前記測定信号を収音するマイクロフォンと、該マイクロフォンで収音された測定音を２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換手段と、該フーリエ変換手段により求められた周波数特性を帯域分割することにより、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成する帯域分割手段と、該帯域分割手段によって生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、前記第２周波数特性におけるノイズを除去する間引き手段と、前記帯域分割手段により生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された第１周波数特性を生成する第１平均化処理手段と、前記間引き手段により間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、平均化処理された第２周波数特性を生成する第２平均化処理手段と、前記第１平均化処理手段により平均化処理された前記第１周波数特性と、前記第２平均化処理手段により平均化処理された前記第２周波数特性とを合成することにより、全帯域の信号成分を備えた音場の周波数特性を求める合成手段とを有し、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００２３】
また、本発明に係る音場測定装置の音場測定方法は、符号長が２ｎ−１の周期関数からなる測定信号をスピーカより出力させるために、外部出力手段が、当該スピーカに対して前記測定信号を出力する外部出力ステップと、該外部出力ステップにおいて前記スピーカより出力された前記測定信号を、マイクロフォンで収音する収音ステップと、該収音ステップにおいて前記マイクロフォンで収音された測定音を、フーリエ変換手段が、２ｍのサンプル長でフーリエ変換することにより、周波数特性を求めるフーリエ変換ステップと、
該フーリエ変換ステップにおいて求められた周波数特性を帯域分割することにより、帯域分割手段が、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成する帯域分割ステップと、該帯域分割ステップにおいて生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理することにより、間引き手段が、前記第２周波数特性におけるノイズを除去する間引きステップと、前記帯域分割ステップにおいて生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、第１平均化処理手段が、平均化処理された第１周波数特性を生成する第１平均化処理ステップと、前記間引きステップにおいて間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出することにより、第２平均化処理手段が、平均化処理された第２周波数特性を生成する第２平均化処理ステップと、前記第１平均化処理ステップにおいて平均化処理された前記第１周波数特性と、前記第２平均化処理ステップにおいて平均化処理された前記第２周波数特性とを合成することにより、合成手段が、全帯域の信号成分を備えた音場の周波数特性を求める合成ステップとを有し、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００２４】
さらに、本発明に係る音場測定装置の音場測定プログラムは、符号長が２ｎ−１の周期関数からなる測定信号を用いて音場の周波数特性を測定する音場測定装置の音場測定プログラムであって、該音場測定装置のコンピュータに、スピーカより前記測定信号を出力させるために当該スピーカに対して前記測定信号を出力させる外部出力機能と、該外部出力機能により前記スピーカより出力された前記測定信号を、マイクロフォンを用いて収音させる収音機能と、該収音機能によって収音された測定音を、２ｍのサンプル長でフーリエ変換させることにより、周波数特性を求めさせるフーリエ変換機能と、該フーリエ変換機能によって求められた周波数特性を帯域分割することにより、高域側成分からなる第１周波数特性と、低域側成分からなる第２周波数特性とを生成させる帯域分割機能と、該帯域分割機能によって生成された第２周波数特性より、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引き処理させることにより、前記第２周波数特性におけるノイズを除去させる間引き機能と、前記帯域分割機能により生成された第１周波数特性に基づいて、所定の第１周波数間隔における信号レベルの平均値を、当該第１周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された第１周波数特性を生成させる第１平均化処理機能と、前記間引き機能により間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、当該第２周波数間隔よりも短い周波数ずつシフトさせながら算出させることにより、平均化処理された第２周波数特性を生成させる第２平均化処理機能と、前記第１平均化処理機能により平均化処理された前記第１周波数特性と、前記第２平均化処理機能により平均化処理された前記第２周波数特性とを合成することにより、全帯域の信号成分を備えた前記音場の周波数特性を求めさせる合成機能とを実現させるための音場測定プログラムであって、前記ｎと前記ｍは、ｍ＞ｎを満たす自然数であって、前記ｋは、ｋ＝０，１，２,・・・であることを特徴とする。
【００２５】
本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、フーリエ変換処理により求められた周波数特性を、高域側成分の第１周波数特性と、低域側成分の第２周波数特性とに分割する。そして、低域側の第２周波数特性に対してのみ間引き処理を施すことにより、間引き処理に伴って発生し得る高域成分の信号レベル低減を避けることが可能となる。
【００２６】
また、フーリエ変換長と測定信号の符号長との非同期に伴って生ずる、変動するレベルの小さな線スペクトルは、低中域の周波数間隔が粗いこともあって、低中域でノイズとして判断される可能性が高い。このため、低域側の第２周波数特性に対して間引き処理を行うことによって、低域成分のノイズを効果的に低減することが可能となる。
【００２７】
さらに、高域側の第１周波数特性には間引き処理が行われないため、間引き処理に伴って発生し得る高域側成分の信号レベル低減を回避することが可能となり、高域側成分の増幅を行う必要がなくなる。そして平均化処理された第１周波数特性および第２周波数特性を合成して全帯域の信号成分を備えた周波数特性を生成することにより、より精度良く音場の周波数特性を求めることが可能になる。
【発明の効果】
【００２８】
本発明に係る音場測定装置、音場測定方法および音場測定プログラムでは、線スペクトルの周波数間隔が粗い場合であっても、間引き処理により変動するレベルの小さな線スペクトルを取り除くことができるため、線スペクトルの周波数間隔が粗い場合であっても、周波数特性における低中域の測定精度を十分に確保することが可能になる。
【図面の簡単な説明】
【００２９】
【図１】実施の形態に係る音場測定装置のハードウエア的な概略構成を示したブロック図である。
【図２】実施の形態に係るＣＰＵが処理プログラムに基づいて、周波数特性の測定処理を行う場合における音場測定装置の機能部における第１の概略構成を示したブロック図である。
【図３】実施の形態に係るＣＰＵによる周波数特性の測定処理の内容を示したフローチャートである。
【図４】（ａ）は、実施の形態における間引き処理部で間引き処理を行う前の周波数特性の第１の例を示し、（ｂ）は、間引き処理を行った後の周波数特性の第１の例を示している。
【図５】（ａ）は、実施の形態における間引き処理部で間引き処理を行う前の周波数特性の第２の例を示し、（ｂ）は、間引き処理を行った後の周波数特性の第２の例を示している。
【図６】実施の形態に係る平均化処理部において、間引き処理された信号を１サンプルずつシフトさせながら、所定サンプル数内の平均値を算出する処理を説明するための図である。
【図７】周波数サンプルに応じて設定される平均化処理のサンプル数幅（平均化幅）の関係を示した図である。
【図８】（ａ）は、ループバック方式により、Ｍ系列符号をフーリエ変換部でフーリエ変換処理した場合における周波数特性を示しており、（ｂ）は、（ａ）に示す周波数特性を用いて間引き処理部で間引き処理を行った場合の周波数特性を示している。
【図９】ループバック方式により、平均化処理部で平均化処理された信号の周波数特性を示しており、（ａ）は、間引き処理部で間引き処理を行った後に平均化処理を行った場合を示し、（ｂ）は、間引き処理を行わずに平均化処理を行った場合を示している。
【図１０】音場測定方式により、平均化処理部で平均化処理された信号の周波数特性を示しており、（ａ）は、間引き処理部で間引き処理を行った場合を示し、（ｂ）は、間引き処理を行わなかった場合を示している。
【図１１】実施の形態に係るＣＰＵが処理プログラムに基づいて、周波数特性の測定処理を行う場合における音場測定装置の機能部の第２の概略構成を示したブロック図である。
【図１２】従来の音場測定装置を用いて求められる周波数特性の一例を示しており、（ａ）は、３２，７６７の符号長を持つＭ系列符号を用いた場合の周波数特性を示し、（ｂ）は、１／３オクターブバンド幅で対数的な平均化処理を行った場合の周波数特性を示している。
【図１３】従来の音場測定装置を用いて求められる周波数特性の一例を示しており、（ａ）は、４,０９６の符号長を持つＭ系列符号を用いた場合の周波数特性を示しており、（ｂ）は、１／３オクターブバンド幅で対数的な平均化処理を行った場合の周波数特性を示している。
【発明を実施するための形態】
【００３０】
以下、本発明に係る音場測定装置について、図面を用いて詳細に説明する。図１は、本発明に係る音場測定装置のハードウエア的な概略構成の一例を示したブロック図である。音場測定装置１は、図１に示すように、ＣＰＵ２と、ＲＯＭ（Read Only Memory）３と、ＲＡＭ（Random Access Memory）４と、記憶部５と，外部出力部（外部出力手段）６と、マイクロフォン７と、表示部８とを有している。また外部出力部６は、スピーカ９に接続されている（図２参照）。
【００３１】
ＲＯＭ３には、音場測定装置１において実行される処理プログラム等が記憶されている。例えば、音場測定装置１の起動時、あるいは、ユーザによる操作に応じてＣＰＵ２がＲＯＭ３の処理プログラム等を読み出すことによって、周波数特性の測定等の各種処理を行うことが可能になっている。ＲＡＭ４は、ＣＰＵ２において行われる処理のワークエリア等として用いられる。
【００３２】
記憶部５は、いわゆる補助記憶装置であって、主としてハードディスク、ＳＳＤ（Solid State Drive）、あるいは不揮発性メモリ（フラッシュＲＯＭ・フラッシュメモリ等）等により構成されている。また、ＳＤカードのように取り外し可能なメモリーカードを記憶部５として利用することも可能である。記憶部５には、ＣＰＵ２における各種処理に用いられる様々なデータ等が記録される。
【００３３】
なお、音場測定装置１としてスマートフォンなどの情報携帯端末を利用する場合には、記憶部５にダウンロード等されたアプリケーションプログラムを記録させ、このアプリケーションプログラムに基づいて、周波数特性の測定処理を行うことも可能である。
【００３４】
外部出力部６は、後述する測定信号を、スピーカ９より出力させる役割を有している。外部出力部６は、測定信号をスピーカ９より出力するために必要な機器等で構成されている。例えば、外部出力部６として、測定信号をアナログ信号に変換するＤ／Ａ変換機器や、測定信号の信号出力を増大させるためのアンプ部が該当する。さらに、外部出力部６は、スピーカ９の入力端子にオーディオケーブルを介して接続するための外部出力端子等にも該当することになる。
【００３５】
また、外部出力部６とスピーカ９との接続はオーディオケーブル等を用いた物理的な接続には限定されない。例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）や無線ＬＡＮなどの無線技術を用いることにより測定信号をスピーカ９より出力させる構成であってもよい。
【００３６】
マイクロフォン７は、スピーカ９より出力される測定音を収音する役割を有している。マイクロフォン７により収音された測定音は、ＲＡＭ４あるいは記憶部５に記録され、後述する周波数特性の測定処理に利用される。表示部８は、一般的な液晶ディスプレイやＣＲＴディスプレイ（ブラウン管ディスプレイ）が該当し、周波数特性の測定処理により求められた音場の周波数特性（例えば、後述する図８〜図１０に示す周波数特性）を、ユーザに視認可能に表示させる役割を有している。
【００３７】
ＣＰＵ２は、ＲＯＭ３に記録される処理プログラム、あるいは、記憶部５に記録される周波数特性の測定用アプリケーションプログラムに従って、スピーカ９とマイクロフォン７との間の周波数特性の測定処理を実行する役割を有している。図２は、ＣＰＵ２が処理プログラムあるいはアプリケーションプログラムに基づいて、周波数特性の測定処理を行う場合における音場測定装置１の機能部の概略構成を示したブロック図である。また、図３は、処理プログラム等に基づくＣＰＵ２の処理内容を示したフローチャートである。
【００３８】
音場測定装置１は、図２に示すように、測定信号生成部１１と、フーリエ変換部（フーリエ変換手段）１２と、間引き処理部（間引き手段）１３と、平均化処理部（平均化処理手段）１４と、高域増幅部１５と、外部出力部６と、マイクロフォン７と表示部８とを有している。また、図２には、外部出力部６に接続されたスピーカ９が示されている。なお、外部出力部６とマイクロフォン７と表示部８とについては、既に図１において説明したので、ここでの説明は省略する。
【００３９】
測定信号生成部１１は、測定信号として、任意の生成多項式にてＭ系列符号を生成する。既に説明したように、Ｍ系列符号は、符号長が２ｎ−１の周期関数に該当する。符号長２ｎ−１におけるｎは、自然数である。
【００４０】
ＣＰＵ２は、処理プログラム等に従って、測定信号生成部１１として機能し、Ｍ系列符号からなる測定信号を生成する処理を行う（図３のＳ１）。そして、ＣＰＵ２は、生成されたＭ系列符号を、外部出力部６を用いて、スピーカ９に出力させる（図３のＳ２、外部出力ステップ、外部出力機能）。その後、ＣＰＵ２は、スピーカ９より出力された測定音を、マイクロフォン７を用いて収音させる（図３のＳ３、収音ステップ、収音機能）。収音された測定音の信号（測定信号）は、フーリエ変換部１２へ出力される。
【００４１】
フーリエ変換部１２は、収音された測定信号に対してフーリエ変換処理（高速ＦＦＴ処理）を行う役割を有している。ＣＰＵ２は、フーリエ変換部１２において、収音した測定信号に対して窓関数で重みづけを行った後、フーリエ変換処理を行うことにより、時間領域の測定信号を周波数領域に変換し、フーリエ変換毎に線スペクトルを出力する（図３のＳ４、フーリエ変換ステップ、フーリエ変換機能）。ここで、線スペクトルはパワースペクトルであり、線スペクトル数はフーリエ変換のサンプル長の半分となる。フーリエ変換処理された測定信号は、間引き処理部１３へ出力される。
【００４２】
間引き処理部１３は、求められた周波数特性の線スペクトルから、ノイズとなる線スペクトルを取り除く（間引きする）役割を有している。上述したように、Ｍ系列符号の符号長は２ｎ−１となる。一方で、フーリエ変換処理に求められる線スペクトル数は、１／２・２ｍ（ｍは自然数）となり、フーリエ変換長（フーリエ変換のサンプル長）は、２ｍとなる。一般に、Ｍ系列符号からなる測定信号を収音してフーリエ変換処理する場合には、フーリエ変換長をＭ系列符号の符号長の２倍以上（つまり、ｍ＞ｎ）に設定するが、Ｍ系列符号の符号長は２ｎ−１であるため、フーリエ変換長がＭ系列符号の符号長の整数倍（例えば２倍、４倍、８倍等）の関係にはならない。このようにフーリエ変換長とＭ系列符号の符号長とが整数倍の関係とはならず、非同期となってしまう場合には、フーリエ変換毎に、変動するレベルの小さな線スペクトルが、一様な線スペクトルの間に発生してしまう。この変動するレベルの小さな線スペクトルは、周波数特性の検出においてノイズとなってしまう。このため、間引き処理部１３は、ノイズとなってしまう線スペクトルを取り除く（間引く）ことによって、周波数特性におけるノイズを除去し、測定精度の向上を図る役割を担っている。
【００４３】
次に、間引き処理部１３が間引き処理を行う場合における処理内容について説明する。図４は、符号長が４，０９５（２ｎ−１においてｎ＝１２）のＭ系列符号を測定信号として用い、フーリエ変換長を８，１９２（２ｍにおいてｍ＝１３）に設定して、ループバック方式により、フーリエ変換部１２で求められる線スペクトル（周波数特性）を示している。なお、図４（ａ）は、間引き処理部１３により間引き処理を行う前の周波数特性を示し、図４（ｂ）は、間引き処理部１３により間引き処理を行った後の周波数特性を示している。
【００４４】
ここで、ループバック方式とは、外部出力部６より出力された測定信号をそのままマイクロフォン７で収音された信号としてフーリエ変換部１２へ出力することによって、周波数特性を測定する方法を意味している。このようなループバック方式を用いることにより、測定信号の周波数特性を、音場の影響を受けずにそのままフーリエ変換した状態で示すことができる。このため、測定信号としてＭ系列符号を用いることにより、理想的にはフラットな周波数特性を得ることができ、測定処理におけるノイズ等を容易に判断することが可能になる。
【００４５】
間引き処理部１３は、フーリエ変換部１２により生成された線スペクトルに基づいて、低域側の線スペクトルから順番に、０×２ｍ−ｎ＋１，１×２ｍ−ｎ＋１，２×２ｍ−ｎ＋１，３×２ｍ−ｎ＋１，・・・，ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引く処理を行う。但し、変数ｋはｋ＝０，１，２，３，・・・と１つずつ増加する整数であり、ｋ×２ｍ−ｎ＋１は、フーリエ変換により生成された最後の線スペクトル（高域側の最後の線スペクトル）が含まれる値（最後の線スペクトルの順番≦ｋ×２ｍ−ｎ＋１）である。
【００４６】
図４（ａ）（ｂ）を用いて説明すると、Ｍ系列符号長のｎが１２（ｎ＝１２）であり、フーリエ変換長のｍが１３（ｍ＝１３）であるため、２ｍ−ｎ＝２１３−１２＝２１＝２となる。従って、間引き処理部１３が間引き処理を行う線スペクトルは、図４（ａ）において低域側から０×２＋１，１×２＋１，２×２＋１，３×２＋１，・・・番目以外の線スペクトル、つまり、１，３，５，７，９，・・・番目以外の線スペクトルである。より詳細には、間引き処理部１３により、２，４，６，８，１０，・・・番目の線スペクトルの間引き処理が行われる。
【００４７】
図４（ａ）において、低域側から１，３，５，７，９，・・・番目の線スペクトルは、ノイズとはならないため、信号レベルが０［ｄＢ］となっている。一方で、低域側から２，４，６，８，１０，・・・番目の線スペクトルは、０［ｄＢ］以外の信号レベルを示し、ノイズとして検出されてしまう。このため、図４（ａ）に示す線スペクトル（周波数特性）から、２，４，６，８，１０，・・・番目（ｋ×２１＋１番目以外）の線スペクトルを間引くことにより、図４（ｂ）に示すように、求められた周波数特性から、０［ｄＢ］以外の値を示す線スペクトル、つまり「変動するレベルの小さな線スペクトル」を除去することが可能になる。
【００４８】
図５は、符号長が４，０９５（２ｎ−１においてｎ＝１２）のＭ系列符号を測定信号として用い、フーリエ変換長を１６，３８４（２ｍにおいてｍ＝１４）に設定して、ループバック方式により、フーリエ変換部１２で求められる線スペクトル（周波数特性）を示している。なお、図５（ａ）は、間引き処理部１３により間引き処理を行う前の周波数特性を示し、図５（ｂ）は、間引き処理部１３により間引き処理を行った後の周波数特性を示している。
【００４９】
図５の場合には、Ｍ系列符号長のｎが１２（ｎ＝１２）であり、フーリエ変換長のｍが１４（ｍ＝１４）であるため、２ｍ−ｎ＝２１４−１２＝２2＝４となる。従って、間引き処理部１３が間引き処理を行う線スペクトルは、図５（ａ）において低域側から１，５，９，１３，１７，・・・番目以外の線スペクトルとなる。つまり、低域側から、２，３，４，６，７，８，１０，１１，１２，１４，１５，１６・・・番目の線スペクトルの間引き処理が行われる。
【００５０】
図５（ａ）において、低域側から１，５，９，１３，１７，・・・番目の線スペクトルは、ノイズとはならないため、信号レベルが０［ｄＢ］となっている。一方で、２，３，４，６，７，８，１０，１１，１２，１４，１５，１６・・・番目の線スペクトルは、０［ｄＢ］以外の信号レベルを示し、ノイズとして検出されてしまう。このため、図５（ａ）に示す線スペクトルから、１，５，９，１３，１７，・・・番目（ｋ×２２＋１番目）以外の線スペクトルを間引くことにより、図５（ｂ）に示すように、求められた周波数特性から０［ｄＢ］以外の値を示す線スペクトル（変動するレベルの小さな線スペクトル）を除去することができる。
【００５１】
ＣＰＵ２は、上述したように、フーリエ変換処理により求められた線スペクトルから、ｋ×２ｍ−ｎ＋１番目以外の線スペクトルを間引きする処理を行う（図３のＳ５、間引きステップ、間引き機能）。間引き処理された信号（周波数特性、線スペクトル）は、平均化処理部１４へと出力される。
【００５２】
平均化処理部１４は、間引き処理された信号の所定サンプル数毎の平均値を算出する役割を有している。平均化処理部１４は、図６に示すように、間引き処理された信号の線スペクトルを、低域側から高域側にわたって、１サンプルずつシフトさせながら、所定サンプル数内の平均値を算出する。
【００５３】
図７は、１サンプルずつシフトさせる時の周波数サンプルに応じて設定される、平均値算出のサンプル数幅（平均化幅、所定の周波数間隔）を示した図である。図７において、フーリエ変換のサンプル長が４，０９６、線スペクトル数が２，０４８となっており、図７の横軸の周波数サンプルは、線スペクトル数に対応する数字となっている。図７に示すように、周波数サンプルに応じて平均値算出時の所定サンプル数（所定の周波数間隔）が変化しており、低域側から高域側にわたって平均化幅が増加するようにして設定される。ＣＰＵ２では、図７のような平均化幅の設定により、１／９オクターブ幅の平均値算出処理を行う。聴覚の分解能は１／３オクターブ程度であることが知られている。このため、図７に示すような平均化幅の設定を平均化処理部１４で行うことにより、十分に大きな分解能で平均化処理を行うことが可能となる。
【００５４】
ＣＰＵ２は、平均化処理部１４において、間引き処理部１３で間引き処理された信号の平均化処理を行い（図３のＳ６、平均化処理ステップ、平均化処理機能）、平均化処理された信号を、高域増幅部１５へと出力する。
【００５５】
高域増幅部１５は、平均化処理された信号の高域成分の信号レベルを増幅する役割を有している。平均化処理が行われた信号は、間引き処理に伴って高域成分の信号レベルが減衰する傾向がある。このため、高域増幅部１５では、減衰する高域成分を考慮した逆フィルタを用いて、高域成分の信号レベルを補完して、求められる周波数特性（線スペクトル）の信号レベルがフラットになる（一様な状態となる）ように増幅処理を行う。このような高域成分の増幅処理により、高域成分における周波数特性の測定精度の向上を図ることが可能となる。
【００５６】
ＣＰＵ２では、平均化処理された信号の高域成分を増幅し（図３のＳ７）、高域成分の増幅処理が行われた信号を、表示部８に出力する。なお、ＣＰＵ２では、間引き処理部１３による間引き処理を行うことなく、フーリエ変換部１２により求められた周波数特性をそのまま高域増幅部１５へ出力し、表示部８で表示させることも可能になっている。表示部８において受信された周波数特性（線スペクトル）は、ＣＰＵ２の指示に基づいて、ユーザに視認可能な状態で表示部８のディスプレイ画面等に表示される（図３のＳ８）。
【００５７】
次に、図８〜図１０に示すように、具体的に測定された周波数特性等を示して、音場測定装置１における処理を説明する。まず、図８（ａ）は、ループバック方式により測定されたＭ系列符号をフーリエ変換部１２でフーリエ変換処理した場合（フーリエ変換処理後）における周波数特性を示している。また、図８（ｂ）は、図８（ａ）に示す周波数特性を用いて間引き処理部１３で間引き処理を行った場合（間引き処理後）の周波数特性を示している。
【００５８】
さらに、図９は、ループバック方式により、平均化処理部１４で平均化処理された信号の周波数特性を示しており、図９（ａ）は、間引き処理部１３で間引き処理を行った後に平均化処理を行った場合を示し、図９（ｂ）は、間引き処理を行わずに平均化処理を行った場合を示している。一方で、図１０は、測定信号をスピーカ９より出力し、マイクロフォン７で収音することにより音場の周波数特性を測定する方法（以下、音場測定方式とする）により、平均化処理部１４で平均化処理された信号の周波数特性を示しており、図１０（ａ）は、間引き処理部１３で間引き処理を行った場合を示し、図１０（ｂ）は、間引き処理を行わなかった場合を示している。
【００５９】
図８〜図１０に示す周波数特性の測定条件は、測定信号としてＭ系列符号を用い、測定信号のサンプリング速度を４４．１ｋＨｚ、Ｍ系列符号長を４，０９５、フーリエ変換部１２によるフーリエ変換のサンプル長は８，１９２、フーリエ変換部１２において使用する窓関数はハミング窓、平均化処理部１４における平均化幅は１／９オクターブに設定されている。
【００６０】
Ｍ系列符号の符号長を４，０９５とし、フーリエ変換のサンプル長を８，１９２に設定した場合には、上述したように、Ｍ系列符号の符号長に対してフーリエ変換のサンプル長が整数倍の関係でなく、非同期となってしまう。このため、図８（ａ）に示すように、フーリエ変換毎に、一様な線スペクトルの間に変動するレベルの小さな線スペクトルが発生してしまい、周波数特性の信号レベルに０［ｄＢ］以外のノイズが検出されてしまっている。さらに、符号長が４，０９５のＭ系列符号は、符号長が短い測定信号に該当する。このため、線スペクトルの周波数間隔が荒くなる傾向があり、特に低域成分において検出された線スペクトルに顕著なばらつきが生じ、線スペクトルの包絡線が必ずしも一様な状態にはなっていない。
【００６１】
しかしながら、図８（ｂ）に示すように、間引き処理部１３で、変動するレベルの小さな線スペクトルを間引くことにより、符号長の短いＭ系列符号を測定信号に用いた場合であっても、線スペクトルの信号レベルのばらつきを抑制し、線スペクトルの低域側の包絡線を一様な状態にすることができる。図８（ｂ）では、３，０００Ｈｚ以下の帯域において、信号レベルの変動が抑制されており、周波数特性が一様な状態となっている。しかしながら、３，０００Ｈｚ以上の帯域では、線スペクトルが変動した状態が示されている。
【００６２】
一方で、図９（ａ）に示すように、図８（ｂ）に示した間引き処理された信号に対して、対数的な平均化処理を行うことにより、低中域はもちろんのこと、３，０００Ｈｚ以上の高域成分においても線スペクトルの変動が抑えられている。
【００６３】
図９（ｂ）は、間引き処理を行わずに対数的な平均化処理を行った場合の周波数特性が示されている。図９（ｂ）に示すように、単に平均化処理を行っても、間引き処理を十分に行わなかった場合には、低中域のレベル変動を抑制することができず、周波数特性の測定精度が大きく劣化してしまう。このため、間引き処理部１３による間引き処理により、低中域において変動する信号レベルの小さな線スペクトルを除くことができ、周波数特性の測定精度向上を図ることができると共に、その後に平均化処理を行うことによって、高域成分の線スペクトルの変動も効果的に抑制することが可能となる。
【００６４】
なお、図９（ａ）に示すように、間引き処理を行うことによって、高域成分の信号レベルが低減した状態となるが、高域増幅部１５において高域成分の増幅処理を行うことによって、高域の減衰量を補完することが可能であり、測定信号の周波数測定を一様にフラットな状態とすることができる。
【００６５】
図１０（ａ）および（ｂ）は、音場測定方式を用いた場合の周波数特性を示している。図１０（ａ）（ｂ）では、スピーカ９から出力された測定音をマイクロフォン７で収音して周波数特性を測定しているため、音場（マイクロフォン７の設置位置における音場）の周波数特性が測定されることになる。図１０（ａ）では、間引き処理を行った後に平均化処理を行っているので、低中域におけるレベル変動が効果的に抑制されているが、図１０（ｂ）では、間引き処理を行わずに平均化処理を行っているので、低中域におけるレベル変動を抑えることができず、音場における周波数特性の測定精度が大きく劣化した状態となっている。
【００６６】
以上説明したように、本実施の形態に係る音場測定装置１では、間引き処理部１３において、フーリエ変換処理により求められた線スペクトルの間引き処理を行う。この間引き処理により、Ｍ系列符号の符号長に対してフーリエ変換のサンプル長が非同期となることにより生ずる「変動するレベルの小さな線スペクトル」を取り除くことができ、周波数特性の測定精度を高めることが可能になる。
【００６７】
特に、間引き処理部１３において間引き処理を行う場合には、測定信号の符号長が２ｎ−１、フーリエ変換処理のサンプル長が２ｍであるときに、ｋ×２ｍ−ｎ＋１番目の以外の線スペクトルを間引くことにより、効果的に変動するレベルの小さな線スペクトルを取り除くことが可能となる。
【００６８】
さらに、間引き処理を行うことにより、測定信号の符号長が短い場合であって、求められた周波数特性における線スペクトル（周波数スペクトル）の周波数間隔が粗い場合であっても、効果的に変動するレベルの小さな線スペクトルを取り除くことができる。このため、符号長の短い測定信号を用いた場合であっても、周波数特性の測定精度を十分に確保することができるとともに、周波数特性の測定に必要とされる測定時間や測定負担を低減し、処理に必要なメモリ量を効果的に減らすことが可能となる。
【００６９】
また、対数的な平均化処理を行うことにより、全帯域において線スペクトルの変動を抑制することができ、音場の周波数特性の測定精度をさらに向上させることが可能となる。
【００７０】
以上、本発明に係る音場測定装置、音場測定方法および音場測定プログラムについて、図面を用いて詳細に説明を行ったが、本発明に係る音場測定装置、音場測定方法および音場測定プログラムは、実施の形態で説明した例には限定されない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても本発明の技術的範囲に属するものと了解される。
【００７１】
例えば、上述した実施の形態では、フーリエ変換処理により求められた周波数特性の全ての帯域に対して間引き処理を行う場合を一例として示して説明した。しかしながら、全ての帯域に対して間引き処理を行うと、高域成分における信号レベルが低減する傾向があり、この低減した高域成分の信号レベルを増幅するために、高域増幅部１５を設ける構成を採用した。
【００７２】
しかしながら、間引き処理を行う帯域を、変動するレベルの小さな線スペクトルの影響が顕著な低中域だけに限定することにより、高域成分における信号レベルの増幅処理を行う必要性が低くなる。
【００７３】
図１１は、フーリエ変換処理により求められる周波数特性の低域成分にだけ間引き処理を行い、高域成分に対しては間引き処理を行わないことを特徴とする音場測定装置１ａの概略構成を示した図である。図１１において、図２に示した構成と同様の処理を行う箇所については同一の符号を附してある。図２に示す音場測定装置１と、図１１に示す音場測定装置１ａとを比較すると、図１１には、帯域分割処理部（帯域分割手段）２０と、ゲイン部２１と、合成処理部（合成手段）２２とが設けられる一方で、図２に示した高域増幅部１５が設けられていない点で相違する。また、図１１における第１平均化処理部（第１平均化処理手段）１４ａおよび第２平均化処理部（第２平均化処理手段）１４ｂは、周波数特性に対して平均化処理を行う点で、平均化処理部１４と同様の処理を行う。
【００７４】
図１１に示す音場測定装置１ａにおいて、帯域分割処理部２０は、フーリエ変換部１２のフーリエ変換処理により求められた周波数特性を、高域側成分からなる周波数特性と、低域側成分からなる周波数特性とに分割する役割を有している。帯域分割処理部２０による分割処理では、フーリエ変換部１２より入力される信号を、予め決定された所定の周波数を境にして、高域側成分の第１周波数特性の信号と、低域側成分の第２周波数特性の信号とに分割する（帯域分割ステップ、帯域分割機能）。この分割処理は、ハイパスフィルタおよびローパスフィルタ等のフィルタを用いる帯域分割ではなく、信号のデジタル処理等を用いることにより、境となる所定の周波数値で２つに分割処理される。従って、帯域分割処理部２０で分割された第１周波数特性の信号は、所定の周波数値以上の周波数における信号レベルしか有しておらず、第２周波数特性の信号は、所定の周波数値以下の周波数における信号レベルしか有していない。
【００７５】
このようにして分割された低域の第２周波数特性に対してのみ、間引き処理部１３で間引き処理を施し、第２平均化処理部１４ｂで平均化処理を行う（第２平均化処理ステップ、第２平均化処理機能）。このように低域の第２周波数特性に対してのみ間引き処理を行うことによって、間引き処理に伴って発生し得る高域成分の信号レベル低減を避けることが可能となる。なお、第２平均化処理部１４ｂは、間引き処理された第２周波数特性に基づいて、所定の第２周波数間隔における信号レベルの平均値を、第２周波数間隔よりも短い周波数ずつ、例えば１サンプルずつ、シフトさせながら算出することにより、平均化処理された第２周波数特性を生成する。
【００７６】
一方で、分割された高域の第１周波数特性に対しては、間引き処理を行うことなく第１平均化処理部１４ａで平均化処理を行う（第１平均化処理ステップ、第１平均化処理機能）。さらに、ゲイン部２１において、第２周波数特性との信号レベル差を考慮したゲイン調整を行う。高域の周波数特性に対しては、間引き処理が行われないので、間引き処理による高域成分の信号レベル低減が発生せず、図２に示すような高域増幅部１５を設ける必要がなくなる。
【００７７】
また、第１平均化処理部１４ａは、間引き処理されていない第１周波数特性の信号に基づいて、所定の第１周波数間隔における信号レベルの平均値を、第１周波数間隔よりも短い周波数ずつ、例えば１サンプルずつ、シフトさせながら算出することにより、平均化処理された第１周波数特性を生成する。
【００７８】
そして、ゲイン部２１によりゲイン調整された高域の第１周波数特性と、第２平均化処理部１４ｂにより平均化処理された低域の第２周波数特性とは、合成処理部２２において合成処理されて、全帯域の信号成分を備えた周波数特性が生成される（合成ステップ、合成機能）。合成処理部２２では、低域側を第２周波数特性により構成され、高域側を第１周波数特性により構成される全帯域の周波数特性を、合成処理によって生成する。
【００７９】
このようにして合成・生成された周波数特性は、低域成分に対してのみ間引き処理が行われて、変動するレベルの小さな線スペクトルが効果的に低減されているので、低域のノイズが抑制された周波数特性を実現することが可能となる。さらに、高域成分については間引き処理が行われていないので、平均化処理後に高域成分の増幅処理を行う必要がなくなり、十分な測定精度を確保した周波数特性を求めることが可能となる。
【００８０】
また、本実施の形態に係る音場測定装置１では、図１に示すようにＲＯＭ３や記憶部５に記録される処理プログラムあるいはアプリケーションプログラムに基づいて、ＣＰＵ２が図２に示すような機能部の機能を実現する旨の説明を行ったが、各機能部の機能を実現するＣＰＵは、１つには限定されない。各機能部の一部の機能を実現するための専用の処理手段（例えば、特定の処理専用のＣＰＵ，チップ等）を設置し、各専用の処理手段が少なくとも１以上の機能を実現する構成であってもよい。このように、複数の専用の処理手段を設定した場合であっても、１つのＣＰＵが処理プログラム等に基づいて音場測定処理を行う場合であっても、間引き処理により変動するレベルの小さな線スペクトルを取り除くことによって、ノイズを効果的に低減することができ、符号長の短い測定信号を用いた場合であっても、精度良く音場環境の周波数特性を測定することが可能である。
【符号の説明】
【００８１】
１，１ａ       …音場測定装置
２       …ＣＰＵ
３       …ＲＯＭ
４       …ＲＡＭ
５       …記憶部
６       …外部出力部（外部出力手段）
７       …マイクロフォン
８       …表示部
９       …スピーカ
１１     …測定信号生成部
１２     …フーリエ変換部（フーリエ変換手段）
１３     …間引き処理部（間引き手段）
１４     …平均化処理部（平均化処理手段）
１４ａ  …第１平均化処理部（第１平均化処理手段）
１４ｂ  …第２平均化処理部（第２平均化処理手段）
１５     …高域増幅部
２０     …帯域分割処理部（帯域分割手段）
２１     …ゲイン部
２２     …合成処理部（合成手段）
【図１】
2015075345.tif 000003
【図２】
2015075345.tif 000004
【図３】
2015075345.tif 000005
【図４】
2015075345.tif 000006
【図５】
2015075345.tif 000007
【図６】
2015075345.tif 000008
【図７】
2015075345.tif 000009
【図８】
2015075345.tif 000010
【図９】
2015075345.tif 000011
【図１０】
2015075345.tif 000012
【図１１】
2015075345.tif 000013
【図１２】
2015075345.tif 000014
【図１３】
2015075345.tif 000015