(word完整版)OFDM-MATLAB仿真程序.doc

(word完整版)OFDMMATLAB仿真程序
(word完整版)OFDMMATLAB仿真程序
(word完整版)OFDMMATLAB仿真程序
HYPERLINK”http：//web。mit。edu/brzezin/www/359/matlab/OFDM.m"OFDM。m：OFDMSimulator(outerfunction)
clearall;

A=［11/exp(1）1/exp(2)］;％powerdelayprofile
N=64；%numberofsymbolsinasingleOFDMsymbol
GI=16；％guardinterval
Mt=1;%numberofTxantennas
Mr=1；％numberofRxantennas
sig2=1e—3;%noisevariance
M=8；%maxconstellationbitnumber
Mgap=10.^(1:（1.7/10）：2。7)；%gap
Btot=100*Mt；%total＃bitsperOFDMsymbol
TransmitIter=50;％#iterationsofsymboltransmissionsforeachchannelinstance
ChannelIter=100；％#iterationsofindependentidenticallydistributedchannelinstances
GapIter=length(Mgap）；

loadENC2.mat
loadENC4.mat
loadENC16.mat
loadENC64.mat
loadENC256.mat

TotEbNo=［]；
Errors=［］；
EbNo=[];
forlGap=1:GapIter
lGap
gap=Mgap（lGap)；
totalErrors=0;
forlChan=1：ChannelIter

％createchannel
[Hh_f]=create_channel（Mt,Mr,A，N+GI）;

%decomposeeachsubchannelinthefrequencydomain
[USV］=svd_decompose_channel（Mt，Mr,h_f，N）;

%bitloading
［bits_alloc，energy_alloc］=BitLoad（S,Btot，Mt＊N,gap，sig2,M）;
%energy_alloc=energy_alloc/(mean(energy_alloc））；
％energy_alloc=ones（1,128);

forlTrans=1：TransmitIter

％bitstotransmit
x=（randn（1，Btot）>0);

%modulate
x_mod=modulate(x,bits_alloc，energy_alloc，s2,s4，s16，s64，s256）；

％precodemodulatedsignal
x_pre=precode（Mt，x_mod,V，N);

%ifft,withcyclicprefixforeachantenna
ofdm_symbol=[];
fori=1:Mt
ofdm_symbol=［ofdm_symbol；ifft_cp_tx_blk（x_pre(i：Mt：Mt＊(N-1)+i）,N，GI)］;
end
ofdm_symbol2=reshape(ofdm_symbol,Mt*(N+GI),1)；

％channel
y=transpose（channel（sig2,Mt，Mr，ofdm_symbol2，H,N+GI)）；

％fft
rec_symbol=[］;
fori=1:Mt
rec_symbol=［rec_symbol；fft_cp_rx_blk（y（i:Mt:Mt*（N+GI—1)+i)，N,GI）］；
end

rec_symbol2=reshape(rec_symbol，1，Mt*N）；

％shapereceivedsignal
shaped_vals=shape（rec_symbol2,Mr,U,N);

%demodulate
y_demod=demodulate(shaped_vals,bits_alloc，energy_alloc，S,s2,s4,s16,s64,s256，c2,c4，c16,c64，c256）;

%comparison
totalErrors=totalErrors+sum(xor(y_demod，x)）；
end
EbNo=［EbNosum(energy_alloc）/Btot/sig2]；
end

Errors=[ErrorstotalErrors/Btot/Channe