OFDM-MATLAB-仿真程序.doc

% - %B6)eu Zx0% 参数定义 %3TNZk9x)OI6xK0% - %kH#Bs4J o0NumLoop = 1000;中国通信网无限分享8C.e/AUyfNumSubc = 128;pK1dA?;O0NumCP = 8;k+y$E X&y6RW0SyncDelay = 0;中国通信网无限分享0d*m;*|.J/aC#i(X0% 子载波数 128!hb%j!MTuO;Y2V0h0% 位数/ 符号 2中国通信网无限分享$k6,N?% 符号数/ 载波 1000p;I3sOzRR0% 训练符号数 0中国通信网无限分享8S3V_!Deej% 循环前缀长度 8 (1/16)*T+H7g/6O2d)pJ0H 0% 调制方式 4-QAMReH3H Cq$y9C0% 多径信道数 3 中国通信网无限分享(HF,N1u;xH% IFFT Size 128 中国通信网无限分享 !vu)fQ% 信道最大时延 2中国通信网无限分享%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)中国通信网无限分享tW6F0Wy/SK% input: BitsTx(1,NumLoop*NumSubc); output: SymQAM(NumLoop,NumSubc/2)中国通信网无限分享2w2x/eB4S2+SymQAMtmp = reshape(BitsTx,2,NumLoop*NumSubc/2).;fM8POZ) j0SymQAMtmptmp = bi2de(SymQAMtmp,2,left-msb);% QAM modulation 中国通信网无限分享 F+? f2H+q4D% 00-1-i,01-1+i,10-1-i,11-1+i5e-dxA!I:w0% 利用查表法进行QAM星座映射GVa#W/9(UJy0QAMTable = -1-i -1+i 1-i 1+i;中国通信网无限分享:Y.j)ozh YSymQAM = QAMTable(SymQAMtmptmp+1);% - %L$T-X8_*nO!r C0% IFFT %中国通信网无限分享4M9$E#Cn% - %中国通信网无限分享hZT$vZ% input: SymQAM(NumLoop,NumSubc/2); output: SymIFFT(NumSubc,NumLoop)中国通信网无限分享+Ek/nj fR,B/ESymIFFT = zeros(NumSubc,NumLoop);:#h&P5r29W H&I8t% O0SymIFFTtmp = reshape(SymQAM,NumSubc/2,NumLoop);中国通信网无限分享wPO!v SymIFFTtmptmp = zeros(NumSubc,NumLoop);)A!_vb?F)Ci0SymIFFTtmptmp(1,:) = real(SymIFFTtmp(1,:); % 实数1Nmkgh0SymIFFTtmptmp(NumSubc/2+1,:) = imag(SymIFFTtmp(1,:); % 实数RX+FD0% 构造共轭对称矩阵中国通信网无限分享7q N m b2G!RCSymIFFTtmptmp(2:NumSubc/2,:) = SymIFFTtmp(2:NumSubc/2,:);中国通信网无限分享CPqE7bSymIFFTtmptmp(NumSubc/2+2):NumSubc,:) = flipdim(conj(SymIFFTtmp(2:NumSubc/2,:),1);SymIFFT = ifft(SymIFFTtmptmp,NumSubc,1);% - %)gB)N Rk*Lls0% Add cyclic prefix %中国通信网无限分享-k9r4Ia g p$HT&L% - %中国通信网无限分享:_1jkbI6eP9% input: SymIFFT(NumSubc,NumLoop); output: SymCP(NumSubc + NumCP,NumLoop)Y9JftM0NumAddPrefix = NumSubc + NumCP;中国通信网无限分享,nK3e,Z$g#b2lnSymCP = zeros(NumAddPrefix,NumLoop);gnB9bXu9p b*a0RowPrefix = (NumSubc - NumCP + 1):NumSubc;中国通信网无限分享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)中国通信网无限分享D-Do5P0h1uSymCh = zeros(1,(NumSubc + NumCP)*NumLoop);中国通信网无限分享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 %中国通信网无限分享 CGP%O3C/ewa% - %:G6k)K8A)O*B.BA0% input: SymCh(1,(NumSubc + NumCP)*NumLoop); output: SymDeCP(NumSubc,NumLoop)中国通信网无限分享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 %中国通信网无限分享m_.Ev Z9O,w/% - %CQ-f%|K1z(k1JJ0% input: SymDeCP(NumSubc,NumLoop); output: SymFFT(NumSubc,NumLoop)MfIqOE!Y0SymFFT = fft(SymDeCP,NumSubc,1);% - %中国通信网无限分享 m.B i9k*O%l% Make Decision(Include DeQAM) %$JY2P(vjz0% - %o7fI?HM z4_0% SymFFT(NumSubc,NumLoop); output: SymDec(NumSubc,NumLoop)中国通信网无限分享 ACt lm5D8SymDec = zeros(NumSubc,NumLoop);+mO5Vgm!A#b0SymEqtmp(1,:) = SymFFT(1,:)+i*SymFFT(NumSubc/2+1,:);中国通信网无限分享#Uxj*dnJ-sdSymEqtmp(2:NumSubc/2,:) = SymFFT(2:NumSubc/2,:);中国通信网无限分享zA G$N:?2Lfor m = 1:NumLoop中国通信网无限分享a3U qUl* for n = 1:NumSubc/2中国通信网无限分享$q _IWL%vTr Real = real(SymEqtmp(n,m);)S:L |3J+Irk0 Imag = imag(SymEqtmp(n,m);中国通信网无限分享BQD+?Sd R x-K6Yar oMA0 if( abs(Real -1) abs(Real +1)中国通信网无限分享d6f!| SG Yd SymDec(2*n-1,m) = 1;中国通信网无限分享v*sA-W/p(AC elseR,aR3l0N(j5x+?0 SymDec(2*n-1,m) = 0;%|;jAQl7i0 end if( abs(Imag -1) abs(Imag +1 ) ) 中国通信网无限分享DfI2B C:KA#dR%qW SymDec(2*n,m) = 1;中国通信网无限分享 JAavVx else中国通信网无限分享;c5I8N? z SymDec(2*n,m) = 0; C T2H7+_0 end中国通信网无限分享$q&k5O Z(mG.x endp2Qzs0hf0a0end% Bit Error中国通信网无限分享de6 trF&CBitsRx = zeros(1,NumSubc*NumLoop);中国通信网无限分享2R4ko9IM%Mq5_BitsRx = SymDec(:).;中国通信网无限分享.3OG/cNum,Ber = symerr(BitsTx,BitsRx)ff eY)kO#S!q0BerSnrTable(snr+1,2) = Num ;中国通信网无限分享5EPyIQMq!?BerSnrTable(snr+1,3) = Ber ;7hnp5F k6|)a)P1dN8|0end,EI+X)G_.Ot0figure(1);中国通信网无限分享Y(S B2o3N Usubplot(2,1,1);中国通信网无限分享N3 R y0%W5f)LEsemilogy(BerSnrTable(:,1),BerSnrTable(:,2),o-);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);title(误比特数曲线);xlabel(SNR/dB);ylabel(误比特数);中国通信网无限分享6wrgc