CN114302487A - Energy efficiency optimization method, device and equipment based on adaptive particle swarm power distribution - Google Patents

Energy efficiency optimization method, device and equipment based on adaptive particle swarm power distribution Download PDF

Info

Publication number
CN114302487A
CN114302487A CN202111628181.9A CN202111628181A CN114302487A CN 114302487 A CN114302487 A CN 114302487A CN 202111628181 A CN202111628181 A CN 202111628181A CN 114302487 A CN114302487 A CN 114302487A
Authority
CN
China
Prior art keywords
energy efficiency
cell
particle
channel
user
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111628181.9A
Other languages
Chinese (zh)
Other versions
CN114302487B (en
Inventor
邓宏贵
张捷
封雨鑫
刘洪梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN202111628181.9A priority Critical patent/CN114302487B/en
Publication of CN114302487A publication Critical patent/CN114302487A/en
Application granted granted Critical
Publication of CN114302487B publication Critical patent/CN114302487B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses an energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution, wherein the method comprises the following steps: initializing the positions of all base stations and the positions of all users in a cell; calculating a large-scale fading factor beta and a covariance matrix R of a channel according to the position of each base station and the position of each user; receiving an uplink pilot signal, and estimating a channel by adopting an MMSE (minimum mean square error) method by combining a large-scale fading factor beta and a covariance matrix R of the channel; deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system; and for the energy efficiency optimization model, according to an optimization objective function, carrying out data power distribution on a user by using a self-adaptive particle swarm algorithm on the basis of fixed pilot frequency power. The invention can effectively improve the energy efficiency of the system and meet the requirement of green communication.

Description

Energy efficiency optimization method, device and equipment based on adaptive particle swarm power distribution
Technical Field
The invention belongs to the technical field of wireless communication of large-scale MIMO systems, and particularly relates to an energy efficiency optimization method, device and equipment based on adaptive particle swarm power distribution.
Background
In a large-scale MIMO system, hundreds of thousands of antennas are configured at a base station, and spatial diversity gain and spatial multiplexing gain are obtained by using multipath scattering, so that the system can greatly improve the system frequency efficiency, energy efficiency and reliability of a wireless link under the condition of an original bandwidth, and the system becomes one of 5G key technologies. In a large-scale MIMO system, a TDD mode is adopted in order to utilize reciprocity of uplink and downlink channels within the same coherence time. The position relationship between the user and the base station and the surrounding environment are different, so that the channel gains between the user and the base station are different, and the fading degrees of signals in the channel transmission process are also different. When uplink data transmission is performed, due to the non-orthogonality of channels between users, signals transmitted between users interfere with each other, so that signals received by a base station end include useful signals and interference signals, and the performance of the whole communication system is affected. Power control has always been one of the hot problems in communication system research. In the energy efficiency research of a large-scale MIMO system, the total power consumption of the system can be reduced under the condition of ensuring the data transmission rate by controlling the power of the pilot frequency and the data signal, the energy efficiency of the system is improved, the performance of the system is further improved, and the requirement of green communication is met.
At present, the power distribution problem has been widely studied, and the specific scheme is as follows:
the first and the second power distribution algorithms are that the total uplink power is distributed to each user in the cell averagely without considering the channel gain. This approach can make the communication quality of the user with good channel gain good, and the communication quality of the user with poor channel gain not so poor, but this algorithm does not consider the user minimum rate, and is not optimal for energy efficiency and frequency efficiency, because the number of antennas configured by the base station cannot approach infinity in reality.
And secondly, a water injection power algorithm from the angle of channel gain, wherein the water injection power algorithm is used for adaptively distributing the transmission power according to the channel state, and usually, the transmission power is distributed more at the moment when the channel state is good, and the transmission power is distributed less at the moment when the channel state is poor, so that the maximum transmission rate is achieved. The power allocation algorithm is simple in implementation, but fairness of user quality is not considered, so that users with good channels allocate more power and have good communication quality, users with poor channels allocate less power and have poor communication quality, and strong interference of other users is received.
According to a power distribution algorithm of a convex optimization theory such as a fractional programming theory, Lagrange multiplication and Dinkelbach algorithm and the like, the energy efficiency optimization target of the large-scale MIMO system is non-convex, generally, the non-convex target function is reconstructed to form a convex target function through relaxation and reforming operations, then, power is distributed through a convex optimization method such as an interior point method and the like, so that the obtained local optimal power distribution vector is the optimal power distribution vector, and the system has good performance for energy efficiency and frequency efficiency.
Disclosure of Invention
Aiming at various problems which are not fully considered in the existing power distribution scheme, the invention aims to provide an energy efficiency optimization method, device and equipment based on adaptive particle swarm power distribution.
In order to achieve the purpose, the invention provides the following technical scheme:
an energy efficiency optimization method based on adaptive particle swarm power distribution comprises the following steps:
s1, initializing each base station position and each user position in the cell;
s2, calculating a large-scale fading factor beta and a covariance matrix R of a channel according to the positions of all base stations and the positions of all users;
s3, receiving an uplink pilot signal, and estimating a channel by adopting an MMSE (minimum mean square error) method by combining a large-scale fading factor beta and a covariance matrix R of the channel;
s4, deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system;
and S5, according to the optimization objective function, the energy efficiency optimization model carries out data power distribution of users by using a self-adaptive particle swarm algorithm on the basis of fixed pilot frequency power.
Preferably, in step S2, the large scale fading factor
Figure BDA0003439331980000031
Expressed as:
Figure BDA0003439331980000032
wherein the content of the first and second substances,
Figure BDA0003439331980000033
representing shadow fading, logarithm thereof
Figure BDA0003439331980000034
Obeying a Gaussian distribution
Figure BDA0003439331980000035
Figure BDA0003439331980000041
Figure BDA0003439331980000042
Denotes the distance r from the ith user in cell l to the base station at the center of cell j0Represents the radius of the cell; α represents a path loss coefficient during signal transmission.
The covariance matrix R of the channel is obtained according to the approximate solution of a Gaussian local scattering model, and the following conditions are met:
Figure BDA0003439331980000043
where β represents a large scale fading factor, dHM is the number of antennas configured for the base station,
Figure BDA0003439331980000044
is the angle of arrival of the user.
Preferably, step S3 specifically includes:
in the pilot transmission phase, the signal received by the base station in cell j is transmitted
Figure BDA0003439331980000045
Expressed as:
Figure BDA0003439331980000046
wherein p isjkRepresenting the power at which the kth user in cell j transmits a pilot,
Figure BDA0003439331980000047
a pilot sequence transmitted for the kth user in cell j;
Figure BDA0003439331980000048
representing additive white Gaussian noise at the receiving end of the base stationFrom independent and identically distributed CN to (0, sigma)p);
Multiplication of both sides of the above formula
Figure BDA0003439331980000049
Obtaining:
Figure BDA00034393319800000410
wherein, since the pilots are orthogonal,
Figure BDA00034393319800000411
a value of 0;
using a minimum mean square error estimation method based on
Figure BDA00034393319800000412
Performing channel estimation, wherein the channel estimation value between the kth user in the cell l and the base station j can be expressed as:
Figure BDA0003439331980000051
Figure BDA0003439331980000052
therefore, in the uplink data transmission phase, the data received by the base station j can be represented as:
Figure BDA0003439331980000053
maximum ratio combining reception is adopted for the received signals, namely the received vector is as follows:
Figure BDA0003439331980000054
therefore, the transmitted signal of user k in cell j is represented by:
Figure BDA0003439331980000055
wherein the first term is a useful signal; the second term is intra-cell interference; the remaining part is inter-cell interference and other uncorrelated noise, sjkIs the transmitted signal of user k in cell j.
Preferably, in step S4, the signal to interference plus noise ratio SINR of the system
Figure BDA0003439331980000056
The lower bound for knowing the frequency efficiency of the system by the shannon capacity theorem is expressed as:
Figure BDA0003439331980000057
the optimization problem is represented as:
Figure BDA0003439331980000061
Figure BDA0003439331980000062
Figure BDA0003439331980000063
wherein the content of the first and second substances,
Figure BDA0003439331980000064
represents the signal-to-interference-and-noise ratio, p, of the kth user signal in cell j received by the base station in cell jliDenotes the data power, V, allocated to user i in cell jjkIndicating that the combined vector is received,
Figure BDA0003439331980000065
representing the estimated value of the channel between the base station of cell j and the kth user in the cell,
Figure BDA0003439331980000066
an estimation error matrix representing a channel matrix,
Figure BDA0003439331980000067
which is indicative of the power of the noise,
Figure BDA0003439331980000068
is MjIdentity matrix of order, EEULRepresenting system energy efficiency, τ, in a massive MIMO uplink transmission processuIndicating the length of the coherent block, τ, for transmitting uplink datacDenotes the total length of the coherent block, M being the number of antennas with which the base station is equipped, prPower consumed by the radio link of the antenna for uplink data transmission, psStatic electrical power consumption, p, consumed during data transmission for large-scale MIMO systemsmaxIs the maximum electrical power during uplink transmission for all users in each cell. RminIs a value set under consideration of the minimum rate constraint.
Preferably, step S5 specifically includes:
s51, initializing the parameters of the particles by using the rand function and the upper and lower limits of the set user data power;
s52, taking the energy efficiency function as a fitness function, calculating the fitness of the initialized particles, and initializing a local optimal solution and a global optimal solution;
s53, judging the iteration termination condition, terminating the iteration when the iteration reaches the initial set maximum iteration number or the fitness value tends to be stable, otherwise, continuing the following steps;
s54, updating the speed and the position of the particle, simultaneously performing boundary processing on the particle beyond the boundary range, and recalculating the fitness value of the updated particle;
and S55, storing the result of each iteration, comparing the results after each iteration, and updating the local optimal solution and the global optimal solution.
And S56, when the iteration termination condition is reached, obtaining the global optimal value, namely the required optimal user data power allocation vector.
Preferably, in step S52, the fitness function is set as the energy efficiency expression:
Figure BDA0003439331980000071
and substituting the position parameter of each particle as an input parameter into the formula to solve a corresponding fitness value, finishing the first iteration process, and setting the initialized particle parameters as a local optimal value and a global optimal value.
Preferably, in step S54, the velocity and position of the particle are updated and boundary processing is performed as follows:
the update formula of the particle speed and the position is as follows:
Vj(t+1)=Vj(t)+c1*rand*(Gbestj-popj(t))+c2*rand*(Zbest-popj(t))
popj(t+1)=popj(t)+ω*Vj(t+1)
Figure BDA0003439331980000072
in the formula, Vj(t +1) represents the search speed, V, of the jth particle at the (t +1) th iterationj(t) represents the search speed of the jth particle in the tth iteration, c1 is an individual learning factor, the iteration representing the speed is only related to the historical position of the particle, c2 is a social learning factor, the relation between the iteration representing the speed and the historical position of the whole particle swarm is represented, omega represents a weight factor, when omega is larger, the local minimum value can be skipped, the global search is facilitated, and when the inertia factor omega is smaller, the local trend can be accurately searched. Thus, when the number of iterations is small, the value of ω is large, and when the number of iterations is largeWhen large, the value of ω is small. GbestjZbest represents the optimal particle position, pop, for the current position of the jth particle for the current iteration number of the entire particle swarmj(t) represents a position parameter representing the jth particle at the tth iteration.
The boundary processing method of the particle speed and the particle position comprises the following steps:
Figure BDA0003439331980000081
Figure BDA0003439331980000082
v (i, j) represents the velocity of the ith particle in the j-dimension direction, the transition speed of the particle in the solution space can be limited by processing the boundary of the particle velocity, and the solution space can be searched sufficiently. pop (i, j) indicates the displacement size of the ith particle in the j-dimension direction, and the position is subjected to boundary processing, so that the particle can be strictly limited to move in the solution space.
The embodiment of the present invention further provides an energy efficiency optimization device based on adaptive particle swarm power distribution, which includes:
the initialization unit is used for initializing the positions of all base stations and the positions of all users in a cell;
a calculating unit, configured to calculate a large-scale fading factor β and a covariance matrix R of a channel according to the positions of the base stations and the positions of the users;
the channel estimation unit is used for receiving the uplink pilot signal, and simultaneously estimating the channel by adopting an MMSE (minimum mean square error) method in combination with the large-scale fading factor beta and the covariance matrix R of the channel;
the model establishing unit is used for deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to the Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system;
and the power distribution unit is used for carrying out data power distribution on the user on the basis of fixed pilot frequency power by utilizing a self-adaptive particle swarm algorithm according to the optimization objective function for the energy efficiency optimization model.
The embodiment of the invention also provides energy efficiency optimization equipment based on adaptive particle swarm power distribution, which comprises a memory and a processor, wherein a computer program is stored in the memory, and the computer program can be executed by the processor so as to realize the energy efficiency optimization method based on adaptive particle swarm power distribution.
In the invention, the communication quality fairness of each user can be fully considered under the condition of low calculation complexity, the energy efficiency of the system is effectively improved, and the requirement of green communication is met.
Drawings
Fig. 1 is a schematic flowchart of a work flow of an energy efficiency optimization method based on adaptive particle swarm power allocation according to a first embodiment of the present invention.
Fig. 2 is a schematic block diagram of a flow of an adaptive particle swarm algorithm provided by an embodiment of the present invention.
Fig. 3 is a schematic diagram illustrating comparison of average frequency efficiency performance of cells based on adaptive particle swarm optimization and other algorithms provided by the present invention.
Fig. 4 is a schematic diagram illustrating comparison of average energy efficiency performance of a cell based on an adaptive particle swarm algorithm and other algorithms according to the present invention.
Fig. 5 is a schematic structural diagram of an energy efficiency optimization device based on adaptive particle swarm power allocation according to a second embodiment of the present invention.
Detailed Description
The invention will be further elucidated with reference to the embodiments and the accompanying drawings.
Referring to fig. 1, a first embodiment of the present invention provides an energy efficiency optimization method based on adaptive particle swarm power allocation, which includes the following steps:
s1, the location of each base station and the location of each user in the cell are initialized.
In one implementation manner, the present embodiment provides 16 cells in total, and the number of the cells is 4 × 4The arrangement is made such that the system employs orthogonal pilots for users within a cell in order to overcome inter-cell interference. Meanwhile, in order to avoid loss of generality, the pilot multiplexing factor of the system is set to be 4. Thus, the pilot sequence employed by the present system may be represented as
Figure BDA0003439331980000101
Figure BDA0003439331980000102
In each cell, users are set in a random distribution mode, the minimum distance from the users to the base station is constrained, and the topological distance of the surrounding base stations is calculated for each user.
S2, calculating a large-scale fading factor β and a covariance matrix R of the channel according to the positions of the base stations and the positions of the users.
In this example, consider a multi-cell, multi-user massive MIMO system, which consists of L square cells. In each cell, a base station equipped with M antennas is set up in the center of the cell, and serves K single-antenna users randomly distributed in the cell (where K is<<M). Without loss of generality, the channel gain from the ith user to cell base station j in cell l
Figure BDA0003439331980000103
Can be expressed as:
Figure BDA0003439331980000104
wherein the content of the first and second substances,
Figure BDA0003439331980000111
representing small scale fading, obeying a circularly symmetric complex Gaussian distribution, i.e. CN — (0, I)M)。
Wherein the content of the first and second substances,
Figure BDA0003439331980000112
is a large scale fading factor, generally expressed as
Figure BDA0003439331980000113
Herein, the
Figure BDA0003439331980000114
Representing shadow fading, logarithm thereof
Figure BDA0003439331980000115
Obeying a Gaussian distribution
Figure BDA0003439331980000116
Figure BDA0003439331980000117
Wherein
Figure BDA0003439331980000118
The distance from the ith user in the cell l to the base station at the center of the cell j is shown, r0 shows the radius of the cell, and alpha shows the path loss coefficient in the signal transmission process; in particular, it is possible to use, for example,
Figure BDA0003439331980000119
slowly varying and easily tracking over the coherence time of several channels, so that within one coherence time
Figure BDA00034393319800001110
Considered as a constant.
In this embodiment, according to the set base station and the user position, the arrival angle between the user and the base station is calculated, and the covariance matrix R of the channel is solved according to the approximate solution of the gaussian local scattering model:
Figure BDA00034393319800001111
where β represents a large scale fading factor, dHFor spatial distance of antenna (generally set as half wavelength), M is allocated to base stationThe number of the antennas to be placed is,
Figure BDA00034393319800001112
is the angle of arrival of the user.
And S3, receiving the uplink pilot signal, and estimating the channel by adopting an MMSE (minimum mean square error) method by combining the large-scale fading factor beta and the covariance matrix R of the channel.
Wherein, in the pilot frequency transmitting phase, the signal received by the base station in the cell j
Figure BDA00034393319800001113
Can be expressed as:
Figure BDA00034393319800001114
where p isjkIndicating the power at which the kth user in cell j transmits a pilot,
Figure BDA00034393319800001115
a pilot sequence transmitted for the kth user in cell j.
Figure BDA0003439331980000121
It is the additive white Gaussian noise at the receiving end of the base station, which is subject to independent same distribution CN (0-sigma)p)。
Multiplying both sides of the equation (4) simultaneously
Figure BDA0003439331980000122
It is possible to obtain:
Figure BDA0003439331980000123
due to the orthogonality of the pilot frequency, the third term on the right side of the above equation (5) is 0, so that the inter-cell interference can be removed, which is beneficial to channel estimation.
Through the above operations, after the base station in the cell receives the pilot signal sent by the user, the base stationCan be based on
Figure BDA0003439331980000124
For channel estimation, the present embodiment adopts a Minimum Mean Square Error (MMSE) estimation method. Thus, the channel estimate between the kth user in cell i and base station j can be expressed as:
Figure BDA0003439331980000125
Figure BDA0003439331980000126
therefore, in the uplink data transmission phase, the data received by the base station j can be represented as:
Figure BDA0003439331980000127
in this embodiment, maximum ratio combining reception is applied to the received signals, that is, the received vector is:
Figure BDA0003439331980000128
therefore, the transmitted signal for user k in cell j can be represented by the following sub-expression:
Figure BDA0003439331980000129
wherein the first term is a useful signal; the second term is intra-cell interference; the remaining portion is inter-cell interference and other uncorrelated noise.
And S4, deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to the Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system.
And deducing SINR of the system and expressions of frequency efficiency SE and energy efficiency EE according to the known process. Wherein:
Figure BDA0003439331980000131
the lower bound of the frequency efficiency of the system known by the shannon capacity theorem can be expressed as:
Figure BDA0003439331980000132
the purpose of this embodiment is to try to improve the average energy efficiency in the system, and to realize the development of green communication once. The optimization problem can therefore be expressed as:
Figure BDA0003439331980000133
Figure BDA0003439331980000134
Figure BDA0003439331980000135
wherein the content of the first and second substances,
Figure BDA0003439331980000136
represents the signal-to-interference-and-noise ratio, p, of the kth user signal in cell j received by the base station in cell jliDenotes the data power, V, allocated to user i in cell jjkIndicating that the combined vector is received,
Figure BDA0003439331980000137
representing the estimated value of the channel between the base station of cell j and the kth user in the cell,
Figure BDA0003439331980000138
an estimation error matrix representing a channel matrix,
Figure BDA0003439331980000139
which is indicative of the power of the noise,
Figure BDA00034393319800001310
is MjIdentity matrix of order, EEULRepresenting system energy efficiency, τ, in a massive MIMO uplink transmission processuIndicating the length of the coherent block, τ, for transmitting uplink datacDenotes the total length of the coherent block, M being the number of antennas with which the base station is equipped, prPower consumed by the radio link of the antenna for uplink data transmission, psStatic electrical power consumption, p, consumed during data transmission for large-scale MIMO systemsmaxIs the maximum electrical power during uplink transmission for all users in each cell. RminIs a value set under consideration of the minimum rate constraint.
And S5, according to the optimization objective function, the energy efficiency optimization model carries out data power distribution of users by using a self-adaptive particle swarm algorithm on the basis of fixed pilot frequency power.
Preferably, step S5 specifically includes:
s51, initializing the particles by using the rand function and the upper and lower limits of the set user data power,
Figure BDA0003439331980000141
wherein is at
Figure BDA0003439331980000142
The upper mark is iteration times, and the lower mark is particle number; p is a radical of12The first numerical subscript represents the cell number and the second numerical subscript represents the user number in the cell.
S52, taking the energy efficiency function EE as a fitness function, calculating the fitness of all initialized particles, and initializing local optimal
Figure BDA0003439331980000143
And global optimum
Figure BDA0003439331980000144
S53, updating the speed and position of the particle and performing boundary processing:
Figure BDA0003439331980000145
Figure BDA0003439331980000146
calculating the updated fitness value EE corresponding to the particles and updating the local optimal solution
Figure BDA0003439331980000147
And global optimal solution
Figure BDA0003439331980000148
And S54, processing boundary conditions, wherein the power constraint condition is required to be used as the boundary condition of the particle position in the iteration process, and the motion of the particle is constrained within a certain range.
The updating formula of the particle speed and the position is as follows:
Vj(t+1)=Vj(t)+c1*rand*(Gbestj-popj(t))+c2*rand*(Zbest-popj(t)) (15)
Figure BDA0003439331980000151
in the formula, Vj(t +1) represents the velocity value of the jth particle at the (t +1) th iteration, Vj(t) represents the speed value of the jth particle in the tth iteration, c1 is an individual learning factor, the iteration representing the speed is only related to the historical position of the particle, c2 is a social learning factor, the iteration representing the speed is related to the historical position of the whole particle swarm, and omega represents a weight factorLocal search is facilitated, and when the inertia factor omega is small, accurate search tends to be conducted on local parts. Therefore, when the number of iterations is small, the value of ω is large, and when the number of iterations is large, the value of ω is small. GbestjThe optimal position of the jth particle's current position, Zbest is the optimal particle position so far, pop, of the entire particle populationj(t) represents a position parameter representing the jth particle at the tth iteration.
The boundary processing method of the particle speed and the particle position comprises the following steps:
Figure BDA0003439331980000152
Figure BDA0003439331980000153
v (i, j) represents the velocity of the ith particle in the j-dimension direction, and boundary processing is performed on the particle velocity, so that the search velocity of the particle in the solution space can be limited, and the solution space can be sufficiently searched. pop (i, j) indicates the displacement size of the ith particle in the j-dimension direction, and the particle motion can be strictly limited in the solution space by performing boundary processing on the position.
In step S54, after the speed and position of each particle are updated, the fitness value of the newly generated particle group is calculated, and according to the comparison between the fitness value of the current particle and the fitness value of the previous particle, the optimal position of each particle is selected as the local optimal position, and the optimal position of the whole particle group is selected as the global optimal position. While the boundary processing for particle velocity and position can be more efficient and reliable.
S55, in order to avoid the particle swarm optimization falling into a local extremum, the embodiment combines the mutation operation of the genetic algorithm, and processes the positions of the particles with a variation probability of twenty percent:
k=ceil(LK*rand) (19)
pop(i,k)=rand*(popmax-popmin)+popmin (20)
and S56, when the iteration termination condition is reached, the global optimal value is the optimal user data power allocation vector to be selected.
Fig. 3 and 4 are graphs of performance simulation results of embodiments of the present invention. It can be seen from the figure that, under the same condition, compared with the conventional power allocation algorithm, the method of the present invention has a certain improvement in frequency efficiency and a great improvement in energy efficiency, and as the number of antennas increases, the performance of several methods finally gradually decreases. Therefore, compared with the traditional allocation algorithm, the method has obvious performance improvement.
Referring to fig. 5, a second embodiment of the present invention further provides an energy efficiency optimization apparatus based on adaptive particle swarm power allocation, including:
an initializing unit 210, configured to initialize positions of each base station and each user in a cell;
a calculating unit 220, configured to calculate a large-scale fading factor β and a covariance matrix R of a channel according to the positions of the base stations and the positions of the users;
a channel estimation unit 230, configured to receive an uplink pilot signal, and estimate a channel by using an MMSE method in combination with a large-scale fading factor β and a covariance matrix R of the channel;
the model establishing unit 240 is configured to derive a signal to interference plus noise ratio SINR expression of the system according to the estimated channel, calculate frequency efficiency according to shannon capacity theorem, and establish an energy efficiency optimization model of the system in combination with a power consumption model of the system;
and the power distribution unit 250 is configured to distribute the data power of the user to the energy efficiency optimization model according to the optimization objective function by using a self-adaptive particle swarm algorithm on the basis of the fixed pilot power.
The third embodiment of the present invention further provides an energy efficiency optimization device based on adaptive particle swarm power allocation, which includes a memory and a processor, where a computer program is stored in the memory, and the computer program can be executed by the processor, so as to implement the energy efficiency optimization method based on adaptive particle swarm power allocation as described above.
The specific embodiments described herein are merely illustrative of the spirit of the invention. Various modifications or additions may be made to the described embodiments or alternatives may be employed by those skilled in the art without departing from the spirit or scope of the invention as defined in the appended claims.

Claims (9)

1. An energy efficiency optimization method based on adaptive particle swarm power distribution is characterized by comprising the following steps:
s1, initializing each base station position and each user position in the cell;
s2, calculating a large-scale fading factor beta and a covariance matrix R of a channel according to the positions of all base stations and the positions of all users;
s3, receiving an uplink pilot signal, and estimating a channel by adopting an MMSE (minimum mean square error) method by combining a large-scale fading factor beta and a covariance matrix R of the channel;
s4, deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system;
and S5, according to the optimization objective function, the energy efficiency optimization model carries out data power distribution of users by using a self-adaptive particle swarm algorithm on the basis of fixed pilot frequency power.
2. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 1, wherein in step S2, the large-scale fading factor
Figure FDA0003439331970000011
Expressed as:
Figure FDA0003439331970000012
wherein the content of the first and second substances,
Figure FDA0003439331970000013
representing shadow fading, logarithm thereof
Figure FDA0003439331970000014
Obeying a Gaussian distribution
Figure FDA0003439331970000015
Figure FDA0003439331970000016
Figure FDA0003439331970000017
Represents the distance from the ith user in cell l to the base station at the center of cell j, and r0 represents the radius of the cell; α represents the number of path losses during signal transmission.
The covariance matrix R of the channel is obtained according to the approximate solution of a Gaussian local scattering model, and the following conditions are met:
Figure FDA0003439331970000021
Figure FDA0003439331970000022
where β represents a large scale fading factor, dHM is the number of antennas configured for the base station,
Figure FDA0003439331970000023
is the angle of arrival of the user.
3. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 2, wherein step S3 specifically includes:
in the pilot transmission phase, the signal received by the base station in cell j is transmitted
Figure FDA0003439331970000024
Expressed as:
Figure FDA0003439331970000025
wherein p isjkRepresenting the power at which the kth user in cell j transmits a pilot,
Figure FDA0003439331970000026
a pilot sequence transmitted for the kth user in cell j;
Figure FDA0003439331970000027
representing additive white Gaussian noise at the receiving end of the base station, which obeys independent same distribution CN to (0, sigma)p);
Multiplication of both sides of the above formula
Figure FDA0003439331970000028
Obtaining:
Figure FDA0003439331970000029
wherein, due to the adoption of the orthogonal pilot frequency,
Figure FDA00034393319700000210
a value of 0;
using a minimum mean square error estimation method based on
Figure FDA00034393319700000211
Performing channel estimation, wherein the cellThe channel estimate between the kth user in l and base station j can be expressed as:
Figure FDA00034393319700000212
Figure FDA0003439331970000031
wherein P isliRepresenting the set of users transmitting the same pilot as user i in cell l.
Therefore, in the uplink data transmission phase, the data received by the base station j is represented as:
Figure FDA0003439331970000032
maximum ratio combining reception is adopted for the received signals, namely the received vector is as follows:
Figure FDA0003439331970000033
therefore, the transmitted signal of user k in cell j is represented by:
Figure FDA0003439331970000034
wherein the first term is a useful signal; the second term is intra-cell interference; the remaining part is inter-cell interference and other non-correlated noise; sjkRepresenting the transmitted signal of user k in cell j.
4. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 3, wherein in step S4,
Figure FDA0003439331970000035
the lower bound of the frequency efficiency of the system can be known by shannon capacity theorem as follows:
Figure FDA0003439331970000036
the optimization problem is represented as:
Figure FDA0003439331970000041
Figure FDA0003439331970000042
Figure FDA0003439331970000043
wherein the content of the first and second substances,
Figure FDA0003439331970000044
represents the signal-to-interference-and-noise ratio, p, of the kth user signal in cell j received by the base station in cell jliDenotes the data power, V, allocated to user i in cell jjkIndicating that the combined vector is received,
Figure FDA0003439331970000045
representing the estimated value of the channel between the base station of cell j and the kth user in the cell,
Figure FDA0003439331970000046
an estimation error matrix representing a channel matrix,
Figure FDA0003439331970000047
which is indicative of the power of the noise,
Figure FDA0003439331970000048
is MjIdentity matrix of order, EEULRepresenting system energy efficiency, τ, in a massive MIMO uplink transmission processcDenotes the total length of the coherent block, τuIndicating the length of the coherent block for transmitting uplink data, M being the number of antennas provided for the base station, prPower consumed by the radio link of the antenna for uplink data transmission, psStatic electrical power consumption, p, consumed during data transmission for large-scale MIMO systemsmaxIs the maximum electric power, R, of all users in each cell during uplink transmissionminIs the minimum user transmission rate limit we consider.
5. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 4, wherein step S5 specifically includes:
s51, initializing the parameters of the particles by using the rand function and the upper and lower limits of the set user data power;
s52, taking the energy efficiency function as a fitness function, calculating the fitness value of the initialized particles, and initializing a local optimal solution and a global optimal solution;
s53, judging the iteration termination condition, terminating the iteration when the iteration reaches the initial set maximum iteration number or the fitness value tends to be stable, otherwise, continuing the following steps;
s54, updating the speed and the position of the particle, simultaneously performing boundary processing on the particle beyond the boundary range, and recalculating the fitness value of the updated particle;
and S55, storing the result of each iteration, comparing the results after each iteration, and updating the local optimal solution and the global optimal solution.
And S56, when the iteration termination condition is reached, obtaining the global optimal value, namely the required optimal user data power allocation vector.
6. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 5,
in step S52, the fitness function is set as an energy efficiency expression:
Figure FDA0003439331970000051
and substituting the position parameter of each particle as an input parameter into the formula to solve a corresponding fitness value, finishing the first iteration process, and setting the initialized particle parameters as a local optimal value and a global optimal value.
7. The energy efficiency optimization method based on adaptive particle swarm power distribution according to claim 5, wherein in step S54, the particle velocity and position are updated and boundary processing is performed as follows:
the update formula of the particle speed and the position is as follows:
Vj(t+1)=Vj(t)+c1*rand*(Gbestj-popj(t))+c2*rand*(Zbest-popj(t))
popj(t+1)=popj(t)+ω*Vj(t+1)
Figure FDA0003439331970000061
in the formula, Vj(t +1) represents the search speed, V, of the jth particle at the (t +1) th iterationj(t) represents the search speed of the jth particle at the tth iteration, c1 is an individual learning factor, the iteration representing the speed is only related to the historical position of the particle, c2 is a social learning factor, the relation between the iteration representing the speed and the historical position of the whole particle swarm, and omega represents a weight factor; gbestjFor the current optimal position of the jth particle, Zbest is the so far optimal particle position, pop, of the entire populationj(t) denotes the time of the tth iteration representing the jth particleA location parameter of (a);
the boundary processing method of the particle speed and the particle position comprises the following steps:
Figure FDA0003439331970000062
Figure FDA0003439331970000063
v (i, j) represents the search speed of the ith particle in the j-dimension direction; pop (i, j) represents the displacement size of the ith particle in the j-dimension direction.
8. An energy efficiency optimization device based on adaptive particle swarm power distribution is characterized by comprising the following components:
the initialization unit is used for initializing the positions of all base stations and the positions of all users in a cell;
a calculating unit, configured to calculate a large-scale fading factor β and a covariance matrix R of a channel according to the positions of the base stations and the positions of the users;
the channel estimation unit is used for receiving the uplink pilot signal, and simultaneously estimating the channel by adopting an MMSE (minimum mean square error) method in combination with the large-scale fading factor beta and the covariance matrix R of the channel;
the model establishing unit is used for deducing a signal to interference plus noise ratio (SINR) expression of the system according to the estimated channel, calculating frequency efficiency according to the Shannon capacity theorem, and establishing an energy efficiency optimization model of the system by combining a power consumption model of the system;
and the power distribution unit is used for carrying out data power distribution on the user on the basis of fixed pilot frequency power by utilizing a self-adaptive particle swarm algorithm according to the optimization objective function for the energy efficiency optimization model.
9. An energy efficiency optimization device based on adaptive particle swarm power distribution, comprising a memory and a processor, wherein the memory stores a computer program, and the computer program can be executed by the processor to implement the energy efficiency optimization method based on adaptive particle swarm power distribution according to any one of claims 1 to 7.
CN202111628181.9A 2021-12-28 2021-12-28 Energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution Active CN114302487B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111628181.9A CN114302487B (en) 2021-12-28 2021-12-28 Energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111628181.9A CN114302487B (en) 2021-12-28 2021-12-28 Energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution

Publications (2)

Publication Number Publication Date
CN114302487A true CN114302487A (en) 2022-04-08
CN114302487B CN114302487B (en) 2024-03-05

Family

ID=80971298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111628181.9A Active CN114302487B (en) 2021-12-28 2021-12-28 Energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution

Country Status (1)

Country Link
CN (1) CN114302487B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116390135A (en) * 2023-04-26 2023-07-04 中南大学 Self-return millimeter wave cellular network communication method based on dynamic time division duplex communication

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103561430A (en) * 2013-11-20 2014-02-05 东南大学 Method for balancing energy efficiency and spectral efficiency
CN105722203A (en) * 2016-02-02 2016-06-29 东南大学 High energy efficiency power distribution method based on particle swarm optimization for large-scale antenna system
CN107733488A (en) * 2017-10-16 2018-02-23 中南大学 Water injection power distribution improved method and system in a kind of extensive mimo system
US10069592B1 (en) * 2015-10-27 2018-09-04 Arizona Board Of Regents On Behalf Of The University Of Arizona Systems and methods for securing wireless communications
CN109474317A (en) * 2019-01-07 2019-03-15 河南理工大学 MR pre-processes the lower extensive MIMO bidirectional relay system power distribution method of hardware damage
CN109890075A (en) * 2019-03-29 2019-06-14 中南大学 A kind of suppressing method of extensive mimo system pilot pollution, system
CN110149130A (en) * 2019-05-16 2019-08-20 杭州电子科技大学 A kind of extensive mimo system upgoing energy efficiency optimization method based on PSO
CN110190879A (en) * 2019-04-30 2019-08-30 杭州电子科技大学 Efficiency optimization method based on the low extensive mimo system of Precision A/D C
CN110505643A (en) * 2019-09-23 2019-11-26 杭州电子科技大学 Extensive mimo system uplink efficiency optimization method based on simulated annealing
CN110808765A (en) * 2019-08-30 2020-02-18 南京航空航天大学 Power distribution method for optimizing spectrum efficiency of large-scale MIMO system based on incomplete channel information
CN110958102A (en) * 2019-12-03 2020-04-03 中南大学 Pilot pollution suppression method based on pilot distribution and power control joint optimization

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103561430A (en) * 2013-11-20 2014-02-05 东南大学 Method for balancing energy efficiency and spectral efficiency
US10069592B1 (en) * 2015-10-27 2018-09-04 Arizona Board Of Regents On Behalf Of The University Of Arizona Systems and methods for securing wireless communications
CN105722203A (en) * 2016-02-02 2016-06-29 东南大学 High energy efficiency power distribution method based on particle swarm optimization for large-scale antenna system
CN107733488A (en) * 2017-10-16 2018-02-23 中南大学 Water injection power distribution improved method and system in a kind of extensive mimo system
CN109474317A (en) * 2019-01-07 2019-03-15 河南理工大学 MR pre-processes the lower extensive MIMO bidirectional relay system power distribution method of hardware damage
CN109890075A (en) * 2019-03-29 2019-06-14 中南大学 A kind of suppressing method of extensive mimo system pilot pollution, system
CN110190879A (en) * 2019-04-30 2019-08-30 杭州电子科技大学 Efficiency optimization method based on the low extensive mimo system of Precision A/D C
CN110149130A (en) * 2019-05-16 2019-08-20 杭州电子科技大学 A kind of extensive mimo system upgoing energy efficiency optimization method based on PSO
CN110808765A (en) * 2019-08-30 2020-02-18 南京航空航天大学 Power distribution method for optimizing spectrum efficiency of large-scale MIMO system based on incomplete channel information
CN110505643A (en) * 2019-09-23 2019-11-26 杭州电子科技大学 Extensive mimo system uplink efficiency optimization method based on simulated annealing
CN110958102A (en) * 2019-12-03 2020-04-03 中南大学 Pilot pollution suppression method based on pilot distribution and power control joint optimization

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
DENG, HG; ZHANG, MN; DENG, HG; MA, SS; LIU, G: "Preamble-Based Channel Estimation Corrected by Equalization in FBMC/OQAM System", 《 2018 IEEE INTERNATIONAL CONFERENCE ON COMPUTER AND COMMUNICATION ENGINEERING TECHNOLOGY (CCET)》, 6 June 2019 (2019-06-06) *
XIAOSHU ZHU;JIE ZHANG ;JUNHONG FENG: "Multiobjective Particle Swarm Optimization Based on PAM and Uniform Design", MATHEMATICAL PROBLEMS IN ENGINEERING, vol. 2015, 22 July 2015 (2015-07-22) *
张继荣;孟繁克;王晟寰: "改进粒子群算法的D2D功率分配", 《西安邮电大学学报》, vol. 26, no. 2, 10 March 2021 (2021-03-10) *
徐进,费少梅,张树有,施岳定: "自适应粒子群求解资源动态分配项目调度问题", 《计算机集成制造系统》, vol. 17, no. 8, 15 August 2011 (2011-08-15) *
杜炜民: "大规模MIMO系统中频带效率和能量效率的优化", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》, 15 March 2017 (2017-03-15) *
路坷平: "大规模MIMO系统导频数据功率比及信道协方差矩阵信息的研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》, 15 February 2019 (2019-02-15) *
邓宏贵,赵俊值: "基于多普勒辅助迭代的高速OFDM系统信道估计方法", 《计算机应用与软件》, vol. 30, no. 1, 15 January 2013 (2013-01-15) *
邓宏贵,赵俊值: "基于多普勒辅助迭代的高速OFDM系统信道估计方法", 计算机应用软件, vol. 30, no. 1, 15 January 2013 (2013-01-15) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116390135A (en) * 2023-04-26 2023-07-04 中南大学 Self-return millimeter wave cellular network communication method based on dynamic time division duplex communication
CN116390135B (en) * 2023-04-26 2024-02-02 中南大学 Self-return millimeter wave cellular network communication method based on dynamic time division duplex communication

Also Published As

Publication number Publication date
CN114302487B (en) 2024-03-05

Similar Documents

Publication Publication Date Title
US8755358B2 (en) Wireless base station device, terminal, and wireless communication method
Ge et al. Deep reinforcement learning for distributed dynamic MISO downlink-beamforming coordination
EP2321912B1 (en) Beamforming systems and method
Seifi et al. Adaptive multicell 3-D beamforming in multiantenna cellular networks
US11057080B2 (en) Software-defined massive multi-input multi-output (MIMO)
Sabbagh et al. Pilot allocation and sum-rate analysis in cell-free massive MIMO systems
US10666329B2 (en) Method and device for beam forming
CN110958102B (en) Pilot pollution suppression method based on pilot distribution and power control joint optimization
CN106160806B (en) Method and apparatus for performing interference coordination in wireless communication system
Banoori et al. Pilot contamination mitigation under smart pilot allocation strategies within massive MIMO-5G system
CN114302487B (en) Energy efficiency optimization method, device and equipment based on self-adaptive particle swarm power distribution
KR20140089890A (en) Method and apparatus for adaptive inter-cell interference canellation for wireless communication system
CN110677858A (en) Transmission power and computing resource allocation method based on task migration period of Internet of things
KR101571103B1 (en) Apparatus and method for transmitting linearly in distributed mimo system
Hawej et al. Iterative weighted nuclear norm minimization-based channel estimation for massive multi-user MIMO systems
CN107872255B (en) Pilot frequency scheduling method suitable for large-scale MIMO cellular mobile communication network
Krunz et al. Online Reinforcement Learning for Beam Tracking and Rate Adaptation in Millimeter-wave Systems
Wang et al. Deep transfer reinforcement learning for beamforming and resource allocation in multi-cell MISO-OFDMA systems
Liu et al. A Reinforcement Learning Approach for Energy Efficient Beamforming in NOMA Systems
Nugraha et al. Block diagonalization precoding and power allocation for clustering small-cell networks
Akbarpour-Kasgari et al. Deep Reinforcement Learning in mmW-NOMA: Joint Power Allocation and Hybrid Beamforming
Kim User scheduling and grouping in massive MIMO broadcast channels with heterogeneous users
Godana et al. Coordinated beamforming in multicell networks with Channel State Information exchange delays
Tekgul et al. Uplink-downlink joint antenna optimization in cellular systems with sample-efficient learning
Sun et al. Predictive Resource Allocation in mmWave Systems with Rotation Detection

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