function?[?iout_qout?]?=?qammod（paradataparandm1?）
%********************?variables?*************************************

%?paradata??:?input?data?（para-by-nd?matrix）
%?iout??????:?output?Ich?data
%?qout??????:?output?Qch?data
%?para??????:?Number?of?parallel?channels
%?nd????????:?Number?of?data
%?m1????????:?Number?of?modulation?levels
%?（?QPSK-2??16QAM-4）

%********************************************************************

%?The?constellation?power
para=1;
nd=1024;
m1=4;
paradata=rand（1nd*m1）>0.5;