CN110535575A - 一种计算并补偿i/q信号线性相位失衡的方法 - Google Patents
一种计算并补偿i/q信号线性相位失衡的方法 Download PDFInfo
- Publication number
- CN110535575A CN110535575A CN201910705992.0A CN201910705992A CN110535575A CN 110535575 A CN110535575 A CN 110535575A CN 201910705992 A CN201910705992 A CN 201910705992A CN 110535575 A CN110535575 A CN 110535575A
- Authority
- CN
- China
- Prior art keywords
- phase
- spectrum
- signal
- formula
- channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/24—Testing correct operation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/24—Testing correct operation
- H04L1/248—Distortion measuring systems
Abstract
本发明公开了一种计算并补偿I/Q信号线性相位失衡的方法,在不考虑噪声影响下,对解调后的I/Q两路基带信号进行ADC采样,以及N点傅里叶变换,再取单边谱计算互功率谱及互功率谱的相位谱,然后利用三点法进行相位解缠,最后根据解缠后的群延时差δ和相位偏置差进行相位补偿,还原出时域信号。
Description
技术领域
本发明属于误差校正技术领域,更为具体地讲,涉及一种计算并补偿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采样,得到采用后的两路离散信号,记为:xi[n]和xq[n],n表示时间的离散时刻;再对xi[n]和xq[n]分别进行N点傅里叶变换,得到Xi[k]和Xq[k],其中,N为偶数,k满足:k=-N/2+1,…,N/2;
(2)、对Xq[k]进行希尔伯特变换,得到
分别取Xi[k]和的正频率一侧的单边谱,得到Xi'[k]和X'q[k],k>0且满足:
其中,δ和分别为Q通道相对于I通道的群延时差和相位偏置差;
(3)、计算变换后的I/Q两路信号的互功率谱
其中,是Xq[k]的共轭;
(4)、计算互功率谱的相位谱,得相位差Δθ[k];
(5)、利用三点法进行相位解缠
(5.1)、在Δθ[k]加上2π的整数倍来进行解缠绕;
其中,是Δθ[ki]解缠后的相位值,mi为整数,M为选取的频点个数;
(5.2)、取相位谱上的信噪比SNR最大的三个频点及其相位值,记为(k1,Δθ[k1])、(k2,Δθ[k2])、(k3,Δθ[k3]),并作为未解缠的样本点;
(5.3)、利用前两个样本点计算Δθ[k]的斜率;
(5.4)、计算I/Q两通道的群延时差δ和相位偏置差
设第一个样本点不需要加2π的整数倍,则m1=0;
根据在同一条直线上,可得:
将(4)式代入(6)式,化简可得:
(k3-k1)m2=C+(k2-k1)m3 (7)
其中,变量
将(7)式变形可得:(k3-k1)m2-C必须是(k2-k1)的整数倍,因此,将(7)式改写为:
(k3-k1)m2=C(mod(k2-k1)) (8)
其中,mod(·)表示求余;
利用线性模方程求解(8)式,可得:
(k3-k1)m2=m2,0+(L(k2-k1)/D) (9)
其中,D是(k3-k1)和(k2-k1)的最大公约数,L=0,1,…,D-1;m2,0=A(C/D)(mod(k2-k1)),A为整数;
根据(9)式得到m2,再将m2代入到(5)式,得到Q通道相对于I通道的群延时差δ;将δ代入到(3)式,得到Q通道相对于I通道的相位偏置差
(6)、延时模块根据群延时差δ和相位偏置差对单边谱X'q[k]进行补偿,得到补偿后的单边谱X'qc[k];
利用镜像恢复模块得到X'qc[k]的共轭再将X'qc[k]和组合成双边谱Xqcd[k];
将双边谱Xqcd[k]进行傅里叶反变换,得到补偿后的Q路时域信号xqcd[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采样,得到采用后的两路离散信号,记为:xi[n]和xq[n],n表示时间的离散时刻;再对xi[n]和xq[n]分别进行N点傅里叶变换,得到Xi[k]和Xq[k],其中,N为偶数,k满足:k=-N/2+1,…,N/2;
S2、计算单边谱信号
对Xq[k]进行希尔伯特变换,得到
分别取Xi[k]和的正频率一侧的单边谱,得到Xi'[k]和X'q[k],k>0,通过上述变换后得到的Q路信号是原Q信号的正交信号,此时Q路信号与I路信号只存在时延误差和相位偏置误差,即满足:
其中,δ和分别为Q通道相对于I通道的群延时差和相位偏置差;
S3、计算变换后的I/Q两路信号的互功率谱
其中,是Xq[k]的共轭;那么,还可以改写为: 为xi[n]的自功率谱,相位差为该参数包括待估计参数δ和
S4、计算互功率谱的相位谱,得相位差Δθ[k];
S5、利用三点法进行相位解缠,具体过程如图3所示;
S5.1、由(3)可知,Δθ[k]是一个斜率为的线性相位方程,群延时差δ是该方程的斜率,相位偏置差是该方程在Δθ[k]轴上的截距。由于正切函数的周期性,Δθ[k]在-π到π之间缠绕,因此在计算斜率和截距之前,需要在Δθ[k]加上2π的整数倍来进行解缠绕;
其中,是Δθ[ki]解缠后的相位值,mi为整数,M为选取的频点个数;
S5.2、为了估计斜率,至少需要三个相位值,取相位谱上的信噪比SNR最大的三个频点及其相位值,记为(k1,Δθ[k1])、(k2,Δθ[k2])、(k3,Δθ[k3]),并作为未解缠的样本点,利用三点法相位解缠计算群延时差δ和相位偏置差
S5.3、利用前两个样本点计算Δθ[k]的斜率;
S5.4、计算I/Q两通道的群延时差δ和相位偏置差
设第一个样本点不需要加2π的整数倍,则m1=0;
由于(k1,Δθ[k1])、(k2,Δθ[k2])已知,m2是唯一的待估计变量;根据 在同一条直线上,可得:
将(4)式代入(6)式,化简可得:
(k3-k1)m2=C+(k2-k1)m3 (7)
其中,变量
将(7)式变形可得:(k3-k1)m2-C必须是(k2-k1)的整数倍,因此,将(7)式改写为:
(k3-k1)m2=C(mod(k2-k1)) (8)
其中,mod(·)表示求余;
利用线性模方程求解(8)式,可得:
(k3-k1)m2=m2,0+(L(k2-k1)/D) (9)
其中,D是(k3-k1)和(k2-k1)的最大公约数,L=0,1,…,D-1;m2,0=A(C/D)(mod(k2-k1)),A为整数,由扩展欧几里得算法得到,且与最大公约数D满足如下关系:
D=A(k3-k1)+B(k2-k1)
其中,B为整数;
根据(9)式得到m2,再将m2代入到(5)式,得到Q通道相对于I通道的群延时差δ;将δ代入到(3)式,得到Q通道相对于I通道的相位偏置差
S6、补偿并还原时域信号
延时模块根据群延时差δ和相位偏置差对单边谱X'q[k]进行补偿,得到补偿后的单边谱X'qc[k];
利用镜像恢复模块得到X'qc[k]的共轭再将X'qc[k]和组合成双边谱Xqcd[k];
将双边谱Xqcd[k]进行傅里叶反变换,得到补偿后的Q路时域信号xqcd[n]。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种计算并补偿I/Q信号线性相位失衡的方法,其特征在于,包括以下步骤:
(1)、在不考虑噪声影响下,对解调后的I/Q两路基带信号进行ADC采样,得到采用后的两路离散信号,记为:xi[n]和xq[n],n表示时间的离散时刻;再对xi[n]和xq[n]分别进行N点傅里叶变换,得到Xi[k]和Xq[k],其中,N为偶数,k满足:k=-N/2+1,...,N/2;
(2)、对Xq[k]进行希尔伯特变换,得到
分别取Xi[k]和的正频率一侧的单边谱,得到X′i[k]和X′q[k],k>0且满足:
其中,δ和分别为Q通道相对于I通道的群延时差和相位偏置差;
(3)、计算变换后的I/Q两路信号的互功率谱谱
其中,是Xq[k]的共轭;
(4)、计算互功率谱的相位谱,得相位差Δθ[k];
(5)、利用三点法进行相位解缠
(5.1)、在Δθ[k]加上2π的整数倍来进行解缠绕;
其中,是Δθ[ki]解缠后的相位值,mi为整数,M为选取的频点个数;
(5.2)、取相位谱上的信噪比SNR最大的三个频点及其相位值,记为(k1,Δθ[k1])、(k2,Δθ[k2])、(k3,Δθ[k3]),并作为未解缠的样本点;
(5.3)、利用前两个样本点计算Δθ[k]的斜率;
(5.4)、计算I/Q两通道的群延时差δ和相位偏置差
设第一个样本点不需要加2π的整数倍,则m1=0;
根据在同一条直线上,可得:
将(4)式代入(6)式,化简可得:
(k3-k1)m2=C+(k2-k1)m3 (7)
其中,变量
将(7)式变形可得:(k3-k1)m2-C必须是(k2-k1)的整数倍,因此,将(7)式改写为:
(k3-k1)m2=C(mod(k2-k1)) (8)
其中,mod(·)表示求余;
利用线性模方程求解(8)式,可得:
(k3-k1)m2=m2,0+(L(k2-k1)/D) (9)
其中,D是(k3-k1)和(k2-k1)的最大公约数,L=0,1,…,D-1;m2,0=A(C/D)(mod(k2-k1)),A为整数;
根据(9)式得到m2,再将m2代入到(5)式,得到Q通道相对于I通道的群延时差δ;将δ代入到(3)式,得到Q通道相对于II通道的相位偏置差
(6)、延时模块根据群延时差δ和相位偏置差对单边谱X′q[k]进行补偿,得到补偿后的单边谱X′qc[k];
利用镜像恢复模块得到X′qc[k]的共轭再将X′qc[k]和组合成双边谱Xqcd[k];
将双边谱Xqcd[k]进行傅里叶反变换,得到补偿后的Q路时域信号xqcd[n]。
2.根据权利要求1所述的一种计算并补偿I/Q信号线性相位失衡的方法,其特征在于,所述整数A和最大公约数D满足关系为:
D=A(k3-k1)+B(k2-k1)
其中,B为整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910705992.0A CN110535575B (zh) | 2019-08-01 | 2019-08-01 | 一种计算并补偿i/q信号线性相位失衡的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910705992.0A CN110535575B (zh) | 2019-08-01 | 2019-08-01 | 一种计算并补偿i/q信号线性相位失衡的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110535575A true CN110535575A (zh) | 2019-12-03 |
CN110535575B CN110535575B (zh) | 2021-05-14 |
Family
ID=68660604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910705992.0A Active CN110535575B (zh) | 2019-08-01 | 2019-08-01 | 一种计算并补偿i/q信号线性相位失衡的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110535575B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113421224A (zh) * | 2021-05-27 | 2021-09-21 | 合肥工业大学 | 一种基于视觉的缆索结构健康监测方法及系统 |
CN113702701A (zh) * | 2021-08-09 | 2021-11-26 | 成都玖锦科技有限公司 | 一种基于梳状波信号的幅相特性测量方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1787500A (zh) * | 2004-12-08 | 2006-06-14 | 上海华虹Nec电子有限公司 | 平衡相位分离器路径信号传递时间差的方法及补偿电路 |
US20070058750A1 (en) * | 2005-09-15 | 2007-03-15 | See Puay Hoe A | Baseband compensation of an offset phase locked loop |
WO2010033345A2 (en) * | 2008-09-18 | 2010-03-25 | Motorola, Inc. | Method and system for radio frequency (rf) group delay compensation in a broadband system |
CN105738696A (zh) * | 2016-04-18 | 2016-07-06 | 天津大学 | 全相位时移相位差频率估计方法及装置 |
CN106254071A (zh) * | 2016-09-27 | 2016-12-21 | 四川航天机电工程研究所 | 一种量子密钥分发的实时相位补偿系统和方法 |
US20170317792A1 (en) * | 2016-05-02 | 2017-11-02 | Keysight Technologies, Inc. | Measuring amplitude and phase response of measurement instrument with binary phase shift keying test signal |
CN107735836A (zh) * | 2015-06-19 | 2018-02-23 | 英特尔公司 | 采用可编程跟踪斜率的数字相位控制 |
CN108011853A (zh) * | 2017-11-27 | 2018-05-08 | 电子科技大学 | 混合滤波器组dac延迟和相位偏移的估计和补偿方法 |
US20180174571A1 (en) * | 2015-09-16 | 2018-06-21 | Kabushiki Kaisha Toshiba | Speech processing device, speech processing method, and computer program product |
-
2019
- 2019-08-01 CN CN201910705992.0A patent/CN110535575B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1787500A (zh) * | 2004-12-08 | 2006-06-14 | 上海华虹Nec电子有限公司 | 平衡相位分离器路径信号传递时间差的方法及补偿电路 |
US20070058750A1 (en) * | 2005-09-15 | 2007-03-15 | See Puay Hoe A | Baseband compensation of an offset phase locked loop |
WO2010033345A2 (en) * | 2008-09-18 | 2010-03-25 | Motorola, Inc. | Method and system for radio frequency (rf) group delay compensation in a broadband system |
CN107735836A (zh) * | 2015-06-19 | 2018-02-23 | 英特尔公司 | 采用可编程跟踪斜率的数字相位控制 |
US20180174571A1 (en) * | 2015-09-16 | 2018-06-21 | Kabushiki Kaisha Toshiba | Speech processing device, speech processing method, and computer program product |
CN105738696A (zh) * | 2016-04-18 | 2016-07-06 | 天津大学 | 全相位时移相位差频率估计方法及装置 |
US20170317792A1 (en) * | 2016-05-02 | 2017-11-02 | Keysight Technologies, Inc. | Measuring amplitude and phase response of measurement instrument with binary phase shift keying test signal |
CN106254071A (zh) * | 2016-09-27 | 2016-12-21 | 四川航天机电工程研究所 | 一种量子密钥分发的实时相位补偿系统和方法 |
CN108011853A (zh) * | 2017-11-27 | 2018-05-08 | 电子科技大学 | 混合滤波器组dac延迟和相位偏移的估计和补偿方法 |
Non-Patent Citations (1)
Title |
---|
MILUTIN PAJOVIC等: "Experimental Demonstration of Multi-Pilot Aided Carrier Phase Estimation for DP-64QAM and DP-256QAM", 《2015 EUROPEAN CONFERENCE ON OPTICAL COMMUNICATION (ECOC)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113421224A (zh) * | 2021-05-27 | 2021-09-21 | 合肥工业大学 | 一种基于视觉的缆索结构健康监测方法及系统 |
CN113702701A (zh) * | 2021-08-09 | 2021-11-26 | 成都玖锦科技有限公司 | 一种基于梳状波信号的幅相特性测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110535575B (zh) | 2021-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4261578B2 (ja) | 無線通信装置及び受信方法 | |
TWI241793B (en) | Method and apparatus for I/Q mismatch calibration in a receiver | |
US7573954B2 (en) | Apparatus and method for compensating IQ imbalance in OFDM system with carrier frequency offset | |
US7782928B2 (en) | Method and apparatus for self-calibration in a mobile transceiver | |
JP4573062B2 (ja) | 位相雑音補正装置及びその方法 | |
US8594241B2 (en) | Estimation and compensation method for IQ imbalance | |
CN101027865B (zh) | 利用自适应模拟正交调制器补偿的通信系统中的差分同相正交延迟补偿的系统及方法 | |
US20020018531A1 (en) | Correction of DC-offset of I/Q modulator | |
US9438463B2 (en) | System for the correction of amplitude and phase errors of in-quadrature signals, corresponding receiver and method | |
CN102377727A (zh) | 利用校准信号补偿i/q失衡的正交接收机及其方法 | |
US20120270516A1 (en) | Receiver for estimating and compensating for in-phase/quadrature mismatch | |
KR100809206B1 (ko) | 직교 주파수 분할 다중 수신기의 시간영역 iq 부정합검출 장치 | |
US9686103B2 (en) | Method for compensating the frequency dependent phase imbalance | |
WO2008018263A1 (fr) | circuit de traitement de signal et procédé de traitement de signal | |
CN114374593A (zh) | 用于WiFi宽带收发通路的IQ失衡补偿方法及应用 | |
TWI241096B (en) | Method and apparatus for I/Q mismatch calibration in a receiver | |
CN110535575A (zh) | 一种计算并补偿i/q信号线性相位失衡的方法 | |
CN107078755A (zh) | 用于发射机校准的电路和系统 | |
CN107547458B (zh) | Iq调制中镜像抑制参数的设置方法、装置及射频拉远单元 | |
TWI280002B (en) | Apparatus and method for calibrating IQ mismatch | |
US20070280380A1 (en) | Method and device for compensating inphase-quadrature (iq) imbalance | |
EP1916764A2 (en) | Method and apparatus for compensating for mismatch occurring in radio frequency quadrature transceiver | |
CN100517986C (zh) | 直接转换接收机 | |
US8982975B2 (en) | Method for signal estimation and compensation and apparatus using the same | |
CN107426121A (zh) | 用于广义频分复用系统的训练序列设计及应用方法和装置 |
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 |