Air-ground channel dynamic modeling method and system for high-mobility unmanned aerial vehicle

Abstract

The invention is suitable for the technical field of wireless communication, and provides an air-ground channel dynamic modeling method and system for a high-mobility unmanned aerial vehicle, and the method comprises the following steps: constructing an unmanned aerial vehicle dynamic model, and quantifying the antenna position offset caused by the attitude change of the unmanned aerial vehicle; based on topographic data and obstacle data, constructing a three-dimensional digital scene model for restoring a real scene, simulating multipath propagation by adopting a ray tracing method, and extracting channel statistical characteristic parameters; constructing a statistical model of the jitter characteristics of the unmanned aerial vehicle, and quantifying the position offset of the unmanned aerial vehicle under the jitter effect; and according to the unmanned aerial vehicle kinetic model, the channel statistical characteristic parameters and the statistical model of the unmanned aerial vehicle jitter characteristics, constructing a channel model fusing the attitude dynamics and the jitter effect by using an analytical method. The method solves the problem of how to fully consider attitude dynamics and jitter effects in unmanned aerial vehicle communication to accurately model channel characteristics and quantify key parameters, and has high applicability to complex and changeable unmanned aerial vehicle communication scenes.

Description

Technical Field

[0001] This invention belongs to the field of communication technology, and in particular relates to a method and system for dynamic modeling of air-to-ground channels for highly maneuverable unmanned aerial vehicles (UAVs). Background Technology

[0002] Currently, the widely deployed 5G technologies are insufficient to meet the demands of emerging services for existing communication performance. Next-generation mobile communication systems will extend to natural spaces such as space, sky, and land, achieving ubiquitous connectivity with full-domain intelligence. Unmanned aerial vehicles (UAVs), as small aircraft, play a crucial role in achieving global, three-dimensional, and deep 6G air-space coverage. The UAV communication network is a vital component of the UAV system; it is the neural network enabling UAVs to communicate externally, facilitating information exchange between UAVs or between UAVs and other receiving stations, thus meeting the information transmission needs of the UAV platform. The UAV air-to-ground channel model reflects the propagation patterns of wireless signals in specific communication scenarios and serves as an important reference for UAV communication network planning and system design. As UAV communication scenarios become increasingly widespread, higher demands are placed on the universality, accuracy, and complexity of UAV channel models.

[0003] To efficiently and reliably design and evaluate the performance of UAV wireless communication systems, it is crucial to accurately characterize the wireless propagation channel using appropriate methods, thereby establishing a reliable, accurate, and easy-to-use UAV communication channel model. When designing, simulating, and optimizing wireless communication systems, a mathematical description of the wireless propagation channel is used, and the complex and ever-changing channel is abstracted to obtain a channel model. However, compared to traditional terrestrial communication channels such as cellular and vehicle-mounted communications, UAV air-to-ground wireless communication channels have their unique characteristics, mainly reflected in:

[0004] 1. Significantly high dynamics: The high-speed movement, attitude angle changes, and random jitter of UAVs cause the distance, relative speed, and propagation path between the transmitter and receiver to change rapidly over time, resulting in strong time-varying characteristics of channel parameters, which poses a significant challenge to channel stability.

[0005] 2. Complex propagation environment: UAVs fly at extremely wide altitudes, covering both low and high altitudes. This means their propagation paths may simultaneously include line-of-sight (LoS), non-line-of-sight (NLoS), and reflection and scattering paths from ground structures, vegetation, water surfaces, etc. Furthermore, the proportion and attenuation characteristics of these paths vary significantly in different flight environments, increasing the complexity of channel modeling.

[0006] 3. Prominent Three-Dimensional Spatial Characteristics: Unlike terrestrial communication, which primarily relies on two-dimensional planar propagation, UAV air-to-ground channels involve three-dimensional spatial issues. Furthermore, changes in UAV attitude can alter the directivity of antenna gain. Therefore, the channel fading characteristics in the three-dimensional angular domain must be fully considered during modeling to accurately characterize the channel.

[0007] Currently, existing models simplify UAVs into point masses for analysis, focusing only on the position and trajectory of their center of mass, ignoring the direct impact of UAV attitude changes and random jitter on the channel. This results in insufficient model accuracy, making it difficult to meet the design and optimization requirements of highly maneuverable UAV communication systems. Summary of the Invention

[0008] The purpose of this invention is to provide a dynamic modeling method for air-to-ground communication channels for highly maneuverable unmanned aerial vehicles (UAVs), aiming to solve the aforementioned technical problems.

[0009] This invention is implemented as follows: a dynamic modeling method for air-to-ground communication channels for highly maneuverable unmanned aerial vehicles (UAVs), comprising the following steps:

[0010] Construct a dynamic model of the UAV and quantify the antenna position offset caused by changes in UAV attitude;

[0011] Based on the collected terrain and obstacle data, a three-dimensional digital scene model that recreates the real scene is constructed.

[0012] Based on a three-dimensional digital scene model, ray tracing is used to simulate multipath propagation and extract channel statistical characteristic parameters.

[0013] Construct a statistical model of the drone's jitter characteristics and quantify the drone's positional shift under jitter.

[0014] Based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, an analytical method is used to construct a channel model that integrates attitude dynamics and jitter effects.

[0015] Furthermore, the steps for constructing a UAV dynamics model and quantifying the antenna position offset caused by changes in UAV attitude include:

[0016] The basic physical parameters of the UAV are set as physical constraints. The physical constraints are input and combined with real-time IMU data. The discretized Euler integral method is used to establish the dynamic equations. The rotation matrix is ​​used to realize the force and torque conversion between the machine system and the inertial frame. The dynamic equations are iteratively updated to complete the complete mapping from control input to global motion trajectory. This realizes the UAV dynamic modeling that combines hard constraints of physical parameters with soft feedback of real-time state, and obtains the UAV dynamic model.

[0017] Based on the attitude data output by the UAV dynamics model, a position offset model caused by the UAV attitude is established. The position offset caused by the attitude angle is converted into pitch and azimuth angles in the ground coordinate system through the attitude matrix, and the antenna position offset is quantified.

[0018] Furthermore, the dynamic equation is:

[0019] ;

[0020] In the formula, This represents the angular velocity of the UAV at time t in the body coordinate system. Let J represent the roll, pitch, and yaw angles of the body coordinate system relative to the inertial coordinate system at time t, where J is the moment of inertia matrix and W is the angular velocity transformation matrix. This represents the velocity of the UAV at time t in the ground coordinate system. This represents the position of the UAV at time t in the ground coordinate system; and These represent the control torque and total tension in the body coordinate system, respectively. The matrix represents the rotation of the body coordinate system, and m and g represent the mass and gravitational acceleration of the UAV.

[0021] The position offset model caused by the drone's attitude is as follows:

[0022] ;

[0023] In the formula, This represents the position offset caused by the drone's attitude. and These are the azimuth and elevation angles, respectively. and This represents the displacement caused by the swaying at that angle; the models for the azimuth and pitch angles are as follows:

[0024] ;

[0025] In the formula, , and These are the pitch angle, roll angle, and yaw angle output by the UAV dynamics model, respectively.

[0026] Furthermore, the steps of constructing a 3D digital scene model that recreates the real scene based on the collected terrain and obstacle data specifically include:

[0027] Digital analysis of real-world scenes is performed; specifically, for elevation maps, key information is extracted through algorithms: elevation values ​​are converted into vertex heights in three-dimensional space, geographic coordinates are mapped into position parameters in a digital coordinate system, and the physical characteristics of real terrain are transformed into computable geometric data to form an initial terrain framework and obtain a terrain grid.

[0028] The parsed real-world scene data is classified, modeled, and fused to obtain obstacle models. For static obstacles, geometric models are constructed using their size and coordinate data. For dynamic obstacles, positional change parameters over time are associated to form digital volumes with temporal information. These obstacle models are aligned with the terrain grid according to their real-world spatial locations, and collision detection algorithms are used to ensure that their spatial relationships conform to physical logic, thus avoiding spatial conflicts in the digital scene.

[0029] Achieving model standardization and scene integration: Converting terrain meshes and obstacle models into the universal OBJ format, unifying their coordinate system and data structure, using vertex coordinates to record spatial positions, using texture mapping to restore surface materials, and using normal vectors to define surface orientation, ultimately achieving digital mapping of the real scene and obtaining a 3D digital scene model.

[0030] Furthermore, based on the three-dimensional digital scene model, the steps of simulating multipath propagation using ray tracing and extracting channel statistical characteristic parameters specifically include:

[0031] Configure ray emission and scene association: Set the three-dimensional coordinates of the signal source and receiver in the three-dimensional digital scene model, set the basic parameters, and define the ray propagation angle range and maximum tracking distance based on the distribution characteristics of scatterers in the scene to ensure that the ray can cover the area where all potential scatterers are located; the basic parameters include the frequency of the transmitted signal, the initial power, and the antenna type.

[0032] The ray tracing algorithm is initiated: Based on the geometric properties of the terrain mesh, obstacles, and scatterers in the 3D digital scene model, and the corresponding physical properties of their materials, the interaction point between the ray and the scatterer is located using a collision detection algorithm. Combining the scattering angle, incident power, and scatterer material characteristics, the power attenuation of each scattering path is derived using scattering theory formulas. The propagation time of the ray from the emission source to the scatterer and then to the receiving point is recorded simultaneously to determine the time delay parameters of the scattering path. The azimuth and elevation angles of the ray incident on the scatterer, as well as the scattering exit angle, are extracted through vector calculations to form a complete set of angle parameters. The geometric properties include dimensions, spatial coordinates, and surface normal vectors; the corresponding physical properties of the material include the scattering coefficient and dielectric constant.

[0033] Channel statistical characteristic parameters are aggregated as follows: First, the loss power, delay, and angle data collected from all scattering paths are systematically classified and organized. After obtaining the statistical sample set, the pitch angle, yaw angle, scatterer loss, and delay parameters in the sample set are first organized, and outliers with too low energy or invalidity are removed to ensure data validity. The mean and variance of each of the above parameters are calculated using statistical tools. For delay, the root mean square delay spread parameter is also calculated. At the same time, frequency histograms of each parameter are plotted. Based on the histograms, the probability density function is fitted using the maximum likelihood estimation method, and the cumulative distribution function is obtained by integrating the probability density function. Finally, the fitted probability density function and cumulative distribution function are compared and verified with typical channel models to determine the optimal fitting parameters, thereby obtaining the final channel statistical characteristics.

[0034] Furthermore, a statistical model of the drone's jitter characteristics is constructed, and the steps for quantifying the drone's positional shift under jitter are detailed, including:

[0035] Statistical models are used to quantify the azimuth and pitch angle shifts caused by UAV jitter, and then the position change caused by the UAV is calculated based on geometric relationships. UAV jitter includes internal jitter caused by rotor vibration and external jitter caused by external factors. The position shift of the UAV under jitter is modeled as follows:

[0036] ;

[0037] In the formula, The positional offset of the drone under jitter. This refers to the positional offset of the drone due to internal jitter. This represents the positional offset of the drone under external jitter.

[0038] Furthermore, regarding internal vibration, based on the structural parameters and dynamic characteristics of the UAV, and combined with historical vibration data, the periodic variation law of the azimuth and pitch angles caused by the rotor mechanical motion is described by the amplitude, frequency, and phase of a sine function. The specific model is as follows:

[0039] ;

[0040] In the formula, and These are the azimuth and pitch angles generated by the drone's shaking due to internal jitter. and This is the displacement caused by the swaying at that angle;

[0041] ;

[0042] In the formula, It is the internal jitter frequency.

[0043] Furthermore, regarding external jitter, influencing factors include wind speed fluctuations and airflow disturbances. These factors exhibit a random distribution, while the real-time compensation mechanism of the UAV flight control system ensures that the changes in azimuth and pitch angles statistically follow a normal distribution with zero mean and stable variance. The specific model is as follows:

[0044] ;

[0045] In the formula, and These are the azimuth and pitch angles generated by the drone's shaking under external vibrations. and This represents the displacement caused by the shaking at that angle.

[0046] Furthermore, based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, the steps for constructing a channel model that integrates attitude dynamics and jitter effects using analytical methods specifically include:

[0047] The analytical expression for the time-varying characteristics of the channel is derived using analytical methods. The channel impulse response formula is as follows:

[0048] ;

[0049] In the formula, H LoS and H NLoS These are the impulse response matrices for the direct and indirect beam paths, respectively. The transmitter uses a P-row, Q-column two-dimensional antenna array containing P*Q transmitting antennas; the receiver uses an A-row, B-column two-dimensional antenna array containing A*B receiving antennas. Therefore, each matrix contains P*Q*A*B elements, as detailed below:

[0050] ;

[0051] ;

[0052] In the formula, and Let be the channel impulse responses of the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna, respectively, expressed by the following formulas:

[0053] ;

[0054] ;

[0055] In the formula, and , respectively, are the antenna gains along the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna; and Let be the direct and indirect large-scale fading from the qp-th transmitting antenna to the ab-th receiving antenna, respectively. and These are the signal phase change factors for the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna, respectively.

[0056] Next, the parameters constituting the channel impulse response are described separately. In the UAV air-to-ground communication system, the modeling of antenna gain needs to comprehensively consider the influence of platform dynamic characteristics and hardware losses. The effective aperture of the antenna is the equivalent area that the antenna can actually receive or transmit electromagnetic waves. Its actual effectiveness will be directionally attenuated due to UAV attitude deviation and high-frequency random jitter. When the antenna main lobe deviates from the target direction, the equivalent effective aperture can be regarded as the projection of the antenna physical aperture on the signal pointing direction, forming a gain attenuation factor related to the attitude error angle. At the same time, hardware losses correct the antenna gain through a loss factor independent of spatial pointing. The attitude deviations include roll, pitch, and yaw angle changes.

[0057] ;

[0058] ;

[0059] The total antenna gain can be expressed as the product of the ideal gain, the effective aperture efficiency correction factor caused by attitude and jitter, and the hardware loss factor. The specific model is as follows:

[0060] ;

[0061] in, Indicates the wavelength of electromagnetic waves. and These are the time delays for direct and indirect beams, respectively. and It is the hardware loss of the transmitting antenna qp and the receiving antenna ab, which follows a chi-square distribution; and The effective apertures of the transmitting and receiving antennas under direct illumination conditions are given by the following formula:

[0062] ;

[0063] ;

[0064] In the formula, It is the position of the transmitter at time t. It is the position of the receiver at time t. Indicates the spacing between the antenna elements at the transmitting end. Indicates the spacing between antenna elements at the receiving end;

[0065] Immediately afterwards, and The effective apertures of the transmitting and receiving antennas under non-direct illumination conditions are given by the following formula:

[0066] ;

[0067] in, The physical area of ​​the transmitting QP antenna. This represents the mutual coupling matrix of the transmitting antenna. The physical area of ​​the receiver's ab antenna. This represents the mutual coupling matrix of the receiving antenna. and They are time t and t+ respectively. Consider the non-direct path distance of antenna qp-ab after jitter offset;

[0068] ;

[0069] ;

[0070] in and Representing t and scattering body at time Location;

[0071] Large-scale fading is used to characterize the slowly varying characteristics of signal strength in wireless communication over a large spatial and environmental range, including path loss and shadow fading (SH). The formulas for direct-path large-scale fading and indirect-path large-scale fading are as follows:

[0072] ;

[0073] ;

[0074] Path loss originates from the signal propagation distance and the propagation environment, causing energy to attenuate with distance or complex scenarios. It determines the basic signal coverage and received strength. Quantized attitude and jitter offsets are incorporated into the path loss calculation, specifically the path loss of the direct path. Path loss of non-direct path The formula is expressed as follows:

[0075] ;

[0076] ;

[0077] Where c=2 is the free space loss model; c=3~5 is the urban environment; and These are the direct and indirect trajectory lengths after incorporating quantized attitude and jitter parameters, respectively, expressed by the following formulas:

[0078] ;

[0079] ;

[0080] In the formula, This represents the position offset caused by the drone's attitude.

[0081] When radio waves encounter obstacles in their propagation path, a shadow is formed, causing a slow change in the median signal strength. This phenomenon is known as the shadowing effect. Statistically, it follows a log-normal distribution, and its probability density function can be expressed as:

[0082] ;

[0083] in, and These are the mean and variance of shadow fading, respectively.

[0084] Next, we model the small-scale parameters of the channel. The small-scale parameters caused by jitter and attitude are represented by the position offsets modeled in the above statistical model. Here, we only describe the influence caused by velocity. Based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, the small-scale parameter model of the UAV is obtained as follows:

[0085] ;

[0086] ;

[0087] The Doppler frequency shift of the channel caused by the direct trajectory velocity is modeled as follows:

[0088] ;

[0089] in, and These are the velocities of the receiving and transmitting ends, respectively, and u is a unit vector representing the signal propagation direction. The specific formula is as follows:

[0090] ;

[0091] The Doppler frequency shift formula for the channel caused by non-direct path velocity is modeled as follows:

[0092] ;

[0093] ;

[0094] in, It is the speed of the transmitting end. The relative velocity between the transmitter and the scatterer It is the velocity of the scatterer n, expressed by the formula: Similarly, It is the speed of the transmitting end. The relative velocity between the transmitter and the scatterer It is the velocity of the scatterer n, expressed by the formula: , The angle between the signal propagation direction and the velocity direction of the transmitting end. It is the angle between the signal propagation direction and the reference direction of the receiving antenna array.

[0095] Another objective of this invention is to provide an air-to-ground channel dynamic modeling system for highly maneuverable unmanned aerial vehicles (UAVs), used to implement the aforementioned air-to-ground channel dynamic modeling method, comprising:

[0096] The dynamics model building module is used to build the dynamics model of the UAV and quantify the antenna position offset caused by changes in the UAV's attitude.

[0097] The scene model building module is used to build a 3D digital scene model that recreates the real scene based on the collected terrain and obstacle data.

[0098] The channel statistical characteristic parameter extraction module is used to simulate multipath propagation based on a three-dimensional digital scene model using the ray tracing method and extract channel statistical characteristic parameters.

[0099] The statistical model building module is used to build a statistical model of the drone's jitter characteristics and quantify the drone's positional offset under jitter.

[0100] The channel model construction module is used to construct a channel model that integrates attitude dynamics and jitter effects using analytical methods, based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics.

[0101] This invention provides a dynamic air-to-ground channel modeling method for highly maneuverable unmanned aerial vehicles (UAVs). This method is based on the UAV attitude-jitter coupling effect. The overall approach establishes a UAV dynamic model and a 3D scene digitization platform, extracts channel statistical characteristics using ray tracing, and statistically models UAV jitter characteristics to construct a channel model combining attitude and jitter effects. It then uses analytical methods to derive the analytical expression for the time-varying channel characteristics and quantifies the key channel parameters including the attitude-jitter coupling effect. This method solves the problem of how to fully consider attitude dynamics and jitter effects in UAV communication to accurately model channel characteristics and quantify key parameters, and it has high applicability to complex and ever-changing UAV communication scenarios. Attached Figure Description

[0102] Figure 1 This is a flowchart illustrating the air-to-ground channel dynamic modeling method for highly maneuverable unmanned aerial vehicles (UAVs) provided in an embodiment of the present invention.

[0103] Figure 2 A schematic diagram of a channel model for highly maneuverable unmanned aerial vehicles provided in an embodiment of the present invention.

[0104] Figure 3 This is a schematic diagram of the air-to-ground channel dynamic modeling system architecture for highly maneuverable unmanned aerial vehicles (UAVs) provided in an embodiment of the present invention.

[0105] Figure 4 This is a channel impulse response diagram obtained through simulation in an embodiment of the present invention. Detailed Implementation

[0106] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0107] like Figure 1 and Figure 2 As shown, in one embodiment of the present invention, a method for dynamic modeling of air-to-ground channels for highly maneuverable unmanned aerial vehicles (UAVs) is provided, which includes the following steps:

[0108] S100. Construct a dynamic model of the UAV and quantify the antenna position offset caused by changes in UAV attitude.

[0109] S200: Based on the collected terrain and obstacle data, construct a three-dimensional digital scene model that recreates the real scene;

[0110] S300: Based on a three-dimensional digital scene model, the ray tracing method is used to simulate multipath propagation and extract channel statistical characteristic parameters.

[0111] S400, construct a statistical model of the drone's shaking characteristics, and quantify the drone's positional shift under shaking effects;

[0112] S500: Based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, an analytical method is used to construct a channel model that integrates attitude dynamics and jitter effects.

[0113] In this embodiment of the invention, the channel modeling process is as follows: First, the basic physical parameters of the UAV are set and the antenna parameters are configured. A UAV dynamic model based on Newtonian mechanics is constructed, and the three-dimensional displacement caused by attitude angle changes is quantified through Euler angle decomposition and coordinate transformation. Then, ray tracing simulation is carried out for typical urban scenarios. Based on geometric optics and uniform diffraction theory, statistical characteristics such as time delay spread and angle spread of multipath components are extracted. Subsequently, the jitter displacement of the UAV is statistically modeled based on stochastic process theory. Finally, the attitude offset output by the dynamic model, the scene statistical characteristics extracted by ray tracing, and the time-varying parameters of jitter modeling are integrated. The analytical expression of the time-varying characteristics of the channel is derived using analytical methods. A channel model that integrates UAV attitude dynamics and jitter effects is constructed. Finally, the full-link mapping from physical scene to mathematical model is completed, realizing the air-to-ground channel modeling of highly maneuverable UAVs.

[0114] The air-to-ground channel modeling method for highly maneuverable UAVs proposed in this invention reduces computational overhead during the modeling process by dynamically modeling attitude changes and statistically modeling the time-varying characteristics of the channel caused by jitter effects, thus providing accurate channel model support for evaluating system transmission performance and optimizing transmission quality.

[0115] In a preferred embodiment of the present invention, the step of constructing a UAV dynamics model and quantifying the antenna position offset caused by changes in UAV attitude, i.e., step S100, specifically includes:

[0116] S110. Set the basic physical parameters of the UAV as physical constraints, input the physical constraints, combine the real-time data collected by the inertial measurement unit (IMU), establish the dynamic equation using the discretized Euler integral method, use the rotation matrix to realize the force and torque conversion between the machine system and the inertial system, iteratively update the dynamic equation, thereby completing the complete mapping from control input to global motion trajectory, realizing the UAV dynamic modeling that combines hard constraints of physical parameters with soft feedback of real-time state, and obtaining the UAV dynamic model;

[0117] S120. Based on the attitude data output by the UAV dynamics model, establish a position offset model caused by the UAV attitude. Convert the position offset caused by the attitude angle into pitch and azimuth angles in the ground coordinate system through the attitude matrix, and quantify the antenna position offset.

[0118] In step S110, the UAV dynamics model is achieved by using physical constraints such as the actual mass of the UAV, the moment of inertia matrix of the airframe, the propeller thrust coefficient and the propeller torque coefficient, as well as data such as the airframe angular velocity and initial attitude angle acquired in real time by the IMU, as input variables for UAV dynamics modeling.

[0119] Specifically, the dynamic equations used in the modeling are:

[0120] ;

[0121] In the formula, This represents the angular velocity of the UAV at time t in the body coordinate system. Let J represent the roll, pitch, and yaw angles of the body coordinate system relative to the inertial coordinate system at time t, where J is the moment of inertia matrix and W is the angular velocity transformation matrix. This represents the velocity of the UAV at time t in the ground coordinate system. This represents the position of the UAV at time t in the ground coordinate system; and These represent the control torque and total tension in the body coordinate system, respectively. The matrix represents the rotation of the body coordinate system, and m and g represent the mass and gravitational acceleration of the UAV.

[0122] In step S120, based on the above UAV dynamics model, the position offset caused by the UAV attitude angle is converted into pitch and azimuth angles in the ground coordinate system using the attitude matrix. The position offset model caused by the UAV attitude is as follows:

[0123] ;

[0124] In the formula, This represents the position offset caused by the drone's attitude. and These are the azimuth and elevation angles, respectively. and This represents the displacement caused by the swaying at that angle; the azimuth and pitch angle models are as follows:

[0125] ;

[0126] In the formula, , and These are the pitch angle, roll angle, and yaw angle output by the UAV dynamics model, respectively.

[0127] In a preferred embodiment of the present invention, the step of constructing a three-dimensional digital scene model that recreates the real scene based on the collected terrain data and obstacle data, namely step S200, specifically includes:

[0128] S210. Perform digital analysis on real-world scenes; among them, for elevation maps of typical urban scenes, extract key information through algorithms: convert elevation values ​​into vertex heights in three-dimensional space, map geographic coordinates into position parameters in a digital coordinate system, transform the physical characteristics of real terrain into computable geometric data, form an initial terrain framework, and obtain a terrain grid.

[0129] S220. Classify, model, and fuse the parsed real-world scene data to obtain obstacle models. For static obstacles, construct geometric models using their size and coordinate data. For dynamic obstacles, associate positional change parameters with time dimension to form digital volumes with temporal information. Align these obstacle models with the terrain grid according to their real-world spatial locations and use collision detection algorithms to ensure that their spatial relationships conform to physical logic and avoid spatial conflicts in the digital scene.

[0130] S230. Achieve model standardization and scene integration: Convert terrain meshes and obstacle models into a universal OBJ format, unify their coordinate system and data structure, use vertex coordinates to record spatial positions, use texture mapping to restore surface materials, use normal vectors to define surface orientation, and finally achieve digital mapping of real scenes to obtain three-dimensional digital scene models.

[0131] In a preferred embodiment of the present invention, the step of simulating multipath propagation using ray tracing based on a three-dimensional digital scene model and extracting channel statistical characteristic parameters, namely step S300, specifically includes:

[0132] S310. Configure ray emission and scene association: Set the three-dimensional coordinates of the signal source and receiver in the three-dimensional digital scene model, set the basic parameters, and define the ray propagation angle range and maximum tracking distance based on the distribution characteristics of scatterers in the scene to ensure that the ray can cover the area where all potential scatterers are located; the above basic parameters include, but are not limited to, the frequency of the transmitted signal, the initial power, and the antenna type.

[0133] S320. Begin ray tracing algorithm: Based on the geometric properties of the terrain mesh, obstacles, and scatterers in the 3D digital scene model, and the physical properties of their corresponding materials, locate the interaction point between the ray and the scatterer using a collision detection algorithm; combine the scattering angle, incident power, and scatterer material characteristics, derive the power attenuation of each scattering path according to scattering theory formulas; simultaneously record the propagation time of the ray from the emission source to the scatterer and then to the receiving point to determine the time delay parameter of the scattering path; extract the azimuth and elevation angles of the ray incident on the scatterer and the exit angle after scattering through vector calculation to form a complete set of angle parameters; the aforementioned geometric properties include, but are not limited to, dimensions, spatial coordinates, and surface normal vectors; the aforementioned physical properties of the material include, but are not limited to, scattering coefficient and dielectric constant;

[0134] S330. Aggregate channel statistical characteristic parameters: First, systematically classify and organize the loss power, delay, and angle data collected from all scattering paths. After obtaining the statistical sample set, first organize the parameters such as pitch angle, yaw angle, scatterer loss, and delay in the sample set, and remove outliers with too low energy or invalid values ​​to ensure data validity. Calculate the mean and variance of each parameter using statistical tools. For delay, the root mean square delay extension parameter also needs to be calculated. At the same time, draw the frequency histogram of each parameter. Based on the histogram, use the maximum likelihood estimation method to fit the probability density function (PDF), and obtain the cumulative distribution function (CDF) by integrating the probability density function. Finally, compare and verify the fitted probability density function and cumulative distribution function with typical channel models to determine the optimal fitting parameters, thereby obtaining the final channel statistical characteristics.

[0135] In this embodiment of the invention, based on the three-dimensional digital scene model constructed in step S200, the electromagnetic wave propagation path in the scene is simulated and calculated using the ray tracing method to construct a ray tracing model (i.e., RT model), and then the statistical characteristics of the channel (i.e., channel statistical features) are extracted; wherein, the channel statistical features include key parameters such as the power distribution of multipath components, time delay spread, and angle spread.

[0136] In practical applications, a 3D digital scene model constructed based on real urban scenarios includes 3D geometric structures such as high-rise buildings, streets, green belts, streetlights, and billboards, as well as dielectric constants, conductivity, and roughness parameters of materials such as concrete walls, glass curtain walls, asphalt pavements, and vegetation. When simulating the propagation path of electromagnetic waves using the ray tracing method, an initial ray beam is first generated from the UAV's antenna at preset angular intervals, while simultaneously setting the termination condition for ray propagation: the maximum number of reflections, diffractions, and scatterings is three. Subsequently, the propagation process of each ray in the scene is tracked: when specular reflection occurs on the high-rise building wall, the power loss is calculated based on the Fresnel reflection coefficient; when diffraction occurs at building corners and streetlight pole edges, the uniform diffraction theory (UTD) is used to calculate diffraction loss; when scattering occurs at green belt vegetation and rough walls, rays are generated in conjunction with the surface roughness model, and their power attenuates according to the scattering angle distribution. Then, based on the effective multipath components obtained from the above simulation, key parameters are extracted to calculate statistical characteristics: the scattered power distribution is calculated by statistically analyzing the received power of all multipaths, and the cumulative power distribution function is plotted; the delay spread is calculated by first calculating the propagation delay of each multipath, and then calculating the average delay and root mean square delay spread with the power of each multipath as the weight; the angle spread is calculated by decomposing the azimuth and elevation angles of each ray for the angle of arrival at the receiver, statistically analyzing the probability density function of the angle samples, and finally obtaining the root mean square angle spread through the second moment calculation, thus completing the extraction of channel statistical characteristics in urban scenarios.

[0137] In a preferred embodiment of the present invention, the step of constructing a statistical model of the drone's jitter characteristics and quantifying the drone's positional offset under jitter, specifically step S400, includes:

[0138] Statistical models are used to quantify the azimuth and pitch angle shifts caused by UAV jitter, and then the position change caused by the UAV is calculated based on geometric relationships. UAV jitter includes internal jitter caused by rotor vibration and external jitter caused by external factors. The position shift of the UAV under jitter is modeled as follows:

[0139] ;

[0140] In the formula, The positional offset of the drone under jitter. This refers to the positional offset of the drone due to internal jitter. This represents the positional offset of the drone under external jitter.

[0141] Specifically, regarding internal vibration, based on the structural parameters and dynamic characteristics of the UAV, and combined with historical vibration data, the periodic variation law of azimuth and pitch angles caused by mechanical motion such as rotors is described by the amplitude, frequency, and phase of a sine function. The specific model is as follows:

[0142] ;

[0143] In the formula, and These are the azimuth and pitch angles generated by the drone's shaking due to internal jitter. and This is the displacement caused by the swaying at that angle;

[0144] ;

[0145] In the formula, It is the internal jitter frequency.

[0146] In addition, external jitter factors include wind speed fluctuations and airflow disturbances. These factors are complex and independent, and wind speed fluctuations and airflow disturbances generally exhibit random distributions. However, the real-time compensation mechanism of the UAV flight control system ensures that the changes in azimuth and pitch angles statistically follow a normal distribution with zero mean and stable variance. The specific model is as follows:

[0147] ;

[0148] In the formula, and These are the azimuth and pitch angles generated by the drone's shaking under external vibrations. and This represents the displacement caused by the shaking at that angle.

[0149] like Figure 3 As shown, in a preferred embodiment of the present invention, the step of constructing a channel model that integrates attitude dynamics and jitter effects using an analytical method based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, i.e., step S500, specifically includes:

[0150] S510. Using analytical methods, derive the analytical expression for the time-varying characteristics of the channel. The channel impulse response formula is as follows:

[0151] ;

[0152] In the formula, H LoS and H NLoS These are the impulse response matrices for the direct and indirect beam paths, respectively. The transmitter uses a P-row, Q-column two-dimensional antenna array containing P*Q transmitting antennas; the receiver uses an A-row, B-column two-dimensional antenna array containing A*B receiving antennas. Therefore, each matrix contains P*Q*A*B elements, as detailed below:

[0153] ;

[0154] ;

[0155] In the formula, and Let be the channel impulse responses of the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna, respectively, expressed by the following formulas:

[0156] ;

[0157] ;

[0158] In the formula, and , respectively, are the antenna gains along the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna; and Let be the direct and indirect large-scale fading from the qp-th transmitting antenna to the ab-th receiving antenna, respectively. and These are the signal phase change factors for the direct and indirect paths from the qp-th transmitting antenna to the ab-th receiving antenna, respectively.

[0159] S520. Next, the parameters constituting the channel impulse response will be described separately. In the UAV air-to-ground communication system, the modeling of antenna gain needs to comprehensively consider the influence of platform dynamic characteristics and hardware losses. The effective aperture of the antenna is the equivalent area that the antenna can actually receive or transmit electromagnetic waves. Its actual effectiveness will be directionally attenuated due to UAV attitude deviations and high-frequency random jitter. When the antenna main lobe deviates from the target direction, the equivalent effective aperture can be regarded as the projection of the antenna physical aperture on the signal pointing direction, forming a gain attenuation factor related to the attitude error angle. At the same time, hardware losses correct the antenna gain through a loss factor independent of spatial pointing. The aforementioned attitude deviations include, but are not limited to, roll, pitch, and yaw angle changes.

[0160] ;

[0161] ;

[0162] The total antenna gain can be expressed as the product of the ideal gain, the effective aperture efficiency correction factor caused by attitude and jitter, and the hardware loss factor. The specific model is as follows:

[0163] ;

[0164] in, Indicates the wavelength of electromagnetic waves. and These are the time delays for direct and indirect beams, respectively. and It is the hardware loss of the transmitting antenna qp and the receiving antenna ab, which follows a chi-square distribution; and The effective apertures of the transmitting and receiving antennas under direct illumination conditions are given by the following formula:

[0165] ;

[0166] ;

[0167] in, and Let qp and ab be the normal directions of the antenna, respectively. Here is the attitude transformation matrix for the UAV. and They are time t and t+ respectively. Consider the distance of antenna qp-ab after jitter offset;

[0168] ;

[0169] ;

[0170] In the formula, It is the position of the transmitter at time t. It is the position of the receiver at time t. Indicates the spacing between the antenna elements at the transmitting end. Indicates the spacing between antenna elements at the receiving end;

[0171] Immediately afterwards, and The effective apertures of the transmitting and receiving antennas under non-direct illumination conditions are given by the following formula:

[0172] ;

[0173] in, The physical area of ​​the transmitting QP antenna. This represents the mutual coupling matrix of the transmitting antenna. The physical area of ​​the receiver's ab antenna. This represents the mutual coupling matrix of the receiving antenna. and They are time t and t+ respectively. Consider the non-direct path distance of antenna qp-ab after jitter offset;

[0174] ;

[0175] ;

[0176] in, and Representing t and scattering body at time Location;

[0177] Large-scale fading is used to characterize the slowly varying characteristics of signal strength in wireless communication over a large spatial and environmental range, including path loss and shadow fading (SH). The formulas for direct-path large-scale fading and indirect-path large-scale fading are as follows:

[0178] ;

[0179] ;

[0180] Path loss originates from the signal propagation distance and the propagation environment, causing energy to attenuate with distance or complex scenarios. It determines the basic signal coverage and received strength. Quantized attitude and jitter offsets are incorporated into the path loss calculation, specifically the path loss of the direct path. Path loss of non-direct path The formula is expressed as follows:

[0181]

[0182] ;

[0183] Where c=2 is the free space loss model; c=3~5 is the urban environment; and These are the direct and indirect trajectory lengths after incorporating quantized attitude and jitter parameters, respectively, expressed by the following formulas:

[0184] ;

[0185] ;

[0186] In the formula, This represents the position offset caused by the drone's attitude.

[0187] When radio waves encounter obstacles such as rolling hills, buildings, and forests, forming a shadow region and causing a slow change in the median signal strength, this phenomenon is called the shadowing effect. Statistically, it follows a log-normal distribution, and its probability density function can be expressed as:

[0188] ;

[0189] in, and These are the mean and variance of shadow fading, respectively.

[0190] S540. Next, the small-scale parameters of the channel are modeled. The small-scale parameters caused by jitter and attitude are represented by the position offsets modeled in the above statistical model. Here, only the influence caused by velocity is described. Based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics, the small-scale parameter model of the UAV is obtained as follows:

[0191] ;

[0192] ;

[0193] The Doppler frequency shift of the channel caused by the direct trajectory velocity is modeled as follows:

[0194] ;

[0195] in, and These are the velocities of the receiving and transmitting ends, respectively, and u is a unit vector representing the signal propagation direction. The specific formula is as follows:

[0196] ;

[0197] The Doppler frequency shift formula for the channel caused by non-direct path velocity is modeled as follows:

[0198] ;

[0199] ;

[0200] in, It is the speed of the transmitting end. The relative velocity between the transmitter and the scatterer It is the velocity of the scatterer n, expressed by the formula: Similarly, It is the speed of the transmitting end. The relative velocity between the transmitter and the scatterer It is the velocity of the scatterer n, expressed by the formula: , The angle between the signal propagation direction and the velocity direction of the transmitting end. It is the angle between the signal propagation direction and the reference direction of the receiving antenna array.

[0201] The time-domain graph of the channel impulse response (CIR) of the channel model obtained in practical applications according to the embodiments of the present invention is shown below. Figure 4 As shown, this channel is a time-varying fading channel, with prominent multipath effects and deep fading problems.

[0202] In another embodiment of the present invention, an air-to-ground channel dynamic modeling system for highly maneuverable unmanned aerial vehicles is also provided, for implementing the above-mentioned air-to-ground channel dynamic modeling method, comprising:

[0203] The dynamics model building module is used to build the dynamics model of the UAV and quantify the antenna position offset caused by changes in the UAV's attitude.

[0204] The scene model building module is used to build a 3D digital scene model that recreates the real scene based on the collected terrain and obstacle data.

[0205] The channel statistical characteristic parameter extraction module is used to simulate multipath propagation based on a three-dimensional digital scene model using the ray tracing method and extract channel statistical characteristic parameters.

[0206] The statistical model building module is used to build a statistical model of the drone's jitter characteristics and quantify the drone's positional offset under jitter.

[0207] The channel model construction module is used to construct a channel model that integrates attitude dynamics and jitter effects using analytical methods, based on the UAV dynamics model, channel statistical characteristic parameters, and statistical model of UAV jitter characteristics.

[0208] It should be noted that each of the above modules can be implemented as a computer program, which can run on a computer device. The computer device's memory can store the computer program that makes up each module, enabling the processor to execute each step of the above method.

[0209] It should be understood that although the steps in the flowcharts of the various embodiments of the present invention are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the various embodiments may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least a portion of the sub-steps or stages of other steps.

[0210] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Furthermore, any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory.

[0211] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of this patent should be determined by the appended claims.

Claims

1. A method for dynamic modeling of air-to-ground channel for high maneuverable unmanned aerial vehicles, characterized in that, The method comprises the following steps: constructing a UAV dynamics model and quantifying the antenna position offset caused by UAV attitude change; based on the collected terrain data and obstacle data, constructing a three-dimensional digital scene model that restores the real scene; based on the three-dimensional digital scene model, simulating multipath propagation by using a ray tracing method and extracting channel statistical characteristic parameters; constructing a statistical model of UAV jitter characteristics and quantifying the position offset of the UAV under the action of jitter; according to the UAV dynamics model, the channel statistical characteristic parameters and the statistical model of the UAV jitter characteristics, constructing a channel model that integrates the effects of attitude dynamics and jitter by using an analytical method.

2. The air-to-ground channel dynamic modeling method for high-maneuverable UAVs of claim 1, wherein, The step of constructing a UAV dynamics model and quantifying the antenna position offset caused by UAV attitude change specifically comprises: setting the basic physical parameters of the UAV as physical constraints, inputting the physical constraints, combining real-time IMU data, establishing a dynamics equation by using a discretized Euler integral method, realizing the conversion of forces and torques between the machine system and the inertial system by using a rotation matrix, iteratively updating the dynamics equation, thereby completing the complete mapping from the control input to the global motion trajectory, realizing the UAV dynamics modeling combining the physical parameter hard constraint and the real-time state soft feedback, and obtaining the UAV dynamics model; according to the attitude data output by the UAV dynamics model, establishing a position offset model caused by the UAV attitude, converting the position offset caused by the attitude angle into the pitch angle and the azimuth angle in the ground coordinate system by using an attitude matrix, and quantifying the antenna position offset.

3. The air-to-ground channel dynamic modeling method for high-maneuverable UAVs of claim 2, wherein, The dynamics equation is: ； wherein, denotes the angular velocity of the UAV at time t in the body coordinate system, denotes the roll angle, pitch angle and yaw angle of the body coordinate system relative to the inertial coordinate system at time t, J is the moment of inertia matrix, and W is the angular velocity conversion matrix, denotes the velocity of the UAV at time t in the ground coordinate system, denotes the position of the UAV at time t in the ground coordinate system; and denote the control torque and total pull force in the body coordinate system, respectively, denotes the rotation matrix of the body coordinate system, and m and g denote the mass of the UAV and the gravitational acceleration, respectively. The position offset model caused by the UAV attitude is: ； wherein is the position offset caused by the UAV attitude, and are the azimuth and the pitch angle, respectively, and is the displacement caused by the sway in the angle; the azimuth and the pitch angle models are as follows: ； wherein, , and are the pitch, roll and yaw angles of the UAV dynamics model output, respectively.

4. The method of claim 1, wherein, The step of constructing a three-dimensional digital scene model that restores the real scene based on the collected terrain data and obstacle data specifically comprises: digitally analyzing the real scene; wherein, for the elevation map, key information is extracted by algorithm: the elevation value is converted into the vertex height in the three-dimensional space, the geographic coordinates are mapped into the position parameters in the digital coordinate system, the physical characteristics of the real terrain are converted into calculable geometric data, an initial terrain framework is formed, and the terrain grid is obtained; classifying and modeling the analyzed real scene data and fusing them to obtain an obstacle model; wherein, for static obstacles, a geometric model is constructed by using their size and coordinate data; for dynamic obstacles, a digital body with time sequence information is formed by associating the position change parameters in the time dimension; the obstacle models are aligned according to the real space position and the terrain grid, and the spatial relationship between them is ensured to meet the physical logic by using a collision detection algorithm, so as to avoid spatial conflicts in the digital scene; standardizing the model and integrating the scene: the terrain grid and the obstacle model are converted into the general OBJ format, the coordinate system and the data structure are unified, the spatial position is recorded by using vertex coordinates, the ground surface material is restored by using texture mapping, and the surface orientation is defined by using a normal vector, so that the digital mapping of the real scene is finally realized, and the three-dimensional digital scene model is obtained.

5. The method of claim 1, wherein, The step of simulating multipath propagation by using a ray tracing method and extracting channel statistical characteristic parameters based on the three-dimensional digital scene model specifically comprises: Performing ray emission and scene association configuration: setting the three-dimensional coordinates of the signal emission source and receiving point in the three-dimensional digital scene model, setting the basic parameters, and based on the distribution characteristics of the scatterers in the scene, the angle range and maximum tracking distance of ray propagation are determined to ensure that the rays can cover all potential scatterer regions; the basic parameters include the frequency, initial power and antenna type of the transmitted signal; Starting the ray tracing algorithm: according to the geometric properties and material properties corresponding to the physical properties of the terrain grid, obstacles and scatterers in the three-dimensional digital scene model, the interaction points of the rays and the scatterers are located through the collision detection algorithm; combined with the scattering angle, incident power and scattering body material characteristics, the power attenuation of each scattering path is derived according to the scattering theory formula; the propagation time from the emission source to the scatterer and then to the receiving point is recorded simultaneously to determine the time delay parameter of the scattering path; the azimuth and elevation angles of the ray incident to the scatterer and the exit angle after scattering are extracted through vector calculation to form a complete set of angle parameters; the geometric properties include size, spatial coordinates and surface normal vector; the material corresponding physical properties include scattering coefficient and dielectric constant; Performing channel statistical characteristic parameter aggregation: first, all the loss power, time delay and angle data collected from the scattering paths are systematically classified and arranged; after obtaining the statistical sample set, first, the elevation angle, azimuth angle, scatterer loss and time delay parameters in the sample set are arranged, and the abnormal values with too low energy or invalidity are removed to ensure data validity; the mean and variance of each parameter are calculated through statistical tools, and the mean square root time delay expansion parameter is also calculated; at the same time, the frequency histogram of each parameter is drawn, and the probability density function is fitted based on the histogram using the maximum likelihood estimation method, and the cumulative distribution function is obtained by integrating the probability density function; finally, the fitted probability density function and cumulative distribution function are compared and verified with typical channel models to determine the optimal fitting parameters, thereby obtaining the final channel statistical characteristics.

6. The method of claim 1, wherein, A statistical model of the shaking characteristics of the unmanned aerial vehicle is constructed, and the steps of quantifying the position deviation of the unmanned aerial vehicle under the shaking effect are as follows, specifically comprising: The azimuth and elevation angle deviations caused by the shaking of the unmanned aerial vehicle are quantified using the statistical model, and the position change caused by the unmanned aerial vehicle is calculated according to the geometric relationship; the shaking of the unmanned aerial vehicle includes internal shaking caused by rotor vibration and external shaking caused by external factors; the position deviation of the unmanned aerial vehicle under the shaking effect is modeled as follows: ； In the formula, is a position offset of the UAV under the internal shaking action, is a position offset of the UAV under the internal shaking action, is a position offset of the UAV under the external shaking action.

7. The air-to-ground channel dynamic modeling method for high-maneuverable UAVs of claim 6, wherein, For internal shaking, based on the structural parameters and dynamic characteristics of the unmanned aerial vehicle, and combined with historical vibration data, the periodic variation law of the azimuth and elevation angle caused by the mechanical movement of the rotor is described through the amplitude, frequency and phase of the sine function, and the specific model is as follows: ； wherein, and are azimuth and pitch angles of the UAV shaking under the internal jitter effect, respectively, and is the displacement of the shaking in the angle. ； In the formula, is the internal dither frequency.

8. The air-to-ground channel dynamic modeling method for high-maneuverable UAVs of claim 6, wherein, For external shaking, the external shaking influencing factors include wind speed fluctuation and air flow disturbance; the external shaking influencing factors are randomly distributed, and the real-time compensation mechanism of the unmanned aerial vehicle flight control system makes the azimuth and elevation angle change obey a normal distribution with a mean of zero and a stable variance in statistics, and the specific model is as follows: ； wherein, and are azimuth and pitch angles of the UAV shaking under the external shaking effect, respectively, and is the displacement of the shaking in the angle.

9. The method of claim 1, wherein, According to the unmanned aerial vehicle dynamics model, channel statistical characteristic parameters and the statistical model of the unmanned aerial vehicle jitter characteristics, the steps of constructing the channel model fusing the attitude dynamics and jitter effects by using the analytical method, specifically comprising: The analytical expression of the channel time-varying characteristics is derived by using the analytical method, and the channel impulse response formula is as follows: ； where H LoS and H NLoS are the direct and non-direct channel impulse response matrices, respectively; the transmitter uses a two-dimensional antenna array with P rows and Q columns, which contains P*Q transmit antennas in total; the receiver uses a two-dimensional antenna array with A rows and B columns, which contains A*B receive antennas in total; therefore, each matrix includes P*Q*A*B elements, which are specifically expressed as follows: ； ； wherein and are the direct and non-direct channel impulse responses from the qpthtransmit antenna to the abrreceive antenna, respectively, and are expressed by the formulas: ； ； wherein and Gqpaband Gqpab, respectively, are the free-space and non-free-space antenna gains of the qpthtransmit antenna to the abrreceive antenna; and and and and Next, the parameters constituting the channel impulse response are described respectively. In the unmanned aerial vehicle air-ground communication system, the modeling of the antenna gain needs to consider the influence of the platform dynamic characteristics and the hardware loss; the antenna effective aperture is the equivalent area that the antenna can actually receive or transmit electromagnetic waves, and its actual performance will be directionally attenuated due to the attitude deviation and high-frequency random jitter of the unmanned aerial vehicle. When the main lobe of the antenna deviates from the target direction, the equivalent effective aperture is regarded as the projection of the physical aperture of the antenna in the signal pointing direction, forming a gain attenuation factor related to the attitude error angle; at the same time, the hardware loss corrects the antenna gain through a loss factor independent of the spatial pointing; the attitude deviation includes roll, pitch and yaw angle changes; ； ； The total antenna gain is expressed as the product of the ideal gain, the effective aperture efficiency correction factor caused by the attitude and jitter, and the hardware loss factor, and the specific model is as follows: ； wherein, denotes the electromagnetic wavelength, and are the direct path and non-direct path delays, respectively, and are the hardware losses of the transmitting antenna qp and the receiving antenna ab, respectively, and are subject to a chi-square distribution; and are the effective apertures of the transmitting antenna and the receiving antenna under the direct path condition, respectively, and are given by the following equations: ； ； wherein, and are the normal directions of the antennas qp and ab, respectively, is the attitude transformation matrix of the drone, and are the distances between the antennas qp and ab at time t and t+1, respectively, considering the jitter offset. ； ； In the formula, is the position offset of the UAV under the jitter effect, is the position of the transmitting end at time t, is the position of the receiving end at time t, represents the spacing between the antenna elements of the transmitting end, represents the spacing between the antenna elements of the receiving end; Next, and The effective aperture of the transmitting antenna and the receiving antenna under the non-orthodiametric condition is as follows: ； wherein, a physical area of the transmit-end antenna qp, denotes a mutual coupling matrix of the transmit-end antenna, a physical area of the receive-end antenna ab, denotes a mutual coupling matrix of the receive-end antenna, and are the non-line-of-sight distances of the antennas qp-ab at time t and t+1, respectively. considering the non-line-of-sight distances of the antennas qp-ab after the jitter offset; ； ； wherein, and respectively represent t and the position of the scatterer at time t. Large-scale fading is used to describe the slow variation characteristics of signal strength in wireless communication with large-scale space and environment, including path loss and shadow fading SH; the direct path large-scale fading formula and the non-direct path large-scale fading formula are as follows: ； ； The path loss is caused by the signal propagation distance extension and the propagation environment, and makes the energy attenuate with the distance or the complex scene, determines the basic coverage range and the receiving strength of the signal, and the quantized attitude and the jitter offset are integrated into the path loss calculation, and the path loss of the direct diameter and the path loss of the non-direct diameter is expressed as follows: ； ； Wherein, c=2 is free space loss model; c=3~5 is urban environment; and Ld and Lnd are the direct path length and non-direct path length after adding the quantized attitude and jitter parameters, and the formula is expressed as follows: ； ； In the formula, is the position offset caused by the UAV attitude; When the radio wave propagation path is blocked by obstacles, a shadow area of radio waves is formed, causing the slow change of the median value of the signal field strength, which is the shadow effect. From the statistical characteristics, it is subject to a lognormal distribution, and its probability density function can be expressed as: ； wherein and are the mean and variance of the shadow fading, respectively; Next, the small-scale parameters of the channel are modeled, wherein the small-scale parameters caused by jitter and attitude are represented by the position offset obtained by the above statistical modeling, and only the influence of the speed is described here. Based on the unmanned aerial vehicle dynamics model, channel statistical characteristic parameters and the statistical model of the unmanned aerial vehicle jitter characteristics, the small-scale parameter model of the unmanned aerial vehicle is obtained as follows: ； ； The Doppler shift formula of the channel caused by the speed of the direct path is modeled as follows: ； where, and are the velocities of the receiving and transmitting ends, respectively, and u is the unit vector in the direction of signal propagation, which is given by ； The Doppler shift formula of the channel caused by the speed of the non-direct path is modeled as follows: ； ； wherein, is the transmitting end moving speed, is the relative speed between the transmitting end and the scatterer, is the moving speed of the scatterer n, which is expressed by the formula: , similarly, is the transmitting end moving speed, is the relative speed between the transmitting end and the scatterer, is the moving speed of the scatterer n, which is expressed by the formula: ; is the angle between the signal propagation direction and the transmitting end speed direction, is the angle between the signal propagation direction and the receiving end antenna array reference direction.

10. A dynamic modeling system for air-to-ground channel for high-maneuverable unmanned aerial vehicles, for implementing the dynamic modeling method for air-to-ground channel according to any one of claims 1-9, characterized in that, Comprise: The dynamics model construction module is used to construct the dynamics model of the unmanned aerial vehicle and quantify the antenna position offset caused by the attitude change of the unmanned aerial vehicle; The scene model construction module is used to construct a three-dimensional digital scene model restoring the real scene based on the collected terrain data and obstacle data; The channel statistical characteristic parameter extraction module is used to simulate multipath propagation by using the ray tracing method based on the three-dimensional digital scene model, and extract the channel statistical characteristic parameters; The statistical model construction module is used to construct the statistical model of the unmanned aerial vehicle jitter characteristics and quantify the position offset of the unmanned aerial vehicle under the action of jitter; The channel model construction module is used to construct the channel model fusing the attitude dynamics and jitter effects by using the analytical method according to the dynamics model of the unmanned aerial vehicle, channel statistical characteristic parameters and the statistical model of the unmanned aerial vehicle jitter characteristics.

Images