(57)【要約】      （修正有）【課題】スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置及びプログラムを提供する。【解決手段】聴取位置で測定された再生周波数応答（振幅周波数特性）とその振幅周波数特性を多重解像度解析した結果から、パラメトリックイコライザ２０のパラメータ（中心周波数、ゲイン、帯域幅）を算出する。再生周波数分析装置１０は、測定信号発生部１０１と、相互相関関数測定部１０３と、遅延部１０２と、コヒーレンス関数測定部１０４と、周波数応答関数測定部１０５と、理想周波数特性データベース１０７と、差分部１０６と、多重解像度分析部１０８と、有効帯域判定部１０９と、イコライザパラメータ計算部１１０とを備える。【選択図】図１
2012028874.tif 000002
【特許請求の範囲】
【請求項１】
スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置であって、
所定の音源から得られる測定信号を発生する測定信号発生部と、
該測定信号発生部から発生した測定信号に対して予め定めた中心周波数、帯域幅及びゲインをパラメータとする周波数特性のイコライジングを行って前記スピーカから発生させた音を予め定めた聴取位置に設けたマイクを介して収音し、マイクアンプを介して増幅して得られる観測信号と、前記測定信号発生部から得られる測定信号との間で相互相関関数に基づく遅延量を測定する相互相関関数測定部と、
該遅延量に従って前記測定信号発生部から得られる測定信号を遅延させ、前記相互相関関数測定部に再入力し、当該相互相関関数の繰返し演算に基づく遅延調整を行い、遅延調整後の測定信号を生成する遅延部と、
該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、コヒーレンス関数を測定するコヒーレンス関数測定部と、
該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、周波数応答関数を測定する周波数応答関数測定部と、
前記音響再生システムに対する聴取位置を含む室空間に応じて理想周波数特性のデータを予め格納しておく理想周波数特性データベースと、
該理想周波数特性データベースに格納された前記理想周波数特性を取り出し、前記周波数応答関数との差分を算出し、差分周波数応答を生成する差分部と、
該差分部より得られた差分周波数応答の多重解像度解析を行い、中心周波数及びスケールパラメータを対応付けた多重解像度の分析値を算出する多重解像度分析部と、
前記コヒーレンス関数で所定値以上を有する周波数に該当する周波数帯域を、最適化対象の有効周波数帯域として判定する有効帯域判定部と、
前記有効周波数帯域に該当する前記多重解像度解析の分析値の絶対値が所定値以上となる中心周波数と、前記スケールパラメータから得られる該中心周波数に対応する帯域幅と、前記周波数応答関数から該中心周波数に対応するゲインとを算出し、イコライジング用のパラメータとして決定するイコライザパラメータ計算部と、
を備えることを特徴とする再生周波数分析装置。
【請求項２】
前記差分周波数応答の値の絶対値が有効周波数帯域における全ての周波数にわたり所定値以下となるまで、前記音響再生システムに対する再生周波数の分析を繰り返す手段を有することを特徴とする、請求項１に記載の再生周波数分析装置。
【請求項３】
前記多重解像度分析部は、前記分析結果のピーク又はディップの絶対値の最大のものから順に所定個数分の調整すべき中心周波数を決定し、該決定した中心周波数におけるイコライジング用のパラメータを決定する手段を有することを特徴とする、請求項２に記載の再生周波数分析装置。
【請求項４】
スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置として構成するコンピュータに、
所定の音源から得られる測定信号を発生するステップと、
該発生した測定信号に対して予め定めた中心周波数、帯域幅及びゲインをパラメータとする周波数特性のイコライジングを行って前記スピーカから発生させた音を予め定めた聴取位置に設けたマイクを介して収音し、マイクアンプを介して増幅して得られる観測信号と、前記測定信号との間で相互相関関数に基づく遅延量を測定するステップと、
該遅延量に従って前記測定信号を遅延させ、前記相互相関関数測定部に再入力し、当該相互相関関数の繰返し演算に基づく遅延調整を行い、遅延調整後の測定信号を生成するステップと、
該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、コヒーレンス関数を測定するステップと、
該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、周波数応答関数を測定するステップと、
前記音響再生システムに対する聴取位置を含む室空間に応じて理想周波数特性のデータを予め格納しておく理想周波数特性データベースに格納された前記理想周波数特性を取り出し、前記周波数応答関数との差分を算出し、差分周波数応答を生成するステップと、
該ステップより得られた差分周波数応答の多重解像度解析を行い、中心周波数及びスケールパラメータを対応付けた多重解像度の分析値を算出するステップと、
前記コヒーレンス関数で所定値以上を有する周波数に該当する周波数帯域を、最適化対象の有効周波数帯域として判定するステップと、
前記有効周波数帯域に該当する前記多重解像度解析の分析値の絶対値が所定値以上となる中心周波数と、前記スケールパラメータから得られる該中心周波数に対応する帯域幅と、前記周波数応答関数から該中心周波数に対応するゲインとを算出し、イコライジング用のパラメータとして決定するステップと、
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音源からの音をスピーカで再生する際に、聴取位置における再生周波数特性を調整する技術に関し、特に、スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置及びそのプログラムに関する。
【背景技術】
【０００２】
従来、室内の各所に聴取者を取り囲むように複数（２個以上の）スピーカを設置し、これらのスピーカを駆動することで、立体的な音響空間を再現できることが知られている。例えば、２２．２マルチチャンネル音響のような多数のスピーカを駆動することで高さ方向も含めた立体的な音響空間の再現が可能となる。
【０００３】
一方、スピーカを利用して音響を再生する音響再生システムにおいて、カラレーションやハウリングが発生しないような音響特性となるように、当該音響の周波数特性を調整するためのイコライザが用いられる技術が知られている（例えば、特許文献１参照）。
【０００４】
イコライザのパラメータは、一般的に周波数特性上のピークを無くすように設定されるものであり、イコライジングを行う場合、このピークのゲインが所定の値（目標レベル）以下となるようにパラメトリックイコライザの特性を決定するにあたり、中心周波数と、中心周波数に対するゲイン及びＱ値（共振の鋭さを表す量であり、「帯域幅」とも称する）がイコライザのパラメータとなる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００７−２０６６３６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、スピーカを設置した室空間の影響により、スピーカの音質が均一に聴取者へ到達することはほとんどありえず、スピーカの位置と部屋の関係により、聴取位置で音圧及び周波数特性がスピーカごとに変化する。
【０００７】
特に、多くのスピーカを駆動する音響再生システムにおいては、聴取位置において、それぞれのスピーカの音質がばらついてしまうと、元の音響信号を正確に再生することができない。例えば、音像移動を再生する場合、音響信号が時々刻々とスピーカを変えて再生されるが、各スピーカの音質が揃っていないと、聴取位置において、再生スピーカが切り替わるごとに異なる音質が再生され、定位感や臨場感に影響を与えてしまう。従って、このような音響再生システムで再生するには全てのスピーカによる再生周波数特性を揃える必要がある。
【０００８】
そこで、本発明の目的は、上述の問題に鑑みて、スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置及びそのプログラムを提供することにある。
【課題を解決するための手段】
【０００９】
音響信号が音響再生システム（帯域により複数のスピーカユニットに分かれているものを含む）から再生され、ある室空間を伝搬し、受聴される場合、再生されるスピーカ位置や空間の形状などにより、聴取位置の再生周波数特性が所望の特性とならないときに、パラメトリックイコライザにより、再生周波数特性を所望の特性になるよう調整する必要がある。そこで、本発明は聴取位置で測定された再生周波数応答（振幅周波数特性）とその振幅周波数特性を多重解像度解析した結果から、パラメトリックイコライザのパラメータ（中心周波数、ゲイン、帯域幅）を算出する。このパラメータを使用して再生周波数特性を所望の特性に調整可能にする。
【００１０】
即ち、本発明の再生周波数分析装置は、スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置であって、所定の音源から得られる測定信号を発生する測定信号発生部と、該測定信号発生部から発生した測定信号に対して予め定めた中心周波数、帯域幅及びゲインをパラメータとする周波数特性のイコライジングを行って前記スピーカから発生させた音を予め定めた聴取位置に設けたマイクを介して収音し、マイクアンプを介して増幅して得られる観測信号と、前記測定信号発生部から得られる測定信号との間で相互相関関数に基づく遅延量を測定する相互相関関数測定部と、該遅延量に従って前記測定信号発生部から得られる測定信号を遅延させ、前記相互相関関数測定部に再入力し、当該相互相関関数の繰返し演算に基づく遅延調整を行い、遅延調整後の測定信号を生成する遅延部と、該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、コヒーレンス関数を測定するコヒーレンス関数測定部と、該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、周波数応答関数を測定する周波数応答関数測定部と、前記音響再生システムに対する聴取位置を含む室空間に応じて理想周波数特性のデータを予め格納しておく理想周波数特性データベースと、該理想周波数特性データベースに格納された前記理想周波数特性を取り出し、前記周波数応答関数との差分を算出し、差分周波数応答を生成する差分部と、該差分部より得られた差分周波数応答の多重解像度解析を行い、中心周波数及びスケールパラメータを対応付けた多重解像度の分析値を算出する多重解像度分析部と、前記コヒーレンス関数で所定値以上を有する周波数に該当する周波数帯域を、最適化対象の有効周波数帯域として判定する有効帯域判定部と、前記有効周波数帯域に該当する前記多重解像度解析の分析値の絶対値が所定値以上となる中心周波数と、前記スケールパラメータから得られる該中心周波数に対応する帯域幅と、前記周波数応答関数から該中心周波数に対応するゲインとを算出し、イコライジング用のパラメータとして決定するイコライザパラメータ計算部と、を備えることを特徴とする。
【００１１】
また、本発明の再生周波数分析装置において、前記差分周波数応答の値の絶対値が有効周波数帯域における全ての周波数にわたり所定値以下となるまで、前記音響再生システムに対する再生周波数の分析を繰り返す手段を有することを特徴とする。
【００１２】
また、本発明の再生周波数分析装置において、前記多重解像度分析部は、前記分析結果のピーク又はディップの絶対値の最大のものから順に所定個数分の調整すべき中心周波数を決定し、該決定した中心周波数におけるイコライジング用のパラメータを決定する手段を有することを特徴とする。
【００１３】
更に、本発明のプログラムは、スピーカを利用して音響を再生する音響再生システムに対して再生周波数を分析し、最適化したイコライジング用のパラメータを決定する再生周波数分析装置として構成するコンピュータに、所定の音源から得られる測定信号を発生するステップと、該発生した測定信号に対して予め定めた中心周波数、帯域幅及びゲインをパラメータとする周波数特性のイコライジングを行って前記スピーカから発生させた音を予め定めた聴取位置に設けたマイクを介して収音し、マイクアンプを介して増幅して得られる観測信号と、前記測定信号との間で相互相関関数に基づく遅延量を測定するステップと、該遅延量に従って前記測定信号を遅延させ、前記相互相関関数測定部に再入力し、当該相互相関関数の繰返し演算に基づく遅延調整を行い、遅延調整後の測定信号を生成するステップと、該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、コヒーレンス関数を測定するステップと、該遅延調整後の測定信号と前記マイクアンプを介して増幅して得られる観測信号から、周波数応答関数を測定するステップと、前記音響再生システムに対する聴取位置を含む室空間に応じて理想周波数特性のデータを予め格納しておく理想周波数特性データベースに格納された前記理想周波数特性を取り出し、前記周波数応答関数との差分を算出し、差分周波数応答を生成するステップと、該ステップより得られた差分周波数応答の多重解像度解析を行い、中心周波数及びスケールパラメータを対応付けた多重解像度の分析値を算出するステップと、前記コヒーレンス関数で所定値以上を有する周波数に該当する周波数帯域を、最適化対象の有効周波数帯域として判定するステップと、前記有効周波数帯域に該当する前記多重解像度解析の分析値の絶対値が所定値以上となる中心周波数と、前記スケールパラメータから得られる該中心周波数に対応する帯域幅と、前記周波数応答関数から該中心周波数に対応するゲインとを算出し、イコライジング用のパラメータとして決定するステップと、を実行させるためのプログラムである。
【発明の効果】
【００１４】
本発明によれば、聴取位置で測定された再生周波数応答（振幅周波数特性）とその振幅周波数特性を多重解像度解析した結果から、パラメトリックイコライザの各パラメータ（中心周波数、帯域幅、ゲイン）を算出することが可能となり、音響再生システムの周波数特性を最適に再生できるようになる。
【図面の簡単な説明】
【００１５】
【図１】本発明による一実施例の再生周波数分析装置と被調整システムを示すブロック図である。
【図２】本発明による一実施例の再生周波数分析装置におけるイコライジングパラメータによるイコライジングを例示する図である。
【図３】本発明による一実施例の再生周波数分析装置におけるウェーブレット（wavelet）多重解像度分析に用いるマザーウェーブレットの算出例を示す図である。
【図４】本発明による一実施例の再生周波数分析装置におけるウェーブレット（wavelet）多重解像度分析の結果を例示する図である。
【発明を実施するための形態】
【００１６】
以下、本発明による一実施例の再生周波数分析装置について説明する。
【００１７】
図１は、本発明による一実施例の再生周波数分析装置と被調整システムを示すブロック図である。図１を参照するに、所定の音響信号の再生周波数特性を調整するためのパラメトリックイコライザ２０と、再生周波数特性を調整した音響信号を増幅するための増幅器３０と、増幅した音響信号を再生するためのスピーカ４０とを備える、室空間で用いる音響再生システム（帯域により複数のスピーカユニットに分かれているものを含む）を、本発明に係る「被調整システム」として説明する。ここで、パラメトリックイコライザ２０は、パラメータ（中心周波数、ゲイン、帯域幅）で周波数特性を調整する装置である。
【００１８】
本実施例の再生周波数分析装置１０は、被調整システムにおけるスピーカ４０から発生する音を調整するために、所定の音源から得られる測定信号を被調整システムにおけるパラメトリックイコライザ２０に供給し、スピーカ４０を介して再生される音の再生周波数を分析し、最適化したイコライジング用のパラメータを決定する装置である。聴取位置Ａには、スピーカ４０を介して再生される音を収音するためのマイク５０が設置され、マイク５０で収音した音声信号は、マイクアンプ６０にて増幅され、再生周波数分析装置１０へと入力される。
【００１９】
再生周波数分析装置１０は、マイクアンプ６０にて増幅した音声信号について再生周波数応答（振幅周波数特性）とその振幅周波数特性を多重解像度解析し、この多重解像度解析した結果から、最適化したパラメトリックイコライザのパラメータ（中心周波数、ゲイン、帯域幅）を算出する。この算出したパラメータを使用してパラメトリックイコライザ２０のパラメータを調整し、所定の音響信号の再生周波数特性を所望の特性に調整することができる。
【００２０】
図１に示すように、再生周波数分析装置１０は、測定信号発生部１０１と、遅延部１０２と、相互相関関数測定部１０３と、コヒーレンス関数測定部１０４と、周波数応答関数測定部１０５と、差分部１０６と、理想周波数特性データベース１０７と、ウェーブレット（wavelet）多重解像度分析部１０８と、有効帯域判定部１０９と、イコライザパラメータ計算部１１０とを備える。
【００２１】
測定信号発生部１０１は、所定の音源から得られる測定信号（例えば、ピンクノイズ）を発生し、パラメトリックイコライザ２０及び遅延部１０２に送出する機能を有する。
【００２２】
パラメトリックイコライザ２０は、測定信号発生部１０１から発生した測定信号を、遅延部１０２を介して入力し、この測定時には、予め定めた初期値の中心周波数、帯域幅及びゲインを初期パラメータとする周波数特性のイコライジングを行って、増幅器３０を介してスピーカ４０から音を発生させる。この測定信号に対応する音は、予め定めた聴取位置Ａに設けたマイク５０を介して収音され、マイクアンプ６０によって増幅される。
【００２３】
遅延部１０２は、測定信号発生部１０１から発生した測定信号について初期遅延量を０として相互相関関数測定部１０３に送出し、さらに、相互相関関数測定部１０３から得られる遅延量を設定して当該測定信号を遅延させる動作を、相互相関関数測定部１０３から得られる遅延量の変化が所定値以下（例えば、ゼロ）となるまで繰り返す機能を有する。従って、相互相関関数測定部１０３から得られる遅延量が安定すると、遅延部１０２は、室空間の影響を遅延調整した測定信号を発生する。
【００２４】
相互相関関数測定部１０３は、マイクアンプ６０を介して増幅して得られる観測信号と、測定信号発生部１０１から得られる測定信号との間で相互相関関数に基づく遅延量を測定する機能を有する。
【００２５】
コヒーレンス関数測定部１０４は、遅延部１０２からの遅延調整後の測定信号と、マイクアンプ６０を介して増幅して得られる観測信号から、コヒーレンス関数を測定する機能を有する。
【００２６】
周波数応答関数測定部１０５は、遅延部１０２からの遅延調整後の測定信号と、マイクアンプ６０を介して増幅して得られる観測信号から、周波数応答関数を測定する機能を有する。
【００２７】
差分部１０６は、理想周波数特性データベース１０７に格納された１つ以上の理想周波数特性を取り出し、周波数応答関数測定部１０５から得られる周波数応答関数との差分を算出し、差分周波数応答を生成する機能を有する。
【００２８】
理想周波数特性データベース１０７は、被調整システムに対する聴取位置Ａを含む室空間に応じて１つ以上の理想周波数特性のデータを予め格納しておく機能を有する。理想周波数特性は、ユーザが所望する周波数特性として予め設定されるものであり、例えば、図２に示すような周波数応答のデータを格納しておく。また、理想周波数特性データベース１０７に格納される１つ以上の理想周波数特性は、実際に各ユーザが聴取する音と、マイク５０及びマイクアンプ６０を介して得られる音との音圧差を吸収するのに用いる。
【００２９】
ウェーブレット（wavelet）多重解像度分析部１０８は、差分部１０６より得られた差分周波数応答の多重解像度解析を行い、中心周波数及びスケールパラメータを対応付けた多重解像度の分析値を算出する機能を有する。
【００３０】
有効帯域判定部１０９は、コヒーレンス関数測定部１０４から得られるコヒーレンス関数で所定値以上を有する周波数に該当する周波数帯域を、最適化対象の有効周波数帯域として判定する機能を有する。
【００３１】
イコライザパラメータ計算部１１０は、有効帯域判定部１０９によって判定した有効周波数帯域に該当する当該多重解像度解析の分析値の絶対値が所定値以上となる中心周波数と、スケールパラメータから得られる該中心周波数に対応する帯域幅を特定し、周波数応答関数測定部１０５から得られる周波数応答関数から該中心周波数に対応するゲインを算出し、パラメトリックイコライザ２０におけるイコライジング用のパラメータとして決定する機能を有する。決定したパラメータは、前述した初期パラメータ（又は前回調整時のパラメータ）を更新するのに用いられ、予め定めた回数（１回の測定でもよい）を繰返し行うこともできる。
【００３２】
これにより、室空間に応じて音響再生システムの周波数特性を最適化することができ、パラメトリックイコライザ２０を介して得られる任意の音響信号を最適に再生できるようになる。
【００３３】
尚、本実施例の再生周波数分析装置１０をコンピュータで構成した場合、各機能を実現する処理内容を記述したプログラムを、当該コンピュータの内部又は外部の記憶部に格納しておき、当該コンピュータの中央演算処理装置（ＣＰＵ）によってこのプログラムを読み出して実行させることで、再生周波数分析装置１０を実現することもできる。また、そのようなプログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラム又はサーバから転送されたプログラムを、一旦、自己の記憶部に記憶することができる。また、このプログラムの別の実施態様として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、更に、このコンピュータにサーバからプログラムが転送される度に、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。
【００３４】
以下、本実施例の再生周波数分析装置１０の各機能ブロックについて更に詳細に説明する。
【００３５】
まず、測定信号発生部１０１は、所定の測定信号を、パラメトリックイコライザ２０、増幅器３０及び１つ以上のスピーカ４０で構成される被調整システムへ入力する。この測定信号は、広帯域で全帯域に対してＳ／Ｎのよい信号であれば、任意の信号を採用することができるが、この要求条件を満足する信号としてピンクノイズを使用することが望ましい。この測定信号は、デジタル信号の測定信号ｘ（ｎ）として出力され、被調整システムの入力段に設けられるＡ／Ｄ変換器（図示せず）によりアナログ信号となる。
【００３６】
被調整システムから再生された測定信号は、聴取位置Ａにあるマイク５０により収音され、マイクアンプ６０を通じ、観測信号ｙ（ｔ）が得られる。この観測信号ｙ（ｔ）は、マイクアンプ６０の直後に設けられるＤ／Ａ変換器（図示せず）によりデジタル変換され、観測信号ｙ（ｎ）となり、相互相関関数測定部１０３、コヒーレンス関数測定部１０４、及び、周波数応答測定部１０５の各々へ入力される。ここでは、このデジタル化された観測信号ｙ（ｎ）は、Ｌ個のサンプリングレート（Ｌサンプル）を１フレームとしてＮフレームごとに相互相関関数測定部１０３、コヒーレンス関数測定部１０４、及び、周波数応答測定部１０５における計算を行うように構成される。
【００３７】
測定信号発生部１０１からの測定信号ｘ（ｎ）は、遅延部１０２を通じて、相互相関関数測定部１０３、コヒーレンス関数測定部１０４、及び周波数応答測定部１０３に入力される。遅延部１０２の動作直後（初期値）は、遅延量を０とする。
【００３８】
相互相関関数測定部１０３は、測定信号ｘ（ｎ）と観測信号ｙ（ｎ）の相互相関関数Ｒｘｙ（ｊ）を算出する。一般的には式（１）で算出することができる。Ｌ個のサンプルについて、ｊ＝０，１，２，・・・，Ｌ−１とする。
【００３９】
【数１】
2012028874.tif 000003
【００４０】
より具体的には、測定信号ｘ（ｎ）と観測信号ｙ（ｎ）のそれぞれのフーリエ変換をＸ（ｆ），Ｙ（ｆ）とすると、式（２）に示すように、Ｘ（ｆ），Ｙ（ｆ）のクロススペクトルＷｘｙについて逆フーリエ変換を実行することにより相互相関関数Ｒｘｙ（ｊ）を算出することができる。
【００４１】
【数２】
2012028874.tif 000004
【００４２】
ここで、〈 〉Nは、Ｎフレーム分の平均を表し、Ｘ＊は、Ｘの複素共役を表す。この相互相関関数の最大値の時刻τが、測定信号ｘ（ｎ）と観測信号のｙ（ｎ）の時間差（遅延量）となる。
【００４３】
相互相関関数の最大値の時刻τは、遅延量として遅延部１０２へ送られる。
【００４４】
遅延部１０２では、相互相関関数測定部１０３から送られた遅延量の時刻τで遅延させた測定信号ｘ（ｎ−τ）を出力する。この測定信号ｘ（ｎ−τ）により、相互相関関数測定部１０３が再度、相互相関関数を算出すると、相互相関関数測定部１０３で算出される相互相関関数の最大値の時刻はほぼ０となる。
【００４５】
尚、相互相関関数測定部１０３は、相互相関関数の最大値の時刻が０となるように繰返し遅延量を算出するように構成され、この測定期間中はいつも最大値の時刻が０となるよう遅延部１０２の遅延量は制御される。以下の説明では、遅延量が制御されたものとして測定信号ｘ（ｎ−τ）→測定信号ｘ（ｎ）とする。
【００４６】
コヒーレンス関数測定部１０４は、測定信号ｘ（ｎ）と観測信号ｙ（ｎ）を入力として式（３）によりコヒーレンス関数γ２（ｆ）を算出する。
【００４７】
【数３】
2012028874.tif 000005
【００４８】
ここで、Ｗｘｙはクロススペクトル、Ｗｘｘ，Ｗｙｙは、それぞれ測定信号ｘ（ｎ）と観測信号ｙ（ｎ）のパワースペクトルである。コヒーレンス関数の算出には、Ｎフレーム分の平均化を必要とする。
【００４９】
コヒーレンス関数測定部１０４で算出されたコヒーレンス関数γ２（ｆ）は、有効帯域判定部１０９へ送られる。
【００５０】
有効帯域判定部１０９では、入力されたコヒーレンスが０から１の値をとり、１に近い値をとる周波数は、系の出力が全て測定入力に起因していることを示しており、ノイズや歪みなどの影響が少ない信号であると考えられる。コヒーレンス値を基準のコヒーレンス値ref_γ２と比較し、基準コヒーレンス値ref_γ２よりも大きな値をとる周波数に対して１（測定の信頼度が高く調整可能な周波数）、小さな値をとる周波数に対して０（測定の信頼度が低く調整不要な周波数）として有効周波数帯域を決定し、この有効周波数帯域の情報を有効周波数帯域判定信号としてイコライザパラメータ計算部１１０に出力する。
【００５１】
イコライザパラメータ計算部１１０には、有効帯域判定部１０９から得られる有効周波数帯域判定信号と、ウェーブレット（wavelet）多重解像度分析部１０８からの分析結果を示す出力が入力される。
【００５２】
周波数応答関数測定部１０５は、測定信号ｘ（ｎ）と観測信号ｙ（ｎ）を入力として式（４）により周波数応答関数Ｈ（ｆ）を算出する。
【００５３】
【数４】
2012028874.tif 000006
【００５４】
差分部１０６は、周波数応答関数測定部１０５によって算出された周波数応答関数Ｈ（ｆ）の振幅周波数特性と、理想周波数特性データベース１０７に予め格納される、ユーザが所望する周波数特性ＨＵ（ｆ）とを、周波数ごとに差分を行った差分周波数応答Ｈｄ（ｆ）を算出する。この差分周波数応答Ｈｄ（ｆ）（＝Ｈ（ｆ）−ＨＵ（ｆ））が、調整すべき周波数応答となる。
【００５５】
差分部１０６で算出した差分周波数応答Ｈｄ（ｆ）と、周波数応答関数測定部１０５によって算出された周波数応答関数Ｈ（ｆ）は、被調整システムを制御するためのパラメトリックイコライザ２０の各パラメータを計算するためのウェーブレット（wavelet）多重解像度分析部１０８に入力される。
【００５６】
パラメトリックイコライザ２０は、図３に示す例のように、中心周波数ｆｃ、帯域幅（Ｑ）及びゲインからなるパラメータを制御可能なフィルタを有する。
【００５７】
ウェーブレット（wavelet）多重解像度分析部１０８は、入力された差分周波数応答Ｈｄ（ｆ）が離散化された周波数の関数であるので、連続ウェーブレット（wavelet）分析を行うため、スプライン関数などにより差分周波数応答Ｈｄ（ｆ）について内挿処理を施し、連続関数化した上で、中心周波数ｆｃとスケールパラメータＳＰとの対応関係を表す多重解像度の分析値を算出する。この分析のために用いるマザーウェーブレットには、図３に例示するように、関数の概形がパラメトリックイコライザ２０におけるフィルタによる周波数領域のフィルタ形に類似するメキシカンハット関数を選択して用いるのが好適である。
【００５８】
ウェーブレット（wavelet）多重解像度分析部１０８による分析例を図４に示す。図４では、Ｘ軸は中心周波数ｆｃ、Ｙ軸はスケールパラメータＳＰ、Ｚ軸は分析値となる。
【００５９】
イコライザパラメータ計算部１１１は、有効帯域判定部１０９から得られる有効周波数帯域判定信号（０又は１のラベル）を用いて、ウェーブレット（wavelet）多重解像度分析部１０８から得られる分析結果の周波数軸に対してラベル付けを行い調整すべき中心周波数ｆｃを決定するとともに、この中心周波数に対応するスケールパラメータＳＰ（帯域幅Ｑの算出に用いるパラメータ）を分析値の値に従って決定し、決定した中心周波数ｆｃ及びスケールパラメータＳＰを用いて帯域幅Ｑを算出するとともに、決定した中心周波数ｆｃと周波数応答関数測定部１０５から得られる周波数応答関数測定部１０５によって算出された周波数応答関数Ｈ（ｆ）とを用いてゲインＧを算出する。つまり、理想周波数帯域データベース１０７から得られる周波数特性ＨＵ（ｆ）と、周波数応答関数測定部１０５によって算出された周波数応答関数Ｈ（ｆ）との間に差がある場合、図４に示したように分析値は、周波数−スケールパラメータ平面の中で一つ以上のピーク又はディップを持つ。それぞれのピークとディップは、測定信号と理想特性との間の差を表す周波数成分を意味する。従って、これらのピーク又はディップを小さくする（周波数特性を平坦化する）ことで、被調整システムの周波数特性を調整することができるようになる。
【００６０】
より具体的には、ピーク又はディップに対応する周波数が、パラメトリックイコライザ２０における調整すべきパラメータの一つである中心周波数ｆｃとなる。ただし、有効帯域判定部１０９で有効として判断された中心周波数ｆｃが調整対象となることに留意する。
【００６１】
また、パラメトリックイコライザ２０における調整すべきパラメータの一つである帯域幅Ｑは、ピーク又はディップに対応するスケールパラメータＳＰから、式（５）で与えられる。
【００６２】
【数５】
2012028874.tif 000007
【００６３】
ここで、中心周波数ｆｃは、ピーク又はディップの周波数であり、半値幅ｗｈは、図３に示すマザーウェーブレットが０．５となる値の幅である。
【００６４】
パラメトリックイコライザ２０における調整すべきパラメータの一つであるゲインＧは、ウェーブレット（wavelet）多重解像度分析部１０８での解析結果で得られた中心周波数ｆｃに対応する周波数応答関数Ｈ（ｆ）の値Ａｆｃと周波数応答関数Ｈ（ｆ）の平均値Ａｍとを用いて、式（６）により得られる。
【００６５】
【数６】
2012028874.tif 000008
【００６６】
イコライザパラメータ計算部１１１は、得られた各パラメータ（ｆｃ，Ｑ，Ｇ）を、被調整システムに含まれるパラメトリックイコライザ２０へ適用し、この測定・分析を繰り返すことにより、被調整システムの周波数特性を所望の特性に調整することができる。
【００６７】
また、本実施例の再生周波数分析装置１０を用いて、自動で被調整システムの周波数を調整することが可能である。この場合、差分周波数応答Ｈｄ（ｆ）の値の絶対値が有効周波数帯域における全ての周波数にわたり所定値以下となった場合に、所望の特性になったと判断することで自動調整することができる。
【００６８】
尚、自動調整の場合には、一つの中心周波数を調整することで、他の中心周波数の振幅特性に影響を及ぼすことがあるので、一度の計測で得られたウェーブレット（wavelet）多重解像度分析部１０８の分析結果のピーク又はディップの絶対値の最大のものから順に所定個数分の調整すべき中心周波数を決定し、パラメトリックイコライザ２０における調整すべきパラメータを決定する。
【００６９】
また、本実施例の再生周波数分析装置１０を用いて、調整を繰り返す際に、パラメトリックイコライザの得られた各パラメータ（ｆｃ，Ｑ，Ｇ）のうち、中心周波数ｆｃと帯域幅Ｑが非常に近い値が得られる場合、上記ピーク又はディップの平坦化を得るべく効果的に収束させるために、前回の調整すべき値である中心周波数ｆｃと帯域幅Ｑを採用して修正を繰り返すようにするのが好適である。
【００７０】
上述の実施例については代表的な例として説明したが、本発明の趣旨及び範囲内で、多くの変形及び置換することができることは当業者に明らかである。例えば、理想周波数帯域データベース１０７から得られる周波数特性ＨＵ（ｆ）を複数種類用意しておき、それぞれについて最適化計算を行って得られるパラメータの平均を算出して最終決定してもよい。また、測定信号には、必ずしもピンクノイズでなくともよく、ホワイトノイズやピンクノイズ及びホワイトノイズの合成のノイズを用いることもできる。従って、本発明は、上述の実施例によって制限するものと解するべきではなく、特許請求の範囲によってのみ制限される。
【産業上の利用可能性】
【００７１】
本発明によれば、聴取位置で測定された再生周波数応答（振幅周波数特性）とその振幅周波数特性を多重解像度解析した結果から、パラメトリックイコライザの各パラメータ（中心周波数、帯域幅、ゲイン）を高精度で算出することが可能となり、音響再生システムの周波数特性の最適化を図る用途に有用である。
【符号の説明】
【００７２】
１０  再生周波数分析装置
２０  パラメトリックイコライザ
３０  増幅器
４０  スピーカ
５０  マイク
６０  マイクアンプ
１０１  測定信号発生部
１０２  遅延部
１０３  相互相関関数測定部
１０４  コヒーレンス関数測定部
１０５  周波数応答関数測定部
１０６  差分部
１０７  理想周波数特性データベース
１０８  ウェーブレット（wavelet）多重解像度分析部
１０９  有効帯域判定部
１１０  イコライザパラメータ計算部
【図１】
2012028874.tif 000009
【図２】
2012028874.tif 000010
【図３】
2012028874.tif 000011
【図４】
2012028874.tif 000012