CN114844584A - 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
CN114844584A
CN114844584A CN202210778010.2A CN202210778010A CN114844584A CN 114844584 A CN114844584 A CN 114844584A CN 202210778010 A CN202210778010 A CN 202210778010A CN 114844584 A CN114844584 A CN 114844584A
Authority
CN
China
Prior art keywords
scattering
cluster
clusters
link
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.)
Granted
Application number
CN202210778010.2A
Other languages
Chinese (zh)
Other versions
CN114844584B (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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Electromagnetism (AREA)
  • Mobile Radio Communication Systems (AREA)

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, 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 grow and extinguish 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 simultaneously do not consider the influence of beams on scattering clustering and extinction, 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 moment when the simulation starts, calculating beam response according to the initialization parameters of the 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 moment;
updating time, the positions of a transmitting and receiving end, beam pointing and scattering cluster positions, calculating to obtain a channel transmission function of each link at the 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 terminal and the terminating terminal,
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
represents the pitch width of the transmit/receive 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 an 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
Is determined by the coordinate of (a) in the space,
Figure 602916DEST_PATH_IMAGE040
representing effective radius of scattering cluster, for scattering clustermAll scattering points within the spectrum 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
Represents any horizontal angle and pitch angle of the transmitting end/receiving end,
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 originmA 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 represents the wavelength of the light emitted by the light source,
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
Distance ofAfter the separation, the water is separated from the water,
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 time of day
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
Originating/terminating scattering clusters
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 time of day
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 variables 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 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
Representing originating/receiving endScattering cluster of
Figure 228204DEST_PATH_IMAGE140
The parameters of the occurrence and extinction 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
representing a 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, if
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 from the poisson process for the mean
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 generation and the extinction of the scattering cluster, makes up for the fact that the random mobility of the scattering cluster and the influence of the wave beam on the generation and the extinction characteristics of the scattering cluster 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 function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
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 at the transmitting and receiving end is determined according to the width of the beam at 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 transmit/receive 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/receiving 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
。
Step S2, calculating the wave beam response according to the wave beam parameters of the receiving and transmitting end, constructing a line-of-sight link, a single-hop link and a double-hop link, and calculating the transmission function corresponding to each link by adopting the propagation diagram theory, thereby obtaining the 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 terminating endnA scattering cluster, and have
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,
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 the horizontal and pitch angles of approach at the closeout,
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
step S3, updating time, updating a position of a transceiver end, a beam direction and a scattering cluster position, filtering scattering clusters outside a beam range according to the beam direction, and calculating the survival probability of the scattering clusters 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 moment of timetThe 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 moment of time
Figure 887514DEST_PATH_IMAGE270
The speed of the vehicle is measured by the speed sensor,
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
to be 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 birth or death 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
a reference distance is indicated which is the distance of reference,is a constant with a constant value, and is,
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.
Step S4, determining the occurrence and the 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 deadline.
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 repeating all the steps corresponding to the formula (1-25) until the simulation deadline time to obtain the channel transmission function of the link after the update time.
In conclusion, the method of the embodiment of the invention makes up the problem that the random mobility of the scattering cluster and the influence of the 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 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.
All 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 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 (5)

1. A simulation method of beam channels based on random twin clusters is characterized by comprising 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 moment when the simulation starts, calculating beam response according to the initialization parameters of the 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 moment;
updating time, the positions of a transmitting and receiving end, beam pointing and scattering cluster positions, calculating to obtain a channel transmission function of each link at the 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.
2. The method of claim 1, wherein initializing transceiving end beam parameters including a beam direction and a beam width, wherein the beam direction is determined according to a transceiving end position, the number of transceiving end scattering point clusters is determined according to the width of the transceiving end beam, the position of the scattering point clusters is initialized by using a Matern hard kernel Poisson cluster process, and each scattering point cluster is given a random speed, comprises:
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 DEST_PATH_IMAGE001
And angle of departure pitch
Figure DEST_PATH_IMAGE002
The calculation formula of (a) is as follows:
Figure DEST_PATH_IMAGE003
Figure DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE005
Figure DEST_PATH_IMAGE006
wherein
Figure DEST_PATH_IMAGE007
And
Figure DEST_PATH_IMAGE008
respectively representing the array center coordinates of the originating terminal and the terminating terminal,
Figure DEST_PATH_IMAGE009
respectively representx、y、zAxial unit vector, sign
Figure DEST_PATH_IMAGE010
It is shown that the operation of the inner product of the vectors,
Figure DEST_PATH_IMAGE011
representing an inverse cosine function; similar horizontal angle of arrival for line-of-sight links
Figure DEST_PATH_IMAGE012
And angle of pitch of arrival
Figure DEST_PATH_IMAGE013
；
The beam pointing direction of the transceiving end is calculated as follows:
Figure DEST_PATH_IMAGE014
Figure DEST_PATH_IMAGE015
Figure DEST_PATH_IMAGE016
Figure DEST_PATH_IMAGE017
wherein
Figure DEST_PATH_IMAGE018
Indicating the transmit/receive beam horizontal pointing angle,
Figure DEST_PATH_IMAGE019
representing the transmit/receive beam elevation pointing angle;
Figure DEST_PATH_IMAGE020
indicating the horizontal angle of departure/arrival of the transceiving end line-of-sight link,
Figure DEST_PATH_IMAGE021
a pitch angle representing departure/arrival of the transceiving end line-of-sight link;
Figure DEST_PATH_IMAGE022
indicating the horizontal width of the originating/receiving beam,
Figure DEST_PATH_IMAGE023
represents the pitch width of the transmit/receive beam;
Figure DEST_PATH_IMAGE024
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 DEST_PATH_IMAGE025
Figure DEST_PATH_IMAGE026
wherein
Figure DEST_PATH_IMAGE027
Which represents the initial moment of time of day,
Figure DEST_PATH_IMAGE028
indicating the number of originating/terminating scattering clusters,
Figure DEST_PATH_IMAGE029
in order to be the density factor, the density of the sample,
Figure DEST_PATH_IMAGE030
representing an effective radius centered at the base station/terminal,
Figure DEST_PATH_IMAGE031
represents the minimum and maximum values of the transmit/receive beam in elevation, and has
Figure DEST_PATH_IMAGE032
；
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 DEST_PATH_IMAGE033
Figure DEST_PATH_IMAGE034
wherein
Figure DEST_PATH_IMAGE035
And
Figure DEST_PATH_IMAGE036
respectively representing the coordinates of the centers of two different scattering clusters,
Figure DEST_PATH_IMAGE037
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 DEST_PATH_IMAGE038
Figure DEST_PATH_IMAGE039
wherein
Figure DEST_PATH_IMAGE040
Representing scattering clustersmInner to the firstiScattering point
Figure DEST_PATH_IMAGE041
Is determined by the coordinate of (a) in the space,
Figure DEST_PATH_IMAGE042
representing effective radius of scattering cluster, for scattering clustermAll scattering points within the image are assigned a random velocity, denoted
Figure DEST_PATH_IMAGE043
。
3. The method of claim 1, wherein the determining the current time when the simulation starts, calculating a beam response according to initialization parameters of a transceiver end beam, 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 comprises:
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 DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE045
wherein
Figure DEST_PATH_IMAGE046
And
Figure DEST_PATH_IMAGE047
represents any horizontal angle and pitch angle of the transmitting end/receiving end,
Figure DEST_PATH_IMAGE048
maximum beam gain representing the transmit/receive beam;
the line-of-sight link refers to
Figure DEST_PATH_IMAGE049
In a direct link of
Figure DEST_PATH_IMAGE050
Refer to the originating endpAn antenna is arranged on the base plate, and the antenna is arranged on the base plate,
Figure DEST_PATH_IMAGE051
refer to the receiving endqAn antenna; the single-hop link is represented as
Figure DEST_PATH_IMAGE052
Indicating the proximity of the originating endmA plurality of scattering clusters, each scattering cluster having a plurality of scattering lines,
Figure DEST_PATH_IMAGE053
indicating the vicinity of the terminating endnA scattering cluster, and have
Figure DEST_PATH_IMAGE054
(ii) a The double-hop link is represented as
Figure DEST_PATH_IMAGE055
；
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_IMAGE056
Figure DEST_PATH_IMAGE057
wherein
Figure DEST_PATH_IMAGE058
Representing line-of-sight links
Figure DEST_PATH_IMAGE059
The distance between the first and second electrodes,
Figure DEST_PATH_IMAGE060
which represents the wavelength of the light emitted by the light source,
Figure DEST_PATH_IMAGE061
in order to be a time-delay variable,
Figure DEST_PATH_IMAGE062
is a variable of the time, and is,
Figure DEST_PATH_IMAGE063
in order to be the speed of light,
Figure DEST_PATH_IMAGE064
is a dirac function; the transmission function of the single-hop link is expressed as:
Figure DEST_PATH_IMAGE065
Figure DEST_PATH_IMAGE066
wherein
Figure DEST_PATH_IMAGE067
Indicating a link
Figure DEST_PATH_IMAGE068
The distance of (a) to (b),
Figure DEST_PATH_IMAGE069
representing scattering clustersnInner firstjScattering point
Figure DEST_PATH_IMAGE070
，
Figure DEST_PATH_IMAGE071
And
Figure DEST_PATH_IMAGE072
representing scattering clusters
Figure DEST_PATH_IMAGE073
The number of internal scattering points is greater than the total number of internal scattering points,
Figure DEST_PATH_IMAGE074
indicating a link
Figure DEST_PATH_IMAGE075
Gain, the expression is as follows:
Figure DEST_PATH_IMAGE076
Figure DEST_PATH_IMAGE077
Figure DEST_PATH_IMAGE078
Figure DEST_PATH_IMAGE079
wherein
Figure DEST_PATH_IMAGE080
Representing scattering points
Figure DEST_PATH_IMAGE081
Relative to the horizontal and pitch angles of departure from the origin,
Figure DEST_PATH_IMAGE082
representing scattering points
Figure DEST_PATH_IMAGE083
Relative to the horizontal and pitch angles of approach at the closeout,
Figure DEST_PATH_IMAGE084
respectively represent
Figure DEST_PATH_IMAGE085
And
Figure DEST_PATH_IMAGE086
the calculation formula of the average transmission delay is as follows:
Figure DEST_PATH_IMAGE087
Figure DEST_PATH_IMAGE088
Figure DEST_PATH_IMAGE089
Figure DEST_PATH_IMAGE090
Figure DEST_PATH_IMAGE091
Figure DEST_PATH_IMAGE092
Figure DEST_PATH_IMAGE093
Figure DEST_PATH_IMAGE094
Figure DEST_PATH_IMAGE095
the transfer function of the dual-hop link is expressed as:
Figure DEST_PATH_IMAGE096
Figure DEST_PATH_IMAGE097
wherein
Figure DEST_PATH_IMAGE098
Indicating a link
Figure DEST_PATH_IMAGE099
The distance of (a) to (b),
Figure DEST_PATH_IMAGE100
indicating a link
Figure DEST_PATH_IMAGE101
The expression is:
Figure DEST_PATH_IMAGE102
Figure DEST_PATH_IMAGE103
wherein
Figure DEST_PATH_IMAGE104
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_IMAGE105
Figure DEST_PATH_IMAGE106
。
4. the method of claim 1, wherein the updating time, the positions of the transceiving ends, 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 probability of survival of the scattering clusters within the beam range is calculated according to the channel transfer function of each link at each time, and the method comprises:
the receiving and transmitting end position is updated as follows:
Figure DEST_PATH_IMAGE107
Figure DEST_PATH_IMAGE108
wherein
Figure DEST_PATH_IMAGE109
Refer to
Figure DEST_PATH_IMAGE110
At the moment of timetIn the position of (a) in the first,
Figure DEST_PATH_IMAGE111
refer to
Figure DEST_PATH_IMAGE112
At the moment of time
Figure DEST_PATH_IMAGE113
In the position of (a) in the first,
Figure DEST_PATH_IMAGE114
representing the velocity vectors of the originating terminal and the terminating terminal; the scattering cluster
Figure DEST_PATH_IMAGE115
The location update is as follows:
Figure DEST_PATH_IMAGE116
Figure DEST_PATH_IMAGE117
wherein
Figure DEST_PATH_IMAGE118
Refers to originating/terminating scattering clusters
Figure DEST_PATH_IMAGE119
Centered at the momenttIn the position of (a) or (b),
Figure DEST_PATH_IMAGE120
the scattering cluster center of the transmitting end/receiving end is in the time
Figure DEST_PATH_IMAGE121
In the position of (a) in the first,
Figure DEST_PATH_IMAGE122
representing originating/receiving scattering clusters
Figure 775753DEST_PATH_IMAGE119
At the moment of timetThe scattering cluster is in a random walk state, the scattering cluster
Figure 240364DEST_PATH_IMAGE119
The speed of (2) is updated as follows:
Figure DEST_PATH_IMAGE123
Figure DEST_PATH_IMAGE124
wherein
Figure DEST_PATH_IMAGE125
Transmitting/receiving end scattering cluster
Figure 659581DEST_PATH_IMAGE119
At the moment of time
Figure 585949DEST_PATH_IMAGE121
The instantaneous speed of the vehicle (c) is,
Figure DEST_PATH_IMAGE126
is a factor of the inertia, and is,
Figure DEST_PATH_IMAGE127
representing originating/receiving scattering clusters
Figure 275688DEST_PATH_IMAGE119
At the moment of time
Figure 629309DEST_PATH_IMAGE121
The instantaneous velocity direction vector of (a), expressed as follows:
Figure DEST_PATH_IMAGE128
Figure DEST_PATH_IMAGE129
wherein
Figure DEST_PATH_IMAGE130
Representing originating/receiving scattering clusters
Figure 538752DEST_PATH_IMAGE119
At the moment of timetIn order to maintain track continuity, the phase is updated as follows:
Figure DEST_PATH_IMAGE131
Figure DEST_PATH_IMAGE132
wherein
Figure DEST_PATH_IMAGE133
Which represents the maximum allowed phase shift and,
Figure DEST_PATH_IMAGE134
to be uniformly distributed in
Figure DEST_PATH_IMAGE135
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 DEST_PATH_IMAGE136
And
Figure DEST_PATH_IMAGE137
determining the center of the scattering cluster
Figure DEST_PATH_IMAGE138
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 DEST_PATH_IMAGE139
Figure DEST_PATH_IMAGE140
wherein
Figure DEST_PATH_IMAGE141
Scattering cluster for representing originating/receiving end
Figure DEST_PATH_IMAGE142
The parameters of the birth or death over time,
Figure DEST_PATH_IMAGE143
scattering cluster for representing originating/receiving end
Figure DEST_PATH_IMAGE144
A birth-to-death parameter on the beam;
Figure DEST_PATH_IMAGE145
which represents the reference distance, is a fixed constant,
Figure DEST_PATH_IMAGE146
the reference angle is represented as a fixed constant;
Figure DEST_PATH_IMAGE147
indicating a link
Figure DEST_PATH_IMAGE148
The angle between the beam direction and the transmitting end/receiving end.
5. The method of claim 1, wherein the occurrence and extinction of the scattering clusters after the update time is determined according to the scattering cluster occurrence probability, the number of surviving scattering clusters is counted, the number of scattering clusters at the moment is determined based on the poisson process, and if the surviving number is lower than the threshold, 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 DEST_PATH_IMAGE149
If, if
Figure DEST_PATH_IMAGE150
Then the scattering cluster of the transmitting/receiving end
Figure DEST_PATH_IMAGE151
At the moment of time
Figure DEST_PATH_IMAGE152
Survival; if it is
Figure DEST_PATH_IMAGE153
Then the scattering cluster of the transmitting/receiving end
Figure 406826DEST_PATH_IMAGE151
At the time
Figure 606864DEST_PATH_IMAGE152
Eliminating;
to be provided with
Figure DEST_PATH_IMAGE154
Generating a random variable from the poisson process for the mean
Figure DEST_PATH_IMAGE155
And for random variables
Figure 944435DEST_PATH_IMAGE155
Rounding is carried out if the number of the scattering clusters surviving at the transmitting end/the receiving end is less than
Figure 495502DEST_PATH_IMAGE155
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.
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 true CN114844584A (en) 2022-08-02
CN114844584B 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
JIACHI ZHANG等: "A Semi-Deterministic MIMO-based Beam Channel Model of Ray-Tracing and Propagation-Graph for mmWave Communications", 《WS24 IEEE ICC 2022 WORKSHOP ON WIRELESS PROPAGATION CHANNELS FOR 5G AND B5G》 *
YU LIU等: "3D Non-Stationary Wideband Tunnel Channel Models for 5G High-Speed Train Wireless Communications", 《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》 *
张嘉驰等: "基于3D 几何信道模型非对称波束信道特性分析", 《北京邮电大学学报》 *
彭懿琳: "基于DMC几何随机簇信道建模与仿真性能评估", 《中国优秀硕士论文全文数据库（电子期刊）》 *
韩柏涛: "面向真空管高速列车的无线通信***关键技术研究", 《中国博士学位论文全文数据库（电子期刊）》 *

Also Published As

Publication number Publication date
CN114844584B (en) 2022-10-21

Similar Documents

Publication Publication Date Title
CN113162679B (en) DDPG algorithm-based IRS (intelligent resilient software) assisted unmanned aerial vehicle communication joint optimization method
CN107248877B (en) MassiveMIMO channel modeling method based on measured data
Jiang et al. Three-dimensional geometry-based stochastic channel modeling for intelligent reflecting surface-assisted UAV MIMO communications
CN107579789B (en) Large-scale unmanned aerial vehicle relay network channel simulation device and GPU real-time simulation method
CN113472419B (en) Safe transmission method and system based on space-based reconfigurable intelligent surface
CN108718454B (en) Cooperative autonomous layout method for communication relay platforms of multiple unmanned aerial vehicles
CN111479239B (en) Sensor emission energy consumption optimization method of multi-antenna unmanned aerial vehicle data acquisition system
Hayashi et al. A study on the variety and size of input data for radio propagation prediction using a deep neural network
CN105553584A (en) 3DMIMO channel modeling method
CN112968743B (en) Time-varying de-cellular large-scale MIMO channel modeling method based on visible region division
Dhekne et al. Extending cell tower coverage through drones
CN110958069A (en) Three-dimensional non-stationary broadband double-cluster channel modeling method based on tunnel environment
Huang et al. A non-stationary 6G V2V channel model with continuously arbitrary trajectory
Lee et al. Path loss exponent prediction for outdoor millimeter wave channels through deep learning
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
Wang et al. Joint power and QoE optimization scheme for multi-UAV assisted offloading in mobile computing
CN113489560A (en) Geometric random modeling method for non-stationary air-ground MIMO channel of unmanned aerial vehicle
CN113784314B (en) Unmanned aerial vehicle data and energy transmission method assisted by intelligent reflection surface
CN114980169A (en) Unmanned aerial vehicle auxiliary ground communication method based on combined optimization of track and phase
Xu et al. A 3D non-stationary mmWave channel model for vacuum tube ultra-high-speed train channels
CN115765899A (en) Unmanned aerial vehicle communication beam domain channel simulation method and device, electronic equipment and medium
CN114844584B (en) Simulation method of beam channel based on random twin cluster
Eldowek et al. A survey of 5G millimeter wave, massive multiple‐input multiple‐output, and vehicle‐to‐vehicle channel measurements and models
CN113162712B (en) Multi-beam channel modeling method based on propagation diagram theory

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