CN103067062A - Base station antenna selecting method based on interference alignment in multi-cell system - Google Patents
Base station antenna selecting method based on interference alignment in multi-cell system Download PDFInfo
- Publication number
- CN103067062A CN103067062A CN2013100256634A CN201310025663A CN103067062A CN 103067062 A CN103067062 A CN 103067062A CN 2013100256634 A CN2013100256634 A CN 2013100256634A CN 201310025663 A CN201310025663 A CN 201310025663A CN 103067062 A CN103067062 A CN 103067062A
- Authority
- CN
- China
- Prior art keywords
- base station
- antenna
- matrix
- omega
- prime
- 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
Images
Abstract
The invention discloses a base station antenna selecting method based on interference alignment in a multi-cell system, and mainly solves the problems that an existing multi-antenna interference alignment technology is relatively high in system hardware cost, and the base station antenna selecting method based on interference alignment can only be used in a network comprising 3 cells. The method comprises the following implementation steps of: estimating state information of channels from various base stations to users in the resident cells; initiating base station antenna selection schemes of the various cells; along a system and capacity increase optimizing direction, under the condition that other base station antenna selection schemes are unchanged, determining the optimal base station selection schemes of the various cells in sequence by using partial iterative interference alignment; and performing multi-cell interference coordination by adopting traditional iterative interference alignment. The base station antenna selecting method is applicable to an LTE (Long Term Evolution) network oriented to commercial use with a plurality of cells, is low in computation complexity, can reduce the system hardware cost, and realizes compromise between system and capacity and computation complexity.
Description
Technical field
The invention belongs to wireless communication technology field, further relate to the antenna for base station system of selection of aliging based on interference in a kind of multi-cell system, can be used for the LTE network that business is used.
Background technology
In order to improve the capacity of cellular system, the Next-Generation Wireless Communication Systems take the LTE network as representative requires to realize as far as possible that the full rate of neighbor cell is multiplexing, i.e. the spectrum reuse factor is 1.Therefore, with respect to noise and decline, presence of intercell interference becomes the principal element that affects systematic function.How effectively carrying out many Inter-Cell Interference Coordination, the interference of particularly eliminating Cell Edge User has caused widely to be paid close attention to.Disturbing alignment is a kind of interference coordination technique that proposes recent years, it can will snap to from the interference signal of other base stations a sub spaces that receives signal space at receiver user, useful signal from respective base station is projected to another interference-free subspace, thereby effectively disturb inhibition, improve power system capacity.
The scholars such as K.Gomadam are at article " Approaching the Capacity of Wireless Networksthrough Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, propose a kind of distributed iterative 2008:1-5) and disturbed alignment schemes (be called traditional iteration and disturb alignment), be applicable to the interference coordination of many residential quarters multiple-input and multiple-output (MIMO) network.The method is the pre-coding matrix of random initializtion base station at first, then according to disturbing the interference of revealing minimum principle Iterative Design user to suppress the pre-coding matrix of matrix and base station until convergence.The weak point that the method exists is: although user's antenna number is generally less in the mimo system, but a fairly large number of antenna of base station possible configuration, the method is not considered the increase along with antenna for base station quantity, the radio frequency link that comprises the elements such as digital-to-analogue/analog to digital converter, low-converter, low noise amplifier also will increase gradually, thereby increases the problem of system hardware cost.
The people such as J.G.Klotz are at article " Antenna Selection Criteria for Interference Alignment " (Proceedings of IEEE International Symposium on Personal, Indoor and Mobile RadioCommunications, Istanbul, Turkey, proposed 2010:527-531) for adopting non-iteration to disturb the antenna for base station system of selection of 3 cell systems of alignment, to reduce the system hardware cost, obtain simultaneously system and volumetric properties and promote.The main performing step of the method is: the first, and all days line options situation that all alternative antenna subsets of utilizing poor searching method to travel through 3 base stations are combined into; The second, for each day line options situation, carry out non-iteration and disturb alignment, and computing system and capacity or each user's useful and interference signal vector space chordal distance with; The 3rd, select the antenna selecting plan of each base station optimum according to maximization system and capacity criterion or maximization chordal distance criterion, the antenna of choosing is connected on the corresponding radio frequency link; The 4th, realize non-iteration interference alignment.The weak point that the method exists is: disturb alignment in view of non-iteration and can only finish the interference coordination of 3 residential quarters, this method is only applicable to comprise the system of 3 residential quarters, therefore uses to have limitation, can not expand to the system of many residential quarters.
Summary of the invention
The object of the invention is to overcome above-mentioned the deficiencies in the prior art, propose in a kind of multi-cell system based on the antenna for base station system of selection of disturbing alignment, reduce the system hardware cost, obtain a day line options gain, simultaneously with range of application from 3 cell extension to a plurality of residential quarters.
For achieving the above object, the thinking of the inventive method is:
Do not fix and more characteristics for carrying out the multiplexing number of cells of full rate in the LTE system, the inventive method adopts iteration to disturb alignment to carry out the interference coordination of a plurality of cell systems; Owing to determining that the antenna for base station selection scheme of each residential quarter optimum is a combinatorial optimization problem, finding the solution this combinatorial optimization problem need to repeatedly search for various days line options situations that are combined into by each cell base station, therefore in each day line options situation search procedure, adopt the traditional iteration that comprises repeatedly iteration to disturb alignment can cause the computation complexity of method higher; The inventive method adopts following two modes to reduce respectively the number of times of various days line options situations of search and the computation complexity in each search procedure: first, optimization direction along system and capacity increase, in the situation that keep other antenna for base station selection schemes constant, successively the antenna for base station selection scheme of each residential quarter is optimized, thereby reduces the number of times of finding the solution various days line options situations of search in the combinatorial optimization problem; The second, when optimizing the antenna for base station selection scheme of each residential quarter, adopt the part iteration of only carrying out the several times iteration to disturb alignment, then computing system and capacity, rather than finish traditional iteration and disturb alignment, thereby reduce the computation complexity of each search procedure; According to the antenna for base station selection scheme of each residential quarter optimum, the inventive method adopts traditional iteration to disturb alignment to realize many Inter-Cell Interference Coordination at last.
Concrete steps of the present invention are as follows:
(1) estimates that each base station is to the channel condition information of this community user
Base station k sends pilot frequency information for this community user k, and user k estimates each antenna of base station k to the channel condition information of self, and feeds back to base station k, k=1 wherein ..., B, B are the residential quarter number; Base station k determines channel matrix according to the channel condition information of this community user feedback
Wherein N and M represent respectively the antenna number of base station k and user k,
The expression complex field;
(2) each cell-site antenna selection scheme of initialization
The antenna selecting plan ω of initialization base station k
kFor having the matrix of maximum F-norm
Corresponding alternative antenna is selected vector, ω
kComputing formula as follows:
Wherein, k=1 ..., B,
Expression line options scheme on the same day is φ
lThe time according to channel matrix H
KkThe base station k that determines is to the channel matrix of this community user, φ
lL alternative antenna of expression base station selected vector,
Select the index variables of vector for alternative antenna, N
fBe the rf chain way of base station,
Expression is chosen N from the N number
fThe number of combinations of number, || ||
FThe F-norm of representing matrix; The antenna selecting plan of each base station is combined into B cell-site antenna selection set omega={ ω
1, ω
2..., ω
B;
(3) along the optimization direction of system and capacity increase, in the situation that keep other antenna for base station selection schemes constant, utilize the part iteration to disturb alignment to determine successively the antenna for base station selection scheme of each residential quarter optimum:
3.1) to select the index variables l of vector be 1 to the initialization alternative antenna, system and the capacity C of simultaneously initialization maximum
MaxBe 0;
3.2) with the antenna for base station selection scheme ω of residential quarter 1
1Be updated to alternative antenna and select vectorial φ
l, keep the antenna selecting plan of other base stations constant, simultaneously the 1st element in B the cell-site antenna selection set omega is updated to φ
l
3.3) according to above-mentioned antenna selecting plan, carry out the part iteration and disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix, are implemented as follows:
3.3a) pre-coding matrix of random initializtion base station k
K=1 wherein ..., B, d represent that the base station sends the dimension of data flow;
3.3b) initialization iterations indicator variable c is 1;
3.3c) determine the interference covariance matrix of user k ' when the c time iteration
K '=1 wherein ..., B;
3.3d) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The interference of design user k ' when the c time iteration suppresses matrix
K '=1 wherein ..., B, υ
dD characteristic value characteristic of correspondence vector of { } representing matrix minimum;
3.3e) suppress matrix according to above-mentioned interference
Determine the interference covariance matrix of base station k when the c time iteration
K '=1 wherein ..., B, k=1 ..., B;
3.3f) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The pre-coding matrix of design base station k when the c time iteration
K=1 wherein ..., B;
3.3g) value of iterations indicator variable c adds 1, if c≤N
Ite, then return step 3.3c, proceed iteration; Otherwise, jump to step 3.4, obtain N
IteThe pre-coding matrix of base station k after the inferior iteration
Suppress matrix with the interference of user k '
K=1 wherein ..., B, k '=1 ..., B, N
IteDisturb the iterations of alignment for the part iteration;
3.4) determine that base station k ' is to the equivalent channel of user k '
And the interference covariance matrix of user k '
Wherein
Be illustrated in antenna selecting plan ω
K 'Base station k ' arrives the channel matrix of user k ' in the situation,
Expression N
IteThe pre-coding matrix of base station k ' after the inferior iteration, k '=1 ..., B;
3.5) equivalent channel of disturbing alignment to obtain according to the part iteration
Disturb and suppress matrix
And interference covariance matrix
Computing system and capacity
Wherein
Be the capacity of user k ', k '=1 wherein ..., B;
3.6) if C
Ω>C
Max, C then is set
Max=C
Ω, and the optimum alternative antenna of order is selected the index variables l of vector
*=l, namely optimum alternative antenna is selected vector
Otherwise, C
MaxRemain unchanged;
3.7) alternative antenna selects the value of the index variables l of vector to add 1, if
Then return step 3.2, continue other alternative antenna of traversal and select vector; Otherwise, skip to step 3.8;
3.9) according to the described method of step 3.1 to 3.8, determine that successively residential quarter 2 is to the antenna for base station selection scheme of residential quarter B optimum
3.10) the antenna for base station selection scheme of above-mentioned each residential quarter optimum is combined into B optimum cell-site antenna selection set
(4) adopt traditional iteration to disturb alignment to carry out many coordinating district interferences:
4.1) select set omega according to B cell-site antenna of above-mentioned optimum
*, carrying out iterations is N
ITraditional iteration disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix;
4.2) disturb the pre-coding matrix of alignment design and interference to suppress matrix according to above-mentioned traditional iteration, each cell base station and user transmit and receive data respectively, carry out many coordinating district interferences, eliminate presence of intercell interference.
The present invention has following advantage:
1. the present invention is owing to utilizing iteration to disturb alignment to carry out many coordinating district interferences, and the antenna for base station of realizing each residential quarter optimum is selected, overcome the existing shortcoming that can only in comprising 3 subzone networks, use based on the antenna for base station system of selection of disturbing alignment, with range of application from 3 cell extension to a plurality of residential quarters, the LTE network of using applicable to the business with a plurality of residential quarters;
2. the present invention is owing to the optimization direction along system and capacity increase, in the constant situation of the antenna for base station selection scheme that keeps other residential quarters, determine successively the antenna for base station selection scheme of each residential quarter optimum, reduced the number of times of searching for various days line options situations, and in the search procedure of each day line options situation, adopt the less part iteration of iterations to disturb alignment, thereby greatly reduced computation complexity;
3. the present invention has reduced the rf chain way of base station because antenna for base station is selected, thereby reduces the system hardware cost, obtains simultaneously a day line options gain;
4. the present invention can the adjustment member iteration disturb the iterations of alignment, thereby reaches the compromise of system and capacity and computation complexity.
Description of drawings
Fig. 1 is the B cell system illustraton of model that the present invention is suitable for;
Fig. 2 is general flow chart of the present invention;
Fig. 3 is the sub-process figure that determines the antenna for base station selection scheme of each residential quarter optimum among the present invention;
Fig. 4 is the sub-process figure that the part iteration is disturbed alignment among the present invention;
Fig. 5 is the present invention and the system of random antenna for base station system of selection and poor search base station antenna selecting method and the simulation comparison figure that capacity changes with signal to noise ratio.
Embodiment
Below in conjunction with accompanying drawing the present invention is done further detailed description.
With reference to Fig. 1, in the B cell system, user of the base station of each residential quarter service, the antenna number of base station and rf chain way are respectively N and N
f, user's antenna number and radio frequency link number average are M.The sequence number of residential quarter represents with k, and k=1 ..., B, then base station k and user k represent respectively base station and the user of residential quarter k.The system that a plurality of users are arranged for each residential quarter, for example can adopting, time-multiplexed mode is reduced to the system that each time slot only has a user.Suppose N>N
f, N need to be chosen in each base station from all N antenna
fIndividual antenna, and be connected on the corresponding radio frequency link, so the alternative antenna of base station set Φ comprises
Individual alternative antenna is selected vector, can be expressed as
Wherein
Expression is chosen N from the N number
fThe number of combinations of number, φ
lRepresent that l alternative antenna select vector, it is the long column vector of N, and its element value is 0 or 1, and l is the index variables that alternative antenna is selected vector.Sequence number is that the element of i is that 1 to represent i antenna of base station selected in the column vector, otherwise not selected.
Expression base station k is to the channel matrix of user k ', wherein ω
kThe antenna selecting plan that expression base station k chooses from alternative antenna is selected to gather,
The expression complex field.
With reference to Fig. 2, performing step of the present invention comprises as follows:
1a) base station k sends to user k
Inferior pilot frequency information is chosen each time some antennas and they is connected to radio frequency link transmission pilot frequency information, wherein chooses N the 1st time from all N antenna
fIndividual antenna, the 2nd time from remaining N-N
fChoose N in the individual antenna
fIndividual antenna, the 3rd time from remaining N-2N
fChoose N in the individual antenna
fIndividual antenna, until the
Inferiorly choose last remaining antenna, thereby guarantee that all antennas all can once be connected to radio frequency link at certain and send pilot frequency information, k=1 wherein ..., B,
Representative rounds up to scalar x;
1b) user k estimates base station k respective antenna to the channel condition information of self each antenna by the pilot frequency information that receives at every turn, and with corresponding channel matrix feedback to base station k, k=1 wherein ..., B;
1c) base station k obtains its all antennas to the channel coefficients of each antenna of user according to the channel matrix of each time feedback, and channel coefficients is formed channel matrix
K=1 wherein ..., B.
Step 2, each cell-site antenna selection scheme of initialization.
2a) according to channel matrix H
KkSelect vectorial φ with alternative antenna
l, determine to work as φ
lBase station k is to the channel matrix of this community user during as antenna selecting plan
Wherein, k=1 ..., B,
Diag{ φ
lRepresent that the element on the leading diagonal is column vector φ
lThe diagonal matrix of corresponding element, function
The submatrix that expression is made of the non-zero row of matrix;
2b) the antenna selecting plan ω of initialization base station k
kMatrix for the F-norm of maximum
Corresponding alternative antenna is selected vector, and computing formula is as follows:
Wherein, k=1 ..., B, || ||
FThe F-norm of representing matrix;
2c) antenna selecting plan with above-mentioned each base station is combined into B cell-site antenna selection set omega={ ω
1, ω
2..., ω
B.
Step 3 along the optimization direction of system and capacity increase, in the situation that keep other antenna for base station selection schemes constant, utilizes the part iteration to disturb alignment to determine successively the antenna for base station selection scheme of each residential quarter optimum.
With reference to Fig. 3, being achieved as follows of this step:
3.1) to select the index variables l of vector be 1 to the initialization alternative antenna, system and the capacity C of simultaneously initialization maximum
MaxBe 0;
3.2) with the antenna for base station selection scheme ω of residential quarter 1
1Be updated to alternative antenna and select vectorial φ
l, keep the antenna selecting plan of other base stations constant, simultaneously the 1st element in B the cell-site antenna selection set omega is updated to φ
l
3.3) according to above-mentioned antenna selecting plan, carry out the part iteration and disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix:
With reference to Fig. 4, being implemented as follows of this step:
3.3a) pre-coding matrix of random initializtion base station k
And satisfy
K=1 wherein ..., B, ()
HThe conjugate transpose of representing matrix, I
nThe unit matrix of expression n * n dimension, d represents that the base station sends the dimension of data flow;
3.3b) initialization iterations indicator variable c is 1;
3.3c) according to scholars such as K.Gomadam at article " Approaching the Capacity of WirelessNetworks through Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, training method 2008:1-5) is determined the interference covariance matrix of user k ' when the c time iteration
Wherein, k '=1 ..., B, P
kThe transmitting power of expression base station k,
The pre-coding matrix of base station k when being illustrated in the c-1 time iteration;
3.3d) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The interference of design user k ' when the c time iteration suppresses matrix
K '=1 wherein ..., B, υ
dD characteristic value characteristic of correspondence vector of { } representing matrix minimum;
3.3e) according to scholars such as K.Gomadam at article " Approaching the Capacity of WirelessNetworks through Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, training method 2008:1-5) is determined the interference covariance matrix of base station k when the c time iteration
Wherein, k=1 ..., B, P
K 'The transmitting power of expression user k ', value equates with the transmitting power of this cell base station;
3.3f) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The pre-coding matrix of design base station k when the c time iteration
K=1 wherein ..., B;
3.3g) value of iterations indicator variable c adds 1, if c≤N
Ite, then return step 3.3c, proceed iteration; Otherwise, jump to step 3.4, obtain N
IteThe pre-coding matrix of base station k after the inferior iteration
Suppress matrix with the interference of user k '
K=1 wherein ..., B, k '=1 ..., B, N
IteDisturb the iterations of alignment for the part iteration, and satisfy N
Ite≤ N
I, N
IThe iterations empirical value that the expression tradition disturbs the alignment convergence to need is preferably 30;
3.4) according to scholars such as K.Gomadam at article " Approaching the Capacity of WirelessNetworks through Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, training method 2008:1-5) determines that base station k ' is to the equivalent channel of user k '
And the interference covariance matrix of user k '
Wherein, k '=1 ..., B, wherein
Be illustrated in antenna selecting plan ω
K 'Base station k ' is to the channel matrix of user k ' in the situation;
3.5) equivalent channel of disturbing alignment to obtain according to the part iteration
Disturb and suppress matrix
And interference covariance matrix
Computing system and capacity C
Ω:
Wherein, k '=1 ..., B, σ
2The variance of white Gaussian noise in the channel, || the determinant of representing matrix, bodge are bit/s/Hz;
3.6) if C
Ω>C
Max, C then is set
Max=C
Ω, and the optimum alternative antenna of order is selected the index variables l of vector
*=l, namely optimum alternative antenna is selected vector
Otherwise, C
MaxRemain unchanged;
3.7) alternative antenna selects the value of the index variables l of vector to add 1, if
Then return step 3.2, other alternative antenna that continue among the basic alternative antenna set of the traversal Φ are selected vector; Otherwise, skip to step 3.8;
3.9) according to the described method of step 3.1 to 3.8, determine that successively residential quarter 2 is to the antenna for base station selection scheme of residential quarter B optimum
3.10) the antenna for base station selection scheme of above-mentioned each residential quarter optimum is combined into B optimum cell-site antenna selection set
Step 4 adopts traditional iteration to disturb alignment to carry out many coordinating district interferences.
4.1) select set omega according to B cell-site antenna of above-mentioned optimum
*, carrying out iterations is N
ITraditional iteration disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix:
4.1a) pre-coding matrix of random initializtion base station k
And satisfy
K=1 wherein ..., B, ()
HThe conjugate transpose of representing matrix, I
nThe unit matrix of expression n * n dimension, d represents that the base station sends the dimension of data flow;
4.1b) initialization iterations indicator variable m is 1;
4.1c) according to scholars such as K.Gomadam at article " Approaching the Capacity of WirelessNetworks through Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, training method 2008:1-5) is determined the interference covariance matrix of user k ' when the m time iteration
Wherein, k '=1 ..., B, P
kThe transmitting power of expression base station k,
The pre-coding matrix of base station k when being illustrated in the m-1 time iteration,
Be illustrated in antenna selecting plan
Base station k is to the channel matrix of user k ' in the situation;
4.1d) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The interference of design user k ' when the m time iteration suppresses matrix
K '=1 wherein ..., B, υ
dD characteristic value characteristic of correspondence vector of { } representing matrix minimum;
4.1e) according to scholars such as K.Gomadam at article " Approaching the Capacity of WirelessNetworks through Distributed Interference Alignment " (Proceedings of the IEEE GlobalTelecommunications Conference, Miami, USA, training method 2008:1-5) is determined the interference covariance matrix of base station k when the m time iteration
Wherein, k=1 ..., B, P
kThe transmitting power of ' expression user k ', its value equates with the transmitting power of this cell base station;
4.1f) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix
The pre-coding matrix of design base station k when the m time iteration
K=1 wherein ..., B;
4.1g) value of iterations variable m adds 1, if m≤N
I, then return step 4.1c, proceed iteration; Otherwise iterative process finishes, and obtains the pre-coding matrix of base station k
Suppress matrix with the interference of user k '
K=1 wherein ..., B, k '=1 ..., B jumps to step 4.2;
4.2) disturb the pre-coding matrix of alignment design and interference to suppress matrix according to above-mentioned traditional iteration, each cell base station and user transmit and receive data respectively, carry out many coordinating district interferences, eliminate presence of intercell interference.
Advantage of the present invention further specifies by following computation complexity analysis and simulation result:
1) computation complexity analysis
Before analyzing computation complexity of the present invention, the computation complexity of the poor search base station antenna selecting method of given first.Poor search base station antenna selecting method is the most direct a kind of method, namely travel through by all alternative antenna of all base stations and select various days synthetic line options situations of Vector Groups, then disturb alignment to select so that the antenna for base station selection scheme of a kind of each residential quarter of system and capacity maximum by traditional iteration.Owing in the system B base station arranged, and each base station having
Individual alternative antenna is selected vector, so the number of times of various days line options situations of the method search is
In each search procedure, carrying out iterations is N
ITraditional iteration disturb alignment, disturb the computation complexity of alignment to be in each search procedure:
Therefore, the poor total computation complexity of search base station antenna selecting method is:
The present invention is along the optimization direction of system and capacity increase, in the situation that keep other antenna for base station selection schemes constant, successively the antenna selecting plan of each base station is optimized, the number of times that the not possible situation of limit institute, so the present invention is searched for various days line options situations is
In each search procedure, the present invention adopts the part iteration to disturb alignment, has reduced the iterations that disturbs alignment, disturbs the computation complexity of alignment to be in each search procedure:
Disturb the computation complexity that aligns to compare in the computation complexity of matrix F-norm and the each search procedure in the step 2 very little, so this part complexity can be ignored.Consideration step 4 is carried out the computation complexity that traditional iteration is disturbed alignment, and the total computation complexity of the present invention can be expressed as:
Further abbreviation is:
Wherein, θ=N
Ite/ N
IBe the part iteration factor, its span is 0<θ≤1.Because in B>1 and
Condition under
Less than
Especially B or
When larger
Much smaller than
Therefore the present invention is lower with respect to poor search base station antenna selecting method computation complexity.
2) emulation experiment
2.1) simulated conditions:
Adopt the MATLAB simulation software, emulation 1000 times.The rf chain way N of the antenna number N=4 of residential quarter number B=4, base station, base station
f=3, the iterations of user's antenna and rf chain way M=3, the dimension d=1 that sends data flow, traditional iteration interference alignment is N
I=30.The transmitting power of supposing each cell base station equates, all is P, and the variances sigma of Gaussian noise
2=1, then power P has represented the signal to noise ratio snr of each community user in the system.
2.2) emulation content and result:
Disturb the iterations N of alignment in the part iteration
Ite=1,5,10, under 30 conditions, system and capacity to the present invention and random antenna for base station system of selection and poor search base station antenna selecting method carry out Monte-Carlo Simulation with the signal to noise ratio snr situation of change, obtain corresponding simulation comparison figure, as shown in Figure 5, wherein abscissa represents signal to noise ratio snr, and ordinate represents system and capacity.
As can be seen from Figure 5, system and the volumetric properties of random antenna for base station system of selection are the poorest, because random device is just chosen at random antenna for base station and accessed the respective radio-frequency link, are difficult to obtain a day line options gain.Poor search base station antenna selecting method has traveled through by all alternative antenna of all base stations and has selected various days synthetic line options situations of Vector Groups, and adopts traditional iteration to disturb the method for alignment, has therefore obtained optimum system and capacity.When signal to noise ratio snr is increased to 30dB from 0dB, system of the present invention and capacity all a little less than poor search base station antenna selecting method apparently higher than system and the capacity of random antenna for base station system of selection.The present invention is at N
IteConcrete system and volumetric properties situation when getting different value: work as N
IteWhen getting minimum value 1, poor-performing of the present invention, the optimal performance that the poor search base station antenna selecting method of distance obtains is far away, because N
IteToo little, the part iteration has been disturbed the iteration of only aliging once, but the performance of this moment still is better than the performance of random antenna for base station system of selection; Work as N
IteWhen getting maximum 30, the present invention obtains the optimal performance that can reach, and this optimal performance is close to the performance of poor search base station antenna selecting method.Work as N
IteGet respectively 5 and at 10 o'clock, performance of the present invention is with respect to N
IteThe performance of getting at 1 o'clock is obviously promoted.N particularly
IteGot 10 o'clock, the present invention can obtain almost and N
IteGet 30 o'clock same performances, and iterations is with respect to N
IteGet 30 o'clock less, show N
Ite=10 can obtain the compromise of system and capacity and computation complexity.
The analysis showed that of simulation result and computation complexity, the present invention can reduce computation complexity effectively, obtain simultaneously close to poor search base station antenna selecting method and volumetric properties.In addition, the present invention disturbs the iterations of alignment, the compromise of realization system and capacity and computation complexity by the adjustment member iteration.
Claims (10)
- In the multi-cell system based on the antenna for base station system of selection of disturbing alignment, comprise the steps:(1) estimates that each base station is to the channel condition information of this community userBase station k sends pilot frequency information for this community user k, and user k estimates each antenna of base station k to the channel condition information of self, and feeds back to base station k, k=1 wherein ..., B, B are the residential quarter number; Base station k determines channel matrix according to the channel condition information of this community user feedback Wherein N and M represent respectively the antenna number of base station k and user k, The expression complex field;(2) each cell-site antenna selection scheme of initializationThe antenna selecting plan ω of initialization base station k kFor having the matrix of maximum F-norm Corresponding alternative antenna is selected vector, ω kComputing formula as follows:Wherein, k=1 ..., B, Expression line options scheme on the same day is φ lThe time according to channel matrix H KkThe base station k that determines is to the channel matrix of this community user, φ lL alternative antenna of expression base station selected vector, Select the index variables of vector for alternative antenna, N fBe the rf chain way of base station, Expression is chosen N from the N number fThe number of combinations of number, || || FThe F-norm of representing matrix; The antenna selecting plan of each base station is combined into B cell-site antenna selection set omega={ ω 1, ω 2..., ω B;(3) along the optimization direction of system and capacity increase, in the situation that keep other antenna for base station selection schemes constant, utilize the part iteration to disturb alignment to determine successively the antenna for base station selection scheme of each residential quarter optimum:3.1) to select the index variables l of vector be 1 to the initialization alternative antenna, system and the capacity C of simultaneously initialization maximum MaxBe 0;3.2) with the antenna for base station selection scheme ω of residential quarter 1 1Be updated to alternative antenna and select vectorial φ l, keep the antenna selecting plan of other base stations constant, simultaneously the 1st element in B the cell-site antenna selection set omega is updated to φ l3.3) according to above-mentioned antenna selecting plan, carry out the part iteration and disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix, are implemented as follows:3.3a) pre-coding matrix of random initializtion base station k K=1 wherein ..., B, d represent that the base station sends the dimension of data flow;3.3b) initialization iterations indicator variable c is 1;3.3c) determine the interference covariance matrix of user k ' when the c time iteration K '=1 wherein ..., B;3.3d) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix The interference of design user k ' when the c time iteration suppresses matrix K '=1 wherein ..., B, υ dD characteristic value characteristic of correspondence vector of { } representing matrix minimum;3.3e) suppress matrix according to above-mentioned interference Determine the interference covariance matrix of base station k when the c time iteration K '=1 wherein ..., B, k=1 ..., B;3.3f) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix The pre-coding matrix of design base station k when the c time iteration K=1 wherein ..., B;3.3g) value of iterations indicator variable c adds 1, if c≤N Ite, then return step 3.3c, proceed iteration; Otherwise, jump to step 3.4, obtain N IteThe pre-coding matrix of base station k after the inferior iteration Suppress matrix with the interference of user k ' K=1 wherein ..., B, k '=1 ..., B, N IteDisturb the iterations of alignment for the part iteration;3.4) determine that base station k ' is to the equivalent channel of user k ' And the interference covariance matrix of user k ' Wherein Be illustrated in antenna selecting plan ω K 'Base station k ' arrives the channel matrix of user k ' in the situation, Expression N IteThe pre-coding matrix of base station k ' after the inferior iteration, k '=1 ..., B;3.5) equivalent channel of disturbing alignment to obtain according to the part iteration Disturb and suppress matrix And interference covariance matrix Computing system and capacity Wherein Be the capacity of user k ', k '=1 wherein ..., B;3.6) if C Ω>C Max, C then is set Max=C Ω, and the optimum alternative antenna of order is selected the index variables l of vector *=l, namely optimum alternative antenna is selected vector Otherwise, C MaxRemain unchanged;3.7) alternative antenna selects the value of the index variables l of vector to add 1, if Then return step 3.2, continue other alternative antenna of traversal and select vector; Otherwise, skip to step 3.8;3.9) according to the described method of step 3.1 to 3.8, determine that successively residential quarter 2 is to the antenna for base station selection scheme of residential quarter B optimum3.10) the antenna for base station selection scheme of above-mentioned each residential quarter optimum is combined into B optimum cell-site antenna selection set(4) adopt traditional iteration to disturb alignment to carry out many coordinating district interferences:4.1) select set omega according to B cell-site antenna of above-mentioned optimum *, carrying out iterations is N ITraditional iteration disturb alignment, the base station pre-coding matrix and the user that design each residential quarter disturb the inhibition matrix;4.2) disturb the pre-coding matrix of alignment design and interference to suppress matrix according to above-mentioned traditional iteration, each cell base station and user transmit and receive data respectively, carry out many coordinating district interferences, eliminate presence of intercell interference.
- 2. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, wherein the described base station k of step (1) sends pilot frequency information to user k, and be to carry out in the following manner: base station k divides The inferior user k of giving sends pilot frequency information, chooses each time some antennas and they are connected to radio frequency link to send pilot frequency information, wherein chooses N the 1st time from all N antenna fIndividual antenna, the 2nd time from remaining N-N fChoose N in the individual antenna fIndividual antenna, the 3rd time from remaining N-2N fChoose N in the individual antenna fIndividual antenna, until the Inferiorly choose last remaining antenna, thereby guarantee that all antennas all can once be connected to radio frequency link at certain and send pilot frequency information, k=1 wherein ..., B, B are the residential quarter number, N and N fRepresent antenna number and the rf chain way of base station, and N>N f, Representative rounds up to scalar x.
- 3. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, the channel matrix in the wherein said step (2) Determine with following formula:Wherein, k=1 ..., B, B are the residential quarter number, H KkBe the channel matrix of base station k to this community user, φ lL alternative antenna of expression base station selected vector, and it is the column vector that N grows, and the element value is whether 1 or 0 not represent respective antenna selected in this vector, For alternative antenna is selected vector index variable, N and N fRepresent antenna number and the rf chain way of base station, and N>N f, Expression is chosen N from the N number fThe number of combinations of number, Diag{ φ lRepresent that the element on the leading diagonal is column vector φ lThe diagonal matrix of corresponding element, function The submatrix that expression is made of the non-zero row of matrix.
- 4. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, the pre-coding matrix in the wherein said step (3.3a) Satisfy condition: K=1 wherein ..., B, B are the residential quarter number, Ω is that the B cell-site antenna is selected set, N fThe rf chain way of expression base station, d represents that the base station sends the dimension of data flow, () HThe conjugate transpose of representing matrix, I nThe unit matrix of expression n * n dimension, The expression complex field.
- In the multi-cell system according to claim 1 based on the antenna for base station system of selection of disturbing alignment, the interference covariance matrix of user k ' when the c time iteration in the wherein said step (3.3c) Its computing formula is as follows:Wherein, k '=1 ..., B, B are the residential quarter number, Ω is that B cell-site antenna selected set, P kThe transmitting power of expression base station k, d represents that the base station sends the dimension of data flow, The pre-coding matrix of base station k when being illustrated in the c-1 time iteration, Be illustrated in antenna selecting plan ω kBase station k is to the channel matrix of user k ', () in the situation HThe conjugate transpose of representing matrix.
- In the multi-cell system according to claim 1 based on the antenna for base station system of selection of disturbing alignment, the interference covariance matrix of base station k when the c time iteration in the wherein said step (3.3e) Its computing formula is as follows:Wherein, k=1 ..., B, B are the residential quarter number, Ω is that B cell-site antenna selected set, P K 'The transmitting power of expression user k ', value equate with the transmitting power of this cell base station, and d represents the dimension of base station transmission data flow, The interference of user k ' suppresses matrix when being illustrated in the c time iteration, Be illustrated in antenna selecting plan ω kBase station k is to the channel matrix of user k ', () in the situation HThe conjugate transpose of representing matrix.
- 7. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, the part iteration in the wherein said step (3.3g) is disturbed the iterations N of alignment IteNeed to satisfy N Ite≤ N I, N wherein IThe iterations empirical value that the expression tradition disturbs the alignment convergence to need is preferably 30.
- 8. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, the interference covariance matrix of the user k ' in the wherein said step (3.4) Its computing formula is as follows:Wherein, k '=1 ..., B, B are the residential quarter number, Ω is that B cell-site antenna selected set, P kThe transmitting power of expression base station k, d represents that the base station sends the dimension of data flow, Be illustrated in N IteThe pre-coding matrix of base station k after the inferior iteration, Be illustrated in antenna selecting plan ω kBase station k is to the channel matrix of user k ', () in the situation HThe conjugate transpose of representing matrix.
- 9. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, the capacity in the wherein said step (3.5) Its computing formula is as follows:Wherein, k '=1 ..., B, B are the residential quarter number, Ω is that B cell-site antenna selected set, P K 'The transmitting power of expression base station k ', d represents that the base station sends the dimension of data flow, With Represent respectively N IteThe interference of the pre-coding matrix of base station k ' and user k ' suppresses matrix after the inferior iteration, Be illustrated in antenna selecting plan φ K 'Base station k ' arrives the channel matrix of user k ' in the situation, The interference covariance matrix of expression user k ', () HThe conjugate transpose of representing matrix, σ 2The variance of white Gaussian noise in the channel, I nThe unit matrix of expression n * n dimension, M is user's antenna number, || the determinant of representing matrix, bodge are bit/s/Hz.
- 10. the antenna for base station system of selection of aliging based on interference in the multi-cell system according to claim 1, wherein described B the cell-site antenna according to optimum of step (4.1) selected set omega *, carrying out iterations is N ITraditional iteration disturb alignment, the base station pre-coding matrix and the user that design in each residential quarter disturb the inhibition matrix, carry out as follows:4.1a) pre-coding matrix of random initializtion base station k And satisfy K=1 wherein ..., B, B represent residential quarter number, N fRepresent the rf chain way of base station, d represents that the base station sends the dimension of data flow, () HThe conjugate transpose of representing matrix, I nThe unit matrix of expression n * n dimension, The expression complex field;4.1b) initialization iterations variable m is 1;4.1c) determine the interference covariance matrix of user k ' when the m time iterationWherein, k '=1 ..., B, P kThe transmitting power of expression base station k, The pre-coding matrix of base station k when being illustrated in the m-1 time iteration, Be illustrated in antenna selecting plan Base station k is to the channel matrix of user k ' in the situation;4.1d) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix The interference of design user k ' when the m time iteration suppresses matrix K '=1 wherein ..., B, d characteristic value characteristic of correspondence vector of υ d{} representing matrix minimum;4.1e) suppress matrix according to above-mentioned interference Determine the interference covariance matrix of base station k when the m time iterationWherein, k=1 ..., B, k '=1 ..., B, P K 'The transmitting power of expression user k ', value equates with the transmitting power of this cell base station;4.1f) reveal minimum principle according to disturbing, according to above-mentioned interference covariance matrix The pre-coding matrix of design base station k when the m time iteration K=1 wherein ..., B;4.1g) value of iterations variable m adds 1, if m≤N I, N wherein IThe iterations empirical value that the expression tradition disturbs the alignment convergence to need is preferably 30, then returns step 4.1c, proceeds iteration; Otherwise iterative process finishes, and obtains the pre-coding matrix of base station k Suppress matrix with the interference of user k ' K=1 wherein ..., B, k '=1 ..., B.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310025663.4A CN103067062B (en) | 2013-01-23 | 2013-01-23 | Base station antenna selecting method based on interference alignment in multi-cell system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310025663.4A CN103067062B (en) | 2013-01-23 | 2013-01-23 | Base station antenna selecting method based on interference alignment in multi-cell system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103067062A true CN103067062A (en) | 2013-04-24 |
CN103067062B CN103067062B (en) | 2015-03-04 |
Family
ID=48109543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310025663.4A Expired - Fee Related CN103067062B (en) | 2013-01-23 | 2013-01-23 | Base station antenna selecting method based on interference alignment in multi-cell system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103067062B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103209059A (en) * | 2013-04-26 | 2013-07-17 | 哈尔滨工业大学 | MIMO (Multiple Input Multiple Output) interference alignment wireless communication method under F-norm (Frobenius norm) bounded channel error model |
CN104683005A (en) * | 2013-12-02 | 2015-06-03 | 华为技术有限公司 | Interference suppression method and interference suppression device |
CN105007108A (en) * | 2015-07-03 | 2015-10-28 | 北京邮电大学 | Distributed interference alignment method based on transmit antenna selection |
CN105262526A (en) * | 2015-10-30 | 2016-01-20 | 哈尔滨工程大学 | Interference alignment method based on minimized interference projection matrix nuclear norm |
CN105978615A (en) * | 2016-05-17 | 2016-09-28 | 电子科技大学 | Multi-cell large-scale MIMO interference elimination method based on antenna selection and interference alignment |
CN106059638A (en) * | 2016-06-21 | 2016-10-26 | 重庆邮电大学 | Interference elimination and antenna optimization method in multi-cell large-scale MIMO system |
CN106603448A (en) * | 2016-12-23 | 2017-04-26 | 西安电子科技大学 | Interference alignment method based on multi-cell multi-user cooperation communication |
CN107104715A (en) * | 2017-06-07 | 2017-08-29 | 南京邮电大学 | A kind of interference alignment schemes based on day line options |
CN108429574A (en) * | 2018-01-24 | 2018-08-21 | 西安科技大学 | Extensive mimo system emitting antenna selecting method |
CN109714093A (en) * | 2018-11-05 | 2019-05-03 | 南京邮电大学 | A kind of joint antenna selection method towards isomery cellular network |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100227637A1 (en) * | 2009-03-03 | 2010-09-09 | Tae Soo Kwon | Signal transmission method and system for transmitting signal by using interference control method and/or transmission power control method |
CN102136857A (en) * | 2010-09-01 | 2011-07-27 | 华为技术有限公司 | Interference processing method for multi-antenna base station network, centralized controller and base station |
CN102710393A (en) * | 2012-05-25 | 2012-10-03 | 中国科学技术大学 | Interference alignment precoding method based on Stiefel manifold |
-
2013
- 2013-01-23 CN CN201310025663.4A patent/CN103067062B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100227637A1 (en) * | 2009-03-03 | 2010-09-09 | Tae Soo Kwon | Signal transmission method and system for transmitting signal by using interference control method and/or transmission power control method |
CN102136857A (en) * | 2010-09-01 | 2011-07-27 | 华为技术有限公司 | Interference processing method for multi-antenna base station network, centralized controller and base station |
CN102710393A (en) * | 2012-05-25 | 2012-10-03 | 中国科学技术大学 | Interference alignment precoding method based on Stiefel manifold |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103209059A (en) * | 2013-04-26 | 2013-07-17 | 哈尔滨工业大学 | MIMO (Multiple Input Multiple Output) interference alignment wireless communication method under F-norm (Frobenius norm) bounded channel error model |
US9887748B2 (en) | 2013-12-02 | 2018-02-06 | Huawei Technologies Co., Ltd. | Interference suppression method and apparatus |
CN104683005A (en) * | 2013-12-02 | 2015-06-03 | 华为技术有限公司 | Interference suppression method and interference suppression device |
WO2015081847A1 (en) * | 2013-12-02 | 2015-06-11 | 华为技术有限公司 | Interference suppressing method and apparatus |
CN104683005B (en) * | 2013-12-02 | 2018-11-30 | 华为技术有限公司 | Disturbance restraining method and device |
CN105007108A (en) * | 2015-07-03 | 2015-10-28 | 北京邮电大学 | Distributed interference alignment method based on transmit antenna selection |
CN105262526A (en) * | 2015-10-30 | 2016-01-20 | 哈尔滨工程大学 | Interference alignment method based on minimized interference projection matrix nuclear norm |
CN105262526B (en) * | 2015-10-30 | 2019-05-21 | 哈尔滨工程大学 | A kind of interference alignment schemes based on minimum interference projection matrix nuclear norm |
CN105978615B (en) * | 2016-05-17 | 2019-08-13 | 电子科技大学 | The extensive MIMO interference elimination method of multiple cell based on day line options and interference alignment |
CN105978615A (en) * | 2016-05-17 | 2016-09-28 | 电子科技大学 | Multi-cell large-scale MIMO interference elimination method based on antenna selection and interference alignment |
CN106059638A (en) * | 2016-06-21 | 2016-10-26 | 重庆邮电大学 | Interference elimination and antenna optimization method in multi-cell large-scale MIMO system |
CN106059638B (en) * | 2016-06-21 | 2019-12-31 | 重庆邮电大学 | Interference elimination and antenna optimization method in multi-cell large-scale MIMO system |
CN106603448A (en) * | 2016-12-23 | 2017-04-26 | 西安电子科技大学 | Interference alignment method based on multi-cell multi-user cooperation communication |
CN106603448B (en) * | 2016-12-23 | 2019-06-21 | 西安电子科技大学 | Interference alignment schemes based on the communication of multiple cell multi-user Cooperation |
CN107104715A (en) * | 2017-06-07 | 2017-08-29 | 南京邮电大学 | A kind of interference alignment schemes based on day line options |
CN107104715B (en) * | 2017-06-07 | 2020-06-12 | 南京邮电大学 | Interference alignment method based on antenna selection |
CN108429574A (en) * | 2018-01-24 | 2018-08-21 | 西安科技大学 | Extensive mimo system emitting antenna selecting method |
CN108429574B (en) * | 2018-01-24 | 2023-04-14 | 西安科技大学 | Method for selecting transmitting antenna of large-scale MIMO system |
CN109714093A (en) * | 2018-11-05 | 2019-05-03 | 南京邮电大学 | A kind of joint antenna selection method towards isomery cellular network |
CN109714093B (en) * | 2018-11-05 | 2022-04-05 | 南京邮电大学 | Combined antenna selection method for heterogeneous cellular network |
Also Published As
Publication number | Publication date |
---|---|
CN103067062B (en) | 2015-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103067062B (en) | Base station antenna selecting method based on interference alignment in multi-cell system | |
CN107453795B (en) | Beam allocation method of multi-user millimeter wave communication system, device and system thereof | |
CN103166685B (en) | A kind of interference alignment schemes based on joint Power distribution in LTE | |
CN103248461B (en) | Based on the multiple cell interference alignment iterative algorithm of beam forming | |
CN104868947A (en) | Method of realizing beam forming and base station | |
CN104702390A (en) | Pilot frequency distribution method in distributed compressive sensing (DCS) channel estimation | |
CN103840870A (en) | Method for lowering limiting feedback expenditure of 3D MIMO channel | |
CN101882952A (en) | Space division multiple address transmission method of statistical characteristic mode | |
CN105007108A (en) | Distributed interference alignment method based on transmit antenna selection | |
CN105450274A (en) | Optimal energy efficiency-based user number optimization method for large-scale and multi-antenna relay system | |
CN103814531B (en) | Improve the method with the transmitting gain at the network element of multiple antenna | |
CN107404375B (en) | Space-time cooperative transmission optimization method in irregular multipoint sea area communication network | |
CN106788631A (en) | A kind of extensive MIMO reciprocities calibration method based on local alignment | |
CN112564746B (en) | Optimal GEE-based power distribution method in CF mmWave mMIMO system | |
CN105207723B (en) | channel correcting method, base station, user equipment and communication system | |
CN104901732A (en) | Pilot frequency multiplexing method in dense node configuration system | |
Almasaoodi et al. | New Quantum Strategy for MIMO System Optimization. | |
CN102882576A (en) | Multi-point coordinated transmission method in LTE-A (long term evolution-advance) downlink system | |
CN1953348B (en) | A estimation method and device for the direction of incoming wave of intelligent multiple sector antenna system | |
CN101989875B (en) | Multi-cell interference suppression method and base station controller | |
CN101321008A (en) | Descending beam forming emission method and device | |
CN111277313B (en) | Bipartite graph-based large-scale MIMO beam selection and transmission method for cellular internet of vehicles | |
CN101953088B (en) | Orthogonal projection precoding and decoding method and equipment by using part CSI feedback | |
Gao et al. | Making wideband channel estimation feasible for mmwave massive MIMO: A doubly sparse approach | |
CN103354658B (en) | Multiple cell multi-input multi-output system considers the optimized for energy efficiency power distribution algorithm that co-channel interference suppresses |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150304 Termination date: 20200123 |
|
CF01 | Termination of patent right due to non-payment of annual fee |