JP2003333697A - 信号補間装置及びその補間方法 - Google Patents
信号補間装置及びその補間方法Info
- Publication number
- JP2003333697A JP2003333697A JP2002136343A JP2002136343A JP2003333697A JP 2003333697 A JP2003333697 A JP 2003333697A JP 2002136343 A JP2002136343 A JP 2002136343A JP 2002136343 A JP2002136343 A JP 2002136343A JP 2003333697 A JP2003333697 A JP 2003333697A
- Authority
- JP
- Japan
- Prior art keywords
- time difference
- signal
- arrival
- signals
- impulse response
- 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.)
- Pending
Links
Landscapes
- Stereophonic System (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
(57)【要約】
【課題】インパルス応答の推定処理において、到来信号
の時間差補正を行ったあと補間処理を行うことにより、
インパルス応答の推定をより正確に行うことができ、高
品質の音響環境をシミュレーションできる信号補間装置
及びその補間方法を提供する。 【解決手段】 異なる音源からの到来信号Ｓ１とＳ２に
応じて、時間差測定回路１０はこれらの到来信号の到来
時間差τを測定し、当該測定結果に基づき、時間差補正
回路２０は到来信号Ｓ１とＳ２が時間軸上揃うように補
正を行い、補間処理回路３０は上記補正結果を用いて補
間処理を行い、さらに時間差調整回路４０は補間結果に
対して測定された時間差τに基づき時間差調整を行い、
仮想音源からの到来信号Ｓ０を推定するので、仮想音源
に対応するインパルス応答の推定を正確に実現できる。
の時間差補正を行ったあと補間処理を行うことにより、
インパルス応答の推定をより正確に行うことができ、高
品質の音響環境をシミュレーションできる信号補間装置
及びその補間方法を提供する。 【解決手段】 異なる音源からの到来信号Ｓ１とＳ２に
応じて、時間差測定回路１０はこれらの到来信号の到来
時間差τを測定し、当該測定結果に基づき、時間差補正
回路２０は到来信号Ｓ１とＳ２が時間軸上揃うように補
正を行い、補間処理回路３０は上記補正結果を用いて補
間処理を行い、さらに時間差調整回路４０は補間結果に
対して測定された時間差τに基づき時間差調整を行い、
仮想音源からの到来信号Ｓ０を推定するので、仮想音源
に対応するインパルス応答の推定を正確に実現できる。
Description
【０００１】
【発明の属する技術分野】本発明は、例えば、所定の音
響環境において、それぞれ異なる方位角にある複数の音
源に対して視聴者の両耳のインパルス応答の時間差を測
定し、当該測定した時間差に応じて音源若しくは視聴者
の移動に伴う両耳のインパルス応答を補間処理により推
定する信号補間装置及びその補間方法に関するものであ
る。
響環境において、それぞれ異なる方位角にある複数の音
源に対して視聴者の両耳のインパルス応答の時間差を測
定し、当該測定した時間差に応じて音源若しくは視聴者
の移動に伴う両耳のインパルス応答を補間処理により推
定する信号補間装置及びその補間方法に関するものであ
る。
【０００２】
【従来の技術】音声を発生する音源を含む実際の音響環
境では、音源若しくは視聴者の移動に伴って、音源と視
聴者の距離が変化する。その結果、視聴者の両耳に到達
する音声信号の波形が変わらず、信号の到来時間が変わ
ることが多い。このような音響環境において、例えば、
所定の数の音源に対応するインパルス応答に応じて、他
の仮想音源に対応するインパルス応答を補間処理などの
方法により推定することが一般的に行われている。この
補間処理において、波形だけではなく、到来時間の補間
も重要となってくる。
境では、音源若しくは視聴者の移動に伴って、音源と視
聴者の距離が変化する。その結果、視聴者の両耳に到達
する音声信号の波形が変わらず、信号の到来時間が変わ
ることが多い。このような音響環境において、例えば、
所定の数の音源に対応するインパルス応答に応じて、他
の仮想音源に対応するインパルス応答を補間処理などの
方法により推定することが一般的に行われている。この
補間処理において、波形だけではなく、到来時間の補間
も重要となってくる。
【０００３】
【発明が解決しようとする課題】ところで、これまでの
補間処理では、複数の音源に対するインパルス応答に基
づいて、線形処理によって仮想音源のインパルス応答を
推定する。このとき、補間されるインパルス応答の到達
時間差は、補間に用いられるインパルス応答の時間軸に
ついての補間からではなく、振幅の補間から決定される
ため、正確に求められなかった。また、補間に用いられ
るインパルス応答を最小位相成分のみからなるインパル
ス応答に変換したあと、元の信号との共分散（Covarian
ce）が最大となる遅延時間を近似的な到来時間として、
その時間を用いて補間されるインパルス応答の到来時間
を求める方法も提案されていたが、この方法では近似的
な到来時間から間接にしか求められないという不利益が
ある。
補間処理では、複数の音源に対するインパルス応答に基
づいて、線形処理によって仮想音源のインパルス応答を
推定する。このとき、補間されるインパルス応答の到達
時間差は、補間に用いられるインパルス応答の時間軸に
ついての補間からではなく、振幅の補間から決定される
ため、正確に求められなかった。また、補間に用いられ
るインパルス応答を最小位相成分のみからなるインパル
ス応答に変換したあと、元の信号との共分散（Covarian
ce）が最大となる遅延時間を近似的な到来時間として、
その時間を用いて補間されるインパルス応答の到来時間
を求める方法も提案されていたが、この方法では近似的
な到来時間から間接にしか求められないという不利益が
ある。
【０００４】本発明は、かかる事情に鑑みてなされたも
のであり、その目的は、インパルス応答の推定処理にお
いて、到来信号の時間差補正を行ったあと補間処理を行
うことにより、インパルス応答の推定をより正確に行う
ことができ、仮想音源をよりリアルに表現でき、高品質
の音響環境をシミュレーションできる信号補間装置及び
その補間方法を提供することにある。
のであり、その目的は、インパルス応答の推定処理にお
いて、到来信号の時間差補正を行ったあと補間処理を行
うことにより、インパルス応答の推定をより正確に行う
ことができ、仮想音源をよりリアルに表現でき、高品質
の音響環境をシミュレーションできる信号補間装置及び
その補間方法を提供することにある。
【０００５】
【課題を解決するための手段】上記目的を達成するた
め、本発明の信号補間装置は、測定点から異なる距離を
離れた第１と第２の位置から発信された信号が、上記測
定点に到達した第１と第２の到来信号の到来時間差を測
定する時間差測定手段と、上記時間差測定手段によって
測定した時間差に応じて、上記測定点で測定された上記
第１と第２の到来信号が時間軸上揃うように補正する時
間差補正手段と、上記時間差補正手段によって補正され
た上記第１と第２の到来信号に応じて補間処理を行い、
上記第１と第２の位置と異なる第３の位置からの到来信
号を推定する補間処理手段とを有する。
め、本発明の信号補間装置は、測定点から異なる距離を
離れた第１と第２の位置から発信された信号が、上記測
定点に到達した第１と第２の到来信号の到来時間差を測
定する時間差測定手段と、上記時間差測定手段によって
測定した時間差に応じて、上記測定点で測定された上記
第１と第２の到来信号が時間軸上揃うように補正する時
間差補正手段と、上記時間差補正手段によって補正され
た上記第１と第２の到来信号に応じて補間処理を行い、
上記第１と第２の位置と異なる第３の位置からの到来信
号を推定する補間処理手段とを有する。
【０００６】また、本発明では、好適には、上記第１と
第２の位置が上記測定点に対して異なる方位角に配置さ
れ、上記第３の位置が上記第１と第２の位置の間に配置
されている。
第２の位置が上記測定点に対して異なる方位角に配置さ
れ、上記第３の位置が上記第１と第２の位置の間に配置
されている。
【０００７】また、本発明では、好適には、上記時間差
測定手段は、上記第１と第２の到来信号の相互相関関数
を演算する相関演算手段と、上記相関演算手段によって
算出した上記相互相関関数の最大値の位置に応じて、上
記第１と第２の到来信号の時間差を計算する時間差計算
手段とを有する。
測定手段は、上記第１と第２の到来信号の相互相関関数
を演算する相関演算手段と、上記相関演算手段によって
算出した上記相互相関関数の最大値の位置に応じて、上
記第１と第２の到来信号の時間差を計算する時間差計算
手段とを有する。
【０００８】また、本発明では、好適には、上記時間差
測定手段は、上記第１と第２の到来信号の一方を所定の
時間差で順次ずらしながら、他方の信号との差の二乗値
の和を計算する二乗差演算手段と、上記二乗差演算手段
によって算出した結果の最小値の位置に応じて、上記第
１と第２の到来信号の時間差を計算する時間差計算手段
とを有する。
測定手段は、上記第１と第２の到来信号の一方を所定の
時間差で順次ずらしながら、他方の信号との差の二乗値
の和を計算する二乗差演算手段と、上記二乗差演算手段
によって算出した結果の最小値の位置に応じて、上記第
１と第２の到来信号の時間差を計算する時間差計算手段
とを有する。
【０００９】また、本発明では、好適には、上記補間処
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号を用いて、線形補間処理によって
上記第３の位置からの到来信号を推定する。
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号を用いて、線形補間処理によって
上記第３の位置からの到来信号を推定する。
【００１０】また、本発明では、好適には、上記補間処
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号を周波数領域に変換する変換手段
と、上記変換手段によって周波数領域に変換された信号
の帯域を変換し、帯域変換された信号を時間領域逆変換
し、当該逆変換の結果に応じて上記第３の位置からの到
来信号を推定する推定手段とを有する。
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号を周波数領域に変換する変換手段
と、上記変換手段によって周波数領域に変換された信号
の帯域を変換し、帯域変換された信号を時間領域逆変換
し、当該逆変換の結果に応じて上記第３の位置からの到
来信号を推定する推定手段とを有する。
【００１１】また、本発明では、好適には、上記補間処
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号に応じて、所定の補間関数を用い
て上記第３の位置からの到来信号を推定する。
理手段は、上記時間差補正手段によって補正された上記
第１と第２の到来信号に応じて、所定の補間関数を用い
て上記第３の位置からの到来信号を推定する。
【００１２】また、本発明では、好適には、上記補間処
理手段によって補間処理された上記第３の位置からの到
来信号に対して、上記時間差測定手段によって測定され
た上記第１と第２の到来信号の時間差に応じて、到来時
間差の補正を行う補正手段をさらに有する。
理手段によって補間処理された上記第３の位置からの到
来信号に対して、上記時間差測定手段によって測定され
た上記第１と第２の到来信号の時間差に応じて、到来時
間差の補正を行う補正手段をさらに有する。
【００１３】また、本発明の信号補間方法は、測定点か
ら異なる距離を離れた第１と第２の位置から発信された
信号が、上記測定点に到達した第１と第２の到来信号の
到来時間差を測定するステップと、上記測定した時間差
に応じて、上記測定点で測定された上記第１と第２の到
来信号が時間軸上揃うように補正するステップと、上記
時間差が補正された上記第１と第２の到来信号に応じて
補間処理を行い、上記第１と第２の位置と異なる第３の
位置からの到来信号を推定するステップとを有する。
ら異なる距離を離れた第１と第２の位置から発信された
信号が、上記測定点に到達した第１と第２の到来信号の
到来時間差を測定するステップと、上記測定した時間差
に応じて、上記測定点で測定された上記第１と第２の到
来信号が時間軸上揃うように補正するステップと、上記
時間差が補正された上記第１と第２の到来信号に応じて
補間処理を行い、上記第１と第２の位置と異なる第３の
位置からの到来信号を推定するステップとを有する。
【００１４】また、本発明では、好適には、上記補間処
理された上記第３の位置からの到来信号に対して、上記
測定された上記第１と第２の到来信号の時間差に応じ
て、到来時間差の補正を行うステップをさらに有する。
理された上記第３の位置からの到来信号に対して、上記
測定された上記第１と第２の到来信号の時間差に応じ
て、到来時間差の補正を行うステップをさらに有する。
【００１５】
【発明の実施の形態】図１は本発明に係る信号補間装置
の一実施形態を示すブロック図である。図示のように、
本実施形態の信号補間装置１００は、時間差測定回路１
０、時間差補正回路２０、補間処理回路３０及び時間差
調整回路４０によって構成されている。
の一実施形態を示すブロック図である。図示のように、
本実施形態の信号補間装置１００は、時間差測定回路１
０、時間差補正回路２０、補間処理回路３０及び時間差
調整回路４０によって構成されている。
【００１６】本実施形態の信号補間装置１００は、例え
ば、所定の音響環境の中異なる音源Ｐ１とＰ２からの到
来信号Ｓ１とＳ２を受信し、受信したこれらの到来信号
に応じて、音源Ｐ１とＰ２以外の仮想音源Ｐ０からの到
来信号Ｓ０を補間処理によって推定する。
ば、所定の音響環境の中異なる音源Ｐ１とＰ２からの到
来信号Ｓ１とＳ２を受信し、受信したこれらの到来信号
に応じて、音源Ｐ１とＰ２以外の仮想音源Ｐ０からの到
来信号Ｓ０を補間処理によって推定する。
【００１７】時間差測定回路１０は、音源Ｐ１とＰ２か
らの到来信号Ｓ１とＳ２の到来時間差τを測定し、測定
結果を時間差補正回路２０及び時間差調整回路４０にそ
れぞれ供給する。時間差補正回路２０は、時間差測定回
路１０によって測定された到来時間差τに応じて、到来
信号Ｓ１とＳ２が時間軸上揃うように時間差を補正す
る。補間処理回路３０は、時間差補正回路２０によって
到来時間差が補正された２つの到来信号に対して、所定
の補間処理方法で補間処理を施し、音源Ｐ１及びＰ２と
異なる仮想音源Ｐ０からの到来信号Ｓ０を推定する。時
間差調整回路４０は、補間処理回路３０によって推定さ
れた到来信号Ｓ０に対して、時間差測定回路１０によっ
て測定した時間差τに基づいて時間差を補正する。時間
差補正の結果、推定信号Ｓ０’を出力する。
らの到来信号Ｓ１とＳ２の到来時間差τを測定し、測定
結果を時間差補正回路２０及び時間差調整回路４０にそ
れぞれ供給する。時間差補正回路２０は、時間差測定回
路１０によって測定された到来時間差τに応じて、到来
信号Ｓ１とＳ２が時間軸上揃うように時間差を補正す
る。補間処理回路３０は、時間差補正回路２０によって
到来時間差が補正された２つの到来信号に対して、所定
の補間処理方法で補間処理を施し、音源Ｐ１及びＰ２と
異なる仮想音源Ｐ０からの到来信号Ｓ０を推定する。時
間差調整回路４０は、補間処理回路３０によって推定さ
れた到来信号Ｓ０に対して、時間差測定回路１０によっ
て測定した時間差τに基づいて時間差を補正する。時間
差補正の結果、推定信号Ｓ０’を出力する。
【００１８】本実施形態の応用例として、例えば、所定
の音響環境における視聴者の両耳のインパルス応答の測
定に用いることができる。図２は、両耳のインパルス応
答を測定するとき、所定の音響環境において音源と視聴
者の空間配置を示す図である。図示のように、音源Ｐは
空間の所定の位置に配置されている。視聴者は、音源か
ら一定の距離を離れた位置にある。
の音響環境における視聴者の両耳のインパルス応答の測
定に用いることができる。図２は、両耳のインパルス応
答を測定するとき、所定の音響環境において音源と視聴
者の空間配置を示す図である。図示のように、音源Ｐは
空間の所定の位置に配置されている。視聴者は、音源か
ら一定の距離を離れた位置にある。
【００１９】図２に示す音響環境において、音源Ｐに対
して視聴者の両耳それぞれのインパルス応答ｈl （ｔ）
及びｈr （ｔ）を求めることができる。さらに、音源Ｐ
が空間において異なる位置に配置される場合、それぞれ
の位置に対して視聴者の両耳のインパルス応答を求める
ことができる。ここで、２次元空間における視聴者の両
耳インパルス応答について考える。この場合、視聴者と
音源は同一の平面上にある。視聴者の両耳インパルス応
答は、この平面上あるゆる方位角に対して求めることが
できる。こうして求められた視聴者の両耳インパルス応
答は、音響環境を高品質で表現する場合に利用される。
例えば、音響環境において有限の数の音源を用いなが
ら、仮想現実感を視聴者に与えたり、空間に仮想の音源
を作ったり、移動する音源を仮想的に発生するために利
用することができる。
して視聴者の両耳それぞれのインパルス応答ｈl （ｔ）
及びｈr （ｔ）を求めることができる。さらに、音源Ｐ
が空間において異なる位置に配置される場合、それぞれ
の位置に対して視聴者の両耳のインパルス応答を求める
ことができる。ここで、２次元空間における視聴者の両
耳インパルス応答について考える。この場合、視聴者と
音源は同一の平面上にある。視聴者の両耳インパルス応
答は、この平面上あるゆる方位角に対して求めることが
できる。こうして求められた視聴者の両耳インパルス応
答は、音響環境を高品質で表現する場合に利用される。
例えば、音響環境において有限の数の音源を用いなが
ら、仮想現実感を視聴者に与えたり、空間に仮想の音源
を作ったり、移動する音源を仮想的に発生するために利
用することができる。
【００２０】図３は、平面上において、視聴者の両耳イ
ンパルス応答を求めるときの一例を示している。図示の
ように、この場合、視聴者の頭部中心軸を中心とした円
弧上に音源を配置し、音源によって所定の音声信号を発
生しながら、それぞれの位置に対して視聴者の両耳イン
パルス応答が測定される。
ンパルス応答を求めるときの一例を示している。図示の
ように、この場合、視聴者の頭部中心軸を中心とした円
弧上に音源を配置し、音源によって所定の音声信号を発
生しながら、それぞれの位置に対して視聴者の両耳イン
パルス応答が測定される。
【００２１】図３において、例えば、視聴者の正面の方
位角を０°として、時計廻りに方位角が増加する方向と
すると、視聴者の右側は方位角９０°にあたり、また、
視聴者の後ろ方向は、方位角１８０°にあたる。図３に
示す平面において、音源Ｐの方位角を変えながら、各方
位角に対して視聴者の両耳インパルス応答を測定するこ
とにより、平面上における全方位の両耳インパルス応答
が求められる。
位角を０°として、時計廻りに方位角が増加する方向と
すると、視聴者の右側は方位角９０°にあたり、また、
視聴者の後ろ方向は、方位角１８０°にあたる。図３に
示す平面において、音源Ｐの方位角を変えながら、各方
位角に対して視聴者の両耳インパルス応答を測定するこ
とにより、平面上における全方位の両耳インパルス応答
が求められる。
【００２２】しかし、実際の測定においては、方位角を
細かい刻みで変えながら、あらゆる方位角において視聴
者の両耳インパルス応答を逐次測定することは、場合に
よって大変困難な作業である。このため、所定の間隔で
方位角度を変えて音源を設置し、各位置に対して視聴者
の両耳インパルス応答を測定し、そして、測定した方位
角以外の角度における両耳インパルス応答を、補間処
理、例えば内挿法によって推定することできる。これに
よって、少ない測定値からより多くのインパルス応答を
求めることができる。
細かい刻みで変えながら、あらゆる方位角において視聴
者の両耳インパルス応答を逐次測定することは、場合に
よって大変困難な作業である。このため、所定の間隔で
方位角度を変えて音源を設置し、各位置に対して視聴者
の両耳インパルス応答を測定し、そして、測定した方位
角以外の角度における両耳インパルス応答を、補間処
理、例えば内挿法によって推定することできる。これに
よって、少ない測定値からより多くのインパルス応答を
求めることができる。
【００２３】例えば、図３において、方位角０°と１５
°における両耳インパルス応答をそれぞれ測定し、例え
ば、方位角０°においてインパルス応答ｈi （ｔ）が得
られ、方位角１５°においてインパルス応答ｈK （ｔ）
が得られるとすると、方位角０°と１５°の間の任意の
角度ａj におけるインパルス応答ｈj （ｔ）が補間処理
によって求められる。
°における両耳インパルス応答をそれぞれ測定し、例え
ば、方位角０°においてインパルス応答ｈi （ｔ）が得
られ、方位角１５°においてインパルス応答ｈK （ｔ）
が得られるとすると、方位角０°と１５°の間の任意の
角度ａj におけるインパルス応答ｈj （ｔ）が補間処理
によって求められる。
【００２４】本発明の信号補間装置１００では、線形補
間、離散フーリエ変換（ＤＦＴ）による補間、及びスプ
ライン補間の３つの方法が取り入れられている。そし
て、補間処理の前に、それぞれの方位角におけるインパ
ルス応答の到来時間差を測定し、測定結果に基づきイン
パルス応答の時間差を補正し、時間差の補正を施したイ
ンパルス応答を用いて補間処理を行うことにより、より
正確な補間処理を実現できる。以下、到来時間差の測定
及び到来時間差の補正について説明し、そして、時間差
補正後のインパルス応答を用いた補間処理について説明
する。
間、離散フーリエ変換（ＤＦＴ）による補間、及びスプ
ライン補間の３つの方法が取り入れられている。そし
て、補間処理の前に、それぞれの方位角におけるインパ
ルス応答の到来時間差を測定し、測定結果に基づきイン
パルス応答の時間差を補正し、時間差の補正を施したイ
ンパルス応答を用いて補間処理を行うことにより、より
正確な補間処理を実現できる。以下、到来時間差の測定
及び到来時間差の補正について説明し、そして、時間差
補正後のインパルス応答を用いた補間処理について説明
する。
【００２５】到来時間差の測定方法
図３に示すように、音源Ｐは、視聴者の頭部を中心とし
た円弧上に配置されている。このため、方位角の異なる
位置にある音源に対するインパルス応答は、音源と視聴
者の耳との間の音声の伝搬経路の距離差によって、視聴
者の耳に到達する到来時間も異なる。本発明の信号補正
装置において、図１に示す時間差測定回路１０におい
て、この到来時間差τを測定し、そして、時間差補正回
路２０において、測定した到来時間差τに基づき異なる
方位角からの到来信号の到来時間差を補正し、補正結果
を補間処理回路３０に供給する。
た円弧上に配置されている。このため、方位角の異なる
位置にある音源に対するインパルス応答は、音源と視聴
者の耳との間の音声の伝搬経路の距離差によって、視聴
者の耳に到達する到来時間も異なる。本発明の信号補正
装置において、図１に示す時間差測定回路１０におい
て、この到来時間差τを測定し、そして、時間差補正回
路２０において、測定した到来時間差τに基づき異なる
方位角からの到来信号の到来時間差を補正し、補正結果
を補間処理回路３０に供給する。
【００２６】時間差測定回路１０は、相互相関演算法及
び二乗差演算法を用いて到来時間差τを測定する。以
下、それぞれの時間差測定方法について説明する。
び二乗差演算法を用いて到来時間差τを測定する。以
下、それぞれの時間差測定方法について説明する。
【００２７】相互相関演算法
相互相関演算法は、２つの到来信号の相互相関関数を演
算し、そして、相互相関関数の最大値の出現位置に応じ
て到来時間差τを求める。図４は、時間差測定回路の一
構成例を示している。図示のように、本例の時間差測定
回路１０ａは、相互相関関数演算部１１と、時間差計算
部１２とによって構成されている。
算し、そして、相互相関関数の最大値の出現位置に応じ
て到来時間差τを求める。図４は、時間差測定回路の一
構成例を示している。図示のように、本例の時間差測定
回路１０ａは、相互相関関数演算部１１と、時間差計算
部１２とによって構成されている。
【００２８】相互相関関数演算部１１は、２つの方位角
にある音源からの到来信号の相互相関関数を求める。時
間差計算部１２は、相互相関関数演算部１１によって求
められた相互相関関数の最大値を検出し、当該最大値の
位置に応じて２つの信号の到来時間差τを求める。
にある音源からの到来信号の相互相関関数を求める。時
間差計算部１２は、相互相関関数演算部１１によって求
められた相互相関関数の最大値を検出し、当該最大値の
位置に応じて２つの信号の到来時間差τを求める。
【００２９】ここで、例えば、方位角方位角ａi におけ
るインパルス応答をｈi （ｔ）とし、方位角ｋにおける
インパルス応答をｈk （ｔ）として、さらに、これらの
インパルス応答の時間サンプルをそれぞれｈi （ｎ）、
ｈk （ｎ）とする。インパルス応答ｈi （ｎ）、ｈk
（ｎ）との相互相関関数Ｒik（τ）は、次式によって計
算される。
るインパルス応答をｈi （ｔ）とし、方位角ｋにおける
インパルス応答をｈk （ｔ）として、さらに、これらの
インパルス応答の時間サンプルをそれぞれｈi （ｎ）、
ｈk （ｎ）とする。インパルス応答ｈi （ｎ）、ｈk
（ｎ）との相互相関関数Ｒik（τ）は、次式によって計
算される。
【００３０】
【数１】
【００３１】図５は、インパルス応答の波形例を示して
いる。図示のように、この波形例は、方位角９０°と１
０５°における音源からのインパルス応答の波形を示し
ている。図６は、これらのインパルス応答に基づいて算
出される相互相関関数Ｒik（τ）の波形を示している。
いる。図示のように、この波形例は、方位角９０°と１
０５°における音源からのインパルス応答の波形を示し
ている。図６は、これらのインパルス応答に基づいて算
出される相互相関関数Ｒik（τ）の波形を示している。
【００３２】時間差計算部１２は、式（１）によって求
めた相互相関関数Ｒik（τ）の最大値Ｒmax を検出し、
当該最大値Ｒmax に対応する時間τを求める。この時間
τがインパルス応答ｈi （ｎ）とｈk （ｎ）の時間差と
して出力する。
めた相互相関関数Ｒik（τ）の最大値Ｒmax を検出し、
当該最大値Ｒmax に対応する時間τを求める。この時間
τがインパルス応答ｈi （ｎ）とｈk （ｎ）の時間差と
して出力する。
【００３３】図６に示すように、相互相関関数Ｒ
ik（τ）は、ほぼτ＝１１３において最大値Ｒmax が検
出される。これに応じて、時間差計算部１２は、インパ
ルス応答ｈ i （ｎ）とｈk （ｎ）の時間差を１１３と推
定する。
ik（τ）は、ほぼτ＝１１３において最大値Ｒmax が検
出される。これに応じて、時間差計算部１２は、インパ
ルス応答ｈ i （ｎ）とｈk （ｎ）の時間差を１１３と推
定する。
【００３４】二乗差演算法
次に、二乗差演算法によるインパルス応答の時間差の求
め方について説明する。二乗差演算法は、２つのインパ
ルス応答をずらしながら、それらの差の二乗和を演算
し、当該演算結果に応じて２つのインパルス応答の時間
差を計算する。
め方について説明する。二乗差演算法は、２つのインパ
ルス応答をずらしながら、それらの差の二乗和を演算
し、当該演算結果に応じて２つのインパルス応答の時間
差を計算する。
【００３５】図７は、二乗差演算法を用いる時間差測定
回路の一構成例を示している。図示のように、本例の時
間差測定回路１０ｂは、二乗差演算部１５と、時間差計
算部１６とによって構成されている。
回路の一構成例を示している。図示のように、本例の時
間差測定回路１０ｂは、二乗差演算部１５と、時間差計
算部１６とによって構成されている。
【００３６】二乗差演算部１５は、２つの方位角にある
音源からの到来信号の時間サンプルｈi （ｎ）とｈk
（ｎ）を用いて、これらの信号をずらしながら、二乗差
関数を演算する。即ち、二乗差演算部１５では、次式に
基づき２つのインパルス応答ｈ i （ｎ）とｈk （ｎ）を
用いて二乗差関数Ｄik（τ）を計算する。
音源からの到来信号の時間サンプルｈi （ｎ）とｈk
（ｎ）を用いて、これらの信号をずらしながら、二乗差
関数を演算する。即ち、二乗差演算部１５では、次式に
基づき２つのインパルス応答ｈ i （ｎ）とｈk （ｎ）を
用いて二乗差関数Ｄik（τ）を計算する。
【００３７】
【数２】
【００３８】図８は、図５に示すインパルス応答の波形
例に対して、式（２）を用いて算出される二乗差関数Ｄ
ik（τ）の波形を示している。
例に対して、式（２）を用いて算出される二乗差関数Ｄ
ik（τ）の波形を示している。
【００３９】時間差計算部１６は、式（２）によって求
めた二乗差関数Ｄik（τ）の最小値Ｄmin を検出し、当
該最小値Ｄmin に対応する時間τを求める。そして、こ
の時間τがインパルス応答ｈi （ｎ）とｈk （ｎ）の時
間差として出力する。図８に示す二乗差関数Ｄik（τ）
は、ほぼτ＝１１３において最小値Ｄmin が検出され
る。このため、時間差計算部１６は、インパルス応答ｈ
i （ｎ）とｈk （ｎ）の時間差を１１３と推定する。
めた二乗差関数Ｄik（τ）の最小値Ｄmin を検出し、当
該最小値Ｄmin に対応する時間τを求める。そして、こ
の時間τがインパルス応答ｈi （ｎ）とｈk （ｎ）の時
間差として出力する。図８に示す二乗差関数Ｄik（τ）
は、ほぼτ＝１１３において最小値Ｄmin が検出され
る。このため、時間差計算部１６は、インパルス応答ｈ
i （ｎ）とｈk （ｎ）の時間差を１１３と推定する。
【００４０】上述したように、時間差計算回路１０は、
相互相関演算または二乗差演算により、２つのインパル
ス応答の時間差τを求めることができる。なお、図５に
示すインパルス応答の波形に対して、図６に示す相互相
関関数Ｒik（τ）の最大値Ｒ max に基づき、時間差τを
求めることができ、また、図８に示す二乗差関数Ｄ
ik（τ）の最小値Ｄmin に基づいても、時間差τを求め
ることができる。即ち、図５に示すインパルス応答の波
形例では、何れかの時間差計算方法でも正しく時間差τ
を計算することができる。
相互相関演算または二乗差演算により、２つのインパル
ス応答の時間差τを求めることができる。なお、図５に
示すインパルス応答の波形に対して、図６に示す相互相
関関数Ｒik（τ）の最大値Ｒ max に基づき、時間差τを
求めることができ、また、図８に示す二乗差関数Ｄ
ik（τ）の最小値Ｄmin に基づいても、時間差τを求め
ることができる。即ち、図５に示すインパルス応答の波
形例では、何れかの時間差計算方法でも正しく時間差τ
を計算することができる。
【００４１】相互相関関数の演算処理は積和演算のた
め、インパルス応答の波形によって、相互相関関数Ｒik
（τ）の最大値Ｒmax が明確に現れにくい場合がある。
この場合、二乗差演算法では、２つの信号の差分の二乗
を計算するので、２つの信号のもつ電力の大きさに関わ
らずその下限が“０”に近づく。また、必要に応じて、
計算された二乗差関数Ｄik（τ）の対数値を求めること
でその最小値Ｄmin をさらに正確に求められる。
め、インパルス応答の波形によって、相互相関関数Ｒik
（τ）の最大値Ｒmax が明確に現れにくい場合がある。
この場合、二乗差演算法では、２つの信号の差分の二乗
を計算するので、２つの信号のもつ電力の大きさに関わ
らずその下限が“０”に近づく。また、必要に応じて、
計算された二乗差関数Ｄik（τ）の対数値を求めること
でその最小値Ｄmin をさらに正確に求められる。
【００４２】上述したように、時間差測定回路１０にお
いて、２つのインパルス応答の相互相関関数または二乗
差関数を求め、その結果に基づき、２つのインパルス応
答の時間差τを算出できる。時間差補正回路２０は、上
述したように算出した時間差τに基づき、２つのインパ
ルス応答の時間差を補正する。
いて、２つのインパルス応答の相互相関関数または二乗
差関数を求め、その結果に基づき、２つのインパルス応
答の時間差τを算出できる。時間差補正回路２０は、上
述したように算出した時間差τに基づき、２つのインパ
ルス応答の時間差を補正する。
【００４３】時間差補正回路２０における補正処理は、
２つのインパルス応答が時間軸上揃うように、算出した
時間差だけ２つのインパルス応答の一方を時間軸上移動
させる。これによって、２つのインパルス応答の時間差
が補正される。図９は、時間差補正回路２０において、
時間差補正処理を行う場合の波形例を示している。図９
（ａ）に示すように、例えば、インパルス応答ｈi
（ｔ）はｈ k （ｔ）より時間τだけ遅れて到着する場
合、時間差補正回路２０において、インパルス応答ｈi
（ｔ）を時間軸上前方に移動させる。時間差補正の結
果、図９（ｂ）に示すように、補正後のインパルス応答
ｈi ’（ｔ）とｈk ’（ｔ）が時間軸上揃えられる。
２つのインパルス応答が時間軸上揃うように、算出した
時間差だけ２つのインパルス応答の一方を時間軸上移動
させる。これによって、２つのインパルス応答の時間差
が補正される。図９は、時間差補正回路２０において、
時間差補正処理を行う場合の波形例を示している。図９
（ａ）に示すように、例えば、インパルス応答ｈi
（ｔ）はｈ k （ｔ）より時間τだけ遅れて到着する場
合、時間差補正回路２０において、インパルス応答ｈi
（ｔ）を時間軸上前方に移動させる。時間差補正の結
果、図９（ｂ）に示すように、補正後のインパルス応答
ｈi ’（ｔ）とｈk ’（ｔ）が時間軸上揃えられる。
【００４４】時間差補正回路２０によって時間差が補正
された２つのインパルス応答ｈi ’（ｎ）とｈk ’
（ｎ）が補間処理回路３０に出力される。補間処理回路
３０では、時間補正後のインパルス応答に基づき、例え
ば、内挿法などで補間処理を行い、インパルス応答ｈi
（ｎ）とｈk （ｎ）に対応する方位角の間の任意の方位
角に対応するインパルス応答を計算する。
された２つのインパルス応答ｈi ’（ｎ）とｈk ’
（ｎ）が補間処理回路３０に出力される。補間処理回路
３０では、時間補正後のインパルス応答に基づき、例え
ば、内挿法などで補間処理を行い、インパルス応答ｈi
（ｎ）とｈk （ｎ）に対応する方位角の間の任意の方位
角に対応するインパルス応答を計算する。
【００４５】補間処理回路３０における補間処理は、線
形補間法、ＤＦＴを用いた補間法及びスプライン法など
に基づいて行われる。以下、それぞれの補間法について
説明する。なお、以下の説明では、インパルス応答を時
間サンプルとして取り扱う。ここで、方位角ａi におけ
るインパルス応答をｈi （ｎ）とし、方位角ｋにおける
インパルス応答をｈk （ｎ）とする。補間処理回路３０
における補間処理は、これらのインパルス応答に基づ
き、方位角ａi とａk の間にある任意の方位角ａj にお
けるインパルス応答ｈj （ｎ）を求める。
形補間法、ＤＦＴを用いた補間法及びスプライン法など
に基づいて行われる。以下、それぞれの補間法について
説明する。なお、以下の説明では、インパルス応答を時
間サンプルとして取り扱う。ここで、方位角ａi におけ
るインパルス応答をｈi （ｎ）とし、方位角ｋにおける
インパルス応答をｈk （ｎ）とする。補間処理回路３０
における補間処理は、これらのインパルス応答に基づ
き、方位角ａi とａk の間にある任意の方位角ａj にお
けるインパルス応答ｈj （ｎ）を求める。
【００４６】線形補間法
線形補間法は、方位角ａi とａk においてそれぞれ測定
したインパルス応答を用いて、線形演算により方位角ａ
i とａk の間にある方位角ａj におけるインパルス応答
を計算する。
したインパルス応答を用いて、線形演算により方位角ａ
i とａk の間にある方位角ａj におけるインパルス応答
を計算する。
【００４７】具体的に、補間処理装置３０は、次式によ
り、方位角ａi とａk においてそれぞれ測定したインパ
ルス応答ｈi （ｎ）とｈk （ｎ）に基づき、方位角ａj
におけるインパルス応答ｈj （ｎ）を演算する。
り、方位角ａi とａk においてそれぞれ測定したインパ
ルス応答ｈi （ｎ）とｈk （ｎ）に基づき、方位角ａj
におけるインパルス応答ｈj （ｎ）を演算する。
【００４８】
【数３】
ｈj （ｎ）＝（ａｈk （ｎ）＋ｂｈi （ｎ））／（ａ＋ｂ）
…（３）
【００４９】式（３）において、ａは方位角ａi とａj
との間の角度を示し、ｂは方位角ａ j とａk との間の角
度を示す。
との間の角度を示し、ｂは方位角ａ j とａk との間の角
度を示す。
【００５０】ＤＦＴを用いた補間法
図１０は、ＤＦＴを用いた補間処理の方法を示してい
る。図１０において、行列Ｈの各列は、所定の角度間隔
で測定したインパルス応答の時間サンプルからなる。こ
のため、インパルス応答を補間処理によって推定するこ
とは、図１０に示す行列の列の数を補間処理によって増
やすことである。なお、行列Ｈの各列を形成するインパ
ルス応答は、上述した時間差補正処理によって、時間差
が補正されたインパルス応答である。
る。図１０において、行列Ｈの各列は、所定の角度間隔
で測定したインパルス応答の時間サンプルからなる。こ
のため、インパルス応答を補間処理によって推定するこ
とは、図１０に示す行列の列の数を補間処理によって増
やすことである。なお、行列Ｈの各列を形成するインパ
ルス応答は、上述した時間差補正処理によって、時間差
が補正されたインパルス応答である。
【００５１】ＤＦＴを用いた補間処理では、まず、行列
Ｈの各行からなる行ベクトルに対して離散フーリエ変換
を行い、行ベクトルを周波数領域に変換する。図１０に
示すように、行ベクトルｓは、各インパルス応答の同じ
時刻のサンプルから構成されている。そして、行ベクト
ルｓに対してＤＦＴ処理の結果、周波数領域の行ベクト
ルｖが得られる。
Ｈの各行からなる行ベクトルに対して離散フーリエ変換
を行い、行ベクトルを周波数領域に変換する。図１０に
示すように、行ベクトルｓは、各インパルス応答の同じ
時刻のサンプルから構成されている。そして、行ベクト
ルｓに対してＤＦＴ処理の結果、周波数領域の行ベクト
ルｖが得られる。
【００５２】次に、行ベクトルｖに対して０を挿入し、
行ベクトルの係数を増やす。図１０に示すように、行ベ
クトルｖの係数が低周波数帯域に移動され、高周波数帯
域に０が埋め込まれる。こうして得られた行ベクトル
ｖ’が逆ＤＦＴを用いて時間領域に逆変換され、行ベク
トルｓ’が得られる。この結果、行ベクトルｓ’は、処
理前の行ベクトルｓに較べて、係数の数が増えた。この
処理を行列Ｈのすべての行ベクトルに対して行われる。
この結果、列の数が増えた新しい行列Ｈ’が得られる。
例えば、処理前の行列Ｈにおいて、隣接した列ｈi
（ｎ）とｈk （ｎ）の間に、新しい列ｈj （ｎ）が生成
される。この列ｈj （ｎ）は、インパルス応答ｈi
（ｎ）とｈk （ｎ）に対応する音源の間に位置する仮想
の音源に対応するインパルス応答である。
行ベクトルの係数を増やす。図１０に示すように、行ベ
クトルｖの係数が低周波数帯域に移動され、高周波数帯
域に０が埋め込まれる。こうして得られた行ベクトル
ｖ’が逆ＤＦＴを用いて時間領域に逆変換され、行ベク
トルｓ’が得られる。この結果、行ベクトルｓ’は、処
理前の行ベクトルｓに較べて、係数の数が増えた。この
処理を行列Ｈのすべての行ベクトルに対して行われる。
この結果、列の数が増えた新しい行列Ｈ’が得られる。
例えば、処理前の行列Ｈにおいて、隣接した列ｈi
（ｎ）とｈk （ｎ）の間に、新しい列ｈj （ｎ）が生成
される。この列ｈj （ｎ）は、インパルス応答ｈi
（ｎ）とｈk （ｎ）に対応する音源の間に位置する仮想
の音源に対応するインパルス応答である。
【００５３】上述したように、インパルス応答によって
構成された行列Ｈの各行ベクトルに対して、ＤＦＴで周
波数領域に変換された行ベクトルに数値０を挿入したベ
クトルを逆ＤＦＴを用いて時間領域に逆変換した結果、
新しいインパルス応答を推定できる。
構成された行列Ｈの各行ベクトルに対して、ＤＦＴで周
波数領域に変換された行ベクトルに数値０を挿入したベ
クトルを逆ＤＦＴを用いて時間領域に逆変換した結果、
新しいインパルス応答を推定できる。
【００５４】スプライン補間法
スプライン補間法は、所定の補間関数を用いて２つ若し
くは２つ以上のインパルス応答の時間サンプルを用い
て、新しいインパルス応答の時間サンプルの値を計算す
る。スプライン補間に使用される補間関数は、例えば、
キュービック関数と呼ばれている３次関数が一般的に使
用されている。キュービック関数を用いて、例えば、隣
接する４つのインパルス応答サンプルに応じて、新しい
インパルス応答のサンプルを計算することができる。こ
れによって、測定されたインパルス応答から新しいイン
パルス応答を推定することができる。
くは２つ以上のインパルス応答の時間サンプルを用い
て、新しいインパルス応答の時間サンプルの値を計算す
る。スプライン補間に使用される補間関数は、例えば、
キュービック関数と呼ばれている３次関数が一般的に使
用されている。キュービック関数を用いて、例えば、隣
接する４つのインパルス応答サンプルに応じて、新しい
インパルス応答のサンプルを計算することができる。こ
れによって、測定されたインパルス応答から新しいイン
パルス応答を推定することができる。
【００５５】本実施形態の信号補間装置１００におい
て、補間処理回路３０において上述した線形補間法、Ｄ
ＦＴ補間法またはスプライン補間法の何れかを用いて補
間処理を行い、仮想の音源に対応するインパルス応答を
求める。そして、時間差調整回路４０において、補間処
理によって求められたインパルス応答に対して、さらに
時間補正を行うことによって、より正確に両耳インパル
ス応答を求める。
て、補間処理回路３０において上述した線形補間法、Ｄ
ＦＴ補間法またはスプライン補間法の何れかを用いて補
間処理を行い、仮想の音源に対応するインパルス応答を
求める。そして、時間差調整回路４０において、補間処
理によって求められたインパルス応答に対して、さらに
時間補正を行うことによって、より正確に両耳インパル
ス応答を求める。
【００５６】補間処理回路３０によって求められたイン
パルス応答の時間補正は、時間差測定回路１０によって
測定されたインパルス応答ｈi （ｎ）とｈk （ｎ）の到
来時間差τ及び補間処理で求められたインパルス応答ｈ
j （ｎ）の方位角に応じて行われる。
パルス応答の時間補正は、時間差測定回路１０によって
測定されたインパルス応答ｈi （ｎ）とｈk （ｎ）の到
来時間差τ及び補間処理で求められたインパルス応答ｈ
j （ｎ）の方位角に応じて行われる。
【００５７】ここで、例えば、インパルス応答ｈi
（ｎ）に対応する音源の方位角をａi 、インパルス応答
ｈk （ｎ）に対応する音源の方位角をａk とし、補間処
理で求められたインパルス応答ｈj （ｎ）に対応する音
源の方位角をａj とする。そして、方位角ａi とａj と
の間の角度をａとし、方位角ａj とａk との間の角度を
ｂとする。インパルス応答ｈj （ｎ）の到来時間差τj
は、次式のように線形演算によって求められる。
（ｎ）に対応する音源の方位角をａi 、インパルス応答
ｈk （ｎ）に対応する音源の方位角をａk とし、補間処
理で求められたインパルス応答ｈj （ｎ）に対応する音
源の方位角をａj とする。そして、方位角ａi とａj と
の間の角度をａとし、方位角ａj とａk との間の角度を
ｂとする。インパルス応答ｈj （ｎ）の到来時間差τj
は、次式のように線形演算によって求められる。
【００５８】
【数４】
τj ＝ｂτ／（ａ＋ｂ） …（４）
【００５９】このように補間処理によって求められたイ
ンパルス応答ｈj （ｎ）に対して、インパルス応答ｈi
（ｎ）とｈk （ｎ）の到来時間差τ及び補間処理で求め
られたインパルス応答ｈj （ｎ）の方位角に応じて、式
（４）で時間差の補正を行うことができる。これによっ
て、より正確にインパルス応答を推定することができ
る。
ンパルス応答ｈj （ｎ）に対して、インパルス応答ｈi
（ｎ）とｈk （ｎ）の到来時間差τ及び補間処理で求め
られたインパルス応答ｈj （ｎ）の方位角に応じて、式
（４）で時間差の補正を行うことができる。これによっ
て、より正確にインパルス応答を推定することができ
る。
【００６０】
【発明の効果】以上説明したように、本発明の信号補間
装置及びその補間方法によれば、所定の音響環境におい
て異なる空間位置にある音源に対応するインパルス応答
に基づいて、補間処理により他の空間位置にある仮想音
源に対応するインパルス応答を求めることができる。ま
た、本発明において、異なる空間位置からの到来信号の
到来時間差を測定し、当該測定結果に応じて到来時間差
を補正し、補正結果に基づいて補間処理を行う。そし
て、補間処理の結果に対して上記測定した時間差に基づ
きさらに時間軸上の補正処理を行うことで、より正確に
インパルス応答を推定することができる利点がある。こ
のため、本発明によって推定された仮想音源のインパル
ス応答を用いることによって、仮想音源をよりリアルに
表現でき、高品質の音響環境をシミュレーションでき
る。
装置及びその補間方法によれば、所定の音響環境におい
て異なる空間位置にある音源に対応するインパルス応答
に基づいて、補間処理により他の空間位置にある仮想音
源に対応するインパルス応答を求めることができる。ま
た、本発明において、異なる空間位置からの到来信号の
到来時間差を測定し、当該測定結果に応じて到来時間差
を補正し、補正結果に基づいて補間処理を行う。そし
て、補間処理の結果に対して上記測定した時間差に基づ
きさらに時間軸上の補正処理を行うことで、より正確に
インパルス応答を推定することができる利点がある。こ
のため、本発明によって推定された仮想音源のインパル
ス応答を用いることによって、仮想音源をよりリアルに
表現でき、高品質の音響環境をシミュレーションでき
る。
【図１】本発明に係る信号補間装置の一実施形態を示す
ブロック図である。
ブロック図である。
【図２】本発明の信号補間装置を用いて視聴者の両耳イ
ンパルス応答を測定する場合、音源と視聴者の空間配置
を示す図である。
ンパルス応答を測定する場合、音源と視聴者の空間配置
を示す図である。
【図３】平面上において視聴者の両耳インパルス応答を
求める場合の一例を示す図である。
求める場合の一例を示す図である。
【図４】時間差測定回路の一構成例を示すブロック図で
ある。
ある。
【図５】インパルス応答の波形例を示す波形図である。
【図６】インパルス応答の相互相関関数Ｒikの波形を示
す図である。
す図である。
【図７】時間差測定回路の他の構成例を示すブロック図
である。
である。
【図８】インパルス応答の二乗差関数Ｄikの波形を示す
図である。
図である。
【図９】時間差補正処理を示す波形図である。
【図１０】ＤＦＴを用いた補間処理法を示す概念図であ
る。
る。
１００…信号補間回路、１０，１０ａ，１０ｂ…時間差
測定回路、１１…相関演算部、１２…時間差計算部、１
５…二乗差演算部、１６…時間差計算部、２０…時間差
補正回路、３０…補間処理回路、４０…時間差調整回
路。
測定回路、１１…相関演算部、１２…時間差計算部、１
５…二乗差演算部、１６…時間差計算部、２０…時間差
補正回路、３０…補間処理回路、４０…時間差調整回
路。
【手続補正書】
【提出日】平成１４年５月２１日（２００２．５．２
１）
１）
【手続補正１】
【補正対象書類名】図面
【補正対象項目名】図２
【補正方法】変更
【補正内容】
【図２】
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図５
【補正方法】変更
【補正内容】
【図５】
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図６
【補正方法】変更
【補正内容】
【図６】
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図８
【補正方法】変更
【補正内容】
【図８】
─────────────────────────────────────────────────────
フロントページの続き
Ｆターム(参考） 5D062 AA65 AA67 AA74 CC16
5J083 AA05 AC07 AC29 AD04 AD18
BE10 BE41 EA08
Claims (11)
- 【請求項１】測定点から異なる距離を離れた第１と第２
の位置から発信された信号が、上記測定点に到達した第
１と第２の到来信号の到来時間差を測定する時間差測定
手段と、 上記時間差測定手段によって測定した時間差に応じて、
上記測定点で測定された上記第１と第２の到来信号が時
間軸上揃うように補正する時間差補正手段と、 上記時間差補正手段によって補正された上記第１と第２
の到来信号に応じて補間処理を行い、上記第１と第２の
位置と異なる第３の位置からの到来信号を推定する補間
処理手段とを有する信号補間装置。 - 【請求項２】上記第１と第２の位置が、上記測定点に対
して異なる方位角に配置されている請求項１記載の信号
補間装置。 - 【請求項３】上記第３の位置が、上記第１と第２の位置
の間に配置されている請求項２記載の信号補間装置。 - 【請求項４】上記時間差測定手段は、上記第１と第２の
到来信号の相互相関関数を演算する相関演算手段と、 上記相関演算手段によって算出した上記相互相関関数の
最大値の位置に応じて、上記第１と第２の到来信号の時
間差を計算する時間差計算手段とを有する請求項１記載
の信号補間装置。 - 【請求項５】上記時間差測定手段は、上記第１と第２の
到来信号の一方を所定の時間差で順次ずらしながら、他
方の信号との差の二乗値の和を計算する二乗差演算手段
と、 上記二乗差演算手段によって算出した結果の最小値の位
置に応じて、上記第１と第２の到来信号の時間差を計算
する時間差計算手段とを有する請求項１記載の信号補間
装置。 - 【請求項６】上記補間処理手段は、上記時間差補正手段
によって補正された上記第１と第２の到来信号を用い
て、線形補間処理によって上記第３の位置からの到来信
号を推定する請求項１記載の信号補間装置。 - 【請求項７】上記補間処理手段は、上記時間差補正手段
によって補正された上記第１と第２の到来信号を周波数
領域に変換する変換手段と、 上記変換手段によって周波数領域に変換された信号の帯
域を変換し、帯域変換された信号を時間領域逆変換し、
当該逆変換の結果に応じて上記第３の位置からの到来信
号を推定する推定手段とを有する請求項１記載の信号補
間装置。 - 【請求項８】上記補間処理手段は、上記時間差補正手段
によって補正された上記第１と第２の到来信号に応じ
て、所定の補間関数を用いて上記第３の位置からの到来
信号を推定する請求項１記載の信号補間装置。 - 【請求項９】上記補間処理手段によって補間処理された
上記第３の位置からの到来信号に対して、上記時間差測
定手段によって測定された上記第１と第２の到来信号の
時間差に応じて、到来時間差の補正を行う補正手段をさ
らに有する請求項１記載の信号補間装置。 - 【請求項１０】測定点から異なる距離を離れた第１と第
２の位置から発信された信号が、上記測定点に到達した
第１と第２の到来信号の到来時間差を測定するステップ
と、 上記測定した時間差に応じて、上記測定点で測定された
上記第１と第２の到来信号が時間軸上揃うように補正す
るステップと、 上記時間差が補正された上記第１と第２の到来信号に応
じて補間処理を行い、上記第１と第２の位置と異なる第
３の位置からの到来信号を推定するステップとを有する
信号補間方法。 - 【請求項１１】上記補間処理された上記第３の位置から
の到来信号に対して、上記測定された上記第１と第２の
到来信号の時間差に応じて、到来時間差の補正を行うス
テップをさらに有する請求項１０記載の信号補間方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002136343A JP2003333697A (ja) | 2002-05-10 | 2002-05-10 | 信号補間装置及びその補間方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002136343A JP2003333697A (ja) | 2002-05-10 | 2002-05-10 | 信号補間装置及びその補間方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2003333697A true JP2003333697A (ja) | 2003-11-21 |
Family
ID=29698377
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002136343A Pending JP2003333697A (ja) | 2002-05-10 | 2002-05-10 | 信号補間装置及びその補間方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2003333697A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006222801A (ja) * | 2005-02-10 | 2006-08-24 | Nec Tokin Corp | 移動音像提示装置 |
| JP2010041425A (ja) * | 2008-08-05 | 2010-02-18 | Yamaha Corp | 音響再生装置およびプログラム |
| US7933421B2 (en) | 2004-05-28 | 2011-04-26 | Sony Corporation | Sound-field correcting apparatus and method therefor |
| JP2015154207A (ja) * | 2014-02-13 | 2015-08-24 | 本田技研工業株式会社 | 音響処理装置、及び音響処理方法 |
-
2002
- 2002-05-10 JP JP2002136343A patent/JP2003333697A/ja active Pending
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7933421B2 (en) | 2004-05-28 | 2011-04-26 | Sony Corporation | Sound-field correcting apparatus and method therefor |
| JP2006222801A (ja) * | 2005-02-10 | 2006-08-24 | Nec Tokin Corp | 移動音像提示装置 |
| JP2010041425A (ja) * | 2008-08-05 | 2010-02-18 | Yamaha Corp | 音響再生装置およびプログラム |
| JP2015154207A (ja) * | 2014-02-13 | 2015-08-24 | 本田技研工業株式会社 | 音響処理装置、及び音響処理方法 |
| US10139470B2 (en) | 2014-02-13 | 2018-11-27 | Honda Motor Co., Ltd. | Sound processing apparatus and sound processing method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10609504B2 (en) | Audio signal processing method and apparatus for binaural rendering using phase response characteristics | |
| Gillette et al. | A linear closed-form algorithm for source localization from time-differences of arrival | |
| US10477337B2 (en) | Audio processing device and method therefor | |
| EP0788723B1 (en) | Method and apparatus for efficient presentation of high-quality three-dimensional audio | |
| KR101591220B1 (ko) | 공간적 전력 밀도에 기초하여 마이크 위치 결정을 위한 장치 및 방법 | |
| Sachar et al. | Microphone position and gain calibration for a large-aperture microphone array | |
| Kotus et al. | Calibration of acoustic vector sensor based on MEMS microphones for DOA estimation | |
| CN109188362A (zh) | 一种麦克风阵列声源定位信号处理方法 | |
| KR20090128221A (ko) | 음원 위치 추정 방법 및 그 방법에 따른 시스템 | |
| Matsumoto et al. | A method of interpolating binaural impulse responses for moving sound images | |
| EP2716074A1 (en) | Method for self - calibrating a set of acoustic sensors, and corresponding system | |
| JP2003333697A (ja) | 信号補間装置及びその補間方法 | |
| KR101103361B1 (ko) | 입체 음향 시스템 및 입체 음향을 구현하는 방법 | |
| KR101815584B1 (ko) | 음원의 위치를 정밀 추정하기 위한 3차원 음향 인텐시티 탐촉장치의 편향오차 보상방법 | |
| RU2623831C1 (ru) | Способ пассивного определения координат движущегося источника излучения | |
| Annibale et al. | Geometric inference of the room geometry under temperature variations | |
| US20230336938A1 (en) | Efficient head-related filter generation | |
| JP2005077205A (ja) | 音源方向推定装置、信号の時間遅延推定装置及びコンピュータプログラム | |
| Keyrouz et al. | A rational HRTF interpolation approach for fast synthesis of moving sound | |
| Koyama | Boundary integral approach to sound field transform and reproduction | |
| CN115699811A (zh) | 头部相关（hr）滤波器 | |
| Mandlik et al. | Real-time sound source localization | |
| Vina et al. | A spherical cross-channel algorithm for binaural sound localization | |
| CN111157949A (zh) | 一种语音识别及声源定位方法 | |
| EP2908125B1 (en) | Ultrasound based measurement apparatus and method |