物理层安全截获概率散点图

clc;
clear all;
x1=[-4.87097 -4.13978 -3.32258 -2.11828 -1 -0.0107527 1.02151 1.96774 3 3.90323 4.89247 5.88172 6.95699 8.03226 9.02151 9.83871 10.914 11.7742 12.6774 13.1935];
y1=[0.387882 0.338799 0.288027 0.213879 0.150452 0.105835 0.0705274 0.0482878 0.0288775 0.01873 0.0109021 0.0063457 0.00322621 0.00168523 0.000880285 0.000512383 0.000253546 0.000139806 7.50312e-005 5.13714e-005];
%scatter(x1,y1,'k');
semilogy(x1,y1,'-<b');
axis([-5,15,0.0001,1]);
hold on
x2=[-4.95699 -4.13978 -3.15054 -1.94624 -0.913978 0.333333 1.4086 2.56989 3.60215 4.50538 5.53763 6.52688 7.55914 8.54839 9.53763 10.3978 11.5161 12.4194 13.4086 14.4839 15];
y2=[0.444075 0.409451 0.367448 0.320951 0.280338 0.231965 0.191938 0.146436 0.11172 0.0875729 0.0668121 0.0469988 0.0330611 0.0214434 0.0139082 0.00926829 0.00539474 0.00322621 0.00177894 0.00098091 0.000768894];
%scatter(x2,y2,'o');
semilogy(x2,y2,'-or');
axis([-5,15,0.0001,1]);
hold on
x3=[-4.87097 -4.05376 -3.02151 -1.94624 -0.956989 -0.0107527 0.978495 2.01075 3 3.98925 4.97849 6.05376 6.95699 8.03226 8.93548 9.92473 11 11.9892 13.0215 13.9677 14.914];
y3=[0.827446 0.805357 0.742564 0.722741 0.631285 0.566526 0.46877 0.398521 0.312383 0.238327 0.163175 0.114785 0.0744494 0.0496123 0.0313195 0.01823 0.0109021 0.0063457 0.00340562 0.00209252 0.00115382];
%scatter(x3,y3,'s');
semilogy(x3,y3,'-sg');
axis([-5,15,0.0001,1]);
hold on
x4=[-5 -4.05172 -3.06034 -2.11207 -1.03448 1.77636e-015 0.991379 2.02586 2.93103 3.92241 5.0431 5.94828 6.98276 8.01724 9.00862 10.0431 11.1207 12.0259 13.1034 14.2241 14.9569];
y4=[0.709494 0.687699 0.646098 0.607013 0.552774 0.503381 0.458402 0.392189 0.357146 0.296172 0.253393 0.216792 0.185478 0.149088 0.123635 0.0993779 0.0824118 0.070508 0.0566745 0.0441557 0.0366172];
%scatter(x4,y4,'d');
semilogy(x4,y4,'-dk');
xlabel('MER(dB)');
ylabel('Intercept Probability');
title('四种方案随MER变化的截获概率曲线');
legend('P-AFbORS','T-AFbORS','T-AFbMRC','常规直接传输（S-D）',3);