2025-11-21T23:10:16.385556

A Scalable MVDR Beamforming Algorithm That is Linear in the Number of Antennas

Herath, Gerami, Wagner et al.

The Minimum Variance Distortionless Response (MVDR) beamforming technique is widely applied in array systems to mitigate interference. However, applying MVDR to large arrays is computationally challenging; its computational complexity scales cubically with the number of antenna elements. In this paper, we introduce a scalable MVDR beamforming method tailored for massive arrays. Our approach, which is specific to scenarios where the signal of interest is below the noise floor (e.g.,~GPS), leverages the Sherman-Morrison formula, low-rank Singular Value Decomposition (SVD) approximations, and algebraic manipulation. Using our approach, we reduce the computational complexity from cubic to linear in the number of antennas. We evaluate the proposed method through simulations, comparing its computational efficiency and beamforming accuracy with the conventional MVDR approach. Our method significantly reduces the computational load while maintaining high beamforming accuracy for large-scale arrays. This solution holds promise for real-time applications of MVDR beamforming in fields like radar, sonar, and wireless communications, where massive antenna arrays are proliferating.

academic

アンテナ数に対して線形なスケーラブルMVDRビームフォーミングアルゴリズム

基本情報

論文ID: 2510.14802
タイトル: A Scalable MVDR Beamforming Algorithm That is Linear in the Number of Antennas
著者: Sanjaya Herath, Armin Gerami, Kevin Wagner, Ramani Duraiswami, Christopher A. Metzler
分類: eess.SP（信号処理）
発表日時: 2025年10月16日（arXiv プレプリント）
論文リンク: https://arxiv.org/abs/2510.14802

要約

最小分散無歪応答（MVDR）ビームフォーミング技術は、アレイシステムにおける干渉抑制に広く応用されている。しかし、大規模アレイへのMVDR適用は計算上の課題があり、計算複雑度はアンテナ素子数に対して3乗の関係にある。本論文は大規模アレイに対応するスケーラブルなMVDRビームフォーミング方法を提案する。本方法は、目的信号がノイズフロア以下の場面（GPS等）を想定し、Sherman-Morrison公式、低ランク特異値分解（SVD）近似、および代数操作を活用している。この方法により、計算複雑度をアンテナ数の3乗から線形関係に削減する。シミュレーション評価を通じて、提案方法の計算効率とビームフォーミング精度を検証し、従来のMVDR方法と比較した。本方法は計算負荷を大幅に削減しながら、大規模アレイにおける高いビームフォーミング精度を維持する。

研究背景と動機

問題定義

本研究が解決する中核的な問題は、従来のMVDRビームフォーミングが大規模アンテナアレイにおいて直面する計算複雑度の問題である。具体的には以下の通りである：

計算複雑度のボトルネック：従来のMVDRは共分散行列の逆行列計算を必要とし、その複雑度はO(M³)である。ここでMはアンテナ数である
リアルタイム性の要求：動的環境では共分散行列の頻繁な更新が必要となり、リアルタイム実装が困難になる
大規模アレイの傾向：現代のレーダ、ソナー、無線通信システムにおいて、アンテナアレイの規模が増大している（数百から数千のアンテナ）

重要性分析

本問題の重要性は以下の点に表れている：

応用の広範性：MVDRはレーダ目標検出、音場景分析等の分野で広く応用されている
技術発展の必要性：大規模アンテナアレイは高い空間分解能と強力な干渉抑制能力を提供する
リアルタイム処理の要求：多くの応用場面ではリアルタイムビームフォーミング処理が要求される

既存方法の限界

文献における既存方法には以下の限界がある：

アルゴリズム的手法：Nyström法に基づく低ランク近似、QR分解等は依然として計算複雑度が高い
分散型手法：複雑な通信プロトコルと同期機構が必要である
深層学習手法：大量の訓練データが必要であり、汎化能力が限定的である

核心的貢献

線形複雑度のMVDRアルゴリズムを提案：計算複雑度をO(M³)からO(MK²)に削減。ここでK≪M
複数の数学的技法を統合：Sherman-Morrison公式、低ランクSVD近似、代数操作を巧妙に融合
特定の場面に最適化：GPS応用等、信号がノイズフロア以下の場面に特化した設計
高いビームフォーミング精度を維持：計算複雑度を大幅に削減しながら、従来のMVDRと同等の性能を保持
完全なアルゴリズムフレームワークを提供：初期化、更新、重み計算の完全なプロセスを提示

方法の詳細

タスク定義

M個の等方性アンテナで構成されるアレイが、目標方向θ₀からの目的信号と方向θ₁,θ₂,...,θLからのL個の干渉信号を受信する場合を考える。時刻tで受信される信号ベクトルは以下の通りである：

$x_t = a(\theta_0)s_0(t) + \sum_{l=1}^{L} a(\theta_l)s_l(t) + s_n(t)$

ここでa(θᵢ)は操舵ベクトル、s₀(t)は目的信号、sₗ(t)は干渉信号、sₙ(t)はノイズである。

従来のMVDRビームフォーマの解は以下の通りである： $w = \frac{R^{-1}a(\theta_0)}{a(\theta_0)^H R^{-1}a(\theta_0)}$

モデルアーキテクチャ

1. 再帰的共分散行列更新

忘却因子αを用いた再帰更新規則： $R_{n+1} = \alpha R_n + (1-\alpha)x_n x_n^H$

2. Sherman-Morrison公式の適用

Sherman-Morrison公式を用いて共分散行列の逆行列を更新： $R_{n+1}^{-1} = \frac{R_n^{-1}}{\alpha} - \frac{(1-\alpha)}{\alpha} \frac{R_n^{-1}x_n x_n^H R_n^{-1}}{\alpha + (1-\alpha)x_n^H R_n^{-1}x_n}$

3. 低ランクSVD近似

共分散行列をK階近似： $R \approx U_K D_K U_K^H$

ここでU_K ∈ C^{M×K}は上位K個の固有ベクトルを含み、D_K ∈ C^{K×K}は上位K個の固有値を含む。

4. 次元削減更新規則

低ランク近似を通じて、更新規則は以下となる： $D_{K,n+1}^{-1} = \frac{D_{K,n}^{-1}}{\alpha} - \frac{(1-\alpha)}{\alpha} \frac{D_{K,n}^{-1}U_K^H x_n x_n^H U_K D_{K,n}^{-1}}{\alpha + (1-\alpha)x_n^H U_K D_{K,n}^{-1}U_K^H x_n}$

技術的革新点

次元削減戦略：K≪Mの低ランク近似により、M×M共分散行列の明示的な形成を回避
漸進的更新機構：Sherman-Morrison公式を利用してO(MK²)複雑度の更新を実現
特定場面への最適化：信号がノイズフロア以下の場面では、Kは通常L+1に設定可能
アルゴリズム統合：SVD、Sherman-Morrison公式、再帰更新を有機的に結合

実験設定

シミュレーションパラメータ

アレイ構成：均一線形アレイ（ULA）、アンテナ間隔は半波長
アンテナ数：M = 50, 75, 100（主要実験）、500まで拡張（複雑度テスト）
信号設定：
- 目標信号：10km距離、500m/s接線速度
- 送信信号：線形周波数変調パルス、帯域幅300MHz、パルス持続時間100ms
- SINR：-10dB
- サンプリングレート：1MHz
アルゴリズムパラメータ：
- 忘却因子α = 0.99
- 低ランク次元K = 10
- スナップショット数：1000