ARQ.zip_ARQ_zip

%% Definicion de variables aa = 8*[1:1:1024]; % instring; %Temporary input string % d; %Link distance in meters % R; %Link data rate in bits per second % L; %Frame length in bits % pr_bit; %Probability of a bit error % protocol_type; %Protocol type (1 = SAW, 2 = SW) % ARQ_type; %ARQ type for SW (1 = SR, 2 = GBN) % N; %Window size for SW % P; %Probability of a frame error % a; %Stallings "a" parameter % U; %Link utilization as a ratio V = 2.0e8; %Copper propagation in meters per second %% INICIALIZACION DE VARIABLES disp('-------------------------------------------- dlc.c -----'); disp('- Program to compute utilization for DLC protocols -'); disp('--------------------------------------------------------'); %Prompt for distance of link in meters instring = input('Length of link in meters =======================> '); d = instring; %Prompt for data rate of link in bits per second instring = input('Data rate of link in bits per second ===========> '); R = instring; %Prompt for size of frames in bytes (converted to bits for L) instring = input('Size of frames in bytes ========================> '); L = 8 * instring; %Prompt for error rate (BER) instring = input('Pr[bit error] (0.0 to 1.0) =====================> '); pr_bit = instring; %Prompt for protocol (SAW or SW) instring = input('Protocol type (1 = SAW, 2 = SW) ================> '); protocol_type = instring; %If SW then prompt for other parameters if (protocol_type == 2) %Prompt for window size instring = input('Window size ====================================> '); N = instring; %Prompt for error control (GBN or SR) instring = input('ARQ type (1 = SR, 2 = GBN) =====================> '); ARQ_type = instring; end for i=1:1024 L=aa(i); %Compute P (probability of frame error) P = 1 - (1 - pr_bit)^L; %Compute "a" value a = (R * d) / ( V * L); %Compute U for SAW if (protocol_type == 1) U(i) = (1 - P) / (1 + (2 * a)); end %Compute U for SW and SR if ((protocol_type == 2) && (ARQ_type == 1)) if (N >= ((2 * a) + 1)) U(i) = (1 - P); else U(i) = (N * (1 - P)) / ((2 * a) + 1); end end %Compute U for SW and GBN if ((protocol_type == 2) && (ARQ_type == 2)) if (N >= ((2 * a) + 1)) U(i) = (1 - P) / (1 + (2 * a * P)); else U(i) = (N * (1 - P)) / (((2 * a) + 1) * (1 - P + N * P)); end end %Output results % disp('-------------------------------------------------------- '); % disp('a = ') % disp(a) % disp('Pr[frame loss] = ') % disp(P) % disp('U = ') % disp(U) end U = U*100; plot (aa,U)