This application claims priority to Chinese Patent Application No. 201811353219.4, filed on Nov. 14, 2018, which is hereby incorporated by reference in its entirety.
The present disclosure relates to the field of a network selection, and in particular to a network selection method and an apparatus for integrated cellular and drone-cell networks.
In order to provide better network dada service and alleviate the congestion of cellular network, a solution of resorting to drone-cells (i.e., low-altitude drones equipped with transceivers) to offload traffic from the congested cellular network can be adopted. To make full use of resources of drone-cells, a significant challenge is how to make an efficient and fair network selection for integrated cellular and drone-cell networks.
An approach for solving a network selection problem is a game-theory-based approach. In the game-theory-based approach, the network selection problem is firstly modeled as a game and then a centralized/distributed approach is exploited to achieve an equilibrium. For example, Cheung et al. formulated the network selection problem as a Bayesian game under the condition that mobility information of users is not complete. Then, they proposed a distributed approach with good convergence properties to achieve the Bayesian Nash equilibrium.
The existing game-theory-based approaches consider the interaction and competition among users. Most of the approaches, however, studied the network selection problem under the situation with quasi-static or predictable network states. However, the integrated cellular and drone-cell networks are highly dynamic, and the network state is hard to predict. Therefore, it is difficult for the existing game-theory-based approaches to solve the network selection problem of the integrated cellular and drone-cell networks.
The present disclosure provides a network selection method and an apparatus for integrated cellular and drone-cell networks, aiming at solving the problem that the existing game-theory-based approaches cannot solve the network selection problem for the integrated cellular and drone-cell networks because the integrated cellular and drone-cell networks are highly dynamic, and the network state is hard to predict.
In a first aspect, the present disclosure provides a network selection method for integrated cellular and drone-cell networks, including:
acquiring a dynamic network model and a dynamic user model, where the dynamic network model includes at least a location model of the drone-cell networks, a capacity model of the cellular network and a capacity model of the drone-cell networks, and the dynamic user model includes at least a location model and a transmission rate model of users;
generating accessible network sets of the users according to the location model of the drone-cell networks and the location model of the users;
generating a random event vector according to the capacity model of the cellular network, the capacity model of the drone-cell networks, the accessible network sets of the users and the transmission rate model, where the accessible network sets of the users include the drone-cell networks and/or the cellular network:
generating an action vector according to the random event vector where the action vector is used for indicating that the users choose to access to the drone-cell networks and/or the cellular network:
obtaining an individual utility of each user according to the action vector and the random event vector;
constructing a first selection model, where the first selection model includes a first objective function and a first constraint, where the first objective function is a proportional fairness function with time average of individual utilities of the users as independent variables, and the first constraint includes at least a first coarse correlated equilibrium constraint, a first minimum individual time average utility constraint and a first action probability constraint, where the first coarse correlated equilibrium constraint is used for constraining the time average of the individual utilities and first auxiliary variables, the first minimum individual time average utility constraint is used for constraining the time average of the individual utilities, and the first action probability constraint is used for constraining the action probability under the condition of the random event vector;
obtaining the time average of the individual utilities according to the individual utilities, the random event probability and the action probability under the condition of the random event vector, where the action probability under the condition of the random event vector is the probability that the users execute the action vector under the condition that the random event vector occurs, and the random event probability is the probability that the random event occurs; and
obtaining a value of the action probability according the first selection model, and determining networks that the users choose to access according to the value of the action probability.
The network selection method provided by the present disclosure includes acquiring the dynamic network model and the dynamic user model, generating the random event vector according to the dynamic network model and the dynamic user model, constructing the first selection model according to the random event vector and the action vector, obtaining the value of the action probability according to the first selection model, and determining the networks that the users choose to access according to the value of the action probability. It solves the problem that the existing game-theory-based approaches cannot solve the network selection problem for the integrated cellular and drone-cell networks because the integrated cellular and drone-cell networks are highly dynamic, and the network state is hard to predict.
In a second aspect, the present disclosure provides a network selection apparatus for integrated cellular and drone-cell networks, including:
a transceiver, configured to collect capacity information of the drone-cell networks, capacity information of the cellular network, set information of accessible networks of the users, data transmission rate information, and transmitting action vector information to the users so that the users can determine the networks to access according to the action vector information;
a processor, configured to generate the action vector information according to the capacity information of the drone-cell networks, the capacity information of the cellular network, the set information of the accessible networks of the users, the data transmission rate information and a fourth selection model.
The fourth selection model is that a difference value between drift of total violation and an utility is less than or equal to a penalty upper bound, where
the drift of the total violation is obtained according to a value of the total violation at a current time slot and a value of the total violation at a next time slot;
the value of the total violation at the current time slot is obtained according to a first virtual value at the current time slot, a second virtual value at the current time slot and a third virtual value at the current time slot;
the first virtual value in the first virtual queue at the current time slot is generated according to violation of second coarse correlated equilibrium constraints at a previous time slot and the first virtual value in the first virtual queue at the previous time slot;
the second virtual value in the second virtual queue at the current time slot is generated according to violation of third auxiliary variable constraints at the previous time slot and the second virtual value in the second virtual queue at the previous time slot;
the third virtual value in the third virtual queue at the current time slot is generated according to violation of a second minimum individual time average utility constraint at the previous time slot and the third virtual value in the third virtual queue at the previous time slot, where the first virtual value at an initial time slot, the second virtual value at the initial time slot, and the third virtual value at the initial time slot are all zero;
a third selection model includes a third objective function and a third constraint, where the third objective function is time average expectation of a proportional fair function with third auxiliary variables as independent variables, and the third constraint includes at least the second coarse correlated equilibrium constraints, the second minimum individual time average utility constraint, a second auxiliary variable constraint, and the third auxiliary variable constraints, where the second coarse correlated equilibrium constraints are used for constraining time average expectation of the individual utilities and time average expectation of the second auxiliary variables, the second minimum individual time average utility constraint is used for constraining the time average expectation of the individual utilities, the second auxiliary variable constraint is used for constraining the second auxiliary variables, and the third auxiliary variable constraints are used for constraining the time average expectation of the third auxiliary variables and the average time expectation of the individual utilities;
the first selection model includes a first objective function and a first constraint, where the first objective function is a proportional fairness function with time average of the individual utilities as independent variables, and the first constraint includes at least a first coarse correlated equilibrium constraint, a first minimum individual time average utility constraint and a first action probability constraint, where the first coarse correlated equilibrium constraint is used for constraining the time average of the individual utilities and first auxiliary variables, the first minimum individual time average utility constraint is used for constraining the time average of the individual utilities, and the first action probability constraint is used for constraining the action probability under the condition of an random event vector;
the time average of the individual utilities is obtained according to the individual utilities, a random event probability and the action probability under the condition of the random event vector, where the action probability under the condition of the random event vector is the probability that the users execute the action vector under the condition that the random event vector occurs;
an individual utility of each user is obtained according to the action vector and the random event vector; and
the action vector is generated according to the random event vector, and the random event vector is generated according to a capacity model of the cellular network, a capacity model of the drone-cell networks, accessible network sets of the users and a transmission rate model, where the accessible network sets of the users is generated according to a location model of the drone-cell networks and a location model of the users.
In the network selection method and the apparatus provided by the present disclosure, the network selection method includes acquiring the dynamic network model and the dynamic user model, generating the random event vector according to the acquired dynamic network model and the dynamic user model, constructing the first selection model according to the random event vector and the action vector, obtaining the value of the action probability according to the first selection model, and determining the networks that the users choose to access according to the value of the action probability.
The present disclosure constructs a dynamic network model and a dynamic user model, simulating the high dynamic of the drone-cell-user connection, the fluctuation of the network capacity and the time-varying of user traffic, etc. The present disclosure formulates the network selection problem as a problem of a repeated stochastic game, which simulates the competition and interaction among the users well. The method can maximize the total user utility while ensuring the fairness among the users. It solves the problem that the existing game-theory-based approaches cannot solve the network selection problem for the integrated cellular and drone-cell networks because the integrated cellular and drone-cell networks are highly dynamic, and the network state is hard to predict.
In order to make purpose, technical solutions and advantages of embodiments of the present disclosure clearer, the technical solutions of the embodiments of the present disclosure will be described hereunder clearly and comprehensively with reference to accompanying drawings. Apparently, the described embodiments are only a part of embodiments of the present disclosure, rather than all of them. Based on the embodiments of the present disclosure, all other embodiments obtained by those of ordinary skilled in the art without any creative effort shall fall into the protection scope of the present disclosure.
S101: acquiring a dynamic network model and a dynamic user model.
More specifically, the dynamic network model includes at least a location model of drone-cell networks, a capacity model of a cellular network and a capacity model of the drone-cell networks, and the dynamic user model includes at least a location model and a transmission rate model of users;
For the network coverage characteristic, this embodiment assumes that the cellular network can cover the whole considered area and a drone-cell can only cover a small region. Specifically, every drone-cell is assumed to have the same and limited coverage radius (denoted by Rd). Let rij denotes a horizontal distance between the i-th user and the j-th drone-cell. The i-th user can access to the j-th drone-cell, if rij≤Rd; otherwise, it cannot.
For the location model of the drone-cell networks, this embodiment introduces a smooth turn mobility model with reflection boundary. In this model, each drone will fly with a smooth and random trajectory without sharp turning. Specifically, a drone is assumed to have a constant forward speed Vd (in m/s) and can change its centripetal acceleration randomly. The duration (in seconds) for the drone to maintain its current centripetal acceleration is subject to an exponential distribution with mean 1/λd. Meanwhile, the reciprocal of the turning radius (in m) of the drone-cell is subject to a Gaussian distribution with zero-mean and variance σd2.
For the capacity model of the networks, this embodiment assumes that the network capacity, denoted by C1(t), of the cellular network at time slot t is subject to a truncated Gaussian distribution N(μb,σb2) within the interval C1(t)∈[μb−2σb,μb+2σb], where 2σb<μb. Meanwhile, this embodiment assumes that for each drone-cell network j∈, its capacity at slot t, denoted by Cj(t), is independent and identically distributed (i.i.d.) and is subject to a truncated Gaussian distribution N(μc,σc2) within the interval Cj(t)∈[μc−2σc,μc+2σc], where 2σc<μc.
For the location model of the users, this embodiment develops a boundary Gauss-Markov mobility model based on the Gauss-Markov mobility model. Specifically, based on the Gauss-Markov mobility model, this embodiment considers a user that moves within a rectangular area
and reflects on the boundary. Thus, the location (denoted by lu(t)=((xu(t),yu(t))) and the velocity (denoted by vu (t)=((vux(t),vuy(t))) of a user in the boundary Gauss-Markov mobility model can meet the following updated formulas:
denotes the Hadamard product, └⋅┘ denotes the floor operation, αu=(αux,αuy) is a two-dimensional (2-D) memory level vector, σu=(σux,σuy) is a 2-D asymptotic standard deviation vector of the velocity, the 2-D memory level vector and the 2-D asymptotic standard deviation vector of the velocity are constants, wu (t)=(wux(t),wuy(t)) represents a 2-D uncorrelated Gaussian process, and wux(t) and wuy(t) are independent and are both zero-mean and unit-variance.
For the transmission rate model of the users, this embodiment assumes that for each user i∈, the required data transmission rate at time slot t (denoted by Ri(t)) is i.i.d. and is subject to a truncated Gaussian distribution N(μR(t),σR(t)2) within the interval Ri(t)∈[μR (t)−2σR (t),μR(t)+2σR(t)], where σR (t)=ρRμR(t) with ρR<½. Furthermore, μR(t) is a Markov process that taking a value in a finite set {μ1, μ2, . . . , μK
S102: generating accessible network sets of users according to a location model of the drone-cell networks and a location model of the users; and generating a random event vector according to a capacity model of the cellular network, a capacity model of the drone-cell networks, the accessible network sets of the users and the transmission rate model.
More specifically, this embodiment assumes that the proposed network selection is a repeated stochastic game problem.
S103: generating an action vector according to the random event vector.
More specifically, the action vector is used for indicating that the users choose to access to the drone-cell networks and/or the cellular network. After observing the random event vector ω(t) at time slot t, the game manager 202 sends a suggestion si(t)∈(ωi(t)) to each player 203-i, where (ωi(t))=(t)∪{0} is the finite set of actions available to the player 203-i. For example, si(t)=j indicates that the game manager 202 suggests the player 203-i choose the network j. Besides, for convenience of description, this embodiment simplifies (ωi(t)) as (t).
For each player 203-i∈, it will choose an action αi(t)∈(t) based on the suggestion si(t). For example, αi(t)=j indicates that the player 203-i chooses to access to the network j. This embodiment lets s(t)=(s1(t), s2(t), . . . , sN(t)) and α(t)=(α1(t), α2(t), . . . , αN(t)) denote the suggestion vector and the action vector, respectively, and defines (t)=1(t)× . . . ×(t).
S104: obtaining an individual utility of each user according to the action vector and the random event vector.
The random event vector ω(t) and the action vector α(t) at time slot t determine the individual utility ui(t) of each player 203-i. Formally, the individual utility ui(t) can take the following formula:
ui(t)=ûi(α(t),ω(t)) (3)
More explicitly, this embodiment utilizes the following expression to define ui(t).
Definition 1: for all i∈, the individual utility ui(t) can be defined as:
where ƒ(x) is an effective transmission ratio function defined as:
where xb is a constant representing the network busyness ratio threshold. Further,
with 1{αk(t)=αi (t)} denoting a 0-1 indicator function. The indicator function equals one, if αk(t)=αi(t); otherwise, it equals zero. Cα
This embodiment assumes that the upper bound of the required data transmission rate Ri(t) for each player 203-i is uimax. Then, according to Definition 1, ûi(α(t),ω(t)) satisfies the following condition:
0≤ûi(α(t),ω(t))≤uimax (6)
S105: constructing a first selection model.
More specifically, the probability density function (denoted by π[ω]) of the random event vector ω(t) can be defined as the following formula:
π[ω]Pr[ω(t)=ω], ∀ω∈Ω (7)
where the notation “” means “defined as equal to”.
This embodiment further defines the action probability Pr[α|ω] as a conditional probability density function based on α∈ and ω∈Ω, where =1× . . . ×N and
According to the probability theory, the action probability satisfies the first action probability constraint:
where =1× . . . ×N, and i(ωi) represents the finite set of actions available to the player 203-i after the random event ωi is observed.
With the definition of Pr[α|ω], this embodiment defines a variable ūi denoting the time average of the individual utility ui(t). According to the law of large numbers, if the action vector α(t) is chosen independently at each time slot t according to the same conditional probability density function Pr[α|ω], it can be guaranteed that for all i∈, ūi, can be written as the following form with probability 1 (w.p. 1):
In addition, considering that the game manager 202 has an objective of formulating the Pr[α|ω] to maximize the total user utility while ensuring the fairness among users, this embodiment designs an increasing and concave proportional fairness function ϕ(ū1, ū2, . . . , ūN) for the game manager 202 as the first objective function of the first selection model. Explicitly, this embodiment supposes that the proportional fairness function ϕ(ū1, ū2, . . . , ūN) is a sum of logarithmic functions:
ϕ(ū1,ū2, . . . ,ūN)=log2(ūi) (12)
Each player 203-i, however, is interested in maximizing its own time average utility ūi. Thus, players 203 may choose whether to accept the suggestions provided by the game manager 202. For each player 203-i, there are two types of selections as presented below.
Participate:
if a player 203-i always chooses to accept the suggestion si(t) at each time slot t∈{0, 1, 2, . . . ,}, it is called as participate. That is, αi (t)=si(t) for all t∈{0, 1, 2, . . . }.
Non-Participate:
if a player 203-i chooses the action αi(t) according to its observed random event ωi(t) at each time slot t∈{0, 1, 2, . . . ,}, it is called as non-participate.
This embodiment assumes that non-participating players 203 will not receive the suggestions si(t).
Definition 2: in order to incentivize all players 203 to participate, the Pr[α|ω] formulated by the game manager 202 needs to be a coarse correlated equilibrium (CCE), which is defined as follows:
For a stochastic game, Pr[α|ω] is a CCE if there are first auxiliary variables φi(υi)∈[0, uimax] for all i∈ meeting the following conditions:
where, αī=α\{αi} represents all entries in the action vector α except for the αi, βi∈ is a preset specific action of a non-participating player 203-i, and υi∈Ωi is a preset specific event of a non-participating player 203-i. Intuitively, φi(υi) represents the largest conditional expected utility earned by the non-participating player 203-i when ωi=υi is observed.
According to Definition 2, the total number of the CCE constraints (13) and (14) is N+|Ωi∥|, which is a linear function over the size of sets Ωi and , where |⋅| represents the number of elements in a set. The value of |Ωi∥|, however, is too large in the system model of this embodiment; thus, the computation of a CCE is complex. This embodiment next describes how to reduce the value of |Ωi∥|.
First, this embodiment simplifies the value space Ωi of the preset event υi. From the description of the network model in this embodiment, all drone-cells are homogenous. Therefore, for the preset event υi, this embodiment only considers the number of the accessible drone-cell networks of a user i∈, rather than the difference in their indexes. Furthermore, the probability that a user i is covered by more than two drone-cell networks simultaneously is small. Thus, the set of (t) can be simplified as {0, 1, 2+}, where, “0”, “1”, and “2+” indicate that the number of the drone-cell networks covering a user i is zero, one and not less than two, respectively. Besides, for preset event υi, this embodiment divides the interval [0, uimax] into KR segments uniformly. If
where ik=1, 2, . . . , KR, then Ri(t) belongs to the ik-th segment. To sum up, this embodiment can simplify the value space Ωi of the preset event υi as Ωis={0, 1, 2+}×{1, 2, . . . , KR}. Note that the value space of ωi is still Ωi. For all ωi∈Ωi and υi∈Ωis, ωi=υi in the constraint (14) indicates that ωi is the original form of υi, and in the constraint (13), φi(ωi)=φi(υi|υi=ωi) with υi=ωi indicating that υi is the simplified form of ωi.
Second, this embodiment simplifies the value space of the preset action βi. Owing to the homogeneity of drone-cells, this embodiment does not identify the index of the drone-cell network that a user i chooses to access for the preset action βi. Meanwhile, since the individual utility ui(t) of a user i at time slot t equals zero when the user i accesses to the empty network according to Definition 1, this embodiment does not consider the preset action of accessing to the empty network, i.e., βi≠0. Therefore, this embodiment simplifies the value space of the preset action βi as ={cellular, drone-cell}, where “cellular” and “drone-cell” indicate that the user i chooses to access to a cellular network and a drone-cell network, respectively. Besides, when choosing to access to a drone-cell network, the user i accesses randomly to an accessible drone-cell network with equal probability. Note that the value space of αi is still .
Therefore, the value of |Ωi∥| is reduced to |Ωis∥|=(|{0,1,2+}|×KR)×||=6KR. Meanwhile, the preset action βi=drone-cell is infeasible when there is no accessible drone-cell networks for a user i. Thus, this embodiment can ignore the following preset event-preset action pairs: {(0, ik;drone-cell), ∀ik∈{1, 2, . . . , KR} }. In this way, the value of |Ωis∥| is further reduced to 6KR−KR=5KR. Finally, the total number of the CCE constraints (13) and (14) is reduced to N+|Ωis∥|=(5KR+1)N.
In actual scenarios, some users have a requirement of minimum time average utilities (MTAU), and this embodiment denotes the set of these users as Su. Therefore, the game manager 202 must guarantee the utilities of these users meet the following MTAU constraints:
ūi≥uic, ∀i∈Su (15)
Based on the above analysis, the constructed first selection model includes: a first objective function and a first constraint. The first constraint includes: a first coarse correlated equilibrium constraint, a first minimum individual time average utility constraint and a first action probability constraint. The constructed first selection model is as follows:
S106: obtaining a value of an action probability according the first selection model, and determining networks that the users choose to access according to the value of the action probability.
The objective of the game manager 202 is to solve the first selection model to obtain the action probability Pr[α|ω], and to choose suggestion vector s(t)=α(t) according to Pr[α|ω] to determine the networks that the users choose to access according to the suggestion vector.
The network selection method provided by this embodiment includes: acquiring the dynamic network model and the dynamic user model, generating the random event vector according to the dynamic network model and the dynamic user model, constructing the first selection model according to the random event vector and the action vector, obtaining the value of the action probability according to the first selection model to determine the networks that the users choose to access according to the value of the action probability. This embodiment solves the problem that the existing game-theory-based approaches cannot solve the network selection problem for the integrated cellular and drone-cell networks because the integrated cellular and drone-cell networks are highly dynamic, and the network state is hard to predict.
Although the above problem (16) is a convex optimization problem, it is highly challenging to solve it owing to the following two reasons: 1) π[ω] is essential to solve this problem (16); however, it may be impossible to obtain π[ω] because π[ω] is influenced by various factors such as network capacity, the mobility of drones and users, and user traffic. 2) The size of the variable Pr[α|ω] is
which increases exponentially with the increasing number of users. To solve these issues, this embodiment converts this challenging problem into a new problem without knowing π[ω], and the size of the new problem is greatly reduced.
This embodiment illustrates a network selection method according to another exemplary embodiment, which differs from the embodiment illustrated in
S1051: constructing a second selection model.
More specifically, for a real-valued stochastic process u(t) over time slots t∈{0, 1, 2, . . . }, this embodiment defines its time average expectation over the first t time slots as:
For all i∈, υi∈Ωis and βi∈, this embodiment defines:
ui,υ
By the stochastic process theory, this embodiment converts the first selection model (16) equivalently into the second selection model, where the second selection model includes a second objective function and a second constraint, where the second objective function is a proportional fairness function with the time average expectation of the individual utilities as independent variables, and the second constraint includes at least second coarse correlated equilibrium constraints, a second minimum individual time average utility constraint and a second auxiliary variable constraint, where the second coarse correlated equilibrium constraints are used for constraining the time average expectation of the individual utilities and the time average expectation of second auxiliary variables, the second minimum individual time average utility constraint is used for constraining the time average expectation of the individual utilities, and the second auxiliary variable constraint is used for constraining the second auxiliary variables. At each time slot t∈{0, 1, 2, . . . }, the game manager 202 observes the random event vector ω(t)∈Ω, and solves an action vector α(t)∈(t) and variables φi,υ
S106: obtaining a value of the action vector according to the second selection model, and determining the networks that the users choose to access according to the value of the action vector.
In the network selection method provided in this embodiment, the constructed second selection model is based on the time average expectation of the individual utilities, the action vector can be obtained when the probability of the random event vector is unknown, and the networks that the users choose to access are determined according to the value of the action vector.
This embodiment illustrates a network selection method according to another exemplary embodiment, which differs from the previous embodiment in that, after S1051, it further includes:
S1052: constructing a third selection model.
More specifically, the objective of the above problem (19) is to maximize a nonlinear function of time average. In order to convert it equivalently into a maximization of the time average of a nonlinear function, this embodiment introduces a third auxiliary vector γ(t)=(γ1(t), . . . , γN(t)) with 0≤γi(t)≤uimax for all i∈, and defines g(t)=ϕ(γ1(t), . . . , γN(t). Then, the Jensen's inequality indicates that
This embodiment considers converting the second selection model (19) into a third selection model (21.1) by leveraging the Jensen's inequality. Where the third selection model includes a third objective function and a third constraint, and the third objective function (21.1) is time average expectation of a proportional fairness function with the third auxiliary variables as independent variables, and the third constraint includes at least the second coarse correlated equilibrium constraints (21.4) and (21.5), the second minimum individual time average utility constraint (21.7), the second auxiliary variable constraint (21.6) and third auxiliary variable constraints (21.2) and (21.3). At each time slot t∈{0, 1, 2, . . . }, the game manager 202 observes the random event vector ω(t)∈Ω and chooses an action vector α(t)∈(t), variables φi,υ
where, the second coarse correlated equilibrium constraints (21.4) and (21.5) are used to constrain the time average expectation of the individual utilities and the time average expectation of the second auxiliary variables, the second minimum individual time average utility constraint (21.7) is used to constrain the time average expectation of the individual utilities, the second auxiliary variable constraint (21.6) is used to constrain the second auxiliary variables, and the third auxiliary variable constraints (21.2) and (21.3) are used to constrain the time average expectation of the third auxiliary variables and the time average expectation of the individual utilities.
S106: obtaining the value of the action vector according to the third selection model, and determining the networks that the users choose to access according to the value of the action vector.
In the network selection method provided in this embodiment, the time average expectation of the proportional fairness function with the third auxiliary variables as independent variables is taken as the third objective function, so that the objective function can be simplified to obtain the action vector according to the third selection model.
This embodiment illustrates a network selection method according to another exemplary embodiment, which differs from the previous embodiment in that, after S1052, it includes:
S1053: converting the third selection model into a fourth selection model by leveraging a drift-plus-penalty technique.
Considering the principle of the drift-plus-penalty technique, as for the constraint (21.1), this embodiment defines the first term Qi(t) of the first virtual queue for all i∈ as the following form:
The constraint (21.4) can be satisfied, if the following mean-rate stability condition is held:
where, the nonnegative operation [x]+=max{x,0}.
Likewise, to enforce the constraints (21.1), (21.1) and (21.1), this embodiment defines other three types of virtual queues, respectively. The second item Di,υ
Di,υ
The second virtual queue Zi(t) is defined for all i∈:
Zi(t+1)=Zi(t)+γi(t)−ui(t) (25)
The third virtual queue Hi(t) is defined for all i∈Su:
Hi(t+1)=Hi+uic−ui(t) (26)
And the constraints (21.1), (21.1) and (21.1) can be satisfied, if the following mean-rate stability conditions are held:
For simplicity, this embodiment assumes that all virtual queues are initialized to zero.
According the formulas (22), (24), (25) and (26), the first virtual value in the first virtual queue at the current time slot is generated according to the violation of the second coarse correlated equilibrium constraints at the previous time slot and the first virtual value in the first virtual queue at the previous time slot. The second virtual value in the second virtual queue at the current time slot is generated according to the violation of the third auxiliary variable constraints at the previous time slot and the second virtual value in the second virtual queue at the previous time slot. The third virtual value in the third virtual queue at the current time slot is generated according to the violation of the second minimum individual time average utility constraint at the previous time slot and the third virtual value in the third virtual queue at the previous time slot.
This embodiment defines a function L(t) as a sum of squares of all four types of queues [Qi(t)]+, [Di,υ
where, Hi(t)=0 for all i∉Su.
Besides, this embodiment defines a drift-plus-penalty expression as Δ(t)−Vg(t), where Δ(t)=L(t+1)−L(t) represents a Lyapunov drift, that is, the drift value of the total violation, −g(t) is a “penalty”, g(t) is the proportional fairness function with the third auxiliary variables as independent variables, and V is a non-negative penalty coefficient that affects a trade-off between the constraint violation and the optimality. The drift-plus-penalty expression satisfies the following condition: minimizing the constraint violation and maximizing the objective. Therefore, the fourth selection model is constructed as follows:
where a penalty upper bound includes: a constant term, a first penalty upper bound term, a second penalty upper bound term and a third penalty upper bound term. The constant term is
the first penalty upper bound term is
the second penalty upper bound term is
and the third penalty upper bound term is
S106: obtaining the value of the action vector according the fourth selection model, and determining the networks that the users choose to access according to the value of the action vector.
In the network selection method provided in this embodiment, the fourth selection model is an inequality form, and this model is simple and easy to obtain the action vector according to the fourth selection model.
This embodiment illustrates a network selection method according to another exemplary embodiment, which differs from the previous embodiment in that: S106: obtaining the value of the action vector according the fourth selection model, and determining the networks that the users choose to access according to the value of the action vector, specifically including the following steps:
S1061: at each time slot t, the game manager 202 observes the first term Qi(t) of the first virtual value in the first virtual queue at the current time slot, the second term Di,υ
This embodiment solves the problem (21.1) by minimizing the upper bound of Δ(t)−Vg(t) at each time slot t greedily. Meanwhile, the upper bound of Δ(t)−Vg(t) can be decomposed into four independent terms. At each time slot t, the first term is a constant, the second term is a function of the third auxiliary vector γ(t), the third term is a function of the second auxiliary variables φi,υ
S1062: choosing the value of the third virtual variables γi(t) for all i∈ according to the second virtual value at the current time slot and the first penalty upper bound term to solve:
The closed-form solution of the problem (32) can take the following form for all i∈:
S1063: choosing the value of the second virtual variables φi,υ
The closed-form solution of the problem (34) can take the following form for all i∈and υi ∈Ωis:
S1064: choosing the value of the action vector α(t) according to the random event vector, the first virtual value at the current time slot, the second virtual value at the current time slot, the third virtual value at the current time slot and the third penalty upper bound term to solve:
S1065: sending αi(t) to each player 203-i, so that the player 203-i determines the network to choose to access according to the action αi(t).
The individual utilities ui(t) and ui,υ
The problem (36) is a non-linear integer programming problem, where ui(t) and ui,υ
which increases exponentially with the increasing number of users. Although heuristic algorithms (e.g., genetic algorithm) can be leveraged to mitigate this problem, it may require a long processing time due to its slow convergence rate. To accelerate the optimization process, this embodiment designs a linear approximation mechanism for the problem (36).
According to Definition 1, the network j∈\{0} may be congested when the following condition is held: 1{αi(t)=j}Rι(t)>xbCj (t). To avoid this situation, the suggested action vector α(t) formulated by the game manager 202 should satisfy the first action vector constraint:
For each participating player 203-i∈, according to Definition 1, if αi(t)=0, ui(t)=0. According to both the constraint (37) and Definition 1, if αi(t)≠0, ui(t)=Ri(t). Therefore, the utility function ui(t) in both cases can be calculated in the following way, thereby forming a mapping table between the individual utilities and the transmission rate of the participating players:
ui(t)=1{αi(t)≠0}Ri(t) (38)
For each non-participating player 203-i∈, and for each υi ∈Ωis and βi∈is, (βi≠0), this embodiment considers the definition (18) of ui,υ
1) If the player 203-i accesses to a network that is just the one suggested by the game manager 202, i.e., βi=αi (t), then ui,υ
2) If the player 203-i accesses to a network that is not the one suggested by the game manager 202, i.e., βi≠αi(t), then this embodiment estimates the effective transmission ratio of the network j=βi at time slot t (denoted by θi,β
where,
Cβ
Therefore, ui,υ
Next, this embodiment discusses how to transform the problem (36) into an integer linear programming problem by introducing a set of auxiliary variables aij, where {aij} is a suggestion matrix. For all i∈ and j∈, this embodiment defines the mapping relationship between the suggestion matrix and the action vector as:
where aij=1 indicates that the game manager 202 suggests the player 203-i access to the network j, and aij=0 indicates that the game manager 202 suggests the player 203-i not access to the network j. According to the definition of aij and the constraint α(t)∈(t), the suggestion matrix constraint can be obtained:
α(t) in (37) is substituted with the variables aij, and the second action vector constraint can be obtained:
α(t) in (37) is substituted with the variables aij, and the individual utilities of the participating players can be obtained:
Furthermore, this embodiment lets
for all i∈, j∈ and βi∈. α(t) in (40) is substituted with the variables aij, and the individual utilities of the non-participating players can be obtained:
ui,υ
According to (41)-(46), this embodiment can transform the problem (36) into the following integer linear programming problem:
where, the weights cij are defined as:
is the fourth penalty upper bound term,
constitutes the second action vector constraint, and
Note that at the initial time slot (t=0), all weights cij will be zero since all virtual queues are initialized to zero. To handle this problem, this embodiment defines the weights cij at time slot t=0 as:
The problem (47) is an integer linear programming problem with respect to the auxiliary variables aij that can be solved by the MOSEK Optimization Tools. Furthermore, in the MOSEK, the branch-and-bound scheme is leveraged to relax the integer variables; thus, the integer linear optimization problem is relaxed to a solvable linear optimization problem.
This embodiment uses the Lyapunov optimization method and the linear approximation mechanism to transform the constructed problem without knowing the state probability π[ω] of the networks and the users, thereby reducing the computational complexity of the problem greatly.
Based on the main framework of solving the problem (21), and combined with the solutions of (32), (34) and (36), this embodiment proposes an network selection apparatus for the integrated cellular and drone-cell networks based on the efficient and fair network selection (EFNS) method shown in
Referring further to
Referring further to
The state information ω(t) and the virtual queue information Qi(t), Di,υ
S201: obtaining an upper bound uimax of data transmission rate, a segment value KR, a penalty coefficient V, and initializing a first virtual queue to a third virtual queue.
More specifically, these parameters are stored in the memory 312. These parameters can be given default values, and the computer operator 340 can modify these parameters through the human-computer interaction module 314. The virtual queues Qi(0)=0, Di,υ
Steps 2-7 are repeated at each time slot, where T is a total number of the time slots.
S202: collecting state information of networks and users to form a random event vector.
More specifically, the processor 313 collects the state information ω(t)∈Ω of the networks and the users through the transceiver 311.
Specifically, ω(t) will be stored in the memory 312 temporarily until the end of Step 6.
S203: obtaining third auxiliary variables γi(t) according to a second virtual value at a current time slot and a first penalty upper bound term.
More specifically, for each i∈, the processor 313 calculates the third auxiliary variables γi(t) according to the formula (33). The third auxiliary variables γi(t) will be temporarily stored in the memory 312 until the end of Step 6.
S204: obtaining second auxiliary variables according to the random event vector at the current time slot, a first virtual value at the current time slot and a second penalty upper bound term.
More specifically, for each i∈ and υi ∈Ωis, the processor 313 calculates the second auxiliary variables φi,υ
S205: obtaining a suggestion matrix according to a fourth penalty upper bound, a second action vector constraint and a suggestion matrix constraint, and obtaining an action vector according to the suggestion matrix.
More specifically, the processor 313 obtains the suggestion matrix {aij} by solving the problem (47), and obtains the action vector α(t) according to the suggestion matrix {aij}. Then the transceiver 311 sends the suggested actions αi(t) to the network access system 331 of each user device 330-i. The suggested action vector α(t) will be temporarily stored in the memory 312 until the end of Step 6.
S206: calculating individual utilities and the first virtual queue to the third virtual queue.
More specifically, the processor 313 calculates ui(t) and ui,υ
S207: judging whether a time slot t reaches a preset time slot value, if not, turning to S208, otherwise, stopping the loop.
S208: updating the time slot t, and turning to S201.
The following context presents a simulation leveraging the EFNS method provided by this embodiment to perform the network selection for the integrated cellular and drone-cell networks.
In order to verify the effectiveness of the network selection method, this embodiment designs three benchmark comparison methods: the cellular-only (CO) method, the random access (RA) method and the on-the-spot offloading (OTSO) method. For the CO method, at each time slot, every user always chooses to access to the cellular network. For the RA method, at each time slot, every user always accesses randomly to an accessible non-empty network with equal probability. For the OTSO method, at each time slot, every user checks whether the drone-cell networks can be access, and if so, the user accesses randomly to an accessible drone-cell network with equal probability; otherwise, the user accesses to the cellular network.
The parameter setting in the simulation is summarized as the following: the size of the considered geographic area is 500×500 m2, i.e., L=500 m and V=500 m. The repeated stochastic game lasts for 1000 seconds and let the duration of a time slot be one second; thus, the simulation runs for 1000 episodes, i.e., T=1000. In the location model of the drone-cell networks, the initial locations of the drones are distributed in the considered area independently and uniformly, and their initial heading angles are independent and subject to a uniform distribution on [0,2π), and parameters (Vd, λd, σd2)=(10, 0.1, 0.02). In the network capacity model, the capacity (in Mb/s) of the cellular network is subject to a truncated Gaussian distribution Ntru(200, 202, ±40), and the capacity (in Mb/s) of each drone-cell network is independent and subject to a truncated Gaussian distribution Ntru(30, 32, ±6). Furthermore, the coverage radius of a drone-cell Rd=100 m.
In the location model of the users, the initial locations of the users are distributed in the considered area independently and uniformly, and the initial velocities of the users are independent and subject to 2-D Gaussian distribution N(0, 0; 22, 22, 0), where parameters αu=(0.73,0.73) and σu=(2,2). In the transmission rate model of the users, a parameter ρR=0.2 and the process μR(t) (in Mb/s) takes a value from the set {μ1, μ2, . . . , μ5 }={2.5, 5, 7.5, 10, 12.5}. Besides, the one-step transition probability matrix P of μR(t) is shown in Table 1:
For the function ƒ(x) in Definition 1, this embodiment sets the network busyness ratio threshold xb=0.9. For the MTAU constraints, this embodiment sets Su={1, . . . , 10} and sets uic=6 for all i∈Su. For the network selection method provided by the present disclosure, set uimax=20 for all i∈, KR=5, and the penalty coefficient V=100.
Meanwhile, this embodiment leverages the following four indexes for performance evaluation of the following proposed method, including:
Queue Stability: this embodiment uses the stability variables defined as
over time slots t=1, 2, . . . , T−1 to measure the stability of the queues of the EFNS method.
Running Time: this is the total time of executing the EFNS method for T=1000 episodes.
Total Utility: this is the total utility of all users during the entire simulation process,
Fairness: this embodiment uses the Jain's fairness index, defined as
to measure the fairness of the network resource allocation, where ūi represents the time average utility of a user i during the entire simulation process, i.e.,
In the simulation, this embodiment tests all comparison methods on one hundred randomly generated data sets. For each comparison method, this embodiment may obtain one hundred results, and the final result is their average.
all stability variables decrease rapidly with the increase of time slot t. After a long period of time, all stability variables tend to be zero. This result indicates that the EFNS method can guarantee that all queues are mean-rate stable, and thus the constraints (21.1), (21.1), (21.1) and (21.1) are satisfied.
the average running time of the EFNS method increases with the increase of N or Md. This is because when N or Md increases, the scale of the problem becomes larger. Meanwhile, the EFNS method can achieve an online network selection.
the EFNS method can always achieve the highest total utility compared with the other three methods by avoiding congestion and making full use of network resources.
For the EFNS method, the total utility thereof increases with the increase of N. However, the rate of increase slows down because network capacity limit the increase of the total utility when N is large. For the other three methods, these total utilities thereof begin quickly to decrease with the increase of N. This is because these three methods do not have a mechanism to avoid network congestion and a large number of users may lead to network congestion.
For all methods except for the CO method, their total utilities increase monotonically with the increase of Md, since users can offload traffic to drone-cell networks.
the EFNS method can achieve a high level of fairness. Specifically, its fairness index is close to 1 because the fair allocation of network resources is considered in the proportional fairness function (12). However, the fairness index of the EFNS method decreases gradually with the increase of N or the decrease of Md, that is because a great N or a small Md may activate the MTAU constraints.
Although the RA and OTSO methods do not consider fairness, they can also achieve a high level of fairness. This is because the time average utilities of users are close to each other after a long time due to the homogeneity of users in the model of this embodiment.
The CO method can achieve the highest level of fairness since users always have the same effective transmission ratio at each time slot. This CO method, however, achieves the lowest total utility.
Finally, it should be noted that each of the above embodiments is merely intended to describe, rather than limit, the technical solutions of the present disclosure. Although the present disclosure is described in detail with reference to the foregoing embodiments, persons of ordinary skill in the art should understand that it is possible to make modifications to the technical solutions described in the foregoing embodiments, or make equivalent substitutions of part or all technical features therein. However, these modifications or substitutions do not make the essence of corresponding technical solutions depart from the scope of the technical solutions in the embodiment solutions of the present disclosure.
Number | Date | Country | Kind |
---|---|---|---|
201811353219.4 | Nov 2018 | CN | national |
Number | Name | Date | Kind |
---|---|---|---|
9866313 | Murphy | Jan 2018 | B1 |
20170111228 | Obaidi | Apr 2017 | A1 |
20180007518 | O'Berry et al. | Jan 2018 | A1 |
20180102831 | Murphy | Apr 2018 | A1 |
20180276995 | Priest | Sep 2018 | A1 |
20190035287 | Priest | Jan 2019 | A1 |
20190295427 | Priest | Sep 2019 | A1 |
20190327712 | Moradi | Oct 2019 | A1 |
20200202284 | Singh | Jun 2020 | A1 |
20200223543 | Singh | Jul 2020 | A1 |
20200260404 | Vrind | Aug 2020 | A1 |
20200314628 | Panchal | Oct 2020 | A1 |
20200355571 | Priest | Nov 2020 | A1 |
20210044613 | Sugarbaker | Feb 2021 | A1 |
Number | Date | Country |
---|---|---|
103332290 | Oct 2013 | CN |
106054928 | Oct 2016 | CN |
106604205 | Apr 2017 | CN |
107517158 | Dec 2017 | CN |
Entry |
---|
Huang, Hailong, “An Algorithm of Efficient Proactive Placement of Autonomous Drones for Maximum Coverage in Cellular Networks” IEEE Wireless Communications Letters; vol. 7, Issue 6; (Jun. 2018); pp. 994-997. |
Xu, Min et al., “Research of UAV Communications Based on Mobile Cellular Networks” Chinese Journal of Mobile Communication; No. 22; (Dec. 2017); pp. 23-28, English Abstract is included. |
The Chinese First Examination Report of correspondence Chinese application No. 201811353219.4, dated Dec. 4, 2019. |
Notice of Allowance of prior Chinese application No. 201811353219.4. |
Number | Date | Country | |
---|---|---|---|
20200154349 A1 | May 2020 | US |