performance_eval.m
资源名称:MIMO-OFDM.rar [点击查看]
上传用户:look542
上传日期:2009-06-04
资源大小:784k
文件大小:2k
源码类别:
传真(Fax)编程
开发平台:
Matlab
- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
- %性能评估
- save saved_data.mat; % 保存数据
- snr = [Eb_NoStart:Eb_NoInterval:Eb_NoEnd];
- % 以下均考虑第u个用户的误比特率性能
- u = 1;
- % 误比特率,误符号率性能曲线
- ber = sum(user_bit_err{u},1)./( length(user_bit{u})*N_frame) ;
- figure(1);
- semilogy(snr,ber,'b-*');
- hold on ;
- grid on ;
- xlabel('Eb/No (dB)','FontSize',12);
- ylabel('BER','FontSize',12);
- title('系统BER~Eb/No性能曲线');
- % figure(2);
- % semilogy(snr,ser,'b-*');
- % hold on ;
- % grid on ;
- % xlabel('Eb/No (dB)','FontSize',12);
- % ylabel('BER','FontSize',12);
- % title('系统SER~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性能曲线');
- % nmse_ce = sum(nmse_ce_acc{u},1)/N_frame;
- % figure(2);
- % semilogy(snr,nmse_ce,'k-v');
- % hold on ;
- % grid on ;
- % xlabel('Eb/No (dB)','FontSize',12);
- % ylabel('NMSE','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;
- NMSE_freq = MSE_freq/delta_fc^2;
- % figure(4);
- % semilogy(snr,MSE_freq,'k-s');
- % hold on ;
- % grid on ;
- % xlabel('Eb/No (dB)','FontSize',12);
- % ylabel('载波同步误差MSE','FontSize',12);
- % title('载波同步误差MSE性能曲线');
- % figure(4);
- % semilogy(snr,NMSE_freq,'k-s');
- % hold on ;
- % grid on ;
- % xlabel('Eb/No (dB)','FontSize',12);
- % ylabel('载波同步误差NMSE','FontSize',12);
- % title('载波同步误差归一化MSE性能曲线');