Network selection method and apparatus for integrated cellular and drone-cell networks

Information

  • Patent Grant
  • 11102717
  • Patent Number
    11,102,717
  • Date Filed
    Tuesday, November 5, 2019
    5 years ago
  • Date Issued
    Tuesday, August 24, 2021
    3 years ago
Abstract
Provided are a network selection method and an apparatus for integrated cellular and drone-cell networks, where the method includes: acquiring a dynamic network model and a dynamic user model, generating a random event vector according to the dynamic network model and the dynamic user model; generating an action vector according to the random event vector; obtaining an individual utility of each user according to the action vector and the random event vector; constructing a first selection model; obtaining the value of the action probability according the first selection model to determine networks that the users choose to access according to the value of the action probability.
Description
CROSS-REFERENCE TO RELATED APPLICATION

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.


TECHNICAL FIELD

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.


BACKGROUND

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.


SUMMARY

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.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a scenario diagram of a network on which a network selection method according to the present disclosure is based;



FIG. 2 is a flowchart illustrating a network selection method according to an exemplary embodiment of the present disclosure;



FIG. 3 is a schematic diagram of a repeated stochastic game structure obeyed by the network selection method according to the embodiment shown in FIG. 2;



FIG. 4 is a schematic structural diagram of a network selection apparatus for the integrated cellular and drone-cell networks according to an exemplary embodiment of the present disclosure;



FIG. 5 is a flowchart of the method performed by a processor in the network selection apparatus for the integrated cellular and drone-cell networks according to the embodiment shown in FIG. 4 of the present disclosure;



FIG. 6 is a schematic diagram illustrating the change of stability variables with time of the network selection method proposed by the present disclosure when the number of users is N=50 and the number of drones is Md=6;



FIG. 7 is a schematic diagram illustrating the effect of the number of users N on the running time of the network selection method proposed by the present disclosure when the number of drones is Md=6;



FIG. 8 is a schematic diagram illustrating the effect of the number of drones Md on the running time of the network selection method proposed by the present disclosure when the number of users is N=50;



FIG. 9 is a schematic diagram illustrating the effect of the number of users N on the total user utilities obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of drones is Md=6;



FIG. 10 is a schematic diagram illustrating the effect of the number of drones Md on the total user utilities obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of users is N=50;



FIG. 11 is a schematic diagram illustrating the effect of the number of users N on the Jain's fairness index obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of drones is Md=6; and



FIG. 12 is a schematic diagram illustrating the effect of the number of drones Md on the Jain's fairness index obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of users is N=50.





DETAILED IMPLEMENTATION OF THE EMBODIMENTS

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.



FIG. 1 is a scenario diagram of a network on which a network selection method according to the present disclosure is based. As shown in FIG. 1, the present disclosure investigates a network selection scenario for the integrated cellular and drone-cell networks. In this scenario, a cellular network 103 is constructed to provide radio access to a collection of users custom character={1, 2, . . . , N} moving randomly and independently in a geographical area of L×W m2 over an infinite sequence of time slots t={0, 1, 2, . . . }. Meanwhile, a set of drone-cells are deployed to alleviate the congestion situation of the cellular network. Meanwhile, the present disclosure assumes that every drone-cell connects directly to a ground station moves independently and is deployed at the same and fixed altitude h. Let custom character={0, 1, 2, . . . , M} be a network set, where j=1 denotes the cellular network, j∈custom characterd={2, . . . , M} denotes the j-th drone-cell network, and j=0 represents an empty network indicating that a user does not access to any network. At each time slot t, a user can select a network to access from the accessible network set. Meanwhile, the present disclosure does not consider the switching cost when a user changes its network access state.



FIG. 2 is a flowchart illustrating a network selection method according to an exemplary embodiment of the present disclosure. As shown in FIG. 2, the network selection method according to the present disclosure includes:


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(μbb2) within the interval C1(t)∈[μb−2σbb+2σb], where 2σbb. Meanwhile, this embodiment assumes that for each drone-cell network j∈custom character, 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(μcc2) within the interval Cj(t)∈[μc−2σcc+2σc], where 2σcc.


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







[


-

L
2


,

L
2


]

×

[


-

W
2


,

W
2


]






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:











l
u



(

t
+
1

)


=



(

-
1

)

k



(



l
u



(
t
)


+


v
u



(
t
)


-

k


(

L
,
W

)



)






(
1
)









v
u



(

t
+
1

)


=




(

-
1

)

k



α
u




v
u



(
t
)



+



1
+

α
u
2





σ
u




w
u



(
t
)










where







k
=


(


k
x

,

k
y


)

=

(





(



x
u



(
t
)


+


v
u
x



(
t
)


+

L
2


)

/
L



,








(



y
u



(
t
)


+


v
u
y



(
t
)


+

W
2


)

/
W




)



,






(
2
)








denotes the Hadamard product, └⋅┘ denotes the floor operation, αu=(αuxuy) is a two-dimensional (2-D) memory level vector, σu=(σuxuy) 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∈custom character, 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μ}, and its one-step transfer probability matrix is defined by this embodiment as P=[piμ,jμ]Kμ×Kμ, where piμ,jμ is the probability that the μR(t) will transfer from the given current value μiμ to the value μjμ at next time slot.


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. FIG. 3 is a schematic diagram of a repeated stochastic game structure obeyed by the network selection method according to the embodiment shown in FIG. 2. As shown in FIG. 3, this structure includes an environment 201, a game manager 202 and players 203. This embodiment regards N users as players 203 and network selection strategies taken by the users as actions. At each time slot t∈{0, 1, 2, . . . ,}, each player 203-icustom character can observe a random event ωi (t)∈Ωi from the environment 201 and the game manager 202 can observe the whole random event vector ω(t)=(ω0(t), ω1(t), . . . , ωN(t))∈Ω from the environment 201, where ω0(t)∈Ω0 represents the random event known only by the game manager 202 and Ω=Ω0×Ω1× . . . ×ΩN. Specifically, in this embodiment, the random event ω0(t) known only by the game manager 202 includes the capacity of the cellular network and the capacity of the drone-cell networks, i.e., ω0 (t)=(C1(t), C2(t), . . . , CM(t)). For all i∈custom character, set ωi(t)={custom character(t),Ri(t)}, where custom character(t) represents the set of accessible non-empty networks of the player 203-i, i.e., custom character(t)={j∈custom character|rij(t)≤Rd}∪{1}.


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)∈custom characteri(t)) to each player 203-i, where custom characteri(t))=custom character(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 custom characteri(t)) as custom character(t).


For each player 203-icustom character, it will choose an action αi(t)∈custom character(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 custom character(t)=custom character1(t)× . . . ×custom character(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∈custom character, the individual utility ui(t) can be defined as:











u
i



(
t
)


=




u
^

i



(


α


(
t
)


,

ω


(
t
)



)


=

{




0
,






α
i



(
t
)


=
0









R
i



(
t
)




f


(


η
i



(
t
)


)



,



otherwise









(
4
)








where ƒ(x) is an effective transmission ratio function defined as:










f


(
x
)


=

{




1
,




0

x


x
b









1


(

1
-

x
b


)

2





(

1
-
x

)

2


,





x
b

<
x

1






0
,



otherwise








(
5
)








where xb is a constant representing the network busyness ratio threshold. Further,








η
i



(
t
)


=






k









1


{



α
k



(
t
)


=


α
i



(
t
)



}




R
k



(
t
)





C


α
i



(
t
)





0






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αi(t) represents the capacity of the network corresponding to αi(t).


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:

π[ω]custom characterPr[ω(t)=ω], ∀ω∈Ω  (7)


where the notation “custom character” means “defined as equal to”.


This embodiment further defines the action probability Pr[α|ω] as a conditional probability density function based on α∈custom character and ω∈Ω, where custom character=custom character1× . . . ×custom characterN and







i

=






ω
i



Ω
i






i



(

ω
i

)



=

𝒥
.







According to the probability theory, the action probability satisfies the first action probability constraint:











Pr


[

α
|
ω

]



0

,



α



,

ω

Ω





(
8
)








Pr


[

α
|
ω

]


=
0

,



α




(
ω
)




,

ω

Ω





(
9
)










Pr


[

α
|
ω

]



=
1

,



ω

Ω






(
10
)







where custom character=custom character1× . . . ×custom characterN, and custom characterii) 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∈custom character, ūi, can be written as the following form with probability 1 (w.p. 1):










u
i

=




ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(

α
,
ω

)








(
11
)







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:

ϕ(ū12, . . . ,ūN)=custom characterlog2(ū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 φii)∈[0, uimax] for all i∈custom character meeting the following conditions:


















ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(

α
,
ω

)









ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]





φ
i



(

ω
i

)









(
13
)













ω

Ω

|

ω
i


=

υ
i










π


[
ω
]




Pr


[

α
|
ω

]





φ
i



(

υ
i

)











ω

Ω

|

ω
i


=

υ
i










π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(


(


β
i

,

α

i
_



)

,
ω

)



















υ
i



Ω
i



,


β
i



i







(
14
)







where, αī=α\{αi} represents all entries in the action vector α except for the αi, βicustom character 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, φii) represents the largest conditional expected utility earned by the non-participating player 203-i when ωii is observed.


According to Definition 2, the total number of the CCE constraints (13) and (14) is N+custom charactericustom character|, which is a linear function over the size of sets Ωi and custom character, where |⋅| represents the number of elements in a set. The value of |Ωicustom character|, 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 |Ωicustom character|.


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∈custom character, 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 custom character(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









R
i



(
t
)




[




(


i
k

-
1

)



u
i

ma





x




K
R


,



i
k



u
i

ma





x




K
R



)


,





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, ωii in the constraint (14) indicates that ωi is the original form of υi, and in the constraint (13), φii)=φiiii) with υii indicating that υi is the simplified form of ωi.


Second, this embodiment simplifies the value space custom character 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 custom character of the preset action βi as custom character={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 custom character.


Therefore, the value of |Ωicustom character| is reduced to |Ωiscustom character|=(|{0,1,2+}|×KR)×|custom character|=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 |Ωiscustom character| is further reduced to 6KR−KR=5KR. Finally, the total number of the CCE constraints (13) and (14) is reduced to N+custom characteriscustom character|=(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:




















Maximize







Pr


[

α
|
ω

]





φ
i



(

υ
i

)










ϕ


(



u
_

1

,





,


u
_

N


)













subject





to


:













u
_

i


=




ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(

α
,
ω

)





,



i




















ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(

α
,
ω

)









ω

Ω









π


[
ω
]




Pr


[

α
|
ω

]





φ
i



(

ω
i

)





,











i

















ω

Ω

|

ω
i


=

υ
i










π


[
ω
]




Pr


[

α
|
ω

]





φ
i



(

υ
i

)











ω

Ω

|

ω
i


=

υ
i










π


[
ω
]




Pr


[

α
|
ω

]






u
^

i



(


(


β
i

,

α

i
_



)

,
ω

)


















i




,


υ
i



Ω
i
s


,


β
i



i
s
















u
_

i



u
i
c


,



i


𝒮
u
















Pr


[

α
|
ω

]



0

,



α



,

ω

Ω














Pr


[

α
|
ω

]


=
0

,



α




(
ω
)




,

ω

Ω
















Pr


[

α
|
ω

]



=
1

,



ω

Ω








(
16
)







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










i












i








i





{
0
}











Ω
i






,





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 FIG. 2 in that, after S104, it includes:


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:











u
_



(
t
)


=


1
t






τ
=
0


t
-
1




E


[

u


(
τ
)


]








(
17
)







For all i∈custom character, υi∈Ωis and βicustom character, this embodiment defines:

ui,υii)(t)custom characterûi((βiī(t)),ω(t))1{ωi(t)=υi}  (18)


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)∈custom character(t) and variables φi,υi(t), where φi,υi(t) represent the second auxiliary variable at time slot t.












Maximize







α


(
t
)


,

φ

i
,

υ
i












lim





inf


t
->









ϕ


(



u
1



(
t
)


,





,



u
_

N



(
t
)



)









subject





to


:












lim





inf


t
->



[




u
_

i



(
t
)


-





υ
i



Ω
i
s










φ
_


i
,

υ
i





(
t
)




]


0

,



i















lim





inf


t
->





[




φ
_


i
,

υ
i





(
t
)


-



u
_


i
,

υ
i



(

β
i

)




(
t
)



]



0

,



i




,


υ
i



Ω
i
s


,


β
i



i
s










0



φ

i
,

υ
i





(
t
)





u
i

ma





x



1


{



ω
1



(
t
)


=

υ
i


}



,



t


{

0
,
1
,
2
,






}



,

i



,






υ
i



Ω
i
s













lim





inf


t
->





[




u
_

i



(
t
)


-


u
i
c



(
t
)



]



0

,



i


𝒮
u












α


(
t
)






(
t
)



,



t


{

0
,
1
,
2
,






}








(
19
)







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∈custom character, and defines g(t)=ϕ(γ1(t), . . . , γN(t). Then, the Jensen's inequality indicates that

g(t)≤ϕ(γ1(t), . . . ,γN(t))  (20)


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)∈custom character(t), variables φi,υi(t) and an auxiliary vector γ(t) to solve:












Maximize







α


(
t
)


,


φ

i
,

υ
i





(
t
)


,

γ


(
t
)











lim





inf


t
->






g
_



(
t
)









subject





to


:






(
21.1
)









lim

t
->










γ
_

i



(
t
)


-



u
_

i



(
t
)






=
0

,



i








(
21.2
)







0



γ
i



(
t
)




u
i

ma





x



,



t


{

0
,
1
,
2
,






}



,

i







(
21.3
)










lim





inf


t
->



[




u
_

i



(
t
)


-





υ
i



Ω
i
s










φ
_


i
,

υ
i





(
t
)




]


0

,



i








(
21.4
)










lim





inf


t
->





[




φ
_


i
,

υ
i





(
t
)


-



u
_


i
,

υ
i



(

β
i

)




(
t
)



]



0

,



i




,


υ
i



Ω
i
s


,


β
i



i
s






(
21.5
)







0



φ

i
,

υ
i





(
t
)





u
i

ma





x



1


{



ω
1



(
t
)


=

υ
i


}



,



t


{

0
,
1
,
2
,






}



,

i



,






υ
i



Ω
i
s






(
21.6
)










lim





inf


t
->





[




u
_

i



(
t
)


-

u
i
c


]



0

,



i


𝒮
u







(
21.7
)








α


(
t
)






(
t
)



,



t


{

0
,
1
,
2
,






}







(
21.8
)







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∈custom character as the following form:











Q
i



(

t
+
1

)


=



Q
i



(
t
)


+





υ
i



Ω
i
s









φ

i
,

υ
i





(
t
)



-


u
i



(
t
)







(
22
)







The constraint (21.4) can be satisfied, if the following mean-rate stability condition is held:











lim

t







E


[


[


Q
i



(
t
)


]

+

]


t


=
0




(
23
)







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,υii)(t) of the first virtual queue is defined for all i∈custom character, υi∈Ωis, and βicustom character:

Di,υii)(t+1)=Di,υii)(t)+ui,υii)(t)−φi,υi(t)  (24)


The second virtual queue Zi(t) is defined for all i∈custom character:

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:











lim

t







E


[


[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+

]


t


=
0




(
27
)








lim

t







E


[


Z
i



(
t
)


]


t


=
0




(
28
)








lim

t







E


[


[


H
i



(
t
)


]

+

]


t


=
0




(
29
)







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,υii)(t)]+, Zi(t) and [Hi(t)]+ (divided by 2 for convenience) at time slot t, which is called a Lyapunov function as the total violation:










L


(
t
)




=






1
2






i










(


[


Q
i



(
t
)


]

+

)

2



+


1
2






i














υ
i



Ω
i
s


,


β
i











(


[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+

)

2




+


1
2






i











Z
i



(
t
)


2



+


1
2






i










(


[


H
i



(
t
)


]

+

)

2








(
30
)







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:











Δ


(
t
)


-

V






g


(
t
)






B
+




l











[


H
i



(
t
)


]

+



u
l
c



-

V






ϕ


(



γ
1



(
t
)


,





,


γ
N



(
t
)



)



+




i











Z
i



(
t
)





γ
i



(
t
)




+




i









{



[


Q
i



(
t
)


]

+







υ
i



Ω
i
s









φ

i
,

υ
i





(
t
)




}


-




i














υ
i



Ω
i
s


,


β
i












[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+




φ

i
,

υ
i





(
t
)





+




i














υ
i



Ω
i
s


,


β
i












[


D

l
,

υ
i



(

β
l

)




(
t
)


]

+




u

i
,

υ
i



(

β
i

)




(
t
)





-




i










{



[


Q
i



(
t
)


]

+

+


Z
i



(
t
)


+


[


H
i



(
t
)


]

+


}




u
i



(
t
)









(
31
)







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







B


=








i










(

u
i
max

)

2


+


1
2






i


S
u









(

u
i
max

)

2



+


1
2






i
















(

u
i
max

)

2






,





the first penalty upper bound term is









-
V







ϕ


(



γ
1



(
t
)


,





,


γ
N



(
t
)



)



+




i








Z
i



(
t
)





γ
i



(
t
)





,





the second penalty upper bound term is










i









{



[


Q
i



(
t
)


]

+







υ
i



Ω
i
s









φ

i
,

υ
i





(
t
)




}


-




l














υ
i



Ω
i
s


,


β
i












[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+




φ

i
,

υ
i





(
t
)










and the third penalty upper bound term is










i














υ
i



Ω
i
s


,


β
i












[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+




u

i
,

υ
i



(

β
i

)




(
t
)





-




i










{



[


Q
i



(
t
)


]

+

+


Z
i



(
t
)


+


[


H
i



(
t
)


]

+


}





u
i



(
t
)


.







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,υii)(t) of the first virtual value in the first virtual queue at the current time slot, the second virtual value Zi(t) in the second virtual queue at the current time slot, the third virtual value Hi(t) in the third virtual queue at the current time slot and the random event vector ω(t)∈Ω.


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,υi(t), and the forth term is a function of the individual utilities ui(t) and ui,υii)(t).


S1062: choosing the value of the third virtual variables γi(t) for all i∈custom character according to the second virtual value at the current time slot and the first penalty upper bound term to solve:










Minimize
-


V

ϕ

(



γ
1



(
t
)


,





,


γ
N



(
t
)



)

+




i











Z
i



(
t
)





γ
i



(
t
)













subject





to


:






0




γ
i



(
t
)




u
i
max


,



i









(
32
)







The closed-form solution of the problem (32) can take the following form for all i∈custom character:











γ
l



(
t
)


=

{





u
i
max

,






Z
i



(
t
)



0







min






{


V



Z
i



(
t
)



ln





2


,

u
i
max


}


,






Z
i



(
t
)


>
0









(
33
)







S1063: choosing the value of the second virtual variables φi,υi(t) for all i∈custom characterand υi ∈Ωis according to the random event vector, the first virtual value at the current time slot and the second penalty upper bound term to solve:











Minimize









i









{



[


Q
i



(
t
)


]

+







υ
i



Ω
i
s









φ

i
,

υ
i





(
t
)




}



-




l














υ
i



Ω
i
s


,


β
i
















[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+




φ

i
,

υ
i





(
t
)


















subject





to


:






0




φ

i
,

υ
i





(
t
)





u
i
max


1


{



ω
i



(
t
)


=

υ
i


}





i






,


υ
i



Ω
i
s







(
34
)







The closed-form solution of the problem (34) can take the following form for all i∈custom characterand υi ∈Ωis:











φ

i
,

υ
i





(
t
)


=

{





u
i
max

,







ω
i



(
t
)


=

υ
i


,



[


Q
i



(
t
)


]

+

<





β
i










[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+









0
,



otherwise








(
35
)







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:










Maximize
-




l














υ
i



Ω
i
s


,


β
i












[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+




u

i
,

υ
i



(

β
i

)




(
t
)





+




i










{



[


Q
i



(
t
)


]

+

+


Z
i



(
t
)


+


[


H
i



(
t
)


]

+


}




u
i



(
t
)
















subject





to


:







α


(
t
)







(
t
)







(
36
)







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,υii)(t) are computed according to formulas (14) and (18), respectively. The virtual queues Qi(t), Di,υii)(t), Zi(t) and Hi(t) are updated according to formulas (22), (24), (25) and (26), respectively.


The problem (36) is a non-linear integer programming problem, where ui(t) and ui,υii)(t) are complicate non-linear functions of α(t). The exhaustive algorithm for solving the problem (36) has a complexity of







O
(




i

















(
t
)





)

,





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∈custom character\{0} may be congested when the following condition is held: custom character1{α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:














i









1


{



α
i



(
t
)


=
j

}




R
i



(
t
)







x
b




C
j



(
t
)




,



j


𝒥


\



{
0
}








(
37
)







For each participating player 203-icustom character, 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-icustom character, and for each υi ∈Ωis and βicustom characteris, (βi≠0), this embodiment considers the definition (18) of ui,υii)(t). If υi≠ωi(t), ui,υii)(t)=0. If υii(t), ui,υii)(t)=ûi((βiī(t)), ω(t)). And this embodiment can estimate ui,υii)(t) according to the following two different cases, thereby forming a mapping table between the individual utilities and the transmission rate of the non-participating players.


1) If the player 203-i accesses to a network that is just the one suggested by the game manager 202, i.e., βii (t), then ui,υii)(t)=Ri(t).


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,βi (t)). Specifically, this embodiment defines the remaining capacity of the network j∈custom character\{0} at time slot t as Cjr(t)=xb Cj(t)custom character1{αi(t)=j}Ri(t), and assumes that Cjr(t)=Cjr(t−1). As a result, ui,υii)(t) can take the following formula:











u

i
,

υ
i



(

β
i

)




(
t
)


=



θ

i
,

β
i





(
t
)





R
i



(
t
)







(
39
)







where,









θ

i
,

β
i





(
t
)


=

f


(




[



x
b




C

β
i




(
t
)



-


C

β
i

r



(

t
-
1

)



]

+

+


R
i



(
t
)





C

β
i




(
t
)



)



,





Cβi(t) represents the capacity of the network corresponding to βi.


Therefore, ui,υii)(t) can be estimated by:











u

i
,

υ
i



(

β
i

)




(
t
)


=

{




0
,





υ
i




ω
i



(
t
)










R
i



(
t
)


,






υ
i

=


ω
i



(
t
)



,


β
i

=


α
i



(
t
)












θ

i
,

β
i





(
t
)





R
i



(
t
)



,






υ
i

=


ω
i



(
t
)



,


β
i




α
i



(
t
)












(
40
)







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∈custom character and j∈custom character, this embodiment defines the mapping relationship between the suggestion matrix and the action vector as:







a
ij

=

{





1
,






α
i



(
t
)


=
j






0
,



otherwise



,







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)∈custom character(t), the suggestion matrix constraint can be obtained:















j

𝒥












a
ij


=
1

,



i








(
41
)








a
ij



{

0
,
1

}


,



i








,

j

𝒥





(
42
)








a
ij

=
0

,



i




,

j








(
t
)







(
43
)







α(t) in (37) is substituted with the variables aij, and the second action vector constraint can be obtained:














i










a
ij




R
i



(
t
)







x
b




C
j



(
t
)




,



j


𝒥


\



{
0
}








(
44
)







α(t) in (37) is substituted with the variables aij, and the individual utilities of the participating players can be obtained:












u
i



(
t
)


=



R
i



(
t
)







j


𝒥


\



{
0
}









a
ij




,



i








(
45
)







Furthermore, this embodiment lets







θ

i
,
j


(

β
i

)


=

{




1
,




j
=

β
i









θ

i
,

β
i





(
t
)


,




j


β
i











for all i∈custom character, j∈custom character and βicustom character. α(t) in (40) is substituted with the variables aij, and the individual utilities of the non-participating players can be obtained:

ui,υii)(t)custom characterûi((βiī(t)),ω(t))1{ωi(t)=υi}, ∀i∈custom characteri∈Ωisicustom characteris  (46)


According to (41)-(46), this embodiment can transform the problem (36) into the following integer linear programming problem:










Maximize

{

a
ij

}











i












j

𝒥








c
ij



a
ij








(
47
)






subject





to


:

















i











R
l



(
t
)




a
ij






x
b




C
j



(
t
)




,



j


𝒥


\



{
0
}




















j

𝒥







a
ij


=
1

,



i

















a
ij



{

0
,
1

}


,



i




,

j

𝒥














a
ij

=
0

,



i




,

j




(
t
)















where, the weights cij are defined as:










c
ij

=

{






(



E
l



(
t
)


-


F
ij



(
t
)



)




R
i



(
t
)



,






j


𝒥


\



{
0
}











-


F
ij



(
t
)






R
i



(
t
)



,




j
=
0









(
48
)





where












E
i



(
t
)


=



[


Q
i



(
t
)


]

+

+


Z
i



(
t
)


+


[


H
i



(
t
)


]

+












and












F
ij



(
t
)


=






υ
i



Ω
i
s


,


β
i



A
i
s











[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+



θ
ij

(

β
i

)



1



{



ω
i



(
t
)


=

υ
i


}

.

















i












j

𝒥








c
ij



a
ij















is the fourth penalty upper bound term,











i











R
i



(
t
)




a
ij






x
b




C
j



(
t
)




,



j


𝒥


\



{
0
}









constitutes the second action vector constraint, and






{











j

𝒥








a
ij


=
1

,



i












a
ij



{

0
,
1

}


,



i








,

j

𝒥









a
ij

=
0

,



i




,

j








(
t
)
















constitutes





the





suggestion





matrix






constraint
.






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:










c
ij

=

{






R
i



(
0
)


,






j


𝒥


\



{
0
}









0
,




j
=
0









(
49
)







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 FIG. 4.



FIG. 4 is a schematic structural diagram of a network selection apparatus for the integrated cellular and drone-cell networks according to an exemplary embodiment of the present disclosure. As shown in FIG. 4, a network selection apparatus 300 includes a transceiver 311 and a memory 312. The transceiver 311 is used for collecting the drone-cell network capacity information Cj(t), j∈custom characterd, the cellular network capacity information C1(t), the users' accessible network set information custom character(t), i∈custom character, and data transmission rate information Ri(t), i∈custom character, to form the random event vector ω(t). Meanwhile, the transceiver 311 is also responsible for transmitting the suggested action information formulated by the game manager 202 to each user device 330. The memory 312 can be any type of computer readable medium, and it is responsible for storing information such as parameters, state data, action data, and virtual queue data. The network selection apparatus 300 further includes a processor 313. The processor 313 can be any type of central processing unit, and it is responsible for processing data in the EFNS method. Specifically, the processor 313 obtains the action vector according to the drone-cell network capacity information, the cellular network capacity information, the users' accessible network set information, and data transmission rate information and the fourth selection model. Each user device 330 includes a network access system 331 for controlling the user device 330 to select a network to access according to the received action. The fourth selection model has been described in detail in the foregoing embodiment, and details are not described herein again.


Referring further to FIG. 4, the network selection apparatus provided by this embodiment also includes a human-computer interaction module 314, which includes a display and an operator input interface. The display can display the results to a computer operator 340. The operator input interface can obtain information inputted by the computer operator 340 from one or more input devices such as a keyboard and a mouse.


Referring further to FIG. 4, the network selection apparatus provided by this embodiment further includes networks 320. FIG. 5 is a flowchart of the method performed by a processor in the network selection apparatus for the integrated cellular and drone-cell networks according to the embodiment shown in FIG. 4 of the present disclosure. As shown in FIG. 5, the processor 313 performs the following actions:


The state information ω(t) and the virtual queue information Qi(t), Di,υii)(t), Zi(t), and Hi (t) at the current time slot t are received. And the output information is the suggested action vector α(t) of the game manager 202 at the current time slot t and the virtual queue information Qi(t+1), Di,υii)(t+1), Zi(t+1), and Hi(t+1) of the next time slot t+1.


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,υii)(0)=0, Zi(0)=0 and Hi(0)=0 are initialized, and these virtual queues are stored in the memory 312.


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∈custom character, 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∈custom character and υi ∈Ωis, the processor 313 calculates the second auxiliary variables φi,υi(t) by the formula (35). The auxiliary variables φi,υi(t) will be temporarily stored in the memory 312 until the end of Step 6.


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,υii)(t) using formulas (4) and (18), respectively, calculates Qi(t+1) Di,υii)(t+1), Zi(t+L) and Hi(t+1) using formulas (22), (24), (25) and (26), and updates the virtual queues in the memory 312.


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:









TABLE 1





One-step transition probability matrix P of μR (t)



















0.8
0.2
0
0
0


0.2
0.6
0.2
0
0


0
0.2
0.6
0.2
0


0
0
0.2
0.6
0.2


0
0
0
0.8
2









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∈custom character, 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









s
Q



(
t
)


=



max

i







[


Q
i



(
t
)


]

+


t


,



s
D



(
t
)


=



max


i



,


υ
i



Ω
i
s


,


β
i








[


D

i
,

υ
i



(

β
i

)




(
t
)


]

+


t


,







s
Z



(
t
)


=



max

i









Z
i



(
t
)





t


,


and







s
H



(
t
)



=



max

i


s
u






[


H
i



(
t
)


]

+


t







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,







i
.
e
.

,




t
=
0


T
-
1








i
















u
i



(
t
)


.







Fairness: this embodiment uses the Jain's fairness index, defined as









(




l










u
_

i


)

2



/



(

N





i










u
_

l
2



)


,





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.,








u
_

i

=


1
T






t
=
0


T
-
1






u
i



(
t
)


.







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.



FIG. 6 illustrates the change of stability variables with time when the network selection method according to the present disclosure is used and when the number of users is N=50 and the number of drones is Md=6. As shown in FIG. 6:


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.



FIG. 7 illustrates the effect of the number of users IN on the running time of the network selection method proposed by the present disclosure when the number of drones is Md=6. FIG. 8 illustrates the effect of the number of drones Md on the running time of the network selection method proposed by the present disclosure when the number of users is N=50. As shown in FIG. 7 and FIG. 8:


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.



FIG. 9 illustrates the effect of the number of users N on the total user utilities obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of drones is Md=6. FIG. 10 illustrates the effect of the number of drones Md on the total user utilities obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of users is N=50. As shown in FIG. 9 and FIG. 10:


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.



FIG. 11 illustrates the effect of the number of users N on the Jain's fairness indices obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of drones is Md=6. FIG. 12 illustrates the effect of the number of drones Md on the Jain's fairness index obtained by the network selection method proposed by the present disclosure and the comparison methods when the number of users is N=50. As shown in FIG. 11 and FIG. 12:


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.

Claims
  • 1. A network selection method for integrated cellular and drone-cell networks, comprising: acquiring a dynamic network model and a dynamic user model, wherein the dynamic network model comprises 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 comprises 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, wherein the accessible network sets of the users comprise the drone-cell networks and/or the cellular network;generating an action vector according to the random event vector, wherein the action vector 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, wherein the first selection model comprises a first objective function and a first constraint, wherein the first objective function is a proportional fairness function with the time average of individual utilities of the users as independent variables, and the first constraint comprises at least a first coarse correlated equilibrium constraint, a first minimum individual time average utility constraint and a first action probability constraint, wherein 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; wherein the time average of the individual utilities is obtained according to the individual utilities, the random event probability and the action probability under the condition of the random event vector, wherein 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 vector occurs; andobtaining 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.
  • 2. The method according to claim 1, wherein after the obtaining an individual utility of each user according to the action vector and the random event vector, the method further comprises: constructing a second selection model, wherein the second selection model comprises a second objective function and a second constraint, wherein the second objective function is a proportional fairness function with time average expectation of the individual utilities as independent variables, and the second constraint comprises at least second coarse correlated equilibrium constraints, a second minimum individual time average utility constraint and a second auxiliary variable constraint, wherein the second coarse correlated equilibrium constraints are used for constraining the 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, and the second auxiliary variable constraint is used for constraining the second auxiliary variables; andobtaining a value of the action vector according the second selection model, and determining networks that the users choose to access according to the value of the action vector.
  • 3. The method according to claim 2, wherein after the obtaining an individual utility of each user according to the action vector and the random event vector, the method further comprises: constructing a third selection model according to the second selection model, wherein the third selection model comprises a third objective function and a third constraint, wherein the third objective function is time average expectation of a proportional fair function with third auxiliary variables as independent variables, and the third constraint comprises at least the second coarse correlated equilibrium constraints, the second minimum individual time average utility constraint, the second auxiliary variable constraint, and third auxiliary variable constraints, wherein the second coarse correlated equilibrium constraints are used for constraining 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 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 time average expectation of the third auxiliary variables and the time average expectation of the individual utilities;obtaining a value of the action vector according the third selection model, and determining networks that the users choose to access according to the value of the action vector.
  • 4. The method according to claim 3, wherein after the constructing a third selection model, the method further comprises: constructing a fourth selection model;obtaining a value of the action vector according the fourth selection model, and determining networks that the users choose to access according to the value of the action vector;wherein the fourth selection model is that a difference value between a drift of total violation and an utility is less than or equal to a penalty upper bound;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 a first virtual queue at the current time slot is generated according to 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 a 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 a third virtual queue at the current time slot is generated according to the 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, wherein the first virtual value at 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;the utility comprises: the proportional fair function and a penalty coefficient, with the third auxiliary variables as independent variables; andthe penalty upper bound comprises: a constant term, a first penalty upper bound term, a second penalty upper bound term and a third penalty upper bound term, wherein the first penalty upper bound term comprises the third auxiliary variables and the second virtual value, the second penalty upper bound term comprises the second auxiliary variables and the first virtual value, and the third penalty upper bound comprises the individual utilities, the first virtual value, the second virtual value and the third virtual value.
  • 5. The method according to claim 4, wherein the obtaining a value of the action vector according the fourth selection model, and determining networks that the users choose to access according to the value of the action vector, specifically comprises: obtaining 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 random event vector at the current time slot;obtaining values of the third auxiliary variables according to the second virtual value at the current time slot and the first penalty upper bound term;obtaining values of the second auxiliary variables according to the random event vector at the current time slot, the first virtual value at the current time slot and the second penalty upper bound term;obtaining the value of the action vector according to the random event vector at the current time slot, 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, and determining the networks that the users choose to access according to the value of the action vector;the obtaining the value of the action vector according to the random event vector at the current time slot, 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, specifically comprises:constructing a suggestion matrix constraint and a first action vector constraint, wherein a suggestion matrix represents that the users are suggested to access to the drone-cell networks and/or the cellular network, and the suggestion matrix constraint is used for constraining the suggestion matrix;using a mapping table between the individual utilities and transmission rate and a mapping relationship between the suggestion matrix and the action vector to process the third penalty upper bound term, to generate a fourth penalty upper bound term;using the mapping relationship between the suggestion matrix and the action vector to process the first action vector constraint, to obtain a second action vector constraint; andobtaining a value of the suggestion matrix according the fourth penalty upper bound term, the suggestion matrix constraint and the second action vector constraint, and determining the networks that the users choose to access according to the value of the suggestion matrix.
  • 6. The method according to claim 1, wherein the first objective function specifically comprises: ƒ1=ϕ(ū1, . . . , ui, . . . , uN), wherein, ƒ1 represents the first objective function, ū1, represents the time average of the individual utility of the i -th user, and ϕ(ū1, . . . , ui, . . . , uN)represents a proportional fair function; the first coarse correlated equilibrium constraint specifically comprises:
  • 7. The method according to claim 3, wherein the second objective function specifically comprises:
  • 8. The method according to claim 4, wherein the fourth selection model specifically comprises:
  • 9. The method according to claim 5, wherein the fourth penalty upper bound term specifically comprises:
Priority Claims (1)
Number Date Country Kind
201811353219.4 Nov 2018 CN national
US Referenced Citations (14)
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
Foreign Referenced Citations (4)
Number Date Country
103332290 Oct 2013 CN
106054928 Oct 2016 CN
106604205 Apr 2017 CN
107517158 Dec 2017 CN
Non-Patent Literature Citations (4)
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.
Related Publications (1)
Number Date Country
20200154349 A1 May 2020 US