16qam仿真及注解

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 信道信噪比