16QAM的星座图

%%%%%%%%% 绘制16QAM的星座图 %%%%%%%%
M=16;
SNRin=-10:1:20; %定义对数信噪比,求出来的是信号的功率谱密度

%%%%% randi函数：生成0 1的单极性序列
%%%% 第一个参数：生成随机序列的最大值 第2个参数：行数 第3个参数：列数
x=randi(M,5000,1)-1; %信源 %产生5000行1列的随机序列，最大值为M-1

y=modulate(modem.qammod(M),x); %调制，modem.qammod(M)确定调制方式
figure(1);%plot(y);
scatterplot(y);
title(‘原信号星座图’);axis([-4 4 -4 4]);

snr=10;ynoisy=awgn(y,snr); %加入高斯白噪声，信噪比SNR=0 y = awgn(x,snr)
figure(2);scatterplot(ynoisy);
title(‘接收信号星座图(SNR=10)’);axis([-4 4 -4 4]);

snr=0;ynoisy=awgn(y,snr); %加入高斯白噪声，信噪比SNR=0 y = awgn(x,snr)
figure(3);scatterplot(ynoisy);
title(‘接收信号星座图(SNR=0)’);axis([-4 4 -4 4]);