CN113364501A - Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel - Google Patents
Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel Download PDFInfo
- Publication number
- CN113364501A CN113364501A CN202110626350.9A CN202110626350A CN113364501A CN 113364501 A CN113364501 A CN 113364501A CN 202110626350 A CN202110626350 A CN 202110626350A CN 113364501 A CN113364501 A CN 113364501A
- Authority
- CN
- China
- Prior art keywords
- power control
- follows
- channel
- user
- low
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 26
- 235000009566 rice Nutrition 0.000 title claims abstract description 26
- 240000007594 Oryza sativa Species 0.000 title 1
- 238000005457 optimization Methods 0.000 claims abstract description 39
- 241000209094 Oryza Species 0.000 claims abstract description 25
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 230000001413 cellular effect Effects 0.000 claims abstract description 11
- 238000012549 training Methods 0.000 claims abstract description 11
- 238000011217 control strategy Methods 0.000 claims abstract description 6
- 238000013139 quantization Methods 0.000 claims description 33
- 239000000654 additive Substances 0.000 claims description 21
- 230000000996 additive effect Effects 0.000 claims description 21
- 238000005562 fading Methods 0.000 claims description 18
- 238000011282 treatment Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000005265 energy consumption Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction 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
- H04B7/0452—Multi-user MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W52/00—Power management, e.g. TPC [Transmission Power Control], power saving or power classes
- H04W52/04—TPC
- H04W52/18—TPC being performed according to specific parameters
- H04W52/24—TPC being performed according to specific parameters using SIR [Signal to Interference Ratio] or other wireless path parameters
- H04W52/241—TPC being performed according to specific parameters using SIR [Signal to Interference Ratio] or other wireless path parameters taking into account channel quality metrics, e.g. SIR, SNR, CIR, Eb/lo
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The invention provides a power control method for a large-scale MIMO system based on low-precision ADC de-cellular under a Rice channel, which comprises the following steps: s1, establishing a model for removing the cellular massive MIMO system based on the low-precision ADC under the Rice channel; s2, performing uplink pilot training based on the constructed model; s3, carrying out downlink data transmission based on the constructed model; s4, analyzing the downlink reachable rate of the user; s5, constructing an optimization problem of maximizing the total rate of the system according to the downlink reachable rate; and S6, designing a power control strategy according to the optimization problem. The technical scheme of the invention further improves the total rate of the large-scale MIMO system based on the low-precision ADC in the Rice channel on the premise of ensuring that the service quality of each user is not less than the set minimum reachable rate and the sending power of each AP is not more than the set maximum sending power.
Description
Technical Field
The invention relates to the technical field of wireless communication, in particular to a power control method for a large-scale MIMO system based on low-precision ADC (analog-to-digital converter) de-cellular under a Rice channel.
Background
In recent years, with the rapid development of the mobile communication industry, the demand for data traffic has increased. In order to meet various diversified communication demands of people, a large-scale cellular Multiple Input Multiple Output (MIMO) system, which has made a thorough change in the architecture of a cellular network, gradually gets into the field of view of people and becomes one of core technologies that people should pay close attention to in the Beyond 5G and 6G era. The large-scale MIMO system breaks the division of cell boundaries on the basis of distributed large-scale MIMO, namely, a large number of Access Points (APs) and users are distributed in a wide area, different APs are connected with a Central Processing Unit (CPU) through backhaul links, and under the same time-frequency resource, the CPU instructions are listened to provide services for all users. Because the large-scale de-cellular MIMO system greatly reduces the distance between the user and the AP, the large-scale de-cellular MIMO system has stronger space macro diversity gain and the capability of resisting path loss, and can greatly improve the service quality of edge users.
However, the cellular massive MIMO system usually needs to deploy a large number of APs and users, which inevitably bring high hardware cost and huge energy consumption when the APs and users are configured with full-precision analog-to-digital converters (ADCs) to perform quantization processing on the received information. To cope with this problem, quantization processing by configuring ADCs of low precision at the AP and the user is undoubtedly a relatively straightforward solution. In addition, in many future communication scenarios, in order to cope with increasingly tight spectrum resources, the cellular massive MIMO system is likely to operate in the millimeter wave frequency band. Because millimeter waves have the characteristics of short wavelength and strong directivity, the line-of-sight component can occupy a dominant position in the whole channel, and therefore, the research of removing a cellular massive MIMO system based on a low-precision ADC under a Rice fading channel becomes a research hotspot of the current academic community.
For a large-scale cellular MIMO system based on a low-precision ADC in a Rice fading channel, the existing research only uses a maximum-minimum power control method to improve the reachable rate of the user with the worst performance, so as to ensure that the service quality of all users is equivalent. For modern communication systems, the total rate of the system is also a very important performance parameter, but the max-min power control method only puts the optimization center on improving the quality of service of the worst user, and is still deficient in improving the total rate of the system.
Disclosure of Invention
In light of the above-mentioned technical problems, the present invention provides a power control method for a large-scale cellular MIMO system based on a low-precision ADC in a rice channel. The invention further improves the total rate of the large-scale MIMO system based on the low-precision ADC in the Rice channel on the premise of ensuring that the service quality of each user is not less than the set minimum reachable rate and the sending power of each AP is not more than the set maximum sending power.
The technical means adopted by the invention are as follows:
a power control method for a large-scale MIMO system based on low-precision ADC de-cellular under a Rice channel comprises the following steps:
s1, establishing a model for removing the cellular massive MIMO system based on the low-precision ADC under the Rice channel;
s2, performing uplink pilot training based on the constructed model;
s3, carrying out downlink data transmission based on the constructed model;
s4, analyzing the downlink reachable rate of the user;
s5, constructing an optimization problem of maximizing the total rate of the system according to the downlink reachable rate;
and S6, designing a power control strategy according to the optimization problem.
Further, the model constructed in step S1 is specifically represented as follows:
wherein, gmnRepresenting the channel coefficient between the mth access point AP and the nth user; gamma raymnRepresenting large scale fading coefficients forReflecting the influence of path loss and shadow fading on the channel coefficient; the part in parentheses represents the small scale fading coefficient, represented by the line-of-sight componentAnd a scattered component hmnCN (0,1) in whichmn~[-π,π]Represents the angle of arrival and CN (0,1) represents a circularly symmetric complex gaussian variable with a mean of 0 and a variance of 1; in addition, kmnRepresenting the Rice K-factor, representing the ratio between the line-of-sight component and the scatter component, is definedAndwherein beta ismn=γmn/(1+kmn) Thus, the channel coefficients are re-represented asAnd obeyThe statistical distribution of (c).
Further, the specific implementation process of step S2 is as follows:
s21, let the length of the correlation interval be T, and let the time occupied by the pilot training in each correlation interval be taup(τp≥N);
S22, supposing that in the uplink pilot training phase, the nth user can send the uplink pilot sequence to all the APAnd that the pilot sequences transmitted between different users are mutually orthogonal, i.e.The uplink pilot information received by the mth access point AP is represented as:
wherein p ispRepresents the normalized signal-to-noise ratio of the uplink pilot signal,representing additive white gaussian noise received by the mth AP;
s23, modeling the quantization process using an additive quantization noise model, and then re-expressing the uplink pilot information received by the mth AP as:
wherein alpha ismIndicating the ADC accuracy at the mth access point AP and the accuracy corresponding to the number of quantization bits ρ used by the mth access point APmRelated, equivalent number of bits ρ m1,2,3,4,5, alphamIs precisely given as alpham0.6366,0.8825,0.96546,0.990503, 0.997501; when rhom>At 5 time, αmAnd rhomSatisfyThe approximate relationship of (a) to (b),representing and receiving a signal ym,pExtraneous additive quantization noise;
s25, for the mth AP to acquireChannel state information from the nth user will beAndmultiplying to obtain:
s26, since the access point AP can obtain the line-of-sight component information and the large-scale fading coefficient from different users, in the channel estimation stage, the system may eliminate the influence of the line-of-sight component, that is, equation (5) may be represented again as:
wherein,representing the quantization noise remaining after eliminating the influence of the line-of-sight component;
s28, obtaining g by using the minimum mean square error channel estimation scheme to the formula (6)mnThe channel estimation information of (2) is as follows:
s29, obtaining the channel estimation scheme according to the properties of the MMSE channel estimation schemeWherein
Further, the specific implementation process of step S3 is as follows:
s31, after the access point AP acquires the CSI from the user, precoding the downlink transmission signal in a conjugate transpose manner, so that the transmission signal of the mth access point AP is represented as:
wherein x isnRepresenting data signals destined for the nth user, satisfiesThe conditions of (a); p is a radical ofdIndicating the normalized SNR of the downlink transmission data signal; etamnRepresents a power control coefficient;
s32, adjusting power control coefficient etamnEach access point AP is satisfiedThe power constraint of (c), namely:
s33, according to equation (9), after all the access points AP have sent the downlink data signal, the received signal of the nth user is represented as:
wherein, wnCN (0,1) represents additive white Gaussian noise received by the nth user;
s34, based on the low-precision ADC used by the user end, representing the quantized received signal of the nth user as:
wherein, munIndicates the ADC precision of the nth user, which is equal to the quantization bit number sigma used by the nth usernIn connection with this, the present invention is,representing and receiving a signal rnUncorrelated additive quantization noise;
further, the specific implementation process of step S4 is as follows:
s41, according to the formula (12), the received signal of the nth user is represented again as follows:
s42, obtaining the expression of the downlink reachable speed of the nth user according to the formula (14) as follows:
s43, deriving the achievable rate closed expression of the nth user according to formula (15):
wherein q isnTo representThe nth column of (1);andrespectively representAnd(N-1). times.N + N1A row; the corresponding elements in the matrices Q, U and V are given as [ Q [ ]]mn=ηmn, δ (m, n) represents a dirac δ function, i.e., when m is equal to n, δ (m, n) is equal to 1; when m ≠ n, δ (m, n) is 0.
Further, the specific implementation process of step S5 is as follows:
s51, ensuring the service quality of each user not less than the set minimum reachable rateAnd the transmission power of each AP is not more than the set maximum transmission power pdOn the premise of (2), the optimization problem of maximizing the total rate of the system is summarized into the following form:
s52, optimizing the problem P1The constraint (17b) in (2) is mathematically transformed to obtain:
wherein the constraint (17b) is equivalent to a second order cone in equation (18), i.e. equation (17b) is converted to a convex constraint in equation (18); both equations (17c) and (17d) are convex constraints;
s53, adopting a continuous convex approximation method to optimize the non-convex problem P1Approximate equivalence to optimization problem P2Form of (2), optimization problem P2The concrete form of (A) is as follows:
wherein l ═ l1,…,lN]And t ═ t1,…,tN];
S54, for N ═ 1, …, N, the constraints given are as follows:
wherein the optimization problem P2The objective function (19a) in (1) is convex, but the constraint (20) is non-convex;
s55, converting the constraint (20) into a convex function form, and equating the constraint (20) to N of 1, …, N as follows:
s56, when N is 1, …, N, re-expressing equation (21) in the form:
wherein Q iskAndrespectively represent Q and lnThe iteration values after k SCA treatments, based on which equation (20) has been successfully converted to the constraint conditions of a second order cone by equation (23), and in addition, by observing the optimization problem P2The resulting constraint (19c) is convex;
s57, according to the method of successive convex approximation, approximating the formula (22) by using a first order Taylor expansion in the formula (24);
s58, according to the method of continuous convex approximation, using inequality ln (x) equal to or more than 1-x-1Converting the constraint (19c) into a constraint having the form of a second order cone as follows:
s59, according to the formula (23) and the formula (25), at the k-th iteration, the problem P is optimized2Second order cone expressed in the formPlanning the problem:
further, the specific implementation process of step S6 is as follows:
s61, initializing and setting the performance parameters: initial iteration number K is 1, maximum iteration number K, tolerance valueMinimum achievable rate for userInitial value Q1、l1And t1;
S62, solving the optimization problem in the formula (26) through a convex optimization solver, and enabling Q*、l*And t*Is the solution of the iteration;
s64, definition k ═ k +1, Qk=Q*,lk=l*And tk=t*The process goes back to step S62.
Compared with the prior art, the invention has the following advantages:
1. compared with the de-cellular massive MIMO system configured with the full-precision ADC under the Rice fading channel, the power control method based on the low-precision ADC de-cellular massive MIMO system under the Rice channel provided by the invention has the advantage that the low-precision ADC configured in the invention can greatly reduce the high hardware cost and the huge energy consumption brought by the full-precision ADC.
By adopting the power control algorithm based on the SCA mode, the power control coefficient with the maximum total system rate can be iterated under the conditions of ensuring the service quality of each user and meeting the limitation of the sending power of each AP.
2. The power control method based on the low-precision ADC de-cellular large-scale MIMO system under the Rice channel provided by the invention designs the power control strategy based on the continuous convex approximation mode, and can iterate the power control coefficient with the maximum total rate of the system under the conditions of ensuring the service quality of each user and meeting the transmission power limit of each AP.
For the above reasons, the present invention can be widely applied to the fields of wireless communication and the like.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a model structure of a massive MIMO system based on low-precision ADC de-cellular under rice channel for which the scheme of the present invention is applicable.
Fig. 3 is a simulation diagram of the system rate performance under different iteration times according to the scheme of the invention.
Fig. 4 is a simulation diagram of a cumulative distribution function of the system rate performance using the power control method in the solution of the present invention and without using the power control algorithm.
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the 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 noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in sequences other than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
As shown in fig. 1, the present invention provides a power control method for a large-scale MIMO system based on low-precision ADC de-cellular under rice channel, comprising the following steps:
s1, establishing a model for removing the cellular massive MIMO system based on the low-precision ADC under the Rice channel;
s2, performing uplink pilot training based on the constructed model;
s3, carrying out downlink data transmission based on the constructed model;
s4, analyzing the downlink reachable rate of the user;
s5, constructing an optimization problem of maximizing the total rate of the system according to the downlink reachable rate;
and S6, designing a power control strategy according to the optimization problem.
In specific implementation, as a preferred embodiment of the present invention, a decellularized massive MIMO system under a leis fading channel is established as shown in fig. 2, wherein M single-antenna APs and N (M > N) single-antenna users are randomly distributed in a wide service area, and all APs are connected to a CPU through a backhaul link and provide services to all users. In addition, in order to effectively reduce the high hardware cost and the huge energy consumption brought by the full-precision ADC quantizer, in the present embodiment, it is assumed that low-precision ADC quantizers are configured at the AP and the user to quantize the received data. The model constructed in step S1 is specifically represented as follows:
wherein, gmnRepresenting the channel coefficient between the mth access point AP and the nth user; gamma raymnRepresenting a large-scale fading coefficient for reflecting the influence of path loss and shadow fading on a channel coefficient; the part in parentheses represents the small scale fading coefficient, represented by the line-of-sight componentAnd a scattered component hmnCN (0,1) in whichmn~[-π,π]Represents the angle of arrival and CN (0,1) represents a circularly symmetric complex gaussian variable with a mean of 0 and a variance of 1; in addition, kmnRepresenting the Rice K-factor, representing the ratio between the line-of-sight component and the scatter component, is definedAndwherein beta ismn=γmn/(1+kmn) Thus, the channel coefficients are re-represented asAnd obeyThe statistical distribution of (c).
In specific implementation, as a preferred embodiment of the present invention, the step S2 performs uplink pilot training based on the constructed model, and the specific implementation process is as follows:
s21, let the length of the correlation interval be T, and let the time occupied by the pilot training in each correlation interval be taup(τp≥N);
S22, supposing that in the uplink pilot training phase, the nth user can send the uplink pilot sequence to all the APAnd that the pilot sequences transmitted between different users are mutually orthogonal, i.e.The uplink pilot information received by the mth access point AP is represented as:
wherein p ispRepresents the normalized signal-to-noise ratio of the uplink pilot signal,representing additive white gaussian noise received by the mth AP;
s23, since the AP is equipped with a low-precision ADC, its quantization error is affected by the received signal strength. In order to explore the influence of the low-precision ADC on the system performance, in this embodiment, an Additive Quantization Noise Model (AQNM) is used to model the quantization process, and then the uplink pilot information received by the mth access point AP is represented again as:
wherein alpha ismIndicating the ADC accuracy at the mth access point AP and the accuracy corresponding to the number of quantization bits ρ used by the mth access point APmRelated, equivalent number of bits ρm1,2,3,4,5, alphamIs precisely given as alpham0.6366,0.8825,0.96546,0.990503, 0.997501; when rhom>At 5 time, αmAnd rhomSatisfyThe approximate relationship of (a) to (b),representing and receiving a signal ym,pExtraneous additive quantization noise;
s25, the mth AP will acquire the Channel State Information (CSI) from the nth userAndmultiplying to obtain:
s26, since the access point AP can obtain the line-of-sight component information and the large-scale fading coefficient from different users, in the channel estimation stage, the system may eliminate the influence of the line-of-sight component, that is, equation (5) may be represented again as:
wherein,representing the quantization noise remaining after eliminating the influence of the line-of-sight component;
s28, obtaining g by using the minimum mean square error channel estimation scheme to the formula (6)mnThe channel estimation information of (2) is as follows:
s29, obtaining the channel estimation scheme according to the properties of the MMSE channel estimation schemeWherein
In specific implementation, as a preferred embodiment of the present invention, the step S3 performs downlink data transmission based on the constructed model, and the specific implementation process is as follows:
s31, after the access point AP acquires the CSI from the user, precoding the downlink transmission signal by using a conjugate transpose (CB) method, so that the transmission signal of the mth access point AP is represented as:
wherein x isnRepresenting data signals destined for the nth user, satisfiesThe conditions of (a); p is a radical ofdIndicating the normalized SNR of the downlink transmission data signal; etamnRepresents a power control coefficient;
s32, adjusting power control coefficient etamnEach access point AP is satisfiedThe power constraint of (c), namely:
s33, according to equation (9), after all the access points AP have sent the downlink data signal, the received signal of the nth user is represented as:
wherein, wnCN (0,1) represents additive white Gaussian noise received by the nth user;
s34, based on the low-precision ADC used by the user end, representing the quantized received signal of the nth user as:
wherein, munIndicates the ADC precision of the nth user, which is equal to the quantization bit number sigma used by the nth usernIn connection with, in particularnAnd σnCan refer to the ADC precision alpha of the mth APmAnd the number of quantization bits ρmThe numerical correspondence between the two is not repeated herein.Representing and receiving a signal rnUncorrelated additive quantization noise;
in specific implementation, as a preferred embodiment of the present invention, the step S4 analyzes the downlink reachable rate of the user, and the specific implementation process is as follows:
s41, according to the formula (12), the received signal of the nth user is represented again as follows:
s42, obtaining the expression of the downlink reachable speed of the nth user according to the formula (14) as follows:
s43, deriving the achievable rate closed expression of the nth user according to formula (15):
wherein q isnTo representThe nth column of (1);andrespectively representAnd(N-1). times.N + N1A row; the corresponding elements in the matrices Q, U and V are given as [ Q [ ]]mn=ηmn, δ (m, n) represents a dirac δ function, i.e., when m is equal to n, δ (m, n) is equal to 1; when m ≠ n, δ (m, n) is 0.
In specific implementation, as a preferred embodiment of the present invention, the step S5 constructs an optimization problem that maximizes the total system rate according to the downlink reachable rate, and the specific implementation process is as follows:
s51, aiming at the modern communication system, the total rate of the system and the reachable rate of each user are important performance indexes for evaluating the quality of a communication system. Therefore, the service quality of each user is ensured not to be less than the set minimum achievable rateAnd the transmission power of each AP is not more than the set maximum transmission power pdOn the premise of (2), the optimization problem of maximizing the total rate of the system is summarized into the following form:
s52, optimizing the problem P1The constraint (17b) in (2) is mathematically transformed to obtain:
wherein the constraint (17b) is equivalent to a second order cone in equation (18), i.e. equation (17b) is converted to a convex constraint in equation (18); both equations (17c) and (17d) are convex constraints;
s53, optimization problem P1The objective function in (1) is non-convex, so the optimization problem cannot directly obtain a global optimal solution through a convex optimization solver. Therefore, in order to solve the above optimization problem, in this embodiment, a continuous convex approximation method is adopted to solve the non-convex optimization problem P1Approximate equivalence to optimization problem P2If the optimization problem P can be solved smoothly2Then a solution to the optimization problem P can be obtained1And this suboptimal solution is close to the global optimal solution. In this embodiment, an optimization problem P is given2The concrete form of (A) is as follows:
wherein l ═ l1,…,lN]And t ═ t1,…,tN];
S54, for N ═ 1, …, N, the constraints given are as follows:
wherein the optimization problem P2The objective function (19a) in (1) is convex, but the constraint (20) is non-convex;
s55, converting the constraint (20) into a convex function form, and equating the constraint (20) to N of 1, …, N as follows:
s56, when N is 1, …, N, re-expressing equation (21) in the form:
wherein Q iskAndrespectively represent Q and lnThe iteration values after k SCA treatments, based on which equation (20) has been successfully converted to the constraint conditions of a second order cone by equation (23), and in addition, by observing the optimization problem P2The resulting constraint (19c) is convex;
s57, according to the method of successive convex approximation, approximating the formula (22) by using a first order Taylor expansion in the formula (24);
s58, according to the method of continuous convex approximation, using inequality ln (x) equal to or more than 1-x-1Converting the constraint (19c) into a constraint having the form of a second order cone as follows:
s59, according to the formula (23) and the formula (25), at the k-th iteration, the problem P is optimized2The second order cone programming problem (convexity optimization problem) is expressed as follows:
in specific implementation, as a preferred embodiment of the present invention, the step S6 designs a power control strategy according to the optimization problem, and the specific implementation process is as follows:
s61, initializing and setting the performance parameters: initial iteration number K is 1, maximum iteration number K, tolerance valueMinimum achievable rate of user to be guaranteedInitial value Q1、l1And t1;
S62, solving the optimization problem in the formula (26) through a convex optimization solver, and enabling Q*、l*And t*Is the solution of the iteration;
s64, definition k ═ k +1, Qk=Q*,lk=l*And tk=t*The process goes back to step S62.
Examples
In order to verify the validity of the scheme of the invention, the following simulation experiment is carried out:
setting a scene: all users and access points AP are assumed to be uniformly and randomly distributed in a square area of 1 x 1km, and the distance between the mth access point AP and the nth user is defined as dmn. The large-scale fading coefficients describe path loss and shadow fading and are given by: gamma raymn=-30.18-26log10(dmn)+Fmn(dB) in whichIndicating a shadow fade. Coefficient of performanceAndare independent of each other and are given bysf0.5 and σ sf8. The Rice K-factor describes the ratio between the line-of-sight component and the scatter component, given as
Other parameter values required in the simulation by the inventive solution are set as in table 1.
TABLE 1 parameter settings
As shown in fig. 3, the simulation results of the total rate of the system varying with the number of iterations when using the inventive scheme for power control are given. As can be seen from fig. 2, as the number of iterations increases, the overall rate performance of the system increases continuously, and after the number of iterations exceeds 8, the system gradually approaches a fixed value, which indicates that the method in the present invention will gradually converge after 8 iterations.
As shown in fig. 4, the cumulative distribution function of the total rate of the system when the power control algorithm in the solution of the present invention is used and when the power control algorithm is not used is shown when the number of iterations K is 10. By observation, it can be seen that without the system using a power control algorithm, the total rate of the system is more than 31.8bits/s/Hz with a 95% probability; however, in the case of a system using the scheme of the present invention for power control, the total rate of the system is more than 42.9bits/s/Hz with a 95% probability, which is 1.35 times that of the system without the power control algorithm. The simulation result verifies the effectiveness of the scheme of the invention in improving the overall rate performance of the system.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (7)
1. A power control method for a large-scale MIMO system based on low-precision ADC de-cellular under a Rice channel is characterized by comprising the following steps:
s1, establishing a model for removing the cellular massive MIMO system based on the low-precision ADC under the Rice channel;
s2, performing uplink pilot training based on the constructed model;
s3, carrying out downlink data transmission based on the constructed model;
s4, analyzing the downlink reachable rate of the user;
s5, constructing an optimization problem of maximizing the total rate of the system according to the downlink reachable rate;
and S6, designing a power control strategy according to the optimization problem.
2. The power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel according to claim 1, wherein the model constructed in the step S1 is specifically expressed as follows:
wherein, gmnRepresenting the channel coefficient between the mth access point AP and the nth user; gamma raymnRepresenting the large-scale fading coefficients of the signal,reflecting the influence of path loss and shadow fading on the channel coefficient; the part in parentheses represents the small scale fading coefficient, represented by the line-of-sight componentAnd a scattered component hmnCN (0,1) in whichmn~[-π,π]Represents the angle of arrival and CN (0,1) represents a circularly symmetric complex gaussian variable with a mean of 0 and a variance of 1; in addition, kmnRepresenting the Rice K-factor, representing the ratio between the line-of-sight component and the scatter component, is definedAndwherein beta ismn=γmn/(1+kmn) Thus, the channel coefficients are re-represented asAnd obeyThe statistical distribution of (c).
3. The power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel as claimed in claim 1, wherein the step S2 is implemented as follows:
s21, let the length of the correlation interval be T, and let the time occupied by the pilot training in each correlation interval be taup(τp≥N);
S22, supposing that in the uplink pilot training phase, the nth user can send the uplink pilot sequence to all the APAnd assumes that the pilot sequences sent between different users are relative to each otherOrthogonal, i.e.The uplink pilot information received by the mth access point AP is represented as:
wherein p ispRepresents the normalized signal-to-noise ratio of the uplink pilot signal,representing additive white gaussian noise received by the mth AP;
s23, modeling the quantization process using an additive quantization noise model, and then re-expressing the uplink pilot information received by the mth AP as:
wherein alpha ismIndicating the ADC accuracy at the mth access point AP and the accuracy corresponding to the number of quantization bits ρ used by the mth access point APmRelated, equivalent number of bits ρm1,2,3,4,5, alphamIs precisely given as alpham0.6366,0.8825,0.96546,0.990503, 0.997501; when rhom>At 5 time, αmAnd rhomSatisfyThe approximate relationship of (a) to (b),representing and receiving a signal ym,pExtraneous additive quantization noise;
s25, the mth AP will acquire the channel state information from the nth userAndmultiplying to obtain:
s26, since the access point AP can obtain the line-of-sight component information and the large-scale fading coefficient from different users, in the channel estimation stage, the system may eliminate the influence of the line-of-sight component, that is, equation (5) may be represented again as:
wherein,representing the quantization noise remaining after eliminating the influence of the line-of-sight component;
s28, obtaining g by using the minimum mean square error channel estimation scheme to the formula (6)mnThe channel estimation information of (2) is as follows:
4. The power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel as claimed in claim 1, wherein the step S3 is implemented as follows:
s31, after the access point AP acquires the CSI from the user, precoding the downlink transmission signal in a conjugate transpose manner, so that the transmission signal of the mth access point AP is represented as:
wherein x isnRepresenting a data signal addressed to an nth user, satisfies E { | xn|2Condition of 1; p is a radical ofdIndicating the normalized SNR of the downlink transmission data signal; etamnRepresents a power control coefficient;
s32, adjusting power control coefficient etamnMake each access point AP satisfy E { | sm|2}≤pdWork ofRate constraints, namely:
s33, according to equation (9), after all the access points AP have sent the downlink data signal, the received signal of the nth user is represented as:
wherein, wnCN (0,1) represents additive white Gaussian noise received by the nth user;
s34, based on the low-precision ADC used by the user end, representing the quantized received signal of the nth user as:
wherein, munIndicates the ADC precision of the nth user, which is equal to the quantization bit number sigma used by the nth usernIn connection with this, the present invention is,representing and receiving a signal rnUncorrelated additive quantization noise;
5. the power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel as claimed in claim 1, wherein the step S4 is implemented as follows:
s41, according to the formula (12), the received signal of the nth user is represented again as follows:
s42, obtaining the expression of the downlink reachable speed of the nth user according to the formula (14) as follows:
s43, deriving the achievable rate closed expression of the nth user according to formula (15):
wherein q isnTo representThe nth column of (1);andrespectively representAnd(N-1). times.N + N1A row; the corresponding elements in the matrices Q, U and V are given as [ Q [ ]]mn=ηmn, δ (m, n) represents a dirac δ function, i.e., when m is equal to n, δ (m, n) is equal to 1; when m ≠ n, δ (m, n) is 0.
6. The power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel as claimed in claim 1, wherein the step S5 is implemented as follows:
s51, ensuring the service quality of each user not less than the set minimum reachable rateAnd the transmission power of each AP is not more than the set maximum transmission power pdOn the premise of (2), the optimization problem of maximizing the total rate of the system is summarized into the following form:
s52, optimizing the problem P1The constraint (17b) in (2) is mathematically transformed to obtain:
wherein the constraint (17b) is equivalent to a second order cone in equation (18), i.e. equation (17b) is converted to a convex constraint in equation (18); equations (17c) and (17d) are also both convex constraints;
s53, adopting a continuous convex approximation method to optimize the non-convex problem P1Approximate equivalence to optimization problem P2Form of (2), optimization problem P2The concrete form of (A) is as follows:
wherein l ═ l1,…,lN]And t ═ t1,…,tN];
S54, for N ═ 1, …, N, the constraints given are as follows:
wherein the optimization problem P2The objective function (19a) in (1) is convex, but the constraint (20) is non-convex;
s55, converting the constraint (20) into a convex function form, and equating the constraint (20) to N of 1, …, N as follows:
s56, when N is 1, …, N, re-expressing equation (21) in the form:
wherein Q iskAndrespectively represent Q and lnThe iteration values after k SCA treatments, based on which equation (20) has been successfully converted to the constraint conditions of a second order cone by equation (23), and in addition, by observing the optimization problem P2The resulting constraint (19c) is convex;
s57, according to the method of successive convex approximation, approximating the formula (22) by using a first order Taylor expansion in the formula (24);
s58, according to the method of continuous convex approximation, using inequality ln (x) equal to or more than 1-x-1Converting the constraint (19c) into a constraint having the form of a second order cone as follows:
s59, according to the formula (23) and the formula (25), at the k-th iteration, the problem P is optimized2A second order cone programming problem expressed in the form:
7. the power control method for the massive MIMO system based on the low-precision ADC de-cellular under rice channel as claimed in claim 1, wherein the step S6 is implemented as follows:
s61, initializing and setting the performance parameters: the initial iteration number K is 1, the maximum iteration number K, the tolerance value theta is 0.01, and the minimum reachable rate of the userInitial value Q1、l1And t1;
S62, solving the optimization problem in the formula (26) through a convex optimization solver, and enabling Q*、l*And t*Is the solution of the iteration;
s64, definition k ═ k +1, Qk=Q*,lk=l*And tk=t*The process goes back to step S62.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110626350.9A CN113364501B (en) | 2021-06-04 | 2021-06-04 | Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110626350.9A CN113364501B (en) | 2021-06-04 | 2021-06-04 | Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113364501A true CN113364501A (en) | 2021-09-07 |
CN113364501B CN113364501B (en) | 2022-07-05 |
Family
ID=77532336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110626350.9A Active CN113364501B (en) | 2021-06-04 | 2021-06-04 | Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113364501B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113890578A (en) * | 2021-09-29 | 2022-01-04 | 浙大城市学院 | Cellular-free system power design method based on RIS and related channels |
CN114665930A (en) * | 2022-03-16 | 2022-06-24 | 南京邮电大学 | Downlink blind channel estimation method of large-scale de-cellular MIMO system |
CN114980332A (en) * | 2022-05-17 | 2022-08-30 | 清华大学 | Downlink power distribution method and device for large-scale cellular MIMO (multiple input multiple output) system |
CN115021846A (en) * | 2022-05-23 | 2022-09-06 | 浙江师范大学 | Balanced optimization method for spectrum efficiency and energy efficiency of large-scale cellular MIMO downlink |
CN116707687A (en) * | 2023-06-27 | 2023-09-05 | 南京盛航海运股份有限公司 | Channel prediction method for de-cellular large-scale MIMO system |
CN116865798A (en) * | 2023-07-06 | 2023-10-10 | 河北大学 | Intelligent super-surface phase shift method for high-speed railway honeycomb removing large-scale MIMO system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105049383A (en) * | 2015-07-01 | 2015-11-11 | 东南大学 | FDD large-scale MIMO system downlink training sequence design method |
CN110166088A (en) * | 2019-05-15 | 2019-08-23 | 南京邮电大学 | The power control algorithm without cell mimo system of customer-centric |
CN111405596A (en) * | 2020-03-24 | 2020-07-10 | 西安电子科技大学 | Resource optimization method for large-scale antenna wireless energy-carrying communication system under Rice channel |
-
2021
- 2021-06-04 CN CN202110626350.9A patent/CN113364501B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105049383A (en) * | 2015-07-01 | 2015-11-11 | 东南大学 | FDD large-scale MIMO system downlink training sequence design method |
CN110166088A (en) * | 2019-05-15 | 2019-08-23 | 南京邮电大学 | The power control algorithm without cell mimo system of customer-centric |
CN111405596A (en) * | 2020-03-24 | 2020-07-10 | 西安电子科技大学 | Resource optimization method for large-scale antenna wireless energy-carrying communication system under Rice channel |
Non-Patent Citations (5)
Title |
---|
SINIAN JIN 等: "On the Energy Efficiency of Multi-Cell Massive MIMO With Beamforming Training", 《IEEE ACCESS》, no. 8, 27 April 2020 (2020-04-27), pages 80739 - 80754, XP011787810, DOI: 10.1109/ACCESS.2020.2990456 * |
SINIAN JIN 等: "Spectral and Energy Efficiency in Cell-Free Massive MIMO Systems Over Correlated Rician Fading", 《IEEE SYSTEMS JOURNAL》, vol. 15, no. 2, 26 May 2020 (2020-05-26), pages 2822 - 2833, XP011859078, DOI: 10.1109/JSYST.2020.2993048 * |
章嘉懿: "去蜂窝大规模MIMO系统研究进展与发展趋势", 《重庆邮电大学学报(自然科学版)》 * |
章嘉懿: "去蜂窝大规模MIMO系统研究进展与发展趋势", 《重庆邮电大学学报(自然科学版)》, no. 03, 15 June 2019 (2019-06-15) * |
谢斌等: "莱斯信道下大规模MIMO系统混合模数转换器的随机向量量化方案", 《高技术通讯》, no. 04, 15 April 2017 (2017-04-15), pages 3 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113890578A (en) * | 2021-09-29 | 2022-01-04 | 浙大城市学院 | Cellular-free system power design method based on RIS and related channels |
CN114665930A (en) * | 2022-03-16 | 2022-06-24 | 南京邮电大学 | Downlink blind channel estimation method of large-scale de-cellular MIMO system |
CN114980332A (en) * | 2022-05-17 | 2022-08-30 | 清华大学 | Downlink power distribution method and device for large-scale cellular MIMO (multiple input multiple output) system |
CN115021846A (en) * | 2022-05-23 | 2022-09-06 | 浙江师范大学 | Balanced optimization method for spectrum efficiency and energy efficiency of large-scale cellular MIMO downlink |
CN116707687A (en) * | 2023-06-27 | 2023-09-05 | 南京盛航海运股份有限公司 | Channel prediction method for de-cellular large-scale MIMO system |
CN116707687B (en) * | 2023-06-27 | 2024-05-10 | 南京盛航海运股份有限公司 | Channel prediction method for de-cellular large-scale MIMO system |
CN116865798A (en) * | 2023-07-06 | 2023-10-10 | 河北大学 | Intelligent super-surface phase shift method for high-speed railway honeycomb removing large-scale MIMO system |
CN116865798B (en) * | 2023-07-06 | 2024-01-05 | 河北大学 | Intelligent super-surface phase shift method for high-speed railway honeycomb removing large-scale MIMO system |
Also Published As
Publication number | Publication date |
---|---|
CN113364501B (en) | 2022-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113364501B (en) | Power control method based on low-precision ADC (analog to digital converter) de-cellular large-scale MIMO (multiple input multiple output) system under Rice channel | |
CN106253956B (en) | Codebook-based modulus mixing method for precoding | |
CN113179109B (en) | Honeycomb-removing large-scale MIMO uplink spectrum efficiency optimization method | |
CN111835406B (en) | Robust precoding method suitable for energy efficiency and spectral efficiency balance of multi-beam satellite communication | |
CN108365873B (en) | Large-scale MIMO self-adaptive transmission method adopting low-precision ADC millimeter waves | |
Liu et al. | Two-timescale hybrid compression and forward for massive MIMO aided C-RAN | |
CN110365388B (en) | Low-complexity millimeter wave multicast beam forming method | |
CN110166088B (en) | Power control algorithm of user-centered cell-free MIMO system | |
CN111970033B (en) | Large-scale MIMO multicast power distribution method based on energy efficiency and spectrum efficiency joint optimization | |
CN113438002B (en) | LSTM-based analog beam switching method, device, equipment and medium | |
CN107171708A (en) | A kind of channel tracking of extensive mimo system is with mixing method for precoding | |
CN112260737A (en) | Multi-beam satellite communication robust precoding method with total energy efficiency and minimum energy efficiency balanced | |
Kaushik et al. | Energy efficiency maximization of millimeter wave hybrid MIMO systems with low resolution DACs | |
CN115941010B (en) | IRS auxiliary honeycomb removing large-scale MIMO system beam forming method based on branch definition | |
Ding et al. | Performance analysis of mixed-ADC massive MIMO systems over spatially correlated channels | |
CN111510188B (en) | Beam searching method and device | |
CN108173575B (en) | Design method of multi-input multi-output relay antenna | |
CN105812032A (en) | Channel estimation method based on beam block structure compressed sensing | |
CN112261728A (en) | Beam selection matrix design method based on lens array | |
CN112312569A (en) | Lens array-based precoding and beam selection matrix joint design method | |
CN108923831B (en) | Method and device for precoding transmission signals | |
CN106533524A (en) | Forming method for beam with maximum energy efficiency in distributed antenna system | |
Choi et al. | User scheduling for millimeter wave MIMO communications with low-resolution ADCs | |
Kaushik et al. | Energy efficiency maximization in millimeter wave hybrid MIMO systems for 5G and beyond | |
Xiao et al. | Proportionally fair robust beamforming for multicast multibeam satellite systems |
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 |