一种增强带外DPD线性化的谱加权SAMP方法与流程

一种增强带外dpd线性化的谱加权samp方法
技术领域
1.本发明属于数字信号处理领域，尤其涉及一种增强带外dpd线性化的谱加权samp方法。


背景技术：

2.当今解决功率放大器(pa)固有线性度与低效率的最常用的方法是数字预失真技术(digital preditortion，dpd)。随着4g通信而兴起的一些技术，如载波聚合或者mimo，继而在5g中，不断增长的信号带宽和峰均功率比让预失真器在处理pa时遇到了困难。而dpd模型中补偿静态非线性和动态记忆效应所需要的参数数量非常大，增加了计算复杂性，这让dpd模型提取和自适应过程变得低效，导致了dpd估计参数时的过度拟合。然而，随着压缩感知(compressive sensing，cs)理论的发展与兴起，可以很好地解决自适应系统中的过拟合和复杂度高的问题，进而减少了所需dpd基函数的数目。
3.dpd模型降阶可以通过cs算法，如匹配追踪(mp)算法、正交匹配追踪(omp)算法、广义正交匹配追踪(gomp)算法、稀疏度自适应匹配追踪(samp)算法、最小二乘(ls)算法、稀疏贝叶斯学习(sbl)算法等。这些算法都能在不同程度上对功放的非线性进行有效改善，但为了无线通信系统的快速发展，开发新的更高效的预失真算法是促进当今行业发展的关键。


技术实现要素：

4.针对现有技术的不足之处，本发明提供一种增强带外dpd线性化的谱加权samp方法，通过基于谱加权策略的稀疏度自适应匹配追踪(sw
‑
samp)算法对功放模型进行降阶处理，进而补偿功放带外失真。
5.本发明的具体技术方案是：一种增强带外dpd线性化的谱加权samp方法，包括以下步骤：
6.步骤1：用一个测试信号来驱动gan器件的ab型超宽带功率放大器，获得pa的输入输出数据；
7.步骤2：采用gmp模型作为全功放行为模型，获得用于pa模型辨识的感知矩阵；
8.步骤3：对pa的输出信号在频域中进行带内陷波处理，得到频域中的由滤波基函数组成的感知矩阵；
9.步骤4：采用基于谱加权策略的稀疏度自适应匹配追踪(sw
‑
samp)算法对目标稀疏模型进行估计，得到一组最相关的基函数来描述pa的特性；
10.步骤4.1：初始化残差e0＝y
f
，初始化信号稀疏度k，迭代稀疏度l＝k，迭代次数m＝1，支撑集y
f
为o
×
1维观测向量，ψ
x,f
为n
×
o维滤波感知矩阵，输出o
×
1维滤波残差
11.步骤4.2：计算选择i中l个最大值，得到与滤波残差最匹配的原子，将这些值对应于ψ
x,f,{j}
的列序号j构成集合λ
k
；
12.步骤4.3：更新支撑集c
k
＝s
(m
‑
1)
∪λ
k
，ψ
x,f,m
＝{υ
x,f,j
|j∈c
k
}；
13.步骤4.4：求解的最小二乘解：
14.步骤4.5：在中选择绝对值最大的l项，记为
15.更新采样矩阵：ψ
x,f,m
中对应的l列记为ψ
x,f,ml
；
16.更新支撑集：选择初始支撑集中的l个元素记为s
ml
，则最终集为f＝s
ml
；
17.步骤4.6：计算滤波残差：
18.步骤4.7：如果残差e
new
＝0，则迭代停止，执行步骤4.8；若||e
new
||2≥||e
m
‑1||2，则l＝l k，返回步骤4.2继续迭代；若不符合前两个条件，则确定最终支撑集为f＝s
ml
，残差为e
m
＝e
new
，如果m≤n，则停止迭代执行步骤4.8，否则返回步骤4.2继续迭代；
19.步骤4.8：输出由最相关基函数组成的最终支撑集f和重构所得的信号。
20.步骤5：将输出得到的最优支撑集中的基函数按照其相关性排序，利用贝叶斯准则(bic)，确定重构所得系数参数的最佳数目，最后采用最小二乘法(ls)准则迭代来估计dpd系数。
21.进一步的，所述采用gmp模型作为全功放行为模型，获得的用于pa模型辨识的感知矩阵，表示如下：
[0022][0023]
式中是在n时刻的o
×
1维基函数。pa行为模型输出信号定义为：其中θ
x
为n
×
o维感知矩阵，包含基向量，为o
×
1维稀疏向量；
[0024]
进一步的，对输出信号进行带内陷波处理得到的频域内的n
×
o滤波基函数矩阵，表示如下：
[0025]
ψ
x,f
＝(υ
x,f
[0],
…
,υ
x,f
[k],
…
,υ
x,f
[n
‑
1])；
[0026]
其中为pa的输出信号y(n)在频域中进行带内陷波处理得到的输出信号，α是控制陷波水平的参数，y[k]是y(n)的离散傅里叶变换，如下式，k＝0,
…
,n
‑
1且n≥l；其中是的离散傅里叶变换，γ{
·
}是过滤算子，为频域中滤波得到的基函数。
[0027]
进一步的，用最小二乘法(ls)准则对功放的行为模型参数进行提取，表示如下：
[0028][0029]
式中是m
×
1(m＜o)维向量，μ为加权系数，g0是pa的线性增益。
[0030]
由上，本发明的增强带外dpd线性化的谱加权samp方法通过参考信号的带内陷波，大幅度地减小了有带外失真所导致的均方误差，因此在补偿相邻信道频谱再生的dpd分量的同时，还对dpd模型进行降阶处理，在高度拟合功放非线性系统的情况下显著地减少了模型系数，降低了系统的复杂度，减小固有的邻信道功率比(aclr)和归一化均方误差(nmse)，同时增强了dpd系统的稳定性和可靠性。
附图说明
[0031]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例的附图作简单地介绍。
[0032]
图1为本发明具体实施方式中基于sw
‑
samp算法的功放预失真系统的模型框图；
[0033]
图2为本发明具体实施方式中基于谱加权策略的稀疏度自适应匹配追踪算法的dpd模型降阶的方法的流程图。
具体实施方式
[0034]
下面结合附图对本发明做进一步说明，其作为本说明书的一部分，通过实施例来说明本发明的原理，本发明的其他方面、特征及其优点通过该详细说明将会变得一目了然。在所参照的附图中，不同的图中相同或相似的部件使用相同的附图标号来表示。
[0035]
本发明的基于sw
‑
samp算法的功放预失真系统如图1所示，其中离线运行模块中包括功率放大器(pa)和sw
‑
samp算法模块；功放逆模型的预失真模块中包括功率放大器(pa)、数字预失真器(dpd)、数模转换器(dac)、上变频模块、模数转换器(adc)、下变频模块、自适应参数提取模块与行为模型求逆模块。用一个测试信号来驱动pa，获得pa的输入输出数据，将测试信号通过所述的离线运行模块，此过程只进行一次；构建闭环dpd系统，将离线运行模块中选择的由最相关的基函数来定义前向路径中的闭环dpd函数，预失真器的输入端与信号源连接，输出端连接数模转换与上变频模块，经过上变频处理的信号与功放连接，由功放输出的信号经过g倍衰减器再进入模数转换与下变频模块，通过自适应参数提取和功放逆模型求取模块，把系数赋值给预失真器。
[0036]
本发明的增强带外dpd线性化的谱加权samp方法包括以下步骤：
[0037]
步骤一，用一个测试信号来驱动gan器件的ab型超宽带功率放大器，获得pa的输入输出数据；
[0038]
步骤二，通过matlab运行提出的基于谱加权策略的稀疏度自适应匹配追踪算法，以获得一组最相关的基函数；
[0039]
步骤三，构建功放逆模型的预失真系统，通过步骤二中得到的最相关基函数来定义前向路径中的预失真器；
[0040]
步骤四，将输出得到的最优支撑集中的基函数按照其相关性排序，利用贝叶斯准则(bic)，确定重构所得系数参数的最佳数目，最后采用最小二乘法(ls)准则迭代来估计dpd系数。
[0041]
基于谱加权策略的稀疏度自适应匹配追踪算法的具体步骤如图2所示。
[0042]
步骤1：初始化残差e0＝y
f
，初始化信号稀疏度k，迭代稀疏度l＝k，迭代次数m＝1，支撑集y
f
为o
×
1维观测向量，ψ
x,f
为n
×
o维滤波感知矩阵，输出o
×
1维滤波残差
[0043]
步骤2：计算选择i中l个最大值，得到与滤波残差最匹配的原子，将这些值对应于ψ
x,f,{j}
的列序号j构成集合λ
k
；
[0044]
步骤3：更新支撑集c
k
＝s
(m
‑
1)
∪λ
k
，ψ
x,f,m
＝{υ
x,f,j
|j∈c
k
}；
[0045]
步骤4：求解的最小二乘解：
[0046]
步骤5：在中选择绝对值最大的l项，记为
[0047]
更新采样矩阵：ψ
x,f,m
中对应的l列记为ψ
x,f,ml
；
[0048]
更新支撑集：选择初始支撑集中的l个元素记为s
ml
，则最终集为f＝s
ml
；
[0049]
步骤6：计算滤波残差：
[0050]
步骤7：如果残差e
new
＝0，则迭代停止，执行步骤8；若||e
new
||2≥||e
m
‑1||2，则l＝l k，返回步骤2继续迭代；若不符合前两个条件，则确定最终支撑集为f＝s
ml
，残差为e
m
＝e
new
，如果m≤n，则停止迭代执行步骤8，否则返回步骤2继续迭代；
[0051]
步骤8：输出由最相关基函数组成的最终支撑集f和重构所得的信号。
[0052]
其中采用gmp模型作为全功放行为模型，获得的用于pa模型辨识的感知矩阵，表示如下：
[0053][0054]
式中是在n时刻的o
×
1维基函数。pa行为模型输出信号定义为：其中θ
x
为n
×
o维感知矩阵，包含基向量，为o
×
1维稀疏向量；
[0055]
其中输出信号进行带内陷波处理得到的频域内的n
×
o滤波基函数矩阵，表示如下：
[0056]
ψ
x,f
＝(υ
x,f
[0],
…
,υ
x,f
[k],
…
,υ
x,f
[n
‑
1])；
[0057]
为pa的输出信号y(n)在频域中进行带内陷波处理得到的输出信号，α是控制陷波水平的参数，y[k]是y(n)的离散傅里叶变换，如下式，k＝0,
…
,n
‑
1且n≥l；其中是的离散傅里叶变换，γ{
·
}是过滤算子，为频域中滤波得到的基函数。
[0058]
本发明采用最小二乘法(ls)准则迭代来估计dpd系数，具体如下：
[0059]
用最小二乘法(ls)准则对功放的行为模型参数进行提取，表示如下：
[0060]
[0061]
式中是m
×
1(m＜o)维向量，μ为加权系数，g0是pa的线性增益。
[0062]
本发明利用基于谱加权策略的稀疏度自适应匹配追踪算法来降低数字预失真行为模型在前向路径中的计算复杂度，通过参考信号的带内陷波，构建了功放逆模型的预失真系统，大幅度地减小了有带外失真所导致的均方误差，因此不仅可以在补偿相邻信道频谱再生的dpd分量时，还可以对dpd模型进行降阶处理。该方法可以减小固有的邻信道功率比(aclr)和归一化均方误差(nmse)，同时降低了系统的复杂度并且提高了整个系统的稳定性。
[0063]
以上所述是本发明的优选实施方式而已，当然不能以此来限定本发明之权利范围，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和变动，这些改进和变动也视为本发明的保护范围。

技术特征：
1.一种增强带外dpd线性化的谱加权samp方法，其特征在于，包括以下步骤：步骤1：用一个测试信号来驱动gan器件的ab型超宽带功率放大器，获得pa的输入输出数据；步骤2：采用gmp模型作为全功放行为模型，获得用于pa模型辨识的感知矩阵；步骤3：对pa的输出信号在频域中进行带内陷波处理，得到频域中的由滤波基函数组成的感知矩阵；步骤4：采用基于谱加权策略的稀疏度自适应匹配追踪算法对目标稀疏模型进行估计，得到一组最相关的基函数来描述pa的特性；步骤5：将输出得到的最优支撑集中的基函数按照其相关性排序，利用贝叶斯准则，确定重构所得系数参数的最佳数目，最后采用最小二乘法准则迭代来估计dpd系数。2.根据权利要求1所述的增强带外dpd线性化的谱加权samp方法，其特征在于，所述采用gmp模型作为全功放行为模型，获得的用于pa模型辨识的感知矩阵，表示如下：式中是在n时刻的o
×
1维基函数，pa行为模型输出信号定义为：其中θ
x
为n
×
o维感知矩阵，包含基向量，为o
×
1维稀疏向量。3.根据权利要求1所述的增强带外dpd线性化的谱加权samp方法，其特征在于，用最小二乘法准则对功放的行为模型参数进行提取，表示如下：式中是m
×
1(m＜o)维向量，μ为加权系数，g0是pa的线性增益。4.根据权利要求1所述的增强带外dpd线性化的谱加权samp方法，其特征在于，步骤4包括以下步骤：步骤4.1：初始化残差e0＝y
f
，初始化信号稀疏度k，迭代稀疏度l＝k，迭代次数m＝1，支撑集y
f
为o
×
1维观测向量，ψ
x,f
为n
×
o维滤波感知矩阵，输出o
×
1维滤波残差步骤4.2：计算选择i中l个最大值，得到与滤波残差最匹配的原子，将这些值对应于ψ
x,f,{j}
的列序号j构成集合λ
k
；步骤4.3：更新支撑集c
k
＝s
(m
‑
1)
∪λ
k
，ψ
x,f,m
＝{υ
x,f,j
|j∈c
k
}；步骤4.4：求解的最小二乘解：步骤4.5：在中选择绝对值最大的l项，记为更新采样矩阵：ψ
x,f,m
中对应的l列记为ψ
x,f,ml
；更新支撑集：选择初始支撑集中的l个元素记为s
ml
，则最终集为f＝s
ml
；
步骤4.6：计算滤波残差：步骤4.7：如果残差e
new
＝0，则迭代停止，执行步骤4.8；若||e
new
||2≥||e
m
‑1||2，则l＝l k，返回步骤4.2继续迭代；若不符合前两个条件，则确定最终支撑集为f＝s
ml
，残差为e
m
＝e
new
，如果m≤n，则停止迭代执行步骤4.8，否则返回步骤4.2继续迭代；步骤4.8：输出由最相关基函数组成的最终支撑集f和重构所得的信号。
技术总结
本发明公开了一种增强带外DPD线性化的谱加权SAMP方法，通过构建一个离线运行模块，先用测试信号驱动功率放大器(PA)得到PA的输入、输出数据，再将PA的输出信号在频域中进行带内陷波处理，然后对稀疏度自适应匹配追踪(SAMP)算法进行改进，运用改进后的算法，即谱加权策略的稀疏度适应匹配追踪(SW


技术研发人员：高明明 边廷玥 刘超 南敬昌 杨会玉 梁琦 孙骥
受保护的技术使用者：辽宁工程技术大学
技术研发日：2021.03.25
技术公布日：2021/6/29