Next Article in Journal
Estimating Global Anthropogenic CO2 Gridded Emissions Using a Data-Driven Stacked Random Forest Regression Model
Next Article in Special Issue
Deep Reinforcement Learning-Based Adaptive Modulation for Underwater Acoustic Communication with Outdated Channel State Information
Previous Article in Journal
Impact of Topographic Correction on PRISMA Sentinel 2 and Landsat 8 Images
Previous Article in Special Issue
Energy-Efficient Cooperative MIMO Formation for Underwater MI-Assisted Acoustic Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Analysis of Relay-Aided NOMA Optical Wireless Communication System in Underwater Turbulence Environment

School of Information and Communication Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 3894; https://doi.org/10.3390/rs14163894
Submission received: 3 July 2022 / Revised: 1 August 2022 / Accepted: 9 August 2022 / Published: 11 August 2022
(This article belongs to the Special Issue Underwater Communication and Networking)

Abstract

:
Non-orthogonal multiple access (NOMA) is a promising technology to improve spectrum utilization effectively for underwater optical wireless communications (UOWC). To exploit the benefits of NOMA in a turbulent environment, cooperative transmission has been introduced in the NOMA–UOWC network. The existing studies on NOMA suggest that relay selection and power optimization are the main factors affecting system performance. In this paper, a general NOMA node pairing method and two power optimization schemes for NOMA–UOWC are proposed, and both schemes are proven to be strictly quasi-convex. The two optimization schemes are solved by the BFGS algorithm and the particle swarm algorithm, respectively. The effectiveness of the proposed schemes are evaluated by our simulations, and the main factors affecting the relay-aided NOMA performance are derived.

Graphical Abstract

1. Introduction

Underwater optical wireless communication (UOWC) will be employed as a vital component of the marine internet of things (IoT) in the future due to its high bit rate, ultra-low latency, and high security [1]. Many researchers have carried out work on UOWC, and various theoretical and experimental studies have been conducted to investigate the behavior of optical beams underwater [2,3,4]. In the meantime, UOWC still confronts many challenges; e.g., the complex noise sources in the ocean environment and the severe channel fading caused by absorption, scattering, turbulence, and bubbles in the seawater [3,4]. On the other hand, most current underwater communication devices use light-emitting diodes (LEDs) as light sources or transmitters, for cost reasons, and the limited modulation bandwidth of LED limits the performance improvement of UOWC systems to some extent [5]. Non-orthogonal multiple access (NOMA) technology could be an attractive alternative to overcome these drawbacks.
There have been notable researches that demonstrated the ability to improve system spectrum utilization of NOMA in radio frequency (RF) networks and free-space optical networks [6,7]. As for underwater networks, most of the research on NOMA focuses on optimizing the receiving rate and improving the communication performance of the system through reasonable resource management, while researches on the exact underwater communication channel are investigated infrequently [8,9,10]. In [8], a practical NOMA scheme for a visible light communication (VLC) system on land was proposed, which improved spectral efficiency by 50% compared to orthogonal frequency division multiple access (OFDMA). An underwater optical wireless cellular network based upon NOMA was presented in [9], and the closed-form expressions of ergodic capacity for the near-end and the far-end receiving nodes in a single NOMA group were deduced. The influence of weak ocean turbulence on the performance of the underwater optical NOMA system was studied in [10], and an optimizing problem joining the user–subcarrier pairing with power allocation was formulated to maximize the system utility. Nevertheless, the fairness of the far-end receiving nodes in the NOMA system was neglected. In addition, the effectiveness and reliability of a terrestrial NOMA communication system can be appropriately balanced through flexible power control and interference suppression; however, there still exist many issues to be explored before this technology can be employed in the underwater wireless optical NOMA communication network in the complex underwater noise and channel environment.
Nevertheless, the superimposed transmission of signals has inevitably contributed to the increase of inter-user interference, which imposes restrictions on the performance of the NOMA communication system [11]. Motivated by the fact that the cooperative communication technology can improve transmission reliability of communication networks markedly by means of relay forwarding [12], the relay-aided NOMA scheme may be more consistent with practical applications. Nevertheless, the research on power allocation algorithms is mostly focused on the improvement of the FPA (fixed power allocation) method based on the characteristics of the channel to date, which distinctly lack flexibility in optimizing the performance of the NOMA–UOWC system.
Simultaneously, many studies in recent years have been implemented for relay-assisted UOWC. Most of the current research has been based on the scattering and absorption characteristics of the underwater channel [13,14,15,16], with hardly any consideration of background ocean noise [17]. Ref. [13] discusses the requirements and design challenges for relay-assisted UOWC. In [14], the attenuation and fading characteristics of UOWC in turbulent channels were systematically analyzed, and the performances of amplified forwarding (AF) and decoded forward (DF) multi-hop relay networks were investigated, respectively. The connectivity of the UOWC network based on the multi-hop transmission was investigated in [15,16]. An optimization strategy for relay deployment in the UOWC network was proposed in [17], and a resource allocation optimization algorithm was designed to uniformly manage the power allocation for relay nodes.
For the underwater optical network, there have been few researches on the combination of NOMA and cooperative communication technology. In [18], a cooperative NOMA system was proposed with the near-end node playing a role of relay node. However, this study ignores the direct transmission path between the source node and the receiver node 2, and it does not discuss the combined reception of this part of the optical signal. The above researches did not consider the combination of NOMA and the independent relay nodes, nor did they take into account the more realistic and complex underwater channel environment, such as solar radiation background noise, etc.
In this paper, the relay selection and power allocation technologies of a relay-aided NOMA–UOWC network are deeply explored, focusing on the combination of NOMA and cooperative communication technologies. In consideration of the influence of underwater solar radiation noise and turbulent channels, the cooperative transmission model of UOWC based on NOMA is established, and a user pairing strategy for joint relay selection is proposed. A global optimal power allocation algorithm (GOPA) and a stepwise sub-optimization power allocation algorithm (SSOPA) are designed for the relay-assisted NOMA communication network based on AF. The effectiveness of the proposed strategy and algorithm is verified by our simulations.

2. System Model

The deployment of nodes is constrained by the communication range in the seawater environment, and an overly long distance between communication nodes will lead to unstable optical communication links. In this paper, we consider the cooperative transmission of UOWC based on AF relay, since such a relaying technique is comparatively easy to implement and its computational complexity is relatively low. The diagram of the relay-aided NOMA–UOWC system is depicted in Figure 1. An underwater optical network is taken into account, in which a transmitter serves N receiving nodes simultaneously with the assistance of M relaying nodes.
According to the principle of NOMA, the user pairing and the relay selection need to be accomplished in real-time before the transmission [19]. Firstly, the channel gains [20,21] between the communication nodes can be calculated by:
h = α 2 attenuation   exp ( c d ) transmission   loss A de η tr η de cos θ 2 π d 2 ( 1 cos θ 0 ) geometric   loss ,
where c = ɑ + b is the attenuation coefficient of seawater, and ɑ and b are the absorption coefficient and the scattering coefficient of seawater, respectively. The optical properties for four typical seawaters are listed in Table 1 [1,22]. d indicates the transmission distance of the optical signal, Ade expresses the aperture area of the optical receiver, ηtr and ηde signify the transmitter optical efficiency and the receiver optical efficiency, respectively. θ is the tilt angle between the optical axis of the transmitter and the receiver, and θ0 represents the light divergence angle of the transmitter. α represents the attenuation of the optical signal caused by the turbulence. Here, the weak turbulent channel is taken into consideration, so the probability density function (PDF) of α can be expressed as [23,24]:
f ( α ) = 1 2 α ( π ln ( σ sci 2 + 1 ) ) / 2 exp ( 2 ( ln α μ s ) 2 ln ( σ sci 2 + 1 ) ) ,
where σ sci 2 indicates the scintillation coefficient, and lnα obeys a Gaussian distribution with a mean value of μs and variance of (1/4)ln ( σ sci 2 + 1).
Without loss of generality, the channel state information (CSI) between communication nodes satisfies h 1 > h 2 > > h n > > h N , and h n (n ∈ [1, N]) indicates the optical channel gain from the transmitter to the nth receiving node. It is assumed that the mth relay node is chosen as the cooperative relay by the pairing NOMA receiving nodes i and j, the superimposed signal transmitted by the transmitter to the relay node at the first time slot can be provided by [25]:
x t m = a i m P t m x i + a j m P t m x j ,
where P t m represents the transmitted power when the mth relay node is chosen, a i m and a j m are the power allocation coefficients of the paired users, and xi and xj represent the modulated signals of the ith node and the jth node. The signal received by the mth relay node can be written as:
y m = h t m ( a i m P t m x i + a j m P t m x j ) + n ¯ t m ,
where h t m (m ∈ [1, M]) is the channel gain between the transmitter and the mth relay node, and n ¯ t m signifies the noise between the transmitter and the mth relay node. Unlike terrestrial communications networks (TCN), there exist multiple noise sources in a marine environment. We mainly consider solar radiation noise, dark current noise, and shot noise as the main noises; therefore, the following expression can be obtained as [26,27]:
n ¯ = 2 B q η de A de π ( θ fov ) 2 Δ λ T F L sol solar   noise + 2 B q η de P t m shot   noise + 4 B k B T e φ F / R lr thermal   noise ,
where B signifies the system bandwidth, q is the elementary charge, Δλ indicates the optical filter bandwidth, TF represents the transmittance of the optical filter, Lsol expresses the solar irradiance, Te is the equivalent temperature, φF represents the system noise figure, Rlr is the load resistance, θfov indicates the field of view for the receiver, and kB denotes the Boltzmann constant.
Utilizing AF mentioned above, the signal forwarded by the mth relay node can be expressed as:
x m = δ η tr y m ,
where δ represents the amplification coefficient of the relay node. In terms of the principle of AF, the amplification coefficient, δ, can be obtained as follows [28]:
δ = P m η t η de ( h m P t m + n m ) ,
where Pm indicates the transmitted power of mth relaying.
As shown schematically in Figure 1, the receiving node receives the transmitted signals from the transmitter and the relay node in two-time slots under the line-of-sight (LOS) links. The maximal ratio combining (MRC) is employed at the receiver [29]. With synchronizing adopted, the received signal at the ith receiving node is:
y i = τ 1 ( h i x t m + n ¯ i ) + τ 2 ( h m , i x m + n ¯ m , i ) ,
where τ1, τ2 represent the combination factors, and hm,i indicates the channel gains between the ith receiving nodes and the mth relay node. The received signal at the jth receiving node can be calculated the same way. The stronger user restores the signal by performing success interference cancellation (SIC) to remove the signal of the weak user. The received signal-to-noise ratios obtained by the ith and the jth nodes can be respectively formulated as:
Γ i m = ( a i m P t m h i / n ¯ i ) 2 ( 1 + a i m P t m h i / n ¯ i ) 2 + a i 2 ( P t m h m / n ¯ t m ) 2 ( P m h m , i / n ¯ m , i ) 2 ( a j m P m h m , i P t m h m / n ¯ m , i n ¯ t m + P t m h m / n ¯ t m + P m h m , i / n ¯ m , i ) 2 ,
Γ j m = ( a j m P t m h j n ¯ j ) 2 + ( a j m P t m h m / n ¯ t m ) 2 ( P m h m , j / n ¯ m , j ) 2 ( P t m h m / n ¯ t m + P m h m , j / n ¯ m , j ) 2 .

3. User Pairing Scheme Jointed with Relay Selection

The first and most important issue in conducting the user pairing study is how many users access a NOMA group, which can implement optimal system performance. More users accessing the NOMA group implies more users sharing power resources. The current mainstream studies are mostly based on two-user paired NOMA, the multi-user pairing is realized by adding new user accesses to the two-user paired group. Hence, the performance of multi-user paired NOMA is highly dependent on the channel quality of the newly added users, according to the NOMA principle. Taking a three-user paired NOMA system as an example, the system performance is improved only when the channel quality of the newly added users is simultaneously better than all users in the original NOMA group. In addition, as the number of pairing users increase, the probability of miscommunication increases geometrically. Therefore, this study focuses on the situation in which two users are paired considering the reliability of the system.
Although the user pairing technique can improve the performance of NOMA systems, the existing user pairing scheme in TCN is not completely applicable to the UOWC network with the combination of NOMA and cooperative communication techniques. In this section, an algorithm of user pairing joint relay selection is proposed, which can maximize the edge user communication rate while ensuring the central user communication requirement based on comprehensively considering the CSIs of the transmitter to receiving NOMA nodes and the relay nodes to receiving nodes.
In accordance with the analysis of the system model above, the transmitter selects an optimal relay node to forward the pairing NOMA signal according to the CSIs of the central user i and the relay nodes. The selection of the relay node can be expressed by:
m = arg max r [ 1 , M ] | h r , i + h i | .
This algorithm preferentially ensures the communication requirement of the NOMA central user. Based on the principle of NOMA, the set of edge nodes, Ced, can be sorted by:
C ed = { j * :   i < j * N , B log 2 ( 1 + Γ j * m ) R th } ,
where Rth indicates the threshold of the data rate. Accordingly, the NOMA pairing edge user j is elected from Ced through (13):
j = arg max k C ed { log 2 ( 1 + Γ k m ) + log 2 ( 1 + Γ i m ) } .
It is obvious that the pairing NOMA group in the proposed scheme can obtain the maximum sum rate while ensuring a certain quality of service.

4. Sum Rate-Based Power Optimization

The deployment of UOWC nodes is one of the most significant current discussions in the construction of marine IoT. Since the complexity of secondary replacement for batteries is very high, energy efficiency is an important consideration for the UOWC network where the nodes are powered with finite energy. In this section, two optimization schemes have been proposed to maximize the system sum rate with limited energy consumption while supporting a certain communication requirement.

4.1. Fixed Power Allocation Method

One of the most commonly used power allocation algorithms in a practical NOMA network is the FPA method [30,31]. The fixed power allocation algorithm’s central mechanism is to assign a fixed power allocation ratio to all NOMA pairing users. Each user is allocated a proportionate share of the transmitted power.
In this case, we set the power allocation ratio as ξ. The optimization problem can be expressed as follows:
max R total
subject   to   a j m = ξ a i m
a j m + ξ a i m = 1
m = 1 M ( P t m + P m ) P t o t a l
P t m > 0 , P m > 0
It is noted that the average allocation with the lowest algorithmic complexity has been employed for the bandwidth allocation between the NOMA groups. For the sake of providing a fair comparison, the results obtained by using different multiple access schemes are evaluated with the same bandwidth allocation method.
Hence, the total sum rates of the system in Equation (14c) can be calculated by the following equation:
R total = m = 1 M B M ( log 2 ( 1 + Γ j m ) + log 2 ( 1 + Γ i m ) ) .
The algorithm is designed to be simple and easy to implement. The fairness of the NOMA system can be adjusted by adjusting the power allocation ratios. The specific algorithm implementation can be referred to in Algorithm 1. However, the algorithm ignores the NOMA group’s users’ channel conditions and cannot traverse all possible power allocation schemes. The step ratio configuration is critical to the optimization solution. The complexity of the system algorithm will increase if the step factor is too small, and the optimal solution of the algorithm cannot be found if the step factor is too large. While solving the optimization problem of (14d), the FPA algorithm, on the other hand, cannot solve the power distribution between the relay node and the source node directly. Hence, we employ the average distribution method to allocate the power between the relay node and the source node.
Algorithm 1 Fixed Power Allocation Algorithm
Initialize Power allocation ratio ξ1 = 0, step factor Δξ, the total transmitted power of system Ptotal, the number of relay nodes M.
Output: The power allocation ratioζ, the sum rates Rtotal.
1: Initialize smax = [1/Δξ], R m * = 0 , Rtotal = 0, ζ[].
2:  for m = 1→M do
3:    Determine NOMA pairing nodes i and j according to (11)~(13)
4:    for s = 1→smax do
5:      ξ#= ξs + Δξ
6:      Compute a i m , a j m according to (14b), (14c)
7:      Compute Γ i m , Γ j m according to (9) and (10)
8:      R sum m = B / M ( log 2 ( 1 + Γ i m ) + B log 2 ( 1 + Γ j m ) )
9:     if  R sum m > = R m * then
10:        R m * = R sum m , ξm = ξ#
11:      else
12:       s = s + 1
13:     end if
14:    end for
15:        R total = R total + R m *
16:        ζ(m)= ξm
17:  end for

4.2. Global Optimal Power Allocation

In this subsection, the optimization objective is set to maximize the data rate of the NOMA group while the total transmitted power is a certain value. As to the design of the sum rates-based power allocation scheme, we have referred to the optimal power control scheme proposed in Ref. [32]. It is assumed that the total transmitted power included the transmitted powers of the transmitter and relay nodes, which also implies that the transmitted power for each NOMA group is identical. The sum rate-based power optimization problem can be formulated as:
max R total
subject   to   0   <   a i m < a j m  
  a i m + a j m 1
m = 1 M ( P t m + P m ) P t o t a l
P t m > 0 , P m > 0
It can be seen from (14d) that the power allocation optimization model of multi-user NOMA joint relay selection is a complex non-deterministic polynomial problem. In general, there is no algorithm with polynomial complexity to solve the optimal solution of this type of problem directly. Nevertheless, in the NOMA system, the pairing users share the same time–frequency resource through superposition in power domains. Additionally, since orthogonal multiple access (OMA) is applied for transmission between NOMA groups, the NOMA group consumes a certain power at each data transmission [33]. The above problem model can be equivalently simplified as follows:
max R m
subject   to   0 < a i m < a j m
  a i m + a j m 1
P t m + P m P total M
P t m > 0 , P m > 0
where Ptotal takes a certain constant value. It can be easily proved that the necessary condition for the objective function to obtain the optimal solution is a i m + a j m = 1 and P t m + P m = P total / M . The optimization problem is demonstrated to be strictly quasi-concave with respect to a i m and P t m . Therefore, the Lagrangian function of this optimization problem can be expressed as:
L g ( P t m , a i m , l 1 , l 2 , l 3 , l 4 ) = R m + l 1 ( 2 a i m 1 ) + l 2 ( 2 P t m P total M ) l 3 a i m l 4 P t m ,
where l1, l2, l3, and l4 represent the Lagrange multipliers, and the Karush–Kuhn–Tucker (KKT) conditions are shown below:
L g P t m = R m P t m + 2 l 2 l 4
L g a i m = R m P t m + 2 l 1 l 3
l 1 ( 2 a i m 1 ) = 0
l 2 ( 2 P t m P total M ) = 0
l 3 a i m = 0 , l 4 P t m = 0
It can be calculated from (18a) and (18c) that l1 = 0, and meanwhile, l3 and l4 are denoted by (18d), (18e) to be 0. The optimization problem is clearly shown to be a complex nonlinear optimization problem. The other parameter values are supposed to be solved by the BFGS quasi-Newton algorithm [34].
As the BFGS algorithm is used to solve the optimization problem, the iterative equation is as follows:
P ( s + 1 ) = P ( s ) + ϖ ( s ) d ( s ) ,
where P = ( P t m , a i m , l 2 ) , ϖ(s) is the step-size factor which can be defined as (20), and d(s) is the search direction, which can be calculated by (21):
ϖ ( s ) = arg min ϖ L g ( P ( s ) + ϖ d ( s ) ) ,
d ( s ) = ( H ( s ) ) 1 ( g ( s ) ) 1 .
The gradient vector of the objective function Rm concerning a i m and P t m , g, can be expressed as:
g = F ( P ) = [ L g P t m , L g a i m , l 2 ( 2 P t m P total M ) ] T ,
where P = ( P t m , a i m , l 2 ) , and Hessian matrix is as follows:
H = F 2 ( P ) = [ 2 L g ( P t m ) 2 2 L g P t m a i m l 2 2 L g a i m P t m 2 L g ( a i m ) 2 l 2 l 2 l 2 2 P t m P total M ] .
Thus, the iteration can be obtained by (24).
F ( P ) F ( P ( s + 1 ) ) + F 2 ( P ( s + 1 ) ) ( P P ( s + 1 ) ) .
When the value of P is set to be P(s), the real-valued vector calculated from the gradient vector is g s = F ( P ( s ) ) , and the formulation can be made out as follows:
F ( P ( s + 1 ) ) F ( P ( s ) ) F 2 ( P ( s + 1 ) ) ( P ( s + 1 ) P ( s ) ) .
The BFGS algorithm takes the second-order approximation when Taylor expansion of the objective function is performed and the gradient operators on both ends of the expansion equation are applied, which simplifies the operation difficulty. For the convenience of subsequent calculations, the Hessian matrix is also approximated in the iterative process, and the correction matrix Δ B ( s ) is added, which means B(s)H(s), and the iterative formula can be obtained by:
B ( s + 1 ) = B ( s ) + Δ B ( s ) ,
where ΔB(s) is the correction matrix, which can be formulated as:
Δ B ( s ) = ( p ( s ) ) ( p ( s ) ) T ( p ( s ) ) T ( q ( s ) ) B ( s ) ( q ( s ) ) ( q ( s ) ) T B ( s ) ( q ( s ) ) T B ( s ) ( q ( s ) ) ,
where q ( s ) = P ( s + 1 ) P ( s ) and p ( s ) = g ( s + 1 ) g ( s ) are the marker vector symbols.
The optimization algorithm of global optimal power allocation is designed as shown in Algorithm 2.
Algorithm 2 Global Optimal Power Allocation Optimization
Initialize Target rate Rth, maximum iterations smax, the precision threshold τ, the total transmitted power of system Ptotal, the number of relay nodes M.
Output: The   allocation   of   transmitted   power   and   the   Lagrange   multiplier ( a i m , P t m , l 2 ) T .
1 : Set   the   initial   power   allocation   coefficient   and   transmitted   power ( a i m * , P t m * , l 2 ) T = ( 1 / 2 , P total / 2 M , 0 ) T .
2: for s = 1→smax do
3:   Determine the search direction
4:      d(s) = −(B(s))−1 × g(s)
5:    Compute the step factor ϖ(s) according to (20)
6:    Compute the mark vector
7:     q(s) = l(s)d(s), P(s+1) = P(s)+ q(s)
6:    Compute ||g(s+1)||2
7:     if ||g(s+1)||2 < τ then
8:        break
9 :     else   compute Γ i m , Γ j m according to (9) and (10)
10 :   R i m = B / M ( log 2 ( 1 + Γ i m ) ) , R j m = B / M ( log 2 ( 1 + Γ j m ) )
11 :   if R i m > R j m & R j m > = R th then
12:        Compute p(s) = g(s+1)g(s) according to (25)
13:        Compute B(s+1) = B(s) + ΔB(s) according to (27)
14:        Let s = s + 1 and return to step 3
15:       else break
16:       end if
17:     end if
18: end for

4.3. Stepwise Sub-Optimization Power Allocation

The optimal solution for the system can be obtained by our proposed GOPA algorithm, proposed in Section 4.1, through traversing the overall system. Nonetheless, the complexity of the system algorithm will be raised to a great degree as the number of nodes accessing the network increases and hence, we further propose a step-by-step sub-optimization algorithm (SSOPA) with a lower computing complexity.
In this problem formulation, the power allocation coefficients in the NOMA group are preferentially determined, with the relay transmitted power assumed to fulfil the QoS (quality of service) of the edge user in the NOMA group. Under this condition, the power allocation issues can be optimized between the transmitter and relay nodes. In particular, the CSIs between the transmitter and the relays, and between the relays and receiving nodes, are comprehensively considered. According to (1), the power allocation coefficients can be conducted as follows:
{ a i m = 1 1 + ( ( h i h m , i ) / ( h j h m , j ) ) a j m = ( h i h m , i ) / ( h j h m , j ) 1 + ( ( h i h m , i ) / ( h j h m , j ) ) .
The problem model can be simplified as follows:
max R m
subject   to   P t m + P m P total / M
P t m > 0 , P m > 0
This optimization model is still an optimizing problem of nonlinear functions. To reduce the complexity of the algorithm, the particle swarm algorithm (PSA) is employed to solve this problem, which is easier to be implemented. As a bionic algorithm, the PSA has natural advantages in solving global optimization problems. Compared with the genetic algorithm, the PSA can memorize all solutions, only needs to update according to the internal speed, and its convergence speed is even faster [35].
The fitness function in this problem is defined as follows:
fit ( P t m , P m ) = R m .
The PSA’s solution procedure is as follows:
(1) The particle swarm is initialized randomly. Each particle carries two types of information: position information c and velocity information V, with the position information corresponding to the value taken in the fit function, which is the allocated transmitted power [ P t m , P m ] . We define the maximum number of iterations smax, the particle swarm size K, and the particle’s maximum speed Vmax.
(2) The initial adaptation value of the particles in each population is calculated according to (30) to obtain the optimal solution pBest in that population.
(3) Searching the initial optimal solution gBest in all swarms.
(4) Update particle velocity and position. During the tth iteration of the kth particle, the value of ck and the speed value of vk, are updated by:
v k ( t ) = χ v k ( t 1 ) + ϑ 1 r 1 ( p b e s t k c k ( t 1 ) ) + ϑ 2 r 2 ( g b e s t ( t 1 ) c k ( t 1 ) ) ,
c k ( t ) = c k ( t 1 ) + c k ( t 1 ) ,
where χ represents the inertia weight. When χ becomes larger, the global optimizing ability of the model is enhanced, and the local optimizing ability is weakened. c1 and c2 denote acceleration coefficients that reflect the intensities of particle self-learning and group learning.
(5) Return to step 2.
The optimizing algorithm for the global optimal power allocation is designed as shown in Algorithm 3.
Algorithm 3 Stepwise Sub-Optimization Power Allocation
Step 1 Compute the power allocation coefficient of pairing nodes according to (26)
Step 2 Optimize the power allocation between transmitter and relay node
Initialize The   number   of   particle   swarm   K ,   the   maximum   iteration   s max ,   the   initial   adaptive   function   value   R m * ←0.
Output: The   allocation   of   transmitted   power ( P t m , P m ) T .
1:  for k ←1 to K
2:   Initialize the initial source node transmission power factor ck and the velocity factor vk.
3 :   compute   Γ i m , k , Γ j m , k according to (9) and (10)
4 :     R m ( k ) = B / M ( log 2 ( 1 + Γ i m , k ) + log 2 ( 1 + Γ j m , k ) )
5 :         if     R m ( k )   >   R m *  then
6:     pbest(k)c(k), gbest   c ( k ) , R m * R m ( k )
7:     else break
8:    end if
9:  end for
10:    while not stop
11:     fors ← 1 to smax
12:     Update the source node transmission power factor c(k) and the velocity factor v(k) according to (31) and (32)
13:     if fit(c(k)) > fit(pbest(k)) then
14:       pbest(k)c(k)
15:     if fit(pbest(k)) > fit(gbest) then
16:       gbestpbest(k)
17:     end for
18 :           [ P t m , P m ] g b e s t
19:   end while

4.4. Analysis of Algorithm Complexity

In this paper, three algorithms are proposed to solve the problem model in different ways. This section analyzes the three algorithms in terms of time complexity and problem size, in order to measure the time and resources consumed by the three algorithms during the calculation process.
The FPA method can complete the algorithm’s solution with a simple nested loop, and the sum rate can be solved in the nested problem using simple multiplication and addition operations. The matrix multiplication operation is required in GOPA. Although we optimize the search direction and step factor solution process by constructing an approximate Hesse matrix rather than solving the Hesse process, the problem scale remains relatively large at this time. The algorithm complexity of SSOPA is proportional to the number of iterations and the initialized population size, and the problem size in a single iteration is small. The problem size of a single mathematical calculation is set to 1 here. To summarize, the complexity statistics for the three algorithms are shown in Table 2. Some of the variables employed in this paper are summarized in Table 3.

5. Simulation and Analysis

In this section, we will present comprehensive simulation studies to verify the effectiveness of the proposed schemes. In this paper, we mainly consider the scenarios for the line-of-sight communication links. The configuration scheme of the transmitter in this paper refers to the quasi-omnidirectional prismatic LED array module that is proposed in [36]. The parameter settings are summarized in Table 4.

5.1. Simulation Scenarios for NOMA–UOWC Network

To study the performance of the relay-aided NOMA–UOWC system, three scenarios are designed as shown in Figure 2. For the convenience of simulation calculation, a rectangular coordinate is established with the source node as the origin. In scenarios 1 and 2, the specific positions of receiving nodes are given as (8.69, 2.33), (11.59, −3.11), (14.94, −1.31), and (17.93,1.57), severally, and the Euclidean distances between the source node and receiving nodes are 9 m, 12 m, 15 m, and 18 m, separately. The tilt angles are set up as 15°, −15°, −5°, and 5°, respectively. In scenario 3, the tilt angle of the receiving nodes remains unchanged, but the Euclidean distances between the source node and receiving nodes increase to 18 m, 24 m, 30 m, and 36 m, separately. During the deployment phase of the relay nodes, we refer to the relay deployment method for single-hop networks of Ref. [17], which has been mentioned in the introduction, and temporarily disregard the geometric loss of the optical communication link. The specific positions of relay nodes are given as (3.94, 0.69) and (6.89, −1.22) in scenarios 1 and 2, while the tilt angles are set up as 10° and −10°. As with the deployment of receiving nodes, in scenario 3 the tilt angle of relay nodes remains unchanged, and the Euclidean distances between the source node and relay nodes increase to 8 m and 14 m, respectively. The specific position of each node is shown in Figure 2.
In this subsection, our main work is to verify the rationality of the proposed system model, and the simulations are performed temporarily, considering only a single NOMA group transmission. The transmission power only occupied half of the total system power. Prior to further system performance analysis, the effectiveness of NOMA for underwater wireless optical communication systems is to be validated. In this case, a comparative simulation is performed to investigate the performances of NOMA-based UOWC and OMA-based UOWC with the same available bandwidth resources, and the receiving nodes in the transmission model are opted for node 1 and node 3 of scenario 1 in Figure 2. The most current 2 × 2 NOMA network is applied in this simulation. The detailed results can be seen in Figure 3. We can easily conclude from Figure 3 that the sum rate of the system increases when the power transmitted increases. It is notable that a NOMA-based system can support almost 1.4 times more effective communication sum rates than an OMA-based systems for the UOWC network when the transmitted power is 1000 mW, due to sharing the same time–frequency resource.
One more important issue that we consider for the research of NOMA–UOWC network is the user pairing. To demonstrate the system performance variation of a three-user paired NOMA relative to a traditional two-user paired NOMA with the same bandwidth and power resources, a simple simulation is conducted with the results provided in Figure 4. Referring to the theoretical analysis in Section 3, the total transmitted power in the two simulations are assumed to be the same, 1000 mW. The transmission distances of user 1 and user 2 are 4 m and 9 m, respectively. The blue line is a control group presented in the simulation, showing the system sum rates of the traditional two-user paired NOMA. In other words, the resources allocated to the third user are zero, so the system sum rates are constant. The orange line shows the system sum rates for three-user paired NOMA. Compared to the control group, it includes the original two-user and the newly increased user 3. As the Euclidean distance between the newly increased user 3 and the transmitter increases, the system sum rates decrease. Furthermore, we can obtain the following conclusions. As the number of paired users increases in a single NOMA group, the system sum rates do not definitely increase for a constant power resource, and the sum rates of three-user paired NOMA are better than that of two-user paired NOMA only when user 3 is the nearest user to the transmitter. Hence, we mainly consider the 2 × 2 NOMA networks for the UOWC network based on cooperative transmission in this paper.
In addition, to highlight the performance of the relay-aided NOMA scheme (denoted as ICNOMA), ICNOMA is compared with two different schemes, which are summarized as follows. (a) ONOMA: a NOMA system with no relay or cooperative communication adopted, seen in Figure 5a. (b) CNOMA: it is assumed that the near-end node has both signal reception and signal forwarding functions in the NOMA system shown in Figure 5b.

5.2. Performance Evaluations

To verify the effectiveness of the proposed relay-aided NOMA under different user pairing conditions, relative to the NOMA scheme proposed in Ref. [37] and the conventional NOMA scheme, a series of simulations were performed, and the detailed results are presented in Figure 6. In Figure 6, 1&2, 1&3, and 1&4 represent user 1 and user 2, 3, and 4 pairings, respectively. The system model of each NOMA scheme can be referred to in Figure 5. To minimize the interference of other components, the simulation environment is set as scenario 1, and the relay node 1 is set to be the relay node. Since the relay nodes have been identified, the FPA algorithm is performed in both NOMA schemes. The results for ONOMA, CNOMA, and ICNOMA are shown in Figure 6. It is evident that the introduction of cooperative communication can significantly improve the performance of a NOMA system. Therein, an independent relay can further enhance the system sum rate; for instance, the sum rates for ICNOMA with different users paired are enhanced by 0.6 bps/Hz to 0.9 bps/Hz compared to the CNOMA when the transmitted power is the same, 1000 mW.
For the purpose of assessing the proposed GOPA and SSOPA algorithms, the sum rates of the NOMA–UOWC network for three scenarios are illustrated in Table 5, where the FPA represents a fixed power allocation scheme, in which the transmitted power for the transmitter and relay nodes are evenly distributed. The attenuation coefficient of seawater is selected as 0.056 m−1 to analyze the impact of solar noise on the UOWC network performance. Therein, the GOPA and SSOPA algorithms perform the prominent ability in improving the system communication capability in scenario 1. Significantly, these two proposed schemes have the absolute advantage over the traditional FPA method, due to the fact that the proposed user pairing scheme jointed with relay selection not only optimizes the cooperative transmission, but also equilibrates the channel fading caused by the path loss and the ocean turbulence. It is noted that in scenario 3, with the increases in depth of nodes and the transmission distance, the GOPA algorithm has a better capability of balancing the path loss and the solar noise. It can be concluded that the GOPA algorithm presents the ability to extend the communication range of the UOWC network.
The results of the relay selection and power allocation for the proposed network in the clear seawater of scenario 1 are summarized in Table 6. To analyze the impact of cooperative transmission on the UOWC network performance, a control group without a cooperative relay was introduced in the performance evaluation. It was noted that the optimal results can be obtained by the GOPA algorithm; meanwhile, the SSOPA algorithm achieves a sub-optimal solution. It can be observed that the GOPA algorithm increases the sum rate by 0.034 bps/Hz in comparison with the SSOPA algorithm. To further evaluate the impact of cooperative transmission on UOWC network performance, the None-relay (a control group without a cooperative relay) scheme was performed [18]. The results of the four schemes show that the relay-assisted approaches prove to be distinct superior compared to the None-relay mode, which verifies the effectiveness of the proposed relay-aided NOMA network. In addition, the GOPA algorithm and SSOPA algorithm proposed in this paper retain the effectiveness in improving the system sum-rates in a conventional NOMA system without cooperative communication.
For the purpose of assessing the proposed scheme in different ocean environments, the comparisons of GOPA, SSOPA, and FPA algorithms are demonstrated in Figure 7, where the scintillation indexes reflect the turbulence intensity. We adopt the FPA algorithm as a contrast to assess the proposed GOPA and SSOPA algorithms. Among these three methods, the GOPA algorithm has the most obvious sum-rate increase compared to the others. Therein, the increment shows obvious changes with different σ sci 2 varying from 0.1 to 1.0. This is due to the fact that the power allocation efficiency has disparate computation complexity in different turbulent environments. It can be summarized from Figure 7a that the sum rate of the GOPA algorithm is nearly the same as that of the SSOPA algorithm, with σ sci 2 = 0.1 . Meanwhile, the sum rate of the GOPA algorithm exceeds those of SSOPA and FPA algorithms by far, with σ sci 2 = 1.0 . We can conclude that the sum rate of GOPA is clearly improved in the three scenarios, and the SSOPA algorithm shows relatively better performance in the weaker turbulence environment.
To further study the impact of solar noise on the results of a relay-aided NOMA–UOWC network, three typical ocean seawaters, which can be referred to in Table 1, were employed to create the depth simulation. It can be observed that the receiving sum rates of NOMA nodes change greatly with the node depth varying from 0 m to 50 m. Significantly, the receiving data rates in the harbor water are far too low compared to the other three conditions. Therefore, the comparisons of these three scenarios are shown in Figure 8. Similar to the results detailed above, the GOPA algorithm obtains the optimized performance in the sum rate of the NOMA–UOWC network under all conditions. However, the gap between GOPA and SSOPA algorithms has narrowed almost to vanishing point in the clear ocean water. In contrast, the GOPA algorithm shows a considerable advantage to the SSOPA and FPA algorithms in the coastal water, with the sum rate of the proposed network showing a maximum at the depth of around 20 m. In the meantime, the sum rates of the proposed network in clear ocean water and pure ocean water increase approximately in line with the depth within 50 m. According to the formulas derived above, the solar noise is positively associated with the depth of node, from which we can deduce that the depth of node has a crucial impact on the performance of the UOWC network. However, the impact of solar noise in the coastal water is relatively weak, which is caused by the fact that the path loss and the noise are balanced at a depth of around 20 m. From a comprehensive perspective, the proposed algorithms of relay selection and power allocation optimization have outstanding capabilities of not only balancing the path loss and solar noise, but also equilibrating the channel fading between the transmitter and the relays, and between the relays and the receiving nodes.

6. Conclusions

A relay-aided NOMA scheme for UOWC networks is proposed in this paper, in which a user pairing method of joint relay selection is suggested. Due to solar noise, two power allocation-optimizing algorithms based on the maximization of the sum rates, i.e., the GOPA and SSOPA algorithms, are designed. Both schemes are proven to be quasi-convex optimizing problems, whose solutions are obtained by the BFGS Newton algorithm and the particle swarm algorithm, respectively. The proposed algorithms are evaluated in AF-based cooperative transmission UOWC networks, and the simulation results verify that the proposed optimization scheme performs more effectively in improving the system communication capability. Furthermore, solar radiation noise has been shown to be one of the main factors affecting NOMA–UOWC network performance.

Author Contributions

Conceptualization, Y.L., H.Y.; formal analysis, Y.L.; methodology, Y.L. and H.Y.; software, Y.L., L.J., X.J. and J.W.; validation, Y.L., X.J. and J.W.; writing—original draft, Y.L.; writing—review and editing, H.Y. and Y.L.; visualization, Y.L.; supervision, H.Y.; project administration, H.Y.; funding acquisition, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the National Natural Science Foundation of China (NSFC) under Grant 61871418 and 61801079, China.

Data Availability Statement

The data presented in this paper are available after contacting the corresponding author.

Acknowledgments

The authors would like to thank the anonymous reviewers for their careful reading and valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaushal, H.; Kaddoum, G. Underwater Optical Wireless Communication. IEEE Access 2016, 4, 1518–1547. [Google Scholar] [CrossRef]
  2. Schirripa Spagnolo, G.; Cozzella, L.; Leccese, F. Underwater Optical Wireless Communications: Overview. Sensors 2020, 20, 2261. [Google Scholar] [CrossRef] [PubMed]
  3. Xue, K.; Ma, R.; Wang, D.; Shen, M. Optical Classification of the Remote Sensing Reflectance and Its Application in Deriving the Specific Phytoplankton Absorption in Optically Complex Lakes. Remote Sens. 2019, 11, 184. [Google Scholar] [CrossRef]
  4. Uudeberg, K.; Aavaste, A.; Kõks, K.-L.; Ansper, A.; Uusõue, M.; Kangro, K.; Ansko, I.; Ligi, M.; Toming, K.; Reinart, A. Optical Water Type Guided Approach to Estimate Optical Water Quality Parameters. Remote Sens. 2020, 12, 931. [Google Scholar] [CrossRef]
  5. Verniers, K.; Nobby, S. Impact of LED Optical Bandwidth Limitation on the Performance of Power Switched Signaling. In Proceedings of the 2020 12th International Symposium on Communication Systems, Networks and Digital Signal Processing (CSNDSP), Porto, Portugal, 20–22 July 2020; pp. 1–6. [Google Scholar]
  6. Elouafadi, R.; Raiss El Fenni, M.; Benjillali, M. NOMA Clustering for Improved Multicast IoT Schemes. J. Sens. Act. Netw. 2022, 11, 26. [Google Scholar] [CrossRef]
  7. Islam, S.M.R.; Avazov, N.; Dobre, O.A.; Kwak, K.-S. Power-Domain Non-Orthogonal Multiple Access (NOMA) in 5G Systems: Potentials and Challenges. IEEE Commun. Surv. Tutor. 2017, 19, 721–742. [Google Scholar] [CrossRef]
  8. Lin, B.; Tang, X.; Zhou, Z.; Lin, C.; Ghassemlooy, Z. Experimental demonstration of SCMA for visible light communications. Opt. Commun. 2018, 419, 36–40. [Google Scholar] [CrossRef]
  9. Jiang, R.; Sun, C.; Tang, X.; Zhang, L.; Wang, H.; Zhang, A. Joint User-Subcarrier Pairing and Power Allocation for Uplink ACO-OFDM-NOMA Underwater Visible Light Communication Systems. J. Lightwave Technol. 2021, 39, 1997–2007. [Google Scholar] [CrossRef]
  10. Geldard, C.; Thompson, J.; Popoola, W.O. A Study of Non-Orthogonal Multiple Access in Underwater Visible Light Communication Systems. In Proceedings of the 2018 IEEE 87th Vehicular Technology Conference (VTC Spring), Porto, Portugal, 3–6 June 2018; pp. 1–6. [Google Scholar]
  11. Akhtar, M.W.; Hassan, S.A.; Jung, H. On the Symbol Error Probability of STBC-NOMA with Timing Offsets and Imperfect Successive Interference Cancellation. Electronics 2021, 10, 1386. [Google Scholar] [CrossRef]
  12. Li, Y.; Yu, Z.; Zhou, H.; Tao, J. To Relay or not to Relay: Open Distance and Optimal Deployment for Linear Underwater Acoustic Networks. IEEE Trans. Commun. 2018, 66, 3797–3808. [Google Scholar] [CrossRef]
  13. Saeed, N.; Celik, A.; Al-Naffouri, T.Y.; Alouini, M.S. Underwater optical wireless communications, networking, and localization: A survey. Ad Hoc Netw. 2019, 94, 101935. [Google Scholar] [CrossRef]
  14. Jamali, M.V.; Chizari, A.; Salehi, J.A. Performance Analysis of Multi-Hop Underwater Wireless Optical Communication Systems. IEEE Photonics Technol. Lett. 2017, 29, 462–465. [Google Scholar] [CrossRef]
  15. Saeed, N.; Celik, A.; Alouini, M.S.; Al-Naffouri, T.Y. Performance Analysis of Connectivity and Localization in Multi-Hop Underwater Optical Wireless Sensor Networks. IEEE. Trans. Mob. Comput. 2019, 18, 2604–2615. [Google Scholar] [CrossRef]
  16. Inoue, Y.; Kodama, T.; Kimura, T. Global Optimization of Relay Placement for Seafloor Optical Wireless Networks. IEEE Trans. Wirel. Comm. 2021, 20, 1801–1815. [Google Scholar] [CrossRef]
  17. Xing, F.; Yin, H.; Shen, Z.; Leung, V.C.M. Joint Relay Assignment and Power Allocation for Multiuser Multirelay Networks over Underwater Wireless Optical Channels. IEEE Internet Things J. 2020, 7, 9688–9701. [Google Scholar] [CrossRef]
  18. Palitharathna, K.W.S.; Suraweera, H.A.; Godaliyadda, R.I.; Herath, V.R.; Thompson, J.S. Average Rate Analysis of Cooperative NOMA Aided Underwater Optical Wireless Systems. IEEE Open J. Comm. Soc. 2021, 2, 2292–2310. [Google Scholar] [CrossRef]
  19. Palitharathna, K.W.S.; Suraweera, H.A.; Godaliyadda, R.I.; Herath, V.R.; Ding, Z. Impact of Receiver Orientation on Full-Duplex Relay Aided NOMA Underwater Optical Wireless Systems. In Proceedings of the 2020 IEEE International Conference on Communications (ICC), Virtual Conference, 7–11 June 2020; pp. 1–7. [Google Scholar]
  20. Miroshnikova, N.E.; Petruchin, G.S.; Titovec, P.A. Estimation of the Effect of Dispersion on the Communication Range in a Wireless Underwater Optical Channel. In Proceedings of the 2021 Systems of Signal Synchronization, Generating and Processing in Telecommunications (SYNCHROINFO), Svetlogorsk, Russia, 30 June–2 July 2021; pp. 1–4. [Google Scholar]
  21. Aravind, J.V.; Kumar, S.; Prince, S. Mathematical Modelling of Underwater Wireless Optical Channel. In Proceedings of the 2018 International Conference on Communication and Signal Processing (ICCSP), Chennai, India, 3–5 April 2018; pp. 776–780. [Google Scholar]
  22. Sticklus, J.; Hoeher, P.A.; Röttgers, R. Optical Underwater Communication: The Potential of Using Converted Green LEDs in Coastal Waters. IEEE J. Ocean Eng. 2019, 44, 535–547. [Google Scholar] [CrossRef]
  23. Wang, K.; Gong, C.; Zou, D.; Xu, Z. Turbulence Channel Modeling and Non-Parametric Estimation for Optical Wireless Scattering Communication. J. Lightwave Technol. 2017, 35, 2746–2756. [Google Scholar] [CrossRef]
  24. Bi, C.; Qing, C.; Wu, P.; Jin, X.; Liu, Q.; Qian, X.; Zhu, W.; Weng, N. Optical Turbulence Profile in Marine Environment with Artificial Neural Network Model. Remote Sens. 2022, 14, 2267. [Google Scholar] [CrossRef]
  25. Zhang, L.; Wang, Z.; Wei, Z.; Chen, C.; Wei, G.; Fu, H.Y.; Dong, Y. Towards a 20 Gbps multi-user bubble turbulent NOMA UOWC system with green and blue polarization multiplexing. Opt. Express 2020, 28, 31796–31807. [Google Scholar] [CrossRef]
  26. Zhang, J.; Gao, G.; Wang, B.; Guan, X.; Yin, L.; Chen, J.; Luo, B. Background noise resistant underwater wireless optical communication using Faraday atomic line laser and filter. J. Lightwave Technol. 2021, 40, 63–73. [Google Scholar] [CrossRef]
  27. Patterson, A.M.; Spence, J.H.; Fischer, R.W. Evaluation of underwater noise from vessels and marine activities. In Proceedings of the IEEE/OES Acoustics in Underwater Geosciences Symposium, Rio de Janeiro, Brazil, 25–27 July 2013; pp. 1–9. [Google Scholar]
  28. Essiambre, R.-J.; Kramer, G.; Winzer, P.J.; Foschini, G.J.; Goebel, B. Capacity Limits of Optical Fiber Networks. J. Lightwave Technol. 2010, 28, 662–701. [Google Scholar] [CrossRef]
  29. Balti, E.; Guizani, M.; Hamdaoui, B. Hybrid Rayleigh and Double-Weibull over impaired RF/FSO system with outdated CSI. In Proceedings of the 2017 IEEE International Conference on Communications (ICC), Paris, France, 21–25 May 2017; pp. 1–6. [Google Scholar]
  30. Benjebbovu, A.; Li, A.; Saito, Y.; Kishiyama, Y.; Harada, A.; Nakamura, T. System-level performance of downlink NOMA for future LTE enhancements. In Proceedings of the 2013 IEEE Globecom Workshops (GC Wkshps), Atlanta, GA, USA, 9–13 December 2013; pp. 66–70. [Google Scholar]
  31. Zhang, H.; Wang, B.; Chen, C.; Cheng, X.; Li, H. Resource Allocation in UAV-NOMA Communication Systems Based on Proportional Fairness. J. Commun. Inf. Netw. 2020, 5, 111–120. [Google Scholar] [CrossRef]
  32. Ghosh, J.; Haci, H.; Kumar, N.; Al-Utaibi, K.A.; Sait, S.M.; So-In, C. A Novel Channel Model and Optimal Power Control Schemes for Mobile mmWave Two-Tier Networks. IEEE Access 2022, 10, 54445–54458. [Google Scholar] [CrossRef]
  33. Gagnon, F.; Dahman, G.; Poitau, G. Extending the ITU-R P.530 Deep-Fading Outage Probability Results to SIMO-MRC and MIMO-MRC Line-of-Sight Systems. IEEE Wirel. Commun. Lett. 2018, 7, 1086–1089. [Google Scholar] [CrossRef]
  34. Aryan, M.; Alejandro, R. RES: Regularized Stochastic BFGS Algorithm. IEEE Trans. Signal Process. 2014, 62, 6089–6104. [Google Scholar]
  35. Grubisic, S.; Carpes, W.P.; Bastos, J.P.A.; Santos, G. Association of a PSO Optimizer with a Quasi-3D Ray-Tracing Propagation Model for Mono and Multi-Criterion Antenna Positioning in Indoor Environments. IEEE Trans. Magn. 2013, 49, 1645–1648. [Google Scholar] [CrossRef]
  36. Tong, Z.; Yang, X.; Chen, X.; Zhang, H.; Zhang, Y.; Zou, H.; Zhao, L.; Xu, J. Quasi-omnidirectional transmitter for underwater wireless optical communication systems using a prismatic array of three high-power blue LED modules. Opt. Express 2021, 29, 20262–20274. [Google Scholar] [CrossRef]
  37. Elamassie, M.; Bariah, L.; Uysal, M.; Muhaidat, S.; Sofotasios, P.C. Capacity Analysis of NOMA-Enabled Underwater VLC Networks. IEEE Access 2021, 9, 153305–153315. [Google Scholar] [CrossRef]
Figure 1. Architecture of the relay-aided NOMA–UOWC system.
Figure 1. Architecture of the relay-aided NOMA–UOWC system.
Remotesensing 14 03894 g001
Figure 2. Diagram of simulation scenarios.
Figure 2. Diagram of simulation scenarios.
Remotesensing 14 03894 g002
Figure 3. Sum rates versus transmitted power for OMA-based UOWC and NOMA-based UOWC.
Figure 3. Sum rates versus transmitted power for OMA-based UOWC and NOMA-based UOWC.
Remotesensing 14 03894 g003
Figure 4. Sum rates of relay-aided UOWC system with two-user paired and three-user paired NOMA.
Figure 4. Sum rates of relay-aided UOWC system with two-user paired and three-user paired NOMA.
Remotesensing 14 03894 g004
Figure 5. Diagram of three kinds of NOMA schemes. (a) Ordinary NOMA, (b) cooperative NOMA without independent relay nodes, (c) independent-relay-aided NOMA.
Figure 5. Diagram of three kinds of NOMA schemes. (a) Ordinary NOMA, (b) cooperative NOMA without independent relay nodes, (c) independent-relay-aided NOMA.
Remotesensing 14 03894 g005
Figure 6. Sum rates versus transmitted power of NOMA-based UOWC system.
Figure 6. Sum rates versus transmitted power of NOMA-based UOWC system.
Remotesensing 14 03894 g006
Figure 7. Sum rates versus transmitted power of ICNOMA–UOWC system for (a) σ sci 2 = 0.1 , (b) σ sci 2 = 0.4 , (c )   σ sci 2 = 1.0 .
Figure 7. Sum rates versus transmitted power of ICNOMA–UOWC system for (a) σ sci 2 = 0.1 , (b) σ sci 2 = 0.4 , (c )   σ sci 2 = 1.0 .
Remotesensing 14 03894 g007
Figure 8. Sum rates versus depth of transmitter for ICNOMA system for: (a) clear ocean water, (b) pure ocean water, and (c) coastal water.
Figure 8. Sum rates versus depth of transmitter for ICNOMA system for: (a) clear ocean water, (b) pure ocean water, and (c) coastal water.
Remotesensing 14 03894 g008
Table 1. Typical values of attenuation coefficient for seawater.
Table 1. Typical values of attenuation coefficient for seawater.
Water Typeɑ (m−1)b (m−1)c (m−1)
Pure seawater0.0530.0030.056
Clean seawater0.1140.0370.151
Coastal seawater0.1790.2190.398
Harbor water0.2951.8752.170
Table 2. Comparison of Algorithm Complexity.
Table 2. Comparison of Algorithm Complexity.
AlgorithmTime ComplexityProblem Size
FPAO(smax)107M
GOPA O ( s max 2 ) 1570M/τ
SSOPAO(smax)278MK
Table 3. Variables commonly used in this paper.
Table 3. Variables commonly used in this paper.
xTransmitted signalyReceived signal
PPowerhChannel gain
nNumber of the receiving nodes, n ∈ [1, N]mNumber of the relay nodes, m ∈ [1, M]
aNOMA power allocation coefficientΖPower allocation ratio
n ¯ Noise powersmaxIterations of algorithm
τStep-size factorkNumber of the particle swarm, k ∈ [1, K]
δAmplification coefficientdTransmission distance
cAttenuation coefficientθBeam divergence angle
θ0Inclination angleαOptical fading amplitude
Table 4. Simulation parameters.
Table 4. Simulation parameters.
ParametersValues
System bandwidth B32 MHz
Total transmitted power Ptotal2000 mW
Divergence angle of the transmitter θ030°
Aperture area of the optical receiver Ade0.01 m2
P-I modulation conversion coefficient ηtr0.9 A/W
Responsibility of the photodetector ηde0.9 W/A
Table 5. Comparisons of three scenarios for ICNOMA–UOWC.
Table 5. Comparisons of three scenarios for ICNOMA–UOWC.
Relay SelectionRelay SelectionRtotal (bps/Hz)
Scenario 1GOPAC1 (i = 1, j = 4)C2(i = 2, j = 3)7.92
SSOPAC1 (i = 1, j = 4)C2(i = 2, j = 3)7.62
FPAC2 (i = 1, j = 3)C2(i = 2, j = 4)6.10
Scenario 2GOPAC1 (i = 1, j = 4)C2(i = 2, j = 3)11.05
SSOPAC1 (i = 1, j = 4)C2(i = 2, j = 3)10.87
FPAC1 (i = 1, j = 3)C2(i = 2, j = 4)8.65
Scenario 3GOPAC1 (i = 1, j = 4)C1 (i = 2, j = 3)2.15
SSOPAC1 (i = 1, j = 4)C1 (i = 2, j = 3)1.50
FPAC1 (i = 1, j = 3)C1 (i = 2, j = 4)0.37
Table 6. Comparisons of three schemes for ICNOMA–UOWC.
Table 6. Comparisons of three schemes for ICNOMA–UOWC.
Relay Selection P t m   ( mW ) Pm (mW) a i m a j m Rtotal (bps/Hz)
GOPAC1 (i = 1, j = 4)421.52578.480.2060.7967.92
C2 (i = 2, j = 3)361.20638.800.3770.623
SSOPAC1 (i = 1, j = 4)445.73554.270.2110.7897.62
C2 (i = 2, j = 3)356.85643.150.3840.616
FPAC1 (i = 1, j = 4)500.00500.000.1910.8096.10
C2 (i = 2, j = 3)500.00500.000.4010.599
N-relay
with GOPA
(i = 1, j = 4)1000.000.000.0510.9491.08
(i = 2, j = 3)1000.000.000.1130.887
N-relay
with SSOPA
(i = 1, j = 4)1000.000.000.0450.9550.97
(i = 2, j = 3)1000.000.000.1010.899
N-relay
with FPA
(i = 1, j = 4)1000.000.000.0340.9660.74
(i = 2, j = 3)1000.000.000.0870.913
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liang, Y.; Yin, H.; Jing, L.; Ji, X.; Wang, J. Performance Analysis of Relay-Aided NOMA Optical Wireless Communication System in Underwater Turbulence Environment. Remote Sens. 2022, 14, 3894. https://doi.org/10.3390/rs14163894

AMA Style

Liang Y, Yin H, Jing L, Ji X, Wang J. Performance Analysis of Relay-Aided NOMA Optical Wireless Communication System in Underwater Turbulence Environment. Remote Sensing. 2022; 14(16):3894. https://doi.org/10.3390/rs14163894

Chicago/Turabian Style

Liang, Yanjun, Hongxi Yin, Lianyou Jing, Xiuyang Ji, and Jianying Wang. 2022. "Performance Analysis of Relay-Aided NOMA Optical Wireless Communication System in Underwater Turbulence Environment" Remote Sensing 14, no. 16: 3894. https://doi.org/10.3390/rs14163894

APA Style

Liang, Y., Yin, H., Jing, L., Ji, X., & Wang, J. (2022). Performance Analysis of Relay-Aided NOMA Optical Wireless Communication System in Underwater Turbulence Environment. Remote Sensing, 14(16), 3894. https://doi.org/10.3390/rs14163894

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop