This study involved only mathematical modeling, algorithm development, and numerical simulations for an OFDM-ISAC system. No human participants, animals, biological specimens, or personally identifiable data were involved. Therefore, institutional ethics approval, including Institutional Review Board (IRB) or Institutional Animal Care and Use Committee (IACUC) approval, was not required.
Design of power allocation model for the Orthogonal Frequency Division Multiplexing-Based Integrated Sensing and Communication System
In ISAC systems, traditional equal PAMs often distribute transmission power evenly across each channel and antenna during simultaneous communication and sensing without considering their different performance requirements. This may allocate excessive power to inefficient transmission links and fail to satisfy the heterogeneous requirements of communication and sensing services, thereby reducing overall power utilization efficiency12,13,14. Therefore, this study considers the dual functional objectives of the system by jointly optimizing communication channel capacity and PDP. A multi-objective optimization framework is constructed to develop a power allocation strategy for the OFDM-based ISAC system, referred to as the OFDM-ISAC-P model. OFDM is an efficient multicarrier modulation technique consisting of a transmitting end and a receiving end. The OFDM process is illustrated in Figure 1.

Figure 1. Signal-processing workflow of the Orthogonal Frequency Division Multiplexing (OFDM) waveform used in the Integrated Sensing and Communication (ISAC) system. The workflow includes data encoding, interleaving, symbol mapping, serial-to-parallel conversion, inverse fast Fourier transform (IFFT), cyclic-prefix (CP) insertion, radio-frequency (RF) transmission, wireless channel propagation, RF reception, CP removal, fast Fourier transform (FFT), channel equalization, demapping, deinterleaving, decoding, and data recovery. Please click here to view a larger version of this figure.
In Figure 1, OFDM signals serve as the shared waveform for communication and sensing because of their high spectrum utilization and strong resistance to multipath interference, thereby enabling unified waveform transmission for both functions15,16,17. Neglecting the spectral shaping effect of the window function, the transmitted OFDM waveform is constructed using the carrier frequency, number of subcarriers, subcarrier complex weights for power allocation, communication signal encoding, and subcarrier spacing. The system bandwidth is determined by the subcarrier spacing and the number of subcarriers18. The communication signal generation process is specified as follows. Binary information bits are first encoded using a low-density parity-check (LDPC) coding scheme to improve transmission reliability and are subsequently mapped to complex-valued symbols using quadrature phase shift keying (QPSK) modulation. The modulated symbols are allocated to OFDM subcarriers according to the optimized power allocation strategy. Following inverse fast Fourier transform (IFFT) processing and cyclic prefix (CP) insertion, the time-domain OFDM communication waveform is generated and transmitted together with the sensing signal. The same modulation and coding configuration is maintained throughout all simulations. The ISAC system configuration used in this study is shown in Figure 2.

Figure 2. System configuration of the Orthogonal Frequency Division Multiplexing (OFDM)-based Integrated Sensing and Communication (ISAC) system. The diagram illustrates the transmitting (Tx) array, receiving (Rx) array, direct line-of-sight (LoS) communication path, target-echo sensing path, ground-reflection path, sky-reflection path, and representative obstacles within the simulated propagation environment. Please click here to view a larger version of this figure.
The system consists of a transmitting antenna array and a receiving antenna array. Signal transmission occurs through two main paths. The first is the Line-of-Sight (LoS) path, in which the signal is transmitted directly from the transmitting array to the receiving array. The second is the reflection path, in which a portion of the signal is reflected from the ground and reaches the receiving array, whereas another portion propagates through a sky-reflection path before reaching the receiving array19,20,21. The system is a LoS-dominant multi-antenna OFDM-ISAC system in which the communication and sensing receivers may adopt different antenna configurations. Although Figure 2 illustrates possible propagation components, including the direct LoS path and reflection-related paths, the mathematical model and simulations consider only the dominant LoS communication link and target echo sensing link. Ground and sky reflections are treated as background propagation effects rather than independently optimized non-LoS (NLoS) channels. The channel response matrix is defined under the LoS-dominant assumption, with multipath effects represented through fading, Doppler shift, and noise terms. Based on this channel response matrix, together with the transmitted signal and channel noise model, the received communication signal is expressed in Equation 1.

Here, nc,m represents Gaussian white noise, Sm denotes the transmitted signal, and Hc,m represents the communication channel response matrix.
In the sensing function of the ISAC system, the signal transmission path consists of two stages. The first stage is transmission from the transmitter to the target, and the second is the echo path, in which the reflected signal propagates from the target back to the receiver through the sensing channel22,23,24. For the m-th subcarrier, the sensing channel response matrix characterizes signal propagation and is expressed as Equation 2:

Here, bs represents the sensing echo-path fading coefficient, fs,m denotes the reflected Doppler frequency shift, and τs represents the propagation delay of the reflection path. Based on the sensing channel response matrix, the received sensing signal is expressed as Equation 3:

Here, ns,m represents Gaussian white noise in the sensing channel.
In wireless communication systems, efficient spectrum planning and resource allocation are essential for improving overall system performance25. Accordingly, this study uses communication channel capacity as the primary design criterion. An optimization model is established to maximize channel capacity by allocating the total transmission power across the OFDM subcarriers according to the proposed resource allocation strategy, as expressed in Equation 4.
The communication optimization problem is formulated as follows (Equation 4):

subject to

Here, σc2 represents the noise variance of the communication Gaussian channel, P denotes the total communication transmission power, Pc represents the optimal communication power allocation vector, am denotes the weight assigned to the m-th subcarrier, and Nc represents the number of subcarriers.
The optimal solution obtained using the WA is expressed as Equation 5:

Here,
denotes the optimal power allocated to the m-th communication subchannel, gm represents the channel-to-noise gain ratio of the m-th communication subchannel, and μc denotes the water level of the WA.
According to the WA illustrated in Figure 3, subchannels with higher channel gains, corresponding to smaller inverse channel gains, receive greater transmission power under the selected water level. This enables the system to exploit favorable channel conditions while satisfying the total power constraint and improving overall communication capacity. Consequently, subchannels with better channel conditions receive more power, whereas those with poorer channel conditions receive less. Through this adaptive power allocation strategy, the WA efficiently utilizes channel resources to improve communication performance.

Figure 3. Principle of the Water-Filling Algorithm (WA) used for communication power allocation. The horizontal axis represents the OFDM subchannels, and the vertical axis represents the inverse channel gain and allocated transmission power. The water level determines the transmit power assigned to each subchannel under the total transmission-power constraint. Please click here to view a larger version of this figure.
For the sensing objective, an appropriate power allocation strategy is essential when the false alarm probability is fixed. An optimized power allocation strategy enhances sensing performance by improving target detection accuracy. The sensing process requires hypothesis testing under the target-present hypothesis (H1) and the target-absent hypothesis (H0). Under these hypotheses, the received signal for each subcarrier during the k-th pulse is expressed as Equation 6:

Here, y denotes the sensing observation vector, A represents the diagonal matrix diag(a0, a1, …, aNc−1), D(v) denotes the diagonal matrix containing the target Doppler frequency shift, β represents the target reflection coefficient vector, and ns denotes Gaussian white noise in the sensing channel.
The probability density functions under H0 and H1 are then evaluated. For a specified false alarm probability, the PDP is calculated based on these two hypotheses. According to likelihood ratio detection theory, a generalized maximum likelihood ratio detector is subsequently derived to maximize the PDP.
The generalized maximum likelihood ratio detector is expressed as Equation 7:

Here, lG and LG denote the generalized likelihood ratio detection statistic and likelihood ratio function, respectively, yH represents the Hermitian transpose of y, and γ0 denotes the detection threshold. During the detection process, the threshold is used to determine whether the target-present hypothesis (H1) or the target-absent hypothesis (H0) is accepted. The relationship between the sensing detection threshold and the target detection probability is illustrated in Figure 4.

Figure 4. Relationship between the sensing detection threshold and target detection probability. The probability density functions under the target-absent hypothesis and target-present hypothesis are shown together with the detection threshold used for hypothesis testing. The shaded regions correspond to the decision regions associated with the false alarm probability and target detection probability. Please click here to view a larger version of this figure.
As shown in Figure 4, the PDP is directly related to the non-central complex chi-square distribution. Under a fixed false alarm probability (i.e., the maximum allowable false detection probability), the non-centrality parameter is positively correlated with the PDP. Consequently, increasing the non-centrality parameter improves the probability of correctly detecting the target signal. Based on this relationship, the total sensing transmission power is treated as the power constraint, whereas maximizing the non-centrality parameter is selected as the optimization objective. The resulting constrained optimization problem is used to derive the optimal sensing power allocation strategy and is expressed as Equation 8:

where α denotes the optimal sensing transmit weight vector obtained under the power constraint and also represents the sensing subcarrier transmit weight vector, B(v, β) denotes the sensing-information vector, and P represents the total sensing transmission power.
The optimization indicates that the optimal sensing power allocation strategy does not distribute power across multiple subchannels. Instead, all available sensing transmission power is allocated to the subchannel identified through channel evaluation as providing the highest sensing performance, thereby maximizing single-channel target detection performance.
In the OFDM-ISAC system, power allocation is the core optimization task, and the optimal strategy must satisfy the design requirements of both communication and sensing. Specifically, the optimization simultaneously considers communication performance metrics, including bit error rate and channel capacity, together with sensing performance metrics, including PDP and parameter estimation accuracy. Overall system performance is improved through coordinated resource allocation. However, directly combining the optimal communication and sensing power allocation strategies using conventional weighted-sum optimization has important limitations. Because the different physical dimensions and relative contributions of the two performance metrics are not adequately considered, the resulting optimization may not accurately reflect the requirements of integrated communication-sensing systems26,27.
Therefore, this study proposes an integrated design framework based on normalized weighting. First, the communication and sensing performance metrics are normalized to eliminate the effects of dimensional inconsistency. The weighting coefficients are then adjusted according to the application scenario to achieve a more balanced power allocation. The mathematical formulation of the normalized OFDM-ISAC-P model is given in Equation 9:

Here, w denotes the communication weighting factor, Fc represents the optimal objective value obtained from the communication optimization problem, pm denotes the power allocated to the m-th OFDM subcarrier, gm represents the equivalent channel-to-noise gain of the m-th communication subcarrier, P denotes the total transmission power, Fs represents the optimal objective value obtained from the sensing optimization problem, and bm denotes the m-th element of the sensing vector B(v, β).
The communication weighting factor w, together with its complementary sensing weighting factor (1 − w), serves as a tunable design parameter that reflects the relative importance of communication channel capacity and PDP under different application scenarios28,29. In this study, w is treated as a scenario-level configuration parameter rather than a time-varying control variable. Accordingly, w is fixed for each simulation scenario (e.g., w = 0.5 for balanced communication and sensing performance and w = 0.8 for communication-priority operation). The improved BOA is then used to optimize the corresponding power allocation. Therefore, the term “dynamic weighting factor” refers to selecting different weighting factors for different operating scenarios rather than implementing real-time automatic adaptation during algorithm execution.
Solution of the Power Optimization Model for the Orthogonal Frequency Division Multiplexing-Based Integrated Sensing and Communication System Using the Improved Butterfly Optimization Algorithm
After constructing the OFDM-ISAC-P model, the improved BOA is used to solve the optimization problem. The BOA has a simple structure, is easy to implement, and possesses strong global search capability. It mimics butterfly foraging behavior and employs a Lévy flight mechanism to efficiently explore the search space. The algorithm requires relatively few parameter settings, making it suitable for a wide range of optimization problems. However, the conventional BOA may converge prematurely to local optima, is sensitive to parameter settings, and often requires many iterations when solving high-dimensional optimization problems30. Therefore, the mutation strategy and optimization procedure are enhanced in this study.
To improve global exploration during the early optimization stage, a DSP strategy is introduced. During the initial iterations, a relatively small switching probability enhances exploration of the search space. As the number of iterations increases, the switching probability is gradually adjusted to strengthen local exploitation and improve convergence. The DSP is defined in Equation 10:

where α denotes the initial value of the DSP, t represents the current iteration number, and β is the control coefficient. This adaptive strategy increases population diversity during the early optimization stage, broadens the search space, and reduces the likelihood of premature convergence to local optima. As the optimization progresses, the search gradually shifts from global exploration to local exploitation.
To further improve optimization performance, DVGM is incorporated into the algorithm. A relatively large mutation variance is applied during the early iterations to expand the search range and improve global exploration. As the optimization proceeds, the mutation variance is gradually reduced to enhance local search accuracy and facilitate convergence. The maximum variance σmax and minimum variance σmin are defined in Equation 11:

The enhanced global and local search strategies are summarized in Equation 12:

Here, normrnd(0, σ(t)) denotes a Gaussian random number with a mean of 0 and variance σ(t), and fi represents the fragrance parameter that controls the butterfly search behavior. After incorporating the DSP and DVGM strategies, the improved BOA is obtained, as illustrated in Figure 5.

Figure 5. Workflow of the improved Butterfly Optimization Algorithm (BOA). The optimization procedure includes initialization, population generation, fitness evaluation, Dynamic Switching Probability (DSP), global and local searches, Dynamic Variance Gaussian Mutation (DVGM), solution updating, convergence assessment, and optimal solution output. Please click here to view a larger version of this figure.
Figure 5 illustrates the two improvement strategies incorporated into the BOA: DSP and DVGM. The DSP strategy adaptively adjusts the switching frequency between global and local search. During the early optimization stage, the algorithm emphasizes global exploration to increase population diversity and improve the probability of identifying promising solution regions. As the optimization progresses, the search gradually shifts toward local exploitation to accelerate convergence. The DVGM strategy perturbs the current best individual, thereby enhancing local search capability, increasing solution diversity, and reducing the likelihood of premature convergence to local optima.
When the improved BOA is applied to solve the OFDM-ISAC-P model, the optimization variables are first encoded, an initial population is generated, and candidate power allocation strategies are assigned to each individual. The overall optimization workflow is illustrated in Figure 6.

Figure 6. Workflow of the improved Butterfly Optimization Algorithm (BOA) for solving the Orthogonal Frequency Division Multiplexing-based Integrated Sensing and Communication (OFDM-ISAC) power allocation model. The procedure includes parameter initialization, encoding of candidate power-allocation vectors, population generation, fitness evaluation, iterative solution updating, convergence assessment, and output of the optimal power-allocation strategy. Please click here to view a larger version of this figure.
As shown in Figure 6, the optimization begins by applying the DSP strategy to perform iterative optimization. During the early iterations, global exploration is emphasized to investigate diverse power allocation combinations for the OFDM-ISAC system. In the later iterations, the search gradually shifts toward local exploitation to refine high-quality candidate solutions while satisfying the total transmission power constraint. The DVGM strategy is subsequently applied to the current best solution to escape local optima and maintain solution diversity. Throughout the optimization process, the communication and sensing objective functions are used to evaluate the fitness of each candidate solution, and higher-quality individuals are retained for the next iteration. The optimization continues until the predefined termination criterion is satisfied. The final optimal solution represents the power allocation strategy that provides the best balance between communication performance and target detection performance.
Reproducibility settings
To improve reproducibility, the simulation parameters, initialization settings, constraints, and termination criteria are specified in detail. All experiments were implemented on Ubuntu 22.04 LTS using Python 3.10, PyTorch 2.1, NumPy 1.26, CVXPY 1.4, and MATLAB R2023b (refer to Table of Materials). The simulation parameters and computational environment used in this study are summarized in Table 1, whereas the detailed OFDM-ISAC simulation parameter configuration is provided in Supplementary Table 1. The simulations used a noise power spectral density of −174 dBm/Hz, 64 subcarriers, 16 OFDM symbols, a carrier frequency of 4 GHz, a subcarrier spacing of 15 kHz, eight transmit antennas, and a maximum of 10,000 iterations. The power allocation ratio was initialized within the range of 0–1 with a resolution of 0.01 and normalized to satisfy the total transmission power constraint. The communication weight, sensing weight, and multi-layer prior accuracy were set to 0.2, 0.3, and 0.5, respectively. The optimization constraints included non-negative user power allocation, limited total transmission power, limited bandwidth, and a false alarm probability on the order of 10⁻3. The optimization terminated when either the maximum number of iterations was reached or the relative change in the objective function between two consecutive iterations fell below 10⁻4. All random initialization and channel perturbation experiments used fixed random seeds and were independently repeated 30 times.
| Parameter | Value | Description |
| Carrier frequency | 4 GHz | OFDM carrier frequency. |
| Communication encoding/modulation scheme | Quadrature Phase Shift Keying (QPSK) modulation with Low-Density Parity-Check (LDPC) channel coding | Communication signal generation method used to reproduce the transmitted OFDM waveform. |
| Communication weight | 0.2 | Weight assigned to the communication objective. |
| False alarm probability | 10⁻³ | False alarm probability used as the sensing optimization constraint. |
| Hardware platform | Intel Core i7-13700 CPU, NVIDIA RTX 4070 GPU, 32 GB RAM | Computational platform used for the simulations. |
| ISAC simulation benchmark dataset | Not applicable. All benchmark data were generated through numerical simulations. | Benchmark dataset used for performance evaluation. |
| Maximum number of iterations | 10,000 | Maximum number of optimization iterations. |
| Multi-layer prior accuracy | 0.5 | Prior-information parameter used during optimization. |
| Noise power spectral density | −174 dBm/Hz | Communication channel noise spectral density. |
| Number of OFDM symbols | 16 | Number of OFDM symbols per frame. |
| Number of OFDM subcarriers | 64 | Number of OFDM subcarriers. |
| Operating system | Ubuntu 22.04 LTS | Computational environment. |
| Optimization software | CVXPY 1.4 | Convex optimization package. |
| Power allocation resolution | 0.01 | Resolution of the initialized power allocation ratios. |
| Programming language | Python 3.10 | Programming language used for simulation implementation. |
| Random seed | Fixed | Fixed random seed used for all simulations. |
| Repeated runs | 30 | Number of independent simulation runs performed for statistical analysis. |
| Scientific computing library | NumPy 1.26 | Numerical computing library used for numerical computations. |
| Sensing weight | 0.3 | Weight assigned to the sensing objective. |
| Simulation platform | MATLAB R2023b | Platform used for numerical simulation and visualization. |
| Stopping criterion | Maximum iterations or relative objective change < 10⁻⁴ | Optimization termination criterion. |
| Subcarrier spacing | 15 kHz | OFDM subcarrier spacing. |
| Tensor computation library | PyTorch 2.1 | Tensor computation framework used for matrix and tensor computations. |
| Total transmit antennas | 8 | Number of transmitting antennas. |
Table 1: Simulation parameters and computational environment used for the Orthogonal Frequency Division Multiplexing (OFDM)-based Integrated Sensing and Communication (ISAC) simulations. The table summarizes the waveform parameters, optimization parameters, simulation settings, computational environment, software, datasets, and other simulation conditions used throughout the study.
Unless otherwise specified, percentage improvements reported as “increased by X%” are calculated as relative percentage gains using Equation 13:

For probability- or rate-based performance metrics, such as the PDP, absolute changes are explicitly reported as percentage-point (pp) increases. For example, an increase from 46% to 63% is reported as +17 pp, whereas +37.0% denotes the corresponding relative percentage increase calculated using 46% as the baseline.