M=16;
k=log2（M）;
n=100000;??????????????????????????%比特序列長度
samp=1;????????????????????????????%過采樣率
x=randint（n1）;????????????????????%生成隨機二進制比特流
stem（x（1:50）‘filled‘）;????????????%畫出相應的二進制比特流信號
title（‘二進制隨機比特流‘）;
xlabel（‘比特序列‘）;ylabel（‘信號幅度‘）;
x4=reshape（xklength（x）/k）;???????%將原始的二進制比特序列每四個一組分組，并排列成k行length（x）/k列的矩陣
xsym=bi2de（x4.‘‘left-msb‘）;???????%將矩陣轉化為相應的十進制信號序列?左邊為有效高位
figure;
stem（xsym（1:50））;??????????????????%畫出相應的十進制信號序列??畫出來是離散函數?桿圖
title（‘10進制隨機信號‘）;
xlabel（‘信號序列‘）;ylabel（‘信號幅度‘）;
y=modulate（modem.qammod（M）xsym）;??%用16QAM調制器對信號進行調制
scatterplot（y）;????????????????????%畫出16QAM信號的星座圖
text（real（y）+0.1imag（y）dec2bin（xsym））;%在指定位置上目標??這里是指在星座點上顯示二進制
axis（[-5?5?-5?5]）;%坐標軸顯示范圍
SNR_in_dB=8:2:24;?%AWGN?信道信噪比?