CN106125547A - 一种基于分段多项式的功放非线性自适应模拟方法及系统 - Google Patents

一种基于分段多项式的功放非线性自适应模拟方法及系统 Download PDF

Info

Publication number
CN106125547A
CN106125547A CN201610435897.XA CN201610435897A CN106125547A CN 106125547 A CN106125547 A CN 106125547A CN 201610435897 A CN201610435897 A CN 201610435897A CN 106125547 A CN106125547 A CN 106125547A
Authority
CN
China
Prior art keywords
signal
input
power amplifier
input signal
amplitude
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
Application number
CN201610435897.XA
Other languages
English (en)
Other versions
CN106125547B (zh
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.)
CETC 38 Research Institute
Original Assignee
CETC 38 Research Institute
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 CETC 38 Research Institute filed Critical CETC 38 Research Institute
Priority to CN201610435897.XA priority Critical patent/CN106125547B/zh
Publication of CN106125547A publication Critical patent/CN106125547A/zh
Application granted granted Critical
Publication of CN106125547B publication Critical patent/CN106125547B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Amplifiers (AREA)

Abstract

本发明涉及一种基于分段多项式的功放非线性自适应模拟方法,包括:根据输入的数字基带信号,获得输入信号的幅度;将功放非线性特性分为三段;将三段多项式系数的初值存入ROM查找表中;查表得到对应段的多项式系数;获得输出信号的幅度和相位,合成获得期望输出信号;比较期望输出信号和反馈的部分模拟输出信号,获得误差信号,不断更新查找表中的多项式系数,最终使模拟输出信号逼近期望输出信号。本发明还公开了一种基于分段多项式的功放非线性自适应模拟系统。本发明采用基于分段多项式的模拟方法,根据输入信号幅度,将功放非线性特性进行分段描述,每一段仅需采用一个低阶多项式进行逼近,这不仅提高了模拟精度,还降低了运算复杂度。

Description

一种基于分段多项式的功放非线性自适应模拟方法及系统
技术领域
本发明涉及飞行器测控系统技术领域,尤其是一种基于分段多项式的功放非线性自适应模拟方法及系统。
背景技术
近年来,全球通信技术发展日新月异,尤其是无线通信呈现如火如荼的发展态势,而射频功率放大器是无线通信系统中的关键器件,也是实现难度最大、价格最昂贵的部分。它用于对已调制信号进行功率放大,使其获得较高的信号功率值,满足远距离传输的需求,一般而言,当输入信号功率较小时,功放具有良好的线性度。但在实际工程中,某些特殊场合下,例如中继卫星通信系统,为提高工作效率,常要求功放工作在饱和点附近,其非线性较为严重,能够引起信号非线性失真。此外,随着无线通信技术的发展,线性数字、多载波调制方式,例如宽带码分多址(WCDMA)、多载波码分多址(MC-CDMA)和正交频分复用(OFDM)等获得了越来越广泛的应用,这些调制方式产生的信号具有较大波动的包络和较高的峰均功率比,通过射频功放时将会产生新的频率分量,引起严重的带内和带外失真,导致接收系统误码率增加,造成邻信道干扰等,影响整个通信系统质量。
信道模拟器能够模拟真实信道场景,包括时延、多普勒、衰落、色散和噪声等特性。在测控和中继卫星通信等系统中,常被用来进行基带和卫星设备入网前的验证试验、系统对接试验等,以检测设备性能和整个系统的通信质量,成为无线通信发展中必不可少的设备之一。考虑射频功放非线性特性对无线信道通信质量产生的影响不容忽略,因此信道模拟器要准确模拟信道特性、正确评估通信系统质量,则需对功放的非线性特性进行模拟。
传统的功放非线性模拟方法一般采用前置法,即事先对要模拟的非线性功放进行离线建模,然后将模型放置在射频线性功放前端即可。这种方法存在局限性,仅能模拟某一种功放非线性特性,同时由于功放非线性容易受到温度和信号等因素的影响,传统模拟方法不具备适应环境变化的能力。对功放特性进行建模是模拟功放非线性的基础。目前较为常用的建模方法主要有:查找表法和多项式法。查找表法能够逼近各种功放特性曲线,查找表越大,量化误差越少,但收敛速度越慢。多项式法收敛速度快,节约存储空间,其阶数和记忆深度越大,精度就越高,但会导致计算量急剧增加,系统实现困难。两种方法都存在一定的局限性,至今还没出现一种能够被广泛接受的建模方法。
发明内容
本发明的首要目的在于提供一种利用自适应学习结构,使其能够模拟不同的功放非线性特性并能适应环境变化,同时采用多个低阶多项式对功放非线性进行分段建模,提高模拟精度,满足信道模拟器对功放非线性特性模拟需求的基于分段多项式的功放非线性自适应模拟方法。
为实现上述目的,本发明采用了以下技术方案:一种基于分段多项式的功放非线性自适应模拟方法,该方法包括下列顺序的步骤:
(1)根据输入的数字基带信号xI(n)和xQ(n),获得输入信号的幅度|x(n)|;
(2)根据输入信号幅度|x(n)|,将功放非线性特性分为三段,每一段采用一个低阶奇次多项式表示,如公式(2)所示:
z ( n ) = Σ m = 0 M Σ k = 1 ( K + 1 ) / 2 a 2 k - 1 , m | x ( n - m ) | 2 ( k - 1 ) x ( n - m ) - - - ( 2 )
其中,a称为多项式系数,K为非线性阶数,M为记忆深度,k为当前非线性阶数,n是输入信号的离散化采样时刻,z(n)是当前输出信号,x(n-m)是输入信号的第m个延迟分量;
(3)将三段多项式系数的初值存入ROM查找表中;
(4)根据输入信号幅值|x(n)|,查表得到对应段的多项式系数;
(5)输入信号经被模拟非线性功放模块,根据输入信号幅值|x(n)|进行查表,获得输出信号的幅度和相位,然后将幅度和相位合成获得期望输出信号
(6)比较期望输出信号和反馈的部分模拟输出信号获得误差信号自适应训练模块根据误差信号利用基于LMS最小均方误差的自适应更新算法不断更新查找表中的多项式系数,最终使模拟输出信号逼近期望输出信号
在步骤(1)中,所述输入信号的幅度|x(n)|的计算公式如下:
| x ( n ) | = x I ( n ) 2 + x Q ( n ) 2 - - - ( 1 ) .
在步骤(2)中,将功放非线性特性分为三段的分段方法取决于功放非线性特性曲线的形状,即AM-AM和AM-PM曲线;AM-AM曲线中横坐标为输入信号幅度,纵坐标为输出信号幅度;AM-PM曲线中纵坐标为输出信号与输入信号的相位差;输入信号幅度较小时,AM-AM曲线近似于线性,AM-PM曲线纵坐标值近似常量,曲线斜率变化慢,分段区间较大,多项式阶数小;当输入信号幅度较大时,AM-AM和AM-PM曲线斜率变化剧烈,分段区间较小,多项式阶数高。
在步骤(3)中,所述多项式系数为复系数,分为实部和虚部,建立ROM查找表的方法如下:
对公式(2)进行转换,得到公式(3):
z ( n ) = Σ m = 0 M B m * x ( n - m ) - - - ( 3 )
其中,m当前记忆深度,M是最大记忆深度,x(n-m)是输入信号第m个延迟分量,Bm为FIR滤波器的系数:
Bm=LUTm(|x(n-m)|)=a1m+a3m|x(n-m)|2+…+a2K-1,m|x(n-m)|2(K-1) (4)。
在步骤(4)中,按照公式(2)对输入信号进行拟合得到信号数字信号经过D/A转换和正交调制为射频信号,然后经射频线性功放放大获得实际模拟输出信号,由衰减耦合器耦合部分信号给自适应训练模块,然后根据输入信号幅度查表获得多项式系数,输入信号与多项式系数进行拟合获得模拟输出信号
在步骤(6)中,自适应更新算法的实现方法如下:
当输入信号为实信号时,LMS最小均方误差算法的递推公式为:
A(n+1)=A(n)+2ue(n)X(n) (5)
其中,A(n)为多项式系数向量,X(n)为实输入信号向量,e(n)为瞬时误差,u为自适应收敛步长;
当输入信号为复信号时,瞬时误差为通过分别求解瞬时平方误差对AI(n)和AQ(n)的偏导数获得递推公式如下:
A I ( n + 1 ) = A I ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A I ( n ) A Q ( n + 1 ) = A Q ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A Q ( n ) - - - ( 6 )
其中,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,化简得到:
A I ( n + 1 ) = A I ( n ) - 2 u × ( e I ( n ) X I ( n ) + e Q ( n ) X Q ( n ) ) A Q ( n + 1 ) = A Q ( n ) - 2 u × ( e Q ( n ) X I ( n ) - e I ( n ) X Q ( n ) ) - - - ( 7 )
其中,XI(n)和XQ(n)分别为输入信号向量的同相和正交分量,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,收敛的充分条件为:
0≤u<(λmax)-1 (8)
其中,λmax为输入信号向量自相关矩阵的最大特征值。
本发明的另一目的在于提供一种基于分段多项式的功放非线性自适应模拟系统,包括非线性预处理模块,其输入端接收输入的数字基带信号xI(n)和xQ(n),其输出端与数模转换器D/A的输入端相连,数模转换器D/A的输出端与正交调制器的输入端相连,正交调制器的输出端与线性射频功放的输入端相连,线性射频功放的输出端与衰减耦合器的输入端相连,衰减耦合器的输出端与正交解调器的输入端相连,正交解调器的输出端与模数转换器A/D的输入端相连,模数转换器A/D的输出端通过加法器与自适应训练模块的输入端相连,自适应训练模块的输出端与ROM查找表的第一输入端相连,ROM查找表输出AI和AQ信号至非线性预处理模块;第一地址计算模块、第一地址计算模块的输入端均接收输入的数字基带信号xI(n)和xQ(n),第一地址计算模块的输出端与ROM查找表的第二输入端相连,第二地址计算模块的输出端与被模拟非线性功放模块的输入端相连,被模拟非线性功放模块的输出端与时延模块的输入端相连,时延模块的输出端通过加法器与自适应训练模块相连。
由上述技术方案可知,本发明的优点在于:第一,本发明采用基于分段多项式的模拟方法,根据输入信号幅度,将功放非线性特性进行分段描述,每一段仅需采用一个低阶多项式进行逼近,这不仅提高了模拟精度,还降低了运算复杂度;第二,本发明结合ROM查找表,将多项式系数存储在ROM查找表中,进一步减少了复乘、复加运算量,提高了收敛速度,硬件实现时有效节约逻辑资源;第三,本发明利用自适应学习结构,能够实现对任何功放非线性特性曲线进行逼近,提高了系统的普适性,同时能够适应功放非线性特性因温度、环境和信号等不定因素的影响而产生的变化,具有很高的稳定性。
附图说明
图1为本发明的系统结构框图;
图2为功放的AM-AM特性曲线;
图3为功放的AM-PM特性曲线;
图4为ROM查找表的技术实现结构框图;
图5为输出信号的功率谱密度曲线对比示意图;
图6为收敛速度曲线对比示意图;
图7为本发明的方法流程图。
具体实施方式
如图7所示,一种基于分段多项式的功放非线性自适应模拟方法,该方法包括下列顺序的步骤:
(1)根据输入的数字基带信号xI(n)和xQ(n),获得输入信号的幅度|x(n)|;
(2)根据输入信号幅度|x(n)|,将功放非线性特性分为三段,每一段采用一个低阶奇次多项式表示,如公式(2)所示:
z ( n ) = Σ m = 0 M Σ k = 1 ( K + 1 ) / 2 a 2 k - 1 , m | x ( n - m ) | 2 ( k - 1 ) x ( n - m ) - - - ( 2 )
其中,a称为多项式系数,K为非线性阶数,M为记忆深度,k为当前非线性阶数,n是输入信号的离散化采样时刻,z(n)是当前输出信号,x(n-m)是输入信号的第m个延迟分量;
(3)将三段多项式系数的初值存入ROM查找表中;
(4)根据输入信号幅值|x(n)|,查表得到对应段的多项式系数;
(5)输入信号经被模拟非线性功放模块,根据输入信号幅值|x(n)|进行查表,获得输出信号的幅度和相位,然后将幅度和相位合成获得期望输出信号
(6)比较期望输出信号和反馈的部分模拟输出信号获得误差信号自适应训练模块根据误差信号利用基于LMS最小均方误差的自适应更新算法不断更新查找表中的多项式系数,最终使模拟输出信号逼近期望输出信号
在步骤(1)中,所述输入信号的幅度|x(n)|的计算公式如下:
| x ( n ) | = x I ( n ) 2 + x Q ( n ) 2 - - - ( 1 ) .
在步骤(2)中,将功放非线性特性分为三段的分段方法取决于功放非线性特性曲线的形状,即AM-AM和AM-PM曲线,如图2和图3所示;AM-AM曲线中横坐标为输入信号幅度,纵坐标为输出信号幅度;AM-PM曲线中纵坐标为输出信号与输入信号的相位差;输入信号幅度较小时,AM-AM曲线近似于线性,AM-PM曲线纵坐标值近似常量,曲线斜率变化慢,分段区间较大即分段区间可以大一些,多项式阶数小;当输入信号幅度较大时,AM-AM和AM-PM曲线斜率变化剧烈,分段区间较小即分段区间应当小一些,多项式阶数高。
在步骤(3)中,所述多项式系数为复系数,分为实部和虚部,建立ROM查找表的方法如下:
对公式(2)进行转换,得到公式(3):
z ( n ) = Σ m = 0 M B m * x ( n - m ) - - - ( 3 )
其中,m当前记忆深度,M是最大记忆深度,x(n-m)是输入信号第m个延迟分量,Bm为FIR滤波器的系数:
Bm=LUTm(|x(n-m)|)=a1m+a3m|x(n-m)|2+…+a2K-1,m|x(n-m)|2(K-1) (4)。
由此可见,Bm的值与输入信号幅度|x(n-m)|有关,因此可建立查找表,输入信号幅度作为查找表索引地址。如图4所示,式(3)可以由多个查找表的级联扩展实现,对记忆多项式中每个输入延迟分量都构建一个相应的查找表,然后根据输入信号的幅度,查表得到多项式系数,进行拟合。
在步骤(4)中,按照公式(2)对输入信号进行拟合得到信号数字信号经过D/A转换和正交调制为射频信号,然后经射频线性功放放大获得实际模拟输出信号,由衰减耦合器耦合部分信号给自适应训练模块,然后根据输入信号幅度查表获得多项式系数,输入信号与多项式系数进行拟合获得模拟输出信号
在步骤(6)中,自适应更新算法的实现方法如下:
当输入信号为实信号时,LMS最小均方误差算法的递推公式为:
A(n+1)=A(n)+2ue(n)X(n) (5)
其中,A(n)为多项式系数向量,X(n)为实输入信号向量,e(n)为瞬时误差,u为自适应收敛步长;
当输入信号为复信号时,则不能简单地利用式(5)进行多项式复系数更新,此时多项式复系数的实部和虚部具有不同的更新算法。当输入信号为复信号时,瞬时误差为通过分别求解瞬时平方误差对AI(n)和AQ(n)的偏导数获得递推公式如下:
A I ( n + 1 ) = A I ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A I ( n ) A Q ( n + 1 ) = A Q ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A Q ( n ) - - - ( 6 )
其中,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,化简得到:
A I ( n + 1 ) = A I ( n ) - 2 u × ( e I ( n ) X I ( n ) + e Q ( n ) X Q ( n ) ) A Q ( n + 1 ) = A Q ( n ) - 2 u × ( e Q ( n ) X I ( n ) - e I ( n ) X Q ( n ) ) - - - ( 7 )
其中,XI(n)和XQ(n)分别为输入信号向量的同相和正交分量,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,收敛的充分条件为:
0≤u<(λmax)-1 (8)
其中,λmax为输入信号向量自相关矩阵的最大特征值。
如图1所示,本系统包括非线性预处理模块1,其输入端接收输入的数字基带信号xI(n)和xQ(n),其输出端与数模转换器D/A的输入端相连,数模转换器D/A的输出端与正交调制器的输入端相连,正交调制器的输出端与线性射频功放的输入端相连,线性射频功放的输出端与衰减耦合器的输入端相连,衰减耦合器的输出端与正交解调器的输入端相连,正交解调器的输出端与模数转换器A/D的输入端相连,模数转换器A/D的输出端通过加法器与自适应训练模块的输入端相连,自适应训练模块的输出端与ROM查找表的第一输入端相连,ROM查找表输出AI和AQ信号至非线性预处理模块1;第一地址计算模块、第一地址计算模块的输入端均接收输入的数字基带信号xI(n)和xQ(n),第一地址计算模块的输出端与ROM查找表的第二输入端相连,第二地址计算模块的输出端与被模拟非线性功放模块的输入端相连,被模拟非线性功放模块的输出端与时延模块的输入端相连,时延模块的输出端通过加法器与自适应训练模块相连。
基带数字输入信号为I和Q两路正交信号,一方面经过非线性预处理模块1,根据输入信号幅度,查找对应段的多项式系数并进行乘加运算,然后经数模转换和正交调制后输出给线性射频功放,放大后获得模拟输出信号,经过衰减耦合器进行功率耦合,耦合回来的部分模拟输出信号经正交解调、模数转换后生成数字基带信号;另一方面输入信号经过被模拟非线性功放模块,该功放模块为用户事先测量得到被模拟非线性功放的AM-AM和AM-PM曲线,存储在ROM查找表中,根据输入信号幅度查表得到获得期望输出信号;期望输出信号进行延时处理,与耦合回来的模拟输出信号同时进入自适应训练模块,比较得到误差信号,然后根据误差信号,利用自适应算法不断更新ROM查找表中的多项式系数,直至误差达到最小,即实现对功放非线性特性的模拟。
以下对本发明进行试验验证:
仿真条件:数字基带系统输入信号为非恒包络双载波宽带信号,并通过带通滤波整形处理,其中输入信号幅度为0.9,采样率为409MHz,载波频率为20.5MHz,码速率为10.25MHz,滤波器为50阶带通滤波器,其工作带宽为10.25MHz;
仿真内容:非线性预处理部分则分别采用无记忆saleh模型、记忆多项式模型,以及本发明提供的分段多项式法进行仿真并对比其模拟性能,分段方法为:根据输入信号幅度将功放非线性分为三段,多项式系数初值存入查找表中:
a)幅度区间,采用M=3、K=1记忆多项式,参数初始化:
a10=0.1+0.1j,a11=0.01+0.01j,a12=0.01+0.01j
a30=0,a31=0,a32=0
a50=0,a51=0,a52=0
b)幅度区间,采用M=3、K=5记忆多项式,参数初始化:
a10=0.1+0.1j,a11=0.01+0.01j,a12=0.01+0.01j
a30=0.1+0.1j,a31=0.001+0.001j,a32=0.001+0.001j
a50=0.1+0.1j,a51=0.0001+0.0001j,a52=0.0001+0.0001j
c)幅度区间,采用M=3、K=1记忆多项式,参数初始化;
a10=0.1+0.1j,a11=0.01+0.01j,a12=0.01+0.01j
a30=0,a31=0,a32=0
a50=0,a51=0,a52=0
仿真结果:图5给出了收敛后,模拟输出信号的功率谱密度曲线。所有非线性模型都产生了新的频率分量,引起频谱扩展。无记忆多项式法因为没有考虑记忆效应,模拟效果最差,记忆多项式法随着阶数和记忆深度的增加,模拟精度增大。但由图6的收敛速度曲线可知,记忆多项式法需要200个点才能收敛,收敛速度下降,而本发明提供的分段多项式法,在保证模拟精度的前提下,因为减少了复乘复加运算量,很大程度提高了收敛速度,在数字电路实现时,能够有效节约硬件逻辑资源。
综上所述,本发明采用基于分段多项式的模拟方法,根据输入信号幅度,将功放非线性特性进行分段描述,每一段仅需采用一个低阶多项式进行逼近,这不仅提高了模拟精度,还降低了运算复杂度;本发明结合ROM查找表,将多项式系数存储在ROM查找表中,进一步减少了复乘、复加运算量,提高了收敛速度,硬件实现时有效节约逻辑资源;本发明利用自适应学习结构,能够实现对任何功放非线性特性曲线进行逼近,提高了系统的普适性,同时能够适应功放非线性特性因温度、环境和信号等不定因素的影响而产生的变化,具有很高的稳定性。

Claims (7)

1.一种基于分段多项式的功放非线性自适应模拟方法,该方法包括下列顺序的步骤:
(1)根据输入的数字基带信号xI(n)和xQ(n),获得输入信号的幅度|x(n)|;
(2)根据输入信号幅值|x(n)|,将功放非线性特性分为三段,每一段采用一个低阶奇次多项式表示,如公式(2)所示:
z ( n ) = Σ m = 0 M Σ k = 1 ( K + 1 ) / 2 a 2 k - 1 , m | x ( n - m ) | 2 ( k - 1 ) x ( n - m ) - - - ( 2 )
其中,a称为多项式系数,K为非线性阶数,M为记忆深度,k为当前非线性阶数,n是输入信号的离散化采样时刻,z(n)是当前输出信号,x(n-m)是输入信号的第m个延迟分量;
(3)将三段多项式系数的初值存入ROM查找表中;
(4)根据输入信号幅值|x(n)|,查表得到对应段的多项式系数;
(5)输入信号经被模拟非线性功放模块,根据输入信号幅值|x(n)|进行查表,获得输出信号的幅度和相位,然后将幅度和相位合成获得期望输出信号
(6)比较期望输出信号和反馈的部分模拟输出信号获得误差信号自适应训练模块根据误差信号利用基于LMS最小均方误差的自适应更新算法不断更新查找表中的多项式系数,最终使模拟输出信号逼近期望输出信号
2.根据权利要求1所述的基于分段多项式的功放非线性自适应模拟方法,其特征在于:在步骤(1)中,所述输入信号的幅度|x(n)|的计算公式如下:
| x ( n ) | = x I ( n ) 2 + x Q ( n ) 2 - - - ( 1 ) .
3.根据权利要求1所述的基于分段多项式的功放非线性自适应模拟方法,其特征在于:在步骤(2)中,将功放非线性特性分为三段的分段方法取决于功放非线性特性曲线的形状,即AM-AM和AM-PM曲线;AM-AM曲线中横坐标为输入信号幅度,纵坐标为输出信号幅度;AM-PM曲线中纵坐标为输出信号与输入信号的相位差;输入信号幅度较小时,AM-AM曲线近似于线性,AM-PM曲线纵坐标值近似常量,曲线斜率变化慢,分段区间较大,多项式阶数小;当输入信号幅度较大时,AM-AM和AM-PM曲线斜率变化剧烈,分段区间较小,多项式阶数高。
4.根据权利要求1所述的基于分段多项式的功放非线性自适应模拟方法,其特征在于:在步骤(3)中,所述多项式系数为复系数,分为实部和虚部,建立ROM查找表的方法如下:
对公式(2)进行转换,得到公式(3):
z ( n ) = Σ m = 0 M B m * x ( n - m ) - - - ( 3 )
其中,m当前记忆深度,M是最大记忆深度,x(n-m)是输入信号第m个延迟分量,Bm为FIR滤波器的系数:
Bm=LUTm(|x(n-m)|)=a1m+a3m|x(n-m)|2+...+a2K-1,m|x(n-m)|2(K-1) (4)。
5.根据权利要求1所述的基于分段多项式的功放非线性自适应模拟方法,其特征在于:在步骤(4)中,按照公式(2)对输入信号进行拟合得到信号数字信号经过D/A转换和正交调制为射频信号,然后经射频线性功放放大获得实际模拟输出信号,由衰减耦合器耦合部分信号给自适应训练模块,然后根据输入信号幅度查表获得多项式系数,输入信号与多项式系数进行拟合获得模拟输出信号
6.根据权利要求1所述的基于分段多项式的功放非线性自适应模拟方法,其特征在于:在步骤(6)中,自适应更新算法的实现方法如下:
当输入信号为实信号时,LMS最小均方误差算法的递推公式为:
A(n+1)=A(n)+2ue(n)X(n) (5)
其中,A(n)为多项式系数向量,X(n)为实输入信号向量,e(n)为瞬时误差,u为自适应收敛步长;
当输入信号为复信号时,瞬时误差为通过分别求解瞬时平方误差对AI(n)和AQ(n)的偏导数获得递推公式如下:
A I ( n + 1 ) = A I ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A I ( n ) A Q ( n + 1 ) = A Q ( n ) - u × ∂ ( e I 2 ( n ) + e Q 2 ( n ) ) ∂ A Q ( n ) - - - ( 6 )
其中,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,化简得到:
A I ( n + 1 ) = A I ( n ) - 2 u × ( e I ( n ) X I ( n ) + e Q ( n ) X Q ( n ) ) A Q ( n + 1 ) = A Q ( n ) - 2 u × ( e Q ( n ) X I ( n ) - e Q ( n ) X Q ( n ) ) - - - ( 7 )
其中,XI(n)和XQ(n)分别为输入信号向量的同相和正交分量,AI(n)为多项式系数向量A(n)的实部分量,AQ(n)为多项式系数向量A(n)的虚部分量,u为自适应收敛步长,eI(n)为瞬时误差的实部分量,eQ(n)为瞬时误差的虚部分量,收敛的充分条件为:
0≤u<(λmax)-1 (8)
其中,λmax为输入信号向量自相关矩阵的最大特征值。
7.实施权利要求1至6中任一项所述的基于分段多项式的功放非线性自适应模拟方法的系统,其特征在于:包括非线性预处理模块,其输入端接收输入的数字基带信号xI(n)和xQ(n),其输出端与数模转换器D/A的输入端相连,数模转换器D/A的输出端与正交调制器的输入端相连,正交调制器的输出端与线性射频功放的输入端相连,线性射频功放的输出端与衰减耦合器的输入端相连,衰减耦合器的输出端与正交解调器的输入端相连,正交解调器的输出端与模数转换器A/D的输入端相连,模数转换器A/D的输出端通过加法器与自适应训练模块的输入端相连,自适应训练模块的输出端与ROM查找表的第一输入端相连,ROM查找表输出AI和AQ信号至非线性预处理模块;第一地址计算模块、第一地址计算模块的输入端均接收输入的数字基带信号xI(n)和xQ(n),第一地址计算模块的输出端与ROM查找表的第二输入端相连,第二地址计算模块的输出端与被模拟非线性功放模块的输入端相连,被模拟非线性功放模块的输出端与时延模块的输入端相连,时延模块的输出端通过加法器与自适应训练模块相连。
CN201610435897.XA 2016-06-16 2016-06-16 一种基于分段多项式的功放非线性自适应模拟方法及系统 Active CN106125547B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610435897.XA CN106125547B (zh) 2016-06-16 2016-06-16 一种基于分段多项式的功放非线性自适应模拟方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610435897.XA CN106125547B (zh) 2016-06-16 2016-06-16 一种基于分段多项式的功放非线性自适应模拟方法及系统

Publications (2)

Publication Number Publication Date
CN106125547A true CN106125547A (zh) 2016-11-16
CN106125547B CN106125547B (zh) 2018-07-03

Family

ID=57469843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610435897.XA Active CN106125547B (zh) 2016-06-16 2016-06-16 一种基于分段多项式的功放非线性自适应模拟方法及系统

Country Status (1)

Country Link
CN (1) CN106125547B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777533A (zh) * 2016-11-24 2017-05-31 北京理工大学 一种基于查找表的功放非线性模拟方法
CN106792866A (zh) * 2016-11-24 2017-05-31 西安电子科技大学 基于分段线性eh模型的wpc系统功率和时间的分配方法
CN108011598A (zh) * 2017-12-26 2018-05-08 中国工程物理研究院电子工程研究所 一种数字预失真自适应处理方法
CN113055323A (zh) * 2021-03-03 2021-06-29 青岛矽昌通信技术有限公司 一种通信系统的数字预失真处理的方法及系统
CN113691267A (zh) * 2021-09-06 2021-11-23 中国工程物理研究院电子工程研究所 一种同时双频带内载波聚合接收端自干扰抑制方法
CN114754844A (zh) * 2022-04-13 2022-07-15 合肥工业大学 流量仪表修正非线性特性是否放大重复性误差的判断方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102111361A (zh) * 2009-12-28 2011-06-29 电子科技大学 基于放大器估计器的查表法自适应预失真器设计
CN105356854A (zh) * 2015-10-13 2016-02-24 重庆临菲电子科技有限公司 组合算法应用于组合学习结构数字预失真系统的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102111361A (zh) * 2009-12-28 2011-06-29 电子科技大学 基于放大器估计器的查表法自适应预失真器设计
CN105356854A (zh) * 2015-10-13 2016-02-24 重庆临菲电子科技有限公司 组合算法应用于组合学习结构数字预失真系统的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HSIN-HUNG CHEN,ET AL.: "Joint Polynomial and Look-Up-Table Predistortion Power Amplifier Linearization", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS》 *
OUALID HAMMI,ET AL.: "Twin Nonlinear Two-Box Models for Power Amplifiers and Transmitters Exhibiting Memory Effects With Application to Digital Predistortion", 《IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS》 *
张鹏等: "基于分段多项式的RF功放预失真线性化新方法", 《北京理工大学学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777533A (zh) * 2016-11-24 2017-05-31 北京理工大学 一种基于查找表的功放非线性模拟方法
CN106792866A (zh) * 2016-11-24 2017-05-31 西安电子科技大学 基于分段线性eh模型的wpc系统功率和时间的分配方法
CN106792866B (zh) * 2016-11-24 2019-11-22 西安电子科技大学 基于分段线性eh模型的wpc系统功率和时间的分配方法
CN106777533B (zh) * 2016-11-24 2020-06-09 北京理工大学 一种基于查找表的功放非线性模拟方法
CN108011598A (zh) * 2017-12-26 2018-05-08 中国工程物理研究院电子工程研究所 一种数字预失真自适应处理方法
CN108011598B (zh) * 2017-12-26 2021-04-27 中国工程物理研究院电子工程研究所 一种数字预失真自适应处理方法
CN113055323A (zh) * 2021-03-03 2021-06-29 青岛矽昌通信技术有限公司 一种通信系统的数字预失真处理的方法及系统
CN113691267A (zh) * 2021-09-06 2021-11-23 中国工程物理研究院电子工程研究所 一种同时双频带内载波聚合接收端自干扰抑制方法
CN113691267B (zh) * 2021-09-06 2023-03-10 中国工程物理研究院电子工程研究所 一种同时双频带内载波聚合接收端自干扰抑制方法
CN114754844A (zh) * 2022-04-13 2022-07-15 合肥工业大学 流量仪表修正非线性特性是否放大重复性误差的判断方法
CN114754844B (zh) * 2022-04-13 2024-05-14 合肥工业大学 流量仪表修正非线性特性是否放大重复性误差的判断方法

Also Published As

Publication number Publication date
CN106125547B (zh) 2018-07-03

Similar Documents

Publication Publication Date Title
CN106125547A (zh) 一种基于分段多项式的功放非线性自适应模拟方法及系统
US8610500B2 (en) Orthogonal basis function set for digital predistorter
EP3089414B1 (en) Digital pre-distortion parameter obtaining method and pre-distortion system
EP2641326B1 (en) Configurable basis-function generation for nonlinear modeling
CN104796091B (zh) 基于分段记忆多项式的功率放大器建模及数字预失真方法
Naraharisetti et al. Efficient least-squares 2-D-cubic spline for concurrent dual-band systems
CN104363191B (zh) 一种跳频通信系统的数字预失真方法
CN103715992A (zh) 基于简化Volterra级数的功放预失真装置及方法
CN101860326A (zh) 基于Volterra级数间接学习型预失真线性化系统
CN104219182A (zh) 基于宽带调制器的射频功放线性化系统
CN104883140A (zh) 基于宽带射频功率放大器的数字预失真装置
Zhu et al. RF power amplifier behavioral modeling using Volterra expansion with Laguerre functions
CN111064439A (zh) 一种改善短波数字预失真性能的系统及方法
CN114189413B (zh) 一种基于fpga的多载波宽带数字预失真装置
CN101350597A (zh) 一种宽带射频功率放大器建模方法
CN102611661B (zh) 基于精确反解记忆多项式模型方程的预失真装置及方法
CN107659273A (zh) 毫米波宽带功率放大器的可重构数字预失真系统及方法
CN103929136A (zh) 宽带移动通信lte系统发射机功率放大器预失真处理方法
CN106093645A (zh) 确定微波部件宽功率变化范围高阶无源互调电平的方法
CN103414671B (zh) 基于单位基函数的2d-dpd迭代削减算法及应用
US10298366B2 (en) Circuit and method for predistortion
Liu et al. Load-mismatch tracking digital predistortion technique for mobile-terminal power amplifiers
Manai et al. Robust digital predistorter for RF power amplifier linearisation
Soury et al. Handling long-term memory effects in X-parameter model
CN102855417B (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
GR01 Patent grant
GR01 Patent grant