Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Service blockage on the downlink in large-scale satellite optical networks: a multi-downlink scheduling for routing selection

Open Access Open Access

Abstract

Over years of space laser communication technology advances, satellite optical networks (SONs) have emerged as a pivotal component in 6 G networks. Satellite services are transmitted from the global view, undergoing transmission through SONs, and being downloaded to the targeted areas. However, the transmission capacity of satellites passing through the areas where users are concentrated may be insufficient to download services transmitted worldwide. This problem exists in various kinds of satellite networks and may cause a large amount of service congestion. In this paper, we propose a multi-downlink delivery routing selection (MDD-RS) strategy to study the total utilization of transmission capacity of SONs. We construct an integer linear programming (ILP) model to establish an optimal case study for minimal network capacity occupation. Also, we design an online option, MDD-RS heuristic algorithm, dynamically calculating path routes, considering bandwidth allocation and resource constraints. A comparative analysis against the conventional single-downlink scheme reveals superior performance of the MDD-RS heuristic algorithm, with a reduction in blocking probability of 0.129 and an improvement in bandwidth utilization of 0.032.

© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

In the past decades, satellites have played a pivotal role in television communications, enabling widespread broadcasts to users in regional areas [13]. By exploiting satellite’s high-speed mobility, nearly all television signals can be forwarded, directly attaching homes [1]. Nowadays, satellite communication is regarded as a promising way to achieve connections for global users. It is a complementary composition of 6 G networks to incorporate the spatial and the temporal networks. It helps to address the digital gap between urban and rural areas through remote communications. To ensure a seamless coverage area, satellites are developed into a large-scale constellation, numbering in tens of thousands, so-called mega-constellations [46]. Some well-known constellations include Starlink [7,8], which consists of over 6,000 small satellites as of March 2024, and OneWeb [9,10], having a constellation exceeding 648 satellites. In the near future, it is expected that these satellite constellations will be applied in various communication scenarios facilitated by some advanced technologies, such as the direct connection between satellites and mobile phones.

Traditionally, satellites have predominantly served as signal relays in satellite-to-ground (SG) communication processes. Through the satellite-to-ground links (SGLs) [11], signals can be forwarded by satellites in a bent-pipe manner, merely forwarding signals among users without modifying the data. This process may introduce latency on the waiting of satellite moving right over the destination. To improve transmission performances, researchers consider the way of forwarding signals from one satellite to its nearby satellites until finding the satellite right over the destination. Such an inter-satellite way is based on the establishment of inter-satellite links (ISLs). While the microwave communications in space have some disadvantages like relatively low data rates, and thus limit the application efficiency. Over years of advances in optical communications, researcher trends to exploit laser communication to establish ISLs and have an enhanced transmission capacity [12,13]. Laser communication is a viable way, which can potentially support gigabits-per-second data rates to match the projected performances of satellite internetwork. The benefits of using laser communication were demonstrated with numerous experiments in various projects [1417], like OneWeb and Starlink. Consequently, ISL achieved through laser communication in satellite networks, often referred to as satellite optical networks (SONs), are emerging as promising and forward-looking trends in the evolution of communication networks.

Nowadays, satellite numbers are rapidly increasing to extend the coverage areas and enable more accessibility regardless of time and position. While there are more numbers of ground stations on the ground. Satellites are enabled to connect with several ground stations in visual range at the same time to perform communications. They establish uplinks and downlinks to upload and download services by occupying a certain bandwidth. The connections can generally maintain a short interval until the satellites fly out of the range. During these connection periods, the Satellite-Ground Link (SGL) connections are generally assumed to have similar transmission capacities. This assumption stems from the situation that satellites in the same orbit theoretically maintain consistent speeds and inclinations relative to ground stations, resulting in nearly equal connection durations. However, user distribution across regions varies due to factors like population density and political strategy. Certain regions become primary recipients of a large number of services. During the download process, the satellite’s downlink transmission capacity may be constrained due to the disparity in transmission capacity between ISLs and SGLs. Moreover, a satellite possesses fixed port numbers and link bandwidth for establishing downlinks with ground stations. It will potentially lead to service blockage. Existing solutions, such as the proposal of double-layer constellations to enhance transmission capacity in specific regions [1820], are constrained by technological maturity and substantial costs. Therefore, the problem of how to use the fixed transmission capacity of downlinks to download more services is inevitable and should have a focus.

In the background of large-scale SONs, this paper focuses on the problem of service blockage caused by limited downlink transmission capacity. We provide case studies to see the network performances of increasing satellite numbers to jointly download services. First, we construct the integer linear programming (ILP) model, to seek the optimal results of minimum network capacity utilization as an optimal case study. Second, we propose an online option, i.e. multi-downlink delivery routing selection (MDD-RS) heuristic algorithm, to dynamically calculate the path routes consisting of ISLs and SGLs. MDD-RS algorithm requires the ISLs and SGLs in multiple paths have the common connection periods. Also, the scheduling is subjected to the following two limitations: i) available transmission capacity of ISLs and SGL, and ii) available port numbers of satellites (to the ground stations) and the port numbers of the ground station (to the satellites) to establish SGLs. With the comparison of one downlink scheme as the benchmark, we analyzed the simulation performances of MDD-RS algorithm. Simulation results show better performances of MDD-RS in the items of blocking probability of service delivery and bandwidth utilization of ISLs and downlinks.

The reminder of this paper is organized into 6 sections. Sec. 2 provides a comprehensive overview of SONs and nowadays popular routing schemes in SONs. Sec. 3 addresses the DD-RS problem with the constructive network model and service model. Sec.4 addresses the DD-RS problem and designs the MDD-RS heuristic algorithm. The numerical results are discussed and presented in Sec. 5. Finally, Sec. 6 draws the paper’s conclusion.

2. Background

2.1 SONs

SONs typically compose the space segment and the ground segment. In the space segment, satellites are divided by the location orbits at different altitudes. Low earth orbit (LEO) is in the range of 400-2000km, medium earth orbit is in the range of 2000-35786 km, and geostationary earth orbit is beyond the altitude of 35786 km [2126]. Satellites perform laser communication to transmit high-speed data between each other in SONs, enabling the exchange of large volumes of data. In the ground segment, it mainly includes the ground stations [27] and ground terminals [28]. Ground stations are mostly considered on their transmission performances of uplinks and downlinks [29], and the connection relationship with the moving satellites is an important topic in existing papers. Within the coverage area, satellite and the ground station can establish bidirectional communications with two independent links affected by visible distance, inclination angle, atmospheric conditions, users’ demand pattern and other factors [30]. They can be achieved by multibeam technology, which are always studied on the issues of power and bandwidth allocation [31,32]. Without losing practicality, this context mainly considers the background of the ground stations with unbalanced service distribution and delivers them to SONs. In SONs, service delivery in SONs can be divided into three processes, as shown in Fig. 1(i) The upload of services via uplinks. ii) The forwarding of services between satellites. iii) The download of services via downlinks. In a service’s delivery path, the satellite responsible for the service upload is called access satellite, and the satellite responsible for the service download is called feeder satellite.

 figure: Fig. 1.

Fig. 1. The presentation of SONs.

Download Full Size | PDF

2.2 Routing schemes in SONs

SONs are attracted significant attention in recent years. In Table 1, it is an overview of relevant literature addressing various service congestion challenges across different types of satellite networks. Additionally, the table summarizes various multipath schemes developed specifically for SONs. Despite lots of existing state-of-art routing schemes were designed in SONs, the problem of how to exploit multiple downlinks for a service delivery was less considered. Motivated and differing from the above strategies, we focused on the service blockage caused by the problem of unbalanced capacity between uplink and downlink in this context.

Tables Icon

Table 1. Summary of Routing Schemes in Satellite Optical Networks

3. Network model and problem statement

3.1 SON model

The modeling of SON can be denoted as a graph at arrival time of service request. In the graph G, we formally take satellites and ground stations in the node set V, respectively marked as S and G. E is link set, including ISLs (${L_S}$) and SGLs (${L_G}$). ${L_G}$ includes uplinks (${L_A}$) and downlinks (${L_F}$). The bandwidth of ${L_A}$, ${L_S}$ and ${L_F}$ are marked as ${B_{as}}$, ${B_{ss}}$ and ${B_{sg}}$, respectively. The occupied bandwidth of ${L_A}$, ${L_S}$ and ${L_F}$ are represented as ${B_{a{s_o}}}$, $\textrm{}{B_{s{s_o}}}$ and $\textrm{}{B_{s{g_o}}}$.

3.2 Service model

Service requests are considered as a tuple $r({{s_r},{g_r},{b_r}} )$, where ${s_r}$ is a content and refers to the source satellite, ${g_r}$ is a ground station (i.e., a content), and ${b_r}$ is the required bandwidth. Suppose source satellite ${s_r}$ and $\textrm{G}({V,E} )$ is given. $r({{s_r},{g_r},{b_r}} )$ can be separated into ${r_s}({{s_r},{d_r},{b_r}} )$ and ${r_g}({{s_r},{d_r},{b_r}} )$ for the path calculation in IS process and download process, as shown in Eq. (1). The separation begins from a common source satellite ${s_r}$ to several satellites ${d_r}$, involving x IS paths and x downlinks for service delivery, as depicted in Eq. (2) and Eq. (3). The total bandwidth occupation across x paths equals ${b_r}$. As for the x IS paths and x downlinks of each service request, suppose they belongs to set $\textrm{R}_s^r$ and set $\textrm{R}_g^r$, respectively.

$$r({{s_r},{g_r},{b_r}} )\to {r_s}({{s_r},{d_r},{b_r}} )\parallel {r_g}({{d_r},{g_r},{b_r}} )$$
$${r_s}({{s_r},{d_r},{b_r}} )\to ({\textrm{r}_s^1({{s_r},d_r^1,b_r^1} )\parallel \cdots \parallel \textrm{r}_s^x({{s_r},d_r^x,b_r^x} )} )$$
$${r_g}({{d_r},{g_r},{b_r}} )\to ({r_g^1({d_r^1,{g_r},b_r^1} )\parallel \cdots \parallel r_g^x({d_r^x,{g_r},b_r^x} )} )$$
$$b_r^1 + \cdots + b_r^x = {b_r},0 < b_r^x \le |{{L_G}} |\le |{{L_S}} |$$
$$\{{d_r^1, \cdots ,d_r^x} \}\in S,\textrm{}0 < x \le {N_{gs}}$$
$$\forall r_s^x \in \textrm{R}_s^r,\forall r_g^x \in \textrm{R}_g^r,r \in R$$
$${s_r} \in S,{g_r} \in G$$

3.3 MDD-RS problem statement

The MDD-RS problem is defined as follows. The above-mentioned models of network topology and service request are taken as an input. Given i) a network graph $\textrm{G}({V,E} )$, ii) the service request model $r({{s_r},{g_r},{b_r}} ),r \in \{{1, \cdots ,|R |} \}$, separated into $r_s^x({{s_r},d_r^x,b_r^x} )$ and $r_g^x({d_r^x,{g_r},b_r^x} )$ . iii) The problem addresses the service blockage caused by a limited downlink transmission capacity, i.e., ${L_{A1}} + {L_{A2}} + {L_{A3}} > {L_F}$, as presented in Fig. 2. Both ILP and MDD-RS heuristic algorithm are constrained by, i) the bandwidth capacity of IS and SG satisfied the service requests at least, and ii) the available port numbers of satellites (to the ground stations) and the port numbers of the ground station (to the satellites).

 figure: Fig. 2.

Fig. 2. Diagram of the MDD-RS problem.

Download Full Size | PDF

4. ILP model

In this section, we develop the MDD-ILP model to allocate available transmission capacity for end-to-end services in a static scenario. The ILP model shows the constraints associated with bandwidth scheduling across multiple downlinks, and enables the optimal results. Table 2 summarizes the model parameters.

Tables Icon

Table 2. Parameters Description

4.1 Supposes illustration

There are two assumptions outlined below. i) The first assumption is that the enabled transmission capacity of every ISL and SGL remains fixed over time. ii) The second assumption is that when the diagonal distances and elevation angles between feeder satellites and ground stations are feasible, they are allowed to establish SGLs in a physical view. As is customary, the connections of every S and G can be identified by coordinates.

4.2 ILP model

The object is to minimize the total capacity occupation ${C_T}$ in the SON, as shown in Eq. (4). It is composed of ${C_S}$ in Eq. (5) and ${C_G}$ in Eq. (6). It is for the constraints of capacity allocation in IS, as shown in Eq. (716). The first constraint, as described in Eq. (7), represents the bandwidth constraint for ISLs. It ensures that the bandwidth occupation of any ISL for a service should not exceed its requirement, which is also less than or equal to the initial bandwidth, as depicted in Eq. (8). Equation (9) ensures the network capacity constraint by ensuring that the total bandwidth demand of all service requests does not exceed the total capacity of the network. Equation (10) states the constraint of occupied numbers for one service in an ISL should not be beyond x. Equation (11) imposes the constraint that the number of services on any ISL should not exceed the total number $\textrm{N}$. Equation (12) constrains the number of ISLs that can be assigned to a service, ensuring it is less than the total number of network links. Equation (13) and Eq. (14) ensures the number of ports that can be occupied by a service cannot exceed the allowable number of ports. The third constraint is the flow conservation in IS, as shown in Eq. (15). It ensures the input and output of every intermediate satellite must be equal. The feeder satellites here are only considered for IS service input and output.

Objective function:

$$Minimum\; {C_T} = {C_S} + {C_G}$$
$$\textrm{}0 \le {C_S} \le {B_{ss}} \cdot |{{L_S}} |,0 \le {C_G} \le {B_{sg}} \cdot |{{L_G}} |$$
$${C_S} = \mathop \sum \nolimits_{r \in R} \mathop \sum \nolimits_{({i,j} )\epsilon {L_S}} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,j} )}^{r_s^x}$$
$${C_G} = \mathop \sum \nolimits_{r \in R} \mathop \sum \nolimits_{({i,m} )\epsilon {L_G}} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,m} )}^{\textrm{r}_g^x}$$

Subjected to,

$$0 \le \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,j} )}^{r_s^x} \le {b_r},\,\forall \textrm{}({i,j} )\in {L_S},\forall r$$
$$0 \le \mathop \sum \nolimits_{r \in R} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,j} )}^{r_s^x} \le {B_{ss}},\textrm{}\forall \textrm{ }({i,j} )\in {L_S}\textrm{}$$
$$\textrm{}0 \le \mathop \sum \nolimits_{r \in R} \mathop \sum \nolimits_{({i,j} )\epsilon Ls} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,j} )}^{\textrm{r}_s^x} \le {B_{ss}} \cdot |{Ls} |\textrm{}$$
$$0 \le \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} K_{({i,j} )}^{r_s^x} \le x,\forall \textrm{}({i,j} )\in {L_S},\forall r$$
$$0 \le \mathop \sum \nolimits_{r \in R} K_{({i,j} )}^{\textrm{r}_s^x} \le |R |,\textrm{}\forall \textrm{}({i,j} )\in {L_S},\forall \textrm{}x$$
$$0 \le \mathop \sum \nolimits_{({i,j} )\epsilon Ls} \mathop \sum \nolimits_{ {\in} ({0,{N_{gs}}} ]} K_{({i,j} )}^{r_s^x} \le |{Ls} |,\; \forall r$$
$$0 \le \mathop \sum \nolimits_{j \in \textrm{S}} K_{({i,j} )}^{r_s^x} \le {N_{ss}},\forall r,\forall \textrm{}i,\forall x$$
$$0 \le \mathop \sum \nolimits_{i \in \textrm{S}} K_{({i,j} )}^{\textrm{r}_s^x} \le {N_{ss}},\forall r,\forall \textrm{}j,\forall x$$
$$\mathop \sum \nolimits_{j \in \textrm{S}} K_{({j,i} )}^{\textrm{r}_s^x} - \mathop \sum \nolimits_{\textrm{j} \in \textrm{S}} K_{({i,j} )}^{\textrm{r}_s^x} = \left\{ {\begin{array}{{cc}} 1&{i = {s_r}}\\ { - 1}&{i = {d_x}}\\ 0&{else} \end{array}} \right.,\; \; \forall r,\forall \textrm{}i,\forall x$$

The following equations are the constraints for the downlink selection. It is for the constraints of capacity allocation in SG, as shown in Eq. (1621). The first constraint is bandwidth constraint of SGLs, as shown in Eq. (16). Equation (17) refers to the constraint of service numbers assigned in any SGL which cannot beyond the total numbers. Equation (18) imposes a constraint on services where the x paths have a total bandwidth of ${b_r}$ if $m = {g_r}$. It is also the flow conservation for a ground station. In this case, Eq. (19) means that the occupied SGL numbers should be larger than x. Equation (20) and Eq. (21) represent the port number constraints for any ground station or satellite, ensuring that they are less than or equal to ${N_{sg}}$.

Subjected to,

$$0 \le \mathop \sum \nolimits_{r \in R} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,m} )}^{\textrm{r}_g^x}{B_{sg}},\,\forall \textrm{}({i,j} )\in {L_G}$$
$$0 \le \mathop \sum \nolimits_{r \in R} K_{({i,m} )}^{\textrm{r}_g^x} \le |R |,\textrm{}\forall \textrm{}({i,j} )\in {L_G},\forall x$$
$$\mathop \sum \nolimits_{i\epsilon S} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} b_r^x \cdot K_{({i,m} )}^{r_g^x} = {b_r},\; \; \; \forall r,m = {g_r}$$
$$\textrm{} < \mathop \sum \nolimits_{i\epsilon S} \mathop \sum \nolimits_{x \in ({0,{N_{gs}}} ]} K_{({i,m} )}^{r_g^x} \le x,\; \; \; \forall r,m = {g_r}$$
$$0 \le \mathop \sum \nolimits_{i \in \textrm{S}} K_{({i,m} )}^{\textrm{r}_g^x} \le {N_{gs}},\; \; \forall r,\forall x,\forall \textrm{}m \in G$$
$$0 \le \mathop \sum \nolimits_{m \in \textrm{G}} K_{({i,m} )}^{r_g^x} \le {N_{sg}},\forall r,\forall x,\forall \textrm{}i \in S$$

5. Online algorithm

ILP models provide optimal results but can be computationally intensive and spend a long running time, especially in large-scale SONs. This section introduces the MDD-RS strategy and its heuristic algorithm. These methods efficiently select feeder satellites and compute multiple routes under the constraints to ground stations.

5.1 MDD-RS strategy

To find multiple paths and downlinks with available resources, we propose an MDD-RS strategy that decomposes one feeder satellite into multiple assistant satellites, collectively delivering services. In MDD-RS strategy, destination ground stations selectively establish SGLs with certain passing satellites. First, we conduct a potential set of satellites for the destination ground station that can establish connections, which can be obtained in the ephemeris table. Then, it is the selection of assistant satellites to calculate paths with available transmission resources. We calculate possible routes by exploiting several times of the Dijkstra algorithm for surplus resources and putting them in a set. Based on the route set, considering the assistant satellites act as destination nodes for service delivery, we go through the set to find whether they have available port numbers and enough bandwidth to fit the demand of requests. Finally, by traversing the above when the total transmission capacity exceeds requirements, determining the number of assistant feeder satellites needed.

5.2 MDD-RS algorithm

In Algorithm 1, as the inputs, $G({V,E} )$, $r({{s_r},{g_r},{b_r}} )$, $r \in R$, $G({V,E} )$, ephemeris table, ${B_{sg}}$, ${B_{ss}}$, ${N_{gs}}$, represents a network having given links in IS and SG with a certain bandwidth and enabled port numbers to be allocated for a series of service requests. As the outputs, we design the MDD-RS algorithm to find x as the number of paths for the content request, ${P_{MDD}}$ as the calculated routes for the request, ${B_{MDD}}$ as the offered bandwidth set for the request.

Algorithm 1 starts by the arrival of requests $r({{s_r},{g_r},{b_r}} )$. It contains two parts, respectively in IS and SG. Considering iteration process for the scheduling in IS is more complex compared to SG, the algorithm first goes to the scheduling in SG and then in IS. First, it is the source and destination definiteness steps. We extract the destination ground station from the request, and construct a possible-connected satellite set (i.e., ${\textrm{D}_{{g_r}}}$) of this ground station (i.e., ${\textrm{D}_{{g_r}}}\{{d_r^1, \cdots ,d_r^x} \}$) by searching the ephemeris table. The algorithm checked the available port numbers of $d_r^x$ in ${\textrm{D}_{{g_r}}}$, which should have at least one. Otherwise, delete the ineligible ones from the ${\textrm{D}_{{g_r}}}$. Then, it goes to the scheduling in IS, i.e., route calculation steps. Based on the given ${s_r}$ and ${\textrm{D}_{{g_r}}}\{{d_r^1, \cdots ,d_r^x} \}$, we use Dijkstra algorithm to calculate the shortest routes between ${s_r}$ and each $d_r^x$, respectively. The x numbers of paths are put into path set of the request ${P_r}$. Next, it is the resource allocation steps. The algorithm sorts the paths according to their number of hops from smallest to largest. Traversing each link of routes in ${P_r}\{{p_r^1, \cdots ,p_r^x} \}$, the algorithm regards the lowest bandwidth of the link in the path as the available bandwidth of the path. The algorithm puts the available bandwidth of each path in set of the request (i.e., ${B_r}\{{b_r^1, \cdots ,b_r^x} \}$), and compare them one by one with ${b_r}$ in the request to judge whether a path having enough bandwidth. If there are $b_r^x$ bigger than ${b_r}$, allocating one of them to the request and ending the scheduling of this request. If there is no $b_r^x$ bigger than ${b_r}$, traversing the ${P_r}\{{p_r^1, \cdots ,p_r^x} \}$ and allocating ${B_r}\{{b_r^1, \cdots ,b_r^x} \}$ for the request until $\sum b_r^x \ge {b_r}$. If bandwidth requirement can be satisfied in advance, delete the surplus routes in the ${P_r}$ and their corresponding bandwidth in ${B_r}$. If the x numbers routes in ${P_r}$ having $\sum b_r^x < {b_r}$, it means the scheduling of the request is failed, otherwise it is successful. Finally, the algorithm output the route set and bandwidth set of successful scheduled requests.

Tables Icon

Table 2. MDD-RS Algorithm

5.3 Complexity analysis

The complexity of MDD-RS algorithm includes three parts, i) complexity of the MDD-RS by traversing satellite set V with O($V$). ii) The complexity of multiple path calculation which is consisted of x times of Dijkstra algorithm which uses two times of “for” loops marked as O(${V^2}$). As for x paths calculated in IS part, the number of iterations needed to transform V satellites for feeder function and converge to a successful route selection solution is in the order of O($V$) times and therefore, since the worst-case complexity of O($x{V^2}$). iii) As for the x downlinks calculation in SG part, the process will go for the bandwidth judgement with O($N$) and N is a constant. It should be noted that will have a relative-high complexity with the increasing numbers of satellites or ground stations.

6. Case study and simulation results

6.1 Simulation settings

In this section, we conduct simulations to evaluate the proposed ILP model and MDD-RS algorithm. As for the basic setting of SON, we consider a Walker constellation (altitude: 1050 km, degree: 89°, planes: 24, phase factor: 9, range: 180°) with 1152 satellites to form a basic topology, whose time-varying connection relationship is shown in Fig. 3(a). We suppose there are 50 ground stations in the terrestrial which have a relative centralized distribution, as shown in Fig. 3(b). The specific connection data between satellites and ground stations was obtained from the simulative constellation data in Satellite Tool Kit. Figure 3(b) is presentation of parts of Ground Stations. The entire holding period was set to 2 hours referred by the practical orbit period, then the whole period was divided into 120 continuous time slices per min by using snapshot. We use part of data as the input of topology to verify the proposed route calculation and bandwidth allocation algorithm. The bandwidth of each ISL was set to 40 Gbps, and the bandwidth of uplinks/downlinks are set to 8 Gbps. We suppose the arrival of services in ground stations followed Poisson distribution. Service numbers (SN) varied in the range of [3000, 4000, 5000, and 6000]. Service request was generated by the program and has the known source satellite, destination ground station, and the required bandwidth. The source and destination satellites were generated by Random distribution within the range of node set, and the generation of required bandwidth follows Normal distribution with 30 standard deviations. Therefore, the services were supposed to transmit in a snapshot topology conducted by the time slices sampled every 1 minute in the algorithm simulation. The designed ILP models are solved by the CPLEX software in version 12.7.1 installed on this computer. The simulations were run on a computer with a 1.8 GHz Intel Core i5-8250U CPU and 8GB RAM. It is noted that ILP may not be suitable for the solution applied in practice due to the complexity of its calculation in time dimension. The simulation results of the ILP model will be referred to refer as an optimal operation diagram, which provides the meaning of model construction and constraint forming. In the dynamic scenario, we developed a C++ event-driven network simulator with a joint constellation with STK To evaluate the performances of the MDD-RS algorithm, with the dynamic route-calculated algorithm always used in existing methods as the benchmark, i.e., Dijkstra-routing selection (D-RS), to calculate the shortest path for services. The numerical performance evaluation for the proposed algorithms with the benchmark is mainly evaluated in terms of blocking probability and bandwidth utilization of downlinks and ISLs.

 figure: Fig. 3.

Fig. 3. (a) 3D presentation of 1152 satellites. (b) The presentation of parts location of 50 ground stations on the Earth.

Download Full Size | PDF

6.2 Simulation results for MDD-ILP and MDD-RS

  • 1. Performance evaluation of the ILP

Figure 4(a) shows the minimal utilization of network capacity when SN is 3000, 4000, 5000, and 6000, respectively. We counted the total utilization of ISLs and SGLs and their partial utilization, as shown in the legend. When SN increases, the network capacity utilization gradually becomes higher. The capacity utilization of ISLs and SGLs show an upward trend, but the increasing of ISLs is slower than that of SGLs. This is because the transmission capacity of ISLs is large and the increase of ISL is not obvious to it. In Fig. 4(b), it is the average number of feeder satellites used for the delivery of services. This performance also indicates the number of IS paths for the services. The results show that the number is irregular as the increase of SNs, and also has a fewer gap between each two of them. This is because that multiple paths can alleviate the bandwidth occupation conflict when the bandwidth of the downlink is tight. Figure 4(c) shows the counted hops for the averaged deliveries of services. The trends are similar to the multiple paths, which are also irregular as SN increases. Also, it shows a fewer gap between each two of them. While it should be noted that the calculation of paths for services will be detoured and thus have lots of extra hops especially when the resources become tight. The results of ILP model can be seen as a case study to seek the optimal results of minimum resource occupation under various SNs at a time point.

  • 2. Performance evaluation of the MDD-RS algorithm

 figure: Fig. 4.

Fig. 4. Performances for the MDD-ILP, i.e., (a) minimal bandwidth utilization of the network respectively when SN = 3000, 4000, 5000 and 6000. (b) The average $x$ under SN = 3000, 4000, 5000 and 6000. (c) The average hops of the minimal results.

Download Full Size | PDF

Figure 5(a) shows the blocking probability of the MDD-RS and D-RS. This performance is calculated by the average blocking probability of SNs in 120 time slices. The figure illustrates a reduction of approximately 0.129 in the blocking probability achieved by MDD-RS compared to D-RS when SN = 3000. With the increase of SNs, the trend of MDD-RS gradually becomes stable, and the D-RS has an increase trend. It means that gap between the two trends changes from a small value to a big value and then small back until the MDD-RS cannot effectively alleviate the blockage. If it occurs to this situation, the current SN is the limit value that can be carried in the SON in the situation of fixed transmission capacity of the downlinks. Therefore, we can draw two conclusions from this performance. One is that using one feeder satellite for service delivery will indeed have the fixed and limited transmission capacity, and thus resulting in blockage. The second is that MDD-RS can support more services delivered in the SONs. Figure 5(b) shows the impact of different time slices on the blocking probability. There are 10 continuous time slices selected from the entire orbit period, and each of them is with different holding times. We can see from the figure that the blocking probability of time slice became bigger as SN increases. Each fixed SN’s blocking probability has a slight and irregular fluctuation in different time slices. This fluctuation might be caused by the connection changes between satellites and ground stations on each time slice. This means that time slices and their holding time have less impact on the SNs that can be carried in SONs. Furthermore, we can also find that the improvement of MDD-RS over SN in SONs has less affection for the time slice changing. Figure 5(c) shows the real-time SN of the MDD-RS and D-RS dealt by each feeder satellite changing in two time slices. The X-axis shows the downlink ID, and the Y-axis shows four sizes of SNs, respectively in two time slices. The Y-axis is composed of two parts, i.e., the front block belonged to the D-RS and the back block belonged to the MDD-RS. As for the SN changing in a time slice, the density of the back part is significantly greater than that of the front part, indicating that more SNs can be carried in the downlink, which is consistent with the conclusion in Fig. 4(b). As for the same value of SN under different time slices, it shows the SN of each downlink fluctuates slightly, which also means that time slices have little effect on this. As for the SN in two time slices, it shows that the downlink ID from 50 to 150 has relatively continuous and stable SNs than the others. It means that ground stations can always establish downlinks with some fixed and specific satellites. That is to say, no matter which constellation it is in, a specific area on the earth can be always covered by the satellites.

 figure: Fig. 5.

Fig. 5. Blocking probability of services under increased SN (i.e., 3000, 4000, 5000, and 6000). A gradual decomposition of blocking probability, i.e., (a) a total result under the comparison of D-RS and MDD-RS with increased SN, (b) results with the averaged blocking probabilities under different service numbers respectively in 1-10 sampled continuous time slices, and (c) results with a real time count of carried service numbers for every feeder satellite in MDD-RS and D-RS algorithms during 2 sampled time slices.

Download Full Size | PDF

In Fig. 6(a), the bandwidth utilization of downlinks is presented for both the MDD-RS and D-RS algorithms. As the SN increases, the MDD-RS curve exhibits a gradual rise, contrasting with the stable curve of the D-RS. The gap between them fluctuates, initially small, then widening before contracting again. Notably, the average downlink utilization ratio can be improved by at least 0.032. Figure 6(b) examines the impact of time slices on downlink transmission capacity utilization. Both curves display a similar trend, indicating that the enhanced bandwidth utilization brought by the MDD-RS is less influenced by time slices. However, the maximum downlink bandwidth utilization ratio approximately reaches 0.4576, while the corresponding blocking probability is about 0.075. This illustrates that there is a room for further improvement in bandwidth utilization. Blockages occur even when there is sufficient transmission capacity of downlinks, hinting at potential issues related to insufficient capacity of ISLs. Figure 6(c) delves into the real-time bandwidth utilization of downlinks, featuring various SNs on two time slices, comparing the MDD-RS and D-RS. The front block corresponds to the D-RS, and the rear block to the MDD-RS. As in Fig. 6(c), the proportional relationship between feeder satellites and downlinks persists, showing trends similar to those observed previously. The real-time data reveals that maximum bandwidth utilization ratio approaches 0.9, a significant contrast to the average time slice in Fig. 6(b), the ratio was less than 0.5. This discrepancy underscores the potential generation of fragmented resources through dynamic service allocation, necessitating further exploration in downlink scheduling.

 figure: Fig. 6.

Fig. 6. Bandwidth utilization of downlinks under an increased SN (i.e., 3000, 4000, 5000, and 6000). Also, a gradual decomposition of it, i.e., (a) a total comparison of MDD-RS and D-RS with increasing service numbers, (b) the averaged bandwidth utilization of downlinks to be compared with different service numbers in 10 continuous time slices, and (c) results with a real-time count of downlink bandwidth under different service numbers respectively in 2 time slices.

Download Full Size | PDF

Figure 7(a) illustrates the bandwidth utilization of ISLs, serving as a reflection of the bandwidth occupation by services. Notably, the bandwidth utilization of ISLs in the MDD-RS scenario exhibits a clear increase with the growing number of services, contrasting with the gradual rise observed in the D-RS. The gap between them widens, reaching a point where the MDD-RS curve stabilizes. Given that ISL bandwidth is generally times larger than that of the downlink, a suitable allocation of bandwidth resources for ISLs becomes imperative. In Fig. 7(b), the influence of time slices on ISL bandwidth utilization is examined. Similarly, the MDD-RS demonstrates a less pronounced sensitivity to variations in time slices, contributing to improved bandwidth utilization of ISLs. Figure 7(c) provides insights into the real-time bandwidth utilization of individual ISLs, with the D-RS occupying the front block and the MDD-RS the rear. Different time slices appear to have a relatively minor impact on real-time bandwidth occupancy. However, both MDD-RS and D-RS show similar uneven utilization of ISLs. The figure depicts concentrated utilization in specific sections, resulting in significantly higher utilization rates, while some ISLs experience lower utilization. This concentration is attributed to the MDD-RS’s ability to support more services, contributing to a relatively focalized utilization pattern. The additional services supported by MDD-RS exacerbate uneven resource utilization among ISLs, potentially leading to partial service blockages, which is an aspect to be explored further in the context of flexible resources within SONs.

 figure: Fig. 7.

Fig. 7. Bandwidth utilization of ISLs under an increased SN (i.e., 3000, 4000, 5000, and 6000), i.e., (a) a total comparison of MDD-RS and D-RS, (b) averaged utilizations changing with different service numbers during 10 continuous time slices, and (c) real-time utilizations of ISLs under different service numbers.

Download Full Size | PDF

7. Conclusions

In this paper, we propose MDD-RS strategy by using possible satellites as the feeder to jointly download services. The main work of this paper is concluded as follows: i) Construct ILP models to select available satellites and schedule multiple routes for the ground stations under the object of a minimum network capacity and subjected to the constraints of bandwidth and port number. ii) The MDD-RS heuristic algorithm is designed to dynamically compute routes with multiple feeder satellites for service delivery. Comparative analysis illustrates the effectiveness of MDD-RS algorithm in enhancing network capacity while maintaining a certain performance guarantee. The results indicate that a higher number of services can be efficiently delivered within the networks, accompanied by improvements in the utilization of satellites as the feeders. In the future, we will explore methods to efficiently utilize dynamic network resources and optimize their utilization. This could leverage machine learning techniques or enhancing the adaptability of routing algorithms.

Funding

National Natural Science Foundation of China (62350001); Funds for Creative Research Groups of China (62021005).

Acknowledgements

This work has been supported in part by the Funds for Creative Research Groups of China (62021005).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. M. Marchese, A. Moheddine, and F. Patrone, “IoT and UAV integration in 5 G hybrid terrestrial-satellite networks,” Sensors 19(17), 3704 (2019). [CrossRef]  

2. R. De Gaudenzi, O. Del Rio Herrero, G. Gallinaro, et al., “Random access schemes for satellite networks, from VSAT to M2M: A survey,” Int. J. Satell. Commun. Netw. 36(1), 66–107 (2018). [CrossRef]  

3. P. Wang, J. Zhang, X. Zhang, et al., “Convergence of satellite and terrestrial networks: A comprehensive survey,” IEEE Access 8, 5550–5588 (2020). [CrossRef]  

4. K. Ahan and I. F. Akyildiz, “Designing large-scale constellations for the internet of space things with cubesats,” IEEE Internet Things J. 8(3), 1749–1768 (2021). [CrossRef]  

5. K. Ahan and I. F. Akyildiz, “Large-scale constellation design for the Internet of Space Things/CubeSats,” 2019 IEEE Globecom (IEE2019), pp. 1–6.

6. S. Liu, Z. Gao, Y. Wu, et al., “LEO satellite constellations for 5 G and beyond: How will they reshape vertical domains?” IEEE Commun. Mag. 59(7), 30–36 (2021). [CrossRef]  

7. J. C. McDowell, “The low earth orbit satellite population and impacts of the SpaceX Starlink constellation,” Astrophys. J. Lett. 892(2), L36 (2020). [CrossRef]  

8. K. Joe, M. Neinavaie, and Z. M. Kassas, “The first carrier phase tracking and positioning results with Starlink LEO satellite signals,” IEEE Trans. Aerosp. Electron. Syst. 58(2), 1487–1491 (2022). [CrossRef]  

9. The OneWeb satellite system, Handbook of Small Satellites: Technology, Design, Manufacture, Applications, Economics and Regulation, Springer, Cham, 1-10 (2020).

10. O. B. Ogutu and E. J. Oughton, “A techno-economic framework for satellite networks applied to low earth orbit constellations: Assessing Starlink, OneWeb and Kuiper,” IEEE Access 9, 141611–141625 (2021). [CrossRef]  

11. G. Lipika, A. Kand, M. Pradhan, et al., “Automated ground station design for an amateur LEO satellite system,” 2019 IEEE Aerospace Conference (IEEE, 2019), pp.1–3.

12. W. Chien, C. Lai, M. S. Hossain, et al., “Heterogeneous space and terrestrial integrated networks for IoT: Architecture and challenges,” IEEE Netw. 33(1), 15–21 (2019). [CrossRef]  

13. G. Marco and M. Zorzi, “Satellite communication at millimeter waves: A key enabler of the 6 G era,” International Conference on Computing, Networking and Communications (IEEE, 2020), pp. 1–6.

14. N. Pachler, I. del Portillo, E. F. Crawley, et al., “An updated comparison of four low earth orbit satellite constellation systems to provide global broadband,” 2021 IEEE International Conference on Communications (IEEE, 2021), pp.1–7.

15. B. Rödiger, C. Menninger, C. Fuchs, et al., “High data-rate optical communication payload for CubeSats,” 2020 Laser Communication and Propagation through the Atmosphere and Oceans IX., (Online Only).

16. A. A. Kazemi, “Intersatellite laser communication systems for harsh environment of space,” 2013 Photonic Applications for Aerospace, Commercial, and Harsh Environments IV.

17. A. Talamante, J. D. Bowman, D. C. Jacobs, et al., “Deployable optical receiver array cubesat,” 2021 Small Satellite Conference, pp. 1–10.

18. P. Wang, B. Di, and L. Song, “Multi-layer LEO satellite constellation design for seamless global coverage,” Global Communications Conference (IEEE, 2021), pp. 01–06.

19. M. Hu, J. Li, C. Cai, et al., “Software defined multicast for large-scale multi-layer LEO satellite networks,” IEEE Trans. Netw. Serv. Manage. 19(3), 2119–2130 (2022). [CrossRef]  

20. H. Dang, Y. Zhao, Y. Jing Y, et al., “Link planning schemes for uninterrupted inter-layer communication in dual-layer LEO optical satellite networks,” International Conference on Optical Communications and Networks (IEEE, 2022), pp. 1–3.

21. T. Hiriart, J. F. Castet, J. M. Lafleur, et al., “Comparative reliability of GEO, LEO, and MEO satellites,” in 2009 Proceedings of the International Astronautical Congress, pp.D.1.

22. N. Hiroki, D. Kudoh, N. Kato, et al., “Load balancing and QoS provisioning based on congestion prediction for GEO/LEO hybrid satellite networks,” Proc. IEEE 99(11), 1998–2007 (2011). [CrossRef]  

23. M. Emilio, “Geocentric spherical surfaces emulating the geostationary orbit at any latitude with zenith links,” Future Internet 12(1), 16 (2020). [CrossRef]  

24. Y. Su, Y. Liu, Y. Zhou, et al., “Broadband LEO satellite communications: Architectures and key technologies,” IEEE Wirel. Commun. 26(2), 55–61 (2019). [CrossRef]  

25. Satellite Orbits for Communications Satellites, Handbook of Satellite Applications. Springer, Cham, pp. 99-120 (2017).

26. P. Gu, R. Li, C. Hua, et al., “Dynamic cooperative spectrum sharing in a multi-beam LEO-GEO co-existing satellite system,” IEEE Trans. Wireless Commun. 21(2), 1170–1182 (2022). [CrossRef]  

27. I. Del Portillo, B. G. Cameron, and E. F. Crawley, “A technical comparison of three low earth orbit satellite constellation systems to provide global broadband,” Acta Astronaut. 159, 123–135 (2019). [CrossRef]  

28. A. Paris, I. Del Portillo, B. Cameron, et al., “A genetic algorithm for joint power and bandwidth allocation in multibeam satellite systems,” 2019 IEEE Aerospace Conference (IEEE, 2019), pp. 1–15.

29. G. Rajeev and N. BenAmmar, “Framework for unifying 5 G and next generation satellite communications,” IEEE Netw. 32(5), 16–24 (2018). [CrossRef]  

30. S. Fu, J. Gao, and L. Zhao, “Integrated resource management for terrestrial-satellite systems,” IEEE Trans. Veh. Technol. 69(3), 3256–3266 (2020). [CrossRef]  

31. I. del Portillo, B. Cameron, and E. Crawley, “Ground segment architectures for large LEO constellations with feeder links in EHF-bands,” 2018 IEEE Aerospace Conference (IEEE, 2018), pp. 1–14.

32. S. Nasir, H. Almorad, H. Dahrouj, et al., “Point-to-point communication in integrated satellite-aerial 6 G networks: State-of-the-art and future challenges,” IEEE Open J. Commun. Soc. 2, 1505–1525 (2021). [CrossRef]  

33. L. Zhang, F. Yan, Y. Zhang, et al., “A routing algorithm based on link state information for LEO satellite networks,” Globecom (IEEE, 2020), pp. 1–6.

34. Y. Dong, S. Zhao, H. Ran, et al., “Routing and wavelength assignment in a satellite optical network based on ant colony optimization with the small window strategy,” J. Opt. Commun. Netw. 7(10), 995–1000 (2015). [CrossRef]  

35. X. Sun and S. Cao, “A routing and wavelength assignment algorithm based on two types of LEO constellations in optical satellite networks,” J. Lightwave Technol. 38(8), 2106–2113 (2020). [CrossRef]  

36. P. Kumar, S. Bhushan, D. Halder, et al., “fybrrLink: Efficient QoS-aware routing in SDN enabled future satellite networks,” IEEE Trans. Netw. Serv. Manage. 19(3), 2107–2118 (2022). [CrossRef]  

37. X. Sun and S. Cao, “Wavelength routing assignment of different topological optical networks based on typical LEO satellite constellations,” International Conference on Wireless and Satellite Systems, pp 612–629.

38. Y. Ning, Y. Zhao, X. Li, et al., “Load-balancing routing algorithm against inter-satellite link congestion in LEO satellite optical networks,” Optical Fiber Communication Conference (IEEE, 2022), pp. 1–3.

39. Y. Jing, Y. Zhao, H. Wang, et al., “Deep-learning-based path computation without routing convergence in optical satellite networks,” J. Opt. Commun. Netw. 15(5), 294–303 (2023). [CrossRef]  

40. Y. Zhang, W. Wang, Y. Zhao, et al., “Access point selection based on regular coding in Walkerdelta satellite optical networks,” Optical Fiber Communication Conference (IEEE, 2023), pp. 1–3.

41. D. Li, X. Li, Y. Liu, et al., “Time window intersection (TWI) based snapshot configuration scheme in terrestrial-satellite optical networks,” Opto-Electronics and Communications Conference (IEEE, 2023), pp 1–5.

42. F. Jiang, Q. Zhang, Z. Yang, et al., “A space–time graph based multipath routing in disruption-tolerant earth-observing satellite networks,” IEEE Trans. Aerosp. Electron. Syst. 55(5), 2592–2603 (2019). [CrossRef]  

43. C. Wang, Z. Ren, W. Cheng, et al., “CDMR: Effective computing-dependent multi-path routing strategies in satellite and terrestrial integrated networks,” IEEE Trans. Netw. Sci. Eng. 9(5), 3715–3730 (2022). [CrossRef]  

44. P. Du, X. Li, Y. Lu, et al., “Multipath TCP over LEO satellite networks,” 2015 International Wireless Communications and Mobile Computing Conference, pp. 1–6.

45. Q. Zhao, P. Du, M. Gerla, et al., “Software defined multi-path tcp solution for mobile wireless tactical networks,” 2018 IEEE Military Communications Conference (IEEE, 2018), pp. 1–9.

46. F. Tang, H. Zhang, and L. T. Yang, “Multipath cooperative routing with efficient acknowledgement for LEO satellite networks,” IEEE Trans. on Mobile Comput. 18(1), 179–192 (2019). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. The presentation of SONs.
Fig. 2.
Fig. 2. Diagram of the MDD-RS problem.
Fig. 3.
Fig. 3. (a) 3D presentation of 1152 satellites. (b) The presentation of parts location of 50 ground stations on the Earth.
Fig. 4.
Fig. 4. Performances for the MDD-ILP, i.e., (a) minimal bandwidth utilization of the network respectively when SN = 3000, 4000, 5000 and 6000. (b) The average $x$ under SN = 3000, 4000, 5000 and 6000. (c) The average hops of the minimal results.
Fig. 5.
Fig. 5. Blocking probability of services under increased SN (i.e., 3000, 4000, 5000, and 6000). A gradual decomposition of blocking probability, i.e., (a) a total result under the comparison of D-RS and MDD-RS with increased SN, (b) results with the averaged blocking probabilities under different service numbers respectively in 1-10 sampled continuous time slices, and (c) results with a real time count of carried service numbers for every feeder satellite in MDD-RS and D-RS algorithms during 2 sampled time slices.
Fig. 6.
Fig. 6. Bandwidth utilization of downlinks under an increased SN (i.e., 3000, 4000, 5000, and 6000). Also, a gradual decomposition of it, i.e., (a) a total comparison of MDD-RS and D-RS with increasing service numbers, (b) the averaged bandwidth utilization of downlinks to be compared with different service numbers in 10 continuous time slices, and (c) results with a real-time count of downlink bandwidth under different service numbers respectively in 2 time slices.
Fig. 7.
Fig. 7. Bandwidth utilization of ISLs under an increased SN (i.e., 3000, 4000, 5000, and 6000), i.e., (a) a total comparison of MDD-RS and D-RS, (b) averaged utilizations changing with different service numbers during 10 continuous time slices, and (c) real-time utilizations of ISLs under different service numbers.

Tables (3)

Tables Icon

Table 1. Summary of Routing Schemes in Satellite Optical Networks

Tables Icon

Table 2. Parameters Description

Tables Icon

Table 3 Table 2. MDD-RS Algorithm

Equations (26)

Equations on this page are rendered with MathJax. Learn more.

r ( s r , g r , b r ) r s ( s r , d r , b r ) r g ( d r , g r , b r )
r s ( s r , d r , b r ) ( r s 1 ( s r , d r 1 , b r 1 ) r s x ( s r , d r x , b r x ) )
r g ( d r , g r , b r ) ( r g 1 ( d r 1 , g r , b r 1 ) r g x ( d r x , g r , b r x ) )
b r 1 + + b r x = b r , 0 < b r x | L G | | L S |
{ d r 1 , , d r x } S , 0 < x N g s
r s x R s r , r g x R g r , r R
s r S , g r G
M i n i m u m C T = C S + C G
0 C S B s s | L S | , 0 C G B s g | L G |
C S = r R ( i , j ) ϵ L S x ( 0 , N g s ] b r x K ( i , j ) r s x
C G = r R ( i , m ) ϵ L G x ( 0 , N g s ] b r x K ( i , m ) r g x
0 x ( 0 , N g s ] b r x K ( i , j ) r s x b r , ( i , j ) L S , r
0 r R x ( 0 , N g s ] b r x K ( i , j ) r s x B s s ,   ( i , j ) L S
0 r R ( i , j ) ϵ L s x ( 0 , N g s ] b r x K ( i , j ) r s x B s s | L s |
0 x ( 0 , N g s ] K ( i , j ) r s x x , ( i , j ) L S , r
0 r R K ( i , j ) r s x | R | , ( i , j ) L S , x
0 ( i , j ) ϵ L s ( 0 , N g s ] K ( i , j ) r s x | L s | , r
0 j S K ( i , j ) r s x N s s , r , i , x
0 i S K ( i , j ) r s x N s s , r , j , x
j S K ( j , i ) r s x j S K ( i , j ) r s x = { 1 i = s r 1 i = d x 0 e l s e , r , i , x
0 r R x ( 0 , N g s ] b r x K ( i , m ) r g x B s g , ( i , j ) L G
0 r R K ( i , m ) r g x | R | , ( i , j ) L G , x
i ϵ S x ( 0 , N g s ] b r x K ( i , m ) r g x = b r , r , m = g r
< i ϵ S x ( 0 , N g s ] K ( i , m ) r g x x , r , m = g r
0 i S K ( i , m ) r g x N g s , r , x , m G
0 m G K ( i , m ) r g x N s g , r , x , i S
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.