Matlab中BPSK的BER计算问题：手动BER无法匹配SNR序列求解

我帮你找到了代码里的几个关键问题，修正后就能得到对应每个SNR值的手动BER序列了：

主要问题分析 #

1. BER赋值错误：你之前的err=error+1;既没有把结果绑定到数组的对应索引k上，也搞错了BER的计算逻辑——BER是错误比特数/总比特数，不是error+1。
2. 比特判决的对比逻辑错误：生成的s是BPSK的±1，但检测后的rd是0/1，直接对比两者的差值会导致错误的计数，需要把s转换回0/1的原始比特再做对比。
3. 噪声生成与判决的问题：BPSK是实值调制，你生成了复噪声，但判决时直接用复数和0比较会触发Matlab的语法错误，应该取接收信号的实部来判决，或者直接生成实噪声。

修正后的完整代码 #

N = 100; % 总比特数
SNR = 0:1:10; % 不同SNR值
err = zeros(1, length(SNR)); % 初始化手动BER数组

for k = 1:length(SNR)
    N0 = 10^(-SNR(k)/10); % 计算噪声功率
    % 生成0/1原始比特，转换为BPSK的±1
    s_bits = randn(1, N) > 0.5; % 0/1比特序列
    s = 2*s_bits - 1; % BPSK调制：1→1，0→-1
    
    % 生成实噪声（BPSK是实调制，复噪声无必要；若用复基带需取实部判决）
    n = sqrt(N0/2) * randn(1, N); 
    r = s + n; % 接收信号
    
    % 判决接收信号：实部大于0判为1，否则0
    rd = real(r) > 0; % 直接生成逻辑数组，简化代码
    
    % 计算错误比特数
    error = sum(rd ~= s_bits); % 直接对比0/1比特，统计不同的数量
    
    % 计算BER并赋值到对应SNR的位置
    err(k) = error / N; 
end

% 计算理论BER
th_err = 0.5 * erfc(sqrt(10.^(SNR/10))); 

% 绘图对比
semilogy(SNR, th_err, '*r', SNR, err, 'bo-'); 
grid on;
title('BPSK比特误码率对比');
xlabel('Eb/N0 (dB)');
ylabel('BER');
legend('理论BER', '手动计算BER');

关键修正点说明 #

• BER计算：用err(k) = error / N替代错误的err=error+1，确保每个SNR对应的BER都存储在数组的对应位置。
• 比特对比：保留原始的0/1比特s_bits，和判决后的rd直接对比，用sum(rd ~= s_bits)快速统计错误数，去掉冗余的嵌套循环，代码更高效。
• 噪声与判决：改为实噪声，判决时直接用real(r) > 0（如果坚持用复噪声，必须取实部，否则Matlab会报错）。
• 代码简化：用逻辑数组直接生成rd，减少不必要的变量初始化和循环，让代码更简洁易读。

现在运行修正后的代码，你就能得到每个SNR值对应的手动BER序列，并且可以和理论BER一起绘图对比了。

内容的提问来源于stack exchange，提问作者Mohamed Osman