一种数字线性预畸变方法和系统的制作方法

文档序号:7514406阅读:608来源:国知局
专利名称:一种数字线性预畸变方法和系统的制作方法
技术领域
本发明属于通信领域,尤其涉及一种数字线性预畸变方法和系统。
技术背景功率放大器(Power Amplifier, PA)是一种在通信系统广泛使用的非线性器 件,只要输入信号的幅度超出其线性区,输出就会产生非线性失真,从而造成 信号带内失真和邻带信号干扰,因此必须将功率放大器的失真控制在一定的范 围之内。从效率角度考虑,不能只通过增大功率放大器线性区来消除非线性失 真,因此线性化技术应运而生。对于目前笫三代移动通信主流标准WCDMA系 统,其射频信号是多电平的,具有4艮高的功率峰均比(peak-to-average-power ratio, PAPR),这对线性化技术提出了新的挑战。功率放大器的线性化技术主要有负反馈技术、前馈技术和线性预畸变技术。 其中,线性预畸变技术的成本相对较低,是一种比较简单、实用的技术。现有 技术提供了一种模拟射频预畸变方法,该方法具有提高功率放大器效率、成本 低等优点,但是在具体实现时需要使用射频非线性有源器件,对这些器件的控 制和调整是一个复杂的过程。发明内容本发明实施例的目的在于提供一种数字线性预畸变方法,旨在解决现有技 术实现线性预畸变时必须使用射频非线性有源器件,控制过程复杂的问题。为解决上迷技术问题,本发明提出了一种数字线性预畸变方法,包括以下 步骤获取功率放大器的输入、输出信号;对所述输出信号进行数字线性预畸变训练,计算数字线性预畸变系数; 根据所述数字线性预畸变系数,进行预畸变滤波;
判断所述功率放大器的输出信号是否满足线性化要求,是则结束数字线性 预畸变训练。
本发明实施例的另一目的在于提供一种数字线性预畸变系统,所述系统包

预畸变滤波器,用于根据预畸变训练器计算出的数字线性预畸变系数,对 预畸变滤波器的输入数字基带信号进行滤波处理,并且根据功率放大器的特征, 动态调整滤波参数以达到最佳滤波效果;
预畸变训练器,用于对预畸变训练器的输入信号进行预畸变训练,并根据 预畸变滤波器的实际输出信号和该输出信号的估计值之间的误差信号,利用最 小二乘法调整预畸变训练的抽头系数直至算法收敛,获取最终的数字线性预畸 变系数。
在本发明的实施例中,基于Volterra级数的精简模型,通过对功率放大器 的输入、输出信号进行预畸变训练和预畸变滤波,可以自适应地改善功率放大 器的线性度。经过仿真验证,此数字线性预畸变方法对单载波信号可以改善功 放线性度高达15db,而对于高峰均比的多载波信号改善也至少有10db。另夕卜, 数字线性预畸变系统星座图指标也有明显改善,经预畸变处理后功率放大器输 出的星座图接近于信源的星座图。经测试平台(基于3GPP25.141协议)验证, 结果表明,该方法是改善功率放大器线性度的一种高效方法。


图1是本发明实施例4是供的数字线性预畸变系统的结构框图; 图2是本发明实施例提供的数字线性预畸变方法的工作流程图; 图3是本发明实施例提供的预畸变浮点仿真功率谱示意图; 图4是本发明实施例提供的24bits定点仿真功率谱示意图;图5是本发明实施例提供的16bits定点仿真功率谱示意图; 图6是本发明实施例提供的12bits定点仿真功率镨示意图; 图7是本发明实施例提供的数据块的F矩阵的生成过程示意图; 图8是本发明实施例提供的DLP滤波器蝶形结构示意图。 图9是本发明实施例提供的单载波PSD图; 图IO是本发明实施例提供的两载波PSD图; 图11是本发明实施例提供的三载波PSD图; 图12是本发明实施例提供的四载波PSD图13是本发明实施例提供的单载波直接经过功放输出的星座示意图,其中 矢量误差率(Error Vector Magnitude, EVM ) EVM =7.6164%;
图14是本发明实施例提供的单载波经过预畸变和功放输出的星座示意图, 其中EVM = 0.0813%;
图15是本发明实施例提供的两载波直接经过功放输出的星座示意图,其中 EVM =4.6440 %;
图16是本发明实施例提供的两载波经过预畸变和功放输出的星座示意图, 其中EVM = 0.1361%;
图17是本发明实施例提供的三载波直接经过功放输出的星座示意图,其中 EVM = 6.3061%;
图18是本发明实施例提供的三载波经过预畸变和功放输出的星座示意图, 其中EVM = 0.5692%;
图19是本发明实施例提供的四载波直接经过功放输出的星座示意图,其中 EVM = 7.8442%;
图20是本发明实施例提供的四载波经过预畸变和功放输出的星座示意图, 其中EVM = 1.5777%;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实 施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅 仅用以解释本发明,并不用于限定本发明。
在本发明实施例中,基于Volterra级数的精简模型,通过对功率放大器的 输出、输入信号进行预畸变训练和预畸变滤波,可以自适应地改善功率放大器 的线性度。
图1示出了本发明实施例提供的数字线性预畸变(Digital Line Processing, DLP)系统的结构,为了便于说明,仅示出了本发明实施例相关的部分。该数 字线性预畸变系统包括预畸变滤波器ll、预畸变训练器12、 D/A转换器13、 上变频器14、功率放大器15、下变频器16、 A/D转换器17、加减器18。
预畸变滤波器11用于根据预畸变训练器12计算出的DLP系数L对预畸 变滤波器的输入数字基带信号;c(w)进行滤波处理,并且根据功率放大器的特征,
可动态调整滤波参^:以达到最佳滤波效果。
其中,预畸变滤波器11根据预畸变训练器12计算出的DLP系数"采用 蝶形结构运算,对输入的基带信号x(")进行滤波处理,并且根据功率放大器的 特征,动态调整滤波参数,输出线性度得到改善的信号z(")。
预畸变训练器12用于对输入信号;K")进行预畸变训练,并根据预畸变滤波 器11的实际输出信号z(")和z(w)的估计值S(")之间的误差信号e("),利用最小二 乘法调整预畸变训练的抽头系数直至算法收敛,得出最终的DLP系数G。
其中,预畸变训练器12对输入信号y(n)进^f亍预畸变训练,先对该信号进行 数据分块,对每个数据块依次构建F矩阵,计算对应的F"F,然后再利用最小 二乘法得到的公式"(F"F)"F 或者"F、/L/1/,计算出预畸变训练器12 的抽头系数。将DLP系数G拷贝到预畸变滤波器11中进行预畸变滤波处理, 如系统输出满足线性度要求,则结束DLP训练,否则继续进行DLP训练。 D/A转换器13用于实现数模转换。
上变频器14用于实现频率转换。
8功率放大器15用于对输入信号进行功率放大。 下变频器16用于实现频率转换。 A/D转换器17用于实现模数转换。
加减器18用于计算预畸变滤波器实际输出信号z(")和该输出信号的估计值 S(")之间的误差信号
图2示出了本发明提供的数字线性预畸变方法的工作流程,在该实现流程 中,DLP系统获取功率放大器的输入、输出信号,对该输出信号进行DLP训练, 得到DLP系数,再将该系数送入DLP滤波器进行滤波。判断DLP滤波后并送 往功率放大器的输出信号是否满足线性化要求,是则结束DLP训练。其具体的 实现流程详述如下
在步骤S201中,获取功率放大器的输入、输出信号。
作为本发明的一个实施例,DLP系统计算功率力文大器的输入信号z(")和输 出信号少("),并对输出信号:K")进行归一化处理,即将少(")归一化到[O, l]间。
此外,作为本发明的一个实施例,在实现DLP训练时,可以采用数字信号 处理器(Digital Signal Processor, DSP),该DSP可以进行浮点运算,但在利 用现场可编程门阵列(Field Programmable Gate Array, FPGA)实现时,对DLP
训练器的数据:K")不可能采用无限精度,必须对:KW进行量化处理,即用有限字 长表示。为了体现定点量化对预畸变性能的影响,在DLP训练和滤波阶段引入 24、 16以及12比特的定点量化,并在计算过程中,进行比特截位以防止乘加 操作带来的扩位溢出。考虑到性能和资源占用的问题,在FPGA实现时采用定 点16比特量化即可实现DLP线性化的效果。图3至图6示出了本发明实施例 提供的数字线性预畸变系统的仿真功率谱,,人中可以看出,在采用浮点运算和 定点量化进行数字线性预畸变均可以高效改善功率放大器线性度,在使用 FPGA实现时,采用定点16比特量化即可实现DLP线性化的效果。
在步骤S202中,构建F矩阵,计算F〃F。
在进行DLP训练时,首先要构建F矩阵。其具体的过程如下第一步根据基于Volterra级数的精简模型,推导得出DLP模型,其模型 可以表示为
= 2] 2] :K"-《)I ><"-《)I"1
本发明基于Volterra级数的精简模型,提出了有记忆多项式的DLP模型及 其正交化结构,具体的推导过程如下
具有记忆效应的非线性系统,可以用Volterra级数表示为
z(O = Z J"…J"4(^)nKMX (1)
其中,7一[f,,…,^r,反(.)是i:阶Volterra级数的实数部分,在离散时域上, 公式(l)可以表示为
z("^SS…Z、"(",…,")Eb("-夂)n,("') (2)
当取Volterra级凄t的对角^t时,即/, =/2 =." = /24+1 =1时,/^+1(1,...,1)为常数, 记为~,,则公式(2)可以简化为
"")-iliv("-WX"")r (3)
其中,p = 2K+i, p为非线性系统模型的最高阶数,y(")和w")分别为该模 型的输入和输出,L为记忆长度,Zv为系数,这一简化有效地降低了模型的复 杂度。
根据上述Volterra级数的简化模型,本发明实施例所提供的DLP模型可表 示为
a")=t S- ) i x"—《)r (4)
其中,j(")和z(")分别为DLP滤波器的输入和输出信号;2为预畸变滤波 器的记忆长度,c^为DLP的抽头系数。
第二步根据所述DLP模型,构建F矩阵,并计算F"F,其中F〃为F的共 辄转置矩阵。
10在DLP训练方法中,对公式(4)的结构进行了基于正交多项式的等效展 开,减少数据之间的相关性,进一步提高了 DLP的性能。 定义
(")=少("—《)广一1 (5) 则 z = Ra (6)
其中定义i^二[ (0),L, (iV —l)f, Rg=[VL, ],R = [R。,L,Re]。
利用最小二乘法的间接训练方法,可以得到DLP系数
"(R"R)-11^ (7) 由于公式(5)中的非线性多项式y , y|7l, W少卩等具有高度相关性,造成RWR 矩阵的条件数特别大,会影响5的准确性和稳定性。为减少这一相关性,改写 么、式(4)为
"")=f fx似"》 (8)
其中,A(力-^:"汰卟r (9)
/=1
(-1 + /)!
"议
(/-1)!(/ + 1)!("1)! (10) 令 f = [R。u,l ,RqU] (11) 则z = Fb (12) 公式(11)所给出的表达式中的F即为F矩阵。
作为本发明的一个实施例,虽然公式(8)中的正交多项式方法,是在复随 机信号幅度均匀分布在[O, l]时得到的,但是此方法对于未归一化分布的幅度同 样适用。为了充分利用正交多项式的优点,本实施例在进行DLP系统工作时先 将X")归一化到[O, l]间。
在步骤S203中,计算DLP系数。
由公式(12),再利用最小二乘法可以得到DLP系数 G = (FwF)-iFwz (13)其中F"是矩阵F的共辄转置矩阵,可以通过计算F〃F以及F〃的值计算出 DLP系数fe。
从公式(8)可以看出,可以釆用移位的方法来简化算法。首先由R。U产生矩 阵F的前*列,则F的其他列即为这前&列的移位,又因为F"F是Hermitian矩阵, 为了提高运算的效率,采取Cholesky分解产生一个下三角矩阵L,该下三角矩 阵L满足
IX^F〃F (14)
将公式(14)代入公式(13),得到
LL〃"F"z (15)
通过简单的置换,得到dlp系数"fWz/l/lw。
作为本发明的一个实施例,例如对长度为38400个码片的每帧WCDMA数 据进行DLP训练,那么矩阵F需要占用38400行和《*(2 +l)列,但其随后得到 的F"F却只需占用^:*(2 + 1)行和《*(。 + 1)列。因此在实现时,中间处理过程中 的F占用了最多的空间。为了降低空间复杂度,即在RAM里存放较少的数据量, 可以对预畸变训练器的输入信号先进行分块,对每个数据块依次构建F矩阵, 计算对应的F"F。
例如,假设将38400个码片分为20个数据块,每块均有1920个码片(在 实际应用时,可以按照RAM的具体占用情况来灵活设定每个块的大小),再 分别对每个数据块执行DLP训练,得到的F都占用了 1920行和夂*(2 + 1)列, 因此此时占用RAM的最大空间从原来的38400*K*(g+l)降低到了
1920*〖*(e + l),大大降低了空间的占用率。
图7示出了本发明实施例提供的数据块的F矩阵生成过程,其中,方块A 表示每行中的第1,L ,K列数据,B表示每行中的第K+1,L ,2K列数据,……。B=0 表示该方块中的所有数据为0 。 A—B中的箭头表示利用 (")=《)I:K"-《)l"得到方块A中的数据后,直接复制便可得到方块B中 的数据。通过对数据进行分块计算DLP系数的具体实现步骤如下第一步计算第一个数据块的f矩阵及对应的f"f。
输入第一个数据块。由~(")=><"-《)l"可以看出,矩阵f每一行的 第[//:+ l,L ,(/+ l)K]列数据等于前一行中的第[(/- 1)《+ l,L jX]列数据,这里 "l,L,g,而第/= l,L ,e行中的第p:+i,L ,(G+1)尺]列数据,全部等于0,因此
对于第一个数据块,只需要计算第1列到第《列的数据,便可得到第一个数据 块的矩阵F,计算对应的Ft- f"f。该数据块其它列的数据可由该/C列数据直接 复制得到,节省了实现的逻辑。
第二步计算第二个数据块的f矩阵及对应的f〃f。
输入第二个数据块,此时仍只需计算第1列到第X列的数据,再利用复制 的方法得到其他列的数据,而其中z: l,L ,g行中的第p:+l,L ,(g+l)《]列数据,
可由第 一个数据块对应f的最后e行数据复制得到,此时利用新的F更新 Ft= Ft+ F 。
第三步计算其它数据块的f矩阵及对应的f"f。
输入其它数据块,得到其对应f的方法类似于第二个lt据块,每次得到f后, 都对Ft进行更新。
第四步如预畸变训练是每帧数据进行一次,那么应对每帧数据都计算f矩 阵及对应的f〃f。
等每帧中的最后一个数据块计算完毕,Ft也得到更新后,此时的Ft就等于 将该帧38400个码片直接计算得到的fwf,大大降低了空间的占用率。利用 ll" = f"f得到Cholesky算法分解后对应的下三角矩阵l 。
第五步计算DLP的系数。
G= f、/:l/lw
在步骤S204中,根据DLP训练所得到的DLP系数,利用公式(8 )进行 预畸变滤波。
根据DLP训练得到的DLP系数,利用公式(8)进行非线性的DLP滤波。作 为本发明的一个实施例,图8示出了本发明实施例提供的蝶形运算结构,详述
13如下
在WCDMA系统中,根据数据速率和DLP系统的性能,对于载波数小于 等于4的情况,可以取[K,Q〗-[5,2],则
z(") = 4,, (_y (")) + 6函(y (" -1)) + 6函(:v (" - 2)) +62(^2 + ^ W" - 0) + W" - 2))
+63,3 + 63iK W" -!))+6函W"-2)) (13)
+640^4 (y ("))+641^4 (_y (" _ l)) + - 2》
+65,5 (, (")) + 651^5 W" - 0) + 一 W" - 2)) 其中,
^tK")H—"(")+仃(")l:K")1 (15)
^W")),(")-20 | + 15,W")|2 (16) ^CK")) = -10^(") + 60;;(")|y(")|-105y(")|y(")|2+56y(")|j(")|3 (17)
K= 15少(")-140+)|少(")|+420>>(")卜(")|2 -504_y(")|_y(")|3 +210y(")|》'(")|4 (18)
公式(10)所对应的正交多项式系数如下表所示
<formula>formula see original document page 14</formula>在步骤S205中,判断功率放大器的输出信号是否满足线性化要求,是则执 行步骤S206,否则执行步骤S201。
判断经DLP滤波处理后并经过功率放大器的输出信号是否满足线性化要 求,是则执行步骤S206,否则执^f亍步骤S201。图9至图12示出了本发明实施 例提供的数字线性预畸变系统的PSD指标,经过仿真-3t〖正,此数字线性预畸变 方法对单载波信号可以改善功率放大器线性度高达15db,而对于高峰均比的多 载波信号改善也至少有10db。图13至图20示出了本发明实施例提供的数字线性预畸变系统的星座图指标,/人中可以看出,星座图指标也有明显改善,经预 畸变后功率》文大器專命出的星座图接近于信源的星座图。
在步骤S206中,结束DLP训练。
在本发明实施例中,基于Volterra级数的精简模型,通过对功率放大器的 输出、输入信号进行预畸变训练和预畸变滤波,可以自适应地改善功率放大器 的线性度。DLP系统在通过FPGA实现时,采用16比特定点量化,从图3至 图6可以看出,采用16比特定点量化进行定点量化即可获得4^好的DLP效果, 而无需采用占用资源更多的24比特量化。另外,在进4亍DLP训练时利用对邀: 据分块计算F矩阵的方法,每次只需反复利用RAM中的1920*《* (0 +1)个空间, 与原来的38400*〖*(2 + 1)相比,大大减少了空间资源的占用率;并且只需计算前 K列数据,再利用数据块复制的方法得到其它列的凄t据,可以大大减少逻辑资 源的使用率,并提高DLP实现的速度。最后,在进行DLP滤波时采用蝶形结 构,大大简化了计算的复杂度,便于硬件实现,提高了 DLP效率。图9至图 12示出了本发明实施例提供的数字线性预畸变系统的PSD指标,经过仿真验 证,此数字线性预畸变方法对单载波信号可以改善功率放大器线性度高达 15db,而对于高峰均比的多载波信号改善也至少有10db。图13至图20示出了 本发明实施例提供的数字线性预畸变系统的星座图指标,从中可以看出,星座 图指标也有明显改善,经预畸变后功率放大器输出的星座图接近于信源的星座 图。经测试平台(基于3GPP 25.141协议)验证,结果表明,该数字线性预畸 变方法是改善功率放大器线性度的一种高效方法。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发 明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明
的保护范围之内。
权利要求
1、一种数字化线性预畸变方法,其特征在于,所述方法包括下述步骤获取功率放大器的基带输入、输出信号;对所述输出信号进行数字线性预畸变训练,计算DLP系数;根据所述数字线性预畸变系数,进行预畸变滤波;判断所述功率放大器的输出信号是否满足线性化要求,是则结束数字线性预畸变训练。
2、 如权利要求1所述的方法,其特征在于,对所述输出信号进行数字线性 预畸变训练,计算数字线性预畸变系数的步骤具体为构建F矩阵,计算F"F,其中F〃为F的共轭转置矩阵'; 计算数字线性预畸变系数。
3、 如权利要求l所述的方法,其特征在于,对所述输出信号进行数字线性 预畸变训练,计算数字线性预畸变系数的步骤之前,所述方法还包括对功率放大器的输出信号进行归一化处理。
4、 如权利要求l所述的方法,其特征在于,对所述输出信号进行数字线性 预畸变训练,计算数字线性预畸变系数的步骤之前,所述方法还包括在利用现场可编程门阵列实现数字线性预畸变时,对功率放大器的输出信 号采用定点16比特量化。
5、 如权利要求2所述的方法,其特征在于,所述构建F矩阵,计算F"F的 步骤具体为对预畸变训练器的输入信号先进行分块,对每个数据块依次构建F矩阵, 计算对应的F"F;其中F矩阵满足下述公式 F= [R0U,l ,RqU]其中,0为预畸变滤波器的记忆长度,R。U是矩阵F的前a列。
6、 如权利要求2所述的方法,其特征在于,所述计算数字线性预畸变系数 可利用最小二乘法获得,其步骤具体满足下述/>式其中,G为数字线性预畸变系数,F"是矩阵F的共轭转置矩阵,z是预畸变 滤波器的输出信号。
7、 如权利要求6所述的方法,其特征在于,所述计算数字线性预畸变系数 的步骤具体满足下述公式<formula>formula see original document page 3</formula>其中,G为数字线性预畸变系数,F"是矩阵F的转置矩阵,z是预畸变滤波 器的输出信号,L是采取Cholesky分解产生 一 个下三角矩阵,该下三角矩阵L满 足下述公式<formula>formula see original document page 3</formula>其中,1/是三角矩阵L的共轭转置矩阵。
8、 如权利要求1所述的方法,其特征在于,根据所述数字线性预畸变系数,进行预畸变滤波的步骤具体满足下述公式 = E ZX A CK" - g))其中,z(")为预畸变滤波器的输出信号,y(")为预畸变训练器的输入信号, 为预畸变训练器的抽头系数,<formula>formula see original document page 3</formula>其中"_(—l)'十K阔!
9、 如权利要求8所述的方法,其特征在于,所述根据所述数字线性预畸变 系数,进行预畸变滤波的步骤具体为采用蝶形运算结构进行预畸变滤波。
10、 一种数字化线性预畸变系统,其特征在于,所述系统包括 预畸变滤波器,用于根据预畸变训练器计算出的数字线性预畸变系数,对预畸变滤波器的输入数字基带信号进行滤波处理,并且根据功率放大器的特征,动态调整滤波参数以达到最佳滤波效果;预畸变训练器,用于对预畸变训练器的输入信号进行预畸变训练,并根据 预畸变滤波器的实际输出信号和该输出信号的估计值之间的误差信号,利用最 小二乘法调整预畸变训练的抽头系数直至算法收敛,获取最终的数字线性预畸 变系数。
11、 如权利要求IO所述的系统,其特征在于,所述系统还包括 D/A转换器,用于实现数模转换;上变频器,用于实现频率转换; 功率放大器,用于对输入信号进行功率放大; 下变频器,用于实现频率转换; A/D转换器,用于实现模数转换;加减器,用于计算预畸变滤波器实际输出信号z(")和该输出信号的估计值 S(")之间的i吴差4言号。
12、 如权利要求IO所述的系统,其特征在于,所述预畸变滤波器对输入数 字基带信号采用蝶形运算结构进行预畸变滤波。
13、 如权利要求IO所述的系统,其特征在于,所述预畸变训练器对预畸变 训练器的输入信号先进行分块,对每个数据块依次构建F矩阵,计算对应的 FWF,再利用最小二乘法得到的公式"(F"F)"F"z或者"F"z/L/L",计算出 数字线性预畸变系数6。
全文摘要
本发明适用于通信领域,提供了一种数字线性预畸变方法和系统,所述方法包括下述步骤获取功率放大器的基带输入、输出信号;对所述输出信号进行数字线性预畸变训练,计算数字线性预畸变系数;根据所述数字线性预畸变系数,进行预畸变滤波;判断所述功率放大器的输出信号是否满足线性化要求,是则结束数字线性预畸变训练。所述系统包括预畸变滤波器和预畸变训练器。本发明可以很好地改善功率放大器线性度,另外,可以极大地减少空间资源的占用率,减少逻辑资源的使用率,简化计算的复杂度,便于硬件实现,提高数字线性预畸变效率。
文档编号H03F1/32GK101505139SQ20081021735
公开日2009年8月12日 申请日期2008年11月17日 优先权日2008年11月17日
发明者曾庆丰, 杨嗣环 申请人:深圳市云海通讯股份有限公司
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1