QAM16.m

clc;
clear
close all;
M=64;
k=log2（M）;
n=120000;
samp=1;%过采样率？？？？？？？？？
%snr=0:1:14;


x=randint（n1）;
stem（x（1:50）‘filled‘）;
title（‘二进制随机比特流‘）;
xlabel（‘比特序列‘）;
ylabel（‘信号幅度‘）;

x6=reshape（xklength（x）/k）;
xsym=bi2de（x6.‘‘left-msb‘）;%2到32进制序列

figure;
stem（xsym（1:50））;   %画出相应的16进制信号序列    32
title（‘16进制随机信号‘）;
xlabel（‘信号序列‘）;
ylabel（‘信号幅度‘）;

y=modulate（modem.qammod（M）xsym）;
scatterplot（y）;

text（real（y）+0.1imag（y）dec2bin（xsym））;%?

axis（[-8 8 -8 8]）;

snr=15+10*log10（k）-10*log10（samp）;
%snrtem=10.^（snr/10）;%信噪比？？？？
%pn=1./snrtem;
%sigma=sqrt（pn）;