资源描述：
《LTE下行链路OFDM系统的MATLAB实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.clearall;%%%%%%%%%参数设置部分%%%%%%%%%%Nsp=52;%系统子载波数（不包括直流载波）Nfft=64;%FFT长度Ncp=16;%循环前缀长度Ns=Nfft+Ncp;%一个完整OFDM符号长度noc=53;%包含直流载波的总的子载波数Nd=6;%每帧包含的完整的OFDM符号数（不包括训练符号)M1=4;%QPSKM2=16;%16-QAM调制sr=250000;%OFDM符号速率EbN0=0:2:30;%归一化信噪比Nfrm=100;%每种信噪比下的仿真帧数ts=1/sr/Ns;%OFDM符号抽样时间间隔t=0
2、:ts:(Ns*(Nd+1)*Nfrm-1)*ts;%抽样时刻fd=100;%最大多普勒频移h=rayleigh(fd,t);%生成单径Rayleigh衰落信道*%训练符号频域数据，采用802.11a中的长训练符号数据Preamble=[11-1-111-11-1111111-1-111-11-11111...1-1-111-11-11-1-1-1-1-111-1-11-11-11111];Preamble1=zeros(1,Nfft);Preamble1(2:27)=Preamble(27:end);%训练符号重排后的数据Preamble
3、1(39:end)=Preamble(1:26);preamble1=ifft(Preamble1);%训练符号时域数据..preamble1=[Preamble1(Nfft-Ncp+1:end)preamble1];%加入循环前缀%%%%%%%%%%仿真循环%%%%%%%%%%forii=1:length(EbN0)%********************发射机部分*******************msg1=randsrc(Nsp,Nd*Nfrm,[0:M1-1]);%QPSK信息数据msg2=randsrc(Nsp,Nd*Nfrm
4、,[0:M2-1]);%16-QAM信息数据data1=pskmod(msg1,M1,pi/4);%QPSK调制data2=qammod(msg2,M2)/sqrt(10);%16-QAM调制并归一化data3=zeros(Nfft,Nd*Nfrm);%根据FFT要求，对数据重排data4=zeros(Nfft,Nd*Nfrm);data3(2:27,:)=data1(27:end,:);data3(39:end,:)=data1(1:26,:);data4(2:27,:)=data2(27:end,:);data4(39:end,:)=d
5、ata2(1:26,:);cleardata1data2;%清除不需要的临时变量data3=ifft(data3);%IFFT变换data4=ifft(data4);data3=[data3(Nfft-Ncp+1:end,:);data3];%加入循环前缀data4=[data4(Nfft-Ncp+1:end,:);data4];..spow1=norm(data3,'fro').^2/(Nsp*Nd*Nfrm);%计算数据符号能量spow2=norm(data4,'fro').^2/(Nsp*Nd*Nfrm);data5=zeros(Ns
6、,(Nd+1)*Nfrm);%加入训练符号data6=data5;forindx=1:Nfrmdata5(:,(indx-1)*(Nd+1)+1)=preamble1.';data5(:,(indx-1)*(Nd+1)+2:indx*(Nd+1))=data3(:,(indx-1)*Nd+1:indx*Nd);data6(:,(indx-1)*(Nd+1)+1)=preamble1.';data6(:,(indx-1)*(Nd+1)+2:indx*(Nd+1))=data4(:,(indx-1)*Nd+1:indx*Nd);endclear
7、data3data4;data5=reshape(data5,1,Ns*(Nd+1)*Nfrm);%并串转换data6=reshape(data6,1,Ns*(Nd+1)*Nfrm);sigma1=sqrt(1/2*spow1/log2(M1)*10.^(-EbN0(ii)/10));%根据EbN0计算噪声标准差sigma2=sqrt(1/2*spow2/log2(M2)*10.^(-EbN0(ii)/10));forindx=1:Nfrmdd1=data5((indx-1)*Ns*(Nd+1)+1:indx*Ns*(Nd+1));%当前帧
8、的发射数据dd2=data6((indx-1)*Ns*(Nd+1)+1:indx*Ns*(Nd+1));hh=h((indx-1)*Ns*(Nd+1)+1:indx*Ns*(Nd
