Abstract

In this paper, we investigate the transmit power allocation problem to minimize the average packet error rate at the access point in the cluster flight spacecraft network, which adopts the CSMA/CA channel access mechanism. First, the node mobility, nodal distance distribution, and probabilistic adjacency matrix were formulated for cluster flight spacecraft network based on twin-satellite mode. Then, the optimization-theoretic model described the optimized transmit power allocation strategy and its implementation algorithm was proposed. And the problem of minimizing the packet error rate of the cluster flight spacecraft network system can be converted into maximizing the expectation of the binary probabilistic adjacency matrix, i.e., maximizing the sum of the nondiagonal elements in the probabilistic adjacency matrix. Due to discreteness of nodal distance distribution, Monte Carlo method was applied to solve the transmit power allocation problem. Yet importantly, the influence of node transmit power on the QoS performance of cluster flight spacecraft network was simulated and analyzed under the assumption of finite overall network transmit power and low traffic load. Finally, the results show that the pocket error rate increases with the provided traffic load, but the pocket error rate hardly changes with the same traffic load in different sequential time slots of any orbital hyperperiod or in the same time slot of different orbital hyperperiods, and by maximizing the sum of the nondiagonal elements in the probabilistic adjacency matrix, the pocket error rate minimum is achieved for a given total network transmit power at any time slot for cluster flight spacecraft network.

1. Introduction

In recent years, fractionated spacecraft with cluster flight model has become a hot topic in the field of distributed space network, due to its advantages of flexibility, rapid response, low cost, strong scalability, and long lifetime. The previous work has made a contribution to earth observation and space exploration [13]. Fractionated spacecraft distributes the functionality of a traditional large monolithic spacecraft into a number of heterogeneous modules. Each module can be regarded as a node through wireless communication, and the nodes construct the cluster flight spacecraft network (CFSN). Cluster flight spacecraft require mutual cooperation between nodes to realize information exchange, navigation communication, and power sharing. These spacecraft constitute a virtual satellite platform with information exchange structure. In addition, like other distributed space systems, the cluster flight spacecraft is a resource-sharing and energy-limited system. Therefore, how to allocate nodal transmit power efficiently and optimize the performance of the cluster flight spacecraft is an important issue needs to be solved [46].

For wireless communication systems, including wireless sensor networks and radar networks, research on optimizing systematic performance by allocating nodal power efficiently has always been a hotspot [7, 8]. For example, in an end-to-end MIMO multihops wireless network with outage probability limited, the power allocation method was studied and the analytic solution of optimal power allocation with the minimum total transmit power of the system was obtained in [9]. Aiming at the problem that it is difficult to obtain an analytic solution for the optimal power allocation in amplification and forwarding (AF) relay selection of cooperative communication system, artificial neural network was adopted to obtain an efficient solution from the target of minimum bit error rate in [10]. For the cooperative network of spectrum sharing, by solving the convex optimization problem, the optimal power allocation strategy with minimum energy consumption under the requirements of QoS was obtained in [11]. For the symbol programming problem of target detection in distributed radar sensor network, the closed-form expression of optimal power allocation by establishing the optimal linear unbiased estimation model was obtained in [12]. In [13], based on the received signal interference plus noise ratio, Markov chain was adopted to obtain the dynamic power control method with the minimum packet error rate (PER).

To the best of our knowledge, there are few research reports on the optimal power allocation in CFSN. Only the solution of the minimum spanning tree by constructing the spatial-temporal network topology to improve energy efficiency was obtained for the CFSN in [14]. In fact, due to the fact that heterogeneous modules have to meet miscellaneous requirements given different tasks, the CFSN faces the requirements of multi-QoS, including the requirements of different delay and bit error [3]. In general, the resource allocation with multi-QoS, including power, bandwidth, and CPU, is an NP-hard problem in the static network [15, 16]. However, because of the high-speed flight of modules, the topology of the CFSN is highly dynamic and nodes are randomly connected. Therefore, it is foreseeable that the resource allocation with QoS is more complicated in CFSN.

Many researchers have studied the problem of power distribution with multi-QoS in static networks. Among them, based on the resource allocation mode with QoS under the constraint of resource [17], the resource allocation problem with multiresources multi-QoS and single-resource multi-QoS by adopting polynomial concave optimal control method was solved in [15]. By introducing a truncated based on slope, an approximated method of concave optimal control was obtained and the resource allocation problem with multiresources multi-QoS was solved efficiently [16]. In addition, as for radar tracking under the constraints of radar bandwidth, processing time, and transmit power, the method of [16] was adopted to solve the optimization problem of radar tracking and speed accuracy efficiently in [18, 19]. In recent years, with the development of large-scale MIMO technology, the optimization of transmit power allocation has also received widespread attention. For example, the issue of the minimum total transmit power of the system from the perspective of satisfying QoS constraints was studied in [20]. An efficient energy allocation algorithm was studied based on binary search for 5G carrier aggregation scenario in [21]. Under the conditions of given QoS requirements, the optimization of transmit power allocation of full-duplex access core network was studied in [22].

The connection between nodes is the fundamental problem of wireless communication networks. It not only reflects link quality, but also determines network performance. In general, for a static and deterministic network composed of n nodes, the nodes can be denoted as a graph , where denotes the vertex set and denotes the edge set. The probabilistic adjacency matrix, denoted by , describes the state of the nodal connection and is a symmetric matrix. If there exists an edge between nodes and , , otherwise 0, that is, nodes and are not connected [23]. From the perspective of wireless communication, the power of transmitter determines the connection between nodes. Therefore, intuitively, the adjacency matrix directly reflects the transmit power of the system. In this way, based on the probabilistic adjacency matrix, the research on transmit power allocation of wireless communication system is plausible. However, as mentioned above, as the topology of the CFSN is highly dynamic and nodes are randomly connected, we proposed a method of transmit power allocation with multi-QoS by establishing the nodal distance distribution model, defining the binary probabilistic adjacency matrix, and adopting the CSMA/CA channel access mechanism. Our transmit power allocation method is proposed on the following premises: (1) low traffic load, (2) finite overall transmit power, and (3) the star topology of the network, that is, one node is responsible for the earth communication, and others are connected to the node by frequency division multiple access (FDMA) with subcarrier binary phase-shift keying (BPSK) modulation [24]. The purpose of optimal transmit power allocation is to minimize the average PER at the access point (AP) based on the probabilistic adjacency matrix.

The structure of this paper is organized as follows: Section 2 introduces and analyzes the basic model of CFSN and describes the definition of the probabilistic adjacency matrix of this model. In Section 3, the optimization-theoretic model is present. And based on the model, the optimized transmit power allocation strategy is proposed. In Section 4, the simulation results of the PER and delay in the network are presented, focusing on the impact of the probabilistic adjacency matrix, the traffic load, and the adopted power allocation strategy. Finally, Section 5 concludes the paper.

2. System Model

2.1. Definition of a Simplified Model in CFSN

The link connection characteristics between nodes in CFSN depend on the relative orbits, transmit power, and receiving sensitivity of the nodes. Therefore, the following definition is given:

Definition 1. In ECI coordinates [24, 25], for each time slot in an orbital hyperperiod of CFSN, a node can be defined as a couple , where is the location of node and is its transmit power.
Given a spacecraft network consisting of L nodes with interference limited, if node can receive the signal transmitted by the node , the bit-signal-to-noise ratio at the receiving node should satisfy the following equation:where is the additive white noise power spectral density of receiver, is Boltzmann’s constant, is the noise temperature, is the transmit power of node , is the Euclidean distance between and , is the rate of data transmission, and are the gains of the transmit and receive antennas, is the working wavelength, and is the receiving sensitivity and depends on the modulation mode, etc. For QPSK modulation, when the bit error rate is less than , . In order to facilitate the calculation, it can be assumed that .
In CFSN, , and are the upper bound and the lower bound of nodal distance, respectively [24].
Therefore, equation (1) can be simplified as follows:According to the following equation, the distance threshold of any successful connection between and can be obtained:That is, if , then and are connected to each other. And we can derive from equation (3):Therefore, the network can be defined as follows:

Definition 2. A CFSN of L nodes, in each time slot of its orbital hyperperiod, can be defined as an order set: , where is the location of the AP.
Based on Definition 2, the definition of binary adjacency matrix in CFSN can be obtained.

Definition 3. For a given CFSN , in each time slot of its orbital hyperperiod, the binary adjacency matrix is given by , whereAnd the complement of corresponds toTo simplify the calculation, the number of ones in the adjacency matrix is denoted by , and the complementary adjacency is given by the number of zeros in adjacency matrix, denoted by .
Moreover, for the sake of analysis, for , two index sets of nodes are defined: and . They represent the node index that node can receive from and transmit to, respectively.

2.2. Nodal Mobility Model and Distance Distributions in CFSN

To accomplish the cluster flight model within bounded distance, twin-satellite mode was adopted to study the nodal mobility model. The node position is uniformly distributed on sphere within radius as shown in Figure 1.

So the mobility model within bounded distance for CFSN can be defined as follows.

Definition 4. In ECI coordinates, if the position sets of nodes in CFSN are at initial time , the position set is , and the positions are uniformly distributed within sphere at time , where and are center and radius of the sphere, respectively. Moreover, positions among all nodes are mutually independent and independent of all previous locations.
Based on the nodal mobility model, the nodal distance distribution can be described in Figure 2. Nodes are assumed to be uniformly located in a circle of the two-dimensional plane.
In Figure 2, the coordinates of the transmitter A and the receiver B are and , respectively, where, and , and D with the value is the distance between the centers of the two circle. The random variables with probability density functions are given bywhere is the uniform distribution over range , and the subscripts denote , and , , respectively.
Therefore, the distance between the transmitter A and the receiver B is given byDespite the simplicity of equation (8), the derivation of the distance density cannot be given in closed form. According to Glivenko–Cantelli Lemma, adopting empirical statistical method and eighth-order polynomial approximation, the probability density function of the distance between nodes with eighth-order polynomial can be denoted byWhen and , the polynomial coefficients of the empirical probability density function are provided in Table 1.

2.3. Determination of Probabilistic Adjacency Matrix in CFSN

According to the foregoing, the definition of the probabilistic adjacency matrix for CFSN is described as follows.

Definition 5. Given the CFSN , in each time slot of an orbital hyperperiod, the probabilistic adjacency matrix is a matrix. is the element, and . is the probability that and successfully connect. The diagonal entries are all equal to 1.
Because the satellite has the capability of storage and forwarding, the diagonal entries are all equal to 1.
If , the connectivity probability between and is given by

3. Packet Error Rate, Transmission Delay, and Power Optimization of the System

In order to analyze the relationship between the system’s PER, transmission delay, and the probabilistic adjacent matrix, the following two counting processes [26] are assumed: in the time interval , the number of times that packet transmission of other nodes has been checked by node i is a homogeneous Poisson process with intensity (dimension: pck/s), and the number of times that a node i has transmitted packet is also a homogeneous Poisson process with intensity . Besides, the processes associated with different nodes are independent of each other, all data can be transmitted successfully within , and the data transmission length is , where is the packet length (dimension: b/pck) and is the transmission data rate (dimension: b/s).

3.1. Packet Error Rate and Delay of the System
3.1.1. Packet Error Rate

For a certain time slot in the orbital hyperperiod of CFSN, if only the adjacency matrix is taken into account, the bidirectional packet error intensity accumulated by node i will bewhere represents the intensity of transmission errors that cannot receive the packet because the nodes are not connected. represents the intensity of transmission errors due to interference from other nodes that cannot receive .

When the system has low traffic load, i.e., , the coefficients and in equation (11) can be expressed, respectively, as follows:

Using the expressions above for coefficients and in (12) and (13), the transmission error intensity of node i is given bywhere denotes the expected value. Therefore, the overall network error intensity is given by

According to the previous assumption, the overall network intensity is, so the probability of packet error, i.e., the ratio between the overall network error intensity and the generation intensity (given by ), is as follows:

3.1.2. Transmission Delay

In order to analyze the relation between transmission delay and probabilistic adjacency matrix, the CSMA/CA mechanism is assumed to be adopted for the node access in CFSN. If the channel idle assessment time is , and the access back-off time of the busy channel is , the transmission delay connecting node i that can receive the data packet is as follows:

Therefore, the average transmission delay of CFSN can be estimated as follows:

3.2. Optimal Transmit Power Allocation

The purpose of the transmit power allocation is to optimize the QoS of the CFSN system. It can be seen from equation (14) that the problem of optimizing the transmit power allocation for minimizing the PER of the CFSN system is equivalent to maximizing , i.e., maximizing the sum of the nondiagonal elements in the probabilistic adjacency matrix.

If , and total power of CFSN system is , the discrete optimization problems can be formulated as follows.

Problem. For a given CFSN , in each time slot of its orbital hyperperiod, each node chooses a transmit power , and the transmit power allocation with the smallest PER needs be optimized as follows:Using the conditional probability method, (10) can be rewritten asSubstituting (20) to (19), because maximizing is equivalent to maximizing the sum of the nondiagonal elements in the probabilistic adjacency matrix, the optimization problem can be converted intoThe constraint C1 and C2 denote that the total network transmit power is finite and the distance between connected nodes is within the threshold of node and . Since , the nonconvexity of constraint C2 can be proved by the second-order partial derivative with respect to the variable.
It is noted that equation (21) is a nonlinear multichoice knapsack problem, which can be solved by Monte Carlo method [27, 28]. Moreover, it is also noted from equations (16) and (18) that the PER and the average delay of the system are linear with .
The outline of the power allocation algorithm is given as follows.Step 1: initialization: set , , , , , , .Step 2: set , and choose a large value for N, where: k = 1, N = total number of trials.Step 3: generate a uniformly distributed random transmit power for each node, (using MATLAB, denotes a matrix where its elements range from 0 to 1. In this paper, the transmit power for each node is , so RN can be computed by mapping ).Step 4: calculate the probability of the AP received from node : compute the transmit distance threshold from node i to AP using equation (3), and compute the connectivity probability of this link using equation (21).Step 5: calculate the probability of the AP transmitted to node : compute the transmit distance threshold from AP to node i using equation (3), and compute the connectivity probability of this link by using equation (21).Step 6: sum the connectivity probability of Step 4 and Step 5.Step 7: add 1 to k, if , calculate the max sum of Step 6 and end; otherwise, go to Step 3.NOTE: Monte Carlo algorithm means that the more the samples, the more approximate the optimal value [29]. As , the objective value tends to the optimal solution.

4. Simulation Analysis

In order to simulate and analyze the influence of node transmit power on the QoS performance of CFSN, packet error rate, average delay, and two allocation strategies are considered: (1) each node has the same transmit power and (2) the optimal transmit power varies from node to node and is allocated using strategy presented in Section 3.2.

The simulations have been carried out referring to star topology, i.e., all nodes transmit (receive) directly to (from) the AP, and using different values of the overall network transmit power and, consequently, different values of the transmit powers allocated to the spacecraft. Simulation parameter settings are as follows: (1) transmit power. It is assumed that the module used to collect solar panels is powered by microwave wireless power transfer to other modules, and each module provides an effective total power of 0.13 W [24]. (2) The number of nodes of the cluster flight spacecraft and , and the orbital elements are derived from [24]. (3) The noise temperature is 300 K. The QPSK modulation is adopted between nodes. The transmission data rate is 100 kbps, the operating frequency is S-band, and f = 2.2 GHz [30]. The gains of the transmit and receive antennas are 1 [31]. Other parameters are listed in Table 2 and one considered topology with is shown in Figure 3.

4.1. Each Node with the Same Transmit Power
4.1.1. Impact of the Traffic Load on PER

According to STK and MATLAB simulation, the nodal distance between AP and other nodes can be listed in Tables 3 and 4. Tables 3 and 4 show that the distance between AP and other nodes changes slowly in different time slots of an orbital hyperperiod, and the distance also has a slight change in corresponding time slots of different orbital hyperperiods.

Since the transmit power of each node is the same, it can be assigned to by equation (4), and the value satisfies all the constraints. The PER is shown as a function of the offered traffic load in CFSN. Comparing the first, the second, and the third time slots in an orbital hyperperiod, the PER hardly changes. The simulation curves are almost overlapping, as shown in Figure 4. Because the change of the distance is not obvious, the probability of connection between nodes has small change in different sequential time slots of any orbital hyperperiod. Comparing the first and the second orbital hyperperiod in a time slot, the phenomenon in Figure 5 is the same as Figure 4. The reason is that the probability of connection between nodes has a slight change in corresponding time slots of different hyperperiods, and the PER is likely to change periodically.

In order to verify the impacts of different transmit power on PER, the same transmit power is allocated to each node with different values in the first time slot. The result is shown in Figure 6. Comparing the curves referring to each node with 0.13 W, each node with 0.08 W with referring to each node with 0.03 W, it can be concluded that if the transmit power is higher, the PER will be lower under the satisfied constraints. As shown in Figures 46, the PER increases with the offered traffic load. In other words, when the traffic load is low, the number of collisions at the AP is likely to be low. Instead, when the traffic load is lager, the probability that two nodes transmit at the same time increases and, subsequently, the PER increases as well.

Keeping other parameters unchanged, the PER under two cases of different numbers of nodes was analyzed. As shown in (b) of Figures 46, the PER also increases with the offered traffic load. Comparing (a) and (b) of Figures 46, it can be further obtained that the number of nodes increases in the network and the PER increases.

4.1.2. Impact of Each Time Slot on PER and Delay

According to the analytical results in Section 3.2, the performance, in terms of PER, depends only on the adjacency matrix. For this scenario, the packet generation rate is set to . We consider the first time slot of the first orbital hyperperiod in CSFN. Figure 7 shows the impact of on PER and delay. When the number of nodes increases, the dimension of the adjacency matrix increases, which in turn affects . The dotted line and solid line indicate 7 nodes and 5 nodes, respectively. The results show that the larger the becomes, the lower the PER will be, and the lager the delay will be. However, the PER and the delay increase with the number of nodes.

4.2. Optimal Transmit Power Allocation Strategy

In this section, we present the impact of the adopted transmit power allocation strategy on the PER of CFSN. In particular, we consider strategies such that optimized and uniformly distributed transmit power is allocated to each node. The optimized transmit power is different at each node and is set according to the power allocation algorithm presented in Section 3.2, where the total transmit power is assigned to each node in order to minimize the PER in the CFSN. This method leads to allocating transmit power reasonably.

Figure 8 shows the impacts of optimized and uniformly distributed transmit power on PER in the first time slot of the first orbital hyperperiod in CSFN. Different values of total network transmit power are considered. Under the transmit power allocation strategy presented in Section 3.2, in Figure 8, a performance comparison between scenarios with and without the use of the proposed transmit power allocation strategy is presented. Comparing the curves referring to with referring to , it can be concluded that if the total transmit power is higher, the sum of the nondiagonal elements in the probabilistic adjacency matrix will be larger and the PER will be lower under the satisfied constraints. Comparing (a) and (b) of Figure 8, it can be concluded that the increment in the number of nodes will increase the PER of the system. For the sake of comparison, in Figure 8, the PER in scenarios where no transmit power allocation strategy is used (solid line) is also shown. In this case, the performance is worse than the case with the optimized transmit power allocation strategy. In fact, given a value of total network transmit power, the proposed transmit power allocation strategy allows to maximize the sum of the nondiagonal elements in the probabilistic adjacency matrix and, therefore, reduce the PER.

5. Conclusion

In this paper, we have presented an optimized transmit power allocation strategy which allows to minimize the PER of CFSN. First, according to the probabilistic adjacency matrix, we have derived a simplified analytical model which describes the performance of CFSN in different time slots, under the assumption of offered traffic load. Then, we have presented optimization-theoretic transmit power allocation algorithm and implemented it under the assumption of finite total network transmit power under two cases of different numbers of nodes. In particular, we have shown the performance depends on the probabilistic adjacency matrix: the sum of the nondiagonal elements and traffic load. Our analytical model has been validated through the Monte Carlo method. This paper has presented the impact of the probabilistic adjacency matrix, the offered traffic load, and transmit power allocation strategy on relevant network performance indicators (PER and delay). Finally, we have verified that the proposed transmit power allocation strategy, by maximizing the sum of the nondiagonal elements in the probabilistic adjacency matrix, allows to minimize the PER for a given total network transmit power at any time slot for CFSN.

Data Availability

The orbital data of CFSN composed of 5 satellites and 7 satellites that used to support the findings of this study are included within the article and from literature [24].

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This research is a project partially supported by the National Natural Science Foundation of China (Grant no. 61362004) and Guizhou Provincial Education Innovation Group Foundation (Grant no. [2017] 031).