CN113162712B - Multi-beam channel modeling method based on propagation diagram theory - Google Patents
Multi-beam channel modeling method based on propagation diagram theory Download PDFInfo
- Publication number
- CN113162712B CN113162712B CN202110103746.5A CN202110103746A CN113162712B CN 113162712 B CN113162712 B CN 113162712B CN 202110103746 A CN202110103746 A CN 202110103746A CN 113162712 B CN113162712 B CN 113162712B
- Authority
- CN
- China
- Prior art keywords
- scattering
- transmitting
- receiving
- point
- antenna array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
- Radio Transmission System (AREA)
Abstract
The invention provides a multi-beam channel modeling method based on a propagation diagram theory, which corrects the propagation path gain related to a transmitting and receiving end by calculating the gain of an antenna array at the transmitting and receiving end at each angle, simultaneously considers the mobility of an obstacle, divides scatterers into a static discrete scattering point set and a dynamic discrete scattering point set, and considers different scattering times for two different scattering point sets to construct a multi-beam channel model based on the propagation diagram theory. The method makes up the defect that beam forming is not considered in the existing propagation diagram channel modeling theory, and improves the accuracy of channel modeling.
Description
Technical Field
The invention relates to the technical field of wireless communication, in particular to a multi-beam channel modeling method based on a propagation diagram theory.
Background
In 2019, 5G commercial license plates will be released to enter the original year of commercial formally in China. As one of the key technologies in 5G, the Multiple Input Multiple Output (MIMO) technology can greatly increase the data rate and greatly increase the bandwidth utilization by using space resources under the condition of limited bandwidth, and thus attracts extensive attention in the academic and industrial fields. The capacity improvement effect of the MIMO technology is mainly related to the channel propagation characteristics, so that an MIMO channel simulation model reflecting the real channel characteristics is accurately and effectively established, and the method is particularly important for the optimization of an MIMO system model, the subsequent signal processing algorithm, the evaluation of system performance and the like. One type of beamforming is a signal preprocessing technique based on an MIMO antenna array, which can generate a highly directional radiation pattern by adjusting parameters of array elements of the antenna array according to a certain rule based on the interference principle of electromagnetic waves. An extremely strong or weak beam pointing in a specific direction can be formed, so that the purpose of enhancing the useful signal and simultaneously suppressing the interference is achieved.
Today, wireless channel models can be divided into two broad categories: empirical models and deterministic models. The empirical model is mainly used for counting the overall variation trend of parameters through a large amount of measured data, such as a double-slope path loss model for describing large-scale fading and a random model based on geometry for describing small-scale fading. As an empirical model for a large number of statistics, the empirical model has a wide range of application scenarios, but it is difficult to accurately characterize the radio propagation characteristics in a specific scenario. The deterministic channel modeling model can just make up for the defects of the empirical model.
Typical deterministic channel modeling theories are, for example, ray tracing-based time domain channel modeling theories and propagation map-based frequency domain channel modeling theories. The propagation diagram channel modeling theory models physical phenomena such as reflectors, scatterers and propagation encountered by radio waves in a propagation process into a relation between a vertex and an edge in computer science, and combines the cascading characteristic of a frequency domain transfer function to obtain the frequency domain channel transfer function of the whole process. The propagation diagram theory has a good effect in representing the scattering phenomenon of radio waves. Compared with other deterministic channel modeling methods (ray tracing-based method, FDTD) and the like, the method has the characteristics of small calculated amount, accurate result and wide applicability.
However, the propagation map channel modeling theory is based on the ray theory, and the rays are independent from each other, so that no interference phenomenon exists, and beam forming cannot be realized. The propagation signals among the multiple antennas are processed independently, and beam forming caused by correlation among the antennas is not considered.
Disclosure of Invention
The embodiment of the invention provides a multi-beam channel modeling method based on a propagation diagram theory, which is used for solving the problems that the existing propagation diagram channel modeling theory is incompatible with beam forming and cannot accurately model a channel under a multi-antenna scene.
In order to achieve the purpose, the invention adopts the following technical scheme.
A multi-beam channel modeling method based on propagation diagram theory comprises the following steps:
according to the parameters of the receiving and transmitting antenna array, the beam direction of the receiving and transmitting antenna array and the beam gain of the receiving and transmitting antenna array at each angle are obtained;
correcting the propagation gain of a propagation edge according to the beam gain of the beam of the transmitting-receiving antenna array at each angle;
and establishing a multi-beam channel model by using a propagation diagram channel modeling theory according to the size and the position of the obstacle and the position relation between the obstacle and the transmitting and receiving antenna array element to obtain the transmission function of the transmitting and receiving antenna array.
Preferably, the obtaining of the beam pointing direction of the transceiving antenna array and the beam gain of the beam of the transceiving antenna array at each angle according to the parameter of the transceiving antenna array comprises:
each provided with a transceiver antenna array
And
a track beam of
And
the direction of the track beam is respectively
And
obtaining the angle omega of each beam of the transmitting end antenna array Tx Has a gain of
Each beam of the receiving end antenna array is at an angle omega Rx A gain of
In the formula, superscripts Tx and Rx respectively represent a transmitting end and a receiving end, superscripts dir represents beam direction, t and r respectively represent the serial numbers of the beams at the transmitting and receiving ends, and theta and phi represent the azimuth angle and the horizontal angle of the beams relative to the antenna array;
setting a certain beam generated by a two-dimensional uniform plane antenna array, and converting the gain of the certain beam at a certain angle into gain
G(θ,φ)＝G 0 (θ,φ)*F (1)；
In formula (1), G 0 (θ, φ) represents the single antenna gain at that angle, and F represents the antenna factor, which can be expressed in more detail as:
in the formula (2), j is an imaginary number, theta 0 And phi 0 Indicating the azimuth angle and the horizontal angle corresponding to the beam direction, M, N respectively indicating the number of array elements along the transverse longitudinal axis of the antenna array, k indicating the beam, d x ,d y Respectively represents the distance between two adjacent array elements in the direction of the horizontal and vertical axes, I m,n The weight coefficients on the mth antenna on the horizontal axis and the nth antenna on the vertical axis are shown.
Preferably, the angle Ω of each beam of the transmitting end antenna array is obtained Tx Has a gain of
The method comprises the following steps:
based on the linear additivity of the influence of the multi-channel wave beams in the angle domain, combining the formulas (1) and (2), obtaining the angle omega of the transmitting-end antenna array Tx ＝(θ Tx ,φ Tx ) Multiple beam gain on
And receiving end antenna array is in certain angle omega Rx Multiple beam gain on
Preferably, the modifying the propagation gain of the propagation edge according to the beam gain of the beam of the transceiving antenna array at each angle comprises:
constructing a point set V and an edge set E according to the propagation diagram theory; in which the set of points is divided into transmissionsSet of points V T Set of scattering points V S And a set of receiving points V R . The edge set E is divided into an edge set E from a transmitting point to a receiving point D Set of edges from emission point to scattering point E T Set of edges from scattering point to scattering point E S And set of scattering point to receiver point edges E R (ii) a Let an arbitrary edge e, which corresponds to a frequency domain transfer function of
In the formula, τ e Which is indicative of the propagation delay time,
is a phase random variable uniformly distributed in [0,2 π), g e Represents the path propagation gain;
for a set of scattering points V S Representing a scattering object by grid-shaped scattering points, and setting the distance between two adjacent scattering points as c/B, wherein c is the light speed, and B is the system bandwidth;
setting the Euclidean distance between two scattering points to satisfy
d th,min ≤||r p -r q ||≤d th,max (6)；
Wherein d is th,min And d th,max Representing two distance thresholds, r p And r q Representing vectors corresponding from the origin of the coordinate system to the scattering points p and q;
dividing a scattering point set into a static discrete point set according to the mobility of the obstacle
And a plurality of dynamic discrete point sets
For the dynamic scattering point, the generation of electric wave in the interior of the scattering point is determined according to the concave-convex property of the obstacle
Limited order scattering;
according to the propagation diagram channel modeling theory, the transmission function of the edge is shown as the formula (5), and the transmission gain g of each type of edge e Is composed of
Where f denotes the frequency, for E ∈ E T ,E R ，
Represents the average propagation delay from the emission point to the set of scattering points, where | represents the potential of the set,
representing an adjustment factor, g represents a scattering point gain factor, and odi (e) represents the number of edges from one scattering point to other scattering points;
according to equation (7), the propagation edges associated with the set of transmission points and the set of reception points are modified to
Wherein omega Tx And Ω Rx Respectively, the angle that the propagation edge makes with the transmit/receive antenna array.
Preferably, the establishing a multi-beam channel model by using a propagation diagram channel modeling theory according to the size and the position of the obstacle and the position relationship between the obstacle and the transceiver antenna array element, and the obtaining of the transmission function of the transceiver antenna array comprises:
the transmission function of the receiving and transmitting antenna array comprises a transmission function matrix D (t, f) from a transmitting end to a receiving end, a transmission function T (t, f) from the transmitting end to a scattering point, a transmission function B (t, f) between the scattering points and a transmission function R (t, f) from the scattering points to the receiving end; a transmission end to receiving end transmission function matrix D (t, f), a transmission end to scattering point transmission function T (t, f) and a scattering point to receiving end transmission function R (t, f) are subjected toG Tx (Ω Tx ) And G Rx (Ω Rx ) The influence of (a);
when the beam generates infinite scattering in the static scattering point set, setting the channel transmission function from the transmitting end to the receiving end through the scattering point
The calculation formula is
Wherein the matrix subscript represents the dimension of the matrix;
when the beams respectively occur in the M dynamic scattering point sets
When the scattering is carried out for a limited time, the matrix of the channel transfer function from the transmitting end to the receiving end through the scattering point is set as
In the formula
Representing the number of scatter points, N, in the ith set of dynamic discrete scatter points T And N R Respectively representing the number of elements in a transmitting point set and a receiving point set, namely the number of transmitting antennas and receiving antennas;
the total transmission function of the system is from the transmitting end to the receiving end
The transmitting terminal is gathered to the receiving terminal through the static discrete scattering points
And the transmitting terminal is gathered to the receiving terminal through a plurality of dynamic discrete scattering points
The sum of the transfer functions of the three is obtained
Each element of the matrix in equation (11) represents a transfer function of one edge.
It can be seen from the technical solutions provided by the embodiments of the present invention that, the present invention provides a multi-beam channel modeling method based on a propagation diagram theory, including: the transmitting and receiving end is an antenna array, and a plurality of beams which can be generated by controlling the phase or the weight of each antenna unit signal; each generated beam has a certain directional characteristic and a certain width; due to the influence of multi-beam superposition, the gain of the antenna array has selectivity on all angles; when a propagation diagram channel modeling theory is used for constructing a scattering point set, the minimum distance between scattering points is required to be larger than the multipath resolution; dividing a scattering point set into a static discrete scattering point set and a plurality of dynamic discrete scattering point sets according to the mobility of the obstacle; radio waves are scattered infinitely times in the static discrete scattering point set, and scattered limitedly times in the dynamic discrete scattering point set; the propagation edge associated with the transmitting end is affected by the originating beam; the propagation edge associated with the receiving end is affected by the beam at the receiving end; giving a propagation function of each type of edge; the final system transfer function is the sum of three parts: namely, a transmission end → a visual distance transmission function of a receiving end, a transmission end → a static scattering point (infinite order reflection) → a transmission function of a receiving point, a transmission point → a dynamic scattering point (finite order reflection) → a transmission function of the receiving point. According to the method, on one hand, the influence of the beam generated by the antenna array is considered, on the other hand, different influences of the dynamic/static discrete scattering point set on the transmission function are considered, the defect that beam forming is not considered in the existing propagation diagram channel modeling theory is overcome, and the accuracy of channel modeling is improved.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a processing flow chart of a multi-beam channel modeling method based on propagation diagram theory according to the present invention;
fig. 2 is a schematic diagram of multi-beam channel modeling based on a propagation diagram theory in the multi-beam channel modeling method based on the propagation diagram theory provided by the present invention;
fig. 3 is a diagram illustrating different types of transfer functions.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
Referring to fig. 1, the present invention provides a multi-beam channel modeling method based on propagation diagram theory, mainly studying a propagation diagram channel modeling method with beam forming capability, including the following steps:
according to the parameters of the receiving and transmitting antenna array, the beam direction of the receiving and transmitting antenna array and the beam gain of the beam of the receiving and transmitting antenna array at each angle are obtained;
correcting the propagation gain of a propagation edge according to the beam gain of the beam of the transmitting-receiving antenna array at each angle;
and establishing a multi-beam channel model by using a propagation diagram channel modeling theory according to the size and the position of the obstacle and the position relation between the obstacle and the transmitting and receiving antenna array element to obtain the transmission function of the transmitting and receiving antenna array.
In the embodiment provided by the invention, the transceiving end adopts a multi-antenna array, and generates a plurality of directional beams by controlling the phase or weight of each array element signal.
Further, each beam has a certain directivity and width, so that the antenna array gain is angularly selective, the gain is larger at an angle close to the center of the beam, and the gain is smaller at an angle far from the center of the beam. Therefore, in some preferred embodiments, the second step specifically includes:
as shown in FIG. 2, the transmit-receive antenna arrays are provided with
And
a track beam of
And
the direction of the track beam is respectively
And
obtaining the angle omega of each beam of the transmitting end antenna array Tx Has a gain of
Each beam of the receiving end antenna array is at an angle omega Rx A gain of
In the formula, superscripts Tx and Rx respectively represent a transmitting end and a receiving end, superscripts dir represents beam direction, t and r respectively represent the serial numbers of the beams at the transmitting and receiving ends, theta and phi represent the azimuth angle and the horizontal angle of the beams relative to the antenna array, and theta and phi represent the azimuth angle and the horizontal angle of the beams relative to the antenna array;
the above
And
the antenna array can be determined by the number of antennas, the distance between array elements and the direction of a beam finger;
setting a certain beam to be generated by a two-dimensional plane antenna array, and converting the gain of the certain beam at a certain angle into gain
G(θ,φ)＝G 0 (θ,φ)*AF (1)；
In formula (1), G 0 (θ, φ) represents the single antenna gain at that angle, and F represents the antenna factor, which can be expressed in more detail as:
in the formula (2), j is an imaginary number, theta 0 And phi 0 Indicating the azimuth angle and the horizontal angle corresponding to the beam direction, M, N respectively indicating the number of array elements along the transverse longitudinal axis of the antenna array, k indicating the beam, d x ,d y Respectively represents the distance between two adjacent array elements in the direction of the horizontal and vertical axes, I m,n The weight coefficients on the mth antenna on the horizontal axis and the nth antenna on the vertical axis are shown.
Further, since the influence of the multi-channel beam is linearly additive in the angular domain, equations (1) to (5) can be combined to obtain the antenna array at the transmitting end at a certain angle Ω Tx ＝(θ Tx ,φ Tx ) Multiple beam gain on
And receiving end antenna array is in certain angle omega Rx Multiple beam gain on
In the preferred embodiment provided by the invention, when the frequency domain channel modeling is performed by adopting the propagation map channel modeling theory, the transmission end, the receiving end and the scattering object are abstracted into points in the graph theory by the theory, edges between the points are regarded as electric wave propagation paths, and each edge corresponds to a corresponding frequency domain transmission function. Specifically, the sub-step of correcting the propagation gain of the propagation edge according to the beam gain of the beam of the transmitting-receiving antenna array at each angle specifically includes the following steps:
constructing a point set V and an edge set epsilon according to a propagation diagram theory; wherein the point set is divided into a set of transmission points V T Set of scattering points V S And a set of receiving points V R . The edge set E is divided into an edge set E from a transmitting point to a receiving point D Set of edges from emission point to scattering point E T Set of edges from scattering point to scattering point E S And set of scattering point to receiver point edges E R (ii) a (ii) a For any edge e, its corresponding frequency domain transfer function can be expressed as:
wherein, tau e Which is indicative of the propagation delay time,
is a phase random variable which is uniformly distributed in [0,2 pi), g e Represents the path propagation gain;
when the scattering points are used for representing scattering objects, grid-shaped scattering points are adopted, one scattering point represents all the scattering points in a small part of area on a scattering surface, and the distance between two adjacent scattering points is c/B, wherein c is the light speed, and B is the system bandwidth;
when the beam is scattered inside the scattering point, only part of the links can occur; if and only if the visual distance between two scattering points is reachable, and the Euclidean distance satisfies:
d th,min ≤||r p -r q ||≤d th,max (6)，
wherein d is th,min And d th,max Is two thresholds, r p And r q Representing vectors corresponding from the origin of the coordinate system to the scattering points p and q;
as shown in fig. 3, the set of scattering points is divided into two categories: a static set of discrete points
And a plurality of dynamic discrete point sets
Wherein radio waves are scattered infinitely times in the static discrete scattering point set, and scattered limitedly times in the dynamic discrete scattering point set;
according to the propagation diagram channel modeling theory, the transmission function of the edge is shown as the formula (5), and the transmission gain g of each type of edge e Is composed of
Where f denotes the frequency, for e ∈ ε t ,ε r ，
Represents the average propagation delay from the emission point to the set of scattering points, where | represents the potential of the set,
the adjustment factor is expressed as a whole, g represents the scattering point gain factor, and odi (e) represents the number of edges from one scattering point to the other.
Further, according to equation (7), propagation edges associated with the set of transmission points and the set of reception points are modified to
Wherein omega Tx And Ω Rx Respectively, the angle that the propagation edge makes with the transmit/receive antenna array.
In the preferred embodiment provided by the invention, the transmission function for obtaining the propagation diagram channel modeling theory according to the size and the position of the obstacle and the position relation between the obstacle and the transmitting and receiving antenna array element comprises four parts: transmitting end → receiving end transfer function matrix D (t, f), transmitting end → scattering point transfer function T (t, f), scattering point → scattering point transferFunction B (t, f), scattering point → receiver transfer function R (t, f); wherein, the transmission terminal → receiving terminal transfer function matrix D (t, f), the transmission terminal → scattering point transfer function T (t, f), the scattering point → receiving terminal transfer function R (t, f) are all received by the wave beam G generated by the antenna array Tx (Ω Tx ) And G Rx (Ω Rx ) The influence of (c).
When the electric wave generates infinite scattering in the static scattering point set, the channel transfer function of the transmitting end → the scattering point → the receiving end
The calculation formula is
When the beams respectively occur in the M dynamic scattering point sets
For the finite scattering, the channel transfer function matrix of the transmitting terminal → the scattering point → the receiving terminal is the superposition of M transfer functions, which can be expressed as:
wherein
Representing the number of scatter points, N, in the ith set of dynamic discrete scatter points T And N R Respectively representing the number of elements in a transmitting point set and a receiving point set, namely the number of transmitting antennas and receiving antennas;
the total transfer function of the system is transmitting end → receiving end
Transmitting end → static discrete scattering point set → receiving end
Transmitting terminal → several dynamic state discrete scattering point sets → receiving terminal
The sum of the three partial transfer functions can be expressed as:
each element of the matrix in equation (11) represents a transfer function of one edge.
In summary, the present invention provides a multi-beam channel modeling method based on propagation diagram theory, including: the receiving and transmitting end is an antenna array, and a plurality of beams can be generated by controlling the phase or weight of each antenna unit signal; each generated beam has a certain directional characteristic and a certain width; due to the influence of multi-beam superposition, the gain of the antenna array is selective at various angles; when a propagation diagram channel modeling theory is used for constructing a scattering point set, the minimum distance between scattering points is required to be larger than the multipath resolution; dividing a scattering point set into a static discrete scattering point set and a plurality of dynamic discrete scattering point sets according to the mobility of the obstacle; radio waves are scattered infinitely times in the static discrete scattering point set, and scattered limitedly times in the dynamic discrete scattering point set; the propagation edge associated with the transmitting end is affected by the originating beam; the propagation edge associated with the receiving end is affected by the beam at the receiving end; giving a propagation function of each type of edge; the final system transfer function is the sum of three parts: namely, a transmission end → a visual distance transmission function of a receiving end, a transmission end → a static scattering point (infinite order reflection) → a transmission function of a receiving point, a transmission point → a dynamic scattering point (finite order reflection) → a transmission function of the receiving point. The method corrects the propagation path gain related to the transmitting and receiving ends by calculating the gain of the antenna array of the transmitting and receiving ends at each angle, simultaneously considers the mobility of the barrier, divides the scattering objects into a static discrete scattering point set and a dynamic discrete scattering point set, and considers different scattering times for the two different scattering point sets to construct a multi-beam channel model based on the propagation diagram theory; the defect that beam forming is not considered in the existing propagation diagram channel modeling theory is overcome, and the accuracy of channel modeling is improved.
Those of ordinary skill in the art will understand that: the figures are merely schematic representations of one embodiment, and the blocks or flow diagrams in the figures are not necessarily required to practice the present invention.
From the above description of the embodiments, it is clear to those skilled in the art that the present invention can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which may be stored in a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the embodiments or some parts of the embodiments.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for apparatus or system embodiments, since they are substantially similar to method embodiments, they are described in relative terms, as long as they are described in partial descriptions of method embodiments. The above-described embodiments of the apparatus and system are merely illustrative, and the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (1)
1. A multi-beam channel modeling method based on propagation diagram theory is characterized by comprising the following steps:
according to the parameters of the receiving and transmitting antenna array, the beam direction of the receiving and transmitting antenna array and the beam gain of the beam of the receiving and transmitting antenna array at each angle are obtained; the method specifically comprises the following steps:
provided with a transmitting-receiving antenna array respectively having
And
a track beam of
And
the direction of the track beam is respectively
And
obtaining the angle omega of each beam of the transmitting end antenna array Tx Has a gain of
Each beam of the receiving end antenna array is at an angle omega Rx A gain of
In the formula, superscripts Tx and Rx respectively represent a transmitting end and a receiving end, superscripts dir represents beam direction, t and r respectively represent the serial numbers of the beams at the transmitting and receiving ends, and theta and phi represent the azimuth angle and the horizontal angle of the beams relative to the antenna array;
setting a certain beam generated by a two-dimensional uniform plane antenna array, and converting the gain of the certain beam at a certain angle into gain
G(θ,φ)＝G 0 (θ,φ)*F (1)；
In formula (1), G 0 (θ, φ) represents the single antenna gain at that angle, and F represents the antenna factor, which can be expressed in more detail as:
in the formula (2), j is an imaginary number, theta 0 And phi 0 Indicating the azimuth angle and the horizontal angle corresponding to the beam direction, M, N respectively indicating the number of array elements along the transverse longitudinal axis of the antenna array, k indicating the beam, d x ,d y Respectively represents the distance between two adjacent array elements in the direction of the horizontal and vertical axes, I m,n Representing the weight coefficient on the mth antenna on the horizontal axis and the nth antenna on the vertical axis;
the obtaining of the angle omega of each wave beam of the transmitting terminal antenna array Tx Has a gain of
The method comprises the following steps:
based on the linear additivity of the influence of the multi-channel wave beams in the angle domain, combining the formulas (1) and (2), obtaining the angle omega of the transmitting-end antenna array Tx ＝(θ Tx ,φ Tx ) Multiple beam gain on
And receiving end antenna array is in certain angle omega Rx Multiple beam gain on
Correcting the propagation gain of a propagation edge according to the beam gain of the beam of the transmitting-receiving antenna array at each angle; the method specifically comprises the following steps:
constructing a point set V and an edge set E according to the propagation diagram theory; wherein the point set is divided into a set of transmission points V T Set of scattering points V S And a set of receiving points V R (ii) a The edge set E is divided into an edge set E from a transmitting point to a receiving point D Set of edges from emission point to scattering point E T Set of edges from scattering point to scattering point E S And set of scattering point to receiver point edges E R (ii) a Let an arbitrary edge e, which corresponds to a frequency domain transfer function of
In the formula, τ e Which is indicative of the propagation delay time,
is a phase random variable which is uniformly distributed in [0,2 pi), g e Represents the path propagation gain;
for a set of scattering points V S Representing a scattering object by grid-shaped scattering points, and setting the distance between two adjacent scattering points as c/B, wherein c is the light speed, and B is the system bandwidth;
setting the Euclidean distance between two scattering points to satisfy
d th,min ≤||r p -r q ||≤d th,max (6)；
Wherein d is th,min And d th,max Representing two distance thresholds, r p And r q Representing vectors corresponding from the origin of the coordinate system to the scattering points p and q;
dividing the set of scattering points into one according to the mobility of the obstacleA static set of discrete points
And a plurality of dynamic discrete point sets
For the dynamic scattering point, the generation of electric wave in the interior of the scattering point is determined according to the concave-convex property of the obstacle
Limited order scattering;
according to the propagation diagram channel modeling theory, the transmission function of the edge is shown as the formula (5), and the transmission gain g of each type of edge e Is composed of
Where f denotes the frequency, for E ∈ E T ,E R ，
Represents the average propagation delay from the emission point to the set of scattering points, where | represents the potential of the set,
representing an adjustment factor, g represents a scattering point gain factor, and odi (e) represents the number of edges from one scattering point to other scattering points;
according to equation (7), the propagation edges associated with the set of transmission points and the set of reception points are modified to
Wherein omega Tx And Ω Rx Respectively representing the angles of the propagation edge and the transmitting/receiving antenna array;
establishing a multi-beam channel model by using a propagation diagram channel modeling theory according to the size and the position of the obstacle and the position relation between the obstacle and the array element of the transmitting and receiving antenna to obtain a transmission function of the transmitting and receiving antenna array; the method specifically comprises the following steps:
the transmission function of the receiving and transmitting antenna array comprises a transmission function matrix D (t, f) from a transmitting end to a receiving end, a transmission function T (t, f) from the transmitting end to a scattering point, a transmission function B (t, f) between scattering points and a transmission function R (t, f) from the scattering points to the receiving end; the matrix D (t, f) of transmission function from the transmitting end to the receiving end, the T (t, f) of transmission function from the transmitting end to the scattering point and the R (t, f) of transmission function from the scattering point to the receiving end are subjected to G Tx (Ω Tx ) And G Rx (Ω Rx ) The influence of (c);
when the beam generates infinite scattering in the static scattering point set, setting the channel transmission function from the transmitting end to the receiving end through the scattering point
The calculation formula is
Wherein the matrix subscript represents the dimension of the matrix;
when the beams respectively occur in the M dynamic scattering point sets
When the scattering is carried out for a limited time, the matrix of the channel transfer function from the transmitting end to the receiving end through the scattering point is set as
In the formula
Representing the ith set of dynamic discrete scattering pointsNumber of inner scattering points, N T And N R Respectively representing the number of elements in a transmitting point set and a receiving point set, namely the number of transmitting antennas and receiving antennas;
the total transmission function of the system is from the transmitting end to the receiving end
The transmitting terminal is gathered to the receiving terminal through the static discrete scattering points
And the transmitting terminal is gathered to the receiving terminal through a plurality of dynamic discrete scattering points
The sum of the transfer functions of the three is obtained
Each element of the matrix in equation (11) represents a transfer function of one edge.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110103746.5A CN113162712B (en) | 2021-01-26 | 2021-01-26 | Multi-beam channel modeling method based on propagation diagram theory |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110103746.5A CN113162712B (en) | 2021-01-26 | 2021-01-26 | Multi-beam channel modeling method based on propagation diagram theory |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113162712A CN113162712A (en) | 2021-07-23 |
| CN113162712B true CN113162712B (en) | 2022-08-16 |
Family
ID=76879143
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110103746.5A Active CN113162712B (en) | 2021-01-26 | 2021-01-26 | Multi-beam channel modeling method based on propagation diagram theory |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113162712B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023133766A1 (en) * | 2022-01-13 | 2023-07-20 | 华为技术有限公司 | Beam control method and apparatus |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111010222A (en) * | 2019-11-13 | 2020-04-14 | 东南大学 | Deep learning-based large-scale MIMO downlink user scheduling method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105721087B (en) * | 2016-03-21 | 2019-01-04 | 北京交通大学 | Theoretical non line of sight class is schemed away from the channel modeling method under tunnel environment based on propagating |
| CN109412723B (en) * | 2017-08-16 | 2020-06-16 | 中兴通讯股份有限公司 | MIMO channel analysis model, modeling method and computer readable storage medium |
| CN109818658A (en) * | 2017-11-21 | 2019-05-28 | 索尼公司 | For the electronic equipment of wireless communication system, method and storage medium |
| CN108833041B (en) * | 2018-04-27 | 2021-05-04 | 东南大学 | Multi-beam low-orbit satellite channel simulation method based on elliptical orbit |
-
2021
- 2021-01-26 CN CN202110103746.5A patent/CN113162712B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111010222A (en) * | 2019-11-13 | 2020-04-14 | 东南大学 | Deep learning-based large-scale MIMO downlink user scheduling method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113162712A (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2347531B1 (en) | Over-the-air test | |
| CN109639378B (en) | Rectangular tunnel wireless propagation channel modeling method | |
| Jiang et al. | Three-dimensional geometry-based stochastic channel modeling for intelligent reflecting surface-assisted UAV MIMO communications | |
| Jiang et al. | Physics-based 3D end-to-end modeling for double-RIS assisted non-stationary UAV-to-ground communication channels | |
| Priebe et al. | AoA, AoD and ToA characteristics of scattered multipath clusters for THz indoor channel modeling | |
| CN106972882B (en) | Self-adaptive beam forming method of co-prime array based on virtual domain space power spectrum estimation | |
| CN109714120B (en) | Method for simulating propagation characteristics of coupling multi-antenna indoor space fading channel | |
| CN108718224B (en) | Wireless communication channel simulation calculation method | |
| CN107104720B (en) | Mutual-prime array self-adaptive beam forming method based on covariance matrix virtual domain discretization reconstruction | |
| Xie et al. | A 3D geometry-based stochastic model for 5G massive MIMO channels | |
| CN110212955B (en) | Ray-based 3D MIMO channel modeling method | |
| CN110958069A (en) | Three-dimensional non-stationary broadband double-cluster channel modeling method based on tunnel environment | |
| Mehrnia et al. | Novel maritime channel models for millimeter radiowaves | |
| CN110620627B (en) | Non-stationary channel modeling method and device for vehicle-to-vehicle multi-antenna system | |
| CN107171703B (en) | Method capable of simulating propagation characteristics of fading signals in indoor multi-antenna communication system | |
| Xiong et al. | A statistical MIMO channel model for reconfigurable intelligent surface assisted wireless communications | |
| CN109951213B (en) | High-altitude platform MIMO three-dimensional geometric random model establishing method and communication method | |
| CN103338094A (en) | Modeling method for multi-input multi-output system channel | |
| CN113162712B (en) | Multi-beam channel modeling method based on propagation diagram theory | |
| Neekzad et al. | Comparison of ray tracing simulations and millimeter wave channel sounding measurements | |
| Jiang et al. | Map-based UAV mmWave channel model and characteristics analysis | |
| Bai et al. | A mixed-bouncing based non-stationary model for 6G massive MIMO mmWave UAV channels | |
| CN113949474A (en) | Unmanned aerial vehicle geometric model establishing method based on assistance of intelligent reflecting surface | |
| Cao et al. | Research on sea surface Ka-band stochastic multipath channel modeling | |
| CN104808178B (en) | A kind of airborne radar transmitting pattern method for designing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |