%(7,4)HAMMING码系统仿真结果与未编码的比较main_ce55_74.m
%仿真序列长度2*10^4,运行时间约2分钟。仿真序列长度10^5,运行时间约20分钟
echo on
SNRindB55741=0:1:10;
SNRindB55742=0:1:10;
smld_err_prb55=zeros(1,length(SNRindB55741));
smld_err_prb74=zeros(1,length(SNRindB55741));
SNR5574=0;
theo_err_prb5574=zeros(1,length(SNRindB55742));%
for i=1:length(SNRindB55741)
smld_err_prb55(i)=smldPe55(SNRindB55741(i));
smld_err_prb74(i)=smldPe55_74(SNRindB55741(i));
end
for i=1:length(SNRindB55742)
SNR5574=exp(SNRindB55742(i)*log(10)/10);
theo_err_prb5574(i)=(1/2)*erfc(sqrt(2*SNR5574)/sqrt(2));
%Qfunct y=(1/2)*erfc(x/sqrt(2));
%theo_err_prb(i)=Qfunct(sqrt(2*SNR));
end
%[p_err_ha,gamma_db_ha]=p_e_hd_a(5,16,4,7,3);
%p_e_hd_a.m 硬判决性能估计函数
%[p_err_sa,gamma_db_sa]=p_e_sd_a(5,16,4,7,3);
%p_e_sd_a.m 软判决性能估计函数
%绘图语句
figure;
semilogy(SNRindB55741,smld_err_prb55,'b*-');
hold on
axis([0,10,1e-5,1]);
xlabel('Eb/N0(dB)');
ylabel('BER');
title('(7,4)HAMMING码,未编码与理论值的仿真结果比较');
semilogy(SNRindB55741,smld_err_prb74,'ro-');
semilogy(SNRindB55742,theo_err_prb5574,'bs:');
%semilogy(gamma_db_ha,p_err_ha,'g');
%semilogy(gamma_db_sa,p_err_sa,'m');
legend('未编码','(7,4)HAMMING码','理论值');
%ce55_74.m 仿真绘图语句