Disclosure of Invention
The invention provides a complicated mobile time-varying wireless channel simulation method based on ray tracing, which solves the problems that when the conventional wireless channel simulation method is applied to an Internet of things scene, the accuracy of simulation test is low, so that a user cannot accurately predict the communication range of wireless communication in the real Internet of things scene according to the result of the simulation test, and the wireless communication layout design cost of the real Internet of things scene is increased in a phase-changing manner.
In order to achieve the above purpose, the invention adopts the following technical scheme:
in a first aspect, the present invention provides a method for simulating a complex mobile time-varying wireless channel based on ray tracing, including:
Responding to the operation of configuring a virtual Internet of things environment by a user, acquiring target data of the virtual Internet of things environment, wherein the virtual Internet of things environment is an environment formed by a plurality of objects, the objects comprise a transmitting end and a receiving end communicated with the transmitting end, the objects comprise dynamic objects, the target data comprise environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual Internet of things environment, and the antenna parameters comprise the number of antennas of the transmitting end and the number of antennas of the receiving end;
determining the ray parameters of the direct rays, the ray parameters of each first reflected ray, the ray parameters of each first diffracted ray and the ray parameters of each mixed ray corresponding to each pair of receiving and transmitting end antennas at each sampling moment by adopting a preset ray tracking algorithm according to the environmental parameters, the antenna parameters, the reflection orders and the diffraction orders; the pair of receiving and transmitting end antennas comprise an antenna of the transmitting end and an antenna of the receiving end, the mixed rays are mixed rays of second reflected rays and second diffracted rays, and the ray parameters comprise paths and receiving power;
adopting a preset scattering model, and determining the ray parameters of each scattered ray corresponding to each pair of transceiver antennas at each sampling moment according to the environment parameters and the ray parameters of each first reflected ray corresponding to each pair of transceiver antennas at each sampling moment;
And determining a channel simulation result according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time, wherein the multipath components comprise direct rays, first reflected rays, first diffracted rays, mixed rays and scattered rays, the channel simulation result comprises Doppler frequency shift, the Doppler frequency shift is used for representing the influence of the movement of an object on the transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end.
In one possible implementation, the channel simulation result further includes a root mean square delay spread corresponding to each sampling instant, where the root mean square delay spread is used to represent the effect of multipath effects on the transmission characteristics of the channel.
In a possible implementation manner, the antenna parameters further include a gain of the transmitting end and a transmitting power of the transmitting end, the scattering model is a single scattering lobe model, and a preset scattering model is adopted to determine, according to the environmental parameters and a radiation parameter of each reflected radiation corresponding to each pair of transmitting-receiving end antennas at each sampling time, a radiation parameter of each scattered radiation corresponding to each pair of transmitting-receiving end antennas at each sampling time, including: determining each scattering surface corresponding to each pair of receiving-transmitting end antennas at each sampling time according to the path of each reflected ray corresponding to each pair of receiving-transmitting end antennas at each sampling time; determining a path of each corresponding scattered ray in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, wherein the first included angle is a maximum included angle between each corresponding scattered ray and a reflected ray in each scattering surface; and determining the receiving power of each corresponding scattered ray in each scattering surface at each sampling moment according to the path of each scattered ray included in each scattering surface at each sampling moment, the gain of the transmitting end, the transmitting power of the transmitting end and the scattering coefficient corresponding to each scattering surface.
In one possible implementation manner, determining a path of each scattered ray corresponding to each scattering surface at each sampling moment according to a preset scattering lobe width factor, a number of scattered rays and a first included angle of each scattering surface at each sampling moment includes: determining each scattering point included in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment; and determining the path of each scattered ray corresponding to each scattering surface at each sampling moment according to each scattering point included in each scattering surface at each sampling moment, and the transmitting end and the receiving end at the sampling moment.
In one possible implementation manner, determining each scattering point included in each scattering surface at each sampling time according to a preset scattering lobe width factor, a number of scattering rays and a first included angle of each scattering surface at each sampling time includes: according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, determining a second included angle corresponding to each scattered ray in each scattering surface at each sampling moment, wherein the second included angle is an included angle between the scattered ray and the reflected ray, and the second included angle is smaller than the first included angle; and determining each scattering point included in each scattering surface at each sampling moment according to the second included angle corresponding to each scattering ray included in each scattering surface at each sampling moment.
In one possible implementation manner, at least one of the transmitting end and the receiving end is a dynamic object, the antenna parameters further include frequencies of simulated radio waves emitted by the transmitting end, the environmental parameters include a motion path and a motion speed of the transmitting end or a motion path and a motion speed of the receiving end, and the determining a channel simulation result according to a ray parameter of each multipath component corresponding to each pair of transmitting and receiving end antennas at each sampling moment includes: determining the relative speed of the transmitting end and the receiving end at each sampling moment according to the motion path and the motion speed of the transmitting end or the motion path and the motion speed of the receiving end; determining normalized channel capacity according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining Doppler frequency shift corresponding to each sampling time according to the path, the frequency and the relative speed of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time.
In one possible implementation manner, determining a channel simulation result according to a ray parameter of each multipath component corresponding to each pair of transceiver antennas at each sampling time includes: determining the total received power corresponding to each sampling moment according to the received power of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; determining the time delay of each multipath component corresponding to each sampling moment according to the path of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining the root mean square delay spread corresponding to each sampling time according to the total received power corresponding to each sampling time and the delay and the received power of each multipath component corresponding to each sampling time.
In one possible implementation manner, in response to an operation of configuring the virtual internet of things environment by a user, obtaining target data of the virtual internet of things environment includes: determining environmental parameters of an environment in which the channel is located in response to an environmental layout operation of a user, the environmental parameters including a shape, a position, a material type of each static object included in the environment, and a shape, an initial position, a material type, a motion path, and a motion speed of each dynamic object; determining antenna parameters in response to antenna configuration operation of a user, wherein the antenna parameters comprise frequency of simulated radio waves sent by a transmitting end, gain of the transmitting end, transmitting power of the transmitting end, moving path and moving speed of a receiving end, antenna types and the number of antennas of the transmitting end, and antenna types and the number of antennas of the receiving end; determining a reflection order and a diffraction order in response to a ray tracing configuration operation of a user; each sampling instant is determined in response to a sampling configuration operation by a user.
In a second aspect, the present invention provides a complex mobile time-varying wireless channel simulation device based on ray tracing, which is characterized by comprising:
The system comprises an acquisition module, a sampling module and a processing module, wherein the acquisition module is used for responding to the operation of configuring a virtual Internet of things environment by a user, the virtual Internet of things environment is an environment formed by a plurality of objects, the objects comprise a transmitting end and a receiving end communicated with the transmitting end, the objects comprise dynamic objects, the target data comprise environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual Internet of things environment, and the antenna parameters comprise the number of antennas of the transmitting end and the number of antennas of the receiving end;
the first determining module is used for determining the ray parameters of the direct rays, the ray parameters of the first reflected rays, the ray parameters of the first diffracted rays and the ray parameters of the mixed rays corresponding to each pair of receiving and transmitting end antennas at each sampling moment according to the environmental parameters, the antenna parameters, the reflection orders and the diffraction orders by adopting a preset ray tracking algorithm; the pair of receiving and transmitting end antennas comprise an antenna of a transmitting end and an antenna of a receiving end, the mixed rays are mixed rays of second reflected rays and second diffracted rays, and the ray parameters comprise paths and receiving power;
the second determining module is used for determining the ray parameters of each scattered ray corresponding to each pair of transceiver antennas at each sampling moment according to the environment parameters and the ray parameters of each first reflected ray corresponding to each pair of transceiver antennas at each sampling moment by adopting a preset scattering model;
The third determining module is configured to determine a channel simulation result according to a ray parameter of each multipath component corresponding to each pair of transceiver antennas at each sampling time, where the multipath components include a direct ray, a first reflected ray, a first diffracted ray, a mixed ray, and a scattered ray, the channel simulation result includes a doppler shift, the doppler shift is used to represent an effect of movement of an object on a transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end.
In one possible implementation, the channel simulation result further includes a root mean square delay spread corresponding to each sampling instant, where the root mean square delay spread is used to represent the effect of multipath effects on the transmission characteristics of the channel.
In one possible implementation manner, the antenna parameter further includes a gain of the transmitting end and a transmitting power of the transmitting end, the scattering model is a single scattering lobe model, and the second determining module is specifically configured to: determining each scattering surface corresponding to each pair of receiving-transmitting end antennas at each sampling time according to the path of each reflected ray corresponding to each pair of receiving-transmitting end antennas at each sampling time; determining a path of each corresponding scattered ray in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, wherein the first included angle is a maximum included angle between each corresponding scattered ray and a reflected ray in each scattering surface; and determining the receiving power of each corresponding scattered ray in each scattering surface at each sampling moment according to the path of each scattered ray included in each scattering surface at each sampling moment, the gain of the transmitting end, the transmitting power of the transmitting end and the scattering coefficient corresponding to each scattering surface.
In one possible implementation manner, the second determining module is specifically configured to: determining each scattering point included in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment; and determining the path of each scattered ray corresponding to each scattering surface at each sampling moment according to each scattering point included in each scattering surface at each sampling moment, and the transmitting end and the receiving end at the sampling moment.
In one possible implementation manner, the second determining module is specifically configured to: according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, determining a second included angle corresponding to each scattered ray in each scattering surface at each sampling moment, wherein the second included angle is an included angle between the scattered ray and the reflected ray, and the second included angle is smaller than the first included angle; and determining each scattering point included in each scattering surface at each sampling moment according to the second included angle corresponding to each scattering ray included in each scattering surface at each sampling moment.
In one possible implementation manner, at least one of the transmitting end and the receiving end is a dynamic object, the antenna parameter further includes a frequency of a simulated radio wave emitted by the transmitting end, the environmental parameter includes a motion path and a motion speed of the transmitting end or a motion path and a motion speed of the receiving end, and the third determining module is specifically configured to: determining the relative speed of the transmitting end and the receiving end at each sampling moment according to the motion path and the motion speed of the transmitting end or the motion path and the motion speed of the receiving end; determining normalized channel capacity according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining Doppler frequency shift corresponding to each sampling time according to the path, the frequency and the relative speed of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time.
In one possible implementation manner, the third determining module is specifically configured to: determining the total received power corresponding to each sampling moment according to the received power of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; determining the time delay of each multipath component corresponding to each sampling moment according to the path of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining the root mean square delay spread corresponding to each sampling time according to the total received power corresponding to each sampling time and the delay and the received power of each multipath component corresponding to each sampling time.
In one possible implementation manner, the acquiring module is specifically configured to: determining environmental parameters of an environment in which the channel is located in response to an environmental layout operation of a user, the environmental parameters including a shape, a position, a material type of each static object included in the environment, and a shape, an initial position, a material type, a motion path, and a motion speed of each dynamic object; determining antenna parameters in response to antenna configuration operation of a user, wherein the antenna parameters comprise frequency of simulated radio waves sent by a transmitting end, gain of the transmitting end, transmitting power of the transmitting end, moving path and moving speed of a receiving end, antenna types and the number of antennas of the transmitting end, and antenna types and the number of antennas of the receiving end; determining a reflection order and a diffraction order in response to a ray tracing configuration operation of a user; each sampling instant is determined in response to a sampling configuration operation by a user.
In a third aspect, the present invention provides a computer device comprising: a processor and a memory; the memory is for storing computer program code, the computer program code comprising computer instructions; when the processor executes the computer instructions, the computer device performs the complex mobile time-varying radio channel emulation method based on ray tracing as described in the first aspect and any one of its possible implementations.
In a fourth aspect, the present invention provides a computer readable storage medium storing computer instructions that, when run on a computer device, cause the computer device to perform the complex mobile time-varying radio channel simulation method based on ray tracing as described in the first aspect and any one of its possible implementations.
The invention provides a complicated mobile time-varying wireless channel simulation method based on ray tracing, which comprises the following steps: responding to the operation of configuring a virtual Internet of things environment by a user, acquiring target data of the virtual Internet of things environment, wherein the virtual Internet of things environment is an environment formed by a plurality of objects, the objects comprise a transmitting end and a receiving end communicated with the transmitting end, the objects comprise dynamic objects, the target data comprise environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual Internet of things environment, and the antenna parameters comprise the number of antennas of the transmitting end and the number of antennas of the receiving end; determining the ray parameters of the direct rays, the ray parameters of each first reflected ray, the ray parameters of each first diffracted ray and the ray parameters of each mixed ray corresponding to each pair of receiving and transmitting end antennas at each sampling moment according to the environmental parameters, the antenna parameters, the reflection orders and the diffraction orders by adopting a preset ray tracking algorithm, wherein the ray parameters comprise paths and receiving power; the pair of receiving and transmitting end antennas comprise an antenna of the transmitting end and an antenna of the receiving end, and the mixed rays are mixed rays of second reflected rays and second diffracted rays; adopting a preset scattering model, and determining the ray parameters of each scattered ray corresponding to each pair of transceiver antennas at each sampling moment according to the environment parameters and the ray parameters of each first reflected ray corresponding to each pair of transceiver antennas at each sampling moment; and determining a channel simulation result according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time, wherein the multipath components comprise direct rays, first reflected rays, first diffracted rays, mixed rays and scattered rays, the channel simulation result comprises Doppler frequency shift, the Doppler frequency shift is used for representing the influence of the movement of an object on the transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end. The invention realizes the simulation support and adaptation aiming at the characteristics of a mobile time-varying channel, multipath components and the like in the virtual Internet of things scene based on the configured virtual Internet of things scene comprising the complex moving objects, has higher simulation test accuracy, and ensures that a user can accurately analyze and predict the communication range of wireless communication in the real Internet of things scene according to the simulation test result, thereby reducing the wireless communication layout design cost of the real Internet of things scene.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The terms "first" and "second" are used below for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the embodiments of the present disclosure, unless otherwise indicated, the meaning of "a plurality" is two or more. In addition, the use of "based on" or "according to" is intended to be open and inclusive in that a process, step, calculation, or other action "based on" or "according to" one or more of the stated conditions or values may in practice be based on additional conditions or beyond the stated values.
The sixth generation mobile communication technology is characterized by wider coverage, higher efficiency, more intelligent application and higher safety performance, and some novel application scenes appear along with the development of the technology. Internet of things scenes, particularly industrial internet of things (Industrial Internet ofThings, IIoT) scenes, are attracting attention as typical scenes therein.
At present, the method for analyzing and evaluating the scene of the industrial Internet of things mainly comprises channel measurement and numerical simulation. The ray tracing method is applied in some papers and verified by the actual channel measurement result to better explain the channel characteristics. In the prior art, channel simulation is performed on path attenuation parameters of an industrial Internet of things channel under a fifth generation mobile communication technology, and coverage of a wireless communication system is obtained. For the sixth generation mobile communication technology, the parameters of the channel need to be re-evaluated due to the higher frequency band used.
However, the industrial internet of things scene has mobile timeliness, the channel composition contains strong multipath, and the current channel simulation software lacks targeted support for the industrial internet of things scene. Taking the CloudRT software as an example, it lacks a motion pattern generation function, so it is difficult to configure a more complex moving scene, and it lacks a speed parameter, resulting in limited analysis of the transmission characteristics of the channel.
In addition, in the existing wireless channel simulation method, the generated scattered rays are based on the assumption that scattering is weak, the difference between the scattered rays and the geometrical propagation paths of the reflected rays depending on the scattered rays is ignored, an offset is added on the basis of the reflected rays only when the ray parameters of the scattered rays are calculated, only the situation that scattering is not strong can be analyzed, multipath effects caused by scattering cannot be well analyzed, simulation test accuracy is low, a user cannot accurately predict the communication range of wireless communication in a real Internet of things scene according to the result of the simulation test, and accordingly the problem of wireless communication layout design cost of the real Internet of things scene is improved in a variable manner.
In order to solve the problem that when the conventional wireless channel simulation method is applied to an Internet of things scene, the accuracy of simulation test is low, so that a user cannot accurately predict the communication range of wireless communication in the real Internet of things scene according to the result of the simulation test, and accordingly the design cost of the wireless communication layout of the real Internet of things scene is increased in a phase-changing manner, the embodiment of the invention provides a complicated mobile time-varying wireless channel simulation method based on ray tracing, which comprises the following steps: responding to the operation of configuring a virtual Internet of things environment by a user, acquiring target data of the virtual Internet of things environment, wherein the virtual Internet of things environment is an environment formed by a plurality of objects, the objects comprise a transmitting end and a receiving end communicated with the transmitting end, the objects comprise dynamic objects, the target data comprise environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual Internet of things environment, and the antenna parameters comprise the number of antennas of the transmitting end and the number of antennas of the receiving end; according to the environment parameters and the antenna parameters, determining the ray parameters of the direct rays corresponding to each sampling time, wherein the ray parameters comprise paths and received power; determining the ray parameters of each first reflected ray, the ray parameters of each first diffracted ray and the ray parameters of each mixed ray corresponding to each pair of receiving and transmitting end antennas at each sampling moment by adopting a preset ray tracking algorithm according to the environmental parameters, the antenna parameters, the reflection orders and the diffraction orders; the pair of receiving and transmitting end antennas comprise an antenna of the transmitting end and an antenna of the receiving end, and the mixed rays are mixed rays of second reflected rays and second diffracted rays; adopting a preset scattering model, and determining the ray parameters of each scattered ray corresponding to each pair of transceiver antennas at each sampling moment according to the environment parameters and the ray parameters of each first reflected ray corresponding to each pair of transceiver antennas at each sampling moment; and determining a channel simulation result according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time, wherein the multipath components comprise direct rays, first reflected rays, first diffracted rays, mixed rays and scattered rays, the channel simulation result comprises Doppler frequency shift, the Doppler frequency shift is used for representing the influence of the movement of an object on the transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end. The embodiment of the invention realizes the simulation support and adaptation aiming at the characteristics of a mobile time-varying channel, multipath components and the like in the virtual Internet of things scene based on the configured virtual Internet of things scene comprising the complex moving object, has higher simulation test accuracy, and ensures that a user can accurately analyze and predict the communication range of wireless communication in the real Internet of things scene according to the simulation test result, thereby reducing the wireless communication layout design cost of the real Internet of things scene.
The execution subject of the complex mobile time-varying wireless channel simulation method based on ray tracing provided by the embodiment of the invention is computer equipment. The computer device may be a terminal device, a server, or a server cluster. The embodiment of the invention is not limited.
Fig. 1 is a schematic flow chart of a complicated mobile time-varying wireless channel simulation method based on ray tracing according to an embodiment of the present invention. As shown in fig. 1, the complex mobile time-varying channel simulation method based on ray tracing may include the following steps S101 to S105.
S101, a computer device responds to operation of configuring a virtual Internet of things environment by a user to acquire target data of the virtual Internet of things environment, wherein the virtual Internet of things environment is an environment formed by a plurality of objects, the objects comprise a transmitting end and a receiving end communicated with the transmitting end, the objects comprise dynamic objects, the target data comprise environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual Internet of things environment, and the antenna parameters comprise the number of antennas of the transmitting end and the number of antennas of the receiving end.
It can be understood that the virtual internet of things environment is configured, that is, the process of building a virtual internet of things scene of wireless channel simulation, where the virtual internet of things scene can be a warehouse, and the virtual internet of things scene includes a plurality of static objects and dynamic objects. The target data includes, but is not limited to, environmental parameters of the virtual internet of things environment, antenna parameters, reflection orders, diffraction orders, sampling moments, and the antenna parameters include the number of antennas of a transmitting end, the number of antennas of a receiving end, and the like.
In one possible implementation manner, the obtaining, by the computer device, target data of the virtual internet of things environment in response to an operation of configuring the virtual internet of things environment by the user may include: the computer equipment responds to the environment layout operation of the user, and determines environment parameters of the environment in which the channel is positioned, wherein the environment parameters comprise the shape, the position and the material type of each static object included in the environment, and the shape, the initial position, the material type, the motion path and the motion speed of each dynamic object; the computer equipment responds to the antenna configuration operation of a user, and determines antenna parameters, wherein the antenna parameters comprise the frequency of simulated radio waves sent by a transmitting end, the gain of the transmitting end, the transmitting power of the transmitting end, the moving path and the moving speed of a receiving end, the antenna types and the number of antennas of the transmitting end, and the antenna types and the number of antennas of the receiving end; the computer device determining a reflection order and a diffraction order in response to a ray tracing configuration operation by a user; the computer device determines each sampling instant in response to a sampling configuration operation by a user.
Fig. 2 is a schematic view of a virtual intelligent warehouse scene constructed in the present embodiment, and as shown in fig. 2, the virtual intelligent warehouse scene is an indoor scene, and is 200 meters long, 100 meters wide, and 10-15 meters high. There are various types of cube-type wood cargo distributions, and there are concrete columns. The warehouse frame is made of metal, and the ground is made of concrete.
For environmental parameters, the embodiment draws the shape of a static object (such as a wall, the ground, static machine equipment and the like) and determines the position of the static object, and for moving objects, the moving parameters such as the path, the speed and the like of the moving object are respectively configured. In the present embodiment, the cart 1 is set to move from (61, 130,0.8) to (61, 105,0.8), and the cart 2 travels along a folding line path of (35, 140,0.8), (60, 140,0.8), (60, 120,0.8). Their rates were 0.83 m/s and 0.9 m/s, respectively, for 30 seconds and 50 seconds, respectively. The transmitting antenna is located at the position of the coordinates (38, 138,8), and the initial position of the receiving antenna is (65, 50, 1). The sender remains stationary and the receiver starts at (65, 50, 1) and then moves 100 meters in the positive y-axis direction at a speed of 1 meter/second. Wherein, the transmitting end is a base station antenna, and the receiving end is an intelligent trolley. The mobile intelligent vehicle receives the task information through the base station antenna or cooperates with other vehicles.
For antenna parameters, in this embodiment, the receiving end antenna and the transmitting end antenna are all omni-directional antennas, and are all antenna arrays. The receiving end antennas and the transmitting end antennas are uniform linear arrays of 1 row and 8 columns, each vertical direction comprises 1 antenna, each horizontal direction comprises 8 antennas, and the antenna spacing in the horizontal direction is 0.05 meter.
For the reflection order and the diffraction order, the reflection order is set to 1 order in the present embodiment, and diffraction is not considered.
For the sampling time, each sampling time interval is set to 10 seconds in this embodiment.
For the frequency of the artificial radio wave emitted from the transmitting end, the frequency is set to 3GHz in this embodiment.
In addition, according to the type of the material contained in the object, the dielectric constant of the material at the corresponding center frequency is introduced into the material library, and accordingly, the dielectric constant parameters of the materials such as metal, concrete, wood and the like at 3GHz are introduced into the material library.
S102, a computer device adopts a preset ray tracing algorithm, and determines the ray parameter of each first reflected ray, the ray parameter of each first diffracted ray and the ray parameter of each mixed ray corresponding to each pair of transceiver antennas at each sampling moment according to environmental parameters, antenna parameters, reflection orders and diffraction orders; the pair of receiving and transmitting end antennas comprise an antenna of the transmitting end and an antenna of the receiving end, the mixed rays are mixed rays of second reflected rays and second diffracted rays, and the ray parameters comprise paths and received power.
It will be appreciated that since a scene includes a plurality of moving objects, the position of each object at each sampling instant is different. In this embodiment, on the premise that the sampling time interval is 10 seconds, the position of each object at each sampling time point is shown in table 1.
TABLE 1
Sampling time t(s)
|
Transmitting antenna coordinates
|
Receiving antenna coordinates
|
Trolley 1 coordinate
|
Trolley 2 coordinates
|
0
|
(38,138,8)
|
(65,50,1)
|
(61,130,0.8)
|
(35,140,0.8)
|
10
|
(38,138,8)
|
(65,60,1)
|
(61,121.7,0.8)
|
(44,140,0.8)
|
20
|
(38,138,8)
|
(65,70,1)
|
(61,113.3,0.8)
|
(53,140,0.8)
|
30
|
(38,138,8)
|
(65,80,1)
|
(61,105,0.8)
|
(53,138,0.8)
|
40
|
(38,138,8)
|
(65,90,1)
|
(61,105,0.8)
|
(53,129,0.8)
|
50
|
(38,138,8)
|
(65,100,1)
|
(61,105,0.8)
|
(53,120,0.8)
|
60
|
(38,138,8)
|
(65,110,1)
|
(61,105,0.8)
|
(53,120,0.8)
|
70
|
(38,138,8)
|
(65,120,1)
|
(61,105,0.8)
|
(53,120,0.8)
|
80
|
(38,138,8)
|
(65,130,1)
|
(61,105,0.8)
|
(53,120,0.8)
|
90
|
(38,138,8)
|
(65,140,1)
|
(61,105,0.8)
|
(53,120,0.8)
|
100
|
(38,138,8)
|
(65,150,1)
|
(61,105,0.8)
|
(53,120,0.8) |
In addition, the second reflected ray and the second diffracted ray forming the mixed ray may not belong to the ray reaching the receiving end, and the first reflected ray and the first diffracted ray are both rays reaching the receiving end, unlike the second reflected ray and the second diffracted ray.
The preset ray tracing algorithm may be a mirror image method.
Specifically, the computer device may determine, according to the environmental parameter, the antenna parameter, the reflection order and the diffraction order, a radiation parameter of each first reflected radiation, a radiation parameter of each first diffracted radiation, and a radiation parameter of each mixed radiation, which correspond to each pair of transceiver antennas at each sampling time, by using a preset radiation tracking algorithm; the pair of receiving and transmitting end antennas comprise an antenna of the transmitting end and an antenna of the receiving end, and the mixed rays are mixed rays of second reflected rays and second diffracted rays.
S103, the computer equipment adopts a preset scattering model, and determines the radiation parameters of each scattered radiation corresponding to each pair of transceiver antennas at each sampling moment according to the environment parameters and the radiation parameters of each first reflected radiation corresponding to each pair of transceiver antennas at each sampling moment.
The preset scattering model may be a quasi-deterministic model (Q-D model) or a single scattering valve model (direct model).
Specifically, the computer device may use a preset scattering model, and determine, according to the environmental parameter and the radiation parameter of each first reflected radiation corresponding to each pair of transceiver antennas at each sampling time, the radiation parameter of each scattered radiation corresponding to each pair of transceiver antennas at each sampling time.
S104, the computer equipment determines a channel simulation result according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time, wherein the multipath components comprise direct rays, first reflected rays, first diffracted rays, mixed rays and scattered rays, the channel simulation result comprises Doppler frequency shift, the Doppler frequency shift is used for representing the influence of the movement of an object on the transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end.
The multipath components refer to electromagnetic wave components with different arrival times after the electromagnetic wave propagates through different paths. In this embodiment, the multipath components may include the direct radiation, the first reflected radiation, the first diffracted radiation, the mixed radiation, and the scattered radiation described above.
Specifically, the computer device may determine a channel simulation result according to a ray parameter of each multipath component corresponding to each pair of antennas at each sampling time, where the channel simulation result includes a doppler shift, where the doppler shift is used to represent an influence of motion of an object on a transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end.
In one possible implementation, the channel simulation results may further include a root mean square delay spread corresponding to each sampling instant, where the root mean square delay spread is used to represent the effect of multipath effects on the transmission characteristics of the channel.
In the embodiment, the computer equipment responds to the configuration operation of the user, builds up a virtual internet of things scene comprising a plurality of complex moving objects, realizes simulation support and adaptation aiming at characteristics such as a mobile time-varying channel, multipath components and the like in the virtual internet of things scene, has higher simulation test accuracy, and enables the user to accurately analyze and predict the communication range of wireless communication in the real internet of things scene according to the simulation test result, thereby reducing the wireless communication layout design cost of the real internet of things scene.
Optionally, the antenna parameters further include a gain of the transmitting end and a transmitting power of the transmitting end, the scattering model is a single scattering lobe model, fig. 1 and fig. 3 show a second schematic flow diagram of a complex mobile time-varying wireless channel simulation method based on ray tracing, and as shown in fig. 3, the step S103 may specifically include the following steps S301 to S303.
S301, the computer equipment determines each scattering surface corresponding to each pair of transceiver antennas at each sampling moment according to the path of each reflected ray corresponding to each pair of transceiver antennas at each sampling moment.
S302, the computer equipment determines a path of each corresponding scattered ray in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, wherein the first included angle is a maximum included angle between each corresponding scattered ray and a reflected ray in each scattering surface.
It can be understood that the preset scattering lobe width factor, the number of scattering rays and the first included angle describe a scattering phenomenon that scattering power is intensively distributed near the scattering surface reflection direction and the maximum value is obtained in the scattering surface reflection direction, so that the distribution of scattering field intensity or scattering power can be more accurately described. Wherein, the larger the value of the scattering lobe width factor, the narrower the scattering lobe width, and the better the direction concentration of the path of the scattered ray.
Specifically, the computer device may determine, according to a preset scatter lobe width factor, a number of scattered rays, and a first included angle of each scattering plane at each sampling time, a path of each corresponding scattered ray in each scattering plane at each sampling time, where the first included angle is a maximum included angle between each corresponding scattered ray and a reflected ray in each scattering plane at each sampling time.
In one possible implementation manner, the computer device may determine, according to a preset scatter lobe width factor, a number of scattered rays, and a first included angle of each scattering surface at each sampling time, a second included angle corresponding to each scattered ray included in each scattering surface at each sampling time, where the second included angle is an included angle between a scattered ray and a reflected ray, the second included angle is smaller than the first included angle, determine, according to a second included angle corresponding to each scattered ray included in each scattering surface at each sampling time, each scattering point included in each scattering surface at each sampling time, and determine, according to each scattering point included in each scattering surface at each sampling time, and a transmitting end and a receiving end of the sampling time, a path of each scattered ray corresponding to each scattering surface at each sampling time.
Illustratively, in this embodiment, a diffusion lobe width factor α is set R The number of scattered rays generated by each scattering surface is N=30, the first included angle is
The second included angle ψ satisfies:
Ψ~U(0,Ψ max ) (1);
for each scattered ray generated, the computer device may generate a curve of values of ψ on the scattering plane from the generated included angles ψ. A point is randomly generated on the curve as a scattering point and scattered radiation is generated from the determined scattering point.
S303, the computer equipment determines the receiving power of each scattering ray corresponding to each scattering surface at each sampling moment according to the path of each scattering ray included in each scattering surface at each sampling moment, the gain of the transmitting end, the transmitting power of the transmitting end and the scattering coefficient corresponding to each scattering surface.
It will be appreciated that the scattering coefficient corresponds to the material of the scattering surface.
Specifically, the computer device may determine, according to the path of each scattered ray included in each scattering surface at each sampling time, and the gain of the transmitting end, the transmitting power of the transmitting end, and the scattering coefficient corresponding to each scattering surface, the receiving power of each scattered ray corresponding to each scattering surface at each sampling time.
Exemplary, the reflected ray field strength is set as E in the present embodiment 0 A central scattered ray field strength E S0 Scattered ray field strength |E S The I can be expressed as:
wherein, the liquid crystal display device comprises a liquid crystal display device,
/>
wherein G is t Representing the gain of the transmitting end, P t Representing the transmitting power of the transmitting end, S is the scattering coefficient of the material of the scattering surface, r i R is the distance from the transmitting antenna to the incident point of the surface of the material s For receiving the distance from the antenna to the incident point on the surface of the material, θ i Is the included angle theta between the incident path of the ray and the normal line s dS represents the scattering source per unit area of the surface of the material, which is the angle between the ray scattering path and the normal.
Of course, the accurate qualitative model may also be used to determine the radiation parameters of the scattered radiation, the specific principles of which are not described in detail herein. Illustratively, in the present embodiment, at the initial time of t=0, the path of the multipath component calculated using the accurate qualitative model is shown in fig. 4, and the path of the multipath component calculated using the single scattering lobe model is shown in fig. 5.
In this embodiment, the single scattering lobe model considers the difference between the geometric propagation paths of the scattered rays and the reflected rays, and compared with the accurate qualitative model, the single scattering lobe model not only supports the simulation under the stronger scattering condition, but also can better consider the multipath effect. In addition, the single scattering valve model also supports the independent configuration of the scattering coefficient of each scattering surface, and the accuracy of channel simulation is further improved.
Optionally, at least one of the transmitting end and the receiving end is a dynamic object, the antenna parameter further includes a frequency of a simulated radio wave emitted by the transmitting end, and the environmental parameter includes a motion path and a motion speed of the transmitting end or a motion path and a motion speed of the receiving end. Based on fig. 3, fig. 6 shows a third flow chart of a complicated mobile time-varying wireless channel simulation method based on ray tracing, and as shown in fig. 6, the above step S104 may specifically include the following steps S601 to S602.
S601, the computer equipment determines the relative speed of the transmitting end and the receiving end at each sampling moment according to the motion path and the motion speed of the transmitting end or the motion path and the motion speed of the receiving end.
S602, the computer equipment determines Doppler frequency shift corresponding to each sampling time according to the path, the frequency and the relative speed of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time.
Specifically, the computer device may determine the doppler shift corresponding to each sampling time according to the path, the frequency, and the relative speed of each multipath component corresponding to each pair of transceiver antennas at each sampling time.
In the present embodiment, doppler shift f D The expression of (2) is:
wherein c is the speed of light, f is the frequency, v is the relative speed, θ is the complement angle of the included angle between the direction of motion of the receiving end and the direction of rays of the multipath component. The maximum doppler shift versus time is shown in fig. 7 and the minimum doppler shift versus time is shown in fig. 8.
Optionally, the computer device further determines the normalized channel capacity according to a ray parameter of each multipath component corresponding to each pair of transceiver antennas at each sampling instant.
Specifically, the computer device may determine a channel transmission matrix according to the ray parameters of each multipath component corresponding to each pair of transceiver antennas at each sampling time, and determine the normalized channel capacity according to the channel transmission matrix and a preset signal-to-noise ratio.
In this embodiment, let the channel transmission matrix be H, element H in H ij The expression is:
wherein h is ij Representing the transfer function of the channel between the antenna of the ith transmitting end and the antenna of the jth receiving end, N R Representing the number of multipath components between the pair of transceivers, P k 、τ k And theta k The power, delay and phase shift of the kth multipath component between the pair of transceivers are represented, respectively.
Let the normalized channel capacity be C, its expression be:
Wherein SNR is a preset signal-to-noise ratio, I NR Is N R ×N R Identity matrix, N Tx Indicating the number of antennas at the transmitting end. Fig. 9 is a schematic diagram showing a relationship between a preset signal-to-noise ratio and a normalized channel capacity under different transceiver antenna configurations, as shown in fig. 9, in which, at an initial time of t=0, the transceiver antenna configurations are changed, and channel capacities of multiple-input multiple-output (MIMO), multiple-input single-output (MISO), single-input multiple-output (SIMO) and single-input single-output (SISO) are compared.
Optionally, based on fig. 6, fig. 10 shows a fourth flow chart of a complicated mobile time-varying wireless channel simulation method based on ray tracing, as shown in fig. 10, step S104 may further include the following steps S1001 to S1003.
S1001, the computer equipment determines the total received power corresponding to each sampling time according to the received power of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time.
In the present embodiment, the total power P is received R The expression of (2) is:
the variation of the total received power with time is shown in fig. 11.
S1002, the computer equipment determines the time delay of each multipath component corresponding to each sampling moment according to the path of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment.
S1003, the computer equipment determines root mean square delay spread corresponding to each sampling time according to the total received power corresponding to each sampling time and the delay and the received power of each multipath component corresponding to each sampling time.
Specifically, the root mean square delay spread σ is calculated as follows:
wherein p is i And t i Representing the received power and time delay of the ith multipath component,representing the average delay of all multipath components.
The embodiment of the invention also provides a complicated mobile time-varying wireless channel simulation device based on ray tracing, fig. 12 is a schematic diagram of a possible complicated mobile time-varying wireless channel simulation device based on ray tracing, and as shown in fig. 12, the device may include an acquisition module 121, a first determination module 122, a second determination module 123, and a third determination module 124.
The obtaining module 121 is configured to obtain target data of a virtual internet of things environment in response to an operation of configuring the virtual internet of things environment by a user, where the virtual internet of things environment is an environment composed of a plurality of objects, the plurality of objects include a transmitting end and a receiving end in communication with the transmitting end, the plurality of objects include dynamic objects, the target data includes environment parameters, antenna parameters, reflection orders, diffraction orders and sampling moments of the virtual internet of things environment, and the antenna parameters include the number of antennas of the transmitting end and the number of antennas of the receiving end;
The first determining module 122 is configured to determine, according to the environmental parameter, the antenna parameter, the reflection order and the diffraction order, a radiation parameter of each first reflected radiation, a radiation parameter of each first diffracted radiation and a radiation parameter of each mixed radiation, which correspond to each pair of transceiver antennas at each sampling moment, by using a preset radiation tracking algorithm; the pair of receiving and transmitting end antennas comprise an antenna of a transmitting end and an antenna of a receiving end, the mixed rays are mixed rays of second reflected rays and second diffracted rays, and the ray parameters comprise paths and receiving power;
a second determining module 123, configured to determine, according to the environmental parameter and the radiation parameter of each first reflected radiation corresponding to each pair of transceiver antennas at each sampling time, the radiation parameter of each scattered radiation corresponding to each pair of transceiver antennas at each sampling time by using a preset scattering model;
the third determining module 124 is configured to determine a channel simulation result according to a radiation parameter of each multipath component corresponding to each pair of transceiver antennas at each sampling time, where the multipath components include a direct radiation, a first reflected radiation, a first diffracted radiation, a mixed radiation, and a scattered radiation, and the channel simulation result includes a doppler shift, where the doppler shift is used to represent an effect of a motion of an object on a transmission characteristic of a channel, and the channel is a signal transmission channel between a transmitting end and a receiving end.
Optionally, the channel simulation result further includes a root mean square delay spread corresponding to each sampling time, where the root mean square delay spread is used to represent an influence of the multipath effect on the transmission characteristic of the channel.
Optionally, the antenna parameters further include a gain of the transmitting end and a transmitting power of the transmitting end, the scattering model is a single scattering lobe model, and the second determining module 123 is specifically configured to: determining each scattering surface corresponding to each pair of receiving-transmitting end antennas at each sampling time according to the path of each reflected ray corresponding to each pair of receiving-transmitting end antennas at each sampling time; determining a path of each corresponding scattered ray in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment, wherein the first included angle is a maximum included angle between each corresponding scattered ray and a reflected ray in each scattering surface; and determining the receiving power of each corresponding scattered ray in each scattering surface at each sampling moment according to the path of each scattered ray included in each scattering surface at each sampling moment, the gain of the transmitting end, the transmitting power of the transmitting end and the scattering coefficient corresponding to each scattering surface.
Optionally, the second determining module 123 is specifically configured to: determining each scattering point included in each scattering surface at each sampling moment according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface at each sampling moment; and determining the path of each scattered ray corresponding to each scattering surface at each sampling moment according to each scattering point included in each scattering surface at each sampling moment, and the transmitting end and the receiving end at the sampling moment.
Optionally, the second determining module 123 is specifically configured to: determining a second included angle corresponding to each scattered ray included in each scattering surface according to a preset scattering valve width factor, the number of scattered rays and a first included angle of each scattering surface, wherein the second included angle is an included angle between the scattered ray and the reflected ray, and the second included angle is smaller than the first included angle; and determining each scattering point included in each scattering surface according to the second included angle corresponding to each scattering ray included in each scattering surface.
Optionally, at least one of the transmitting end and the receiving end is a dynamic object, the antenna parameter further includes a frequency of a simulated radio wave emitted by the transmitting end, the environmental parameter includes a motion path and a motion speed of the transmitting end or a motion path and a motion speed of the receiving end, and the third determining module 124 is specifically configured to: determining the relative speed of the transmitting end and the receiving end at each sampling moment according to the motion path and the motion speed of the transmitting end or the motion path and the motion speed of the receiving end; determining normalized channel capacity according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining Doppler frequency shift corresponding to each sampling time according to the path, the frequency and the relative speed of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling time.
Optionally, the third determining module 124 is further configured to: determining a channel transmission matrix according to the ray parameters of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining the normalized channel capacity according to the channel transmission matrix and a preset signal-to-noise ratio.
Optionally, the third determining module 124 is specifically configured to: determining the total received power corresponding to each sampling moment according to the received power of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; determining the time delay of each multipath component corresponding to each sampling moment according to the path of each multipath component corresponding to each pair of receiving and transmitting end antennas at each sampling moment; and determining the root mean square delay spread corresponding to each sampling time according to the total received power corresponding to each sampling time and the delay and the received power of each multipath component corresponding to each sampling time.
Optionally, the obtaining module 121 is specifically configured to: determining environmental parameters of an environment in which the channel is located in response to an environmental layout operation of a user, the environmental parameters including a shape, a position, a material type of each static object included in the environment, and a shape, an initial position, a material type, a motion path, and a motion speed of each dynamic object; determining antenna parameters in response to antenna configuration operation of a user, wherein the antenna parameters comprise frequency of simulated radio waves sent by a transmitting end, gain of the transmitting end, transmitting power of the transmitting end, moving path and moving speed of a receiving end, antenna types and the number of antennas of the transmitting end, and antenna types and the number of antennas of the receiving end; determining a reflection order and a diffraction order in response to a ray tracing configuration operation of a user; each sampling instant is determined in response to a sampling configuration operation by a user.
Of course, the complex mobile time-varying wireless channel simulation device based on ray tracing provided by the embodiment of the invention comprises but is not limited to the above modules.
The device for simulating the complicated mobile time-varying wireless channel based on the ray tracing provided by the embodiment of the invention is used for executing the complicated mobile time-varying wireless channel simulation method based on the ray tracing in the embodiment, so that the same technical effect as that of the complicated mobile time-varying wireless channel simulation method based on the ray tracing in the embodiment can be achieved.
Another embodiment of the present invention also provides a computer apparatus including: a processor and a memory; the memory is for storing computer program code, the computer program code comprising computer instructions; when the processor executes the computer instructions, the computer device performs the complex mobile time-varying wireless channel simulation method based on ray tracing as shown in the above method embodiments.
Another embodiment of the present invention further provides a computer readable storage medium having stored thereon computer instructions that, when executed on a computer device, cause the computer device to perform the complex mobile time-varying wireless channel simulation method based on ray tracing as shown in the above method embodiment.
Another embodiment of the present invention also provides a computer program product comprising computer instructions which, when run on a computer device, cause the computer device to perform the complex mobile time-varying radio channel simulation method based on ray tracing as shown in the above method embodiment.
The foregoing is merely illustrative of specific embodiments of the present invention, and the scope of the present invention is not limited thereto, but any changes or substitutions within the technical scope of the present invention should be covered by the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.