ofdm.m

是一个根据最新LTE标准写的一个链路实现程序和信道估计程序

function [TrData_p1,TrData_p2]=ofdm(data1,daopin1,data2,daopin2)
N_count=13;
para=180;
fl = 1024;
gl=32;
SNR=10;
%规定相关参数
send_data1=data1+daopin1;
send_data2=data2+daopin2;
%产生发射天线的数据

% 进行IFFT运算
y1_ifft = ifft(send_data1,fl); 
y2_ifft = ifft(send_data2,fl);
%通过傅立叶反变换,将频域数据转换为时域数据
ich2_p1 = real(y1_ifft);      %I信道取变换后的实部
qch2_p1 = imag(y1_ifft);      %Q信道取变换后的虚部
ich2_p2 = real(y2_ifft);
qch2_p2 = imag(y2_ifft);
%插入保护间隔
ich3_p1 = [ich2_p1(fl-gl+1:fl,:);ich2_p1];
qch3_p1 = [qch2_p1(fl-gl+1:fl,:);qch2_p1];
ich3_p2 = [ich2_p2(fl-gl+1:fl,:);ich2_p2];
qch3_p2 = [qch2_p2(fl-gl+1:fl,:);qch2_p2];
%并串转换
ich4_p1 = reshape(ich3_p1,1,(fl+gl)*N_count);
qch4_p1 = reshape(qch3_p1,1,(fl+gl)*N_count);
ich4_p2 = reshape(ich3_p2,1,(fl+gl)*N_count);
qch4_p2 = reshape(qch3_p2,1,(fl+gl)*N_count);
%形成复数发射数据
TrData_p1 = ich4_p1+qch4_p1.*sqrt(-1);
TrData_p2 = ich4_p2+qch4_p2.*sqrt(-1);