clc;clear;
close all;
% Beleg_Nr = 3; % Christoph
Beleg_Nr = 5; % Zoe
% Beleg_Nr = 7; % Samuel
% Beleg_Nr = 9; % Sebastian
if Beleg_Nr == 3;
    %Parameter Christoph
    z_p = 1;
    R = 0.18;
    L = 0.35e-3;
    Psi = 0.1528;
    J = 0.25;
    m_l01 = 0;
    m_l02 = 1;
    I = 5;
    Phi = 0;
    Phi01 = 0;
    omega01 = 10;
    t_s_m = 15;%nicht relevant
    t_s_j = 3.5;
    Phi02 = 0;
    omega02 = -20;
    t_end1 = 5;
    t_end2 = 10;
elseif Beleg_Nr == 5;
    %Parameter Zoe
    z_p = 1;
    R = 0.07;
    L = 0.1e-3;
    Psi = 0.1528;
    J = 0.6;
    m_l01 = 2;
    m_l02 = 2;
    I = 10;
    Phi = 0;
    Phi01 = 0;
    omega01 = 0;
    t_s_m = 3;
    t_s_j = 15;%nicht relevant
    Phi02 = 0;
    omega02 = -40;
    t_end1 = 5;
    t_end2 = 10;
elseif Beleg_Nr == 7;
    %Parameter Samuel
    z_p = 1;
    R = 0.18;
    L = 0.35e-3;
    Psi = 0.1528;
    J = 0.2;
    m_l01 = 0;
    m_l02 = 1;
    I = 4;
    Phi = 0;
    Phi01 = 0;
    omega01 = -30;
    t_s_m = 15;%nicht relevant
    t_s_j = 3;
    Phi02 = 0;
    omega02 = 10;
    t_end1 = 5;
    t_end2 = 10;
elseif Beleg_Nr == 9;
    %Parameter Sebastian
    z_p = 1;
    R = 0.07;
    L = 0.1e-3;
    Psi = 0.1528;
    J = 0.6;
    m_l01 = 0;
    m_l02 = 2;
    I = 10;
    Phi = 0;
    Phi01 = 0;
    omega01 = 0;
    t_s_m = 15;%nicht relevant
    t_s_j = 3;
    Phi02 = 0;
    omega02 = -40;
    t_end1 = 5;
    t_end2 = 10;
end

%% Aufgabe 2 Zoe/Sebastian
sim('mdl_a2_20160513');

% Plots Aufgabe 2
% 1.Winkelgeschwindigkeit
figure;
plot(omega_out(:,1),omega_out(:,2),'linewidth',2)
set(gca,'FontSize',16);
grid on;
title('Winkelgeschwindigkeit','FontSize',20)
xlabel('Zeit in s','FontSize',14)
ylabel('Winkelgeschwindigkeit in 1/s','FontSize',14)
fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe2_1.csv'];
csvwrite(fname,omega_out);

%Spannungsortskurve
figure;
plot(U1_out(:,2),U2_out(:,2),'linewidth',2)
set(gca,'FontSize',16);
grid on;
title('Spannungsortskurve','FontSize',20)
xlabel('','FontSize',14)
ylabel('','FontSize',14)
fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe2_2.csv'];
csvwrite(fname,[U1_out(:,2),U2_out(:,2)]);

%Mechanische/elektrische Leistung
fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe2_3.csv'];
if (Beleg_Nr == 3)||(Beleg_Nr == 7)
    figure;
    plot(Pmech_out(:,1),Pmech_out(:,2),'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Mechanische Leistung','FontSize',20)
    xlabel('Zeit in s','FontSize',14)
    ylabel('Leistung in W','FontSize',14)
    csvwrite(fname,Pmech_out);
else
    figure;
    plot(Pel_out(:,1),Pel_out(:,2),'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Elektrische Leistung','FontSize',20)
    xlabel('Zeit in s','FontSize',14)
    ylabel('Leistung in kW','FontSize',14);
    csvwrite(fname,Pel_out);
end

%% Aufgabe 3
sim('mdl_a3_20160519')

% Plots Aufgabe 3
% 1.Winkelgeschwindigkeit
figure;
plot(a3_t,a3_omega,'linewidth',2)
set(gca,'FontSize',16);
grid on;
title('Winkelgeschwindigkeit','FontSize',20)
xlabel('Zeit in s','FontSize',14)
ylabel('Winkelgeschwindigkeit in 1/s','FontSize',14)
fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe3_1.csv'];
csvwrite(fname,[a3_t,a3_omega]);

%Stromortskurve
figure;
if (Beleg_Nr == 7)||(Beleg_Nr == 9)
    plot(a3_t,a3_I12(:,1),a3_t,a3_I12(:,2),'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Strangströme','FontSize',20)
    xlabel('Strom in A','FontSize',14)
    ylabel('Zeit in s','FontSize',14)
    fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe3_2.csv'];
    csvwrite(fname,[a3_t, a3_I12(:,1),a3_I12(:,2)]);
else
    plot(a3_I12(:,1),a3_I12(:,2),'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Stromortskurve','FontSize',20)
    xlabel('','FontSize',14)
    ylabel('','FontSize',14)
    fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe3_2.csv'];
    csvwrite(fname,[a3_I12(:,1),a3_I12(:,2)]);
end

%Mechanische/elektrische Leistung
fname = ['Beleg' num2str(Beleg_Nr) '_Aufgabe3_3.csv'];
if (Beleg_Nr == 3)||(Beleg_Nr == 7)
    figure;
    plot(a3_t,a3_mi,'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Inneres Moment','FontSize',20)
    xlabel('Zeit in s','FontSize',14)
    ylabel('Inneres Moment in Nm','FontSize',14)
    csvwrite(fname,[a3_t, a3_mi]);
else
    a3_Pmech = a3_mi.*a3_omega;
    figure;
    plot(a3_t,a3_Pmech,'linewidth',2)
    set(gca,'FontSize',16);
    grid on;
    title('Mechanische Leistung','FontSize',20)
    xlabel('Zeit in s','FontSize',14)
    ylabel('Leistung in W','FontSize',14);
    csvwrite(fname,[a3_t,a3_Pmech]);
end