JP3859594B2 - 適応フィルタリングのための装置 - Google Patents
適応フィルタリングのための装置 Download PDFInfo
- Publication number
- JP3859594B2 JP3859594B2 JP2002524089A JP2002524089A JP3859594B2 JP 3859594 B2 JP3859594 B2 JP 3859594B2 JP 2002524089 A JP2002524089 A JP 2002524089A JP 2002524089 A JP2002524089 A JP 2002524089A JP 3859594 B2 JP3859594 B2 JP 3859594B2
- Authority
- JP
- Japan
- Prior art keywords
- vector
- stage
- group
- inner product
- correlation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0025—Particular filtering methods
- H03H21/0029—Particular filtering methods based on statistics
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Radio Transmission System (AREA)
- Complex Calculations (AREA)
- Filters That Use Time-Delay Elements (AREA)
Description
(関連する出願に関するクロスリファレンス)
本願発明は、いずれもが本文においてその全文に亘って参考として織り込まれている、出願人によって2000年8月29日付けで出願された「適応アルゴリズムを実行するための有効アーキテクチャー」という表題にて現在出願中の特許出願第60/228,408号および2000年10月27日付けで出願された「非定常データを処理するための方法およびシステム」という表題の特許出願第60/243,480号に関連し、これらに基づいて優先権を主張するものである。
【0002】
(連邦政府によって支援される研究に関する記述)
本願発明は、海軍調査課によって裁定された契約N00014−00−C−0068として政府のサポートを受けてなされたものである。政府は、本願発明に対する一定の権利を有するものである。米政府は、本願発明に対して支払済ライセンスを有し、海軍調査課によって裁定された契約N00014−00−C−0068の項目によって規定される穏当な条件において特許権者が第3者に対してライセンスを付与させる権利を限られた環境下において有するものである。
【0003】
(発明の分野)
本発明は、概してデータ群の離散表示を処理するための装置に関連するものである。より詳しくは、本発明は、データ群を統計的相関によって関連付けられた直交成分へと分解するための装置(システム)に関するものである。この装置は、観測離散データから対象となるデータを回復するために用いられ、この場合、これらの観測データは、例えばベクトルデータからスカラー信号を推定する場合において対象となるデータとその他のデータとの両方を含むものであり得る。
【0004】
(背景情報)
1.計算コストに関する注釈
演算の計算コストは、通常「浮動小数点演算」あるいは「FLOP」の回数を数えることによって評価される。この評価されたコストは、「FLOPカウント」と称される。例えば、マトリックスMと縦ベクトルvがあるとして、MおよびvはいずれもがN個の行を有し、MはK個の列を有するものとする。この積Mvを演算するには、スカラー・スカラーの積の合計演算回数はNKとなり、スカラー・スカラーの和の合計演算回数はN(K−1)となる。したがって、FLOPカウントの合計はN(2K−1)あるいはおよそ2NKとなる。実際においてはその他の操作(例えばメモリにおけるデータの移動など)が必要となるため、このFLOPカウントは、実際の計算コストの概算に過ぎないが、FLOPカウントは一般的に認められているコストのインジケータである。FLOPカウントを定義する一つの参考文献としてはG.H.ゴルブおよびC.F.ヴァン・ローン著「マトリックス計算」第3版、ジョン・ホプキンス・ユニバーシティ出版、1996(特に18乃至19ページ)が挙げられ、各例が254、263および270ページにおいて一覧にされている。
【0005】
通常FLOPカウントは、処理されるべきデータの数とともにカウントがどのように増加するかを見積もる「成長の度合い」を表す簡単な方式によって評価される。上述の例において、FLOPカウントの増加の程度はO(NK)であって、ここにおける記号「O( )」は「オーダー」を表すものである。「O( )」とはデータ入力の増加に伴う成長あるいはスケーリングの法則を表すものであるため、この表示O(NK)では係数の2を省略している。すなわち、係数2はNおよびKの直線的成長と比較した場合無視できる項目である。「成長のオーダー」を定義する一つの参考文献としてはハロルド・アベルソンおよびジェラルド・サスマン著「コンピュータプログラムの構造および解釈」、MIT出版およびマグローヒル、39乃至40ページが挙げられる。
【0006】
このような状況において、我々の目標は、O(N 3 )のオーダーで増加する計算方法を避け、O(N 2 )あるいはO(N)のオーダーで増加する方法を見出すことにある。
【0007】
2.例示的な問題ドメイン
本願において提示されている例示的な問題ドメインは、離散データサンプル群に対する適応信号処理である。本願においては、このドメインに関連して本願発明による好ましい実施例が開示されている。対象となる一般的なブロック図が、図1において示されている。この図において、スカラー出力ε0(k)群を得るためにフィルタがスカラー入力d0(k)群とベクトル入力x0(k)群とに対して作用する。スカラー入力とベクトル入力とは、ともに入力データを特徴づけるものである。整数指数kは、このデータが離散サンプルからなることを示す。フィルタウェイトw0は、例えばコスト関数を最小化するため、あるいは例えば信号対雑音比などのその他の方法を最大化するためなど、ある種の性能測度を最適化するために選択される。一例として、二次コスト関数を最小化するためにフィルタウェイトが選択される場合、このフィルタは「ウィーナ」あるいは「最小自乗」フィルタとして知られる(サイモン・ヘイキン著「適応フィルタ理論」第3版、特に194および483ページ、プレンティスホールを参照のこと)。ここにおいて、ウィーナフィルタという用語を使用する。
【0008】
より具体的な例として、図2においてウィーナフィルタを用いて離散データ入力サンプル群x(k)を処理してステアリングベクトルs(あるいはレプリカベクトルまたはフォーカシングベクトルとも称される)に最もマッチするデータの部分を得るケースが示されている。例えば、単周波適応ビーム形成において、sとは信号がアンテナアレイに沿ってトレース可能な空間パターンのモデルであり、sは典型的には複素量である。無線通信においては、sは符号分割多重アクセス(CDMA)システムにおける個人のユーザの拡張コードであり得る。時空適応処理(STAP)へのアプリケーションにおいては、sは時間および空間における(周波数および角度へと変換される)パターンであり得る。これら全てのアプリケーションを共通の表記法にて表記するために、ステアリングベクトルを正規化することにより、単位ノルムを有する無次元量とした。
【数1】
【0009】
これらいずれのアプリケーションにおいても、ブロッキングマトリックスBにより、x(k)のうちsに対して直交する成分が求められる。したがって、d0(k)=sHx(k)がsの方向におけるデータの部分を表し、x0(k)=Bx(k)がsに対して直交するデータの残りの部分を表す。
【0010】
ここにおいて図2におけるフィルタが図3において示されるデータベクトルのブロックにどのように適用されるかを示す。x(k)がN個のエントリーを有する観測データの縦ベクトルとする。これら観測データは複素量であってもよい。指数kとは、x(k)がデータベクトルのブロックにおける1つのデータベクトルであることを示しており、1≦k≦Kである。増加するkの値は通常増加する時間に対応する。データベクトルのブロックが以下のN×Kマトリックスであるとする。
【数2】
【0011】
この例において、ステアリングベクトルまたはレプリカベクトルsに最もマッチする、観測データの部分を抽出するためにXをフィルタリングすることが目的である。
【0012】
3.観測データから対象となるデータを検索するためのアプローチ
非適応型アプローチ
対象となるデータを検索するための簡単なアプローチとしては、入力データベクトルをステアリングベクトルsに積算することであり、よってフィルタリングの結果は積sHXとなる。このアプローチの一般化は、ステアリングベクトルsのある関数であるウェイトベクトルwを定義することにあり、この場合のフィルタリングの結果は積wHXとなる。例えば、wはサイドローブを低減するために「テーパリング」を適応したsのコピーであり得る。この種のフィルタリングは、ウェイトベクトルwがデータXに依存しないため「非適応型」と称される。ビーム形成において用いられる場合、このアプローチは、従来のビーム形成と称される。(例えばCDMAなど)別のアプリケーションにおいて用いられる場合、これはマッチドフィルタリングと称される。
【0013】
データ適応アプローチ全般について
対象となるデータを検索するためのより高度なアプローチとしては、ステアリングベクトルsと処理されるデータXとに依存するウェイトベクトルwを計算することである。このフィルタリングの結果は積wHXとなり、この種のフィルタリングは、ウェイトベクトルwがデータXに依存するため「データ適応型」あるいは「適応型」と称される。多くの適応型の方法が存在する。一つの概論本としてはサイモン・ヘイキン著「適応フィルタ理論」第3版、プレンティスホール、1996年があり、また別の概論本としてはバーナード・ウィドロウおよびサミュエル・スターンズ著「適応信号処理」、プレンティスホール、1985年がある。このアプローチにおいて、フィルタリングはそのデータベクトル群にカスタマイズされる。
【0014】
「最適な」ウェイトベクトルwを定義するために、適応型アプローチにおいては「性能関数」を定義し、次にこの性能関数を最適化するウェイトベクトルを見つける。おそらく最も一般的である性能関数とは最小化されるべき二次式あるいはエネルギー状のコスト関数であるが、その他の性能関数であってもよく、ここにおいて詳述される発明はある特定の性能関数に限定されるものではない。
【0015】
最小分散無ひずみレスポンス(MVDR)
離散観測データ群Xから対象となるデータを検索するための一つの簡単でポピュラーである適応型アプローチとしては、最小分散無ひずみレスポンス(MVDR)の解を求めるものがある。このアプローチにおいて、ウェイトベクトルwがレプリカベクトルsにマッチするデータをパスするという制約の条件下において二次のコスト関数、すなわちエネルギー状のコスト関数を最小化するためのウェイトベクトルwに対してデータベクトルx(k)が積算される。ここにおける本発明の好ましい実施例においては、以下のようなK個の有限である群の「ブロック平均」である平均(全体平均ではない)を用いる。
【数3】
【0016】
最小化は以下の制約wHS=1に基づき以下のように表される。
【数4】
【0017】
4.MVDRの解を得るための伝統的な適応型フィルタアプローチ
MVDRの解を得るための伝統的な適応型フィルタアプローチにおいては、解を以下のデータ相関あるいは共分散マトリックスによって表す。
【数5】
【0018】
Rを用いて最小値(4)は以下のように書き換えられ得る。
【数6】
【0019】
この伝統的なアプローチにおいては、Rが最大階数である(すなわちRは逆行列R −1 を有する)と仮定する。(4)または(6)の伝統的な解は以下となる。
【数7】
【0020】
Rを演算するために、1ブロック分の(5)を評価するにはN自乗×K個のオーダー、つまりO(N2K)の積算処理が必要となる。データベクトルT全ての数についてマトリックスを演算するためにはO(N2T)の処理が必要となる。R-1を演算するために、1ブロック分のマトリックス反転を行うためにはO(N3)の処理が必要となり、データベクトルT全ての数についてマトリックスを反転するためにはO(N3T/K)の処理が必要となる。各ブロックにおけるR-1に対してSを積算するためには、O(N2)の処理が必要となり、データベクトルT全ての数に対してはO(N2T/K)の処理が必要となる。したがって、伝統的なアプローチにおいてはO(N2)あるいはO(N3)のオーダーのコストが必要となり、多くの実際的なアプリケーションにおいてNが大きい場合にはシステムに対して面倒な処理負荷となる。
【0021】
5.MVDRのためのマルチステージウィーナフィルタ(MWF)による解決
マルチステージウィーナフィルタは、信号の検出、推定および分類のために用いられ得るデータ適応型フィルタである。例えば、ステアリングベクトルsとともに用いられる場合、MWFはそのデータがどの程度レプリカベクトルsに類似するかを推定するために一つ以上のデータベクトルx(k)のブロックストリームを処理することができる。より高度なアプリケーションにおいては、MWFは(単一のステアリングベクトルsがステアリングベクトルのマトリックスによって置き換えられ得る場合の)複数の制約を満足させることが可能である。
【0022】
J.S.ゴールドスタイン、I.S.リードおよびL.L.シャーフによって、本文においてその全体に亘って参考として織り込まれている「直交射影に基づくウィーナフィルタのマルチステージ的表示」(「情報理論に対するIEEE処理」第44版、7巻、2943−2959ページ、1998年11月[ゴールドスタイン]、においてマルチステージ的分解を用いたウィーナフィルタ構造が開示されている。[ゴールドスタイン]においては概括サイドロープキャンセラ(GSC)MDVR−制約付きウィーナフィルタとこのようなフィルタをMWFとして実施するための方法が開示されている。このような2つの段階を有するフィルタの例が図2において示されている。
【0023】
[ゴールドスタイン]におけるMWFは、例えばRが逆行列を有しない場合などにおいてMWFが機能するなど、ここにおいて数式7として示されている伝統的なMVDRの解に対して有利であるMVDRに基づく解を提供するものである。データが定常であってMFWが最大階数まで演算されるという特別なケースにおいて、データ共分散マトリックスに基づいてMVDRを求める伝統的な解と同等の性能を付与するものである。
【0024】
より具体的には、実世界における
・ 無線通信−使用者が受信機をON/OFF切り替えする場合
・ 受動的ビーム形成−波源が受信機に対して移動する場合
・ 能動的レーダーおよびソナー−クラッタ統計が地底あるいは海底に沿って幅広く変動する場合
などのアプリケーションにおいては非定常データが用いられる。このような非定常データを鑑みた場合、MWFはデータ共分散マトリックス反転に基づく伝統的な解よりもよりよい性能を示す。MWFが低い階数で演算される場合であってサンプルサポートが低い場合にはMVDRによる解より効率がよい場合がある。MWFは共分散マトリックスよりもより直接的且つ効率的に干渉を表す統計に基づいている。MWFフィルタは、無線通信のシミュレーションにおいて調査されて、実際のレーダーデータを用いて実験された。各ケースにおいてその技術的性能は従来の方法のそれに匹敵するものあるいはそれを超越するものであり得る。
【0025】
しかしながら図2において示されているように、[ゴールドスタイン]におけるMWFは、sに対して直交するデータを得るためにブロッキングマトリックスBとデータx(k)とのマトリックス積を演算することを求めている。忠実に演算された場合、マトリックス積Bx(k)には、1≦k≦Kのブロックに対してO(N2K)のコストが必要となる。その後の適応段階においては、相関方向ベクトルhiに対して直交であるデータを得るためにより多くのブロッキングマトリックスBiが必要となる。各段については、O(N2K)のコストが必要となる。
【0026】
ブロッキングマトリックスをあからさまに用いずに直交分解を実施することが求められている。
【0027】
(発明の概要)
スカラー群とベクトル群とからなるデータを分析し、統計的相関によって関連付けられた成分に分析するための発明である。好ましい実施例において、本発明は、スカラー群とベクトル群とを入力として受け取り;スカラーおよびベクトル入力に関連付けられた相関方向ベクトルを演算し;相関方向ベクトルを用いて入力ベクトルの内積を演算し;内積を相関方向ベクトルに積算してスケーリングされた相関方向ベクトル群を形成し;スケーリングされた相関方向ベクトルを入力ベクトルから減算して相関方向ベクトルに対して直交であるまたは実質的に直交である入力ベクトルの射影を得るステップあるいは手段を有するものである。その出力は、スカラー内積群および相関ベクトルに対して直交であるあるいは実質的に直交であるベクトル群である。これらステップあるいは手段は、カスケード状に繰り返されてデータのマルチステージ的分析を形成することが可能である。本発明は、また適応型分析段階に先立ってステアリングベクトルとともに用いることも可能である。好ましい実施例において、本発明は従来の方法よりもより少ない演算的リソースを必要とするものである。
【0028】
ここにおいて説明される好ましい実施例は、信号処理問題に対する最適な解を演算するための基礎を得るためにデータを分析するための方法に関するものである。伝統的には、データ平均された共分散マトリックスを演算してその後マトリックスを分析し、各固有値および固有ベクトルを求めるアプローチが一般的であった。より最近のアプローチは、データを統計的相関に基づいて識別された直交成分へとマルチステージ的に分解することである。ここにおける発明は、統計的相関によって識別されたマルチステージ的分解である後者のアプローチを適応するものである。
【0029】
発明の詳細な説明
通例
図4において、本発明の好ましい実施例が概略的に示されている。上側のチェーンは、(左から右へと流れる)データを分析するものである。下側のチェーンは、フィルタウェイトとフィルタ出力(右から左へと流れる)とを合成するものである。
【0030】
ステアリングベクトルsを用いる例においては、処理は非適応型演算から始まる。最初に、データベクトルがステアリングベクトルに対して射影され、以下のような所望された信号の初期推定を得る。
【数8】
【0031】
例えば、d0(k)は、従来のビーム形成器あるいはマッチドフィルタの出力であり得る。一般的に、d0(k)は、sのサイドローブによってもたらされる干渉を含有するものである。この干渉を識別して減算するための統計的分析を得るために、以下のように記載して残りのデータを絶縁することが可能である:
【数9】
ここにおけるBはsに対して直交であるデータの射影、すなわちsの零空間に対する射影を得る「ブロッキングマトリックス」である。
【0032】
(9)における射影処理は、一意的に定義される。しかしながら、射影は少なくとも2通りに説明することができる。[ゴールドスタイン]において説明されるものを含む以前のアプローチにおいては、ブロッキングマトリックスは長方形であり、結果はN−1次元のベクトル(すなわちその長さがN−1であるベクトル)である。もう一方の可能性としては以下の正方ブロッキングマトリックスを考える。
【数10】
【0033】
この正方マトリックスによる射影は、上記と同一のものであるが、この射影は元のN次元の(すなわち全てのベクトルが長さNを有する)座標システムの固定された視点においては減算としてみなされ得る。
【0034】
これら2つのブロッキングマトリックス間の選択は、本発明において実現される数値的効率に関係するものである。本発明の好ましい実施例において(10)における正方ブロッキングマトリックスを選択することによって(9)を以下のように書き換えることが可能である。
【数11】
【0035】
(11)において示される減算を用いて本発明の好ましい実施例における射影を得ることによってブロックのための計算コストはO(NK)である。これにより、(9)で表される、O(N2K)のコストを必要とするマトリックス積算と比べて明らかな節約が可能となる。
【0036】
d0(k)およびx0(k)が入力として与えられた場合、処理は次に適応段階の再帰をもって続行される。指数iは適応段階数を識別するものである(ここにおいてi=1は、第1の適応段階を示す、など)。段階iをサポートするために、この段階に対する2つの入力間のブロック平均された相関r i-1 を演算する。
【数12】
【0037】
段階iにおいては、相関の大きさδiと相関方向ベクトルhiとにおいてこの相関を用いる。
【数13】
【数14】
【0038】
本発明の好ましい実施例においては、この相関方向ベクトルをブロックの均一加重をすることによるブロック平均化を用いて演算する一方、特定のアプリケーション(すなわち時間的拘束、精度に関する要件など)および処理リソースの利用可能性などに基づいて、当業者にとって明らかである他の平均化技術を用いることも可能である。
【0039】
(11)において得られた数値的効率を念頭において、この方向に沿う、またこの方向に対して直交であるデータの射影を得るために、本発明の好ましい実施例においてこの射影は減算としてみなされる。
【数15】
【数16】
【0040】
この再帰的分析は、様々な態様で終了しうる。データXのブロックが最大階数である場合、これら射影は全てのデータが直交単位ベクトル群[s,h1,h2,...,hN−1]に対して射影されるまで継続され得る。この分析においてはN−1段階が必要となる(なぜならばsが既に1つの方向に該当するからである)。あるいは、データXのブロックが最大階数でない場合(例えば意図的にK<Nを選択した場合)データの階数はこの分析において適切に「アンダーフロー」する。この場合、xi(k)は最も高い各段階においてゼロのみを含有するものとなる。あるいは、フィルタリングを意図的にここではS(但し1≦S≦N−1)と称するより少ない数の段階へとトランケートすることも可能である。トランケーションは、図2においてグラフの形で示されており、「ターミネータ」がεS(k)=dS(k)と設定する。
【0041】
分析が終了した後、初期化εS(k)=dS(k)によって図4における下側のチェーンにおいて右から左への合成が始まる。この合成は再帰的である。各適応段階i(i=Sからi=lまで)について以下の演算を行う。
【数17】
【数18】
【数19】
【0042】
すなわち、直後の段階からの出力ε i (k) の自乗平均ξ i が求められ、相関の大きさδ i を前記自乗平均(ξ i )で除してウエイトw i が求められる。さらに、直後の段階からの出力ε i (k) およびウエイトw i を乗じて、スケーリングされた出力ε i ・w i が求められ、これと入力スカラー群d i-1 (k) との間の差d i-1 (k) −ε i ・w i を求めて、現在の段階のスカラー出力ε i-1 (k) (ε0(k))を得る。
【0043】
陽関数表示にてウェイトベクトルを使用するその他のフィルタと比較しつつ、本発明に関する理解をさらに深めるため、以下の説明を行なう。合成ε0(k)は、以下のように書き換えられることに注意されたい。
【数20】
【0044】
さらに、(ここで詳述される好ましい実施例と同様)直交ベクトル[s,h1,h2,h3,・・・]のうちいずれもがN個のエントリーを有する場合、x(k)に対する従属は直接以下のように表され得る。
【数21】
【0045】
これによって以下のマルチステージ式ウィーナフィルタのための同等のウェイトベクトルが識別される。
【数22】
【0046】
アプリケーションによっては、図2におけるフィルタ出力ε0(k)に対して直接アクセス可能であるため、このウェイトベクトルを数値的に実行する必要はない。また別のアプリケーションにおいて、ウェイトベクトルは利用価値のある結果であり得る。
【0047】
本発明の好ましい実施例の例として、図5において水中音響問題のシミュレーションが示されている。N=100とし、水平ラインアレイがノイズの背景において3つの空間的離散音源を分解している。離散的音源およびノイズは全て擬似乱数複素ガウス変数として計算されて単一の高速フーリエ変換(FFT)周波数ビンの内容を表す。データベクトルは一秒毎に収集される。図面における単位は1μPa2/HzにつきdBであって(μPa2にて変数を与えるために周波数スペクトラムが積分するような)受動的ソナー問題における周波数ビンの単位を表す。要素レベルにおいて、周辺音響ノイズレベルは(1μPa2/Hzにつき)65dBであり、この例を簡素化するためにこの周辺ノイズは等方性のものである。処理される周波数ビンは、空間における偽信号が生じないようにアレイが処理できる最大周波数の1/3の周波数に集中しており、等方性周辺ノイズに対して要素間の相関を付与する。相関されない要素ノイズも(1μPa2/Hzにつき)25dBの効果的なレベルにおいてモデル化される。ここでは3つの空間的離散ソースはエンドファイヤからの到来角がそれぞれ0°、30°、60°である局地的平面波としてモデル化されている。要素レベルにおけるこれらのRMS振幅はそれぞれ(1μPa2/Hzにつき)85、65、75dBである。ビーム形成は、そのインクレメントがΔφ=0.25°であるステアリング角度0≦φ≦180°の範囲について演算される(ここにおいてこの極めて細かい解が用いられるのは、幅狭であるピークのピークレベルがTIFFグラフィック出力において正しく表示されるようにするためである)。これら結果は、K=200データベクトルである一つのブロックについて示されている(200秒をカバーする)。
【0048】
図5(a)において振幅および角度分解能のリファレンスとして従来のビーム形成(CBF)の結果が示されている。CBFのため、ここにおける固定のウェイトベクトルはテーパリングなしのステアリングベクトルである(wCBF=s)。示された実際の出力を演算するには複素データベクトルx(k)がウェイトベクトルに積算されてその単位がμPa2/Hzであるパワー状の量を得る。
【数23】
【0049】
K=200のデータベクトルのコストは、各ステアリング角度につきO(NK)である。1μPa2/HzにつきdBの単位を示すため図面においては10log10[PCBF/(1μPa2/Hz)]がプロットされている。
【0050】
図5(b)において(7)におけるマトリックス反転によって計算された伝統的な適応型の解が示されている。マトリックスRは、O(N2K)のコストにてK=200のデータベクトルの全ブロックについて計算され、反転はO(N3)のコストにて計算される。各ステアリング角度について適応ウェイトベクトルwTraditionalが(7)にて全ブロックについて計算され、その後データベクトルが以下のパワー状の量を計算するためにウェイトベクトルに対して積算される。
【数24】
【0051】
図5(c)において最大階数(S=N−1段階)にて計算されたMWFの解が示されており、ここにおいて各ステアリング角度に対する計算コストはO(N2K)である。プロットから明らかであるように、最大階数のMWFは伝統的な解を回復するものである。パワー状出力を(22)においてウェイトベクトルを用いて計算することが可能であるものの、最も効率的なかたちで以下のようにスカラー出力ε0(k)から計算するものである。
【数25】
【0052】
図5(d)において低減された階数(S=3段階)にて計算されたMWFの解が示されており、ここにおいて各ステアリング角度に対する計算コストはO(NKS)である。これは顕著な結果である。適応型ビーム形成の超分解を従来のビーム形成と比較できる計算コストにて得ることが可能である。
【0053】
図6において本発明における仕様の分析チェーンが概略的に示されている。この例示的な実施例は、1つの非適応型段階100と1つの適応型段階200を含むものである。この例において「前段階」という記述子を使用して考えられる最初の適応型段階200の後に次の適応型段階を加えることが可能である。200において用いられる表記は、第1の適応段階を表す。指数iは以下の説明において段階を示すために用いられるものであって、i=0は、i=1である適応段階の「前段階」である非適応段階に相当する。
【0054】
非適応型分析段階100は、データベクトル群x(k)とステアリングベクトルsとを入力として受け取り、各データベクトルx(k)とステアリングベクトルsとの内積群d0(k)、すなわちd0(k)=sHx(k)を形成する内積ロジック装置120を有する。ロジック装置120と連携しているベクトルスケーリングロジック装置130は、射影群d0(k)とステアリングベクトルsを受けるものであり、全てsに対して平行であるがd0(k)によってスケールされている、すなわちsd0(k)であるベクトル群を形成する。ロジック装置130と連携しており、データベクトル群x(k)とスケーリングされたベクトル群sd0(k)を受けるためのベクトル減算ロジック装置140は、sの零空間に対するx(k)の射影に対応するベクトル差群x0(k)=x(k)−sd0(k)を形成するために用いられる。
【0055】
各適応分析段階200は、そのベクトル群xi−1(k)とスカラー群di−1(k)とからなる入力をその直前の段階から受け取る。相関方向ロジック装置210がこれら入力間の相関を計算して以下の単位ノルムを得るようにその相関方向ベクトルhiを正規化する。
【数26】
【0056】
本発明の好ましい実施例においては、この相関方向ベクトルをブロックの均一加重をすることによるブロック平均化を用いて演算する一方、特定のアプリケーション(すなわち時間的拘束、精度に関する要件など)および処理リソースの利用可能性などに基づいて、その他の当業者にとって明らかである他の平均化技術を用いることも可能である。
【0057】
その操作がロジック装置120に類似する内積ロジック装置220は、その直前の段階と現段階のロジック装置210と連携している。ロジック装置220は、ベクトル入力群xi−1(k)と現在の適応段階の相関方向ベクトルhiを受け取り、各入力ベクトルxi−1(k)の相関方向ベクトルに対する内積群di(k)、すなわちdi(k)=hi Hxi−1(k)を形成するものである。現段階におけるロジック装置210および220と連携しているベクトルスケーリングロジック装置230は、現段階の射影群di(k)と現段階の相関方向ベクトルhiを受けるものであり、全てhiに対して平行であるがdi(k)によってスケールされている、すなわちhidi(k)であるベクトル群を形成する。
【0058】
その操作がロジック装置140に類似するベクトル減算ロジック装置240は、入力ベクトル入力群xi−1(k)とスケーリングされたベクトル群hidi(k)を受けるものであり、hiに対して直交であるxi−1(k)の直交射影に対応するベクトル差群xi(k)=xi−1(k)−hidi(k)を形成するのに使用される。
【0059】
本発明の好ましい実施例において、MWFの計算コストは、データベクトルTの総数をより小さなブロックに分割するKの選択に左右されない。これに対して、伝統的なマトリックス反転において、より小さなKを選択することによって与えられたTについて反転すべきマトリックスが増えるためコストが上昇することになる。
【0060】
好ましい実施例において、入力データのチャネルあるいはデータサンプルが悪い場合には本発明の性能は卒直に劣化する。図3における手段を未だ用いながらも「悪い」チャネルあるいはサンプルにゼロを埋め込み、ステアリングベクトルsを再プログラムしてゼロを含むようにして単位ノルム(sHs=1)を有するようにすることが可能である。本発明においてゼロは影響を与えることなく流れていく。これによってデータXのブロックの階数が不足する場合、ゼロによる除法を避けるためにwi=δi/ξiにおける分母に対して小さなゼロでない数値が付与されるとフィルタはこの階数の「アンダーフロー」を適切に扱う。この変形によって重大な計算上のコストを加えるものではない。
【0061】
こうした階数における適切な「アンダーフロー」により、データベクトルKの数を意図的にセンサNの数よりも小さく選択したことによる階数不足のブロックを処理することが容易となる。センサNの数は、元の座標システムにおける次元の数に対応するものである。このようなケース(K<N)の場合、(4)または(6)における最小化は不定である(それらを決定するための式より多くの未知数がwに存在する)。これによってwには、(4)または(6)におけるコスト関数に影響を与えない余分の自由度が与えられる。このような場合、MWFは(4)または(6)において述べられたMVDRまたは「最小自乗」問題を解決する最小ノルムウェイトベクトルmin(‖w‖)を見出す。
【0062】
適応型フィルタの効率的な実行によって性能(例えば無線通信におけるビットエラー率あるいは適応型ビーム形成器におけるゲインなど)とコスト(金銭あるいは例えばウェイト、空間、パワー、メモリ、FLOPなどのその他のリソース)との間のトレードオフを可能とする。この方法を用いてコストあるいは性能のうちいずれを向上させるのかはこの技術的優位をどう費やすべきかを選択するデザイナーによる。
【0063】
例えば、無線通信において、本発明を電話に組み込むことによって1個の携帯電話の製造コストが増加するものでないと仮定せよ。あるシステムデザイナーは、(システムに対する与えられたユーザの数に対して)各ユーザにとっての性能を改善することを選択することもできるし、あるいはそのデザイナーはシステムに対してより多くのユーザを充てることを選択することもできる(性能を現在のレベルに保持しながらシステムをインストールして操作するためのユーザ当たりのコストを下げる)。
【0064】
まとめ
まとめると、適応型フィルタにおいて用いられる部分空間射影は、ブロッキングマトリックスに対する積算の式として表すことができるものの、この射影の効率的な数値的具現化は、固定された元のN−次元空間における減算として実施されるべきである。空間を固定することによって計算上の要求が低減される。
【0065】
なお、この発明の好ましい実施例の開示は、抑制されるべき干渉の定常性あるいはガウス性の仮定に依存するものではない。唯一なされた仮定は、シミュレーションされたデータの生成と固定ステアリングベクトルsの選択における利便さのためだけである。この開示において記載された全ての式は、図3において示されたデータブロックのための二次コスト関数の最小化に関するものである。データブロックXのソースあるいは構成に関する仮定はなされていない。
【図面の簡単な説明】
【図1】 対象となる一般的なプロセッサの例示的実施例の図である。
【図2】 最小分散無ひずみレスポンス(MVDR)を決定するために概括サイドローブキャンセラ(GSC)として用いられるプロセッサの例を示す図である。
【図3】 図2においてフィルタにて処理されるべきデータの概略図である。
【図4】 本発明の好ましい実施例の概略図である。
【図5】 単周波適応型ビーム形成に適応される本発明の例を図示するための一連のプロットである。
【図6】 本発明による分析チェーン装置の概略図である。
Claims (4)
- マルチステージ的適応フィルタのための分析チェーン装置であって、
この分析チェーンは、非適応分析段階および少なくとも1つの適応分析段階を含み、
前記非適応分析段階は、
データベクトル群と、このデータベクトル群から対象となるデータを抽出するための固定ベクトルであるステアリングベクトルとを受け、ステアリングベクトルと各データベクトルとの第1の内積群を形成する第1の内積ロジック装置と、
第1の内積ロジック装置と連携し、ステアリングベクトルと第1の内積群を受け、ステアリングベクトルと第1の内積群の各内積との第1のスケーリングされた方向ベクトル群を形成する第1のベクトルスケーリングロジック装置と、
第1のベクトルスケーリングロジック装置と連携し、データベクトル群と第1のスケーリングされたベクトル群とを受け、データベクトル群と第1のスケーリングされたベクトル群との対応する要素間の第1のベクトル差群を形成する第1のベクトル差ロジック装置とを有するものであり、
前記少なくとも1つの適応分析段階は、
直前の段階と連携し、直前の段階のベクトル差群と直前の段階の内積群とを受け、直前の段階のベクトル差群と対応する直前の段階の内積群との間の現在の段階の相関方向ベクトルを形成する相関方向ベクトルロジック装置と、
直前の段階と現在の段階の相関方向ベクトルロジック装置と連携し、直前の段階のベクトル差群と現在の段階の相関方向ベクトルを受け、直前の段階の各ベクトル差と現在の段階の相関方向ベクトルとの現在の段階の内積群を形成する適応段階内積ロジック装置と、
現在の段階の相関方向ベクトルロジック装置と現在の段階の適応段階内積ロジック装置と連携し、現在の段階の内積群と現在の段階の相関方向ベクトルとを受け、現在の段階の内積群の各内積と現在の段階の相関方向ベクトルとを乗じて現在の段階のスケーリングされた方向ベクトル群を形成する適応段階ベクトルスケーリングロジック装置と、
直前の段階のベクトル差ロジック装置と現在の段階の適応段階ベクトルスケーリングロジック装置と連携し、直前の段階のベクトル差群と現在の段階のスケーリングされた方向ベクトル群とを受け、直前の段階のベクトル差群と現在の段階のスケーリングされた方向ベクトル群との対応する要素の間の現在の段階のベクトル差群を形成する適応段階ベクトル差ロジック装置とを有する、
分析チェーン装置。 - マルチステージ的適応フィルタのための分析チェーン装置の適応段階であって、
直前の段階と連携し、直前の段階のベクトル差群と直前の段階の内積群とを受け、直前の段階のベクトル差群と対応する直前の段階の内積群との間の現在の段階の相関方向ベクトルを形成する相関方向ベクトルロジック装置と、
直前の段階と現在の段階の相関方向ベクトルロジック装置と連携し、直前の段階のベクトル差群と現在の段階の相関方向ベクトルを受け、直前の段階の各ベクトル差と現在の段階の相関方向ベクトルとの現在の段階の内積群を形成する適応段階内積ロジック装置と、
現在の段階の相関方向ベクトルロジック装置と現在の段階の適応段階内積ロジック装置と連携し、現在の段階の内積群と現在の段階の相関方向ベクトルとを受け、現在の段階の内積群の各内積と現在の段階の相関方向ベクトルとを乗じて現在の段階のスケーリングされた方向ベクトル群を形成する適応段階ベクトルスケーリングロジック装置と、
直前の段階のベクトル差ロジック装置と現在の段階の適応段階ベクトルスケーリングロジック装置と連携し、直前の段階のベクトル差群と現在の段階のスケーリングされた方向ベクトル群とを受け、直前の段階のベクトル差群と現在の段階のスケーリングされた方向ベクトル群との対応する要素の間の現在の段階のベクトル差群を形成する適応段階ベクトル差ロジック装置と
を有する分析チェーン装置の適応段階。 - 複数の段階を含む適応フィルタであって、
前記複数の段階のうちの少なくとも1つの段階は、適応分析部を含み、
前記適応分析部は、
直前の段階からの入力ベクトル群および入力スカラー群の間の相関方向ベクトルを形成する相関方向ベクトルロジック装置と、
各入力ベクトルと前記相関方向ベクトルとの内積群を形成する内積ロジック装置と、
各内積および前記相関方向ベクトルからスケーリングされた方向ベクトル群を形成するベクトルスケーリングロジック装置と、
前記入力ベクトル群と前記スケーリングされた方向ベクトル群との対応する要素間のベクトル差群を形成するベクトル差ロジック装置とを含む、
適応フィルタ。 - 前記相関方向ベクトルロジック装置は、直前の段階からの入力ベクトル群(x i-1 )および入力スカラー群(d i-1 )の相関(r i-1 )およびその相関の大きさ(δ i )を計算し、当該相関(r i-1 )の単位ノルムを求めることによって正規化された相関方向ベクトル(h i )を形成するものであり、
前記適応分析部を含む少なくとも1つの段階は、さらに、合成部を含み、
前記合成部は、
直後の段階からの出力(ε i )の自乗平均(ξ i )を形成する手段と、
前記相関の大きさ(δ i )を前記自乗平均(ξ i )で除して、ウエイト(w i )を形成する手段と、
直後の段階からの出力(ε i )および前記ウエイト(w i )を乗じて、スケーリングされた出力(ε i ・w i )を形成する手段と、
前記入力スカラー群(d i-1 )と前記スケーリングされた出力(ε i ・w i )との間の差(d i-1 −ε i ・w i )を、現在の段階の出力(ε i-1 )として形成する手段と
を含む、請求項3記載の適応フィルタ。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US22840800P | 2000-08-29 | 2000-08-29 | |
US24348000P | 2000-10-27 | 2000-10-27 | |
PCT/US2001/026592 WO2002019258A1 (en) | 2000-08-29 | 2001-08-27 | System and method for adaptive filtering |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2004507992A JP2004507992A (ja) | 2004-03-11 |
JP2004507992A5 JP2004507992A5 (ja) | 2005-09-15 |
JP3859594B2 true JP3859594B2 (ja) | 2006-12-20 |
Family
ID=26922347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002524089A Expired - Fee Related JP3859594B2 (ja) | 2000-08-29 | 2001-08-27 | 適応フィルタリングのための装置 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7120657B2 (ja) |
EP (1) | EP1323121A1 (ja) |
JP (1) | JP3859594B2 (ja) |
KR (1) | KR100829485B1 (ja) |
CN (1) | CN1211755C (ja) |
AU (1) | AU2001288409A1 (ja) |
WO (1) | WO2002019258A1 (ja) |
Families Citing this family (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002211881A1 (en) * | 2000-10-13 | 2002-04-22 | Science Applications International Corporation | System and method for linear prediction |
FR2825856B1 (fr) * | 2001-06-06 | 2003-09-12 | Nortel Networks Ltd | Procede et dispositif de traitement de signal dans un recepteur de radiocommunication a etalement de spectre |
US6787119B2 (en) * | 2002-03-13 | 2004-09-07 | Council Of Scientific And Industrial Research | Process for the preparation of dinitramidic acid and salts thereof |
US7181085B1 (en) * | 2002-04-04 | 2007-02-20 | Acorn Technologies, Inc. | Adaptive multistage wiener filter |
US8082286B1 (en) | 2002-04-22 | 2011-12-20 | Science Applications International Corporation | Method and system for soft-weighting a reiterative adaptive signal processor |
US20040081131A1 (en) | 2002-10-25 | 2004-04-29 | Walton Jay Rod | OFDM communication system with multiple OFDM symbol sizes |
US8208364B2 (en) * | 2002-10-25 | 2012-06-26 | Qualcomm Incorporated | MIMO system with multiple spatial multiplexing modes |
AU2003286602A1 (en) * | 2002-10-25 | 2004-05-25 | Science Applications International Corporation | Adaptive filtering in the presence of multipath |
US8320301B2 (en) | 2002-10-25 | 2012-11-27 | Qualcomm Incorporated | MIMO WLAN system |
US7986742B2 (en) | 2002-10-25 | 2011-07-26 | Qualcomm Incorporated | Pilots for MIMO communication system |
US7280627B2 (en) * | 2002-12-09 | 2007-10-09 | The Johns Hopkins University | Constrained data-adaptive signal rejector |
GB2402308B (en) * | 2003-05-28 | 2006-01-04 | Nokia Corp | Par reduction for edge clipper |
US7193558B1 (en) * | 2003-09-03 | 2007-03-20 | The United States Of America As Represented By The Secretary Of The Navy | Radar processor system and method |
US9473269B2 (en) | 2003-12-01 | 2016-10-18 | Qualcomm Incorporated | Method and apparatus for providing an efficient control channel structure in a wireless communication system |
US7606319B2 (en) * | 2004-07-15 | 2009-10-20 | Nokia Corporation | Method and detector for a novel channel quality indicator for space-time encoded MIMO spread spectrum systems in frequency selective channels |
US7599426B2 (en) * | 2004-07-19 | 2009-10-06 | Acorn Technologies, Inc. | Use of adaptive filters in multiple access wireless systems employing predictable signals |
US7684761B2 (en) * | 2004-11-04 | 2010-03-23 | Nokia Corporation | Closed-loop signalling method for controlling multiple transmit beams and correspondingly adapted transceiver devices |
US7756196B1 (en) * | 2005-04-04 | 2010-07-13 | Acorn Technologies, Inc. | Efficient adaptive filters for CDMA wireless systems |
US20070098088A1 (en) * | 2005-10-31 | 2007-05-03 | Chih-Yuan Lin | Equalizer applied in mimo-ofdm system and related method |
US7773032B2 (en) * | 2006-05-10 | 2010-08-10 | Research Foundation Of The City University Of New York | Methods and applications utilizing signal source memory space compression and signal processor computational time compression |
GB0712270D0 (en) * | 2007-06-22 | 2007-08-01 | Nokia Corp | Wiener filtering arrangement |
US7868819B2 (en) * | 2007-09-07 | 2011-01-11 | The Board Of Trustees Of The Leland Stanford Junior University | Arrangements for satellite-based navigation and methods therefor |
US20090122853A1 (en) * | 2007-11-12 | 2009-05-14 | Acorn Technologies, Inc. | Channel tracking methods for subspace equalizers |
US8386549B1 (en) * | 2008-03-04 | 2013-02-26 | Acorn Technologies, Inc. | Reduced complexity adaptive multistage wiener filter |
US7700416B1 (en) | 2008-04-25 | 2010-04-20 | Acorn Technologies, Inc. | Tensile strained semiconductor on insulator using elastic edge relaxation and a sacrificial stressor layer |
US7851325B1 (en) | 2008-09-12 | 2010-12-14 | Acorn Technologies, Inc. | Strained semiconductor using elastic edge relaxation, a buried stressor layer and a sacrificial stressor layer |
CN101534177B (zh) * | 2009-04-03 | 2012-05-23 | 西安交通大学 | Mimo系统中基于rumswf的低复杂度降秩均衡方法 |
US9559417B1 (en) * | 2010-10-29 | 2017-01-31 | The Boeing Company | Signal processing |
EP2842384A4 (en) * | 2012-04-26 | 2015-12-16 | Propagation Res Associates Inc | METHOD AND SYSTEM FOR USING ORTHOGONAL SPATIAL PROJECTIONS TO MITIGATE PARASITING |
US9401741B2 (en) | 2012-04-26 | 2016-07-26 | Propagation Research Associates, Inc. | Methods and systems for mitigating signal interference |
US10101445B2 (en) * | 2014-04-29 | 2018-10-16 | Research Foundation Of The City University Of New York | Power centroid radar |
US9525934B2 (en) * | 2014-12-31 | 2016-12-20 | Stmicroelectronics Asia Pacific Pte Ltd. | Steering vector estimation for minimum variance distortionless response (MVDR) beamforming circuits, systems, and methods |
US10571224B2 (en) | 2015-05-04 | 2020-02-25 | Propagation Research Associates, Inc. | Systems, methods and computer-readable media for improving platform guidance or navigation using uniquely coded signals |
CN109741401A (zh) * | 2018-12-25 | 2019-05-10 | 南通蓝岛海洋工程有限公司 | 一种基于图像恢复的用于导管架水下安装在线测量方法 |
US11018705B1 (en) | 2020-07-17 | 2021-05-25 | Propagation Research Associates, Inc. | Interference mitigation, target detection, location and measurement using separable waveforms transmitted from spatially separated antennas |
US11329705B1 (en) * | 2021-07-27 | 2022-05-10 | King Abdulaziz University | Low-complexity robust beamforming for a moving source |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US425966A (en) | 1890-04-15 | weems | ||
US3892051A (en) * | 1973-10-31 | 1975-07-01 | Gen Electric | Simulated collimation of computer generated images |
US4259661A (en) * | 1978-09-01 | 1981-03-31 | Burroughs Corporation | Apparatus and method for recognizing a pattern |
US5151919A (en) * | 1990-12-17 | 1992-09-29 | Ericsson-Ge Mobile Communications Holding Inc. | Cdma subtractive demodulation |
US5285291A (en) * | 1991-02-08 | 1994-02-08 | Adobe Systems Incorporated | Methods of assigning pixels to cells of a halftone grid |
JPH08335899A (ja) * | 1995-06-07 | 1996-12-17 | N T T Ido Tsushinmo Kk | Cdma復調回路 |
US6122628A (en) * | 1997-10-31 | 2000-09-19 | International Business Machines Corporation | Multidimensional data clustering and dimension reduction for indexing and searching |
IL123782A (en) * | 1998-03-22 | 2001-05-20 | Eci Telecom Ltd | Signal Equality |
-
2001
- 2001-08-21 US US09/933,004 patent/US7120657B2/en not_active Expired - Lifetime
- 2001-08-27 KR KR1020037001848A patent/KR100829485B1/ko not_active IP Right Cessation
- 2001-08-27 JP JP2002524089A patent/JP3859594B2/ja not_active Expired - Fee Related
- 2001-08-27 CN CNB018142303A patent/CN1211755C/zh not_active Expired - Fee Related
- 2001-08-27 AU AU2001288409A patent/AU2001288409A1/en not_active Abandoned
- 2001-08-27 WO PCT/US2001/026592 patent/WO2002019258A1/en active Application Filing
- 2001-08-27 EP EP01968135A patent/EP1323121A1/en not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
CN1447953A (zh) | 2003-10-08 |
KR20030038691A (ko) | 2003-05-16 |
AU2001288409A1 (en) | 2002-03-13 |
JP2004507992A (ja) | 2004-03-11 |
CN1211755C (zh) | 2005-07-20 |
EP1323121A1 (en) | 2003-07-02 |
WO2002019258A1 (en) | 2002-03-07 |
KR100829485B1 (ko) | 2008-05-16 |
US20020152253A1 (en) | 2002-10-17 |
US7120657B2 (en) | 2006-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3859594B2 (ja) | 適応フィルタリングのための装置 | |
US7426463B2 (en) | System and method for linear prediction | |
US5694474A (en) | Adaptive filter for signal processing and method therefor | |
Huang et al. | Reduced-rank MDL method for source enumeration in high-resolution array processing | |
Jakobsson et al. | Computationally efficient two-dimensional Capon spectrum analysis | |
Harmanci et al. | Relationships between adaptive minimum variance beamforming and optimal source localization | |
WO1997011538A9 (en) | An adaptive filter for signal processing and method therefor | |
US20030198340A1 (en) | Multistage median cascaded canceller | |
De Lamare | Adaptive reduced-rank LCMV beamforming algorithms based on joint iterative optimisation of filters | |
Zhang et al. | Spatial averaging of time-frequency distributions for signal recovery in uniform linear arrays | |
Sharman | Adaptive algorithms for estimating the complete covariance eigenstructure | |
Zhang et al. | The fractional lower order moments based ESPRIT algorithm for noncircular signals in impulsive noise environments | |
Yu et al. | Efficient eigenspace-based array signal processing using multiple shift-invariant subarrays | |
Ergul et al. | Enhancing the accuracy of the interpolations and anterpolations in MLFMA | |
Gantayat et al. | An efficient direction‐of‐arrival estimation of multipath signals with impulsive noise using satin bowerbird optimization‐based deep learning neural network | |
Pan et al. | Simplified spatial smoothing for DOA estimation of coherent signals | |
Ganage et al. | PCA assisted DTCWT denoising for improved DOA estimation of closely spaced and coherent signals | |
Viberg et al. | Two decades of statistical array processing | |
Owsley et al. | Dominant mode power spectrum estimation | |
Pinto et al. | Direction Finding Based on Multi-Step Knowledge-Aided Iterative Conjugate Gradient Algorithms | |
Guo et al. | Matched field localization based on CS-MUSIC algorithm | |
Bashar et al. | Minimum array elements and snapshots for performance analysis based on direction of Arrival (DOA) estimation methods | |
Tian et al. | Beamspace mode | |
KR20090078416A (ko) | 적응적 빔포밍 장치 및 그 방법 | |
RU2014681C1 (ru) | Адаптивная антенная решетка |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A524 | Written submission of copy of amendment under article 19 pct |
Free format text: JAPANESE INTERMEDIATE CODE: A524 Effective date: 20031225 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050912 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051129 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20060224 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20060303 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060524 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060905 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060919 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100929 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |