% HS 2009
% Milestone 3
% 
% Anastasiades Bryan
% Schneider Melanie
% Weber Marco
% -------------------------------
% -------------------------------
% 
% Plots Kombiniert mit hydro plant
% 
% -------------------------------
% -------------------------------

clc

 if (exist('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;

if (exist('daten')==0)
   
    sim('hydroplant_MS3_091110');
end

ind.h_W=1;
ind.V_flw_T=2;
ind.V_flw_F=3;
ind.p_F=4;
ind.omega=5;
ind.U=6;
ind.I=7;
ind.P=8;
ind.x=9;


% Plots wie in HydroPlant angeordnet 3:3

figure('name','Milestone 3- KOMBINATION')


% Plot Ventilposition
subplot(3,3,1)
plot(sim_daten.time, sim_daten.signals.values(:,gem_x),'LineWidth',2);
hold on
plot(daten.time, daten.signals.values(:,ind.x),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.h_W),'r');
grid on
xlabel('Zeit [s]');
ylabel('Wasserstand [m]');
title({'Milestone 3: PLOT-VERGLEICH';...
    'Anastasiades Bryan / Schneider Melanie / Weber Marco';...
    '';'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);
hold on
plot(daten.time, daten.signals.values(:,ind.p_F),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.V_flw_T),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.V_flw_F),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.omega),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.U),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.I),'r');
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);
hold on
plot(daten.time, daten.signals.values(:,ind.P),'r');
grid on
xlabel('Zeit [s]');
ylabel('Leistung [W]');
title('Leistung');
legend('US','Hydro Plant')