%Uebungsbeispiel zur Bluckdruck Messwerterfassung
%Aufgabe: Blutdruck.m
%Christian Widmer, 11.06.2002

clear;

klassen = 8;



%Messwerte

Atacand_16 = [
    126 89 82;
    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;
    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;
    124 87 68;
    132 79 64;
    124 84 65;
    148 92 68;
    141 90 69;
    128 89 68;
    140 86 68;
    136 91 68];


summe = 29;                     %Anzahl einträge (Zeilen)

%Berechnungen
MaxWerte=max(Atacand_16);
MinWerte=min(Atacand_16);
Mittelwert = mean(Atacand_16);
stdAbw = std(Atacand_16);

deltaX=(MaxWerte-MinWerte)/klassen;
Flaeche=deltaX*summe;


%Berechnung Systolischer Druck
x=MinWerte(:,1):0.1:MaxWerte(:,1);
u1=(x-Mittelwert(:,1))/stdAbw(:,1);
SysDruck=1/(stdAbw(:,1)*sqrt(2*pi))*exp(-0.5*u1.*u1);
FlaecheGauss1 = 0.1*sum(SysDruck);
Gauss1 = SysDruck*Flaeche(1)/FlaecheGauss1;

%Berechnung Diastolischer Druck
y=MinWerte(:,2):0.1:MaxWerte(:,2);
u2=(y-Mittelwert(:,2))/stdAbw(:,2);
DiasDruck=1/(stdAbw(:,2)*sqrt(2*pi))*exp(-0.5*u2.*u2);
FlaecheGauss2 = 0.1*sum(DiasDruck);
Gauss2 = DiasDruck*Flaeche(2)/FlaecheGauss2;

%Berechnung Puls
z=MinWerte(:,3):0.1:MaxWerte(:,3);
u3=(z-Mittelwert(:,3))/stdAbw(:,3);
Puls=1/(stdAbw(:,3)*sqrt(2*pi))*exp(-0.5*u3.*u3);
FlaecheGauss3 = 0.1*sum(Puls);
Gauss3 = Puls*Flaeche(3)/FlaecheGauss3;

 
%Ausgabe

% Ausgabe in Figure 1

subplot(3,1,1)
hold on;
title ('Systolischer Druck');
hist(Atacand_16(:,1),klassen)
plot(x,Gauss1, 'r')
text (138,9, ['Mittelwert = ',num2str(Mittelwert(:,1))]);
text (138,6, ['Standartabw. = ',num2str(stdAbw(:,1))]);



% Ausgabe in Figure 2

subplot(3,1,2)
hold on;
title ('Diastolischer Druck');
hist(Atacand_16(:,2),klassen)
plot(y,Gauss2, 'r')
text (79,9, ['Mittelwert = ',num2str(Mittelwert(:,2))]);
text (79,6, ['Standartabw. = ',num2str(stdAbw(:,2))]);


% Ausgabe in Figure 3

subplot(3,1,3)
hold on;
title ('Puls');
hist(Atacand_16(:,3),klassen)
plot(z,Gauss3, 'r')
text (57,7, ['Mittelwert = ',num2str(Mittelwert(:,3))]);
text (76,7, ['Standartabw. = ',num2str(stdAbw(:,3))]);

