CN112083227B - I/Q unbalance phase error compensation method based on FPA - Google Patents

I/Q unbalance phase error compensation method based on FPA Download PDF

Info

Publication number
CN112083227B
CN112083227B CN202010877152.5A CN202010877152A CN112083227B CN 112083227 B CN112083227 B CN 112083227B CN 202010877152 A CN202010877152 A CN 202010877152A CN 112083227 B CN112083227 B CN 112083227B
Authority
CN
China
Prior art keywords
phase error
frequency
fpa
phase
compensation module
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.)
Expired - Fee Related
Application number
CN202010877152.5A
Other languages
Chinese (zh)
Other versions
CN112083227A (en
Inventor
孟婕
曾浩
郭连平
叶芃
赵禹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202010877152.5A priority Critical patent/CN112083227B/en
Publication of CN112083227A publication Critical patent/CN112083227A/en
Application granted granted Critical
Publication of CN112083227B publication Critical patent/CN112083227B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R25/00Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • G01R23/165Spectrum analysis; Fourier analysis using filters

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Tests Of Electronic Circuits (AREA)

Abstract

The invention discloses an I/Q unbalance phase error compensation method based on an FPA (field programmable gate array), which comprises the steps of firstly obtaining an I/Q unbalance phase error theta (omega) to be compensated through inputting a continuous single tone frequency sweep signal, then designing an I/Q unbalance phase compensation module by using the FPA based on the unbalance phase error theta (omega), equating the I/Q unbalance phase compensation module to be a nonlinear optimization problem, calculating the solution of the nonlinear optimization problem by using the FPA to obtain the optimal parameter of the compensation module, and finally performing I/Q phase error compensation by using the designed I/Q unbalance phase compensation module.

Description

I/Q unbalance phase error compensation method based on FPA
Technical Field
The invention belongs to the technical field of error correction, and particularly relates to an I/Q unbalance phase error compensation method based on an FPA (field programmable Gate array).
Background
The choice of front-end architecture greatly affects receiver integratability and flexibility, with zero-if architectures being favored over the past decade for their advantages of ease of monolithic integration, simplicity of construction, and low cost and power consumption. However, this architecture employs in-phase and quadrature (I/Q) mixing, which is particularly sensitive to I/Q imbalance. Ideally, the I and Q signals should have an exact 90 ° phase difference and equal amplitude. However, in practical situations, errors in the used devices, circuit design, and PCB layout, and differences in the frequency responses of the I-path and Q-path low-pass filters all cause amplitude and phase mismatch I/Q imbalance of the I-path and Q-path signals, which results in incomplete attenuation of the image signal and increases the error rate of baseband signal processing. The effects of I/Q imbalance become more pronounced with the use of higher order modulated waveforms or large bandwidth multichannel signals, which must be improved using additional analog or digital signal processing methods.
The phase unbalance of I/Q unbalance is a key problem, and compared with amplitude errors, the calculation and compensation difficulty of the phase errors is higher; the Flower Pollination Algorithm (FPA) is a natural heuristic that simulates the principle of natural pollination to represent the solution to the problem as pollen of the flower and the process of iterating new solutions as the pollination process between flowers. In recent years, FPA receives wide attention due to the effectiveness of solving practical problems, and compared with a traditional optimization algorithm based on gradient, the FPA is not easy to fall into local optimization and can better find an approximate optimal solution closer to global optimization. By setting proper individual number and iteration number, the solution obtained by FPA can well solve the problem of I/Q unbalance phase error compensation.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an I/Q unbalance phase error compensation method based on an FPA (field programmable gate array), wherein a compensation module of the I/Q unbalance phase error is designed to be equivalent to a nonlinear optimization problem, the nonlinear optimization problem is calculated by the FPA, the optimal parameters of the compensation module are obtained, and then phase compensation is completed.
In order to achieve the above object, the present invention provides an I/Q imbalance phase error compensation method based on FPA, comprising the following steps:
(1) acquiring an I/Q imbalance phase error to be compensated by inputting a single-tone frequency sweeping signal;
(1.1) carrying out equal-interval frequency sweeping on the bandwidth to be tested by using a signal generator, so that single-tone frequency sweeping signals with N frequency values are output at equal intervals within the frequency range of the bandwidth to be tested and are sequentially input to an I/Q receiver;
(1.2) the monophonic frequency sweep signal ri(t) has phase error in the two paths with I/Q, respectively
Figure BDA0002652940480000021
Local oscillator signal cos (ω)LOt) and
Figure BDA0002652940480000022
multiplication of ω whereinLOIs the local oscillator signal frequency;
then, obtaining an unbalanced I/Q test signal after low-pass filtering processing;
xI(t)=cos(ωit)
xQ(t)=sin(ωit+θ(ωi))
where i is 1,2, … N, θ (ω)i) Representing frequency omegaiThe phase error is detected, and the phase error is detected,
Figure BDA0002652940480000023
φ(ωi) Representing a frequency dependent phase error;
(1.3) sampling the unbalanced I/Q test signal by an ADC to obtain an unbalanced I/Q test signal x after samplingI(n) and xQ(n);
Then inputting the test signal into a frequency domain conversion processor, performing discrete Fourier transform to complete the conversion of the I/Q test signal from a time domain to a frequency domain, wherein the converted frequency domain I/Q test signal is as follows:
Figure BDA0002652940480000024
Figure BDA0002652940480000025
(1.4) testing the frequency domain I/Q test signal XI(ω)、XQ(ω) inputting to the imbalance phase estimator to estimate the imbalance phase error;
Figure BDA0002652940480000026
(1.5) continuously processing the N single tone frequency sweep signals according to the method of the steps (1.2) - (1.4) to obtain an imbalance phase error theta (omega), wherein omega is [ omega ], [ omega ]12,…,ωN]T
(2) Designing an I/Q unbalance phase compensation module by using an FPA (field programmable gate array) based on the unbalance phase error theta (omega);
(2.1) building an I/Q unbalance phase compensation module, which specifically comprises a delay module, an FD compensation module and an all-pass filter;
(2.2) designing the total phase-frequency response function phi of the I/Q unbalance phase compensation modulec(ω,R);
φc(ω,R)=-δd·ω+φAPF(ω,R)
Wherein R is [ delta ]d,M,ψ],δdFor fractional delay of FD compensation module, M ═ Mn]T=[M1,M2,…,MN/2]T,MnIs the amplitude of the nth second-order junction pole of the all-pass filter, [ psi ═ phi [ [ phi ]n]T=[ψ12,…,ψN/2]T,ψnThe phase angle of the nth second-order junction pole of the all-pass filter; phi is aAPF(ω, R) is the phase frequency response function of the all-pass filter;
(2.3) defining an error matrix E;
E=[e(ω,R)]T=[e(ω1,R),e(ω2,R),…,e(ωi,R),…,e(ωN,R)]T
wherein e (ω)i,R)=φci,R)-θ(ωi);
(2.4) defining a cost function J (R);
Figure BDA0002652940480000031
(2.5) optimizing J (R) by using a flower pollination algorithm FPA to minimize J (R) and satisfy Mn<1,n=1,2,…,N/2;
(3) Compensating the I/Q phase error by utilizing an I/Q unbalance phase compensation module;
testing the unbalanced I/Q signal x after samplingI(n) and xQ(n) input to the I/Q imbalance phase compensation module to output a compensated I/Q signal xIC[n]、xQC[n]。
The invention aims to realize the following steps:
the invention relates to an I/Q unbalance phase error compensation method based on an FPA (field programmable gate array), which comprises the steps of firstly obtaining an I/Q unbalance phase error theta (omega) to be compensated through inputting a continuous single tone frequency sweep signal, then designing an I/Q unbalance phase compensation module by utilizing the FPA based on the unbalance phase error theta (omega), equivalently designing the I/Q unbalance phase compensation module into a nonlinear optimization problem, then calculating the solution of the nonlinear optimization problem by utilizing the FPA to obtain the optimal parameter of the compensation module, and finally compensating the I/Q phase error by utilizing the designed I/Q unbalance phase compensation module.
Meanwhile, the I/Q unbalance phase error compensation method based on the FPA also has the following beneficial effects:
(1) the method comprises the steps of obtaining the I/Q unbalance phase error to be compensated by inputting a single-tone frequency sweeping signal, and equating the phase error in the whole bandwidth to be tested to the phase error of each frequency sweeping frequency point by utilizing the principle that the phase error in a small bandwidth of a frequency sweeping interval can be approximated to be a constant, so that the novel method for testing the I/Q unbalance phase error is provided.
(2) The I/Q phase unbalance compensation module based on the all-pass filter and the fractional delay filter is designed, and a new thought for compensating I/Q phase unbalance is provided for innovatively converting the design problem of the I/Q phase unbalance compensation module into a constrained nonlinear optimization problem in order to ensure the stability of the all-pass filter realized by adopting an IIR filter structure.
(3) Compared with the traditional optimization algorithm based on gradient, the FPA is not easy to fall into local optimization and can better find an approximate optimal solution closer to global optimization, thereby better completing the compensation work of I/Q (input/output) unbalance phase.
Drawings
FIG. 1 is a flow chart of an FPA-based I/Q imbalance phase error compensation method of the present invention;
FIG. 2 is a schematic diagram of I/Q imbalance phase error;
FIG. 3 is a block diagram of an I/Q imbalance phase compensation module;
fig. 4 is a flow chart of FPA optimization.
Detailed Description
The following description of the embodiments of the present invention is provided in order to better understand the present invention for those skilled in the art with reference to the accompanying drawings. It is to be expressly noted that in the following description, a detailed description of known functions and designs will be omitted when it may obscure the subject matter of the present invention.
Examples
FIG. 1 is a flow chart of the I/Q imbalance phase error compensation method based on FPA according to the present invention.
In this embodiment, as shown in fig. 1, the method for compensating an I/Q imbalance phase error based on FPA of the present invention mainly includes three stages S1-S3:
s1, acquiring an I/Q imbalance phase error to be compensated by inputting a single tone frequency sweep signal;
assuming that there is no amplitude imbalance error in the path, the I/Q imbalance model is shown in the dashed box of fig. 2. OmegaLO=2πfLOIn order to be the angular frequency of the local oscillator signal,
Figure BDA0002652940480000051
the phase difference between two local oscillation signals is irrelevant to the frequency. And the frequency-dependent phase error is partly caused by the frequency response differences introduced by the different components between the I/Q paths, denoted as (ω). The low pass filter LPF is used to remove high frequency components and the analog-to-digital converter ADC is used to sample the I/Q signal. Without loss of generality, the total phase error between the I/Q channels is all equivalent to the Q channel, denoted as
Figure BDA0002652940480000052
Where ω is the baseband signal angular frequency.
In order to correctly compensate the phase error, the I/Q phase imbalance error to be compensated needs to be accurately measured, and the embodiment estimates the I/Q phase imbalance error by using a single-tone test signal, and the specific process is as follows:
s1.1, performing equal-interval frequency sweeping on a bandwidth to be tested by using a signal generator, so that single-tone frequency sweeping signals with N frequency values are output at equal intervals within the frequency range of the bandwidth to be tested and are sequentially input to an I/Q receiver;
in this embodiment, the N single tone frequency sweep signals cover the entire test bandwidth, i.e., the entire test bandwidth is divided into N equal narrow frequency bands, where the phase error parameter in each frequency band can be considered as a constant. Therefore, the phase error in the whole test bandwidth can be equivalent to the phase error of the N frequency points.
S1.2, tone frequency sweep signal ri(t) has phase error in the two paths with I/Q, respectively
Figure BDA0002652940480000053
Local oscillator signal cos (ω)LOt) and
Figure BDA0002652940480000054
multiplication of ω whereinLOIs the local oscillator signal frequency;
then, obtaining an unbalanced I/Q test signal after low-pass filtering processing;
xI(t)=cos(ωit)
xQ(t)=sin(ωit+θ(ωi))
where i is 1,2, … N, θ (ω)i) Representing frequency omegaiThe phase error is detected, and the phase error is detected,
Figure BDA0002652940480000055
φ(ωi) A phase error portion, which represents a frequency dependence, is caused by differences in the frequency response of different components between the I/Q paths;
s1.3, sampling the unbalanced I/Q test signal by an ADC (analog to digital converter) to obtain an unbalanced I/Q test signal x after samplingI(n) and xQ(n);
Then inputting the test signal into a frequency domain conversion processor, performing discrete Fourier transform to complete the conversion of the I/Q test signal from a time domain to a frequency domain, wherein the converted frequency domain I/Q test signal is as follows:
Figure BDA0002652940480000056
Figure BDA0002652940480000061
s1.4, testing the frequency domain I/Q test signal XI(ω)、XQ(omega) is input to an imbalance phase estimator to obtain an imbalance phase estimation result;
Figure BDA0002652940480000062
so as to estimate the unbalanced phase error;
Figure BDA0002652940480000063
s1.5, continuously processing the N monophonic frequency sweep signals according to the method of steps S1.2-S1.4 to obtain an imbalance phase error θ (ω), where ω is [ ω ═ ω12,…,ωN]T
S2, designing an I/Q unbalance phase compensation module by using an FPA (field programmable gate array) based on the unbalance phase error theta (omega);
s2.1, building an I/Q unbalance phase compensation module, as shown in fig. 3, specifically comprising a Fractional Delay (FD) compensation module, an all-pass filter and a delay module. Since the FD compensation filter introduces a delay associated with the order of the FD filter on the channel to be compensated, the delay module adds the delay value to the other channel to keep the delay consistent.
In this embodiment, the FD compensation module of the Q channel is implemented by an FD filter, and the order of D (D is an even number) will bring D/2 delay to the Q channel, so in order to make the I/Q channel delay consistent, the delay module is used to delay the I channel by D/2.
S2.2, the implementation structure of the all-pass filter adopts a cascade structure of an Infinite Impulse Response (IIR) filter, can be expressed in a form of multiplication of a plurality of second-order nodes, and the expression in a z domain is shown as the following formula
Figure BDA0002652940480000064
Wherein N is an even number and is the order of the all-pass filter,
Figure BDA0002652940480000065
is the pole of the nth second-order junction, and has amplitude and phase angle of MnAnd psinThe phase frequency response of the all-pass filter can be expressed as:
Figure BDA0002652940480000066
assuming fractional delay delta of FD compensation moduledThus, the total phase-frequency response of the phase compensation module can be expressed as:
φc(ω,R)=-δd·ω+φAPF(ω,R)
wherein R is [ delta ]d,M,ψ],M=[Mn]T=[M1,M2,…,MN/2]T,MnIs the amplitude of the nth second-order junction pole of the all-pass filter, [ psi ═ phi [ [ phi ]n]T=[ψ12,…,ψN/2]T,ψnThe phase angle of the nth second-order junction pole of the all-pass filter; phi is aAPF(ω, R) is the phase frequency response function of the all-pass filter;
s2.3, in order to compensate for the previously obtained theta (omega), phic(ω, R) should be equal to θ (ω), thus defining an error matrix E;
E=[e(ω,R)]T=[e(ω1,R),e(ω2,R),…,e(ωi,R),…,e(ωN,R)]T
wherein e (ω)i,R)=φci,R)-θ(ωi);
S2.4, the design of the I/Q unbalance phase compensation module can be equivalent to a nonlinear optimization problem, so that a cost function J (R) is defined;
Figure BDA0002652940480000071
s2.5, in the embodiment, an all-pass filter in the I/Q imbalance phase compensation module is realized by adopting an IIR filter structure, and the stability of the all-pass filter is guaranteed to be crucial, so that constraint is added in a nonlinear optimization problem to guarantee a stability condition, the stability condition is further converted into a least square problem with constraint, and finally J (R) is optimized by using a pollination algorithm FPA to enable J (R) to be minimum and meet Mn< 1, N ═ 1,2, …, N/2; compared with the traditional optimization algorithm based on gradient, the FPA is not easy to fall into local optimum, and an approximate optimal solution closer to global optimum can be better found, so that I/Q imbalance phase compensation is better completed.
As shown in fig. 4, we describe in detail the specific process of FPA optimizing j (r), specifically:
1) and setting the pollen individual number S of the FPA, wherein the maximum iteration number maximum and the conversion probability p are 0.8. Randomly generating an initial population containing S individuals, taking a cost function J (R) as a fitness function and calculating the fitness value of each individual of the initial population;
2) and starting to circularly perform population propagation: during each circulation, sequentially traversing each individual and generating a new individual corresponding to the next generation; wherein, when the t-th cycle is performed, the mth individual is traversed
Figure BDA0002652940480000072
Then, a number r is randomly generatedm,t∈[0,1]Then, a new individual of the next generation is generated according to the following formula
Figure BDA0002652940480000073
Figure BDA0002652940480000081
Wherein ε is [0,1 ]]The random parameters are uniformly distributed on the surface of the substrate,
Figure BDA0002652940480000082
is the current optimal solution and gamma is a scale factor, set to 0.01. The global pollination step length L follows the levy distribution:
Figure BDA0002652940480000083
where Γ (λ) is the standard gamma function, s is the Levy flight step, s is0Is the step minimum, step s can be given by:
Figure BDA0002652940480000084
Figure BDA0002652940480000085
where U and V satisfy a gaussian distribution and the constant λ is set to 1.5.
Then will be
Figure BDA0002652940480000086
Substituting into fitness function to calculate its fitness value
Figure BDA0002652940480000087
3) After each circulation is finished, the fitness values of S new individuals generated by the circulation are compared, and the new individual with the minimum fitness value is selected as the optimal individual generated by the circulation
Figure BDA0002652940480000088
4) Comparing the optimal individuals generated by the circulation
Figure BDA0002652940480000089
And the last timeCyclically generated optimal individuals
Figure BDA00026529404800000810
The magnitude of the corresponding fitness value if
Figure BDA00026529404800000811
Then use
Figure BDA00026529404800000812
Is replaced by
Figure BDA00026529404800000813
And is used as the current optimal solution; otherwise, the optimal individuals generated by the last cycle are kept
Figure BDA00026529404800000814
Unchanged and used as the optimal individual generated by the current cycle;
5) judging whether the current cycle time t reaches the maximum iteration number maximum, if not, continuing the cycle; if so, selecting the optimal solution recorded after the maximum cycles
Figure BDA00026529404800000815
As the final optimal solution.
S3, I/Q phase error compensation is carried out by utilizing an I/Q unbalance phase compensation module;
testing the unbalanced I/Q signal x after samplingI(n) and xQ(n) input to the I/Q imbalance phase compensation module to output a compensated I/Q signal xIC[n]、xQC[n]。
Although illustrative embodiments of the present invention have been described above to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, and various changes may be made apparent to those skilled in the art as long as they are within the spirit and scope of the present invention as defined and defined by the appended claims, and all matters of the invention which utilize the inventive concepts are protected.

Claims (3)

1. An I/Q imbalance phase error compensation method based on FPA is characterized by comprising the following steps:
(1) acquiring an I/Q imbalance phase error to be compensated by inputting a single-tone frequency sweeping signal;
(1.1) carrying out equal-interval frequency sweeping on the bandwidth to be tested by using a signal generator, so that single-tone frequency sweeping signals with N frequency values are output at equal intervals within the frequency range of the bandwidth to be tested and are sequentially input to an I/Q receiver;
(1.2) the monophonic frequency sweep signal ri(t) has phase error in the two paths with I/Q, respectively
Figure FDA0002652940470000014
Local oscillator signal cos (ω)LOt) and
Figure FDA0002652940470000015
multiplication of ω whereinLOIs the local oscillator signal frequency;
then, obtaining an unbalanced I/Q test signal after low-pass filtering processing;
xI(t)=cos(ωit)
xQ(t)=sin(ωit+θ(ωi))
where i is 1,2, … N, θ (ω)i) Representing frequency omegaiThe phase error is detected, and the phase error is detected,
Figure FDA0002652940470000016
(1.3) sampling the unbalanced I/Q test signal by an ADC to obtain an unbalanced I/Q test signal x after samplingI(n) and xQ(n);
Then inputting the test signal into a frequency domain conversion processor, performing discrete Fourier transform to complete the conversion of the I/Q test signal from a time domain to a frequency domain, wherein the converted frequency domain I/Q test signal is as follows:
Figure FDA0002652940470000011
Figure FDA0002652940470000012
(1.4) testing the frequency domain I/Q test signal XI(ω)、XQ(ω) inputting to the imbalance phase estimator to estimate the imbalance phase error;
Figure FDA0002652940470000013
(1.5) continuously processing the N single tone frequency sweep signals according to the method of the steps (1.2) - (1.4) to obtain an imbalance phase error theta (omega), wherein omega is [ omega ], [ omega ]12,…,ωN]T
(2) Designing an I/Q unbalance phase compensation module by utilizing an FPA algorithm based on the unbalance phase error theta (omega);
(2.1) building an I/Q unbalance phase compensation module, which specifically comprises a delay module, an FD compensation module and an all-pass filter;
(2.2) designing the total phase-frequency response function phi of the I/Q unbalance phase compensation modulec(ω,R);
φc(ω,R)=-δd·ω+φAPF(ω,R)
Wherein R is [ delta ]d,M,ψ],δdFor fractional delay of FD compensation module, M ═ Mn]T=[M1,M2,…,MN/2]T,MnIs the amplitude of the nth second-order junction pole of the all-pass filter, [ psi ═ phi [ [ phi ]n]T=[ψ12,…,ψN/2]T,ψnThe phase angle of the nth second-order junction pole of the all-pass filter; phi is aAPF(ω, R) is the phase frequency response function of the all-pass filter;
(2.3) defining an error matrix E;
E=[e(ω,R)]T=[e(ω1,R),e(ω2,R),…,e(ωi,R),…,e(ωN,R)]T
wherein e (ω)i,R)=φci,R)-θ(ωi);
(2.4) defining a cost function J (R);
Figure FDA0002652940470000021
(2.5) optimizing J (R) by using a flower pollination algorithm FPA to minimize J (R) and satisfy Mn<1,n=1,2,…,N/2;
(3) Compensating the I/Q phase error by utilizing an I/Q unbalance phase compensation module;
testing the unbalanced I/Q signal x after samplingI(n) and xQ(n) input to the I/Q imbalance phase compensation module to output a compensated I/Q signal xIC[n]、xQC[n]。
2. The FPA based I/Q imbalance phase error compensation method of claim 1, wherein the phase-frequency response of the all-pass filter is expressed as:
Figure FDA0002652940470000022
3. the FPA-based I/Q imbalance phase error compensation method according to claim 1, wherein the specific process of optimizing j (r) by using FPA algorithm in step (2.5) is as follows:
1) setting the pollen individual number S, the maximum iteration number Maxiter and the conversion probability p of the FPA algorithm; randomly generating an initial population containing S individuals, taking a cost function J (R) as a fitness function and calculating the fitness value of each individual of the initial population;
2) and starting to circularly perform population propagation: during each circulation, sequentially traversing each individual and generating a new individual corresponding to the next generation; wherein, when the t-th circulation,go to the m-th individual
Figure FDA0002652940470000031
Then, a number r is randomly generatedm,t∈[0,1]Then, a new individual of the next generation is generated according to the following formula
Figure FDA0002652940470000032
Figure FDA0002652940470000033
Wherein ε is [0,1 ]]The random parameters are uniformly distributed on the surface of the substrate,
Figure FDA0002652940470000034
is the current optimal solution, gamma is a scale factor, L is the global pollination step length, m is not equal to j is not equal to k, m, j, k belongs to [1, S ]];
Then will be
Figure FDA0002652940470000035
Substituting into fitness function to calculate its fitness value
Figure FDA0002652940470000036
3) After each circulation is finished, the fitness values of S new individuals generated by the circulation are compared, and the new individual with the minimum fitness value is selected as the optimal individual generated by the circulation
Figure FDA0002652940470000037
4) Comparing the optimal individuals generated by the circulation
Figure FDA0002652940470000038
And the optimal individuals generated in the last circulation
Figure FDA0002652940470000039
The magnitude of the corresponding fitness value if
Figure FDA00026529404700000310
Then use
Figure FDA00026529404700000311
Is replaced by
Figure FDA00026529404700000312
And is used as the current optimal solution; otherwise, the optimal individuals generated by the last cycle are kept
Figure FDA00026529404700000313
Unchanged and used as the optimal individual generated by the current cycle;
5) judging whether the current cycle time t reaches the maximum iteration number maximum, if not, continuing the cycle; if so, selecting the optimal solution recorded after the maximum cycles
Figure FDA00026529404700000314
As the final optimal solution.
CN202010877152.5A 2020-08-27 2020-08-27 I/Q unbalance phase error compensation method based on FPA Expired - Fee Related CN112083227B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010877152.5A CN112083227B (en) 2020-08-27 2020-08-27 I/Q unbalance phase error compensation method based on FPA

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010877152.5A CN112083227B (en) 2020-08-27 2020-08-27 I/Q unbalance phase error compensation method based on FPA

Publications (2)

Publication Number Publication Date
CN112083227A CN112083227A (en) 2020-12-15
CN112083227B true CN112083227B (en) 2021-05-14

Family

ID=73728656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010877152.5A Expired - Fee Related CN112083227B (en) 2020-08-27 2020-08-27 I/Q unbalance phase error compensation method based on FPA

Country Status (1)

Country Link
CN (1) CN112083227B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113132031B (en) * 2021-04-25 2022-07-12 成都天奥测控技术有限公司 IQ correction method for receiving end
CN113259286B (en) * 2021-07-07 2021-10-15 易兆微电子(杭州)股份有限公司 Transmitter and IQ imbalance and direct current offset compensation method and device thereof
CN115276683B (en) * 2022-09-27 2022-12-20 北京力通通信有限公司 IQ imbalance blind estimation tracking correction method for zero intermediate frequency broadband radio frequency transmitter
CN115765016B (en) * 2022-11-30 2024-03-19 国网四川省电力公司电力科学研究院 Distributed photovoltaic output fluctuation suppression method and system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101461201A (en) * 2006-06-06 2009-06-17 高通股份有限公司 Fast in-phase and quadrature imbalance calibration
CN101917726A (en) * 2010-07-15 2010-12-15 新邮通信设备有限公司 Method for reducing out-of-band radiation of user equipment and user equipment
CN102449165A (en) * 2009-04-03 2012-05-09 解码遗传学私营有限责任公司 Genetic markers for risk management of atrial fibrillation and stroke
CN103425051A (en) * 2013-08-21 2013-12-04 北京航空航天大学 Unbalance identification and vibration suppression control system for magnetic suspension rotating machinery
CN104320204A (en) * 2014-11-11 2015-01-28 中国电子科技集团公司第五十四研究所 Local oscillator IQ signal phase and amplitude calibration device
CN108140232A (en) * 2015-06-10 2018-06-08 无比视视觉技术有限公司 For handling the image processor of image and method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7280696B2 (en) * 2002-05-20 2007-10-09 Simmonds Precision Products, Inc. Video detection/verification system
US10401499B2 (en) * 2016-03-16 2019-09-03 Raytheon Company Laser beam projection system with dynamic phase compensation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101461201A (en) * 2006-06-06 2009-06-17 高通股份有限公司 Fast in-phase and quadrature imbalance calibration
CN102449165A (en) * 2009-04-03 2012-05-09 解码遗传学私营有限责任公司 Genetic markers for risk management of atrial fibrillation and stroke
CN101917726A (en) * 2010-07-15 2010-12-15 新邮通信设备有限公司 Method for reducing out-of-band radiation of user equipment and user equipment
CN103425051A (en) * 2013-08-21 2013-12-04 北京航空航天大学 Unbalance identification and vibration suppression control system for magnetic suspension rotating machinery
CN104320204A (en) * 2014-11-11 2015-01-28 中国电子科技集团公司第五十四研究所 Local oscillator IQ signal phase and amplitude calibration device
CN108140232A (en) * 2015-06-10 2018-06-08 无比视视觉技术有限公司 For handling the image processor of image and method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Analysis of Load Allocation Problem Using Flower Pollination Algorithm with Constraints;Deblina Maity等;《2019 16th International Conference on Electrical Engineering/Electronics, Computer, Telecommunications and Information Technology (ECTI-CON)》;20200113;全文 *
OFDM联合信道均衡与I/Q失衡补偿算法;陈楚等;《通信技术》;20080930;全文 *

Also Published As

Publication number Publication date
CN112083227A (en) 2020-12-15

Similar Documents

Publication Publication Date Title
CN112083227B (en) I/Q unbalance phase error compensation method based on FPA
US8218687B2 (en) Frequency dependent I/Q imbalance estimation
CN102870389B (en) For the digital IQ imbalance compensation of the two conversion receivers of two carrier waves
US8842771B2 (en) Amplitude flatness and phase linearity calibration for RF sources
US7088765B1 (en) Vector calibration system
US8224269B2 (en) Vector modulator calibration system
JP6334132B2 (en) Integrated calibration apparatus and runtime calibration method
US20110199096A1 (en) Correction of non-linearities in adcs
US8976914B2 (en) Multi-tap IQ imbalance estimation and correction circuit and method
US9985812B1 (en) Systems and methods for IQ demodulation with error correction
TW200832946A (en) Frequency response correction for a receiver having a frequency translation device
CN107797099A (en) A kind of real-time internal calibration processing method and processing device of Multichannel Digital Receiver
CN109583575B (en) Processing method for improving instrument vector signal analysis performance based on deep learning
Song et al. Time-and frequency-interleaving: Distinctions and connections
Park et al. A calibration for the modulated wideband converter using sinusoids with unknown phases
CN110557122A (en) method for correcting frequency response non-uniformity error of TIADC system
US4809203A (en) Hybrid analog-digital filter
CN109308453A (en) Undersampled signal frequency estimating methods and device based on pattern clustering and spectrum correction
Ivanisevic et al. Impedance spectroscopy systems: Review and an all-digital adaptive IIR filtering approach
CN104065598B (en) Broadband IQ disequilibrium regulatings method, apparatus and system
CN115882970A (en) Method and system for correcting received IQ imbalance
JP6756491B2 (en) Methods and Devices for Generating Normalized Phase Spectrum
CN104040982B (en) The mechanism of correction is detracted for I/Q and detracts measurement using the transmitter for offseting local oscillator
US20050228843A1 (en) Method and apparatus for demodulating square root
CN111614586B (en) Estimation method and device of transmit-end IQ imbalance parameters

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210514