%%
clear all
close all
clc
% Ausgangswerte zur Messdatenerzeugung
% K   = 5;
% w0  = 50;
% zeta= 0.08;
% koeff = [K w0 zeta];

% w = 1:0.01:300;
% A = 20*log10(koeff(1)*koeff(2).^2./sqrt((koeff(2).^2-w.^2).^2+(2*koeff(3)*koeff(2)*w).^2));
% semilogx(w,A);
% hold on
% bode(tf([K*w0^2],[1 2*zeta*w0 w0^2]))

% opt = optimset('Display','iter','TolFun',850,'MaxIter',20);
opt = optimset('TolFun',850);
x0  = [11 0.1 50];
x = fminsearch(@(x) berechneEGlobal(x),x0,opt)