QAM误码率与SNR的MATLAB仿真

仿真在LTE中16QAM和64QAM信号的SNR与误码率间的关系。即实现了误码率与SNR间的关系，又实现了误比特率与SNR间的关系，还实现了添加高斯白噪声时的实际效果与理论值间的差别比较

close all
clear all

%M=64;
M=16;
k=log2（M）;
L=sqrt（M）;
n=120000;                     %比特序列长度
samp=1;                       %过采样率

x=randi（[01]1n）;           %生成随机二进制比特流
stem（x（1:50）‘filled‘）;       %画出相应的二进制比特流信号
title（‘二进制随机比特流‘）;
xlabel（‘比特序列‘）;
ylabel（‘信号幅度‘）;

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

%y=modulate（qammod（xsymM）xsym）;    %用64QAM调制器对信号进行
y=qammod（xsymM）.‘;             %用64QAM调制器对信号进行
scatterplot（y）;                 %画出64QAM信号的星座图
%text（real（y）imag（y）dec2bin（xsym））;
axis（[-L L -L L]）;

a = 3*log2（L）;
b = L^（2）-1;
c = 1-L^（-1）;

EbNo=0:0.1:14;
for n=1:length（EbNo）
  snr（n）=EbNo（n）+10*log10（k）-10*log10（samp）; %信噪比
  yn=awgn（ysnr（n）‘measured‘）;              %加入高斯白噪声
  yd=qamdemod（ynM）;                         %此时解