matlab

Development Platform:

CHM

performance_eval.m：Code Content

							%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%性能评估
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性能曲线');

						