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（‘二進(jìn)制隨機(jī)比特流‘）;
xlabel（‘比特序列‘）;
ylabel（‘信號幅度‘）;

x6=reshape（xklength（x）/k）;
xsym=bi2de（x6.‘‘left-msb‘）;%2到32進(jìn)制序列

figure;
stem（xsym（1:50））;???%畫出相應(yīng)的16進(jìn)制信號序列????32
title（‘16進(jìn)制隨機(jī)信號‘）;
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）;