CN113328770A - Large-scale MIMO channel state estimation method and device - Google Patents
Large-scale MIMO channel state estimation method and device Download PDFInfo
- Publication number
- CN113328770A CN113328770A CN202110616813.3A CN202110616813A CN113328770A CN 113328770 A CN113328770 A CN 113328770A CN 202110616813 A CN202110616813 A CN 202110616813A CN 113328770 A CN113328770 A CN 113328770A
- Authority
- CN
- China
- Prior art keywords
- matrix
- channel
- channel state
- millimeter wave
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims abstract description 128
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 17
- 230000006870 function Effects 0.000 claims description 31
- 238000011084 recovery Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 10
- 239000013598 vector Substances 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000003190 augmentative effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- WYROLENTHWJFLR-ACLDMZEESA-N queuine Chemical compound C1=2C(=O)NC(N)=NC=2NC=C1CN[C@H]1C=C[C@H](O)[C@@H]1O WYROLENTHWJFLR-ACLDMZEESA-N 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
-
- 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
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
Abstract
The invention discloses a large-scale MIMO channel state estimation method and a device, wherein the method comprises the following steps: determining sparsity and low rank of millimeter wave channel transmission according to millimeter wave transmission in a large-scale MIMO system; under a matrix complete model, determining sparsity by adopting a beam domain space representation matrix, and determining low rank by adopting a matrix kernel norm; constructing an objective function according to sparsity and low rank of a millimeter wave channel, an unknown sparse channel gain matrix and an unknown channel state matrix; and optimizing and solving the objective function to determine the state information of the millimeter wave channel, wherein iterative updating variable solving is carried out by adopting an alternative direction multiplier method, and a random singular value decomposition method is introduced to solve a channel state matrix in the alternative direction multiplier method. By means of the low-rank characteristic and the sparse characteristic of the millimeter wave channel, a random singular value decomposition method of the matrix is introduced in the process of iteratively updating the variable, and the efficiency of obtaining channel state information is improved.
Description
Technical Field
The present invention relates to the field of wireless communication technologies, and in particular, to a method and an apparatus for estimating a large-scale MIMO channel state.
Background
The development of mobile communication technology has promoted the historical progress of human, but the development of communication technology is restricted by the complicated and changeable wireless Channel of mobile communication, and the acquisition of accurate Channel State Information (CSI) is crucial to the acquisition of diversity gain and the utilization of Channel fading. The large-scale multi-input multi-output technology is a key technology in a 5G communication system, the number of antennas in the technology can reach 64 or more than 128, compared with the number of antennas with the number not exceeding 4 of a 4G system terminal, the number is increased by one order of magnitude, the data transmission rate is improved, meanwhile, the calculation complexity is high, so that more antenna units need higher signal frequency than a microwave frequency band, therefore, millimeter waves with frequency bands distributed at 30-300 GHz can meet the requirements of people, researchers combine the millimeter waves with the large-scale multi-input multi-output technology, and the millimeter wave large-scale MIMO technology is provided.
The millimeter wave has the characteristics of large bandwidth, high frequency and short wavelength, so that the scattering of the millimeter wave is weak when the millimeter wave is transmitted in the air, and the characteristic of cluster transmission is presented, so that a millimeter wave channel looks sparse, except the sparse characteristic, the millimeter wave channel can show angular expansion in the arrival angle, the departure angle and the elevation angle domain, the angular expansion can generate a useful low-rank structure, but the millimeter wave generates very large loss when the millimeter wave is transmitted in space, is easy to be absorbed by the atmosphere and is seriously attenuated by rainfall.
In order to compensate path loss of millimeter waves, a millimeter wave communication system realizes a high-directivity antenna through a super-large-scale antenna array and a beam forming technology to obtain gain and improve transmission quality, the beam forming technology is a design that a user carries out precoding according to information such as the direction of a transmitting end or a receiving end, in a millimeter wave large-scale MIMO system, researchers usually design a mixed precoding matrix according to the obtained CSI, and can also carry out fading matching on a transmitting signal by utilizing the CSI so as to achieve efficient transmission, meanwhile, the transmitting end obtains accurate and effective channel state information and is beneficial to assisting in realizing mixed beam forming, but in a large-scale antenna array scene, the environment with low signal-to-noise ratio is usually adopted before beam forming, so that the accurate and effective channel state information is more difficult to obtain.
Disclosure of Invention
The invention aims to provide a large-scale MIMO channel state estimation method to solve the problem that channel state information is difficult to estimate.
In order to achieve the above object, the present invention provides a large-scale MIMO channel state estimation method, including:
determining sparsity and low rank of millimeter wave channel transmission according to millimeter wave transmission in a large-scale MIMO system; wherein, under a matrix-complete model,
determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm;
constructing a target function according to the sparsity and low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix and a preset unknown channel state recovery matrix;
and optimizing and solving the objective function to determine the state information of the millimeter wave channel, wherein the optimizing and solving comprises iteratively updating variables by adopting an alternative direction multiplier method, and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
Preferably, the first and second electrodes are formed of a metal,
obtaining a beam domain space representation matrix S according to the beam domain space decomposition of the low rank matrix H, as follows:
wherein D isRAnd DTEach represents a discrete fourier matrix that is,representing the discrete Fourier matrix DTThe conjugate matrix of (2).
Preferably, the objective function is as follows:
minH,S τH||H||*+τS||S||1;
wherein | H | Y*Represents the nuclear norm of the low-rank matrix H, | S | | luminance1Representing the 1 norm of the beam domain space representation matrix S, omega representing the sampling matrix, o representing the Hadamard product operation, HΩRepresenting the result of sampling the channel matrix, τHAnd τSEach representing a weighting factor related to the number of millimeter-wave channel paths.
Preferably, the objective function further comprises:
introducing two auxiliary variables to eliminate the noise of the millimeter wave channel, wherein the two auxiliary variables are as follows:
updating the objective function as follows:
The present invention also provides a large-scale MIMO channel state estimation apparatus, comprising:
the acquisition module is used for determining the sparsity and low rank of millimeter wave channel transmission according to the transmission of millimeter waves in the large-scale MIMO system; wherein, under a matrix-complete model,
determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm;
the constructing module is used for constructing a target function according to the sparsity and the low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix and a preset unknown channel state recovery matrix;
and the optimization module is used for optimizing and solving the objective function to determine the state information of the millimeter wave channel, wherein the optimization and solution comprises the steps of iteratively updating variables by adopting an alternative direction multiplier method and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
Preferably, the acquisition module is configured to acquire, from the user,
and is further configured to obtain the beam domain spatial representation matrix S according to a beam domain spatial decomposition of the low rank matrix H, as follows:
wherein D isRAnd DTEach represents a discrete fourier matrix that is,representing the discrete Fourier matrix DTThe conjugate matrix of (2).
Preferably, the constructing module is further configured to construct the objective function as follows:
minH,S τH||H||*+τS||S||1;
wherein | H | Y*Represents the nuclear norm of the low-rank matrix H, | S | | luminance1Representing the beam domain spatial representation1 norm of matrix S, Ω denotes sampling matrix, o denotes Hadamard product operation, HΩRepresenting the result of sampling the channel matrix, τHAnd τSEach representing a weighting factor related to the number of millimeter-wave channel paths.
Preferably, the constructing module is further configured to construct the objective function as follows:
introducing two auxiliary variables to eliminate the noise of the millimeter wave channel, wherein the two auxiliary variables are as follows:
updating the objective function as follows:
The invention also provides a computer terminal device comprising one or more processors and a memory. A memory coupled to the processor for storing one or more programs; when executed by the one or more processors, cause the one or more processors to implement a massive MIMO channel state estimation method as in any one of the above embodiments.
The present invention also provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the massive MIMO channel state estimation method according to any of the embodiments described above.
In the background of a complete matrix model, the method for introducing the random singular value decomposition of the matrix in the process of iteratively updating the variable improves the efficiency of acquiring the state information of the channel by the low-rank characteristic and the sparse characteristic of the millimeter wave channel.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the embodiments will be briefly described 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 that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flowchart of a massive MIMO channel state estimation method according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating normalized mean square error as a function of signal-to-noise ratio according to another embodiment of the present invention;
FIG. 3 is a schematic diagram of achievable spectral efficiency as a function of signal-to-noise ratio according to yet another embodiment of the present invention;
FIG. 4 is a diagram illustrating normalized mean square error as a function of iteration count for an algorithm according to an embodiment of the present invention;
fig. 5 is a schematic structural diagram of a massive MIMO channel state estimation apparatus according to an embodiment of the present invention.
Detailed Description
The technical solutions in the present invention will be described clearly and completely with reference to the accompanying drawings, and it is obvious that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be understood that the step numbers used herein are for convenience of description only and are not used as limitations on the order in which the steps are performed.
It is to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the specification of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of the described 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.
The term "and/or" refers to and includes any and all possible combinations of one or more of the associated listed items.
Referring to fig. 1, the present invention provides a large-scale MIMO channel state estimation method, which includes:
s101, determining sparsity and low-rank property of millimeter wave channel transmission according to millimeter wave transmission in a large-scale MIMO system; and under a complete matrix model, determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm.
Specifically, under a complete matrix model, channel estimation can be modeled as a recovery problem of a low-rank matrix, and meanwhile, beam space representation of the channel matrix can be introduced to be used as auxiliary information for joint recovery. Single-user multi-antenna time division duplex uplink scene, analog beam forming structure of sampling switch mode of base station end and user end, and user end configuration NTRoot antenna, base station side equipped with NRThe method comprises the following steps of root-fire antennas, wherein the antenna arrangement mode is a uniform linear array, the antenna spacing is d, the antenna departure angle is phi, the antenna arrival angle is theta, the transmitted training sequence is unit power, a high-frequency millimeter wave channel is adopted, and a classical S-V channel model is used for modeling the millimeter wave channel, and the method comprises the following steps:
wherein alpha islDenotes the complex gain of the l-th path, aR(θl) And aT(φl) Array response vectors for the ue and the bs respectively,representing a vector of response to an array aT(φl) And (4) making a conjugation device.
The problem of recovery of the low-rank matrix H under the complete matrix model can be modeled as follows:
min||H||*
s.t.ρΩ(H)=ρΩ(D)
wherein the nuclear normRepresenting the sum of the singular values of the matrix, Ω being the set of observation samples, ρΩ():CM×N→CM×NRepresents a projection matrix, defined as:
the beam domain spatial decomposition of the low rank matrix H is as follows:
wherein D isRAnd DTEach represents a discrete fourier matrix that is,representing the discrete Fourier matrix DTThe beam space representation matrix S may reflect the sparse characteristics of the millimeter wave channel.
Under the background of a complete matrix model, modeling low-rank characteristics generated by millimeter wave channel angle expansion into a recovery problem of a low-rank matrix, and establishing a constrained target function by using a matrix kernel norm as optimal convex approximation of a rank function.
S102, constructing a target function according to the sparsity and low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix and a preset unknown channel state recovery matrix.
In the design of the millimeter wave channel estimation algorithm, because the millimeter wave channel has the sparse characteristic, the millimeter wave channel matrix is expressed by the beam domain space as follows:
wherein D isRAnd DTEach represents a discrete fourier matrix that is,representing the discrete Fourier matrix DTThe beam space representation matrix S can reflect the sparse characteristic of a millimeter wave channel and only contains a small number of high channel gains, because the millimeter wave channel angular expansion can represent a useful low-rank structure, the recovery problem of the low-rank matrix can be established under a matrix complete model, the beam domain space representation of the matrix is used as auxiliary information of the recovery problem of the low-rank matrix, an unknown channel state matrix H and the beam space representation thereof are jointly recovered through the unknown sparse channel gain matrix S, weighting factors related to the number of millimeter wave channel paths are introduced, and a target convex optimization problem is formulated as follows:
minH,S τH||H||*+τS||S||1;
wherein | H | Y*Represents the nuclear norm of the low-rank matrix H, | S | | luminance1Representing the 1 norm of the beam domain space representation matrix S, omega representing the sampling matrix, o representing the Hadamard product operation, HΩRepresenting the result of sampling the channel matrix, τHAnd τSEach representing a weighting factor related to the number of millimeter-wave channel paths.
The convex optimization problem is a dual-target convex optimization problem, so that the convex optimization problem has a global optimal solution, the solution is carried out by adopting an alternating direction multiplier method, and two auxiliary variables are introduced to eliminate the noise of a millimeter wave channel, wherein the two auxiliary variables are as follows:
meanwhile, in order to make the target problem separable in variables, the target problem is updated as follows:
S103, optimizing and solving the objective function to determine state information of the millimeter wave channel, wherein the optimizing and solving include iteratively updating variables by adopting an alternative direction multiplier method, and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
Specifically, the target problem in step S103 is converted into an augmented lagrange function as follows:
wherein,the lagrange multiplier is represented by a number of lagrange multipliers,representation matrix Z1The conjugate transpose of (a) is performed,representation matrix Z2Tr (·) represents the tracing of the matrix, and ρ represents the step size of the lagrange multiplier method.
Splitting an augmented Lagrange function of a target problem into a plurality of independent subproblems according to an alternating direction multiplier method, and solving the problems when the L-th is 0,1, 2.
For the variables H, H(L+1)=Udiag({sign(ξi)max(ξi,0)}1≤i≤r)VHWherein a matrix of random singular value decompositions is usedObtaining approximate i truncated singular vector matrixes U and V and the first i larger singular values xiiThe specific calculation flow is as follows:
1) a random matrix omega is generated.
2) And taking the matrix to be decomposed as A, projecting the A onto a random matrix omega, calculating a subspace A omega of the matrix A, and generating an orthogonal basis of the subspace.
3) The iterative update uses the orthogonal bases to generate a new subspace of matrix a and a new orthogonal base, and adds together the subspaces of each iteration into a larger subspace.
4) The iteration is completed, and the orthogonal base Q of the maximum subspace is calculated
5) Let final subspace B be QTA, then calculating the singular vector matrix and singular values of the matrix B. For the variables Y, Y(L+1)=unvec(y(L+1)),
"vec.)" means that the columns of a matrix are connected as a vector, and "unwec.)" means the inverse operation thereof.
For the variable S, S(L+1)=unvec(s(L+1)),
s(L+1)=sign(Re(υ(L+1)))°max(|Re(υ(L+1))|-τ′s,0)+sign(Im(υ(L+1)))°max(|Im(υ(L+1))|-τ′s0), "Im (. -)" means taking the imaginary part.
According to experimental simulation verification, an environment is configured, the number of antennas is set to be 64 multiplied by 64, the number of channel clusters is set to be 2, the number of paths in each cluster is set to be 1, and path gain a is set to beilThe array response vector is generated as a uniform linear array, subject to a rayleigh distribution with a mean of 0 and a variance of 1.
Normalized mean square error:"ε.)" means the mean value obtained by averaging the results of multiple Monte Carlo iterations in an experiment.
Referring to fig. 2, fig. 3 and fig. 4, compared to other methods for acquiring channel information, the method of the present invention can estimate the state information of the channel more accurately and obtain good spectrum efficiency.
The invention adopts the sparse characteristic of a millimeter wave channel and a low-rank structure generated by angular expansion, establishes a dual-target convex problem represented by an unknown channel state matrix H and a beam space thereof through joint recovery of an unknown sparse channel gain matrix S under the background of a complete matrix model, obtains a global optimal solution by using an alternative direction multiplier method, and uses random singular value decomposition in the process of solving the sub-problem. The invention adopts a random singular value decomposition method to solve variables in the subproblems, avoids directly decomposing a large matrix by solving a subspace with lower dimensionality, effectively reduces the complexity in a large-scale multi-antenna scene, and improves the efficiency of obtaining channel state information by jointly utilizing the low-rank characteristic and the sparse characteristic of a millimeter wave channel and combining a matrix complete model.
Referring to fig. 4, the present invention provides a massive MIMO channel state estimation apparatus, including:
the acquisition module 11 is configured to determine sparsity and low rank of millimeter wave channel transmission according to transmission of millimeter waves in the massive MIMO system; wherein, under a matrix-complete model,
and determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm.
And the constructing module 12 is configured to construct an objective function according to the sparsity and low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix, and a preset unknown channel state recovery matrix.
And the optimization module 13 is configured to optimize and solve the objective function to determine the state information of the millimeter wave channel, where the optimizing and solving includes iteratively updating variables by using an alternative direction multiplier method, and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
For specific limitations of the massive MIMO channel state estimation apparatus, reference may be made to the above limitations, which are not described herein again. The modules in the massive MIMO channel state estimation apparatus may be wholly or partially implemented by software, hardware, and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
The invention provides a computer terminal device comprising one or more processors and a memory. A memory is coupled to the processor for storing one or more programs, which when executed by the one or more processors, cause the one or more processors to implement the massive MIMO channel state estimation method as in any one of the embodiments above.
The processor is used for controlling the overall operation of the computer terminal equipment so as to complete all or part of the steps of the massive MIMO channel state estimation method. The memory is used to store various types of data to support the operation at the computer terminal device, which data may include, for example, instructions for any application or method operating on the computer terminal device, as well as application-related data. The Memory may be implemented by any type of volatile or non-volatile Memory device or combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk, or optical disk.
In an exemplary embodiment, the computer terminal Device may be implemented by one or more Application Specific 1 integrated circuits (AS 1C), a Digital Signal Processor (DSP), a Digital Signal Processing Device (DSPD), a Programmable Logic Device (PLD), a Field Programmable Gate Array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, for performing the above-mentioned large-scale MIMO channel state estimation method, and achieving technical effects consistent with the above-mentioned methods.
In another exemplary embodiment, a computer readable storage medium comprising program instructions which, when executed by a processor, implement the steps of the massive MIMO channel state estimation method in any one of the above embodiments is also provided. For example, the computer-readable storage medium may be the above-mentioned memory including program instructions executable by a processor of a computer terminal device to perform the above-mentioned massive MIMO channel state estimation method and achieve the technical effects consistent with the above-mentioned method.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (10)
1. A massive MIMO channel state estimation method is characterized by comprising the following steps:
determining sparsity and low rank of millimeter wave channel transmission according to millimeter wave transmission in a large-scale MIMO system; wherein, under a matrix-complete model,
determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm;
constructing a target function according to the sparsity and low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix and a preset unknown channel state recovery matrix;
and optimizing and solving the objective function to determine the state information of the millimeter wave channel, wherein the optimizing and solving comprises iteratively updating variables by adopting an alternative direction multiplier method, and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
2. The massive MIMO channel state estimation method of claim 1, wherein the beam domain spatial representation matrix S is obtained according to a beam domain spatial decomposition of a low rank matrix H, as follows:
3. The massive MIMO channel state estimation method of claim 2, wherein the objective function is as follows:
minH,S τH||H||*+τS||S||1;
wherein | H | Y*Represents the nuclear norm of the low-rank matrix H, | S | | luminance1Representing the 1 norm of the beam domain space representation matrix S, omega representing the sampling matrix, o representing the Hadamard product operation, HΩRepresenting the result of sampling the channel matrix, τHAnd τSEach representing a weighting factor related to the number of millimeter-wave channel paths.
4. The massive MIMO channel state estimation method of claim 3, wherein the objective function further comprises:
introducing two auxiliary variables to eliminate the noise of the millimeter wave channel, wherein the two auxiliary variables are as follows:
updating the objective function as follows:
5. A massive MIMO channel state estimation apparatus, comprising:
the acquisition module is used for determining the sparsity and low rank of millimeter wave channel transmission according to the transmission of millimeter waves in the large-scale MIMO system; wherein, under a matrix-complete model,
determining the sparsity by adopting a preset beam domain space representation matrix, and determining the low rank by adopting a preset matrix kernel norm;
the constructing module is used for constructing a target function according to the sparsity and the low rank of the millimeter wave channel, a preset unknown sparse channel gain matrix and a preset unknown channel state recovery matrix;
and the optimization module is used for optimizing and solving the objective function to determine the state information of the millimeter wave channel, wherein the optimization and solution comprises the steps of iteratively updating variables by adopting an alternative direction multiplier method and solving a channel state matrix in the alternative direction multiplier method by introducing a random singular value decomposition method.
6. The massive MIMO channel state estimation apparatus of claim 5, wherein the obtaining module is further configured to obtain the beam-domain spatial representation matrix S according to a beam-domain spatial decomposition of a low rank matrix H, as follows:
7. The massive MIMO channel state estimation apparatus of claim 6, wherein the constructing module is further configured to construct the objective function as follows:
minH,SτH||H||*+τS||S||1;
wherein | H | Y*Represents the nuclear norm of the low-rank matrix H, | S | | luminance1Representing the 1 norm of the beam domain space representation matrix S, omega representing the sampling matrix, o representing the Hadamard product operation, HΩRepresenting the result of sampling the channel matrix, τHAnd τSEach representing a weighting factor related to the number of millimeter-wave channel paths.
8. The massive MIMO channel state estimation apparatus of claim 7, wherein the constructing module is further configured to construct the objective function as follows:
introducing two auxiliary variables to eliminate the noise of the millimeter wave channel, wherein the two auxiliary variables are as follows:
updating the objective function as follows:
9. A computer terminal device, comprising:
one or more processors;
a memory coupled to the processor for storing one or more programs;
when executed by the one or more processors, cause the one or more processors to implement the massive MIMO channel state estimation method of any one of claims 1 to 4.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the massive MIMO channel state estimation method according to any one of claims 1 to 4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110616813.3A CN113328770A (en) | 2021-06-02 | 2021-06-02 | Large-scale MIMO channel state estimation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110616813.3A CN113328770A (en) | 2021-06-02 | 2021-06-02 | Large-scale MIMO channel state estimation method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113328770A true CN113328770A (en) | 2021-08-31 |
Family
ID=77421524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110616813.3A Pending CN113328770A (en) | 2021-06-02 | 2021-06-02 | Large-scale MIMO channel state estimation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113328770A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114726686A (en) * | 2022-03-24 | 2022-07-08 | 南京邮电大学 | Uniform area array millimeter wave large-scale MIMO channel estimation method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110518946A (en) * | 2019-08-30 | 2019-11-29 | 电子科技大学 | Based on the sparse time-varying channel estimation method with low-rank of millimeter wave time varying channel block |
CN112436872A (en) * | 2020-11-02 | 2021-03-02 | 北京邮电大学 | Multi-user large-scale MIMO channel estimation method and device |
CN112737649A (en) * | 2020-12-25 | 2021-04-30 | 杭州电子科技大学 | Millimeter wave channel estimation method based on angle grid optimization and norm constraint |
-
2021
- 2021-06-02 CN CN202110616813.3A patent/CN113328770A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110518946A (en) * | 2019-08-30 | 2019-11-29 | 电子科技大学 | Based on the sparse time-varying channel estimation method with low-rank of millimeter wave time varying channel block |
CN112436872A (en) * | 2020-11-02 | 2021-03-02 | 北京邮电大学 | Multi-user large-scale MIMO channel estimation method and device |
CN112737649A (en) * | 2020-12-25 | 2021-04-30 | 杭州电子科技大学 | Millimeter wave channel estimation method based on angle grid optimization and norm constraint |
Non-Patent Citations (1)
Title |
---|
邱佳锋: "基于矩阵完备的低复杂度毫米波大规模MIMO信道估计", 《电信科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114726686A (en) * | 2022-03-24 | 2022-07-08 | 南京邮电大学 | Uniform area array millimeter wave large-scale MIMO channel estimation method |
CN114726686B (en) * | 2022-03-24 | 2023-07-14 | 南京邮电大学 | Uniform area array millimeter wave large-scale MIMO channel estimation method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | Uplink cascaded channel estimation for intelligent reflecting surface assisted multiuser MISO systems | |
Qi et al. | Integrating sensing, computing, and communication in 6G wireless networks: Design and optimization | |
Bellili et al. | Generalized approximate message passing for massive MIMO mmWave channel estimation with Laplacian prior | |
Kaushik et al. | Dynamic RF chain selection for energy efficient and low complexity hybrid beamforming in millimeter wave MIMO systems | |
Tsai et al. | Efficient compressive channel estimation for millimeter-wave large-scale antenna systems | |
Eliasi et al. | Low-rank spatial channel estimation for millimeter wave cellular systems | |
CN109104225B (en) | Large-scale MIMO beam domain multicast transmission method with optimal energy efficiency | |
Kaushik et al. | Joint bit allocation and hybrid beamforming optimization for energy efficient millimeter wave MIMO systems | |
CN110838859B (en) | High-energy-efficiency robust precoding method suitable for multi-beam satellite communication system | |
Wang et al. | Channel parameter estimation of mmWave MIMO system in urban traffic scene: A training channel-based method | |
Song et al. | Coordinated hybrid beamforming for millimeter wave multi-user massive MIMO systems | |
CN112260737B (en) | Multi-beam satellite communication robust precoding method with total energy efficiency and minimum energy efficiency balanced | |
Yang et al. | Active 3D double-RIS-aided multi-user communications: Two-timescale-based separate channel estimation via Bayesian learning | |
CN109831233B (en) | Multi-cell coordinated large-scale MIMO beam domain multicast power distribution method | |
Zhang et al. | Robust beamforming for coherent signals based on the spatial-smoothing technique | |
US11252045B1 (en) | Processing blind beamforming for multi-user multiple-input multiple-output (MU-MIMO) systems | |
Zhang et al. | Hybrid precoder and combiner design for single-user mmWave MIMO systems | |
CN112469119A (en) | Positioning method, positioning device, computer equipment and storage medium | |
Hung | Modified particle swarm optimization structure approach to direction of arrival estimation | |
CN106130938B (en) | Multi-user joint channel estimation method for TDD large-scale MIMO system | |
CN113328770A (en) | Large-scale MIMO channel state estimation method and device | |
CN112235022B (en) | Low-complexity large-scale MIMO low-orbit satellite precoding method | |
Hu et al. | PRINCE: A pruned AMP integrated deep CNN method for efficient channel estimation of millimeter-wave and terahertz ultra-massive MIMO systems | |
CN110086734B (en) | Channel tracking and recovery method in time-varying large-scale MIMO system | |
Antreich et al. | Two-dimensional channel parameter estimation for millimeter-wave systems using Butler matrices |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210831 |
|
RJ01 | Rejection of invention patent application after publication |