CN108303602A - A kind of test method of supersparsity aerial array transmitting-receiving beam pattern - Google Patents
A kind of test method of supersparsity aerial array transmitting-receiving beam pattern Download PDFInfo
- Publication number
- CN108303602A CN108303602A CN201711468663.6A CN201711468663A CN108303602A CN 108303602 A CN108303602 A CN 108303602A CN 201711468663 A CN201711468663 A CN 201711468663A CN 108303602 A CN108303602 A CN 108303602A
- Authority
- CN
- China
- Prior art keywords
- antenna
- delay
- signal
- phase difference
- difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/10—Radiation diagrams of antennas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention discloses the test method that a kind of supersparsity aerial array receives and dispatches beam pattern, mainly solve the problems, such as that existing measurement method and system can not carry out ground Pattern measurement to supersparsity aerial array.The present invention receives the echo-signal of same satellite by multiple antennas, and each antenna delay diference and phase difference, the zero as antenna delay diference and phase difference are obtained using FX correlation synthesizers；Rasterizing is carried out to the arrival bearing in single antenna beam angle on this basis, the satellite-signal that each antenna receives accurately is compensated using corresponding delay difference and phase difference, to realize the test of supersparsity aerial array transmit-receive position figure.The present invention has many advantages, such as that delay/phase measurement accuracy is high, method is simple and reliable, is relatively specific for the far-field pattern test of the supersparsities such as area antennas group battle array, networking radar, distributed aerial array.
Description
Technical field
Antenna array signal synthetic technology, Distributed Network Radar in being communicated the present invention relates to deep space TT＆C are mutually referred and synthesized
At technology, interferometry data treatment technology etc., especially suitable for extensive reception antenna group battle array, networking radar distributed
The related applications fields such as the test of transmitting-receiving Beam synthesis, combined coefficient analysis and the assessment of array antenna.
Background technology
In deep space TT＆C communication system, remote away due to deep-space spacecraft target and earth survey station, from several ten million to
More than one hundred million kilometers cause the received signal strength of existing telemetry communication equipment very faint, and distributed dual-mode antenna group battle array is to realize
It improves transmitting ERIP and receives the effective means of signal SNR；In networking radar system, the target of the array antenna of supersparsity
Detection, AF panel, coherent performance etc. are required to antenna radiation pattern poised for battle and carry out test verification.Traditional phase array antenna beam side
It cannot be satisfied the directional diagram far field test condition of supersparsity aerial array to the ground test method of figure, it is therefore desirable to by in-orbit
The approximation far field objects such as satellite are tested；In addition, different from traditional phased array width phase weighted antenna, sparse array antenna is also
It needs to solve array antenna delay diference, the problems such as phase difference and antenna gain are inconsistent, it is therefore desirable to a kind of new supersparsity
Aerial array receives and dispatches the test method of beam pattern.
Invention content
The invention discloses the test method that a kind of supersparsity aerial array receives and dispatches beam pattern, this method passes through to big
Scale array antenna received satellite-signal carries out the synthesis of FX correlations, high-precision array delay difference/phase difference estimation and compensation, array
Antenna SNR methods of estimation etc. realize the test of supersparsity aerial array transmit-receive position figure.This method sets battle array reference antenna first,
It is right using orbit prediction and reference antenna constructing antennas delay diference and offset model based on aerial position and the thick rail of satellite
The thick alignment of each aerial signal delay difference and phase difference is realized after each antenna reception signal progress coarse delay and phase difference compensation, then
Complex base band signal is obtained after quadrature frequency conversion and integral cleaning and carries out FFT operations, carries out composing certainly in integration period
The cross-spectrum calculating for calculating and carrying out each aerial signal with " total reference signal " in addition to own signal using FX correlations synthesizer, profit
The aerial signal is obtained relative to the residual retardance difference and phase information with total reference signal, for repairing with cross-spectrum phase result
Positive frequency difference model between antenna delay diference and antenna guides the signal of the lower integration period of each antenna to realize delay difference and phase difference
Essence alignment；The signal power of each aerial signal is obtained by iterative calculation using mutual spectral amplitude result, is believed in conjunction with from spectral amplitude
The noise power of each antenna can be calculated in breath, to realize estimating for each antenna snr of received signal and composite signal signal-to-noise ratio
Meter, so far obtains each aerial signal delay diference, phase difference and the estimation with interior signal-to-noise ratio, by all aerial signals by essence
It is true to prolong delay difference, frequency difference and phase difference compensation and be directly added to obtain total composite signal after amplitude weighting；Utilize each antenna SNR
Estimated result and composite signal SNR estimated results obtain array gain value.The result of the delay difference, phase difference and amplitude weighting is made
In system null fill-in to delay difference and offset model, to obtain antenna according to rasterizing difference arrival bearing on this basis
Delay diference and phase difference are added in compensation model, and the battle array to obtain all arrival bearings in single antenna wave beam synthesizes gain
Figure in conjunction with single antenna directional diagram and then obtains the transmitting-receiving beam pattern test result of aerial array.The present invention is relatively specific for
The far-field pattern test of the supersparsities such as area antennas group battle array, networking radar, distributed aerial array, overcomes such array
It can not carry out the problem of ground Pattern measurement.
The object of the present invention is achieved like this：
A kind of test method of supersparsity aerial array transmitting-receiving beam pattern, it is characterised in that include the following steps：
1. antenna delay diference/offset model unit is forecast according to externally input 1~antenna of antenna N delays/phase
Model, selects antenna 1 as antenna is referred to, 1~antenna of antenna N delays/phase forecasting model is postponed with reference antenna respectively/
Phase forecasting model is made the difference to obtain the corresponding antenna delay diference/phase difference multinomial models of 1~antenna of antenna N；
It is respectively right to 1~antenna of antenna N according to the residual retardance difference and residual phase difference initial value from spectrum and the output of cross-spectrum computing unit
The antenna delay diference answered/phase difference multinomial model one-to-one correspondence is modified, and obtains postponing between N number of revised antenna
Difference/phase difference multinomial model, as the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N；Root
According to externally input single antenna beam angle rasterizing information, the two of K*K are carried out to the arrival bearing in single antenna beam angle
Rasterizing is tieed up, delay difference/phase difference multinomial model of K*K different arrival bearings is generated；By K*K delay difference/phase difference
Multinomial model is separately summed to obtain day with the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N
The corresponding antenna delay diference/phase difference compensation models of 1~antenna of line N；By the corresponding antennas of 1~antenna of antenna N
Delay diference compensation model one-to-one correspondence is sent to the corresponding coarse delay compensating units of 1~antenna of antenna N and smart delay compensation
Unit send the corresponding antenna delay diferences of 1~antenna of antenna N/phase difference compensation model one-to-one correspondence to 1~day of antenna
The corresponding carrier phase difference compensating units of line N；Wherein, N is the natural number more than 1, and K is the natural number more than 1；
2. coarse delay compensating unit according to the time scale information of antenna delay diference compensation model and its respective antenna signal into
Row calculates, and obtains delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration, and the delay difference is pre-
Report value is decomposed into integer bit delay compensation value and decimal bit delays offset；Coarse delay compensating unit to aerial signal when
Integer bit delay compensation in the corresponding time of integration is completed in domain, and the signal after compensation is sent to carrier phase difference compensating unit；
3. carrier phase difference compensating unit is according to antenna delay diference/phase difference compensation model and its respective antenna signal
Time scale information is calculated, and delay difference predicted value and phase of the intermediate time spacecraft relative to antenna in the corresponding time of integration are obtained
The delay difference predicted value is separated into integer bit delay compensation value and decimal bit delays offset by potential difference predicted value；Carrier wave
Phase difference compensation unit completes phase difference compensation to the signal that coarse delay compensating unit exports in time domain and by integer bit
Carrier phase difference caused by delay compensation compensates, and the signal after compensation is sent to quadrature frequency conversion and integral cleaning unit；
4. quadrature frequency conversion and integral cleaning unit construct down-conversion signal according to externally input IF carrier nominal value,
Quadrature frequency conversion is carried out to the signal that carrier phase difference compensating unit exports using the down-conversion signal of construction, obtains complex baseband
Signal will integrate I, Q after cleaning after carrying out integral cleaning reduction data rate respectively to I, Q branch of complex base band signal
Baseband signal is sent into FFT unit；
5. FFT unit constructs complex base band signal and simultaneously carries out plural FFT fortune according to I, Q baseband after integral cleaning
It calculates, and FFT result of calculations is sent into smart delay compensation unit；
6. smart delay compensation unit according to the time scale information of antenna delay diference compensation model and its respective antenna signal into
Row calculates, and obtains delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration, and the delay difference is pre-
Report value is decomposed into integer bit delay compensation value and decimal bit delays offset；According to decimal bit delays value and FFT frequencies
The corresponding phase compensation weighted vector of decimal bit delays is constructed, it is small to carry out striped rotation realization to the result of FFT unit output
The essence compensation of bit delays is counted, the result after essence is compensated is respectively fed to amplitude weighting unit and certainly spectrum and cross-spectrum computing unit；
7. corresponding amplitude weighting value is calculated according to externally input antenna radiation pattern yield value in amplitude weighting unit,
Then amplitude weighting is carried out to the result of smart delay compensation unit output, and result after amplitude weighting is sent to synthesis unit；
8. the signal that the corresponding amplitude weighting units of 1~antenna of antenna N export directly is added by synthesis unit, obtain
On the one hand total composite signal utilizes the corresponding baseband signal of each antenna and noise power, the noise of total composite signal is calculated
Power and signal power, and then the interior SNR of band that total composite signal is calculated, and by the interior SNR value of the band of total composite signal send to
Array gain computing unit；On the other hand it is corresponding from spectrum and cross-spectrum calculating total composite signal to be fed back to 1~antenna of antenna N
Unit；
9. the sum total of the signal and synthesis unit output exported from spectrum and cross-spectrum computing unit according to smart delay compensation unit
At signal calculate from spectrum and cross-spectrum, wherein it is each aerial signal and total composite signal respectively to participate in the signal that cross-spectrum calculates
Subtract the reference signal that own antenna signal obtains；By the cross-spectrum phase information of cross-spectrum result by least square linear fit into
Row calculates, and residual retardance difference and residual phase difference is obtained, by residual retardance difference and residual phase difference antennas delay diference/phase
Poor model unit is updated for residual retardance difference and residual phase difference iteration；It is iteratively solved using the amplitude spectrum of cross-spectrum result
To the signal amplitude of each antenna, corresponding noise power is calculated in conjunction with from spectral amplitude information, and then each antenna is calculated
The interior SNR of band of signal, and the interior SNR value of the band of each aerial signal is sent to array gain computing unit；
10. array gain computing unit is according in the band of SNR value and total composite signal in the corresponding bands of 1~antenna of antenna N
The array gain value on the arrival bearing is calculated in SNR value, which is worth corresponding arrival bearing exports together；
For the reception data in the same amount of time of record, generated not according to single antenna beam angle rasterizing
K*K arrival bearing and its corresponding array gain value are obtained 2. to 10. repeating K*K-1 times from step with arrival bearing, it will
Arrival bearing and its corresponding array gain value carry out three-dimensional output of drawing to get to the transmitting-receiving beam gain of sparse array antenna
Directional diagram；
Complete the test of supersparsity aerial array transmitting-receiving beam pattern.
It has the following advantages that compared with the background technology, the present invention：
1. the present invention solves the problems, such as that current sparse array antenna radiation pattern can not carry out ground test；
2. the present invention can realize delay difference, phase difference and the amplitude between each antenna reception signal using FX correlations synthesizer
The accurate estimation of weighting and automatic compensation realize high-precision delay, the three-dimensional weighting of phase and amplitude, it is phased to overcome tradition
Battle array only can phase and amplitude method of weighting sparse array antenna radiation pattern discomfort.
3. supersparsity antenna array pattern can be rapidly completed using one section of satellite-signal for receiving storage in the present invention
Test and generation, the testing time is substantially reduced while improving measuring accuracy.
Description of the drawings
Fig. 1 is the principle of the present invention flow chart；
Specific implementation mode
Referring to Fig.1, the invention will be further described.
A kind of test method of supersparsity aerial array transmitting-receiving beam pattern, it is characterised in that include the following steps：
1. antenna delay diference/offset model unit is forecast according to externally input 1~antenna of antenna N delays/phase
Model, selects antenna 1 as antenna is referred to, 1~antenna of antenna N delays/phase forecasting model is postponed with reference antenna respectively/
Phase forecasting model is made the difference to obtain the corresponding antenna delay diference/phase difference multinomial models of 1~antenna of antenna N；
It is respectively right to 1~antenna of antenna N according to the residual retardance difference and residual phase difference initial value from spectrum and the output of cross-spectrum computing unit
The antenna delay diference answered/phase difference multinomial model one-to-one correspondence is modified, and obtains postponing between N number of revised antenna
Difference/phase difference multinomial model, as the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N；Root
According to externally input single antenna beam angle rasterizing information, the two of K*K are carried out to the arrival bearing in single antenna beam angle
Rasterizing is tieed up, delay difference/phase difference multinomial model of K*K different arrival bearings is generated；By K*K delay difference/phase difference
Multinomial model is separately summed to obtain day with the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N
The corresponding antenna delay diference/phase difference compensation models of 1~antenna of line N；By the corresponding antennas of 1~antenna of antenna N
Delay diference compensation model one-to-one correspondence is sent to the corresponding coarse delay compensating units of 1~antenna of antenna N and smart delay compensation
Unit send the corresponding antenna delay diferences of 1~antenna of antenna N/phase difference compensation model one-to-one correspondence to 1~day of antenna
The corresponding carrier phase difference compensating units of line N；Wherein, N is the natural number more than 1, and K is the natural number more than 1；
Delay/phase forecasting model refers to delay/phase forecast data based on antenna site and the thick rail forecast of satellite, number
It is contained in including caused by the factors such as geometric position, satellite clock correction, troposphere/ionospheric refraction and receiver clock-offsets
Delay/phase information.
Buffer path difference/phase difference multinomial model refers to passing through least square linear according to delay/phase forecast data
It is fitted obtained unitary quartic polynomial f (t), as follows：
τ (t)=a4t4+a3t3+a2t2+a1t+a0
The multinomial model is the function of time t, according to the time t of setting, can be calculated t moment delay difference and
Phase difference.
The delay difference of arrival bearing/phase difference multinomial model refers to carrying out two-dimensional grid according to single antenna beam angle
Metaplasia at arrival bearing and the delay difference that is calculated of aerial array location information and phase difference.
2. coarse delay compensating unit according to the time scale information of antenna delay diference compensation model and its respective antenna signal into
Row calculates, and obtains delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration, and the delay difference is pre-
Report value is decomposed into integer bit delay compensation value and decimal bit delays offset；Coarse delay compensating unit to aerial signal when
Integer bit delay compensation in the corresponding time of integration is completed in domain, and the signal after compensation is sent to carrier phase difference compensating unit；
3. carrier phase difference compensating unit is according to antenna delay diference/phase difference compensation model and its respective antenna signal
Time scale information is calculated, and delay difference predicted value and phase of the intermediate time spacecraft relative to antenna in the corresponding time of integration are obtained
The delay difference predicted value is separated into integer bit delay compensation value and decimal bit delays offset by potential difference predicted value；Carrier wave
Phase difference compensation unit completes phase difference compensation to the signal that coarse delay compensating unit exports in time domain and by integer bit
Carrier phase difference caused by delay compensation compensates, and the signal after compensation is sent to quadrature frequency conversion and integral cleaning unit；
4. quadrature frequency conversion and integral cleaning unit construct down-conversion signal according to externally input IF carrier nominal value,
Quadrature frequency conversion is carried out to the signal that carrier phase difference compensating unit exports using the down-conversion signal of construction, obtains complex baseband
Signal will integrate I, Q after cleaning after carrying out integral cleaning reduction data rate respectively to I, Q branch of complex base band signal
Baseband signal is sent into FFT unit；
5. FFT unit constructs complex base band signal and simultaneously carries out plural FFT fortune according to I, Q baseband after integral cleaning
It calculates, and FFT result of calculations is sent into smart delay compensation unit；
6. smart delay compensation unit according to the time scale information of antenna delay diference compensation model and its respective antenna signal into
Row calculates, and obtains delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration, and the delay difference is pre-
Report value is decomposed into integer bit delay compensation value and decimal bit delays offset；According to decimal bit delays value and FFT frequencies
The corresponding phase compensation weighted vector of decimal bit delays is constructed, it is small to carry out striped rotation realization to the result of FFT unit output
The essence compensation of bit delays is counted, the result after essence is compensated is respectively fed to amplitude weighting unit and certainly spectrum and cross-spectrum computing unit；
7. corresponding amplitude weighting value is calculated according to externally input antenna radiation pattern yield value in amplitude weighting unit,
Then amplitude weighting is carried out to the result of smart delay compensation unit output, and result after amplitude weighting is sent to synthesis unit；
8. the signal that the corresponding amplitude weighting units of 1~antenna of antenna N export directly is added by synthesis unit, obtain
On the one hand total composite signal utilizes the corresponding baseband signal of each antenna and noise power, the noise of total composite signal is calculated
Power and signal power, and then the interior SNR of band that total composite signal is calculated, and by the interior SNR value of the band of total composite signal send to
Array gain computing unit；On the other hand it is corresponding from spectrum and cross-spectrum calculating total composite signal to be fed back to 1~antenna of antenna N
Unit；
9. the sum total of the signal and synthesis unit output exported from spectrum and cross-spectrum computing unit according to smart delay compensation unit
At signal calculate from spectrum and cross-spectrum, wherein it is each aerial signal and total composite signal respectively to participate in the signal that cross-spectrum calculates
Subtract the reference signal that own antenna signal obtains；By the cross-spectrum phase information of cross-spectrum result by least square linear fit into
Row calculates, and residual retardance difference and residual phase difference is obtained, by residual retardance difference and residual phase difference antennas delay diference/phase
Poor model unit is updated for residual retardance difference and residual phase difference iteration；It is iteratively solved using the amplitude spectrum of cross-spectrum result
To the signal amplitude of each antenna, corresponding noise power is calculated in conjunction with from spectral amplitude information, and then each antenna is calculated
The interior SNR of band of signal, and the interior SNR value of the band of each aerial signal is sent to array gain computing unit；
10. array gain computing unit is according in the band of SNR value and total composite signal in the corresponding bands of 1~antenna of antenna N
The array gain value on the arrival bearing is calculated in SNR value, which is worth corresponding arrival bearing exports together；
For the reception data in the same amount of time of record, generated not according to single antenna beam angle rasterizing
K*K arrival bearing and its corresponding array gain value are obtained 2. to 10. repeating K*K-1 times from step with arrival bearing, it will
Arrival bearing and its corresponding array gain value carry out three-dimensional output of drawing to get to the transmitting-receiving beam gain of sparse array antenna
Directional diagram；
Complete the test of supersparsity aerial array transmitting-receiving beam pattern.
The brief operation principle of the present invention is as follows：
A kind of test method of supersparsity aerial array transmitting-receiving beam pattern is relatively specific for area antennas group battle array, network
Change the far-field pattern test of the supersparsities such as radar, distributed aerial array, overcoming such array can not be in ground progress side
To the problem of figure test.
The invention discloses the test method that a kind of supersparsity aerial array receives and dispatches beam pattern, this method passes through to big
Scale array antenna received satellite-signal carries out the synthesis of FX correlations, high-precision array delay difference/phase difference estimation and compensation, array
Antenna SNR methods of estimation etc. realize the test of supersparsity aerial array transmit-receive position figure.This method sets battle array reference antenna first,
It is right using orbit prediction and reference antenna constructing antennas delay diference and offset model based on aerial position and the thick rail of satellite
The thick alignment of each aerial signal delay difference and phase difference is realized after each antenna reception signal progress coarse delay and phase difference compensation, then
Complex base band signal is obtained after quadrature frequency conversion and integral cleaning and carries out FFT operations, carries out composing certainly in integration period
The cross-spectrum calculating for calculating and carrying out each aerial signal with " total reference signal " in addition to own signal using FX correlations synthesizer, profit
The aerial signal is obtained relative to the residual retardance difference and phase information with total reference signal, for repairing with cross-spectrum phase result
Positive frequency difference model between antenna delay diference and antenna guides the signal of the lower integration period of each antenna to realize delay difference and phase difference
Essence alignment；The signal power of each aerial signal is obtained by iterative calculation using mutual spectral amplitude result, is believed in conjunction with from spectral amplitude
The noise power of each antenna can be calculated in breath, to realize estimating for each antenna snr of received signal and composite signal signal-to-noise ratio
Meter, so far obtains each aerial signal delay diference, phase difference and the estimation with interior signal-to-noise ratio, by all aerial signals by essence
It is true to prolong delay difference, frequency difference and phase difference compensation and be directly added to obtain total composite signal after amplitude weighting；Utilize each antenna SNR
Estimated result and composite signal SNR estimated results obtain array gain value.The result of the delay difference, phase difference and amplitude weighting is made
In system null fill-in to delay difference and offset model, to obtain antenna according to rasterizing difference arrival bearing on this basis
Delay diference and phase difference are added in compensation model, and the battle array to obtain all arrival bearings in single antenna wave beam synthesizes gain
Figure in conjunction with single antenna directional diagram and then obtains the transmitting-receiving beam pattern test result of aerial array.
Claims (1)
1. a kind of test method of supersparsity aerial array transmitting-receiving beam pattern, it is characterised in that include the following steps：
1. antenna delay diference/offset model unit according to externally input 1~antenna of antenna N delays/phase forecasting model,
Select antenna 1 as refer to antenna, by 1~antenna of antenna N delays/phase forecasting model respectively with reference antenna delay/phase
Forecasting model is made the difference to obtain the corresponding antenna delay diference/phase difference multinomial models of 1~antenna of antenna N；According to
It is corresponding to 1~antenna of antenna N from the residual retardance difference and residual phase difference initial value of spectrum and the output of cross-spectrum computing unit
Antenna delay diference/phase difference multinomial model one-to-one correspondence is modified, and obtains N number of revised antenna delay diference/phase
Potential difference multinomial model, as the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N；According to outside
The single antenna beam angle rasterizing information of input carries out the arrival bearing in single antenna beam angle the two-dimensional grid of K*K
Change, generates delay difference/phase difference multinomial model of K*K different arrival bearings；By K*K delay difference/phase difference multinomial
Model be separately summed to obtain with the corresponding buffer path difference/phase difference multinomial models of 1~antenna of antenna N antenna 1~
The corresponding antenna delay diference/phase difference compensation models of antenna N；It will postpone between the corresponding antennas of 1~antenna of antenna N
Poor compensation model one-to-one correspondence is sent to the corresponding coarse delay compensating units of 1~antenna of antenna N and smart delay compensation unit, will
The corresponding antenna delay diferences of 1~antenna of antenna N/phase difference compensation model one-to-one correspondence is sent to 1~antenna of antenna N respectively
Corresponding carrier phase difference compensating unit；Wherein, N is the natural number more than 1, and K is the natural number more than 1；
2. coarse delay compensating unit is counted according to the time scale information of antenna delay diference compensation model and its respective antenna signal
It calculates, delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration is obtained, by the delay difference predicted value
It is decomposed into integer bit delay compensation value and decimal bit delays offset；Coarse delay compensating unit is complete in time domain to aerial signal
At integer bit delay compensation in the corresponding time of integration, and the signal after compensation is sent to carrier phase difference compensating unit；
3. carrier phase difference compensating unit is according to the markers of antenna delay diference/phase difference compensation model and its respective antenna signal
Information is calculated, and delay difference predicted value and phase difference of the intermediate time spacecraft relative to antenna in the corresponding time of integration are obtained
The delay difference predicted value is separated into integer bit delay compensation value and decimal bit delays offset by predicted value；Carrier phase
Poor compensating unit completes phase difference compensation to the signal that coarse delay compensating unit exports in time domain and by integer bit delays
Carrier phase difference compensation caused by compensation, and the signal after compensation is sent to quadrature frequency conversion and integral cleaning unit；
4. quadrature frequency conversion and integral cleaning unit construct down-conversion signal according to externally input IF carrier nominal value, utilize
The down-conversion signal of construction carries out quadrature frequency conversion to the signal that carrier phase difference compensating unit exports, and obtains complex baseband letter
Number, after carrying out integral cleaning reduction data rate respectively to I, Q branch of complex base band signal, I, Q base after cleaning will be integrated
Band signal is sent into FFT unit；
5. FFT unit constructs complex base band signal and simultaneously carries out plural FFT operations according to I, Q baseband after integral cleaning, and
FFT result of calculations are sent into smart delay compensation unit；
6. smart delay compensation unit is counted according to the time scale information of antenna delay diference compensation model and its respective antenna signal
It calculates, delay difference predicted value of the intermediate time spacecraft relative to antenna in the corresponding time of integration is obtained, by the delay difference predicted value
It is decomposed into integer bit delay compensation value and decimal bit delays offset；It is constructed according to decimal bit delays value and FFT frequencies
The corresponding phase compensation weighted vector of decimal bit delays carries out striped rotation to the result of FFT unit output and realizes decimal ratio
The essence compensation that spy postpones, the result after essence is compensated are respectively fed to amplitude weighting unit and certainly spectrum and cross-spectrum computing unit；
7. corresponding amplitude weighting value is calculated according to externally input antenna radiation pattern yield value in amplitude weighting unit, then
Amplitude weighting is carried out to the result of smart delay compensation unit output, and result after amplitude weighting is sent to synthesis unit；
8. the signal that the corresponding amplitude weighting units of 1~antenna of antenna N export directly is added by synthesis unit, summed up
At signal, the corresponding baseband signal of each antenna and noise power are on the one hand utilized, the noise power of total composite signal is calculated
And signal power, and then the interior SNR of band that total composite signal is calculated, and the interior SNR value of the band of total composite signal is sent to battle array and is increased
Beneficial computing unit；On the other hand it is corresponding single from spectrum and cross-spectrum calculating total composite signal to be fed back to 1~antenna of antenna N
Member；
9. the sum total of signal and the synthesis unit output exported from spectrum and cross-spectrum computing unit according to smart delay compensation unit is at letter
Number carry out from spectrum and cross-spectrum calculate, wherein participate in cross-spectrum calculate signal be that each aerial signal and total composite signal subtract respectively
The reference signal that own antenna signal obtains；The cross-spectrum phase information of cross-spectrum result is counted by least square linear fit
It calculates, residual retardance difference and residual phase difference is obtained, by residual retardance difference and residual phase difference antennas delay diference/phase differential mode
Type unit is updated for residual retardance difference and residual phase difference iteration；It iteratively solves to obtain using the amplitude spectrum of cross-spectrum result each
The signal amplitude of antenna is calculated corresponding noise power in conjunction with from spectral amplitude information, and then each aerial signal is calculated
The interior SNR of band, and the interior SNR value of the band of each aerial signal is sent to array gain computing unit；
10. array gain computing unit is according to the interior SNR of the band of SNR value and total composite signal in the corresponding bands of 1~antenna of antenna N
Value, is calculated the array gain value on the arrival bearing, which is worth corresponding arrival bearing exports together；
For the reception data in the same amount of time of record, the different incoming waves generated according to single antenna beam angle rasterizing
Direction, 2. to 10. repeating K*K-1 times, obtains K*K arrival bearing and its corresponding array gain value, by incoming wave side from step
To and its corresponding array gain value carry out three-dimensional output of drawing to get to the transmitting-receiving beam gain direction of sparse array antenna
Figure；
Complete the test of supersparsity aerial array transmitting-receiving beam pattern.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711468663.6A CN108303602B (en) | 2017-12-29 | 2017-12-29 | Test method for receiving and transmitting beam pattern of ultra-sparse antenna array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711468663.6A CN108303602B (en) | 2017-12-29 | 2017-12-29 | Test method for receiving and transmitting beam pattern of ultra-sparse antenna array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108303602A true CN108303602A (en) | 2018-07-20 |
| CN108303602B CN108303602B (en) | 2020-05-08 |
Family
ID=62868029
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711468663.6A Active CN108303602B (en) | 2017-12-29 | 2017-12-29 | Test method for receiving and transmitting beam pattern of ultra-sparse antenna array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108303602B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018364A (en) * | 2019-05-07 | 2019-07-16 | 中国人民解放军32039部队 | Antenna radiation pattern Orbital detection method, system and electronic equipment |
| CN110109048A (en) * | 2019-05-23 | 2019-08-09 | 北京航空航天大学 | A kind of invasion signal arrival bearing's angular range estimation method based on phase difference |
| CN110708105A (en) * | 2019-08-15 | 2020-01-17 | 北京航天驭星科技有限公司 | Method for developing satellite measurement, operation and control by using ground station antenna array |
| CN110927751A (en) * | 2019-12-09 | 2020-03-27 | 中国电子科技集团公司第五十四研究所 | Array antenna self-adaptive correction implementation method based on carrier phase measurement |
| CN113203897A (en) * | 2021-04-14 | 2021-08-03 | 之江实验室 | Side lobe suppression and beam forming method of arbitrary two-dimensional antenna array |
| CN113381786A (en) * | 2021-06-11 | 2021-09-10 | 军事科学院系统工程研究院网络信息研究所 | Multi-satellite co-location distributed constellation antenna array cooperative transmission method |
| RU2789466C1 (en) * | 2022-03-29 | 2023-02-03 | Публичное акционерное общество "Радиофизика" | Method for measuring the characteristics of the radiation pattern of a digital phased array antenna |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012052856A1 (en) * | 2010-10-21 | 2012-04-26 | Reutech Radar Systems (Proprietary) Limited | Floodlight radar system for detecting and locating moving targets in three dimensions |
| CN102957472A (en) * | 2012-09-18 | 2013-03-06 | 西安电子科技大学 | Method and device for receiving communication signal of thinned array antenna system |
-
2017
- 2017-12-29 CN CN201711468663.6A patent/CN108303602B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012052856A1 (en) * | 2010-10-21 | 2012-04-26 | Reutech Radar Systems (Proprietary) Limited | Floodlight radar system for detecting and locating moving targets in three dimensions |
| CN102957472A (en) * | 2012-09-18 | 2013-03-06 | 西安电子科技大学 | Method and device for receiving communication signal of thinned array antenna system |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018364A (en) * | 2019-05-07 | 2019-07-16 | 中国人民解放军32039部队 | Antenna radiation pattern Orbital detection method, system and electronic equipment |
| CN110018364B (en) * | 2019-05-07 | 2020-07-31 | 中国人民解放军32039部队 | Antenna directional pattern on-orbit testing method and system and electronic equipment |
| CN110109048A (en) * | 2019-05-23 | 2019-08-09 | 北京航空航天大学 | A kind of invasion signal arrival bearing's angular range estimation method based on phase difference |
| CN110708105A (en) * | 2019-08-15 | 2020-01-17 | 北京航天驭星科技有限公司 | Method for developing satellite measurement, operation and control by using ground station antenna array |
| CN110927751A (en) * | 2019-12-09 | 2020-03-27 | 中国电子科技集团公司第五十四研究所 | Array antenna self-adaptive correction implementation method based on carrier phase measurement |
| CN113203897A (en) * | 2021-04-14 | 2021-08-03 | 之江实验室 | Side lobe suppression and beam forming method of arbitrary two-dimensional antenna array |
| CN113381786A (en) * | 2021-06-11 | 2021-09-10 | 军事科学院系统工程研究院网络信息研究所 | Multi-satellite co-location distributed constellation antenna array cooperative transmission method |
| RU2789466C1 (en) * | 2022-03-29 | 2023-02-03 | Публичное акционерное общество "Радиофизика" | Method for measuring the characteristics of the radiation pattern of a digital phased array antenna |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108303602B (en) | 2020-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108303602A (en) | A kind of test method of supersparsity aerial array transmitting-receiving beam pattern | |
| CN108365872B (en) | A kind of signal synthesis method suitable for strange land multiple antennas group battle array | |
| CN101980048B (en) | Antenna array forming technology-based ground-based radar system for space debris | |
| RU2565386C2 (en) | Method, apparatus and system for determining position of object, having global navigation satellite system receiver, by processing non-differential data, similar to carrier phase measurements, and external data similar to ionospheric data | |
| US9304184B1 (en) | System and method for mitigating severe multipath interference for geolocation and navigation | |
| US8184042B2 (en) | Self calibrating conformal phased array | |
| US7855681B2 (en) | Systems and methods for determining element phase center locations for an array of antenna elements | |
| CN104137447B (en) | Estimate the method and system of the echo signal progressive error of spacecraft or vehicle launch | |
| US9099776B2 (en) | Method for iterative estimation of global parameters | |
| CN104714244A (en) | Multi-system dynamic PPP resolving method based on robust self-adaption Kalman smoothing | |
| EP3106897A1 (en) | Gnss receiver with an on-board capability to implement an optimal error correction mode | |
| JP5759676B2 (en) | Propagation path estimation system and propagation path estimation method | |
| Shah et al. | Precision of Ku-band reflected signals of opportunity altimetry | |
| CN107132423B (en) | A kind of method and device of ionization detection layer electron density sum | |
| Sweeney et al. | Enabling Mars radio occultation by smallsats | |
| Han et al. | Angle of arrival and time of flight estimation as an ising energy minimization problem | |
| Joseph et al. | Development of VHF (240–270 MHz) antennas for SoOp (Signal of Opportunity) receiver for 6U cubesat platforms | |
| CN113671535B (en) | Tri-frequency beacon receiver separation layer TEC calculation method based on channel simulator | |
| Bhardwaj et al. | Investigation of Ionospheric Vertical Delay at S1 and L5 Frequencies, Based on Thick-Shell Model Using NavIC System, for Mid Latitude Region of India. | |
| Kim et al. | Development of an end-to-end mission simulator for land remote sensing with signals of opportunity | |
| Lanyi et al. | Determination of angular separation between spacecraft and quasars with the very long baseline array | |
| Weiss | Getting accurate time from GNSS receivers: Considerations to approach nanosecond time | |
| RU2309425C2 (en) | Method of forming calibration data for radio direction finder/ range finder (versions) | |
| Li et al. | Extraction of electron density profiles with geostationary satellite-based GPS side lobe occultation signals | |
| Rodon et al. | Results on spatial-temporal atmospheric phase screen retrieval from long-term GEOSAR acquisition |
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 |