Wireless channel simulation method and system based on time slice differential trajectory

Abstract

The invention discloses a wireless channel simulation method and system based on a time slice differential trajectory. The method comprises the following steps: acquiring motion time of a receiver and a transmitter, initial coordinates of the receiver and the transmitter, and initial speeds of the receiver and the transmitter; splitting the movement time of the receiver and the transmitter into a plurality of micro-segments, and configuring parameters of movement tracks at the two ends of the receiver and the transmitter based on the micro-segment time and the initial coordinates of the receiver and the transmitter; and calculating motion coordinates of the receiver and the transmitter, a transmitting and receiving angle change parameter and a multipath channel coefficient between transmitting and receiving links. According to the method and the device, the track parameter configuration efficiency of irregular motion scenes at the receiver and the transmitter and the subsequent channel generation efficiency based on differential parameter configuration are improved, the space-time continuity of a wireless channel in a complex motion process is ensured, and convenience is provided for multi-scene complex track channel data generation.

Description

Technical Field

[0001] The present application belongs to the technical field of mobile communication wireless channel simulation, and specifically relates to a wireless channel simulation method and system based on time slice differential trajectory. Background Art

[0002] With the successful commercialization of 5G technology and the increasing number of connected devices, mobile communications have officially entered the era of the Internet of Everything. In the future, the development of 6G technology will continue to evolve from low frequency to millimeter wave, and then to the high frequency band and large bandwidth of terahertz. By then, multiple functions such as communication, perception, and intelligent computing will achieve technical integration and network collaboration. Whether it is for the industrial development of 5G technology upgrades or for 6G technology research, the theory and method of multi-band broadband and multi-scenario wireless channel modeling is an important research direction in the field of mobile communications. Accurate channel modeling can effectively build a performance test environment for wireless communication equipment by simulating the wireless environment through a channel simulator, thereby realizing the performance verification of wireless equipment such as base stations, terminals, automobiles, and satellites, and helping equipment to upgrade technology.

[0003] The evolution of 5G technology will be oriented to more vertical industries. In the future, the millimeter wave to terahertz frequency band for 6G will also generate new communication scenarios. The motion trajectories at both ends of the transceiver will also become increasingly complex and irregular, such as the motion trajectory of a car under complex road conditions, and the three-dimensional dynamic flight trajectory of a drone. The trajectory configuration method of the existing wireless channel simulator is based on a combination of regular trajectories such as straight lines, vertical and horizontal arcs. Even if there is a customized segmented configuration method, it is based on the input of the coordinate values ​​of the trajectory turning points, which makes it difficult to configure irregular complex trajectories. The calculation complexity of the dotted coordinates is high during the parameter configuration process, and this coordinate dot-dot matching method will cause discontinuity in the subsequent channel coefficient calculation results. If the discontinuity problem is handled, additional complexity will be introduced. Therefore, a new method is needed to configure complex trajectory parameters and calculate and generate wireless channel data. Summary of the invention

[0004] In order to solve the deficiencies in the prior art, this application proposes a wireless channel simulation method and system based on time-sliced ​​differential trajectory. This application targets the complex motion trajectory of both ends of the transceiver, has the characteristics of both configuration arbitrariness and quickness, and can achieve fast channel model data calculation while completing parameter configuration and trajectory calculation and drawing.

[0005] In a first aspect, the present application proposes a wireless channel simulation method based on time slice differential trajectory, comprising:

[0006] Obtain the movement time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial speed of the receiver and transmitter, and the distance between the receiver and transmitter;

[0007] The motion time of the receiver and the transmitter is divided into multiple micro-segments, and the parameters of the motion trajectory of the receiver and the transmitter are configured based on the micro-segment time and the starting coordinates of the receiver and the transmitter;

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

[0009] Based on the coordinates of each sample point, the main path departure angle vector and the main path arrival angle vector of each sample point are obtained by using the intra-segment sample point calculation method;

[0010] Based on the main path departure angle vector and the main path arrival angle vector, 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 the transmitter are obtained;

[0011] Based on the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter, the distance between the receiver and the transmitter, the main path direction vectors of the receiver and the transmitter, the receiver velocity vector and the 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, the main path channel coefficient of the trajectory between the receiver and the transmitter link is determined.

[0012] Optionally, the parameters for completing the motion trajectory of both ends of the receiver and the transmitter include:

[0013] The respective micro-segment times of the receiver and the transmitter, the respective micro-segment motion orientation vectors of the receiver and the transmitter, the respective micro-segment attitude orientation vectors of the receiver and the transmitter, and the respective micro-segment acceleration vectors of the receiver and the transmitter.

[0014] Optionally, the motion orientation vectors of the micro-segments of the receiver and the transmitter in the parameters for completing the motion trajectories of both ends of the receiver and the transmitter include:

[0015] The starting azimuth angle of the receiver, the starting elevation angle of the receiver, the ending azimuth angle of the receiver, and the ending elevation angle of the receiver.

[0016] The starting azimuth angle of the transmitter, the starting elevation angle of the transmitter, the ending azimuth angle of the transmitter, and the ending elevation angle of the transmitter.

[0017] Optionally, the micro-segmented attitude orientation vectors of the receiver and the transmitter in the parameters for completing the motion trajectory of both ends of the receiver and the transmitter include:

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

[0019] Optionally, the adjusted antenna array coordinate vector at both ends of the receiver and the transmitter is calculated by the antenna array coordinate vector before transformation at both ends of the receiver and the transmitter, the rotation angle of the three-dimensional coordinate axis of each sample point of the receiver and the rotation angle of the three-dimensional coordinate axis of each sample point of the transmitter.

[0020] Optionally, 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 required to calculate the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter are obtained by linear interpolation of the starting rotation angle values ​​and the ending rotation angle values ​​within the respective micro-segment times of the receiver and the transmitter.

[0021] Optionally, the receiver velocity vector of each sample point is calculated by the initial velocity of the receiver within the micro-segment time, the movement time of the receiver within the segment, the acceleration vector of the receiver within the micro-segment time, and the unit vector of the micro-segment movement direction of the receiver;

[0022] The transmitter velocity vector of each sample point is calculated by the transmitter's initial velocity within the micro-segment time, the transmitter's movement time within the segment, the transmitter's acceleration vector within the micro-segment time, and the transmitter's micro-segment movement direction unit vector.

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

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

[0025] In the second aspect, a wireless channel simulation system based on time-sliced ​​differential trajectory is proposed, including:

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

[0027] A calculation module is used to combine the data in the acquisition module, and by splitting the movement time of the receiver and the transmitter into multiple micro-segments, calculate the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter, the distance between the two ends of the receiver and the transmitter, the main path direction vector of the receiver and the transmitter, the receiver velocity vector and the transmitter velocity vector of each sample point, the antenna gain vector of each sample point on the two polarizations of the main path transmitter end, and the antenna gain vector of each sample point on the two polarizations of the main path receiver end;

[0028] The determination module is used to determine the main path channel coefficient of the trajectory between the receiver and the transmitter link according to the calculation result of the calculation module.

[0029] According to a third aspect, a computer storage medium is provided, on which a computer program is stored, and when the program is executed by a processor, the steps of any one of the wireless channel simulation methods based on time-sliced ​​differential trajectories in the first aspect are implemented.

[0030] In a fourth aspect, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the steps of any one of the wireless channel simulation methods based on time-sliced ​​differential trajectories in the first aspect are implemented.

[0031] The beneficial effects brought about by the technical solutions provided by some embodiments of the present application include at least:

[0032] The present application provides a wireless channel simulation method and system based on time-sliced ​​differential trajectories, which configures the motion time vectors, direction, attitude, acceleration and other easily obtainable parameters of each differential segment; then combines the acquired motion time of the receiver and the transmitter, the starting coordinates of the receiver and the transmitter, and the initial speed of the receiver and the transmitter to complete the coordinate calculation of each sample point in the motion trajectory of the receiver and the transmitter, and adopts intra-segment interpolation and segment splicing to complete the rapid calculation of the parameters of the full trajectory of each transceiver link, and directly completes the calculation of the wireless channel coefficient of the entire trajectory in the form of differential trajectories; improves the efficiency of trajectory parameter configuration for irregular motion scenes at both ends of the receiver and the transmitter, and the efficiency of subsequent channel generation based on differential parameter configuration, and ensures the space-time continuity of the wireless channel during complex motion, which provides convenience for the generation of complex trajectory channel data for multiple scenes of Internet of Vehicles, drones, and satellites. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments of the present application or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0034] Figure 1 A flow chart of a wireless channel simulation method based on time-slice differential trajectory shown in an embodiment of the present application;

[0035] Figure 2 A schematic diagram of a 3D scene from a base station to a vehicle and a time-sliced ​​differential trajectory parameter configuration for a vehicle networking scenario shown in an embodiment of the present application;

[0036] Figure 3 A schematic diagram of a time-varying Doppler curve of a channel coefficient after parameter configuration based on a time-slice differential trajectory shown in an embodiment of the present application;

[0037] Figure 4 A schematic diagram of an arc trajectory of a geostationary satellite motion suitable for a satellite communication scenario as shown in an embodiment of the present application;

[0038] Figure 5 A schematic diagram of the motion trajectory of a receiver suitable for a millimeter wave multi-user communication scenario disclosed in an embodiment of the present application. DETAILED DESCRIPTION

[0039] In order to make the objectives, technical solutions and advantages of the present application more clear, the embodiments of the present application will be further described in detail below with reference to the accompanying drawings.

[0040] When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Instead, they are only examples of methods and devices consistent with some aspects of the present application as detailed in the attached claims.

[0041] In the description of the present application, it should be understood that the terms "first", "second", etc. are used for descriptive purposes only and should not be understood as indicating or implying relative importance. For those of ordinary skill in the art, the specific meanings of the above terms in the present application can be understood according to the specific circumstances. In addition, in the description of the present application, unless otherwise specified, "multiple" refers to two or more. "And / or" describes the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. The character " / " generally indicates that the objects associated before and after are in an "or" relationship.

[0042] The present application proposes a wireless channel simulation method, device and storage medium based on time slice differential trajectory. By analyzing and decomposing a task, a security resource configuration model is established, and the model is solved using an improved marginal benefit algorithm, thereby achieving configuration optimization of product security resources.

[0043] The following is further described in conjunction with the embodiments and drawings.

[0044] Embodiment 1

[0045] The following will be combined Figure 1 , a wireless channel simulation method based on time slice differential trajectory provided in an embodiment of the present application is introduced in detail.

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

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

[0048] Step S2: Split the motion time of the receiver and the transmitter into multiple micro-segments, and configure the parameters of the motion trajectory of both ends of the receiver and the transmitter based on the micro-segment time and the starting coordinates of the receiver and the transmitter;

[0049] In the specific implementation process, the number, frequency and sampling rate, initial position, posture and antenna configuration of the transceivers (i.e., the receivers and transmitters hereinafter) are configured through the differential trajectory configurator. Different from the two mainstream configuration methods of the existing channel simulation configurator, that is, the regular patterning of the trajectory by adding straight line segments and circular arc segments, or the configuration method of inputting the position coordinates of the turning points of each segment, this embodiment adopts a time-based micro-segmentation method to split the complex motion trajectory of each transceiver into many small micro-segments, and determine the number of micro-segments S n , the total motion time T of the trajectory, and the motion time vector t within the differential segment time i ;in,

[0050]

[0051] Specifically, the configuration parameters include: the initial velocity scalar vr of the receiver starting point, the initial velocity scalar v of the transmitter starting point t ; Receiver micro-segment acceleration vector a r , transmitter micro-segment acceleration vector at; receiver motion direction vector micro-segment starting azimuth Transmitter motion direction vector micro-segment start azimuth Receiver movement towards vector micro segment start pitch angle Transmitter movement towards vector micro segment starting pitch angle Receiver motion heading vector micro-segment end azimuth Transmitter motion direction vector micro segment end azimuth Receiver motion towards vector micro segment end pitch angle Transmitter movement towards vector micro segment end pitch angle The receiver micro-segment attitude heading vector starts the Z-axis rotation angle Transmitter micro-segment attitude heading vector starting Z-axis rotation angle The attitude vector of the receiver micro-segment starts the Y-axis rotation angle The transmitter micro-segment attitude heading vector starts the Y-axis rotation angle The receiver micro-segment attitude heading vector starts the X-axis rotation angle Transmitter micro-segment attitude heading vector starting X-axis rotation angle The receiver micro-segment attitude heading vector ends the Z-axis rotation angle Transmitter micro-segment attitude heading vector end Z axis rotation angle The receiver micro-segment attitude heading vector ends the Y-axis rotation angle Transmitter micro-segment attitude heading vector ends Y-axis rotation angle The receiver micro-segment attitude heading vector ends the X-axis rotation angle Transmitter micro-segment attitude heading vector ends 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 According to the above configuration parameters, the total number of sample points of the channel coefficient 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 Nyquist sampling theorem;

[0055] Then the number of time domain sample points of the channel coefficient of the i-th micro-segment is P 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 a configuration panel, a software interface, etc. The configuration form includes but is not limited to software interface configuration, parameter table import, etc.

[0058] Step S3: based on the configured parameters, the initial speeds of the receiver and the transmitter, and each micro-segment time, the coordinates of each sample point, the receiver speed vector and the transmitter speed vector of each sample point are obtained;

[0059] First, after the configuration is executed, the vector coordinates [xr ,y r , z r ], transmitter vector coordinates [x t ,y t , z t ] and the receiver motion distance vector d r , the transmitter motion distance vector d t Perform calculations.

[0060] In this embodiment, the coordinates of the transceiver are calculated by a differential trajectory wireless channel calculator, wherein the coordinates of each sample point j in each differential segment i of the receiver are iteratively calculated as follows:

[0061]

[0062] where v ri is the initial speed of the receiver in each micro-segment time, which is determined by the initial speed v r1 , the acceleration and motion time of each micro-segment are iterated to obtain the initial velocity v of the second micro-segment r2 The calculation formula is:

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

[0064] v ri is the velocity vector of the receiver in each micro-segment time, which is composed of the unit vector of the motion direction of each micro-segment Initial velocity v in each micro-segment ri and the acceleration a in each micro-segment time ri , and the intra-segment movement time jΔt;

[0065] And so on:

[0066]

[0067] Sample point identification j=1,2,…,P in micro-segment time i , where θ ri (j) Yes arrive The monotonic linear difference within the interval, by default

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

[0069]

[0070] where v ti is the initial speed of the receiver in each micro-segment time, which is determined by the initial speed v t1It is obtained by iteration with the acceleration and motion time of each micro-segment;

[0071] The initial velocity v of the second microsegment t2 The calculation formula is:

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

[0073] where v ti is the velocity vector of the transmitter in each micro-segment time, and is the unit vector of the motion direction of each micro-segment Initial velocity v in each micro-segment ti and the acceleration a in each micro-segment time ti , and the movement time jΔt within the segment, so:

[0074]

[0075] Sample point identification j=1,2……P within micro-segment time i , where θ ti (j) Yes arrive The monotonic linear difference within the interval, by default yes arrive The monotonic linear difference within the interval, by default

[0076] The starting coordinates of the receiver and transmitter can be used to configure the motion trajectory of the two ends of the transceiver and draw the trajectory lines. In addition, for the convenience of configuration, there is a linkage operation relationship between the motion time and the motion distance in the micro-segment time through the micro-segment acceleration vector. The configuration can be performed by adjusting the motion time, and it is also compatible with the distance-based configuration method. That is, the transceiver motion distance is calculated, and the motion distance d of each micro-segment transmitter is ti The calculation formula is as follows:

[0077]

[0078] The moving distance d of each micro-segment receiver ri The calculation formula is as follows:

[0079]

[0080] During the calculation process, the moving distance d ri , d ti With movement time t i The values ​​can be converted to each other, and the distance calculation of the micro-segment start and end point coordinates can also be used to reversely calculate the movement time.

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

[0082] Specifically, the differential trajectory wireless channel calculator performs the wireless channel parameter calculation based on the differential trajectory parameters at both ends of the transceiver, including completing the departure angle vector of the main path of the full trajectory of each transceiver link by using the intra-segment interpolation and segment splicing method based on the configuration parameters of step S2 within each differential segment time. The arrival angle vector of the main path

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

[0084]

[0085] in, is the departure pitch angle of each sampling point j in the i-th segment, is the departure azimuth of each sampling point j in the i-th segment, is the arrival pitch angle of each sampling point j in the i-th segment, is the arrival azimuth of each sampling point j in the i-th segment.

[0086] Step S5: Based on the main path departure angle vector and the main path arrival angle vector, 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 the transmitter are obtained; specifically including:

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

[0088]

[0089] Among them, F r,θ,gcs , is the antenna pattern of the receiver end perpendicular to the horizontal polarization global coordinate system.

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

[0091]

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

[0093] Main radial direction vector They are:

[0094]

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

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

[0097] The rotation angles of the Z-axis, Y-axis, and X-axis of the j-th sample point in the i-th segment of the receiver are α r , β r , γ r The value of each sample point can be determined by the starting angle value within the i segment. and the final rotation angle value Perform linear interpolation to obtain;

[0098] The rotation angles of the Z-axis, Y-axis, and X-axis of the j-th sample point in the i-th section of the transmitter are α t , β t , γ t The value of each sample point can be determined by the starting angle value within the i segment. and the final rotation angle value Perform linear interpolation respectively.

[0099] The coordinate vector of the antenna array at the receiver end after adjustment based on the attitude change parameters of the receiver and transmitter, including the X-axis, Y-axis and Z-axis rotation parameters for:

[0100]

[0101] in Transform the antenna array coordinate vector before the receiver;

[0102] The coordinate vector of the antenna array at the transmitter end after adjustment based on the attitude change parameters of the receiver and transmitter, including the X-axis, Y-axis and Z-axis rotation parameters for:

[0103]

[0104] in Antenna element coordinate vector before transformation at the transmitter.

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

[0106] Specifically, the receiver velocity vector v of each sample point obtained in step S3 is 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] Finally, the calculation of the main path channel coefficient of the complex trajectory between the link from the receiver antenna u to the transmitter antenna s is realized by the following formula.

[0111]

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

[0113] In the specific calculation process, based on the acquired data and configuration parameters in S1 and S2, the calculation of the micro-segmented motion distance and the coordinates of the end point of the transceiver, as well as the drawing of the trajectory line are completed. It can be understood by those skilled in the art that the trajectory types drawn by this method include but are not limited to regular trajectory types, such as arcs, polygons, continuous broken lines 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 complex trajectory multipath channel coefficients between the transceiver links in step S5 is performed, wherein the calculation form includes but is not limited to software background calculation, code script calculation, etc.;

[0114] During the specific implementation process, in order to obtain a data format file that can be used for channel simulator playback, the data format of the complex trajectory multipath channel coefficient data between the transceiver links calculated by micro-segmentation is converted through the differential trajectory wireless channel data generator, including but not limited to the conversion of computer memory storage data to .npy or .bin files, .npy output files to .bin output files and other file formats. Finally, a data format file that can be used for channel simulator playback at the transceiver end under any trajectory is generated by segment splicing, and the size of the segment storage may depend on the memory of the running hardware device.

[0115] It can be understood that the differential trajectory wireless channel data generator performs micro-segment calculation data format conversion based on the channel coefficient calculation result obtained by the differential trajectory wireless channel calculator executing the above step S2; the conversion form includes direct conversion and generation by the software background, and conversion again through a data processing tool after generating a certain data form, etc.; finally, a channel file suitable for playback by a channel simulator is generated by segment splicing, and the size of the segment storage may depend on the memory of the running hardware device.

[0116] The application method disclosed in this embodiment directly completes the calculation of the configuration vector parameters related to the entire trajectory by splicing differential trajectories, and finally completes the wireless channel calculation under the complex trajectory by a segmented acceleration algorithm based on the continuous configuration parameters of the entire trajectory, and the size of the segmented storage can depend on the memory of the running hardware device. The problem of difficult configuration of the dot position parameters of complex trajectories in the existing wireless channel simulator is solved, and the efficiency of dynamic wireless channel modeling and channel data generation of complex trajectory receivers and transmitters is improved.

[0117] In order to further describe the content of this embodiment, three specific examples are listed below for illustration:

[0118] Example 1

[0119] In the scenario of the Internet of Vehicles, when a car communicates with a base station while driving, when the car turns at an elevated intersection, it drives from north to south, then changes lanes to enter the ramp, and then winds up from the ramp to drive from west to east. The specific technical solution is as follows:

[0120] Step 1: After completing the configuration of the conventional parameters such as the number of transceivers, frequency and sampling rate, initial position, attitude and antenna configuration through the differential trajectory configurator, such as Figure 2 As shown in the figure, the time-based micro-segmentation method is used to split the complex motion trajectory of the car receiver into many small micro-segments, and the motion time vector t of each micro-segment of the receiver is configured through the differential track configurator, in seconds; the motion orientation vector of the receiver micro-segment includes the starting azimuth Pitch Angle End Azimuth Pitch Angle The receiver micro-segment attitude heading vector includes the starting Z-axis rotation angle Y-axis rotation X-axis rotation End Z-axis rotation Y-axis rotation X-axis rotation The units of the above angles are all degrees; the receiver acceleration vector a r , in meters per second squared. It can be seen that since the differential trajectory configurator divides the irregular curve into many micro-segments, and the angles of each micro-segment are continuous, the final configured parameters meet the space-time consistency of the wireless channel. The main parameter configuration table is as follows:

[0121] Table 1 Time slice trajectory differential segmentation configuration parameter table

[0122]

[0123] Step 2: After the configuration is executed, the differential trajectory wireless channel calculator calculates the transceiver's differential segment motion distance vector D and the transceiver's endpoint vector coordinates [x, y, z] in meters. The calculation results can be found in Figure 2 , and finally realize the configuration of the motion trajectory of the automobile terminal and the drawing of the trajectory line, and the coordinate position of the base station remains unchanged. The differential trajectory wireless channel calculator performs the wireless channel parameter calculation of the differential trajectory parameters, including the configuration parameters of step 2 in each differential segment time, using the intra-segment interpolation and segment splicing method based on formulas (13) to (16) to complete the calculation of the main path departure angle vector AOD, ZOD, and the main path arrival angle vector AOA, ZOA of each transceiver link; based on formulas (21) to (22), the calculation of the receiver and transmitter direction vectors is completed. Based on formulas (17) to (20), the main path antenna gain vector parameters on the two polarizations of the receiver and transmitter are quickly calculated. Based on the formulas of the specific embodiment, all parameters of each sample point are calculated, 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 space-time consistency, no additional processing is required when the wireless channel coefficient calculation is realized through parameters, thereby speeding up the simulation time.

[0124] Step 3: Through the differential trajectory wireless channel data generator, the data format of the complex trajectory multipath channel coefficient data between the transceiver links calculated by the differential segment is converted, including but not limited to the conversion of computer memory storage data to .npy or .bin files, .npy files to .bin files, etc., to generate a data format file that can be used for channel simulator playback at the transceiver end under any trajectory. At the same time, the generator can generate data based on the channel and realize the verification of the main channel parameters through digital signal processing, such as Figure 3 As shown in the figure, the Doppler of the generated channel link is verified by short-time Fourier transform. Finally, the change trend of the time-varying Doppler frequency deviation is in line with expectations, and its change process satisfies continuous consistency.

[0125] Example 2

[0126] This example is for the communication between satellite and land terminal in satellite communication scenario, where the communication satellite is a geostationary orbit satellite, 35786.043 kilometers from the ground, and the projected position is 101.4 degrees east longitude and 0 degrees north latitude. The satellite terminal is at 22 degrees 17 minutes 37 seconds north latitude (22.2936 degrees) and 114 degrees 10 minutes 31 seconds east longitude (114.1753 degrees) on the earth. The specific technical solution is implemented as follows:

[0127] Step 1: After the number of satellites and ground terminals, frequency and sampling rate, initial position, attitude and antenna configuration and other conventional related parameters are configured through the differential trajectory configurator, the motion trajectories of the ground terminal and the satellite are divided into many small micro-segments by using the time-based micro-segmentation method, and the motion time vector t of the transceiver in each micro-segment time is configured through the differential trajectory configurator. r , t t ; The motion direction vectors of the ground terminal and the satellite for each micro-segment include the starting azimuth Pitch Angle End Azimuth Pitch Angle The attitude heading vectors of the ground terminal and the satellite for each micro-segment include the starting Z-axis rotation angle Y-axis rotation X-axis rotation End Z-axis rotation Y-axis rotation X-axis rotation The acceleration vector a of the ground terminal and the satellite in their respective micro-segment time r , a t .

[0128] Step 2: After the configuration is executed, the differential trajectory wireless calculator is used to complete the motion distance vector d of each differential segment between the ground terminal and the satellite. r , d t , the end point ground terminal and satellite vector coordinates [x r ,y r , z r ]、[x t ,y t , z t ] calculation to realize the configuration of the motion trajectory of the ground terminal and the satellite and the drawing of the trajectory lines. Figure 4A schematic diagram is drawn for the arc trajectory of the geostationary satellite after 31 minutes and 40 seconds of movement. The wireless channel parameter calculation based on the differential trajectory parameters at both ends of the transceiver is performed by the differential trajectory wireless channel calculator, including the mobile coordinate configuration parameters of step 2 in each differential segment time, and the calculation of the main path departure angle vector AOD, ZOD, and the main path arrival angle vector AOA, ZOA of the full trajectory of each transceiver link is completed based on formulas (13) to (16) by using intra-segment interpolation and segment splicing methods; the calculation of the receiver and transmitter direction vectors is completed based on formulas (21) to (22). The main path antenna gain vector parameters on both polarizations of the receiver and transmitter are quickly calculated based on formulas (17) to (20). All parameters of each sample point are calculated based on the various formulas of the specific embodiment, and finally the calculation of the complex trajectory main path channel coefficient between the transceiver links is realized through formula (27).

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

[0130] Example 3

[0131] This example is for the millimeter wave multi-user mobile communication scenario. The communication signal frequency is 26GHz, and the millimeter wave terminal is stationary. The millimeter wave terminal moves around the millimeter wave base station in an irregular manner. In this use case, the number of base stations is 1 and the number of terminals is 2. The specific technical solution is implemented as follows:

[0132] Step 1: After completing the configuration of conventional related parameters such as the number of transceivers, millimeter wave frequency and sampling rate, initial position, attitude and antenna configuration through the differential trajectory configurator, the complex motion trajectory of each receiver is divided into many small micro-segments by using a time-based micro-segmentation method, and the motion time vector t of each receiver in each micro-segment time is configured through the differential trajectory configurator, in seconds; the motion orientation vector of the receiver micro-segment includes the starting azimuth Pitch angle θ s , Ending azimuth Pitch angle θ e , the receiver micro-segment attitude heading vector includes the starting Z-axis rotation angle α s , Y-axis rotation angle β s , X-axis rotation angle γ s 、End Z-axis rotation angle α e , Y-axis rotation angle β e , X-axis rotation angle γ e , the units of the above angles are all degrees; the receiver acceleration vector a, the unit is meters per second squared.

[0133] Step 2: After the configuration is executed, the differential trajectory wireless channel calculator completes the calculation of the motion distance vector D of each micro-segment of each receiver and the coordinates [x, y, z] of the end point vector of each receiver segment, in meters. The final implementation is as follows Figure 5 The configuration of the receiver motion trajectory and the drawing of the trajectory line are shown. The coordinate position of the transmitter base station remains unchanged, and the motion trajectories experienced by different receivers in the air and on the ground are different, but the motion time is the same. It can be seen that since the differential trajectory configurator divides the irregular curve into many micro-segments, and the angles of each micro-segment are continuous, the parameters finally configured meet the space-time consistency of the wireless channel. The differential trajectory parameters are calculated by the differential trajectory wireless channel calculator, including configuring the parameters of the moving coordinates of step 2 within each micro-segment time, using segment interpolation and segment splicing to complete the calculation of the main path departure angle vector AOD, ZOD, and the main path arrival angle vector AOA, ZOA of each transceiver link based on formula (13) to formula (16); based on formula (21) to formula (22), the calculation of the receiver and transmitter direction vector is completed. Based on formula (17) to formula (20), the fast calculation of the main path antenna gain vector parameters on the two polarizations of the receiver and transmitter is realized. Based on the formulas of the specific embodiment, all parameters of each sample point are calculated, and finally the calculation of the main path channel coefficient of the complex trajectory between the transmitting and receiving links is realized through formula (27). Since the parameters configured by the differential trajectory configurator already have space-time consistency, no additional processing is required when the wireless channel coefficient is calculated through the parameters, thereby speeding up the simulation time. In addition, the parallel link processing method between multiple users can also speed up the simulation speed.

[0134] Step 3: Through the differential trajectory wireless channel data generator, the data format of the complex trajectory multipath channel coefficient data between the transceiver links calculated by micro-segmentation is converted, including but not limited to the conversion of computer memory storage data to .npy or .bin files, .npy files to .bin files, etc., and finally, through segmented splicing, a data format file that can be used for channel simulator playback for each transceiver link under any trajectory is generated.

[0135] Embodiment 2

[0136] The present application proposes a wireless channel simulation system based on time slice differential trajectory, comprising:

[0137] include:

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

[0139] A calculation module is used to combine the data in the acquisition module, and by splitting the movement time of the receiver and the transmitter into multiple micro-segments, calculate the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter, the distance between the two ends of the receiver and the transmitter, the main path direction vector of the receiver and the transmitter, the receiver velocity vector and the transmitter velocity vector of each sample point, the antenna gain vector of each sample point on the two polarizations of the main path transmitter end, and the antenna gain vector of each sample point on the two polarizations of the main path receiver end;

[0140] The determination module is used to determine the main path channel coefficient of the trajectory between the receiver and the transmitter link according to the calculation result of the calculation module.

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

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

[0143] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented 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] Embodiment 3

[0145] The present application proposes a computer storage medium on which a computer program is stored. When the program is executed by a processor, the steps of the wireless channel simulation method based on time-slice differential trajectory in any embodiment are implemented.

[0146] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0147] Therefore, the present application also proposes a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the method described in any embodiment of the present application is implemented. The computer-readable storage medium can be configured in any device of the present application.

[0148] Embodiment 4

[0149] The present application also proposes an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method described in any embodiment of the present application when executing the computer program.

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

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

[0152] As a computer-readable storage medium, the storage device can be used to store software programs, computer executable programs and module units, such as program instructions corresponding to the method in the embodiment of the present application. The storage device may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, at least one application required for a function; the storage data area may store data created according to the use of the terminal, etc. In addition, the storage device may include a high-speed random access memory, and may also include a non-volatile memory, such as at least one disk storage device, a flash memory device, or other non-volatile solid-state storage device. In some instances, the storage device may further include a memory remotely arranged relative to the processor, and these remote memories may be connected via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network and a combination thereof. The various embodiments in the present application are described in a progressive manner, and the same and similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments.

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

[0154] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. The embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interface, indirect coupling or communication connection of the device or unit, which can be electrical, mechanical or other forms.

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

[0156] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0157] It should be noted that the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the systems, methods and computer program products according to multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of a code, and the module, a program segment or a part of a code contains 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 box can also occur in a different order from the order marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of the boxes in the block diagram and / or flowchart can be implemented with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.

[0158] Finally, it should be noted that the above-described embodiments are only specific implementation methods of the present application, which are used to illustrate the technical solutions of the present application, rather than to limit them. The protection scope of the present application is not limited thereto. Although the present application is described in detail with reference to the above-mentioned embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-mentioned embodiments within the technical scope disclosed in the present application, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be included in the protection scope of the present application. Therefore, the protection scope of the present application shall be based on the protection scope of the claims.

[0159] The applicant of this application has made a detailed explanation and description of the implementation examples of this application in conjunction with the drawings in the specification. However, those skilled in the art should understand that the above implementation examples are only preferred implementation schemes of this application, and the detailed description is only to help readers better understand the spirit of this application, and it is not a limitation on the scope of protection of this application. On the contrary, any improvements or modifications 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 trajectory, characterized in that: include: Obtain the movement time of the receiver and transmitter, the starting coordinates of the receiver and transmitter, the initial speed of the receiver and transmitter, and the distance between the receiver and transmitter; The motion time of the receiver and the transmitter is divided into multiple micro-segments, and the parameters of the motion trajectory of the receiver and the transmitter are configured based on the micro-segment time and the starting coordinates of the receiver and the transmitter; Based on the configured parameters, the initial velocities of the receiver and the transmitter, and the time of each micro-segment, the coordinates of each sample point, the receiver velocity vector and the transmitter velocity vector of each sample point are obtained; Based on the coordinates of each sample point, the main path departure angle vector and the main path arrival angle vector of each sample point are obtained by using the intra-segment sample point calculation method; Based on the main path departure angle vector and the main path arrival angle vector, 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 the transmitter are obtained; Based on the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter, the distance between the receiver and the transmitter, the main path direction vectors of the receiver and the transmitter, the receiver velocity vector and the 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, the main path channel coefficient of the trajectory between the receiver and the transmitter link is determined.

2. The wireless channel simulation method based on time slice differential trajectory according to claim 1 is characterized in that: The parameters for completing the motion trajectory of the receiver and the transmitter include: The respective micro-segment times of the receiver and the transmitter, the respective micro-segment motion orientation vectors of the receiver and the transmitter, the respective micro-segment attitude orientation vectors of the receiver and the transmitter, and the respective micro-segment acceleration vectors of the receiver and the transmitter.

3. The wireless channel simulation method based on time slice differential trajectory according to claim 1 is characterized in that: The motion orientation vectors of the micro-segments of the receiver and the transmitter in the parameters for completing the motion trajectories at both ends of the receiver and the transmitter include: The starting azimuth angle of the receiver, the starting elevation angle of the receiver, the ending azimuth angle of the receiver, and the ending elevation angle of the receiver. The starting azimuth angle of the transmitter, the starting elevation angle of the transmitter, the ending azimuth angle of the transmitter, and the ending elevation angle of the transmitter.

4. The wireless channel simulation method based on time slice differential trajectory according to claim 1 is characterized in that: The parameters for completing the motion trajectories of the receiver and the transmitter include the micro-segmented attitude orientation vectors of the receiver and the transmitter, respectively: The rotation angle of the starting three-dimensional coordinate axis of the receiver, the rotation angle of the ending three-dimensional coordinate axis of the receiver, the rotation angle of the starting three-dimensional coordinate axis of the transmitter, and the rotation angle of the ending three-dimensional coordinate axis of the transmitter.

5. The wireless channel simulation method based on time slice differential trajectory according to claim 1 is characterized in that: The adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter are calculated by the antenna array coordinate vectors before transformation at both ends of the receiver and the 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 trajectory according to claim 1 is characterized in that: 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 required to calculate the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter are obtained by linear interpolation of the starting rotation angle values ​​and the ending rotation angle values ​​within the respective differential segment times of the receiver and the transmitter.

7. The wireless channel simulation method based on time slice differential trajectory according to claim 1 is characterized in that: The receiver velocity vector of each sample point is calculated by the initial velocity of the receiver within the micro-segment time, the movement time of the receiver within the segment, the acceleration vector of the receiver within the micro-segment time, and the unit vector of the receiver micro-segment movement direction; The transmitter velocity vector of each sample point is calculated by the transmitter's initial velocity within the micro-segment time, the transmitter's movement time within the segment, the transmitter's acceleration vector within the micro-segment time, and the transmitter's micro-segment movement direction unit vector.

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

9. A wireless channel simulation system based on time slice differential trajectory, characterized in that: include: An acquisition module is used to acquire the movement time of the receiver and the transmitter, the starting coordinates of the receiver and the transmitter, the initial speed of the receiver and the transmitter, and the distance between the two ends of the receiver and the transmitter; A calculation module is used to combine the data in the acquisition module, and by splitting the movement time of the receiver and the transmitter into multiple micro-segments, calculate the adjusted antenna array coordinate vectors at both ends of the receiver and the transmitter, the distance between the two ends of the receiver and the transmitter, the main path direction vector of the receiver and the transmitter, the receiver velocity vector and the transmitter velocity vector of each sample point, the antenna gain vector of each sample point on the two polarizations of the main path transmitter end, and the antenna gain vector of each sample point on the two polarizations of the main path receiver end; The determination module is used to determine the main path channel coefficient of the trajectory between the receiver and the transmitter link according to the calculation result of the calculation module.

10. A computer storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the wireless channel simulation method based on time-sliced ​​differential trajectory in any one of claims 1 to 8 are implemented.

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, the method according to any one of claims 1 to 8 is implemented.

Images