JPH08184488A - 音響特性測定装置 - Google Patents
音響特性測定装置Info
- Publication number
- JPH08184488A JPH08184488A JP34058194A JP34058194A JPH08184488A JP H08184488 A JPH08184488 A JP H08184488A JP 34058194 A JP34058194 A JP 34058194A JP 34058194 A JP34058194 A JP 34058194A JP H08184488 A JPH08184488 A JP H08184488A
- Authority
- JP
- Japan
- Prior art keywords
- spectrum
- average value
- input
- signal
- input signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【構成】 測定用マイクＭＣ1 〜ＭＣ5 にそれぞれ接続
されたマイクアンプ４７1 〜４７5 からの出力を、音響
特性測定器４５からの制御信号で信号切換器４６で切換
選択して、各測定用マイク毎に、入力端子４１からの入
力信号とマイクアンプの出力とを用いて各周波数ポイン
ト毎の伝達関数の平均値を求め、各測定用マイクの周波
数ポイント毎の伝達関数の平均値を平均化して、１つの
音響特性を得る。 【効果】 再生する音声信号の調整を容易に行うことが
できる。
されたマイクアンプ４７1 〜４７5 からの出力を、音響
特性測定器４５からの制御信号で信号切換器４６で切換
選択して、各測定用マイク毎に、入力端子４１からの入
力信号とマイクアンプの出力とを用いて各周波数ポイン
ト毎の伝達関数の平均値を求め、各測定用マイクの周波
数ポイント毎の伝達関数の平均値を平均化して、１つの
音響特性を得る。 【効果】 再生する音声信号の調整を容易に行うことが
できる。
Description
【０００１】
【産業上の利用分野】本発明は、被測定物に供給する音
声信号と被測定物から出力される音声信号とから被測定
物の伝達関数を測定して音響特性を得る音響特性測定装
置に関する。
声信号と被測定物から出力される音声信号とから被測定
物の伝達関数を測定して音響特性を得る音響特性測定装
置に関する。
【０００２】
【従来の技術】従来、複数の測定用マイク及びマイクア
ンプを用いて音響特性の測定を行う音響特性測定装置に
おいては、パワーアンプ及びスピーカ等から成る音声再
生装置に供給される入力信号及び各マイクアンプからの
出力信号に高速フーリエ変換（ＦＦＴ）処理を施して入
力信号及び出力信号のパワースペクトル及びクロススペ
クトルを求め、このクロススペクトルの平均値を入力信
号のパワースペクトルの平均値で除算することにより、
マイクの位置の音響特性を測定している。
ンプを用いて音響特性の測定を行う音響特性測定装置に
おいては、パワーアンプ及びスピーカ等から成る音声再
生装置に供給される入力信号及び各マイクアンプからの
出力信号に高速フーリエ変換（ＦＦＴ）処理を施して入
力信号及び出力信号のパワースペクトル及びクロススペ
クトルを求め、このクロススペクトルの平均値を入力信
号のパワースペクトルの平均値で除算することにより、
マイクの位置の音響特性を測定している。
【０００３】具体的には、図５に示す音響特性の測定手
順のフローチャートを用いて説明する。
順のフローチャートを用いて説明する。
【０００４】ここで、図６は音響特性測定装置に用いる
信号を示すものであり、例えば音声データである入力信
号ｘ（ｔ）が被測定物１に供給されると共に、入力ＩＮ
1 の信号として入力端子２に出力される。また、上記被
測定物１から出力される、上記入力信号ｘ（ｔ）に対応
する出力信号ｙ（ｔ）は、入力ＩＮ2 の信号として入力
端子３に出力される。音響特性測定装置では、上記入力
端子２、３に出力される入力信号ｘ（ｔ）及び出力信号
ｙ（ｔ）を用いて、被測定物の伝達関数を測定して音響
特性を求める。
信号を示すものであり、例えば音声データである入力信
号ｘ（ｔ）が被測定物１に供給されると共に、入力ＩＮ
1 の信号として入力端子２に出力される。また、上記被
測定物１から出力される、上記入力信号ｘ（ｔ）に対応
する出力信号ｙ（ｔ）は、入力ＩＮ2 の信号として入力
端子３に出力される。音響特性測定装置では、上記入力
端子２、３に出力される入力信号ｘ（ｔ）及び出力信号
ｙ（ｔ）を用いて、被測定物の伝達関数を測定して音響
特性を求める。
【０００５】先ず、図５のステップＳ１で、入力Ｉ
Ｎ1 、ＩＮ2 の音声データを取り込み、ステップＳ２
で、入力ＩＮ1 の音声データに対して入力ＩＮ2 の音声
データが高速フーリエ変換のポイント数の範囲を越えて
遅延している場合を考慮して、入力ＩＮ2 の音声データ
としてはそのまま先頭から取り出し、入力ＩＮ1 の音声
データとしては補正分の音声データの先頭から遅延させ
て取り出す。
Ｎ1 、ＩＮ2 の音声データを取り込み、ステップＳ２
で、入力ＩＮ1 の音声データに対して入力ＩＮ2 の音声
データが高速フーリエ変換のポイント数の範囲を越えて
遅延している場合を考慮して、入力ＩＮ2 の音声データ
としてはそのまま先頭から取り出し、入力ＩＮ1 の音声
データとしては補正分の音声データの先頭から遅延させ
て取り出す。
【０００６】次に、ステップＳ３で入力ＩＮ1 の元の波
形に時間窓を乗算した後にＦＦＴ処理を行い、また、ス
テップＳ４で入力ＩＮ2 の元の波形に時間窓を乗算した
後にＦＦＴ処理を行うことにより、入力ＩＮ1 、ＩＮ2
の周波数スペクトルを求める。
形に時間窓を乗算した後にＦＦＴ処理を行い、また、ス
テップＳ４で入力ＩＮ2 の元の波形に時間窓を乗算した
後にＦＦＴ処理を行うことにより、入力ＩＮ1 、ＩＮ2
の周波数スペクトルを求める。
【０００７】これにより、ステップＳ５で入力ＩＮ1 の
周波数ポイント毎のパワースペクトルを算出し、ステッ
プＳ６で入力ＩＮ2 の周波数ポイント毎のパワースペク
トルを算出する。さらに、ステップＳ７で上記入力ＩＮ
1 、ＩＮ2 のパワースペクトルを用いてクロススペクト
ルを算出する。
周波数ポイント毎のパワースペクトルを算出し、ステッ
プＳ６で入力ＩＮ2 の周波数ポイント毎のパワースペク
トルを算出する。さらに、ステップＳ７で上記入力ＩＮ
1 、ＩＮ2 のパワースペクトルを用いてクロススペクト
ルを算出する。
【０００８】この後、ステップＳ８で、上記ステップＳ
３〜Ｓ７までの処理動作を、各周波数ポイント毎に平均
する回数Ｎ回分行ったか否かを判別し、Ｎ回行っていな
いならばステップＳ３に戻り、ステップＳ３〜Ｓ７まで
の処理動作を行う。また、Ｎ回行ったならばステップＳ
９に進んで、各周波数ポイント毎に入力ＩＮ1 、ＩＮ2
のパワースペクトルの平均値、及びクロススペクトルの
平均値を求める。
３〜Ｓ７までの処理動作を、各周波数ポイント毎に平均
する回数Ｎ回分行ったか否かを判別し、Ｎ回行っていな
いならばステップＳ３に戻り、ステップＳ３〜Ｓ７まで
の処理動作を行う。また、Ｎ回行ったならばステップＳ
９に進んで、各周波数ポイント毎に入力ＩＮ1 、ＩＮ2
のパワースペクトルの平均値、及びクロススペクトルの
平均値を求める。
【０００９】さらに、ステップＳ１０で、入力ＩＮ1 、
ＩＮ2 のパワースペクトルの平均値及びクロススペクト
ルの平均値を用いて、互いに干渉する光波の性質である
干渉性いわゆるコヒーレンスを計算し、ステップＳ１１
で伝達関数を計算する。
ＩＮ2 のパワースペクトルの平均値及びクロススペクト
ルの平均値を用いて、互いに干渉する光波の性質である
干渉性いわゆるコヒーレンスを計算し、ステップＳ１１
で伝達関数を計算する。
【００１０】この後、ステップＳ１２で、コヒーレンス
の値が一定値以上であるか否かを判別する信頼性判別処
理を各周波数ポイント毎に行い、コヒーレンスの値が一
定値以上である伝達関数のみを測定結果とする。
の値が一定値以上であるか否かを判別する信頼性判別処
理を各周波数ポイント毎に行い、コヒーレンスの値が一
定値以上である伝達関数のみを測定結果とする。
【００１１】さらに、ステップＳ１３で、各周波数ポイ
ント毎の伝達関数の平均化いわゆるアベレージングが設
定されているか否かを判別する。アベレージングの設定
がされていないならば音響特性の測定処理を終了する。
しかし、アベレージングの設定がされているならば、そ
の周波数ポイントの伝達関数を音響特性として表示す
る。
ント毎の伝達関数の平均化いわゆるアベレージングが設
定されているか否かを判別する。アベレージングの設定
がされていないならば音響特性の測定処理を終了する。
しかし、アベレージングの設定がされているならば、そ
の周波数ポイントの伝達関数を音響特性として表示す
る。
【００１２】また、ステップＳ１４に進んで、全ての測
定が終了したか否かを判別する。これにより、全ての測
定動作が終了していると判別されるならば、音響特性の
測定処理を終了する。また、全ての測定が終了していな
いと判別されるならば、ステップＳ１に戻り、入力ＩＮ
1 、ＩＮ2 の信号を取り込んで、ステップＳ２〜Ｓ１３
までの音響特性の測定処理を繰り返して行う。
定が終了したか否かを判別する。これにより、全ての測
定動作が終了していると判別されるならば、音響特性の
測定処理を終了する。また、全ての測定が終了していな
いと判別されるならば、ステップＳ１に戻り、入力ＩＮ
1 、ＩＮ2 の信号を取り込んで、ステップＳ２〜Ｓ１３
までの音響特性の測定処理を繰り返して行う。
【００１３】
【発明が解決しようとする課題】ところで、従来の音響
特性測定装置においては、複数の測定用マイクを用い、
これらの測定用マイクの位置でそれぞれ音響特性の測定
を行い、各々の測定用マイクの伝達関数を別々の音響特
性として測定を行っている。即ち、上記音響特性測定装
置では、各々の測定用マイクの音響特性を測定するのみ
であり、結果的に複数の異なる音響特性を得ることにな
る。
特性測定装置においては、複数の測定用マイクを用い、
これらの測定用マイクの位置でそれぞれ音響特性の測定
を行い、各々の測定用マイクの伝達関数を別々の音響特
性として測定を行っている。即ち、上記音響特性測定装
置では、各々の測定用マイクの音響特性を測定するのみ
であり、結果的に複数の異なる音響特性を得ることにな
る。
【００１４】しかし、複数の測定用マイクに対して音声
再生装置が１系統しかない場合に、複数の音響特性、即
ち測定結果を元にして１系統のイコライザを動作させ、
音響特性を調整することは容易ではない。
再生装置が１系統しかない場合に、複数の音響特性、即
ち測定結果を元にして１系統のイコライザを動作させ、
音響特性を調整することは容易ではない。
【００１５】そこで、本発明は上述の実情に鑑み、複数
の測定用マイクのそれぞれの音響特性を用いて音声の再
生を容易に調整することができる音響特性測定装置を提
供するものである。
の測定用マイクのそれぞれの音響特性を用いて音声の再
生を容易に調整することができる音響特性測定装置を提
供するものである。
【００１６】
【課題を解決するための手段】本発明に係る音響特性測
定装置は、入力信号及び出力信号に直交変換処理を施す
直交変換手段と、上記直交変換手段からの入力信号及び
該入力信号に対応する出力信号のスペクトルを用いて入
力信号及び出力信号の各パワースペクトルを算出し、ま
た、上記入力信号及び出力信号の各スペクトルからクロ
ススペクトルを算出し、上記入力信号及び出力信号の各
パワースペクトル及びクロススペクトルを用いて入力信
号及び出力信号のパワースペクトルの平均値及びクロス
スペクトルの平均値を求めるスペクトル算出手段と、上
記スペクトル算出手段からの入力信号及び出力信号のパ
ワースペクトルの平均値、及びクロススペクトルの平均
値を用いて各周波数成分の伝達特性を算出する伝達特性
算出手段と、上記伝達特性算出手段からの各周波数成分
の伝達特性が一定値以上であるか否かを判別し、一定値
以上の周波数成分の伝達特性のみを平均化する平均化手
段と、上記平均化手段からの複数の測定用マイク毎の伝
達特性の平均値を平均化して音響特性を得る音響特性算
出手段とを有して成ることにより上述した課題を解決す
る。
定装置は、入力信号及び出力信号に直交変換処理を施す
直交変換手段と、上記直交変換手段からの入力信号及び
該入力信号に対応する出力信号のスペクトルを用いて入
力信号及び出力信号の各パワースペクトルを算出し、ま
た、上記入力信号及び出力信号の各スペクトルからクロ
ススペクトルを算出し、上記入力信号及び出力信号の各
パワースペクトル及びクロススペクトルを用いて入力信
号及び出力信号のパワースペクトルの平均値及びクロス
スペクトルの平均値を求めるスペクトル算出手段と、上
記スペクトル算出手段からの入力信号及び出力信号のパ
ワースペクトルの平均値、及びクロススペクトルの平均
値を用いて各周波数成分の伝達特性を算出する伝達特性
算出手段と、上記伝達特性算出手段からの各周波数成分
の伝達特性が一定値以上であるか否かを判別し、一定値
以上の周波数成分の伝達特性のみを平均化する平均化手
段と、上記平均化手段からの複数の測定用マイク毎の伝
達特性の平均値を平均化して音響特性を得る音響特性算
出手段とを有して成ることにより上述した課題を解決す
る。
【００１７】ここで、上記入力信号及び出力信号のパワ
ースペクトルの平均値、及びクロススペクトルの平均値
を用いてコヒーレンスの値を計算するコヒーレンス算出
手段を有し、上記平均化手段では上記入力信号のパワー
スペクトルの平均値のレベルと上記コヒーレンス算出手
段からのコヒーレンスの値とが一定値以上である周波数
成分の伝達特性のみを平均化することを特徴とする。
ースペクトルの平均値、及びクロススペクトルの平均値
を用いてコヒーレンスの値を計算するコヒーレンス算出
手段を有し、上記平均化手段では上記入力信号のパワー
スペクトルの平均値のレベルと上記コヒーレンス算出手
段からのコヒーレンスの値とが一定値以上である周波数
成分の伝達特性のみを平均化することを特徴とする。
【００１８】また、本発明に係る音響特性測定装置は、
入力信号及び出力信号に直交変換処理を施す直交変換手
段と、上記直交変換手段からの入力信号及び該入力信号
に対応する出力信号のスペクトルを用いて入力信号及び
出力信号の各パワースペクトルを算出し、また、上記入
力信号及び出力信号の各スペクトルからクロススペクト
ルを算出し、上記入力信号及び出力信号の各パワースペ
クトル及びクロススペクトルを用いて入力信号及び出力
信号のパワースペクトルの平均値及びクロススペクトル
の平均値を求めるスペクトル算出手段と、上記スペクト
ル算出手段からの入力信号及び出力信号のパワースペク
トルの平均値、及びクロススペクトルの平均値を用いて
各周波数成分の伝達特性を算出する伝達特性算出手段
と、上記伝達特性算出手段からの各周波数成分の伝達特
性が一定値以上であるか否かを判別し、一定値以上の周
波数成分の伝達特性のみを平均化する平均化手段と、上
記平均化手段からの複数の測定用マイク毎の伝達特性の
平均値から少なくとも１つの測定用マイクの伝達特性の
平均値を独立に取り出す音響特性算出手段とを有して成
ることを特徴とする。
入力信号及び出力信号に直交変換処理を施す直交変換手
段と、上記直交変換手段からの入力信号及び該入力信号
に対応する出力信号のスペクトルを用いて入力信号及び
出力信号の各パワースペクトルを算出し、また、上記入
力信号及び出力信号の各スペクトルからクロススペクト
ルを算出し、上記入力信号及び出力信号の各パワースペ
クトル及びクロススペクトルを用いて入力信号及び出力
信号のパワースペクトルの平均値及びクロススペクトル
の平均値を求めるスペクトル算出手段と、上記スペクト
ル算出手段からの入力信号及び出力信号のパワースペク
トルの平均値、及びクロススペクトルの平均値を用いて
各周波数成分の伝達特性を算出する伝達特性算出手段
と、上記伝達特性算出手段からの各周波数成分の伝達特
性が一定値以上であるか否かを判別し、一定値以上の周
波数成分の伝達特性のみを平均化する平均化手段と、上
記平均化手段からの複数の測定用マイク毎の伝達特性の
平均値から少なくとも１つの測定用マイクの伝達特性の
平均値を独立に取り出す音響特性算出手段とを有して成
ることを特徴とする。
【００１９】
【作用】本発明においては、入力信号及びこの入力信号
が音響信号に変換され、複数の測定用マイクにより収音
されて得られた出力信号を用いて、入力信号及び出力信
号のスペクトルから周波数成分毎の伝達特性を算出し、
一定値以上の周波数成分の伝達特性のみを平均化して１
つの測定用マイクの各周波数成分の伝達特性の平均値を
得、さらに、複数の測定用マイクの各周波数成分の伝達
特性の平均値を平均化して音響特性を得ることにより、
測定用マイクの位置による音響特性のばらつきを抑える
ことができる。
が音響信号に変換され、複数の測定用マイクにより収音
されて得られた出力信号を用いて、入力信号及び出力信
号のスペクトルから周波数成分毎の伝達特性を算出し、
一定値以上の周波数成分の伝達特性のみを平均化して１
つの測定用マイクの各周波数成分の伝達特性の平均値を
得、さらに、複数の測定用マイクの各周波数成分の伝達
特性の平均値を平均化して音響特性を得ることにより、
測定用マイクの位置による音響特性のばらつきを抑える
ことができる。
【００２０】また、入力信号及び出力信号のパワースペ
クトルの平均値、及びクロススペクトルの平均値を用い
てコヒーレンスの値を計算し、上記入力信号のパワース
ペクトルの平均値のレベルとコヒーレンスの値とが一定
値以上である周波数成分の伝達特性のみを平均化するこ
とにより、求める伝達関数の信頼性を高めることができ
る。
クトルの平均値、及びクロススペクトルの平均値を用い
てコヒーレンスの値を計算し、上記入力信号のパワース
ペクトルの平均値のレベルとコヒーレンスの値とが一定
値以上である周波数成分の伝達特性のみを平均化するこ
とにより、求める伝達関数の信頼性を高めることができ
る。
【００２１】また、複数の測定用マイクの伝達特性の平
均値から少なくとも１つの測定用マイクの伝達特性の平
均値を独立に取り出すことにより、任意の測定用マイク
の測定結果を音響特性として得ることができる。
均値から少なくとも１つの測定用マイクの伝達特性の平
均値を独立に取り出すことにより、任意の測定用マイク
の測定結果を音響特性として得ることができる。
【００２２】
【実施例】以下、本発明の好ましい実施例について、図
面を参照しながら説明する。図１には、本発明に係る音
響特性測定装置を用いたシステムの概略的な構成を示
す。
面を参照しながら説明する。図１には、本発明に係る音
響特性測定装置を用いたシステムの概略的な構成を示
す。
【００２３】この図１のシステムの入力端子４１からの
音声データは、イコライザ４１に入力されると共に、音
響特性測定器４５に入力される。
音声データは、イコライザ４１に入力されると共に、音
響特性測定器４５に入力される。
【００２４】上記イコライザ４１で調整された音声デー
タは、パワーアンプ４３を介してスピーカ４４から出力
される。
タは、パワーアンプ４３を介してスピーカ４４から出力
される。
【００２５】また、このシステムは５本の測定用マイク
を装備している。これらの測定用マイクＭＣ1 〜ＭＣ5
はそれぞれマイクアンプ４７1 〜４７5 に接続されてお
り、さらに、マイクアンプ４７1 〜４７5 は信号切換器
４６内の端子ｂ〜ｆにそれぞれ接続されている。この信
号切換器４６内の端子ｂ〜ｆのいずれかの端子が切換接
続されることにより、上記測定用マイクＭＣ1 〜ＭＣ5
から入力されるいずれかの音声データが、対応するマイ
クアンプで増幅され、信号切換器４６を介して音響特性
測定器４５に入力される。
を装備している。これらの測定用マイクＭＣ1 〜ＭＣ5
はそれぞれマイクアンプ４７1 〜４７5 に接続されてお
り、さらに、マイクアンプ４７1 〜４７5 は信号切換器
４６内の端子ｂ〜ｆにそれぞれ接続されている。この信
号切換器４６内の端子ｂ〜ｆのいずれかの端子が切換接
続されることにより、上記測定用マイクＭＣ1 〜ＭＣ5
から入力されるいずれかの音声データが、対応するマイ
クアンプで増幅され、信号切換器４６を介して音響特性
測定器４５に入力される。
【００２６】これにより、上記入力端子４１からの音声
データ及び測定用マイクＭＣ1 〜ＭＣ5 の内のいずれか
の音声データの２チャンネルのアナログ音声信号を用い
て伝達関数を計算し、振幅及び位相等の伝達特性を得
て、それらの伝達特性を各周波数成分毎に加算して平均
化いわゆるアベレージングを行う。さらに、測定用マイ
クＭＣ1 〜ＭＣ5 毎の伝達特性のアベレージングを行う
ことによって１つの音響特性を得て、パワーアンプ４３
及びスピーカ４４から成る１系統の音響再生装置から出
力される信号を調整する。
データ及び測定用マイクＭＣ1 〜ＭＣ5 の内のいずれか
の音声データの２チャンネルのアナログ音声信号を用い
て伝達関数を計算し、振幅及び位相等の伝達特性を得
て、それらの伝達特性を各周波数成分毎に加算して平均
化いわゆるアベレージングを行う。さらに、測定用マイ
クＭＣ1 〜ＭＣ5 毎の伝達特性のアベレージングを行う
ことによって１つの音響特性を得て、パワーアンプ４３
及びスピーカ４４から成る１系統の音響再生装置から出
力される信号を調整する。
【００２７】また、信号切換器４６が端子ａに切り換え
られたときには、入力端子４１からの音声データとイコ
ライザ４２により調整された音声データとの間の音響特
性も測定することができる。
られたときには、入力端子４１からの音声データとイコ
ライザ４２により調整された音声データとの間の音響特
性も測定することができる。
【００２８】次に、音響特性測定器の概略的な構成を図
２に示し、音響特性の測定について以下に説明する。
２に示し、音響特性の測定について以下に説明する。
【００２９】ここで、伝達関数の測定に必要な音声信号
のサンプル数は、後述する高速フーリエ変換（ＦＦＴ）
のポイント数、及び後述するパワースペクトル及びクロ
ススペクトルの平均する回数を用いて算出される。例え
ば、上記伝達関数の測定に必要なサンプル数をＳ、高速
フーリエ変換のポイント数をＰ、平均する回数をＮとす
る場合には、上記伝達関数の測定に必要なサンプル数Ｓ
は以下の（１）式で表される。
のサンプル数は、後述する高速フーリエ変換（ＦＦＴ）
のポイント数、及び後述するパワースペクトル及びクロ
ススペクトルの平均する回数を用いて算出される。例え
ば、上記伝達関数の測定に必要なサンプル数をＳ、高速
フーリエ変換のポイント数をＰ、平均する回数をＮとす
る場合には、上記伝達関数の測定に必要なサンプル数Ｓ
は以下の（１）式で表される。
【００３０】Ｓ＝Ｐ×Ｎ ・・・（１） 具体的には、高速フーリエ変換のポイント数Ｐが２０４
８、平均する回数Ｎが１０のときには、サンプル数Ｓ
は、 ２０４８×１０＝２０４８０ の値となる。
８、平均する回数Ｎが１０のときには、サンプル数Ｓ
は、 ２０４８×１０＝２０４８０ の値となる。
【００３１】例えば、図２の入力端子１１に入力される
入力ＩＮ1 の信号は、図１のマイクＭＣ1 〜ＭＣ5 の内
で信号切換器４６によって切換選択された１つのマイク
に接続されるマイクアンプ４７からの入力信号がディジ
タル変換された信号であり、図２の入力端子２１に入力
される入力ＩＮ2 の信号は、図１のパワーアンプ４３に
出力される信号がディジタル変換された信号であるとす
る。
入力ＩＮ1 の信号は、図１のマイクＭＣ1 〜ＭＣ5 の内
で信号切換器４６によって切換選択された１つのマイク
に接続されるマイクアンプ４７からの入力信号がディジ
タル変換された信号であり、図２の入力端子２１に入力
される入力ＩＮ2 の信号は、図１のパワーアンプ４３に
出力される信号がディジタル変換された信号であるとす
る。
【００３２】ここで、入力ＩＮ1 の信号に対して入力Ｉ
Ｎ2 の信号が高速フーリエ変換のポイント数の範囲を越
えて遅延している場合を考慮して、入力ＩＮ2 の信号と
してはそのまま信号データの先頭から取り出し、入力Ｉ
Ｎ1 の信号としては補正分の信号データの先頭から遅延
させて取り出す。
Ｎ2 の信号が高速フーリエ変換のポイント数の範囲を越
えて遅延している場合を考慮して、入力ＩＮ2 の信号と
してはそのまま信号データの先頭から取り出し、入力Ｉ
Ｎ1 の信号としては補正分の信号データの先頭から遅延
させて取り出す。
【００３３】上記入力ＩＮ1 の信号は入力端子１１から
乗算器１３に出力されると共に、この入力ＩＮ1 の信号
に対応する入力ＩＮ2 の信号が入力端子２１から乗算器
２３に出力される。
乗算器１３に出力されると共に、この入力ＩＮ1 の信号
に対応する入力ＩＮ2 の信号が入力端子２１から乗算器
２３に出力される。
【００３４】ここで、所定のサンプル数の音声波形の切
り出しの際には、切り出し区間の両端に急激な変化が起
こらないようにすると共に、スペクトル領域では信号の
スペクトルに窓関数のフーリエ変換の畳み込み、即ち重
みつき移動平均を行うために、元の波形に時間窓を乗算
する処理が必要である。よって、窓関数発生器１２、２
２からは、上記乗算器１３、２３に送られる信号の帯域
に対応する窓関数が発生されて乗算器１３、２３に供給
される。これにより、上記乗算器１３、２３では帯域の
信号と窓関数とが乗算され、この乗算された信号は、そ
れぞれＦＦＴ解析器１４、２４に送られる。
り出しの際には、切り出し区間の両端に急激な変化が起
こらないようにすると共に、スペクトル領域では信号の
スペクトルに窓関数のフーリエ変換の畳み込み、即ち重
みつき移動平均を行うために、元の波形に時間窓を乗算
する処理が必要である。よって、窓関数発生器１２、２
２からは、上記乗算器１３、２３に送られる信号の帯域
に対応する窓関数が発生されて乗算器１３、２３に供給
される。これにより、上記乗算器１３、２３では帯域の
信号と窓関数とが乗算され、この乗算された信号は、そ
れぞれＦＦＴ解析器１４、２４に送られる。
【００３５】上記ＦＦＴ解析器１４、２４では、送られ
た信号データに高速フーリエ変換処理を施すことによ
り、入力ＩＮ1 の信号の周波数スペクトル及び入力ＩＮ
2 の信号の周波数スペクトルが求められる。
た信号データに高速フーリエ変換処理を施すことによ
り、入力ＩＮ1 の信号の周波数スペクトル及び入力ＩＮ
2 の信号の周波数スペクトルが求められる。
【００３６】上記入力ＩＮ1 のスペクトルの複素データ
をＸ（ｋ）とするとき、この複素データＸ（ｋ）は乗算
器１６に送られると共に、複素共役変換器１５にも送ら
れる。この複素共役変換器１５では送られた複素データ
Ｘ（ｋ）が複素共役データＸ* ( ｋ）に変換されて乗算
器１６に送られる。この乗算器１６では上記ＦＦＴ解析
器１４からの複素データＸ（ｋ）と上記複素共役データ
Ｘ* ( ｋ）とが乗算されて入力ＩＮ1 のパワースペクト
ルＸ* ( ｋ）Ｘ（ｋ）が求められる。このパワースペク
トルＸ* ( ｋ）Ｘ（ｋ）はレジスタ３１ａに記憶され
る。
をＸ（ｋ）とするとき、この複素データＸ（ｋ）は乗算
器１６に送られると共に、複素共役変換器１５にも送ら
れる。この複素共役変換器１５では送られた複素データ
Ｘ（ｋ）が複素共役データＸ* ( ｋ）に変換されて乗算
器１６に送られる。この乗算器１６では上記ＦＦＴ解析
器１４からの複素データＸ（ｋ）と上記複素共役データ
Ｘ* ( ｋ）とが乗算されて入力ＩＮ1 のパワースペクト
ルＸ* ( ｋ）Ｘ（ｋ）が求められる。このパワースペク
トルＸ* ( ｋ）Ｘ（ｋ）はレジスタ３１ａに記憶され
る。
【００３７】同様にして、上記入力ＩＮ2 のスペクトル
の複素データをＹ（ｋ）とすると、この複素データＹ
（ｋ）は乗算器２６に送られると共に、複素共役変換器
２５に送られる。この複素共役変換器２５では送られた
複素データＹ（ｋ）が複素共役データＹ* ( ｋ）に変換
されて、この複素共役データＹ* ( ｋ）は乗算器２６に
送られる。この乗算器２６では上記ＦＦＴ解析器２４か
らの複素データＹ（ｋ）と上記複素共役データＹ* (
ｋ）とが乗算されて入力ＩＮ2 のパワースペクトルＹ*
( ｋ）Ｙ（ｋ）が求められる。このパワースペクトルＹ
* ( ｋ）Ｙ（ｋ）はレジスタ３１ｃに記憶される。
の複素データをＹ（ｋ）とすると、この複素データＹ
（ｋ）は乗算器２６に送られると共に、複素共役変換器
２５に送られる。この複素共役変換器２５では送られた
複素データＹ（ｋ）が複素共役データＹ* ( ｋ）に変換
されて、この複素共役データＹ* ( ｋ）は乗算器２６に
送られる。この乗算器２６では上記ＦＦＴ解析器２４か
らの複素データＹ（ｋ）と上記複素共役データＹ* (
ｋ）とが乗算されて入力ＩＮ2 のパワースペクトルＹ*
( ｋ）Ｙ（ｋ）が求められる。このパワースペクトルＹ
* ( ｋ）Ｙ（ｋ）はレジスタ３１ｃに記憶される。
【００３８】また、上記複素共役変換器１５から出力さ
れる入力ＩＮ1 のスペクトルの複素共役データＸ* (
ｋ）と上記ＦＦＴ解析器２４から出力される入力ＩＮ2
のスペクトルの複素データＹ（ｋ）とは乗算器１７で乗
算されて、クロススペクトルＸ* ( ｋ）Ｙ（ｋ）が求め
られ、このクロススペクトルＸ* ( ｋ）Ｙ（ｋ）はレジ
スタ３１ｂに記憶される。
れる入力ＩＮ1 のスペクトルの複素共役データＸ* (
ｋ）と上記ＦＦＴ解析器２４から出力される入力ＩＮ2
のスペクトルの複素データＹ（ｋ）とは乗算器１７で乗
算されて、クロススペクトルＸ* ( ｋ）Ｙ（ｋ）が求め
られ、このクロススペクトルＸ* ( ｋ）Ｙ（ｋ）はレジ
スタ３１ｂに記憶される。
【００３９】ここで、上記レジスタ３１ａ、３１ｂ、３
１ｃは、上記平均する回数Ｎに対応してＮ個から成るも
のであり、Ｎ個の入力ＩＮ1 のパワースペクトルＸ* (
ｋ）Ｘ（ｋ）、入力ＩＮ2 のパワースペクトルＹ* (
ｋ）Ｙ（ｋ）、及びクロススペクトルＸ* ( ｋ）Ｙ
（ｋ）をそれぞれ記憶するものである。
１ｃは、上記平均する回数Ｎに対応してＮ個から成るも
のであり、Ｎ個の入力ＩＮ1 のパワースペクトルＸ* (
ｋ）Ｘ（ｋ）、入力ＩＮ2 のパワースペクトルＹ* (
ｋ）Ｙ（ｋ）、及びクロススペクトルＸ* ( ｋ）Ｙ
（ｋ）をそれぞれ記憶するものである。
【００４０】この後、上記レジスタ３１ａ、３１ｂ、３
１ｃにそれぞれ記憶されたＮ個の入力ＩＮ1 のパワース
ペクトルＸ* ( ｋ）Ｘ（ｋ）、クロススペクトルＸ* (
ｋ）Ｙ（ｋ）、及び入力ＩＮ2 のパワースペクトルＹ*
( ｋ）Ｙ（ｋ）は、それぞれ対応する平均値化回路３２
ａ、３２ｂ、３２ｃに送られて、入力ＩＮ1 のパワース
ペクトルの平均値Ｐ1 （ｋ）、クロススペクトルの平均
値Ｃ（ｋ）、及び入力ＩＮ2 のパワースペクトルの平均
値Ｐ2 （ｋ）が計算される。
１ｃにそれぞれ記憶されたＮ個の入力ＩＮ1 のパワース
ペクトルＸ* ( ｋ）Ｘ（ｋ）、クロススペクトルＸ* (
ｋ）Ｙ（ｋ）、及び入力ＩＮ2 のパワースペクトルＹ*
( ｋ）Ｙ（ｋ）は、それぞれ対応する平均値化回路３２
ａ、３２ｂ、３２ｃに送られて、入力ＩＮ1 のパワース
ペクトルの平均値Ｐ1 （ｋ）、クロススペクトルの平均
値Ｃ（ｋ）、及び入力ＩＮ2 のパワースペクトルの平均
値Ｐ2 （ｋ）が計算される。
【００４１】ここで、上記入力ＩＮ1 のパワースペクト
ルの平均値Ｐ1 （ｋ）は（２）式で表され、入力ＩＮ2
のパワースペクトルの平均値Ｐ2 （ｋ）は（３）式で表
され、クロススペクトルの平均値Ｃ（ｋ）は（４）式で
表される。
ルの平均値Ｐ1 （ｋ）は（２）式で表され、入力ＩＮ2
のパワースペクトルの平均値Ｐ2 （ｋ）は（３）式で表
され、クロススペクトルの平均値Ｃ（ｋ）は（４）式で
表される。
【００４２】
【数１】
【００４３】
【数２】
【００４４】
【数３】
【００４５】上記入力ＩＮ1 のパワースペクトルの平均
値Ｐ1 （ｋ）及びクロススペクトルの平均値Ｃ（ｋ）は
伝達関数演算器３５に送られ、この伝達関数演算器３５
において被測定物の伝達関数Ｈ（ｋ）が計算される。具
体的には、振幅及び位相が伝達関数として計算される。
この伝達関数の値は出力端子３９から出力される。
値Ｐ1 （ｋ）及びクロススペクトルの平均値Ｃ（ｋ）は
伝達関数演算器３５に送られ、この伝達関数演算器３５
において被測定物の伝達関数Ｈ（ｋ）が計算される。具
体的には、振幅及び位相が伝達関数として計算される。
この伝達関数の値は出力端子３９から出力される。
【００４６】振幅及び位相を含めた伝達関数Ｈ（ｋ）は
以下の（５）式で表され、振幅のみの伝達関数Ｈ（ｋ）
は以下の（６）式で表される。
以下の（５）式で表され、振幅のみの伝達関数Ｈ（ｋ）
は以下の（６）式で表される。
【００４７】
【数４】
【００４８】
【数５】
【００４９】また、上記入力ＩＮ1 のパワースペクトル
の平均値Ｐ1 （ｋ）、入力ＩＮ2 のパワースペクトルの
平均値Ｐ2 （ｋ）、及びクロススペクトルの平均値Ｃ
（ｋ）は、コヒーレンス演算器３４に送られる。
の平均値Ｐ1 （ｋ）、入力ＩＮ2 のパワースペクトルの
平均値Ｐ2 （ｋ）、及びクロススペクトルの平均値Ｃ
（ｋ）は、コヒーレンス演算器３４に送られる。
【００５０】尚、上記平均値化回路３２ｂで求められた
クロススペクトルの平均値Ｃ（ｋ）は、複素共役変換器
３３に送られてクロススペクトルの平均値Ｃ（ｋ）の複
素共役データＣ* ( ｋ）が求められており、このクロス
スペクトルの平均値Ｃ（ｋ）の複素共役データＣ* (
ｋ）もコヒーレンス演算器３４に送られる。
クロススペクトルの平均値Ｃ（ｋ）は、複素共役変換器
３３に送られてクロススペクトルの平均値Ｃ（ｋ）の複
素共役データＣ* ( ｋ）が求められており、このクロス
スペクトルの平均値Ｃ（ｋ）の複素共役データＣ* (
ｋ）もコヒーレンス演算器３４に送られる。
【００５１】上記コヒーレンス演算器３４では、上記入
力ＩＮ1 のパワースペクトルの平均値Ｐ1 （ｋ）、入力
ＩＮ2 のパワースペクトルの平均値Ｐ2 （ｋ）、クロス
スペクトルの平均値Ｃ（ｋ）、及びクロススペクトルの
平均値Ｃ（ｋ）の複素共役データＣ* ( ｋ）を用いて、
互いに干渉する光波の性質である干渉性いわゆるコヒー
レンスが求められる。このコヒーレンスをｒとすると、
コヒーレンスｒは以下の（７）式で表される。
力ＩＮ1 のパワースペクトルの平均値Ｐ1 （ｋ）、入力
ＩＮ2 のパワースペクトルの平均値Ｐ2 （ｋ）、クロス
スペクトルの平均値Ｃ（ｋ）、及びクロススペクトルの
平均値Ｃ（ｋ）の複素共役データＣ* ( ｋ）を用いて、
互いに干渉する光波の性質である干渉性いわゆるコヒー
レンスが求められる。このコヒーレンスをｒとすると、
コヒーレンスｒは以下の（７）式で表される。
【００５２】
【数６】
【００５３】ここで、各データの周波数ポイントにおい
て、入力ＩＮ1 の信号のパワースペクトルの平均値のレ
ベルとコヒーレンスの値とが一定値以上である場合の周
波数成分の伝達特性のみを測定結果としてデータを更新
する。
て、入力ＩＮ1 の信号のパワースペクトルの平均値のレ
ベルとコヒーレンスの値とが一定値以上である場合の周
波数成分の伝達特性のみを測定結果としてデータを更新
する。
【００５４】また、さらに次の音響特性の測定が継続さ
れるならば、再び入力ＩＮ1 及び入力ＩＮ2 の信号を取
り込んで、上述の処理を繰り返す。
れるならば、再び入力ＩＮ1 及び入力ＩＮ2 の信号を取
り込んで、上述の処理を繰り返す。
【００５５】図１では、上記音響特性測定器４５からの
コントロール信号によって信号切換器４６の切換を制御
することにより、先ず、マイクアンプ４７1 からの出力
を音響特性測定器４５に入力し、測定用マイクＭＣ1 の
伝達特性を測定する。この後、信号切換器４６を切換制
御することにより、マイクアンプ４７2 、４７3 、４７
4 、４７5 からの出力を順次選択して、測定用マイクＭ
Ｃ2 、ＭＣ3 、ＭＣ4、ＭＣ5 の伝達特性をそれぞれ測
定する。
コントロール信号によって信号切換器４６の切換を制御
することにより、先ず、マイクアンプ４７1 からの出力
を音響特性測定器４５に入力し、測定用マイクＭＣ1 の
伝達特性を測定する。この後、信号切換器４６を切換制
御することにより、マイクアンプ４７2 、４７3 、４７
4 、４７5 からの出力を順次選択して、測定用マイクＭ
Ｃ2 、ＭＣ3 、ＭＣ4、ＭＣ5 の伝達特性をそれぞれ測
定する。
【００５６】また、この音響特性測定器４５では、１つ
の測定用マイクにおける各周波数ポイント毎の測定結果
の値を加算してアベレージングを行い、さらに、全ての
測定用マイクの測定結果を加算してアベレージングを行
い、１つの測定結果を音響特性として表示する。
の測定用マイクにおける各周波数ポイント毎の測定結果
の値を加算してアベレージングを行い、さらに、全ての
測定用マイクの測定結果を加算してアベレージングを行
い、１つの測定結果を音響特性として表示する。
【００５７】次に、図３に示す音響特性測定装置の測定
手順のフローチャートを用いて、音響特性の測定手順を
以下に説明する。
手順のフローチャートを用いて、音響特性の測定手順を
以下に説明する。
【００５８】先ず、ステップＳ２１でマイクアンプから
の出力の切換選択を行い、ステップＳ２２で、入力ＩＮ
1 の信号に対して入力ＩＮ2 の信号が高速フーリエ変換
のポイント数の範囲を越えて遅延している場合を考慮し
て、入力ＩＮ1 の信号は補正分の信号データの先頭から
遅延させて取り出すように遅延時間を設定する。
の出力の切換選択を行い、ステップＳ２２で、入力ＩＮ
1 の信号に対して入力ＩＮ2 の信号が高速フーリエ変換
のポイント数の範囲を越えて遅延している場合を考慮し
て、入力ＩＮ1 の信号は補正分の信号データの先頭から
遅延させて取り出すように遅延時間を設定する。
【００５９】この後、ステップＳ２３で、入力ＩＮ1 及
び入力ＩＮ2 からのアナログ音声信号を取り込んでディ
ジタル変換し、このディジタル信号から音響特性の測定
に必要なサンプル数を取り出す。
び入力ＩＮ2 からのアナログ音声信号を取り込んでディ
ジタル変換し、このディジタル信号から音響特性の測定
に必要なサンプル数を取り出す。
【００６０】そして、ステップＳ２４で、入力ＩＮ1 の
信号に窓関数を乗算し、この出力にＦＦＴ処理を施して
入力ＩＮ1 の信号の周波数スペクトルを求め、また、入
力ＩＮ2 の信号に窓関数を乗算し、この出力にＦＦＴ処
理を施して入力ＩＮ2 の信号の周波数スペクトルを求め
る。さらに、入力ＩＮ1 の信号及び入力ＩＮ2 の信号の
周波数スペクトルを用いてそれぞれのパワースペクトル
を算出し、上記算出された入力ＩＮ1 のパワースペクト
ルと入力ＩＮ2 のパワースペクトルとを用いてクロスス
ペクトルを算出する。
信号に窓関数を乗算し、この出力にＦＦＴ処理を施して
入力ＩＮ1 の信号の周波数スペクトルを求め、また、入
力ＩＮ2 の信号に窓関数を乗算し、この出力にＦＦＴ処
理を施して入力ＩＮ2 の信号の周波数スペクトルを求め
る。さらに、入力ＩＮ1 の信号及び入力ＩＮ2 の信号の
周波数スペクトルを用いてそれぞれのパワースペクトル
を算出し、上記算出された入力ＩＮ1 のパワースペクト
ルと入力ＩＮ2 のパワースペクトルとを用いてクロスス
ペクトルを算出する。
【００６１】また、入力ＩＮ1 のパワースペクトル、入
力ＩＮ2 のパワースペクトル、及びクロススペクトルを
平均する回数Ｎ個分だけ求めるために、ステップＳ２５
において、ステップＳ２４の処理操作をＮ回行ったか否
かを判別する。これにより、入力ＩＮ1 及び入力ＩＮ2
のパワースペクトル、及びクロススペクトルのそれぞれ
の値を平均する回数分だけ求める。
力ＩＮ2 のパワースペクトル、及びクロススペクトルを
平均する回数Ｎ個分だけ求めるために、ステップＳ２５
において、ステップＳ２４の処理操作をＮ回行ったか否
かを判別する。これにより、入力ＩＮ1 及び入力ＩＮ2
のパワースペクトル、及びクロススペクトルのそれぞれ
の値を平均する回数分だけ求める。
【００６２】この後、ステップＳ２６において、上記求
められた平均する回数分の入力ＩＮ1 のパワースペクト
ル、入力ＩＮ2 のパワースペクトル、及びクロススペク
トルの値を用いて、入力ＩＮ1 のパワースペクトル、入
力ＩＮ2 のパワースペクトル、及びクロススペクトルの
平均値をそれぞれ算出する。
められた平均する回数分の入力ＩＮ1 のパワースペクト
ル、入力ＩＮ2 のパワースペクトル、及びクロススペク
トルの値を用いて、入力ＩＮ1 のパワースペクトル、入
力ＩＮ2 のパワースペクトル、及びクロススペクトルの
平均値をそれぞれ算出する。
【００６３】そして、ステップＳ２７で、入力ＩＮ1 の
パワースペクトル、入力ＩＮ2 のパワースペクトル、及
びクロススペクトルの平均値を用いてコヒーレンスの値
を求める。また、入力ＩＮ1 のパワースペクトル及びク
ロススペクトルの平均値を用いて伝達関数を算出する。
パワースペクトル、入力ＩＮ2 のパワースペクトル、及
びクロススペクトルの平均値を用いてコヒーレンスの値
を求める。また、入力ＩＮ1 のパワースペクトル及びク
ロススペクトルの平均値を用いて伝達関数を算出する。
【００６４】さらに、ステップＳ２８で、各周波数ポイ
ント毎に求めたコヒーレンス値及び伝達関数を用いて、
後述する信頼性判別処理により、各周波数ポイント毎に
伝達関数のアベレージングを行って平均値を求め、１つ
の測定用マイクの音響特性を外部の表示装置上に表示す
る。
ント毎に求めたコヒーレンス値及び伝達関数を用いて、
後述する信頼性判別処理により、各周波数ポイント毎に
伝達関数のアベレージングを行って平均値を求め、１つ
の測定用マイクの音響特性を外部の表示装置上に表示す
る。
【００６５】この後、ステップＳ２１に戻り、マイクア
ンプの出力を切り換えて上記ステップＳ２２〜Ｓ２８ま
での操作を行うことにより、図１の測定用マイクＭＣ1
〜ＭＣ5 の各伝達関数を求めて表示する。
ンプの出力を切り換えて上記ステップＳ２２〜Ｓ２８ま
での操作を行うことにより、図１の測定用マイクＭＣ1
〜ＭＣ5 の各伝達関数を求めて表示する。
【００６６】そして、ステップＳ２９に進んで、測定用
マイクＭＣ1 〜ＭＣ5 の各周波数ポイント毎の伝達関数
の平均値を用いてアベレージングを行い、各周波数ポイ
ント毎に１つの伝達関数の平均値を求め、音響特性とし
て表示する。
マイクＭＣ1 〜ＭＣ5 の各周波数ポイント毎の伝達関数
の平均値を用いてアベレージングを行い、各周波数ポイ
ント毎に１つの伝達関数の平均値を求め、音響特性とし
て表示する。
【００６７】この後、ステップＳ３０で、アベレージン
グ動作をＯＦＦするならば、音響特性の測定動作を終了
するが、さらに音響特性を測定するならば、ステップＳ
２１に戻って測定用マイクを切換選択してステップＳ２
２〜Ｓ２８までの動作を行い、ステップＳ２９で各測定
用マイクの伝達関数の平均値のアベレージングを行う。
グ動作をＯＦＦするならば、音響特性の測定動作を終了
するが、さらに音響特性を測定するならば、ステップＳ
２１に戻って測定用マイクを切換選択してステップＳ２
２〜Ｓ２８までの動作を行い、ステップＳ２９で各測定
用マイクの伝達関数の平均値のアベレージングを行う。
【００６８】次に、各周波数ポイント毎の信頼性判別処
理の手順のフローチャートを図４に示し、以下に説明す
る。
理の手順のフローチャートを図４に示し、以下に説明す
る。
【００６９】先ず、ステップＳ４１で、得られたコヒー
レンスの値が一定値以上であるか否かを判別し、一定値
以上であるならば 'ＯＫ’とし、一定値未満であるなら
ば 'ＮＧ’とする。この判別により、 'ＯＫ’とされた
ならば、ステップＳ４２に進んで、アベレージングを行
うことが設定されているか否かを判別し、また、 'Ｎ
Ｇ’とされたならば、ステップＳ４６に進んで、アベレ
ージングを行うことが設定されているか否かを判別す
る。
レンスの値が一定値以上であるか否かを判別し、一定値
以上であるならば 'ＯＫ’とし、一定値未満であるなら
ば 'ＮＧ’とする。この判別により、 'ＯＫ’とされた
ならば、ステップＳ４２に進んで、アベレージングを行
うことが設定されているか否かを判別し、また、 'Ｎ
Ｇ’とされたならば、ステップＳ４６に進んで、アベレ
ージングを行うことが設定されているか否かを判別す
る。
【００７０】上記ステップＳ４６において、アベレージ
ングを行うことが設定されていないと判別されるなら
ば、何も行わずにこの処理を終了する。
ングを行うことが設定されていないと判別されるなら
ば、何も行わずにこの処理を終了する。
【００７１】また、上記ステップＳ４２において、アベ
レージングを行うことが設定されていないならば、ステ
ップＳ５０に進んで、現在得られている周波数ポイント
の伝達関数を音響特性として外部の表示装置上に表示す
る。
レージングを行うことが設定されていないならば、ステ
ップＳ５０に進んで、現在得られている周波数ポイント
の伝達関数を音響特性として外部の表示装置上に表示す
る。
【００７２】また、上記ステップＳ４２において、アベ
レージングを行うことが設定されているならば、得られ
ている伝達関数を記憶しておき、ステップＳ４３に進ん
で、'ＯＫ’とされた伝達関数のデータ数とアベレージ
ングを行うために設定されたデータ数とを比較する。こ
の比較により、 'ＯＫ’とされたデータ数がアベレージ
ングを行うために設定されたデータ数に満たないときに
は、ステップＳ４４に進んでデータを追加し、ステップ
Ｓ４５でデータの値の合計を求め、データ数によるアベ
レージングを行う。そして、このアベレージングにより
得られた伝達関数データの平均値は、ステップＳ５０で
音響特性として表示される。
レージングを行うことが設定されているならば、得られ
ている伝達関数を記憶しておき、ステップＳ４３に進ん
で、'ＯＫ’とされた伝達関数のデータ数とアベレージ
ングを行うために設定されたデータ数とを比較する。こ
の比較により、 'ＯＫ’とされたデータ数がアベレージ
ングを行うために設定されたデータ数に満たないときに
は、ステップＳ４４に進んでデータを追加し、ステップ
Ｓ４５でデータの値の合計を求め、データ数によるアベ
レージングを行う。そして、このアベレージングにより
得られた伝達関数データの平均値は、ステップＳ５０で
音響特性として表示される。
【００７３】また、ステップＳ４３で、 'ＯＫ’とされ
たデータ数がアベレージングを行うために設定されたデ
ータ数と同じであるときには、ステップＳ４８に進ん
で、一番古いデータを捨てて最後に測定した新しいデー
タを記憶し、ステップＳ４９でデータの値の合計を求
め、アベレージング回数によるアベレージングを行う。
そして、このアベレージングにより得られた伝達関数デ
ータの平均値は、ステップＳ５０で音響特性として表示
される。
たデータ数がアベレージングを行うために設定されたデ
ータ数と同じであるときには、ステップＳ４８に進ん
で、一番古いデータを捨てて最後に測定した新しいデー
タを記憶し、ステップＳ４９でデータの値の合計を求
め、アベレージング回数によるアベレージングを行う。
そして、このアベレージングにより得られた伝達関数デ
ータの平均値は、ステップＳ５０で音響特性として表示
される。
【００７４】尚、上記ステップＳ４６で、アベレージン
グを行うことが設定されていると判別されるならば、ス
テップＳ４７に進んで、今まで 'ＯＫ’とされたデータ
があるか否かを判別し、 'ＯＫ’とされたデータがない
と判別されるときには、何も行わずにこの処理を終了す
る。また、 'ＯＫ’とされたデータがあると判別される
ときには、ステップＳ４５に進んでデータの値の合計を
求め、データ数によるアベレージングを行う。そして、
このアベレージングにより得られたデータの平均値は、
ステップＳ５０で音響特性として表示される。
グを行うことが設定されていると判別されるならば、ス
テップＳ４７に進んで、今まで 'ＯＫ’とされたデータ
があるか否かを判別し、 'ＯＫ’とされたデータがない
と判別されるときには、何も行わずにこの処理を終了す
る。また、 'ＯＫ’とされたデータがあると判別される
ときには、ステップＳ４５に進んでデータの値の合計を
求め、データ数によるアベレージングを行う。そして、
このアベレージングにより得られたデータの平均値は、
ステップＳ５０で音響特性として表示される。
【００７５】上述の信頼性判別処理は、各周波数ポイン
ト毎に行う。
ト毎に行う。
【００７６】尚、上記音響特性測定装置においては、複
数の測定用マイクの音響特性のアベレージングを行うこ
とにより１つの音響特性を得ているが、任意の測定用マ
イクの測定結果を音響特性測定装置における音響特性と
して選択することも可能である。
数の測定用マイクの音響特性のアベレージングを行うこ
とにより１つの音響特性を得ているが、任意の測定用マ
イクの測定結果を音響特性測定装置における音響特性と
して選択することも可能である。
【００７７】
【発明の効果】以上の説明からも明らかなように、本発
明に係る音響特性測定装置は、入力信号及び出力信号に
直交変換処理を施す直交変換手段と、上記直交変換手段
からの入力信号及び該入力信号に対応する出力信号のス
ペクトルを用いて入力信号及び出力信号の各パワースペ
クトルを算出し、また、上記入力信号及び出力信号の各
スペクトルからクロススペクトルを算出し、上記入力信
号及び出力信号の各パワースペクトル及びクロススペク
トルを用いて入力信号及び出力信号のパワースペクトル
の平均値及びクロススペクトルの平均値を求めるスペク
トル算出手段と、上記スペクトル算出手段からの入力信
号及び出力信号のパワースペクトルの平均値、及びクロ
ススペクトルの平均値を用いて各周波数成分の伝達特性
を算出する伝達特性算出手段と、上記伝達特性算出手段
からの各周波数成分の伝達特性が一定値以上であるか否
かを判別し、一定値以上の周波数成分の伝達特性のみを
平均化する平均化手段と、上記平均化手段からの複数の
測定用マイク毎の伝達特性の平均値を平均化して音響特
性を得る音響特性算出手段とを有して成ることにより、
１系統の音声再生装置の調整を容易に行うことができ
る。
明に係る音響特性測定装置は、入力信号及び出力信号に
直交変換処理を施す直交変換手段と、上記直交変換手段
からの入力信号及び該入力信号に対応する出力信号のス
ペクトルを用いて入力信号及び出力信号の各パワースペ
クトルを算出し、また、上記入力信号及び出力信号の各
スペクトルからクロススペクトルを算出し、上記入力信
号及び出力信号の各パワースペクトル及びクロススペク
トルを用いて入力信号及び出力信号のパワースペクトル
の平均値及びクロススペクトルの平均値を求めるスペク
トル算出手段と、上記スペクトル算出手段からの入力信
号及び出力信号のパワースペクトルの平均値、及びクロ
ススペクトルの平均値を用いて各周波数成分の伝達特性
を算出する伝達特性算出手段と、上記伝達特性算出手段
からの各周波数成分の伝達特性が一定値以上であるか否
かを判別し、一定値以上の周波数成分の伝達特性のみを
平均化する平均化手段と、上記平均化手段からの複数の
測定用マイク毎の伝達特性の平均値を平均化して音響特
性を得る音響特性算出手段とを有して成ることにより、
１系統の音声再生装置の調整を容易に行うことができ
る。
【００７８】ここで、上記入力信号及び出力信号のパワ
ースペクトルの平均値、及びクロススペクトルの平均値
を用いてコヒーレンスの値を計算するコヒーレンス算出
手段を有し、上記平均化手段では上記入力信号のパワー
スペクトルの平均値のレベルと上記コヒーレンス算出手
段からのコヒーレンスの値とが一定値以上である周波数
成分の伝達特性のみを平均化することにより、求める伝
達関数の信頼性を高めることができる。
ースペクトルの平均値、及びクロススペクトルの平均値
を用いてコヒーレンスの値を計算するコヒーレンス算出
手段を有し、上記平均化手段では上記入力信号のパワー
スペクトルの平均値のレベルと上記コヒーレンス算出手
段からのコヒーレンスの値とが一定値以上である周波数
成分の伝達特性のみを平均化することにより、求める伝
達関数の信頼性を高めることができる。
【００７９】また、本発明に係る音響特性測定装置は、
入力信号及び出力信号に直交変換処理を施す直交変換手
段と、上記直交変換手段からの入力信号及び該入力信号
に対応する出力信号のスペクトルを用いて入力信号及び
出力信号の各パワースペクトルを算出し、また、上記入
力信号及び出力信号の各スペクトルからクロススペクト
ルを算出し、上記入力信号及び出力信号の各パワースペ
クトル及びクロススペクトルを用いて入力信号及び出力
信号のパワースペクトルの平均値及びクロススペクトル
の平均値を求めるスペクトル算出手段と、上記スペクト
ル算出手段からの入力信号及び出力信号のパワースペク
トルの平均値、及びクロススペクトルの平均値を用いて
各周波数成分の伝達特性を算出する伝達特性算出手段
と、上記伝達特性算出手段からの各周波数成分の伝達特
性が一定値以上であるか否かを判別し、一定値以上の周
波数成分の伝達特性のみを平均化する平均化手段と、上
記平均化手段からの複数の測定用マイク毎の伝達特性の
平均値から少なくとも１つの測定用マイクの伝達特性の
平均値を独立に取り出す音響特性算出手段とを有して成
ることにより、任意の測定用マイクの測定結果を音響特
性として得ることができる。
入力信号及び出力信号に直交変換処理を施す直交変換手
段と、上記直交変換手段からの入力信号及び該入力信号
に対応する出力信号のスペクトルを用いて入力信号及び
出力信号の各パワースペクトルを算出し、また、上記入
力信号及び出力信号の各スペクトルからクロススペクト
ルを算出し、上記入力信号及び出力信号の各パワースペ
クトル及びクロススペクトルを用いて入力信号及び出力
信号のパワースペクトルの平均値及びクロススペクトル
の平均値を求めるスペクトル算出手段と、上記スペクト
ル算出手段からの入力信号及び出力信号のパワースペク
トルの平均値、及びクロススペクトルの平均値を用いて
各周波数成分の伝達特性を算出する伝達特性算出手段
と、上記伝達特性算出手段からの各周波数成分の伝達特
性が一定値以上であるか否かを判別し、一定値以上の周
波数成分の伝達特性のみを平均化する平均化手段と、上
記平均化手段からの複数の測定用マイク毎の伝達特性の
平均値から少なくとも１つの測定用マイクの伝達特性の
平均値を独立に取り出す音響特性算出手段とを有して成
ることにより、任意の測定用マイクの測定結果を音響特
性として得ることができる。
【図１】本発明に係る音響特性測定装置の概略的な構成
を示す図である。
を示す図である。
【図２】音響特性測定器の概略的な構成を示す図であ
る。
る。
【図３】音響特性測定の手順のフローチャートである。
【図４】信頼性判別処理手順のフローチャートである。
【図５】従来の音響特性測定の手順のフローチャートで
ある。
ある。
【図６】音響特性測定用信号を説明するための図であ
る。
る。
１２、２２ 窓関数発生器 １４、２４ ＦＦＴ解析器 １５、２５、３３ 複素共役変換器 ３１ａ、３１ｂ、３１ｃ レジスタ ３２ａ、３２ｂ、３２ｃ 平均値化回路 ３４ コヒーレンス演算器 ３５ 伝達関数演算器 ４２ イコライザ ４３ パワーアンプ ４４ スピーカ ４５ 音響特性測定器 ４６ 信号切換器 ４７ マイクアンプ
Claims (3)
- 【請求項１】 入力信号及び上記入力信号が音響信号に
変換され、複数の測定用マイクにより収音されて得られ
た出力信号を用いて音響特性を算出する音響特性測定装
置において、 上記入力信号及び出力信号に直交変換処理を施す直交変
換手段と、 上記直交変換手段からの入力信号及び該入力信号に対応
する出力信号のスペクトルを用いて入力信号及び出力信
号の各パワースペクトルを算出し、また、上記入力信号
及び出力信号の各スペクトルからクロススペクトルを算
出し、上記入力信号及び出力信号の各パワースペクトル
及びクロススペクトルを用いて入力信号及び出力信号の
パワースペクトルの平均値及びクロススペクトルの平均
値を求めるスペクトル算出手段と、 上記スペクトル算出手段からの入力信号及び出力信号の
パワースペクトルの平均値、及びクロススペクトルの平
均値を用いて各周波数成分の伝達特性を算出する伝達特
性算出手段と、 上記伝達特性算出手段からの各周波数成分の伝達特性が
一定値以上であるか否かを判別し、一定値以上の周波数
成分の伝達特性のみを平均化する平均化手段と、 上記平均化手段からの複数の測定用マイク毎の伝達特性
の平均値を平均化して音響特性を得る音響特性算出手段
とを有して成ることを特徴とする音響特性測定装置。 - 【請求項２】 上記入力信号及び出力信号のパワースペ
クトルの平均値、及びクロススペクトルの平均値を用い
てコヒーレンスの値を計算するコヒーレンス算出手段を
有し、上記平均化手段では上記入力信号のパワースペク
トルの平均値のレベルと上記コヒーレンス算出手段から
のコヒーレンスの値とが一定値以上である周波数成分の
伝達特性のみを平均化することを特徴とする請求項１記
載の音響特性測定装置。 - 【請求項３】 入力信号及び上記入力信号が音響信号に
変換され、複数の測定用マイクにより収音されて得られ
た出力信号を用いて音響特性を算出する音響特性測定装
置において、 上記入力信号及び出力信号に直交変換処理を施す直交変
換手段と、 上記直交変換手段からの入力信号及び該入力信号に対応
する出力信号のスペクトルを用いて入力信号及び出力信
号の各パワースペクトルを算出し、また、上記入力信号
及び出力信号の各スペクトルからクロススペクトルを算
出し、上記入力信号及び出力信号の各パワースペクトル
及びクロススペクトルを用いて入力信号及び出力信号の
パワースペクトルの平均値及びクロススペクトルの平均
値を求めるスペクトル算出手段と、 上記スペクトル算出手段からの入力信号及び出力信号の
パワースペクトルの平均値、及びクロススペクトルの平
均値を用いて各周波数成分の伝達特性を算出する伝達特
性算出手段と、 上記伝達特性算出手段からの各周波数成分の伝達特性が
一定値以上であるか否かを判別し、一定値以上の周波数
成分の伝達特性のみを平均化する平均化手段と、 上記平均化手段からの複数の測定用マイク毎の伝達特性
の平均値から少なくとも１つの測定用マイクの伝達特性
の平均値を独立に取り出す音響特性算出手段とを有して
成ることを特徴とする音響特性測定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34058194A JP3277440B2 (ja) | 1994-12-29 | 1994-12-29 | 音響特性測定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP34058194A JP3277440B2 (ja) | 1994-12-29 | 1994-12-29 | 音響特性測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08184488A true JPH08184488A (ja) | 1996-07-16 |
| JP3277440B2 JP3277440B2 (ja) | 2002-04-22 |
Family
ID=18338373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP34058194A Expired - Fee Related JP3277440B2 (ja) | 1994-12-29 | 1994-12-29 | 音響特性測定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3277440B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09258749A (ja) * | 1996-03-21 | 1997-10-03 | Nippon Hoso Kyokai <Nhk> | 音場におけるインパルス応答測定装置 |
| WO2003094576A1 (fr) * | 2002-04-30 | 2003-11-13 | Sony Corporation | Dispositif et procede de mesure de caracteristique d'emission, et amplificateur |
| US20110150230A1 (en) * | 2009-12-17 | 2011-06-23 | Canon Kabushiki Kaisha | Sound processing apparatus and method |
| JP2012141214A (ja) * | 2010-12-28 | 2012-07-26 | Mitsubishi Heavy Ind Ltd | 音量誤差測定装置及び音量誤差測定方法 |
| JP2016197030A (ja) * | 2015-04-02 | 2016-11-24 | 川崎重工業株式会社 | 音信号の信号処理装置及び方法 |
-
1994
- 1994-12-29 JP JP34058194A patent/JP3277440B2/ja not_active Expired - Fee Related
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09258749A (ja) * | 1996-03-21 | 1997-10-03 | Nippon Hoso Kyokai <Nhk> | 音場におけるインパルス応答測定装置 |
| WO2003094576A1 (fr) * | 2002-04-30 | 2003-11-13 | Sony Corporation | Dispositif et procede de mesure de caracteristique d'emission, et amplificateur |
| US7286946B2 (en) | 2002-04-30 | 2007-10-23 | Sony Corporation | Transmission characteristic measuring device transmission characteristic measuring method, and amplifier |
| US20110150230A1 (en) * | 2009-12-17 | 2011-06-23 | Canon Kabushiki Kaisha | Sound processing apparatus and method |
| US8401201B2 (en) * | 2009-12-17 | 2013-03-19 | Canon Kabushiki Kaisha | Sound processing apparatus and method |
| JP2012141214A (ja) * | 2010-12-28 | 2012-07-26 | Mitsubishi Heavy Ind Ltd | 音量誤差測定装置及び音量誤差測定方法 |
| JP2016197030A (ja) * | 2015-04-02 | 2016-11-24 | 川崎重工業株式会社 | 音信号の信号処理装置及び方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3277440B2 (ja) | 2002-04-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1954096B1 (en) | Apparatus and method for measuring loudspeaker transfer function with enhanced frequency resolution | |
| EP0159546B1 (en) | Digital graphic equalizer | |
| JPH07143034A (ja) | ハウリング抑制装置 | |
| WO2003094576A1 (fr) | Dispositif et procede de mesure de caracteristique d'emission, et amplificateur | |
| JP4440865B2 (ja) | 計測装置 | |
| JPH0795684A (ja) | 音響特性補正装置 | |
| US6453253B1 (en) | Impulse response measuring method | |
| JP3277440B2 (ja) | 音響特性測定装置 | |
| JP3381109B2 (ja) | 伝達関数測定装置 | |
| JPS63234699A (ja) | 音場補正装置 | |
| JP2867769B2 (ja) | 音響測定方法およびその装置 | |
| JP3139803B2 (ja) | インパルス応答測定装置 | |
| CN114545078A (zh) | 中频信号处理系统及处理方法 | |
| JPH055760Y2 (ja) | ||
| US5420934A (en) | Electronic sound processing system | |
| JP2589864B2 (ja) | 信号比較装置 | |
| JPH09215085A (ja) | 音響再生装置および受聴周波数特性補正方法 | |
| JP3612764B2 (ja) | 遅延時間測定方法及び装置 | |
| JP4132693B2 (ja) | イコライザ | |
| JP3455843B2 (ja) | 搬送波周波数の検出回路、検出方法及び記録媒体 | |
| JPH08294195A (ja) | 信号処理方法 | |
| KR900008161B1 (ko) | 음질조정장치 | |
| JPH0526921A (ja) | スペクトラム・アナライザ | |
| JPH0894690A (ja) | サーボアナライザ | |
| JPH09258749A (ja) | 音場におけるインパルス応答測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20020108 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080215 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090215 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100215 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |