througput.m

主要是针对现在的合作通信中

eb_n0_sd=0:0.25:6.5;
y=[0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.0000,0.000,0.000,0.000,0.000,0.000,0.000000,0.006464,0.011614,0.025755,0.076923,0.119275,0.157116,0.235984,0.317500,0.387117,0.459371];
eb_n0_sd0=0:0.25:6.00;
z=[0.000,0.000,0.000,0.000,0.0000,0.000,0.158730,0.625000,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286,0.714286];
plot(eb_n0_sd0,z,'g');
hold on;
eb_n0_sd1=0:0.25:6.5;
w=[0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.087108,0.084746,0.191638,0.357143,0.435464,0.569581,0.644716,0.695489,0.707786,0.713004];
plot(eb_n0_sd1,w,'b');
hold on;
plot(eb_n0_sd,y,'r');
hold on;
xlabel('Eb/N0(dB,s-d)');
ylabel('throughput');
title('the curve of throughput relating to Eb/N0');
legend('co-EbN0-sr=5dB','co-EbN0-sr=1dB','pure truncated ARQ');
hold;