CN102970261A - 一种提高数字预失真处理速度的方法 - Google Patents
一种提高数字预失真处理速度的方法 Download PDFInfo
- Publication number
- CN102970261A CN102970261A CN2012104293483A CN201210429348A CN102970261A CN 102970261 A CN102970261 A CN 102970261A CN 2012104293483 A CN2012104293483 A CN 2012104293483A CN 201210429348 A CN201210429348 A CN 201210429348A CN 102970261 A CN102970261 A CN 102970261A
- Authority
- CN
- China
- Prior art keywords
- matrix
- distortion
- expression
- power amplifier
- sigma
- 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.)
- Pending
Links
Images
Landscapes
- Amplifiers (AREA)
Abstract
本发明公开了一种提高数字预失真处理速度的方法,包括以下步骤:(1)并行处理器件连续抓取功放输入信号和功放输出信号,并传给串行处理器件;(2)串行处理器件根据功放模型和LS算法原理计算自相关矩阵Ryy和互相关矩阵Ryz;(3)得到Ryy和互相关矩阵Ryz后,计算预失真系数矩阵;(4)并行处理器件接收预失真系数矩阵后,对输入信号进行预失真处理,再传给发射链路。本发明在实际应用中能够快速跟踪信号的变化,减少带外干扰。
Description
技术领域
本发明涉及功率放大器(功放)的数字预失真技术,特别涉及一种提高数字预失真处理速度的方法。
背景技术
随着无线通信的发展,信号带宽越来越宽,其包络起伏和峰均比也越来越大,这对宽带发信通道,尤其是混频器和功放提出了更高的线性要求。这主要是因为半导体器件在大信号下具有不可避免的非线性特性,当宽带发信机工作在非线性区时,会产生严重的互调分量,在工作频带内就可以造成信号间的相互干扰,在工作频带外产生不需要的频谱泄漏。
功率放大器是无线通信系统中必不可少的部分,非线性是其固有特性。然而,非线性会带来失真,导致系统误码率增加,产生邻道干扰。为了保证功率放大器的线性度与效率,人们已经提出了各种各样的方法对功放进行线性化处理,例如功率回退法(Back-off),负反馈法(Feedback),前馈法(Feed-forward),数字预失真法(Predistortion)等方法。其中,功率回退法是一种简单、可靠地线性化措施,它将输入功率降低使功率放大器工作在线性区域,但是,从经济方面来考虑,大功率的放大管价格昂贵,用牺牲功率放大器效率的方法来提高线性度是并不可取的;反馈法是采用外部电路来改善功率放大器的线性度的,但它并不适用于宽带信号;前馈法相比于反馈法,克服了延迟带来的影响,更加适宜用于宽带信号,但随着器件特性的变化其性能将变坏;数字预失真就是在功率放大器之前对输入信号先做预处理,它的传输特性刚好是放大器特性的反函数。由于数字预失真法具有电路形式简单、调整方便、效率高、造价低等优点,因此,它是目前使用得最多,最有发展前景的一种方法。
数字预失真技术的出现是线性化技术发展中非常重要的一步,其具有稳定性高,适用带宽高,功放成本低,散热量小,实现难度低等优点。同时,随着数字信号处理技术的高速发展,使得数字预失真系统的实现成本越来越低,而性能却不断得到提高。数字预失真技术由于自身具备的各种优势,正逐渐成为功率放大器的主要线性化技术。
在工程实践中,数字预失真处理的步骤通常是,在数字基带采集功放的输入和输出信号,基于预先设定的功放模型计算预失真系数,再将该系数用于预失真处理。数字预失真处理的一般硬件实现方法是:首先采集功放输入和输出信号,将数据传给DSP或软核等串行处理器件计算预失真系数,而使用预失真系数进行预失真处理在FPGA或专用芯片中进行,如申请号为CN201010558073.4的专利所述。这种实现方法在工程应用中存在的主要问题是,计算预失真系数的方法常采用最小均方算法(least square,i.e.LS),该算法中包括自相关矩阵计算步骤,该步骤的计算量非常大,在300MHz的DSP时钟下耗费时间达几十秒,占到整个DSP处理流程耗费时间的80%以上。而现代移动通信信号的变化是非常快的,随着用户的接入和退出,移动通信信号的带宽和幅度变化非常大,同时移动通信系统常采用跳频技术,即信号所处的频点也会在不同帧的时间段内发生跳变,通常变化的速度小于1秒。所有上述因素都对传统预失真技术实现方式的效果造成了影响,使得预失真技术在快变信号的条件下难以产生较好的效果,对相邻频段产生较大干扰,从而降低了移动通信系统的容量和覆盖范围。
发明内容
为了克服现有技术的上述缺点与不足,本发明的目的在于提供一种提高数字预失真处理速度的方法,在实际应用中能够快速跟踪信号的变化,减少带外干扰。
本发明的目的通过以下技术方案实现:
一种提高数字预失真处理速度的方法,包括以下步骤:
(1)并行处理器件连续抓取功放输入信号和功放输出信号,并传给串行处理器件;
(2)串行处理器件计算自相关矩阵Ryy和互相关矩阵Ryz;
设
其中,z(n)为功放输入信号,y(n)为功放输出信号,n=1,2,…,N,N为采集的功放输入信号和功放输出信号的总点数,akl为预失真系数,k为多项式阶数,l为记忆深度;K为最大多项式阶数,L为最大记忆深度;
Ryy为(K*L)*(K*L)的二维矩阵,设Ryy(p,q)为Ryy的第p行,第q列的元素,则
其中,y*(n+lp)表示y(n+lp)的复数共轭,lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
Ryz为(K*L)*1的一维矩阵,设Ryz(p)为Ryz的第p行的元素,则
其中lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
所述Ryy由以下方法计算:
(2-1)计算并保存矩阵polyMat={pn,k},其中
pn,k=y(n)*|y(n)|2k-1
(2-2)利用pn,k计算Ryy的上三角矩阵中的元素Ryy(p,q),
p=1,2,…,(K*L),对于第p行q=p,p+1,…,(K*L),lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
(2-3)由Ryy的上三角矩阵获得Ryy的下三角矩阵:
其中p=1,2,…,(K*L),对于第p行q=1,2,…,p-1;
(3)得到Ryy和互相关矩阵Ryz后,根据下式计算预失真系数矩阵:
将预失真系数矩阵传给并行处理器件;
(4)并行处理器件接收预失真系数矩阵后,对输入信号进行预失真处理,再传给发射链路。
与现有技术相比,本发明具有以下优点和有益效果:
本发明通过对自相关矩阵的矩阵元计算步骤和方法的简化,使得其计算量为传统方法的1/10以下,由于自相关矩阵的计算在DSP所执行的预失真系数计算步骤所耗费的时间中所占比重为80%以上,并且是周期性的运行,所以该发明能大大提高了预失真系数的更新速度,使得预失真系数的更新周期缩短到秒级时间单位,同时也减少了DSP的功耗,在实际系统的应用中,由于经过功放的信号随着用户的接入和退出,其幅度和带宽的变化非常快,自相关矩阵计算速度的提高使得预失真系统能够快速跟踪信号的变化,从而减少工作频带内干扰和带外的泄漏,提高移动通信系统的容量和覆盖范围。
附图说明
图1为本发明用于实现数字预失真处理的整体硬件框图。
图2为本发明提高数字预失真处理速度的方法的流程图。
具体实施方式
下面结合实施例及附图,对本发明作进一步地详细说明,但本发明的实施方式不限于此。
实施例
图1为本实施例用于实现数字预失真处理的整体硬件框图,包括串行处理器件、并行处理器件、发射链路、反馈链路、功率放大器和天线。
如图2所示,本实施例的提高数字预失真处理速度的方法,包括以下步骤:
(1)并行处理器件连续抓取功放输入信号和功放输出信号,并传给串行处理器件;
(2)串行处理器件根据功放模型和LS算法原理计算自相关矩阵Ryy和互相关矩阵Ryz;
设
其中,z(n)为功放输入信号,y(n)为功放输出信号,n=1,2,…,N,N为采集的功放输入信号和功放输出信号的总点数,akl为预失真系数,k为多项式阶数,l为记忆深度;K为最大多项式阶数,L为最大记忆深度;
Ryy为(K*L)*(K*L)的二维矩阵,设Ryy(p,q)为Ryy的第p行,第q列的元素,则
其中,y*(n+lp)表示y(n+lp)的复数共轭,lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
Ryz为(K*L)*1的一维矩阵,设Ryz(p)为Ryz的第p行的元素,则
其中lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
所述Ryy由以下方法计算:
(2-1)计算并保存矩阵polyMat={pn,k},其中
pn,k=y(n)*|y(n)|2k-1
(2-2)利用pn,k计算Ryy的上三角矩阵中的元素Ryy(p,q),
p=1,2,…,(K*L),对于第p行q=p,p+1,…,(K*L),lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
(2-3)由Ryy的上三角矩阵获得Ryy的下三角矩阵:
其中p=1,2,…,(K*L),对于第p行q=1,2,…,p-1;
(3)得到Ryy和互相关矩阵Ryz后,根据下式计算预失真系数矩阵:
将预失真系数矩阵传给并行处理器件;
(4)并行处理器件接收预失真系数矩阵后,对输入信号进行预失真处理,再传给发射链路。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (1)
1.一种提高数字预失真处理速度的方法,其特征在于,包括以下步骤:
(1)并行处理器件连续抓取功放输入信号和功放输出信号,并传给串行处理器件;
(2)串行处理器件计算自相关矩阵Ryy和互相关矩阵Ryz;
设
其中,z(n)为功放输入信号,y(n)为功放输出信号,n=1,2,…,N,N为采集的功放输入信号和功放输出信号的总点数,akl为预失真系数,k为多项式阶数,l为记忆深度;K为最大多项式阶数,L为最大记忆深度;
Ryy为(K*L)*(K*L)的二维矩阵,设Ryy(p,q)为Ryy的第p行,第q列的元素,则
其中,y*(n+lp)表示y(n+lp)的复数共轭;lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
Ryz为(K*L)*1的一维矩阵,设Ryz(p)为Ryz的第p行的元素,则
其中,lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
所述Ryy由以下方法计算:
(2-1)计算并保存矩阵polyMat={pn,k},其中
pn,k=y(n)*|y(n)|2k-1
(2-2)利用pn,k计算Ryy的上三角矩阵中的元素Ryy(p,q),
p=1,2,…,(K*L),对于第p行q=p,p+1,…,(K*L);lp=[p/K],表示p/K的整数部分;kp=mod(p,K),表示p/K的余数;lq=[q/K],kq=mod(q,K);
(2-3)由Ryy的上三角矩阵获得Ryy的下三角矩阵:
其中p=1,2,…,(K*L),对于第p行q=1,2,…,p-1;
(3)得到Ryy和互相关矩阵Ryz后,根据下式计算预失真系数矩阵:
将预失真系数矩阵传给并行处理器件;
(4)并行处理器件接收预失真系数矩阵后,对输入信号进行预失真处理,再传给发射链路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104293483A CN102970261A (zh) | 2012-10-31 | 2012-10-31 | 一种提高数字预失真处理速度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012104293483A CN102970261A (zh) | 2012-10-31 | 2012-10-31 | 一种提高数字预失真处理速度的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102970261A true CN102970261A (zh) | 2013-03-13 |
Family
ID=47800148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012104293483A Pending CN102970261A (zh) | 2012-10-31 | 2012-10-31 | 一种提高数字预失真处理速度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102970261A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103618684A (zh) * | 2013-11-25 | 2014-03-05 | 京信通信系统(中国)有限公司 | 一种估计dpd系数的方法及装置 |
CN103685110A (zh) * | 2013-12-17 | 2014-03-26 | 京信通信系统(中国)有限公司 | 一种预失真处理的方法、系统及预失真系数运算器 |
CN103701738A (zh) * | 2013-12-09 | 2014-04-02 | 京信通信系统(中国)有限公司 | 一种预失真参数确定方法、设备及数字预失真处理系统 |
CN103957179A (zh) * | 2014-04-04 | 2014-07-30 | 京信通信系统(中国)有限公司 | Dpd实现方法和系统 |
CN106031027A (zh) * | 2014-03-10 | 2016-10-12 | 德州仪器公司 | 用于开关模式功率放大器的数字预失真的方法及设备 |
CN107241070A (zh) * | 2016-03-29 | 2017-10-10 | 大唐移动通信设备有限公司 | 一种确定dpd系数的方法、fpga及dpd处理系统 |
CN110598271A (zh) * | 2019-08-22 | 2019-12-20 | 中国电子科技集团公司第二十九研究所 | 一种基于fpga实现4辅助天线slc功能的系统及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101286963A (zh) * | 2008-05-30 | 2008-10-15 | 北京北方烽火科技有限公司 | 一种基于可编程器件的宽带自适应数字预失真引擎装置 |
CN101771383A (zh) * | 2008-12-31 | 2010-07-07 | 大唐移动通信设备有限公司 | 一种实现信号预失真处理的方法和装置 |
US7822146B2 (en) * | 2006-12-01 | 2010-10-26 | Texas Instruments Incorporated | System and method for digitally correcting a non-linear element |
CN101902258A (zh) * | 2009-05-31 | 2010-12-01 | 大唐移动通信设备有限公司 | 一种数字预失真处理参数求取的方法及装置 |
CN102780665A (zh) * | 2012-07-25 | 2012-11-14 | 华南理工大学 | 一种数字预失真处理方法 |
-
2012
- 2012-10-31 CN CN2012104293483A patent/CN102970261A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7822146B2 (en) * | 2006-12-01 | 2010-10-26 | Texas Instruments Incorporated | System and method for digitally correcting a non-linear element |
CN101286963A (zh) * | 2008-05-30 | 2008-10-15 | 北京北方烽火科技有限公司 | 一种基于可编程器件的宽带自适应数字预失真引擎装置 |
CN101771383A (zh) * | 2008-12-31 | 2010-07-07 | 大唐移动通信设备有限公司 | 一种实现信号预失真处理的方法和装置 |
CN101902258A (zh) * | 2009-05-31 | 2010-12-01 | 大唐移动通信设备有限公司 | 一种数字预失真处理参数求取的方法及装置 |
CN102780665A (zh) * | 2012-07-25 | 2012-11-14 | 华南理工大学 | 一种数字预失真处理方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103618684A (zh) * | 2013-11-25 | 2014-03-05 | 京信通信系统(中国)有限公司 | 一种估计dpd系数的方法及装置 |
CN103618684B (zh) * | 2013-11-25 | 2017-09-29 | 京信通信系统(中国)有限公司 | 一种估计dpd系数的方法及装置 |
CN103701738A (zh) * | 2013-12-09 | 2014-04-02 | 京信通信系统(中国)有限公司 | 一种预失真参数确定方法、设备及数字预失真处理系统 |
CN103685110A (zh) * | 2013-12-17 | 2014-03-26 | 京信通信系统(中国)有限公司 | 一种预失真处理的方法、系统及预失真系数运算器 |
CN103685110B (zh) * | 2013-12-17 | 2017-03-22 | 京信通信系统(中国)有限公司 | 一种预失真处理的方法、系统及预失真系数运算器 |
CN106031027A (zh) * | 2014-03-10 | 2016-10-12 | 德州仪器公司 | 用于开关模式功率放大器的数字预失真的方法及设备 |
CN106031027B (zh) * | 2014-03-10 | 2020-02-28 | 德州仪器公司 | 用于开关模式功率放大器的数字预失真的方法及设备 |
CN103957179A (zh) * | 2014-04-04 | 2014-07-30 | 京信通信系统(中国)有限公司 | Dpd实现方法和系统 |
CN103957179B (zh) * | 2014-04-04 | 2018-04-20 | 京信通信系统(中国)有限公司 | Dpd实现方法和系统 |
CN107241070A (zh) * | 2016-03-29 | 2017-10-10 | 大唐移动通信设备有限公司 | 一种确定dpd系数的方法、fpga及dpd处理系统 |
CN107241070B (zh) * | 2016-03-29 | 2020-07-10 | 大唐移动通信设备有限公司 | 一种确定dpd系数的方法、fpga及dpd处理系统 |
CN110598271A (zh) * | 2019-08-22 | 2019-12-20 | 中国电子科技集团公司第二十九研究所 | 一种基于fpga实现4辅助天线slc功能的系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102970261A (zh) | 一种提高数字预失真处理速度的方法 | |
US20230370937A1 (en) | Method and system for baseband predistortion linearization in multi-channel wideband communication systems | |
CN102893399B (zh) | 预失真校正方法、预失真校正装置、发射机及基站 | |
CN104796091B (zh) | 基于分段记忆多项式的功率放大器建模及数字预失真方法 | |
CN102195912B (zh) | 数字预失真处理设备和方法 | |
CN106506417A (zh) | 一种窄带反馈的数字预失真系统与方法 | |
CN102075469B (zh) | 用于数字预失真系统的信号延迟时间的估计方法 | |
CN102143114A (zh) | 使用峰值加窗降低峰均功率比的方法和装置 | |
CN102970262B (zh) | 一种提高数字预失真稳定性的方法 | |
CN103296978B (zh) | 一种数字预失真方法及装置 | |
CN103618684A (zh) | 一种估计dpd系数的方法及装置 | |
CN102780665A (zh) | 一种数字预失真处理方法 | |
CN108011598A (zh) | 一种数字预失真自适应处理方法 | |
CN107659273A (zh) | 毫米波宽带功率放大器的可重构数字预失真系统及方法 | |
CN107565914A (zh) | 共发双带包络跟踪功率放大器的数字预失真系统及其方法 | |
CN103179074A (zh) | 基于正交多项式的自适应预失真系统及方法 | |
CN102904846B (zh) | 一种适应快变信号的数字预失真处理方法 | |
CN110086438B (zh) | 一种针对无源多波束发射机的数字预失真系统及方法 | |
CN103701738A (zh) | 一种预失真参数确定方法、设备及数字预失真处理系统 | |
CN103957179B (zh) | Dpd实现方法和系统 | |
CN111131104B (zh) | 一种预失真处理方法、装置、存储介质和设备 | |
CN109039969A (zh) | 宽带数字预失真器的实现方法 | |
CN104270330B (zh) | 一种数字预失真系统中的数据采集方法 | |
CN102522957A (zh) | 一种改善射频功率放大器预失真性能的方法 | |
CN102378336A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130313 |