A wireless channel simulation method and system based on time-slice differential trajectories

Abstract

This application discloses a wireless channel simulation method and system based on time-slice differential trajectories. The method includes: acquiring the motion time, initial coordinates, and initial velocities of the receiver and transmitter; dividing the motion time of the receiver and transmitter into multiple differential segments; configuring the parameters of the motion trajectories at both ends of the receiver and transmitter based on the differential segment time and the initial coordinates of the receiver and transmitter; and calculating the motion coordinates, transmit / receive angle change parameters, and multipath channel coefficients between the transmit and receive links of the receiver and transmitter. This application improves the efficiency of trajectory parameter configuration and subsequent channel generation based on differential parameter configuration in scenarios where both the receiver and transmitter are in irregular motion, and ensures the spatiotemporal continuity of the wireless channel during complex motion processes, providing convenience for the generation of complex trajectory channel data in multiple scenarios.

Description

Technical Field

[0001] This application belongs to the field of mobile communication wireless channel simulation technology, specifically relating to a wireless channel simulation method and system based on time-slice differential trajectories. Background Technology

[0002] With the successful commercialization of 5G technology and the increasing number of connected devices, mobile communication has officially entered the era of the Internet of Everything. The future development of 6G technology will continue the evolution from low frequency to millimeter wave, and then to high frequency bands with large bandwidth in terahertz. At that time, communication, sensing, intelligent computing, and other functions will achieve technological integration and network collaboration. Whether for industrial development upgrading from 5G technology or for 6G technology research, the theory and methods of multi-band broadband multi-scenario wireless channel modeling are important research directions in the field of mobile communication. Accurate channel modeling, through channel simulators, can effectively construct performance testing environments for wireless communication equipment, thereby enabling performance verification of wireless devices such as base stations, terminals, automobiles, and satellites, and helping to upgrade these devices.

[0003] The evolution of 5G technology will extend to more vertical industries, and the millimeter-wave to terahertz frequency bands of 6G will also create entirely new communication scenarios. The motion trajectories at both ends of the transceiver will become increasingly complex and irregular, such as the motion trajectory of a car in complex road conditions and the three-dimensional dynamic flight trajectory of a drone. Existing trajectory configuration methods for wireless channel simulators are based on combinations of regular trajectories such as straight lines, vertical and horizontal arcs. Even custom segmented configuration methods are based on inputting the coordinate values ​​of trajectory turning points, which makes the configuration of irregular and complex trajectories difficult. The calculation complexity of point coordinates during parameter configuration is high, and this point-based method will cause discontinuities in the subsequent channel coefficient calculation results. Handling discontinuities introduces additional complexity. Therefore, a new method is needed for complex trajectory parameter configuration and the calculation and generation of wireless channel data. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this application proposes a wireless channel simulation method and system based on time-slice differential trajectories. This application addresses the complex motion trajectories at both ends of the transceiver, offering both arbitrary configuration and speed. Simultaneously, it enables rapid channel model data calculation while configuring parameters and calculating and plotting the trajectory.

[0005] Firstly, this application proposes a wireless channel simulation method based on time-slice differential trajectories, including:

[0006] Obtain the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter;

[0007] The motion time of the receiver and transmitter is divided into multiple differential segments, and the parameters of the motion trajectories at both ends of the receiver and transmitter are configured based on the differential segment time and the starting coordinates of the receiver and transmitter.

[0008] Based on the configured parameters, the initial velocities of the receiver and transmitter, and the time of each differential segment, the coordinates of each sample point, the receiver velocity vector of each sample point, and the transmitter velocity vector are obtained.

[0009] Based on the coordinates of each sample point, the method of calculating within the segment sample points is used to obtain the principal diameter departure angle vector and principal diameter arrival angle vector of each sample point.

[0010] Based on the departure angle vector and arrival angle vector of the main path, the antenna gain vectors of each sample point on the two polarizations of the main path receiver end, the antenna gain vectors of each sample point on the two polarizations of the main path transmitter end, and the main path direction vectors of the receiver and transmitter are obtained.

[0011] Based on the adjusted antenna array coordinate vectors at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main path direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector at each sample point, the antenna gain vectors at each sample point on the two polarizations at the main path transmitter end, and the antenna gain vectors at each sample point on the two polarizations at the main path receiver end, the main path channel coefficients of the trajectory between the receiver and transmitter are determined.

[0012] Optionally, the parameters for determining the motion trajectories at both ends of the receiver and transmitter include:

[0013] The respective differential time intervals of the receiver and transmitter, the motion orientation vectors of the respective differential time intervals of the receiver and transmitter, the attitude orientation vectors of the respective differential time intervals of the receiver and transmitter, and the acceleration vectors of the respective differential time intervals of the receiver and transmitter.

[0014] Optionally, the motion orientation vectors of the receiver and transmitter, which are the differential segments of the parameters used to complete the motion trajectories at both ends of the receiver and transmitter, include:

[0015] The receiver's starting azimuth angle, starting elevation angle, ending azimuth angle, and ending elevation angle.

[0016] The transmitter's starting azimuth angle, starting elevation angle, ending azimuth angle, and ending elevation angle.

[0017] Optionally, the attitude orientation vectors of the receiver and transmitter, which are the differential segments of the parameters used to complete the motion trajectories at both ends of the receiver and transmitter, include:

[0018] The rotation angles of the receiver's starting three-dimensional coordinate axis, the rotation angles of the receiver's ending three-dimensional coordinate axis, the rotation angles of the transmitter's starting three-dimensional coordinate axis, and the rotation angles of the transmitter's ending three-dimensional coordinate axis.

[0019] Optionally, the adjusted antenna element coordinate vectors at both ends of the receiver and transmitter are calculated using the antenna element coordinate vectors before the transformation at both ends of the receiver and transmitter, the rotation angles of the three-dimensional coordinate axes of each sample point of the receiver, and the rotation angles of the three-dimensional coordinate axes of each sample point of the transmitter.

[0020] Optionally, the values ​​of the rotation angles of the three-dimensional coordinate axes of each sample point of the receiver and the three-dimensional coordinate axes of each sample point of the transmitter required to calculate the coordinate vectors of the antenna array after adjustment at both ends of the receiver and the transmitter are obtained by linear interpolation of the initial rotation angle value and the final rotation angle value within the respective differential time interval of the receiver and the transmitter.

[0021] Optionally, the receiver velocity vector of each sample point is calculated using the receiver's initial velocity within the differential time interval, the receiver's motion time within the differential time interval, the receiver's acceleration vector within the differential time interval, and the unit vector of the receiver's differential motion direction.

[0022] The transmitter velocity vector at each sample point is calculated using the transmitter's initial velocity within the differential time interval, the transmitter's motion time within the differential time interval, the transmitter's acceleration vector within the differential time interval, and the unit vector of the transmitter's differential motion direction.

[0023] Optionally, the antenna gain vectors of each sample point on the two polarizations at the main diameter transmitter end are calculated based on the antenna radiation pattern of the vertical and horizontal polarization global coordinate system at the transmitter end and the main diameter departure angle vector;

[0024] The antenna gain vectors at each sample point on the two polarizations of the main diameter receiver are calculated based on the antenna pattern in the vertical and horizontal polarization global coordinate system of the receiver and the main diameter arrival angle vector.

[0025] Secondly, a wireless channel simulation system based on time-slice differential trajectories is proposed, including:

[0026] The acquisition module is used to acquire the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter;

[0027] The calculation module is used to combine the data in the acquisition module and calculate the antenna array coordinate vectors after adjustment at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main diameter direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector at each sample point, the antenna gain vectors at each sample point on the two polarizations at the main diameter transmitter end, and the antenna gain vectors at each sample point on the two polarizations at the main diameter receiver end.

[0028] The determination module is used to determine the main path channel coefficients of the trajectory between the receiver and transmitter based on the calculation results of the calculation module.

[0029] Thirdly, a computer storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the wireless channel simulation methods based on time-slice differential trajectories in the first aspect.

[0030] Fourthly, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the steps of any of the wireless channel simulation methods based on time-slice differential trajectories in the first aspect.

[0031] The beneficial effects of the technical solutions provided in some embodiments of this application include at least the following:

[0032] This application provides a wireless channel simulation method and system based on time-slice differential trajectories. By configuring the motion time vectors of each differential segment, along with easily obtainable parameters such as direction, attitude, and acceleration, and combining the obtained motion time, initial coordinates, and initial velocities of the receiver and transmitter, the system calculates the coordinates of various points in the receiver and transmitter's motion trajectory. It then uses intra-segment interpolation and segmented splicing to rapidly calculate the parameters of the entire trajectory for each transceiver link, directly calculating the wireless channel coefficients of the entire trajectory using differential trajectories. This improves the efficiency of trajectory parameter configuration in scenarios where both the receiver and transmitter are in irregular motion, as well as the efficiency of subsequent channel generation based on differential parameter configuration. Furthermore, it ensures the spatiotemporal continuity of the wireless channel during complex motion processes, facilitating the generation of complex trajectory channel data for various scenarios such as vehicle-to-everything (V2X), drones, and satellites. Attached Figure Description

[0033] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0034] Figure 1 This is a flowchart of a wireless channel simulation method based on time-slice differential trajectories, as shown in an embodiment of this application.

[0035] Figure 2 This is a schematic diagram of the vehicle time-slice differential trajectory parameter configuration and base station to vehicle 3D scene applicable to the Internet of Vehicles scenario, as shown in the embodiments of this application.

[0036] Figure 3 This is a schematic diagram of the time-varying Doppler curves of the channel coefficients after parameter configuration based on the time-slice differential trajectory, as shown in an embodiment of this application.

[0037] Figure 4 This is a schematic diagram of the arc-shaped trajectory of a geostationary satellite in a satellite communication scenario, as shown in the embodiments of this application.

[0038] Figure 5 This is a schematic diagram of the receiver motion trajectory for millimeter-wave multi-user communication scenarios disclosed in the embodiments of this application. Detailed Implementation

[0039] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0040] In the following description, when referring to the accompanying drawings, the same numbers in different drawings denote the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of methods and apparatus consistent with some aspects of this application as detailed in the appended claims.

[0041] In the description of this application, it should be understood that the terms "first," "second," etc., are used for descriptive purposes only and should not be construed as indicating or implying relative importance. Those skilled in the art can understand the specific meaning of the above terms in this application based on the specific circumstances. Furthermore, in the description of this application, unless otherwise stated, "multiple" refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship.

[0042] This application proposes a wireless channel simulation method, device, and storage medium based on time-slice differential trajectories. By analyzing and decomposing a task, a resource allocation model for its support is established, and an improved marginal benefit algorithm is used to solve the model, thereby optimizing the allocation of product support resources.

[0043] The following description, in conjunction with the embodiments and accompanying drawings, provides further details.

[0044] Example 1

[0045] The following will combine Figure 1 This paper provides a detailed description of a wireless channel simulation method based on time-slice differential trajectories provided in the embodiments of this application.

[0046] A wireless channel simulation method based on time-slice differential trajectories includes the following steps:

[0047] Step S1: Obtain the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter;

[0048] Step S2: Divide the motion time of the receiver and transmitter into multiple differential segments, and configure the parameters of the motion trajectories at both ends of the receiver and transmitter based on the differential segment time and the starting coordinates of the receiver and transmitter;

[0049] In the specific implementation process, the differential trajectory configurator configures the conventional related parameters such as the number, frequency, sampling rate, initial position, attitude, and antenna configuration of the transceivers (i.e., receivers and transmitters in the following text). Unlike the two mainstream configuration methods of existing channel simulation configurators, which are either regular trajectory patterns using straight line segments plus arc segments, or configuration methods that input the coordinates of the turning points of each segment, this embodiment adopts a time-based differential segmentation method, which breaks down the complex motion trajectory of each transceiver into numerous small differential segments, and determines the number of differential segments S. n The total trajectory motion time T, and the motion time vector t within the differential segment. i ;in,

[0050]

[0051] Specifically, the configuration parameters include: the receiver's initial velocity scalar vr and the transmitter's initial velocity scalar v. t ; Receiver differential piecewise acceleration vector a r Transmitter differential acceleration vector at; receiver motion orientation vector differential segment start azimuth angle. Transmitter motion orientation vector differential segment start azimuth angle Receiver motion orientation vector differential segment start pitch angle Transmitter motion orientation vector differential segment start pitch angle Receiver motion orientation vector differential segment ending azimuth angle Transmitter motion orientation vector differential segment ending azimuth angle Receiver motion orientation vector differential segment end pitch angle Transmitter motion orientation vector differential segment end pitch angle The receiver differential segment's attitude orientation vector starts at the Z-axis rotation angle. The transmitter's differential segment attitude orientation vector starts at the Z-axis rotation angle. The attitude orientation vector of the receiver differential segment starts at the Y-axis rotation angle. Transmitter differential segment attitude orientation vector start Y-axis rotation angle The receiver differential segment's attitude orientation vector starts at the X-axis rotation angle. Transmitter differential segment attitude orientation vector start X-axis rotation angle The attitude orientation vector of the receiver differential segment ends at the Z-axis rotation angle. The transmitter's differential segment attitude orientation vector ends at the Z-axis rotation angle. Attitude orientation vector ending Y-axis rotation of receiver differential segment Transmitter differential segment attitude orientation vector end Y-axis rotation angle Attitude orientation vector ending X-axis rotation of receiver differential segment Transmitter differential segment attitude orientation vector end X-axis rotation angle The dimension of each of the above vectors is equal to S. n Dimension, each element is labeled i = 1, 2, ..., S n Based on the above configuration parameters, the total number of channel coefficient samples P can be obtained. n ;

[0052]

[0053] Δt=1 / f s (3)

[0054] Where Δt is the time sampling interval, f s To satisfy the update rate of the Nyquist sampling theorem;

[0055] Then the number of time-domain samples P of the channel coefficients of the i-th differential segment i The calculation formula is as follows:

[0056]

[0057] In the specific implementation process, the configuration of the differential trajectory parameters of each time slice in step S1 is completed through configuration panels, software interfaces, etc. The configuration methods include, but are not limited to, software interface configuration and parameter table import.

[0058] Step S3: Based on the configured parameters, the initial velocities of the receiver and transmitter, and the time of each differential segment, obtain the coordinates of each sample point, the receiver velocity vector of each sample point, and the transmitter velocity vector.

[0059] First, after configuration execution, the receiver vector coordinates [x] at the endpoint of each differential segment are set.r y r , z r ], transmitter vector coordinates [x t y t , z t ] and receiver motion distance vector d r The transmitter's motion distance vector d t Perform the calculation.

[0060] This embodiment uses a differential trajectory wireless channel calculator to calculate the transceiver coordinates. The iterative calculation of the coordinates of each sample point j within each differential segment i of the receiver is as follows:

[0061]

[0062] Where v ri Let v be the initial velocity of the receiver during each differential time interval, which is determined by the initial velocity v at the starting point. r1 The initial velocity v of the second differential segment is obtained by iteratively calculating the acceleration and motion time of each differential segment. r2 The calculation formula is:

[0063] v r2 =v r1 +a r1 t1 (6)

[0064] v ri The velocity vector of the receiver in each differential time interval is given by the unit vector of the motion direction in each differential segment. Initial velocity v in each differential time interval ri and acceleration a during each differential time interval ri The time of motion within the segment, jΔt, is jointly determined;

[0065] And so on:

[0066]

[0067] Sample point labels j = 1, 2, ..., P within the differential time interval i , where θ ri (j) is arrive The monotonically linear difference within the interval, by default

[0068] The coordinates of each sample point j within each differential segment i of the transmitter are calculated iteratively as follows:

[0069]

[0070] Where v ti Let v be the initial velocity of the receiver during each differential time interval, which is determined by the initial velocity v at the starting point. t1It is obtained by iterating over the acceleration and motion time of each differential segment;

[0071] The initial velocity v of the second differential segment t2 The formula for calculation is:

[0072] v t2 =v t1 +a t1 t1 (9)

[0073] Where v ti The velocity vector of the transmitter during each differential time interval is represented by the unit vector of the motion direction during each differential segment. Initial velocity v in each differential time interval ti and acceleration a during each differential time interval ti The motion time jΔt within the segment is jointly determined, therefore:

[0074]

[0075] The sample point labels j = 1, 2, ..., P within the differential time interval i , where θ ti (j) is arrive The monotonically linear difference within the interval, by default yes arrive The monotonically linear difference within the interval, by default

[0076] The starting coordinates of the receiver and transmitter allow for the configuration and plotting of their motion trajectories. Furthermore, for ease of configuration, a linkage relationship exists between the motion time and distance within each differential segment through the differential acceleration vector. Configuration can be achieved by adjusting the motion time, and it is also compatible with distance-based configuration methods. Specifically, the transceiver's motion distance is calculated, and the transmitter's motion distance d in each differential segment is... ti The calculation formula is as follows:

[0077]

[0078] The movement distance d of each differential segment receiver ri The calculation formula is as follows:

[0079]

[0080] The distance d during the calculation process ri d ti With motion time t i The values ​​can be converted to each other, and the distance between the starting and ending coordinates of the differential segment can be used to deduce the motion time.

[0081] Step S4: Based on the coordinates of each sample point, the main path departure angle vector and main path arrival angle vector of the entire trajectory of each link are obtained by using the intra-segment sample point calculation method;

[0082] Specifically, the wireless channel parameters based on the differential trajectory parameters at both ends of the transceiver are calculated using a differential trajectory wireless channel calculator. This includes calculating the departure angle vector of the principal path of the entire trajectory of each transceiver link based on the configuration parameters of step S2 within each differential segment time, using intra-segment interpolation and segmented splicing methods. Arrival angle vector of the principal path

[0083] The angle update method for each sampling point j within the i-th segment is as follows:

[0084]

[0085] in, Let j be the pitch angle of each sampling point j within the i-th segment. Let be the azimuth angle of each sampling point j within the i-th segment. Let be the elevation angle of each sampling point j within the i-th segment. Let be the azimuth angle of arrival of each sampling point j within the i-th segment.

[0086] Step S5: Based on the main path departure angle vector and the main path arrival angle vector, obtain the antenna gain vectors of each sample point on the two polarizations at the main path receiver end, the antenna gain vectors of each sample point on the two polarizations at the main path transmitter end, and the main path direction vectors of the receiver and transmitter; specifically including:

[0087] The gain vector F of each sample antenna on the two polarizations at the main path receiver end rx,u,θ ,

[0088]

[0089] Where F r,θ,gcs , This is the antenna pattern at the receiver end, perpendicular to the horizontally polarized global coordinate system.

[0090] The gain vector F of each sample point antenna on the two polarizations at the main transmitter end tx,s,θ ,

[0091]

[0092] Where F tx,θ,gcs , This is the antenna pattern in the vertical and horizontal polarization global coordinate system at the transmitter end.

[0093] Principal diameter direction vector They are respectively:

[0094]

[0095] The attitude orientation vector of the transceiver's differential segment includes the rotation angles of the starting three-dimensional coordinate axis and the ending three-dimensional coordinate axis; specifically:

[0096] Start Z-axis rotation Y-axis rotation X-axis rotation End of Z-axis rotation Y-axis rotation X-axis rotation

[0097] The rotation angles of the receiver's i-th segment and j-th sample point along the Z-axis, Y-axis, and X-axis are respectively α. r ,β r γ r The value of each sample point can be derived from the initial rotation angle value within segment i. and termination rotation angle value The results were obtained by performing linear interpolation separately.

[0098] The rotation angles of the i-th segment and the j-th sample point on the Z-axis, Y-axis, and X-axis are respectively α. t ,β t γ t The value of each sample point can be derived from the initial rotation angle value within segment i. and termination rotation angle value The results were obtained by performing linear interpolation separately.

[0099] The receiver-side antenna array coordinate vector is adjusted based on the receiver and transmitter attitude change parameters, including the X-axis, Y-axis, and Z-axis rotation parameters. for:

[0100]

[0101] in This is the coordinate vector of the antenna array before transformation at the receiver end;

[0102] The transmitter-side antenna array coordinate vector is adjusted based on the receiver and transmitter attitude change parameters, including the X-axis, Y-axis, and Z-axis rotation parameters. for:

[0103]

[0104] in This is the coordinate vector of the antenna array before the transmitter transformation.

[0105] Step S6: Based on the adjusted antenna array coordinate vectors at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main path direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector at each sample point, the antenna gain vectors at each sample point on the two polarizations at the transmitter end of the main path, and the antenna gain vectors at each sample point on the two polarizations at the receiver end of the main path, determine the main path channel coefficients of the trajectory between the receiver and transmitter link.

[0106] Specifically, the receiver velocity vector v of each sampling point obtained from step S3 r for

[0107] v r =v ri (j) (25)

[0108] The transmitter velocity vector vt at each sample point is:

[0109] v t =v ti (j) (26)

[0110] The main path channel coefficient of the complex trajectory link between the receiver antenna u and the transmitter antenna s is finally calculated using the following formula.

[0111]

[0112] in d3D is the wavenumber. Here, d3D represents the distance between the two ends of the transceiver, calculated based on the Euclidean distance between the current sample coordinates of the transmitter and receiver.

[0113] In the specific calculation process, based on the data and configuration parameters obtained in S1 and S2, the calculation of the differential segment motion distance and the transceiver endpoint coordinates, as well as the drawing of the trajectory lines, are completed. Those skilled in the art will understand that the trajectory types drawn by this method include, but are not limited to, regular trajectory types, such as arcs, polygons, continuous polylines, and irregular trajectory types. Based on the parameter configuration in step S2 and the calculation of the receiver and transmitter coordinates in step S3, the calculation of the wireless channel parameters and the multipath channel coefficients of the complex trajectory between the transceiver links described in step S5 is performed. The calculation methods include, but are not limited to, software background calculations, code script calculations, etc.

[0114] In the specific implementation process, in order to obtain a data format file that can be played by the channel simulator, the multipath channel coefficient data of complex trajectory between the transmit and receive links calculated by the differential trajectory wireless channel data generator is converted into a data format. This includes, but is not limited to, converting data stored in computer memory to .npy or .bin files, converting .npy output files to .bin output files, and other file formats. Finally, a data format file that can be played by the channel simulator at the transmit and receive ends under any trajectory is generated by segment splicing. The size of its segment storage can depend on the memory of the running hardware device.

[0115] It is understandable that the differential trajectory wireless channel data generator performs differential segment calculation data format conversion based on the channel coefficient calculation results obtained by the differential trajectory wireless channel calculator in step S2 above. The conversion methods include direct conversion and generation in the software background, and conversion again through data processing tools after generating a certain data format. Finally, a channel file suitable for playback by the channel simulator is generated by segment splicing. The size of its segment storage can depend on the memory of the running hardware device.

[0116] The method disclosed in this embodiment directly calculates the configuration vector parameters related to the entire trajectory through differential trajectory stitching. Finally, based on the continuous configuration parameters of the entire trajectory, a segmented acceleration algorithm is used to complete the wireless channel calculation under complex trajectories. The size of the segmented storage depends on the memory of the running hardware device. This solves the problem of difficult configuration of point position parameters for complex trajectories in existing wireless channel simulators, and improves the efficiency of dynamic wireless channel modeling and channel data generation for receivers and transmitters with complex trajectories.

[0117] To illustrate this embodiment in more detail, three specific examples are provided below:

[0118] Example 1

[0119] For scenarios involving communication between vehicles and base stations during vehicle-to-everything (V2X) driving, the specific technical solution for a vehicle turning at an overpass intersection, from north to south, changing lanes onto a ramp, and then spiraling uphill from the ramp to travel from west to east, is as follows:

[0120] Step 1: After configuring the transceiver quantity, frequency and sampling rate, initial position, attitude and antenna configuration, etc., through the differential trajectory configurator, as follows: Figure 2 As shown, a time-based differential segmentation method is used to break down the complex motion trajectory of the car receiver into numerous small differential segments. A differential track configurator is then used to configure the motion time vector t (in seconds) for each differential segment. The motion orientation vector of each differential segment includes the starting azimuth angle. Pitch angle End azimuth Pitch angle The attitude orientation vector of the receiver differential segment includes the starting Z-axis rotation angle. Y-axis rotation X-axis rotation End of Z-axis rotation Y-axis rotation X-axis rotation All angles mentioned above are in degrees; receiver acceleration vector a r The unit is meters per second squared. It can be seen that because the differential trajectory configurator divides the irregular curve into numerous differential segments, and the angles of each differential segment are continuous, the final configured parameters satisfy the spatiotemporal consistency of the wireless channel. The main parameter configuration table is as follows:

[0121] Table 1. Configuration Parameters for Time Slice Trajectory Differentiation Segmentation

[0122]

[0123] Step 2: After configuration execution, use the differential trajectory wireless channel calculator to calculate the motion distance vector D of each differential segment of the transceiver and the coordinates [x, y, z] of the transceiver's endpoint vector. The unit is meters. The calculation results can be found in [link to calculation]. Figure 2 Ultimately, the configuration of the vehicle terminal's motion trajectory and the drawing of the trajectory lines are realized, while the base station coordinate position remains unchanged. The wireless channel parameter calculation of the differential trajectory parameters is performed by the differential trajectory wireless channel calculator, including the calculation of the main path departure angle vectors AOD and ZOD, and the main path arrival angle vectors AOA and ZOA of the full trajectory of each transceiver link based on the configuration parameters of step 2 within each differential segment time, using intra-segment interpolation and segment splicing methods based on formulas (13) to (16); and the calculation of the receiver and transmitter direction vectors based on formulas (21) to (22). The rapid calculation of the main path antenna gain vector parameters of the two polarizations of the receiver and transmitter is realized based on formulas (17) to (20). All parameters of each sample point are calculated based on the formulas of the specific embodiment, and finally the calculation of the main path channel coefficient of the complex trajectory between the transceiver links is realized through formula (27). Since the parameters configured by the differential trajectory configurator already have spatiotemporal consistency, no additional processing is required when calculating the wireless channel coefficient through the parameters, thus speeding up the simulation time.

[0124] Step 3: Using a differential trajectory wireless channel data generator, the data format of the complex trajectory multipath channel coefficient data between the transmit and receive links calculated by differential segmentation is converted, including but not limited to converting data stored in computer memory to .npy or .bin files, and .npy files to .bin files, generating data format files that can be played by a channel simulator at the transceiver end under any trajectory. Simultaneously, the generator can generate data based on the channel and verify the main channel parameters through digital signal processing, such as... Figure 3 As shown, the Doppler of the generated channel link is verified by short-time Fourier transform. The frequency offset trend of the time-varying Doppler is consistent with expectations, and its change process is continuous and consistent.

[0125] Example 2

[0126] This example addresses communication between a satellite and a land-based terminal in a satellite communication scenario. The communication satellite is a geostationary orbit satellite, located 35,786.043 kilometers above the Earth, with a projected position of 101.4 degrees east longitude and 0 degrees north latitude. The satellite terminal is located at 22 degrees 17 minutes 37 seconds north latitude (22.2936 degrees) and 114 degrees 10 minutes 31 seconds east longitude (114.1753 degrees). The specific technical solution is implemented as follows:

[0127] Step 1: After configuring the conventional parameters such as the number of satellites and ground terminals, frequency and sampling rate, initial position, attitude and antenna configuration through the differential trajectory configurator, the motion trajectories of the ground terminals and satellites are divided into numerous small differential segments using a time-based differential segmentation method. The motion time vector t of the transceiver within each differential segment is then configured using the differential orbit configurator. r , t t The motion orientation vectors of the ground terminal and the satellite, each with its own differential segment, include the starting azimuth angle. Pitch angle End azimuth Pitch angle The attitude orientation vectors of the ground terminal and the satellite, each in its differential segment, include the starting Z-axis rotation angle. Y-axis rotation X-axis rotation End of Z-axis rotation Y-axis rotation X-axis rotation The acceleration vector 'a' of the ground terminal and the satellite within their respective differential time intervals r a t .

[0128] Step 2: After configuration execution, use the differential trajectory wireless calculator to complete the calculation of the motion distance vector d of each differential segment of the ground terminal and the satellite. r d t The endpoint ground terminal and satellite vector coordinates [x r y r , z r ]、[x t y t , z t The calculations are used to configure the motion trajectories and draw the trajectory lines at both the ground terminal and the satellite. For example... Figure 4A schematic diagram of the arc trajectory of a geostationary satellite after 31 minutes and 40 seconds of motion is drawn. The wireless channel parameters based on the differential trajectory parameters at both ends of the transceiver are calculated using a differential trajectory wireless channel calculator. This includes configuring the motion coordinates in step 2 based on each differential segment time. Intra-segment interpolation and segmented splicing are used to calculate the main path departure angle vectors AOD and ZOD, and the main path arrival angle vectors AOA and ZOA for each transceiver link, based on formulas (13) to (16). The direction vectors of the receiver and transmitter are calculated based on formulas (21) to (22). The gain vector parameters of the main path antennas on the two polarizations of the receiver and transmitter are quickly calculated based on formulas (17) to (20). All parameters for each sample point are calculated based on the formulas in the specific embodiment. Finally, the main path channel coefficients for the complex trajectory between the transceiver links are calculated using formula (27).

[0129] Step 3: Using the differential trajectory wireless channel data generator, the data format of the complex trajectory multipath channel coefficient data between each transceiver link calculated by differential segment is converted, including but not limited to converting data stored in computer memory to .npy or .bin files, and converting .npy files to .bin files, to generate a data format file that can be played by the channel simulator at the transceiver end under any trajectory.

[0130] Example 3

[0131] This example targets a millimeter-wave multi-user mobile communication scenario. The communication signal frequency is 26GHz and the terminal is stationary. The millimeter-wave terminal moves irregularly around the millimeter-wave base station. In this use case, there is one base station and two terminals. The specific technical solution is implemented as follows:

[0132] Step 1: After configuring the routine parameters such as the number of transceivers, millimeter-wave frequency and sampling rate, initial position, attitude, and antenna configuration through the differential trajectory configurator, a time-based differential segmentation method is used to break down the complex motion trajectory of each receiver into numerous small differential segments. The motion time vector t, in seconds, is then configured for each receiver within each differential segment using the differential trajectory configurator. The motion orientation vector of each receiver differential segment includes the starting azimuth angle. Pitch angle θ s Ending azimuth Pitch angle θ e The attitude orientation vector of the receiver differential segment includes the initial Z-axis rotation angle α. s Y-axis rotation angle β s X-axis rotation angle γ s End of Z-axis rotation angle α e Y-axis rotation angle β e X-axis rotation angle γ e All angles are in degrees; the receiver acceleration vector a is in meters per second squared.

[0133] Step 2: After configuration execution, the differential trajectory wireless channel calculator calculates the motion distance vector D for each differential segment of each receiver, as well as the coordinates [x, y, z] of the endpoint vector of each receiver segment, in meters. The final result is as follows: Figure 5 The configuration of the receiver's motion trajectory and the drawing of the trajectory lines are shown. The coordinate position of the transmitter base station remains unchanged. The motion trajector of different receivers in the air and on the ground is different, but the motion time is the same. It can be seen that since the differential trajectory configurator divides the irregular curve into many differential segments and the angle of each differential segment is continuous, the final configured parameters meet the spatiotemporal consistency of the wireless channel. The wireless channel parameters of the differential trajectory parameters are calculated by the differential trajectory wireless channel calculator, including the configuration parameters of the moving coordinates in step 2 based on the time of each differential segment. The calculation of the main path departure angle vectors AOD and ZOD, and the main path arrival angle vectors AOA and ZOA of the full trajectory of each transceiver link is completed based on formulas (13) to (16) using the segment interpolation and segment splicing method. The direction vectors of the receiver and transmitter are calculated based on formulas (21) to (22). The gain vector parameters of the main path antennas on the two polarizations of the receiver and transmitter are calculated quickly based on formulas (17) to (20). Based on the formulas of the specific embodiments, all parameters of each sample point are calculated, and finally the main path channel coefficient of the complex trajectory between the transmit and receive links is calculated through formula (27). Since the parameters configured by the differential trajectory configurator already have spatiotemporal consistency, no additional processing is required when calculating the wireless channel coefficients through the parameters, thus speeding up the simulation time. In addition, the parallel link processing method between multiple users can also speed up the simulation.

[0134] Step 3: Using the differential trajectory wireless channel data generator, the data format of the complex trajectory multipath channel coefficient data between the transmit and receive links calculated by differential segmentation is converted, including but not limited to converting data stored in computer memory to .npy or .bin files, and .npy files to .bin files. Finally, by segmenting and splicing, a data format file that can be played by the channel simulator is generated for each transmit and receive link under any trajectory.

[0135] Example 2

[0136] This application proposes a wireless channel simulation system based on time-slice differential trajectories, including:

[0137] include:

[0138] The acquisition module is used to acquire the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter;

[0139] The calculation module is used to combine the data in the acquisition module and calculate the antenna array coordinate vectors after adjustment at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main diameter direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector at each sample point, the antenna gain vectors at each sample point on the two polarizations at the main diameter transmitter end, and the antenna gain vectors at each sample point on the two polarizations at the main diameter receiver end.

[0140] The determination module is used to determine the main path channel coefficients of the trajectory between the receiver and transmitter based on the calculation results of the calculation module.

[0141] The wireless channel simulation system based on time-slice differential trajectories provided in this embodiment is implemented by the wireless channel simulation method based on time-slice differential trajectories provided in Embodiment 1. It has the same technical features as the wireless channel simulation method based on time-slice differential trajectories provided in Embodiment 1, so it can solve the same technical problems and achieve the same technical effects.

[0142] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0143] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0144] Example 3

[0145] This application proposes a computer storage medium storing a computer program that, when executed by a processor, implements the steps of the wireless channel simulation method based on time-slice differential trajectories in any embodiment.

[0146] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0147] Therefore, this application also proposes a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the methods described in any embodiment of this application. The computer-readable storage medium can be configured in any device of this application.

[0148] Example 4

[0149] This application also proposes an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method described in any embodiment of this application.

[0150] For example, it includes: one or more processors; and a storage device for storing one or more programs, which, when executed by the one or more processors, cause the one or more processors to implement the method provided in the embodiments of this application. The methods described are included in the functional descriptions above and will not be repeated here.

[0151] The electronic device also includes input and output devices; the processor, storage device, input and output devices in the electronic device can be connected by a bus or other means.

[0152] A storage device, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, and module units, such as program instructions corresponding to the methods in the embodiments of this application. The storage device may mainly include a program storage area and a data storage area. The program storage area may store the operating system and at least one application program required for a function; the data storage area may store data created based on the use of the terminal, etc. Furthermore, the storage device may include high-speed random access memory and non-volatile memory, such as at least one disk storage device, flash memory device, or other non-volatile solid-state storage device. In some instances, the storage device may further include memory remotely located relative to the processor, and these remote memories can be connected via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof. The various embodiments in this application are described in a progressive manner; similar or identical parts between embodiments can be referred to mutually. Each embodiment focuses on describing the differences from other embodiments.

[0153] Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or part of the technical solution, can be embodied in the form of a computer program product.

[0154] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0155] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0156] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0157] It should be noted that the flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions marked in the blocks may occur in a different order than those marked in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.

[0158] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The scope of protection of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the scope of the technology disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

[0159] The applicant has provided a detailed description of the implementation examples of this application in conjunction with the accompanying drawings. However, those skilled in the art should understand that the above implementation examples are merely preferred embodiments of this application. The detailed description is only intended to help readers better understand the spirit of this application and is not intended to limit the scope of protection of this application. On the contrary, any improvements or modifications made based on the inventive spirit of this application should fall within the scope of protection of this application.

Claims

1. A wireless channel simulation method based on time-slice differential trajectories, characterized in that, include: Obtain the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter; The motion time of the receiver and transmitter is divided into multiple differential segments, and the parameters of the motion trajectories at both ends of the receiver and transmitter are configured based on the differential segment time and the starting coordinates of the receiver and transmitter. Based on the configured parameters, the initial velocities of the receiver and transmitter, and the time of each differential segment, the coordinates of each sample point, the receiver velocity vector of each sample point, and the transmitter velocity vector are obtained. Based on the coordinates of each sample point, the method of calculating within the segment sample points is used to obtain the principal diameter departure angle vector and principal diameter arrival angle vector of each sample point. Based on the departure angle vector and arrival angle vector of the main path, the antenna gain vectors of each sample point on the two polarizations of the main path receiver end, the antenna gain vectors of each sample point on the two polarizations of the main path transmitter end, and the main path direction vectors of the receiver and transmitter are obtained. Based on the adjusted antenna array coordinate vectors at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main path direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector at each sample point, the antenna gain vectors at each sample point on the two polarizations at the main path transmitter end, and the antenna gain vectors at each sample point on the two polarizations at the main path receiver end, the main path channel coefficients of the trajectory between the receiver and transmitter are determined.

2. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The parameters for completing the motion trajectories at both ends of the receiver and transmitter include: The respective differential time intervals of the receiver and transmitter, the motion orientation vectors of the respective differential time intervals of the receiver and transmitter, the attitude orientation vectors of the respective differential time intervals of the receiver and transmitter, and the acceleration vectors of the respective differential time intervals of the receiver and transmitter.

3. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The motion orientation vectors of the receiver and transmitter, which are the differential segments of the parameters used to complete the motion trajectories at both ends of the receiver and transmitter, include: The receiver's starting azimuth angle, starting elevation angle, ending azimuth angle, and ending elevation angle. The transmitter's starting azimuth angle, starting elevation angle, ending azimuth angle, and ending elevation angle.

4. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The attitude orientation vectors of the receiver and transmitter, which are the differential segments of the parameters used to complete the motion trajectories at both ends of the receiver and transmitter, include: The rotation angles of the receiver's starting three-dimensional coordinate axis, the rotation angles of the receiver's ending three-dimensional coordinate axis, the rotation angles of the transmitter's starting three-dimensional coordinate axis, and the rotation angles of the transmitter's ending three-dimensional coordinate axis.

5. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The adjusted antenna array coordinate vectors at both ends of the receiver and transmitter are calculated using the original antenna array coordinate vectors at both ends of the receiver and transmitter, the rotation angles of the three-dimensional coordinate axes of each sample point of the receiver, and the rotation angles of the three-dimensional coordinate axes of each sample point of the transmitter.

6. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The values ​​of the rotation angles of the three-dimensional coordinate axes of each sample point of the receiver and the three-dimensional coordinate axes of each sample point of the transmitter required to calculate the coordinate vectors of the antenna array after adjustment at both ends of the receiver and the transmitter are obtained by linear interpolation of the initial rotation angle value and the final rotation angle value of each differential time segment of the receiver and the transmitter.

7. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The receiver velocity vector at each sample point is calculated using the receiver's initial velocity within the differential time interval, the receiver's motion time within the differential time interval, the receiver's acceleration vector within the differential time interval, and the unit vector of the receiver's differential motion direction. The transmitter velocity vector at each sample point is calculated using the transmitter's initial velocity within the differential time interval, the transmitter's motion time within the differential time interval, the transmitter's acceleration vector within the differential time interval, and the unit vector of the transmitter's differential motion direction.

8. The wireless channel simulation method based on time-slice differential trajectories according to claim 1, characterized in that, The antenna gain vectors of each sample point on the two polarizations of the main diameter transmitter end are calculated based on the antenna radiation pattern of the vertical and horizontal polarization global coordinate system of the transmitter end and the main diameter departure angle vector. The antenna gain vectors at each sample point on the two polarizations of the main diameter receiver are calculated based on the antenna pattern in the vertical and horizontal polarization global coordinate system of the receiver and the angle of arrival vector of the main diameter.

9. A wireless channel simulation system based on time-slice differential trajectories, characterized in that, include: The acquisition module is used to acquire the motion time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial velocity of the receiver and transmitter, and the distance between the two ends of the receiver and transmitter; The calculation module is used to divide the motion time of the receiver and transmitter into multiple differential segments. Based on the differential segment time and the initial coordinates of the receiver and transmitter, it configures the parameters of the motion trajectory at both ends of the receiver and transmitter. Based on the configured parameters, the initial velocities of the receiver and transmitter, and the time of each differential segment, it obtains the coordinates of each sample point, the receiver velocity vector, and the transmitter velocity vector of each sample point. Based on the coordinates of each sample point, it uses an intra-segment sample point calculation method to obtain the principal diameter departure angle vector and the principal diameter arrival angle vector of each sample point. Based on the principal diameter departure angle vector and the principal diameter arrival angle vector, it obtains the antenna gain vectors of each sample point on the two polarizations at the receiver end of the principal diameter, the antenna gain vectors of each sample point on the two polarizations at the transmitter end of the principal diameter, and the principal diameter direction vectors of the receiver and transmitter. The determination module is used to determine the main path channel coefficients of the trajectory between the receiver and transmitter based on the adjusted antenna array coordinate vectors at both ends of the receiver and transmitter, the distance between the two ends of the receiver and transmitter, the main path direction vectors of the receiver and transmitter, the receiver velocity vector and transmitter velocity vector of each sample point, the antenna gain vectors of each sample point on the two polarizations of the main path transmitter end and the antenna gain vectors of each sample point on the two polarizations of the main path receiver end.

10. A computer storage medium having a computer program stored thereon, characterized in that, When executed by the processor, the program implements the steps of any one of the wireless channel simulation methods based on time-slice differential trajectories as claimed in claims 1 to 8.

11. An electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method as described in any one of claims 1 to 8.

Images