Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide an accurate unmanned aerial vehicle channel model building method based on a large-scale intelligent reflection unit, and the model building method can provide powerful support for the exploration of the key technology of a 6G communication system.
The technical scheme is as follows: the invention discloses an unmanned aerial vehicle channel model building method based on a large-scale intelligent reflection unit, which comprises the following steps:
s1, configuring the intelligent reflecting surface IRS on the surface of a building at the edge of a served cell, simulating a vertical building structure around a receiving end by using a three-dimensional cylinder, simulating the near field effect of the large-scale intelligent reflecting surface IRS by using second-order approximation of spherical wavefront, assuming that a scatterer is positioned on the surface of the three-dimensional cylinder, and establishing an unmanned aerial vehicle channel model based on the large-scale intelligent reflecting unit, wherein the intelligent reflecting surface IRS comprises intelligent reflecting units which are uniformly arranged; obtaining the complex channel gain of the channel according to the model;
s2, obtaining the complex channel gain according to the unmanned aerial vehicle channel model based on the large-scale intelligent reflection unit, wherein the complex channel gain comprises two components: the complex channel gain of the Unmanned Aerial Vehicle (UAV) which is directly transmitted with the receiving end without the Intelligent Reflector (IRS) and the complex channel gain of the UAV which is transmitted with the receiving end by the Intelligent Reflector (IRS);
s3, designing an optimization problem by maximizing the received signal power according to the unmanned aerial vehicle channel model based on the large-scale intelligent reflection unit;
s4, simplifying optimization problems: the optimization problem proposed in S3 is computationally very complex, so to reduce complexity, the problem needs to be further simplified; when the scale of the intelligent reflection unit in the intelligent reflection surface IRS is larger, the power of the received signal is mainly controlled by the reflection signal passing through the intelligent reflection surface IRS, and the complex channel gain of the reflection signal is mainly the direct component of the reflection signal, so that the process of solving the reflection phase is simplified;
s5, according to the simplified optimization problem, considering the time-varying Doppler frequency shift of the multipath component between the UAV antenna unit p and the receiving end antenna unit q after passing through the (m, n) -th intelligent reflection unit
When the optimal IRS reflection phase of the intelligent reflection surface is solved, the Doppler frequency shift of the direct component is subtracted to enhance the power of the received signal;
s6, solving a space-time correlation function based on the intelligent reflector IRS assistance through the complex channel gain obtained in the step S2 and the optimal intelligent reflector IRS reflection phase obtained in the step S5, and determining the influence of different parameters on the unmanned aerial vehicle channel characteristics through correlation analysis.
Further, the complex channel gain h obtained in step S1pq(t, τ), which is expressed as follows:
where t represents a time variable, L represents the number of taps, L represents the total number of taps, c
lDenotes the gain of the first tap, h
l,pq(t) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q directly without the smart reflector IRS,
represents the complex channel gain, τ, between the UAV antenna unit p and the receiving-end antenna unit q through the intelligent reflector IRS
l(t) represents the propagation delay of the ith tap, and δ (·) represents the impulse function.
Further, in step S2, the unmanned aerial vehicle UAV antenna unit p directly connects to the receiving-end antenna without passing through the intelligent reflection surface IRSComplex channel gain h between units ql,pq(t) is represented as follows:
wherein ,
representing the complex channel gain of the direct component between the unmanned aerial vehicle UAV antenna element p and the receiving end antenna element q,
a complex channel gain representing a scattering component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q;
wherein ,G
tDenotes the transmit antenna gain, G
rRepresenting the receiving-end antenna gain, gamma
TRRepresents the path loss, K, of the UAV to the receiving end
1Representing the rice factor, λ the carrier wavelength, t the time variable, π the circumference ratio, ξ
pq(t) represents the time-varying distance between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q,
the time-varying Doppler frequency shift of a direct component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q is shown, delta (l-1) represents a delay impact function after l taps, and N
lRepresenting scatterers
Number of (1), xi
pnl(t) denotes unmanned aerial vehicle UAV antenna unit p and scatterers
The time-varying distance between them,
representing scatterers
And the time-varying distance between the receiving-end antenna unit q,
representing the time-varying doppler shift of the scattered component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q.
Further, in step S2, the complex channel gain of the unmanned aerial vehicle UAV transmitted through the intelligent reflector IRS and the receiving end
Is represented as follows:
wherein ,
represents the complex channel gain of the direct component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q after being reflected by the intelligent reflecting surface IRS,
the complex channel gain of a scattering component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after being scattered by an intelligent reflecting surface IRS and a scatterer is represented;
wherein M represents a row position index of the intelligent reflection unit, N represents a column position index of the intelligent reflection unit, M represents a row reflection unit number of the intelligent reflection surface, N represents a column reflection unit number of the intelligent reflection surface, G
tRepresenting the gain of the transmitting antenna, G representing the gain of the IRS reflecting element, G
rRepresenting the receiving-end antenna gain, gamma
TIRRepresents the path loss, K, from the UAV to the IRS and then to the receiving end
2Representing the rice factor, λ the carrier wavelength, t the time variable, π the circumference ratio, ξ
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, theta
mn(t) denotes the reflection phase of the intelligent reflective surface IRS at time t,
the time-varying Doppler frequency shift of multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after passing through an (m, n) -th intelligent reflection unit is represented; n is a radical of
lRepresenting scatterers
The number of taps, l represents the number of taps,
expressing (m, n) -th intelligent reflection unit and scatterer
The time-varying distance between them,
representing scatterers
And the time-varying distance between the receiving-end antenna unit q,
representing multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q via an intelligent reflecting surface IRS and a scatterer
The latter time-varying doppler shift.
Further, the optimization problem in step S3 is expressed as follows:
where t represents a time variable, θ
mn(t) denotes the reflection phase of the intelligent reflective surface IRS at time t,
represents a statistical mean operation, h
pq(t) represents the complex channel gain of the multipath component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q.
Further, step S4 includes the following steps:
s41, since the power of the received signal is mainly concentrated on the direct component reflected by the intelligent reflective surface IRS, the optimization problem in step S3 is simplified as follows:
wherein, t represents a time variable,
the method comprises the steps that complex channel gain of a direct component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after being reflected by an intelligent reflecting surface IRS is represented, and | represents an absolute value function;
s42, handle
In (1)The phase relationship is substituted into the formula (5), and the optimization problem is further simplified as follows:
wherein ,
the method is characterized in that the time-varying Doppler frequency shift of multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after passing through an (M, N) -th intelligent reflection unit is represented, M represents the number of row reflection units of an intelligent reflection surface, M represents the row position index of the intelligent reflection unit, N represents the number of column reflection units of the intelligent reflection surface, N represents the column position index of the intelligent reflection unit, lambda represents the carrier wavelength, pi represents the circumferential ratio, xi represents the circumferential ratio
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, theta
mn(t) represents the reflection phase of the intelligent reflective surface IRS at time t.
Further, the step of solving the optimal IRS reflection phase of the intelligent reflection surface in step S5 includes the following steps:
s51, solving the optimal IRS reflection phase of the intelligent reflection surface according to the optimal IRS reflection phase optimization problem obtained by the formula (6)
The expression of (a) is as follows:
wherein ,
to represent
And 2 pi two numberThe remainder of the division;
s52, because formula (7) does not consider the time-varying Doppler frequency shift of the multipath component between the UAV antenna unit p and the receiving end antenna unit q after passing through the (m, n) -th intelligent reflection unit
Therefore, will
Further rewritten as:
wherein ,
the time-varying Doppler frequency shift of multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after passing through an (m, n) -th intelligent reflection unit is shown.
Further, step S6 includes the following steps:
s61, solving a space-time correlation function based on the intelligent reflecting surface IRS assistance according to the definition by using the complex channel gain obtained in the step S2 and the optimal reflecting phase of the intelligent reflecting surface IRS obtained in the step S5, and calculating the formula as follows:
firstly, according to the definition formula of the space-time correlation function:
wherein ,
representing a spatio-temporal correlation function, δ, between two time-varying transfer functions
TRepresenting the antenna spacing, δ, between the unmanned aerial vehicle UAV antenna units
RRepresents the antenna spacing between the antenna units at the user end, tau represents the propagation delay, t represents the time variable,
representing statistical mean operations (·)
*Denotes a complex conjugate operation, h
pq(t) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p and the user side antenna unit q, h
p′q′(t + τ) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p 'and the user-side antenna unit q' after a time delay τ; | represents an absolute value function;
the expressions of the complex channel gain functions obtained in step S2 are respectively substituted to obtain specific spatio-temporal correlation functions as follows:
wherein ,
represents the space-time correlation of the direct component between the UAV antenna unit and the receiving end antenna unit, wherein lambda represents the carrier wavelength, pi represents the circumferential ratio, xi
pq(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna unit p and the user side antenna unit q
p′q′(t + τ) represents the time-varying distance between the human-machine UAV antenna unit p 'and the user-side antenna unit q' after a time delay τ,
indicating unmanned aerial vehicleTime-varying doppler shift of the direct component between UAV antenna unit p and user side antenna unit q;
represents the space-time correlation of the scattering component between the UAV antenna unit and the receiving-end antenna unit, N
lRepresenting scatterers
The number of taps, l, xi
pnl(t) denotes unmanned aerial vehicle UAV antenna unit p and scatterers
Time-varying distance between, xi
p′nl(t + τ) denotes the unmanned aerial vehicle UAV antenna unit p' and scatterers after a time delay τ
The time-varying distance between them,
representing scatterers
And the time-varying distance between the subscriber-side antenna unit q,
showing scatterers after a time delay of tau
And the time varying distance between the subscriber side antenna unit q',
a time-varying doppler frequency shift representing a scattering component between the unmanned aerial vehicle UAV antenna unit p and the user-side antenna unit q;
representation of Unmanned Aerial Vehicle (UAV) skyThe space-time correlation of the direct component reflected by the intelligent reflector IRS between the line unit and the receiving end antenna unit, wherein M represents the number of row reflecting units of the intelligent reflector IRS, N represents the number of column reflecting units of the intelligent reflector IRS, and xi
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
p′m′n′(t + τ) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p ' and the (m ', n ') -th smart reflector unit after a time delay τ has elapsed
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the user terminal antenna unit q, xi
m′n′q′(t + τ) represents the time-varying distance between the (m ', n ') -th intelligent reflection unit and the user side antenna unit q ' after the time delay τ;
indicating that the space between the UAV antenna unit and the receiving end antenna unit passes through the intelligent reflector IRS and the scatterer
Space-time dependence of the reflected scattered component, δ
NIndicating the distance, δ, between the column reflective elements of the intelligent reflective surface IRS
MRepresents the distance between the column reflective elements of the intelligent reflective surface IRS;
and S62, when the number of reflection units of the intelligent reflection IRS, the reflection phase of the intelligent reflection surface IRS and the flight track of the unmanned aerial vehicle UAV are changed, analyzing the influence of the parameter changes on the channel characteristics of the unmanned aerial vehicle through the correlation according to the space-time correlation function obtained above.
In one embodiment of the invention: an apparatus comprising a memory and a processor, wherein:
a memory for storing a computer program capable of running on the processor;
and the processor is used for executing the steps of the unmanned aerial vehicle channel model building method based on the large-scale intelligent reflection unit when the computer program is run.
In another embodiment of the invention: a storage medium having stored thereon a computer program that, when executed by at least one processor, performs the steps of the above-described method for channel model building for a large-scale intelligent reflection unit-based drone.
Has the advantages that: compared with the prior art, the unmanned aerial vehicle channel model building method based on the large-scale intelligent reflection unit adopts the intelligent reflection surface IRS to assist the unmanned aerial vehicle communication compared with the traditional unmanned aerial vehicle communication technology, and proves that the signal power of the received signal can be improved by adopting the intelligent reflection surface IRS. The reflection phase of the intelligent reflection surface IRS is the key to increase the received signal power. Therefore, the main contribution of the invention is to provide the optimal reflection phase of the intelligent reflecting surface, and the power of the received signal can be maximized. Meanwhile, the invention also solves the space-time correlation function based on the assistance of the intelligent reflector IRS, and researches the influence of the parameter changes on the channel characteristics of the unmanned aerial vehicle by analyzing the changes of the space-time correlation function when the number of reflection units of the intelligent reflector IRS, the reflection phase of the intelligent reflector IRS and the flight track of the unmanned aerial vehicle UAV are changed. Therefore, the method provides certain help for exploring the influence of the intelligent reflector IRS on the channel statistical characteristics of the unmanned aerial vehicle.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
According to the invention, the intelligent reflector IRS is adopted to control the propagation environment of the UAV channel, the near field effect of the large intelligent reflector IRS is considered, the second-order approximation of spherical wavefront is utilized for modeling, the optimal received signal power is considered, and the optimal intelligent reflector IRS reflection phase is obtained. The invention considers the capability of the intelligent reflector IRS to change the unmanned aerial vehicle channel propagation environment, namely the influence of the number of the intelligent reflector IRS reflection units and the size of the reflection units on the unmanned aerial vehicle channel statistical characteristics. The invention considers the influence of different flight trajectories of the unmanned aerial vehicle UAV on the channel statistical characteristics of the unmanned aerial vehicle. The invention considers the influence of the intelligent reflector IRS on the channel statistical characteristics of the unmanned aerial vehicle under different time-varying reflection phases. The invention considers the unmanned aerial vehicle channel based on the intelligent reflector IRS, explores the influence of the intelligent reflector IRS on the unmanned aerial vehicle channel statistical characteristics, and better provides a basis for system performance analysis and precoding algorithm design in the future.
The invention discloses an unmanned aerial vehicle channel model building method based on a large-scale intelligent reflection unit, which comprises the following steps:
s1, according to an actual ground communication scene of the unmanned aerial vehicle UAV, determining the position relation among the unmanned aerial vehicle UAV, the intelligent reflector IRS and the receiving end, establishing an unmanned aerial vehicle channel model based on a large-scale intelligent reflector unit, and obtaining the complex channel gain of a channel;
the established unmanned aerial vehicle channel model based on the large-scale intelligent reflection unit is used for serving all users in the cell, and the intelligent reflection surface is configured on the surface of a building at the edge of the served cell. Then using a three-dimensional cylinder toAnd simulating a vertical building structure around the receiving end, and simulating the near-field effect of the large-scale intelligent reflecting surface by utilizing second-order approximation of spherical wavefront. The scatterers are assumed to be positioned on the surfaces of the three-dimensional cylinders, and the intelligent reflecting units on the intelligent reflecting surface are uniformly arranged; considering the change of the intelligent reflector IRS to the channel propagation environment, the established large-scale intelligent reflector-unit-based unmanned aerial vehicle channel model is shown in fig. 1, the invention adopts a three-dimensional cylinder to simulate the intelligent reflector IRS, the unmanned aerial vehicle UAV, and scatterers around the receiving end, the intelligent reflector adopts a uniform planar reflection array unit, the number of reflection units in each row is assumed to be M, the number of reflection units in each column is assumed to be N, and the intelligent reflector IRS is assumed to be configured on the surface of the building, so that all users in the local cell can be served. The height of the Unmanned Aerial Vehicle (UAV) is obviously higher than that of a ground building, no building is shielded between the UAV and the intelligent reflecting surface IRS, and a direct link is assumed between the UAV and the intelligent reflecting surface IRS. Wherein, in FIG. 1, H
IRSIndicating the height of the intelligent reflecting surface, H
TIndicating the altitude, R, of the unmanned aerial vehicle UAV
lRadius, xi, of the ith three-dimensional cylinder
pq(t) represents the time-varying distance between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q,
representing unmanned aerial vehicle UAV antenna unit p and scatterers
The time-varying distance between them,
representing scatterers
And time-varying distance xi between the antenna unit q and the receiving end
TRRepresents the distance, ξ, between the unmanned aerial vehicle UAV and the receiving end
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) denotes (m, n) -th intelligenceThe time-varying distance between the reflecting unit and the receiving-end antenna unit q,
expressing (m, n) -th intelligent reflection unit and scatterer
Time-varying distance between, xi
IRSRIndicating the distance, ξ, between the intelligent reflector IRS and the receiving end
IRSTRepresenting the distance, θ, between the unmanned aerial vehicle UAV and the Intelligent Reflector IRS
IRSIndicating the direction of the intelligent reflecting surface IRS in the x-y plane, alpha
IRSTRepresents the direction, θ, of the unmanned aerial vehicle UAV relative to the Intelligent reflective surface IRS in the x-y plane
RRepresenting the direction of the receiving end in the x-y plane, alpha
TRRepresenting the orientation of the unmanned UAV relative to the receiver in the x-y plane, v
TRepresenting the speed of movement, gamma, of the unmanned aerial vehicle UAV
TAzimuth, theta, representing unmanned aerial vehicle UAV motion, elevation, representing unmanned aerial vehicle UAV motion
TRepresenting the direction, v, of the receiving end of the unmanned UAV in the x-y plane
RIndicating the speed of movement, gamma, of the receiving end
RIndicating the azimuth angle of the receiving end,
representing scatterers
Is reaching an azimuth angle (AAoA),
representing scatterers
Is in the direction of departure (AAoD),
representing scatterers
Is the angle of elevation of arrival (EAoA),
representing scatterers
Elevation of departure (EAoD).
Obtaining the complex channel gain h of the channel according to the unmanned aerial vehicle channel model based on the large-scale intelligent reflection unitpqThe formula for the calculation of (t, τ) is:
where t represents a time variable, L represents the number of taps, L represents the total number of taps, c
lDenotes the gain of the first tap, h
l,pq(t) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q directly without the smart reflector IRS,
represents the complex channel gain, τ, between the UAV antenna unit p and the receiving-end antenna unit q through the intelligent reflector IRS
l(t) represents the propagation delay of the ith tap, and δ (·) represents the impulse function.
S2, obtaining the complex channel gain according to the unmanned aerial vehicle channel model based on the large-scale intelligent reflection unit, wherein the complex channel gain comprises two components: the complex channel gain of the Unmanned Aerial Vehicle (UAV) which is directly transmitted with the receiving end without the Intelligent Reflector (IRS) and the complex channel gain of the UAV which is transmitted with the receiving end by the Intelligent Reflector (IRS);
s21, complex channel gain h between Unmanned Aerial Vehicle (UAV) antenna unit p and receiving end antenna unit q directly without Intelligent Reflector (IRS)l,pq(t) is represented as follows:
wherein ,
representing the complex channel gain of the direct component between the unmanned aerial vehicle UAV antenna element p and the receiving end antenna element q,
a complex channel gain representing a scattering component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q;
wherein ,G
tDenotes the transmit antenna gain, G
rRepresenting the receiving-end antenna gain, gamma
TRRepresents the path loss, K, of the UAV to the receiving end
1Representing the rice factor, λ the carrier wavelength, t the time variable, π the circumference ratio, ξ
pq(t) represents the time-varying distance between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q,
the time-varying Doppler frequency shift of a direct component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q is shown, delta (l-1) represents a delay impact function after l taps, and N
lRepresenting scatterers
The number of the (c) component (a),
representing unmanned aerial vehicle UAV antenna unit p and scatterers
The time-varying distance between them,
representing scatterers
And the time-varying distance between the receiving-end antenna unit q,
representing the time-varying doppler shift of the scattered component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q.
S22, complex channel gain of unmanned aerial vehicle UAV through transmission between intelligent reflector IRS and receiving end
Is represented as follows:
wherein ,
represents the complex channel gain of the direct component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q after being reflected by the intelligent reflecting surface IRS,
the complex channel gain of a scattering component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after being scattered by an intelligent reflecting surface IRS and a scatterer is represented;
wherein m represents the row position index of the intelligent reflection unit, and n represents the column bit of the intelligent reflection unitIndex, M represents the number of row reflective elements of the intelligent reflective surface, N represents the number of column reflective elements of the intelligent reflective surface, G
tRepresenting the gain of the transmitting antenna, G representing the gain of the IRS reflecting element, G
rRepresenting the receiving-end antenna gain, gamma
TIRRepresents the path loss, K, from the UAV to the IRS and then to the receiving end
2Representing the rice factor, λ the carrier wavelength, t the time variable, π the circumference ratio, ξ
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, theta
mn(t) denotes the reflection phase of the intelligent reflective surface IRS at time t,
represents the time-varying Doppler shift between the multipath components of the UAV antenna unit p and the receiving end antenna unit q after passing through the (m, N) -th intelligent reflection unit, N
lRepresenting scatterers
The number of taps, l represents the number of taps,
expressing (m, n) -th intelligent reflection unit and scatterer
The time-varying distance between them,
representing scatterers
And the time-varying distance between the receiving-end antenna unit q,
representing multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q via an intelligent reflecting surface IRS and a scatterer
The latter time-varying doppler shift.
S3, designing an optimization problem according to the maximum power of the received signal according to the unmanned aerial vehicle channel model based on the large-scale intelligent reflection unit, wherein the design purpose is to maximize the power of the received signal;
thus, the optimization problem is represented as follows:
where t represents a time variable, θ
mn(t) denotes the reflection phase of the intelligent reflective surface IRS at time t,
represents a statistical mean operation, h
pq(t) represents the complex channel gain of the multipath component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q.
S4, simplifying optimization problems: the received signal is composed of a direct component, a scattered component, a direct component via the intelligent reflecting surface IRS and a scattered component via the intelligent reflecting surface, which results in a computationally significant complexity of the optimization problem proposed in step S3, so that further simplification of the problem is required in order to reduce the complexity. The invention discovers that when the scale of the intelligent reflection unit in the intelligent reflection surface IRS is larger, the power of the received signal is mainly controlled by the reflection signal passing through the intelligent reflection surface IRS, and the complex channel gain of the reflection signal is mainly determined by the direct component, so that the process of solving the reflection phase can be simplified;
the simplifying process comprises the following steps:
s41, since the power of the received signal is mainly concentrated on the direct component reflected by the intelligent reflective surface IRS, the optimization problem in step S3 is simplified as follows:
wherein, t represents a time variable,
the method comprises the steps that complex channel gain of a direct component between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after being reflected by an intelligent reflecting surface IRS is represented, and | represents an absolute value function;
s42, handle
The phase relation in (3) is substituted into the formula (5), and the optimization problem is further simplified into:
wherein ,
the time-varying Doppler frequency shift of multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after passing through an (M, N) -th intelligent reflection unit is represented, M represents the number of row reflection units of an intelligent reflection surface, M represents the row position index of the intelligent reflection unit, N represents the number of column reflection units of the intelligent reflection surface, N represents the column position index of the intelligent reflection unit, lambda represents the carrier wavelength, and xi
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, theta
mn(t) represents the reflection phase of the intelligent reflective surface IRS at time t.
And S5, solving the optimal IRS reflection phase of the intelligent reflection surface according to the simplified optimization problem. It should be noted that when solving the reflected phase, the doppler shift of the direct component needs to be subtracted, and only the doppler shift is removed to enhance the received signal power, and the existence of the doppler shift is often ignored in the previous research.
Which comprises the following steps:
S51、according to the optimization problem of the optimal IRS reflection phase of the intelligent reflecting surface obtained by the formula (6), the optimal IRS reflection phase of the intelligent reflecting surface can be solved
The expression of (a) is as follows:
where t represents a time variable, λ represents the carrier wavelength, ξpmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unitmnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, pi represents the circumferential ratio, thetamn(t) denotes the reflection phase of the intelligent reflective surface IRS at time t, and mod (a, b) denotes the remainder of the division of the two numbers a and b.
S52, because formula (7) does not consider the time-varying Doppler frequency shift of the multipath component between the UAV antenna unit p and the receiving end antenna unit q after passing through the (m, n) -th intelligent reflection unit
Therefore, the invention will
Further rewritten as:
where t represents a time variable, λ represents the carrier wavelength, ξ
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
mnq(t) represents the time-varying distance between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q, pi represents the circumferential ratio, theta
mn(t) denotes the reflection bit of the intelligent reflective surface IRS at time t,
the time-varying Doppler frequency shift of multipath components between an Unmanned Aerial Vehicle (UAV) antenna unit p and a receiving end antenna unit q after passing through an (m, n) -th intelligent reflection unit is shown, and mod (a, b) shows a remainder of division of a number and a number b.
S6, solving a space-time correlation function based on the assistance of the intelligent reflecting surface IRS according to the definition by the complex channel gain obtained in the step S2 and the optimal reflecting phase of the intelligent reflecting surface IRS obtained in the step S5, and determining the influence of different parameters on the channel characteristics of the unmanned aerial vehicle through correlation analysis.
Which comprises the following steps:
s61, firstly, according to the definition formula of the space-time correlation function:
wherein ,
representing a spatio-temporal correlation function, δ, between two time-varying transfer functions
TRepresenting the antenna spacing, δ, between the unmanned aerial vehicle UAV antenna units
RDenotes the antenna spacing between antenna elements at the receiving end, τ denotes the propagation delay, t denotes a time variable,
representing statistical mean operations (·)
*Denotes a complex conjugate operation, h
pq(t) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q, h
p′q′(t + τ) represents the complex channel gain between the unmanned aerial vehicle UAV antenna unit p 'and the receiving end antenna unit q' after a time delay τ; | represents an absolute value function;
the expressions of the complex channel gain functions obtained in step S2 are respectively substituted to obtain specific spatio-temporal correlation functions as follows:
wherein ,
represents the space-time correlation of the direct component between the UAV antenna unit and the receiving end antenna unit, wherein lambda represents the carrier wavelength, pi represents the circumferential ratio, xi
pq(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q
p′q′(t + τ) represents the time-varying distance between the unmanned aerial vehicle UAV antenna unit p 'and the receiving end antenna unit q' after a time delay τ,
a time-varying doppler frequency shift representing a direct component between an unmanned aerial vehicle UAV antenna unit p and a receiving end antenna unit q;
represents the space-time correlation of the scattering component between the UAV antenna unit and the receiving-end antenna unit, N
lRepresenting scatterers
The number of taps, l, xi
pnl(t) denotes unmanned aerial vehicle UAV antenna unit p and scatterers
Time-varying distance between, xi
p′nl(t + τ) denotes the unmanned aerial vehicle UAV antenna unit p' and scatterers after a time delay τ
The time-varying distance between them,
representing scatterers
And the time-varying distance between the receiving-end antenna unit q,
showing scatterers after a time delay of tau
And the time-varying distance between the receiving-end antenna unit q',
a time-varying doppler frequency shift representing a scattering component between the unmanned aerial vehicle UAV antenna unit p and the receiving end antenna unit q;
the space-time correlation of the direct component reflected by the intelligent reflecting surface IRS between the unmanned aerial vehicle UAV antenna unit and the receiving end antenna unit is represented, M represents the number of row reflecting units of the intelligent reflecting surface IRS, N represents the number of column reflecting units of the intelligent reflecting surface IRS, and xi represents
pmn(t) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p and the (m, n) -th smart reflector unit
p′m′n′(t + τ) represents the time-varying distance, ξ, between the unmanned aerial vehicle UAV antenna units p ' and the (m ', n ') -th smart reflector unit after a time delay τ has elapsed
mnq(t) represents the time-varying distance, ξ, between the (m, n) -th intelligent reflection unit and the receiving-end antenna unit q
m′n′q′(t + τ) represents a time-varying distance between the (m ', n ') -th intelligent reflection unit and the receiving-end antenna unit q ' after a time delay τ;
indicating that the space between the UAV antenna unit and the receiving end antenna unit passes through the intelligent reflector IRS and the scatterer
Space-time dependence of the reflected scattered component, δ
NIndicating the distance, δ, between the column reflective elements of the intelligent reflective surface IRS
MThe distance between the column reflection elements of the intelligent reflection surface IRS is indicated and τ represents the propagation delay.
And S62, when the number of reflection units of the intelligent reflection IRS, the reflection phase of the intelligent reflection surface IRS and the flight track of the unmanned aerial vehicle UAV are changed, analyzing the influence of the parameter changes on the channel characteristics of the unmanned aerial vehicle through the correlation according to the space-time correlation function obtained above.
The invention also provides an apparatus comprising a memory and a processor, wherein:
a memory for storing a computer program capable of running on the processor;
and the processor is used for executing the steps of the unmanned aerial vehicle channel model building method based on the large-scale intelligent reflection unit when the computer program is run.
The invention also provides a storage medium, wherein the storage medium is stored with a computer program, and the computer program is executed by at least one processor to realize the steps of the unmanned aerial vehicle channel model building method based on the large-scale intelligent reflection unit.
Fig. 2 is a schematic diagram of three different flight trajectories of an unmanned aerial vehicle UAV. In fig. 2, the drone under path I is far from the intelligent reflector IRS and the receiving end, parameters
and γ
TSet to 0 and 90, respectively, the drone under path II is in 0 to 5.2 seconds
and γ
TMoving towards the receiving end in the-30 ° direction, the drone under path III is in 0 to 5.2 seconds
and γ
TMoves towards the intelligent reflection IRS reflection unit in the direction of-120 °. Wherein,
elevation angle, gamma, representing the direction of movement of an Unmanned Aerial Vehicle (UAV)
TRepresenting the unmanned aerial vehicle UAV direction of motion azimuth.
Fig. 3 is a comparison graph of absolute envelope amplitudes of a traditional unmanned aerial vehicle channel model and an unmanned aerial vehicle channel model based on the assistance of an intelligent reflector IRS in different unmanned aerial vehicle flight trajectories. As can be seen from fig. 3, as the unmanned aerial vehicle UAV moves towards the intelligent reflective surface IRS, i.e., path III, the envelope size received by the unmanned aerial vehicle model employing the intelligent reflective surface IRS gradually increases. Conversely, when the unmanned aerial vehicle UAV is far away from the intelligent reflective surface IRS, the size of the envelope of the unmanned aerial vehicle model using the intelligent reflective surface IRS gradually decreases, i.e., the path I and the path II. And when the intelligent reflector IRS is not adopted, when the unmanned aerial vehicle UAV gradually moves to the receiving end, the absolute envelope amplitude of the unmanned aerial vehicle model is not obviously increased in the path II. This shows that the propagation environment between unmanned aerial vehicle and the receiving terminal can be effectively changed to intelligent plane of reflection IRS.
FIG. 4 shows a conventional UAV channel model and an intelligent reflector IRS-based broadband UAV channel model in different IRS reflection phases θmn(t) absolute envelope size comparison of the broadband UAV-MIMO model. In phase 1, the IRS reflection phase is aligned with the time-varying phase and the doppler shift of the LoS component, and in phase 2, the IRS reflection phase, i.e., the optimal reflection phase of the intelligent reflective surface IRS, is set according to equation (8). As can be seen from fig. 6, when the size of the reflection unit of the intelligent reflector IRS is increased, the performance and the multipath fading phenomenon of the UAV-MIMO communication system employing the intelligent reflector IRS can be significantly improved. Meanwhile, when the number of the intelligent reflecting surface IRS reflecting units is greater than 100 × 100, the performance of method 2 is the same as that of method 1. The results show that the method has the advantages of high yield,the provided IRS reflection phase is suitable for research of a broadband unmanned aerial vehicle communication system based on intelligent reflector IRS assistance.
Fig. 5 is a comparison graph of the emission spatial correlation curves of the conventional unmanned aerial vehicle channel model and the unmanned aerial vehicle channel model based on the assistance of the intelligent reflector IRS in three different unmanned aerial vehicle trajectories. It can be seen from fig. 5 that the spatial correlation of the drone-MIMO model without the assistance of the intelligent reflector IRS is affected by the drone trajectory, and that the normalized antenna spacing δTThe spatial correlation gradually decreases as/λ increases. While the spatial correlation of the unmanned aerial vehicle channel model based on the assistance of the intelligent reflector IRS is constant. This also means that the intelligent reflective surface IRS can reduce the non-stationarity of the transmitting antenna elements in the spatial domain.
Fig. 6 is a comparison graph of emission spatial correlation curves of a traditional unmanned aerial vehicle channel model and a broadband unmanned aerial vehicle channel model based on the assistance of an intelligent reflector IRS under three different unmanned aerial vehicle flight trajectories. As can be seen in FIG. 6, the spatial correlation of the broadband drone channel model without the assistance of the intelligent reflector IRS is subject to the drone trajectory and the normalized antenna spacing δTThe effect of/λ, which is consistent with the results obtained in fig. 5.
To sum up, the unmanned aerial vehicle channel model building method based on the large-scale intelligent reflection unit is a non-stable three-dimensional broadband channel model building method of an unmanned aerial vehicle multi-input multi-output communication system based on a large-scale intelligent reflection surface, and comprises the steps of intelligent reflection surface time-varying reflection phase design: the received signal power maximization is taken as a target design optimization problem, and the optimization problem is solved to obtain an optimal time-varying reflection phase; designing time-varying Doppler frequency shift parameters: obtaining time-varying Doppler frequency shift parameters among the unmanned aerial vehicle, the receiving end and the intelligent reflecting surface according to the intelligent reflecting surface assisted geometric model; analyzing the statistical characteristics of the channel: and analyzing the statistical characteristics of the unmanned aerial vehicle MIMO channel model based on the assistance of the intelligent reflector according to the time-varying reflection phase and the time-varying parameters of the intelligent reflector. In the invention, the communication system adopting the intelligent reflecting surface IRS has better signal propagation environment, can obviously improve the power of received signals and reduce the influence of multipath fading and Doppler frequency shift on the received signals. Therefore, the model building method can provide powerful support for the exploration of the key technology of the 6G communication system.