qammod2.m

OFDM系统仿真

%qammod2.m
%16QAM调制程序
function[iout,qout]=qammod2(paradata,para,nd,ml)
%******************variables*************************
%paradata:输入数据
%iout:I信道输出数据
%qout:Q信道输出数据
%para:并行信道数
%nd:数据数
%ml:调制数
%(QPSK->2 16QAM->4)
%*****************************************************
%The constellation power
k=sqrt(10);
iv=[-3-1 3 1];
m2=ml/2;
count2=0;
for ii=1:nd
isi=zeros(para,1);
isq=zeros(para,1);
for jj=1:ml
if jj<=m2
  isi=isi+2.^(m2-jj).*paradata((1:para),count2+jj);
else
  isq=isq+2.^(ml-jj).*paradata((1:para),count2+jj);
 end
 end
iout((1:para),ii)=reshape(iv(isi+1)./k,para,1);
qout((1:para),ii)=reshape(iv(isq+1)./k,para,1);
count2=count2+ml;
end
%********************end of file***************************