用MATLAB实现OFDM仿真分析

para=128;

fftlen=128;

noc=128;   

nd=6;      

ml=2;   

sr=2500;

br=sr.*ml;

gilen=18;  

ebn0=3;  

%************************** 主循环部分 **************************

nloop=100;

noe = 0;  

nod = 0;  

eop=0;    

nop=0;    

for iii=1:nloop

%****************** 发射部分 ******************************************************

%****************** 产生数据 **********************

seldata=rand(1,para*nd*ml)>0.5;        

%****************** 串并转换 ***********************

paradata=reshape(seldata,para,nd*ml); 

%****************** QPSK调制 ***********************

[ich,qch]=qpskmod(paradata,para,nd,ml);

kmod=1/sqrt(2);                      

ich1=ich.*kmod;

qch1=qch.*kmod;                      

%****************** IFFT(离散傅里叶反变换) *********

x=ich1+qch1.*i;

y=ifft(x);                       

ich2=real(y);                           

qch2=imag(y);                          

%****************** 插入保护间隔 *******************

[ich3,qch3]= giins(ich2,qch2,fftlen,gilen,nd);

fftlen2=fftlen+gilen;

%****************** 衰减计算 ************************

spow=sum(ich3.^2+qch3.^2)/nd./para;    

attn=0.5*spow*sr/br*10.^(-ebn0/10);

attn=sqrt(attn);

%****************** 高斯白信道 *****************************************************

[ich4,qch4]=comb(ich3,qch3,attn);

%****************** 接收部分  ******************************************************

%****************** 去除保护间隔 ********************

[ich5,qch5]= girem(ich4,qch4,fftlen2,gilen,nd);

%****************** FFT(离散傅里叶变换)  ************

rx=ich5+qch5.*i;

ry=fft(rx); 

ich6=real(ry);      

qch6=imag(ry);   

%***************** QPSK解调 **************************

ich7=ich6./kmod;

qch7=qch6./kmod;

[demodata]=qpskdemod(ich7,qch7,para,nd,ml);  

%**************  并串变换  ***************************

demodata1=reshape(demodata,1,para*nd*ml);

%*************** 比特误码率 (BER) *********************

% instantaneous number of error and data

noe2=sum(abs(demodata1-seldata));     

nod2=length(seldata);                  

% cumulative the number of error and data in noe and nod

noe=noe+noe2;

nod=nod+nod2;

% calculating PER

if noe2~=0 

   eop=eop+1;

else

   eop=eop;

end  

   eop;

   nop=nop+1;

fprintf('%dt%et%dn',iii,noe2/nod2,eop);  

end

%****************** 结果输出 ***************************

per=eop/nop;

ber=noe/nod;

fprintf('%ft%et%et%dtn',ebn0,ber,per,nloop);

fid = fopen('BERofdm.dat','a');

fclose(fid);

%****************** end of file ******************************************************

最后

以上就是清秀水杯为你收集整理的用MATLAB实现OFDM仿真分析的全部内容，希望文章能够帮你解决用MATLAB实现OFDM仿真分析所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。