误码率计算用MATLAB仿真

误码率是衡量数字通信系统性能的一个重要指标，可以通过MATLAB进行仿真计算。下面以二进制基带传输为例，介绍如何用MATLAB计算BPSK调制下的误码率。

假设发送端发送的二进制序列为`s`，接收端接收到的二进制序列为`s_hat`，则误码率可以计算为：

BER = sum(s~=s_hat)/length(s);

其中，`~= `表示不等于运算符，`sum`函数用于计算误码数，除以序列长度即可得到误码率。

下面是一个简单的BPSK调制的MATLAB代码：

clc; clear; close all;

% 参数设置
N = 10^5;         % 发送的二进制序列长度
Eb_No_dB = 0:1:10;  % 信噪比范围
SNRdB = Eb_No_dB + 10*log10(1);  % 信噪比
Eb = 1;           % 符号能量
Es = Eb;          % 符号能量等于比特能量
BER = zeros(1, length(SNRdB));  % 初始化误码率

% 生成发送端二进制序列
s = randi([0 1], 1, N);

% BPSK调制
s_mod = 2*s - 1;

% 发送信号加噪声并接收
for i = 1:length(SNRdB)
    noise = sqrt(Eb/2/10^(SNRdB(i)/10))*randn(1, N);  % 生成高斯噪声
    r = s_mod + noise;                                 % 发送信号加噪声
    r_demod = sign(r);                                 % BPSK解调
    s_hat = (r_demod + 1)/2;                           % 二进制解调
    BER(i) = sum(s ~= s_hat)/N;                        % 计算误码率
end

% 绘制误码率曲线
semilogy(Eb_No_dB, BER, 'b-o');
xlabel('信噪比（dB）');
ylabel('误码率');
title('BPSK调制下的误码率曲线');
grid on;

代码中，首先设置参数，包括发送的二进制序列长度、信噪比范围、符号能量、误码率等。然后生成一个随机的二进制序列，进行BPSK调制，并添加高斯噪声。接收端进行BPSK解调和二进制解调，计算误码率。最后用`semilogy`函数绘制了BPSK调制下的误码率曲线。

需要注意的是，误码率计算的准确性与序列长度有关，通常需要进行多次实验取平均值。