source: ETALON/Filters/Experiment results/Experiment3.m @ 713

%Experiment N3
close all; clear all; clc
nu_f=[20 23 24 25 26 27 30 40 50];
nu_f0=[37 53 75 124 186 368];
F=[416 368 364 390 420 422 382 352 328]/438;
F0=[350 230 180 104 60 24]/438;
figure
hold on
plot(nu_f,F,'-ob')
plot(nu_f0,F0,'-or')
%%%

F37=[74.9 74.4 81.33 83.06 75.6 83.3 90.74 80.42 79.78 75.7];
F53=[63.8 58.6 61.2 57.2 69.7 50.2 94.4 52.01 88.2 54.73];
F75=[51.5 35.7 47.13 43.57 53.47 30.16 64.12 41.9 52.5 35.53];
F124=[33.1 10.79 11.06 24.49 11.35 17.22 12.3 25 14.7 21.6];
F186=[11.38 4.3 9.15 21.64 7.89 6.7 10.3 24.4 13 25.3];
F368=[2.6 0.5 0.9 2.12 2.24 2.3 2.8 1.26 7.17 1.68];
Frame=[74.5 72.6 73.1 72.8 73.56 74.2 74.1];
Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)];
eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)];
figure
errorbar(nu_f0,Fx,eFx)


F37=[79.2 119.3 85.7 93.2 79.4 88 89.7 81.8];
F53=[40.7 38.5 51.2 48.3 47.1 34.6 47.5 53.5];
F75=[42.3 23.6 45.6 22.2 30.9 15.2 29.8 25.4];
F124=[18.5 8.18 17.4 16.7 6.15 9.84 7.0 7.5];
F186=[7.8 11.4 9.5 6.4 4.8 4.1 4.9 4.5 ];
F368=[11.03 3.57 5.36 5.7 7.1 2.53 4.5 2.3];
Frame=[95.8 92.5 92.7 92.05 92 92.1];
Fx=[mean(F37) mean(F53) mean(F75) mean(F124) mean(F186) mean(F368)];
eFx=[std(F37) std(F53) std(F75) std(F124) std(F186) std(F368)];
figure
errorbar(nu_f0,Fx,eFx)
%%%%
F20=[89.1 84.8 83.2 81.2];
F23=[72 71.4 70.3 72];
F24=[67.7 71 70.5 72.7 ];
F25=[74.5 74.6 74.2 74.2 ];
F26=[81.7 82 81 82.4];
F27=[88 84.2 87.8 88.5];
F30=[77.6 77.4 76.6 79.7];
F40=[61.3 60.21 60.9 62];
F50=[66.2 66 67.3 66];
Frame=[94.3 93.8 93 92.3 93.7 93.2 93.1];
Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame);
eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame);
figure
errorbar(nu_f,Fx,eFx)

F20=[94.7 95.4 96.2 94.2 ];
F23=[ 90.3 93.3 87 86.6 ];
F24=[ 87.6 91.6 88.4 89.5 ];
F25=[87.7 88 88.3 87.5];
F26=[85.2 83.9 83.1 84.6];
F27=[84.9 83 80.5 77.3];
F30=[82.6 82 83.3 85.5];
F40=[87 89.5 89.2 88];
F50=[91.4 94.2 94.8 90.7];
Frame=[96.3 93.3 92.7 92 91.8 91.6 91 91.2];
Fx=[mean(F20) mean(F23) mean(F24) mean(F25) mean(F26) mean(F27) mean(F30) mean(F40) mean(F50)]./mean(Frame);
eFx=[std(F20) std(F23) std(F24) std(F25) std(F26) std(F27) std(F30) std(F40) std(F50)]./mean(Frame);
figure
errorbar(nu_f,Fx,eFx)