WO2011007706A1 - インパルス応答測定方法およびインパルス応答測定装置 - Google Patents
インパルス応答測定方法およびインパルス応答測定装置 Download PDFInfo
- Publication number
- WO2011007706A1 WO2011007706A1 PCT/JP2010/061567 JP2010061567W WO2011007706A1 WO 2011007706 A1 WO2011007706 A1 WO 2011007706A1 JP 2010061567 W JP2010061567 W JP 2010061567W WO 2011007706 A1 WO2011007706 A1 WO 2011007706A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- frequency
- impulse response
- phase
- inverse filter
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/28—Measuring attenuation, gain, phase shift or derived characteristics of electric four pole networks, i.e. two-port networks; Measuring transient response
Abstract
【課題】被測定系のインパルス応答を計測する際に、送信側と受信側のサンプリングクロックが非同期であっても、簡単な装置あるいは信号処理によって精度の高いインパルス応答計測を可能にする。 【解決手段】インパルス応答測定方法は、第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成ステップと、第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換ステップと、前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の周波数特性を示す関数の逆関数である逆フィルタの位相を少なくとも補正する逆フィルタ補正ステップと、を有する。そして、補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測する。
Description
本発明は、音響機器・音響空間・電気信号の伝送線路等の被測定系の伝達特性を計測するためのインパルス応答測定方法およびインパルス応答測定装置に関し、特に、送信側と受信側の処理装置間で同期信号周波数が異なっていても正確に伝達特性を計測することができるインパルス応答測定方法およびインパルス応答測定装置に関する。
従来、音響信号を再生する音響再生装置から再生されたインパルス応答信号によって、被測定系の周波数特性等の様々な音響特性を把握する技術が広く知られており、このインパルス応答信号を精度良く測定しようとするいくつかの方法が提案されている。
例えば、従来技術の一つとして、複数(I個)の同一信号(繰返し信号)を用いて音響機器、音響空間の伝達特性を測定し、この信号の繰返し周期でそれぞれの信号を同期加算平均することによって、ランダムな雑音を持つ信号に対する信号対雑音(S/N)比を向上(10logI dBアップ)させることが行われている。この際、繰返し信号の繰返し周期に対して正確に同期加算するために、図2に示すように送信側と受信側で共通の同期信号(クロックパルス)を用いるのが一般的である。しかしながら、送信側と受信側でクロックパルスを共通にするには、送信側と受信側の測定装置の一体化、もしくはクロックパルスを送信側と受信側で共通にするための伝送路が必要である。
例えば、ホールや競技場のような大規模空間の音響測定において、空間そのものの伝達特性を測定するときには、試験用信号発生器(送信側)からの信号を音源用増幅器へ供給し、スピーカによって音を空間に放射し、受音点で受信して測定するが、広い空間では、場所による違いを把握するために、複数の測定点に亘って音を測定するのが普通である。この際、測定装置が同一のクロックパルスを用いるために、音を受信するマイクロホンにケーブルを引き回して測定するか、もしくは測定装置そのものをマイクロホンと共に測定点に移動させ、ケーブルによって、音源用増幅器を含むスピーカ側に接続する必要がある。
このため、クロックパルスを送信側と受信側で共通にしなくても正確に応答測定ができる非同期の測定システムが求められている。
本発明の発明者らは、非同期システムの受信側でA/D変換された連続する複数の各信号サンプリングデータ列を切り出し、その各データ列を離散フーリエ変換(DFT/Discrete Fourier Transform)処理をして、時間領域で同期加算するのではなく、周波数領域で同期加算する方法を見出した(特許文献1、非特許文献1参照)。
また、非同期システムにTSP(Time Stretched Pulses)法を適用し、受信側でA/D変換された複数の同一信号のサンプリングデータ列を時間領域で同期加算するために、受信側で先頭にある信号データ列と、それに続く信号データ列を切り出し、切り出した先頭にある信号データ列を基準に、それに続く信号列データ間との相互相関処理を行って相互相関値のピークを求めることによって同期加算時における各信号データ列の位置を合わせて加算をする方法が提案されている(例えば、特許文献2参照)。
この方法では、前記した相互相関処理により、受信側のクロックでサンプリングして得られた各データ列の長さ(ビット数)と送信側の各データ列の長さ(ビット数)との違いを求める。そして、インパルス応答波形を送出するために同期加算後の信号と逆TSP信号との畳み込み処理を時間領域で行うとき、逆TSP信号を再サンプリング(生成信号数を補償)して、受信側のクロックでサンプリングすることによって得られる各データ列の長さ(クロックパルス数)に対応したデータ列の長さを有する補正逆TSP信号を求めている。
また、低域の外来ノイズによる影響を最小限にするために各データ列間隔をランダムあるいは増減させる方法が提案されている(例えば、特許文献4、特許文献5参照)。
また、TSP法(例えば、特許文献3、非特許文献2参照)は、インパルス応答のデジタルによる測定に多く使われている。TSP信号は、高い周波数から低い周波数へ変化する信号、あるいは低い周波数から高い周波数へ変化する信号(周波数をスイープ(掃引)する信号)であって、インパルスを時間軸上で引き延ばすことによりエネルギーを増大させている。したがって、高いS/N比が得られるインパルス応答測定ができる。近年、TSP法の改善がさらに試みられ、低域のパワーをさらに上げられるようなTSPフィルタが提案されている(例えば、特許文献6、特許文献7、特許文献8参照)。
日野、土屋、遠藤,"非同期な測定系で同期加算を行う方法の検討",第10回AESリージョナルコンベンション東京 予稿集,2001年6月
大賀寿郎、山崎芳雄、金田豊,"音響システムとデジタル処理",電子情報通信学会誌p.158-159,1995年
ところで、たとえばCDプレーヤを使った計測方式あるいは装置では、CDプレーヤは、ディスクに記録されたTSP情報を検出し、D/A変換器によりアナログ変換し、アンプやスピーカを介して被測定室内で放音する。そして、受信側のコンピュータでは、マイクを介してスピーカから放音された音を取り込み、A/D変換器によってデジタル変換し、信号処理を実行し、インパルス応答測定をすることになる。
しかしながら、上記CDプレーヤ側(送信側)のD/A変換器の同期信号であるサンプリングクロックと、コンピュータ(受信側)のA/D変換器のサンプリングクロックの周波数(あるいは位相)が異なっているとき、つまり、送信側と受信側でサンプリングクロックが非同期である場合には、同じ波形を表すサンプリング数(標本化数)が異なり、測定される信号に誤差が生じるので、正確なインパルス応答を求めることができない。
また、上述の従来技術(例えば、特許文献2に記載の技術)では、送信時のTSP信号データ列に対する逆TSP信号から再サンプリングすることによって補正逆TSP信号データ列を生成することは容易でない。たとえば、大容量の波形メモリに時間軸に対して高分解能な逆TSP信号データ列を記憶しておき、送信側と受信側の同期信号周波数比(サンプリングクロックパルス周波数比)に応じて、順次、必要な逆TSP信号データ列を取り出して補正逆TSP信号データ列を生成するといった複雑な処理が必要である。
本発明は、このような課題を解決するためになされたものであって、被測定系のインパルス応答を計測する際に、送信側と受信側のサンプリングクロックが非同期であっても、簡単な装置あるいは信号処理によって精度の高いインパルス応答計測ができるインパルス応答測定方法およびインパルス応答測定装置を提供することを目的とする。
本発明は、第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成ステップと、第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換ステップと、前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の周波数特性で表現される関数の逆関数である逆フィルタの位相を少なくとも補正する逆フィルタ補正ステップと、を有し、補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測することを特徴とする、インパルス応答測定方法である。
また、本発明は、第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成手段と、第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換手段と、前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の生成フィルタの逆フィルタの位相を少なくとも補正する逆フィルタ補正手段と、を有し、補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測することを特徴とする、インパルス応答測定装置である。
本発明に係るインパルス応答測定方法およびインパルス応答測定装置によれば、被測定系のインパルス応答を計測する際に、送信側と受信側のサンプリングクロックが非同期であっても、簡単な装置あるいは信号処理によって精度の高いインパルス応答計測ができる。
また、ホールや競技場のような大規模空間の音響測定において、試験用信号発生器からの信号を音源用増幅器へ供給しスピーカによって音を空間に放射し、受音点で測定装置によって受信して測定する場合、敢えて試験用信号発生器と測定装置に同一の同期信号(クロックパルス)を使わなくても正しいインパルス応答計測ができる。このため、音を受信するマイクロホンにケーブルを引き回して音を測定したり、もしくは測定装置そのものをマイクロホンと共に測定点に移動させ、ケーブルによって音源用増幅器を含むスピーカ側に接続したりする必要がない。
また、すでに車に搭載された音響機器から、人が聴く位置(運転位置)までのインパルス応答特性を、その音響機器とはケーブルで接続されていない計測器によって計測することができる。たとえば、CDに同期加算平均する繰り返しのTSPデータを記録し、それを音響機器であるCDプレーヤで再生することによって音のインパルス応答特性を計測器で容易かつ高精度に測定することが可能となる。
以下、図面を用いて、本実施形態に係るインパルス応答測定方法およびインパルス応答測定装置について詳細に説明する。
<1.全体構成>
まず、図1を用いて、非同期系システムにおけるTSP法によるインパルス応答測定の処理の流れについて説明する。なお、非同期系システムにおける複数の同一の信号(繰返し信号)の高精度な周期の検出法と、周波数領域での同期加算平均である同期ベクトル加算平均については後述する。
まず、図1を用いて、非同期系システムにおけるTSP法によるインパルス応答測定の処理の流れについて説明する。なお、非同期系システムにおける複数の同一の信号(繰返し信号)の高精度な周期の検出法と、周波数領域での同期加算平均である同期ベクトル加算平均については後述する。
同図に示すように、非同期系システム1は、互いに異なるサンプリングクロック周波数(fsとf's)で動作する、送信側ユニット2と受信側ユニット3を備える。
<1-1.送信側ユニット>
送信側ユニット2では、周波数領域でのインパルス入力I(k)を、伝達関数H(k)のTSPフィルタ4に入力することによって、TSPフィルタ4の出力であるH(k)I(k)が得られる。次に、このTSPフィルタ4の出力であるH(k)I(k)に対して、逆離散フーリエ変換処理部5でIDFT(Inverse Discrete Fourier Transform)を実行することによって、送信側ユニット2のデータのサンプリングクロック周波数fsに同期したサンプル列であるTSP信号h(n)が得られる。次に、このTSP信号h(n)をD/A変換器6に通すことによってTSP信号h(t)が得られる。そして、このTSP信号h(t)を、インパルス応答特性g(t)である空間などからなる測定対象(被測定系)に入力することによって、被測定信号の出力x(t)=h(t)*g(t)(*は畳み込み演算を表す)が得られる。
送信側ユニット2では、周波数領域でのインパルス入力I(k)を、伝達関数H(k)のTSPフィルタ4に入力することによって、TSPフィルタ4の出力であるH(k)I(k)が得られる。次に、このTSPフィルタ4の出力であるH(k)I(k)に対して、逆離散フーリエ変換処理部5でIDFT(Inverse Discrete Fourier Transform)を実行することによって、送信側ユニット2のデータのサンプリングクロック周波数fsに同期したサンプル列であるTSP信号h(n)が得られる。次に、このTSP信号h(n)をD/A変換器6に通すことによってTSP信号h(t)が得られる。そして、このTSP信号h(t)を、インパルス応答特性g(t)である空間などからなる測定対象(被測定系)に入力することによって、被測定信号の出力x(t)=h(t)*g(t)(*は畳み込み演算を表す)が得られる。
<1-2.受信側ユニット>
受信側ユニット3では、被測定系を通して得られた被測定信号x(t)(=h(t)*g(t))を、サンプラー(あるいはA/D変換器)7において、受信側ユニット3のデータのサンプリングクロック周波数f's=fs/α(αは周波数比を表す)でサンプリングすることによって、受信側ユニット3のサンプリングクロック周波数f'sに同期したサンプル列である出力信号x(m)が得られる。次に、この出力信号x(m)に対して、DFT処理部8でDFTを実行することによって、X(l)が得られ、このX(l)とTSP逆フィルタH-1(l)9の伝達特性との乗算H-1(l)X(l)を行うことによってG'(l)が得られる。次に、このG'(l)に対してIDFT処理部10でIDFTを実行することによって、サンプル列であるg(m)が得られ、このg(m)をD/A変換器11に通すことによってインパルス応答波形g(t)が得られる。なお、この例で示したD/A変換器6、11には、入力される信号をD/A変換し、アナログゲート等でパルス列を作り、アナログ低域フィルタで補間を行う機能が含まれている。
受信側ユニット3では、被測定系を通して得られた被測定信号x(t)(=h(t)*g(t))を、サンプラー(あるいはA/D変換器)7において、受信側ユニット3のデータのサンプリングクロック周波数f's=fs/α(αは周波数比を表す)でサンプリングすることによって、受信側ユニット3のサンプリングクロック周波数f'sに同期したサンプル列である出力信号x(m)が得られる。次に、この出力信号x(m)に対して、DFT処理部8でDFTを実行することによって、X(l)が得られ、このX(l)とTSP逆フィルタH-1(l)9の伝達特性との乗算H-1(l)X(l)を行うことによってG'(l)が得られる。次に、このG'(l)に対してIDFT処理部10でIDFTを実行することによって、サンプル列であるg(m)が得られ、このg(m)をD/A変換器11に通すことによってインパルス応答波形g(t)が得られる。なお、この例で示したD/A変換器6、11には、入力される信号をD/A変換し、アナログゲート等でパルス列を作り、アナログ低域フィルタで補間を行う機能が含まれている。
このように、受信側ユニット3では、複数の各TSP信号によるインパルス応答出力x(m)をDFT処理し、周波数領域での同期ベクトル加算平均を実行した後、TSPフィルタの逆関数であるTSP逆フィルタH-1(l)の伝達特性を乗ずるのであるが、この逆関数に対して、サンプリング周波数ずれの影響を考慮するとより高精度のインパルス応答測定ができる。以下、周波数領域での同期ベクトル加算平均と、TSP逆フィルタH-1(k)の補正原理について詳細に説明する。
<2.周波数領域での同期ベクトル加算平均>
まず、周波数領域での同期ベクトル加算平均について詳細に説明する。
まず、周波数領域での同期ベクトル加算平均について詳細に説明する。
(ステップ1)受信側ユニット3のA/D変換器7より、受信側ユニット3のサンプリングクロック周波数f'sに同期して出力される同期加算平均対象の全信号データ列x(m)を読み出す。このとき、送信側ユニット2で決定されているデータフォーマットに基づいて、同期加算平均対象のデータ列の前に設けた同期データを基準にサンプリングデータ列が連続しているとして複数の各信号データ列を切り出す。
(ステップ2)先頭信号データ列とそれに続く信号データ間の相互相関を取り、相関値がピークとなる時間軸上の各信号データ列のサンプリングクロック数とサンプリングクロック間の位相で表現される位置誤差(位相差)あるいは同期位置を推定する。このような相互相関による各データ列の位置計測は、受信側(受信側ユニット3)のデータ列を長くできるTSP法によるインパルス応答計測には好適であり、データ列を長くすれば相互相関の値はノイズに影響されにくくなる。これは、信号とノイズ間及び一方のデータ列にのっているノイズと他方に乗っているノイズ間に相関がないからである。
(ステップ3)相互相関値の正確なピーク点のもつ位置を求めるために相関形状の対称性を利用して、二次関数による補間を行う。これにより、同期位置の位相情報を正確に知ることができる。
(ステップ4)切り出された各データ列に対してDFT処理部8においてDFTを実行する。受信側ユニット3のサンプリングクロック周波数が送信側ユニット2のサンプリングクロック周波数と異なるときは、切り出した各信号データ列の位置の位相はばらついている。このため、これら位置誤差(時間ずれ)に相当する周波数領域の各周波数に対する位相ずれ分を先に求めたDFTデータに対して補正し、補正したDFTデータについて周波数ごとの複素ベクトル加算平均を行う。これが同期ベクトル加算平均である。
ここで、位置誤差(時間ずれ)は周波数領域では位相の回りに対応している。時間ずれτと周波数領域での位相θとの関係について数式的に表現しておく。周波数領域での角周波数をω、各角周波数での位相回転量をθ(ω)とすると以下の式が成り立つ。
<数式1>
θ(ω)=-ωτ、ω=2πf 、τ=[CN+{Cφ/(2π)}]/f's ・・・ (1)
ただし、f's:サンプリングクロック周波数
CN:位置誤差時間τ内に存在するサンプリングクロック数
Cφ:位置誤差時間τ間位相を2πCNで測定しきれない分のサンプリングクロック間位相
この値は、数値の表現範囲での分解能を持つことができるので、データ量を増加せずに精度を確保できる。この特長により、切り出した各信号データ列間の時間ずれ量に見合うようDFTした信号データの周波数毎の位相回転量を加減した上で、信号データ列全体における周波数毎にベクトル加算と平均処理がされたDFT(伝達関数)の結果X(l)を得ることができる。
<数式1>
θ(ω)=-ωτ、ω=2πf 、τ=[CN+{Cφ/(2π)}]/f's ・・・ (1)
ただし、f's:サンプリングクロック周波数
CN:位置誤差時間τ内に存在するサンプリングクロック数
Cφ:位置誤差時間τ間位相を2πCNで測定しきれない分のサンプリングクロック間位相
この値は、数値の表現範囲での分解能を持つことができるので、データ量を増加せずに精度を確保できる。この特長により、切り出した各信号データ列間の時間ずれ量に見合うようDFTした信号データの周波数毎の位相回転量を加減した上で、信号データ列全体における周波数毎にベクトル加算と平均処理がされたDFT(伝達関数)の結果X(l)を得ることができる。
(ステップ5)同期加算した結果をTSP逆フィルタH-1(l)9と乗算を行い、その結果をIDFT処理部10においてIDFTすることによって、インパルス応答波形g(m)を得る。
<3.TSP逆フィルタH-1(k)の補正原理>
次に、TSP逆フィルタH-1(k)の補正原理について詳細に説明する。
次に、TSP逆フィルタH-1(k)の補正原理について詳細に説明する。
受信側ユニット3の応答波形のサンプル列x(m)は、送信側ユニット2のサンプリング周波数とは異なる周波数でサンプリングされているので、DFT変換したX(l)にTSP逆フィルタH-1(l)を乗じてその結果をIDFTしてインパルス応答のパルス列を求めるためには、TSP逆フィルタH-1(l)が受信側ユニット3のサンプリング周波数に対応して生成されていなければならない。
送信側ユニット2のサンプリングクロック周波数fsと受信側ユニット3のサンプリングクロック周波数f'sとの比をα(=fs/f's)とする。
図3において受信側ユニット3では、α>1のときは、送信側ユニット2と同一のサンプル総数Nで処理しようとした場合、検出信号幅が送信側ユニット2ではサンプル総数がWN個であったものが受信側ユニット3ではWNS=〔WN/α〕と縮小されてしまう(ここで、〔A〕はAを整数に丸め処理したことを表す)。つまり、受信側ユニット3を基準とする測定では波形が時間軸方向に縮小したように見える。
応答波形のサンプル列x(m)にDFTを行いX(l)を得て、これにTSPフィルタH(k)4の逆関数H-1(l)を乗ずるとき、TSPフィルタH(k)4をサンプリング周波数の変化分だけ補正して得られたTSP逆フィルタH-1(l)9を使う必要がある。言い換えれば、受信側ユニット3では、自身のサンプリングクロック周波数f'sに基づいて処理するので、あたかも検出波形が縮小したように見える。したがって、ここでは検出側での絶対時間は変わらないとし、サンプリング周波数の変化によって波形が伸縮したとして扱う。
このときの周波数領域との関係を見るために、受信側ユニット3で得られたサンプル列x(m)を連続なx(αt)とみなして、x(t)のフーリエ変換X(f)における時間のスケールファクタαの関係を示す次の(2)式の公式(E.Oran Brigham,"高速フーリエ変換" Prentice-Hall Inc,第3章(1974)参照)を利用して考察する。スケールファクタαは、先に述べた送信側のサンプリング周波数fsと受信側のサンプリング周波数f'sとの比である。
<数式2>
x(αt)⇔(1/|α|)X(f/α) ・・・ (2)
これは定性的に言えば、α>1とすると、時間領域の波形は時間軸側に縮小した形になるが、周波数領域では周波数軸上で拡大し、振幅は小さくなる。
<数式2>
x(αt)⇔(1/|α|)X(f/α) ・・・ (2)
これは定性的に言えば、α>1とすると、時間領域の波形は時間軸側に縮小した形になるが、周波数領域では周波数軸上で拡大し、振幅は小さくなる。
図4は、連続系と離散系における時間のスケールファクタαに関する関係を模式的に描いた図であり、上段と中段には連続系に関する式(2)の関係を示し、下段には離散系DFTの関係を示している。
連続系において測定時間Tで時間長Wt及び周波数帯域-Fm~+Fmをもつ信号x(t)の振幅周波数特性と位相周波数特性を上段に示し、スケールファクタα>1の信号x(αt)の振幅周波数特性と位相周波数特性を中段に示す。後者のスケールファクタα>1の信号x(αt)では、測定時間Tで時間長がWt/αとなる波形となり、その周波数帯域が-αFm~+αFmで、振幅周波数特性の振幅が1/αとなる。
このスケールファクタαの連続系の信号x(αt)をDFTするとき、信号のx(αt)が周期Tで繰り返され、周波数領域の帯域幅-Fx~+Fxは、Fx>αmaxFm(ただし、αmaxはα>0の最大値)の条件で設定する。DFTであるので周波数領域で基底周波数帯域が-Fx~+Fxとなる周期2Fxの周期関数となる。
このとき、時間軸のサンプル総数Nは高速フーリエ変換(FFT/Fast Fourier Transform)演算を考慮してN=2i(ただしiは整数)とする。時間軸のサンプリング周波数f'sは2Fx以上であれば良いが、ここではf's=2Fxと設定するとともに、時間周期TをT=N/f'sと設定する。このとき周波数領域では、信号の周波数特性として1/Tごとの線スペクトルが現れ、帯域-αFm~+αFmの特性が周波数軸上において周期2Fxで繰り替えされる。2Fx/(1/T)=Nとなるので、周波数領域では、基底離散周波数帯域が-N/2~+N/2-1となる周期Nの周期関数である。離散信号列が周期Tで繰り返すとして扱っているので、DFT後の周波数領域では1/Tごとの線スペクトルが現れる。
ここまでは、スケールファクタαがあるときの離散時間の取り扱いを示した。また、図4は模式的に示した図でありTSP信号による実際の応答波形x(t)をかなり簡略して描いている。また、TSPフィルタH(k)4の出力H(k)I(k)をIDFT処理部5でIDFTした信号列h(n)は、そのままでは先頭部分には信号が現れないので、図5に示す環状シフトを行ったことにしている。つまり先頭の信号列を後方にシフトしているが、DFTは周期性を前提としているのでこのようにしても結果は変わらない。
次に、受信側ユニット3の応答波形のサンプル列x(m)が、送信側ユニット2のサンプリングクロック周波数とは異なる周波数でサンプリングされているときの補正されたTSP逆フィルタH-1(l)を導出する。
ここで、TSPフィルタ(伝達関数)あるいは離散的フーリエ変換(DFT)H(k)は以下の式が提案されている(例えば、特許文献3参照)。
<数式3>
H(k)=exp(jβk2) 0≦k≦M/2(Mは整数)
H*(M-k)=exp{-jβ(M-k)2} (M/2)+1<k<M
β(M/2)2=2Lπ⇔β=8Lπ/M2(Lは整数)
・・・ (3)
この式においてN≧αmaxMの条件が成り立つN点DFTを考慮し、かつ式(3)を参考にして、以下の式(4)に変形する。図4のように周期Nの繰り返し特性となる。
<数式4>
H(k)=exp(jβk2) 0≦k≦M/2-1
H*(k)=exp(-jβk2) (M/2)≦k<0
H(k)=0 -N/2≦k≦-(M/2+1)及び(M/2)≦k≦N/2-1
β=8Lπ/M2(Lは整数)
・・・ (4)
上記フーリエ変換対の式x(αt)⇔(1/|α|)X(f/α)より検出波形の時間軸がα倍されたとするとフーリエ変換した周波数軸は1/α倍となる。図4を用いて考察したスケールファクタαのDFTへの展開を踏まえると、送信側ユニット2のTSPフィルタH(k)と受信側ユニット3のTSPフィルタH(l)の位相θの特性は図6のようになる。なお、同図はα>1として描いている。
ここで、TSPフィルタ(伝達関数)あるいは離散的フーリエ変換(DFT)H(k)は以下の式が提案されている(例えば、特許文献3参照)。
<数式3>
H(k)=exp(jβk2) 0≦k≦M/2(Mは整数)
H*(M-k)=exp{-jβ(M-k)2} (M/2)+1<k<M
β(M/2)2=2Lπ⇔β=8Lπ/M2(Lは整数)
・・・ (3)
この式においてN≧αmaxMの条件が成り立つN点DFTを考慮し、かつ式(3)を参考にして、以下の式(4)に変形する。図4のように周期Nの繰り返し特性となる。
<数式4>
H(k)=exp(jβk2) 0≦k≦M/2-1
H*(k)=exp(-jβk2) (M/2)≦k<0
H(k)=0 -N/2≦k≦-(M/2+1)及び(M/2)≦k≦N/2-1
β=8Lπ/M2(Lは整数)
・・・ (4)
上記フーリエ変換対の式x(αt)⇔(1/|α|)X(f/α)より検出波形の時間軸がα倍されたとするとフーリエ変換した周波数軸は1/α倍となる。図4を用いて考察したスケールファクタαのDFTへの展開を踏まえると、送信側ユニット2のTSPフィルタH(k)と受信側ユニット3のTSPフィルタH(l)の位相θの特性は図6のようになる。なお、同図はα>1として描いている。
式(4)においてkにl/αを代入し、振幅を1/|α|としたDFTのH(l)を求める。
<数式5>
H(l)=(1/|α|)exp{jβ(l/α)2}
=(1/|α|)exp{j(β/α2 )l2}, 0≦l≦〔αM/2〕-1
H*(l)=(1/|α|)exp{-j(β/α2)l2}, -〔αM/2〕≦l<0
H(l)=0,-N/2≦l≦-〔αM/2〕-1及び〔αM/2〕≦l≦N/2-1
β=8Lπ/M2
・・・ (5)
ただし、〔αM/2〕はαM/2を計算し整数に丸めた値を示す。丸めはたとえば四捨五入をする。
<数式5>
H(l)=(1/|α|)exp{jβ(l/α)2}
=(1/|α|)exp{j(β/α2 )l2}, 0≦l≦〔αM/2〕-1
H*(l)=(1/|α|)exp{-j(β/α2)l2}, -〔αM/2〕≦l<0
H(l)=0,-N/2≦l≦-〔αM/2〕-1及び〔αM/2〕≦l≦N/2-1
β=8Lπ/M2
・・・ (5)
ただし、〔αM/2〕はαM/2を計算し整数に丸めた値を示す。丸めはたとえば四捨五入をする。
l=-〔αM/2〕のときβ=8Lπ/M2をH*(l)に代入すると、H*(l)=(1/|α|)exp{-j(2Lπ)}となり、図6における位相特性は、l=-〔αM/2〕での位相が-2Lπとなることがわかる。
H(l)の逆関数あるいは逆フィルタH-1(l)は以下となる。
<数式6>
H-1(l)=|α|exp{-j(β/α2)l2}, 0≦l≦〔αM/2〕-1
H*-1(l)=|α|exp{j(β/α2)l2}, -〔αM/2〕≦l<0
H-1(l)=0,-N/2≦l≦-〔αM/2〕-1及び〔αM/2〕≦l≦N/2-1
β=8Lπ/M2
α=fs/f's
fs:送信側のサンプリング周波数
f's:受信側のサンプリング周波数
・・・ (6)
なお、ここでは被測定系のインパルス応答の時間長TLは、図4における受信側の周期T以下であるとしている。
<数式6>
H-1(l)=|α|exp{-j(β/α2)l2}, 0≦l≦〔αM/2〕-1
H*-1(l)=|α|exp{j(β/α2)l2}, -〔αM/2〕≦l<0
H-1(l)=0,-N/2≦l≦-〔αM/2〕-1及び〔αM/2〕≦l≦N/2-1
β=8Lπ/M2
α=fs/f's
fs:送信側のサンプリング周波数
f's:受信側のサンプリング周波数
・・・ (6)
なお、ここでは被測定系のインパルス応答の時間長TLは、図4における受信側の周期T以下であるとしている。
このように受信側のサンプリング周波数がずれたときでも、式(6)のTSP逆フィルタH-1(l)を用いれば正確なインパルス応答列g(m)を得ることができる。また、補正TSP逆フィルタH-1(l)は、数値の表現範囲での分解能で補正できるので、データ量を増加せずに精度を確保できる。
<4.複数データ列の構成>
ここで、複数の同期加算データ列が連続している場合と、各データ列間隔が一定でない場合について検討する。特に後者は、低域の外来ノイズに対して強くするために提案されている(特許文献4.特許文献5参照)。
ここで、複数の同期加算データ列が連続している場合と、各データ列間隔が一定でない場合について検討する。特に後者は、低域の外来ノイズに対して強くするために提案されている(特許文献4.特許文献5参照)。
(1)複数のデータ列が連続している場合(図7a参照)
このときのスケールファクタαは、データ列全体の自己相関を取るかあるいは隣接するデータ列の相互相関を取れば周期が求まるので、これにより送信側の予め決められた周期と比較して求められる。前記相互相関処理によって周期を求めるとき補間処理をして精度を上げる方法は、上記した同期加算平均を取るときと同様である。同期加算する先頭のデータ列は回路系の過渡応答等の影響で波形が安定していない場合は、同期加算平均データ列としては使わない。
このときのスケールファクタαは、データ列全体の自己相関を取るかあるいは隣接するデータ列の相互相関を取れば周期が求まるので、これにより送信側の予め決められた周期と比較して求められる。前記相互相関処理によって周期を求めるとき補間処理をして精度を上げる方法は、上記した同期加算平均を取るときと同様である。同期加算する先頭のデータ列は回路系の過渡応答等の影響で波形が安定していない場合は、同期加算平均データ列としては使わない。
(2)各データ列が不等間隔の場合(図7b参照)
この場合は、同期加算平均するデータの先端部が過渡応答の影響を受けないようにするため、あるいはその後端部が切り出したときにデータを失わないようにする必要がある。そのために、2つ以上のN点信号データブロックを一塊としたデータブロックとし、このデータブロックの間隔を変えてデータ列を構成する。過渡応答あるいは切り出し時のデータ損失を伴わないのであれば、1つのN点信号データブロックを一塊としたデータブロック構成でも良い。
この場合は、同期加算平均するデータの先端部が過渡応答の影響を受けないようにするため、あるいはその後端部が切り出したときにデータを失わないようにする必要がある。そのために、2つ以上のN点信号データブロックを一塊としたデータブロックとし、このデータブロックの間隔を変えてデータ列を構成する。過渡応答あるいは切り出し時のデータ損失を伴わないのであれば、1つのN点信号データブロックを一塊としたデータブロック構成でも良い。
データの切り出しは、同期加算対象のデータ列の前に設けた同期データを基準にサンプリングデータ列が予め決められたフォーマットで連続しているとして各信号データ列を切り出す。そして相互相関によって切り出した各データ列の位相誤差を求める。また隣接するデータ列の相互相関を取ることにより、データ列の周期を測定し、その周期と送信側で予め決められたデータ列の周期と比較して、スケールファクタαを求める。
このように、本実施形態に係るインパルス応答計測方法およびインパルス応答計測装置によれば、非同期システムにおいても、同期加算するデータ列が連続でも、あるいはそれらの間隔が一定でなくても精度の高い同期加算が実行できる。
以上説明したように、本実施形態に係るインパルス応答測定装置は、第1のサンプリングクロック周波数(例えば、fs)を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成手段(例えば、送信ユニット2のD/A変換器6など)と、第2のサンプリングクロック周波数(例えば、f's)を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換手段(例えば、受信ユニット3のA/D変換器7やDFT変換部8)と、前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比(例えば、α)に応じて前記入力信号の生成フィルタの逆フィルタ(例えば、TSP逆フィルタH-1(k))の少なくとも位相を補正する逆フィルタ補正手段と、を有し、補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測することを特徴とする、インパルス応答測定装置である。なお、「入力信号の生成フィルタの逆フィルタ」とは、入力信号の離散的周波数特性で表現される関数の逆関数となる伝達関数のことである。
本実施形態に係るインパルス応答測定装置(またはインパルス応答測定方法)によれば、被測定系のインパルス応答を計測する際に、送信側と受信側のサンプリングクロックが非同期であっても、簡単な装置あるいは信号処理によって精度の高いインパルス応答計測ができる。
また、ホールや競技場のような大規模空間の音響測定において、試験用信号発生器からの信号を音源用増幅器へ供給しスピーカによって音を空間に放射し、受音点で測定装置によって受信して測定する場合、敢えて試験用信号発生器と測定装置に同一の同期信号(クロックパルス)を使わなくても正しいインパルス応答計測ができる。このため、音を受信するマイクロホンにケーブルを引き回して音を測定したり、もしくは測定装置そのものをマイクロホンと共に測定点に移動させ、ケーブルによって音源用増幅器を含むスピーカ側に接続したりする必要がない。
さらに、すでに車に搭載された音響機器から、人が聴く位置(運転位置)までのインパルス応答特性を、その音響機器とはケーブルで接続されていない計測器によって計測することができる。たとえば、CDに同期加算平均する繰り返しのTSPデータを記録し、それを音響機器であるCDプレーヤで再生することによって音のインパルス応答特性を計測器で容易かつ高精度に測定することが可能となる。
なお、前記入力信号生成手段は、同じ任意波形の前記入力信号を等間隔あるいは不等間隔で繰返し発生する信号発生器、または前記入力信号と同じ信号が記録された媒体(例えば、図8のCD12)と前記入力信号を繰返し再生する再生器(例えば、CDプレーヤ13)を測定用信号源として用い、前記測定用信号源から発生される繰返し信号を前記入力信号として前記被測定系に入力し、前記信号変換手段は、前記被測定信号を受信点で受信し、前記測定用信号源と前記受信点の間で共通の同期信号を使用せずに、前記被測定信号の波形から得られた各周期の波形情報を使用して前記被測定信号を切り分け、切り分けた各周期の前記波形情報の時間のずれ量は基準とする周期と他の周期の相互相関の相関値が真の最大値となる時間のずれ量により求め、前記時間ずれの補正には各周期の波形を周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、先に変換し位相が補正された各周期の振幅と位相の情報をベクトル量として加算平均し、前記逆フィルタ補正手段は、前記被測定信号の自己相関あるいは繰り返す信号の隣接する信号間の相互相関により求められる前記被測定信号の信号周期と前記繰返し信号の周期とを比較して前記周波数比を求め、前記測定用信号源の繰返し信号の周波数領域における前記逆フィルタの位相を前記周波数比に応じて補正し、前記信号変換ステップによる加算平均の結果と、前記逆フィルタ補正ステップによる補正後の逆フィルタとの積を算出し、この算出結果を時間領域に変換して前記インパルス応答を計測してもよい。
このような構成とすれば、送信側と受信側のサンプリングクロックが非同期であっても、簡易で精度の高い同期加算平均による雑音低減手法を利用することができるので、インパルス応答測定において、現場の作業性を機動性のあるものにし、測定方法や測定装置を軽便かつ簡易にすることができる上に、インパルス応答測定の精度をさらに向上させることができる。
また、前記信号変換手段は、前記時間ずれの補正に各周期を離散的フーリエ変換(DFT)して前記周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、各周期の前記波形情報の和は、各周波数の複素ベクトル量を加算平均してもよい。
また、前記逆フィルタ補正手段に用いられるスケールファクタαは、前記自己相関または前記相互相関による相関値の最大値を補間して求めてもよい。
本発明に係るインパルス応答計測方法の一例として、すでに車に搭載された音響機器から、人が聴く位置(運転位置)までの応答特性を、音響機器と計測器とをケーブル接続しないで計測する方法について説明する。
図8に示すように、このインパルス応答計測方法では、同期加算平均する繰り返しのTSPデータをCDに記録し、このCDを音響機器であるCDプレーヤで再生し、この再生音を計測器(オーディオサウンドアナライザ16(以下、ASA16と言う場合がある))で測定する。
上記図1を併せて参照し、本実施例をより詳細に説明する。CD12には、複数個の同一のTSPデータ列h(n)が連続して記録されている。本実施例では、CD12に記録されたデータをCDプレーヤ13によって読み出し、読み出した信号h(t)をスピーカ14に供給する。
スピーカ14より出力された音は、マイクロホン15によって電気信号に変換された後に、ASA16に入力される。これにより、ASA16では、検出アンプ出力(応答波形)x(t)を得ることになる。なお、ASA16が得た応答波形x(t)には、空間の特性とスピーカ14の特性が含まれている。ASA16は、応答波形の連続した複数のデータ列x(m)を、USBなどのインターフェースを介してPC(パーソナルコンピュータ)17へ転送する。
PC17は、CD12に記録されたデータフォーマットに基づいて、同期加算対象のデータ列の前に設けた同期データを基準にサンプリングデータ列が連続しているとして複数の各信号データ列を切り出し、相互相関処理によって各データ列x(m)間の位置誤差(時間差、および自己相関によって繰り返しの信号周期を求めて入力側とのサンプリング周波数比であるスケールファクタαを求める。続いて、PC17は、各データ列x(m)に対してDFT処理部8でDFTを実行して各X(l)を求め、先に求めたデータ位置情報によって補正位相量を求めて各X(l)の周波数ごとの位相を補正する。続いて、PC17は、複数ある位相補正された各X(l)に対して同期ベクトル加算を行う。最後に、PC17は、スケールファクタαで補正されたTSP逆フィルタH-1(l)関数9を求め、同期ベクトル加算平均結果のX(l)との乗算G'(l)=H-1(l)X(l)を実行し、その結果を、IDFT処理部10で高速逆フーリエ変換(IFFT/Inverse Fast Fourier Transform)処理し、インパルス応答列g(m)を生成し、表示部にこの結果を示す。
本実施例では、送信側のデータフォーマットとして、同期加算対象のデータ列の前に同期データ列を設けたので、受信側では、受信信号のデータ中の同期データ列と同期を取ることによって同期位置が判定できる。つまりデータの存在場所を知ることができる。しかしながら、同期データを設けないで最初にデータが現れた位置を基準としてデータ処理をしても良い。
また、同期データ列と受信側で得られた同期加算対象の全信号データ列を最初に同期加算する単位に切り出すために、本実施例では、送信側のデータフォーマットが既知であることを利用して連続するデータ列を切り出した。しかしながら、同期加算平均対象の同一の繰り返しデータ列が連続しているのであれば、自己相関を実行してデータの繰り返し周期を求めて各データ列が切り出せる。
次に、本実施形態に係るTSP逆フィルタH-1(l)の有効性を音響空間など被測定対象のない状態でのシミュレーションによって、その効果を確認する。そのために送信側のサンプリング周波数がずれて、受信側の波形が時間方法にずれたとしたときに、TSP逆フィルタを補正したときとしないときで、インパルス波形がどうなるかを調べた。
なお、シミュレーションは以下の条件で実施した。同期加算回数は5回、サンプリング周波数はfs=44.1kHz、式(4)および式(6)に使った定数は、Nが65536、Lが8192、Mが65536×α、αが0.999(-0.1%)、Tが約1.5秒である。このシミュレーション結果を図9に示す。また、このときのTSP波形h(t)を図10に示す。
図9aに示すように、TSP逆フィルタH-1(l)を補正しないとインパルス波形とは程遠い波形となってしまうが、図9bに示すように、TSP逆フィルタH-1(l)を補正すれば、非同期システムの送信側と受信側のサンプリング周波数が異なっても精度の高いインパルス波形が得られることがわかる。
なお、本発明に係るインパルス応答測定方法およびインパルス応答測定装置の構成は、上記実施形態で示した構成に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
したがって、例えば、上記実施形態では、インパルス応答測定の処理において、周波数領域での同期ベクトル加算平均を実行する例を示したが、例えば、インパルス応答測定にそれほど高い精度が要求されない場合には、周波数領域での同期ベクトル加算平均を実行せずにインパルス応答測定を行ってもよい。
すなわち、本発明に係るインパルス応答測定方法およびインパルス応答測定装置は、第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成し、第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行い、前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の逆フィルタの位相を少なくとも補正し、補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測するものであればよい。
本発明に係るインパルス応答測定方法およびインパルス応答測定装置は、音響機器・音響空間・電気信号の伝送線路等の被測定系の伝達特性を計測するために利用することができる。
1 非同期系システム
2 送信側ユニット
3 受信側ユニット
4 TSPフィルタ
5 IDFT処理部
6、11 D/A変換器
7 サンプラー(あるいはA/D変換器)
8 DFT処理部
9 TSP逆フィルタH-1(l)
10 IDFT処理部
12 CD
13 CDプレーヤ
14 スピーカ
15 マイクロホン
16 オーディオサウンドアナライザ(ASA)
17 PC
2 送信側ユニット
3 受信側ユニット
4 TSPフィルタ
5 IDFT処理部
6、11 D/A変換器
7 サンプラー(あるいはA/D変換器)
8 DFT処理部
9 TSP逆フィルタH-1(l)
10 IDFT処理部
12 CD
13 CDプレーヤ
14 スピーカ
15 マイクロホン
16 オーディオサウンドアナライザ(ASA)
17 PC
Claims (8)
- 第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成ステップと、
第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換ステップと、
前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の周波数特性を示す関数の逆関数である逆フィルタの少なくとも位相を補正する逆フィルタ補正ステップと、を有し、
補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測することを特徴とする、
インパルス応答測定方法。 - 前記入力信号生成ステップは、
同じ任意波形の前記入力信号を等間隔あるいは不等間隔で繰返し発生する信号発生器、または前記入力信号と同じ信号が記録された媒体と前記入力信号を繰返し再生する再生器を測定用信号源として用い、前記測定用信号源から発生される繰返し信号を前記入力信号として前記被測定系に入力し、
前記信号変換ステップは、
前記被測定信号を受信点で受信し、前記測定用信号源と前記受信点の間で共通の同期信号を使用せずに、前記被測定信号の波形から得られた各周期の波形情報を使用して前記被測定信号を切り分け、切り分けた各周期の前記波形情報の時間のずれ量は基準とする周期と他の周期の相互相関の相関値が真の最大値となる時間のずれ量により求め、前記時間ずれの補正には各周期の波形を周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、先に変換し位相が補正された各周期の周波数ごとの振幅と位相の情報をベクトル量として加算平均し、
前記逆フィルタ補正ステップは、
前記被測定信号の自己相関あるいは繰り返す信号の隣接する信号間の相互相関により求められる前記被測定信号の信号周期と前記繰返し信号の周期とを比較して前記周波数比を求め、前記繰返し信号の周波数領域における前記逆フィルタの位相を前記周波数比に応じて補正し、
前記信号変換ステップによる加算平均の結果と、前記逆フィルタ補正ステップによる補正後の逆フィルタとの積を算出し、この算出結果を時間領域に変換して前記インパルス応答を計測することを特徴とする、
請求項1に記載のインパルス応答測定方法。 - 前記信号変換ステップは、
前記時間ずれの補正に各周期を離散的フーリエ変換(DFT)して前記周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、各周期の前記波形情報の和は、各周波数の複素ベクトル量を加算平均することを特徴とする、
請求項2に記載のインパルス応答測定方法。 - 前記逆フィルタ補正ステップは、
前記自己相関または前記相互相関による相関値の最大値を補間して求めることを特徴とする、
請求項2または3に記載のインパルス応答測定方法。 - 第1のサンプリングクロック周波数を有する同期信号を用いて、被測定系に入力する任意波形の入力信号を生成する入力信号生成手段と、
第2のサンプリングクロック周波数を有する同期信号を用いて、前記被測定系から出力される被測定信号の離散値系への変換を行う信号変換手段と、
前記第1のサンプリングクロック周波数と前記第2のサンプリングクロック周波数との周波数比に応じて前記入力信号の周波数特性を示す関数の逆関数である逆フィルタの位相を少なくとも補正する逆フィルタ補正手段と、を有し、
補正後の前記逆フィルタを用いて前記被測定系のインパルス応答を計測することを特徴とする、
インパルス応答測定装置。 - 前記入力信号生成手段は、
同じ任意波形の前記入力信号を等間隔あるいは不等間隔で繰返し発生する信号発生器、または前記入力信号と同じ信号が記録された媒体と前記入力信号を繰返し再生する再生器を測定用信号源として用い、前記測定用信号源から発生される繰返し信号を前記入力信号として前記被測定系に入力し、
前記信号変換手段は、
前記被測定信号を受信点で受信し、前記測定用信号源と前記受信点の間で共通の同期信号を使用せずに、前記被測定信号の波形から得られた各周期の波形情報を使用して前記被測定信号を切り分け、切り分けた各周期の時間のずれ量は基準とする周期と他の周期の相互相関の相関値が真の最大値となる時間のずれ量により求め、前記時間ずれの補正には各周期の波形を周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、先に変換し位相が補正された各周期の周波数ごとの振幅と位相の情報をベクトル量として加算平均し、
前記逆フィルタ補正手段は、
前記被測定信号の自己相関あるいは繰り返す信号の隣接する信号間の相互相関により求められる前記被測定信号の信号周期と前記繰返し信号の周期とを比較して前記周波数比を求め、前記繰返し信号の周波数領域における前記逆フィルタの位相を前記周波数比に応じて補正し、
前記信号変換ステップによる加算平均の結果と、前記逆フィルタ補正ステップによる補正後の逆フィルタとの積を算出し、この算出結果を時間領域に変換して前記インパルス応答を計測することを特徴とする、
請求項5に記載のインパルス応答測定装置。 - 前記信号変換手段は、
前記時間ずれの補正に各周期を離散的フーリエ変換(DFT)して前記周波数領域の振幅と位相の情報に変換した後、前記時間ずれに相当する周波数ごとの位相ずれ情報に基づいて位相を補正し、前記波形情報の各周期の和は、各周波数の複素ベクトル量を加算平均することを特徴とする、
請求項6に記載のインパルス応答測定装置。 - 前記逆フィルタ補正手段は、
前記自己相関または前記相互相関による相関値の最大値を補間して求めることを特徴とする、
請求項6または7に記載のインパルス応答測定装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/384,381 US20120166123A1 (en) | 2009-07-17 | 2010-07-07 | Impulse response measuring method and impulse response measuring device |
EP10799765.2A EP2455768A4 (en) | 2009-07-17 | 2010-07-07 | Impulse response measuring method and impulse response measuring device |
CN201080032402.XA CN102472784B (zh) | 2009-07-17 | 2010-07-07 | 脉冲响应测定方法和脉冲响应测定装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009-168427 | 2009-07-17 | ||
JP2009168427A JP5540224B2 (ja) | 2009-07-17 | 2009-07-17 | インパルス応答測定方法およびインパルス応答測定装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011007706A1 true WO2011007706A1 (ja) | 2011-01-20 |
Family
ID=43449316
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2010/061567 WO2011007706A1 (ja) | 2009-07-17 | 2010-07-07 | インパルス応答測定方法およびインパルス応答測定装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20120166123A1 (ja) |
EP (1) | EP2455768A4 (ja) |
JP (1) | JP5540224B2 (ja) |
CN (1) | CN102472784B (ja) |
WO (1) | WO2011007706A1 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015182608A1 (ja) * | 2014-05-27 | 2015-12-03 | 独立行政法人石油天然ガス・金属鉱物資源機構 | 振動検出システム、信号処理装置及び信号処理方法 |
JP2018133682A (ja) * | 2017-02-15 | 2018-08-23 | 株式会社Jvcケンウッド | フィルタ生成装置、及びフィルタ生成方法 |
RU2719325C2 (ru) * | 2015-09-08 | 2020-04-17 | Сэн-Гобэн Изовер | Способ и система получения по меньшей мере одного акустического параметра окружающей среды |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6436616B2 (ja) * | 2013-06-12 | 2018-12-12 | キヤノン株式会社 | 計測装置、計測方法、および処理装置 |
US9596553B2 (en) * | 2013-07-18 | 2017-03-14 | Harman International Industries, Inc. | Apparatus and method for performing an audio measurement sweep |
EP3627494B1 (en) | 2017-05-17 | 2021-06-23 | Panasonic Intellectual Property Management Co., Ltd. | Playback system, control device, control method, and program |
CN107144343A (zh) * | 2017-07-10 | 2017-09-08 | 薛天 | 低频振动位移传感器组网方法、系统及装置 |
JP6935370B2 (ja) * | 2018-07-24 | 2021-09-15 | アンリツ株式会社 | 信号発生装置および該装置を用いた周波数特性表示方法 |
WO2021152918A1 (ja) * | 2020-01-31 | 2021-08-05 | 住友電気工業株式会社 | 検知システム、検知装置および検知方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06265400A (ja) | 1993-03-11 | 1994-09-20 | Sony Corp | インパルス応答測定装置 |
JP2725838B2 (ja) | 1989-06-05 | 1998-03-11 | 株式会社小野測器 | インパルス応答の測定方法 |
JP2867769B2 (ja) | 1991-10-24 | 1999-03-10 | ヤマハ株式会社 | 音響測定方法およびその装置 |
JP2004193782A (ja) * | 2002-12-09 | 2004-07-08 | Toa Corp | スピーカとマイクロホン間の音波伝搬時間測定方法およびその装置 |
JP3718642B2 (ja) | 2001-06-12 | 2005-11-24 | エタニ電機株式会社 | 音響機器、音響空間、電気信号伝送線路等の伝達特性測定方法 |
WO2006011356A1 (ja) | 2004-07-29 | 2006-02-02 | Wakayama University | インパルス応答測定方法及び装置 |
JP3766975B1 (ja) | 2004-08-25 | 2006-04-19 | 辰哉 加藤 | パラメトリックな時間引き延ばしパルス生成装置 |
JP2007156393A (ja) | 2005-11-09 | 2007-06-21 | Pioneer Electronic Corp | インパルス応答検出装置及びインパルス応答検出プログラム |
JP2007232492A (ja) | 2006-02-28 | 2007-09-13 | Kawai Musical Instr Mfg Co Ltd | 伝達特性測定方法および装置 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5404379A (en) * | 1991-01-28 | 1995-04-04 | Industrial Technology Research Institute | Timing recovery method and system |
US6307868B1 (en) * | 1995-08-25 | 2001-10-23 | Terayon Communication Systems, Inc. | Apparatus and method for SCDMA digital data transmission using orthogonal codes and a head end modem with no tracking loops |
US7046694B2 (en) * | 1996-06-19 | 2006-05-16 | Digital Radio Express, Inc. | In-band on-channel digital broadcasting method and system |
US6005901A (en) * | 1997-02-27 | 1999-12-21 | Advanced Micro Devices | Arrangement for asynchronous decimation using a frequency ratio estimator and method thereof |
AU5009399A (en) * | 1998-09-24 | 2000-05-04 | Sony Corporation | Impulse response collecting method, sound effect adding apparatus, and recording medium |
US6507606B2 (en) * | 2000-03-29 | 2003-01-14 | Symmetrican, Inc. | Asymmetric digital subscriber line methods suitable for long subscriber loops |
US20030112861A1 (en) * | 2001-12-18 | 2003-06-19 | Erdogan Alper Tunga | Method and system for adaptively training time domain equalizers |
US20080144839A1 (en) * | 2005-02-28 | 2008-06-19 | Pioneer Corporation | Characteristics Measurement Device and Characteristics Measurement Program |
JP4285469B2 (ja) * | 2005-10-18 | 2009-06-24 | ソニー株式会社 | 計測装置、計測方法、音声信号処理装置 |
JP4466658B2 (ja) * | 2007-02-05 | 2010-05-26 | ソニー株式会社 | 信号処理装置、信号処理方法、プログラム |
-
2009
- 2009-07-17 JP JP2009168427A patent/JP5540224B2/ja active Active
-
2010
- 2010-07-07 EP EP10799765.2A patent/EP2455768A4/en not_active Withdrawn
- 2010-07-07 US US13/384,381 patent/US20120166123A1/en not_active Abandoned
- 2010-07-07 WO PCT/JP2010/061567 patent/WO2011007706A1/ja active Application Filing
- 2010-07-07 CN CN201080032402.XA patent/CN102472784B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2725838B2 (ja) | 1989-06-05 | 1998-03-11 | 株式会社小野測器 | インパルス応答の測定方法 |
JP2867769B2 (ja) | 1991-10-24 | 1999-03-10 | ヤマハ株式会社 | 音響測定方法およびその装置 |
JPH06265400A (ja) | 1993-03-11 | 1994-09-20 | Sony Corp | インパルス応答測定装置 |
JP3718642B2 (ja) | 2001-06-12 | 2005-11-24 | エタニ電機株式会社 | 音響機器、音響空間、電気信号伝送線路等の伝達特性測定方法 |
JP2004193782A (ja) * | 2002-12-09 | 2004-07-08 | Toa Corp | スピーカとマイクロホン間の音波伝搬時間測定方法およびその装置 |
WO2006011356A1 (ja) | 2004-07-29 | 2006-02-02 | Wakayama University | インパルス応答測定方法及び装置 |
JP3766975B1 (ja) | 2004-08-25 | 2006-04-19 | 辰哉 加藤 | パラメトリックな時間引き延ばしパルス生成装置 |
JP2007156393A (ja) | 2005-11-09 | 2007-06-21 | Pioneer Electronic Corp | インパルス応答検出装置及びインパルス応答検出プログラム |
JP2007232492A (ja) | 2006-02-28 | 2007-09-13 | Kawai Musical Instr Mfg Co Ltd | 伝達特性測定方法および装置 |
Non-Patent Citations (3)
Title |
---|
E. ORAN BRIGHAM: "The Fast Fourier Transform", 1974, PRENTICE-HALL INC. |
HINO; TSUCHIYA; ENDO: "An Examination of the Synchronous Addition Method under an Asynchronous Measurement System", AES 10TH REGIONAL CONVENTION, June 2001 (2001-06-01) |
JURO OHGA; YOSHIO YAMASAKI; YUTAKA KANEDA: "Acoustic Systems and Digital Processing", THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS, 1995, pages 158 - 159 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015182608A1 (ja) * | 2014-05-27 | 2015-12-03 | 独立行政法人石油天然ガス・金属鉱物資源機構 | 振動検出システム、信号処理装置及び信号処理方法 |
JP2015224916A (ja) * | 2014-05-27 | 2015-12-14 | 独立行政法人石油天然ガス・金属鉱物資源機構 | 振動検出システム、信号処理装置及び信号処理方法 |
US10281601B2 (en) | 2014-05-27 | 2019-05-07 | Japan Oil, Gas And Metals National Corporation | Vibration detecting system, signal processing apparatus and signal processing method |
RU2719325C2 (ru) * | 2015-09-08 | 2020-04-17 | Сэн-Гобэн Изовер | Способ и система получения по меньшей мере одного акустического параметра окружающей среды |
JP2018133682A (ja) * | 2017-02-15 | 2018-08-23 | 株式会社Jvcケンウッド | フィルタ生成装置、及びフィルタ生成方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2455768A1 (en) | 2012-05-23 |
JP2011022055A (ja) | 2011-02-03 |
EP2455768A4 (en) | 2017-10-25 |
JP5540224B2 (ja) | 2014-07-02 |
CN102472784B (zh) | 2014-07-16 |
US20120166123A1 (en) | 2012-06-28 |
CN102472784A (zh) | 2012-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5540224B2 (ja) | インパルス応答測定方法およびインパルス応答測定装置 | |
JP4552016B2 (ja) | インパルス応答測定方法及び装置 | |
JP6508539B2 (ja) | 音場収音装置および方法、音場再生装置および方法、並びにプログラム | |
WO2006126672A1 (ja) | アナログデジタル変換装置、プログラム、及び記録媒体 | |
JP2011022055A5 (ja) | ||
JP2020036271A (ja) | 伝達関数生成装置、伝達関数生成方法、およびプログラム | |
JP2002365320A (ja) | 音響用ならびに電気電子伝送路用伝達特性の測定方法 | |
JP2019052938A (ja) | 光反射測定装置及びその方法 | |
JP5035815B2 (ja) | 周波数測定装置 | |
Hafez et al. | Automatic detection of geomagnetic sudden commencement via time–frequency clusters | |
JP2007017415A (ja) | インパルス応答の時間差測定方法 | |
JP3139803B2 (ja) | インパルス応答測定装置 | |
JP2010171785A (ja) | 頭部伝達関数補間用係数算出装置、音像定位装置、頭部伝達関数補間用係数算出方法、及びプログラム | |
JP4247168B2 (ja) | 音場特性測定システム | |
CN110231515B (zh) | 一种绝缘寿命试验介质损耗因数的实时在线测量方法 | |
JP3795446B2 (ja) | Mtie試験信号発生装置 | |
JP5802573B2 (ja) | 計測器およびその制御方法 | |
Gamper | Clock drift estimation and compensation for asynchronous impulse response measurements | |
JP2008224306A (ja) | スペクトラム解析装置 | |
TWI792779B (zh) | 用於同步聲音播放及記錄時脈的系統及其使用之方法 | |
Mori et al. | Robustness of pure white pseudonoise signal to temporal fluctuation in impulse response measurement | |
JP2007156393A (ja) | インパルス応答検出装置及びインパルス応答検出プログラム | |
CN115792446A (zh) | 一种单相输电线路可听噪声声源的测试方法及装置 | |
JP3455778B2 (ja) | 信号合成装置 | |
JP4299953B2 (ja) | 振動試験制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201080032402.X Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 10799765 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2010799765 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13384381 Country of ref document: US |