%%%%%%%%%%%%%%%%%%%%%% Mean Rate 2 %%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%  Assumptions: 
%%%%%  --- n = even
%%%%%  --- Linear Time Scales: tau_{2k-1}=tau_{2k}=tau_{1,2}+deltatau (k-1)
%%%%%  --- average symbol time T_n=\tau_{1,2}(4+(alpha-1)(n-2))/4
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
D=20; % discrete points for curve resolution 
ainit=1.2;aend=3.5;ainc=(aend-ainit)/D;  % Range of parameter alpha
alpha=ainit:ainc:aend;
%s=a-1; % s=deltatau/tau_1
tau_1=1;% Normalization. Value is not need for simulation. 
N=10; % number of channels included for comparison
figure(201)
hold off
for n=2:2:N
    plot(alpha, 4*log(n)./(4+(alpha-1).*(n-2))/log(2),'linewidth',1)
    hold on
    %R_n=[R_n;2*log(n)./(2+(alpha-1).*(n-1))];
end
legend('\itn = 2','\itn = 4','\itn = 6','\itn = 8','\itn = 10')
%axis([ainit aend min([Ceven(12/2,:) Ceven(4/2,:) Ceven(6/2,:) Ceven(8/2,:)]) max([Ceven(2/2,:) Ceven(4/2,:) Ceven(6/2,:) Ceven(8/2,:)])])
xlabel('\alpha = \tau_{3,4} / \tau_{1,2}','fontsize',18)
ylabel('\tau_{1,2}R_{n}','fontsize',18)
grid on
