(57)【要約】  最適化されたセンサーアレイ用のモーダルビームフォーマが開示されている。ビームフォーマのビームパターンを形成する方法では、 ビームフォーマにおいて、センサーアレイから入力信号が受信され、これらの入力信号が球面調和関数ドメインに分解され、 球面調和関数に重み付け係数が付与され、 これらを組み合わせて出力信号が形成され、また、重み付け係数が、ある与えられた一組の入力パラメータに対して凸最適化演算を用いて最適化される。多重主ローブ生成、均一および不均一副ローブ制御、自動ヌルステアリング、ロバストネスおよびホワイトノイズ利得に対する二次錐計画法制約を形成するための公式が提供されている。【選択図】図１８
2012523731.tif 000002
【特許請求の範囲】
【請求項１】
ビームフォーマのビームパターンを形成する方法であって、
前記ビームフォーマが、
センサーアレイから入力信号を受信し、
前記入力信号を球面調和関数ドメインに分解し、
前記球面調和関数に重み付け係数を付与し、
これらを組み合わせて出力信号を形成し、
これらの重み付け係数が、凸最適化演算によって、ある与えられた一組の入力パラメータに対して最適化される、方法。
【請求項２】
前記センサーアレイが、センサー位置が概念上の球面に配置される球面アレイである、請求項１に記載の方法。
【請求項３】
前記センサーアレイが、開球アレイ、剛球アレイ、半球アレイ、２重開球アレイ、球殻アレイ、およびカーディオイドマイクロホンを備えた単一開球アレイからなる群から選択される形態を有している、請求項２に記載の方法。
【請求項４】
前記センサーアレイが、音声帯域用途に設計され、約８ｃｍ〜約３０ｃｍの最大寸法を有している、請求項１，２または３に記載の方法。
【請求項５】
前記センサーアレイがマイクロホンアレイである、請求項１乃至４のうちのいずれか一項に記載の方法。
【請求項６】
最適化問題、および任意選択的に制約が、前記センサーアレイの出力を最小限に抑えること、副ローブレベルを最小限に抑えること、主ローブ領域のひずみを最小限に抑えること、およびホワイトノイズ利得を最大化することのうちの１つ以上として定式化されている、請求項１乃至５のうちのいずれか一項に記載の方法。
【請求項７】
前記最適化問題が、前記センサーアレイの出力を最小限に抑えることとして定式化されている、請求項１乃至６のうちのいずれか一項に記載の方法。
【請求項８】
前記一組の入力パラメータが、前記ビームパターンに主ローブを形成するように、指定された方向の前記アレイ利得をある与えられたレベルに維持するという要件を含んでいる、請求項１乃至７のうちのいずれか一項に記載の方法。
【請求項９】
前記一組の入力パラメータが、前記ビームパターンに複数の主ローブを形成するように、複数の指定された方向の前記アレイ利得をある与えられたレベルに維持するという要件を含んでいる、請求項８に記載の方法。
【請求項１０】
前記ビームパターンに異なるレベルの複数の主ローブを形成するように、前記複数の指定された方向に対して必要な利得レベルが提供される、請求項９に記載の方法。
【請求項１１】
前記ビームフォーマが、凸制約として前記要件または各要件を定式化する、請求項８、９または１０に記載の方法。
【請求項１２】
前記ビームフォーマが、線形等式制約として前記要件または各要件を定式化する、請求項１１に記載の方法。
【請求項１３】
前記ビームフォーマが、前記指定された方向から前記センサーアレイに入射する単位強度平面波に対する前記アレイ出力が前もって決められている定数と等しいという要件として、前記要件または各要件を定式化する、請求項１２に記載の方法。
【請求項１４】
前記一組の入力パラメータが、前記ビームパターンにヌルを形成するように、指定された方向の前記アレイ利得をある与えられたレベル未満にするという要求を含んでいる、請求項１乃至１３のうちのいずれか一項に記載の方法。
【請求項１５】
前記一組の入力パラメータが、前記ビームパターンに複数のヌルを形成するように、複数の指定された方向の前記アレイ利得をある与えられたレベル未満にするという要件を含んでいる、請求項１４に記載の方法。
【請求項１６】
前記ビームパターンに異なる深さの複数のヌルを形成するように、前記複数の指定された方向の各々に対して最大利得レベルが提供される、請求項１５に記載の方法。
【請求項１７】
前記ビームフォーマが、凸制約として前記要件または各要件を定式化する、請求項１４、１５または１６に記載の方法。
【請求項１８】
前記ビームフォーマが、二次錘制約として前記要件または各要件を定式化する、請求項１７に記載の方法。
【請求項１９】
前記ビームフォーマが、前記指定された方向から前記センサーアレイに入射する単位強度平面波に対する前記アレイ出力の強度が前もって決められている定数未満であるという要件として、前記要件または各要件を定式化する、請求項１８に記載の方法。
【請求項２０】
前記一組の入力パラメータが、指定されたレベルのロバストネスを前記ビームパターンが有しているという要求を含んでいる、請求項１乃至１９のうちのいずれか一項に記載の方法。
【請求項２１】
前記指定されたレベルのロバストネスが、前記重み付け係数を含むベクトルのノルムに対する制限として指定される、請求項２０に記載の方法。
【請求項２２】
前記ノルムがユークリッドノルムである、請求項２１に記載の方法。
【請求項２３】
前記重み付け係数が、二次錐計画法によって最適化される、請求項１乃至２２のうちのいずれか一項に記載の方法。
【請求項２４】
位数ｎの各球面調和関数に対して１つ以上の重み付け係数が最適化されるが、各位数の球面調和関数内で、前記１つ以上の重み付け係数が、前記位数ｎのすべての次数ｍ＝−ｎ〜次数ｍ＝ｎまでにおいて共通である、請求項１乃至２３のうちのいずれか一項に記載の方法。
【請求項２５】
前記一組の入力信号が、前記球面調和関数ドメインへ分解される前に、周波数ドメインに変換される、請求項１乃至２４のうちのいずれか一項に記載の方法。
【請求項２６】
前記ビームフォーマが、周波数ドメイン信号が狭帯域周波数ビンに分割される広帯域ビームフォーマであり、各狭帯域周波数ビンが別々に最適化さ、重み付けされてから、これらの狭帯域周波数ビンが再結合されて広帯域出力を形成する、請求項２５に記載の方法。
【請求項２７】
前記一組の入力信号が、時間ドメインで処理され、また、前記重み付け係数が、前記球面調和関数信号に適用される有限インパルス応答フィルタのタップ重みである、請求項１乃至２４のうちのいずれか一項に記載の方法。
【請求項２８】
ビームフォーマであって、
各センサーが信号を生成するように構成されたアレイ状に並べられているセンサーと、
入力信号を球面調和関数ドメインへ分解し、分解された前記信号を出力するように構成された球面調和関数分解手段と、
１組の入力パラメータに基づいて、前記分解された信号に付与される重み付け係数を凸最適化演算を用いて計算するように構成された重み付け係数計算手段と、
計算された前記重み付け係数と前記分解された信号を組み合わせて出力信号を形成する出力形成手段とを備えてなる、ビームフォーマ。
【請求項２９】
信号追跡手段をさらに備えており、該信号追跡手段が、前記センサーからの前記信号を評価して所望の信号ソースの方向および不要な干渉ソースの方向を判断するように構成されてなる、請求項２８に記載のビームフォーマ。
【請求項３０】
ビームフォーマのビームパターンを形成する方法であって、
前記ビームフォーマが、
センサーアレイから入力信号を受信し、
前記信号に重み付け係数を付与し、
これらを組み合わせて出力信号を形成し、
前記重み付け係数が、前記ビームパターンに複数の主ローブを形成するように、複数の指定された方向のアレイ利得をある与えられたレベルに維持するという制約条件下で、ある与えられた一組の入力パラメータに対して凸最適化演算によって最適化され、
各要件が、指定された方向から前記センサーアレイに入射する単位強度平面波に対する前記センーアレイの出力が前もって決められている定数と等しいという要件として定式化される、方法。
【請求項３１】
コンピュータ上で実行されると、請求項１乃至２７または３０のうちのいずれか一項に記載のステップを実行するように構成されてなる、ソフトウェア製品。
【請求項３２】
前記ソフトウェア製品がデータキャリアである、請求項３１に記載のソフトウェア製品。
【請求項３３】
前記ソフトウェア製品が、リモート位置から送信される信号を含む、請求項３１に記載のソフトウェア製品。
【請求項３４】
物理的キャリアの形態を有しているソフトウェア製品を製造する方法であって、前記データキャリアにインストラクションを格納することを含み、該インストラクションがコンピュータにより実行されると、該コンピュータが請求項１乃至２７または３０のうちのいずれか一項に記載の方法を実行する、方法。
【請求項３５】
ソフトウェア製品をリモート位置へ、該リモート位置のコンピュータへデータを送信することによって送る方法であって、前記データがインストラクションを含み、該インストラクションが前記コンピュータによって実行されると、該コンピュータが請求項１乃至２７または３０のうちのいずれか一項に記載の方法を実行する、方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はビームフォーミングに関するものである。
【背景技術】
【０００２】
ビームフォーミングとは、アレイ状に配置された複数のセンサーからの入力を組み合わせるための技術のことである。アレイ状に配置されたセンサーの各々は、その場所に応じて異なる信号を生成し、これらの信号はシーン全体を表すものである。これらの信号を異なる方法で組み合わせることによって、たとえば各受信信号に対し異なる重み付け係数を付与するまたは異なるフィルタを使用することによって、シーンの異なる側面を強調するおよび／または抑制することができる。具体的にいえば、アレイの指向性については、ある方向に対する重み付けを大きくすることにより選択された方向に対してアレイの感度をより高くすることができる。
【０００３】
ビームフォーミングは、電磁波および音波に適用することができ、たとえばレーダおよびソナーに用いられている。センサーアレイは、用途および使用される波長に応じて実質的にいかなるサイズまたは形状をとることもできる。簡単な用途では、１次元の線形アレイで十分である場合もある。複雑な用途では、２次元または３次元のアレイが必要となる場合もある。最近では、ビームフォーミングは、３次元（３Ｄ）での音響受信、室内音響学のための音場解析、ビデオ・テレコンファレンスでの音声受信、到来方向推定ならびに騒音制御の用途の分野で用いられている。これらの用途の場合、十分な３Ｄ音響解析を実現するために、３次元のアレイ状に並べられた複数のマイクロホンが必要となる。
【０００４】
３次元アレイの可能な配置のうち、球面アレイがとくに興味深い。というのは、他の標準的なアレイの幾何学構造に比べてより柔軟な３次元ビームパターンの合成が可能となるとともに、球面調和関数ドメインの数学的フレームワークを用いてアレイ処理を行うことができるからである。通常、球面アレイは、その表面に複数のセンサーが分配されている球体の形態をとっている。最も一般的な実施形態としては、複数のセンサーが物理的な球体の表面上に配置されている「剛球」、および表面が単に概念上のものでしかなく、複数のセンサーがこの概念上の表面の位置に他の手段によって保持されている「開球」が含まれている。２重開球（一方が他方の内側に位置する同心円状に並べられた概念上の２つの球面にセンサーが配置されている）、球殻アレイ（概念上の２つの球面の間に、すなわち２つの球面により形成されるシェル内にセンサーが配置されている）、カーディオイドマイクロホンを備えた単一開球および半球の如き他の構成も適切な実施形態である。これらはすべて、音場を分解して球面調和関数を形成するために用いることができる。
【０００５】
与えられたアレイ（たとえば、音響用途の場合、マイクロホンもしくはハイドロホン、または無線通信用途の場合、アンテナからなるアレイ）の場合、アレイ内のセンサーの各々に付与される重み付けにより、そのアレイの「ビームパターン」が定義されることになる。しかしながら、通常、アレイの１つ以上の部分が他の部分よりも大きな重み付けがなされる場合、ビームパターンには、受信強度が高く、信号利得が高い領域を示す「ローブ」と、入射波が著しく減衰される、受信強度の低い領域を示す「ヌル」とが形成される。ローブおよびヌルの配置は、センサーに付与された重みおよびセンサーの物理的な配置に依存する。しかしながら、通常、ビームパターンには、最も強く信号を受信する方向の「主」ローブ（すなわち、ビームパターンの一次マキシマム）およびビームパターンの二次（または、他の）マキシマムである１つ以上の「副」ローブが含まれる。ローブとローブとの間にはヌルが形成される。
【０００６】
音響用途では、聴覚情景解析のこと考えれば、この場合の問題は、あるソース（たとえば、あなたに話しかけている友達）の話を聞こうとするとともに、ある干渉ソース（たとえば、あなたの隣でなされている他の会話）からの音を無視するまたは遮ろうとすることが望まれるカクテルパーティ問題にたとえることができる。それと同時に、一般的にパーティの暗騒音を無視するまたは遮ることも望まれる。同様に、マイクロホンアレイにおけるビームフォーミング問題は、所望のソースにアレイの受信力を集中させるとともに、干渉ソースおよび暗騒音の影響を最小限に抑えることである。
【０００７】
これらの問題は、２つの部屋がマイクロホンアレイおよびラウドスピーカにより通信可能にリンクされている、すなわち各部屋が音を拾い上げて他の部屋へ音声信号として送信するマイクロホンアレイと、他の部屋から受信される信号を変換して音を発生するラウドスピーカとを備えているテレコンファレンスの如き用途において非常に重要である。任意の与えられた時間において、一方の部屋（近い方）には、音声を拾い上げなければならない一人以上の話者、通信の相手（遠い方）からの音を生じるラウドスピーカの如き理想的には遮断されなければならない干渉ソース、ならびにたとえば空調機のノイズまたはエコーおよび話者および／またはラウドスピーカに起因する反響音のような暗騒音が存在しうる。
【０００８】
一般的に、この問題は、「ビームステアリング」として知られているプロセスによって解決しようとされている。「ビームステアリング」では、ビームパターンの主ローブが目的信号の方向に向けられるとともに、ビームパターンのヌル（ノッチとしても知られている）が、干渉信号の方向に向けられる（「ヌルステアリング」）ようになっている。
【０００９】
一般的に、副ローブは、所望の信号より強いものを受信するビームパターンの領域を表わしている。すなわち副ローブは、ビームパターンの望まれないローカルマキシマムのことである。副ローブは避けられないが、重み付け係数の適切な選択によって、副ローブのサイズを調製することができる。
【００１０】
さらに、１を超える目的信号方向がある場合、ビームパターンに複数の主ローブを形成することが可能である。調製することが望ましいビームパターンの他の特徴には、主ローブのビーム幅、ロバストネス、すなわち異常な入力または予期しない入力に持ちこたえるシステムの能力、およびアレイ信号利得（すなわち、信号対ノイズ比（ＳＮＲ）の利得）が含まれる。
【００１１】
ほとんどの環境では、聴覚情景は常に変化している。目的の信号が行き来し、干渉ソースからの信号が行き来し、信号の方向が変化し、振幅ノイズレベルが高くなる。これらの状況では、理想的には変化する状況にセンサーアレイが適応できる必要がある。たとえば、センサーアレイは、移動する目的信号を追尾するようにビームパターンの主ローブを移動させる必要がある。または、センサーアレイは、新規の干渉ソースを打ち消すための新規のヌルを形成する必要がある。同様に、干渉ソースが消滅した場合、システムの制約を変更してさらに最適な解を得ることが可能となる。したがって、これらの状況では、センサーアレイは適応性を有したものである必要がある。すなわち、センサーアレイは、制約を再評価し、最適化問題を解きなおし、新しい最適解を見つけることができるようになしてある必要がある。さらに、聴覚情景が急速に変わるような状況、たとえばテレコンファレンスのような場合、その間ずっと人々が会話を始めたり止めたりしており、目的信号ソースおよび干渉ソースが絶えず数と方向において変化しているので、理想的には、ビームフォーマは実時間で動作可能である必要がある。
【００１２】
この分野では、複数の研究が成されてきている。いくつかの具体例を挙げると、マイヤーおよびエルコは、下記の文献において、ルック方向に対して対称となっており、ビームパターンの形状を変更することなく３Ｄスペースで向きを変えることが可能な球面マイクロホンアレイビームパターン設計を用いた音場球面調和関数分解の応用および解析について報告している。
［ジェー・マイヤー（Ｊ．Ｍｅｙｅｒ）およびジー・エルコ（Ｇ．Ｅｉｋｏ）、「音場の正規直交分解に基づいた高度にスケーリング可能な球面マイクロホンアレイ（Ａ  ｈｉｇｈｌｙ  ｓｃａｌａｂｌｅ  ｓｐｈｅｒｉｃａｌ  ｍｉｃｒｏｐｈｏｎｅ  ａｒｒａｙ  ｂａｓｅｄ  ｏｎ  ａｎ  ｏｒｔｈｏｎｏｍａｌ  ｄｅｃｏｍｐｏｓｉｔｉｏｎ  ｏｆ  ｔｈｅ  ｓｏｕｎｄｆｉｅｌｄ）」、ＩＣＡＳＳＰ会報、第２巻、ページ１７８１〜１７８４、２００２年５月］
また、ＷＯ２００６／１１０２３０をさらに参照されたい。これらの研究に対する展開として、ラファリーは、下記の文献において、一般的に用いられている遅延和型ビームパターン設計法を球面マイクロホンアレイに適用した。
［ビー・ラファリー（Ｂ．Ｒａｆａｅｌｙ）、「位相モード対遅延和型の球面マイクロホンアレイ処理（Ｐｈａｓｅ−ｍｏｄｅ  ｖｅｒｓｕｓ  ｄｅｌａｙ−ａｎｄ−ｓｕｍ  ｓｐｈｅｒｉｃａｌｍｉｃｒｏｐｈｏｎｅ  ａｒｒａｙ  ｐｒｏｃｅｓｓｉｎｇ）」、ＩＥＥＥ信号プロセスレター（ＩＥＥＥ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓ．Ｌｅｔｔ．）、第１２巻、１０号、ページ７１３〜７１６、２００５年１０月］
すなわちアレイ重みを用いて、単一の平面波に起因する自由音場マイクロホンでの遅れを補償した。このアプローチでは、高いロバストネスが得られたが、低周波数での指向性の減少が犠牲になった。他の研究では、ラファリーらは、下記の文献において、音場の指向性の解析を向上させるために古典的なドルフ・チェビシェフ（Ｄｏｌｐｈ−Ｃｈｅｂｙｓｈｅｖ）パターン設計アプローチ用いて、与えられた主ローブ幅およびアレイ位数（ａｒｒａｙ  ｏｒｄｅｒ）において副ローブの調製も達成した。
［ビー・ラファリー（Ｂ．Ｒａｆａｅｌｙ）、エー・コレッツ（Ａ．Ｋｏｒｅｔｚ）、アール・ウィニック（Ｒ．Ｗｉｎｉｋ）、エム・アグモン（Ｍ．Ａｇｍｏｎ）、「室内音響解析を向上させるための球面マイクロホンアレイビームパターン設計」、室内音響学国際シンポジウム会報、ページＳ４２、２００７年９月］
リーおよびデュラスワミは、下記の文献において、ビームパターン合成にホワイトノイズ利得（ＷＮＧ）制約を課することにより、実用化に役立つビームフォーミング指向性とロバストネスとの間のバランスを見出すためのアレイ重み付け最適化方法を報告している。
［ジー・ワイ・リー（Ｚ．Ｙ．Ｌｉ）およびアール・デュラスワミ（Ｒ．Ｄｕｒａｉｓｗａｍｉ）、「ビームフォーミングのための球面マイクロホンアレイの柔軟で最適なデザイン」、オーディオ・スピーチ・言語プロセスＩＥＥＥ会議録、第１５巻、２号、ページ７０２〜７１４、２００７年２月］
上述の研究では対称なビームパターンだけが想定されているが、ラファリーは、下記の文献において、ビームパターン設計法を球面マイクロホンアレイの非対称なケースまで適用範囲を広げた。
［ビー・ラファリー、「指向性のあるルームインパルス応答の解析のための複数のヌルを備えた球面マイクロホンアレイ」、ＩＣＡＳＳＰ会報、ページ２８１〜２８４、２００８年４月］
このアプローチは、スペースドメインおよび球面調和関数ドメインの両方で定式化されていることに加えて、多重ヌルステアリング手法を含んでいる。多重ヌルステアリング手法では、固定ヌルが、ビームパターンに形成され、信号対ノイズ比を向上させるために、既知の外部ビーム方向から入ってくる干渉に向けて方向付けされるようになっている。
【００１３】
アージェンチエリら（Ａｒｇｅｎｔｉｅｒ  ｅｔ  ａｌ)による以下の文献では、凸最適化技法が用いられ、その問題が球面調和関数フレームワークを用いて解析されたが、波動場は球面調和関数へ分解されなかった。
［「ロボット工学における近接場または遠距離場のスピーカ配置のためのモーダル解析に基づくビームフォーミング」、２００６年ＩＥＥＥ／ＲＳＪ知能ロボットおよびシステムに関する国際会議会報、ページ８６６〜８７１］
しかしながら、球面調和関数ドメインビームフォーミングに関する上述の研究では、任意の外部ビーム方向から入ってくる動的な干渉を抑制するために、複数の深いヌルが、それに適応して、ビームパターン内に形成され、方向付けできるようにはなっていなかった。このような干渉抑制は、ビデオコンファレンスまたはテレコンファレンスの用途におけるスピーチ強化および多重チャネル音響エコーキャンセレーションにおいて、ならびに指向性ルームインパルス応答解析（すなわち、インパルスの生成および反射解析による部屋の音響解析）において望まれることが多い。これに加えて、上述の研究は、副ローブ制御制約およびロバストネス制約の如き多重ビームフォーミングの性能パラメータを単一の最適化アルゴリズムに有効に含めることができなかったので、これらの相互に関連するパラメータすべてに対する大域的（グローバル）最適解を得ることは今日までできなかった。
【００１４】
主な困難は、最適化アルゴリズムが集中的に多数の計算をすることを必要としているということである。テレコンファレンスのような上述の用途が消費者向けの用途であるので、アルゴリズムは、容易に入手可能な消費者レベルの計算能力で合理的な時間内に実行可能でなければならない。さらに特筆すべきことは、これらの用途が実時間に基づくものであるので、実時間において適応性があることが必要となることである。したがって、実時間動作を維持しながら、所望のパラメータすべてを最適化するのは非常に困難なことである。実時間動作のための要件はセンサーアレイの用途に応じて変わりうる。しかしながら、テレコンファレンスのような音声を拾い上げる用途では、センサーアレイは、聴覚情景の動的変化に同一の速度で適合して行くことができなければならない。人々が一度に数秒間話す傾向があるので、ビームパターンを再最適化するために数秒（最大約５秒）しか必要とならないビームフォーマが有用である。しかしながら、好ましくは、いわれたことをなにも逃さないように、システムが、一秒程度の時間スケールでビームパターンを再最適化する（すなわち、最適な重み付けを再計算する）ことができる必要がある。最も好ましくは、システムは、新規な信号ソース（たとえば、新規なスピーカ）が検出されると可及的速やかに、ビームフォーマがその方向に対する適切なアレイ利得を提供することを担保するように、一秒間に数回重み付けを再最適化することができる必要がある。
【００１５】
いうまでもなく、ムーアの法則によれば計算能力がいまだ指数関数的に増大していっているので、計算能力の進歩が必要な計算を行なう時間を急速に減少させており、将来的には、実時間に基づく用途では、再最適化が著しく速い速度で実行されるであろうことが期待される。
【００１６】
与えられたシナリオにおいてビームパターンの選択に影響を与えるパラメータが複数存在するので、これらのパラメータのうちの１つに対する最適解が必ずしもその他のパラメータにとって最適なものであるとは限らない。したがって、それらの間で妥協点を見出さなければならない。これらの要因の間で最良の（最適な）妥協点を見つけることは、システムの要件に左右される。これらは、最適化問題における制約として定式化することができる。たとえば、システムがある指向性を有することが必要とされる場合もあれば、または選択されたしきい値を超える利得を得ることを必要とされる場合もある。これに代えて、副ローブがあるしきい値未満であることを必要とされる場合もあれば、システムがあるロバストネスを有していることを必要とされる場合もある。上述のように、最適化は集中的に多数の計算をするプロセスであり、制約を加えるたびに、集中的にさらに多くの数の計算をすることになる。したがって、実際問題として、最適解を合理的な時間内に見つけようとする場合、１以上の制約をシステムに課すことは通常実行不可能なことである。
【発明の概要】
【発明が解決しようとする課題】
【００１７】
今日まで行なわれてきた研究では、最適化アルゴリズムは１つまたは２つの制約だけに制限されていた。場合によっては、複数の制約がそれぞれ別個に１つずつ個々のステージで解かれるようになっていたが、大域的な最適解を得ることはできなかった。
【００１８】
複数の制約をシステムに課して、球面アレイにとっての大域的に最適なビームパターンを見出す方法を提供する必要性が依然としてある。
【課題を解決するための手段】
【００１９】
本発明の第一の態様によれば、ビームフォーマのビームパターンを形成する方法では、  ビームフォーマが、センサーアレイから入力信号を受信し、これらの入力信号を球面調和関数ドメインに分解し、球面調和関数に重み付け係数を付与し、これらを組み合わせて出力信号を形成し、これらの重み付け係数が、ある与えられた一組の入力パラメータに対して凸最適化演算を用いて最適化される。
【００２０】
目的関数および制約を凸関数として表すことによって、凸最適化技術を用いることが可能となる。凸最適化は、グローバルミニマムが存在するならばグローバルミニマムを見つけることができ、さらに、数値法を用いてグローバルミニマムを迅速にかつ効率的に見つけることができることを担保するという利点を有している。
【００２１】
従来の研究では、周波数に依存しない規則的なまたは不規則なビームパターンを容易に形成するため、アレイ重み付け設計アプローチでは、周波数に依存するコンポーネントを分離するために、球面調和関数ドメインにおいて（下記に詳細に説明されている）モード振幅ｂn（ｋａ）反転が用いられる。しかしながら、ｂn（ｋａ）は、あるｋａ値およびｎ値が小さな値であり、その反転は、実際の実施にあたって、ビームフォーマのロバストネスを劣化させる恐れがある。本発明では、直接的により一般的な重みを最適化フレームのターゲットとすることにより、最適化問題を、凸最適化問題として、すなわち目的関数および制約がすべて凸関数である最滴化問題として定式化することができるようになる。 上述のように、凸最適化演算の利点は、最適化変数の最適値を素早く見つけることができる高速（すなわち、高速計算可能な）数値ソルバーが存在するということにある。さらに、上述のように、凸最適化演算の結果は、常に局地的（ローカル）な最適解ではなく大域的（グローバル）な最適解をもたらす。したがって、上述のように定式化することにより、本発明にかかるビームフォーマは、複数の制約が課された場合であっても、実時間で、アレイビームパターンを適応的に最適化することができる。
【００２２】
凸最適化技術は従来から知られているものである。凸最適化問題を解くためのさまざまな数値法およびソフトウェアツールも従来から知られている。しかしながら、凸最適化は、目的関数および最適化する制約がすべて凸関数である場合にしか用いることができない。ある関数ｆが凸関数であるためには次の条件が満たされたければならない。すなわち、すべてのｘ、ｙおよびａ、ｂにおいて、ｆ（ａｘ＋ｂｙ）≦ａｆ（ｘ）＋ｂｆ（ｙ）である。この式で、ａ＋ｂ＝１であり、ａ≧０であり、ｂ≧０である。したがって、与えられた最適化問題を凸最適化技術を用いて解くことは必ず可能であるとは限らない。まず、最適化問題は、凸最適化演算を適用することができるように定式化されなければならない。換言すれば、最小化する必要のあるシステム特性を考え、それを凸関数として定式化しなければならない。さらに、最適化問題における制約すべてを、凸の等式／不等式または線形等式として定式化しなければならない。ビームフォーミング問題を凸最適化演算問題として定式化することによって、本発明において、複数の制約が課されたビームフォーミング問題の実時間の解を高速で計算することを可能とする複数の非常に効率的なアルゴリズムを用いることができるようになる。
【００２３】
好ましくは、センサーアレイは、センサーが概念上の球面上に配置されている球面アレイである。このような配置が対称であると処理が簡単なものとなる。本発明では、複数の異なる球面センサーアレイ配置を用いることが可能である。好ましくは、センサーアレイは、開球アレイ、剛球アレイ、半球アレイ、２重開球アレイ、球殻アレイおよびカーディオイドマイクロホンを備えた単一開球アレイからなる群から選択される１形態を有している。
【００２４】
アレイサイズは、用途および用いられる波長に応じて大きく変わりうる。しかしながら、音声を拾い上げる用途で用いられるマイクロホンアレイの場合、センサーアレイは、約８ｃｍと約３０ｃｍとの間の大きな寸法を有していることが好ましい。球面アレイの場合、最大寸法は直径のことである。大きな球体は、低周波をうまく扱うことができる利点を有しているが、高周波での空間エイリアシングを回避するためには、２つのマイクロホンの間の距離が最高周波数の波長の半分よりも小さい必要がある。したがって、マイクロホンの数が有限である場合、球体が小さいということは、マイクロホンとマイクロホンとの間の距離が短く、空間エイリアシングの問題が少ないということを意味する。いうまでもなく、５〜１００ＭＨｚの周波数を想定しうる超音波イメージングの如き高周波用途では、センサーアレイサイズは著しく小さくなる。同様に、ソナー用途では、アレイサイズは著しく大きくなる場合もある。
【００２５】
好ましくは、センサーアレイはアレイ状に並べられるマイクロホン（マイクロホンアレイ）のことである。マイクロホンアレイは、複数の音声を拾い上げる用途、テレコンファレンス用途およびテレプレゼンス用途において、異なるスピーカの声を他の干渉ノイズおよび暗騒音から分離して選択的に増幅するために用いることができる。本明細書に記載の具体例はテレコンファレンスにおけるマイクロホンアレイに関するものであるが、いうまでもなく、本発明は、ビームフォーミングの基本技術に関するものであり、音楽のレコーディングの如き他のオーディオ分野や他の分野、たとえば位置検出または通信のための水面下ハイドロホンアレイなどのソナー、およびセンサー用のアンテナを備えたレーダの如き無線周波数用途にも同様に適用可能である。
【００２６】
好ましい実施形態では、最適化問題、および任意選択的に制約は、センサーアレイの出力を最小限に抑えること、副ローブレベルを最小限に抑えること、主ローブ領域のひずみを最小限に抑えること、およびホワイトノイズ利得を最大化することのうちの１つ以上を満たすように定式化される。これらの要件のうちの１つ以上をビームフォーマの入力パラメータとして選択することができる。さらに、これらの要件のうちのいずれかを最適化問題として定式化することができる。また、要件のうちのいずれかを最適化問題のさらなる制約として定式化することができる。たとえば、最適化問題を、副ローブレベルを最小限に抑えることを条件として、アレイの出力電力を最小限に抑えるように定式化することができ、または、最適化問題を、主ローブ領域のひずみを最小限に抑えることを条件として、副ローブレベルを最小限に抑えるように定式化することができる。所望ならば、個々のビームフォーミング問題に応じて、複数の制約を課すようにしてもよい。
【００２７】
いくつかの好ましい実施形態では、最適化問題は、アレイ出力を最小限に抑えることとして定式化される。これは、システムに課される任意の制約を条件として全域的に最小限に抑えられるパラメータである。したがって、これとは反対に、ビームパターンのある与えられた領域（方向）に制約がない場合、最適化アルゴリズムは、その領域におけるアレイ利得の出力をそのアレイ利得を削減することにより削減しようとする。このことは、利得が必要とされる領域以外のすべての領域で可能なだけ利得を最小限に抑えるという一般的な利点を有している。
【００２８】
好ましくは、入力パラメータは、ビームパターンに主ローブを形成するように、指定された方向のアレイ利得をある与えられたレベルに維持するという要件を含んでいる。上述のように利得を削減する最適化アルゴリズムの一般的な傾向として、指定方向において利得をある与えられたレベルに維持するという要件により、ビームパターンに主ローブ（すなわち、高い利得の領域、すなわち信号を減衰するのではなく信号を増幅する領域）が存在することが担保される。
【００２９】
さらに好ましくは、入力パラメータは、ビームパターンに複数の主ローブを形成するように、複数の指定された方向のアレイ利得をある与えられたレベルに維持するという要件を含んでいる。換言すれば、複数の方向のアレイの利得を選択されたレベルに維持するように、複数の制約を課することよりアレイの指向性が最適化される。このようにして、センサーアレイのビームパターンに複数の主ローブを形成することができ、また、複数の信号ソース方向に対して、その他の方向よりも高い利得が提供することができるようになる。
【００３０】
さらに好ましくは、ビームパターンに異なるレベルの複数の主ローブを形成するように、複数の指定された方向に対して必要な利得レベルが提供される。換言すれば、最適化制約は、異なる方向に異なるレベルの信号メンテナンス（すなわち、アレイ利得）を課すようなものである。たとえば、一方の方向に対して他方の方向よりも高いレベルまたは低いレベルにアレイ利得を維持することができる。このようにして、ビームフォーマは、複数の信号ソースに集中し、それと同時に、それらの信号のレベルを等しくすることができるようになる。たとえば、拾い上げることが必要な３つの信号ソースがあり、３つの信号のうちの２つの信号が３番目の信号よりも強度が高い場合、システムは、ビームパターンに３つの主ローブを形成し、弱い信号に向けられているローブが、強い信号に向けられているローブよりも高い利得を有するようにすることにより、弱いソースを増幅して３つのソースの信号強度を等しくするようにできる。
【００３１】
好ましくは、ビームフォーマは、凸制約として１つ以上の要件を定式化する。さらに好ましくは、ビームフォーマは、線形等式制約として１つ以上の要件を定式化する。このようにして制約を定式化すると、最適化問題は、凸最適化問題のサブセットである二次錘計画法問題となる。二次錘計画法問題の数値解法は詳細に研究されており、また、複数の高速かつ効率的なアルゴリズムが凸二次錘計画法問題の解くのに利用可能となっている。
【００３２】
好ましくは、ビームフォーマは、指定された方向からセンサーアレイに入射する単位強度平面波に対するセンサーアレイ出力が前もって決められている定数と等しいという要件として、１つ以上の主ローブ要件を定式化する。換言すれば、ビームフォーミングパターンは、アレイ出力が指定の方向から入射する平面波に対して特定の利得を提供するように抑制される、  この制約の形態は、線形の等式であるので、上記のように二次錘計画法問題に適用することができる。
【００３３】
本発明の好ましい実施形態では、入力パラメータは、ビームパターンにヌルを形成するように、指定された方向のアレイ利得がある与えられたレベル未満であるという要件を含んでいる。換言すれば、ビームフォーマ最適化問題は、少なくとも一つの方向のアレイ利得が選択されたしきい値未満であるという最適化制約を条件としている。このことにより、ビームパターンの副ローブ領域の最小化を可能とするので、システムの二次マキシマムのサイズが限定される。さらに、このことにより、ビームパターンに「ノッチ」を形成することが可能となり、干渉信号を遮断するために選択された方向に対する利得がとくに低くされる。
【００３４】
さらに好ましくは、入力パラメータは、ビームパターンに複数のヌルを形成するように、複数の指定された方向のアレイ利得がある与えられたレベル未満であるという要件を含んでいる。換言すれば、ビームフォーマ最適化問題は、複数の方向のアレイ利得が対応するしきい値未満であるという最適化制約を条件とするものである。このようにして、複数のヌルをビームパターンに形成することができるので、複数の干渉ソースを抑制することが可能となる。
【００３５】
さらに好ましくは、ビームパターンに異なる深さの複数のヌルを形成するように、最大利得レベルが複数の指定された方向の各々に対して課される。このようにして、異なるレベルの制約をビームパターンの異なる領域に課することができるようになる。たとえば、副ローブをおおむねあるレベル未満に維持することはできるものの、干渉信号を遮断するためにノッチまたはヌルを望む領域には、もっと厳しい制約が課せられる。必要なところだけに最も厳しい制約を課すことにより、ビームパターンの自由度に影響を与えることを少なくし、ビームパターンのその他の部分をより均一に最小化する。
【００３６】
好ましくは、ビームフォーマは、１つ以上の副ローブ要件を凸制約として定式化する。さらに好ましくは、ビームフォーマは、１つ以上副ローブ要件を二次錘制約として定式化する。上記のように、このようにして制約を定式化すると、最適化問題は、凸最適化問題のサブセットである二次錘計画法問題となる。二次錘計画法問題の数値解法は詳細に研究されており、また、複数の高速かつ効率的なアルゴリズムが凸二次錘計画法問題を解くのに利用可能となっている。
【００３７】
最も好ましくは、ビームフォーマは、指定された方向からセンサーアレイに入射する単位強度平面波に対するアレイ出力の強度が前もって決められている定数未満であるという要件として１つ以上の副ローブ要件を定式化する。上記のように、この制約の形態は、凸不等式であるので、上記のような二次錘計画法問題に適用することができる。
【００３８】
好ましくは、入力パラメータは、ビームパターンが、指定されたレベルのロバストネスを有しているという要件を含んでいる。所望のソース信号を拾い上げるということが重大な用途では、単なる軽微なズレ、ランダムノイズまたは他の予期しない干渉によりシステムが故障することのないように担保されていることが望ましい。換言すれば、システムがある程度までエラーに対する回復力があることが望まれる。好ましくは、ロバストネスのレベルは、重み付け係数を有するベクトルのノルム上の制限として指定される。さらに好ましくは、ノルムはユークリッドノルムである。下記にさらに詳細に記載されているように、重み付け係数ベクトルのノルムを最小限に抑えることは、アレイのホワイトノイズ利得を最大化するので、システムのロバストネスを向上させることになる。
【００３９】
好ましくは、重み付け係数は二次錐計画法によって最適化される。上述のように、二次錐計画法は、詳細に研究されてきている凸最適化のサブセットであり、そのような問題を素早く解くための高速かつ効率的なアルゴリズムが利用可能となっている。複数の制約がシステムに課されたとしても、このような数値アルゴリズムは最適化問題のグローバルミニマムに非常に速く収束することができる。
【００４０】
好ましくは、１つ以上の重み付け係数が、位数（ｏｒｄｅｒ）ｎの各球面調和関数に対して最適化されるが、各位数の球面調和関数内では、重み付け係数が、位数ｎのすべての次数（ｄｅｇｒｅｅ）ｍ＝−ｎ〜次数ｍ＝ｎにおいて共通である。このようにして重み付け係数の数を減らすことによって、ビームパターンは、ルック方向に対して回転対称であることに限定される。しかしながら、このようなビームパターンは複数の状況において有益であり、また、係数の数を削減することにより、最適化問題が単純化され、より速く解を算出することが可能となる。
【００４１】
好ましい実施形態によっては、入力信号は、球面調和関数ドメインへ分解される前に、周波数ドメインに変換される場合もある。好ましい実施形態によっては、ビームフォーマは、周波数ドメイン信号が狭帯域周波数ビンに分割される広帯域ビームフォーマであり、各ビンは、これらの周波数ビンが再結合されて広帯域出力を形成する前に、別々に最適化さて、重み付けされる。他の好ましい実施形態では、入力信号は時間ドメインで処理され、また、重み付け係数は、球面調和関数信号に適用される有限インパルス応答フィルタのタップ重みである。
【００４２】
処理ドメインの選択は、個々のシナリオ（すなわち、個々のビームフォーミング問題）の状況に依存する。たとえば、受信されて処理されるであろう予測周波数スペクトルにより、時間ドメインと周波数ドメインとの間の選択が影響を受けることになる。というのは、一方のドメインの方がよりよい解を与え、計算がより効率的であるからである。
【００４３】
時間ドメインでの処理は状況によっては非常に有効である。というのは、時間ドメインは、本来的に広帯域であるからである。したがって、このように実施すると、最適化前に周波数ドメインへフーリエ変換するための集中的な高速計算や、最適化後に時間ドメインへ戻すための逆のフーリエ変換をするための集中的な高速計算を実行する必要がなくなる。また、広帯域の解を得るために入力を複数の狭帯域周波数ビンに分割する必要もなくなる。もっと正確にいえば、すべての重み付け係数を求めるために単一の最適化問題を解くだけでよくなる。いくつかの実施形態では、重み付け係数は、有限インパルス応答（ＦＩＲ）フィルタのタップ重みの形態を取っている。
【００４４】
原則として、ビームフォーミングの性能の観点からすると、時間ドメインにおける実施および周波数ドメインにおける実施は、ＦＩＲ長さがＦＦＴ長さと等しい場合には、同一のビームフォーミング性能を与えることができる。時間ドメインは、実際の実施においては、ＦＦＴおよび逆ＦＦＴを必要としないという、周波数ドメインよりも重要な利点を有している。しかしながら、最適化の複雑さの観点からすると、ＦＩＲとＦＦＴとが同一の長さＬを有していると仮定した場合、単一の最適化によって、１組のＦＩＲ（各チャネル毎にＬ個のＦＩＲ係数）を最適化する計算の複雑さは、Ｌ個のサブバンド最適化によって、１組のアレイ重み（すなわち、各チャネル毎に単一の重み）を最適化することに比べればはるかに大きいと考えられる。したがって、各アプローチは異なる状況において長所を有しうる。
【００４５】
第二の態様によれば、本発明は、各センサーが信号を生成するように構成されているアレイ状に並べられた複数のセンサーと、入力信号を球面調和関数ドメインへ分解し、分解された前記信号を出力するように構成された球面調和関数分解手段と、分解された信号に付与される重み付け係数を１組の入力パラメータに基づいて凸最適化を用いて計算するように構成されている重み付け係数計算手段と、計算された重み付け係数と分解された信号を組み合わせて出力信号を形成する出力形成手段とを備えている。
【００４６】
このようなビームフォーマは、上述のビームフォーミング方法の利点をすべて実現する。さらに、ビームフォーミング方法に関して上述された好ましい特徴はすべて、このビームフォーマの実施形態にも当てはまる。上述のように、時間ドメインでの実施では、出力形成手段は、複数の有限インパルス応答フィルタを有するものとしてもよい。
【００４７】
好ましくは、ビームフォーマは、信号追跡手段をさらに備えており、当該信号追跡手段は、センサーからの信号を評価して所望の信号ソースの方向および不要な干渉ソースの方向を判断するように構成されている。このようなアルゴリズムは、同一のデータを用いて、ビームフォーミング最適化アルゴリズムと並行して動作することができる。位置推定アルゴリズムは、目的信号の方向および干渉ソースの方向を見つけ、ビームフォーマは、ソース信号を増幅し干渉信号を減衰するための適切なビームパターンを形成する。
【００４８】
上述のように、この記載は、主に球面調和関数ドメインでの信号処理に関するものである。しかしながら、本明細書に記載の技術は、他のドメイン、とくにスペースドメインにも適用可能である。一部の用途において、凸最適化がスペースドメイン処理に用いられているが、球面アレイの最適化問題を定式化することもさらなる発明概念であると考えられる。したがって、本発明のさらなる態様によれば、球面センサーアレイのためのビームフォーマのビームパターンを形成する方法であって、ビームフォーマが、センサーアレイから入力信号を受信し、これらの信号に重み付け係数を付与し、これらを組み合わせて出力信号を形成し、これらの重み付け係数が、ある与えられた一組の入力パラメータに対して凸最適化を用いて最適化される。本発明者らは、球面調和関数ドメインに関連して開発した技術および式をスペースドメイン内における球面アレイの処理にも適用可能である点、したがって、本発明により、スペースドメインにおいて実時間で複数の制約含む最適化を実行することも可能である点を認識している。
【００４９】
本発明のさらなる態様によれば、ビームフォーマのビームパターンを形成する方法であって、ビームフォーマが、センサーアレイから入力信号を受信し、これらの入力信号に重み付け係数を付与し、これらを組み合わせて出力信号を形成し、これらの重み付け係数が、ビームパターンに複数の主ローブを形成するように、複数の指定された方向のアレイ利得をある与えられたレベルに維持するという制約を条件として、ある与えられた一組の入力パラメータに対して凸最適化を用いて最適化され、各要件が、指定された方向からセンサーアレイに入射する単位強度平面波に対するセンサーアレイの出力が前もって決められている定数と等しいという要件として定式化される。
【００５０】
上述のように、この記載により導かれるこれらの方法の適用可能性についていえば、システムを実用化できなくなるほど遅くすることなく、最適化問題に複数の制約を課することが可能となる。したがって、本発明の技術および式を用いると、複数のヌル形成およびステアリング制約、ロバストネス制約および主ローブビーム幅制約を課すのと同時に、複数の主ローブの形成および指向性制約を課すことが可能となる。
【００５１】
好ましくは、ビームフォーマは実時間または疑似実時間で動作することができる。いうまでもなく、環境（たとえば、オーディオ用途における音響環境）が固定されている場合、実行時間中にアレイ重みを更新する必要はない。もっと正確にいえば、単一の組の最適化された重みを、前もって（たとえば、システム起動時または校正指令時に）算出しておくことができ、また、その算出された重みを動作中に変更する必要がない。しかしながら、この構成では、本発明のすべての特徴が利用されているわけではない。したがって、センサーアレイは、環境および制約の変化に合わせて最適化問題を解きなおすことにより、最適な重みを動的に変更するようになしてあることが好ましい。上述のように、システムは、実時間または疑似実時間でアレイ重みを再最適化できることが好ましい。実時間の定義は用途によってさまざまである。しかしながら、この記載では、実時間とは、アレイが、１秒間で、アレイ重みを再最適化して新規の最適化されたビームパターンを形成することができることを意味する。疑似実時間とは、最適化にかかる時間が最大約５秒までであることを意味する。このような疑似実時間は、環境の動力学がそれほど急速に変わらない状況、たとえばソースの数および方向がまれにしか変わらない講義中の音響のような場合には有用となりうる。
【００５２】
実時間または疑似実時間における運用では、最適化演算は、重み付けを徐々にかつ連続的に更新するためにバックグラウンド内で実行されることが好ましい。それに代えて、いくつかの状況に対する複数の組の重みを前もって算出してメモリーに格納するようにしておいてもよい。そして、環境の変化にともなって、最も適切な組の重みをシステムに搭載するようにすればよい。しかしながら、いうまでもなく、この実施形態は、実時間における実際の最適化のために本発明の能力および速度を十分に活用しているわけではない。
【００５３】
本発明にかかるビームフォーマは、スペースドメインおよび球面調和関数ドメインにおいて良好に動作することができる。ドメインの選択は、アレイの個々の用途、アレイの幾何学形状、アレイが取り扱うと考えられる信号の特性、および必要となる処理タイプに依存する。一般的にスペースドメインおよび球面調和関数ドメインが最も有用であるが、他のドメイン（たとえば、円筒調和関数ドメイン）も同様に用いられてもよい。それに加えて、処理を周波数ドメインで行ってもよいしまたは時間ドメインで行ってもよい。具体的にいえば、球面調和関数分解をともなう時間ドメイン処理も同様に有用である。したがって、センサー信号が、一組の正規直交基底関数へ分解されて、さらに処理されることが好ましい。最も好ましくは、正規直交基底関数は、球面調和関数、すなわち球面座標系の波動方程式の解であり、波動場分解が球面フーリエ変換によって行なわれる。球面調和関数ドメインは、球面アレイまたは略球面アレイにとくに適している。
【００５４】
さらなる態様によれば、本発明は、センサーアレイのビームフォーマのビームパターンを最適化する方法であって、複数のセンサーからの複数の入力信号が、重み付けされ、組み合わされて、アレイ出力信号を形成し、これらのセンサーの重みが、これらのセンサーの重みの凸関数としてアレイ出力を表現し、１つ以上の制約を条件として出力を最小限に抑えることによって、最適化され、１つ以上の制約が、センサーの重みの凸関数の等式および／または不平等として表現される。
【００５５】
明らかなように、本発明の方法がビームフォーミング問題の一般解を提供する。複数の制約を単一の最適化問題に同時に課し、１つの大域的最適解を得ることができる。しかしながら、課される制約が少ない場合、上術の従来の研究の結果を繰り返すことができる。したがって、本発明は、最適化問題に対するより一般的な解と見なすことができる。
【００５６】
下記に、本システムの好ましい形態についてのさらに詳細な解析を説明する。
【００５７】
通常、空間オーバサンプリング（ｓｐａｔｉａｌ  ｏｖｅｒ−ｓａｍｐｌｉｎｇ）が現実的には用いられるので、次の解析では、より効率的と考えられる球面調和関数ドメイン処理にフォーカスが与えられる。しかしながら、いうまでもなく、球面調和関数ドメイン重み付け関数に関して記載される技術をスペースドメイン解析に同様に適用し、類似の凸最適化問題としてもよい。
【００５８】
背景資料の導出および有用な結果を本出願の付属資料として添付しておく。下記の説明に記載の式番号は付録資料の式番号からの連続番号である。
【００５９】
従来の研究から、周波数に依存しない規則的または不規則なビームパターンを容易に形成するために、アレイ重み付けアプローチでは、球面調和関数ドメインにおけるｂn（ｋａ）の反転が利用され、周波数に依存するコンポーネントが分離される。しかしながら、ｂn（ｋａ）があるｋａ値およびｎ値において小さな値を有し、その反転が実際の実施においてロバストネスを劣化させるので、より一般的な重みであるｗ*（ｋ）ｔが最適化フレームワークの直接目標となる。
【００６０】
次のセクションでは、付属資料で導出された結果を行列式を用いて発展させ、本発明の凸最適化問題および対応する制約が導出される。
【００６１】
次の表記法が用いられる。
【００６２】
【数１】
2012523731.tif 000003
この表記法で、ｖｅｃ（・）は括弧内の項目を積み重ねて（Ｎ＋１）2×１列ベクトルを得ることを表し、（・）Tは転置を表している。
【００６３】
この表記法を用いて、次の式を得ることができる。
【００６４】
【数２】
2012523731.tif 000004
【００６５】
【数３】
2012523731.tif 000005
【００６６】
【数４】
2012523731.tif 000006
【００６７】
【数５】
2012523731.tif 000007
（１８）は、ｂが（ｎ2＋１）から（ｎ＋１）2までのｂnの反復を有していることに留意されたい。（９）から、ｐがモーダルアレイマニホルドベクトル（ｍｏｄａｌ  ａｒｒａｙ  ｍａｎｉｆｏｌｄ  ｖｅｃｔｏｒ）であると考えることができる。
【００６８】
（１４）は次のようにベクトル表記法で書くことができる：
【００６９】
【数６】
2012523731.tif 000008
この式で、（・）Hはエルミート転換を表す。
【００７０】
次の記載では、最適化問題が、ビーム方向以外から入ってくる干渉を抑えるためにアレイ出力を最小限に抑えるとともに、主ローブ方向からの信号を維持し、副ローブを制御するものとして定式化されている。さらに、ビームフォーマのロバストネスを向上させるために、アレイ重みのノルムを指定された定数に制限するために、ホワイトノイズ利得制約が課される。
【００７１】
アレイ出力は次の式で与えられる：
【００７２】
【数７】
2012523731.tif 000009
この式で、Ｅ［・］はブラケット内の量の統計的期待値を表し、Ｒ（ω）はｘの共分散行列（スペクトル行列）である。
【００７３】
Ｈ（ｋａ、Ω）により表される指向性パターンは、すべての目標角度からの単位入力信号に対するアレイ応答の関数である。したがって、次の式で表される：
【００７４】
【数８】
2012523731.tif 000010
信号ソースが相関関係を有していないと仮定すると、ｘの共分散行列は次の式で表される：
【００７５】
【数９】
2012523731.tif 000011
この式で、
【数１０】
2012523731.tif 000012
は、Ｄ＋１個の相関関係のない信号のパワー（ｐｏｗｅｒ）であり、Ｑ（ω）＝Ｅ［Ｎ（ω）ＮH（ω）］は、
【数１１】
2012523731.tif 000013
でのノイズ共分散行列である。
【００７６】
ここで、ノイズ場の特別なケースである等方性ノイズ、すなわち球面上に均一に分配されるノイズを想定する。パワースペクトル密度
【数１２】
2012523731.tif 000014
を有している等方性ノイズというのは、全方向から球面に入ってくる相関関係のない均一なパワー密度
【数１３】
2012523731.tif 000015
を有している無数の平面波が存在しているものと考えることができる。したがって、すべての方向に対して共分散行列を積分することにより、等方性ノイズ共分散行列は次の式で与えられる：
【００７７】
【数１４】
2012523731.tif 000016
（７）、（１８）および（１９）を用いて、（２５）を次のように書き直すことができる。
【００７８】
【数１５】
2012523731.tif 000017
この式で、
【数１６】
2012523731.tif 000018
は２つのベクトルのアダマール（すなわち、各要素毎の）積を表している。 上記の導出において、球面調和関数正規直交特性（４）が用いられていることに留意されたい。
【００７９】
実際の用途では、正確な共分散行列Ｒ（ω）は入手不可能である。したがって、式（２４）の代わりに、標本共分散行列が用いられる。標本共分散行列は次の式で与えられる：
【００８０】
【数１７】
2012523731.tif 000019
この式で、Ｉはスナップショットの数である。
【００８１】
アレイ利得Ｇ（ｋ）は、アレイ出力での信号対ノイズ比（ＳＮＲ）の、入力センサーでのＳＮＲに対する比率であると定義されている：
【００８２】
【数１８】
2012523731.tif 000020
この式で、
【数１９】
2012523731.tif 000021
は正規化されたノイズ共分散行列である。
【００８３】
アレイの性能の一般的な尺度は指向性である。指向係数Ｄ（ｋ）、すなわち指向利得は、等方性ノイズに対するアレイ利得として解釈することができる。（２７）の中のＱをＱisoで交換することにより次の指向係数が与えられる。
【００８４】
【数２０】
2012523731.tif 000022
次いで、指向性指数（ＤＩ）がＤＩ（ｋ）＝１０ｌｏｇ10Ｄ（ｋ）ｄＢとして定義される。
【００８５】
ビームフォーマの能力を評価しうる性能尺度が複数ある。一般的に用いられているアレイの性能尺度は、指向性、アレイ利得、ビーム幅、副ローブレベルおよびロバストネスである。
【００８６】
これらの相反する性能尺度間のトレードオフはビームフォーマ設計最適化問題を意味する。本発明にかかる方法では、最適化問題は、目的信号（ＳＯＩ）の無歪制約（すなわち、ビームパターンに主ローブを形成すること）と共に、任意の数の他の所望の制約、たとえば副ローブ制約およびロバストネス制約を条件として、出力を最小限に抑えることに向けられる。アレイ重みベクトルｗ（ｋ）を最適化変数とすると、多重制約ビームフォーミング最適化問題を次のように定式化することが可能となる：
【００８７】
【数２１】
2012523731.tif 000023
この式で、ΩSLは副ローブ領域であり、εおよびζは、それぞれ、副ローブおよびホワイトノイズ利得（すなわち、ホワイトノイズに対するアレイ利得）ＷＮＧを制御するためのユーザパラメータである。ホワイトノイズ利得制約はビームフォーマのロバストネスを向上させるために一般的に用いられている。ルック方向（すなわち、主ローブの方向）は、Ω0であり、ＳＯＩの到来方向（入ってくる方向または入射方向）である。
【００８８】
ワイトノイズ利得（ＷＮＧ）は次の式で与えられる：
【００８９】
【数２２】
2012523731.tif 000024
（１５）を用いると、ＷＮＧを次のように書き直すことができる：
【００９０】
【数２３】
2012523731.tif 000025
ホワイトノイズ利得が重みベクトルのノルムに反比例することが分かる。ビームフォーマのロバストネスを向上させるために、分母、すなわちアレイ重みのノルムをあるしきい値に制限するようにしてもよい。
【００９１】
隣接している方向の応答と応答との間の相関関係に起因して、Ωl∈ΘSL、ｌ＝１、．．．、Ｌの方向の有限数の格子点を用いて副ローブ領域ΩSLを近似することができる。Ｌの選択は必要とされる近似の精度によって決まる。
【００９２】
（２３）および（３１）を用いると、（２９）は次のように表される：
【００９３】
【数２４】
2012523731.tif 000026
この式で、‖・‖はユークリッドノルムを表している。
【００９４】
二次錐計画法は、一組の二次錘制約および一組の線形等式制約を条件として、線形関数を最小化する一般的な凸計画法問題の下位分類である。この問題は次のように記載することができる：
【００９５】
【数２５】
2012523731.tif 000027
この式で、
【数２６】
2012523731.tif 000028
であり、
【数２７】
2012523731.tif 000029
であり、
【数２８】
2012523731.tif 000030
であり、
【数２９】
2012523731.tif 000031
であり、
【数３０】
2012523731.tif 000032
であり、
【数３１】
2012523731.tif 000033
であり、
【数３２】
2012523731.tif 000034
であり、
【数３３】
2012523731.tif 000035
であり、
【数３４】
2012523731.tif 000036
であり、
【数３５】
2012523731.tif 000037
およびＣは一組の実数および複素数（または、行列）である。
【００９６】
上述の（３２）に記載の最適化問題において、便宜上、一時的に独立変数ω、ｋを除外し、
【００９７】
【数３６】
2012523731.tif 000038
をＲのコレスキー分解とすると次の式が得られる：
【００９８】
【数３７】
2012523731.tif 000039
新しい非負スカラー変数ｙlを導入し、ｙ＝［ｙl、ｗT］Tであり、ｂ＝［ｌ、０T］Tであると定義し、０が所望の次元のゼロベクトルであるとすると、最適化問題（３２）を次のように書き直すことができる：
【００９９】
【数３８】
2012523731.tif 000040
この式で、Ｉは単位行列である。このようにして、最適化問題（３２）は、二次錘計画法問題の形式に書き直される。したがって、数値法を用いてこの問題の解を効率的に見つけることができる。最適化問題の解を求めた後、変数ｙのベクトルの目的パラメータはそのサブベクトルｗにより与えられる。
【０１００】
したがって、この最適化問題は、１組の二次錘制約および一組の線形等式制約を条件として、線形関数を最小限に抑える凸二次錐計画法（ＳＯＣＰ）問題として定式化されることが分かる。これは、より一般的な凸計画法問題の下位分類である。ＳＯＣＰ問題は、計算可能であり、既知の数値ソルバーを用いて効率的に解を求めることができる問題である。このような数値ソルバーの一例がＭＡＴＬＡＢで使用可能なＳｅＤｕＭｉソルバー（ｈｔｔｐ：／／ｓｅｄｕｍｉ．ｉｅ．ｌｅｈｉｇｈ．ｅｄｕ／）である。存在するならば、ＳＯＣＰ問題のグローバル最適数値解が保証される。すなわち、ＳＯＣＰ問題にグローバル最小値が存在するならば、数値ソルバーのアルゴリズムはそれを確実に見つける。さらに、これらの技術が高速計算可能であるので、実時間での最適化を維持しつつ、最適化問題に複数の制約を含めることができる。ＳＯＣＰは、一般的な凸最適化より計算がより効率的であので、実時間用途に非常に好適である。
【０１０１】
計算の複雑さに関していえば、上述の（３２．３）で導出されたＳＯＣＰ問題を解決するために内点法が用いられる場合、双対性ギャップ（ｄｕａｌｉｔｙ  ｇａｐ）をそれ自体の一定の割合に削減させるための反復回数は、
【数３９】
2012523731.tif 000041
（この式で、項「１」は等式制約（ｅｑｕａｌｉｔｙ  ｃｏｎｓｒａｉｎｔ）によるものである）により上方で境界され、また、反復毎の計算量は
【数４０】
2012523731.tif 000042
である。
【０１０２】
最適化問題（３２．２）の場合、反復毎の計算量は
【数４１】
2012523731.tif 000043
であり、反復回数は
【数４２】
2012523731.tif 000044
である。通常、アルゴリズムは反復回数１０未満で収束する（最適化分野において広く認められている事実）。
【０１０３】
本発明の好ましい実施形態について説明する前に特筆すべき点は、上記の解析はすべて、信号ソースが、遠距離場にあり、そのためアレイに入射する平面波として近似可能であるという仮定に基づいているということである。
【０１０４】
さらに特筆すべき点は、解析が狭帯域ビームフォーマ設計に基づいているということである。広帯域ビームフォーマは、周波数帯域をより狭い周波数ビンに分解し、狭帯域ビームフォーマにより各ビンを処理することによって簡単に実現することができる。
【０１０５】
時間ドメインで実施される場合、広帯域ビームフォーマを達成するためには、適切な時間遅延および重みを各副帯域（ｓｕｂ−ｂａｎｄ）のセンサーの各々に課してビームパターンを形成するが、またはそれに代えて、ＦＩＲおよび重み付け方法（ＦＩＲ−ａｎｄ−ｗｅｉｇｈｔ  ｍｅｔｈｏｄ）を用いて時間ドメインの広帯域ビームフォーミングを達成することができる。しかしながら、周波数ドメインで実施される場合、それぞれの狭い周波数ビンについて、複素数重み（ｃｏｍｐｌｅｘ  ｗｅｉｇｈｔ）が各センサーに付与される。上記の記載は、周波数ドメインでの実施にフォーカスを与え、各周波数ついて複素数重みを最適化している。時間ドメインでの実施についての詳細な説明は下記になされている。
【０１０６】
上記のアプローチは、周波数ドメインにおける信号モデルに基づいており、このモデルでは、複素数値のモーダル変換およびアレイ処理が用いられている。スピーチ用途およびオーディオ用途において非常に重要な広帯域ビームフォーマを実現するために、広帯域アレイ信号が、離散的フーリエ変換（ＤＦＴ）を用いてより狭い周波数ビンへ分解され、次いで各周波数ビンが狭帯域ビームフォーミングアルゴリズムを用いて独立して処理され、次いで、広帯域出力信号が離散的逆フーリエ変換を用いて生成される。周波数ドメインでの実施がブロック化処理でなされるので、それに付随する時間遅れのため、時間が重要となる（ｔｉｍｅ−ｃｒｉｔｉｃａｌ）スピーチ用途およびオーディオ用途には適さない場合もある。
【０１０７】
古典的エレメント−スペースアレー処理（ｅｌｅｍｅｎｔ  ｓｐａｃｅ  ａｒｒａｙ  ｐｒｏｃｅｓｓｉｎｇ）では、フィルタ・アンド・サム構造（ｆｉｌｔｅｒ−ａｎｄ−ｓｕｍ  ｓｔｒｕｃｔｕｒｅ）を用いて広帯域ビームフォーマを時間ドメインで具象化することができることは周知のことである、フィルタ・アンド・サム構造では、ひとそろいの有限インパルス応答（ＦＩＲ）フィルタをセンサーの出力部に設け、フィルタ出力を合計して最終出力時系列を形成するようになしてある。時間ドメインでフィルタ・アンド・サムを具象化する主な利点は、新規のスナップショットが到着する毎に、実行時にビームフォーマを更新することができるということである。フィルタ・アンド・サムビームフォーマ設計のキーポイントは、所望のビームフォーミング性能を達成するためにＦＩＲフィルタのタップ重みをどのように計算するかということである。
【０１０８】
また、球面アレイモーダルビームフォーミングを実数値モーダル変換およびフィルタ・アンド・サムビームフォーミング構造を用いて時間ドメインで具象化することができる。ＷＯ０３／０６１３３６は、球面アレイモーダルビームフォーマを時間ドメインで具象化するための新規な球面調和関数フレームワーク内の構造を提案している。この具体例では、信号処理チャネルの数が著しく削減され、球面調和関数の実数部分および虚数部分が、時間ドメイン広帯域信号を実数値球面調和関数ドメインへ変換するための球面フーリエ変換の基礎として用いられ、ビームフォーマのルック方向がビームパターン形状から巧妙に分離可能となっている。周波数に依存しないビームパターンを達成するために、ＷＯ０３／０６１３３６では、各信号チャンネルの周波数依存コンポーネントを分離するために逆フィルタを用いることが提案されているが、このような種類の逆フィルタリングは、システムのロバストネスを損傷する恐れがある（ジェー・マイヤーおよびジー・エルコ、「音場の正規直交分解に基づいた高度にスケーリング可能な球面マイクロホンアレイ」、ＩＣＡＳＳＰ会報、第２巻、ページ１７８１〜１７８４、２００２年５月）。さらに、このようなフィルタ・アンド・サムモーダルビームフォーミング構造に関する系統的な性能解析フレームワークが定式化されていないので、指向性係数、副ローブレベル、ロバストネスなどの如きすべての相反する広帯域ビームフォーミング性能尺度を効果的に制御することができない。
【０１０９】
ここで、時間ドメインで具象化される広帯域モーダルビームフォーミングフレームワークを説明する。この技術は、フィルタ・アンドサムモーダルビームフォーミングス構造を修正したものに基づいている。アレイ応答、等方性ノイズおよび空間ホワイトノイズの両方に対するビームフォーマ出力および主ローブ空間応答変化（ＭＳＲＶ）のための表現を、ＦＩＲフィルタタップ重みについて導き出す。複数の相反する性能尺度（たとえば、指向性指数、ロバストネス、副ローブレベル、主ローブ応答変化など）間の適切なトレードオフを実現するために、ＦＩＲフィルタタップ重み設計問題が、複数の制約を有する計算可能な最適化問題として定式化される。
【０１１０】
それに加えて、ここに記載の構成で、ステアリングユニットについて説明する。ステアリングユニットにより信号処理チャネルの数が削減され、また、古典的エレメントスペースアレイ処理と比較して、モーダルビームフォーミングアプローチは計算という点においてより効率的である。ステアリングユニットは、ルック方向に対して回転対称となるビームパターンを形成することにより、計算の複雑さを減少させている。上述の非対称ビームパターンほど一般的ではないものの、このような構造はいまだに有用であることが多い。しかしながら、いうまでもなく、ステアリングユニットが下記の時間ドメインビームフォーマの主要なコンポーネントではなく、また、より一般的なビームパターン形成が望まれる場合には、ステアリングユニットが省略されてもよい。
【０１１１】
下記では、周波数ドメインアプローチのために先に導出した結果のうちのいくつかが改めて定式化され、ビームステアリングユニットが加えられる。ｓ番目のマイクロホンで受信される時系列がｘs（ｔ）であり、周波数ドメイン表記がｘ（ｆ、Ω）であると仮定する。ｘ（ｆ、Ω）の離散的球面フーリエ変換（球面フーリエ係数）が次のように与えられる：
【０１１２】
【数４３】
2012523731.tif 000045
（Ｔ５）を用いて、音場が、時間ドメインまたは周波数ドメインから球面調和関数ドメインに変換される。
【０１１３】
各マイクロホンがｗ*（ｆ、Ωs）で表される重みを有していると仮定する。 ｙ（ｆ）で表されるアレイ出力を次のように計算することができる：
【０１１４】
【数４４】
2012523731.tif 000046
この式で、
【数４５】
2012523731.tif 000047
はｗ*（ｆ、Ωs）の球面フーリエ係数である。（Ｔ６）内の第二の合計項は球面調和関数ドメイン内の重みであると考えることができる。
【０１１５】
先の場合と同様に、次の表記が用いられる：
【０１１６】
【数４６】
2012523731.tif 000048
この表記で、ｖｅｃ（・）は括弧内の項目を積み重ねて（Ｎ＋１）2×１列ベクトルを得ることを表し、（・）Tは転置を表している。
【０１１７】
（Ｔ６）をベクトル表記法で書き直すことができる。
【０１１８】
【数４７】
2012523731.tif 000049
この式で、
【数４８】
2012523731.tif 000050
である。
【０１１９】
アレイ出力は次の式で与えられる：
【０１２０】
【数４９】
2012523731.tif 000051
この式で、Ｅ［・］はブラケット内の量の統計的期待値を表し、Ｒb（ｆ）はｘbの共分散行列（スペクトル行列）である。
【０１２１】
Ｂ（ｆ、Ω）により表される指向性パターンは、すべての目標角度からの単位入力信号に対するアレイ応答の関数である。したがって、
【０１２２】
【数５０】
2012523731.tif 000052
重み付けに対する球面フーリエ変換に関するパーセバルの関係を適用することによって、次の式が得られる：
【０１２３】
【数５１】
2012523731.tif 000053
直観的に、マイクロホンが球面上に一様に分配されることが望ましい。しかしながら、真の等距離での空間サンプリングは、５つの正多面体幾何学形状、すなわち四面体、立方体、八面体、十二面体および２０面体に従って構築される構造に対してのみ可能である。ほぼ均一なサンプリングスキームを提供する配置が用いられている。この配置では、３２のマイクロホンが、切頭２０面体の面の中心に配置される。球面アレイにとって良好に働くと考えられている簡単でかつほぼ均一な格子の他の具体例は、Ｆｌｉｅｇｅ格子である。これらのほぼ均一なケースでは、
【数５２】
2012523731.tif 000054
である。
【０１２４】
ルック方向Ω0に対して回転対称であるビームパターンを形成するために、アレイ重みは次の式で表される：
【０１２５】
【数５３】
2012523731.tif 000055
この式で、
【数５４】
2012523731.tif 000056
は、ルック方向をΩ0だけ移動させる役割を有するステアリングユニットとして働き、Ｃn（ｆ）はパターン生成の役割を有している。
【０１２６】
（Ｔ６）において（Ｔ１２）を用いると、次の式が得られる：
【０１２７】
【数５５】
2012523731.tif 000057
（Ｔ５）および（Ｔ１３）に従って、図２０に記載のようなモーダルビームフォーマ構造が得られる。まず、音場データｘ（ｆ、Ω）が時間ドメインまたは周波数ドメインから球面調和関数ドメインデータｘnm（ｆ）に変換される。次いで、ハーモニックスドメインデータｘnm（ｆ）がモーダルビームフォーマへ直接送られる（ステアリング、重み付け、および合計）。これは、マイヤーおよびエルコによって次の文献に記載されたものとの違いである：「音場の正規直交分解に基づいた高度にスケーリング可能な球面マイクロホンアレイ」、ＩＣＡＳＳＰ会報、第２巻、ページ１７８１〜１７８４、２００２年５月。この文献では、代わりに、ｂnについて補償されている球面調和関数がモーダルビームフォーマに送られる。この変更は、補償ユニットによって引き起こされるビームフォーマの劣悪なロバストネスを回避するためになされている。
【０１２８】
（Ｔ１０）において（Ｔ１２）、（５）および（７）を用いると、次の式が与えられる：
【０１２９】
【数５６】
2012523731.tif 000058
この式で、Ｐnはルジャンドル多項式であり、ΘはΩとΩ0との間の角度である。
【０１３０】
ロバストネスは、アレイの性能の重要な尺度であり、ホワイトノイズ利得（ＷＮＧ）によって定量化される。すなわちホワイトノイズに対するアレイ利得のことである。（Ｔ１１）を用い、
【数５７】
2012523731.tif 000059
であると仮定すると、ＷＮＧは次の式で与えられる：
【０１３１】
【数５８】
2012523731.tif 000060
この式で、ｃ＝［Ｃ0、．．．、Ｃn、．．．、ＣN］Tは（Ｎ＋１）×１の列ベクトルのことである。
【０１３２】
最大ＤＩモーダルビームフォーマおよび最大ＷＮＧモーダルビームフォーマの場合、次の式が得られる：
【０１３３】
【数５９】
2012523731.tif 000061
【０１３４】
【数６０】
2012523731.tif 000062
この式で、下付き文字ＭＤＩおよびＭＷＮＧはそれぞれ最大ＤＩビームフォーマおよび最大ＷＮＧビームフォーマを表している。
【０１３５】
ここまでのところ、モーダル変換モおよびビームフォーミングの数理解析が複素数球面調和関数（ｃｏｍｐｌｅｘ  ｓｐｈｅｒｉｃａｌ  ｈａｒｍｏｎｉｃｓ）について説明されている。次に、広帯域モーダルビームフォーマの時間ドメインでの具象化について検討する。時間ドメインで実施する場合、実数値係数がより適切であるので、球面調和関数ドメインデータの実数部分および虚数部分を用いることができる。
【０１３６】
ｓ番目のマイクロホンで受信されるサンプリングされた広帯域時系列がｘs（ｌ）＝ｘs（ｔ）｜t=lTであり、この式で、Ｔsがサンプリング間隔であると仮定する。（Ｔ５）と同様に、
【数６１】
2012523731.tif 000063
が周波数に依存しないと想定すると、広帯域球面調和関数ドメインデータは次の式で与えられる：
【０１３７】
【数６２】
2012523731.tif 000064
この式で、ｘnm（ｌ）は（Ｔ５）内のｘnm（ｆ）の時間ドメイン表記、すなわちｘnm（ｆ）のフーリエ逆変換であり、
【数６３】
2012523731.tif 000065
は入力データの長さである。
【０１３８】
フィルタ・アンド・サム構造は、古典的エレメントスペースアレイ処理における広帯域ビームフォーミングで用いられている。古典的エレメントスペースアレイ処理では、各センサーがＦＩＲフィルタに供給し、フィルタ出力が合計されビームフォーマ出力時系列を生じる。古典的アレイ処理の類似性を用いて、フィルタ・アンド・サム構造をモーダルビームフォーマに適用することができる。すなわち、ステアリングユニットの出力部にひとそろいの実数値ＦＩＲフィルタが設けられ、これらのフィルタが、広帯域周波数バンドの複素数重み（ｃｏｍｐｌｅｘ  ｗｒｉｇｈｔｉｎｇ）ｃn（ｆ）の役割をなす。ステアリングユニットを備えたモーダルビームフォーマの利点は、Ｍ個のフィルタを必要とする古典的エレメントスペースビームフォーマとは対照的に、Ｎ＋１個のＦＩＲフィルタしか必要としないので、計算が効率的であるということである。Ｍ≧（Ｎ＋１）2であることに留意されたい。特筆すべき点は、ステアリングユニットは本発明の任意選択な構成要素であり、用いられない場合、（Ｎ＋１）2個の球面調和関数（
【数６４】
2012523731.tif 000066
）の各々についてＦＩＲフィルタが用いられるということである。
【０１３９】
ｈnを、位数ｎの球面調和関数に対応するＦＩＲフィルタのインパルス応答であるとする、すなわちｈn＝［ｈn1、ｈn2、．．．、ｈnL］T、ｎ＝０、．．．、Ｎであるとする。ここで、ＬはＦＩＲフィルタの長さである。（Ｔ１３）に対してフーリエ逆変換を行ない、動作周波数バンド上のフィルタの応答をｃn（ｆ）にほぼ等しいと想定すると、
【数６５】
2012523731.tif 000067
により表される時間ドメインビームフォーマ出力を次の２つの式で与えることができる：
【０１４０】
【数６６】
2012523731.tif 000068
この式で、＊は畳み込み（コンボリューション）を表し、
【０１４１】
【数６７】
2012523731.tif 000069
であり、
ここで、Ｒｅ（・）およびＩｍ（・）はそれぞれ実数部および虚数部を表しており、
【数６８】
2012523731.tif 000070
であり、
【数６９】
2012523731.tif 000071
である。上述の導出には特性
【数７０】
2012523731.tif 000072
が利用されていることに留意されたい。
【０１４２】
（Ｔ２０）において（３）を用いると、次の式が与えられる：
【０１４３】
【数７１】
2012523731.tif 000073
（Ｔ１９）および（Ｔ２１）に従って、広帯域モーダルビームフォーマの時間ドメインでの具象化を図２１のように示すことができる。プリディレイＴ0が各高調波（ｈａｒｍｏｎｉｃｓ）に対してＦＩＲフィルタの前に取り付けられることに留意されたい。このプリディレイは、通常Ｔ0＝−（Ｌ−１）Ｔs／２として選択されるＦＩＲフィルタの固有の群遅延を償うために用いられる。この目的は、これらのＦＩＲフィルタのインパルス応答（または、タップ重み）を選択してモーダルビームフォーマの所望の周波数−波数応答を達成することである。
【０１４４】
インパルス応答ｈnを備えたＦＩＲフィルタの複素周波数応答は次の式で与えられる：
【０１４５】
【数７２】
2012523731.tif 000074
この式で、ｅ（ｆ）＝［１、ｅ-j2πfTs、．．．、ｅ-j(L-1)2πfTs］Tである。
【０１４６】
η＝ｅ-j2πfT0とする。周波数ｆでの次数ｎの球面調和関数に対応するパターン生成ユニットの総合重み関数は、次の式で与えられる：
【０１４７】
【数７３】
2012523731.tif 000075
（Ｔ１４）においてｃn（ｋ）に代えて（Ｔ２３）の
【数７４】
2012523731.tif 000076
を用いると、次の式が得られる：
【０１４８】
【数７５】
2012523731.tif 000077
【０１４９】
【数７６】
2012523731.tif 000078
であり、ａ＝「ａ0、．．．、ａn、．．．、ａN］Tであるとし、（Ｎ＋１）Ｌ×１複合ベクトル
【数７７】
2012523731.tif 000079
を定義する。式（Ｔ２４）を次のように書き直すことができる。
【０１５０】
【数７８】
2012523731.tif 000080
この式で、
【数７９】
2012523731.tif 000081
はクロネッカー積を表し、
【数８０】
2012523731.tif 000082
である。
【０１５１】
αs＝４π／Ｍである場合、（Ｔ６）のアレイ出力振幅は、古典的アレイ処理
【数８１】
2012523731.tif 000083
よりも４π／Ｍ倍だけ大きいことに留意されたい。したがって、球面調和関数ドメインの無歪制約は次のようになる：
【０１５２】
【数８２】
2012523731.tif 000084
ここで、球面等方性ノイズ、すなわちノイズが球面上に均一に分布する、ノイズフィールドの特別なケースを想定する。スペクトル密度
【数８３】
2012523731.tif 000085
を有する等方性ノイズは、全方向Ωから均一な出力密度
【数８４】
2012523731.tif 000086
で球面に到来する相関関係のない平面波が無数に存在しているものである想定することができる。したがって、すべての方向の共分散行列を積分することによって、等方性ノイズ共分散行列は次の式で与えられる：
【０１５３】
【数８５】
2012523731.tif 000087
この式で
【数８６】
2012523731.tif 000088
であり、
【数８７】
2012523731.tif 000089
であり、
【数８８】
2012523731.tif 000090
であり、
【数８９】
2012523731.tif 000091
は２つのベクトルのアダマール（すなわち、エレメント毎の）積であることを表し、ｄｉａｇ｛・｝は対角線上に独立変数（ａｒｇｕｍｅｎｔ）エレメントを有する正方行列を表している。上述の導出において球面調和関数正規直交特性が用いられていることに留意されたい。
【０１５４】
マイクロホンアレイに等方性ノイズだけが入射する特別なケースを想定する。Ｒb（ｆ）が等方性ノイズ共分散行列Ｑbiso（ｆ）と交換された（Ｔ９）が用いられ、Ｐisoout（ω）で表される等方性ノイズのみのビームフォーマ出力が得られる。
【０１５５】
【数９０】
2012523731.tif 000092
この式で、
【０１５６】
【数９１】
2012523731.tif 000093
である。
ここで、ｂc（ｋａ）＝［ｂ0（ｋａ）、ｂ1（ｋａ）、ｂ2（ｋａ）、．．．、ｂN（ｋａ）］Tである。
【０１５７】
（Ｔ２３）を用い、
【数９２】
2012523731.tif 000094
であると表すと、次の式が与えられる：
【０１５８】
【数９３】
2012523731.tif 000095
（Ｔ２９）においてｃ（ｋ）に代えて
【数９４】
2012523731.tif 000096
を用いると、次の式が与えられる：
【０１５９】
【数９５】
2012523731.tif 000097
この式で、
【数９６】
2012523731.tif 000098
はｈに関する等方性ノイズ共分散行列である。
【０１６０】
ｆLおよびｆUがそれぞれ下限周波数および上限周波数である周波数バンド［ｆL、ｆU］を占める広帯域等方性ノイズの場合、その広帯域共分散行列
【数９７】
2012523731.tif 000099
を、領域［ｆL、ｆU］においてｆに対して積分することにより得ることができる：
【０１６１】
【数９８】
2012523731.tif 000100
この式で、積分を、加算を行なうことによって近似させることができる。
【０１６２】
空間ホワイトノイズが周波数バンド［ｆL、ｆU］にわたって平坦なスペクトル
【数９９】
2012523731.tif 000101
を有していると仮定する。広帯域等方性ノイズのみのビームフォーマ出力は次のようになる：
【０１６３】
【数１００】
2012523731.tif 000102
パワースペクトル密度
【数１０１】
2012523731.tif 000103
を有する空間ホワイトノイズのみがマイクロホンアレイに衝突する特別なケースを想定する。
【数１０２】
2012523731.tif 000104
の場合、Ｐwoutで表される空間ホワイトノイズのみのビームフォーマ出力は次の式で与えられる：
【０１６４】
【数１０３】
2012523731.tif 000105
空間ホワイトノイズが、周波数バンド［０、ｆs／２」全域にわたって平坦なスペクトル
【数１０４】
2012523731.tif 000106
を有していると仮定する。
【数１０５】
2012523731.tif 000107
で表される広帯域ビームフォーマ出力は次の式で与えられる：
【０１６５】
【数１０６】
2012523731.tif 000108
次いで、ＢＷＮＧにより表される広帯域ホワイトノイズ利得は次のように定義される：
【０１６６】
【数１０７】
2012523731.tif 000109
アレイの性能の一般的な尺度は指向性である。指向係数Ｄ（ｆ）、すなわち指向利得は、等方性ノイズに対するアレイ利得として解釈することができ、次の式で与えられる：
【０１６７】
【数１０８】
2012523731.tif 000110
ほとんどの場合、指向係数は、ｄＢで表され、指向性指数（ＤＩ）と呼ばれる。ＤＩ（ｆ）＝１０ｌｇＦ（ｆ）、ここで、ｌｇ（・）＝ｌｏｇ10（・）である。
【０１６８】
主ローブ空間応答変化（ＭＳＲＶ）は、次のように定義される：
【０１６９】
【数１０９】
2012523731.tif 000111
この式で、ｆ0は選択された基準周波数である。
【０１７０】
ｆk∈［ｆL、ｆU］（ｋ＝１，２、．．．、Ｋ）、Θj∈ΘML（ｊ＝１，２、．．．、ＮML）およびΘi∈ΘSL（ｉ＝１，２、．．．、ＮSL）をそれぞれ周波数バンド［ｆL、ｆU］、主ローブ領域ΘMLおよび副ローブ領域ΘSLを近似する選択された（均一または非均一）格子であるとする。ＮMLＫ×１列ベクトルＢγMSRVおよびＮSLＫ×１列ベクトルＢSLが定義され、それぞれの項目が次のように与えられる：
【０１７１】
【数１１０】
2012523731.tif 000112
【０１７２】
【数１１１】
2012523731.tif 000113
次いで、γMSRVのノルム、すなわち‖γMSRV‖qを合成広帯域ビームパターンの周波数に依存しない近似の尺度として用いることができる。下付き文字ｑ∈｛２、∞｝は、それぞれｌ2（ユークリッド）ノルムおよびｌ∞（Ｃｈｅｂｙｓｈｅｖ）ノルム表している。同様に、‖ＢSL‖qは副ローブ挙動の尺度である。
【０１７３】
ビームフォーマの能力を評価するための性能尺度が複数存在する。一般的に用いられているアレイの性能尺度は指向性、ＭＳＲＶ、副ローブレベルおよびロバストネスである。これらの相反する性能尺度間のトレードオフはビームフォーマ設計最適化問題を意味する。広帯域球面調和関数ドメインビームパターンＢ（ｆ、Ω）（Ｔ２５）、広帯域等方性ノイズのみのビームフォーマ出力
【数１１２】
2012523731.tif 000114
（Ｔ３４）、広帯域ホワイトノイズ利得ＢＷＮＧ（Ｔ３７）、主ローブ空間応答変化ベクトルγMSRV（Ｔ４０）および副ローブ挙動ベクトルＢSL（Ｔ４１）を定式化した後、広帯域モーダルビームフォーマのための最適アレイパターン合成問題を次のように定式化することができる：
【０１７４】
【数１１３】
2012523731.tif 000115
この式で、ｑ1、ｑ2∈｛２、∞｝および
【数１１４】
2012523731.tif 000116
は、コスト関数と、３つのユーザパラメータとを有している。上述の周波数ドメイン問題と同様に、最適化問題（Ｔ４２）を、凸形態で表すこととができ、いわゆる二次錘計画法（ＳＯＣＰ）として定式化することができる。二次錘計画法（ＳＯＣＰ）は、ＳｅＤｕＭｉの如きＳＯＣＰソルバーを用いて効率的に解くことができる。
【０１７５】
（Ｔ４２）は、ビームフォーミングの目的に応じて適切な最適化問題を定式化するために用いることができる一般的な式として与えられている。たとえば、４つの関数（ｌ＝１、２、３、４）のうちのいずれかを目的関数として用い、その他の関数のうちのいずれかを制約（関数）として用いることができる。ｌ＝１の場合、最適化問題は、アレイの出力を最小限に抑えることとして定式化される。ｌ＝２の場合、最適化問題は、主ローブ領域のひずみを最小限に抑えることである。ｌ＝３の場合、最適化問題は、副ローブレベルを最小限に抑えることであり、ｌ＝４の場合、最適化問題は、ホワイトノイズ利得（ロバストネス）を最大化させることである。それぞれの場合において、最適化問題を他の制約のうちのいずれかまたはすべてを条件として定式化することができる。たとえば、ｌ＝２を目的関数とし、ｌ＝１、ｌ＝３およびｌ＝４を最適化問題のさらなる制約として、最適化問題を定式化することができる。このようにして、このビームフォーマを非常に柔軟なものにすることができることが分かる。
【０１７６】
この構成では、フィルタタップ重みが、与えられた一組の入力パラメータについて、凸最適化により最適化される。センサーアレイからの入力信号が球面調和関数ドメインへ分解され、次いで、分解された球面調和関数コンポーネントが、ＦＩＲタップ重みにより重み付けされた後、組み合わせられて出力信号を形成する。
【０１７７】
特筆すべき点は、この記載では、ほとんど電話コンファレンスに関する具体例が提供されているが、本発明は、電話コンファレンス用途に制限されるわけではないということである。もっと正確にいえば、本発明は、他の技術分野にも同様に適用可能なビームフォーミング方法に関するものである。これらには、非常に複雑な聴覚情景のうちのある領域を強調または非強調することが望まれる、高性能周辺音響システムおよび音楽録音システムのためのアンビソニックス（ａｍｂｉｓｏｎｉｃｓ）が含まれる。このような用途には、とくに本発明にかかる多重主ローブ指向性およびレベル制御ならびにそれと同時の選択肢である複数の副ローブ制約が適用可能である。
【０１７８】
同様に、本発明のビームフォーマを、音声帯域用途よりも著しく高いまたは低い周波数に適用することができる。たとえば、通信および位置確認のためのハイドロホンアレイを備えたソナーシステムは、低周波数で動作するようになっている傾向があるものの、超音波ドランスデューサが通常５〜３０ＭＨｚの周波数範囲で動作するような超音波用途であっても本発明にかかるビームフォーマから恩恵を受ける。超音波ビームフォーミングは、多重指向性の迅速な選択および干渉の迅速な抑制により高い品質のイメージが提供されうる医療イメージング、断層撮影法などの用途において用いることができる。超音波は、患者のイメージングが呼吸および心臓の鼓動からの定常的な動きならびに不随意運動から影響を受ける実時間スピードおいて非常に効果的である。
【０１７９】
また、本発明は縦波音波の解析に限定されるものではない。ビームフォーミングは、センサーがアンテナである電磁波放射に対しても適用可能である。具体的にいえば、無線周波数用途において、レーダシステムはビームフォーミングから非常に恩恵を受けることができる。いうまでもなく、これらのシステムは、ビームパターンの実時間適応をさらに必要とする。たとえば相当の速度で移動している複数の航空機を追跡する場合、実時間において複数の主ローブを形成することは非常に有益である。
【０１８０】
さらに、本発明の用途には、地震探査、たとえば石油探査が含まれる。この分野では、ルック方向が非常に具体的であり、かつ正確であることが重要である。したがって、広大な土地を網羅しなければならないような場合、主ローブ幅および指向性制約を速く課することにより、このようなシステムが高速で動作することが可能となる。
【０１８１】
したがって、１つの好ましい実施形態では、本発明は、上述のようなビームフォーマを備えており、センサーアレイはアレイ状に並べられているハイドロホンである。
【０１８２】
他の好ましい実施形態では、本発明は、上述のようなビームフォーマを備えており、センサーアレイはアレイ状に並べられている超音波トランスデューサである。
【０１８３】
他の好ましい実施形態では、本発明は、上述のようなビームフォーマを備えており、センサーアレイはアレイ状に並べられているアンテナである。好ましい実施形態によっては、これらのアンテナは無線周波数アンテナである。
【０１８４】
いうまでもなく、本発明にかかるビームフォーマは、大部分がソフトウェとして実現され、このソフトウェアがコンピューティングデバイス（これは、たとえば一般的なパーソナルコンピュータ（ＰＣ）であってもよいし、もしくはメインフレームコンピュータであってもよい）により実行されるようになっていてもよいし、これが特別に設計されプログラムされたＲＯＭ（読み取り専用メモリ）であってもよいし、または、これがフィールドプログラマブルゲートアレイ（ＦＰＧＡ）で実現されてもよい。このようなデバイスでは、ソフトウェアは、あらかじめ搭載されていてもよいし、データキャリアによってもしくはネットワークを通じた転送によってシステムへ転送されてもよい。インターネットの如き広域ネットワークに接続されているシステムは、ソフトウェアの新バージョンをダウンロードして、それを更新するように構成されていてもよい。
【０１８５】
したがって、さらなる態様によれば、本発明は、ソフトウェア製品を提供し、このソフトウェア製品がコンピュータにより実行されると、コンピュータに上述の方法のステップを実行させる。ソフトウェア製品はデータキャリアであってもよい。それに代えて、ソフトウェア製品は、リモート位置から送信される信号を含んでいてもよい。
【０１８６】
本発明の他の態様によれば、物理的キャリアの形態を有しているソフトウェア製品を製造する方法は、データキャリアにインストラクションを格納することを含み、当該インストラクションがコンピュータにより実行されると、当該コンピュータが上述の方法を実行する。
【０１８７】
本発明のさらに他の態様によれば、ソフトウェア製品をリモート位置へ、当該リモート位置のコンピュータへデータを送信することによって送る方法は、データがインストラクションを含み、当該インストラクションがコンピュータによって実行されると、当該コンピュータが上述の方法を実行する。
【０１８８】
下記に、添付の図面を参照しながら、本発明の好ましい実施形態を例示のみを意図して記載する。
【図面の簡単な説明】
【０１８９】
【図１】選択値ζに関して、第一の実施形態にかかる位数Ｎ＝４のノルム制約球面アレイビームフォーマの指向性指数をｋａの関数として表すグラフである。
【図２】選択値ζに関して、第一の実施形態にかかる位数Ｎ＝４のノルム制約球面アレイビームフォーマのホワイトノイズ利得をｋａの関数として表すグラフである。
【図３】選択値ｋａに関して、第一の実施形態にかかる位数Ｎ＝４のノルム制約球面アレイビームフォーマの指向性指数をホワイトノイズ利得の関数として表すグラフである。
【図４Ａ】ｋａ＝３であり、アレイの位数がＮ＝４であり、マイクロホンの使用数が２５個である、遅延和型ビームフォーマの指向性パターンを示す図である。
【図４Ｂ】ｋａ＝３であり、アレイの位数がＮ＝４であり、マイクロホンの使用数が２５個である、純粋位相モードビームフォーマの指向性パターンを示す図である。
【図４Ｃ】ｋａ＝３であり、アレイの位数がＮ＝４であり、マイクロホンの使用数が２５個である、ロバスト最大ＤＩノルム制約ビームフォーマの指向性パターンを示す図である。
【図５】ζ＝Ｍ／４の場合の、ｋａ＝１、２および４に対応する周波数における、第一の実施形態にかかる遅延和型ビームフォーマおよびノルム制約型ビームフォーマの指向性パターンを高さの関数として示す図である。
【図６】ζ＝Ｍ／４であり、ｋａ＝３である、第二の実施形態にかかるノルム制約ビームフォーマの指向性パターンを示す図である。
【図７Ａ】ｋａ＝３であり、ＤＩが最大化されている、第三の実施形態にかかる副ローブ制御を備えたロバストビームフォーマの指向性パターンを示す図である。
【図７Ｂ】ｋａ＝３であり、深さが−４０ｄＢで幅が３０°のノッチが（６０°２７０°）の方向に形成されている、第三の実施形態にかかる副ローブ制御を備えたロバストビームフォーマの指向性パターンを示す図である。
【図７Ｃ】ｋａ＝３であり、ＳＮＲ出力が最大化され、（６０°、２７０°）の干渉の到来方向に対してヌルが形成されている、第三の実施形態にかかる副ローブ制御を備えたロバストビームフォーマの指向性パターンを示す図である。
【図８Ａ】副ローブ制御が均一であるロバストビームフォーミング用のビームパターンを示す図である。
【図８Ｂ】副ローブ制御が不均一であり、ノッチが形成されているロバストビームフォーミングのビームパターンを示す図である。
【図９Ａ】副ローブ制御および自動多重ヌルステアリングを備えたロバストビームフォーミングのビームパターンを示す図である。
【図９Ｂ】副ローブ制御、多重主ローブおよび自動多重ヌルステアリングを備えたロバストビームフォーミングのビームパターンを示す図である。
【図１０Ａ】副ローブ制御を備えていない単一ビームのビームパターンを示す図である。
【図１０Ｂ】副ローブ制御が不均一である単一ビーム用のビームパターンを示す図である。
【図１１Ａ】均一な副ローブ制御および適応ヌルステアリングを備えた単一ビームのビームパターンを示す図である。
【図１１Ｂ】副ローブ制御のない多重ビームのビームパターンを示す図である。
【図１２Ａ】副ローブ制御および適応ヌルステアリングを備えた多重ビームフォーミングのビームパターンを示す図である。
【図１２Ｂ】主ローブレベル制御を備えた多重ビームフォーミングのビームパターンを示す図である。
【図１３】ロバストネス制約が課され、副ローブ制御が行われていない、位数４の標準ビームパターンを示す図である。
【図１４】ロバストネス制約が課され、副ローブ制御制約が課されている、位数４の最適ビームパターンを示す図である。
【図１５】ロバストネス制約が課され、副ローブ制御が行われ、方向（５０、９０）から来る干渉対して深いヌルステアリングが施されている、位数４の最適ビームパターンを示す図である。
【図１６】該当する信号の方向に対する６つの無歪制約が課されている、最適多重主ローブビームパターンを示す図である。
【図１７】ヌルが（０、０）で形成され、た下半球に対して副ローブが制御される、該当する信号の方向に６つの無歪の制約を備えた最適多重主ローブビームパターンを示す図である。
【図１８Ａ】本発明にかかる方法とその方法を実行するための装置とを概略的に示すフローチャートである。
【図１８Ｂ】本発明にかかる方法とその方法を実行するための装置とを概略的に示すフローチャートである。
【図１９】テレコンファレンスというシナリオで本発明を実施したものを示す図である。
【図２０】周波数ドメインで動作し、ステアリングユニットを備えているモーダルビームフォーマの構造を概略的に示す図である。
【図２１】ステアリングユニットおよび複数のＦＩＲフィルタを備えている広帯域モーダルビームフォーマを時間ドメインで具象化したものを概略的に示す図である。
【図２２Ａ】ＦＩＲフィルタ係数を示している、最大ロバストネス設計を用いたモーダルビームフォーマの性能を示す図である。
【図２２Ｂ】最大ロバストネス設計を用いて時間ドメインビームフォーマおよび周波数ドメインビームフォーマの重み付け関数を周波数の関数として示している、最大ロバストネス設計を用いたモーダルビームフォーマの性能を示す図である。
【図２２Ｃ】周波数とアングルとの関数としてビームパターンを示している、最大ロバストネス設計を用いたモーダルビームフォーマの性能を示す図である。
【図２２Ｄ】さまざまな周波数でＤＩとＷＮＧを示している、最大ロバストネス設計を用いたモーダルビームフォーマの性能を示す図である。
【図２３Ａ】ＦＩＲフィルタの係数を示している、最大指向性設計を用いた時間ドメインモーダルビームフォーマの性能を示す図である。
【図２３Ｂ】重み付け関数を示している、最大指向性設計を用いた時間ドメインモーダルビームフォーマの性能を示す図である。
【図２３Ｃ】ビームパターンを示している、最大指向性設計を用いた時間ドメインモーダルビームフォーマの性能を示す図である。
【図２３Ｄ】さまざまな周波数でのＤＩおよびＷＮＧを示している、最大指向性設計を用いた時間ドメインモーダルビームフォーマの性能を示す図である。
【図２４Ａ】ロバスト最大指向性設計を用いてビームフォーマの性能を示す図である。
【図２４Ｂ】ロバスト最大指向性設計を用いてビームフォーマの性能を示す図である。
【図２４Ｃ】ロバスト最大指向性設計を用いてビームフォーマの性能を示す図である。
【図２４Ｄ】ロバスト最大指向性設計を用いてビームフォーマの性能を示す図である。
【図２５Ａ】２オクターブにわたって周波数が変動しないパターンを備えたビームフォーマの性能を示す図である。
【図２５Ｂ】２オクターブにわたって周波数が変動しないパターンを備えたビームフォーマの性能を示す図である。
【図２５Ｃ】２オクターブにわたって周波数が変動しないパターンを備えたビームフォーマの性能を示す図である。
【図２６Ａ】多重制約最適化を用いるビームフォーマの性能を示す図である。
【図２６Ｂ】多重制約最適化を用いるビームフォーマの性能を示す図である。
【図２６Ｃ】多重制約最適化を用いるビームフォーマの性能を示す図である。
【図２７Ａ】２つの典型的なマイクロホンで受信した時系列データおよびそのうちの一方のスペクトログラムからなる実験結果を示す図である。
【図２７Ｂ】ＴＤＭＲモーダルビームフォーマの２つの異なるステアリング方向の時系列出力およびそのうちの一方のスペクトログラムからなる実験結果を示す図である。
【図２７Ｃ】ＴＤＭＤモーダルビームフォーマの２つの異なるステアリング方向の時系列出力およびそのうちの一方のスペクトログラムからなる実験結果を示す図である。
【図２７Ｄ】ＴＤＲＭＤモーダルビームフォーマの２つの異なるステアリング方向の時系列出力およびそのうちの一方のスペクトログラムからなる実験結果を示す図である。
【発明を実施するための形態】
【０１９０】
まず図１８を参照すると、Ｍ個のマイクロホンからなる球状マイクロホンアレイ用のビームフォーミングシステムとして、本発明にかかるシステムの好ましい実施形態が概略的に示されている。
【０１９１】
複数のマイクロホン１０（図１８において概略的に示されているが、実際には球面アレイ状に配置される）は、それぞれアレイのまわりの環境から音波を受信し、これらを電気信号へ変換するようになしてある。Ｍ個のマイクロホンの各々からの信号は、ステージ１１において、まずＭ個のプリアンプ、Ｍ個のＡＤＣ（アナログ−ディジタル変換器）およびＭ個の校正フィルタによって処理される。次いで、これらの信号はすべてステージ２０へ送信され、ここで、データが高速フーリエ変換アルゴリズムによりＭ個の周波数ビンチャネルへと分割される。次いで、これらはステージ１２へ送信され、ここで、球面フーリエ変換（ｓｐｈｅｒｉｃａｌ  Ｆｏｕｒｉｅ  ｔｒａｎｓｆｏｒｍ）が行われる。ここでは、信号が位数Ｎの球面調和関数ドメインに変換される。すなわち、位数（ｏｒｄｅｒ）ｎ＝０、．．．、Ｎおよび次数（ｄｅｇｒｅｅ）ｍ＝−ｎ、．．．、ｎの（Ｎ＋１）2個の球面調和関数の各々に対して球面調和関数係数が生成される。
【０１９２】
球面調和関数ドメイン情報が、制約定式化のためにステージ１３へ送られ、最適化後ビームパターン合成のためにステージ１６へと送られる。ステージ１３では、システムの所望のパラメータが、調整可能パラメータステージ１４から入力される。図示されているように、入力することができる所望のパラメータには、信号のルック方向、主ローブ幅（１４ａ）、ロバストネス（１４ｂ）、所望の副ローブレベルおよび副ローブ領域（１４ｃ）ならびに所望のヌルロケーションおよびヌル深さ（１４ｄ）が含まれる。
【０１９３】
ステージ１３では、ステージ１２から球面調和関数ドメイン信号情報と組み合わせられた、ビームパターンのための所望の入力パラメータが受け取られ、これらを、凸最適化に適した凸二次最適化制約式（ｃｏｎｖｅｘ  ｑｕａｔｒａｔｉｃ  ｏｐｔｉｍｉｚａｔｉｏｎ  ｃｏｎｓｔｒａｉｎｔｓ）へと定式化される。制約は、自動ヌルステアリング、主ローブ制御、副ローブ制御およびロバストネスに関して定式化される。次いで、これらの制約は、ステージ１５へ送られる。ステージ１５は、内点法または二次錐計画法の如き数値最適化アルゴリズムを実行するための凸最適化ソルバーであり、ここで、入力される制約の下で最適なビームパターンを提供するために球面調和関数係数に適用される最適重み付け係数が求められる。特筆すべき点は、スペースドメインでは、球面調和関数ドメインへの変換が行われず、最適化された重み付け係数が入力信号に直接適用されるということである。
【０１９４】
次いで、これらの求められた重み付け係数がステージ１６へ送られる。ステージ１６では、これらの係数がステージ１２からのデータと加重和（ｗｅｉｇｈｔｅｄ  ｓｕｍ）として組み合わせられ、ステージ１７で、最終的に、単一チャネル逆高速フーリエ変換が実行され、アレイ出力信号が形成される。
【０１９５】
ここで、本発明の現実的な実施を説明する。図１９では、本発明がテレコンファレンスを想定して実施されている。２つの会議室３０ａおよび３０ｂが示されている。各部屋には、三次元で音声を取り込むための球面マイクロホンアレイ３２ａ、３２ｂと、一組のラウドスピーカ３４ａ、３４ｂとを備えているテレコンファレンス用のシステムが装備されている。各部屋には、その角に４つのスピーカが設けられているように示されているが、いうまでもなく、他の構成も同様に有効である。また、各部屋は、マイクロホンアレイのまわりのさまざまな位置に、３人の話者３６ａ、３６ｂが存在していることが示されている。マイクロホンアレイは、ビームフォーマおよびそれに付随するコントローラ３８ａ、３８ｂに接続されている。ビームフォーマおよびそれに付随するコントローラは、マイクロホンアレイ３２ａ、ｂのために最適なビームパターンを生成するための最適化アルゴリズムを実行するようになしてある。
【０１９６】
動作時、３人の話者３４ａのうちの一人が話し、その他の人が黙っていると想定する。コントローラ３８ａがソース信号を検出し、ビームフォーマを制御して部屋３０ａ内のマイクロホンアレイ３２ａのビームフォーミングパターンを生成することにより話者３６ａの方向に主ローブ（すなわち、高利得領域）を形成し、かつ、その他の方向に対するアレイ利得を最小限に抑える。
【０１９７】
部屋３０ｂでは、ビームフォーマ３８ｂは、ラウドスピーカ３４ｂの各々からの音響ソースを干渉ソースとして検出する。これらの方向からの音を最小限に抑えることは、２部屋間のフィードバックループの回避に望ましいことである。
【０１９８】
ここで、部屋３０ｂ内の話者３６ｂが部屋３０ａ内の人と話し始めると、部屋３０ｂ内のビームフォーマは、直ちにその話者の方向に沿って主ローブを形成して彼または彼女の声が部屋３０ａまで安全に送られることを担保するようにしなければならない。同様に、部屋３０ａ内のビームフォーマ３８ａは、部屋３０ｂのフィードバックを回避するために、ビームパターンにおいて、ラウドスピーカ３４ａの方向に対して深いヌルを直ちに形成しなければならない。
【０１９９】
ビームフォーマ３８ａ、３８ｂが複数の主ローブおよび複数の深いヌルを形成することができ、実時間中にこれらの指向性を制御することができるので、話者のうちの一人が会話中に部屋の中を歩き回り始めたとしても、このシステムが故障してしまうことはない。オフィスを通り過ぎていく警察のサイレンの如き予期しない干渉であっても、実時間中に深いヌルの指向性を制御することにより対処することができる。同時に、ビームフォーマ３８ａ、３８ｂは、ビルディングの空調機の送風機の如き一般的な背景ノイズの影響を最小限に抑えるよう、適用される制約条件の境界内にアレイ出力を抑えるようになしてある。
【０２００】
このシステムは、全二重送信、騒音削減、残響除去、およびアコースティックエコーキャンセレーションなどの機能を備えた高品質な３Ｄ空間オーディオ機器を提供する。
Ａ．特別なケース
次に上記の最適化問題（３２）のいくつかの特別なケースを考え、これらを従来の研究結果と比較する。
特別なケース１
指向性が最大であり、ＷＮＧまたは副ローブ制御がない。このことは、（２４）においてε＝０、ζ＝０、
【数１１５】
2012523731.tif 000117
、Ｑ（ω）＝Ｑiso（ω）とすることにより定式化される。このことにより、Ｒ（ω）＝Ｑiso（ω）となり、また、（３２）内の２つの不等式制約が常にアクティブではない（制約がない）ので無視することができる。
【０２０１】
指向係数を等方性ノイズに対するアレイ利得として解釈することができるので、この場合の最適化問題は指向係数を最大化することになる。
【０２０２】
この場合の最適化問題は、古典的なアレイ処理におけるカポン（Ｃａｐｏｎ）ビームフォーマに類似し、（３２）の解は容易に導出される。
【０２０３】
【数１１６】
2012523731.tif 000118
（７）および（２６）を用いて、さらに下記の式
【０２０４】
【数１１７】
2012523731.tif 000119
を用いると、式（３３）を以下の式にさらに変換することができる。
【０２０５】
【数１１８】
2012523731.tif 000120
この式で、
【数１１９】
2012523731.tif 000121
は要素毎の割り算を行うことを表している。すなわち
【数１２０】
2012523731.tif 000122
のことを表している。（３５）内の重み付けは、アレイ利得には影響を与えることのないスカラー乗数を除いて、純粋位相モード球面マイクロホンアレイの重み付けと同一である（たとえば下記の文献を参照されたい：、ビー・ラファリー（Ｂ．ラファリー）、「位相モード対遅延和型球面マイクロホンアレイ処理（Ｐｈａｓｅ−ｍｏｄｅ  ｖｅｒｓｕｓ  ｄｅｌａｙ−ａｎｄ−ｓｕｍ  ｓｐｈｅｒｉｃａｌ  ｍｉｃｒｏｐｈｏｎｅ  ａｒｒａｙ  ｐｒｏｃｅｓｓｉｎｇ）」、ＩＥＥＥ信号プロセスレター（ＩＥＥＥ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓ．Ｌｅｔｔ．）、第１２巻、１０号、ページ７１３〜７１６、２００５年１０月（本明細書の最初のページでも引用されている）。
【０２０６】
（３１）および（２８）において（３５）を用いると下記の式が得られる：
【０２０７】
【数１２１】
2012523731.tif 000123
【０２０８】
【数１２２】
2012523731.tif 000124
（特筆すべき点は、これらが、先に引用されたラファリーの文献においてｄn≡１である場合の（１１）および（１２）と同一であることである。）この結果は、位数Ｎの純粋位相モード球面マイクロホンアレイが周波数に依存しない２０ｌｏｇ10（Ｎ＋１）ｄＢの最大ＤＩを有することを裏付けることになる。
特別なケース２  
ＷＮＧが最大であり、指向性おおび副ローブ制御がない。これは、次のように定式化される。Ｒ（ω）＝Ｉ、この式で、Ｉは恒等行列であり、ε＝∞であり、ζ＝０である。
【０２０９】
明らかに、このケースの最適化問題は、重み付けベクトルのノルムを最小化すること、またはホワイトノイズ利得を最大化することになる。
【０２１０】
（３３）内のＱisoがＩにより交換されると、このケースの解は下記の式となる：
【０２１１】
【数１２３】
2012523731.tif 000125
【０２１２】
【数１２４】
2012523731.tif 000126
これは、開球構造の場合、スカラー乗数を除いて、遅延和型の球面マイクロホンアレイの重み付けと同一である。
【０２１３】
さらに、（３１）および（２８）において（３８）を用いると、下記の式が得られる。
【０２１４】
【数１２５】
2012523731.tif 000127
【０２１５】
【数１２６】
2012523731.tif 000128
（これが上述のラファリーの文献内の（１７）および（１８）と同一の結果であることに留意されたい）。
【０２１６】
Ｎ→∞の場合、（４０）内の総和が（４π）2に接近するので、遅延和型アレイは、周波数に依存しない定数である、Ｍに等しいＷＮＧを達成する。これは、古典的アレイ処理において周知の結果である。
特別なケース３  
指向性およびＷＮＧが制御され、副ローブが制御されない。このケースは、ε＝０の基準によって定式化される。
【０２１７】
この場合最適化問題は、ホワイトノイズ利得制約（または、ノルム制約）ロバストカポンビームフォーミング問題に類似した形態を有する。
【０２１８】
ζ＝ＷＮＧ2である場合、対応する解は特別なケース２で説明されたような遅延和型アレイであるということを確認するのは容易なことである。さらに、Ｒ（ω）＝Ｑiso（ω）であり、ζの値を（０，ＷＮＧ2］の範囲で調製すると、純粋な位相モードと遅延和型球面アレイ処理との間のトレードオフを得ることができる。
【０２１９】
本発明にかかる下記の好ましい実施形態は、上述のビームフォーマのシミュレーションであり、その性能を説明および評価するために用いられる。図１〜図７のシミュレーションでは、位数Ｎ＝４の開球アレイを想定しており、また、マイクロホンの数がＭ＝（Ｎ＋１）2の数であると仮定している。
【０２２０】
本明細書に記載のシミュレーションはすべて、消費者等級のコンピュータ設備、たとえばＣＰＵ速度が２．４ＧＨｚで、ＲＡＭが２ＧＢであるノート型ＰＣを用いて行われたものである。これらのシミュレーションは、ＭＡＴＬＡＢを用いて行われ、各狭帯域シミュレーションついて約２〜５秒を要した。いうまでもなく、ＭＡＴＬＡＢコードは、数理解析およびシミュレーションのために設計されたハイレベルプログラミング言語であり、最適化アルゴリズムが、Ｃまたはアセンブリー言語の如き低レベルプログラミング言語に実装された場合、またはこれらがフィールドプログラマブルゲートアレイに実装された場合、その速度が著しく高まることが期待できる。
Ｂ．純粋な位相モードと遅延和型アレイとの間のトレードオフ
Ｒ（ω）＝Ｑiso（ω）であり、ε＝∞であるとする。最適化問題（３２）はノルム制約最大ＤＩビームフォーミング問題となる。球面アレイ構造は三次元対称性を提供する。一般性を失うことなく、ルック方向がΩ0＝［０°、０°］であると仮定する。ある与えられたζの値に対して、ｋａの関数としてこの最適化問題を解いて重みベクトルｗ（ｋ）を求め、それらを（２８）および（３１）に代入してＤＩおよびＷＮＧを得る。図１および図２は、それぞれＤＩおよびＷＮＧを、ζが０，Ｍ／２，Ｍ／４およびＷＮＧ2である場合ついて、ｋａの関数として示している。ζ＝０の場合およびζ＝ＷＮＧ2の場合は、それぞれ純粋位相モードアレイおよび遅延和型アレイに対応している。ζ＝Ｍ／２の場合およびζ＝Ｍ／４の場合は、理想的な最大ＷＮＧであるＭと比較して、それぞれＷＮＧが３ｄＢおよび６ｄＢだけ低下したロバストビームフォーマに対応している。
【０２２１】
図２は、ノルム制約ビームフォーマが、ある与えられたしきい値を超えるＷＮＧを生じ、良好なロバストネスを実現できることを示している。２つのノルム制約ビームフォーマのＤＩ、すなわちζ＝Ｍ／２およびζ＝Ｍ／４が、遅延和型ビームフォーマのものよりもはるかに高い。これらのＤＩが、純粋位相モードビームフォーマのものより小さいが、取得可能である。しかしながら、後者のものは、実社会での適用の際に遭遇するどのような小さなランダムアレイ誤差（ｒａｎｄｏｍ  ａｒｒａｙ  ｅｒｒｏｒｓ）に対しても非常に敏感であるために、取得可能でないことが多い。これに加えて、純粋位相モードビームフォーマに関する図２のｋａ＝３．１４およびｋａ＝４．５０において観察される２つのＷＮＧの値が非常に低いことについては、開球アレイにおける周知の問題であって、剛球アレイを用いることによりこの問題は回避される。要約すると、この具体例では、ノルム制約ビームフォーミングが純粋位相モードと遅延和型アレイとの間の有用なトレードオフを提供しうることが明示されている。
【０２２２】
さらに、ζ＝Ｍ／２およびζ＝Ｍ／４の場合、重み付けベクトルノルム制約がｋａ＝４およびｋａ＝５のあたりでインアクティブとなることが分かる。このことは、これらの領域あたりでは、純粋位相モードビームフォーマが既に相当なＷＮＧを提供しているという事実に起因している。したがって、これらの領域あたりでは、これらの２つのビームフォーマは純粋位相モードビームフォーマと同一となる。
【０２２３】
図３には、ノルム制約ビームフォーマのＤＩが、ｋａ＝１、２、３および４に対応する周波数について、ＷＮＧの関数として示されている。高周波数において、アレイが良好なＷＮＧ−ＤＩ性能を有していることが分かる。低周波数では、アレイのＷＮＧ−ＤＩ性能が著しく劣化している。
【０２２４】
ｋａ＝３に対応する周波数に関して、３つのビームフォーマ、すなわち遅延和型ビームフォーマ、純粋位相モードビームフォーマ、ノルム制約ビームフォーマの三次元アレイパターンが、ζ＝Ｍ／４の場合についで、（２３）により計算されている。図４に、これらの結果が示されており、この図には、ルック方向のパターンの振幅が１（または、０ｄＢ）に等しくなるように、正規化係数Ｍ／４πが含まれている。このケースの場合のアレイパターンがルック方向に対して対称となっていることが分かる。また、ノルム制約ビームフォーマが遅延和型ビームフォーマに比べて狭い主ローブを形成することも分かる。また、これらのビームフォーマのＤＩおよびＷＮＧの値も複数の図に示されている。図４（ｃ）内のＷＮＧは正確に１０ｌｏｇ10（Ｍ／４）＝７．９６ｄＢとなっている。
【０２２５】
図５では、ζ＝Ｍ／４の場合に関して、遅延和型（ＤＡＳ）ビームフォーマおよびノルム制約ビームフォーマの指向性パターンが、高さの関数として、ｋａ＝１、２および４に対応する周波数について比較されている。 純粋位相モードビームフォーマの指向性パターンが周波数に依存しないこと、図２により示唆されているように、ｋａ＝４においてζ＝Ｍ／４である場合に、ノルム制約ビームフォーマのものと同一であることは注目に値する。
Ｃ干渉が除去されるロバストビームフォーミング
上述の特別なケース３を考える。ノイズが等方性ノイズであると仮定する。信号および干渉が（０°、０°）および（−９０°、６０°）からアレイに入り、各センサーの信号（干渉）対ノイズ比がそれぞれ０ｄＢおよび３０ｄＢであと仮定する。また、正確な共分散が、知られており、理論的なアレイ共分散行列Ｒ（ω）（２４）によって表されると仮定する。
【０２２６】
このケースの場合、最適化問題は、ノルム制約ロバストカポンビームフォーミング問題となり、指向性の低下を犠牲にして高アレイ利得を得るビームフォーマに帰着する。
【０２２７】
図６は、ζ＝Ｍ／４およびｋａ＝３という値の場合に得られるアレイパターンが示されている   予想通り、アレイパターンは、干渉の到着方向に対してヌルが深くなっている。このケースの場合のアレイパターンは、図４に示されている純粋位相モードビームフォーマおよび遅延和型ビームフォーマによるものとは異なり、ルック方向に対してもはや対称ではなくなっている。
Ｄ．副ローブが制御され、干渉が除去されるロバストビームフォーミング
図４および図６には、ｋａ＝３におけるこれらのアレイパターンの副ローブレベルがほぼ−１３．２ｄＢから−１６．３ｄＢであることが示されている。このような値は、ほとんどの用途では高すぎて、予期しないまたは急に生じる干渉に対する性能が著しく劣化してしまう恐れがある。このような状況での用途のために副ローブが制御されるビームフォーマについでの具体例を検討する。
【０２２８】
Ｒ（ω）＝Ｑiso（ω）の場合の等方性ノイズをまず仮定し、ｋａ＝３であり、ζ＝Ｍ／４であり、ε＝０．１であるケース、すなわち所望の副ローブレベルが−２０ｄＢであるケースを考える。 副ローブ領域が、ΩSL＝｛（θ、φ）｜θ≧４５°｝として定義される。（３２）の最適化問題の解は、副ローブが制御されるノルム制約最大ＤＩビームフォーマである。図７（ａ）には、得られたアレイパターンが示されている。指定されている副ローブレベルは−２０ｄＢ未満である。
【０２２９】
ここで、副ローブ制御に加えて、深さが−４０ｄＢでありかつ幅が３０°である、（６０°、２７０°）方向のノッチを形成することを考える。このケース場合、所望の副ローブ構造は方向に依存する。図７（ｂ）には、所望のノッチ領域にε＝０．０１を設定し、それと同時にその他の副ローブ領域ではε＝０．１を維持して、最適化問題を解くことによって得られたアレイパターンが示されている。所定のノッチが形成され、−２０ｄＢという低い副ローブレベルが維持されていることが分かる。
【０２３０】
上述のＣに記載されているシナリオを検討する。副ローブを−２０ｄＢ未満に、すなわちε＝０．１となるように制御したいと仮定する。その他のパラメータをＣで用いられたものと同一のままにする。ビームフォーマ重み付けベクトルを最適化問題（３２）を解くことにより求める。図７（ｃ）には、得られたアレイパターンが示されている。図４（ａ）と比較して、干渉の入ってくる方向に対してヌルであることに加えて、この方法による副ローブが厳密に−２０ｄＢ未満であることが分かる。
【０２３１】
剛球アレイの下記のシミュレーションでは、位数がＮ＝４であり、複数の主ローブ制約が課され、非均一副ローブ制約が課されている。ビームパターンに複数の主ローブを形成するためには、該当する方向に、非ひずみ制約が課されなければならない。非均一副ローブ制御の場合、副ローブ領域のすべての標本点（ｓａｍｐｌｉｎｇ  ｐｏｉｎｔｓ）がある与えられたしきい値未満であることを必要とすることに代えて、副ローブ方向の各々に異なるしきい値を課すことができるようにすることができる。たとえば、干渉方向に対しては強い制約を課すようにし、その他の方向にはそれよりも弱いしきい値を課すようにすることができる。これらのさらなる制約（Ｋ個の主ローブ制約およびＬ個の副ローブ制約）を用いると、最適化問題（３２）は次のよう書き換えることができる：
【０２３２】
【数１２７】
2012523731.tif 000129
先の場合と同様に、この最適化問題の特性により、凸最適化を適用することができ、とくにこの最適化問題が凸二次錐計画法問題であるので、ＳＯＣＰ技術を用いてその問題を解くことができる。これらの技術を用いると、多数の制約が含まれることになったとしても、この問題を今でどおり効率的に実時間で最適化することができる。
【０２３３】
このビームフォーマの性能を評価するために、さらなるシミュレーションが用いられる。位数がＮ＝４であり、Ｍ＝（Ｎ＋１）2である剛球アレイを考える。単一の主ローブの場合についてのルック方向が［０°、０°］であり、ｋａ＝３であり、各センサーの信号対ノイズ比および干渉対ノイズ比が０ｄＢおよび３０ｄＢであり、ＷＮＧ制約が８ｄＢに設定されていると仮定する。図８（ａ）には、副ローブ領域がΩSL＝｛（θ、φ）｜θ≧４５°｝として定義され、副ローブレベルが−２０ｄＢ未満であるアレイパターンが示されている。図８（ｂ）には、非均一副ローブ制御の性能が示されている。（６０°、２７０°）の方向に、深さが−４０ｄＢで、幅が３０°のノッチが形成され、また、その他の副ローブレベルが今までと同様に−２０ｄＢで維持されている。
【０２３４】
図９（ａ）では、２つの干渉が（６０°、１９０°）および（９０°、２６０°）からアレイに入り、次いで、副ローブが−２０ｄＢ未満に厳密に抑えられた状態で、ヌルが、干渉の入ってくる方向に対して自動的に形成・方向付けされていることが分かる。図９（ｂ）には、複数の主ローブ形成の性能、副ローブ−が２０ｄＢに制御された状態での自動ヌルステアリングの性能が示されており、この図では、３つの干渉が（０°、０°）、（４５°、９０°）および（５０°、２７０°）からアレイに入ってくる状態で、２つの所望の信号が（４０°、０°）および（４０°、１８０°）からアレイに入ってくることが仮定されている  また、図８および図９に関して、実際の指向性指数（ＤＩ）およびＷＮＧ値が計算されている。
【０２３５】
次の解析では、部屋に設けられたコンパクトな球面マイクロホンアレイが想定されている。信号ソースがすべてアパーチャの遠方界に設けられていると仮定し（その結果、これらの信号ソースをアレイ上に入射する平面波で近似しうる）、部屋内での初期段階での反射は点光源としてモデル化され、遅い段階ｄでの反射は等方性ノイズとしてモデル化される。ここで、Ｌ＋Ｄ個のソース信号がΩ1Ω2、．．．、ΩL、ΩL+1、．．．、ΩL+Dの方向からの球面に入って来ること、また、ノイズが存在することを仮定する。次いで、各マイクロホン位置のスペースドメイン音圧は次のように書くことができる：
【０２３６】
【数１２８】
2012523731.tif 000130
この式で、
【数１２９】
2012523731.tif 000131
がＬ＋Ｄ個の信号スペクトルであり、
【数１３０】
2012523731.tif 000132
および
【数１３１】
2012523731.tif 000133
がＲ個の初期段階での反射であり、αおよびτが初期段階での反射の減衰時間および伝播時間を表し、Ｎ（ω、Ωs）が追加のノイズスペクトルである。（４３）内の第一の項は、取り込むことが望まれるＬ個の所望の信号に相当し、（４３）内の第二の項はＤ個の干渉に相当する。
【０２３７】
ｘ（ｋａ、Ωs）の球面フーリエ変換は次の式で与えられる：
【０２３８】
【数１３２】
2012523731.tif 000134
この式で、Ｎnm（ω）はノイズの球面フーリエ変換であり、先の場合のように、Ｎは、Ｍ≧（Ｎ＋１）2を満たす球面調和関数の位数（ｏｒｄｅｒ）である。
【０２３９】
次いで、アレイ処理をスペースドメイン内または球面調和関数ドメイン内で行なうことができ、またアレイ出力ｙ（ｋａ）が以下の式で計算される：
【０２４０】
【数１３３】
2012523731.tif 000135
先の場合のように、αsはサンプリング手法によって異なる。サンプリングが均一な場合、αs＝４π／Ｍである。
【０２４１】
一部の実施形態の場合と同様に、次の実施形態にかかるビームフォーマでは、複数の主ローブ方向が維持され、副ローブレベルが制御されるが、ビーム方向外からくる干渉を順応的に抑制するためにアレイ出力が最小限に抑えられるようになしてある。さらに、システムのロバストネスを向上させる目的で、重み付けノルム制約（すなわち、ホワイトノイズ利得制御）も用いられ、アレイの重み付けノルムをある選択されたしきい値に制限するようになしてある。
【０２４２】
Ωl＝Ω1、Ω2、．．．、ΩLの方向からくるＬ個の所望の信号が良好に取り込まれ、均一にされることを担保するために、以下のＬ×（Ｎ＋１）2マニホルド行列と、
【０２４３】
【数１３４】
2012523731.tif 000136
Ｌ個の所望の主ローブレベルを含んでいる以下のＬｘ１ベクトル列とが定義されている。
【０２４４】
【数１３５】
2012523731.tif 000137
この式で、４π／Ｍは正規化係数である。次いで、主ローブレベルが扱い易い多重ビームを形成する問題は、次のような単一の一次等式制約として定式化することができ
【０２４５】
【数１３６】
2012523731.tif 000138
Ｌ個の主ローブ応答レベルは異なるＡの値を設定することにより制御することができる。このことは、スピーチレベルが異なるＬ個の所望のスピーカの音声振幅を等化させるような単純な用途においてとくに有益なものとなる。このことは、主として部屋の異なる位置に話者が座るという事実に起因している。
【０２４６】
上述の実施形態と同様に、すべての副ローブをしきい値未満に厳格に抑えることを担保するために、次のような一組の二次不等制約を定式化することができる。
【０２４７】
【数１３７】
2012523731.tif 000139
この式で、ΩSL、jは副ローブ領域を表しており、これらは、複数の主ローブのビーム幅を制御するために用いられる。
【０２４８】
上記の実施形態でのように、適応する主ローブの形成および多重のヌルステアリングは、さまざまな制約を適用しながら、実行時間においてアレイ出力を最小限に抑えることにより達成される。（２２）において先に述べたように、アレイ出力は次の式で求められる：
【０２４９】
【数１３８】
2012523731.tif 000140
この式で、Ｅ［・］は統計的期待値を表し、Ｒ（ω）は共分散行列ｘを表している。 単純化するため、ルーム内の初期の反射が直接音よりもはるかに低いと仮定すると、Ｒ（ω）は次の式で表される：
【０２５０】
【数１３９】
2012523731.tif 000141
この式で、Ｒa（ω）はａ番目の信号に対応する信号分散行列であり、Ｒn（ω）はノイズ共分散行列である。
【０２５１】
ここで、変数ξの導入によって、最適化問題を次のように再定式化することができる：
【０２５２】
【数１４０】
2012523731.tif 000142
単一の主ローブについて（３１）で先に導出された重み付けベクトルノルム制約が多重主ローブの場合にも当てはまる。というのは、それがアレイの重み付けのダイナミックレンジを調節してアレイ出力においてノイズを大きく増幅することを回避するからである。
【０２５３】
これを（４６）、（４７）および（５０）と組み合わせることにより、（３２）の最適化問題を次のように表すことができる：
【０２５４】
【数１４１】
2012523731.tif 000143
このようにして、主ローブレベルが異なり、多重ヌル形成および多重ヌルステアリングを用いて副ローブが制御され、ロバストネス制約が課される多重主ローブ形成を実現する単一最適化問題が定式化されている。さらに、この最適化問題は、凸二次錐最適化問題であるので、実時間で、二次錐計画法を用いて効率的に解くことができるものである。
【０２５５】
上述の式から分かるように、重み付けベクトルノルム制約は、ζを分母に置いているのではなく、しきい値δを分子に置いて表されている。次のシミュレーションは、用いられているδの値が示されている。
【０２５６】
次のシミュレーションでは、ｒ＝５ｃｍの剛球がＭ＝（Ｎ＋１）2個のマイクロホンによりサンプリングされ、また、ｋａ＝３であると想定する。各マイクロホンの信号対ノイズ比および干渉対ノイズ比はそれぞれ０ｄＢおよび３０ｄＢである。副ローブ領域を離散化するために、５°の均一なグリッドが用いられる。他に指定がない限り、便宜上、適応ビームフォーミングの具体例では、理論的なデータ共分散行列Ｒ（ω）が用いられる。
【０２５７】
単一ビーム（Ｌ＝１）の場合、位数がＮ＝４であり、ルック方向が［０°、０°］であり、ＷＮＧ制約が８ｄＢ（δ＝０．１５９）に設定されていると仮定する。 図１０（ａ）には、副ローブが制御されず、適応ヌルステアリング制約がない、（５１）を用いた標準的な単一ビームパターンの合成が示されている。図１０（ｂ）には、不均一副ローブ制御の性能が示されている。主要な副ローブ領域が、副ローブレベルが均一に−２０ｄＢ）未満であり（εj＝０．０１）、ΩSL＝｛（θ、φ）｜θ≧４５°｝であるとして定義されている。ノッチが、（６０°、２７０°）の方向であり、深さが−４０ｄＢ（εj＝０．０００１）であり、幅が３０°であると定義されている。図１１（ａ）では、ノッチが取り除かれ、２つの干渉が［６０°、１９０°］および［９０°、２６０°］からアレイに入ってくることが仮定され、ヌルが自動的に形成され、干渉の入ってくる方向に向けられ、副ローブが−２０ｄＢ未満に厳密に維持されていることが分かる。すべての単一ビームのケースについて実際のＷＮＧおよび指向性指数（ＤＩ）の値が計算されていることに留意されたい。
【０２５８】
図１０（ｂ）では、主ローブが少し幅広くなっており、また、ＤＩが、副ローブが制御されないものよりも０．３ｄＢだけ低くなっている。しかしながら、これらの犠牲は実用化の際には受け入れ可能である。劣化の理由は、ビームフォーミング性能パラメータ、すなわちビーム幅、副ローブレベル、ＤＩおよびロバストネスがすべて相互に関連しているからである。本明細書に記載のアルゴリズムは、これらの相反する複数の目的間において適切な妥協点を提供する。
【０２５９】
多重ビームの例（Ｌ＝３）の場合、より多くの自由度を得るためにＮ＝５のアレイ位数（ａｒｒａｙ  ｏｒｄｅｒ）が用いられている。３つの所望の信号が［６０°、０°］［６０°、１２０°］および［６０°、２４０°］からアレイに入ってくると仮定されている。 図１１（ｂ）には、Ａ1、2、3＝１であり、δ＝４である多重ビームの形成性能が示されている。図１２（ａ）には、［０°、０°］［６５°、６０°］、［６５°、１８０°］および［６５°、３００°］から干渉が入ってくると仮定し、適応ヌルステアリングが行われ、−２０ｄＢの副ローブ制御が行われる多重ビーム受理可能性能が示されている。次に、第二の所望の信号の振幅が他の２つの信号よりも６ｄＢだけ低いと想定し、音レベルを同等にするために、Ａ2＝２に設定し、δ＝１に設定する。図１２（ｂ）には、ビームパターンが示されており、第二の主ローブ方向から入ってくる信号に対して約６ｄＢの振幅増強が得られることが示されている。
【０２６０】
図１３〜図１７には、本発明にかかる最適ビームフォーマの利点を示すさらなるシミュレーションが示されている。図１３は、ロバストネス制約は行われているが副ローブ制御は行われていない四位数（４ｔｈ  ｏｒｄｅｒ）標準ビームパターンが示されている。対照的に、図１４には、ロバストネス制約および副ローブ制御制約が行われている本発明にかかる四位数最適ビームパターンが示されている。主ローブは＋Ｚ軸から４５°の領域にある。図１５には、ロバストネス制約が行われ、副ローブ制御が行われ、方向（５０、９０）から入ってくる干渉に対して深いヌルステアリングが行われている、本発明に従って形成された四位数最適ビームパターンが示されている。
【０２６１】
図１６には、該当する信号の方向に対する６つの無歪制約が行われてビームパターンに６つの主ローブが形成されている、本発明に従って形成される最適多重主ローブビームパターンが示されている。図１７には、該当する信号の方向に対する６つの無歪制約が行われ、（０、０）でヌルが形成され、下半球側にある副ローブが制御されている、本発明に従って形成された最適多重主ローブビームパターンが示されている。
時間ドメインの例
下記には、広帯域モーダルビームフォーマのためのアレイパターン合成に対する時間ドメインアプローチの性能を例示するためにいくつかの数値の例が説明されている。
【０２６２】
下記で検討する具体例では、半径が４．２ｃｍであり、Ｍ＝３２個のマイクロホンが切頭２０面体の面の中央に設けられている剛球アレイが想定されている。音場分解のために位数Ｎ＝４が用いられ、αs≡４π／Mである。サンプリング周波数はｆs＝１４７００Ｈｚである。周波数バンド［ｆL、ｆU］は、Ｋ＝５１個の周波数グリッド、ｆk＝ｆL・１０lg(fU/fL)*(k-1)/(K-1)、ｋ＝１，２，．．．、Ｋを用いて離散化される。ＦＩＲフィルタの長さはＬ＝６５である。とくに明記されていない限り、ΘML＝［０°：２°：４０°］およびΘSL＝［４８°：２°：１８０°］であることが仮定され、このことは、方向を離散化するために２°の均一なグリッドが用いられることを意味する。
Ｔ．Ａ．最大ロバストネス設計
式（Ｔ４２）を参照し、ｆL＝５００Ｈｚ、ｆU＝５０００Ｈｚであると仮定する。ｌ＝４、μ1＝∞、μ2＝∞、μ3＝∞とする。最適化問題は次の式で表される。
【０２６３】
【数１４２】
2012523731.tif 000144
この問題の解は時間ドメイン最大ロバスト（ＴＤＭＲ）モーダルビームフォーマと呼ばれる。ＦＩＲフィルタｈは、最適化問題（Ｔ４３）を解くことにより求められ、そのサブベクトルｈ0、ｈ1、．．．、ｈNは、図２２（ａ）に示されている。ｈが（Ｔ２３）に代入され、
【数１４３】
2012523731.tif 000145
が求められ、それらが図２２（ｂ）に示されている。この図には、比較のため、（Ｔ１７）を用いて計算された［ｃn（ｆk）］MWNGも示されている。時間ドメイン最大ロバストモーダルビームフォーマの重み付け
【数１４４】
2012523731.tif 000146
は、周波数バンド［ｆL、ｆU］内において、周波数ドメイン最大ＷＮＧモーダルビームフォーマの重み付けに近似する。
【０２６４】
（Ｔ２５）を用いて、ビームパターンが、周波数および角度の関数として計算され、グリッド上の周波数および角度の点で表されている。図２２（ｃ）には、得られたビームパターンが示されている。この図には正規化係数Ｍ／４πが含まれているので、ルック方向のビームパターンの振幅は１（または０ｄＢ）に等しい。
【０２６５】
ＤＩおよびＷＮＧがそれぞれ（Ｔ３８）および（Ｔ１５）を用いて計算される。比較のために、周波数ドメイン最大ＷＮＧモーダルビームフォーマのＤＩおよびＷＮＧも計算される。図２２（ｄ）は、さまざまな周波数で得られた結果が示されている。
Ｔ．Ｂ．最大指向性設計
ｌ＝１，μ2＝∞、μ3＝∞、μ4＝∞とする。最適化問題（Ｔ４２）は最大指向性設計問題となる。得られたビームフォーマは、時間ドメイン最大指向性（ＴＤＭＤ）モーダルビームフォーマと呼ばれる。
【０２６６】
ｆL＝５００Ｈｚ、ｆU＝５０００Ｈｚであると仮定する。図２３（ａ）、図２３（ｂ）、図２３（ｃ）および図２３（ｄ）には、それぞれ得られたＦＩＲフィルタｈ0、ｈ1、．．．、ｈN、重み付け関数
【数１４５】
2012523731.tif 000147
、ビームパターン、ならびにＤＩおよびＷＮＧが示されている。これらの図には、比較のため、周波数ドメイン最大ＤＩモーダルビームフォーマの重み付け関数［ｃn（ｆk）］MDI（Ｔ１６）、ＤＩおよびＷＮＧがさらに示されている。最大指向性設計を用いた時間ドメインモーダルビームフォーマの重み付けが、周波数バンド［ｆL，ｆU］内において、周波数ドメインのものに近似することが分かる。
【０２６７】
図２２（ａ）、図２２（ｂ）、図２２（ｄ）と比較すると、ＦＩＲフィルタの係数、ひいてはＴＤＭＤビームフォーマの得られた重み付け関数が非常に大きく、また、低周波数におけるＷＮＧが小さすぎることが分かる。これらはすべて、このビームフォーマがロバストネスを欠いていることを示唆している。
Ｔ．Ｃ．ロバストネス制御を伴う最大指向性
ビームフォーマのロバストネスを向上させるためには、広帯域ホワイトノイズ利得制約を課す必要がある。このことは、ｌ＝１，μ2＝∞、μ3＝∞、μ4がユーザパラメータであるとして定式化することができる。得られたビームフォーマは、時間ドメインロバスト最大指向性（ＴＤＲＭＤ）モーダルビームフォーマと呼ばれる。
【０２６８】
ｆL＝５００Ｈｚ、ｆU＝５０００Ｈｚ、μ4＝４π／Ｍであると仮定する。 図２４（ａ）、図２４（ｂ）、図２４（ｃ）および図２４（ｄ）には、それぞれ、得られたＦＩＲフィルタｈ0、ｈ1、．．．、ｈN、重み付け関数
【数１４６】
2012523731.tif 000148
、ビームパターン、ならびにＤＩおよびＷＮＧが示されている。
【０２６９】
図２４（ｄ）から、このビームフォーマのＷＮＧが−３ｄＢよりも高く、低周波では、図２３に示されるような最大指向性設計のものよりもはるかに高いことが分かる。このビームフォーマのＤＩは、図２２に示されるような最大ロバストネス設計のものよりもはるかに高い。したがって、これらの結果から、この設計が指向性とロバストネスとの間の良好なトレードオフを提供していることが分かる。
Ｔ．Ｄ．周波数に依存しないビームフォーマ
周波数に依存しない広帯域ビームパターンの合成を仮定する。帯域幅を２オクターブ削減するので、ｆL＝１２５０Ｈｚ、ｆU＝５０００Ｈｚとなる。ｌ＝１、μ2＝１０-1.5・４π／Ｍ、ｑ1＝２，μ3＝∞、μ4＝２π／Ｍ、ΘML＝［０°：２°：１８０°］とする。図２５には、結果が示されている。期待された、周波数に依存しないビームパターンが得られ、ＷＮＧが適度であることが分かる。
Ｔ．Ｅ．多数の制約を伴う最適ビームフォーマ
ｆL＝１２５０Ｈｚ、ｆU＝５０００Ｈｚであると仮定する。ｌ＝１、μ2＝０．１・４π／Ｍ、ｑ1＝２，μ3＝１０-14/20・４π／Ｍ、ｑ2＝∞、μ4＝１０-4/10・４π／Ｍ、ΘML＝［０°：２°：４０°］、ΘSL＝［４８°：２°：１８０°］とする。図２６には得られた結果が示されている。制約がすべて担保され、また、複数の性能尺度間のトレードオフが得られている。
実験結果
ＭＨアコースティックスからのＥｉｇｅｎｍｉｋｅ（登録商標）マイクロホンアレイは、半径が４．２ｃｍの剛球アレイであり、３２個のマイクロホンが、切頭２０面体の面の中心に設けられている。実験は、７５Ｈｚまで無響になっている無響室で行なわれた。また、Ｅｉｇｅｎｍｉｋｅ（登録商標）は録音用のために無響室の中心に配置された。ラウドスピーカが、ほぼ（２０°、１８０°）の方向に向かってＥｉｇｅｎｍｉｋｅ（登録商標）マイクロホンアレイから１．５メートル離して設けられ、掃引周波数コサイン信号（１００Ｈｚ〜５ｋＨｚまでの範囲）を生じさせるために用いられた。音は、１４．７ｋＨｚのサンプリング周波数で、サンプル当たり１６ビットで、Ｅｉｇｅｎｍｉｋｅ（登録商標）マイクロホンアレイにより録音された。
【０２７０】
図２７（ａ）の上側プロットおよび下側プロットには、それぞれ２つの典型的なマイクロホン（すなわち、日の当たる側にある１３番のマイクロホンおよび日の当たらない側にある３１番のマイクロホン）で受信された信号が示されている。上側プロットに示されている信号の短時間フーリエ変換を用いたスペクトログラムが、中央プロットに示されている。
【０２７１】
サブセクションＴ．Ａ．で示されるＴＤＭＲモーダルビームフォーマが用いられている。到来（入ってくる）方向、すなわち（２０°、１８０°）にビームが向けられる場合のビームフォーマの時系列出力およびスペクトログラムが、それぞれ、図２７（ｂ）の上側プロットおよび中央プロットに示されている。図２７（ｂ）の下側プロットには、到来方向から６０°ズラされている他の方向（８０°１８０°）にビームが向けられる場合の時系列出力が示されている。
【０２７２】
サブセクションＴ．Ｂ．およびＴ．Ｃ．に記載のＴＤＭＤモーダルビームフォーマおよびＴＤＲＭＤモーダルビームフォーマを同一のマイクロホンアレイデータに適用する。上述のプロセスが繰り返えされる。これら２つの方法を用いた場合の図２７（ｂ）と同様の結果が、それぞれ図２７（ｃ）および図２７（ｄ）に示されている。
【０２７３】
図２７（ｂ）、図２７（ｃ）および図２７（ｄ）の上側プロットをについて説明する。ＴＤＭＲＤビームフォーマの出力がＴＤＭＲビームフォーマのものと類似していることが分かる。しかしながら、ＴＤＭＤビームフォーマの場合、低周波数における大きさがはるかに大きい。この理由は、低周波数における重み付けノルムが非常に大きく、期待されるアレイ応答ベクトルと実際のアレイ応答ベクトルと間のミスマッチが僅かであっても出力が大きくなるからである。換言すれば、このビームフォーマは少しのミスマッチに対しても非常に敏感であるということである。
【０２７４】
図２７（ｂ）の下側プロットを図２７（ｄ）の下側プロットと比較すると、ＴＤＭＲビームフォーマの時系列の大きさが、とくに低周波数において、ＴＤＲＭＤビームフォーマのものよりはるかに大きく、前者のビーム幅が後者のビーム幅よりも広いことを意味していることが分かる。このことは、図２２および図２４に示されているビームパターンからも見出すことができる。したがって、図２７に記載の結果から、ＴＤＲＭＤビームフォーマが指向性とロバストネスとの間の良好なトレードオフを提供していることが分かる。
【０２７５】
上記の実施例は、球面調和関数ドメインの広帯域モーダルビームフォーマを実数値時間ドメイン（ｒｅａｌ−ｖａｌｕｅｄ  ｔｉｍｅ−ｄｏｍａｉｎ）で具象化したものである。これらの実施例に記載の広帯域モーダルビームフォーマはモーダル変換ユニットと、ステアリングユニットと、パターン作成ユニットとを備えている。ステアリングユニットは、任意選択的なものであり、ルック方向に対して回転対称ではないビームパターンの作成が必要な場合には省略することができる。パターン作成ユニットは、ステアリング方向に依存せず、フィルタ和構造（ｆｉｌｔｅｒ−ａｎｄ−ｓｕｍ  ｓｔｒｕｃｔｕｒｅ）を用いて実現される。優れた球面調和関数フレームワークは、従来のエレメントスペース基づいたアプローチ（ｅｌｅｍｅｎｔ−ｓｐａｃｅ  ｂａｓｅｄ  ａｐｐｒｏａｃｈ）に比べて演算がより効率的な最適化アルゴリズムおよび実現スキームに結びつく。広帯域アレイ応答、等方性ノイズおよび空間ホワイトノイズに対するビームフォーマ出力および主ローブ空間応答変化はすべて、ＦＩＲフィルタのタップ重みの関数として表されている。ＦＩＲフィルタ設計問題は、多重制約問題として定式化されており、このことにより、得られるビームフォーマが、指向性、主ローブ空間応答変化、副ローブレベルおよびロバストネスの如き相反する複数のアレイ性能尺度間の適切なトレードオフを提供することができることが担保されることになる。
【０２７６】
上述のすべてから理解することができるように、二次錘計画法ソルバーを用いて効率的に解くことできる多重制約凸最適化問題として最適化問題を定式化することによって、球面マイクロホンアレイの最適ビームフォーマ設計の問題が取り組まれている。得られたビームフォーマが、指向性指数、ロバストネス、アレイ利得、副ローブレベル、主ローブ幅などの如き複数の性能尺度間の適切なトレードオフを提供し、異なるローブ／領域の利得制約を変えながら、複数の主ローブの形成および干渉の回避のために、適応する複数のヌルを形成することができることが実証されている。このアプローチによって柔軟な設計ツールが提供されることは明白である。というのは、このアプローチが以前に学習された遅延和型ビームフォーマ、および特別なケースとして純粋位相モードビームフォーマを包含するとともに、もっと複雑な最適化問題を許容可能な時間枠内で解決することを可能とするからである。
添付資料
下記のセクションでは、球面フーリエ変換および球面調和関数に基づいたビームフォーミングの背景が説明してあり、この明細書において用いられているいくつかの結果が導出されている。
【０２７７】
標準直交座標系（ｘ、ｙ、ｚ）および球面座標系（ｒ、θ、φ）が用いられている。
【０２７８】
ここでは、仰角および方位角（アジマス）は、＋ｚ軸および＋ｘ軸から測定されたラジアン単位で表される角度変位をｚ＝０の面に投影したものである。Ω0＝（θ0、φ0）の方向から半径ａの球体に入ってくる単位強度の平面波を想定する。この用途では、時間係数ｅｘｐ（ｉωｔ）が抑制される。ここで、ｉは
【数１４７】
2012523731.tif 000149
であり、ωは時間角周波数（ｔｅｍｐｏｒａｌ  ｒａｄｉａｎ  ｆｒｅｑｕｅｎｃｙ）である。
【０２７９】
波数ｋについての球体表面における観測点での総合音圧波は、球面調和関数を用いて次のように書くことができる：
【０２８０】
【数１４８】
2012523731.tif 000150
この式で、ｋ＝‖ｋ‖＝ｗ／ｃであり、ｃは音速であり、
【数１４９】
2012523731.tif 000151
は位数がｎであり、次数がｍである球面調和関数であり、上付き文字＊は複素共役を表しており。ｂn（ｂａ）は、球体構成、たとえば、剛球（ｒｉｇｉｄ  ｓｐｈｅｒｅ）、開球（ｏｐｅｎ  ｓｐｈｅｒｅ）などに応じて次の式で求められる：
【０２８１】
【数１５０】
2012523731.tif 000152
この式で、ｊnおよびｈnは、それぞれｎ位数の球状ベッセル関数およびハンケル関数であり、
【数１５１】
2012523731.tif 000153
および
【数１５２】
2012523731.tif 000154
は、それぞれｊnおよびｈnの独立変数（ａｒｇｕｍｅｎｔ）に対する導関数（ｄｅｒｉｖａｔｉｖｅ）である。
【０２８２】
球面調和関数は、波動方程式、または球座標のヘルムホルツの式の解である。これらは、次の式で求められる：
【０２８３】
【数１５３】
2012523731.tif 000155
この式で、
【数１５４】
2012523731.tif 000156
はルジャンドル陪関数を表している。球面調和関数は、正規直交であり、次の式を満たす：
【０２８４】
【数１５５】
2012523731.tif 000157
この式で、δn-n'およびδm-m'は、クロネッカー（Ｋｒｏｎｅｃｋｅｒ）デルタ関数であり、積分
【数１５６】
2012523731.tif 000158
は単位球面Ｓ2の全表面をカバーしている。
【０２８５】
単位球面上で２乗積分可能な関数（ｓｑｕａｒｅｄ  ｉｎｔｅｇｒａｂｌｅ  ｆｕｎｃｔｉｏｎ）ｐの球面調和関数分解または球面フーリエ変換は、ｐnmで表され、その逆変換と共に次のそれぞれ対応する式で求められる：
【０２８６】
【数１５７】
2012523731.tif 000159
【０２８７】
【数１５８】
2012523731.tif 000160
球面フーリエ変換（５）を（１）で表されるような平面波に適用すると、球面調和関数ドメイン表現ｐ（ｋａ、Ω0、Ω）が求められる：
【０２８８】
【数１５９】
2012523731.tif 000161
ここで、球面アレイの性能を解析するために、球面に入ってくる、Ω0方向からの目的信号（ＳＯＩ）平面波と、Ω1、．．．、Ωd、．．．、ΩD方向からのＤ干渉平面波とを想定する。相関関係のないノイズを加えると、球面上の音圧を次のように書くことができる：
【０２８９】
【数１６０】
2012523731.tif 000162
この式で、
【数１６１】
2012523731.tif 000163
はＤ＋１個のソース信号スペクトルであり、Ｎ（ω）は付加ノイズスペクトルであり、βはＳＯＩが存在するか否かを示すバイナリパラメータ（ｂｉｎａｒｙ  ｐａｒａｍｅｔｅｒ）である。
【０２９０】
ｘ（ｋａ、Ωs）の球面フーリエ変換は、次の式で求められる。
【０２９１】
【数１６２】
2012523731.tif 000164
この式で、
【数１６３】
2012523731.tif 000165
はノイズの球面フーリエ変換を表している。
【０２９２】
アレイ処理は、アレイ入力信号とアレイ重み付け関数との積を全球面にわたって積分することにより、または球面調和関数ドメインにおいて同様の重み付けおよび総和を行うことにより、スペースドメインまたは球面調和関数ドメインのどちらにおいてでも処理することができる。アパーチャ重み付け関数をｗで表すことにより、アレイ出力は、アレイ入力信号と複素共役重み関数（ｃｏｍｐｌｅｘ  ｃｏｎｊｕｇａｔｅｄ  ｗｅｉｇｈｔｉｎｇ  ｆｕｎｃｔｉｏｎ）ｗ*との積の球面全体にわたる積分として求められる。
【０２９３】
【数１６４】
2012523731.tif 000166
この式で、ｗnmはｗの球面フーリエ変換係数である。（１０）内の総和項は、球面調和関数ドメインにおける重み付けと考えることができ、位相モード処理とも呼ばれることに留意されたい。
【０２９４】
実際問題として、音圧は、マイクロホン位置Ωsで空間的にサンプリングされる。ここで、ｓ＝１、．．．、Ｍであり、Ｍはマイクロホンの数である。マイクロホン位置は次の離散的正規直交性条件（ｄｉｓｃｒｅｔｅ  ｏｒｔｈｎｏｒｍａｌｉｔｙ  ｃｏｎｄｉｔｉｏｎ）を満たす必要がある：
【０２９５】
【数１６５】
2012523731.tif 000167
この式で、αsはサンプリングスキームによって異なる。サンプリングが均一（ｕｎｉｆｏｒｍ）な場合、
【数１６６】
2012523731.tif 000168
であるためには、αs≡４π／Ｍである。 いうまでもなく、球面上にマイクロホンを配設する他の空間的なサンプリングスキームにおいても同様に有効である。
【０２９６】
有限数の数のマイクロホンが球面をサンプリングする場合、位数Ｎの球面調和関数は、空間エイリアシング を回避するために、Ｍ≧（Ｎ＋１）2を満たす必要があることに留意されたい。換言すれば、位数Ｎが与えられた場合、マイクロホンの数Ｍは少なくとも（Ｎ＋１）2でなければならない。
【０２９７】
ｘ（ｋａ、Ωｓ）の離散球面フーリエ変換（球面フーリエ係数）およびその逆変換はそれぞれ対応する次の式で求められる：
【０２９８】
【数１６７】
2012523731.tif 000169
【０２９９】
【数１６８】
2012523731.tif 000170
解析を単純化するために、この明細書では、マイクロホンによる空間サンプリングは完全であり、空間エイリアシングは無視できるので、αs≡４π／Ｍであると仮定する。
【０３００】
これに対応するアレイ出力は次の式により計算することができる：
【０３０１】
【数１６９】
2012523731.tif 000171
この式で、ｗ*（ｋ、Ωs）はアレイ重みである。
【数１７０】
2012523731.tif 000172
はそれらの球面フーリエ係数である。 理想的に均一にサンプリングが行われる場合、（１４）内のアレイ出力の振幅は
【数１７１】
2012523731.tif 000173
である古典的なアレイ処理よりも４π／Ｍ倍だけ大きくなることに留意されたい。
【０３０２】
球面フーリエ変換に対する重み付けのパルゼバル関係を用いると、次の式が得られる：      
【０３０３】
【数１７２】
2012523731.tif 000174
この式は因数αsを示す。
【図１】
2012523731.tif 000175
【図２】
2012523731.tif 000176
【図３】
2012523731.tif 000177
【図４Ａ−４Ｃ】
2012523731.tif 000178
【図５】
2012523731.tif 000179
【図６】
2012523731.tif 000180
【図７Ａ−７Ｃ】
2012523731.tif 000181
【図８Ａ】
2012523731.tif 000182
【図８Ｂ】
2012523731.tif 000183
【図９Ａ】
2012523731.tif 000184
【図９Ｂ】
2012523731.tif 000185
【図１０Ａ−１０Ｂ】
2012523731.tif 000186
【図１１Ａ−１１Ｂ】
2012523731.tif 000187
【図１２Ａ】
2012523731.tif 000188
【図１２Ｂ】
2012523731.tif 000189
【図１３】
2012523731.tif 000190
【図１４】
2012523731.tif 000191
【図１５】
2012523731.tif 000192
【図１６】
2012523731.tif 000193
【図１７】
2012523731.tif 000194
【図１８Ａ】
2012523731.tif 000195
【図１８Ｂ】
2012523731.tif 000196
【図１９】
2012523731.tif 000197
【図２０】
2012523731.tif 000198
【図２１】
2012523731.tif 000199
【図２２Ａ】
2012523731.tif 000200
【図２２Ｂ】
2012523731.tif 000201
【図２３Ａ】
2012523731.tif 000202
【図２３Ｂ】
2012523731.tif 000203
【図２４Ａ】
2012523731.tif 000204
【図２４Ｂ】
2012523731.tif 000205
【図２４Ｄ】
2012523731.tif 000206
【図２５Ｂ】
2012523731.tif 000207
【図２６Ｂ】
2012523731.tif 000208
【図２２Ｃ】
2012523731.tif 000209
【図２２Ｄ】
2012523731.tif 000210
【図２３Ｃ】
2012523731.tif 000211
【図２３Ｄ】
2012523731.tif 000212
【図２４Ｃ】
2012523731.tif 000213
【図２５Ａ】
2012523731.tif 000214
【図２５Ｃ】
2012523731.tif 000215
【図２６Ａ】
2012523731.tif 000216
【図２６Ｃ】
2012523731.tif 000217
【図２７Ａ】
2012523731.tif 000218
【図２７Ｂ】
2012523731.tif 000219
【図２７Ｃ】
2012523731.tif 000220
【図２７Ｄ】
2012523731.tif 000221
【国際調査報告】
2012523731.tif 000222
2012523731.tif 000223
2012523731.tif 000224
2012523731.tif 000225