clear all ; close all ; clc ;

% 시뮬레이션을 통한 결과 도출
T = 10000; % 타임슬롯 T를 10000으로 설정
N = 100; % 노드 N을 100으로 설정
A = rand(N,T); % N행 T열에 0~1까지의 랜덤값들을 생성
S1 = zeros(21,T); % 모든 노드에서 한 슬롯 타임에 전송되는 모든 프레임을 카운트할 변수 설정
S2 = zeros(1,21); % 전송에 성공한 프레임을 카운트할 변수 설정
for pp = 1 : 21 % 각 노드에서 프레임 전송이 발생할 확률을 0%~20%까지 변화
    for tt = 1 : T % 타임슬롯을 변화
        for nn = 1 : N % 각 노드를 변화
            if A(nn,tt) <= (pp-1)/(N*10) % 시간의 변화에 따라 각 노드들에서 프레임 전송이 발생할 확률에 따른 프레임 발생을 감지
                S1(pp,tt) = S1(pp,tt)+1; % 한번 감지될때마다 s1의 값을 1씩 증가
            end
        end
        if S1(pp,tt) == 1 % 각 슬롯에 충돌이 일어났는지를 판별
            S2(pp) = S2(pp)+1; % 충돌없이 전송될 경우 1씩 카운트
        end
    end
end
S_F = S2./T; % throughput을 계산

% 이론적인 결과 도출
G = [0:10:200]/100; % (offered load) G = N p
S = G.*exp(-G); % throughput = s / T

% 그래프 설정
figure;
plot(G,S_F,'-bo');
hold on;
plot(G,S,'-k');
hold off;
title('Slotted - ALOHA ');
xlabel('G = offered load = Np');
ylabel('S = Throughput');
legend('Simulation','Theory');

-----------------------------------------------------------------------------------------------------------------------------------------------------------------------

공식을 통해서 구한 이론값과 시뮬레이션 값을 통한 그래프를 그려서 비교해봤습니다.