Как построить кривую BER против SNR для двоичных данных

Question

У меня есть два вектора измерения (130,1). Один для переданных данных, а другой для полученных данных.Используя следующий код, я получил BER и количество ошибок в векторе.Проблема в том, что я понятия не имею, как приступить к построению кривой BER против SNR.Хотя SNR может изменяться от 1 до 30 с шагом 1, BER, в то время как плоттер, кажется, не меняется, и я получаю прямую линию вместо кривой водопада.Пожалуйста, помогите, спасибо заранее.

clc
clear all
close all

A = 'sample-25-4-18.xlsx';
sheet1 = 1;
sheet2 = 2
Range = 'A1:AN140';
B = xlsread (A ,sheet1, Range);
C = xlsread (A ,sheet2, Range);

N = 10^4 % number of bits or symbols
rand('state',100); % initializing the rand() function
randn('state',200); % initializing the randn() function

% Transmitter
Eb_N0_dB = [-3:10]; % multiple Eb/N0 values

for ii = 1:length(Eb_N0_dB)

   % counting the errors
   nErr(ii) = size(find([B - C]),1);

end
simBer = nErr/N; % simulated ber
theoryBer = 0.5*erfc(sqrt(10.^(Eb_N0_dB/10))); % theoretical ber

% plot
close all
figure
semilogy(Eb_N0_dB,theoryBer,'b.-');
hold on
semilogy(Eb_N0_dB,simBer,'b*');
axis([-3 10 10^-5 0.5])
grid on
legend('theory', 'simulation');
xlabel('Eb/No, dB');
ylabel('Bit Error Rate');
title('Bit error probability curve for BPSK modulation');

Сюжет выглядит следующим образом.

Пожалуйста, помогите.Заранее спасибо.