OFDM技術(shù)仿真(MATLAB代碼)..

1、第一章緒論1.1簡述OFDM是一種特殊的多載波傳輸方案，它可以被看作是一種調(diào)制技術(shù)，也可以被當作一種復(fù)用技術(shù)。多載波傳輸把數(shù)據(jù)流分解成若干子比特流，這樣每個子數(shù)據(jù)流將具有低得多的比特速率，用這樣的低比特率形成的低速率多狀態(tài)符號再去調(diào)制相應(yīng)的子載波，就構(gòu)成多個低速率符號并行發(fā)送的傳輸系統(tǒng)。正交頻分復(fù)用是對多載波調(diào)制(MCM,Multi-CarrierModulation)的一種改進。它的特點是各子載波相互正交，所以擴頻調(diào)制后的頻譜可以相互重疊，不但減小了子載波間的干擾,還大大提高了頻譜利用率。符號間干擾是多徑衰落信道寬帶傳輸?shù)闹饕獑栴},多載波調(diào)制技術(shù)包括正交頻分復(fù)用(OFDM)是解決這一難題中最

2、具前景的方法和技術(shù)。利用OFDM技術(shù)和IFFT方式的數(shù)字實現(xiàn)更適宜于多徑影響較為顯著的環(huán)境,如高速WLAN和數(shù)字視頻廣播DVB等。OFDM作為一種高效傳輸技術(shù)備受關(guān)注，并已成為第4代移動通信的核心技術(shù)。如果進行OFDM系統(tǒng)的研究，建立一個完整的OFDM系統(tǒng)是必要的。本文在簡要介紹了OFDM基本原理后，基于MATLAB構(gòu)建了一個完整的OFDM動態(tài)仿真系統(tǒng)。1.2OFDM基本原理概述1.2.1 OFDM的產(chǎn)生和發(fā)展OFDM的思想早在20世紀60年代就已經(jīng)提出，由于使用模擬濾波器實現(xiàn)起來的系統(tǒng)復(fù)雜度較高，所以一直沒有發(fā)展起來。在20世紀70年代，提出用離散傅里葉變換(DFT)實現(xiàn)多載波調(diào)制，為OFD

3、M的實用化奠定了理論基礎(chǔ)；從此以后，OFDM在移動通信中的應(yīng)用得到了迅猛的發(fā)展。OFDM系統(tǒng)收發(fā)機的典型框圖如圖1.1所示，發(fā)送端將被傳輸?shù)臄?shù)字信號轉(zhuǎn)換成子載波幅度和相位的映射，并進行離散傅里葉變換(IDFT)將數(shù)據(jù)的頻譜表達式變換到時域上。IFFT變換與IDFT變換的作用相同，只是有更高的計算效率，所以適用于所有的應(yīng)用系統(tǒng)。其中，上半部分對應(yīng)于發(fā)射機鏈路，下半部分對應(yīng)于接收機鏈路。由于FFT操作類似于IFFT，因此發(fā)射機和接收機可以使用同一硬件設(shè)備。當然，這種復(fù)雜性的節(jié)約則意味著接收發(fā)機不能同時進行發(fā)送和接收操作。接收端進行發(fā)送相反的操作，將射頻（RF，RadioFrequency）信號與基

4、帶信號進行混頻處理，并用FFT變換分解頻域信號。子載波幅度和相位被采集出來并轉(zhuǎn)換回數(shù)字信號。IFFT和FFT互為反變換，選擇適當?shù)淖儞Q將信號接收或發(fā)送。但信號獨立于系統(tǒng)時，F(xiàn)FT變換和IFFT變換可以被交替使用。1.2.2 串并變換數(shù)據(jù)傳輸?shù)牡湫托问绞谴袛?shù)據(jù)流，符號被連續(xù)傳輸，每一個數(shù)據(jù)符號的頻譜可占據(jù)整個可以利用的帶寬。但在并行數(shù)據(jù)傳輸系統(tǒng)中，許多符號同時傳輸，減少了那些在串行系統(tǒng)中出現(xiàn)的問題。在OFDM系統(tǒng)中，每個傳輸符號速率的大小大約在幾十bit/s到幾十kbit/s之間，所以必須進行串并變換，將輸入串行比特流轉(zhuǎn)換成為可以傳輸?shù)腛FDM符號。由于調(diào)試模式可以自適應(yīng)調(diào)節(jié)，所以每個子載波的

5、調(diào)制模式是可以變化的，因為而每個子載波可傳輸?shù)谋忍財?shù)也是可以變化的，所以串并變換需要分配給每個子載波數(shù)據(jù)段的長度是不一樣的。在接收端執(zhí)行相反的過程，從各個子載波出來的數(shù)據(jù)長度不一樣。在接收端執(zhí)行相反的過程，從各個子載波處來的數(shù)據(jù)被轉(zhuǎn)換回原來的串行數(shù)據(jù)。當一個OFDM符號在多徑無線信道中傳輸時，頻率選擇性衰落會導(dǎo)致某幾組子載波收到相當大的的衰減，從而引起比特錯誤。這些在信道頻率響應(yīng)的零點會造成在鄰近的子載波上發(fā)射的信息受到破壞，導(dǎo)致在每個符號中出現(xiàn)一連串的比特錯誤。與一大串錯誤連續(xù)出現(xiàn)的情況相比較，大多數(shù)前向糾錯編碼(FEC,ForwardErrorCorrection)在錯誤分布均與的情況下會

6、工作得更有效。所以，為了提高系統(tǒng)的性能，大多數(shù)系統(tǒng)采用數(shù)據(jù)加擾作為串并變換工作的一部分。這可以通過把每個連續(xù)的數(shù)據(jù)比特隨機地分配到各個子載波上來實現(xiàn)。在接收機端，進行一個對應(yīng)的逆過程解出信號。這樣，不僅可以還原出數(shù)據(jù)比特原來的順序，同時還可以分散由于信號衰落引起的連串的比特錯誤使其在時間上近似均勻分布。這種將比特錯誤位置的隨機化可以提高前向糾錯編碼(FEC)的性能，并且系統(tǒng)的總的性能也得到改善。1.2.3 子載波調(diào)制正交頻分復(fù)用(OFDM)技術(shù)就是在頻域內(nèi)將給定信道分成許多正交子信道，在每個子信道上使用一個子載波進行調(diào)制，并且各子載波并行傳輸。盡管總的信道是非平坦的，具有頻率選擇性，但是每個子

7、信道是相對平坦的，在每個子信道上進行的是窄帶傳輸，信號帶寬小于信道的相應(yīng)帶寬，因此大大消除了信號波形間的干擾。而且子信道的載波相互正交，一個OFDM符號包括多個經(jīng)過PSK調(diào)制或QAM調(diào)制的子載波的合成信號，每個子載波的頻譜相互重疊，從而又提高了頻譜利用率。用N表示子載波個數(shù)，T表示OFDM符號的持續(xù)時間，di(i=0,1,N-1)為分配給每個子信道的數(shù)據(jù)符號，為第i個子載波的載波頻率，從t=ts開始的OFDM符號的等效基帶信號可表示為(模擬信號表示式)：s(t)二才direct(t一ts一T/2)expj2兀i/T(t一ts)t=o(1-1)ttt+Tsss(t)二0tt十tT+ts或ss(t

8、)的實部和虛部分別對應(yīng)于OFDM符號的同相分量和正交分量，在實際系統(tǒng)中可分別與相應(yīng)子載波的余弦分量和正弦分量相乘，構(gòu)成最終的子信道。其相應(yīng)的數(shù)字表示式如下：令ts=0，采樣速率為N/T，則發(fā)送速率的第k(k=：0,1,N-1)個采樣表示為：s(k)二s(kT/N)二亍diexp(j2兀ik/N)(i_2)t=0顯然式上式恰好為IDFT的表達式，可知OFDM的調(diào)制和解調(diào)可以通過IDFT和DFT或(IFFT和FFT)來實現(xiàn)。如圖1.2所示，在一個OFDM符號內(nèi)包含四個載波的實例。其中，所有的子載波都具有相同的幅度和相位，但在實際應(yīng)用中，根據(jù)數(shù)據(jù)符號的調(diào)制方式，每個子載波都有相同的幅度和相位是不可能

9、的。從圖1.2可以看出每個子載波在一個OFDM符號周期內(nèi)都包含整數(shù)倍個周期，而且各個相鄰的子載波之間相差1個周期。這一特性可以用來解釋子載波之間的正交性，即：(1-3)1/TJTexp(j叫t)exp(j叫)dto如對式1-3中的第j個子載波進行調(diào)制，然后在時間長度T內(nèi)進行積分，即:bj1/T藝djexpj2兀i=0(t一t)dts丿=dj(1-4)=1/T寸expCj2兀j/T(t-t：diexpj2兀i/T(t一ts)t根據(jù)對式1-4可以看到，對第J個子載波進行解調(diào)可以恢復(fù)出期望的符號。而對其他載波來說，由于積分間隔內(nèi)，頻率差別(I-J)/T可以產(chǎn)生整數(shù)倍個周期，所以積分結(jié)果為零。這種正交

10、性還可以從頻率角度來解釋。根據(jù)式1-2，每個OFDM符號在其周期T內(nèi)包含多個非零子載波。因此其頻譜可以看作是周期為T的矩形脈沖的頻譜與一組位于各個子載波頻率上的8函數(shù)的卷積。矩形脈沖頻譜幅度值為sinc(/T)函數(shù)，這種函數(shù)的零點出現(xiàn)在頻率為1/T整數(shù)倍的位置上。圖1.2OFDM載波10.80.60.40.20-0.2-0.4-3000-6000-4000-200002000400060008000頻率圖1.3OFDM子載波頻譜這種現(xiàn)象可以參見圖1.3,圖中給出了相互覆蓋的各個子信道內(nèi)經(jīng)過矩形波形成型得到的符號的sine函數(shù)頻譜。在每個子載波頻率最大值處，所有其他子信道的頻譜值恰好為零。因為在

11、對OFDM符號進行解調(diào)過程中，需要計算這些點上所對應(yīng)的每個子載波頻率的最大值，所以可以從多個相互重疊的子信道中提取每一個子信道的符號，而不會受到其他子信道的干擾。從圖1.3可以看出，OFDM符號頻譜實際上可以滿足奈奎斯特準則，即多個子信道頻譜之間不存在相互干擾。因此這種一個子信道頻譜出現(xiàn)最大值而其他信道頻譜為零點的特點可以避免載波間的干擾(ICI)的出現(xiàn)。1.2.4 DFT的實現(xiàn)傅里葉變換將時域與頻域聯(lián)系在一起，傅里葉變換的形式有幾種，選擇哪種形式的傅里葉也變化由工作的具體環(huán)境決定。大多數(shù)信號處理使用DFT。DFT是常規(guī)變換的一種變化形式，信號在時域和頻域上均抽樣。由DFT的定義，時間上波形連

12、續(xù)重復(fù)，因此導(dǎo)致頻域上頻譜的連續(xù)重復(fù)?？焖俑道锶~變換(FFT)僅是計算應(yīng)用的一種快速數(shù)學方法，由于其高效性，使OFDM技術(shù)發(fā)展迅速。對于N比較大的系統(tǒng)來說，式1-1中的OFDM復(fù)等效基帶信號可以采用離散傅里葉逆變換(IDFT)方法來實現(xiàn)。為了敘述的簡潔，可以令式1-1中的ts=0，并且忽略矩形函數(shù)，對于信號s(t)以T/N的速率進行抽樣，即令t=kT/N(k=0,l,.,N-l),則得到：sk=s(kT/N)=diexp(j2兀ik/N)(0kN一1)(1-5)i=0可以看到Sk等效為對dj進行IDFT運算。同樣在接收端，為了恢復(fù)出原來的數(shù)據(jù)ki符號dj可以對sk進行逆變換，即DFT得到：di

13、二Lskexp(-j2兀ik/N)(0iN-1)(1-6)i=0根據(jù)以上分析可以看到，OFDM系統(tǒng)的調(diào)制和解調(diào)可以分別由IDFT和DFT來代替。通過N點的IDFT運算，把頻域數(shù)據(jù)符號q變換為時域數(shù)據(jù)符號Sk，經(jīng)過射頻載波調(diào)制之后，發(fā)送到無線信道中。其中每個IDFT輸出的數(shù)據(jù)符號$11111111保護間隔1111保護間隔IFFTFFT輸出；111FFTTFFT時間L丄1Jtel111111Ts11符號N-1111符號N111符號N-1F圖1.5OFDM符號形成過程符號的總長度為Ts=Tg+Tfft其中Ts為OFDM符號的總長度，Tg為抽樣的保護間隔長度，tff1為FFT變換產(chǎn)生的無保護間隔的OF

14、DM符號長度，則在接收端抽樣開始的時刻Tg應(yīng)滿足下式：ITimax的要求，則可以完全克服ISI的影響。同時由于OFDM延時副本內(nèi)所有包含的子載波的周期個數(shù)也為整數(shù)，時延信號就不會在解調(diào)過程中產(chǎn)生ICI。17第二章OFDM仿真結(jié)構(gòu)2.1OFDM傳輸系統(tǒng)一個完整的OFDM系統(tǒng)原理框圖如圖2.1所示，在發(fā)射端，輸入的高速比特流通過調(diào)制映射產(chǎn)生調(diào)制符號，經(jīng)過串并變換變成N條并行的低速子數(shù)據(jù)流，每N個并行數(shù)據(jù)構(gòu)成一個OFDM符號。插入導(dǎo)頻信號后經(jīng)快速傅立葉反變換(IFFT)對每個OFDMM符號的N個數(shù)據(jù)進行調(diào)制，變成的時域信號為：x(n)=IF=2X(m)ej2rann/Nn=0,1,N12-1)其中m

15、為頻域上的離散點，n為時域上的離散點，N為載波數(shù)目，為了在接收端有效抑ISI,通常在每一時域OFDM符號前要附加上長度為NG1個采樣的保護間隔(在OFDM中保護間隔一般選循環(huán)前綴CP)。加保護間隔后的信號可表示為公式(2-2)最后信號經(jīng)并/串變換及D/A轉(zhuǎn)換，由發(fā)送天線發(fā)送出去。=Jx(N+n)n=NCi，NCi+1,1、x(n)n=0,1,2,N1xCI(n)2-2)接收端將接收的信號進行處理，完成定時同步和載波同步。經(jīng)A/D轉(zhuǎn)換，串并轉(zhuǎn)換后的信號可表示為公式(2-3)：y(n)=x(n)*h(n)+z(n)+w(n)CICI2-3)然后,去CP后進行FFT解調(diào),同時進行信道估計(依據(jù)插入的

16、導(dǎo)頻信號),接著將信道估計值和FFT解調(diào)值一同送入檢測器進行相干檢測,檢測出每個子載波上的信息符號，最后通過反映射及信道譯碼恢復(fù)出原始比特流。移除cP，經(jīng)FFT變換后的信號可表示為式(2-4)：Y(m)=FFT(n)=1/N21y(n)ej2兀mn/Nn=0m=0,1,N1Y(m)=X(m)H(m)+Z(m)+W(m)m=0,1,N1(2-4)n圖2.1OFDM系統(tǒng)原理框圖其中H(m)為信道h(n)的傅立葉轉(zhuǎn)換，Z(m)為符號問干擾和載波問干擾z(n)的傅立葉轉(zhuǎn)換，W(m)是加性高斯白噪聲w(n)的傅立葉轉(zhuǎn)換。2.2OFDM仿真構(gòu)建OFDM系統(tǒng)編譯碼的數(shù)據(jù)處理量很大，利用矩陣對信息序列進行編碼

17、，譯碼等大量的運算都涉及到了矩陣運算，因此采用MATLAB來進行仿真。根據(jù)OFDM系統(tǒng)原理，下面以數(shù)字廣播電視(DVB)為例進行仿真。數(shù)字視頻廣播(DVB)通過兩種模式利用OFDM，這兩種模式的子載波個數(shù)分別為1705和6817，根據(jù)這兩種不同的子載波數(shù)量選擇所需要的FFT/IFFT的規(guī)模，因此這兩種模式也分別被稱為2K模式和8K模式。2K系統(tǒng)的子載波數(shù)量僅為8K的1/4，被稱為8K的簡化版本。本論文仿真的是2K模式的DVB，由于保護間隔也縮小到8K的1/4，因此在單頻網(wǎng)絡(luò)內(nèi)，2K系統(tǒng)處理時延擴展以及發(fā)射機之間的傳輸能力要下降。8K系統(tǒng)的FFT長度為896us,而保護間隔可以介于28us到22

18、4us之間。而2K系統(tǒng)的取值只為前者的1/4，圖2.4和圖2.5分別為DVB系統(tǒng)的發(fā)射機和接收機框圖。圖2.4DVB系統(tǒng)的發(fā)射機框圖在發(fā)射端，數(shù)據(jù)被分為若干組，每組內(nèi)包含188B，它們通過加擾碼和外碼R-S編碼，能夠在204B幀內(nèi)糾正8個錯誤字節(jié)。然后，對經(jīng)過編碼的比特由交織器在12B深度內(nèi)進行交織。并在按編碼效率為1/2，約束長度為7,生成多項式（171，133）的卷積碼進行編碼。通過打孔，編碼效率可以提高到2/3，3/4，5/6以及刀8。最后，經(jīng)卷積編碼的比特再經(jīng)過內(nèi)交織器的交織，被映射為4QAM。圖2.5DVB系統(tǒng)的接收框圖在接收端，要執(zhí)行相干QAM解調(diào)，就必須得到參考幅度、相位，這就要

19、求發(fā)送導(dǎo)頻子載波。對8K模式來說，每個OFDM符號內(nèi)包含768個導(dǎo)頻，剩余6048個子載波用于數(shù)據(jù)傳輸，對于2K模式來說，每個OFDM符號內(nèi)包含192導(dǎo)頻，剩余1512個子載波提供數(shù)據(jù)使用。導(dǎo)頻位置圖樣在每4個OFDM符號中重復(fù)一次，但是符號和符號之間是不同的。第三章OFDM仿真實現(xiàn)及結(jié)果3.1OFDM發(fā)送模塊一個從ts時刻開始的OFDM符號可以表示為:st)二0其中，Re(3-1)di為復(fù)合調(diào)制符號，率，標準的DVB(數(shù)字視頻廣播)表示如下:Ns為載波數(shù)，T為符號持續(xù)時間，fc為載波頻s(t)=Reexpj2兀fct藝工67見Cm,l,km,l,k(t)m二0l二0k二kmin(3-2)ex

20、p0j2兀kTU(tAl-Ts-68-m幾)(i+gg-m)-Ts其它其中:k為載波數(shù)；l為OFDM符號數(shù)；m為傳輸幀數(shù)；K為已傳輸載波數(shù)；Ts為符號持續(xù)時間；Tu為時延載波間隔時間;A為保護間隔；fc為射頻信號中心頻率；k為載波相對中心頻率，k=k(K+K.)/2；maxminm,0,k為復(fù)合符號表示m幁中第1個數(shù)據(jù)符號的第k個載波;m,1,k為復(fù)合符號表示m幁中第2個數(shù)據(jù)符號的第k個載波;m,63,k為復(fù)合符號表示m幁中第64個數(shù)據(jù)符號的第k個載波;在此采用傳輸速率為2K的數(shù)字廣播發(fā)送標準，這種模式在數(shù)字廣播電視（DTV）中被定義為移動接收標準。傳送的OFDM符號由很多幀結(jié)構(gòu)組成，每一幀持

21、續(xù)時間為Tf共包含68個OFDM符號。四個幀組成一個大幀結(jié)構(gòu).每一個符號是由2K模式下1705個子載波構(gòu)成并且其傳輸持續(xù)時間為Ts.在符號持續(xù)時間耳中有效符號持續(xù)時間為Tu,保護間隔時間為A。2K模式的具體參數(shù)參見表3.1：表3.12K模式OFDM參數(shù)參數(shù)2K模型載波數(shù)目K1705最小載波數(shù)Kmin0取大載波數(shù)Kman1704持續(xù)時間Tu224Ms載波間隔1/T“4464Hz最小載波K皿血與最大載波Kman（K-1）/兒間隔7.61MHz允許保護間隔時間A/Tu1/41/81/161/32有效符號持續(xù)時間Tu2048XT224Ms保護見個持續(xù)時間A512XT56Ms256XT28Ms128XT

22、14Ms64XT56Ms基本周期T7/64MsOFDM符號持續(xù)時間Ts=A+Tu2560XT280Ms2304XT252Ms2176XT238Ms2112XT231Ms從t=0到t=Ts對式3-2進行分析可以得到式3-3：maxs(t)=ReJej2兀fct冬maxcej2兀k(tA)/TuI0,0,k(3-3)(3-4)Ik=Kmink、k-(K.+K)/2minmax很明顯上式與反傅里葉變換（IDF）有相似之處:j2knq=l/NXeNqq=0有很多不同的FFT算法可以實現(xiàn)離散傅里葉變換(DFT)及離散傅里葉反變換(IDFT)這樣就很方便實際應(yīng)用中形成N個樣本xn使其對應(yīng)的每個符號有用部分

23、的持續(xù)時間為T。在時間保護間隔內(nèi)將后面NA/幾個樣點復(fù)制到前面，然后經(jīng)過集成上行轉(zhuǎn)換使信s(t)的中心頻率為fc。3.2OFDM符號的產(chǎn)生OFDM頻譜主要集中在fc附近，一種比較方便的實現(xiàn)方法是利用2-FFT和2-IFFT并且以T/2作為其基本周期。從表格2.1可以看出，OFDM符號持續(xù)時間為Tu，其為2048點的IFFT變換；因此要進行4096點的IFFT。圖3.2給出了OFDM符號產(chǎn)生方框圖，其中部分變量已標示出其用于Matlab代碼中以方便分析。T定義為信號的基本周期，既然模擬的是一個帶通信號就必須考慮其時間周期(1/Rs)其至少為載波頻率的兩倍。更一般地，用其整數(shù)倍Rs=40/T。這樣

24、一個關(guān)系式使載波頻率接近于902MHz，其描述如圖3.2所示。首先，隨機產(chǎn)生一個長為3412的二進制序列。然后，采用QAM映射，每兩位二進制比特映射成1土j中的一個。之后，進行4096點IFFT變換，先變?yōu)槟M值，再通過一個巴特沃斯低通濾波器，最后在發(fā)射端上變頻到射頻段以s(t)發(fā)送出去。CU3.3.2莫模擬生生PFD符符號在信源符號A中加入4906-1708=2391個零使其取樣為原來的兩倍并達到預(yù)期的中心頻率。從圖3.3和3.4可以看出這樣做的效果使得載波以T/2作為其時間周期。同時也注意到載波為離散時間的基帶信號，用發(fā)送濾波器產(chǎn)生一個連續(xù)時間信號g(t)作為復(fù)信號載波。其脈沖響應(yīng)和脈沖形

25、狀如圖3.5所示。500-600.20.460.61-6150100500時間(s)10-10000.20.40.60.811:2時間（s）*10CaiioGFFTDS1D5D.10圖3.5g(t)脈沖信號i_i|丨1|j02Q4:JJSCbI1.4I.-J1J8:f-:ir蛭皺趣扇功率諸豈度-打11:2ea：o1?：.i16ienhxcds冊餉謝屈圖3.3信號載波在B處時域響應(yīng)圖3.4在點B處載波信號的頻率響應(yīng)這個發(fā)射濾波器在時域和頻域的輸出顯示在圖3.7和圖3.8中。圖3.8的頻率響應(yīng)是周期的，這是因為離散時間信號在頻域是周期的，其頻譜帶寬取決于Rs。U(t)的周期是T/2,重建濾波器將會

26、有(T/2=18.286)-7.61=10.675MHz的過渡帶寬可以利用。如果用N點IFFT,過度帶寬只有(1/T=9.143)-7.61=1.533MHz,19因此為了避免混淆需要一個非常尖銳的滾降來較少重建濾波器的復(fù)雜程度。IIJJ-2斗6310I214圖3.6D/A濾波器響應(yīng)182D衰-:-5D-6070-8010UJ06M：It間杪)J.處正交廿覽直i02i0!繪FFT11F239刼率CHz疑訓(xùn)Bic演換的功率諾密崖051a2aalffi團?00.2D：015a2DISI：015殖耳(GHz)D，RHmE.暫囂餌吞60JJ16.S時間：松10X.圖3.7信號U在點C處時域響應(yīng)3.8信

27、號U在點D處頻域響應(yīng)33圖3.6給出了相對理想的DA濾波器器的頻率響應(yīng)。它是一個13階的截止頻率為1/T的巴特沃斯濾波器。該濾波器的時域和頻域響應(yīng)分別為如圖3.9和圖3.10。首先值得注意的是在濾波過程中在延遲產(chǎn)生在2x10-7附近，除了這一時刻其將按照預(yù)期進行濾波。這時從子載波853到1705其位置都為位于中心頻率(0Hz)的右邊，而1號子載波到852號在中心頻率(0Hz)以左4fc范圍內(nèi)。下一步要執(zhí)行多重雙正交單邊帶幅度調(diào)制uoft(t)。在這一調(diào)制中，存在一個同相信號m【(T)和一個正交信號mQ(T)其滿足式(3-5)：s(t)=叫(t)cos(2兀fct)+mQ(t)sin(2兀fct

28、)(3-5)式2-3可以展開為式2-6：k=Kmin(2兀KxRe(c)cos0,0,k+ft仝CT丿丿uIK+K.kmaxminKxIm(c)sin0,0,k(kminK+K.maxmin2+ft人TucTu丿丿(3-6)其中將同相信號和正交信號分別作為C和4-QAM的實部和虛部。對應(yīng)的IFFTm,l,k處理過程為：s(t)=uoftTcos(2兀ft)uofsin(2兀ft)(37)TcQc(3-7)信號s(t)的時域和頻域響應(yīng)如圖3.11和圖3.12。UOFT正交分量圖3.9信號在D點處的時域響應(yīng)UOFT經(jīng)FFT6040InES-E-:20000.511.522.533.54頻率糾垃K1

29、0s經(jīng)訕變換功率譜密度圖3.10信號在在點D處頻率響應(yīng)圖3.11信號s(t)在點E處時域響應(yīng)302D*IQDD.?i：F?蛭Welc茨換的功率沽密度25154:.；出005O?-10QI丄02D04Di方108II12I14HIEU10頻率(OR-15：-圖3.12信號s(t)在點E處頻率響應(yīng)3.3OFDM接收部分圖3.20是一個基本的OFDM接收機結(jié)構(gòu)。OFDM系統(tǒng)對時間和頻率偏移非常敏感。即使在理想的模擬環(huán)境下也要考慮濾波過程產(chǎn)生的延時。重建濾波和解調(diào)濾波所產(chǎn)生的延時td=64/Rs。這種延時影響了對信號的接收，從輸入（圖3.4）輸出（圖3.26）信號可以看出其細微的差別。處理好延時問題后

30、，接收部分的其他過程將很順利地進行。和發(fā)送部分一樣，我們再接收部分定義了各個階段的變量名稱（如圖3.20所示）以方便Matlab仿真對各個部分進行分析，其各部分仿真結(jié)果如圖3.20到圖3.29。r(t)Ja-hatfc圖3.20OFDM接收模擬r-tilde正交分皋圖3.21信號r-tilde在點F處的時域響應(yīng)可見發(fā)送端的信號s（t）經(jīng)信道傳輸?shù)竭_接收端r（t）。首先，在接收端下變頻到基帶形成信號r-tilde,其各部分分量響應(yīng)如圖3.21所示，再通過與發(fā)送端匹配的濾波器進行濾波量得到信號r-info,如圖3.23所示，可以發(fā)現(xiàn)于圖3.21的包絡(luò)相同這是因為發(fā)送信號為調(diào)頻信號而前面所做的工作只

31、不過是對原始發(fā)送信號的解調(diào)即恢復(fù)出來了信號UOFT,在不考慮信道的噪聲的影響所以其與發(fā)送端完全一樣。然后，通過A/D轉(zhuǎn)換抽樣成數(shù)字信號得到0FDM符號數(shù)據(jù)流r-data，如圖3.23所示。接著進行4096點FFT變換。最后，進行檢測與譯碼得到與發(fā)送端相同的二進制序列。醫(yī)0-E0Df-n同相廿擔Q20.4DEO.eI時間秒卜inf口正交分量1.2圖3.22信號r-tilde在點F處的頻域響應(yīng)-ZH/lap)姿陰凰K1EOim-50a20.4060.8時闔帖I3.23信號r-info在點G處的時域響應(yīng)101.2H10pinfoFFTED1111do:;20?:：qiiii5CsoD0:1P一，甲6

32、pir-dats同相分雖J.2J.4L.BL.t吋間（秒）r-dalsI亠00.511.522.533.5冊率甲9CK-50-IDGI：.12KlO5I.6I.l：m間（秒）圖3.25信號r-data在點H處的時域響應(yīng)圖3.24信號r-info在點G處的頻域響應(yīng)隨機輸入一組長為3412的二進制序列，從中截取前20個顯示出來，每個比特間隔為T/2,T為基本時間間隔，T=7/64s。將隨機輸入的二進制序列映射到QPSK星座圖上，結(jié)果如圖3.29所示。經(jīng)編碼和映射后的二進制序列變?yōu)橄鄳?yīng)的復(fù)數(shù)，再經(jīng)過IFFT、D/A轉(zhuǎn)換器和LPF后輸出的調(diào)制后的信號波形圖如圖3.11所示。圖3.11為經(jīng)理想信道傳輸?shù)?/p>

33、OFDM信號，圖3.29為經(jīng)高斯白噪聲信道傳輸?shù)腛FDM信號。從圖中可以看出在時域內(nèi)信號幅度值變化近似相等，但是在下圖中由于受噪聲的干擾信號幅度在短時間內(nèi)起伏變化很快。圖3.10和圖3.24為調(diào)制輸出信號的頻譜圖，在頻譜中可以發(fā)現(xiàn)，兩者的幅度值變化近似相同，只是受噪聲的信號圖3.24在接近零幅度處幅度起伏嚴重，從兩者的功率密度譜中對比發(fā)現(xiàn)理想信號的邊帶功率要比存在噪聲信號的邊帶功率大約低-36dB。圖3.26為接收信號經(jīng)A/D轉(zhuǎn)換后的功率譜密度，可以看出理想信號（圖3.4）和存在噪聲的信號兩者在幅度上有微小的變化，只是受噪聲影響后信號幅度在接近零點和接近最大值時旁邊有很多微小的起伏變化，這就導(dǎo)

34、致兩者的功率譜密度的值相差大約30dB。r-dataFFT1-2Di:r宣囂惡-eo邛口III1Iiiiirv-n:j呷-丨”I”丄111lJiiiiiiiii0.20.40.60.611.21.41.61.02頻!率Hz：K1O經(jīng)WmIE變換啟對率諸會康_-40-6045曰口12141E1曰圖3.26信號r-data在點H處的頻域響應(yīng)二圖326信號r-data在點H處的頻域響應(yīng)弓工工?1、:c-y11.510.501.5-0.500.511.5實部0.80.J0JI0.20.6-0.00.5a0.5實部圖3.28接收端a-hat星座圖圖3.27接收端info-h星座圖圖3.27為經(jīng)過檢測與譯

35、碼后的信號的星座圖，可以看出在理想情況下接收到的QPSK星座圖與發(fā)送端的QPSK星座圖十分接近，在理想情況下接收到的二進制序列與發(fā)送端的二進制序列相比在幅度上有微小的差距，相差大約為0.04%，可以認為接收到的數(shù)據(jù)是正確的；通過高斯白噪聲信道接收到的QPSK星座圖與發(fā)送端的QPSK星座圖相比，信號的星座圖展寬，但是基本能量還是集中在一點。第四章全文總結(jié)通過建立基于OFDM技術(shù)的DVB調(diào)制、解調(diào)模型，仿真分析了信號通過理想信道和加性高斯白噪聲兩種信道下的情況，可以證明OFDM系統(tǒng)具有內(nèi)在的適應(yīng)性，能夠適應(yīng)高斯信道，各子載波相互正交，所以擴頻調(diào)制后的頻譜可以相互重疊，不但減小了子載波間的干擾，還大

36、大提高了頻譜利用率。所以，作為第四代移動通信的主流技術(shù)，應(yīng)該廣泛的應(yīng)用到各種通信技術(shù)當中，提高數(shù)據(jù)傳輸速率和傳輸?shù)目煽啃?。附?：OFDM發(fā)送部分代碼clearall;closeall;%DVB傳送參數(shù)Tu=224e-6;%可利用OFDM符號時間T=Tu/2048;%原始基帶周期G=0;%允許保護時間間隔選擇1/4,1/8,1/16或1/32delta=G*Tu;%保護間隔持續(xù)時間Ts=delta+Tu;%整個OFDM持續(xù)時間Kmax=1705;%子載波數(shù)目Kmin=0;FS=4096;%IFFT/FFT長度q=10;%載波周期與原始基帶周期比fc=q*l/T;%載波頻率Rs=4*fc;%模擬

37、周期t=0:l/Rs:Tu;%數(shù)據(jù)產(chǎn)生程序(A)M=Kmax+l;rand(state,0);a=-l+2*round(rand(M,l).+i*(-l+2*round(rand(M,l).;A=length(a);info=zeros(FS,l);info(1:(A/2)=a(1:(A/2).;%補充零info(FS-(A/2)-l):FS)=a(A/2)+l):A).;%子載波產(chǎn)生程序(B)carriers=FS.*ifft(info,FS);tt=0:T/2:Tu;figure(1);subplot(211);stem(tt(1:20),real(carriers(1:20);subpl

38、ot(212);stem(tt(1:20),imag(carriers(1:20);figure(2);f=(2/T)*(1:(FS)/(FS);subplot(211);plot(f,abs(fft(carriers,FS)/FS);subplot(212);pwelch(carriers,2/T);%D/A轉(zhuǎn)換模擬L=length(carriers);chips=carriers.;zeros(2*q)-1,L);p=1/Rs:1/Rs:T/2;g=ones(length(p),1);%脈沖形成figure(3);stem(p,g);dummy=conv(g,chips(:);u=dumm

39、y(1:length(t);%(C)figure(4);subplot(211);plot(t(1:400),real(u(1:400);subplot(212);plot(t(1:400),imag(u(1:400);figure(5);ff=(Rs)*(1:(q*FS)/(q*FS);subplot(211);plot(ff,abs(fft(u,q*FS)/FS);subplot(212);pwelch(u,Rs);b,a=butter(13,1/20);%重構(gòu)濾波器H,F=FREQZ(b,a,FS,Rs);figure(6);plot(F,20*log10(abs(H);uoft=fil

40、ter(b,a,u);%基帶信號(D)figure(7);subplot(211);plot(t(80:480),real(uoft(80:480);subplot(212);plot(t(80:480),imag(uoft(80:480);figure(8);subplot(211);plot(ff,abs(fft(uoft,q*FS)/FS);subplot(212);pwelch(uoft,Rs);%Upconverters_tilde=(uoft.).*exp(1i*2*pi*fc*t);s=real(s_tilde);%通頻帶信號(E)figure(9);plot(t(80:480)

41、,s(80:480);figure(10);subplot(211);%plot(ff,abs(fft(real(uoft).).*cos(2*pi*fc*t),q*FS)/FS);%plot(ff,abs(fft(imag(uoft).).*sin(2*pi*fc*t),q*FS)/FS);plot(ff,abs(fft(s,q*FS)/FS);subplot(212);%pwelch(real(uoft).).*cos(2*pi*fc*t),Rs);%pwelch(imag(uoft).).*sin(2*pi*fc*t),Rs);pwelch(s,Rs);附錄2：OFDM接收部分代碼%DV

42、B2K模式接收部分clearall;closeall;Tu=224e-6;%有用OFDM符號持續(xù)時間T=Tu/2048;%原始基帶周期G=0;%允許保護時間間隔選擇1/4,1/8,1/16或1/32delta=G*Tu;%保護間隔持續(xù)時間Ts=delta+Tu;%整個OFDM符號持續(xù)時間Kmax=1705;%子載波數(shù)Kmin=0;FS=4096;%IFFT/FFT長度q=10;%載波周期與原始基帶周期比fc=q*l/T;%載波頻率Rs=4*fc;%模擬周期t=0:l/Rs:Tu;tt=0:T/2:Tu;%數(shù)據(jù)產(chǎn)生程序sM=2;x,y=meshgrid(-sM+l):2:(sM-l),(-sM+

43、l):2:(sM-l);alphabet=x(:)+li*y(:);N=Kmax+l;rand(state,0);a=-l+2*round(rand(N,l).+i*(-l+2*round(rand(N,l).;A=length(a);info=zeros(FS,l);info(l:(A/2)=a(l:(A/2).;info(FS-(A/2)-l):FS)=a(A/2)+l):A).;carriers=FS.*ifft(info,FS);%UpconverterL=length(carriers);chips=carriers.;zeros(2*q)-1,L);p=1/Rs:1/Rs:T/2;g=ones(length(p),1);dummy=conv(g,chips(:);u=dummy;zeros(46,1);b,aa=butter(13,1/20);uoft=filter(b,aa,u