threeber.m

来自「主要是针对现在的合作通信中」· M 代码 · 共 17 行

eb_n0_sd=0:0.25:6.5;
y=[1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.0000,1.000,1.000,1.000,1.000,1.000,1.000000,0.947368,0.906977,0.800000,0.513043,0.329502,0.235546,0.084912,0.030111,0.008989,0.002051];
eb_n0_sd0=0:0.25:6.00;
z=[1.000,1.000,1.000,1.000,1.0000,1.000,0.777778,0.125000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0000];
plot(eb_n0_sd0,z,'g');
hold on;
eb_n0_sd1=0:0.25:6.5;
w=[1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,0.878049,0.881356,0.731707,0.500000,0.390351,0.202586,0.097397,0.026316,0.009099,0.001795];
plot(eb_n0_sd1,w,'b');
hold on;
plot(eb_n0_sd,y,'r');
hold on;
xlabel('Eb/N0(dB,s-d)');
ylabel('PER');
title('the curve of PER relating to Eb/N0');
legend('co-EbN0-sr=5dB','co-EbN0-sr=1dB','pure truncated ARQ');
hold;