close?all
clear?all

%M=64;
M=16;
k=log2（M）;
L=sqrt（M）;
n=120000;?????????????????????%比特序列長度
samp=1;???????????????????????%過采樣率

x=randi（[01]1n）;???????????%生成隨機二進制比特流
stem（x（1:50）‘filled‘）;???????%畫出相應的二進制比特流信號
title（‘二進制隨機比特流‘）;
xlabel（‘比特序列‘）;
ylabel（‘信號幅度‘）;

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

%y=modulate（qammod（xsymM）xsym）;????%用64QAM調制器對信號進行
y=qammod（xsymM）.‘;?????????????%用64QAM調制器對信號進行
scatterplot（y）;?????????????????%畫出64QAM信號的星座圖
%text（real（y）imag（y）dec2bin（xsym））;
axis（[-L?L?-L?L]）;

a?=?3*log2（L）;
b?=?L^（2）-1;
c?=?1-L^（-1）;

EbNo=0:0.1:14;
for?n=1:length（EbNo）
??snr（n）=EbNo（n）+10*log10（k）-10*log10（samp）;?%信噪比
??yn=awgn（ysnr（n）‘measured‘）;??????????????%加入高斯白噪聲
??yd=qamdemod（ynM）;?????????????????????????%此時解