performance_eval.m

OFDM的详细仿真代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%性能评估


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