JP2870440B2 - 立体音場再生方式 - Google Patents
立体音場再生方式Info
- Publication number
- JP2870440B2 JP2870440B2 JP7025210A JP2521095A JP2870440B2 JP 2870440 B2 JP2870440 B2 JP 2870440B2 JP 7025210 A JP7025210 A JP 7025210A JP 2521095 A JP2521095 A JP 2521095A JP 2870440 B2 JP2870440 B2 JP 2870440B2
- Authority
- JP
- Japan
- Prior art keywords
- filter
- transfer function
- output
- error
- sound field
- 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.)
- Expired - Lifetime
Links
Landscapes
- Stereophonic System (AREA)
- Circuit For Audible Band Transducer (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Stereophonic Arrangements (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、２チャンネルスピーカ
により立体音場を再生する立体音場再生方式に関する。
により立体音場を再生する立体音場再生方式に関する。
【０００２】
【従来の技術】従来、この種の立体音場再生方式は、右
側スピーカから聴取者の左の耳に伝わる信号および左側
スピーカから聴取者の右の耳に伝わる信号（以下、クロ
ストークパス信号という）をキャンセルして、右側スピ
ーカから聴取者の右の耳に伝わる信号および左側スピー
カから聴取者の左の耳に伝わる信号（以下、メインパス
信号という）のみを聴取者に伝送することにより、通常
のオーディオ信号により再生される音場に比べて臨場感
のある音場を再生するものである。
側スピーカから聴取者の左の耳に伝わる信号および左側
スピーカから聴取者の右の耳に伝わる信号（以下、クロ
ストークパス信号という）をキャンセルして、右側スピ
ーカから聴取者の右の耳に伝わる信号および左側スピー
カから聴取者の左の耳に伝わる信号（以下、メインパス
信号という）のみを聴取者に伝送することにより、通常
のオーディオ信号により再生される音場に比べて臨場感
のある音場を再生するものである。
【０００３】図３は、従来の立体音場再生方式を説明す
るための図である。
るための図である。
【０００４】従来の立体音場再生方式に用いられるフィ
ルタ構成は、４つのＦＩＲフィルタ（インパルス応答が
有限個の系列で表されるデジタルフィルタ）よりなる次
のような構成となっている。左側スピーカ１０９用のオ
ーディオデジタル信号Ｌをそれぞれ入力とするＦＩＲフ
ィルタＸLL１０１，ＸLR１０２と、右側スピーカ１１０
用オーディオデジタル信号Ｒをそれぞれ入力とするＦＩ
ＲフィルタＸRL１０３，ＸRR１０４とにより構成され、
ＦＩＲフィルタＸLL１０１およびＦＩＲフィルタＸRL１
０３の各出力が、Ｄ／Ａコンバータ１０５、アンプ１０
７を順次介して左側スピーカ１０９の入力となってお
り、ＦＩＲフィルタＸLR１０２およびＦＩＲフィルタＸ
RR１０４の各出力が、Ｄ／Ａコンバータ１０６、アンプ
１０８を順次介して右側スピーカ１１０の入力となって
いる。
ルタ構成は、４つのＦＩＲフィルタ（インパルス応答が
有限個の系列で表されるデジタルフィルタ）よりなる次
のような構成となっている。左側スピーカ１０９用のオ
ーディオデジタル信号Ｌをそれぞれ入力とするＦＩＲフ
ィルタＸLL１０１，ＸLR１０２と、右側スピーカ１１０
用オーディオデジタル信号Ｒをそれぞれ入力とするＦＩ
ＲフィルタＸRL１０３，ＸRR１０４とにより構成され、
ＦＩＲフィルタＸLL１０１およびＦＩＲフィルタＸRL１
０３の各出力が、Ｄ／Ａコンバータ１０５、アンプ１０
７を順次介して左側スピーカ１０９の入力となってお
り、ＦＩＲフィルタＸLR１０２およびＦＩＲフィルタＸ
RR１０４の各出力が、Ｄ／Ａコンバータ１０６、アンプ
１０８を順次介して右側スピーカ１１０の入力となって
いる。
【０００５】この立体音場再生方式では、オーディオ信
号Ｌ，Ｒがそれぞれ入力されると、各ＦＩＲフィルタＸ
LL１０１，ＸLR１０２，ＸRL１０３，ＸRR１０４によっ
て、左右のスピーカ１０９，１１０から信号が出力され
る際のクロストークパス信号がキャンセルされ、聴取者
にメインパス信号のみが伝送される。このときの各ＦＩ
ＲフィルタＸLL１０１，ＸLR１０２，ＸRL１０３，ＸRR
１０４のフィルタ係数は、以下のようにして決定されて
いる。
号Ｌ，Ｒがそれぞれ入力されると、各ＦＩＲフィルタＸ
LL１０１，ＸLR１０２，ＸRL１０３，ＸRR１０４によっ
て、左右のスピーカ１０９，１１０から信号が出力され
る際のクロストークパス信号がキャンセルされ、聴取者
にメインパス信号のみが伝送される。このときの各ＦＩ
ＲフィルタＸLL１０１，ＸLR１０２，ＸRL１０３，ＸRR
１０４のフィルタ係数は、以下のようにして決定されて
いる。
【０００６】各ＦＩＲフィルタＸLL１０１，ＸLR１０
２，ＸRL１０３，ＸRR１０４のフィルタ係数は、左右の
スピーカ１０９，１１０からそれぞれ聴取者の左右の耳
に伝わる伝達関数を用いて表すことができる。すなわ
ち、右側スピーカ１０９および左側スピーカ１１０から
左耳に伝わる伝達関数をそれぞれＨRL、ＨLLとし、右側
スピーカ１０９および左側スピーカ１１０から右耳に伝
わる伝達関数をそれぞれＨ RR、ＨLRとすると、各ＦＩＲ
フィルタＸLL１０１，ＸLR１０２，ＸRL１０３，Ｘ RR１
０４のフィルタ係数は以下のような式で表される。
２，ＸRL１０３，ＸRR１０４のフィルタ係数は、左右の
スピーカ１０９，１１０からそれぞれ聴取者の左右の耳
に伝わる伝達関数を用いて表すことができる。すなわ
ち、右側スピーカ１０９および左側スピーカ１１０から
左耳に伝わる伝達関数をそれぞれＨRL、ＨLLとし、右側
スピーカ１０９および左側スピーカ１１０から右耳に伝
わる伝達関数をそれぞれＨ RR、ＨLRとすると、各ＦＩＲ
フィルタＸLL１０１，ＸLR１０２，ＸRL１０３，Ｘ RR１
０４のフィルタ係数は以下のような式で表される。
【０００７】 ＸLL＝ＨRR÷（ＨLL×ＨRR−ＨLR×ＨRL） （１） ＸLR＝−ＨLR÷（ＨLL×ＨRR−ＨLR×ＨRL） （２） ＸRR＝ＨLL÷（ＨLL×ＨRR−ＨLR×ＨRL） （３） ＸRL＝−ＨRL÷（ＨLL×ＨRR−ＨLR×ＨRL） （４） 各フィルタ係数ＸLL，ＸLR，ＸRL，ＸRRを決定する各伝
達関数ＨLL，ＨRR，Ｈ LR，ＨRLは、無響室等にて、所望
波と相関の強い信号（参照信号）を系のなかで発生し、
所定の聴取位置に設置されたダミーヘッドマイクから入
力される信号（音圧に応じた信号）と参照信号の差の２
乗平均誤差を最小化することにより求められる。この無
響室等にて求められた各伝達関数を基に、ＦＦＴ（Fast
FourieTransfor の略）による畳み込み演算処理を行う
ことにより上記式（１）〜（４）における各ＦＩＲフィ
ルタＸLL１０１，ＸLR１０２，ＸRL１０３，ＸRR１０４
のフィルタ係数が決定される。
達関数ＨLL，ＨRR，Ｈ LR，ＨRLは、無響室等にて、所望
波と相関の強い信号（参照信号）を系のなかで発生し、
所定の聴取位置に設置されたダミーヘッドマイクから入
力される信号（音圧に応じた信号）と参照信号の差の２
乗平均誤差を最小化することにより求められる。この無
響室等にて求められた各伝達関数を基に、ＦＦＴ（Fast
FourieTransfor の略）による畳み込み演算処理を行う
ことにより上記式（１）〜（４）における各ＦＩＲフィ
ルタＸLL１０１，ＸLR１０２，ＸRL１０３，ＸRR１０４
のフィルタ係数が決定される。
【０００８】
【発明が解決しようとする課題】しかしながら、上述の
従来の立体音場再生方式には、以下のような問題があ
る。
従来の立体音場再生方式には、以下のような問題があ
る。
【０００９】（１）伝達関数を求める際、フィルタの収
束性が十分でないためにダミーヘッドマイクを用いた音
圧の測定に遅延誤差が生じる。そのため、従来例により
求められた伝達関数は、理想伝達関数に対してその遅延
誤差分の誤差を含んだものとなる。このように遅延誤差
を含んだ伝達関数を用いてＦＩＲフィルタのフィルタ係
数の決定が行われる従来の方式では、正確にクロストー
クパス信号をキャンセルすることができず、正確にメイ
ンパス信号を聴取者に伝達できないという問題点があ
る。
束性が十分でないためにダミーヘッドマイクを用いた音
圧の測定に遅延誤差が生じる。そのため、従来例により
求められた伝達関数は、理想伝達関数に対してその遅延
誤差分の誤差を含んだものとなる。このように遅延誤差
を含んだ伝達関数を用いてＦＩＲフィルタのフィルタ係
数の決定が行われる従来の方式では、正確にクロストー
クパス信号をキャンセルすることができず、正確にメイ
ンパス信号を聴取者に伝達できないという問題点があ
る。
【００１０】（２）従来の方式では、各ＦＩＲフィルタ
のフィルタ係数は、無響室等での測定により求められた
伝達関数を基に決定される。しかし、室内音響特性は設
置環境によって異なるため、無響室等での測定により求
められた伝達関数を基にフィルタ係数が決定されたもの
では、無響室と室内音響特性が大きく異なる環境におい
ては、十分にクロストークパス信号をキャンセルするこ
とはできず、所望の臨場感を得ることができないという
問題点がある。
のフィルタ係数は、無響室等での測定により求められた
伝達関数を基に決定される。しかし、室内音響特性は設
置環境によって異なるため、無響室等での測定により求
められた伝達関数を基にフィルタ係数が決定されたもの
では、無響室と室内音響特性が大きく異なる環境におい
ては、十分にクロストークパス信号をキャンセルするこ
とはできず、所望の臨場感を得ることができないという
問題点がある。
【００１１】本発明の目的は、遅延誤差を考慮した伝達
関数を基にＦＩＲフィルタのフィルタ係数を求めること
ができる立体音場再生方式を提供することにある。
関数を基にＦＩＲフィルタのフィルタ係数を求めること
ができる立体音場再生方式を提供することにある。
【００１２】さらには、上記フィルタ特性を設置環境に
よって変化する室内音響特性に最適なものとして正確に
メインパス信号を聴取者に伝達することのできる立体音
場再生方式を提供することにある。
よって変化する室内音響特性に最適なものとして正確に
メインパス信号を聴取者に伝達することのできる立体音
場再生方式を提供することにある。
【００１３】
【課題を解決するための手段】上記目的を達成するた
め、本発明の立体音場再生方式は、立体音場を生成する
ための第１および第２のオーディオ信号と、前記第１の
オーディオ信号入力とする第１および第２のＦＩＲフィ
ルタと、前記第２のオーディオ信号を入力とする第３お
よび第４のＦＩＲフィルタと、前記第１および第３のＦ
ＩＲフィルタの出力により動作する第１のスピーカと、
前記第２および第４のＦＩＲフィルタの出力により動作
する第２のスピーカとを有し、前記第１および第２のス
ピーカの各出力によって立体音場が再生されるように前
記第１乃至第４のＦＩＲフィルタのフィルタ係数を設定
する立体音場再生方式であって、所定の聴取位置に設置
されるダミーヘッドマイクと、所望波と相関の強い参照
信号を発生し、該参照信号に基づいて前記第１および第
２のスピーカをそれぞれ動作させ、該第１および第２の
スピーカの出力を測定した前記ダミーヘッドマイクの出
力と前記参照信号との誤差が最小となるような、メイン
パス信号に関する伝達関数とクロストークパス信号に関
する伝達関数を求める伝達関数測定手段と、前記伝達関
数測定手段によって求められたメインパス信号に関する
伝達関数およびクロストークパス信号に関する伝達関数
と、該メインパス信号に関する伝達関数の逆関数とをパ
ラメータとして、ＦＦＴによる畳み込み演算処理を用い
て前記第１乃至第４のＦＩＲフィルタのフィルタ係数を
設定するフィルタ係数計算手段とを有することを特徴と
する。
め、本発明の立体音場再生方式は、立体音場を生成する
ための第１および第２のオーディオ信号と、前記第１の
オーディオ信号入力とする第１および第２のＦＩＲフィ
ルタと、前記第２のオーディオ信号を入力とする第３お
よび第４のＦＩＲフィルタと、前記第１および第３のＦ
ＩＲフィルタの出力により動作する第１のスピーカと、
前記第２および第４のＦＩＲフィルタの出力により動作
する第２のスピーカとを有し、前記第１および第２のス
ピーカの各出力によって立体音場が再生されるように前
記第１乃至第４のＦＩＲフィルタのフィルタ係数を設定
する立体音場再生方式であって、所定の聴取位置に設置
されるダミーヘッドマイクと、所望波と相関の強い参照
信号を発生し、該参照信号に基づいて前記第１および第
２のスピーカをそれぞれ動作させ、該第１および第２の
スピーカの出力を測定した前記ダミーヘッドマイクの出
力と前記参照信号との誤差が最小となるような、メイン
パス信号に関する伝達関数とクロストークパス信号に関
する伝達関数を求める伝達関数測定手段と、前記伝達関
数測定手段によって求められたメインパス信号に関する
伝達関数およびクロストークパス信号に関する伝達関数
と、該メインパス信号に関する伝達関数の逆関数とをパ
ラメータとして、ＦＦＴによる畳み込み演算処理を用い
て前記第１乃至第４のＦＩＲフィルタのフィルタ係数を
設定するフィルタ係数計算手段とを有することを特徴と
する。
【００１４】上記の場合、前記ダミーヘッドマイクが、
聴取者の左右の耳にかかる音圧をそれぞれ測定する第１
および第２のダミーヘッドマイクからなり、前記伝達関
数測定手段が、前記参照信号を発生する参照信号発生手
段と、前記参照信号発生手段から発生された参照信号を
入力とする第１および第２のアダプティブフィルタと、
前記第１のダミーヘッドマイクの出力を一方の入力と
し、前記第１のアダプティブフィルタの出力を他方の入
力とし、これら入力の誤差を出力する第１の加算手段
と、前記第２のダミーヘッドマイクの出力を一方の入力
とし、前記第２のアダプティブフィルタの出力を他方の
入力とし、これら入力の誤差を出力する第２の加算手段
と、前記第１の加算手段から出力される誤差の２乗平均
誤差を求め、該２乗平均誤差が最小となるように前記第
１のアダプティブフィルタのフィルタ係数を更新する第
１の最小２乗法処理手段と、前記第２の加算手段から出
力される誤差の２乗平均誤差を求め、該２乗平均誤差が
最小となるように前記第２のアダプティブフィルタのフ
ィルタ係数を更新する第２の最小２乗法処理手段と、を
有し、前記フィルタ係数計算手段が、前記第１および第
２のアダプティブフィルタのフィルタ係数に基づいて前
記第１乃至第４のＦＩＲフィルタのフィルタ係数を設定
するように構成してもよい。
聴取者の左右の耳にかかる音圧をそれぞれ測定する第１
および第２のダミーヘッドマイクからなり、前記伝達関
数測定手段が、前記参照信号を発生する参照信号発生手
段と、前記参照信号発生手段から発生された参照信号を
入力とする第１および第２のアダプティブフィルタと、
前記第１のダミーヘッドマイクの出力を一方の入力と
し、前記第１のアダプティブフィルタの出力を他方の入
力とし、これら入力の誤差を出力する第１の加算手段
と、前記第２のダミーヘッドマイクの出力を一方の入力
とし、前記第２のアダプティブフィルタの出力を他方の
入力とし、これら入力の誤差を出力する第２の加算手段
と、前記第１の加算手段から出力される誤差の２乗平均
誤差を求め、該２乗平均誤差が最小となるように前記第
１のアダプティブフィルタのフィルタ係数を更新する第
１の最小２乗法処理手段と、前記第２の加算手段から出
力される誤差の２乗平均誤差を求め、該２乗平均誤差が
最小となるように前記第２のアダプティブフィルタのフ
ィルタ係数を更新する第２の最小２乗法処理手段と、を
有し、前記フィルタ係数計算手段が、前記第１および第
２のアダプティブフィルタのフィルタ係数に基づいて前
記第１乃至第４のＦＩＲフィルタのフィルタ係数を設定
するように構成してもよい。
【００１５】
【作用】従来の立体音場再生方式では、ＦＦＴによる畳
み込み演算処理を行う際のパラメータには、メインパス
信号に関する遅延伝達関数が用いられてないため、これ
により算出されたＦＩＲフィルタのフィルタ係数は遅延
誤差を含んだものとなっていた。
み込み演算処理を行う際のパラメータには、メインパス
信号に関する遅延伝達関数が用いられてないため、これ
により算出されたＦＩＲフィルタのフィルタ係数は遅延
誤差を含んだものとなっていた。
【００１６】上記遅延誤差について検討した結果、後に
詳述するようにメインパス信号に関する伝達関数の逆数
を用いることにより、遅延誤差を補正できることがわか
った。
詳述するようにメインパス信号に関する伝達関数の逆数
を用いることにより、遅延誤差を補正できることがわか
った。
【００１７】本発明の立体音場再生方式では、ＦＦＴに
よる畳み込み演算処理を行う際のパラメータとして、メ
インパス信号に関する伝達関数およびクロストークパス
信号に関する伝達関数に加えて、メインパス信号に関す
る伝達関数の逆関数が用いられる。このメインパス信号
に関する伝達関数の逆関数をパラメータに加えることに
より、フィルタ係数を求める際の遅延誤差が補正され、
ＦＩＲフィルタのフィルタ係数として遅延誤差が補正さ
れたフィルタ係数が設定される。このようにＦＩＲフィ
ルタのフィルタ係数として遅延誤差が補正されたフィル
タ係数が設定される本発明の方式においては、従来の方
式のものに比べて、正確にクロストーク信号をキャンセ
ルでき、正確にパス信号を聴取者に伝達することができ
る。
よる畳み込み演算処理を行う際のパラメータとして、メ
インパス信号に関する伝達関数およびクロストークパス
信号に関する伝達関数に加えて、メインパス信号に関す
る伝達関数の逆関数が用いられる。このメインパス信号
に関する伝達関数の逆関数をパラメータに加えることに
より、フィルタ係数を求める際の遅延誤差が補正され、
ＦＩＲフィルタのフィルタ係数として遅延誤差が補正さ
れたフィルタ係数が設定される。このようにＦＩＲフィ
ルタのフィルタ係数として遅延誤差が補正されたフィル
タ係数が設定される本発明の方式においては、従来の方
式のものに比べて、正確にクロストーク信号をキャンセ
ルでき、正確にパス信号を聴取者に伝達することができ
る。
【００１８】また、本発明では、伝達関数測定手段にお
いて求められるメインパス信号に関する伝達関数とクロ
ストークパス信号に関する伝達関数が、設置環境によっ
て変化する室内音響特性に最適なものとなるように求め
られるので、ＦＩＲフィルタのフィルタ係数を設置環境
に応じた値に設定することができる。よって、フィルタ
特性を設置環境によって変化する室内音響特性に最適な
ものとすることができ、正確にメインパス信号を聴取者
に伝達することができる。
いて求められるメインパス信号に関する伝達関数とクロ
ストークパス信号に関する伝達関数が、設置環境によっ
て変化する室内音響特性に最適なものとなるように求め
られるので、ＦＩＲフィルタのフィルタ係数を設置環境
に応じた値に設定することができる。よって、フィルタ
特性を設置環境によって変化する室内音響特性に最適な
ものとすることができ、正確にメインパス信号を聴取者
に伝達することができる。
【００１９】
【実施例】次に、本発明の実施例について図面を参照し
て説明する。
て説明する。
【００２０】図１は、本発明の一実施例の立体音場再生
方式を説明するための図である。
方式を説明するための図である。
【００２１】本実施例の立体音場再生方式は２チャンネ
ルスピーカ立体音場再生方式であって、伝達関数を測定
する際に生じる遅延誤差を補正できるようにしたもので
ある。さらには、室内音響特性が設置環境により変化し
ても、その変化に応じて最適なＦＩＲフィルタのフィル
タ係数を決定することをできるようにしたものである。
以下、この方式のシステム構成について説明する。
ルスピーカ立体音場再生方式であって、伝達関数を測定
する際に生じる遅延誤差を補正できるようにしたもので
ある。さらには、室内音響特性が設置環境により変化し
ても、その変化に応じて最適なＦＩＲフィルタのフィル
タ係数を決定することをできるようにしたものである。
以下、この方式のシステム構成について説明する。
【００２２】図１において、ＨLLは左側スピーカ１０か
ら左耳に伝わる伝達関数、ＨRLは右側スピーカ９から左
耳に伝わる伝達関数、ＨRRは右側スピーカ９から右耳に
伝わる伝達関数、ＨLRは左側スピーカ１０から右耳に伝
わる伝達関数を示す。
ら左耳に伝わる伝達関数、ＨRLは右側スピーカ９から左
耳に伝わる伝達関数、ＨRRは右側スピーカ９から右耳に
伝わる伝達関数、ＨLRは左側スピーカ１０から右耳に伝
わる伝達関数を示す。
【００２３】フィルタ構成は従来の方式のものと同様、
左側スピーカ９用の入力伝達関数（オーディオデジタル
信号）Ｌをそれぞれ入力とするＦＩＲフィルタＸLL１，
ＸLR２と、右側スピーカ１０用の入力伝達関数（オーデ
ィオデジタル信号）Ｒをそれぞれ入力とするＦＩＲフィ
ルタＸRL３，ＸRR４とにより構成され、ＦＩＲフィルタ
ＸLL１およびＦＩＲフィルタＸRL３の各出力が、Ｄ／Ａ
コンバータ５、アンプ７を順次介して左側スピーカ９の
入力となっており、ＦＩＲフィルタＸLR２およびＦＩＲ
フィルタＸRR４の各出力が、Ｄ／Ａコンバータ６、アン
プ８を順次介して右側スピーカ１０の入力となってい
る。
左側スピーカ９用の入力伝達関数（オーディオデジタル
信号）Ｌをそれぞれ入力とするＦＩＲフィルタＸLL１，
ＸLR２と、右側スピーカ１０用の入力伝達関数（オーデ
ィオデジタル信号）Ｒをそれぞれ入力とするＦＩＲフィ
ルタＸRL３，ＸRR４とにより構成され、ＦＩＲフィルタ
ＸLL１およびＦＩＲフィルタＸRL３の各出力が、Ｄ／Ａ
コンバータ５、アンプ７を順次介して左側スピーカ９の
入力となっており、ＦＩＲフィルタＸLR２およびＦＩＲ
フィルタＸRR４の各出力が、Ｄ／Ａコンバータ６、アン
プ８を順次介して右側スピーカ１０の入力となってい
る。
【００２４】上記フィルタ構成部における各ＦＩＲフィ
ルタＸLL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数を
決定するための手段として、信号検出部１２、伝達関数
測定部１３およびフィルタ計算部１４が設けられてい
る。信号検出部１２は、所定の聴取位置における聴取者
の左右の耳にかかる音圧ＰL，ＰRを測定するものであ
る。伝達関数測定部１３は、信号検出部１２にて検出さ
れた音圧を基に伝達関数を求めるものである。フィルタ
計算部１４は、伝達関数測定部１３にて求められた伝達
関数を基に各ＦＩＲフィルタＸLL１，ＸLR２，ＸRL３，
ＸRR４のフィルタ係数を計算するものである。
ルタＸLL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数を
決定するための手段として、信号検出部１２、伝達関数
測定部１３およびフィルタ計算部１４が設けられてい
る。信号検出部１２は、所定の聴取位置における聴取者
の左右の耳にかかる音圧ＰL，ＰRを測定するものであ
る。伝達関数測定部１３は、信号検出部１２にて検出さ
れた音圧を基に伝達関数を求めるものである。フィルタ
計算部１４は、伝達関数測定部１３にて求められた伝達
関数を基に各ＦＩＲフィルタＸLL１，ＸLR２，ＸRL３，
ＸRR４のフィルタ係数を計算するものである。
【００２５】上記のようなシステム構成をとる本実施例
の立体音場再生方式では、設置場所において信号検出部
１２、伝達関数測定部１３およびフィルタ計算部１４に
よって各ＦＩＲフィルタＸLL１，ＸLR２，ＸRL３，ＸRR
４のフィルタ係数が決定れる。
の立体音場再生方式では、設置場所において信号検出部
１２、伝達関数測定部１３およびフィルタ計算部１４に
よって各ＦＩＲフィルタＸLL１，ＸLR２，ＸRL３，ＸRR
４のフィルタ係数が決定れる。
【００２６】以下、各ＦＩＲフィルタＸLL１，ＸLR２，
ＸRL３，ＸRR４のフィルタ係数の求め方について詳しく
説明する。
ＸRL３，ＸRR４のフィルタ係数の求め方について詳しく
説明する。
【００２７】聴取者の左右の耳にかかるＰL，ＰRは、フ
ィルタ係数ＸLL，ＸLR，ＸRL，ＸRRがクロストークパス
信号をキャンセルしてメインパス信号を伝えることがで
きる係数とすると、入力伝達関数Ｌ，Ｒ、伝達関数
ＨLL，ＨRL，ＨRR，ＨLRを用いて以下の式で表すことが
できる。
ィルタ係数ＸLL，ＸLR，ＸRL，ＸRRがクロストークパス
信号をキャンセルしてメインパス信号を伝えることがで
きる係数とすると、入力伝達関数Ｌ，Ｒ、伝達関数
ＨLL，ＨRL，ＨRR，ＨLRを用いて以下の式で表すことが
できる。
【００２８】ＰL＝（ＸLL×ＨLL＋ＸLR×ＨRL）Ｌ ＋（ＸRR×ＨRL＋ＸRL×ＨLL）Ｒ （５） ＰR＝（ＸLL×ＨLR＋ＸLR×ＨRR）Ｌ ＋（ＸRR×ＨRR＋ＸRL×ＨLR）Ｒ （６） ここで、メインパス信号のみが両耳に伝わるように上記
式（５），（６）の係数項を以下の式のように置く。
式（５），（６）の係数項を以下の式のように置く。
【００２９】 ＸLL×ＨLL＋ＸLR×ＨRL＝１ （７） ＸLL×ＨLR＋ＸLR×ＨRR＝０ （８） ＸRR×ＨRL＋ＸRL×ＨLL＝０ （９） ＸRR×ＨRR＋ＸRL×ＨLR＝１ （１０） 式（７），（１０）の左辺「１」は聴取者に伝わるメイ
ンパス信号の理想伝達関数、式（８），（９）の左辺
「０」はクロストークパス信号の理想伝達関数である。
ＸLL，ＸLRとＸRL，ＸRRとは同様にして求めることとな
るので、以下ＸLL，ＸLRについての算出について説明す
る。
ンパス信号の理想伝達関数、式（８），（９）の左辺
「０」はクロストークパス信号の理想伝達関数である。
ＸLL，ＸLRとＸRL，ＸRRとは同様にして求めることとな
るので、以下ＸLL，ＸLRについての算出について説明す
る。
【００３０】上記式（７），（８）を行列により表す
と、
と、
【００３１】
【数１】 となる。上記式（１１）を以下の式のように置く。
【００３２】 ＨＸ＝１ （１２） 上記式（１２）の両辺に行列Ｈの逆行列Ｈ-1をかけて Ｈ-1ＨＸ＝Ｈ-1 （１３） を得る。この式（１３）より遅延パラメータＨ-1Ｈが発
生したことがわかる。よって、ＸLL，ＸLRは以下の式に
より求めることができる（ＸRL，ＸRRも同様にして求め
ることができる）。
生したことがわかる。よって、ＸLL，ＸLRは以下の式に
より求めることができる（ＸRL，ＸRRも同様にして求め
ることができる）。
【００３３】 ＸLL＝（１÷ｄｅｔ Ｈ）×（ＨRR÷（ＨLL×ＨRR）） （１４） ＸLR＝（１÷ｄｅｔ Ｈ）×（−ＨLR÷（ＨLL×ＨRR）） （１５） ｄｅｔ Ｈ＝（ＨLL×ＨRR）÷（ＨLL×ＨRR） −（ＨLR×ＨRL）÷（ＨLL×ＨRR） （１６） 上記式（１４）〜（１６）を上述の式（７）に代入し、
メインパス信号のスピーカから聴取者までの伝達関数を
求める。
メインパス信号のスピーカから聴取者までの伝達関数を
求める。
【００３４】 ｄｅｔ Ｈ÷ｄｅｔ Ｈ＝１ （１７） ここで、各伝達関数の逆関数は最小２乗解であることか
ら上記式（１７）より上述の式（７）は以下のように表
すことができる。
ら上記式（１７）より上述の式（７）は以下のように表
すことができる。
【００３５】 ＸLL×ＨLL＋ＸLR×ＨRL＝ｄｅｔ Ｈ÷ｄｅｔ Ｈ （１８） 上記式（１８）と上述の式（８）よりＸLL，ＸLRを求め
ることができる（ＸRL，ＸRRについても同様にして求め
ることができる）。
ることができる（ＸRL，ＸRRについても同様にして求め
ることができる）。
【００３６】各フィルタ係数ＸLL，ＸLR，ＸRL，Ｘ
RRは、伝達関数ＨLL，ＨRL，ＨRR，ＨLRと伝達関数
ＨLL，ＨRRの逆関数（ＨLL）-1，（ＨRR）-1をＦＦＴに
よる畳み込み演算処理を行うことにより求められること
から、以下の式で表される。
RRは、伝達関数ＨLL，ＨRL，ＨRR，ＨLRと伝達関数
ＨLL，ＨRRの逆関数（ＨLL）-1，（ＨRR）-1をＦＦＴに
よる畳み込み演算処理を行うことにより求められること
から、以下の式で表される。
【００３７】 ＸLL＝Ｘ×（Ｘ2）-1×（ＨLL×ＨRR）-1×ＨRR （１９） ＸLR＝Ｘ×（Ｘ2）-1×（ＨLL×ＨRR）-1×（−ＨLR） （２０） ＸRR＝Ｘ×（Ｘ2）-1×（ＨLL×ＨRR）-1×ＨLL （２１） ＸRL＝Ｘ×（Ｘ2）-1×（ＨLL×ＨRR）-1×（−ＨRL） （２２） Ｘ＝ＨLL×ＨRR×（ＨLL×ＨRR）-1 −ＨLR×ＨRL×（ＨLL×ＨRR）-1 （２３） （Ｘ2）-1＝｛〔ＨLL×ＨRR×（ＨLL×ＨRR）-1 −ＨLR×ＨRL×（ＨLL×ＨRR）-1〕2｝-1 （２４） 以上のようにして求められた式（１９）〜（２４）を用
いて各フィルタ係数Ｘ LL，ＸLR，ＸRL，ＸRRを求めれ
ば、従来方式に比べてクロストークパス信号を減衰さ
せ、メインパス信号を正確に伝送することが可能とな
る。
いて各フィルタ係数Ｘ LL，ＸLR，ＸRL，ＸRRを求めれ
ば、従来方式に比べてクロストークパス信号を減衰さ
せ、メインパス信号を正確に伝送することが可能とな
る。
【００３８】なお、逆関数（Ｘ2）-1，（ＨLL×ＨRR）
-1は、逆フィルタの定理により、レビンソンアルゴリズ
ムなどを用いて次のようにして求めることができる。理
想インパルス応答Ｙ、伝達関数Ｈ、求める逆関数Ｚの関
係をマトリックスで表現すると、
-1は、逆フィルタの定理により、レビンソンアルゴリズ
ムなどを用いて次のようにして求めることができる。理
想インパルス応答Ｙ、伝達関数Ｈ、求める逆関数Ｚの関
係をマトリックスで表現すると、
【００３９】
【数２】 のようになる。上記式（２５）の両辺に転置マトリクス
をかけ、逆関数Ｚの最小２乗解を求めると、 ＨtＹ＝ＨtＨＺ （２６） となる。相互相関ベクトル Ｐ＝ＨtＹ ，自己相関マト
リクス Ｒ＝ＨtＨと置くと Ｐ＝ＲＺ と表すことがで
き、逆関数Ｚは以下の式により求めることができる。
をかけ、逆関数Ｚの最小２乗解を求めると、 ＨtＹ＝ＨtＨＺ （２６） となる。相互相関ベクトル Ｐ＝ＨtＹ ，自己相関マト
リクス Ｒ＝ＨtＨと置くと Ｐ＝ＲＺ と表すことがで
き、逆関数Ｚは以下の式により求めることができる。
【００４０】 Ｚ＝Ｒ-1Ｐ （２７） この逆関数Ｚの最小２乗解が逆関数（Ｘ2）-1，（ＨLL
×ＨRR）-1に相当し、伝達関数ＨがＸ2、ＨLL×ＨRRに
相当する。
×ＨRR）-1に相当し、伝達関数ＨがＸ2、ＨLL×ＨRRに
相当する。
【００４１】次に、上述の式（１９）〜（２４）を用い
て各フィルタ係数ＸLL，ＸLR，ＸRL，ＸRRを決定する手
段である、信号検出部１２、伝達関数測定部１３および
フィルタ計算部１４について詳しく説明する。
て各フィルタ係数ＸLL，ＸLR，ＸRL，ＸRRを決定する手
段である、信号検出部１２、伝達関数測定部１３および
フィルタ計算部１４について詳しく説明する。
【００４２】図２は、図１に示した信号検出部１２、伝
達関数測定部１３およびフィルタ計算部１４の内部構成
を示す図である。
達関数測定部１３およびフィルタ計算部１４の内部構成
を示す図である。
【００４３】信号検出部１２は、所定の聴取位置におけ
る聴取者の左右の耳にかかる音圧Ｐ L，ＰRを測定するた
めのダミーヘッドマイク１５L，１５Rと、これらダミー
ヘッドマイク１５L，１５Rによって検出されたアナログ
信号をＡ／Ｄ変換するＡ／Ｄコンバータ１６，１７とに
より構成されている。
る聴取者の左右の耳にかかる音圧Ｐ L，ＰRを測定するた
めのダミーヘッドマイク１５L，１５Rと、これらダミー
ヘッドマイク１５L，１５Rによって検出されたアナログ
信号をＡ／Ｄ変換するＡ／Ｄコンバータ１６，１７とに
より構成されている。
【００４４】伝達関数測定部１３は、系のなかで参照信
号を発生するＭ系列（Maximum Period sequence の略）
ノイズ部２２と、アダプティブフィルタ１８，２０と、
最小２乗法処理部１９，２１と、加算器２３，２４とに
より構成されている。アダプティブフィルタ１８，２０
は、それぞれＭ系列ノイズ部２２からの参照信号を入力
として出力信号ｏを出力する。加算器２３，２４は、Ａ
／Ｄコンバータ１６，１７のそれぞれの出力を理想出力
ｄとして一方の入力とし、アダプティブフィルタ１８，
２０の各出力信号ｏを他方の入力として、理想出力ｄと
出力信号ｏとの誤差ｅをそれぞれ出力する。最小２乗法
処理部１９，２１は、加算器２３，２４から出力される
誤差ｅの２乗平均誤差がそれぞれ最小となるように最小
２乗法を用いてアダプティブフィルタ１８，２０のフィ
ルタ係数をそれぞれ更新する。
号を発生するＭ系列（Maximum Period sequence の略）
ノイズ部２２と、アダプティブフィルタ１８，２０と、
最小２乗法処理部１９，２１と、加算器２３，２４とに
より構成されている。アダプティブフィルタ１８，２０
は、それぞれＭ系列ノイズ部２２からの参照信号を入力
として出力信号ｏを出力する。加算器２３，２４は、Ａ
／Ｄコンバータ１６，１７のそれぞれの出力を理想出力
ｄとして一方の入力とし、アダプティブフィルタ１８，
２０の各出力信号ｏを他方の入力として、理想出力ｄと
出力信号ｏとの誤差ｅをそれぞれ出力する。最小２乗法
処理部１９，２１は、加算器２３，２４から出力される
誤差ｅの２乗平均誤差がそれぞれ最小となるように最小
２乗法を用いてアダプティブフィルタ１８，２０のフィ
ルタ係数をそれぞれ更新する。
【００４５】フィルタ計算部１４は、最小２乗法処理部
１９，２１によって更新されたアダプティブフィルタ１
８，２０のフィルタ係数から各伝達関数（ＨLL，ＨRL，
ＨRR，ＨLR）を得、ＦＦＴによる畳み込み演算処理、お
よび逆フィルタを求めるためのアルゴリズムによりＦＩ
ＲフィルタＸLL１，ＸLR２，ＸRL３，ＸRR４に設定する
フィルタ係数の計算を上述した式（１９）〜（２４）を
用いて行う。
１９，２１によって更新されたアダプティブフィルタ１
８，２０のフィルタ係数から各伝達関数（ＨLL，ＨRL，
ＨRR，ＨLR）を得、ＦＦＴによる畳み込み演算処理、お
よび逆フィルタを求めるためのアルゴリズムによりＦＩ
ＲフィルタＸLL１，ＸLR２，ＸRL３，ＸRR４に設定する
フィルタ係数の計算を上述した式（１９）〜（２４）を
用いて行う。
【００４６】なお、図中各スピーカ９，１０の入力につ
いては、Ｄ／Ａコンバータ５およびアンプ７を介して、
またはＤ／Ａコンバータ６およびアンプ８を介して、伝
達関数ＨLL，ＨLRを求めるときには左側スピーカ９へ、
伝達関数ＨRR，ＨRLを求めるときには右側スピーカ１０
に入力されるよう構成されている。また、信号検出部１
２に聴取者の両耳に対応するダミーヘッドマイク１
５L，１５Rを設け、伝達関数測定部１３にアダプティブ
フィルタ１８，２０、最小２乗法処理部１９，２１、加
算器２３，２４を設けることにより、伝達関数ＨRRとＨ
RLを、または伝達関数ＨLLとＨLRを同時に求めることが
できるものとなっている。
いては、Ｄ／Ａコンバータ５およびアンプ７を介して、
またはＤ／Ａコンバータ６およびアンプ８を介して、伝
達関数ＨLL，ＨLRを求めるときには左側スピーカ９へ、
伝達関数ＨRR，ＨRLを求めるときには右側スピーカ１０
に入力されるよう構成されている。また、信号検出部１
２に聴取者の両耳に対応するダミーヘッドマイク１
５L，１５Rを設け、伝達関数測定部１３にアダプティブ
フィルタ１８，２０、最小２乗法処理部１９，２１、加
算器２３，２４を設けることにより、伝達関数ＨRRとＨ
RLを、または伝達関数ＨLLとＨLRを同時に求めることが
できるものとなっている。
【００４７】以下、信号検出部１２、伝達関数測定部１
３およびフィルタ計算部１４の具体的な動作について説
明する。ここでは、伝達関数ＨRRとＨRL、および伝達関
数Ｈ LLとＨLRは同様にして求めれることから、ここで
は、伝達関数ＨLLとＨLRを求める場合について説明す
る。
３およびフィルタ計算部１４の具体的な動作について説
明する。ここでは、伝達関数ＨRRとＨRL、および伝達関
数Ｈ LLとＨLRは同様にして求めれることから、ここで
は、伝達関数ＨLLとＨLRを求める場合について説明す
る。
【００４８】伝達関数ＨLLとＨLRを求める場合は、Ｍ系
列ノイズ部２２から参照信号がＤ／Ａコンバータ５およ
びアンプ７を介して左側スピーカ９に入力されるととも
に、アダティブフィルタ１８，２０のそれぞれに入力さ
れる。
列ノイズ部２２から参照信号がＤ／Ａコンバータ５およ
びアンプ７を介して左側スピーカ９に入力されるととも
に、アダティブフィルタ１８，２０のそれぞれに入力さ
れる。
【００４９】右側スピーカ９に参照信号が入力される
と、左側スピーカ９から信号が出力され、これが信号検
出部１２のダミーヘッドマイク１５L，１５Rによってそ
れぞれ検出される。すなわち、ダミーヘッドマイク１５
Lによって左側スピーカ９から聴取者の左の耳に伝わる
信号が検出され、ダミーヘッドマイク１５Rによって左
側スピーカ９から聴取者の右の耳に伝わる信号が検出さ
れる。
と、左側スピーカ９から信号が出力され、これが信号検
出部１２のダミーヘッドマイク１５L，１５Rによってそ
れぞれ検出される。すなわち、ダミーヘッドマイク１５
Lによって左側スピーカ９から聴取者の左の耳に伝わる
信号が検出され、ダミーヘッドマイク１５Rによって左
側スピーカ９から聴取者の右の耳に伝わる信号が検出さ
れる。
【００５０】右側スピーカ９からの信号がダミーヘッド
マイク１５L，１５Rによってそれぞれ検出されると、ダ
ミーヘッドマイク１５L，１５Rの出力アナログ信号がＡ
／Ｄコンバータ１６，１７によってＡ／Ｄ変換され、理
想出力ｄとして伝達関数測定部１３の加算器２３，２４
へそれぞれ入力される。
マイク１５L，１５Rによってそれぞれ検出されると、ダ
ミーヘッドマイク１５L，１５Rの出力アナログ信号がＡ
／Ｄコンバータ１６，１７によってＡ／Ｄ変換され、理
想出力ｄとして伝達関数測定部１３の加算器２３，２４
へそれぞれ入力される。
【００５１】各Ａ／Ｄコンバータ１６，１７から理想出
力ｄが加算器２３，２４へそれぞれ入力されると、加算
器２３では、Ａ／Ｄコンバータ１６から出力された理想
出力ｄと、Ｍ系列ノイズ２２からの参照信号を入力とし
てアダプティブフィルタ１８から出力される出力信号ｏ
との差分が求められ、これが誤差ｅとして最小２乗法処
理部１９に入力される。同様にして、加算器２４では、
Ａ／Ｄコンバータ１６からの理想出力ｄとアダプティブ
フィルタ１８からの出力信号ｏとの差分が求められ、こ
れが誤差ｅとして最小２乗法処理部２１に入力される。
力ｄが加算器２３，２４へそれぞれ入力されると、加算
器２３では、Ａ／Ｄコンバータ１６から出力された理想
出力ｄと、Ｍ系列ノイズ２２からの参照信号を入力とし
てアダプティブフィルタ１８から出力される出力信号ｏ
との差分が求められ、これが誤差ｅとして最小２乗法処
理部１９に入力される。同様にして、加算器２４では、
Ａ／Ｄコンバータ１６からの理想出力ｄとアダプティブ
フィルタ１８からの出力信号ｏとの差分が求められ、こ
れが誤差ｅとして最小２乗法処理部２１に入力される。
【００５２】上述のようにして各加算器２３，２４から
誤差ｅが最小２乗法処理部１９，２１に入力されると、
最小２乗法処理部１９は、その入力される誤差ｅの２乗
平均誤差が最小となるように最小２乗法を用いてアダプ
ティブフィルタ１８のフィルタ係数を更新させる。同様
にして、最小２乗法処理部２１は、その入力される誤差
ｅが最小となるように最小２乗法を用いてアダプティブ
フィルタ２０のフィルタ係数を更新させる。このように
して、更新されたアダプティブフィルタ１８のフィルタ
係数から伝達関数ＨLRを得、更新されたアダプティブフ
ィルタ２０のフィルタ係数から伝達関数ＨLLを得る。
誤差ｅが最小２乗法処理部１９，２１に入力されると、
最小２乗法処理部１９は、その入力される誤差ｅの２乗
平均誤差が最小となるように最小２乗法を用いてアダプ
ティブフィルタ１８のフィルタ係数を更新させる。同様
にして、最小２乗法処理部２１は、その入力される誤差
ｅが最小となるように最小２乗法を用いてアダプティブ
フィルタ２０のフィルタ係数を更新させる。このように
して、更新されたアダプティブフィルタ１８のフィルタ
係数から伝達関数ＨLRを得、更新されたアダプティブフ
ィルタ２０のフィルタ係数から伝達関数ＨLLを得る。
【００５３】上述のような処理を伝達関数ＨRRとＨRLに
ついても行い、伝達関数ＨLL，ＨLR，ＨRR，ＨRLが全て
求められると、フィルタ計算部１４が、その求められた
伝達関数ＨLL，ＨLR，ＨRR，ＨRLを基に、ＦＦＴによる
畳み込み演算処理、および逆フィルタを求めるためのア
ルゴリズムによりＦＩＲフィルタＸLL１，ＸLR２，Ｘ RL
３，ＸRR４に設定するフィルタ係数の計算を上述した式
（１９）〜（２４）を用いて行う。各ＦＩＲフィルタＸ
LL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数を計算す
ると、フィルタ計算部１４は、その計算されたフィルタ
係数をそれぞれ図１に示す各ＦＩＲフィルタＸLL１，Ｘ
LR２，ＸRL３，ＸRR４へ転送して設定する。
ついても行い、伝達関数ＨLL，ＨLR，ＨRR，ＨRLが全て
求められると、フィルタ計算部１４が、その求められた
伝達関数ＨLL，ＨLR，ＨRR，ＨRLを基に、ＦＦＴによる
畳み込み演算処理、および逆フィルタを求めるためのア
ルゴリズムによりＦＩＲフィルタＸLL１，ＸLR２，Ｘ RL
３，ＸRR４に設定するフィルタ係数の計算を上述した式
（１９）〜（２４）を用いて行う。各ＦＩＲフィルタＸ
LL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数を計算す
ると、フィルタ計算部１４は、その計算されたフィルタ
係数をそれぞれ図１に示す各ＦＩＲフィルタＸLL１，Ｘ
LR２，ＸRL３，ＸRR４へ転送して設定する。
【００５４】以上のようにして各ＦＩＲフィルタＸ
LL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数が設定さ
れることにより、従来方式に比べてクロストークパス信
号を減衰させ、メインパス信号を正確に伝送することが
可能となる。
LL１，ＸLR２，ＸRL３，ＸRR４のフィルタ係数が設定さ
れることにより、従来方式に比べてクロストークパス信
号を減衰させ、メインパス信号を正確に伝送することが
可能となる。
【００５５】なお、上述した実施例では、信号検出部１
２に聴取者の両耳に対応するダミーヘッドマイク１
５L，１５Rを設け、伝達関数測定部１３にアダプティブ
フィルタ１８，２０、最小２乗法処理部１９，２１、加
算器２３，２４を設けることにより、伝達関数ＨRRとＨ
RLを、または伝達関数ＨLLとＨLRを同時に求めるように
した構成のものを説明したが、ダミーヘッドマイク、ア
ダプティブフィルタ、最小２乗法処理部、加算器は１つ
であってもよく、この場合には、各伝達関数はそれぞれ
別々に求められることとなる。
２に聴取者の両耳に対応するダミーヘッドマイク１
５L，１５Rを設け、伝達関数測定部１３にアダプティブ
フィルタ１８，２０、最小２乗法処理部１９，２１、加
算器２３，２４を設けることにより、伝達関数ＨRRとＨ
RLを、または伝達関数ＨLLとＨLRを同時に求めるように
した構成のものを説明したが、ダミーヘッドマイク、ア
ダプティブフィルタ、最小２乗法処理部、加算器は１つ
であってもよく、この場合には、各伝達関数はそれぞれ
別々に求められることとなる。
【００５６】
【発明の効果】本発明は、以上説明したように構成され
ているので、以下に記載するような効果を奏する。
ているので、以下に記載するような効果を奏する。
【００５７】本発明の立体音場再生方式においては、設
定されたＦＩＲフィルタのフィルタ係数は遅延誤差を補
正したものとなっているので、従来の方式のものに比べ
て、正確にクロストーク信号をキャンセルでき、臨場感
が向上するという効果がある。
定されたＦＩＲフィルタのフィルタ係数は遅延誤差を補
正したものとなっているので、従来の方式のものに比べ
て、正確にクロストーク信号をキャンセルでき、臨場感
が向上するという効果がある。
【００５８】さらには、フィルタ特性を設置環境によっ
て変化する室内音響特性に最適なものとすることができ
るので、室内音響特性が大きく異なる場合においても、
所望の臨場感を得ることができるという効果がある。
て変化する室内音響特性に最適なものとすることができ
るので、室内音響特性が大きく異なる場合においても、
所望の臨場感を得ることができるという効果がある。
【図１】本発明の一実施例の立体音場再生方式を説明す
るための図である。
るための図である。
【図２】図１に示した信号検出部１２、伝達関数測定部
１３およびフィルタ計算部１４の内部構成を示す図であ
る。
１３およびフィルタ計算部１４の内部構成を示す図であ
る。
【図３】従来の立体音場再生方式を説明するための図で
ある。
ある。
１〜４ ＦＩＲフィルタ ５，６ Ｄ／Ａコンバータ ７，８ アンプ ９，１０ スピーカ １２ 信号検出部 １３ 伝達関数測定部 １４ フィルタ計算部 １５L,１５R ダミーヘッドマイク １６，１７ Ａ／Ｄコンバータ １８，２０ アダプティブフィルタ １９，２１ 最小２乗法処理部 ２２ Ｍ系列ノイズ部
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 ＦＩ Ｈ０４Ｒ 5/033 Ｇ１０Ｋ 15/00 Ｌ
Claims (2)
- 【請求項１】 立体音場を生成するための第１および第
２のオーディオ信号と、前記第１のオーディオ信号入力
とする第１および第２のＦＩＲフィルタと、前記第２の
オーディオ信号を入力とする第３および第４のＦＩＲフ
ィルタと、前記第１および第３のＦＩＲフィルタの出力
により動作する第１のスピーカと、前記第２および第４
のＦＩＲフィルタの出力により動作する第２のスピーカ
とを有し、前記第１および第２のスピーカの各出力によ
って立体音場が再生されるように前記第１乃至第４のＦ
ＩＲフィルタのフィルタ係数を設定する立体音場再生方
式であって、 所定の聴取位置に設置されるダミーヘッドマイクと、 所望波と相関の強い参照信号を発生し、該参照信号に基
づいて前記第１および第２のスピーカをそれぞれ動作さ
せ、該第１および第２のスピーカの出力を測定した前記
ダミーヘッドマイクの出力と前記参照信号との誤差が最
小となるような、メインパス信号に関する伝達関数とク
ロストークパス信号に関する伝達関数を求める伝達関数
測定手段と、 前記伝達関数測定手段によって求められたメインパス信
号に関する伝達関数およびクロストークパス信号に関す
る伝達関数と、該メインパス信号に関する伝達関数の逆
関数とをパラメータとして、ＦＦＴによる畳み込み演算
処理を用いて前記第１乃至第４のＦＩＲフィルタのフィ
ルタ係数を設定するフィルタ係数計算手段とを有するこ
とを特徴とする立体音場再生方式。 - 【請求項２】 請求項１に記載の立体音場再生方式にお
いて、 前記ダミーヘッドマイクが、 聴取者の左右の耳にかかる音圧をそれぞれ測定する第１
および第２のダミーヘッドマイクからなり、 前記伝達関数測定手段が、 前記参照信号を発生する参照信号発生手段と、 前記参照信号発生手段から発生された参照信号を入力と
する第１および第２のアダプティブフィルタと、 前記第１のダミーヘッドマイクの出力を一方の入力と
し、前記第１のアダプテ ィブフィルタの出力を他方の入
力とし、これら入力の誤差を出力する第１の加算手段
と、 前記第２のダミーヘッドマイクの出力を一方の入力と
し、前記第２のアダプティブフィルタの出力を他方の入
力とし、これら入力の誤差を出力する第２の加算手段
と、 前記第１の加算手段から出力される誤差の２乗平均誤差
を求め、該２乗平均誤差が最小となるように前記第１の
アダプティブフィルタのフィルタ係数を更新する第１の
最小２乗法処理手段と、 前記第２の加算手段から出力される誤差の２乗平均誤差
を求め、該２乗平均誤差が最小となるように前記第２の
アダプティブフィルタのフィルタ係数を更新する第２の
最小２乗法処理手段と、を有し、 前記フィルタ係数計算手段が、前記第１および第２のア
ダプティブフィルタのフィルタ係数に基づいて前記第１
乃至第４のＦＩＲフィルタのフィルタ係数を設定するこ
とを特徴とする立体音場再生方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7025210A JP2870440B2 (ja) | 1995-02-14 | 1995-02-14 | 立体音場再生方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7025210A JP2870440B2 (ja) | 1995-02-14 | 1995-02-14 | 立体音場再生方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08223700A JPH08223700A (ja) | 1996-08-30 |
| JP2870440B2 true JP2870440B2 (ja) | 1999-03-17 |
Family
ID=12159604
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7025210A Expired - Lifetime JP2870440B2 (ja) | 1995-02-14 | 1995-02-14 | 立体音場再生方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2870440B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9031268B2 (en) * | 2011-05-09 | 2015-05-12 | Dts, Inc. | Room characterization and correction for multi-channel audio |
| JP6565709B2 (ja) * | 2016-01-26 | 2019-08-28 | 株式会社Ｊｖｃケンウッド | 音像定位処理装置、及び音像定位処理方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS583640B2 (ja) * | 1977-08-29 | 1983-01-22 | 日本ビクター株式会社 | 音質補正回路 |
| JPH06165298A (ja) * | 1992-11-24 | 1994-06-10 | Nissan Motor Co Ltd | 音響再生装置 |
-
1995
- 1995-02-14 JP JP7025210A patent/JP2870440B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08223700A (ja) | 1996-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100739798B1 (ko) | 청취 위치를 고려한 2채널 입체음향 재생 방법 및 장치 | |
| US5018202A (en) | Electronic noise attenuation system | |
| JP3264489B2 (ja) | 音響再生装置 | |
| JP4780119B2 (ja) | 頭部伝達関数測定方法、頭部伝達関数畳み込み方法および頭部伝達関数畳み込み装置 | |
| EP1341399B1 (en) | Sound field control method and sound field control system | |
| US20070147636A1 (en) | Acoustics correcting apparatus | |
| KR20060026730A (ko) | 청취 위치를 고려한 2채널 가상 음향 재생 방법 및 장치 | |
| JP2001507879A (ja) | ステレオ音響エキスパンダ | |
| US20030076973A1 (en) | Sound signal processing method and sound reproduction apparatus | |
| JP3167259B2 (ja) | 音響再生装置 | |
| JPH08182100A (ja) | 音像定位方法および音像定位装置 | |
| JP4835151B2 (ja) | オーディオシステム | |
| JP2001285998A (ja) | 頭外音像定位装置 | |
| JP2006279863A (ja) | 頭部伝達関数の補正方法 | |
| EP0955789A2 (en) | Method and device for synthesizing a virtual sound source | |
| JP4171675B2 (ja) | 音場制御システム、および音場制御方法 | |
| JP2870440B2 (ja) | 立体音場再生方式 | |
| JP2714098B2 (ja) | 音響周波数特性の補正方法 | |
| JP3739438B2 (ja) | 音像定位方法及びその装置 | |
| JP4306815B2 (ja) | 線形予測係数を用いた立体音響処理装置 | |
| JP4845407B2 (ja) | リファレンスフィルタの生成方法 | |
| JPH08102999A (ja) | 立体音響再生装置 | |
| JP3254461B2 (ja) | 音響伝達特性予測方法およびその装置、音響装置 | |
| JP3393738B2 (ja) | 適応信号処理方法 | |
| JP5163685B2 (ja) | 頭部伝達関数測定方法、頭部伝達関数畳み込み方法および頭部伝達関数畳み込み装置 |