clc;
clear?all

M=16;???????????????%設置16QAM參量
k=log2（M）;
n=100000;???????????%比特序列長度
p=0.5;??????????????%二進制比特流產生1的概率
samp=1;?????????????%過采樣率
x=randsrc（1n[10;p1-p]）;?????%生成隨機二進制比特流

figure（‘Name‘‘二進制隨機比特流‘‘Numbertitle‘‘Off‘）
stem（x（11:50）‘filled‘）;?????????%畫出相應的二進制比特流信號
title（‘二進制隨機比特流‘）;
xlabel（‘比特序列‘）;ylabel（‘信號幅度‘）;

x4=reshape（xklength（x）/k）;??????%將原始的二進制比特序列每四個一組分組，并排列成k行length（x）/k列的矩陣
xsym=bi2de（x4.‘‘left-msb‘）;??????%將矩陣轉化為相應的16進制信號序列
figure（‘Name‘‘16進制信號序列‘‘Numbertitle‘‘Off‘）
stem（xsym（1:50））;?????????????????%畫出相應的16進制信號序列
title（‘16進制隨機信號‘）;
xlabel（‘信號序列‘）;ylabel（‘信號幅度‘）;

y=qammod（xsymM）;?????????????????????%用16QAM調制器對信號進行調制

scatterplot（y）;???????????????????????%畫出16QAM信號的星座圖