通信原理QAM调制和星座图

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