performance_eval.m
上传用户:look542
上传日期:2009-06-04
资源大小:784k
文件大小:2k
源码类别:

传真(Fax)编程

开发平台:

Matlab

  1. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  2. %性能评估
  3. save saved_data.mat;    % 保存数据
  4. snr = [Eb_NoStart:Eb_NoInterval:Eb_NoEnd];
  5. % 以下均考虑第u个用户的误比特率性能
  6. u = 1;
  7. % 误比特率,误符号率性能曲线
  8. ber = sum(user_bit_err{u},1)./( length(user_bit{u})*N_frame) ;
  9. figure(1);
  10. semilogy(snr,ber,'b-*'); 
  11. hold on ;
  12. grid on ;
  13. xlabel('Eb/No (dB)','FontSize',12);
  14. ylabel('BER','FontSize',12);
  15. title('系统BER~Eb/No性能曲线');
  16. % figure(2);
  17. % semilogy(snr,ser,'b-*'); 
  18. % hold on ;
  19. % grid on ;
  20. % xlabel('Eb/No (dB)','FontSize',12);
  21. % ylabel('BER','FontSize',12);
  22. % title('系统SER~Eb/No性能曲线');
  23. % 信道估计器的MSE性能曲线
  24. mse_ce = sum(mse_ce_acc{u},1)/N_frame;
  25. figure(2);
  26. semilogy(snr,mse_ce,'k-v');
  27. hold on ;
  28. grid on ;
  29. xlabel('Eb/No (dB)','FontSize',12);
  30. ylabel('MSE','FontSize',12);
  31. title('信道估计器的MSE性能曲线');
  32. % nmse_ce = sum(nmse_ce_acc{u},1)/N_frame;
  33. % figure(2);
  34. % semilogy(snr,nmse_ce,'k-v');
  35. % hold on ;
  36. % grid on ;
  37. % xlabel('Eb/No (dB)','FontSize',12);
  38. % ylabel('NMSE','FontSize',12);
  39. % title('信道估计器的归一化MSE性能曲线');
  40. % 帧定时的误差性能曲线
  41. err_timing1 = sum( abs(frame_timing_acc{u}),1 )/N_frame;
  42. % figure(3);
  43. % stem(snr,err_timing1,'k-s');
  44. % hold on ;
  45. % grid on ;
  46. % axis([0 20 0 100]);
  47. % xlabel('Eb/No (dB)','FontSize',12);
  48. % ylabel('定时误差(样点数)','FontSize',12);
  49. % title('帧定时的误差性能曲线');
  50. % 载波同步误差MSE性能曲线
  51. MSE_freq = sum( abs(freq_syn_acc{u}).^2 ,1 )/N_frame;
  52. NMSE_freq = MSE_freq/delta_fc^2;
  53. % figure(4);
  54. % semilogy(snr,MSE_freq,'k-s');
  55. % hold on ;
  56. % grid on ;
  57. % xlabel('Eb/No (dB)','FontSize',12);
  58. % ylabel('载波同步误差MSE','FontSize',12);
  59. % title('载波同步误差MSE性能曲线');
  60. % figure(4);
  61. % semilogy(snr,NMSE_freq,'k-s');
  62. % hold on ;
  63. % grid on ;
  64. % xlabel('Eb/No (dB)','FontSize',12);
  65. % ylabel('载波同步误差NMSE','FontSize',12);
  66. % title('载波同步误差归一化MSE性能曲线');