G1 = [
          100            0            0      0.25379      0.23057            0
          100            3            0       0.1718      0.15755            0
          100            6            0       0.1108     0.087356            0
          100            9            0     0.052356     0.024713            0
          100           12            0     0.015718    0.0023851            0
          100           15            0    0.0088649   0.00020115            0
];

G2 = [
          100            0            0      0.22519      0.20494            0
          100            3            0      0.16433      0.14147            0
          100            6            0     0.091299     0.071169            0
          100            9            0     0.031169     0.017857            0
          100           12            0     0.010292    0.0019156            0
          100           15            0    0.0037662   0.00012987            0
];

hold off
subplot(2,2,1)
semilogy(G1(:,2),G1(:,5),'k-')
hold on
semilogy(G1(:,2),G1(:,4),'k-.')
xlabel 'SNR (dB)'
ylabel 'BER'
hold off

% a = axis;
% 
% axis( [0 18 a(3) a(4)] );

x = 0.5; % 11; % a(1) + (a(2)-a(1))/20;
y = 0.3e-2; % 1.2e-1; % a(3) + (a(4)-a(3))/20;

text(x,y, ['.-.-.- EM               '; ...
           '\_\_\_\_\_ known channel']);


title 'pilot period = 10'

hold off
subplot(2,2,2)
semilogy(G2(:,2),G2(:,5),'k-')
hold on
semilogy(G2(:,2),G2(:,4),'k-.')
xlabel 'SNR (dB)'
ylabel 'BER'
hold off

% a = axis;
% 
% axis( [0 18 a(3) a(4)] );

x = 0.5; % 11; % a(1) + (a(2)-a(1))/20;
y = 0.3e-2; % 1.2e-1; % a(3) + (a(4)-a(3))/20;

text(x,y, ['.-.-.- EM               '; ...
           '\_\_\_\_\_ known channel']);


title 'pilot period = 5'
