plot.m

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

clc
clear all
tic
SNR_in_DB=-4:0.3:10;
snr_in_DB=-4:0.3: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)
  [p1(i),p5(i)]=GPSmain(snr_in_DB(i));%py为没有经过信道编码的误码率，pb为经过信道编码的误码率
   p2(i)=GPSmain1(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,'-',snr_in_DB,p1,'-*',snr_in_DB,p2,'-+',snr_in_DB,p5,'k-o');
%semilogy(SNR_in_DB,y,'-',snr_in_DB,p1,'-*',snr_in_DB,p5,'-+');
xlabel('SNR');
ylabel('BER(db)');
legend('theory','GPSL1 encoding','GPSL2C encoding','GPSL5 encoding','Location','SouthWest');
grid on;
toc
% py
% pb

%  plot(SNR_in_DB,log10(y));%误码率理想特性曲线
% plot(SNR_in_DB,y);