Mein MATLAB Forum

%%Berechnung des polytropischen Exponent der Messungen

clc;

Y_Time = S_druck(:,1);

%% get Variable from Workspace
%% Saugtemparatur

T1.D = S_temperaturD(:,2);
T1.C = S_temperaturC(:,2);
T1.B = S_temperaturB(:,2);
T1.A = S_temperaturA(:,2);

%% Endtemparatur

T2.D = E_temperaturD(:,2);
T2.C = E_temperaturC(:,1);
T2.B = E_temperaturB(:,2);
T2.A = E_temperaturA(:,2);


%% Saugdruck

P1.D = S_druck(:,5);
P1.C = S_druck(:,4);
P1.B = S_druck(:,3);
P1.A = S_druck(:,2);

%% Enddruck

P2.D = E_druck(:,5);
P2.C = E_druck(:,4);
P2.B = E_druck(:,3);
P2.A = E_druck(:,2);

%%% Funktion zur Berechnung des polytropischen Exponent

P_A = log(P2.A/P1.A);
P_B = log(P2.B/P1.B);
P_C = log(P2.C/P1.C);
P_D = log(P2.D/P1.D);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

T_A = log(T2.A/T1.A);
T_B = log(T2.B/T1.B);
T_C = log(T2.C/T1.C);
T_D = log(T2.D/T1.D);

%%______________________________________________________________

N_A = P_A/(P_A-T_A);
N_B = P_B/(P_B-T_B);
N_C = P_C/(P_C-T_C);
N_D = P_D/(P_D-T_D);

%Create figure
H.fig = figure;
H.axs = axes;
hold on;
grid on;
box on;

%Plot N_i

plot(N_A, Y_Time, 'r-');
plot(N_B, Y_Time, 'r--');
plot(N_C, Y_Time, 'b-');
plot(N_D, Y_Time, 'k.');

hodl off