CN110958041A - Dynamic user hybrid precoding method for millimeter wave system - Google Patents

Dynamic user hybrid precoding method for millimeter wave system Download PDF

Info

Publication number
CN110958041A
CN110958041A CN201910965145.8A CN201910965145A CN110958041A CN 110958041 A CN110958041 A CN 110958041A CN 201910965145 A CN201910965145 A CN 201910965145A CN 110958041 A CN110958041 A CN 110958041A
Authority
CN
China
Prior art keywords
user
storage unit
path
matrix
mth
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
CN201910965145.8A
Other languages
Chinese (zh)
Other versions
CN110958041B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201910965145.8A priority Critical patent/CN110958041B/en
Publication of CN110958041A publication Critical patent/CN110958041A/en
Application granted granted Critical
Publication of CN110958041B publication Critical patent/CN110958041B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters

Landscapes

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

Abstract

The invention relates to a dynamic user hybrid precoding method of a millimeter wave system, which comprises the following steps: acquiring channel state information of each user, and storing the channel state information of adjacent time slots of each user in a storage unit; obtaining the optimal propagation path of each user according to the stored channel state information; judging whether the optimal propagation path is a direct path or not to obtain a judgment result of each user; obtaining the state of each user according to the judgment result; designing a simulation pre-coding matrix of the user according to the state of the user; designing a digital pre-coding matrix of a user according to the analog pre-coding matrix; and carrying out data transmission according to the analog precoding matrix and the digital precoding matrix. The method of the invention can obtain the state of each user only by comparing and judging the channel state information of the adjacent time slots of the users, thereby improving the efficiency of user state detection and reducing the complexity of detection, and provides a self-adaptive mixed pre-coding method according to the real-time state of the detected user.

Description

Dynamic user hybrid precoding method for millimeter wave system
Technical Field
The invention belongs to the technical field of wireless communication, and particularly relates to a dynamic user hybrid precoding method for a millimeter wave system.
Background
As one of the key technologies of the fifth generation mobile communication system, the millimeter wave can effectively solve the problem of global frequency band shortage and greatly improve the transmission speed and the spectrum efficiency of the system. However, the sparse channels and narrow beams in mmwave systems make their performance extremely sensitive to the user state in the system. For example, when the millimeter wave beam is blocked by a human body, about 35dB of beam gain will be lost. In addition, the millimeter wave beam gain also suffers from different degrees of attenuation when the user is in motion.
Aiming at the precoding problem in a millimeter wave partial connection antenna array system, a plurality of scholars have conducted intensive research at present, for example, a simulation precoding matrix is designed by utilizing a simulation precoding matrix block diagonalization structure and a semi-positive definite relaxation algorithm so as to reduce the Euclidean distance between a hybrid precoding matrix and a full digital precoding matrix; obtaining an optimal transmission direction through finite search in a beam space, and taking an array steering vector corresponding to the direction as a corresponding column of a simulation pre-coding matrix to simplify the design complexity of the hybrid pre-coding matrix; the calculation complexity in the optimization process is reduced by converting the optimization problem of the traditional analog precoding matrix into sub-problems of respectively optimizing the amplitude and the phase.
However, the conventional hybrid precoding algorithm only considers optimizing the sum-rate performance of the static users and simplifying the computational complexity in the design process of the static user hybrid precoding matrix. When dynamic users exist in the system, the sum rate performance of the traditional hybrid precoding algorithm is extremely easily influenced by the user states in the system. The above conventional algorithms are therefore not suitable for dynamic user scenarios or it is difficult to maintain their expected performance.
Disclosure of Invention
In order to solve the above problems in the prior art, the present invention provides a dynamic user hybrid precoding method for millimeter wave system. The technical problem to be solved by the invention is realized by the following technical scheme, the invention provides a dynamic user hybrid precoding method of a millimeter wave system, which is used for a millimeter wave part connection antenna array system and comprises the following steps:
acquiring channel state information of each user, and storing the channel state information of adjacent time slots of each user in a storage unit;
obtaining the optimal propagation path of each user according to the stored channel state information;
judging whether the optimal propagation path is a direct path or not to obtain a judgment result of each user;
obtaining the state of each user according to the judgment result;
designing a simulation pre-coding matrix of the user according to the state of the user;
designing a digital pre-coding matrix of a user according to the analog pre-coding matrix;
and carrying out data transmission according to the analog precoding matrix and the digital precoding matrix.
In an embodiment of the present invention, acquiring channel state information of each user and storing the channel state information of each user's adjacent time slot in a storage unit includes:
the base station obtains the channel state information H scanned by each user in the beam training time slotm，
Figure BDA0002230244020000021
Wherein HmA channel state information matrix representing the mth user, M representing the number of users,
Figure BDA0002230244020000022
a channel state matrix representing the direct path of the mth user,
Figure BDA0002230244020000023
a channel state matrix representing the ith non-direct path of the mth user, K represents the Rice factor, PmIndicates the number of propagation paths of the mth user,
Figure BDA0002230244020000024
a complex path gain representing an ith non-direct path of an mth user;
and respectively storing the channel state information of the t-1 time slot and the t time slot of each user in a first storage unit and a second storage unit, wherein t is more than or equal to 1.
In an embodiment of the present invention, obtaining an optimal propagation path of each user according to the stored channel state information includes:
obtaining an optimal propagation path of each user in the first storage unit and the second storage unit according to the following formula,
Figure BDA0002230244020000025
wherein g represents a memory cell,
Figure BDA0002230244020000026
indicating the departure angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure BDA0002230244020000027
indicating the arrival angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure BDA0002230244020000028
represents the complex path gain of the xth propagation path of the user in the g storage unit, | | | represents taking the absolute value, max represents taking the maximum value,
Figure BDA0002230244020000029
the number of propagation paths of the mth user in the g-th storage unit is shown.
In an embodiment of the present invention, determining whether the optimal propagation path is a direct path to obtain a determination result of each user includes:
if the optimal propagation path satisfies
Figure BDA0002230244020000031
The optimal propagation path of the user in the g-th storage unit is a direct propagation path, and if the optimal propagation path meets the requirement
Figure BDA0002230244020000032
The optimal propagation path of the user at the g-th storage unit is a non-direct path,
wherein epsilon represents the threshold values of the direct path and the indirect path, and epsilon is more than or equal to 0 and less than or equal to muAS，μASRepresents the average of the angular spread of the propagating beam.
In an embodiment of the present invention, obtaining the state of each user according to the determination result includes:
if the optimal propagation paths of the user in the first storage unit and the second storage unit are direct paths, executing the following steps:
if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit satisfies the requirement
Figure BDA0002230244020000033
The user is in a static state, and if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit meets the requirement
Figure BDA0002230244020000034
The user is in a state of motion and,
wherein, delta represents the judgment threshold value of the motion and the static state of the user,
Figure BDA0002230244020000035
μASan average value representing the angular spread of the propagation beam, BW representing the half-power beamwidth of the base station propagation beam;
if the optimal propagation path of the user in the first storage unit is a direct path and the optimal propagation path in the second storage unit is a non-direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtainTo
Figure BDA0002230244020000036
The number of the path pairs is one,
Figure BDA0002230244020000037
if it is
Figure BDA0002230244020000038
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA0002230244020000039
The user is in a blocking + movement state, otherwise the user is in a blocking state;
if the optimal propagation path of the user in the first storage unit is a non-direct path and the optimal propagation path in the second storage unit is a direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure BDA00022302440200000310
The number of the path pairs is one,
Figure BDA00022302440200000311
if it is
Figure BDA0002230244020000041
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA0002230244020000042
The user is in a motion state, otherwise, the user is in a static state;
if the optimal propagation paths of the user in the first storage unit and the second storage unit are both non-direct paths, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure BDA0002230244020000043
The number of the path pairs is one,
Figure BDA0002230244020000044
if it is
Figure BDA0002230244020000045
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA0002230244020000046
The user is in the block + move state otherwise the user is in the block state.
In an embodiment of the present invention, designing an analog precoding matrix of a user according to a state of the user includes:
prioritizing a status of the user, wherein,
stationary > occlusion > motion > occlusion + motion;
designing a simulation pre-coding matrix of the first user according to the priority sequence;
and designing analog precoding matrixes of other users to obtain the analog precoding matrixes of all the users.
In an embodiment of the present invention, designing an analog precoding matrix of the first user according to the priority ranking includes:
i. let m equal to 1, wmFor all-zero vectors, a precoding matrix F is simulatedRFJudging the state of the first user for an all-zero matrix, wherein the simulation precoding matrix FRFIn order to realize the purpose,
Figure BDA0002230244020000047
wherein, simulatingPrecoding matrix FRFHas the dimension of
Figure BDA0002230244020000048
NtIndicating the number of transmit antennas of the base station,
Figure BDA0002230244020000049
indicating the number of transmit sub-arrays of the base station, flAn analog precoding vector representing the ith transmit sub-array of the base station,
Figure BDA00022302440200000410
wmrepresenting the mth user dimension as NrX 1 received combined vector, NrIndicating the number of receive antennas per user;
if the status of the first user is static or blocked, then
Figure BDA0002230244020000051
f1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝1，
Figure BDA0002230244020000052
Wherein phi issRepresenting the angle of departure, theta, of the user's propagating beamsThe arrival angle of the user's propagation beam is shown, the propagation beam of the user is the departure angle and the arrival angle corresponding to the path gain maximum path, wm(1: N) represents the 1 st to Nth elements of the m-th user reception combining vector, N represents the number of transmission sub-arrays used in the base station,
Figure BDA0002230244020000053
aN(x) Indicating the array steering vector, N indicating the number of linear antennas per transmit sub-array of the base station,
Figure BDA0002230244020000054
λ represents the carrier wavelength, d represents the spacing between antenna elements,
Figure BDA0002230244020000055
if the state of the first user is a moving or moving + blocking state, passing through the adjacent N in the base stationcCooperation of transmission sub-arrays, and N of userscThe receiving sub-arrays cooperate to carry out beam combination, and the angle difference of the transmission beams of the adjacent transmitting sub-arrays and the adjacent receiving sub-arrays is
Figure BDA0002230244020000056
Then it is determined that,
Figure BDA0002230244020000057
if the tangential direction to the first user's direction of motion is counterclockwise, then,
Figure BDA0002230244020000058
Figure BDA0002230244020000059
Figure BDA00022302440200000510
Figure BDA0002230244020000061
Figure BDA0002230244020000062
n＝Nc，
if the tangential direction to the first user's direction of motion is clockwise, then,
Figure BDA0002230244020000063
Figure BDA0002230244020000064
Figure BDA0002230244020000065
Figure BDA0002230244020000066
Figure BDA0002230244020000067
n＝Nc；
wherein, wm(N +1:2N) denotes that the mth user receives the N +1 th to 2N-th elements of the merged vector, wm[(Nc-1)N+1:NcN](N) th user receiving combined vectorc-1) N +1 to NthcN elements;
updating the analog precoding matrix FRFA value of (d);
in an embodiment of the present invention, designing the analog precoding matrices of other users to obtain the analog precoding matrices of all users includes:
i. selecting a propagation beam for the mth user, determining a departure angle and an arrival angle of the propagation beam for the mth user,
Figure BDA0002230244020000068
Figure BDA0002230244020000069
wherein | | | purple hair2The norm of 2 is taken, and min is the minimum value;
judging the status of the mth user, if the status of the mth user is static or blocking, then,
Figure BDA0002230244020000071
fn+1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝n+1，
Figure BDA0002230244020000072
if the mth user state is a moving or moving + blocking state, passing through the adjacent N in the base stationcCooperation of the transmitting sub-arrays, and N of userscThe receiving sub-arrays cooperate to carry out beam combination, and the angle difference of the transmission beams of the adjacent transmitting sub-arrays and the adjacent receiving sub-arrays is
Figure BDA0002230244020000073
Then it is determined that,
Figure BDA0002230244020000074
if the tangential direction of the mth user's motion direction is counterclockwise, then,
Figure BDA0002230244020000075
Figure BDA0002230244020000076
Figure BDA0002230244020000077
Figure BDA0002230244020000078
Figure BDA0002230244020000079
n＝n+Nc，
if the tangential direction of the mth user's motion direction is clockwise, then,
Figure BDA00022302440200000710
Figure BDA00022302440200000711
Figure BDA0002230244020000081
Figure BDA0002230244020000082
Figure BDA0002230244020000083
n＝n+Nc；
updating the analog precoding matrix FRFIs given as M +1, repeating steps i-iv, and iterating M-1 times to obtain the analog precoding matrix F of the userRF。
In one embodiment of the present invention, designing a digital precoding matrix of a user according to the analog precoding matrix comprises:
constructing a combined interference matrix for the mth user
Figure BDA0002230244020000084
Figure BDA0002230244020000085
Figure BDA0002230244020000086
Wherein the content of the first and second substances,
Figure BDA0002230244020000087
the dimension of (a) is M-1 xM,
Figure BDA0002230244020000088
a valid channel vector representing the mth user;
combined interference matrix to the mth user
Figure BDA0002230244020000089
The singular value decomposition is carried out, and the singular value decomposition,
Figure BDA00022302440200000810
wherein, UmRepresenting combined interference matrices
Figure BDA00022302440200000811
Left singular matrix of, VmRepresenting combined interference matrices
Figure BDA00022302440200000812
Right singular matrix ofmRepresenting a singular value diagonal matrix;
according to the right singular matrix VmObtaining the digital precoding vector F of the mth userBB(:,m)，
FBB(:,m)＝Vm(:,M)，
Wherein, FBB(: m) denotes a digital precoding matrix FBBRow m of (5), Vm(M) represents the right singular matrix VmThe M-th column of (1);
repeating the steps by making M equal to M +1, and iterating for M-1 times to obtain the analog precoding matrix F of the userBB。
In an embodiment of the present invention, performing data transmission according to the analog precoding matrix and the digital precoding matrix includes:
the base station transmits data according to the channel state information, the analog pre-coding matrix, the digital pre-coding matrix and the receiving merging vector of the user, and the receiving signal y of the mth usermIn order to realize the purpose,
Figure BDA0002230244020000091
wherein n ismRepresenting the received noise vector for the mth user and x representing the transmitted signal vector.
Compared with the prior art, the invention has the beneficial effects that:
1. in the millimeter wave system dynamic user hybrid precoding method, the state of each user can be obtained only by comparing and judging the channel state information of the adjacent time slots of the users without a user positioning system or a gyroscope, so that the user state detection efficiency is improved, and the detection complexity is reduced;
2. the millimeter wave system dynamic user hybrid precoding method of the invention utilizes the cooperation of adjacent antenna sub-arrays and enables the interval between the departure angle and the arrival angle of the adjacent sub-arrays to transmit the wave beams to be
Figure BDA0002230244020000092
Under the condition of the system, stable and wide beams are synthesized, so that the stability of the motion user in the millimeter wave system is improved, the time for the motion user to move out of the beams is prolonged, and the training period of the system beams is effectively prolonged.
The foregoing description is only an overview of the technical solutions of the present invention, and in order to make the technical means of the present invention more clearly understood, the present invention may be implemented in accordance with the content of the description, and in order to make the above and other objects, features, and advantages of the present invention more clearly understood, the following preferred embodiments are described in detail with reference to the accompanying drawings.
Drawings
FIG. 1 is a system model diagram according to an embodiment of the present invention;
fig. 2 is a flowchart of a dynamic user hybrid precoding method for a millimeter wave system according to an embodiment of the present invention;
FIG. 3 is a diagram comparing beam patterns provided by embodiments of the present invention;
FIG. 4 is a diagram comparing beam patterns provided by embodiments of the present invention;
FIG. 5 is a simulation graph of user movement time and rate over time provided by an embodiment of the present invention;
FIG. 6 is a simulation graph of time and rate of user movement provided by an embodiment of the present invention;
FIG. 7 is a simulation graph of time and rate of user movement provided by an embodiment of the present invention;
FIG. 8 is a simulation graph of time and rate of movement of a user according to another embodiment of the present invention;
FIG. 9 is a simulation diagram of random time and rate variation with time for a user state according to an embodiment of the present invention;
FIG. 10 is a diagram of another simulation of random user state and rate over time provided by an embodiment of the present invention;
FIG. 11 is a graph of a simulation of random time and rate over time for a user state provided by an embodiment of the present invention;
fig. 12 is a simulation diagram of random time and rate variation with time for a user state according to another embodiment of the present invention.
Detailed Description
In order to further explain the technical means and effects of the present invention adopted to achieve the predetermined object, the following describes in detail a dynamic user hybrid precoding method for millimeter wave system according to the present invention with reference to the accompanying drawings and the detailed description.
The foregoing and other technical matters, features and effects of the present invention will be apparent from the following detailed description of the embodiments, which is to be read in connection with the accompanying drawings. The technical means and effects of the present invention adopted to achieve the predetermined purpose can be more deeply and specifically understood through the description of the specific embodiments, however, the attached drawings are provided for reference and description only and are not used for limiting the technical scheme of the present invention.
Example one
Referring to fig. 1, fig. 1 is a system model diagram according to an embodiment of the present invention, and as shown in the diagram, a millimeter wave communication system used in this embodiment is formed by connecting a millimeter wave section with an antenna array base station, a transmission channel, and a user. Base station has
Figure BDA0002230244020000101
A plurality of transmitting sub-arrays, each transmitting sub-array having N antennas, the base station having a common structure
Figure BDA0002230244020000102
A transmitting antenna, M users in the system, each user having
Figure BDA0002230244020000103
A receiving sub-array, each transmitting sub-array has N antennas, each user has
Figure BDA0002230244020000104
A plurality of receiving antennas, wherein each antenna is independent and the adjacent antennas are spaced from each other by a distance of
Figure BDA0002230244020000105
Each user state is random and includes four states of rest, motion, occlusion, and occlusion + motion.
Referring to fig. 2, fig. 2 is a flowchart of a dynamic user hybrid precoding method for a millimeter wave system according to an embodiment of the present invention, where as shown in the figure, the method according to the embodiment includes:
s1: acquiring channel state information of each user, and storing the channel state information of adjacent time slots of each user in a storage unit;
specifically, the method comprises the following steps:
s11: the base station obtains the channel state information H scanned by each user in the beam training time slotm，
Figure BDA0002230244020000106
Wherein HmA channel state information matrix representing the mth user, M representing the number of users,
Figure BDA0002230244020000107
a channel state matrix representing the direct path of the mth user,
Figure BDA0002230244020000108
a channel state matrix representing the ith non-direct path of the mth user, K represents the Rice factor, PmIndicates the number of propagation paths of the mth user,
Figure BDA0002230244020000109
a complex path gain representing an ith non-direct path of an mth user;
in this embodiment, the channel state information includes departure angle, arrival angle and path gain between each user and a uniform linear antenna array base station,
Figure BDA0002230244020000111
and
Figure BDA0002230244020000112
the expression of (a) is as follows,
Figure BDA0002230244020000113
Figure BDA0002230244020000114
Figure BDA0002230244020000115
Figure BDA0002230244020000116
wherein the content of the first and second substances,
Figure BDA0002230244020000117
λ represents the carrier wavelength, d represents the spacing between antenna elements,
Figure BDA0002230244020000118
Figure BDA0002230244020000119
representing the angle of arrival of the ith propagation path of the mth user,
Figure BDA00022302440200001110
denotes the departure angle of the ith propagation path of the mth user, i is 0,1,2, …, PmAnd-1, when i is equal to 0, the path is a direct path, and when i is equal to 0, the path is a non-direct path.
S12: and respectively storing the channel state information of the t-1 time slot and the t time slot of each user in a first storage unit and a second storage unit, wherein t is more than or equal to 1.
In this embodiment, the first storage unit and the second storage unit only store the channel state information of the adjacent time slot of each user, and when the channel state information of the next time slot needs to be stored, the channel state information of the previous time slot is overwritten or deleted.
S2: obtaining the optimal propagation path of each user according to the stored channel state information;
specifically, the method comprises the following steps:
obtaining an optimal propagation path of each user in the first storage unit and the second storage unit according to the following formula,
Figure BDA00022302440200001111
wherein g represents a memory cell,
Figure BDA00022302440200001112
indicating the departure angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure BDA00022302440200001113
indicating the arrival angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure BDA00022302440200001114
represents the complex path gain of the xth propagation path of the user in the g storage unit, | | | represents taking the absolute value, max represents taking the maximum value,
Figure BDA00022302440200001115
the number of propagation paths of the mth user in the g-th storage unit is shown.
S3: judging whether the optimal propagation path is a direct path or not to obtain a judgment result of each user;
specifically, the method comprises the following steps:
if the optimal propagation path satisfies
Figure BDA0002230244020000121
The optimal propagation path of the user in the g storage unit is a direct path;
if the optimal propagation path satisfies
Figure BDA0002230244020000122
The optimal propagation path of the user in the g storage unit is a non-direct path;
wherein epsilon represents the threshold values of the direct path and the indirect path, and epsilon is more than or equal to 0 and less than or equal to muAS，μASThe average value of the angle spread of the propagation beam is represented, and in this embodiment, the threshold value is set according to the angle spread value of the system and the system requirement. As shown in the above equation, if the absolute value of the sum of the departure angle and the reception angle of the optimal propagation path in the g-th storage unit is smaller than ∈, the optimal propagation path is a direct path, and otherwise, the optimal propagation path is a non-direct path.
S4: obtaining the state of each user according to the judgment result;
in this embodiment, the states of the user include a still state, a moving state, a blocking state and a blocking + moving state, where the still state represents that the user is not blocked and does not move, the moving state represents that the user is in a moving state, the blocking state represents that the millimeter wave beam is blocked by an obstacle, and the blocking + moving state represents that the user has both moving and blocking.
Specifically, the method comprises the following steps:
s41: if the optimal propagation paths of the user in the first storage unit and the second storage unit are direct paths, executing the following steps:
if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit satisfies the requirement
Figure BDA0002230244020000123
The user is in a stationary state;
if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit satisfies the requirement
Figure BDA0002230244020000124
The user is in motion;
wherein, δ represents the judgment threshold of the user's motion and static state, and the value thereof is determined according to the system requirement, in the embodiment,
Figure BDA0002230244020000125
μASan average value representing the angular spread of the propagation beam, BW representing the half-power beamwidth of the base station propagation beam;
s42: if the optimal propagation path of the user in the first storage unit is a direct path and the optimal propagation path in the second storage unit is a non-direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure BDA0002230244020000126
The number of the path pairs is one,
Figure BDA0002230244020000131
if it is
Figure BDA0002230244020000132
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA0002230244020000133
The user is in a blocking + movement state, otherwise the user is in a blocking state;
s43: if the optimal propagation path of the user in the first storage unit is a non-direct path and the optimal propagation path in the second storage unit is a direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure BDA0002230244020000134
The number of the path pairs is one,
Figure BDA0002230244020000135
if it is
Figure BDA0002230244020000136
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA0002230244020000137
The user is in a motion state, otherwise, the user is in a static state;
s44: if the optimal propagation paths of the user in the first storage unit and the second storage unit are both non-direct paths, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure BDA0002230244020000138
The number of the path pairs is one,
Figure BDA0002230244020000139
if it is
Figure BDA00022302440200001310
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure BDA00022302440200001311
The user is in the block + move state otherwise the user is in the block state.
In the millimeter wave system dynamic user hybrid precoding method of the embodiment, the state of each user can be obtained only by comparing and judging the channel state information of the adjacent time slots of the users without a user positioning system or a gyroscope, so that the user state detection efficiency is improved, and the detection complexity is reduced.
S5: designing a simulation pre-coding matrix of the user according to the state of the user;
specifically, the method comprises the following steps:
s51: prioritizing a status of the user, wherein,
stationary > occlusion > motion > occlusion + motion;
in this embodiment, the processing priority of the user in the static state is the highest, the data transfer processing is performed first, and the processing priority of the user in the blocking + moving state is the lowest, and finally, the data transfer processing is performed, so that if a plurality of users exist in the same state and the processing priorities among the plurality of users are not in sequence, any one of the users can be selected to start the data transfer processing.
S52: designing a simulation pre-coding matrix of the first user according to the priority sequence
Specifically, the method comprises the following steps:
i. let m equal to 1, wmFor all-zero vectors, a precoding matrix F is simulatedRFFor all-zero matrix, the state of the first user is enteredRow judgment, wherein the analog precoding matrix FRFIn order to realize the purpose,
Figure BDA0002230244020000141
wherein, the analog precoding matrix FRFHas the dimension of
Figure BDA0002230244020000142
NtIndicating the number of transmit antennas of the base station,
Figure BDA0002230244020000143
representing the number of transmit subarrays of a base station, in this embodiment, the millimeter wave base station antenna and the user antenna are partially connected antenna arrays, flAn analog precoding vector representing the ith transmit sub-array of the base station,
Figure BDA0002230244020000144
wmrepresenting the mth user dimension as NrX 1 received combined vector, NrIndicating the number of receive antennas per user;
if the status of the first user is static or blocked, then
Figure BDA0002230244020000145
f1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝1，
Figure BDA0002230244020000146
Wherein phi issRepresenting the angle of departure, theta, of the user's propagating beamsThe arrival angle of the user's propagation beam is shown, the propagation beam of the user is the departure angle and the arrival angle corresponding to the path gain maximum path, wm(1: N) denotes that the mth user receives the 1 st of the merged vectorsTo the nth element, N represents the number of transmission sub-arrays used in the base station,
Figure BDA0002230244020000147
aN(x) Indicating the array steering vector, N indicating the number of linear antennas per transmit sub-array of the base station,
Figure BDA0002230244020000148
λ represents the carrier wavelength, d represents the spacing between antenna elements,
Figure BDA0002230244020000151
if the state of the first user is a moving or moving + blocking state, passing through the adjacent N in the base stationcCooperation of transmission sub-arrays, and N of userscThe receiving sub-arrays cooperate to perform beam combination to increase the beam width of the transmission beam, and the angle difference between the transmission beam of the adjacent transmitting sub-array and the transmission beam of the adjacent receiving sub-array is
Figure BDA0002230244020000152
Then it is determined that,
Figure BDA0002230244020000153
if the tangential direction to the first user's direction of motion is counterclockwise, then,
Figure BDA0002230244020000154
Figure BDA0002230244020000155
Figure BDA0002230244020000156
Figure BDA0002230244020000157
Figure BDA0002230244020000158
n＝Nc，
if the tangential direction to the first user's direction of motion is clockwise, then,
Figure BDA0002230244020000159
Figure BDA00022302440200001510
Figure BDA00022302440200001511
Figure BDA00022302440200001512
Figure BDA00022302440200001513
n＝Nc；
wherein, wm(N +1:2N) denotes that the mth user receives the N +1 th to 2N-th elements of the merged vector, wm[(Nc-1)N+1:NcN](N) th user receiving combined vectorc-1) N +1 to NthcN elements;
in the present embodiment, the adjacent subarrays N in cooperationcThe number of the transmitting subarrays is set according to the requirements of the system, the number n of the used transmitting subarrays in the base station is accumulated and calculated in the process of transmitting data to users by the base station, and the base station is shared
Figure BDA0002230244020000161
And transmitting the sub-arrays.
Updating the analog precoding matrix FRFA value of (d);
specifically, f of the first user is obtained in step ii or step iii1Or f1To fNcIs updated to the analog precoding matrix FRFIn (1).
S53: and designing analog precoding matrixes of other users to obtain the analog precoding matrixes of all the users.
Specifically, the method comprises the following steps:
i. selecting a propagation beam for the mth user, determining a departure angle and an arrival angle of the propagation beam for the mth user,
Figure BDA0002230244020000162
Figure BDA0002230244020000163
wherein | | | purple hair2The norm of 2 is taken, and min is the minimum value;
judging the status of the mth user, if the status of the mth user is static or blocking, then,
Figure BDA0002230244020000164
fn+1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝n+1，
Figure BDA0002230244020000165
if the mth user state is a moving or moving + blocking state, passing through the adjacent N in the base stationcCooperation of the transmitting sub-arrays, and N of userscThe receiving sub-arrays cooperate to carry out beam combination, and the angle difference of the transmission beams of the adjacent transmitting sub-arrays and the adjacent receiving sub-arrays is
Figure BDA0002230244020000171
Then it is determined that,
Figure BDA0002230244020000172
if the tangential direction of the mth user's motion direction is counterclockwise, then,
Figure BDA0002230244020000173
Figure BDA0002230244020000174
Figure BDA0002230244020000175
Figure BDA0002230244020000176
Figure BDA0002230244020000177
n＝n+Nc，
if the tangential direction of the mth user's motion direction is clockwise, then,
Figure BDA0002230244020000178
Figure BDA0002230244020000179
Figure BDA00022302440200001710
Figure BDA00022302440200001711
Figure BDA00022302440200001712
n＝n+Nc；
updating the analog precoding matrix FRFIs given as M +1, repeating steps i-iv, and iterating M-1 times to obtain the analog precoding matrix F for all usersRF。
S6: designing a digital pre-coding matrix of a user according to the analog pre-coding matrix;
specifically, the method comprises the following steps:
s61: constructing a combined interference matrix for the mth user
Figure BDA0002230244020000181
Figure BDA0002230244020000182
Figure BDA0002230244020000183
Wherein the content of the first and second substances,
Figure BDA0002230244020000184
the dimension of (a) is M-1 xM,
Figure BDA0002230244020000185
a valid channel vector representing the mth user;
s62: combined interference matrix to the mth user
Figure BDA0002230244020000186
The singular value decomposition is carried out, and the singular value decomposition,
Figure BDA0002230244020000187
wherein, UmRepresenting combined interference matrices
Figure BDA0002230244020000188
Left singular matrix of, VmRepresenting combined interference matrices
Figure BDA0002230244020000189
Right singular matrix ofmRepresenting a singular value diagonal matrix;
s63: according to the right singular matrix VmObtaining the digital precoding vector F of the mth userBB(:,m)，
FBB(:,m)＝Vm(:,M)，
Wherein, FBB(: m) denotes a digital precoding matrix FBBRow m of (5), Vm(M) represents the right singular matrix VmThe M-th column of (1);
s64: repeating S61-S64 by making M equal to M +1, and iterating M-1 times to obtain the analog precoding matrix F of the userBB。
S7: and carrying out data transmission according to the analog precoding matrix and the digital precoding matrix.
Specifically, the method comprises the following steps:
the base station transmits data according to the channel state information, the analog pre-coding matrix, the digital pre-coding matrix and the receiving merging vector of the user, and the receiving signal y of the mth usermIn order to realize the purpose,
Figure BDA00022302440200001810
wherein n ismRepresenting the received noise vector for the mth user and x representing the transmitted signal vector.
The millimeter wave system dynamic user hybrid precoding method of the embodiment utilizes the cooperation of the adjacent antenna subarrays and makes the interval between the departure angle and the arrival angle of the propagation beam of the adjacent subarrays be
Figure BDA00022302440200001811
In the case of (2), a smooth wide wave is synthesizedAnd the beams are formed, so that the stability of the mobile user in the millimeter wave system is improved, the time for the mobile user to move out of the beams is prolonged, and the training period of the system beams is effectively prolonged.
Example two
The present embodiment is a simulation experiment of the dynamic user hybrid precoding method of the millimeter wave system in the first embodiment. The simulation conditions of this embodiment are: the base station has 15 transmitting subarrays, each transmitting subarray has 20 antennas, the system has 5 users, each user has 3 receiving subarrays, each receiving subarray has 20 antennas, the Rice factor K is 7dB, and the transmitting power of each signal is unit power.
Referring to fig. 3 and 4, fig. 3 shows a method according to the present invention and the adjacent cooperative sub-array is NcA comparison of 2 and a beam pattern map without the method of the present invention; FIG. 4 shows a method of the present invention with N adjacent cooperative sub-arraysc2 and NcA beam pattern map comparison plot of 3. As shown in FIG. 3, N is usedcThe beam of the inventive method of 2 is much wider than the beam without the inventive method, as shown in fig. 4 using NcThe process of the present invention using N as compared with 3cThe beam of the method of the invention is wider and the top end is smoother, and by comparing fig. 3 with fig. 4, the angle difference between adjacent subarray beams is shown to be
Figure BDA0002230244020000191
In this case, a wide beam with a smoother top end can be formed by a plurality of sub-arrays.
Fig. 5-8 are graphs showing the simulation of the velocity change with time of the user in different exercise states according to the present invention. Wherein, 5 users in the system are all in motion state, and each user has only direct path, the distances from the base station for the users in the system are respectively 5 meters, 15 meters and 35 meters, the moving speed of the user in fig. 5 is 2 meters/second, the moving speed of the user in fig. 6 is 5 meters/second, the moving speed of the user in fig. 7 is 15 meters/second, and the moving speed of the user in fig. 8 is 35 meters/second. In FIG. 5, when the user moves at a speed2 m/s and 15 m from the base station, NcThe 3 beam combining algorithm (beam combining of 3 sub-arrays) can be maintained and the rate performance is kept unchanged within 1000 ms, NcThe beam combining algorithm 2 (beam combining of 2 sub-arrays) drops the sum rate to 13 bits/s/hz after 1000 ms, while the sum rate performance without beam combining algorithm drops to 5 bits/s/hz after 750 ms. As shown in fig. 6, using NcThe method of the invention of 3 can maintain the sum speed unchanged within 500 milliseconds under the condition of the distance of 15 meters and the speed of 5 meters/second. Without the method of the invention, at a distance of 15 meters and a speed of 5 m/s, and after 300 milliseconds, the sum rate decays to 5 bits/s/hz. As shown in FIG. 7, when the user is 35 m away and 15 m/s at speed, NcThe sum rate of the method of the invention decays to 5 bits/s/hz after 200 ms 2, NcThe sum rate of the method of the invention decays to 5 bits/s/hz after 300 ms, while the sum rate without the method of the invention decays to 5 bits/s/hz after 100 ms. As shown in FIG. 8, when the user speed is 35 m/s, N is usedc2 and NcThe inventive method of 3 has better sum rate performance at distances of 5 meters, 15 meters and 35 meters than the inventive method, and in particular, when the distance is 35 meters and the time is 500 milliseconds, N iscThe sum rate of the inventive method of 2 is 4.5 bits/sec/hz, NcThe sum rate for the inventive method of 3 is 15 bits/sec/hz sum, while the sum rate without the inventive method is 1.5 bits/sec/hz. As can be seen from fig. 5, 6, 7 and 8, the method of the present invention can effectively maintain the sum rate performance and slow down the rate of decrease of the sum rate, compared to the curves without the method of the present invention.
Fig. 9-12 are graphs showing the simulation of the sum rate over time when the user uses the method of the present invention in a random state and does not use the method of the present invention. Wherein, all users in the system are in random states, the states include stationary, moving, blocking and moving + blocking, each user has 5 propagation paths, the distances from the base station to the users are respectively 5 meters, 15 meters and 35 meters, the moving speed of the user in fig. 9 is 2 meters/second, the moving speed of the user in fig. 10 is 5 meters/second, the moving speed of the user in fig. 11 is 15 meters/second, and the moving speed of the user in fig. 12 is 35 meters/second. The propagation delay of the ith propagation path from the transmitting end to the receiving end follows an exponential delay profile, which is expressed by the following formula,
τ′i＝-γτDSln(Xi)，
wherein, γτWith 3 representing the delay scaling factor, DS representing the delay spread, and lgDS having an average of-0.28 log10(1+fc)-7.173，fc28GHz is the carrier frequency. Then, the delay of the ith path can be obtained by subtracting the minimum delay to normalize the delay and arranging the delays in descending order,
τi＝sort{τ′i-min(τ′i)}，
wherein sort represents a descending order function, and min represents taking a minimum value.
As can be seen from fig. 9, 10, 11 and 12, the method of the present invention can improve the sum rate performance of dynamic users compared to NcProcess according to the invention when 2, NcThe method of the present invention at 3 maintains the sum rate performance for a longer period of time. As shown in FIG. 9, when the user moves at 2 m/s and 15 m from the base station, Nc2 and NcThe adaptive hybrid precoding algorithm of 3 can maintain the system sum rate above 20 bits/s/hz in 1000 ms, while the sum rate without the adaptive hybrid precoding algorithm drops to 16 bits/s/hz after 1000 ms. As shown in FIG. 10, NcThe present method can maintain the sum rate performance of a user with a distance of 5 meters and a speed of 5 meters/second for 500 milliseconds at 3, but the sum rate without the present method is rapidly reduced. As shown in FIG. 11, when the distance is 15 m, N is used without using the method of the present inventioncMethod of the invention when 2 and using NcWhen the sum rate of the method of the present invention is 3, the sum rate decays to 14 bits/s/hz after 200 ms, 325 ms and 500 ms, respectively, it can be seen that the method of the present invention can effectively delay the decay of the sum rate performance of the dynamic users, and improve the system stability. As shown in fig. 12As shown, when the distance is 15 m, N is used without using the method of the present inventioncMethod of the invention when 2 and using NcWhen the sum rate decays to 14 bits/s/hz after 200 ms, 300 ms and 500 ms, respectively, the method of the present invention can be seen to also effectively mitigate the sum rate performance fading when the user speed is 35 m/s and the distance is 15 m. Therefore, the method of the invention can effectively delay the sum rate performance of dynamic users and improve the system stability.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (10)

1. A dynamic user hybrid precoding method for a millimeter wave system is characterized by comprising the following steps:
acquiring channel state information of each user, and storing the channel state information of adjacent time slots of each user in a storage unit;
obtaining the optimal propagation path of each user according to the stored channel state information;
judging whether the optimal propagation path is a direct path or not to obtain a judgment result of each user;
obtaining the state of each user according to the judgment result;
designing a simulation pre-coding matrix of the user according to the state of the user;
designing a digital pre-coding matrix of a user according to the analog pre-coding matrix;
and carrying out data transmission according to the analog precoding matrix and the digital precoding matrix.
2. The method of claim 1, wherein obtaining channel state information of each user and storing the channel state information of each user's adjacent time slot in a storage unit comprises:
the base station obtains the channel state information H scanned by each user in the beam training time slotm，
Figure FDA0002230244010000011
Wherein HmA channel state information matrix representing the mth user, M representing the number of users,
Figure FDA0002230244010000012
a channel state matrix representing the direct path of the mth user,
Figure FDA0002230244010000013
a channel state matrix representing the ith non-direct path of the mth user, K represents the Rice factor, PmIndicates the number of propagation paths of the mth user,
Figure FDA0002230244010000014
a complex path gain representing an ith non-direct path of an mth user;
and respectively storing the channel state information of the t-1 time slot and the t time slot of each user in a first storage unit and a second storage unit, wherein t is more than or equal to 1.
3. The method of claim 2, wherein obtaining an optimal propagation path for each user according to the stored channel state information comprises:
obtaining an optimal propagation path of each user in the first storage unit and the second storage unit according to the following formula,
Figure FDA0002230244010000021
wherein g represents a memory cell,
Figure FDA0002230244010000022
indicating the departure angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure FDA0002230244010000023
indicating the arrival angle corresponding to the optimal propagation path of the user in the g-th storage unit,
Figure FDA0002230244010000024
represents the complex path gain of the xth propagation path of the user in the g storage unit, | | | represents taking the absolute value, max represents taking the maximum value,
Figure FDA0002230244010000025
the number of propagation paths of the mth user in the g-th storage unit is shown.
4. The method of claim 3, wherein determining whether the optimal propagation path is a direct path to obtain a determination result for each user comprises:
if the optimal propagation path satisfies
Figure FDA0002230244010000026
The optimal propagation path of the user in the g-th storage unit is a direct propagation path, and if the optimal propagation path meets the requirement
Figure FDA0002230244010000027
The optimal propagation path of the user at the g-th storage unit is a non-direct path,
wherein epsilon represents the threshold values of the direct path and the indirect path, and epsilon is more than or equal to 0 and less than or equal to muAS，μASRepresents the average of the angular spread of the propagating beam.
5. The method of claim 4, wherein obtaining the status of each user according to the determination result comprises:
if the optimal propagation paths of the user in the first storage unit and the second storage unit are direct paths, executing the following steps:
if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit satisfies the requirement
Figure FDA0002230244010000028
The user is in a static state, and if the departure angle of the optimal propagation path of the user in the first storage unit and the second storage unit meets the requirement
Figure FDA0002230244010000031
The user is in a state of motion and,
wherein, delta represents the judgment threshold value of the motion and the static state of the user,
Figure FDA0002230244010000032
μASan average value representing the angular spread of the propagation beam, BW representing the half-power beamwidth of the base station propagation beam;
if the optimal propagation path of the user in the first storage unit is a direct path and the optimal propagation path in the second storage unit is a non-direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure FDA0002230244010000033
The number of the path pairs is one,
Figure FDA0002230244010000034
if it is
Figure FDA0002230244010000035
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure FDA0002230244010000036
The user is in a blocking + movement state, otherwise the user is in a blocking state;
if the optimal propagation path of the user in the first storage unit is a non-direct path and the optimal propagation path in the second storage unit is a direct path, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure FDA0002230244010000037
The number of the path pairs is one,
Figure FDA0002230244010000038
if it is
Figure FDA0002230244010000039
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure FDA00022302440100000310
The user is in a motion state, otherwise, the user is in a static state;
if the optimal propagation paths of the user in the first storage unit and the second storage unit are both non-direct paths, executing the following steps:
matching the nearest non-direct path of the user in the first storage unit and the second storage unit according to the following formula to obtain
Figure FDA0002230244010000041
The number of the path pairs is one,
Figure FDA0002230244010000042
if it is
Figure FDA0002230244010000043
The angular difference of the departure angles of at least 2 path pairs in each path pair satisfies
Figure FDA0002230244010000044
The user is in the block + move state otherwise the user is in the block state.
6. The method of claim 5, wherein designing an analog precoding matrix for a user based on the state of the user comprises:
prioritizing a status of the user, wherein,
stationary > occlusion > motion > occlusion + motion;
designing a simulation pre-coding matrix of the first user according to the priority sequence;
and designing analog precoding matrixes of other users to obtain the analog precoding matrixes of all the users.
7. The method of claim 6, wherein designing the first user's analog precoding matrix according to the prioritization comprises:
i. let m equal to 1, wmFor all-zero vectors, a precoding matrix F is simulatedRFJudging the state of the first user for an all-zero matrix, wherein the simulation precoding matrix FRFIn order to realize the purpose,
Figure FDA0002230244010000045
wherein, the analog precoding matrix FRFHas the dimension of
Figure FDA0002230244010000046
NtIndicating the number of transmit antennas of the base station,
Figure FDA0002230244010000047
transmitting subarrays representing base stationsNumber flAn analog precoding vector representing the ith transmit sub-array of the base station,
Figure FDA0002230244010000048
wmrepresenting the mth user dimension as NrX 1 received combined vector, NrIndicating the number of receive antennas per user;
if the status of the first user is static or blocked, then
Figure FDA0002230244010000051
f1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝1，
Figure FDA0002230244010000052
Wherein phi issRepresenting the angle of departure, theta, of the user's propagating beamsThe arrival angle of the user's propagation beam is shown, the propagation beam of the user is the departure angle and the arrival angle corresponding to the path gain maximum path, wm(1: N) represents the 1 st to Nth elements of the m-th user reception combining vector, N represents the number of transmission sub-arrays used in the base station,
Figure FDA0002230244010000053
aN(x) Indicating the array steering vector, N indicating the number of linear antennas per transmit sub-array of the base station,
Figure FDA0002230244010000054
λ represents the carrier wavelength, d represents the spacing between antenna elements,
Figure FDA0002230244010000055
if said first one is usedIf the state of the user is moving or moving + blocking state, the user passes through the adjacent N in the base stationcCooperation of transmission sub-arrays, and N of userscThe receiving sub-arrays cooperate to carry out beam combination, and the angle difference of the transmission beams of the adjacent transmitting sub-arrays and the adjacent receiving sub-arrays is
Figure FDA0002230244010000056
Then it is determined that,
Figure FDA0002230244010000057
if the tangential direction to the first user's direction of motion is counterclockwise, then,
Figure FDA0002230244010000058
Figure FDA0002230244010000059
Figure FDA00022302440100000510
Figure FDA0002230244010000061
Figure FDA0002230244010000062
n＝Nc，
if the tangential direction to the first user's direction of motion is clockwise, then,
Figure FDA0002230244010000063
Figure FDA0002230244010000064
Figure FDA0002230244010000065
Figure FDA0002230244010000066
Figure FDA0002230244010000067
n＝Nc；
wherein, wm(N +1:2N) denotes that the mth user receives the N +1 th to 2N-th elements of the merged vector, wm[(Nc-1)N+1:NcN](N) th user receiving combined vectorc-1) N +1 to NthcN elements;
updating the analog precoding matrix FRFThe value of (c).
8. The method of claim 7, wherein designing the analog precoding matrix of other users to obtain the analog precoding matrix of all users comprises:
i. selecting a propagation beam for the mth user, determining a departure angle and an arrival angle of the propagation beam for the mth user,
Figure FDA0002230244010000068
Figure FDA0002230244010000069
wherein | | | purple hair2The norm of 2 is taken, and min is the minimum value;
judging the status of the mth user, if the status of the mth user is static or blocking, then,
Figure FDA0002230244010000071
fn+1＝aN(φs)，
wm(1:N)＝aN(θs)，
n＝n+1，
Figure FDA0002230244010000072
if the mth user state is a moving or moving + blocking state, passing through the adjacent N in the base stationcCooperation of the transmitting sub-arrays, and N of userscThe receiving sub-arrays cooperate to carry out beam combination, and the angle difference of the transmission beams of the adjacent transmitting sub-arrays and the adjacent receiving sub-arrays is
Figure FDA0002230244010000073
Then it is determined that,
Figure FDA0002230244010000074
if the tangential direction of the mth user's motion direction is counterclockwise, then,
Figure FDA0002230244010000075
Figure FDA0002230244010000076
Figure FDA0002230244010000077
Figure FDA0002230244010000078
Figure FDA0002230244010000079
n＝n+Nc，
if the tangential direction of the mth user's motion direction is clockwise, then,
Figure FDA0002230244010000081
Figure FDA0002230244010000082
Figure FDA0002230244010000083
Figure FDA0002230244010000084
Figure FDA0002230244010000085
n＝n+Nc；
updating the analog precoding matrix FRFLet M be M +1, repeat steps i-iv, iterate M-1 times to obtain the analog precoding matrix F for all usersRF。
9. The method of claim 8, wherein designing a user's digital precoding matrix based on the analog precoding matrix comprises:
constructing a combined interference matrix for the mth user
Figure FDA0002230244010000086
Figure FDA0002230244010000087
Figure FDA0002230244010000088
Wherein the content of the first and second substances,
Figure FDA0002230244010000089
the dimension of (a) is M-1 xM,
Figure FDA00022302440100000810
a valid channel vector representing the mth user;
combined interference matrix to the mth user
Figure FDA00022302440100000811
The singular value decomposition is carried out, and the singular value decomposition,
Figure FDA00022302440100000812
wherein, UmRepresenting combined interference matrices
Figure FDA00022302440100000813
Left singular matrix of, VmRepresenting combined interference matrices
Figure FDA00022302440100000814
Right singular matrix ofmRepresenting a singular value diagonal matrix;
according to the right singular matrix VmObtaining the digital precoding vector F of the mth userBB(:,m)，
FBB(:,m)＝Vm(:,M)，
Wherein, FBB(: m) denotes a digital precoding matrix FBBRow m of (5), Vm(M) represents the right singular matrix VmThe M-th column of (1);
making m equal to m +1, repeating the above steps and iteratingObtaining the simulation pre-coding matrix F of the user M-1 timesBB。
10. The method of claim 9, wherein transmitting data according to the analog precoding matrix and the digital precoding matrix comprises:
the base station transmits data according to the channel state information, the analog pre-coding matrix, the digital pre-coding matrix and the receiving merging vector of the user, and the receiving signal y of the mth usermIn order to realize the purpose,
Figure FDA0002230244010000091
wherein n ismRepresenting the received noise vector for the mth user and x representing the transmitted signal vector.
CN201910965145.8A 2019-10-11 2019-10-11 Dynamic user hybrid precoding method for millimeter wave system Active CN110958041B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910965145.8A CN110958041B (en) 2019-10-11 2019-10-11 Dynamic user hybrid precoding method for millimeter wave system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910965145.8A CN110958041B (en) 2019-10-11 2019-10-11 Dynamic user hybrid precoding method for millimeter wave system

Publications (2)

Publication Number Publication Date
CN110958041A true CN110958041A (en) 2020-04-03
CN110958041B CN110958041B (en) 2021-08-24

Family

ID=69976354

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910965145.8A Active CN110958041B (en) 2019-10-11 2019-10-11 Dynamic user hybrid precoding method for millimeter wave system

Country Status (1)

Country Link
CN (1) CN110958041B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105162556A (en) * 2015-08-19 2015-12-16 南京邮电大学 Large-scale MIMO system channel feedback method based on spatial- temporal correlation
CN107005289A (en) * 2014-11-25 2017-08-01 高通股份有限公司 Technology for detecting line-of-sight transmission using millimetre-wave attenuator equipment
CN107294585A (en) * 2016-04-01 2017-10-24 中兴通讯股份有限公司 The feedback method and device of channel condition information
CN108111199A (en) * 2017-05-05 2018-06-01 中兴通讯股份有限公司 Feedback, method of reseptance and device, equipment, the storage medium of channel state information
US20180198511A1 (en) * 2016-11-03 2018-07-12 Futurewei Technologies, Inc. System and Method for Hierarchal Beamforming and Rank Adaptation for Hybrid Antenna Architecture
CN108471325A (en) * 2018-03-23 2018-08-31 北京理工大学 A kind of sparse radio frequency/base band mixing method for precoding
CN109617585A (en) * 2019-01-18 2019-04-12 杭州电子科技大学 Mixing method for precoding based on part connection in the extensive MIMO of millimeter wave
WO2019173775A1 (en) * 2018-03-08 2019-09-12 Cohere Technologies, Inc. Scheduling multi-user mimo transmissions in fixed wireless access systems

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107005289A (en) * 2014-11-25 2017-08-01 高通股份有限公司 Technology for detecting line-of-sight transmission using millimetre-wave attenuator equipment
CN105162556A (en) * 2015-08-19 2015-12-16 南京邮电大学 Large-scale MIMO system channel feedback method based on spatial- temporal correlation
CN107294585A (en) * 2016-04-01 2017-10-24 中兴通讯股份有限公司 The feedback method and device of channel condition information
US20180198511A1 (en) * 2016-11-03 2018-07-12 Futurewei Technologies, Inc. System and Method for Hierarchal Beamforming and Rank Adaptation for Hybrid Antenna Architecture
CN108111199A (en) * 2017-05-05 2018-06-01 中兴通讯股份有限公司 Feedback, method of reseptance and device, equipment, the storage medium of channel state information
WO2019173775A1 (en) * 2018-03-08 2019-09-12 Cohere Technologies, Inc. Scheduling multi-user mimo transmissions in fixed wireless access systems
CN108471325A (en) * 2018-03-23 2018-08-31 北京理工大学 A kind of sparse radio frequency/base band mixing method for precoding
CN109617585A (en) * 2019-01-18 2019-04-12 杭州电子科技大学 Mixing method for precoding based on part connection in the extensive MIMO of millimeter wave

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
尤若楠: "基于有限反馈的毫米波MIMO系统的混合预编码方法", 《电信科学》 *

Also Published As

Publication number Publication date
CN110958041B (en) 2021-08-24

Similar Documents

Publication Publication Date Title
CN100574144C (en) The method and apparatus that arrival direction is followed the tracks of and in the application of the beam shaping of base station
CN113411110B (en) Millimeter wave communication beam training method based on deep reinforcement learning
US6590532B1 (en) Radio device
Booth et al. Multi-armed bandit beam alignment and tracking for mobile millimeter wave communications
KR101199569B1 (en) Opportunistic interference aligned user selection in multiuser mimo interference channels
CN111245481B (en) Large-scale MIMO satellite mobile communication downlink transmission method and system
KR20050004605A (en) Combined beamforming-diversity wireless fading channel de-modulator using sub-array grouped adaptive array antennas, portable telecommunication receiving system comprising it and method thereof
CN100492935C (en) Wave beam shaping method
JP2007159130A (en) Uplink receiving method and device in distributed antenna mobile communication system
Shen et al. Design and implementation for deep learning based adjustable beamforming training for millimeter wave communication systems
CN113438002B (en) LSTM-based analog beam switching method, device, equipment and medium
CN109951216B (en) Codebook-assisted large-scale MIMO DOA estimation method
Shen et al. Multi-beam design for near-field extremely large-scale RIS-aided wireless communications
Abdallah et al. Multi-agent deep reinforcement learning for beam codebook design in RIS-aided systems
CN110958041B (en) Dynamic user hybrid precoding method for millimeter wave system
CN115549745B (en) RIS phase shift design method, RIS phase shift design device, RIS phase shift design computer equipment and RIS phase shift storage medium
CN115361043B (en) Communication control method and control system for high-speed rail millimeter wave communication system
Liu et al. Hybrid RIS-assisted MIMO dual-function radar-communication system
Chen et al. A deep learning driven hybrid beamforming method for millimeter wave MIMO system
CN114499615B (en) Near-far field unified transmitting beam forming method in terahertz communication system
Cheng et al. Lightweight CNNs-based interleaved sparse array design of phased-MIMO radar
CN113437999B (en) Adaptive beam width modulation method for inhibiting beam drift effect in millimeter wave communication system
Loadman et al. An overview of adaptive antenna technologies for wireless communications
CN115865160A (en) Beam forming method and system of large-scale MIMO-NOMA system in low-orbit satellite communication scene
CN115021779A (en) RIS assisted MIMO system discrete phase shift design method and apparatus

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