CN114567400A - Unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness - Google Patents
Unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness Download PDFInfo
- Publication number
- CN114567400A CN114567400A CN202210199505.XA CN202210199505A CN114567400A CN 114567400 A CN114567400 A CN 114567400A CN 202210199505 A CN202210199505 A CN 202210199505A CN 114567400 A CN114567400 A CN 114567400A
- Authority
- CN
- China
- Prior art keywords
- unmanned aerial
- aerial vehicle
- antenna
- time
- angle
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000009826 distribution Methods 0.000 claims description 15
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 14
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/18502—Airborne stations
- H04B7/18506—Communications with or from aircraft, i.e. aeronautical mobile service
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention provides an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness. Compared with the prior art, the model provided by the modeling method is a three-dimensional double-cylinder model, the model comprehensively considers the scattering of the unmanned aerial vehicle side, the reflection and scattering of the ground station side, the mobility of the unmanned aerial vehicle transmitting end and the ground station receiving end, particularly the influence of the rotation of the unmanned aerial vehicle side and the three-dimensional moving track of the unmanned aerial vehicle, and accords with various actual communication scenes; the time-varying distance and the time-varying angle are calculated, the complex envelope signal is received, the actual communication situation and the non-stable statistical characteristics of the unmanned aerial vehicle space MIMO channel can be accurately described, based on the modeling method, some general rules of the unmanned aerial vehicle communication scheme can be summarized through actual simulation, and an idea is provided for building a robust unmanned aerial vehicle wireless communication system.
Description
Technical Field
The invention belongs to the field of wireless channel modeling, and particularly relates to an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness.
Background
Along with the development of unmanned aerial vehicles in military and civil fields, the unmanned aerial vehicle communication network that the information exchange formed between unmanned aerial vehicle and the unmanned aerial vehicle or between unmanned aerial vehicle and other receiving stations constitutes the important component part of unmanned aerial vehicle system. Meanwhile, the unmanned aerial vehicle plays an important role in the application of the fifth generation (5G) mobile communication system, the unmanned aerial vehicle communication technology has become one of the key points of the research of the wireless communication technology, and the establishment of an accurate and reliable unmanned aerial vehicle communication channel model is particularly important.
At present, an unmanned aerial vehicle air-ground channel modeling method mainly comprises a deterministic model, a statistical model and a geometric random channel model, but the deterministic model modeling process is complex, completely depends on detailed information of a channel environment, and has no universality; the random parameters obtained by the statistical model according to experience have certain difference with the actual scene, the accuracy is low, some unique characteristics of unmanned aerial vehicle communication cannot be captured, and the geometric random channel model balances the complexity and the accuracy, and becomes the mainstream unmanned aerial vehicle air-ground channel modeling method. However, most of the existing unmanned aerial vehicle air-ground channel models based on geometric randomness assume that an unmanned aerial vehicle does uniform linear motion in a three-dimensional space and neglects the influence caused by the position change of an array antenna, so that the unmanned aerial vehicle three-dimensional space model does not conform to the characteristics of any flight track and rotating antenna, and meanwhile neglects channel non-stationarity caused by the movement of the unmanned aerial vehicle and a ground station.
Disclosure of Invention
The invention aims to provide an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness, which aims to solve the technical problems that the existing unmanned aerial vehicle air-ground channel model based on geometric randomness mostly assumes that an unmanned aerial vehicle does uniform linear motion in a three-dimensional space and ignores the influence caused by the position change of an array antenna, the unmanned aerial vehicle air-ground channel model does not conform to the characteristics that the actual unmanned aerial vehicle three-dimensional space has any flight track and rotating antenna, and meanwhile, the channel instability caused by the movement of the unmanned aerial vehicle and a ground station is ignored.
In order to solve the technical problems, the specific technical scheme of the invention is as follows:
an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness comprises the following steps:
step 1, constructing a three-dimensional model of an unmanned aerial vehicle air-ground MIMO channel, wherein the three-dimensional model of the unmanned aerial vehicle air-ground MIMO channel comprises an unmanned aerial vehicle which moves at a high speed in a low-altitude three-dimensional space and serves as a transmitting end, a ground station which moves in a ground two-dimensional space and serves as a receiving end, and two right cylinders which are positioned on the unmanned aerial vehicle side and the ground station side.
The three-dimensional movement of the unmanned aerial vehicle is determined by the speed vTHorizontal angle of velocity gammaTAnd speed pitch angle ξTCharacterizing, and satisfying the Gauss Markov motion model with an initial position of OT,tiTime shift to OT’;
The moving speed of the ground station is vRHorizontal angle of velocity of gammaRInitial position is OR,tiTime shift to OR’;
The radius of the right cylinder at the unmanned aerial vehicle side is RTSurface distribution N1An effective scatterer, wherein the n-th1An effective scatterer can be represented asThe radius of the right cylinder at the side of the ground station is RRSurface distribution N2An effective scatterer, wherein the n-th2An effective scatterer can be represented as Bottom surface distribution N3A ground reflection effective scatterer, wherein3An effective scatterer can be expressed as
Wherein, the initial horizontal distance between unmanned aerial vehicle and the ground station is D, the vertical distance is H, and the pitch angle of unmanned aerial vehicle position satisfies beta0H/D, the ground station is at a height H0。
Step 2, use of LT×LRThe antenna array matrix describes an unmanned aerial vehicle air-ground MIMO channel, the antenna array is a uniform linear antenna array, and an unmanned aerial vehicle side antenna array TXAntenna rotation pitch angle psiTAnd the rotation azimuth angle thetaTAre all time-varying angles, and the ground station side antenna array RXAntenna pitch angle psiRAnd azimuth angle thetaRFixing;
step 3, deducing nonstationary time-varying distance and time-varying angle of the unmanned aerial vehicle air-ground MIMO channel caused by unmanned aerial vehicle antenna rotation and unmanned aerial vehicle three-dimensional movement according to the geometric relationship;
and 4, calculating the total receiving complex envelope signal of the unmanned aerial vehicle space-ground MIMO channel according to different path component signals of the unmanned aerial vehicle side antenna propagation signal propagated to the ground station side antenna through different scatterers.
Further, the gaussian markov motion model satisfied by the three-dimensional motion of the unmanned aerial vehicle in step 1 is:
wherein the velocity is t0The initial value of the time is vT(t0) At a pitch angle t0The initial value of the time is xiT(t0) Horizontal angle at t0The initial value of time is gammaT(t0) (ii) a Velocity magnitude at tiThe value of the time is vT(ti) Elevation angle at tiThe value of time is xiT(ti) Horizontal angle at tiThe value of time is gammaT(ti) (ii) a When i → infinity, the asymptotic mean of the velocity magnitudes isAsymptotic mean value of pitch angle ofAsymptotic mean of horizontal angle isρv、 ρξAnd ργIs [0,1 ]]Inner tuning value, pvRandomness, p, characterizing the magnitude of the velocityξCharacterizing the randomness of the pitch angle, pγCharacterizing randomness of horizontal angles; l, M and N are variables that follow a Gaussian distribution; in order to ensure that the speed, the horizontal angle and the pitch angle are in reasonable ranges, the maximum value of the speed is set as vTmaxThe maximum value of the pitch angle is set to xiTmaxThe maximum horizontal angle is set to gammaTmax。
Further, in step 2, the unmanned aerial vehicle side antenna array TXThe rotation time-varying angle of the antenna meets the cosine-remaining process:
ψT(ti)=ψT(t0)+ψTm cos(πti)
θT(ti)=θT(t0)+θTm cos(πti)
wherein the antenna rotates at a pitch angle t0The initial value of the time is psiT(t0) The antenna is rotated at a horizontal angle t0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated by a horizontal angle tiValue of time of dayIs thetaT(ti) (ii) a The maximum amplitude of the variation of the antenna rotation pitch angle is psiTmThe maximum amplitude of the variation of the horizontal angle of rotation of the antenna is thetaTm。
Further, unmanned aerial vehicle side antenna array TXThe time-varying angle of rotation of the antenna satisfies a Gaussian Markov process:
wherein the antenna has a rotation pitch angle of t0The initial value of the time is psiT(t0) At the antenna rotation pitch angle t0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated by a horizontal angle tiThe value of time is thetaT(ti) (ii) a When i → ∞, the asymptotic mean of the antenna rotation pitch angle isThe asymptotic mean value of the antenna rotation horizontal angle isρψAnd ρθIs [0,1 ]]Inner tuning value, pψCharacterizing the randomness, rho, of the rotational pitch angle of the antennaθCharacterizing randomness of a horizontal angle of rotation of the antenna; x and Y are variables subject to Gaussian distribution; in order to ensure that the rotation pitch angle and the horizontal angle of the antenna are within a reasonable range, the maximum value of the rotation pitch angle of the antenna is set to psiTmaxThe maximum value of the horizontal angle of the line rotation is set to thetaTmax。
Further, in step 3, the non-stationary time-varying distance is a distance from the p-th antenna in the antenna array of the unmanned aerial vehicle side to the q-th antenna in the antenna array of the ground station side through different effective scatterersFormed linkAndtime-varying distance epsilonpq(ti)、 Andthe path components formed by the links are LOS, SB1, SB2, SB3 and DB, and the time-varying distance is tiThe position coordinates of the unmanned aerial vehicle and the ground station at the moment are calculated by a distance formula, tiTime of dayAndrepresenting the three-dimensional position coordinates of the transmitting end of the unmanned aerial vehicle and the receiving end of the ground station, wherein the expression of the coordinates is
Wherein, DeltaTIs the distance, delta, between the p-th antenna of the drone and the center of the drone antenna arrayRIs the distance between the q-th antenna of the ground station and the center of the ground station antenna array, and meets the requirement
δTIndicating the antenna spacing of the drone, deltaRRepresenting the antenna spacing of the ground station.
Further, in step 3, the non-stationary time-varying angle includes tiThe departure azimuth, departure pitch, arrival azimuth and arrival pitch of the LOS path at time tiThe departure azimuth angle of the LOS path at the momentDenotes, tiThe angle of departure and pitch of the LOS path at the momentDenotes, tiAzimuth of arrival of time LOS pathDenotes, tiThe angle of elevation of arrival of the LOS path at timeIt is shown that,azimuth of departure of linkIt is shown that,the departure and pitch angle of the link is set byIt is shown that,link arrival azimuth angle ofIt is shown that,angle of arrival of link fromIs represented by, k is 1, 2, 3, wherein
Further, in step 4, an expression of each path component signal of the receiving complex packet is obtained by considering the non-stationary time-varying angle and the time-varying distance:
the complex envelope expression of the LoS path component is
The complex envelope of the SB1 path component is expressed as
The complex envelope of the SB2 path component is expressed as
The complex envelope of the SB3 path component is expressed as
The complex envelope expression of the DB path component is
Wherein omegapqRepresenting the total received power, K representing the Rice factor, etaSB1、ηSB2、ηSB3And ηDBRespectively representing the total scattered power omega of each path componentpqA ratio of/(K +1) and satisfies ηSB1+ηSB2+ ηSB2+ηDB=1,Andrepresenting the phase produced by each path component through the scatterer as an independent random variable obeying a uniform distribution over [ - π, π ]TmRepresenting the maximum Doppler frequency, f, of the droneRmRepresenting the maximum doppler frequency of the ground station.
The invention discloses an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness, which has the following advantages:
the invention provides an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness, wherein the model is a three-dimensional double-cylinder model, and compared with the prior art, the model comprehensively considers the scattering of an unmanned aerial vehicle side, the reflection and scattering of a ground station side, the mobility of an unmanned aerial vehicle transmitting end and a ground station receiving end, particularly the rotation of the unmanned aerial vehicle side and the influence of the three-dimensional moving track of the unmanned aerial vehicle, and accords with various actual communication scenes; the time-varying distance and the time-varying angle are calculated, the complex envelope signal is received, the actual communication situation and the non-stable statistical characteristic of the unmanned aerial vehicle space MIMO channel can be accurately described, based on the modeling method, some general rules of the unmanned aerial vehicle communication scheme can be summarized through actual simulation, and an idea is provided for building a robust unmanned aerial vehicle communication wireless communication system.
Drawings
FIG. 1 is a schematic diagram of an unmanned aerial vehicle space-ground MIMO channel model based on geometric randomness according to the present invention;
fig. 2 is a schematic diagram of the rotation of the drone antenna array of the present invention;
fig. 3 is a diagram illustrating the components of each path formed in the channel model of the present invention.
Detailed Description
In order to better understand the purpose, structure and function of the present invention, the following describes the unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness in further detail with reference to the accompanying drawings.
The technical scheme adopted by the invention is as follows: an unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness comprises the following steps:
step 1, constructing a three-dimensional model of the unmanned aerial vehicle air-ground MIMO channel, wherein the three-dimensional model of the unmanned aerial vehicle air-ground MIMO channel comprises an unmanned aerial vehicle which moves at a high speed in a low-altitude three-dimensional space and serves as a transmitting end, a ground station which moves in a ground two-dimensional space and serves as a receiving end, and two right cylinders which are positioned on the side of the unmanned aerial vehicle and the side of the ground station, as shown in FIG. 1.
The three-dimensional movement of the unmanned aerial vehicle is controlled by the speed vTHorizontal angle of velocity gammaTAnd speed pitch angle ξTCharacterizing, and satisfying the Gaussian Markov motion model with an initial position of OT,tiTime of day shifts to OT’;
The moving speed of the ground station is vRHorizontal angle of velocity of gammaRInitial position is OR,tiTime shift to OR’;
The radius of the right cylinder at the unmanned aerial vehicle side is RTSurface distribution N1An effective scatterer, wherein the n-th1An effective scatterer can be represented asThe radius of the right cylinder at the side of the ground station is RRSurface distribution N2An effective scatterer, wherein2An effective scatterer can be represented as Bottom surface distribution N3A ground reflection effective scatterer, wherein3Effective powderThe projectile may be represented as
Wherein, the initial horizontal distance between unmanned aerial vehicle and the ground station is D, the vertical distance is H, and the pitch angle of unmanned aerial vehicle position satisfies beta0H/D, the ground station is at a height H0。
In order to better represent the influence of the three-dimensional movement of the unmanned aerial vehicle on the non-stationary characteristics of the channel, the three-dimensional movement modeling of the unmanned aerial vehicle is a Gaussian Markov movement model and is expressed as follows:
wherein the velocity is t0The initial value of the time is vT(t0) At a pitch angle t0The initial value of the time is xiT(t0) Horizontal angle at t0The initial value of time is gammaT(t0) The velocity is tiThe value of the time is vT(ti) Elevation angle at tiThe value of time is xiT(ti) Horizontal angle at tiThe value of time is gammaT(ti) (ii) a When i → ∞ is reached, the asymptotic mean value of the velocity magnitude is μvTThe asymptotic mean value of the pitch angle isAsymptotic mean of horizontal angle ofρv、ρξAnd ργIs [0,1 ]]Inner tuning value, pvRandomness, p, characterizing the magnitude of the velocityξCharacterizing randomness of pitch angle, pγCharacterizing randomness of horizontal angles; l, M and N are variables that follow a Gaussian distribution; in order to ensure that the speed, the horizontal angle and the pitch angle are in reasonable ranges, the maximum value of the speed is set as vTmaxThe maximum value of the pitch angle is set to xiTmaxThe maximum horizontal angle is set to gammaTmax。
Step 2, use LT×LRThe antenna array matrix describes the unmanned aerial vehicle air-ground MIMO channel, which is a uniform linear antenna array, assuming LT=LRAs shown in fig. 2, the drone side antenna array T takes into account the roll and pitch motion of the droneXAntenna rotation pitch angle psiTAnd rotation azimuth angle thetaTAre all time-varying angles, and the ground station side antenna array RXAntenna pitch angle psiRAnd azimuth angle thetaRFixing;
the unmanned aerial vehicle side antenna array TXThe antenna rotation time-varying angle is modeled as a cosine process:
ψT(ti)=ψT(t0)+ψTm cos(πti)
θT(ti)=θT(t0)+θTm cos(πti)
wherein the antenna rotates at a pitch angle t0The initial value of the time is psiT(t0) The antenna is rotated at a horizontal angle t0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated at a horizontal angle tiThe value of time is thetaT(ti) (ii) a The maximum amplitude of the variation of the antenna rotation pitch angle is psiTmThe maximum amplitude of the variation of the horizontal angle of rotation of the antenna is thetaTm。
When the rotation of the unmanned aerial vehicle antenna is random and unpredictable, the unmanned aerial vehicle side antenna array T is arrangedXThe antenna rotation time-varying angle is modeled as a Gauss horseThe Erkov process:
wherein the antenna rotates at a pitch angle t0The initial value of the time is psiT(t0) At the antenna rotation pitch angle t0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated at a horizontal angle tiThe value of time is thetaT(ti) (ii) a When i → ∞, the asymptotic mean of the antenna rotation pitch angle isThe asymptotic mean value of the antenna rotation horizontal angle isρψAnd ρθIs [0,1 ]]Inner tuning value, pψCharacterizing the randomness, rho, of the rotational pitch angle of the antennaθCharacterizing randomness of a horizontal angle of rotation of the antenna; x and Y are variables subject to Gaussian distribution; in order to ensure that the rotation pitch angle and the horizontal angle of the antenna are in a reasonable range, the maximum value of the rotation pitch angle of the antenna is set to be psiTmaxThe maximum value of the horizontal angle of linear rotation is set to thetaTmax。
And 3, deducing non-stationary time-varying distance and time-varying angle of the unmanned aerial vehicle air-ground MIMO channel caused by rotation of an unmanned aerial vehicle antenna and three-dimensional movement of the unmanned aerial vehicle according to the geometric relation.
The non-stationary time-varying distance is a link T formed by a pth antenna in the antenna array of the unmanned aerial vehicle side to a qth antenna in the antenna array of the ground station side through different effective scatterersp-Rq、 Andtime-varying distance epsilonpq(ti)、 Andas shown in FIG. 3, the path components formed by the links are LOS, SB1, SB2, SB3 and DB, and the time-varying distance is tiThe position coordinates of the unmanned aerial vehicle and the ground station at the moment are calculated by a distance formula, tiTime of dayAndrepresenting the three-dimensional position coordinates of the transmitting end of the unmanned aerial vehicle and the receiving end of the ground station, wherein the expression of the coordinates is
Wherein, DeltaTIs the distance, delta, between the p-th antenna of the drone and the center of the drone antenna arrayRIs the distance between the qth antenna of the ground station and the center of the antenna array of the ground station, and satisfies
δTIndicating the antenna spacing, delta, of the droneRIndicating the antenna spacing of the ground station.
The non-stationary time-varying angles due to the movement of the drone and the ground station include tiThe departure azimuth, departure pitch angle, arrival azimuth and arrival pitch angle, t, of the LOS path at timeiAzimuth of departure of time LOS pathDenotes tiThe departure-pitch angle of the LOS path at that momentDenotes, tiAzimuth of arrival of time LOS pathDenotes, tiThe angle of elevation of arrival of the LOS path at timeIt is shown that,the departure azimuth angle of the link is determined byIt is shown that the process of the present invention,the departure and pitch angle of the link is set byIt is shown that,link arrival azimuth angle ofIt is shown that,angle of arrival pitch of linkIs represented by, k is 1, 2, 3, wherein
And 4, calculating the total receiving complex envelope signal of the unmanned aerial vehicle space-ground MIMO channel according to different path component signals of the unmanned aerial vehicle side antenna propagation signal propagated to the ground station side antenna through different scatterers.
And obtaining an expression of each path component signal of the received complex envelope by considering the non-stationary time-varying angle and the time-varying distance:
(1) the complex envelope expression of the LoS path component is
(2) The complex envelope of the SB1 path component is expressed as
(3) The complex envelope of the SB2 path component is expressed as
(4) The complex envelope of the SB3 path component is expressed as
(5) The complex envelope expression of the DB path component is
Wherein omegapqRepresenting the total received power, K representing the Rice factor, etaSB1、ηSB2、ηSB3And ηDBRespectively representing the total scattered power omega of each path componentpqA ratio of (K +1) and satisfies ηSB1+ηSB2+ ηSB2+ηDB=1,Andrepresenting the phase generated by each path component through the scatterer as an independent random variable subject to uniform distribution over [ - π, π [ ], fTmRepresenting the maximum Doppler frequency, f, of the droneRmRepresenting the maximum doppler frequency of the ground station.
It is to be understood that the present invention has been described with reference to certain embodiments, and that various changes in the features and embodiments, or equivalent substitutions may be made therein by those skilled in the art without departing from the spirit and scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.
Claims (7)
1. An unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness is characterized by comprising the following steps:
step 1, constructing a three-dimensional model of an unmanned aerial vehicle air-ground MIMO channel, wherein the three-dimensional model of the unmanned aerial vehicle air-ground MIMO channel comprises an unmanned aerial vehicle which moves at a high speed in a low-altitude three-dimensional space and serves as a transmitting end, a ground station which moves in a ground two-dimensional space and serves as a receiving end, and two right cylinders which are positioned at the side of the unmanned aerial vehicle and the side of the ground station;
the three-dimensional movement of the unmanned aerial vehicle is controlled by the speed vTHorizontal angle of velocity gammaTAnd velocity pitch angle ξTCharacterizing, and satisfying the Gaussian Markov motion model with an initial position of OT,tiTime shift to OT’;
The moving speed of the ground station is vRHorizontal angle of velocity of gammaRInitial position is OR,tiTime shift to OR’;
The radius of the right cylinder at the unmanned aerial vehicle side is RTSurface distribution N1An effective scatterer, wherein1An effective scatterer is represented asThe radius of the right cylinder at the side of the ground station is RRSurface distribution N2An effective scatterer, wherein2An effective scatterer is represented as Bottom surface distribution N3A ground reflection effective scatterer, wherein3An effective scatterer is represented as
Wherein, the initial horizontal distance between unmanned aerial vehicle and the ground station is D, the vertical distance is H, and the pitch angle of unmanned aerial vehicle position satisfies beta0H/D, the ground station is at a height H0;
Step 2, use of LT×LRAntenna array matrix descriptionAn unmanned aerial vehicle air-ground MIMO channel, wherein the antenna array is a uniform linear antenna array, and an unmanned aerial vehicle side antenna array TXAntenna rotation pitch angle psiTAnd rotation azimuth angle thetaTAll of which are time-varying angles, ground station side antenna array RXAntenna pitch angle psiRAnd azimuth angle thetaRFixing;
step 3, deducing nonstationary time-varying distance and time-varying angle of the unmanned aerial vehicle air-ground MIMO channel caused by unmanned aerial vehicle antenna rotation and unmanned aerial vehicle three-dimensional movement according to the geometric relationship;
and 4, calculating the total receiving complex envelope signal of the unmanned aerial vehicle space-ground MIMO channel according to different path component signals of the unmanned aerial vehicle side antenna propagation signal propagated to the ground station side antenna through different scatterers.
2. The method for modeling the unmanned aerial vehicle air-ground MIMO channel based on geometric randomness according to claim 1, wherein the Gaussian Markov motion model satisfied by the three-dimensional motion of the unmanned aerial vehicle in step 1 is:
wherein the velocity is t0The initial value of the time is vT(t0) At a pitch angle t0The initial value of the time is xiT(t0) (ii) a Horizontal angle at t0The initial value of time is gammaT(t0) The velocity is tiThe value of the time is vT(ti) At a pitch angle tiThe value of time is xiT(ti) Water, waterStraight angle is at tiThe value of time is gammaT(ti) (ii) a When i → ∞ the asymptotic mean value of the velocity magnitude isAsymptotic mean value of pitch angle ofAsymptotic mean of horizontal angle isρv、ρξAnd ργIs [0,1 ]]Inner tuning value, pvRandomness, p, characterizing the magnitude of the velocityξCharacterizing the randomness of the pitch angle, pγCharacterizing randomness of horizontal angles; l, M and N are variables that follow a Gaussian distribution; in order to ensure that the speed, the horizontal angle and the pitch angle are in reasonable ranges, the maximum value of the speed is set as vTmaxThe maximum value of the pitch angle is set to xiTmaxThe maximum horizontal angle is set to gammaTmax。
3. The method of claim 2, wherein in step 2, the drone-side antenna array TXThe antenna rotation time variation angle satisfies the cosine process:
ψT(ti)=ψT(t0)+ψTmcos(πti)
θT(ti)=θT(t0)+θTmcos(πti)
wherein the antenna rotates at a pitch angle t0The initial value of the time is psiT(t0) The antenna is rotated at a horizontal angle t0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated at a horizontal angle tiThe value of time is thetaT(ti) (ii) a Maximum change in antenna rotation pitch angleAmplitude of psiTmThe maximum amplitude of the change of the antenna rotation horizontal angle is thetaTm。
4. The method of claim 2, wherein the drone side antenna array T is a space MIMO channel modeling method based on geometric randomnessXSatisfies the gaussian markov process:
wherein the antenna rotates at a pitch angle t0The initial value of the time is psiT(t0) At a rotational pitch angle t of the antenna0The initial value of the time is thetaT(t0) (ii) a The antenna rotates at a pitch angle tiThe value of time is psiT(ti) The antenna is rotated at a horizontal angle tiThe value of time is thetaT(ti) (ii) a When i → ∞, the asymptotic mean of the antenna rotation pitch angle isAsymptotic mean of the horizontal angle of rotation of the antenna isρψAnd ρθIs [0,1 ]]Inner tuning value, pψCharacterizing the randomness, rho, of the rotational pitch angle of the antennaθCharacterizing randomness of a horizontal angle of rotation of the antenna; x and Y are variables subject to a Gaussian distribution; in order to ensure that the rotation pitch angle and the horizontal angle of the antenna are in a reasonable range, the maximum value of the rotation pitch angle of the antenna is set to be psiTmaxThe maximum value of the horizontal angle of linear rotation is set to thetaTmax。
5. Root of herbaceous plantThe method as claimed in claim 3 or 4, wherein in step 3, the non-stationary time-varying distance is a link T formed by a p-th antenna in the unmanned aerial vehicle side antenna array to a q-th antenna in the ground station side antenna array via different effective scatterersp-Rq、Andtime-varying distance epsilonpq(ti)、Andthe path components formed by the links are LOS, SB1, SB2, SB3 and DB, and the time-varying distance is tiThe position coordinates of the unmanned aerial vehicle and the ground station at the moment are calculated by a distance formula, tiTime of dayAndrepresenting the three-dimensional position coordinates of the transmitting end of the unmanned aerial vehicle and the receiving end of the ground station, wherein the expression of the coordinates is
Wherein, DeltaTIs the distance, delta, between the p-th antenna of the drone and the center of the drone antenna arrayRIs the distance between the q-th antenna of the ground station and the center of the ground station antenna array, and meets the requirement
δTIndicating the antenna spacing, delta, of the droneRRepresenting the antenna spacing of the ground station.
6. The method of claim 5, wherein in step 3, the non-stationary time-varying angle comprises tiThe departure azimuth, departure pitch angle, arrival azimuth and arrival pitch angle, t, of the LOS path at timeiAzimuth of departure of time LOS pathDenotes, tiThe angle of departure and pitch of the LOS path at the momentDenotes, tiOf the LOS path at that momentArrival azimuth angle ofDenotes, tiThe angle of elevation of arrival of the LOS path at timeIt is shown that,azimuth of departure of linkIt is shown that,the departure and pitch angle of the link is set byIt is shown that the process of the present invention,link arrival azimuth angle ofIt is shown that,angle of arrival pitch of linkIs represented by, k is 1, 2, 3, wherein
7. The unmanned aerial vehicle space-earth MIMO channel modeling method based on geometric randomness as claimed in claim 6, wherein in step 4, an expression of each path component signal of the received complex envelope is obtained by considering non-stationary time-varying angle and time-varying distance:
the complex envelope expression of the LoS path component is
The complex envelope of the SB1 path component is expressed as
The complex envelope of the SB2 path component is expressed as
The complex envelope of the SB3 path component is expressed as
The complex envelope expression of the DB path component is
Wherein omegapqRepresenting the total received power, K representing the Rice factor, etaSB1、ηSB2、ηSB3And ηDBRespectively representing the total scattered power omega of each path componentpqA ratio of (K +1) and satisfies ηSB1+ηSB2+ηSB2+ηDB=1,Andrepresenting the phase generated by each path component through the scatterer as an independent random variable subject to uniform distribution over [ - π, π [ ], fTmRepresenting the maximum Doppler frequency, f, of the droneRmRepresenting the maximum doppler frequency of the ground station.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210199505.XA CN114567400B (en) | 2022-03-02 | 2022-03-02 | Unmanned aerial vehicle space MIMO channel modeling method based on geometric random |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210199505.XA CN114567400B (en) | 2022-03-02 | 2022-03-02 | Unmanned aerial vehicle space MIMO channel modeling method based on geometric random |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114567400A true CN114567400A (en) | 2022-05-31 |
CN114567400B CN114567400B (en) | 2024-06-07 |
Family
ID=81715524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210199505.XA Active CN114567400B (en) | 2022-03-02 | 2022-03-02 | Unmanned aerial vehicle space MIMO channel modeling method based on geometric random |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114567400B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150270885A1 (en) * | 2015-06-08 | 2015-09-24 | Donald C.D. Chang | MIMO Systems with Active Scatters and their Performance Evaluation |
CN109951213A (en) * | 2017-12-21 | 2019-06-28 | 上海交通大学 | High altitude platform MIMO three-dimensional geometry stochastic model method for building up and communication means |
CN113489560A (en) * | 2021-05-12 | 2021-10-08 | 东南大学 | Geometric random modeling method for non-stationary air-ground MIMO channel of unmanned aerial vehicle |
CN113938233A (en) * | 2021-11-18 | 2022-01-14 | 重庆邮电大学 | Geometric random modeling method for non-stationary air-to-air MIMO channel of unmanned aerial vehicle |
-
2022
- 2022-03-02 CN CN202210199505.XA patent/CN114567400B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150270885A1 (en) * | 2015-06-08 | 2015-09-24 | Donald C.D. Chang | MIMO Systems with Active Scatters and their Performance Evaluation |
CN109951213A (en) * | 2017-12-21 | 2019-06-28 | 上海交通大学 | High altitude platform MIMO three-dimensional geometry stochastic model method for building up and communication means |
CN113489560A (en) * | 2021-05-12 | 2021-10-08 | 东南大学 | Geometric random modeling method for non-stationary air-ground MIMO channel of unmanned aerial vehicle |
CN113938233A (en) * | 2021-11-18 | 2022-01-14 | 重庆邮电大学 | Geometric random modeling method for non-stationary air-to-air MIMO channel of unmanned aerial vehicle |
Non-Patent Citations (1)
Title |
---|
燕锋;周天香;李昊;庞井明;丁凯;夏玮玮;沈连丰: "一种无人机空地MIMO信道的三维几何随机模型", 东南大学学报(英文版), vol. 38, no. 004, 31 December 2022 (2022-12-31) * |
Also Published As
Publication number | Publication date |
---|---|
CN114567400B (en) | 2024-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112865897B (en) | Non-stationary channel simulation method and system for ground scene by unmanned aerial vehicle | |
CN108768565B (en) | Dynamic MIMO channel modeling and parameter calculation method for Internet of vehicles | |
CN100561250C (en) | The clutter base-band analog signal production method of spaceborne bistatic radar | |
CN113489560B (en) | Geometric random modeling method for non-stationary air-ground MIMO channel of unmanned aerial vehicle | |
CN113204015B (en) | Multi-unmanned aerial vehicle track planning method and system for distributed radar target simulation | |
CN114257293B (en) | Modeling method for interference of NGSO constellation on GEO satellite communication earth station based on space-time cut set | |
CN112070894A (en) | Real environment navigation multipath real-time simulation method, device, medium and electronic equipment | |
CN114970213B (en) | Sea clutter dynamic simulation method under soliton internal wave fluctuation effect | |
CN114124266B (en) | Channel modeling method based on IRS (intelligent resilient system) for assisting communication between unmanned aerial vehicle and unmanned ship | |
CN105044667A (en) | Double-satellite tracking method, device and system for moving target | |
CN113225711A (en) | Tunnel scene vehicle-to-vehicle MIMO wireless channel capacity estimation method | |
CN108614248B (en) | Frequency modulation detection target echo simulation method based on fractal theory | |
CN105388467A (en) | Method of correcting Doppler weather radar echo attenuation | |
CN113784314A (en) | Unmanned aerial vehicle data and energy transmission method assisted by intelligent reflection surface | |
CN114039684A (en) | Broadband UAV MIMO channel modeling method based on geometry | |
CN111551968B (en) | Unmanned aerial vehicle alignment system and method based on deep learning motion prediction | |
CN114567400A (en) | Unmanned aerial vehicle air-ground MIMO channel modeling method based on geometric randomness | |
CN104101869A (en) | Ground wave radar moving object simulation modeling method in polar coordinate | |
CN106019242A (en) | Space-based bistatic radar flight state configuration method | |
CN106405515B (en) | A kind of space based radar Clutter simulation method | |
CN114422060B (en) | Method and system for constructing unmanned aerial vehicle communication channel model | |
CN113938233B (en) | Geometric random modeling method for non-stationary space MIMO channel of unmanned aerial vehicle | |
CN113189585A (en) | Motion error compensation algorithm based on unmanned aerial vehicle bistatic SAR system | |
Ahmed et al. | 3D geometry‐based UAV to vehicle multiple‐input multiple‐output channel model incorporating Unmanned aerial vehicle heaving motion | |
CN109633580A (en) | A kind of space based radar Clutter simulation method and device considering Terrain Elevation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |