CN118381570A - Automatic rapid calibration system for ultra-wideband multichannel - Google Patents
Automatic rapid calibration system for ultra-wideband multichannel Download PDFInfo
- Publication number
- CN118381570A CN118381570A CN202410817963.4A CN202410817963A CN118381570A CN 118381570 A CN118381570 A CN 118381570A CN 202410817963 A CN202410817963 A CN 202410817963A CN 118381570 A CN118381570 A CN 118381570A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- channel
- calibrated
- amplitude
- 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
- 230000001360 synchronised effect Effects 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 108091006146 Channels Proteins 0.000 claims description 158
- 238000000034 method Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 11
- 238000004088 simulation Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The invention provides an automatic rapid calibration system of ultra-wideband multichannel, which comprises: the system comprises a clock and control module, a calibration source generation module, an amplitude and phase calculation module, a compensation coefficient calculation module and a filter. The clock and control module provided by the invention generates the synchronous pulse signals of each stage through the second pulse signals, and can conveniently realize the automatic calibration function of the amplitude and phase consistency of the ultra-wideband multichannel according to the set time sequence as the synchronous signals of each stage of the whole automatic rapid calibration system. The automatic rapid calibration system provided by the invention is convenient and rapid to use, easy to realize and low in resource consumption, and is suitable for calibrating the amplitude and the phase of a large-scale array antenna.
Description
Technical Field
The invention belongs to the technical field of calibration, and particularly relates to an automatic rapid calibration system for ultra-wideband multiple channels.
Background
For ultra-wideband multichannel receiving systems, there are unavoidable amplitude and phase differences in the array antennas, LNAs (Low Noise Amplifier, low noise amplifiers), filters, frequency converters, digital-to-analog converters, and other analog devices involved. And, factors such as working environment, temperature variation and the like can influence the consistency of the amplitude and the phase of the device. Quantization errors introduced by analog signal digitization, and distortion errors caused by nonlinearity such as finite word length effect of an FPGA (Field Programmable GATE ARRAY ) chip in signal processing can influence the consistency of amplitude and phase among channels. The amplitude and phase inconsistencies between channels can have serious impact on subsequent signal processing, such as wideband DBF (Digital Beam Forming ), and even render the system inoperable.
Therefore, in engineering implementation engineering, a convenient, fast, easy-to-implement, and low-resource-consumption fast and automated calibration method is needed to calibrate amplitude and phase consistency between wideband signal channels.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an automatic rapid calibration system for ultra-wideband multiple channels. The technical problems to be solved by the invention are realized by the following technical scheme:
The invention provides an automatic rapid calibration system of ultra-wideband multichannel, comprising:
the clock and control module is used for acquiring a second pulse signal, generating an mth synchronous pulse signal Ts after receiving the second pulse signal, generating an mth synchronous pulse signal Td after a first preset time interval, and generating an mth synchronous pulse signal Te after a second preset time interval, wherein M is more than or equal to 1 and less than or equal to M, and M is the total number of the synchronous pulse signals Ts in one second pulse signal;
the calibration source generation module is used for switching the frequency to the mth frequency after receiving the mth synchronous pulse signal Ts so as to generate N paths of signals to be calibrated of the mth frequency;
The amplitude-phase calculation module is used for receiving an mth synchronous pulse signal Td, calculating the amplitude and the phase of each path of signal to be calibrated of the mth frequency at the same moment after receiving the mth synchronous pulse signal Td, completing the calculation of the amplitude and the phase of each path of signal to be calibrated of the mth frequency after receiving the mth synchronous pulse signal Te, and circularly executing the generation process of the amplitude and the phase of the signal to be calibrated until the calculation of the amplitude and the phase of all signals to be calibrated of all frequencies is completed after receiving the M synchronous pulse signals Te;
the compensation coefficient calculation module is used for obtaining a compensation coefficient according to the amplitude and the phase of the signal to be calibrated;
And the filter is used for updating the coefficient of the filter by utilizing the compensation coefficient to obtain an updated filter so as to calibrate each path of signal to be calibrated of each frequency by utilizing the updated filter to obtain a calibrated signal.
In one embodiment of the present invention, the clock and control module is further configured to set a start frequency, an end frequency, and a frequency interval;
The mth frequency is the sum of the initial frequency and m-1 frequency intervals;
The method for generating the N paths of signals to be calibrated of the mth frequency comprises the following steps of:
judging whether the mth frequency is smaller than or equal to the termination frequency, if so, switching the frequency to the mth frequency to generate an N-path signal to be calibrated of the mth frequency, and if not, stopping generating the signal to be calibrated.
In one embodiment of the present invention, calculating the amplitude and phase of each path of the signal to be calibrated for the mth frequency includes:
carrying out Fourier transform on the nth path of signals to be calibrated to obtain the nth path of signals to be calibrated in the frequency domain, wherein N is more than or equal to 1 and less than or equal to N;
And obtaining the phase and the amplitude of the nth signal to be calibrated according to the maximum amplitude of the nth signal to be calibrated.
In one embodiment of the present invention, obtaining the phase and the amplitude of the nth signal to be calibrated according to the maximum amplitude of the nth signal to be calibrated includes:
determining the maximum amplitude value of the nth path of the frequency domain signal to be calibrated, and obtaining the complex value of the maximum amplitude value;
and according to the complex value of the amplitude maximum value, the phase and the amplitude of the signal to be calibrated in the nth path are obtained.
In one embodiment of the present invention, obtaining the compensation coefficient according to the amplitude and the phase of the signal to be calibrated includes:
Obtaining time domain signals of M frequencies according to the amplitude and the phase of signals to be calibrated of the M frequencies, and obtaining input signals of N channels according to the time domain signals of the M frequencies;
selecting the first channel from the input signals of the N channels The input signal of the channel is used as the signal of the reference channel;
obtaining a filter frequency response of an e-th channel according to a frequency response of a reference channel obtained based on the signal of the reference channel and a frequency response of the e-th channel obtained based on the input signal of the e-th channel, ；
And obtaining the compensation coefficient of the e-th channel according to the filter frequency response of the e-th channel based on the least square method.
In one embodiment of the present invention, obtaining a time domain signal of M frequencies according to the amplitude and phase of signals to be calibrated of M frequencies, and obtaining an input signal of N channels according to the time domain signal of M frequencies, including:
Obtaining an nth time domain signal of an mth frequency according to the amplitude and the phase of the nth signal to be calibrated of the mth frequency, wherein the nth time domain signal of the mth frequency is expressed as:
；
Wherein, An nth time domain signal of an mth frequency,The amplitude of the nth signal to be calibrated for the mth frequency,The phase of the nth signal to be calibrated for the mth frequency,Is a natural index of the Chinese characters,In units of imaginary numbers,As a function of the discrete time variable,Is the mth frequency;
Obtaining an input signal of an nth channel according to an nth time domain signal of M frequencies, wherein the input signal of the nth channel is expressed as:
；
Wherein, An input signal for an nth channel;
And obtaining the input signals of the N paths of channels according to the input signals of the 1 st path of channels to the input signals of the N paths of channels.
In one embodiment of the present invention, the first channel is selected from the N-channel input signalsAn input signal of a channel as a signal of a reference channel, comprising:
selecting the first with minimum in-band fluctuation in the input signals of the N channels The input signal of the channel is used as the signal of the reference channel.
In one embodiment of the present invention, obtaining a filter frequency response of an e-th channel according to a frequency response of a reference channel obtained based on a signal of the reference channel and a frequency response of the e-th channel obtained based on an input signal of the e-th channel includes:
obtaining a frequency response of the reference channel according to the signal of the reference channel, wherein the frequency response of the reference channel is expressed as:
；
Wherein, For the frequency response of the reference channel, T is the FFT length,K is a discrete time variable, k is a FFT length variable, 1.ltoreq.k.ltoreq.T,Is the firstThe input signal of the channel is provided with a signal,Is a natural index of the Chinese characters,Is an imaginary unit;
obtaining the frequency response of the e-th channel according to the input signal of the e-th channel, wherein the frequency response of the e-th channel is expressed as:
；
Wherein, For the frequency response of the path e channel,An input signal for the e-th channel;
Obtaining a filter frequency response of the e-th channel according to the frequency response of the reference channel and the frequency response of the e-th channel, wherein the filter frequency response of the e-th channel is expressed as:
；
；
Wherein, For the filter frequency response of the e-th channel,Is the length of the filter.
In one embodiment of the present invention, the compensation coefficient of the e-th channel is expressed as:
；
；
Wherein, Is the compensation coefficient of the e-th channel,Is the conjugate transpose of the matrix I,For the filter frequency response of the e-th channel, T is the length of the FFT, k is the variable of the length of the FFT, k is more than or equal to 1 and less than or equal to T,For the length of the filter it is,In imaginary units.
In one embodiment of the present invention, the automated fast calibration system further includes a multi-channel data acquisition module, configured to acquire N paths of signals to be calibrated of the mth frequency, and transmit the N paths of signals to be calibrated to the amplitude-phase calculation module, and further configured to transmit the signals to be calibrated to the filter after updating coefficients of the filter with the compensation coefficients obtained by the compensation coefficient calculation module.
The invention has the beneficial effects that: the automatic rapid calibration system provided by the invention is provided with the synchronous pulse signal Ts, the synchronous pulse signal Td and the synchronous pulse signal Te, so that the broadband automatic calibration can follow the working time sequence set by the system. The clock and control module generates synchronous pulse signals of each stage through the second pulse signals, and the synchronous pulse signals are used as synchronous signals for working of each stage of the whole automatic rapid calibration system, and according to the set time sequence, the automatic calibration function of the ultra-wideband multichannel amplitude and phase consistency can be conveniently realized. And the automatic rapid calibration system provided by the invention calculates the amplitude and the phase of all signals to be calibrated through the amplitude and phase calculation module, then obtains the compensation coefficient of the filter based on the obtained amplitude and phase to update the coefficient of the filter, and calibrates the signals to be calibrated through the updated filter to obtain calibrated signals, thereby rapidly and accurately calibrating the amplitude and phase consistency among ultra-wideband signal channels. In engineering implementation, the form is convenient and quick to use, easy to implement and low in resource consumption, and is suitable for calibrating the amplitude and the phase of the large-scale array antenna.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
FIG. 1 is a schematic diagram of an automatic rapid calibration system for ultra wideband multiple channels according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of calibration operation time sequence of an ultra-wideband multichannel automatic rapid calibration system according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of an amplitude error before calibration according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a pre-calibration phase error according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of an amplitude error after calibration when the filter provided in the embodiment of the present invention is 15 th order;
Fig. 6 is a schematic diagram of a phase error after calibration when the filter provided in the embodiment of the invention is 15 th order.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but embodiments of the present invention are not limited thereto.
Example 1
Referring to fig. 1, fig. 1 is a schematic structural diagram of an automatic rapid calibration system for ultra-wideband multiple channels according to an embodiment of the present invention, where the automatic rapid calibration system for ultra-wideband multiple channels includes:
the clock and control module is used for acquiring a second pulse signal, generating an mth synchronous pulse signal Ts after receiving the second pulse signal, generating an mth synchronous pulse signal Td after a first preset time interval, and generating an mth synchronous pulse signal Te after a second preset time interval, wherein M is more than or equal to 1 and less than or equal to M, and M is the total number of the synchronous pulse signals Ts in one second pulse signal;
the calibration source generation module is used for switching the frequency to the mth frequency after receiving the mth synchronous pulse signal Ts so as to generate N paths of signals to be calibrated of the mth frequency;
The amplitude phase calculation module is used for receiving an mth synchronous pulse signal Td, calculating the amplitude and the phase of each path of signal to be calibrated of the mth frequency at the same moment after receiving the mth synchronous pulse signal Td, completing the calculation of the amplitude and the phase of each path of signal to be calibrated of the mth frequency after receiving the mth synchronous pulse signal Te, and circularly executing the generation process of the amplitude and the phase of the signal to be calibrated until the calculation of the amplitude and the phase of all signals to be calibrated of all frequencies is completed after receiving the M synchronous pulse signals Te;
the compensation coefficient calculation module is used for obtaining a compensation coefficient according to the amplitude and the phase of the signal to be calibrated;
And the filter is used for updating the coefficient of the filter by using the compensation coefficient to obtain an updated filter so as to calibrate the signal to be calibrated by using the updated filter to obtain a calibrated signal.
In this embodiment, referring to fig. 2, a second pulse signal (i.e. 1 pps) is provided from the outside of the system, then a synchronization pulse signal Ts, a synchronization pulse signal Td and a synchronization pulse signal Te are set based on the second pulse signal, and M synchronization pulse signals Ts, M synchronization pulse signals Td and M synchronization pulse signals Te are provided in the second pulse signal, where the synchronization pulse signal Ts is a trigger signal for generating a switching frequency of the module by the calibration source, the synchronization pulse signal Td is a trigger signal for starting to calculate the amplitude and phase of the signal to be calibrated currently, and the synchronization pulse signal Te is a signal for completing the amplitude and phase calculation of the current frequency.
In this embodiment, an initial frequency, a termination frequency, a frequency interval, and a residence time are set in the clock and control module, where the initial frequency is a frequency corresponding to the first synchronization pulse signal Ts, the termination frequency is a frequency for generating the synchronization pulse signal Ts by cutting off, the frequency interval is a frequency difference between two frequencies, the mth frequency corresponding to the mth synchronization pulse signal Ts is a sum of the initial frequency and M-1 frequency intervals, and the residence time is a holding time of one synchronization pulse signal, and the initial frequency, the termination frequency, and the frequency interval may be set according to practical requirements, and in this embodiment, the initial frequency is 1000MHz, the termination frequency is 2000MHz, the frequency interval is 10MHz, and the residence time is 20ms or 50ms, for example, when the total number M of frequencies is 50, the total calibration time is a calibration time of 1GHz bandwidth, and 1s is required.
In this embodiment, after the clock and control module sends the calibration command, the calibration program is started, the clock and control module will first receive the externally sent second pulse signal, when the clock and control module receives the second pulse signal, it will generate the first synchronization pulse signal Ts, then after the first preset time interval, it will generate the first synchronization pulse signal Td, when the amplitude and phase of the signal to be calibrated need to be calculated, the clock and control module will transmit the first synchronization pulse signal Td to the amplitude and phase calculation module to start to calculate the amplitude and phase of each signal to be calibrated of the first frequency, after the second preset time interval, it will generate the first synchronization pulse signal Te, the clock and control module will transmit the first synchronization pulse signal Te to the amplitude and phase calculation module, it will complete the calculation of the amplitude and phase of each signal to be calibrated of the first frequency, generate the corresponding amplitude matrix and phase matrix, and store them. After the amplitude and phase of each signal to be calibrated of the first frequency are generated, the second synchronous pulse signal Ts, the second synchronous pulse signal Td and the second synchronous pulse signal Te are sequentially generated, then the amplitude and phase of each signal to be calibrated of the second frequency are generated, and so on until the amplitude and phase of each signal to be calibrated of the Mth frequency are generated. The second preset time interval is greater than the first preset time interval, for example, 5ms, and is less than the time interval for generating the two adjacent synchronization pulse signals Ts, for example, 30ms.
In this embodiment, after receiving the mth synchronization pulse signal Ts sent by the clock and control module, the calibration source generating module switches the frequency to the mth frequency so as to generate N channels of signals to be calibrated at the mth frequency, where the number of channels N may be set according to the actual requirement, and in this embodiment, N is not limited specifically, and is, for example, 16 or 32.
In an alternative embodiment, switching the frequency to the mth frequency to generate the N-way signal to be calibrated for the mth frequency includes:
judging whether the mth frequency is smaller than or equal to the termination frequency, if so, switching the frequency to the mth frequency to generate N paths of signals to be calibrated of the mth frequency, and if not, stopping generating the signals to be calibrated.
That is, after receiving the mth synchronization pulse signal Ts, the calibration source generating module needs to determine whether the mth frequency corresponding to the mth synchronization pulse signal Ts has reached the termination frequency, if not, switch the frequency to the mth frequency to generate the signal to be calibrated, and if so, indicate that all the frequencies in the second pulse signal have been processed, and no further signal to be calibrated is needed.
In an alternative embodiment, the automated rapid calibration system may further include a multi-channel data acquisition module, where the multi-channel data acquisition module is configured to acquire N-th signals to be calibrated at the mth frequency, and transmit the N-th signals to be calibrated to the amplitude-phase calculation module.
That is, the present embodiment collects the N paths of signals to be calibrated generated by the calibration source generating module by using the multi-channel data collecting module, and then after the amplitude-phase calculating module receives the synchronization pulse signal Td, the collected N paths of signals to be calibrated are simultaneously sent to the amplitude-phase calculating module, so as to generate the amplitude and phase of the signals to be calibrated.
In an alternative embodiment, calculating the amplitude and phase of each signal to be calibrated for the mth frequency includes:
s1.1, carrying out Fourier transform on the nth path of signals to be calibrated to obtain the nth path of signals to be calibrated in the frequency domain, wherein N is more than or equal to 1 and less than or equal to N.
S1.2, obtaining the phase and the amplitude of the nth signal to be calibrated according to the maximum amplitude of the nth signal to be calibrated.
S1.21, determining the maximum amplitude value of the nth frequency domain signal to be calibrated, and obtaining the complex value of the maximum amplitude value, wherein the complex value of the maximum amplitude value is marked as a+bj, a is a real part, b is an imaginary part, and j is an imaginary unit.
S1.22, according to the complex value of the maximum amplitude value, the phase and the amplitude of the nth path of signal to be calibrated are obtained.
Here, the phase of the nth signal to be calibrated is expressed as:，， for the phase of the nth signal to be calibrated, An angle is the amplitude of the nth signal to be calibrated) Abs @ as a function of angle calculation) Taking the modulus value.
In this embodiment, after receiving the M synchronization pulse signals Te, the amplitude-phase calculation module calculates the amplitudes and phases of all the signals to be calibrated at all the frequencies, then generates the amplitudes and phases of all the signals to be calibrated at all the M frequencies, establishes an amplitude matrix from the amplitudes of all the signals to be calibrated at all the M frequencies, establishes a phase matrix from the phases of all the signals to be calibrated at all the M frequencies, and establishes an amplitude matrix a and a phase matrixExpressed as:
；
；
Wherein, The amplitude of the nth signal to be calibrated for the mth frequency,The phase of the nth signal to be calibrated for the mth frequency.
In an alternative embodiment, deriving the compensation factor from the amplitude and phase of the signal to be calibrated comprises:
s2.1, obtaining time domain signals of M frequencies according to the amplitude and the phase of signals to be calibrated of the M frequencies, and obtaining input signals of N channels according to the time domain signals of the M frequencies.
S2.11, obtaining an nth time domain signal of an mth frequency according to the amplitude and the phase of the nth signal to be calibrated of the mth frequency, wherein the nth time domain signal of the mth frequency is expressed as:
；
Wherein, An nth time domain signal of an mth frequency,Is a natural index of the Chinese characters,As a function of the discrete time variable,Is the mth frequency.
S2.12, obtaining an input signal of an nth channel according to an nth time domain signal of M frequencies, wherein the input signal of the nth channel is expressed as:
；
Wherein, Is the input signal of the nth channel.
S2.13, obtaining the input signals of the N paths of channels according to the input signals of the 1 st path of channels to the input signals of the N paths of channels.
Specifically, the input signals of the 1 st channel to the input signals of the N th channel can be sequentially obtained by circularly executing S2.11 to S2.12, so as to obtain the input signals of the N channels.
S2.2, selecting the first channel in the input signals of the N channelsThe input signal of the channel serves as the signal of the reference channel.
Specifically, the first with minimum in-band fluctuation is selected from the input signals of the N channelsThe input signal of the channel serves as the signal of the reference channel.
Here, the minimum in-band fluctuation is the minimum deviation of the amplitude, that is, the minimum difference between the maximum value and the minimum value of the amplitude. The present embodiment assumes the firstIn-band fluctuations of the input signal of the path channel are minimized.
S2.3, obtaining the filter frequency response of the e-th channel according to the frequency response of the reference channel obtained based on the signal of the reference channel and the frequency response of the e-th channel obtained based on the input signal of the e-th channel,。
S2.31, obtaining the frequency response of the reference channel according to the signal of the reference channel, wherein the frequency response of the reference channel is expressed as:
；
Wherein, For the frequency response of the reference channel, T is the length of the FFT (Fast Fourier Transfor, fast Fourier transform), k is the FFT length variable, 1.ltoreq.k.ltoreq.T,Is the firstInput signals of the channel.
S2.32, obtaining the frequency response of the e-th channel according to the input signal of the e-th channel, wherein the frequency response of the e-th channel is expressed as:
；
Wherein, For the frequency response of the path e channel,Is the input signal for the e-th channel.
S2.33, obtaining the filter frequency response of the e-th channel according to the frequency response of the reference channel and the frequency response of the e-th channel, wherein the filter frequency response of the e-th channel is expressed as:
；
；
Wherein, For the filter frequency response of the e-th channel,Is the length of the filter.
S2.4, based on a least square method, obtaining a compensation coefficient of the e-th channel according to the filter frequency response of the e-th channel.
Specifically, fitting is performed by adopting a least square method, so that the compensation coefficient of the e-th channel is obtained through the filter frequency response fitting of the e-th channel, and the compensation coefficient of the e-th channel is expressed as:
；
；
Wherein, Is the compensation coefficient of the e-th channel,Is the conjugate transpose of matrix I.
In this embodiment, the multi-channel data acquisition module is further configured to transmit the signal to be calibrated to the filter after the filter is updated by the compensation coefficient obtained by the compensation coefficient calculation module.
That is, after the compensation coefficient calculation module obtains the compensation coefficient of the corresponding channel, the obtained compensation coefficient is used to update the filter coefficient of the corresponding channel, and at this time, the multi-channel data acquisition module can transmit the signal to be calibrated to the channel corresponding to the updated filter, so as to calibrate the signal to be calibrated, and obtain the calibrated signal.
Optionally, the multi-channel data acquisition module, the amplitude and phase calculation module, the compensation coefficient calculation module and the filter are integrated on the same FPGA.
Optionally, the filter comprises a FIR (Finite Impulse Response, finite length unit impulse response) filter.
Simulation and result analysis:
Simulation conditions: through an antenna radiation mode, through radio frequency channels, ADC (Analog-to-Digital Converter) data acquisition is used, the scanning center frequency is 4.9 GHz-5.9 GHz, the frequency interval is 10MHz, amplitude phase information of different frequency points of each channel is extracted through acquired data, filter fitting is carried out, a least square method is adopted in a fitting algorithm, simulation results are shown in fig. 3, fig. 4, fig. 5 and fig. 6, fig. 3 is an amplitude error schematic diagram before calibration provided by the embodiment of the invention, fig. 4 is an amplitude error schematic diagram before calibration provided by the embodiment of the invention, fig. 5 is an amplitude error schematic diagram after calibration when a filter provided by the embodiment of the invention is 15 th order, fig. 6 is an amplitude error schematic diagram after calibration when the filter provided by the embodiment of the invention is 15 th order, and the phase error in a channel band is comprehensively known when the filter order number is 15 th order: not more than minus or plus ; Amplitude error in the channel band: not more than minus or plus。
The automatic rapid calibration system provided by the invention relates to broadband synchronous acquisition, broadband channel equalization and rapid realization of multichannel amplitude and phase test technology, and can rapidly realize the consistency of amplitude and phase among channels within a certain large bandwidth.
The automatic rapid calibration system provided by the invention is provided with the synchronous pulse signal Ts, the synchronous pulse signal Td and the synchronous pulse signal Te, so that the broadband automatic calibration can follow the working time sequence set by the system, the clock and control module generates synchronous pulse signals of each stage through the second pulse signal to serve as the synchronous signals working in each stage of the whole automatic rapid calibration system, and the automatic calibration function of ultra-broadband multichannel amplitude and phase consistency can be conveniently realized according to the set time sequence.
The method mainly adopts a pulse synchronization mode to carry out an automatic calibration flow, can effectively prescribe the working time of each functional module, ensures that each functional module is carried out regularly and in a beat, can effectively improve the working efficiency of the system, reduces the automatic calibration time, has great time advantage when calibrating a large-scale array antenna, can ensure that the system is orderly carried out according to the time sequence of the system, effectively avoids the time sequence blockage of the system in the working process, and reduces the probability of faults.
The automatic rapid calibration system provided by the invention calculates the amplitude and the phase of all signals to be calibrated through the amplitude and phase calculation module, then obtains the compensation coefficient of the filter based on the obtained amplitude and phase to update the coefficient of the filter, and calibrates the signals to be calibrated through the updated filter to obtain calibrated signals, thereby rapidly and accurately calibrating the amplitude and phase consistency among ultra-wideband signal channels. In engineering implementation, the form is convenient and quick to use, easy to implement and low in resource consumption, and is suitable for calibrating the amplitude and the phase of the large-scale array antenna.
In the description of the present invention, it should be understood that the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present invention, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, one skilled in the art can combine and combine the different embodiments or examples described in this specification.
The foregoing is a further detailed description of the invention in connection with the preferred embodiments, and it is not intended that the invention be limited to the specific embodiments described. It will be apparent to those skilled in the art that several simple deductions or substitutions may be made without departing from the spirit of the invention, and these should be considered to be within the scope of the invention.
Claims (10)
1. An automated rapid calibration system for ultra-wideband multiple channels, the automated rapid calibration system comprising:
the clock and control module is used for acquiring a second pulse signal, generating an mth synchronous pulse signal Ts after receiving the second pulse signal, generating an mth synchronous pulse signal Td after a first preset time interval, and generating an mth synchronous pulse signal Te after a second preset time interval, wherein M is more than or equal to 1 and less than or equal to M, and M is the total number of the synchronous pulse signals Ts in one second pulse signal;
the calibration source generation module is used for switching the frequency to the mth frequency after receiving the mth synchronous pulse signal Ts so as to generate N paths of signals to be calibrated of the mth frequency;
The amplitude-phase calculation module is used for receiving an mth synchronous pulse signal Td, calculating the amplitude and the phase of each path of signal to be calibrated of the mth frequency at the same moment after receiving the mth synchronous pulse signal Td, completing the calculation of the amplitude and the phase of each path of signal to be calibrated of the mth frequency after receiving the mth synchronous pulse signal Te, and circularly executing the generation process of the amplitude and the phase of the signal to be calibrated until the calculation of the amplitude and the phase of all signals to be calibrated of all frequencies is completed after receiving the M synchronous pulse signals Te;
the compensation coefficient calculation module is used for obtaining a compensation coefficient according to the amplitude and the phase of the signal to be calibrated;
And the filter is used for updating the coefficient of the filter by utilizing the compensation coefficient to obtain an updated filter so as to calibrate each path of signal to be calibrated of each frequency by utilizing the updated filter to obtain a calibrated signal.
2. The automated rapid calibration system of claim 1, wherein the clock and control module is further configured to set a start frequency, an end frequency, a frequency interval;
The mth frequency is the sum of the initial frequency and m-1 frequency intervals;
The method for generating the N paths of signals to be calibrated of the mth frequency comprises the following steps of:
judging whether the mth frequency is smaller than or equal to the termination frequency, if so, switching the frequency to the mth frequency to generate an N-path signal to be calibrated of the mth frequency, and if not, stopping generating the signal to be calibrated.
3. The automated rapid calibration system of claim 1, wherein calculating the amplitude and phase of each of the signals to be calibrated for the mth frequency comprises:
carrying out Fourier transform on the nth path of signals to be calibrated to obtain the nth path of signals to be calibrated in the frequency domain, wherein N is more than or equal to 1 and less than or equal to N;
And obtaining the phase and the amplitude of the nth signal to be calibrated according to the maximum amplitude of the nth signal to be calibrated.
4. An automated fast calibration system according to claim 3, wherein deriving the phase and amplitude of the nth signal to be calibrated from the maximum amplitude of the nth signal to be calibrated comprises:
determining the maximum amplitude value of the nth path of the frequency domain signal to be calibrated, and obtaining the complex value of the maximum amplitude value;
and according to the complex value of the amplitude maximum value, the phase and the amplitude of the signal to be calibrated in the nth path are obtained.
5. The automated rapid calibration system of claim 1, wherein deriving compensation coefficients from the amplitude and phase of the signal to be calibrated comprises:
Obtaining time domain signals of M frequencies according to the amplitude and the phase of signals to be calibrated of the M frequencies, and obtaining input signals of N channels according to the time domain signals of the M frequencies;
selecting the first channel from the input signals of the N channels The input signal of the channel is used as the signal of the reference channel;
obtaining a filter frequency response of an e-th channel according to a frequency response of a reference channel obtained based on the signal of the reference channel and a frequency response of the e-th channel obtained based on the input signal of the e-th channel, ；
And obtaining the compensation coefficient of the e-th channel according to the filter frequency response of the e-th channel based on the least square method.
6. The automated rapid calibration system of claim 5, wherein obtaining M frequency time domain signals from the amplitude and phase of M frequency signals to be calibrated, and obtaining N channel input signals from the M frequency time domain signals, comprises:
Obtaining an nth time domain signal of an mth frequency according to the amplitude and the phase of the nth signal to be calibrated of the mth frequency, wherein the nth time domain signal of the mth frequency is expressed as:
；
Wherein, An nth time domain signal of an mth frequency,The amplitude of the nth signal to be calibrated for the mth frequency,The phase of the nth signal to be calibrated for the mth frequency,Is a natural index of the Chinese characters,In units of imaginary numbers,As a function of the discrete time variable,Is the mth frequency;
Obtaining an input signal of an nth channel according to an nth time domain signal of M frequencies, wherein the input signal of the nth channel is expressed as:
；
Wherein, An input signal for an nth channel;
And obtaining the input signals of the N paths of channels according to the input signals of the 1 st path of channels to the input signals of the N paths of channels.
7. The automated rapid calibration system of claim 5, wherein a first selected among the N-way input signals isAn input signal of a channel as a signal of a reference channel, comprising:
selecting the first with minimum in-band fluctuation in the input signals of the N channels The input signal of the channel is used as the signal of the reference channel.
8. The automated rapid calibration system of claim 5, wherein deriving the filter frequency response for the e-th channel from the frequency response for the reference channel derived based on the signal for the reference channel and the frequency response for the e-th channel derived based on the input signal for the e-th channel comprises:
obtaining a frequency response of the reference channel according to the signal of the reference channel, wherein the frequency response of the reference channel is expressed as:
；
Wherein, For the frequency response of the reference channel, T is the FFT length,K is a discrete time variable, k is a FFT length variable, 1.ltoreq.k.ltoreq.T,Is the firstThe input signal of the channel is provided with a signal,Is a natural index of the Chinese characters,Is an imaginary unit;
obtaining the frequency response of the e-th channel according to the input signal of the e-th channel, wherein the frequency response of the e-th channel is expressed as:
；
Wherein, For the frequency response of the path e channel,An input signal for the e-th channel;
Obtaining a filter frequency response of the e-th channel according to the frequency response of the reference channel and the frequency response of the e-th channel, wherein the filter frequency response of the e-th channel is expressed as:
；
；
Wherein, For the filter frequency response of the e-th channel,Is the length of the filter.
9. The automated rapid calibration system of claim 5, wherein the compensation coefficient of the e-th path is expressed as:
；
；
Wherein, Is the compensation coefficient of the e-th channel,Is the conjugate transpose of the matrix I,For the filter frequency response of the e-th channel, T is the length of the FFT, k is the variable of the length of the FFT, k is more than or equal to 1 and less than or equal to T,For the length of the filter it is,In imaginary units.
10. The automated rapid calibration system of claim 1, further comprising a multi-channel data acquisition module for acquiring N-way signals to be calibrated for the mth frequency and transmitting the N-way signals to be calibrated to the amplitude phase calculation module, and further for transmitting the signals to be calibrated to the filter after updating coefficients of the filter with the compensation coefficients obtained by the compensation coefficient calculation module.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410817963.4A CN118381570B (en) | 2024-06-24 | 2024-06-24 | Automatic rapid calibration system for ultra-wideband multichannel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410817963.4A CN118381570B (en) | 2024-06-24 | 2024-06-24 | Automatic rapid calibration system for ultra-wideband multichannel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118381570A true CN118381570A (en) | 2024-07-23 |
| CN118381570B CN118381570B (en) | 2024-09-10 |
Family
ID=91912682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410817963.4A Active CN118381570B (en) | 2024-06-24 | 2024-06-24 | Automatic rapid calibration system for ultra-wideband multichannel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118381570B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130329841A1 (en) * | 2012-06-08 | 2013-12-12 | Deere & Company | High Frequency Signal Receiver with Self-Calibrated Group Delay Compensation |
| US20150030104A1 (en) * | 2013-07-26 | 2015-01-29 | Research & Business Foundation Sungkyunkwan University | Apparatus and method for compensating output signals of magnetic encoder using digital phase-locked loop |
| CN104779989A (en) * | 2015-05-11 | 2015-07-15 | 重庆大学 | Boardband array correcting filter coefficient calculation method |
| KR101554287B1 (en) * | 2014-05-23 | 2015-09-18 | 국방과학연구소 | Apparatus for correcting amplitude and phase for multi-mode mono pulse antenna using pilot signals |
| CN110031809A (en) * | 2019-04-15 | 2019-07-19 | 北京电子工程总体研究所 | A kind of more receiving channel amplitude-phase consistency Calibration Methods of radar and readable storage medium storing program for executing |
| US20190238176A1 (en) * | 2018-02-01 | 2019-08-01 | Thales | Method and system for calibrating a radiofrequency multichannel subsystem of a telecommunications payload |
| CN110221268A (en) * | 2019-06-21 | 2019-09-10 | 西北核技术研究院 | A kind of method phase-magnitude calibration filter and promote the non-coherent radar signal processing performance of short pulse |
| CN111211797A (en) * | 2020-01-14 | 2020-05-29 | 安徽大学 | Receiver compensation method and compensation system and electronic equipment thereof |
| CN111257883A (en) * | 2020-02-13 | 2020-06-09 | 航天新气象科技有限公司 | Calibration device and method for phased array weather radar system |
| CN113884996A (en) * | 2021-09-18 | 2022-01-04 | 北京遥感设备研究所 | Multi-time-width pulse signal correction method and system of special test equipment for radar |
| CN114185008A (en) * | 2021-09-30 | 2022-03-15 | 西安电子科技大学 | System and method for compensating amplitude-phase error of receiving channel of narrow-band digital array radar system |
| CN116436537A (en) * | 2022-11-17 | 2023-07-14 | 北京理工大学前沿技术研究院 | Digital TR (transmitter and receiver) component amplitude and phase calibration method |
-
2024
- 2024-06-24 CN CN202410817963.4A patent/CN118381570B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130329841A1 (en) * | 2012-06-08 | 2013-12-12 | Deere & Company | High Frequency Signal Receiver with Self-Calibrated Group Delay Compensation |
| US20150030104A1 (en) * | 2013-07-26 | 2015-01-29 | Research & Business Foundation Sungkyunkwan University | Apparatus and method for compensating output signals of magnetic encoder using digital phase-locked loop |
| KR101554287B1 (en) * | 2014-05-23 | 2015-09-18 | 국방과학연구소 | Apparatus for correcting amplitude and phase for multi-mode mono pulse antenna using pilot signals |
| CN104779989A (en) * | 2015-05-11 | 2015-07-15 | 重庆大学 | Boardband array correcting filter coefficient calculation method |
| US20190238176A1 (en) * | 2018-02-01 | 2019-08-01 | Thales | Method and system for calibrating a radiofrequency multichannel subsystem of a telecommunications payload |
| CN110031809A (en) * | 2019-04-15 | 2019-07-19 | 北京电子工程总体研究所 | A kind of more receiving channel amplitude-phase consistency Calibration Methods of radar and readable storage medium storing program for executing |
| CN110221268A (en) * | 2019-06-21 | 2019-09-10 | 西北核技术研究院 | A kind of method phase-magnitude calibration filter and promote the non-coherent radar signal processing performance of short pulse |
| CN111211797A (en) * | 2020-01-14 | 2020-05-29 | 安徽大学 | Receiver compensation method and compensation system and electronic equipment thereof |
| CN111257883A (en) * | 2020-02-13 | 2020-06-09 | 航天新气象科技有限公司 | Calibration device and method for phased array weather radar system |
| CN113884996A (en) * | 2021-09-18 | 2022-01-04 | 北京遥感设备研究所 | Multi-time-width pulse signal correction method and system of special test equipment for radar |
| CN114185008A (en) * | 2021-09-30 | 2022-03-15 | 西安电子科技大学 | System and method for compensating amplitude-phase error of receiving channel of narrow-band digital array radar system |
| CN116436537A (en) * | 2022-11-17 | 2023-07-14 | 北京理工大学前沿技术研究院 | Digital TR (transmitter and receiver) component amplitude and phase calibration method |
Non-Patent Citations (5)
| Title |
|---|
| 张治西;刘长青;王震;: "RCS动态测量雷达幅度相位校准技术", 现代雷达, no. 01, 15 January 2010 (2010-01-15) * |
| 张诚;陈励军;: "多路径信号建模在换能器校准中的应用", 电声技术, no. 04, 17 April 2007 (2007-04-17) * |
| 彭涛;: "自适应抗干扰系统中多通道幅相校准的工程实现", 舰船电子工程, no. 05, 20 May 2019 (2019-05-20) * |
| 杜俊蛟;曹旭东;李宝森;孟华东;: "基于FPGA的多通道雷达接收机幅相不一致校正", 电光与控制, no. 01, 1 January 2011 (2011-01-01) * |
| 黄蕾;王萍;周哲帅;董航;王磊;: "相控阵快速校准系统研究", 电子科学技术, no. 04, 10 July 2016 (2016-07-10) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118381570B (en) | 2024-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114185008B (en) | System and method for compensating amplitude and phase errors of receiving channel of narrow-band digital array radar system | |
| CN105548978A (en) | Active reflection coefficient test system and method of phased array antenna | |
| CN108614247A (en) | A kind of secondary radar channel calibration method | |
| EP3675365A1 (en) | Linear and non-linear calibration for time interleaved digital-to-analog converter | |
| US7668509B2 (en) | Frequency selective leveling loop for multi-signal phased array transmitters | |
| CN110824466A (en) | Multi-target tracking system and DBF channel calibration FPGA implementation method thereof | |
| US8788234B2 (en) | Method of calibrating interleaved digitizer channels | |
| CN107796991A (en) | The active standing wave automatic testing equipment of phased array antenna and method | |
| CN118381570B (en) | Automatic rapid calibration system for ultra-wideband multichannel | |
| US20040108953A1 (en) | Method and device for the calibration-equalization of a reception system | |
| CN106256044B (en) | A kind of phased array calibration method and phased array calibrate circuit | |
| US7649492B2 (en) | Systems and methods for providing delayed signals | |
| CN117388811A (en) | Multi-channel calibration equipment and calibration method of semi-physical injection simulation system | |
| CN117713960A (en) | Channel calibration and equalization method for broadband zero intermediate frequency receiving array | |
| Naldi et al. | Development of a new digital signal processing platform for the Square Kilometre Array | |
| CN114545078B (en) | Intermediate frequency signal processing system and processing method | |
| CN110958062A (en) | Radio frequency multichannel calibration system and method | |
| CN114639957A (en) | Digital-analog mixed multi-beam forming vehicle-mounted device and phase-shifting wave control method thereof | |
| Nuteson et al. | Digital beamforming and calibration for smart antennas using real-time FPGA processing | |
| EP3384312A1 (en) | Array and module calibration with delay line | |
| Fischer et al. | X-Band HRWS demonstrator: Digital Beamforming test results | |
| CN107918070B (en) | Digital T/R assembly test system and transmitting and receiving state test method thereof | |
| CN217212975U (en) | Digital source single receiving network vector analysis system | |
| CN111913160B (en) | Correction coefficient acquisition method and array receiving system | |
| Tidd et al. | Sequential beamspace DOA estimation |
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 |