CN107181519A - A kind of null based on mobile target DOA extends 3D MIMO beam form-endowing methods - Google Patents
A kind of null based on mobile target DOA extends 3D MIMO beam form-endowing methods Download PDFInfo
- Publication number
- CN107181519A CN107181519A CN201710556742.6A CN201710556742A CN107181519A CN 107181519 A CN107181519 A CN 107181519A CN 201710556742 A CN201710556742 A CN 201710556742A CN 107181519 A CN107181519 A CN 107181519A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- array
- array factor
- doa
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Abstract
3D MIMO beam form-endowing methods are extended the invention discloses a kind of null based on mobile target DOA, belong to AF panel field.First on the basis of existing 3D MIMO wave beam formings, for base station and certain mobile terminal of communication, the antenna array of M × N number of array element is constructed on x/y plane;And the weights of each array factor are calculated respectively;Determine the DOA information between destination mobile terminal and interfering mobile terminal and base station;Then the weights of each array factor are readjusted according to DOA information;Constrained finally by applying to desired orientation and the peripherad angular regions of unexpected side, null extension is carried out on unexpected direction, while undistorted in the desired direction.The problem of interference rejection capability declines when solving null narrower width under user mobility scene, improve the robustness of system, so as to effectively suppress the interference of change in location, and make desired signal obtain, preferably close to undistorted corresponding, effectively increasing spectrum efficiency to meet higher rate requirement.
Description
Technical field
The invention belongs to AF panel field, specifically a kind of null based on mobile target DOA extends 3D-MIMO wave beams
Shaping method.
Background technology
According to《5G visions and demand white paper》Description, the wide area covering scene that 5G is supported will be provided the user up to
10Mbit/s experience density, supports the data traffic density of hot spot region to experience speed and 10Tbit/km2 up to 10Gbit/s
Packing density.But, not only quantity is few for traditional mimo antenna, and control ability is poor, and its spatial reuse, space diversity and ripple
The low performance of the functions such as beam figuration and the application of stressing property, can not also receive for 5G.
3D-MIMO technologies are one of 5G key technologies, and the sky of vertical direction is added on the basis of 2D-MIMO technologies
Between the free degree, it is possible to increase power system capacity, increase spectrum efficiency, effectively reduce inter-cell interference.From LTE Rel-11 standardization
Starting stage just there is company to start to promote 3D-MIMO beamforming techniques, and had more next in Rel-12 preparatory stage
More operators and enterprise of equipment manufacturers show the enthusiasm to 3D-MIMO beamforming techniques.
The operation principle of smart antenna mainly includes two processes, is antenna system first to launching from mobile terminal
The direction of arrival (DOA) of multipath signal is estimated, calculates and is carried out behind distance, angle of declination and the azimuth of mobile terminal and base station
Space filtering, so as to suppress the interference of other terminal-pair base stations;Next to that base station is according to DOA information, adjust on each oscillator of antenna
The amplitude of signal and the weights of phase, carry out digital beam forming to the signal that smart antenna is sent, make Base Transmitter signal
Main lobe can be with less field angle and higher power density, and the direction of arrival for moving along terminal electric wave signal sends movement back to
Terminal, so that the main lobe direction of the electromagnetic wave of smart antenna transmitting is to quasiexpectation user, zero valve direction alignment interference source.
Obviously, the core of smart antenna is the generation of wave beam forming and the positioning of shaped-beam, should produce power density
The farther shaped-beam of high, radiation length, will also enable shaped-beam to be accurately positioned user terminal.
The wave beam forming of mimo antenna is a kind of Signal Pretreatment technology based on aerial array, and it is by adjusting antenna
The weight coefficient of each array element produces the wave beam with directive property in array, so as to obtain obvious array gain.If can basis
The weight coefficient of each array element of channel condition suitable control, can equally be obtained while desired orientation signal intensity is strengthened,
The interference to unexpected direction is reduced as far as possible.If weight coefficient is obtained by the DOA parameters pretreatment of terminal, wave beam
Figuration can just point to terminal user.But existing beamforming technique is all based on what signal during user's inactive state was estimated
DOA information is adjusted, and in user mobility scene, wave beam forming performance will be reduced seriously.
3D-MIMO beam-forming systems need a sets of plan to determine weight, so as to adjust beam position so that wave beam refers to
To ideal user terminal, and form null to mitigate interference in the direction of non-ideal user.Signal to Interference plus Noise Ratio can so be increased, and
And mitigated the interference from other secondary lobes.Therefore a set of weight determining method is needed, it is horizontal and vertical to build respectively
The weights of array, minimize the mean square error of array pattern vector sum unit vector, and unit vector represents ideal orientation here
Array pattern, null vector represents the array pattern in non-ideal direction.Because rectangle plane array can regard the linear array of M N-dimensional as
Row, the weight that the weight that such M-N ties up array can be respectively from horizontal array and orthogonal array is obtained.
Beam steering and null method rely primarily on angle domain information and carry out shape-endowing weight value generation, when user moves
Wait, because the uncertainty of user terminal location may influence the accuracy that direction of arrival DOA estimates, so as to reduce wave beam forming
The AF panel performance of technology.The null of traditional beam form-endowing method formation is narrower, when signal Mutual coupling is accurate
Interference signal can be suppressed well.But in practice due to signal DOA change, can cause interference with to null positional misalignment with
And desired signal is oriented to appropriate mismatch so that contain not repressed interference signal in the signal of wave beam forming output.It is oriented to arrow
Amount mismatch can cause wave beam forming hydraulic performance decline.
The content of the invention
The present invention is on the basis of existing 3D-MIMO wave beam formings, for null narrower width under user mobility scene
When interference rejection capability decline, in order to constrain null width, it is more stable that side lobe gain and weight vector modulus value etc. obtain performance
Wave beam forming, it is proposed that a kind of null based on mobile target DOA extends 3D-MIMO beam form-endowing methods.
Comprise the following steps that:
Step 1: base station and certain mobile terminal for communication, construct the antenna array of M × N number of array element on x/y plane;
Level ground is chosen as x-axis, and x-axis is provided with M array element, the level ground vertical with x-axis is y-axis, and y-axis
It is provided with N number of array element;
Step 2: calculating the weights of each array factor in antenna array respectively;
M-n array factor weight wmnIt is calculated as follows:
wmn=am·bn
amRefer to m-n weights of the array factor in x-axis direction;bnRefer to m-n weights of the array factor in y-axis direction.
Step 3: according to the position of destination mobile terminal and interfering mobile terminal, determining respectively between each terminal and base station
DOA information;
Destination mobile terminal refers to the mobile terminal with base station communication;
DOA between destination mobile terminal and base station is set toDOA between interfering mobile terminal and base station is set toI=1,2,3 ...;Because the angle estimation deviation that the mobility of terminal is produced is set to
All DOA angles play initial line by angle of x-axis;
Step 4: according to the DOA information between each terminal and base station, readjusting the weights of each array factor;
It is specific as follows:
Step 401, by the every value AF of row array factor in the horizontal direction in M × N number of array elementxWith each column array factor in Vertical Square
To value AFy, utilization orientation figure product principle obtains the expression formula of each array factor;
dxFor the spacing between two neighboring array factor in x-axis, dyFor the spacing between two neighboring array factor in y-axis.
K is regulation coefficient, is integer;βxRepresent the phase delay in horizontal direction, βyRepresent that the phase in vertical direction is prolonged
Late, (θ, φ) is the angle value between each terminal and base station;
Step 402, the DOA information brought between each terminal and base station, the corresponding array factor of destination mobile terminal is expressed
The value of formula is set to 1, and the value of the corresponding array factor expression formula of interfering mobile terminal is set into 0;
As a result it is as follows:
Step 403, the value according to each array factor expression formula, utilize the value AF of often row array factor in the horizontal directionxCalculate
The weight W of horizontal directionx;
First, the value AF according to every row array factor in the horizontal directionxFormula, calculates weight WxExpression formula;
Wherein, Wx=[a1 a2…aM]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
Finally, to weight WxExpression formula is solved;
IfB=[1 0 ... 0]T, above formula abbreviation is
AWx=B;
Wx=A-1B, A are M × M invertible matrix;Wx=pinv (A) * B, pinv (A) are Moore-Penrose generalized inverses.
Step 404, the value according to each array factor expression formula, using each column array factor vertical direction value AFyCalculate
The weight W of vertical directiony;
First, according to each column array factor vertical direction value AFyFormula, calculates weight WyExpression formula;
Wherein, Wy=[b1 b2…bN]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
Finally, to weight WyExpression formula is solved;
IfB'=[1 0 ... 0]T, above formula abbreviation is A'Wy
=B';
Wy=A'-1B', A' are N × N invertible matrix;Wy=pinv (A') * B', pinv (A') are Moore-Penrose
Generalized inverse.
Step 405, the weight W using horizontal directionxWith the weight W of vertical directionyThe weight matrix W of array antenna is calculated, is entered
One step obtains the weights after each array factor adjustment;
W=Wx(Wy)T=(wmn)M×N
Step 5: with the weights of each array factor after adjustment, by desired orientation and the peripherad angle of unexpected side
Spend region and apply constraint, null extension is carried out on unexpected direction, while undistorted in the desired direction;
Null extension includes two parts:DOA signals θ in desired orientation0Expand to [θ0l,θ0h] and unexpected direction
DOA signals θiExpand to [θil,θih]。
Specially:To extend the DOA signals θ in unexpected directioni, work as θi∈[θil,θih] when array meet:|wHs(θi)|≤
δ, then θiExpand to [θil,θih];δ is null response constraint threshold value;
For the DOA signals θ in extension desired orientation0, as θ ∈ [θ0l,θ0h] when θ0Close to undistorted, and meet wHs(θ)-1
≤ μ, then θ0Expand to [θ0l,θ0h];μ is a minimum.
The advantage of the invention is that:
1), a kind of null based on mobile target DOA extends 3D-MIMO beam form-endowing methods, compared to traditional wave beam
Figuration scheme, the present invention effectively increases spectrum efficiency to meet higher rate requirement in terms of wireless technology.
2), a kind of null based on mobile target DOA extends 3D-MIMO beam form-endowing methods, in moving target, works as letter
Number estimation is not punctual, by effective null-broadening, extends effective range of receiving to signal, and applies to expecting and interference radiating way
Array response constraint, and the mode such as Sidelobe control constraint, to moving target and when there is larger signal DOA estimation deviation,
The robustness of raising system.
Brief description of the drawings
Fig. 1 is end-user mobility scene graph in 3D-MIMO of the present invention;
Fig. 2 is that null of the present invention based on mobile target DOA extends 3D-MIMO beam form-endowing method flow charts;
Fig. 3 is active face array factor geometry distribution map of the present invention;
Fig. 4 is intelligent antenna beam shaping 3D antenna radiation patterns of the present invention.
Embodiment
Below in conjunction with drawings and examples, the present invention is described in further detail.
The present invention relates to terminal in 3D-MIMO (3-Dimensional Multiple-Input Multiple-Output)
User mobility scene, as shown in figure 1, θ0Represent mobile terminal to BS initial level angular separation, θ1Represent that movement is later
Horizontal direction angle;Δ θ represents angle difference;d0Represent base station and mobile terminal initial distance;d1After representing that mobile terminal is moved
With the distance of base station.t0Represent initial time, t1The mobile later time is represented, Δ t represents the time difference.
For the base station BS and certain mobile terminal MES of communication, when the target terminal enters interference region, incoming signal
DOA change causes wave beam forming hydraulic performance decline when there is evaluated error, forms wider at interference by null broadening
Null, so as to effectively suppress the interference of change in location, and makes desired signal obtain preferably close to undistorted corresponding.
As shown in Fig. 2 comprising the following steps that:
Step 1: base station and certain mobile terminal for communication, construct the antenna array of M × N number of array element on x/y plane;
To the rectangle plane battle array in an x-y plane, as shown in figure 3, having in x-axis has N number of array element in M array element, y-axis,
Constitute the antenna array of M × N number of array element;Rectangle plane battle array can be considered M line array of N number of array element, or M array element it is N number of directly
Linear array.The horizontal azimuth of array antenna is represented, θ represents the vertical direction elevation angle, dxAnd dyThe level of antenna element is represented respectively
With vertical direction spacing;U represents to receive signal.
Step 2: calculating the weights of each array factor in antenna array respectively;
The weights of m-n array elements are wmnIt is calculated as follows:
wmn=am·bn
amRefer to m-n weights of the array factor in x-axis direction;bnRefer to m-n weights of the array factor in y-axis direction.
Step 3: according to the position of destination mobile terminal and interfering mobile terminal, calculating respectively between each terminal and base station
DOA information;
Destination mobile terminal refers to the mobile terminal with base station communication;
DOA in the present embodiment between destination mobile terminal and base station is set toBetween interfering mobile terminal and base station
DOA be set toWithBecause the angle estimation deviation that the mobility of terminal is produced is set to
All DOA angles play initial line by angle of x-axis;
Step 4: according to the DOA information between each terminal and base station, readjusting the weights of each array factor;
It is specific as follows:
Step 401, by the every value AF of row array factor in the horizontal direction in M × N number of array elementxWith each column array factor in Vertical Square
To value AFy, utilization orientation figure product principle obtains the expression formula of each array factor;
The array factor of the antenna array independent role of M array element or N number of array element is had learned that, so being obtained by directional diagram product principle
To the Antenna Array Pattern of whole M × N number of array element, obtained by directional diagram product principle:
dxFor the spacing between two neighboring array factor in x-axis, dyFor the spacing between two neighboring array factor in y-axis.
K is regulation coefficient, is integer;βxRepresent the phase delay in horizontal direction, βyRepresent that the phase in vertical direction is prolonged
Late, when needing wave beam phase modulation,Wherein θ0WithRepresent phase modulation to θ0WithAngle.(θ, φ) is the angle value between each terminal and base station;
Step 402, the DOA information brought between each terminal and base station, the corresponding array factor of destination mobile terminal is expressed
The value of formula is set to 1, and the value of the corresponding array factor expression formula of interfering mobile terminal is set into 0;
In order to receive the signal of desired orientation, and it is minimized the influence in unexpected direction.Need to make expectation
Direction produces unit gain so that AF=1;And it is upwardly formed zero point in unexpected side so that AF=0.I.e. in the desired direction
Have:AFx=1and AFy=1;There is AF on unexpected directionx=0and AFy=0.
In order to receive the signal of desired orientation, and it is minimized the influence in unexpected direction;Need to make expectation
Direction produces unit gain, i.e.,And it is upwardly formed zero point i.e. in unexpected side
Step 403, the value according to each array factor expression formula, utilize the value AF of often row array factor in the horizontal directionxCalculate
The weight W of horizontal directionx;
First, the value AF according to every row array factor in the horizontal directionxFormula, calculates weight WxExpression formula;
Wherein, Wx=[a1 a2…aM]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
There are four wave beams in a region in the present invention, there is θ=θ in the desired direction0,The other three is unexpected
Direction has respectively
On unexpected direction, it can obtain:
Finally, to weight WxExpression formula is solved;
IfB=[1 0 ... 0]T, above formula abbreviation is AWx=B;
Wx=A-1B, A are 4 × M invertible matrix;The weight W of level dimensionx=pinv (A) * B, pinv (A) are Moore-
Penrose generalized inverses.
Step 404, the value according to each array factor expression formula, using each column array factor vertical direction value AFyCalculate
The weight W of vertical directiony;
First, according to each column array factor vertical direction value AFyFormula, calculates weight WyExpression formula;
Wherein, Wy=[b1 b2…bN]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
Finally, to weight WyExpression formula is solved;
IfB'=[1 0 ... 0]T, above formula abbreviation is A'Wy
=B';
Wy=A'-1B', A' are N × N invertible matrix;Wy=pinv (A') * B', pinv (A') are Moore-Penrose
Generalized inverse.
Step 405, the weight W using horizontal directionxWith the weight W of vertical directionyThe weight matrix W of array antenna is calculated, is entered
One step obtains the weights after each array factor adjustment;
W=Wx(Wy)T=(wmn)M×N
Because the mobility of user may influence the estimation in angle of arrival:
θ1=θ0+Δθ
Wherein, θtmp=θ0+(π-θv), Δ d=v × Δ t, 0≤Δ t < Tud。
Step 5: with the weights of each array factor after adjustment, by desired orientation and the peripherad angle of unexpected side
Spend region and apply constraint, null extension is carried out on unexpected direction, while undistorted in the desired direction;
The null of traditional Beamforming Method formation is narrower, can press down well when signal Mutual coupling is accurate
Interference signal processed.But it can cause interference with and be lost to null position due to reasons such as the motor message DOA changes of user in practice
Accurate and desired signal is oriented to mismatch;When null expansion technique can effectively solve interference signal motion, under wave beam forming performance
The problem of drop;
Null extension includes two parts:DOA signals θ in desired orientation0Expand to [θ0l,θ0h] and unexpected direction
DOA signals θiExpand to [θil,θih]。
Specially:Interference signal is moved
Assuming that desired signal DOA is θ0, interference signal DOA is in [θil,θih] interior change, in order to preferably receive desired signal
And suppressing interference signal, array accordingly needs to meet:
H(θ0)=wHs(θ0)=1
H(θi)=| wHs(θi)|≤δ
θi∈[θil,θih], δ is null response constraint threshold value;
Beam pattern side lobe gain is controlled to be less than a certain threshold value η, then the array response in secondary lobe region needs to meet
|wHs(θΘ)|≤η
Wherein, Θ represents the angular range of secondary lobe region overlay.
Desired signal is moved
When desired signal is in certain limit [θ0l,θ0h] in change when, to obtain preferable expected response, then need meet
In [θ0l,θ0h] in the range of close to undistorted, i.e.,
wHs(θ)-1≤μ
Wherein, θ ∈ [θ0l,θ0h], μ is a minimum.
Discrete angular information is equally spaced chosen in desired signal DOA angular ranges, then
wHs(θ0t) -1≤μ, t=1,2 ..., T
Wherein, T is the number of selected discrete angular;Assuming that K fixation carrys out self-interference signal respectively from θik, k=1,
2 ..., K direction are incident, then are in the constraints of interference radiating way
H(θik)≤δ, k=1,2 ..., K
It is broadly divided into two steps:
Suitable pre-coding matrix is selected according to DOE information, and forms null in this direction.
Receiving signal can be write as
Y=HPx
Wherein, H is channel vector, and P is pre-coding matrix
Select DFT code books as candidate wave beam formed matrix, the precoding based on DFT code books in linear array have compared with
Good beam direction, and can be with design processes simplified, pre-coding matrix can be write as
G represents codebook size, P(g)(m, n) represents the m rows n row of pre-coding matrix.Represent g
Individual pre-coding matrix.
Assuming that linear array shows NtIndividual antenna element, spacing is half-wavelengthThe corresponding factor of array can be write as
Order
Assuming thatForM rows, such zero point can obtain by following formula
Obtained by pre-coding matrix.
The null based on mobile target DOA that the present invention is ultimately formed extends intelligent antenna beam shaping 3D antenna directions
Figure, as shown in Figure 4.
Claims (4)
1. a kind of null based on mobile target DOA extends 3D-MIMO beam form-endowing methods, it is characterised in that specific steps are such as
Under:
Step 1: base station and certain mobile terminal for communication, construct the antenna array of M × N number of array element on x/y plane;
Level ground is chosen as x-axis, and x-axis, provided with M array element, the level ground vertical with x-axis is to be set in y-axis, and y-axis
There is N number of array element;
Step 2: calculating the weights of each array factor in antenna array respectively;
M-n array factor weight wmnIt is calculated as follows:
wmn=am·bn
amRefer to m-n weights of the array factor in x-axis direction;bnRefer to m-n weights of the array factor in y-axis direction;
Step 3: according to the position of destination mobile terminal and interfering mobile terminal, the DOA between each terminal and base station is determined respectively
Information;
DOA between destination mobile terminal and base station is set toDOA between interfering mobile terminal and base station is set toBecause the angle estimation deviation that the mobility of terminal is produced is set toAll DOA angles are equal
Initial line is played by angle of x-axis;
Step 4: according to the DOA information between each terminal and base station, readjusting the weights of each array factor;
Implement step as follows:
Step 401, by the every value AF of row array factor in the horizontal direction in M × N number of array elementxWith each column array factor in vertical direction
Value AFy, utilization orientation figure product principle obtains each array factor AF expression formula;
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mi>F</mi>
<mo>=</mo>
<msub>
<mi>AF</mi>
<mi>x</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>AF</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>m</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>kd</mi>
<mi>x</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>kd</mi>
<mi>y</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>kd</mi>
<mi>x</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>+</mo>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>kd</mi>
<mi>y</mi>
</msub>
<mi>sin</mi>
<mi>&theta;</mi>
<mi>cos</mi>
<mi>&phi;</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</msup>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Wherein, dxFor the spacing between two neighboring array factor in x-axis, dyFor the spacing between two neighboring array factor in y-axis;k
For regulation coefficient, positive integer is taken;βxRepresent the phase delay in horizontal direction, βyThe phase delay in vertical direction is represented, (θ,
It is φ) angle value between each terminal and base station;
Step 402, the DOA information brought between each terminal and base station, by the corresponding array factor expression formula of destination mobile terminal
Value is set to 1, and the value of the corresponding array factor expression formula of interfering mobile terminal is set into 0;
As a result it is as follows:
Step 403, the value according to each array factor expression formula, utilize the value AF of often row array factor in the horizontal directionxCalculated level side
To weight Wx;
Step 404, the value according to each array factor expression formula, using each column array factor vertical direction value AFyCalculate Vertical Square
To weight Wy;
Step 405, the weight W using horizontal directionxWith the weight W of vertical directionyThe weight matrix W of array antenna is calculated, further
Obtain the weights after each array factor adjustment;
W=Wx(Wy)T=(wmn)M×N
Step 5: with the weights of each array factor after adjustment, by desired orientation and the peripherad angular area of unexpected side
Domain applies constraint, null extension is carried out on unexpected direction, while undistorted in the desired direction.
2. a kind of null based on mobile target DOA as claimed in claim 1 extends 3D-MIMO beam form-endowing methods, it is special
Levy and be, described step 403 is specific as follows:
First, the value AF according to every row array factor in the horizontal directionxFormula, calculates weight WxExpression formula;
Wherein, Wx=[a1 a2 … aM]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
Finally, to weight WxExpression formula is solved;
IfB=[1 0 ... 0]T, above formula abbreviation is AWx
=B;
Wx=A-1B, A are M × M invertible matrix;Wx=pinv (A) * B, pinv (A) are Moore-Penrose generalized inverses.
3. a kind of null based on mobile target DOA as claimed in claim 1 extends 3D-MIMO beam form-endowing methods, it is special
Levy and be, described step 404 is specific as follows:
First, according to each column array factor vertical direction value AFyFormula, calculates weight WyExpression formula;
Wherein, Wy=[b1 b2 … bN]T;
Then, in the desired direction, it is 1 according to the value of the corresponding array factor expression formula of destination mobile terminal, obtains:
On unexpected direction, it is 0 according to the value of the corresponding array factor expression formula of interfering mobile terminal, obtains:
Finally, to weight WyExpression formula is solved;
IfB'=[1 0 ... 0]T, above formula abbreviation is A'Wy=
B';
Wy=A'-1B', A' are N × N invertible matrix;Wy=pinv (A') * B', pinv (A') are Moore-Penrose broad sense
It is inverse.
4. a kind of null based on mobile target DOA as claimed in claim 1 extends 3D-MIMO beam form-endowing methods, it is special
Levy and be, null extension includes two parts in described step five:DOA signals θ in desired orientation0Expand to [θ0l,θ0h] with
And the DOA signals θ in unexpected directioniExpand to [θil,θih];Specially:
To extend the DOA signals θ in unexpected directioni, work as θi∈[θil,θih] when array meet:|wHs(θi) |≤δ, then θiExtension
For [θil,θih];δ is null response constraint threshold value;
For the DOA signals θ in extension desired orientation0, as θ ∈ [θ0l,θ0h] when θ0Close to undistorted, and meet wHS (θ) -1≤μ,
Then θ0Expand to [θ0l,θ0h];μ is a minimum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710556742.6A CN107181519B (en) | 2017-07-10 | 2017-07-10 | Null steering expansion 3D-MIMO beam forming method based on moving target DOA |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710556742.6A CN107181519B (en) | 2017-07-10 | 2017-07-10 | Null steering expansion 3D-MIMO beam forming method based on moving target DOA |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107181519A true CN107181519A (en) | 2017-09-19 |
CN107181519B CN107181519B (en) | 2020-08-28 |
Family
ID=59845158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710556742.6A Active CN107181519B (en) | 2017-07-10 | 2017-07-10 | Null steering expansion 3D-MIMO beam forming method based on moving target DOA |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107181519B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109039405A (en) * | 2018-07-04 | 2018-12-18 | 西安电子科技大学 | A kind of skip communication means based on interference space rotation |
CN109039413A (en) * | 2018-08-01 | 2018-12-18 | 北京航空航天大学 | It is a kind of applied to the 3D flexibility covering method for facing empty communication |
CN109120322A (en) * | 2018-06-06 | 2019-01-01 | 西安科技大学 | Multiple cell multi-user's 3D-MIMO beam form-endowing method |
CN109444809A (en) * | 2018-10-30 | 2019-03-08 | 天津津航计算技术研究所 | A kind of UAV TT & C's method based on smart antenna |
CN113543144A (en) * | 2020-05-26 | 2021-10-22 | 中兴通讯股份有限公司 | Wireless communication method, terminal, base station, system, electronic device, and medium |
WO2023284571A1 (en) * | 2021-07-14 | 2023-01-19 | 中兴通讯股份有限公司 | Beamforming method and apparatus, base station, and storage medium |
WO2023124959A1 (en) * | 2021-12-31 | 2023-07-06 | 华为技术有限公司 | Interference avoidance method, apparatus, and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080170554A1 (en) * | 2005-03-02 | 2008-07-17 | Zte Corporation | Method and Equipment for Realizing Smart Antenna in Wcdma System |
CN101425832A (en) * | 2008-11-10 | 2009-05-06 | 中兴通讯股份有限公司 | Adaptive multiuser wave beam forming method based on wide nulling |
US7570211B1 (en) * | 2008-03-25 | 2009-08-04 | Rockwell Collins, Inc. | Digital beamforming method and apparatus for pointing and null steering without calibration or calculation of covariance matrix |
CN101924586A (en) * | 2009-06-10 | 2010-12-22 | 中兴通讯股份有限公司 | Multi-user beamforming method and base station |
-
2017
- 2017-07-10 CN CN201710556742.6A patent/CN107181519B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080170554A1 (en) * | 2005-03-02 | 2008-07-17 | Zte Corporation | Method and Equipment for Realizing Smart Antenna in Wcdma System |
US7570211B1 (en) * | 2008-03-25 | 2009-08-04 | Rockwell Collins, Inc. | Digital beamforming method and apparatus for pointing and null steering without calibration or calculation of covariance matrix |
CN101425832A (en) * | 2008-11-10 | 2009-05-06 | 中兴通讯股份有限公司 | Adaptive multiuser wave beam forming method based on wide nulling |
CN101924586A (en) * | 2009-06-10 | 2010-12-22 | 中兴通讯股份有限公司 | Multi-user beamforming method and base station |
Non-Patent Citations (2)
Title |
---|
EHAB ALI 等: "Beamforming techniques for massive MIMO systems in 5G: overview, classification, and trends for future research", 《FRONTIERS OF INFORMATION TECHNOLOGY & ELECTRONIC ENGINEERING》 * |
徐佳祥: "3D MIMO场景下的波束赋形技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑(月刊 )》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109120322A (en) * | 2018-06-06 | 2019-01-01 | 西安科技大学 | Multiple cell multi-user's 3D-MIMO beam form-endowing method |
CN109039405A (en) * | 2018-07-04 | 2018-12-18 | 西安电子科技大学 | A kind of skip communication means based on interference space rotation |
CN109039405B (en) * | 2018-07-04 | 2021-03-23 | 西安电子科技大学 | Air-skipping communication method based on interference subspace rotation |
CN109039413A (en) * | 2018-08-01 | 2018-12-18 | 北京航空航天大学 | It is a kind of applied to the 3D flexibility covering method for facing empty communication |
CN109039413B (en) * | 2018-08-01 | 2020-07-03 | 北京航空航天大学 | 3D flexible covering method applied to air-to-air communication |
CN109444809A (en) * | 2018-10-30 | 2019-03-08 | 天津津航计算技术研究所 | A kind of UAV TT & C's method based on smart antenna |
CN109444809B (en) * | 2018-10-30 | 2023-02-03 | 天津津航计算技术研究所 | Unmanned aerial vehicle measurement and control method based on intelligent antenna |
CN113543144A (en) * | 2020-05-26 | 2021-10-22 | 中兴通讯股份有限公司 | Wireless communication method, terminal, base station, system, electronic device, and medium |
WO2021238323A1 (en) * | 2020-05-26 | 2021-12-02 | 中兴通讯股份有限公司 | Wireless communication method, terminal, base station, system, electronic device and medium |
CN113543144B (en) * | 2020-05-26 | 2023-01-10 | 中兴通讯股份有限公司 | Wireless communication method, terminal, base station, system, electronic device, and medium |
WO2023284571A1 (en) * | 2021-07-14 | 2023-01-19 | 中兴通讯股份有限公司 | Beamforming method and apparatus, base station, and storage medium |
WO2023124959A1 (en) * | 2021-12-31 | 2023-07-06 | 华为技术有限公司 | Interference avoidance method, apparatus, and system |
Also Published As
Publication number | Publication date |
---|---|
CN107181519B (en) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107181519A (en) | A kind of null based on mobile target DOA extends 3D MIMO beam form-endowing methods | |
CN109347525B (en) | Adaptive multi-beam forming and beam switching method for millimeter-wave high-speed rail train-ground communication | |
US20170332249A1 (en) | Methods and Apparatus for Generating Beam Pattern with Wider Beam Width in Phased Antenna Array | |
JP5345690B2 (en) | Method for multi-antenna signal processing in antenna element arrangement, corresponding transceiver and corresponding antenna element arrangement | |
RU2515503C2 (en) | Communication method and system for directional transmission of digital data between aircraft and ground station | |
Luo et al. | Efficient multiple-group multiple-antenna (MGMA) scheme for high-speed railway viaducts | |
CN109444809B (en) | Unmanned aerial vehicle measurement and control method based on intelligent antenna | |
CN101364828A (en) | Downlink beam forming method | |
CN101848021B (en) | Method and device for generating broadcast beam weight of intelligent antenna array | |
CN106031068A (en) | Beamforming-based communication method and apparatus | |
CN112039567B (en) | Beam forming design method of multi-intelligent reflector system | |
US20160127920A1 (en) | Methods and Apparatus for Antenna Elevation Design | |
CN104320169B (en) | Three-dimensional wave beam Shape design method in multi-user's 3D mimo systems | |
CN101924586A (en) | Multi-user beamforming method and base station | |
CN106034307A (en) | Method and device for estimating direction of arrival (DOA) beam forming (BF) weight. | |
CN102457319A (en) | Business data distribution method and device | |
Alluhaibi et al. | Capacity maximisation for hybrid digital-to-analog beamforming mm-wave systems | |
Su et al. | Wideband precoding for RIS-aided THz communications | |
CN105356924A (en) | Method and device for estimating DOA-BF weight of dual-polarized antenna system | |
CN116436502A (en) | Intelligent reflecting surface selection and phase matrix adjustment method | |
CN107171710A (en) | beam-forming method and system | |
CN103248412B (en) | A kind of method of satellite navigation Multibeam synthesis | |
CN106027131A (en) | Base station and beam forming method | |
CN114222310B (en) | Method for optimizing reflection of combined 3D wave beam forming and intelligent reflecting surface | |
CN105510872A (en) | Two-dimensional arrival direction measuring method suitable for large-scale MIMO system |
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 |