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.
- ID Articolo: 2510.14802
- Titolo: A Scalable MVDR Beamforming Algorithm That is Linear in the Number of Antennas
- Autori: Sanjaya Herath, Armin Gerami, Kevin Wagner, Ramani Duraiswami, Christopher A. Metzler
- Classificazione: eess.SP (Elaborazione dei Segnali)
- Data di Pubblicazione: 16 ottobre 2025 (preprint arXiv)
- Link dell'Articolo: https://arxiv.org/abs/2510.14802
La tecnica di beamforming a Risposta Minima Varianza Senza Distorsione (MVDR) è ampiamente utilizzata nei sistemi di array per sopprimere le interferenze. Tuttavia, l'applicazione dell'MVDR a array di grandi dimensioni presenta sfide computazionali significative, con una complessità computazionale che cresce cubicamente rispetto al numero di elementi antenna. Questo articolo propone un metodo di beamforming MVDR scalabile per array di grandi dimensioni. Il metodo è specificamente progettato per scenari in cui il segnale di interesse è al di sotto del floor di rumore (come il GPS), sfruttando la formula di Sherman-Morrison, l'approssimazione di Decomposizione ai Valori Singolari (SVD) a basso rango e operazioni algebriche. Attraverso questo approccio, la complessità computazionale viene ridotta da una relazione cubica al numero di antenne a una relazione lineare. L'efficienza computazionale e la precisione del beamforming del metodo proposto sono state valutate mediante simulazioni e confrontate con il metodo MVDR tradizionale. Il metodo riduce significativamente il carico computazionale mantenendo al contempo un'elevata precisione di beamforming per array di grandi dimensioni.
Il problema fondamentale affrontato da questa ricerca è la complessità computazionale del beamforming MVDR tradizionale negli array di antenne di grandi dimensioni. Nello specifico:
- Collo di bottiglia della complessità computazionale: L'MVDR tradizionale richiede il calcolo dell'inversa della matrice di covarianza, con complessità O(M³), dove M è il numero di antenne
- Requisiti di tempo reale: In ambienti dinamici è necessario aggiornare frequentemente la matrice di covarianza, rendendo l'implementazione in tempo reale difficile
- Tendenza verso array di grandi dimensioni: Nei moderni sistemi radar, sonar e comunicazioni wireless, la scala degli array di antenne continua ad aumentare (da centinaia a migliaia di antenne)
L'importanza di questo problema si manifesta in:
- Ampia applicabilità: L'MVDR è ampiamente applicato nella rilevazione di bersagli radar, nell'analisi di scene acustiche e in altri campi
- Esigenze di sviluppo tecnologico: Gli array di antenne di grandi dimensioni forniscono un'elevata risoluzione spaziale e una forte capacità di soppressione delle interferenze
- Requisiti di elaborazione in tempo reale: Molti scenari applicativi richiedono l'elaborazione del beamforming in tempo reale
I metodi esistenti in letteratura presentano le seguenti limitazioni:
- Approcci algoritmici: L'approssimazione a basso rango basata su Nyström, la decomposizione QR e altri metodi presentano ancora una complessità computazionale elevata
- Metodi distribuiti: Richiedono protocolli di comunicazione complessi e meccanismi di sincronizzazione
- Metodi di apprendimento profondo: Richiedono grandi quantità di dati di addestramento con capacità di generalizzazione limitata
- Propone un algoritmo MVDR a complessità lineare: Riduce la complessità computazionale da O(M³) a O(MK²), dove K≪M
- Integra molteplici tecniche matematiche: Combina abilmente la formula di Sherman-Morrison, l'approssimazione SVD a basso rango e operazioni algebriche
- Ottimizzazione per scenari specifici: Progettato specificamente per scenari in cui il segnale è al di sotto del floor di rumore, come le applicazioni GPS
- Mantiene elevata precisione di beamforming: Riduce significativamente la complessità computazionale mantenendo prestazioni comparabili all'MVDR tradizionale
- Fornisce un framework algoritmico completo: Presenta una procedura completa di inizializzazione, aggiornamento e calcolo dei pesi
Si consideri un array composto da M antenne isotropiche che ricevono il segnale di interesse dalla direzione θ₀ e L segnali di interferenza dalle direzioni θ₁, θ₂, ..., θL. Il vettore di segnale ricevuto al tempo t è:
xt=a(θ0)s0(t)+∑l=1La(θl)sl(t)+sn(t)
dove a(θᵢ) è il vettore di steering, s₀(t) è il segnale di interesse, sₗ(t) è il segnale di interferenza e sₙ(t) è il rumore.
La soluzione del beamformer MVDR tradizionale è:
w=a(θ0)HR−1a(θ0)R−1a(θ0)
Utilizzo della regola di aggiornamento ricorsivo con fattore di dimenticanza α:
Rn+1=αRn+(1−α)xnxnH
Utilizzo della formula di Sherman-Morrison per aggiornare l'inversa della matrice di covarianza:
Rn+1−1=αRn−1−α(1−α)α+(1−α)xnHRn−1xnRn−1xnxnHRn−1
Approssimazione di rango K della matrice di covarianza:
R≈UKDKUKH
dove U_K ∈ C^{M×K} contiene i primi K autovettori e D_K ∈ C^{K×K} contiene i primi K autovalori.
Attraverso l'approssimazione a basso rango, la regola di aggiornamento diventa:
DK,n+1−1=αDK,n−1−α(1−α)α+(1−α)xnHUKDK,n−1UKHxnDK,n−1UKHxnxnHUKDK,n−1
- Strategia di riduzione dimensionale: Attraverso l'approssimazione a basso rango con K≪M si evita la formazione esplicita della matrice di covarianza M×M
- Meccanismo di aggiornamento incrementale: Utilizza la formula di Sherman-Morrison per realizzare aggiornamenti con complessità O(MK²)
- Ottimizzazione per scenari specifici: Per scenari in cui il segnale è al di sotto del floor di rumore, K può tipicamente essere impostato a L+1
- Integrazione algoritmica: Combina organicamente SVD, formula di Sherman-Morrison e aggiornamento ricorsivo
- Configurazione dell'array: Array lineare uniforme (ULA), spaziatura tra antenne di mezza lunghezza d'onda
- Numero di antenne: M = 50, 75, 100 (esperimenti principali), esteso a 500 (test di complessità)
- Impostazione del segnale:
- Segnale bersaglio: distanza 10 km, velocità tangenziale 500 m/s
- Segnale trasmesso: impulso chirp lineare, larghezza di banda 300 MHz, durata dell'impulso 100 ms
- SINR: -10 dB
- Frequenza di campionamento: 1 MHz
- Parametri dell'algoritmo:
- Fattore di dimenticanza α = 0,99
- Dimensionalità a basso rango K = 10
- Numero di snapshot: 1000
- Larghezza del lobo principale (MLW): Larghezza angolare del lobo principale del diagramma di radiazione
- Livello dei lobi laterali (SLL): Livello di potenza dei lobi laterali relativo al lobo principale
- Profondità della tacca nulla: Livello di potenza della tacca nulla relativo al lobo principale
- Tempo di calcolo: Tempo totale per eseguire 10.000 passi temporali
- Guadagno SINR: Rapporto tra SINR di uscita e SINR di ingresso
- MVDR tradizionale: Metodo standard di inversione della matrice di covarianza
- Confronto del tempo di esecuzione: Testato su processore AMD EPYC 7443 a 24 core
| M | L | MLW (°) | | Profondità Tacca Nulla (dB) | | SLL (dB) | |
|---|
| | MVDR | Prop. | MVDR | Prop. | MVDR | Prop. |
| 50 | 1 | 2,27 | 2,32 | -49,56 | -42,44 | -13,02 | -9,14 |
| 75 | 2 | 1,36 | 1,33 | -41,02 | -34,84 | -12,75 | -11,05 |
| 100 | 3 | 1,06 | 1,08 | -45,83 | -41,78 | -11,37 | -12,47 |
- MVDR tradizionale: Complessità O(M³), il tempo di esecuzione cresce con M³
- Metodo proposto: Complessità O(MK²), il tempo di esecuzione cresce linearmente con M
- Miglioramento delle prestazioni: Per array con M=500, il tempo di calcolo si riduce di diversi ordini di grandezza
I risultati sperimentali mostrano che il metodo proposto è comparabile al metodo MVDR tradizionale nei seguenti aspetti:
- Puntamento del lobo principale: Indirizza correttamente il lobo principale verso la direzione del bersaglio
- Formazione della tacca nulla: Forma tacche nulle efficaci nelle direzioni dei segnali di interferenza
- Forma complessiva del diagramma di radiazione: Altamente coerente con l'MVDR tradizionale
Attraverso simulazioni con 100.000 passi temporali si è scoperto che:
- Decadimento del SINR: L'uso prolungato comporta una diminuzione del guadagno SINR
- Effetto della reinizializzazione: Dopo la reinizializzazione al passo 50.000, il guadagno SINR si recupera
- Costo della reinizializzazione: La complessità O(M²) della reinizializzazione rimane inferiore al metodo tradizionale O(M³)
- Approssimazione a basso rango di Nyström: Utilizza approssimazione a basso rango per ridurre il carico computazionale e di memoria
- Metodo di decomposizione QR: Tracciamento dinamico dei cambiamenti di voce e rumore
- MVDR-SMI: Implementazione ricorsiva utilizzando decomposizione di Cholesky e trasformazioni di Householder
- Algoritmi di passaggio di messaggi: Realizzazione del beamforming distribuito attraverso la comunicazione tra nodi locali
- Metodo ADMM: Distribuzione del carico computazionale su più processori
- Apprendimento per rinforzo avversariale profondo: Miglioramento delle capacità di beamforming per MIMO di grandi dimensioni
- Reti neurali convoluzionali: Riduzione della complessità di addestramento e calibrazione di array di antenne di grandi dimensioni
- Riduzione significativa della complessità computazionale: Riduzione da O(M³) a O(MK²), realizzando scalabilità lineare
- Mantenimento di elevata precisione di beamforming: Prestazioni comparabili all'MVDR tradizionale in metriche chiave come larghezza del lobo principale e livello dei lobi laterali
- Applicabilità a applicazioni in tempo reale: Fornisce una soluzione pratica per il beamforming MVDR in tempo reale su array di grandi dimensioni
- Restrizioni sugli scenari applicabili: Specificamente progettato per scenari in cui il segnale è al di sotto del floor di rumore, con prestazioni scadenti quando il segnale è al di sopra del floor di rumore
- Decadimento delle prestazioni a lungo termine: Richiede reinizializzazione periodica per mantenere il guadagno SINR
- Costo di inizializzazione: Il calcolo SVD iniziale richiede ancora complessità O(M³)
- Estensione a scenari in cui il segnale è al di sopra del floor di rumore
- Metodi di inizializzazione più efficienti: Come metodi SVD randomizzati
- Strategie di reinizializzazione adattive: Attivazione automatica della reinizializzazione in base al decadimento delle prestazioni
- Forte innovazione teorica: Combina abilmente molteplici strumenti matematici per risolvere problemi pratici
- Elevato valore pratico: Risolve il collo di bottiglia critico dell'MVDR per array di grandi dimensioni
- Esperimenti completi: Verifica da molteplici prospettive incluse precisione, complessità e prestazioni a lungo termine
- Completezza algoritmica: Fornisce procedura algoritmica completa e dettagli di implementazione
- Scenari applicativi limitati: Applicabile solo a condizioni di segnale specifiche
- Analisi teorica insufficiente: Mancanza di garanzie teoriche su convergenza e stabilità
- Guida sulla scelta dei parametri: Mancanza di guida teorica sulla scelta del valore K
- Assenza di verifica su hardware reale: Solo risultati di simulazione, mancanza di verifica su piattaforme hardware reali
- Contributo accademico: Fornisce nuovi approcci risolutivi per il beamforming su array di grandi dimensioni
- Valore ingegneristico: Potenziale applicazione in radar, sonar, comunicazioni wireless e altri campi
- Riproducibilità: Descrizione algoritmica chiara, facile da riprodurre e migliorare
- Ricezione GPS: Scenario tipico di segnale al di sotto del floor di rumore
- Rilevamento di segnali deboli: Applicazioni che richiedono forte soppressione delle interferenze
- Array di antenne di grandi dimensioni: Sistemi con centinaia o migliaia di antenne
- Applicazioni sensibili alla latenza computazionale: Applicazioni sensibili al ritardo di calcolo
L'articolo cita 21 articoli correlati, coprendo teoria fondamentale del beamforming, elaborazione di array di grandi dimensioni, algoritmi SVD e altri aspetti, fornendo una base teorica solida per la ricerca.
Valutazione Complessiva: Questo è un articolo di notevole valore pratico nel campo dell'elaborazione dei segnali. Attraverso tecniche matematiche ingegnose risolve il collo di bottiglia computazionale del beamforming MVDR per array di grandi dimensioni. Sebbene presenti limitazioni come restrizioni sugli scenari applicativi, fornisce contributi preziosi allo sviluppo del campo.