《OFDM-matlab实现.ppt》由会员分享,可在线阅读,更多相关《OFDM-matlab实现.ppt(9页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、串/并变换QPSK映射FFT移去保护时隙QPSK逆映射IFFT上变频插入保护时隙并/串转换串/并变换并/串转换下变频信道串行二进制数据输入串行二进制数据输出SNR=10; %信噪比取值信噪比取值FL=128; %fft信号长度信号长度NS=6; %一个帧结构中一个帧结构中ofdm信号的个数信号的个数PARA=128; %子载波个数子载波个数GL=32; %保护时隙的长度保护时隙的长度%*signal generation*%Signal=rand(1,PARA*NS*2)0.5;%产生产生0,1随即序列随即序列for i=1:PARA%串并变换串并变换 for j=1:NS*2 SigPara
2、(i,j)=Signal(i*j); endendfigure(1)subplot(211);stem(Signal(1:50),.);051015202530354000.20.40.60.81前50个序列% *QPSK映射映射 *%for j=1:NS ichannal(:,j)=SigPara(:,2*j-1); qchannal(:,j)=SigPara(:,2*j);endk=1./sqrt(2);ichannal1=ichannal.*k;qchannal1=qchannal.*k;x=ichannal1+qchannal1.*sqrt(-1);y=ifft(x);ichannal
3、2=real(y);qchannal2=imag(y);05010015020025030035040045050000.10.20.30.40.50.60.7用用ifft/fft实现的多载波调制系统是在频域上进行的实现的多载波调制系统是在频域上进行的00.20.40.60.800.20.40.60.8数据x的星座图020406080100120140-0.100.10.20.30.4020406080100120140160-0.100.10.20.30.4%* 插入保护时隙插入保护时隙 *%ichannal3=ichannal2(FL-GL+1:FL,:);ichannal2;qchann
4、al3=qchannal2(FL-GL+1:FL,:);qchannal2;%* 并串变换并串变换*%ichannal4=reshape(ichannal3,1,(FL+GL)*NS);qchannal4=reshape(qchannal3,1,(FL+GL)*NS);Tdata=ichannal4+qchannal4.*sqrt(-1);插入循环前缀减小码间串扰,如图所示,后20个序列搬到序列的前面%* recieved data *%Rdata=awgn(Tdata,SNR,measured);%* remove pretection interval *%idata=real(Rdata
5、);qdata=imag(Rdata);idata1=reshape(idata,FL+GL,NS);qdata1=reshape(qdata,FL+GL,NS);idata2=idata1(GL+1:FL+GL,:);qdata2=qdata1(GL+1:FL+GL,:);%* FFT *%Rex=idata2+qdata2*sqrt(-1);ry=fft(Rex);Reichannal=real(ry)./k;Reqchannal=imag(ry)./k;%*QPSK reverse reflection *%for j=1:NS Repara(:,2*j-1)=Reichannal(:,j); Repara(:,2*j)=Reqchannal(:,j);endRes=reshape(Repara,1,);Resig=Res0.5;figure(3)stem(Resig(1:40),.);解调部分是调制逆过程,不再演示解调部分是调制逆过程,不再演示右边给出解调部分的代码:右边给出解调部分的代码:上变频部分只要将载波乘以cos(wct)将整个频谱搬移到wc处下变频是与之相反的过程。该过程与Am调制类似051015202530354000.20.40.60.81最终得到的接收端的前50个序列与发送端相同