%% VP - T3 - 3-Faktoren-Plan
% Authoren: andreas.voigt@ovgu.de
%% Aufräumen
clear all;close all;clc;
%% Daten
T = [120 140 120 140 120 140 120 140];
Z = [ 2 2 4 4 2 2 4 4];
K = [0.1 0.1 0.1 0.1 0.5 0.5 0.5 0.5];
M1 = [52.8 61.5 56.7 67.9 53.6 62.2 56.5 68.5];
M2 = [54.1 61.1 55.2 70.2 54.1 62.9 54.6 67.2];
%% Mittelwert und Fehler
M = zeros(1,8); 
for i=1:8
M(i) = (M1(i)+M2(i))/2; % einfacher mittelwert 
end
%% Berechnen von Effekten
ET = abs((M(1)-M(2)+M(3)-M(4)+M(5)-M(6)+M(7)-M(8))/4); % Temperatureffekt
EK = abs((M(1)-M(5)+M(2)-M(6)+M(3)-M(7)+M(4)-M(8))/4); % Katalysatoreffekt
%% Grafik
figure(1)
subplot(3,1,1)
plot(T,M,'rx-')
xlabel(' Temperatur ')
ylabel(' Ausbeute ')
text(122,65,[' Effekt ',num2str(ET,'%6.2f'),' % '])
subplot(3,1,2)
plot(Z,M,'gx-')
xlabel(' Zeitdauer ')
ylabel(' Ausbeute ')
subplot(3,1,3)
plot(K,M,'bx-')
xlabel(' Katalysator ')
ylabel(' Ausbeute ')
text(0.2,55,[' Effekt ',num2str(EK,'%6.2f'),' % '])
