5 files changed, 137 insertions, 0 deletions

diff --git a/analysis/e1anl/src/__pycache__/agregate.cpython-313.pyc b/analysis/e1anl/src/__pycache__/agregate.cpython-313.pyc
new file mode 100644
index 0000000..3a0d6cb
--- /dev/null
+++ b/analysis/e1anl/src/__pycache__/agregate.cpython-313.pyc
diff --git a/analysis/e1anl/src/__pycache__/analyze.cpython-313.pyc b/analysis/e1anl/src/__pycache__/analyze.cpython-313.pyc
new file mode 100644
index 0000000..a891763
--- /dev/null
+++ b/analysis/e1anl/src/__pycache__/analyze.cpython-313.pyc
diff --git a/analysis/e1anl/src/agregate.py b/analysis/e1anl/src/agregate.py
new file mode 100644
index 0000000..1054288
--- /dev/null
+++ b/analysis/e1anl/src/agregate.py
@@ -0,0 +1,71 @@
+# mqtt_receiver.py
+
+import paho.mqtt.client as mqtt
+import csv
+import time
+from datetime import datetime
+from collections import deque
+import os
+os.makedirs("output", exist_ok=True)
+
+BROKER = "192.168.1.101"
+SUB_TOPIC = "device/echo/in"
+CSV_FILE = "output/rtt_throughput_log.csv"
+WINDOW_SEC = 1
+
+msg_times = deque()
+should_stop = False  # flag to stop loop externally
+
+
+def on_connect(client, userdata, flags, rc):
+    if rc == 0:
+        print("Connected to broker")
+        client.subscribe(SUB_TOPIC)
+    else:
+        print("Failed to connect, return code:", rc)
+
+
+def on_message(client, userdata, msg):
+    try:
+        payload_str = msg.payload.decode()
+        sent_ms = int(payload_str)
+        received_ms = int(time.time() * 1000)
+        rtt = received_ms - sent_ms
+        timestamp = datetime.now().isoformat()
+
+        now_sec = time.time()
+        msg_times.append(now_sec)
+        while msg_times and msg_times[0] < now_sec - WINDOW_SEC:
+            msg_times.popleft()
+        throughput = len(msg_times) / WINDOW_SEC
+
+        with open(CSV_FILE, "a", newline="") as f:
+            writer = csv.writer(f)
+            writer.writerow([timestamp, sent_ms, received_ms, rtt, throughput])
+
+    except Exception as e:
+        print("Error processing message:", e)
+
+
+def collect_data(duration_sec):
+    global should_stop
+    should_stop = False
+
+    # Start CSV log
+    with open(CSV_FILE, "w", newline="") as f:
+        writer = csv.writer(f)
+        writer.writerow(["timestamp", "sent_ms", "received_ms", "rtt_ms", "throughput_msg_per_s"])
+
+    client = mqtt.Client()
+    client.on_connect = on_connect
+    client.on_message = on_message
+    client.connect(BROKER, 1883, 60)
+
+    client.loop_start()
+    print(f"Zbieranie danych przez {duration_sec / 60:.1f} minut...")
+
+    time.sleep(duration_sec)
+
+    print("Zbieranie danych zakończone.")
+    client.loop_stop()
+    client.disconnect()
diff --git a/analysis/e1anl/src/analyze.py b/analysis/e1anl/src/analyze.py
new file mode 100644
index 0000000..6e73df3
--- /dev/null
+++ b/analysis/e1anl/src/analyze.py
@@ -0,0 +1,49 @@
+# analyze_metrics.py
+
+import pandas as pd
+import matplotlib.pyplot as plt
+import seaborn as sns
+
+import os
+os.makedirs("output", exist_ok=True)
+
+def analyze():
+    CSV_FILE = "output/rtt_throughput_log.csv"
+    df = pd.read_csv(CSV_FILE)
+
+    df['sent_ms'] = df['sent_ms'].astype(int)
+    df['received_ms'] = df['received_ms'].astype(int)
+    df['rtt_ms'] = df['rtt_ms'].astype(int)
+    df['throughput_msg_per_s'] = df['throughput_msg_per_s'].astype(float)
+
+    mean_rtt = df['rtt_ms'].mean()
+    p95_rtt = df['rtt_ms'].quantile(0.95)
+    mean_throughput = df['throughput_msg_per_s'].mean()
+
+    print(f"Średnie RTT: {mean_rtt:.2f} ms")
+    print(f"RTT 95-percentyl: {p95_rtt:.2f} ms")
+    print(f"Średni throughput: {mean_throughput:.2f} msg/s")
+
+    plt.figure(figsize=(6,4))
+    sns.barplot(x=['Variant A'], y=[mean_rtt], palette="Set2")
+    plt.title("Średnie RTT - wariant A")
+    plt.ylabel("RTT [ms]")
+    plt.tight_layout()
+    plt.savefig("output/rtt_mean_a.png")
+    plt.show()
+
+    plt.figure(figsize=(6,4))
+    sns.barplot(x=['Variant A'], y=[p95_rtt], palette="Set3")
+    plt.title("RTT 95-percentyl - wariant A")
+    plt.ylabel("RTT [ms]")
+    plt.tight_layout()
+    plt.savefig("output/rtt_p95_a.png")
+    plt.show()
+
+    plt.figure(figsize=(6,4))
+    sns.barplot(x=['Variant A'], y=[mean_throughput], palette="Set1")
+    plt.title("Średni throughput - wariant A")
+    plt.ylabel("Messages per second")
+    plt.tight_layout()
+    plt.savefig("output/throughput_a.png")
+    plt.show()
diff --git a/analysis/e1anl/src/main.py b/analysis/e1anl/src/main.py
new file mode 100644
index 0000000..b35eead
--- /dev/null
+++ b/analysis/e1anl/src/main.py
@@ -0,0 +1,17 @@
+# main.py
+
+from agregate import collect_data
+from analyze import analyze
+
+def main():
+    DURATION_HOURS = 2
+    DURATION_SECONDS = DURATION_HOURS * 3600
+
+    print("Collecting data...")
+    collect_data(DURATION_SECONDS)
+
+    print("Analyzing data...")
+    analyze()
+
+if __name__ == "__main__":
+    main()