CN103414436B - 一种放大器的弱非线性失真计算方法 - Google Patents
一种放大器的弱非线性失真计算方法 Download PDFInfo
- Publication number
- CN103414436B CN103414436B CN201310339185.4A CN201310339185A CN103414436B CN 103414436 B CN103414436 B CN 103414436B CN 201310339185 A CN201310339185 A CN 201310339185A CN 103414436 B CN103414436 B CN 103414436B
- Authority
- CN
- China
- Prior art keywords
- voltage
- omega
- sigma
- order
- place
- 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
Landscapes
- Amplifiers (AREA)
Abstract
本发明涉及一种放大器的弱非线性失真计算方法,包括:设定电路元件的初始值,同时将每一个非线性单元描述为线性单元并联一个电流源;由低到高依次计算一阶、二阶和三阶的电流和电压并保存,高阶的电流和电压的计算使用已经保存的低价数据;利用一阶、二阶和三阶的电流和电压计算Volterra泛函级数序列的失真。这种弱非线性失真计算方法,在保证高速通讯信号品质的同时,大幅度缩短了二极管桥式开关结构放大器弱非线性失真计算耗时,改善超高带宽无线接收器性能。
Description
技术领域
本发明涉及半导体集成电路无线通讯信号处理领域,具体涉及一种放大器的弱非线性失真计算方法。
背景技术
近年来,随着超高速移动通讯系统的不断发展,对低成本超高带宽无线接收器的需求与日剧增。该无线接收器包括置于模数转换器前端的采用二极管桥式开关结构的跟踪保持放大器。为了获得大的输入带宽,该放大器采用差分开环二极管桥式开关结构。失真是限制放大器动态范围的主要因素之一,是放大器性能进一步提升发展受到限制的瓶颈。因此为了提高放大器的性能必须要减小失真效应。一个可行的方法就采取有效的方法把失真计算出来,根据计算结果对失真采取抑制措施。采用高频肖特基二极管模型对二极管桥式开关放大器的失真进行分析,分析显示失真主要由频率无关的非线性成分所主导。提出了一个依据维他里(Volterra)泛函级数序列分析的用于计算该放大器弱非线性失真的计算机程序。本发明对改善超高带宽无线接收器性能和促进半导体集成电路发展有积极作用。
发明内容
本发明需要解决的技术问题是,如何提供一种放大器的弱非线性失真计算方法,能在改善超高带宽无线接收器性能的基础上降低失真计算复杂度,从而在保证高速通讯信号品质的同时,大幅度缩短了弱非线性失真计算耗时。
本发明的上述技术问题这样解决:构建一种放大器的弱非线性失真计算方法,其特征在于,包括以下步骤:
(一)建立电路模型:设定电路元件的初始值,同时将每一个非线性单元描述为线性单元并联一个电流源;
(二)计算电流和电压:由低到高依次计算一阶、二阶和三阶的电流和电压并保存,高阶的电流和电压的计算使用已经保存的低价数据,具体公式包括:
一阶电流:
i=g1υ+g2υ2+g3υ3+...+gnvn
其中,i是一阶小信号增量电流,g1,g2,g3…gn是泰勒序列系数,v是一阶电压;n是自然数,等于电路模型中的电容总数;
一阶电压:
q=C1v+C2v2+C3v3+...;
其中,Q是总的电量,Vs,q是q处的电压,ωq是q处的频率,t是时间变量,C1,C2,C3。。。Cn是电容;
二阶电流:
其中,V1(t)是一阶电压,g2是泰勒序列系数,v1,q1,v1,q2分别是q1,q2处的电压,ωq1,ωq2是q1,q2处的频率;
二阶电压:
其中,R是系统输入阻抗,g1,g2是泰勒序列系数,vs,q1,vs,q2分别是q1,q2处的源电压,ωq1,ωq2是q1,q2处的频率;
三阶电流:
其中,R是系统输入阻抗,g1,g2,g3是泰勒序列系数,vs,q1,vs,q2,vs,q3分别是q1,q2处的源电压,ωq1,ωq2,ωq3是q1,q2,q3处的频率;
三阶电压:
其中,R是系统输入阻抗,g1是泰勒序列系数;
(三)计算失真特性:利用一阶、二阶和三阶的电流和电压计算Volterra泛函级数序列的失真,具体公式包括:
其中,Vin为输入电压,H1(jω1)和H3(jω1,jω2,jω3)分别是一阶和三阶Volterra核,运算符"o"表示在Vin中每一项的幅值和相位都随着的幅值和相位而改变。
按照本发明提供的弱非线性失真计算方法,所述弱非线性失真计算方法还包括对计算出来的放大器失真采取抑制。
按照本发明提供的弱非线性失真计算方法,所述放大器是二极管桥式开关结构放大器,所述弱非线性失真计算方法应用在超高带宽无线接收器的控制程序中。
本发明提供的放大器的弱非线性失真计算方法,基于弱非线性电路的Volterra序列分析提出一种非线性电流算法来降低计算复杂度,在保证高速通讯信号品质的同时,大幅度缩短了弱非线性失真计算耗时。
附图说明
下面结合附图和具体实施例进一步对本发明进行详细说明。
图1是本发明针对二极管桥式开关结构放大器的弱非线性失真的非线性电流算法示意图;
图2是“SpectreRF”工具仿真与本发明程序计算二极管桥型跟踪保持放大器三次谐波失真(HD3)值的比较;
图3是“SpectreRF”工具仿真与本发明程序计算二极管桥型跟踪保持放大器三阶交调点(IIP3)值的比较;
图4是用于“SpectreRF”工具仿真和本发明程序计算比较的二极管桥型跟踪保持放大器;
图5是“SpectreRF”工具仿真时间与本发明非线性电流程序计算所需时间比较。
具体实施方式
首先,说明本发明二极管桥式开关结构放大器的弱非线性失真计算方法,该方法通过以下步骤实现:
101)首先对弱非线性电路的一阶电流响应进行计算,一个简化压控电导的小信号I/V特性的计算公式如下:
i=g1υ+g2υ2+g3υ3+...(1)
其中第一阶电压υ(t):
102)根据一阶电流响应的结果计算二阶电流响应、三阶电流响应和更高阶电流响应,二阶电流的计算公式如下:
103)通过详细的计算流程,在每一阶均计算出混合频率向量、相应的电压、激励电流和最终的节点电压。
第二,说明本发明的程序具体实现:
该程序基于非线性电流方法和基尔霍夫电流定律(KCL)采用易于移植的“Python”计算机程序设计语言编写。程序采用线性电阻,二极管,电容,电感,可控源,传输线和非线性电阻,二极管,电容,电感,可控源作为基本单元来计算Volterra序列的失真,程序的计算流程图如图1所示。
程序启动后首先初始化环境变量,包括设定电阻,电容,二极管,电感和压控电流源等的初始值,同时用线性等效元件替换非线性元件。然后程序开始用线性等效元件计算一阶电压。输出的数据被保存。接下来程序计算二阶和三阶电流和电压,高阶的电流和电压的计算要用到已经保存的低价数据。正如流程图所示,每一阶相对应的电流和电压值求出来后,期望的相应响应即可从Volterra序列得出。
第三,下面结合具体实施实例和附图对本发明作进一步阐述,但本发明并不限于以下实施实例。
如图1所示,程序启动后首先初始化环境变量,包括设定电阻,电容,二极管,电感和压控电流源等的初始值,同时用线性等效元件替换非线性元件。然后程序开始用线性等效元件计算一阶电压。输出的数据被保存。接下来程序计算二阶和三阶电流和电压,高阶的电流和电压的计算要用到已经保存的低价数据。正如流程图所示,每一阶相对应的电流和电压值求出来后,期望的相应响应即可从Volterra序列得出。
可用SpectreRF仿真器对该二极管桥型放大器进行小信号仿真分析就可以直接测出该放大器的失真特性。
在该仿真中,工艺采用商用的130nmCMOS工艺,输入频率从1GHz变化到10GHz,时钟源CLKP设为4GHz,CLKN是CLKP的互补信号。用于仿真的电路如图4所示,是一个带有采样电容的差分二极管桥型开关放大器。仿真所采用的相关参数为CH=350fF和Id=5mA.在SpectreRF仿真器中运行周期稳态分析(PSS),周期AC分析(PAC)和瞬态分析即可得到失真分析所需的数据。非线性电流算法和SpectreRF仿真的IIP3结果如图3所示。由图可以看出,两者的结果很相似,因此该IIP3分析证实了我们程序计算的正确性。非线性电流算法和SpectreRF仿真的HD3结果如图2所示。同样两者的结果很相似,从而证实了我们分析的正确性。
在相同的硬件条件下,非线性电流算法程序计算的时间和SpectreRF仿真所需的时间如下面的如图5所示。显然在某些特定条件下,程序计算的时间远小于SpectreRF仿真,其执行效率比较高。
以上所述仅为本发明的较佳实施例,凡依本发明权利要求范围所做的均等变化与修饰,皆应属本发明权利要求的涵盖范围。
Claims (3)
1.一种放大器的弱非线性失真计算方法,其特征在于,包括以下步骤:
㈠建立电路模型:设定电路元件的初始值,同时将每一个非线性单元描述为线性单元并联一个电流源;
㈡计算电流和电压:由低到高依次计算一阶、二阶和三阶的电流和电压并保存,高阶的电流和电压的计算使用已经保存的低价数据,具体公式包括:
一阶电流:
i=g1v+g2v2+g3v3+...+gnvn
其中,i是一阶小信号增量电流,g1,g2,g3…gn是泰勒序列系数,v是一阶电压;n是自然数,等于电路模型中的电容总数;
一阶电压:
q=C1v+C2v2+C3v3+...+Cnvn;
其中,Q是总的电量,Vs,q是q处的电压,ωq是q处的频率,t是时间变量,C1,C2,C3…Cn是电容;
二阶电流:
其中,v1(t)是一阶电压,g2是泰勒序列系数,v1,q1,v1,q2分别是q1,q2处的电压,ωq1,ωq2是q1,q2处的频率;
二阶电压:
其中,R是系统输入阻抗,g1,g2是泰勒序列系数,vs,q1,vs,q2分别是q1,q2处的源电压,ωq1,ωq2是q1,q2处的频率;
三阶电流:
其中,R是系统输入阻抗,g1,g2,g3是泰勒序列系数,vs,q1,vs,q2,vs,q3分别是q1,q2,q3处的源电压,ωq1,ωq2,ωq3是q1,q2,q3处的频率;
三阶电压:
其中,R是系统输入阻抗,g1是泰勒序列系数;
㈢计算失真特性:利用一阶、二阶和三阶的电流和电压计算Volterra泛函级数序列的失真,具体公式包括:
其中,Vin为输入电压,H1(jω1)和H3(jω1,jω2,jω3)分别是一阶和三阶Volterra核,运算符"o"表示在Vin中每一项的幅值和相位都随着的幅值和相位而改变。
2.根据权利要求1所述弱非线性失真计算方法,其特征在于,所述弱非线性失真计算方法还包括对计算出来的放大器失真采取抑制。
3.根据权利要求1或2所述弱非线性失真计算方法,其特征在于,所述放大器是二极管桥式开关结构放大器,所述弱非线性失真计算方法应用在超高带宽无线接收器的控制程序中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310339185.4A CN103414436B (zh) | 2013-08-06 | 2013-08-06 | 一种放大器的弱非线性失真计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310339185.4A CN103414436B (zh) | 2013-08-06 | 2013-08-06 | 一种放大器的弱非线性失真计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103414436A CN103414436A (zh) | 2013-11-27 |
CN103414436B true CN103414436B (zh) | 2016-03-30 |
Family
ID=49607426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310339185.4A Active CN103414436B (zh) | 2013-08-06 | 2013-08-06 | 一种放大器的弱非线性失真计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103414436B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1557046A (zh) * | 2002-11-14 | 2004-12-22 | 株式会社日立国际电气 | 失真补偿电路、失真补偿信号生成方法及功率放大器 |
CN1842223A (zh) * | 2005-03-29 | 2006-10-04 | 索尼株式会社 | 音频信号放大设备和失真校正方法 |
CN102111361A (zh) * | 2009-12-28 | 2011-06-29 | 电子科技大学 | 基于放大器估计器的查表法自适应预失真器设计 |
-
2013
- 2013-08-06 CN CN201310339185.4A patent/CN103414436B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1557046A (zh) * | 2002-11-14 | 2004-12-22 | 株式会社日立国际电气 | 失真补偿电路、失真补偿信号生成方法及功率放大器 |
CN1842223A (zh) * | 2005-03-29 | 2006-10-04 | 索尼株式会社 | 音频信号放大设备和失真校正方法 |
CN102111361A (zh) * | 2009-12-28 | 2011-06-29 | 电子科技大学 | 基于放大器估计器的查表法自适应预失真器设计 |
Non-Patent Citations (1)
Title |
---|
射频功率放大器的基带自适应预失真技术;王晓峰;《中国优秀硕士论文全文数据库》;20100515;I135-161 * |
Also Published As
Publication number | Publication date |
---|---|
CN103414436A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | An evolutionary power spectrum model of fully nonstationary seismic ground motion | |
CN106874548A (zh) | 一种基于双重傅里叶变换分析逆变器的方法 | |
Mudunkotuwa et al. | Development of a hybrid simulator by interfacing dynamic phasors with electromagnetic transient simulation | |
Garanayak et al. | Fast and accurate measurement of harmonic parameters employing hybrid adaptive linear neural network and filtered‐x least mean square algorithm | |
CN102623062A (zh) | 一种忆阻器仿真模型 | |
CN102983819A (zh) | 一种功率放大器的模拟方法及功率放大器模拟装置 | |
Wang et al. | Adaptive notch filter design under multiple identical bandwidths | |
Zhu et al. | RF power amplifier behavioral modeling using Volterra expansion with Laguerre functions | |
CN106484962A (zh) | 一种基于e类逆变的谐振式无线输电系统的符号分析方法 | |
Wang et al. | Robust dynamic state estimation of power systems with model uncertainties based on adaptive unscented filter | |
Cruz Nunez-Perez et al. | Flexible test bed for the behavioural modelling of power amplifiers | |
CN110414118A (zh) | 一种基于分离式建模的Boost变换器建模方法及应用 | |
Da Costa et al. | Digital filtering of oscillations intrinsic to transmission line modeling based on lumped parameters | |
CN103414436B (zh) | 一种放大器的弱非线性失真计算方法 | |
CN104217043A (zh) | 分析电源网络的方法和装置 | |
Yamanaka et al. | Theoretical study on operational-amplifier-based negative impedance converters with symmetrically allocated impedance elements | |
CN103051293A (zh) | 射频放大器预失真处理方法及其系统 | |
Sotiriadis et al. | RF switched-capacitor power amplifier modeling | |
US20070033000A1 (en) | Efficient non-iterative frequency domain method and system for nonlinear analysis | |
Comberiate et al. | Using the latency insertion method (LIM) to generate X parameters | |
CN103425836A (zh) | 高速cmos电路失真的程序快速算法 | |
Mrad et al. | Conducted EMI prediction for integrated Class D audio amplifier | |
Su et al. | Evaluation study for the integration of electromagnetic transients simulator and transient stability simulator | |
US20030065496A1 (en) | Effective value impedance simulation method and apparatus and effective value impedance simulation program | |
CN108959757B (zh) | 一种线性分组排序方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |