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

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

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

% 清空环境变量和命令窗口
clc;
clear;

% 设置参数
M = 2; % 调制阶数
k = log2(M); % 每个符号的比特数
EbNoVec = (0:10); % 信噪比范围
numBits = 1e6; % 发送的总比特数
numSym = numBits/k; % 发送的总符号数
m = 1; % 码元时间长度

% 初始化误码率向量
berVec = zeros(size(EbNoVec));

% 生成BPSK调制器和QPSK调制器
bpskMod = comm.BPSKModulator;
qpskMod = comm.QPSKModulator;

% 生成AWGN信道
awgnChan = comm.AWGNChannel('NoiseMethod', 'Signal to noise ratio (SNR)');

for i = 1:length(EbNoVec)
    % 计算信噪比
    EbNo = EbNoVec(i);
    snr = EbNo + 10*log10(k) - 10*log10(m);
    
    % 设置AWGN信道的信噪比
    awgnChan.SNR = snr;
    
    % 生成随机比特流
    dataIn = randi([0 1], numBits, 1);
    
    % BPSK调制
    bpskTx = bpskMod(dataIn);
    
    % QPSK调制
    qpskTx = qpskMod(dataIn);
    
    % 传输信号通过AWGN信道
    bpskRx = awgnChan(bpskTx);
    qpskRx = awgnChan(qpskTx);
    
    % BPSK解调
    bpskDemod = comm.BPSKDemodulator;
    dataOut = bpskDemod(bpskRx);

    % 计算误码率
    [~, berVec(i)] = biterr(dataIn, dataOut);
    
    % QPSK解调
    qpskDemod = comm.QPSKDemodulator;
    dataOut = qpskDemod(qpskRx);

    % 计算误码率
    [~, berVec(i)] = biterr(dataIn, dataOut);
end

% 绘制误码率-信噪比曲线
semilogy(EbNoVec, berVec, 'b.-');
grid on;
xlabel('Eb/No (dB)');
ylabel('BER');
title('BPSK and QPSK Modulation in AWGN');
legend('BPSK', 'QPSK');