CN103117781A - Method and device for antenna array calibration under complex electromagnetic environment - Google Patents

Method and device for antenna array calibration under complex electromagnetic environment Download PDF

Info

Publication number
CN103117781A
CN103117781A CN2013100410054A CN201310041005A CN103117781A CN 103117781 A CN103117781 A CN 103117781A CN 2013100410054 A CN2013100410054 A CN 2013100410054A CN 201310041005 A CN201310041005 A CN 201310041005A CN 103117781 A CN103117781 A CN 103117781A
Authority
CN
China
Prior art keywords
signal
channel
array
frequency
passage
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
Application number
CN2013100410054A
Other languages
Chinese (zh)
Other versions
CN103117781B (en
Inventor
马洪
饶蓉
李攀
余德
程莉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201310041005.4A priority Critical patent/CN103117781B/en
Publication of CN103117781A publication Critical patent/CN103117781A/en
Application granted granted Critical
Publication of CN103117781B publication Critical patent/CN103117781B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention discloses a method and a device for antenna array calibration under complex electromagnetic environment. The device transmits calibration signals from a center point of an antenna array, and extracts array antennas and receives front-end amplitude inconsistency and phase inconsistency by means of wiener solution according to the principle of adaptive filter. The influence of small-amplitude same frequency interference or adjacent channel interference upon the result of frequency array calibration can be inhibited automatically in the complex electromagnetic environment with intense interference. When high-amplitude same frequency interference or adjacent channel interference has great influence on channel parameters, the extracted channel parameters can be processed to eliminate or greatly weaken the influence of the interference; received array signals are then subjected to broadband compensation or multiple narrowband compensation, and accordingly errors caused by channel inconsistency are eliminated. The method and device are applicable to real-time channel parameter extraction and automatic compensation for circular antenna arrays and receiving channels operating in any frequency range under the complex electromagnetic environment.

Description

Antenna array calibration method and device thereof under a kind of complex electromagnetic environment
Technical field
The present invention relates to the Array Signal Processing field, be specially antenna array calibration method and device thereof under a kind of complex electromagnetic environment.Complex electromagnetic environment refers to that there is the environment of intensive interference signal in the external world.Circular array refers to that antenna is evenly distributed on circumference.The present invention extracts and compensates circular array channel amplitude and phase place inconsistency error under complex electromagnetic environment.
Background technology
Develop rapidly along with communication theory and technology, array antenna has been played the part of more and more important role in the radio systems such as mobile communication, broadcasting station, radar system and navigation, and the quality of Array Signal Processing performance will directly affect the performance of whole communication system.Greatly reduce the performance of array acceptor system due to the inconsistency of each passage of array antenna, these parameters mainly comprise amplitude and the phase place of each reception antenna and receiver front end equipment, therefore, must adopt various technology and means to solve as much as possible the inconsistency of array antenna passage, each passage be carried out the calibration of amplitude inconsistency and phase place inconsistency.
At present, most of array calibration document is to the calibration in narrow band signal, or the calibration of carrying out under very little and interference signal amplitude circlet border in the external interference closeness, realizes that background is fairly simple." Sensor array calibration in the presence of mutual coupling and unknown sensor gains and phases " (C.M.S See, Electron Lett, 1994, 30 (5): 373-374) propose to send the calibration auxiliary signal, ask the covariance matrix of the array signal of reception, and associating Estimation of Spatial Spectrum, try to achieve channel amplitude inconsistency and phase place inconsistency, but this applied environment is restricted, external interference will be less and is faint, and in the situation that intensive external interference and interference signal are stronger, the method is just inapplicable.
Generally, the assisted calibration signal of transmission is all easy identification, and the simple signal that can not change of the waveform on frequency domain.Exactly because the intrinsic attribute of simple signal, the external interference pulse will affect the reception of this auxiliary signal to be differentiated, and affects the result of array calibration but simultaneously.Characteristics of the present invention are to solve existing array calibration method not consider external interference and can only be in the problem of arrowband internal calibration, Acquisition channel parameter (comprising channel amplitude inconsistency parameter and channel phases inconsistency parameter) accurately, and compensate in broadband.
Summary of the invention
The invention provides antenna array calibration method and device thereof under a kind of complex electromagnetic environment, the present invention can and exist under the complex electromagnetic environment of intensive interference in any frequency band range, automatically suppresses by a small margin same frequency or the impact of adjacent interference frequently on this frequency array calibration result.
Antenna array calibration method under a kind of complex electromagnetic environment provided by the invention comprises the following steps:
The 1st step sent calibrating signal;
The 2nd step utilized circular array to receive the calibrating signal that sends, and each road signal that receives is carried out analog-to-digital conversion, obtained the array digital signal;
The 3rd step was utilized the array digital signal that obtains, extract amplitude and the phase place of each passage, again with the amplitude of a certain passage and phase place as benchmark, ask the amplitude of other each passage and the relative value of phase place and benchmark, i.e. channel amplitude inconsistency and channel phases inconsistency;
The 4th step judge whether all frequencies are calibrated completely, if so, enter for the 5th step, otherwise, switch to next frequency, repeated for the 1st step to go on foot to the 3rd, until complete the extraction of the channel parameters of all frequencies;
The 5th step stopped sending calibrating signal;
The 6th step processed channel parameters, eliminated the impact of disturbing;
On the 7th step circular array, each reception antenna receives aerial signal, and each road signal that will receive carries out analog-to-digital conversion, obtains the array digital signal;
The 8th step compensated the array digital signal that obtains, and was compensated output signal.
Antennas array calibration arrangement under a kind of complex electromagnetic environment provided by the invention, it is characterized in that, this device comprises signal generator unit, circular antenna array element, array acceptor front end unit, array channel parameter extraction unit, channel parameters processing unit and channel compensation unit;
Described circular antenna array element is used for receiving signal and transmitting calibration signal, a calibrating signal transmitting antenna that is used for sending calibrating signal is placed at the center of circle battle array, also evenly place many reception antennas on the circumference of described circular array, distance between every adjacent reception antenna all equates, this reception antenna is extracting the array channel parameter stage, receives the calibrating signal that is mixed with external interference; In the compensated array signal phase, receive required array signal to be calibrated, described signal generator unit for generation of simple signal as the calibrating signal in the array calibration process;
Described array acceptor front end unit be used for to the signal that receives amplify, filtering, analog-to-digital conversion process, and offer described array channel parameter extraction unit;
The array signal that described array channel parameter extraction unit be used for to receive and the Wei Na solution form of reference signal according to sef-adapting filter carry out computing, extract channel parameters, and offer described channel parameters processing unit;
Process if described channel parameters processing unit is used for calibration frequency or adjacent place are frequently existed in strongly disturbing situation, and the data after processing offer the channel compensation unit; Somely disturb littlely, signal interference ratio is larger, does not process, and obtains final channel parameters;
Described channel compensation unit comprises broadband compensation unit and/or many narrow-band compensations unit, eliminates the array channel error that the signal that receives comprises.
The present invention can and exist under the complex electromagnetic environment of intensive interference in any frequency band range, automatically suppress by a small margin same frequency or the adjacent impact of frequently disturbing this frequency array calibration result, when if the Amplitude Ratio that disturbs is larger, can this frequency channel parameters of automatic rejection, and interpolation replaces it, obtain accurately final array channel parameter, and carry out many narrow-band compensations, perhaps disposable compensation in broad frequency range.Compare with other array calibration technology, the present invention has following characteristics:
(1) the present invention according to the principle of sef-adapting filter, uses Wei Na solution form to carry out the extraction of channel parameters, and operand is little, calculates simply, and easily realizes.
(2) can namely exist in the situation of intensive interference at complex electromagnetic environment, according to the interference magnitude size at calibration frequency or adjacent frequency place, channel parameters is processed.Avoided conventional array calibration method easily to be subjected to the impact of external interference.
(3) the array channel compensation can be carried out many narrow-band compensations or broadband compensation, and broadband compensation is the structure calibration filters, realizes the disposable compensation in required frequency band; Many narrow-band compensations are the signal channellizations that will receive, and are divided into a plurality of arrowbands, then compensate in each arrowband.Having avoided the compensation of conventional array channel to require the signal of compensation to be calibrated is signal in the arrowband.
(4) apparatus of the present invention can be based upon on the basis of any band limits, have avoided the general array calibration can only be in the arrowband scope, and disturb in the situation of wanting smaller just feasible in this frequency range.
Description of drawings
Fig. 1 is the method flow schematic diagram of example of the present invention;
Fig. 2 is channel parameters extracting method block diagram;
Fig. 3 is channel parameters processing method block diagram;
Fig. 4 is broadband compensation method block diagram;
Fig. 5 is many narrow-band compensations method block diagram;
Fig. 6 is the apparatus structure schematic diagram of example of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described further.Need to prove at this, understand the present invention for the explanation of these execution modes for helping, but do not consist of limitation of the invention.In addition, below in each execution mode of described the present invention involved technical characterictic just can mutually not make up as long as consist of each other conflict.
As shown in Figure 1, the array calibration method that example of the present invention provides comprises the following steps:
(1) send calibrating signal
Signal generator, is sent by the calibrating signal transmitting antenna through after digital-to-analogue conversion according to the parameter generating simple signal of the calibrating signal of waiting to send out.The parameter of calibrating signal comprises type, amplitude and the frequency etc. of signal.
(2) data acquisition
After sending calibrating signal, utilize circular array to receive each road signal, obtain digital signal after analog-to-digital conversion, the passage ginseng that is used for the back is extracted calculating.
(3) carry out the channel parameters extraction algorithm
Each railway digital signal that utilization obtains carries out computing, extract amplitude and the phase place of each passage, again with the amplitude of a certain passage and phase place as benchmark, ask the amplitude of other each passage and the relative value of phase place and benchmark, i.e. channel amplitude inconsistency and channel phases inconsistency.As shown in Figure 3, concrete steps are as follows:
(3.1) each reception antenna of circular array receives the calibrating signal that sends from the calibrating signal transmitting antenna, and the interference signals such as air broadcast radio station, air communications, and it is s that the signal that i reception antenna receives passes through the analog-to-digital conversion postscript i(n), n is the sequence number of sampled point.
(3.2) produce single-frequency sequence cos2 π fn with the identical frequency of calibrating signal, as reference signal 1, and sin2 π fn, as reference signal 2, the calibration signal frequency of f for launching.
(3.3) according to the signal of reference signal and reception, carry out related operation, extract channel parameters, specific as follows described:
Supposing to be positioned at the calibrating signal that the calibrating signal transmitting antenna of circular array center launches is cos2 π ft, and t is the time of emission.I reception antenna receives this calibrating signal and after amplification, filtering and analog-to-digital conversion, its output signal is
Figure BDA00002809571000051
A iFor the amplitude of this passage output signal,
Figure BDA00002809571000052
Be the phase shift of this passage output signal, they have reflected the amplitude of this receive path introducing and the inconsistency of phase place indirectly, and its value can based on the lowest mean square criterion, with the Wei Na solution form of sef-adapting filter, namely be calculated in the steps below:
Calculate the ASSOCIATE STATISTICS amount w of this passage opt:
w opt = I i Q i = 2 N Σ n = 1 N { s i ( n ) cos 2 πfn } - Σ n = 1 N { s i ( n ) sin 2 πfn
In following formula, f is the frequency of the calibrating signal of launching, and N is for being used for calculating ASSOCIATE STATISTICS amount w optTotal number of sample points, obtain weights I iAnd Q i, and utilize following formula to calculate the amplitude A of i passage i, phase place
Figure BDA00002809571000061
A i = I i 2 + Q i 2 ,
Figure BDA00002809571000063
(3.4) according to A iWith
Figure BDA00002809571000064
Ask for the amplitude of each passage and the relative value of phase place and benchmark, i.e. the amplitude inconsistency of passage and phase place inconsistency are as channel parameters.Can be with amplitude and the phase parameter A of first passage 1,
Figure BDA00002809571000065
As benchmark, ask for the amplitude of each passage and the deviation of phase place and this benchmark, obtain channel amplitude inconsistency parameter A i'=A i/ A 1With phase place inconsistency parameter
Figure BDA00002809571000066
Wherein, i is channel position, and 1≤i≤P, P are the array received total number of channels.
(4) the judgement frequency has been calibrated no
The present invention is applicable in any frequency range, if need to ask the channel parameters in the broad frequency band scope, need send the calibrating signal of a plurality of frequencies.After sending the calibrating signal of first frequency and extracting the channel parameters of this frequency, the judgement frequency has been calibrated no, do not calibrated when frequency and switched to next frequency, repeated (1)~(3) step, until complete the extraction of the channel parameters of all frequencies.
(5) stop sending calibrating signal
After the channel parameters extraction was complete, each reception antenna did not need to receive calibrating signal, therefore, stop sending calibrating signal, and the shutdown signal generator makes the calibrating signal transmitting antenna no longer send calibrating signal.
(6) treatment channel parameter
The present invention processes channel parameters, eliminates the impact of disturbing, and processing procedure as shown in Figure 4.The concrete processing procedure of this step is:
(6.1) suppose that a plurality of calibration frequencies that arrange are f in the working frequency range of required calibration k, k=1,2 ..., K, the calibration frequency total number of K for arranging, the channel amplitude relative value of extracting at each calibration frequency is A ' I, k, and the phase place relative value is
Figure BDA00002809571000067
I is channel position.If at this calibration frequency place or the adjacent frequency place external interference is arranged, but interference magnitude is smaller, the channel parameters at this frequency place directly is taken as A ' I, kWith
Figure BDA00002809571000071
If the external interference Amplitude Ratio is larger, need to process.
Example of the present invention is according to resulting phase place-frequency curve The burr situation determine whether to process the channel parameters at this frequency place.Wherein, A ' I, kWith
Figure BDA00002809571000073
The amplitude, the phase place inconsistency parameter that represent i passage, k frequency place.
At phase place-frequency curve
Figure BDA00002809571000074
In, judge and process the method for burr point: a threshold value T is set, when
Figure BDA00002809571000075
Figure BDA00002809571000076
And Judge this frequency place Be burr, at this moment with the channel amplitude relative value A ' at this frequency place I, kWith the phase place relative value
Figure BDA00002809571000079
Deletion.T is by adjacent calibration frequency interval, passage length, and signaling rate determines, because
Figure BDA000028095710000710
The unit of T is degree, and Δ f is the difference on the frequency of two calibrating signals, and c is the transmission speed of signal, and l is the length of passage.
(6.2) at the deleted enterprising row interpolation of frequency of channel parameters, replace deleted channel parameters with interpolation result, as the channel parameters at this frequency place.
Concrete interpolation method can be the methods such as three rank polynomial interopolations, with the phase place-frequency curve after the assurance interpolation
Figure BDA000028095710000711
Level and smooth in whole working frequency range.Can also double according to demand or many times of interpolation on the basis of calibrated frequency points, namely increase the total number K of calibration frequency, obtain channel parameters corresponding to more intensive frequency.
(6.3) repeating step (6.1)~(6.2) obtain the channel parameters after all passages are processed.
(7) data acquisition
On circular array, each reception antenna receives required aerial signal, each road signal of synchronous acquisition.
(8) channel compensation
Channel compensation has broadband compensation method and many narrow-band compensations method.
The broadband compensation process the invention provides the method that the signal that each passage is received carries out amplitude and phase compensation in broadband as shown in Figure 5, and concrete steps are as follows:
(a1), the structural wideband compensating filter, according to the symmetry of FIR filter parameter, the channel parameters A ' that each calibration frequency place of i paths is obtained I, kWith
Figure BDA00002809571000081
Frequency domain response H as the broadband compensation filter i(f) parameter on positive frequency band is with A ' J, kWith
Figure BDA00002809571000082
As H i(f) parameter on negative frequency band, then to H i(f) carry out Fourier inversion, can obtain the time domain impulse response h of broadband compensation filter i(n).
(a2), the signal x to be calibrated that the i passage is received i(n) enter in digital compensation filter, with the time domain impulse response h of broadband compensation filter i(n) carry out convolution, that is:
S i ′ = h i ( n ) ⊗ x i ( n )
s′ I (n)For through the signal after the i channel compensation.
(a3), repeating step (a1)~(a2), obtain the reception signal after all channel compensations.
But the present invention also signal that receives of each passage carries out many narrow-band compensations, process as shown in Figure 6, concrete steps are as follows:
(b1), the wideband digital signal x that the i passage is received i(n), process through channelizing, namely dedicated bandwidth is that Δ f, centre frequency are respectively f mSeveral band pass filters to the reception signal x of Whole frequency band i(n) carry out filtering, obtain the channelizing output signal:
y i , m ( n ) = x i ( n ) ⊗ g i , m ( n )
g I, m(n) be the impulse response of m channel band-pass filter of i passage, m=1,2 ..., M, M are the channel number.
After processing by channelizing, with wideband received signal x i(n) become a plurality of arrowband complex signal y I, m(n) (channelizing signal).
(b2), according to channel parameters, structure i passage, the interior compensating factor r of a m channel I, m:
Figure BDA00002809571000085
J is empty unit, and e is the truth of a matter of natural logrithm, A ' I, mBe the amplitude inconsistency parameter in i passage, a m channel,
Figure BDA00002809571000086
It is the phase place inconsistency parameter in i passage, a m channel.
(b3), the channel parameters in i passage, a m channel is compensated in the arrowband complex signal of respective channels, obtain the compensation output signal s at this passage, this channel place " I, m(n):
s″ i,m(n)=r i,my i,m(n)
The reception signal that utilizes following formula to obtain this passage is positioned at the compensation output signal at all narrow band channels place.
(b4), repeat (b1)~(b3), try to achieve all passages in the compensation output signal at each narrow band channel place.
The present invention is based on circular antenna array, existing under the complex electromagnetic environment of intensive interference, to the calibration of each passage.A calibrating signal transmitting antenna 10 that is used for sending calibrating signal is placed at the circular array center, and the simple signal of signal generator unit 11 generations is as the calibrating signal in the array calibration process.Evenly place P root reception antenna 12 on circular array, extracting the array channel parameter stage, receive the calibrating signal that is mixed with external interference; In the compensated array signal phase, receive the array signal of required calibration.
Array calibration device of the present invention as shown in Figure 1, array acceptor front end 13 can comprise the modules such as amplification, filtering, frequency conversion, analog-to-digital conversion.14 li of array channel parameter extraction unit, the array signal of reception and reference signal are separated form according to the Wei Na of sef-adapting filter, carry out related operation, extract channel parameters.If channel parameters processing unit 15 is calibration frequency or adjacent place frequently to be existed in strongly disturbing situation process, somely to disturb littlely, signal interference ratio is larger, does not process, and obtains final channel parameters.Channel compensation unit 16 comprises broadband compensation unit 160 and many narrow-band compensations unit 161, i.e. the array channel error that the signal that receives comprises is eliminated in two kinds of compensation methodes.
The above is preferred embodiment of the present invention, but the present invention should not be confined to the disclosed content of this embodiment and accompanying drawing.So everyly do not break away from the equivalence of completing under spirit disclosed in this invention or revise, all falling into the scope of protection of the invention.

Claims (9)

1. the antenna array calibration method under a complex electromagnetic environment comprises the following steps:
The 1st step sent calibrating signal;
The 2nd step utilized circular array to receive the calibrating signal that sends, and each road signal that receives is carried out analog-to-digital conversion, obtained the array digital signal;
The 3rd step was utilized the array digital signal that obtains, extract amplitude and the phase place of each passage, again with the amplitude of a certain passage and phase place as benchmark, ask the amplitude of other each passage and the relative value of phase place and benchmark, i.e. channel amplitude inconsistency and channel phases inconsistency;
The 4th step judge whether all frequencies are calibrated completely, if so, enter for the 5th step, otherwise, switch to next frequency, repeated for the 1st step to go on foot to the 3rd, until complete the extraction of the channel parameters of all frequencies;
The 5th step stopped sending calibrating signal;
The 6th step processed channel parameters, eliminated the impact of disturbing;
On the 7th step circular array, each reception antenna receives aerial signal, and each road signal that will receive carries out analog-to-digital conversion, obtains the array digital signal;
The 8th step compensated the array digital signal that obtains, and was compensated output signal.
2. antenna array calibration method according to claim 1, is characterized in that, the 3rd step specifically comprised following process:
(3.1) each reception antenna of circular array receives the calibrating signal that sends from the calibrating signal transmitting antenna, and interference signal, and the signal of i reception antenna reception is s through the analog-to-digital conversion postscript i(n), n is the sequence number of sampled point;
(3.2) produce single-frequency sequence cos2 π fn with the identical frequency of calibrating signal, as reference signal 1, and sin2 π fn, as reference signal 2, the calibration signal frequency of f for launching;
(3.3) according to the signal of reference signal and reception, carry out related operation, extract channel parameters, comprise amplitude A i, phase place
Figure FDA00002809570900021
(3.4) with the amplitude of one of them passage and phase parameter as benchmark, ask for the amplitude of each passage and the deviation of phase value and this benchmark, obtain channel amplitude inconsistency parameter and phase place inconsistency parameter.
3. antenna array calibration method according to claim 2, it is characterized in that, the 3rd step comprised following process: step (3.3) is: supposing to be positioned at the calibrating signal that the calibrating signal transmitting antenna of circular array center launches is cos2 π ft, t is the time of emission, i reception antenna receives this calibrating signal and after amplification, filtering and analog-to-digital conversion, its output signal is
Figure FDA00002809570900022
A iBe i passage output signal amplitude,
Figure FDA00002809570900023
Be the phase shift of this passage output signal, calculate the ASSOCIATE STATISTICS amount w of this passage opt:
w opt = I i Q i = 2 N Σ n = 1 N { s i ( n ) cos 2 πfn } - Σ n = 1 N { s i ( n ) sin 2 πfn }
In following formula, f is the frequency of the calibrating signal of launching, and N is for being used for calculating ASSOCIATE STATISTICS amount w optTotal number of sample points, obtain weights I iAnd Q i, and utilize
Figure FDA00002809570900025
Figure FDA00002809570900026
Calculate the amplitude A of i passage i, phase place
4. according to claim 1,2 or 3 described antenna array calibration methods, is characterized in that, the 6th step specifically comprised following process:
Suppose that a plurality of calibration frequencies that arrange are f in the working frequency range of required calibration k, k=1,2 ..., K, the calibration frequency total number of K for arranging, the channel amplitude relative value of extracting at each calibration frequency is A ' i,k, and the phase place relative value is I is channel position; If at certain calibration frequency place or adjacent frequency place, external interference is arranged, but interference magnitude is smaller, the channel parameters at this frequency place directly is taken as A ' i,kWith
Figure FDA00002809570900029
Otherwise, to reject the channel parameters of this frequency, then carry out interpolation processing, the channel parameters before replacing is eliminated the impact of disturbing.
5. antenna array calibration method according to claim 4, is characterized in that, the 6th the step according to resulting phase place-frequency curve
Figure FDA00002809570900031
The burr situation determine whether to process the channel parameters at this frequency place.
6. antenna array calibration method according to claim 5, is characterized in that, at phase place-frequency curve
Figure FDA00002809570900032
In, judge and process the method for burr point:
(6.1) when And
Figure FDA00002809570900035
Judge this frequency place
Figure FDA00002809570900036
Be burr, at this moment with the channel amplitude relative value A ' at this frequency place I, kWith the phase place relative value Deletion, wherein, T is predefined threshold value;
(6.2) at the deleted enterprising row interpolation of frequency of channel parameters, replace deleted channel parameters with interpolation result, as the channel parameters at this frequency place;
(6.3) repeating step (6.1) ~ (6.2) obtain the channel parameters after all passages are processed.
7. arbitrary described linear array calibration steps according to claim 1 to 6, is characterized in that, the 8th step was carried out broadband compensation according to following process:
(a1), the structural wideband compensating filter, according to the symmetry of FIR filter parameter, the channel parameters A ' that each calibration frequency place of i paths is obtained i,kWith
Figure FDA00002809570900038
Frequency domain response H as the broadband compensation filter i(f) parameter on positive frequency band is with A ' i,kWith
Figure FDA00002809570900039
As H i(f) parameter on negative frequency band, then to H i(f) carry out Fourier inversion, namely obtain the time domain impulse response h of broadband compensation filter i(n);
(a2), the signal x to be calibrated that the i passage is received i(n) enter in digital compensation filter, with the time domain impulse response h of broadband compensation filter i(n) carry out convolution, that is:
s i ′ ( n ) = h i ( n ) ⊗ x i ( n )
s′ I (n)For through the signal after the i channel compensation;
(a3), repeating step (a1) ~ (a2), obtain the reception signal after all channel compensations.
8. arbitrary described antenna array calibration method in 6 according to claim 1, is characterized in that, the 8th step was carried out many narrow-band compensations according to following process:
(b1), the signal x to be calibrated that the i passage is received i(n), process through channelizing, obtain the channelizing output signal:
y i , m ( n ) = x i ( n ) ⊗ g i , m ( n )
g I, m(n) be the impulse response of m channel band-pass filter of i passage, m=1,2 ..., M, M are the channel number,
After processing by channelizing, with wideband received signal x i(n) become a plurality of arrowband complex signal y i,m(n);
(b2), according to channel parameters, structure i passage, the interior compensating factor r of a m channel i,m:
J is empty unit, and e is the truth of a matter of natural logrithm, A ' I, mBe the amplitude inconsistency parameter in i passage, a m channel,
Figure FDA00002809570900043
It is the phase place inconsistency parameter in i passage, a m channel;
(b3), the channel parameters in i passage, a m channel is compensated in the arrowband complex signal of respective channels, obtain the compensation output signal s at this passage, this channel place " I, m(n):
s″ i,m(n)=r i,my i,m(n)
The reception signal that utilizes following formula to obtain this passage is positioned at the compensation output signal at all narrow band channels place;
(b4), repeat (b1) ~ (b3), try to achieve all passages in the compensation output signal at each narrow band channel place.
9. the antennas array calibration arrangement under a complex electromagnetic environment, it is characterized in that, this device comprises signal generator unit, circular antenna array element, array acceptor front end unit, array channel parameter extraction unit, channel parameters processing unit and channel compensation unit;
Described circular antenna array element is used for receiving signal and transmitting calibration signal, a calibrating signal transmitting antenna that is used for sending calibrating signal is placed at the center of circle battle array, evenly place many reception antennas on the circumference of circular array, distance between every adjacent reception antenna all equates, this reception antenna is used for being mixed with in the array channel parameter stage reception of extraction the calibrating signal of external interference, in the compensated array signal phase, receive required array signal to be calibrated; Described signal generator unit for generation of simple signal as the calibrating signal in the array calibration process;
Described array acceptor front end unit be used for to the signal that receives amplify, filtering, analog-to-digital conversion process, and offer described array channel parameter extraction unit;
Described array channel parameter extraction unit is used for the array signal and the Wei Na solution form of reference signal according to sef-adapting filter that will receive, and carries out computing, extracts channel parameters, and offers described channel parameters processing unit;
Process if described channel parameters processing unit is used for calibration frequency or adjacent place are frequently existed in strongly disturbing situation, and the data after processing offer the channel compensation unit; Somely disturb littlely, signal interference ratio is larger, does not process, and obtains final channel parameters;
Described channel compensation unit comprises broadband compensation unit and/or many narrow-band compensations unit, eliminates the array channel error that the signal that receives comprises.
CN201310041005.4A 2013-02-01 2013-02-01 A kind of antenna array calibration method under complex electromagnetic environment and device thereof Expired - Fee Related CN103117781B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310041005.4A CN103117781B (en) 2013-02-01 2013-02-01 A kind of antenna array calibration method under complex electromagnetic environment and device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310041005.4A CN103117781B (en) 2013-02-01 2013-02-01 A kind of antenna array calibration method under complex electromagnetic environment and device thereof

Publications (2)

Publication Number Publication Date
CN103117781A true CN103117781A (en) 2013-05-22
CN103117781B CN103117781B (en) 2016-04-13

Family

ID=48416064

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310041005.4A Expired - Fee Related CN103117781B (en) 2013-02-01 2013-02-01 A kind of antenna array calibration method under complex electromagnetic environment and device thereof

Country Status (1)

Country Link
CN (1) CN103117781B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104506253A (en) * 2015-01-13 2015-04-08 重庆大学 Amplitude phase error correction system and method for transmitting channel of phased-array antenna
CN106019212A (en) * 2016-05-16 2016-10-12 中国电子科技集团公司第十研究所 Method for correcting non-ideal factor among multiple channels
WO2017201765A1 (en) * 2016-05-27 2017-11-30 Huawei Technologies Co., Ltd. Method and apparatus for beam forming calibration in point to multi point communication systems
CN108549058A (en) * 2018-02-28 2018-09-18 四川九洲电器集团有限责任公司 Secondary radar broadband active phased array system and dynamic calibration method thereof
CN110149157A (en) * 2018-02-11 2019-08-20 西南电子技术研究所(中国电子科技集团公司第十研究所) Array antenna wideband channel parallel calibration method
CN110940955A (en) * 2019-10-24 2020-03-31 湖北三江航天险峰电子信息有限公司 Radar polarization self-adaptive identification method and system
CN111562445A (en) * 2020-05-06 2020-08-21 中国人民解放军63892部队 Real-time monitoring method for angular simulation precision of radio frequency simulation test system
CN111711500A (en) * 2020-05-06 2020-09-25 中国人民解放军63892部队 Simulation antenna array calibration and radio frequency signal monitoring system
CN111999782A (en) * 2020-07-30 2020-11-27 北京遥测技术研究所 Automatic calibration method of rotary scanning imaging system
CN112083389A (en) * 2020-06-09 2020-12-15 成都微联云智科技有限责任公司 Real-time calibration method for circular array antenna
CN112350751A (en) * 2020-11-04 2021-02-09 中国电子科技集团公司第五十四研究所 Multi-channel wide-band amplitude and phase calibration device for satellite communication
CN112346023A (en) * 2019-08-06 2021-02-09 北京华航无线电测量研究所 Amplitude self-calibration module of passive radar
CN113075701A (en) * 2021-03-25 2021-07-06 陕西凌云电器集团有限公司 Satellite signal receiving method and module with differential positioning and anti-interference functions
CN114513228A (en) * 2022-04-19 2022-05-17 中国人民解放军海军工程大学 L-band high-speed frequency hopping data link non-cooperative interference cancellation device and method
CN117169801A (en) * 2023-11-02 2023-12-05 成都德辰博睿科技有限公司 Electromagnetic environment monitoring and calibrating system, method, device and medium

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424730A (en) * 2007-11-02 2009-05-06 杰脉通信技术(上海)有限公司 Multichannel positioning system calibrating method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424730A (en) * 2007-11-02 2009-05-06 杰脉通信技术(上海)有限公司 Multichannel positioning system calibrating method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张帆: "VHF均匀圆阵有源校准方案研究与实现", 《硕士学位论文》, 31 July 2012 (2012-07-31) *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104506253A (en) * 2015-01-13 2015-04-08 重庆大学 Amplitude phase error correction system and method for transmitting channel of phased-array antenna
CN106019212A (en) * 2016-05-16 2016-10-12 中国电子科技集团公司第十研究所 Method for correcting non-ideal factor among multiple channels
CN106019212B (en) * 2016-05-16 2018-08-31 中国电子科技集团公司第十研究所 A kind of bearing calibration for non-ideal factor between multichannel
WO2017201765A1 (en) * 2016-05-27 2017-11-30 Huawei Technologies Co., Ltd. Method and apparatus for beam forming calibration in point to multi point communication systems
US9948407B2 (en) 2016-05-27 2018-04-17 Huawei Technologies Co., Ltd. Method and apparatus for beamforming calibration in point to multipoint communication systems
CN110149157A (en) * 2018-02-11 2019-08-20 西南电子技术研究所(中国电子科技集团公司第十研究所) Array antenna wideband channel parallel calibration method
CN108549058A (en) * 2018-02-28 2018-09-18 四川九洲电器集团有限责任公司 Secondary radar broadband active phased array system and dynamic calibration method thereof
CN112346023A (en) * 2019-08-06 2021-02-09 北京华航无线电测量研究所 Amplitude self-calibration module of passive radar
CN110940955A (en) * 2019-10-24 2020-03-31 湖北三江航天险峰电子信息有限公司 Radar polarization self-adaptive identification method and system
CN111711500B (en) * 2020-05-06 2022-08-19 中国人民解放军63892部队 Simulation antenna array calibration and radio frequency signal monitoring system
CN111562445B (en) * 2020-05-06 2022-08-19 中国人民解放军63892部队 Real-time monitoring method for angular simulation precision of radio frequency simulation test system
CN111562445A (en) * 2020-05-06 2020-08-21 中国人民解放军63892部队 Real-time monitoring method for angular simulation precision of radio frequency simulation test system
CN111711500A (en) * 2020-05-06 2020-09-25 中国人民解放军63892部队 Simulation antenna array calibration and radio frequency signal monitoring system
CN112083389A (en) * 2020-06-09 2020-12-15 成都微联云智科技有限责任公司 Real-time calibration method for circular array antenna
CN112083389B (en) * 2020-06-09 2024-03-15 成都明思数源科技有限公司 Real-time calibration method for circular array antenna
CN111999782A (en) * 2020-07-30 2020-11-27 北京遥测技术研究所 Automatic calibration method of rotary scanning imaging system
CN111999782B (en) * 2020-07-30 2023-02-03 北京遥测技术研究所 Automatic calibration method of rotary scanning imaging system
CN112350751A (en) * 2020-11-04 2021-02-09 中国电子科技集团公司第五十四研究所 Multi-channel wide-band amplitude and phase calibration device for satellite communication
CN112350751B (en) * 2020-11-04 2022-07-08 中国电子科技集团公司第五十四研究所 Multi-channel wide-band amplitude and phase calibration device for satellite communication
CN113075701A (en) * 2021-03-25 2021-07-06 陕西凌云电器集团有限公司 Satellite signal receiving method and module with differential positioning and anti-interference functions
CN113075701B (en) * 2021-03-25 2023-04-18 陕西凌云电器集团有限公司 Satellite signal receiving method and module with differential positioning and anti-interference functions
CN114513228A (en) * 2022-04-19 2022-05-17 中国人民解放军海军工程大学 L-band high-speed frequency hopping data link non-cooperative interference cancellation device and method
CN114513228B (en) * 2022-04-19 2022-07-15 中国人民解放军海军工程大学 L-band high-speed frequency hopping data link non-cooperative interference cancellation device and method
CN117169801A (en) * 2023-11-02 2023-12-05 成都德辰博睿科技有限公司 Electromagnetic environment monitoring and calibrating system, method, device and medium
CN117169801B (en) * 2023-11-02 2024-03-19 成都德辰博睿科技有限公司 Electromagnetic environment monitoring and calibrating system, method, device and medium

Also Published As

Publication number Publication date
CN103117781B (en) 2016-04-13

Similar Documents

Publication Publication Date Title
CN103117781B (en) A kind of antenna array calibration method under complex electromagnetic environment and device thereof
EP2041897B1 (en) Method and arrangement for calibrating transmit paths of an antenna array
US10298338B2 (en) Method for evaluating quality of radio frequency signals for stellite navigation system
CN103901407B (en) C-band frequency-agile radar signal intercept receiving method
CN104601259B (en) Wireless communication receiver with i/q imbalance estimation and correction techniques
CN103869334B (en) The automatic identification of GNSS spacing wave distortion and disposal route
CN105656834A (en) Digital correction method for IQ channel mismatch of novel broadband receiver
WO2014006523A3 (en) Low sampling rate adaptation scheme for dual-band linearization
EP3014831B1 (en) Double transmission detection
CN104158777A (en) Range finder pulse interference suppression method combining orthogonal projection and circular beam forming
CN101881821B (en) Fractional order Fourier domain channelized receiving method
CN104901754B (en) Channel monitoring system based on channelization frequency spectrum perception
CN103338024B (en) The complementary Kalman filtering apparatus and method of time delay in antenna array
CN110213010A (en) A kind of unmanned plane detection system and method based on multi-channel radio frequency signal
CN105791194A (en) Synchronization method and system for anti-narrowband interference
CN103399301B (en) The receiving trap of a kind of broadband SAR signal and method of reseptance
CN105721078A (en) Automatic test method for general group delay ripple
CN103399030B (en) A kind of system and method for realizing answer signal phase difference detection under blank pipe 3/A pattern
CN101133655A (en) Tunable filter for tone detection
CN105549037A (en) High-precision satellite navigation broadband array signal generation method
CN102608663A (en) Interference canceller applied to detecting core quadrupole moment resonance signal
CN104967491B (en) Multichannel width tests system signal reception processing method mutually
CN202837534U (en) Distributed receiver and external active calibration device using same
CN108141242A (en) High-frequency front-end circuit and impedance matching methods
CN103701741B (en) Method for achieving signal modulation quality measurement and constellation diagram display

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160413

Termination date: 20220201

CF01 Termination of patent right due to non-payment of annual fee