劉素娟 張仲侯

(北京工業(yè)大學(xué)微電子學(xué)院 北京 100124)

1 引言

模擬數(shù)字轉(zhuǎn)換器(Analog-to-Digital Converter,ADC)在現(xiàn)代電子系統(tǒng)中扮演著至關(guān)重要的角色，在各個(gè)領(lǐng)域都獲得了廣泛的應(yīng)用[1-3]，高速、高分辨率是當(dāng)今ADC的發(fā)展方向。然而，對(duì)于ADC來(lái)說(shuō)，采樣率和分辨率是一對(duì)相互制約的指標(biāo)，由于制造工藝的限制，單個(gè)ADC很難同時(shí)滿足高采樣率和高分辨率的要求。文獻(xiàn)[4]在1980年提出通過(guò)時(shí)間交織(time-interleaved)的方式對(duì)信號(hào)進(jìn)行采樣，即時(shí)間交織模數(shù)轉(zhuǎn)換器(Time-Interleaved Analog-to-Digital Converter, TIADC)。從而可以在現(xiàn)有制造工藝的條件下，成倍地提高采樣率。

由于器件制造工藝的非理想性，不可避免地會(huì)在通道之間引入失配誤差，這些誤差會(huì)降低TIADC的有效分辨率等動(dòng)態(tài)性能[5,6]。在過(guò)去的幾十年間，大量的學(xué)者對(duì)TIADC通道之間如增益失配誤差、采樣時(shí)刻失配誤差、頻率響應(yīng)失配誤差等線性失配誤差的估計(jì)和校正方法進(jìn)行了深入的研究[7-10]。近些年來(lái)，學(xué)者逐漸開(kāi)始對(duì)TIADC中的非線性失配誤差(nonlinear mismatch error)產(chǎn)生興趣，非線性誤差是由于模擬前端電路的非理想性、子ADC中的積分非線性以及微分非線性引入[11]，TIADC通道之間非線性的差異又會(huì)引起非線性失配誤差。

很多學(xué)者對(duì)TIADC非線性失配誤差進(jìn)行分析并進(jìn)行估計(jì)和校正[12-18]。文獻(xiàn)[12,13]提出了一種通道隨機(jī)化的校正策略，但是引入了額外的子ADC，增加了硬件資源消耗。文獻(xiàn)[14,15]提出了一種基于Volterra級(jí)數(shù)的TIADC行為級(jí)模型，對(duì)TIADC的動(dòng)態(tài)非線性失配誤差進(jìn)行建模并提出了一種適用于硬件實(shí)現(xiàn)的補(bǔ)償結(jié)構(gòu)。文獻(xiàn)[16,17]提出了TIADC靜態(tài)非線性失配誤差的自適應(yīng)校正算法，文獻(xiàn)[16,17]均利用多項(xiàng)式模型對(duì)TIADC靜態(tài)非線性失配誤差進(jìn)行建模，分別采用了帶寬有效利用(Band Width Efficiency, BWE)和哈達(dá)瑪調(diào)制的方法進(jìn)行校正，并采取了濾波最小二乘的方法進(jìn)行非線性失配誤差的估計(jì)。文獻(xiàn)[18]采用了簡(jiǎn)化的Volterra級(jí)數(shù)對(duì)雙通道TIADC靜態(tài)非線性失配誤差進(jìn)行建模并利用自相關(guān)函數(shù)和最小均方(Least Mean Square, LMS)迭代算法進(jìn)行了失配誤差的估計(jì)。

目前只有少量的文章對(duì)TIADC帶有記憶效應(yīng)的非線性失配誤差的估計(jì)和校正方法進(jìn)行研究。文獻(xiàn)[19]提出了聯(lián)合盲校正算法對(duì)雙通道TIADC的帶有記憶效應(yīng)的非線性失配誤差進(jìn)行校正，使用歸一化最小均方(normalized-LMS)算法對(duì)非線性失配誤差進(jìn)行估計(jì)和補(bǔ)償，但是該方法只適用于兩通道非線性失配誤差的校正，難以擴(kuò)展到多通道的情況；文獻(xiàn)[20]提出了多相非線性均衡器(polyphase NonLinear EQualizer, pNLEQ)，可以同時(shí)對(duì)線性失配誤差和高階非線性失配誤差進(jìn)行校正，但是因?yàn)樾Ｕ^(guò)程涉及復(fù)雜的復(fù)數(shù)運(yùn)算，會(huì)消耗大量的硬件資源，所以不便在實(shí)際的硬件結(jié)構(gòu)中實(shí)現(xiàn)。

本文針對(duì)M通道的TIADC帶有記憶效應(yīng)的非線性失配誤差提出了一種基于子通道重構(gòu)結(jié)構(gòu)的盲校正方法。相比于BWE和哈達(dá)瑪變換的方法：首先，不僅能對(duì)靜態(tài)非線性失配誤差進(jìn)行校正，還能對(duì)帶有記憶效應(yīng)的非線性失配誤差進(jìn)行校正；其次，不需要產(chǎn)生正弦調(diào)制序列，有更少的硬件資源消耗；最后，LMS算法模塊以及誤差重構(gòu)中乘法器模塊的工作頻率僅為T(mén)IADC總采樣頻率的1/M，降低了系統(tǒng)功耗。

2 系統(tǒng)建模

理想的TIADC結(jié)構(gòu)如圖1所示，本文假設(shè)只存在帶有記憶效應(yīng)的非線性失配誤差，偏置失配、增益失配以及時(shí)間失配誤差均已對(duì)其進(jìn)行補(bǔ)償。

圖1 TIADC結(jié)構(gòu)示意圖

2.1 帶有記憶效應(yīng)的非線性失配誤差的子ADC的離散模型

本文考慮一個(gè)帶限的模擬輸入信號(hào)xa(t)，其最高頻率為fmax。xa(t)的模擬頻譜(傅里葉變換)表示為Xa(jΩ)，其中Ω=2 πf表示模擬角頻率(rad/s)，f表示模擬頻率(Hz)。

采樣周期為T(mén)s理想ADC的離散時(shí)間輸出信號(hào)表示為

xd[n]的數(shù)字頻譜離散時(shí)間傅里葉變換(Discrete-Time Fourier Transform, DTFT)表示為Xd(ejω)，其中ω表示采樣率為fs下的數(shù)字角頻率(弧度/采樣周期，rad/sample)，關(guān)系為

因此，模擬頻譜和數(shù)字頻譜之間的關(guān)系可表示為[21]

圖2是帶有記憶效應(yīng)的非線性失配誤差的子ADC的離散模型，模擬輸入信號(hào)xa(t)首先經(jīng)過(guò)采樣得到了離散信號(hào)xd[n]，然后經(jīng)過(guò)P 階的離散時(shí)間Volterra級(jí)數(shù)(Discrete Time Volterra Series, DTVS)后得到了輸出信號(hào)xv[n]，因?yàn)楸疚膬H考慮非線性失配誤差，所以2≤p≤P，DTVS的表達(dá)式為

圖2 帶有記憶效應(yīng)的非線性失配誤差的子ADC的離散模型

2.2 帶有記憶效應(yīng)的非線性失配誤差的TIADC的離散模型

對(duì)于一個(gè)M通道的TIADC，每個(gè)子ADC都具有不同的非線性誤差，整個(gè)系統(tǒng)的離散等效模型如圖3所示。為了方便后續(xù)表示，2.1節(jié)中的模擬輸入信號(hào)xa(t)以及第m通道采樣后的值分別用x(t)以及xm[nM]表示。此外，不同于2.1節(jié)的單通道ADC，本節(jié)以及之后的內(nèi)容，TIADC的采樣率和采樣周期為fs以及Ts，而每個(gè)子ADC的采樣周期為T(mén)IADC的M倍，即MTs。

圖3 帶有記憶效應(yīng)的非線性失配誤差的TIADC的離散模型

由2.1節(jié)中的理論可知，第m個(gè)子ADC的降采樣序列xm[v]可以表示為

其傅里葉變換表示為

cm,p,q表示第m個(gè)子ADC的p階非線性項(xiàng)記憶長(zhǎng)度為q時(shí)的誤差系數(shù)。

ym[v]經(jīng)過(guò)升采樣并相加之后得到了TIADC的輸出y[n]，其表達(dá)式為

3 校正方法

基于圖3中的TIADC等效模型，本節(jié)提出了一種針對(duì)M通道中帶有記憶效應(yīng)的非線性失配誤差的自適應(yīng)盲校正方法。校正方法的基本原理是利用估計(jì)的非線性誤差系數(shù)和子通道重構(gòu)結(jié)構(gòu)(Sub-Channel Reconstruction, SCR)重構(gòu)誤差信號(hào)，然后從TIADC的輸出中減去誤差信號(hào)，采用了濾波降采樣最小二乘(Filtered-Down-sampled Least Mean Square, FDLMS)算法在輸入無(wú)關(guān)頻帶(Input Free Band, IFB)中得到誤差信號(hào)以此來(lái)估計(jì)非線性誤差系數(shù)。

3.1 校正結(jié)構(gòu)

誤差信號(hào)e[n]可以由輸入信號(hào)x[n]進(jìn)行完美重構(gòu)，但是實(shí)際情況中，輸入信號(hào)是未知的，可以用y[n]來(lái)近似x[n]進(jìn)行重構(gòu)[17]。重構(gòu)信號(hào)可以表示為

c?m表示誤差估計(jì)矩陣。

圖4展示了M通道TIADC的校正結(jié)構(gòu)，TIADC的輸出y[n]經(jīng)過(guò)M倍降采樣率的降采樣模塊，數(shù)據(jù)的頻率降到了TIADC總采樣率的1/M，也就是每個(gè)子ADC的采樣率fs/M。后續(xù)模塊的工作頻率也均為fs/M。降采樣模塊子通道輸出ym[v]通過(guò)子通道誤差重構(gòu)(Sub-Channel Reconstruction, SCR)模塊后，再經(jīng)過(guò)升采樣模塊以及一個(gè)單位的延遲得到了重構(gòu)后的誤差信號(hào)e?[n]。利用ym[v]進(jìn)行乘方操作依次得到2階項(xiàng)ym2[v]到最高階P階項(xiàng)ymP[v]，然后分別乘以對(duì)應(yīng)的系數(shù)cm,2,0～cm,P,0，再累加得到了一個(gè)基本誤差重構(gòu)單元(Basic Error Reconstruction Unit, BERU)，利用ym[v]得到Q-1個(gè)延時(shí)分量ym[v-1]～ym[v-Q+1]，這Q-1個(gè)延時(shí)分量再經(jīng)過(guò)Q-1個(gè)基本誤差重構(gòu)單元，最后將Q個(gè)基本誤差重構(gòu)單元的輸出合并得到重構(gòu)后的信號(hào)yRS_m[v]，其公式為

圖4 M通道TIADC帶有記憶效應(yīng)的非線性誤差校正原理結(jié)構(gòu)圖

3.2 估計(jì)結(jié)構(gòu)和整體框架

為了估計(jì)系數(shù)矩陣cm的值，考慮一個(gè)帶限輸入信號(hào)，帶寬范圍為[0, βπ]，其中0<β<1。若TIADC內(nèi)不存在任何誤差，那么在TIADC的輸出頻譜僅在[0,βπ]內(nèi)有信號(hào)能量，當(dāng)非線性誤差存在時(shí)，在[βπ,π]頻帶范圍內(nèi)將會(huì)產(chǎn)生誤差能量，但是沒(méi)有輸入信號(hào)的能量，這個(gè)部分稱為輸入無(wú)關(guān)頻帶(Input Free Band, IFB)或者失配頻帶(Mismatch Band)。圖5給出了含有2階非線性誤差的四通道TIADC的IFB示意圖，其中X(jΩ)和E(jΩ)分別表示x[n]和e[n]的連續(xù)時(shí)間傅里葉變換。因此在[βπ,π]范圍內(nèi)的這部分能量可以通過(guò)一個(gè)高通濾波器濾出并結(jié)合FDLMS算法對(duì)誤差系數(shù)矩陣cm進(jìn)行估計(jì)。誤差信號(hào)e[n]經(jīng)過(guò)高通濾波器f[n]濾波后的誤差濾波信號(hào)ε[n]表示為

圖5 含有2階非線性誤差的四通道TIADC輸入無(wú)關(guān)頻帶示意圖[7]

表示信號(hào)向量X和Y經(jīng)過(guò)高通濾波器之后的部分。圖6為基于FDLMS迭代算法的M通道TIADC的非線性誤差自適應(yīng)盲校正的結(jié)構(gòu)示意圖，D代表高通濾波器f[n]的延遲，ε[n]經(jīng)過(guò)另外一個(gè)降采樣模塊得到εm[v],m=0, 1, 2, ··· , M-1。LMS為最小均方算法模塊，估計(jì)誤差系數(shù)矩陣c?m中元素c?m(q,p-2)=c?m,p,q基于FDLMS迭代算法的迭代公式為

其中，μ為迭代步長(zhǎng)，T=-(M+D)/M。迭代式(31)通過(guò)不斷迭代減小ε[n]的值，當(dāng)式(31)收斂時(shí)，ε[n]的值被減到最小，從而估計(jì)出系數(shù)c?m,p,q的值，進(jìn)而得到估計(jì)的系數(shù)矩陣c?m。 圖6中采用了向量表達(dá)方式使得圖示更簡(jiǎn)潔，其中

圖6 基于FDLMS算法的M通道自適應(yīng)盲校正算法結(jié)構(gòu)示意圖

4 仿真結(jié)與實(shí)驗(yàn)

本節(jié)展示了所提出校正結(jié)構(gòu)的仿真結(jié)果。利用MATLAB Simulink建立了四通道、16 bit的TIADC模型，高通濾波器利用FADTOOL工具進(jìn)行設(shè)計(jì)，階數(shù)為40階，通帶左側(cè)歸一化頻率設(shè)置為0.8π。迭代步長(zhǎng)μ設(shè)置為0.001，估計(jì)系數(shù)矩陣c?m(m=0, 1, ··· , M-1)初始設(shè)為零矩陣，非線性最高階設(shè)置為3，記憶長(zhǎng)度為2，系數(shù)矩陣c設(shè)置為

4.1 多頻點(diǎn)實(shí)驗(yàn)

本實(shí)驗(yàn)選取了如下輸入信號(hào)

由DC到0.79π(歸一化頻率)均勻分布的正弦信號(hào)。圖7比較了校正前后的頻譜。從結(jié)果可以看出，本文所提出的校正方法有效地提高了TIADC系統(tǒng)的無(wú)雜散動(dòng)態(tài)范圍(Spurious Free Dynamic Range,SFDR)。校正之前由于非線性失配誤差，TIADC的SFDR只有51.1083 dB，校正之后，TIADC的SFDR增加到了92.3539 dB。圖8給出了每個(gè)通道分支對(duì)應(yīng)的誤差收斂曲線，其中黑色虛線表示真實(shí)值。

圖7 多頻點(diǎn)輸入信號(hào)校正效果示意圖

圖8 多頻點(diǎn)輸入信號(hào)情況下誤差收斂曲線

4.2 白噪聲實(shí)驗(yàn)

本實(shí)驗(yàn)選取了一個(gè)低通信號(hào)，仍然使用4.1節(jié)中的四通道TIADC。一個(gè)均值為0，方差為1的高斯白噪聲信號(hào)經(jīng)過(guò)一個(gè)截止頻率為0.78π的低通濾波器得到所需的低通信號(hào)。為了量化在通帶中的校正效果，本文計(jì)算了誤差向量幅度(Error Vector Magnitude, EVM)，其公式為[16]

其中，X(q)表示期望信號(hào)的頻譜，而Y(q)則表示校正信號(hào)的頻譜。Q是FFT的總點(diǎn)數(shù)(圖9中采用了點(diǎn)數(shù)為3.2×104的FFT)。圖9是校正前后的頻譜對(duì)比，可以看到，TIADC系統(tǒng)在輸入無(wú)關(guān)帶中的最大失真由-49.74 dB降低到了-89.85 dB，圖10給出了4個(gè)子通道所對(duì)應(yīng)的誤差收斂曲線，黑色虛線代表系數(shù)的真實(shí)值。TIADC校正前后的EVMdB分別為-47.77 dB和-88.11 dB。

圖9 低通白噪聲輸入信號(hào)校正效果示意圖

圖10 低通白噪聲輸入信號(hào)情況下誤差收斂曲線

4.3 計(jì)算復(fù)雜度

本節(jié)分析了所提出方法的計(jì)算復(fù)雜度。本方法對(duì)比其他論文主要的優(yōu)勢(shì)在于對(duì)硬件資源消耗更少，避免了正弦調(diào)制信號(hào)。此外，本文的帶有記憶效應(yīng)的校正結(jié)構(gòu)也適用于靜態(tài)非線性失配誤差校正，當(dāng)記憶長(zhǎng)度為1的時(shí)候，所提出的結(jié)構(gòu)即等效于靜態(tài)非線性模型。從圖6可以看到，一個(gè)BERU里面的乘法器數(shù)目為2(P-1)，因此每一個(gè)降采樣后的SCR模塊包含乘法器數(shù)目為2(P-1)Q，總的乘法器數(shù)目為2M(P-1)Q，但是這些乘法器工作頻率均為Fs/M，因此就等效于在Fs的工作頻率下使用了2(P-1)Q個(gè)乘法器[16]，加上高通濾波器的乘法器數(shù)目Nf/2+1，本文校正算法每個(gè)采樣周期使用的乘法器數(shù)目為2(P-1)Q+Nf/2+1。表1給出了本文和其他方法的對(duì)比。對(duì)比文獻(xiàn)[16]和文獻(xiàn)[17]，在只考慮靜態(tài)非線性的情況下，本文采用方法的乘法器較少，且不需要調(diào)制器，避免了生成正弦序列的復(fù)雜電路。

表1 計(jì)算復(fù)雜度的對(duì)比

圖11展示了3種方法所消耗的乘法器數(shù)量以及乘法器數(shù)目的比率。高通濾波器的階數(shù)為44，因此系數(shù)對(duì)稱的有限沖激響應(yīng)(Finite Impulse Response,FIR)濾波器乘法器數(shù)目為22，可以看到，隨著通道數(shù)目M以及非線性最高階數(shù)P的增加，BWE方法[16]以及哈達(dá)瑪矩陣調(diào)制方法[17]的計(jì)算復(fù)雜度不斷增加，而本文的方法僅隨著非線性誤差最高階P以及記憶長(zhǎng)度Q增長(zhǎng)，和通道數(shù)目M無(wú)關(guān)。在M=16,P=7, Q=2時(shí)3種方法消耗乘法器數(shù)目分別(BWE、哈達(dá)瑪、本文)為227, 125和47，消耗乘法器比率分別為6.31和3.57。本文所提結(jié)構(gòu)硬件消耗遠(yuǎn)小于另外兩種方法。其中比率公式為

圖11 計(jì)算復(fù)雜度對(duì)比

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

本文介紹了一種自適應(yīng)盲校正方法來(lái)校正M通道TIADC中帶有記憶效應(yīng)的非線性失配誤差。通過(guò)子通道重構(gòu)結(jié)構(gòu)(SCR)重構(gòu)非線性引起的誤差信號(hào)，并通過(guò)FDLMS算法估計(jì)了各通道的非線性失配誤差系數(shù)。本文所提出的方法首先不僅能校正靜態(tài)非線性失配誤差，也能校正帶有記憶效應(yīng)的非線性失配誤差，其次避免了調(diào)制器的使用，與基于正弦序列和哈達(dá)瑪變換的方法相比，硬件資源和功耗大量減少。