(57)【特許請求の範囲】
【請求項１】
１あるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理装置であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段と、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段と、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３の手段と、
を具備することを特徴とする信号処理装置。
【請求項２】
前記第１の手段は、それぞれ異なる周波数帯域の信号を通過する複数の帯域通過フィルタで構成されることを特徴とする請求項１に記載の信号処理装置。
【請求項３】
前記信号源は音源であることを特徴とする請求項１に記載の信号処理装置。
【請求項４】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理装置であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段と、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段と、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３の手段と、
を具備することを特徴とする信号処理装置。
【請求項５】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理装置であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段と、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段と、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３の手段と、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４の手段と、
を具備することを特徴とする信号処理装置。
【請求項６】
前記第１の手段は、それぞれ異なる周波数帯域の信号を通過する複数の帯域通過フィルタで構成されることを特徴とする請求項５に記載の信号処理装置。
【請求項７】
前記信号源は音源であることを特徴とする請求項５に記載の信号処理装置。
【請求項８】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理装置であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段と、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段と、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３の手段と、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４の手段と、
を具備することを特徴とする信号処理装置。
【請求項９】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理方法であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１のステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２のステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３のステップと、
を具備することを特徴とする信号処理方法。
【請求項１０】
前記信号源は音源であることを特徴とする請求項９に記載の信号処理方法。
【請求項１１】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理方法であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１のステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２のステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３のステップと、
を具備することを特徴とする信号処理方法。
【請求項１２】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理方法であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１のステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２のステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３のステップと、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４のステップと、
を具備することを特徴とする信号処理方法。
【請求項１３】
前記信号源は音源であることを特徴とする請求項１２に記載の信号処理方法。
【請求項１４】
１つあるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理方法であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１のステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２のステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３のステップと、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４のステップと、
を具備することを特徴とする信号処理方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、音源などの信号源が存在する方向や位置を推定する信号源方向又は位置推定方法及び装置に係り、特に、複数の受信装置で受信される各受信信号間の相互相関関数に基づいて１つあるいは複数の信号源の方向や位置を推定する信号源方向又は位置推定方法及び装置に関する。
【０００２】
また、本発明は、複数の受信装置で受信される各受信信号間の相互相関関数に基づいて１つあるいは複数の信号源の信号を強調する信号強調方法及び装置に関する。
【０００３】
【従来の技術】
雑音の多い環境や、複数の信号源が同時に信号を発生しているような環境下で、信号源の方向や位置を推定する技術は、信号源からの信号をより良く受信するよう受信装置を適応させるシステムや、信号源の方向にビデオ・カメラを自動的に向けて監視を行うシステム、発言者の方向にビデオ・カメラを自動的に向けて発言者の映像を送信するテレビ会議システムなどに利用される。
【０００４】
このような信号処理の技術分野においては、発信信号を複数の受信装置を用いて受信して、各受信信号間の相互相関関数に基づいて信号源の方向や位置を推定することが従来から行われている。
【０００５】
例えば、特開平１１−８３９８２号公報には、複数のマイクロホンのうち選択された２個のマイクロホンの出力間における、仮定した音源方向角度に対応する遅延時間を与えて相互相関係数をそれぞれ求めた上で、これを加算した値を音源方向角度に対して表示する音源方向検知装置について開示されている。該公報に記載の音源方向検知装置によれば、ＳＮ比が悪い場合でも、音波の種類の如何を問わず、音の到来方向を的確に知ることができる。
【０００６】
また、特開平１１−３０４９０６号公報には、複数のマイクロホンで受信した信号の相互相関関数を、すべてのマイクロホンの組について計算し、これら相互相関関数の最大値を与える時間差を予備推定時間差とし、続いて、すべてのマイクロホンについての遅延和のパワーを最大にする時間差を予備推定時間差の近傍で探索してこれを推定時間差とし、この推定時間差に基づいて音源の位置を計算する音源位置推定方法について開示されている。時間差は音源の方向角度に相当する。各マイクに遅延を付けて足し合わせることにより、特定方向の受信感度を上げることができる。該公報に記載の音源位置推定方法は、耐雑音性に優れ、且つ、演算量が少なくて済む。
【０００７】
また、特許第２９８２７６６号公報には、複数のマイクロホンで獲得した複数のオーディオ信号から、信号そのもの、あるいは信号から計算した自己回帰係数を係数とする逆フィルタにより白色化した信号の極性のみを抽出した符号時間系列を生成し、これらの符号時間系列の相互相関関数を計算し、これら相互相関関数から正規化電力を計算してその時間平均を計算し、この時間平均を基に音源方向を推定する方法及びその装置について開示されている。
【０００８】
しかしながら、上記の従来技術はいずれも、雑音が多い場合や同時に信号を発生する信号源の数が多い場合には、充分に信号源の方向や位置を推定することができない。
【０００９】
特開平１１−３０４９０６号公報では、相互相関関数の最大値を与える時間差を予備推定時間差とし、この近傍で遅延和のパワーを最大にする時間差を探索して、その時間差から音源位置を計算するようにしている（上述）。この方法を複数の信号源が存在する環境下で、それぞれの信号源の方向や位置の推定に適用することを想定した場合、それぞれの信号源に対応する予備推定時間を相互相関関数から求め、各々の予備推定時間の近傍で遅延和のパワーを細大にする時間差をそれぞれ探索する必要がある。このため、探索に要する計算量が信号源の数に比例して増大してしまうことになる。
【００１０】
また、特許第２９８２７６６号公報では、装置のハードウェア規模を小さくするために、信号の相互相関関数を計算するのではなく、信号そのもの、あるいは白色化した信号の極性のみを抽出した符号時間系列の相互相関関数を基に音源方向を推定するようにしている。
【００１１】
信号そのものの極性を抽出する方法においては、比較的レベルの高い低周波の雑音が受信信号に含まれている場合、抽出された符号時間系列は−１又は＋１がその雑音の周期の２分の１程度の期間にわたって連続することになる。したがって、符号時間系列は音源の信号ではなくこの低周波雑音に対応することになるので、その相互相関関数から音源方向を求めることができない。
【００１２】
また、白色化した信号の極性を抽出する方法を利用する場合、白色化の過程において受信信号に含まれる音源からの符号に特有の性質が失われてしまう。このため、その相互相関関数は雑音の影響を強く受けるようになるので、音源方向の推定性能が劣化してしまう。なお、白色化の方法については、例えば"Improvement of the performance of cross correlation method for identifying aircraft noise with pre-whitening of signals"（The Journal of the Acoustical Society of Japan, vol.13, no.4, pp.241-252, July 1992）に記載されているが、これはそもそも飛行場近くで測定した飛行機の騒音を対象としたものである。
【００１３】
また、特許第２９８５９８２号公報には、２つのマイクロホンの出力を、それぞれ帯域分割した上で、各周波数帯域ごとに信号のパワーを求め、ピーク・ホールドし、その対数を計算して、時間差分処理した信号の対応する周波数帯域ごとの相互相関関数を求め、これらの帯域ごとの相互相関関数を重み付け平均化し、その最大値を取る時間差から音源方向を演算する音源方向推定方法について開示されている。
【００１４】
該公報に記載の方法によれば、反射音が多い環境下においても、直接音に基づいて音源方向を推定することができるとされている。しかしながら、この方法による場合、各マイクロホンからの入力信号の各周波数帯域ごとの対数計算に要する計算量が多く、これを実現するためのハードウェアが大規模化してしまう。また、信号のパワーが比較的小さい、若しくは雑音のパワーが比較的大きい場合には、対数処理によって反射音の影響を小さくすることができないこともある。例えば、暗騒音のパワーが１で、そこにパワーが２の直接音が到達し、続いてパワーが３の反射音が到達したとすると、対数化処理後の数値では、暗騒音は３．０ｄＢ、反射音は４．７ｄＢとなる。したがって、対数化処理前では反射音の大きさは直接音の１．５倍であるが、対数化処理後では１．５７倍となり、反射音の影響は数値的に軽減されない/*/。
【００１５】
雑音の多い環境や、複数の信号源が同時に信号を発生しているような環境において、雑音や、他の信号源からの信号の影響を抑制し、ある１つの信号源の信号を強調したり分離する技術は、信号が音声の場合には、音声認識装置の認識性能を高めたり、予想される信号源の種類ごとにあらかじめ測定しておいた信号と受信信号とを比較して信号源の特定を行う信号源識別装置の識別性能を高めたりする目的で利用される。
【００１６】
このような信号強調・分離技術の分野においては、複数の受信装置により信号を受信し、その受信信号間の相互相関関数などに基づいて信号源の方向又は位置と受信装置の位置によって定まる各受信装置ごとの遅延時間を推定して、この推定遅延時間を用いて各受信信号を遅延させた上で加算し、信号源の信号を強調したり、あるいは分離することが行われている。
【００１７】
例えば、特開平５−９５５９６号公報には、複数のマイクロホンにより受信された各音声信号を、それぞれ帯域通過フィルタにより周波数帯域別に分解して、周波数帯域別信号間の相互相関関数を求め、この相互相関関数から各音声信号の時間差を検出して、各音声信号をこの検出時間差に基づいて遅延した後に加算する雑音低減装置について開示されている。該雑音低減装置によれば、入力音声信号を合成することで雑音を抑えつつ音声信号を強調して取り出すことにによってＳＮ比を向上させることができる。
【００１８】
また、特開平９−２５１２９９号公報には、複数のマイクロホンからなるマイクロホン・アレイからの入力信号を、帯域通過フィルタ・バンクにより周波数帯域ごとのバンド・パス波形をマイクロホン・チャンネル別に求め、仮定する音源位置又は方向ごとのバンド・パス・パワー分布を、最小分散法などにより求め、この周波数帯域別のバンド・パス・パワー分布をすべての周波数帯域について統合して、音源位置又は方向を推定するとともに、この推定した音源位置又は方向に基づいて、周波数帯域ごとのバンド・パス・パワー分布から該当するバンド・パス・パワーを音声パラメータとして抽出して、音声認識を行うマイクロホン・アレイ入力型音声認識装置及び方法について開示されている。
【００１９】
また、特許第２９２８８７３号公報には、複数の波動収集回路により検出された波動を、周波数帯域ごとの信号又はその時間的レベル変化を信号成分として求め、これらの各周波数帯域ごとの信号成分の相互相関関数を計算して、相関値が設定された閾値を越える信号成分間の時間差を求め、この時間差が所定の設定遅れ時間に含まれる周波数帯域の信号成分を取り出し、この所定の設定遅れ時間に対応する特定の位置から到来した波動成分を各周波数帯域ごとに出力したり、あるいはこれらの信号成分を加算して、特定の位置から到来した波動成分を出力する信号処理装置について開示されている。
【００２０】
しかしながら、特開平５−９５５９６号公報に記載の雑音低減装置は、主に自動車電話に利用することが想定されており、その信号を強調すべき音源は運転者の声であって、信号源は１つと仮定されている。すなわち、複数の音源からの信号をそれぞれ強調することは考慮されておらず、また、音源のおおよその位置をあらかじめ想定しているので、任意の位置にある音源に関する処理が考慮されていない。また、強調信号は周波数帯域別に分解する前の信号をそのまま遅延して加算しているので、ある周波数帯域に比較的大きな雑音がある場合には、これを充分に抑制することができない。
【００２１】
また、特開平９−２５１２９９号公報に記載のマイクロホン・アレイ入力型音声認識装置及び方法では、仮定する音源位置又は方向ごとのバンド・パス・パワー分布を最小分散法などにより求めている。ところが、最小分散法は、音源の数がマイクロホン・アレイを構成するマイクロホンの総数より少ない場合に適用できる方法であり、またその計算量は多い。また、最小分散法に代えて遅延和法を用いれば、計算量をある程度削減することができるが、音源位置又は方向の推定精度は低くなってしまう。したがって、信号源の数、及びその位置又は方向が未知である場合、充分と考えられるだけのマイクロホンを用意するとともに仮定する信号源位置又は方向を変化させながら最小分散法などを適用するか、あるいは、仮定する信号源位置又は方向を複数設定して並列に最小分散法を適用することが必要となる。このような処理を実現するハードウェアは高価なものになるか、あるいは極めて大規模なものになる。
【００２２】
また、特許第２９２８８７３号公報に記載の信号処理装置は、所定の位置から到来する波動成分を抽出・分離することを目的として構成されたものであり、信号源の方向又は位置が未知であって信号の到来方向が未知である場合に、信号源の方向又は位置を推定して、その方向からの信号を抽出・分離したり強調する点については考慮されていない。
【００２３】
以上を総括すれば、従来の技術においては、あらかじめ信号源の方向や位置が知られておらず、且つ、雑音が多い場合や、同時に信号を発生する信号源の数が多い場合、特に受信装置の数より信号源の数が多い場合に、１つあるいは複数の信号源の方向又は位置を充分な精度で推定することができず、また、１つあるいは複数の信号源からの信号を充分に強調することができなかった。
【００２４】
《注釈》
＊：ここで言う「反射音」は、持続している直接音とそれに加えて到来する反射音を含むものとする。通常、反射音単独のパワーが直接音のパワーより大きくなることはない。指向性マイクロホンを用いている場合は、指向性の高い方向が反射音の到来方向近くを向いていて、指向性の低い方向が直接音の到来方向に向いている場合を仮定すると、反射音単独のパワーが直接音のそれよりも大きくなることもある。
【００２５】
【発明が解決しようとする課題】
本発明の目的は、音源などの信号源が存在する方向や位置を推定することができる、優れた信号源方向又は位置推定方法及び装置を提供することにある。
【００２６】
本発明の更なる目的は、存在する方向や位置があらかじめ知られていない１あるいは複数の信号源からの信号を複数の受信装置で受信した信号を用いて、信号源の方向又は位置を推定することができる、優れた信号源方向又は位置推定方法及び装置を提供することにある。
【００２７】
本発明の更なる目的は、複数の受信装置で受信される各受信信号間の相互相関関数に基づいて１つあるいは複数の信号源の信号を強調することができる、優れた信号強調方法及び装置を提供することにある。
【００２８】
本発明の更なる目的は、存在する方向や位置があらかじめ知られていない１あるいは複数の信号源からの信号を、複数の受信装置で受信した信号を用いて強調することができる、優れた信号強調方法及び装置を提供することにある。
【００２９】
【課題を解決するための手段】
本発明は、上記課題を参酌してなされたものであり、その第１の側面は、１あるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理装置又は方法であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段又はステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段又はステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３の手段又はステップと、
を具備することを特徴とする信号処理装置又は方法である。
【００３０】
前記第１の手段は、それぞれ異なる周波数帯域の信号を通過する複数の帯域通過フィルタで構成することができる。
【００３５】
また、前記信号源は音源であってもよい。
【００３６】
また、本発明の第２の側面は、１あるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理装置又は装置であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段又はステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段又はステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定する第３の手段又はステップと、
を具備することを特徴とする信号処理装置又は方法である。
【００３７】
また、本発明の第３の側面は、１あるいは複数の信号源からの信号を複数の受信装置を介して受信して得た複数の受信信号を処理する信号処理装置又は方法であって、
それぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段又はステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段又はステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３の手段又はステップと、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４の手段又はステップと、
を具備することを特徴とする信号処理装置又は方法である。
【００３８】
前記第１の手段は、それぞれ異なる周波数帯域の信号を通過する複数の帯域通過フィルタで構成することができる。
【００４４】
また、前記信号源は音源であってもよい。
【００４５】
また、本発明の第４の側面は、１あるいは複数の信号源からの信号を複数の受信装置を介して受信した信号を記録した記録媒体から再生される信号を処理する信号処理装置又は方法であって、
再生されたそれぞれの受信信号を複数の周波数帯域別信号に分解する第１の手段又はステップと、
２つの異なる受信装置に由来する周波数帯域別信号間の相互相関関数を、対応する周波数帯域ごとに、受信装置の組み合わせごとに計算する第２の手段又はステップと、
周波数帯域ごとに前記相互相関関数を規格化して、受信装置の組み合わせごとに規格化された相互相関関数をすべての周波数帯域、または事前に信号や雑音が含まれる周波数帯域が知られている場合には信号が含まれる若しくは雑音が含まれない一部の複数の周波数帯域にわたって加算し、加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求めて候補遅延時間とし、各候補遅延時間を最も整合的に与える信号源の方向又は位置を推定し、該推定した信号源の方向又は位置により各受信装置ごとの遅延時間を推定する第３の手段又はステップと、
該推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、該遅延した周波数帯域別信号を加算する第４の手段又はステップと、
を具備することを特徴とする信号処理装置又は方法である。
【００４６】
【作用】
１あるいは複数の信号源からの信号を、複数の受信装置で受信して、各受信信号を複数の帯域通過フィルタにより周波数帯域別の信号に分解する。そして、２つの異なる受信装置に由来するこれらの周波数帯域別信号間の相互相関関数を、相当する周波数帯域ごとで且つ受信装置の組合せごとに計算する。方向性のない雑音のパワーがある一部の周波数帯域において大きい場合には、その周波数帯域における相互相関関数は極大値を示さないので、信号源の方向又は位置、あるいは信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定する際に、この雑音の影響を効果的に抑制することができる。
【００４７】
また、方向性のない雑音のパワーが、すべての周波数帯域において大きい場合であっても、信号源からの信号が存在する周波数帯域においては、その周波数帯域内で信号が雑音により打ち消されてしまわない限り、その周波数帯域の相互相関関数は、信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間に極大値を持つので、信号源の方向又は位置、あるいは信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定することができる。
【００４８】
また、方向性のある雑音の場合は、これも信号源の１つとして、その方向又は位置、あるいはその方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定することができる。
【００４９】
また、相互相関関数から信号源の方向又は位置を推定する方法、あるいは信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定する方法において、周波数帯域ごとに相互相関関数を規格化し、受信装置の組合せごとにこれらの規格化した相互相関関数を、すべてのあるいは一部の複数の周波数帯域にわたって加算する方法を用いることによって、周波数帯域によってパワー・レベルが異なる信号の場合に、パワー・レベルが低い周波数帯域の情報も、パワー・レベルが高い他の周波数帯域の信号に妨害されることなく、有効に利用することができるので、信号源の方向又は位置、あるいは信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定する際の精度を高めることができる。
【００５０】
また、受信装置の組合せごとにこれらの規格化した相互相関関数を複数の周波数帯域にわたって加算する方法を用いることによって、これら加算された相互相関関数が極大となる遅延時間を、受信装置の組合せごとに求め、これらを候補遅延時間とするので、加算を用いない方法と比較して候補遅延時間の数を削減することができる。すなわち、これらの候補遅延時間を最も整合的に与える信号源の数及び方向又は位置を推定する際の計算量を減らすこと、又はその計算に必要なハードウェアの規模を小さくすることができる。
【００５１】
また、上述した相互相関関数の加算においては、事前に信号や雑音に関する性質が知られている場合などには、信号による成分を強調して、雑音による成分を抑制するように重み付き加算することによって、信号源の方向又は位置、あるいは信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定する際の精度をいっそう高めることができる。
【００５２】
また、本発明に係る信号強調方法及び装置において、遅延した周波数帯域別信号を加算する際に、推定した遅延時間における周波数帯域別信号間の相互相関関数の値に応じて、周波数帯域ごとに重み係数を定め、それらの重み係数を用いて重み付き加算する場合は、歪みが生じる欠点はあるが、雑音や他の信号源からの影響を小さくすることができ、強調をより明確に行うことができる。
【００５３】
本発明のさらに他の目的、特徴や利点は、後述する本発明の実施例や添付する図面に基づくより詳細な説明によって明らかになるであろう。
【００５４】
【発明の実施の形態】
以下、図面を参照しながら本発明の実施例を詳解する。
【００５５】
［第１の実施例］
図１には、本発明の第１の実施例に係る信号処理装置の機能構成を模式的に示している。
【００５６】
同図に示すように、この信号処理装置は、１以上の信号源１−１，１−２…から発生する信号を受信する複数（Ｍ個）のマイクロホン２−１，２−２，…，２−Ｍと、各マイクロホン２−１…の受信信号を増幅する増幅器３−１，３−２，…，３−Ｍと、各増幅信号をデジタル変換するＡ／Ｄ変換器４−１，４−２，…，４−Ｍと、帯域通過フィルタ・バンク５−１，５−２，…，５−Ｍとからなる複数（図示の例ではＭ個）の受信系統と、対応する周波数帯域ごとに受信系統の組み合わせごとに周波数帯域別信号間の相互相関関数を計算する相互相関関数計算部６と、これら算出された相互相関関数を基に信号源の方向又は位置を推定する信号源方向又は位置推定部７とで構成される。
【００５７】
マイクロホン２−１，２−２，…，２−Ｍでそれぞれ受信され、増幅器３−１，３−２，…，３−Ｍでそれぞれ増幅された受信信号は、Ａ／Ｄ変換器４−１，４−２，…，４−Ｍによりそれぞれデジタル変換されて、デジタル受信信号ｘ（ｍ；ｔ）を得る。ここで、ｍ（＝１，２，…，Ｍ）は、マイクロホンを識別するための参照番号であり、また、ｔ（＝１，２，…）は信号が劣化しないように選ばれたサンプリング周期を単位とする時刻である。
【００５８】
次いで、デジタル変換された各受信信号ｘ（ｍ，ｔ）は、それぞれ対応する帯域通過フィルタ・バンク部５−ｍ（ｍ＝１，２，…，Ｍ）に入力される。帯域通過フィルタ・バンク部５−ｍは、それぞれ複数の帯域通過フィルタで構成されて、増幅及びデジタル変換された受信信号を複数の周波数帯域別信号に分解することができる。ここでは、Ｆ（＞１）個の帯域別信号ｘ（ｍ，１；ｔ），ｘ（ｍ，２；ｔ），…，ｘ（ｍ，Ｆ；ｔ）に分解されて出力されるものとする。Ｆは、各フィルタ・バンク部のバンク数である。
【００５９】
次いで、Ｍ×Ｆ個の帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）は相互相関関数計算部６に入力される。相互相関関数計算部６では、異なる２つのマイクロホンの組み合わせ（ｍ，ｎ）ごとに、対応する周波数帯域別の相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ；ｄ＝−Ｔ，−Ｔ＋１，…，Ｔ）を計算して出力する。ここで、ｄ（＝−Ｔ，−Ｔ＋１，…，Ｔ）は、ｘ（ｎ，ｆ；ｔ）を基準としてｘ（ｍ，ｆ，ｔ）を遅延させたと考える遅延時間であり、Ｔは以下で説明するように、マイクロホンの間隔と信号の速さとサンプリング周期から定める最大遅延時間である。
【００６０】
次いで、算出されたＭ×（Ｍ−１）×Ｆ／２個の帯域別相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｔ＝１，２，…；ｄ＝−Ｔ，−Ｔ＋１，…，Ｔ）は信号源方向又は位置推定部７に入力されて、信号源の方向又は位置が推定されて出力される。
【００６１】
図２には、帯域通過フィルタ・バンク部５−ｍの内部構成を図解している。同図に示すように、帯域通過フィルタ・バンク部５−ｍは、それぞれ異なる周波数帯域を受け持つＦ個の帯域通過フィルタ８−１，８−２，…，８−Ｆで構成される。各帯域通過フィルタ８−１，８−２，…，８−Ｆは、例えばデジタル・フィルタで構成される。
【００６２】
ｍ番目のマイクロホン２−ｍからの受信信号に由来する信号ｘ（ｍ；ｔ）は、それぞれ異なる周波数帯域を受け持つＦ個の帯域通過フィルタ８−１，８−２，…，８−Ｆによりフィルタリングされることにより、帯域別のＦ個の信号ｘ（ｍ，ｆ；ｔ）（ｆ＝１，２，…，Ｆ）に分解されて出力される。
【００６３】
図３には、相互相関関数計算部６の内部構成を図解している。同図では、相互相関関数係数部６のうち、ｍ番目のマイクロホン５−ｍとｎ番目のマイクロホン５−ｎの各々からの受信信号間の相互相関関数を計算する部分を抽出して図解している。
【００６４】
ｍ番目のマイクロホン２−ｍとｎ番目のマイクロホン２−ｎの各々からの受信信号におけるｆ番目の周波数帯域の各帯域別信号ｘ（ｍ，ｆ；ｔ）及びｘ（ｎ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）は、それぞれＴ個の遅延器９１，９−２，…，９−Ｔ、及び１０−１，１０−２，…，１０−Ｔに順次入力される。
【００６５】
特定量の遅延を受けたこれら２つの信号は、積算器１１−１，１１−２，…，１１−（２Ｔ＋１）により掛け合わされ、さらに低域通過フィルタ１２−１，１２−２，…，１２−（２Ｔ＋１）に入力され、重み付き時間平均され、その結果がｆ番目の周波数帯域のｍ番目とｎ番目のマイクロホン２−ｍ及び２−ｎからの帯域別信号の各時間後とに２Ｔ＋１個の値からなる相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｔ＝１，２，…，；ｄ＝−Ｔ−ｔ＋１，…，Ｔ）として出力される。
【００６６】
遅延器の個数Ｔは、ｍ番目とｎ番目の各マイクロホン２−ｍ，２−ｎ間の距離をＤ（ｍ，ｎ）、信号の速さをｃとおいたとき、サンプリング周期を時間の単位とすれば、Ｔ〜Ｄ（ｍ，ｎ）／ｃとなるように選ばれる。後の過程における処理を単純化するなどの目的のために、Ｔより多くの遅延器を設けてもよい。但し、その場合はマイクロホン間の距離Ｄ（ｍ，ｎ）を信号が進む時間Ｄ（ｍ，ｎ）／ｃより長い遅延を与えることになり、物理的には冗長である。
【００６７】
次いで、信号源方向又は位置推定部７の内部構成及び動作特性について説明する。但し、後述するようにマイクロホンの個数Ｍが２の場合と３以上の場合とでは内部構成と動作は相違する。
【００６８】
図４には、マイクロホンの個数Ｍが２の場合における信号源方向又は位置推定部７の内部構成を図解している。
【００６９】
同図に示す場合、２つの異なるマイクロホンの組は参照番号１−１と１−２のマイクロホンによる一組だけであるので、信号源方向又は位置推定部７への入力はＦ個の相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）である。これらの相互相関関数は、規格化計算部１３により周波数帯域ごとに規格化され、規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）として出力される。
【００７０】
規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）は、相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）の遅延時間ｄに関して最大値又は分散などの値を求め、これらの値によりその相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）を割り算して計算することができる。
【００７１】
次いで、規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）は、加算部１４に入力され、それらのフィルタの番号ｆについて加算した値が加算相互相関関数Ｒ（１，２；ｔ；ｄ）として出力される。
【００７２】
この加算においては、通常はすべてのフィルタ番号について加算することが望ましい。あるいは、事前に信号や雑音が含まれる周波数帯域が知られている場合には、信号が含まれる若しくは雑音が含まれない、一部の複数の周波数帯域に対応する番号についてのみ加算してもよいし、あるいは信号による成分を強調して雑音による成分を抑制するように重み付き加算してもよい。
【００７３】
次いで、加算相互相関関数Ｒ（１，２；ｔ；ｄ）は、極大検出部１５に入力され、遅延時間ｄに関して極大値を与える遅延時間ｄ（１，２；１），ｄ（１，２；２），…が検出され、出力される。
【００７４】
この極大検出においては、信号源に対応しない意味のない極大を検出することを抑制するため、あらかじめ設定されるか、あるいは加算相互相関関数Ｒ（１，２；ｔ；ｄ）の、平均や分散などの統計的数値によって設定される閾値に、その極大値が達している、あるいはその閾値を越えていることを条件として、極大の検出を行うようにしてもよい。
【００７５】
次いで、上記の遅延時間ｄ（１，２；１），ｄ（１，２；２），…は信号源方向又は位置推定部１６に入力され、それぞれの遅延時間に対応する信号源の方向θ1，θ2，…が以下の式により計算され、出力される。
【００７６】
【数１】
0004815661.tif 000002
【００７７】
図５には、信号源方向推定に用いられる上式［数１］を図解している。但し、２つのマイクロホン２−１と２−２と信号源１−ｉを通る平面内で説明する。
【００７８】
マイクロホンが２つの場合、信号源の位置を遅延時間ｄ（１，２；ｉ）から推定することは数学的に不可能であり、その方向を推定することだけができる。この場合、信号源は２つのマイクロホン間の距離Ｄ（１，２）に比べて充分に離れた位置にあり、信号は平面波で記述されるものと近似される。
【００７９】
図５に示すようなマイクロホン２−１と２−２の作る直線に垂直な方向を基準として信号源１−ｉがθiの方向にあるものとする。このとき、θiの符号は信号源がマイクロホン２−１側にある場合を正方向、マイクロホン２−２側にある場合を負方向と定義する。したがって、図５に示す例ではθi＞０である。
【００８０】
図５に示す例では、信号の速さをｃと置くと、信号源１−ｉからの信号は、マイクロホン２−１に先に到達し、マイクロホン２−２にはＤ（１，２）ｓｉｎθi／ｃだけ遅れて到達する。
【００８１】
したがって、マイクロホン２−１及び２−２によって受信した信号ｘ（１；ｔ）及びｘ（２；ｔ）のうち信号源１−ｉに由来する成分をそれぞれｘｓ（１，ｉ；ｔ）及びｘｓ（２，ｉ；ｔ）と定義すると、時間遅れに着目し振幅の違いを無視すると、以下の式が成立する。
【００８２】
【数２】
0004815661.tif 000003
【００８３】
加算相互相関関数Ｒ（１，２；ｔ；ｄ）の極大として検出された遅延時間ｄ（１，２；ｉ）が信号源１−ｉに対応するものとすれば、
【００８４】
【数３】
0004815661.tif 000004
【００８５】
となるから、これを変形して上式［数１］を得ることができる。
【００８６】
なお、雑音による影響やサンプリングによる影響のためｃ・ｄ（１，２；ｉ）／Ｄ（１，２）の絶対値|ｃ・ｄ（１，２，ｉ）／Ｄ（１，２）|が１を越え、そのａｒｃｓｉｎが実数の範囲内では求められない事態が起こり得る。この場合は、ｄ（１，２；ｉ）が正のとき、θi＝π／２＝９０゜、ｄ（１，２；ｉ）が負のときθi＝−π／２＝−９０゜と推定して差し支えない。
【００８７】
次に、マイクロホンの個数Ｍが３以上の場合について、信号源方向又は位置推定部７の内部構成及び動作について、図６を参照しながら説明する。
【００８８】
マイクロホンの個数が２の場合（図４を参照のこと）と同様に、相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）は規格化計算部１３に入力され、周波数帯域ごとに規格化され、規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）として出力される。
【００８９】
規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）は、相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）の遅延時間ｄに関して最大値又は分散等の値を求め、これらの値によりその相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）を割り算して計算する。
【００９０】
次いで、規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）は、加算部１４に入力され、それらをフィルタの参照番号ｆについて加算した値が加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）として出力される。
【００９１】
加算部１４での加算においては、通常は、すべてのフィルタ番号について加算することが望ましい。あるいは、事前に信号や雑音が含まれる周波数帯域が知られている場合には、信号が含まれる若しくは雑音が含まれない、一部の複数の周波数帯域に対応する番号についてのみ加算してもよいし、あるいは信号による成分を強調し雑音による成分を抑制するように重み付き加算してもよい。
【００９２】
次いで、加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）は、極大検出部１５に入力され、遅延時間ｄに関して極大値を与える遅延時間ｄ（ｍ，ｎ；１），ｄ（ｍ，ｎ；２），…（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）が検出され、出力される。
【００９３】
この極大検出においては、信号源に対応しない意味のない極大を検出することを抑制するため、あらかじめ設定される、あるいは加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）の平均や分散などの統計的数値によって設定される閾値を、その極大値が達しているあるいは越えていることを条件として、極大値の検出を行うようにしてもよい。
【００９４】
次いで、上記の遅延時間ｄ（ｍ，ｎ；ｉ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｉ＝１，２，…）は、信号源位置推定部１７に入力され、それぞれの遅延時間に対応する信号源のＳ（１）＝（Ｓx（１），Ｓy（１），Ｓz（１）），Ｓ（２）＝（Ｓx（２），Ｓy（２），Ｓz（２）），…が計算され、出力される。
【００９５】
図７には、信号源位置推定部１７の内部構成を示している。
【００９６】
上記の遅延時間ｄ（ｍ，ｎ；ｉ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｉ＝１，２，…）は、遅延時間分析部１８に入力され、推定される信号源ごとに集合Ｅ（ｋ）＝｛ｄ（ｍ1，ｍ2；ｋ），ｄ（ｍ1，ｍ3；ｋ），…，ｄ（ｍ2，ｍ3；ｋ），…，ｄ（ｍi，ｍj；ｋ）｝（ｋ＝１，２，…）に分類される。
【００９７】
このとき、遅延時間ｄ（ｍ，ｎ；ｉ）の異なる遅延時間を区別する番号ｉを適宜入れ替え、異なる上記推定される信号源を区別する番号ｋに一致するようにした。この分類は、ある１つの信号源に対応する集合Ｅ（ｋ）から、相異なる３つのマイクロホンの組（ｎ1，ｎ2，ｎ3）（但し、ｎ1＜ｎ2＜ｎ3とする）の要素から作られる３つの相異なる２つのマイクロホンの組（ｎ1，ｎ2），（ｎ1，ｎ3），（ｎ2，ｎ3）の受信信号間の遅延時間ｄ（ｎ1，ｎ2；ｋ），ｄ（ｎ1，ｎ3；ｋ），ｄ（ｎ2，ｎ3；ｋ）を任意に選んだときに、次の条件式が満たされるように行われる。
【００９８】
【数４】
0004815661.tif 000005
【００９９】
ここで、|ａ|はａの絶対値を表し、εは１，２，３程度の小さな数に選ばれる。
【０１００】
上記の分類を行うには、例えば、すべての相異なるマイクロホンの番号の組（ｎ1，ｎ2，ｎ3）（１≦ｎ1＜ｎ2＜ｎ3≦Ｍ）について、上式［数４］を満たす遅延時間の組を順次探索し、さらに、共通する遅延時間を含む組を併合して、その共通する遅延時間を取り除く手続を繰り返せばよい。
【０１０１】
図８には、上式［数４］を説明するための図を示している。
【０１０２】
信号源１−ｋがＳ（ｋ）に位置するものとし、そこから発生する信号が、位置Ｍｎ1，Ｍｎ2，Ｍｎ3にそれぞれ位置するマイクロホン２−ｎ1，２−ｎ2，２−ｎ3の各々に到達すると、これらに要する時間はそれぞれ、ｔｎ1（ｋ）＝|Ｓ（ｋ）−Ｍｎ1|／ｃ，ｔｎ2（ｋ）＝|Ｓ（ｋ）−Ｍｎ2|／ｃ，ｔｎ3（ｋ）＝|Ｓ（ｋ）−Ｍｎ3|／ｃ，となる。ここで、|Ａ|は、ベクトルＡの大きさである。
【０１０３】
遅延時間ｄ（ｎi，ｎj：ｋ）（ｎi，ｎj＝ｎ1，ｎ2，ｎ3）は、ｄ（ｎi，ｎj；ｋ）＝ｔｎj（ｋ）−ｔｎi（ｋ）となる。これを上式［数４］の左辺に代入すると０になることが分かる。
【０１０４】
遅延時間ｄ（ｎi，ｎj：ｋ）（ｎi，ｎj＝ｎ1，ｎ2，ｎ3）の推定において、雑音による影響やサンプリングによる影響のために生じる誤差により、上式［数４］の左辺の値が厳密には０にならない場合があることを考慮して、ある小さな値εよりも小さければ、同一の信号源ｋによるものであると判定することにしたものが［数４］である。
【０１０５】
次に、上記の信号源ごとに分類された遅延時間の集合Ｅ（ｋ）（ｋ＝１，２，…）は、信号源位置計算部１９に入力され、信号源の位置Ｓ（ｋ）＝（Ｓｘ（ｋ），Ｓｙ（ｋ），Ｓｚ（ｋ））（ｋ＝１，２，…）が計算され、出力される。
【０１０６】
信号源源位置の計算方法の１つは、ｉ番目のマイクロホン１−ｉの位置をＭi＝（Ｍｘi，Ｍｙi，Ｍｚi）とおいて、Ｅ（ｋ）の各要素ｄ（ｍi，ｍj；ｋ）について、それぞれ［|Ｍｍj−Ｓ（ｋ）|／ｃ−|Ｍｍi−Ｓ（ｋ）|／ｃ−ｄ（ｍi，ｍj；ｋ）］の２乗を計算して、これらをすべて足し合わせた量を最小化するＳ（ｋ）を求めることである。
【０１０７】
次いで、このような信号源位置の計算方法によりＳ（ｋ）が正しく計算される理由について説明する。
【０１０８】
信号源１−ｋからの信号が、マイクロホン２−ｍi,２−ｍjの各々に到達するのに、それぞれｔｍi（ｋ），ｔｍj（ｋ）だけの時間を要するとすると、ｔｍi（ｋ）＝|Ｍｍi−Ｓ（ｋ）|／ｃ，ｔｍj（ｋ）＝|Ｍｍj−Ｓ（ｋ）|／ｃである。
【０１０９】
マイクロホン２−ｍi,２−ｍj間の遅延時間ｄ（ｍi，ｍj；ｋ）は、雑音による影響やサンプリングによる影響のために生じる誤差がなければ、ｄ（ｍi，ｍj；ｋ）＝ｔｍj（ｋ）−ｔｍi（ｋ）＝|Ｍｍj−Ｓ（ｋ）|／ｃ−|Ｍｍi−Ｓ（ｋ）|／ｃとなる。
【０１１０】
雑音による影響やもサンプリングによる影響のために生じる誤差を考慮して、このような誤差の２乗をすべてのＥ（ｋ）の要素である遅延時間に関して足し合わせた量をＳ（ｋ）に関して最小化すれば、最小自乗誤差基準に基づいて最適なＳ（ｋ）の推定値が得られる。
【０１１１】
なお、上述した信号源位置推定において、マイクロホンの個数Ｍが３の場合、空間的に任意の位置に存在する信号源の位置を推定することは数学的に不可能である。このため、信号源が存在する平面を仮定して、これを信号源の位置に関する拘束条件として位置を推定すればよい。この場合、計算が正しく行われるためには、この信号源が存在すると仮定する平面は、３つのマイクロホンがその上に存在する平面に対して垂直に交じらないものを選ぶ必要があり、例えば、３つのマイクロホンがその上に存在する平面を選べばよい。
【０１１２】
［第２の実施例］
図９には、本発明の第２の実施例に係る信号処理装置の構成を模式的に示している。
【０１１３】
同図に示すように、この信号処理装置は、１以上の信号源１−１，１−２…から発生する信号を受信する複数（Ｍ個）のマイクロホン２−１，２−２，…，２−Ｍと、各マイクロホン２−１…の受信信号を増幅する増幅器３−１，３−２，…，３−Ｍと、各増幅信号をデジタル変換するＡ／Ｄ変換器４−１，４−２，…，４−Ｍと、帯域通過フィルタ・バンク５−１，５−２，…，５−Ｍとからなる複数（Ｍ系統）の受信系統と、対応する周波数帯域ごとに受信系統の組み合わせごとに周波数帯域別信号間の相互相関関数を計算する相互相関関数計算部６と、信号源からの信号を強調する信号強調部２０とで構成される。
【０１１４】
マイクロホン２−１，２−２，…，２−Ｍでそれぞれ受信され、増幅器３−１，３−２，…，３−Ｍでそれぞれ増幅された受信信号は、Ａ／Ｄ変換器４−１，４−２，…，４−Ｍによりデジタル変換されて、デジタル受信信号ｘ（ｍ；ｔ）を得る。ここで、ｍ（＝１，２，…，Ｍ）は、マイクロホンを識別するための参照番号であり、また、ｔ（＝１，２，…）は信号が劣化しないように選ばれたサンプリング周期を単位とする時刻である。
【０１１５】
次いで、デジタル変換された各受信信号ｘ（ｍ，ｔ）は、それぞれ対応する帯域通過フィルタ・バンク部５−ｍ（ｍ＝１，２，…，Ｍ）に入力される。帯域通過フィルタ・バンク部５−ｍは、それぞれ複数の帯域通過フィルタで構成されて、増幅及びデジタル変換された受信信号を複数の周波数帯域別信号に分解することができる。ここでは、Ｆ（＞１）個の帯域別信号ｘ（ｍ，１；ｔ），ｘ（ｍ，２；ｔ），…，ｘ（ｍ，Ｆ；ｔ）に分解されて出力されるものとする。Ｆは、各フィルタ・バンク部のバンク数である。
【０１１６】
次いで、Ｍ×Ｆ個の帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）は相互相関関数計算部６に入力される。相互相関関数計算部６では、異なる２つのマイクロホンの組み合わせ（ｍ，ｎ）ごとに、対応する周波数帯域別の相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ；ｄ＝−Ｔ，−Ｔ＋１，…，Ｔ）を計算して出力する。ここで、ｄ（＝−Ｔ，−Ｔ＋１，…，Ｔ）は、ｘ（ｎ，ｆ；ｔ）を基準としてｘ（ｍ，ｆ，ｔ）を遅延させたと考える遅延時間であり、Ｔは［００６］で既に説明したように、マイクロホンの間隔と信号の速さとサンプリング周期から定める最大遅延時間である。
【０１１７】
次いで、算出されたＭ×（Ｍ−１）×Ｆ／２個の帯域別相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｔ＝１，２，…；ｄ＝−Ｔ，−Ｔ＋１，…，Ｔ）は信号強調部２０に入力されて、１つあるいは複数の信号源ごとの信号が強調されて出力される。
【０１１８】
図９において、図１と同じ参照番号が付された各機能モジュールは略同一の内部構成及び動作特性を備えているので、ここでは説明を省略する。以下では、信号強調部２０について、その内部構成と動作特性について説明する。但し、マイクロホンの個数Ｍが２の場合と３以上の場合とでは内部構成と動作は相違する。
【０１１９】
図１０には、マイクロホンの個数Ｍが２の場合における信号強調部２０の内部構成を図解している。
【０１２０】
同図に示す場合、２つの異なるマイクロホンの組は参照番号１−１と１−２のマイクロホンによる一組だけであるので、信号源方向又は位置推定部７への入力はＦ個の相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）である。これらの相互相関関数は、規格化計算部１３により周波数帯域ごとに規格化され、規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）として出力される。
【０１２１】
規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）は、相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）の遅延時間ｄに関して最大値又は分散などの値を求め、これらの値によりその相互相関関数ｒ（１，２，ｆ；ｔ；ｄ）を割り算して計算することができる。
【０１２２】
次いで、規格化相互相関関数ｒｎ（１，２，ｆ；ｔ；ｄ）（ｆ＝１，２，…，Ｆ）は、加算部１４に入力され、それらのフィルタの番号ｆについて加算した値が加算相互相関関数Ｒ（１，２；ｔ；ｄ）として出力される。
【０１２３】
この加算においては、通常はすべてのフィルタ番号について加算することが望ましい。あるいは、事前に信号や雑音が含まれる周波数帯域が知られている場合には、信号が含まれる若しくは雑音が含まれない、一部の複数の周波数帯域に対応する番号についてのみ加算してもよいし、あるいは信号による成分を強調して雑音による成分を抑制するように重み付き加算してもよい。
【０１２４】
次いで、加算相互相関関数Ｒ（１，２；ｔ；ｄ）は、極大検出部１５に入力され、遅延時間ｄに関して極大値を与える遅延時間ｄ（１，２；１），ｄ（１，２；２），…が検出され、出力される。
【０１２５】
この極大検出においては、信号源に対応しない意味のない極大を検出することを抑制するため、あらかじめ設定されるか、あるいは加算相互相関関数Ｒ（１，２；ｔ；ｄ）の、平均や分散などの統計的数値によって設定される閾値に、その極大値が達している、あるいはその閾値を越えていることを条件として、極大の検出を行うようにしてもよい。
【０１２６】
極大検出部１５で検出された、遅延時間ｄ（１，２；１），ｄ（１，２；２），…と、帯域通過フィルタ・バンク部５−１，５−２で分解された帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｆ）とは、強調信号生成部２１に入力され、強調信号ｙ（１；ｔ），ｙ（２；ｔ），…が、遅延器と加算器を用いた遅延加算により計算されて出力される。
【０１２７】
上記の強調信号ｙ（ｉ；ｔ）（ｉ＝１，２，…）は、例えば、図示しない遅延器により遅延時間ｄ（１，２；ｋ）に応じて、帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２；ｆ＝１，２，…，Ｆ）をそれぞれｘ（１，ｆ；ｔ−Ｔ／２−ｄ（１，２；ｉ）／２），ｘ（２，ｆ；ｔ−Ｔ／２＋ｄ（１，２；ｉ）／２）のように遅延し、さらに、図示しない加算器により、以下の加算を計算すればよい。すなわち、
【０１２８】
【数５】
0004815661.tif 000006
【０１２９】
上記の加算においては、遅延させた帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２；ｆ＝１，２，…，Ｆ）を、同じ重みで加算しているが、例えばｄ（１，２；ｉ）やｒ（１，２，ｆ；ｔ；ｄ（１，２；ｉ））の値に応じて適当な重みを乗じて加算したりしてもよい。
【０１３０】
また、図１１には、マイクロホンの個数Ｍが３以上の場合における信号強調部２０の内部構成を図解している。
【０１３１】
マイクロホンの個数が２の場合（図１０を参照のこと）と同様に、相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）は規格化計算部１３に入力され、周波数帯域ごとに規格化され、規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）として出力される。
【０１３２】
規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）は、相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）の遅延時間ｄに関して最大値又は分散等の値を求め、これらの値によりその相互相関関数ｒ（ｍ，ｎ，ｆ；ｔ；ｄ）を割り算して計算する。
【０１３３】
次いで、規格化相互相関関数ｒｎ（ｍ，ｎ，ｆ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｆ＝１，２，…，Ｆ）は、加算部１４に入力され、それらをフィルタの参照番号ｆについて加算した値が加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）として出力される。
【０１３４】
加算部１４での加算においては、通常は、すべてのフィルタ番号について加算することが望ましい。例えば、事前に信号や雑音が含まれる周波数帯域が知られている場合には、信号が含まれる若しくは雑音が含まれない、一部の複数の周波数帯域に対応する番号についてのみ加算してもよいし、あるいは信号による成分を強調し雑音による成分を抑制するように重み付き加算してもよい。
【０１３５】
次いで、加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）は、極大検出部１５に入力され、遅延時間ｄに関して極大値を与える遅延時間ｄ（ｍ，ｎ；１），ｄ（ｍ，ｎ；２），…（ｍ＝１，２，…，Ｍ；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ）が検出され、出力される。
【０１３６】
この極大検出においては、信号源に対応しない意味のない極大を検出することを抑制するため、あらかじめ設定される、あるいは加算相互相関関数Ｒ（ｍ，ｎ；ｔ；ｄ）の平均や分散などの統計的数値によって設定される閾値を、その極大値が達しているあるいは越えていることを条件として、極大値の検出を行うようにしてもよい。
【０１３７】
極大検出部１５で検出された遅延時間ｄ（１，２；１），ｄ（１，２；２），…と、帯域通過フィルタ・バンク部５−１，５−２，…，５−Ｍで分解された帯域別信号ｘ（ｍ，ｔ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）とは、強調信号生成部２２に入力され、強調信号ｙ（１；ｔ），ｙ（２；ｔ），…が、置換器と加算器を用いた遅延加算により計算されて出力される。
【０１３８】
図１２には、この強調信号生成部２２の内部構成を図解している。
【０１３９】
上記の遅延時間ｄ（ｍ，ｎ；ｉ）（ｍ＝１，２，…，Ｍ−１；ｎ＝ｍ＋１，ｍ＋２，…，Ｍ；ｉ＝１，２，…）は、遅延時間分析部１８に入力され、推定される信号源ごとに集合Ｅ（ｋ）＝ｄ（ｍ1，ｍ2；ｋ），ｄ（ｍ1，ｍ3；ｋ），…，ｄ（ｍ2，ｍ3；ｋ），…，ｄ（ｍi，ｍj；ｋ）（ｋ＝１，２，…）に分類される。
【０１４０】
このとき、遅延時間ｄ（ｍ，ｎ；ｉ）の異なる遅延時間を区別する番号ｉを適宜入れ替え、異なる上記推定される信号源を区別する番号ｋに一致するようにした。この分類は、ある１つの信号源に対応する集合Ｅ（ｋ）から、相異なる３つのマイクロホンの組（ｎ1，ｎ2，ｎ3）（但し、ｎ1＜ｎ2＜ｎ3とする）の要素から作られる３つの相異なる２つのマイクロホンの組（ｎ1，ｎ2），（ｎ1，ｎ3），（ｎ2，ｎ3）の受信信号間の遅延時間ｄ（ｎ1，ｎ2；ｋ），ｄ（ｎ1，ｎ3；ｋ），ｄ（ｎ2，ｎ3；ｋ）を任意に選んだときに、上述した条件式［数４］が満たされるように行われる。
【０１４１】
次いで、信号源ごとに分類された遅延時間の集合Ｅ（ｋ）（ｋ＝１，２，…）と、各帯域通過フィルタ・バンク部５−１，５−２，…，５−Ｍで分解された帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）とは、強調信号計算部２２に入力され、強調信号ｙ（１；ｔ），ｙ（２；ｔ），…が計算され、出力される。
【０１４２】
図１３には、この強調信号計算の方法の一例を図解している。同図に示す例では、Ｅ（ｋ）＝｛ｄ（ｍ1，ｍ2；ｋ），ｄ（ｍ2，ｍ3；ｋ），ｄ（ｍ2，ｍ4；ｋ），ｄ（ｍ3，ｍ4；ｋ）｝（ｍ1＜ｍ2＜ｍ3＜ｍ4）である。このとき、マイクロホン２−ｍiの番号ｍiを頂点ｍiに割り当て、異なる２つのマイクロホン番号ｉ＜ｊに対して、ｄ（ｉ，ｊ；ｋ）がＥ（ｋ）の要素であれば、ｉを始点とし、ｊを終点とする有向線分を表す矢印を描いたものが、図１３に示すグラフである。
【０１４３】
第１の実施例において説明した信号源ごとに分類された遅延時間の集合Ｅ（ｋ）（ｋ＝１，２，…）を求める手続に従えば、グラフは必ず連結したもの、すなわちグラフの任意の２つの頂点はいくつかの有向成分により繋がったものになる。
【０１４４】
頂点を１つ選び（図示の例ではｍ1を選択したものとする）、各頂点ｍi（ｉ＝１，２，３，４）に対して、遅延量ｚｍi（ｋ）を次のように計算する。
【０１４５】
すなわち、ｍ1から有向線分ｄ（ｉ，ｊ；ｋ）をたどって、ｍiに至る経路を1つ選び、その経路に従ってｄ（ｉ，ｊ；ｋ）の値を加算する。この加算の際、経路を進む向きと有向線分の向きが同じであれば−１を、反対であれば＋１を重みとして乗じて加算を行う。図示の例では、以下に示す式の通りとなる。すなわち、
【０１４６】
【数６】
0004815661.tif 000007
【０１４７】
上述の遅延量ｚｍi（ｋ）の計算においては、１つの経路を選んで計算を行っているが、複数の経路について計算を行い、それらの平均を適当に計算して、これを遅延量ｚｍi（ｋ）として用いてもよい。
【０１４８】
次いで、ｚｍi（ｋ）が正又は０になるように、すべてのｚｍi（ｋ）に等しくある数Ｚ（ｋ）を加えたものを改めてｚｍi（ｋ）とする。図示の場合、
【０１４９】
【数７】
0004815661.tif 000008
【０１５０】
最後に、これらの遅延量を用いて、図示しない遅延器によりｘ（ｍi，ｆ；ｔ）をそれぞれ遅延し、さらに図示しない加算器により、
【０１５１】
【数８】
0004815661.tif 000009
【０１５２】
と加算され、遅延加算された強調信号ｙ（ｋ；ｔ）が計算される。
【０１５３】
上記の遅延加算により、信号源１−ｋからの信号がマイクロホン２−ｍ1，２−ｍ2，２−ｍ3，２−ｍ4に到達するのに要する時間の差が、遅延量ｚｍ1（ｋ），ｚｍ2（ｋ），ｚｍ3（ｋ），ｚｍ4（ｋ）による遅延により補償され、同時刻に発せられた信号として加算されることは、その計算方法から明らかである。
【０１５４】
また、上記の強調信号計算方法の他の例では、第１の実施例において説明したものと同様の方法により、信号源位置推定法により各信号源の位置Ｓ（ｋ）（ｋ＝１，２，…）を推定し、この推定信号源位置Ｓ（ｋ）と各マイクロホンの位置Ｍi（ｉ＝１，２，…，Ｍ）とから計算される遅延時間ｚi（ｋ）＝|Ｓ（ｋ）−Ｍi|／ｃを用いて、帯域通過フィルタ・バンク部５−１，５−２，…，５−Ｍで分解された帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）を、以下の計算式に従い、遅延加算により強調信号ｙ（ｋ；ｔ）を計算することができる。すなわち、
【０１５５】
【数９】
0004815661.tif 000010
【０１５６】
上記の２通りの強調信号計算の方法における遅延加算においては、加算される遅延させた帯域別信号ｘ（ｍ，ｆ；ｔ）（ｍ＝１，２，…，Ｍ；ｆ＝１，２，…，Ｆ）を、すべて同じ重みで加算しているが、例えばｄ（ｍi，ｍj；ｋ）の値に応じて適当な重みを乗じて計算したり、あるいはｒ（ｍi，ｍj，ｆ；；ｔ；ｄ（ｍi，ｍj；ｋ））の値に応じて適当な重みを乗じて加算したりしてもよい。
【０１５７】
以下では、本実施例に係る信号源方向又は位置推定装置、並びに信号強調装置の適用例を挙げておく。
【０１５８】
（適用例１）
信号源方向又は位置推定装置のＴＶ会議システム、監視システムへの適用：
複数のマイクを含んだ信号源方向又は位置推定装置の出力をカメラ制御装置に接続する（図１４を参照のこと）ことで、カメラを常時信号源に向けることができる。あるいは、カメラ制御装置に複数台のカメラが接続されている場合には、信号源方向又は位置推定装置の出力をモニタすべきカメラを逐次切り替えるようにしてもよい。または、信号源の位置推定とともに、信号源に向かってカメラをズームインさせるようにしてもよい。
【０１５９】
（適用例２）
信号強調装置の入力装置への適用：
本発明に係る信号強調装置を、コンピュータや携帯端末（ＰＤＡ）などの情報処理装置におけるキーボードその他のユーザ入力装置の代替とすることができる。すなわち、複数のマイクを含んだ信号強調装置の出力を音声認識装置で音声認識処理することによって（図１５を参照のこと）、情報処理装置に対するキャラクタ形式のコマンド入力とすることができる。また、音声ベースでのユーザインターフェースにすることもできる。
例えば、話者ごとの発言内容を記録装置に記録することによって議事録作成システムを構築することができる。
【０１６０】
（適用例３）
信号強調装置のコマンダ／コントローラへの適用：
本発明に係る信号強調装置を、ロボット、自動車、コンピュータ、携帯端末（ＰＤＡ）、各種の情報家電機器に対して、ユーザ操作を行うコマンダ／コントローラとして適用することができる。すなわち、複数のマイクを含んだ信号強調装置の出力を音声認識装置で音声認識処理してキャラクタに変換し、さらに入力キャラクタをコマンド解析装置により解析処理することによって（図１６を参照のこと）、ユーザのコマンド内容を特定することができる。
例えば、このようなコマンダ／コントローラを介して、ロボットに対して「歩行」やその他の動作指示を行うことができる。また、コンピュータに対して、「ワープロ」と話しかけることにより、該当するアプリケーションを起動させたり、さらにキャラクタ入力を行うことができる。また、ＰＤＡのような個人情報を扱う携帯端末に対して、「今日の予定は？」と話しかけることによりスケジュール帳を起動して該当する日付や時刻の予定を引き出すことができる。また、電話に対して、「○△さんに電話をかけて」と話しかけることにより、電話帳を検索して該当する人物に対して発呼することができる。また、情報家電に対して、「部屋の温度を上げて」とか「××を録画して」などの操作を指示することができる。また、カーナビゲーション・システムに対して「○△宅まで行きたい」と話しかけることにより、道案内を起動することができる。
【０１６１】
（適用例４）
信号強調装置、並びに信号源方向又は位置推定装置のＴＶ会議システムへの適用：
複数のマイクを含んだ信号強調装置の出力を音声認識装置に入力し、さらに音声入力結果を記録・表示装置に出力するとともに、複数のマイクを含んだ信号源方向又は位置推定装置による信号源の方向又は位置の推定結果を記録・表示装置に出力する（図１７を参照のこと）。
【０１６２】
記録・表示装置では、何処の席にいる人がどのような発言を行ったかを逐次的に記録したり、現在発言中の人を特定して、カメラで追跡して画面表示したりすることができる。
【０１６３】
（適用例５）
スマート・アンテナへの適用：
複数のアンテナを含んだ信号方向又は位置推定装置の出力をアンテナ指向性制御装置に出力するとともに、このアンテナ指向性制御装置の制御指令を基にアンテナの位置や向きなどを調整することができる（図１８を参照のこと）。このようなシステムによれば、アンテナの受信感度を自律的且つダイナミックに調整することができる。また、基地局などに配備することによって、携帯電話からの接続状況を向上させることができる。
【０１６４】
［追補］
以上、特定の実施例を参照しながら、本発明について詳解してきた。しかしながら、本発明の要旨を逸脱しない範囲で当業者が該実施例の修正や代用を成し得ることは自明である。
【０１６５】
上記の第１の実施例においては、相互相関関数から信号源の方向又は位置を推定する際に、また、第２の実施例においては、相互相関関数と周波数帯域別信号とから信号強調する際に、周波数帯域ごとに相互相関関数を規格化し、受信装置の組み合わせごとにこれらの規格化した相互相関関数を、すべてあるいは一部の複数の周波数帯域にわたって加算し、これらの加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求め、これらを候補遅延時間としている。但し、上記の規格化を必ずしも必要とする訳ではなく、規格化を行わずに以降の処理を行ってもよい。また、相互相関関数をすべてあるいは一部の複数の周波数帯域にわたって加算することなく、それぞれの相互相関関数が極大となる遅延時間を、受信装置の組み合わせごとに求め、これらを候補遅延時間としてもよいし、あるいは複数の異なる周波数帯域を加算した、複数の加算された相互相関関数が極大となる遅延時間を受信装置の組み合わせごとに求め、これらを候補遅延時間としてもよい。
【０１６６】
また、上記の第１及び第２の各実施例においては、信号を音として、信号源は音源、受信装置はマイクロホンとしているが、信号は他のもの、例えば光、電磁波、水面波などであっても本発明を好適に適用することができる。また、これらの信号の場合は、信号源はそれぞれの信号を発する物又は装置であり、受信装置はそれぞれの信号を受信可能な装置を用いた構成とすればよい。
【０１６７】
また、上記の第１及び第２の各実施例においては、ディジタル信号処理技術を用いて実現しているが、信号源の種類、用いることのできるハードウェアなどの条件に応じて、その一部あるいは全部をアナログ信号処理技術や、計算機による処理技術によって置き換えることができることは言うまでもない。
【０１６８】
要するに、例示という形態で本発明を開示してきたのであり、限定的に解釈されるべきではない。本発明の要旨を判断するためには、冒頭に記載した特許請求の範囲の欄を参酌すべきである。
【０１６９】
【発明の効果】
以上詳記したように、本発明によれば、存在する方向や位置があらかじめ知られていない１あるいは複数の信号源からの信号を複数の受信装置で受信した信号を用いて、信号源の方向又は位置を推定することができる、優れた信号源方向又は位置推定方法及び装置を提供することができる。
【０１７０】
また、本発明によれば、存在する方向や位置があらかじめ知られていない１あるいは複数の信号源からの信号を、複数の受信装置で受信した信号を用いて強調することができる、優れた信号強調方法及び装置を提供することができる。
【０１７１】
本発明によれば、１あるいは複数の信号源からの信号を複数の受信装置で受信した各受信信号を、複数の帯域通過フィルタによりそれぞれ複数の周波数帯域別信号に分解して、これら周波数帯域別信号間の相互相関関数を相当する周波数帯域ごとに受信装置の組合せごとに計算して、これらの相互相関関数を基に信号源の方向又は位置を推定するようにしたので、雑音の多い環境や複数の信号源が同時に信号を発生している場合においても信号源の方向又は位置を充分に推定することができる。
【０１７２】
また、本発明では、１あるいは複数の信号源からの信号を複数の受信装置で受信した各受信信号を、複数の帯域通過フィルタによりそれぞれ複数の周波数帯域別信号に分解し、これら周波数帯域別信号間の相互相関関数を相当する周波数帯域ごとに受信装置の組合せごとに計算して、これらの相互相関関数を基に信号源の方向又は位置と受信装置の配置とによって定まる各受信装置ごとの遅延時間を推定し、この推定遅延時間を用いて各受信信号の周波数帯域別信号をそれぞれ遅延し、これらの遅延した周波数帯域別信号を加算して、信号源の信号を強調するようにしたので、雑音の多い環境や複数の信号源が同時に信号を発生しているような場合においても、１あるいは複数の信号源ごとの信号を強調することができる。
【図面の簡単な説明】
【図１】本発明の第１の実施例に係る信号処理装置の構成を模式的に示したブロック図である。
【図２】帯域通過フィルタ・バンク部５−ｍの内部構成を示したブロック図である。
【図３】相互相関関数計算部６の内部構成を示したブロック図である。
【図４】マイクロホンの個数Ｍが２の場合における信号源方向又は位置推定部７の内部構成（第１の実施例）を示したブロック図である。
【図５】信号源方向推定に用いられる式［数１］を解説するための図である。
【図６】マイクロホンの個数Ｍが３以上の場合における信号源方向又は位置推定部７の内部構成（第１の実施例）を示したブロック図である。
【図７】信号源位置推定部の内部構成を示したブロック図である。
【図８】［数４］を説明するための図である。
【図９】本発明の第２の実施例に係る信号処理装置の構成を模式的に示したブロック図である。
【図１０】マイクロホンの個数Ｍが２の場合における信号源方向又は信号強調部２０の内部構成（第２の実施例）を示したブロック図である。
【図１１】マイクロホンの個数Ｍが３以上の場合における信号強調部２０の内部構成（第２の実施例）を示したブロック図である。
【図１２】強調信号生成部２２の内部構成を示したブロック図である。
【図１３】強調信号計算の方法の一例を説明するための図である。
【図１４】信号源方向又は位置推定装置をＴＶ会議システムや監視システムに適用した例を説明するための図である。
【図１５】信号強調装置を入力装置に適用した例を説明するための図である。
【図１６】信号強調装置をコマンダ／コントローラに適用した例を説明するための図である。
【図１７】信号源方向又は位置推定装置と信号強調装置をＴＶ会議システムに適用した例を説明するための図である。
【図１８】信号源方向または位置推定装置をスマート・アンテナに適用した例を説明するための図である。
【符号の説明】
１…音源，２…マイクロホン
３…増幅器，４…Ａ／Ｄ変換器
５…帯域通過フィルタ・バンク
６…相互相関関数計算部
７…信号源方向又は位置推定部
８…帯域通過フィルタ
９，１０…遅延器
１１…積算器，１２…重み付き平均低域通過フィルタ
１３…規格化計算部，１４…加算部
１５…極大検出部，１６…信号源方向推定部
１７…信号源位置推定部
１８…遅延時間分析部，１９…信号源位置計算部
２０…信号強調部
２１，２２…強調信号生成部
２３…強調信号計算部
【図１】
0004815661.tif 000011
【図２】
0004815661.tif 000012
【図３】
0004815661.tif 000013
【図４】
0004815661.tif 000014
【図５】
0004815661.tif 000015
【図６】
0004815661.tif 000016
【図７】
0004815661.tif 000017
【図８】
0004815661.tif 000018
【図９】
0004815661.tif 000019
【図１０】
0004815661.tif 000020
【図１１】
0004815661.tif 000021
【図１２】
0004815661.tif 000022
【図１３】
0004815661.tif 000023
【図１４】
0004815661.tif 000024
【図１５】
0004815661.tif 000025
【図１６】
0004815661.tif 000026
【図１７】
0004815661.tif 000027
【図１８】
0004815661.tif 000028