Le cadre proposé pour la modélisation des canaux en champ proche et le beamforming dans les systèmes 6G dotés d'IMS étendus (XL-IMS) combine une modélisation basée sur la physique avec une optimisation par apprentissage automatique afin de surmonter les limites des approches conventionnelles en champ lointain. Un modèle de canal en champ proche intégrant la propagation d'ondes sphériques et la non-stationnarité spatiale est développé et validé à l'aide du jeu de données DeepMIMO. L'information sur le canal est représentée sous forme de graphe afin de capturer les dépendances spatiales entre les éléments de la surface méta et les utilisateurs, et un réseau d'attention GNN de bout en bout est utilisé pour apprendre les poids de beamforming nécessaires à la focalisation du faisceau en champ proche. Contrairement aux méthodes d'optimisation itératives classiques, ce cadre fondé sur l'apprentissage prédit directement les solutions de beamforming, réduisant ainsi la complexité computationnelle tout en conservant une grande précision et une bonne adaptabilité aux conditions dynamiques du canal. Les performances sont évaluées en termes de débit atteignable, de gain de beamforming et d'efficacité énergétique, montrant des résultats supérieurs par rapport aux méthodes traditionnelles de beamforming en champ lointain. Figure 1 illustre le cadre global.
Les étapes proposées sont les suivantes :
Étape 1 : Configurer l'environnement DeepMIMO
La plateforme DeepMIMO v3.x est d'abord configurée à l'aide du paramétrage extérieur O1, qui utilise une fréquence porteuse de 28 GHz et une bande passante de 100 MHz. Dans la simulation, l'XL-IMS est représenté par un réseau plan uniforme de 256 éléments réfléchissants disposés en configuration 16×16, espacés de λ/2. Le modèle de propagation choisi est une onde sphérique comprenant des trajets en visibilité directe (LoS) et en visibilité non directe (NLoS).
Étape 2 : Définir la géométrie XL-IMS et les emplacements des utilisateurs
Concevoir le XL-IMS comme une métasurface plane de 16×16 comprenant 256 éléments réfléchissants reconfigurables. Les dispositifs terminaux ont été répartis aléatoirement dans la zone de champ proche, entre 0,5 m et 10 m du XL-IMS, selon le critère de la distance de Rayleigh.
Étape 3 : Générer des canaux en champ proche en utilisant la propagation d'ondes sphériques
Modélisez les canaux entre chaque élément de la métasurface reconfigurable et l'utilisateur en utilisant la propagation d'ondes sphériques, qui tient compte de la variation de phase induite par la distance.
Étape 4 : Générer des données de canal en champ proche
Les terminaux utilisateur étaient répartis de manière uniforme dans la zone en champ proche (0,5-10 m) par rapport au panneau XL-IMS. Les trajets de propagation en vue directe (LoS) et en vue non directe (NLoS), générés par le scénario extérieur DeepMIMO O1, ont été utilisés pour générer les canaux.
Étape 5 : Préparer le jeu de données pour l'entraînement et les tests
L'ensemble de données généré comprenait 20 000 échantillons, dont 16 000 (80 %) ont été utilisés pour l'apprentissage et 4 000 (20 %) pour les tests.
Étape 6 : Normaliser les caractéristiques d'entrée et préparer les entrées du modèle
La dernière étape du processus de préparation consistait à appliquer une normalisation Min-Max à toutes les caractéristiques d'entrée, telles que l'amplitude du canal, la phase, la distance et les coordonnées de l'utilisateur. Ce processus a permis de normaliser les données et de les convertir en tenseurs, qui ont ensuite été injectés dans le réseau GNN avec mécanisme d'attention. Le modèle a été optimisé à l'aide de l'optimiseur Adam avec un taux d'apprentissage de 0,001 et sur 50 époques.
Étape 7 : Appliquer le mécanisme d'attention
Ajouter un mécanisme d'attention afin de prendre en compte les informations globales du graphe et de capturer les dépendances à longue portée entre les nœuds du graphe. Cela aide le modèle à se concentrer sur les caractéristiques spatiales importantes.
Étape 8 : Prédire les poids de formation de faisceau
Prédire les poids de formation de faisceau et les valeurs de décalage de phase à l'aide d'un embedding basé sur l'attention des GNN. Les procédures complexes, fastidieuses et itératives de formation de faisceau sont remplacées par une méthode de prédiction.
Étape 9 : Évaluer les métriques de performance
Pour évaluer l'efficacité de la communication et du beamforming, calculez des indicateurs clés de performance tels que le débit réalisable, le gain de beamforming, l'efficacité spectrale, l'efficacité énergétique et le SINR.
Étape 10 : Comparer aux méthodes de référence
Comparer l'algorithme proposé aux algorithmes traditionnels, notamment ZF, MMSE, SDR et AO.
Définition du modèle système
Le système proposé considère un réseau sans fil 6G rendu possible par un XL-IMS composé de 256 éléments réfléchissants reconfigurables qui contrôlent la propagation des ondes électromagnétiques grâce à des décalages de phase ajustables. Étant donné que les utilisateurs fonctionnent dans la région du champ proche, l'hypothèse classique d'onde plane en champ lointain devient invalide, ce qui impose de prendre en compte la propagation d'ondes sphériques, l'atténuation dépendante de la distance et la non-stationnarité spatiale. La limite entre les régions du champ proche et du champ lointain est déterminée par la distance de Rayleigh, qui dépend de la taille de la métasurface et de la longueur d'onde porteuse, et marque la transition entre les effets de front d'onde sphérique et la propagation d'onde plane. Ceci est décrit comme suit :
(1)
Ici, D est la taille de la métasurface, λ est la longueur d'onde, et R est la distance de Rayleigh. Un phénomène de champ proche se produit dans le canal de communication lorsqu'un utilisateur se trouve à une position plus proche de R que sa valeur réelle. Lorsque le signal traverse différents composants de la métasurface, il subit des variations distinctes de phase et d'amplitude.
Le canal de la n-ième méta-surface vers l'utilisateur utilise un modèle d'onde sphérique pour la propagation en champ proche, représenté par l'équation du coefficient de canal suivante :
(2)
Ici, α désigne le paramètre de perte de parcours, λ désigne la longueur d'onde porteuse, et dn est la distance de l'utilisateur par rapport au composant.
Dans ce cas, βn est le facteur de perte de trajet pour l'élément XL-IMS nème et peut s'exprimer par βn = dn−α, où α représente l'exposant de perte de trajet. Il convient de noter que dn est la distance euclidienne entre la position de l'utilisateur et l'élément nème de la métasurface. Les variations du canal en fonction de l'espace sont dues aux différences significatives de dn entre les différentes métasurfaces, contrairement à la situation rencontrée avec les canaux en champ lointain. En général, le vecteur de canal peut être expliqué comme suit :
(3)
Le signal d'entrée reçu par l'utilisateur dépend du réglage de phase du système XL-IMS. Considérons θ = [θ1, θ2, …., θN] pour représenter le déphasage réalisé à l'aide de la métasurface. Le signal d'entrée peut alors être exprimé comme suit :
(4)
où Φ = diag(ejθ1, ….., ejθN) est le signal transmis, w est le vecteur de précodage de base, et n représente le bruit. Le vecteur de canal en champ proche est défini comme h ∈ ℂN×1, la matrice diagonale de décalage de phase du XL-IMS comme Φ ∈ ℂN×N, le vecteur de formation de faisceau comme x ∈ ℂN×1, et le signal reçu comme y ∈ ℂ. Le terme de bruit additif est défini comme n ∼ CN (0, σ2). Ces dimensions préservent la rigueur mathématique de la formulation.
Dans ce cas, ws est le vecteur de précodage de la station de base utilisé pour pré-coder le symbole transmis (s). Lors de l'évaluation des performances, des approches courantes de précodage numérique telles que ZF et MMSE sont mises en œuvre au niveau de la station de base à l'aide de réalisations similaires du canal en champ proche. Toutefois, la nouvelle approche vise à optimiser la matrice de déphasage XL-IMS (θ) grâce au mécanisme d'apprentissage par attention GNN.
Afin de concentrer l'énergie dans la région du champ proche, un angle de conception θ permettant une interférence constructive pour l'utilisateur doit être réalisable. La stratégie de formation de faisceau vise à maximiser l'efficacité spectrale ou la puissance du signal reçu. La figure 2 décrit l'architecture du modèle de système de communication en champ proche basé sur l'approche XL-IMS.
Modélisation des canaux en champ proche
Le modèle de réseau de communication 6G basé sur XL-IMS proposé a été spécifiquement conçu pour imiter les conditions de propagation proches de la distance de Rayleigh par rapport à la métasurface. Dans la région du champ proche, la propagation est sphérique, avec des valeurs de phase variables et des paramètres spatiaux non homogènes à travers la métasurface, contrairement aux modèles conventionnels de propagation en champ lointain qui supposent des ondes planes. Dans la présente étude, un modèle d'onde sphérique fondé sur la physique a été utilisé pour concevoir une matrice de réponse impulsionnelle de canal précise, tenant compte des positions physiques des éléments de la métasurface ainsi que des terminaux utilisateur.
Représentation des canaux basée sur la distance
La distance entre chaque composant de la métasurface et l'utilisateur est un facteur crucial dans l'élaboration du modèle de canal en régime proche. La longueur pn entre les positions de l'utilisateur pu et le n-ième composant du XL-IMS est donnée par :
(5)
Cette valeur varie considérablement sur l'ensemble de la métasurface en raison de sa surface extrêmement grande, ce qui entraîne des temps de propagation et des décalages de phase inégaux qui doivent être pris en compte dans le modèle.
Modèle de canal à onde sphérique
Le coefficient de canal entre l'élément n-ième et l'utilisateur est déterminé à l'aide du modèle d'onde sphérique décrit dans l'équation (2), qui prend en compte les variations de phase et d'amplitude dépendant de la distance. Pour la génération du jeu de données et l'entraînement du modèle, les coefficients de canal sont obtenus à partir du cadre de travail DeepMIMO basé sur le tracé de rayons, qui intègre à la fois la propagation multipath en vue directe (LoS) et hors vue directe (NLoS), selon des structures environnementales réalistes. À partir de ce qui précède, il est possible de construire la matrice de réponse de canal à partir des réponses de canal pour les N éléments de la surface métamatérielle, comme indiqué dans l'équation (3).
Non-stationnarité spatiale des canaux en champ proche et importance du modèle
Contrairement aux canaux en champ lointain, les canaux en champ proche présentent une non-stationnarité spatiale, où les signaux reçus varient à travers la surface métamatérielle, permettant une focalisation de faisceau dépendante de la position. Cette information spatiale est intégrée dans la matrice de réponse du canal et est capturée par le modèle d'attention GNN. Figure supplémentaire 1 illustre le cadre proposé de modélisation du canal en champ proche pour le système de communication 6G basé sur XL-IMS et la génération de données de canal réalistes pour le beamforming basé sur l'apprentissage profond.
Génération du jeu de données
La génération des données a été réalisée à l'aide du cadre DeepMIMO, qui fournit des simulations de canal réalistes basées sur le tracé de rayons, intégrant la géométrie de l'environnement, la cohérence spatiale et la mobilité des utilisateurs. Les réponses de canal pour plusieurs emplacements d'utilisateurs en champ proche ont été générées à l'aide du scénario extérieur O1, à une fréquence porteuse de 28 GHz et une largeur de bande de 100 MHz. Le XL-IMS était composé de 256 éléments réfléchissants uniformément répartis, et les emplacements des utilisateurs ont été choisis dans la région de champ proche selon le critère de Rayleigh. Les coefficients de canal, les gains de trajet et les informations spatiales ont été extraits à l'aide de la bibliothèque Python DeepMIMO. Tableau supplémentaire 1 résume les paramètres de simulation et de reproductibilité. Le modèle a été entraîné selon un apprentissage supervisé, les caractéristiques d'état de canal servant d'entrées et les poids de formation de faisceau optimisés selon des principes physiques servant de labels cibles.
Création des caractéristiques d'entrée et structure du jeu de données
Ensuite, les données brutes des canaux ont été converties en caractéristiques en vue de leur entrée dans des modèles d'apprentissage profond. Ces caractéristiques comprenaient les amplitudes et les phases des réponses des canaux, prises séparément. D'autres paramètres d'entrée pouvaient également être inclus selon les exigences du problème. Par exemple, les coordonnées des utilisateurs, les distances et certains paramètres spatiaux pouvaient également être considérés comme des caractéristiques. L'ensemble de données a été divisé en ensembles d'apprentissage et de test selon un ratio de 80:20.
Préparation du jeu de données pour l'apprentissage
La dernière étape consistait à normaliser le jeu de données préparé et à le convertir en tenseurs adaptés en entrée pour l'architecture GNN avec mécanisme d'attention lors de l'apprentissage. Le jeu de données d'apprentissage a été utilisé pour apprendre l'application de l'état du canal vers les poids optimaux de formation de faisceau, tandis que le jeu de données de test a permis d'évaluer l'efficacité du modèle. Ce jeu de données soigneusement préparé a facilité et rendu plus efficace le processus d'apprentissage. Supplementary Table 2 présente les détails de la description du jeu de données DeepMIMO.
Représentation des caractéristiques
L'étape de représentation des caractéristiques a transformé les données de canal en champ proche à valeurs complexes en un format adapté à l'apprentissage profond. Les réponses du canal ont été décomposées en caractéristiques capturant les propriétés spatiales et de propagation, tout en préservant les informations essentielles du champ proche pour l'apprentissage. Figure supplémentaire 2 illustre le jeu de données DeepMIMO utilisé dans ce processus. Les données de canal provenant de l'environnement en champ proche ont été représentées sous forme de nombres complexes hn, avec des composantes réelle et imaginaire qui représentent l'information de phase et d'amplitude des coefficients du canal. Afin d'être efficacement utilisées dans des algorithmes d'apprentissage automatique, les coefficients ont été décomposés en leurs parties réelle (Rehn) et imaginaire (Imhn). Les données peuvent également être représentées par l'amplitude hn et la phase des coefficients ∠hn.
Extraction de caractéristiques spatiales
En plus du coefficient de canal, des caractéristiques spatiales ont été incluses afin de tenir compte de la nature géométrique de l'environnement en champ proche. Cela comprenait la distance dn entre chaque élément de la métasurface et le récepteur, ainsi que les coordonnées du récepteur (x, y, z). De cette manière, le modèle prenait en compte le fait que les caractéristiques du signal peuvent dépendre de la position spatiale de la paire émetteur-récepteur, une propriété importante pour la propagation en champ proche.
Normalisation et préparation des entrées du modèle
Comme étape préliminaire avant de fournir les entrées au modèle, un processus de normalisation a été effectué afin d'éviter l'instabilité numérique et d'accélérer la convergence. Les vecteurs structurés ont été divisés en ensembles de données d'entraînement et de test, notés (Xtrain, Xtest). De cette manière, les entrées ont ensuite été fournies aux modèles GNN et aux modèles d'attention.
Sur la base de la représentation des caractéristiques présentée dans la Figure supplémentaire 2, différents types de caractéristiques ont été utilisés dans ce jeu de données, notamment les coefficients de canal hn, qui indiquaient la réponse complexe du signal de tous les éléments de la métasurface. Une autre caractéristique était la distance dn, qui, dans ce cas, indiquait les effets de propagation. Bien que les composantes imaginaire et réelle fournissent une description numérique du signal, les parties en amplitude et en phase décrivent explicitement l'intensité du signal et le déphasage, respectivement. La coordonnée de position de l'utilisateur (x, y, z) constituait une autre caractéristique importante qui améliorait la perception spatiale du modèle.
Conception par apprentissage profond utilisant un GNN pour la modélisation des dépendances spatiales, intégrée à un mécanisme d'attention hybride pour l'apprentissage du contexte global
Cet algorithme d'apprentissage a été conçu pour entraîner efficacement le canal XL-IMS en champ proche complexe, avec ses caractéristiques spatiales et de propagation, en utilisant des réseaux neuronaux graphiques (GNN) et un mécanisme d'attention. Étant donné que les éléments de la métasurface sont fortement corrélés entre eux et que leurs interactions sont non uniformes dans la zone de champ proche, un réseau entièrement connecté n'est pas adapté pour représenter correctement le canal. Ainsi, dans cette approche, la métasurface est représentée sous forme de graphe, chaque élément jouant le rôle d'un nœud et les interactions entre eux constituant les arêtes. Le GNN a capturé les dépendances spatiales localisées entre les nœuds voisins, tandis que le mécanisme d'attention a amélioré l'apprentissage des dépendances globales.
Représentation graphique de la SM-IPL
Le XL-IMS est créé sous la forme d'un graphe G = (V, E), chaque sommet vn ∈ V représentant un élément de la méta-surface, tandis que les arêtes eij ∈ E représentent les relations spatiales entre les éléments. Chaque sommet possède un vecteur de caractéristiques xn, obtenu lors de la phase de représentation des caractéristiques. La matrice d'adjacence A décrit les connexions entre les sommets proches ou situés dans certaines limites de distance. Une telle formulation sous forme de graphe permet d'intégrer l'information spatiale présente dans la propagation en champ proche.
Pour assurer la reproductibilité, la connectivité du graphe est déterminée à l'aide de l'approche des K-plus proches voisins avec k fixé à 8. Soient pi et pj les vecteurs de position respectifs des cellules XL-IMS de rang i et j. La distance entre deux sommets quelconques est mesurée à l'aide de la métrique euclidienne selon la formule suivante : dij = ||pi − pj||2. Une arête existera entre les sommets i et j si le sommet j appartient au groupe des huit voisins les plus proches du sommet i. Par conséquent, la matrice d'adjacence est construite selon la règle suivante : Aij = 1 si j ∈ Nk(i) et Aij = 0 sinon.
Dans le GNN proposé, les mises à jour de chaque nœud sont calculées à partir des informations provenant des autres nœuds voisins. En particulier, la propagation de la couche l peut être exprimée comme suit :
(6)
où hi(l) représente le vecteur de caractéristiques du sommet i dans la couche l, N(i) est le voisinage du sommet i, W(l) est la matrice de poids apprenable, et σ(·) symbolise une fonction d'activation non linéaire.
Le modèle d'attention GNN comprenait trois couches de convolution sur graphes avec des dimensions cachées fixées à 128, et l'unité rectifiée linéaire (ReLU) était utilisée comme fonction d'activation. La génération du graphe a été réalisée à l'aide de la méthode des k-plus proches voisins (k = 8). Chaque caractéristique XL-IMS agit comme un nœud, tandis que les arêtes sont formées en fonction de la distance euclidienne entre les nœuds adjacents. Pour extraire les dépendances spatiales globales, nous avons utilisé un réseau d'auto-attention multi-têtes avec 4 têtes. La probabilité de dropout a été fixée à 0,3 afin d'éviter le surapprentissage.
Mécanisme d'attention
Pour modéliser les dépendances à longue distance existant au-delà de la portée du voisinage local, l'attention a été introduite dans le modèle, qui fonctionne comme suit :
(7)
où eij = LeakyReLU(aT[Whi | Whj]) représente l'importance des nœuds i et j. Grâce à ce modèle, il est possible d'attribuer différents niveaux d'importance aux nœuds en fonction de leurs rôles, facilitant ainsi la compréhension du contexte global à travers la métasurface. Le cadre hybride capture à la fois les informations spatiales locales et globales, et les caractéristiques apprises sont traitées par des couches entièrement connectées afin de générer les poids de formation de faisceau. Comme illustré dans les figures supplémentaires 3 et 4, le réseau de neurones graphique (GNN) modélise les dépendances spatiales locales par le biais de la transmission de messages et de l'agrégation de caractéristiques, tandis que le mécanisme d'attention capture les interactions globales entre les éléments de la métasurface, permettant une formation efficace de faisceau en champ proche. Des poids d'attention sont calculés pour quantifier la pertinence des paires de nœuds, en attribuant une importance plus grande aux interactions les plus informatives. À l'aide de ces poids, chaque nœud agrège l'information provenant de tous les autres nœuds, permettant au modèle de capturer des dépendances globales allant au-delà des voisinages locaux. La fusion pondérée des caractéristiques qui en résulte produit une représentation optimisée des caractéristiques, mettant en évidence les caractéristiques spatiales importantes tout en supprimant les informations moins pertinentes, ce qui soutient une optimisation efficace de la formation de faisceau et un focalisation précise de l'énergie.
Optimisation du beamforming
L'étape d'optimisation du beamforming a déterminé les décalages de phase XL-IMS et les poids de beamforming afin de focaliser l'énergie sur les utilisateurs en champ proche. Contrairement au pointage de faisceau classique en champ lointain, le beamforming en champ proche concentre l'énergie à des emplacements spatiaux spécifiques. Dans le cadre proposé, un modèle à attention basé sur un GNN a appris directement les paramètres optimaux de beamforming à partir des caractéristiques du canal, réduisant ainsi la complexité des méthodes d'optimisation itératives. Le signal reçu par l'utilisateur est donné par l'équation (4).
Les données d'entraînement ont été obtenues à partir de DeepMIMO pour chaque réalisation de canal, et les décalages de phase XL-IMS correspondants ont été déterminés à l'aide de l'algorithme AO. La matrice de canal en champ proche h et la matrice de décalage de phase XL-IMS ont été utilisées pour calculer le vecteur de décalage de phase θ en maximisant la puissance du signal reçu :
(8)
(9)
L'optimiseur AO a été exécuté pendant 100 itérations jusqu'à ce qu'il atteigne le seuil de convergence souhaité de 10−4. Ces valeurs optimisées ont servi de labels d'apprentissage pour notre modèle GNN avec attention. L'objectif était d'obtenir une interférence constructive afin de concentrer efficacement le faisceau en champ proche et d'améliorer la puissance du signal reçu. Plutôt que d'utiliser des méthodes d'optimisation classiques telles que SDR ou AO, le problème a été résolu par une approche fondée sur l'apprentissage. Plus précisément, le problème d'optimisation du beamforming a été approché par un modèle de deep learning fΘ(⋅).
Un apprentissage supervisé a été utilisé pour entraîner le modèle GNN avec attention, les vecteurs de déphasage optimaux issus de l'algorithme AO servant de cibles d'entraînement. Plus précisément, notons le vecteur de déphasage optimal généré par AO comme étant θAO,i et les prédictions du réseau comme
. Le critère d'entraînement a alors été formulé comme suit :
(10)
où N est le nombre d'échantillons d'apprentissage.
Algorithme 1 : La formation de faisceau en champ proche à l’aide d’un modèle avec attention basée sur les GNN est présentée dans le Fichier Supplémentaire 1. Dans cet algorithme, la première étape consiste à configurer les paramètres du système et la base de données DeepMIMO afin de créer des conditions de canal réalistes. La caractéristique de canal en champ proche pour chaque échantillon est calculée à l’aide du modèle d’onde sphérique, puis les caractéristiques importantes, telles que les coefficients de canal et les caractéristiques spatiales, sont extraites de la base de données. Les caractéristiques extraites sont converties sous une forme structurée, puis représentées sous forme de graphe, les éléments de la métasurface constituant les nœuds. Le modèle d’apprentissage profond est construit à l’aide de couches GNN afin de capturer les dépendances spatiales locales, et un mécanisme d’attention est introduit pour apprendre les dépendances globales. Le cadre d’apprentissage profond est entraîné selon une approche supervisée, en utilisant une fonction de perte MSE entre les décalages de phase prédits et les étiquettes de décalages de phase optimaux générées par AO.
Cette implémentation a été réalisée avec Python et PyTorch. Les échantillons pour le canal DeepMIMO ont été obtenus à partir du scénario extérieur O1. Ces coefficients de canal ont été convertis en graphes et traités à l'aide d'un GNN comportant trois couches et un mécanisme d'attention multi-têtes. Le modèle a été entraîné à l'aide de l'optimiseur Adam avec une fonction de perte MSE. Ce modèle a généré des décalages de phase et des poids pour le beamforming auprès d'utilisateurs en champ proche, et les résultats ont été évalués en termes de débit, d'efficacité et de SINR.