CN114844584B - Simulation method of beam channel based on random twin cluster - Google Patents

Simulation method of beam channel based on random twin cluster Download PDF

Info

Publication number
CN114844584B
CN114844584B CN202210778010.2A CN202210778010A CN114844584B CN 114844584 B CN114844584 B CN 114844584B CN 202210778010 A CN202210778010 A CN 202210778010A CN 114844584 B CN114844584 B CN 114844584B
Authority
CN
China
Prior art keywords
scattering
cluster
clusters
originating
representing
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
CN202210778010.2A
Other languages
Chinese (zh)
Other versions
CN114844584A (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.)
Beijing Jiaotong University
Original Assignee
Beijing Jiaotong University
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 Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CN202210778010.2A priority Critical patent/CN114844584B/en
Publication of CN114844584A publication Critical patent/CN114844584A/en
Application granted granted Critical
Publication of CN114844584B publication Critical patent/CN114844584B/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
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
    • 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
    • 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/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842—Weighted combining
    • H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention provides a simulation method of a beam channel based on a random twin cluster. Determining the number of scattering point clusters at a transmitting and receiving end according to the width of a beam at the transmitting and receiving end, and initializing the positions of the scattering point clusters to obtain a beam parameter at the transmitting and receiving end; calculating beam response according to the beam parameters of the transceiving end, and calculating the channel transmission function of each link at the current moment; updating time, the positions of a transceiving end, a beam direction and a scattering cluster position, calculating a channel transmission function of each link at the updating time, calculating the survival probability of the scattering cluster in a beam range, determining the occurrence and the extinction of the scattering cluster after the updating time, counting the number of the surviving scattering clusters, determining the number of the scattering clusters based on a Poisson process, and generating a new scattering cluster if the survival number is lower than the threshold. The method of the invention considers the random mobility of the scattering cluster and the influence of the wave beam on the scattering cluster generation and extinction, and makes up for the fact that the random mobility of the scattering cluster and the influence of the wave beam on the scattering cluster generation and extinction characteristics are not considered in the conventional random channel modeling theory.

Description

Simulation method of beam channel based on random twin cluster
Technical Field
The invention relates to the technical field of wireless communication, in particular to a beam channel simulation method based on a random twin cluster.
Background
As a key technology in 5G, MIMO (Multiple Input Multiple Output) technology and beamforming technology, not only can the frequency band utilization rate be greatly improved, but also the serious path loss caused by a high frequency band, especially a millimeter wave frequency band, can be overcome, thereby attracting wide attention in academic and industrial fields. In recent years, with the continuous development of car networking, unmanned aerial vehicles, high-speed railways, urban rail transit and the like, the beamforming technology shows huge application advantages and potentials in wireless communication of high-mobility scenes, and therefore, the accurate and efficient beam channel model can provide important physical layer technical support for research and analysis of radio wave propagation characteristics, communication system optimization and the like in the high-mobility scenes.
The wireless channel model methods are mainly classified into deterministic channel modeling and stochastic channel modeling. The typical method of the deterministic channel model is a ray tracing method, which calculates parameters of multipath, such as power, angle, time delay and the like, by constructing an accurate propagation environment. Representative models in the stochastic channel modeling method are a geometric-based stochastic channel model and a non-geometric stochastic modeling method, typical non-geometric stochastic channel models include a tapped delay line model and a clustered delay line model, and the geometric-stochastic-based model assumes that scattering points are distributed on a regular geometric shape. Although some work is carried out on random channel models based on geometry to carry out research on beam channel modeling, the models do not take the mobility of scattering points into consideration, and do not take the influence of beams on scattering cluster extinction into consideration, so that the non-stationary characteristics of channels under beam channels are difficult to characterize.
Disclosure of Invention
The invention provides a beam channel simulation method based on a random twin cluster, which is used for solving the problem that the movement characteristic of a scattering cluster is not considered in the conventional random channel modeling theory.
In order to achieve the purpose, the invention adopts the following technical scheme.
A simulation method of beam channels based on random twin clusters comprises the following steps:
initializing a receiving and transmitting end beam parameter, wherein the receiving and transmitting end beam parameter comprises a beam direction and a beam width, the beam direction is determined according to the position of a receiving and transmitting end, the number of scattering point clusters of the receiving and transmitting end is determined according to the width of the receiving and transmitting end beam, the position of each scattering point cluster is initialized by adopting a Matern hard kernel Poisson cluster process, and a random speed is given to each scattering point cluster;
determining the current time when simulation starts, calculating beam response according to initialization parameters of receiving and transmitting end beams, constructing a line-of-sight link, a single-hop link and a double-hop link, calculating a transmission function corresponding to each link by adopting a propagation diagram theory, and calculating to obtain a channel transmission function of each link at the current time;
updating time, the positions of a transceiving end, a beam direction and a scattering cluster position, calculating to obtain a channel transmission function of each link at an updating moment, and calculating the survival probability of each scattering cluster according to the position and the speed of the scattering cluster at each moment;
determining the occurrence and the extinction of the scattering clusters after the updating time according to the scattering cluster survival probability, counting the number of the survival scattering clusters, determining the number of the scattering clusters at the moment based on the poisson process, and if the survival number is lower than the threshold value, generating a new scattering cluster based on the Matern hard-core poisson cluster process; the above process is repeated until the simulation deadline.
Preferably, the initializing a transceiving end beam parameter includes a beam direction and a beam width, where the beam direction is determined according to a transceiving end position, the number of scattering point clusters at the transceiving end is determined according to the width of a transceiving end beam, the position of each scattering point cluster is initialized by using a matern hard kernel poisson cluster process, and a random speed is given to each scattering point cluster, including:
the beam pointing direction of the transmitting and receiving ends is determined by the angle of the line-of-sight link, the horizontal angle of departure of the line-of-sight link
Figure 266240DEST_PATH_IMAGE001
And away from pitch
Figure 888852DEST_PATH_IMAGE002
The calculation formula is as follows:
Figure 722815DEST_PATH_IMAGE003
Figure 538368DEST_PATH_IMAGE004
Figure 483190DEST_PATH_IMAGE005
Figure 225887DEST_PATH_IMAGE006
wherein
Figure 496331DEST_PATH_IMAGE007
Respectively representing the array center coordinates of the originating and terminating ends,
Figure 793320DEST_PATH_IMAGE008
respectively representx、 y、zAxial unit vector, sign
Figure 682779DEST_PATH_IMAGE009
It is shown that the operation of the inner product of the vectors,
Figure 279983DEST_PATH_IMAGE010
representing an inverse cosine function; similar horizontal angle of arrival for line-of-sight links
Figure 721328DEST_PATH_IMAGE011
And angle of pitch of arrival
Figure 771193DEST_PATH_IMAGE012
；
The beam pointing direction of the transceiving end is calculated as follows:
Figure 57817DEST_PATH_IMAGE013
Figure 916052DEST_PATH_IMAGE014
Figure 856195DEST_PATH_IMAGE015
Figure 737563DEST_PATH_IMAGE016
wherein
Figure 618757DEST_PATH_IMAGE017
Indicating the transmit/receive beam horizontal pointing angle,
Figure 269182DEST_PATH_IMAGE018
representing the transmit/receive beam elevation pointing angle;
Figure 255592DEST_PATH_IMAGE019
indicating the horizontal angle of departure/arrival of the transceiving end line-of-sight link,
Figure 748890DEST_PATH_IMAGE020
a pitch angle representing departure/arrival of the transceiving end line-of-sight link;
Figure 642897DEST_PATH_IMAGE021
indicating the horizontal width of the originating/receiving beam,
Figure 475724DEST_PATH_IMAGE022
representing the pitch width of the originating/receiving beam;
Figure 429773DEST_PATH_IMAGE023
represents a rounding function;
determining the number of scattering point clusters at the transmitting and receiving ends according to the width of the wave beams at the transmitting and receiving ends, wherein the calculation formula is as follows:
Figure 613630DEST_PATH_IMAGE024
Figure 780169DEST_PATH_IMAGE025
wherein
Figure 264240DEST_PATH_IMAGE026
Which represents the initial moment of time of day,
Figure 123612DEST_PATH_IMAGE027
indicating the number of originating/terminating scattering clusters,
Figure 201289DEST_PATH_IMAGE028
in order to be the density factor, the density of the sample,
Figure 171519DEST_PATH_IMAGE029
representing the effective radius centered at the base station/terminal,
Figure 775676DEST_PATH_IMAGE030
represents the minimum and maximum values of the transmit/receive beam in elevation, and has
Figure 478053DEST_PATH_IMAGE031
；
Generating scattering clusters based on a Matern hard kernel Poisson cluster process, firstly generating the positions of scattering point cluster centers by adopting the Matern hard kernel point process, wherein the positions of any two scattering cluster centers need to meet the following requirements:
Figure 167660DEST_PATH_IMAGE032
Figure 348106DEST_PATH_IMAGE033
wherein
Figure 275610DEST_PATH_IMAGE034
Respectively representing the coordinates of the centers of two different scattering clusters,
Figure 7943DEST_PATH_IMAGE035
expressing the minimum interval between clusters, and generating cluster scattering points based on the Poisson point process for the generated scattering cluster center, wherein the cluster scattering points need to satisfy the following conditions:
Figure 60213DEST_PATH_IMAGE036
Figure 106666DEST_PATH_IMAGE037
wherein
Figure 956854DEST_PATH_IMAGE038
Representing scattering clustersmInner to the firstiScattering point
Figure 735454DEST_PATH_IMAGE039
The coordinates of (a) are calculated,
Figure 602916DEST_PATH_IMAGE040
representing effective radius of scattering cluster, for scattering clustermAll scattering points within the image are assigned a random velocity, denoted
Figure 249798DEST_PATH_IMAGE041
。
Preferably, the determining the current time when the simulation starts, calculating a beam response according to an initialization parameter of a beam at a transceiving end, constructing a line-of-sight link, a single-hop link, and a double-hop link, calculating a transmission function corresponding to each link by using a propagation map theory, and calculating a channel transmission function of each link at the current time includes:
determining the current time when the simulation starts, wherein the beam parameters at the transceiving end comprise beam pointing direction and beam width, and the beam response calculation formula is as follows:
Figure 558419DEST_PATH_IMAGE042
Figure 101396DEST_PATH_IMAGE043
wherein
Figure 987312DEST_PATH_IMAGE044
Any horizontal angle and pitch angle of the originating/terminating end are shown,
Figure 578831DEST_PATH_IMAGE045
maximum beam gain representing the transmit/receive beam;
the line-of-sight link refers to
Figure 69855DEST_PATH_IMAGE046
In a direct link of
Figure 580471DEST_PATH_IMAGE047
Refer to the originating endpAn antenna is arranged on the base plate, and the antenna is arranged on the base plate,
Figure 94629DEST_PATH_IMAGE048
refer to the receiving endqAn antenna; the single-hop link is represented as
Figure 83313DEST_PATH_IMAGE049
Indicating the proximity of the originating endmA plurality of scattering clusters, each scattering cluster having a plurality of scattering lines,
Figure 428844DEST_PATH_IMAGE050
indicating the vicinity of the terminating endnA scattering cluster, and have
Figure 251306DEST_PATH_IMAGE051
(ii) a The double-hop link is represented as
Figure 111815DEST_PATH_IMAGE052
；
The transmission function corresponding to each link is calculated by adopting a propagation diagram theory, wherein the transmission function expression of the direct link is as follows:
Figure 779557DEST_PATH_IMAGE053
Figure 979594DEST_PATH_IMAGE054
wherein
Figure 832012DEST_PATH_IMAGE055
Representing line-of-sight links
Figure 320762DEST_PATH_IMAGE056
The distance between the first and second electrodes,
Figure 854512DEST_PATH_IMAGE057
which is indicative of the wavelength of the light,
Figure 971373DEST_PATH_IMAGE058
in order to be a time-delay variable,
Figure 135638DEST_PATH_IMAGE059
is a variable of the time, and is,
Figure 908422DEST_PATH_IMAGE060
in order to be the speed of light,
Figure 42600DEST_PATH_IMAGE061
is a dirac function; the transmission function of the single-hop link is expressed as:
Figure 154912DEST_PATH_IMAGE063
Figure 366711DEST_PATH_IMAGE064
wherein
Figure 892371DEST_PATH_IMAGE065
Indicating a link
Figure 440027DEST_PATH_IMAGE066
The distance of (a) to (b),
Figure 321DEST_PATH_IMAGE067
representing scattering clustersnInner firstjScattering point
Figure 99864DEST_PATH_IMAGE068
，
Figure 50502DEST_PATH_IMAGE069
Representing scattering clusters
Figure 526483DEST_PATH_IMAGE070
The number of internal scattering points is greater than the total number of internal scattering points,
Figure 675705DEST_PATH_IMAGE071
indicating a link
Figure 352674DEST_PATH_IMAGE072
Gain, the expression is as follows:
Figure 384084DEST_PATH_IMAGE073
Figure 601438DEST_PATH_IMAGE074
Figure 542849DEST_PATH_IMAGE075
Figure 249774DEST_PATH_IMAGE076
wherein
Figure 643847DEST_PATH_IMAGE077
Representing scattering points
Figure 664892DEST_PATH_IMAGE078
Relative to the horizontal and pitch angles of departure from the origin,
Figure 585444DEST_PATH_IMAGE079
representing scattering points
Figure 604215DEST_PATH_IMAGE080
Relative to the horizontal and pitch angles of approach at the closeout,
Figure 610217DEST_PATH_IMAGE081
respectively represent
Figure 841479DEST_PATH_IMAGE082
Figure 819799DEST_PATH_IMAGE083
The calculation formula of the average transmission delay is as follows:
Figure 868526DEST_PATH_IMAGE084
Figure 502770DEST_PATH_IMAGE085
Figure 662356DEST_PATH_IMAGE086
Figure 495183DEST_PATH_IMAGE087
Figure 590178DEST_PATH_IMAGE088
Figure 842211DEST_PATH_IMAGE089
Figure 743170DEST_PATH_IMAGE090
Figure 368187DEST_PATH_IMAGE091
Figure 758717DEST_PATH_IMAGE092
the transfer function of the dual-hop link is expressed as:
Figure 101974DEST_PATH_IMAGE093
Figure 72204DEST_PATH_IMAGE094
wherein
Figure 676360DEST_PATH_IMAGE095
Indicating a link
Figure 378737DEST_PATH_IMAGE096
The distance of (a) to (b),
Figure 68344DEST_PATH_IMAGE097
indicating a link
Figure 983211DEST_PATH_IMAGE098
The expression is:
Figure 910716DEST_PATH_IMAGE099
Figure 643048DEST_PATH_IMAGE100
Figure 960897DEST_PATH_IMAGE101
wherein
Figure 804088DEST_PATH_IMAGE102
Is a scattering attenuation factor, the value of which depends on the material of the scattering cluster;
the channel transfer function is the sum of the transfer functions of the links, and is expressed as:
Figure 258203DEST_PATH_IMAGE103
Figure 364699DEST_PATH_IMAGE104
。
preferably, the updating time, the positions of the transceiving end, the beam direction, and the positions of the scattering clusters are calculated to obtain a channel transfer function of each link at the updating time, and the survival probability of the scattering clusters in the beam range is calculated according to the channel transfer function of each link at each time, including:
the receiving and transmitting end position is updated as follows:
Figure 294478DEST_PATH_IMAGE105
Figure 144623DEST_PATH_IMAGE106
wherein
Figure 718823DEST_PATH_IMAGE107
Refer to
Figure 996221DEST_PATH_IMAGE108
At the moment of timetIn the position of (a) in the first,
Figure 147717DEST_PATH_IMAGE109
refer to
Figure 473656DEST_PATH_IMAGE110
At the moment of time
Figure 964680DEST_PATH_IMAGE111
In the position of (a) in the first,
Figure 209716DEST_PATH_IMAGE112
representing the velocity vectors of the originating and the terminating ends; the scattering cluster
Figure 723874DEST_PATH_IMAGE113
The location update is as follows:
Figure 972279DEST_PATH_IMAGE114
Figure 317810DEST_PATH_IMAGE115
wherein
Figure 874693DEST_PATH_IMAGE116
Refers to originating/terminating scattering clusters
Figure 781DEST_PATH_IMAGE117
Centered at the momenttIn the position of (a) in the first,
Figure 402943DEST_PATH_IMAGE118
the scattering cluster center of the transmitting end/receiving end is in the time
Figure 665297DEST_PATH_IMAGE119
In the position of (a) in the first,
Figure 393082DEST_PATH_IMAGE120
representing originating/receiving scattering clusters
Figure 944149DEST_PATH_IMAGE117
At the moment of timetThe scattering cluster is in a random walk state, the scattering cluster
Figure 540215DEST_PATH_IMAGE117
The speed of (2) is updated as follows:
Figure 860338DEST_PATH_IMAGE121
Figure 493445DEST_PATH_IMAGE122
wherein
Figure 594125DEST_PATH_IMAGE123
Transmitting/receiving end scattering cluster
Figure 931565DEST_PATH_IMAGE117
At the moment of time
Figure 778299DEST_PATH_IMAGE119
The instantaneous speed of the vehicle (c) is,
Figure 238099DEST_PATH_IMAGE124
is a factor of the inertia, and is,
Figure 170283DEST_PATH_IMAGE125
representing originating/receiving scattering clusters
Figure 311414DEST_PATH_IMAGE117
At the moment of time
Figure 402867DEST_PATH_IMAGE119
Is expressed as follows:
Figure 643355DEST_PATH_IMAGE126
Figure 453048DEST_PATH_IMAGE127
wherein
Figure 804395DEST_PATH_IMAGE128
Representing originating/receiving scattering clusters
Figure 953617DEST_PATH_IMAGE117
At the moment of timetIn order to maintain track continuity, the phase is updated as follows:
Figure 224061DEST_PATH_IMAGE129
Figure 396416DEST_PATH_IMAGE130
wherein
Figure 676088DEST_PATH_IMAGE131
Which represents the maximum allowed phase shift and,
Figure 685676DEST_PATH_IMAGE132
to be uniformly distributed in
Figure 267967DEST_PATH_IMAGE133
Random variable in between;
the updated beam pointing angle is specifically as follows: according to the updated coordinates of the transmitting and receiving ends, the departure angle and the arrival angle of the line-of-sight link are calculated, and the beam direction is updated based on a formula (1-4);
the filtering of the scattering clusters outside the beam range according to the beam direction specifically includes: for single hop links
Figure 52252DEST_PATH_IMAGE134
And
Figure 745401DEST_PATH_IMAGE135
determining the center of the scattering cluster
Figure 869215DEST_PATH_IMAGE136
Whether the scattering cluster is positioned in the wave beam range of the transmitting and receiving ends at the same time or not, if not, removing the scattering cluster;
the scattering cluster survival probability calculation formula is as follows:
Figure 747041DEST_PATH_IMAGE137
Figure 628410DEST_PATH_IMAGE138
wherein
Figure 249884DEST_PATH_IMAGE139
Scattering cluster for representing originating/receiving end
Figure 228204DEST_PATH_IMAGE140
The parameters of the birth or death over time,
Figure 417877DEST_PATH_IMAGE141
scattering cluster for representing originating/receiving end
Figure 911175DEST_PATH_IMAGE140
A birth-to-death parameter on the beam;
Figure 946127DEST_PATH_IMAGE142
which represents the reference distance, is a fixed constant,
Figure 778954DEST_PATH_IMAGE143
the reference angle is represented as a fixed constant;
Figure 998583DEST_PATH_IMAGE144
indicating a link
Figure 854544DEST_PATH_IMAGE145
The angle between the beam direction and the transmitting end/receiving end.
Preferably, the occurrence and the extinction of the scattering clusters after the updating time are determined according to the scattering cluster survival probability, the number of the survival scattering clusters is counted, the number of the scattering clusters at the moment is determined based on the poisson process, and if the survival number is lower than the threshold value, a new scattering cluster is generated based on the Matern hard-core poisson cluster process; repeating the above process until the simulation deadline time, comprising:
randomly generating a random number uniformly distributed among (0,1)
Figure 817820DEST_PATH_IMAGE146
If at all
Figure 770733DEST_PATH_IMAGE147
Then the scattering cluster of the transmitting/receiving end
Figure 36629DEST_PATH_IMAGE140
At the moment of time
Figure 504520DEST_PATH_IMAGE119
Survival; if it is
Figure 881274DEST_PATH_IMAGE148
Then the scattering cluster of the transmitting/receiving end
Figure 688693DEST_PATH_IMAGE140
At the moment of time
Figure 718966DEST_PATH_IMAGE119
Eliminating;
to be provided with
Figure 408573DEST_PATH_IMAGE149
Generating a random variable of a slave Poisson process for a mean value
Figure 589019DEST_PATH_IMAGE150
And for random variables
Figure 313261DEST_PATH_IMAGE150
Rounding is carried out if the number of the scattering clusters surviving at the transmitting end/the receiving end is less than
Figure 186540DEST_PATH_IMAGE150
Generating a new scattering cluster based on a Matern hard-core Poisson cluster process; all steps corresponding to equations (1-25) are repeated until the simulation deadline.
According to the technical scheme provided by the embodiment of the invention, the method considers the random mobility of the scattering cluster and the influence of the wave beam on the scattering cluster generation and extinction, makes up the problem that the random mobility of the scattering cluster and the influence of the wave beam on the scattering cluster generation and extinction characteristics are not considered in the conventional random channel modeling theory, and provides a new method for random channel simulation in a wave beam mode.
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 process flow diagram of a simulation method of a beam channel based on a random twin cluster according to the present invention.
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 functions 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.
The processing flow of the simulation method based on the random twin cluster beam channel provided by the invention is shown in fig. 1, and comprises the following processing steps:
step S1, beam pointing is determined according to the position of a transmitting and receiving end of wireless communication, the number of scattering point clusters of the transmitting and receiving end is determined according to the width of the beam of the transmitting and receiving end, the position of the scattering point clusters is initialized by adopting a Matern hard kernel Poisson cluster process, and random speed is given to each scattering point cluster.
The beam pointing direction of the transmitting and receiving end is determined by the angle of a line-of-sight link, and the horizontal angle of the line-of-sight link
Figure 115442DEST_PATH_IMAGE151
And angle of departure pitch
Figure 427474DEST_PATH_IMAGE152
The calculation formula of (a) is as follows:
Figure 147169DEST_PATH_IMAGE153
Figure 50403DEST_PATH_IMAGE154
Figure 121127DEST_PATH_IMAGE155
Figure 33588DEST_PATH_IMAGE156
wherein
Figure 342210DEST_PATH_IMAGE157
Respectively representing the array center coordinates of the originating terminal and the terminating terminal,
Figure 619607DEST_PATH_IMAGE158
respectively representx、y、zAxial unit vector, sign
Figure 771103DEST_PATH_IMAGE159
It is shown that the operation of the inner product of the vectors,
Figure 97042DEST_PATH_IMAGE160
representing an inverse cosine function; similar horizontal angle of arrival for line-of-sight links
Figure 650383DEST_PATH_IMAGE161
And angle of pitch of arrival
Figure 98682DEST_PATH_IMAGE162
；
The beam pointing is calculated as follows:
Figure 612840DEST_PATH_IMAGE163
Figure 601525DEST_PATH_IMAGE164
Figure 884738DEST_PATH_IMAGE165
Figure 503939DEST_PATH_IMAGE166
wherein
Figure 364447DEST_PATH_IMAGE167
Indicating the originating/terminating beam horizontal pointing angle,
Figure 297768DEST_PATH_IMAGE168
representing the transmit/receive beam elevation pointing angle;
Figure 560122DEST_PATH_IMAGE169
indicating the horizontal angle of departure/arrival of the transceiving end line-of-sight link,
Figure 22328DEST_PATH_IMAGE170
a pitch angle representing departure/arrival of the transceiving end line-of-sight link;
Figure 573395DEST_PATH_IMAGE171
indicating the horizontal width of the originating/terminating beam,
Figure 169461DEST_PATH_IMAGE172
represents the pitch width of the transmit/receive beam;
Figure 427267DEST_PATH_IMAGE173
represents a rounding function;
determining the number of scattering point clusters at the transmitting and receiving ends according to the width of the wave beams at the transmitting and receiving ends, wherein the calculation formula is as follows:
Figure 185007DEST_PATH_IMAGE174
Figure 895474DEST_PATH_IMAGE175
wherein
Figure 238774DEST_PATH_IMAGE176
Which represents the initial moment of time of day,
Figure 210141DEST_PATH_IMAGE177
indicating the number of originating/terminating scattering clusters,
Figure 545308DEST_PATH_IMAGE178
in order to be the density factor, the density of the sample,
Figure 602125DEST_PATH_IMAGE179
representing an effective radius centered at the base station/terminal,
Figure 743257DEST_PATH_IMAGE180
represents the minimum and maximum values of the transmit/receive beam in elevation, and has
Figure 444497DEST_PATH_IMAGE181
。
Generating scattering clusters based on a Matern hard kernel Poisson cluster process, firstly generating the positions of scattering point cluster centers by adopting the Matern hard kernel point process, wherein the positions of any two scattering cluster centers need to meet the following requirements:
Figure 75198DEST_PATH_IMAGE182
Figure 494678DEST_PATH_IMAGE183
wherein
Figure 173921DEST_PATH_IMAGE184
And
Figure 119880DEST_PATH_IMAGE185
respectively representing the coordinates of the centers of two different scattering clusters,
Figure 796849DEST_PATH_IMAGE186
indicating the minimum spacing between clusters. For the generated scattering cluster center, generating an in-cluster scattering point based on a poisson point process, wherein the in-cluster scattering point needs to satisfy the following conditions:
Figure 297101DEST_PATH_IMAGE187
Figure 576772DEST_PATH_IMAGE188
wherein
Figure 252604DEST_PATH_IMAGE189
Representing scattering clustersmInner to the firstiScattering point
Figure 21846DEST_PATH_IMAGE190
Is determined by the coordinate of (a) in the space,
Figure 743815DEST_PATH_IMAGE191
representing the effective radius of the scattering cluster. For scattering clustermAll scattering points within the image are assigned a random velocity, denoted
Figure 436964DEST_PATH_IMAGE192
。
And S2, calculating beam response according to the beam parameters of the transmitting and receiving ends, constructing a line-of-sight link, a single-hop link and a double-hop link, and calculating a transmission function corresponding to each link by adopting a propagation diagram theory so as to obtain a channel transmission function at the moment.
The receiving and transmitting end beam parameters refer to beam pointing direction and beam width, and the beam response calculation formula is as follows:
Figure 623095DEST_PATH_IMAGE193
Figure 110708DEST_PATH_IMAGE194
wherein
Figure 319972DEST_PATH_IMAGE195
And
Figure 207026DEST_PATH_IMAGE196
represents any horizontal angle and pitch angle of the transmitting end/receiving end,
Figure 857450DEST_PATH_IMAGE197
representing the maximum beam gain of the originating/receiving beam.
The line-of-sight link refers to
Figure 109440DEST_PATH_IMAGE198
In a direct link of
Figure 331299DEST_PATH_IMAGE199
Refer to the originating endpAn antenna is arranged on the base plate, and the antenna is arranged on the base plate,
Figure 366251DEST_PATH_IMAGE200
refer to the receiving endqAn antenna; the single-hop link is represented as
Figure 526974DEST_PATH_IMAGE201
，
Figure 215445DEST_PATH_IMAGE202
Indicating the proximity of the originating endmA plurality of scattering clusters, each scattering cluster having a plurality of scattering lines,
Figure 399301DEST_PATH_IMAGE203
indicating the vicinity of the receiving endnA scattering cluster, and has
Figure 237944DEST_PATH_IMAGE204
(ii) a The double-hop link is represented as
Figure 987595DEST_PATH_IMAGE205
。
The transmission function corresponding to each link is calculated by adopting a propagation diagram theory, wherein the transmission function expression of the direct link is as follows:
Figure 581387DEST_PATH_IMAGE206
Figure 659064DEST_PATH_IMAGE207
wherein
Figure 426032DEST_PATH_IMAGE208
Representing line-of-sight links
Figure 905555DEST_PATH_IMAGE209
The distance between the first and second electrodes is less than the predetermined distance,
Figure 670249DEST_PATH_IMAGE210
which represents the wavelength of the light emitted by the light source,
Figure 625435DEST_PATH_IMAGE211
in order to be a time-delay variable,
Figure 805881DEST_PATH_IMAGE212
is a variable of the time, and is,
Figure 264544DEST_PATH_IMAGE213
in order to be the speed of light,
Figure 200139DEST_PATH_IMAGE214
is a dirac function; the single-hop link transfer function can be expressed as:
Figure 252409DEST_PATH_IMAGE215
Figure 626758DEST_PATH_IMAGE216
wherein
Figure 815294DEST_PATH_IMAGE217
Indicating a link
Figure 921790DEST_PATH_IMAGE218
The distance of (a) to (b),
Figure 585990DEST_PATH_IMAGE219
representing scattering clustersnInner firstjScattering point
Figure 373817DEST_PATH_IMAGE220
，
Figure 807073DEST_PATH_IMAGE221
Representing scattering clusters
Figure 756574DEST_PATH_IMAGE222
The number of internal scattering points is greater than the total number of internal scattering points,
Figure 111332DEST_PATH_IMAGE223
and
Figure 302185DEST_PATH_IMAGE224
indicating a link
Figure 730892DEST_PATH_IMAGE225
Gain, the expression is as follows:
Figure 241508DEST_PATH_IMAGE226
Figure 817983DEST_PATH_IMAGE227
Figure 682034DEST_PATH_IMAGE228
Figure 355461DEST_PATH_IMAGE229
wherein
Figure 646765DEST_PATH_IMAGE230
Representing scattering points
Figure 710536DEST_PATH_IMAGE231
Relative to the horizontal and pitch angles of departure from the origin,
Figure 502911DEST_PATH_IMAGE232
representing scattering points
Figure 906211DEST_PATH_IMAGE233
Relative to arrival at the receiving endThe horizontal angle and the pitch angle of arrival,
Figure 493050DEST_PATH_IMAGE234
respectively represent
Figure 716221DEST_PATH_IMAGE235
Figure 249970DEST_PATH_IMAGE236
The calculation formula of the average transmission delay is as follows:
Figure 366831DEST_PATH_IMAGE237
Figure 265517DEST_PATH_IMAGE238
Figure 366197DEST_PATH_IMAGE239
Figure 438058DEST_PATH_IMAGE240
Figure 550370DEST_PATH_IMAGE241
Figure 744591DEST_PATH_IMAGE242
Figure 942355DEST_PATH_IMAGE243
Figure 614644DEST_PATH_IMAGE244
Figure 174939DEST_PATH_IMAGE245
the transmission function of the double-hop link can be expressed as:
Figure 415427DEST_PATH_IMAGE246
Figure 162803DEST_PATH_IMAGE247
wherein
Figure 904363DEST_PATH_IMAGE248
Indicating a link
Figure 725689DEST_PATH_IMAGE249
The distance of (a) to (b),
Figure 279290DEST_PATH_IMAGE250
indicating a link
Figure 513963DEST_PATH_IMAGE251
The expression is:
Figure 403421DEST_PATH_IMAGE252
Figure 735045DEST_PATH_IMAGE253
wherein
Figure 317336DEST_PATH_IMAGE254
The value of the scattering attenuation factor depends on the material of the scattering cluster.
The channel transfer function at the time is calculated as:
Figure 39305DEST_PATH_IMAGE255
Figure 122667DEST_PATH_IMAGE256
s3, updating time, updating the position of the transceiver end, the beam direction and the position of a scattering cluster, filtering the scattering cluster outside the beam range according to the beam direction, and calculating the survival probability of the scattering cluster in the beam range;
the receiving and transmitting end position is updated as follows:
Figure 653006DEST_PATH_IMAGE257
Figure 734094DEST_PATH_IMAGE258
wherein
Figure 740096DEST_PATH_IMAGE259
Refer to
Figure 236937DEST_PATH_IMAGE260
At the moment of timetIn the position of (a) in the first,
Figure 277574DEST_PATH_IMAGE261
refer to
Figure 936088DEST_PATH_IMAGE260
At the moment of time
Figure 632649DEST_PATH_IMAGE262
In the position of (a) in the first,
Figure 792235DEST_PATH_IMAGE263
representing the velocity vectors of the originating terminal and the terminating terminal; the scattering cluster
Figure 828324DEST_PATH_IMAGE264
The location update is as follows:
Figure 782374DEST_PATH_IMAGE265
Figure 966230DEST_PATH_IMAGE266
wherein
Figure 804873DEST_PATH_IMAGE267
Refers to originating/terminating scattering clusters
Figure 554523DEST_PATH_IMAGE268
Centered at the momenttIn the position of (a) in the first,
Figure 820420DEST_PATH_IMAGE269
the scattering cluster center of the transmitting end/receiving end is in the time
Figure 225993DEST_PATH_IMAGE270
In the position of (a) in the first,
Figure 992961DEST_PATH_IMAGE271
representing originating/receiving scattering clusters
Figure 738063DEST_PATH_IMAGE272
At the time of daytThe scattering cluster is in a random walk state, the scattering cluster
Figure 299494DEST_PATH_IMAGE272
The speed of (2) is updated as follows:
Figure 864468DEST_PATH_IMAGE273
Figure 378669DEST_PATH_IMAGE274
wherein
Figure 102911DEST_PATH_IMAGE275
Transmitting/receiving end scattering cluster
Figure 710610DEST_PATH_IMAGE272
At the time of day
Figure 887514DEST_PATH_IMAGE270
The instantaneous speed of the vehicle (c) is,
Figure 933967DEST_PATH_IMAGE276
is a factor of the inertia, and is,
Figure 388082DEST_PATH_IMAGE277
representing originating/receiving scattering clusters
Figure 556895DEST_PATH_IMAGE272
At the moment of time
Figure 96461DEST_PATH_IMAGE270
Is expressed as follows:
Figure 946605DEST_PATH_IMAGE278
Figure 645440DEST_PATH_IMAGE279
wherein
Figure 391679DEST_PATH_IMAGE280
Representing originating/receiving scattering clusters
Figure 746437DEST_PATH_IMAGE272
At the moment of timetIn order to maintain track continuity, the phase is updated as follows:
Figure 197010DEST_PATH_IMAGE281
Figure 625717DEST_PATH_IMAGE282
wherein
Figure 870754DEST_PATH_IMAGE283
Which represents the maximum allowed phase shift and,
Figure 978387DEST_PATH_IMAGE284
is uniformly distributed in
Figure 842438DEST_PATH_IMAGE285
Random variable in between.
The updating of the beam pointing angle is specifically as follows: according to the updated coordinates of the transmitting and receiving ends, the departure angle and the arrival angle of the line-of-sight link are calculated, and the beam direction is updated based on a formula (1-4);
the method for filtering the scattering clusters outside the beam range according to the beam direction specifically comprises the following steps: for single hop links
Figure 187969DEST_PATH_IMAGE286
And
Figure 869486DEST_PATH_IMAGE287
determining the center of the scattering cluster
Figure 870940DEST_PATH_IMAGE288
Whether the scattering cluster is positioned in the wave beam range of the transmitting and receiving ends at the same time or not, if not, removing the scattering cluster;
the scattering cluster survival probability calculation formula is as follows:
Figure 866578DEST_PATH_IMAGE289
Figure 863352DEST_PATH_IMAGE290
wherein
Figure 591137DEST_PATH_IMAGE291
Scattering cluster for representing originating/receiving end
Figure 198662DEST_PATH_IMAGE272
The parameters of the occurrence and extinction over time,
Figure 404515DEST_PATH_IMAGE292
scattering cluster for representing originating/receiving end
Figure 724638DEST_PATH_IMAGE272
A birth-to-death parameter on the beam;
Figure DEST_PATH_IMAGE293
which represents the reference distance, is a fixed constant,
Figure 216799DEST_PATH_IMAGE294
the reference angle is represented as a fixed constant;
Figure 51900DEST_PATH_IMAGE295
indicating a link
Figure 61444DEST_PATH_IMAGE296
The angle between the beam direction and the transmitting end/receiving end.
S4, determining the occurrence and extinction of the scattering clusters after the updating time according to the survival probability, counting the number of the survival scattering clusters, determining the number of the scattering clusters at the moment based on the poisson process, and generating new scattering clusters based on the Matern hard-core poisson cluster process if the survival number is lower than the threshold; the above process is repeated until the simulation cutoff time.
The method for determining the occurrence and the extinction of the scattering clusters after the updating time according to the survival probability comprises the following steps: randomly generating a random number uniformly distributed among (0,1)
Figure 501653DEST_PATH_IMAGE297
If, if
Figure 695874DEST_PATH_IMAGE298
Then the scattering cluster of the transmitting/receiving end
Figure 893637DEST_PATH_IMAGE272
At the moment of time
Figure 565927DEST_PATH_IMAGE270
Survival, otherwise, death; to be provided with
Figure 860642DEST_PATH_IMAGE299
Generating a random variable from the poisson process for the mean
Figure DEST_PATH_IMAGE300
And rounding the scattering clusters if the number of the scattering clusters surviving at the transmitting end/receiving end is less than
Figure 694606DEST_PATH_IMAGE301
Generating a new scattering cluster based on a Matern hard-core Poisson cluster process; and (4) repeating all the steps corresponding to the formulas (1-25) until the simulation deadline time, and obtaining the channel transmission function of the link after the update time.
In summary, the method of the embodiment of the invention makes up for the fact that the influence of random mobility of the scattering cluster and the generation and extinction characteristics of the beam on the scattering cluster is not considered in the existing random channel modeling theory, and provides a new method for random channel modeling in a beam mode.
Compared with the traditional random geometric channel modeling method, the random mobility of the scattering cluster is considered on one hand, and the influence of the beam on channel simulation is considered on the other hand, so that the method is more suitable for simulating the random channel in a multi-antenna beam forming mode moving scene.
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 (3)

1. A simulation method of beam channels based on random twin clusters is characterized by comprising the following steps:
initializing a transmitting end wave beam parameter and a receiving end wave beam parameter, wherein the transmitting end wave beam parameter and the receiving end wave beam parameter comprise wave beam direction and wave beam width, the wave beam direction is determined according to the transmitting end position and the receiving end position, the number of transmitting end scattering clusters and the receiving end scattering clusters is determined according to the transmitting end wave beam width and the receiving end wave beam width, the position of the scattering clusters is initialized by adopting a Matern hard core point and Poisson point process, and random speed is given to each scattering cluster;
step 2, determining the initial moment of the simulation start
Figure DEST_PATH_IMAGE001
According to the beginning of the originating end and the receiving endThe method comprises the steps of calculating beam response by using an initialized beam parameter, constructing a line-of-sight link, a single-hop link and a double-hop link, calculating transmission functions corresponding to all links by adopting a propagation diagram theory, and calculating the sum of the transmission functions corresponding to all links as an initialized channel transmission function;
updating time, updating positions of the transmitting end and the receiving end, beam direction and speed and positions of scattering clusters, filtering the scattering clusters outside intersection of the transmitting end beam range and the receiving end beam range according to the transmitting end and the receiving end beam direction, calculating to obtain a channel transmission function of the current updating time, and calculating survival probability of each scattering cluster according to the position and speed of the scattering cluster of the current updating time;
determining the occurrence and the extinction of the scattering clusters at the current updating time according to the scattering cluster occurrence probability, counting the number of the survival scattering clusters, calculating the number of the scattering clusters at the current updating time based on the Poisson process and setting the number as a threshold, and if the number of the survival scattering clusters is lower than the threshold, generating new scattering clusters based on the Matern hard nuclear point and the Poisson point process; repeating the processing processes corresponding to the steps 3-4 until the simulation deadline;
the step 2 comprises the following steps:
determining an initial time at which simulation begins
Figure DEST_PATH_IMAGE002
The beam parameters of the transmitting end and the receiving end comprise beam pointing direction and beam width, and the beam response calculation formula is as follows:
Figure DEST_PATH_IMAGE003
Figure DEST_PATH_IMAGE004
wherein
Figure DEST_PATH_IMAGE005
Indicating any horizontal angle of the originating/terminating end,
Figure DEST_PATH_IMAGE006
indicating any pitch angle of the originating/terminating end,
Figure DEST_PATH_IMAGE007
represents the maximum beam gain of the originating/receiving beam,
Figure DEST_PATH_IMAGE008
representing the horizontal pointing angle of the originating/terminating beam,
Figure DEST_PATH_IMAGE009
indicating the elevation pointing angle of the originating/receiving beam,
Figure DEST_PATH_IMAGE010
indicating the horizontal width of the originating/receiving beam,
Figure DEST_PATH_IMAGE011
represents the pitch width of the transmit/receive beam;
the line-of-sight link refers to
Figure DEST_PATH_IMAGE012
In a direct link of
Figure DEST_PATH_IMAGE013
Refer to the originating endpAn antenna is arranged on the base plate, and the antenna is arranged on the base plate,
Figure DEST_PATH_IMAGE014
refer to the receiving endqAn antenna; the single-hop link is represented as
Figure DEST_PATH_IMAGE015
，
Figure DEST_PATH_IMAGE016
Indicating the proximity of the originating endmA plurality of scattering clusters, each scattering cluster having a plurality of scattering lines,
Figure DEST_PATH_IMAGE017
indicating the vicinity of the receiving endnA scattering cluster, and have
Figure DEST_PATH_IMAGE018
，
Figure DEST_PATH_IMAGE019
Indicating the number of scattering clusters near the originating end,
Figure DEST_PATH_IMAGE020
the number of scattering clusters near the receiving end is shown, and the double-hop link is shown as
Figure DEST_PATH_IMAGE021
；
The transmission function corresponding to each link is calculated by adopting a propagation diagram theory, wherein the transmission function expression of the direct link is as follows:
Figure DEST_PATH_IMAGE022
Figure DEST_PATH_IMAGE023
wherein
Figure DEST_PATH_IMAGE024
Is the number of the imaginary numbers,
Figure DEST_PATH_IMAGE025
representing line-of-sight links
Figure DEST_PATH_IMAGE026
The distance of (a) to (b),
Figure DEST_PATH_IMAGE027
which represents the wavelength of the light emitted by the light source,
Figure DEST_PATH_IMAGE028
in order to be a time-delay variable,
Figure DEST_PATH_IMAGE030
is a variable of the time, and is,
Figure DEST_PATH_IMAGE031
in order to be the speed of light,
Figure DEST_PATH_IMAGE032
is a dirac function; the transmission function of the single-hop link is expressed as:
Figure DEST_PATH_IMAGE033
Figure DEST_PATH_IMAGE034
Figure DEST_PATH_IMAGE035
wherein
Figure DEST_PATH_IMAGE036
Indicating a link
Figure DEST_PATH_IMAGE037
The distance of (a) to (b),
Figure DEST_PATH_IMAGE038
representing scattering clusters
Figure DEST_PATH_IMAGE039
Inner firstiThe number of scattering points is such that,
Figure DEST_PATH_IMAGE040
indicating a link
Figure DEST_PATH_IMAGE041
The distance of (a) to (b),
Figure DEST_PATH_IMAGE042
indicating a link
Figure DEST_PATH_IMAGE043
The distance of (a) to (b),
Figure DEST_PATH_IMAGE044
indicating a link
Figure DEST_PATH_IMAGE045
The distance of (a) to (b),
Figure DEST_PATH_IMAGE046
representing scattering clusters
Figure DEST_PATH_IMAGE047
Inner firstjThe number of the scattering points is,
Figure DEST_PATH_IMAGE048
representing scattering clusters
Figure DEST_PATH_IMAGE049
The number of internal scattering points is greater than the total number of internal scattering points,
Figure DEST_PATH_IMAGE050
indicating a link
Figure DEST_PATH_IMAGE051
The gain of (a) is obtained by the gain of (b),
Figure DEST_PATH_IMAGE052
indicating a link
Figure DEST_PATH_IMAGE053
The gain of (a) is obtained,
Figure DEST_PATH_IMAGE054
the expression is as follows:
Figure DEST_PATH_IMAGE055
Figure DEST_PATH_IMAGE056
Figure DEST_PATH_IMAGE057
Figure DEST_PATH_IMAGE058
wherein
Figure DEST_PATH_IMAGE059
Representing scattering points
Figure DEST_PATH_IMAGE060
Relative to the horizontal and pitch angles of departure from the origin,
Figure DEST_PATH_IMAGE061
representing scattering points
Figure DEST_PATH_IMAGE062
Relative to the horizontal and pitch angles of approach at the closeout,
Figure DEST_PATH_IMAGE063
respectively representing links
Figure DEST_PATH_IMAGE064
Figure DEST_PATH_IMAGE065
The calculation formula of the average transmission delay is as follows:
Figure DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE067
Figure DEST_PATH_IMAGE068
Figure DEST_PATH_IMAGE069
Figure DEST_PATH_IMAGE070
the calculation formula is as follows:
Figure DEST_PATH_IMAGE071
Figure DEST_PATH_IMAGE072
Figure DEST_PATH_IMAGE073
Figure DEST_PATH_IMAGE074
the transfer function of the dual-hop link is expressed as:
Figure DEST_PATH_IMAGE075
Figure DEST_PATH_IMAGE076
wherein
Figure DEST_PATH_IMAGE077
Indicating a link
Figure DEST_PATH_IMAGE078
The distance of (a) to (b),
Figure DEST_PATH_IMAGE079
indicating a link
Figure DEST_PATH_IMAGE080
The expression is:
Figure DEST_PATH_IMAGE081
Figure DEST_PATH_IMAGE082
wherein
Figure DEST_PATH_IMAGE083
Is a scattering attenuation factor, the value of which depends on the material of the scattering cluster;
the channel transfer function is the sum of the transfer functions of the links, and is expressed as:
Figure DEST_PATH_IMAGE084
Figure DEST_PATH_IMAGE085
the step 4 comprises the following steps:
randomly generating a random number uniformly distributed among (0,1)
Figure DEST_PATH_IMAGE086
If, if
Figure DEST_PATH_IMAGE087
In which
Figure DEST_PATH_IMAGE088
Representing originating scattering clusters
Figure DEST_PATH_IMAGE089
End scattering cluster
Figure DEST_PATH_IMAGE090
At the time of day
Figure DEST_PATH_IMAGE091
Probability of survival, the originating scattering cluster
Figure DEST_PATH_IMAGE092
End scattering cluster
Figure DEST_PATH_IMAGE093
At the moment of time
Figure DEST_PATH_IMAGE094
Survival; if it is
Figure DEST_PATH_IMAGE095
Then the originating scattering cluster
Figure DEST_PATH_IMAGE096
End scattering cluster
Figure DEST_PATH_IMAGE097
At the moment of time
Figure DEST_PATH_IMAGE098
Eliminating;
the number of scattering clusters at the current update time is calculated based on the poisson process and is set as a threshold, specifically: to be provided with
Figure DEST_PATH_IMAGE099
Generating a random variable from the originating/receiving end of the poisson process for the mean
Figure DEST_PATH_IMAGE100
Wherein
Figure DEST_PATH_IMAGE101
Representing the initial time and for random variables
Figure DEST_PATH_IMAGE102
Rounding, using the rounded value as the set threshold value,
if the number of the scattering clusters which survive at the transmitting end/the receiving end is smaller than a threshold value, generating a new scattering cluster based on Matern hard core point and Poisson point processes; and repeating the processing procedures of the steps 3 to 4 until the simulation deadline.
2. The method of claim 1, wherein step 1 comprises:
the beam pointing direction of the originating and receiving end is determined by the angle of the line-of-sight link, which is away from the horizontal
Figure DEST_PATH_IMAGE103
And an angle of departure pitch
Figure DEST_PATH_IMAGE104
The calculation formula of (a) is as follows:
Figure DEST_PATH_IMAGE105
Figure DEST_PATH_IMAGE106
Figure DEST_PATH_IMAGE107
Figure DEST_PATH_IMAGE108
wherein
Figure DEST_PATH_IMAGE109
Respectively representing the center coordinates of the antenna array at the transmitting end and the receiving end,
Figure DEST_PATH_IMAGE110
respectively representxAndzaxial unit vector, sign
Figure DEST_PATH_IMAGE111
It is shown that the operation of the inner product of the vectors,
Figure DEST_PATH_IMAGE112
representing an inverse cosine function; solving for the horizontal angle of arrival of the line-of-sight link
Figure DEST_PATH_IMAGE113
And angle of pitch of arrival
Figure DEST_PATH_IMAGE114
；
The beam pointing direction of the originating/receiving end is calculated as follows:
Figure DEST_PATH_IMAGE115
Figure DEST_PATH_IMAGE116
Figure DEST_PATH_IMAGE117
Figure DEST_PATH_IMAGE118
wherein
Figure DEST_PATH_IMAGE119
Indicating the originating/terminating beam horizontal pointing angle,
Figure DEST_PATH_IMAGE120
indicating the transmit/receive beam elevation pointing angle,
Figure DEST_PATH_IMAGE121
representing the horizontal angle of arrival at the originating/terminating end of the line-of-sight link,
Figure DEST_PATH_IMAGE122
representing the elevation angle of arrival of elevation departure/convergence of the line-of-sight link origination;
Figure DEST_PATH_IMAGE123
indicating the horizontal width of the originating/receiving beam,
Figure DEST_PATH_IMAGE124
represents the pitch width of the transmit/receive beam;
Figure DEST_PATH_IMAGE125
represents a rounding function;
determining the number of the transmitting end/receiving end scattering point clusters according to the transmitting end/receiving end beam width, wherein the calculation formula is as follows:
Figure DEST_PATH_IMAGE126
Figure DEST_PATH_IMAGE127
wherein
Figure DEST_PATH_IMAGE128
Which represents the initial moment of time of day,
Figure DEST_PATH_IMAGE129
represent
Figure DEST_PATH_IMAGE130
The number of time-to-make/receive scattering clusters,
Figure DEST_PATH_IMAGE131
in order to be the density factor, the density of the sample,
Figure DEST_PATH_IMAGE132
indicating an effective radius centered on the originating/terminating end,
Figure DEST_PATH_IMAGE133
represents the minimum and maximum values of the transmit/receive beam in elevation, and has
Figure DEST_PATH_IMAGE134
；
Generating a scattering cluster based on a Matern hard core point process and a Poisson cluster process, firstly generating the center position of the scattering cluster by adopting the Matern hard core point process, wherein the center positions of any two scattering clusters need to meet the following requirements:
Figure DEST_PATH_IMAGE135
Figure DEST_PATH_IMAGE136
wherein
Figure DEST_PATH_IMAGE137
Respectively representing the coordinates of the centers of two different scattering clusters,
Figure DEST_PATH_IMAGE138
representing the minimum interval between the scattering clusters, and generating scattering points in the scattering clusters based on the Poisson point process for the generated scattering cluster center, wherein the scattering points in the scattering clusters need to meet the requirement：
Figure DEST_PATH_IMAGE139
Figure DEST_PATH_IMAGE140
Wherein
Figure DEST_PATH_IMAGE141
Representing scattering clusters
Figure DEST_PATH_IMAGE142
Inner to the firstiScattering point
Figure DEST_PATH_IMAGE143
Is determined by the coordinate of (a) in the space,
Figure DEST_PATH_IMAGE144
representing effective radius of scattering cluster, for scattering cluster
Figure DEST_PATH_IMAGE145
All scattering points within the image are assigned a random velocity, denoted
Figure DEST_PATH_IMAGE146
。
3. The method of claim 2, wherein step 3 comprises:
the transmitting end/receiving end position is updated as follows:
Figure DEST_PATH_IMAGE147
Figure DEST_PATH_IMAGE148
wherein
Figure DEST_PATH_IMAGE149
Is referred to as the first of originpAn antenna
Figure DEST_PATH_IMAGE150
At the end ofqAn antenna
Figure DEST_PATH_IMAGE151
At the time of daytIn the position of (a) in the first,
Figure DEST_PATH_IMAGE152
is referred to as the first of originpAn antenna
Figure DEST_PATH_IMAGE153
At the end ofqAn antenna
Figure DEST_PATH_IMAGE154
At the moment of time
Figure DEST_PATH_IMAGE155
In the position of (a) in the first,
Figure DEST_PATH_IMAGE156
a velocity vector representing the originating/receiving end; the originating scattering cluster
Figure DEST_PATH_IMAGE157
End scattering cluster
Figure DEST_PATH_IMAGE158
The location update is as follows:
Figure DEST_PATH_IMAGE159
Figure DEST_PATH_IMAGE160
wherein
Figure DEST_PATH_IMAGE161
Refers to originating scattering clusters
Figure DEST_PATH_IMAGE162
Center/terminus scattering cluster
Figure DEST_PATH_IMAGE163
Centered at the momenttIn the position of (a) or (b),
Figure DEST_PATH_IMAGE164
is an originating scattering cluster
Figure DEST_PATH_IMAGE165
End scattering cluster
Figure DEST_PATH_IMAGE166
Centered at the moment
Figure DEST_PATH_IMAGE167
In the position of (a) in the first,
Figure DEST_PATH_IMAGE168
representing originating scattering clusters
Figure DEST_PATH_IMAGE169
End scattering cluster
Figure DEST_PATH_IMAGE170
At the moment of timetThe scattering cluster is in a random walk state, and the scattering cluster is initiated
Figure DEST_PATH_IMAGE171
End scattering cluster
Figure DEST_PATH_IMAGE172
The speed of (2) is updated as follows:
Figure DEST_PATH_IMAGE173
Figure DEST_PATH_IMAGE174
wherein
Figure DEST_PATH_IMAGE175
Originating scattering cluster
Figure DEST_PATH_IMAGE176
End scattering cluster
Figure DEST_PATH_IMAGE177
At the moment of time
Figure DEST_PATH_IMAGE178
The speed of the vehicle is measured by the speed sensor,
Figure DEST_PATH_IMAGE179
is a factor of the inertia of the object,
Figure DEST_PATH_IMAGE180
representing originating scattering clusters
Figure DEST_PATH_IMAGE181
End scattering cluster
Figure DEST_PATH_IMAGE182
At the time of day
Figure DEST_PATH_IMAGE183
Is expressed as follows:
Figure DEST_PATH_IMAGE184
Figure DEST_PATH_IMAGE185
wherein
Figure DEST_PATH_IMAGE186
Representing originating scattering clusters
Figure DEST_PATH_IMAGE187
End scattering cluster
Figure DEST_PATH_IMAGE188
At the moment of timetIn order to maintain track continuity, the phase is updated as follows:
Figure DEST_PATH_IMAGE189
Figure DEST_PATH_IMAGE190
wherein
Figure DEST_PATH_IMAGE191
Representing originating scattering clusters
Figure DEST_PATH_IMAGE192
End scattering cluster
Figure DEST_PATH_IMAGE193
At the moment of time
Figure DEST_PATH_IMAGE194
The horizontal phase of the instantaneous speed of the rotor,
Figure DEST_PATH_IMAGE195
which represents the maximum allowed phase shift and,
Figure DEST_PATH_IMAGE196
to be uniformly distributed in
Figure DEST_PATH_IMAGE197
Random variable in between;
the update beam pointing is specifically as follows: according to the updated positions of the transmitting end and the receiving end, calculating a departure angle and an arrival angle of the line-of-sight link, and updating the beam direction based on formulas (15) - (16);
according to the beam direction, scattering clusters outside the beam range are filtered, specifically as follows: for single hop links
Figure DEST_PATH_IMAGE198
And
Figure DEST_PATH_IMAGE199
Figure 489239DEST_PATH_IMAGE200
determining originating scattering clusters
Figure DEST_PATH_IMAGE201
Centered at the moment
Figure 468696DEST_PATH_IMAGE202
Position of
Figure DEST_PATH_IMAGE203
Whether the scattering cluster is positioned in the range of the transmitting end wave beam and the receiving end wave beam simultaneously, if not, removing the transmitting end scattering cluster
Figure 419072DEST_PATH_IMAGE204
(ii) a Judging receiving end scattering cluster
Figure DEST_PATH_IMAGE205
Centered at the moment
Figure 361752DEST_PATH_IMAGE206
In the position of
Figure DEST_PATH_IMAGE207
Whether the scattering cluster is located in the transmitting end wave beam range and the receiving end wave beam range at the same time, if not, removing the receiving end scattering cluster
Figure 972862DEST_PATH_IMAGE208
；
Calculating the survival probability of the scattering clusters in the intersection of the originating beam range and the terminating beam range according to the position and the speed of the scattering clusters at the current updating time, wherein the calculation formula is as follows:
Figure DEST_PATH_IMAGE209
Figure 841329DEST_PATH_IMAGE210
wherein
Figure DEST_PATH_IMAGE211
Representing originating scattering clusters
Figure 780467DEST_PATH_IMAGE212
End scattering cluster
Figure DEST_PATH_IMAGE213
The parameters of the birth or death over time,
Figure 776104DEST_PATH_IMAGE214
representing originating scattering clusters
Figure DEST_PATH_IMAGE215
End scattering cluster
Figure 490988DEST_PATH_IMAGE216
A birth-to-death parameter on the beam;
Figure DEST_PATH_IMAGE217
representing a reference distance, as a fixed constant，
Figure 281090DEST_PATH_IMAGE218
The reference angle is represented as a fixed constant;
Figure DEST_PATH_IMAGE219
indicating a link
Figure 848469DEST_PATH_IMAGE220
Pointing with the originating beam/link
Figure DEST_PATH_IMAGE221
And the beam direction of the receiving end.
CN202210778010.2A 2022-07-04 2022-07-04 Simulation method of beam channel based on random twin cluster Active CN114844584B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210778010.2A CN114844584B (en) 2022-07-04 2022-07-04 Simulation method of beam channel based on random twin cluster

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210778010.2A CN114844584B (en) 2022-07-04 2022-07-04 Simulation method of beam channel based on random twin cluster

Publications (2)

Publication Number Publication Date
CN114844584A CN114844584A (en) 2022-08-02
CN114844584B true CN114844584B (en) 2022-10-21

Family

ID=82574646

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210778010.2A Active CN114844584B (en) 2022-07-04 2022-07-04 Simulation method of beam channel based on random twin cluster

Country Status (1)

Country Link
CN (1) CN114844584B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014178287A1 (en) * 2013-05-03 2014-11-06 Nec Corporation A method used for a mobile communications system, a user equipment used in a mobile communications system and a simulator used for a mobile communications system
CN113949439A (en) * 2021-10-22 2022-01-18 东南大学 Geometric random channel modeling method for unmanned aerial vehicle air-to-air communication
CN114598408A (en) * 2022-03-10 2022-06-07 东南大学 6G universal channel modeling method suitable for full-band full scene

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014178287A1 (en) * 2013-05-03 2014-11-06 Nec Corporation A method used for a mobile communications system, a user equipment used in a mobile communications system and a simulator used for a mobile communications system
CN113949439A (en) * 2021-10-22 2022-01-18 东南大学 Geometric random channel modeling method for unmanned aerial vehicle air-to-air communication
CN114598408A (en) * 2022-03-10 2022-06-07 东南大学 6G universal channel modeling method suitable for full-band full scene

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
3D Non-Stationary Wideband Tunnel Channel Models for 5G High-Speed Train Wireless Communications;Yu Liu等;《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》;20200131;第21卷(第1期);全文 *
A Semi-Deterministic MIMO-based Beam Channel Model of Ray-Tracing and Propagation-Graph for mmWave Communications;Jiachi Zhang等;《WS24 IEEE ICC 2022 Workshop on Wireless Propagation Channels for 5G and B5G》;20220520;全文 *
基于3D 几何信道模型非对称波束信道特性分析;张嘉驰等;《北京邮电大学学报》;20220430;第45卷(第2期);全文 *
基于DMC几何随机簇信道建模与仿真性能评估;彭懿琳;《中国优秀硕士论文全文数据库（电子期刊）》;20150215;全文 *
面向真空管高速列车的无线通信系统关键技术研究;韩柏涛;《中国博士学位论文全文数据库（电子期刊）》;20220215;全文 *

Also Published As

Publication number Publication date
CN114844584A (en) 2022-08-02

Similar Documents

Publication Publication Date Title
CN107248877B (en) MassiveMIMO channel modeling method based on measured data
CN111245480B (en) Double-cluster geometric channel modeling method based on three-dimensional space
CN113162679A (en) DDPG algorithm-based IRS (inter-Range instrumentation System) auxiliary unmanned aerial vehicle communication joint optimization method
CN108718454B (en) Cooperative autonomous layout method for communication relay platforms of multiple unmanned aerial vehicles
CN105553584A (en) 3DMIMO channel modeling method
Jiang et al. Three-dimensional geometry-based stochastic channel modeling for intelligent reflecting surface-assisted UAV MIMO communications
Hayashi et al. A study on the variety and size of input data for radio propagation prediction using a deep neural network
CN111479239B (en) Sensor emission energy consumption optimization method of multi-antenna unmanned aerial vehicle data acquisition system
CN111147163A (en) Wireless communication link loss prediction method based on DNN neural network
CN113472419B (en) Safe transmission method and system based on space-based reconfigurable intelligent surface
CN110958069A (en) Three-dimensional non-stationary broadband double-cluster channel modeling method based on tunnel environment
CN110620627B (en) Non-stationary channel modeling method and device for vehicle-to-vehicle multi-antenna system
CN114268903B (en) Geographic information assisted unmanned aerial vehicle relay position deployment and power distribution method
CN112968743A (en) Time-varying de-cellular large-scale MIMO channel modeling method based on visible region division
CN113489560A (en) Geometric random modeling method for non-stationary air-ground MIMO channel of unmanned aerial vehicle
Huang et al. A non-stationary 6G V2V channel model with continuously arbitrary trajectory
Xu et al. A 3D non-stationary mmWave channel model for vacuum tube ultra-high-speed train channels
CN114980169A (en) Unmanned aerial vehicle auxiliary ground communication method based on combined optimization of track and phase
CN111479216A (en) Unmanned aerial vehicle cargo conveying method based on UWB positioning
CN114844584B (en) Simulation method of beam channel based on random twin cluster
CN113162712B (en) Multi-beam channel modeling method based on propagation diagram theory
CN111917498B (en) Simulation method for realizing millimeter wave channel space consistency
CN115765899A (en) Unmanned aerial vehicle communication beam domain channel simulation method and device, electronic equipment and medium
CN114584237A (en) Method for calculating capacity of non-stationary wireless channel in airspace in large-scale antenna array communication
Pattanayak A genetically trained neural network for prediction of path loss in outdoor microcell

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