logs = logspace(0.1,2.301, 15);

magnitudes = [];
phases = [];
for i = 1:15
    data = load('meas' + string(i));
    magnitudes(i) = 20*log10(max(data.simout.signals.values(:,5))/5);
    phases(i) = givepd(data.simout.signals.values(:,5), ...
                       data.simout.signals.values(:,1)) * (180/pi);
end
figure(2)
semilogx(logs,magnitudes);
semilogx(logs,phases);

kc = 0.1586;
magnitudes_kc = magnitudes + (20*log10(kc));
semilogx(logs,magnitudes_kc);
semilogx(logs,phases);