CN1262130C - Method of estimating space reaching direction of fixed wave beam - Google Patents

Method of estimating space reaching direction of fixed wave beam Download PDF

Info

Publication number
CN1262130C
CN1262130C CN 03121293 CN03121293A CN1262130C CN 1262130 C CN1262130 C CN 1262130C CN 03121293 CN03121293 CN 03121293 CN 03121293 A CN03121293 A CN 03121293A CN 1262130 C CN1262130 C CN 1262130C
Authority
CN
China
Prior art keywords
search
angle
user
estimating
signal
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.)
Expired - Lifetime
Application number
CN 03121293
Other languages
Chinese (zh)
Other versions
CN1535048A (en
Inventor
李峰
邹素玲
康绍莉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
Datang Mobile Communications Equipment Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Datang Mobile Communications Equipment Co Ltd filed Critical Datang Mobile Communications Equipment Co Ltd
Priority to CN 03121293 priority Critical patent/CN1262130C/en
Priority to PCT/CN2004/000296 priority patent/WO2004088347A1/en
Publication of CN1535048A publication Critical patent/CN1535048A/en
Application granted granted Critical
Publication of CN1262130C publication Critical patent/CN1262130C/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • G01S3/38—Systems for determining direction or deviation from predetermined direction using adjustment of real or effective orientation of directivity characteristic of an antenna or an antenna system to give a desired condition of signal derived from that antenna or antenna system, e.g. to give a maximum or minimum signal

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The present invention relates to a method for estimating direction of arrival (DOA) of a fixing beam in space, which forms a gradually hierarchical thinning estimation method for calculating DOA by single-stage search substituted by multi-stage search. The method comprises that by adopting a classical beam forming algorithm, a user carries out first-stage search in the entire beam space range of an aerial array and total searches L1 times, the maximum value is taken from output power obtained by each search so as to obtain that a user arrival angle corresponding to the maximum value is a first search angle; then, in the appointed space range of the first search angle, the user carries out a second-stage search by adopting the classical beam forming algorithm and total searches L2 times, the maximum value is taken from the output power obtained by each search so as to obtain that a careful search angle of the user arrival angle corresponding to the maximum, and the user can perform the steps to carry out multi-stage search again and again, until the angle finally needed by a system is obtained, wherein L1 and L2 are positive integers. The careful search method can obtain needed angular resolution (precision), and can also reduce the amount of calculation of the system.

Description

A kind of method of estimating the fixed beam space arrival direction
Technical field
The present invention relates to a kind of mobile communication technology, relate to a kind ofly in mobile radio system or rather, the method for terminal fixed beam space arrival direction is estimated in the base station.
Background technology
Fixed beam space arrival direction (DOA:Direction of Arrival) estimation technique is the product that Estimation of Spatial Spectrum combines with array signal process technique, is widely used in radar, sonar, marine charting, seismology, radio astronomy and the digital mobile communication technical field.
DOA estimates that problem to be solved is how to utilize the output of antenna (or transducer, transducer etc.) array, detects useful signal, and estimates the direction parameter of this useful signal from background noise.
DOA is estimated vital two performance index of weighing the method for estimation quality are spatial resolution and computational complexity.Spatial resolution reflection be the minimum angles interval of each source signal that can distinguish of method of estimation, the resolution of method of estimation is high more, each source signal that can differentiate is just approaching more; The computational complexity of method then embodies the realizability of method of estimation or claims practicality.In the development history of DOA estimation technique, these two technical indicators are determining the direction of advance in this field just.
It also is one of major issue of studying in the intelligent antenna technology field that DOA estimates.How obtaining higher angular resolution and reducing operand as far as possible is the problem that researchers generally are concerned about.In numerous DOA algorithm for estimating, ESPRIT and various variant algorithm thereof can obtain higher relatively precision.Yet this class algorithm mainly is applicable to the large-scale antenna array, and computational complexity is very high.For present mobile communication system, smart antenna is mainly used in the base station.Because be subjected to the influence of digital signal processor (DSP) disposal ability, at present, the array element number in the smart antenna array can not be done very greatly, about 10, will influence the estimated accuracy of ESPRIT algorithm so usually.Therefore, in present mobile communication system, seldom adopt the ESPRIT algorithm to realize that DOA estimates.
Classical Beam-former algorithm, also claiming fixed beam figuration algorithm, is a kind of DOA algorithm for estimating of low complex degree, and it is by searching for (or claiming scanning) to discrete sector region, obtain each regional spatial spectrum parameter, arrive angle according to these parameters then and estimate.Fixed beam figuration algorithm is applicable to various environment, and the DOA resolution that obtains is not subjected to the influence of bay number, therefore is widely used in the present mobile communication system.
Yet the DOA resolution that fixed beam figuration algorithm obtains is relevant with the discrete fanned partition number in space.Usually, number of partitions big more (the search angle is more little), the resolution of DOA increases, however the amount of calculation of system also increases.
Fig. 1 illustrates the method for estimating DOA with classical beamforming algorithm, the antenna array of forming by N bay, N radio frequency (RF) circuit correspondence receives the signal of bay output separately, deliver in the digital baseband signal processor after analog-to-digital conversion (A/D) and handle, its spatial spectrum parameter can be expressed as with formula (1):
P(θ)＝a H(θ)Ra(θ)，θ∈[0，θ cell] (1)
Wherein, P (θ) represents power output, and a (θ) represents direction vector, R representation space correlation matrix.H is a conjugate transpose symbol, and promptly power output is made up of the direction vector and the spatial correlation matrix that arrive angle, and its space search zone is 0 to θ CellCoverage.According to the spatial spectrum parameter of trying to achieve, just can carry out DOA and estimate.Direction vector a (θ) in the formula is different with linear array difference (as even linear array, uniform circular array).
From formula (1) as can be seen, when finding the solution spatial spectrum parameter P (θ), θ is a variable quantity with Δ θ, 0 to θ CellCoverage in, θ of every change, find the solution once corresponding power output P (θ), wherein the pairing arrival angle of maximum P (θ) is the angle of final needs, the more little then DOA of Δ θ estimated accuracy is high more, corresponding spatial spectrum searching times (available L represents) is just big more, yet the system-computed amount is also big more.
Summary of the invention
The objective of the invention is to design a kind of method of estimating the fixed beam space arrival direction,, should obtain certain angular resolution and promptly reach certain precision, reduce the system-computed amount simultaneously again based on fixed beam figuration algorithm.
Method of the present invention is that a kind ofization single-stage is searched for the method for estimation (" the progressively thin search method of classification ") that DOA is asked in a kind of progressively classification refinement that forms for multi-pass decoding.
The technical scheme that realizes the object of the invention is such: a kind of method of estimating the fixed beam space arrival direction is characterized in that comprising:
A. adopt classical beamforming algorithm, the user is carried out the 1st grade of search in the whole beam space scope of aerial array, search for L1 time altogether, from the power output that each time search obtains, take out maximum, obtain the first angle θ that searches that the pairing user of this maximum arrives angle 1, L1 is a positive integer;
B. adopt classical beamforming algorithm, search angle θ just 1Designated space scope [θ 1-θ Cell/ L 1, θ 1+ θ Cell/ L 1] in, this user is carried out the 2nd grade of search, search for L2 time altogether, from the power output that each time search obtains, take out maximum, what obtain that the pairing user of this maximum arrives angle carefully searches angle, wherein, θ CellBe whole beam space scope, L2 is a positive integer.
Described step B repeats, and the next stage search is carried out in the designated space scope of carefully searching angle that the upper level search obtains.
Whole beam space scope in the described steps A is 180 ° when aerial array is even linear array; When being even circular array, aerial array is 360 °.
The power output that described search obtains, the direction vector and the spatial correlation matrix that are arrived angle by signal are formed.
Described spatial correlation matrix can be expressed as user's signal space correlation matrix, user's signal to noise ratio spatial correlation matrix and user's signal interference ratio spatial correlation matrix respectively by power maximal criterion, signal to noise ratio maximal criterion and signal interference ratio maximal criterion.
Described signal space correlation matrix can directly be made up of received signal, or directly is made up of channel response, or directly is made up of certain or some taps of channel response.
The 1st grade of search in the described steps A comprises that also determining time big power output arrives angle with the big pairing user of power output of this time; Judge and the pairing user of peak power output arrive angle and time big pairing user of power output whether arrive angle be between adjacent region; Be between adjacent region the time, will carry out the 2nd grade of search of step B between this adjacent region as the designated space scope of just searching angle; Not between adjacent region the time, intervally carry out the 2nd grade of search of step B respectively as the designated space scope of just searching angle with two.
Described search is carried out separately each user, comprising: one of the search subscriber main diameter direction arrives angle, and perhaps search comprises a plurality of arrival angles of main diameter direction and other footpath directions.
" progressively classification carefully search DOA estimate " of the present invention method, branch two-stage or multistage whole spatial spectrum is searched for.Wherein, the 1st grade is slightly search just, determines the approximate range of user DOA from whole beam space scope; The 2nd grade or the follow-up multistage precise search more of being, according to needed precision, from appointment first search in the spatial dimension (1 grade of search or near the angle that once search obtains) determine the accurate position of user DOA.For every grade of search, we adopt the classical beamforming algorithm shown in the formula (1) to find the solution.
The present invention is directed to the estimation of the DOA under the intelligent antenna technology in the mobile communication system, " the progressively classification refinement DOA method of estimation " that is proposed, have simple, amount of calculation is little and effective advantage, and provides prerequisite for the realization of up-downgoing wave beam forming technology and mobile subscriber's location simultaneously.
The method of classical wave beam (or claiming fixed beam) DOA arrival direction is asked in the employing that the present invention proposes progressively classification refinement, is earlier whole sub-district (perhaps sector) coverage tentatively to be searched for, and obtains the DOA of rough estimate; Precision is as required further searched near the DOA that above-mentioned rough estimate obtains then, obtains more accurate DOA.This " the progressively classification refinement DOA method of estimation " that the present invention proposes can obtain higher DOA estimated accuracy by enough less operands.
Description of drawings
The classical beamforming algorithm of Fig. 1 is estimated the schematic diagram of DOA;
The DOA of Fig. 2 even linear array estimates schematic diagram;
The DOA of Fig. 3 uniform circular array estimates schematic diagram;
The flow chart of Fig. 4 " progressively classification refinement DOA estimates " method;
Fig. 5 " the progressively classification refinement DOA estimates " applicating flow chart of method in the TD-SCDMA system.
Embodiment
The present invention still adopts the fixed beam figuration algorithm of formula (1), obtains the spatial spectrum parameter, estimates DOA:
P(θ)＝a H(θ)Ra(θ)，θ∈[0，θ cell]
Spatial spectrum parameter P (θ) comprises direction vector a (θ) and the spatial correlation matrix R that arrives angle.The direction vector a (θ) that wherein arrives angle is relevant with the type of antenna array.
Referring to Fig. 2, the finding the solution of direction vector a (θ) in the even linear array supposes that with the 1st antenna be reference point shown in the figure, and the direction vector of information source can be expressed as formula (2):
a ( θ ) = f 1 ( θ ) f 1 ( θ ) e j 2 πd λ sin θ · · · f N ( θ ) e j 2 πd λ ( N - 1 ) sin θ T - - - ( 2 )
θ∈[0，180°]
Wherein, N represents bay number (N equals 5 among the figure), and d represents the distance between the antenna, and T is a transpose operator.When each bay is omnidirectional antenna, f 1(θ)=f 2(θ)=...=f N(θ)=1.
Referring to Fig. 3, the finding the solution of direction vector a (θ) in the uniform circular array in uniform circular array, supposes that with center of circle O be reference point shown in the figure, is reference line (0x) with the line of the center of circle and the 1st antenna, and the direction vector of information source can be expressed as formula (3):
a ( θ ) = f 1 ( θ ) e j 2 πr λ cos ( θ ) f 2 ( θ ) e j 2 πr λ cos ( θ - 2 π N ) · · · f N ( θ ) e j 2 πr λ cos ( θ - 2 π ( N - 1 ) N ) T - - - ( 3 )
θ∈[0，360°)
Wherein, r represents radius of circle, and N is the number (N is 8 among the figure) of bay, and T is a transpose operator, and the among the figure is the elevation angle (also being one of parameter of DOA estimation).In like manner, when each bay is omnidirectional antenna, f is arranged 1(θ)=f 2(θ)=...=f N(θ)=1.
Spatial correlation matrix R, different criterions according to the spatial spectrum parameter, comprise power maximal criterion, signal to noise ratio maximal criterion and signal interference ratio maximal criterion, both can represent each user's signal space correlation matrix, also can represent each user's signal to noise ratio spatial correlation matrix, can also represent each user's signal interference ratio spatial correlation matrix.Suppose each user's signal space correlation matrix R XxExpression, noise correlation matrix R nExpression, so, for a certain terminal use k, its spatial correlation matrix R can be expressed as formula (4):
Referring to Fig. 4,, be the flow chart of a secondary refinement DOA method of estimation for the flow chart of the present invention's " progressively classification refinement DOA estimates " method.
Step 41, determine spatial correlation matrix R, according to the different criterions of selecting the spatial spectrum parameter, comprise the power maximal criterion, signal to noise ratio maximal criterion and signal interference ratio maximal criterion, spatial correlation matrix R both can represent each user's signal space correlation matrix, also can represent each user's signal to noise ratio spatial correlation matrix, can also represent each user's signal interference ratio spatial correlation matrix, spatial correlation matrix R and each user's signal space correlation matrix, each user's signal to noise ratio spatial correlation matrix, should satisfy relation between each user's the signal interference ratio spatial correlation matrix by formula (4) decision.
Suppose θ CellRepresent whole space search scope, this scope is relevant with the type of aerial array, for linear array, and θ Cell=180 °; For circle battle array, θ Cell=360 °, use P 1(θ) represent the 1st grade of power output that search obtains, search for L1 time altogether, represent arbitrary time search in the 1st grade of search with i; Use P 2(θ) represent the 2nd grade of power output that search obtains, search for L2 time altogether, represent arbitrary time search in the 2nd grade of search with j, so, the concrete calculation procedure of " progressively classification carefully search DOA estimate " method is as follows:
Step 42 adopts classical beamforming algorithm, to user's (information source) at whole spatial dimension (θ Cell) carrying out the 1st grade of search, searching times is L 1(relevant) with the size of Δ θ, the power output that each search obtains is:
P 1i(θ 1i)＝a H(θ 1i)Ra(θ 1i)，θ 1i∈[0，θ cell]，i＝1…L 1 (5)
Step 43 finds a maximum from the 1st grade of power output that L1 search obtains altogether, the pairing arrival angle of this maximum is just searches angle θ 1, represent with formula (6):
θ 1 ← P ( θ 1 ) = max θ 1 i ( P 1 i ( θ 1 i ) ) , i = 1 · · · L 1 - - - ( 6 )
Step 44 adopts classical beamforming algorithm, at the space specified scope (θ of step 43 acquisition 1Near) in carry out the 2nd grade of search, searching times is L 2(relevant) with the size of Δ θ, the power output that each search obtains is:
P 2j(θ 2j)＝a H(θ 2j)Ra(θ 2j)，θ 2j∈[θ 1-Δθ 1，θ 1+Δθ 1]，j＝1…L 2 (7)
Δθ 1＝θ cell/L 1 (8)
Step 45 finds a maximum from the 2nd grade of power output that L2 search obtains altogether, the pairing arrival angle of this maximum is carefully searches angle θ 2, represent with formula (9):
θ 2 ← P ( θ 2 ) = max θ 2 j ( P 2 j ( θ 2 j ) ) , j = 1 · · · L 2 - - - ( 9 )
If DOA is asked in the two-stage refinement, the θ of this moment 2It is the angle that system finally needs.Ask DOA if adopt the above refinement of secondary, then also need set by step 44,45 identical processes progressively to obtain, when power maximum and pairing arrival angle thereof are asked in back one-level search, be that power maximum and the pairing arrival angle thereof of utilizing the previous stage search to obtain is the hunting zone, promptly with this arrive angle (near) be the space specified hunting zone, carry out next stage search.The designated space scope of each grade search obtains by the corresponding entry that changes in the formula (7) and (8), as: the 3rd level search is to use θ 2Replace θ 1, use Δ θ 2Replace Δ θ 1, and Δ θ 2Be expressed as Δ θ 2=2 Δ θ 1/ L 2The 4th grade of search is to use θ 3Replace θ 1, use Δ θ 3Replace Δ θ 1, and Δ θ 3Be expressed as Δ θ 3=2 Δ θ 2/ L 3And the like, promptly be the process that a region of search reduces refinement step by step.
" progressively the DOA estimation is carefully searched in classification " method is in specific implementation process, can also adopt following transform method: arrival angle maximum and time big power output correspondence is determined in the 1st grade of search, judge whether both are between adjacent region, if (whether angle difference is smaller or equal to Δ θ between adjacent region for both 1), then directly in this interval, search for; If (angle difference is greater than Δ θ not between adjacent region for both 1), then in two intervals, search for respectively, from the power output that two interval search obtain, pick out peak power output, then with the angle of peak power output value correspondence, as the arrival angle of being asked.
Referring to Fig. 5, be example with the TD-SCDMA system below, the application of " progressively the DOA estimation is carefully searched in classification " method is described.
TD-SCDMA (Time Division-Synchronous Code Division Multiple Access) system is based on the mobile communication system of smart antenna, and wherein, the base station adopts array antenna to transmit and receive.For up link, smart antenna plays the micro diversity effect that maximum power ratio merges; For down link, smart antenna carries out wave beam forming, suppresses cochannel and disturbs.In wave beam forming, adopted fixed beam to form algorithm based on DOA, DOA then obtains by " progressively the DOA estimation is carefully searched in classification " of the present invention method.
Step 51,52, by the antenna array that N bay formed, N radio frequency (RF) circuit correspondence receives the signal of bay output separately, delivers in the digital baseband signal processor after analog-to-digital conversion (A/D) and handles.
Step 53,54 is carried out channel estimating and is got spatial correlation matrix.In the TD-SCDMA system, channel estimating is as a key technology, and it reflects the wireless channel of signal experience.Under the situation that sends signal normalization, the spatial correlation matrix of received signal just is equal to the channel space correlation matrix.Therefore, under the indissociable situation of each subscriber signal, directly use signal space correlation matrix R XxCarry out user's differentiation, signal space correlation matrix R XxCan directly form by received signal, also can directly form, also can directly form by certain or some taps of channel response by channel response.Like this, the spatial correlation matrix of user k (total K user) is expressed as:
R xx ( k ) = E { H ( k ) · H ( k ) H } , k = 1 · · · K - - - ( 10 )
E is for averaging H (k)The channel response of expression user k, H is the conjugate transpose symbol, is dot product, it is embodied as:
H (k)＝(h (k、l)…h (k、N)) T (11)
Wherein, h (k, n)The channel response of expression user k respective antenna n, T is a transpose operator in the formula.
h ( k , n ) = h 1 ( k , n ) h 2 ( k , n ) · · · h w ( k , n ) , n = 1 · · · N - - - ( 12 )
Wherein, h (k, n)The channel response tap of expression user k respective antenna n, W represents that the maximal window of channel response is long, it is the sum of tap, what following formula was represented is the stack of the spatial correlation matrix of all channel response taps, also can be the stack of the spatial correlation matrix of certain several channel response tap (a few items in the formula 12).Also can utilize the spatial correlation matrix of each channel response tap (a certain item in the formula 12) to try to achieve each 1 arrival angle during enforcement respectively.
Step 55 utilizes formula (1) to ask the spatial spectrum parameter.For single sub-district, the TD-SCDMA system has adopted the even circular array of omnidirectional's array element, its array element sum N=8.Therefore, the spectrum search volume of its correspondence is 360 °.In the wave beam forming process, because only interested in each user's main footpath usually, like this, when asking DOA, for each user, a demand gets 1 and arrives angle (the arrival angle of main diameter direction).Adopt this " classification " searching method can also try to achieve a plurality of arrival angles of each user (main diameter direction and other footpath directions), can directly utilize whole channel response spatial correlation matrix to try to achieve a plurality of arrival angles when specifically implementing
Step 56, in order to reach 1 ° resolution, system has carried out two-stage search (by step shown in Figure 4), and every grade of search is all carried out even subregion to given sector.Decompose by numeral 360 is carried out approximate number, try to achieve the optimum search number of times of two-stage search: L respectively 1=36, L 2=18.
Step 57,58, the k user's that utilization is finally tried to achieve arrival angle θ (k), the location that can carry out this user respectively reaches by the simplest formula w (k)=a *(θ (k)) carrying out downlink or uplink beamforming algorithm design (can also adopt other down beam shaping formula), * is a conjugate of symbol in the formula.The DOA that promptly adopts the inventive method to obtain can be used for carrying out the user location, also can be used for carrying out the design of downlink or uplink beamforming algorithm.
The present invention can be applicable to and adopts in the communication system that array antenna receives, this " progressively the DOA of classification refinement estimates " method, adopt classical beamforming algorithm, be classified to whole beam space is searched for, progressively obtain the DOA of refinement, can when reaching certain precision, reduce the system-computed amount.

Claims (8)

1. method of estimating the fixed beam space arrival direction is characterized in that comprising:
A. adopt classical beamforming algorithm, the user is carried out the 1st grade of search in the whole beam space scope of aerial array, search for L1 time altogether, from the power output that each time search obtains, take out maximum, obtain the first angle θ that searches that the pairing user of this maximum arrives angle 1, L1 is a positive integer;
B. adopt classical beamforming algorithm, search angle θ just 1Designated space scope [θ 1-θ Cell/ L 1, θ 1+ θ Cell/ L 1] in, this user is carried out the 2nd grade of search, search for L2 time altogether, from the power output that each time search obtains, take out maximum, what obtain that the pairing user of this maximum arrives angle carefully searches angle, wherein, θ CellBe the whole beam space scope of aerial array, L2 is a positive integer.
2. a kind of method of estimating the fixed beam space arrival direction according to claim 1, it is characterized in that: described step B repeats, and the next stage search is carried out in the designated space scope of carefully searching angle that the upper level search obtains.
3. a kind of method of estimating the fixed beam space arrival direction according to claim 1 and 2 is characterized in that: the whole beam space scope in the described steps A is 180 ° when aerial array is even linear array; When being even circular array, aerial array is 360 °.
4. a kind of method of estimating the fixed beam space arrival direction according to claim 1 and 2 is characterized in that: the power output that described search obtains, the direction vector and the spatial correlation matrix that are arrived angle by signal are formed.
5. a kind of method of estimating the fixed beam space arrival direction according to claim 4 is characterized in that: described spatial correlation matrix can be expressed as user's signal space correlation matrix, user's signal to noise ratio spatial correlation matrix and user's signal interference ratio spatial correlation matrix respectively by power maximal criterion, signal to noise ratio maximal criterion and signal interference ratio maximal criterion.
6. a kind of method of estimating the fixed beam space arrival direction according to claim 5, it is characterized in that: described user's signal space correlation matrix can directly be made up of received signal, or the channel response of antenna is directly formed, or certain or some taps of the channel response of antenna are directly formed by the user by the user.
7. a kind of method of estimating the fixed beam space arrival direction according to claim 1 and 2 is characterized in that: the 1st grade of search in the described steps A comprises that also determining time big power output arrives angle with the big pairing user of power output of this time; Judge and the pairing user of peak power output arrive angle and time big pairing user of power output whether arrive angle be between adjacent region; Be between adjacent region the time, will carry out the 2nd grade of search of step B between this adjacent region as the designated space scope of just searching angle; Not between adjacent region the time, intervally carry out the 2nd grade of search of step B respectively as the designated space scope of just searching angle with two.
8. a kind of method of estimating the fixed beam space arrival direction according to claim 1 and 2, it is characterized in that: described search is carried out separately each user, comprise: one of the search subscriber main diameter direction arrives angle, or search comprises a plurality of arrival angles of main diameter direction and other footpath directions.
CN 03121293 2003-03-31 2003-03-31 Method of estimating space reaching direction of fixed wave beam Expired - Lifetime CN1262130C (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN 03121293 CN1262130C (en) 2003-03-31 2003-03-31 Method of estimating space reaching direction of fixed wave beam
PCT/CN2004/000296 WO2004088347A1 (en) 2003-03-31 2004-03-31 Method for estimating the direction of arrival of fixation beam in space

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 03121293 CN1262130C (en) 2003-03-31 2003-03-31 Method of estimating space reaching direction of fixed wave beam

Publications (2)

Publication Number Publication Date
CN1535048A CN1535048A (en) 2004-10-06
CN1262130C true CN1262130C (en) 2006-06-28

Family

ID=33102896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 03121293 Expired - Lifetime CN1262130C (en) 2003-03-31 2003-03-31 Method of estimating space reaching direction of fixed wave beam

Country Status (2)

Country Link
CN (1) CN1262130C (en)
WO (1) WO2004088347A1 (en)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102006114B (en) * 2004-12-17 2012-10-03 华为技术有限公司 Intelligent antenna system and spread-spectrum signal receiver
CN101084633B (en) * 2004-12-17 2011-05-18 华为技术有限公司 Multipath searcher, refinement treatment unit and method for direction of arrival estimation using multipath searcher
CN100388846C (en) * 2005-07-08 2008-05-14 中兴通讯股份有限公司 Downlink beamforming method in mobile communication intelligent antenna system
CN101330304B (en) * 2007-06-22 2012-07-11 中兴通讯股份有限公司 Method and device for calculating arrival direction in an intelligent antenna system
CN101252382B (en) * 2008-03-17 2012-01-25 成都国恒空间技术工程有限公司 Wide frequency range signal polarizing and DOA estimating method and apparatus
CN101615943B (en) * 2008-06-27 2014-04-30 中兴通讯股份有限公司 Method for estimating elevation angle of intelligent antenna multi-subarray system
CN101349742B (en) * 2008-08-29 2010-12-29 西安电子科技大学 Method for optimizing space between broad band phased array elements and measuring frequency and direction of frequency domain multiple targets
CN101487888B (en) * 2009-02-23 2011-01-26 重庆大学 Spacing spectrum peak search method
CN101860387B (en) * 2009-04-09 2013-06-12 中兴通讯股份有限公司 DOA estimation method and system thereof
CN102164374B (en) * 2011-05-10 2013-09-11 华为技术有限公司 Method, device and system for processing beam searching
CN102385048A (en) * 2011-08-10 2012-03-21 西安交通大学 Mixed signal direction estimation method based on even linear array
CN103701515B (en) * 2013-12-11 2017-05-10 北京遥测技术研究所 Digital multi-beam forming method
CN104168047B (en) * 2014-08-13 2017-04-12 电子科技大学 Single-ended time domain beam searching method based on compressed sensing
CN106559118B (en) * 2015-09-24 2020-03-20 中国电信股份有限公司 Method and device for estimating azimuth angle of user terminal under large-scale antenna
CN105891771B (en) * 2016-04-06 2018-01-02 北京邮电大学 It is a kind of improve estimated accuracy based on continuously distributed angle estimating method and equipment
CN106209705B (en) * 2016-09-07 2019-06-14 江苏中兴微通信息科技有限公司 A kind of main diameter estimation method and device of millimeter wave condition of sparse channel
CN106291453B (en) * 2016-09-18 2018-12-18 电子科技大学 A kind of two-dimentional direction of arrival measuring method based on beam space transformation
CN108093360B (en) 2016-11-20 2019-12-24 上海朗帛通信技术有限公司 Method and equipment in UE, base station and service center
US10571557B2 (en) * 2017-06-12 2020-02-25 GM Global Technology Operations LLC Two-stage beamforming
CN109116297B (en) * 2017-06-23 2022-09-27 北京遥感设备研究所 Passive radar space spectrum estimation and beam synthesis combined direction finding method
CN108169710A (en) * 2017-11-16 2018-06-15 捷开通讯(深圳)有限公司 Localization method and alignment system based on reconfigurable antenna
CN108535698B (en) * 2018-04-04 2021-07-20 西安电子科技大学 Meter-wave radar low elevation angle estimation method based on beam space
US11710887B2 (en) * 2018-05-31 2023-07-25 Kymeta Corporation Satellite signal acquisition
CN109188345B (en) * 2018-08-27 2023-03-10 电子科技大学 Coherent signal source DOA estimation method based on pre-delay space-time structure removal
CN109633520B (en) * 2019-01-21 2023-03-24 重庆邮电大学 Uniform circular array super-resolution spatial spectrum estimation method
CN112578354B (en) * 2020-02-28 2024-02-23 加特兰微电子科技(上海)有限公司 Method for determining azimuth angle of target object, computer device and storage medium
CN111679266B (en) * 2020-06-18 2021-09-03 成都纳雷科技有限公司 Automobile millimeter wave radar sparse array grating lobe false target identification method and system
CN112946564B (en) * 2021-04-12 2024-02-02 西北大学 DOA estimation method and device based on DNN (digital optical network) beam space and computer storage medium
CN117278144B (en) * 2023-11-22 2024-02-13 西安迅尔电子有限责任公司 Detection method for low signal-to-noise ratio signal of reconnaissance receiver

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE514830C2 (en) * 1993-11-17 2001-04-30 Ericsson Telefon Ab L M Method and apparatus for directional determination
US5786791A (en) * 1997-02-24 1998-07-28 Motorola, Inc. Method for determining an angle of arrival of a signal transmitted by a remote unit in a communication system
WO1998043106A2 (en) * 1997-03-25 1998-10-01 Siemens Aktiengesellschaft Method for evaluating direction
GB2337171A (en) * 1998-05-06 1999-11-10 Motorola Ltd Direction finder
JP2000147083A (en) * 1998-11-06 2000-05-26 Ntt Mobil Communication Network Inc Method and apparatus for measuring arrival angle
JP4726306B2 (en) * 2001-01-31 2011-07-20 パナソニック株式会社 Wireless communication system, mobile terminal station and direction determination method
CN1199052C (en) * 2001-05-26 2005-04-27 华为技术有限公司 Method of raising accuracy estimation of received signal wave arrival direction

Also Published As

Publication number Publication date
CN1535048A (en) 2004-10-06
WO2004088347A1 (en) 2004-10-14

Similar Documents

Publication Publication Date Title
CN1262130C (en) Method of estimating space reaching direction of fixed wave beam
CN105912791B (en) DOA estimation method based on local search in virtual relatively prime array
CN1242566C (en) Method and appts. for estimating downlink beamforming weights in communications system
CN107505602A (en) DOA estimation method based on DFT under nested battle array
CN106802402B (en) DOA estimation method based on dual-layer Parallel circular array antenna
CN1535046A (en) Wave beam shaping method based on broad band antenna and its device
CN1317891A (en) Method and appts. for reckoning optimum weighting of mobile communication closed transmitting diversity
CN1129237C (en) Digital wave beam forming method and module in radio communication system and its array receiver
CN102045837A (en) Mobile node positioning method and device
CN113747360B (en) Indoor positioning method based on mixed visible light and BLE
CN109116295A (en) The passive direction finding algorithm of baseline is chosen based on phased array
CN110244273A (en) It is a kind of based on the target angle estimation method for being uniformly distributed formula array
CN1574683A (en) Method and device for searching distributed signal sources of base station
CN1207859C (en) Method of time space solution for estimating wave diretion of maultiple paths signals in correlative CDMA and its device
CN113960594A (en) Phased array radar searching method and system
CN1731706A (en) Device for forming annular array beam in VAN
CN104656055A (en) Single signal arrival angle estimating method based on large-scale multi-antenna system
CN1622490A (en) Method and apparatus for implementing omnidirectional coverage of array antennas
CN1199052C (en) Method of raising accuracy estimation of received signal wave arrival direction
CN1475054A (en) Method and apparatus for determining operation condition in communications system
CN113917388A (en) DOA estimation and positioning method using hybrid ADC
CN116405348A (en) Communication perception integrated system array antenna and 5G new air interface angle estimation method
CN1148834C (en) Communication method between base station with N antennae and mobile phone and base for implementing same
CN1186895C (en) Method for generating 3D wave beams in intelligent antenna
CN116232391A (en) Beam training method for ultra-large-scale antenna array, electronic equipment and storage medium

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
PE01 Entry into force of the registration of the contract for pledge of patent right

Effective date of registration: 20070510

Pledge (preservation): Pledge

PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20100413

Granted publication date: 20060628

Pledgee: CHINA DEVELOPMENT BANK

Pledgor: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.|SHANGHAI DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.|CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Registration number: 2007110000354

ASS Succession or assignment of patent right

Owner name: INST OF TELECOMMUNICATION SCIENCE AND TECHNOLGOY

Free format text: FORMER OWNER: DATANG MOBILE COMMUNICATION EQUIPMENT CO., LTD.

Effective date: 20110707

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100083 NO. 40, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING TO: 100191 NO. 40, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING

TR01 Transfer of patent right

Effective date of registration: 20110707

Address after: 100191 Haidian District, Xueyuan Road, No. 40,

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100083 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210607

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

CX01 Expiry of patent term

Granted publication date: 20060628

CX01 Expiry of patent term