CN105471521B - 一种多音测试信号失真产物抑制方法 - Google Patents
一种多音测试信号失真产物抑制方法 Download PDFInfo
- Publication number
- CN105471521B CN105471521B CN201510802568.XA CN201510802568A CN105471521B CN 105471521 B CN105471521 B CN 105471521B CN 201510802568 A CN201510802568 A CN 201510802568A CN 105471521 B CN105471521 B CN 105471521B
- Authority
- CN
- China
- Prior art keywords
- signal
- distorted
- phase
- value
- power
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/0082—Monitoring; Testing using service channels; using auxiliary channels
- H04B17/0085—Monitoring; Testing using service channels; using auxiliary channels using test signal generators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/38—Synchronous or start-stop systems, e.g. for Baudot code
- H04L25/40—Transmitting circuits; Receiving circuits
- H04L25/49—Transmitting circuits; Receiving circuits using code conversion at the transmitter; using predistortion; using insertion of idle bits for obtaining a desired frequency spectrum; using three or more amplitude levels ; Baseband coding techniques specific to data transmission systems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Electromagnetism (AREA)
- Measuring Phase Differences (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明涉及一种多音测试信号失真产物抑制方法,其特征在于:在基带失真信号频点处,分四次产生幅度与失真信号幅度相等,相位分别为θ′y、θ′y+π、的预失真信号,测量叠加后信号的功率值,根据公式直接计算出失真信号的相位估计值,本发明不需要进行多次迭代即可得到基带失真信号相位估计值,且运行稳定,任何情况下均可通过以上计算方法直接获得相位估计值。
Description
技术领域
本发明涉及一种多音测试信号失真产物抑制方法,适用于电子测试测量领域。
背景技术
双音信号常用于测试设备的三阶互调失真产物大小,以获取三阶互调失真截断点(THD)等非线性指标参数;多音信号凭借其带限连续谱特性,常常用于通信设备性能评估和测试。在卫星测试领域,测试通道的三阶互调特性时需要用到双音信号;在卫星进行EMC试验时,需要多音信号作为上行激励信号,使多路转发器工作在额定状态。
目前,一般使用矢量信号源产生多音信号,其内部集成任意波发生器,可通过播放多音信号基带波形文件来产生多音信号。但是,由于信号源内部器件的非线性,得到的多音信号多带有失真产物(本振泄露,镜像泄露,杂散),这些信号严重影响了多音信号作为测试信号时的动态范围。
针对这些失真产物,一般采用预失真方法进行抑制,以改善信号质量。预失真方法采用等幅反相抵消原理,即在失真信号频点处产生一个幅度相同、相位相反的信号对原失真信号进行抵消。估计失真信号的幅度和相位是预失真方法的核心。
幅度估计一般直接测量失真信号电平和主音信号电平的相对值,在假定系统为线性系统的条件下,可得到基带失真信号相对于主音信号的幅度,进而得到基带失真信号的幅度估计。由于一般测试设备(频谱仪)无法测试相位,因此不能直接测得失真信号的相位值。
《电信技术研究》2011年第4期中的论文《多音测试信号失真产物的预失真抑制方法与实现》给出了一种基带失真信号相位的估计方法,其思想为:分三次在基带失真信号频率处产生有一定相位间隔的信号,通过测量叠加后失真信号的大小,取三次测量中失真产物最小时的对应相位,该相位与失真信号的相反相位最接近;以上次所得信号相位为基准,减小相位间隔,再次进行搜索,得到基带失真信号相位估计值;通过多次迭代,可不断缩小范围,提升估计精度。
该方法主要存在的问题是:(1)该算法的稳定性与初始相位密切相关,当初始相位与基带失真信号相位相隔较小时,可较快定位基带失真信号相位;(2)进行相位估计时需要多次迭代(一般需要7次以上),效率较低;(3)搜索方法可能出现迭代不收敛问题。
发明内容
本发明解决的技术问题是:克服以上缺点,提供一种多音测试信号失真产物抑制方法,该方法稳定不需要迭代即可估计出失真信号相位,进而完成失真产物抑制。
本发明的技术方案是:一种多音测试信号失真产物抑制方法,如图1所示,实现步骤如下:
(1)取得有用信号Ss功率值Ps与失真信号Sd的功率值Pd,结合基带有用信号S′s的功率设置值P′s,估算出基带失真信号S′d的功率值P′d;
(2)在基带失真信号频率处分四次产生初相分别为θ′y、θ′y+π、功率均为P′y=P′d的基带预失真信号S′y,S′y对应的预失真信号Sy与失真信号Sd进行叠加,得到叠加后信号的线性功率值Pd1、Pd2、Pd3、Pd4,根据下式即可估算出基带失真信号S′d的相位值θ′d;
θ′d=θ′y-atan2(Pd4-Pd2,Pd1-Pd3)
atan2是matlab中用于计算四象限反正切的函数;
(3)在基带失真信号频率处,产生初相为θ′d+π、功率为P′d的信号S,即可对原失真信号进行等幅反相抵消,以实现失真产物抑制。
θ′y的取值为任意值,取0方便计算。
P′d=Pd-Ps+P′s,Ps、Pd、P′s、P′d均为对数功率值。
所述的功率值Ps、Pd、Pd1、Pd2、Pd3、Pd4均通过频谱仪测量获得,Ps、Pd为对数功率值,Pd1、Pd2、Pd3、Pd4为线性功率值。
本发明的原理:在基带失真信号频点处,分四次产生幅度与失真信号幅度相等,相位分别为θ′y、θ′y+π、的预失真信号,测量叠加后信号的功率值,根据公式直接计算出失真信号的相位估计值。
本发明具体理论支持推导过程如下:
假设原多音信号中失真信号Sd的功率为Pd,幅度为Ad,相位为θd;该失真信号对应基带信号S′d的功率为P′d,幅度为A′d,相位为θ′d。
首先,使用频谱仪测量取得有用信号Ss的功率值Ps与失真信号Sd的功率值Pd,结合基带有用信号的功率设置值P′s,根据公式1估算出基带失真信号的功率值P′d;
P′d=Pd-Ps+P′s (公式1)
公式1中功率值均为对数功率值。
然后,在基带失真信号频点处叠加功率为P′y=P′d,幅度为A′y,相位为θ′y的预失真信号S′y,该基带预失真信号S′y对应射频预失真信号Sy的功率为Py,幅度为Ay,相位为θy。预失真信号Sy与失真信号Sd叠加后得到信号的线性功率值Pd1满足公式2,矢量叠加示意图见图2。
其中常数k是幅度与功率进行换算的系数。
保持基带预失真信号S′y的幅度不变,相位在θ′y的基础上分别增加π、则射频预失真信号Sy的相位θy也分别旋转π、得到叠加后信号的线性功率值见公式3、公式4、公式5:
功率值Pd1、Pd2、Pd3、Pd4均是线性功率值,并满足公式6、公式7
Pd4-Pd2=4kAdAy sin(θy-θd) (公式6)
Pd1-Pd3=4kAdAy cos(θy-θd) (公式7)
则原多音信号中的失真信号Sd的相位θd可采用公式8进行计算。
θd=θy-atan2(Pd4-Pd2,Pd1-Pd3) (公式8)
atan2是matlab中用于计算四象限反正切的函数,其根据角度的正弦和余弦值来计算反正切,计算结果范围为(-π,π]。
失真信号Sd与预失真信号Sy的相位差等于基带失真信号S′d和基带预失真信号S′y的相位差,即
θd-θy=θ′d-θ′y
则基带失真信号S′d的相位θ′d可由公式9进行计算:
θ′d=θ′y-atan2(Pd4-Pd2,Pd1-Pd3) (公式9)
θ′y是首次叠加预失真信号的初相值,可以为任意实数,为了简化计算,常取θ′y=0。
至此,已经得到基带失真信号S′d的相位值θ′d。在基带失真信号频率处,产生初相为θ′d+π、功率为P′d的信号S,即可对基带失真信号进行等幅反相抵消,进而以实现失真产物抑制。
本发明与现有技术相比的优点在于:
(1)本发明不需要进行多次迭代即可得到基带失真信号相位估计值;
(2)本发明运行稳定,任何情况下均可通过以上计算方法直接获得相位估计值。
附图说明
图1为本发明方法的实现流程图;
图2分别给出了基带预失真信号与基带失真信号叠加示意图和预失真信号与失真信号叠加示意图。
图3是具体实施时设备连接示意图;
图4是未进行失真抑制时的双音信号对应频谱仪截图;
图5是进行失真抑制后的双音信号对应频谱仪截图。
具体实施方式
下面以产生双音信号为例,结合失真信号抑制流程,对本发明实施方式做进一步说明。
1)系统连接方式
系统连接方式见图3。本示例中使用任意波形发生器N8241A产生基带信号,使用矢量信号源E8267D对信号进行调制得到最终的输出信号,使用频谱分析仪N9030A对输出信号进行回采分析。
由于本算法对频率精度要求较高,因此三台设备需要共用10MHz参考,具体连接方式为N9030A的10MHz输出端接至E8267D的10MHz输入端,E8267D的10MHz输出端接至N8241A的10MHz输入端(AUX 10MHz REF IN)。
2)双音信号产生
编程控制N8241A和E8267D产生等幅双音信号,信号频率为1999MHz和2000MHz。对应基带信号频率为-1MHz和0MHz,基带功率均为0dBm,初相位分别为124°和17°(均为随机值);将矢量信号源调制载波频率设置为2000MHz,输出电平为0dBm。频谱仪采集到的频谱见图4,由于采用非对称安排(信号频率关于调制载波频率非对称),此时出现镜像失真信号(1998MHz,2001MHz)。以下内容以1998MHz失真信号抑制过程为例。
3)基带失真信号功率估计
记1999MHz信号为有用信号Ss,1998MHz信号为失真信号Sd。
计算有用信号与失真信号的功率差。有用信号功率Ps取-4.60dBm,失真信号功率Pd为-79.80dBm。基带信号中有用信号功率设置值P′s为0dBm,则基带失真信号的功率P′d为-75.2dBm。
4)基带失真信号相位估计
在基带信号中叠加基带预失真信号,频率为-2MHz、功率P′d为-75.2dBm、初相位θ′y为0°,进行下载。为了更精确的进行功率测量,频谱仪设置为ZERO SPAN模式,RB设置为100Hz,对测量结果进行平均。
基带预失真信号初相位θ′y为0°时,测得叠加后失真信号的功率Pd1为-73.70dBm。
再分别设置基带预失真信号初相位θ′y为90°,180°,270°,进行下载、测量,得到功率结果如下:
Pd2=-75.48dBm,Pd3=-85.81dBm,Pd4=-79.52dBm
将上述对数功率值转换为线性功率值,带入公式9即可计算出基带失真信号相位值θ′d为23.18°。为了进行抵消,预失真信号相位为θ′d+π,即203.18°。
5)失真信号抑制
在产生基带波形时,叠加基带预失真信号,频率为-2MHz,功率为-75.2dBm,相位为203.18°,即可对1998MH处失真信号进行抑制。
使用同样方法也可对2001MHz处信号进行抑制,图5给出了同时对1998MHz和2001MHz失真信号进行抑制后的频谱仪截图。1998MHz处失真信号电平降低8dB,2001MHz处降低39dB。
Claims (3)
1.一种多音测试信号失真产物抑制方法,其特征在于实现步骤如下:
(1)取得有用信号Ss功率值Ps与失真信号Sd的功率值Pd,结合基带有用信号S′s的功率设置值P′s,估算出基带失真信号S′d的功率值P′d;其中,P′d=Pd-Ps+P′s,Ps、Pd、P′s、P′d均为对数功率值;
(2)在基带失真信号频率处分四次产生初相分别为θ′y、θ′y+π、功率均为P′y=P′d的基带预失真信号S′y,S′y对应的预失真信号Sy与失真信号Sd进行叠加,得到叠加后信号的线性功率值Pd1、Pd2、Pd3、Pd4,根据下式即可估算出基带失真信号S′d的相位值θ′d;
θ′d=θ′y-atan2(Pd4-Pd2,Pd1-Pd3)
atan2是matlab中用于计算四象限反正切的函数;
(3)在基带失真信号频率处,产生初相为θ′d+π、功率为P′d的信号S,即可对原失真信号进行等幅反相抵消,以实现失真产物抑制。
2.根据权利要求1所述的多音测试信号失真产物抑制方法,其特征在于:θ′y的取值为任意值。
3.根据权利要求1所述的多音测试信号失真产物抑制方法,其特征在于:所述的功率值PS、Pd、Pd1、Pd2、Pd3、Pd4均通过频谱仪测量获得,Ps、Pd为对数功率值,Pd1、Pd2、Pd3、Pd4为线性功率值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510802568.XA CN105471521B (zh) | 2015-11-19 | 2015-11-19 | 一种多音测试信号失真产物抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510802568.XA CN105471521B (zh) | 2015-11-19 | 2015-11-19 | 一种多音测试信号失真产物抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105471521A CN105471521A (zh) | 2016-04-06 |
CN105471521B true CN105471521B (zh) | 2018-07-24 |
Family
ID=55608873
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510802568.XA Active CN105471521B (zh) | 2015-11-19 | 2015-11-19 | 一种多音测试信号失真产物抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105471521B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1599266A (zh) * | 2003-09-16 | 2005-03-23 | 安德鲁公司 | 无线发射机中滤波器的补偿 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6567762B2 (en) * | 2000-12-22 | 2003-05-20 | Agilent Technologies, Inc. | Dynamic range extension apparatus and method |
US7940198B1 (en) * | 2008-04-30 | 2011-05-10 | V Corp Technologies, Inc. | Amplifier linearizer |
-
2015
- 2015-11-19 CN CN201510802568.XA patent/CN105471521B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1599266A (zh) * | 2003-09-16 | 2005-03-23 | 安德鲁公司 | 无线发射机中滤波器的补偿 |
Also Published As
Publication number | Publication date |
---|---|
CN105471521A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9794006B2 (en) | Envelope tracking RF transmitter calibration | |
US8175549B2 (en) | Closed loop transmitter IQ calibration | |
TWI241793B (en) | Method and apparatus for I/Q mismatch calibration in a receiver | |
CN103297155B (zh) | 失真度估算设备和方法以及移动通信装置 | |
JP2010283589A5 (zh) | ||
US9485598B2 (en) | Method and a device for measuring the amplitude noise and/or phase noise of a signal | |
US9596120B2 (en) | Signal transmission apparatus, distortion compensation apparatus, and signal transmission method | |
US20070237258A1 (en) | Measurement method and arrangement for amplitude and phase synchronization in a polar transmitter | |
JP2011182068A (ja) | べき級数型ディジタルプリディストータとその歪補償制御方法 | |
US8798191B2 (en) | Transmission apparatus, radio communication apparatus, and transmission method | |
US20150241494A1 (en) | Nonlinear distortion detection device and distortion compensation power amplifier | |
CN105471521B (zh) | 一种多音测试信号失真产物抑制方法 | |
KR101664797B1 (ko) | Ofdm신호의 papr를 감쇄하기 위한 방법 및 이를 이용한 ofdm 송신기 | |
US11126693B2 (en) | Method and device for reducing noise in a modulated signal | |
JP2014179900A (ja) | 送信装置 | |
Rawat et al. | Characterization and modeling scheme for harmonics at power amplifier output | |
US10122391B2 (en) | Radio frequency systems and methods for polar phase distortion calibration | |
TW201528693A (zh) | 信號處理的方法及裝置 | |
Haojie et al. | Enhanced Imbalance-Distortion Mitigation and Noise Suppression in I/Q-Based Phase Demodulation Systems | |
US10666303B1 (en) | Spurious signal cancellation in radio frequency transmitters | |
CN105610762B (zh) | 一种可校准的宽载波正交调制装置 | |
Simoneau et al. | Multitone feedback through demodulating log detector for detection of spurious emissions in software radio | |
TW201914273A (zh) | 處理訊號干擾的裝置及方法 | |
RU2726184C1 (ru) | Способ формирования предварительно искаженного сигнала | |
JP2016119601A (ja) | 歪み補償回路 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |