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）;????????????????????