qam_rayleigh.m

clc;
clear all

M=16;               %设置16QAM参量
k=log2（M）;
n=100000;           %比特序列长度
p=0.5;              %二进制比特流产生1的概率
samp=1;             %过采样率
x=randsrc（1n[10;p1-p]）;     %生成随机二进制比特流

figure（‘Name‘‘二进制随机比特流‘‘Numbertitle‘‘Off‘）
stem（x（11:50）‘filled‘）;         %画出相应的二进制比特流信号
title（‘二进制随机比特流‘）;
xlabel（‘比特序列‘）;ylabel（‘信号幅度‘）;

x4=reshape（xklength（x）/k）;      %将原始的二进制比特序列每四个一组分组，并排列成k行length（x）/k列的矩阵
xsym=bi2de（x4.‘‘left-msb‘）;      %将矩阵转化为相应的16进制信号序列
figure（‘Name‘‘16进制信号序列‘‘Numbertitle‘‘Off‘）
stem（xsym（1:50））;                 %画出相应的16进制信号序列
title（‘16进制随机信号‘）;
xlabel（‘信号序列‘）;ylabel（‘信号幅度‘）;

y=qammod（xsymM）;                     %用16QAM调制器对信号进行调制

scatterplot（y）;                       %画出16QAM信号的星座图