function eGlobal = berechneEGlobal(x)

f       =   [1,2,3,4,5,6,7,8,9,10,...
            20,30,40,50,60,70,80,90,100,200,300];
AMess05 =   [0 0 0 0 0 0 0 0 0 0 0 0 0.7 1 1.1 1.45 1.5 1.5 1.6 1.55 1.1 0 -1.5 ...
            -2.9 -10.5];
AMess25 =   [0 0 0 0 0 0 0 0 0 0 0 0 0.9 1.1 1 0.9 -0.1 -1.1 -1.9 -2.7...
            -4.6 -7 -9.5 -12.8 -22.9];
AMess50 =   [0 0 0 0 0 0 0 0 0 0 0 0 0.9 0 -1 -2.2 ...
            -3.2 -4.5 -5.5 -7 -9.5 -13.2 -16 -19 -29.5];
        
A = 1;      % Amplitude des Eingangssignals
abweichung1 = zeros(1,length(f));
amplitude1 = zeros(1,length(f));
AdB1     = zeros(1,length(f));
phase1 = zeros(1,length(f));

for i = 1:length(f)
    % Filtern --> in der Filterfunktion muss simuliert werden!!!
    [simulationszeit,eingangssignal,ausgangssignalGefiltert] = berechneGefiltertesSignal(f(i),x,A);
    
    % Amplitude und Phase bestimmen
    [amplitude1(i),phase1(i),~] = berechneAmpPhs(simulationszeit,f(i),eingangssignal,ausgangssignalGefiltert);
    AdB1(i) = 20*log10(amplitude1(i)/A);
    abweichung1(1,i) = AdB1(i)-AMess05(i);
disp(['Die Abweichung in dB bei i = ',num2str(i),' ist ',num2str(abweichung1(1,i))]);
end
eSquared1 = abweichung1.^2;
eSquaredSum1 = sum(eSquared1);


abweichung2 = zeros(1,length(f));
amplitude2 = zeros(1,length(f));
AdB2     = zeros(1,length(f));
phase2 = zeros(1,length(f));

for i = 1:length(f)
    % Filtern --> in der Filterfunktion muss simuliert werden!!!
    [simulationszeit,eingangssignal,ausgangssignalGefiltert] = berechneGefiltertesSignal(f(i),x,A);
    % Amplitude bestimmen
    [amplitude2(i),phase2(i),~] = berechneAmpPhs(simulationszeit,f(i),eingangssignal,ausgangssignalGefiltert);
    AdB2(i) = 20*log10(amplitude2(i)/A);
    abweichung2(1,i) = AdB2(i)-AMess25(i);
%     disp(['Die Abweichung in dB bei i = ',num2str(i),' ist ',num2str(abweichung1(1,i))])
end
eSquared2 = abweichung2.^2;
eSquaredSum2 = sum(eSquared2);


abweichung3 = zeros(1,length(f));
amplitude3 = zeros(1,length(f));
AdB3     = zeros(1,length(f));
phase3 = zeros(1,length(f));

for i = 1:length(f)
    % Filtern --> in der Filterfunktion muss simuliert werden!!!
    [simulationszeit,eingangssignal,ausgangssignalGefiltert] = berechneGefiltertesSignal(f(i),x,A);
    % Amplitude bestimmen
    [amplitude3(i),phase3(i),~] = berechneAmpPhs(simulationszeit,f(i),eingangssignal,ausgangssignalGefiltert);
    AdB3(i) = 20*log10(amplitude3(i)/A);
    abweichung3(1,i) = AdB3(i)-AMess50(i);
%     disp(['Die Abweichung in dB bei i = ',num2str(i),' ist ',num2str(abweichung1(1,i))])
end
eSquared3 = abweichung3.^2;
eSquaredSum3 = sum(eSquared3);


eGlobal = 4*eSquaredSum1+eSquaredSum2+eSquaredSum3;
disp(['eGlobal bei i = ',num2str(i),' ist ',num2str(eGlobal)])