请写出完成BPSK和QPSK调制在AWGN信道下的SNR-BER误码率性能评估的matlab代码

以下是BPSK和QPSK调制在AWGN信道下的SNR-BER误码率性能评估的matlab代码：

% BPSK和QPSK调制在AWGN信道下的SNR-BER误码率性能评估

clear all; close all;

% 信号参数
M = 2; % 2PSK
% M = 4; % 4QPSK
k = log2(M);

% 构建调制映射表
if M == 2
    modTable = 1/sqrt(2)*[-1 1];
elseif M == 4
    modTable = 1/sqrt(2)*[-1-1i 1-1i -1+1i 1+1i];
else
    error('不支持的调制方式');
end

% 仿真参数
SNR_dB = 0:2:20; % 信噪比范围
numSym = 1e4; % 发送符号数

% 循环执行
for idxSNR = 1:length(SNR_dB)
    % 生成随机发送符号
    txSym = randi([0 M-1], 1, numSym);
    % 映射到调制符号
    modSym = modTable(txSym+1);
    % 加噪声
    SNR = 10^(SNR_dB(idxSNR)/10);
    noiseVar = 1/(2*k*SNR); % AWGN信道的噪声方差
    noise = sqrt(noiseVar)*(randn(size(modSym))+1i*randn(size(modSym)));
    rxSym = modSym + noise;
    % 解调
    if M == 2
        rxSymDec = real(rxSym)>0;
    elseif M == 4
        rxSymDec = zeros(size(rxSym));
        rxSymDec(real(rxSym)>0 & imag(rxSym)>0) = 3;
        rxSymDec(real(rxSym)>0 & imag(rxSym)<0) = 1;
        rxSymDec(real(rxSym)<0 & imag(rxSym)>0) = 2;
        rxSymDec(real(rxSym)<0 & imag(rxSym)<0) = 0;
    end
    % 统计误码率
    numErr = sum(txSym~=rxSymDec);
    ber(idxSNR) = numErr/numSym;
end

% 绘图
semilogy(SNR_dB, ber, '-o');
grid on;
xlabel('SNR (dB)');
ylabel('BER');
title('BPSK/QPSK in AWGN Channel');
if M == 2
    legend('BPSK');
elseif M == 4
    legend('QPSK');
end

注意事项：

- 代码中有两个注释需要根据调制方式选择注释掉一个，另一个保留。
- 代码中假设发送符号是等概率随机的，如果需要其他分布的符号，需要修改代码。
- 代码中假设AWGN信道的噪声是复高斯分布的，如果需要其他分布的噪声，需要修改代码。
- 代码中假设发送符号和接收符号的能量都是1，如果需要其他能量，需要修改代码。