%-------------------------------------------------------------------------
%Plot über Frequenz für Parameterstudie 2 Massenschwinger
%Dominik Marquardt/ TI-F/ 15.11.2010
%zwei/drei Radmassen werden Simuliert
%-------------------------------------------------------------------------
%--------------------------LOAD-DATA--------------------------------------

close all;

%%
for x=1:1  % 3mal wird der Plotfile durchlaufen einmal für jede Anregung

    %Funktion zur Entscheidung welche Ergebnismatrix geladen wird. Z wird
    %im Steuerfile definiert.
    %     if z==6
%     [s_ca,s_auf,s_rad,s_er,a_auf,a_rad,f_ca,f_ka,F_N,v_er,leg,sww,tm,TS]=erg_matrix_2(x);

    %     elseif z==9
            [s_ca,s_auf,s_rad,s_er,a_auf,a_rad,f_ca,f_ka,F_N,v_er,leg,sww,tm,TS]=erg_matrix_3(x);

%
    a=3;            %Anfangswert für fortlaufende nummerierung der Figures

    %################## BILDE-BASISTITEL-JE-NACH-ANREGUNG ##################

    % x wird aus Steuerfile gelesen
    if x==1
        basistitel = ('Variation Radmasse Sweep');%bilde HAUPTTITEL

%     elseif x==2
%         basistitel= ('Variation Radmasse StepUp');%bilde HAUPTTITEL

%     elseif x==3
%         basistitel = ('Variation Radmasse FB');%bilde HAUPTTITEL
    end
    %-----------------------PLOT-FIGURES----------------------------------

%    %########################### Bodediagramme #############################
% 
%     for b=1:3  % Alle Bodediagramme haben den selben Subplot(212). Deshalb die Schleife
%         % mit a werden die wenigen Veränderungen in Subplot(211)
%         % definiert.
%         a = a+1;figure(a);
% 
%         if a==1
%             [f,Yamp,Ypha] = fft_ausein(s_auf, s_er, TS); %AUFBAUBEWEGUNG Bode Diagramme (Übertragungsverhalten)
%         elseif a==2
%             [f,Yamp,Ypha] = fft_ausein(s_rad, s_er, TS); %RADWEG-> Bode Diagramme (Übertragungsverhalten)
%         elseif a==3
%             [f,Yamp,Ypha] = fft_ausein(s_ca, s_er, TS); %FEDERWEG--> Bode Diagramme (Übertragungsverhalten)
%         end
% 
%         ax(1) = subplot(2,1,1);
%         semilogx(f,Yamp);hold on;box on;grid on;
%         set(gca,'Fontsize',16);
%         if a==1
%             title({basistitel,'Bode-Diagramm s auf/s er'},'Fontsize',18);
%         elseif a==2
%             title({basistitel,'Bode-Diagramm s rad / s er'},'Fontsize',18);
%         else
%             title({basistitel,'Bode-Diagramm s ca/ s er'},'Fontsize',18);
%         end
%         ylabel('Amplitudenverhältnis','FontSize',16);
%         legend(leg,'Location', 'NorthWest');
%         axis([-inf 25 -inf inf]);          %0.5/TS ist die Grenzfrequenz
% 
%         ax(2) = subplot(2,1,2);
%         semilogx(f,Ypha);hold on;box on;grid on;
%         set(gca,'Fontsize',16);
%         xlabel(gca,'Frequenz [Hz]','FontSize',16);
%         ylabel(gca,'Phase [Grad]','FontSize',16);
%         legend(leg,'Location', 'NorthWest');
%         axis([-inf 25 -360 360]);
%         linkaxes(ax,'x');
% %         set(figure(a),'Position', [0,0, 1280, 1024]); %Vergrößert figure auf screen size
%     end

    %###################### VERGRÖSSERUNGSFUNKTIONEN ######################
    for b=4:6
        a = a+1;figure(a);
            
        if a==4
            %Amplitudenspektrum der Anregung
            [f,Yamp] = fft_amp(s_er,TS);  %ANREGUNG        %fft_amp(Data,Abtastzeitintervall)
        elseif a==5
            %BODEDIAGRAMM
            [f,Yamp,Ypha] = fft_ausein(F_N,s_er,TS); %RADLAST
        elseif a==6
            %BODEDIAGRAMM
            [f,Yamp,Ypha] = fft_ausein(a_auf,s_er,TS);  %AUFBAUBESCHLEUNIGUNG
        end

        if (a==4)||(a==5)||(a==6)
            semilogx(f,Yamp);grid on;hold on;box on;
            % Ein Peak in der ersten Zeile ist der Grund, dass erst ab der
            % zweiten geplottet wird.
            set(gca,'Fontsize',16);
            if a==4
                title({basistitel,'Amplitudenspektrum von s er'},'FontSize',18);
                ylabel(gca,'Amplitude (fft/n)','FontSize',16);
                legend('Anregung','Location', 'NorthWest')
            elseif a==5
                title({basistitel,'Vergrößerungsfunktion der Radlast F N/ s er'},'FontSize',18);
                ylabel(gca,'Amplitudenverhältnis N/m','FontSize',16);
                legend(leg,'Location', 'NorthWest');
            else
                title({basistitel,'Vergrößerungsfunktion der Aufbaubeschleunigung a auf/ s er'},'FontSize',18);
                ylabel(gca,'Amplitudenverhältnis 1/s^2','FontSize',16);
                legend(leg,'Location', 'NorthWest');
            end
            xlabel(gca,'Frequenz [Hz]','FontSize',16);
%             axis([-inf 25 -inf inf]);
            XLIM([-inf 25]);
            
        end

%          %########################## LEISTUNGSDICHTE ########################
%  
%         %Leistungsdichtespektrum der Anregung
%         if a==7
%             [f,Ypsd] = fft_psd(s_er,TS);  %ANREGUNG   %fft_psd(Data,Abtastzeitintervall)
% 
%             %Leistungsdichtespektrum der Radlast
%         elseif a==8
%             [f,Ypsd] = fft_psd(F_N,TS);   %RADLAST   %fft_psd(Data,Abtastzeitintervall)
% 
%             %Leistungsdichtespektrum der Aufbaubeschleunigung
%         elseif a==9
%             [f,Ypsd] = fft_psd(a_auf,TS); %AUFBAUBESCHLEUNIGUNG   %fft_psd(Data,Abtastzeitintervall)
%         end
% 
%         if (a==7)||(a==8)||(a==9)       % a sorgt für die richtige Titelwahl
% 
%             loglog(f,Ypsd);grid on;hold on;box on;
%             set(gca,'Fontsize',16);
%             if a==7
%                 title({basistitel,'Leistungsdichtespektrum von s er'},'FontSize',18);
%                 legend('Anregung','Location', 'NorthWest');
%             elseif a==8
%                 title({basistitel,'Leistungsdichtespektrum der Radlast N^2/Hz'},'FontSize',18);
%                 legend(leg,'Location', 'NorthWest');
%             elseif a==9
%                 title({basistitel,'Leistungsdichtespektrum der Aufbaubeschleunigung m^2/s^4*Hz'},'FontSize',18);
%                 legend(leg,'Location', 'NorthWest');
%             end
%             xlabel(gca,'Frequenz [Hz]','FontSize',16);
%             ylabel(gca,'spektrale Dichte m^2/Hz','FontSize',16);
%             axis([-inf 25 -inf inf]);
%         end

         set(figure(a),'Position', [0,0, 1280, 1024]); %Vergrößert figure auf screen size
    end


         save_plots_freq;                % Files zum Speichern der Plots aufrufen
%         clear all;clc;clf;close all;
end