本發(fā)明涉及一種基于單位基函數(shù)的2D-DPD迭代削減算法及應(yīng)用。

背景技術(shù)：
在無(wú)線通信系統(tǒng)中，發(fā)射機(jī)的線性度（Linearity）是評(píng)價(jià)發(fā)射機(jī)性能的重要指標(biāo)。而在所有用于提高發(fā)射機(jī)線性度的技術(shù)中，數(shù)字預(yù)失真（DPD,DigitalPredistortion）是當(dāng)前最受關(guān)注，也是使用最多的一種技術(shù)。數(shù)字預(yù)失真技術(shù)通過(guò)在數(shù)字基帶對(duì)發(fā)射機(jī)進(jìn)行測(cè)量，并建立發(fā)射機(jī)的行為級(jí)反向模型作為預(yù)失真器，由此對(duì)輸入信號(hào)進(jìn)行預(yù)矯正，用以抵消發(fā)射機(jī)的非線性，進(jìn)而獲取整體的高線性度。這其中，高精度的反向模型是整個(gè)數(shù)字預(yù)失真技術(shù)能夠生效的關(guān)鍵，當(dāng)前已經(jīng)有多種行為級(jí)模型可以對(duì)發(fā)射機(jī)建模，如：神經(jīng)網(wǎng)絡(luò)模型、頻譜選擇模型、Volterra模型等。這其中，Volterra模型最適合用于嵌入式系統(tǒng)的自適應(yīng)提取和實(shí)時(shí)實(shí)現(xiàn)。隨著無(wú)線通信技術(shù)的不斷發(fā)展，新一代的移動(dòng)通信系統(tǒng)要求能夠同時(shí)支持不同的標(biāo)準(zhǔn)在不同的頻段上進(jìn)行傳輸。然而，當(dāng)信號(hào)在不同頻帶同時(shí)傳輸時(shí)，其產(chǎn)生的行為效果要遠(yuǎn)比傳統(tǒng)的單頻段的射頻發(fā)射機(jī)更復(fù)雜，在每個(gè)頻帶上對(duì)信號(hào)產(chǎn)生干擾的不止有頻帶內(nèi)的信號(hào)本身的交調(diào)，還包含兩個(gè)頻帶信號(hào)之間的交調(diào)。因此，傳統(tǒng)的發(fā)射機(jī)行為模型并不能滿足對(duì)每個(gè)頻帶單獨(dú)進(jìn)行建模的需求。這預(yù)示著傳統(tǒng)數(shù)字預(yù)失真要經(jīng)過(guò)更仔細(xì)的擴(kuò)展才能在新型的協(xié)同式雙頻發(fā)射機(jī)中使用。目前，很多單頻發(fā)射機(jī)行為模型已經(jīng)被擴(kuò)展為可以用于對(duì)協(xié)同式雙頻發(fā)射機(jī)進(jìn)行建模的二維模型，這其中，最有潛力的當(dāng)屬二維數(shù)字預(yù)失真模型(2D-DPD),該模型通過(guò)在記憶多項(xiàng)式(MP,MemoryPolynomials)中引入兩個(gè)頻率之間的交調(diào)項(xiàng)實(shí)現(xiàn)對(duì)協(xié)同式雙頻發(fā)射機(jī)的精確建模。然而，由于考慮了兩個(gè)頻率之間的交調(diào)，其復(fù)雜度遠(yuǎn)高于傳統(tǒng)的記憶多項(xiàng)式模型。要將其在實(shí)際系統(tǒng)中進(jìn)行實(shí)現(xiàn)將消耗大量的硬件資源，且系統(tǒng)在高頻時(shí)種下的同步性很難保證。為了簡(jiǎn)化二維記憶多項(xiàng)式模型，有一種改進(jìn)的模型是改良的二維記憶多項(xiàng)式（twodimensionmodifiedmemorypolynomials,2D-MMP），這種模型通過(guò)引入一個(gè)自適應(yīng)的耦合因子，大大減少了所用基的數(shù)量。2D-MMP的缺點(diǎn)包括：自適應(yīng)選擇耦合因子的過(guò)程增加了模型提取的負(fù)擔(dān)；模型的基不是簡(jiǎn)單地乘積形式，因此模型本身也不再是簡(jiǎn)單地Volterra級(jí)數(shù)形式；最終模型中基的數(shù)量對(duì)于在線實(shí)現(xiàn)來(lái)說(shuō)仍然太高了。

技術(shù)實(shí)現(xiàn)要素：
為了克服上述現(xiàn)有技術(shù)的缺點(diǎn)，本發(fā)明的目的在于提供一種基于單位基函數(shù)的2D-DPD迭代削減算法及應(yīng)用，利用對(duì)發(fā)射機(jī)進(jìn)行測(cè)量采集到的信號(hào)反復(fù)訓(xùn)練模型,在此過(guò)程中可以觀察得到對(duì)系統(tǒng)影響較強(qiáng)的基，并削減掉對(duì)系統(tǒng)影響較弱的基，從而大大簡(jiǎn)化模型，該方法可以在線或離線完成，從而降低系統(tǒng)復(fù)雜度，節(jié)省硬件資源，降低系統(tǒng)延時(shí)，提高模型穩(wěn)定性，具有靈活度高、效果好的特點(diǎn)。為了實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案是：一種基于單位基函數(shù)的2D-DPD迭代削減算法，包括如下兩步：第一步：先對(duì)2D-DPD的基函數(shù)通過(guò)如下公式進(jìn)行歸一化得到單位基函數(shù)：（1）其中x1(n)和x2(n)表示發(fā)射機(jī)在高頻和低頻的輸入信號(hào)，y1(n)和y2(n)表示發(fā)射機(jī)在高頻和低頻的輸出信號(hào)，F(xiàn)m,k,j[]為模型的基，為對(duì)應(yīng)于Fm,k,j[y1(n),y2(n)]的標(biāo)準(zhǔn)差，為對(duì)應(yīng)于Fm,k,j[y2(n),y1(n)]的標(biāo)準(zhǔn)差，和為模型中需要估計(jì)的系數(shù)，稱之為對(duì)應(yīng)Fm,k,j[y1(n),y2(n)]或Fm,k,j[y2(n),y1(n)]的核，通過(guò)將基函數(shù)除以其標(biāo)準(zhǔn)差，對(duì)基函數(shù)進(jìn)行歸一化，此時(shí)和作為衡量該基重要性的標(biāo)準(zhǔn)，||表示幅度；第二步：通過(guò)迭代削減法逐一將對(duì)系統(tǒng)影響最小的模型削減掉直到模型復(fù)雜度符合要求。所述第二步采用迭代削減法，在每一次迭代的過(guò)程中利用第一步的公式對(duì)系統(tǒng)進(jìn)行建模，并削去最小的項(xiàng).該算法中用以評(píng)價(jià)基對(duì)系統(tǒng)影響強(qiáng)弱的標(biāo)準(zhǔn)是在使用歸一化基進(jìn)行建模時(shí)基對(duì)應(yīng)系數(shù)的模。所述基于單位基函數(shù)的2D-DPD迭代削減算法可以應(yīng)用于協(xié)同式雙頻發(fā)射機(jī)數(shù)字預(yù)失真系統(tǒng)中，包括如下過(guò)程：（1）在模型提取過(guò)程中選用單位基函數(shù)進(jìn)行訓(xùn)練；（2）在系統(tǒng)設(shè)計(jì)和發(fā)射機(jī)生產(chǎn)過(guò)程中，利用離線的迭代削減算法，減少資源分配，并保留一些用于保證模型兼容性的系數(shù)；（3）在基帶系統(tǒng)自適應(yīng)建模時(shí)，利用在線的迭代削減算法，減少基函數(shù)個(gè)數(shù)，提高建模算法穩(wěn)定性，減少系統(tǒng)延時(shí)?？傮w而言，該算法在發(fā)射機(jī)生產(chǎn)過(guò)程中離線執(zhí)行，在發(fā)射機(jī)基帶系統(tǒng)建模過(guò)程中在線執(zhí)行。與現(xiàn)有技術(shù)相比，本發(fā)明的算法具有靈活度高、效果好的特點(diǎn)。在使用時(shí)，該算法既可以在生產(chǎn)過(guò)程中離線執(zhí)行，也可以在系統(tǒng)建模過(guò)程中自適應(yīng)實(shí)現(xiàn)。本發(fā)明的算法若只是在生產(chǎn)過(guò)程中離線完成，可以保證最小的硬件資源消耗；但是因?yàn)槿鄙僮赃m應(yīng)訓(xùn)練環(huán)節(jié)，不能保證達(dá)到最好效果，因此應(yīng)當(dāng)接受一些冗余基函數(shù)，保證模型兼容性。本發(fā)明的算法若在系統(tǒng)建模過(guò)程中在線執(zhí)行，則可以最大程度降低在模型運(yùn)行過(guò)程中的延時(shí)消耗，提高模型穩(wěn)定性；但是因?yàn)橄到y(tǒng)硬件資源分配要照顧最大可能消耗情況，所以仍然需要占用大量硬件資源。因此，本發(fā)明的算法最好是首先在生產(chǎn)過(guò)程中離線執(zhí)行，選擇足夠保證模型兼容性的少量基函數(shù)，使系統(tǒng)的硬件資源消耗降低；在線建模時(shí)再次執(zhí)行，選取最少基函數(shù)個(gè)數(shù)。則可以同時(shí)獲得前面描述兩部分的全部?jī)?yōu)點(diǎn)。附圖說(shuō)明圖1為本發(fā)明基于單位基函數(shù)的二維數(shù)字預(yù)失真模型迭代削減算法流程圖。圖2為本發(fā)明具體實(shí)施方法流程圖。圖3為本發(fā)明算法測(cè)試平臺(tái)示意圖。圖4為本發(fā)明迭代削減過(guò)程中模型精度隨模型項(xiàng)數(shù)變化示意圖。圖5為本發(fā)明削減前后的模型用于DPD效果比較，其中A為高頻頻譜，B為低頻頻譜。具體實(shí)施方式下面結(jié)合附圖和實(shí)施例詳細(xì)說(shuō)明本發(fā)明的實(shí)施方式。如圖1所示，一種基于單位基函數(shù)的2D-DPD迭代削減算法，包括如下步驟：首先，進(jìn)行初始化，選定2D-DPD中的所有基函數(shù)；然后，利用單位基函數(shù)建模并計(jì)算模型誤差；如果模型誤差大于預(yù)設(shè)值或者基函數(shù)個(gè)數(shù)少于預(yù)設(shè)值，則裝載緩存的模型和基函數(shù)，算法結(jié)束；否則，將緩存模型和選定的基函數(shù)，削減掉對(duì)應(yīng)系數(shù)模最小的基函數(shù)，再次計(jì)算模型誤差，直到模型復(fù)雜度符合要求其中，單位基函數(shù)是通過(guò)對(duì)2D-DPD的基函數(shù)根據(jù)如下公式進(jìn)行歸一化得到的：（1）其中x1(n)和x2(n)表示發(fā)射機(jī)在高頻和低頻的輸入信號(hào)，y1(n)和y2(n)表示發(fā)射機(jī)在高頻和低頻的輸出信號(hào)，F(xiàn)m,k,j[]為模型的基，為對(duì)應(yīng)于Fm,k,j[y1(n),y2(n)]的標(biāo)準(zhǔn)差，為對(duì)應(yīng)于Fm,k,j[y2(n),y1(n)]的標(biāo)準(zhǔn)差，和為模型中需要估計(jì)的系數(shù)，稱之為對(duì)應(yīng)Fm,k,j[y1(n),y2(n)]或Fm,k,j[y2(n),y1(n)]的核，通過(guò)將基函數(shù)除以其標(biāo)準(zhǔn)差，對(duì)基函數(shù)進(jìn)行歸一化，此時(shí)和作為衡量該基重要性的標(biāo)準(zhǔn)，||表示幅度。具體實(shí)施該算法時(shí)可以參考附圖2所示流程圖，首先對(duì)模型進(jìn)行初始化，此例中假設(shè)初始模型為原始的2D-DPD模型，并將其格式化為向量Phi，其具體格式為：Phi={(0,0,0),…,(m,k,j),…,(M,N-1,0)}(2)這里，對(duì)應(yīng)2D-DPD模型的單位基函數(shù)：其中m表示基函數(shù)的記憶深度，k表示本頻帶的非線性階數(shù)，j表示交調(diào)頻帶的非線性階數(shù)。定義另一個(gè)向量Phi_Buf用來(lái)緩沖選定的基函數(shù)，首先對(duì)Phi進(jìn)行緩沖，即：Phi_Buf=Phi(4)隨后對(duì)選定的基函數(shù)Phi做關(guān)于(1)和(3)的建模，建模方法可以選用最小二乘法(LS,LeastSquareMethod)或共軛梯度法(CG,ConjugateGradientMethod)等任意用于解線性方程組的方法。所得的模型的基存儲(chǔ)在向量Kernel中。為了確定當(dāng)前模型的精度，下面計(jì)算模型輸出與期望輸出之間的歸一化均方誤差(NMSE,NormalizedMeanSquareError)，其定義式如下：其中P為用于檢測(cè)的采樣點(diǎn)的個(gè)數(shù)，xexp(n)為期望的模型輸出，由于DPD中為反向模型的訓(xùn)練，所以xexp(n)即為實(shí)際的輸入信號(hào)，xmod(n)為測(cè)試模型時(shí)得到的模型輸出；(5)計(jì)算得到的結(jié)果單位為dB。模型的NMSE計(jì)算結(jié)果被存儲(chǔ)在變量Error中。接下來(lái)對(duì)循環(huán)終止條件進(jìn)行判斷，該條件設(shè)置為：Error>Error_T或LENGTH(Phi)<LENGTH_T(6)其中LENGTH(Phi)表示向量Phi中元素的個(gè)數(shù)，Error_T和LENGTH_T為預(yù)設(shè)的兩個(gè)閾值。該條件既要求模型滿足精度要求，同時(shí)期望模型中基函數(shù)個(gè)數(shù)足夠少。若循環(huán)終止條件不滿足，循環(huán)繼續(xù)進(jìn)行，基函數(shù)向量Phi被緩存在Phi_Buf中，隨后在Kernel中找到幅度最小項(xiàng)，記下其坐標(biāo)n，并從Phi中刪除對(duì)應(yīng)項(xiàng)，即：Phi=Phi–Phi(n)(8)若循環(huán)條滿足，則讀取已經(jīng)緩存的模型，并結(jié)束程序。最終選取的基函數(shù)為：Phi_Buf=Phi(9)該算法的執(zhí)行效果可以通過(guò)儀器搭成的實(shí)驗(yàn)平臺(tái)進(jìn)行測(cè)試，測(cè)試平臺(tái)示意圖如附圖3所示。該平臺(tái)由兩臺(tái)信號(hào)發(fā)生器ESGE4438C，一臺(tái)頻譜儀PSAE4440，功率合成器和功率放大器（PA,PowerAmplifier）構(gòu)成；計(jì)算機(jī)（PC,PersonalComputer）在其中實(shí)現(xiàn)系統(tǒng)控制和信號(hào)處理的功能。雖然該算法是在儀器測(cè)試平臺(tái)上進(jìn)行測(cè)試的，但是它也可以在嵌入式系統(tǒng)中在線實(shí)現(xiàn)。整個(gè)基函數(shù)削減的過(guò)程如附圖4所示，從中可見(jiàn)，當(dāng)少量削減基函數(shù)時(shí)，模型的精度幾乎不受到影響；僅當(dāng)基函數(shù)數(shù)量很少的時(shí)候，模型的精度開(kāi)始對(duì)基函數(shù)數(shù)量很敏感，即是說(shuō)，最終剩余的基函數(shù)對(duì)于模型是很重要的。利用削減前后的模型進(jìn)行DPD的效果比較如附圖5所示，從中可見(jiàn)對(duì)基函數(shù)進(jìn)行削減前后的模型在DPD效果上幾乎沒(méi)有變化，然而基函數(shù)數(shù)量從24個(gè)削減到了9個(gè)。如前問(wèn)所述，當(dāng)該算法在線實(shí)現(xiàn)時(shí)，初始化向量Phi無(wú)需使用所有2D-DPD模型的基，而應(yīng)當(dāng)根據(jù)硬件設(shè)置恰當(dāng)選擇。