
% need to import f-exp, p_exp, a_exp
f = logspace(1, 4, 100);
cutoff = 1590;
figure(1); clf
subplot(2,1,1)
plot(log10(f), -atan(f/cutoff)*57.3, 'LineWidth', 2)
hold on
plot(log10(f_exp), p_exp, '+', 'MarkerSize',10)
ylabel('phase')
title('Phase -- RC 10Kradians/sec, averaging 100000 samples')
subplot(2,1,2)
plot(log10(f), log10( 1./sqrt(1+f.^2/cutoff.^2)), 'LineWidth', 2)
hold on
plot(log10(f_exp), log10(a_exp), '+', 'MarkerSize',10)
ylabel('log amp')
xlabel('log freq')
title('Amp -- RC 10Kradians/sec, averaging 100000 samples')
