用MATLAB实现OFDM仿真分析

#附MATLAB仿真程序%**********************备部^分***************************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;foriii=1:nloop*******************发射部分*************************************************************************J产牛数**********************seldata二rand(l,para*nd*ml)>0.5;*******************串并转***********************paradata二reshape(seldata,para,nd*ml);%******************qpsk调制***********************[ich,qch]二qpskmod(paradata,para,nd,ml);kmod=l/sqrt(2);ich1=ich.*kmod;qch1=qch.*kmod;*******************IFFT(离散傅里叶反变换)*********x二ichl+qchl.*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);***************并串变***************************demodatai二reshape(demodata,l,para*nd*ml);****************比特误码率(BER)*********************%instantaneousnumberoferroranddatanoe2=sum(abs(demodatal-seldata));nod2=length(seldata);%cumulativethenumberoferroranddatainnoeandnodnoe=noe+noe2;nod二nod+nod2;%calculatingPERifnoe2~=0eop二eop+1;elseeop二eop;endeop;nop二nop+1;fprintf('%d\t%e\t%d\n',iii,noe2/nod2,eop);end*******************结果输出***************************per二eop/nop;ber二noe/nod;fprintf('%f\t%e\t%e\t%d\t\n',ebnO,ber,per,nloop);fid=fopen('BERofdm.