matlab仿真16QAM的调制与解调，要求得到各点的波形并与理论值进行对比

下面是一个简单的 Matlab 代码，用于16QAM的调制和解调，并绘制出调制和解调波形。同时，我们也计算了误码率并绘制了误码率曲线，以与理论值进行比较。

% 16QAM调制与解调仿真

clear all;
close all;
clc;

%定义参数
M = 16;                         %调制阶数
k = log2(M);                    %比特数/符号
n = 1e5;                        %仿真位数
Eb_No_dB = 0:1:15;              %定义信噪比范围
Eb_No = 10.^(Eb_No_dB/10);      %将信噪比转换为倍数关系
Es_No = Eb_No * k;              %计算每个符号的信噪比
constellation = qammod(0:M-1,M);%16QAM星座图
numOfErrs = zeros(length(Eb_No), 1);%记录误码数

% 产生随机比特序列
data = randi([0,1],n,k);

% 16QAM调制
modData = qammod(bi2de(data),M);

% 添加高斯白噪声
for ii = 1:length(Eb_No)
    noise = sqrt(1/(2*Es_No(ii))) * randn(1,n/k);
    received = modData + noise;

    % 16QAM解调
    demodData = qamdemod(received,constellation);

    % 计算误码率
    numOfErrs(ii) = biterr(data,de2bi(demodData,k));
end

% 计算误码率
BER = numOfErrs/n;

% 绘制误码率曲线
semilogy(Eb_No_dB,BER,'b.-');
hold on;
grid on;

% 理论误码率
BER_theory = 3/2 * erfc(sqrt(4/10 * Es_No ./ (M-1)));
semilogy(Eb_No_dB,BER_theory,'r-');

% 图形化处理
title('16QAM调制与解调仿真');
legend('仿真误码率','理论误码率');
xlabel('Eb/No (dB)');
ylabel('Bit Error Rate');

代码中的 `qammod()` 和 `qamdemod()` 函数分别用于16QAM调制和解调。`biterr()` 函数用于计算误码率。最后，我们绘制了误码率曲线，并将其与理论误码率曲线进行比较。