OFDM調(diào)制解調(diào)仿真_第1頁
OFDM調(diào)制解調(diào)仿真_第2頁
OFDM調(diào)制解調(diào)仿真_第3頁
OFDM調(diào)制解調(diào)仿真_第4頁
OFDM調(diào)制解調(diào)仿真_第5頁
已閱讀5頁,還剩4頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、實驗四 OFDM調(diào)制解調(diào)仿真一、 實驗?zāi)康模?) 了解OFDM調(diào)制解調(diào)的原理。(2) 學(xué)會用星座圖分析系統(tǒng)性能。二、 實驗內(nèi)容(1) 編寫MATLAB程序,實現(xiàn)OFDM系統(tǒng)調(diào)制解調(diào)。(2) 繪制各步驟圖形并分析系統(tǒng)特性。三、 實驗代碼1、主代碼如下:clear all;close all;N=input('請輸入碼元數(shù)');SNR=input('請輸入信噪比');xx=randint(1,4*N); %原序列figure(1),stem(xx,'.k'); %原序列圖形title('原序列');B=0;for m=1:4:4*N

2、A=xx(m)*8+xx(m+1)*4+xx(m+2)*2+xx(m+3); B=B+1; ee(B)=A;endfigure(2),stem(ee,'.b');title('化為015的碼元');yy=star(ee,N);figure(3),plot(yy,'.r'); %映射后的星座圖title('映射后的星座圖'); ff=ifft(yy,N); %傅里葉反變換N1=floor(N*1/4);N3=floor(N*3/4);N5=floor(N*5/4);figure(4),stem(ff,'.m');ti

3、tle('傅里葉反變換后');for j=1:N1 %加循環(huán)前綴 ss(j)=ff(N3+j);endfor j=1:N %變成長度為N*5/4的序列 ss(N1+j)=ff(j);endfigure(5),stem(ss,'.k'); %畫出圖形title('加N/4循環(huán)前后綴');%ss=wgn(1,N5,0,10,'dBW','complex');ss=awgn(ss,SNR); %加入噪聲figure(6),stem(ss,'.m'); %加入噪聲后的圖形title('加入噪聲后&#

4、39;); zz=fft(ss(N1+1):N5),N); %傅里葉變換figure(7),plot(zz,'.b'); %畫圖title('傅里葉變換后');rr=istar(zz,N); %星座圖糾錯figure(8),plot(rr,'.r'); %畫圖title('糾錯后的星座圖');dd=decode(rr,N); %解碼figure(9),stem(dd,'.m'); %畫圖title('星座圖糾錯并解碼后');%bb=d2b(dd,N);bb=d2bb(dd,N); %轉(zhuǎn)化為0/1比特

5、流figure(10),stem(bb,'.b');title('轉(zhuǎn)化為0/1比特流后');2、調(diào)用的函數(shù):(1) 星座圖映射%星座圖映射function yy=star(xx,N)B=-3-3*i,-3-i,-1-3*i,-1-i,-3+3*i,-3+i,-1+3*i,-1+i,3-3*i,3-i,1-3*i,1-i,3+3*i,3+i,1+3*i,1+i;for j=1:N yy(j)=B(xx(j)+1);end(2) 星座圖逆映射%星座圖逆映射function rr=istar(zz,N)for j=1:N if(mod(floor(real(zz(j)

6、,2)=0) zz1(j)=ceil(real(zz(j); else zz1(j)=floor(real(zz(j); end if(mod(floor(imag(zz(j),2)=0) zz1(j)=zz1(j)+ceil(imag(zz(j)*i; else zz1(j)=zz1(j)+floor(imag(zz(j)*i; end rr(j)=zz1(j);end(3) 十進制轉(zhuǎn)二進制function bb=d2b(dd,N)%十進制轉(zhuǎn)二進制for j=1:N*4 bb(j)=;endfor j=1:4:N*4 bb1=dec2bin(dd(floor(j/4)+1),4); for

7、k=1:4 bb(4*(j-1)+k)=bb1(k); endend(4) 十進制轉(zhuǎn)化為01比特流%十進制轉(zhuǎn)化為01比特流function bb=d2bb(dd,N)for j=1:N*4 bb(j)=1;endj=1;while(j<=N*4) N1=ceil(j/4) a4=mod(dd(N1),2); dd(N1)=floor(dd(N1)/2); a3=mod(dd(N1),2); dd(N1)=floor(dd(N1)/2); a2=mod(dd(N1),2); dd(N1)=floor(dd(N1)/2); a1=mod(dd(N1),2); bb(j)=a1; j=j+1;

8、 bb(j)=a2;j=j+1;bb(j)=a3;j=j+1;bb(j)=a4;j=j+1;end(5) 星座圖逆映射function yy=decode(rr,N)for j=1:N switch(rr(j)%星座圖逆映射 case -3-3*i yy(j)=0; case -3-i yy(j)=1; case -1-3*i yy(j)=2; case -1-i yy(j)=3; case -3+3*i yy(j)=4; case -3+i yy(j)=5; case -1+3*i yy(j)=6; case -1+i yy(j)=7; case 3-3*i yy(j)=8; case 3-i yy(j)=9; case 1-3*i yy(j)=10; cas

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論