在这段代码后面加一个计算最后输出的信号和最开始基带信号的曲麻吕

好的，以下是修改后的代码，添加了计算输出信号和基带信号的曲麻吕的部分：

M = 4; % 调制阶数
k = log2(M); % 每个符号的比特数
numBits = 1e6; % 发送的比特数
SNR_dB = 0:2:16; % 信噪比范围
numSymbols = numBits/k; % 发送的符号数

% 构造发送符号
bits = randi([0 1], numBits, 1);
x = qammod(bits, M);

% 构造加性高斯白噪声
for i = 1:length(SNR_dB)
    SNR = 10^(SNR_dB(i)/10);
    noise_var = 1/SNR;
    noise = sqrt(noise_var/2)*(randn(numSymbols, 1) + 1i*randn(numSymbols, 1));

    % 接收信号
    y = x + noise;

    % 检测符号，并计算误码率
    bits_hat = qamdemod(y, M);
    numErr(i) = sum(bits ~= bits_hat);
    ber(i) = numErr(i)/numBits;

    % 计算输出信号和基带信号的曲麻吕
    t = (0:numSymbols-1).';
    s = exp(-1i*2*pi*t/numSymbols);
    c = x.*s;
    d = y.*s;
    OF = abs(sum(c))^2/numSymbols;
    AF = abs(sum(d))^2/numSymbols;
    SER(i) = AF/OF;
end

% 画图
semilogy(SNR_dB, ber, 'o-', 'LineWidth', 2);
hold on;
semilogy(SNR_dB, SER, 'x-', 'LineWidth', 2);
xlabel('信噪比 (dB)');
ylabel('误码率/曲麻吕');
legend('误码率', '曲麻吕');
grid on;

在这个代码中，我们添加了计算输出信号和基带信号的曲麻吕的部分。我们首先计算了发送符号的离散时间傅里叶变换（DTFT）系数 `s`。然后，我们分别计算了发送信号和接收信号的 DTFT 系数（分别为 `c` 和 `d`）。最后，我们计算了输出信号和基带信号的曲麻吕，即 `OF` 和 `AF`，并计算了曲麻吕误码率（SER）。最后，我们使用 `hold on` 命令，同时绘制误码率和曲麻吕之间的关系图。