阅读说明：本技术 一种计算并补偿i/q信号线性相位失衡的方法 (A method of it calculates and compensates i/q signal linear phase imbalance ) 是由 孟婕 叶芃 王厚军 郭连平 赵禹 曾浩 田雨 蒋俊 王猛 于 2019-08-01 设计创作，主要内容包括：本发明公开了一种计算并补偿I/Q信号线性相位失衡的方法,在不考虑噪声影响下,对解调后的I/Q两路基带信号进行ADC采样,以及N点傅里叶变换,再取单边谱计算互功率谱及互功率谱的相位谱,然后利用三点法进行相位解缠,最后根据解缠后的群延时差δ和相位偏置差<Image he="53" wi="48" file="DDA0002152130370000011.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>进行相位补偿,还原出时域信号。(The invention discloses a kind of methods for calculating and compensating i/q signal linear phase imbalance, in the case where not considering influence of noise, ADC sampling is carried out to the I/Q two-way baseband signal after demodulation, and N point Fourier transformation, unilateral spectrum is taken to calculate the phase spectrum of crosspower spectrum and crosspower spectrum again, then phase unwrapping is carried out using line-of-sight course, the group delay time difference δ and phase bias after finally twining according to solution are poor Phase compensation is carried out, time-domain signal is restored.)

一种计算并补偿I/Q信号线性相位失衡的方法

技术领域

本发明属于误差校正技术领域，更为具体地讲，涉及一种计算并补偿I/Q信号线性相位失衡的方法。

背景技术

前端架构的选择处于影响接收器可集成性和灵活性的重要位置，而传统的超外差架构有很多缺点，在过去的十几年，零中频结构由于具有易于单片集成、结构简单和成本功耗低等优点而备受青睐，与外差式接收结构不同的是，零中频接收结构中的本振频率等同于或者接近于输入射频信号的频率，从而得到直流中心或低中频信号。然而该结构采用同相和正交(I/Q)正交混频，对射频失真尤为敏感，其中I/Q信号正交失衡是其中的核心影响因素。理想情况下，I路和Q路信号应该具有精确的90°相位差和相等的幅度。然而，在实际情况中，所用器件的不同、电路设计的不同、PCB布局的误差都会使I路和Q路信号的信号产生相位和幅度的差异。此外，I路和Q路的低通滤波器，数据转换器和放大器的频率响应的差异也很大程度的影响到的整体幅度和相位失配。以上因素都会造成I/Q失衡，导致不完全的镜像信号衰减，且在基带处理信号时提高误码率。除镜像频率干扰问题外，I/Q失衡会严重影响基于预失真的功率放大器线性化技术的性能。随着使用更高阶调制波形或更宽带的多信道信号，这些失衡的影响变得更加明显。在许多系统中必须使用额外的模拟或数字信号处理来改善。

I/Q失衡的相位失衡是其中的关键难题，其主要由模拟解调本振的不理想以及I/Q信号两条通路的群延时不一致造成。其中不同的相位偏置一般由本振信号的不理想造成，其导致的相位误差与频率无关，不受基带信号频率影响，可理解为常数，是窄带通信考虑的主要误差，但是对于宽带通信而言，还要考虑群延时偏差。群延时偏差由I、Q两路路径引入的不同延时导致，这部分导致的相位失衡与频率有关，会随着基带信号频率的变化而改变，较难计算及消除，因此本发明提供了一种基于计算互功率谱的同时消除I/Q信号两种相位失衡误差，即总的线性相位失衡的方法。

发明内容

本发明的目的在于克服现有技术的不足，提供计算并补偿I/Q信号线性相位失衡的方法，通过三点法进行相位解缠，进而补偿I/Q信号的线性相位失衡。

为实现上述发明目的，本发明一种计算并补偿I/Q信号线性相位失衡的方法，其特征在于，包括以下步骤：

(1)、在不考虑噪声影响下，对解调后的I/Q两路基带信号进行ADC采样，得到采用后的两路离散信号，记为：x_i[n]和x_q[n]，n表示时间的离散时刻；再对x_i[n]和x_q[n]分别进行N点傅里叶变换，得到X_i[k]和X_q[k]，其中，N为偶数，k满足：k＝-N/2+1,…,N/2；

(2)、对X_q[k]进行希尔伯特变换，得到

分别取X_i[k]和的正频率一侧的单边谱，得到X_i'[k]和X'_q[k]，k＞0且满足：

其中，δ和分别为Q通道相对于I通道的群延时差和相位偏置差；

(3)、计算变换后的I/Q两路信号的互功率谱

其中，是X_q[k]的共轭；

(4)、计算互功率谱的相位谱，得相位差Δθ[k]；

(5)、利用三点法进行相位解缠

(5.1)、在Δθ[k]加上2π的整数倍来进行解缠绕；

其中，是Δθ[k_i]解缠后的相位值，m_i为整数，M为选取的频点个数；

(5.2)、取相位谱上的信噪比SNR最大的三个频点及其相位值，记为(k₁,Δθ[k₁])、(k₂,Δθ[k₂])、(k₃,Δθ[k₃])，并作为未解缠的样本点；

(5.3)、利用前两个样本点计算Δθ[k]的斜率；

(5.4)、计算I/Q两通道的群延时差δ和相位偏置差

设第一个样本点不需要加2π的整数倍，则m₁＝0；

根据在同一条直线上，可得：

将(4)式代入(6)式，化简可得：

(k₃-k₁)m₂＝C+(k₂-k₁)m₃ (7)

其中，变量

将(7)式变形可得：(k₃-k₁)m₂-C必须是(k₂-k₁)的整数倍，因此，将(7)式改写为：

(k₃-k₁)m₂＝C(mod(k₂-k₁)) (8)

其中，mod(·)表示求余；

利用线性模方程求解(8)式，可得：

(k₃-k₁)m₂＝m_2,0+(L(k₂-k₁)/D) (9)

其中，D是(k₃-k₁)和(k₂-k₁)的最大公约数，L＝0,1,…,D-1；m_2,0＝A(C/D)(mod(k₂-k₁))，A为整数；

根据(9)式得到m₂，再将m₂代入到(5)式，得到Q通道相对于I通道的群延时差δ；将δ代入到(3)式，得到Q通道相对于I通道的相位偏置差

(6)、延时模块根据群延时差δ和相位偏置差对单边谱X'_q[k]进行补偿，得到补偿后的单边谱X'_qc[k]；

利用镜像恢复模块得到X'_qc[k]的共轭再将X'_qc[k]和组合成双边谱X_qcd[k]；

将双边谱X_qcd[k]进行傅里叶反变换，得到补偿后的Q路时域信号x_qcd[n]。

本发明的发明目的是这样实现的：

本发明一种计算并补偿I/Q信号线性相位失衡的方法，在不考虑噪声影响下，对解调后的I/Q两路基带信号进行ADC采样，以及N点傅里叶变换，再取单边谱计算互功率谱及互功率谱的相位谱，然后利用三点法进行相位解缠，最后根据解缠后的群延时差δ和相位偏置差进行相位补偿，还原出时域信号。

同时，本发明一种计算并补偿I/Q信号线性相位失衡的方法还具有以下有益效果：

(1)、本发明只需要I/Q信号的三个频点就可以实现估计，不需要导频信号及额外的硬件。

(2)、本发明可以实现群延时差和相位偏置差的同时计算及补偿。

(3)、存在I/Q幅度失衡的情况下，本发明还可以完成相位失衡校正。

附图说明

图1是本发明一种计算并补偿I/Q信号线性相位失衡的方法流程图；

图2是正交接收机I/Q通道不平衡模型框图；

图3是计算线性相位失衡的原理图；

图4是线性相位失衡补偿框图。

具体实施方式

下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。

实施例

图1是本发明一种计算并补偿I/Q信号线性相位失衡的方法流程图。

在本实施例中，如图1所示，本发明一种计算并补偿I/Q信号线性相位失衡的方法，包括以下步骤：

S1、基带信号预处理

在不考虑噪声影响下，如图2所示，接收到的射频模拟信号为r(t)，本振(LO)输出的两路信号分别为cos(ω_ct)和 代表Q路本振信号相对I路的相位偏置(理想情况)，ω_c为载波频率，设基带信号频率为ω，低通滤波器用于滤除混频后的高频成分；

对解调后的I/Q两路基带信号进行ADC采样，得到采用后的两路离散信号，记为：x_i[n]和x_q[n]，n表示时间的离散时刻；再对x_i[n]和x_q[n]分别进行N点傅里叶变换，得到X_i[k]和X_q[k]，其中，N为偶数，k满足：k＝-N/2+1,…,N/2；

S2、计算单边谱信号

对X_q[k]进行希尔伯特变换，得到

分别取X_i[k]和的正频率一侧的单边谱，得到X_i'[k]和X'_q[k]，k＞0，通过上述变换后得到的Q路信号是原Q信号的正交信号，此时Q路信号与I路信号只存在时延误差和相位偏置误差，即满足：

其中，δ和分别为Q通道相对于I通道的群延时差和相位偏置差；

S3、计算变换后的I/Q两路信号的互功率谱

其中，是X_q[k]的共轭；那么，还可以改写为： 为x_i[n]的自功率谱，相位差为该参数包括待估计参数δ和

S4、计算互功率谱的相位谱，得相位差Δθ[k]；

S5、利用三点法进行相位解缠，具体过程如图3所示；

S5.1、由(3)可知，Δθ[k]是一个斜率为的线性相位方程，群延时差δ是该方程的斜率，相位偏置差是该方程在Δθ[k]轴上的截距。由于正切函数的周期性，Δθ[k]在-π到π之间缠绕，因此在计算斜率和截距之前，需要在Δθ[k]加上2π的整数倍来进行解缠绕；

其中，是Δθ[k_i]解缠后的相位值，m_i为整数，M为选取的频点个数；

S5.2、为了估计斜率，至少需要三个相位值，取相位谱上的信噪比SNR最大的三个频点及其相位值，记为(k₁,Δθ[k₁])、(k₂,Δθ[k₂])、(k₃,Δθ[k₃])，并作为未解缠的样本点，利用三点法相位解缠计算群延时差δ和相位偏置差

S5.3、利用前两个样本点计算Δθ[k]的斜率；

S5.4、计算I/Q两通道的群延时差δ和相位偏置差

设第一个样本点不需要加2π的整数倍，则m₁＝0；

由于(k₁,Δθ[k₁])、(k₂,Δθ[k₂])已知，m₂是唯一的待估计变量；根据 在同一条直线上，可得：

将(4)式代入(6)式，化简可得：

(k₃-k₁)m₂＝C+(k₂-k₁)m₃ (7)

其中，变量

将(7)式变形可得：(k₃-k₁)m₂-C必须是(k₂-k₁)的整数倍，因此，将(7)式改写为：

(k₃-k₁)m₂＝C(mod(k₂-k₁)) (8)

其中，mod(·)表示求余；

利用线性模方程求解(8)式，可得：

(k₃-k₁)m₂＝m_2,0+(L(k₂-k₁)/D) (9)

其中，D是(k₃-k₁)和(k₂-k₁)的最大公约数，L＝0,1,…,D-1；m2,0＝A(C/D)(mod(k₂-k₁))，A为整数，由扩展欧几里得算法得到，且与最大公约数D满足如下关系：

D＝A(k₃-k₁)+B(k₂-k₁)

其中，B为整数；

根据(9)式得到m₂，再将m₂代入到(5)式，得到Q通道相对于I通道的群延时差δ；将δ代入到(3)式，得到Q通道相对于I通道的相位偏置差

S6、补偿并还原时域信号

延时模块根据群延时差δ和相位偏置差对单边谱X'_q[k]进行补偿，得到补偿后的单边谱X'_qc[k]；

利用镜像恢复模块得到X'_qc[k]的共轭再将X'_qc[k]和组合成双边谱X_qcd[k]；

将双边谱X_qcd[k]进行傅里叶反变换，得到补偿后的Q路时域信号x_qcd[n]。

尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。