%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%性能評估


save saved_data.mat;    % 保存數據


snr = [Eb_NoStart:Eb_NoInterval:Eb_NoEnd];

% 以下均考慮第u個用戶的誤比特率性能
u = 1;

% 誤比特率性能曲線

err_rate_user = sum(user_bit_err{u},1)./( length(user_bit{u})*N_frame) ;

figure(1);
semilogy(snr,err_rate_user,'b-*'); 
hold on ;
grid on ;
xlabel('Eb/No (dB)','FontSize',12);
ylabel('BER','FontSize',12);
title('系統BER~Eb/No性能曲線');


% 信道估計器的MSE性能曲線

mse_ce = sum(mse_ce_acc{u},1)/N_frame;

figure(2);
semilogy(snr,mse_ce,'k-v');
hold on ;
grid on ;
xlabel('Eb/No (dB)','FontSize',12);
ylabel('MSE','FontSize',12);
title('信道估計器的MSE性能曲線');


% 幀定時的誤差性能曲線
err_timing1 = sum( abs(frame_timing_acc{u}),1 )/N_frame;

figure(3);
stem(snr,err_timing1,'k-s');
hold on ;
grid on ;
axis([0 20 0 100]);
xlabel('Eb/No (dB)','FontSize',12);
ylabel('定時誤差(樣點數)','FontSize',12);
title('幀定時的誤差性能曲線');

% 載波同步誤差MSE性能曲線

MSE_freq = sum( abs(freq_syn_acc{u}).^2 ,1 )/N_frame;

figure(4);
semilogy(snr,MSE_freq,'k-s');
hold on ;
grid on ;
xlabel('Eb/No (dB)','FontSize',12);
ylabel('載波同步誤差MSE','FontSize',12);
title('載波同步誤差MSE性能曲線');