基于等效采樣的數(shù)字預(yù)失真方法

蔡順燕， 庹先國(guó)， 楊劍波

（1.成都理工大學(xué)地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，四川 成都 610059；2.成都師范學(xué)院，四川 成都 611130）

0 引 言

隨著現(xiàn)代無(wú)線通信技術(shù)的發(fā)展，頻帶資源變得越來(lái)越緊張，為此業(yè)界提出了WCDMA（寬帶碼分多址）、OFDM（正交頻分復(fù)用）等高頻譜利用率的傳輸和調(diào)制技術(shù)。采用這類技術(shù)所傳輸?shù)男盘?hào)具有寬頻帶、高峰均比等特點(diǎn)，這些特點(diǎn)決定了必須采用高線性度的射頻功率放大器，否則會(huì)產(chǎn)生嚴(yán)重的帶內(nèi)和帶外失真，增大通信系統(tǒng)誤碼率并干擾鄰近信道。目前還無(wú)法通過電路設(shè)計(jì)的方式從根本上解決功率放大器的非線性問題，于是人們提出了許多功放線性化技術(shù)，常用的功率放大器線性化技術(shù)有功率回退、前饋、預(yù)失真等，其中數(shù)字預(yù)失真技術(shù)[1-3]具有穩(wěn)定、高效、自適應(yīng)等優(yōu)勢(shì)，能達(dá)到中等程度的線性化，是目前廣泛應(yīng)用的一種功放線性化技術(shù)。

為獲取功放輸出信號(hào)的失真信息，數(shù)字預(yù)失真系統(tǒng)需將功放輸出信號(hào)反饋回來(lái)。由于功放的非線性會(huì)導(dǎo)致輸出信號(hào)頻譜擴(kuò)展，當(dāng)信號(hào)頻帶較寬時(shí)，經(jīng)過功放后反饋回來(lái)的信號(hào)帶寬較寬，ADC需要高采樣率，這對(duì)電路設(shè)計(jì)及硬件實(shí)現(xiàn)較高要求。為降低對(duì)反饋回路ADC采樣率的要求，文獻(xiàn)[4]提出了欠采樣預(yù)失真方法，使預(yù)失真反饋系統(tǒng)可用2倍基帶信號(hào)帶寬的采樣率對(duì)反饋信號(hào)進(jìn)行采樣，大大降低了預(yù)失真系統(tǒng)的實(shí)現(xiàn)難度。但隨著無(wú)線通信技術(shù)的發(fā)展，信號(hào)帶寬進(jìn)一步提高，使用欠采樣方法仍然需要較高的采樣速率。為此，文中提出一種等效采樣預(yù)失真方法，可進(jìn)一步降低要求，降低系統(tǒng)成本以及硬件設(shè)計(jì)和實(shí)現(xiàn)的難度。

1 等效采樣基本原理

在某些特定應(yīng)用中，當(dāng)信號(hào)頻帶較寬時(shí)，為降低對(duì)ADC采樣率的要求，人們提出等效采樣技術(shù)[5-7]，可利用低成本、低速率的ADC器件實(shí)現(xiàn)對(duì)高速周期信號(hào)的采樣。等效采樣利用信號(hào)的周期性，以增加采集時(shí)間為代價(jià)，通過對(duì)信號(hào)進(jìn)行多周期的采樣后重組原信號(hào)（即實(shí)時(shí)采樣得到的信號(hào)）。常用的等效采樣方法包括順序等效采樣和隨機(jī)等效采樣。

圖1 等效采樣基本原理

順序等效采樣的基本原理如圖1所示，每次采樣值依次對(duì)應(yīng)信號(hào)在一個(gè)周期內(nèi)某個(gè)采樣點(diǎn)的值，每次采樣相對(duì)于周期信號(hào)來(lái)說(shuō)都有相同的時(shí)間步進(jìn)，經(jīng)過多次采樣后，把采樣得到的數(shù)據(jù)按順序進(jìn)行重組，可以得到與實(shí)時(shí)采樣等效的采樣數(shù)據(jù)。圖1（a）中經(jīng)過20次等效采樣后的數(shù)據(jù)按順序重組后可得到圖1（b）中的實(shí)時(shí)采樣數(shù)據(jù)。圖1是每個(gè)周期內(nèi)只采集一個(gè)數(shù)據(jù)點(diǎn)的等效采樣實(shí)例，實(shí)際上也可以根據(jù)具體的工程需要，在一個(gè)周期中采集多個(gè)等間隔的數(shù)據(jù)點(diǎn)，從而縮短等效采樣時(shí)間（但所需采樣頻率大）。

2 等效采樣數(shù)字預(yù)失真方法

2.1 數(shù)字預(yù)失真系統(tǒng)結(jié)構(gòu)

典型的數(shù)字預(yù)失真系統(tǒng)結(jié)構(gòu)如圖2所示。數(shù)字域中的基帶IQ信號(hào)u（n）首先經(jīng)預(yù)失真器進(jìn)行預(yù)處理，預(yù)處理后的輸出信號(hào)x（n）經(jīng)過D/A轉(zhuǎn)換、正交上變頻后送入射頻功放，經(jīng)功率放大后的射頻信號(hào)s（t）由天線輻射出去。當(dāng)預(yù)失真器與功放的非線性特性相逆時(shí)，功放輸出信號(hào)為線性放大的射頻信號(hào)，所以預(yù)失真器的設(shè)計(jì)非常關(guān)鍵。為設(shè)計(jì)出比較理想的預(yù)失真器，需要將功放輸出信號(hào)s（t）反饋回來(lái)，反饋信號(hào)經(jīng)過耦合衰減、正交下變頻解調(diào)、A/D轉(zhuǎn)換后得到含有失真信息的基帶IQ信號(hào)y（n），基于輸出和反饋回來(lái)的的數(shù)字信號(hào)，通過相應(yīng)的算法就可設(shè)計(jì)出預(yù)失真器（在設(shè)計(jì)預(yù)失真器的過程中，通常將D/A變換、上變頻、功率放大、耦合衰減、下變頻及A/D變換整個(gè)通路看作一個(gè)整體的非線性系統(tǒng)進(jìn)行研究）。

圖2 數(shù)字預(yù)失真系統(tǒng)結(jié)構(gòu)

2.2 等效采樣預(yù)失真原理

實(shí)現(xiàn)數(shù)字預(yù)失真的重點(diǎn)是通過適當(dāng)?shù)乃惴ǐ@取預(yù)失真器的模型參數(shù)，為獲取預(yù)失真器的參數(shù)，常用的預(yù)失真學(xué)習(xí)結(jié)構(gòu)有直接學(xué)習(xí)結(jié)構(gòu)[8]、間接學(xué)習(xí)結(jié)構(gòu)[9]和基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)。其中，基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)通過使用兩次自適應(yīng)學(xué)習(xí)算法獲取預(yù)失真器模型參數(shù)，其優(yōu)點(diǎn)是能克服功放非線性系統(tǒng)加性噪聲的影響。考慮到等效采樣是利用多個(gè)周期采樣得到的數(shù)據(jù)重構(gòu)出實(shí)時(shí)采樣數(shù)據(jù)，不同周期采樣時(shí)刻可能會(huì)出現(xiàn)誤差，導(dǎo)致重構(gòu)出的采樣數(shù)據(jù)產(chǎn)生新的噪聲。為克服等效采樣引入噪聲帶來(lái)的影響，文中采用基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)。

結(jié)合等效采樣技術(shù)和基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)，文中提出的等效采樣預(yù)失真系統(tǒng)框圖如圖3所示。采用等效采樣技術(shù)可極大地降低對(duì)ADC采樣率的要求，但要求待采樣信號(hào)具有一定周期性，而實(shí)際的無(wú)線通信系統(tǒng)發(fā)送的是非周期信號(hào)，為使用等效采樣技術(shù)對(duì)反饋信號(hào)進(jìn)行采樣，采用周期性訓(xùn)練序列的方式來(lái)獲取失真器的參數(shù)。

圖3 等效采樣預(yù)失真系統(tǒng)框圖

文中提出的等效采樣預(yù)失真方法分設(shè)計(jì)和實(shí)時(shí)兩個(gè)階段。在設(shè)計(jì)階段，預(yù)失真器的輸出用周期性訓(xùn)練序列取代（開關(guān)S置向訓(xùn)練序列端），訓(xùn)練序列發(fā)生器產(chǎn)生的周期性訓(xùn)練序列經(jīng)過功放非線性系統(tǒng)后，在反饋回路中用等效采樣的方式采集回來(lái)，然后采用適當(dāng)?shù)乃惴ň涂稍O(shè)計(jì)出預(yù)失真器在預(yù)失真器實(shí)時(shí)工作階段，預(yù)失真器的輸出送到功放非線性系統(tǒng)（開關(guān)S置向預(yù)失真器端），若預(yù)失真器設(shè)計(jì)合理，則功放非線性系統(tǒng)的輸出相對(duì)于預(yù)失真器的輸入是線性放大的。

下面分析基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)的抗噪聲特性，基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)分兩個(gè)步驟：首先識(shí)別功放的非線性模型，然后通過識(shí)別到的模型得到逆模型。如圖3所示的系統(tǒng)框圖，x（n）表示D/A轉(zhuǎn)換前的數(shù)字信號(hào)，功放非線性系統(tǒng)包含D/A變換、上變頻、功率放大、耦合反饋、下變頻、A/D變換整個(gè)過程，功放非線性系統(tǒng)的輸出y（n）表示反饋信號(hào)進(jìn)行A/D變換后的數(shù)字信號(hào)（y（n）是一個(gè)關(guān)于 x（n）的非線性函數(shù)），n0（n）表示加性噪聲（包括A/D量化噪聲、外部干擾噪聲、等效采樣所引入的新噪聲），功放模型的輸出可以表示為

式中：X（n）——由輸入信號(hào)x（n）構(gòu)成的向量（ 具體形式與功放模型有關(guān)）；

w——功放模型參數(shù)向量；

wH——共軛轉(zhuǎn)置。

功放模型識(shí)別的代價(jià)函數(shù)可以表示為

其中

該代價(jià)函數(shù)是一個(gè)關(guān)于模型參數(shù)向量w的二次性能函數(shù)，故存在唯一的w0使該代價(jià)函數(shù)的值最小。對(duì)該代價(jià)函數(shù)求偏導(dǎo)并令偏導(dǎo)數(shù)為0，則可得功放模型參數(shù)的最佳解（即維納解）為

從式（2）、式（5）可以看出，加性噪聲n0（n）會(huì)影響代價(jià)函數(shù)的J（w）極值，但是不會(huì)影響最佳解，即采用該方法得到的功放系統(tǒng)模型參數(shù)不受加性噪聲影響，相應(yīng)地，采用逆模型算法所得到的預(yù)失真器模型參數(shù)也不會(huì)受此加性噪聲影響。

2.3 等效采樣的實(shí)現(xiàn)方案

順序等效采樣在具體實(shí)現(xiàn)時(shí)，控制多周期采樣中的采樣點(diǎn)與觸發(fā)點(diǎn)的時(shí)間間隔是關(guān)鍵，為此需要設(shè)計(jì)一個(gè)精確的延時(shí)單元，使下一次采樣相對(duì)于本次采樣延遲一個(gè)采樣點(diǎn)的時(shí)間。此外，在實(shí)際的預(yù)失真系統(tǒng)中，D/A輸出的信號(hào)經(jīng)過功放系統(tǒng)反饋回來(lái)后會(huì)存在延時(shí)，這也是在等效采樣時(shí)需要考慮的問題。若在預(yù)失真系統(tǒng)中按照常規(guī)的等效采樣方法來(lái)實(shí)現(xiàn)則會(huì)使設(shè)計(jì)復(fù)雜化，為此文中提出了一種簡(jiǎn)易的等效采樣實(shí)現(xiàn)方案。

等效采樣實(shí)現(xiàn)方案框圖如圖4所示。該方案包含3個(gè)步驟，首先，利用周期性的訓(xùn)練序列送到D/A，信號(hào)經(jīng)過功放系統(tǒng)后反饋至A/D，其中A/D的工作頻率為D/A工作頻率的1/（N+1），A/D采樣得到的數(shù)據(jù)按順序組合后即得到等效采樣的數(shù)據(jù)；然后，利用等效采樣得到的一組反饋數(shù)據(jù)序列y（n）（采樣點(diǎn)數(shù)至少為周期性訓(xùn)練序列一個(gè)周期點(diǎn)數(shù)的2倍）與訓(xùn)練序列x（n）采用相關(guān)算法進(jìn)行延時(shí)估計(jì)，利用延時(shí)估計(jì)得到的延時(shí)將兩組數(shù)據(jù)序列（ x（ n）和 y（ n））在時(shí)域上對(duì)齊；最后，使用基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)設(shè)計(jì)出預(yù)失真器。

圖4 等效采樣實(shí)現(xiàn)方案框圖

如圖4所示，設(shè)D/A的工作頻率為f，周期性訓(xùn)練序列一個(gè)周期中包含的采樣點(diǎn)數(shù)為N，所以輸出訓(xùn)練序列的周期為N/f。根據(jù)訓(xùn)練序列的周期性X（ N+n）=X（ n），其中 n=1，2，…，N，若 A/D 采樣頻率為f/（N+1）（即周期為（N+1）/f），則采樣周期與訓(xùn)練序列的周期相比，增加了一個(gè)采樣點(diǎn)的時(shí)間，這樣經(jīng)過連續(xù)N次采樣后可得到如下式所示的采樣序列值：

由于采樣周期為（N+1）/f，所以相對(duì)于周期為N/f的訓(xùn)練序列，每次得到的采樣點(diǎn)都后推一個(gè)采樣步進(jìn)（步進(jìn)為1/f），這樣經(jīng)過N次采樣后得到的采樣值（y1，y2，y3，…，yN）進(jìn)行順序組合后就可以得到周期信號(hào)完整的等時(shí)采樣周期。在數(shù)字預(yù)失真系統(tǒng)中，A/D和D/A的工作時(shí)鐘可以由統(tǒng)一的參考時(shí)鐘提供，使用FPGA等可編程器件可以非常精確、方便地完成A/D和D/A采樣時(shí)鐘的控制和分配。由以上分析得出的A/D工作頻率為D/A工作頻率的1/（N+1）。式（6）中是連續(xù)N個(gè)采樣點(diǎn)的示例，若連續(xù)采取M×N個(gè)采樣點(diǎn)，則順序組合后得到的采樣序列為連續(xù)M個(gè)周期的訓(xùn)練序列。

在預(yù)失真系統(tǒng)中，信號(hào)經(jīng)過功放系統(tǒng)后反饋回來(lái)存在時(shí)間延遲，所以，等效采樣得到的采樣序列與原周期性序列相比有一個(gè)固定的延時(shí)，在使用自適應(yīng)算法設(shè)計(jì)預(yù)失真器之前，需要將等效采樣得到的反饋信號(hào)與輸出信號(hào)在時(shí)間上對(duì)齊，這里采用相關(guān)算法對(duì)該延時(shí)進(jìn)行估計(jì)。此外，考慮到等效采樣的誤差問題，為降低采樣誤差，可以將多次采樣得到的多周期序列取算術(shù)平均。

3 實(shí)驗(yàn)與結(jié)論

為驗(yàn)證文中提出的等效采樣預(yù)失真方法的正確性，在半實(shí)物實(shí)驗(yàn)平臺(tái)上進(jìn)行驗(yàn)證實(shí)驗(yàn)。Matlab中產(chǎn)生的數(shù)字基帶信號(hào)通過以太網(wǎng)下載到矢量信號(hào)發(fā)生器E4438C中，完成對(duì)數(shù)字基帶IQ信號(hào)的D/A轉(zhuǎn)換和正交上變頻調(diào)制（本實(shí)驗(yàn)中射頻載波為2.14GHz），得到上變頻后的射頻信號(hào)；經(jīng)過射頻功率放大器后，其中一路信號(hào)經(jīng)耦合器后送入實(shí)時(shí)頻譜分析儀RSA3408A，完成對(duì)射頻信號(hào)的正交下變頻解調(diào)、基帶信號(hào)的A/D采樣，并將解調(diào)后所得的數(shù)字基帶IQ信號(hào)反饋到PC機(jī)中，PC機(jī)Matlab中完成數(shù)字預(yù)失真的處理。另一路射頻信號(hào)經(jīng)衰減器后直接送入頻譜分析儀E4447A，實(shí)時(shí)顯示輸出信號(hào)的頻譜。

實(shí)驗(yàn)采用單載波WCDMA信號(hào)作為測(cè)試信號(hào)（碼片速率為3.84Mc/s），采用滾降因子為0.22的根升余弦成型濾波器，WCDMA信號(hào)的采樣點(diǎn)數(shù)為20000個(gè)，Matlab將信號(hào)下載到E4438C中后將循環(huán)產(chǎn)生以這20000個(gè)采樣點(diǎn)為周期的信號(hào)（E4438C的D/A變換速率為51.2 MHz）。實(shí)驗(yàn)中RSA3408A充當(dāng)反饋通道的功能，其實(shí)時(shí)采樣率為51.2 MHz。實(shí)驗(yàn)中等效采樣速率設(shè)置為實(shí)時(shí)采樣率的1/20（即2.56 MHz），每個(gè)周期中可以采樣到1 000個(gè)數(shù)據(jù)點(diǎn)（由于訓(xùn)練序列的周期較長(zhǎng)，若每個(gè)周期只采樣一個(gè)數(shù)據(jù)，則RSA3408A需要采集大量的數(shù)據(jù)，所以采用在一個(gè)周期內(nèi)采集多個(gè)等間隔數(shù)據(jù)的方式）。實(shí)驗(yàn)使用RSA3408A采集的20個(gè)周期的數(shù)據(jù)，每周期中按照2.56 MHz的等效采樣速率抽取出1 000個(gè)采樣點(diǎn)，20個(gè)周期即可構(gòu)成一個(gè)完整的采樣周期。

得到的實(shí)時(shí)采樣與等效采樣數(shù)據(jù)的時(shí)域?qū)Ρ热鐖D5所示，可以看出等效采樣可以很好地重現(xiàn)實(shí)時(shí)采樣的數(shù)據(jù)，等效采樣和實(shí)時(shí)采樣得到的數(shù)據(jù)誤差較小。圖6是實(shí)時(shí)采樣和等效采樣信號(hào)的頻域?qū)Ρ葓D，可以看出，等效采樣比實(shí)時(shí)采樣明顯多了噪聲信號(hào)，可知這些噪聲是由等效采樣引入的。

圖5 等效采樣與實(shí)時(shí)采樣數(shù)據(jù)時(shí)域?qū)Ρ?/p>

圖6 等效采樣與實(shí)時(shí)采樣數(shù)據(jù)頻域?qū)Ρ?/p>

基于等效采樣和實(shí)時(shí)采樣分別獲得的20000個(gè)采樣點(diǎn)，在Matlab中采用RLS自適應(yīng)算法首先得到功放非線性模型參數(shù)，然后再得到預(yù)失真器的模型參數(shù)（實(shí)驗(yàn)中功放非線性失真模型和預(yù)失真器模型均采用階次為7、記憶深度都為3記憶多項(xiàng)式模式）。圖7為采用實(shí)時(shí)采樣和等效采樣方法得到的預(yù)失真效果對(duì)比，可以看出，采用等效采樣預(yù)失真方法可以達(dá)到與采用實(shí)時(shí)采樣預(yù)失真方法相當(dāng)?shù)木€性化效果，等效采樣重構(gòu)數(shù)據(jù)時(shí)產(chǎn)生的誤差是導(dǎo)致其預(yù)失真效果略差的主要原因。

圖7 不同方法預(yù)失真效果對(duì)比

4 結(jié)束語(yǔ)

為進(jìn)一步降低對(duì)反饋回路ADC采樣率的要求，文中提出了等效采樣預(yù)失真方法。首先送出周期性訓(xùn)練序列，反饋回路采用等效采樣方式得到反饋信號(hào)，采用基于模型識(shí)別的學(xué)習(xí)結(jié)構(gòu)獲取預(yù)失真器的模型參數(shù)。然后再利用得到的預(yù)失真器對(duì)真實(shí)發(fā)送的信號(hào)進(jìn)行預(yù)失真處理，從而在功放輸出端得到線性放大的射頻信號(hào)。

文中提出的等效采樣預(yù)失真方法可以大大降低對(duì)反饋回路ADC采樣率的要求，從而降低系統(tǒng)實(shí)現(xiàn)難度和系統(tǒng)成本，有利于功放線性化系統(tǒng)的小型化。經(jīng)過實(shí)驗(yàn)證明，等效采樣預(yù)失真方法正確可行，該預(yù)失真方法可用于小型化、低功耗場(chǎng)合。

[1]Padmanabhan S， Peng L， Andreas S.Optimizing the identification of digital predistorters for improved power amplifier linearization performance[J].IEEE Transactions on Circuits and Systems-II：Express Briefs，2014，61（ 9）：671-675.

[2]Hu X， Wang G， Wang Z C， et al.Predistortion linearization of an X-Band TWTA for communications applications[J].IEEE Transactions on Electron Devices，2011，58（ 6）：1768-1774.

[3]Deleu T，Dervin M，Kasai K，et al.Iterative Predistortion ofthe NonlinearSatellite Channel[J].IEEE Transactions on Communications，2014，62（ 8）：2916-2926.

[4]詹鵬，秦開宇，陳長(zhǎng)偉，等.欠采樣方法在預(yù)失真線性化技術(shù)中的應(yīng)用[J].計(jì)算機(jī)工程與應(yīng)用，2011，47（ 16）：62-64.

[5]車?yán)Y留兵.偽隨機(jī)等效采樣法在超寬帶接收機(jī)中的應(yīng)用[J].現(xiàn)代雷達(dá)，2014，36（ 7）：62-44.

[6]余慧敏.基于FPGA的低功耗等效數(shù)據(jù)采樣裝置[J].儀表技術(shù)與傳感器，2013（ 6）：96-98.

[7]吳志剛，黃在朝，吳鵬.基于FPGA的等效采樣在光時(shí)域反射計(jì)中的應(yīng)用[J].光通信技術(shù)，2012（ 10）：44-46.

[8]Yang G H，Wang L L.One-step model extraction method for direct learning digital predistortion[J].Electronics Letters，2014，50（ 16）：1148-1150.

[9]Meenakshi R，Karun R，F(xiàn)adhel M，et al.Generalized rationalfunctions for reduced-complexity behavioral modeling and digital predistortion of broadband wireless transmitters[J].IEEE Transactions on Instrumentation and Measurement，2014，63（ 2）：485-498.