JP3724726B2

JP3724726B2 - 反響消去装置、反響消去方法、および反響消去プログラム

Info

Publication number: JP3724726B2
Application number: JP2002060892A
Authority: JP
Inventors: 朗中川
Original assignee: Nippon Telegraph and Telephone Corp; NTT Inc USA
Current assignee: NTT Inc; NTT Inc USA
Priority date: 2002-03-06
Filing date: 2002-03-06
Publication date: 2005-12-07
Anticipated expiration: 2022-03-06
Also published as: JP2003258688A

Description

【０００１】
【発明の属する技術分野】
この発明は、反響消去装置、反響消去方法、および反響消去プログラムに関し特に、多チャネル再生経路を有する通信会議装置の如き拡声通話装置において、ハウリングの原因および聴覚上の障害となる室内反響信号を消去する反響消去装置、反響消去方法、および反響消去プログラムに関する。
【０００２】
【従来の技術】
拡声通話装置にはハウリングの原因および聴覚上の障害となる室内反響信号を消去する反響消去装置が使用されているが、この従来例を図３を参照して説明する。
拡声通話において、相手の発話により得られる音声信号は再生信号として受話端子１１に受信されるが、この再生信号はそのままスピーカ１２により再生される場合と、スピーカ１２に入力されるに先だって、再生信号の振幅、或はパワーの大きさに応じて自動的に利得を調節するという様な何らかの加工を施された後にスピーカ１２から再生される場合とがある。従って、この明細書における再生信号ｘ₁（ｋ）とは、相手から送信された再生信号そのものの他に、何等かの処理が施された後の再生信号をも意味するものとする。ｋは離散時間を表す。反響消去装置１４は、再生信号ｘ₁（ｋ）がスピーカ１２から反響路１５を経て、マイクロホン１６に収音されて得られる反響信号ｙ₁（ｋ）を消去する装置である。反響信号ｙ₁（ｋ）は、時刻ｋにおける反響路１５のインパルス応答をｈ₁₁（ｋ、ｎ）として、
ｙ₁（ｋ）＝Σｈ₁₁（ｋ，ｎ）ｘ₁（ｋ−ｎ）・・・・・・（１）
Σはｎ＝０から（Ｌ−１）までの総和
の様な畳み込み演算で得られるもの、とモデル化することができる。Ｌはタップ数であり、反響路１５の残響時間に対応させて、予め設定しておく定数である。
【０００３】
先ず、再生信号蓄積ベクトル生成部１７において、再生信号ｘ₁（ｋ）を（Ｌ−１）時刻過去のものまでを蓄積しておく。蓄積された信号は、再生信号ベクトルｘ ₁（ｋ）、即ち

として出力される。但し、Ｔはベクトルの転置を表す。擬似反響信号生成部１８では、式（２）の再生信号ベクトルｘ ₁（ｋ）と、反響路推定部１９から得
られる擬似反響ベクトルｈ ₁₁＾（ｋ）との内積演算
ｙ₁＾（ｋ）＝ｈ ₁₁＾^T（ｋ）ｘ ₁（ｋ）・・・・・・（３）
を行い、その結果として、擬似反響信号ｙ₁＾（ｋ）を生成する。この内積演算は式（１）のような畳み込み演算と等価である。反響路推定部１９では、擬似反響信号生成部１８で用いる擬似反響路ベクトルｈ ₁₁＾（ｋ）を生成する。この
反響路推定に用いるもっとも一般的なアルゴリズムは、ＮＬＭＳアルゴリズム（学習同定法）である。ＮＬＭＳアルゴリズムでは時刻ｋにおける再生信号ベクトルｘ ₁（ｋ）と、残留反響信号ｅ₁（ｋ）、即ち差回路２１でマイクロホン１６
の出力ｙ₁（ｋ）から擬似反響信号ｙ₁＾（ｋ）を差し引いた信号
ｅ₁（ｋ）＝ｙ₁（ｋ）−ｙ₁＾（ｋ）・・・・・・（４）
とから、時刻ｋ＋１において用いる擬似反響路ベクトルｈ ₁₁＾（ｋ＋１）を次
式のように求める。
【０００４】

但し、μはステップサイズパラメータと呼ばれ０＜μ＜２の範囲の定数で適応動作の調整に用いる。以上のような処理を繰り返すことにより、反響路推定部１９は、次第に擬似反響路ベクトルｈ ₁₁＾（ｋ）を真の反響路１５のインパルス
応答ｈ₁₁（ｋ、ｎ）の時系列を各要素として持つ反響路ベクトルｈ ₁₁（ｋ）、
即ち

と一致させることが可能となり、その結果、式（４）の残留反響信号ｅ１（ｋ）を小さくすることができる。
【０００５】
一般にＮ（≧２）チャネルの再生系とＭ（≧１）チャネルの収音系とで構成される通信会議システムの場合の反響の消去は、図４に示すような構成により行われる。即ち再生側の全Ｎチャネルと収音側の各１チャネルとの間にＮ入力１出力時系列信号を処理するＮチャネル反響消去装置２２₁ 、２２₂ 、…、２２_M をそれぞれ接続した反響消去システム２３として実現される。この場合システム全体でＮ×Ｍ個の反響路１５_nm（１≦ｎ≦Ｎ、１≦ｍ≦Ｍ）が存在する。このシステムの構成単位である再生側の全Ｎチャネルと収音側の各１チャネルとの間に接続されるＮチャネル反響消去装置２２₁ 、２２₂ 、…、２２_M については、図３に示した反響消去装置１４の構成を拡張して、図５に示すように構成される。これは例えば電子情報通信学会誌、８６／１０．Ｖｏｌ．Ｊ６９−ＡＮｏ．１０「多チャネル適応ディジタルフィルタ」に詳しく述べられている。ここで、収音側が第ｍ収音チャネル（１≦ｍ≦Ｍ）に接続されているＮチャネル反響消去装置２２_m を考える。第ｍチャネルのマイクロホン１６_m で収音される反響信号は、各再生チャネルの再生信号がそれぞれの反響路１５_1m〜１５_Nmを経て収音側ですべて加算されることにより得られるために、反響路推定をどの再生チャネルについても、統一的に同じ１つの残留反響信号ｅ_m（ｋ）のみを評価して行うための工夫が必要となる。まず、各再生チャネルの再生信号について、再生信号蓄積・ベクトル生成部（１７₁ ，１７₂ ，…，１７_N ）により、再生信号ベクトル

を生成する。但し、Ｌ₁，Ｌ₂ ，…，Ｌ_Nはタップ数で、各反響路１５₁，１５₂，…，１５_Mの残響時間に対応させて、あらかじめ設定する定数である。これらのベクトルをベクトル結合部２４によって、

と結合する。また、反響路推定部１９_Mにおいても、各再生チャネルと第ｍ収音チャネルとの間のＮ個の反響路を模擬するための、各擬似反響路ベクトルｈ _1m
＾（ｋ），ｈ _2m＾（ｋ），…，ｈ _Nm＾（ｋ）を結合して、

として扱う。擬似反響路結合ベクトルｈ _m＾（ｋ）の更新は、ＮＬＭＳアルゴ
リズムを用いた場合、

のように行われる。擬似反響信号生成部１８_mでは、内積演算
ｙ_m＾（ｋ）＝ｈ _m＾^T（ｋ）Ｘ（ｋ）・・・・・・（13）
により、第ｍ収音チャネルで収音された反響信号ｙ_m（ｋ）に対する擬似反響信号ｙ_m＾（ｋ）を生成する。このように、各チャネル毎のベクトルを結合して１つのベクトルとして扱うことにより、基本的な処理の流れは、図３に示した１チャネル反響消去装置と同様となる。
【０００６】
この多チャネル反響消去方法において、チャネル間の信号の相互相関が強いと適応アルゴリズムは，真値に収束できなくなる。そこで、各チャネルの再生信号に非線形前処理を用い、相互相関を変動させることにより、真値を求められることが知られている。非線形前処理後の信号をｘ ₁ ^-（ｋ），ｘ ₂ ^-（ｋ），…，
ｘ _N ^-（ｋ）とし、非線形処理をｆ₁（），ｆ₂（），…，ｆ_N（）で表すものと
すると、例えば、以下の様に表すことができる。
【０００７】
ｘ ₁ ^-（ｋ）＝ｘ ₁（ｋ）＋α₁ｆ₁（ｘ ₁（ｋ））・・・・・・（14）
ｘ ₂ ^-（ｋ）＝ｘ ₂（ｋ）＋α₂ｆ₂（ｘ ₂（ｋ））・・・・・・（15）
…
ｘ _N ^-（ｋ）＝ｘ _N（ｋ）＋α_Nｆ_N（ｘ _N（ｋ））・・・・・・（15'）
即ち、元の再生信号にα_i （ｉ＝１，２，…，Ｎ）倍した非線形成分を加えることにより、チャネル間の信号の相互相関を変動させることにより適応フィルタによる真値への収束が期待できる。
【０００８】
ここで、各再生信号に加える非線形成分はα_i の大きさに依存するが、この値は再生音質と収束速度のトレードオフとなる。即ち、加算する非線形成分が大きくなると、再生音質が劣化するが収束速度が速くなり、小さくなると再生音質の劣化が小さく押さえられる反面、収束速度が遅くなる。
【０００９】
【発明が解決しようとする課題】
上述した通り、Ｎチャネル再生信号の相互相関を変動させるために加算する非線形成分の大きさは、再生音質の面から小さく抑える必要があり、その結果適応フィルタの収束速度が遅くなる。従来、非線形成分は各チャネルの信号強度に関して所定の一定の比率で各チャネルの信号に加算されていた。
この発明は、再生信号の各チャネル毎のパワーの比が異なる時に、再生信号の小さい方のチャネルの音は大きい方のチャネルの音にマスキングされて聞こえづらいことを利用し、強度が最大のチャネルに対しては所定の非線形成分比率を、強度が最大のチャネル以外のチャネルに対しては強度が最大のチャネルで用いられる非線形成分比率よりも大きい非線形成分比率をそれぞれ用いる構成を採用して、受聴者には再生音質劣化を感じさせずに、適応フィルタの収束速度を改善する反響消去装置、反響消去方法、および反響消去プログラムを提供するものである。
【００１０】
【課題を解決するための手段】
Ｎ個（Ｎは２以上の整数）のチャネルの再生信号をそれぞれＮ個の擬似反響路に通してＮ個の擬似反響信号を求め、Ｎ個の擬似反響信号の和を総合擬似反響信号とし、前記Ｎ個の再生信号を同時に再生して収音された反響信号から前記総合擬似反響信号を差し引くことにより残差信号を求め、前記Ｎ個の再生信号と前記残差信号とを用いてＮ個の疑似反響路を逐次推定する反響消去装置において、前記Ｎチャネル再生信号に各々非線形処理を行って非線形成分信号を生成する非線形前処理部９と、前記Ｎチャネル再生信号の強度をそれぞれ測定する信号測定部７と、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出する非線形成分比率演算部８とを有し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算して、再生信号を再生する反響消去装置を構成した。
【００１１】
そして、先の反響消去装置において、前記非線形成分演算部８は、前記強度が最大のチャネルに対しては、所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては、前記強度が最大のチャネルで用いられる非線形成分比率よりも大きい非線形成分比率をそれぞれ用いる反響消去装置を構成した。また、Ｎ個（Ｎは２以上の整数）のチャネルの再生信号をそれぞれＮ個の擬似反響路に通してＮ個の擬似反響信号を求め、Ｎ個の擬似反響信号の和を総合擬似反響信号とし、前記Ｎ個の再生信号を同時に再生して収音された反響信号から前記総合擬似反響信号を差し引くことにより残差信号を求め、前記Ｎ個の再生信号と前記残差信号とを用いてＮ個の疑似反響路を逐次推定する反響消去方法において、前記Ｎチャネル再生信号に各々非線形処理を行い非線形成分信号を生成し、前記Ｎチャネル再生信号の強度をそれぞれ測定し、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算し、再生信号を再生する反響消去方法を構成した。
【００１２】
更に、Ｎチャネル再生信号に各々非線形処理を行い非線形成分信号を生成し、前記Ｎチャネル再生信号の強度をそれぞれ測定し、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算し、再生信号を再生することを実行させる反響消去プログラムを構成した。
【００１３】
【発明の実施の形態】
この発明の実施の形態を図１の実施例を参照して説明する。
図１においては、説明を簡単にするために、受話２チャネル、収音１チャネルとした場合の例が説明される。各チャネルの周波数領域に変換された再生信号は信号測定部７に入力されて各周波数の信号の強度が測定される。非線形成分演算部８は、その信号強度を入力として、各チャネルに加えることのできる非線形成分を求める。非線形成分演算部８により求められた非線形成分は、非線形前処理部９₁および９₂において再生信号に加えられる。この場合、各チャネル間において、強度が最大のチャネルに対しては、所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては、前記強度が最大のチャネルで用いられる非線形成分比率よりも大きい非線形成分比率をそれぞれ用いて、再生信号の小さいチャネルにより大きい非線形成分を加えることにより、信号の相互相関を大きく変動させる、即ち、小さくすることができる。非線形処理された再生信号はスピーカ２を介して反響路１５に送出されると共に、各Ｎチャネル反響消去装置１０にも送出され、各反響消去装置１０により反響路から収音チャネルに混入する反響信号を低減することができる。
【００１４】
実際に加える非線形成分に関して、聴覚実験により求める方法を説明する。先ず、使用した非線形前処理は以下の式で得られる。非線形成分量は式（16）および（17）のα₁、α₂により制御される。
ｘ₁（ｋ）＝ｕ₁（ｋ）＋α₁（（ｕ₁（ｋ）＋｜ｕ₁（ｋ）｜）／２）・・・（16）
ｘ₂（ｋ）＝ｕ₂（ｋ）＋α₂（（ｕ₂（ｋ）−｜ｕ₂（ｋ）｜）／２）・・・（17）
｜ｕ（ｋ）｜は、ｕ（ｋ）の絶対値を示す。α_i（ｉ＝１、２）は、元の再生信号に加える非線形成分の割合を示す。信号測定部７で信号強度Ｐを求める式を、例えば、以下の通りにして求める。
【００１５】
Ｐx_1(k)＝Σｘ₁ ²（ｋ-i）（ｉ＝０〜適応フィルタタップ数−１）・・・（18）
Ｐx_2(k)＝Σｘ₂ ²（ｋ-i）（ｉ＝０〜適応フィルタタップ数−１）・・・（19）
式（18）、（19）で求められたパワーの差が１５ｄＢあった場合、即ち、
１０ｌｏｇ₁₀（Ｐx_1(k)）−１０ｌｏｇ₁₀（Ｐx_2(k)）＝１５・・・・・・・・（20）
であり、ｘ₁（ｋ）＞ｘ₂（ｋ）であり、α₁＝０.５である場合、α₂＝１.０としても音の劣化はα₂＝０.５としたときと同等であった。同じく、パワー比が１０ｄＢであった場合、ｘ₁（ｋ）＞ｘ₂（ｋ）、α₁＝０.５、α₂＝０.８３で音の劣化が同程度であった。パワー比が０ｄＢ、即ち、等パワーステレオ信号の場合、α₁＝α₂＝０.５と考えられるので、増加させることのできる非線形成分付加量とパワー比が比例していると考えられる。ｘ₁（ｋ）＞ｘ₂（ｋ）の場合の非線形成分加算パラメータは、
α₂＝α₁＋β・・・・・・・・（21）
で与えられ、βは、信号の比がγｄＢとすると、
β＝γ／３０・・・・・・・・（22）
で求められる。
【００１６】
ｘ₁（ｋ）＞ｘ₂（ｋ）で、式（20）を満たしたとき、即ち、信号比が１５ｄＢあったときの従来例の式（16）、式（17）でα₁＝α₂＝０.５とした場合と、同様に、α₁＝０.５、α₂＝１.０としたときの実施例の場合の適応フィルタの推定精度を図２に示す。波線により示される従来例と比較して、実線により示される実施例の方が、速く収束していることが分かる。また、信号のパワー比がないときはγ＝０となり、式（21）でβ＝０となるので、従来例と比較して性能が低下するということはない。
【００１７】
以上においては、受話２チャネル、収音１チャネルとした場合を例として説明したが、以下、チャネル数をＮチャネルに一般化した説明をする。
この発明は、再生信号に加算する非線形成分の再生信号に対する比率を各チャネル毎に設定する。Ｎチャネル再生信号の信号強度を各々測定し、各受話チャネルの内の再生信号強度の最大のチャネルを決定する。再生信号強度最大チャネル以外のチャネルについて、強度最大チャネルの再生信号強度との間の強度差或いは強度比を求める。ここで、各再生信号をｘ_i（ｋ）（ｉ＝１〜Ｎ，Ｎ：チャネル数）とし、加算する非線形成分をｆ₁（ｘ_i（ｋ））、非線形成分比率をα_i（ｋ）とすると、非線形処理されスピーカから出力される信号ｕ_i（ｋ）は、
ｕ_i（ｋ）＝ｘ_i（ｋ）＋α_i（ｋ）ｆ₁（ｘ_i（ｋ））
となる。ここで、非線形処理には、例えば、半波整流付加を用いる。半波整流付加は、以下の様に与えられる。
ｆ₁（ｘ_i（ｋ））＝ｘ_i（ｋ）＋０.５（ｘ_i（ｋ）＋（−１）＾_i｜ｘ_i（ｋ）｜）
ここで、強度最大チャネルがチャネルｊであると判定されたとき、非線形成分比率α_i（ｋ）は以下のように得られる。
【００１８】

ここで、ｋは現在のサンプル、ｔはパワーを求める時定数、αは強度最大チャネルに使用される非線形成分比率、ｂは信号強度差もしくは信号強度比から非線形成分比率を決定する定数である。α、ｂは再生音を受聴し、その再生品質の優劣を評価する主観評価を用いて決定されたものを用いてもよい。
【００１９】
例えば、Ｎ＝２チャネルの場合は、α＝０.５、ｂ＝３０となる。一般に、非線形成分比率を多くするほど再生音の音質は劣化するが、時不変で各チャネルに対して同じ値の非線形成分比率を用いる従来法において、α₁（ｋ）＝α₂（ｋ）＝０.５（固定）を用いた場合と同等な音質が得られる。例えば、チャネル間のレベル差が１５ｄＢであったとすると、再生信号と非線形成分の比は、従来例は５０％、実施例では約７０％となり、実施例の方が加算する非線形成分の量が多い。
【００２０】
Ｎ＝２、チャネル間レベル差が１５ｄＢである場合の適応フィルタの推定精度の計算機シミュレーション結果を図２に示す。反響路を推定する適応アルゴリズムには、学習同定法（ＮＬＭＳ）、適応フィルタ長は５１２タップ、サンプリング周波数は１６ｋＨｚ、再生信号は白色雑音に、実測したインパルス応答を畳み込みステレオ信号を生成したものを用いている。図２は、３２、０００サンプルのときに、反響路を変動させ、その時点からの収束速度を示している。従来例（破線）は、推定精度を−１０ｄＢまで推定するのに、２６８、０００サンプル必要であるのと比較して、実施例（実線）は同様に推定精度を−１０ｄＢまで推定するのに、１４８、０００サンプル必要となり、同じ推定精度（−１０ｄＢ）まで推定するのに、約１．８倍高速化されていることが分かる。
【００２１】
一般的な通信会議の強度（受聴位置で約７０ｄＢspl）付近（＋−６ｄＢ）においては、絶対強度を変化（受聴位置で６４ｄＢspl〜７６ｄＢspl程度）させても非線形成分比率を求める定数αとｂに、上述した値α＝０．５、ｂ＝３０を用いても、以上と同様に再生音を受聴させた場合、実施例と従来例とで再生音の音質は同等となる。即ち、受聴位置において絶対強度が変化しても、実施例の効果は同様に得られると言える。
なお、以上の反響消去装置を電子計算機を主要な構成部材として構成してもよい。そして、この発明を、ＣＤその他の記憶媒体からダウンロードし或いは通信回線を介してダウンロードしたプログラムをこの電子計算機にインストールして実施することができる。
【００２２】
【発明の効果】
以上の通りであって、この発明によれば、各チャネルの再生信号の強度を測定し、従来の如く一定の比率ではなしに信号強度差或いは信号強度比に応じた量の非線形成分を再生信号に加算することにより、受聴者に対して再生音質の劣化を感じさせることなく、反響消去装置の収束速度を改善することができる。
【図面の簡単な説明】
【図１】実施例を説明する図。
【図２】適応フィルタの推定精度を説明する図。
【図３】従来例を説明する図。
【図４】他の従来例を説明する図。
【図５】反響消去装置の詳細を説明する図。
【符号の説明】
７信号測定部
８非線形成分演算部
９非線形前処理部

Claims

Ｎ個（Ｎは２以上の整数）のチャネルの再生信号をそれぞれＮ個の擬似反響路に通してＮ個の擬似反響信号を求め、Ｎ個の擬似反響信号の和を総合擬似反響信号とし、前記Ｎ個の再生信号を同時に再生して収音された反響信号から前記総合擬似反響信号を差し引くことにより残差信号を求め、前記Ｎ個の再生信号と前記残差信号とを用いてＮ個の疑似反響路を逐次推定する反響消去装置において、
前記Ｎチャネル再生信号に各々非線形処理を行って非線形成分信号を生成する非線形前処理部と、前記Ｎチャネル再生信号の強度をそれぞれ測定する信号測定部と、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出する非線形成分比率演算部とを有し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算して、再生信号を再生することを特徴とする反響消去装置。
請求項１に記載される反響消去装置において、
前記非線形成分演算部８は、前記強度が最大のチャネルに対しては、所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては、前記強度が最大のチャネルで用いられる非線形成分比率よりも大きい非線形成分比率をそれぞれ用いることを特徴とする反響消去装置。
Ｎ個（Ｎは２以上の整数）のチャネルの再生信号をそれぞれＮ個の擬似反響路に通してＮ個の擬似反響信号を求め、Ｎ個の擬似反響信号の和を総合擬似反響信号とし、前記Ｎ個の再生信号を同時に再生して収音された反響信号から前記総合擬似反響信号を差し引くことにより残差信号を求め、前記Ｎ個の再生信号と前記残差信号とを用いてＮ個の疑似反響路を逐次推定する反響消去方法において、
前記Ｎチャネル再生信号に各々非線形処理を行って非線形成分信号を生成し、
前記Ｎチャネル再生信号の強度をそれぞれ測定し、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算し、再生信号を再生することを特徴とする反響消去方法。
Ｎ個（Ｎは２以上の整数）のチャネルの再生信号をそれぞれＮ個の擬似反響路に通してＮ個の擬似反響信号を求め、Ｎ個の擬似反響信号の和を総合擬似反響信号とし、前記Ｎ個の再生信号を同時に再生して収音された反響信号から前記総合擬似反響信号を差し引くことにより残差信号を求め、前記Ｎ個の再生信号と前記残差信号とを用いてＮ個の疑似反響路を逐次推定する反響消去プログラムであって、
前記Ｎチャネル再生信号に各々非線形処理を行って非線形成分信号を生成し、
前記Ｎチャネル再生信号の強度をそれぞれ測定し、前記強度が最大のチャネルに対しては所定の非線形成分比率を、前記強度が最大のチャネル以外のチャネルに対しては前記強度が最大のチャネルの信号の強度との強度差または強度比に基づく所定の非線形成分比率を算出し、前記非線形成分信号に各々対応するチャネルの非線形成分比率で前記Ｎチャネル再生信号に加算し、再生信号を再生することを実行させる反響消去プログラム。