系研一下學期-移動通信技術講義

2008：p81-同步技術分ej

SNR

10.5947SNRsin2

同步誤差對系統(tǒng)性能影xn

NX

kej2kn Nk

,NGnN2為調制數(shù)據(jù)信號的能量，根據(jù)中心極限定理，當N 過假設信道最大多徑時延為L個樣點（L hn,lAlln2

2

n2Al為第l徑的能量，l

1,

為第l En,

h*n,

A

n2

fd ll

0 f hn,hn,lAllnl lNN

1Nl

k0OFDM系統(tǒng)中ICI分頻偏：頻率同步算法更快速準 定時偏差：接收端符號同步在無ISI的CP部－時變信道產(chǎn)生的接收信號經(jīng)過FFTr[k]

NN1y[n]ej2knN1k

N[k,k]d[k] [m,kL1

Aej2ml

Njj[m,

l[n]e

nmkNl0 第k2N

dNd

[m,k

[m,[m, E[m,k]2E[m,k]2E[m,k]*[m,kN1N

2n1n2fdej2mknnN2Nn10n2

J0[m,k]

1

N Nn

2nfdcos2nmk N2

0 f

子載波m對子載波k矩陣 h[n]ej2nN,是用子載波間隔歸一化的頻率N[m,k]1ej2nmkN[m,k E[m,k

sin2N2sin2mkN－定時偏差產(chǎn)生的若FFT窗偏離無ISI區(qū)域M個樣點，用x'ddLICI干擾信號：yICI[n] LlLM

A[n]x'[n

[n]

1M Nn0lLM

A[n]x'[nl]ej2kn 22M [k]E [k]Y*[k] NN

n0lLM1n每個子載波受到的均值干擾能量都相同，與子載波序號k同步誤差對OFDM系統(tǒng)性能影FFT窗定時在無ISI的CP區(qū)FFT窗定時在有ISI的CP區(qū)FFT窗定時在當前符號的數(shù)據(jù)區(qū)域：FFT窗內將包 同步算一般先作定時估計，再作頻偏估基于CP的同步算法：不需要額外的開基于CP和數(shù)據(jù)輔助的同步算xn

1XNNkN

L-hnlnlL-lL-ynlxnlL-l Snynej2nN和就是需要估計的信號時基于CP的同步算 r[k]skej2nN觀察2NG個連續(xù)抽樣點時，其中必包含一個完整的OFDMFFT窗的起始位置設為I,...,G1和IN,...,NG將2NG個觀察點作為一個向量rr1r2NGT22,m sEkr*ks

2ej

,m

2

,

對數(shù)似然函數(shù),定義為概率密度函數(shù)fr(k,,lgfr,只有屬于集合I和集合I'中的對應元素存在相關性，其它2N

lg

r(k),r(kN)

fr(k

kI'lgfr(k),r(kN)fr(k' fr(k)fr(kN) 因為fr(k)是對所有的2NL個點求乘積，故其值與k所以fr(k)的值也與k最大似然估計是要求出使,最大的和因此，省去fr(k)并不影響對和k,lgfr(k),r(kN fr(k)fr(kN) 分布，由CP的性質s(k)s(kNfr(k),r(kN)

swswexp

r(k)22Reej2r(k)r*(kN

r(kN)2wws s

2ej

2

2 r(k

r(kNwwss wwss

2

2

SNR是r(k)和r(kN),中的分母部分為2個一維 分fr(k)

exp

r(k)

2

2 wfr(kN)

exp

r(kN)2sw sw

2 ws,c Reej2r(k)r*(kN r(k)2r(kN)2wsG1 k c1c2cos2 clg12,c

21sw由于c1，c2為常數(shù)，對最大似 不產(chǎn)生影響sw,cos2

1 m

r(k)r(kN),m2

r(k)

r(kNk k,

因子與頻偏有是能量部分，與頻偏無關使似然函數(shù)最大化應該分2步來完成：max,max,max,max,?ML cos項為1，即2n2n為整數(shù)，可

1 由于三角函數(shù)的周期性，只能取其主值區(qū)間-,內的?ML12，只與有關，令,?ML最大化即可得 =argmax,?1

置的ML估計 12為對系統(tǒng)小數(shù)頻在估計過程中只用了N+G個信號可以利用多個OFDM符號數(shù)據(jù)來提高同步參數(shù)估計性接收信號r[kskej2kN發(fā)送信號sk的協(xié)方 M1m(NG)G式中，N r(k)r(kN*02

km(NGM1m(NG)G km(NG

2r(kN)2SNRSNR定時和頻偏的最佳估?argmaxN0?1 基于CP的同步算法（續(xù)集相關：一個OFDM塊中的一個采樣點的共軛與其相隔N從判定有用數(shù)據(jù)部分起始點來看，集相關算法檢測定時量度kMr(km(NG))r*(km(NG)NMrkm(NG)MMrkm(NG)Mrkm(NG)N2m0r(km(NG))r*(km(NG)Nm0rkm(N2m0rkm(NG)N2

0kN相關值的計算，灰色為CP部 基于CP的同步算法（續(xù)徑和頻偏時，接收信號的CP部分與它后移子載L-經(jīng)過多徑信道后的接收信號：rklksklwkL-lL-假設最大多徑時延擴展max小于保護間隔TG,將接收信號延時N個采樣點后減去當前L-lwdifkwkNwk，假設NG為CP長度，m是整由CP的性質，sklskNlmNlkmNl1k

kN

k,mN

kmN

kl

L-

l l當最

k2

kM

MGG

km2,M

從上面討論可知，當L1kNG-M1時，ravek應近似為在k＝NG-M2時，ravek變?yōu)橐粋€加大的值。于是，可以通過檢測這個突變可以通過計算ravek前后兩個樣值的

ravekSN=512，NG=100，QPSK調制

vSNR=5dB,基于CP的同步算法（續(xù)獲得較好的性能，但在多徑信道，由于CP中一C

自適應調整為無ISI污染的CP部n1avendifn,SNRSNR1了頻偏對符號定時的影響，且在一定范圍內隨頻移的數(shù)據(jù)輔助同步算s(n)

11

j2

,k,n0,1,...,NNkNk

Xk為在子載波k上調制的數(shù)據(jù)符號h,tltll對收到的信號在采樣之后得到的基帶r(k)expj2f0dTsexpj2f0kTslkdTssknlw(kl式中，nld Ts,f0是頻率偏移，d是定時偏移，Ts是采樣周定義f0NTs為相對頻率偏移，i和f分別為子載波間的整數(shù)倍頻偏和小數(shù)倍頻連續(xù)發(fā)送2個OFDM符號塊，在AWGN下接收到2Nr(n)s(n)expj2nNw(n),k0,1,...,2N前N點信號經(jīng)過FFT變換后的第k個元

N

rnexpjNkn,k0,1,...,N 后半部分N點信號經(jīng)過FFT變換后的第n個元素可表

2N

N

rnexpjNknrnNexpjNkn,n0,1,...,Nn

rnNr(n)expj2R(2)R(1)expj2 如果不考慮噪聲，R(2)和R(1)只差一相位expj2，如果加上白 Y(1)R(1)W

,W為對應的頻域噪Y(2)R(2)W 根據(jù)概率論中條件概率知識，經(jīng)過推導得到頻偏的MLN

* Im

k N

k 由于指數(shù)函數(shù)的周期為2，這種算法只能估計出12子載波間隔的頻 法是縮短FFT長度，f0Tf,TNTs,假設初始頻偏不超過max，那么最小的初始子載波間隔fminmax,FFT長度就確定了。如果縮短后符號的平均方差因為能量減小而變大，捕獲范圍增大是以犧牲同步數(shù)據(jù)輔助同步算法（續(xù))

Nr*

,r為接收到的時域信

N

dmdmNP(dP(d)rdmN2 在系統(tǒng)存在頻偏情況下，第一個訓練符號前后半個符號中具有 頻偏可估計為：1P(d P(d?)表示主值區(qū)間(-,內的值，f為小數(shù)頻偏，i為整數(shù)頻?1P(d?), i2zz為整第二個訓練符號偶數(shù)子載波上傳送的數(shù)據(jù)與第一個訓練符號中偶數(shù)子載波具有某種特定的PN序列關系，該已知PN序列關系將用于用戶整數(shù)頻偏估計。第二個訓練符號奇數(shù)子載波上傳送任意數(shù)據(jù)。偶數(shù)子載波位置，第二個訓練符號與第一個訓練符號具有一樣的X2X序列關系，Cm 2m,m0,1,...,N2,選擇的PN序列Cm應具有較2X的自相關特性，同時具有較低的PAPR特整數(shù)頻偏將使子載波位置整體發(fā)生循環(huán)移位，不會產(chǎn)生ICI，因可在頻域進行檢C，在變換到頻域前必須進行小數(shù)倍頻偏補償。用序列epfnN對訓練符號進行小數(shù)頻偏補償，之后變換到頻域，且用R(1)和R(2)分別表示變換到頻域的第一和第二 由于訓練符號1和符號2子載波相對位置并不?m

2RR ,?m與Cm只存在一定的相差，即整數(shù)倍頻偏的大小，因此2RRm可通過尋找互相關函數(shù)的最大值得到i的估

N?C*

,2zzN8,...,0,N8,估計范圍 4,Ni

數(shù)據(jù)輔助同步算法（續(xù)符號定時和小數(shù)頻偏估計與Cox算法相 Xk

,k0,2,...,N XN2,kN2 X2在符號同步和小數(shù)頻偏f估計完成后，首先需要用序列2k有Rejkk

k2z k2z

2

k2z 為初始相位，n為定時偏差， 與Cox算法一樣，小數(shù)頻偏補償后系統(tǒng)殘余整數(shù)頻偏為i2z整數(shù)頻偏的影響，除了移位2z外k Rk2的值近似等于Ck當存在定時偏差n0時，對所有信號對Rk，Rk2的值與Ck之間將引入一相位4n0 N。由于信道響應Hk的變化，也將引入一定的相位差在系統(tǒng)子載波數(shù)比較大的FM?max

C*

,2zzN8,...,0,N 0,2N2,估計范圍N4N

j40在用估計出的整數(shù)頻偏對信號糾正后仍有Yk Yk2Xk Xk2e0相鄰2個子載波上的信號的相位差與定時殘余偏移量n0有關，用去掉信號的影響，可求出 NYY* k 聯(lián)合CP和數(shù)據(jù)輔助的同步算 導頻信號自相關 設一個OFDM塊中的N個子載波上包含Np個導頻符s(n)

1

j2knNN 2 22Xk為在子載波k上調制的數(shù)據(jù)符號，其平均功率為2

p(n)

1Pexpj2knN x xPk為導頻符號，假設其平均功率為

2r(n)s(np(nw(n),w w假設時域信號s(n)是方差為N N） 當OFDM系統(tǒng)中Np 加CP后，s(ns(nNm(nm(nNn0,G1。Np如果很小，連續(xù)OFDM塊之間的相關性可忽略，只需要考慮因為噪聲是0均值過程，導頻符號在接收端確知，r(n)x為時變均值m(n)、方差為 x 1,kmck,m,kmN,k,G ,kmN,m,G x x

,SNR SNR 對應的值就是定時偏移的ML估計，即?

argmax經(jīng)過相關推導，CP1P ReG1r*(k)r(kN)

G

r(k)2

r(kN)2

2 k k 本質上是利用CP與其相隔N個采樣點的相關關 1Re r*(k)p(kN)ReG1r(k)r(kN)*p(k)

k 本質上是利用導頻符號的相關關系，由于導頻符號均勻分布，導致存多個相利用了CP的唯一性和P的尖對于高SNR，估計主要依賴于CP，反之考慮頻偏時：r(ns(np(n)expj2nN如果CP和P還取實部，由于復數(shù)乘以復指數(shù)函數(shù)?

SNR

1

(k)r(kN)

2

r(k

r(kN)2k

k

GP 1r*(k)p(kN)Pk

k

r(k)r(kN)

p(kSNR=9dB。參考定時估計為僅用CP的ML估計盲同步算不需要CP和訓練12子載波寬度的限盲同步算法（續(xù)K設噪聲干擾后的正弦信號：z(nxknw(nn1Kkxnexpjn,w(n)是方差為2 白噪聲 k，k和k分別是正弦波的幅度、頻率和初始相位令MNm1,KmNKm的目的是使采樣點的個數(shù) y(n)z(n),...,z(nm w(n)w(n),...,w(nm x(n)x(n),...,x(nm1) Aa(),...,a( y(n)Ax(n)w(n)，n1,...,kPEx(n)x(nk)k

,REy(n)y(n)H

APA

0 00PP 2 k根據(jù)協(xié)方差矩陣定義得到：Pk D

expjk因為A是范德蒙矩陣并且P0是正定的，可 k 2,k1,...,K, ... k 與k1m相對應的正交特征矩陣向量：Gg1gmk1 0 1 令

,

2 k R性質1AP0AHS性質2SAC，CP0AHSΛ1，且C1SH設A1Im10AA20Im1A,則A2A1D因為D是酉矩陣，A2A1D等效于空間旋轉對S作同樣處理，S1Im10S,S20Im1S ACADCSC1DCS S左邊乘以

H,再取逆，ΦSH

當使用SIT算法估計DIN是NN的單位陣，ei是IN的第i1 N12expj2mnN是NN的FFT矩N NDNh自變量是矢量表示對角陣，對角線上的元素DNhm,nhn對于對角陣，當對角線上元素是范德蒙矢量時僅用標量DN

DNfN第n個OFDM塊中有K個數(shù)據(jù)信息符s(n)s(nK),s(nK1),...,s(nKK經(jīng)過串并和插入虛子載波后得到N1矢量,N用i1...iNK作為插入虛子載波的索引，余下的K個索引i1...ik使

k k插入的虛子載波可用虛子載波插入矩陣表示：T

當TSC左乘s(n)時得到N1矢量，虛子載波的位置索引為用FH左乘 s(n)就相當于實現(xiàn)了IFFT變

Nk為了補償信道的時間彌散效應需要加CP，用 :

TIT N表示加CP操作，G是CP長度，信道階數(shù)LG，IGn是IN的最G行，加CP后符號塊長度為PN u(n) FH u(n是串行輸入的，經(jīng)過多徑信道h(t,)

L1(t)tl l 接收端的采樣周期為T，定義 0T為歸一化的頻偏，l l x(nexpj0nPDP0H0u(nH1u(n1)wnH0，H1分別是PP上、下三角Toeplitz矩陣T然 去掉CP，取CP矩陣表示為：RCP:0NGINy(n)

x(n)exp

nPG

s(n)CP

RCPH10NP,HRCPH0TCP,v(n)RCPwkH是循環(huán)矩陣 hklmodN,循環(huán)矩陣可由FFT矩陣kFHFH hNh(0),...,h(2N

N),h(2nN)

l

h(l)expj2lnN 表示子載波上的頻率響應，因為FFHI，利用H循環(huán)矩陣 y(n)expjnPGFH hTs(n 上式是去掉CP后的時域信號，如果00，不考慮噪聲，那DN0IN,FNy(n)DNhN這個等式證明了OFDM把卷積信道轉化為每個子載波上是平 信道如果0，F(xiàn) FH不是對角化矩陣，子載波間不再正交導

所以，在作FFT之前一定要消除00的影響，即要補償頻定義協(xié)方差矩陣

FH h DH

THFs(n)DH

和2I分別是數(shù)據(jù)信息和噪聲協(xié)方差矩陣

h N

K 實際中，集協(xié)方差矩陣（統(tǒng)計意義上）Ryy是由Nb個塊的采樣點的來近似?

y(n)y(n)00

APA

2I比較，A FH P h DH k2k1NΦ和D是相似矩陣，在OFDM系統(tǒng)expj1 D expj2k1N0將trDtrΦ代入上式，得到頻偏的閉式解 tr 0 k1expj2k1N當子載波上的頻率響應有一個為0，P0就會發(fā)生虧秩現(xiàn)令每個OFDM塊中連續(xù)的虛子載波位置以OFDM盲同步算法（續(xù)不存在頻偏時系統(tǒng)模N個子載波的OFDM系統(tǒng),Nc個子載波用作傳輸經(jīng)過調制的原始數(shù)假設子載波p0到p0Nc1用作數(shù)據(jù)子載波部分，p0是第一個數(shù)據(jù)子載 的索引，G為CP長度。設第k個OFDM符號頻域信息數(shù) S(k) S0(k),....,SP1T定義NN的IDFT u(N N , exp NN 1u(N

由IFFT變換進行多載波調制，產(chǎn)生的時域信號矢量x(k),....,

FS(kF為NN維部分IDFT矩0Fu

N

W(:,1:N up0(N1) 信號加入CP后，得到J1JNG維OFDMx(k) (k),...., (k),x(k),

(k),...,

N N

N FNG1:N,:S(k) 矢量x(k)中的每個元素或碼片經(jīng)過脈沖成型gtr(t)產(chǎn)生連續(xù)的時域信號送入信 Jx(tx(k tpkjT,T為碼片周期，設qpkJkqk和pqmod

x(txq(k)gtrtqTxq(k)為與xp(k)等價的傳輸碼片序列。接收信號r(t)xq(k)htqTw(t),htqTgtr(t)c(t)grxw(t)n(t)grx假設ht存在有限邊界0LT，并且GLJ，經(jīng)過采樣速率為1T采樣L的信號為：r(ir(t0iTxilht0lTw(it0為采樣相位，w(iw(t0iTLl設hlhtlT,則信道矢量hhLhL1h0T,由h構造一JL