% Capacity vs. arrival rate given QoS a = 0:0.01:1; %arrival rate [p/us] la = length(a); z = 10; %longest time to wait [us] L = 1; % average packet length [kb] C = 0; % link speed [Gbps or kbpus] Pz_obj = [0.01 0.1 0.2]; lPz = length(Pz_obj); LS = zeros(la,lPz); % variable to store C(lambda) for j=1:lPz for l = 1:la c = C; lambda = a(l); Pz = Wait(c,lambda,L,z); while Pz > Pz_obj(j) c = c + 0.001; Pz = Wait(c,lambda,L,z); end LS(l,j) = c; end end figure f = plot(a,LS); f(1).LineWidth = 1.5; f(2).LineWidth = 1.5; f(3).LineWidth = 1.5; grid on title('Capacity vs. Arrival rate') xlabel('Arrival rate (\lambda) [packets/\mus]') ylabel({'Link speed (C)';'[Gbps]'}) legend({['Pz_{obj} = ',num2str(Pz_obj(1))],['Pz_{obj} = ',num2str(Pz_obj(2))],... ['Pz_{obj} = ',num2str(Pz_obj(3))]})