16PSK眼图，星座图，误码率，噪声特性的研究与仿真

clear all;
clc;
close all;

M=16;
k=log2（M）;
n=100000;                          %比特序列长度
samp=1;                            %过采样率
 x=randi（[01][1n]）;                    %生成随机二进制比特流
 stem（x（1:50）‘filled‘）;            %画出相应的二进制比特流信号
 title（‘二进制随机比特流‘）;
 xlabel（‘比特序列‘）;ylabel（‘信号幅度‘）;

x4=reshape（xklength（x）/k）;       %将原始的二进制比特序列每四个一组分组，并排列成k行length（x）/k列的矩阵
 xsym=bi2de（x4.‘‘left-msb‘）;       %将矩阵转化为相应的16进制信号序列
figure;
 stem（xsym（1:50））;                  %画出相应的16进制信号序列
xlabel（‘信号序列‘）;ylabel（‘信号幅度‘）;
y=modulate（modem.pskmod（M）xsym）;  %用16PSK调制器对信号进行调制
scatterplot（y）;                    %画出16PSK信号的星座图
text