Untitled

% 1.3

% 1. Load synchronization data from pre stmpSync.dat
syncData = load('pre_stmpSync.dat');
sensorTimeSync = syncData(:, 2); % Sensor timestamps from sync points
gatewayTimeSync = syncData(:, 1); % Gateway timestamps from sync points

% 2. Plot sync points on two subplots
figure;

% 2.1 First subplot: Gateway time vs Sensor time (scatter plot)
subplot(2,1,1);
plot(sensorTimeSync, gatewayTimeSync, '-');
xlabel('Szenzor lokális ideje');
ylabel('Gateway lokális ideje');
title('Szinkronizációs pontok: Szenzor vs Gateway lokális ideje');
grid on;

% 2.2 Second subplot: Plot sync points with the same timescale (index-based)
subplot(2,1,2);
plot(sensorTimeSync, 'r', 'DisplayName', 'Szenzor lokális ideje');
hold on;
plot(gatewayTimeSync, 'b', 'DisplayName', 'Gateway lokális ideje');
xlabel('Index');
ylabel('Idő (s)');
legend;
title('Szinkronizációs pontok egy időskálán');
grid on;

% 3. Fit a linear model to find transformation parameters
p = polyfit(sensorTimeSync, gatewayTimeSync, 1);
slope = p(1);
offset = p(2);

% 4. Load signal data from pre_gwySampFull.dat and pre_sensorDataFull.dat
gatewayData = load('pre_gwySampFull.dat');
sensorData = load('pre_sensorDataFull.dat');

gatewayTime = gatewayData(:,1); % Time from gateway samples
gatewaySignal = gatewayData(:,2); % Signal from gateway samples

sensorTime = sensorData(:,1); % Time from sensor samples
sensorSignal = sensorData(:,2); % Signal from sensor samples

% 5. Apply timestamp transformation to sensor time
transformedSensorTime = slope * sensorTime + offset;

% 6. Plot both the gateway and sensor signals before and after transformation
figure;

% 6.1 Plot original sensor signal vs gateway signal
subplot(2,1,1);
plot(gatewayTime, gatewaySignal, 'b', 'DisplayName', 'Gateway Signal');
hold on;
plot(sensorTime, sensorSignal, 'r', 'DisplayName', 'Sensor Signal (Original)');
xlabel('Time (s)');
ylabel('Signal');
legend;
title('Original Sensor vs Gateway Signals');
grid on;

% 6.2 Plot transformed sensor signal vs gateway signal
subplot(2,1,2);
plot(gatewayTime, gatewaySignal, 'b', 'DisplayName', 'Gateway Signal');
hold on;
plot(transformedSensorTime, sensorSignal, 'r', 'DisplayName', 'Sensor Signal (Transformed)');
xlabel('Time (s)');
ylabel('Signal');
legend;
title('Transformed Sensor vs Gateway Signals');
grid on;