% plots for various N
%
hold off
subplot(2,2,1)
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 '(a) gm.07: N = 10'

x = 1;
y = 0.5e-2;

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

hold off
subplot(2,2,2)
load gm.15.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.15: N = 20'

hold off
subplot(2,2,3)
load gm.16.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.16: N = 30'

hold off
subplot(2,2,4)
load gm.17.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.17: N = 40'

print -deps gmN4.eps
