CN1761171A - Signal processing system and method based on estimating signal direction of anived wave - Google Patents
Signal processing system and method based on estimating signal direction of anived wave Download PDFInfo
- Publication number
- CN1761171A CN1761171A CN 200410009679 CN200410009679A CN1761171A CN 1761171 A CN1761171 A CN 1761171A CN 200410009679 CN200410009679 CN 200410009679 CN 200410009679 A CN200410009679 A CN 200410009679A CN 1761171 A CN1761171 A CN 1761171A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- antenna
- subarray
- receiver
- 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
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The system includes following parts: array antenna part, antenna array element for inductive receiving signal of remote arrived wave; switching matrix part, connecting or disconnecting node according to instruction of control circuit; signal receiver part, magnifying or attenuating amplitude of signal output from antenna array, carrying out AGC process, carrying out demodulating and decoding so as to convert signal of arrived wave to digital baseband signal; control part, controlling connection and disconnection of switching matrix, controlling operation for accessing data memory; covariance count part for calculating covariance of signal output from array element; subarray composition part (adder), composing equivalent covariance; data memory part, storing intermediate result of data operation and final result; space feature detection part, detecting space feature parameters of signal of arrived wave in order to realize estimating direction of arrived wave.
Description
Technical field
The present invention relates to satellite communication field, relate in particular to the processing unit and the method for the DOA estimation of signal in a kind of satellite communication, this method can be applicable to because system complexity and power consumption constraints processor number are less than the occasion of bay number.
Background technology
In spectrum is estimated, typically compose method of estimation and can be divided into three kinds: the period map method (periodogram) that Schuster in 1898 proposes is based on the direct spectrum method of estimation of Fourier conversion; Nineteen twenty-seven Yule propose based on AR, the model parameter The Power Spectrum Estimation Method of MA and arma modeling; The super-resolution spectrum method of estimation of Schmidt proposition in 1986, the MUSIC method is typical case's representative, is to use one of more method.MUSIC super-resolution spectrum method of estimation is based on the arrived wave signal Subspace Decomposition, and is closely related with array antenna array element spatial distribution characteristic (array manifold).The ESPRIT algorithm that people such as Roy propose also belongs to the super-resolution spectrum algorithm for estimating that decomposes based on proper subspace, this method has been utilized the symmetry of array antenna array element in spatial distribution, require the received signal vector space to have rotational invariance on the invariant subspace of spatial domain, the MUSIC method distributes to array element does not have the symmetry requirement.
DOA method of estimation based on the characteristic vector Subspace Decomposition is: the additive noise in arrived wave signal and array element adds under the mutual incoherent condition of interference signal, the full signal space of array antenna output signal can be decomposed into the subspace of two quadratures, the noise subspace that signal subspace that signal phasor is dominant and noise vector are dominant.The projection of vector in the signal subspace on noise subspace is zero vector.So the inner product of signal phasor and noise vector equals zero, it is infinitely great that its inverse is tending towards.Based on this, search just can obtain arrived wave signal direction (DOA) angle with the pairing wave angle that reaches of the null vector of noise vector inner product.
In actual applications, the data that we can observe are data of signal and noise stack, can be expressed as the output random data sequence x of array antenna
i(n), i=1,2 ..., M; N=1,2 ..., N.Wherein i represents the locus numbering of array antenna array element with respect to reference array element, and M is an array number; N represents the snap sampling constantly, and N is fast umber of beats.Be without loss of generality, can suppose that array dateout random sequence is that the interior broad sense of spatial domain and time-domain is steady, the random process of traversal.At M>K, during N → ∞, can obtain array antenna output signal sequence x
i(n) ensemble average value, i.e. the valuation of covariance matrix,
R wherein
IjBe array antenna i, the assessment of covariation of j delayed output signals, K are that information source reaches wave number.Carry out DOA and estimate, need calculate the valuation of covariance matrix earlier
Each information source reach ripple be on spatial domain and the time-domain broad sense steadily, during ergodic random process, the ensemble average covariance matrix
Be equivalent to the statistical average covariance matrix of delayed output signals in the spatial domain.Angular resolution based on Subspace Decomposition DOA algorithm for estimating depends on the covariance matrix valuation
Levels of precision.Valuation
More little with the error of assembly average R, the angular resolution that DOA estimates is high more.The covariance matrix valuation of array output signal
Calculating be to carry out the important content that DOA estimates.
In the DOA algorithm for estimating, according to different receiving system structures, the computational methods of covariance matrix can be divided into two classes, a class is based on the multi-receiver system configuration, and second class is based on the single-receiver system configuration.In existing system based on multi-receiver, each array element of array antenna is corresponding one by one with receiver (signal processor): logically, the array element output is directly connected to receiver inlet, in each snap sampling constantly, receiver needs to receive, handle the data of whole array element outputs, therefore, the quantity of receiver (signal processor) equates with bay quantity.Because the quantity of system receiver equals bay quantity, in the bigger application scenario of bay number, the corresponding increase of receiver quantity, on the one hand, system cost increases, and on the other hand, the total power consumption of receiver (signal processor) will increase rapidly.This type systematic can not be used in equipment power dissipation, the strict limited spaceborne communication system of weight and volume.Second type systematic has only a receiver, the output of each array element of array antenna links to each other with the phase shifter/delayer of a series of complex value weight coefficients, each phase shifter/delayer output links to each other with the adder input, and the output of adder links to each other with receiver (signal processor) input.The hardware cost of this system is lower, is equivalent to a time domain serial received system.Its shortcoming is when the bay number is bigger, and required weighting phase shift phase device and delayer quantity also increase, and calculates complex value amplitude and the required operand of phase weighting coefficient and increases sharply, and is unfavorable for the data in real time processing.These two kinds of systems all can not satisfy the requirement of satellite communication, are not suitable for using in satellite communication.
Summary of the invention
Based on the problems referred to above, be that common array signal processing system requires array antenna array element quantity to equal the quantity of signal receiver usually, and in some applications, perhaps can't satisfy this condition, and perhaps needing to make the two unequal---receiver quantity is less than array element quantity.The present invention discloses a kind of signal processing system, a kind of ripple DOA method of estimation and a kind of best subarray quantity computation method of reaching.
In satellite communication, spaceborne receiver quantity is lacked than array antenna array number M.Estimate resolution in order to improve DOA, must make full use of the array antenna resource, make its effective aperture reach maximum.This requirement must all use the ripple snap that reaches on whole array elements of antenna.On the other hand, in each snap sampling constantly, L
rIndividual receiver at most can only receiving array antenna L
rOutput snap on the individual array element.
Therefore, an object of the present invention is to provide a kind of array signal processing system, according in the array signal processing system of the present invention, a kind of time division multiplexing receiver structure scheme is disclosed, operation principle is as follows: different with traditional reception programme, in this reception structure, array antenna is not directly to link to each other with signal receiver, but links to each other with receiver by a programme-controlled exchange matrix.The programme-controlled exchange matrix is equivalent to electronic switch, and each node has connection and disconnects two operating states.In any snap sampling constantly, the programme-controlled exchange matrix makes always has L in the array antenna
rIndividual array element and L
rThe input of individual receiver links to each other; In different snap samplings constantly, the programme-controlled exchange matrix may switch connected node, makes another group L of receiver and array antenna
rIndividual array element links to each other.In the time of need switching, the L in the programme-controlled exchange matrix
rIndividual node switches simultaneously.Represent the programme-controlled exchange battle array with sw, then have:
Line display receiver numbering in formula (1) the sw matrix, tabulation are shown the array element numbering of main array antenna, k
I, jExpression node connection status, k
I, j∈ (0,1), k
I, j=0 expression does not connect k
I, j=1 expression receiver i links to each other with bay j, and then array element is exported fast beat of data and is input in the receiver.i=1,2,…L
r,j=0,1,…,L-1。
In foundation array signal processing system of the present invention, each array element of antenna realizes that with receiver grouping, timesharing are connected, and each snap sampling always has L constantly
rIndividual array element links to each other with receiver; Receiver constantly can receive L in each snap sampling
rThe fast beat of data of the output of individual array element, receiver is in running order always.
According in the array signal processing system of the present invention, can receive the fast beat of data of output of any amount bay with a small amount of receiver.Adopt the working method of time division multiplexing receiver, improved the utilance of receiver, reduced the quantity of required receiver, make the total power consumption of receiver controlled.
According in the array signal processing system of the present invention, signal processing method can both be connected all array elements of main array antenna with receiver, and the array element resource is fully used, thereby effective aperture that can the expanded matrix antenna, improves the DOA estimated resolution.
According in the array signal processing system of the present invention, the dateout of receiver structure on can the whole array elements of receiving array antenna obtains the full detail of arrived wave signal, reaches ripple DOA estimated resolution and reaches maximum thereby can make.What but scheme adopted is the mode of time division multiplexing receiver, and receiver receives only L constantly in each snap sampling
rThe output snap of individual array element.Array antenna is actual to be made up of L array element, the bay sum is greater than the receiver sum, therefore receiver can not receive data on the whole array elements of main array antenna constantly a snap sampling, therefore, can not directly calculate the covariance matrix of array antenna output snap in a conventional manner.For this reason, the present invention disclose a kind of array signal processing method-based on subarray synthetic reach the ripple direction determining method
The receiver of signal processing system involved in the present invention (signal processor) quantity is between the receiver quantity of single-receiver and multi-receiver system, promptly greater than 1, less than array antenna array number (M), concrete quantity is according to applied environment condition and application requirements adjustment.
According in the array signal processing system of the present invention, to the array antenna type without limits, both can be uniform line-array, also can be the linear array of irregular alignment, or the array of plane distribution.
According in the array signal processing system of the present invention, the quantity of bay not to be strict with, array antenna is minimum can to comprise 2 array elements, and array element quantity for example reaches the wave number amount, processing speed etc. by the practical application conditional decision.
According in the array signal processing system of the present invention, the relative populations of signal receiver and bay is not had strict the restriction, but receiver quantity generally is less than bay quantity.When equaling bay quantity, receiver quantity just has been degenerated to common array signal processing system.
First aspect of the present invention is a kind of array antenna signal processing system, comprise: a kind of treatment system and corresponding processing method of estimating based on the synthetic DOA of subarray, the arrived wave signal direction (DOA) that this treatment system realization receiver quantity is less than under the bay quantity situation is estimated.This treatment system comprises: the array antenna part, and bay (1-1 to 1-N) induction receives the far field arrived wave signal; The switching matrix part is according to the instruction connection or the disconnected node of control circuit; The signal receiver part is amplified or decay bay amplitude output signal, carries out AGC and handles, and carries out demodulation, decoding is transformed to digital baseband signal with arrived wave signal; Control section, the connection or the disconnection of control switching matrix, the accessing operation of control data memory; The covariance calculating section, the covariance of calculating delayed output signals; Subarray composition part (adder), synthetic equivalent covariance; Data memory part storing intermediate result of data operation and final result; The space characteristics test section, the space characteristic parameter of detection arrived wave signal is realized reaching ripple direction (DOA) and is estimated.
Another aspect of the present invention provides a kind of signal processing method, and this method realizes reaching the super-resolution estimation of ripple direction (DOA) based on the arrived wave signal snap.Described method comprises: based on the time division multiplexing processing method, ask the covariance matrix of sub antenna array output; Ask the equivalent covariance matrix of main antenna array output; Reaching ripple DOA estimates.
The 3rd aspect of the present invention discloses to be optimized criterion, method, the scheme of choosing sub-array antenna and counts computing formula.
According to the present invention, can realize in the communication system, signal processing system of bay number that the super-resolution of arrived wave signal direction (DOA) estimates because system complexity and power consumption constraints, signal receiver (signal processor) number are less than.
System involved in the present invention can be applicable to a plurality of occasions, comprising: satellite communication system, onboard radar system, multiple-input and multiple-output (MIMO) mobile communication system, seismic signal detection system, subterranean resource exploration etc.
In foundation array signal processing system of the present invention, main array antenna 1-1 to 1-L is divided into several subarray antennas, the array number of each subarray antenna is less than the array number of main array antenna, and the array number of each subarray antenna equals the quantity of signal receiver.The array element of subarray antenna does not directly link to each other with receiver, but by switching matrix, adopts time division multiplexing mode to be connected with the input of receiver, realizes the time division multiplexing of receiver with this.The array number of each subarray equates.
In according to array signal processing system of the present invention and respective handling method, according to the array element numbering, main array antenna is divided into M sub-array antenna, each subarray antenna is by L
rIndividual array element is formed, and equals receiver quantity.The principle of dividing is: (1) exports fast beat of data by whole subarrays can synthesize the array antenna output correlation matrix of exporting the correlation matrix equivalence with main array, is called equivalent main array output correlation matrix; This requires any two array elements of main array antenna to occur once simultaneously in some subarray antennas at least.(2) minimum number of subarray antenna.
In according to array signal processing system of the present invention and respective handling method, because receiver quantity L
rBe less than main array antenna array number L, receiver constantly can only be received L in each snap sampling
rThe output snap of individual array element can't receive all the other L-L constantly in same sampling
rOutput snap on the individual array element.So, can not obtain to calculate the output snap that main array antenna is exported the required whole array elements of covariance matrix constantly a snap sampling.Yet, when arrived wave signal and the characteristic of channel satisfy certain condition, can obtain to calculate the required whole snaps of covariance by the mode of time-division processing.
In according to array signal processing system of the present invention and respective handling method, it is as follows to propose processing method at this situation: (1) is in each snap sampling moment, to the L of a sub-array antenna
rDateout on the individual array element is carried out synchronized sampling, obtains L
rThe output of individual array element reach wave datum; (2) constantly sampling is implemented in the output of different subarray bays in different snap samplings, through obtaining the snap of whole subarray bay outputs after several snap samplings constantly.(3) according to the signal processing method among the present invention the whole fast beat of data that receives is handled, can be calculated required covariance matrix valuation.
In foundation array signal processing system of the present invention, receiver is implemented the data of the whole array elements outputs of sub-array antenna constantly in each snap sampling and is taken a sample.By time-sharing multiplex, can receive the fast beat of data of output of handling any amount bay, realize reception and processing to a large amount of delayed output signals.Compare with other scheme, under certain condition, this programme can significantly improve systematic function, improves the precision that reaches ripple DOA resolving power and spectrum estimation, reduces the main lobe width that the digital wave number form becomes.
Description of drawings
Below in conjunction with accompanying drawing the present invention is illustrated in further detail:
Fig. 1 is according to array antenna signal processing system function composition frame chart of the present invention;
Fig. 2-Fig. 5 is the DOA estimation curve under different condition, and the estimation curve parameter among each figure is respectively:
Fig. 2: it is 15,45 ,-45 degree that fast umber of beats N=40, signal to noise ratio snr=20dB, ripple reach the angle.
Fig. 3: it is 15,19 ,-60 degree that fast umber of beats N=400, signal to noise ratio snr=20dB, ripple reach the angle.
Fig. 4: it is 60,64 ,-30 degree that fast umber of beats N=50, signal to noise ratio snr=100dB, ripple reach the angle.
Fig. 5: it is 60,64 ,-30 degree that fast umber of beats N=10000, signal to noise ratio snr=30dB, ripple reach the angle.
Most preferred embodiment is described in detail
Below with reference to accompanying drawing of the present invention, more detailed description goes out most preferred embodiment of the present invention.
Fig. 1 forms structured flowchart according to the array antenna signal processing system function of the specific embodiment of the invention.Array antenna signal processing system shown in Figure 1 can be installed on the satellite, also can be installed in other array signal processing system, for example in the mobile communication equipment.
Explain the specific embodiment of the invention with reference to Fig. 1 below:
Array antenna in the signal processing system involved in the present invention (1) receives the signal that sends from one or more far fields equipment (for example, ground satellite station, base station etc.); The on off operating mode of control section (4) control switching matrix (2) node; Signal receiver (3) receives the arrived wave signal of switching matrix (2) output, becomes digital baseband signal after treatment; Covariance calculating section (5) calculates the estimated value of arrived wave signal correlation matrix; Subarray composition part (6) is based on the synthetic equivalent main array covariance matrix of valuation of subarray output covariance matrix, and the result is stored in the data memory (7); Wave angle is estimated to reach based on the covariance matrix valuation in space characteristics test section (8), realizes that super-resolution reaches ripple direction (DOA) and estimates.Specific implementation method is as follows:
Among Fig. 1, array antenna part (1) comprises bay 1-1 to 1-L, and L is greater than 1.In foundation array signal processing system of the present invention, the array antenna of being made up of whole array element 1-1 to 1-L is defined as main array antenna; The array antenna of being made up of part array element is defined as the subarray antenna.As required, different subarray antennas can comprise identical array element.
In according to array signal processing system of the present invention and corresponding signal process method, the quantity of basis signal receiver, main array antenna is divided into several subarray antennas, and each subarray antenna comprises the array element of equal number, and array element quantity equals signal receiver quantity.
In according to array signal processing system of the present invention and corresponding signal process method, according to sub-array partition criterion disclosed by the invention main array antenna array element to be divided into groups, each group array element is formed a sub-array antenna.Suppose that each subarray antenna is by L
rIndividual array element is formed, and promptly hypothesis has L
rIndividual receiver, 2≤L
r≤ L.The subarray antenna is serial number successively.
In foundation array signal processing system of the present invention, array antenna part (1) receives the wireless signal that comes from far-end, and the arrived wave signal complex amplitude of each array element output has different phase places.The signal of array element 1-1 to 1-N output becomes the input signal 2-i1 to 2-iL of switching matrix (2).Signal expression is as follows:
Be without loss of generality, establish: array antenna (1) is a uniform line-array, adjacent array element distance
λ is an arrived wave signal centre wavelength; Each subarray is by L
rIndividual array element is formed; K arrived wave signal, far field arrived wave signal complex amplitude is S
k(n), the angle of pitch is θ
k, uncorrelated mutually between signal and the signal; N (n) is an additive noise, obeys multiple Gaussian Profile, and average is zero, and variance is σ
2, uncorrelated mutually between the noise on each array element, uncorrelated mutually between noise and the arrived wave signal.Then array element l at the arrived wave signal complex amplitude of sampling n induction constantly is
k=1,2,...,K,L=1,2,…,L。y
l(n) n corresponding to l array element exports snap.All L array element at the arrived wave signal complex amplitude vector of sampling n induction constantly is
Y(n)
T=[y
1(n)y
2(n)…y
L(n)]+N(n) (3)
Y (n)
TCorresponding to the output snap of array antenna at sampling moment n.
Order
Then have:
A(θ)=[a(θ
1),a(θ
2),…,a(θ
K)] (4)
A (θ
k) be reach wave angle be θ k reach the ripple direction vector.
s
k(n) be to reach the arrived wave signal base band complex envelope that wave angle is θ k.
N (n)
TIt is array antenna additive noise vector.
Under the effect of control section (4), switching matrix (2) is selected L from the signal that input 2-i1 to 2-iL receives
rIndividual, 2-1 to 2-Lr exports by output.Node K in the switching matrix (2)
11To K
LrLOperating state determine by the output signal 4-1 of control section (4).Control section (4) is by port 4-1 output control signal, and in return the node of matrix (2) is selected signal, is input to the control input end 2-c1 of switching matrix (2).Node K when switching matrix (2)
MnWhen not selected, K
Mn=0, the input 2-in of switching matrix (2) and output 2-m disconnect, input signal y
n(n) can not be sent to the output of switching matrix (2), this snap is rejected; Choose the node K of switching matrix (2) when control section (4)
MnThe time, node K
Mn=1, the input 2-in of switching matrix (2) is communicated with output 2-m, input signal y
n(n) the output port 2-m by switching matrix (2) exports.
In foundation array signal processing system of the present invention, according to array signal processing method disclosed by the invention, in each snap sampling constantly, control section (4) will be controlled the operating state of each node in the switching matrix (2), the whole output 2-1 to 2-Lr that make switching matrix (2) always can with the L among switching matrix (2) the input 2-i1 to 2-iL
rIndividual input is communicated with, with remaining L-L
rIndividual input disconnects and connecting.The L of switching matrix (2)
rThe L of individual input and array antenna
rIndividual array element is corresponding one by one.Therefore, under the effect of control section (4), switching matrix (2) is L in the sampling of each snap constantly all can be array antenna
rOutput signal on the individual array element is sent to output 2-1 to 2-Lr undampedly.
According to array signal processing method disclosed by the invention, in each snap sampling constantly, control section (4) is selected L from L input arrived wave signal of switching matrix (2)
rIndividual arrived wave signal sends signal receiver (3) to by switching matrix (2).
Switching matrix (2) passes through output 2-1 to 2-Lr with L
rThe arrived wave signal of individual array element is sent to the input 3-il to 3-iL of signal receiver part (3) respectively
r
In foundation array signal processing system of the present invention, signal receiver (3) comprises automatic gain control (AGC), demodulating unit, decoding unit, analog to digital conversion (A/D) circuit.Agc circuit carries out AGC to received signal to be handled, and makes the arrived wave signal amplitude of reception satisfy the job requirement of subsequent conditioning circuit; Demodulating unit carries out demodulation to the signal of handling through AGC, becomes baseband complex signal, and baseband complex signal contains the relative amplitude and the phase information of arrived wave signal; The A/D translation circuit is transformed to the digital baseband complex signal with the Analog Baseband complex signal, so that system applies DSP technology.Signal receiver (3) is by output 3-01 to 3-0Lr output arrived wave signal snap y
1(n) to y
Lr(n).Wherein, y
l(n), l=1 ..., Lr is corresponding to any L in L the arrived wave signal
rIndividual arrived wave signal.
The output signal of signal receiver (3) is the digital baseband complex signal, is input to covariance calculating section (5) by input 5-1 to 5-Lr.Covariance calculating section (5) is finished auto-correlation and the cross-correlation calculation to input signal, and the matrix that its assessment of covariation that is output as each input signal is formed is defined as covariance matrix, R
YYThe form of covariance matrix:
In the formula, r
IjBe bay i, the covariance between the j dateout, y represents array antenna output signal vector, i, j ∈ (1 ..., Lr), Lr is an array number.
By formula (8) as seen, covariance matrix R
YYBe the covariance matrix of array antenna output signal vector on spatial domain, form by the covariance between each delayed output signals of antenna.
In foundation array signal processing system of the present invention, according to array signal processing method disclosed by the invention, constantly the data on the whole Lr array elements in the subarray are taken a sample simultaneously in each snap sampling, obtain an output snap of Lr array element, according to the error precision requirement, each sub antenna array is repeated the sampling of N correlations, obtain N snap of Lr array element, can calculate according to formula (8).Introduce concrete processing method below.
The covariance matrix R of array antenna output signal
YYBe defined as:
R
YY=E{[A(θ)S(n)+N(n)][A(θ)S(n)+N(n)]
H} (9)
Wherein E{} represents to get the statistical expectation of bracket interior element.Formula (9) can further be written as:
R
YY=E{[A(θ)S(n)+N(n)][S(n)
HA(θ)H+N(n)
H]}
=E{A(θ)S(n)S(n)
HA(θ)
H+N(n)S(n)
HA(θ)
H+A(θ)S(n)N(n)
H+N(n)N(n)
H} (10)
Under the separate condition of arrived wave signal and noise, the cross-product term in the formula (10) equals zero, so have
R
YY=E{A(θ)S(n)S(n)
HA(θ)
H}+E{N(n)N(n)
H}
=A(θ)R
S(n)A(θ)
H+σ
2I (11)
Wherein, I is Lr * Lr rank unit matrix.
By formula (11) as can be known, need ask source signal covariance matrix R earlier
S, just can obtain the covariance matrix R of array antenna output signal again according to formula (11)
YY
R
S=E{[s
1(n)…s
K(n)]
T[s
1(n)…s
K(n)]
*} (12)
When disregarding noise, array antenna output covariance matrix is R
XX:
R
XX=A(θ)R
S(l)A(θ)
H (13)
When considering additive noise, the covariance matrix of array antenna output signal is following form:
R
YY=R
XX+σ
2I (14)
σ wherein
2Be noise variance, I is a unit matrix.
Can obtain the valuation of covariance matrix with N snap of array output at the mean value of time domain
Can obtain the covariance matrix valuation of any one subarray antenna output according to formula (15).
Can determine required whole subarrays according to subarray antenna selection criterion disclosed by the invention and method.
According to signal processing method disclosed by the invention,, successively whole subarray antenna output signals are taken a sample according to the sequential of determining.Behind the snap sampling number of times through regulation, can obtain the output snap on whole subarray antennas, the whole array element.
According to signal processing method disclosed by the invention, under the coordination of control section (4), covariance calculating section (5) is obtained the covariance matrix valuation of each subarray antenna output snap one by one according to formula (15).Result of calculation is sent to data memory part (7) by the output 5-0 output of covariance matrix calculating section (5), and as the intermediate object program of calculating, each subarray output covariance matrix valuation is temporary among the data memory (7).
The output of each sub antenna array produces a covariance matrix.Generally speaking, the quantity of subarray antenna equals
Wherein, L is main array antenna array number, and Lr is a subarray bay number, and it equals signal receiver quantity, C
L LrIt is the number of combinations of selecting Lr element in L the element.
Signal processing method disclosed by the invention is: the covariance matrix that calculates each subarray array element output; According to signal processing criterion disclosed by the invention,, construct main array antenna output equivalent covariance matrix by subarray output covariance matrix based on the subarray synthetic technology; Use based on the algorithm of eigenmatrix singular value decomposition arrived wave signal space characteristics parameter is estimated, realize reaching ripple DOA and estimate.The detailed construction method of main array antenna equivalence covariance matrix is described with an embodiment below.
With reference to Fig. 1, to establish main array antenna array number and equal L=4, the subarray bay is counted Lr=2.Sequence number is 1 array element array element for referencial use, and the array element number consecutively is 1,2,3,4.Then the subarray antenna number equals
Six subarrays are respectively by array element (1,2), and (2,3), (3,4) (4,1), (1,3), (2,4) are formed, and subarray stream shape is respectively:
The covariance matrix of subarray antenna output is R
YYi:
R
YYi=A
i(θ)R
SiA
i(θ)
H+σ
2I (22)
Wherein I is 2 * 2 rank unit matrix, i=1, and 2 ..., 6.
In form, R
SiWith R
SBeing identical, all is K * kth moment battle array, is 3 * 3 rank here, is the cross-correlation matrix between the source signal, and it is mutual when uncorrelated to reach ripple in information source, R
Si, R
SBe the diagonal angle square formation.In fact the corresponding constantly arrived wave signal S of different snap samplings
i(n) and additive noise N
i(n) incomplete same, so R
Si≠ R
S, i represents i subarray.Construct equivalent main array in order to use based on the subarray synthetic method, require ripple reach signal be broad sense steadily, in the random process of spatial domain and time domain traversal.At this moment, can think R
Si=R
S,
σ wherein
i 2Be the variance of subarray antenna i additive noise, σ
2It is the variance of main array antenna additive noise.
In arrived wave signal processing method disclosed by the invention, the sampling method and the process of regulation subarray output snap are as follows: control section (4) sends a control command to switching matrix (2) constantly in each sampling, choose a wherein group node, Lr output node of switching matrix (2) is communicated with Lr input node, the output snap on Lr the array element is sent to Lr signal receiver (3).
Control section (4) is chosen Lr the input node that links to each other with first subarray bay in the switching matrix (2) constantly in first snap sampling, constantly choose Lr the input node that links to each other with second sub-array antenna array element in the switching matrix (2) second snap sampling, and the like, circulation is carried out.
Like this, what signal receiver (3) constantly received signal in the sampling of first snap is the output snap of first subarray antenna, and what constantly receive signal second snap sampling is the output snap of second sub-array antenna.。。, what constantly receive signal in the sampling of n snap is the output snap of n sub-array antenna.
Signal receiver (3) continuous operation is taken a sample to the whole Lr of first subarray array element outputs constantly in first snap sampling; Second snap sampling taken a sample to second whole Lr of subarray array element outputs constantly ..., n snap sampling is constantly to n the whole Lr of subarray array element output samplings.So successively Lr array element dateout of each subarray antenna sampled.Finish after the sampling process to Lr array element output of last subarray antenna, restart Lr array element output sampling to first subarray antenna again, so circulation is carried out.Can require to determine hits according to different valuation error precisions.Each sampling sampling snap constantly is temporary in the data memory (7).
In arrived wave signal processing method disclosed by the invention, definition: whole subarray antennas required fast umber of beats of once sampling respectively is called sampling period of subarray antenna.Therefore, equal based on the C of the direct computing method of main array output doubly based on the synthetic sampling period of subarray, wherein C is the subarray antenna amount.
Based on the output of signal receiver (3), covariance calculating section (5) can go out to calculate the output signal covariance matrix of subarray antenna in proper order, is respectively R
YY1, R
YY2, R
YY3, R
YY4, R
YY5, R
YY6:
Unification is write as following form:
In the formula, r
Ij(k) be covariance between array element i and the array element j output signal, k represents r
Ij(k) this number of times that in a sampling period of subarray, occurs, for example r
11(2) be illustrated in the sampling period and occur covariance item r for the second time
11In the present embodiment, a sampling period of subarray equals 6 snap sampling periods, supposes that sampling frequency equals 1.
Covariance calculating section (5) passes through data/address bus 5-0 with result of calculation R
YY1, R
YY2, R
YY3, R
YY4, R
YY5, R
YY6Send to the input/output port 7-1 of data memory (7).
Control section (4) is read subarray covariance matrix R from data memory (7)
YY1, R
YY2, R
YY3, R
YY4, R
YY5, R
YY6, data are sent to the port 6-1 of subarray composition part (6) by bus 7-2.
In foundation signal processing system of the present invention, subarray composition part (6) is based on the covariance matrix R of subarray
YY1, R
YY2, R
XX3, R
YY4, R
YY5And R
YY6, according to signal processing method disclosed by the invention, the synthetic covariance matrix valuation R that is equivalent to main array output
YY
In the formula, r
Ij(n) n in represents r
IjThe number of times that this repeats in the subarray covariance is calculated.
According to the control command that control section (4) sends to subarray composition part (6) by port 4-3, subarray composition part (6) is with equivalent covariance matrix R
YYBy port 6-1 output, the data/address bus 7-2 input by memory part (7) is temporary in the memory (7).
In according to signal processing system of the present invention and signal processing method, the covariance item that repeats is taked following processing method:
1) keep the value that occurs for the first time, the value of abandoning repeating:
r
ij=r
ij(1) (31)
r
Ij(1) is the value that occurs first.
2) or, get the value of the arithmetic mean of each time valuation as this:
N is r
Ij(n) number of times of Chu Xianing.
3) or, the mean value of N valuation of getting equivalent covariance matrix is as covariance matrix:
N is the sampling period number.
When calculating main array equivalence covariance matrix based on the subarray synthetic method, need C fast beat of data, wherein, first fast beat of data is the sampling value of first subarray antenna output, second fast beat of data is the sampling value of second sub-array antenna output,, C fast beat of data is the sampling value of C sub-array antenna output.
Compare with other signal processing system, by increasing the fast umber of beats of sampling, proposed by the invention has following advantage based on subarray synthetic signal processing system and signal processing method: make the quantity of desired signal receiver (3) be reduced to Lr from L, Lr<L, simultaneously, do not reduce the main performance of signal processing system, as resolution, beamwidth etc.
Subarray composition part (6) outputs to data memory (7) afterwards with equivalent covariance matrix data, and control section (4) sends instruction by port 4-4 to data memory (7), and 7-c1 is input to data memory by port.This instruction makes data memory (7) from port 7-3 equivalent covariance matrix be exported, and is sent to space characteristics test section (8) by port 8-2.
Space characteristics test section (8) is by the instruction of port 8-1 reception from control section (4), read in based on the synthetic equivalent covariance matrix data of subarray from port 8-2, according to DOA method of estimation based on the eigenmatrix singular value decomposition, equivalent covariance matrix is carried out singular value decomposition, extract the matrix exgenvalue valuation, carry out parameter Estimation, realize the estimation of arrived wave signal DOA super-resolution.Ripple reaches the port 8-3 output of angular data by space characteristics test section (8).
In according to signal processing system of the present invention and signal processing method, reach ripple DOA estimation approach and be not limited to the MUSIC algorithm that the embodiment of the invention adopts.
Signal processing system disclosed in this invention and signal processing method are applicable to any method for parameter estimation based on Singular Value Decomposition Using, estimate wave beam formation etc. comprising reaching ripple DOA.
Optimize the scheme number of choosing the subarray antenna:
In the embodiment of above-mentioned foundation signal processing system of the present invention and signal processing method, the subarray antenna amount is C,
Discover that during the array number Lr that comprises when subarray>2, when calculating equivalent covariance matrix, required subarray antenna amount can be less than C
L LrIn the content below, the present invention will disclose the computational methods formula of required subarray antenna minimum number; Also will openly choose the criterion of required subarray antenna, and be called to optimize and choose criterion.
Optimize and choose criterion:
1. the number of times that occurs in different subarray antennas of any one array element is few as far as possible.
2. any two array elements occur once in a sub-array antenna at least simultaneously.
3. count under the selectable situation of Lr at the subarray bay, get the selectable maximum of Lr.
4. under the condition that main array array number L and subarray array number Lr determine, get the minimum scheme of choosing of subarray antenna number.
Choose criterion according to this optimization, can optimize the scheme of choosing and count K and be:
For reducing the disturbance of random noise, can limit minimum number P that array element occurs greater than certain threshold value in subarray, for example P 〉=50, the i.e. minimum number of coherent sampling to covariance.
Embodiment performance simulation result:
Accompanying drawing 2-5 compares with the ripple DOA estimation performance simulation result that reaches that adopts main array output based on the synthetic ripple DOA estimated result that reaches of subarray for employing.
Adopt the MUSIC algorithm, to carrying out emulation, pseudo-spectral function P based on subarray synthetic arrived wave signal treatment system and signal processing method according to of the present invention
M U(θ) be
Simulated conditions: 3 information sources reach ripple, and main array antenna is made up of four array elements, and the subarray antenna comprises two array elements.
Each picture group comprises two sets of curves respectively.Wherein, last figure is based on the synthetic DOA estimated result of subarray, and figure below is the DOA estimated result based on main array output.
Fig. 2: it is 15,45 ,-45 degree that fast umber of beats N=40, signal to noise ratio snr=20dB, ripple reach the angle.
Fig. 3: it is 15,19 ,-60 degree that fast umber of beats N=400, signal to noise ratio snr=20dB, ripple reach the angle.
Fig. 4: it is 60,64 ,-30 degree that fast umber of beats N=50, signal to noise ratio snr=100dB, ripple reach the angle.
Fig. 5: it is 60,64 ,-30 degree that fast umber of beats N=10000, signal to noise ratio snr=30dB, ripple reach the angle.
The embodiment simulation result shows:
Among Fig. 2, the DOA angle of arrived wave signal differs bigger.The result shows, under medium signal to noise ratio condition, use less sampling snap just can estimate exactly that efferent echo reaches the angle, from two figure up and down as can be known, the DOA estimated result of exporting fast beat of data based on synthetic DOA estimated result of subarray and the direct main array of employing is identical.
Among Fig. 3, the DOA angle of arrived wave signal differs smaller.The result shows, can tell different arrived wave signals based on the DOA algorithm for estimating that subarray is synthetic.From two figure up and down as can be known, under identical signal to noise ratio condition, needing increases fast umber of beats, just can make the DOA estimated result more accurate.
Among Fig. 4, it is identical with Fig. 3 that signal reaches wave angle, also gets smaller value.Compare with Fig. 3, simulated conditions changes into: reduce the fast umber of beats of sampling, increase the arrived wave signal signal to noise ratio.The result shows, under the high s/n ratio condition, what very little fast umber of beats just can be told close proximity reaches the ripple direction.
Among Fig. 5, it is bigger that signal reaches wave angle, and adjacent to reach wave angle degree difference less.Signal to noise ratio is got representative value, increases fast umber of beats.The result shows, and is effective equally based on the DOA algorithm for estimating that subarray is synthetic.
Comparison diagram 5 and Fig. 3 as can be known, it is influential to DOA resolution that signal reaches wave angle.Reach wave angle and increase, DOA resolution will reduce.
Fig. 2 to Fig. 5 shows that disclosed in this invention is effective based on subarray synthetic arrived wave signal treatment system and signal processing method.
Although disclose most preferred embodiment of the present invention and accompanying drawing for the purpose of illustration, it will be appreciated by those skilled in the art that: without departing from the spirit and scope of the invention and the appended claims, various replacements, variation and modification all are possible.Therefore, the present invention should not be limited to most preferred embodiment and the disclosed content of accompanying drawing.
Claims (9)
1. an array antenna signal processing system based on the synthetic arrived wave signal direction determining method of subarray, realizes that receiver quantity is less than the arrived wave signal direction estimation under the bay quantity situation, and this treatment system comprises:
The array antenna part, the bay induction receives the far field arrived wave signal;
The switching matrix part is according to the instruction connection or the disconnected node of control circuit;
The signal receiver part is amplified or decay bay amplitude output signal, carries out AGC and handles, and carries out demodulation, decoding is transformed to digital baseband signal with arrived wave signal;
Control section, the connection or the disconnection of control switching matrix, the accessing operation of control data memory;
The covariance calculating section, the covariance of calculating delayed output signals;
Subarray composition part, synthetic equivalent covariance;
Data memory part, storing intermediate result of data operation and final result;
The space characteristics test section, the space characteristic parameter of detection arrived wave signal is realized the estimation of arrived wave signal direction.
2. according to claim 1, it is characterized in that: to the array antenna type without limits, both can be uniform line-array, also can be the linear array of irregular alignment, or the array of plane distribution.。
3. according to claim 1, it is characterized in that: the quantity of bay is not strict with, and array antenna is minimum can to comprise 2 array elements, array element quantity can by reach the wave number amount, processing speed decides.
4. according to claim 1, it is characterized in that: receiver quantity also can equal bay quantity.
5. according to claim 1, it is characterized in that: the array element of subarray antenna does not directly link to each other with receiver, but by switching matrix, adopts time division multiplexing mode to be connected with the input of receiver.
6. realize reaching ripple direction super-resolution estimated signals processing method based on the arrived wave signal snap for one kind, described method comprises:
According to the array element numbering, main array antenna is divided into the subarray antenna;
Based on the time division multiplexing processing method, ask the covariance matrix of sub antenna array output;
Ask the equivalent covariance matrix of main antenna array output;
Carrying out the arrived wave signal direction estimates;
7. according to claim 6, it is characterized in that: receiver is implemented sampling to the data of the whole array element outputs of sub-array antenna constantly in each snap sampling.
8. according to claim 6, it is characterized in that: by time-sharing multiplex, can receive the fast beat of data of output of handling any amount bay, realize reception and processing a large amount of delayed output signals.
9. the method for sub-array antenna is chosen in an optimization, be divided into the subarray antenna according to described method master's array antenna, the array number of each subarray antenna is less than the array number of main array antenna, the array number of each subarray antenna equals the quantity of signal receiver, the array element of subarray antenna does not directly link to each other with receiver, but by switching matrix, adopts time division multiplexing mode to be connected with the input of receiver, it is characterized in that the following principle of the employing of division:
Export fast beat of data by whole subarrays and can synthesize the array antenna output correlation matrix of exporting the correlation matrix equivalence with main array, be called equivalent main array output correlation matrix;
Any two array elements of main array antenna occur once in some subarray antennas at least simultaneously;
The minimum number of subarray antenna.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100096797A CN100373808C (en) | 2004-10-15 | 2004-10-15 | Signal processing system and method based on estimating signal direction of anived wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100096797A CN100373808C (en) | 2004-10-15 | 2004-10-15 | Signal processing system and method based on estimating signal direction of anived wave |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1761171A true CN1761171A (en) | 2006-04-19 |
CN100373808C CN100373808C (en) | 2008-03-05 |
Family
ID=36707137
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100096797A Expired - Fee Related CN100373808C (en) | 2004-10-15 | 2004-10-15 | Signal processing system and method based on estimating signal direction of anived wave |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100373808C (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103901394A (en) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | Quantum universal gravitation searching dynamic DOA estimation method under impulsive noise environment |
CN108430090A (en) * | 2017-02-15 | 2018-08-21 | 华为技术有限公司 | The determining and method for determining DOA and access network equipment and terminal |
CN109061620A (en) * | 2013-05-24 | 2018-12-21 | 罗伯特·博世有限公司 | Method for running MIMO radar |
CN110515033A (en) * | 2019-08-09 | 2019-11-29 | 南京航空航天大学 | A kind of deficient channel direction-finding system and method restored based on Toeplitz matrix |
CN111650553A (en) * | 2020-06-02 | 2020-09-11 | 斯凯瑞利(北京)科技有限公司 | Signal processing system and method for time division multiplexing-based direction estimation of arriving signals |
CN111830460A (en) * | 2020-07-23 | 2020-10-27 | 哈尔滨工业大学(威海) | DOA estimation method based on sequential MUSIC |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5262789A (en) * | 1992-04-30 | 1993-11-16 | General Electric Company | Source identification system for closely separated spatial sources |
FR2805614B1 (en) * | 2000-02-25 | 2003-08-22 | Thomson Csf | METHOD FOR LOCATING RADIO SOURCES USING A TWO-CHANNEL HIGH RESOLUTION RADIOGONIOMETER |
US6567034B1 (en) * | 2001-09-05 | 2003-05-20 | Lockheed Martin Corporation | Digital beamforming radar system and method with super-resolution multiple jammer location |
JP4098026B2 (en) * | 2002-07-31 | 2008-06-11 | 富士通株式会社 | Method of estimating direction of arrival of periodic stationary signal in multipath propagation environment and reception beam forming apparatus using the same |
CN1297822C (en) * | 2003-02-21 | 2007-01-31 | 重庆邮电学院 | Estimation method for radio orientation incoming wave direction based on TD-SCMA |
-
2004
- 2004-10-15 CN CNB2004100096797A patent/CN100373808C/en not_active Expired - Fee Related
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085566B (en) * | 2013-05-24 | 2023-02-28 | 罗伯特·博世有限公司 | Method for operating a multiple-input multiple-output radar |
CN109061620A (en) * | 2013-05-24 | 2018-12-21 | 罗伯特·博世有限公司 | Method for running MIMO radar |
CN109085566A (en) * | 2013-05-24 | 2018-12-25 | 罗伯特·博世有限公司 | Method for running MIMO radar |
CN103901394B (en) * | 2014-03-21 | 2016-07-06 | 哈尔滨工程大学 | Quantum gravitation under a kind of impulsive noise environment searches for dynamic DOA estimation method |
CN103901394A (en) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | Quantum universal gravitation searching dynamic DOA estimation method under impulsive noise environment |
CN108430090A (en) * | 2017-02-15 | 2018-08-21 | 华为技术有限公司 | The determining and method for determining DOA and access network equipment and terminal |
WO2018149346A1 (en) * | 2017-02-15 | 2018-08-23 | 华为技术有限公司 | Doa determining method, method for use in doa determining, access network device, and terminal |
CN108430090B (en) * | 2017-02-15 | 2020-07-28 | 华为技术有限公司 | Method for determining and determining DOA, access network equipment and terminal |
CN110515033A (en) * | 2019-08-09 | 2019-11-29 | 南京航空航天大学 | A kind of deficient channel direction-finding system and method restored based on Toeplitz matrix |
CN110515033B (en) * | 2019-08-09 | 2022-06-28 | 南京航空航天大学 | Toeplitz matrix recovery-based under-channel direction finding system and method |
CN111650553A (en) * | 2020-06-02 | 2020-09-11 | 斯凯瑞利(北京)科技有限公司 | Signal processing system and method for time division multiplexing-based direction estimation of arriving signals |
CN111830460A (en) * | 2020-07-23 | 2020-10-27 | 哈尔滨工业大学(威海) | DOA estimation method based on sequential MUSIC |
CN111830460B (en) * | 2020-07-23 | 2022-11-01 | 哈尔滨工业大学(威海) | DOA estimation method based on sequential MUSIC |
Also Published As
Publication number | Publication date |
---|---|
CN100373808C (en) | 2008-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112350759B (en) | RIS-assisted multi-user multi-antenna communication and radar spectrum sharing method | |
Yu et al. | Constrained waveform design for colocated MIMO radar with uncertain steering matrices | |
CN1106086C (en) | Signal processing method utilizing eigenvector corresponding to maximum eigenvalue of autocorrelation matrix of received signals for antenna array system | |
CN1835416A (en) | Method and apparatus for direction-of-arrival tracking | |
CN1187924C (en) | Intelligence antenna for IMI-2000 CDMA | |
CN1194494C (en) | Adaptive antenna receiving equipment | |
CN1211755C (en) | System and method for adaptive filtering | |
CN1808949A (en) | Non-physical modeling and emulation method for channels in multi-input and multi-output communication system | |
CN1445885A (en) | Directional diagram forming method and device for self-adaptive antenna array in base station | |
KR101199569B1 (en) | Opportunistic interference aligned user selection in multiuser mimo interference channels | |
CN1384680A (en) | Signal receiving method and device | |
KR20050004605A (en) | Combined beamforming-diversity wireless fading channel de-modulator using sub-array grouped adaptive array antennas, portable telecommunication receiving system comprising it and method thereof | |
CN1878026A (en) | Intelligent antenna implementing method and system based on software radio | |
KR101088602B1 (en) | Channel information prediction system and channel information prediction method | |
JP2007159130A (en) | Uplink receiving method and device in distributed antenna mobile communication system | |
CN115085774B (en) | Common sense fusion hybrid beamforming method based on Keramelteon | |
KR20100107063A (en) | Channel information prediction system and channel information prediction method | |
Deng et al. | Co-design of waveform correlation matrix and antenna positions for MIMO radar transmit beampattern formation | |
CN1761171A (en) | Signal processing system and method based on estimating signal direction of anived wave | |
CN1802771A (en) | Adaptive antenna reception method and device | |
CN1819574A (en) | Carrier frequency bias estimation with OFDMA up link system intersection | |
Wang et al. | Average rate maximization for mobile STAR-RIS-aided NOMA system | |
Liu et al. | An efficient method to determine the diagonal loading factor using the constant modulus feature | |
CN1163015C (en) | Multi-user disturbance eliminator | |
CN117498900A (en) | Resource allocation device and method for honeycomb-removing large-scale MIMO (multiple input multiple output) general sense integrated system |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080305 Termination date: 20111015 |