发明名称 组合算法应用于组合学结构数字预失真系统的方法
摘要 本发明公开了一种组合算法应用于组合学结构数字预失真系统的方法,采用经典的峰均比抑制CFR技术和数字预失真DPD技术联合方案,同时针对直接学结构和间接学结构的优缺点,提出间接学结构与直接学结构结合的组合学结构设计数字预失真系统,然后结合LMS算法、RLS算法以及变步长LMS即NLMS算法的优势,提出RLS和NLMS联合的组合算法对数字预失真器参数估计和更新,组合算法具有RLS算法收敛迅速的特点,同时简化了计算量。对输入的高峰均比信号先进行CFR处理,再经DPD处理之后送给功放,得到比较理想的线性输出,有效改善了功放非线性引起的带内失真和带外频谱扩展。本发明不仅降低了宽带功放因记忆效应造成的性能下降,而且改善了宽带功放的线性度和效率。
申请公布号 CN105356854A 申请公布日期 2016.02.24
申请号 CN201510677011.8 申请日期 2015.10.13
申请人 重庆临菲电子科技有限公司 发明人 谷林海;葛利嘉;毛键;张振宇
分类号 H03F1/32(2006.01)I;H04L25/49(2006.01)I 主分类号 H03F1/32(2006.01)I
代理机构 重庆大学专利中心 50201 代理人 王翔
主权项 组合算法应用于组合学习结构数字预失真系统的方法,其特征在于,包括两个部分,组合学习结构数字预失真系统,以及组合算法基于搭建的组合学习结构数字预失真系统对信号处理方法;所述组合学习结构数字预失真系统包括前向通路和反馈回路;所述前向通路包括信号产生模块(100)、峰均值抑制模块即CFR模块(101)、上变频模块即DUC模块(102)、数字预失真模块即DPD模块(103)、DSP模块(104)、数模转换模块即DAC模块(105)、带通滤波器(106)、功率放大器即PA(107)以及天线调谐器和天线(109);所述反馈回路包括耦合器(108)、衰减器(110)和模数转换模块即ADC(111);所述DPD模块(103)采用多项式预失真技术,该技术主要包括功放模型、学习结构和自适应算法三个部分;DPD内核是基于记忆多项式模型的数字预失真;DPD学习结构采用组合学习结构,DPD的自适应算法采用组合算法;所述学习结构采用组合学习结构,由间接学习结构和直接学习结构结合的组合学习结构设计数字预失真系统;系统在间接学习结构与直接学习结构之间进行切换;所述组合算法是在间接学习结构下采用RLS和NLMS算法对DPD内核系数进行估计;而在直接学习结构下采取LMS算法对DPD内核系数进行估计和更新;所述组合算法基于所搭建的组合学习结构数字预失真系统对信号处理方法包括以下步骤:1)将基带信号模块生成的基带信号进行数字上变频和峰值削波处理产生x(n);2)采用记忆多项式模型构造数字预失真内核即DPD内核,其数学表述参见公式1;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>z</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>C</mi><mrow><mi>k</mi><mi>q</mi></mrow></msub><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>q</mi><mo>)</mo></mrow><mo>|</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>q</mi><mo>)</mo></mrow><msup><mo>|</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000011.GIF" wi="1301" he="159" /></maths>其中,z(n)为预失真输出信号,K表示模型的最高非线性阶数,Q为记忆深度,C<sub>kq</sub>是DPD系数;3)设误差判断门限|e<sub>1</sub>(n)|=|e(n)|=1,则|e(n)|>0.01,开关掷1,工作在自适应算法1,此时的学习结构为间接学习结构;且|e<sub>1</sub>(n)|>0.05,则RLS算法进行DPD内核系数<img file="FDA0000824189280000021.GIF" wi="103" he="68" />快速估计来其初始值;基于RLS算法包括以下过程;a)设<img file="FDA0000824189280000022.GIF" wi="103" he="64" />为DPD内核系数,初始化<img file="FDA0000824189280000023.GIF" wi="127" he="67" />b)通过公式1计算出信号x(n)经过预失真输出信号z(n);c)功放模型采用Wiener‑Hammerstein模型,其数学表达式参见公式2‑4;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>L</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>a</mi><mi>l</mi></msub><mi>z</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000024.GIF" wi="1037" he="140" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>v</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>b</mi><mi>k</mi></msub><mi>u</mi><msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>k</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000025.GIF" wi="1038" he="134" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>c</mi><mi>q</mi></msub><mi>v</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>q</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000026.GIF" wi="1037" he="150" /></maths>式中,u(n)为输入向量,v(n)为功放输出衰减得到的信号,y(n)为功放的输出采样信号,a<sub>l</sub>是线性时不变系统的冲值响应值,L代表模型记忆深度;b<sub>k</sub>为无记忆非线性系统多项式的系数,K表示模型的最高非线性阶数;c<sub>q</sub>为线性时不变系统的冲值响应值,Q是模型记忆深度;通过公式2、公式3和公式4,计算得到信号z(n)经过功放的输出采样信号为y(n);d)y(n)经尺度变换后y(n)/G,作为“后失真器”的输入,其中G代表功放的增益;令v(n)=y(n)/G,u(n)为输入向量,其中:u(n)=[v(n),v(n)|v(n)|<sup>2</sup>,v(n)|v(n)|<sup>4</sup>,v(n‑1),v(n‑1)|v(n‑1)|<sup>2</sup>,v(n‑1)|v(n‑1)|<sup>4</sup>,v(n‑2),v(n‑2)|v(n‑2)|<sup>2</sup>,v(n‑2)|v(n‑2)|<sup>4</sup>];   (5)e)v(n)通过公式1的“后失真器”处理得到输出信号z<sub>1</sub>(n);f)通过公式6计算预失真输出信号z(n)与“后失真器”输出信号z<sub>1</sub>(n)之间的估计误差e<sub>1</sub>(n):e<sub>1</sub>(n)=z(n)‑z<sub>1</sub>(n)   (6)g)通过公式7、8、9和10,计算得到DPD内核系数<img file="FDA0000824189280000027.GIF" wi="127" he="70" />π(n)=P(n‑1)u(n)   (7)k(n)=π(n)/(λ+u<sup>H</sup>(n)π(n))   (8)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>k</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><msub><mi>e</mi><mn>1</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000031.GIF" wi="1245" he="79" /></maths>P(n)=λ<sup>‑1</sup>P(n‑1)‑λ<sup>‑1</sup>k(n)u<sup>H</sup>(n)P(n‑1)   (10)上式中,P(n)为预失真器最佳相位系数向量,π(n)为预失真器最佳幅度系数向量,k(n)为增益向量,u<sup>H</sup>(n)是u(n)的复共轭转置,e<sub>1</sub><sup>*</sup>(n)为e<sub>1</sub>(n)复共轭,λ(0<λ<1)为遗忘因子;h)判断|e<sub>1</sub>(n)|大小,若|e<sub>1</sub>(n)|>0.01,则重复b到g;4)判断|e<sub>1</sub>(n)|大小,若0.01<|e<sub>1</sub>(n)|<0.05,切换至NLMS算法模块,同时把RLS迭代得到的系数送入所述模块,作为NLMS权系数的初始值,通过公式11和公式12进行进一步迭代处理,计算得到DPD内核系数<img file="FDA0000824189280000032.GIF" wi="132" he="79" />并且将<img file="FDA0000824189280000033.GIF" wi="103" he="76" />值复制给DPD系数C<sub>kq</sub>;<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;e</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>*</mo></msup><mi>u</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000034.GIF" wi="1182" he="79" /></maths>μ(n+1)=αμ(n)+βe<sub>1</sub><sup>2</sup>(n)   (12)上式中,e<sub>1</sub><sup>*</sup>(n)为e<sub>1</sub>(n)复共轭,μ为迭代步长,0<α<1,β>0;5)判断|e<sub>1</sub>(n)|大小,若|e<sub>1</sub>(n)|≤0.01,则|e(n)|=|e<sub>1</sub>(n)|,即|e(n)|≤0.01,故开关掷向2,此时的学习结构为直接学习结构,并用自适应算法1收敛后的权值向量作为自适应算法2权值向量的初始值,采用一种适合于硬件实现的自适应算法即LMS算法来计算DPD内核系数<img file="FDA0000824189280000035.GIF" wi="133" he="70" />a)用自适应算法1收敛后的权值向量作为自适应算法2权值向量的初始值<img file="FDA0000824189280000036.GIF" wi="133" he="68" />b)通过公式1计算得到x(n)经过预失真输出信号z(n);c)通过公式2、3和4计算得到信号z(n)经过功放的输出采样信号为y(n),y(n)经尺度变换后y(n)/G,其中G代表功放的增益,令z<sub>2</sub>=y(n)/G;d)通过公式13计算期望信号x(n)与功放估计输入z<sub>2</sub>之间的估计误差;e<sub>2</sub>(n)=z(n)‑z<sub>2</sub>(n)   (13)e)通过公式14计算得到DPD内核系数<img file="FDA0000824189280000037.GIF" wi="127" he="70" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;e</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>*</mo></msup><mi>u</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000824189280000038.GIF" wi="1301" he="78" /></maths>上式中,e<sub>2</sub><sup>*</sup>(n)为e<sub>2</sub>(n)复共轭,μ为迭代步长;f)判断|e<sub>2</sub>(n)|大小,若|e<sub>2</sub>(n)|>0.005,则重复b到e;若|e<sub>2</sub>(n)|≤0.005,DPD内核系数<img file="FDA0000824189280000041.GIF" wi="110" he="71" />值复制给DPD系数C<sub>kq</sub>。
地址 400041 重庆市九龙坡区科园三路7号2-4-4号