%Ersteller:     D'Agostino Peter
%Aufgabe:       Blutdruck
%Dateiname:     blutdruck.m
%-------------------------------------------------------

clear;clf;
% Geg: Reihenfolge der Werte : systolischer Druck, diastolischer Druck und Puls
% Das Medikament Atacand_16 wurde verschrieben.

Atacand_16 = [  126 89 82;      %13.05.01
                138 86 72;
                147 92 77;
                146 93 78;
                128 94 68;
                145 88 73;
                130 92 66;
                132 89 64;
                131 89 70;
                135 96 66;
                132 87 75;
                136 96 75;      %06.06.01
                131 94 68;
                133 94 59;
                132 95 67;
                136 88 75;
                132 92 63;
                133 92 75;
                136 94 67;
                139 83 74;
                132 87 72;      %06.07.01
                124 87 68;
                132 79 64;      %03.08.01
                124 84 65;
                148 92 68;
                141 90 69;
                128 89 68;
                140 86 68;
                136 91 68];     %20.08.01
            
 
    
     Mittelwerte   = mean(Atacand_16)  % Berechnen des Mittelwerte
     MaxWerte      = max(Atacand_16)   % Berechnen des Maximalwerte      
     MinWerte      = min(Atacand_16)   % Berechnen des Minimalwerte 
     StdAbweichung = std(Atacand_16)    % Berechnen des Standardabweichungen
     summe = 29;
     
     
     AnzahlKlassen = 8;                             % Anzahl Klassen definieren Wurzel aus Werte                
     deltaX=(MaxWerte-MinWerte)/AnzahlKlassen;
     Flaeche= deltaX*summe;
     SystolischerDruck = Atacand_16(:,1);
     stdabSystDruck =std(SystolischerDruck);        % Standardabweichung
     MittelWertSystDruck = mean(Atacand_16(:,1));
     MittelWertDiastDruck = mean(Atacand_16(:,2));    
     MittelWertPuls = mean(Atacand_16(:,3));
     x=MinWerte(1):0.1:MaxWerte(1);                 % Werte Bereich x
     
     u=(x-MittelWertSystDruck)/stdabSystDruck;      % u Einsetzen in Formel Normalverteilung
     
     
     % Berechnen Normalverteilung 
     NormSystDruck = 1/(sqrt(2*pi*stdabSystDruck))*exp(-0.5*(u.*u));
     FlaecheGauss = 0.1*sum(NormSystDruck);
     Gauss = NormSystDruck*Flaeche(1)/FlaecheGauss;
     
     subplot(3,1,1);
     hold on; 
     % Häufigkeitsverteilung der Werte systolischer Druck
              
     hist(SystolischerDruck,AnzahlKlassen);     % Erstellen Histogramm
     plot(x,Gauss,'b','LineWidth',2);     
     h = findobj(gca,'Type','patch');
     set(h,'FaceColor','r','EdgeColor','w');    % Histogramm in der Farbe rot mit weissen Rand
     sa = sprintf ('%5.4f',stdabSystDruck);
     m = sprintf ('%5.4f' ,MittelWertSystDruck);
     text(135,7,['Standardabw. = ',sa],'FontSize',12)
     text(135,9,['Mittelwert = ',m],'FontSize',12)
     grid;
     % Ende Aufzeichnung Systolischer Druck
     % -----------------------------------------------------------------------------
     
     subplot(3,1,2);
     hold on;
     x=MinWerte(2):0.1:MaxWerte(2);  
     DiastolischerDruck = Atacand_16(:,2)             % Häufigkeitsverteilung der Werte systolischer Druck
     stdabDiastDruck = std(DiastolischerDruck);
     u=(x-MittelWertDiastDruck)/stdabDiastDruck;      % u Einsetzen in Formel Normalverteilung
     % Berechnen Normalverteilung 
     NormDiastDruck = 1/(sqrt(2*pi*stdabDiastDruck))*exp(-0.5*(u.*u));
     FlaecheGauss = 0.1*sum(NormDiastDruck);
     Gauss = NormDiastDruck*Flaeche(2)/FlaecheGauss;
     hist(DiastolischerDruck,AnzahlKlassen);
     plot(x,Gauss,'r','LineWidth',2);  
     h = findobj(gca,'Type','patch');
     set(h,'FaceColor','g','EdgeColor','w')     % Histogramm in der Farbe gruen mit blauen Rand
     stdDD = sprintf ('%5.4f',stdabDiastDruck);
     m = sprintf ('%5.4f' ,MittelWertDiastDruck);
     text(78,5,['Standardabw. = ',stdDD],'FontSize',12)
     text(85,8,['Mittelwert = ',m],'FontSize',12)
     grid;
     % Ende Aufzeichnung Diastolischer Druck
     % -----------------------------------------------------------------------------   
     
     
     subplot(3,1,3);
     hold on;
     x=MinWerte(3):0.1:MaxWerte(3);  
     Puls = Atacand_16(:,3)                      % Häufigkeitsverteilung der Werte Puls
     stdabPuls = std(Puls);
     u=(x-MittelWertPuls)/stdabPuls; % u Einsetzen in Formel Normalverteilung
     % Berechnen Normalverteilung 
     NormPuls = 1/(sqrt(2*pi*stdabPuls))*exp(-0.5*(u.*u));
     FlaecheGauss = 0.1*sum(NormPuls);
     Gauss = NormPuls*Flaeche(3)/FlaecheGauss;
     hist(Puls,AnzahlKlassen);
     plot(x,Gauss,'g','LineWidth',2);  
     h = findobj(gca,'Type','patch');
     set(h,'FaceColor','c','EdgeColor','w'); % Histogramm in der Farbe gruen mit blauen Rand
     stdP = sprintf ('%5.4f',stdabPuls);
     m =sprintf ('%5.4f' ,MittelWertPuls);
     text(55,7,['Standardabw. = ',stdP],'FontSize',12)
     text(75,7,['Mittelwert = ',m],'FontSize',12)
     grid;
     % Ende Aufzeichnung Puls
     % -----------------------------------------------------------------------------   