% plots for various alpha
%
hold off
subplot(2,2,1)
load gm.10.txt
G=gm;
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
% semilogy(G(:,2),G(:,6),'k--')
axis([0 20 1e-3 1e0]);
xlabel 'SNR (dB)'
ylabel 'BER'
title '(a) gm.10: \alpha = 0.9'

x = 1;
y = 0.5e-2;

text(x,y, ['.......  deterministic h'; ...
           '-.-.-. Gauss-Markov h   '; ...
           '\_\_\_\_ known channel  ']);

hold off
subplot(2,2,2)
load gm.07.txt
G=gm;
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
% semilogy(G(:,2),G(:,6),'k--')
axis([0 20 1e-3 1e0]);
xlabel 'SNR (dB)'
ylabel 'BER'
title '(b) gm.07: \alpha = 0.8'

hold off
subplot(2,2,3)
load gm.11.txt
G=gm;
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
% semilogy(G(:,2),G(:,6),'k--')
axis([0 20 1e-3 1e0]);
xlabel 'SNR (dB)'
ylabel 'BER'
title '(c) gm.11: \alpha = 0.7'

hold off
subplot(2,2,4)
load gm.12.txt
G=gm;
semilogy(G(:,2),G(:,3),'k:')
hold on
semilogy(G(:,2),G(:,4),'k-.')
semilogy(G(:,2),G(:,5),'k-')
% semilogy(G(:,2),G(:,6),'k--')
axis([0 20 1e-3 1e0]);
xlabel 'SNR (dB)'
ylabel 'BER'
title '(d) gm.12: \alpha = 0.6'

print -deps gm4a.eps
