Il framework proposto per la modellizzazione del canale in campo prossimale e il beamforming nei sistemi 6G abilitati da XL-IMS combina la modellizzazione basata sulla fisica con l'ottimizzazione mediante apprendimento automatico, superando i limiti dei tradizionali approcci in campo lontano. Viene sviluppato e convalidato utilizzando il dataset DeepMIMO un modello di canale in campo prossimale che incorpora la propagazione dell'onda sferica e la non stazionarietà spaziale. L'informazione del canale è rappresentata come un grafo per catturare le dipendenze spaziali tra gli elementi della metasuperficie e gli utenti, e viene impiegata una rete end-to-end basata su GNN con attenzione per apprendere i pesi di beamforming necessari al focalizzamento del fascio in campo prossimale. A differenza dei metodi tradizionali di ottimizzazione iterativa, il framework basato sull'apprendimento prevede direttamente le soluzioni di beamforming, riducendo la complessità computazionale pur mantenendo un'elevata accuratezza e adattabilità a condizioni di canale dinamiche. Le prestazioni vengono valutate in termini di velocità raggiungibile, guadagno di beamforming ed efficienza energetica, mostrando risultati superiori rispetto ai metodi tradizionali di beamforming in campo lontano. Figura 1 illustra il framework complessivo.
Di seguito sono riportati i passaggi proposti:
Passaggio 1: Configurare l'ambiente DeepMIMO
La piattaforma DeepMIMO v3.x viene inizialmente configurata tramite l'ambiente esterno O1, che utilizza una frequenza portante di 28 GHz e una larghezza di banda di 100 MHz. Nella simulazione, l'XL-IMS è rappresentato da un array planare uniforme 16×16 composto da 256 elementi riflettenti, distanziati tra loro di λ/2. Il modello di propagazione scelto è un'onda sferica con percorsi LoS e NLoS.
Passaggio 2: Definire la geometria XL-IMS e le posizioni degli utenti
Progettare l'XL-IMS come una metasuperficie planare 16×16 che comprende 256 elementi riflettenti riconfigurabili. I dispositivi terminali utente sono stati distribuiti casualmente entro il campo vicino, da 0,5 m a 10 m dall'XL-IMS, in base al criterio della distanza di Rayleigh.
Passaggio 3: Generare canali in campo prossimale utilizzando la propagazione di onde sferiche
Modellare i canali tra ciascun elemento della metasuperficie riconfigurabile e l'utente utilizzando la propagazione di onde sferiche, che tiene conto della variazione di fase indotta dalla distanza.
Passaggio 4: Generare dati del canale in campo prossimale
I terminali utente sono stati distribuiti in modo uniforme all'interno della zona in prossimità (0,5-10 m) dal pannello XL-IMS. Sono stati utilizzati sia i percorsi di propagazione LoS che NLoS generati dallo scenario esterno DeepMIMO O1 per generare i canali.
Passaggio 5: Preparare il set di dati per l'addestramento e il test
Il set di dati generato comprendeva 20.000 campioni, di cui 16.000 (80%) sono stati utilizzati per l'addestramento e 4.000 (20%) per i test.
Passaggio 6: Normalizzare le caratteristiche di input e preparare gli input del modello
L'ultima fase del processo di preparazione ha comportato l'applicazione della normalizzazione Min-Max a tutte le caratteristiche in ingresso, come l'ampiezza del canale, la fase, la distanza e le coordinate dell'utente. Questo processo ha normalizzato i dati e li ha convertiti in tensori, che sono stati inseriti nella rete GNN con attenzione. Il modello è stato ottimizzato utilizzando l'ottimizzatore Adam con un tasso di apprendimento di 0,001 e 50 epoche.
Passaggio 7: Applicare il meccanismo di attenzione
Aggiungere un meccanismo di attenzione per considerare le informazioni globali del grafo e catturare le dipendenze a lungo raggio tra i nodi del grafo. Ciò consente al modello di concentrarsi su caratteristiche spaziali importanti.
Passaggio 8: Prevedere i pesi della formazione del fascio
Predire i pesi del beamforming e i valori di sfasamento utilizzando l'embedding basato sull'attenzione delle GNN. Le complesse, laboriose e iterative procedure di beamforming vengono sostituite da un metodo predittivo.
Passaggio 9: Valutare le metriche di prestazione
Per valutare l'efficacia della comunicazione e del beamforming, calcolare indicatori chiave di prestazione (KPI) come la velocità raggiungibile, il guadagno del beamforming, l'efficienza spettrale, l'efficienza energetica e il SINR.
Passo 10: Confronto con i metodi di riferimento
Confrontare l'algoritmo proposto con gli algoritmi tradizionali, inclusi ZF, MMSE, SDR e AO.
Definizione del modello di sistema
Il sistema proposto considera una rete wireless 6G abilitata da un XL-IMS costituito da 256 elementi riflettenti riconfigurabili che controllano la propagazione dell'onda elettromagnetica mediante spostamenti di fase regolabili. Poiché gli utenti operano nella regione del campo vicino, l'ipotesi convenzionale di onda piana del campo lontano diventa non valida, rendendo necessaria la considerazione della propagazione dell'onda sferica, dell'attenuazione dipendente dalla distanza e della non stazionarietà spaziale. Il confine tra le regioni del campo vicino e del campo lontano è determinato dalla distanza di Rayleigh, che dipende dalle dimensioni della metasuperficie e dalla lunghezza d'onda del segnale portante, e indica la transizione dagli effetti del fronte d'onda sferico alla propagazione dell'onda piana. Ciò è descritto come segue:
(1)
Qui, D è la dimensione della metasuperficie, λ è la lunghezza d'onda e R è la distanza di Rayleigh. Un fenomeno in campo prossimale si verifica nel canale di comunicazione quando un utente si trova in una posizione più vicina a R rispetto al suo valore reale. Mentre il segnale attraversa diversi componenti della metasuperficie, subisce variazioni distinte di fase e di ampiezza.
Il canale dalla n-esima meta-superficie all'utente utilizza un modello di onda sferica per la propagazione in campo vicino, rappresentato dall'equazione del coefficiente del canale come segue:
(2)
In questo caso, α indica il parametro di attenuazione del percorso, λ indica la lunghezza d'onda della portante e dn è la distanza dell'utente dal componente.
In questo caso, βn è il fattore di attenuazione del percorso per l'elemento XL-IMS nesimo e può essere espresso come βn = dn−α, dove α rappresenta l'esponente di attenuazione del percorso. Si noti che dn è la distanza euclidea tra la posizione dell'utente e l'elemento nesimo della metasuperficie. Le variazioni del canale rispetto allo spazio sono causate dalle differenze significative di dn tra diversi elementi della metasuperficie, a differenza della situazione riscontrata nei canali in campo lontano. In generale, il vettore del canale può essere descritto come segue:
(3)
Il segnale in ingresso ricevuto dall'utente dipende dalla sintonizzazione di fase del sistema XL-IMS. Si consideri θ = [θ1, θ2, …., θN] per rappresentare lo sfasamento realizzato utilizzando la metasuperficie. In tal caso, il segnale in ingresso può essere espresso come:
(4)
dove Φ = diag(ejθ1, ….., ejθN) è il segnale trasmesso, w è il vettore di precodifica a codifica base e n rappresenta il rumore. Il vettore del canale in campo prossimo è definito come h ∈ ℂN×1, la matrice diagonale dei cambi di fase dell'XL-IMS come Φ ∈ ℂN×N, il vettore di beamforming come x ∈ ℂN×1 e il segnale ricevuto come y ∈ ℂ. Il termine di rumore additivo è definito come n ∼ CN (0, σ2). Queste dimensioni mantengono la correttezza matematica della formulazione.
In questo caso, ws è il vettore di precodifica della stazione base utilizzato per precodificare il simbolo trasmesso (s). Durante il confronto delle prestazioni, presso la stazione base vengono implementati comuni approcci di precodifica digitale come ZF e MMSE, utilizzando realizzazioni simili del canale in campo prossimo. Tuttavia, il nuovo approccio mira a ottimizzare la matrice di sfasamento XL-IMS (θ) attraverso il meccanismo di apprendimento con attenzione GNN.
Per concentrare l'energia nella regione del campo vicino, deve essere realizzabile una progettazione di θ in cui si verifica un'interferenza costruttiva per l'utente. La strategia di beamforming mira a massimizzare l'efficienza spettrale o la potenza del segnale ricevuto. La Figura 2 descrive l'architettura del modello di sistema di comunicazione nel campo vicino basato sull'approccio XL-IMS.
Modellizzazione del canale in campo prossimo
Il modello proposto per le reti di comunicazione 6G basate su XL-IMS nella regione del campo vicino è stato progettato specificamente per imitare le condizioni di propagazione prossime alla distanza di Rayleigh rispetto alla metasuperficie. Nella regione del campo vicino, la propagazione è sferica, con valori di fase variabili e parametri spaziali non omogenei lungo la metasuperficie, a differenza dei modelli convenzionali di propagazione nel campo lontano che assumono onde piane. Nello studio attuale, è stato utilizzato un modello fisico di onda sferica per progettare una matrice di risposta all'impulso del canale accurata, che tiene conto delle posizioni fisiche sia degli elementi della metasuperficie sia dei terminali utente.
Rappresentazione del canale basata sulla distanza
La distanza tra ciascun componente della metasuperficie e l'utente è un fattore cruciale nella costruzione del modello del canale nel regime del campo vicino. La lunghezza pn tra le posizioni dell'utente pu e il n-il componente XL-IMS è fornito da:
(5)
Questo valore varia notevolmente lungo la metasuperficie a causa della sua area estremamente grande, determinando tempi di propagazione e spostamenti di fase disuguali che devono essere considerati nel modello.
Modello del canale a onda sferica
Il coefficiente del canale tra l'elemento n-esimo e l'utente è determinato utilizzando il modello ad onda sferica descritto nell'Equazione (2), che tiene conto delle variazioni di fase e di ampiezza dipendenti dalla distanza. Per la generazione del dataset e l'addestramento del modello, i coefficienti del canale sono ottenuti dal framework di ray-tracing DeepMIMO, che incorpora sia la propagazione in vista diretta (LoS) sia quella in vista indiretta (NLoS) basata su strutture ambientali realistiche. A partire da quanto sopra, è possibile costruire la matrice di risposta del canale a partire dalle risposte del canale per tutti gli N elementi della superficie metamatematica, come indicato nell'Equazione (3).
Non stazionarietà spaziale dei canali in campo prossimale e importanza del modello
A differenza dei canali in campo lontano, i canali in campo vicino mostrano una non stazionarietà spaziale, in cui i segnali ricevuti variano lungo la metasuperficie, permettendo un'illuminazione del fascio dipendente dalla posizione. Questa informazione spaziale è incorporata nella matrice di risposta del canale ed è catturata dal modello GNN-attention. Supplementary Figure 1 illustra il framework proposto per la modellizzazione del canale in campo vicino per il sistema di comunicazione 6G basato su XL-IMS e la generazione di dati di canale realistici per la formazione del fascio basata su apprendimento profondo.
Generazione del set di dati
La generazione dei dati è stata eseguita utilizzando il framework DeepMIMO, che fornisce simulazioni di canale realistiche basate sul ray-tracing, incorporando la geometria ambientale, la coerenza spaziale e la mobilità dell'utente. Le risposte del canale per più posizioni utente in campo prossimo sono state generate utilizzando lo scenario esterno O1 a una frequenza portante di 28 GHz e una larghezza di banda di 100 MHz. L'XL-IMS era composto da 256 elementi riflettenti distribuiti uniformemente, e le posizioni degli utenti sono state selezionate all'interno della regione di campo prossimo in base al criterio di Rayleigh. I coefficienti del canale, i guadagni di percorso e le informazioni spaziali sono stati estratti utilizzando la libreria Python DeepMIMO. Tabella supplementare 1 riassume i parametri di simulazione e riproducibilità. Il modello è stato addestrato mediante apprendimento supervisionato, utilizzando caratteristiche dello stato del canale come ingressi e pesi di beamforming ottimizzati basati sulla fisica come etichette obiettivo.
Creazione delle caratteristiche di input e struttura del set di dati
Successivamente, i dati grezzi dei canali sono stati convertiti in caratteristiche da utilizzare come input per modelli di apprendimento profondo. Queste caratteristiche includevano le ampiezze e le fasi delle risposte dei canali come componenti separati. Altri parametri di input potevano essere inclusi in base alle esigenze del problema. Ad esempio, le coordinate degli utenti, le distanze e alcuni parametri spaziali potevano essere considerati anch'essi come caratteristiche. Il dataset è stato suddiviso in un insieme di addestramento e un insieme di test in un rapporto 80:20.
Preparazione del set di dati per l'addestramento
L'ultimo passaggio prevedeva la normalizzazione del dataset preparato e la sua conversione in tensori adatti all'ingresso nell'architettura GNN-attention per l'addestramento. Il dataset di addestramento è stato utilizzato per apprendere la mappatura dallo stato del canale ai pesi ottimali di beamforming, mentre il dataset di test è stato impiegato per valutare l'efficacia del modello. Questo dataset accuratamente preparato ha reso il processo di apprendimento più semplice ed efficace. Supplementary Table 2 mostra i dettagli descrittivi del dataset di DeepMIMO.
Rappresentazione delle caratteristiche
La fase di rappresentazione delle caratteristiche ha trasformato i dati del canale a campo vicino di valore complesso in un formato adatto all'apprendimento profondo. Le risposte del canale sono state scomposte in caratteristiche che catturano le proprietà spaziali e di propagazione, preservando al contempo le informazioni essenziali del campo vicino per l'apprendimento. Figura supplementare 2 illustra il dataset DeepMIMO utilizzato in questo processo. I dati del canale provenienti dall'ambiente a campo vicino sono stati rappresentati come numeri complessi hn, con componenti reale e immaginaria che rappresentano l'informazione di fase e ampiezza dei coefficienti del canale. Per consentirne un uso efficace negli algoritmi di apprendimento automatico, i coefficienti sono stati scomposti nelle loro parti reale (Rehn) e immaginaria (Imhn). I dati possono inoltre essere rappresentati in termini di modulo hn e fase dei coefficienti ∠hn.
Estrazione di caratteristiche spaziali
Oltre al coefficiente del canale, sono state incluse caratteristiche spaziali per tenere conto della natura geometrica dell'ambiente in campo prossimo. Queste comprendevano la distanza dn da ciascun elemento della metasuperficie al ricevitore, insieme alle coordinate del ricevitore (x, y, z). In questo modo, il modello ha considerato il fatto che le caratteristiche del segnale possono dipendere dalla posizione spaziale della coppia trasmettitore-ricevitore, una proprietà importante per la propagazione in campo prossimo.
Normalizzazione e preparazione dell'input del modello
Come passo preliminare prima di fornire al modello i suoi input, è stato effettuato un processo di normalizzazione per evitare instabilità numerica e accelerare la convergenza. I vettori strutturati sono stati suddivisi in set di dati di addestramento e di test, indicati rispettivamente con (Xtrain, Xtest). In questo modo, gli input sono stati successivamente forniti ai modelli GNN e ai modelli con attenzione.
Sulla base della rappresentazione delle caratteristiche mostrata nella Figura Supplementare 2, sono stati utilizzati diversi tipi di caratteristiche in questo dataset, inclusi i coefficienti di canale hn, che indicavano la risposta del segnale complesso di tutti gli elementi della metasuperficie. Un'altra caratteristica era la distanza dn, che in questo caso indicava gli effetti di propagazione. Sebbene le componenti immaginarie e reali forniscano una descrizione numerica del segnale, le parti di ampiezza e fase forniscono rispettivamente una descrizione esplicita dell'intensità del segnale e della variazione di fase. La coordinata della posizione dell'utente (x, y, z) era un'altra caratteristica importante che migliorava la consapevolezza spaziale del modello.
Progettazione basata sul deep learning mediante GNN per la modellizzazione della dipendenza spaziale integrata con un meccanismo ibrido di attenzione per l'apprendimento del contesto globale
Questo algoritmo di apprendimento è stato sviluppato per addestrare in modo efficiente il complesso canale XL-IMS nel campo vicino, con le sue caratteristiche spaziali e di propagazione, utilizzando reti neurali grafiche (GNN) e un meccanismo di attenzione. Poiché gli elementi della metasuperficie sono fortemente correlati tra loro e le loro interazioni sono non uniformi all'interno della zona di campo vicino, una rete completamente connessa non è adatta a rappresentare correttamente il canale. Pertanto, in questo approccio, la metasuperficie è rappresentata come un grafo, in cui ciascun elemento funge da nodo e le interazioni tra essi costituiscono gli archi. La GNN ha catturato le dipendenze spaziali localizzate tra nodi adiacenti, mentre il meccanismo di attenzione ha migliorato l'apprendimento delle dipendenze globali.
Rappresentazione grafica di XL-IMS
L'XL-IMS è creato come un grafo G = (V, E), in cui ogni vertice vn ∈ V rappresenta un elemento della meta-superficie, mentre i lati eij ∈ E rappresentano le relazioni spaziali tra gli elementi. Ogni vertice possiede un vettore di caratteristiche xn, ottenuto tramite la fase di rappresentazione delle caratteristiche. La matrice di adiacenza A descrive le connessioni tra i vertici che sono vicini o entro determinati limiti di distanza. Tale formulazione a grafo viene utilizzata per incorporare le informazioni spaziali presenti nella propagazione in campo prossimale.
Per garantire la riproducibilità, la connettività del grafo viene determinata utilizzando l'approccio dei K-Vicini più Prossimi con k impostato a 8. Lasciare pi e pj rappresentano i vettori di posizione di i-e e j-th cellule XL-IMS, rispettivamente. La distanza tra due vertici qualsiasi è misurata mediante la metrica euclidea come segue: dij = ||pi − pj||2. ci sarà un arco tra i vertici i e j se vertice j è incluso nel gruppo degli otto vicini più prossimi al vertice i. Pertanto, la matrice delle adiacenze viene costruita secondo la seguente regola: Aij = 1 se j ∈ Nk(i) e Aij = 0 altrimenti.
Nella GNN proposta, gli aggiornamenti per ogni nodo vengono calcolati sulla base delle informazioni provenienti dagli altri nodi adiacenti. In particolare, la propagazione del livello l può essere espressa come segue:
(6)
dove hi(l) rappresenta il vettore di caratteristiche del vertice i nel livello l, N(i) è il vicinato del vertice i, W(l) è la matrice di pesi apprendibile e σ(·) indica una funzione di attivazione non lineare.
Il modello GNN-attention era composto da tre livelli di convoluzione su grafi con dimensioni nascoste impostate a 128, e come funzione di attivazione è stata utilizzata l'unità lineare rettificata (ReLU). La generazione del grafo è stata eseguita mediante il metodo dei k-vicini più vicini (k=8). Ogni caratteristica XL-IMS funge da nodo, mentre i collegamenti sono creati in base alla distanza euclidea tra nodi adiacenti. Per estrarre dipendenze spaziali globali, abbiamo utilizzato una rete di auto-attenzione multi-testa con 4 teste. La probabilità di dropout è stata impostata a 0,3 per prevenire l'overfitting.
Meccanismo di attenzione
Per modellare le dipendenze a lungo raggio che esistevano al di fuori del raggio di vicinato locale, è stata introdotta l'attenzione nel modello, il cui funzionamento è il seguente:
(7)
dove eij = LeakyReLU(aT[Whi | Whj]) rappresenta l'importanza dei nodi i e j. Con questo modello, è possibile assegnare diversi livelli di importanza ai nodi in base ai loro ruoli, facilitando così la comprensione del contesto globale sull'intera metasuperficie. Il framework ibrido cattura sia informazioni spaziali locali che globali, e le caratteristiche apprese vengono elaborate attraverso strati completamente connessi per generare i pesi di formazione del fascio. Come mostrato nelle Figure Supplementari 3 e 4, la GNN modella le dipendenze spaziali locali mediante il passaggio di messaggi e l'aggregazione delle caratteristiche, mentre il meccanismo di attenzione cattura le interazioni globali tra gli elementi della metasuperficie, consentendo una formazione efficace del fascio in campo prossimale. I pesi di attenzione vengono calcolati per quantificare la rilevanza delle coppie di nodi, assegnando un'importanza maggiore alle interazioni più informative. Utilizzando questi pesi, ogni nodo aggrega informazioni da tutti gli altri nodi, consentendo al modello di catturare dipendenze globali al di là dei vicinati locali. La fusione ponderata risultante delle caratteristiche genera una rappresentazione ottimizzata delle caratteristiche che enfatizza le caratteristiche spaziali importanti, sopprimendo al contempo le informazioni meno rilevanti, sostenendo così un'ottimizzazione efficace della formazione del fascio e della focalizzazione dell'energia.
Ottimizzazione del beamforming
La fase di ottimizzazione del beamforming ha determinato gli scostamenti di fase XL-IMS e i pesi del beamforming per focalizzare l'energia sugli utenti in campo prossimale. A differenza del classico puntamento del fascio in campo lontano, il beamforming in campo prossimale concentra l'energia in specifiche posizioni spaziali. Nel framework proposto, un modello con attenzione basato su GNN ha appreso direttamente i parametri ottimali di beamforming a partire dalle caratteristiche del canale, riducendo la complessità dei metodi di ottimizzazione iterativa. Il segnale ricevuto dall'utente è espresso nell'Equazione (4).
I dati di addestramento sono stati ottenuti da DeepMIMO per ogni realizzazione del canale, e gli spostamenti di fase XL-IMS corrispondenti sono stati determinati utilizzando l'algoritmo AO. La matrice del canale in campo prossimale h e la matrice degli spostamenti di fase XL-IMS sono state utilizzate per calcolare il vettore di spostamento di fase θ massimizzando la potenza del segnale ricevuto:
(8)
(9)
L'ottimizzatore AO è stato eseguito per 100 iterazioni fino al raggiungimento della soglia di convergenza desiderata di 10−4. Questi valori ottimizzati sono stati utilizzati come etichette di addestramento per il nostro modello GNN con attenzione. L'obiettivo era ottenere un'interferenza costruttiva per una focalizzazione efficace del fascio in campo prossimo e un potenziamento della potenza del segnale ricevuto. Invece di ricorrere a metodi di ottimizzazione convenzionali come SDR o AO, il problema è stato risolto mediante un approccio basato sull'apprendimento. In particolare, il problema di ottimizzazione della formazione del fascio è stato approssimato da un modello di apprendimento profondo fΘ(⋅).
È stato adottato un apprendimento supervisionato per addestrare il modello GNN con attenzione, utilizzando come obiettivi di addestramento i vettori di sfasamento ottimali per il beamforming forniti dall'algoritmo AO. In particolare, indichiamo il vettore di sfasamento ottimale generato dall'AO come θAO,i e le previsioni del modello come
. Quindi, il criterio di addestramento è stato formulato come segue:
(10)
dove N è il numero di campioni di addestramento.
Algoritmo 1: La formazione del fascio in campo prossimo mediante modello GNN con meccanismo di attenzione è presentata nel File Supplementare 1. In questo algoritmo, il primo passo consiste nell'impostare i parametri del sistema e il database DeepMIMO per creare condizioni di canale realistiche. La caratteristica del canale in campo prossimo per ogni campione è calcolata utilizzando il modello d'onda sferica, e quindi vengono estratte dal database le caratteristiche rilevanti, come i coefficienti del canale e le caratteristiche spaziali. Le caratteristiche estratte vengono convertite in una forma strutturata e successivamente rappresentate come un grafo, con gli elementi della metasuperficie come nodi. Il modello di apprendimento profondo è generato mediante strati GNN per catturare le dipendenze spaziali locali, e viene introdotto un meccanismo di attenzione per apprendere le dipendenze globali. Il framework di apprendimento profondo viene addestrato mediante apprendimento supervisionato, utilizzando la perdita MSE tra gli spostamenti di fase previsti e le etichette di spostamento di fase ottimali generate tramite AO.
Questa implementazione è stata eseguita con Python e PyTorch. I campioni per il canale DeepMIMO sono stati ottenuti dallo scenario esterno O1. Questi coefficienti di canale sono stati convertiti in grafi ed elaborati utilizzando una GNN con tre strati e un meccanismo di attenzione multitestina. Il modello è stato addestrato utilizzando l'ottimizzatore Adam con una funzione di perdita MSE. Questo modello ha generato gli scostamenti di fase e i pesi per la formazione del fascio in prossimità di utenti in campo vicino, e i risultati sono stati valutati in termini di velocità, efficienza e SINR.