Предлагаемая концепция моделирования каналов ближней зоны и формирования диаграммы направленности в системах 6G с поддержкой XL-IMS сочетает в себе физическое моделирование с оптимизацией на основе машинного обучения для преодоления ограничений традиционных подходов для дальней зоны. Разработана модель канала ближней зоны, учитывающая распространение сферических волн и пространственную нестационарность, которая была валидирована с использованием набора данных DeepMIMO. Информация о канале представлена в виде графа для фиксации пространственных зависимостей между элементами метаповерхности и пользователями; при этом для обучения весам формирования луча при фокусировке в ближней зоне применяется сквозная сеть GNN-attention. В отличие от традиционных итерационных методов оптимизации, данная концепция на основе обучения напрямую предсказывает решения для формирования луча, что снижает вычислительную сложность при сохранении высокой точности и адаптивности к динамическим условиям канала. Эффективность оценивается по достижимой скорости передачи данных, коэффициенту усиления при формировании луча и энергоэффективности, что демонстрирует превосходные результаты по сравнению с традиционными методами формирования луча для дальней зоны. На Рисунке 1 представлена общая схема концепции.
Ниже приведены предлагаемые этапы:
Шаг 1: Конфигурация среды DeepMIMO
Платформа DeepMIMO v3.x сначала настраивается с помощью конфигурации O1 Outdoor, в которой используется несущая частота 28 GHz и ширина полосы пропускания 100 MHz. В симуляции XL-IMS представлена в виде равномерной плоской антенной решетки 16×16 из 256 отражающих элементов с интервалом λ/2 между ними. В качестве модели распространения выбрана сферическая волна с путями прямой (LoS) и непрямой (NLoS) видимости.
Шаг 2: Определение геометрии XL-IMS и местоположения пользователей
Разработайте XL-IMS в виде плоской метаповерхности размером 16×16, состоящей из 256 реконфигурируемых отражающих элементов. Пользовательские терминальные устройства были случайным образом распределены в ближней зоне, на расстоянии от 0.5m до 10m от XL-IMS, согласно критерию расстояния Рэлея.
Шаг 3: Генерация ближних полей с помощью распространения сферических волн
Моделируйте каналы между каждым элементом перенастраиваемой метаповерхности и пользователем с помощью распространения сферических волн, что позволяет учесть изменение фазы, вызванное расстоянием.
Шаг 4: Генерация данных ближнего поля для канала
Пользовательские терминалы были равномерно распределены в ближней зоне (0,5–10 м) от панели XL-IMS. Для формирования каналов использовались как пути распространения с прямой видимостью (LoS), так и без прямой видимости (NLoS), созданные с помощью сценария DeepMIMO O1 Outdoor.
Шаг 5: Подготовка набора данных для обучения и тестирования
Сгенерированный набор данных состоял из 20 000 образцов, из которых 16 000 (80%) были использованы для обучения и 4 000 (20%) — для тестирования.
Шаг 6: Нормализация входных признаков и подготовка входных данных для модели
Заключительный этап процесса подготовки включал выполнение масштабирования Min-Max для всех входных признаков, таких как амплитуда канала, фаза, расстояние и координаты пользователя. В ходе этого процесса данные были масштабированы и преобразованы в тензоры, которые затем подавались в сеть GNN-attention. Оптимизация модели проводилась с использованием оптимизатора Adam при скорости обучения 0.001 и количестве эпох 50.
Шаг 7: Применение механизма внимания
Добавьте механизм внимания для учета глобальной информации графа и захвата дальних зависимостей между узлами графа. Это помогает модели сосредоточиться на важных пространственных признаках.
Шаг 8: Предсказание весовых коэффициентов диаграммообразования
Предскажите весовые коэффициенты диаграммы направленности и значения фазового сдвига с помощью вложения на основе GNN-attention. Сложные, трудоемкие и итерационные процедуры формирования луча заменяются методом прогнозирования.
Шаг 9: Оценка показателей эффективности
Для оценки эффективности связи и формирования луча рассчитайте такие ключевые показатели эффективности (KPI), как достижимая скорость, коэффициент усиления формирования луча, спектральная эффективность, энергоэффективность и SINR.
Шаг 10: Сравнение с эталонными методами
Сравните предлагаемый алгоритм с традиционными алгоритмами, включая ZF, MMSE, SDR и AO.
Определение модели системы
Предлагаемая система рассматривает беспроводную сеть 6G с поддержкой XL-IMS, состоящей из 256 реконфигурируемых отражающих элементов, которые управляют распространением электромагнитных волн посредством регулируемых фазовых сдвигов. Поскольку пользователи находятся в области ближнего поля, традиционное допущение о плоской волне дальнего поля становится недействительным, что требует учета распространения сферических волн, затухания в зависимости от расстояния и пространственной нестационарности. Граница между областями ближнего и дальнего поля определяется расстоянием Рэлея, которое зависит от размера метаповерхности и длины волны несущей частоты, и отмечает переход от эффектов сферического волнового фронта к распространению плоской волны. Это описывается следующим образом:
(1)
Здесь D — размер метаповерхности, λ — длина волны, а R — расстояние Рэлея. В канале связи возникает эффект ближнего поля, когда пользователь находится в точке, расстояние до которой меньше истинного значения R. По мере прохождения сигнала через различные компоненты метаповерхности он подвергается определенным изменениям фазы и амплитуды.
Канал от n-й метаповерхности к пользователю использует модель сферической волны для распространения в ближней зоне, которая представляется следующим уравнением коэффициента канала:
(2)
Здесь α обозначает параметр потерь при распространении сигнала, λ — длину волны несущей, а dn — расстояние от пользователя до компонента.
В данном случае, βn является ли коэффициент потерь при распространении для nth элемент XL-IMS, который можно представить как βn = дн−α, где α обозначает показатель потерь при распространении сигнала. Следует отметить, что дн представляет собой евклидово расстояние от положения пользователя до нth элемент метаповерхности. Изменения канала в пространстве обусловлены значительными различиями в дn между различными метаповерхностями, в отличие от ситуации с каналами в дальней зоне. В общем виде вектор канала можно представить следующим образом:
(3)
Входной сигнал, получаемый пользователем, зависит от фазовой настройки системы XL-IMS. Пусть θ = [θ1, θ2, …., θN] представляет фазовый сдвиг, реализованный с помощью метаповерхности. Тогда входной сигнал можно выразить следующим образом:
(4)
где Φ = diag(ejθ1, ….., ejθN) — передаваемый сигнал, w — вектор предварительного кодирования базового кода, а n обозначает шум. Вектор канала в ближней зоне определяется как h ∈ ℂN×1, диагональная матрица фазового сдвига XL-IMS как Φ ∈ ℂN×N, вектор диаграммообразования как x ∈ ℂN×1, а принимаемый сигнал как y ∈ ℂ. Аддитивный член шума определяется как n ∼ CN (0, σ2). Данные размерности обеспечивают математическую корректность формулировки.
В данном случае ws представляет собой вектор предкодирования базовой станции, который используется для предкодирования передаваемого символа (s). При проведении сравнительного анализа производительности на базовой станции реализуются общепринятые методы цифрового предкодирования, такие как ZF и MMSE, с использованием аналогичных реализаций каналов ближнего поля. Однако новый подход направлен на оптимизацию матрицы фазовых сдвигов XL-IMS (θ) с помощью механизма обучения GNN-attention.
Для концентрации энергии в ближней зоне необходимо обеспечить достижимость для пользователя такого дизайна θ, при котором возникает конструктивная интерференция. Стратегия формирования луча направлена на максимизацию спектральной эффективности или мощности принимаемого сигнала. На рисунке 2 представлена архитектура модели системы связи в ближней зоне, основанной на подходе XL-IMS.
Моделирование каналов в ближней зоне
Предложенная модель ближнего поля для сетей связи 6G на базе XL-IMS была специально разработана для имитации условий распространения сигнала вблизи расстояния Рэлея от метаповерхности. В области ближнего поля распространение является сферическим, с изменяющимися значениями фазы и неоднородными пространственными параметрами по всей метаповерхности, в отличие от традиционных моделей распространения в дальнем поле, в которых предполагаются плоские волны. В настоящем исследовании была использована физическая модель сферической волны для разработки точной матрицы импульсного отклика канала, учитывающей физическое расположение как элементов метаповерхности, так и пользовательских терминалов.
Представление каналов на основе расстояния
Расстояние между каждым компонентом метаповерхности и пользователем является решающим фактором при построении модели канала в режиме ближнего поля. Расстояние pn между координатами пользователя pu и n-м компонентом XL-IMS определяется следующим образом:
(5)
Это значение значительно варьируется по всей метаповерхности из-за её чрезвычайно большой площади, что приводит к неравномерному времени распространения и фазовым сдвигам, которые необходимо учитывать в модели.
Модель канала со сферическими волнами
Коэффициент канала между n-м элементом и пользователем определяется с использованием модели сферической волны, описанной в уравнении (2), которая учитывает зависимости фазы и амплитуды от расстояния. Для генерации набора данных и обучения модели коэффициенты канала извлекаются из среды трассировки лучей DeepMIMO, которая учитывает многолучевое распространение сигнала как в условиях прямой видимости (LoS), так и в условиях отсутствия прямой видимости (NLoS) на основе реалистичных структур окружающей среды. Исходя из вышеизложенного, можно построить матрицу отклика канала на основе откликов всех N элементов метаповерхности, как показано в уравнении (3).
Пространственная нестационарность каналов ближней зоны и значимость моделей
В отличие от каналов дальней зоны, каналы ближней зоны характеризуются пространственной нестационарностью, при которой принимаемые сигналы варьируются по поверхности метаповерхности, что позволяет осуществлять фокусировку луча в зависимости от местоположения. Эта пространственная информация заложена в матрице отклика канала и фиксируется моделью GNN-attention. На Дополнительном рисунке 1 представлена предлагаемая структура моделирования каналов ближней зоны для системы связи 6G на базе XL-IMS и процесс генерации реалистичных данных о канале для формирования луча на основе глубокого обучения.
Создание набора данных
Генерация данных осуществлялась с использованием платформы DeepMIMO, которая обеспечивает реалистичное моделирование каналов на основе трассировки лучей с учетом геометрии окружающей среды, пространственной согласованности и мобильности пользователей. Импульсные характеристики каналов для нескольких положений пользователей в ближней зоне были сформированы с использованием сценария O1 (открытая местность) при несущей частоте 28 GHz и ширине полосы 100 MHz. XL-IMS состоял из 256 равномерно распределенных отражающих элементов; положения пользователей в ближней зоне выбирались на основе критерия Рэлея. Коэффициенты канала, усиления путей и пространственная информация были извлечены с помощью библиотеки DeepMIMO для Python. В Дополнительной таблице 1 обобщены параметры моделирования и воспроизводимости. Модель обучалась методом обучения с учителем, где в качестве входных данных использовались характеристики состояния канала, а в качестве целевых меток — оптимизированные на основе физических принципов веса диаграммы направленности.
Создание входных признаков и структура набора данных
Затем необработанные данные каналов были преобразованы в признаки для подачи в модели глубокого обучения. Эти признаки включали амплитуды и фазы отклика каналов в качестве отдельных компонентов. В зависимости от требований задачи могли быть включены и другие входные параметры. Например, в качестве признаков могли быть рассмотрены координаты пользователей, расстояния и некоторые пространственные параметры. Набор данных был разделен на обучающую и тестовую выборки в соотношении 80:20.
Подготовка набора данных для обучения
Заключительный этап включал нормализацию подготовленного набора данных и его преобразование в тензоры, подходящие для ввода в архитектуру GNN-attention для обучения. Обучающий набор данных использовался для изучения отображения состояния канала на оптимальные веса диаграммообразования, в то время как тестовый набор данных служил для оценки эффективности модели. Тщательная подготовка данных упростила и повысила эффективность процесса обучения. В Дополнительной таблице 2 приведены подробные сведения о наборе данных DeepMIMO.
Представление признаков
Этап представления признаков позволил преобразовать данные о ближнепольных каналах с комплексными значениями в формат, подходящий для глубокого обучения. Отклики каналов были разложены на признаки, которые отражают пространственные и пропогационные характеристики, сохраняя при этом важную информацию о ближнем поле для процесса обучения. Дополнительный рисунок 2 иллюстрирует набор данных DeepMIMO, использованный в этом процессе. Данные о канале связи в ближней зоне были представлены в виде комплексных чисел чn, где действительная и мнимая составляющие представляют информацию о фазе и амплитуде коэффициентов канала. Для обеспечения эффективного использования в алгоритмах машинного обучения коэффициенты были разложены на действительную (Reчн) и мнимую (Imчн) частей. Данные также могут быть представлены в виде величины чн и фаза коэффициентов ∠чн.
Извлечение пространственных признаков
Помимо коэффициента канала, были включены пространственные признаки, чтобы учесть геометрическую природу среды ближнего поля. К ним относятся расстояние dn от каждого элемента метаповерхности до приемника, а также координаты приемника (x, y, z). Таким образом, модель учитывала тот факт, что характеристики сигнала могут зависеть от пространственного положения пары передатчик-приемник, что является важным свойством распространения сигнала в ближнем поле.
Нормализация и подготовка входных данных модели
В качестве предварительного этапа перед подачей входных данных в модель была выполнена процедура нормализации для предотвращения численной нестабильности и ускорения сходимости. Структурированные векторы были разделены на обучающую и тестовую выборки, обозначенные как (Xtrain, Xtest). После этого входные данные были поданы в модели GNN и attention.
Исходя из представления признаков на Дополнительном рисунке 2, в данном наборе данных использовались признаки различных типов, включая коэффициенты каналов hn, которые указывали на комплексный отклик сигнала всех элементов метаповерхности. Другим признаком было расстояние dn, которое в данном случае указывало на эффекты распространения. Хотя мнимые и вещественные компоненты обеспечивают численное описание сигнала, амплитудные и фазовые составляющие дают явное описание силы сигнала и фазового сдвига соответственно. Координаты положения пользователя (x, y, z) были еще одним важным признаком, который расширил пространственное восприятие модели.
Архитектура глубокого обучения с использованием GNN для моделирования пространственных зависимостей, интегрированная с гибридным механизмом внимания для изучения глобального контекста
Данный алгоритм обучения был разработан для эффективного обучения сложного канала XL-IMS в ближней зоне с учетом его пространственных характеристик и характеристик распространения с использованием GNN и механизма внимания. Поскольку элементы метаповерхности сильно коррелируют друг с другом, а их взаимодействия в зоне ближнего поля неоднородны, полносвязная сеть не подходит для адекватного представления канала. Следовательно, в данном подходе метаповерхность представлена в виде графа, где каждый элемент выступает в качестве узла, а взаимодействия между ними — в качестве ребер. GNN позволила зафиксировать локальные пространственные зависимости между соседними узлами, в то время как механизм внимания улучшил процесс обучения глобальным зависимостям.
Графическое представление XL-IMS
XL-IMS создается в виде графа G = (V, E), где каждая вершина vn ∈ V представляет собой один элемент метаповерхности, а ребра eij ∈ E представляют пространственные связи между элементами. Каждая вершина имеет вектор признаков xn, который получается на этапе представления признаков. Матрица смежности A отражает связность между вершинами, которые расположены близко друг к другу или находятся в пределах определенных ограничений по расстоянию. Такая формулировка графа используется для учета пространственной информации, существующей при распространении в ближнем поле.
Для обеспечения воспроизводимости связность графа определяется с использованием метода k-ближайших соседей, где k равно 8. Пусть pi и pj представляют собой векторы координат i-й и j-й клеток XL-IMS соответственно. Расстояние между любыми двумя вершинами измеряется с помощью евклидовой метрики следующим образом: dij = ||pi − pj||2. Ребро между вершинами i и j будет существовать, если вершина j входит в группу восьми ближайших соседей вершины i. Таким образом, матрица смежности строится по следующему правилу: Aij = 1, если j ∈ Nk(i), и Aij = 0 в противном случае.
В предлагаемой GNN обновления для каждого узла рассчитываются на основе информации из других соседних узлов. В частности, распространение слоя l может быть представлено следующим образом:
(6)
где hi(l) представляет вектор признаков вершины i в слое l, N(i) — окрестность вершины i, W(l) — обучаемую матрицу весов, а σ(·) обозначает нелинейную функцию активации.
Модель GNN-attention состояла из трех графовых сверточных слоев со скрытыми размерностями 128; в качестве функции активации использовалась линейная функция с выпрямлением (ReLU). Генерация графа осуществлялась методом k-ближайших соседей (k=8). Каждый признак XL-IMS выступает в роли узла, а ребра формируются на основе евклидова расстояния между соседними узлами. Для извлечения глобальных пространственных зависимостей мы использовали многоголовую сеть самовнимания с 4 головами. Вероятность исключения (dropout) была установлена на уровне 0,3 для предотвращения переобучения.
Механизм внимания
Для моделирования дальних зависимостей, выходящих за пределы локальной окрестности, в модель был введен механизм внимания, который работает следующим образом:
(7)
где eij = LeakyReLU(aT[Whi | Whj]) — это значимость узлов i и j. С помощью этой модели можно присваивать узлам разный уровень значимости в зависимости от их ролей, что облегчает понимание моделью глобального контекста по всей метаповерхности. Гибридная архитектура фиксирует как локальную, так и глобальную пространственную информацию, а извлеченные признаки обрабатываются через полносвязные слои для генерации весов диаграммообразования. Как показано на дополнительных рисунках 3 и 4, GNN моделирует локальные пространственные зависимости посредством передачи сообщений и агрегации признаков, в то время как механизм внимания фиксирует глобальные взаимодействия между элементами метаповерхности, обеспечивая эффективное формирование луча в ближнем поле. Веса внимания вычисляются для количественной оценки релевантности пар узлов, при этом более информативным взаимодействиям присваивается большая значимость. Используя эти веса, каждый узел агрегирует информацию от всех остальных узлов, что позволяет модели улавливать глобальные зависимости за пределами локальных окрестностей. Результирующее взвешенное слияние признаков формирует оптимизированное представление признаков, которое подчеркивает важные пространственные характеристики и подавляет менее релевантную информацию, что способствует эффективной оптимизации формирования луча и фокусировке энергии.
Оптимизация формирования луча
На этапе оптимизации диаграммообразования были определены фазовые сдвиги XL-IMS и веса диаграммообразования для фокусировки энергии на пользователях в ближней зоне. В отличие от традиционного управления лучом в дальней зоне, диаграммообразование в ближней зоне концентрирует энергию в конкретных пространственных точках. В предложенной структуре модель GNN-attention позволила определить оптимальные параметры диаграммообразования непосредственно на основе характеристик канала, что снизило сложность итерационных методов оптимизации. Сигнал, принимаемый пользователем, представлен в уравнении (4).
Обучающие данные для каждой реализации канала были получены из DeepMIMO, а соответствующие фазовые сдвиги XL-IMS были определены с помощью алгоритма AO. Матрица канала ближнего поля h и матрица фазовых сдвигов XL-IMS использовались для вычисления вектора фазового сдвига θ путем максимизации мощности принимаемого сигнала:
(8)
(9)
Оптимизатор AO запускался на 100 итераций до достижения требуемого порога сходимости 10−4. Эти оптимизированные значения служили обучающими метками для нашей модели GNN-attention. Целью было достижение конструктивной интерференции для эффективной фокусировки луча в ближней зоне и увеличения мощности принимаемого сигнала. Вместо традиционных методов оптимизации, таких как SDR или AO, задача была решена с использованием подхода на основе обучения. В частности, задача оптимизации диаграммы направленности была аппроксимирована моделью глубокого обучения fΘ(⋅).
Для обучения модели GNN-attention был применен метод контролируемого обучения, в котором в качестве целевых значений для обучения использовались оптимальные векторы фазового сдвига диаграммообразования, полученные с помощью алгоритма AO. В частности, обозначим оптимальный вектор фазового сдвига, сгенерированный AO, как θAO,i, а предсказания сети — как
. Тогда критерий обучения был сформулирован следующим образом:
(10)
где N — количество обучающих выборок.
Алгоритм 1: Формирование луча в ближней зоне с использованием модели GNN-attention представлен в Дополнительном файле 1. В данном алгоритме первым шагом была настройка параметров системы и базы данных DeepMIMO для создания реалистичных условий канала. Характеристики канала ближней зоны для каждого образца рассчитывались с использованием модели сферической волны, после чего из базы данных извлекались важные признаки, такие как коэффициенты канала и пространственные характеристики. Извлеченные признаки были преобразованы в структурированную форму и представлены в виде графа, где узлами являются элементы метаповерхности. Модель глубокого обучения была создана с использованием слоев GNN для захвата локальных пространственных зависимостей, а также был введен механизм внимания (attention mechanism) для изучения глобальных зависимостей. Фреймворк глубокого обучения обучался с помощью контролируемого обучения с использованием функции потерь MSE между предсказанными фазовыми сдвигами и оптимальными метками фазовых сдвигов, сгенерированными методом AO.
Данная реализация была выполнена с использованием Python и PyTorch. Образцы для канала DeepMIMO были получены из сценария O1 для открытого пространства. Эти коэффициенты канала были преобразованы в графы и обработаны с помощью GNN с тремя слоями и механизмом многоголового внимания (multihead attention). Модель была обучена с использованием оптимизатора Adam и функции потерь MSE. С помощью этой модели были сгенерированы фазовые сдвиги и веса для формирования луча для пользователей в ближней зоне, а результаты были оценены по показателям скорости, эффективности и SINR.