数字预失真分数时延估计与信号对齐算法及系统技术方案

本发明专利技术公开了一种数字预失真分数时延估计与信号对齐算法及系统，该算法包括如下步骤：首先，对DPD输入信号与PA反馈输出信号进行等速率采样，同时分别截取N点送存储器；对截取的数据做互相关运算；根据幅度相关函数法估计两路信号的整数倍采样间隔时延；然后，根据“三点二次”插值算法由幅度相关函数的最大值与次大值估计分数倍采样间隔时延；并将估出的整数倍采样间隔时延与分数倍采样间隔时延相加输出；最后，根据已知的信号时延，结合窗函数法与Nyquist插值滤波器设计“前插”滤波器，并作用于功率放大器的反馈信号，实现DPD输入信号与PA反馈输出信号的对齐。本发明专利技术DPD的信号时延估计时运算量低。

【技术实现步骤摘要】

本专利技术涉及通信信号数字处理
，具体涉及一种数字预失真分数时延估计与信号对齐算法及系统。
技术介绍
数字预失真(DPD)是功率放大器(PA)线性化校正的主流方法，其性能一直受困于信号对齐精度。由于，DPD系统需考虑PA带宽几个以上的失真，所以采样速率非常高，需要有性能良好的时延估计算法与信号插值对齐算法保障该系统的实现。目前，关于DPD系统中信号的时延估计算法主要有五种。一，迭代法时延估计算法。利用信号间存在一定线性关系，通过迭代方法实现DPD输出与PA反馈输出对齐。实践证明：该方法稳定性受迭代步长影响，延时估计精度差。二，延时锁定环(DLL)算法。利用环路鉴相器给出的误差信号控制环路中的压控振荡器VCO的输出时钟，实现PA反馈回路的采样脉冲的脉位变化达到信号对齐的目的。实践证明：该算法的实现复杂、收敛速度慢。三，相关函数法。又分幅度差相关函数法和幅度相关函数法两种。前者对输入信号和反馈信号的对应样点取幅度差，然后对差求互相关。再通过最大值估计信号间的整数(倍Ts)时延；与下面要讨论的幅度相关函数信号时延估计方法相比，增加了求插值过程。后者是最流行的信号同步方法，也是利用相关函数的最大值位置实现信号时延的估计，同样也实现整数(倍Ts)时延估计。实践中，为了采用该方法实现对分数时延的估计，需要对DPD的输入信号或者对PA的反馈输出信号进行高倍插值，再按插前间隔抽取，构成多路后，再与未进行插值操作的信号分别相关，然后取各路相关结果的峰值，再将这些峰值放在一起比较，取其中取最大的峰值，其对应的抽取信号即认为与另一相关信号是对齐的。该方法虽具有很多优点，但是缺点也相当明显：(1)不能从根本上消除时延的模糊问题；(2)运算量过大，DPD系统不能跟随PA参数的变化速度。为了保证DPD的性能，一般需要32或64倍甚至更高倍数的插值再抽取。四，基于小波分析时延估计方法。通过对信号的幅度相关函数做小波分析后，实现对信号时延的估计的。这种基于全局的寻找最优解过程运算量也不小。五，基于黄金分割法实现时延的估计。首先假设时延范围为[a，b]，采用黄金分割法在[a，b]中寻优。具体过程为：利用黄金分割方法，令利[a，b]中一点τ1，并用其对升余弦补偿器的系数进行估计，再讲估计出系数的滤波器与PA的反馈输出卷积得到输出v(n)，再将v(n)与PA反馈输出信号相关运算，重复上述这个过程，直到|τk-τk-1|小于设定的一个小常数ε，此时即得到同步输出vk(n)。可见对DPD系统的时延范围有先验知识决定了系统的运算复杂性。其实质是一种穷举方法，且在这一系列的运算过程中，难免出现算法的数值计算问题。实现信号的对齐方法主要有三种方法：一，基于重采滤波器法。因速度问题或硬件消耗问题不适合于DPD系统。二，基于多项式近似法。利用了一个低阶的多项式在一定区域对插值滤波器系数进行估计，以适应信号的变化。其实质为一个时变滤波器，该方法所实现的滤波器的新系数的误差取决于插值算法与老系数的量化精度；三。基于FRROW滤波器法。该方法是将多项式系数作为滤波器系数与输入卷积后，再与时间偏差卷积实现信号对齐，运算过程中也需要对滤波器权值修正，存在的数值计算问题也会影响到DPD系统的性能。
技术实现思路
为解决上述问题，本专利技术提供了一种数字预失真分数时延估计与信号对齐算法及系统，DPD系统中存在采样模糊问题，在详细分析了基于DPD的PA信道模型后认为：对DPD的高速采样而言，可利用互相关函数在峰值点附近具有二次函数特征，取其最大值与次大值三个点，采用本专利技术的分数倍采样间隔时延估计算法估计出小于采样间隔的信号时延，彻底解决离散互相关函数与连续互相关函数峰值间因采样造成的偏差问题；在获得了信号的精确时延信息后，根据多速率数字信号处理技术，采用窗函数方法对PA反馈输出进行“前插”递推，实现DPD系统中的信号的对齐。为实现上述目的，本专利技术采取的技术方案为：一种数字预失真分数时延估计与信号对齐算法，包括如下步骤：(1)对PA反馈信号y(n)与DPD输入信号d(n)进行相同速率采样，并同时分别取N点(L＝2M≥N，M∈Z+)，将其暂存到双口RAM中；(2)采样2L点FFT计算y(n)与d(n)的互相关；并提取r(n)中取最大值与两个次大值；(3)按照本专利技术所述时延估计y(n)与d(n)之间的分数间隔时延τ；(4)按照本专利技术所述插值滤波器系数计算方法计算hw(n)；(5)按照实现信号前插。结束。具体为：S1、同时读取DPD输入信号d(n)与PA反馈输出信号y(n)各N个点，并将其分别放到具有L＝2m，其中(m∈Z，L＞2N)个单元存储器中，d(n)与y(n)分别按照式(1)，d(n)=d1(n)+j·dQ(n)y(n)=y1(n)+j·yQ(n)---(1)]]>操作；S2、计算上述信号的互相关函数r(n)，取{|r(n)|本文档来自技高网...

【技术保护点】
数字预失真分数时延估计与信号对齐算法，其特征在于，包括如下步骤：首先，对DPD输入信号与PA反馈输出信号进行等速率采样，同时分别截取N点送存储器；对截取的数据做互相关运算；根据幅度相关函数法估计两路信号的整数倍采样间隔时延；然后，根据“三点二次”插值算法由幅度相关函数的最大值与次大值估计分数倍采样间隔时延；并将估出的整数倍采样间隔时延与分数倍采样间隔时延相加输出；最后，根据已知的信号时延，结合窗函数法与Nyquist插值滤波器设计“前插”滤波器，并作用于功率放大器的反馈信号，实现DPD输入信号与PA反馈输出信号的对齐；具体为：S1、同时读取DPD输入信号d(n)与PA反馈输出信号y(n)各N个点，并将其分别放到具有L＝2m，其中(m∈Z，L>2N)个单元存储器中，d(n)与y(n)分别按照式(1)，d(n)=dl(n)+j·dQ(n)y(n)=yl(n)+j·yQ(n)---(1)]]>操作；S2、计算上述信号的互相关函数r(n)，取{|r(n)|}的最大值记为f2，f2的左右两边的次大值分别记为f1、f3，本专利技术所述分数时延估计多项式如式(4)所示，τ=(f3-f1)4(f2-ηΔδ)-2·(f1+f3)·Ts---(4)]]>进行分数倍采样间隔估计；设信号的整数采样间隔时延为Γ，按照式(5)，Υ＝Γ+τ    (5)估计信号的总时延值Υ；S3、采用窗函数法计算“前”插对齐滤波器系数hw(n)；插值滤波器单位冲击响应按式(6)，hw(n)＝W(k)×sin[π·(Υ‑kTs)]/[π·(Υ‑kTs)]，k＝0，1，…M   (6)进行，其中W(n)=(1+Δδ/2)(n-M-12),(M-1)/2∈Z+(1+Δδ/2)(n-M2),M/2∈Z+]]>为窗函数，M为窗函数的长。...

【技术特征摘要】
1.数字预失真分数时延估计与信号对齐算法，其特征在于，包括如下步骤：
首先，对DPD输入信号与PA反馈输出信号进行等速率采样，同时分别截取N点送存储器；
对截取的数据做互相关运算；根据幅度相关函数法估计两路信号的整数倍采样间隔时延；
然后，根据“三点二次”插值算法由幅度相关函数的最大值与次大值估计分数倍采样间
隔时延；并将估出的整数倍采样间隔时延与分数倍采样间隔时延相加输出；
最后，根据已知的信号时延，结合窗函数法与Nyquist插值滤波器设计“前插”滤波器，

【专利技术属性】
技术研发人员：赵林军，张海林，刘龙伟，杨栋，武德斌，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61