OFDM_matlab实现.ppt

OFDM简易模型matlabsimulation串/并变换SNR=10; %信噪比取值
FL=128; %fft信号长度
NS=6; %一个帧结构中ofdm信号的个数
PARA=128;%子载波个数
GL=32; %保护时隙的长度

%****signalgeneration****************%

Signal=rand(1,PARA*NS*2)>0.5;
%产生0，1随即序列
fori=1:PARA %串并变换
forj=1:NS*2
SigPara(i,j)=Signal(i*j);
end
end
figure(1)
subplot(211);stem(Signal(1:50),'.');
%**********QPSK映射*********%
forj=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);%******插入保护时隙*******%
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);%**********recieveddata***********%
Rdata=awgn(Tdata,SNR,'measured');

%**********removepretectioninterval*********%
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;

%***********QPSKreversereflection************%
forj=1:NS
Repara(:,2*j-1)=Reichannal(:,j);
Repara(:,2*j)=Reqchannal(:,j);
end
Res=reshape(Repara,1,[]);
Resig=Res>0.5;
figure(3)
stem(Resig(1:40),'.');最终得到的接收端的前50个序列与发送端相同