clc;
clear?all;
close?all;
M?=16;?%M進(jìn)制
SNRin=-10:1:20;
x?=?randi（21001）-1;%信源
k=length（x）/4;
x1=zeros（k1）;
for?i=1:k
????x1（i）=8*x（4*（i-1）+1）+4*x（4*（i-1）+2）+2*x（4*（i-1）+3）+x（4*i）;
end
y=modulate（modem.qammod（M）x1）;???%調(diào)制??modem.qammod（M）確定調(diào)制方式
ynoisy?