JP3634490B2

JP3634490B2 - Impulse response measurement device in sound field

Info

Publication number: JP3634490B2
Application number: JP06492696A
Authority: JP
Inventors: 洋幸大久保; 隆司西
Original assignee: Japan Broadcasting Corp
Current assignee: Japan Broadcasting Corp
Priority date: 1996-03-21
Filing date: 1996-03-21
Publication date: 2005-03-30
Anticipated expiration: 2016-03-21
Also published as: JPH09258749A

Description

【０００１】
【発明の属する技術分野】
本発明は、音場におけるインパルス応答を測定する装置に係わり、特に高精度にインパルス応答を自動測定することのできる装置に関する。
【０００２】
［発明の概要］
本発明は、室内の音響伝達特性を表すインパルス応答をクロススペクトル法によって求める場合、測定されるインパルス応答のノイズ成分を考慮することにより高精度なインパルス応答を自動的に求める装置に関するものである。
【０００３】
従来のインパルス応答を求める自動測定法は広い測定周波数帯域に渡って必要十分な信号対雑音比（以下Ｓ／Ｎ比と称する）が取れているかの判断が難しく自動測定が困難であった。
【０００４】
本発明では、算出されるインパルス応答を平均化していく段階で前後のインパルス応答波形についての差のノルムを基に応答信号のノイズ成分を、監視する評価値として採用することにより高精度なインパルス応答を自動的に求めることを可能にする装置である。
【０００５】
また、他の発明では、測定音源として使用される楽音信号と、得られる応答信号とをいくつかの周波数帯域に分割し、帯域毎にコヒーレンス関数の値を基に十分なＳ／Ｎ比を持った応答かどうかを判断し、不十分なＳ／Ｎ比を持つ応答については平均化の操作を行わないことによって有効な平均化操作を行い、高精度なインパルス応答を自動的に短時間で求めることを可能にする装置である。
【０００６】
【従来の技術】
従来より、室内の音響伝達特性を測定する場合には、クロススペクトル法を用いたものが知られている（参考文献１：城戸健一著『デジタル信号処理入門』、Ｐ４７−Ｐ４９、丸善（１９８５））。
【０００７】
クロススペクトル法によってインパルス応答を求める場合、測定結果として得られる応答のうち真のインパルス応答を信号Ｓとし、その他の成分を雑音Ｎとしたとき、Ｓ／Ｎ比は次の式１のように表される。
【０００８】
【数１】

【０００９】
従来のインパルス応答測定において、広帯域に渡って十分なＳ／Ｎ比を持つインパルス応答を得るには、予め必要と思われる時間以上の長時間の測定を行い、得られた応答波形に平均化処理を行った上で、インパルス応答のＳ／Ｎ比を目で確認する必要があった。
【００１０】
しかし、実際の音響測定では、真のインパルス応答が不明であるために、Ｓ／Ｎ比の算出が不可能であり、このＳ／Ｎ比を用いた測定の自動化も困難であった。
【００１１】
ところで、インパルス応答を求めるための測定用の信号源に広帯域の測定用信号を用いた場合は、精度良く室内の音響伝達特性が測定できる反面、ホールやスタジオ等の室内の音響伝達特性を知るには聴衆が入るなど、実際の運用状態における特性を知ることが最も重要となる。
【００１２】
運用状態の特性を知るために室内に聴衆がいる状態の特性が必要な場合、従来のインパルス応答の測定においは、特殊な測定信号を使用している。このため、聴衆が聞き苦しかったり、聴覚に悪い影響を与えてしまう等の心配があり、運用状態での測定は困難であった。
【００１３】
そこで、測定用の信号源として聴衆の聞き慣れた楽音を用いる方法が採用されている（参考文献２：西隆司、古川宣一 ”楽音信号を用いた室のインパルス応答推定について” 日本音響学会講演論文集（１９９８／１０））。
【００１４】
しかし、楽音を測定音源として使用した場合に得られる応答信号の周波数成分は定常的に広帯域なものではなく、時々刻々に変化している。このため、Ｓ／Ｎ比が広帯域に渡って十分な応答を得るためには、やはり予め必要と思われる時間以上の長時間の測定を行わねばならない。また、得られた応答波形に多くの平均化処理を行った上で、インパルス応答のＳ／Ｎ比を目視により確認しなければならなかった。
【００１５】
【発明が解決しようとする課題】
このように、上述した従来のインパルス応答測定においては、真のインパルス応答が不明であることから、Ｓ／Ｎ比の算出が不可能であるという問題がある。一方、測定信号源として楽音を用いる従来例では、楽音信号の性質上、得られる応答信号の周波数特性は定常的でなく時々刻々に変化するため、広い測定周波数帯域に渡って必要十分なＳ／Ｎ比が取れているかの判断が難しく精度の高いインパルス応答を得るのは困難であるという問題があった。
【００１６】
本発明は上記の事情に鑑みてなされたものであり、その目的は、必要以上に時間をかけることなく、かつ聴衆に苦痛を与えること無く自動的に十分なＳ／Ｎ比を持つインパルス応答を得ることを可能にし、高精度な音響伝達特性を得ることのできる音場におけるインパルス応答測定装置を提供することにある。
【００１７】
【課題を解決するための手段】
上記の目的を達成するために請求項１に記載の発明は、音源信号と、当該音源信号を音場に放射・収音して得た応答信号を入力するとともに、第１〜第４処理を行うインパルス応答測定装置であって、第１処理は、同一時刻における音源信号と応答信号を一定長の時間窓で切り出し、第２処理は、第１処理で切り出した音源信号と応答信号からクロススペクトル法を用いてインパルス応答を算出し、第３処理は、第２処理で算出したインパルス応答とそれまでに算出したインパルス応答の平均である第１平均インパルス応答と、それまでに算出したインパルス応答の平均である第２平均インパルス応答の差のノルムに基づく評価値を算出し、第４処理は、第３処理で求めた評価値がしきい値を下回った場合は、第１平均インパルス応答を真のインパルス応答として出力し、そうでない場合は、時刻をシフトした時間窓にて第１処理に戻ることを特徴としている。
【００１８】
請求項１の発明では、逐次平均化処理を行う中で、平均化操作を行う前後のインパルス応答の差のノルムを基にした評価値を導入している。一般的に平均化が行われるに従いインパルス応答のノイズ成分が減少し、それと共に評価値は０％に漸近していく。評価値が十分小さくなった時点で測定を終了するという判断を行う。これにより、必要以上に時間をかけることなく、自動的に十分なＳ／Ｎ比を持つインパルス応答を得ることが可能となる。
【００１９】
請求項２の発明は、楽音信号と、当該楽音信号を音場に放射・収音して得た応答信号を入力するとともに、第１〜第５処理を行うインパルス応答測定装置であって、第１処理は、同一時刻における楽音信号と応答信号を一定長の時間窓で切り出して、複数の帯域に分割し、第２処理は、各帯域の楽音信号と応答信号に基づいて、帯域毎にコヒーレンス関数値を算出し、コヒーレンス関数値がしきい値を上回る帯域を有効とし、第３処理は、第２処理で有効とした帯域の楽音信号と応答信号からクロススペクトル法を用いて帯域別インパルス応答を算出し、第４処理は、第３処理で算出した帯域別インパルス応答を帯域毎に平均化して帯域別平均インパルス応答を算出し、第５処理は、全ての帯域において平均化した回数が設定された回数に達した場合は、各帯域別平均インパルス応答を合成して全帯域インパルス応答を生成し、そうでない場合には、時刻をシフトした時間窓にて第１処理に戻ることを特徴としている。
【００２０】
請求項２の発明では、音源となる楽音信号とその応答信号とから、ある窓長で波形を切り出した後、周波数軸上で狭帯域に分割し、帯域毎に音源信号と応答信号とのコヒーレンス関数値を求め、このコヒーレンス関数値に基づいてＳ／Ｎ比の評価を行う。楽音信号は時間経過と共に周波数特性が変動しているために、切り出した時刻と帯域によってＳ／Ｎ比が変動する。
【００２１】
Ｓ／Ｎ比とコヒーレンス関数とは１対１の対応関係にあるため、Ｓ／Ｎ比の変動に伴ってコヒーレンス関数も増減する。コヒーレンス関数の値が予め設定したしきい値よりも大きい場合には、Ｓ／Ｎ比が十分とれていると考えられるので、クロススペクトル法の処理を続行して平均化処理を行う。一方、コヒーレンス関数の値が前記しきい値よりも小さい場合には、Ｓ／Ｎ比が十分でないと考えられるので、その後の処理は行わないという判断をする。従って、Ｓ／Ｎ比が小さい帯域の伝達特性は平均化操作に採用しないので、収束判定が容易になり、少ない平均回数で処理を終了することができ、測定時間の短縮化、測定精度の向上が可能となる。
【００２２】
請求項３の発明は、請求項２に記載のインパルス応答測定装置において、第２処理は、ある帯域のコヒーレンス関数値が、設定期間しきい値を下回る場合には、楽音信号を切り替える旨を促す警告信号を生成することを特徴としている。
【００２３】
請求項３の発明では、楽音信号が全ての周波数成分を含んでおらず、ある帯域についてコヒーレンス関数値が長時間に渡ってしきい値を越えない場合には、警告を発して楽音切り替えを促すようにしているので、信号処理上は無駄の楽音を長時間に渡って放射し続けることが防止され、測定時間の短縮化、測定精度の向上が可能となる。
【００２４】
【発明の実施の形態】
＜第１の実施の形態＞
《第１の実施の形態の構成》
図１は本発明の第１の実施の形態を示すブロック図である。
【００２５】
同図に示すように、このインパルス応答測定装置は、信号源１Ａの音源信号Ｓａをアンプ２により増幅して室内３に放射する一方、この放射された音源信号Ｓａを収音し応答信号ＳｂとしてＡ／Ｄ変換器４によりディジタル信号に変換した後、信号処理部５Ａに供給するように構成されている。
【００２６】
この場合、信号源１Ａはテープ等で構成され、この信号源１Ａの音源信号Ｓａはアンプ２を介して室内３に設けられたスピーカ６から出力され、同じく室内３に設けられたマイクロホン７によって収音される。室内３は、ホールやスタジオ、あるいは一般家庭等、測定対象となる種々の環境が考えられる。
【００２７】
この実施の形態における信号処理部５Ａは、図２に示すように、２つのゲート回路６Ａ，７Ａと、クロススペクトル法処理部８Ａと、平均化処理部９Ａと、計算結果格納部１０Ａと、評価値演算部１１と、比較演算部１２とを備えている。
【００２８】
ゲート回路６Ａには信号源１Ａから出力された音源信号Ｓａが供給され、また、ゲート回路７Ａにはマイクロホン７で収音された応答信号Ｓｂが供給されている。各ゲート回路６Ａ，７Ａにおける信号の取り込み開始時刻ｔｓｔは、ｔ１，ｔ２，ｔ３，…，ｔｎに設定されており、その取り込み時間はＴに設定されている。すなわち、このゲート回路６Ａ，７Ａでは、応答波形は設定された時間長単位で切り出される。
【００２９】
クロススペクトル法処理部８Ａは、ゲート回路６Ａ，７Ａから取り込まれた信号に対してクロススペクトル法によりインパルス応答を算出する。
【００３０】
平均化処理部９Ａは、クロススペクトル法処理部８Ａで算出されたインパルス応答を平均化して計算結果格納部１０Ａに供給する。
【００３１】
計算結果格納部１０Ａは、平均化処理部９Ａで算出されたインパルス応答の平均化値を格納するとともに、後述する比較演算部１２から設定終了トリガ信号が出力された時、格納されているインパルス応答の平均化値を真のインパルス応答として出力する。
【００３２】
評価値演算部１１は、平均化処理部９Ａによってインパルス応答の平均化演算がされる度に、インパルス応答の差のノルムの評価値（ＤＲ値）を後述する第２式に基づいて算出する。
【００３３】
比較演算部１２は、評価値演算部１１で算出されたＤＲ値を予め設定されたしきい値と比較し、ＤＲ値がしきい値より小さい場合には、設定終了トリガ信号を出力する。この場合、しきい値としては、“０”に近い値が設定される。
【００３４】
《第１の実施の形態の作用》
次に本実施の形態の作用を図３のフローチャートを参照して説明する。
【００３５】
信号源１Ａの音源信号Ｓａはアンプ２で増幅された後、スピーカ６に供給されると共に、Ａ／Ｄ変換器４に供給されてディジタル化される。スピーカ６から室内に放射された音源信号Ｓａはマイクロホン７で収音され応答信号ＳｂとなりＡ／Ｄ変換器４に供給されてディジタル化される。
【００３６】
前記ディジタル化された音源信号Ｓａは信号処理部５Ａのゲート回路６Ａに供給され、また、前記ディジタル化された応答信号Ｓｂは信号処理部５Ａのゲート回路７Ａに供給される。
【００３７】
ゲート回路６Ａ、及び７Ａでは入力された音源信号Ｓａ、及び応答信号Ｓｂに対して時刻ｔ１，ｔ２，ｔ３，…，ｔｎにおける一定区間（取り込み時間Ｔ）の波形が切り出される。クロススペクトル法処理部８Ａでは最初の時刻ｔ１における（すなわち、１つめの）インパルス応答が算出される（ステップＳＴ１）。次いで、切り出された部分が時間Ｔだけシフトされ、時刻ｔ２におけるインパルス応答が算出される（ステップＳＴ２）。
【００３８】
平均化処理部では、時刻ｔ１と時刻ｔ２で得られたインパルス応答の平均が演算される。計算結果は計算結果格納部１０Ａに格納される（ステップＳＴ３）。
【００３９】
次いで、評価値演算部１１では、平均化を実行する都度、次の式２に示すようなインパルス応答の差のノルムの評価値（ＤＲ値）が算出される（ステップＳＴ４）。
【００４０】
【数２】

ここで、ＤＲ_ｋ：平均化をｋ回行ったときのインパルス応答の差のノルム評価値
Ｅ_ｋ（ｎ）：平均化をｋ回行ったときのインパルス応答波形をＮ個のサンプルに分割したときのｎサンプル目の値
【００４１】
比較演算部１２では、算出されたＤＲ値が設定値（しきい値）と比較して小さいか否かが判定される（ステップＳＴ５）。平均化処理部９Ａで実行される平均化が十分行われていない段階では、得られるインパルス応答にはノイズ成分が大量に含まれており、得られる平均化前後のＤＲ値は大きな値となる。この場合、処理はステップＳＴ２に戻り、次の時刻ｔ３におけるインパルス応答が算出され、次いで、ｔ１，ｔ２，ｔ３で求められたインパルス応答の平均が算出される。そして、平均化前後のＤＲ値が再び算出され、このＤＲ値が設定値より小さいか否かが判定される（ステップＳＴ２〜ＳＴ５）。
【００４２】
このように、時刻ｔ１，ｔ２，ｔ３，…，ｔｎにおけるインパルス応答が逐次算出され（ステップＳＴ２）、得られたインパルス応答が逐次平均化され（ステップＳＴ３）、そのＤＲ値が算出される（ステップＳＴ４）。
【００４３】
インパルス応答の平均化が十分行われると差のノルムは０へ漸近する。これに伴って、ＤＲ値も０（％）に漸近していく。そして、ＤＲ値が設定値を下回ったときに比較演算部１２は測定終了トリガ信号を出力する。計算結果格納部１０Ａは、この測定終了トリガ信号を入力すると、このとき格納されている最新のインパルス応答平均値を真のインパルス応答として出力する。
【００４４】
《第１の実施の形態の効果》
このように、この実施の形態によれば、クロススペクトル法によって測定される室内の音響伝達特性を表すインパルス応答のＳ／Ｎ比を大きくするために、必要以上に平均化回数を増やすことなく、ノイズ成分をＤＲ値を基にして定量的に評価できる。従って、平均化回数を減らすことができるため、測定時間の短縮が可能となるばかりでなく、時間経過に伴って変化してしまう伝達特性に対して追随性の良い測定が可能となる。
【００４５】
＜第２の実施の形態＞
《第２の実施の形態の構成》
図４は本発明の第２の実施の形態を示すブロック図である。
【００４６】
同図に示すように、このインパルス応答測定装置は、信号源１Ｂの音源信号Ｓｃをアンプ２により増幅して室内３に放射する一方、この放射された音源信号Ｓｃを収音し応答信号ＳｂとしてＡ／Ｄ変換器４によりディジタル信号に変換した後、信号処理部５Ｂに供給するように構成されている。
【００４７】
この場合、信号源１Ｂは楽音を収録したテープ等で構成される。また、この音源として使用する楽音は高い周波数成分まで含んでいる楽音が望ましい。例えば、大編成のオーケストラ、電気楽器による音楽、打楽器中心の音楽など、複数の楽音テープが用意されている。この信号源１Ｂの音源信号Ｓｃはアンプ２を介して室内３に設けられたスピーカ６から出力され、同じく室内３に設けられたマイクロホン７によって収音される。室内３は、ホールやスタジオ、あるい一般家庭等、測定対象となる種々の環境が考えられる。
【００４８】
この実施の形態における信号処理部５Ｂは、図５に示すように、その機能上、２つのゲート回路６Ｂ，７Ｂと、各ゲート回路６Ｂ，７Ｂ別に設けられた帯域分割フィルタ部１３，１４と、コヒーレンス関数値演算部１５と、比較演算部１６と、クロススペクトル法処理部８Ｂと、平均化処理部９Ｂと、計算結果格納部１０Ｂとを備えている。
【００４９】
ゲート回路６Ｂには、信号源１Ｂから出力された音源信号Ｓｃが供給され、また、ゲート回路７Ｂにはマイクロホン７で収音された応答信号Ｓｄが供給されている。各ゲート回路６Ｂ，７Ｂにおける信号の取り込み開始時刻ｔｓｔは、ｔ１，ｔ２，ｔ３，…，ｔｎであり、取り込み時間はＴに設定されている。
【００５０】
帯域分割フィルタ部１３は、供給された音源信号Ｓｃを１／１オクターブ帯域に分割するフィルタを備えており、分割された各帯域毎の音源信号Ｓｃをコヒーレンス関数値演算部１５に供給する。
【００５１】
帯域分割フィルタ部１４は、供給された応答信号Ｓｄを１／１オクターブ帯域に分割するフィルタを備えており、分割された各帯域毎の応答信号Ｓｄをコヒーレンス関数値演算部１５に供給する。
【００５２】
コヒーレンス関数値演算部１５は、各帯域毎の音源信号Ｓｃと応答信号Ｓｄとを入力し帯域毎のコヒーレンス関数値を算出する。
【００５３】
比較演算部１６は、コヒーレンス関数値演算部１５により求められたコヒーレンス関数値が予め設定されたしきい値を下回った場合に、クロススペクトル処理を行わないという制御信号と、音源切り替え警告信号を生成する。制御信号はクロススペクトル法処理部８Ｂに供給されている。
【００５４】
クロススペクトル法処理部８Ｂは、帯域分割フィルタ部１３，１４を介して帯域分割された各信号成分に対してクロススペクトル法によりインパルス応答を算出する処理を、比較演算部１６から前記制御信号が供給されるまで実行する。
【００５５】
平均化処理部９Ｂは、前述した時刻ｔ１、ｔ２，…，ｔｎをシフトしながら、逐次インパルス応答を平均化する演算を実行する。
【００５６】
計算結果格納部１０Ｂは、平均化処理部９Ｂで演算されたインパルス応答の平均化値を格納するとともに、後述する比較演算部１２から設定終了トリガ信号が出力された時、格納されている各帯域インパルス応答の平均化値を合成し、広帯域のインパルス応答を算出する。
【００５７】
《第２の実施の形態の作用》
次に本実施の形態の作用を図６のフローチャートを参照して説明する。なお、図６のフローチャートにおいては、原則として各算出処理や各判定処理は帯域別に並行して実行されているものとする。
【００５８】
第１の実施の形態で説明したと同様、信号源１Ｂの音源信号（楽音信号）ＳＣはアンプ２で増幅された後、スピーカ６に供給されると共に、Ａ／Ｄ変換器４に供給されてディジタル化される。スピーカ６から室内に放射された音源信号Ｓｃはマイクロホン７で収音され応答信号ＳｄとなりＡ／Ｄ変換器４に供給されてディジタル化される。
【００５９】
前記ディジタル化された音源信号Ｓｃは信号処理部５Ｂのゲート回路６Ｂに供給され、また、前記ディジタル化された応答信号Ｓｄは信号処理部５Ｂのゲート回路７Ｂに供給される。
【００６０】
ゲート回路６Ｂ、及び７Ｂでは入力された音源信号Ｓｃ、及び応答信号Ｓｄに対して時刻ｔ１，ｔ２，ｔ３，…，ｔｎにおける一定時間窓（取り込み時間Ｔ）の波形が切り出される。次いで、帯域分割フィルタ部１３、及び１４では、それぞれ切り出された波形が１／１オクターブ帯域に分割され、それぞれコヒーレンス関数値演算部１５に供給される（ステップＳＴ１１）。
【００６１】
コヒーレンス関数演算部１５では、帯域毎のコヒーレンス関数値が次の式３に基づいて算出される（ステップＳＴ１２）。
【００６２】
【数３】

【００６３】
比較演算部１６では、算出されたコヒーレンス関数値が予め設定されたしきい値と比較され、コヒーレンス関数値がしきい値よりも小さい場合には、その後の処理をスキップさせる制御信号をクロススペクトル法処理部８Ｂに供給するとともに、時間窓をシフトさせ、次の時刻ｔ２における切り出された波形に対する前記一連の処理を実行する（ステップＳＴ１５）。また、コヒーレンス関数値が設定値よりも小さい場合において、音源として用いられている楽音信号が全ての周波数成分を含んでおらず、ある帯域についてコヒーレンス関数値が長時間に渡ってしきい値を越えない場合には（ステップＳＴ１４Ｙｅｓ）、楽音信号を別の楽音信号に切り替える警告信号を出力する。これにより、信号源１Ｂの音源信号は別の楽音信号に切り替わった後、ステップＳＴ１１からの処理を開始する（ステップＳＴ１６）。
【００６４】
コヒーレンス関数値がしきい値を越えた帯域のデータは有効なものとし、クロススペクトル法処理部では、クロススペクトル法によってインパルス応答が算出される（ステップＳＴ１７）。
【００６５】
そして、平均化処理部９Ｂでは、時間窓をシフト（開始時間ｔ３，…，ｔｎ、切り出し時間Ｔ）させつつ逐次求められたインパルス応答が平均化されて計算結果格納部１０Ｂに格納される（ステップＳＴ１８）。
【００６６】
平均化処理部９Ｂでは、平均回数が設定された回数よりも大きいか否かが判定され、設定された平均回数を終了すると（ステップＳＴ１９ＹＥＳ）、測定終了トリガ信号が出力されて、平均化処理が終了する。計算結果格納部１０Ｂは、測定終了トリガ信号が入力されると、各帯域のインパルス応答を合成し、全帯域のインパルス応答を算出して測定を終了する（ステップＳＴ２０）。
【００６７】
《第２の実施の形態の効果》
このようにこの実施の形態によれば、実際に室内で使用される状態での音響伝達特性を測定できることにより、聴衆がいる等の現実により近い状態の室内の音響伝達特性を精度良く知ることができる。ここで得られる知見を発展させるとこにより、人が居ない状態から人が居る状態を精度良く予測することが可能となり、設計段階で完成後の運用状態の音響伝達特性を予測することができ、音響設計の精度も向上する。
【００６８】
また、ホール等の特殊な施設でなく一般家庭のおいても音響機器等のこの実施の形態を適用して、楽音による測定を行うことにより、室のユーザー（聴取者）が意識することなく測定をすることができる。これにより、室内、音響機器伝達特性を知ることができれば、その特性を等化することも可能となり、仮に聴取環境が異なっても音質、定位を初めとする音声制作上の細かい効果を同一条件とすることができる。
【００６９】
＜変形例＞
なお、上述した第２の実施の形態では、帯域分割フィルタ部１３，１４により、供給された音源信号Ｓｃ、応答信号Ｓｄを狭帯域の信号に分割するようにしたが、コヒーレンス関数値演算部１５において帯域別の関数値を算出するようにしても良い。このように構成することにより、帯域分割フィルタ部１３，１４を省略でき、装置構成が簡素化される。
【００７０】
【発明の効果】
以上説明したように請求項１の発明によれば、必要以上に時間をかけることなく、自動的に十分なＳ／Ｎ比を持つインパルス応答を得ることが可能となる。
【００７１】
請求項２の発明によれば、Ｓ／Ｎ比が小さい帯域の伝達特性は平均化操作に採用しないので、収束判定が容易になり、少ない平均回数で処理を終了することができ、測定時間の短縮化、測定精度の向上が可能となる。
【００７２】
請求項３の発明によれば、楽音信号が全ての周波数成分を含んでおらず、ある帯域についてコヒーレンス関数値が長時間に渡ってしきい値を越えない場合には、警告を発して楽音切り替えを促すようにしているので、信号処理上は無駄の楽音を長時間に渡って放射し続けることが防止され、測定時間の短縮化、測定精度の向上が可能となる。
【図面の簡単な説明】
【図１】本発明の第１の実施の形態を示す構成図である。
【図２】図１に示す第１の実施の形態における信号処理部の構成を示すブロック図である。
【図３】第１の実施の形態における処理手順を説明するフローチャートである。
【図４】本発明の第２の実施の形態を示す構成図である。
【図５】図４に示す第２の実施の形態における信号処理部の構成を示すブロック図である。
【図６】第２の実施の形態における処理手順を説明するフローチャートである。
【符号の説明】
１Ａ，１Ｂ信号源
２アンプ
３室内
４Ａ／Ｄ変換器
５Ａ，５Ｂ信号処理部
６Ａ，６Ｂ，７Ａ，７Ｂゲート回路
８Ａ，８Ｂクロススペクトル法処理部
９Ａ，９Ｂ平均化処理部
１０Ａ，１０Ｂ計算結果格納部
１１評価値演算部
１２，１６比較演算部
１３，１４帯域分割フィルタ部
１５コヒーレンス関数値演算部[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an apparatus for measuring an impulse response in a sound field, and more particularly to an apparatus capable of automatically measuring an impulse response with high accuracy.
[0002]
[Summary of Invention]
The present invention relates to an apparatus for automatically obtaining a highly accurate impulse response by considering a noise component of a measured impulse response when an impulse response representing an indoor acoustic transfer characteristic is obtained by a cross spectrum method.
[0003]
In the conventional automatic measurement method for obtaining the impulse response, it is difficult to determine whether a necessary and sufficient signal-to-noise ratio (hereinafter referred to as S / N ratio) is obtained over a wide measurement frequency band, and automatic measurement is difficult.
[0004]
In the present invention, the noise component of the response signal is adopted as an evaluation value to be monitored based on the norm of the difference between the preceding and succeeding impulse response waveforms at the stage of averaging the calculated impulse responses, thereby providing a highly accurate impulse response. It is a device that makes it possible to automatically obtain.
[0005]
In another invention, the musical sound signal used as a measurement sound source and the obtained response signal are divided into several frequency bands, and each band has a sufficient S / N ratio based on the value of the coherence function. If the response has an insufficient S / N ratio, an effective averaging operation is performed by not performing the averaging operation, and a highly accurate impulse response is automatically obtained in a short time. It is a device that makes it possible.
[0006]
[Prior art]
Conventionally, in order to measure the acoustic transfer characteristics in a room, one using a cross spectrum method is known (Reference 1: Kenichi Kido, “Introduction to Digital Signal Processing”, P47-P49, Maruzen (1985). ).
[0007]
When the impulse response is obtained by the cross spectrum method, the S / N ratio is expressed by the following equation 1 when the true impulse response of the response obtained as a measurement result is the signal S and the other component is the noise N. Is done.
[0008]
[Expression 1]

[0009]
In the conventional impulse response measurement, in order to obtain an impulse response with a sufficient S / N ratio over a wide band, measurement is performed for a long time that is longer than necessary, and the obtained response waveform is averaged. It was necessary to visually confirm the S / N ratio of the impulse response.
[0010]
However, in actual acoustic measurement, since the true impulse response is unknown, it is impossible to calculate the S / N ratio, and it is difficult to automate the measurement using this S / N ratio.
[0011]
By the way, when a wide-band measurement signal is used as a measurement signal source for obtaining an impulse response, it is possible to accurately measure the acoustic transfer characteristics in a room, but to know the acoustic transfer characteristics in a room such as a hall or a studio. It is most important to know the characteristics in the actual operational state, such as the audience entering.
[0012]
When the characteristic of the state where the audience is in the room is necessary to know the characteristic of the operation state, a special measurement signal is used in the measurement of the conventional impulse response. For this reason, there were concerns that the audience was hard to hear or had a bad influence on hearing, and it was difficult to measure in an operational state.
[0013]
Therefore, the method of using the musical sound familiar to the audience is adopted as the signal source for the measurement (Reference 2: Takashi Nishi, Nobuichi Furukawa "About Impulse Response Estimation in the Room Using Musical Signal" Lecture by the Acoustical Society of Japan Papers (1998/10)).
[0014]
However, the frequency component of the response signal obtained when the musical sound is used as a measurement sound source is not constantly a wide band, but changes every moment. For this reason, in order to obtain a sufficient response over an S / N ratio in a wide band, it is necessary to perform measurement for a long time that is longer than necessary. In addition, it was necessary to visually check the S / N ratio of the impulse response after performing many averaging processes on the obtained response waveform.
[0015]
[Problems to be solved by the invention]
Thus, in the conventional impulse response measurement described above, there is a problem that the S / N ratio cannot be calculated because the true impulse response is unknown. On the other hand, in the conventional example using the musical sound as the measurement signal source, the frequency characteristic of the response signal obtained is not steady due to the nature of the musical sound signal, and changes every moment, so that the necessary S / S is sufficient over a wide measurement frequency band. There is a problem that it is difficult to determine whether the N ratio is obtained, and it is difficult to obtain a highly accurate impulse response.
[0016]
The present invention has been made in view of the above circumstances, and its purpose is to automatically generate an impulse response having a sufficient S / N ratio without taking more time than necessary and without causing pain to the audience. It is an object of the present invention to provide an impulse response measuring apparatus in a sound field that can be obtained and can obtain a highly accurate acoustic transmission characteristic.
[0017]
[Means for Solving the Problems]
In order to achieve the above object, the invention described in claim 1 inputs a sound source signal and a response signal obtained by radiating and collecting the sound source signal in a sound field, and performs the first to fourth processes. An impulse response measuring device to perform, wherein the first process cuts out a sound source signal and a response signal at the same time in a fixed time window, and the second process cross-spectrums from the sound source signal and the response signal cut out in the first process. An impulse response is calculated using the method, and the third process includes a first average impulse response that is an average of the impulse response calculated in the second process and the impulse response calculated so far, and the impulse response calculated so far. An evaluation value based on the norm of the difference between the second average impulse responses, which is an average, is calculated. In the fourth process, if the evaluation value obtained in the third process falls below the threshold value, the first average impulse response is true. of Output as impulse response, otherwise, is characterized in that the time at shifted time window to return to the first processing.
[0018]
In the first aspect of the present invention, an evaluation value based on the norm of the difference between impulse responses before and after performing the averaging operation is introduced during the sequential averaging process. In general, the noise component of the impulse response decreases as the averaging is performed, and the evaluation value gradually approaches 0%. It is determined that the measurement is terminated when the evaluation value becomes sufficiently small. This makes it possible to automatically obtain an impulse response having a sufficient S / N ratio without taking more time than necessary.
[0019]
The invention of claim 2 is an impulse response measuring apparatus for inputting a musical sound signal and a response signal obtained by radiating and collecting the musical sound signal to a sound field, and performing first to fifth processing, In the first process, the musical tone signal and the response signal at the same time are cut out by a fixed time window and divided into a plurality of bands. The second process is a coherence for each band based on the musical tone signal and the response signal in each band. The function value is calculated and the band in which the coherence function value exceeds the threshold value is valid. The third process is the impulse response for each band using the cross spectrum method from the musical sound signal and the response signal in the band that is valid in the second process. In the fourth process, the impulse response for each band calculated in the third process is averaged for each band to calculate the average impulse response for each band. In the fifth process, the number of times of averaging in all the bands is set. Reached the number of times If the, each band-specific average impulse responses to generate a synthesized to full band impulse response, otherwise, is characterized in that the time at shifted time window to return to the first processing.
[0020]
According to the second aspect of the present invention, after a waveform is cut out with a certain window length from a musical sound signal as a sound source and its response signal, it is divided into narrow bands on the frequency axis, and the coherence between the sound source signal and the response signal for each band A function value is obtained, and the S / N ratio is evaluated based on the coherence function value. Since the frequency characteristic of the musical sound signal varies with time, the S / N ratio varies depending on the cut time and band.
[0021]
Since the S / N ratio and the coherence function have a one-to-one correspondence, the coherence function also increases or decreases as the S / N ratio varies. When the value of the coherence function is larger than a preset threshold value, it is considered that the S / N ratio is sufficient, so the cross spectrum method process is continued and the averaging process is performed. On the other hand, when the value of the coherence function is smaller than the threshold value, it is considered that the S / N ratio is not sufficient, so that the subsequent processing is not performed. Therefore, transfer characteristics in a band with a small S / N ratio are not adopted in the averaging operation, so that it is easy to determine convergence, processing can be completed with a small number of averages, shortening measurement time, and improving measurement accuracy. Is possible.
[0022]
According to a third aspect of the present invention, in the impulse response measuring apparatus according to the second aspect , the second process prompts the user to switch the musical tone signal when the coherence function value of a certain band falls below a set period threshold value. It is characterized by generating a warning signal .
[0023]
According to the third aspect of the present invention, when the tone signal does not include all frequency components and the coherence function value does not exceed the threshold value for a long time in a certain band, a warning is issued to prompt tone switching. As a result, it is possible to prevent unnecessary musical sounds from being radiated for a long time in signal processing, and it is possible to shorten the measurement time and improve the measurement accuracy.
[0024]
DETAILED DESCRIPTION OF THE INVENTION
<First Embodiment>
<< Configuration of the First Embodiment >>
FIG. 1 is a block diagram showing a first embodiment of the present invention.
[0025]
As shown in the figure, this impulse response measuring apparatus amplifies the sound source signal Sa of the signal source 1A by the amplifier 2 and radiates it into the room 3, while collecting the radiated sound source signal Sa as the response signal Sb. The digital signal is converted into a digital signal by the A / D converter 4 and then supplied to the signal processing unit 5A.
[0026]
In this case, the signal source 1A is composed of a tape or the like, and the sound source signal Sa of the signal source 1A is output from the speaker 6 provided in the room 3 via the amplifier 2 and collected by the microphone 7 also provided in the room 3. Sounded. The room 3 may be various environments to be measured, such as a hall, studio, or ordinary home.
[0027]
As shown in FIG. 2, the signal processing unit 5A in this embodiment includes two gate circuits 6A and 7A, a cross spectrum method processing unit 8A, an averaging processing unit 9A, a calculation result storage unit 10A, and an evaluation. A value calculation unit 11 and a comparison calculation unit 12 are provided.
[0028]
The sound source signal Sa 1 output from the signal source 1A is supplied to the gate circuit 6A, and the response signal Sb collected by the microphone 7 is supplied to the gate circuit 7A. The signal capture start time tst in each of the gate circuits 6A and 7A is set to t1, t2, t3,..., Tn, and the capture time is set to T. That is, in the gate circuits 6A and 7A, the response waveform is cut out in units of the set time length.
[0029]
The cross spectrum method processing unit 8A calculates an impulse response by the cross spectrum method for the signals taken from the gate circuits 6A and 7A.
[0030]
The averaging processing unit 9A averages the impulse response calculated by the cross spectrum method processing unit 8A and supplies the averaged impulse response to the calculation result storage unit 10A.
[0031]
The calculation result storage unit 10A stores the average value of the impulse response calculated by the averaging processing unit 9A, and stores a stored impulse response when a setting end trigger signal is output from the comparison calculation unit 12 described later. Is output as a true impulse response.
[0032]
The evaluation value calculation unit 11 calculates an evaluation value (DR value) of a norm of a difference between impulse responses based on the second formula described later each time the averaging processing unit 9A performs an impulse response averaging calculation.
[0033]
The comparison calculation unit 12 compares the DR value calculated by the evaluation value calculation unit 11 with a preset threshold value, and outputs a setting end trigger signal when the DR value is smaller than the threshold value. In this case, a value close to “0” is set as the threshold value.
[0034]
<< Operation of First Embodiment >>
Next, the operation of the present embodiment will be described with reference to the flowchart of FIG.
[0035]
The sound source signal Sa of the signal source 1A is amplified by the amplifier 2 and then supplied to the speaker 6 and also supplied to the A / D converter 4 and digitized. The sound source signal Sa radiated indoors from the speaker 6 is picked up by the microphone 7 and becomes a response signal Sb, which is supplied to the A / D converter 4 and digitized.
[0036]
The digitized sound source signal Sa is supplied to the gate circuit 6A of the signal processing unit 5A, and the digitized response signal Sb is supplied to the gate circuit 7A of the signal processing unit 5A.
[0037]
In the gate circuits 6A and 7A, a waveform in a certain section (capture time T) at time t1, t2, t3,..., Tn is cut out from the input sound source signal Sa and response signal Sb. The cross spectrum method processing unit 8A calculates the impulse response at the first time t1 (that is, the first) (step ST1). Next, the cut out portion is shifted by time T, and the impulse response at time t2 is calculated (step ST2).
[0038]
In the averaging processing unit, an average of impulse responses obtained at time t1 and time t2 is calculated. The calculation result is stored in the calculation result storage unit 10A (step ST3).
[0039]
Next, every time averaging is performed, the evaluation value calculation unit 11 calculates an evaluation value (DR value) of a norm of a difference in impulse response as shown in the following equation 2 (step ST4).
[0040]
[Expression 2]

Here, DR _k : Norm of difference in impulse response when averaging is performed k times Evaluation value E _k (n): When the impulse response waveform when averaging is performed k times is divided into N samples The value of the nth sample of
The comparison calculation unit 12 determines whether or not the calculated DR value is smaller than the set value (threshold value) (step ST5). In the stage where the averaging performed by the averaging processing unit 9A is not sufficiently performed, the obtained impulse response contains a large amount of noise components, and the obtained DR value before and after averaging becomes a large value. In this case, the process returns to step ST2, the impulse response at the next time t3 is calculated, and then the average of the impulse responses obtained at t1, t2, and t3 is calculated. Then, the DR value before and after averaging is calculated again, and it is determined whether or not this DR value is smaller than the set value (steps ST2 to ST5).
[0042]
In this way, impulse responses at times t1, t2, t3,..., Tn are sequentially calculated (step ST2), the obtained impulse responses are sequentially averaged (step ST3), and the DR value is calculated (step ST2). ST4).
[0043]
When the impulse response is sufficiently averaged, the norm of the difference approaches 0. Along with this, the DR value gradually approaches 0 (%). When the DR value falls below the set value, the comparison calculation unit 12 outputs a measurement end trigger signal. When the calculation result storage unit 10A receives the measurement end trigger signal, the calculation result storage unit 10A outputs the latest impulse response average value stored at this time as a true impulse response.
[0044]
<< Effects of First Embodiment >>
Thus, according to this embodiment, in order to increase the S / N ratio of the impulse response representing the acoustic transfer characteristics in the room measured by the cross spectrum method, without increasing the number of averaging more than necessary, The noise component can be quantitatively evaluated based on the DR value. Therefore, since the number of times of averaging can be reduced, not only the measurement time can be shortened, but also a measurement with good followability can be performed with respect to the transfer characteristics that change with the passage of time.
[0045]
<Second Embodiment>
<< Configuration of Second Embodiment >>
FIG. 4 is a block diagram showing a second embodiment of the present invention.
[0046]
As shown in the figure, this impulse response measuring apparatus amplifies the sound source signal Sc of the signal source 1B by the amplifier 2 and radiates it into the room 3, while collecting the radiated sound source signal Sc as the response signal Sb. The digital signal is converted into a digital signal by the A / D converter 4 and then supplied to the signal processing unit 5B.
[0047]
In this case, the signal source 1B is composed of a tape or the like on which music is recorded. In addition, it is desirable that the musical sound used as the sound source includes a high frequency component. For example, a plurality of musical sound tapes are prepared, such as large orchestras, music using electric instruments, and music mainly percussion instruments. The sound source signal Sc of the signal source 1B is output from the speaker 6 provided in the room 3 via the amplifier 2, and is picked up by the microphone 7 also provided in the room 3. The room 3 can be various environments to be measured, such as a hall, studio, or ordinary home.
[0048]
As shown in FIG. 5, the signal processing unit 5B in this embodiment is functionally divided into two gate circuits 6B and 7B, and band

division filter units

13 and 14 provided for the respective gate circuits 6B and 7B. A coherence function value calculation unit 15, a comparison calculation unit 16, a cross spectrum method processing unit 8B, an averaging processing unit 9B, and a calculation result storage unit 10B are provided.
[0049]
The gate circuit 6B is supplied with the sound source signal Sc output from the signal source 1B, and the gate circuit 7B is supplied with the response signal Sd collected by the microphone 7. Signal capture start times tst in the gate circuits 6B and 7B are t1, t2, t3,..., Tn, and the capture time is set to T.
[0050]
The band division filter unit 13 includes a filter that divides the supplied sound source signal Sc into 1/1 octave bands, and supplies the divided sound source signal Sc for each band to the coherence function value calculation unit 15.
[0051]
The band division filter unit 14 includes a filter that divides the supplied response signal Sd into 1/1 octave bands, and supplies the divided response signal Sd for each band to the coherence function value calculation unit 15.
[0052]
The coherence function value calculation unit 15 receives the sound source signal Sc and the response signal Sd for each band and calculates a coherence function value for each band.
[0053]
When the coherence function value obtained by the coherence function value calculation unit 15 falls below a preset threshold value, the comparison calculation unit 16 generates a control signal not to perform cross spectrum processing and a sound source switching warning signal To do. The control signal is supplied to the cross spectrum method processing unit 8B.
[0054]
The cross spectrum method processing unit 8B supplies the control signal from the comparison operation unit 16 to calculate the impulse response by the cross spectrum method for each signal component band-divided through the band

division filter units

13 and 14. Run until
[0055]
The averaging processing unit 9B performs an operation of sequentially averaging the impulse responses while shifting the above-described times t1, t2,..., Tn.
[0056]
The calculation result storage unit 10B stores the average value of the impulse response calculated by the averaging processing unit 9B, and stores each band stored when a setting end trigger signal is output from the comparison calculation unit 12 described later. An impulse response average value is synthesized to calculate a broadband impulse response.
[0057]
<< Operation of Second Embodiment >>
Next, the operation of the present embodiment will be described with reference to the flowchart of FIG. In the flowchart of FIG. 6, in principle, each calculation process and each determination process are executed in parallel for each band.
[0058]
As described in the first embodiment, the sound source signal (musical signal) SC of the signal source 1B is amplified by the amplifier 2 and then supplied to the speaker 6 and supplied to the A / D converter 4. Digitized. The sound source signal Sc radiated into the room from the speaker 6 is picked up by the microphone 7 and becomes a response signal Sd, which is supplied to the A / D converter 4 and digitized.
[0059]
The digitized sound source signal Sc is supplied to the gate circuit 6B of the signal processing unit 5B, and the digitized response signal Sd is supplied to the gate circuit 7B of the signal processing unit 5B.
[0060]
In the gate circuits 6B and 7B, a waveform of a certain time window (capture time T) at time t1, t2, t3,..., Tn is cut out from the input sound source signal Sc and response signal Sd. Next, in the band

division filter units

13 and 14, the cut out waveforms are divided into 1/1 octave bands and supplied to the coherence function value calculation unit 15 (step ST11).
[0061]
The coherence function calculation unit 15 calculates a coherence function value for each band based on the following Equation 3 (step ST12).
[0062]
[Equation 3]

[0063]
In the comparison operation unit 16, the calculated coherence function value is compared with a preset threshold value. When the coherence function value is smaller than the threshold value, a control signal for skipping the subsequent processing is applied to the cross spectrum method. While supplying to the process part 8B, a time window is shifted and the said series of processes with respect to the cut-out waveform in the following time t2 are performed (step ST15). In addition, when the coherence function value is smaller than the set value, the tone signal used as a sound source does not contain all frequency components, and the coherence function value exceeds a threshold for a long time for a certain band. If not (step ST14 Yes), a warning signal for switching the tone signal to another tone signal is output. Thereby, after the sound source signal of the signal source 1B is switched to another musical sound signal, the processing from step ST11 is started (step ST16).
[0064]
Data in a band where the coherence function value exceeds the threshold value is valid, and the cross spectrum method processing unit calculates an impulse response by the cross spectrum method (step ST17).
[0065]
Then, the averaging processor 9B averages the impulse responses obtained sequentially while shifting the time window (start time t3,..., Tn, cutout time T) and stores the averaged impulse responses in the calculation result storage unit 10B (step S1). ST18).
[0066]
In the averaging processing unit 9B, it is determined whether or not the average number is larger than the set number. When the set average number is finished (YES in step ST19), a measurement end trigger signal is output and the averaging process is performed. finish. When the measurement end trigger signal is input, the calculation result storage unit 10B combines the impulse responses of the respective bands, calculates the impulse responses of all the bands, and ends the measurement (step ST20).
[0067]
<< Effects of Second Embodiment >>
As described above, according to this embodiment, it is possible to accurately measure the acoustic transfer characteristics in a room that is closer to reality, such as the presence of an audience, by measuring the acoustic transfer characteristics in a state where they are actually used indoors. it can. By developing the knowledge obtained here, it becomes possible to accurately predict the state where there is a person from the state where there is no person, and it is possible to predict the acoustic transfer characteristics of the operational state after completion at the design stage, The accuracy of acoustic design is also improved.
[0068]
In addition, it is possible to measure without sound by the user (listener) of the room by applying this embodiment of the audio equipment, etc., in a general household as well as in a special facility such as a hall, etc. Can do. As a result, if it is possible to know the transfer characteristics of indoor and acoustic equipment, it is also possible to equalize the characteristics, and even if the listening environment is different, fine effects on sound production including sound quality and localization are the same conditions. can do.
[0069]
<Modification>
In the second embodiment described above, the supplied sound source signal Sc and response signal Sd are divided into narrowband signals by the band

division filter units

13 and 14, but the coherence function value calculation unit 15 is used. The function value for each band may be calculated. With this configuration, the band

division filter units

13 and 14 can be omitted, and the apparatus configuration is simplified.
[0070]
【The invention's effect】
As described above, according to the first aspect of the present invention, it is possible to automatically obtain an impulse response having a sufficient S / N ratio without taking more time than necessary.
[0071]
According to the second aspect of the present invention, since the transfer characteristic of the band having a small S / N ratio is not adopted in the averaging operation, the convergence determination is facilitated, and the processing can be completed with a small average number of times. Shortening and improvement of measurement accuracy are possible.
[0072]
According to the invention of claim 3, when the musical sound signal does not include all frequency components and the coherence function value does not exceed the threshold value for a long time in a certain band, a warning is issued and the musical sound is switched. Therefore, it is possible to prevent unnecessary musical sounds from being continuously emitted for a long time in signal processing, and it is possible to shorten the measurement time and improve the measurement accuracy.
[Brief description of the drawings]
FIG. 1 is a configuration diagram showing a first embodiment of the present invention.
FIG. 2 is a block diagram showing a configuration of a signal processing unit in the first embodiment shown in FIG.
FIG. 3 is a flowchart illustrating a processing procedure according to the first embodiment.
FIG. 4 is a configuration diagram showing a second embodiment of the present invention.
FIG. 5 is a block diagram showing a configuration of a signal processing unit in the second embodiment shown in FIG. 4;
FIG. 6 is a flowchart illustrating a processing procedure in the second embodiment.
[Explanation of symbols]
1A, 1B Signal source 2 Amplifier 3 Indoor 4 A /

D converters

5A, 5B

Signal processing units

6A, 6B, 7A, 7B Gate circuits 8A, 8B Cross spectrum

method processing units

9A, 9B

Averaging processing units

10A, 10B Calculation results Storage unit 11 Evaluation value calculation unit 12, 16

Comparison calculation unit

13, 14 Band division filter unit 15 Coherence function value calculation unit

Claims

An impulse response measuring device that inputs a sound source signal and a response signal obtained by radiating and collecting the sound source signal in a sound field, and that performs first to fourth processing,
In the first process, the sound source signal and the response signal at the same time are cut out with a fixed time window,
The second process calculates an impulse response using the cross spectrum method from the sound source signal and the response signal cut out in the first process,
The third process is a difference between the first average impulse response that is the average of the impulse responses calculated in the second process and the impulse responses calculated so far, and the second average impulse response that is the average of the impulse responses calculated so far. Calculate an evaluation value based on the norm of
The fourth process outputs the first average impulse response as a true impulse response when the evaluation value obtained in the third process falls below the threshold value, otherwise, the time window is shifted in time. Returning to the first process, the impulse response measuring apparatus.

An impulse response measuring device for inputting a musical sound signal and a response signal obtained by radiating and collecting the musical sound signal in a sound field, and performing first to fifth processing,
In the first process, the musical sound signal and the response signal at the same time are cut out by a fixed length time window and divided into a plurality of bands.
The second process calculates a coherence function value for each band based on the musical tone signal and response signal of each band, and validates the band where the coherence function value exceeds the threshold value.
The third process calculates the impulse response for each band using the cross spectrum method from the musical sound signal and the response signal in the band enabled in the second process,
The fourth process calculates the average impulse response for each band by averaging the impulse response for each band calculated in the third process for each band,
In the fifth process, when the number of times of averaging in all the bands reaches the set number of times, an average impulse response for each band is synthesized to generate an all-band impulse response. An impulse response measuring apparatus which returns to the first processing in the shifted time window .

The impulse response measurement device according to claim 2, wherein the second process generates a warning signal that prompts the user to switch the musical tone signal when the coherence function value of a certain band is lower than a set period threshold value. .