OFDM_matlab实现.ppt

OFDM简易模型matlab simulation * * 串/并变换 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(i,j)=Signal(i*j); end end figure(1) subplot(211);stem(Signal(1:50),.); 前50个序列 % **********QPSK映射 *********% for j=1:NS ichannal(:,j)=SigPara(:,2*j-1); qchannal(:,j)=SigPara(:,2*j); end k=1./sqrt(2); ichannal1=ichannal.*k; qchannal1=qchannal.*k; x=ichannal1+qchannal1.*sqrt(-1); y=ifft(x); ichannal2=real(y); qchannal2=imag(y); 用ifft/fft实现的多载波调制系统是在频域上进行的 数据x的星座图 %****** 插入保护时隙 *******% ichannal3=[ichannal2(FL-GL+1:FL,:);ichannal2]; qchannal3=[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); 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); end Res=reshape(Repara,1,[]); Resig=Res0.5; figure(3) stem(Resig(1:40),.); 解调部分是调制逆过程，不再演示 右边给出解调部分的代码： 上变频部分只要将载波乘以cos(wct)将整个频谱搬移到wc处 下变频是与之相反的过程。该过程与Am调制类似 *