OFDM-MATLAB-仿真程序.doc

% - %B6)eu Zx0% 參數(shù)定義 %3TNZk9x)OI6xK0% - %kH#Bs4J o0NumLoop = 1000;中國通信網(wǎng)無限分享8C.e/AUyfNumSubc = 128;pK1dA?;O0NumCP = 8;k+y$E X&y6RW0SyncDelay = 0;中國通信網(wǎng)無限分享0d*m;*|.J/aC#i(X0% 子載波數(shù) 128!hb%j!MTuO;Y2V0h0% 位數(shù)/ 符號(hào) 2中國通信網(wǎng)無限分享$k6,N?% 符號(hào)數(shù)/ 載波 1000p;I3sOzRR0% 訓(xùn)練符號(hào)數(shù) 0中國通信網(wǎng)無限分享8S3V_!Deej% 循環(huán)前綴長度 8 (1/16)*T+H7g/6O2d)pJ0H 0% 調(diào)制方式 4-QAMReH3H Cq$y9C0% 多徑信道數(shù) 3 中國通信網(wǎng)無限分享(HF,N1u;xH% IFFT Size 128 中國通信網(wǎng)無限分享 !vu)fQ% 信道最大時(shí)延 2中國通信網(wǎng)無限分享%F-CW9 W|.I% - %:g?7&P$u3J-P(U1W0% QAM MODULATION %;r r asd6z|0% - %nS0KR2j0YSi|Md0% Generate the random binary stream for transmit testpk.j,FqU;z-y0BitsTx = floor(rand(1,NumLoop*NumSubc)*2);% Modulate (Generates QAM symbols)中國通信網(wǎng)無限分享tW6F0Wy/SK% input: BitsTx(1,NumLoop*NumSubc); output: SymQAM(NumLoop,NumSubc/2)中國通信網(wǎng)無限分享2w2x/eB4S2+SymQAMtmp = reshape(BitsTx,2,NumLoop*NumSubc/2).;fM8POZ) j0SymQAMtmptmp = bi2de(SymQAMtmp,2,left-msb);% QAM modulation 中國通信網(wǎng)無限分享 F+? f2H+q4D% 00-1-i,01-1+i,10-1-i,11-1+i5e-dxA!I:w0% 利用查表法進(jìn)行QAM星座映射GVa#W/9(UJy0QAMTable = -1-i -1+i 1-i 1+i;中國通信網(wǎng)無限分享:Y.j)ozh YSymQAM = QAMTable(SymQAMtmptmp+1);% - %L$T-X8_*nO!r C0% IFFT %中國通信網(wǎng)無限分享4M9$E#Cn% - %中國通信網(wǎng)無限分享hZT$vZ% input: SymQAM(NumLoop,NumSubc/2); output: SymIFFT(NumSubc,NumLoop)中國通信網(wǎng)無限分享+Ek/nj fR,B/ESymIFFT = zeros(NumSubc,NumLoop);:#h&P5r29W H&I8t% O0SymIFFTtmp = reshape(SymQAM,NumSubc/2,NumLoop);中國通信網(wǎng)無限分享wPO!v SymIFFTtmptmp = zeros(NumSubc,NumLoop);)A!_vb?F)Ci0SymIFFTtmptmp(1,:) = real(SymIFFTtmp(1,:); % 實(shí)數(shù)1Nmkgh0SymIFFTtmptmp(NumSubc/2+1,:) = imag(SymIFFTtmp(1,:); % 實(shí)數(shù)RX+FD0% 構(gòu)造共軛對稱矩陣中國通信網(wǎng)無限分享7q N m b2G!RCSymIFFTtmptmp(2:NumSubc/2,:) = SymIFFTtmp(2:NumSubc/2,:);中國通信網(wǎng)無限分享CPqE7bSymIFFTtmptmp(NumSubc/2+2):NumSubc,:) = flipdim(conj(SymIFFTtmp(2:NumSubc/2,:),1);SymIFFT = ifft(SymIFFTtmptmp,NumSubc,1);% - %)gB)N Rk*Lls0% Add cyclic prefix %中國通信網(wǎng)無限分享-k9r4Ia g p$HT&L% - %中國通信網(wǎng)無限分享:_1jkbI6eP9% input: SymIFFT(NumSubc,NumLoop); output: SymCP(NumSubc + NumCP,NumLoop)Y9JftM0NumAddPrefix = NumSubc + NumCP;中國通信網(wǎng)無限分享,nK3e,Z$g#b2lnSymCP = zeros(NumAddPrefix,NumLoop);gnB9bXu9p b*a0RowPrefix = (NumSubc - NumCP + 1):NumSubc;中國通信網(wǎng)無限分享V_W;9m1QNSymCP = SymIFFT(RowPrefix,:);SymIFFT;% - %ikG(Cp,_$v*db0% Go through the channel %0 ai0% - %0w!jv| gJ Q0% input: SymCP(NumSubc + NumCP,NumLoop); output: SymCh(1,(NumSubc + NumCP)*NumLoop)中國通信網(wǎng)無限分享D-Do5P0h1uSymCh = zeros(1,(NumSubc + NumCP)*NumLoop);中國通信網(wǎng)無限分享rSn0WSymChtmp = SymCP(:).;$m B&x$lw0Ch = 1 1/2 1/4;9a&u4XLI0I0SymChtmptmp = filter(Ch,1,SymChtmp);% Add the AWGNGd(S K0BerSnrTable = zeros(20,3);+a P#c$y0for snr=0:19; % = SNR + 10*log10(log2(2);6PJ81_Sej,x0yV)Y0 BerSnrTable(snr+1,1) = snr;-e2 u0n%qWM&|ka0SymCh = awgn(SymChtmptmp,snr,measured);% - %3tn&e0% Remove Guard Intervals %中國通信網(wǎng)無限分享 CGP%O3C/ewa% - %:G6k)K8A)O*B.BA0% input: SymCh(1,(NumSubc + NumCP)*NumLoop); output: SymDeCP(NumSubc,NumLoop)中國通信網(wǎng)無限分享r5hw*scx%c5?SymDeCP = zeros(NumSubc,NumLoop);&EbZ TG H&ZD8rlYq:r0SymDeCPtmp = reshape(SymCh,NumSubc + NumCP,NumLoop);9pv:n!Z&gw#J40SymDeCP = SymDeCPtmp(NumCP+1+SyncDelay):NumAddPrefix+SyncDelay,:);% - %V3t5E0N$0% FFT %中國通信網(wǎng)無限分享m_.Ev Z9O,w/% - %CQ-f%|K1z(k1JJ0% input: SymDeCP(NumSubc,NumLoop); output: SymFFT(NumSubc,NumLoop)MfIqOE!Y0SymFFT = fft(SymDeCP,NumSubc,1);% - %中國通信網(wǎng)無限分享 m.B i9k*O%l% Make Decision(Include DeQAM) %$JY2P(vjz0% - %o7fI?HM z4_0% SymFFT(NumSubc,NumLoop); output: SymDec(NumSubc,NumLoop)中國通信網(wǎng)無限分享 ACt lm5D8SymDec = zeros(NumSubc,NumLoop);+mO5Vgm!A#b0SymEqtmp(1,:) = SymFFT(1,:)+i*SymFFT(NumSubc/2+1,:);中國通信網(wǎng)無限分享#Uxj*dnJ-sdSymEqtmp(2:NumSubc/2,:) = SymFFT(2:NumSubc/2,:);中國通信網(wǎng)無限分享zA G$N:?2Lfor m = 1:NumLoop中國通信網(wǎng)無限分享a3U qUl* for n = 1:NumSubc/2中國通信網(wǎng)無限分享$q _IWL%vTr Real = real(SymEqtmp(n,m);)S:L |3J+Irk0 Imag = imag(SymEqtmp(n,m);中國通信網(wǎng)無限分享BQD+?Sd R x-K6Yar oMA0 if( abs(Real -1) abs(Real +1)中國通信網(wǎng)無限分享d6f!| SG Yd SymDec(2*n-1,m) = 1;中國通信網(wǎng)無限分享v*sA-W/p(AC elseR,aR3l0N(j5x+?0 SymDec(2*n-1,m) = 0;%|;jAQl7i0 end if( abs(Imag -1) abs(Imag +1 ) ) 中國通信網(wǎng)無限分享DfI2B C:KA#dR%qW SymDec(2*n,m) = 1;中國通信網(wǎng)無限分享 JAavVx else中國通信網(wǎng)無限分享;c5I8N? z SymDec(2*n,m) = 0; C T2H7+_0 end中國通信網(wǎng)無限分享$q&k5O Z(mG.x endp2Qzs0hf0a0end% Bit Error中國通信網(wǎng)無限分享de6 trF&CBitsRx = zeros(1,NumSubc*NumLoop);中國通信網(wǎng)無限分享2R4ko9IM%Mq5_BitsRx = SymDec(:).;中國通信網(wǎng)無限分享.3OG/cNum,Ber = symerr(BitsTx,BitsRx)ff eY)kO#S!q0BerSnrTable(snr+1,2) = Num ;中國通信網(wǎng)無限分享5EPyIQMq!?BerSnrTable(snr+1,3) = Ber ;7hnp5F k6|)a)P1dN8|0end,EI+X)G_.Ot0figure(1);中國通信網(wǎng)無限分享Y(S B2o3N Usubplot(2,1,1);中國通信網(wǎng)無限分享N3 R y0%W5f)LEsemilogy(BerSnrTable(:,1),BerSnrTable(:,2),o-);title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));title(誤比特?cái)?shù)曲線);xlabel(SNR/dB);ylabel(誤比特?cái)?shù));中國通信網(wǎng)無限分享6wrgc