Proponowany schemat modelowania kanałów w polu bliskim i formowania wiązki w systemach 6G obsługujących XL-IMS łączy modelowanie oparte na fizyce z optymalizacją za pomocą uczenia maszynowego, aby przezwyciężyć ograniczenia konwencjonalnych podejść dla pola dalekiego. Opracowano i zwalidowano za pomocą zbioru danych DeepMIMO model kanału w polu bliskim, uwzględniający propagację fal sferycznych oraz przestrzenną niestacjonarność. Informacje o kanale są reprezentowane jako graf, aby uchwycić zależności przestrzenne między elementami metasuperficyjnymi a użytkownikami, a do wyznaczania wag formowania wiązki dla ogniskowania wiązki w polu bliskim wykorzystano end-to-end sieć GNN-attention. W przeciwieństwie do konwencjonalnych iteracyjnych metod optymalizacji, schemat oparty na uczeniu bezpośrednio przewiduje rozwiązania formowania wiązki, co zmniejsza złożoność obliczeniową przy zachowaniu wysokiej dokładności i adaptacyjności do dynamicznych warunków kanału. Wydajność oceniono pod kątem osiągalnej przepływności, wzmocnienia formowania wiązki i efektywności energetycznej, wykazując wyniki lepsze w porównaniu z tradycyjnymi metodami formowania wiązki w polu dalekim. Rysunek 1 przedstawia ogólny schemat.
Poniżej przedstawiono proponowane kroki:
Krok 1: Konfiguracja środowiska DeepMIMO
Platforma DeepMIMO v3.x jest w pierwszej kolejności konfigurowana za pomocą ustawień O1 Outdoor, w których zastosowano częstotliwość przenoszenia 28 GHz oraz szerokość pasma 100 MHz. W symulacji XL-IMS jest reprezentowany przez jednorodną macierz planarna 16×16 składającą się z 256 elementów odbijających, rozmieszczonych w odstępach λ/2. Wybrany model propagacji to fala sferyczna z uwzględnieniem ścieżek LoS i NLoS.
Krok 2: Definiowanie geometrii XL-IMS i lokalizacji użytkowników
Zaprojektuj XL-IMS jako płaską metapowierzchnię o wymiarach 16×16, składającą się z 256 rekonfigurowalnych elementów refleksyjnych. Urządzenia terminalne użytkownika rozmieszczono losowo w obszarze pola bliskiego, od 0,5m do 10m od XL-IMS, w oparciu o kryterium odległości Rayleigha.
Krok 3: Generowanie kanałów pola bliskiego przy użyciu propagacji fal sferycznych
Modeluj kanały pomiędzy każdym elementem rekonfigurowalnej metasurface a użytkownikiem przy użyciu propagacji fal sferycznych, uwzględniając zmianę fazy wywołaną odległością.
Krok 4: Generowanie danych kanału bliskiego pola
Terminalne urządzenia użytkowników zostały rozmieszczone jednorodnie w strefie bliskiego pola (0,5-10 m) od panelu XL-IMS. Do wygenerowania kanałów wykorzystano zarówno ścieżki propagacji LoS, jak i NLoS, stworzone w ramach scenariusza DeepMIMO O1 Outdoor Scenario.
Krok 5: Przygotowanie zbioru danych do uczenia i testowania
Wygenerowany zestaw danych obejmował 20 000 próbek, z których 16 000 (80%) wykorzystano do uczenia, a 4 000 (20%) do testowania.
Krok 6: Normalizacja cech wejściowych i przygotowanie danych wejściowych modelu
Ostatni etap procesu przygotowania polegał na przeprowadzeniu skalowania Min-Max wszystkich cech wejściowych, takich jak amplituda kanału, faza, odległość oraz współrzędne użytkownika. Proces ten przeskalował dane i przekształcił je w tensory, które zostały wprowadzone do sieci GNN-attention. Model został zoptymalizowany przy użyciu optymalizatora Adam z szybkością uczenia 0,001 i 50 epokami.
Krok 7: Zastosowanie mechanizmu uwagi
Dodaj mechanizm uwagi, aby uwzględnić globalne informacje o grafie i uchwycić długodystansowe zależności między węzłami grafu. Pomaga to modelowi skupić się na istotnych cechach przestrzennych.
Krok 8: Przewidywanie wag formowania wiązki
Przewidź wagi formowania wiązki i wartości przesunięcia fazowego, wykorzystując osadzanie z mechanizmem uwagi GNN (GNN-attention embedding). Złożone, żmudne i iteracyjne procedury formowania wiązki zostają zastąpione metodą predykcyjną.
Krok 9: Ocena wskaźników wydajności
Aby ocenić skuteczność komunikacji i formowania wiązki, należy obliczyć wskaźniki KPI, takie jak osiągalna przepływność, zysk formowania wiązki, efektywność widmowa, efektywność energetyczna oraz SINR.
Krok 10: Porównanie z metodami referencyjnymi
Porównaj proponowany algorytm z algorytmami tradycyjnymi, w tym ZF, MMSE, SDR oraz AO.
Definicja modelu systemu
Proponowany system rozważa bezprzewodową sieć 6G wspieraną przez XL-IMS składający się z 256 rekonfigurowalnych elementów odbijających, które kontrolują propagację fal elektromagnetycznych poprzez regulowane przesunięcia fazowe. Ponieważ użytkownicy operują w obszarze pola bliskiego, konwencjonalne założenie o płaskiej fali pola dalekiego staje się nieprawidłowe, co wymaga uwzględnienia propagacji fal sferycznych, tłumienia zależnego od odległości oraz przestrzennej niestacjonarności. Granica między obszarami pola bliskiego i dalekiego jest określona przez odległość Rayleigha, która zależy od rozmiaru metasurfacy oraz długości fali nośnej i wyznacza przejście od efektów sferycznych frontów fal do propagacji fali płaskiej. Jest to opisane w następujący sposób:
(1)
Tutaj D oznacza rozmiar metasuperfacy, λ to długość fali, a R to odległość Rayleigha. Zjawisko pola bliskiego występuje w kanale komunikacyjnym, gdy użytkownik znajduje się w pozycji bliższej niż rzeczywista wartość R. W miarę przemieszczania się sygnału przez poszczególne komponenty metasuperfacy, ulega on wyraźnym zmianom fazy i amplitudy.
Kanał od n-tej metasurface do użytkownika wykorzystuje model fali sferycznej dla propagacji w polu bliskim, co jest reprezentowane przez równanie współczynnika kanału w następujący sposób:
(2)
W tym przypadku α oznacza parametr straty propagacyjnej, λ oznacza długość fali nośnej, a dn to odległość użytkownika od komponentu.
W tym przypadku βn jest współczynnikiem strat propagacyjnych dla n-tego elementu XL-IMS i może być wyrażony jako βn = dn−α, gdzie α oznacza wykładnik strat propagacyjnych. Należy zauważyć, że dn to odległość euklidesowa od pozycji użytkownika do n-tego elementu metapowierzchni. Zmiany kanału w odniesieniu do przestrzeni są powodowane znacznymi różnicami w dn pomiędzy różnymi metapowierzchniami, w przeciwieństwie do sytuacji w przypadku kanałów pola dalekiego. W ujęciu ogólnym wektor kanału można wyjaśnić w następujący sposób:
(3)
Sygnał wejściowy odbierany przez użytkownika zależy od strojenia fazowego systemu XL-IMS. Przyjmijmy, że θ = [θ1, θ2, …., θN] reprezentuje przesunięcie fazowe zrealizowane za pomocą metasurfacy. Wówczas sygnał wejściowy można zapisać jako:
(4)
gdzie Φ = diag(ejθ1, ….., ejθN) jest sygnałem nadawanym, w to wektor prekodowania kodowania bazowego, a n oznacza szum. Wektor kanału w polu bliskim jest zdefiniowany jako h ∈ ℂN×1, diagonalna macierz przesunięcia fazy XL-IMS jako Φ ∈ ℂN×N, wektor formowania wiązki jako x ∈ ℂN×1, a sygnał odebrany jako y ∈ ℂ. Składnik szumu addytywnego jest zdefiniowany jako n ∼ CN (0, σ2). Wymiary te zapewniają poprawność matematyczną sformułowania.
W tym przypadku ws jest wektorem prek either kodowania stacji bazowej, który jest wykorzystywany do prek either kodowania transmitowanego symbolu (s). Podczas benchmarkingu wydajności, w stacji bazowej implementowane są powszechne podejścia do cyfrowego prek either kodowania, takie jak ZF i MMSE, przy użyciu podobnych realizacji kanałów w polu bliskim. Jednak nowe podejście dąży do optymalizacji macierzy przesunięć fazowych XL-IMS (θ) poprzez mechanizm uczenia z wykorzystaniem sieci GNN-attention.
Aby skoncentrować energię w obszarze bliskiego pola, dla użytkownika powinno być możliwe osiągnięcie konfiguracji θ, w której występuje interferencja konstruktywna. Strategia formowania wiązki ma na celu maksymalizację efektywności widmowej lub mocy sygnału odbieranego. Rycina 2 przedstawia architekturę modelu systemu komunikacji w bliskim polu opartego na podejściu XL-IMS.
Modelowanie kanału w polu bliskim
Zaproponowany model pola bliskiego dla sieci komunikacyjnych 6G oparty na XL-IMS został specjalnie zaprojektowany tak, aby naśladować warunki propagacji w pobliżu odległości Rayleigha od metapowierzchni. W obszarze pola bliskiego propagacja ma charakter sferyczny, z zmieniającymi się wartościami fazy i niejednorodnymi parametrami przestrzennymi w obrębie metapowierzchni, w przeciwieństwie do konwencjonalnych modeli propagacji w polu dalekim, które zakładają fale płaskie. W niniejszym badaniu wykorzystano fizyczny model fali sferycznej do zaprojektowania dokładnej macierzy odpowiedzi impulsowej kanału, która uwzględnia fizyczne położenie zarówno elementów metapowierzchni, jak i terminali użytkowników.
Reprezentacja kanałów w oparciu o odległość
Odległość między każdym elementem metapowierzchni a użytkownikiem jest kluczowym czynnikiem podczas konstruowania modelu kanału w obszarze pola bliskiego. Długość pn między lokalizacjami użytkownika pu a n-tym elementem XL-IMS jest określona wzorem:
(5)
Wartość ta różni się znacznie w obrębie metasurfacy ze względu na jej niezwykle dużą powierzchnię, co prowadzi do nierównych czasów propagacji i przesunięć fazowych, które muszą zostać uwzględnione w modelu.
Model kanału fali sferycznej
Współczynnik kanału pomiędzy n-tym elementem a użytkownikiem jest określany przy użyciu modelu fali sferycznej opisanego w równaniu (2), który uwzględnia zależne od odległości zmiany fazy i amplitudy. W celu generowania zbiorów danych i trenowania modelu współczynniki kanałów są pobierane z frameworka ray-tracingu DeepMIMO, który uwzględnia propagację wielodrożną zarówno w warunkach widoczności bezpośredniej (LoS), jak i braku widoczności bezpośredniej (NLoS), w oparciu o realistyczne struktury środowiskowe. Na podstawie powyższych danych można skonstruować macierz odpowiedzi kanału z odpowiedzi kanałowych dla wszystkich N elementów metapowierzchni, zgodnie z równaniem (3).
Przestrzenna niestacjonarność kanałów pola bliskiego i istotność modelu
W przeciwieństwie do kanałów pola dalekiego, kanały pola bliskiego wykazują niestacjonarność przestrzenną, w której sygnały odbióru różnią się w obrębie metapowierzchni, co umożliwia ogniskowanie wiązki zależne od lokalizacji. Ta informacja przestrzenna jest zawarta w macierzy odpowiedzi kanału i jest przechwytywana przez model GNN-attention. Rycina uzupełniająca 1 przedstawia proponowany model kanału pola bliskiego dla systemu komunikacji 6G opartego na XL-IMS oraz generowanie realistycznych danych kanałowych dla formowania wiązki w oparciu o głębokie uczenie.
Generowanie zestawu danych
Generowanie danych przeprowadzono przy użyciu frameworka DeepMIMO, który zapewnia realistyczne symulacje kanałów oparte na śledzeniu promieni (ray-tracing), uwzględniające geometrię otoczenia, spójność przestrzenną oraz mobilność użytkowników. Odpowiedzi kanałowe dla wielu lokalizacji użytkowników w polu bliskim wygenerowano dla scenariusza zewnętrznego O1 przy częstotliwości nośnej 28 GHz i szerokości pasma 100 MHz. XL-IMS składał się z 256 równomiernie rozłożonych elementów odbijających, a lokalizacje użytkowników wybrano w obszarze pola bliskiego na podstawie kryterium Rayleigha. Współczynniki kanałowe, wzmocnienia ścieżek oraz informacje przestrzenne wyekstrahowano za pomocą biblioteki Python DeepMIMO. Tabela uzupełniająca 1 podsumowuje parametry symulacji i odtwarzalności. Model został wytrenowany przy użyciu uczenia nadzorowanego, gdzie wejściami były cechy stanu kanału, a etykietami docelowymi zoptymalizowane wagi formowania wiązki oparte na fizyce.
Tworzenie cech wejściowych i struktura zbioru danych
Następnie surowe dane z kanałów przekształcono w cechy służące jako dane wejściowe dla modeli głębokiego uczenia. Cechy te obejmowały amplitudy i fazy odpowiedzi kanałów jako oddzielne komponenty. W zależności od wymagań problemu można było uwzględnić inne parametry wejściowe. Na przykład jako cechy można było rozważyć współrzędne użytkowników, odległości oraz niektóre parametry przestrzenne. Zbiór danych podzielono na zbiór treningowy i testowy w stosunku 80:20.
Przygotowanie zestawu danych do uczenia
Ostatni krok polegał na normalizacji przygotowanego zbioru danych i przekształceniu go w tensory odpowiednie do wprowadzenia do architektury GNN-attention w celu przeprowadzenia uczenia. Zbiór treningowy został wykorzystany do nauczenia odwzorowania stanu kanału na optymalne wagi formowania wiązki, natomiast zbiór testowy posłużył do oceny skuteczności modelu. Tak przygotowany zbiór danych sprawił, że proces uczenia był łatwiejszy i bardziej efektywny. Tabela uzupełniająca 2 przedstawia szczegółowy opis zbioru danych DeepMIMO.
Reprezentacja cech
Etap reprezentacji cech przekształcił złożone dane kanału pola bliskiego w format odpowiedni dla głębokiego uczenia. Odpowiedzi kanału zostały rozłożone na cechy, które oddają charakterystykę przestrzenną i propagacyjną, zachowując jednocześnie kluczowe informacje o polu bliskim niezbędne do procesu uczenia. Rycina uzupełniająca 2 ilustruje zbiór danych DeepMIMO wykorzystany w tym procesie. Dane kanału z otoczenia pola bliskiego zostały przedstawione jako liczby zespolone hn, których części rzeczywista i urojona reprezentują informacje o fazie i amplitudzie współczynników kanału. Aby umożliwić efektywne wykorzystanie w algorytmach uczenia maszynowego, współczynniki zostały rozłożone na części rzeczywistą (Rehn) i urojoną (Imhn). Dane można było również przedstawić jako moduł hn oraz fazę współczynników ∠hn.
Ekstrakcja cech przestrzennych
Oprócz współczynnika kanału uwzględniono cechy przestrzenne, aby odzwierciedlić geometryczny charakter środowiska pola bliskiego. Obejmowały one odległość dn od każdego elementu metapowierzchni do odbiornika, a także współrzędne odbiornika (x, y, z). W ten sposób model uwzględniał fakt, że charakterystyka sygnału może zależeć od położenia przestrzennego pary nadajnik-odbiornik, co jest właściwością istotną dla propagacji w polu bliskim.
Normalizacja i przygotowanie danych wejściowych modelu
Jako krok wstępny przed wprowadzeniem danych do modelu przeprowadzono proces normalizacji, aby uniknąć niestabilności numerycznej i przyspieszyć zbieżność. Strukturyzowane wektory podzielono na zbiory danych treningowych i testowych, oznaczone jako (Xtrain, Xtest). W ten sposób dane wejściowe zostały następnie przekazane do modeli GNN i attention.
Na podstawie reprezentacji cech przedstawionej na Rysunku uzupełniającym 2, w zbiorze danych zastosowano różne rodzaje cech, w tym współczynniki kanałów hn, które określały złożoną odpowiedź sygnału wszystkich elementów metasurface. Inną cechą była odległość dn, która w tym przypadku wskazywała na efekty propagacji. Chociaż składowe rzeczywiste i urojone zapewniają numeryczny opis sygnału, to moduł i faza dają odpowiednio jawny opis natężenia sygnału i przesunięcia fazowego. Współrzędne pozycji użytkownika (x, y, z) stanowiły kolejną istotną cechę, która zwiększyła orientację przestrzenną modelu.
Projekt głębokiego uczenia z wykorzystaniem GNN do modelowania zależności przestrzennych zintegrowany z hybrydowym mechanizmem uwagi do uczenia kontekstu globalnego
Ten algorytm uczenia został opracowany w celu efektywnego trenowania złożonego kanału XL-IMS w polu bliskim, z uwzględnieniem jego charakterystyk przestrzennych i propagacyjnych, przy użyciu GNN oraz mechanizmu uwagi. Ponieważ elementy metasurface są ze sobą silnie skorelowane, a ich oddziaływania w strefie pola bliskiego są niejednorodne, sieć w pełni połączona nie jest odpowiednia do prawidłowej reprezentacji kanału. W związku z tym, w tym podejściu metasurface jest reprezentowana jako graf, w którym każdy element pełni rolę węzła, a oddziaływania między nimi stanowią krawędzie. GNN uchwyciła lokalne zależności przestrzenne między sąsiednimi węzłami, natomiast mechanizm uwagi poprawił proces uczenia się zależności globalnych.
Reprezentacja graficzna XL-IMS
XL-IMS jest tworzony jako graf G = (V, E), w którym każdy wierzchołek vn ∈ V reprezentuje jeden element meta-powierzchni, natomiast krawędzie eij ∈ E reprezentują relacje przestrzenne między elementami. Każdy wierzchołek posiada wektor cech xn, który jest uzyskiwany w fazie reprezentacji cech. Macierz sąsiedztwa A oddaje łączność między wierzchołkami, które znajdują się blisko siebie lub w ramach określonych ograniczeń odległości. Taka sformułowanie grafowe jest wykorzystywane w celu uwzględnienia informacji przestrzennych istniejących w propagacji pola bliskiego.
Aby zapewnić powtarzalność, łączność grafu określono przy użyciu podejścia K-najbliższych sąsiadów (K-Nearest Neighbors) z wartością k równą 8. Niech pi oraz pj reprezentują wektory położenia odpowiednio i-tej i j-tej komórki XL-IMS. Odległość między dowolnymi dwoma wierzchołkami mierzono za pomocą metryki euklidesowej zgodnie z wzorem: dij = ||pi − pj||2. Krawędź zostanie utworzona między wierzchołkami i oraz j, jeśli wierzchołek j znajduje się w grupie ośmiu najbliższych sąsiadów wierzchołka i. W związku z tym macierz sąsiedztwa skonstruowano według następującej reguły: Aij = 1, jeśli j ∈ Nk(i), a w przeciwnym razie Aij = 0.
W proponowanej sieci GNN aktualizacje dla każdego węzła są obliczane na podstawie informacji z innych sąsiednich węzłów. W szczególności propagację warstwy l można przedstawić w następujący sposób:
(6)
gdzie hi(l) oznacza wektor cech wierzchołka i w warstwie l, N(i) to sąsiedztwo wierzchołka i, W(l) jest macierzą wag podlegającą uczeniu, a σ(·) symbolizuje nieliniową funkcję aktywacji.
Model GNN-attention składał się z trzech splotowych warstw grafowych z wymiarami ukrytymi ustawionymi na 128, a jako funkcja aktywacji została zastosowana jednostka liniowa z prostowaniem (ReLU). Generowanie grafów przeprowadzono metodą k-najbliższych sąsiadów (k=8). Każda cecha XL-IMS pełni rolę węzła, natomiast krawędzie są tworzone na podstawie odległości euklidesowej między sąsiednimi węzłami. W celu wyodrębnienia globalnych zależności przestrzennych wykorzystano wielogłowicową sieć samoatencji z 4 głowicami. Prawdopodobieństwo dropout ustawiono na 0,3, aby zapobiec przeuczeniu modelu.
Mechanizm uwagi
Aby zamodelować zależności długodystansowe występujące poza zakresem lokalnego sąsiedztwa, w modelu wprowadzono mechanizm uwagi (attention), który działa w następujący sposób:
(7)
gdzie eij = LeakyReLU(aT[Whi | Whj]) określa istotność węzłów i oraz j. Dzięki temu modelowi można przypisać różnym węzłom różny poziom istotności w zależności od ich roli, co ułatwia modelowi zrozumienie kontekstu globalnego w obrębie metapowierzchni. Hybrydowe ramy przechwytują zarówno lokalne, jak i globalne informacje przestrzenne, a wyuczone cechy są przetwarzane przez warstwy w pełni połączone w celu wygenerowania wag formowania wiązki. Jak pokazano na rysunkach uzupełniających 3 i 4, GNN modeluje lokalne zależności przestrzenne poprzez przekazywanie komunikatów i agregację cech, podczas gdy mechanizm uwagi przechwytuje globalne interakcje między elementami metapowierzchni, umożliwiając efektywne formowanie wiązki w polu bliskim. Wagi uwagi są obliczane w celu ilościowego określenia istotności par węzłów, przypisując większą wagę bardziej informacyjnym interakcjom. Wykorzystując te wagi, każdy węzeł agreguje informacje ze wszystkich pozostałych węzłów, co pozwala modelowi uchwycić globalne zależności wykraczające poza lokalne sąsiedztwa. Wynikowa ważona fuzja cech generuje zoptymalizowaną reprezentację cech, która podkreśla istotne charakterystyki przestrzenne, jednocześnie tłumiąc mniej relewantne informacje, co wspiera efektywną optymalizację formowania wiązki oraz ogniskowanie energii.
Optymalizacja kształtowania wiązki
Etap optymalizacji formowania wiązki określił przesunięcia fazowe XL-IMS oraz wagę formowania wiązki w celu skupienia energii na użytkownikach w polu bliskim. W przeciwieństwie do konwencjonalnego sterowania wiązką w polu dalekim, formowanie wiązki w polu bliskim koncentruje energię w konkretnych lokalizacjach przestrzennych. W proponowanej strukturze model GNN-attention nauczył się optymalnych parametrów formowania wiązki bezpośrednio z charakterystyk kanału, co zmniejszyło złożoność iteracyjnych metod optymalizacji. Sygnał odebrany przez użytkownika przedstawiono w równaniu (4).
Dane treningowe dla każdej realizacji kanału pobrano z DeepMIMO, a odpowiadające im przesunięcia fazowe XL-IMS wyznaczono przy użyciu algorytmu AO. Macierz kanału pola bliskiego h oraz macierz przesunięć fazowych XL-IMS zostały wykorzystane do obliczenia wektora przesunięcia fazowego θ poprzez maksymalizację mocy sygnału odbiorczego:
(8)
(9)
Optymalizator AO był uruchamiany przez 100 iteracji, aż do osiągnięcia pożąnego progu zbieżności wynoszącego 10−4. Zoptymalizowane wartości te posłużyły jako etykiety treningowe dla naszego modelu GNN-attention. Celem było uzyskanie interferencji konstruktywnej dla efektywnego ogniskowania wiązki w polu bliskim oraz zwiększenia mocy sygnału odbieranego. Zamiast konwencjonalnych metod optymalizacji, takich jak SDR lub AO, problem rozwiązano przy użyciu podejścia opartego na uczeniu. W szczególności problem optymalizacji kształtowania wiązki został przybliżony przez model głębokiego uczenia fΘ(⋅).
Do przeszkolenia modelu GNN-attention zastosowano uczenie nadzorowane, przyjmując optymalne wektory przesunięcia fazowego formowania wiązki, wyznaczone przez algorytm AO, jako cele treningowe. W szczególności, niech wektor optymalnego przesunięcia fazowego wygenerowany przez AO będzie oznaczony jako θAO,i, a predykcje sieci jako
. Wówczas kryterium uczenia sformułowano w następujący sposób:
(10)
gdzie N to liczba próbek treningowych.
Algorytm 1: Formowanie wiązki w polu bliskim z wykorzystaniem modelu GNN-attention został przedstawiony w Pliku uzupełniającym 1. W tym algorytmie pierwszym krokiem było skonfigurowanie parametrów systemu oraz bazy danych DeepMIMO w celu stworzenia realistycznych warunków kanału. Cechę kanału w polu bliskim dla każdej próbki obliczono przy użyciu modelu fali sferycznej, a następnie z bazy danych wyekstrahowano istotne cechy, takie jak współczynniki kanału i cechy przestrzenne. Wyekstrahowane cechy przekształcono w formę ustrukturyzowaną, a następnie przedstawiono jako graf, w którym elementami węzłów były elementy metapowierzchni. Model głębokiego uczenia wygenerowano z wykorzystaniem warstw GNN w celu uchwycenia lokalnych zależności przestrzennych, a następnie wprowadzono mechanizm uwagi (attention mechanism) do nauki zależności globalnych. Framework głębokiego uczenia wytrenowano metodą uczenia nadzorowanego z zastosowaniem funkcji straty MSE pomiędzy przewidywanymi przesunięciami fazy a etykietami optymalnych przesunięć fazy wygenerowanymi przez AO.
Implementacja została wykonana przy użyciu języka Python i biblioteki PyTorch. Próbki dla kanału DeepMIMO pobrano ze scenariusza zewnętrznego O1. Współczynniki tych kanałów przekształcono w grafy, a następnie przetworzono za pomocą GNN z trzema warstwami i mechanizmem wielogłowicowej uwagi (multihead attention). Model trenowano przy użyciu optymalizatora Adam z funkcją straty MSE. Model ten wygenerował przesunięcia fazowe i wagi dla formowania wiązki (beamforming) dla użytkowników w polu bliskim, a wyniki oceniano pod kątem przepływności, wydajności i SINR.