%% 測試 ofdm，畫測試曲線
%%
%% 本測試為  有無均衡對系統性能的影響
%%
close all;clear all; clc; tic ; %計時開始
global CP_len;
global nSubC;
global ifft_len;
global bitPerSymbol;
global symbolPerCarrier;
global s_len;
global fade;        %信道衰弱

%% 測試參數
%SNR = 40; % dB
CP_len = 0;
nSubC = 48;
ifft_len = 256;
bitPerSymbol = 2;
symbolPerCarrier = 50;
fade = [0, 0, 0, 0.4, 0.3]; % 信道參數，下標為延時點數，值為衰減系數

% 開始測試
s_len = nSubC*bitPerSymbol* symbolPerCarrier;

SNR = 0:2:20;
ber1 = zeros(1,length(SNR));
ber2 = zeros(1,length(SNR));
ber3 = zeros(1,length(SNR));

for kk = 1:length(SNR) 
    max = floor((2^kk)/10)+1;   
    SNR()
    fade = [0]; % 信道參數，下標為延時點數，值為衰減系數
    ber = 0;
    for m = 1:max
        % 生成用于測試的隨機比特
        s_in = floor( rand(1, s_len )*2 );    
        sch = ofdm_mod_eq(s_in);
        sch = channel(sch,SNR(kk));
        s_out = ofdm_demod_eq(sch);    
        %計算誤比特率
        err_count = sum(abs(s_in-s_out));
        ber = ber + err_count/s_len;
    end
    ber1(kk) = ber/max;    
    
    fade = [0, 0, 0, 0, 0.4]; % 信道參數，下標為延時點數，值為衰減系數
    ber = 0;
    for m = 1:max
        % 生成用于測試的隨機比特
        s_in = floor( rand(1, s_len )*2 );    
        sch = ofdm_mod_eq(s_in);
        sch = channel(sch,SNR(kk));
        s_out = ofdm_demod_eq(sch);    
        %計算誤比特率
        err_count = sum(abs(s_in-s_out));
        ber = ber + err_count/s_len;
    end
    ber2(kk) = ber/max ;
    
    fade = [0, 0, 0, 0, 0.4, 0.3]; % 信道參數，下標為延時點數，值為衰減系數
    ber = 0;
    for m = 1:max
        % 生成用于測試的隨機比特
        s_in = floor( rand(1, s_len )*2 );    
        sch = ofdm_mod_eq(s_in);
        sch = channel(sch,SNR(kk));
        s_out = ofdm_demod_eq(sch);    
        %計算誤比特率
        err_count = sum(abs(s_in-s_out));
        ber = ber + err_count/s_len;
    end
    ber3(kk) = ber/max ;    
end
figure(2);
semilogy(SNR,ber1,'o-');
hold on;
semilogy(SNR,ber2,'r*-');
hold on;
semilogy(SNR,ber3,'gv-');

xlabel('SNR in dB');
ylabel('Bit Error Rate');
legend('AWGN','AWGN + 2 path','AWGN + 3 path');
grid;
toc