67 lines
1.4 KiB
Matlab
67 lines
1.4 KiB
Matlab
close all;
|
|
events= readtable("001rec_SI.csv");
|
|
result = readmatrix("001rec_ADC_data.csv");
|
|
%result = result(14430:end, : );
|
|
result(:,14) = result(:,14)./2;
|
|
time = result(:,1)./1000;
|
|
|
|
temp = 0.000111 .* result(:,4:9) + 2.31991;
|
|
|
|
pressure = result(:,10:17).*(0.000015)-6.565;
|
|
|
|
load_cell = result(:,2:3);
|
|
weight = -1.166759307845543e-04 .* 0.5.*(load_cell(:,1) + load_cell(:,2)) + 4.971416323340051e+02;
|
|
|
|
figure()
|
|
scatter(time,weight,10)
|
|
title("load cell in Kg")
|
|
xlabel("Time in s")
|
|
ylabel("KG")
|
|
xlim([878.207 920.111+5.000])
|
|
|
|
|
|
pressureNames = { ...
|
|
'LOX tank pressure', ...
|
|
'Ethanol tank pressure', ...
|
|
'Ethanol pressure in injection plate', ...
|
|
'Chamber pressure', ...
|
|
'Ethanol pressure at inlet', ...
|
|
'Pressure channel 6', ...
|
|
'Pressure channel 7'};
|
|
|
|
for x = 1:7
|
|
figure()
|
|
hold on;
|
|
|
|
% Left axis → weight/thrust
|
|
yyaxis left
|
|
scatter(time, weight, 'DisplayName', 'Thrust');
|
|
ylabel("Thrust / Weight")
|
|
|
|
% Right axis → pressure
|
|
yyaxis right
|
|
scatter(time, pressure(:,x), 10, 'DisplayName', 'Pressure');
|
|
ylabel("Pressure")
|
|
|
|
% Title using your descriptive spreadsheet names
|
|
title(pressureNames{x})
|
|
|
|
xlabel("Time in s")
|
|
legend()
|
|
xlim([878.207 920.111 + 5.000])
|
|
end
|
|
|
|
%{
|
|
for x = 1:6
|
|
figure()
|
|
scatter(time,temp(:,x),10)
|
|
title(sprintf("temperature channel %d",x))
|
|
xlabel("Time in s")
|
|
ylabel("Temprature in C")
|
|
xlim([878.207 920.111+5.000])
|
|
end
|
|
|
|
|
|
%}
|
|
|