16QAM研究，眼图，星座图，误码率，噪声的研究。

clear all;
close all;
clc;
M=16;
k=log2（M）;
n=100000;                          %比特序列长度
samp=1;                            %过采样率
x=randi（[12][1n]）-1;  
figure;%生成随机二进制比特流
stem（x（1:50）‘filled‘）;            %画出相应的二进制比特流信号
title（‘二进制随机比特流‘）;
xlabel（‘比特序列‘）;ylabel（‘信号幅度‘）;
x4=reshape（xklength（x）/k）;       %将原始的二进制比特序列每四个一组分组，并排列成k行length（x）/k列的矩阵
xsym=bi2de（x4.‘‘left-msb‘）;
figure%将矩阵转化为相应的16进制信号序列
stem（xsym（1:50））;                  %画出相应的16进制信号序列
title（‘16进制随机信号‘）;
xlabel（‘信号序列‘）;ylabel（‘信号幅度‘）;
 
y=modulate（modem.qammod（M）xsym）;  %用16QAM调制器对信号进行调制
scatterplot（y）;