Данное исследование включало только математическое моделирование, разработку алгоритмов и численное моделирование для системы OFDM-ISAC. В работе не принимали участия люди, животные, не использовались биологические образцы или персональные данные. Таким образом, одобрение институционального комитета по этике, включая одобрение Институционального наблюдательного совета (IRB) или Институционального комитета по уходу и использованию животных (IACUC), не требовалось.
Разработка модели распределения мощности для интегрированной системы зондирования и связи на базе ортогонального частотного разделения каналов
В системах ISAC традиционные равные PAM часто равномерно распределяют мощность передачи по каждому каналу и антенне в процессе одновременной связи и зондирования, не учитывая их различные требования к производительности. Это может привести к выделению избыточной мощности неэффективным линиям передачи и невозможности удовлетворить гетерогенные требования служб связи и зондирования, что снижает общую эффективность использования мощности12,13,14. Поэтому в данном исследовании рассматриваются двойные функциональные цели системы путем совместной оптимизации пропускной способности канала связи и PDP. Для разработки стратегии распределения мощности для системы ISAC на базе OFDM, называемой моделью OFDM-ISAC-P, построена структура многокритериальной оптимизации. OFDM представляет собой эффективный метод многочастотной модуляции, состоящий из передающей и приемной частей. Процесс OFDM показан на Рисунке 1.

Рисунок 1. Рабочий процесс обработки сигналов формы волны с ортогональным частотным разделением каналов (OFDM), используемой в системе интегрированного зондирования и связи (ISAC). Рабочий процесс включает кодирование данных, перемежение, отображение символов, последовательно-параллельное преобразование, обратное быстрое преобразование Фурье (IFFT), вставку циклического префикса (CP), радиочастотную (РЧ) передачу, распространение по беспроводному каналу, РЧ-прием, удаление CP, быстрое преобразование Фурье (FFT), выравнивание канала, обратное отображение, обратное перемежение, декодирование и восстановление данных. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
На Рисунке 1 сигналы OFDM используются в качестве общего рабочего сигнала для связи и зондирования благодаря их высокой спектральной эффективности и высокой устойчивости к многолучевым помехам, что позволяет обеспечить единую передачу сигналов для обеих функций15,16,17. Без учета эффекта спектрального формирования оконной функции передаваемая форма сигнала OFDM формируется с использованием несущей частоты, количества поднесущих, комплексных весов поднесущих для распределения мощности, кодирования сигнала связи и интервала между поднесущими. Полоса пропускания системы определяется интервалом между поднесущими и их количеством18. Процесс генерации сигнала связи определен следующим образом. Бинарные информационные биты сначала кодируются с использованием схемы кодирования с малой плотностью проверок на четность (LDPC) для повышения надежности передачи, а затем отображаются в символы с комплексными значениями с использованием квадратурной фазовой манипуляции (QPSK). Модулированные символы распределяются по поднесущим OFDM в соответствии с оптимизированной стратегией распределения мощности. После применения быстрого преобразования Фурье (БПФ) и добавления циклического префикса (CP) формируется временная форма сигнала связи OFDM, которая передается вместе с сигналом зондирования. Одна и та же конфигурация модуляции и кодирования сохраняется во всех симуляциях. Конфигурация системы ISAC, использованная в данном исследовании, показана на Рисунке 2.

Рисунок 2. Конфигурация системы интегрированного зондирования и связи (ISAC) на базе ортогонального частотного разделения каналов (OFDM). На диаграмме показаны передающая (Tx) антенная решетка, приемная (Rx) антенная решетка, прямой путь связи в пределах прямой видимости (LoS), путь зондирования по эху от цели, путь отражения от земли, путь отражения от неба и типичные препятствия в моделируемой среде распространения сигнала. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Система состоит из передающей антенной решетки и приемной антенной решетки. Передача сигнала происходит по двум основным путям. Первый — это путь прямой видимости (LoS), при котором сигнал передается непосредственно от передающей решетки к приемной. Второй — путь отражения, при котором часть сигнала отражается от земли и достигает приемной решетки, в то время как другая часть проходит через путь отражения от неба перед попаданием на приемную решетку19,20,21. Система представляет собой многоантенную систему OFDM-ISAC с доминирующим путем LoS, в которой приемники связи и зондирования могут иметь разные конфигурации антенн. Хотя на Рисунке 2 показаны возможные компоненты распространения, включая прямой путь LoS и пути, связанные с отражением, в математической модели и при моделировании учитываются только доминирующий канал связи LoS и канал зондирования эхо-сигнала от цели. Отражения от земли и неба рассматриваются как фоновые эффекты распространения, а не как независимо оптимизируемые каналы без прямой видимости (NLoS). Матрица отклика канала определена с учетом доминирования пути LoS, при этом многолучевые эффекты представлены через замирания, доплеровский сдвиг и шумовые члены. На основе этой матрицы отклика канала, совместно с передаваемым сигналом и моделью шума в канале, принимаемый сигнал связи выражен в Уравнении 1.

Здесь nc,m представляет собой гауссовский белый шум, Sm обозначает передаваемый сигнал, а Hc,m представляет собой матрицу отклика канала связи.
В функции зондирования системы ISAC путь передачи сигнала состоит из двух этапов. Первый этап — это передача от передатчика к цели, а второй — путь эха, по которому отраженный сигнал распространяется от цели обратно к приемнику через канал зондирования22,23,24. Для m-й поднесущей матрица отклика канала зондирования характеризует распространение сигнала и выражается Уравнением 2:

Здесь bs представляет собой коэффициент замирания эхо-сигнала зондирования, fs,m обозначает отраженный доплеровский сдвиг частоты, а τs представляет собой задержку распространения сигнала по пути отражения. На основе матрицы отклика канала зондирования принимаемый сигнал зондирования выражается Уравнением 3:

Здесь ns,m представляет собой гауссов белый шум в канале измерения.
В системах беспроводной связи эффективное планирование спектра и распределение ресурсов имеют решающее значение для повышения общей производительности системы25. Соответственно, в данном исследовании в качестве основного критерия проектирования используется пропускная способность канала связи. Для максимизации пропускной способности канала разработана модель оптимизации, которая распределяет общую мощность передачи по поднесущим OFDM согласно предлагаемой стратегии распределения ресурсов, как показано в Уравнении 4.
Задача оптимизации связи формулируется следующим образом (Уравнение 4):

подвержен

Здесь σc2 представляет собой дисперсию шума коммуникационного гауссовского канала, P обозначает общую мощность передачи данных, Pc представляет вектор оптимального распределения мощности связи, am обозначает вес, присвоенный m-й поднесущей, а Nc представляет количество поднесущих.
Оптимальное решение, полученное с помощью WA, выражается следующим образом: Уравнение 5:

Здесь
обозначает оптимальную мощность, выделенную для m-го коммуникационного подканала, gm представляет собой отношение коэффициента усиления канала к шуму m-го коммуникационного подканала, а μc обозначает уровень воды WA.
Согласно алгоритму WA, проиллюстрированному на Рисунке 3, подканалы с более высокими коэффициентами усиления канала, что соответствует меньшим обратным коэффициентам усиления, получают большую мощность передачи при выбранном уровне заполнения (water level). Это позволяет системе использовать благоприятные условия канала при соблюдении ограничения по общей мощности и повышении общей пропускной способности связи. Следовательно, подканалы с лучшими условиями распространения получают больше мощности, в то время как подканалы с худшими условиями — меньше. Благодаря этой стратегии адаптивного распределения мощности алгоритм WA эффективно использует ресурсы канала для повышения качества связи.

Рисунок 3. Принцип алгоритма заполнения водой (WA), используемого для распределения мощности связи. Горизонтальная ось представляет подканалы OFDM, а вертикальная ось — обратный коэффициент усиления канала и выделяемую мощность передачи. Уровень воды определяет мощность передачи, назначенную каждому подканалу при ограничении общей мощности передачи. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Для достижения целей зондирования при фиксированной вероятности ложной тревоги крайне важна соответствующая стратегия распределения мощности. Оптимизированная стратегия распределения мощности повышает эффективность зондирования за счет увеличения точности обнаружения цели. Процесс зондирования требует проверки гипотез в условиях гипотезы о наличии цели (H1) и гипотезы об отсутствии цели (H0). Согласно этим гипотезам, принимаемый сигнал для каждой поднесущей во время k-го импульса выражается формулой 6:

Здесь y обозначает вектор наблюдений датчика, A представляет собой диагональную матрицу diag(a0, a1, …, aNc−1), D(v) обозначает диагональную матрицу, содержащую доплеровский сдвиг частоты цели, β представляет собой вектор коэффициентов отражения цели, а ns обозначает белый гауссовый шум в канале зондирования.
Затем оцениваются функции плотности вероятности для H0 и H1. Для заданной вероятности ложного срабатывания рассчитывается PDP на основе этих двух гипотез. Согласно теории детектирования по отношению правдоподобия, впоследствии выводится обобщенный детектор максимального отношения правдоподобия для максимизации PDP.
Обобщенный детектор отношения максимального правдоподобия выражается Уравнением 7:

Здесь lG и LG обозначают статистику обнаружения на основе обобщенного отношения правдоподобия и функцию отношения правдоподобия соответственно, yH представляет собой эрмитово-сопряженное транспонирование y, а γ0 обозначает порог обнаружения. В процессе обнаружения порог используется для определения того, какая гипотеза принимается: гипотеза о наличии цели (H1) или гипотеза об отсутствии цели (H0). Зависимость между порогом обнаружения при зондировании и вероятностью обнаружения цели показана на рисунке 4.

Рисунок 4. Взаимосвязь между порогом обнаружения датчиком и вероятностью обнаружения цели. Функции плотности вероятности для гипотезы отсутствия цели и гипотезы присутствия цели показаны вместе с порогом обнаружения, используемым для проверки гипотез. Заштрихованные области соответствуют областям принятия решения, связанным с вероятностью ложной тревоги и вероятностью обнаружения цели. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Как показано на рисунке 4, вероятность обнаружения (PDP) напрямую связана с нецентральным комплексным распределением хи-квадрат. При фиксированной вероятности ложной тревоги (т. е. максимально допустимой вероятности ложного обнаружения) параметр нецентральности положительно коррелирует с PDP. Следовательно, увеличение параметра нецентральности повышает вероятность правильного обнаружения целевого сигнала. Основываясь на этой взаимосвязи, общая мощность передачи при зондировании рассматривается как ограничение по мощности, в то время как максимизация параметра нецентральности выбирается в качестве цели оптимизации. Результирующая задача условной оптимизации используется для вывода оптимальной стратегии распределения мощности зондирования и выражается в уравнении 8:

где α обозначает оптимальный вектор весовых коэффициентов передачи сигнала зондирования, полученный с учетом ограничения по мощности, и также представляет собой вектор весовых коэффициентов передачи для поднесущих зондирования, B(v, β) обозначает вектор информационной составляющей зондирования, а P представляет общую мощность передачи сигнала зондирования.
Оптимизация показывает, что оптимальная стратегия распределения мощности зондирования не предполагает распределение мощности по нескольким подканалам. Вместо этого вся доступная мощность передачи зондирующего сигнала выделяется подканалу, который по результатам оценки канала обеспечивает наилучшие характеристики зондирования, что позволяет максимизировать эффективность обнаружения цели в одном канале.
В системе OFDM-ISAC распределение мощности является основной задачей оптимизации, и оптимальная стратегия должна удовлетворять проектным требованиям как связи, так и зондирования. В частности, оптимизация одновременно учитывает показатели эффективности связи, включая вероятность битовой ошибки и пропускную способность канала, а также показатели эффективности зондирования, включая PDP и точность оценки параметров. Общая производительность системы повышается за счет скоординированного распределения ресурсов. Однако прямое объединение оптимальных стратегий распределения мощности для связи и зондирования с использованием традиционной оптимизации взвешенной суммы имеет существенные ограничения. Поскольку различные физические размерности и относительные вклады двух показателей эффективности учитываются недостаточно, полученная оптимизация может неточно отражать требования интегрированных систем связи и зондирования26,27.
Таким образом, в данном исследовании предлагается интегрированная основа проектирования, основанная на нормализованном взвешивании. Сначала метрики производительности связи и зондирования нормализуются для устранения эффектов размерной несогласованности. Затем весовые коэффициенты корректируются в соответствии со сценарием применения для достижения более сбалансированного распределения мощности. Математическая формулировка нормализованной модели OFDM-ISAC-P представлена в Уравнении 9:

Здесь w обозначает весовой коэффициент связи, Fc представляет собой оптимальное значение целевой функции, полученное в результате решения задачи оптимизации связи, pm обозначает мощность, выделяемую m-й поднесущей OFDM, gm представляет собой эквивалентный коэффициент усиления отношение сигнал/шум m-й коммуникационной поднесущей, P обозначает общую мощность передачи, Fs представляет собой оптимальное значение целевой функции, полученное в результате решения задачи оптимизации зондирования, а bm обозначает m-й элемент вектора зондирования B(v, β).
Весовой коэффициент связи w вместе с дополнительным весовым коэффициентом зондирования (1 − w) служит настраиваемым параметром проектирования, который отражает относительную значимость пропускной способности канала связи и PDP в различных сценариях применения28,29. В данном исследовании w рассматривается как параметр конфигурации на уровне сценария, а не как изменяющаяся во времени управляющая переменная. Соответственно, w фиксируется для каждого сценария моделирования (например, w = 0.5 для сбалансированной производительности связи и зондирования и w = 0.8 для работы с приоритетом связи). Затем улучшенный алгоритм BOA используется для оптимизации соответствующего распределения мощности. Таким образом, термин «динамический весовой коэффициент» относится к выбору различных весовых коэффициентов для разных рабочих сценариев, а не к реализации автоматической адаптации в режиме реального времени во время выполнения алгоритма.
Решение модели оптимизации мощности для интегрированной системы зондирования и связи на основе ортогонального частотного разделения с использованием улучшенного алгоритма оптимизации бабочек
После построения модели OFDM-ISAC-P для решения задачи оптимизации используется улучшенный алгоритм BOA. BOA имеет простую структуру, легко внедряется и обладает высокой способностью к глобальному поиску. Он имитирует кормовое поведение бабочек и использует механизм левитовского полета для эффективного исследования пространства поиска. Алгоритм требует настройки относительно небольшого количества параметров, что делает его подходящим для широкого спектра задач оптимизации. Однако традиционный BOA может преждевременно сходиться к локальным оптимумам, чувствителен к настройкам параметров и часто требует большого количества итераций при решении многомерных задач оптимизации30. Поэтому в данной работе усовершенствованы стратегия мутации и процедура оптимизации.
Для улучшения глобального поиска на начальном этапе оптимизации вводится стратегия DSP. В ходе первых итераций относительно низкая вероятность переключения способствует более активному исследованию пространства поиска. По мере увеличения числа итераций вероятность переключения постепенно корректируется для усиления локальной эксплуатации и улучшения сходимости. DSP определяется в Уравнении 10:

где α обозначает начальное значение DSP, t представляет номер текущей итерации, а β является коэффициентом управления. Эта адаптивная стратегия увеличивает разнообразие популяции на раннем этапе оптимизации, расширяет пространство поиска и снижает вероятность преждевременной сходимости к локальным оптимумам. По мере продвижения оптимизации поиск постепенно переключается с глобального исследования на локальную эксплуатацию.
Для дальнейшего повышения эффективности оптимизации в алгоритм включен метод DVGM. На начальных итерациях применяется относительно большая дисперсия мутации для расширения диапазона поиска и улучшения глобального исследования. По мере протекания оптимизации дисперсия мутации постепенно снижается для повышения точности локального поиска и облегчения сходимости. Максимальная дисперсия σmax и минимальная дисперсия σmin определены в Уравнении 11:

Улучшенные стратегии глобального и локального поиска обобщены в уравнении 12:

Здесь normrnd(0, σ(t)) обозначает гауссово случайное число со средним значением 0 и дисперсией σ(t), а fi представляет собой параметр аромата, который определяет поведение бабочек при поиске. После внедрения стратегий DSP и DVGM получается улучшенный алгоритм BOA, как показано на Рисунке 5.

Рисунок 5Рабочий процесс усовершенствованного алгоритма оптимизации бабочек (BOA). Процедура оптимизации включает инициализацию, генерацию популяции, оценку приспособленности, динамическую вероятность переключения (DSP), глобальный и локальный поиск, динамическую гауссову мутацию с переменной дисперсией (DVGM), обновление решения, оценку сходимости и вывод оптимального решения. Пожалуйста, нажмите здесь, чтобы просмотреть эту фигуру в увеличенном размере.
Рисунок 5 иллюстрирует две стратегии улучшения, внедренные в BOA: DSP и DVGM. Стратегия DSP адаптивно регулирует частоту переключения между глобальным и локальным поиском. На ранней стадии оптимизации алгоритм делает упор на глобальное исследование для повышения разнообразия популяции и увеличения вероятности обнаружения перспективных областей решения. По мере продвижения оптимизации поиск постепенно смещается в сторону локальной эксплуатации для ускорения сходимости. Стратегия DVGM вносит возмущения в текущую лучшую особь, тем самым расширяя возможности локального поиска, увеличивая разнообразие решений и снижая вероятность преждевременной сходимости к локальным оптимумам.
При применении улучшенного алгоритма BOA для решения модели OFDM-ISAC-P переменные оптимизации сначала кодируются, создается начальная популяция, и каждому индивидууму назначаются возможные стратегии распределения мощности. Общий рабочий процесс оптимизации показан на Рисунке 6.

Рисунок 6. Алгоритм работы улучшенного алгоритма оптимизации «Бабочка» (BOA) для решения модели распределения мощности в интегрированной системе связи и зондирования на базе мультиплексирования с ортогональным частотным разделением каналов (OFDM-ISAC). Процедура включает инициализацию параметров, кодирование векторов-кандидатов для распределения мощности, генерацию популяции, оценку приспособленности, итеративное обновление решения, оценку сходимости и вывод оптимальной стратегии распределения мощности. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Как показано на Рисунке 6, оптимизация начинается с применения стратегии DSP для проведения итерационной оптимизации. На ранних итерациях основной акцент делается на глобальном исследовании для изучения различных комбинаций распределения мощности в системе OFDM-ISAC. На более поздних итерациях поиск постепенно смещается в сторону локальной эксплуатации для уточнения высококачественных кандидатов в решения при соблюдении ограничения на общую мощность передачи. Впоследствии к текущему лучшему решению применяется стратегия DVGM для выхода из локальных оптимумов и поддержания разнообразия решений. На протяжении всего процесса оптимизации целевые функции связи и зондирования используются для оценки приспособленности каждого потенциального решения, при этом особи более высокого качества сохраняются для следующей итерации. Оптимизация продолжается до выполнения заданного критерия завершения. Итоговое оптимальное решение представляет собой стратегию распределения мощности, которая обеспечивает наилучший баланс между характеристиками связи и эффективностью обнаружения цели.
Параметры воспроизводимости
Для повышения воспроизводимости подробно указаны параметры симуляции, настройки инициализации, ограничения и критерии завершения. Все эксперименты были реализованы в Ubuntu 22.04 LTS с использованием Python 3.10, PyTorch 2.1, NumPy 1.26, CVXPY 1.4 и MATLAB R2023b (см. Таблицу материалов). Параметры симуляции и вычислительная среда, использованные в данной работе, обобщены в Таблице 1, в то время как подробная конфигурация параметров симуляции OFDM-ISAC представлена в Дополнительной таблице 1. В симуляциях использовались спектральная плотность мощности шума −174 dBm/Hz, 64 поднесущие, 16 OFDM-символов, несущая частота 4 GHz, шаг поднесущей 15 kHz, восемь передающих антенн и максимум 10 000 итераций. Коэффициент распределения мощности был инициализирован в диапазоне 0–1 с разрешением 0,01 и нормализован для соблюдения ограничения по общей мощности передачи. Вес связи, вес зондирования и точность многослойного априорного распределения были установлены на уровне 0,2, 0,3 и 0,5 соответственно. Ограничения оптимизации включали неотрицательное распределение мощности пользователей, ограниченную общую мощность передачи, ограниченную полосу пропускания и вероятность ложной тревоги порядка 10⁻3. Оптимизация завершалась либо при достижении максимального количества итераций, либо когда относительное изменение целевой функции между двумя последовательными итерациями опускалось ниже 10⁻4. Все эксперименты с рандомной инициализацией и возмущениями канала проводились с использованием фиксированных случайных чисел (seeds) и независимо повторялись 30 раз.
| Параметр | Значение | Описание |
| Несущая частота | 4 GHz | Несущая частота OFDM. |
| Схема кодирования/модуляции связи | Квадратурная фазовая манипуляция (QPSK) с канальным кодированием с малой плотностью проверок на четность (LDPC) | Метод генерации сигнала связи, используемый для воспроизведения передаваемой формы сигнала OFDM. |
| Вес связи | 0.2 | Вес, присвоенный задаче связи. |
| Вероятность ложной тревоги | 10⁻³ | Вероятность ложной тревоги, используемая в качестве ограничения при оптимизации зондирования. |
| Аппаратная платформа | Intel Core i7-13700 CPU, NVIDIA RTX 4070 GPU, 32 GB RAM | Вычислительная платформа, использованная для моделирования. |
| Эталонный набор данных для моделирования ISAC | Неприменимо. Все эталонные данные были получены с помощью численного моделирования. | Эталонный набор данных, использованный для оценки эффективности. |
| Максимальное количество итераций | 10,000 | Максимальное количество итераций оптимизации. |
| Точность многослойного априорного распределения | 0.5 | Параметр априорной информации, используемый при оптимизации. |
| Спектральная плотность мощности шума | −174 dBm/Hz | Спектральная плотность шума канала связи. |
| Количество OFDM-символов | 16 | Количество OFDM-символов в одном кадре. |
| Количество OFDM-поднесущих | 64 | Количество OFDM-поднесущих. |
| Операционная система | Ubuntu 22.04 LTS | Вычислительная среда. |
| Программное обеспечение для оптимизации | CVXPY 1.4 | Пакет для выпуклой оптимизации. |
| Разрешение распределения мощности | 0.01 | Разрешение инициализированных коэффициентов распределения мощности. |
| Язык программирования | Python 3.10 | Язык программирования, использованный для реализации моделирования. |
| Случайное число (seed) | Фиксированное | Фиксированное случайное число, использованное для всех серий моделирования. |
| Повторные запуски | 30 | Количество независимых запусков моделирования, выполненных для статистического анализа. |
| Библиотека для научных вычислений | NumPy 1.26 | Библиотека численных вычислений, использованная для проведения расчетов. |
| Вес зондирования | 0.3 | Вес, присвоенный задаче зондирования. |
| Платформа моделирования | MATLAB R2023b | Платформа, использованная для численного моделирования и визуализации. |
| Критерий остановки | Максимальное число итераций или относительное изменение целевой функции < 10⁻⁴ | Критерий завершения оптимизации. |
| Интервал между поднесущими | 15 kHz | Интервал между OFDM-поднесущими. |
| Библиотека тензорных вычислений | PyTorch 2.1 | Фреймворк для тензорных вычислений, использованный для матричных и тензорных операций. |
| Общее количество передающих антенн | 8 | Количество передающих антенн. |
Таблица 1: Параметры моделирования и вычислительная среда, использованные для моделирования интегрированного зондирования и связи (ISAC) на основе ортогонального частотного разделения каналов (OFDM).В таблице приведены параметры формы сигнала, параметры оптимизации, настройки моделирования, описание вычислительной среды, программного обеспечения, наборов данных и другие условия моделирования, использованные в ходе исследования.
Если не указано иное, процентные улучшения, представленные как «увеличение на X%», рассчитываются как относительный процентный прирост с использованием Уравнения 13:

Для показателей эффективности, основанных на вероятности или скорости, таких как PDP, абсолютные изменения эксплицитно указываются в процентных пунктах (pp). Например, увеличение с 46% до 63% указывается как +17 pp, тогда как +37,0% обозначает соответствующее относительное процентное увеличение, рассчитанное с использованием 46% в качестве базового значения.