CN109120305B - A kind of digital tracking method of broadband beams, system and equipment - Google Patents
A kind of digital tracking method of broadband beams, system and equipment Download PDFInfo
- Publication number
- CN109120305B CN109120305B CN201810724910.2A CN201810724910A CN109120305B CN 109120305 B CN109120305 B CN 109120305B CN 201810724910 A CN201810724910 A CN 201810724910A CN 109120305 B CN109120305 B CN 109120305B
- Authority
- CN
- China
- Prior art keywords
- maximum value
- beam synthesis
- arrival bearing
- value
- parameter
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7075—Synchronisation aspects with code phase acquisition
- H04B1/70754—Setting of search window, i.e. range of code offsets to be searched
Abstract
The invention discloses a kind of digital tracking methods of broadband beams, comprising: carries out the calculating of frequency domain Beam synthesis to each circuit-switched data and arrival bearing's vector of finishing FFT processing, obtains Beam synthesis amount；Obtain the maximum value of Beam synthesis amount；Signal detection and normalized are done to the maximum value of Beam synthesis amount, obtain parameter value；When parameter value is greater than or equal to preset useful signal reference value, and the maximum value of Beam synthesis amount is less than or equal to the useful signal reference value, determine that there are useful signals in current this subchannel of minimum basis；The maximum value of parameter value is obtained, and the final bearing and final pitch angle of wave beam are obtained according to the maximum value of parameter value.The embodiment of the invention also discloses the figure tracking equipment of a kind of digital servomechanism of broadband beams and a kind of broadband beams.Using the embodiment of the present invention, quickly estimate the arrival bearing of current time satellite-signal, achieve the effect that phased array antenna it is real-time to broadband beams, it is quick, accurately track.
Description
Technical field
The present invention relates to technical field of satellite communication more particularly to a kind of digital tracking method of broadband beams, system and
Equipment.
Background technique
Satellite communication has the characteristics that covering is wide, transmission is fast, capacity is not limited by landform greatly, receives in recent years extensively
Concern.Now, as commercial market autonomous satellite phone is formally distributed telephone numbers, satellite communication is to spread comprehensively in dual-use field
It opens.Related with wave frequency is received due to arranging spacing between each array element of array antenna, frequency is higher, between the arrangement needed
Away from smaller, thus for the electromagnetic wave of VHF band, array antenna spacing needs very small, this is to the technique of hardware and heat dissipation etc.
It is proposed very high requirement.
In current satellite signal receiving system, the common satellite earth antenna of fixed station is parabola antenna, movement station
The mode that satellite antenna is often used mechanical scanning tracks satellite received signal.Mechanical scanning satellite antenna is generally side
Position, the driving of two axis of pitching carry out real-time control to mechanical device using mechanical servo.Mechanical scanning antennas according to sweeping in real time
The reception level retouched, to judge azimuth and the pitch angle of satellite-signal, so that it is right in real time by way of machinery rotation
Satelloid sense.But existing mechanical scanning antennas is for the ICBM SHF satellite terminal of high-speed mobile that (mountainous region high speed is mobile
Automobile, train, aircraft etc.), since its automatic reaction time is slower, causes it that cannot reach and real-time, quickly and accurately trace into
Signal.
Summary of the invention
The purpose of the embodiment of the present invention is that a kind of digital tracking method of broadband beams, system and equipment are provided, it can be quickly
Estimate the arrival bearing of current time satellite-signal, with reach phased array antenna it is real-time to broadband beams, it is quick, accurately track
Effect.
To achieve the above object, the embodiment of the invention provides a kind of digital tracking methods of broadband beams, comprising:
Obtain arrival bearing's vector of two-dimensional planar array；
The data for obtaining ADC acquisition, do FFT processing to each circuit-switched data respectively；
The calculating of frequency domain Beam synthesis is carried out to each circuit-switched data and arrival bearing's vector of finishing FFT processing, is obtained
Beam synthesis amount；
The maximum value of the Beam synthesis amount described in same this subchannel of minimum basis of acquisition；
Signal detection and normalized are done to the maximum value of the Beam synthesis amount, obtain parameter value；
When the parameter value is greater than or equal to preset useful signal reference value, and the maximum value of the Beam synthesis amount is small
When the useful signal reference value, determine that there are useful signals in current this subchannel of minimum basis；
When judgement is there are the maximum value of the parameter value when useful signal, is obtained, the maximum value of the parameter value is corresponding
Scanning space is final arrival bearing, and according to the final arrival bearing obtain wave beam final bearing and final pitching
Angle.
Compared with prior art, the digital tracking method of broadband beams disclosed by the invention, first by finishing FFT
Each circuit-switched data of processing and arrival bearing's vector carry out the calculating of frequency domain Beam synthesis, obtain Beam synthesis amount；Then right
The maximum value of the Beam synthesis amount does signal detection and normalized, obtains parameter value；When the parameter value is greater than or waits
In preset useful signal reference value, and the maximum value of the Beam synthesis amount is less than or equal to the useful signal reference value
When, determine that there are useful signals in current this subchannel of minimum basis；The maximum value of the parameter value is finally obtained, and according to described
The maximum value of parameter value obtains the final bearing and final pitch angle of wave beam.Solves the existing mechanical scanning day of the prior art
Line causes it that cannot reach real-time, quick, accurate the ICBM SHF satellite terminal of high-speed mobile since its automatic reaction time is slower
Ground traces into the problem of signal.The arrival bearing of current time satellite-signal can be estimated, quickly to reach phased array antenna
Effect that is real-time to broadband beams, quick, accurately tracking.
As an improvement of the above scheme, the arrival bearing's vector for obtaining two-dimensional planar array, comprising:
According to the azimuth and pitch angle of the array manifold of bay number composition and reception wave beam, it is flat to obtain two dimension
Arrival bearing's vector of face array.
It is as an improvement of the above scheme, described to do FFT processing respectively to per data all the way, comprising:
FFT operation is done respectively to each circuit-switched data；
An operation is taken to finish that each circuit-switched data after the FFT operation takes predetermined number point.
As an improvement of the above scheme, acquisition Beam synthesis amount described in same this subchannel of minimum basis is most
Big value, comprising:
Signal search and the judgement of two dimensional beam direction are carried out to the Beam synthesis amount, obtained in the same minimum basis book
The maximum value of Beam synthesis amount described in channel.
As an improvement of the above scheme, before the arrival bearing's vector for obtaining two-dimensional planar array further include:
Obtain conventional parameter, frequency range including scanning, number, the ADC sampling rate, FFT of minimum basis this subchannel
The points of processing, the frequency of arrival bearing and useful signal reference value；
The parameter for obtaining phase battle array control aerial array, bows including bay number, the dividing condition of attitude, space
The dividing condition at angle and the dividing condition of scanning space.
To achieve the above object, the embodiment of the invention also provides a kind of digital servomechanisms of broadband beams, comprising:
Arrival bearing's vector acquiring unit, for obtaining arrival bearing's vector of two-dimensional planar array；
Fft processing unit does FFT processing to each circuit-switched data for obtaining the data of ADC acquisition respectively；
Frequency domain Beam synthesis computing unit, for each circuit-switched data and arrival bearing's vector for finishing FFT processing
The calculating of frequency domain Beam synthesis is carried out, Beam synthesis amount is obtained；
Beam synthesis amount maximum value acquiring unit, for obtaining the Beam synthesis described in same this subchannel of minimum basis
The maximum value of amount；
Signal detection and normalized unit do signal detection and normalizing for the maximum value to the Beam synthesis amount
Change processing, obtains parameter value；
Useful signal judging unit, for being greater than or equal to preset useful signal reference value, and institute when the parameter value
When stating the maximum value of Beam synthesis amount less than or equal to the useful signal reference value, determine to deposit in current this subchannel of minimum basis
In useful signal；
Final result acquiring unit, it is described for when determining there are when useful signal, obtaining the maximum value of the parameter value
The corresponding scanning space of the maximum value of parameter value is final arrival bearing, and obtains wave beam according to the final arrival bearing
Final bearing and final pitch angle.
Compared with prior art, the digital servomechanism of broadband beams disclosed by the invention is closed by frequency domain wave beam first
The calculating of frequency domain Beam synthesis is carried out to each circuit-switched data and arrival bearing's vector of finishing FFT processing at computing unit, is obtained
To Beam synthesis amount；Then signal detection and normalized unit to the maximum value of the Beam synthesis amount do signal detection and
Normalized obtains parameter value；When the parameter value is greater than or equal to preset useful signal reference value, and wave beam conjunction
When being less than or equal to the useful signal reference value at the maximum value of amount, useful signal judging unit determines current minimum basis book
There are useful signals in channel；Last final result acquiring unit obtains the maximum value of the parameter value, and according to the parameter
The maximum value of value obtains the final bearing and final pitch angle of wave beam.Solves the existing mechanical scanning antennas pair of the prior art
In the ICBM SHF satellite terminal of high-speed mobile, since its automatic reaction time is slower, cause its cannot reach real-time, quickly and accurately with
The problem of track is to signal.The arrival bearing of current time satellite-signal can be estimated, quickly to reach phased array antenna to width
Effect that is real-time with wave beam, quick, accurately tracking.
As an improvement of the above scheme, arrival bearing's vector acquiring unit is used to be formed according to bay number
Array manifold and azimuth and the pitch angle for receiving wave beam, obtain arrival bearing's vector of two-dimensional planar array.
As an improvement of the above scheme, the fft processing unit is used to do FFT operation respectively to each circuit-switched data, and right
What each circuit-switched data after finishing the FFT operation took predetermined number point takes operation.
As an improvement of the above scheme, the Beam synthesis amount maximum value acquiring unit be used for the Beam synthesis amount into
Row signal search and the judgement of two dimensional beam direction, obtain the maximum of the Beam synthesis amount described in same this subchannel of minimum basis
Value；
The system also includes:
Parameter acquiring unit, for obtaining the parameter of conventional parameter and phase battle array control aerial array；Wherein, the conventional parameter
Frequency range including scanning, the number of minimum basis this subchannel, ADC sampling rate, the points of FFT processing, arrival bearing
The frequency and useful signal reference value；The parameter of the phase battle array control aerial array includes stroke of bay number, attitude
Divide situation, the dividing condition of the dividing condition at the space angle of depression and scanning space.
To achieve the above object, the embodiment of the present invention also provides a kind of figure tracking equipment of broadband beams, including processing
Device, memory and storage in the memory and are configured as the computer program executed by the processor, the place
Reason device realizes the digital tracking method of broadband beams described in any of the above-described embodiment when executing the computer program.
Detailed description of the invention
Fig. 1 is a kind of flow chart of the digital tracking method of broadband beams provided in an embodiment of the present invention；
Fig. 2 is a kind of another flow chart of the digital tracking method of broadband beams provided in an embodiment of the present invention；
Fig. 3 is a kind of structural block diagram of the digital servomechanism 10 of broadband beams provided in an embodiment of the present invention；
Fig. 4 is a kind of structural block diagram of the figure tracking equipment 20 of broadband beams provided in an embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Embodiment one
It is a kind of flow chart of the digital tracking method of broadband beams provided in an embodiment of the present invention referring to Fig. 1, Fig. 1；Packet
It includes:
S11, the arrival bearing's vector for obtaining two-dimensional planar array；
S12, the data for obtaining ADC acquisition, do FFT processing to each circuit-switched data respectively；
S13, the calculating of frequency domain Beam synthesis is carried out to each circuit-switched data and arrival bearing's vector of finishing FFT processing,
Obtain Beam synthesis amount；
S14, the maximum value for obtaining the Beam synthesis amount described in same this subchannel of minimum basis；
S15, signal detection and normalized are done to the maximum value of the Beam synthesis amount, obtains parameter value；
S16, it is greater than or equal to preset useful signal reference value, and the maximum of the Beam synthesis amount when the parameter value
When value is less than or equal to the useful signal reference value, determine that there are useful signals in current this subchannel of minimum basis；
S17, when determine there are when useful signal, obtain the maximum value of the parameter value, the maximum value pair of the parameter value
The scanning space answered is final arrival bearing, and obtains the final bearing of wave beam and final according to the final arrival bearing
Pitch angle.
Preferably, before executing step S11, step S10 is first carried out: obtaining conventional parameter and phase battle array control aerial array
Parameter.Wherein, the conventional parameter include the frequency range of scanning, the number of minimum basis this subchannel, ADC sampling rate,
Points, the frequency of arrival bearing and the useful signal reference value of FFT processing；The parameter of the phase battle array control aerial array includes antenna
Element number of array, the dividing condition of attitude, the dividing condition of the dividing condition at the space angle of depression and scanning space.
Specifically, the frequency range of the system of reception scanning is S frequency range, then the frequency range of the scanning is f1~f2 (Hz), frequency
Band difference is denoted as fd=f2-f1.Preferably, f1 is typically chosen 10MHz, and f2 is typically chosen 30MHz；S frequency range refers generally to frequency model
It is trapped among the electromagnetic wave band limits of 1.55-3.4GHz, is mapped to the array element spacing of array antenna in decimetre rank, this is for the modern times
Technique and heat dissipation are not all problem, thus are come into being for the various applications of S frequency range array antenna.It is fast in order to solve ground
Fast mobile terminal is the broadband beams tracking problem and antenna tracking of S band satellite communication are slow, height is high, volume is big, weight
Weight, power consumption be high, cannot conformal design the problem of, need the number using S band satellite communication two dimensional phased array antenna broadband beams
Word tracking designs the internal S band satellite communication two dimensional phased array antenna for using Digital Signal Processing, to realize S frequency range
Satellite communication broadband beams quick and precisely track target.
Specifically, the bandwidth of this subchannel of minimum basis is Bw (Hz) in current reception system, then entire scanning band can wrap
The number of described this subchannel of minimum basis contained meets Nbas=fd/Bw.Common this subchannel bandwidth of reception system minimum basis
Generally 2.4KHz, namely can satisfy the transmission of speech all the way, certainly, it can choose as needed smaller or larger.With ground
For the Cellular Networks of face, the bandwidth of a cellular cell may be several million bandwidth, however each channel is not to have taken entire bee
The bandwidth of nest cell, and each channel only may occupy the bandwidth of several K, satellite is also the band in this way, satellite beams
It is wide very wide, but apply channel that can not all take bandwidth, at this moment just will appear extreme case is exactly a satellite beams band
It is wide by tens, sometime may only have several K datas to transmit above, and other bandwidth are idle state entirely, the present invention is implemented
The method that example proposes also can fast and accurately trace into wave beam in the state that satellite bandwidth is almost comprehensively idle.
Specifically, the ADC sampling rate of the system of reception meets fs=2xfd.The base band that reception system samples ADC
The points that signal does FFT processing are NF point, and NF generally chooses 2 power side here, general often selection 8192 or 16384.At this point,
It defines each this subchannel of minimum basis to be made of NB point after FFT transform, meets NB=Bw/ (fs/NF), ensure that NB is most here
Amount is close to integer, to reduce the error of subsequent calculating.
Specifically, the frequency that the system of reception obtains the arrival bearing is f_refresh (Hz), namely every T_
Arrival bearing estimation of starting in refresh=1/f_refresh seconds, this parameter can sets itself according to actual needs, generally
Per second refresh 50~150 times of car terminals setting in mountain area is sufficient for communication performance.
Specifically, the useful signal reference value is obtained, wherein the useful signal reference value is εmin~εmax, εminFor
The lower limit of the useful signal reference value, εmaxFor the upper limit of the useful signal reference value, the two are fixed values, can be artificial
It need to be set according to practical reception system.
Specifically, the bay number of phased-array antenna array is M, the bay number M is generally chosen
Tens.
Specifically, according to actual needs, the division number for receiving system to 360 degree of attitude is Nazi, to space
90 degree of pitch angle of division number is Npit.Here attitude and space pitch angle can be divided equidistantly, it can also not
It is equidistant to divide, sets itself can be needed according to reception system.Preferably, the value range of Nazi is 30~60 degree, Npit's
Value range is 10~20 degree.In addition, according to different regions service condition difference, such as region of the equator uses, basic satellite is just
In 90 degree of positions of pitch angle, polar regions, satellite, can be according to actual use region, weight just in the position that pitch angle is almost 0 degree
Point divides common position.The dividing condition of the scanning space meets Nspa=Nazi × Npit.
Specifically, in step s 11, the array manifold formed according to M array element and the reception all possible side of wave beam
Parallactic angle and pitch angle, precalculate or microwave dark room measurement in advance obtains arrival bearing's arrow of Nspa M rank two-dimensional planar array
Amount, is denoted as Al(i), wherein I=1 ..., Nspa；I=1 ..., M
It is worth noting that selecting different array manifolds, the direction vector of two-dimensional surface battle array is different, such as plane
The direction vector of circle battle array meets:
Wherein i=1 ... M, θlFor the pitch angle of the i-th space block divided,For the i-th space block that has divided
Azimuth, R are the radius of circular array, and λ is wavelength, and j is pure imaginary number, if becoming square matrix, direction vector can change again, are needed
Recalculate formula.But it has to be noticed that above formula is all theoretical calculation formula, due to actual hardware material and technique
Difference, cause array antenna anti-line vector can with theoretical value gap, so preferably being gone on the basis of theoretical calculation micro-
Wave darkroom carries out correction in kind, to obtain optimal direction vector.
Specifically, in step s 12, timer is detected first, when the pulse that timer was generated every T_refresh seconds is arrived
When coming, the data of current time M road ADC acquisition are saved, are NF point per the points for saving base band data all the way.Since ADC is one
Directly work, thus ADC can output data always, it is every when being not per but due to the limitation of antenna beam tracking refreshing frequency
Quarter is all calculated in the data acquired using ADC, is only taken the data that ADC samples out in fixed time, so when most of
It carves, the data of ADC acquisition are that discarding is unused.Preferably, general ground-plane antenna selection takes an ADC every 10ms to 50ms
The data sampled out are proper.
Then, the FFT operation of NF point is done respectively to the base band sample data on the road M, it is assumed that the road M base band sample data is xm
(n), then the base band sample data after FFT operation meets:
Xm(K)=FFT (xm(n)) formula (2)；
Wherein, XmIt (K) is the base band sample data after FFT operation；M=1 ..., M；N=1 ..., NF；K=1 ..., NF.
It is worth noting that the spectral range that the result after above-mentioned FFT operation indicates is [- fs/2, fs/2], the general (scene FPGA
Programmable gate array) result that the IP kernel of FFT operation calculates is done in the inside is not this range, so should be noted that when realizing pair
Operation result is corrected.
Finally, to the transformed road M Xm(K) operation is taken to obtain Xm(K '), every road only takes NF/2 point, described default
Number point is NF/2, and just corresponding frequency spectrum is exactly the signal spectrum for needing to scan to these points on every road, specifically takes a formula
It is as follows:
Xm(K ')=Xm(K '+NF/4) formula (3)；
Wherein, m=1 ..., M；K '=1 ..., NF/2.Since the echo signal bandwidth to be detected only has fd (Hz), and
The expression range of signal is [- fd, fd] after FFT operation, thus takes Xm(K) NF/2 point come obtain echo signal frequency spectrum [-
Fd/2, fd/2].
Specifically, in step s 13, by Xm(K ') and arrival bearing's vector Al(i) frequency domain Beam synthesis meter is carried out
It calculates, specific formula for calculation is as follows:
Wherein, Y (p, l) is Beam synthesis amount；P=1 ..., Nbas；I=1 ..., Nspa.Frequency domain wave beam in formula (4)
Synthesis is exactly with the X receivedmArrival bearing's vector A of (K ') signal and divided all spacesl(i) traversal synthesis one
Time, namely exhaustion retrieval.
Specifically, carrying out signal search in step S14 to the Beam synthesis amount Y (p, l) and two dimensional beam direction being added to sentence
Certainly, the maximum value in Nbas Y (p, l) for belonging to same this subchannel of minimum basis is found first here, maximum value Ymax
(p) it indicates, the maximum value Y of the Beam synthesis amountmax(p) corresponding I (number of scanning space) uses lmaxIt indicates.This operation is
It is that most possible arrival bearing is found in all possible arrival bearing, the value of the Beam synthesis amount Y (p, l) is bigger,
Then show that this segment data signal is most likely to be to come from I direction.
Specifically, the maximum value Y in step S15, first to the Beam synthesis amountmax(p) signal detection is done, it is described
The process of signal detection includes: will be in the same p and Ymax(p) all Beam synthesis amount Y (p, l) of 1dB are theoretically differed
It adds up and is denoted as Ymax1dB(p), Ymax1dBIt (p) is signal energy detection limit, the gross energy of all Y (p, l) is denoted as in the same p
Ysum(p).Due to FFT operation signal energy can be made to reveal, thus the point of maximum value can not representation signal whole energy, this
In selection and maximum value difference 1dB size signal, be construed as the energy part of signal leakage, thus it mutually added up
As the detection limit of signal energy, and the signal of entire frequency spectrum is added the reduced value as signal detection amount, seeks the two
Variable is to do normalized to signal energy detection limit for next, the benefit of normalized be exactly detection threshold not by
The influence of ADC sampled signal amplitude variation.
Then, to Ymax1dB(p) normalized is done, specific formula is as follows:
Wherein, p=1 ..., Nbas；R (p) is to signal energy detection limit Ymax1dB(p) normalized done as a result,
That is parameter value.
Specifically, in step s 16, judging whether the parameter value R (p) is greater than or equal to the useful signal reference value
Lower limit εmin, while judging the maximum value Y of the Beam synthesis amountmax(p) whether it is less than or equal to the useful signal to refer to
The upper limit ε of valuemax。εminAnd εmaxFor judging in this subchannel bandwidth of p-th of current time minimum basis with the presence or absence of effectively letter
Number.
If R (p) >=εminAnd Ymax(p)≤εmax, then it is assumed that p-th of minimum basis this subchannel bandwidth memory is in useful signal.
If R (p) < εmin, then it is assumed that without useful signal in this subchannel bandwidth of p-th of current time minimum basis；If Ymax
(p) > εmax, then it is assumed that there is interference signal in this subchannel bandwidth of p-th of current time minimum basis, joined interference processing herein
Partial judgement, so that method of the invention has certain anti-interference ability in practical applications.As R (p) < εminOr Ymax
(p) > εmaxWhen, the last moment calculated azimuth and pitch angle of output storage storage, at the same return step S12, weight
FFT processing newly is done to the data of the road M ADC acquisition, waits the arrival bearing of subsequent time to estimate, to reach target frequency bands broadband
Wave beam in real time, quickly, the effect that accurately tracks.The estimation frequency can be improved to reach intimate in the terminal relatively high for mobility
The effect of uninterrupted real-time estimation beam direction.
Specifically, in step S17, when determining that there are when useful signal, find all useful signals that are determined to have
Maximum value R in the parameter value R (p)max(p), Rmax(p)=maxP=1 ..., Nbas{R(p)}.The maximum value of the parameter value
Rmax(p) corresponding lmaxFor the final arrival bearing.The above process is to find signal from all this subchannels of minimum basis
The strongest sub-channels of energy, using the direction of the sub-channel signal as the final beam direction of output.
Specifically, according to lmaxThe final bearing and final pitch angle of the final arrival bearing corresponding can be found,
This moment calculated final bearing and the final pitch angle are stored in memory, and export the final orientation
Angle and the final pitch angle.Each scanning space I can correspond to an azimuth and pitch angle, by tabling look-up, can look for from I
To its corresponding azimuth and pitch angle.Certainly also the azimuth successfully estimated every time and pitch angle moment are updated here
It stores, to prevent sometime because of the factors such as signal interference, use when cannot accurately estimate azimuth and pitch angle.
Preferably, the process of above-mentioned steps S11~S17 can also refer to Fig. 2.It is worth noting that when Fig. 2 indicates current
It carves to the figure tracking processes of broadband beams, on final bearing and the pitch angle, or output described in the output current time
Behind the azimuth at one moment and pitch angle, start the azimuth for recalculating subsequent time and pitch angle, at this point, return step
S12 is recalculated.
When it is implemented, first by carrying out frequency to each circuit-switched data and arrival bearing's vector of finishing FFT processing
Domain Beam synthesis calculates, and obtains Beam synthesis amount；Then signal detection and normalization are done to the maximum value of the Beam synthesis amount
Processing, obtains parameter value；When the parameter value is greater than or equal to preset useful signal reference value, and the Beam synthesis amount
When maximum value is less than or equal to the useful signal reference value, determine that there are useful signals in current this subchannel of minimum basis；Most
It obtains the maximum value of the parameter value afterwards, and the final bearing of wave beam is obtained according to the maximum value of the parameter value and is finally bowed
The elevation angle.
Compared with prior art, the digital tracking method of broadband beams disclosed by the invention, it is existing to solve the prior art
Mechanical scanning antennas cause it that cannot reach real since its automatic reaction time is slower the ICBM SHF satellite terminal of high-speed mobile
When, quickly and accurately trace into the problem of signal.The arrival bearing of current time satellite-signal can be estimated, quickly to reach
Effect that is real-time to broadband beams, quick to phased array antenna, accurately tracking.
The present invention overcomes general antenna arrival bearing using digital Mutual coupling technology and estimates slower problem, has
The advantages of quickly estimating and tracking arrival bearing.The problem of also overcoming the frequent losing lock of general antenna tracking satellite beams simultaneously,
Has the advantages that tenacious tracking beam direction.The technology for refreshing the estimation frequency is set using number, and overcoming general antenna cannot
The problem of adapting to high dynamic tracking terminal satellite beams, and then have the advantages that adapt to various terminals tracking satellite wave beam.Using
The technology of digital signal detection amount normalized, overcome general antenna cannot adaptive adjustment signal detection threshold ask
Topic has the advantages that the wave beam tracking for adapting to various satellite-signal sizes.The technology detected using digital signal upper energy limit, gram
The problem of general antenna cannot exclude interference signal has been taken, has had the advantages that certain anti-interference ability.
Embodiment two
It is a kind of structural frames of the digital servomechanism 10 of broadband beams provided in an embodiment of the present invention referring to Fig. 3, Fig. 3
Figure；Include:
Arrival bearing's vector acquiring unit 11, for obtaining arrival bearing's vector of two-dimensional planar array；
Fft processing unit 12 does FFT processing to each circuit-switched data for obtaining the data of ADC acquisition respectively；
Frequency domain Beam synthesis computing unit 13, for being sweared to each circuit-switched data and the arrival bearing of finishing FFT processing
Amount carries out the calculating of frequency domain Beam synthesis, obtains Beam synthesis amount；
Beam synthesis amount maximum value acquiring unit 14 is closed for obtaining the wave beam described in same this subchannel of minimum basis
At the maximum value of amount；
Signal detection and normalized unit 15 do signal detection for the maximum value to the Beam synthesis amount and return
One change processing, obtains parameter value；
Useful signal judging unit 16, for being greater than or equal to preset useful signal reference value when the parameter value, and
When the maximum value of the Beam synthesis amount is less than or equal to the useful signal reference value, determine in current this subchannel of minimum basis
There are useful signals；
Final result acquiring unit 17, for when there are when useful signal, obtain the maximum value of the parameter value, institute for judgement
The corresponding scanning space of maximum value for stating parameter value is final arrival bearing, and obtains wave beam according to the final arrival bearing
Final bearing and final pitch angle.
Preferably, the digital servomechanism 10 of the broadband beams further includes parameter acquiring unit 18, and the parameter obtains
Unit 18 is used to obtain the parameter of conventional parameter and phase battle array control aerial array.Wherein, the conventional parameter includes the frequency of scanning
Range, the number of minimum basis this subchannel, ADC sampling rate, FFT processing points, the frequency of arrival bearing and useful signal
Reference value；The parameter of the phase battle array control aerial array includes bay number, the dividing condition of attitude, the space angle of depression
Dividing condition and scanning space dividing condition.
Specifically, the frequency range of the system of reception scanning is S frequency range, then the frequency range of the scanning is f1~f2 (Hz), frequency
Band difference is denoted as fd=f2-f1.Preferably, f1 is typically chosen 10MHz, and f2 is typically chosen 30MHz；S frequency range refers generally to frequency model
It is trapped among the electromagnetic wave band limits of 1.55-3.4GHz, is mapped to the array element spacing of array antenna in decimetre rank, this is for the modern times
Technique and heat dissipation are not all problem, thus are come into being for the various applications of S frequency range array antenna.It is fast in order to solve ground
Fast mobile terminal is the broadband beams tracking problem and antenna tracking of S band satellite communication are slow, height is high, volume is big, weight
Weight, power consumption be high, cannot conformal design the problem of, need the number using S band satellite communication two dimensional phased array antenna broadband beams
Word tracking designs the internal S band satellite communication two dimensional phased array antenna for using Digital Signal Processing, to realize S frequency range
Satellite communication broadband beams quick and precisely track target.
Specifically, the bandwidth of this subchannel of minimum basis is Bw (Hz) in current reception system, then entire scanning band can wrap
The number of described this subchannel of minimum basis contained meets Nbas=fd/Bw.Common this subchannel bandwidth of reception system minimum basis
Generally 2.4KHz, namely can satisfy the transmission of speech all the way, certainly, it can choose as needed smaller or larger.With ground
For the Cellular Networks of face, the bandwidth of a cellular cell may be several million bandwidth, however each channel is not to have taken entire bee
The bandwidth of nest cell, and each channel only may occupy the bandwidth of several K, satellite is also the band in this way, satellite beams
It is wide very wide, but apply channel that can not all take bandwidth, at this moment just will appear extreme case is exactly a satellite beams band
It is wide by tens, sometime may only have several K datas to transmit above, and other bandwidth are idle state entirely, the present invention is implemented
The method that example proposes also can fast and accurately trace into wave beam in the state that satellite bandwidth is almost comprehensively idle.
Specifically, the ADC sampling rate of the system of reception meets fs=2 × fd.The base band that reception system samples ADC
The points that signal does FFT processing are NF point, and NF generally chooses 2 power side here, general often selection 8192 or 16384.At this point,
It defines each this subchannel of minimum basis to be made of NB point after FFT transform, meets NB=Bw/ (fs/NF), ensure that NB is most here
Amount is close to integer, to reduce the error of subsequent calculating.
Specifically, the frequency that the system of reception obtains the arrival bearing is f_refresh (Hz), namely every T_
Arrival bearing estimation of starting in refresh=1/f_refresh seconds, this parameter can sets itself according to actual needs, generally
Per second refresh 50~150 times of car terminals setting in mountain area is sufficient for communication performance.
Specifically, the useful signal reference value is obtained, wherein the useful signal reference value is εmin~εmax, εminFor
The lower limit of the useful signal reference value, εmaxFor the upper limit of the useful signal reference value, the two are fixed values, can be artificial
It need to be set according to practical reception system.
Specifically, the bay number of phased-array antenna array is M, the bay number M is generally chosen
Tens.
Specifically, according to actual needs, the division number for receiving system to 360 degree of attitude is Nazi, to space
90 degree of pitch angle of division number is Npit.Here attitude and space pitch angle can be divided equidistantly, it can also not
It is equidistant to divide, sets itself can be needed according to reception system.Preferably, the value range of Nazi is 30~60 degree, Npit's
Value range is 10~20 degree.In addition, according to different regions service condition difference, such as region of the equator uses, basic satellite is just
In 90 degree of positions of pitch angle, polar regions, satellite, can be according to actual use region, weight just in the position that pitch angle is almost 0 degree
Point divides common position.The dividing condition of the scanning space meets Nspa=Nazi × Npit.
Specifically, array manifold and received wave that arrival bearing's vector acquiring unit 11 is formed according to M array element
The all possible azimuth of beam and pitch angle, precalculate or microwave dark room measurement in advance obtains Nspa M rank two-dimensional surface battle array
Arrival bearing's vector of column, is denoted as Al(i), wherein I=1 ..., Nspa；I=1 ..., M
It is worth noting that selecting different array manifolds, the direction vector of two-dimensional surface battle array is different, such as plane
The direction vector of circle battle array meets:
Wherein i=1 ... M, θlFor the pitch angle of the i-th space block divided,For the i-th space block that has divided
Azimuth, R are the radius of circular array, and λ is wavelength, and j is pure imaginary number, if becoming square matrix, direction vector can change again, are needed
Recalculate formula.But it has to be noticed that above formula is all theoretical calculation formula, due to actual hardware material and technique
Difference, cause array antenna anti-line vector can with theoretical value gap, so preferably being gone on the basis of theoretical calculation micro-
Wave darkroom carries out correction in kind, to obtain optimal direction vector.
Specifically, detecting timer first, when the pulse that timer was generated every T_refresh seconds arrives, preservation is worked as
The data of preceding moment M road ADC acquisition are NF point per the points for saving base band data all the way.Since ADC works always, institute
With ADC can output data always but due to the limitation of antenna beam tracking refreshing frequency be not all to use all the time
The data of ADC acquisition calculate, and only take the data that ADC samples out in fixed time, so most of moment, ADC acquisition
Data be abandon it is unused.Preferably, general ground-plane antenna selects the number for taking an ADC to sample out every 10ms to 50ms
According to proper.
Then, the fft processing unit 12 does the FFT operation of NF point to the base band sample data on the road M respectively, it is assumed that the road M
Base band sample data is xm(n), then the base band sample data after FFT operation meets:
Xm(K)=FFT (xm(n)) formula (2)；
Wherein, XmIt (K) is the base band sample data after FFT operation；M=1 ..., M；N=1 ..., NF；K=1 ..., NF.
It is worth noting that the spectral range that the result after above-mentioned FFT operation indicates is [- fs/2, fs/2], the general (scene FPGA
Programmable gate array) result that the IP kernel of FFT operation calculates is done in the inside is not this range, so should be noted that when realizing pair
Operation result is corrected.
Finally, the fft processing unit 12 is to the transformed road M Xm(K) operation is taken to obtain Xm(K '), every road only take
NF/2 point, the predetermined number point is NF/2, and just corresponding frequency spectrum is exactly the signal for needing to scan to these points on every road
Frequency spectrum specifically takes a formula as follows:
Xm(K ')=Xm(K '+NF/4) formula (3)；
Wherein, m=1 ..., M；K '=1 ..., NF/2.Since the echo signal bandwidth to be detected only has fd (Hz), and
The expression range of signal is [- fd, fd] after FFT operation, thus takes Xm(K) NF/2 point come obtain echo signal frequency spectrum [-
Fd/2, fd/2].
Specifically, the frequency domain Beam synthesis computing unit 13 is by Xm(K ') and arrival bearing's vector Al(i) it carries out
Frequency domain Beam synthesis calculates, and specific formula for calculation is as follows:
Wherein, Y (p, l) is Beam synthesis amount；P=1 ..., Nbas；I=1 ..., Nspa.Frequency domain wave beam in formula (4)
Synthesis is exactly with the X receivedmArrival bearing's vector A of (K ') signal and divided all spacesl(i) traversal synthesis one
Time, namely exhaustion retrieval.
Specifically, the Beam synthesis amount maximum value acquiring unit 14 searches the Beam synthesis amount Y (p, l) progress signal
Rope adds two dimensional beam direction to adjudicate, and is found first in Nbas Y (p, l) for belonging to same this subchannel of minimum basis here
Maximum value, maximum value Ymax(p) it indicates, the maximum value Y of the Beam synthesis amountmax(p) corresponding I (of scanning space
Number) use lmaxIt indicates.The Beam synthesis amount maximum value acquiring unit 14 is to find most to have in all possible arrival bearing
The value of possible arrival bearing, the Beam synthesis amount Y (p, l) are bigger, then show that this segment data signal is most likely to be from I
What a direction came.
Specifically, the maximum value Y of the signal detection and normalized unit 15 first to the Beam synthesis amountmax
(p) signal detection is done, the process of the signal detection includes: will be in the same p and Ymax(p) all institutes of 1dB are theoretically differed
It states Beam synthesis amount Y (p, l) and adds up and be denoted as Ymax1dB(p), Ymax1dBIt (p) is signal energy detection limit, all Y in the same p
The gross energy of (p, l) is denoted as Ysum(p).Since FFT operation can be such that signal energy reveals, thus the point of maximum value can not represent
Whole energy of signal select to be construed as the energy portion of signal leakage with the signal of maximum value difference 1dB size here
Divide, thus it is mutually added up to the detection limit as signal energy, and the signal of entire frequency spectrum is added and is used as signal detection amount
Reduced value, seeking the two variables is to do normalized to signal energy detection limit for next, normalized
Benefit is exactly that detection threshold is not influenced by the variation of ADC sampled signal amplitude.
Then, the signal detection and normalized unit 15 are to Ymax1dB(p) normalized is done, specific formula is such as
Under:
Wherein, p=1 ..., Nbas；R (p) is to signal energy detection limit Ymax1dB(p) normalized done as a result,
That is parameter value.
Specifically, the useful signal judging unit 16 judges whether the parameter value R (p) is effective more than or equal to described
The lower limit ε of signal reference valuemin, while judging the maximum value Y of the Beam synthesis amountmax(p) have described in whether being less than or equal to
Imitate the upper limit ε of signal reference valuemax。εminAnd εmaxFor judging whether deposit in this subchannel bandwidth of p-th of current time minimum basis
In useful signal.
If R (p) >=εminAnd Ymax(p)≤εmax, then the useful signal judging unit 16 determines p-th of minimum basis book
There are useful signals in channel width.
If R (p) < εmin, then the useful signal judging unit 16 determines this sub-channel of p-th of current time minimum basis
Without useful signal in wide；If Ymax(p) > εmax, then the useful signal judging unit 16 determines p-th of minimum basis of current time
There is interference signal in this subchannel bandwidth, joined the judgement of interference processing part herein, so that method of the invention is in reality
There is certain anti-interference ability in.As R (p) < εminOr Ymax(p) > εmaxWhen, upper a period of time of output storage storage
Calculated azimuth and pitch angle are carved, FFT processing is done to the data of the road M ADC acquisition again, waits the incoming wave of subsequent time
Direction estimation, thus achieve the effect that target frequency bands broadband beams in real time, quickly, accurately track.It is relatively high for mobility
The estimation frequency can be improved to achieve the effect that almost uninterrupted real-time estimation beam direction in terminal.
Specifically, when the useful signal judging unit 16 determines that the final result obtains single there are when useful signal
Member 17 finds the maximum value R in all parameter value R (p) for being determined to have useful signalmax(p), Rmax(p)=
maxP=1 ..., Nbas{R(p)}.The maximum value R of the parameter valuemax(p) corresponding lmaxFor the final arrival bearing.Above-mentioned mistake
Journey is that the strongest sub-channels of signal energy are found from all this subchannels of minimum basis, with the direction of the sub-channel signal
Final beam direction as output.
Specifically, the final result acquiring unit 17 is according to lmaxIt corresponding can find the final arrival bearing's
Final bearing and final pitch angle, by this moment calculated final bearing and the final pitch angle deposit storage
In device, and export the final bearing and the final pitch angle.Each scanning space I can correspond to an azimuth and
Pitch angle can find its corresponding azimuth and pitch angle by tabling look-up from I.Certainly will also successfully estimate every time here
Azimuth and pitch angle moment out updates storage, cannot be accurate to prevent sometime because of the factors such as signal interference
Use when estimating azimuth and pitch angle.
When it is implemented, first by frequency domain Beam synthesis computing unit 13 to finish FFT processing each circuit-switched data with
Arrival bearing's vector carries out the calculating of frequency domain Beam synthesis, obtains Beam synthesis amount；Then signal detection and normalized
Unit 15 does signal detection and normalized to the maximum value of the Beam synthesis amount, obtains parameter value；When the parameter value
More than or equal to preset useful signal reference value, and the maximum value of the Beam synthesis amount is less than or equal to the useful signal
When reference value, useful signal judging unit 16 determines that there are useful signals in current this subchannel of minimum basis；Last final result
Acquiring unit 17 obtains the maximum value of the parameter value, and obtains the final bearing of wave beam according to the maximum value of the parameter value
With final pitch angle.
Compared with prior art, it is existing to solve the prior art for the digital servomechanism 10 of broadband beams disclosed by the invention
Some mechanical scanning antennas cause it that cannot reach the ICBM SHF satellite terminal of high-speed mobile since its automatic reaction time is slower
The problem of real-time, quickly and accurately tracing into signal.The arrival bearing of current time satellite-signal can be quickly estimated, with
Achieve the effect that phased array antenna it is real-time to broadband beams, it is quick, accurately track.
Embodiment three
Referring to fig. 4, Fig. 4 is a kind of structural frames of the figure tracking equipment 20 of broadband beams provided in an embodiment of the present invention
Figure.The figure tracking equipment 20 of broadband beams described in the embodiment includes: processor 21, memory 22 and is stored in institute
State the computer program that can be run in memory 22 and on the processor 21.The processor 21 executes the computer journey
The step in above-mentioned each screen control method embodiment, such as step S11~S17 shown in FIG. 1 are realized when sequence.Alternatively, institute
State the function that each unit in above-mentioned each Installation practice is realized when processor 21 executes the computer program, such as the incoming wave
The function of direction vector acquiring unit 11.
Illustratively, the computer program can be divided into one or more module/units, one or more
A module/unit is stored in the memory, and is executed by the processor 21, to complete the present invention.It is one or
Multiple module/units can be the series of computation machine program instruction section that can complete specific function, and the instruction segment is for describing
Implementation procedure of the computer program in the figure tracking equipment 20 of the broadband beams.For example, the computer program
Arrival bearing's vector acquiring unit 11, fft processing unit 12, frequency domain Beam synthesis computing unit 13, wave beam can be divided into
Synthetic quantity maximum value acquiring unit 14, signal detection and normalized unit 15, useful signal judging unit 16, final result
Acquiring unit 17 and parameter acquiring unit 18, each module concrete function refer to the figure tracking of above-described embodiment middle width strip wave beam
The function of each unit in system 10, details are not described herein.
The figure tracking equipment 20 of the broadband beams can be desktop PC, notebook, palm PC and cloud
Server etc. calculates equipment.The figure tracking equipment 20 of the broadband beams may include, but be not limited only to, processor 21, storage
Device 22.It will be understood by those skilled in the art that the schematic diagram is only the example of the figure tracking equipment 20 of broadband beams, and
The restriction to the figure tracking equipment 20 of broadband beams is not constituted, may include than illustrating more or fewer components or group
Certain components or different components are closed, such as the figure tracking equipment 20 of the broadband beams can also include input and output
Equipment, network access equipment, bus etc..
The processor 21 can be central processing unit (Central Processing Unit, CPU), can also be
Other general processors, digital signal processor (Digital Signal Processor, DSP), specific integrated circuit
(Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor is also possible to any conventional processor
Deng the processor 21 is the control centre of the figure tracking equipment 20 of the broadband beams, and various interfaces and route is utilized to connect
Connect the various pieces of the figure tracking equipment 20 of entire broadband beams.
The memory 22 can be used for storing the computer program and/or module, the processor 22 by operation or
Computer program and/or the module stored in the memory is executed, and calls the data being stored in memory, is realized
The various functions of the figure tracking equipment 20 of the broadband beams.The memory 22 can mainly include storing program area and deposit
Store up data field, wherein application program needed for storing program area can store operating device, at least one function (for example sound is broadcast
Playing function, image player function etc.) etc.；Storage data area, which can be stored, uses created data (such as audio according to mobile phone
Data, phone directory etc.) etc..In addition, the memory 22 may include high-speed random access memory, it can also include non-volatile
Property memory, such as hard disk, memory, plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital
(Secure Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or other
Volatile solid-state part.
Wherein, if the integrated module/unit of the figure tracking equipment 20 of the broadband beams is with SFU software functional unit
Form realize and when sold or used as an independent product, can store in a computer readable storage medium.Base
In such understanding, the present invention realizes all or part of the process in above-described embodiment method, can also pass through computer program
It is completed to instruct relevant hardware, the computer program can be stored in a computer readable storage medium, the calculating
Machine program is when being executed by processor, it can be achieved that the step of above-mentioned each embodiment of the method.Wherein, the computer program includes
Computer program code, the computer program code can for source code form, object identification code form, executable file or certain
A little intermediate forms etc..The computer-readable medium may include: any entity that can carry the computer program code
Or device, recording medium, USB flash disk, mobile hard disk, magnetic disk, CD, computer storage, read-only memory (ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software
Distribution medium etc..It should be noted that the content that the computer-readable medium includes can be according to making laws in jurisdiction
Requirement with patent practice carries out increase and decrease appropriate, such as in certain jurisdictions, according to legislation and patent practice, computer
Readable medium does not include electric carrier signal and telecommunication signal.
It should be noted that the apparatus embodiments described above are merely exemplary, wherein described be used as separation unit
The unit of explanation may or may not be physically separated, and component shown as a unit can be or can also be with
It is not physical unit, it can it is in one place, or may be distributed over multiple network units.It can be according to actual
It needs that some or all of the modules therein is selected to achieve the purpose of the solution of this embodiment.In addition, device provided by the invention
In embodiment attached drawing, the connection relationship between module indicate between them have communication connection, specifically can be implemented as one or
A plurality of communication bus or signal wire.Those of ordinary skill in the art are without creative efforts, it can understand
And implement.
The above is a preferred embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (8)
1. a kind of digital tracking method of broadband beams characterized by comprising
Obtain arrival bearing's vector of two-dimensional planar array；
The data for obtaining ADC acquisition, do FFT processing to each circuit-switched data respectively；
The calculating of frequency domain Beam synthesis is carried out to each circuit-switched data and arrival bearing's vector of finishing FFT processing, obtains wave beam
Synthetic quantity；
The maximum value of the Beam synthesis amount described in same this subchannel of minimum basis of acquisition；
Signal detection and normalized are done to the maximum value of the Beam synthesis amount, obtain parameter value；
When the parameter value be greater than or equal to preset useful signal reference value, and the maximum value of the Beam synthesis amount be less than or
When equal to the useful signal reference value, determine that there are useful signals in current this subchannel of minimum basis；
When there are when useful signal, obtain the maximum value of the parameter value, the corresponding scanning of the maximum value of the parameter value for judgement
Space is final arrival bearing, and obtains the final bearing and final pitch angle of wave beam according to the final arrival bearing；
It is wherein, described to do FFT processing respectively to per data all the way, comprising:
FFT operation is done respectively to each circuit-switched data；
An operation is taken to finish that each circuit-switched data after the FFT operation takes predetermined number point.
2. the digital tracking method of broadband beams as described in claim 1, which is characterized in that the acquisition two-dimensional planar array
Arrival bearing's vector, comprising:
According to the azimuth and pitch angle of the array manifold of bay number composition and reception wave beam, two-dimensional surface battle array is obtained
Arrival bearing's vector of column.
3. the digital tracking method of broadband beams as described in claim 1, which is characterized in that the acquisition is in the same minimum
The maximum value of Beam synthesis amount described in basic subchannel, comprising:
Signal search and the judgement of two dimensional beam direction are carried out to the Beam synthesis amount, obtained in same this subchannel of minimum basis
Described in Beam synthesis amount maximum value.
4. the digital tracking method of broadband beams as described in claim 1, which is characterized in that the acquisition two-dimensional planar array
Arrival bearing's vector before further include:
Obtain conventional parameter, frequency range including scanning, number, the ADC sampling rate, FFT processing of minimum basis this subchannel
Points, the frequency of arrival bearing and useful signal reference value；
The parameter for obtaining phase battle array control aerial array, including bay number, the dividing condition of attitude, the space angle of depression
The dividing condition of dividing condition and scanning space.
5. a kind of digital servomechanism of broadband beams characterized by comprising
Arrival bearing's vector acquiring unit, for obtaining arrival bearing's vector of two-dimensional planar array；
Fft processing unit does FFT processing to each circuit-switched data for obtaining the data of ADC acquisition respectively；
Frequency domain Beam synthesis computing unit, for being carried out to each circuit-switched data and arrival bearing's vector of finishing FFT processing
Frequency domain Beam synthesis calculates, and obtains Beam synthesis amount；
Beam synthesis amount maximum value acquiring unit, for obtaining the Beam synthesis amount described in same this subchannel of minimum basis
Maximum value；
Signal detection and normalized unit are done at signal detection and normalization for the maximum value to the Beam synthesis amount
Reason, obtains parameter value；
Useful signal judging unit, for being greater than or equal to preset useful signal reference value, and the wave when the parameter value
The maximum value of beam synthetic quantity be less than or equal to the useful signal reference value when, determine in current this subchannel of minimum basis there are
Imitate signal；
Final result acquiring unit, for when there are when useful signal, obtain the maximum value of the parameter value, the parameter for judgement
The corresponding scanning space of the maximum value of value is final arrival bearing, and obtains the final of wave beam according to the final arrival bearing
Azimuth and final pitch angle；
Wherein, the fft processing unit is for doing FFT operation to each circuit-switched data respectively, and to the FFT operation is finished after
What each circuit-switched data took predetermined number point takes operation.
6. the digital servomechanism of broadband beams as claimed in claim 5, which is characterized in that arrival bearing's vector obtains
Azimuth and the pitch angle of array manifold and reception wave beam that unit is used to be formed according to bay number, it is flat to obtain two dimension
Arrival bearing's vector of face array.
7. the digital servomechanism of broadband beams as claimed in claim 5, which is characterized in that the Beam synthesis amount maximum value
Acquiring unit is used to carry out signal search to the Beam synthesis amount and two dimensional beam direction is adjudicated, and obtains in the same minimum basis
The maximum value of Beam synthesis amount described in this subchannel；
The system also includes:
Parameter acquiring unit, for obtaining the parameter of conventional parameter and phase battle array control aerial array；Wherein, the conventional parameter includes
The frequency range of scanning, number, ADC sampling rate, the points, the frequency of arrival bearing of FFT processing of minimum basis this subchannel
With useful signal reference value；The parameter of the phase battle array control aerial array includes bay number, the division feelings of attitude
The dividing condition of condition, the dividing condition at the space angle of depression and scanning space.
8. a kind of figure tracking equipment of broadband beams, including processor, memory and storage in the memory and by
It is configured to the computer program executed by the processor, the processor realizes aforesaid right when executing the computer program
It is required that the digital tracking method of broadband beams described in any one of 1~4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810724910.2A CN109120305B (en) | 2018-07-04 | 2018-07-04 | A kind of digital tracking method of broadband beams, system and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810724910.2A CN109120305B (en) | 2018-07-04 | 2018-07-04 | A kind of digital tracking method of broadband beams, system and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109120305A CN109120305A (en) | 2019-01-01 |
| CN109120305B true CN109120305B (en) | 2019-12-03 |
Family
ID=64821980
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810724910.2A Active CN109120305B (en) | 2018-07-04 | 2018-07-04 | A kind of digital tracking method of broadband beams, system and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109120305B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111798869B (en) * | 2020-09-10 | 2020-11-17 | 成都启英泰伦科技有限公司 | Sound source positioning method based on double microphone arrays |
| CN116582877A (en) * | 2022-01-29 | 2023-08-11 | 维沃移动通信有限公司 | Incoming wave direction estimation method, terminal and network side equipment |
| CN114745301B (en) * | 2022-03-29 | 2023-06-30 | 中国电子科技集团公司第三十八研究所 | Method and system for detecting single-frequency signal with high dynamic low signal-to-noise ratio |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6816112B1 (en) * | 2003-05-30 | 2004-11-09 | Lockheed Martin Corporation | Hybrid RF/optical acquisition and tracking system and method |
| CN101995566A (en) * | 2010-10-15 | 2011-03-30 | 西安电子科技大学 | System and method for forming digital wave beams of two-dimensional digital array radar |
| CN102522632A (en) * | 2011-12-19 | 2012-06-27 | 东南大学 | Digital beam tracking method for molecular matrix of satellite mobile communication phased-array antenna |
| CN104101786A (en) * | 2014-06-24 | 2014-10-15 | 中国电子科技集团公司第十研究所 | All airspace active multi beam spherical phased array antenna direction diagram measurement system |
| CN104391276A (en) * | 2014-10-08 | 2015-03-04 | 西安电子工程研究所 | Transmit-receive split planar array phased-array radar antenna array and beam former |
| CN107785663A (en) * | 2017-10-10 | 2018-03-09 | 深圳市华讯方舟空间信息产业科技有限公司 | antenna beam attitude control method and system |
| CN107885233A (en) * | 2017-10-19 | 2018-04-06 | 安徽工程大学 | A kind of automatic tracking system of satellite-signal |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103441788B (en) * | 2013-07-04 | 2016-02-24 | 中国电子科技集团公司第十研究所 | A kind of feedback deterministic Larger Dynamic adaptive array antenna beam synthesizing method |
| CN104536017B (en) * | 2015-01-06 | 2016-09-21 | 中国人民解放军国防科学技术大学 | A kind of navigation neceiver STAP method of Beam synthesis after first subspace projection |
| CN105186102B (en) * | 2015-09-15 | 2018-01-19 | 西安星通通信科技有限公司 | Antenna for satellite communication in motion system and tracking based on digital beam tracking |
| CN106712827B (en) * | 2016-11-17 | 2020-07-28 | 上海卫星工程研究所 | Dynamic beam tracking test device and method for satellite-borne digital multi-beam receiving antenna |
| CN106888044A (en) * | 2017-03-28 | 2017-06-23 | 中国电子科技集团公司第三十八研究所 | A kind of optimum synthesis method of round symmetrical antenna Oriented Graphics with Assigned Form |
-
2018
- 2018-07-04 CN CN201810724910.2A patent/CN109120305B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6816112B1 (en) * | 2003-05-30 | 2004-11-09 | Lockheed Martin Corporation | Hybrid RF/optical acquisition and tracking system and method |
| CN101995566A (en) * | 2010-10-15 | 2011-03-30 | 西安电子科技大学 | System and method for forming digital wave beams of two-dimensional digital array radar |
| CN102522632A (en) * | 2011-12-19 | 2012-06-27 | 东南大学 | Digital beam tracking method for molecular matrix of satellite mobile communication phased-array antenna |
| CN104101786A (en) * | 2014-06-24 | 2014-10-15 | 中国电子科技集团公司第十研究所 | All airspace active multi beam spherical phased array antenna direction diagram measurement system |
| CN104391276A (en) * | 2014-10-08 | 2015-03-04 | 西安电子工程研究所 | Transmit-receive split planar array phased-array radar antenna array and beam former |
| CN107785663A (en) * | 2017-10-10 | 2018-03-09 | 深圳市华讯方舟空间信息产业科技有限公司 | antenna beam attitude control method and system |
| CN107885233A (en) * | 2017-10-19 | 2018-04-06 | 安徽工程大学 | A kind of automatic tracking system of satellite-signal |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109120305A (en) | 2019-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109120305B (en) | A kind of digital tracking method of broadband beams, system and equipment | |
| CN111050276B (en) | Method and device for optimizing IRS (intelligent resilient service) auxiliary communication system based on angle information | |
| CN106054123B (en) | A kind of sparse L battle arrays and its arrival direction estimation method | |
| CN110839204B (en) | Communication optimization method and device of IRS auxiliary communication system | |
| CN102175989B (en) | Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival) | |
| CN108646213B (en) | Direct wave AOA (automatic optical inspection) judgment method in indoor multipath environment | |
| CN108445486A (en) | It is rebuild and the modified Beamforming Method of steering vector based on covariance matrix | |
| CN110113088B (en) | Intelligent estimation method for wave arrival angle of separated digital-analog hybrid antenna system | |
| CN109164418B (en) | LPI-based networking radar steady OFDM waveform design method under spectrum sharing environment | |
| CN107124216A (en) | A kind of Capon robust adaptive beamforming method and system for array error | |
| CN111856402B (en) | Signal processing method and device, storage medium and electronic device | |
| CN109143275A (en) | A kind of miniaturization array antenna anti-interference realization method based on population | |
| CN110290542B (en) | Communication coverage optimization method and system for offshore unmanned aerial vehicle | |
| CN110083175B (en) | Unmanned aerial vehicle formation network cooperative scheduling method and device | |
| CN104852758B (en) | Vertical beam shaping method and device under three-dimensional extensive aerial network | |
| CN105612779A (en) | Transmit beamforming sounding with traveling pilots | |
| CN108233984B (en) | Interference suppression method and device | |
| CN108900232A (en) | A kind of method, apparatus and electronic equipment of adaptive beamforming | |
| CN113917949B (en) | Ship-borne three-axis antenna satellite automatic tracking method, system, equipment and storage medium | |
| CN113890798B (en) | Structured sparse estimation method and device for RIS cascade channel multi-user combination | |
| CN106199550B (en) | A kind of wave digression and two-dimentional direction of arrival automatic matching combined estimation method | |
| CN109561381A (en) | Two-frequency signal localization method, device, equipment and the storage medium of mobile terminal | |
| CN112422458B (en) | Channel estimation method, apparatus and computer storage medium | |
| WO2017152504A1 (en) | Beam search method and apparatus | |
| CN1589054A (en) | Up self adaption wave beam shorming method based on signal arrival angle detection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |