帮我写MATLAB代码：仿真双极性不归零码通过高斯白噪声信道后的误码率性能，画出误码率曲线。

好的，以下是MATLAB代码：

clear all;
close all;

%% 参数设置
N = 100000; % 发送比特数
Eb = 1; % 每个比特的能量
EbN0dB_min = 0; % 最小信噪比(dB)
EbN0dB_max = 10; % 最大信噪比(dB)
EbN0dB_step = 1; % 信噪比步长(dB)
sigma = 1; % 高斯白噪声标准差

%% 构造双极性不归零码
bit_stream = randn(1,N)>0; % 生成随机比特流
bp_data = 2*bit_stream - 1; % 双极性不归零码

%% 误码率计算
for i=1:length(EbN0dB_min:EbN0dB_step:EbN0dB_max)
    %% 发送信号
    EbN0dB = EbN0dB_min + (i-1)*EbN0dB_step; % 当前信噪比(dB)
    % 计算噪声能量
    N0 = Eb / (10^(EbN0dB/10));
    % 生成高斯白噪声
    noise = sigma*randn(1,N);
    % 接收信号
    r = bp_data + noise;
    
    %% 解码
    bit_hat = r>0; % 硬判决
    % 统计错误比特数
    num_err(i) = sum(bit_stream~=bit_hat);
    
    %% 绘制误码率曲线
    BER(i) = num_err(i)/N; % 当前信噪比下的误码率
end

%% 绘图
semilogy(EbN0dB_min:EbN0dB_step:EbN0dB_max,BER,'bo-');
grid on;
axis([EbN0dB_min EbN0dB_max 1e-6 1]);
xlabel('Eb/N0(dB)');
ylabel('Bit Error Rate');
title('Bipolar NRZ-I Code over AWGN Channel');

希望能对你有所帮助！