CN103428747A - Aviation self-organizing network wireless link stability prediction method - Google Patents

Aviation self-organizing network wireless link stability prediction method Download PDF

Info

Publication number
CN103428747A
CN103428747A CN2013103978882A CN201310397888A CN103428747A CN 103428747 A CN103428747 A CN 103428747A CN 2013103978882 A CN2013103978882 A CN 2013103978882A CN 201310397888 A CN201310397888 A CN 201310397888A CN 103428747 A CN103428747 A CN 103428747A
Authority
CN
China
Prior art keywords
prime
node
stage
max
aircraft
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013103978882A
Other languages
Chinese (zh)
Other versions
CN103428747B (en
Inventor
雷磊
王丹
罗诚
王成华
张雅静
朱马君
朱明�
李晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanxi Aerospace Industry Research Institute Group Co ltd
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201310397888.2A priority Critical patent/CN103428747B/en
Publication of CN103428747A publication Critical patent/CN103428747A/en
Application granted granted Critical
Publication of CN103428747B publication Critical patent/CN103428747B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

The invention discloses an aviation self-organizing network wireless link stability prediction method. The method includes the steps of firstly, establishing an aviation self-organizing network node mobile model; secondly, determining a node movement rate probability density function in each state and a probability density function of a relative speed between every two nodes according to the movement features of the nodes in the mobile model in each state; thirdly, determining a distribution function of the duration time of a wireless link through combination of the distance between every two nodes and node relative movement direction information on the basis; finally, obtaining stability factors of the multi-hop wireless link according to the distribution function of the duration time of the wireless link, wherein the mobile model is composed of seven states, namely, the acceleration taking off state, the constant-speed rising state, the smooth and steady flying state, the steering state, the constant-speed descending speed, the deceleration stopping state and the static state. The stability factors of the link can be used as the importation foundation for route selection through the airplane nodes, and therefore the wireless transmission link with the high stability is set up for an aviation self-organizing network. The validity of the method is proved through a simulation experiment in the QualNet simulation environment.

Description

A kind of aviation self-organizing network wireless link stability prediction method
Technical field
The invention belongs to field of wireless, particularly aviation self-organizing network wireless link stability prediction method.
Background technology
Wireless self-organization network (Wireless Ad hoc Networks) is a kind of novel wireless communication network that does not rely on static infrastructure, and it is comprised of one group of mobile node with wireless transmitter, and networking fast and flexible, reliability are high.Each node in network is terminal and router, can forward the packet of other node in automatic network.Can be under the procotol of layering be controlled after node start with the mode quickly networking of multi-hop ad hoc, when part of nodes breaks down or destroyed and after quitting work, the operation of whole network can not be affected, thereby network has very strong survivability and self-healing ability.
In recent years, wireless self-organization network has become a kind of important wireless networking mode, in the civilian and military field, all is widely used.Aviation self-organizing network (Aeronautical Ad hoc Networks, AANET) is the product that radio self organizing network technology is applied in the multi-platform networking of aviation.At civil area, the replacement of aviation self-organizing network is expensive, the broadband connections satellite of long time delay, becomes the new selection of following wideband multimedia air communications.The researcher proposes, and between the passenger plane of flight, builds self-organizing network, and the mode that can forward with multi-hop be set up the communication link between aircraft and ground base station, for the passenger provides the aerial broadband the Internet access service of cheap, low time delay on the road.
Aircraft node high-speed mobile (cruising speed is generally between 700km/h to 1000km/h) in the aviation self-organizing network, network topology changes fast, thereby the internodal wireless transmission link of aircraft frequently ruptures.When any hop link in the multi-hop transmission path ruptures, the aircraft node must start route finding process, finds other available transmission paths.A large amount of route finding process has consumed network bandwidth resources greatly, causes that network congestion and data-message transmission time delay increase, and finally causes the degradation of network in general performance.Therefore, in the design of aviation self-organizing network Routing Protocol, must consider the duration of wireless link between node, thereby find by route discovery the wireless transmission link that stability is high, alleviate the impact of network topology change on the network in general performance.
As the above analysis, between the aircraft node, the prediction of wireless link stability is of great significance the overall performance tool that promotes the aviation self-organizing network.Existing wireless link stability prediction method roughly can be divided three classes: the link stability Forecasting Methodology based on received signal power, based on node location, move the link stability Forecasting Methodology of (geometric properties) and the link stability Forecasting Methodology of Based on Probability estimation theory.The basic principle of the link stability Forecasting Methodology based on received signal power is: when node sends packet, radio signal power is the constantly decay with the increase of transmission range in transmitting procedure.After neighbor node receives this signal, according to the distance between the estimation of received signal power and path loss model and sending node, and according to this distance upgrade and sending node between link-state information.The method simple, intuitive, but in actual environment, the propagation of wireless signal may be subject to the impact of the many factors such as atmosphere, sexual intercourse, complex-terrain, deviation inevitably appears in the distance between the sending node obtained by said method estimation and receiving node, and then has influence on final routing strategy.
The basic principle of the link stability prediction of moving based on node location is: according to movable informations such as existing node location, speed, directions, derive and obtain the information of link-state change.The link stability prediction of moving based on node location mainly is divided into again two kinds.The first is based on link Duration Prediction (Prediction based, PBR), and this agreement is as follows to the Forecasting Methodology of link existent time: given two mobile node i and j, define between them the distance be | d Ij|, corresponding translational speed is respectively v iAnd v j, communication distance is R, and definition s is used for distinguishing motion or counter motion in the same way, two internodal vertical width are w, derive thus and obtain the predicted value of link duration between two mobile node i and j.What this Forecasting Methodology was only considered is the situation that movement direction of nodes is identical or contrary, can't solve the problem of two movement direction of nodes not on the same straight line time.Therefore and what this algorithm was only considered is the linear uniform motion of node, can't solve the problems such as acceleration, deceleration, turning of node in the aviation self-organizing network.The second is based on the prediction (Multipath Doppler Routing Algorithm, MUDOR) of many Pood frequency displacement, and Doppler frequency shift is the key factor that the MUDOR agreement is considered.Pass between relative velocity and Doppler frequency shift is v=c (flf 0-1), wherein c is the light velocity, and f is expected frequency, f 0It is observed frequency.When the expectation frequency is less than observing frequency, for be defined as-c of doppler values (Doppler Value, the DV) (flf that reflects link load 0-1), otherwise be defined as+2c (flf 0-1).Doppler values is less, and link stability is better.Whether the MUDOR algorithm neither needs the support of GPS equipment also reliable without the attenuation model of considering the reception signal.But, due to the impact of only considering relative speed, and ignored link establishment internodal initial distance constantly fully.Thereby it can only carry out qualitative forecasting to link property, and can't obtain the predicted value of link duration accurately.
The link stability Forecasting Methodology of Based on Probability estimation theory is: for any two mobile nodes 1 and 2, the communication distance of node 1 is divided into to the n section, every section represents respectively a state S i, writ state S N+1Mean absorbing state, the zone beyond representation node 1 communication range, p IjMean that node 2 is by state S iMove to state S after the unit interval jProbability, define a step transition probability matrix P=|p Ij|, i, j ∈ [1, n+1].
Figure BSA0000094734960000031
Mean that communication initial time node 2 is in state S iProbability, the initial distance distribution matrix of defined node 2
Figure BSA0000094734960000032
I ∈ [1, n+1], can draw linkage availability and the probability density function of link duration thus, analyzes accordingly link stability.Because a step transition probability matrix in the method is changeless, therefore, the link stability that the method analysis obtains is only relevant with the distribution of node initial distance, and irrelevant with the initial relative speed of node.
Summary of the invention
The objective of the invention is to propose a kind of wireless link stability prediction method that is applicable to the aviation self-organizing network, for Design of Routing Protocol provides important evidence, alleviate the impact of network topology change on the network in general performance.In order to realize this purpose, step of the present invention is:
Step 1: build aviation self-organized network nodes mobility model, mobility model is comprised of seven states: accelerate to take off, constant speed rising, smooth flight, turning, constant speed descend, slowing down stops with static;
Step 2: according to node in above-mentioned mobility model, at the motion feature of each state, determine the node movement rate probability density function in each state;
Step 3: the speed probability density function according to node in each motion state, determine the probability density function of relative speed between node;
Step 4: according to the probability density function of relative speed between node, the distance between node, and node direction of relative movement, determine the distribution function of link duration;
Step 5: according to link duration distribution function, calculate the multi-hop link stability factor.
The wireless link stability prediction method of the aviation self-organizing network that the present invention proposes is realized in the QualNet4.5 network simulation environment.Simulating area is 10000 * 10000 * 1000m 3, network node adds up to 100, in the simulating area random distribution.Physical layer adopts the DSSS model, and simulation time is 300s, and all the other simulation parameters are as shown in table 1.Relevant parameter value in the node motion model is as shown in table 2
Table 1 simulation parameter
Figure BSA0000094734960000041
Relevant parameter value in table 2 node motion model
Figure BSA0000094734960000042
Accompanying drawing 4 has provided the theoretical value of the link duration distribution function that the present invention calculates and the contrast of simulation value.The consistency of simulation value and theoretical value has illustrated that the present invention determines the validity of link duration distribution function method between node.Accompanying drawing 5 and accompanying drawing 6 have provided the network performance simulation result contrast before and after the link stability Forecasting Methodology of introducing the present invention's proposition in AODV (ad hoc on-demand distance vector) Routing Protocol.Simulation results show, the link stability Forecasting Methodology of introducing the present invention's proposition in Routing Protocol can effectively provide the saturation throughput performance of aviation self-organizing network, and reduces Packet Delay simultaneously.
The accompanying drawing explanation
Fig. 1 is aviation self-organized network nodes mobility model state transition graph;
Fig. 2 is node relative velocity schematic diagram;
Fig. 3 is that node link duration distribution function calculates schematic diagram;
Fig. 4 is theoretical value and the simulation value comparison diagram distributed link duration of calculating of the present invention;
Fig. 5 is the network throughput performance simulation comparison diagram of introducing the link stability Forecasting Methodology front and back of the present invention's proposition in the AODV Routing Protocol;
Fig. 6 is the network delay performance simulation comparison diagram of introducing the link stability Forecasting Methodology front and back of the present invention's proposition in the AODV Routing Protocol.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
The aviation self-organizing network link stability Forecasting Methodology that the present invention proposes realizes in Wireless Network Simulation environment QualNet4.5.
Below provide the specific embodiment of the invention step:
Step 1: build aviation self-organized network nodes mobility model.
As shown in Figure 1, to be divided into be seven states to this mobility model: accelerate to take off, constant speed rising, smooth flight, turning, constant speed descend, slowing down stops with static.The aircraft node, from accelerating takeoff phase, after speed accelerates to set-point, enters the constant speed ascent stage.After rising to the regulation flying height, aircraft will start cruising flight (smooth flight)., if meet with special circumstances, as prominent, meet thunder and lightning weather awing, aircraft need to be revised course line, and now, aircraft will enter the turning stage.Turn after end, aircraft is got back to regular shipping lines, i.e. the smooth flight state.When being about to arrive destination, aircraft enters the constant speed decline stage.When dropping to ground proximity, aircraft enters the deceleration stop phase, slows down gradually until finally stop, and enters quiescent phase.
The design parameter of each motion state is described as follows:
1, accelerate takeoff phase (α stage)
Targeted rate v α∈ U[v α, min, v α, max], continue timeslot number α ∈ U[α min, α max], the horizontal movement direction
Figure BSA0000094734960000061
Vertical motion direction θ α=pi/2.Wherein U means to be uniformly distributed.So acceleration a in this stage αFor
a α = v α - v t 0 t α - t 0 = v α αΔt - - - ( 1 )
Δ t means time step, and the whole period of motion take Δ t and divided as time slot.
2, constant speed ascent stage (β stage)
After accelerating to finish takeoff phase, aircraft starts to climb at a certain angle, speed v β=v α, continue timeslot number β ∈ U[β min, β max], the horizontal movement direction
Figure BSA0000094734960000063
Vertical motion direction θ β∈ U[0, pi/2].
3, smooth flight (γ stage)
Cruising phase in this stage simulated aircraft motion, aircraft keeps flying speed constant, along the preferential direction level, flies at a constant speed.V γ=v β, continue timeslot number γ ∈ U[γ min, γ max], θ γ=pi/2.
4, turn (ζ stage)
In flight course, may need to change original course line because of certain reason, at this moment aircraft will enter the turning stage, and its turning can be approximately uniform circular motion.Before aircraft is turned, according to random selected horizontal movement direction
Figure BSA0000094734960000071
With vertical motion direction θ ζ, make uniform circular motion, in this stage, the speed v of aircraft ζ=v γ.The turning stage continues timeslot number ζ ∈ U[ζ min, ζ max].After the γ stage finishes, aircraft may enter the turning stage, also may enter the constant speed decline η stage because being about to arrive destination.After supposing that each smooth flight finishes, the probability that aircraft enters the turning stage is P γ-ζ, the probability that enters the constant speed decline stage is P γ-η, obvious P γ-ζ+ P γ-η=1.Suppose that aircraft this process before fall-retarding will carry out K time, K=0,1,2 ...In this case, the γ stage carries out K+1 time, and the ζ stage carries out K time.The desired value of K is as follows
E [ K ] = Σ K = 0 ∞ K ( P γ - ξ ) K ( 1 - P γ - ξ ) = P γ - ξ 1 - P γ - ξ - - - ( 2 )
5, constant speed descends (η stage)
When being about to arrive destination, enter this stage.In this stage, aircraft is with at the uniform velocity straight line decline of fixed speed, and initial velocity magnitude and direction are respectively: v η=v γ,
Figure BSA0000094734960000073
θ η∈ U[pi/2, π].Continue timeslot number η ∈ U[η min, η max].
6, slow down and stop (d stage)
After arriving ground, aircraft will be decelerated to and stop, and can regard the motion in this stage as uniformly retarded motion.Initial velocity magnitude and direction are respectively: v d=v η, θ d=pi/2.Duration d ∈ U[d min, d max], can obtain accordingly acceleration and be
a d = 0 - v η dΔt = - v η dΔt - - - ( 3 )
7, static (p stage)
The duration p ∈ U[p that aircraft remains static min, p max], then start flight next time.
Step 2: the probability density function of determining each stage speed of node.
For constant speed rising, smooth flight, turning and constant speed this four-stage that descends, because the node movement rate all remains unchanged, so the node movement rate distributes to meet and is uniformly distributed, as shown in the formula
f ( v ) = 1 v max - v min v ∈ ( v min , v max ) 0 else - - - ( 4 )
For accelerating takeoff phase, node is with initial velocity v 0=0, accelerate to given target velocity v α, v αEvenly choose from [v min, v max].Simultaneously, due to v α∈ [v min, v max], therefore no matter v αValue is how many, v ∈ [0, v min] value will occur with probability 1 takeoff phase accelerating, due to even acceleration, so their probability density is also identical, be made as k.And v ∈ [v min, v max] can occur in this stage will be by v αValue determine.V αValue is in v minWith v minProbability between+Δ v is Δ v/ (v max-v min), to accelerate in take-off process, speed v comprises v minThe probability of+Δ v is 1-Δ v/ (v max-v min).According to above analysis, can obtain accelerating takeoff phase, the probability density function of node movement rate is as follows
f &alpha; ( v ) = k , 0 &le; v &le; v min k ( 1 - v - v min v max - v min ) v min < v &le; v max - - - ( 5 )
Due to
&Integral; - &infin; + &infin; f &alpha; ( v ) dv = 1 - - - ( 6 )
Can obtain
k = 2 ( v max + v min ) - - - ( 7 )
Thereby formula (5) can be expressed as
f &alpha; ( v ) = 2 v max + v min , 0 &le; v &le; v min 2 v max + v min ( 1 - v - v min v max - v min ) , v min < v &le; v max - - - ( 8 )
Equally, for the deceleration stop phase, also can obtain the conclusion identical with formula (8) according to being similar to top analysis.
Step 3: the probability density function of determining relative speed between the aircraft node.
The speed of supposing optional two mobile nodes is respectively v 1, v 2, as shown in Figure 2, the relative velocity v ' of these two mobile nodes=v 1-v 2, angle ω ∈ [0, π], establish v 1, v 2, v ' mould be respectively v 1, v 2, v ', have
v &prime; = v 1 2 + v 2 2 - 2 v 1 v 2 cos w w = arccos v 1 2 + v 2 2 - v &prime; 2 2 v 1 v 2 - - - ( 9 )
Due to v 1, v 2, ω is separate, so, v 1, v 2, ω joint probability density function can be written as
f w , v 1 , v 2 ( w , v 1 , v 2 ) = 1 &pi; f v 1 ( v 1 ) f v 2 ( v 2 ) - - - ( 10 )
Convert by Jacobi, can obtain v 1, v 2, v ' joint probability density function
f v &prime; , v 1 , v 2 ( v &prime; , v 1 , v 2 ) = f w , v 1 , v 2 ( w , v 1 , v 2 ) | &PartialD; w &PartialD; v &prime; |
= 2 v &prime; f w , v 1 , v 2 ( w , v 1 , v 2 ) 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2 - - - ( 11 )
= 2 v &prime; f v 1 ( v 1 ) f v 2 ( v 2 ) &pi; 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2
The probability density function that can be drawn relative speed v ' by formula (11) is
f v &prime; ( v &prime; ) = &Integral; v 2 , min v 2 , max &Integral; v 1 , min v 1 , max 2 v &prime; f v 1 ( v 1 ) f v 2 ( v 2 ) &pi; 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2 d v 1 d v 2 - - - ( 12 )
Definition status S set={ α, beta, gamma, ζ, η, d, p}.Then, the rate distribution of each state obtained according to formula (4) and (8), in conjunction with formula (12), the probability density function that can obtain the relative speed under any two states is
f v &prime; E , G ( v &prime; ) = &Integral; v G min v G max &Integral; v E min v E max 2 v &prime; f ( v E ) f ( v G ) &pi; 4 v E 2 v G 2 - ( v E 2 + v G 2 - v &prime; 2 ) 2 d v E d v G - - - ( 13 )
Wherein, E, G all belongs to S.
Definition p (E), p (G) is respectively node constantly in E, the probability in G stage, the relative speed under aviation self-organizing network nodes mobility model is distributed as
f v &prime; ( v &prime; ) = &Sigma; E &Element; S &Sigma; G &Element; S p ( E ) p ( G ) f v &prime; E , G ( v &prime; ) - - - ( 14 )
Step 4: the distribution function that calculates the link duration.
As shown in Figure 3, consider the wireless link between node 1 and node 2, using node 2 as with reference to node, node 1 and node 2 apart from being d 0, direction and the d of two node speed of related movement v ' 0Between angle theta obey and be uniformly distributed at interval [0, π], R is the node transmission range; Use geometric knowledge, obtain the distribution function of link duration
F ( t ) = P ( T &le; t ) = P ( d 1 v &prime; &le; t ) = P ( d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t )
( 15 )
= &Integral; &Integral; &Integral; d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t f ( &theta; , d 0 , v &prime; ) d&theta; dd 0 dv &prime;
θ, d 0, tri-variablees of v ' are separate, therefore formula (15) can be expressed as
F ( t ) = &Integral; &Integral; &Integral; d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t f ( &theta; ) f ( d 0 ) f v &prime; ( v &prime; ) d&theta; dd 0 dv &prime; - - - ( 16 )
The probability density function of θ is
f ( &theta; ) = 1 &pi; - - - ( 17 )
By
F ( d 0 ) = P ( d &le; d 0 ) = 4 3 &pi;d 0 3 4 3 &pi;R 3 = d 0 3 R 3 - - - ( 18 )
Can obtain
f ( d 0 ) = 3 d 0 2 R 3 - - - ( 19 )
By qualifications
d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t - - - ( 20 )
Can obtain
cos &theta; &le; v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; td 0 - - - ( 21 )
Order
k=∫f(θ)dθ (22)
{。##.##1},
k = 0 v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; td 0 &le; - 1 &pi; - arccos v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; td 0 &pi; - 1 < v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; td 0 < 1 1 v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; td 0 &GreaterEqual; 1 - - - ( 23 )
Therefore, formula (16) can be expressed as
F ( t ) = &Integral; 0 2 v &prime; max &Integral; 0 R kf ( d 0 ) f v &prime; ( v &prime; ) dd 0 dv &prime; - - - ( 24 )
By formula (14), (17), (19), (23) substitution formula (24), can try to achieve the expression formula of link duration distribution function F (t).
Step 5: calculate the link stability factor.
For the stability to wireless link between node, predicted, the present invention has introduced the link stability factor and computational methods thereof.As shown in Figure 3, the present invention adjusts the distance according to the initial phase between node 1 and node 2 and initial velocity is derived and to be drawn link duration T between two nodes 12
T 12 = d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; - - - ( 25 )
Then by T 12Be to compare 0.9 o'clock corresponding link duration with link duration distribution function F (t) value, obtain
Figure BSA0000094734960000116
T ^ 12 = T 12 T F ( t ) = 0.9 - - - ( 26 )
The present invention is defined as the wireless link stability factor of 2 of node 1 and nodes
S 12 = min ( T ^ 12 , 1 ) - - - ( 27 )
For the multi-hop wireless link be comprised of (N-1) bar one hop link, its wireless link stability factor is defined as the product of many one hop link stability factors,
S 1N=S 12S 23S 34…S (N-1)N (28)
The wireless link stability factor can be used for predicting the stability state of current wireless link.In actual applications, determine F (t) according to the actual value of each motion state parameters in the aircraft mobility model, then calculated the stability factor of wireless link between node by formula (28), and the important evidence using it as Route Selection, thereby set up the wireless transmission link that stability is high for the aviation self-organized network nodes.
The content be not described in detail in the present patent application book belongs to the known prior art of professional and technical personnel in the field.

Claims (6)

1. an aviation self-organizing network wireless link stability prediction method, the step adopted is:
Step 1: build aviation self-organized network nodes mobility model, mobility model is comprised of seven states: accelerate to take off, constant speed rising, smooth flight, turning, constant speed descend, slowing down stops with static;
Step 2: according to node in above-mentioned mobility model, at the motion feature of each state, determine the node movement rate probability density function in each state;
Step 3: the speed probability density function according to node in each motion state, determine the probability density function of relative speed between node;
Step 4: according to the probability density function of relative speed between node, the distance between node, and node direction of relative movement, determine the distribution function of link duration;
Step 5: according to link duration distribution function, calculate the multi-hop link stability factor.
2. a kind of aviation self-organizing network wireless link stability prediction method according to claim 1 is characterized in that the concrete construction method of aviation self-organized network nodes mobility model is:
The aircraft node, from accelerating takeoff phase, after speed accelerates to set-point, enters the constant speed ascent stage, and after rising to the regulation flying height, aircraft will start cruising flight (smooth flight); , if meet with special circumstances, as prominent, meet thunder and lightning weather awing, aircraft need to be revised course line, and now, aircraft will enter the turning stage; Turn after end, aircraft is got back to regular shipping lines, i.e. the smooth flight state; When being about to arrive destination, aircraft enters the constant speed decline stage, and when dropping to ground proximity, aircraft enters the deceleration stop phase, slows down gradually until finally stop, and enters quiescent phase;
The design parameter of each motion state is described as follows:
(1) accelerate takeoff phase (α stage)
Targeted rate v α∈ U[v α, min, v α, max], continue timeslot number α ∈ U[α min, α max], the horizontal movement direction
Figure FSA0000094734950000011
Vertical motion direction θ α=pi/2; Wherein U means to be uniformly distributed; So acceleration a in this stage αFor:
a &alpha; = v &alpha; - v t 0 t &alpha; - t 0 = v &alpha; &alpha;&Delta;t - - - ( 1 )
Δ t means time step, and the whole period of motion take Δ t and divided as time slot;
(2) constant speed ascent stage (β stage)
After accelerating to finish takeoff phase, aircraft starts to climb at a certain angle, speed v β=v α, continue timeslot number β ∈ U[β min, β max], the horizontal movement direction
Figure FSA0000094734950000024
Vertical motion direction θ β∈ U[O, pi/2];
(3) smooth flight (γ stage)
Cruising phase in this stage simulated aircraft motion, aircraft keeps flying speed constant, along the preferential direction level, flies at a constant speed; v γ=v β, continue timeslot number γ ∈ U[γ min, γ max], θ γ=pi/2;
(4) turn (ζ stage)
In flight course, may need to change original course line because of certain reason, at this moment aircraft will enter the turning stage, and its turning can be approximately uniform circular motion; Before aircraft is turned, according to random selected horizontal movement direction
Figure FSA0000094734950000025
With vertical motion direction θ ζ, make uniform circular motion, in this stage, the speed v of aircraft ζ=v γ, the turning stage continues timeslot number ζ ∈ U[ζ min, ζ max]; After the γ stage finishes, aircraft may enter the turning stage, also may enter the constant speed decline η stage because being about to arrive destination; After supposing that each smooth flight finishes, the probability that aircraft enters the turning stage is P γ-ζ, the probability that enters the constant speed decline stage is P γ-η, obvious P γ-ζ+ P γ-η=1; Suppose that aircraft this process before fall-retarding will carry out K time, K=0,1,2 ..., in this case, the γ stage carries out K+1 time, and the ζ stage carries out K time; The desired value of K is as follows:
E [ K ] = &Sigma; K = 0 &infin; K ( P &gamma; - &xi; ) K ( 1 - P &gamma; - &xi; ) = P &gamma; - &xi; 1 - P &gamma; - &xi; - - - ( 2 )
(5) constant speed descends (η stage)
When being about to arrive destination, enter this stage; In this stage, aircraft is with at the uniform velocity straight line decline of fixed speed, and initial velocity magnitude and direction are respectively: v η: v γ,
Figure FSA0000094734950000023
θ η∈ U[pi/2, π]; Continue timeslot number η ∈ U[η min, η max];
(6) slow down and stop (d stage)
After arriving ground, aircraft will be decelerated to and stop, and can regard the motion in this stage as uniformly retarded motion; Initial velocity magnitude and direction are respectively: v d=v η,
Figure FSA0000094734950000031
θ d=pi/2; Duration d ∈ U[d min, d max], can obtain accordingly acceleration and be
a d = 0 - v &eta; d&Delta;t = - v &eta; d&Delta;t - - - ( 3 )
(7) static (p stage)
The duration p ∈ U[p that aircraft remains static min, p max], then start flight next time.
3. a kind of aviation self-organizing network wireless link stability prediction method according to claim 1 is characterized in that in each state, the concrete of node movement rate probability density function determines that method is:
For constant speed rising, smooth flight, turning and constant speed this four-stage that descends, because the node movement rate all remains unchanged, so the node movement rate distributes to meet and is uniformly distributed, as shown in the formula:
f ( v ) = 1 v max - v min v &Element; ( v min , v max ) 0 else - - - ( 4 )
Accelerating takeoff phase, node is with initial velocity v 0=O, evenly accelerate to given target velocity v α, v αEvenly choose from [v min, v max]; At the deceleration stop phase, node is with initial velocity v α, evenly decelerate to O, v αEvenly choose from [v min, v max]; Therefore, accelerating takeoff phase and deceleration stop phase, the probability density function of node movement rate is:
f &alpha; ( v ) = 2 v max + v min , 0 &le; v &le; v min 2 v max + v min ( 1 - v - v min v max - v min ) , v min < v &le; v max - - - ( 5 )
Can obtain node movement rate probability density function in each state by said process.
4. a kind of aviation self-organizing network wireless link stability prediction method according to claim 1 is characterized in that the concrete of probability density function of relative speed between node determines that method is:
The speed of supposing optional two mobile nodes is respectively v 1, v 2, the relative velocity v ' of these two mobile nodes=v 1-v 2, angle ω ∈ [0, π], establish v 1, v 2, v ' mould be respectively v 1, v 2, v ', have
v &prime; = v 1 2 + v 2 2 - 2 v 1 v 2 cos w , w = arccos v 1 2 + v 2 2 - v &prime; 2 2 v 1 v 2 &CenterDot; - - - ( 6 )
Due to v 1, v 2, ω is separate, so, v 1, v 2, ω joint probability density function can be written as
f w , v 1 , v 2 ( w , v 1 , v 2 ) = 1 &pi; f v 1 ( v 1 ) f v 2 ( v 2 ) &CenterDot; - - - ( 7 )
Convert by Jacobi, can obtain v 1, v 2, v ' joint probability density function
f v &prime; , v 1 , v 2 ( v &prime; , v 1 , v 2 ) = f w , v 1 , v 2 ( w , v 1 , v 2 ) | &PartialD; w &PartialD; v &prime; |
= 2 v &prime; f w , v 1 , v 2 ( w , v 1 , v 2 ) 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2 - - - ( 8 )
= 2 v &prime; f v 1 ( v 1 ) f v 2 ( v 2 ) &pi; 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2 &CenterDot;
The probability density function that can be drawn relative speed v ' by formula (8) is
f v &prime; ( v &prime; ) = &Integral; v 2 , min v 2 , max &Integral; v 1 , min v 1 , max 2 v &prime; f v 1 ( v 1 ) f v 2 ( v 2 ) &pi; 4 v 1 2 v 2 2 - ( v 1 2 + v 2 2 - v &prime; 2 ) 2 dv 1 dv 2 - - - ( 9 )
Definition status S set={ α, beta, gamma, ζ, η, d, p}; Then, the rate distribution of each state obtained according to formula (4) and (5), in conjunction with formula (9), the probability density function that can obtain the relative speed under any two states is
f v &prime; E , G ( v &prime; ) = &Integral; v G min v G max &Integral; v E min v E max 2 v &prime; f ( v E ) f ( v G ) &pi; 4 v E 2 v G 2 - ( v E 2 + v G 2 - v &prime; 2 ) 2 dv E dv G - - - ( 10 )
Wherein, E, G all belongs to S;
Definition p (E), p (G) is respectively node constantly in E, the probability in G stage, the node relative speed distributes and can be expressed as:
f v &prime; ( v &prime; ) = &Sigma; E &Element; S &Sigma; G &Element; S p ( E ) p ( G ) f v &prime; E , G ( v &prime; ) - - - ( 11 )
Can obtain the probability density function of node movement rate in each state by said process.
5. a kind of aviation self-organizing network wireless link stability prediction method according to claim 1 is characterized in that between node, the concrete of wireless link duration distribution function determines that method is:
Consider the wireless link between node 1 and node 2, using node 2 as with reference to node, node 1 and node 2 apart from being d 0, direction and the d of two node speed of related movement v ' 0Between angle be θ, R is the node transmission range; Use geometric knowledge, obtain the distribution function of link duration
F ( t ) = P ( T &le; t ) = P ( d 1 v &prime; &le; t ) = P ( d 0 cos &theta; + R 2 + d 0 2 sin 2 &theta; v &prime; &le; t )
= &Integral; &Integral; &Integral; d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t f ( &theta; , d 0 , v &prime; ) d&theta;d d 0 d v &prime; - - - ( 12 )
θ, d 0, tri-variablees of v ' are separate, therefore formula (12) can be expressed as
F ( t ) = &Integral; &Integral; &Integral; d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t f ( &theta; ) f ( d 0 ) f v &prime; ( v &prime; ) d&theta;d d 0 d v &prime; - - - ( 13 )
θ obeys and is uniformly distributed at interval [0, π], and its probability density function is
f ( &theta; ) = 1 &pi; - - - ( 14 )
By
F ( d 0 ) = P ( d &le; d 0 ) = 4 3 &pi; d 0 3 4 3 &pi; R 3 = d 0 2 R 3 - - - ( 15 )
Can obtain
f ( d 0 ) = 3 d 0 2 R 3 - - - ( 16 )
By qualifications
d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; &le; t - - - ( 17 )
Can obtain
cos &theta; &le; v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; t d 0 - - - ( 18 )
Order
k=∫f(θ)dθ (19)
{。##.##1},
k = 0 v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; t d 0 &le; - 1 &pi; - arccos v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; t d 0 &pi; - 1 < v &prime; 2 t 2 + d 0 2 - R 2 2 v &prime; t d 0 < 1 1 v &prime; 2 + t 2 + d 0 2 - R 2 2 v &prime; t d 0 &GreaterEqual; 1 - - - ( 20 )
Therefore, formula (13) can be expressed as
F ( t ) = &Integral; 0 2 v &prime; max &Integral; 0 R kf ( d 0 ) f v &prime; ( v &prime; ) d d 0 d v &prime; - - - ( 21 )
By formula (11), (14), (16), (20) substitution formula (21), can try to achieve the expression formula of link duration distribution function F (t).
6. a kind of aviation self-organizing network wireless link stability prediction method according to claim 1 is characterized in that between node, the concrete of wireless link duration distribution function determines that method is:
Consider the wireless link between node 1 and node 2, using node 2 as with reference to node, node 1 and node 2 apart from being d 0, direction and the d of two node speed of related movement v ' 0Between angle be θ, R is the node transmission range; The duration T of wireless link between node 1 and node 2 12For
T 12 = d 0 cos &theta; + R 2 - d 0 2 sin 2 &theta; v &prime; - - - ( 22 )
By T 12Be to compare 0.9 o'clock corresponding link duration with link duration distribution function F (t) value, obtain
Figure FSA0000094734950000072
T ^ 12 = T 12 T F ( t ) = 0.9 - - - ( 23 )
The present invention is defined as the wireless link stability factor of 2 of node 1 and nodes
S 12 = min ( T ^ 12 , 1 ) - - - ( 24 )
For the multi-hop wireless link be comprised of (N-1) bar one hop link, its wireless link stability factor is defined as the product of many one hop link stability factors,
S 1N=S 12S 23S 34…S (N-1)N (25)
The wireless link stability factor can be used for predicting the stability state of current wireless link, in actual applications, determine F (t) according to the actual value of each motion state parameters in the aircraft mobility model, then calculated the stability factor of wireless link between node by formula (25), and the important evidence using it as Route Selection, thereby set up the wireless transmission link that stability is high for the aviation self-organized network nodes.
CN201310397888.2A 2013-09-05 2013-09-05 A kind of aviation self-organizing network wireless link stability prediction method Active CN103428747B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310397888.2A CN103428747B (en) 2013-09-05 2013-09-05 A kind of aviation self-organizing network wireless link stability prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310397888.2A CN103428747B (en) 2013-09-05 2013-09-05 A kind of aviation self-organizing network wireless link stability prediction method

Publications (2)

Publication Number Publication Date
CN103428747A true CN103428747A (en) 2013-12-04
CN103428747B CN103428747B (en) 2016-02-24

Family

ID=49652759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310397888.2A Active CN103428747B (en) 2013-09-05 2013-09-05 A kind of aviation self-organizing network wireless link stability prediction method

Country Status (1)

Country Link
CN (1) CN103428747B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105320787A (en) * 2014-07-22 2016-02-10 张琳 Method for constructing dynamic evaluation model of network link reliability
CN105912749A (en) * 2016-03-31 2016-08-31 北京润科通用技术有限公司 Simulation method and device
CN106209202A (en) * 2016-06-27 2016-12-07 东南大学 A kind of aviation self-organized network topology builds and the method for linking Internet
CN106658537A (en) * 2016-12-07 2017-05-10 北京大学深圳研究生院 Method and system for establishing 3-dimension spatio temporally correlated group mobility model of unmanned aerial vehicle mobile ad hoc network node
CN109104366A (en) * 2018-08-29 2018-12-28 重庆大学 A kind of link duration calculation method, MPR node selecting method and Routing Protocol
CN109808915A (en) * 2018-12-28 2019-05-28 中国人民解放军63653部队 A kind of a wide range of flight test multistation relay TT & C system of long endurance and method
CN110417588A (en) * 2019-07-19 2019-11-05 北京科技大学 A kind of aviation dynamic network paths planning method based on Game with Coalitions

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101848557A (en) * 2009-05-31 2010-09-29 北京航空航天大学 Method for constructing mobile scale-free self-organized network model
US20110057830A1 (en) * 2009-09-10 2011-03-10 The Boeing Company Method for validating aircraft traffic control data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101848557A (en) * 2009-05-31 2010-09-29 北京航空航天大学 Method for constructing mobile scale-free self-organized network model
US20110057830A1 (en) * 2009-09-10 2011-03-10 The Boeing Company Method for validating aircraft traffic control data

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JING LI, LEI LEI, WEIKANG LIU, YUPING SHEN AND GANG ZHU: "An Improved Semi-Markov Smooth Mobility Model for Aeronautical Ad hoc Networks", 《IEEE》 *
郑 博,张衡阳,黄国策: "基于3-D SGM 移动模型的MANET链路动态性", 《吉林大学学报》 *
郑 博,张衡阳,黄国策: "航空自组网链路动态性分析", 《系统工程与电子技术》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105320787A (en) * 2014-07-22 2016-02-10 张琳 Method for constructing dynamic evaluation model of network link reliability
CN105912749A (en) * 2016-03-31 2016-08-31 北京润科通用技术有限公司 Simulation method and device
CN105912749B (en) * 2016-03-31 2019-06-04 北京润科通用技术有限公司 Emulation mode and device
CN106209202A (en) * 2016-06-27 2016-12-07 东南大学 A kind of aviation self-organized network topology builds and the method for linking Internet
CN106209202B (en) * 2016-06-27 2018-11-06 东南大学 A kind of method of aviation self-organized network topology structure and linking Internet
CN106658537A (en) * 2016-12-07 2017-05-10 北京大学深圳研究生院 Method and system for establishing 3-dimension spatio temporally correlated group mobility model of unmanned aerial vehicle mobile ad hoc network node
CN109104366A (en) * 2018-08-29 2018-12-28 重庆大学 A kind of link duration calculation method, MPR node selecting method and Routing Protocol
CN109104366B (en) * 2018-08-29 2020-10-16 重庆大学 Link duration calculation method, MPR node selection method and routing protocol
CN109808915A (en) * 2018-12-28 2019-05-28 中国人民解放军63653部队 A kind of a wide range of flight test multistation relay TT & C system of long endurance and method
CN109808915B (en) * 2018-12-28 2022-02-01 中国人民解放军63660部队 Multi-station relay measurement and control system and method for long-endurance large-range flight test
CN110417588A (en) * 2019-07-19 2019-11-05 北京科技大学 A kind of aviation dynamic network paths planning method based on Game with Coalitions
CN110417588B (en) * 2019-07-19 2020-10-20 北京科技大学 Aviation dynamic network path planning method based on alliance game

Also Published As

Publication number Publication date
CN103428747B (en) 2016-02-24

Similar Documents

Publication Publication Date Title
CN103428747B (en) A kind of aviation self-organizing network wireless link stability prediction method
Alalewi et al. On 5G-V2X use cases and enabling technologies: A comprehensive survey
CN111786713B (en) Unmanned aerial vehicle network hovering position optimization method based on multi-agent deep reinforcement learning
Lei et al. Impact of packet loss on CACC string stability performance
CN109286913B (en) Energy consumption optimization method of unmanned aerial vehicle mobile edge computing system based on cellular network connection
CN110650039B (en) Multimodal optimization-based network cooperative communication model for unmanned aerial vehicle cluster auxiliary vehicle
CN105979558B (en) A kind of unmanned plane trunked communication system based on cloud
Prathiba et al. SDN-assisted safety message dissemination framework for vehicular critical energy infrastructure
Lin et al. A 3D smooth random walk mobility model for FANETs
Dai et al. UAV-aided wireless communication design with energy constraint in space-air-ground integrated green IoT networks
CN113179116A (en) Joint optimization method for transmission power and flight trajectory of unmanned aerial vehicle relay system
Seliem et al. Delay analysis for drone-based vehicular ad-hoc networks
Wang et al. A-GR: A novel geographical routing protocol for AANETs
CN100536619C (en) Clustering method for air mobile node and air-ground integrated clustering system
Husain et al. Simulated analysis of location and distance based routing in VANET with IEEE802. 11p
Wu et al. FM-MAC: A multi-channel MAC protocol for FANETs with directional antenna
Zhou et al. Research on multi-UAV networks in disaster emergency communication
Zhou et al. A topology aware routing protocol based ADS-B system for aeronautical ad hoc networks
CN116187021A (en) Unmanned airport surface operation scheduling method and application
Han et al. Deep reinforcement learning assisted spectrum management in cellular based urban air mobility
Yue et al. A Channel Knowledge Map-Aided Personalized Resource Allocation Strategy in Air-Ground Integrated Mobility
Jiang et al. Research on OLSR adaptive routing strategy based on dynamic topology of UANET
Newton et al. Analysis of topology algorithms for commercial airborne networks
Taya et al. Coverage expansion through dynamic relay vehicle deployment in mmwave V2I communications
CN112367632A (en) Networked measurement and control system suitable for stratospheric unmanned aerial vehicle and measurement and control method thereof

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20170118

Address after: 210000 floor, block B, LIAN building, No. 18, Hui Hui Road, Gulou District, Jiangsu, Nanjing, China

Patentee after: CHINA COMMUNICATION TECHNOLOGY (NANJING) Co.,Ltd.

Address before: Yudaojie Nanjing 210016 Jiangsu province No. 29

Patentee before: Nanjing University of Aeronautics and Astronautics

DD01 Delivery of document by public notice

Addressee: Yuan Dong

Document name: Notification of Passing Examination on Formalities

DD01 Delivery of document by public notice
TR01 Transfer of patent right

Effective date of registration: 20211229

Address after: 430000 room 02, 4 / F, building A3, Wuhan International Enterprise Center, No. 1, Guanshan Second Road, Donghu New Technology Development Zone, Wuhan, Hubei Province

Patentee after: WUHAN HUAXUN GUORONG TECHNOLOGY CO.,LTD.

Address before: 210000 floor 10, block B, Lianchuang building, No. 18, Jihui Road, Gulou District, Nanjing, Jiangsu Province

Patentee before: CHINA COMMUNICATION TECHNOLOGY (NANJING) Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240416

Address after: No. 200 Nanneihuan Street, Xiaodian District, Taiyuan City, Shanxi Province, 030000

Patentee after: Shanxi Aerospace Industry Research Institute Group Co.,Ltd.

Country or region after: China

Address before: 430000 room 02, 4 / F, building A3, Wuhan International Enterprise Center, No. 1, Guanshan Second Road, Donghu New Technology Development Zone, Wuhan, Hubei Province

Patentee before: WUHAN HUAXUN GUORONG TECHNOLOGY CO.,LTD.

Country or region before: China

TR01 Transfer of patent right