plot.m

来自「AWGN信道下高斯噪声干扰下GALILEO编码译码系统性能仿真」· M 代码 · 共 33 行

clc 
clear all
tic
SNR_in_DB=-4:0.5:10;
snr_in_DB=-4:0.5:10;

for i=1:length(SNR_in_DB),
   snr=10^(SNR_in_DB(i)/10);%计算理想误码率
   x=sqrt(2*snr);
   y(i)=(1/2)*erfc(x/sqrt(2));
end; 

for i=1:length(snr_in_DB)
   [py(i),pb(i)]=main(snr_in_DB(i));%py为没有经过信道编码的误码率，pb为经过信道编码的误码率
end;

echo on;
% semilogy(SNR_in_DB,y);%误码率理想特性曲线
% hold on
% semilogy(snr_in_DB,py,'b');
% semilogy(snr_in_DB,pb,'r');
semilogy(SNR_in_DB,y,'-o',snr_in_DB,py,'-+',snr_in_DB,pb,'-*');
xlabel('SNR');
ylabel('BER(db)');
% set(gca,'Ytick',[1*10^(-5)]);
legend('theory','noencoding','encoding');
grid on;
% py
% pb
toc
%  plot(SNR_in_DB,log10(y));%误码率理想特性曲线
% plot(SNR_in_DB,y);