clc;
clear;
close?all;
%%初始化
n=1920;?%二進制比特數量
x=randi（[0:1]n1）;
snr=20;?%%dB

%%QPSK
x0=reshape（x[2n/2]）;????????%%每個symbol代表3個bit?3x640
x0=x0‘;??%%640x3
xsym_QPSK=bi2de（x0）;?%%640x1?將每3個比特轉換為對應的十進制，為一個symbol
y0=pskmod（xsym_QPSK4）;?%%qpsk調制
y0_noisy=awgn（y0snr‘measured‘）;

scatterplot（y010‘*r‘）;title（‘QPSK?Ideal‘）;
scatterplot（y0_noisy10‘o‘）;title（‘QPSK?noisy‘）;

%%8QAM
x1=reshape（x[3n/3]）;????????%%每個symbol代表3個bit?3x640
x1=x1‘;??%%640x3
xsym_qam8=bi2de（x1）;?%%640x1?將每3個比特轉換為對應的十進制，為一個symbol
y1=qammod（xsym_qam88）;?%%8QAM調制
y11=（[real（y1）imag（y1）]）;
y1_noisy=awgn（y1snr‘measured‘）;
y11_noisy=（[real（y1_noisy）imag（y1_noisy）]）;
scatterplot（y110‘*r‘）;title（‘8QAM?Ideal‘）;
hold?on
scatterplot（y1_noisy10‘o‘）;title（‘8QAM?noisy‘）;

%%16QAM
x2=reshape（x[4n/4]）;????????%%每個symbol代表3個bit?3x640
x2=x2‘;??%%640x3