CN110212955B - Ray-based 3D MIMO channel modeling method - Google Patents

Ray-based 3D MIMO channel modeling method Download PDF

Info

Publication number
CN110212955B
CN110212955B CN201910500886.9A CN201910500886A CN110212955B CN 110212955 B CN110212955 B CN 110212955B CN 201910500886 A CN201910500886 A CN 201910500886A CN 110212955 B CN110212955 B CN 110212955B
Authority
CN
China
Prior art keywords
angle
antenna
horizontal
path
sub
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
Application number
CN201910500886.9A
Other languages
Chinese (zh)
Other versions
CN110212955A (en
Inventor
武畅
魏兆林
李�杰
游远
阿天仁
杨帆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910500886.9A priority Critical patent/CN110212955B/en
Publication of CN110212955A publication Critical patent/CN110212955A/en
Application granted granted Critical
Publication of CN110212955B publication Critical patent/CN110212955B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/0413—MIMO systems

Abstract

The invention discloses a ray-based 3D MIMO channel modeling method, which realizes a ray-based channel model by combining sub paths and obtains a 3D MIMO channel transmission function Hu，s，n(t), the invention not only considers the factors including channel power characteristic, transmission delay, horizontal angle and pitching angle, but also considers the factors of antenna polarization and correlation, and has the functions of simple algorithm and easy scene switching. The invention further adopts a uniform power sub-path method in the process of solving the horizontal emission angle, the vertical emission angle, the horizontal receiving angle and the vertical receiving angle of the sub-path, thereby greatly reducing the calculation amount and the system complexity.

Description

Ray-based 3D MIMO channel modeling method
Technical Field
The invention relates to the field of multi-antenna wireless communication, in particular to a ray-based 3D MIMO channel modeling method.
Background
In the field of wireless communication, signal transmission is often affected by various environmental noises, human interference, multipath effects and doppler shifts, and all of the above effects have uncertainty, so before designing software and hardware of a communication system, researchers need to apply radio wave theory, information theory, random process theory and statistical signal processing theory to theoretically analyze a wireless communication channel and establish a strict mathematical model to determine the transmission characteristics of signals, and design a communication system or a communication unit and deduce the performance of the communication system or the communication unit on the basis.
There is an important branch in the field of wireless communication, MIMO (Multiple-Input Multiple output) technology. The MIMO technology is to use a plurality of transmitting antennas and receiving antennas at a transmitting end and a receiving end, respectively, so that signals are transmitted and received through the plurality of antennas of the transmitting end and the receiving end, thereby improving communication quality. Early MIMO (multiple input multiple output) channel modeling mainly considers factors such as channel loss, transmission delay, horizontal angle and the like, neglects the pitch angle problem of the channel, and is a typical two-dimensional mathematical model.
With the dramatic increase of mobile data traffic, the research of 5G mobile communication has become popular at present, and 5G communication has higher requirements on information transmission rate, so that a multi-antenna transmission MIMO technology, which is an important link of 5G communication, needs a channel model with higher precision to support the transmission characteristic research and system performance optimization.
At present, on the basis of a traditional MIMO channel two-dimensional mathematical model, the pitching angle factor of a channel is added in the academic world, and a 3D MIMO modeling technology is constructed. However, at present, the academic community mainly uses a spatial correlation matrix to realize the model construction based on the correlation, and on this basis, a specific spectrum function must be used to simulate the time correlation. Once the channel scene of the communication system is changed, the channel model must be reconstructed, so the universality is poor, the computation amount is huge, and the hardware resource of a digital device is consumed greatly.
Disclosure of Invention
Aiming at the defects in the prior art, the method for ray-based 3D MIMO channel modeling provided by the invention solves the problems that the existing 3D MIMO channel modeling technology does not support or is not easy to support communication system channel scene transformation, the universality is poor, the operation amount is huge, and a large amount of digital device hardware resources are consumed, can flexibly realize the channel transformation under various scenes, and supports L OS (line-of-sight transmission, which means that wireless signals are transmitted between a transmitting end and a receiving end without shielding) and N L OS (non-line-of-sight transmission, which means that wireless signals are transmitted in a reflection, scattering and diffraction mode under the condition of obstacles).
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a method of ray-based 3D MIMO channel modeling, comprising the steps of:
s1, setting the total number S of transmitting antennas and the total number U of receiving antennas of the 3D MIMO channel model according to the arrangement of the transmitting antennas and the receiving antennas of the MIMO communication system, and carrying out linear approximate processing on the antenna layout;
s2, determining the type of a communication scene according to the communication environment of the MIMO communication system;
s3, determining the number N of propagation paths between antenna arrays consisting of transmitting antennas and receiving antennas and the number M of sub-paths contained in each path according to the type of a communication scene;
s4, determining the path delay of each path according to the channel environment in the communication scene according to the CASE II model method;
s5, according to the CASE II model method, determining the path power P of each path n according to the path delay of each path and the statistical characteristics of the channel environment path power in the communication scenenN is an integer, N is more than or equal to 1 and less than or equal to N;
s6, according to the uniform power sub-diameter method, calculating the horizontal emission angle α of each sub-diameter M in each diameter n according to the number M of sub-diameters contained in each diametern,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000031
M is an integer, M is more than or equal to 1 and less than or equal to M;
s7, calculating the polarization parameter k of each sub-diameter m in each diameter n according to the statistical characteristics of the polarization rate cross parametersn,m；
S8, obtaining the initial phase of the polarized antenna of each sub-diameter m in each diameter n according to the statistical characteristics of the initial phase of the polarized antenna
Figure BDA0002090183460000032
And
Figure BDA0002090183460000033
is a vertical-vertical cross-polarization initial phase,
Figure BDA0002090183460000034
for the vertical-horizontal cross polarization initial phase,
Figure BDA0002090183460000035
for the horizontal-horizontal cross-polarization initial phase,
Figure BDA0002090183460000036
is a horizontal-vertical cross polarization initial phase;
s9, establishing a 3D MIMO channel time domain expression according to the wireless communication channel principle and the multi-antenna characteristics, and obtaining the path power P of each path nnPolarization parameter k for each sub-diameter m in each diameter nn,mHorizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000037
And its polarization antenna initial phase
Figure BDA0002090183460000038
Determining an nth path communication time domain transmission function H received by an s antenna transmitting an u antenna of a 3D MIMO channelu,s,n() And finishing channel modeling, wherein S and U are positive integers, S is more than or equal to 1 and less than or equal to S, and U is more than or equal to 1 and less than or equal to U.
Further, the method for determining the communication scene type in step S2 is to determine that the communication scene is L OS if the wireless signal is transmitted between the transmitting-end antenna and the receiving-end antenna without being blocked, and otherwise determine that the communication scene is N L OS.
Further, the method of determining the number N of paths of the propagation path between the transmitting antenna and the receiving antenna according to the communication scene type in step S3 is that if the communication scene is N L OS, the number value of the paths of the propagation path is set to L and N is L according to the number of obstacles and the physical electrical characteristics by using the electromagnetic field and electromagnetic wave theory, and if the communication scene is L OS, the number of the paths of the propagation path in L OS scene is L +1 and N is L +1, that is, a direct path is added on the basis of the propagation path in N L OS scene.
Further: the statistical characteristic of the path power in the step S5 is an exponential distribution.
Further: the step S6 specifically includes the following steps:
s61, horizontal emission angle α according to sub-diametern,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000041
The angle power statistics of (a) yields the horizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000042
Respective angular power spectral density PAS function pα(θ)、pβ(θ)、pγ(theta) and
Figure BDA0002090183460000043
wherein the horizontal emission angle αn,mThe angle power statistic characteristic of (1) is Laplace distribution and horizontal receiving angle gamman,mHas an angular power statistic of 35 DEG Laplacian distribution or 360 DEG uniform distribution, vertical emission angle βn,mAnd vertical acceptance angle
Figure BDA0002090183460000044
The statistical characteristics of the angle power of the power sensor are all uniformly distributed at 30 degrees;
s62, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000045
All angles theta that holdiWherein i is 1, 2, … and M, and the angles are recorded as a sub-path horizontal emission angle αn,m；
S63, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000046
All angles theta that holdiWherein i is 1, 2, … and M, and the angles are recorded as a sub-path vertical emission angle βn,m；
S64, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000047
All angles theta that holdiWhere i is 1, 2, …, M, and these angles are taken as the horizontal reception angle γ of the minor diametern,m；
S65, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000048
All angles theta that holdiWherein i is 1, 2, …, M, and the angles are perpendicular receiving angles of the sub-paths
Figure BDA0002090183460000049
Further: the statistical characteristic of the polarization rate cross parameter in step S7 is lognormal distribution, and the expression of the polarization parameter is:
Figure BDA00020901834600000410
where are random variables that follow a gaussian distribution.
Further: the statistical characteristic of the initial phase of the polarized antenna in the step S8 is a uniform distribution from-pi to pi.
Further: the 3D MIMO channel time domain expression established in step S9 according to the wireless communication channel principle and the multi-antenna characteristics is as follows:
Figure BDA0002090183460000051
wherein G istx,s,VFor power gain in the direction of vertical polarization of the s-th transmitting antenna, Gtx,s,HFor power gain in horizontal polarization direction of the s-th transmitting antenna, Grx,u,VFor the power gain in the vertical polarization direction of the u-th receiving antenna, Grx,u,HPower gain, x, for the horizontal polarization direction of the u-th receiving antennas、ysAnd zsIs the spatial coordinate, x, of the s-th transmitting antennau、yuAnd zuFor the spatial coordinates of the u-th receiving antenna,
Figure BDA0002090183460000052
is the inverse of the wavelength of the carrier wave,
Figure BDA0002090183460000053
is the velocity vector of the movable antenna end, t is the current time, [ ·]HIs a matrix conjugate transpose operation.
The invention has the beneficial effects that: the wireless transmission link among the multiple antennas can be regarded as a communication channel consisting of a plurality of ray clusters (multipath for short) with path resolution, each path comprises a plurality of sub-paths which can not be resolved, and the invention realizes a channel model based on rays by combining the sub-paths to obtain a 3D MIMO channel transmission function Hu,s,n(t) of (d). In the parameter establishment process, the number of paths and the number of sub-paths are set by a communication scene, and the path power P of each path nnXPR (polarization parameter) k for each sub-diameter m of each diameter nn,mHorizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000054
And its polarization antenna initial phase
Figure BDA0002090183460000055
Are established by the statistical properties of their device characteristics and the environment in which they are located. The invention not only considers the factors including the channel power characteristic, the transmission delay, the horizontal angle, the pitching angle and the like, but also considers the factors of the antenna polarization and the correlation, and has the functions of simple algorithm and easy scene switching. The invention further adopts a uniform power sub-path method in the process of solving the horizontal emission angle, the vertical emission angle, the horizontal receiving angle and the vertical receiving angle of the sub-path, thereby greatly reducing the calculation amount and the system complexity.
Drawings
FIG. 1 is a schematic diagram of a ray-based 3D MIMO channel modeling process;
FIG. 2 is a schematic diagram of a 3D MIMO antenna and multipath;
FIG. 3 is a schematic diagram of resolvable paths and non-resolvable sub-paths contained therein;
FIG. 4 is a waveform diagram of a 3D MIMO multi-path communication time domain transfer function at a mobile terminal speed of 120 km/h;
FIG. 5 is a waveform diagram of a 3D MIMO multi-path communication time domain transfer function at a mobile terminal speed of 240 km/h;
FIG. 6 is a 3D MIMO channel signal amplitude distribution law;
FIG. 7 is a 3D MIMO channel phase distribution law;
fig. 8 is a comparison graph of spatial correlation of a 3D MIMO channel.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 1, in one embodiment of the present invention, a method for ray-based 3D MIMO channel modeling comprises the steps of:
s1, setting the total number S of transmitting antennas and the total number U of receiving antennas of the 3D MIMO channel model according to the arrangement of the transmitting antennas and the receiving antennas of the MIMO communication system, and carrying out linear approximate processing on the antenna layout;
s2, determining the type of a communication scene according to the communication environment of the MIMO communication system, if a wireless signal is transmitted between a transmitting terminal antenna and a receiving terminal antenna without shielding, judging that the communication scene is L OS, otherwise, judging that the communication scene is N L OS, wherein the communication scene of the embodiment is N L OS;
s3, determining the number N of the propagation paths between the antenna array composed of the transmitting antenna and the receiving antenna and the number M of the sub-paths contained in each path according to the communication scene type, wherein if the communication scene is N L OS, the number value of the propagation paths is set to be L according to the number of barriers and physical and electrical characteristics by electromagnetic field and electromagnetic wave theory, the number value is L, if the communication scene is L OS, the number of the propagation paths in L OS scene is L +1, and N is L +1, namely, a direct path is added on the basis of the propagation paths in N L OS scene, as shown in FIG. 2, in N L OS communication scene, due to the reflection of buildings, a plurality of distinguishable propagation paths exist in a communication link from a specific mobile terminal transmitting antenna to a specific receiving antenna of a base station, as shown in FIG. 3, a plurality of indistinguishable sub-paths are contained in one path;
s4, determining the path delay of each path according to the CASE II model method and the channel environment in the communication scene, wherein the model delay is [ 0310710109017302510 ] ns in the embodiment;
s5, according to CASE II model method, determining path power P of each path n according to path delay of each path and statistical characteristics of channel environment path power in communication scene, namely exponential distributionnN is an integer, N is more than or equal to 1 and less than or equal to N;
s6, according to the uniform power sub-diameter method, calculating the horizontal emission angle α of each sub-diameter M in each diameter n according to the number M of sub-diameters contained in each diametern,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000071
M is an integer, M is more than or equal to 1 and less than or equal to M, and the method comprises the following steps:
s61 horizontal emission angle α from minor diametern,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000072
The angle power statistics of (a) yields the horizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000073
Respective PAS (angular Power spectral Density) function pα(θ)、pβ(θ)、pγ(theta) and
Figure BDA0002090183460000074
wherein the horizontal emission angle αn,mThe angle power statistic characteristic of (1) is Laplace distribution and horizontal receiving angle gamman,mHas an angular power statistic of 35 DEG Laplacian distribution or 360 DEG uniform distribution, vertical emission angle βn,mAnd vertical acceptance angle
Figure BDA0002090183460000075
The statistical characteristics of the angle power of the power sensor are all uniformly distributed at 30 degrees;
s62, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000081
All angles theta that holdiWherein i is 1, 2, … and M, and the angles are recorded as a sub-path horizontal emission angle αn,m；
S63, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000082
All angles theta that holdiWherein i is 1, 2, … and M, and the angles are recorded as a sub-path vertical emission angle βn,m；
S64, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000083
All angles theta that holdiWhere i is 1, 2, …, M, and these angles are taken as the horizontal reception angle γ of the minor diametern,m；
S65, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure BDA0002090183460000084
All angles theta that holdiWherein i is 1, 2, …, M, and the angles are perpendicular receiving angles of the sub-paths
Figure BDA0002090183460000085
S7, calculating XPR (polarization parameter) k of each sub-diameter m in each diameter n according to the statistical characteristics of the polarization rate cross parametersn,mThe statistical characteristic is lognormal distribution, and the expression is as follows:
Figure BDA0002090183460000086
wherein is a random variable subject to a gaussian distribution;
s8, obtaining the initial phase of the polarized antenna of each sub-diameter m in each diameter n according to the statistical characteristic of the initial phase of the polarized antenna, namely the uniform distribution from-pi to pi
Figure BDA0002090183460000087
And
Figure BDA0002090183460000088
is a vertical-vertical cross-polarization initial phase,
Figure BDA0002090183460000089
for the vertical-horizontal cross polarization initial phase,
Figure BDA00020901834600000810
for the horizontal-horizontal cross-polarization initial phase,
Figure BDA00020901834600000811
is a horizontal-vertical cross polarization initial phase;
s9, power P of each path nnXPR (polarization parameter) k for each sub-diameter m of each diameter nn,mHorizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA00020901834600000812
And its polarization antenna initial phaseAccording to the principle of wireless communication channel andestablishing a 3D MIMO channel time domain expression by multi-antenna characteristics, and determining an nth path communication time domain transmission function H received by an s antenna transmitting a u antenna of the 3D MIMO channelu,s,n(t), channel modeling is completed, wherein S and U are positive integers, S is more than or equal to 1 and less than or equal to S, U is more than or equal to 1 and less than or equal to U, and the expression is as follows:
Figure BDA0002090183460000091
wherein G istx,s,VFor power gain in the direction of vertical polarization of the s-th transmitting antenna, Gtx,s,HFor power gain in horizontal polarization direction of the s-th transmitting antenna, Grx,u,VFor the power gain in the vertical polarization direction of the u-th receiving antenna, Grx,u,HPower gain, x, for the horizontal polarization direction of the u-th receiving antennas、ysAnd zsIs the spatial coordinate, x, of the s-th transmitting antennau、yuAnd zuFor the spatial coordinates of the u-th receiving antenna,
Figure BDA0002090183460000092
is the inverse of the wavelength of the carrier wave,
Figure BDA0002090183460000093
is the velocity vector of the movable antenna end, t is the current time, [ ·]HIs a matrix conjugate transpose operation.
The invention has the beneficial effects that: the wireless transmission link among the multiple antennas can be regarded as a communication channel consisting of a plurality of ray clusters (multipath for short) with path resolution, each path comprises a plurality of sub-paths which can not be resolved, and the invention realizes a channel model based on rays by combining the sub-paths to obtain a 3D MIMO channel transmission function Hu,s,n(t) of (d). In the parameter establishment process, the number of paths and the number of sub-paths are set by a communication scene, and the path power P of each path nnXPR (polarization parameter) k for each sub-diameter m of each diameter nn,mHorizontal emission angle αn,mVertical emission angle βn,mHorizontal reception angle gamman,mAnd vertical acceptance angle
Figure BDA0002090183460000094
And its polarization antenna initial phase
Figure BDA0002090183460000095
Are established by the statistical properties of their device characteristics and the environment in which they are located. The invention not only considers the factors including the channel power characteristic, the transmission delay, the horizontal angle, the pitching angle and the like, but also considers the factors of the antenna polarization and the correlation, and has the functions of simple algorithm and easy scene switching. The invention further adopts a uniform power sub-path method in the process of solving the horizontal emission angle, the vertical emission angle, the horizontal receiving angle and the vertical receiving angle of the sub-path, thereby greatly reducing the calculation amount and the system complexity.
The channel transfer function H of each path at a carrier frequency of 6GHz and a movable antenna end speed of 120km/Hu,s,n(t) as shown in FIG. 4; the channel transfer function H of each path when the speed is converted to 240km/Hu,s,n(t) is shown in FIG. 5. It can be seen from the comparison that the larger the moving speed of the mobile terminal is, the larger the doppler frequency shift is, the smaller the channel coherence time is, and the channel coherence time is reflected on the channel, that is, the faster the channel change is, and the modeling effect conforms to the common knowledge of wireless communication.
As known in the prior art, the amplitude of the MIMO channel transfer function in the N L OS communication scenario should be subject to rayleigh distribution, and the phase should be subject to uniform distribution, and the amplitude distribution law and the phase distribution law of the MIMO channel signal obtained in this embodiment are as shown in fig. 6 and fig. 7, and both conform to theoretical expectations.
As can be seen from fig. 8, the 3D MIMO channel model constructed by the present embodiment is closer to the theoretical expectation than the conventional MIMO channel model, and is more consistent with the actual physical channel.

Claims (8)

1. A method for ray-based 3D MIMO channel modeling, comprising the steps of:
s1, setting the total number S of transmitting antennas and the total number U of receiving antennas of the 3D MIMO channel model according to the arrangement of the transmitting antennas and the receiving antennas of the MIMO communication system, and carrying out linear approximate processing on the antenna layout;
s2, determining the type of a communication scene according to the communication environment of the MIMO communication system;
s3, determining the number N of propagation paths between antenna arrays consisting of transmitting antennas and receiving antennas and the number M of sub-paths contained in each path according to the type of a communication scene;
s4, determining the path delay of each path according to the channel environment in the communication scene according to the CASE II model method;
s5, according to the CASE II model method, determining the path power P of each path n according to the path delay of each path and the statistical characteristics of the channel environment path power in the communication scenenN is an integer, N is more than or equal to 1 and less than or equal to N;
s6, according to the uniform power sub-diameter method, calculating the horizontal emission angle α of each sub-diameter M in each diameter n according to the number M of sub-diameters contained in each diametern，mVertical emission angle βn，mHorizontal reception angle gamman，mAnd vertical acceptance angle
Figure FDA0002496012840000011
M is an integer, M is more than or equal to 1 and less than or equal to M;
s7, calculating the polarization parameter k of each sub-diameter m in each diameter n according to the statistical characteristics of the polarization rate cross parametersn，m；
S8, obtaining the initial phase of the polarized antenna of each sub-diameter m in each diameter n according to the statistical characteristics of the initial phase of the polarized antenna
Figure FDA0002496012840000012
And
Figure FDA0002496012840000013
is a vertical-vertical cross-polarization initial phase,
Figure FDA0002496012840000014
for the vertical-horizontal cross polarization initial phase,
Figure FDA0002496012840000015
for the horizontal-horizontal cross-polarization initial phase,
Figure FDA0002496012840000016
is a horizontal-vertical cross polarization initial phase;
s9, establishing a 3D MIMO channel time domain expression according to the wireless communication channel principle and the multi-antenna characteristics, and obtaining the path power P of each path nnPolarization parameter k for each sub-diameter m in each diameter nn，mHorizontal emission angle αn，mVertical emission angle βn，mHorizontal reception angle gamman，mAnd vertical acceptance angle
Figure FDA0002496012840000021
And its polarization antenna initial phase
Figure FDA0002496012840000022
Determining an nth path communication time domain transmission function H received by an s antenna transmitting an u antenna of a 3D MIMO channelu，s，nAnd (t) finishing channel modeling when t is the current moment, wherein S and U are positive integers, S is more than or equal to 1 and less than or equal to S, and U is more than or equal to 1 and less than or equal to U.
2. The method of ray-based 3D MIMO channel modeling according to claim 1, wherein the method of determining the communication scenario type in step S2 is to decide the communication scenario to be L OS if the wireless signal is transmitted between the transmitting end antenna and the receiving end antenna without occlusion, and to decide the communication scenario to be N L OS otherwise.
3. The method for ray-based 3D MIMO channel modeling according to claim 2, wherein the method for determining the number N of propagation paths between the transmitting antenna and the receiving antenna according to the communication scenario type in step S3 is that if the communication scenario is N L OS, the number N of propagation paths is L according to the number of barriers and physical electrical characteristics by electromagnetic field and electromagnetic wave theory, and if the communication scenario is L OS, the number N of propagation paths in L OS scenario is L +1, i.e. a direct path is added on the basis of the propagation paths in N L OS scenario.
4. The method for ray-based 3D MIMO channel modeling according to claim 1, wherein the statistical property of the path power in step S5 is exponential distribution.
5. The method for ray-based 3D MIMO channel modeling according to claim 1, wherein the step S6 specifically includes the steps of:
s61, horizontal emission angle α according to sub-diametern，mVertical emission angle βn，mHorizontal reception angle gamman，mAnd vertical acceptance angle
Figure FDA0002496012840000023
The angle power statistics of (a) yields the horizontal emission angle αn，mVertical emission angle βn，mHorizontal reception angle gamman，mAnd vertical acceptance angle
Figure FDA0002496012840000024
Respective angular power spectral density PAS function pα(θ)、pβ(θ)、pγ(theta) and
Figure FDA0002496012840000025
wherein theta is an angle falling within the range of 0-360 DEG, and the horizontal emission angle is αn，mThe angle power statistic characteristic of (1) is Laplace distribution and horizontal receiving angle gamman，mHas an angular power statistic of 35 DEG Laplacian distribution or 360 DEG uniform distribution, vertical emission angle βn，mAnd vertical acceptance angle
Figure FDA0002496012840000031
The statistical characteristics of the angle power of the power sensor are all uniformly distributed at 30 degrees;
s62, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure FDA0002496012840000032
All angles theta that holdiWherein i is 1, 2, 1, M, and the angles are recorded as the horizontal emission angle α of the sub-pathn，m；
S63, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure FDA0002496012840000033
All angles theta that holdiWherein i is 1, 2, 1, M, and the angles are recorded as the perpendicular emission angles β of the sub-pathsn，m；
S64, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure FDA0002496012840000034
All angles theta that holdiWherein, i is 1, 2, and M, and the angles are regarded as the horizontal receiving angle γ of the sub-pathn，m；
S65, traversing the angle of 0-360 degrees to find out the angle of the energy-saving device
Figure FDA0002496012840000035
All angles theta that holdiWherein, i is 1, 2, 1, M, and the angles are perpendicular receiving angles of the sub-paths
Figure FDA0002496012840000036
6. The method for ray-based 3D MIMO channel modeling according to claim 1, wherein the statistical property of the polarization rate cross-over parameter in the step S7 is lognormal distribution, and the expression of the polarization parameter is:
Figure FDA0002496012840000037
where are random variables that follow a gaussian distribution.
7. The method for ray-based 3D MIMO channel modeling according to claim 1, wherein the statistical property of the initial phase of the polarized antennas in the step S8 is a uniform distribution from-pi to pi.
8. The method for ray-based 3D MIMO channel modeling according to claim 1, wherein the 3D MIMO channel time domain expression established according to the wireless communication channel principle and the multi-antenna characteristics in step S9 is:
Figure FDA0002496012840000041
wherein k isn，mFor each sub-diameter m of each diameter n, a polarization parameter Gtx，s，VFor power gain in the direction of vertical polarization of the s-th transmitting antenna, Gtx，s，HFor power gain in horizontal polarization direction of the s-th transmitting antenna, Grx，u，VFor the power gain in the vertical polarization direction of the u-th receiving antenna, Grx，u，HPower gain, x, for the horizontal polarization direction of the u-th receiving antennas、ysAnd zsIs the spatial coordinate, x, of the s-th transmitting antennau、yuAnd zuFor the spatial coordinates of the u-th receiving antenna,
Figure FDA0002496012840000042
is the inverse of the wavelength of the carrier wave,
Figure FDA0002496012840000043
is the velocity vector of the movable antenna end, t is the current time, [ ·]HIs a matrix conjugate transpose operation.
CN201910500886.9A 2019-06-11 2019-06-11 Ray-based 3D MIMO channel modeling method Active CN110212955B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910500886.9A CN110212955B (en) 2019-06-11 2019-06-11 Ray-based 3D MIMO channel modeling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910500886.9A CN110212955B (en) 2019-06-11 2019-06-11 Ray-based 3D MIMO channel modeling method

Publications (2)

Publication Number Publication Date
CN110212955A CN110212955A (en) 2019-09-06
CN110212955B true CN110212955B (en) 2020-07-14

Family

ID=67791924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910500886.9A Active CN110212955B (en) 2019-06-11 2019-06-11 Ray-based 3D MIMO channel modeling method

Country Status (1)

Country Link
CN (1) CN110212955B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112104579B (en) * 2020-08-25 2021-11-19 西安交通大学 Channel estimation method, system, device and storage medium based on model constraint
CN112511199B (en) * 2020-11-10 2022-04-08 河海大学 Low-complexity 3D-MIMO related channel fast generation method
CN113534204B (en) * 2021-07-16 2023-03-10 湖南北迪教育科技有限公司 Navigation multipath signal simulation method and system based on three-dimensional scene
CN114665998B (en) * 2022-03-22 2023-02-21 北京大学 Triple non-stationary wireless communication channel modeling method under space-time consistency

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103984836A (en) * 2014-05-30 2014-08-13 电子科技大学 Combined storage method for 3D-MIMO antenna polarization gains
CN105553584A (en) * 2015-12-10 2016-05-04 国网山东省电力公司烟台供电公司 3DMIMO channel modeling method
CN107302387A (en) * 2017-06-30 2017-10-27 西安电子科技大学 A kind of high-speed aircraft relays dual polarization mimo channel modeling method
CN107425895A (en) * 2017-06-21 2017-12-01 西安电子科技大学 A kind of 3D MIMO statistical channel modeling methods based on actual measurement
CN108365903A (en) * 2018-01-29 2018-08-03 哈尔滨工程大学 A kind of three-dimensional Massive mimo channels modeling method based on random scatter cluster
CN109588064A (en) * 2016-05-13 2019-04-05 瑞典爱立信有限公司 The network architecture, method and apparatus for cordless communication network

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107171757B (en) * 2017-05-22 2020-10-23 电子科技大学 MIMO wireless channel modeling method fusing intelligent antenna
EP3685468A4 (en) * 2017-09-22 2021-06-16 Duke University Enhanced mimo communication systems using reconfigurable metasurface antennas and methods of using same
US10313034B2 (en) * 2017-10-12 2019-06-04 Spirent Communications, Inc. Massive MIMO array testing using a programmable phase matrix and channel emulator
CN109802737A (en) * 2019-01-23 2019-05-24 电子科技大学 A kind of root mean square angle spread acquisition methods of 3D mimo channel modeling

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103984836A (en) * 2014-05-30 2014-08-13 电子科技大学 Combined storage method for 3D-MIMO antenna polarization gains
CN105553584A (en) * 2015-12-10 2016-05-04 国网山东省电力公司烟台供电公司 3DMIMO channel modeling method
CN109588064A (en) * 2016-05-13 2019-04-05 瑞典爱立信有限公司 The network architecture, method and apparatus for cordless communication network
CN107425895A (en) * 2017-06-21 2017-12-01 西安电子科技大学 A kind of 3D MIMO statistical channel modeling methods based on actual measurement
CN107302387A (en) * 2017-06-30 2017-10-27 西安电子科技大学 A kind of high-speed aircraft relays dual polarization mimo channel modeling method
CN108365903A (en) * 2018-01-29 2018-08-03 哈尔滨工程大学 A kind of three-dimensional Massive mimo channels modeling method based on random scatter cluster

Also Published As

Publication number Publication date
CN110212955A (en) 2019-09-06

Similar Documents

Publication Publication Date Title
CN110212955B (en) Ray-based 3D MIMO channel modeling method
US9786999B2 (en) Over-the-air test
US8880002B2 (en) Over-the-air test
CN104698430B (en) It is a kind of for carrying the high-precision angle estimating method based on virtual antenna array
CN109714120B (en) Method for simulating propagation characteristics of coupling multi-antenna indoor space fading channel
CN107171757B (en) MIMO wireless channel modeling method fusing intelligent antenna
CN108494512B (en) Millimeter wave large-scale MIMO channel simulation system and method
CN107171703B (en) Method capable of simulating propagation characteristics of fading signals in indoor multi-antenna communication system
Tushar et al. Distributed transmit beamforming based on a 3-bit feedback system
Zhou et al. Path loss model based on cluster at 28 GHz in the indoor and outdoor environments
CN114665998B (en) Triple non-stationary wireless communication channel modeling method under space-time consistency
Shaddad et al. Channel estimation for intelligent reflecting surface in 6G wireless network via deep learning technique
CN105244634A (en) Large scale MIMO antenna array dimension reduction method and system using the method
CN113965881A (en) Millimeter wave integrated communication and sensing method under shielding effect
CN107426119B (en) Cross-shaped channel estimation method, device and system for millimeter wave communication
CN113162712B (en) Multi-beam channel modeling method based on propagation diagram theory
Fedosov et al. Adaptive algorithm for MIMO-system of the wireless access for the receiver of the mobile station
CN109861776A (en) A method of the simulation multiple antennas exterior space couples fading propagation characteristic
CN115915157A (en) Urban environment terahertz wireless communication base station site selection method based on ray tracing
CN111181879B (en) Anti-interference time-varying channel equalization method based on adaptive basis expansion
Zhao et al. Analysis of a distributed MIMO channel capacity under a special scenario
CN113726390A (en) Energy efficiency maximization transmission system and method based on reconfigurable reflection super surface
CN113726389A (en) Cell-free networking system and method based on reconfigurable reflection super surface
Wang et al. Research on MIMO Channel Capacity in Complex Indoor Environment Based on Deterministic Channel Model
Vilhar et al. Efficient open-source ray-tracing methods for rural environment

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