qam调制.m

QAM的调制与解调程序

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