% Systemmodellierung 
% HS 2009
% Milestone 3
% 
% Anastasiades Bryan
% Schneider Melanie
% Weber Marco
% -------------------------------
% -------------------------------
% 
% Plots
% 
% -------------------------------
% -------------------------------

close all

if (exist('sim_daten')==0)
    sim('wasserkraftwerk');
end



gem_x=1;
gem_h_W=2;
gem_p_F=3;
gem_V_flw_T=4;
gem_V_flw_F=5;
gem_omega=6;
gem_U=7;
gem_I=8;
gem_P=9;



% Plots wie in HydroPlant angeordnet 3:3

figure('name','Milestone 3')

% Plot Ventilposition
subplot(3,3,1)
plot(sim_daten.time, sim_daten.signals.values(:,gem_x),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Position der Ventilnadel [-]');
title('Ventil Position');

% Plot Wasserstand im Wasserschloss
subplot(3,3,2)
plot(sim_daten.time, sim_daten.signals.values(:,gem_h_W),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Wasserstand [m]');
title('Wasserstand im Wasserschloss');

% Plot Druck vor dem Ventil
subplot(3,3,3)
plot(sim_daten.time, sim_daten.signals.values(:,gem_p_F),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Druck [Pa]');
title('Druck vor Ventil');

% Plot Volumenstrom Tunnel
subplot(3,3,4)
plot(sim_daten.time, sim_daten.signals.values(:,gem_V_flw_T),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Volumenstrom Tunnel [m^3/s]');
title('Volumenstrom Tunnel');

% Plot Volumenstrom Fallrohr
subplot(3,3,5)
plot(sim_daten.time, sim_daten.signals.values(:,gem_V_flw_F),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Volumenstrom Fallrohr [m^3/s]');
title('Volumenstrom Fallrohr');

% Plot Drehgeschwindigkeit der Turbine
subplot(3,3,6)
plot(sim_daten.time, sim_daten.signals.values(:,gem_omega),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('\omega [rad/s]');
title('Drehgeschwindigkeit der Turbine');

% Plot Spannung
subplot(3,3,7)
plot(sim_daten.time, sim_daten.signals.values(:,gem_U),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Spannung [V]');
title('Spannung');

% Plot Strom
subplot(3,3,8)
plot(sim_daten.time, sim_daten.signals.values(:,gem_I),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Strom [A]');
title('Strom');

% Plot Leistung
subplot(3,3,9)
plot(sim_daten.time, sim_daten.signals.values(:,gem_P),'LineWidth',2);
grid on
xlabel('Zeit [s]');
ylabel('Leistung [W]');
title('Leistung');