Calibration method, apparatus, medium, and device for ultra-wideband unevenness and iq imbalance

Abstract

The application relates to a method, device, medium and equipment for calibrating ultra-wideband unevenness and IQ imbalance, wherein the method comprises the following steps: receiving a single-tone signal at equal frequency intervals, and performing filtering processing on the single-tone signal to obtain a first single-tone signal; performing windowing processing and frequency domain conversion on the first single-tone signal to obtain a second single-tone signal; obtaining peak points of the second single-tone signal IQ two-way in the corresponding frequency domain, determining IQ phase imbalance information and IQ amplitude imbalance information according to the phases and amplitudes of the peak points; establishing an imbalance signal according to the IQ phase imbalance information and the IQ amplitude imbalance information; establishing a local standard signal, calculating the coefficient of an equivalent filter based on the imbalance signal and the local standard signal by using a least square method, and performing filtering on the single-tone signal according to the coefficient of the equivalent filter, so as to calibrate the ultra-wideband unevenness and the IQ imbalance. The application has the effect of effectively calibrating the ultra-wideband unevenness and the IQ imbalance of a receiver.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of communication technology, in particular to a method, device, medium and equipment for calibrating ultra-wideband unevenness and IQ imbalance. BACKGROUND

[0002] Ultra-wideband is a new wireless communication technology, which modulates impulse with very steep rising and falling time, so that the signal has a bandwidth of GHz. IQ imbalance is due to the limitations of analog devices in the up-conversion and down-conversion receivers, which causes the amplitude and phase imbalance of the signal in the (I, In-phase) and (Q, Quadrature) paths. IQ imbalance varies with frequency in the receive band. Since the wideband of ultra-wideband signal is large, IQ imbalance and unevenness will be further deteriorated, so there is an urgent need for a method for calibrating ultra-wideband unevenness and IQ imbalance of the receiver.

[0003] Currently, the commonly used calibration method is to inject from the radio frequency end of the receiver through a power divider and a directional coupler, or a far / near field source can also be used. The advantage of this method is easy to implement, but the disadvantage is that it does not include correction of the antenna array and the feed line, and the power divider and directional coupler will introduce additional amplitude and phase errors, which will ultimately affect the improvement of the correction performance. Another calibration method is to use a far field source method, which can overcome the above-mentioned shortcomings, but at the same time, due to the existence of surrounding objects, scattering and multipath propagation effects will occur, resulting in ineffective calibration of ultra-wideband unevenness and IQ imbalance. SUMMARY

[0004] In order to effectively calibrate the ultra-wideband unevenness and IQ imbalance of the receiver, the present application provides a method, device, medium and equipment for calibrating ultra-wideband unevenness and IQ imbalance.

[0005] In a first aspect of the present application, a method for calibrating ultra-wideband unevenness and IQ imbalance is provided, which specifically includes: receiving a single tone signal at equal frequency intervals, and filtering the single tone signal to obtain a first single tone signal, the single tone signal being a signal corresponding to a sub-band, and a plurality of sub-bands forming an ultra-wideband;

[0006] The first single tone signal is windowed and frequency domain transformed to obtain a second single tone signal;

[0007] The peak points of the IQ two paths of the second single tone signal in the corresponding frequency domain are obtained, and the IQ phase imbalance information and the IQ amplitude imbalance information are determined according to the phases and amplitudes of the peak points;

[0008] The imbalance signal is established according to the IQ phase imbalance information and the IQ amplitude imbalance information.

[0009] establishing a local standard signal, calculating coefficients of an equivalent filter by a least square method based on the unbalance signal and the local standard signal, and filtering the single tone signal according to the coefficients of the equivalent filter to calibrate the ultra-wideband unevenness and IQ imbalance.

[0010] By adopting the technical scheme, after receiving the multiple single tone signals, the first single tone signal is obtained through filtering processing, so that the noise other than the single tone signal is filtered out, and the first single tone signal is subjected to windowing processing and frequency domain conversion to obtain the second single tone signal, so that the first single tone signal has higher definition and is easier to show the characteristics, thereby facilitating subsequent signal analysis. Then, the IQ two-way phase and amplitude of the peak points in the corresponding frequency domain of the second single tone signal are compared to obtain the IQ phase imbalance information and the IQ amplitude imbalance information, that is, the difference information of the IQ two-way phase and the difference information of the IQ two-way amplitude. Then, the errors of the IQ phase and amplitude are determined according to the IQ phase imbalance information and the IQ amplitude imbalance information, and the unbalance signal is determined according to the errors of the IQ phase and amplitude. Finally, the coefficients of the equivalent filter are calculated by the least square method, so that the frequency response of the channel corresponding to the unbalance signal is equalized to be consistent with the local standard signal under the action of the equivalent filter, thereby suppressing the image interference generated by the single tone signal, and effectively calibrating the ultra-wideband unevenness and IQ imbalance.

[0011] Optionally, the windowing processing and the frequency domain conversion on the first single tone signal to obtain the second single tone signal specifically include:

[0012] The real part signal and the imaginary part signal of the first single tone signal are subjected to windowing processing respectively to obtain a windowed real part signal and a windowed imaginary part signal.

[0013] The windowed real part signal and the windowed imaginary part signal are subjected to frequency domain conversion respectively to obtain the second single tone signal.

[0014] By adopting the technical scheme, since the first single tone signal has a corresponding complex signal, in order to facilitate the processing of the first single tone signal, the real part signal and the imaginary part signal of the complex signal corresponding to the first single tone signal are subjected to windowing processing and frequency domain conversion respectively, so as to finally obtain the second single tone signal, thereby facilitating the subsequent analysis of the unbalance characteristics of the IQ two-way signal of the second single tone signal.

[0015] Optionally, the obtaining of the peak points of the IQ two-way of the second single tone signal in the corresponding frequency domain and the determination of the IQ phase imbalance information and the IQ amplitude imbalance information according to the phase and amplitude of each peak point specifically include:

[0016] The peak points of the IQ two-way of the second single tone signal in the corresponding frequency domain are obtained.

[0017] Subtracting the phase of the peak point corresponding to the I path from the phase of the peak point corresponding to the Q path and subtracting pi / 2, the IQ phase imbalance information is obtained; the ratio of the amplitude of the peak point corresponding to the I path to the amplitude of the peak point corresponding to the Q path is determined as the IQ amplitude imbalance information.

[0018] By adopting the technical scheme, after the second single-tone signal is determined, the peak point of the I path signal of the second single-tone signal and the peak point of the Q path signal are determined respectively. Since the phase of the Q path signal itself is pi / 2 different from the phase of the I path signal, when the phase difference between the Q path signal and the I path signal is calculated, pi / 2 needs to be subtracted additionally, and finally the IQ phase imbalance information is obtained. In addition, the amplitude of the peak point of the I path signal is divided by the amplitude of the peak point corresponding to the Q path, and the ratio of the two is used to represent the IQ amplitude imbalance, and is determined as the IQ amplitude imbalance information, thereby facilitating the establishment of the subsequent imbalance signal.

[0019] Optionally, the imbalance signal is established according to the IQ phase imbalance information and the IQ amplitude imbalance information, and specifically includes:

[0020] The IQ phase imbalance information is recorded in a preset ERR_P vector list to obtain a phase error table, and the IQ amplitude imbalance information is recorded in a preset ERR_A vector list to obtain an amplitude error table, the ERR_P vector list includes standard IQ phase imbalance information, and the ERR_A vector list includes standard IQ amplitude imbalance information; and the imbalance signal is obtained according to the phase error table and the amplitude error table.

[0021] A plurality of groups of the imbalance signal are merged to obtain an error signal amplitude and phase vector table, and a frequency domain corresponding to the imbalance signal is generated according to the error signal amplitude and phase vector table.

[0022] By adopting the technical scheme, after the IQ phase imbalance information and the IQ phase imbalance information are determined, the IQ phase imbalance information is recorded in the ERR_P vector list, so that the IQ phase imbalance information and the standard IQ phase imbalance information constitute a phase error table of the second single-tone signal; similarly, the IQ amplitude imbalance information is recorded in the ERR_A vector list to determine an amplitude error table. Then, the amplitude error table and the phase error table are input into a variable function to restore the imbalance signal and the corresponding frequency domain under the phase and amplitude error, so that the imbalance characteristics of the single-tone signal can be more purely reflected, and further analysis of the imbalance signal is facilitated.

[0023] Optionally, the local standard signal is established, and specifically includes:

[0024] A sampling rate and a frequency point list of the single-tone signal are obtained.

[0025] According to the sampling rate and the frequency point list, a local standard signal is established by matlab;

[0026] A plurality of local standard signals are merged to obtain an original reference signal amplitude and phase vector table, and a frequency domain corresponding to each local standard signal is generated according to the original reference signal amplitude and phase vector table.

[0027] According to the sampling rate and the frequency point list, a local standard signal is established by matlab;

[0028] Optionally, the coefficient of the equivalent filter is calculated by the least square method based on the unbalanced signal and the local standard signal, and specifically includes:

[0029] The channel corresponding to the local standard signal is determined as a reference channel, and the channel corresponding to the unbalanced signal is determined as a mismatch channel;

[0030] The all-pass linear phase shift network frequency response H ref is determined, and the first formula H i (jω) = C ref (jω) * H ref / C i (jω) is used to calculate the frequency response of the i-th equalizer to be solved, wherein C ref (jω) is the inherent frequency response of the reference channel, C i (jω) is the inherent frequency response of any mismatch channel, i is the number of subbands, H i (jω) is the frequency response of the i-th equalizer to be solved.

[0031] The frequency factor matrix is obtained, and the conjugate transpose of the frequency factor matrix is obtained to obtain a weighting matrix, the weighting matrix including the fitting error weight of each frequency point in the frequency factor matrix;

[0032] E(jω) is set, E(jω) is approximated to H i (jω) by the least square method, and the coefficient of the equivalent filter is calculated according to the frequency factor matrix, H i (jω) and the weighting matrix through a preset second formula, wherein E(jω) is the frequency response of the equivalent filter.

[0033] By adopting the technical scheme, after the local standard signal and the unbalanced signal are determined, a channel of the local standard signal is taken as a reference channel, a channel of the unbalanced signal is taken as a mismatch channel, then the inherent frequency response of the reference channel is divided by the inherent frequency response of the mismatch channel, and then multiplied by the all-pass linear phase shift network frequency response, so that the frequency response of the i-th equalizer to be solved is obtained, then the frequency response of the equivalent filter approximates the frequency response of the i-th equalizer to be solved, the interference of the measurement noise is avoided, finally, the frequency factor matrix corresponding to the single tone signal, the frequency response of the i-th equalizer to be solved and the weighting matrix are input into the second formula to calculate the coefficient of the equivalent filter, so that the equivalent filter plays the role of the equalizer, the frequency responses of other channels are equalized to be consistent with the frequency response of the pre-reference channel, and then the ultra-wideband unevenness and IQ imbalance are effectively calibrated.

[0034] Optionally, the second formula is: FIR_h=inv(W*A)*(W*b), wherein b is the value of H i (jω), FIR_h is the coefficient of the equivalent filter, W is the weighting matrix, and A is the frequency factor matrix.

[0035] By adopting the technical scheme, based on the second formula, the product of the weighting matrix and the frequency factor matrix is matrix-inverted to obtain inv(W*A), then the weighting matrix and the frequency response of the i-th equalizer to be solved are multiplied, finally, the result after multiplication is multiplied by inv(W*A), the final matrix calculation fitting of the coefficient of the equivalent filter is completed, and the coefficient FIR_h of the equivalent filter capable of calibrating the ultra-wideband unevenness and IQ imbalance is obtained.

[0036] In a second aspect of the present application, an ultra-wideband unevenness and IQ imbalance calibration device is provided, which specifically comprises: a single tone acquisition module configured to receive single tone signals at equal frequency intervals and perform filtering processing on the single tone signals to obtain first single tone signals, wherein the single tone signals are signals corresponding to subbands, and a plurality of the subbands form an ultra-wideband;

[0037] A single tone processing module is configured to perform windowing processing and frequency domain transformation on the first single tone signals to obtain second single tone signals; an information determination module is configured to obtain peak points of IQ two paths of the second single tone signals in a corresponding frequency domain, and determine IQ phase imbalance information and IQ amplitude imbalance information according to phases and amplitudes of the peak points;

[0038] A signal establishment module is configured to establish unbalanced signals according to the IQ phase imbalance information and the IQ amplitude imbalance information;

[0039] The coefficient determination module is configured to establish a local standard signal, calculate the coefficient of the equivalent filter by the least square method based on the unbalanced signal and the local standard signal, and filter the single-tone signal according to the coefficient of the equivalent filter to calibrate the ultra-wideband unevenness and IQ imbalance.

[0040] By adopting the technical scheme, the single-tone acquisition module acquires the received single-tone signal, filters the single-tone signal to obtain a first single-tone signal, and then the single-tone processing module performs windowing processing and frequency domain conversion on the first single-tone signal to obtain a second single-tone signal. Then, the information determination module acquires the peak points of the IQ two paths of the second single-tone signal in the corresponding frequency domain, determines the IQ phase imbalance information and the IQ amplitude imbalance information according to the phases and amplitudes of the peak points. The signal establishment module establishes the unbalanced signal according to the IQ phase imbalance information and the IQ amplitude imbalance information. Finally, the coefficient determination module calculates the coefficient of the equivalent filter according to the local standard signal and the unbalanced signal by the least square method, so that the coefficient of the equivalent filter is used to effectively calibrate the ultra-wideband unevenness and IQ imbalance.

[0041] In summary, the present application has at least one of the following beneficial technical effects:

[0042] After receiving a plurality of single-tone signals, the first single-tone signal is obtained by filtering processing, so as to filter out noise other than the single-tone signal, and the first single-tone signal is windowed and frequency domain converted to obtain a second single-tone signal, so that the first single-tone signal has high clarity and is easy to see the characteristics, facilitating subsequent signal analysis. Then, the phases and amplitudes of the peak points of the IQ two paths of the second single-tone signal in the corresponding frequency domain are compared to obtain the IQ phase imbalance information and the IQ amplitude imbalance information. Then, the errors of the IQ phase and amplitude are determined according to the IQ phase imbalance information and the IQ amplitude imbalance information, and then the unbalanced signal is determined according to the errors of the IQ phase and amplitude. Finally, the coefficient of the equivalent filter is calculated by the least square method, so that under the action of the equivalent filter, the frequency response of the channel corresponding to the unbalanced signal is equalized to be consistent with the local standard signal, thereby suppressing the image interference generated by the single-tone signal, and effectively calibrating the ultra-wideband unevenness and IQ imbalance. BRIEF DESCRIPTION OF DRAWINGS

[0043] Figure 1 is a flowchart of an ultra-wideband unevenness and IQ imbalance calibration method provided by an embodiment of the present application;

[0044] Figure 2 is a spectrum diagram of 25 left single-tone signals superimposed by a receiver according to an embodiment of the present application;

[0045] Figure 3is a spectrum diagram of 25 right single tone signals superimposed after being received by a receiver provided by an embodiment of the present application.

[0046] Figure 4 is a schematic diagram of an equalizer calibration principle provided by an embodiment of the present application.

[0047] Figure 5 is a flowchart of another calibration method of ultra-wideband unevenness and IQ imbalance provided by an embodiment of the present application.

[0048] Figure 6 is a phase error table provided by an embodiment of the present application.

[0049] Figure 7 is an amplitude error table provided by an embodiment of the present application.

[0050] Figure 8 is a flowchart of another calibration method of ultra-wideband unevenness and IQ imbalance provided by an embodiment of the present application.

[0051] Figure 9 is a schematic diagram of real and imaginary parts of a coefficient corresponding to an equivalent filter provided by an embodiment of the present application.

[0052] Figure 10 is a spectrum diagram of 25 left single tone signals filtered by an equivalent filter provided by an embodiment of the present application.

[0053] Figure 11 is a spectrum diagram of 25 right single tone signals filtered by an equivalent filter provided by an embodiment of the present application.

[0054] Figure 12 is a structural schematic diagram of a calibration device of ultra-wideband unevenness and IQ imbalance provided by an embodiment of the present application.

[0055] Legend: 11, single tone acquisition module; 12, single tone processing module; 13, information determination module; 14, signal establishment module; 15, coefficient determination module. DETAILED DESCRIPTION

[0056] In order to enable persons skilled in the art to better understand the technical solutions in the present specification, the technical solutions in the present specification will be described clearly and completely in conjunction with the drawings in the embodiments of the present specification. Obviously, the described embodiments are only some of the embodiments of the present application, not all the embodiments.

[0057] In the description of the embodiments of the present application, the words "exemplary", "for example", or "e.g." are used to mean "an example of" or "an example, only". Any embodiment or design solution described as "exemplary", "for example", or "e.g." in the embodiments of the present application should not be construed as being more advantageous or superior than other embodiments or design solutions. In fact, the use of the words "exemplary", "for example", or "e.g." is intended to present related concepts in a specific manner.

[0058] In the description of the embodiments of the present application, the term "and / or" is only used to describe the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B, which can represent three cases of existence of A alone, existence of B alone, and existence of A and B at the same time. In addition, unless otherwise specified, the term "multiple" means two or more. For example, multiple systems refer to two or more systems, and multiple screen terminals refer to two or more screen terminals. In addition, the terms "first" and "second" are only used for description purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the indicated technical features. Therefore, the features defined with "first" and "second" can explicitly or implicitly include one or more features. The terms "include", "contain", "have" and their variants mean "include but are not limited to", unless otherwise specifically emphasized.

[0059] Referring to Figure 1 The embodiments of the present application disclose a flowchart of a method for calibrating ultra-wideband unevenness and IQ imbalance, which can be implemented by relying on a computer program and can also run on an ultra-wideband unevenness and IQ imbalance calibration device based on the von Neumann system. The computer program can be integrated in an application or run as an independent tool application, and specifically includes:

[0060] S101: Receive a single-tone signal with equal frequency intervals, and perform filtering processing on the single-tone signal to obtain a first single-tone signal. The single-tone signal is a signal corresponding to a sub-band, and multiple sub-bands form an ultra-wideband.

[0061] Specifically, frequency spacing refers to the interval between two adjacent frequencies, represented by df. A single-tone signal refers to a sinusoidal signal of a single frequency; multiple independent sinusoidal signal waveforms superimposed together form a multi-tone signal. The receiver receives single-tone signals transmitted at equal frequency intervals from the test path, which includes a power divider and a directional coupler. The power divider is used to split the energy of one input signal into two or more outputs of equal or unequal energy. In this embodiment, the ultra-wideband signal is divided into multiple sub-bands by the power divider. Each sub-band is flat within itself, but not between different sub-bands. A sub-band, also called a sub-frequency band, is mainly converted from the time domain to the frequency domain using sub-band coding technology, and then divided into several sub-frequency bands. It should be noted that after the ultra-wideband signal is divided into multiple sub-bands, the number of mismatch ripples within each sub-band is greatly reduced. If the number of ripple mismatches within each sub-band can be guaranteed to not exceed 1, that is, each sub-band is relatively flat. Each sub-band corresponds to a single-tone signal transmitted each time.

[0062] In addition, after the receiver receives multiple monotone signals, it filters the monotone signals to remove noise and spurious signals, and finally obtains the first monotone signal.

[0063] For example, within a relatively flat channel of df = 10MHz, a 500MHz ultra-wideband signal only needs to be divided into K = 50 subbands. If the local oscillator of the power divider is set to f0 = 2GHz, then the single-tone signals emitted by the power divider are sequentially (f0)…f2(f0+1df), f3(f0+2df)…fK(f0+(K-1)df), which are (2GHz, 2.01GHz, 2.02GHz…2.5GHz). Therefore, when the receiver's local oscillator is 2.25GHz, the single-tone signals received by the receiver are sequentially (-250MHz, -240MHz,…220MHz, 230MHz, 240MHz, 250MHz), which is essentially the difference frequency signal between the receiver's local oscillator and the emitted single-tone signal.

[0064] like Figure 2 and Figure 3 As shown, Figure 2 This is a schematic diagram of the spectrum after the left 25 single-tone signals received by the receiver are superimposed, as shown below. Figure 3 This is a schematic diagram of the spectrum of the 25 single-tone signals received by the receiver from the right. The 25 single-tone signals on the left are (-250MHz, -240MHz, -230MHz, ... -10MHz), and the 25 single-tone signals on the right are (10MHz, ... 230MHz, 240MHz, 250MHz). As can be seen from the diagram, image interference occurs to the right of the 25 single-tone signals on the left, and similarly, image interference occurs to the left of the 25 single-tone signals on the right. The spectrum refers to the frequency spectral density, which is the frequency distribution curve.

[0065] S102: windowing processing and frequency domain transformation are performed on the first single-tone signal to obtain a second single-tone signal.

[0066] Specifically, the windowing is a process of cutting a signal of a limited length, which is also called time domain truncation. The windowing is equivalent to observing a signal through a time window of a limited length, and is also called windowing. The main function of the windowing is to make the first single-tone signal present a continuous waveform and increase the clarity of the first single-tone signal. The frequency domain transformation is a transformation of a complex time signal or a spatial signal into a structure identified by frequency components. A common way of the frequency domain transformation is fast Fourier transform (FFT). Because some signals are difficult to identify characteristics in the time domain, the characteristics can be easily identified in the frequency domain, and the signal analysis is facilitated.

[0067] After the plurality of first single-tone signals are obtained according to step S101, windowing processing is performed on each first single-tone signal. A feasible windowing processing manner is that a window function in matlab is used to perform windowing processing on the first single-tone signal. The windowing type can be a Hamming window, a Kaiser window, or a triangular window in other embodiments. It should be noted that when the windowing processing and the frequency domain transformation are performed on each first single-tone signal, the I channel signal and the Q channel signal of each first single-tone signal are subjected to the windowing processing and the frequency domain transformation, that is, the I channel signal and the Q channel signal each have a corresponding frequency domain.

[0068] S103: peak points of the second single-tone signal IQ two channels in the corresponding frequency domain are obtained, and IQ phase imbalance information and IQ amplitude imbalance information are determined according to phases and amplitudes of the peak points.

[0069] In an implementable embodiment, the peak points of the second single-tone signal IQ two channels in the corresponding frequency domain are obtained.

[0070] The phase of the corresponding peak point of the Q channel is subtracted from the phase of the corresponding peak point of the I channel and π / 2 to obtain the IQ phase imbalance information. The ratio of the amplitude of the corresponding peak point of the I channel to the amplitude of the corresponding peak point of the Q channel is determined as the IQ amplitude imbalance information.

[0071] Specifically, after the second single-tone signal is determined, a findpeaks function in matlab is used to detect and locate the peak values in the data. In the embodiment of the present application, the findpeaks function is used to select the peak value points of the I-channel signal in the waveform of the corresponding frequency domain and the peak value points of the Q-channel signal in the waveform of the corresponding frequency domain in the second single-tone signal. Since the signal in the frequency domain is a complex number, i.e., each value is a phasor with an amplitude and a phase, the peak value points of the signal can be converted into a complex number, i.e., a complex signal, by Hilbert transform, and the corresponding phase and amplitude can be extracted therefrom after the conversion.

[0072] After the phase and amplitude of the peak value points of the I-channel signal and the phase and amplitude of the peak value points of the Q-channel signal are determined, the phase of the peak value points of the Q-channel signal in each second single-tone signal is subtracted from the phase of the peak value points of the I-channel signal, and then subtracted by π / 2 to obtain a phase difference value. The multiple phase difference values are recorded to determine the IQ phase imbalance information. The amplitude ratio of the amplitude of the peak value points of the I-channel signal to the amplitude of the peak value points of the Q-channel signal in each second single-tone signal is calculated, and the multiple amplitude ratios are recorded to determine the IQ amplitude imbalance information.

[0073] S104: Establish an imbalance signal according to the IQ phase imbalance information and the IQ amplitude imbalance information.

[0074] Specifically, after the IQ phase imbalance information and the IQ amplitude imbalance information are determined, the IQ phase imbalance information is recorded to a preset standard phase imbalance table, the standard phase imbalance table includes standard phase difference values and corresponding frequencies, an IQ two-channel phase error table is obtained, and a phase error function ERR_P(i) is automatically fitted according to the data in the IQ two-channel phase error table by using the nlinfit function in matlab. Similarly, the IQ amplitude imbalance information is recorded to a preset standard amplitude imbalance table, the standard amplitude imbalance table includes standard amplitude ratios and corresponding frequencies, an IQ two-channel amplitude error table is obtained, and an amplitude error function ERR_A(i) is fitted according to the IQ two-channel amplitude error table. Wherein, i is the number of subbands.

[0075] In addition, the function expression of the IQ two-channel standard signal is:

[0076] St_I(i, :) = cos(2*pi*Freq(i) / Fs*t);

[0077] St_Q(i, :) = sin(2*pi*Freq(i) / Fs*t);

[0078] Wherein, pi is the circular constant π, Freq is the frequency point list of the single-tone signal, that is, (-250MHZ, -240M,... 10MHZ, 240MHZ, 250MHZ) in step S01, and Fs is the sampling rate of the single-tone signal.

[0079] Finally, the phase error function ERR_P(i) and the amplitude error function ERR_A(i) are integrated into the expression of the IQ two-way standard signal, so that the unbalanced signal is obtained, which better reflects the unbalanced characteristics of the single-tone signal and avoids the influence of noise and various analog devices on the single-tone signal. The unbalanced signal is specifically expressed as follows:

[0080] St_I(i, :) = 1 / ERR_A(i)*cos(2*pi*Freq(i) / Fs*t-ERR_P(i));

[0081] St_Q(i, :) = 1 / ERR_A(i)*sin(2*pi*Freq(i) / Fs*t-ERR_P(i)).

[0082] S105: Establishing a local standard signal, based on the unbalanced signal and the local standard signal, calculating the coefficients of the equivalent filter by the least square method, and filtering the single-tone signal according to the coefficients of the equivalent filter to calibrate the ultra-wideband unevenness and IQ imbalance.

[0083] In another implementable embodiment, establishing a local standard signal specifically includes:

[0084] Obtaining the sampling rate and the frequency point list of the single-tone signal;

[0085] According to the sampling rate and the frequency point list, the local standard signal is established by matlab;

[0086] Merging a plurality of local standard signals to obtain an original reference signal amplitude and phase vector table, and generating a frequency domain corresponding to each local standard signal according to the original reference signal amplitude and phase vector table.

[0087] Specifically, after the unbalanced signal is established, the local standard signal is generated by the following function,

[0088] St_I(i, :) = cos(2*pi*Freq(i) / Fs*t);

[0089] St_Q(i, :) = sin(2*pi*Freq(i) / Fs*t); Wherein, Fs is the sampling rate, and Freq is the frequency point list. The sampling rate can be 625M, and the frequency point list is [-250MHZ, -240MHZ,... 240MHZ, 250MHZ], which can be input externally.

[0090] Next, the channel corresponding to the local standard signal is used as the reference channel, and the channel corresponding to the unbalanced signal is used as the mismatched channel. Further explanation of the reference and mismatched channels is needed. To ensure consistent frequency response across channels and minimize imbalance in individual tone signals, achieving channel equalization, the channel with the flattest and least distortion within a sub-band (or channel) can be selected as the reference channel, with the remaining channels compared to it. However, in practical applications, determining the optimal sub-band is difficult; instead, a single channel is used as the reference, and equalizers are inserted into the remaining channels to ensure their frequency response matches the reference channel. In this embodiment, the equalizer can be fitted using an Nth-order equivalent filter. Therefore, the coefficients of the equivalent filter need to be determined, which can be understood as the unit impulse response. Frequency response, in electronics, describes the difference in an instrument's processing capability for signals of different frequencies, specifically referring to the gain versus frequency curve. The equivalent filter is a finite-length unit impulse response (FIR) filter.

[0091] Specifically, such as Figure 4 As shown, Figure 4 This diagram illustrates the equalizer calibration principle. The `frestimate` function in MATLAB is used to determine the fixed frequency response of the reference channel corresponding to the local standard signal and the inherent frequency response of the mismatched channel corresponding to the unbalanced signal. Based on the fixed frequency response of the reference channel and the inherent frequency response of the mismatched channel, the frequency response of the i-th equalizer is determined. Finally, the least squares method is used to make the frequency response of the Nth-order equivalent filter approximate the frequency response of the i-th equalizer, thereby determining the optimal coefficients of the equivalent filter. This allows the equivalent filter to effectively replace the equalizer, ensuring that the frequency responses of other channels are consistent with the reference channel, thus completing the joint calibration of ultra-wideband unevenness and IQ imbalance.

[0092] It should be noted that various errors exist in the manufacturing process of each component in the receiver channel. Component aging, orthogonal errors in downconversion, etc., can all cause changes in the amplitude and phase characteristics of the channel. This change is called channel mismatch, and the corresponding channel is called a mismatched channel.

[0093] See Figure 5 This application discloses a flowchart illustrating another calibration method for ultra-wideband unevenness and IQ imbalance, which can be implemented using a computer program or run on a calibration device for ultra-wideband unevenness and IQ imbalance based on the von Neumann architecture. The computer program can be integrated into an application or run as a standalone utility application, specifically including:

[0094] S201: receive a single tone signal with equal frequency intervals, and perform filtering processing on the single tone signal to obtain a first single tone signal, the single tone signal being a signal corresponding to a subband, and a plurality of subbands forming an ultra-wide band.

[0095] Specifically, reference can be made to step S101, which will not be repeated here.

[0096] S202: windowing processing is performed on the real part signal and the imaginary part signal of the first single tone signal respectively to obtain a windowed real part signal and a windowed imaginary part signal.

[0097] S203: the windowed real part signal and the windowed imaginary part signal are subjected to frequency domain transformation respectively to obtain a second single tone signal.

[0098] Specifically, the first single tone signal is determined, and the first single tone signal is represented in a complex form to be converted into a corresponding complex signal, which facilitates further signal processing and analysis. The first single tone signal is a real signal, which is a signal that is real and can be measured in reality. The complex signal is a signal with a function value being a complex number. The complex signal includes a real part signal and an imaginary part signal. The real part is the real part of the amplitude, which represents the amplitude of the signal, while the imaginary part is the imaginary part of the amplitude, which represents the phase of the signal with respect to the frequency. For example, the complex signal C=a+jb, j is the imaginary unit, a is the real part, and b is the imaginary part.

[0099] After the real part signal and the imaginary part signal are determined, the real part signal and the imaginary part signal of the first single tone signal are subjected to windowing processing respectively to obtain a windowed real part signal and a windowed imaginary part signal. The windowing processing process can be specifically referred to step S102. Finally, the windowed real part signal and the windowed imaginary part signal are subjected to frequency domain transformation respectively to obtain a second single tone signal. The frequency domain transformation can be specifically referred to step S102, which will not be repeated here.

[0100] S204: obtain the peak points of the second single tone signal IQ two-way in the corresponding frequency domain, and determine the IQ phase imbalance information and the IQ amplitude imbalance information according to the phases and amplitudes of the peak points.

[0101] Specifically, reference can be made to step S103, which will not be repeated here.

[0102] S205: record the IQ phase imbalance information to a preset ERR_P vector list to obtain a phase error table, and record the IQ amplitude imbalance information to a preset ERR_A vector list to obtain an amplitude error table.

[0103] S206: obtain an imbalance signal according to the phase error table and the amplitude error table.

[0104] S207: combine a plurality of imbalance signals to obtain an error signal amplitude and phase vector table, and generate a frequency domain corresponding to the imbalance signal according to the error signal amplitude and phase vector table.

[0105] Specifically, after the IQ phase imbalance information and the IQ amplitude imbalance information are determined, the IQ phase imbalance information is recorded to a preset ERR_P vector list, wherein the ERR_P vector list comprises a standard variation curve of the IQ phase difference value and the frequency, i.e., standard IQ phase imbalance information, and finally a phase error table is obtained to record error variation of the IQ phase difference value. For details, refer to Figure 6 , Figure 6 The phase error table provided by the embodiment of the present application has a horizontal coordinate of a frequency of a single-tone signal received by a receiver, in units of MHZ, and a vertical coordinate of a phase difference value, in units of degrees.

[0106] Similarly, the IQ amplitude imbalance information is also recorded to a preset ERR_A vector list, wherein the ERR_A vector list comprises a standard variation curve of the IQ amplitude ratio and the frequency, i.e., standard IQ amplitude imbalance information, and finally an amplitude error table is obtained to record error variation of the IQ amplitude ratio. For details, refer to Figure 7 , Figure 7 The amplitude error table provided by the embodiment of the present application has a horizontal coordinate of a frequency of a single-tone signal received by a receiver and a vertical coordinate of an amplitude.

[0107] After the phase error table and the amplitude error table are determined, the imbalance signal is determined according to the phase error table and the amplitude error table, for details, refer to step S104, which will not be repeated here. Then, a plurality of groups of imbalance signals are combined, the imbalance signal comprises an IQ two-way signal St_I(:,i) and St_Q(:,i), and an error signal amplitude-phase vector table is obtained, specifically as follows:

[0108] for i=1:Len

[0109] Err_real(i)=sum(St_I(:,i));

[0110] Err_imag(i)=sum(St_Q(:,i));

[0111] End

[0112] Wherein, Len is a sampling point number, after the error signal amplitude-phase vector table is determined, the frequency domain corresponding to the imbalance signal is further determined, the error signal amplitude-phase vector table is input into a variable function, and the specific implementation is as follows:

[0113] ERR=complex(Err_real,Err_imag);

[0114] C=fft(ERR);

[0115] Wherein, C is a frequency domain corresponding to the imbalance signal.

[0116] S208: Establish a local standard signal, calculate the coefficients of the equivalent filter based on the unbalanced signal and the local standard signal through the least square method, and filter the single tone signal according to the coefficients of the equivalent filter to calibrate the ultra-wideband unevenness and IQ imbalance.

[0117] Specifically, reference can be made to step S105, which will not be repeated here.

[0118] Referring to Figure 8 The embodiment of the present application discloses another flowchart of a method for calibrating ultra-wideband unevenness and IQ imbalance, which can be implemented by a computer program and can also run on an ultra-wideband unevenness and IQ imbalance calibration device based on the von Neumann system. The computer program can be integrated in an application or run as an independent tool application, and specifically includes:

[0119] S301: Receive a single tone signal with equal frequency intervals and filter the single tone signal to obtain a first single tone signal. The single tone signal is a signal corresponding to a subband, and a plurality of subbands form an ultra-wideband.

[0120] S302: Perform windowing processing and frequency domain transformation on the first single tone signal to obtain a second single tone signal.

[0121] S303: Obtain peak points of the IQ two paths of the second single tone signal in the corresponding frequency domain, and determine IQ phase imbalance information and IQ amplitude imbalance information according to the phases and amplitudes of the peak points.

[0122] S304: Establish an unbalanced signal according to the IQ phase imbalance information and the IQ amplitude imbalance information.

[0123] Specifically, reference can be made to steps S101-S104, which will not be repeated here.

[0124] S305: Establish a local standard signal, determine a reference channel corresponding to the local standard signal, and determine a mismatch channel corresponding to the unbalanced signal.

[0125] S306: Determine the frequency response H ref of the all-pass linear phase shift network, and calculate the frequency response of the i-th equalizer to be solved through the first formula H i (jω) = C ref (jω) * H ref / C i (jω), wherein C ref (jω) is the inherent frequency response of the reference channel, C i (jω) is the inherent frequency response of any mismatch channel, i is the number of subbands, H i (jω) is the frequency response of the i-th equalizer to be solved.

[0126] Specifically, after establishing the local standard signal, the reference channel and mismatch channel are determined sequentially based on the local standard signal and the unbalanced signal, using the following formula: H ref =e -j(N-1)Δ / 2 Determine the frequency response of the all-pass linear phase-shift network;

[0127] Among them, H ref The frequency response of the all-pass linear phase-shift network is given, which ensures that all channels have the same tap delay. j is the imaginary unit, N is the order of the equivalent filter in the receiver that functions identically to the equalizer, and Δ is the tap delay. It should be noted that simulations show that H... ref Setting it to the value 1 is also feasible.

[0128] Frequency response H of all-pass linear phase-shift network ref Once determined, according to the first formula: H i (jω)=C ref (jω)*H ref / C i (jω) is used to calculate the frequency response of the i-th equalizer. The fixed frequency response C of the reference channel corresponding to the local standard signal is determined using the frestimate function in MATLAB. ref (jω), the inherent frequency response C of any mismatched channel corresponding to the unbalanced signal. i (jω).

[0129] S307: Obtain the frequency factor matrix, and perform conjugate transpose on the frequency factor matrix to obtain the weighting matrix. The weighting matrix includes the fitting error weights for each frequency point in the frequency factor matrix.

[0130] Specifically, using the formula: A = M * N;

[0131] A represents the size of the frequency factor matrix corresponding to the single-tone signal, and M represents the number of sampling points. After calculating the frequency factor matrix, its conjugate transpose is used to obtain A_H. Finally, A_H is determined as the weighting matrix. The purpose of the weighting matrix is ​​to weight the fitting error at each frequency point in the frequency factor matrix, so that the proportion of the fitting error at different frequency points in the total error is different. Smaller values ​​are applied to less important frequency points to relax the fitting accuracy at these frequency points, while larger values ​​are applied to important frequency points to strengthen the fitting accuracy, thus enabling the equivalent filter to have a better fit in the important frequency bands. The frequency factor matrix is ​​represented as follows:

[0132] S308: Set E(jω), and approximate H using the least squares method. i (jω), and based on the frequency factor matrix, H i(jω) and the weighting matrix are calculated by a preset second formula to obtain the coefficients of the equivalent filter, where E(jω) is the frequency response of the equivalent filter.

[0133] Specifically, after the frequency factor matrix and the weighting matrix are determined, the frequency response E(jω) of the equivalent filter is set as aT(ω)h, where a T is the transpose of the phase shift vector, and h is the weight coefficient vector of the equivalent filter, that is, the coefficients of the equivalent filter. Considering that there is measurement noise in the single-tone signal, the least square method is used to approximate E(jω) to H i (jω). Then, the coefficients of the equivalent filter are calculated by a second formula, and the second formula is as follows:

[0134] FIR_h = inv(W*A)*(W*b);

[0135] where FIR_h is the coefficient of the equivalent filter, b is the value of H i (jω), which has been obtained, W is the weighting matrix, A is the frequency factor matrix, and inv is the matrix inversion. The three parameters b, W and A have been determined, and are directly substituted into the second formula. Finally, the coefficients of the equivalent filter are set as FIR_h, so that all the other channels of the receiver are consistent with the frequency response of the reference channel, thereby effectively calibrating the ultra-wideband unevenness and IQ imbalance of the receiver.

[0136] In addition, as Figure 9 shown is a schematic diagram of the real part and the imaginary part of the FIR corresponding to the coefficients of the equivalent filter. After the received single-tone signals are processed based on the coefficients of the equivalent filter, as Figure 10 and Figure 11 shown, Figure 10 is a schematic diagram of the frequency spectrum of the left 25 single-tone signals after being filtered by the equivalent filter, Figure 11 is a schematic diagram of the frequency spectrum of the right 25 single-tone signals after being filtered by the equivalent filter. As shown in the figure, the mirror image of the left 25 single-tone signals on the right is well suppressed, and the mirror image of the right 25 single-tone signals on the left is well suppressed.

[0137] The implementation principle of the method for calibrating the ultra-wideband unevenness and IQ imbalance in the embodiment of the application is as follows: after receiving a plurality of single-tone signals, the first single-tone signal is obtained through filtering processing, so that the noise other than the single-tone signal is filtered out, and the first single-tone signal is subjected to windowing processing and frequency domain conversion to obtain the second single-tone signal, so that the first single-tone signal has higher definition and is easier to show the characteristics, thereby facilitating subsequent signal analysis. Then, the phase and amplitude of the peak points of the second single-tone signal in the corresponding frequency domain are compared to obtain the IQ phase imbalance information and the IQ amplitude imbalance information, that is, the difference information of the IQ phase and the difference information of the IQ amplitude. Then, the errors of the IQ phase and the IQ amplitude are determined according to the IQ phase imbalance information and the IQ amplitude imbalance information, and then the imbalance signal is determined according to the errors of the IQ phase and the IQ amplitude. Finally, the coefficients of the equivalent filter are calculated by the least square method, so that the frequency response of the channel corresponding to the imbalance signal is equalized to be consistent with the local standard signal under the action of the equivalent filter, thereby suppressing the image interference generated by the single-tone signal, and then effectively calibrating the ultra-wideband unevenness and IQ imbalance.

[0138] The following is an apparatus embodiment of the application, which can be used to execute the method embodiments of the application. For details not disclosed in the apparatus embodiments of the application, refer to the method embodiments of the application.

[0139] Please refer to Figure 12 The structure schematic diagram of the calibration apparatus for the ultra-wideband unevenness and IQ imbalance provided by the embodiment of the application. The calibration apparatus for the ultra-wideband unevenness and IQ imbalance can be realized as all or part of the apparatus by software, hardware or a combination of both. The apparatus 1 comprises a single-tone acquisition module 11, a single-tone processing module 12, an information determination module 13, a signal establishment module 14 and a coefficient determination module 15.

[0140] The single-tone acquisition module 11 is configured to receive the single-tone signals at equal frequency intervals and filter the single-tone signals to obtain the first single-tone signal. The single-tone signal is the signal corresponding to the sub-band, and a plurality of sub-bands constitute the ultra-wideband.

[0141] The single-tone processing module 12 is configured to perform windowing processing and frequency domain conversion on the first single-tone signal to obtain the second single-tone signal.

[0142] The information determination module 13 is configured to obtain the peak points of the second single-tone signal in the corresponding frequency domain, and determine the IQ phase imbalance information and the IQ amplitude imbalance information according to the phase and amplitude of each peak point.

[0143] The signal establishing module 14 is configured to establish an imbalance signal according to the IQ phase imbalance information and the IQ amplitude imbalance information; the coefficient determining module 15 is configured to establish a local standard signal, calculate the coefficient of an equivalent filter based on the imbalance signal and the local standard signal by using a least square method, and filter the single-tone signal according to the coefficient of the equivalent filter, so as to calibrate the ultra-wideband unevenness and the IQ imbalance.

[0144] Optionally, the single-tone processing module 12 is specifically configured to:

[0145] perform windowing processing on the real part signal and the imaginary part signal of the first single-tone signal respectively to obtain a windowed real part signal and a windowed imaginary part signal;

[0146] perform frequency domain transformation on the windowed real part signal and the windowed imaginary part signal respectively to obtain a second single-tone signal.

[0147] Optionally, the information determining module 13 is specifically configured to:

[0148] obtain the peak points of the IQ two paths of the second single-tone signal in the corresponding frequency domain;

[0149] subtract the phase of the peak point of the I path from the phase of the peak point of the Q path and subtract π / 2 to obtain the IQ phase imbalance information; and determine the ratio of the amplitude of the peak point of the I path to the amplitude of the peak point of the Q path as the IQ amplitude imbalance information.

[0150] Optionally, the signal establishing module 14 is specifically configured to:

[0151] record the IQ phase imbalance information into a preset ERR_P vector list to obtain a phase error table, and record the IQ amplitude imbalance information into a preset ERR_A vector list to obtain an amplitude error table, the ERR_P vector list including standard IQ phase imbalance information, and the ERR_A vector list including standard IQ amplitude imbalance information;

[0152] obtain the imbalance signal according to the phase error table and the amplitude error table;

[0153] merge a plurality of groups of imbalance signals to obtain an error signal amplitude and phase vector table, and generate a frequency domain corresponding to the imbalance signal according to the error signal amplitude and phase vector table.

[0154] Optionally, the coefficient determining module 15 is specifically configured to:

[0155] obtain a sampling rate and a frequency point list of the single-tone signal;

[0156] establish the local standard signal by using matlab according to the sampling rate and the frequency point list;

[0157] Multiple sets of local standard signals are merged to obtain the original reference signal amplitude and phase vector table, and the frequency domain corresponding to each local standard signal is generated based on the original reference signal amplitude and phase vector table.

[0158] Optionally, the coefficient determination module 15 is also used for:

[0159] The channel corresponding to the local standard signal is determined as the reference channel, and the channel corresponding to the unbalanced signal is determined as the mismatch channel; the frequency response H of the all-pass linear phase shift network is determined. ref And through the first formula H i (jω)=C ref (jω)*H ref / C i (jω) Calculate the frequency response of the i-th equalizer to be determined, where C ref (jω) is the natural frequency response of the reference channel, C i (jω) represents the intrinsic frequency response of any mismatched channel, i is the number of sub-bands, and H... i (jω) is the frequency response of the i-th equalizer to be determined; obtain the frequency factor matrix, and perform conjugate transpose on the frequency factor matrix to obtain the weighting matrix, which includes the fitting error weight of each frequency point in the frequency factor matrix.

[0160] Let E(jω), and approximate H using the least squares method. i (jω), and based on the frequency factor matrix, H i The coefficients of the equivalent filter are calculated using a pre-defined second formula based on (jω) and the weighting matrix, where E(jω) is the frequency response of the equivalent filter.

[0161] It should be noted that the calibration device for ultra-wideband unevenness and IQ imbalance provided in the above embodiments is only illustrated by the division of the functional modules described above when performing the calibration method for ultra-wideband unevenness and IQ imbalance. In practical applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the calibration device for ultra-wideband unevenness and IQ imbalance provided in the above embodiments and the calibration method embodiment for ultra-wideband unevenness and IQ imbalance belong to the same concept, and the implementation process is detailed in the method embodiment, which will not be repeated here.

[0162] This application also discloses a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, it employs a calibration method for ultra-wideband unevenness and IQ imbalance as described in the above embodiments.

[0163] The computer program can be stored in a computer-readable medium. The computer program includes computer program code, which can be in the form of source code, object code, executable file, or certain middleware. The computer-readable medium includes any entity or device capable of carrying computer program code, recording media, USB flash drive, portable hard drive, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc. It should be noted that the computer-readable medium includes, but is not limited to, the above-mentioned components.

[0164] The above-described ultrawideband unevenness and IQ imbalance calibration method is stored in the computer-readable storage medium and loaded and executed on the processor to facilitate the storage and application of the above method.

[0165] This application also discloses an electronic device in which a computer program is stored in a computer-readable storage medium. When the computer program is loaded and executed by a processor, it employs the aforementioned calibration method for ultra-wideband unevenness and IQ imbalance.

[0166] The electronic device can be a desktop computer, a laptop computer, or a cloud server, and the electronic device includes, but is not limited to, a processor and a memory. For example, the electronic device may also include input / output devices, network access devices, and buses.

[0167] The processor can be a central processing unit (CPU). Of course, depending on the actual use, it can also be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), off-the-shelf programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or any conventional processor, etc., and this application does not limit it in this regard.

[0168] The memory can be an internal storage unit of an electronic device, such as a hard disk or RAM, or an external storage device, such as a plug-in hard disk, smart memory card (SMC), secure digital card (SD), or flash memory card (FC) equipped on the electronic device. Furthermore, the memory can be a combination of an internal storage unit and an external storage device of the electronic device. The memory is used to store computer programs and other programs and data required by the electronic device. The memory can also be used to temporarily store data that has been output or will be output. This application does not limit this.

[0169] In this electronic device, a calibration method for ultra-wideband unevenness and IQ imbalance according to the above embodiment is stored in the memory of the electronic device and loaded and executed on the processor of the electronic device for convenient use.

[0170] The foregoing description is merely an exemplary embodiment of this disclosure and should not be construed as limiting the scope of this disclosure. Any equivalent changes and modifications made in accordance with the teachings of this disclosure shall still fall within the scope of this disclosure. Other embodiments of this disclosure will be readily apparent to those skilled in the art upon consideration of the specification and practice of the disclosure herein. This application is intended to cover any variations, uses, or adaptations of this disclosure that follow the general principles of this disclosure and include common knowledge or customary techniques in the art not described herein. The specification and embodiments are to be considered exemplary only, and the scope and spirit of this disclosure are defined by the claims.

Claims

1. A method of calibration of ultra-wideband (UWB) unevenness and IQ imbalance, the method comprising: The method comprises: Equal frequency interval receives a single tone signal, and filters the single tone signal to obtain a first single tone signal, the single tone signal being a signal corresponding to a subband, and a plurality of the subbands forming an ultra-wide band; Window processing and frequency domain transformation are performed on the first single tone signal to obtain a second single tone signal; IQ phase imbalance information and IQ amplitude imbalance information are determined according to phases and amplitudes of peak points of the second single tone signal in a corresponding frequency domain; An imbalance signal is established according to the IQ phase imbalance information and the IQ amplitude imbalance information, including recording the IQ phase imbalance information into a preset ERR_P vector list to obtain a phase error table, and recording the IQ amplitude imbalance information into a preset ERR_A vector list to obtain an amplitude error table, the ERR_P vector list including standard IQ phase imbalance information, and the ERR_A vector list including standard IQ amplitude imbalance information; The imbalance signal is obtained according to the phase error table and the amplitude error table, including fitting a phase error function according to the phase error table, fitting an amplitude error function according to the amplitude error table, and incorporating the phase error function and the amplitude error function into an expression of a standard IQ two-way signal to obtain the imbalance signal; A plurality of groups of the imbalance signal are combined to obtain an error signal amplitude and phase vector table, and a frequency domain corresponding to the imbalance signal is generated according to the error signal amplitude and phase vector table; A local standard signal is established, a coefficient of an equivalent filter is calculated by a least square method based on the imbalance signal and the local standard signal, and the single tone signal is filtered according to the coefficient of the equivalent filter to calibrate ultra-wide band unevenness and IQ imbalance.

2. The method of calibrating for ultra-wideband unflatness and IQ imbalance according to claim 1, wherein, The window processing and frequency domain transformation on the first single tone signal to obtain the second single tone signal specifically comprises: Window processing is performed on real part signals and imaginary part signals of the first single tone signal to obtain windowed real part signals and windowed imaginary part signals; The windowed real part signals and the windowed imaginary part signals are subjected to frequency domain transformation to obtain the second single tone signal.

3. The method of calibrating for ultra-wideband unflatness and IQ imbalance according to claim 1, wherein, The peak points of the second single tone signal in the corresponding frequency domain are obtained, and the IQ phase imbalance information and the IQ amplitude imbalance information are determined according to phases and amplitudes of the peak points, specifically comprising: The peak points of the second single tone signal in the corresponding frequency domain are obtained; The IQ phase imbalance information is obtained by subtracting a phase of an I-path corresponding peak point from a phase of a Q-path corresponding peak point and subtracting π / 2; The IQ amplitude imbalance information is determined by a ratio of an amplitude of the I-path corresponding peak point to an amplitude of the Q-path corresponding peak point.

4. The method of calibrating for ultra-wideband unflatness and IQ imbalance according to claim 2, wherein, The local standard signal is established, specifically comprising: A sampling rate and a frequency point list of the single tone signal are obtained; The local standard signal is established by matlab according to the sampling rate and the frequency point list; A plurality of groups of the local standard signal are combined to obtain an original reference signal amplitude and phase vector table, and a frequency domain corresponding to each of the local standard signals is generated according to the original reference signal amplitude and phase vector table.

5. The method of calibrating for ultra-wideband unflatness and IQ imbalance according to claim 1, wherein, The method comprises the following steps: calculating the coefficients of an equivalent filter by least square method based on the unbalanced signal and the local standard signal, specifically comprising: determining the channel corresponding to the local standard signal as a reference channel and determining the channel corresponding to the unbalanced signal as a mismatch channel; Determining the frequency response H of an all-pass linear phase shift network ref and by the first equation H i (jω) = C ref (jω) * H ref / C i (jω) is calculated, wherein C ref (jω) is the intrinsic frequency response of the reference channel, C i (jω) is the intrinsic frequency response of any of the mismatch channels, i is the number of subbands, H i (jω) is the frequency response of the i-th equalizer to be determined obtaining a frequency factor matrix, and obtaining a weighting matrix by conjugate transposition of the frequency factor matrix, wherein the weighting matrix comprises the fitting error weight of each frequency point in the frequency factor matrix; Setting E(jω), approximating E(jω) to H i (jω) according to the frequency factor matrix H i (jω) and the weighting matrix are calculated by a preset second formula to obtain the coefficient of the equivalent filter, wherein E(jω) is the frequency response of the equivalent filter.

6. The method of calibrating for ultra-wideband unflatness and IQ imbalance according to claim 5, wherein, The second formula is: FIR_h = inv(W*A)*(W*b), wherein b is H i (jω), FIR_h is the coefficient of the equivalent filter, W is a weighting matrix, and A is a frequency factor matrix.

7. A calibration device for ultra-wideband unevenness and IQ imbalance, characterized in that, The method for calibrating ultra-wideband unevenness and IQ imbalance comprises the following steps: a single tone acquisition module (11) is configured to receive a single tone signal at equal frequency intervals and filter the single tone signal to obtain a first single tone signal, wherein the single tone signal is a signal corresponding to a sub-band, and a plurality of sub-bands form an ultra-wideband; a single tone processing module (12) is configured to perform windowing processing and frequency domain transformation on the first single tone signal to obtain a second single tone signal; an information determination module (13) is configured to obtain the peak points of the IQ two-way of the second single tone signal in the corresponding frequency domain, and determine IQ phase imbalance information and IQ amplitude imbalance information according to the phase and amplitude of each peak point; a signal establishment module (14) is configured to establish an unbalanced signal according to the IQ phase imbalance information and the IQ amplitude imbalance information; 8. A computer-readable storage medium having stored therein a computer program, characterized in that, a coefficient determination module (15) is configured to establish a local standard signal, calculate the coefficients of an equivalent filter by least square method based on the unbalanced signal and the local standard signal, and filter the single tone signal according to the coefficients of the equivalent filter to calibrate the ultra-wideband unevenness and IQ imbalance.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, characterized in that, The computer program is loaded and executed by the processor, and the method of any one of claims 1-6 is adopted. The processor loads and executes the computer program, and the method of any one of claims 1-6 is adopted.

Images