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 degli approcci convenzionali in campo lontano. Viene sviluppato e convalidato un modello del canale in campo prossimale che incorpora la propagazione delle onde sferiche e la non stazionarietà spaziale, utilizzando il dataset DeepMIMO. Le informazioni sul canale sono rappresentate 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 del beamforming necessari per il focalizzazione del fascio in campo prossimale. A differenza dei metodi convenzionali 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 sono 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 configurata inizialmente tramite l'allestimento 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 di XL-IMS e le posizioni degli utenti
Progettare l'XL-IMS come una metasuperficie planare 16×16 che comprende 256 elementi riflettenti reconfigurabili. I dispositivi terminali utente sono stati distribuiti in modo casuale all'interno del 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 in vista diretta (LoS) che quelli non in vista diretta (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 forniti alla 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ò aiuta il modello a concentrarsi sulle 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 con attenzione GNN. Le complesse, laboriose procedure iterative di beamforming sono 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 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 composto 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, richiedendo 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, α si riferisce al parametro di attenuazione del segnale, λ si riferisce alla 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 ed è esprimibile come βn = dn−α, dove α rappresenta l'esponente di attenuazione del percorso. Va notato che dn è la distanza euclidea tra la posizione dell'utente e l'nesimo elemento della metasuperficie. Le variazioni del canale rispetto allo spazio sono causate dalle significative differenze in dn tra diverse metasuperfici, a differenza della situazione riscontrata nei canali in campo lontano. In generale, il vettore del canale può essere spiegato 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. Quindi, il segnale in ingresso può essere espresso come:
(4)
dove Φ = diag(ejθ1, ….., ejθN) è il segnale trasmesso, w è il vettore di precodifica a base di codifica e n rappresenta il rumore. Il vettore del canale in campo prossimale è 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 che viene utilizzato per precodificare il simbolo trasmesso (s). Durante il confronto delle prestazioni, approcci comuni di precodifica digitale come ZF e MMSE vengono implementati presso la stazione base 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 basato su GNN.
Per concentrare l'energia nella regione del campo vicino, una progettazione di θ in cui si verifica un'interferenza costruttiva, per l'utente, deve essere realizzabile. La strategia di formazione del fascio mira a massimizzare l'efficienza spettrale o la potenza del segnale ricevuto. La Figura 2 descrive l'architettura del modello di sistema di comunicazione in campo vicino basato sull'approccio XL-IMS.
Modellizzazione del canale in campo prossimale
Il modello proposto per reti di comunicazione 6G basate su XL-IMS nella regione prossimale è stato progettato specificamente per imitare le condizioni di propagazione vicine alla distanza di Rayleigh dalla metasuperficie. Nella regione prossimale, la propagazione è sferica, con valori di fase variabili e parametri spaziali non omogenei lungo la metasuperficie, a differenza dei modelli convenzionali di propagazione in campo lontano che assumono onde piane. Nello studio attuale, è stato utilizzato un modello fisico di onda sferica per progettare una matrice precisa della risposta all'impulso del canale, 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 l'n-esimo componente XL-IMS è data da:
(5)
Questo valore varia notevolmente attraverso la metasuperficie a causa della sua area estremamente grande, causando 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 a onda sferica descritto nell'Equazione (2), che tiene conto delle variazioni di fase e 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 la propagazione multipath sia in condizioni di visibilità diretta (LoS) che non in visibilità diretta (NLoS), basandosi 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 meta-superficie, 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, consentendo un'ottimizzazione del fascio dipendente dalla posizione. Questa informazione spaziale è incorporata nella matrice di risposta del canale ed è catturata dal modello di attenzione GNN. 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 dataset
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 a seconda delle 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 insiemi di addestramento e di test in un rapporto 80:20.
Preparazione del set di dati per l'addestramento
L'ultimo passaggio ha comportato la normalizzazione del dataset preparato e la sua conversione in tensori adatti all'input 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 con valori complessi in un formato adatto al deep learning. 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. Supplementary Figure 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 consentire un utilizzo 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 di ciascun elemento della metasuperficie dal ricevitore, insieme alle coordinate del ricevitore (x, y, z). In questo modo, il modello ha tenuto conto del fatto che le caratteristiche del segnale possono dipendere dalla posizione spaziale della coppia trasmettitore-ricevitore, una proprietà importante nella propagazione in campo prossimo.
Normalizzazione e preparazione dell'input del modello
Come passo preliminare prima di fornire al modello i suoi input, è stato eseguito un processo di normalizzazione per evitare instabilità numerica e accelerare la convergenza. I vettori strutturati sono stati suddivisi in insiemi di dati di addestramento e di test, indicati con (Xtrain, Xtest). In questo modo, gli input sono stati successivamente forniti ai modelli GNN e di attenzione.
Sulla base della rappresentazione delle caratteristiche mostrata nella Figura Suppletiva 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 forniscono una descrizione numerica del segnale, le parti in modulo e in fase forniscono rispettivamente una descrizione esplicita dell'intensità del segnale e dello sfasamento. 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 in campo prossimo, 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 prossimo, una rete completamente connessa non è adatta a rappresentare correttamente il canale. Pertanto, in questo approccio, la metasuperficie è rappresentata come un grafo, in cui ogni elemento funge da nodo e le interazioni tra di 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 gli archi 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 la connettività tra i vertici che sono vicini tra loro o entro determinati vincoli 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 su 8. Lasciare pi e pj rappresentano i vettori di posizione dei i-e e j-celle XL-IMS, rispettivamente. La distanza tra due vertici qualsiasi è misurata mediante la metrica euclidea come segue: dij = ||pi − pj||2. sarà presente un arco tra i vertici i e j se vertice j è incluso nel gruppo degli otto vicini più prossimi al vertice iPertanto, la matrice di adiacenza 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 delle caratteristiche del vertice i nel livello l, N(i) è il vicinato del vertice i, W(l) è la matrice dei pesi apprendibile e σ(·) simboleggia una funzione di attivazione non lineare.
Il modello GNN con attenzione era composto da tre strati di convoluzione su grafi con dimensioni nascoste impostate a 128, e come funzione di attivazione è stata utilizzata l'unità lineare retificata (ReLU). La generazione del grafo è stata eseguita mediante il metodo dei k-vicini più prossimi (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 esistenti al di fuori del raggio di vicinato locale, è stata introdotta l'attenzione nel modello, che funziona come segue:
(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 la comprensione del contesto globale attraverso la metasuperficie. Il framework ibrido cattura informazioni spaziali sia locali che globali, e le caratteristiche apprese vengono elaborate attraverso strati completamente connessi per generare i pesi della formazione del fascio. Come mostrato nelle Figure Supplettive 3 e 4, la GNN modella le dipendenze spaziali locali attraverso 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 pesata risultante delle caratteristiche genera una rappresentazione ottimizzata delle caratteristiche che enfatizza le caratteristiche spaziali importanti mentre sopprime le informazioni meno rilevanti, supportando un'ottimizzazione efficace della formazione del fascio e il concentramento dell'energia.
Ottimizzazione del beamforming
La fase di ottimizzazione del beamforming ha determinato gli sfasamenti 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 GNN con attenzione 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 corrispondenti di XL-IMS sono stati determinati utilizzando l'algoritmo AO. La matrice del canale in campo prossimale h e la matrice degli spostamenti di fase di 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 nel campo prossimale e un aumento 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 a fasatura ottimale forniti dall'algoritmo AO. In particolare, indichiamo con θAO,i il vettore di sfasamento ottimale generato da AO e con
le previsioni del modello. Il criterio di addestramento è stato quindi formulato come segue:
(10)
dove N è il numero di campioni di addestramento.
Algoritmo 1: La formazione del fascio in campo prossimale 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 prossimale per ogni campione è stata calcolata utilizzando il modello d'onda sferica, e quindi sono state estratte dal database le caratteristiche rilevanti, come i coefficienti del canale e le caratteristiche spaziali. Le caratteristiche estratte sono state convertite in una forma strutturata e successivamente rappresentate come un grafo, con gli elementi della metasuperficie come nodi. Il modello di apprendimento profondo è stato sviluppato utilizzando strati GNN per catturare le dipendenze spaziali locali, ed è stato introdotto un meccanismo di attenzione per apprendere le dipendenze globali. Il framework di apprendimento profondo è stato addestrato mediante apprendimento supervisionato, utilizzando la perdita MSE tra gli spostamenti di fase predetti e le etichette di spostamento di fase ottimali generate dall'AO.
Questa implementazione è stata eseguita utilizzando Python e PyTorch. I campioni per il canale DeepMIMO sono stati ottenuti dallo scenario esterno O1. Questi coefficienti del canale sono stati convertiti in grafi ed elaborati mediante una GNN con tre strati e un meccanismo di attenzione multi-testa. Il modello è stato addestrato utilizzando l'ottimizzatore Adam con una funzione di perdita MSE. Questo modello ha generato spostamenti di fase e pesi per la formazione del fascio in utenti in campo prossimale, e i risultati sono stati valutati in termini di velocità, efficienza e SINR.