"kは時間を表し、fは周波数を表し、jはフレーム番号を表し、iは適応フィルタのタップ番号を表し、mはスピーカの番号を表し、nはマイクロホンの番号を表し、SetNfは周波数fにおいて空間サンプリング間隔に基づいて選択されるマイクロホンの番号の集合を表すとして、受話信号ベクトル変換部が、M1≦m≦M個のスピーカから出力されるMチャネルの受話信号xmkを、チャネルmごとに、周波数領域に変換して、受話信号Xmf,jを生成する受話信号ベクトル変換ステップと、エコーレプリカ生成部が、前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnについて、タップ数がI1≦i≦I、Iは1以上のフィルタ係数Hm,n,if,jを用いて、空間的に間引かれたエコーレプリカY^SetNff,jを生成するエコーレプリカ生成ステップと、エコーレプリカ空間補間部が、前記空間的に間引かれたエコーレプリカY^SetNff,jから、周波数fごとに、空間補間を行い、エコーレプリカY^nf,jを生成するエコーレプリカ空間補間ステップと、逆FFT部が、前記エコーレプリカY^nf,jを、時間領域に変換して、エコーレプリカy^njを生成する逆FFTステップと、エコー消去部が、直線上に等間隔に配置されたN3≦n≦N個のマイクロホンから収音されるNチャネルの送話信号ynjと前記エコーレプリカy^njから、誤差信号enkを生成するエコー消去ステップと、FFT部が、前記誤差信号enkを、周波数領域に変換して、誤差信号Enf,jを生成するFFTステップと、前記エコーレプリカ生成部が、前記誤差信号Enf,jと前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnについて、修正量dHm,n,if,jを求め、当該修正量dHm,n,if,jを用いて、前記フィルタ係数Hm,n,if,jを更新するフィルタ係数更新ステップと、を含むことを特徴とするマルチチャネルエコー消去方法。kは時間を表し、fは周波数を表し、jはフレーム番号を表し、iは適応フィルタのタップ番号を表し、mはスピーカの番号を表し、SetMfは周波数fにおいて空間サンプリング間隔に基づいて選択されるスピーカの番号の集合を表し、nはマイクロホンの番号を表し、SetNfは周波数fにおいて空間サンプリング間隔に基づいて選択されるマイクロホンの番号の集合を表すとして、受話信号ベクトル変換部が、直線上に等間隔に配置されたM3≦m≦M個のスピーカから出力されるMチャネルの受話信号xmkを、チャネルmごとに、周波数領域に変換して、受話信号Xmf,jを生成する受話信号ベクトル変換ステップと、エコーレプリカ生成部が、前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnと前記SetMfに含まれるスピーカmの組み合わせについて、タップ数がI1≦i≦I、Iは1以上のフィルタ係数Hm,n,if,jを用いて、空間的に間引かれたエコーレプリカY^SetNff,jを生成するエコーレプリカ生成ステップと、エコーレプリカ空間補間部が、前記空間的に間引かれたエコーレプリカY^SetNff,jから、周波数fごとに、空間補間を行い、エコーレプリカY^nf,jを生成するエコーレプリカ空間補間ステップと、逆FFT部が、前記エコーレプリカY^nf,jを、時間領域に変換して、エコーレプリカy^njを生成する逆FFTステップと、エコー消去部が、直線上に等間隔に配置されたN3≦n≦N個のマイクロホンから収音されるNチャネルの送話信号ynjと前記エコーレプリカy^njから、誤差信号enkを生成するエコー消去ステップと、FFT部が、前記誤差信号enkを、周波数領域に変換して、誤差信号Enf,jを生成するFFTステップと、前記エコーレプリカ生成部が、前記誤差信号Enf,jと前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnと前記SetMfに含まれるスピーカmの組み合わせについて、修正量dHm,n,if,jを求め、当該修正量dHm,n,if,jを用いて、前記フィルタ係数Hm,n,if,jを更新するフィルタ係数更新ステップと、を含むことを特徴とするマルチチャネルエコー消去方法。請求項1または2に記載のマルチチャネルエコー消去方法であって、Sはマイクロホン間隔を表し、λfは周波数fの波長を表すとして、前記SetNfは、S×Q≦λf/2を満たす最大の自然数Qを求め、当該Q以下である間引き間隔Qfを決定し、当該間引き間隔Qfに基づいて選択されることを特徴とするマルチチャネルエコー消去方法。請求項1から3のいずれかに記載のマルチチャネルエコー消去方法であって、Sはマイクロホン間隔を表し、λfは周波数fの波長を表すとして、前記エコーレプリカ空間補間ステップは、周波数fごとに、S×Q≦λf/2を満たす最大の自然数Qを求め、当該Q以下である間引き間隔Qfを決定し、のように、空間補間を行うことを特徴とするマルチチャネルエコー消去方法。請求項1から4のいずれかに記載のマルチチャネルエコー消去方法であって、前記フィルタ係数Hm,n,if,jは、周波数fごとに可変であり、所定の周波数を超える周波数に対応するフィルタ係数のタップ数が、所定の周波数以下の周波数に対応するフィルタ係数のタップ数よりも少なく設定されていることを特徴とするマルチチャネルエコー消去方法。請求項1から5のいずれかに記載のマルチチャネルエコー消去方法であって、帯域分割音声スイッチ部が、前記受話信号xmkを、チャネルmごとに、高域成分を減衰する受話側高域減衰ステップと、前記帯域分割音声スイッチ部が、前記誤差信号enkを、チャネルnごとに、高域成分を減衰する送話側高域減衰ステップと、前記帯域分割音声スイッチ部が、前記受話信号xmkの低域成分と前記誤差信号enkの低域成分を用いて、送受話判定を行う送受話判定ステップを、さらに有し、前記送受話判定ステップにおいて、受話状態と判定した場合には、前記受話側高域減衰ステップを実行し、送話状態と判定した場合には、前記送話側高域減衰ステップを実行することを特徴とするマルチチャネルエコー消去方法。kは時間を表し、fは周波数を表し、jはフレーム番号を表し、iは適応フィルタのタップ番号を表し、mはスピーカの番号を表し、nはマイクロホンの番号を表し、SetNfは周波数fにおいて空間サンプリング間隔に基づいて選択されるマイクロホンの番号の集合を表すとして、M1≦m≦M個のスピーカと、直線上に等間隔に配置されたN3≦n≦N個のマイクロホンと、前記スピーカから出力されるMチャネルの受話信号xmkを、チャネルmごとに、周波数領域に変換して、受話信号Xmf,jを生成する受話信号ベクトル変換部と、前記受話信号Xmf,jが入力されると、当該受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnについて、タップ数がI1≦i≦I、Iは1以上のフィルタ係数Hm,n,if,jを用いて、空間的に間引かれたエコーレプリカY^SetNff,jを生成し、誤差信号Enf,jが入力されると、当該誤差信号Enf,jと前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnについて、修正量dHm,n,if,jを求め、当該修正量dHm,n,if,jを用いて、前記フィルタ係数Hm,n,if,jを更新するエコーレプリカ生成部と、前記空間的に間引かれたエコーレプリカY^SetNff,jから、周波数fごとに、空間補間を行い、エコーレプリカY^nf,jを生成するエコーレプリカ空間補間部と、前記エコーレプリカY^nf,jを、時間領域に変換して、エコーレプリカy^njを生成する逆FFT部と、前記マイクロホンから収音されるNチャネルの送話信号ynjと前記エコーレプリカy^njから、誤差信号enkを生成するエコー消去部と、前記誤差信号enkを、周波数領域に変換して、前記誤差信号Enf,jを生成するFFT部と、を備えることを特徴とするマルチチャネルエコー消去装置。kは時間を表し、fは周波数を表し、jはフレーム番号を表し、iは適応フィルタのタップ番号を表し、mはスピーカの番号を表し、SetMfは周波数fにおいて空間サンプリング間隔に基づいて選択されるスピーカの番号の集合を表し、nはマイクロホンの番号を表し、SetNfは周波数fにおいて空間サンプリング間隔に基づいて選択されるマイクロホンの番号の集合を表すとして、直線上に等間隔に配置されたM3≦m≦M個のスピーカと、直線上に等間隔に配置されたN3≦n≦N個のマイクロホンと、前記スピーカから出力されるMチャネルの受話信号xmkを、チャネルmごとに、周波数領域に変換して、受話信号Xmf,jを生成する受話信号ベクトル変換部と、前記受話信号Xmf,jが入力されると、当該受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnと前記SetMfに含まれるスピーカmの組み合わせについて、タップ数がI1≦i≦I、Iは1以上のフィルタ係数Hm,n,if,jを用いて、空間的に間引かれたエコーレプリカY^SetNff,jを生成し、誤差信号Enf,jが入力されると、当該誤差信号Enf,jと前記受話信号Xmf,jから、周波数fごとに、前記SetNfに含まれるマイクロホンnと前記SetMfに含まれるスピーカmの組み合わせについて、修正量dHm,n,if,jを求め、当該修正量dHm,n,if,jを用いて、前記フィルタ係数Hm,n,if,jを更新するエコーレプリカ生成部と、前記空間的に間引かれたエコーレプリカY^SetNff,jから、周波数fごとに、空間補間を行い、エコーレプリカY^nf,jを生成するエコーレプリカ空間補間部と、前記エコーレプリカY^nf,jを、時間領域に変換して、エコーレプリカy^njを生成する逆FFT部と、前記マイクロホンから収音されるNチャネルの送話信号ynjと前記エコーレプリカy^njから、誤差信号enkを生成するエコー消去部と、前記誤差信号enkを、周波数領域に変換して、前記誤差信号Enf,jを生成するFFT部と、を備えることを特徴とするマルチチャネルエコー消去装置。請求項1から6のいずれかに記載されたマルチチャネルエコー消去方法の各ステップをコンピュータに実行させるためのプログラム。しかしながら、非特許文献1に記載されたWFSは、ある地点での音波面を取得し別の地点で再合成するために、数十以上のマイクロホンと数十以上のスピーカを必要とする。そのため、WFSを双方向テレビ会議に導入しようとする場合、スピーカマイクロホンの間の音響パスはチャネル数の二乗になるため、エコーキャンセラが推定する音響パス数はチャネル数の二乗で急激に増大する。また、エコーキャンセラの演算量を下げるために、非特許文献2に記載された周波数領域でエコー計算とフィルタ係数更新を行う方法を用いても、エコーレプリカ生成とフィルタ係数更新の処理は、入力チャネル数と出力チャネル数の積、すなわちM×Nで増大してしまう。本発明はこのような点に鑑みてなされたものであり、マルチチャネル拡声通話系において音響エコーを消去するマルチチャネルエコー消去処理における演算量を低減することができるマルチチャネルエコー消去装置を提供することを目的とする。上記の課題を解決するために、本発明のマルチチャネルエコー消去装置は、M1≦m≦M個のスピーカと直線上に等間隔に配置されたN3≦n≦N個のマイクロホンと受話信号ベクトル変換部とエコーレプリカ生成部とエコーレプリカ空間補間部と逆FFT部とエコー消去部とFFT部を備える。kは時間を表し、fは周波数を表し、jはフレーム番号を表し、iは適応フィルタのタップ番号を表し、mはスピーカの番号を表し、nはマイクロホンの番号を表し、SetNfは周波数fにおいて空間サンプリング間隔に基づいて選択されるマイクロホンの番号の集合を表すとする。受話信号ベクトル変換部は、スピーカから出力されるMチャネルの受話信号xmkを、チャネルmごとに、周波数領域に変換して、受話信号Xmf,jを生成する。エコーレプリカ生成部は、受話信号Xmf,jが入力されると、当該受話信号Xmf,jから、周波数fごとに、SetNfに含まれるマイクロホンnについて、タップ数がI1≦i≦I、Iは1以上のフィルタ係数Hm,n,if,jを用いて、空間的に間引かれたエコーレプリカY^SetNff,jを生成する。エコーレプリカ空間補間部は、空間的に間引かれたエコーレプリカY^SetNff,jから、周波数fごとに、空間補間を行い、エコーレプリカY^nf,jを生成する。逆FFT部は、エコーレプリカY^nf,jを、時間領域に変換して、エコーレプリカy^njを生成する。エコー消去部は、マイクロホンから収音されるNチャネルの送話信号ynjとエコーレプリカy^njから、誤差信号enkを生成する。FFT部は、誤差信号enkを、周波数領域に変換して、誤差信号Enf,jを生成する。エコーレプリカ生成部は、誤差信号Enf,jが入力されると、当該誤差信号Enf,jと受話信号Xmf,jから、周波数fごとに、SetNfに含まれるマイクロホンnについて、修正量dHm,n,if,jを求め、当該修正量dHm,n,if,jを用いて、フィルタ係数Hm,n,if,jを更新する。本発明のマルチチャネルエコー消去装置によれば、周波数ごとにエコー信号の空間方向形状に着目し、低い周波数では空間補間を用いてエコーレプリカを生成することで、エコーレプリカの生成とフィルタ係数の更新における演算回数を低減し、マルチチャネル拡声通話系におけるエコー消去処理全体での演算量を低減することができる。以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。本発明の原理実施例の説明に先立ち、本発明のマルチチャネルエコー消去方法の原理を説明する。本発明は、周波数ごとのエコー信号の空間方向の形状に着目し、低い周波数では空間補間を用いてエコーレプリカを生成することで、エコーキャンセラ演算量の急激な増大を抑えることができる。マイクロホンが等間隔に直線状に並んでいる場合、各周波数でのエコー信号の空間形状は、周波数が低いほど滑らかになる。それは、音波の波長λと周波数fの間に、波長λ×周波数f=音速c一定の関係があるために、周波数fが低いほど波長λは長くなるためである。さらにナイキスト定理によれば、低い周波数ではサンプリング間隔をマイクロホン間隔よりもずっと大きくすることができる。これは、低い周波数では空間サンプリングの間隔を広げられること、すなわちサンプリングを間引けること、を意味する。ここではWFS収音再生系の一例として、直線状スピーカアレーと直線状マイクロホンアレーをサンプリング周波数16kHzで動作させることを考える。ナイキスト定理より、マイクロホン間隔すなわち空間波のサンプリング間隔はλ/2以下である必要がある。帯域上限8kHzの波をサンプリングするには、音速を330[m/s]として、マイクロホン間隔を2cmにする必要がある。しかし周波数が低くなるほど波形は滑らかになる。図4にこの状況を示す。33本のマイクロホンで収音したエコー信号1フレーム分を高速フーリエ変換により周波数分解して得られた、4kHz,2kHz,1kHz,500Hz,250Hz,125Hzのエコー信号の空間形状をそれぞれプロットしている。なお、実線は実数成分を表し、点線は虚数成分を表している。滑らかな波形に対して空間サンプリングを間引きすることができる。すなわち低い周波数については、全マイクロホンではなく一部のマイクロホンでエコー消去処理を行い、残りのマイクロホンについては適応フィルタを用いずに一部のマイクロホンのエコー消去済み信号から、空間補間によって全マイクロホンのエコー消去済み信号を求めることができる。図5に、マイクロホン数9、周波数帯域分割数8の場合について、周波数fごとの送話信号の間引きパターンSetNfを示す。横軸はマイク位置であり、縦軸は周波数である。fの値が大きいほど周波数が高いことを表している。f=3,4では、空間サンプリングが1/2に間引かれている。f=1,2では空間サンプリングが1/4に間引かれている。送話信号の間引きパターンSetNfは、周波数fの波長λf=c音速/fから求められる。マイクロホンが等間隔に並んでいる場合、SNはマイクロホン間隔を表すとして、SN×Q≦λf/2を満たす最大の自然数Qを求め、間引き間隔QfをQ以下に設定する。例えば、マイクロホン数N=9、間引き間隔Qf=4の場合、SetNf={1,5,9}と設定することができる。また、間引き間隔Qf=2の場合、SetNf={1,3,5,7,9}と設定することができる。上記の原理に基づいて、低い周波数でエコーレプリカの生成処理とフィルタ係数の更新処理を間引くことで演算量を下げることが可能になる。図6、図7を参照して、本発明の実施例1に係るマルチチャネルエコー消去装置20の動作を詳細に説明する。図6は本発明の実施例1に係るマルチチャネルエコー消去装置20の構成を示すブロック図である。図7は本発明の実施例1に係るマルチチャネルエコー消去装置20の動作を示すフローチャートである。本実施例では、周波数領域適応アルゴリズムとして、EMoulinesOAAmraneandYGrenierTheGeneralizedMultidelayAdaptiveFilterStructureandConvergenceAnalysisIEEETransonSPvol43no1に記載された方法を用いた場合について説明する。以下では、適応フィルタは時間方向I個Iは1以上に分割されているものとして説明する。以下、実際に行われる手続きの順に説明してゆく。本実施例のマルチチャネルエコー消去装置20は、M1≦m≦M個のスピーカ21~2Mと、N3≦n≦N個のマイクロホン31~3Nと、受話信号ベクトル変換部110と、N個のエコーレプリカ生成部2201~220Nと、エコーレプリカ空間補間部300と、逆FFT部400と、エコー消去部510と、FFT部600を備える。受話信号ベクトル変換部110は、Mチャネルの受話信号xmkをLサンプルごとにブロック化し、1フレーム=2Lサンプルとして、1フレーム分を高速フーリエ変換により周波数領域に変換し、上記式1のように受話信号Xmf,jを生成するS110。エコーレプリカ生成部220nは、周波数f≦L1について、周波数fごとに、SetNfに含まれるnについて、式8のように、受話信号Xmf,jとフィルタ係数Hm,n,if,jを掛けることで、受話信号Xmf,jをフィルタ処理し、Mチャネル分加算する。これにより空間的に間引きされたエコーレプリカY^SetNff,jを求めるS220。エコーレプリカ空間補間部300は、周波数f≦L1について、周波数fごとに、エコーレプリカ生成部210nの生成する空間的に間引きされたエコーレプリカY^SetNff,jを用いて、空間補間によりエコーレプリカY^nf,jを求めるS300。より具体的には、間引きによりスキップされたチャネルn、すなわちSetNfに含まれないnについて、SetNfに含まれるチャネルのエコー消去済み信号を空間補間することでエコーレプリカY^nf,jを求める。マイクロホンが等間隔で直線状に配置されているときは、空間的に間引きされたエコーレプリカY^SetNff,jとsinc関数を使って、式9のように空間補間が可能である。なお、実際の計算では、sinc関数を有限長で打ち切る必要がある。sinc関数を打ち切る範囲としては、間引き間隔Qfに基づいて、式10もしくは式11のように設定することができる。この場合には、sinc関数と乗算するエコーレプリカは4点もしくは6点で済む。続いて、エコーレプリカ空間補間部300は、周波数f>L1について、周波数fごとに、実数信号の高速フーリエ変換結果に関する対称性を用いて、式12のようにエコーレプリカY^nf,jを求めるS300。ここで、conjは複素共役を取ることを表している。逆FFT部400は、エコーレプリカY^nf,jを逆高速フーリエ変換により時間領域に変換し、上記式3のようにエコーレプリカy^njを求めるS400。エコー消去部510は、時間領域でN個のマイクロホン31~3Nから収音されるNチャネルの送話信号ynjとエコーレプリカy^njから誤差信号enjを求めるS510。FFT部600は、誤差信号enjを高速フーリエ変換により周波数領域に変換し、上記式4のように誤差信号Enf,jを求めるS600。受信信号xmkをフレーム化するときフレーム分割数DをD≧2に設定した場合には、フレーム番号jで求めた誤差信号Enk,jと、一つ前のフレーム番号j-1で求めた誤差信号Enk,j-1を窓かけ処理を経て合成して出力する。図9はD=2の場合の合成処理を示す図である。具体的には、jフレーム目で求めた誤差信号enk0t,jを式13のように表し、WHは長さ2L/Dのハニング窓を表すとして、合成後の誤差信号en’k0tは、式14のように表すことができる。エコーレプリカ生成部220nは、SetNfに含まれるnについてのみ、誤差信号Enf,jと受話信号Xmf,jから、式15のように、フィルタ係数の修正量dHm,n,if,jを求める。次に、SetNfに含まれるnについてのみ、各チャネルのフィルタ係数Hm,n,if,jを、式16のように更新するS710。ここで、pf,jは、周波数成分ごとに、Nチャネル分の送話信号パワーの総和を、式17のように計算することで求めたものであり、修正量dHm,n,if,jを補正している。ここで、μは0~1の値をとるステップサイズであり、δは分母が0になることを防止するための微小な正定数であり、βは0~1の値をとるパワー計算で短時間平均をとるための平滑化定数である。上記から明らかなように、本実施例では、SetNfに含まれないnについては、エコーレプリカy^njも修正量dHm,n,if,jもフィルタ係数Hm,n,if,jも計算されないため、マルチチャネルエコー消去処理全体での演算量を低減することができる。実施例1の実験結果本実施例の効果を確認するために、シミュレーションを行った。残響時間200msの部屋で、直線状スピーカアレー33素子、間隔2cmと直線状マイクロホンアレー33素子、間隔2cmを2m離して平行に配置し、スピーカマイク間の全エコー経路インパルス応答をシミュレータにより生成した。サンプリング周波数は16kHzに設定した。空間補間でエコーレプリカを求める際に、間引き間隔Qfに基づいて上記式10を使用した。受話信号は、別途ピンクノイズを音源として、33マイクロホンによる収音を模擬して生成した。図10に本実施例の構成によるエコー消去処理結果を示す。図10では33チャネル中の第1,5,11,17,19チャネルについて、送話信号と誤差信号のレベルを、受話信号を実線で、誤差信号を破線でプロットしている。実線と破線の差だけエコーが消去されていることを表している。いずれもエコーを良好に消去していることが分かる。変形例本実施例では、周波数領域のフィルタ係数Hm,n,if,jが、どの周波数fでもIタップを持つ場合を説明している。本実施例の変形例として、フィルタ係数の高域部分でタップ数をIタップより減らすことが考えられる。例えば、人の声の成分が相対的に小さい4kHzを超える周波数ではタップ数を1とすることなどが考えられる。これにより、全体の演算量をさらに下げることが可能になる。演算量をさらに下げるために、適応フィルタに入力する受話信号についても低い周波数で間引きを行ってもよい。周波数fごとの受話信号の間引きパターンSetMfは、SetNfと同様に、周波数fの波長λf=c/fから求められる。スピーカが等間隔に並んでいる場合、SMはスピーカ間隔を表すとして、SM×Q≦λf/2を満たす最大の自然数Qを求め、間引き間隔QMfをQ以下に設定する。図11、図12を参照して、本発明の実施例2に係るマルチチャネルエコー消去装置30の動作を詳細に説明する。図11は本発明の実施例2に係るマルチチャネルエコー消去装置30の構成を示すブロック図である。図12は本発明の実施例2に係るマルチチャネルエコー消去装置30の動作を示すフローチャートである。実施例2のマルチチャネルエコー消去装置30は、実施例1のマルチチャネルエコー消去装置20と比較して、エコーレプリカ生成部220nの替わりに、エコーレプリカ生成部230nを備える点が相違する。エコーレプリカ生成部230nは、周波数f≦L1について、周波数fごとに、mがSetMfに含まれる受話信号Xmf,jを抽出するS231。続いて、エコーレプリカ生成部230nは、周波数f≦L1について、周波数fごとに、SetNfに含まれるnについて、式18のように、受話信号Xmf,jとフィルタ係数Hm,n,if,jを掛けることで、受話信号Xmf,jをフィルタ処理し、Mチャネル分加算する。これにより空間的に間引きされたエコーレプリカY^SetNff,jを求めるS232。また、エコーレプリカ生成部230nは、mがSetMfに含まれ、かつnがSetNfに含まれる場合にのみ、誤差信号Enf,jと受話信号Xmf,jから、上記式15のように、フィルタ係数の修正量dHm,n,if,jを求める。次に、mがSetMfに含まれ、かつnがSetNfに含まれる場合にのみ、各チャネルのフィルタ係数Hm,n,if,jを、上記式16のように更新するS720。ただし、修正量dHm,n,if,jを補正するpf,jは、式19のように計算して求める。ここで、μは0~1の値をとるステップサイズであり、δは分母が0になることを防止するための微小な正定数であり、βは0~1の値をとるパワー計算で短時間平均をとるための平滑化定数である。上記から明らかなように、本実施例では、m∈SetMfかつn∈SetNfのm,nについてのみ、エコーレプリカY^nf,jと修正量dHm,n,if,jが計算され、フィルタ係数Hm,n,if,jが更新される。このため、上記以外のm,nについては、エコーレプリカY^nf,jも修正量dHm,n,if,jもフィルタ係数Hm,n,if,jも計算されないため、演算量をさらに低減することができる。なお、本実施例では、受話側の信号が間引かれているため、フィルタ係数と特定スピーカ特定マイク間のエコー経路特性は1対1に対応しなくなる。実施例2の実験結果本実施例の効果を確認するために、実施例1と同一の設定でシミュレーションを行った。図13に本実施例の構成によるエコー消去処理結果を示す。図13では33チャネル中の第1、5、11、17、19チャネルについて、送話信号と誤差信号のレベルを、受話信号を実線で、誤差信号を破線でプロットしている。実線と破線の差だけエコーが消去されていることを表している。いずれもエコーを良好に消去していることが分かる。図14、図15、図16を参照して、本発明の実施例3に係るマルチチャネルエコー消去装置40の動作を詳細に説明する。図14は本発明の実施例3に係るマルチチャネルエコー消去装置40の構成を示すブロック図である。図16は本発明の実施例3に係るマルチチャネルエコー消去装置40の動作を示すフローチャートである。実施例3のマルチチャネルエコー消去装置40は、実施例1のマルチチャネルエコー消去装置20と比較して、さらに帯域分割音声スイッチ部800を備える点が相違する。図15は帯域分割音声スイッチ部800の構成を示すブロック図である。受話信号および送話信号を帯域分割音声スイッチ部800で低域成分と高域成分に分け、高域成分のエコーを送受話状態に応じて減衰させることで制御し、低域成分のエコーを適応フィルタで消去して制御する。帯域分割音声スイッチ部800は、送受話判定部810と、M個の受話側高域減衰部8201~820Mと、N個の受話側高域減衰部8301~830Nを備える。受話側高域減衰部820mは、ハイパスフィルタ以下、HPFと略す。8210と、ローパスフィルタ以下、LPFと略す。8220と、信号減衰器8230と、信号加算器8240を備える。送話側高域減衰部830nは、HPF8310と、LPF8320と、信号減衰器8330と、信号加算器8340を備える。受話側高域減衰部820mは、受話信号にHPF8210とLPF8220を適用し、低域信号と高域信号に分割する。信号減衰器8230は、HPF8210の出力する高域信号を指定分減衰する。信号加算器8240は、信号減衰器8230の出力する信号とLPF8220の出力する信号を加算する。送話側高域減衰部820nは、送話信号にHPF8310とLPF8320を適用し、低域信号と高域信号に分割する。信号減衰器8330は、HPF8310の出力する高域信号を指定分減衰する。信号加算器8340は、信号減衰器8330の出力する信号と、LPF8320の出力する信号を加算する。送受話判定部810は、LPF8220の出力するMチャネルの信号とLPF8320の出力するNチャネルの信号を用いて、送受話判定を行うS810。受話状態と判定したときには受話側高域信号のみを減衰させるS820。送話状態と判定したときには送話側高域信号のみを減衰させるS830。信号減衰器8220および信号減衰器8320に指定する減衰量は、残留エコーが気にならないレベルになるように、3~40dBの範囲で設定する。本実施例によると、エコーレプリカ生成部220mが高域成分に対応するフィルタ係数Hm,nf,jをエコーレプリカ生成に使用しないため、全体での演算量を低減することができる。ただし、フレーム間で異音が生じないことを保証するために受話信号ベクトル変換部110が受話信号をフレーム化する際に、フレーム分割数DをD≧2に設定する必要がある。<プログラム、記録媒体>上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。また、このプログラムの流通は、例えば、そのプログラムを記録したDVD、CD-ROM等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるASPApplicationServiceProvider型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるものコンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等を含むものとする。また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。"